An advanced process-based distributed model for the investigation of rainfall-induced landslides: The effect of process representation and boundary conditions

Extreme rainfall events are the major driver of shallow landslide occurrences in mountainous and steep terrain regions around the world. Subsurface hydrology has a dominant role on the initiation of rainfall-induced shallow landslides, since changes in the soil water content affect significantly the soil shear strength. Rainfall infiltration produces an increase of soil water potential, which is followed by a rapid drop in apparent cohesion. Especially on steep slopes of shallow soils, this loss of shear strength can lead to failure even in unsaturated conditions before positive water pressures are developed. We present HYDROlisthisis, a process-based model, fully distributed in space with fine time resolution, in order to investigate the interactions between surface and subsurface hydrology and shallow landslides initiation. Fundamental elements of the approach are the dependence of shear strength on the three-dimensional (3-D) field of soil water potential, as well as the temporal evolution of soil water potential during the wetting and drying phases. Specifically, 3-D variably saturated flow conditions, including soil hydraulic hysteresis and preferential flow phenomena, are simulated for the subsurface flow, coupled with a surface runoff routine based on the kinematic wave approximation. The geotechnical component of the model is based on a multidimensional limit equilibrium analysis, which takes into account the basic principles of unsaturated soil mechanics. A series of numerical simulations were carried out with various boundary conditions and using different hydrological and geotechnical components. Boundary conditions in terms of distributed soil depth were generated using both empirical and process-based models. The effect of including preferential flow and soil hydraulic hysteresis was tested together with the replacement of the infinite slope assumption with the multidimensional limit equilibrium analysis. The results show that boundary conditions play a crucial role in the model performance and that the introduced hydrological (preferential flow and soil hydraulic hysteresis) and geotechnical components (multidimensional limit equilibrium analysis) significantly improve predictive capabilities in the presented case study.