Mapping Landslide Prediction through a GIS-Based Model: A Case Study in a Catchment in Southern Italy

Shallow landslides are an increasing concern in Italy and worldwide because of the frequent association with vegetation management. As vegetation cover plays a fundamental role in slope stability, we developed a GIS-based model to evaluate the influence of plant roots on slope safety, and also included a landslide susceptibility map. The GIS-based model, 4SLIDE, is a physically based predictor for shallow landslides that combines geological, topographical, and hydrogeological data. The 4SLIDE combines the infinite slope model, TOPMODEL (for the estimation of the saturated water level), and a vegetation root strength model, which facilitates prediction of locations that are more susceptible for shallow landslides as a function of forest cover. The aim is to define the spatial distribution of Factor of Safety (FS) in steep-forested areas. The GIS-based model 4SLIDE was tested in a forest mountain watershed located in the Sila Greca (Cosenza, Calabria, South Italy) where almost 93% of the area is covered by forest. The sensitive ROC analysis (Receiver Operating Characteristic) indicates that the model has good predictive capability in identifying the areas sensitive to shallow landslides. The localization of areas at risk of landslides plays an important role in land management activities because landslides are among the most costly and dangerous hazards.


Introduction
Landslides are defined as mass movements of soil and rocks along a slope. Landslides occur when slopes undergo a decrease in the shear strength of the hillside material due to an increase in the shear stress, or due to a combination of natural ecosystem processes and anthropogenic activities [1,2]. In some cases, human activities can trigger or, more often, they can accelerate the dynamics of natural processes by modifying the threshold of occurrence of landslides, thus potentially increasing risks for people. A frequent determinant of these processes consists of extreme weather events in geomorphologically unstable areas. Furthermore, land use changes and different land management practices can interact with those factors due to their influence on soil hydrological response, thus

The 4SLIDE
The 4SLIDE GIS-based model was set up as an add-on toolbox component to ESRI's ArcGIS© 10.3 for Desktop for creating a toolbox. The 4SLIDE tool uses the infinite slope stability equation [25] for calculating the Factor of Safety Index and it applies the hydrological model "TOPMODEL" (a TOPography based hydrological MODEL) [26] for evaluating the landslide risk at the catchment scale. In addition, this model considers soil reinforcement due to the presence of plants roots through a vegetation root strength model according to the theory presented in Wu et al. [27]. The main output of the tool is the Factor of Safety (FS) raster map called a "soil stability map". An overview of the workflow methodology of 4SLIDE in ArcGIS environment for landslide hazard prediction is shown in Figure 1.
Geosciences 2020, 10, x FOR PEER REVIEW 3 of 22 workflow methodology of 4SLIDE in ArcGIS environment for landslide hazard prediction is shown in Figure 1. The soil stability map provides, for each pixel, a value used to classify the terrain in risk classes subdivided in terms of the Factor of Safety Index (Table 1). In Figure S1, the GIS-based model development procedure was described.  The soil stability map provides, for each pixel, a value used to classify the terrain in risk classes subdivided in terms of the Factor of Safety Index (Table 1). In Figure A1, the GIS-based model development procedure was described.

Factor of Safety Index
The FS index is defined as the ratio of stresses resisting failure (Sr) to the stress required to bring the slope into a state of limit equilibrium (SL) along a failure surface according to the infinite slope stability equation (ISSE). This limit equilibrium analysis can be applied for all types of failure slip surfaces. The most useful analytical method for forested slopes is the planar infinite slope analysis [28]. This analysis assumes that the thickness of the soil mantle is small compared to the slope length and the failure plane is parallel to the slope [27,29,30]. For forested slopes with water tables less than or equal to soil depth (in meters), the FS index is [3]: where c R is the cohesion attributed to root systems (kPa), γ m is the unit weight of soil at field moisture content (kN m −3 ), γ sat is the unit weight of saturated soil (kN m −3 ), Z the vertical soil thickness (m), h is the vertical height of the water table (m), β is the slope angle (in radians), and W the vegetation surcharge (kPa). In areas with FS < 1.3, there is a probability of failure.

