Search

Abstract We considered various non‐uniformities such as branch faults, rotation of stress field directions, and changes in tectonic environments to simulate the dynamic rupture process of the 6 February 2023 Mw 7.8 Kahramanmaraş earthquake in SE Türkiye. We utilized near‐fault waveform data, GNSS static displacements, and surface rupture to constrain the dynamic model. The results indicate that the high initial stress accumulated in the Kahramanmaraş‐Çelikhan seismic gap leads to the successful triggering of the East Anatolian Fault (EAF) and the supershear rupture in the northeast segment. Due to the complexity of fault geometry, the rupture speed along the southeastern segment of the EAF varied repeatedly between supershear and subshear, which contributed to the unexpectedly strong ground motion. Furthermore, the triggering of the EAF reminds us to be aware of the risk of seismic gaps on major faults being triggered by secondary faults, which is crucial to prevent significant disasters.

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 leads to localized up- or downdip propagating supershear rupture depending on hypocentral depth. Albeit dynamic earthquakes differ (rupture speed, peak slip-rate, fault slip, bimaterial effects), the effects of hypocentral depth (25–40 km) on tsunami dynamics are negligible. Lateral hypocenter variations lead to small effects such as delayed wave arrival of up to 100 s and differences in tsunami amplitude of up to 0.4 m at the coast. We next analyse inundation on a coastline with complex topo-bathymetry which increases tsunami wave amplitudes up to ≈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 ν = 0.3 which results in close to double the amount of shallow fault slip, ≈1.5 m higher vertical seafloor displacement, and a difference of up to ≈1.5 m in coastal tsunami amplitudes. Lastly, we model a dynamic rupture “tsunami earthquake” with low rupture velocity and low peak slip rates but 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 (≈7.5 m) and inundation area compared to all other scenarios. Our mechanically consistent analysis for a generic megathrust setting can provide building blocks toward using physics-based dynamic rupture modeling in Probabilistic Tsunami Hazard Analysis.

Abstract Elastic dislocation models require large near‐trench slip to explain large tsunamigenesis, which is probably best exemplified in the 2011 M9 Tohoku earthquake. However, it is puzzling that the largest Tohoku tsunami heights occurred about 100 km north of the largest slip zone, where bathymetric surveys indicate no large slip at the trench or submarine landslides. Here we show that coseismic yielding of plentiful sediments in the northern Japan Trench margin can induce large inelastic uplift landward from the trench and diminish slip near the trench. The scarcity of sediments in the south leads to nearly elastic response with large slip at the trench and mostly horizontal seafloor displacement. Thus, the variations of sediments along the Japan Trench and sediment yielding can explain the puzzling variations of tsunamigenesis and near‐trench slip in this earthquake. Inelastic wedge deformation can be an important mechanism of tsunamigenesis in accretionary and other sediment‐filled plate margins.

Fracture mechanics theory and seismological observations suggest that slip-rate is constantly changing during earthquake rupture, including dramatic acceleration from static conditions to high velocity sliding followed by deceleration and arrest. This slip history is partly determined by a complex frictional evolution, including overcoming peak friction, rapid weakening, and re-strengthening (or healing). Recent experimental developments have allowed friction evolution measurements under realistic slip histories reaching high co-seismic slip-rates of meters per second. Theoretical work has focused on describing the observed steady-state weakening at these high-velocities, but the transient behavior has only been fit by direct parameterizations without state variable dependence, needed to simulate arbitrary slip-histories. Commonly used forms of rate-state friction (RSF) are based on low-velocity, step-change experiments and have been shown to not fit the entire frictional evolution using a single set of realistic parameters. Their logarithmic form precludes zero fault slip-rate, assuming it is never truly static, thus does not capture slip initiation phenomena that might contribute to nucleation behavior. Inverting high slip-rate and friction data from different types of experiments, we show that RSF can work by using parameter ranges far from typical low-velocity values. In comparison, we introduce “bristle-state” friction (BSF) models, developed by control-system engineers to predict the transient frictional evolution during arbitrary stressing, especially reversals through static conditions. Although BSF models were also designed for low-velocities, we show that their form provides advantages for fitting frictional evolution measurements under high slip-rate, long-displacement, non-trivial slip histories, especially during the initial strengthening stage.