Search

Gully erosion is one of the advanced forms of water erosion. Identifying the effective factors and gully erosion predicting is one of the important tools to control and manage such phenomenon. The main purpose of this study is to evaluate the effect of four different resampling algorithms including cross-validation (5-fold and 10-fold) and bootstrapping (Bootstrap and Optimism bootstrap) on boosted regression tree (BRT), support vector machine (SVM), and random forest (RF) models in spatial modeling and evaluation of head-cut gully erosion in Konduran watershed. For this purpose, based on an extensive field survey, the points of the head-cut of the gully erosion were identified first, and a map of the distribution of head-cut gully erosion in the study area was prepared. Then 18 variable identify and prepare as factors affecting the occurrence of head-cut gully erosion. To assess the efficiency of the models, receiver operating characteristics (ROC) and area under the curve (AUC) were used. The results of the assessment indicated that the use of resampling algorithms increases the efficiency of the models. The integrated optimism-bootstrap-BRT, optimism-bootstrap-SVM, and Optimism-Bootstrap-RF models with AUC 0.85, 0.823 and 0.89 respectively, outperformed the cross-validation 5fold (BRT, SVM, RF), Cross-validation 10fold (BRT, SVM, RF) and Bootstrap (BRT, SVM, RF) integrated algorithms.

