Search

Rapid and accurate mapping of soil properties in farmlands is crucial for guiding agricultural production and maintaining food security. Traditional methods using spectral features from remote sensing prove valuable for estimating soil properties, but are restricted to short periods of bare soil occurrence within agricultural settings. Addressing the challenge of predicting soil properties under crop cover, this study proposed an improved soil modeling framework that integrates dynamic crop growth information with machine learning techniques. The methodology's robustness was tested on six key soil properties in an agricultural region of China, including soil organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), dissolved organic carbon (DOC), dissolved organic nitrogen (DON), and pH. Four experimental scenarios were established to assess the impact of crop growth information, represented by the normalized difference vegetation index (NDVI) and phenological parameters. Specifically, Scenario I utilized only natural factors (terrain and climate data); Scenario II added phenological parameters based on Scenario I; Scenario III incorporated time-series NDVI based on Scenario I; and Scenario IV combined all variables (traditional natural factors and crop growth information). These were evaluated using three advanced machine learning models: random forest (RF), Cubist, and Extreme Gradient Boosting (XGBoost). Results demonstrated that incorporating phenological parameters and time-series NDVI significantly improved model accuracy, enhancing predictions by up to 36% over models using only natural factors. Moreover, although both are crop growth factors, the contribution of the time-series NDVI variable to model accuracy surpassed that of the phenological variable for most soil properties. Relative importance analysis suggested that the crop growth information, derived from time-series NDVI and phenology data, collectively explained 14–45% of the spatial variation in soil properties. This study highlights the significant benefits of integrating remote sensing-based crop growth factors into soil property inversion under crop-covered conditions, providing valuable insights for digital soil mapping. [ABSTRACT FROM AUTHOR]

Soils play a central role in ecosystem functioning, and thus, mapped soil property information is indispensable to supporting sustainable land management. Digital Soil Mapping (DSM) provides a framework to spatially estimate soil properties. However, broad-scale DSM remains challenging because of non-purposively sampled soil data, large data volumes for processing extensive soil covariates, and high model complexities due to spatially varying soil–landscape relationships. This study presents a three-dimensional DSM framework for Switzerland, targeting the soil properties of clay content (Clay), organic carbon content (SOC), pH value (pH), and potential cation exchange capacity (CECpot). The DSM approach is based on machine learning and a comprehensive exploitation of soil and remote sensing data archives. Quantile Regression Forest was applied to link the soil sample data from a national soil data base with covariates derived from a LiDAR-based elevation model, from climate raster data, and from multispectral raster time series based on satellite imagery. The covariate set comprises spatially multiscale terrain attributes, climate patterns and their temporal variation, temporarily multiscale land use features, and spectral bare soil signatures. Soil data and predictions were evaluated with respect to different landcovers and depth intervals. All reference soil data sets were found to be spatially clustered towards croplands, showing an increasing sample density from lower to upper depth intervals. According to the R2 value derived from independent data, the overall model accuracy amounts to 0.69 for Clay, 0.64 for SOC, 0.76 for pH, and 0.72 for CECpot. Reduced model accuracies were found to be accompanied by soil data sets showing limited sample sizes (e.g., CECpot), uneven statistical distributions (e.g., SOC), and low spatial sample densities (e.g., woodland subsoils). Multiscale terrain covariates were highly influential for all models; climate covariates were particularly important for the Clay model; multiscale land use covariates showed enhanced importance for modeling pH; and bare soil reflectance was a major driver in the SOC and CECpot models. [ABSTRACT FROM AUTHOR]

This paper quantifies the annual dynamics of the shortwave radiation reflected from bare arable land as a result of smoothing previously plowed land located in three different agricultural subregions of the European Union and associated countries. This estimate takes into account the annual variation of the bare arable land area, obtained from Sentinel 2 satellite imagery; the spatial variability of soil units within croplands, obtained from digital soil and land-cover maps; and the laboratory spectral reflectance characteristics of these units, obtained from soil samples stored in the LUCAS soil database. The properties of the soil units, which cover an area of at least 4% of each subregion, were characterized. The highest amounts of shortwave radiation reflected under clear-sky conditions from air-dried, bare arable land surfaces—approximately 850 PJ day−1 and 1.10 EJ day−1 for land shaped by a plow (Pd) and smoothing harrow (Hs), respectively—were found in the summer around 8 August in the western subregion. However, the lowest radiation occurred in the spring on 10 April at 340 PJ day−1 for Pd and 430 PJ day−1 for Hs in the central subregion. The largest and the smallest amounts of this radiation throughout the year—only as a result of smoothing, by Hs, land that was previously treated by Pd—was estimated at 42 EJ for the western and southern subregions and 19 EJ for the central subregion, respectively. [ABSTRACT FROM AUTHOR]

The deployment of remote sensing platforms has facilitated the mapping of soil properties to a great extent. However, the accuracy of these soil property estimates is compromised by the presence of non-soil cover, which introduces interference with the acquired reflectance spectra over pixels. Therefore, current soil property estimation by remote sensing is limited to bare soil pixels, which are identified based on spectral indices of vegetation. Our study proposes a composite mapping approach to extend the soil properties mapping beyond bare soil pixels, associated with an uncertainty map. The proposed approach first classified the pixels based on their bare soil fractional cover by spectral unmixing. Then, a specific regression model was built and applied to each bare soil fractional cover class to estimate clay content. Finally, the clay content maps created for each bare soil fractional cover class were mosaicked to create a composite map of clay content estimations. A bootstrap procedure was used to estimate the standard deviation of clay content predictions per bare soil fractional cover dataset, which represented the uncertainty of estimations. This study used a hyperspectral image acquired by the Airborne Visible/Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) sensor over cultivated fields in South India. The proposed approach provided modest performances in prediction ( R v a l 2 ranging from 0.53 to 0.63) depending on the bare soil fractional cover class and showed a correct spatial pattern, regardless of the bare soil fraction classes. The model's performance was observed to increase with the adoption of higher bare soil fractional cover thresholds. The mapped area ranged from 10.4% for pixels with bare soil fractional cover >0.7 to 52.7% for pixels with bare soil fractional cover >0.3. The approach thus extended the mapped surface by 42.4%, while maintaining acceptable prediction performances. Finally, the proposed approach could be adopted to extend the mapping capability of planned and current hyperspectral satellite missions. [ABSTRACT FROM AUTHOR]

Soil particle size fractions (PSFs) are important properties for understanding the physical and chemical processes in soil systems. Knowledge about the distribution of soil PSFs is critical for sustainable soil management. Although log-ratio transformations have been widely applied to soil PSFs prediction, the statistical distribution of original data and the transformed data given by log-ratio transformations is different, resulting in biased estimates of soil PSFs. Therefore, multivariate random forest (MRF) was utilized for the simultaneous prediction of soil PSFs, as it is able to capture dependencies and internal relations among the three components. Specifically, 243 soil samples collected across the Loess Plateau were used. Meanwhile, Landsat data, terrain attributes, and climatic variables were employed as environmental variables for spatial prediction of soil PSFs. The results depicted that MRF gave satisfactory soil PSF prediction performance, where the R2 values were 0.62, 0.53, and 0.73 for sand, silt, and clay, respectively. Among the environmental variables, nighttime land surface temperature (LST_N) presented the highest importance in predicting soil PSFs in the Loess Plateau, China. Maps of soil PSFs and texture were generated at a 30 m resolution, which can be utilized as alternative data for soil erosion management and ecosystem conservation. [ABSTRACT FROM AUTHOR]

Soil organic carbon (SOC) is generally thought to act as a carbon sink; however, in areas with high spatial heterogeneity, using a single model to estimate the SOC of the whole study area will greatly reduce the simulation accuracy. The earth surface unit division is important to consider in building different models. Here, we divided the research area into different habitat patches using partitioning around a medoids clustering (PAM) algorithm; then, we built an SOC simulation model using machine learning algorithms. The results showed that three habitat patches were created. The simulation accuracy for Habitat Patch 1 (R2 = 0.55; RMSE = 2.89) and Habitat Patch 3 (R2 = 0.47; RMSE = 3.94) using the XGBoost model was higher than that for the whole study area (R2 = 0.44; RMSE = 4.35); although the R2 increased by 25% and 6.8%, the RMSE decreased by 33.6% and 9.4%, and the field sample points significantly declined by 70% and 74%. The R2 of Habitat Patch 2 using the RF model increased by 17.1%, and the RMSE also decreased by 10.5%; however, the sample points significantly declined by 58%. Therefore, using different models for corresponding patches will significantly increase the SOC simulation accuracy over using one model for the whole study area. This will provide scientific guidance for SOC or soil property monitoring with low field survey costs and high simulation accuracy. [ABSTRACT FROM AUTHOR]

Soil organic carbon (SOC) contents and stocks provide valuable insights into soil health, nutrient cycling, greenhouse gas emissions, and overall ecosystem productivity. Given this, remote sensing data coupled with advanced machine learning (ML) techniques have eased SOC level estimation while revealing its patterns across different ecosystems. However, despite these advances, the intricacies of training reliable and yet certain SOC models for specific end-users remain a great challenge. To address this, we need robust SOC uncertainty quantification techniques. Here, we introduce a methodology that leverages conformal prediction to address the uncertainty in estimating SOC contents while using remote sensing data. Conformal prediction generates statistically reliable uncertainty intervals for predictions made by ML models. Our analysis, performed on the LUCAS dataset in Europe and incorporating a suite of relevant environmental covariates, underscores the efficacy of integrating conformal prediction with another ML model, specifically random forest. In addition, we conducted a comparative assessment of our results against prevalent uncertainty quantification methods for SOC prediction, employing different evaluation metrics to assess both model uncertainty and accuracy. Our methodology showcases the utility of the generated prediction sets as informative indicators of uncertainty. These sets accurately identify samples that pose prediction challenges, providing valuable insights for end-users seeking reliable predictions in the complexities of SOC estimation. [ABSTRACT FROM AUTHOR]