Hydrological Model
The TOPMODEL approach has become widely used for hydrological catchment modeling [26,31,32]. Different than other models, TOPMODEL does not divide the catchment into homogeneous units, but rather employs distribution functions to represent the catchment's natural heterogeneity and its effects on hydrological processes. Therefore, the distribution of wetness throughout a catchment can be simulated easily and with low computational demand. These developments follow field observations that show that higher soil moisture or areas of surface saturation tend to occur in convergent hollows. It has also been reported that landslides most commonly originate in areas of topographic convergence [17,26]. As per Beven and Kirkby [32], the index of hydrological similarity (I) is based on the topographic index of Kirkby and Statham [33] and defined as: where, a is the specific catchment area, defined as upslope area per unit contour length (m), and β is the slope angle of the ground surface ( Figure A6). This index describes the tendency of water to accumulate (a) and to be moved downslope by gravitational forces (β). TOPMODEL can be used for determination of the groundwater levels (z) ( Figure A7).; this theory assumes that transmissivity decreases exponentially with increasing depth to the groundwater table z (meter below surface), the hydraulic gradient is equal to the surface gradient (β), and lateral flow in the unsaturated zone can be neglected. A mean depth to the water table, z at each location i, is given by the equation: where, Ti (m 2 t −1 ) is the transmissivity if the groundwater level is at the ground surface, f (m −1 ) is a shape factor describing the exponential decrease of conductivity with soil depth, R is recharge, and A is the catchment area.

The Vegetation Root Strength Model
Plants represent an effective means of defense against surface erosion. The increase in soil shear strength due to the presence of roots was assessed according to the model presented by Wu et al. [27]. This model allows for simple and quick calculation of soil reinforcement by roots using tensile strength and root distribution information, resulting in additional cohesion [34]. If the soil is rooted, Geosciences 2020, 10, 309 5 of 22 the increased soil resistance to slides can be expressed as an additional cohesion c R (kPa) [35,36] that increases shear strength due to the presence of roots: where, T R is the tensile strength and a r is the Root Area Ratio (RAR), both specified per tree diameter class i, N is the number of classes considered, k is the factor accounting for decomposition of root tensile strength according to the bending angle of roots with respect to the shear plane, and k" is a reduction factor accounting for the non-simultaneous breaking of roots [20].

4SLIDE Tool Implementation
4SLIDE GIS-based model was implemented in ModelBuilder application. This is an easy-to-use application for creating, editing, and managing models in ArcGIS Desktop 10© (Esri). The data processing workflow for defining the FS index map is shown in Figure 1. In the first phase, the fundamental data are collected and inserted within the Model Builder. Subsequently, the internal tools of ArcGIS are used to determine the morphological and hydrological parameters needed to run the TOPMODEL equations. Lastly, the root cohesion parameters are inserted into Model Builder and, through the Raster Calculator tool (ArcGIS tool), the infinite slope stability equation is employed ( Figure A1). At this stage, the model computes the Factor of Safety index and a map is created.

Input and Output Tool Data
The 4SLIDE tool requires a digital terrain model (DTM) for calculating the topographic factors like slope angle in radians (rad) and the drainage network. The hydrological data required by the model are: (i) transmissivity (m 2 hr −1 ) (the rate at which groundwater flows horizontally through an aquifer or soil), (ii) the effective drainable porosity (the fraction of voids per unit volume), (iii) the discharge of watershed outlet (m 3 s −1 ) (the amount of fluid that crosses a section of the area A in the time unit), and (iv) basin area (m 2 ). Moreover, the model requires the following soil data: (a) soil cohesion (in kPa), (b) internal angle of friction (radians), (c) dry unit soil weight (kN m −3 ), (d) wet unit soil weight (kN m −3 ), and (e) vertical soil depth (m). The last input parameter is root cohesion (kPa) for plants. Here, a homogeneous value for the entire study area or a raster map with different values can be used ( Figure A2). The model uses either numerical values (N.V.) (across all cells) or a raster map (a value that varies across different cells) for each parameter. The values that can be changed by the user to generate different simulations in the dialog box are reported in Table 2.  The output of 4SLIDE tool is a FS index map in the form of a raster file. FS values < 1 indicate slope failure (i.e., a cell where a landslide is highly likely), while values < 1.3 indicate unstable conditions. For display purposes, in this study, the FS index was grouped into five classes according to Pack et al. [37] (Table 1).

