Previous geodetic and teleseismic observations of the 2021 Mw 7.4 Maduo earthquake imply surprising but difficult-to-constrain complexity, including rupture across multiple fault segments and supershear rupture. Here, we present an integrated analysis of multi-fault 3D dynamic rupture models, high-resolution optical correlation analysis, and joint optical-InSAR slip inversion. Our preferred model, validated by the teleseismic multi-peak moment rate release, includes unilateral eastward double-onset supershear speeds and cascading rupture dynamically triggering two adjacent fault branches. We propose that pronounced along-strike variation in fracture energy, complex fault geometries, and multi-scale variable prestress drives this event's complex rupture dynamics. We illustrate how supershear transition has signatures in modeled and observed off-fault deformation. Our study opens new avenues to combine observations and models to better understand complex earthquake dynamics, including local and potentially repeating supershear episodes across immature faults or under heterogeneous stress and strength conditions, which are potentially not unusual.
more »
« less
3D Dynamic Rupture Modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 Earthquake Doublet Using Early Observations
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
- Award ID(s):
- 2121568
- PAR ID:
- 10535834
- Publisher / Repository:
- SSA
- Date Published:
- Journal Name:
- The Seismic Record
- Volume:
- 3
- Issue:
- 4
- ISSN:
- 2694-4006
- Page Range / eLocation ID:
- 342 to 356
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The destructive 2023 moment magnitude ( M w ) 7.8-7.7 earthquake doublet ruptured multiple segments of the East Anatolian Fault system in Turkey. We integrate multi-scale seismic and space-geodetic observations with multi-fault kinematic inversions and dynamic rupture modeling to unravel the events’ complex rupture history and stress-mediated fault interactions. Our analysis reveals three sub-shear slip episodes during the initial M w 7.8 earthquake with delayed rupture initiation to the southwest. The M w 7.7 event occurred 9 hours later with larger slip and supershear rupture on its western branch. Mechanically consistent dynamic models accounting for fault interactions can explain the unexpected rupture paths, and require a heterogeneous background stress. Our results highlight the importance of combining near- and far-field observations with data-driven and physics-based models for seismic hazard assessment.more » « less
-
Abstract Previous geodetic and teleseismic observations of the 2021 7.4 Maduo earthquake imply surprising but difficult‐to‐constrain complexity, including rupture across multiple fault segments and supershear rupture. Here, we present an integrated analysis of multi‐fault 3D dynamic rupture models, high‐resolution optical correlation analysis, and joint optical‐InSAR slip inversion. Our preferred model, validated by the teleseismic multi‐peak moment rate release, includes unilateral eastward double‐onset supershear speeds and cascading rupture dynamically triggering two adjacent fault branches. We propose that pronounced along‐strike variation in fracture energy, complex fault geometries, and multi‐scale variable prestress drives this event's complex rupture dynamics. We illustrate how supershear transition has signatures in modeled and observed off‐fault deformation. Our study opens new avenues to combine observations and models to better understand complex earthquake dynamics, including local and potentially repeating supershear episodes across immature faults or under heterogeneous stress and strength conditions, which are potentially not unusual.more » « less
-
ABSTRACT Fault zones exhibit geometrical complexity and are often surrounded by multiscale fracture networks within their damage zones, potentially influencing rupture dynamics and near-field ground motions. In this study, we investigate the ground-motion characteristics of cascading ruptures across damage zone fracture networks of moderate-size earthquakes (Mw 5.5–6.0) using high-resolution 3D dynamic rupture simulations. Our models feature a listric normal fault surrounded by more than 800 fractures, emulating a major fault and its associated damage zone. We analyze three cases: a cascading rupture propagating within the fracture network (Mw 5.5), a non-cascading main-fault rupture with off-fault fracture slip (Mw 6.0), and a main-fault rupture without a fracture network (Mw 6.0). Cascading ruptures within the fracture network produce distinct ground-motion signatures with enriched high-frequency content, arising from simultaneous slip of multiple fractures and parts of the main fault, resembling source coda-wave-like signatures. This case shows elevated near-field characteristic frequency (fc) and stress drop, approximately an order of magnitude higher than the estimation directly on the fault of the dynamic rupture simulation. The inferred fc of the modeled vertical ground-motion components reflects the complexity of the radiation pattern and rupture directivity of fracture-network cascading earthquakes. We show that this is consistent with observations of strong azimuthal dependence of corner frequency in the 2009–2016 central Apennines, Italy, earthquake, sequence. Simulated ground motions from fracture-network cascading ruptures also show pronounced azimuthal variations in peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration, with average PGA nearly double that of the non-cascading cases. Cascading ruptures radiate high-frequency seismic energy, yield nontypical ground-motion characteristics including coda-wave-like signatures, and may result in a significantly higher seismologically inferred stress drop and PGA. Such outcomes emphasize the critical role of fault-zone complexity in affecting rupture dynamics and seismic radiation and have important implications for physics-based seismic hazard assessment.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
An official website of the United States government

