Multi-property digital soil mapping at 30-m spatial resolution down to 1 m using extreme gradient boosting tree model and environmental covariates

Information on the variability of soil properties including clay, calcium carbonate equivalent (CCE) and acidity (pH) in space and depth play major role in sustainable agricultural development but scarce in many regions in the world, specifically in developing countries. Digital soil mapping technique has shown strong success in the prediction of soil properties both in space and depth from their relationship with environmental covariates. The aim of the study was to digitally map soil properties including clay, CCE and pH at multiple depths and identify the dominant factors controlling their variability using Minadoab region of Iran as a case study area due to its significance for agricultural production and economic development in the country. A total of 104 soil profiles were randomly selected and described in field during 15–30 July 2017, and 386 samples were collected from each genetic horizon with variable depths. A continuous spline function was fitted to the horizon data and properties were estimated at five standardized depths following GlobalSoilMap project. The soil layers were then grouped into three classes: i) H1: top (0–5 cm); ii) H2–H4: middle (5–15, 15–30 and 30–60 cm); and iii) H5: bottom (60–100 cm) for easy interpretation. A total of 17 static environmental covariates were derived from terrain-related attributes. Additionally, matching with the soil sampling date, a total of eight band ratios in addition to six individual bands were derived from Landsat-8 OLI (Product ID: LC08_L1TP_168034_20170721_20170728_01_T1) for the study area. Finally, all covariates were categorized in nine groups representing various factors controlling soil variations. The digital maps associated with uncertainties were prepared using extreme gradient boosting (XGBoost) Tree model and bootstrapping method, respectively. Overall, the highest accuracy using XGBoost Tree model was found for CCE (R2 = 0.73, on average), followed by pH (R2 = 0.72, on average) and clay (R2 = 0.70, on average). The highest amount of clay was observed for the bottom layer, followed by middle and top layers across the study area. H1 to H3 layers were subjected to tillage operations over the time and may have contributed to the redistribution of clays. A similar trend was also observed for CCE and pH as the location of the study area is underlain by massive limestone. The uncertainties helped quantify over- and under-estimation of the predictions and the factors with the highest variable importance was identified as the dominant factor. The information generated for the study area provides soil variability information that is critical for developing management strategies for sustainable production in the area.