Search USGSSearch

Geology topics

Robert E. Kayen

Publications and source records attributed to Robert E. Kayen.

At least 19 recordsLinked to original sources

Point clouds of bridge generated by terrestrial laser scanner and images via Structure from Motion technique: Comparison study

A collection of points representing an object in space is commonly called a point cloud. There are several techniques for collecting point clouds. This research is focused on a comparison study of two approaches: (1) collecting point clouds with a surveying grade terrestrial laser scanner (TLS) and (2) generating point clouds from drone-taken still images by utilizing Structure-from-Motion (SFM) technique. The paper's main objective is to compare the generated point clouds to each other and show the advantages and disadvantages of both techniques for structural health monitoring of bridges. To achieve the goal of the paper, we selected a pedestrian bridge for a comparison study. For simplicity, the study was limited to a single-span bridge. The bridge we studied is a composite steel and reinforced concrete (RC) bridge with a curved deck. A terrestrial laser scanner was used to scan the bridge from several positions, and point clouds were registered based on point-cloud to point-cloud matching. No targets were used during the collection of the point clouds by the laser scanner. In addition, still images of both bridges were taken by a drone. Based on a quantitative comparison of the results, the paper discusses the pros and cons of both approaches.

California

Shoreline retreat of the Corte Madera marshes, 1853 to 2016, Marin County, California

The greater San Francisco Bay estuary, prior to human intervention, encompassed about 2,200 km 2 of tidal and salt marshes. Over time, these areas became increasingly diked, developed, and altered from their natural state. In addition, natural forces are always driving a continually shifting equilibrium. This study area, the Corte Madera marshes, is a tidal marsh or wetland located in southeastern Marin County, and it borders an embayment of central San Francisco Bay along about 2.8 km of shoreline. Most of this shoreline is located within the Corte Madera Marsh Ecological Reserve, managed by the California Department of Fish and Wildlife. Other areas within the marsh include (1) unincorporated Greenbrae (at the boardwalk), (2) diked land (that is, isolated from tidal action) owned by the Golden Gate Bridge Highway and Transportation District, and (3) urbanized areas such as in the Mariner Cove subdivision of Corte Madera. The present tidal marsh area was historically subdivided into the following informally named tracts, listed from north to south: Heerdt marsh, north Muzzi marsh, inner and outer Muzzi marshes, Marta’s marsh, and Triangle marsh. The purpose of this study is to derive the magnitudes and rates of shoreline change (both erosion and accretion) for the Corte Madera shoreline, with particular emphasis on the time period from 1931 to 2016. The rates of change are then related to different shoreline types (that is, natural or diked) and (or) locations on the shoreline.

California

Examination of differences between three SPT-based seismic soil liquefaction triggering relationships

The preceding companion paper presented the updating of the seismic soil liquefaction triggering relationship of Cetin et al. [1], and compared the resulting updated relationship with the earlier version. In this second paper, a detailed cross-comparison is made between three triggering relationships: (1) Seed et al. [2], as slightly updated by the NCEER Working Group (Youd et al. [3]), (2) Boulanger and Idriss [4], and (3) Cetin et al. [5]. Differences between these three triggering relationships, and the apparent causes of them are examined. Also studied are the impacts of these differences on levels of conservatism with regard to evaluation of liquefaction triggering hazard, and the resulting risks for engineering projects.

Soil Dynamics and Earthquake Engineering

Seismic velocity site characterization of 10 Arizona strong-motion recording stations by spectral analysis of surface wave dispersion

Vertical one-dimensional shear wave velocity ( V S ) profiles are presented for strong-motion sites in Arizona for a suite of stations surrounding the Palo Verde Nuclear Generating Station. The purpose of the study is to determine the detailed site velocity profile, the average velocity in the upper 30 meters of the profile ( V S30 ), the average velocity for the entire profile ( V SZ ), and the National Earthquake Hazards Reduction Program (NEHRP) site classification. The V S profiles are estimated using a non-invasive continuous-sine-wave method for gathering the dispersion characteristics of surface waves. Shear wave velocity profiles were inverted from the averaged dispersion curves using three independent methods for comparison, and the root-mean-square combined coefficient of variation ( COV ) of the dispersion and inversion calculations are estimated for each site.

Arizona

Seismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading

