Next Article in Journal
Application of the Coastal and Marine Ecological Classification Standard (CMECS) to Create Benthic Geologic Habitat Maps for Portions of Acadia National Park, Maine, USA
Previous Article in Journal
The Relationships between Greenstone Belts and the Kryvyi Rih–Kremenchuk Basin in the Middle Dnieper Domain of the Ukrainian Shield Revealed by Detrital Zircon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Deterministic Physically Based Distributed Models for Rainfall-Induced Shallow Landslides

by
Giada Sannino
1,2,*,
Massimiliano Bordoni
3,
Marco Bittelli
4,
Claudia Meisina
3,
Fausto Tomei
2 and
Roberto Valentino
1
1
Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, 43121 Parma, Italy
2
Regional Agency for Environmental Protection and Energy, Emilia-Romagna Region, 40122 Bologna, Italy
3
Department of Earth and Environmental Sciences, University of Pavia, 27100 Pavia, Italy
4
Department of Agricultural and Food Sciences, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Submission received: 9 August 2024 / Revised: 14 September 2024 / Accepted: 24 September 2024 / Published: 27 September 2024
(This article belongs to the Section Natural Hazards)

Abstract

:
Facing global warming’s consequences is a major issue in the present times. Regarding the climate, projections say that heavy rainfalls are going to increase with high probability together with temperature rise; thus, the hazard related to rainfall-induced shallow landslides will likely increase in density over susceptible territories. Different modeling approaches exist, and many of them are forced to make simplifications in order to reproduce landslide occurrences over space and time. Process-based models can help in quantifying the consequences of heavy rainfall in terms of slope instability at a territory scale. In this study, a narrative review of physically based deterministic distributed models (PBDDMs) is presented. Models were selected based on the adoption of the infinite slope scheme (ISS), the use of a deterministic approach (i.e., input and output are treated as absolute values), and the inclusion of new approaches in modeling slope stability through the ISS. The models are presented in chronological order with the aim of drawing a timeline of the evolution of PBDDMs and providing researchers and practitioners with basic knowledge of what scholars have proposed so far. The results indicate that including vegetation’s effects on slope stability has raised in importance over time but that there is still a need to find an efficient way to include them. In recent years, the literature production seems to be more focused on probabilistic approaches.

1. Introduction

Among landslides, those involving shallow soil (up to a 2–3 m depth) are most frequently induced by rainfall and are highly dangerous, as no premonitory signs are present over territories [1]. The landslide occurrence is the final step of a chain of processes that starts with rainfall infiltration and leads to slope collapse [2]. Therefore, a tool including process quantification to prepare for potential landslide hazards is necessary, especially under climate change conditions; in fact, the temperature rise will allow the atmosphere to hold more moisture so that a greater magnitude of heavy rainfall will likely occur [3]. However, little is known about the effects of climate and its variation on slope stability, landslides, landslide hazards, and the related risk of climate change [4].
Rainfall is recognized as the major shallow landslide triggering factor [5], and over time, different techniques have been employed to predict rainfall-induced shallow landslides in order to constitute Landslide Early Warning System (LEWS) tools [6]. To model landslide occurrences, two main approaches are recognized: physically based modeling and statistical modeling. These models can help in assessing not only the landslide susceptibility of large areas but also the related hazard, which is defined as the detection of the triggering time of a landslide (i.e., “when” or “how frequently” will occur) [7]. Physically based models are also referred to as process-based models, while statistical approaches are further named “data-driven models” as they use past events data based on the assumption that environmental conditions leading to landslides in the past are likely to provoke new instability phenomena in the future [8]. Although they represent interesting techniques for large scale susceptibility analysis, they highly depend on the resolution of past landslides inventories, which can propagate high levels of uncertainty to the outputs if they are not sufficiently detailed. Moreover, under climate change conditions, data-driven approaches have to be undertaken carefully since new environmental and meteorological conditions may not be represented by the past [9].
The use of physical–mathematical models has recently gained large consensus among engineers and the scientific community, not only for applications at the local or slope scale but also at larger extensions, as they can take into account the dynamic variability of the system [10]. Physically based model applications over large areas are also supported by remote sensing, which allows extensive observation of landslides, especially where it would not be possible with ground-based techniques, thus providing robust landslide inventories upon which model output accuracy can be assessed [11,12,13,14]. Whether the aim is to design emergency plans to warn people or inform policy makers about the consequences of extreme climatic events, physically based distributed models represent valid tools as they can detect landslide occurrences in advance. When dealing with future landslide prediction and climate change impacts on slope stability, this modeling approach is considerably noteworthy, especially when vegetation presence is considered [15], thanks to the possibility of assessing single aspects’ contribution to the overall stability.
Physically based deterministic distributed models (PBDDMs) adopt deterministic approaches for both input and output and can be used to back-analyze real landslides in order to derive soil hydraulic and geotechnical parameters. Deterministic models can be highly sensitive to the input parameters’ variability. To reduce uncertainty in assigning unknown soil parameters, some physical models include probabilistic treatments of input data or are coupled with other approaches to obtain probability maps of failure as outputs, e.g., [16].
When developing new models for shallow landslide detection, the processes to consider must be decided. Normally, PBDDMs for shallow landslide prediction are composed of a hydrological module interfaced with a geotechnical module, typically based on the estimation of the slope safety factor through the limit equilibrium method, as it is less computationally demanding than finite element numerical methods, which are often not practically useable over large areas [17]. At the same time, the hydrological processes involved in landslide initiation are many, and including all of them in a single model would make it impossible to run in a reasonable computational time. Because of this complexity, the majority of the existing models that consider hydrological aspects in detail tend to simplify the mechanical aspects. The simplified approach based on the infinite slope scheme (ISS) is adopted in the majority of PBBDMs for rainfall-induced shallow landslides. In this approach, the assumption of a planar failure geometry is considered consistent with the shape, size, and failure mechanisms of rainfall-induced shallow landslides, although this assumption has some applicability limitations [18].
In recent years, models including antecedent hydrological information have raised importance in landslide hazard assessment [19], and including vegetation-related processes has become fundamental in shallow landslides’ prediction. This is because canopies and roots modify the stability condition of slopes under both hydrological and mechanical viewpoints [20,21,22]. The rooted portion of soil has peculiar characteristics, and its behavior, together with proper consideration of unsaturated soil mechanics, should be taken into account in shallow landslide prediction. Furthermore, the interest in nature-based solutions (NBSs) and bioengineering techniques has raised the importance of preventing landslides over territories, not only because of environmental purposes but also because of economic aspects [23,24,25]. Quantitatively assessing the role of vegetation on slope stability is of importance, although a comprehensive methodology is still lacking [26].
In this work, a non-systematic narrative review of the main physically based deterministic distributed models developed over time is presented, highlighting which aspects or techniques were used and introduced over time by the different models. The review is restricted to the models that adopt the infinite slope scheme and are based on the limit equilibrium method. The aim is to draw a chronological memory of rainfall-induced shallow landslides deterministic modeling and also to address new research efforts in directions that have not yet been explored. This paper is organized as follows. In the second section, a brief discussion of shallow landslide modeling methods is presented, stretching some conceptual aspects of hydrology and plants’ contribution. The third section presents the materials and methods adopted for this work, while the fourth section describes the selected models in detail. It was chosen to include only models that do not use any kind of probabilistic approach. This choice is meant to be helpful for practitioners and non-academic, but also for new researchers, in order to classify physically based models for rainfall-induced shallow landslides into two different groups: “purely” physically based deterministic distributed models (PBDDMs), i.e., the specific topic of this review, and hybrid forms. In the end, a conclusive section considers aspects that still have to be properly explored in PBDDMs. Some suggestions are also provided.

2. Shallow Landslide Modeling Methods

2.1. Hydrology