In digital soil mapping, machine learning models have been widely applied. However, the accuracy of machine learning models can be limited by the use of a single model and a small number of soil samples. This study introduces a novel method, semi-supervised classification combined with stacking learning (SSC-SL), to enhance soil classification mapping in hilly and low-mountain areas of Northern Jurong City, Jiangsu Province, China. This study incorporated Gaofen-2 (GF-2) remote sensing imagery along with its associated remote sensing indices, the ALOS Digital Elevation Model (DEM) and their derived topographic factors, and soil parent material data in its modelling process. We first used three base learners, Ranger, Rpart, and XGBoost, to construct the SL model. In addition, we employed the fuzzy c-means clustering algorithm (FCM) to construct a clustering map. To fully leverage the information from a multitude of environmental variables, understand the distribution of data, and enhance the effectiveness of the classification, we selected unlabelled samples near the boundaries of the patches on the clustering map. The SSC-SL model demonstrated superior stability and performance, with optimal accuracy at a 0.9 confidence level, achieving an overall accuracy of 0.77 and a kappa coefficient of 0.73. These metrics exceeded those of the highest performing base learner (Ranger model) by 10.4% and 12.3%, respectively, and they outperformed the least effective base learner (Rpart model) by 27.3% and 32.9%. It notably improves the spatial distribution accuracy of soil types. Key environmental variables influencing soil type distribution include soil parent material (SPM), land use (LU), the multi-resolution valley bottom flatness index (MRVBF), and Elevation (Ele). In conclusion, the SSC-SL model offers a novel and effective approach for enhancing the predictive accuracy of soil classification mapping. [ABSTRACT FROM AUTHOR]

Conventional maps of soil parent material (SPM) types obtained by field survey and manual mapping or predictions from other map data have limited accuracy. Digital soil mapping of SPM types necessitates accurate acquisition of SPM distribution information, which is still a challenge in hilly areas. This study developed a high-accuracy method for SPM identification in hilly areas at the county scale. Based on geographic information system technology, seven feature variables were extracted from the geological map, geomorphic map, digital elevation model, and remote sensing image data of Shanggao County, Jiangxi Province, China. Different feature combination schemes were designed to develop SPM identification models based on random forest (RF), support vector machine (SVM), and maximum likelihood classification (MLC) algorithms. The best SPM identification results were obtained from the RF algorithm using the combination of geological type, geomorphic type, elevation, and slope. Confusion matrices were constructed based on a field survey of 586 validation samples, and the results were evaluated in terms of overall accuracy, precision, recall, F1 score, and Kappa coefficient. The overall accuracy and Kappa coefficient of the results from the optimal RF model were 83.11% and 0.79, respectively, which were 26.11% and 0.31 higher than those of the conventional map, respectively. Its precision and recall for various SPM types were greater than 75%. A comprehensive comparison of the accuracy, uncertainty, and plotting performance of the SPM recognition results reveals that the RF algorithm outperforms the SVM algorithm and the MLC algorithm. Geological type was the largest contributor to SPM identification, followed by geomorphic type, elevation, and slope. The importance of different feature variables varied for distinct SPM types. The accuracy of SPM identification was not improved by selecting more feature variables, such as land use type, normalised difference vegetation index, and topographic wetness index. This study demonstrates the feasibility of high-accuracy county-level SPM mapping in hilly areas based on the RF algorithm using geological type, geomorphic type, elevation, and slope as feature variables. As hilly areas have typical topographic features and SPM types, the proposed method of SPM mapping can be useful for application in other similar areas. There are a few limitations in this study with regard to data quality and resolution, feature variable selection, classification algorithm generalisation, and study area representativeness, which may affect the outcomes and need to be solved. [ABSTRACT FROM AUTHOR]

Despite optical remote sensing (and the spectral vegetation indices) contributions to digital soil-mapping studies of soil organic carbon (SOC), few studies have used active radar remote sensing mission data like that from synthetic aperture radar (SAR) sensors to predict SOC. Bearing in mind the importance of SOC mapping for agricultural, ecological, and climate interests and also the recently developed methods for vegetation monitoring using Sentinel-1 SAR data, in this work, we aimed to take advantage of the high operationality of Sentinel-1 imaging to test the accuracy of SOC prediction at different soil depths using machine learning systems. Using linear, nonlinear, and tree regression-based methods, it was possible to predict the SOC content of soils from western Bahia, Brazil, a region with predominantly sandy soils, using as explanatory variables the SAR vegetation indices. The models fed with SAR sensor polarizations and vegetation indices produced more accurate results for the topsoil layers (0–5 cm and 5–10 cm in depth). In these superficial layers, the models achieved an RMSE in the order of 5.0 g kg−1 and an R2 ranging from 0.16 to 0.24, therefore explaining about 20% of SOC variability using only Sentinel-1 predictors. [ABSTRACT FROM AUTHOR]

Soil texture is a key physical property that affects the soil's ability to retain moisture and nutrients. As a result, it is of extreme importance to conduct remote sensing monitoring of soil texture. Songnen Plain is located in the black soil belt of Northeast China. The development of satellite imagery in remote sensing technology enables the rapid monitoring of large areas. This study aimed to map the surface soil texture of cultivated land in Songnen Plain using Sentinel-2 images and Random Forest (RF) algorithm. We conducted this study by collecting 354 topsoil (0–20 cm) samples in Songnen Plain and evaluating the effectiveness of the bands and spectral indices of Sentinel-2 images and RF algorithm in predicting soil texture (sand, silt, and clay fractions). The results demonstrated that the 16 covariates were moderately and highly correlated with soil texture. And, Band11 of Sentinel-2 images could be used as the corresponding band of soil texture. For sand fraction, the Sentinel-2 images and RF algorithm's Coefficient of Determination (R2) and Root Mean Square Error (RMSE) were 0.77 and 10.48%, respectively, and for silt fraction, they were 0.75 and 9.38%. Sand fraction decreased from southwest to northeast in Songnen Plain, while silt and clay fractions increased. We found that the Songnen Plain was affected by water erosion and wind erosion, in the northeast and southwest, respectively, providing reference for the implementation of Conservation Tillage policies. The outcome of the study can provide reference for future soil texture mapping with a high resolution. [ABSTRACT FROM AUTHOR]

There is a growing realization among policymakers that in order to pave the way for the development of evidence-based conservation recommendations for policy, it is essential to improve the capacity for soil-health monitoring by adopting multidimensional and integrated approaches. However, the existing ready-to-use maps are characterized mainly by a coarse spatial resolution (>200 m) and information that is not up to date, making their use insufficient for the EU's policy requirements, such as the common agricultural policy. This work, by utilizing the Soil Data Cube, which is a self-hosted custom tool, provides yearly estimations of soil thematic maps (e.g., exposed soil, soil organic carbon, clay content) covering all the agricultural area in Lithuania. The pipeline exploits various Earth observation data such as a time series of Sentinel-2 satellite imagery (2018–2022), the LUCAS (Land Use/Cover Area Frame Statistical Survey) topsoil database, the European Integrated Administration and Control System (IACS) and artificial intelligence (AI) architectures to improve the prediction accuracy as well as the spatial resolution (10 m), enabling discrimination at the parcel level. Five different prediction models were tested with the convolutional neural network (CNN) model to achieve the best accuracy for both targeted indicators (SOC and clay) related to the R 2 metric (0.51 for SOC and 0.57 for clay). The model predictions supported by the prediction uncertainties based on the PIR formula (average PIR 0.48 for SOC and 0.61 for clay) provide valuable information on the model's interpretation and stability. The model application and the final predictions of the soil indicators were carried out based on national bare-soil-reflectance composite layers, generated by employing a pixel-based composite approach to the overlaid annual bare-soil maps and by using a combination of a series of vegetation indices such as NDVI, NBR2, and SCL. The findings of this work provide new insights for the generation of soil thematic maps on a large scale, leading to more efficient and sustainable soil management, supporting policymakers and the agri-food private sector. [ABSTRACT FROM AUTHOR]

Readily available environmental covariates in current digital soil mapping usually do not indicate the spatial differences between deep soil attributes. This, to a large extent, leads to a decrease in the accuracy of 3D soil mapping with depth, which seriously affects the quality of soil information generated. This study tested the hypothesis that spatialized laboratory soil spectral information can be used as environmental covariates to improve the accuracy of 3D soil attribute mapping and proposed a new type of environmental covariable. In the first step, with soil-forming environmental covariates and independent soil profiles, laboratory vis-NIR spectral data of soil samples resampled into six bands in Anhui province, China, were spatially interpolated to generate spatial distributions of soil spectral measurements at multiple depths. In the second step, we constructed three sets of covariates using the laboratory soil spectral distribution maps at multiple depths: conventional soil-forming variables (C), conventional soil-forming variables plus satellite remote sensing wavebands (C+SRS) and conventional soil-forming variables plus spatialized laboratory soil spectral information (C+LSS). In the third step, we used the three sets of environmental covariates to develop random forest models for predicting soil attributes (pH; CEC, cation exchange capacity; Silt; SOC, soil organic carbon; TP, total phosphorus) at multiple depths. We compared the 3D soil mapping accuracies between these three sets of covariates based on another dataset of 132 soil profiles (collected in the 1980s). The results show that the use of spatialized laboratory soil spectral information as additional environmental covariates has a 50% improvement in prediction accuracy compared with that of only conventional covariates, and a 30% improvement in prediction accuracy compared with that of the satellite remote sensing wavebands as additional covariates. This indicates that spatialized laboratory soil spectral information can improve the accuracy of 3D digital soil mapping. [ABSTRACT FROM AUTHOR]

