Impact Factor 2.689 | CiteScore 3.3
More on impact ›


Front. Earth Sci. | doi: 10.3389/feart.2021.626844

3D linked subduction, dynamic rupture, tsunami and inundation modeling: dynamic effects of supershear and tsunami earthquakes, hypocenter location and shallow fault slip Provisionally accepted The final, formatted version of the article will be published soon. Notify me

  • 1Ludwig Maximilian University of Munich, Germany
  • 2University of Brasilia, Brazil
  • 3Technical University of Munich, Germany
  • 4University of Leeds, United Kingdom
  • 5Utrecht University, Netherlands

Physics-based dynamic rupture models capture the variability of earthquake slip in space and time and can account for the structural complexity inherent to subduction zones. Here we link tsunami generation, propagation, and coastal inundation with 3D earthquake dynamic rupture (DR) models initialized using a 2D seismo-thermo-mechanical geodynamic (SC) model simulating both subduction dynamics and seismic cycles.
We analyze a total of 15 subduction-initialized 3D dynamic rupture-tsunami scenarios in which the tsunami source arises from the time-dependent co-seismic seafloor displacements with flat bathymetry and inundation on a linearly sloping beach.
We first vary the location of the hypocenter to generate 12 distinct unilateral and bilateral propagating earthquake scenarios. Large-scale fault topography (and the associated rheological complexity) leads to localized up- or downdip propagating supershear rupture depending on hypocentral depth.
Albeit dynamic earthquake differences (rupture speed, peak slip-rate, fault slip, bimaterial effects), the effects of hypocentral depth (25 to 40~km, all consistent with the SC model) on tsunami dynamics are negligible. Lateral hypocenter variations lead to small effects such as delayed wave arrival of up to 80~s and differences in tsunami amplitude of up to 0.5~m at the coast.
We next analyse inundation on a coastline with complex topo-bathymetry which increases tsunami wave amplitudes up to $\approx$1.5~m compared to a linearly sloping beach.
Motivated by structural heterogeneity in subduction zones, we analyse a scenario with increased Poisson's ratio of $\nu=0.3$ which results in close to double the amount of shallow fault slip, $\approx$1.5~m higher vertical seafloor displacement, and a difference of up to $\approx$0.6~m in coastal tsunami amplitudes.
Lastly, we model a ``tsunami earthquake'' with low rupture velocity and peak slip rates and twice as high tsunami potential energy. We triple fracture energy which again doubles the amount of shallow fault slip, but also causes a 2~m higher vertical seafloor uplift and the highest coastal tsunami amplitude ($\approx$8~m) and inundation area compared to all other scenarios.
Our sensitivity analysis for a generic megathrust setting can provide building blocks towards dynamic rupture modeling complementing Probabilistic Tsunami Hazard Analysis.

Keywords: Earthquake dynamics, High Performance Computing (HPC), Probabilistic tsunami hazard assessment (PTHA), Megathrust earthquakes, Supershear propagation

Received: 06 Nov 2020; Accepted: 26 Feb 2021.

Copyright: © 2021 Wirp, Gabriel, Madden, Schmeller, van Zelst, Krenz, van Dinther and Rannabauer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

* Correspondence:
Ms. Sara A. Wirp, Ludwig Maximilian University of Munich, Munich, 80539, Bavaria, Germany,
Dr. Alice-Agnes Gabriel, Ludwig Maximilian University of Munich, Munich, 80539, Bavaria, Germany,