As already mentioned, the majority of the existing PBDDMs for rainfall-induced shallow landslides prediction are composed of two interfaced modules: one computes the hydrological phenomena, and the other computes the slope stability.
A recognized mechanism for shallow landslide activation is the rapid formation of nil or positive pore water pressure because of soil saturation [27,28,29]. Saturated lenses develop in response to rainfall events because of infiltration and seepage processes, and they can trigger landslides when enlarged enough. Water can accumulate in the soil because of the presence of an impermeable or low-permeability layer, either during vertical downward seepage or upward movements [30]. Some landslides can then completely or partially mobilize into destructive earth flows because of liquefaction [31,32,33,34], while others only translate the detached material up to a certain distance.
Because of these observations, the first slope stability models were developed based on saturated soil conditions [35]. Under saturation conditions, soil properties and hydrological parameters can be assumed as constant, and the role of matric suction (i.e., negative pore water pressure) is neglected. Real meteorological conditions can even be ignored in models based on saturated soil mechanics. In the most simplified cases where the soil is considered to be partially saturated, real rainfall amounts are considered as contributing to the uprising of a pre-existing groundwater table.
Some limitations of these simplified approaches have been recognized, and practical applications related to slope stability problems also involve unsaturated soil mechanics [36]. In fact, since shallow landslides have been observed to occur even under negative pore water pressure, i.e., when the soil is in partially saturated condition [37,38], it is important to consider the evolution in time of pore-water pressures and unsaturated soil mechanics in PBDDMs.
Unsaturated soil mechanics are normally based on non-linear soil property functions, including the constitutive relationship between soil water content and soil suction (also defined as the water potential or hydraulic head), namely, the Soil Water Retention Curve (SWRC, or Soil Water Characteristic Curve, SWCC), and the one that relates soil suction or soil water content to soil hydraulic conductivity, namely, the Hydraulic Conductivity Function (HCF) [39,40]. The relationships are unique for a certain pore size distribution of soil and are empirically defined. However, because of hysteretic behavior, these properties follow different patterns during wetting or drying phases [41]. Also, the presence of vegetation was demonstrated to alter the pore size distribution because of root growth, leading to an evolution of the SWRC and of the soil permeability over time with respect to SWRC or HCF of bare soil [42].
When considering unsaturated conditions, transient or stationary input rainfall can be assumed for hydrological balance computation, and water movements can be assessed through a transient analysis. The most commonly used approach to determine the transient values of soil water content is the use of partial differential equations—in particular, numerical solutions of Richards’ equation for unsaturated seepage process [43,44] at different spatial domains (i.e., 1-D, 2-D, or 3-D). This approach allows to consideration of different hydrological processes into a single balance equation, although in the 3-D domain, its computation over large areas can be consuming both in terms of energy and time [45].
When dealing with rainfall-induced landslide hydrology, several processes should be quantified. In any case, the major water input is represented by rainfall infiltration, which enters the soil matrix and undergoes gradient and gravity-driven movements. Then, subsurface flows and redistribution processes take place, and soil hydraulic conductivity anisotropy should be considered [46]. In most cases, different hydraulic behaviors are only considered along soil vertical directions [47].
When a numerical problem is involved, it is essential to properly define boundary conditions, that is, defining how the considered physical system acts at the borders of its spatial and temporal domains [48].
As already mentioned, defining the antecedent condition of soil moisture prior to rainfall is fundamental to correctly simulate water movements during a precipitation event [49]. The global soil water balance equation should be solved for a certain time span, which should be antecedent with respect to the period of interest, in order to assess consistent soil moisture as an initial condition. In this sense, as a key hydrological process, plant growth and evapotranspiration activity should not be neglected [50]. In fact, vegetation alters the soil water content through different mechanisms, including the absorption of water over the rooted portion of the soil; the modification of SWRC; and parameters such as the saturated hydraulic conductivity, canopy rainfall interception, and the creation of preferential flow patterns through roots and stems. This latter aspect is very difficult to consider simultaneously with matrix flow equations if the model allows the consideration of two different domains at the same time [51]. Although preferential flow paths can have different natures and can originate in different ways, any of them is replicable by models based on matrix flow equations only if approaches such as dual-permeability (or double porosity) are adopted [52,53].
Notwithstanding the effectiveness of Richards’ equation in soil water balance computation, this approach may not be the most efficient in specific cases. For example, when a landslide is triggered because of an in-depth wetting front propagation, the hydrological mechanisms can be approximated through simplified approaches for infiltration (e.g., [54,55]).

2.2. Slope Stability

With regard to the geotechnical modules of PBDDMs, the most commonly used approach for slope stability is the computation of a Factor of Safety (FoS). This approach is based on the limit equilibrium method and considers the relation between stabilizing and destabilizing actions on a slope. When the FoS is equal to 1, the slope is in a critical equilibrium state; as the FoS drops below 1, the entire slope is estimated as unstable.
Several methods to compute FoS exist [56]. In a general sense, FoS is computed over one or more potential failure surfaces in order to detect under which conditions and/or at which depth the landslide is triggered. PBDDMs should adopt an automatic procedure for analyzing different slip surfaces, especially when the soil strength and pressure profiles differ along with depth and the model can consider multi-layered soils with different parametrization.
Since shallow landslides are normally translational and their length/depth ratio is generally low, the infinite slope method (ISM), which assumes a ground-parallel planar failure surface, represents the most commonly used approach, although it has limitations [18]. Among others, Lu and Godt developed an equation for FoS that allows the classic saturated soil mechanics theory based on the effective stress concept to be easily extended to the unsaturated regime [57]. However, models that aim to simulate rotational movement also exist, as well as models that can approximate the landslide runout [58].
It is known that vegetation contributes to soil reinforcement in different ways and at different spatial scales [59,60]. With regard to plant roots, it is known that the rooting system can extend the soil shear strength either through water absorption or mechanical reinforcement, normally quantified as an additional cohesive term extending the soil’s effective cohesion [61] and normally applied on the sliding surface. The mechanical reinforcement seems to be more effective in the shallower portion of the soil with respect to the hydrological reinforcement, which is effective down to a 1–2 m depth [62]. More in detail, the mechanical improvement in soil shear strength is exerted through the root network tensile strength and its interaction with soil and bedrock anchoring, as well as the ability of roots to cross the slip surface [63]. Both large and fine roots contribute to the global exerted reinforcement [64]. The overall root cohesion can be derived through different methods, namely, the Wu and Waldron model (WWM), the Fiber Bundle Model (FBM), the Analytical Fiber Bundle Model (AFBM), the Root Bundle Model (RBM), and the root bundle model with root-failure Weibull survival function (RBMw) [25,65,66,67,68,69,70]. The root reinforcement can be either basal or lateral when applied at the base of a failure plane or on the lateral sides of the landslide body, respectively. Nevertheless, in soil conditions close to saturation, it is not totally clear if these effects are still present or not [59,60]. Few of the existing models consider the real spatial variability of root distribution to assess the root reinforcement [71], even if the vegetation stand characteristics affect the magnitude of stabilization effects, especially if gaps are present [72].
PBDDMs can be applied at the slope scale, basin scale, or regional scale. In the first case, it may be possible to carry on field campaigns to obtain specific soil hydraulic and geotechnical parameters. In the second case, especially when the analyzed area is larger than a single slope or a small catchment, field work can become very expensive, and assigning reliable soil parameters can be challenging. For models that do not involve probabilistic approaches, parameters can be assigned based on texture or pedotransfer functions applied to geological units if extended field campaign data are not available [73].

3. Materials and Methods

