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: Controls on Fore‐Arc Deformation and Stress Switching After the Great 2011 Tohoku‐Oki Earthquake From Discrete Numerical Simulations
Abstract The large magnitude of the 2011Mw9.0 Tohoku‐Oki earthquake, which occurred off the east coast of Japan, was not expected or predicted by any previous studies. One surprising observation was the sudden change in stress state; local earthquakes confirmed a compressional stress state before the main shock, whereas an extensional stress state was evident after the main shock. Using discrete element method modeling, this project attempts to reproduce the stress change after the main shock, and explores the conditions that can cause stress switching both onshore and offshore Tohoku. Our simulations demonstrate that rapid fault weakening can produce stress change and predominant normal‐fault earthquake mechanisms in the upper plate of Tohoku‐Oki. Several specific conditions seem to favor such stress switching; the megathrust fault must have been frictionally strong before the main shock, and comparable values of internal (μinternal) and basal friction (μbasal) are necessary to cause the formation of widespread normal faults within the wedge. Furthermore, dynamic extension during elastic unloading appears to play an important role in accommodating stress changes and wedge deformation in the Tohoku area; these cannot be explained solely based on Critical Coulomb Wedge models, but instead require dynamic unloading processes.  more » « less
Award ID(s):
1723249
PAR ID:
10455117
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
124
Issue:
8
ISSN:
2169-9313
Page Range / eLocation ID:
p. 9265-9279
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We determine the trace element stratigraphy of Site C0019, drilled during the Japan Fast Trench Drilling Project (JFAST) International Ocean Discovery Program (IODP) Expedition 343, to illuminate the structure of the plate boundary following the Tohoku-Oki earthquake of 2011. The stratigraphic units at the JFAST site are compared to undeformed Western Pacific sediments from two reference sites (Ocean Drilling Program (ODP) Site 1149 and Deep Sea Drilling Project (DSDP) Site 436). The trace element fingerprints in these reference sedimentary units can be correlated to individual JFAST samples. At the JFAST site, we find that the accretionary wedge and downgoing plate sediments in the core are composed primarily of Holocene to Eocene sediments. There are several age reversals and gaps within the sequence, consistent with multiple faults in the bottom 15 m of the JFAST core. Our results point to several candidate faults that could have slipped during the 2011 Tohoku-Oki earthquake, in addition to the pelagic clay layer that has been proposed as the main décollement fault. 
    more » « less
  2. Abstract Fault‐damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale‐dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major listric fault and its damage zone. We investigate 10 distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited off‐fault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (a) the fracture energy scales with fracture and fault size, (b) favorable relative pre‐stress of fractures within the ambient stress field, and (c) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that fractures with favorable relative pre‐stress, embedded within a fault damage zone, may lead to cascading earthquake rupture that shadows main fault slip. We find that such off‐fault events may reach moment magnitudes up toMw ≈ 5.5, comparable to magnitudes that can be otherwise hosted by the main fault. Our findings offer insights into physical processes governing cascading earthquake dynamic rupture within multiscale fracture networks. 
    more » « less
  3. Abstract Large earthquakes rupture faults over hundreds of kilometers within minutes. Finite‐fault models image these processes and provide observational constraints for understanding earthquake physics. However, finite‐fault inversions are subject to non‐uniqueness and uncertainties. The diverse range of published models for the well‐recorded 2011 9.0 Tohoku‐Oki earthquake illustrates this challenge, and its rupture process remains under debate. Here, we comprehensively compare 32 published finite‐fault models of the Tohoku‐Oki earthquake. We aim to identify the most coherent slip features of the Tohoku‐Oki earthquake from these slip models and develop a new method for quantitatively analyzing their variations. We find that the models correlate poorly at 1‐km subfault size, irrespective of the data type. In contrast, model agreement improves significantly with increasing subfault sizes, consistently showing that the largest slip occurs up‐dip of the hypocenter near the trench. We use the set of models to test the sensitivity of available teleseismic, regional seismic, and geodetic observations. For the large Tohoku‐Oki earthquake, we find that the analyzed finite‐fault models are less sensitive to slip features smaller than 64 km. When we use the models to compute synthetic seafloor deformation, we observe strong variations in the synthetics, suggesting their sensitivity to small‐scale slip features. Our newly developed approach offers a quantitative framework to identify common features in distinct finite‐fault slip models and to analyze their robustness using regional and global geophysical observations for megathrust earthquakes. Our results indicate that dense offshore instrumentation is critical for resolving the rupture complexities of megathrust earthquakes. 
    more » « less
  4. The 11 March 2011 M 9.0 Tohoku-oki earthquake was one of the largest earthquakes ever recorded and was accompanied by a devastating tsunami. Slip during the earthquake was exceptionally large at shallow depth on the plate boundary fault, which was one of the primary factors that contributed to the extreme tsunami amplitudes that inundated the coast of Japan. International Ocean Discovery Program Expedition 405 aims to investigate the conditions and processes that facilitated the extremely shallow slip on the subduction interface in the 2011 Tohoku-oki earthquake. Proposed work includes coring and logging operations at two sites in a transect across the trench. The first site, located within the overriding plate, will access the fault zone in the region of large shallow slip, targeting the plate boundary décollement, overlying frontal prism, and subducted units cut by the décollement. The second site, located on the Pacific plate, will access the undisturbed sedimentary and volcanic inputs to the subduction zone. A borehole observatory will be installed into the décollement and surrounding rocks to provide measurements of the temperature in and around the fault over the following several years. Sampling, geophysical logs, and the observatory temperature time series will document the compositional, structural, mechanical, and frictional properties of the rocks in the décollement and adjacent country rock, as well as the hydrogeologic structure and pore fluid pressure of the fault zone and frontal prism—key properties that influence the effective stress to facilitate earthquake slip and potential for large slip. Results from Expedition 405 will address fundamental questions about earthquake slip on subduction zones that may directly inform earthquake and tsunami hazard assessments around the world. 
    more » « less
  5. Abstract We study the effects of pore fluid pressure (Pf) on the pre‐earthquake, near‐fault stress state, and 3‐D earthquake rupture dynamics through six scenarios utilizing a structural model based on the 2004Mw9.1 Sumatra‐Andaman earthquake. As pre‐earthquakePfmagnitude increases, effective normal stress and fault shear strength decrease. As a result, magnitude, slip, peak slip rate, stress drop, and rupture velocity of the scenario earthquakes decrease. Comparison of results with observations of the 2004 earthquake support that pre‐earthquakePfaverages near 97% of lithostatic pressure, leading to pre‐earthquake average shear and effective normal tractions of 4–5 and 22 MPa. The megathrust in these scenarios is weak, in terms of low mean shear traction at static failure and low dynamic friction coefficient during rupture. Apparent co‐seismic principal stress rotations and absolute post‐seismic stresses in these scenarios are consistent with the variety of observed aftershock focal mechanisms. In all scenarios, the mean apparent stress rotations are larger above than below the megathrust. Scenarios with largerPfmagnitudes exhibit lower mean apparent principal stress rotations. We further evaluate pre‐earthquakePfdepth distribution. IfPffollows a sublithostatic gradient, pre‐earthquake effective normal stress increases with depth. IfPffollows the lithostatic gradient exactly, then this normal stress is constant, shifting peak slip and peak slip rate updip. This renders constraints on near‐trench strength and constitutive behavior crucial for mitigating hazard. These scenarios provide opportunity for future calibration with site‐specific measurements to constrain dynamically plausible megathrust strength andPfgradients. 
    more » « less