Assessing the Optimal Tsunami Inundation Modeling Strategy for Large Earthquakes in Subduction Zones

Tsunamis are rare events involving several complex physical phenomena. Due to this complexity and the relative scarcity of observations, tsunami research makes extensive use of numerical simulations. For seismogenic tsunamis, the source is often modeled as an instantaneous sea-floor displacement (IS), while the tsunami propagation and inundation is computed through a shallow water approximation (SW). Here, we investigate what is the best tsunami inundation modeling strategy for different realistic earthquake source size and duration. We use 1D earthquake-tsunami coupled simulations of large M > 8 earthquakes in Tohoku-like subduction zone to test for which conditions the IS and/or the SW approximations can simulate with enough accuracy the tsunami evolution. We use as a reference a time-dependent (TD), multi-layer, non-hydrostatic (NH) 1D model. Source duration, and size, are based on 1D dynamic rupture simulations with realistic stress drop and rigidity. We show that slow ruptures, generating slip in the shallow part of subduction zones (e.g., tsunami earthquakes), and very large events, with an along-dip extent comparable with the trench-coast distance (as occurs for megathrust events) require a TD-NH modeling, especially for regions with steep coastal bathymetry. Conversely, deeper, higher stress-drop events can be modeled through an IS-SW approximation. We finally show that: (a) steeper bathymetries generate larger runups and, (b) a resonant mechanism emerges with runup amplifications associated with larger source size on flatter bathymetries. These results, obtained with 1D modeling, can serve as a guide for the appropriate 2/3D simulation approach for applications ranging from fundamental tsunami science to computational-intensive hazard assessments.