This study builds upon previous review papers such as [20,25,74] and follows a non-systematic approach, mainly focusing on chronology and some specific aspects of different models. The flow chart of the methodology is represented in Figure 1.
A total of 43 models considered by other reviews were screened, and only 18 were detected as “purely deterministic”. A total of 11 duplicates considered by more than one review were removed, while other 13 models presented some probabilistic approaches, and then they were removed. The resulting 18 deterministic models were analyzed, and 4 of them were discarded as they do not use the infinite slope scheme (ISS) as geotechnical model. In the last screening step, 2 models were removed because they were originally presented as non-distributed models, and 2 models considered by other reviews were an extension of a pre-existing model. The models which constitute extension of their previous version were, however, cited in the following section.

4. Discussion on Models

The 12 selected models are discussed in detail.
SHALSTAB [75] is one of the first GIS-based PBDMs published in 1994. It computes slope stability for cohesionless soils over an infinite slope under a steady-state flow condition, and it is intended to be used with ESRI-ArcGIS. By assuming a steady saturated flow, the model is able to calculate the rainfall amount that is necessary to trigger a landslide over a specified area using a contour-based Digital Elevation Model (DEM) methodology. Its output is an estimation of the critical saturated soil height computed through a hydrological model called TOPOG [76]. The model neglects the effects of the degree of saturation in the vadose zone. Because of its structure and mathematical formulations, SHALSTAB is not suitable to forecast the timing of landslide triggering.
In 1995, the dSLAM model was published [77]. The model aims to quantify the slope instability in steep and forested areas. With dSLAM, the mechanical root reinforcement is introduced in PBDDMs for shallow landslides, while hydraulic vegetation effects are ignored. The model is designed for translational slide overlying a lithic contact. dSLAM is a contour-based model and does not account for rainfall spatial distribution. The mechanical root reinforcement is considered a constant additive term to the overall soil cohesion. Hydrology is based on a kinematic wave groundwater model [78], and it can be run either with actual hyetographs or long-term sequences of precipitation. The authors demonstrate that the spatial distribution of slope instability is controlled by topography and forest harvesting. During major rainstorms, groundwater flow greately affects the factor of safety.
In 1998, the SHETRAN model was extended with a shallow landslide erosion and sediment yield component [79,80]. It is a basin-scale model that considers spatial variability of rainfall input and hydrological responses. It considers snowmelt as a triggering factor as well. The model computes the landslides’ impact on sediment yields at the basin outlet, deriving the volume of dislocated material. Two levels of resolution are comprised: first, a finer spatial assessment of shallow landslide susceptibility is carried out, and a critical soil saturation condition is obtained through GISLIP, a GIS analysis tool; subsequently, a time-varying simulation based on the hydrological grid-based physical SHE model [81] is conducted at a coarser resolution, at which the geotechnical stability analysis also is conducted Therefore, the model can be computationally sustainable for basin-scale simulations (areas around 500 km2). The SHE model comprises evapotranspiration activity and canopy interception.
In 2000, Iverson proposed a model intending to assess the effects of transient rainfall on the timing, rates, and locations of landslides [82]. The hydrology is computed considering the soil infiltration capacity equal to the saturated hydraulic conductivity in order to derive analytical values of the pressure head, assuming a pre-existing steady-state pore pressure. The aim of the model is to derive rainfall thresholds for triggering shallow landslides, corresponding to the peak value of pore pressure, by deriving transient pore pressure distributions that are added to the pre-existing one. Although a simplified transient analysis is included, Iverson’s model is better suitable for shallow landslides related to short-duration rainfall (from approximately 1 h to 70 h), as the model neglects lateral water flow [83]. The model also derives post-failure motion.
TRIGRS [84] is a grid-based PBDDM that aims at locating the timing and size of landslides. At the time, most of the landslide prediction models were producing only susceptibility maps, not involving a complete transient and distributed analysis. In its first version, TRIGRS solves on a pixel-by-pixel basis the one-dimensional vertical version of Richards’ equation by [85], assuming differentiated bare soil areas, which are characterized by a unique homogeneous isotropic layer. In general, TRIGRS is the most commonly used PBDDM for shallow landslides at the regional scale, even if no graphical interface is provided. Its use is linked to a Geographical Information System (GIS), where input data can be prepared and outputs can be visualized. Many versions of TRIGRS have been developed over time, including the unsaturated version, the accounting of a vegetation effect [86], and a parallelized version [87]. TRIGRS provides the minimum Factor of Safety (FoS) calculated at selected time steps of a rainstorm, the pore water pressures, and the depth of the minimum FoS at a certain time step. The soil is described by a SWRC. The unique consideration of 1-D water movements, although accurate under certain conditions, is not appropriate when complex topography is present over a large area. Moreover, the input DEM spatial resolution seems to strongly affect the results [88].
In the same year, 2008, the SLIP model was developed [89,90]. It is a PBDDM that adopts a simplified hydrological approach to simulate soil water balance, assuming that the saturated portion of the soil increases because of infiltrating rainfall. At the same time, the soil desaturates through a percolation process. To be quickly applied to large areas with low computational effort, the model avoids the need to be provided with complex hydrological formulation and approximates a transient seepage computation. In its recent version, the SLIP model also includes a simplified approach for plant interception and root cohesion, changing the model into G-SLIP [91]. SLIP represents a valid example of a simplified although effective solution for Early Warning Systems over large areas, which is an urgent need due to climate change.
In 2011, the SUSHI model was proposed by Capparelli and Versace [92]. It describes water movements in a bi-dimensional domain, allowing the consideration of irregular soil stratigraphy and different soil parameters. A 2-D Richards’ equation is involved in the model, following the assumption of isotropic soil through the adoption of a specific “capillarity coefficient”. This coefficient represents the rate at which water is absorbed or released because of pressure head changes in the soil. A fully implicit method and the finite difference method are employed for hydrology computation in the SUSHI model. Evapotranspiration is accounted for by adopting a uniform root distribution and deriving a sink transpiration term. The vegetation is modeled considering a fixed Leaf Area Index value (LAI) throughout the year.
In 2013, Lepore et al. [93] published a model called tRIBS-VEGGIE-Landslide, then modified by Arnone et al. [94] in 2016 to include probabilistic treatment of uncertainties. The methodology is based on a Triangular Irregular Network (TIN) mesh and accounts for post-failure movement by considering selected slope angles as thresholds for determining whether the landslide body will move down to a run-out distance. Landslide movement is assumed to follow the same flow directions evaluated by the tRIBS hydrological component (the deepest descend), which is based on a transient computation of infiltration and redistribution processes. The basin morphology can be modified by landslide deposits, with consequent impact on most of the simulated processes. Concerning vegetation, the model considers root mechanical reinforcement together with general soil cohesion, while the transpiration is estimated through the methodology provided by [95], which is based on vapor pressure deficit, soil moisture levels, rooting profile, leaf area, and available energy.
In 2014, Milledge et al. [96] pointed out that the existing models were highly computationally demanding and were not practically applicable across landscapes. The proposed model, called MD-STAB, simulates lateral resistances acting on landslide margins using earth pressure theory and the lateral root distribution, which is modeled as an exponential function of soil depth in a three-dimensional limit equilibrium force balance. This assumption allows the model to consider roots that laterally cross the shearing surfaces. The possibility of considering forces that act on the lateral sides of rigid blocks contradicts the infinite slope assumption, for which the inter-slice interactions are ignored. This model has also been extended, including the derivation of root reinforcement from field-measured forest stand characteristics [97]. The model ignores infiltration, soil suction, and capillary rises, and the groundwater level is steady and parallel to the slope surface.
In 2017, SOSlope was published by Cohen and Schwarz [98]. The model focuses on the effect of root and soil strength on slope stability in forests, accounting for single-tree contribution. The hydrological aspects are considered through a simplified and empirical dual-porosity model. Through this approach, SOSlope can approximate the water dynamics in both the soil matrix and the preferential flow domains. The model is suitable for assessing fundamental aspects such as the role of the forest structure (e.g., tree size, tree spacing), root distribution, and root mechanical properties on the triggering mechanisms of shallow landslides. SOSlope considers both lateral and basal root reinforcement and is able to reproduce the self-organized redistribution of forces on a slope during rainfall-triggered shallow landslides. The model is particularly suitable for highly detailed forest management purposes, and outputs can be used in GIS environments.
Lizàrraga and Buscarnera, in 2018, developed a model that uses suction-dependent plasticity and limit equilibrium theories to derive the slope Factor of Safety (FoS) in unsaturated soils [99,100]. The model considers mechanical aspects usually not included in other models, such as suction-hardening and liquefaction potential [101]. The aim is to simultaneously quantify the susceptibility to frictional slips and liquefaction-induced flow slides of shallow soil slopes in order to incorporate these considerations in regional-scale landslide hazard mapping. The underlying hydrology is a transient computation based on Richards’ equation. Laboratory data are used to determine input parameters. The model application points out the strong interplay between infiltration mechanisms (i.e., slow or fast) and the mode and depth of slope instability. Vegetation effects are neglected.
In 2023, Abdollhai et al. [102] proposed a model specifically designed for the estimation of hillslope post-wildfire stability against rainfall-induced landslides. The proposed model constitutes a physically based yet practical slope stability framework capable of capturing the interplay of key driving factors and wildfire-induced alterations. The aim is to derive post-wildfire temporal changes in the Factor of Safety in response to rainfall. The transient hydrological analysis is conducted through an analytical approximation of 1D Richards’ equation. The soil is considered isotropic and homogeneous. Changes in saturated soil water content are determined as a function of the soil elastic modulus, and the method proposed by [103] is then used to derive saturated water content in deformable soils. The post-wildfire effects on vegetation are quantified as a reduction in the transpiration rate due to the decrease in the number of roots. The infiltration capacity is also affected by wildfire, although in the first application, the authors decided not to account for it.
Figure 2 reports the normalized number of citations of the different considered models. The normalized number of citations is calculated as the ratio between the total number of citations (as derived by Google Scholar) and the number of years since publication (i.e., the difference between the current year—2024—and the year of publication). It can be seen that the most cited models are the Iverson model, SHALSTAB, and TRIGRS. Out of them, only the dSLAM model has a normalized citation number higher than 25. These four models can thus be considered as milestones of PBDDM knowledge and development over time.
Table 1 summarizes some features of the 12 analyzed models. Notwithstanding the adoption of the infinite slope scheme, the hydrology modules, as highlighted, show that the use of the 1D Richards equation in deterministic models for shallow landslides is the most commonly used method.