[EN] Floods are a consequence of extreme rainfall events. Although surface runoff generation is the origin of discharge, flood research usually focuses on lowlands where the impact is higher. Runoff and sediment delivery at slope and pedon scale receiving much less attention in the effort to understand flood behaviour in time and space. This is especially relevant in areas where, due to climatic and hydrogeological conditions, streams are ephemeral, so-called dry rivers (¿wadis¿, "ramblas" or ¿barrancos¿) that are widespread throughout the Mediterranean. This paper researches the relationship between water delivery at pedon and slope scale with dry river floods in Macizo del Caroig, Eastern Iberian Peninsula. Plots of 1x1, 1x2, 1x4, and 2x8 m located in the ¿El Teularet¿ Soil Erosion and Degradation Research Station were monitored from 2004 to 2014 to measure soil and water delivery. Rainfall and flow at the dry river Barranco de Benacancil were also monitored. Results show that runoff and sediment discharge were concentrated in few events during the 11 years of research. A single flood event was registered in the channel on September 28, 2009, however, the runoff was registered 160 times at the plots. Runoff discharge was dependent on the size of the plots, with larger plots yielding lower runoff discharge per unit area, suggesting short runoff-travel distance and duration. Three rainfall events contributed with 26% of the whole runoff discharge, and five achieved 56% of the runoff. We conclude that the runoff generated at the plot scale is disconnected from the main channel. From a spatial point of view, there is a decrease in runoff coefficient along the slope. From a temporal point of view, the runoff is concentrated in a few rainfall events. These results show that the runoff generated at plot and slope scale does not contribute to the floods except for rainfall events with more than 100 mm day-1. The disconnection of the runoff and sediment delivery is confirmed by the reduction in the runoff delivery at plot scale due to the control of the length of the plot (slope) on the runoff and sediment delivery., [ES] . Lasinundacionesson consecuencia de lluvias extremas. Aunque la generación de escorrentía superficial es el origen de la descarga, la investigación de inundaciones generalmente se enfoca en las tierras bajas donde el impacto es mayor. La escorrentía y la distribución de sedimentos a escala de pendiente y pedón reciben mucha menos atención en la comprensión del comportamiento de las inundaciones en el tiempo y el espacio. Esto es especialmente relevante en zonas donde, debido a las condiciones climáticas e hidrogeológicas, los cauces son efímeros. Son los llamados ríos secos (¿wadis¿, ¿ramblas¿ o ¿barrancos¿) muy extendidos por todo el Mediterráneo. Este artículo investiga la relación entre el suministro de agua a escala de pedón y ladera con las crecidas de ríos secos en Macizo del Caroig, este de la Península Ibérica. Las parcelas de 1x1, 1x2, 1x4 y 2x8 m localizadas en la Estación de Investigación de Erosión y Degradación de Suelos ¿El Teularet¿ fueron monitoreadas de 2004 a 2014 para medir la producción de suelo y agua. También se monitorearon las precipitaciones y el caudal en el río seco Barranco de Benacancil. Los resultados muestran que la escorrentía y la descarga de sedimentos se concentraron en pocos eventos durante los 11 años de investigación. Se registró un solo evento de inundación en el canal el 28 de septiembre de 2009, sin embargo, la escorrentía se registró 160 veces en las parcelas. La descarga de escorrentía dependió del tamaño de las parcelas. Las parcelas más grandes produjeron una menor descarga de escorrentía por unidad de área, lo que sugiere una corta distancia y duración del recorrido de escorrentía. Tres eventos de lluvia contribuyeron con el 26% de la descarga total de la escorrentía y cinco lograron el 56% de la escorrentía. Se concluye que la escorrentía generada a escala de la parcela está desconectada del canal principal. Desde un punto de vista espacial, hay una disminución en el coeficiente de escorrentía a lo largo de la pendiente. Desde un punto de vista temporal, la escorrentía se concentra en unos pocos eventos de lluvia. Estos resultados muestran que la escorrentía generada a escala de parcela y pendiente no contribuyen a las inundaciones excepto para eventos de lluvia con más de 100 mm día-1 . La desconexión de la escorrentía y la entrega de sedimentos se confirma por la reducción de la escorrentía a escala de parcela debido al control de la longitud (pendiente) sobre la escorrentía y la entrega de sedimentos., Artemi Cerda thanks the Co-operative Research program from the OECD (Biological Resource Management for Sustainable Agricultural Systems) for its support with the 2016 CRP fellowship (OCDE TAD/CRP JA00088807), POSTFIRE Project (CGL2013-47862-C2-1 and 2-R), and POSTFIRE_CARE Project (CGL2016-75178-C2-2-R) sponsored by the Spanish Ministry of Economy and Competitiveness and AEI/FEDER, UE. This paper was written as a result of the collaboration that was initiated due to the COST ActionES1306: Connecting European Connectivity research and COST CA18135 FIRElinks: Fire in the Earth System. Science and Society. We wish to thank the Department of Geography secretariat team (Nieves Gomez, Nieves Dominguez, and Susana Tomas) for their support for three decades to our research at the Soil Erosion and Degradation Research team (SEDER), with special thanks to the scientific researchers that as visitors from other research teams contributed to the SEDER research. And we also thank the Laboratory for Geomorphology technicians (Leon Navarro) for the key contribution to our research. The collaboration of the Geography and Environmental Sciences students was fruitful and enjoyable.

Gully erosion is the most active hydro-geomorphological phenomenon in the continental areas due to the high erosion rates triggered by the gully system. Monitoring and modelling gully development and gully distribution will contribute to understand landforms evolution and risk assessment. The purpose of the current research is to model head-cut gully erosion susceptibility (HCGES) using support vector machine (SVM), random forest (RF) and novel ensemble model of random forest with four Oblique methods (Logistic Regression, Ridge Regression, Partial least squares (PLS) and Support vector machine (SVM)) (hereafter called ensemble ORF) data mining models in Fareghan watershed, Hormozghan province, Iran. For this purpose, 14 variables influencing the gully development, were prepared and 145 head-cut gully erosion locations were identified in the study area. The efficiency of SVM, RF, ensemble ORF were evaluated based on receiver operating characteristic (ROC), the results have shown that all these three models are highly accurate and robust in predicting the head-cut gully erosion susceptibility zones. The results of the models were evaluated based on the area under the receiver operatic characteristic curve (AUC) in the validation stage presented that the efficiency of these models are 0.91, 0.94, and 0.96, respectively. Altitude and distance from the road in all three models were more important than other variables. The findings of this research will contribute to develop gully control strategies and to prevent the gully initiation where gully erosion is more susceptible.

