An official website of the United States government
Abstract Understanding the dynamics of microearthquakes is a timely challenge with the potential to address current paradoxes in earthquake mechanics, and to better understand earthquake ruptures induced by fluid injection. We perform fully 3D dynamic rupture simulations caused by fluid injection on a target fault for Fault Activation and Earthquake Ruptures experiments generatingMw ≤ 1 earthquakes. We investigate the dynamics of rupture propagation with spatially variable stress drop caused by pore pressure changes and assuming different slip‐weakening constitutive parameters. We show that the spontaneous arrest of propagating ruptures is possible by assuming a high fault strength parameter S, that is, a high ratio between strength excess and dynamic stress drop. In faults with high S values (low rupturing potential), even minor variations inDc(from 0.45 to 0.6 mm) have a substantial effect on the rupture propagation and the ultimate earthquake size. Modest spatial variations of dynamic stress drop determine the rupture mode, distinguishing self‐arresting from run‐away ruptures. Our results suggest that several characteristics inferred for accelerating dynamic ruptures differ from those observed during rupture deceleration of a self‐arresting earthquake. During deceleration, a decrease of peak slip velocity is associated with a nearly constant cohesive zone size. Moreover, the residual slip velocity value (asymptotic value for a crack‐like rupture) decreases to nearly zero. This means that an initially crack‐like rupture becomes a pulse‐like rupture during spontaneous arrest. These findings highlight the complex dynamics of small induced earthquakes, which differ from solutions obtained from conventional crack‐like models of earthquake rupture.
more »
« less
Quake-DFN: A Software for Simulating Sequences of Induced Earthquakes in a Discrete Fault Network
ABSTRACT We present an earthquake simulator, Quake-DFN, which allows simulating sequences of earthquakes in a 3D discrete fault network governed by rate and state friction. The simulator is quasi-dynamic, with inertial effects being approximated by radiation damping and a lumped mass. The lumped mass term allows for accounting for inertial overshoot and, in addition, makes the computation more effective. Quake-DFN is compared against three publicly available simulation results: (1) the rupture of a planar fault with uniform prestress (SEAS BP5-QD), (2) the propagation of a rupture across a stepover separating two parallel planar faults (RSQSim and FaultMod), and (3) a branch fault system with a secondary fault splaying from a main fault (FaultMod). Examples of injection-induced earthquake simulations are shown for three different fault geometries: (1) a planar fault with a wide range of initial stresses, (2) a branching fault system with varying fault angles and principal stress orientations, and (3) a fault network similar to the one that was activated during the 2011 Prague, Oklahoma, earthquake sequence. The simulations produce realistic earthquake sequences. The time and magnitude of the induced earthquakes observed in these simulations depend on the difference between the initial friction and the residual friction μi−μf, the value of which quantifies the potential for runaway ruptures (ruptures that can extend beyond the zone of stress perturbation due to the injection). The discrete fault simulations show that our simulator correctly accounts for the effect of fault geometry and regional stress tensor orientation and shape. These examples show that Quake-DFN can be used to simulate earthquake sequences and, most importantly, magnitudes, possibly induced or triggered by a fluid injection near a known fault system.
more »
« less
- Award ID(s):
- 1822214
- PAR ID:
- 10536544
- Publisher / Repository:
- Seismological Society of America
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- ISSN:
- 0037-1106
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Numerical and laboratory models of earthquake cycles on faults governed by rate‐and‐state friction often show cycle‐invariant behavior, while natural faults exhibit considerable variability in slip history. Possible explanations include heterogeneities in fault stress and frictional properties. We investigate how various types of heterogeneity in simulations of quasi‐dynamic sequences of seismic and aseismic slip affect rupture complexity, hypocenter location, and slow slip events (SSEs). We model a 2D vertical strike‐slip fault and study the roles of self‐affine fractal heterogeneities in normal stress, rate‐and‐state parameter , and characteristic slip‐weakening distance, as well as the effects of a low‐rigidity fault zone. We find that only a combination of heterogeneous parameters introduces variability in the modeled rupture extent, hypocenter depth, and recurrence interval. In particular, variable hypocenter depths require creeping patches within the velocity‐weakening seismogenic zone. A low‐rigidity fault zone adds little to slip complexity. Slip law simulations produce fewer partial ruptures, smaller stress drops, and lower peak slip rates compared to aging law simulations. We show that the ratio of the seismogenic zone thickness to nucleation size does not entirely predict slip complexity. The most complex aging law model, combining multiple heterogeneities, features system‐size earthquakes preceded by cascades of partial ruptures and spontaneous SSEs. For such models, a transition from aperiodic to quasi‐regular regimes requires more cycles than is typically needed to erase the effect of initial conditions. These results highlight the importance of heterogeneity in reproducing natural fault slip complexity in numerical models of earthquake sequences.more » « less
-
Abstract Earthquake stress drops are inferred to be independent of source depth, contradicting linear scaling predictions for earthquakes as frictional stick‐slip instabilities that assume increasing fault normal stress due to overburden. Here, we examine the scaling between averaged stress drops and increasing normal stress for simulated earthquake sequences in continuum rate‐and‐state fault models. The models produce a weaker dependence of stress drop on normal stress than the linearity of simple friction, which can be well‐fit by a sublinear power‐law. This result is more prominent when the fault dimension is much larger than nucleation scales. In such cases, the averaged behavior of ruptures is dominated by rupture propagation conditions, reflecting more heterogeneous shear stress conditions. As natural faults can be considerably larger than the smallest earthquakes they host, such a weaker scaling between averaged rupture conditions and normal stress may partially explain the lack of an inferred depth‐dependence of earthquake stress drops.more » « less
-
null (Ed.)SUMMARY Earthquake ruptures are generally considered to be cracks that propagate as fracture or frictional slip on pre-existing faults. Crack models have been used to describe the spatial distribution of fault offset and the associated static stress changes along a fault, and have implications for friction evolution and the underlying physics of rupture processes. However, field measurements that could help refine idealized crack models are rare. Here, we describe large-scale laboratory earthquake experiments, where all rupture processes were contained within a 3-m long saw-cut granite fault, and we propose an analytical crack model that fits our measurements. Similar to natural earthquakes, laboratory measurements show coseismic slip that gradually tapers near the rupture tips. Measured stress changes show roughly constant stress drop in the centre of the ruptured region, a maximum stress increase near the rupture tips and a smooth transition in between, in a region we describe as the earthquake arrest zone. The proposed model generalizes the widely used elliptical crack model by adding gradually tapered slip at the ends of the rupture. Different from the cohesive zone described by fracture mechanics, we propose that the transition in stress changes and the corresponding linear taper observed in the earthquake arrest zone are the result of rupture termination conditions primarily controlled by the initial stress distribution. It is the heterogeneous initial stress distribution that controls the arrest of laboratory earthquakes, and the features of static stress changes. We also performed dynamic rupture simulations that confirm how arrest conditions can affect slip taper and static stress changes. If applicable to larger natural earthquakes, this distinction between an earthquake arrest zone (that depends on stress conditions) and a cohesive zone (that depends primarily on strength evolution) has important implications for how seismic observations of earthquake fracture energy should be interpreted.more » « less
-
Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1785/0120230299
