Published March 16, 2026 | Version Published
Journal Article Open

Large megathrust earthquakes in cold mantle wedge corners under lawsonite blueschist facies

  • 1. ROR icon University of Southern California
  • 2. ROR icon California Institute of Technology
  • 3. Institute of Geology, China Earthquake Administration, Beijing, China

Abstract

Giant megathrust earthquakes typically occur above the upper-plate Moho or the 350°C isotherm in warm subduction zones. However, large earthquakes also occur within the mantle wedge corner at 35-55 km depth at cold subduction zones, such as the Kermadec, Japan Trench, and Chilean margins. Here we investigate the frictional behavior of lawsonite blueschist metagreywacke, potentially found in cold subduction channels, to better understand the control of metamorphic grade on such first-order variations in seismogenic depth. We perform velocity-step experiments on metagreywacke gouge from room temperature to 500°C and effective normal stresses from 50 MPa to 320 MPa, capturing the unstable friction regime and the brittle-to-flow transition. The constitutive behavior of metagreywacke indicates a potential seismogenic behavior at high temperatures below the Moho. Large megathrust earthquakes in the mantle wedge corner may develop in the lawsonite blueschist metasediment channel of cold subduction slabs.

Copyright and License

© 2026, The Author(s). Open Access .This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. 

Acknowledgement

We are grateful for the technical support of Wenming Yao and Changrong He at the Institute of Geology, China Earthquake Administration. This study is funded by the National Natural Science Foundation of China under grant 42274234, the National Key Research and Development Program of China (2022YFC3003704), and the National Science Foundation under award number EAR-1848192. H.Z. acknowledges support from the National Science Foundation grant EAR-2041892 awarded to John E. Vidale.

Data Availability

The laboratory data used in this study are available in a zenodo.org repository.

Code Availability

A demonstration code for simulating velocity-step experiments using the constitutive law proposed in this study is available at https://doi.org/10.5281/zenodo.1852974983. The “Unicycle” software package, along with the input files for the seismic cycle simulation, can be accessed at: https://bitbucket.org/sbarbot/unicycle.

Supplemental Material

Supplementary Information (download PDF )

Files

s41467-026-70315-4_reference.pdf

Files (43.7 MB)

Name Size
md5:bfa743abeef03a428f3cbff50df84799
7.3 MB Preview Download
md5:505fa22028114b4b654f890cbb5adde2
36.4 MB Preview Download

Additional details

Related works

Describes
Journal Article: https://rdcu.be/e9edA (URL)

Funding

National Natural Science Foundation of China
42274234
National Science Foundation
EAR-1848192

Caltech Custom Metadata