Gentle sediment-laden slopes are typical of the onshore coastal zone and offshore continental shelf and slope. Coastal sediment are commonly young weakly consolidated materials that are well stratified, have low strength, and can mobilize shear displacements at low levels of stress. Seismically-driven plastic displacements of these sediment pose a hazard to coastal cities, buried onshore utilities, and offshore infrastructure like harbor protection and outfalls. One-dimensional rigid downslope-directed Newmark sliding block analyses have been used to predict earthquake deformations generally on steeper slopes that are modeled as frictional materials. This study probes the effect of multidirectional earthquake motions on inertial displacements of gently sloping ground of the coastal and offshore condition where soft-compliant soil is expected. Toward that objective, this investigation seeks to understand the effect on Newmark-type displacements of [1] multidirectional earthquake shaking and [2] soil compliance. In order to model multidirectional effects, the earthquake motions are rotated into the local slope strike- and dip-components. On gently sloping ground, including the strike component of motion always results in a larger and more accurate shear stress vector. Strike motions are found to contribute to downslope deformations on any declivity. Compliant response of the soil mass also influences the plastic displacements. The magnitude of seismic displacements can be estimated with a simplified model using only the estimated soil yield-acceleration ( k y ) and the peak ground velocity ( V max ) of the earthquake motions. Compliance effects can be effectively mapped using the concept of Plastic Displacement Response Spectra (PDRS).

Engineering Geology

Summary of SPT based field case history data of CETIN (2016) database

This report provides documentation of the Cetin et al. (2016) field performance case histories, probabilistic maximum likelihood assessment and the sources of differences between the liquefaction triggering resistance estimations (CRR values) of the widely used liquefaction triggering relationships of Seed et al. (1985), Cetin et al. (2004, 2016) and Boulanger and Idriss (2012). Cetin et al. (2016-a) presented a concise summary of the improved database and the updated triggering relationships. For the sake of completeness, as part of Chapter 2 of this report, Cetin et al. (2016-a) manuscript is re-presented, as edited to include the electronic supplements, and a broader and more detailed documentation of all of case histories. This broader and more detailed documentation along with the interpretations of other researchers is presented in Appendix A. Similarly, site response analyses details and results from Cetin (2000) is re-presented in Appendix B. To enable the readers to quickly refer to Cetin and Seed (2004), it is also included in Appendix C. The sources of differences between the liquefaction triggering resistance estimations (CRR values) of the widely used liquefaction triggering relationships of Seed et al. (1985), Cetin et al. (2004) and Boulanger and Idriss (2012) were discussed in Cetin et al. (2016-b). Again, for the sake of completeness, as part of Chapter 3 of this report, Cetin et al. (2016-b) manuscript is re-presented, as edited to the include of electronic supplements, and a detailed response to issues addressed by Idriss and Boulanger (2012). Also, the technical response letter (Cetin, 2014) prepared at the request of the members of the NRC Committee on the State of the Art and Practice in Earthquake Induced Soil Liquefaction Assessment, is given in Appendix D, which is helpful to follow the technical discussions. This report is intended only as a concise summary of a vast amount of data. The interpretations presented are those of the research team. A more detailed description of some of the details of the methods and procedures used to evaluate and analyze these field performance case histories is also presented in Cetin (2000) and Cetin et al. (2000) though the final evaluations presented in this report are the most recent interpretations undertaken under these current studies.

Report

An overview of the geotechnical damage brought by the 2016 Kumamoto Earthquake, Japan

The 2016 Kumamoto earthquake with a moment magnitude of 7.0 (Japanese intensity = 7) that struck on April 16 brought devastation in many areas of Kumamoto Prefecture and partly in Oita Prefecture in Kyushu Region, Japan. The earthquake succeeds a foreshock of magnitude 6.5 (Japanese intensity = 7) on April 14. The authors conducted two surveys on the devastated areas: one during April 16-17, and the other during May 11-14. This report summarizes the damage brought to geotechnical structures by the two consecutive earthquakes within a span of twenty-eight hours. This report highlights some of the observed damage and identifies reasons for such damage. The geotechnical challenges towards mitigation of losses from such earthquakes are also suggested.

Report

A comparative study of surface waves inversion techniques at strong motion recording sites in Greece

