Paul Bodin

Strong Ground Motion from the Michoacan, Mexico, Earthquake

The network of strong motion accelerographs in Mexico includes instruments that were installed, under an international cooperative research program, in sites selected for the high potenial of a large earthquake. The 19 September 1985 earthquake (magnitude 8.1) occurred in a seismic gap where an earthquake was expected. As a result, there is an excellent descripton of the ground motions that caused the disaster.

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δt (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t0, and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t0. The approximate solutions underpredict the full solutions, although the difference decreases as t0 approaches the end of the earthquake cycle. The relationship between Δt and t0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δt, transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t0. We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δt equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δt than the ΔCFS model does for static stress steps at small t0, and smaller Δt than the ΔCFS model for stress steps at large t0. The ΔCFS model predicts nonzero Δt only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δt for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Non-volcanic tremor driven by large transient shear stresses

Non-impulsive seismic radiation or ‘tremor’ has long been observed at volcanoes and more recently around subduction zones. Although the number of observations of non-volcanic tremor is steadily increasing, the causative mechanism remains unclear. Some have attributed non-volcanic tremor to the movement of fluids, while its coincidence with geodetically observed slow-slip events at regular intervals has led others to consider slip on the plate interface as its cause. Low-frequency earthquakes in Japan, which are believed to make up at least part of non-volcanic tremor, have focal mechanisms and locations that are consistent with tremor being generated by shear slip on the subduction interface. In Cascadia, however, tremor locations appear to be more distributed in depth than in Japan, making them harder to reconcile with a plate interface shear-slip model. Here we identify bursts of tremor that radiated from the Cascadia subduction zone near Vancouver Island, Canada, during the strongest shaking from the moment magnitude Mw = 7.8, 2002 Denali, Alaska, earthquake. Tremor occurs when the Love wave displacements are to the southwest (the direction of plate convergence of the overriding plate), implying that the Love waves trigger the tremor. We show that these displacements correspond to shear stresses of approximately 40 kPa on the plate interface, which suggests that the effective stress on the plate interface is very low. These observations indicate that tremor and possibly slow slip can be instantaneously induced by shear stress increases on the subduction interface—effectively a frictional failure response to the driving stress.

Observing earthquakes triggered in the near field by dynamic deformations

We examine the hypothesis that dynamic deformations associated with seismic waves trigger earthquakes in many tectonic environments. Our analysis focuses on seismicity at close range (within the aftershock zone), complementing published studies of long-range triggering. Our results suggest that dynamic triggering is not confined to remote distances or to geothermal and volcanic regions. Long unilaterally propagating ruptures may focus radiated dynamic deformations in the propagation direction. Therefore, we expect seismicity triggered dynamically by a directive rupture to occur asymmetrically, with a majority of triggered earthquakes in the direction of rupture propagation. Bilaterally propagating ruptures also may be directive, and we propose simple criteria for assessing their directivity. We compare the inferred rupture direction and observed seismicity rate change following 15 earthquakes ( M 5.7 to M 8.1) that occurred in California and Idaho in the United States, the Gulf of Aqaba, Syria, Guatemala, China, New Guinea, Turkey, Japan, Mexico, and Antarctica. Nine of these mainshocks had clearly directive, unilateral ruptures. Of these nine, seven apparently induced an asymmetric increase in seismicity rate that correlates with the rupture direction. The two exceptions include an earthquake preceded by a comparable-magnitude event on a conjugate fault and another for which data limitations prohibited conclusive results. Similar (but weaker) correlations were found for the bilaterally rupturing earthquakes we studied. Although the static stress change also may trigger seismicity, it and the seismicity it triggers are expected to be similarly asymmetric only if the final slip is skewed toward the rupture terminus. For several of the directive earthquakes, we suggest that the seismicity rate change correlates better with the dynamic stress field than the static stress change. Manuscript received 1 March 2002.

Ground-Motion Scaling in the Kachchh Basin, India, Deduced from Aftershocks of the 2001 Mw 7.6 Bhuj Earthquake

We studied the excitation, propagation, and site effects in the Kachchh basin of India by using ground-motion recordings from a temporary seismograph network deployed to study aftershocks of the M w 7.6 Bhuj earthquake of 26 January 2001. The Kachchh basin has been proposed as a useful analog region for studying hazard in other earthquake-prone but slowly deforming regions, such as the central United States. The earthquakes we studied ranged in size from about M 2 to M 5.2, and travel paths ranged from a few kilometers to about a hundred kilometers. There was a broad range of focal depths among the aftershocks, so the data were divided into two overlapping subsets to test the sensitivity of the derived propagation and source parameters to focal depth. Parameters we constrained include the source excitation terms (related to stress drop), a frequency-dependent attenuation operator, a geometric spreading function, and an operator to account for site effects. Our results indicate that seismic-wave attenuation in Kachchh crust is very low, similar to other continental intraplate areas such as central and eastern North America. We also estimated seismic moments and stress drops for the earthquakes by fitting single-corner-frequency source-model spectra to the observed spectra, corrected for propagation by using our derived parameters. Stress drops were found to scale with seismic moment and to be rather high overall. By using a stochastic point-source model to estimate mainshock ground motions, we found that the distance decay of expected peak ground motions, assuming a stress drop of 15-20 MPa, compare well with the scant observations for the Bhuj earthquake. Ground-motion predictions for Kachchh, based on Bhuj aftershock data, support the idea that the region may have similar hazard to proposed analog areas in North America.

Coastal uplift and mortality of intertidal organisms caused by the september 1985 Mexico earthquakes.

Coastal uplift associated with the great Mexican earthquake of 19 September 1985 and its principal aftershock produced widespread mortality of intertidal organisms along the coast of the states of Michoac�n and Guerrero, Mexico. Measurements of the vertical extent of mortality at ten sites provided estimates of the magnitude of the vertical component of deformation along the coast. Within the affected area, uplift ranged from about 12 centimeters to about 1 meter, and no subsidence was observed. The observations are consistent with models of the tectonic deformation that results from buried slip on a shallow-dipping underthrust fault.

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Noise to signal: A microtremor study at liquefaction sites in the New Madrid Seismic Zone

Understanding how sedimentary basins respond to seismic-wave energy generated by large earthquake events is a significant concern for seismic-hazard estimation. This study explores the use of microtremors, or ambient noise, for evaluating strong-motion site effects. The study focuses on the Mississippi Embayment in the New Madrid Seismic Zone, where widespread liquefaction and ground failure occurred during the 1811–1812 earthquake sequence. Spectral analyses of microtremor data at sites representing different environments of deposition (and sedimentary facies), different embayment thicknesses, and varying liquefaction susceptibility show correlations between (1) calculated vulnerability indices and evidence of liquefaction, (2) sediment thickness and fundamental resonant frequency, and (3) subsurface stratigraphic boundaries and observed peaks in horizontal-to-vertical spectral ratios. Results of the study suggest that the microtremor method could be helpful in identifying those areas most vulnerable to ...

Induced Dynamic Nonlinear Ground Response at Garner Valley, California

We present results from a prototype experiment in which we actively induce, observe, and quantify in situ nonlinear sediment response in the near surface. This experiment was part of a suite of experiments conducted during August 2004 in Garner Valley, California, using a large mobile shaker truck from the Network for Earthquake Engineering Simulation (NEES) facility. We deployed a dense accel- erometer array within meters of the mobile shaker truck to replicate a controlled, laboratory-style soil dynamics experiment in order to observe wave-amplitude- dependent sediment properties. Ground motion exceeding 1g acceleration was pro- duced near the shaker truck. Thewave field was dominated by Rayleigh surface waves and ground motions were strong enough to produce observable nonlinear changes in wave velocity. We found that as the force load of the shaker increased, the Rayleigh- wave phase velocity decreased by as much as ∼30% at the highest frequencies used (up to 30 Hz). Phase velocity dispersion curves were inverted for S-wave velocity as a function of depth using a simple isotropic elastic model to estimate the depth depen- dence of changes to the velocity structure. The greatest change in velocity occurred nearest the surface, within the upper 4 m. These estimated S-wave velocity values were used with estimates of surface strain to compare with laboratory-based shear modulus reduction measurements from the same site. Our results suggest that it may be possible to characterize nonlinear soil properties in situ using a noninvasive field technique.

Slidequake Generation versus Viscous Creep at Softrock-landslides: Synopsis of Three Different Scenarios at Slumgullion Landslide, Heumoes Slope, and Super-Sauze Mudslide

In this study, we describe conditions for slidequake generation at three different creeping softrock landslides: the Slumgullion landslide in the San Juan Mountains, Colorado, U.S., the Heumoes slope in the Austrian Alps, and the mudslide in Super-Sauze, French Alps. From a geomorphologic point of view, all three landslides are classified as creeping landslides with average velocities between centimeters to meters per year. Associating creep with viscous flow, and considering the largely saturated, clayey consistency of the slope body, one would not expect any brittle behavior. Thus, it came as a surprise that impulsive seismic signals indicative of shear fracture could be discovered by sensitive passive monitoring methods at all three slopes. These fracture signals occur in episodes, have similar signatures as small earthquakes, and could be located within the slide bodies, i.e., are evidence of slidequakes. Our investigations identified seismic and aseismic slip in each slide, with slidequakes focusing at significant bedrock structures or at lateral boundaries. Synoptic comparison of three scenarios underlines the importance of landslide-bedrock and landslide-lateral boundary interactions under gravitational loading and Mohr-Coulomb-type failure. Comparison to frictional and asperity models of crustal- and plate-scale boundaries may pave the way to a comprehensive understanding of slidequake generation, and future slope failure prediction.