Delayed Seismicity Rate Changes Controlled by Static Stress Transfer

On 15 June 2010, a M_w5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M_w7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity

Delayed Seismicity Rate Changes Controlled by Static Stress Transfer

On 15 June 2010, a M_w5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M_w7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity

A physics-based earthquake simulator replicates seismic hazard statistics across California

Seismic hazard models are important for society, feeding into building codes and hazard mitigation efforts. These models, however, rest on many uncertain assumptions and are difficult to test observationally because of the long recurrence times of large earthquakes. Physics-based earthquake simulators offer a potentially helpful tool, but they face a vast range of fundamental scientific uncertainties. We compare a physics-based earthquake simulator against the latest seismic hazard model for California. Using only uniform parameters in the simulator, we find strikingly good agreement of the long-term shaking hazard compared with the California model. This ability to replicate statistically based seismic hazard estimates by a physics-based model cross-validates standard methods and provides a new alternative approach needing fewer inputs and assumptions for estimating hazard

Earthquake Recurrence in Simulated Fault Systems

We employ a computationally efficient fault system earthquake simulator, RSQSim, to explore effects of earthquake nucleation and fault system geometry on earthquake occurrence. The simulations incorporate rate- and state-dependent friction, high-resolution representations of fault systems, and quasi-dynamic rupture propagation. Faults are represented as continuous planar surfaces, surfaces with a random fractal roughness, and discontinuous fractally segmented faults. Simulated earthquake catalogs have up to 106 earthquakes that span a magnitude range from ∼M4.5 to M8. The seismicity has strong temporal and spatial clustering in the form of foreshocks and aftershocks and occasional large-earthquake pairs. Fault system geometry plays the primary role in establishing the characteristics of stress evolution that control earthquake recurrence statistics. Empirical density distributions of earthquake recurrence times at a specific point on a fault depend strongly on magnitude and take a variety of complex forms that change with position within the fault system. Because fault system geometry is an observable that greatly impacts recurrence statistics, we propose using fault system earthquake simulators to define the empirical probability density distributions for use in regional assessments of earthquake probabilities

Complex characteristics of slow slip events in subduction zones reproduced in multi-cycle simulations

[1] Since the discovery of slow slip events along subduction zone interfaces worldwide, dense geodetic and seismic networks have illuminated detailed characteristics of these events and associated tremor. High‐resolution observations of tremor, where the spatial‐temporal evolution is presumed to reflect that of the underlying slow slip events, show highly complex patterns in which the origins remain poorly understood. We present a new, computationally efficient modeling technique that reproduces many features of observed slow slip events, including slow initiation, coalescence of separate events, and rapid back‐propagation of renewed slip over previously slipped regions. Rapid back propagation speeds are explained as a consequence of rate‐ and state‐dependent frictional healing, consistent with analytical solutions developed in support of the simulations