Coupled, Physics-Based Modeling Reveals Earthquake Displacements are Critical to the 2018 Palu, Sulawesi Tsunami

• Published in: Pure and applied geophysics 2019-10, Vol.176 (10), p.4069-4109
• Main Authors: Ulrich, T., Vater, S., Madden, E. H., Behrens, J., van Dinther, Y., van Zelst, I., Fielding, E. J., Liang, C., Gabriel, A.-A.
• Format: Article
• Language: English
• Subjects: Bathymeters; Bathymetry; Displacement; Duration; Earth and Environmental Science; Earth Sciences; Earthquakes; Fault lines; Flooding; Geophysics/Geodesy; Hazard assessment; Modelling; Ocean floor; Physics; Propagation; Rupture; Rupturing; Seismic activity; Seismic engineering; Seismic propagation; Seismic response; Seismic wave propagation; Seismic waves; Slip; Sulawesi/Palu-2018 and Anak/Krakatau-2018; Surveys; Tectonics; Time dependence; Tsunami hazard; Tsunamis; Uplift; Wave propagation; Waveforms; Weather hazards
• ISSN: 0033-4553; 1420-9136
• DOI: 10.1007/s00024-019-02290-5

The September 2018, M w  7.5 Sulawesi earthquake occurring on the Palu-Koro strike-slip fault system was followed by an unexpected localized tsunami. We show that direct earthquake-induced uplift and subsidence could have sourced the observed tsunami within Palu Bay. To this end, we use a physics-based, coupled earthquake–tsunami modeling framework tightly constrained by observations. The model combines rupture dynamics, seismic wave propagation, tsunami propagation and inundation. The earthquake scenario, featuring sustained supershear rupture propagation, matches key observed earthquake characteristics, including the moment magnitude, rupture duration, fault plane solution, teleseismic waveforms and inferred horizontal ground displacements. The remote stress regime reflecting regional transtension applied in the model produces a combination of up to 6 m left-lateral slip and up to 2 m normal slip on the straight fault segment dipping 65 ∘ East beneath Palu Bay. The time-dependent, 3D seafloor displacements are translated into bathymetry perturbations with a mean vertical offset of 1.5 m across the submarine fault segment. This sources a tsunami with wave amplitudes and periods that match those measured at the Pantoloan wave gauge and inundation that reproduces observations from field surveys. We conclude that a source related to earthquake displacements is probable and that landsliding may not have been the primary source of the tsunami. These results have important implications for submarine strike-slip fault systems worldwide. Physics-based modeling offers rapid response specifically in tectonic settings that are currently underrepresented in operational tsunami hazard assessment.