Case Study
Landslides are very diffused in large areas of the Italian territory [38], in particular, in the Calabria region (Southern Italy). In fact, in this territory, a great number of sites are quite prone to landsliding, due to the combination of their geological, morphological, climatic, and land management features [39,40]. This region is also characterized by a complex orography and generally small river basins that respond rapidly to rainfall. Intense meteorological events combined with these attributes of the territory give rise to violent shallow landslide phenomena characterized by rapid dynamics. About 70% of the rainfall occurs from October to March, with average annual precipitation ranging from 1000 to 2000 mm y −1 [41]. The study area is a forest catchment (coordinate of the center of the catchment: WGS 84/UTM zone 33N-EPSG: 32633) called "Bacino del Bonis", located in the border area between the municipalities of Acri and Longobucco ( Figure 2) [42].
The basin was instrumented in 1986 with three mechanical rain gauges (with tipping buckets, data compiled in 20-min intervals) located at the basin outlet (975 m a.s.l.) and at representative sites within the north-eastern (Petrarella: 1258 m a.s.l.) and southwestern (Don Bruno: 1175 m a.s.l.) parts of the catchment. Runoff was measured at the catchment outlet using a Thomson weir (capable of measuring discharges up to 17 m 3 s −1 ) equipped with a mechanical stage recorder. The basin area is 1385 km 2 , most of which is forested with an incised stream network. The area is covered by Calabrian pine (Pinus laricio Poiret) with small areas of mixed stands of chestnut (Castanea sativa Mill.) and alder riparian forest (Alnus glutinosa (L.) Gaertn) ( Figure 3). small river basins that respond rapidly to rainfall. Intense meteorological events combined with these attributes of the territory give rise to violent shallow landslide phenomena characterized by rapid dynamics. About 70% of the rainfall occurs from October to March, with average annual precipitation ranging from 1000 to 2000 mm y −1 [41]. The study area is a forest catchment (coordinate of the center of the catchment: WGS 84/UTM zone 33N-EPSG: 32633) called "Bacino del Bonis", located in the border area between the municipalities of Acri and Longobucco ( Figure 2) [42].   The main stream has long steep slopes varying from 35 to 40% with the presence of several waterfalls and some isolated pools. The geology throughout the territory consists of four major complexes of different ages. The oldest complex is comprised of an acid igneous and metamorphic unit (Sila unit) of Paleozoic age and calcareous unit of Mesozoic age [43]. These rocks are covered by a Paludi formation of lower Miocene age [44], followed by a sedimentary succession of mio-pliocene [45] and; Quaternary sediments [46,47]. A complex of acidic intrusive plutonic rocks mainly characterizes the study basin. Soils are predominantly coarse-textured colluvium derived from volcanic rocks and classified as ultic haploxeralfs according to soil taxonomy [47]. Soil depth is shallow, with an average depth of about 1 m, and highly permeable to water. The water table rapidly responds during rainfall, rising and falling in response to rain intensity. The average annual precipitation is 1250 mm year −1 and average air temperature is 8.8 °C [47], with a typical mountain-Mediterranean climate regime, characterized by cool and wet winters and warm and dry summers. For this study, the soil properties considered are: soil cohesion, internal friction angle, unconfined compressive The main stream has long steep slopes varying from 35 to 40% with the presence of several waterfalls and some isolated pools. The geology throughout the territory consists of four major complexes of different ages. The oldest complex is comprised of an acid igneous and metamorphic unit (Sila unit) of Paleozoic age and calcareous unit of Mesozoic age [43]. These rocks are covered by a Paludi formation of lower Miocene age [44], followed by a sedimentary succession of mio-pliocene [45] and; Quaternary sediments [46,47]. A complex of acidic intrusive plutonic rocks mainly characterizes the study basin. Soils are predominantly coarse-textured colluvium derived from volcanic rocks and classified as ultic haploxeralfs according to soil taxonomy [47]. Soil depth is shallow, with an average depth of about 1 m, and highly permeable to water. The water table rapidly responds during rainfall, rising and falling in response to rain intensity. The average annual precipitation is 1250 mm year −1 and average air temperature is 8.8 • C [47], with a typical mountain-Mediterranean climate regime, characterized by cool and wet winters and warm and dry summers. For this study, the soil properties considered are: soil cohesion, internal friction angle, unconfined compressive strength, unit weight, and saturated hydraulic conductivity (Table 3). These soil data were determinate by 135 soil samples taken within the study area [48] ( Figure A3).