5. Discussion and Conclusions

Based on the publishing date, several models were described in this narrative review, trying to draw a comprehensive overview of physically based deterministic models for shallow landslides.
By analyzing the literature, purely deterministic approaches appear to have recently become less explored than they were in recent decades. However, it appears clear that the attention given to the stabilizing effects of vegetation has grown over the years, as testified by the rising complexity involved when trying to assess these effects in a consistent way. Big challenges are still open about this topic since most of the root reinforcement models are based on variables and parameters not easily derivable such as root architecture or root diameters. Most of the root reinforcement models have a strong empirical basis that is not easily applicable at large scale, particularly when different kinds of plant species cohabit [104]. New paradigms and expedients should be explored for root effect quantification over large areas, and purely deterministic models can help in this task, as the processes and the related parameters can be quantified and studied singularly in a specific way. An example is provided by [105], where root-induced modifications of soil hydraulic properties (namely, the saturated hydraulic conductivity and the Soil Water Retention Curve) are included in a physically based model. From this point of view, involving vegetation as a dynamic variable, accounting for growing over time and the space of roots and canopy, based on real meteorological data and intraseasonal dynamic variability, can improve the performance of slope stability models [106,107].
More studies based on remote sensing linking canopy development with root architectures, rooting depth, and root spatial extensions could help in applying root development and reinforcement models over large areas, as vegetation parameters can be more easily derived when related to aboveground biomass [61]. On the contrary, non-invasive techniques for estimating roots’ morphometric characteristics are still difficult to use.
An important aspect that constitutes a strong point since the first PBDMs were published is the interoperability with GIS environments. This aspect can help decision makers to cross different sources of information about territories. As it is difficult to have a unique, comprehensive PBDDM for all the processes that should be accounted for in landslide hazard assessment, overlaying different model outputs or different spatial information can provide more insights on a large scale. In fact, it should not be overlooked that PBDDMs pretend to give a single, specific output based on single, selected input parameters that may not represent reality, especially when large periods are considered. Comparing different sources of information through pre- or post-processing techniques and procedures may help raise the reliability of landslide risk assessment analyses, maintaining an acceptable operational time.
In distributed models, sensitivity analysis to different DEM (or mesh) spatial resolutions should be included. It is known that accuracy changes according to the considered spatial resolution, either when computing hydrological phenomena or slope stability, especially when this latter is based on the infinite slope method. This aspect appears not properly considered when new models are developed, although the problem is crucial for proper application of the model itself. In fact, spatial discretization of the domain can lead to different results.
Although rarely discussed when new models are developed, an important aspect for practical applications is the required running time of the different algorithms at the scale of interest. This problem is crucial, especially if the model should be used in early warning systems for civil protection purposes at the regional scale.
It is worth remembering that a very important expedient to overcome the spatial uncertainty of parameters is including probabilistic approaches in PBDMs, thus leading to hybrid solutions. There are several valid examples in the literature, such as SINMAP [108], GEOTOP-Fs [109], HIRESSS [110], SlideForMap [111], and FSLAM [112] These models include probabilistic assessments at different complexity rates, but in this paper, it was preferred to focus on physically based distributed models characterized by only deterministic input parameters and deterministic output results.

Author Contributions

Conceptualization, G.S.; methodology, G.S. and R.V.; resources, G.S., M.B. (Massimiliano Bordoni) and R.V.; writing—original draft preparation, G.S.; writing—review and editing, G.S., M.B. (Marco Bittelli), C.M. and F.T.; supervision, R.V. and M.B. (Massimiliano Bordoni). All authors have read and agreed to the published version of the manuscript.

Funding

