An official website of the United States government
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
Scale Dependence of Earthquake Rupture Prestress in Models With Enhanced Weakening: Implications for Event Statistics and Inferences of Fault Stress
Abstract Determining conditions for earthquake slip on faults is a key goal of fault mechanics highly relevant to seismic hazard. Previous studies have demonstrated that enhanced dynamic weakening (EDW) can lead to dynamic rupture of faults with much lower shear stress than required for rupture nucleation. We study the stress conditions before earthquake ruptures of different sizes that spontaneously evolve in numerical simulations of earthquake sequences on rate‐and‐state faults with EDW due to thermal pressurization of pore fluids. We find that average shear stress right before dynamic rupture (aka shear prestress) systematically varies with the rupture size. The smallest ruptures have prestress comparable to the local shear stress required for nucleation. Larger ruptures weaken the fault more, propagate over increasingly under‐stressed areas due to dynamic stress concentration, and result in progressively lower average prestress over the entire rupture. The effect is more significant in fault models with more efficient EDW. We find that, as a result, fault models with more efficient weakening produce fewer small events and result in systematically lower b‐values of the frequency‐magnitude event distributions. The findings (a) illustrate that large earthquakes can occur on faults that appear not to be critically stressed compared to stresses required for slip nucleation; (b) highlight the importance of finite‐fault modeling in relating the local friction behavior determined in the lab to the field scale; and (c) suggest that paucity of small events or seismic quiescence may be the observational indication of mature faults that operate under low shear stress due to EDW.
more »
« less
- Award ID(s):
- 1724686
- PAR ID:
- 10371826
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 10
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Mediterranean Hellenic Arc subduction zone (HASZ) has generated several 8 earthquakes and tsunamis. Seismic‐probabilistic tsunami hazard assessment typically utilizes uniform or stochastic earthquake models, which may not represent dynamic rupture and tsunami generation complexity. We present an ensemble of ten 3D dynamic rupture earthquake scenarios for the HASZ, utilizing a realistic slab geometry. Our simplest models use uniform along‐arc pre‐stresses or a single circular initial stress asperity. We then introduce progressively more complex models varying initial shear stress along‐arc, multiple asperities based on scale‐dependent critical slip weakening distance, and a most complex model blending all aforementioned heterogeneities. Thereby, regional initial conditions are constrained without relying on detailed geodetic locking models. Varying epicentral locations in the simplest, homogeneous model leads to different rupture speeds and moment magnitudes. We observe dynamic fault slip penetrating the shallow slip‐strengthening region and affecting seafloor uplift. Off‐fault plastic deformation can double vertical seafloor uplift. A single‐asperity model generates a 8 scenario resembling the 1303 Crete earthquake. Using along‐strike varying initial stresses results in 8.0–8.5 dynamic rupture scenarios with diverse slip rates and uplift patterns. The model with the most heterogeneous initial conditions yields a 7.5 scenario. Dynamic rupture complexity in prestress and fracture energy tends to lower earthquake magnitude but enhances tsunamigenic displacements. Our results offer insights into the dynamics of potential large Hellenic Arc megathrust earthquakes and may inform future physics‐based joint seismic and tsunami hazard assessments.more » « less
-
Abstract Finite-element models of neotectonics require transform faults to rupture seismically even where preseismic shear stresses are low, presumably by dynamic-weakening mechanisms. A long-standing objection is that, if a rupture initiated at an asperity with high static friction stresses, which then transitioned to low dynamic-weakening stresses, local stress drop would be near total and on the order of 80 MPa, which is 4×–40× greater than observed. But the 5 Mw ≥ 7.8 transform earthquakes since 2000 initially ruptured on the branch faults of small net slip (Stein and Bird, 2024). If the slip initiates on a branch fault with different slip physics and no dynamic weakening, this solves the stress-drop problem. We propose that most large shallow earthquakes are hybrid ruptures, which begin on branch faults of small slip with high shear stresses, and then continue propagating on a connected dynamically weakened fault of large slip, even where shear stresses are low. One prediction of this model is that most large shallow ruptures should be unilateral. We test this prediction against the 100 largest (m ≥ 6.49) shallow continental strike-slip earthquakes 1977–2022, using information from the Global Centroid Moment Tensor and International Seismological Centre catalogs. The differences in time and location between the epicenter and the epicentroid define a horizontal “migration” velocity vector for the evolving centroid of each rupture. Early aftershock locations are summarized by a five-parameter elliptical model. Using the geometric relations between these (and mapped traces of active faults) and guided by a symmetrical decision table, we classified 55 ruptures as apparently unilateral, 30 as bilateral, and 15 as ambiguous. Our finding that a majority (55%–70%) of these ruptures are unilateral permits the interpretation that a majority of ruptures are hybrids, both in terms of geometry (branch fault to transform) and in terms of the physics of their fault slip.more » « less
-
Abstract Fault regions inferred to be slowly slipping are interpreted to accommodate much of tectonic plate motion aseismically and potentially serve as barriers to earthquake rupture. Here, we build on prior work using simulations of earthquake sequences with enhanced dynamic fault weakening to show how fault regions that exhibit decades of steady creep or transient slow‐slip events can be driven to dynamically fail by incoming earthquake ruptures. Following substantial earthquake slip, such regions can be under‐stressed and locked for centuries prior to slowly slipping again. Our simulations illustrate that slow fault slip indicates that a region is sufficiently loaded to be failing about its quasi‐static strength. Hence, if a fault region is susceptible to failing dynamically, then observations of slow slip could serve as an indication that the region is critically stressed and ready to fail in a future earthquake, posing a qualitatively different interpretation of slow slip for seismic hazard.more » « less
-
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