Calibration and Sensitivity Analysis
The 4SLIDE tool was applied in the study area and several input data for generating the FS index map were used. Firstly, high-resolution Digital Terrein Model (DTM) ( Figure A5) derived from a LiDAR survey ( Figure A4) conducted in the ALForLab project (www.alforlab.eu) were used. For vegetation, root cohesion data is described in Moresi et al. [36], while the geotechnical parameters were derived from Belloti & Selleri [49] ( Table 3). The soil map was made after a detailed geological survey inside Bonis basin (Figure 4). This map is a product of the ALFORLAB project created by a kriging geostatistics algorithm in QGIS software. These soil data were determinate by 135 soil samples taken within the study area [48] ( Figure  S3).

Calibration and Sensitivity Analysis
The 4SLIDE tool was applied in the study area and several input data for generating the FS index map were used. Firstly, high-resolution Digital Terrein Model (DTM) ( Figure S5) derived from a LiDAR survey ( Figure S4) conducted in the ALForLab project (www.alforlab.eu) were used. For vegetation, root cohesion data is described in Moresi et al. [36], while the geotechnical parameters were derived from Belloti & Selleri [49] ( Table 3). The soil map was made after a detailed geological survey inside Bonis basin (Figure 4). This map is a product of the ALFORLAB project created by a kriging geostatistics algorithm in QGIS software. Rainfall data were collected from the Cecita (CS) weather station, the nearest station to the study area (~10 km). A sensitivity analysis of the 4SLIDE parameters was also conducted. The 4SLIDE sensitivity was assessed by analyzing the changes with respect to a standard value resulting from the variation of the generic i-th independent variable (Xi) with respect to its basic value (XiR), while keeping all the other parameters (Xj) constant for j ≠ i. Therefore, to provide perspectives on the key drivers generating the FS map for the Bonis catchment, the sensitivity was assessed by changing the following parameters: Rainfall data were collected from the Cecita (CS) weather station, the nearest station to the study area (~10 km). A sensitivity analysis of the 4SLIDE parameters was also conducted. The 4SLIDE Geosciences 2020, 10, 309 9 of 22 sensitivity was assessed by analyzing the changes with respect to a standard value resulting from the variation of the generic i-th independent variable (X i ) with respect to its basic value (X iR ), while keeping all the other parameters (X j ) constant for j i. Therefore, to provide perspectives on the key drivers generating the FS map for the Bonis catchment, the sensitivity was assessed by changing the following parameters: 1. Cohesion 2. Friction angle 3. Soil unit weight 4. Root cohesion 5. Water Unit Weight By applying the "one-factor-at-a-time" method OAT [50,51], each parameter was changed by ±5% and ±10% from the standard value since these ranges of variation were consistent with observed variability in both the study area and in the literature [51,52]. Moreover, the changing of parameter values potentially reveals a linear (or non-linear) relationship among the parameter changes and model output, as well as the sensitivity of each parameter at different times during the simulation. This methodology shows the relative change of FS (∆FS) for the variation of each parameter with respect to the reference value (FS R ). FS R value was plotted against the relative variation of each independent parameter X i accounted for (∆X i ) with respect to the corresponding reference value (X iR ) to assess its slope mi as defined by: where, mi is the partial derivate of FS compared to its reference value (FS R ) with respect to the relative variation of the i-th parameter (X i ) compared to its reference value (X iR ). As such, the intrinsic linearity of Equation (1) was immediately checked. Thus, for higher values of mi, the sensitivity of FS towards the relative variation of X i is also higher.

4SLIDE GIS-Based Model Validation
The accuracy of the landslide prediction map produced by the 4SLIDE was evaluated by Relative Operating Characteristics (ROC) analysis [53]. The method provides a curve given by a confusion matrix of binary classification according to four possible outcomes: true positive (TP), true negative (TN), false positive (FP), and false negative (FN). The outcomes are derived by comparing results of the model with the ground truth survey (GTS), which are established by surveys. ROC Curve is a method that compares true positive rates against false positive rates. The Areas Under the Curve (A.U.C., which range from 0.5 to 1.0) are indices of accuracy assessment (Table 4); a nonparametric method for estimating A.U.C. was used for comparing two ROC curves [54]. Statistical analyses were conducted using R software program (www.r-project.org).  For each random point, a sample area of 20 m in diameter was investigated ( Figure 6). Each area was surveyed to check for the presence/absence of landslides using two methods. Firstly, GTS that was conducted in August 2017 for the entire basin. Secondly, a DTM was generated from LiDAR point clouds. "Hillshade" viewing allowed visualization of the terrain such that distinct landslide characteristics were identified within the sample areas. The database created included information on the location, dimensions, and other features of landslide occurrence and the locations of non-landslide areas to use these during the training and predictive phases. The database was used also for the ROC analysis. Fifty sample areas ( Figure 5) were defined and the two relief methods (GTS and LiDAR data) were superimposed to check for the presence/absence of landslides. Based on this overlap of both relief methods, each pixel was assigned 1 (presence) or 0 (absence) for landslide occurrence.