Surface wave method was used for the estimation of Vs vs depth profile at 10 strong motion stations in Greece. The dispersion data were obtained by SASW method, utilizing a pair of electromechanical harmonic-wave source (shakers) or a random source (drop weight). In this study, three inversion techniques were used a) a recently proposed Simplified Inversion Method (SIM), b) an inversion technique based on a neighborhood algorithm (NA) which allows the incorporation of a priori information regarding the subsurface structure parameters, and c) Occam's inversion algorithm. For each site constant value of Poisson's ratio was assumed (ν=0.4) since the objective of the current study is the comparison of the three inversion schemes regardless the uncertainties resulting due to the lack of geotechnical data. A penalty function was introduced to quantify the deviations of the derived Vs profiles. The Vs models are compared as of Vs(z), Vs30 and EC8 soil category, in order to show the insignificance of the existing variations. The comparison results showed that the average variation of SIM profiles is 9% and 4.9% comparing with NA and Occam's profiles respectively whilst the average difference of Vs30 values obtained from SIM is 7.4% and 5.0% compared with NA and Occam's.

Conference Paper

Shear-wave velocity and site-amplification factors for 50 Australian sites determined by the spectral analysis of surface waves method

One-dimensional shear-wave velocity ( V S ) profiles are presented at 50 strong motion sites in New South Wales and Victoria, Australia. The V S profiles are estimated with the spectral analysis of surface waves (SASW) method. The SASW method is a noninvasive method that indirectly estimates the V S at depth from variations in the Rayleigh wave phase velocity at the surface.

Open-File Report

Shear Wave Velocity and Site Amplification Factors for 25 Strong-Motion Instrument Stations Affected by the M5.8 Mineral, Virginia, Earthquake of August 23, 2011

Vertical one-dimensional shear wave velocity (Vs) profiles are presented for 25 strong-motion instrument sites along the Mid-Atlantic eastern seaboard, Piedmont region, and Appalachian region, which surround the epicenter of the M5.8 Mineral, Virginia, Earthquake of August 23, 2011. Testing was performed at sites in Pennsylvania, Maryland, West Virginia, Virginia, the District of Columbia, North Carolina, and Tennessee. The purpose of the study is to determine the detailed site velocity profile, the average velocity in the upper 30 meters of the profile (VS,30), the average velocity for the entire profile (VS,Z), and the National Earthquake Hazards Reduction Program (NEHRP) site classification. The Vs profiles are estimated using a non-invasive continuous-sine-wave method for gathering the dispersion characteristics of surface waves. A large trailer-mounted active source was used to shake the ground during the testing and produce the surface waves. Shear wave velocity profiles were inverted from the averaged dispersion curves using three independent methods for comparison, and the root-mean square combined coefficient of variation (COV) of the dispersion and inversion calculations are estimated for each site.

North Carolina, Pennslyvannia, Tennesee, West Virg

Surface wave site characterization at 27 locations near Boston, Massachusetts, including 2 strong-motion stations

The geotechnical properties of the soils in and around Boston, Massachusetts, have been extensively studied. This is partly due to the importance of the Boston Blue Clay and the extent of landfill in the Boston area. Although New England is not a region that is typically associated with seismic hazards, there have been several historical earthquakes that have caused significant ground shaking (for example, see Street and Lacroix, 1979; Ebel, 1996; Ebel, 2006). The possibility of strong ground shaking, along with heightened vulnerability from unreinforced masonry buildings, motivates further investigation of seismic hazards throughout New England. Important studies that are pertinent to seismic hazards in New England include source-parameter studies (Somerville and others, 1987; Boore and others, 2010), wave-propagation studies (Frankel, 1991; Viegas and others, 2010), empirical ground-motion prediction equations (GMPE) for computing ground-motion intensity (Tavakoli and Pezeshk, 2005; Atkinson and Boore, 2006), site-response studies (Hayles and others, 2001; Ebel and Kim, 2006), and liquefaction studies (Brankman and Baise, 2008). The shear-wave velocity (VS) profiles collected for this report are pertinent to the GMPE, site response, and liquefaction aspects of seismic hazards in the greater Boston area. Besides the application of these data for the Boston region, the data may be applicable throughout New England, through correlations with geologic units (similar to Ebel and Kim, 2006) or correlations with topographic slope (Wald and Allen, 2007), because few VS measurements are available in stable tectonic regions. Ebel and Hart (2001) used felt earthquake reports to infer amplification patterns throughout the greater Boston region and noted spatial correspondence with the dominant period and amplification factors obtained from ambient noise (horizontal-to-vertical ratios) by Kummer (1998). Britton (2003) compiled geotechnical borings in the area and produced a microzonation map based on generalized velocity profiles, where the amplifications were computed using Shake (Schnable and others, 1972), along with an assumed input ground motion. The velocities were constrained by only a few local measurements associated with the Central Artery/Tunnel project. The additional VS measurements presented in this report provide a number of benefits. First, these measurements provide improved spatial coverage. Second, the larger sample size provides better constraints on the mean and variance of the VS distribution for each layer, which may be paired with a three-dimensional (3D) model of the stratigraphy to generate one-dimensional (1D) profiles for use in a standard site-response analysis (for example, Britton, 2003). Third, the velocity profiles may also be used, along with a 3D model of the stratigraphy, as input into a 3D simulation of the ground motion to investigate the effects of basin-generated surface waves and the potential focusing of seismic waves. This report begins with a short review of the geology of the study area and the field methods that we used to estimate the velocity profiles. The raw data, processed data, and the interpreted VS profiles are given in appendix 1. Photographs and descriptions of the sites are provided in appendix 2.