Large damages and losses resulting from floods are widely reported across the globe. Thus, the identification of the flood-prone zones on a flood susceptibility map is very essential. To do so, 13 conditioning factors influencing the flood occurrence in Brisbane river catchment in Australia (i.e., topographic, water-related, geological, and land use factors) were acquired for further processing and modeling. In this study, artificial neural networks (ANN), deep learning neural networks (DLNN), and optimized DLNN using particle swarm optimization (PSO) were exploited to predict and estimate the susceptible areas to the future floods. The significance of the conditioning factors analysis for the region highlighted that altitude, distance from river, sediment transport index (STI), and slope played the most important roles, whereas stream power index (SPI) did not contribute to the hazardous situation. The performance of the models was evaluated against the statistical tests such as sensitivity, specificity, the area under curve (AUC), and true skill statistic (TSS). DLNN and PSO-DLNN models obtained the highest values of sensitivity (0.99) for the training stage to compare with ANN. Moreover, the validations of specificity and TSS for PSO-DLNN recorded the highest values of 0.98 and 0.90, respectively, compared with those obtained by ANN and DLNN. The best accuracies by AUC were evaluated in PSO-DLNN (0.99 in training and 0.98 in testing datasets), followed by DLNN and ANN. Therefore, the optimized PSO-DLNN proved its robustness to compare with other methods.

Complex interrelationships between landscape-level geoenvironmental factors and natural phenomena have rendered land degradation control measures ineffective. For control to be effective, this study argues that the interactions between different geoenvironmental factors and gully erosion (as an indicator of land degradation) should be more fully investigated and spatially mapped. To do so, gully locations of the Konduran watershed, Iran, were detected in the field and modeled in response to seventeen geoenvironmental factors using three machine learning methods, i.e., multivariate adaptive regression splines (MARS), random forest (RF), regularized random forest (RRF), and Bayesian generalized linear model (Bayesian GLM). The models' performance was validated, the relationship of gully occurrence with each factor was quantified, the probability of gully erosion (i.e., land degradation) was retrospectively estimated, and the spatially explicit maps of land degradation susceptibility were produced. Based on the area under the receiver operating characteristic curve (AUC), the RRF and MARS models with AUC = 0.98 achieved the greatest goodness-of-fit with the training dataset, whereas the RF model with AUC = 0.83 showed the greatest ability in predicting future gully occurrences. Further scrutinization using the sensitivity and specificity metrics demonstrated the efficiency of the RF model for correctly classifying the gully (sensitivity-training = 92%; sensitivity-validation = 90%) and non-gully (specificity-training = 95%; specificity-validation = 68%) pixels. Nearly 13% of the study area ended up being the hardest hit region due to their general characteristics of distance from roads and rives, altitude, and normalized difference vegetation index (NDVI) that were identified as the most influential factors in gully erosion occurrence. Given the resolution quality and reliable predictive accuracy, our spatially explicit maps of land susceptibility to gully erosion can be used by authorities and urban planners for identifying the target areas for rehabilitation and making more informed decisions for infrastructure development. Although our study was strictly focused on a certain region, our recommendations and implications are of global significance.

