Soil erodibility mapping using remote sensing and in situ soil data with random forest model in a mountainous...

Land degradation is accelerating in the Himalayan ecosystem, resulting in the loss of soil nutrients due to severe erosion. Soil erosion presents a significant environmental challenge, resulting in both on-site and off-site consequences, such as reduced soil productivity and siltation in reservoirs. Soil erodibility (K factor), an inherent soil property, determines the susceptibility of soils to erosion. Sampling across hilly and mountainous terrain pose challenges due to its complex landscape. Despite these challenges, it is essential to study K factor variations in different land use/land cover types to comprehend the threat of erosion. Digital soil mapping offers an opportunity to overcome this limitation by providing spatial predictions of soil properties. The objective of our study is to map the spatial distribution of soil erodibility using the Random Forest (RF) model, a machine learning method based on sampled in situ soil data and environmental covariates. We collected 556 surface soil samples from the mountainous catchment (Tehri dam catchment) using the stratified random sampling approach. The model performed satisfactorily in both training (r² = 0.91; RMSE = 0.00185) and testing (r² = 0.45; RMSE = 0.00318) phases. Subsequently, we generated a digital map with a resolution of 12.5 m to depict the distribution of the K factor. Our analysis revealed that key environmental variables influencing the prediction of the K factor included geology, mean NDVI, and climatic factors. The average K factor value was estimated at 0.0304 and ranging from 0.0251 to 0.0400 t ha h ha^-1 MJ^-1 mm^-1. A higher K factor was observed in the barren land (0.0344) primarily located in the higher and trans-Himalayan region of seasonally snow-covered areas. These areas typically feature young soils with weak soil formation and unstable soil aggregates. Subsequently cropland/cultivated soils (0.0307) exhibited higher K factor values due to the breakdown of soil aggregates by ploughing activities and exposing carbon to decomposition. The average K factor value of evergreen (0.0294) and deciduous (0.0295) forests were the lowest compared to other land use/land cover types indicating the role of forests in resisting soil erosion. By assessing and predicting soil erodibility, land planners and farmers can implement erosion control measures to protect soil health, prevent sedimentation in water bodies, and sustain agricultural productivity in the Himalayas.