An official website of the United States government
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
This content will become publicly available on April 27, 2027
Dynamic restrengthening and fault heterogeneity explain megathrust earthquake complexity
Abstract Megathrusts host Earth’s largest earthquakes. Understanding the physical conditions controlling their rupture dynamics is critical for assessing seismic and tsunami hazards. These earthquakes often display complex rupture dynamics, exemplified by the 2011 Tohoku-Oki earthquake, which exhibited multiple rupture episodes, depth-dependent seismic radiation, and substantial tsunamigenic slip near the trench. However, how such complexity arises from preexisting physical conditions remains uncertain. Here, we demonstrate that the observed rupture complexity of the Tohoku-Oki earthquake can spontaneously and self-consistently emerge, driven by rapid coseismic frictional restrengthening and data-informed fault heterogeneity. We use an ensemble of 3D dynamic rupture simulations to identify that mixed downdip pulse-like and updip crack-like rupture are driven by dynamic stress redistribution with episodic rupture reactivation. By featuring low fault strength compared to its dynamic stress drop, a preferred model can consistently reproduce the observed complex depth-dependent propagation speeds, multiple rupture fronts as imaged by back-projection, and large tsunamigenic slip at the trench. Our findings demonstrate that preexisting fault heterogeneity conjointly with dynamic frictional weakening and restrengthening drives seemingly unexpected megathrust rupture complexity, highlighting the need to include dynamic effects into physics-based seismic and tsunami hazard assessments of future earthquakes.
more »
« less
- PAR ID:
- 10690639
- Publisher / Repository:
- Nature Communications
- Date Published:
- Journal Name:
- Nature Communications
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Extreme slip at shallow depths on subduction zone faults is a primary contributor to tsunami generation by earthquakes. Improving earthquake and tsunami risk assessment requires understanding the material and structural conditions that favor earthquake propagation to the trench. We use new biomarker thermal maturity indicators to identify seismic faults in drill core recovered from the Japan Trench subduction zone, which hosted 50 m of shallow slip during theMw9.1 2011 Tohoku-Oki earthquake. Our results show that multiple faults have hosted earthquakes with displacement ≥ 10 m, and each could have hosted many great earthquakes, illustrating an extensive history of great earthquake seismicity that caused large shallow slip. We find that lithologic contrasts in frictional properties do not necessarily determine the likelihood of large shallow slip or seismic hazard.more » « less
-
Abstract The Mediterranean Hellenic Arc subduction zone (HASZ) has generated several 8 earthquakes and tsunamis. Seismic‐probabilistic tsunami hazard assessment typically utilizes uniform or stochastic earthquake models, which may not represent dynamic rupture and tsunami generation complexity. We present an ensemble of ten 3D dynamic rupture earthquake scenarios for the HASZ, utilizing a realistic slab geometry. Our simplest models use uniform along‐arc pre‐stresses or a single circular initial stress asperity. We then introduce progressively more complex models varying initial shear stress along‐arc, multiple asperities based on scale‐dependent critical slip weakening distance, and a most complex model blending all aforementioned heterogeneities. Thereby, regional initial conditions are constrained without relying on detailed geodetic locking models. Varying epicentral locations in the simplest, homogeneous model leads to different rupture speeds and moment magnitudes. We observe dynamic fault slip penetrating the shallow slip‐strengthening region and affecting seafloor uplift. Off‐fault plastic deformation can double vertical seafloor uplift. A single‐asperity model generates a 8 scenario resembling the 1303 Crete earthquake. Using along‐strike varying initial stresses results in 8.0–8.5 dynamic rupture scenarios with diverse slip rates and uplift patterns. The model with the most heterogeneous initial conditions yields a 7.5 scenario. Dynamic rupture complexity in prestress and fracture energy tends to lower earthquake magnitude but enhances tsunamigenic displacements. Our results offer insights into the dynamics of potential large Hellenic Arc megathrust earthquakes and may inform future physics‐based joint seismic and tsunami hazard assessments.more » « less
-
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
-
Tsunamigenic earthquakes pose considerable risks, both economically and socially, yet earthquake and tsunami hazard assessments are typically conducted separately. Earthquakes associated with unexpected tsunamis, such as the 2018 Mw 7.5 strike-slip Sulawesi earthquake, emphasize the need to study the tsunami potential of active submarine faults in different tectonic settings. Here, we investigate physics-based scenarios combining simulations of 3D earthquake dynamic rupture and seismic wave propagation with tsunami generation and propagation. We present time-dependent modeling of one-way linked and 3D fully coupled earthquakes and tsunamis for the ∼ 100 km long Húsavík–Flatey Fault Zone (HFFZ) in North Iceland. Our analysis shows that the HFFZ has the potential to generate sizable tsunamis. The six dynamic rupture models sourcing our tsunami scenarios vary regarding hypocenter location, spatiotemporal evolution, fault slip, and fault structure complexity but coincide with historical earthquake magnitudes. Earthquake dynamic rupture scenarios on a less segmented fault system, particularly with a hypocenter location in the eastern part of the fault system, have a larger potential for local tsunami generation. Here, dynamically evolving large shallow fault slip (∼ 8 m), near-surface rake rotation (± 20∘), and significant coseismic vertical displacements of the local bathymetry (± 1 m) facilitate strike-slip faulting tsunami generation. We model tsunami crest to trough differences (total wave heights) of up to ∼ 0.9 m near the town Ólafsfjörður. In contrast, none of our scenarios endanger the town of Akureyri, which is shielded by multiple reflections within the narrow Eyjafjörður bay and by Hrísey island. We compare the modeled one-way linked tsunami waveforms with simulation results using a 3D fully coupled approach. We find good agreement in the tsunami arrival times and location of maximum tsunami heights. While seismic waves result in transient motions of the sea surface and affect the ocean response, they do not appear to contribute to tsunami generation. However, complex source effects arise in the fully coupled simulations, such as tsunami dispersion effects and the complex superposition of seismic and acoustic waves within the shallow continental shelf of North Iceland. We find that the vertical velocity amplitudes of near-source acoustic waves are unexpectedly high – larger than those corresponding to the actual tsunami – which may serve as a rapid indicator of surface dynamic rupture. Our results have important implications for understanding the tsunamigenic potential of strike-slip fault systems worldwide and the coseismic acoustic wave excitation during tsunami generation and may help to inform future tsunami early warning systems.more » « less