The predicts current and future flood risk in the Kalvan watershed of northwestern Markazi Province, Iran. To do this, 512 flood and non-flood locations were identified and mapped. Twenty flood-risk factors were selected to model flood risk using several machine learning techniques: conditional inference random forest (CIRF), the gradient boosting model (GBM), extreme gradient boosting (XGB) and their ensembles. To investigate the future (year 2050) effects of changing climates and changing land use on future flood risk, a general circulation model (GCM) with representative concentration pathways (RCPs) of the 2.6 and 8.5 scenarios by 2050 was tested for impacts on 8 precipitation variables. In addition, future land uses in 2050 was prepared using a CA-Markov model. The performances of the flood risk models were validated with Receiver Operating Characteristic-Area Under Curve (ROC-AUC) and other statistical analyses. The AUC value of the ROC curve indicates that the ensemble model had the highest predictive power (AUC = 0.83) and was followed by GBM (AUC = 0.80), XGB (AUC = 0.79), and CIRF (AUC = 0.78). The results of climate and land use changes on future flood-prone areas showed that the areas classified as having moderate to very high flood risk will increase by 2050. Due to the changes occurring with land uses and in climates, the area classified as moderate to very high risk increased in the predictions from all four models. The areal proportion classes of the risk zones in 2050 under the RCP 2.6 scenario using the ensemble model have changed of the following proportions from the current distribution Very Low = −12.04 %, Low = −8.56 %, Moderate = +1.56 %, High = +11.55 %, and Very High = +7.49 %. The RCP 8.5 scenario has caused the following changes from the present percentages: Very Low = −14.48 %, Low = −6.35 %, Moderate = +4.54 %, High = +10.61 %, and Very High = +5.67 %. The results of current and future flood risk mapping can aid planners and flood hazard managers in their efforts to mitigate impacts.

Identifying areas prone to flooding is a key step in flood hazard management. The purpose of this study is to develop and present a novel flood susceptibility model based on Bayesian Additive Regression Tree (BART) methodology. The predictive performance of new model is assessed via comparison with the Naïve Bayes (NB) and Random Forest (RF) based methods that were previously published in the literature. All models were tested on a real case study based in the Kan watershed in Iran. The following fifteen climatic and geo-environmental variables were used as inputs into all flood susceptibility models: altitude, aspect, slope, plan curvature, profile curvature, drainage density, distance from river distance from road, stream power index (SPI), topographic wetness index (TPI), topographic position index (TPI), curve number (CN), land use, lithology and rainfall. Based on the existing flood field survey and other information available for the analyzed area, a total of 118 flood locations were identified as potentially prone to flooding. The data available were divided into two groups with 70% used for training and 30% for validation of all models. The receiver operating characteristic (ROC) curve parameters were used to evaluate the predictive accuracy of the new and existing models. Based on the area under curve (AUC) the new BART (86%) model outperformed the NB (80%) and RF (85%) models. Regarding the importance of input variables, the results obtained showed that the altitude and distance from the river are the most important variables for assessing flooding susceptibility.

This study predicts forest fire susceptibility in Chaloos Rood watershed in Iran using three machine learning (ML) models—multivariate adaptive regression splines (MARS), support vector machine (SVM), and boosted regression tree (BRT). The study utilizes 14 set of fire predictors derived from vegetation indices, climatic variables, environmental factors, and topographical features. To assess the suitability of the models and estimating the variance and bias of estimation, the training dataset obtained from the Natural Resources Directorate of Mazandaran province was subjected to resampling using cross validation (CV), bootstrap, and optimism bootstrap techniques. Using variance inflation factor (VIF), weight indicating the strength of the spatial relationship of the predictors to fire occurrence was assigned to each contributing variable. Subsequently, the models were trained and validated using the receiver operating characteristics (ROC) area under the curve (AUC) curve. Results of the model validation based on the resampling techniques (non, 5- and 10-fold CV, bootstrap and optimism bootstrap) produced AUC values of 0.78, 0.88, 0.90, 0.86 and 0.83 for the MARS model; 0.82, 0.82, 0.89, 0.87, 0.84 for the SVM and 0.87, 0.90, 0.90, 0.90, 0.91 for the BRT model. Across the individual model, the 10-fold CV performed best in MARS and SVM with AUC values of 0.90 and 0.89. Overall, the BRT outperformed the other models in all ramification with highest AUC value of 0.91 using optimism bootstrap resampling algorithm. Generally, the resampling process enhanced the prediction performance of all the models.