Performance of 4SLIDE (Calibration/Validation) using the Physical Similarity Method
The 4SLIDE considers soil features, vegetation effects, and hydraulic factors for determining slope stability. The resulting map is characterized by five stability classes ( Figure 6) according to the FS index value (Table 1).

Performance of 4SLIDE (Calibration/Validation) Using the Physical Similarity Method
The 4SLIDE considers soil features, vegetation effects, and hydraulic factors for determining slope stability. The resulting map is characterized by five stability classes ( Figure 6) according to the FS index value ( Table 1).
These results highlight that the areas with highest landslide susceptibility are mainly located on steep slopes and unconsolidated lithologies (e.g., colluvial material), which are generally adjacent to the river ( Figure 6). In contrast, areas upslope of the stream are generally characterized by low landslide susceptibility.

Performance of 4SLIDE (Calibration/Validation) using the Physical Similarity Method
The 4SLIDE considers soil features, vegetation effects, and hydraulic factors for determining slope stability. The resulting map is characterized by five stability classes ( Figure 6) according to the FS index value (Table 1).  The classification results show that 15% of the study area is classified in the "upper threshold of instability" (0.0 < FS < 0.5), 36% medium susceptibility (0.5 < FS < 1.3), and 20% in low to very low class of stability (FS > 1.3) (Figure 7). These results highlight that the areas with highest landslide susceptibility are mainly located on steep slopes and unconsolidated lithologies (e.g., colluvial material), which are generally adjacent to the river (Figure 6). In contrast, areas upslope of the stream are generally characterized by low landslide susceptibility.
According to A.U.C. method, the model performance between the GTS and DTM-LiDAR was 0.76 and 0.70, respectively (Figure 8). According to A.U.C. method, the model performance between the GTS and DTM-LiDAR was 0.76 and 0.70, respectively (Figure 8).
These results show that the model has good predictive capability for identifying the areas at risk of shallow landslides and the testing method using DTM-LIDAR data gave similar results compared to GTS. According to the obtained A.U.C. values and classification [54], the model predictive accuracy was assessed as good for both methods. Overall, the result showed good agreement between the susceptibility map and the landslide location data collected in the field. The 4SLIDE, like other predictive models, requires a large amount of information for initialization and parameterization. To further test the effect of the root systems, a simulation with 4SLIDE was performed, eliminating the effect of root cohesion. The main result of this simulation was a greater increase of unstable areas. The new FS index map was reported in Figure 6 and is characterized by five classes (Figure 8) according to the FS index value ( Table 1). The classification results show that 6% of the study area is classified in the "upper threshold of instability" (0.0 < FS < 0.5), 40% medium susceptibility (0.5 < FS < 1.3), and 28% in low to very low class of stability (FS > 1.3) (Figure 8). These results highlight that the areas with highest landslide susceptibility are mainly located on steep slopes. These results highlight that the areas with highest landslide susceptibility are mainly located on steep slopes and unconsolidated lithologies (e.g., colluvial material), which are generally adjacent to the river (Figure 6). In contrast, areas upslope of the stream are generally characterized by low landslide susceptibility.
According to A.U.C. method, the model performance between the GTS and DTM-LiDAR was 0.76 and 0.70, respectively (Figure 8). These results show that the model has good predictive capability for identifying the areas at risk of shallow landslides and the testing method using DTM-LIDAR data gave similar results compared to GTS. According to the obtained A.U.C. values and classification [54], the model predictive accuracy was assessed as good for both methods. Overall, the result showed good agreement between the susceptibility map and the landslide location data collected in the field. The 4SLIDE, like other We had tested the 4SLIDE tool's in the Bonis catchment without the benefit of vegetation. The simulation showed a significant increase in landslide risk areas, not only in the more sloping areas but also in the flat areas ( Figure 9). predictive models, requires a large amount of information for initialization and parameterization. To further test the effect of the root systems, a simulation with 4SLIDE was performed, eliminating the effect of root cohesion. The main result of this simulation was a greater increase of unstable areas.
The new FS index map was reported in Figure 6 and is characterized by five classes (Figure 8) according to the FS index value ( Table 1). The classification results show that 6% of the study area is classified in the "upper threshold of instability" (0.0 < FS < 0.5), 40% medium susceptibility (0.5 < FS < 1.3), and 28% in low to very low class of stability (FS > 1.3) (Figure 8). These results highlight that the areas with highest landslide susceptibility are mainly located on steep slopes. We had tested the 4SLIDE tool's in the Bonis catchment without the benefit of vegetation. The simulation showed a significant increase in landslide risk areas, not only in the more sloping areas but also in the flat areas ( Figure 9).