Massachusetts

Settlement of the USS Arizona, Pearl Harbor, Hawaii

The U.S. Geological Survey, in collaboration with the National Park Service Submerged Resources Center, undertook investigations at the USS Arizona Memorial at Pearl Harbor, Hawaii, in 2002, 2003, and 2005 to characterize geological factors affecting the deterioration and movement of the hull of the USS Arizona. Since sinking on the morning of December 7, 1941, the hull of the USS Arizona has been slowly but steadily disappearing below the surface of Pearl Harbor. Continuous sediment coring at three of four locations around the hull of the Arizona was only partially successful, but it was sufficient to identify a varied sedimentary substrate beneath the hull. A boring near the stern reveals a thick, continuous sequence of soft, gray clay to the bottom of the boring. In contrast, borings near the bow and starboard side, below about 5 meters subbottom depth, indicate the presence of very stiff, brown clay and coral debris and an absence of soft clay. Multisensor core logger scanning of the recovered cores distinguishes the lower density of the soft, gray clay at the stern from the higher density of the stiff, brown clays and coral debris at the bow and starboard side. Uniaxial consolidation testing of the soft gray clay indicates a normally consolidated sequence, whereas the stiff, brown clay and coral debris are overconsolidated. Profiles of shear wave velocity vs. depth obtained through spectral analysis of interface wave testing around the perimeter of the hull in 2005 identified areas of higher velocity, stiffer sediment at the bow and starboard side, which correspond to the dense, stiff clay recovered near the bow and starboard borings. Low shear-wave velocities at the port midship and quarter of the hull correlate with the lower density, softer sediment recovered from the boring at the stern. Cross sections of the subbottom of the Memorial combine results from the sediment borings and geophysical surveys and depict a wedge of soft clay unconformably overlying the stiff clays and coral debris beneath the aft half of the USS Arizona and thickening toward the stern. The 2008 position of the hull has been documented using both tide-based and differential Global Positioning System (GPS) measuring systems. Analysis of historical and recent photographs was done to create a record of settlement from the time of sinking in 1941 to the present. By examining shadows in suitable photos, the sun azimuth, local time of day, and tide levels were determined to derive tide-adjusted and sea-level-rise-corrected elevations for structures on the hull and from these elevations to obtain settlement and tilt trends. The settlement trends, most complete for barbette 3, have two components. An early, nonlinear component ends on December 9, 1941, and represents the initial penetration and displacement of the bottom sediment by the hull. A linear, long-term trend of normal consolidation continues to the present day. Long-term settlement rates are greatest at the stern and decrease linearly to the midship, showing that the aft half of the hull is moving as an intact, rigid body. The recent rate of settlement at the stern is about 3.5 mm/year; rates at the starboard midship and forward part of the hull are less than one-third of the stern rate. The aft half of the USS Arizona hull presently tilts about 2 degrees to port, an increase of at least 1.5 degrees since the initial sinking of the ship. The results of this study identify differential settlement of the Arizona hull, due to the wedge of soft clay underlying the aft half of the hull, as the cause of the movement of the hull beneath the surface of Pearl Harbor. Calculation of sediment consolidation using lab-determined properties of the soft clay demonstrates that the observed settlements can be reproduced by projecting appropriate clay thicknesses beneath the hull. Several of the high-quality photographs analyzed for the historical settlement analysis highlight some of the limitations of this retrospective technique for determining tide-based elevations. In these cases, calculated structure elevations do not conform to the settlement trend, indicating that there can be complicating factors affecting the interpretation of the photos. Conflicting dates for events during the salvage operations were also encountered.