Unsustainable human management has negative effects on cropland soil organic carbon (SOC), causing a decrease in soil health and the emission of greenhouse gas. Due to contiguous fields, large-scale mechanized operations are widely used in the Northeast China Plain, which greatly improves production efficiency while decreasing the soil quality, especially for SOC. Therefore, an up-to-date SOC map is needed to estimate soil health after long-term cultivation to inform better land management. Using Quantile Regression Forest, a total of 396 soil samples from 132 sampling sites at three soil depth intervals and 40 environmental covariates (e.g., Landsat 8 spectral indices, and WorldClim 2 and MODIS products) selected by the Boruta feature selection algorithm were used to map the spatial distribution of SOC in the cropland of the Northeast Plain at a 90 m spatial resolution. The results showed that SOC increased overall from the southern area to the northern area, with an average of 17.34 g kg−1 in the plough layer (PL) and 13.92 g kg−1 in the compacted layer (CL). At the vertical scale, SOC decreased, with depths getting deeper. The average decrease in SOC from PL to CL was 3.41 g kg−1. Climate (i.e., average temperature, daytime and nighttime land surface temperature, and mean temperature of driest quarter) was the dominant controlling factor, followed by position (i.e., oblique geographic coordinate at 105°), and organism (i.e., the average and variance of net primary productivity in the non-crop period). The average uncertainty was 1.04 in the PL and 1.07 in the CL. The high uncertainty appeared in the area with relatively scattered fields, high altitudes, and complex landforms. This study updated the 90 m resolution cropland SOC maps at spatial and vertical scales, which clarifies the influence of mechanized operations and provides a reference for soil conservation policy-making. [ABSTRACT FROM AUTHOR]

Accurately mapping the spatial distribution and variation of soil organic matter (SOM) is of great significance for guiding regional soil management. However, the applicability and prediction performance of machine learning techniques in dryland agroecosystems still needs to be further studied. In this study, we collected a total of 733 topsoil samples from the farmland in Xiao County, Anhui Province, which is a typical dryland agroecosystem in the Huang-Huai-Hai Plain. Then, the environmental covariates were selected, and the ordinary kriging (OK), multiple linear stepwise regression (MLR), regression kriging (RK), radial basis function neural network (RBFNN), and random forest (RF) models were conducted to map the SOM content, and the optimal model was ascertained. The results demonstrated that the alkali-hydrolyzable nitrogen (26.11%), available potassium (17.73%), mean annual precipitation (13.26%), and pH (11.80%) were the main controlling factors affecting the spatial distribution of SOM in the study area. Meanwhile, the introduction of environmental covariates can effectively improve the SOM prediction accuracy, and the RF model (R2 = 0.48, MAE = 2.38 g kg−1, MRE = 12.99%, RMSE = 3.14 g kg−1) has a better performance than the RFBNN, MLR, RK, and OK methods. Although there are local differences in the spatial distribution of SOM predicted by the five methods, the overall spatial distribution of SOM was characterized by the low concentration area (13.44–20.00 g kg−1) distributed in the central and northwest of study area, and the high concentration area (24.00–28.95 g kg−1) distributed in the southeast. Overall, our study demonstrated that machine learning-based models could accurately predict the SOM content in dryland agroecosystem, and the produced maps function as baseline maps for sustainable agricultural management. [ABSTRACT FROM AUTHOR]

For most of the arable land in Russia (132–137 million ha), the dominant and accurate soil information is stored in the form of map archives on paper without coordinate reference. The last traditional soil map(s) (TSM, TSMs) were created over 30 years ago. Traditional and/or archival soil map(s) (ASM, ASMs) are outdated in terms of storage formats, dates, and methods of production. The technology of constructing a multitemporal soil line (MSL) makes it possible to update ASMs and TSMs based on the processing of big remote-sensing data (RSD). To construct an MSL, the spectral characteristics of the bare soil surface (BSS) are used. The BSS on RSD is distinguished within the framework of the conceptual apparatus of the spectral neighborhood of the soil line. The filtering of big RSD is based on deep machine learning. In the course of the work, a vector georeferenced version of the ASM and an updated soil map were created based on the coefficient "C" of the MSL. The maps were verified based on field surveys (76 soil pits). The updated map is called the map of soil interpretation of the coefficient "C" (SIC "C"). The SIC "C" map has a more detailed legend compared to the ASM (7 sections/chapters instead of 5), greater accuracy (smaller errors of the first and second kind), and potential suitability for calculating soil organic matter/carbon (SOM/SOC) reserves (soil types/areals in the SIC "C" map are statistically significant are divided according to the thickness of the organomineral horizon and the content of SOM in the plowed layer). When updating, a systematic underestimation of the numbers of contours and areas of soils with manifestations of negative/degradation soil processes (slitization and erosion) on the TSM was established. In the process of updating, all three shortcomings of the ASMs/TSMs (archaic storage, dates, and methods of creation) were eliminated. The SIC "C" map is digital (thematic raster), modern, and created based on big data processing methods. For the first time, the actualization of the soil map was carried out based on the MSL characteristics (coefficient "C"). [ABSTRACT FROM AUTHOR]

The North China agro–pastoral zone is a large, ecologically fragile zone in the arid and semi-arid regions. Quantitative remote sensing inversion of soil organic carbon (SOC) in this region can facilitate understanding of the current status of degraded land restoration and provide data support for carbon cycling research in the region. Deep learning (DNN) for SOC inversion has been W.a hot topic over the past decade, but there have been few studies at the regional scale in the arid and semi-arid zones. In this study, a DNN model with five hidden layers and five skip connections was established using 644 spatially distributed SOC samples and Landsat 8 OLI imagery. The model was compared with the random forest algorithm in terms of generalization ability. The main conclusions were as follows: 1. The DNN algorithm can establish a high-precision SOC inversion model (R2 = 0.52, RMSE = 0.7), with 90% of errors concentrated in the range of −2.5 to 2.5 kg·C/m2; 2. the Boruta variable-screening algorithm can effectively improve the model accuracy of the random forest algorithm, but due to the DNN's better ability to mine hidden information in the data, the improvement effect on the DNN model accuracy is limited; 3. the SOC samples in arid and semi-arid areas are highly positively skewed, with a significant impact on the modeling accuracy of DNN, and conversion is required to obtain a model with better generalization ability; and 4. in arid and semi-arid regions, SOC has a weak correlation with vegetation indices but a stronger correlation with temperature, elevation, and aridity. This study established a reliable deep learning model for SOC density in a large arid and semi-arid region, providing a reference and framework for the establishment of SOC inversion models in other regions. [ABSTRACT FROM AUTHOR]

Airborne geophysical data (AGD) have great potential to represent soil-forming factors. Because of that, the objective of this study was to evaluate the importance of AGD in predicting soil attributes such as aluminum saturation (ASat), base saturation (BS), cation exchange capacity (CEC), clay, and organic carbon (OC). The AGD predictor variables include total count (μR/h), K (potassium), eU (uranium equivalent), and eTh (thorium equivalent), ratios between these elements (eTh/K, eU/K, and eU/eTh), factor F or F-parameter, anomalous potassium (Kd), anomalous uranium (Ud), anomalous magnetic field (AMF), vertical derivative (GZ), horizontal derivatives (GX and GY), and mafic index (MI). The approach was based on applying predictive modeling techniques using (1) digital elevation model (DEM) covariates and Sentinel-2 images with AGD; and (2) DEM covariates and Sentinel-2 images without the AGD. The study was conducted in Bom Jardim, a county in Rio de Janeiro-Brazil with an area of 382,430 km², with a database of 208 soil samples to a predefined depth (0–30 cm). Non-explanatory covariates for the selected soil attributes were excluded. Through the selected covariables, the random forest (RF) and support vector machine (SVM) models were applied with separate samples for training (75%) and validation (25%). The model's performance was evaluated through the R-squared (R2), root mean square error (RMSE), and mean absolute error (MAE), as well as null model values and coefficient of variation (CV%). The RF algorithm showed better performance with AGD (R2 values ranging from 0.15 to 0.23), as well as the SVM model (R2 values ranging from 0.08 to 0.23) when compared to RF (R2 values ranging from 0.10 to 0.20) and SVM (R2 values ranging from 0.04 to 0.10) models without AGD. Overall, the results suggest that AGD can be helpful for soil mapping. Nevertheless, it is crucial to acknowledge that the accuracy of AGD in predicting soil properties could vary depending on various common factors in DSM, such as the quality and resolution of the covariates and available soil data. Further research is needed to determine the optimal approach for using AGD in soil mapping. [ABSTRACT FROM AUTHOR]