Sensitive Parameters and their Calibration to Simulate Future Hydrogeological Risk
Five cases were simulated and evaluated (Figure 10) forcing the model with five different initialization datasets of model sensitivity analysis.

Sensitive Parameters and Their Calibration to Simulate Future Hydrogeological Risk
Five cases were simulated and evaluated ( Figure 10) forcing the model with five different initialization datasets of model sensitivity analysis. FS was mainly sensitive to the friction angle, its slope being equal to 0.844 ± 0.005. Soil unit weight was the second most sensitive parameter (mi = 0.242 ± 0.012), followed by soil cohesion (mi = 0.186 ± 0.012) and finally root cohesion (mi = 0.063 ± 0.000) ( Table 5).  FS was mainly sensitive to the friction angle, its slope being equal to 0.844 ± 0.005. Soil unit weight was the second most sensitive parameter (m i = 0.242 ± 0.012), followed by soil cohesion (m i = 0.186 ± 0.012) and finally root cohesion (m i = 0.063 ± 0.000) ( Table 5).

Discussion
Landslides are among the most costly natural hazards, and, although large catastrophic events are rare, small mass movements are very frequent and produce a large amount of economic losses in the long term [1]. Therefore, defining measures to mitigate risks of landslide hazards is an urgent need, since it could have a great impact on current and future societies. In this context, predictive models can highlight areas which are predisposed to landslides and can be useful for preventive actions. In geomorphology, predictive modelling is used to provide relevant and functional information on earth surface processes and landforms throughout extensive areas because of the scarcity of more conventional data surveys [7,55]. Identification of areas at risk of landslide, determination of landslide locations, and assessment of susceptibility play important roles in regional management activities and mitigation measures. In this context, the 4SLIDE model develops landslide susceptibility maps in forest catchments using several parameters among which are root cohesion, soil characteristics, and terrain morphology. In fact, as reported in Formetta [56], the landslide susceptibility evaluation is formed by three main components: (1) hydrological model estimation of soil moisture and water table depth; (2) computing the factor of safety index (FS) based on the infinite slope model; and (3) GIS for visualization and calculation of the outputs. Landslide analyses at different scales can be conducted by different deterministic models [16]. One of the most common is "Distributed Shallow Landslide Analysis Models" (dSLAM), which employs physical numerical variables [16], but some other models (e.g., 4Slide) use the infinite slope stability theory, including TRIGRS, SINMAP, CHASM, and GEOtop-FS [57]. These models are useful for calculating the susceptibility of shallow landslides. The outputs of these models are the water table depth maps and landslide susceptibility maps. But these models have all different features and their own advantages and disadvantages; in fact, some of these are commercial software, and some are simply plug-ins for CAD or GIS software. Independent software like CHASM has an easy interface useful for the final user, while the model script-based ones (e.g., TRIGRS in Fortran) are less user-friendly, and these latter models are often open source licensed.
GEOtop is an independent model for investigating landslide prediction using DTMs and geotechnical and moisture parameters [58]. 4SLIDE, like SHALSTAB [59], analyzes landslide susceptibility at the catchment scale and uses steady-state hydrological processes and the infinite slope approach [60], while SINMAP [37], like 4SLIDE, uses numerical parameters and raster-based GIS images.
The 4SLIDE model differs from existing models because it allows for both the use of uniform numerical data throughout the catchment as well as spatially distributed data using raster maps. The latter method facilitates a much more detailed study of the prevention of shallow landslides. The main function of 4SLIDE is to calculate the FS index by analyzing the runoff process using slope values extracted from the DTM. This feature is useful for landslide study over time and over a large area, and 4SLIDE is an open-source model. In addition, the model uses root cohesion as an input parameter and, like the other input data, this parameter can also be spatially distributed allowing for inclusion of a mix of species. Furthermore, the model presents a simple-to-use graphical input interface. 4SLIDE was tested within the Bonis forest catchment where input data could be sampled and spatially represented on raster maps of soil and vegetation data. The model predicted unstable areas adjacent to the stream mainly due to steep terrain and less vegetation cover. Instead, the most stable areas were in zones with gentler slopes and more vegetation cover.
Accuracy assessment showed that the 4SLIDE model estimated the areas susceptible to landslides correctly, as confirmed by both ground surveys and DTM-LiDAR analysis. The performance of the model was evaluated and compared by ROC curves and A.U.C. values. In fact, A.U.C. values showed acceptable results, while the ground-truth sampling method produced slightly better results than the DTM-LiDAR method. The two methods are useful for determination of ground movement and deformation, characteristic of landslide development. In both cases, the performance of the landslide prediction depends on accurate spatial sampling of the ground surface. While both techniques are important, the GTS method can also detect inactive landslides where there are not surface changes; DTM-LIDAR surveys cannot detect this difference [61]. In addition, to test model uncertainty, a sensitivity analysis was performed by varying soil parameters in the ISSE. This analysis showed that a 5 or 10% decrease in the friction angle generates a high increase in slip areas. This is mainly due to the sandy soil that undergoes a greater effect on slope stability by varying the friction angle [7]. Other parameters of the model were much less sensitive to changes. For other soil types, sensitivity results may differ.
4SLIDE accurately determined areas that are susceptible to failure exclusively within the forest basin. The main advantage of the model is the ability to input spatially distributed data, including modifying vegetation cover within the basin and its effects on landslide susceptibility. Such assessments are useful for regional hazard and forest management planning, particularly in cases of removal of forest cover (e.g., fire, large storms, insect outbreaks, poor harvesting, and logging practices). For this reason, a simulation was performed, excluding the contribution of the all vegetation in the Bonis Catchment, which showed a high increase in the areas susceptible to landslides. The new simulation confirmed the importance of vegetation for slope stabilization and how poor forest management can lead to an increase of the susceptibility of landslides.
The main limitation of the model is that it cannot be applied in urban areas, where water drainage structures exist and where buildings and fills may influence slope stability [1]. Another limit of 4SLIDE is that it must interface within commercial GIS software; thus, although the model is freely available, users must have access to the software. Thus, the next steps will be to adapt the model for application in urban areas and utilize an open source software GIS like QGIS.

