Simulated effect of soil depth and bedrock topography on near-surface hydrologic response and slope stability

This paper explores the effect of hillslope hydrological behavior on slope stability in the context of transient subsurface saturation development and landslide triggering. We perform a series of virtual experiments to address how subsurface topography affects the location and spatial pattern of slip surface development and pore pressure dynamics. We use a 3D Darcy–Richards equation solver (Hydrus 3-D) combined with a cellular automata slope stability model to simulate the spatial propagation of the destabilized area. Our results showed that the soil–bedrock interface and in particular, bedrock depressions, played a key role in pore pressure dynamics, acting as an impedance for the downslope drainage of perched water. Filling and spilling of depressions in the bedrock surface microtopography induced localized zones of increased pressure head such that the development of porepressure fields—not predictable by surface topography—lead to rapid landslide propagation. Our work suggests that landslide models should consider the subsurface topography in order to include a connectivity component in the mathematical description of hydrological processes operating at the hillslope scale.