Submarine Landslide Source Modeling using the 3D Slope Stability Analysis Method for the 2018 Palu-Sulawesi Tsunami

Abstract

Abstract. Studies have indicated that submarine landslides played an important role in the 2018 Sulawesi tsunami event, damaging the coast of Palu Bay in addition to the earthquake source. Most of these studies relied on visible landslides to reproduce tsunamis but could still not fully explain the observational data. Recently, several numerical models included hypothesized submarine landslides that were taken into account to obtain a better explanation of the event. In this study, for the first time, submarine landslides were simulated by applying a numerical model based on Hovland’s 3D slope stability analysis for cohesion-frictional soils. To specify landslide volume and location, the model assumed an elliptical slip surface on a vertical slope of 27 m of mesh-divided terrain and evaluated the minimum safety factor in each mesh area based on the surveyed soil property data extracted from the literature. The landslide output was then substituted into a two-layer numerical model based on a shallow-water equation to simulate tsunami propagation. The results were combined with the other tsunami sources, i.e., earthquakes and observed coastal collapses, and validated with various postevent field observational data, including tsunami runup heights and flow depths around the bay, the inundation area around Palu city, waveforms recorded by the Pantoloan tide gauge, and video-inferred waveforms. The model generated several submarine landslides, with lengths of 0.2–2.0 km throughout Palu Bay. The results confirmed the existence of submarine landslide sources in the southern part of the bay and showed agreement with the observed tsunami data, including runups and flow depths. Furthermore, the simulated landslides also reproduced the video-inferred waveforms in 3 out of 6 locations. Although these calculated submarine landslides still cannot fully explain some of the observed tsunami data, they emphasize the possible submarine landslide locations in southern Palu Bay that should be studied and surveyed in the future.

How Small Slip Surfaces Evolve Into Large Submarine Landslides—Insight From 3D Numerical Modeling

“…Such a multi‐directional propagation mechanism can provide additional driving force because of the reduction in strength during slip surface growth, but this physical failure mechanism cannot be considered in a 2D scenario. A few LEM studies have been able to consider 3D submarine slope stability analysis with rotational sliding mechanisms in the past two decades (Somphong et al., 2022; Sultan et al., 2007). Zhang et al.…”