The deployment of remote sensing platforms has facilitated the mapping of soil properties to a great extent. However, the accuracy of these soil property estimates is compromised by the presence of non-soil cover, which introduces interference with the acquired reflectance spectra over pixels. Therefore, current soil property estimation by remote sensing is limited to bare soil pixels, which are identified based on spectral indices of vegetation. Our study proposes a composite mapping approach to extend the soil properties mapping beyond bare soil pixels, associated with an uncertainty map. The proposed approach first classified the pixels based on their bare soil fractional cover by spectral unmixing. Then, a specific regression model was built and applied to each bare soil fractional cover class to estimate clay content. Finally, the clay content maps created for each bare soil fractional cover class were mosaicked to create a composite map of clay content estimations. A bootstrap procedure was used to estimate the standard deviation of clay content predictions per bare soil fractional cover dataset, which represented the uncertainty of estimations. This study used a hyperspectral image acquired by the Airborne Visible/Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) sensor over cultivated fields in South India. The proposed approach provided modest performances in prediction (Rval2 ranging from 0.53 to 0.63) depending on the bare soil fractional cover class and showed a correct spatial pattern, regardless of the bare soil fraction classes. The model’s performance was observed to increase with the adoption of higher bare soil fractional cover thresholds. The mapped area ranged from 10.4% for pixels with bare soil fractional cover >0.7 to 52.7% for pixels with bare soil fractional cover >0.3. The approach thus extended the mapped surface by 42.4%, while maintaining acceptable prediction performances. Finally, the proposed approach could be adopted to extend the mapping capability of planned and current hyperspectral satellite missions.

Soil particle size fractions (PSFs) are important properties for understanding the physical and chemical processes in soil systems. Knowledge about the distribution of soil PSFs is critical for sustainable soil management. Although log-ratio transformations have been widely applied to soil PSFs prediction, the statistical distribution of original data and the transformed data given by log-ratio transformations is different, resulting in biased estimates of soil PSFs. Therefore, multivariate random forest (MRF) was utilized for the simultaneous prediction of soil PSFs, as it is able to capture dependencies and internal relations among the three components. Specifically, 243 soil samples collected across the Loess Plateau were used. Meanwhile, Landsat data, terrain attributes, and climatic variables were employed as environmental variables for spatial prediction of soil PSFs. The results depicted that MRF gave satisfactory soil PSF prediction performance, where the R2 values were 0.62, 0.53, and 0.73 for sand, silt, and clay, respectively. Among the environmental variables, nighttime land surface temperature (LST_N) presented the highest importance in predicting soil PSFs in the Loess Plateau, China. Maps of soil PSFs and texture were generated at a 30 m resolution, which can be utilized as alternative data for soil erosion management and ecosystem conservation.

Soil organic carbon (SOC) is generally thought to act as a carbon sink; however, in areas with high spatial heterogeneity, using a single model to estimate the SOC of the whole study area will greatly reduce the simulation accuracy. The earth surface unit division is important to consider in building different models. Here, we divided the research area into different habitat patches using partitioning around a medoids clustering (PAM) algorithm; then, we built an SOC simulation model using machine learning algorithms. The results showed that three habitat patches were created. The simulation accuracy for Habitat Patch 1 (R2 = 0.55; RMSE = 2.89) and Habitat Patch 3 (R2 = 0.47; RMSE = 3.94) using the XGBoost model was higher than that for the whole study area (R2 = 0.44; RMSE = 4.35); although the R2 increased by 25% and 6.8%, the RMSE decreased by 33.6% and 9.4%, and the field sample points significantly declined by 70% and 74%. The R2 of Habitat Patch 2 using the RF model increased by 17.1%, and the RMSE also decreased by 10.5%; however, the sample points significantly declined by 58%. Therefore, using different models for corresponding patches will significantly increase the SOC simulation accuracy over using one model for the whole study area. This will provide scientific guidance for SOC or soil property monitoring with low field survey costs and high simulation accuracy.

Soil organic carbon (SOC) contents and stocks provide valuable insights into soil health, nutrient cycling, greenhouse gas emissions, and overall ecosystem productivity. Given this, remote sensing data coupled with advanced machine learning (ML) techniques have eased SOC level estimation while revealing its patterns across different ecosystems. However, despite these advances, the intricacies of training reliable and yet certain SOC models for specific end-users remain a great challenge. To address this, we need robust SOC uncertainty quantification techniques. Here, we introduce a methodology that leverages conformal prediction to address the uncertainty in estimating SOC contents while using remote sensing data. Conformal prediction generates statistically reliable uncertainty intervals for predictions made by ML models. Our analysis, performed on the LUCAS dataset in Europe and incorporating a suite of relevant environmental covariates, underscores the efficacy of integrating conformal prediction with another ML model, specifically random forest. In addition, we conducted a comparative assessment of our results against prevalent uncertainty quantification methods for SOC prediction, employing different evaluation metrics to assess both model uncertainty and accuracy. Our methodology showcases the utility of the generated prediction sets as informative indicators of uncertainty. These sets accurately identify samples that pose prediction challenges, providing valuable insights for end-users seeking reliable predictions in the complexities of SOC estimation.

In digital soil mapping, machine learning models have been widely applied. However, the accuracy of machine learning models can be limited by the use of a single model and a small number of soil samples. This study introduces a novel method, semi-supervised classification combined with stacking learning (SSC-SL), to enhance soil classification mapping in hilly and low-mountain areas of Northern Jurong City, Jiangsu Province, China. This study incorporated Gaofen-2 (GF-2) remote sensing imagery along with its associated remote sensing indices, the ALOS Digital Elevation Model (DEM) and their derived topographic factors, and soil parent material data in its modelling process. We first used three base learners, Ranger, Rpart, and XGBoost, to construct the SL model. In addition, we employed the fuzzy c-means clustering algorithm (FCM) to construct a clustering map. To fully leverage the information from a multitude of environmental variables, understand the distribution of data, and enhance the effectiveness of the classification, we selected unlabelled samples near the boundaries of the patches on the clustering map. The SSC-SL model demonstrated superior stability and performance, with optimal accuracy at a 0.9 confidence level, achieving an overall accuracy of 0.77 and a kappa coefficient of 0.73. These metrics exceeded those of the highest performing base learner (Ranger model) by 10.4% and 12.3%, respectively, and they outperformed the least effective base learner (Rpart model) by 27.3% and 32.9%. It notably improves the spatial distribution accuracy of soil types. Key environmental variables influencing soil type distribution include soil parent material (SPM), land use (LU), the multi-resolution valley bottom flatness index (MRVBF), and Elevation (Ele). In conclusion, the SSC-SL model offers a novel and effective approach for enhancing the predictive accuracy of soil classification mapping.

The spatial heterogeneity of soil properties has a significant impact on crop growth, making it difficult to adopt site-specific crop management practices. Traditional laboratory-based analyses are costly, and data extrapolation for mapping soil properties using high-resolution imagery becomes a computationally expensive procedure, taking days or weeks to obtain accurate results using a desktop workstation. To overcome these challenges, cloud-based solutions such as Google Earth Engine (GEE) have been used to analyze complex data with machine learning algorithms. In this study, we explored the feasibility of designing and implementing a digital soil mapping approach in the GEE platform using high-resolution reflectance imagery derived from a thermal infrared and multispectral camera Altum (MicaSense, Seattle, WA, USA). We compared a suite of multispectral-derived soil and vegetation indices with in situ measurements of physical-chemical soil properties in agricultural lands in the Peruvian Mantaro Valley. The prediction ability of several machine learning algorithms (CART, XGBoost, and Random Forest) was evaluated using R2, to select the best predicted maps (R2 > 0.80), for ten soil properties, including Lime, Clay, Sand, N, P, K, OM, Al, EC, and pH, using multispectral imagery and derived products such as spectral indices and a digital surface model (DSM). Our results indicate that the predictions based on spectral indices, most notably, SRI, GNDWI, NDWI, and ExG, in combination with CART and RF algorithms are superior to those based on individual spectral bands. Additionally, the DSM improves the model prediction accuracy, especially for K and Al. We demonstrate that high-resolution multispectral imagery processed in the GEE platform has the potential to develop soil properties prediction models essential in establishing adaptive soil monitoring programs for agricultural regions. [ABSTRACT FROM AUTHOR]

This article aims to explore the use of machine learning (ML) methods for mapping the distribution of mercury (Hg) content in topsoil, using the city of Ufa (Russia) and adjacent areas as an example. For this purpose, a soil dataset of 250 points sampled from a 0–20 cm depth on different land uses, including residential, industrial and undisturbed (forests and parks), was used. Random Forest (RF), Extreme Gradient Boosting (XGboost), Cubist and k-Nearest Neighbor (kNN) ML techniques were employed to model and map the Hg concentrations. We used remote sensing data (RSD) and topographic attributes as explanatory variables. ML models were calibrated and validated using the leave-one-out cross-validation approach. The Hg content varied from 0.005 to 0.58 mg/kg and was characterized by very high variability. According to the MAE and RMSE metrics, the RF method resulted in the most accurate spatial prediction for the Hg content (0.029 and 0.065 mg/kg, respectively), while the XGBoost approach showed the lowest prediction efficiency (0.032 and 0.073 mg/kg, respectively). The results showed that the slope map, spectral index MSI and Sentinel-2A band B11 were the key variables in explaining the variability of Hg content. We found that higher uncertainty values of soil Hg were found in croplands, urban residential and industrial areas, which supports the view that spatial modelling of HM in urban landscapes is challenging. The present study provides insights into the potential of digital soil mapping techniques in combination with RSD and terrain variables for identifying areas at risk of Hg contamination in urban areas, which can inform land-use planning and management strategies to protect human health and the environment. [ABSTRACT FROM AUTHOR]

