Four-Year Monitoring Study of Shallow Landslide Hazards Based on Hydrological Measurements in a Weathered Granite Soil Slope in South Korea

To build a comprehensive understanding of long-term hydro-mechanical processes that lead to shallow landslide hazards, this study explicitly monitored the volumetric water content (VWC) and rainfall amount for a weathered granite soil slope over a four year period. From the 12 operational landslide monitoring stations installed across South Korea, the Songnisan station was selected as the study site. VWC sensors were placed in the subsurface with a grid-like arrangement at depths of 0.5 and 1.0 m. Shallow landslide hazards were evaluated by applying an infinite slope stability model that adopted a previously proposed unified effective stress concept. By analyzing the variations in the monitored VWC values, the derived matric suctions and suction stresses, and the calculated factor of safety values, we were able to obtain numerous valuable insights. In particular, the seasonal effects of drainage and evapotranspiration on the slope moisture conditions and slope stability were addressed. Preliminary test results indicated that continuous rainfall successfully represented the derived matric suction conditions at a depth of 1.0 m in the lower slope, although this was not the case for the upper and middle slopes. The significance of a future study on cumulative field monitoring data from various sites in different geological conditions is highlighted.

Shallow Landslides Physically Based Susceptibility Assessment Improvement Using InSAR. Case Study: Carpathian and Subcarpathian Prahova Valley, Romania

A multi-temporal satellite radar interferometry technique is used for deriving the actual surface displacement patterns in a slope environment in Romania, in order to validate and improve a landslide susceptibility map. The probability the occurrence of future events is established using a deterministic approach based on a classical one-dimension infinite slope stability model. The most important geotechnical parameters for slope failure in the proposed study area are cohesion, unit weight and friction angle, and the triggering factor is a rapid rise in groundwater table under wetting conditions. Employing a susceptibility analysis using the physically based model under completely saturated conditions proved to be the most suitable scenario for identifying unstable areas. The kinematic characteristics are assessed by the Small BAseline Subsets (SBAS) interferometry technique applied to C-band synthetic aperture radar (SAR) Sentinel-1 imagery. The analysis was carried out mainly for inhabited areas which present a better backscatter return. The validation revealed that more than 22% of the active landslides identified by InSAR were predicted as unstable areas by the infinite slope model. We propose a refinement of the susceptibility map using the InSAR results for unravelling the danger of the worst-case scenario.

Short and long term probabilistic slope stability analyses of a large area of unsaturated pyroclastic soils

<p>The analysis of slope stability over large areas is a demanding task for several reasons, such as the need for extensive datasets, the uncertainty of collected data, the difficulty of accounting for site-specific factors, and the considerable computation time required due to the size of investigated areas, which can pose major barriers, particularly in civil protection contexts where rapid analysis and forecasts are essential. However, as the identification of zones of higher failure probability is very useful for stakeholders and decision-makers, the scientific community has attempted to improve capabilities to provide physically based assessments. This study combined a transient seepage analysis of an unsaturated-saturated condition with an infinite slope stability model and probabilistic analysis through the use of a high-computing capacity parallelized platform. Both short- and long-term analyses were performed for a study area, and roles of evapotranspiration, vegetation interception, and the root increment of soil strength were considered. A model was first calibrated based on hourly rainfall data recorded over a 4-day event (December 14–17, 1999) causing destructive landslides to compare the results of model simulations to actual landslide events. Then, the calibrated model was applied for a long-term simulation where daily rainfall data recorded over a 4-year period (January 1, 2005–December 31, 2008) were considered to study the behavior of the area in response to a long period of rainfall. The calibration shows that the model can correctly identify higher failure probability within the time range of the observed landslides as well as the extents and locations of zones computed as the most prone ones. The long-term analysis allowed for the identification of a number of days when the slope factor of safety was lower than 1.2 over a significant number of cells. In all of these cases, zones approaching slope instability were concentrated in specific sectors and catchments of the study area. In addition, some subbasins were found to be the most recurrently prone to possible slope instability. Interestingly, the application of the adopted methodology provided clear indications of both weekly and seasonal fluctuations of overall slope stability conditions.</p>

