The Lituya Bay landslide-generated mega-tsunami – numerical simulation and sensitivity analysis

Abstract

Abstract. The 1958 Lituya Bay landslide-generated mega-tsunami is simulated using the Landslide-HySEA model, a recently developed finite-volume Savage–Hutter shallow water coupled numerical model. Two factors are crucial if the main objective of the numerical simulation is to reproduce the maximal run-up with an accurate simulation of the inundated area and a precise recreation of the known trimline of the 1958 mega-tsunami of Lituya Bay: first, the accurate reconstruction of the initial slide and then the choice of a suitable coupled landslide–fluid model able to reproduce how the energy released by the landslide is transmitted to the water and then propagated. Given the numerical model, the choice of parameters appears to be a point of major importance, which leads us to perform a sensitivity analysis. Based on public domain topo-bathymetric data, and on information extracted from the work of Miller (1960), an approximation of Gilbert Inlet topo-bathymetry was set up and used for the numerical simulation of the mega-event. Once optimal model parameters were set, comparisons with observational data were performed in order to validate the numerical results. In the present work, we demonstrate that a shallow water type of model is able to accurately reproduce such an extreme event as the Lituya Bay mega-tsunami. The resulting numerical simulation is one of the first successful attempts (if not the first) at numerically reproducing, in detail, the main features of this event in a realistic 3-D basin geometry, where no smoothing or other stabilizing factors in the bathymetric data are applied.

Seismic observations, numerical modeling, and geomorphic analysis of a glacier lake outburst flood in the Himalayas

“…If a future event were to cause the same 16 to 23 m of lake surface lowering as occurred during the 1994 GLOF, then this would translate to 22 to 32 million m 3 of water released or approximately 30% more than the 1994 flood volume. While satellite imagery indicates that the lake boundaries (excluding the retreating glacier terminus) are nearly unchanged since the GLOF, the ongoing glacier retreat has also exposed unstable steep valley walls and lateral moraines above the lake ( 37 , 38 ). A large mass movement into the water could result in a sudden increase in hydrostatic pressure and subsequent overtopping or structural failure of the Lugge Tsho moraine dam.…”

Multisource Bayesian Probabilistic Tsunami Hazard Analysis for the Gulf of Naples (Italy)

“…Even if more recent and updated numerical codes exist (e.g., González-Vida et al, 2019;Macías et al, 2016), the COMCOT model fits ours needs, for the possibility of homogeneously covering all the considered sources, as well as for its relatively easy implementation and its availability online. The COMCOT model is able to simulate tsunami waves with an explicit leap-frog finite difference scheme and multisize nested grids.…”

Ocean Observations Required to Minimize Uncertainty in Global Tsunami Forecasts, Warnings, and Emergency Response

“…In 2019, the United States tsunami warning centers will deploy a new version of the forecast system, which is based on a version of the Method of Splitting Tsunamis (MOST) numerical model (Titov and González, 1997;De la Asunción et al, 2013;Titov et al, 2016) that takes advantage of the massive parallel architecture of modern graphics processing units (GPUs). Other tsunami warning centers in Japan and Europe are in the process of adopting similar GPU technology for tsunami forecasting (González-Vida et al, 2019).…”