Soils are at the crossroads of many existential issues that humanity is currently facing. Soils are a finite resource that is under threat, mainly due to human pressure. There is an urgent need to map and monitor them at field, regional, and global scales in order to improve their management and prevent their degradation. This remains a challenge due to the high and often complex spatial variability inherent to soils. Over the last four decades, major research efforts in the field of pedometrics have led to the development of methods allowing to capture the complex nature of soils. As a result, digital soil mapping (DSM) approaches have been developed for quantifying soils in space and time. DSM and monitoring have become operational thanks to the harmonization of soil databases, advances in spatial modeling and machine learning, and the increasing availability of spatiotemporal covariates, including the exponential increase in freely available remote sensing (RS) data. The latter boosted research in DSM, allowing the mapping of soils at high resolution and assessing the changes through time. We present a review of the main contributions and developments of French (inter)national research, which has a long history in both RS and DSM. Thanks to the French SPOT satellite constellation that started in the early 1980s, the French RS and soil research communities have pioneered DSM using remote sensing. This review describes the data, tools, and methods using RS imagery to support the spatial predictions of a wide range of soil properties and discusses their pros and cons. The review demonstrates that RS data are frequently used in soil mapping (i) by considering them as a substitute for analytical measurements, or (ii) by considering them as covariates related to the controlling factors of soil formation and evolution. It further highlights the great potential of RS imagery to improve DSM, and provides an overview of the main challenges and prospects related to digital soil mapping and future sensors. This opens up broad prospects for the use of RS for DSM and natural resource monitoring. [ABSTRACT FROM AUTHOR]

Soil total nitrogen (STN) is a crucial component of the ecosystem's nitrogen pool, and accurate prediction of STN content is essential for understanding global nitrogen cycling processes. This study utilized the measured STN content of 126 sample points and 40 extracted remote sensing variables to predict the STN content and map its spatial distribution in the northeastern coastal region of Hebei Province, China, employing the random forest (RF), gradient boosting machine (GBM), and extreme gradient boosting (XGBoost) methods. The purpose was to compare the ability of remote sensing images (Landsat-8, Sentinel-1, and Sentinel-2) with different machine learning methods for predicting STN content. The research results show the following: (1) The three machine learning methods accurately predicted the STN content and the optimal model provided by the XGBoost method, with an R2 of 0.627, RMSE of 0.127 g·kg−1, and MAE of 0.092 g·kg−1. (2) The combination of optical and synthetic aperture radar (SAR) images improved prediction accuracy, with the R2 improving by 45.5%. (3) The importance of optical images is higher than that of SAR images in the RF, GBM, and XGBoost methods, with optical images accounting for 87%, 76%, and 77% importance, respectively. (4) The spatial distribution of STN content predicted by the three methods is similar. Higher STN contents are distributed in the northern part of the study area, while lower STN contents are distributed in coastal areas. The results of this study can be very useful for inventories of soil nitrogen and provide data support and method references for revealing nitrogen cycling. [ABSTRACT FROM AUTHOR]

Satellite-based soil organic carbon content (SOC) mapping over wide regions is generally hampered by the low soil sampling density and the diversity of soil sampling periods. Some unfavorable topsoil conditions, such as high moisture, rugosity, the presence of crop residues, the limited amplitude of SOC values and the limited area of bare soil when a single image is used, are also among the influencing factors. To generate a reliable SOC map, this study addresses the use of Sentinel-2 (S2) temporal mosaics of bare soil (S2Bsoil) over 6 years jointly with soil moisture products (SMPs) derived from Sentinel 1 and 2 images, SOC measurement data and other environmental covariates derived from digital elevation models, lithology maps and airborne gamma-ray data. In this study, we explore (i) the dates and periods that are preferable to construct temporal mosaics of bare soils while accounting for soil moisture and soil management; (ii) which set of covariates is more relevant to explain the SOC variability. From four sets of covariates, the best contributing set was selected, and the median SOC content along with uncertainty at 90% prediction intervals were mapped at a 25-m resolution from quantile regression forest models. The accuracy of predictions was assessed by 10-fold cross-validation, repeated five times. The models using all the covariates had the best model performance. Airborne gamma-ray thorium, slope and S2 bands (e.g., bands 6, 7, 8, 8a) and indices (e.g., calcareous sedimentary rocks, "calcl") from the "late winter–spring" time series were the most important covariates in this model. Our results also indicated the important role of neighboring topographic distances and oblique geographic coordinates between remote sensing data and parent material. These data contributed not only to optimizing SOC mapping performance but also provided information related to long-range gradients of SOC spatial variability, which makes sense from a pedological point of view. [ABSTRACT FROM AUTHOR]

To overcome spatial, spectral and temporal constraints of different remote sensing products, data fusion is a good technique to improve the prediction capability of soil prediction models. However, few studies have analyzed the effects of image fusion on digital soil mapping (DSM) models. This research fused multispectral (MS) and panchromatic Landsat 8 (L8) bands, and MS Sentinel 2 (S2) and panchromatic L8 bands using the Brovey, Intensity–Hue–Saturation and Gram–Schmidt methods in an agricultural area in Yellow River Basin, China. To analyze the effects of image fusion on DSM models, various SOC prediction models derived from remote sensing image datasets were established by the random forest method. Soil salinity indices and spectral reflectance from all the remote sensing data had relatively strong negative correlations with SOC, and vegetation indices and water indices from all the remote sensing data had relatively strong positive correlations with SOC. Soil moisture and vegetation were the main controlling factors of the SOC spatial pattern in the study area. More spectral indices derived from pansharpened L8 and fused S2–L8 images by all three image fusion methods had stronger relationships with SOC compared with those from MS L8 and MS S2, respectively. All the SOC models established by pansharpened L8 and fused S2–L8 images had higher prediction accuracy than those established by MS L8 and MS S2, respectively. The fusion between S2 and L8 bands had stronger effects on enhancing the prediction accuracy of SOC models compared with the fusion between panchromatic and MS L8 bands. It is concluded that digital soil mapping and image fusion can be utilized to increase the prediction performance of SOC spatial prediction models. [ABSTRACT FROM AUTHOR]

Thermo-erosion gullies (TGs) are typical thermokarst features in upland permafrost; the soil organic carbon (SOC) of TGs has an important influence on soil quality in cold regions. The objectives of this study were to estimate the spatial distribution of SOC content in a typical TG on the northeastern Tibetan Plateau in China by using soil properties from seven different TGs and covariates from unmanned aerial vehicle (UAV) images, and to characterize the SOC content changes in four representative landscape regions (NO-Slumping, Slumping1, Slumping2, and Slumped) within this typical TG. The support vector machine (SVM) was the optimal machine learning algorithm for SOC content prediction, which explained 53.06% (R2) of the SOC content variation. Silt content was the most influential factor which demonstrated a positive relationship with SOC content in different TGs. In addition, the SOC content in the TGs was related to the landscapes. Severe Slumping (Slumping2: 150.79 g·kg−1) had a lower SOC content than NO-Slumped (163.29 g·kg−1) and the initial slumping stage (Slumping1: 169.08 g·kg−1). The results suggested that SVM was an effective algorithm to obtain a profound understanding of the SOC content over space, while future research needs to pay more attention to the SOC content distribution in the different TGs. [ABSTRACT FROM AUTHOR]

The rapid acquisition of high-resolution spatial distribution of soil organic matter (SOM) at the field scale is essential for precision agriculture. The UAV imaging hyperspectral technology, with its high spatial resolution and timeliness, can fill the research gap between ground-based monitoring and remote sensing. This study aimed to test the feasibility of using UAV hyperspectral data (400–1000 nm) with a small-sized calibration sample set for mapping SOM at a 1 m resolution in typical low-relief black soil areas of Northeast China. The experiment was conducted in an approximately 20 ha field. For calibration, 20 samples were collected using a 100 × 100 m grid sampling strategy, while 20 samples were randomly collected for independent validation. UAV captured hyperspectral images with a spatial resolution of 0.05 × 0.05 m. The extracted spectra within every 1 × 1 m were then averaged to represent the spectra of that grid; this procedure was also performed across the whole field. Upon applying various spectral pretreatments, including absorbance conversion, multiple scattering correction, Savitzky–Golay smoothing filtering, and first-order differentiation, the absolute maximum values of the correlation coefficients of the spectra for SOM increased from 0.41 to 0.58. Importance analysis from the optimal random forest (RF) model showed that the characterized bands of SOM were located in the 450–600 and 750–900 nm regions. When the RF model was used, the UAV hyperspectra data (UAV-RF) were able to successfully predict SOM, with an R2 of 0.53 and RMSE of 1.48 g kg−1. The prediction accuracy was then compared with that obtained using ordinary kriging (OK) and the RF model based on proximal sensing (PS-RF) with the same number of calibration samples. However, the OK method failed to predict the SOM accuracy (RMSE = 2.17 g kg−1; R2 = 0.02) due to a low sampling density. The semi-covariance function was unable to describe the spatial variability of SOM effectively. When the sampling density was increased to 50 × 50 m, OK successfully predicted SOM, with RMSE = 1.37 g kg−1 and R2 = 0.59, and its results were comparable to those of UAV-RF. The prediction accuracy of PS-RF was generally consistent with that of UAV-RF, with RMSE values of 1.41 g kg−1 and 1.48 g kg−1 and R2 values of 0.57 and 0.53, respectively, which indicated that SOM prediction based on UAV-RF is feasible. Additionally, compared with the PS platforms, the UAV hyperspectral technology could simultaneously provide spectral information of tens or even hundreds of continuous bands and spatial information at the same time. This study provides a reference for further research and development of UAV hyperspectral techniques for fine-scale SOM mapping using a small number of samples. [ABSTRACT FROM AUTHOR]

