Transformation of Geospatial Modelling of Soil Erosion Susceptibility Using Machine Learning

Transformation of Geospatial Modelling of Soil Erosion Susceptibility Using Machine Learning | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Transformation of Geospatial Modelling of Soil Erosion Susceptibility Using Machine Learning Muhammad Ramdhan Olii, Sartan Nento, Nurhayati Doda, Rizky Selly Nazarina Olii, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4933265/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Soil erosion creates substantial environmental and economic challenges, especially in areas vulnerable to land degradation. This study investigates the use of machine learning (ML) techniques—namely Support Vector Machines (SVM) and Generalized Linear Models (GLM)—for geospatial modeling of soil erosion susceptibility (SES). By leveraging geospatial data and incorporating a range of factors including hydrological, topographical, and environmental variables, the research aims to improve the accuracy and reliability of SES predictions. Results show that the SVM model predominantly identifies areas as having moderate (40.59%) or low (38.50%) susceptibility, whereas the GLM model allocates a higher proportion to very low (24.55%) and low (38.59%) susceptibility. Both models exhibit high performance, with SVM and GLM achieving accuracies of 87.4% and 87.2%, respectively, though GLM slightly surpasses AUC (0.939 vs. 0.916). GLM places greater emphasis on hydrological factors such as distance to rivers and drainage density, while SVM provides a more balanced assessment across various variables. This study demonstrates that ML-based models can significantly enhance SES assessments, offering a more nuanced and accurate approach than traditional methods. The findings highlight the value of adopting innovative, data-driven techniques in environmental modeling and offer practical insights for land management and conservation practices. Soil Erosion Susceptibility (SES) Geospatial Modelling Machine Learning (ML) Support Vector Machines (SVM) Generalized Linear Models (GLM) Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 1 Introduction Soil erosion remains a critical environmental challenge that impacts ecosystems, agriculture, and infrastructure (Olii et al., 2023 ). The degradation of fertile topsoil, sedimentation of waterways, and loss of vegetation cover are just some of the detrimental effects caused by erosion, which can lead to long-term ecological damage (Arabameri et al., 2019 ) and economic losses (Almouctar et al., 2021 ). Predicting and managing soil erosion susceptibility (SES) is therefore essential for sustainable land use and environmental conservation (Kucuker & Cedano Giraldo, 2022 ). To achieve this, accurate and reliable models that can predict the susceptibility of different areas to erosion are needed. Traditionally, modeling SES has relied on empirical methods that use historical data and simple statistical relationships to predict future erosion patterns (Saini et al., 2015 ). However, these methods often fail to capture the complexity of the interactions between the many environmental variables that influence erosion processes, such as rainfall intensity, soil type, land use, and topography (Olii, Olii, et al., 2024 ). As a result, there is a growing need for more sophisticated modeling approaches that can better account for these complexities and provide more accurate predictions (Golijanin et al., 2022 ; Kucuker & Cedano Giraldo, 2022 ). Geospatial modeling, when combined with machine learning techniques, significantly enhances the capability to predict and analyze SES. Machine learning (ML) models such as Random Forest (RF), Decision Tree (DT), Artificial Neural Networks (ANN), Support Vector Machines (SVM), Generalized Linear Models (GLM), etc can handle complex, non-linear relationships between environmental factors and SES, which traditional geospatial models might not fully capture (Al-Bawi et al., 2021 ; Gayen et al., 2019 ). By integrating these models with GIS and remote sensing data, the spatial patterns of SES can be more accurately mapped and predicted (Olii, Olii, et al., 2024 ). This combination allows for a more data-driven approach, where the models can learn from large datasets, adjust to various geographical contexts, and improve prediction accuracy by using classified and weighted factors tailored to local environmental conditions. The synergy between geospatial modeling and ML offers powerful tools for more effective land management and soil erosion prevention strategies. Most studies utilizing ML for SES modeling rely on raw or normalized continuous data, without prior classification into discrete classes or the assignment of weights based on expert judgment (Golkarian et al., 2023 ; Huang et al., 2023 ; Phinzi & Szabó, 2024 ). This approach can lead to less interpretable models, as the continuous nature of the data might obscure important distinctions between different categories of environmental factors. Additionally, the lack of expert-informed weights may result in the model underestimating or overestimating the significance of certain variables, potentially compromising the accuracy and robustness of predictions. This limitation highlights the need for more sophisticated methods that incorporate domain expertise into the ML modeling process to improve the reliability and practical applicability of SES assessments. This leaves a gap in understanding how the integration of this traditional approach with advanced ML models could enhance prediction accuracy and model interpretability. This study introduces a novel approach by integrating traditional classification and weighting of environmental factors with advanced ML models like SVM and GLM for SES mapping. Unlike previous studies that often apply these models individually, this research emphasizes the innovative combination of pre-classifying factors into discrete classes and assigning weights based on expert knowledge. This methodology not only enhances model interpretability, making the results more accessible to practitioners and decision-makers, but also addresses common machine learning challenges such as data complexity, overfitting, and multicollinearity. Furthermore, the adaptable classification system allows the model to be customized to various geographic settings, increasing its applicability and robustness. This study offers a significant advancement over traditional modeling approaches by improving prediction accuracy and model stability, particularly in regions with unique environmental conditions. This study aims to address this research gap by exploring the potential for integrating GLM and SVM in the geospatial modeling of SES. The research will involve a systematic comparison of these models to evaluate their effectiveness in predicting SES across different spatial scales and environmental conditions based on prior classification into discrete classes. Additionally, the study will investigate how these models can be optimized to handle the complexities of spatial data, including the influence of diverse and non-linear environmental variables such as climate, topography, and land use. By integrating GLM and SVM, this research seeks to develop a more robust and comprehensive framework for predicting erosion susceptibility, which could enhance the accuracy and reliability of SES assessments. The findings of this study have the potential to contribute significantly to the field of environmental modeling, offering new insights into the strengths and limitations of GLM and SVM and providing a basis for future research on the integration of statistical and ML approaches in geospatial modeling. 2. Materials and Method 2.1 Study Area The Saddang Watershed is located in the southwestern part of Sulawesi Island, within coordinates ranging from 2° 43' 42.4992" S to 3° 34' 51.4992" S latitude and 119° 14' 49.4988" E to 120° 3' 43.4988" E longitude. Covering an area of 4,909 km², as shown in Fig. 1 , this watershed spans across the South Sulawesi and West Sulawesi provinces. The primary river within the watershed is the Saddang River, which traverses the Enrekang, Tana Toraja, and North Toraja districts in South Sulawesi, and extends into Polewali in West Sulawesi. The river flows into the Makassar Strait through two estuaries: Barbana and Paria. The Saddang Watershed is crucial for both irrigation and energy production. The Benteng Dam provides irrigation for over 94,000 hectares of agricultural land, while the Bakaru Hydroelectric Power Plant, located downstream, has a capacity of 128 MW and plays a key role in meeting the region's energy needs. Additionally, the watershed has significant groundwater potential, estimated at around 1.354 million m 3 /year, which further supports local water supply and agriculture. The region's geomorphology features a variety of landforms, including fallen deposits, mountain/hill volcanoes, strongly incised folded mountains and hills, and karst hills. Land use within the watershed is varied, comprising settlements, rice fields, plantations, grasslands, swamps, water bodies, and areas designated for mixed dryland agriculture and forests. This mix of land uses supports a wide range of economic activities, including agriculture and forestry, which are vital to the local communities. The topography of the study area is notably diverse, with elevations ranging from 44 m to 2,880 m and an average elevation of 1,277 m. The climate in the Saddang Watershed is characterized by an average annual temperature of around 23°C. October is the warmest month with an average temperature of 26°C, while June is the coolest, averaging 22°C. Rainfall is substantial, with an annual average of 2,500 mm. The month of May receives the highest average rainfall of 387 mm, whereas September receives the least at 68 mm. 2.2 Overview of methodological framework The methodological framework of this study involves several key steps to assess SES (Fig. 2 ). It begins with data collection, followed by the mapping of soil erosion inventory. Key factors influencing soil erosion are then selected for analysis. ML models are used to predict SES areas, and their performance is rigorously evaluated. The results are then normalized to ensure consistency and comparability. This structured approach ensures a thorough and accurate assessment of SES. 2.2.1 Collection Data The data utilized in this study includes a diverse range of geospatial and environmental datasets, each sourced from specific platforms and offering distinct resolutions. The SRTM data, with a resolution of 30 x 30 m 2 , provides detailed elevation information crucial for topographic analysis and is accessible through the USGS Earth Explorer website ( https://earthexplorer.usgs.gov/ ). Similarly, Landsat 9 OLI/TIRS imagery, also with a 30 x 30 m 2 resolution, offers high-quality satellite imagery useful for quantifying vegetation greenness, available from the same USGS platform ( https://earthexplorer.usgs.gov/ ). For soil characteristics, the SoilGrids website ( https://soilgrids.org/ ) supplies maps with a 250 x 250 m 2 resolution detailing soil texture, organic carbon content, and bulk density, essential for understanding soil properties and their implications on land use and agriculture. Rainfall data, with a finer resolution of 0.25° x 0.25°, is provided by the NASA POWER Data Access Viewer ( https://power.larc.nasa.gov/data-access-viewer/ ), which helps in analyzing precipitation patterns and their impacts. SAS Planet and Google Earth Images are highly valuable for identifying and monitoring soil erosion. The detailed imagery enables precise visualization of land surface changes, allowing researchers to detect soil erosion patterns. By providing current and historical views of the landscape, these images support a comprehensive analysis of soil erosion processes. Lastly, administrative boundary data, provided in shapefile format, is available from the GADM website ( https://gadm.org/ ), which is instrumental for geographic and spatial analysis, delineating various administrative regions. Each dataset plays a critical role in providing comprehensive insights into the study area, supporting a wide range of analyses from elevation and land cover to soil properties and precipitation. 2.3.2 Soil Erosion Inventory Mapping The soil erosion inventory map is a key component in developing the SES model, serving as the dependent variable for this study. Accurately mapping the SES of the Saddang watershed required identifying both eroded and non-eroded areas. To facilitate this, the coordinates of 1992 locations—993 with soil erosion and 999 without—were collected through field surveys and analyzed using SAS Planet and Google Earth. These data points were then used to create a binary SES model, categorizing locations by soil erosion occurrence or non-occurrence. For model development, 1195 samples (60% of the total) were randomly selected, with the remaining 797 samples (40%) set aside for validation (Fig. 1 ). The types of soil erosion identified in the study included sheet, rill, gully, and mass movements. 2.3.3 Selection of the SES Factors The selection of factors for this study was carefully guided by several criteria: the availability of reliable data, insights from existing literature and prior research, the connectivity and variability of the data, and the specific geo-environmental characteristics of the study area. Based on these considerations, a comprehensive set of 11 key SES factors was identified and compiled. This set includes various hydrological factors such as rainfall erosivity, which measures the potential of rainfall to cause soil erosion, and the Topographical Wetness Index (TWI), which indicates areas of water accumulation. Additional hydrological factors include the distance to the river, the Stream Power Index (SPI), which assesses the energy of water flow, and drainage density, which reflects the network of water channels in the area. Topographic factors include the slope-length factor, which quantifies the impact of slope length on soil erosion, and the Topographic Roughness Index (TRI), which measures the variability of the terrain. Environmental factors encompass bulk density, which affects soil cohesion, clay ratio, which influences soil texture, soil organic carbon, which contributes to soil health, and the Normalized Difference Vegetation Index (NDVI), which indicates vegetation cover. The spatial distribution of SES factors can be seen in Fig. 3 . 2.3.3.1 Rainfall Erosivity Rainfall erosivity measures the potential of rainfall to cause soil erosion based on the intensity and kinetic energy of rainfall events. High rainfall erosivity indicates a greater potential for soil detachment and transport. The most common equation used to calculate rainfall erosivity is the R-factor in the Universal Soil Loss Equation (USLE): $$R=\sum\limits_{{i=1}}^{{12}} {1.735 \times {{10}^{\left( {1.5{{\log }_{10}}\left( {\frac{{{P_m}^{2}}}{{{P_a}}}} \right) - 0.018188} \right)}}}$$ 1 where R is the rainfall erosivity factor (MJ mm ha − 1 h − 1 year − 1 ), P m is the monthly rainfall (mm), and P a is the annual rainfall (mm). 2.3.3.2 Topographical Wetness Index (TWI) The Topographical Wetness Index indicates an area's susceptibility to soil saturation and water accumulation, which can influence soil erosion by increasing soil moisture and reducing its stability. High TWI values are often associated with greater soil erosion potential due to waterlogged conditions. It is calculated using the formula: $${\text{TWI}}=\ln \left( {\frac{{{A_s}}}{{\tan \beta }}} \right)$$ 2 where A s is the upstream contributing area and β is the slope gradient (in radians). 2.3.3.3 Stream Power Index (SPI) The Stream Power Index represents the erosive power of flowing water and its capacity to transport sediment. Higher SPI values suggest a higher potential for soil erosion due to the increased force exerted by flowing water. The value can be computed with the following formula: $${\text{SPI}}=\ln \left( {{A_s}\tan \beta } \right)$$ 3 where A s is the upstream contributing area and β is the slope gradient (in radians). 2.3.3.4 Distance to River The distance to a river significantly influences soil erosion by determining the likelihood of sediment transport into water bodies. Areas closer to rivers are more susceptible to soil erosion due to the increased potential for sediment movement and higher water flow, which can lead to both surface soil erosion and riverbank instability. As proximity to the river increases, the SES and subsequent sedimentation in the water body also rise, exacerbating soil erosion processes in these vulnerable zones. 2.3.3.5 Drainage Density Drainage density reflects the total length of streams and rivers per unit area in a watershed. Higher drainage density indicates a more dissected landscape, which can enhance surface runoff and increase soil erosion potential. It is calculated using the following expression: $${D_d}=\frac{{{l_s}}}{{{W_A}}}$$ 4 where l S is the total length of the river (km) and W A is a watershed area (km 2 ). 2.3.3.6 Slope Length Factor The Slope Length Factor represents the effect of slope length and steepness on soil erosion. Longer and steeper slopes generally contribute to greater soil erosion due to increased velocity and volume of surface runoff. It is commonly calculated using the LS factor in the USLE: \({\text{LS}}={\left( {\frac{\lambda }{{22.13}}} \right)^m}10.8\sin \beta +0.03\) if tan β < 0.09 (5) \({\text{LS}}={\left( {\frac{\lambda }{{22.13}}} \right)^m}16.8\sin \beta - 0.5\) if tan β ≥ 0.09 (6) $${\text{m}}=\frac{{\text{F}}}{{1+{\text{F}}}}$$ 7 $$F=\frac{{{{\sin \beta } \mathord{\left/ {\vphantom {{\sin \beta } {0.0896}}} \right. \kern-0pt} {0.0896}}}}{{3{{\left( {\sin \beta } \right)}^{0.8}}+0.56}}$$ 8 where λ is slope length (m), β is the slope gradient (in radians), m is the slope length exponent, and F is the ratio between rill soil erosion and interrill soil erosion. 2.3.3.7 Topographic Roughness Index (TRI) The Topographic Roughness Index measures the variability of terrain elevation, with higher values indicating a more rugged landscape. Rugged terrains often have more intense soil erosion due to increased surface runoff and reduced vegetation cover. It is calculated by the standard deviation of elevation within a specified window or grid cell using the equation: $${\text{TRI}}=Y{\left( {\sum {{{\left( {{x_{ij}} - {x_{00}}} \right)}^2}} } \right)^{0.5}}$$ 9 where x ij is the elevation of the neighbor grid (0,0). 2.3.3.8 Bulk Density Bulk density reflects the soil's compactness and porosity, influencing water infiltration and root penetration. High bulk density indicates compacted soil with lower infiltration rates, leading to increased surface runoff and potential soil erosion. Bulk density is typically measured directly using soil samples. 2.3.3.9 Clay Ratio The clay ratio indicates the proportion of clay particles in the soil, which affects soil structure and its susceptibility to soil erosion. Soils with higher clay content can be more cohesive, reducing soil erosion, but under certain conditions, they may also be susceptibility to crusting and soil erosion. The clay ratio is calculated as the percentage of clay particles relative to other soil particles. It is measured by applying: 10 where %clay is the percentage of clay content, %sand is the percentage of sand content, and %silt is the percentage of silt content. 2.3.3.10 Carbon Organic Soil organic carbon content influences soil structure, stability, and resistance to soil erosion. Higher organic carbon levels improve soil aggregation, reducing soil erosion by increasing infiltration and decreasing runoff. It is measured in units of decigrams per kilogram (dg/kg) of soil. 2.3.3.11 Normalized Difference Vegetation Index (NDVI) NDVI is a measure of vegetation cover and health, with higher values indicating more dense and vigorous vegetation. Dense vegetation protects soil from soil erosion by reducing the impact of raindrops and slowing surface runoff. NDVI is calculated using satellite imagery as: $${\text{NDVI}}=\frac{{{\text{NI}}{{\text{R}}_{{\text{band}}}} - {\text{Re}}{{\text{d}}_{{\text{band}}}}}}{{{\text{NI}}{{\text{R}}_{{\text{band}}}}+{\text{Re}}{{\text{d}}_{{\text{band}}}}}}$$ 11 Where NIR is light reflected in the near-infrared spectrum and RED band is light reflected in the red range of the spectrum. 2.3.4 Soil Erosion Modeling Using Machine Learning Models Soil erosion modeling using ML involves predicting the likelihood and extent of soil erosion based on various environmental factors. ML models can capture complex relationships between these factors and soil erosion processes, offering a flexible and data-driven approach to SES assessment. 2.3.4.1 Support Vector Machines (SVM) Support Vector Machines (SVM) are applied in soil erosion modeling to classify and predict areas at SES based on various hydrological, environmental, and topographical factors. SVM works by identifying the optimal hyperplane that separates different classes of SES with the maximum margin, effectively distinguishing between high-susceptibility and low- low-susceptibility zones. The model handles complex, non-linear relationships between factors by using kernel functions to map the input data into a higher-dimensional space where the classes become linearly separable. A common equation used in SVM classification is: $$f\left( x \right)=sign\left( {w\phi \left( x \right)+b} \right)$$ 12 where w is the weight vector, ϕ( x ) represents the transformation of the input data into a higher-dimensional space, and b is the bias term. The function f ( x ) determines the class label based on the sign of the output. By training on known soil erosion location and no-soil erosion location data, SVM can effectively predict SES in new, untested areas, making it a valuable tool in soil conservation and land management planning. 2.3.4.2 Generalized Linear Models (GLM) Generalized Linear Models (GLM) are used in soil erosion studies to analyze the relationship between soil erosion factors and soil erosion location or no-soil erosion location data by extending linear regression to handle non-normal distributions of the response variable, such as binary or count data. GLMs link predictors to soil erosion outcomes through a specified link function, allowing for the modeling of complex and non-linear relationships. The flexibility of GLM makes it suitable for predicting SES under varying hydrological, environmental, and topographical conditions. The general equation for a GLM is: $$g\left( \mu \right)={\beta _0}+{\beta _1}{X_1}+{\beta _2}{X_2}+ \cdots +{\beta _p}{X_p}$$ 13 where g ( µ ) is the link function that relates the mean of the response variable µ to the linear predictors β 0 , β 1 , …, β p and X 1 , X 2 , …, X p are the independent variables (SES factors). GLMs are valuable in understanding the probabilistic nature of soil erosion and predicting its occurrence under various scenarios. 2.3.5 Evaluating The Models’ Performance Several metrics and methods are employed to evaluate the performance of SVM and GLM in predicting SES, assessing the models' accuracy, reliability, and ability to generalize to new data, each serving to assess different aspects of model performance. These evaluations are crucial for determining how well the models perform, ensuring they provide reliable predictions across various scenarios and conditions. These evaluations are essential for determining the accuracy of the models in correctly identifying areas at SES, their reliability in providing consistent results across different datasets, and their generalization ability to make accurate predictions in new, unseen data. Metrics such as accuracy, precision, recall, and the F Mesure provide insights into the model’s ability to balance true positive (TP) and false positive rates (FPR), while the ROC curve and AUC measure the overall discriminative power of the model. Cross-validation techniques further ensure that the models are not overfitting to the training data, enhancing their robustness and reliability when applied to different environmental conditions. Together, these methods create a comprehensive framework for evaluating and refining the performance of SVM and GLM in SES modeling. 2.3.6 Normalization of Results Normalization of soil erosion susceptibility (SES) results involves scaling the output predictions to a common range to ensure consistent interpretation and comparison across different datasets or scenarios. This process adjusts the SES values ​​to a standardized scale, typically [0, 1], to make the results more comparable and to avoid biases introduced by varying magnitudes in the raw data. It can be measured by applying the equation: $${X_{norm}}=\frac{{X - {X_{\hbox{min} }}}}{{{X_{\hbox{max} }} - {X_{\hbox{min} }}}}$$ 14 where X is the original value, X min ​ is the minimum value in the dataset, and X max is the maximum value in the dataset. 2.3.6 Weighting and Scoring Modeling SES using weighted forms of SVM and GLM involves assigning weights to various factors to reflect their relative significance in predicting SES. In this approach, each factor is assigned a score based on its influence on SES (Table 2 ). These weights are then multiplied by factor-specific scores for different factor classes in Table 2 , ensuring that the contributions of each factor are accurately represented in the model's predictions. Table 2 Weights, Classes, and Scores of SES Factors No. SES Factors Categories Weight Classes of Factors Area (km 2 ) Area (%) Scores SVM GLM 1 Rainfall erosivity (MJ mm ha − 1 h − 1 year − 1 ) Hydrological Data 0.076 0.085 2,500 - - 5 2 Topographical Wetness Index (TWI) 0.093 0.063 20 2 0.0 5 3 Stream Power Index (SPI) 0.038 0.032 15 5 0.1 5 4 Distance to River (m) 0.107 0.376 > 1,600 2,777 56.6 1 1,200-1,600 454 9.3 2 800-1,200 492 10.0 3 400–800 540 11.0 4 < 400 645 13.1 5 5 Drainage Density (km/km 2 ) 0.134 0.282 0.0-0.2 2,011 41.0 1 0.2–0.4 1,815 37.0 2 0.4–0.6 925 18.8 3 0.6–0.8 150 3.1 4 0.8-1.0 8 0.2 5 6 Slope Length Factor Topographic Data 0.121 0.052 6.8 3,041 61.9 5 7 Topographic Roughness Index (TRI) 0.064 0.052 0.0-0.2 949 19.3 1 0.2–0.4 95 1.9 2 0.4–0.6 210 4.3 3 0.6–0.8 614 12.5 4 0.8-1.0 3,041 61.9 5 8 Bulk Density (cg/cm 3 ) Environmental Data 0.022 0.001 125 - - 5 9 Clay Ratio 0.099 0.143 0.0-0.2 - - 1 0.2–0.4 129 2.6 2 0.5–0.6 2,837 57.8 3 0.7–0.8 1,915 39.0 4 0.8-1.0 28 0.6 5 10 Carbon Organic (dg/kg) 0.088 0.098 > 125 62 1.3 1 100–125 1,865 38.0 2 75–100 2,175 44.3 3 50–75 788 16.0 4 0.7 - - 1 0.5–0.7 313 6.4 2 0.3–0.5 3,889 79.2 3 0.2–0.3 305 6.2 4 < 0.2 402 8.2 5 3 Results 3.1 Model Performance Evaluation The comparative analysis of SVM and GLM performance in predicting SES reveals nuanced differences that can guide model selection based on specific research goals. Table 3 , SVM shows a slight advantage in accuracy (87.4%) over GLM (87.2%), a crucial metric for overall model performance in correctly predicting both positive and negative outcomes. Accuracy remains a widely used measure in SES studies, indicating the proportion of correctly classified instances across the total cases. On the other hand, GLM excels in several critical performance areas. Figure 4 shows the higher Area AUC for GLM (0.939 vs. SVM's 0.916) suggesting that GLM is better at distinguishing between SES areas and those that are not, which is vital in developing reliable SES maps. Furthermore, GLM's higher recall and sensitivity (both at 89.4%) compared to SVM (both at 86.5%) indicate that GLM is more effective in identifying true positives, reducing the susceptibility of underestimating areas susceptibility to soil erosion—a key factor in environmental conservation and land management (Rahmati et al., 2017 ). Additionally, GLM's higher specificity (88.3% vs. SVM's 85%) underscores its effectiveness in correctly identifying non-susceptible areas, thereby minimizing false positives, which can lead to more targeted and cost-effective soil conservation strategies (Bui et al., 2020 ). Although SVM shows better precision (87.6% vs. GLM's 86.1%), indicating fewer false positives, the comprehensive performance of GLM across multiple metrics suggests that it may offer a more balanced approach for SES modeling, particularly when the goal is to minimize both false negatives and false positives. In conclusion, while SVM's slightly higher accuracy and precision may be advantageous in certain contexts, GLM's superior performance in AUC, recall, sensitivity, and specificity makes it a more robust choice for SES modeling, especially in scenarios where accurate identification of both high-susceptibility and low- susceptibility areas is crucial for sustainable land management and environmental protection. Table 3 Accuracy metric of the SES Model Accuracy Metrics Unit ML Algorithm SVM GLM Accuracy % 87.4 87.2 Classification Error % 12.6 12.8 AUC 0.916 0.939 Precision % 87.6 86.1 Recall % 86.5 89.4 F Measure % 87.0 87.7 Sensitivity % 86.5 89.4 Specificity % 85.0 88.3 3.2 Spatial Distribution of Soil Erosion Susceptibility (SES) Table 4 compares SES across different classes using the SVM and GLM presents interesting contrasts in their predictions. GLM categorizes a notably larger area under "Very Low" susceptibility (24.55%) compared to SVM (4.37%), indicating that GLM tends to classify a larger portion of the landscape as having minimal SES. This could be attributed to GLM's generalization tendencies, possibly smoothing over finer variations in the data, resulting in broader classifications. On the other hand, SVM's higher sensitivity to data nuances leads it to identify smaller, more concentrated areas of "Very Low" susceptibility, suggesting a more conservative approach. Additionally, the "Moderate" susceptibility category shows a significant difference, with SVM assigning 40.59% of the area to this class, while GLM assigns only 21.84%. This discrepancy further underscores SVM's tendency to distribute SES across a wider range of moderate susceptibility areas, possibly indicating its precision in capturing gradual variations in SES. In contrast, both models show remarkable alignment in the "Low" susceptibility category, with almost identical percentages (38.50% for SVM and 38.59% for GLM). This similarity suggests a consensus between the models in identifying regions with low SES, which might be attributed to clearer patterns or more robust data. However, the differences reemerge in the "High" and "Very High" classes, where SVM again shows a higher area in the "High" susceptibility category (15.56% compared to GLM's 13.32%), while GLM slightly exceeds SVM in the "Very High" class (1.23% versus 0.98%). This pattern could indicate that SVM is more sensitive to identifying areas that transition from moderate to high susceptibility, whereas GLM might be better at distinguishing extreme cases of susceptibility. Figure 5 of each model depicts that these highly vulnerable zones are primarily found along rivers and river branches. This pattern emerges because water movement in these areas frequently causes greater soil separation and transport. Because of their proximity to water bodies, soil particles are more likely to be swept away by surface runoff, resulting in increased soil erosion. Furthermore, the characteristics of river flow, such as meandering and the power of water during high discharge events, make these places more SES. The overall distribution suggests that SVM provides a more nuanced and detailed classification of SES, which could be useful for targeted management interventions, while GLM offers a broader, more generalized assessment that might be advantageous for broader, strategic planning and prioritization of soil erosion control measures. These differences have important implications for land management, where the choice of model might depend on the specific goals of the soil erosion assessment—whether it is for detailed site-specific management or general planning and resource allocation. Table 4 The SES class area of ​​each Machine Learning Algorithm No. Normalization Range Classes of Soil Erosion Susceptibility Support Vector Machines (SVM) Generalized Linear Models (GLM) Grid Total Area (km 2 ) % Grid Total Area (km 2 ) % 1 0.0–0.2 Very Low 238,436 215 4.4 1,339,038 1,205 24.6 2 0.2–0.4 Low 2,099,817 1,890 38.5 2,104,967 1,894 38.6 3 0.4–0.6 Moderate 2,213,906 1,993 40.6 1,191,423 1,072 21.8 4 0.6–0.8 High 848,573 764 15.6 726,729 654 13.3 5 0.8–1.0 Very High 53,583 48 1.0 67,158 60 1.2 Total 5,454,315 4,909 100 5,429,315 4,886 100 3.3 Feature Importance and Interpretation Table 2 provided highlights the weights assigned to various factors by each model, revealing insights into their relative significance. For SVM, the highest weights are assigned to distance to the river (0.107), drainage density (0.134), and slope length factor (0.121). These weights suggest that SVM places substantial emphasis on topographic and hydrological features. Recent studies support this prioritization. For example, Band et al. ( 2020 ) emphasized the importance of slope length and drainage density in predicting soil erosion due to their direct influence on runoff and soil erosion processes. Slope length affects the velocity of water flow, which is crucial in soil erosion, while drainage density impacts how water is distributed across the landscape (Vu Dinh et al., 2021 ). Additionally, the distance to rivers, which influences sediment transport and water flow patterns, is critical for SVM models, aligning with findings by Pourghasemi et al. ( 2020 ) that proximity to water bodies significantly affects soil erosion dynamics. In contrast, GLM assigns the highest weights to distance to river (0.376), drainage density (0.282), and NDVI (0.283). This prioritization indicates GLM's greater sensitivity to proximity to water bodies and land cover changes. The importance of distance to rivers and drainage density in GLM reflects their significant roles in soil erosion processes, as supported by Band et al. ( 2020 ), who found that proximity to rivers affects sediment transport and soil erosion, while drainage density influences runoff patterns. The Normalized Difference weight in GLM underscores the role of vegetation cover in SES. Recent research by Bai et al. (2020) highlights that NDVI is a critical factor in soil erosion modeling because it reflects vegetation cover, which stabilizes soil and reduces soil erosion rates. 4 Discussion 4.1 Strengths and Weaknesses of Each Model For SES modelling, the GLM offers several advantages due to its flexibility and interpretability. GLM allows for the inclusion of various types of predictor variables and can model different types of SES responses, such as binary or continuous outcomes. This interpretability is crucial in understanding the relationship between environmental factors and SES, which aids in making informed land management decisions. For instance, the GLM's higher emphasis on factors like distance to river and NDVI (with parameter weights of 0.376 and 0.283, respectively) aligns with their significant roles in soil erosion processes. These factors are influence sediment transport and land cover, as supported by studies, which emphasize their importance in soil erosion dynamics (Arabameri et al., 2019 ; Gayen et al., 2019 ; Igwe et al., 2020 ; Rahmati et al., 2017 ). However, GLM’s assumption of a linear relationship between predictors and response can be a limitation, especially in the context of the complex and non-linear processes governing soil erosion (Jiang et al., 2021 ). Moreover, GLM may struggle with high-dimensional data or significant interactions between variables, as seen in the moderate performance it shows in the accuracy metrics with an AUC of 0.939 compared to the SVM's 0.916. On the other hand, SVM are particularly well-suited for SES modelling due to their ability to handle non-linear relationships between predictors and response variables. SVM can effectively manage complex interactions and higher-dimensional data, making it advantageous in modelling the intricate processes involved in soil erosion. The SVM model’s emphasis on drainage density and slope length factor (with parameter weights of 0.134 and 0.121, respectively) reflects its strength in capturing topographic influences on soil erosion, which is consistent with the findings of (Olii et al., 2023 ). SVM's flexibility with different kernel functions allows it to adapt to various types of data distributions, enhancing its predictive accuracy in SES applications (Mustafa et al., 2018 ). However, SVM also has drawbacks, including being computationally intensive and less interpretable (Devos et al., 2009 ), which can be challenging for practitioners who need to understand the factors influencing SES. Despite these challenges, SVM's ability to provide accurate predictions is reflected in its performance metrics, such as a precision of 87.6% and recall of 86.5%, indicating its robustness in SES classification. 4.2 Implications for Spatial Susceptibility Assessment and Environmental Management The findings of SES assessments have significant implications for land management practices, particularly when integrating ML models like SVM and GLM with geospatial data. The data from the SES models, illustrated in the provided Table 4 , shows that both SVM and GLM effectively identify areas with varying degrees of SES, with GLM highlighting a larger area of very low susceptibility (24.55%) compared to SVM (4.37%). This precision in identifying high-susceptibility areas allows for more targeted interventions, such as reforestation, terracing, and the construction of check dams, which are crucial in preventing soil erosion and enhancing land productivity (Mosavi et al., 2020 ). Moreover, the data suggests that different models may emphasize different areas, as seen with the distribution of high and very high susceptibility areas, where SVM shows a higher percentage of moderate susceptibility (40.59%) compared to GLM (21.84%). These insights support findings that targeted soil conservation measures can lead to improved agricultural yields and reduced environmental degradation, demonstrating the practical benefits of using ML models in SES assessments (Folharini et al., 2023 ; Ghorbanzadeh et al., 2020 ; Nguyen et al., 2023 ; Wang et al., 2023 ). The application of such data-driven insights in land management optimizes resource allocation and aligns with sustainable development goals by mitigating the impact of soil erosion on ecosystems and human livelihoods (Arabameri et al., 2018 ; Rahmati et al., 2017 ). The broader environmental implications of spatial SES assessments extend beyond direct land management applications, particularly when considering the different weights and distributions identified by SVM and GLM models. For instance, SVM’s higher sensitivity in identifying moderately susceptible areas (40.59%) suggests that it may be more effective in predicting and managing soil erosion in regions where SES is not immediately apparent but still significant. In contrast, GLM’s emphasis on larger very low susceptibility areas (24.55%) can help prioritize areas for preventive measures. Accurate SES mapping is crucial for watershed management (Farhan et al., 2013 ), as it helps predict sediment loads in rivers and reservoirs, which is essential for maintaining water quality and preventing downstream flooding (Olii, et al., 2024 ). Additionally, the adaptability of these models to various geographic and environmental contexts, as indicated by the correlation between different model outputs and environmental factors, highlights their potential in global environmental management efforts. This adaptability is vital for integrating SES models into climate change adaptation strategies, predicting shifts in soil erosion patterns due to changing precipitation regimes and land use dynamics (Eekhout & de Vente, 2022 ). Ultimately, SES models, with their nuanced understanding of different susceptibility levels, provide valuable tools for combating desertification, preserving biodiversity, and supporting sustainable development on a global scale (Allafta & Opp, 2022 ; Das et al., 2020 ; François et al., 2024 ). 4.3 Limitations and Future Research Directions Despite advancements in modelling techniques, the study faces several limitations that could impact the accuracy and generalizability of the results. One significant limitation is the dependency on the quality and resolution of input data, which can directly affect the precision of model outputs (Olii et al., 2021 ). For instance, the models' ability to accurately classify areas into different SES classes, as shown in the provided Tables 3 , relies heavily on the quality of geospatial data. The SVM model, which identified 40.59% of the area as moderately susceptible, and the GLM model, which identified 21.84% in the same category, demonstrate how variations in data quality can influence model outcomes. The generalizability of these findings may also be constrained by the specific environmental and geographic characteristics of the study area, making it challenging to apply the results to regions with different climates, soil types, and land use patterns (Mosavi et al., 2020 ). Additionally, while SVM and GLM are robust models, their application may introduce biases if not properly calibrated, or if the classification and weighting systems, such as those seen in the different parameter weights for factors like distance to river and drainage density, are overly simplistic (Senanayake et al., 2020 ). These limitations highlight the importance of carefully interpreting the results and validating the models with independent datasets to ensure their reliability (Conoscenti et al., 2014 ). Future research should focus on addressing these limitations by exploring new avenues to enhance the accuracy and applicability of SES models. For instance, the use of higher-resolution spatial data could improve model performance by providing more detailed inputs, which would be particularly beneficial in refining the classification of areas with very high or very low SES, as shown by the differences in areas classified by SVM and GLM. Research could also investigate the effectiveness of more advanced ML models, such as Deep Learning or Neural Networks, which may better capture the complex interactions between environmental variables and soil erosion processes (Khosravi et al., 2023 ). Moreover, future studies should explore integrating real-time data from remote sensing technologies to enable dynamic monitoring of SES (Musasa et al., 2024 ). This approach could significantly enhance the responsiveness of land management strategies to ongoing environmental changes, aligning with the need for more adaptive and resilient frameworks for SES assessment and environmental management (Ejegu & Yegizaw, 2021 ). By tackling these areas, future research can contribute to developing more robust and comprehensive models that not only improve predictive accuracy but also broaden the applicability of SES assessments across diverse environmental contexts. Conclusions The integration of Machine Learning (ML) techniques such as Support Vector Machines (SVM) and Generalized Linear Models (GLM) in the geospatial modelling of Soil Erosion Susceptibility (SES) provides a significant enhancement in predictive accuracy and reliability over traditional methods. The analysis shows that the SVM model predominantly classifies areas as moderate (40.59%) and low (38.50%) susceptibility, while the GLM model identifies a larger portion of the area as very low (24.55%) and low (38.59%) susceptibility. Both models exhibit high accuracy metrics, with SVM achieving an accuracy of 87.4% and GLM close behind at 87.2%. Furthermore, the GLM model demonstrates a slightly better Area Under the Curve (AUC) score of 0.939, compared to 0.916 for the SVM model, indicating superior model performance in distinguishing between different classes of SES. The analysis of contributing factors highlights distinct parameter weights, with the GLM model placing greater emphasis on hydrological factors such as distance to rivers and drainage density, while the SVM model gives more balanced attention across multiple environmental variables. These differences suggest that GLM might offer a more nuanced understanding of hydrological influences on soil erosion. The robust performance of these ML models, validated through cross-validation techniques, underscores their potential for dynamic and accurate SES assessments. This is crucial for guiding effective land management and conservation strategies, particularly in regions vulnerable to soil erosion. These findings emphasize the value of adopting innovative, data-driven approaches in environmental modelling to address the complex challenges of land degradation and soil conservation. Declarations Acknowledgements We sincerely thank the Faculty of Engineering at Universitas Gorontalo for their invaluable support and resources, which were crucial to the success of this research. We also appreciate the contributions of our colleagues, research assistants, and collaborators, whose expertise and dedication significantly enhanced the quality of this study. Author Contributions Olii, M.R. led the conceptualization, methodology design, and overall supervision of the research project, and played a significant role in drafting and revising the manuscript. Nento, S. was responsible for data collection, statistical analysis, and interpretation of results, contributing to the writing of the results and discussion sections. Doda, N. conducted the literature review and developed the theoretical framework, assisting in writing the introduction and literature review sections. Djafar, H. handled fieldwork, data acquisition, and processing, and contributed to the development of environmental variables for the SES model. Olii, R.Z.N. supported software development, model simulation, and validation, and was involved in visualizing the results. Pakaya, P. contributed to writing, formatting the manuscript, and preparing tables and figures for publication. All authors reviewed and approved the final manuscript. Data availability The data supporting this study's findings are available from the corresponding author upon request, for academic and research purposes, subject to confidentiality agreements and data use policies. Conflicts of Interest The authors declare that there are no conflicts of interest regarding the publication of this paper References Al-Bawi, A. J., Al-Abadi, A. M., Pradhan, B., & Alamri, A. M. (2021). Assessing gully erosion susceptibility using topographic derived attributes, multi-criteria decision-making, and machine learning classifiers. Geomatics, Natural Hazards and Risk , 12 (1), 3035–3062. https://doi.org/10.1080/19475705.2021.1994024 Allafta, H., & Opp, C. (2022). Soil Erosion Assessment Using the RUSLE Model, Remote Sensing, and GIS in the Shatt Al-Arab Basin (Iraq-Iran). Applied Sciences (Switzerland) , 12 (15), 1–17. https://doi.org/10.3390/app12157776 Almouctar, M. A. S., Wu, Y., Zhao, F., & Dossou, J. F. (2021). Soil erosion assessment using the rusle model and geospatial techniques (Remote sensing and gis) in south-central niger (maradi region). Water (Switzerland) , 13 (24). https://doi.org/10.3390/w13243511 Arabameri, A., Cerda, A., & Tiefenbacher, J. P. (2019). Spatial pattern analysis and prediction of gully erosion using novel hybrid model of entropy-weight of evidence. Water (Switzerland) , 11 (6), 1–23. https://doi.org/10.3390/w11061129 Arabameri, A., Pradhan, B., Pourghasemi, H. R., & Rezaei, K. (2018). Identification of erosion-prone areas using different multi-criteria decision-making techniques and gis. Geomatics, Natural Hazards and Risk , 9 (1), 1129–1155. https://doi.org/10.1080/19475705.2018.1513084 Band, S. S., Janizadeh, S., Pal, S. C., Saha, A., Chakrabortty, R., Shokri, M., & Mosavi, A. (2020). Novel ensemble approach of deep learning neural network (Dlnn) model and particle swarm optimization (pso) algorithm for prediction of gully erosion susceptibility. Sensors (Switzerland) , 20 (19), 1–28. https://doi.org/10.3390/s20195609 Bui, D. T., Tsangaratos, P., Nguyen, V. T., Liem, N. Van, & Trinh, P. T. (2020). Comparing the prediction performance of a Deep Learning Neural Network model with conventional machine learning models in landslide susceptibility assessment. Catena , 188 (December 2019), 104426. https://doi.org/10.1016/j.catena.2019.104426 Conoscenti, C., Angileri, S., Cappadonia, C., Rotigliano, E., Agnesi, V., & Märker, M. (2014). Gully erosion susceptibility assessment by means of GIS-based logistic regression: A case of Sicily (Italy). Geomorphology , 204 , 399–411. https://doi.org/10.1016/j.geomorph.2013.08.021 Das, B., Bordoloi, R., Thungon, L. T., Paul, A., Pandey, P. K., Mishra, M., & Tripathi, O. P. (2020). An integrated approach of GIS, RUSLE and AHP to model soil erosion in West Kameng watershed, Arunachal Pradesh. Journal of Earth System Science , 129 (1), 1–18. https://doi.org/10.1007/s12040-020-1356-6 Devos, O., Ruckebusch, C., Durand, A., Duponchel, L., & Huvenne, J. P. (2009). Support vector machines (SVM) in near infrared (NIR) spectroscopy: Focus on parameters optimization and model interpretation. Chemometrics and Intelligent Laboratory Systems , 96 (1), 27–33. https://doi.org/10.1016/j.chemolab.2008.11.005 Eekhout, J. P. C., & de Vente, J. (2022). Global impact of climate change on soil erosion and potential for adaptation through soil conservation. Earth-Science Reviews , 226 , 103921. https://doi.org/10.1016/j.earscirev.2022.103921 Ejegu, M. A., & Yegizaw, E. S. (2021). Modeling soil erosion susceptibility and LULC dynamics for land degradation management using geoinformation technology in Debre Tabor district, Northwestern highlands of Ethiopia. J. Degrade. Min. Land Manage , 8 (2), 2623–2633. https://doi.org/10.15243/jdmlm.2021.082.2623 Farhan, Y., Zregat, D., & Farhan, I. (2013). Spatial Estimation of Soil Erosion Risk Using RUSLE Approach, RS, and GIS Techniques: A Case Study of Kufranja Watershed, Northern Jordan. In Journal of Water Resource and Protection (Vol. 05, Issue 12, pp. 1247–1261). https://doi.org/10.4236/jwarp.2013.512134 Folharini, S., Vieira, A., Bento-Gonçalves, A., Silva, S., Marques, T., & Novais, J. (2023). Soil Erosion Quantification using Machine Learning in Sub-Watersheds of Northern Portugal. Hydrology , 10 (1). https://doi.org/10.3390/hydrology10010007 François, M., Gonçalves Pontes, M. C., de Vasconcelos, R. N., de Oliveira, U. C., Peixoto da Silva, H., Faria, D., & Mariano-Neto, E. (2024). Assessing soil erosion and its drivers in agricultural landscapes: a case study in southern Bahia, Brazil. Journal of Water and Climate Change , 00 (0), 1–16. https://doi.org/10.2166/wcc.2024.147 Gayen, A., Pourghasemi, H. R., Saha, S., Keesstra, S., & Bai, S. (2019). Gully erosion susceptibility assessment and management of hazard-prone areas in India using different machine learning algorithms. Science of the Total Environment , 668 , 124–138. https://doi.org/10.1016/j.scitotenv.2019.02.436 Ghorbanzadeh, O., Shahabi, H., Mirchooli, F., Valizadeh Kamran, K., Lim, S., Aryal, J., Jarihani, B., & Blaschke, T. (2020). Gully erosion susceptibility mapping (GESM) using machine learning methods optimized by the multi‑collinearity analysis and K-fold cross-validation. Geomatics, Natural Hazards and Risk , 11 (1), 1653–1678. https://doi.org/10.1080/19475705.2020.1810138 Golijanin, J., Nikolić, G., Valjarević, A., Ivanović, R., Tunguz, V., Bojić, S., Grmuša, M., Lukić Tanović, M., Perić, M., Hrelja, E., & Stankov, S. (2022). Estimation of potential soil erosion reduction using GIS-based RUSLE under different land cover management models: A case study of Pale Municipality, B&H. Frontiers in Environmental Science , 10 (August), 1–13. https://doi.org/10.3389/fenvs.2022.945789 Golkarian, A., Khosravi, K., Panahi, M., & Clague, J. J. (2023). Spatial variability of soil water erosion: Comparing empirical and intelligent techniques. Geoscience Frontiers , 14 (1), 101456. https://doi.org/10.1016/j.gsf.2022.101456 Huang, D., Su, L., Zhou, L., Tian, Y., & Fan, H. (2023). Assessment of gully erosion susceptibility using different DEM-derived topographic factors in the black soil region of Northeast China. International Soil and Water Conservation Research , 11 (1), 97–111. https://doi.org/10.1016/j.iswcr.2022.04.001 Igwe, O., John, U. I., Solomon, O., & Obinna, O. (2020). GIS-based gully erosion susceptibility modeling, adapting bivariate statistical method and AHP approach in Gombe town and environs Northeast Nigeria. Geoenvironmental Disasters , 7 (1). https://doi.org/10.1186/s40677-020-00166-8 Jiang, C., Yang, Z., Wang, X., Dong, X., Li, Z., & Li, C. (2021). Examining the reversal of soil erosion decline in the hotspots of sandstorms: A non-linear ecosystem dynamic perspective. Journal of Arid Environments , 186 (October 2020), 104421. https://doi.org/10.1016/j.jaridenv.2020.104421 Khosravi, K., Rezaie, F., Cooper, J. R., Kalantari, Z., Abolfathi, S., & Hatamiafkoueieh, J. (2023). Soil water erosion susceptibility assessment using deep learning algorithms. Journal of Hydrology , 618 . https://doi.org/10.1016/j.jhydrol.2023.129229 Kucuker, D. M., & Cedano Giraldo, D. (2022). Assessment of soil erosion risk using an integrated approach of GIS and Analytic Hierarchy Process (AHP) in Erzurum, Turkiye. Ecological Informatics , 71 , 101788. https://doi.org/https://doi.org/10.1016/j.ecoinf.2022.101788 Mosavi, A., Sajedi-Hosseini, F., Choubin, B., Taromideh, F., Rahi, G., & Dineva, A. A. (2020). Susceptibility mapping of soil water erosion using machine learning models. Water (Switzerland) , 12 (7). https://doi.org/10.3390/w12071995 Musasa, T., Dube, T., & Marambanyika, T. (2024). Landsat satellite programme potential for soil erosion assessment and monitoring in arid environments: A review of applications and challenges. International Soil and Water Conservation Research , 12 (2), 267–278. https://doi.org/10.1016/j.iswcr.2023.10.003 Mustafa, M. R. U., Sholagberu, A. T., Yusof, K. W., Hashim, A. M., Khan, M. W. A., & Shahbaz, M. (2018). SVM-Based Geospatial Prediction of Soil Erosion under Static and Dynamic Conditioning Factors. MATEC Web of Conferences , 203 . https://doi.org/10.1051/matecconf/201820304004 Nguyen, C. Q., Tran, T. T., Nguyen, T. T. T., Nguyen, T. T., Astarkhanova, T. S., Vu, L. V., Dau, K. T., Nguyen, H. N., Pham, G. H., Nguyen, D. D., Prakash, I., & Pham, B. (2023). Mapping of soil erosion susceptibility using advanced machine learning models at Nghe An, Vietnam. Journal of Hydroinformatics , 26 (1). https://doi.org/10.2166/hydro.2023.327 Olii, M. R., Kironoto, B. A., Olii, A., Pakaya, R., & Olii, A. K. Z. (2024). Advancing Soil Erosion Assessment: Application of Remote Sensing and Geospatial Techniques in Bulango Ulu Reservoir Basin. E3S Web of Conferences , 476 , 1–15. https://doi.org/10.1051/e3sconf/202447601041 Olii, M. R., Olii, A. K. Z., Olii, A., Pakaya, R., & Kironoto, B. A. (2024). Spatial modeling of soil erosion risk a multi‑criteria decision‑making (MCDM) approach in the paguyaman watershed, gorontalo, Indo.pdf. Arabian Journal of Geosciences , 17 (226), 1–13. https://doi.org/https://doi.org/10.1007/s12517-024-12032-0 Olii, M. R., Olii, A., & Pakaya, R. (2021). The Integrated Spatial Assessment of The Flood Hazard Using AHP-GIS: The Case Study of Gorontalo Regency. Indonesian Journal of Geography , 53 (1), 126–135. https://doi.org/10.22146/ijg.59999 Olii, M. R., Olii, A., Pakaya, R., & Olii, M. Y. U. P. (2023). GIS-based analytic hierarchy process (AHP) for soil erosion-prone areas mapping in the Bone Watershed, Gorontalo, Indonesia. Environmental Earth Sciences , 82 (9), 1–14. https://doi.org/10.1007/s12665-023-10913-3 Phinzi, K., & Szabó, S. (2024). Predictive machine learning for gully susceptibility modeling with geo-environmental covariates: main drivers, model performance, and computational efficiency. In Natural Hazards (Vol. 120, Issue 8). Springer Netherlands. https://doi.org/10.1007/s11069-024-06481-9 Pourghasemi, H. R., Sadhasivam, N., Kariminejad, N., & Collins, A. L. (2020). Gully erosion spatial modelling: Role of machine learning algorithms in selection of the best controlling factors and modelling process. Geoscience Frontiers , 11 (6), 2207–2219. https://doi.org/10.1016/j.gsf.2020.03.005 Rahmati, O., Tahmasebipour, N., Haghizadeh, A., Pourghasemi, H. R., & Feizizadeh, B. (2017). Evaluating the influence of geo-environmental factors on gully erosion in a semi-arid region of Iran: An integrated framework. Science of the Total Environment , 579 , 913–927. https://doi.org/10.1016/j.scitotenv.2016.10.176 Saini, S. S., Jangra, R., & Kaushik, S. P. (2015). Vulnerability Assessment of Soil Erosion Using Geospatial Techniques-a Pilot Study of Upper Catchment of Markanda River. Research Article International Journal of Advancement in Remote Sensing , 3 (1), 9–21. Senanayake, S., Pradhan, B., Huete, A., & Brennan, J. (2020). A review on assessing and mapping soil erosion hazard using geo-informatics technology for farming system management. Remote Sensing , 12 (24), 1–25. https://doi.org/10.3390/rs12244063 Vu Dinh, T., Hoang, N. D., & Tran, X. L. (2021). Evaluation of Different Machine Learning Models for Predicting Soil Erosion in Tropical Sloping Lands of Northeast Vietnam. Applied and Environmental Soil Science , 2021 . https://doi.org/10.1155/2021/6665485 Wang, Y., Zhang, Y., & Chen, H. (2023). Gully erosion susceptibility prediction in Mollisols using machine learning models. Journal of Soil and Water Conservation , 78 (5), 385–396. https://doi.org/10.2489/jswc.2023.00019 Additional Declarations No competing interests reported. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4933265","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":347516244,"identity":"5b41ce27-633b-4870-934c-170b01402ae7","order_by":0,"name":"Muhammad Ramdhan Olii","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABAklEQVRIiWNgGAWjYBACCQY2MM24AUQ+MJCQA9EHHhCtJcFAwhisJYF4LQwMiQ1QBk4g2X4s8XNFzTbZ7exnDz5IKLBInx92+CHQFjs53QbsWqR50g5Lnjl223hnT16yAdBhuRtvpxkAtSQbmx3ArkWOIb1BsoHtduKGAzlmEmAtsxNAWg4kbsOlhf9588+Gf0At59+AtaQbzk7/gFeLtETaMcnGNqCWGxBbEuSlc/DbIjnjWZplY99t4w033hiD/GK4QTqn4ECCAW6/SJxPM77Z8O227IbzOYYPPvypk5efnb75w4cKOzlcWjCBAVilAbHKQUC+gRTVo2AUjIJRMBIAAKjLZz/86vdEAAAAAElFTkSuQmCC","orcid":"","institution":"Universitas Gorontalo","correspondingAuthor":true,"prefix":"","firstName":"Muhammad","middleName":"Ramdhan","lastName":"Olii","suffix":""},{"id":347516245,"identity":"56f39394-1d1c-4274-94d9-882d22ab2b60","order_by":1,"name":"Sartan Nento","email":"","orcid":"","institution":"Universitas Gorontalo","correspondingAuthor":false,"prefix":"","firstName":"Sartan","middleName":"","lastName":"Nento","suffix":""},{"id":347516246,"identity":"ed367cb2-b40c-43b1-8dbe-94ffee3bbbf4","order_by":2,"name":"Nurhayati Doda","email":"","orcid":"","institution":"Universitas Gorontalo","correspondingAuthor":false,"prefix":"","firstName":"Nurhayati","middleName":"","lastName":"Doda","suffix":""},{"id":347516247,"identity":"c2407efb-f6b7-4fec-9a3c-9b8c330a7a49","order_by":3,"name":"Rizky Selly Nazarina Olii","email":"","orcid":"","institution":"Universitas Gorontalo","correspondingAuthor":false,"prefix":"","firstName":"Rizky","middleName":"Selly Nazarina","lastName":"Olii","suffix":""},{"id":347516248,"identity":"b1909233-abf8-4e2a-a25c-14d939777ac3","order_by":4,"name":"Haris Djafar","email":"","orcid":"","institution":"Sulawesi II River Basin Center","correspondingAuthor":false,"prefix":"","firstName":"Haris","middleName":"","lastName":"Djafar","suffix":""},{"id":347516249,"identity":"742cb05f-0b9c-44ce-9102-7f1d04735e89","order_by":5,"name":"Ririn Pakaya","email":"","orcid":"","institution":"Universitas Gorontalo","correspondingAuthor":false,"prefix":"","firstName":"Ririn","middleName":"","lastName":"Pakaya","suffix":""}],"badges":[],"createdAt":"2024-08-18 12:11:12","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4933265/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4933265/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":64768945,"identity":"cb5f7798-49b1-4f48-afcc-33cb9e01b6ba","added_by":"auto","created_at":"2024-09-18 14:33:08","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1935548,"visible":true,"origin":"","legend":"\u003cp\u003eLocation of the study area\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-4933265/v1/a5a17b4a1ca40e6f4c503b53.png"},{"id":64768944,"identity":"3b5aab1e-aa2f-427b-82cf-a49d862df13e","added_by":"auto","created_at":"2024-09-18 14:33:08","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":978575,"visible":true,"origin":"","legend":"\u003cp\u003eResearch flowchart.\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-4933265/v1/2be3678c8e3207ccf99c6cb5.png"},{"id":64768948,"identity":"6676985f-88ed-4ccf-9459-04aa8a5e9662","added_by":"auto","created_at":"2024-09-18 14:33:08","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":15990310,"visible":true,"origin":"","legend":"\u003cp\u003eSES factors in Saddang Watershed\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-4933265/v1/414f6eae748bc3ffbfeb87de.png"},{"id":64770161,"identity":"9ba79ad1-5440-4c5d-893d-01a4a514b464","added_by":"auto","created_at":"2024-09-18 14:41:08","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":99758,"visible":true,"origin":"","legend":"\u003cp\u003eAUC-ROC Graph\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-4933265/v1/8c81f8da43c381700523b9c9.png"},{"id":64768947,"identity":"fe924a8c-ec59-4264-89ea-493eaafe42b2","added_by":"auto","created_at":"2024-09-18 14:33:08","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":2972880,"visible":true,"origin":"","legend":"\u003cp\u003eSoil erosion Susceptibility Maps Based on SVM [a] and GLM [b]\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-4933265/v1/54309983e28e1a46d36b3635.png"},{"id":73909638,"identity":"ab5bccff-880f-4693-931b-ffdce32a3d75","added_by":"auto","created_at":"2025-01-15 21:46:32","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":25714322,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4933265/v1/b709f05e-f8d4-4c45-8b9f-def3b92ff0fc.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Transformation of Geospatial Modelling of Soil Erosion Susceptibility Using Machine Learning","fulltext":[{"header":"1 Introduction","content":"\u003cp\u003eSoil erosion remains a critical environmental challenge that impacts ecosystems, agriculture, and infrastructure (Olii et al., \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The degradation of fertile topsoil, sedimentation of waterways, and loss of vegetation cover are just some of the detrimental effects caused by erosion, which can lead to long-term ecological damage (Arabameri et al., \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) and economic losses (Almouctar et al., \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Predicting and managing soil erosion susceptibility (SES) is therefore essential for sustainable land use and environmental conservation (Kucuker \u0026amp; Cedano Giraldo, \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). To achieve this, accurate and reliable models that can predict the susceptibility of different areas to erosion are needed. Traditionally, modeling SES has relied on empirical methods that use historical data and simple statistical relationships to predict future erosion patterns (Saini et al., \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). However, these methods often fail to capture the complexity of the interactions between the many environmental variables that influence erosion processes, such as rainfall intensity, soil type, land use, and topography (Olii, Olii, et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). As a result, there is a growing need for more sophisticated modeling approaches that can better account for these complexities and provide more accurate predictions (Golijanin et al., \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Kucuker \u0026amp; Cedano Giraldo, \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eGeospatial modeling, when combined with machine learning techniques, significantly enhances the capability to predict and analyze SES. Machine learning (ML) models such as Random Forest (RF), Decision Tree (DT), Artificial Neural Networks (ANN), Support Vector Machines (SVM), Generalized Linear Models (GLM), etc can handle complex, non-linear relationships between environmental factors and SES, which traditional geospatial models might not fully capture (Al-Bawi et al., \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Gayen et al., \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). By integrating these models with GIS and remote sensing data, the spatial patterns of SES can be more accurately mapped and predicted (Olii, Olii, et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). This combination allows for a more data-driven approach, where the models can learn from large datasets, adjust to various geographical contexts, and improve prediction accuracy by using classified and weighted factors tailored to local environmental conditions. The synergy between geospatial modeling and ML offers powerful tools for more effective land management and soil erosion prevention strategies.\u003c/p\u003e \u003cp\u003eMost studies utilizing ML for SES modeling rely on raw or normalized continuous data, without prior classification into discrete classes or the assignment of weights based on expert judgment (Golkarian et al., \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Huang et al., \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Phinzi \u0026amp; Szab\u0026oacute;, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). This approach can lead to less interpretable models, as the continuous nature of the data might obscure important distinctions between different categories of environmental factors. Additionally, the lack of expert-informed weights may result in the model underestimating or overestimating the significance of certain variables, potentially compromising the accuracy and robustness of predictions. This limitation highlights the need for more sophisticated methods that incorporate domain expertise into the ML modeling process to improve the reliability and practical applicability of SES assessments. This leaves a gap in understanding how the integration of this traditional approach with advanced ML models could enhance prediction accuracy and model interpretability. This study introduces a novel approach by integrating traditional classification and weighting of environmental factors with advanced ML models like SVM and GLM for SES mapping. Unlike previous studies that often apply these models individually, this research emphasizes the innovative combination of pre-classifying factors into discrete classes and assigning weights based on expert knowledge. This methodology not only enhances model interpretability, making the results more accessible to practitioners and decision-makers, but also addresses common machine learning challenges such as data complexity, overfitting, and multicollinearity. Furthermore, the adaptable classification system allows the model to be customized to various geographic settings, increasing its applicability and robustness. This study offers a significant advancement over traditional modeling approaches by improving prediction accuracy and model stability, particularly in regions with unique environmental conditions.\u003c/p\u003e \u003cp\u003eThis study aims to address this research gap by exploring the potential for integrating GLM and SVM in the geospatial modeling of SES. The research will involve a systematic comparison of these models to evaluate their effectiveness in predicting SES across different spatial scales and environmental conditions based on prior classification into discrete classes. Additionally, the study will investigate how these models can be optimized to handle the complexities of spatial data, including the influence of diverse and non-linear environmental variables such as climate, topography, and land use. By integrating GLM and SVM, this research seeks to develop a more robust and comprehensive framework for predicting erosion susceptibility, which could enhance the accuracy and reliability of SES assessments. The findings of this study have the potential to contribute significantly to the field of environmental modeling, offering new insights into the strengths and limitations of GLM and SVM and providing a basis for future research on the integration of statistical and ML approaches in geospatial modeling.\u003c/p\u003e"},{"header":"2. Materials and Method","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Study Area\u003c/h2\u003e \u003cp\u003eThe Saddang Watershed is located in the southwestern part of Sulawesi Island, within coordinates ranging from 2\u0026deg; 43' 42.4992\" S to 3\u0026deg; 34' 51.4992\" S latitude and 119\u0026deg; 14' 49.4988\" E to 120\u0026deg; 3' 43.4988\" E longitude. Covering an area of 4,909 km\u0026sup2;, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, this watershed spans across the South Sulawesi and West Sulawesi provinces. The primary river within the watershed is the Saddang River, which traverses the Enrekang, Tana Toraja, and North Toraja districts in South Sulawesi, and extends into Polewali in West Sulawesi. The river flows into the Makassar Strait through two estuaries: Barbana and Paria. The Saddang Watershed is crucial for both irrigation and energy production. The Benteng Dam provides irrigation for over 94,000 hectares of agricultural land, while the Bakaru Hydroelectric Power Plant, located downstream, has a capacity of 128 MW and plays a key role in meeting the region's energy needs. Additionally, the watershed has significant groundwater potential, estimated at around 1.354\u0026nbsp;million m\u003csup\u003e3\u003c/sup\u003e/year, which further supports local water supply and agriculture. The region's geomorphology features a variety of landforms, including fallen deposits, mountain/hill volcanoes, strongly incised folded mountains and hills, and karst hills. Land use within the watershed is varied, comprising settlements, rice fields, plantations, grasslands, swamps, water bodies, and areas designated for mixed dryland agriculture and forests. This mix of land uses supports a wide range of economic activities, including agriculture and forestry, which are vital to the local communities. The topography of the study area is notably diverse, with elevations ranging from 44 m to 2,880 m and an average elevation of 1,277 m. The climate in the Saddang Watershed is characterized by an average annual temperature of around 23\u0026deg;C. October is the warmest month with an average temperature of 26\u0026deg;C, while June is the coolest, averaging 22\u0026deg;C. Rainfall is substantial, with an annual average of 2,500 mm. The month of May receives the highest average rainfall of 387 mm, whereas September receives the least at 68 mm.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 Overview of methodological framework\u003c/h2\u003e \u003cp\u003eThe methodological framework of this study involves several key steps to assess SES (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). It begins with data collection, followed by the mapping of soil erosion inventory. Key factors influencing soil erosion are then selected for analysis. ML models are used to predict SES areas, and their performance is rigorously evaluated. The results are then normalized to ensure consistency and comparability. This structured approach ensures a thorough and accurate assessment of SES.\u003c/p\u003e \u003cdiv id=\"Sec5\" class=\"Section3\"\u003e \u003ch2\u003e2.2.1 Collection Data\u003c/h2\u003e \u003cp\u003eThe data utilized in this study includes a diverse range of geospatial and environmental datasets, each sourced from specific platforms and offering distinct resolutions. The SRTM data, with a resolution of 30 x 30 m\u003csup\u003e2\u003c/sup\u003e, provides detailed elevation information crucial for topographic analysis and is accessible through the USGS Earth Explorer website (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://earthexplorer.usgs.gov/\u003c/span\u003e\u003cspan address=\"https://earthexplorer.usgs.gov/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Similarly, Landsat 9 OLI/TIRS imagery, also with a 30 x 30 m\u003csup\u003e2\u003c/sup\u003e resolution, offers high-quality satellite imagery useful for quantifying vegetation greenness, available from the same USGS platform (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://earthexplorer.usgs.gov/\u003c/span\u003e\u003cspan address=\"https://earthexplorer.usgs.gov/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). For soil characteristics, the SoilGrids website (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://soilgrids.org/\u003c/span\u003e\u003cspan address=\"https://soilgrids.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) supplies maps with a 250 x 250 m\u003csup\u003e2\u003c/sup\u003e resolution detailing soil texture, organic carbon content, and bulk density, essential for understanding soil properties and their implications on land use and agriculture. Rainfall data, with a finer resolution of 0.25\u0026deg; x 0.25\u0026deg;, is provided by the NASA POWER Data Access Viewer (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://power.larc.nasa.gov/data-access-viewer/\u003c/span\u003e\u003cspan address=\"https://power.larc.nasa.gov/data-access-viewer/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), which helps in analyzing precipitation patterns and their impacts. SAS Planet and Google Earth Images are highly valuable for identifying and monitoring soil erosion. The detailed imagery enables precise visualization of land surface changes, allowing researchers to detect soil erosion patterns. By providing current and historical views of the landscape, these images support a comprehensive analysis of soil erosion processes. Lastly, administrative boundary data, provided in shapefile format, is available from the GADM website (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://gadm.org/\u003c/span\u003e\u003cspan address=\"https://gadm.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), which is instrumental for geographic and spatial analysis, delineating various administrative regions. Each dataset plays a critical role in providing comprehensive insights into the study area, supporting a wide range of analyses from elevation and land cover to soil properties and precipitation.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section3\"\u003e \u003ch2\u003e2.3.2 Soil Erosion Inventory Mapping\u003c/h2\u003e \u003cp\u003eThe soil erosion inventory map is a key component in developing the SES model, serving as the dependent variable for this study. Accurately mapping the SES of the Saddang watershed required identifying both eroded and non-eroded areas. To facilitate this, the coordinates of 1992 locations\u0026mdash;993 with soil erosion and 999 without\u0026mdash;were collected through field surveys and analyzed using SAS Planet and Google Earth. These data points were then used to create a binary SES model, categorizing locations by soil erosion occurrence or non-occurrence. For model development, 1195 samples (60% of the total) were randomly selected, with the remaining 797 samples (40%) set aside for validation (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). The types of soil erosion identified in the study included sheet, rill, gully, and mass movements.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section3\"\u003e \u003ch2\u003e2.3.3 Selection of the SES Factors\u003c/h2\u003e \u003cp\u003eThe selection of factors for this study was carefully guided by several criteria: the availability of reliable data, insights from existing literature and prior research, the connectivity and variability of the data, and the specific geo-environmental characteristics of the study area. Based on these considerations, a comprehensive set of 11 key SES factors was identified and compiled. This set includes various hydrological factors such as rainfall erosivity, which measures the potential of rainfall to cause soil erosion, and the Topographical Wetness Index (TWI), which indicates areas of water accumulation. Additional hydrological factors include the distance to the river, the Stream Power Index (SPI), which assesses the energy of water flow, and drainage density, which reflects the network of water channels in the area. Topographic factors include the slope-length factor, which quantifies the impact of slope length on soil erosion, and the Topographic Roughness Index (TRI), which measures the variability of the terrain. Environmental factors encompass bulk density, which affects soil cohesion, clay ratio, which influences soil texture, soil organic carbon, which contributes to soil health, and the Normalized Difference Vegetation Index (NDVI), which indicates vegetation cover. The spatial distribution of SES factors can be seen in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e.\u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.1 Rainfall Erosivity\u003c/h2\u003e \u003cp\u003eRainfall erosivity measures the potential of rainfall to cause soil erosion based on the intensity and kinetic energy of rainfall events. High rainfall erosivity indicates a greater potential for soil detachment and transport. The most common equation used to calculate rainfall erosivity is the R-factor in the Universal Soil Loss Equation (USLE):\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$$R=\\sum\\limits_{{i=1}}^{{12}} {1.735 \\times {{10}^{\\left( {1.5{{\\log }_{10}}\\left( {\\frac{{{P_m}^{2}}}{{{P_a}}}} \\right) - 0.018188} \\right)}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere R is the rainfall erosivity factor (MJ mm ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e h\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e year\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e), P\u003csub\u003em\u003c/sub\u003e is the monthly rainfall (mm), and P\u003csub\u003e\u003cem\u003ea\u003c/em\u003e\u003c/sub\u003e is the annual rainfall (mm).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.2 Topographical Wetness Index (TWI)\u003c/h2\u003e \u003cp\u003eThe Topographical Wetness Index indicates an area's susceptibility to soil saturation and water accumulation, which can influence soil erosion by increasing soil moisture and reducing its stability. High TWI values are often associated with greater soil erosion potential due to waterlogged conditions. It is calculated using the formula:\u003cdiv id=\"Equ2\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e\n$${\\text{TWI}}=\\ln \\left( {\\frac{{{A_s}}}{{\\tan \\beta }}} \\right)$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere A\u003csub\u003es\u003c/sub\u003e is the upstream contributing area and β is the slope gradient (in radians).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.3 Stream Power Index (SPI)\u003c/h2\u003e \u003cp\u003eThe Stream Power Index represents the erosive power of flowing water and its capacity to transport sediment. Higher SPI values suggest a higher potential for soil erosion due to the increased force exerted by flowing water. The value can be computed with the following formula:\u003cdiv id=\"Equ3\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ3\" name=\"EquationSource\"\u003e\n$${\\text{SPI}}=\\ln \\left( {{A_s}\\tan \\beta } \\right)$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e3\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere A\u003csub\u003es\u003c/sub\u003e is the upstream contributing area and β is the slope gradient (in radians).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.4 Distance to River\u003c/h2\u003e \u003cp\u003eThe distance to a river significantly influences soil erosion by determining the likelihood of sediment transport into water bodies. Areas closer to rivers are more susceptible to soil erosion due to the increased potential for sediment movement and higher water flow, which can lead to both surface soil erosion and riverbank instability. As proximity to the river increases, the SES and subsequent sedimentation in the water body also rise, exacerbating soil erosion processes in these vulnerable zones.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.5 Drainage Density\u003c/h2\u003e \u003cp\u003eDrainage density reflects the total length of streams and rivers per unit area in a watershed. Higher drainage density indicates a more dissected landscape, which can enhance surface runoff and increase soil erosion potential. It is calculated using the following expression:\u003cdiv id=\"Equ4\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ4\" name=\"EquationSource\"\u003e\n$${D_d}=\\frac{{{l_s}}}{{{W_A}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e4\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere l\u003csub\u003eS\u003c/sub\u003e is the total length of the river (km) and W\u003csub\u003eA\u003c/sub\u003e is a watershed area (km\u003csup\u003e2\u003c/sup\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.6 Slope Length Factor\u003c/h2\u003e \u003cp\u003eThe Slope Length Factor represents the effect of slope length and steepness on soil erosion. Longer and steeper slopes generally contribute to greater soil erosion due to increased velocity and volume of surface runoff. It is commonly calculated using the LS factor in the USLE:\u003c/p\u003e \u003cp\u003e \u003cspan class=\"InlineEquation\"\u003e \u003cspan class=\"mathinline\"\u003e\\({\\text{LS}}={\\left( {\\frac{\\lambda }{{22.13}}} \\right)^m}10.8\\sin \\beta +0.03\\)\u003c/span\u003e \u003c/span\u003e if tan β\u0026thinsp;\u0026lt;\u0026thinsp;0.09 (5)\u003c/p\u003e \u003cp\u003e \u003cspan class=\"InlineEquation\"\u003e \u003cspan class=\"mathinline\"\u003e\\({\\text{LS}}={\\left( {\\frac{\\lambda }{{22.13}}} \\right)^m}16.8\\sin \\beta - 0.5\\)\u003c/span\u003e \u003c/span\u003e if tan β\u0026thinsp;\u0026ge;\u0026thinsp;0.09 (6)\u003cdiv id=\"Equ5\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ5\" name=\"EquationSource\"\u003e\n$${\\text{m}}=\\frac{{\\text{F}}}{{1+{\\text{F}}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e7\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equ6\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ6\" name=\"EquationSource\"\u003e\n$$F=\\frac{{{{\\sin \\beta } \\mathord{\\left/ {\\vphantom {{\\sin \\beta } {0.0896}}} \\right. \\kern-0pt} {0.0896}}}}{{3{{\\left( {\\sin \\beta } \\right)}^{0.8}}+0.56}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e8\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere λ is slope length (m), β is the slope gradient (in radians), m is the slope length exponent, and F is the ratio between rill soil erosion and interrill soil erosion.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.7 Topographic Roughness Index (TRI)\u003c/h2\u003e \u003cp\u003eThe Topographic Roughness Index measures the variability of terrain elevation, with higher values indicating a more rugged landscape. Rugged terrains often have more intense soil erosion due to increased surface runoff and reduced vegetation cover. It is calculated by the standard deviation of elevation within a specified window or grid cell using the equation:\u003cdiv id=\"Equ7\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ7\" name=\"EquationSource\"\u003e\n$${\\text{TRI}}=Y{\\left( {\\sum {{{\\left( {{x_{ij}} - {x_{00}}} \\right)}^2}} } \\right)^{0.5}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e9\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003ex\u003c/em\u003e\u003csub\u003e\u003cem\u003eij\u003c/em\u003e\u003c/sub\u003e is the elevation of the neighbor grid (0,0).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.8 Bulk Density\u003c/h2\u003e \u003cp\u003eBulk density reflects the soil's compactness and porosity, influencing water infiltration and root penetration. High bulk density indicates compacted soil with lower infiltration rates, leading to increased surface runoff and potential soil erosion. Bulk density is typically measured directly using soil samples.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.9 Clay Ratio\u003c/h2\u003e \u003cp\u003eThe clay ratio indicates the proportion of clay particles in the soil, which affects soil structure and its susceptibility to soil erosion. Soils with higher clay content can be more cohesive, reducing soil erosion, but under certain conditions, they may also be susceptibility to crusting and soil erosion. The clay ratio is calculated as the percentage of clay particles relative to other soil particles. It is measured by applying:\u003cimg src=\"https://myfiles.space/user_files/122228_c8a1650c59388082/122228_custom_files/img1726669001.png\"\u003e\u003cdiv class=\"EquationNumber\"\u003e10\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere %clay is the percentage of clay content, %sand is the percentage of sand content, and %silt is the percentage of silt content.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.10 Carbon Organic\u003c/h2\u003e \u003cp\u003eSoil organic carbon content influences soil structure, stability, and resistance to soil erosion. Higher organic carbon levels improve soil aggregation, reducing soil erosion by increasing infiltration and decreasing runoff. It is measured in units of decigrams per kilogram (dg/kg) of soil.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section4\"\u003e \u003ch2\u003e2.3.3.11 Normalized Difference Vegetation Index (NDVI)\u003c/h2\u003e \u003cp\u003eNDVI is a measure of vegetation cover and health, with higher values indicating more dense and vigorous vegetation. Dense vegetation protects soil from soil erosion by reducing the impact of raindrops and slowing surface runoff. NDVI is calculated using satellite imagery as:\u003cdiv id=\"Equ9\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ9\" name=\"EquationSource\"\u003e\n$${\\text{NDVI}}=\\frac{{{\\text{NI}}{{\\text{R}}_{{\\text{band}}}} - {\\text{Re}}{{\\text{d}}_{{\\text{band}}}}}}{{{\\text{NI}}{{\\text{R}}_{{\\text{band}}}}+{\\text{Re}}{{\\text{d}}_{{\\text{band}}}}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e11\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eWhere NIR is light reflected in the near-infrared spectrum and RED\u003csub\u003eband\u003c/sub\u003e is light reflected in the red range of the spectrum.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section3\"\u003e \u003ch2\u003e2.3.4 Soil Erosion Modeling Using Machine Learning Models\u003c/h2\u003e \u003cp\u003eSoil erosion modeling using ML involves predicting the likelihood and extent of soil erosion based on various environmental factors. ML models can capture complex relationships between these factors and soil erosion processes, offering a flexible and data-driven approach to SES assessment.\u003c/p\u003e \u003cdiv id=\"Sec20\" class=\"Section4\"\u003e \u003ch2\u003e2.3.4.1 Support Vector Machines (SVM)\u003c/h2\u003e \u003cp\u003eSupport Vector Machines (SVM) are applied in soil erosion modeling to classify and predict areas at SES based on various hydrological, environmental, and topographical factors. SVM works by identifying the optimal hyperplane that separates different classes of SES with the maximum margin, effectively distinguishing between high-susceptibility and low- low-susceptibility zones. The model handles complex, non-linear relationships between factors by using kernel functions to map the input data into a higher-dimensional space where the classes become linearly separable. A common equation used in SVM classification is:\u003cdiv id=\"Equ10\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ10\" name=\"EquationSource\"\u003e\n$$f\\left( x \\right)=sign\\left( {w\\phi \\left( x \\right)+b} \\right)$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e12\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003ew\u003c/em\u003e is the weight vector, ϕ(\u003cem\u003ex\u003c/em\u003e) represents the transformation of the input data into a higher-dimensional space, and \u003cem\u003eb\u003c/em\u003e is the bias term. The function \u003cem\u003ef\u003c/em\u003e(\u003cem\u003ex\u003c/em\u003e) determines the class label based on the sign of the output. By training on known soil erosion location and no-soil erosion location data, SVM can effectively predict SES in new, untested areas, making it a valuable tool in soil conservation and land management planning.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec21\" class=\"Section4\"\u003e \u003ch2\u003e2.3.4.2 Generalized Linear Models (GLM)\u003c/h2\u003e \u003cp\u003eGeneralized Linear Models (GLM) are used in soil erosion studies to analyze the relationship between soil erosion factors and soil erosion location or no-soil erosion location data by extending linear regression to handle non-normal distributions of the response variable, such as binary or count data. GLMs link predictors to soil erosion outcomes through a specified link function, allowing for the modeling of complex and non-linear relationships. The flexibility of GLM makes it suitable for predicting SES under varying hydrological, environmental, and topographical conditions. The general equation for a GLM is:\u003cdiv id=\"Equ11\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ11\" name=\"EquationSource\"\u003e\n$$g\\left( \\mu \\right)={\\beta _0}+{\\beta _1}{X_1}+{\\beta _2}{X_2}+ \\cdots +{\\beta _p}{X_p}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e13\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003eg\u003c/em\u003e(\u003cem\u003e\u0026micro;\u003c/em\u003e) is the link function that relates the mean of the response variable \u003cem\u003e\u0026micro;\u003c/em\u003e to the linear predictors \u003cem\u003eβ\u003c/em\u003e\u003csub\u003e\u003cem\u003e0\u003c/em\u003e\u003c/sub\u003e, \u003cem\u003eβ\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e, \u0026hellip;, \u003cem\u003eβ\u003c/em\u003e\u003csub\u003ep\u003c/sub\u003e and \u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e, \u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e, \u0026hellip;, \u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003ep\u003c/em\u003e\u003c/sub\u003e are the independent variables (SES factors). GLMs are valuable in understanding the probabilistic nature of soil erosion and predicting its occurrence under various scenarios.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec22\" class=\"Section3\"\u003e \u003ch2\u003e2.3.5 Evaluating The Models\u0026rsquo; Performance\u003c/h2\u003e \u003cp\u003eSeveral metrics and methods are employed to evaluate the performance of SVM and GLM in predicting SES, assessing the models' accuracy, reliability, and ability to generalize to new data, each serving to assess different aspects of model performance. These evaluations are crucial for determining how well the models perform, ensuring they provide reliable predictions across various scenarios and conditions. These evaluations are essential for determining the accuracy of the models in correctly identifying areas at SES, their reliability in providing consistent results across different datasets, and their generalization ability to make accurate predictions in new, unseen data. Metrics such as accuracy, precision, recall, and the F Mesure provide insights into the model\u0026rsquo;s ability to balance true positive (TP) and false positive rates (FPR), while the ROC curve and AUC measure the overall discriminative power of the model. Cross-validation techniques further ensure that the models are not overfitting to the training data, enhancing their robustness and reliability when applied to different environmental conditions. Together, these methods create a comprehensive framework for evaluating and refining the performance of SVM and GLM in SES modeling.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec23\" class=\"Section3\"\u003e \u003ch2\u003e2.3.6 Normalization of Results\u003c/h2\u003e \u003cp\u003eNormalization of soil erosion susceptibility (SES) results involves scaling the output predictions to a common range to ensure consistent interpretation and comparison across different datasets or scenarios. This process adjusts the SES values ​​to a standardized scale, typically [0, 1], to make the results more comparable and to avoid biases introduced by varying magnitudes in the raw data. It can be measured by applying the equation:\u003cdiv id=\"Equ12\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ12\" name=\"EquationSource\"\u003e\n$${X_{norm}}=\\frac{{X - {X_{\\hbox{min} }}}}{{{X_{\\hbox{max} }} - {X_{\\hbox{min} }}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e14\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003eX\u003c/em\u003e is the original value, \u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003emin\u003c/em\u003e\u003c/sub\u003e ​ is the minimum value in the dataset, and \u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003emax\u003c/em\u003e\u003c/sub\u003e is the maximum value in the dataset.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec24\" class=\"Section3\"\u003e \u003ch2\u003e2.3.6 Weighting and Scoring\u003c/h2\u003e \u003cp\u003eModeling SES using weighted forms of SVM and GLM involves assigning weights to various factors to reflect their relative significance in predicting SES. In this approach, each factor is assigned a score based on its influence on SES (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e2\u003c/span\u003e). These weights are then multiplied by factor-specific scores for different factor classes in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e2\u003c/span\u003e, ensuring that the contributions of each factor are accurately represented in the model's predictions.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eWeights, Classes, and Scores of SES Factors\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"9\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c9\" colnum=\"9\"\u003e\u003c/div\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eNo.\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eSES\u003c/p\u003e \u003cp\u003eFactors\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eCategories\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colspan=\"2\" nameend=\"c5\" namest=\"c4\"\u003e \u003cp\u003eWeight\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eClasses of\u003c/p\u003e \u003cp\u003eFactors\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eArea\u003c/p\u003e \u003cp\u003e(km\u003csup\u003e2\u003c/sup\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eArea\u003c/p\u003e \u003cp\u003e(%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eScores\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eSVM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eGLM\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eRainfall erosivity\u003c/p\u003e \u003cp\u003e(MJ mm ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e h\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e year\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\" morerows=\"24\" rowspan=\"25\"\u003e \u003cp\u003eHydrological\u003c/p\u003e \u003cp\u003eData\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.076\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.085\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;1,750\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1,750-2,000\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e1,430\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e29.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e2,000\u0026ndash;2,250\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e1,951\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e39.7\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2,250-2,500\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e1,528\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e31.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026gt;\u0026thinsp;2,500\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eTopographical Wetness\u003c/p\u003e \u003cp\u003eIndex (TWI)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.093\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.063\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e1,960\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e39.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e5\u0026ndash;10\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e2,692\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e54.8\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e10\u0026ndash;15\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e227\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e4.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e15\u0026ndash;20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e28\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026gt;\u0026thinsp;20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eStream Power Index (SPI)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.038\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.032\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e19\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e0\u0026ndash;5\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e4,014\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e81.8\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e5\u0026ndash;10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e830\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e16.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e10\u0026ndash;15\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e41\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026gt;\u0026thinsp;15\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eDistance to River\u003c/p\u003e \u003cp\u003e(m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.107\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e\u003cb\u003e0.376\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e\u0026gt;\u0026thinsp;1,600\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e2,777\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e56.6\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1,200-1,600\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e454\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e9.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e800-1,200\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e492\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e10.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e400\u0026ndash;800\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e540\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e11.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;400\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e645\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e13.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eDrainage Density\u003c/p\u003e \u003cp\u003e(km/km\u003csup\u003e2\u003c/sup\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e\u003cb\u003e0.134\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.282\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e0.0-0.2\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e2,011\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e41.0\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.2\u0026ndash;0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e1,815\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e37.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.4\u0026ndash;0.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e925\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e18.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.6\u0026ndash;0.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e150\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e3.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.8-1.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eSlope Length Factor\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\" morerows=\"9\" rowspan=\"10\"\u003e \u003cp\u003eTopographic\u003c/p\u003e \u003cp\u003eData\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.121\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.052\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e949\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e19.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.4\u0026ndash;1.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e95\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e1.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1.4\u0026ndash;3.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e210\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e4.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e3.1\u0026ndash;6.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e614\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e12.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e\u0026gt;\u0026thinsp;6.8\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e3,041\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e61.9\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eTopographic Roughness\u003c/p\u003e \u003cp\u003eIndex (TRI)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.064\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.052\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.0-0.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e949\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e19.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.2\u0026ndash;0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e95\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e1.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.4\u0026ndash;0.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e210\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e4.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.6\u0026ndash;0.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e614\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e12.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e0.8-1.0\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e3,041\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e61.9\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eBulk Density\u003c/p\u003e \u003cp\u003e(cg/cm\u003csup\u003e3\u003c/sup\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\" morerows=\"19\" rowspan=\"20\"\u003e \u003cp\u003eEnvironmental\u003c/p\u003e \u003cp\u003eData\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.022\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.001\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e50\u0026ndash;75\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e18\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e75\u0026ndash;100\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e2,898\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e59.0\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e100\u0026ndash;125\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e1,993\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e40.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026gt;\u0026thinsp;125\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eClay Ratio\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.099\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.143\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.0-0.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.2\u0026ndash;0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e129\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e2.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e0.5\u0026ndash;0.6\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e2,837\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e57.8\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.7\u0026ndash;0.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e1,915\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e39.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.8-1.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e28\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eCarbon Organic\u003c/p\u003e \u003cp\u003e(dg/kg)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.088\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.098\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026gt;\u0026thinsp;125\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e62\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e1.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e100\u0026ndash;125\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e1,865\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e38.0\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e75\u0026ndash;100\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e2,175\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e44.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e50\u0026ndash;75\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e788\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e16.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e19\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e11\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003eNormalized Difference\u003c/p\u003e \u003cp\u003eVegetation Index (NDVI)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.108\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\" morerows=\"4\" rowspan=\"5\"\u003e \u003cp\u003e0.283\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026gt;\u0026thinsp;0.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.5\u0026ndash;0.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e313\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e6.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e0.3\u0026ndash;0.5\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e3,889\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e79.2\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.2\u0026ndash;0.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e305\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e6.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026lt;\u0026thinsp;0.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e402\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e8.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003c/div\u003e"},{"header":"3 Results","content":"\u003cdiv id=\"Sec26\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Model Performance Evaluation\u003c/h2\u003e \u003cp\u003eThe comparative analysis of SVM and GLM performance in predicting SES reveals nuanced differences that can guide model selection based on specific research goals. Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e3\u003c/span\u003e, SVM shows a slight advantage in accuracy (87.4%) over GLM (87.2%), a crucial metric for overall model performance in correctly predicting both positive and negative outcomes. Accuracy remains a widely used measure in SES studies, indicating the proportion of correctly classified instances across the total cases. On the other hand, GLM excels in several critical performance areas. Figure\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e shows the higher Area AUC for GLM (0.939 vs. SVM's 0.916) suggesting that GLM is better at distinguishing between SES areas and those that are not, which is vital in developing reliable SES maps. Furthermore, GLM's higher recall and sensitivity (both at 89.4%) compared to SVM (both at 86.5%) indicate that GLM is more effective in identifying true positives, reducing the susceptibility of underestimating areas susceptibility to soil erosion\u0026mdash;a key factor in environmental conservation and land management (Rahmati et al., \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Additionally, GLM's higher specificity (88.3% vs. SVM's 85%) underscores its effectiveness in correctly identifying non-susceptible areas, thereby minimizing false positives, which can lead to more targeted and cost-effective soil conservation strategies (Bui et al., \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Although SVM shows better precision (87.6% vs. GLM's 86.1%), indicating fewer false positives, the comprehensive performance of GLM across multiple metrics suggests that it may offer a more balanced approach for SES modeling, particularly when the goal is to minimize both false negatives and false positives.\u003c/p\u003e \u003cp\u003eIn conclusion, while SVM's slightly higher accuracy and precision may be advantageous in certain contexts, GLM's superior performance in AUC, recall, sensitivity, and specificity makes it a more robust choice for SES modeling, especially in scenarios where accurate identification of both high-susceptibility and low- susceptibility areas is crucial for sustainable land management and environmental protection.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eAccuracy metric of the SES Model\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eAccuracy Metrics\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eUnit\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"2\" nameend=\"c4\" namest=\"c3\"\u003e \u003cp\u003eML Algorithm\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSVM\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eGLM\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAccuracy\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e87.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e87.2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eClassification Error\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e12.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e12.8\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAUC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.916\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0.939\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePrecision\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e87.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e86.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRecall\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e86.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e89.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eF Measure\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e87.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e87.7\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSensitivity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e86.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e89.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSpecificity\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e85.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e88.3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec27\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Spatial Distribution of Soil Erosion Susceptibility (SES)\u003c/h2\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e4\u003c/span\u003e compares SES across different classes using the SVM and GLM presents interesting contrasts in their predictions. GLM categorizes a notably larger area under \"Very Low\" susceptibility (24.55%) compared to SVM (4.37%), indicating that GLM tends to classify a larger portion of the landscape as having minimal SES. This could be attributed to GLM's generalization tendencies, possibly smoothing over finer variations in the data, resulting in broader classifications. On the other hand, SVM's higher sensitivity to data nuances leads it to identify smaller, more concentrated areas of \"Very Low\" susceptibility, suggesting a more conservative approach. Additionally, the \"Moderate\" susceptibility category shows a significant difference, with SVM assigning 40.59% of the area to this class, while GLM assigns only 21.84%. This discrepancy further underscores SVM's tendency to distribute SES across a wider range of moderate susceptibility areas, possibly indicating its precision in capturing gradual variations in SES.\u003c/p\u003e \u003cp\u003eIn contrast, both models show remarkable alignment in the \"Low\" susceptibility category, with almost identical percentages (38.50% for SVM and 38.59% for GLM). This similarity suggests a consensus between the models in identifying regions with low SES, which might be attributed to clearer patterns or more robust data. However, the differences reemerge in the \"High\" and \"Very High\" classes, where SVM again shows a higher area in the \"High\" susceptibility category (15.56% compared to GLM's 13.32%), while GLM slightly exceeds SVM in the \"Very High\" class (1.23% versus 0.98%). This pattern could indicate that SVM is more sensitive to identifying areas that transition from moderate to high susceptibility, whereas GLM might be better at distinguishing extreme cases of susceptibility. Figure\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e of each model depicts that these highly vulnerable zones are primarily found along rivers and river branches. This pattern emerges because water movement in these areas frequently causes greater soil separation and transport. Because of their proximity to water bodies, soil particles are more likely to be swept away by surface runoff, resulting in increased soil erosion. Furthermore, the characteristics of river flow, such as meandering and the power of water during high discharge events, make these places more SES.\u003c/p\u003e \u003cp\u003eThe overall distribution suggests that SVM provides a more nuanced and detailed classification of SES, which could be useful for targeted management interventions, while GLM offers a broader, more generalized assessment that might be advantageous for broader, strategic planning and prioritization of soil erosion control measures. These differences have important implications for land management, where the choice of model might depend on the specific goals of the soil erosion assessment\u0026mdash;whether it is for detailed site-specific management or general planning and resource allocation.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 4\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eThe SES class area of ​​each Machine Learning Algorithm\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"9\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c9\" colnum=\"9\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eNo.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eNormalization Range\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eClasses of Soil Erosion Susceptibility\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"3\" nameend=\"c6\" namest=\"c4\"\u003e \u003cp\u003eSupport Vector Machines\u003c/p\u003e \u003cp\u003e(SVM)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"3\" nameend=\"c9\" namest=\"c7\"\u003e \u003cp\u003eGeneralized Linear Models\u003c/p\u003e \u003cp\u003e(GLM)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eGrid\u003c/p\u003e \u003cp\u003eTotal\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eArea\u003c/p\u003e \u003cp\u003e(km\u003csup\u003e2\u003c/sup\u003e)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eGrid\u003c/p\u003e \u003cp\u003eTotal\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003eArea\u003c/p\u003e \u003cp\u003e(km\u003csup\u003e2\u003c/sup\u003e)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c9\"\u003e \u003cp\u003e%\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.0\u0026ndash;0.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eVery Low\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e238,436\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e215\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e4.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e1,339,038\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e1,205\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e24.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.2\u0026ndash;0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eLow\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e2,099,817\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1,890\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e38.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e2,104,967\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e1,894\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e38.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.4\u0026ndash;0.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eModerate\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e2,213,906\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1,993\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e40.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e1,191,423\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e1,072\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e21.8\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.6\u0026ndash;0.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eHigh\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e848,573\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e764\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e15.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e726,729\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e654\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e13.3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.8\u0026ndash;1.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eVery High\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e53,583\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e48\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e67,158\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e60\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1.2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003eTotal\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e\u003cb\u003e5,454,315\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e4,909\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cb\u003e100\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e5,429,315\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e\u003cb\u003e4,886\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e\u003cb\u003e100\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec28\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Feature Importance and Interpretation\u003c/h2\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e2\u003c/span\u003e provided highlights the weights assigned to various factors by each model, revealing insights into their relative significance. For SVM, the highest weights are assigned to distance to the river (0.107), drainage density (0.134), and slope length factor (0.121). These weights suggest that SVM places substantial emphasis on topographic and hydrological features. Recent studies support this prioritization. For example, Band et al. (\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) emphasized the importance of slope length and drainage density in predicting soil erosion due to their direct influence on runoff and soil erosion processes. Slope length affects the velocity of water flow, which is crucial in soil erosion, while drainage density impacts how water is distributed across the landscape (Vu Dinh et al., \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Additionally, the distance to rivers, which influences sediment transport and water flow patterns, is critical for SVM models, aligning with findings by Pourghasemi et al. (\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) that proximity to water bodies significantly affects soil erosion dynamics.\u003c/p\u003e \u003cp\u003eIn contrast, GLM assigns the highest weights to distance to river (0.376), drainage density (0.282), and NDVI (0.283). This prioritization indicates GLM's greater sensitivity to proximity to water bodies and land cover changes. The importance of distance to rivers and drainage density in GLM reflects their significant roles in soil erosion processes, as supported by Band et al. (\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2020\u003c/span\u003e), who found that proximity to rivers affects sediment transport and soil erosion, while drainage density influences runoff patterns. The Normalized Difference weight in GLM underscores the role of vegetation cover in SES. Recent research by Bai et al. (2020) highlights that NDVI is a critical factor in soil erosion modeling because it reflects vegetation cover, which stabilizes soil and reduces soil erosion rates.\u003c/p\u003e \u003c/div\u003e"},{"header":"4 Discussion","content":"\u003cdiv id=\"Sec30\" class=\"Section2\"\u003e \u003ch2\u003e4.1 Strengths and Weaknesses of Each Model\u003c/h2\u003e \u003cp\u003eFor SES modelling, the GLM offers several advantages due to its flexibility and interpretability. GLM allows for the inclusion of various types of predictor variables and can model different types of SES responses, such as binary or continuous outcomes. This interpretability is crucial in understanding the relationship between environmental factors and SES, which aids in making informed land management decisions. For instance, the GLM's higher emphasis on factors like distance to river and NDVI (with parameter weights of 0.376 and 0.283, respectively) aligns with their significant roles in soil erosion processes. These factors are influence sediment transport and land cover, as supported by studies, which emphasize their importance in soil erosion dynamics (Arabameri et al., \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Gayen et al., \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Igwe et al., \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Rahmati et al., \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). However, GLM\u0026rsquo;s assumption of a linear relationship between predictors and response can be a limitation, especially in the context of the complex and non-linear processes governing soil erosion (Jiang et al., \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Moreover, GLM may struggle with high-dimensional data or significant interactions between variables, as seen in the moderate performance it shows in the accuracy metrics with an AUC of 0.939 compared to the SVM's 0.916.\u003c/p\u003e \u003cp\u003eOn the other hand, SVM are particularly well-suited for SES modelling due to their ability to handle non-linear relationships between predictors and response variables. SVM can effectively manage complex interactions and higher-dimensional data, making it advantageous in modelling the intricate processes involved in soil erosion. The SVM model\u0026rsquo;s emphasis on drainage density and slope length factor (with parameter weights of 0.134 and 0.121, respectively) reflects its strength in capturing topographic influences on soil erosion, which is consistent with the findings of (Olii et al., \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). SVM's flexibility with different kernel functions allows it to adapt to various types of data distributions, enhancing its predictive accuracy in SES applications (Mustafa et al., \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). However, SVM also has drawbacks, including being computationally intensive and less interpretable (Devos et al., \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2009\u003c/span\u003e), which can be challenging for practitioners who need to understand the factors influencing SES. Despite these challenges, SVM's ability to provide accurate predictions is reflected in its performance metrics, such as a precision of 87.6% and recall of 86.5%, indicating its robustness in SES classification.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec31\" class=\"Section2\"\u003e \u003ch2\u003e4.2 Implications for Spatial Susceptibility Assessment and Environmental Management\u003c/h2\u003e \u003cp\u003eThe findings of SES assessments have significant implications for land management practices, particularly when integrating ML models like SVM and GLM with geospatial data. The data from the SES models, illustrated in the provided Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e4\u003c/span\u003e, shows that both SVM and GLM effectively identify areas with varying degrees of SES, with GLM highlighting a larger area of very low susceptibility (24.55%) compared to SVM (4.37%). This precision in identifying high-susceptibility areas allows for more targeted interventions, such as reforestation, terracing, and the construction of check dams, which are crucial in preventing soil erosion and enhancing land productivity (Mosavi et al., \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Moreover, the data suggests that different models may emphasize different areas, as seen with the distribution of high and very high susceptibility areas, where SVM shows a higher percentage of moderate susceptibility (40.59%) compared to GLM (21.84%). These insights support findings that targeted soil conservation measures can lead to improved agricultural yields and reduced environmental degradation, demonstrating the practical benefits of using ML models in SES assessments (Folharini et al., \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Ghorbanzadeh et al., \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Nguyen et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Wang et al., \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The application of such data-driven insights in land management optimizes resource allocation and aligns with sustainable development goals by mitigating the impact of soil erosion on ecosystems and human livelihoods (Arabameri et al., \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Rahmati et al., \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2017\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe broader environmental implications of spatial SES assessments extend beyond direct land management applications, particularly when considering the different weights and distributions identified by SVM and GLM models. For instance, SVM\u0026rsquo;s higher sensitivity in identifying moderately susceptible areas (40.59%) suggests that it may be more effective in predicting and managing soil erosion in regions where SES is not immediately apparent but still significant. In contrast, GLM\u0026rsquo;s emphasis on larger very low susceptibility areas (24.55%) can help prioritize areas for preventive measures. Accurate SES mapping is crucial for watershed management (Farhan et al., \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2013\u003c/span\u003e), as it helps predict sediment loads in rivers and reservoirs, which is essential for maintaining water quality and preventing downstream flooding (Olii, et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). Additionally, the adaptability of these models to various geographic and environmental contexts, as indicated by the correlation between different model outputs and environmental factors, highlights their potential in global environmental management efforts. This adaptability is vital for integrating SES models into climate change adaptation strategies, predicting shifts in soil erosion patterns due to changing precipitation regimes and land use dynamics (Eekhout \u0026amp; de Vente, \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Ultimately, SES models, with their nuanced understanding of different susceptibility levels, provide valuable tools for combating desertification, preserving biodiversity, and supporting sustainable development on a global scale (Allafta \u0026amp; Opp, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Das et al., \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Fran\u0026ccedil;ois et al., \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec32\" class=\"Section2\"\u003e \u003ch2\u003e4.3 Limitations and Future Research Directions\u003c/h2\u003e \u003cp\u003eDespite advancements in modelling techniques, the study faces several limitations that could impact the accuracy and generalizability of the results. One significant limitation is the dependency on the quality and resolution of input data, which can directly affect the precision of model outputs (Olii et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). For instance, the models' ability to accurately classify areas into different SES classes, as shown in the provided Tables\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e3\u003c/span\u003e, relies heavily on the quality of geospatial data. The SVM model, which identified 40.59% of the area as moderately susceptible, and the GLM model, which identified 21.84% in the same category, demonstrate how variations in data quality can influence model outcomes. The generalizability of these findings may also be constrained by the specific environmental and geographic characteristics of the study area, making it challenging to apply the results to regions with different climates, soil types, and land use patterns (Mosavi et al., \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Additionally, while SVM and GLM are robust models, their application may introduce biases if not properly calibrated, or if the classification and weighting systems, such as those seen in the different parameter weights for factors like distance to river and drainage density, are overly simplistic (Senanayake et al., \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). These limitations highlight the importance of carefully interpreting the results and validating the models with independent datasets to ensure their reliability (Conoscenti et al., \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2014\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFuture research should focus on addressing these limitations by exploring new avenues to enhance the accuracy and applicability of SES models. For instance, the use of higher-resolution spatial data could improve model performance by providing more detailed inputs, which would be particularly beneficial in refining the classification of areas with very high or very low SES, as shown by the differences in areas classified by SVM and GLM. Research could also investigate the effectiveness of more advanced ML models, such as Deep Learning or Neural Networks, which may better capture the complex interactions between environmental variables and soil erosion processes (Khosravi et al., \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Moreover, future studies should explore integrating real-time data from remote sensing technologies to enable dynamic monitoring of SES (Musasa et al., \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). This approach could significantly enhance the responsiveness of land management strategies to ongoing environmental changes, aligning with the need for more adaptive and resilient frameworks for SES assessment and environmental management (Ejegu \u0026amp; Yegizaw, \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). By tackling these areas, future research can contribute to developing more robust and comprehensive models that not only improve predictive accuracy but also broaden the applicability of SES assessments across diverse environmental contexts.\u003c/p\u003e \u003c/div\u003e"},{"header":"Conclusions","content":"\u003cp\u003eThe integration of Machine Learning (ML) techniques such as Support Vector Machines (SVM) and Generalized Linear Models (GLM) in the geospatial modelling of Soil Erosion Susceptibility (SES) provides a significant enhancement in predictive accuracy and reliability over traditional methods. The analysis shows that the SVM model predominantly classifies areas as moderate (40.59%) and low (38.50%) susceptibility, while the GLM model identifies a larger portion of the area as very low (24.55%) and low (38.59%) susceptibility. Both models exhibit high accuracy metrics, with SVM achieving an accuracy of 87.4% and GLM close behind at 87.2%. Furthermore, the GLM model demonstrates a slightly better Area Under the Curve (AUC) score of 0.939, compared to 0.916 for the SVM model, indicating superior model performance in distinguishing between different classes of SES. The analysis of contributing factors highlights distinct parameter weights, with the GLM model placing greater emphasis on hydrological factors such as distance to rivers and drainage density, while the SVM model gives more balanced attention across multiple environmental variables. These differences suggest that GLM might offer a more nuanced understanding of hydrological influences on soil erosion. The robust performance of these ML models, validated through cross-validation techniques, underscores their potential for dynamic and accurate SES assessments. This is crucial for guiding effective land management and conservation strategies, particularly in regions vulnerable to soil erosion. These findings emphasize the value of adopting innovative, data-driven approaches in environmental modelling to address the complex challenges of land degradation and soil conservation.\u003c/p\u003e"},{"header":"Declarations","content":"\u003ch1\u003eAcknowledgements\u003c/h1\u003e\n\u003cp\u003eWe sincerely thank the Faculty of Engineering at Universitas Gorontalo for their invaluable support and resources, which were crucial to the success of this research. We also appreciate the contributions of our colleagues, research assistants, and collaborators, whose expertise and dedication significantly enhanced the quality of this study.\u003c/p\u003e\n\u003ch1\u003eAuthor Contributions\u003c/h1\u003e\n\u003cp\u003eOlii, M.R. led the conceptualization, methodology design, and overall supervision of the research project, and played a significant role in drafting and revising the manuscript. Nento, S. was responsible for data collection, statistical analysis, and interpretation of results, contributing to the writing of the results and discussion sections. Doda, N. conducted the literature review and developed the theoretical framework, assisting in writing the introduction and literature review sections. Djafar, H. handled fieldwork, data acquisition, and processing, and contributed to the development of environmental variables for the SES model. Olii, R.Z.N. supported software development, model simulation, and validation, and was involved in visualizing the results. Pakaya, P. contributed to writing, formatting the manuscript, and preparing tables and figures for publication. All authors reviewed and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe data supporting this study\u0026apos;s findings are available from the corresponding author upon request, for academic and research purposes, subject to confidentiality agreements and data use policies.\u003c/p\u003e\n\u003ch1\u003eConflicts of Interest\u003c/h1\u003e\n\u003cp\u003eThe authors declare that there are no conflicts of interest regarding the publication of this paper\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAl-Bawi, A. J., Al-Abadi, A. M., Pradhan, B., \u0026amp; Alamri, A. M. (2021). Assessing gully erosion susceptibility using topographic derived attributes, multi-criteria decision-making, and machine learning classifiers. \u003cem\u003eGeomatics, Natural Hazards and Risk\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(1), 3035\u0026ndash;3062. https://doi.org/10.1080/19475705.2021.1994024\u003c/li\u003e\n\u003cli\u003eAllafta, H., \u0026amp; Opp, C. (2022). Soil Erosion Assessment Using the RUSLE Model, Remote Sensing, and GIS in the Shatt Al-Arab Basin (Iraq-Iran). \u003cem\u003eApplied Sciences (Switzerland)\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(15), 1\u0026ndash;17. https://doi.org/10.3390/app12157776\u003c/li\u003e\n\u003cli\u003eAlmouctar, M. A. S., Wu, Y., Zhao, F., \u0026amp; Dossou, J. F. (2021). Soil erosion assessment using the rusle model and geospatial techniques (Remote sensing and gis) in south-central niger (maradi region). \u003cem\u003eWater (Switzerland)\u003c/em\u003e, \u003cem\u003e13\u003c/em\u003e(24). https://doi.org/10.3390/w13243511\u003c/li\u003e\n\u003cli\u003eArabameri, A., Cerda, A., \u0026amp; Tiefenbacher, J. P. (2019). Spatial pattern analysis and prediction of gully erosion using novel hybrid model of entropy-weight of evidence. \u003cem\u003eWater (Switzerland)\u003c/em\u003e, \u003cem\u003e11\u003c/em\u003e(6), 1\u0026ndash;23. https://doi.org/10.3390/w11061129\u003c/li\u003e\n\u003cli\u003eArabameri, A., Pradhan, B., Pourghasemi, H. R., \u0026amp; Rezaei, K. (2018). Identification of erosion-prone areas using different multi-criteria decision-making techniques and gis. \u003cem\u003eGeomatics, Natural Hazards and Risk\u003c/em\u003e, \u003cem\u003e9\u003c/em\u003e(1), 1129\u0026ndash;1155. https://doi.org/10.1080/19475705.2018.1513084\u003c/li\u003e\n\u003cli\u003eBand, S. S., Janizadeh, S., Pal, S. C., Saha, A., Chakrabortty, R., Shokri, M., \u0026amp; Mosavi, A. (2020). Novel ensemble approach of deep learning neural network (Dlnn) model and particle swarm optimization (pso) algorithm for prediction of gully erosion susceptibility. \u003cem\u003eSensors (Switzerland)\u003c/em\u003e, \u003cem\u003e20\u003c/em\u003e(19), 1\u0026ndash;28. https://doi.org/10.3390/s20195609\u003c/li\u003e\n\u003cli\u003eBui, D. T., Tsangaratos, P., Nguyen, V. T., Liem, N. Van, \u0026amp; Trinh, P. T. (2020). Comparing the prediction performance of a Deep Learning Neural Network model with conventional machine learning models in landslide susceptibility assessment. \u003cem\u003eCatena\u003c/em\u003e, \u003cem\u003e188\u003c/em\u003e(December 2019), 104426. https://doi.org/10.1016/j.catena.2019.104426\u003c/li\u003e\n\u003cli\u003eConoscenti, C., Angileri, S., Cappadonia, C., Rotigliano, E., Agnesi, V., \u0026amp; M\u0026auml;rker, M. (2014). Gully erosion susceptibility assessment by means of GIS-based logistic regression: A case of Sicily (Italy). \u003cem\u003eGeomorphology\u003c/em\u003e, \u003cem\u003e204\u003c/em\u003e, 399\u0026ndash;411. https://doi.org/10.1016/j.geomorph.2013.08.021\u003c/li\u003e\n\u003cli\u003eDas, B., Bordoloi, R., Thungon, L. T., Paul, A., Pandey, P. K., Mishra, M., \u0026amp; Tripathi, O. P. (2020). An integrated approach of GIS, RUSLE and AHP to model soil erosion in West Kameng watershed, Arunachal Pradesh. \u003cem\u003eJournal of Earth System Science\u003c/em\u003e, \u003cem\u003e129\u003c/em\u003e(1), 1\u0026ndash;18. https://doi.org/10.1007/s12040-020-1356-6\u003c/li\u003e\n\u003cli\u003eDevos, O., Ruckebusch, C., Durand, A., Duponchel, L., \u0026amp; Huvenne, J. P. (2009). Support vector machines (SVM) in near infrared (NIR) spectroscopy: Focus on parameters optimization and model interpretation. \u003cem\u003eChemometrics and Intelligent Laboratory Systems\u003c/em\u003e, \u003cem\u003e96\u003c/em\u003e(1), 27\u0026ndash;33. https://doi.org/10.1016/j.chemolab.2008.11.005\u003c/li\u003e\n\u003cli\u003eEekhout, J. P. C., \u0026amp; de Vente, J. (2022). Global impact of climate change on soil erosion and potential for adaptation through soil conservation. \u003cem\u003eEarth-Science Reviews\u003c/em\u003e, \u003cem\u003e226\u003c/em\u003e, 103921. https://doi.org/10.1016/j.earscirev.2022.103921\u003c/li\u003e\n\u003cli\u003eEjegu, M. A., \u0026amp; Yegizaw, E. S. (2021). Modeling soil erosion susceptibility and LULC dynamics for land degradation management using geoinformation technology in Debre Tabor district, Northwestern highlands of Ethiopia. \u003cem\u003eJ. Degrade. Min. Land Manage\u003c/em\u003e, \u003cem\u003e8\u003c/em\u003e(2), 2623\u0026ndash;2633. https://doi.org/10.15243/jdmlm.2021.082.2623\u003c/li\u003e\n\u003cli\u003eFarhan, Y., Zregat, D., \u0026amp; Farhan, I. (2013). Spatial Estimation of Soil Erosion Risk Using RUSLE Approach, RS, and GIS Techniques: A Case Study of Kufranja Watershed, Northern Jordan. In \u003cem\u003eJournal of Water Resource and Protection\u003c/em\u003e (Vol. 05, Issue 12, pp. 1247\u0026ndash;1261). https://doi.org/10.4236/jwarp.2013.512134\u003c/li\u003e\n\u003cli\u003eFolharini, S., Vieira, A., Bento-Gon\u0026ccedil;alves, A., Silva, S., Marques, T., \u0026amp; Novais, J. (2023). Soil Erosion Quantification using Machine Learning in Sub-Watersheds of Northern Portugal. \u003cem\u003eHydrology\u003c/em\u003e, \u003cem\u003e10\u003c/em\u003e(1). https://doi.org/10.3390/hydrology10010007\u003c/li\u003e\n\u003cli\u003eFran\u0026ccedil;ois, M., Gon\u0026ccedil;alves Pontes, M. C., de Vasconcelos, R. N., de Oliveira, U. C., Peixoto da Silva, H., Faria, D., \u0026amp; Mariano-Neto, E. (2024). Assessing soil erosion and its drivers in agricultural landscapes: a case study in southern Bahia, Brazil. \u003cem\u003eJournal of Water and Climate Change\u003c/em\u003e, \u003cem\u003e00\u003c/em\u003e(0), 1\u0026ndash;16. https://doi.org/10.2166/wcc.2024.147\u003c/li\u003e\n\u003cli\u003eGayen, A., Pourghasemi, H. R., Saha, S., Keesstra, S., \u0026amp; Bai, S. (2019). Gully erosion susceptibility assessment and management of hazard-prone areas in India using different machine learning algorithms. \u003cem\u003eScience of the Total Environment\u003c/em\u003e, \u003cem\u003e668\u003c/em\u003e, 124\u0026ndash;138. https://doi.org/10.1016/j.scitotenv.2019.02.436\u003c/li\u003e\n\u003cli\u003eGhorbanzadeh, O., Shahabi, H., Mirchooli, F., Valizadeh Kamran, K., Lim, S., Aryal, J., Jarihani, B., \u0026amp; Blaschke, T. (2020). Gully erosion susceptibility mapping (GESM) using machine learning methods optimized by the multi‑collinearity analysis and K-fold cross-validation. \u003cem\u003eGeomatics, Natural Hazards and Risk\u003c/em\u003e, \u003cem\u003e11\u003c/em\u003e(1), 1653\u0026ndash;1678. https://doi.org/10.1080/19475705.2020.1810138\u003c/li\u003e\n\u003cli\u003eGolijanin, J., Nikolić, G., Valjarević, A., Ivanović, R., Tunguz, V., Bojić, S., Grmu\u0026scaron;a, M., Lukić Tanović, M., Perić, M., Hrelja, E., \u0026amp; Stankov, S. (2022). Estimation of potential soil erosion reduction using GIS-based RUSLE under different land cover management models: A case study of Pale Municipality, B\u0026amp;H. \u003cem\u003eFrontiers in Environmental Science\u003c/em\u003e, \u003cem\u003e10\u003c/em\u003e(August), 1\u0026ndash;13. https://doi.org/10.3389/fenvs.2022.945789\u003c/li\u003e\n\u003cli\u003eGolkarian, A., Khosravi, K., Panahi, M., \u0026amp; Clague, J. J. (2023). Spatial variability of soil water erosion: Comparing empirical and intelligent techniques. \u003cem\u003eGeoscience Frontiers\u003c/em\u003e, \u003cem\u003e14\u003c/em\u003e(1), 101456. https://doi.org/10.1016/j.gsf.2022.101456\u003c/li\u003e\n\u003cli\u003eHuang, D., Su, L., Zhou, L., Tian, Y., \u0026amp; Fan, H. (2023). Assessment of gully erosion susceptibility using different DEM-derived topographic factors in the black soil region of Northeast China. \u003cem\u003eInternational Soil and Water Conservation Research\u003c/em\u003e, \u003cem\u003e11\u003c/em\u003e(1), 97\u0026ndash;111. https://doi.org/10.1016/j.iswcr.2022.04.001\u003c/li\u003e\n\u003cli\u003eIgwe, O., John, U. I., Solomon, O., \u0026amp; Obinna, O. (2020). GIS-based gully erosion susceptibility modeling, adapting bivariate statistical method and AHP approach in Gombe town and environs Northeast Nigeria. \u003cem\u003eGeoenvironmental Disasters\u003c/em\u003e, \u003cem\u003e7\u003c/em\u003e(1). https://doi.org/10.1186/s40677-020-00166-8\u003c/li\u003e\n\u003cli\u003eJiang, C., Yang, Z., Wang, X., Dong, X., Li, Z., \u0026amp; Li, C. (2021). Examining the reversal of soil erosion decline in the hotspots of sandstorms: A non-linear ecosystem dynamic perspective. \u003cem\u003eJournal of Arid Environments\u003c/em\u003e, \u003cem\u003e186\u003c/em\u003e(October 2020), 104421. https://doi.org/10.1016/j.jaridenv.2020.104421\u003c/li\u003e\n\u003cli\u003eKhosravi, K., Rezaie, F., Cooper, J. R., Kalantari, Z., Abolfathi, S., \u0026amp; Hatamiafkoueieh, J. (2023). Soil water erosion susceptibility assessment using deep learning algorithms. \u003cem\u003eJournal of Hydrology\u003c/em\u003e, \u003cem\u003e618\u003c/em\u003e. https://doi.org/10.1016/j.jhydrol.2023.129229\u003c/li\u003e\n\u003cli\u003eKucuker, D. M., \u0026amp; Cedano Giraldo, D. (2022). Assessment of soil erosion risk using an integrated approach of GIS and Analytic Hierarchy Process (AHP) in Erzurum, Turkiye. \u003cem\u003eEcological Informatics\u003c/em\u003e, \u003cem\u003e71\u003c/em\u003e, 101788. https://doi.org/https://doi.org/10.1016/j.ecoinf.2022.101788\u003c/li\u003e\n\u003cli\u003eMosavi, A., Sajedi-Hosseini, F., Choubin, B., Taromideh, F., Rahi, G., \u0026amp; Dineva, A. A. (2020). Susceptibility mapping of soil water erosion using machine learning models. \u003cem\u003eWater (Switzerland)\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(7). https://doi.org/10.3390/w12071995\u003c/li\u003e\n\u003cli\u003eMusasa, T., Dube, T., \u0026amp; Marambanyika, T. (2024). Landsat satellite programme potential for soil erosion assessment and monitoring in arid environments: A review of applications and challenges. \u003cem\u003eInternational Soil and Water Conservation Research\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(2), 267\u0026ndash;278. https://doi.org/10.1016/j.iswcr.2023.10.003\u003c/li\u003e\n\u003cli\u003eMustafa, M. R. U., Sholagberu, A. T., Yusof, K. W., Hashim, A. M., Khan, M. W. A., \u0026amp; Shahbaz, M. (2018). SVM-Based Geospatial Prediction of Soil Erosion under Static and Dynamic Conditioning Factors. \u003cem\u003eMATEC Web of Conferences\u003c/em\u003e, \u003cem\u003e203\u003c/em\u003e. https://doi.org/10.1051/matecconf/201820304004\u003c/li\u003e\n\u003cli\u003eNguyen, C. Q., Tran, T. T., Nguyen, T. T. T., Nguyen, T. T., Astarkhanova, T. S., Vu, L. V., Dau, K. T., Nguyen, H. N., Pham, G. H., Nguyen, D. D., Prakash, I., \u0026amp; Pham, B. (2023). Mapping of soil erosion susceptibility using advanced machine learning models at Nghe An, Vietnam. \u003cem\u003eJournal of Hydroinformatics\u003c/em\u003e, \u003cem\u003e26\u003c/em\u003e(1). https://doi.org/10.2166/hydro.2023.327\u003c/li\u003e\n\u003cli\u003eOlii, M. R., Kironoto, B. A., Olii, A., Pakaya, R., \u0026amp; Olii, A. K. Z. (2024). Advancing Soil Erosion Assessment: Application of Remote Sensing and Geospatial Techniques in Bulango Ulu Reservoir Basin. \u003cem\u003eE3S Web of Conferences\u003c/em\u003e, \u003cem\u003e476\u003c/em\u003e, 1\u0026ndash;15. https://doi.org/10.1051/e3sconf/202447601041\u003c/li\u003e\n\u003cli\u003eOlii, M. R., Olii, A. K. Z., Olii, A., Pakaya, R., \u0026amp; Kironoto, B. A. (2024). Spatial modeling of soil erosion risk a multi‑criteria decision‑making (MCDM) approach in the paguyaman watershed, gorontalo, Indo.pdf. \u003cem\u003eArabian Journal of Geosciences\u003c/em\u003e, \u003cem\u003e17\u003c/em\u003e(226), 1\u0026ndash;13. https://doi.org/https://doi.org/10.1007/s12517-024-12032-0\u003c/li\u003e\n\u003cli\u003eOlii, M. R., Olii, A., \u0026amp; Pakaya, R. (2021). The Integrated Spatial Assessment of The Flood Hazard Using AHP-GIS: The Case Study of Gorontalo Regency. \u003cem\u003eIndonesian Journal of Geography\u003c/em\u003e, \u003cem\u003e53\u003c/em\u003e(1), 126\u0026ndash;135. https://doi.org/10.22146/ijg.59999\u003c/li\u003e\n\u003cli\u003eOlii, M. R., Olii, A., Pakaya, R., \u0026amp; Olii, M. Y. U. P. (2023). GIS-based analytic hierarchy process (AHP) for soil erosion-prone areas mapping in the Bone Watershed, Gorontalo, Indonesia. \u003cem\u003eEnvironmental Earth Sciences\u003c/em\u003e, \u003cem\u003e82\u003c/em\u003e(9), 1\u0026ndash;14. https://doi.org/10.1007/s12665-023-10913-3\u003c/li\u003e\n\u003cli\u003ePhinzi, K., \u0026amp; Szab\u0026oacute;, S. (2024). Predictive machine learning for gully susceptibility modeling with geo-environmental covariates: main drivers, model performance, and computational efficiency. In \u003cem\u003eNatural Hazards\u003c/em\u003e (Vol. 120, Issue 8). Springer Netherlands. https://doi.org/10.1007/s11069-024-06481-9\u003c/li\u003e\n\u003cli\u003ePourghasemi, H. R., Sadhasivam, N., Kariminejad, N., \u0026amp; Collins, A. L. (2020). Gully erosion spatial modelling: Role of machine learning algorithms in selection of the best controlling factors and modelling process. \u003cem\u003eGeoscience Frontiers\u003c/em\u003e, \u003cem\u003e11\u003c/em\u003e(6), 2207\u0026ndash;2219. https://doi.org/10.1016/j.gsf.2020.03.005\u003c/li\u003e\n\u003cli\u003eRahmati, O., Tahmasebipour, N., Haghizadeh, A., Pourghasemi, H. R., \u0026amp; Feizizadeh, B. (2017). Evaluating the influence of geo-environmental factors on gully erosion in a semi-arid region of Iran: An integrated framework. \u003cem\u003eScience of the Total Environment\u003c/em\u003e, \u003cem\u003e579\u003c/em\u003e, 913\u0026ndash;927. https://doi.org/10.1016/j.scitotenv.2016.10.176\u003c/li\u003e\n\u003cli\u003eSaini, S. S., Jangra, R., \u0026amp; Kaushik, S. P. (2015). Vulnerability Assessment of Soil Erosion Using Geospatial Techniques-a Pilot Study of Upper Catchment of Markanda River. \u003cem\u003eResearch Article International Journal of Advancement in Remote Sensing\u003c/em\u003e, \u003cem\u003e3\u003c/em\u003e(1), 9\u0026ndash;21.\u003c/li\u003e\n\u003cli\u003eSenanayake, S., Pradhan, B., Huete, A., \u0026amp; Brennan, J. (2020). A review on assessing and mapping soil erosion hazard using geo-informatics technology for farming system management. \u003cem\u003eRemote Sensing\u003c/em\u003e, \u003cem\u003e12\u003c/em\u003e(24), 1\u0026ndash;25. https://doi.org/10.3390/rs12244063\u003c/li\u003e\n\u003cli\u003eVu Dinh, T., Hoang, N. D., \u0026amp; Tran, X. L. (2021). Evaluation of Different Machine Learning Models for Predicting Soil Erosion in Tropical Sloping Lands of Northeast Vietnam. \u003cem\u003eApplied and Environmental Soil Science\u003c/em\u003e, \u003cem\u003e2021\u003c/em\u003e. https://doi.org/10.1155/2021/6665485\u003c/li\u003e\n\u003cli\u003eWang, Y., Zhang, Y., \u0026amp; Chen, H. (2023). Gully erosion susceptibility prediction in Mollisols using machine learning models. \u003cem\u003eJournal of Soil and Water Conservation\u003c/em\u003e, \u003cem\u003e78\u003c/em\u003e(5), 385\u0026ndash;396. https://doi.org/10.2489/jswc.2023.00019\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Soil Erosion Susceptibility (SES), Geospatial Modelling, Machine Learning (ML), Support Vector Machines (SVM), Generalized Linear Models (GLM)","lastPublishedDoi":"10.21203/rs.3.rs-4933265/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4933265/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eSoil erosion creates substantial environmental and economic challenges, especially in areas vulnerable to land degradation. This study investigates the use of machine learning (ML) techniques\u0026mdash;namely Support Vector Machines (SVM) and Generalized Linear Models (GLM)\u0026mdash;for geospatial modeling of soil erosion susceptibility (SES). By leveraging geospatial data and incorporating a range of factors including hydrological, topographical, and environmental variables, the research aims to improve the accuracy and reliability of SES predictions. Results show that the SVM model predominantly identifies areas as having moderate (40.59%) or low (38.50%) susceptibility, whereas the GLM model allocates a higher proportion to very low (24.55%) and low (38.59%) susceptibility. Both models exhibit high performance, with SVM and GLM achieving accuracies of 87.4% and 87.2%, respectively, though GLM slightly surpasses AUC (0.939 vs. 0.916). GLM places greater emphasis on hydrological factors such as distance to rivers and drainage density, while SVM provides a more balanced assessment across various variables. This study demonstrates that ML-based models can significantly enhance SES assessments, offering a more nuanced and accurate approach than traditional methods. The findings highlight the value of adopting innovative, data-driven techniques in environmental modeling and offer practical insights for land management and conservation practices.\u003c/p\u003e","manuscriptTitle":"Transformation of Geospatial Modelling of Soil Erosion Susceptibility Using Machine Learning","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-18 14:33:03","doi":"10.21203/rs.3.rs-4933265/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"54acf0e9-8cf0-4536-bf05-d84eb56e40e2","owner":[],"postedDate":"September 18th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2025-01-15T21:38:10+00:00","versionOfRecord":[],"versionCreatedAt":"2024-09-18 14:33:03","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4933265","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4933265","identity":"rs-4933265","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}