Rainfall-induced shallow landslides and soil wetness: comparison of physically-based and probabilistic predictions

Landslides are an impacting natural hazard in alpine regions, calling for effective forecasting and warning systems. Here we compare two methods (physically-based and probabilistic) for the prediction of shallow rainfall-induced landslides in an application to Switzerland, with a specific focus on the value of antecedent soil wetness. First, we show that landslide susceptibility predicted by the factor of safety in the infinite slope model is strongly dependent on soil data inputs, limiting the hydrologically active range where landslides can occur to only ~20 % of the area with typical soil parameters and soil depth models. Second, the physically-based approach with a coarse resolution model setup (TerrSysMP) 12.5 km × 12.5 km downscaled to 25 m × 25 m with the TopographicWetness Index to provide water table simulations for the infinite slope stability model did not succeed in predicting local scale landsliding satisfactorily, despite spatial downscaling. We argue that this is due to inadequacies of the infinite slope model, soil parameter uncertainty, and the coarse resolution of the hydrological model. Third, soil saturation estimates provided by a higher resolution 500 m × 500 m conceptual hydrological model (PREVAH) provided added value to rainfall threshold curves for landslide prediction in the probabilistic approach, with potential to reduce false alarms and misses. We conclude that although combined physically-based hydrological-geotechnical modelling is the desired goal, we still need to overcome problems of model resolution, parameter constraints, and landslide validation for successful prediction at regional scales.

Multiseasonal probabilistic slope stability analysis of a large area of unsaturated pyroclastic soils

The analysis of slope stability over large areas is a demanding task for several reasons, such as the need for extensive datasets, the uncertainty of collected data, the difficulty of accounting for site-specific factors, and the considerable computation time required due to the size of investigated areas, which can pose major barriers, particularly in civil protection contexts where rapid analysis and forecasts are essential. However, as the identification of zones of higher failure probability is very useful for stakeholders and decision-makers, the scientific community has attempted to improve capabilities to provide physically based assessments. This study combined a transient seepage analysis of an unsaturated-saturated condition with an infinite slope stability model and probabilistic analysis through the use of a high-computing capacity parallelized platform. Both short- and long-term analyses were performed for a study area, and roles of evapotranspiration, vegetation interception, and the root increment of soil strength were considered. A model was first calibrated based on hourly rainfall data recorded over a 4-day event (December 14–17, 1999) causing destructive landslides to compare the results of model simulations to actual landslide events. Then, the calibrated model was applied for a long-term simulation where daily rainfall data recorded over a 4-year period (January 1, 2005–December 31, 2008) were considered to study the behavior of the area in response to a long period of rainfall. The calibration shows that the model can correctly identify higher failure probability within the time range of the observed landslides as well as the extents and locations of zones computed as the most prone ones. The long-term analysis allowed for the identification of a number of days (9) when the slope factor of safety was lower than 1.2 over a significant number of cells. In all of these cases, zones approaching slope instability were concentrated in specific sectors and catchments of the study area. In addition, some subbasins were found to be the most recurrently prone to possible slope instability. Interestingly, the application of the adopted methodology provided clear indications of both weekly and seasonal fluctuations of overall slope stability conditions. Limitations of the present study and future developments are finally discussed.

Probabilistic landslide susceptibility analysis in tropical mountainous terrain using the physically based r.slope.stability model

Landslides triggered by rainfall are very common phenomena in complex tropical environments such as the Colombian Andes, one of the regions of South America most affected by landslides every year. Currently in Colombia, physically based methods for landslide hazard mapping are mandatory for land use planning in urban areas. In this work, we perform probabilistic analyses with r.slope.stability, a spatially distributed, physically based model for landslide susceptibility analysis, available as an open-source tool coupled to GRASS GIS. This model considers alternatively the infinite slope stability model or the 2.5-D geometry of shallow planar and deep-seated landslides with ellipsoidal or truncated failure surfaces. We test the model in the La Arenosa catchment, northern Colombian Andes. The results are compared to those yielded with the corresponding deterministic analyses and with other physically based models applied in the same catchment. Finally, the model results are evaluated against a landslide inventory using a confusion matrix and receiver operating characteristic (ROC) analysis. The model performs reasonably well, the infinite slope stability model showing a better performance. The outcomes are, however, rather conservative, pointing to possible challenges with regard to the geotechnical and geo-hydraulic parameterization. The results also highlight the importance to perform probabilistic instead of – or in addition to – deterministic slope stability analyses.