A methodology is presented for the quantitative assessment of soil biomass productivity at 100 m spatial resolution on a national scale. The traditional land evaluation approach—where crop yield is the dependent variable—was followed using measured yield and net primary productivity data derived from satellite images, together with digital soil and climate maps. In addition to characterizing of soil biomass productivity based on measured data, the weight of soil properties on productivity was also quantified to provide measured soil health and soil quality indicators as an information base for designing sustainable land management practices. To produce these results, we used only the Random Forest method for our calculations. The study considers high-input agriculture, which is predominant in the country. Biomass productivity indices for the main crops (wheat, maize and sunflowers) and general productivity indices were calculated for the whole agricultural area of Hungary. Results can be implemented in cadastral systems, in applied in agricultural and rural development programs. The assessment can be repeated for monitoring purposes to support general monitoring objectives as well as for reporting in relation to the United Nations Sustainable Development Goals. However, on the basis of the results, we also propose a method for periodically updating the assessment, which can also be used for monitoring biomass productivity in the context of climate change, land degradation and the development of cultivation technology. [ABSTRACT FROM AUTHOR]

Soil organic matter (SOM), as the greatest carbon storage in the terrestrial environment, is inextricably related to the global carbon cycle and global climate change. Accurate estimation and mapping of SOM content are crucial for guiding agricultural output and management, as well as controlling the climate issue. Traditional chemical analysis is unable to satisfy the dynamic estimation of SOM due to its low timeliness. Remote and proximal sensing have significant advantages in terms of ease of use, estimation accuracy, and geographical resolution. In this study, we developed a framework based on machine learning to estimate SOM with high accuracy and resolution using Fourier mid-infrared attenuation total reflectance spectroscopy (FTIR-ATR), Sentinel-2 images, and DEM derivatives. This framework's performance was evaluated on a regional scale using 245 soil samples from northeast China. Results indicated that the calibration size could be shrunk to 50% while achieving a fair prediction performance for SOM content. The Lasso, partial least squares (PLS), support vector regression (SVR), and convolutional neural networks (CNN) performed well in predicting SOM from FTIR-ATR spectra, and the performance was enhanced further by using Sentinel-2 images and DEM derivates. The PLS, SVR, and CNN models created SOM maps with higher spatial resolution and variation than the Kriging approach. The PLS and SVR models provided enough variety and were more realistic in the local SOM map, making them usable at the field scale, and the suggested framework took a fresh look at high-resolution SOM mapping. [ABSTRACT FROM AUTHOR]

Precise knowledge about the soil organic carbon (SOC) content in cropland soils is one requirement to design and execute effective climate and food policies. In digital soil mapping (DSM), machine learning algorithms are used to predict soil properties from covariates derived from traditional soil mapping, digital elevation models, land use, and Earth observation (EO). However, such DSM models are trained for a specific dataset and region and have so far only allowed limited general statements to be made that would enable the models to be transferred to different regions. In this study, we test the transferability of SOC models for cropland soils using five different covariate groups: multispectral soil reflectance composites (satellite), soil legacy data (soil), digital elevation model derivatives (terrain), climate parameters (climate), and combined models (combined). The transferability was analyzed using data from two federal states in southern Germany: Bavaria and Baden-Wuerttemberg. First, baseline models were trained for each state with combined models performing best in both cases (R2 = 0.68/0.48). Next, the models were transferred and tested with soil samples from the other state whose data were not used during model calibration. Only satellite and combined models were transferable, but accuracy declined in both cases. In the final step, models were trained with samples from both states (mixed-data models) and applied to each state separately. This process significantly improved the accuracies of satellite, terrain, and combined models, while it showed no effect on climate models and decreased the models based on soil covariates. The experiment underlines the importance of EO for the transfer and extrapolation of DSM models. [ABSTRACT FROM AUTHOR]

Accurate estimation of root-zone soil moisture (SM) is of great significance for accurate irrigation management. This study was purposed to identify planted-by-planted mapping of root-zone SM on three critical fruit growth periods based on UAV multispectral images using three machine learning (ML) algorithms in a kiwifruit orchard in Shaanxi, China. Several spectral variables were selected based on variable importance (VIP) rankings, including reflectance Ri at wavelengths 560, 668, 740, and 842 nm. Results indicated that the VIP method effectively reduced 42 vegetation indexes (VIs) to less than 7 with an evaluation accuracy of root-zone SM models. Compared with deep root-zone SM models (SM40 and SM60), shallow root-zone SM models (SM10, SM20, and SM30) have better performance (R2 from 0.65 to 0.82, RRMSE from 0.02 to 0.03, MAE from 0.20 to 0.54) in the three fruit growth stages. Among three ML algorithms, random forest models were recommended for simulating kiwi root-zone SM during the critical fruit growth period. Overall, the proposed planted-by-planted root-zone SM estimation approach can be considered a great tool to upgrade the toolbox of the growers in site-specific field management for the high spatiotemporal resolution of SM maps. [ABSTRACT FROM AUTHOR]

Soil organic carbon (SOC) is a critical indicator for the global carbon cycle and the overall carbon pool balance. Obtaining soil maps of surface SOC is fundamental to evaluating soil quality, regulating climate change, and global carbon cycle modeling. However, efficient approaches for obtaining accurate SOC information remain challenging, especially in remote or inaccessible regions of the Qinghai–Tibet Plateau (QTP), which is influenced by complex terrains, climate change, and human activities. This study employed field measurements, SoilGrids250m (SOC_250m, a spatial resolution of 250 m × 250 m), and Sentinel-2 images with different machine learning methods to map SOC content in the QTP. Four machine learning methods including partial least squares regression (PLSR), support vector machines (SVM), random forest (RF), and artificial neural network (ANN) were used to construct spatial prediction models based on 396 field-collected sampling points and various covariates from remote sensing images. Our results revealed that the RF model outperformed the PLSR, SVM, and ANN models, with a higher determination coefficient (R2 of 0.82 is from the training datasets) and the ratio of performance to deviation (RPD = 2.54). The selected covariates according to the variable importance in projection (VIP) were: SOC_250m, B2, B11, Soil-Adjusted Vegetation Index (SAVI), Normalized Difference Vegetation Index (NDVI), B5, and Soil-Adjusted Total Vegetation Index (SATVI). The predicted SOC map showed an overall decrease in SOC content ranging from 69.30 g·kg−1 in the southeast to 1.47 g·kg−1 in the northwest. Our prediction showed spatial heterogeneity of SOC content, indicating that Sentinel-2 images were acceptable for characterizing the variability of SOC. The findings provide a scientific basis for carbon neutrality in the QTP and a reference for the digital mapping of SOC in the alpine region. [ABSTRACT FROM AUTHOR]

Conventional maps of soil parent material (SPM) types obtained by field survey and manual mapping or predictions from other map data have limited accuracy. Digital soil mapping of SPM types necessitates accurate acquisition of SPM distribution information, which is still a challenge in hilly areas. This study developed a high-accuracy method for SPM identification in hilly areas at the county scale. Based on geographic information system technology, seven feature variables were extracted from the geological map, geomorphic map, digital elevation model, and remote sensing image data of Shanggao County, Jiangxi Province, China. Different feature combination schemes were designed to develop SPM identification models based on random forest (RF), support vector machine (SVM), and maximum likelihood classification (MLC) algorithms. The best SPM identification results were obtained from the RF algorithm using the combination of geological type, geomorphic type, elevation, and slope. Confusion matrices were constructed based on a field survey of 586 validation samples, and the results were evaluated in terms of overall accuracy, precision, recall, F1 score, and Kappa coefficient. The overall accuracy and Kappa coefficient of the results from the optimal RF model were 83.11% and 0.79, respectively, which were 26.11% and 0.31 higher than those of the conventional map, respectively. Its precision and recall for various SPM types were greater than 75%. A comprehensive comparison of the accuracy, uncertainty, and plotting performance of the SPM recognition results reveals that the RF algorithm outperforms the SVM algorithm and the MLC algorithm. Geological type was the largest contributor to SPM identification, followed by geomorphic type, elevation, and slope. The importance of different feature variables varied for distinct SPM types. The accuracy of SPM identification was not improved by selecting more feature variables, such as land use type, normalised difference vegetation index, and topographic wetness index. This study demonstrates the feasibility of high-accuracy county-level SPM mapping in hilly areas based on the RF algorithm using geological type, geomorphic type, elevation, and slope as feature variables. As hilly areas have typical topographic features and SPM types, the proposed method of SPM mapping can be useful for application in other similar areas. There are a few limitations in this study with regard to data quality and resolution, feature variable selection, classification algorithm generalisation, and study area representativeness, which may affect the outcomes and need to be solved.

Despite optical remote sensing (and the spectral vegetation indices) contributions to digital soil-mapping studies of soil organic carbon (SOC), few studies have used active radar remote sensing mission data like that from synthetic aperture radar (SAR) sensors to predict SOC. Bearing in mind the importance of SOC mapping for agricultural, ecological, and climate interests and also the recently developed methods for vegetation monitoring using Sentinel-1 SAR data, in this work, we aimed to take advantage of the high operationality of Sentinel-1 imaging to test the accuracy of SOC prediction at different soil depths using machine learning systems. Using linear, nonlinear, and tree regression-based methods, it was possible to predict the SOC content of soils from western Bahia, Brazil, a region with predominantly sandy soils, using as explanatory variables the SAR vegetation indices. The models fed with SAR sensor polarizations and vegetation indices produced more accurate results for the topsoil layers (0–5 cm and 5–10 cm in depth). In these superficial layers, the models achieved an RMSE in the order of 5.0 g kg−1 and an R2 ranging from 0.16 to 0.24, therefore explaining about 20% of SOC variability using only Sentinel-1 predictors.

