Abstract. Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need to reexamine the assumption of predetermined frictional breakdown energy common in dynamic rupture modeling and to better understand the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes.
more »
« less
The earthquake arrest zone
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
- Award ID(s):
- 1763499
- NSF-PAR ID:
- 10218729
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 224
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 581 to 589
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Understanding the dynamics of microearthquakes is a timely challengewith 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 fluidinjection on a target fault for FEAR experiments generating Mw ≤ 1earthquakes. We investigate the dynamics of rupture propagation withspatially variable stress drop caused by pore pressure changes andassuming different constitutive parameters. We show that the spontaneousarrest of propagating ruptures is possible by assuming a high faultstrength parameter S, that is, a high ratio between strength excess anddynamic stress drop. In faults with high S values (low rupturingpotential), even minor variations in Dc (from 0.45 to 0.6 mm) have asubstantial effect on the rupture propagation and the ultimateearthquake size. Our results show that modest spatial variations ofdynamic stress drop determine the rupture mode, distinguishingself-arresting from run-away ruptures. Our results suggest that severalcharacteristics inferred for accelerating dynamic ruptures differ fromthose observed during rupture deceleration of a self-arrestingearthquake. During deceleration, a decrease of peak slip velocity isassociated with a nearly constant cohesive zone size. Moreover, theresidual slip velocity value (asymptotic value for a crack-like rupture)decreases to nearly zero. This means that an initially crack-likerupture becomes a pulse-like rupture during spontaneous arrest. Insummary, our findings highlight the complex dynamics of smallearthquakes, which are partially contrasting with established crack-likemodels of earthquake rupture.
-
more » « less
Abstract To better quantify how injection, prior seismicity, and fault properties control rupture growth and propagation of induced earthquakes, we perform a finite‐fault slip inversion on a
M w4.0 earthquake that occurred in April 2015, the largest earthquake in an induced sequence near Guthrie, Oklahoma. The slip inversion reveals a rupture with slip patches that are anti‐correlated to the locations of prior seismicity. The prior seismicity driven by low pore pressure changes and static stress changes occurred on weaker portions of the fault, while theM w4.0 earthquake likely ruptured relatively stronger portions of the fault. To resolve if pore pressure changes or the initial underlying stress distribution and fault strength controlled the final slip distribution of the GuthrieM w4.0 earthquake, we compare strike‐slip events of similar magnitude from tectonically active regions and previously inactive regions. Earthquakes on reactivated faults exhibit different slip distributions than active regions, they have more prominent and well separated slip patches, a behavior often associated with faults of lower fault maturity. Pore pressure shows little effect on the distributions. These observations suggest that the initial underlying stress distribution and fault strength of reactivated faults in low deformation regions is the primary controlling factor of the slip distribution with pore pressure perturbations and earthquake interactions being secondary. Therefore, GuthrieM w4.0 earthquakes slip distribution was enhanced by pore‐pressure perturbations and earthquake interactions by creating an optimal stress state for its failure, but the slip distribution itself is controlled by its fault's initial stress and strength state. -
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.
-
null (Ed.)The fault damage zone is a well-known structure of localized deformation around faults. Its material properties evolve over earthquake cycles due to coseismic damage accumulation and interseismic healing. We will present fully dynamic earthquake cycle simulations to show how the styles of earthquake nucleation and rupture propagation change as fault zone material properties vary temporally. First, we will focus on the influence of fault zone structural maturity quantified by near-fault seismic wave velocities in simulations. The simulations show that immature fault zones promote small and moderate subsurface earthquakes with irregular recurrence intervals, whereas mature fault zones host pulse-like earthquake rupture that can propagate to the surface, extend throughout the seismogenic zone, and occur at regular intervals. The interseismic healing in immature fault zones plays a key role in allowing the development of aseismic slip episodes including slow-slip events and creep, which can propagate into the seismogenic zone, and thus limit the sizes of subsequent earthquakes by releasing fault stress. In the second part, we will discuss how the precursory changes of seismic wave velocities of fault damage zones may affect earthquake nucleation process. Both laboratory experiments and seismic observations show that the abrupt earthquake failure can be preceded by accelerated fault deformation and the accompanying velocity reduction of near-fault rocks. We will use earthquake cycle simulations to systematically test the effects of timing and amplitudes of such precursory velocity changes. Our simulations will provide new insights into the interplay between fault zone structure and earthquake nucleation process, which can be used to guide future real-time monitoring of major fault zones.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggaa386
