Search USGS⌕ Search

Geology topics

N. Beeler

Publications and source records attributed to N. Beeler.

3 recordsLinked to original sources

Laboratory-based maximum slip rates in earthquake rupture zones and radiated energy

Laboratory stick-slip friction experiments indicate that peak slip rates increase with the stresses loading the fault to cause rupture. If this applies also to earthquake fault zones, then the analysis of rupture processes is simplified inasmuch as the slip rates depend only on the local yield stress and are independent of factors specific to a particular event, including the distribution of slip in space and time. We test this hypothesis by first using it to develop an expression for radiated energy that depends primarily on the seismic moment and the maximum slip rate. From laboratory results, the maximum slip rate for any crustal earthquake, as well as various stress parameters including the yield stress, can be determined based on its seismic moment and the maximum slip within its rupture zone. After finding that our new equation for radiated energy works well for laboratory stick-slip friction experiments, we used it to estimate radiated energies for five earthquakes with magnitudes near 2 that were induced in a deep gold mine, an M 2.1 repeating earthquake near the San Andreas Fault Observatory at Depth (SAFOD) site and seven major earthquakes in California and found good agreement with energies estimated independently from spectra of local and regional ground-motion data. Estimates of yield stress for the earthquakes in our study range from 12 MPa to 122 MPa with a median of 64 MPa. The lowest value was estimated for the 2004 M 6 Parkfield, California, earthquake whereas the nearby M 2.1 repeating earthquake, as recorded in the SAFOD pilot hole, showed a more typical yield stress of 64 MPa.

Bulletin of the Seismological Society of America↗

Improved constraints on the estimated size and volatile content of the Mount St. Helens magma system from the 2004–2008 history of dome growth and deformation

The history of dome growth and geodetic deflation during the 2004–2008 Mount St. Helens eruption can be fit to theoretical curves with parameters such as reservoir volume, bubble content, initial overpressure, and magma rheology, here assumed to be Newtonian viscous, with or without a solid plug in the conduit center. Data from 2004–2008 are consistent with eruption from a 10–25 km 3 reservoir containing 0.5–2% bubbles, an initial overpressure of 10–20 MPa, and no significant, sustained recharge. During the eruption we used curve fits to project the eruption's final duration and volume. Early projections predicted a final volume only about half of the actual value; but projections increased with each measurement, implying a temporal increase in reservoir volume or compressibility. A simple interpretation is that early effusion was driven by a 5–10 km 3 , integrated core of fluid magma. This core expanded with time through creep of semi‐solid magma and host rock.

Washington↗

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth↗