This research is co-funded by the Italian Ministry of University and Research (MUR) and the European Union through the REACT-EU program.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Persichillo, M.G.; Bordoni, M.; Meisina, C. The role of land use changes in the distribution of shallow landslides. Sci. Total Environ. 2017, 574, 924–937. [Google Scholar] [CrossRef]
  2. Cotecchia, F.; Santaloia, F.; Lollino, P.; Vitone, C.; Pedone, G.; Bottiglieri, O. From a phenomenological to a geomechanical approach to landslide hazard analysis. Eur. J. Environ. Civ. Eng. 2016, 20, 1004–1031. [Google Scholar] [CrossRef]
  3. Crozier, M.J. Deciphering the effect of climate change on landslide activity: A review. Geomorphology 2010, 124, 260–267. [Google Scholar] [CrossRef]
  4. Gariano, S.L.; Guzzetti, F. Landslides in a changing climate. Earth-Sci. Rev. 2016, 162, 227–252. [Google Scholar] [CrossRef]
  5. Polemio, M.; Petrucci, O. Rainfall as a Landslide Triggering Factor an Overview of Recent International Research; Thomas Telford Ltd.: London, UK, 1999; pp. 1219–1226. [Google Scholar]
  6. Guzzetti, F.; Gariano, S.L.; Peruccacci, S.; Brunetti, M.T.; Marchesini, I.; Rossi, M.; Melillo, M. Geographical landslide early warning systems. Earth-Sci. Rev. 2020, 200, 102973. [Google Scholar] [CrossRef]
  7. Reichenbach, P.; Rossi, M.; Malamud, B.D.; Mihir, M.; Guzzetti, F. A review of statistically-based landslide susceptibility models. Earth-Sci. Rev. 2018, 180, 60–91. [Google Scholar] [CrossRef]
  8. Guzzetti, F.; Carrara, A.; Cardinali, M.; Reichenbach, P. Landslide hazard evaluation: A review of current techniques and their application in a multi-scale study, Central Italy. Geomorphology 1999, 31, 181–216. [Google Scholar] [CrossRef]
  9. Fang, Z.; Wang, Y.; van Westen, C.; Lombardo, L. Landslide hazard spatiotemporal prediction based on data-driven models: Estimating where, when and how large landslide may be. Int. J. Appl. Earth Obs. Geoinf. 2024, 126, 103631. [Google Scholar] [CrossRef]
  10. Crosta, G.B.; Frattini, P.J.N.H. Distributed modelling of shallow landslides triggered by intense rainfall. Nat. Hazards Earth Syst. Sci. 2003, 3, 81–93. [Google Scholar] [CrossRef]
  11. Kyriou, A.; Nikolakopoulos, K.G.; Koukouvelas, I.K. Timely and low-cost remote sensing practices for the assessment of landslide activity in the service of hazard management. Remote Sens. 2022, 14, 4745. [Google Scholar] [CrossRef]
  12. Nikolakopoulos, K.G.; Kyriou, A.; Koukouvelas, I.K.; Tomaras, N.; Lyros, E. UAV, GNSS, and InSAR Data Analyses for Landslide Monitoring in a Mountainous Village in Western Greece. Remote Sens. 2023, 15, 2870. [Google Scholar] [CrossRef]
  13. Nikolakopoulos, K.G.; Kyriou, A.; Koukouvelas, I.K. UAV, GNSS, and GIS for the Rapid Assessment of Multi-Occurrence Landslides. Geosciences 2024, 14, 160. [Google Scholar] [CrossRef]
  14. Tsironi, V.; Ganas, A.; Karamitros, I.; Efstathiou, E.; Koukouvelas, I.; Sokos, E. Kinematics of active landslides in Achaia (Peloponnese, Greece) through InSAR time series analysis and relation to rainfall patterns. Remote Sens. 2022, 14, 844. [Google Scholar] [CrossRef]
  15. Li, Y.; Duan, W. Decoding vegetation’s role in landslide susceptibility mapping: An integrated review of techniques and future directions. Biogeotechnics 2023, 2, 100056. [Google Scholar] [CrossRef]
  16. Raia, S.; Alvioli, M.; Rossi, M.; Baum, R.L.; Godt, J.W.; Guzzetti, F. Improving predictive power of physically based rainfall-induced shallow landslide models: A probabilistic approach. Geosci. Model Dev. 2014, 7, 495–514. [Google Scholar] [CrossRef]
  17. Scaringi, G.; Fan, X.; Xu, Q.; Liu, C.; Ouyang, C.; Domènech, G.; Yang, F.; Dai, L. Some considerations on the use of numerical methods to simulate past landslides and possible new failures: The case of the recent Xinmo landslide (Sichuan, China). Landslides 2018, 15, 1359–1375. [Google Scholar] [CrossRef]
  18. Milledge, D.G.; Griffiths, D.V.; Lane, S.N.; Warburton, J. Limits on the validity of infinite length assumptions for modelling shallow landslides. Earth Surf. Process. Landf. 2012, 37, 1158–1166. [Google Scholar] [CrossRef]
  19. Greco, R.; Marino, P.; Bogaard, T.A. Recent advancements of landslide hydrology. Wiley Interdiscip. Rev. Water 2023, 10, e1675. [Google Scholar] [CrossRef]
  20. Murgia, I.; Giadrossich, F.; Mao, Z.; Cohen, D.; Capra, G.F.; Schwarz, M. Modeling shallow landslides and root reinforcement: A review. Ecol. Eng. 2022, 181, 106671. [Google Scholar] [CrossRef]
  21. Lu, J.; Zhang, Q.; Werner, A.D.; Li, Y.; Jiang, S.; Tan, Z. Root-induced changes of soil hydraulic properties—A review. J. Hydrol. 2020, 589, 125203. [Google Scholar] [CrossRef]
  22. Ni, J.J.; Leung, A.K.; Ng, C.W.W. Modelling effects of root growth and decay on soil water retention and permeability. Can. Geotech. J. 2019, 56, 1049–1055. [Google Scholar] [CrossRef]
  23. de Jesús Arce-Mojica, T.; Nehren, U.; Sudmeier-Rieux, K.; Miranda, P.J.; Anhuf, D. Nature-based solutions (NbS) for reducing the risk of shallow landslides: Where do we stand? Int. J. Disaster Risk Reduct. 2019, 41, 101293. [Google Scholar] [CrossRef]
  24. Vicarelli, M.; Sudmeier-Rieux, K.; Alsadadi, A.; Shrestha, A.; Schütze, S.; Kang, M.M.; Leue, M.; Wasielewski, D.; Mysiak, J. On the cost-effectiveness of Nature-based Solutions for reducing disaster risk. Sci. Total Environ. 2024, 947, 174524. [Google Scholar] [CrossRef] [PubMed]
  25. DiBiagio, A.; Capobianco, V.; Oen, A.; Tallaksen, L.M. State-of-the-art: Parametrization of hydrological and mechanical reinforcement effects of vegetation in slope stability models for shallow landslides. Landslides 2024, 21, 2417–2446. [Google Scholar] [CrossRef]
  26. Mao, Z. Root reinforcement models: Classification, criticism and perspectives. Plant Soil 2022, 472, 17–28. [Google Scholar] [CrossRef]
  27. Vanapalli, S.K.; Fredlund, D.G.; Pufahl, D.E.; Clifton, A.W. Model for the prediction of shear strength with respect to soil suction. Can. Geotech. J. 1996, 33, 379–392. [Google Scholar] [CrossRef]
  28. Godt, J.W.; Baum, R.L.; Savage, W.Z.; Salciarini, D.; Schulz, W.H.; Harp, E.L. Transient deterministic shallow landslide modeling: Requirements for susceptibility and hazard assessments in a GIS framework. Eng. Geol. 2008, 102, 214–226. [Google Scholar] [CrossRef]
  29. Kuriakose, S.L.; Jetten, V.G.; Van Westen, C.J.; Sankar, G.; Van Beek, L.P.H. Pore water pressure as a trigger of shallow landslides in the Western Ghats of Kerala, India: Some preliminary observations from an experimental catchment. Phys. Geogr. 2008, 29, 374–386. [Google Scholar] [CrossRef]
  30. Ran, Q.; Hong, Y.; Li, W.; Gao, J. A modelling study of rainfall-induced shallow landslide mechanisms under different rainfall characteristics. J. Hydrol. 2018, 563, 790–801. [Google Scholar] [CrossRef]
  31. McKenna, J.P.; Santi, P.M.; Amblard, X.; Negri, J. Effects of soil-engineering properties on the failure mode of shallow landslides. Landslides 2012, 9, 215–228. [Google Scholar] [CrossRef]
  32. Gofar, N.; Rahardjo, H. Saturated and unsaturated stability analysis of slope subjected to rainfall infiltration. MATEC Web Conf. 2017, 101, 05004. [Google Scholar] [CrossRef]
  33. Gabet, E.J.; Mudd, S.M. The mobilization of debris flows from shallow landslides. Geomorphology 2006, 74, 207–218. [Google Scholar] [CrossRef]
  34. Iverson, R.M.; Reid, M.E.; LaHusen, R.G. Debris-flow mobilization from landslides. Annu. Rev. Earth Planet. Sci. 1997, 25, 85–138. [Google Scholar] [CrossRef]
  35. Ridley, A.M. Soil suction—What it is and how to successfully measure it. In Proceedings of the FMGM 2015: Ninth Symposium on Field Measurements in Geomechanics, Sydney, Australia, 9–11 September 2015; Australian Centre for Geomechanics: Crawley, Australia, 2015; pp. 27–46. [Google Scholar]
  36. Fredlund, D.G. Unsaturated soil mechanics in engineering practice. J. Geotech. Geoenvironmental Eng. 2006, 132, 286–321. [Google Scholar] [CrossRef]
  37. Bittelli, M.; Valentino, R.; Salvatorelli, F.; Pisa, P.R. Monitoring soil-water and displacement conditions leading to landslide occurrence in partially saturated clays. Geomorphology 2012, 173–174, 161–173. [Google Scholar] [CrossRef]
  38. Godt, J.W.; Baum, R.L.; Lu, N. Landsliding in partially saturated materials. Geophys. Res. Lett. 2009, 36, L02403. [Google Scholar] [CrossRef]
  39. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef]
  40. Van Genuchten, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef]
  41. Bordoni, M.; Meisina, C.; Valentino, R.; Lu, N.; Bittelli, M.; Chersich, S. Hydrological factors affecting rainfall-induced shallow landslides: From the field monitoring to a simplified slope stability analysis. Eng. Geol. 2015, 193, 19–37. [Google Scholar] [CrossRef]
  42. Leung, A.K.; Garg, A.; Ng, C.W.W. Effects of plant roots on soil-water retention and induced suction in vegetated soil. Eng. Geol. 2015, 193, 183–197. [Google Scholar] [CrossRef]
  43. Farthing, M.W.; Ogden, F.L. Numerical solution of Richards’ equation: A review of advances and challenges. Soil Sci. Soc. Am. J. 2017, 81, 1257–1269. [Google Scholar] [CrossRef]
  44. Zha, Y.; Yang, J.; Zeng, J.; Tso, C.H.M.; Zeng, W.; Shi, L. Review of numerical solution of Richardson–Richards equation for variably saturated flow in soils. Wiley Interdiscip. Rev. Water 2019, 6, e1364. [Google Scholar] [CrossRef]
  45. Cockett, R.; Heagy, L.J.; Haber, E. Efficient 3D inversions using the Richards equation. Comput. Geosci. 2018, 116, 91–102. [Google Scholar] [CrossRef]
  46. Priono; Rahardjo, H.; Chatterjea, K.; Leong, E.C.; Wang, J.Y. Effect of hydraulic anisotropy on soil–water characteristic curve. Soils Found. 2016, 56, 228–239. [Google Scholar] [CrossRef]
  47. Jing, Y.S.; Zhang, B.; Thimm, A.; Zepp, H. Anisotropy of soil hydraulic properties along arable slopes. Pedosphere 2008, 18, 353–362. [Google Scholar] [CrossRef]
  48. Ali, A.; Huang, J.; Lyamin, A.V.; Sloan, S.W.; Cassidy, M.J. Boundary effects of rainfall-induced landslides. Comput. Geotech. 2014, 61, 341–354. [Google Scholar] [CrossRef]
  49. Liu, H.; Lei, T.W.; Zhao, J.; Yuan, C.P.; Fan, Y.T.; Qu, L.Q. Effects of rainfall intensity and antecedent soil water content on soil infiltrability under rainfall conditions using the run off-on-out method. J. Hydrol. 2011, 396, 24–32. [Google Scholar] [CrossRef]
  50. Tabacchi, E.; Lambs, L.; Guilloy, H.; Planty-Tabacchi, A.M.; Muller, E.; Decamps, H. Impacts of riparian vegetation on hydrological processes. Hydrol. Process. 2000, 14, 2959–2976. [Google Scholar] [CrossRef]
  51. Gerke, H.H. Preferential flow descriptions for structured soils. J. Plant Nutr. Soil Sci. 2006, 169, 382–400. [Google Scholar] [CrossRef]
  52. Montrasio, L.; Valentino, R.; Losi, G.L. Rainfall infiltration in a shallow soil: A numerical simulation of the double-porosity effect. Electron. J. Geotechnol. Eng. 2011, 16, 1387–1403. [Google Scholar]
  53. Shao, W.; Yang, Z.; Ni, J.; Su, Y.; Nie, W.; Ma, X. Comparison of single-and dual-permeability models in simulating the unsaturated hydro-mechanical behavior in a rainfall-triggered landslide. Landslides 2018, 15, 2449–2464. [Google Scholar] [CrossRef]
  54. Rawls, W.J.; Brakensiek, D.L.; Miller, N. Green-Ampt infiltration parameters from soils data. J. Hydraul. Eng. 1983, 109, 62–70. [Google Scholar] [CrossRef]
  55. Cho, S.E. Prediction of shallow landslide by surficial stability analysis considering rainfall infiltration. Eng. Geol. 2017, 231, 126–138. [Google Scholar] [CrossRef]
  56. Xiao, S.; Dai, T.; Li, S. Review and Comparative Analysis of Factor of Safety Definitions in Slope Stability. Geotech. Geol. Eng. 2024, 42, 4263–4283. [Google Scholar] [CrossRef]
  57. Lu, N.; Godt, J. Infinite slope stability under steady unsaturated seepage conditions. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef]
  58. Liu, W.; Yang, Z.; He, S. Modeling the landslide-generated debris flow from formation to propagation and run-out by considering the effect of vegetation. Landslides 2021, 18, 43–58. [Google Scholar] [CrossRef]
  59. Masi, E.B.; Tofani, V.; Rossi, G.; Cuomo, S.; Wu, W.; Salciarini, D.; Caporali, E.; Catani, F. Effects of roots cohesion on regional distributed slope stability modelling. Catena 2023, 222, 106853. [Google Scholar] [CrossRef]
  60. Masi, E.B.; Segoni, S.; Tofani, V. Root reinforcement in slope stability models: A review. Geosciences 2021, 11, 212. [Google Scholar] [CrossRef]
  61. Bischetti, G.B.; Chiaradia, E.A.; Epis, T.; Morlotti, E. Root cohesion of forest species in the Italian Alps. Plant Soil 2009, 324, 71–89. [Google Scholar] [CrossRef]
  62. Ni, J.J.; Leung, A.K.; Ng, C.W.W.; Shao, W. Modelling hydro-mechanical reinforcements of plants to slope stability. Comput. Geotech. 2018, 95, 99–109. [Google Scholar] [CrossRef]
  63. Schmaltz, E.M.; Mergili, M. Integration of root systems into a GIS-based slip surface model: Computational experiments in a generic hillslope environment. Landslides 2018, 15, 1561–1575. [Google Scholar] [CrossRef]
  64. Giadrossich, F.; Cohen, D.; Schwarz, M.; Ganga, A.; Marrosu, R.; Pirastru, M.; Capra, G.F. Large roots dominate the contribution of trees to slope stability. Earth Surf. Process. Landf. 2019, 44, 1602–1609. [Google Scholar] [CrossRef]
  65. Waldron, L.J. The shear resistance of root-permeated homogeneous and stratified soil. Soil Sci. Soc. Am. J. 1977, 41, 843–849. [Google Scholar] [CrossRef]
  66. Wu, T.H.; McKinnell, W.P., III; Swanston, D.N. Strength of tree roots and landslides on Prince of Wales Island, Alaska. Can. Geotech. J. 1979, 16, 19–33. [Google Scholar] [CrossRef]
  67. Pollen, N.; Simon, A. Estimating the mechanical effects of riparian vegetation on stream bank stability using a fiber bundle model. Water Resour. Res. 2005, 41. [Google Scholar] [CrossRef]
  68. Cohen, D.; Schwarz, M.; Or, D. An analytical fiber bundle model for pullout mechanics of root bundles. J. Geophys. Res. Earth Surf. 2011, 116. [Google Scholar] [CrossRef]
  69. Schwarz, M.; Cohen, D.; Or, D. Root-soil mechanical interactions during pullout and failure of root bundles. J. Geophys. Res. Earth Surf. 2010, 115. [Google Scholar] [CrossRef]
  70. Schwarz, M.; Giadrossich, F.; Cohen, D. Modeling root reinforcement using a root-failure Weibull survival function. Hydrol. Earth Syst. Sci. 2013, 17, 4367–4377. [Google Scholar] [CrossRef]
  71. Schwarz, M.; Lehmann, P.; Or, D. Quantifying lateral root reinforcement in steep slopes–from a bundle of roots to tree stands. Earth Surf. Process. Landf. 2010, 35, 354–367. [Google Scholar] [CrossRef]
  72. Moos, C.; Bebi, P.; Graf, F.; Mattli, J.; Rickli, C.; Schwarz, M. How does forest structure affect root reinforcement and susceptibility to shallow landslides? Earth Surf. Process. Landf. 2016, 41, 951–960. [Google Scholar] [CrossRef]
  73. Patil, N.G.; Singh, S.K. Pedotransfer functions for estimating soil hydraulic properties: A review. Pedosphere 2016, 26, 417–430. [Google Scholar] [CrossRef]
  74. Meena, V.; Kumari, S.P.; Shankar, V. Physically based modelling techniques for landslide susceptibility analysis: A comparison. IOP Conf. Ser.: Earth Environ. Sci. 2022, 1032, 012033. [Google Scholar] [CrossRef]
  75. Montgomery, D.R.; Dietrich, W.E. A physically based model for the topographic control on shallow landsliding. Water Resour. Res. 1994, 30, 1153–1171. [Google Scholar] [CrossRef]
  76. O’Loughlin, E.M. Prediction of surface saturation zones in natural catchments by topographic analysis. Water Resour. Res. 1986, 22, 794–804. [Google Scholar] [CrossRef]
  77. Wu, W.; Sidle, R.C. A distributed slope stability model for steep forested basins. Water Resour. Res. 1995, 31, 2097–2110. [Google Scholar] [CrossRef]
  78. Takasao, T.; Shiiba, M. Incorporation of the effect of concentration of flow into the kinematic wave equations and its applications to runoff system lumping. J. Hydrol. 1988, 102, 301–322. [Google Scholar] [CrossRef]
  79. Burton, A.; Bathurst, J.C. Physically based modelling of shallow landslide sediment yield at a catchment scale. Environ. Geol. 1998, 35, 89–99. [Google Scholar] [CrossRef]
  80. Ewen, J.; Parkin, G.; O’Connell, P.E. SHETRAN: Distributed river basin flow and transport modeling system. J. Hydrol. Eng. 2000, 5, 250–258. [Google Scholar] [CrossRef]
  81. Abbott, M.B.; Bathurst, J.C.; Cunge, J.A.; O’Connell, P.E.; Rasmussen, J. An introduction to the European Hydrological System—Systeme Hydrologique Europeen,“SHE”, 1: History and philosophy of a physically-based, distributed modelling system. J. Hydrol. 1986, 87, 45–59. [Google Scholar] [CrossRef]
  82. Iverson, R.M. Landslide triggering by rain infiltration. Water Resour. Res. 2000, 36, 1897–1910. [Google Scholar] [CrossRef]
  83. Frattini, P.; Crosta, G.; Sosio, R. Approaches for defining thresholds and return periods for rainfall-triggered shallow landslides. Hydrol. Process. Int. J. 2009, 23, 1444–1460. [Google Scholar] [CrossRef]
  84. Baum, R.L.; Savage, W.Z.; Godt, J.W. TRIGRS: A Fortran Program for Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Analysis, Version 2.0, Open File Report 2008-1159; US Geological Survey: Reston, VA, USA, 2008; p. 75.
  85. Srivastava, R.; Yeh, T.C.J. Analytical solutions for one-dimensional, transient infiltration toward the water table in homogeneous and layered soils. Water Resour. Res. 1991, 27, 753–762. [Google Scholar] [CrossRef]
  86. Saadatkhah, N.; Tehrani, M.H.; Mansor, S.; Khuzaimah, Z.; Kassim, A.; Saadatkhah, R. Impact assessment of land cover changes on the runoff changes on the extreme flood events in the Kelantan River basin. Arab. J. Geosci. 2016, 9, 687. [Google Scholar] [CrossRef]
  87. Alvioli, M.; Baum, R.L. Parallelization of the TRIGRS model for rainfall-induced landslides using the message passing interface. Environ. Model. Softw. 2016, 81, 122–135. [Google Scholar] [CrossRef]
  88. Viet, T.T.; Lee, G.; Thu, T.M.; An, H.U. Effect of digital elevation model resolution on shallow landslide modeling using TRIGRS. Nat. Hazards Rev. 2017, 18, 04016011. [Google Scholar] [CrossRef]
  89. Montrasio, L.; Valentino, R. A model for triggering mechanisms of shallow landslides. Nat. Hazards Earth Syst. Sci. 2008, 8, 1149–1159. [Google Scholar] [CrossRef]
  90. Montrasio, L.; Valentino, R.; Corina, A.; Rossi, L.; Rudari, R. A prototype system for space–time assessment of rainfall-induced shallow landslides in Italy. Nat. Hazards 2014, 74, 1263–1290. [Google Scholar] [CrossRef]
  91. Montrasio, L.; Gatto, M.P.A.; Miodini, C. The role of plants in the prevention of soil-slip: The G-SLIP model and its application on territorial scale through G-XSLIP platform. Landslides 2023, 20, 1149–1165. [Google Scholar] [CrossRef]
  92. Capparelli, G.; Versace, P. FLaIR and SUSHI: Two mathematical models for early warning of landslides induced by rainfall. Landslides 2011, 8, 67–79. [Google Scholar] [CrossRef]
  93. Lepore, C.; Arnone, E.; Noto, L.V.; Sivandran, G.; Bras, R.L. Physically based modeling of rainfall-triggered landslides: A case study in the Luquillo forest, Puerto Rico. Hydrol. Earth Syst. Sci. 2013, 17, 3371–3387. [Google Scholar] [CrossRef]
  94. Arnone, E.; Caracciolo, D.; Noto, L.V.; Preti, F.; Bras, R.L. Modeling the hydrological and mechanical effect of roots on shallow landslides. Water Resour. Res. 2016, 52, 8590–8612. [Google Scholar] [CrossRef]
  95. Ivanov, V.Y.; Bras, R.L.; Vivoni, E.R. Vegetation-hydrology dynamics in complex terrain of semiarid areas: 2. Energy-water controls of vegetation spatiotemporal dynamics and topographic niches of favorability. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef]
  96. Milledge, D.G.; Bellugi, D.; McKean, J.A.; Densmore, A.L.; Dietrich, W.E. A multidimensional stability model for predicting shallow landslide size and shape across landscapes. J. Geophys. Res. Earth Surf. 2014, 119, 2481–2504. [Google Scholar] [CrossRef] [PubMed]
  97. Cislaghi, A.; Chiaradia, E.A.; Bischetti, G.B. Including root reinforcement variability in a probabilistic 3D stability model. Earth Surf. Process. Landf. 2017, 42, 1789–1806. [Google Scholar] [CrossRef]
  98. Cohen, D.; Schwarz, M. Tree-root control of shallow landslides. Earth Surf. Dyn. 2017, 5, 451–477. [Google Scholar] [CrossRef]
  99. Lizárraga, J.J.; Frattini, P.; Crosta, G.B.; Buscarnera, G. Regional-scale modelling of shallow landslides with different initiation mechanisms: Sliding versus liquefaction. Eng. Geol. 2018, 228, 346–356. [Google Scholar] [CrossRef]
  100. Buscarnera, G.; Di Prisco, C. Soil stability and flow slides in unsaturated shallow slopes: Can saturation events trigger liquefaction processes? Géotechnique 2013, 63, 801–817. [Google Scholar] [CrossRef]
  101. Lizárraga, J.J.; Buscarnera, G. Safety factors to detect flowslides and slips in unsaturated shallow slopes. Géotechnique 2018, 68, 442–450. [Google Scholar] [CrossRef]
  102. Abdollahi, M.; Vahedifard, F.; Tracy, F.T. Post-wildfire stability of unsaturated hillslopes against rainfall-triggered landslides. Earth’s Future 2023, 11, e2022EF003213. [Google Scholar] [CrossRef]
  103. Tracy, F.T.; Vahedifard, F. Analytical solution for coupled hydro-mechanical modeling of infiltration in unsaturated soils. J. Hydrol. 2022, 612, 128198. [Google Scholar] [CrossRef]
  104. Glenn, E.P.; Huete, A.R.; Nagler, P.L.; Nelson, S.G. Relationship between remotely-sensed vegetation indices, canopy attributes and plant physiological processes: What vegetation indices can and cannot tell us about the landscape. Sensors 2008, 8, 2136–2160. [Google Scholar] [CrossRef]
  105. Lann, T.; Bao, H.; Lan, H.; Zheng, H.; Yan, C.; Peng, J. Hydro-mechanical effects of vegetation on slope stability: A review. Sci. Total Environ. 2024, 926, 171691. [Google Scholar] [CrossRef] [PubMed]
  106. Schmaltz, E.M.; Van Beek, L.P.H.; Bogaard, T.A.; Kraushaar, S.; Steger, S.; Glade, T. Strategies to improve the explanatory power of a dynamic slope stability model by enhancing land cover parameterisation and model complexity. Earth Surf. Process. Landf. 2019, 44, 1259–1273. [Google Scholar] [CrossRef]
  107. Sannino, G.; Tomei, F.; Bittelli, M.; Bordoni, M.; Meisina, C.; Valentino, R. Implementation of a slope stability method in the CRITERIA-1D agro-hydrological modeling scheme. Landslides 2024, 1–20. [Google Scholar] [CrossRef]
  108. Pack, R.T.; Tarboton, D.G.; Goodwin, C.N. The SINMAP approach to terrain stability mapping. In Proceedings of the 8th Congress of the International Association of Engineering Geology, Vancouver, BC, Canada, 21–25 September 1998. [Google Scholar]
  109. Simoni, S.; Zanotti, F.; Bertoldi, G.; Rigon, R. Modelling the probability of occurrence of shallow landslides and channelized debris flows using GEOtop-FS. Hydrol. Process. 2008, 22, 532–545. [Google Scholar] [CrossRef]
  110. Rossi, G.; Catani, F.; Leoni, L.; Segoni, S.; Tofani, V. HIRESSS: A physically based slope stability simulator for HPC applications. Nat. Hazards Earth Syst. Sci. 2013, 13, 151–166. [Google Scholar] [CrossRef]
  111. van Zadelhoff, F.B.; Albaba, A.; Cohen, D.; Phillips, C.; Schaefli, B.; Dorren, L.K.A.; Schwarz, M. Introducing SlideforMAP; a probabilistic finite slope approach for modelling shallow landslide probability in forested situations. Nat. Hazards Earth Syst. Sci. 2021, 22, 2611–2635. [Google Scholar] [CrossRef]
  112. Medina, V.; Hürlimann, M.; Guo, Z.; Lloret, A.; Vaunat, J. Fast physically-based model for rainfall-induced landslide susceptibility assessment at regional scale. Catena 2021, 201, 105213. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the methodology. * One of the reviewed models, SOSlope, uses the Discrete Element Method, but the ISS was adopted for the geometry definition. ** With “similar approaches”, it is intended that the model was extending its previous version.
