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 the
Stochastic slip rupture modeling is a computationally efficient, reduced‐physics approximation that has the capability to create large numbers of unique ruptures based only on a few statistical assumptions. Yet one fundamental question pertaining to this approach is whether the slip distributions calculated in this way are “realistic.” Rather, can stochastic modeling reproduce slip distributions that match what is seen in
- Award ID(s):
- 1835661
- PAR ID:
- 10490259
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract M w9.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 Dynamic wedge failure produces short‐wavelength seafloor uplift efficiently with diminishing shallow slip on the plate interface, generating impulsive tsunami. For ria coasts with prevalent small‐wavelength bathymetric and geomorphologic features, such as the Sanriku coast of Japan, impulsive tsunami can be amplified to produce large runup. We model tsunami propagation and runup of the 1896 Sanriku tsunami by using the seafloor deformation from dynamic rupture models of Ma and Nie (2019) for a
MW 8 earthquake with inelastic wedge deformation. The nonlinear Boussinesq equation is solved by a nested‐grid finite‐difference method with high‐resolution bathymetry data. We show that an inelastic deformation model with extensive wedge failure produces impulsive tsunami similar to those observed offshore the Sanriku coast in the 2011 Tohoku earthquake and generates large runup remarkably consistent with the 1896 Sanriku tsunami. As an alternative to previous models based solely on fault slip, we suggest that the impulsive tsunami and large runup along the Sanriku coast observed in the 2011 Tohoku earthquake can be explained by inelastic wedge deformation north of 38.5°N. -
more » « less
Abstract We explore the potential of the adjoint‐state tsunami inversion method for rapid and accurate near‐field tsunami source characterization using S‐net, an array of ocean bottom pressure gauges. Compared to earthquake‐based methods, this method can obtain more accurate predictions for the initial water elevation of the tsunami source, including potential secondary sources, leading to accurate water height and wave run‐up predictions. Unlike finite‐fault tsunami source inversions, the adjoint method achieves high‐resolution results without requiring densely gridded Green's functions, reducing computation time. However, optimal results require a dense instrument network with sufficient azimuthal coverage. S‐net meets these requirements and reduces data collection time, facilitating the inversion and timely issuance of tsunami warnings. Since the method has not yet been applied to dense, near‐field data, we test it on synthetic waveforms of the 2011
M w 9.0 Tohoku earthquake and tsunami, including triggered secondary sources. The results indicate that with a static source model without noise, using the first 5 min of the waveforms yields a favorable performance with an average accuracy score of 93%, and the largest error of predicted wave amplitudes ranges between −5.6 and 1.9 m. Using the first 20 min, secondary sources were clearly resolved. We also demonstrate the method's applicability using S‐net recordings of the 2016M w 6.9 Fukushima earthquake. The findings suggest that lower‐magnitude events require a longer waveform duration for accurate adjoint inversion. Moreover, the estimated stress drop obtained from inverting our obtained tsunami source, assuming uniform slip, aligns with estimations from recent studies. -
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
-
more » « less
Abstract Physics‐based numerical modeling of earthquake source processes strives to predict quantities of interest for seismic hazard, such as the probability of an earthquake rupture jumping between fault segments. How to assess the predictive power of numerical models remains a topic of ongoing debate. Here, we investigate how sensitive the outcomes of numerical simulations of sequences of earthquakes and aseismic slip are to choices in numerical discretization and treatment of inertial effects, using a simplified 2‐D crustal fault model with two co‐planar segments separated by a creeping barrier. Our simulations demonstrate that simplifying inertial effects and using oversized cells significantly affects the resulting earthquake sequences, including the rate of two‐segment ruptures. We find that fault models with different properties and modeling assumptions can produce similar frequency‐magnitude statistics and static stress drops but have different rates of two‐segment ruptures. For sufficiently long faults, we find that long‐term sequences of events can substantially differ even among simulations that are well resolved by standard considerations. In such simulations, some outcomes, such as static stress drops, are similar among adequately resolved simulations, whereas others, such as the rate of two‐segment ruptures, can be highly sensitive to numerical procedures and modeling assumptions. While it is possible that the response of models with additional ingredients ‐Realistic fault geometry, fluid effects, etc. ‐Would be less sensitive to numerical procedures, our results emphasize the need to examine the potential dependence of simulation outcomes on the modeling procedures and resolution, particularly when assessing their predictive value for seismic hazard assessment.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2022jb025716
