Soil water erosion (SWE) is an important global hazard that affects food availability through soil degradation, a reduction in crop yield, and agricultural land abandonment. A map of soil erosion susceptibility is a first and vital step in land management and soil conservation. Several machine learning (ML) algorithms optimized using the Grey Wolf Optimizer (GWO) metaheuristic algorithm can be used to accurately map SWE susceptibility. These optimized algorithms include Convolutional Neural Networks (CNN and CNN-GWO), Support Vector Machine (SVM and SVM-GWO), and Group Method of Data Handling (GMDH and GMDH-GWO). Results obtained using these algorithms can be compared with the well-known Revised Universal Soil Loss Equation (RUSLE) empirical model and Extreme Gradient Boosting (XGBoost) ML tree-based models. We apply these methods together with the frequency ratio (FR) model and the Information Gain Ratio (IGR) to determine the relationship between historical SWE data and controlling geo-environmental factors at 116 sites in the Noor-Rood watershed in northern Iran. Fourteen SWE geo-environmental factors are classified in topographical, hydro-climatic, land cover, and geological groups. We next divided the SWE sites into two datasets, one for model training (70% of the samples = 81 locations) and the other for model validation (30% of the samples = 35 locations). Finally the model-generated maps were evaluated using the Area under the Receiver Operating Characteristic (AU-ROC) curve. Our results show that elevation and rainfall erosivity have the greatest influence on SWE, while soil texture and hydrology are less important. The CNN-GWO model (AU-ROC = 0.85) outperformed other models, specifically, and in order, SVR-GWO = GMDH-GWO (AUC = 0.82), CNN = GMDH (AUC = 0.81), SVR = XGBoost (AUC = 0.80), and RULSE. Based on the RUSLE model, soil loss in the Noor-Rood watershed ranges from 0 to 2644 t ha–1yr1.

Keywords