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: Dynamic Emergence of Plate Motions and Great Megathrust Earthquakes Across Length and Time Scales
Abstract The slow motion of tectonic plates over thousands of kilometers is intermittently interrupted by great earthquakes with sudden slips localized near convergent plate boundaries. We developed a subduction model that self‐consistently integrates buoyancy forces, diffusion and dislocation creep, and inter‐plate friction. From the nonlinear dynamics emerge long‐term plate motions that achieve velocities of cm/year, effective viscosities of Pas below plates, and sudden slips up to m repeating every several hundred years. Along‐strike resistance arising from long‐wavelength variation of coseismic slip is naturally incorporated with a rupture length scale, . Computations with km generate events with . When decreases, there is a commensurate decrease in the effective moment of rupture events. Predicted long‐term plate velocities, mantle viscosities, cycles of stress loading and release, and rupture event size and magnitude all show good agreement with observations.  more » « less
Award ID(s):
2343864
PAR ID:
10556454
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
51
Issue:
22
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The importance of a low-viscosity asthenosphere underlying mobile plates has been highlighted since the earliest days of the plate tectonics revolution. However, absolute asthenospheric viscosities are still poorly constrained, with estimates spanning up to 3 orders of magnitude. Here we follow a new approach using analytic solutions for Poiseuille-Couette channel flow to compute asthenospheric viscosities under the Caribbean. We estimate Caribbean dynamic topography and the associated pressure gradient, which, combined with flow velocities estimated from geologic markers and tomographic structure, yield our best-estimate asthenospheric viscosity of (3.0 ± 1.5)*10 18 Pa s. This value is consistent with independent estimates for non-cratonic and oceanic regions, and challenges the hypothesis that higher-viscosity asthenosphere inferred from postglacial rebound is globally-representative. The active flow driven by Galapagos plume overpressure shown here contradicts the traditional view that the asthenosphere is only a passive lubricating layer for Earth’s tectonic plates. 
    more » « less
  2. In subduction zones, along‐strike and downdip variations in megathrust slip behavior are linked to changes in the properties of the subducting and overriding plates. Although marine geophysical methods provide insights into subduction zone structures, most surveys consist of sparse 2D profiles, limiting our understanding of first‐order controls. Here, we use active‐source seismic data to derive a 3D crustal‐scale P‐wave velocity model of the Alaska Peninsula subduction zone that encompasses both plates and spans the Semidi segment and SW Kodiak asperity. Our results reveal modest variations within the incoming plate, attributed to a series of fracture zones, seamounts and their associated basement swell, collectively contributing to plate hydration. Basement swell appears to modulate the distribution and type of sediment entering the trench, likely impacting observed variations in slip behavior. The overriding plate exhibits significant heterogeneity. The updip limit and width of the dynamic backstop are similar between the SW Kodiak asperity and eastern Semidi segment but differ significantly from the Western Semidi segment. These distinctions likely account for differences in earthquake rupture patterns and interseismic coupling among these segments. Additionally, high velocities in the mid‐lower forearc crust coincide with the location of the megathrust slip during the Mw 8.2 2021 Chignik event. We interpret these velocities as intracrustal intrusions that contributed to the deep rupture of the 2021 event. Our findings suggest that the contrasting structural and material properties of both the incoming and overriding plates influence the spatially complex and semi‐persistent segmentation of the megathrust offshore the Alaska Peninsula. 
    more » « less
  3. Abstract While many large earthquakes are preceded by observable foreshocks, the mechanism responsible for the occurrence of these smaller‐scale seismic events remains uncertain. One physical explanation of foreshocks with growing support is that they are produced by the interaction of slow slip with fault heterogeneity. Inspired by the suggestion from laboratory experiments that foreshocks occur on fault asperities (bumps), we explore rate‐and‐state fault models with patches of higher normal stress embedded in a larger seismogenic region by conducting 3‐D numerical simulations of their behavior over long‐term sequences of aseismic and seismic slips. The models do produce smaller‐scale seismicity during the aseismic nucleation of larger‐scale seismic events. These smaller‐scale events have reasonable stress drops, despite the highly elevated compression assigned to the source patches. We find that the two main factors contributing to the reasonable stress drops are the significant extent of the rupture into the region surrounding the patches and the aseismic stress release just prior to the seismic events. The smaller‐scale seismicity can only occur if a sufficient separation in nucleation scales between the foreshock‐like events and mainshocks is achieved. Our modeling provides insight into the conditions conducive for generating foreshocks on both natural and laboratory faults. 
    more » « less
  4. 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
  5. Abstract A devastating magnitude 7.2 earthquake struck Southern Haiti on 14 August 2021. The earthquake caused severe damage and over 2000 casualties. Resolving the earthquake rupture process can provide critical insights into hazard mitigation. Here we use integrated seismological analyses to obtain the rupture history of the 2021 earthquake. We find the earthquake first broke a blind thrust fault and then jumped to a disconnected strike‐slip fault. Neither of the fault configurations aligns with the left‐lateral tectonic boundary between the Caribbean and North American plates. The complex multi‐fault rupture may result from the oblique plate convergence in the region, so that the initial thrust rupture is due to the boundary‐normal compression and the following strike‐slip faulting originates from the Gonâve microplate block movement, orienting SW‐NE direction. The complex rupture development of the earthquake suggests that the regional deformation is accommodated by a network of segmented faults with diverse faulting conditions. 
    more » « less