Antonio Scala

and 9 more

Tsunamis are rare, destructive events, whose generation, propagation and coastal impact processes involve several complex physical phenomena. Most tsunami applications, like probabilistic tsunami hazard assessment, make extensive use of large sets of numerical simulations, facing a systematic trade-off between the computational costs and the modelling accuracy. For seismogenic tsunami, the source is often modelled as an instantaneous sea-floor displacement due to the fault static slip distribution, while the propagation in open-sea is computed through a shallow water approximation. Here, through 1D earthquake-tsunami coupled simulations of large M>8 earthquakes in Tohoku-like subduction zone, we tested for which conditions the instantaneous source (IS) and/or the shallow water (SW) approximations can be used to simulate with enough accuracy the whole tsunami evolution. We used as a reference a time-dependent (TD), multi-layer, non-hydrostatic (NH) model whose source features, duration, and size, are based on seismic rupture dynamic simulations with realistic stress drop and rigidity, within a Tohoku-like environment. We showed that slow ruptures, generating slip in shallow part of subduction slabs (e.g. tsunami earthquakes), and very large events, with an along-dip extension comparable with the trench-coast distance (e.g. mega-thrust) require a TD-NH modelling, in particular when the bathymetry close to the coast features sharp depth gradients. Conversely, deeper, higher stress-drop events can be accurately modelled through an IS-SW approximation. We finally showed to what extent inundation depend on bathymetric geometrical features: (i) steeper bathymetries generate larger inundations and (ii) a resonant mechanism emerges with run-up amplifications associated with larger source size on flatter bathymetries.

Manel Prada

and 6 more

Megathrust earthquakes are strongly influenced by the elastic properties of rocks surrounding the fault. However, these properties are often overestimated in numerical simulations, particularly in the shallow megathrust. Here we explore the influence that realistic depth-varying upper-plate elastic properties along the megathrust have on earthquake rupture dynamics and tsunamigenesis using 3D dynamic rupture and tsunami simulations. We compare results from three subduction zone scenarios with homogeneous and heterogeneous elastic media, and bimaterial fault. Elastic properties in the heterogeneous model follow a realistic depth-distribution derived from controlled-source tomography models of subduction zones. We assume the same friction properties for all scenarios. Simulations in the heterogeneous and homogeneous models show that rigidity variation of the country rock determines the depth-varying behavior of slip, slip rate, frequency content, and rupture time. Fault friction may provide additional constraints, but to a lesser extent. The depth-varying behavior of slip, frequency content, and rupture duration quantitatively agree with previous predictions based on worldwide data compilations, explaining the main depth-dependent traits of tsunami earthquakes and large shallow megathrust earthquakes. Large slip, slow rupture and slip rate amplification in bimaterial simulations are largely controlled by the elastic rock properties of the most compliant side of the fault, which in subduction zones is the upper plate. Large shallow slip and trenchward increasing upper-plate compliance of the heterogeneous model lead to the largest co-seismic seafloor deformation and tsunami amplitude. This highlights the importance of considering realistic variations in upper-plate rigidity to properly assess the tsunamigenic potential of megathrust earthquakes.