An official website of the United States government
Abstract Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene M w > 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai’iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.
more »
« less
Subduction Zone Geometry Modulates the Megathrust Earthquake Cycle: Magnitude, Recurrence, and Variability
Abstract Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably‐dipping planar and variably‐curved nonplanar megathrusts using the volumetric, high‐order accurate codetandemto account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly‐dipping megathrusts to host larger earthquakes than steeply‐dipping ones. Under elevated pore pressure and less strongly velocity‐weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly‐dipping and sharply‐curved megathrusts host multi‐period supercycles of slow‐to‐fast, small‐to‐large slip events under higher effective stresses and more strongly velocity‐weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second‐order effects from structurally‐derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co‐ and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.
more »
« less
- PAR ID:
- 10535823
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 8
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Tsunamis generated by seafloor displacements accompanying large submarine earthquakes provide sensitivity to absolute slip position and distribution for offshore faulting analogous to that of geodetic observations for landward faulting. Tsunami recordings at deep‐water and near‐shore ocean bottom pressure sensors and tide gauges, along with runup and inundation measurements, can now be reliably modeled using detailed bathymetric structures and robust numerical codes. As a result, tsunami observations now play an important role in quantifying coseismic slip distributions for large submarine earthquakes in subduction zones and other tectonic environments. Applications of joint modeling or inversion of seismic, geodetic and tsunami observations for recent major earthquakes are described, highlighting the specific contributions of the tsunami observations to source model resolution. Tsunami observations provide unique information on the up‐dip extent of earthquake coseismic slip on subduction zone megathrust faults and occurrence of near‐trench slip, which are usually not well constrained by seismic and land‐based geodetic signals. Tsunami signals also help to detect offshore slow slip that is not evident in seismic or land‐based geodetic data and to balance geophysical constraints on ruptures that extend from on‐shore to off‐shore. Tsunami runup measurements and stratigraphic deposits further provide unique constraints on large earthquake ruptures that occurred prior to modern geophysical instrumentation.more » « less
-
Fluids influence fault zone strength and the occurrence of earthquakes, slow slip events, and aseismic slip. We introduce an earthquake sequence model with fault zone fluid transport, accounting for elastic, viscous, and plastic porosity evolution, with permeability having a power‐law dependence on porosity. Fluids, sourced at a constant rate below the seismogenic zone, ascend along the fault. While the modeling is done for a vertical strike‐slip fault with 2D antiplane shear deformation, the general behavior and processes are anticipated to apply also to subduction zones. The model produces large earthquakes in the seismogenic zone, whose recurrence interval is controlled in part by compaction‐driven pressurization and weakening. The model also produces a complex sequence of slow slip events (SSEs) beneath the seismogenic zone. The SSEs are initiated by compaction‐driven pressurization and weakening and stalled by dilatant suctions. Modeled SSE sequences include long‐term events lasting from a few months to years and very rapid short‐term events lasting for only a few days; slip is ∼1–10 cm. Despite ∼1–10 MPa pore pressure changes, porosity and permeability changes are small and hence fluid flux is relatively constant except in the immediate vicinity of slip fronts. This contrasts with alternative fault valving models that feature much larger changes in permeability from the evolution of pore connectivity. Our model demonstrates the important role that compaction and dilatancy have on fluid pressure and fault slip, with possible relevance to slow slip events in subduction zones and elsewhere.more » « less
-
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
-
Faults are usually surrounded by damage zones associated with localized deformation. Here we use fully dynamic earthquake cycle simulations to quantify the behaviors of earthquakes in fault damage zones. We show that fault damage zones can make a significant contribution to the spatial and temporal seismicity distribution. Fault stress heterogeneities generated by fault zone waves persist over multiple earthquake cycles that, in turn, produce small earthquakes that are absent in homogeneous simulations with the same friction conditions. Shallow fault zones can produce a bimodal depth distribution of earthquakes with clustering of seismicity at both shallower and deeper depths. Fault zone healing during the interseismic period also promotes the penetration of aseismic slip into the locked region and reduces the sizes of fault asperities that host earthquakes. Hence, small and moderate subsurface earthquakes with irregular recurrence intervals are commonly observed in immature fault zone simulations with interseismic healing. To link our simulation results to geological observations, we will use simulated fault slip at different depths to infer the timing and recurrence intervals of earthquakes and discuss how such measurements can affect our understanding of earthquake behaviors. We will also show that the maturity and material properties of fault damage zones have strong influence on whether long-term earthquake characteristics are represented by single events.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2024JB029191
