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.


This content will become publicly available on April 2, 2026

Title: The 2021 M W 8.1 Kermadec Earthquake Sequence: Great Earthquake Rupture Along the Mantle/Slab Contact
Abstract Most great earthquakes on subduction zone plate boundaries have large coseismic slip concentrated along the contact between the subducting slab and the upper plate crust. On 4 March 2021, a magnitude 7.4 foreshock struck 1 hr 47 min before a magnitude 8.1 earthquake along the northern Kermadec island arc. The mainshock is the largest well‐documented underthrusting event along the ∼2,500‐km long Tonga‐Kermadec subduction zone. Using teleseismic, geodetic, and tsunami data, we find that all substantial coseismic slip in the mainshock is located along the mantle/slab interface at depths from 20 to 55 km, with the large foreshock nucleating near the down‐dip edge. Smaller foreshocks and most aftershocks are located up‐dip of the mainshock, where substantial prior moderate thrust earthquake activity had occurred. The upper plate crust is ∼17 km thick in northern Kermadec with only moderate‐size events along the crust/slab interface. A 1976 sequence withMWvalues of 7.9, 7.8, 7.3, 7.0, and 7.0 that spanned the 2021 rupture zone also involved deep megathrust rupture along the mantle/slab contact, but distinct waveforms exclude repeating ruptures. Variable waveforms for eight deep M6.9+ thrusting earthquakes since 1990 suggest discrete slip patches distributed throughout the region. The ∼300‐km long plate boundary in northern Kermadec is the only documented subduction zone region where the largest modeled interplate earthquakes have ruptured along the mantle/slab interface, suggesting that local frictional properties of the putatively hydrated mantle wedge may involve a dense distribution of Antigorite‐rich patches with high slip rate velocity weakening behavior in this locale.  more » « less
Award ID(s):
1802364
PAR ID:
10581148
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
130
Issue:
4
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract On 4 and 6 July 2019, two large strike‐slip earthquakes withW‐phase moment magnitudesMWW6.5 (foreshock) andMWW7.1 (mainshock) struck the Eastern California Shear Zone, northeast of Ridgecrest. The faulting geometry and kinematic coseismic slip distribution of both events are determined by jointly inverting seismological and geodetic observations guided by aftershock and surface rupture locations. The foreshock ruptured two orthogonal faults with a prominent L‐shaped geometry with maximum slip of ~1.1 m on the NE‐SW segment. The mainshock faulting extended NW‐SE along several primary fault segments that straddle the foreshock slip. The surface rupture and slip model indicate mostly near‐horizontal strike‐slip motion with maximum slip of ~3.7 m, but there is a localized vertical dip‐slip motion. Both the foreshock and mainshock ruptures terminate in regions of complex surface offsets. High aftershock productivity and low rupture velocity may be the result of rupture of a relatively immature fault system. 
    more » « less
  2. Abstract On 19 September 2022, a major earthquake struck the northwestern Michoacán segment along the Mexican subduction zone. A slip model is obtained that satisfactorily explains geodetic, teleseismic, and tsunami observations of the 2022 event. The preferred model has a compact large-slip patch that extends up-dip and northwestward from the hypocenter and directly overlaps a 1973 Mw 7.6 rupture. Slip is concentrated offshore and below the coast at depths from 10 to 30 km with a peak value of ∼2.9 m, and there is no detected coseismic slip near the trench. The total seismic moment is 3.1×1020  N·m (Mw 7.6), 72% of which is concentrated in the first 30 s. Most aftershocks are distributed in an up-dip area of the mainshock that has small coseismic slip, suggesting near-complete strain release in the large-slip patch. Teleseismic P waveforms of the 2022 and 1973 earthquakes are similar in duration and complexity with high cross-correlation coefficients of 0.68–0.98 for long P to PP signal time windows, indicating that the 2022 earthquake is a quasi-repeat of the 1973 earthquake, possibly indicating persistent frictional properties. Both the events produced more complex P waveforms than comparable size events along Guerrero and Oaxaca, reflecting differences in patchy locking of the Mexican megathrust. 
    more » « less
  3. Abstract The 2021MW6.0 Yangbi, Yunnan strike‐slip earthquake occurred on an unmapped crustal fault near the Weixi‐Qiaoho‐Weishan Fault along the southeast margin of the Tibetan Plateau. Using near‐source broadband seismic data from ChinArray, we investigate the spatial and temporal rupture evolution of the mainshock using apparent moment‐rate functions (AMRFs) determined by the empirical Green's function (EGF) method. Assuming a 1D line source on the fault plane, the rupture propagated unilaterally southeastward (∼144°) over a rupture length of ∼8.0 km with an estimated rupture speed of 2.1 km/s to 2.4 km/s. A 2D coseismic slip distribution for an assumed maximum rupture propagation speed of 2.2 km/s indicates that the rupture propagated to the southeast ∼8.0 km along strike and ∼5.0 km downdip with a peak slip of ∼2.1 m before stopping near the largest foreshock, where three bifurcating subfaults intersect. Using the AMRFs, the radiated energy of the mainshock is estimated as ∼. The relatively low moment scaled radiated energyof 1.5 × 10−5and intense foreshock and aftershock activity might indicate reactivation of an immature fault. The earthquake sequence is mainly distributed along a northwest‐southeast trend, and aftershocks and foreshocks are distributed near the periphery of the mainshock large‐slip area, suggesting that the stress in the mainshock slip zone is significantly reduced to below the level for more than a few overlapping aftershock to occur. 
    more » « less
  4. Abstract A great earthquake struck the Semidi segment of the plate boundary along the Alaska Peninsula on 29 July 2021, re‐rupturing part of the 1938 rupture zone. The 2021MW8.2 Chignik earthquake occurred just northeast of the 22 July 2020MW7.8 Simeonof earthquake, with little slip overlap. Analysis of teleseismicPandSHwaves, regional Global Navigation Satellite System (GNSS) displacements, and near‐field and far‐field tsunami observations provides a good resolution of the 2021 rupture process. During ∼60‐s long faulting, the slip was nonuniformly distributed along the megathrust over depths from 32 to 40 km, with up to ∼12.9‐m slip in an ∼170‐km‐long patch. The 40–45 km down‐dip limit of slip is well constrained by GNSS observations along the Alaska Peninsula. Tsunami observations preclude significant slip from extending to depths <25 km, confining all coseismic slip to beneath the shallow continental shelf. Most aftershocks locate seaward of the large‐slip zones, with a concentration of activity up‐dip of the deeper southwestern slip zone. Some localized aftershock patches locate beneath the continental slope. The surface‐wave magnitudeMSof 8.1 for the 2021 earthquake is smaller thanMS = 8.3–8.4 for the 1938 event. Seismic and tsunami data indicate that slip in 1938 was concentrated in the eastern region of its aftershock zone, extending beyond the Semidi Islands, where the 2021 event did not rupture. 
    more » « less
  5. Abstract Subduction zones host some of Earth's most damaging natural hazards, including megathrust earthquakes and earthquake‐induced tsunamis. A major control on the initiation and rupture characteristics of subduction megathrust earthquakes is how the coupled zone along the subduction interface accumulates elastic strain between events. We present results from observations of slow slip events (SSEs) in Cascadia occurring during the interseismic period downdip of the fully coupled zone, which imply that the orientation of strain accumulation within the coupled zone can vary with depth. Interseismic GPS motions suggest that forces derived from relative plate motions across a shallow, offshore locked plate interface dominate over decadal timescales. Deeper on the plate interface, below the locked (seismogenic) patch, slip during SSEs dominantly occurs in the updip direction, reflecting a dip‐parallel force acting on the slab, such as slab pull. This implies that in subduction zones with obliquely convergent plate motions, the seismogenic zone of the megathrust is loaded by forces acting in two discrete directions, leading to a depth‐varying orientation of strain accumulation on the plate interface. 
    more » « less