An official website of the United States government
Abstract Megathrust earthquakes exhibit a ubiquitous seismic radiation style: low‐frequency (LF) seismic energy is efficiently emitted from the shallowest portion of the fault, whereas high‐frequency (HF) seismic energy is efficiently emitted from the deepest part of the fault. Although this is observed in many case‐specific studies, we show that it is ubiquitous in global megathrust earthquakes between 1995 and 2021. Previous studies have interpreted this as an effect of systematic depth variation in either the plate interface frictional properties (Lay et al., 2012) or the P wavespeeds (Sallarès & Ranero, 2019). This work suggests an alternative hypothesis: the interaction between waves and ruptures due to the Earth's free surface is the leading mechanism that generates this behavior. Two‐dimensional dynamic rupture simulations of subduction zone earthquakes support this hypothesis. Our simulations show that the interaction between the seismic waves reflected at the Earth's free surface and the updip propagating rupture results in LF radiation at the source. In contrast, the downdip propagation of rupture is less affected by the free surface and is thus dominated by HF radiation typical of buried faults. To a second degree, the presence of a realistic Earth structure derived from P‐wave velocity (VP) tomographic images and realistic VP/VSratio estimated in boreholes further enhances the contrast in source radiation. We conclude that the Earth's free surface is necessary to explain the observed megathrust earthquake radiation style, and the realistic structure of subduction zone is necessary to better predict earthquake ground motion and tsunami potential.
more »
« less
Three-dimensional wave propagation and earthquake dynamic rupture simulations in complex poroelastic media
SUMMARY Numerical simulations of earthquakes and seismic wave propagation require accurate material models of the solid Earth. In contrast to purely elastic rheology, poroelasticity accounts for pore fluid pressure and fluid flow in porous media. Poroelastic effects can alter both the seismic wave field and the dynamic rupture characteristics of earthquakes. For example, the presence of fluids may affect cascading multifault ruptures, potentially leading to larger-than-expected earthquakes. However, incorporating poroelastic coupling into the elastodynamic wave equations increases the computational complexity of numerical simulations compared to elastic or viscoelastic material models, as the underlying partial differential equations become stiff. In this study, we use a Discontinuous Galerkin solver with Arbitrary High-Order DERivative time stepping of the poroelastic wave equations implemented in the open-source software SeisSol to simulate 3-D complex seismic wave propagation and 3-D dynamic rupture in poroelastic media. We verify our approach for double-couple point sources using independent methods including a semi-analytical solution and a finite-difference scheme and a homogeneous full-space and a poroelastic layer-over-half-space model, respectively. In a realistic carbon capture and storage reservoir scenario at the Sleipner site in the Utsira Formation, Norway, we model 3-D wave propagation through poroelastic sandstone layers separated by impermeable shale. Our results show a sudden change in the pressure field across material interfaces, which manifests as a discontinuity when viewed at the length scale of the dominant wavelengths of S or fast P waves. Accurately resolving the resulting steep pressure gradient dramatically increases the computational demands, requiring high-resolution modelling. We show that the Gassmann elastic equivalent model yields almost identical results to the fully poroelastic model when focusing solely on solid particle velocities. We extend this approach using suitable numerical fluxes to 3-D dynamic rupture simulations in complex fault systems, presenting the first 3-D scenarios that combine poroelastic media with geometrically complex, multifault rupture dynamics and tetrahedral meshes. Our findings reveal that, in contrast to modelling wave propagation only, poroelastic materials significantly alter rupture characteristics compared to using elastic equivalent media since the elastic equivalent fails to capture the evolution of pore pressure. Particularly in fault branching scenarios, the Biot coefficient plays a key role in either promoting or inhibiting fault activation. In some cases, ruptures are diverted to secondary faults, while in others, poroelastic effects induce rupture arrest. In a fault zone dynamic rupture model, we find poroelasticity aiding pulse-like rupture. A healing front is induced by the reduced pore pressure due to reflected waves from the boundaries of the poroelastic damage zone. Our results highlight that poroelastic effects are important for realistic simulations of seismic waves and earthquake rupture dynamics. In particular, our poroelastic simulations may offer new insights on the complexity of multifault rupture dynamics, fault-to-fault interaction and seismic wave propagation in realistic models of the Earth’s subsurface.
more »
« less
- PAR ID:
- 10613469
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 242
- Issue:
- 2
- ISSN:
- 0956-540X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Detailed imaging of accretionary wedges reveals splay fault networks that could pose a significant tsunami hazard. However, the dynamics of multiple splay fault activation during megathrust earthquakes and the consequent effects on tsunami generation are not well understood. We use a 2‐D dynamic rupture model with complex topo‐bathymetry and six curved splay fault geometries constrained from realistic tectonic loading modeled by a geodynamic seismic cycle model with consistent initial stress and strength conditions. We find that all splay faults rupture coseismically. While the largest splay fault slips due to a complex rupture branching process from the megathrust, all other splay faults are activated either top down or bottom up by dynamic stress transfer induced by trapped seismic waves. We ascribe these differences to local non‐optimal fault orientations and variable along‐dip strength excess. Generally, rupture on splay faults is facilitated by their favorable stress orientations and low strength excess as a result of high pore‐fluid pressures. The ensuing tsunami modeled with non‐linear 1‐D shallow water equations consists of one high‐amplitude crest related to rupture on the longest splay fault and a second broader wave packet resulting from slip on the other faults. This results in two episodes of flooding and a larger run‐up distance than the single long‐wavelength (300 km) tsunami sourced by the megathrust‐only rupture. Since splay fault activation is determined by both variable stress and strength conditions and dynamic activation, considering both tectonic and earthquake processes is relevant for understanding tsunamigenesis.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
-
The seismic potential of active low-angle normal faults (LANFs, <30° dip) remains enigmatic under Andersonian faulting theory, which predicts that normal faults dipping less than 30° should be inactive. The Alto Tiberina fault (ATF) in the northern Apennines, a partly creeping 17°-dipping LANF, has not been associated with any historical earthquakes but could potentially generate earthquakes up to Mw~7. We investigate the mechanical preconditions and dynamic plausibility of large ATF earthquakes using 3D dynamic rupture and seismic wave propagation simulations constrained by multidisciplinary data from the Alto Tiberina Near Fault Observatory (TABOO-NFO). Our models incorporate the complex non-planar ATF fault geometry, including hanging wall secondary faults and a recent geodetic coupling model. We show that potential large earthquakes (up to Mw~7.4) are mechanically viable under Andersonian extensional stress conditions if the ATF is statically relatively weak (μs=0.37). Large earthquakes only nucleate on favorably oriented, steeper fault sections (dip ≥30°), and remain confined to the coupled portion, limiting earthquake magnitude. These ruptures may dynamically trigger an intersecting synthetic branch but are unlikely to affect more distant antithetic faults. Jointly integrating fault geometry and geodetic coupling is crucial for forecasting dynamic rupture nucleation and propagation.more » « less