Figure 1. Flow chart of the methodology. * One of the reviewed models, SOSlope, uses the Discrete Element Method, but the ISS was adopted for the geometry definition. ** With “similar approaches”, it is intended that the model was extending its previous version.
Geosciences 14 00255 g001
Figure 2. PBDDMs’ normalized number of citations obtained (derived by Google Scholar). Models are in chronological order [99,102].
Figure 2. PBDDMs’ normalized number of citations obtained (derived by Google Scholar). Models are in chronological order [99,102].
Geosciences 14 00255 g002
Table 1. A summary of the 12 selected models. * The columns “Country of application” and “Real case study” refer to the related cited paper in this review. “WWM” stands for “Wu and Waldron root reinforcement model”; “AFBM” stands for “Analytical Fiber Bundle Model” for root reinforcement.
Table 1. A summary of the 12 selected models. * The columns “Country of application” and “Real case study” refer to the related cited paper in this review. “WWM” stands for “Wu and Waldron root reinforcement model”; “AFBM” stands for “Analytical Fiber Bundle Model” for root reinforcement.
ModelCountry of Application *Real Case StudyHydrologyVegetation
SHALSTAB [75]United StatesTennessee Valley (colluvial soils), Mettman Ridge (colluvial soils), and Split Creek (silty sands)Derivation of the critical saturated soil height through TOPOG [76]-
dSLAM [77]United StatesCedar Creek basin, Oregon Coast Range—real landslides not well documentedKinematic wave groundwater model [78]Mechanical basal root reinforcement (WWM) as constant Cr
SHETRAN [79] United Kingdom of Great Britain and Northern Ireland Kirkton research catchment in Balquhidder, Scotland—real landslides not well documentedLateral flow: Boussinesq equation; vertical flow: 1D
Richards’ equation;
snowmelt
ET, interception, spatially variable land use
Iverson [82]--Reduced form of Richards’ equation-
TRIGRS [84]--Richards’ equation per [85]Extension by [86] comprises root cohesion and tree surcharge
SLIP [90]ItalyDifferent areas of ItalyIncrease in the saturated portion of the soil with rainfallExtension by [91] comprises interception and root cohesion
SUSHI [92]ItalyMay 1998 Sarno landslides on pyroclastic soils from Campania (Italy)1D Richards’ equation (Hydro-SUSHI module)ET as function of LAI, interception
tRIBS-VEGGIE-Landslide [93]Puerto RicoLuquillo forest (Puerto Rico)—no real landslides1D Richards’ equationHydrological effects: transpiration [95], evaporation reduction by canopies;
vegetation variable with time, space, and depth (root distribution and root water uptake);
root cohesion
MD-STAB [96]United StatesNovember 1996, Coos Bay, Mettnam ridge debris flow (gravelly sand)Hydrostatic condition Only mechanical reinforcement: basal and lateral, variable with depth
SOSlope [98]--Macropore water pressure and matrix suctionMechanical lateral root reinforcement (RBM)
Lizárraga et al. [99]ItalyMay 1998 landslides on pyroclastic soils from Campania (Italy)Authors used TRIGRS to compute hydrology-
Abdollahi et al. [102]United States2019 Las Lomas watershed (California) shallow landslide (sandy loam and loam soils)1D Richards’ equation ET
Mechanical basal root reinforcement (AFBM)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sannino, G.; Bordoni, M.; Bittelli, M.; Meisina, C.; Tomei, F.; Valentino, R. Deterministic Physically Based Distributed Models for Rainfall-Induced Shallow Landslides. Geosciences 2024, 14, 255. https://doi.org/10.3390/geosciences14100255

AMA Style

Sannino G, Bordoni M, Bittelli M, Meisina C, Tomei F, Valentino R. Deterministic Physically Based Distributed Models for Rainfall-Induced Shallow Landslides. Geosciences. 2024; 14(10):255. https://doi.org/10.3390/geosciences14100255

Chicago/Turabian Style

Sannino, Giada, Massimiliano Bordoni, Marco Bittelli, Claudia Meisina, Fausto Tomei, and Roberto Valentino. 2024. "Deterministic Physically Based Distributed Models for Rainfall-Induced Shallow Landslides" Geosciences 14, no. 10: 255. https://doi.org/10.3390/geosciences14100255

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop