skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Linking 3D Long‐Term Slow‐Slip Cycle Models With Rupture Dynamics: The Nucleation of the 2014 M w 7.3 Guerrero, Mexico Earthquake
Abstract Slow slip events (SSEs) have been observed in spatial and temporal proximity to megathrust earthquakes in various subduction zones, including the 2014Mw7.3 Guerrero, Mexico earthquake which was preceded by aMw7.6 SSE. However, the underlying physics connecting SSEs to earthquakes remains elusive. Here, we link 3D slow‐slip cycle models with dynamic rupture simulations across the geometrically complex flat‐slab Cocos plate boundary. Our physics‐based models reproduce key regional geodetic and teleseismic fault slip observations on timescales from decades to seconds. We find that accelerating SSE fronts transiently increase shear stress at the down‐dip end of the seismogenic zone, modulated by the complex geometry beneath the Guerrero segment. The shear stresses cast by the migrating fronts of the 2014Mw7.6 SSE are significantly larger than those during the three previous episodic SSEs that occurred along the same portion of the megathrust. We show that the SSE transient stresses are large enough to nucleate earthquake dynamic rupture and affect rupture dynamics. However, additional frictional asperities in the seismogenic part of the megathrust are required to explain the observed complexities in the coseismic energy release and static surface displacements of the Guerrero earthquake. We conclude that it is crucial to jointly analyze the long‐ and short‐term interactions and complexities of SSEs and megathrust earthquakes across several (a)seismic cycles accounting for megathrust geometry. Our study has important implications for identifying earthquake precursors and understanding the link between transient and sudden megathrust faulting processes.  more » « less
Award ID(s):
2121568 2311206 2311208
PAR ID:
10535841
Author(s) / Creator(s):
;
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
AGU Advances
Volume:
5
Issue:
2
ISSN:
2576-604X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Shallow slow-slip events (SSEs) contribute to strain release near the shallow portions of subduction interfaces and may contribute to promoting shallow subduction earthquakes. Recent efforts in offshore monitoring of shallow SSEs have provided evidence of possible interactions between shallow SSEs and megathrust earthquakes. In this study, we use a dynamic earthquake simulator that captures both quasi-static (for SSEs) and dynamic (for megathrust earthquakes) slip to explore their interactions and implications for seismic and tsunami hazards. We model slip behaviors of a shallow-dipping subduction interface on which two locally locked patches (asperities) with different strengths are embedded within a conditionally stable zone. We find that both SSEs and earthquakes can occur, and they interact over multiple earthquake cycles in the model. Dynamic ruptures can nucleate on the asperities and propagate into the surrounding conditionally stable zone at slow speeds, generating tsunami earthquakes. A clear correlation emerges between the size of an earthquake and SSE activities preceding it. Small earthquakes rupture only the low-strength asperity, whereas large earthquakes rupture both. Before a large earthquake, periodic SSEs occur around the high-strength asperity, gradually loading stress into its interior. The critically stressed high-strength asperity can be ruptured together with the low-strength one in the large earthquake, followed by a relatively quiet interseismic period with very few SSEs and then a small earthquake. An SSE may or may not directly lead to nucleation of an earthquake, depending on whether a nearby asperity is ready for spontaneously dynamic failure. In addition, because of different SSE activities, the coupling degree may change dramatically between different interseismic periods, suggesting that its estimate based on a short period of observation may be biased. 
    more » « less
  2. 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
  3. Abstract On 18 November 2022, a large earthquake struck offshore southern Sumatra, generating a tsunami with 25 cm peak amplitude recorded at tide gauge station SBLT. OurW‐phase solution indicates a shallow dip of 6.2°, compatible with long‐period surface wave radiation patterns. Inversion of teleseismic body waves indicates a shallow slip distribution extending from about 10 km deep to near the trench with maximum slip of ∼4.1 m and seismic moment of  Nm (MW7.3). Joint modeling of seismic and tsunami data indicates a shallow rigidity of ∼23 GPa. We find a low moment‐scaled radiated energy of , similar to that of the 2010MW7.8 Mentawai event () and other tsunami earthquakes. These characteristics indicate that the 2022 event should be designated as a smaller moment magnitude tsunami earthquake compared to the other 12 well‐documented global occurrences since 1896. The 2022 event ruptured up‐dip of the 2007MW8.4 Bengkulu earthquake, demonstrating shallow seismogenic capability of a megathrust that had experienced both a deeper seismic event and adjacent shallow aseismic afterslip. We consider seismogenic behavior of shallow megathrusts and concern for future tsunami earthquakes in subduction zones globally, noting a correlation between tsunami earthquake occurrence and subducting seafloor covered with siliceous pelagic sediments. We suggest that the combination of pelagic clay and siliceous sediments and rough seafloor topography near the trench play important roles in controlling the genesis of tsunami earthquakes along Sumatra and other regions, rather than the subduction tectonic framework of accretionary or erosive margin. 
    more » « less
  4. 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
  5. null (Ed.)
    Abstract Either the triggering of large earthquakes on a fault hosting aseismic slip or the triggering of slow slip events (SSE) by passing seismic waves involve seismological questions with important hazard implications. Just a few observations plausibly suggest that such interactions actually happen in nature. In this study we show that three recent devastating earthquakes in Mexico are likely related to SSEs, describing a cascade of events interacting with each other on a regional scale via quasi-static and/or dynamic perturbations across the states of Guerrero and Oaxaca. Such interaction seems to be conditioned by the transient memory of Earth materials subject to the “traumatic” stress produced by seismic waves of the great 2017 (Mw8.2) Tehuantepec earthquake, which strongly disturbed the SSE cycles over a 650 km long segment of the subduction plate interface. Our results imply that seismic hazard in large populated areas is a short-term evolving function of seismotectonic processes that are often observable. 
    more » « less