Conclusions
Natural hazards are those processes that occur naturally and when they damage or negatively affect people, property, or the environment, they result in disasters. Natural hazards and associated disasters are of increasing importance in today's globalized societies. In the 20th century, more than 4.5 million people died and 200 million were damaged by natural hazards. Landslides are one of the costliest risks and although large catastrophic events are rare, small movements are very frequent and produce a large amount of long-term economic losses. As such, it is necessary to develop innovative measures to mitigate risks because of the large impacts they can inflict on society.
Mapping inventories of these phenomena show the spatial attributes of landslides. In geomorphology, predictive modelling of landslides is used to provide relevant and functional information on earth surface process and landforms over extensive areas because these typologies are generally unavailable in conventional surveys [58]. In this study, a theoretical model for shallow landslides was developed and implemented. The 4SLIDE model is based on physical and local parameters, such as slope gradient, lithology, elevation, aspect, root cohesion, and land use/land cover. The model delineates the areas most prone to shallow landslides and is intended to for general planning and rather simple predictive applications. The model predicted spatial landslide occurrence with a moderate level accuracy in our catchment based on two validation methods. As such, the application of this model reduces the difficulties associated with data acquisition in harsh mountainous environments. Of course, better ground-based data will produce more reliable predictions. In terms of practical applications, we note that landslide susceptibility is necessary for a risk assessment and additional efforts should focus on the importance of both natural and urban vegetation and practices for reducing landslide hazards.