Digital soil mapping for soil types using machine learning approaches at the landscape scale in the arid regions of Iran

Abstract

The aims of this research are (i) to compare random forest (RF), boosted regression tree (BRT), and multinomial logistic regression (MnLR) models to prepare the prediction maps of soil great group and subgroup levels, (ii) determination of the most important environmental covariates influencing the production of digital soil mapping (DSM) in an arid climate, (iii) to evaluate the efficiency of spectra indices extracted from Sentinel-2A digital images and data capability of ALOS-PALSAR radar data, and (iv) investigating the effect of sub-surface genetic horizons in the modeling of different types of soil map classes distribution. The principal component analysis method was employed to select the best set from the pool of environmental covariates (n = 46) such as geomorphometric parameters (GPs), RS indices, and diagnostic soil properties (DSP). The relative importance results indicate that Gypsic (GYP) subsurface horizon, standardized height (StH), slope length (SL), and normalized different vegetation index (NDVI) had an important role in the prediction of soil classes compared to the other selected covariates. DSM methodology was used in this research by incorporating of RF model and representative soil-forming factors that can be used for preparing the maps of soil classes in low-relief areas with a similar soil-landscape relationship. Totally, in this study places a spotlight on the profound impact of sub-surface genetic horizons, shedding light on their pivotal role in accurately modeling soil class distributions. These findings not only advance our comprehension of soil variability in arid regions but also hold immense implications for the burgeoning field of pedometrics.