Flash flooding is considered one of the most dynamic natural disasters for which measures need to be taken to minimize economic damages, adverse effects, and consequences by mapping flood susceptibility. Identifying areas prone to flash flooding is a crucial step in flash flood hazard management. In the present study, the Kalvan watershed in Markazi Province, Iran, was chosen to evaluate the flash flood susceptibility modeling. Thus, to detect flash flood-prone zones in this study area, five machine learning (ML) algorithms were tested. These included boosted regression tree (BRT), random forest (RF), parallel random forest (PRF), regularized random forest (RRF), and extremely randomized trees (ERT). Fifteen climatic and geo-environmental variables were used as inputs of the flash flood susceptibility models. The results showed that ERT was the most optimal model with an area under curve (AUC) value of 0.82. The rest of the models’ AUC values, i.e., RRF, PRF, RF, and BRT, were 0.80, 0.79, 0.78, and 0.75, respectively. In the ERT model, the areal coverage for very high to moderate flash flood susceptible area was 582.56 km2 (28.33%), and the rest of the portion was associated with very low to low susceptibility zones. It is concluded that topographical and hydrological parameters, e.g., altitude, slope, rainfall, and the river’s distance, were the most effective parameters. The results of this study will play a vital role in the planning and implementation of flood mitigation strategies in the region.

This study aims to evaluate a new approach in modeling gully erosion susceptibility based on deep learning neural network (DLNN) model, ensemble Particle swarm optimization (PSO) algorithm with DLNN (PSO-DLNN) and comparing these approaches with common artificial neural network (ANN) and support vector machine (SVM) models in Shiran watershed, Iran. For this purpose, 13 independent variables affecting gully erosion susceptibility in the study area, including altitude, slope, aspect, plan curvature, profile curvature, drainage density, distance from river, land use, soil, lithology, rainfall, , stream power index (SPI), topographic wetness index (TWI), were prepared. Also, 132 gully erosion locations were identified during field visits. Data for modeling were divided into two categories of training (70%) and testing (30%). Receiver operating characteristic (ROC) parameters including sensitivity, specificity, negative predictive value (NPV), positive predictive value (PPV) and area under curve (AUC) were used to evaluate the performance of the models. The results showed that, the AUC values from ROC with considering testing datasets of PSO-DLNN is 0.89 and which is associated with superb accuracy. Rest of the models also associated with optimal accuracy and near about PSO-DLNN model; the AUC values from ROC of DLNN, SVM and ANN for testing datasets are 0.87, 0.85 and 0.84 respectively. The PSO algorithm has updated and optimized the weights of DLNN model, and as a result, the efficiency of this model in predicting gully erosion susceptibility has increased. Therefore, it can be concluded that the use of DLNN model and its ensemble with PSO algorithm can be used as a novel and practical method in predicting the susceptibility of gully erosion that helps planners and managers in managing and reducing the risk of this phenomenon.

The estimation and mapping of forest stand characteristics are vital because this information is necessary for sustainable forest management. The present study considers the use of a Bayesian additive regression trees (BART) algorithm as a non-parametric classifier using Sentinel-2A data and topographic variables to estimate the forest stand characteristics, namely the basal area (m2/ha), stem volume (m3/ha), and stem density (number/ha). These results were compared with those of three other popular machine learning (ML) algorithms, such as generalised linear model (GLM), K-nearest neighbours (KNN), and support vector machine (SVM). A feature selection was done on 28 variables including the multi-spectral bands on Sentinel-2 satellite, related vegetation indices, and ancillary data (elevation, slope, and topographic solar-radiation index derived from digital elevation model (DEM)) and then the most insignificant variables were removed from the datasets by recursive feature elimination (RFE). The study area was a mountainous forest with high biodiversity and an elevation gradient from 26 to 1636 m. An inventory dataset of 1200 sample plots was provided for training and testing the algorithms, and the predictors were fed into the ML models to compute and predict the forest stand characteristics. The accuracies and certainties of the ML models were assessed by their root mean square error (RMSE), mean absolute error (MAE), and R-squared (R2) values. The results demonstrated that BART generated the best basal area and stem volume predictions, followed by GLM, SVM, and KNN. The best RMSE values for both basal area (8.12 m2/ha) and stem volume (29.28 m3/ha) estimation were obtained by BART. Thus, the ability of the BART model for forestry application was established. On the other hand, KNN exhibited the highest RMSE values for all stand variable predictions, thereby exhibiting the least accuracy for this specific application. Moreover, the effectiveness of the narrow Sentinel-2 bands around the red edge and elevation was highlighted for predicting the forest stand characteristics. Therefore, we concluded that the combination of the Sentinel-2 products and topographic variables derived from the PALSAR data used in this study improved the estimation of the forest attributes in temperate forests.

