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

A comparison of machine learning models for the mapping of groundwater spring potential

Groundwater resources are vitally important in arid and semi-arid areas meaning that spatial planning tools are required for their exploration and mapping. Accordingly, this research compared the predictive powers of five machine learning models for groundwater potential spatial mapping in Wadi az-Zarqa watershed in Jordan. The five models were random forest (RF), boosted regression tree (BRT), support vector machine (SVM), mixture discriminant analysis (MDA), and multivariate adaptive regression spline (MARS). These algorithms explored spatial distributions of 12 hydrological-geological-physiographical (HGP) conditioning factors (slope, altitude, profile curvature, plan curvature, slope aspect, slope length (SL), lithology, soil texture, average annual rainfall, topographic wetness index (TWI), distance to drainage network, and distance to faults) that determine where groundwater springs are located. The area under the curve (AUC) of the receiver operating characteristic (ROC) curve was employed to evaluate the prediction accuracies of the five individual models. Here the results were ranked in descending order as MDA (83.2%), RF (80.6%), SVM (80.2%), BRT (78.0%), and MARS (75.5%).The results show good potential for further use of machine learning techniques for mapping groundwater spring potential in other places where the use and management of groundwater resources is essential for sustaining rural or urban life

Effects of hydrological events on morphological evolution of a fluvial system.

This study quantifies morphological evolution of the Dez River, Iran, from 1955 to 2016. The approach uses a sequence of Landsat images, aerial photos, and topographic maps. In addition, the hydrological data including average daily discharge and yearly maximum discharge at the Dezful hydrological station for the period (1955-2016) were used. The study reach was divided into 48 meander loops from upstream to downstream. Active channel width (w) was determined at 10 m intervals and changes assessed along the study reach of the Dez River. Morphological indices including sinuosity index; straight meander length; centerline flow length; erosion area; erodible length channel migration; centerline elongation; and radius of curvature were calculated in the reach. Results showed that the study reach of the Dez River changed dramatically in response to major floods, although the general trend is towards a narrowing of active channel width by 38% in the period 1955-2016. Results show that most of the meander loops in the study area extended and expanded. Between 1989 and 1995, all types of meander change were observed. There was a direct correlation between the frequency of hydrological events (flood days) bigger than 2-years return period and elongation of bends

Assessment of fractal dimension and geometrical characteristics of the landslides identified in North of Tehran, Iran

The aim of the presented study is to assess the fractal dimension (D) and the geometrical characteristics (length and width) of the landslides identified in North of Tehran, Iran. At first, the landslide locations (528 landslides) were identified by interpretation of aerial photographs, satellite images and field surveys, and then to calculate the fractal dimension (D), we used the computer programming named as FRACEK. In the next step, geometrical characteristics of each landslide such as length (L) and width (W) were calculated by ArcGIS software. The landslide polygons were digitized from the mentioned landslide inventory map and rotated based on movement direction. The fractal dimension for all landslides varied between 1.665 and 1.968. Subsequently, the relationship between the length/width ratios and theirs fractal D values for 528 landslides was calculated. The results showed that correlation coefficients (R), which are different regression models such as exponential, linear, logarithmic, polynomial, and power, between D and L/W ratio are relatively high, respectively (0.75, 0.75, 0.76, 0.78, and 0.75). It can be concluded that the fractal dimension values and geometry characteristics of landslides would be useful indices for the management of hazardous areas, susceptible slopes, land use planning, and landslide hazard mitigation

Crossing the Phantom Divide Line in a DGP-Inspired F(R,ϕ)F(R,\phi)-Gravity

We study possible crossing of the phantom divide line in a DGP-inspired $F(R,\phi)$ braneworld scenario where scalar field and curvature quintessence are treated in a unified framework. With some specific form of $F(R,\phi)$ and by adopting a suitable ansatz, we show that there are appropriate regions of the parameters space which account for late-time acceleration and admit crossing of the phantom divide line.Comment: 23 Pages, 10 figs, Submitted to JCA

Regional prediction of landslide hazard using probability analysis of intense rainfall in the Hoa Binh province, Vietnam.

The main objective of this study is to assess regional landslide hazards in the Hoa Binh province of Vietnam. A landslide inventory map was constructed from various sources with data mainly for a period of 21 years from 1990 to 2010. The historic inventory of these failures shows that rainfall is the main triggering factor in this region. The probability of the occurrence of episodes of rainfall and the rainfall threshold were deduced from records of rainfall for the aforementioned period. The rainfall threshold model was generated based on daily and cumulative values of antecedent rainfall of the landslide events. The result shows that 15-day antecedent rainfall gives the best fit for the existing landslides in the inventory. The rainfall threshold model was validated using the rainfall and landslide events that occurred in 2010 that were not considered in building the threshold model. The result was used for estimating temporal probability of a landslide to occur using a Poisson probability model. Prior to this work, five landslide susceptibility maps were constructed for the study area using support vector machines, logistic regression, evidential belief functions, Bayesian-regularized neural networks, and neuro-fuzzy models. These susceptibility maps provide information on the spatial prediction probability of landslide occurrence in the area. Finally, landslide hazard maps were generated by integrating the spatial and the temporal probability of landslide. A total of 15 specific landslide hazard maps were generated considering three time periods of 1, 3, and 5 years