Comparison of the performance of spatial landslide prediction with TRIGRS1D and SCOOPS3D models and parameter optimization: application to the Oltrepò Pavese

<p>Physically based models based on the combination of hydrological and slope stability models are important tools in spatial and temporal prediction of landslides, since they can be used for hazard mapping as an aid for land planning. In many applications, hydrological models are combined with very simple infinite slope stability analysis, given that multi-dimensional analysis is more computationally demanding. Only a few studies have attempted to apply such algorithms to the catchment scale. Thus, there is a need for more studies on this issue, also to understand the real advantages of applying multi-dimensional slope stability analysis in comparison with the one-dimensional. </p><p>This study aims to compare the performance of two different forecasting models, namely the infinite slope and the three-dimensional stability analysis by SCOOPS3D (Software to analyze three-dimensional slope stability throughout a digital landscape), a very efficient model proposed by USGS to be applied to the catchment scale, which has seldom been applied so far in the literature. In particular, TRIGRS (Transient Rainfall Infiltration and Grid-Based Regional Slope-stability Model) is used for hydrological analysis.  Then the resulting pressure head field is used first as input to the infinite slope stability model embedded into TRIGRS program itself and then as input to SCOOPS3D. To calibrate the terrain stability-related parameters of either piece of software, a multi-objective optimization is proposed in this work to maximize the model predictability performance, in an attempt to optimize ROC performance statistics, i.e. to maximize the true positive rate while simultaneously minimizing the false positive rate.</p><p>The approach was applied to a real case study, a catchment in the Oltrepò Pavese (northern Italy), in which the areas of triggered landslides were accurately monitored during an extreme rainfall on 27-28 April, 2009, featuring 160 mm in 48 h. Compared to other works in the scientific literature, in which only a generic point of location of landslides was known, the present work benefits from the availability of a detailed landslide inventory containing observed landslide shapes.</p><p>The results point out the significantly better performance of  SCOOPS3D, in comparison with the infinite slope stability. Though SCOOPS3D seems to overestimate landslide prone areas, the 3D method is more realistic than the 1D method as far as the slip surface definition is concerned. Therefore, the proposed methodology, lying in the use of SCOOPS 3D with optimized parameters, can be a helpful tool for providing multiple landslide hazard maps for planning.</p>

A new approach to mapping landslide hazards: a probabilistic integration of empirical and physically based models in the North Cascades of Washington, USA

We developed a new approach for mapping landslide hazards by combining probabilities of landslide impacts derived from a data-driven statistical approach and a physically based model of shallow landsliding. Our statistical approach integrates the influence of seven site attributes (SAs) on observed landslides using a frequency ratio (FR) method. Influential attributes and resulting susceptibility maps depend on the observations of landslides considered: all types of landslides, debris avalanches only, or source areas of debris avalanches. These observational datasets reflect the detection of different landslide processes or components, which relate to different landslide-inducing factors. For each landslide dataset, a stability index (SI) is calculated as a multiplicative result of the frequency ratios for all attributes and is mapped across our study domain in the North Cascades National Park Complex (NOCA), Washington, USA. A continuous function is developed to relate local SI values to landslide probability based on a ratio of landslide and non-landslide grid cells. The empirical model probability derived from the debris avalanche source area dataset is combined probabilistically with a previously developed physically based probabilistic model. A two-dimensional binning method employs empirical and physically based probabilities as indices and calculates a joint probability of landsliding at the intersections of probability bins. A ratio of the joint probability and the physically based model bin probability is used as a weight to adjust the original physically based probability at each grid cell given empirical evidence. The resulting integrated probability of landslide initiation hazard includes mechanisms not captured by the infinite-slope stability model alone. Improvements in distinguishing potentially unstable areas with the proposed integrated model are statistically quantified. We provide multiple landslide hazard maps that land managers can use for planning and decision-making, as well as for educating the public about hazards from landslides in this remote high-relief terrain.