This study aims to test the capacity of Flow-3D regarding the simulation of a rockslide-generated impulse wave by evaluating the influences of the extent of the computational domain, the grid resolution, and the corresponding computation times on the accuracy of modelling results. A detailed analysis of the Lituya Bay tsunami event (1958, Alaska, maximum recorded run-up of 524 m a.s.l.) is presented. A focus is put on the tsunami formation and run-up in the impact area. Several simulations with a simplified bay geometry are performed in order to test the concept of a “denser fluid”, compared to the seawater in the bay, for the impacting rockslide material. Further, topographic and bathymetric surfaces of the impact area are set up. The observed maximum run-up can be reproduced using a uniform grid resolution of 5 m, where the wave overtops the hill crest facing the slide source and then flows diagonally down the slope. The model is extended along the entire bay to simulate the wave propagation. The tsunami trimline is well recreated when using (a) a uniform mesh size of 20 m or (b) a non-uniform mesh size of 15 m × 15 m × 10 m with a relative roughness of 2 m for the topographic surface. The trimline mainly results from the primary wave, and in some locations it also results from reflected waves. The denser fluid is a suitable and simple concept to recreate a sliding mass impacting a waterbody, in this case with maximum impact speed of ∼93 m s⁻¹. The tsunami event and the related trimline are well reproduced using the 3D modelling approach with the density evaluation model available in Flow-3D.

In this work, a simple methodology for preliminarily assessing the magnitude of potential landslide-induced impulse waves’ attenuation in mountain lakes is presented. A set of metrics is used to define the geometries of theoretical mountain lakes of different sizes and shapes and to simulate impulse waves in them using the hydrodynamic software Flow-3D. The modeling results provide the ‘wave decay potential’, a ratio between the maximum wave amplitude and the flow depth at the shoreline. Wave decay potential is highly correlated with what is defined as the ‘shape product’, a metric that represents lake geometry. The relation between these two parameters can be used to evaluate wave dissipation in a natural lake given its geometric properties, and thus estimate expected flow depth at the shoreline. This novel approach is tested by applying it to a real-world event, the 2007 landslide-generated wave in Chehalis Lake (Canada), where the results match well with those obtained using the empirical equation provided by ETH Zurich (2019 Edition). This work represents the initial stage in the development of this method, and it encourages additional research and modeling in which the influence of the impacting characteristics on the resulting waves and flow depths is investigated.

On 17 October 2015, a large-scale subaerial landslide occurred in Taan Fiord, Alaska, which released about 50 Mm${}^3$ of rock. This entered the water body and triggered a tsunami with a runup of up to 193 m. This paper aims to simulate the possible formation of a weak layer in this mountainous slope until collapse, and to analyze the possible triggering factors of this landslide event from a geotechnical engineering perspective so that a deeper understanding of this large landslide event can be gained. We analyzed different remote-sensing datasets to characterize the evolution of the coastal landslide process. Based on the acquired remote-sensing data, Digital Elevation Models were derived, on which we employed a 2D limit equilibrium method in this study to calculate the safety factor and compare the location of the associated sliding surface with the most probable actual location at which this landslide occurred. The calculation results reflect the development process of this slope collapse. In this case study, past earthquakes, rainfall before this landslide event, and glacial melting at the toe may have influenced the stability of this slope. The glacial retreat is likely to be the most significant direct triggering factor for this slope failure. This research work illustrates the applicability of multi-temporal remote sensing data of slope morphology to constrain preliminary slope stability analyses, aiming to investigate large-scale landslide processes. This interdisciplinary approach confirms the effectiveness of the combination of aerial data acquisition and traditional slope stability analyses. This case study also demonstrates the significance of a climate change for landslide hazard assessment, and that the interaction of natural hazards in terms of multi-hazards cannot be ignored.