Application of frequency ratio, statistical index, and weights-of-evidence models and their comparison in landslide susceptibility mapping in Central Nepal Himalaya

The Mugling–Narayanghat road section falls within the Lesser Himalaya and Siwalik zones of Central Nepal Himalaya and is highly deformed by the presence of numerous faults and folds. Over the years, this road section and its surrounding area have experienced repeated landslide activities. For that reason, landslide susceptibility zonation is essential for roadside slope disaster management and for planning further development activities. The main goal of this study was to investigate the application of the frequency ratio (FR), statistical index (SI), and weights-of-evidence (WoE) approaches for landslide susceptibility mapping of this road section and its surrounding area. For this purpose, the input layers of the landslide conditioning factors were prepared in the first stage. A landslide inventory map was prepared using earlier reports, aerial photographs interpretation, and multiple field surveys. A total of 438 landslide locations were detected. Out these, 295 (67 %) landslides were randomly selected as training data for the modeling using FR, SI, and WoE models and the remaining 143 (33 %) were used for the validation purposes. The landslide conditioning factors considered for the study area are slope gradient, slope aspect, plan curvature, altitude, stream power index, topographic wetness index, lithology, land use, distance from faults, distance from rivers, and distance from highway. The results were validated using area under the curve (AUC) analysis. From the analysis, it is seen that the FR model with a success rate of 76.8 % and predictive accuracy of 75.4 % performs better than WoE (success rate, 75.6 %; predictive accuracy, 74.9 %) and SI (success rate, 75.5 %; predictive accuracy, 74.6 %) models. Overall, all the models showed almost similar results. The resultant susceptibility maps can be useful for general land use planning

Landslide susceptibility mapping at VAZ watershed (Iran) using an artificial neural network model: a comparison between multilayer perceptron (MLP) and radial basic function (RBF) algorithms

Landslide susceptibility and hazard assessments are the most important steps in landslide risk mapping. The main objective of this study was to investigate and compare the results of two artificial neural network (ANN) algorithms, i.e., multilayer perceptron (MLP) and radial basic function (RBF) for spatial prediction of landslide susceptibility in Vaz Watershed, Iran. At first, landslide locations were identified by aerial photographs and field surveys, and a total of 136 landside locations were constructed from various sources. Then the landslide inventory map was randomly split into a training dataset 70 % (95 landslide locations) for training the ANN model and the remaining 30 % (41 landslides locations) was used for validation purpose. Nine landslide conditioning factors such as slope, slope aspect, altitude, land use, lithology, distance from rivers, distance from roads, distance from faults, and rainfall were constructed in geographical information system. In this study, both MLP and RBF algorithms were used in artificial neural network model. The results showed that MLP with Broyden–Fletcher–Goldfarb–Shanno learning algorithm is more efficient than RBF in landslide susceptibility mapping for the study area. Finally the landslide susceptibility maps were validated using the validation data (i.e., 30 % landslide location data that was not used during the model construction) using area under the curve (AUC) method. The success rate curve showed that the area under the curve for RBF and MLP was 0.9085 (90.85 %) and 0.9193 (91.93 %) accuracy, respectively. Similarly, the validation result showed that the area under the curve for MLP and RBF models were 0.881 (88.1 %) and 0.8724 (87.24 %), respectively. The results of this study showed that landslide susceptibility mapping in the Vaz Watershed of Iran using the ANN approach is viable and can be used for land use planning

Control of style-of-faulting on spatial pattern of earthquake-triggered landslides

Predictive mapping of susceptibility to earthquake-triggered landslides (ETLs) commonly uses distance to fault as spatial predictor, regardless of style-of-faulting. Here, we examined the hypothesis that the spatial pattern of ETLs is influenced by style-of-faulting based on distance distribution analysis and Fry analysis. The Yingxiu–Beichuan fault (YBF) in China and a huge number of landslides that ruptured and occurred, respectively, during the 2008 Wenchuan earthquake permitted this study because the style-of-faulting along the YBF varied from its southern to northern parts (i.e. mainly thrust-slip in the southern part, oblique-slip in the central part and mainly strike-slip in the northern part). On the YBF hanging-wall, ETLs at 4.4–4.7 and 10.3–11.5 km from the YBF are likely associated with strike- and thrust-slips, respectively. On the southern and central parts of the hanging-wall, ETLs at 7.5–8 km from the YBF are likely associated with oblique-slips. These findings indicate that the spatial pattern of ETLs is influenced by style-of-faulting. Based on knowledge about the style-of-faulting and by using evidential belief functions to create a predictor map based on proximity to faults, we obtained higher landslide prediction accuracy than by using unclassified faults. When distance from unclassified parts of the YBF is used as predictor, the prediction accuracy is 80%; when distance from parts of the YBF, classified according to style-of-faulting, is used as predictor, the prediction accuracy is 93%. Therefore, mapping and classification of faults and proper spatial representation of fault control on occurrence of ETLs are important in predictive mapping of susceptibility to ETLs