The objective of this study is to assess the gully head-cut erosion susceptibility and identify gully erosion prone areas in the Meimand watershed, Iran. In recent years, this study area has been greatly influenced by several head-cut gullies due to unusual climatic factors and human induced activity. The present study is therefore intended to address this issue by developing head-cut gully erosion prediction maps using boosting ensemble machine learning algorithms, namely Boosted Tree (BT), Boosted Generalized Linear Models (BGLM), Boosted Regression Tree (BRT), Extreme Gradient Boosting (XGB), and Deep Boost (DB). Initially, we produced a gully erosion inventory map using a variety of resources, including published reports, Google Earth images, and field records of the Global Positioning System (GPS). Subsequently, we distributed this information randomly and choose 70% (102) of the test gullies and the remaining 30% (43) for validation. The methodology was designed using morphometric and thematic determinants, including 14 head-cut gully erosion conditioning features. We have also investigated the following: (a) Multi-collinearity analysis to determine the linearity of the independent variables, (b) Predictive capability of piping models using train and test dataset and (c) Variables importance affecting head-cut gully erosion. The study reveals that altitude, land use, distances from road and soil characteristics influenced the method with the greatest impact on head-cut gully erosion susceptibility. We presented five head-cut gully erosion susceptibility maps and investigated their predictive accuracy through area under curve (AUC). The AUC test reveals that the DB machine learning method demonstrated significantly higher accuracy (AUC = 0.95) than the BT (AUC = 0.93), BGLM (AUC = 0.91), BRT (AUC = 0.94) and XGB (AUC = 0.92) approaches. The predicted head-cut gully erosion susceptibility maps can be used by policy makers and local authorities for soil conservation and to prevent threats to human activities.

In the present study, gully erosion susceptibility was evaluated for the area of the Robat Turk Watershed in Iran. The assessment of gully erosion susceptibility was performed using four state-of-the-art data mining techniques: random forest (RF), credal decision trees (CDTree), kernel logistic regression (KLR), and best-first decision tree (BFTree). To the best of our knowledge, the KLR and CDTree algorithms have been rarely applied to gully erosion modeling. In the first step, from the 242 gully erosion locations that were identified, 70% (170 gullies) were selected as the training dataset, and the other 30% (72 gullies) were considered for the result validation process. In the next step, twelve gully erosion conditioning factors, including topographic, geomorphological, environmental, and hydrologic factors, were selected to estimate gully erosion susceptibility. The area under the ROC curve (AUC) was used to estimate the performance of the models. The results revealed that the RF model had the best performance (AUC = 0.893), followed by the KLR (AUC = 0.825), the CDTree (AUC = 0.808), and the BFTree (AUC = 0.789) models. Overall, the RF model performed significantly better than the others, which may support the application of this method to a transferable susceptibility model in other areas. Therefore, we suggest using the RF, KLR, and CDT models for gully erosion susceptibility mapping in other prone areas to assess their reproducibility.

Shallow landslides through land degrading not only lead to threat the properly and life of human but they also may produce huge ecosystem damages. The aim of this study was to compare the performance of two decision tree machine learning algorithms including classification and regression tree (CART) and reduced error pruning tree (REPTree) for shallow landslide susceptibility mapping in Bijar, Kurdistan province, Iran. We first used 20 conditioning factors and then they were tested by information gain ratio (IGR) technique to select the most important ones. We then constructed a geodatabase based on the selected factors along with a total of 111 landslide locations with a ratio of 80/20 (for calibration/validation). The performance of the models was checked by the true positive rate (TP Rate), false positive rate (FP Rate), precision, recall, F1-Measure, Kappa, mean absolute error, and area under the receiver operatic curve (AUC). Results of IGR specified that the slope angle and TWI had the most contribution to shallow landslide occurrence in the study area. Moreover, results concluded that although these models had a high goodness-of-fit and prediction accuracy, the CART model (AUC=0.856) outperformed the REPTree model (AUC=0.837). Therefore, the CART model can be used as a promising tool and also as a base classifier to hybrid with optimization algorithms and Meta classifiers for spatial prediction of shallow landslide-prone areas.