Hawai'i

Non-destructive measurement of soil liquefaction density change by crosshole radar tomography, Treasure Island, California

A ground penetrating radar (GPR) experiment at the Treasure Island Test Site [TILT] was performed to non-destructively image the soil column for changes in density prior to, and following, a liquefaction event. The intervening liquefaction was achieved by controlled blasting. A geotechnical borehole radar technique was used to acquire high-resolution 2-D radar velocity data. This method of non-destructive site characterization uses radar trans-illumination surveys through the soil column and tomographic data manipulation techniques to construct radar velocity tomograms, from which averaged void ratios can be derived at 0.25 - 0.5m pixel footprints. Tomograms of void ratio were constructed through the relation between soil porosity and dielectric constant. Both pre- and post-blast tomograms were collected and indicate that liquefaction related densification occurred at the site. Volumetric strains estimated from the tomograms correlate well with the observed settlement at the site. The 2-D imagery of void ratio can serve as high-resolution data layers for numerical site response analysis.

Book

Spatial distribution of landslides triggered from the 2007 Niigata Chuetsu–Oki Japan Earthquake

Understanding the spatial distribution of earthquake-induced landslides from specific earthquakes provides an opportunity to recognize what to expect from future events. The July 16, 2007 M w 6.6 (M JMA 6.8) Niigata Chuetsu–Oki Japan earthquake triggered hundreds of landslides in the area surrounding the coastal city of Kashiwazaki and provides one such opportunity to evaluate the impacts of an offshore, magnitude 6 + earthquake on a steep coastal region. As part of a larger effort to document all forms of geotechnical damage from this earthquake, we performed landslide inventory mapping throughout the epicentral area and analyzed the resulting data for spatial, seismic-motion, and geologic correlations to describe the pattern of landsliding. Coupled with examination of a third-party, aerial-photo-based landslide inventory, our analyses reveal several areas of high landslide concentration that are not readily explained by either traditional epicentral and fault–plane-distance metrics or by recorded and inferred ground-motions. Whereas average landslide concentrations averaged less than 1 landslide per square kilometer (LS/km 2 ), some areas reached up to 2 LS/km 2 in the Nishiyama Hills to the northeast of Kashiwazaki and between 2 and 11 LS/km 2 in coastal areas to the north and south of the city. Correlation with seismometer-based and monument overturning back-calculated ground motions suggests that a minimum peak ground acceleration (PGA) of approximately 0.2 g was necessary for landsliding throughout the region, but does not explain the subregional areas of high landslide concentration. However, analysis of topographic slope and the distribution of generally weak, dip-slope, geologic units does sufficiently explain why, on a sub-regional scale, high landslide concentrations occurred where they did. These include: (1) an inland region of steep, dip-slope, anticlinal sedimentary strata with associated fold belt compression and uplift of the anticline and (2) coastal areas with generally weaker, weathered outcrop lithology and steeper slopes resulting from active marine and terrestrial cliff processes. The results offer lessons for understanding the effects of earthquakes on both regional and subregional scales with regard to the spatial distribution of landsliding.

Engineering Geology

Volcano collapse promoted by progressive strength reduction: New data from Mount St. Helens

Rock shear strength plays a fundamental role in volcano flank collapse, yet pertinent data from modern collapse surfaces are rare. Using samples collected from the inferred failure surface of the massive 1980 collapse of Mount St. Helens (MSH), we determined rock shear strength via laboratory tests designed to mimic conditions in the pre-collapse edifice. We observed that the 1980 failure shear surfaces formed primarily in pervasively shattered older dome rocks; failure was not localized in sloping volcanic strata or in weak, hydrothermally altered rocks. Our test results show that rock shear strength under large confining stresses is reduced ∼20% as a result of large quasi-static shear strain, as preceded the 1980 collapse of MSH. Using quasi-3D slope-stability modeling, we demonstrate that this mechanical weakening could have provoked edifice collapse, even in the absence of transiently elevated pore-fluid pressures or earthquake ground shaking. Progressive strength reduction could promote collapses at other volcanic edifices.