There is a growing realization among policymakers that in order to pave the way for the development of evidence-based conservation recommendations for policy, it is essential to improve the capacity for soil-health monitoring by adopting multidimensional and integrated approaches. However, the existing ready-to-use maps are characterized mainly by a coarse spatial resolution (>200 m) and information that is not up to date, making their use insufficient for the EU’s policy requirements, such as the common agricultural policy. This work, by utilizing the Soil Data Cube, which is a self-hosted custom tool, provides yearly estimations of soil thematic maps (e.g., exposed soil, soil organic carbon, clay content) covering all the agricultural area in Lithuania. The pipeline exploits various Earth observation data such as a time series of Sentinel-2 satellite imagery (2018–2022), the LUCAS (Land Use/Cover Area Frame Statistical Survey) topsoil database, the European Integrated Administration and Control System (IACS) and artificial intelligence (AI) architectures to improve the prediction accuracy as well as the spatial resolution (10 m), enabling discrimination at the parcel level. Five different prediction models were tested with the convolutional neural network (CNN) model to achieve the best accuracy for both targeted indicators (SOC and clay) related to the R2 metric (0.51 for SOC and 0.57 for clay). The model predictions supported by the prediction uncertainties based on the PIR formula (average PIR 0.48 for SOC and 0.61 for clay) provide valuable information on the model’s interpretation and stability. The model application and the final predictions of the soil indicators were carried out based on national bare-soil-reflectance composite layers, generated by employing a pixel-based composite approach to the overlaid annual bare-soil maps and by using a combination of a series of vegetation indices such as NDVI, NBR2, and SCL. The findings of this work provide new insights for the generation of soil thematic maps on a large scale, leading to more efficient and sustainable soil management, supporting policymakers and the agri-food private sector.

Soil texture is a key physical property that affects the soil’s ability to retain moisture and nutrients. As a result, it is of extreme importance to conduct remote sensing monitoring of soil texture. Songnen Plain is located in the black soil belt of Northeast China. The development of satellite imagery in remote sensing technology enables the rapid monitoring of large areas. This study aimed to map the surface soil texture of cultivated land in Songnen Plain using Sentinel-2 images and Random Forest (RF) algorithm. We conducted this study by collecting 354 topsoil (0–20 cm) samples in Songnen Plain and evaluating the effectiveness of the bands and spectral indices of Sentinel-2 images and RF algorithm in predicting soil texture (sand, silt, and clay fractions). The results demonstrated that the 16 covariates were moderately and highly correlated with soil texture. And, Band11 of Sentinel-2 images could be used as the corresponding band of soil texture. For sand fraction, the Sentinel-2 images and RF algorithm’s Coefficient of Determination (R2) and Root Mean Square Error (RMSE) were 0.77 and 10.48%, respectively, and for silt fraction, they were 0.75 and 9.38%. Sand fraction decreased from southwest to northeast in Songnen Plain, while silt and clay fractions increased. We found that the Songnen Plain was affected by water erosion and wind erosion, in the northeast and southwest, respectively, providing reference for the implementation of Conservation Tillage policies. The outcome of the study can provide reference for future soil texture mapping with a high resolution.

Accurately mapping the spatial distribution and variation of soil organic matter (SOM) is of great significance for guiding regional soil management. However, the applicability and prediction performance of machine learning techniques in dryland agroecosystems still needs to be further studied. In this study, we collected a total of 733 topsoil samples from the farmland in Xiao County, Anhui Province, which is a typical dryland agroecosystem in the Huang-Huai-Hai Plain. Then, the environmental covariates were selected, and the ordinary kriging (OK), multiple linear stepwise regression (MLR), regression kriging (RK), radial basis function neural network (RBFNN), and random forest (RF) models were conducted to map the SOM content, and the optimal model was ascertained. The results demonstrated that the alkali-hydrolyzable nitrogen (26.11%), available potassium (17.73%), mean annual precipitation (13.26%), and pH (11.80%) were the main controlling factors affecting the spatial distribution of SOM in the study area. Meanwhile, the introduction of environmental covariates can effectively improve the SOM prediction accuracy, and the RF model (R2 = 0.48, MAE = 2.38 g kg−1, MRE = 12.99%, RMSE = 3.14 g kg−1) has a better performance than the RFBNN, MLR, RK, and OK methods. Although there are local differences in the spatial distribution of SOM predicted by the five methods, the overall spatial distribution of SOM was characterized by the low concentration area (13.44–20.00 g kg−1) distributed in the central and northwest of study area, and the high concentration area (24.00–28.95 g kg−1) distributed in the southeast. Overall, our study demonstrated that machine learning-based models could accurately predict the SOM content in dryland agroecosystem, and the produced maps function as baseline maps for sustainable agricultural management.

Unsustainable human management has negative effects on cropland soil organic carbon (SOC), causing a decrease in soil health and the emission of greenhouse gas. Due to contiguous fields, large-scale mechanized operations are widely used in the Northeast China Plain, which greatly improves production efficiency while decreasing the soil quality, especially for SOC. Therefore, an up-to-date SOC map is needed to estimate soil health after long-term cultivation to inform better land management. Using Quantile Regression Forest, a total of 396 soil samples from 132 sampling sites at three soil depth intervals and 40 environmental covariates (e.g., Landsat 8 spectral indices, and WorldClim 2 and MODIS products) selected by the Boruta feature selection algorithm were used to map the spatial distribution of SOC in the cropland of the Northeast Plain at a 90 m spatial resolution. The results showed that SOC increased overall from the southern area to the northern area, with an average of 17.34 g kg−1 in the plough layer (PL) and 13.92 g kg−1 in the compacted layer (CL). At the vertical scale, SOC decreased, with depths getting deeper. The average decrease in SOC from PL to CL was 3.41 g kg−1. Climate (i.e., average temperature, daytime and nighttime land surface temperature, and mean temperature of driest quarter) was the dominant controlling factor, followed by position (i.e., oblique geographic coordinate at 105°), and organism (i.e., the average and variance of net primary productivity in the non-crop period). The average uncertainty was 1.04 in the PL and 1.07 in the CL. The high uncertainty appeared in the area with relatively scattered fields, high altitudes, and complex landforms. This study updated the 90 m resolution cropland SOC maps at spatial and vertical scales, which clarifies the influence of mechanized operations and provides a reference for soil conservation policy-making.

Remote sensing soil properties in a coherent manner is now feasible from regional to global scales. Xu et al. [[10]] showed that soil-probe contact measurements resulted in the optimum spectral quality and estimation accuracy; that sieving the soil sample into particle sizes below 1 mm and drying before spectral measurement effectively enhanced spectral quality and estimation accuracy; and that the variation in soil temperature did not have a distinct influence on spectral quality. As the change in soil organic carbon (SOC) over time is an important indicator of soil health and for CO SB 2 sb sequestration, the determination and mapping of SOC content from optical Earth observation sensors is currently an area of rigorous research. [Extracted from the article]

Soil texture is an important property that controls the mobility of the water and nutrients in soil. This study examined the capability of machine learning (ML) models in estimating soil texture fractions using different combinations of remotely sensed data from Sentinel-1 (S1), Sentinel-2 (S2), and terrain-derived covariates (TDC) across two contrasting agroecological regions in Southwest Germany, Kraichgau and the Swabian Alb. Importantly, we tested the predictive power of three different ML models: the random forest (RF), the support vector machine (SVM), and extreme gradient boosting (XGB) coupled with the remote sensing data covariates. As expected, ML model performance was not consistent regarding the input covariates, soil texture fractions, and study regions. For example, in the Swabian Alb, the SVM model performed the best for the sand content with S2 + TDC (RMSE = 3.63%, R2 = 0.42), and XGB best predicted the clay content with S1 + S2 + TDC (RMSE = 6.84%, R2 = 0.64). In Kraichgau, the best models for sand (RMSE = 7.54%, R2 = 0.79) and clay contents (RMSE = 6.14%, R2 = 0.48) were obtained using XGB and SVM, respectively. Moreover, the results indicated that TDC were critical in estimating soil texture fractions, especially in Kraichgau, which indicated that topography plays an important role in defining the spatial distribution of soil properties. In contrast, the contribution of remote sensing data better predicted the silt and clay content in the Swabian Alb. The transferability of a region-specific model to the other region was low as indicated by poor predictive performance. The resulting soil-texture-fraction maps could be a significant source of information for efficient land resource management and environmental monitoring. Nonetheless, further research to evaluate the added value of the Sentinel imagery and to better analyze the spatial transferability of machine learning models is highly recommended. [ABSTRACT FROM AUTHOR]