Soil erosion is threatening land sustainability. Piping erosion is one of the land degradation processes that lead to significant landscape and environmental changes, and request a proper mapping survey. The purpose of this study is to survey piping erosion susceptibility maps in Zarandeieh watershed of Markazi province using Random Forest (RF), Support Vector Machine (SVM), and Bayesian Generalized Linear Models (Bayesian GLM) machine learning methods. For this purpose, due to the influence of different physiographic, environmental and soil conditions on the development and formation of piping, 18 variables were considered for modeling the piping erosion sensitivity in Zarandieh watershed. Based on field surveys and aerial photographs, 152 points of piping erosion were identified in the studied area, 70% of which was used for modeling, and 30% for model validation. The area under curve (AUC) was used to evaluate the performance of the models used. The results of the pipping erosion susceptibility showed that all three RF, SVM and Bayesian GLM models, have a good performance in the validation stage such as the AUC shown with values of 0.9 for RF, 0.88 for SVM, and 0.87 for Bayesian GLM. Altitude, PH and Bulk density are the variables that had the most impact on the pipping erosion sensitivity in the study area. This result shows that topographical and soil chemical factors are responsible for the piping distribution in the Zarandieh watershed.

Flash floods are one of the most devastating natural hazards, they occur within a catchment (region) where the response time of the drainage basin is short. Identification of probable flash flood locations and development of accurate flash flood susceptibility maps are important for proper flash flood management of a region. With this objective, we proposed and compared several novel hybrid computational approaches of machine learning methods for flash flood susceptibility mapping, namely AdaBoostM1 based Credal Decision Tree (ABM-CDT), Bagging based Credal Decision Tree (Bag-CDT), Dagging based Credal Decision Tree (Dag-CDT), MultiBoostAB based Credal Decision Tree (MBAB-CDT), and single Credal Decision Tree (CDT). These models were applied at a catchment of Markazi state in Iran. About 320 past flash flood events and nine flash flood influencing factors, namely distance from rivers, aspect, elevation, slope, rainfall, distance from faults, soil, land use, and lithology were considered and analyzed for the development of flash flood susceptibility maps. Correlation based feature selection method was used to validate and select the important factors for modeling of flash floods. Based on this feature selection analysis, only eight factors (distance from rivers, aspect, elevation, slope, rainfall, soil, land use, and lithology) were selected for the modeling, where distance to rivers is the most important factor for modeling of flash flood in this area. Performance of the models was validated and compared by using several robust metrics such as statistical measures and Area Under the Receiver Operating Characteristic (AUC) curve. The results of this study suggested that ABM-CDT (AUC = 0.957) has the best predictive capability in terms of accuracy, followed by Dag-CDT (AUC = 0.947), MBAB-CDT (AUC = 0.933), Bag-CDT (AUC = 0.932), and CDT (0.900), respectively. The proposed methods presented in this study would help in the development of accurate flash flood susceptible maps of watershed areas not only in Iran but also other parts of the world.