Washington

Historic and paleo-submarine landslide deposits imaged beneath Port Valdez, Alaska: Implications for tsunami generation in a glacial fiord

During the 1964 M9.2 great Alaskan earthquake, submarine-slope failures resulted in the generation of highly destructive tsunamis at Port Valdez, Alaska. A high-resolution, mini-sparker reflection profiler was used to image debris lobes, which we attribute to slope failures that occurred both during and prior to the 1964 megathrust event. In these reflection profiles, debris lobe deposits are indicated by acoustically opaque units that are separated by undisturbed parallel-layered reflectors. Near-surface debris lobes attributed to the 1964 earthquake include: (1) a debris lobe over 30 m thick that emanates from the fiord-head delta in eastern Port Valdez; and (2) debris flow lobes incorporating large, intact blocks up to 40 m high in western Port Valdez, off the Shoup Glacier moraine. In addition to the near-surface debris lobes, we imaged at least five additional debris lobe deposits buried beneath the 1964 deposit. The debris lobe directly beneath the 1964 deposit has a similar thickness and spatial distribution as the 1964 deposit. However, the older, deeper, debris lobes are thinner, less extensive, and separated by thinner sequences of parallel-layered reflectors. Glacier retreat and concomitant build-up of the fiord-head delta combined with longer time intervals between megathrust events may have resulted in more extensive delta failures and thus thicker debris lobes through time.

Alaska

Surface-wave site characterization at 52 strong-motion recording stations affected by the Parkfield, California, M6.0 earthquake of 28 September 2004

We present one-dimensional shear-wave velocity ( V S ) profiles at 52 strong-motion sites that recorded the 28 September 2004 Magnitude 6.0 Parkfield, Calif., earthquake. We estimate the V S profiles with the Spectral Analysis of Surface-Waves (SASW) method. The SASW method is a noninvasive method that indirectly estimates the V S at depth from variations in the Rayleigh wave phase velocity at the surface. To address the uncertainty associated with these measurements, we compare the SASW profiles to surface-source downhole-receiver (SSDR) profiles at four sites. Three of the four SSDR sites are in close agreement with the adjacent SASW site, while the SASW profile is considerably slower than the SSDR profile at one site.

California

Three-dimensional site response at KiK-net downhole arrays

Ground motions at two Kiban-Kyoshin Network (KiK-net) strong motion downhole array sites in Hokkaido, Japan (TKCH08 in Taiki and TKCH05 in Honbetsu) illustrate the importance of three-dimensional (3D) site effects. These sites recorded the M8.0 2003 Tokachi-Oki earthquake, with recorded accelerations above 0.4 g at both sites as well as numerous ground motions from smaller events. Weak ground motions indicate that site TKCH08 is well modeled with the assumption of plane SH waves traveling through a 1D medium (SH1D), while TKCH05 is characteristic of a poor fit to the SH1D theoretical response. We hypothesized that the misfit at TKCH05results from the heterogeneity of the subsurface. To test this hypothesis, we measured four S-wave velocity profiles in the vicinity (< 300 m) of each site with the spectral analysis of surface waves (SASW) method. This KiK-net site pair is ideal for assessing the relative importance of 3D site effects and nonlinear site effects. The linear ground motions at TKCH05 isolate the 3D site effects, as we hypothesized from the linear ground motions and confirmed with our subsequent SASW surveys. The Tokachi-Oki time history at TKCH08 isolates the effects of nonlinearity from spatial heterogeneity because the 3D effects are negligible. The Tokachi-Oki time history at TKCH05 includes both nonlinear and 3D site effects. Comparisons of the accuracy of the SH1D model predictions of these surface time histories from the downhole time histories indicates that the 3D site effects are at least as important as nonlinear effects in this case. The errors associated with the assumption of a 1D medium and 1D wave propagation will be carried into a nonlinear analysis that relies on these same assumptions. Thus, the presence of 3D effects should be ruled out prior to a 1D nonlinear analysis. The SH1D residuals show that 3D effects can be mistaken for nonlinear effects.

Hokkaido