Automatic subdivision of landscapes into terrain units remains a challenge. Slope units are terrain units bounded by drainage and divide lines, but their use in hydrological and geomorphological studies is limited because of the lack of reliable software for their automatic delineation. We present the r.slopeunits software for the automatic delineation of slope units, given a digital elevation model and a few input parameters. We further propose an approach for the selection of optimal parameters controlling the terrain subdivision for landslide susceptibility modeling. We tested the software and the optimization approach in central Italy, where terrain, landslide, and geo-environmental information was available. The software was capable of capturing the variability of the landscape and partitioning the study area into slope units suited for landslide susceptibility modeling and zonation. We expect r.slopeunits to be used in different physiographical settings for the production of reliable and reproducible landslide susceptibility zonations.

The relation between climate change and its potential effects on the stability of slopes remains an open issue. For rainfall induced landslides, the point consists in determining the effects of the projected changes in the duration and amounts of rainfall that can initiate slope failures. We investigated the relationship between fine-scale climate projections obtained by downscaling and the expected modifications in landslide occurrence in Central Italy. We used rainfall measurements taken by 56 rain gauges in the 9-year period 2003-2011, and the RainFARM technique to generate downscaled synthetic rainfall fields from regional climate model projections for the 14-year calibration period 2002-2015, and for the 40-year projection period 2010-2049. Using a specific algorithm, we extracted a number of rainfall events, i.e. rainfall periods separated by dry periods of no or negligible amount of rain, from the measured and the synthetic rainfall series. Then, we used the selected rainfall events to forcethe Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Model TRIGRS v. 2.1. We analyzed the results in terms of variations (or lack of variations) in the rainfall thresholds for the possible initiation of landslides, in the probability distribution of landslide size (area), and in landslide hazard. Results showed that the downscaled rainfall fields obtained by RainFARM can be used to single out rainfall events, and to force the slope stability model. Results further showed that while the rainfall thresholds for landslide occurrence are expected to change in future scenarios, the probability distribution of landslide areas are not. We infer that landslide hazard in the study area is expected to change in response to the projected variations in the rainfall conditions. We expect our results to contribute to regional investigations of the expected impact of projected climate variations on slope stability conditions and on landslide hazards.

Distributed models to forecast the spatial and temporal occurrence of rainfall-induced shallow landslides are based on deterministic laws. These models extend spatially the static stability models adopted in geotechnical engineering, and adopt an infinite-slope geometry to balance the resisting and the driving forces acting on the sliding mass. An infiltration model is used to determine how rainfall changes pore-water conditions, modulating the local stability/instability conditions. A problem with the operation of the existing models lays in the difficulty in obtaining accurate values for the several variables that describe the material properties of the slopes. The problem is particularly severe when the models are applied over large areas, for which sufficient information on the geotechnical and hydrological conditions of the slopes is not generally available. To help solve the problem, we propose a probabilistic Monte Carlo approach to the distributed modeling of rainfall-induced shallow landslides. For this purpose, we have modified the transient rainfall infiltration and grid-based regional slope-stability analysis (TRIGRS) code. The new code (TRIGRS-P) adopts a probabilistic approach to compute, on a cell-by-cell basis, transient pore-pressure changes and related changes in the factor of safety due to rainfall infiltration. Infiltration is modeled using analytical solutions of partial differential equations describing one-dimensional vertical flow in isotropic, homogeneous materials. Both saturated and unsaturated soil conditions can be considered. TRIGRS-P copes with the natural variability inherent to the mechanical and hydrological properties of the slope materials by allowing values of the TRIGRS model input parameters to be sampled randomly from a given probability distribution. The range of variation and the mean value of the parameters can be determined by the usual methods used for preparing the TRIGRS input parameters. The outputs of several model runs obtained varying the input parameters are analyzed statistically, and compared to the original (deterministic) model output. The comparison suggests an improvement of the predictive power of the model of about 10% and 16% in two small test areas, that is, the Frontignano (Italy) and the Mukilteo (USA) areas. We discuss the computational requirements of TRIGRS-P to determine the potential use of the numerical model to forecast the spatial and temporal occurrence of rainfall-induced shallow landslides in very large areas, extending for several hundreds or thousands of square kilometers. Parallel execution of the code using a simple process distribution and the message passing interface (MPI) on multi-processor machines was successful, opening the possibly of testing the use of TRIGRS-P for the operational forecasting of rainfall-induced shallow landslides over large regions.