There is increasing demand for more detailed soil maps to support fine-scale land use planning, soil carbon management, and precision agriculture in Saskatchewan. Predictive soil mapping that incorporates a combination of environmental covariates provides a cost-effective tool for generating finer resolution soil maps. This study focused on mapping soil properties for multiple soil horizons in Saskatchewan using historical legacy soil data in combination with remote sensing band indices, bare soil composite imagery, climate data, and terrain attributes. Mapped soil properties included soil organic carbon content (SOC), total nitrogen, cation exchange capacity (CEC), electrical conductivity (EC), inorganic carbon (IOC), sand and clay content, and total profile soil organic carbon stocks. For each of these soil properties, a recursive feature elimination was undertaken to reduce the number of features in the overall model. This process involved iteratively removing features such that random forest out-of-bag error was minimized. Final random forest models were built for each property and evaluated using an independent test dataset. Overall, predictive models were successful for SOC (R2 = 0.71), total nitrogen (R2 = 0.65), CEC (R2 = 0.46), sand content (R2 = 0.44) and clay content (R2 = 0.55). The methods used in this study enable mapping of a greater geographic region of Saskatchewan compared to those previously established that relied solely on bare soil composite imagery. [ABSTRACT FROM AUTHOR]

Soil texture has a great influence on the physical–hydric and chemical behavior of soils. In the Amazon regions, due to the presence of dense forest cover and limited access to roads, carrying out surveys and mapping of soils is challenging. When data exist, they are relatively sparse and the distribution is quite uneven. In this context, machine learning algorithms (ML) associated with remote sensor covariates offer a framework to derive digital maps of soil attributes. The objective of this study was to produce maps of surface and subsurface soil clay, silt, and sand contents in a 13.440 km2 area in the Amazon. The specific objectives were to (a) evaluate the gain in prediction accuracy when using the P-band of airborne radar as a covariate; (b) evaluate two sampling approaches (Reference Area—RA and Total Area—TA); and (c) evaluate the transferability and performance of three ML algorithms: regression tree (RT), random forest (RF), and support vector machine (SVM). The study site was divided into three blocks, called Urucu, Araracanga, and Juruá, respectively. The soil dataset consisted of 151 surface and subsurface sand, silt, and clay observations and 21 covariates (20 relief variables and the backscattering coefficient from the P-band). Both the RA and TA sampling approach used 114 observations for training the prediction models (75%) and 37 for validation (25%). The RA approach was better for the development of sand and silt models. Overall, RF derived the most accurate predictions for all variables. The effect of introducing the P-band backscattering coefficient improved the sand prediction accuracy at the surface and subsurface in Araracanga, which had the highest sand content, with relative improvements (RI) of the R2, root mean square error (RMSE), and mean absolute error (MAE) of 46%, 3%, and 4% at the surface, respectively, and 66.7%, 4.4%, and 5.2% at the subsurface, respectively. For silt, the P-band improved the predictions at the surface in Araracanga, which had the lowest silt contents among the blocks. For clay, adding the P-band improved the RF predictions at the subsurface, with RI of the R2, RMSE, and MAE of 29%, 5%, and 5%, respectively. Despite the low observation density, inherently hindered by the low accessibility of the area and high costs of sampling thereof, the results showed the potential of ML algorithms boosted by airborne radar P-band to map soil clay, silt, and sand contents in the Amazon. [ABSTRACT FROM AUTHOR]

Soil texture is a key soil property driving physical, chemical, biological, and hydrological processes in soils. The rapid development of remote sensing techniques shows great potential for mapping soil properties. This study highlights the effectiveness of multitemporal remote sensing data in identifying soil textural class by using retrieved vegetation properties as proxies of soil properties. The impacts of sensors, modeling resolutions, and modeling techniques on the accuracy of soil texture classification were explored. Multitemporal Landsat-8 and Sentinel-2 images were individually acquired at the same time periods. Three satellite-based experiments with different inputs, i.e., Landsat-8 data, Sentinel-2 data (excluding red-edge parameters), and Sentinel-2 data (including red-edge parameters) were conducted. Modeling was carried out at three spatial resolutions (10, 30, 60 m) using five machine-learning (ML) methods: random forest, support vector machine, gradient-boosting decision tree, categorical boosting, and super learner that combined the four former classifiers based on the stacking concept. In addition, a novel SHapley Addictive Explanation (SHAP) technique was introduced to explain the outputs of the ML model. The results showed that the sensors, modeling resolutions, and modeling techniques significantly affected the prediction accuracy. The models using Sentinel-2 data with red-edge parameters performed consistently best. The models usually gave better results at fine (10 m) and medium (30 m) modeling resolutions than at a coarse (60 m) resolution. The super learner provided higher accuracies than other modeling techniques and gave the highest values of overall accuracy (0.8429), kappa (0.7611), precision (0.8378), recall rate (0.8393), and F1-score (0.8398) at 30 m with Sentinel-2 data involving red-edge parameters. The SHAP technique quantified the contribution of each variable for different soil textural classes, revealing the critical roles of red-edge parameters in separating loamy soils. This study provides comprehensive insights into the effective modeling of soil properties on various scales using multitemporal optical images. [ABSTRACT FROM AUTHOR]

The geosciences suffer from a lack of large georeferenced datasets that can be used to assess and monitor the role of soil organic carbon (SOC) in plant growth, soil fertility, and CO2 sequestration. Publicly available, large field-scale georeferenced datasets are often limited in number and design to serve these purposes. This study provides the first publicly accessible dataset of georeferenced topsoil SOC measurements (n = 840) over a 26-hectare (ha) agricultural field located in southern Ontario, Canada, with a sampling density of ~32 points per ha. As SOC is usually influenced by site topography (i.e., slope and landscape position), each point of the database is associated with a wide range of remote sensing topographic derivatives; as well as with normalized difference vegetation index (NDVI) based value. The NDVI data were extracted from remote sensing Sentinel-2 imagery from over a five-year period (2017–2021). In this paper, the methodology for topsoil sampling, SOC measurement in the lab, as well as producing the suite of topographic derivatives is described. We discuss the opportunities that the database offers in terms of spatially explicit and continuous soil information to support international efforts in digital soil mapping (i.e., SoilGrids250m) as well as other potential applications detailed in the discussion section. We believe that the database with very dense point location measurements can help in conducting carbon stocks and sequestration studies. Such information can be used to help bridge the gap between ground data and remotely sensed datasets or data-derived products from modeling approaches intended to evaluate field-scale rates of agricultural carbon accumulation. The generated topsoil database in this study is archived and publicly available on the Zenodo open-access repository. [ABSTRACT FROM AUTHOR]

Rapid and accurate mapping of soil organic carbon (SOC) is of great significance to understanding the spatial patterns of soil fertility and conducting soil carbon cycle research. Previous studies have dedicated considerable efforts to the spatial prediction of SOC content, but few have systematically quantified the effects of environmental covariates selection, the spatial scales and the model types on SOC prediction accuracy. Here, we spatially predicted SOC content through digital soil mapping (DSM) based on 186 topsoil (0–20 cm) samples in a typical hilly red soil region of southern China. Specifically, we first determined an optimal covariate set from different combinations of multiple environmental variables, including multi-sensor remote sensing images (Sentinel-1 and Sentinel-2), climate variables and DEM derivatives. Furthermore, we evaluated the impacts of spatial resolution (10 m, 30 m, 90 m, 250 m and 1000 m) of covariates and the model types (three linear and three non-linear machine learning techniques) on the SOC prediction. The results of the performance analysis showed that a combination of Sentinel-1/2-derived variables, climate and topographic predictors generated the best predictive performance. Among all variables, remote sensing covariates, especially Sentinel-2-derived predictors, were identified as the most important explanatory variables controlling the variability of SOC content. Moreover, the prediction accuracy declined significantly with the increased spatial scales and achieved the highest using the XGBoost model at 10 m resolution. Notably, non-linear machine learners yielded superior predictive capability in contrast with linear models in predicting SOC. Overall, our findings revealed that the optimal combination of predictor variables, spatial resolution and modeling techniques could considerably improve the prediction accuracy of the SOC content. Particularly, freely accessible Sentinel series satellites showed great potential in high-resolution digital mapping of soil properties. [ABSTRACT FROM AUTHOR]

For most of the arable land in Russia (132–137 million ha), the dominant and accurate soil information is stored in the form of map archives on paper without coordinate reference. The last traditional soil map(s) (TSM, TSMs) were created over 30 years ago. Traditional and/or archival soil map(s) (ASM, ASMs) are outdated in terms of storage formats, dates, and methods of production. The technology of constructing a multitemporal soil line (MSL) makes it possible to update ASMs and TSMs based on the processing of big remote-sensing data (RSD). To construct an MSL, the spectral characteristics of the bare soil surface (BSS) are used. The BSS on RSD is distinguished within the framework of the conceptual apparatus of the spectral neighborhood of the soil line. The filtering of big RSD is based on deep machine learning. In the course of the work, a vector georeferenced version of the ASM and an updated soil map were created based on the coefficient “C” of the MSL. The maps were verified based on field surveys (76 soil pits). The updated map is called the map of soil interpretation of the coefficient “C” (SIC “C”). The SIC “C” map has a more detailed legend compared to the ASM (7 sections/chapters instead of 5), greater accuracy (smaller errors of the first and second kind), and potential suitability for calculating soil organic matter/carbon (SOM/SOC) reserves (soil types/areals in the SIC “C” map are statistically significant are divided according to the thickness of the organomineral horizon and the content of SOM in the plowed layer). When updating, a systematic underestimation of the numbers of contours and areas of soils with manifestations of negative/degradation soil processes (slitization and erosion) on the TSM was established. In the process of updating, all three shortcomings of the ASMs/TSMs (archaic storage, dates, and methods of creation) were eliminated. The SIC “C” map is digital (thematic raster), modern, and created based on big data processing methods. For the first time, the actualization of the soil map was carried out based on the MSL characteristics (coefficient “C”).