Cultural heritage sites, particularly those in mountainous regions face serious threats as mountains are hazardous places and many of them are located on shifting tectonic plates and live under the threat of earthquakes and related activities. Alongside, mountains are also exposed to atmospheric interactions and rainfall plays an important role in land movement activities through seepage of rainwater in fragile structures. Moreover, gravity pushing along with land use and climate change induced changing rainfall patterns and modification of slope lead to several hazards in mountainous regions. Environmental degradation in terms of soil erosion, loss of forests, and agricultural products are common phenomena in mountainous regions due to various multi-hazard threats. Therefore, it is necessary to conserve and management of our mountain environments as it is essential to the survival of the global ecosystem. Thus, the current research article focused on multi-hazard susceptibility mapping and evaluation of its risk assessment in some of the famous cultural heritage sites in the eastern Himalayan region of Sikkim state, India. Multi-hazard susceptibility mapping was carried out using boosted regression tree (BRT), Bayesian additive regression tree (BART) and Bayesian generalized linear model (BGLM) considering twenty-two conditioning factors and seismic activity, as this region is highly susceptible to earthquakes. The future climate and land use change were estimated using four representative concentration pathway (RCP) scenarios and Dynamic Conversion of Land-Use and its Effects (Dyna-CLUE) model respectively to identify future multi-hazard susceptibility areas and vulnerable cultural heritage sites. The novelty of this study is to a combination of machine learning, RCPs derived future climate and Dyna-CLUE induced future land use change estimation for multi-hazard modelling and identification of vulnerable cultural heritage sites. The result of this study will help land use planners and archaeologists to adopt proper management strategies for protecting the cultural heritage sites and maintaining environmental sustainability for the proper management of mountain resources.

This research was conducted to determine which areas in the Robat Turk watershed in Iran are sensitive to gully erosion, and to define the relationship between gully erosion and geo-environmental factors by two data mining techniques, namely, Random Forest (RF) and k-Nearest Neighbors (KNN). First, 242 gully locations we determined in field surveys and mapped in ArcGIS software. Then, twelve gully-related conditioning factors were selected. Our results showed that, for both the RF and KNN models, altitude, distance to roads, and distance from the river had the highest influence upon gully erosion sensitivity. We assessed the gully erosion susceptibility maps using the Receiver Operating Characteristic (ROC) curve. Validation results showed that the RF and KNN models had Area Under the Curve (AUC) 87.4 and 80.9%, respectively. As a result, the RF method has better performance compared with the KNN method for mapping gully erosion susceptibility. Rainfall, altitude, and distance from a river were identified as the most important factors affecting gully erosion in this area. The methodology used in this research is transferable to other regions to determine which areas are prone to gully erosion and to explicitly delineate high-risk zones within these areas.

Soil erosion determines landforms, soil formation and distribution, soil fertility, and land degradation processes. In arid and semiarid ecosystems, soil erosion is a key process to understand, foresee, and prevent desertification. Addressing soil erosion throughout watersheds scales requires basic information to develop soil erosion control strategies and to reduce land degradation. To assess and remediate the non-sustainable soil erosion rates, restoration programs benefit from the knowledge of the spatial distribution of the soil losses to develop maps of soil erosion. This study presents Support Vector Machine (SVM), Random Forest (RF), and adaptive boosting (AdaBoost) data mining models to map soil erosion susceptibility in Kozetopraghi watershed, Iran. A soil erosion inventory map was prepared from field rainfall simulation experiments on 174 randomly selected points along the Kozetopraghi watershed. In previous studies, this map has been prepared using indirect methods such as the Universal Soil Loss Equation to assess soil erosion. Direct field measurements for mapping soil erosion susceptibility have so far not been carried out in our study site in the past. The soil erosion rate data generated by simulated rainfall in 1 m2 plots at rainfall rate of 40 mmh&minus, 1 was used to develop the soil erosion map. Of the available data, 70% and 30% were randomly classified to calibrate and validate the models, respectively. As a result, the RF model with the highest area under the curve (AUC) value in a receiver operating characteristics (ROC) curve (0.91), and the lowest mean square error (MSE) value (0.09), has the most concordance and spatial differentiation. Sensitivity analysis by Jackknife and IncNodePurity methods indicates that the slope angle is the most important factor within the soil erosion susceptibility map. The RF susceptibility map showed that the areas located in the center and near the watershed outlet have the most susceptibility to soil erosion. This information can be used to support the development of sustainable restoration plans with more accuracy. Our methodology has been evaluated and can be also applied in other regions.