Predicting gully occurrence at watershed scale: Comparing topographic indices and multivariate statistical models

In this study, the ability of five topographic indices to predict the gully trajectories observed in two adjacent watersheds located in Sicily (Italy) was evaluated. Two of these indices, named MSPI and MTWI, as far as we know, have never been employed to this aim. They were obtained by multiplying the stream power index (SPI) and the topographic wetness index (TWI), respectively, by the convergence index (CI). The predictive ability of the topographic indices was measured by using both cut-off independent (AUC: area under the receiver operating characteristic curve) and dependent statistics (Cohen's kappa index κ, sensitivity, specificity). These statistics were calculated also for 100 MARS (multivariate adaptive regression splines) and 100 LR (logistic regression) model runs, which used as predictors the topographic variables (i.e. contributing area, slope steepness, plan curvature and convergence index) combined into the five indices. Performance statistics of both topographic indices and statistical models were calculated using 100 random samples of 2 m grid cells, which were extracted only from flow concentration lines. This was done in order to focus the validation process on where gully erosion is more likely to occur. MSPI achieved the best predictive skill (AUC > 0.93; κ > 0.71) among the topographic indices and exhibited similar and better accuracy than local (i.e. trained and validated in the same watershed) and transferred (i.e. trained in one watershed and tested in the other one) LR models, respectively. On the other hand, MSPI performed similarly to transferred MARS runs (AUC > 0.92; κ > 0.71) but slightly worse than local MARS runs (AUC > 0.95; κ > 0.77). Based on the results of this experiment, it can be inferred that (i) including CI helps in detecting hollow areas where gullies are more likely to occur and (ii) MPSI can be a valid alternative to a data driven approach for mapping gully erosion susceptibility in areas where a gully inventory is not available, which is necessary to calibrate statistical models

Forward logistic regression for earth-flow landslide susceptibility assessment in the Platani river basin (southern Sicily, Italy)

Forward logistic regression has allowed us to derive an earth-flow susceptibility model for the Tumarrano river basin, which was defined by modeling the statistical relationships between an archive of 760 events and a set of 20 predictors. For each landslide in the inventory, a landslide identification point (LIP) was automatically produced as corresponding to the highest point along the boundary of the landslide polygons, and unstable conditions were assigned to cells at a distance up to 8m. An equal number of stable cells (out of landslides) was then randomly extracted and appended to the LIPs to prepare the dataset for logistic regression. A model building strategy was applied to enlarge the area included in training the model and to verify the sensitivity of the regressed models with respect to the locations of the selected stable cells. A suite of 16 models was prepared by randomly extracting different unoverlapping stable cell subsets that have been appended to the unstable ones. Models were finally submitted to forward logistic regression and validated. The results showed satisfying and stable error rates (0.236 on average, with a standard deviation of 0.007) and areas under the receiver operating characteristic (ROC) curve (AUCs) (0.839 for training and 0.817 for test datasets) as well as factor selections (ranks and coefficients). As regards the predictors, steepness and large-profile and local-plan topographic curvatures were systematically selected. Clayey outcropping lithology, midslope drainage, local and midslope ridges, and canyon landforms were also very frequently (from eight to 15 times) included in the models by the forward selection procedures. The model-building strategy allowed us to produce a performing earth-flow susceptibility model, whose model fitting, prediction skill, and robustness were estimated on the basis of validation procedures, demonstrating the independence of the regressed model on the specific selection of the stable cells

Gully erosion susceptibility mapping using multivariate adaptive regression splines-replications and sample size scenarios

Soil erosion is a serious problem affecting numerous countries, especially, gully erosion. In the current research, GIS techniques and MARS (Multivariate Adaptive Regression Splines) algorithm were considered to evaluate gully erosion susceptibility mapping among others. The study was conducted in a specific section of the Gorganroud Watershed in Golestan Province (Northern Iran), covering 2142.64 km2 which is intensely influenced by gully erosion. First, Google Earth images, field surveys, and national reports were used to provide a gully-hedcut evaluation map consisting of 307 gully-hedcut points. Eighteen gully erosion conditioning factors including significant geoenvironmental and morphometric variables were selected as predictors. To model sensitivity of gully erosion, Multivariate Adaptive Regression Splines (MARS) was used while the Area Under the Receiver Operating Characteristic (ROC) Curve (AUC), drawing ROC curves, efficiency percent, Yuden index, and kappa were used to evaluate model efficiency. We used two different scenarios of the combination of the number of replications, and sample size, including 90%/10% and 80%/20% with 10 replications, and 70%/30% with 5, 10, and 15 replications for preparing gully erosion susceptibility mapping (GESM). Each one involves a various subset of both positive (presence), and negative (absence) cases. Absences were extracted as randomly distributed individual cells. Therefore, the predictive competency of the gully erosion susceptibility model and the robustness of the procedure were evaluated through these datasets. Results did not show considerable variation in the accuracy of the model, with altering the percentage of calibration to validation samples and number of model replications. Given the accuracy, the MARS algorithm performed excellently in predictive performance. The combination of 80%/20% using all statistical measures including SST (0.88), SPF (0.83), E (0.79), Kappa (0.58), Robustness (0.01), and AUC (0.84) had the highest performance compared to the other combinations. Consequently, it was found that the performance of MARS for modelling gully erosion susceptibility is quite consistent while changes in the testing and validation specimens are executed. The intense acceptable prediction capability of the MARS model verifies the reliability of the method employed for use of this model elsewhere and gully erosion studies since they are qualified to quickly generating precise and exact GESMs (gully erosion sensitivity maps) to make decisions and management edaphic and hydrologic features

Evaluation of multi-hazard map produced using MaxEnt machine learning technique

Natural hazards are diverse and uneven in time and space, therefore, understanding its complexity is key to save human lives and conserve natural ecosystems. Reducing the outputs obtained after each modelling analysis is key to present the results for stakeholders, land managers and policymakers. So, the main goal of this survey was to present a method to synthesize three natural hazards in one multi-hazard map and its evaluation for hazard management and land use planning. To test this methodology, we took as study area the Gorganrood Watershed, located in the Golestan Province (Iran). First, an inventory map of three different types of hazards including flood, landslides, and gullies was prepared using field surveys and different official reports. To generate the susceptibility maps, a total of 17 geo-environmental factors were selected as predictors using the MaxEnt (Maximum Entropy) machine learning technique. The accuracy of the predictive models was evaluated by drawing receiver operating characteristic-ROC curves and calculating the area under the ROC curve-AUCROC. The MaxEnt model not only implemented superbly in the degree of fitting, but also obtained significant results in predictive performance. Variables importance of the three studied types of hazards showed that river density, distance from streams, and elevation were the most important factors for flood, respectively. Lithological units, elevation, and annual mean rainfall were relevant for detecting landslides. On the other hand, annual mean rainfall, elevation, and lithological units were used for gully erosion mapping in this study area. Finally, by combining the flood, landslides, and gully erosion susceptibility maps, an integrated multi-hazard map was created. The results demonstrated that 60% of the area is subjected to hazards, reaching a proportion of landslides up to 21.2% in the whole territory. We conclude that using this type of multi-hazard map may be a useful tool for local administrators to identify areas susceptible to hazards at large scales as we demonstrated in this research

Predicting Earthquake-Induced Landslides by Using a Stochastic Modeling Approach: A Case Study of the 2001 El Salvador Coseismic Landslides

In January and February 2001, El Salvador was hit by two strong earthquakes that triggered thousands of landslides, causing 1259 fatalities and extensive damage. The analysis of aerial and SPOT-4 satellite images allowed us to map 6491 coseismic landslides, mainly debris slides and flows that occurred in volcanic epiclastites and pyroclastites. Four different multivariate adaptive regression splines (MARS) models were produced using different predictors and landslide inventories which contain slope failures triggered by an extreme rainfall event in 2009 and those induced by the earthquakes of 2001. In a predictive analysis, three validation scenarios were employed: the first and the second included 25% and 95% of the landslides, respectively, while the third was based on a k-fold spatial cross-validation. The results of our analysis revealed that: (i) the MARS algorithm provides reliable predictions of coseismic landslides; (ii) a better ability to predict coseismic slope failures was observed when including susceptibility to rainfall-triggered landslides as an independent variable; (iii) the best accuracy is achieved by models trained with both preparatory and trigger variables; (iv) an incomplete inventory of coseismic slope failures built just after the earthquake event can be used to identify potential locations of yet unreported landslides

Geospatial analysis of drought tendencies in the carpathians as reflected in a 50-year time series

Climate change is one of the most important issues of anthropogenic activities. The increasing drought conditions can cause water shortage and heat waves and can influence the agricultural production or the water supply of cities. The Carpathian region is also affected by this phenomenon; thus, we aimed at identifying the tendencies between 1960 and 2010 applying the CarpatClim (CC) database. We calculated the trends for each grid point of CC, plotted the results on maps, and applied statistical analysis on annual and seasonal level. We revealed that monthly average temperature, maximum temperature and evapotranspiration had similar patterns and had positive trends in all seasons except autumn. Precipitation also had a positive trend, but it had negative values in winter. The geospatial analysis disclosed an increasing trend from West to East and from north to west. A simple binary approach (value of 1 above the upper quartile in case of temperature and evapotranspiration, value of 1 below the lower quartile; 0 for the rest of the data) helped to identify the most sensitive areas where all the involved climatic variables exceeded the threshold: Western Hungary and Eastern Croatia. Results can help to prepare possible mitigation strategies to climate change and both landowners and planners can draw the conclusions

A multi-scale regional landslide susceptibility assessment approach: the SUFRA_SICILIA (SUscettibilit\ue0 da FRAna in Sicilia) project

The SUFRA project is based on a three level susceptibility mapping. According to the availability of more detailed data, the three scale for susceptibility mapping are increased respect to the ones suggested by the TIER group to 1:100,000, 1:50,000 and 1:25,000/1:10,000. The mapping levels exploit climatic, soil use (CORINE2009) and seismic informative layers, differentiating in the details of the core data (geology and topography), in the quality and resolution of the landslide inventory and in the modelling approach (Tab. 1). SUFRA_100 is based on a heuristic approach which is applied by processing a geologic layer (produced by ARTA integrating pre-CARG 1:100,000 geologic maps); the DEM exploited are IGMI 250m and the mapping units are 1km side square cells. Models are validated with respect to the PAI LIPs (Landslide Identification Points) which are reclassified adopting a simplified scheme. Output cuts of SUFRA100 will be referred to administrative boundaries (provinces). SUFRA50 is based on statistical analysis of new CARG geologic maps and 20m (ITA2000) - 2m (ATA2007) DEM. The mapping units are 500m and 50m cells, hydrographic and hydro-morphometric units. The landslide inventory is the IFFI2012_LIPs (first level) which is the result of the conversion in IFFI format of the PAI archive, which will be supported by remote landslide mapping (exploiting the ATA2007 aerial photos), according to the IFFI first level approach. Validation of the models will be performed exploiting both random spatial partition and temporal partition methods. Output cuts of SUFRA50 will be based on physiographic (basin) and administrative (municipalities) boundaries. SUFRA10/25 is based on statistical analysis of new CARG geologic maps (remotely and field adapted) and 2m (ATA2007) DEM. The mapping units are the slope units (SLUs) which are derived by further partitioning the hydro-morphometric units so to obtain closed morphodynamic units. The landslide inventories is the IFFI2012 which is the results of a field supported (on focus) landslide remote systematic mapping, according to the IFFI full level approach. Examples of SUFRA_100, SUFRA_50 and SUFRA_10 are presented for some representative key sector of Sicily. First results attest for the feasibility and goodness of the proposed protocol. The SUFRA program aims at enabling the regional governmental administration to cope with landslide prevision, which is the required operational concept in land management and planning. PAI has been a great advance with respect to the \u201cpre-SARNO\u201d conditions, but it is very exposed to fail: it is a blind approach for new activations; it is critically dependent on the quality of the landslide inventories; it cannot project the susceptibility outside the landslide area

Identification of debris flow initiation zones using topographic model and airborne laser scanning data

Empirical multivariate predictive models represent an important tool to estimate debris flow initiation areas. Most of the approaches used in modelling debris flows propagation and deposit phases required identifying release (starting point) area or source area. Initiation areas offer a good overview to point out where field investigation should be conducted to establish a detailed hazard map. These zones, usually, are arbitrarily chosen which affect the model outputs; hence, there is a need to have accurate and automated means of identifying the release area. In addition to this, the resolution of the terrain dataset also affects the results of the detection of source areas. In this study, airborne laser scanning (ALS) data was used because of its robustness in providing detailed terrain attributes at high resolution. Primary and secondary conditioning parameters were derived from digital elevation model (DEM) as input into the modelling process. Three models were executed at different spatial resolution scales: 5, 10 and 15 m, respectively. MARSpline multivariate data mining predictive approach was implemented using morphometric indices and topographical derived parameter as independent variables. A statistics validation was calculated to estimate the optimal pixel size, 1200 randomly sample data were generated from existing inventory data. Debris flows and no-debris flows were categorized, and the transform to continuous integer (1 and 0), respectively. To achieve this, the data set was divided into two, 70% (840) for the training dataset and 30% (360) for validation. The best model was selected based on the model performance using the generalized cross validation (GCV) and the receiver operating characteristic (ROC) curve/area under curve (AUC) values. Conditioning parameters were numerically optimized to identify the arbitrarily maximum model basis function for eleven variables, using MARSplines analysis (algorithm). The three most influencing topographic parameters identified are topographic roughness index (TRI), slope angle, and specific catchment area (SCA) with the percentage values of participation in the model of 100, 93, and 86%, respectively. The chosen function appeared to describe the analysed correlation sufficiently well. Consequently, three stages of optimization were made to determine the optimized source areas is possible with 10 m pixel size, 200 maximum basis functions and 3 maximum interactions, resulting into 82% ROC train and 80% test, GCV 0.189 and 85% correlation coefficient. The result will be of great contribution to the advancement of a broad understanding of the dynamics of debris flows hazard and mitigations at regional level which; that is resourceful for comprehensive slope management for safe urban planning decision-making process and debris flow disaster management