Search

The scientific field of precision agriculture employs increasingly innovative techniques to optimize inputs, maximize profitability, and reduce environmental impact. However, obtaining a high number of soil samples is challenging in order to make precision agriculture viable. There is a trade-off between the amount of data needed and the time and resources spent to obtain these data compared to the accuracy of the maps produced with more or fewer points. In the present study, the research was based on an exhaustive dataset of apparent electrical conductivity (aEC) containing 3906 points distributed along 26 transects with spacing between each of up to 40 m, measured by the proximal soil sensor EM38-MK2, for a grain-producing area of 72 ha in São Paulo, Brazil. A second sparse dataset was simulated, showing only four transects with a 400 m distance and, in the end, only 162 aEC points. The aEC map via ordinary kriging (OK) from the grid with 26 transects was considered the reference, and two other mapping approaches were used to map aEC via sparse grid: kriging with external drift (KED) and geographically weighted regression (GWR). These last two methods allow the increment of auxiliary variables, such as those obtained by remote sensors that present spatial resolution compatible with the pivot scale, such as data from the Landsat-8, Aster, and Sentinel-2 satellites, as well as ten terrain covariates derived from the Alos Palsar digital elevation model. The KED method, when used with the sparse dataset, showed a relatively good fit to the aEC data (R2 = 0.78), with moderate prediction accuracy (MAE = 1.26, RMSE = 1.62) and reasonable predictability (RPD = 1.76), outperforming the GWR method, which had the weakest performance (R2 = 0.57, MAE = 1.78, RMSE = 2.30, RPD = 0.81). The reference aEC map using the exhaustive dataset and OK showed the highest accuracy with an R2 of 0.97, no systematic bias (ME = 0), and excellent precision (RMSE = 0.56, RPD = 5.86). Management zones (MZs) derived from these maps were validated using soil texture data from clay samples measured at 0–10 cm depth in a grid of 72 points. The KED method demonstrated the highest potential for accurately defining MZs for irrigation, producing a map that closely resembled the reference MZ map, thereby providing reliable guidance for irrigation management. [ABSTRACT FROM AUTHOR]

Peatlands, which only cover 3 to 5 percent of the global land area, can store up to twice the amount of carbon as the world's forests. Although recognised for their significant role in the global carbon cycle, discovering the global extent of peatlands and their carbon stock remains challenging. Referring to the UNEP's global peatland map, here we present PEATGRIDS, a data product containing global maps of peat thickness and carbon stock created created using the digital soil mapping approach. We compiled over 25,000 observations of peatland thickness, bulk density (BD) and carbon content (CC), globally. Using the Random Forest (RF) algorithm, we estimated peat thickness and peat BD and CC at ~1 km resolution at multiple depths (0–2 m) globally. The estimates were generated using 19 land surface covariates from digital maps and remote sensing images of land use, soil characteristics, topographical features, and climate parameters. The RF models for peat thickness were trained on 25,200 points grouped into six geographic regions. Validation of the peat thickness estimates showed a good performance, with the coefficient of determination (R2)ranging from 0.15 to 0.72. The prediction for peat BD and CC followed the same model architecture and were trained on 17,000 and 7,000 points, respectively. Overall, BD and CC models performed well and consistently across soil layers with average R2 values of 0.61 for BD and 0.48 for CC. Based on the estimated peat thickness, BD and CC, the carbon stock of global peatland was estimated to be 1,029 Pg C for peat dominated area of 6.57 million km2. PEATGRIDS is made available at https://doi.org/10.5281/zenodo.12559239 (Widyastuti et al., 2024) to support further analyses and modelling of peatlands across the globe. [ABSTRACT FROM AUTHOR]

Accurate and high-resolution spatial soil information is crucial for efficient and sustainable land use, management, and conservation. Since the establishment of digital soil mapping (DSM) and the GlobalSoilMap working group, significant advances have been made in spatial soil information globally. However, accurately predicting soil variation over large and complex areas with limited samples remains a challenge, especially for China, which has diverse soil landscapes. To address this challenge, we utilized 11,209 representative multi-source legacy soil profiles (including the Second National Soil Survey of China, World Soil Information Service, First National Soil Survey of China, and regional databases) and high-resolution soil-forming environment characterization. Using advanced Quantile Regression Forest algorithms and a high-performance parallel computing strategy, we developed comprehensive maps of 23 soil physical, chemical and fertility properties at six standard depth layers from 0 to 2 meters in China with a 90 m spatial resolution (China dataset of soil properties for land surface modeling version 2, CSDLv2). Data-splitting and independent samples validation strategies were employed to evaluate the accuracy of the predicted maps quality. The results showed that the predicted maps were significantly more accurate and detailed compared to traditional soil type linkage methods (i.e., CSDLv1, the first version of the dataset), SoilGrids 2.0, and HWSD 2.0 products, effectively representing the spatial variation of soil properties across China. The prediction accuracy of most soil properties at the 0–5 cm depth interval ranged from good to moderate, with Model Efficiency Coefficients for most soil properties ranging from 0.75 to 0.32 during data-splitting validation and from 0.88 to 0.25 during independent sample validation. The wide range between the 5 % lower and 95 % upper prediction limits may indicate substantial room for improvement in current predictions. The relative importance of environmental covariates in predictions varied with soil properties and depth, indicating the complexity of interactions among multiple factors in the soil formation processes. As the soil profiles used in this study mainly originate from the Second National Soil Survey of China during 1970s and 1980s, they could provide new perspectives of soil changes together with existing maps based on 2010s soil profiles. The findings make important contributions to the GlobalSoilMap project and can also be used for regional Earth system modeling and land surface modeling to better represent the role of soil in hydrological and biogeochemical cycles in China. This dataset is freely available and can be accessed at https://doi.org/10.11888/Terre.tpdc.301235 (Shi et al, 2024). [ABSTRACT FROM AUTHOR]

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]

(1) Background: The use of multiscale prediction or the optimal scaling of predictors can enhance soil maps by applying pixel size in digital soil mapping (DSM). (2) Methods: A total of 200, 50, and 129 surface soil samples (0–30 cm) were collected by the CLHS method in three different areas, namely, the Marvdasht, Bandamir, and Lapuee plains in southwest Iran. Then, four soil properties—soil organic matter (SOM), bulk density (BD), soil shear strength (SS), and mean weighted diameter (MWD)—were measured at each sampling point as representative attributes of soil physical and chemical quality. This study examined different-scale scenarios ranging from resampling the original 30 m digital elevation model and remote sensing indices to various pixel sizes, including 60 × 60, 90 × 90, 120 × 120, and up to 2100 × 2100 m. (3) Results: After evaluating 22 environmental covariates, 11 of them were identified as the most suitable candidates for predicting soil properties based on recursive feature elimination (RFE) and expert opinion methods. Furthermore, among different pixel size scenarios for SOM, BD, SS, and MWD, the highest accuracy was achieved at 1200 × 1200 m (R2 = 0.35), 180 × 180 m (R2 = 0.67), 1200 × 1200 m (R2 = 0.42), and 2100 × 2100 m (R2 = 0.34), respectively, in Marvdasht plain. (4) Conclusions: Adjusting the pixel size improves the capture of soil property variability, enhancing mapping precision and supporting effective decision making for crop management, irrigation, and land use planning. [ABSTRACT FROM AUTHOR]

Soil organic matter (SOM) and the ratio of soil organic carbon to total nitrogen (C/N ratio) are fundamental to the ecosystem services provided by soils. Therefore, understanding the spatial distribution and relationships between the SOM components mineral-associated organic matter (MAOM), particulate organic matter (POM), and C/N ratio is crucial. Three ensemble machine learning models were trained to obtain spatial predictions of the C/N ratio, MAOM, and POM in German agricultural topsoil (0–10 cm). Parameter optimization and model evaluation were performed using nested cross-validation. Additionally, a modification to the regressor chain was applied to capture and interpret the interactions among the C/N ratio, MAOM, and POM. The ensemble models yielded mean absolute percent errors (MAPEs) of 8.2% for the C/N ratio, 14.8% for MAOM, and 28.6% for POM. Soil type, pedo-climatic region, hydrological unit, and soilscapes were found to explain 75% of the variance in MAOM and POM, and 50% in the C/N ratio. The modified regressor chain indicated a nonlinear relationship between the C/N ratio and SOM due to the different decomposition rates of SOM as a result of variety in its nutrient quality. These spatial predictions enhance the understanding of soil properties' distribution in Germany. [ABSTRACT FROM AUTHOR]

Many soil functions and processes are controlled by the soil particle size distribution. Accordingly, nationwide geoinformation on this soil property is required to enable climate-smart and resilient land management. This study presents a new deep learning approach to simultaneously model the contents of the three particle sizes of sand, silt, and clay and their variations with depth throughout the landscape. The approach allows for the consideration of the natural soil horizon boundaries and the inclusion of the surrounding landscape context of each soil profile to investigate the soil–landscape relation. Applied to the agricultural soil landscape of Germany, the approach generated a three-dimensional continuous data product with a resolution of 100 m in geographic space and a depth resolution of 1 cm. The approach relies on a patch-wise multi-target convolutional neural network (CNN) model. Genetic algorithm optimization was applied for CNN parameter tuning. Overall, the effectiveness of the CNN algorithm in generating multidimensional, multivariate, national-scale soil data products was demonstrated. The predictive performance resulted in a median root mean square error of 17.8 mass-% for the sand, 14.4 mass-% for the silt, and 9.3 mass-% for the clay content in the top ten centimeters. This increased to 20.9, 16.5, and 11.8 mass-% at a 40 cm depth. The generated data product is the first of its kind. However, even though the potential of this deep learning approach to understand and model the complex soil–landscape relation is virtually limitless, its limitations are data driven concerning the approximation of the soil-forming factors and the available soil profile data. [ABSTRACT FROM AUTHOR]

Sustainable agriculture is becoming increasingly important in the face of growing environmental challenges. One key aspect of sustainable agriculture is managing soil resources effectively. In this context, digital soil mapping (DSM) has emerged as a powerful tool to understand soil variability better and inform land management decisions. This paper proposes a comprehensive data warehouse for DSM that supports climate change analysis. Our architecture integrates frequent itemset mining (FMI) and association rules mining (ARM) to extract insights from large‐scale soil data. We review related studies in soil data warehousing and ARM, identify gaps, and propose a data warehouse architecture leveraging the galaxy multidimensional model for DSM based on the SCORPAN model, which incorporates all relevant soil forming factors. We employ and compare A‐priori, FP‐growth, and ECLAT algorithms to efficiently mine frequent itemsets and generate association rules. Our intensive experiments evaluation demonstrates that FP‐growth outperforms the other algorithms in accuracy, scalability, and speed and requires less memory. Additionally, we utilized correlation metrics for ARM, such as lift, cosine, kulc, and Imbalance ratio, to obtain the most significant and relevant association rules. These rules provide valuable insights into the complex relationships between soil properties and environmental factors, which can inform land management decisions and improve sustainable agriculture practices. This work contributes to the growing body of research on DSM and data‐driven approaches to sustainable agriculture. [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 aims of this research are (i) to compare random forest (RF), boosted regression tree (BRT), and multinomial logistic regression (MnLR) models to prepare the prediction maps of soil great group and subgroup levels, (ii) determination of the most important environmental covariates influencing the production of digital soil mapping (DSM) in an arid climate, (iii) to evaluate the efficiency of spectra indices extracted from Sentinel-2A digital images and data capability of ALOS-PALSAR radar data, and (iv) investigating the effect of sub-surface genetic horizons in the modeling of different types of soil map classes distribution. The principal component analysis method was employed to select the best set from the pool of environmental covariates (n = 46) such as geomorphometric parameters (GPs), RS indices, and diagnostic soil properties (DSP). The relative importance results indicate that Gypsic (GYP) subsurface horizon, standardized height (StH), slope length (SL), and normalized different vegetation index (NDVI) had an important role in the prediction of soil classes compared to the other selected covariates. DSM methodology was used in this research by incorporating of RF model and representative soil-forming factors that can be used for preparing the maps of soil classes in low-relief areas with a similar soil-landscape relationship. Totally, in this study places a spotlight on the profound impact of sub-surface genetic horizons, shedding light on their pivotal role in accurately modeling soil class distributions. These findings not only advance our comprehension of soil variability in arid regions but also hold immense implications for the burgeoning field of pedometrics. [ABSTRACT FROM AUTHOR]

Soil compaction is one of the most important and readily mitigated threats to soil health. Digital Soil Mapping (DSM) has emerged as an efficient method to provide broad-scale maps by combining soil information with environmental covariates. Until now, soil information input to DSM has been mainly composed of point-based quantitative measurements of soil properties and/or of soil type/horizon classes derived from laboratory analysis, point observations, or soil maps. In this study, we used field estimates of soil compaction to map soil behavior to compaction at a national scale. The results from a previous study enabled clustering of six different behaviors using the in situ field observations. Mapping potential responses to soil compaction is an effective land management tool for preventing future compaction. Random forest was used to make spatial predictions of soil behavior to compaction over cultivated soils of mainland France (about 210,000 km2). Modeling was performed at 90 m resolution. The map enabled us to spatially identify clusters of possible responses to compaction. Most clusters were consistent with known geographic distributions of some soil types and properties. This consistency was checked by comparing maps with both national and local-scale external sources of soil information. The best spatial predictors were available digital maps of soil properties (clay, silt, sand, organic carbon (SOC) content, and pH), some indicators of soil structural quality using SOC and clay content, and environmental covariates (T °C and relief-related covariates). Predicted maps were interpretable to support management recommendations to mitigate soil compactness at the soil–scape scale. Simple observational field data that are usually collected by soil surveyors, then stored and available in soil databases, provide valuable input data for digital mapping of soil behavior to compaction and assessment of inherent soil sensitivity to compaction. [ABSTRACT FROM AUTHOR]

Intensive agricultural practices lead to a deterioration in soil quality, causing a decline in farm productivity and quality, and disturbing the ecosystem balance in command areas. To achieve sustainable production and implement effective soil management strategies, understanding the state and spatial variability of soil quality is essential. The study aims to enhance the understanding of soil quality variability and provide actionable insights for sustainable soil management. In this regard, principal component analysis (PCA) and digital soil mapping were used to assess and map the spatial variability of the soil quality index (SQI) in the Cauvery command area, Mandya district, Karnataka, India. A total of 145 georeferenced soil samples were drawn at 0–15 cm depth and analyzed for physico-chemical properties. PCA was used to reduce the dataset into a minimum dataset as eight important soil indicators and to determine relative weightage factors, which were used for assessing SQI with linear and non-linear scoring methods. For spatial assessment of SQI, the random forest algorithm with environmental covariates was used to map eight soil indicators selected in the minimum dataset. The soil property maps were subjected to linear and non-linear scoring, followed by multiplying with corresponding weightage factors and summation to produce SQI maps. Results reveal that values of SQI calculated using linear scoring, range from 0.10 to 0.64, with a mean of 0.39, while non-linear scoring exhibits a wider range of 0.12 to 0.78 and a mean of 0.48. With a slight higher sensitivity index of 6.5, non-linear scoring proved to be the better scoring method compared to linear scoring. Spatial assessment shows that the R2 and LCC between the calculated and predicted SQI were higher for non-linear scoring (0.66 and 0.66) compared to linear scoring (0.60 and 0.65). The SQI maps reveal high spatial variability with more than 40 percent of soils classified as moderate-to-low index. The soils with low SQI were distributed in eastern parts, whereas western parts exhibited high-to-very-high soil quality. To achieve production goals and improve soil quality in the eastern region, sustainable soil and crop management strategies must be developed, and their effects on soil quality should be assessed. [ABSTRACT FROM AUTHOR]

Synthesizing bare soil pictures in regions with complex vegetation is challenging, which hinders the accuracy of predicting soil organic carbon (SOC) in specific areas. An SOC prediction model was developed in this study by integrating the convolutional neural network and long and short-term memory network (CNN-LSTM) algorithms, taking into consideration soil-forming factors such as climate, vegetation, and topography in Hainan. Compared with common algorithmic models (random forest, CNN, LSTM), the SOC prediction model based on the CNN-LSTM algorithm achieved high accuracy (R2 = 0.69, RMSE = 6.06 g kg−1, RPIQ = 1.96). The model predicted that the SOC content ranged from 5.49 to 36.68 g kg−1, with Hainan in the central and southern parts of the region with high SOC values and the surrounding areas with low SOC values, and that the SOC was roughly distributed as follows: high in the mountainous areas and low in the flat areas. Among the four models, CNN-LSTM outperformed LSTM, CNN, and random forest models in terms of R2 accuracy by 11.3%, 23.2%, and 53.3%, respectively. The CNN-LSTM model demonstrates its applicability in predicting SOC content and shows great potential in complex areas where obtaining sample data is challenging and where SOC is influenced by multiple interacting factors. Furthermore, it shows significant potential for advancing the broader field of digital soil mapping. [ABSTRACT FROM AUTHOR]

This study assessed the perception and use of digital applications for soil fertility management and conservation strategies among small-scale crop farmers in southwest Nigeria. A total of 376 farmers were randomly selected across the six southwest states. The data collected were analyzed using descriptive statistics. The majority of the farmers relied on perception and other non-scientific approaches such as the appearance of weeds and performance of crops in the previous season to assess soil fertility. Only 1.1% and 0.3% of the farmers assessed soil fertility through soil tests and digital applications, respectively. Most farmers adopted bush fallowing and the use of inorganic fertilizers to improve soil fertility. Although 4.8% of the farmers indicated that they had digital applications on their mobile phones, only 2.9% claimed to have used these. More than half (56.4%) of the farmers stated that a lack of awareness of the existence of digital applications and internet-enabled telephones were the reasons they have not been able to use digital applications. The majority of the farmers (97.3%) indicated their willingness to embrace the use of new farm decision digital applications which could provide more information, especially on soil fertility, if introduced. More extensive services focusing on older, less literate farmers and farmers who hitherto did not belong to any farmers' association are advocated for in order to encourage the use of digital applications and soil fertility management and conservation practices. [ABSTRACT FROM AUTHOR]

This review focuses on digital soil organic carbon (SOC) mapping at regional or national scales in spatial resolutions up to 1 km using open data remote sensing sources, emphasizing its importance in achieving United Nations' Sustainable Development Goals (SDGs) related to hunger, climate action, and land conservation. The literature review was performed according to scientific studies indexed in the Web of Science Core Collection database since 2000. The analysis reveals a steady rise in total digital soil mapping studies since 2000, with digital SOC mapping studies accounting for over 20% of these studies in 2023, among which SDGs 2 (Zero Hunger) and 13 (Climate Action) were the most represented. Notably, countries like the United States, China, Germany, and Iran lead in digital SOC mapping research. The shift towards machine and deep learning methods in digital SOC mapping has surged post-2010, necessitating environmental covariates like topography, climate, and spectral data, which are cornerstones of machine and deep learning prediction methods. It was noted that the available climate data primarily restrict the spatial resolution of digital SOC mapping to 1 km, which typically requires downscaling to harmonize with topography (up to 30 m) and multispectral data (up to 10–30 m). Future directions include the integration of diverse remote sensing data sources, the development of advanced algorithms leveraging machine learning, and the utilization of high-resolution remote sensing for more precise SOC mapping. [ABSTRACT FROM AUTHOR]

Faced with increasing global soil degradation, spatially explicit data on cropland soil organic matter (SOM) provides crucial data for soil carbon pool accounting, cropland quality assessment and the formulation of effective management policies. As a spatial information prediction technique, digital soil mapping (DSM) has been widely used to spatially map soil information at different scales. However, the accuracy of digital SOM maps for cropland is typically lower than for other land cover types due to the inherent difficulty in precisely quantifying human disturbance. To overcome this limitation, this study systematically assessed a framework of “information extraction-feature selection-model averaging” for improving model performance in mapping cropland SOM using 462 cropland soil samples collected in Guangzhou, China in 2021. The results showed that using the framework of dynamic information extraction, feature selection and model averaging could efficiently improve the accuracy of the final predictions (R2: 0.48 to 0.53) without having obviously negative impacts on uncertainty. Quantifying the dynamic information of the environment was an efficient way to generate covariates that are linearly and nonlinearly related to SOM, which improved the R2 of random forest from 0.44 to 0.48 and the R2 of extreme gradient boosting from 0.37 to 0.43. Forward recursive feature selection (FRFS) is recommended when there are relatively few environmental covariates (500). The Granger-Ramanathan model averaging approach could improve the prediction accuracy and average uncertainty. When the structures of initial prediction models are similar, increasing in the number of averaging models did not have significantly positive effects on the final predictions. Given the advantages of these selected strategies over information extraction, feature selection and model averaging have a great potential for high-accuracy soil mapping at any scales, so this approach can provide more reliable references for soil conservation policy-making.

Abstract Obtaining accurate spatial maps of soil organic carbon (SOC) in farmlands is crucial for assessing soil quality and achieving precision agriculture. The cropping system is an important factor that affects the soil carbon cycle in farmlands, and different agricultural managements under different cropping systems lead to spatial heterogeneity of SOC. However, current research often ignores differences in the main controlling factors of SOC under different cropping systems, especially when the cropping pattern is complex, which is not conducive to farmland zoning management. This study aims to (i) obtain the spatial distribution map of six cropping systems by using multi-phase HJ-CCD satellite images; (ii) explore the stratified heterogeneous relationship between SOC and environmental variables under different cropping systems by using the Cubist model; and (iii) predict the spatial map of SOC. The Xiantao, Tianmen, and Qianjiang cities, which are the core agricultural areas of the Jianghan Plain, were selected as the study area. Results showed that the SOC content in rice–wheat rotation was the highest among the six cropping systems. The Cubist model outperformed random forest, ordinary kriging, and multiple linear regression in SOC mapping. The results of the Cubist model showed that cropping system, climate, soil attributes, and vegetation index were important influencing factors of SOC in farmlands. The main controlling factors of SOC under different cropping systems were different. Specifically, summer crop types had a greater influence on spatial variations in SOC than winter crops. Paddy–upland rotation was more affected by river distance and NDVI, while upland–upland rotation was more affected by irrigation-related factors. This work highlights the differentiated main controlling factors of SOC under different cropping systems and provides data support for farmland zoning management. The Cubist model can improve the prediction accuracy of SOC under complex cropping systems.

ABSTRACTThis research focuses on understanding the spatial variation of Soil Organic Matter (SOM) and pH levels in the North of Morocco. The study employs a comprehensive approach to enhance predictive modelling, incorporating the Boruta algorithm for effective environmental covariates selection and optimizing model parameters through hyperparameter optimization. Utilizing a Random Forest (RF) model with remote sensing indices and topographic features, the research predicts SOM and pH to identify key contributors to their spatial variability. SOM prediction saw significant success, with a notable correlation to remote sensing indices such as the RVI, NDVI, and TNDVI. These indices, indicative of vegetation health and productivity, emerged as primary influencers of SOM. In comparison, the influence of topographic features like elevation, slope, and aspect was found to be less significant. Conversely, predicting pH was challenging due to the minimal spatial variability within the dataset. Addressing this limitation could involve dataset expansion or alternative models for low-correlated data handling. Despite the RF model’s limited efficacy in pH prediction, an observable correlation between SOM and pH was identified, consistent with prior research. Areas with higher SOM exhibited lower pH values, indicating relative soil acidification from organic matter decomposition. The study’s RF model demonstrated potential in SOM prediction using remote sensing indices, but enhancing pH prediction is essential. Future research may explore dataset expansion, diverse sampling, or testing alternative predictive models for better performance with low-correlated datasets. The study offers valuable insights for advanced predictive model development and enriches understanding of soil management practices.

ABSTRACTSpatial prediction is essential for obtaining the spatial distribution of geographic variables and selecting appropriate covariates for this process can be challenging, especially for non-expert users. For easing the burden of selecting the appropriate covariates, two case-based reasoning strategies, namely the most-similar-case and covariate-classification strategies, have been proposed for automated covariate selection. The former may suggest nonessential covariates due to its case-level reasoning way. And the latter with covariate-level reasoning may overlook related covariates and recommend fewer covariates than the case-level reasoning. In this study, we propose a new strategy of integrating case-level and covariate-level reasoning to effectively leverage the strengths of both previous strategies while also addressing their limitations. The proposed strategy is validated through a case study of automatically selecting covariates for digital soil mapping under reasoning with a case base containing 189 cases. The leave-one-out evaluation demonstrated that our proposed strategy outperformed the previous two strategies.

Mapping soil properties in sub-watersheds is critical for agricultural productivity, land management, and ecological security. Machine learning has been widely applied to digital soil mapping due to a rapidly increasing number of environmental covariates. However, the inclusion of many environmental covariates in machine learning models leads to the problem of multicollinearity, with poorly understood consequences for prediction performance. Here, we explored the effects of variable selection on the prediction performance of two machine learning models for multiple soil properties in the Haihun River sub-watershed, Jiangxi Province, China. Surface soils (0–20 cm) were collected from a total of 180 sample points in 2022. The optimal covariates were selected from 40 environmental covariates using a recursive feature elimination algorithm. Compared to all-variable models, the random forest (RF) and extreme gradient boosting (XGBoost) models with variable selection improved in prediction accuracy. The R2 values of the RF and XGBoost models increased by 0.34 and 0.47 for the soil organic carbon, by 0.67 and 0.62 for the total phosphorus, and by 0.43 and 0.62 for the available phosphorus, respectively. The models with variable selection presented reduced global uncertainty, and the overall uncertainty of the RF model was lower than that of the XGBoost model. The soil properties showed high spatial heterogeneity based on the models with variable selection. Remote sensing covariates (particularly principal component 2) were the major factors controlling the distribution of the soil organic carbon. Human activity covariates (mainly land use) and organism covariates (mainly potential evapotranspiration) played a predominant role in driving the distribution of the soil total and soil available phosphorus, respectively. This study indicates the importance of variable selection for predicting multiple soil properties and mapping their spatial distribution in sub-watersheds. [ABSTRACT FROM AUTHOR]

Texture, soil organic matter (SOM), and soil depth (SoD) are crucial properties in forest management because they can supply spatial information on forest site productivity and guide fertilizer applications. However, soil properties possess an inherent uncertainty that must be mapped to enhance decision making in management applications. Most digital soil mapping predictions primarily concentrate on the mean of the distribution, often neglecting the estimation of local uncertainty in soil properties. Additionally, there is a noticeable scarcity of practical soil examples to demonstrate the prediction uncertainty for the benefit of forest managers. In this study, following a digital soil mapping (DSM) approach, a Quantile Regression Forest (QRF) model was developed to generate high-resolution maps and their uncertainty regarding the texture, SoD, and SOM, which were expressed as standard deviation (Sd) values. The results showed that the SOM (R2 = 0.61, RMSE = 2.03% and with an average Sd = 50%), SoD (R2 = 0.74 and RMSE = 19.4 cm), clay (R2 = 0.63, RMSE = 10.5% and average Sd = 29%), silt (R2 = 0.59, RMSE = 6.26% and average Sd = 33%), and sand content (R2 = 0.55, RMSE = 9.49% and average Sd = 35%) were accurately estimated for forest plantations in central south Chile. A practical demonstration of precision fertilizer application, utilizing the predictive distribution of SOM, effectively showcased how uncertainty in soil attributes can be leveraged to benefit forest managers. This approach holds potential for optimizing resource allocation and maximizing economic benefits. [ABSTRACT FROM AUTHOR]

In response to the growing societal awareness of the critical role of healthy soils, there has been an increasing demand for accurate and high-resolution soil information to inform national policies and support sustainable land management decisions. Despite advancements in digital soil mapping and initiatives like GlobalSoilMap, quantifying soil variability and its uncertainty across space, depth and time remains a challenge. Therefore, maps of key soil properties are often still missing on a national scale, which is also the case in the Netherlands. To meet this challenge and fill this data gap, we introduce BIS-4D, a high-resolution soil modeling and mapping platform for the Netherlands. BIS-4D delivers maps of soil texture (clay, silt and sand content), bulk density, pH, total nitrogen, oxalate-extractable phosphorus, cation exchange capacity and their uncertainties at 25 m resolution between 0 and 2 m depth in 3D space. Additionally, it provides maps of soil organic matter and its uncertainty in 3D space and time between 1953 and 2023 at the same resolution and depth range. The statistical model uses machine learning informed by soil observations amounting to between 3815 and 855 950, depending on the soil property, and 366 environmental covariates. We assess the accuracy of mean and median predictions using design-based statistical inference of a probability sample and location-grouped 10-fold cross validation (CV) and prediction uncertainty using the prediction interval coverage probability. We found that the accuracy of clay, sand and pH maps was the highest, with the model efficiency coefficient (MEC) ranging between 0.6 and 0.92 depending on depth. Silt, bulk density, soil organic matter, total nitrogen and cation exchange capacity (MEC of 0.27 to 0.78), and especially oxalate-extractable phosphorus (MEC of -0.11 to 0.38) were more difficult to predict. One of the main limitations of BIS-4D is that prediction maps cannot be used to quantify the uncertainty in spatial aggregates. We provide an example of good practice to help users decide whether BIS-4D is suitable for their intended purpose. An overview of all maps and their uncertainties can be found in the Supplement. Openly available code and input data enhance reproducibility and help with future updates. BIS-4D prediction maps can be readily downloaded at 10.4121/0c934ac6-2e95-4422-8360-d3a802766c71. BIS-4D fills the previous data gap of the national-scale GlobalSoilMap product in the Netherlands and will hopefully facilitate the inclusion of soil spatial variability as a routine and integral part of decision support systems. [ABSTRACT FROM AUTHOR]

Spatial prediction is essential for obtaining the spatial distribution of geographic variables and selecting appropriate covariates for this process can be challenging, especially for non-expert users. For easing the burden of selecting the appropriate covariates, two case-based reasoning strategies, namely the most-similar-case and covariate-classification strategies, have been proposed for automated covariate selection. The former may suggest nonessential covariates due to its case-level reasoning way. And the latter with covariate-level reasoning may overlook related covariates and recommend fewer covariates than the case-level reasoning. In this study, we propose a new strategy of integrating case-level and covariate-level reasoning to effectively leverage the strengths of both previous strategies while also addressing their limitations. The proposed strategy is validated through a case study of automatically selecting covariates for digital soil mapping under reasoning with a case base containing 189 cases. The leave-one-out evaluation demonstrated that our proposed strategy outperformed the previous two strategies. [ABSTRACT FROM AUTHOR]

This research focuses on understanding the spatial variation of Soil Organic Matter (SOM) and pH levels in the North of Morocco. The study employs a comprehensive approach to enhance predictive modelling, incorporating the Boruta algorithm for effective environmental covariates selection and optimizing model parameters through hyperparameter optimization. Utilizing a Random Forest (RF) model with remote sensing indices and topographic features, the research predicts SOM and pH to identify key contributors to their spatial variability. SOM prediction saw significant success, with a notable correlation to remote sensing indices such as the RVI, NDVI, and TNDVI. These indices, indicative of vegetation health and productivity, emerged as primary influencers of SOM. In comparison, the influence of topographic features like elevation, slope, and aspect was found to be less significant. Conversely, predicting pH was challenging due to the minimal spatial variability within the dataset. Addressing this limitation could involve dataset expansion or alternative models for low-correlated data handling. Despite the RF model's limited efficacy in pH prediction, an observable correlation between SOM and pH was identified, consistent with prior research. Areas with higher SOM exhibited lower pH values, indicating relative soil acidification from organic matter decomposition. The study's RF model demonstrated potential in SOM prediction using remote sensing indices, but enhancing pH prediction is essential. Future research may explore dataset expansion, diverse sampling, or testing alternative predictive models for better performance with low-correlated datasets. The study offers valuable insights for advanced predictive model development and enriches understanding of soil management practices. [ABSTRACT FROM AUTHOR]

This document is an editorial from a special issue of the journal "Land" focused on digital soil mapping, decision support tools, and soil monitoring systems in the Mediterranean region. The editorial highlights the evolving role of soil surveyors in the digital era and the increased collaboration between scientists in advancing soil mapping efforts. It emphasizes the importance of integrating legacy and new soil data to enhance our understanding of soil properties and support sustainable soil management. The editorial also provides an overview of seven abstracts from recent studies that showcase the innovative use of machine learning in digital soil mapping and the assessment of soil properties in the Mediterranean region. These studies demonstrate the potential of machine learning models to improve soil management and environmental sustainability. [Extracted from the article]

An analysis of the current potential range of the Pleistocene relict plant species Schoenus ferrugineus and modeling of changes in its future range under moderate (RCP4.5) and strong (RCP8.5) climate change in the middle and second half of the 21st century were carried out. The MaxEnt program was used for modeling. Climate variables from CHELSA Bioclim, the global digital soil mapping system SoilGrids, and a digital elevation model were used as predictors. Modeling has shown that climate change will lead to a significant reduction in the suitability of S. ferrugineus habitat conditions by the mid-21st century. The predicted changes in the distribution of habitats of S. ferrugineus, a diagnostic species of calcareous mires and an indicator of their ecological state, indicate a possible strong transformation of wetland complexes in the Southern Urals region even under moderate climate change. A reduction in the distribution of S. ferrugineus at the eastern limit of its range will also be facilitated by more frequent extreme droughts. To maintain the distribution of S. ferrugineus on the eastern border of its range, a number of measures are proposed to mitigate the negative consequences of climate change, contributing to the preservation of the hydrological regime of calcareous mires. [ABSTRACT FROM AUTHOR]

Obtaining accurate spatial maps of soil organic carbon (SOC) in farmlands is crucial for assessing soil quality and achieving precision agriculture. The cropping system is an important factor that affects the soil carbon cycle in farmlands, and different agricultural managements under different cropping systems lead to spatial heterogeneity of SOC. However, current research often ignores differences in the main controlling factors of SOC under different cropping systems, especially when the cropping pattern is complex, which is not conducive to farmland zoning management. This study aims to (i) obtain the spatial distribution map of six cropping systems by using multi-phase HJ-CCD satellite images; (ii) explore the stratified heterogeneous relationship between SOC and environmental variables under different cropping systems by using the Cubist model; and (iii) predict the spatial map of SOC. The Xiantao, Tianmen, and Qianjiang cities, which are the core agricultural areas of the Jianghan Plain, were selected as the study area. Results showed that the SOC content in rice–wheat rotation was the highest among the six cropping systems. The Cubist model outperformed random forest, ordinary kriging, and multiple linear regression in SOC mapping. The results of the Cubist model showed that cropping system, climate, soil attributes, and vegetation index were important influencing factors of SOC in farmlands. The main controlling factors of SOC under different cropping systems were different. Specifically, summer crop types had a greater influence on spatial variations in SOC than winter crops. Paddy–upland rotation was more affected by river distance and NDVI, while upland–upland rotation was more affected by irrigation-related factors. This work highlights the differentiated main controlling factors of SOC under different cropping systems and provides data support for farmland zoning management. The Cubist model can improve the prediction accuracy of SOC under complex cropping systems. [ABSTRACT FROM AUTHOR]

The presence of the petrogypsic horizon is an impediment to developing agriculture in the Sahara. It hinders the soil's ability to store water and root development of crops. The petrogypsic horizon is commonly difficult to map due to its location either on the surface or at depth. This study used logistic regression-kriging and logistic regression models to map the petrogypsic horizon occurrence probability using 466 observations over an area of 22 573 ha in the Sahara Desert of Algeria. The models included remote sensing indices and topographic variables as environmental covariates. The accuracy of models was verified by the area under the curve (AUC). A binary map was produced by applying a threshold of 0.7 on the most performant probability map. Our results showed that logistic regression-kriging performed the best (AUC = 0.88), due to the consideration of residual spatial correlation in the model. The grain size index covariate was the most relevant compared to topographic variables, which showed the usefulness of spectral indices. Based on the binary map, the risk associated with the presence of the petrogypsic horizon was limited, representing 26% of the study area. In the Sahara Desert, though the petrogypsic horizon was weakly correlated with the tested environmental covariates, the use of satellite images and residual autocorrelation in a predictive modelling approach improved the mapping and thus risk assessment of the petrogypsic horizon. [ABSTRACT FROM AUTHOR]

Recently, numerous spatial prediction methods with diverse characteristics have been developed. Selecting an appropriate spatial prediction method, along with its data preprocessing and parameter settings, presents a challenging task for many users, especially for non-experts. This paper addresses this challenge by exploring the potential of automated machine learning method proposed in artificial intelligent domain to automatically determine the most suitable method among various machine learning methods. As a case study, the automated machine learning method was applied to predict the spatial distribution of soil subgroups in Heshan farm. A total of 110 soil samples and 10 terrain variables were utilized in the designed experiments. To evaluate the performance, the proposed method was compared to each machine learning method with default parameters values or parameters determined by expert knowledge. The results showed that the proposed method typically achieved higher accuracy scores than the two alternative methods. This suggests that automated machine learning performs effectively in scenarios where numerous machine learning methods are available and offers practical utility in reducing the dependence on users' expertise in spatial prediction. However, a more robust automated framework should be developed to encompass a broader range of spatial prediction methods, such as spatial statistic methods, rather than only focusing on machine learning methods. [ABSTRACT FROM AUTHOR]

Soil organic carbon (SOC) is central to the functioning of terrestrial ecosystems, has climate mitigation potential and provides several benefits for soil health. Understanding the spatial distribution of SOC can help formulate sustainable soil management practices. Digital soil mapping (DSM) uses advanced statistical and geostatistical methods to estimate soil properties across large areas. DSM integrates climate data, topographic features, geology, legacy soil maps, land management and remote sensing data. Bare soil spectra may reflect the presence of particular soil components, making satellite derived spectra suitable predictors of SOC. Bare soil spectra derived from Sentinel-2 were used to estimate SOC concentration (SOC%) and granulometric fractions in the plough layer (0–30 cm) of agricultural parcels in northern Belgium. Thereafter, the estimation performance of SOC% was compared for three DSM models: one with bare soil spectra, one with environmental covariates (topography, granulometry and vegetation), and a combined model with bare soil spectra and environmental covariates. The estimation performance of sand, silt and clay fractions using bare soil spectra from the spring seedbed (R2: 0.53–0.74; RPD: 1.49–2.05; RPIQ: 1.52–2.39) was higher than that of SOC% (R2: 0.16; RPD: 1.08; RPIQ: 1.32). The highest estimation performance of SOC% was obtained for a DSM model including all covariates (R2: 0.28; RPD: 1.18; RPIQ: 1.44), but the contribution of spring seedbed spectra to a model containing environmental covariates was small. The results provide valuable insights for refining soil property estimation using DSM with spectral and environmental covariates.

Fine-scale spatial distribution of soil physicochemical properties is crucial for soil quality management, agriculture planning and geotechnical engineering. Existing soil map databases are usually developed in national scale, potentially leading to issues of coarse resolution and restricted applicability in fine-scaled studies. For broad-scale area, conventional digital soil mapping methods are challenging due to the lack of representative soil profiles and their uneven distribution. This study addresses these challenges by developing a downscaling-based framework to map 11 soil physicochemical properties at a 30-m resolution. The methodology involves constructing regression models using soil properties derived from coarse national soil maps and soil-forming covariates. Predictions are subsequently refined using fine-resolution covariates. To capture fine-scale spatial variability across diverse landscapes, 30-m resolution remote sensing data, relief, climate, and spatial covariates were prepared. The fine-gridding framework employs tree-based machine learning models, particularly random forest (RF), to enhance prediction accuracy. Evaluation of the fine-gridded soil maps demonstrated commendable results, with high consistency in summary statistics, semi-variograms, and evaluation metrics compared to reference maps. Mass preservation of RF-predicted maps exhibit high performances with both concordance correlation coefficients (CCCs) and coefficients of determination (R2s) exceeding 0.9 across most scenarios. This study provides a robust approach for enhancing the spatial resolution of soil maps, facilitating their use in local fine-scale applications. The proposed strategy offers a valuable solution for broad-scale areas that require finer soil maps but lack sufficient qualified soil profiles.

This study aims to identify a key soil property influencing farmland prices as a soil production capital indicator for soil security. There is generally a strong relationship between farmland price and production capital, as well as with biomass production. From a soil-centric point of view, under dryland agriculture in Australia, we hypothesise that soil available water capacity (AWC) is a key soil property influencing farmland prices and consequently can act as a soil production capital indicator. Here, we relate soil AWC of Australia to median farmland prices using digital soil maps of AWC. The results revealed that AWC from 0 to 100 cm is significantly associated with farmland prices (r = 0.31, p-value

The crucial role of soil in global food production and its multifaceted contributions to ecosystem services underscore the need for a comprehensive evaluation framework. This study presents a novel approach to soil evaluation by integrating the assessment of soil functions, services, and threats into a unified metric which quantifies the capital dimension which is termed ''soil management capital''. To achieve this, we focused on three key roles: a threat, evaluated by the pH regulation costs, factoring in different buffering capacities and liming/sulphuring treatments to attain an optimal reference pH, a function, representing soil's role as a reservoir and regulator of nutrients; and a service, considering available water holding capacity (AWC) for water security. Integrating these activities culminates in a soil capital metric, providing a comprehensive understanding of the soil's economic value in US$/ha units on a pixel-by-pixel basis across Australia. To achieve this, we analysed spatial information from the Soil and Landscape Grid of Australia (SLGA) with Google Earth Engine (GEE). Our study generated detailed maps illustrating pH regulation costs, OC, nutrients (N, P, K), and AWC across Australia. The integration of these roles revealed regions of large soil capital in the Australian eastern coast and Tasmania.

Soil organic carbon (SOC) plays a crucial role in the global carbon cycle and in maintaining soil functions in the context of land use and climate change. Understanding the spatial distribution of SOC is essential for the management of agricultural land to optimize soil health and carbon storage. In this study, we investigated the spatial distribution of SOC in an agricultural lowland area of the Lombardy region, Italy, using machine learning (ML) techniques combined with residual kriging. ML models, including the artificial neural network (ANN), extreme learning machine (ELM), and random forest (RF), were trained on 120 SOC observations and eight environmental variables to predict SOC values across the study area. The performance of this ML approach was assessed using a ten-fold nested cross-validation process. The ELM and RF models showed better predictive performances based on the concordance correlation coefficient and root mean square error (RMSE), with RF slightly outperforming ELM based on the RMSE. The residuals of each iteration from the ML models were interpolated by ordinary kriging (OK) and added to the ML-based trend model in a hybrid regression-kriging approach. This approach which accounted for the spatial autocorrelation of the prediction residuals, resulting in a marginally improved prediction accuracy in the ML models. In addition, we found that vertical distance to the channel network and channel network base level are important predictor variables that should be considered in future digital soil models for SOC in lowland areas, given their importance in this study. Furthermore, this study highlights that predicted SOC values were low, particularly in Luvisols, which can be explained by the long history of agricultural land use depleting SOC due to agricultural management and loss of organic plant residues. The prediction maps depicted spatial variation and patterns of SOC in the study area. Our findings may help to refine soil management practices and contribute to improving soil health and carbon sequestration in agricultural lowland areas.

For the international digital soil mapping (DSM) community, adequate spatial estimates of nitrogen (N) mineralization have yet to be generated. This is due, in part, to an inability to capture critical N controls at the regional and provincial scales. While the influence of climate, vegetation, and relief are accessible predictors in DSM, the effect of soil management is known for its important influence on N dynamics, but has hitherto been elusive for soil mappers. For the purpose of producing N maps to inform N fertilizer management, the intention of this study was to determine the importance of novel crop frequency layers, as a proxy for soil management, through the development of provincial scale DSMs of total nitrogen (TN), biological nitrogen availability (BNA) and the estimate of N mineralization over a growing season (GSN) as calculated from TN and BNA results. Crop frequency covariates were developed that estimated the frequency a particular crop type was planted over a 10-year period, thus capturing cropping system and tillage intensity. TN results were 27% higher using crop frequency layers and the support vector machine learner, with a Lin’s concordance correlation coefficient (concordance) of 0.45. BNA predictions increased by 24% using crop frequency layers and the stochastic gradient boosting learner for a final concordance of 0.45. GSN showed the least improvement using the crop frequency layers (6%) but resulted in the highest concordance (0.47) using crop frequency layers with the stochastic gradient boosting learner. The stable N pool, represented by TN, showed climate covariates with the highest importance; whereas, the labile pool, based on BNA measures, was best predicted and controlled by organism covariates. The successful inclusion of crop frequency layers into N maps indicated that the number of times in forages and potatoes over a 10-year period was of the greatest importance. As tillage intensity is most pronounced in potatoes, and as forages contribute to increased biomass and building soil organic matter levels, the results showed that increasing the number of years in forages had a positive correlation with GSN and the stable and labile N pools.

Accurately quantifying high-resolution field-scale soil organic carbon (SOC) stocks is challenging yet crucial for improving site-specific land management and carbon accounting. This challenge is even greater when the study units are large heterogenous ranches. This study utilizes a digital soil mapping (DSM) approach and a U.S. legacy soil dataset, combined with soil, climate, biotic, and topographic covariate datasets, to design a targeted soil sampling plan for acquiring local samples. The resulting local samples were then used in combination with the legacy data to build optimal ranch-scale SOC stock models. We provide an example of this approach using ranch units in the western U.S. as a case study. In our approach we first applied a clustering analysis to generate spatial clusters. This was followed by adopting a conditioned Latin hypercube sampling scheme within each cluster, to generate sets of strategically selected local sampling points. The local samples required for improved estimates of SOC stocks were determined to have a sample size of 15 and 40 soil cores, within the respective 13 and 36 km2 parcels. While our modeling results for SOC concentrations at a relatively homogeneous site in eastern Montana showed a significant two-fold improvement in model fit when using individually selected calibration datasets for each point, as opposed to selecting the calibration dataset as a whole at the ranch level, the disparity between the pixel- and ranch-based models was inconsequential for the other two ranch sites in Colorado that were more spatially diverse in terms of land management and vegetation cover. Compared to SOC concentration models (R2 between 0.3 and 0.7), the performance of models for bulk density (BD) (R2

Estimations of the patterns and controls of soil organic carbon (SOC) could provide instructive insights into the potential impact of future global change on soil carbon (C). In this work, we combined GeoDetector and random forest (RF) to estimate SOC stocks in Chinese grassland ecosystems with uncertainty assessments, and identified a network of cross-correlated environmental covariates for determining SOC based on a dataset collected from 813 sampling sites (0–20 cm) collected from 2000 to 2014. We predicted that 17.50 Pg C was stored in Chinese grasslands to a depth of 20 cm and that the average SOC density was 4.69 kg C/m−2(−|−). The effectiveness of using RF to predict SOC was demonstrated by an accuracy assessment based on 10-fold cross-validation, with a ratio of performance to deviation (RPD) of 2.89. The SOC stocks in southern China were lower than those in northern China. A high SOC density was found in northeastern China and on the Qinghai–Tibet Plateau. Soil properties had the strongest direct effect on SOC. Climate was significantly negatively associated with SOC and indirectly affected SOC via its effect on soil properties. Topography had a significant direct impact on vegetation, but its direct effect on SOC was relatively weak. This study emphasizes the patterns and heterogeneity of SOC stocks as well as the relative significance of climate, vegetation, soil characteristics, topography, and their complex interrelationships in controlling SOC. These results may provide a theoretical foundation for developing sustainable management systems and calibrating C process models.

The advent of cloud computing platforms (e.g., Google Earth Engine (GEE)) and the massive amounts of optical and radar Earth Observation (EO) data hosted by these platforms present new opportunities for mapping soil pH at large scales. However, existing studies generally lack consensus on the effects of satellite and radar sensor parameters on GEE-based soil pH prediction models. In this study, we assessed the suitability of long-term radar (C-band Sentinel-1 and L-band PALSAR-1/2) and optical (Sentinel-2) EO data on GEE for the digital mapping of soil pH on a continental (Europe) scale and determined the most appropriate radar sensor parameters. Thirteen scenarios with different data configurations were simulated and combined with the 2018 LUCAS soil database and two machine learners (boosted regression trees and extreme gradient boosting) to develop soil prediction models. Results showed that the selection of modeling techniques, satellite sensors and radar system parameters largely affected the model output. Models involving a single polarization mode of PALSAR-1/2 data performed the worst (RPD = 1.24). Models based on Sentinel-1 data performed better than those built using PALSAR-1/2 data. The model performance was improved when a model involved more polarization bands, orbital directions, and band frequencies. The multiband model built using the two radar datasets achieved a comparable accuracy to the model based on optical data. Moreover, the model that fused radar-optical data achieved better results, with RPD values of 1.56 and 1.46 for the models with and without radar data, respectively; its performance was comparable to that of models built with commonly used variables (topography and climate). The analysis of importance indicated that long-term optical and radar EO data on GEE were important in our model. The modelling of soil pH at the continental scale largely benefits from GEE. The predicted maps exhibited strong spatial heterogeneity among different biogeographic regions, with similar spatial patterns under different modelling scenarios.

Digital soil mapping relies on statistical relationships between soil profile observations and environmental covariates at the sample locations. However, inherent limitations of legacy soil profiles, such as inaccurate georeferencing, could frequently introduce location errors into these soil profiles that affect the quality of digital soil mapping. To address this challenge, this study focuses on reducing the location error of legacy soil profiles and evaluating the resulting impact on digital soil mapping. We improved the agreement between the detailed descriptive information of legacy soil profiles and the relatively accurate environmental covariates (such as elevation, slope, and land use) to reduce the location errors of the legacy soil profiles. Quantile regression forest models were constructed to predict soil properties and their uncertainty using legacy soil profiles before and after location error correction. Our results demonstrate that for the majority of soil variables, correcting positional errors in legacy soil profiles to some extent enhances the accuracy of the digital soil mapping. The largest improvement was found for soil organic carbon at 0–5 cm depth interval, with 21 % increase of MEC. The impact of reduced location error is particularly noteworthy in regions characterized by complex terrain. In addition, the details of the predicted maps of legacy soil profiles were improved after correcting for positional errors, which better represent the spatial variation of soil properties across China. Besides, we also found that elevation was the primary controlling factor for correcting location error of legacy soil profiles. This research presents a step towards producing high-resolution and high-quality spatial soil datasets, which can provide essential support for soil management and ensure future soil security.

Imbalanced datasets are one of the main challenges in digital soil mapping. For these datasets, machine learning techniques commonly overestimate the majority classes and underestimate the minority ones. In general, this generates maps with poor precision and unrealistic results. Considering these maps for land use decision-making can have dire consequences. This is the case of acid sulfate (AS) soils, a type of harmful soil that can generate serious environmental damage when drained in agricultural or forestry activities. Therefore, it is necessary to create high-precision maps to avoid environmental damage. Although most soil class datasets in nature are imbalanced, this problem has hardly been studied. One of the main objectives of this work is the evaluation of different techniques to address the problem of imbalanced datasets. The methods considered to balance the dataset are an undersampling technique, the addition of more samples, and the combination of both. For increasing the number of samples from the minority class, we develop a new technique by creating artificial samples from the quaternary geological map. The method used for the modeling is Random Forest, one of the best methods for the classification of AS soils. Balancing the dataset improves the performance of the model in all the studied cases, where the values of the metrics for both classes are above 80%. The consideration of artificial non-AS soil samples improves the prediction of the model for the AS soils. Furthermore, we create AS soil probability maps for the four balanced datasets and the imbalanced dataset. The modeled AS soil probability maps created from the balanced datasets have high precision. A detailed comparison between the maps is made. The predictions of some of these maps match between 75%–80% of the study area. In addition, the extent of the AS soils obtained in all the cases is compared with the extent of the AS soils in the conventionally produced occurrence map. The good results of this study confirm the importance of balancing the dataset to improve the prediction and classification of AS soils.

Although agricultural intensification has generally increased crop yields, it also resulted in a range of environmental issues. These include increased erosion rates and declined soil organic carbon (SOC) stocks. In order to improve our understanding on how erosion impacts the overall SOC storage capacity of croplands, this study analyses the 3D distribution of SOC as a function of water and tillage erosion in a conventionally ploughed field in the Belgian silt loam region. We present a novel methodological framework to integrate the output of an advanced erosion model as a co-variate within a Digital Soil Mapping (DSM) approach with the objective to create detailed SOC maps. More precisely, we combined (i) the Water and Tillage Erosion Model and Sediment Delivery Model (WaTEM/SEDEM), simulating spatial patterns of soil erosion and sediment deposition due to water and tillage erosions, with (ii) a SOC sampling campaign, resulting in a SOC spatial distribution model that considers both types of erosion. The results show that, as compared to plateaus, SOC stocks are nearly half as large along eroding slopes (i.e. convex slope positions affected by tillage erosion and steep slopes affected by water erosion). Yet, they are up to twice as large in areas characterized by sediment deposition (i.e. concave positions due to tillage erosion and foot slope positions due to water erosion). Our results further show that tillage erosion has a significant influence on the SOC stocks in the top 0.7 m and in particular on the top 0.4 m. The influence of water erosion is less strong but mostly significant along the entire depth profile. Overall, this work demonstrates the relevance of considering different erosion processes when aiming to predict spatial patterns of SOC. Considering the top 1 m and a WaTEM/SEDEM application at a resolution of 10 m our 3D SOC modelling approach obtained a coefficient of determination (R2) of 0.62, a relative root mean square error (RRMSE) of 30.5 % and a relative mean absolute error (RMAE) of 26.0 %. While future work may likely lead to further improvements e.g. a more detailed SOC sampling network along the foot slopes and thalwegs in order to obtain a more spatially detailed prediction of the area characterized by depositions due to water erosion, we demonstrate the great potential of existing erosion and deposition models in doing so.

A new method for constructing detailed maps of the soil organic matter (SOM) distribution in the top layer of arable land has been developed and proposed. The method is based on the theory of spectral neighborhood of the soil line (SNSL) and the technology of constructing a multitemporal soil line (MSL). The method is based on the processing of big remote sensing data (BRSD) from 1984 to 2023. Filtering of BRSD and detection of bare soil surface (BSS) is carried out on the basis of neural networks. The method was implemented for BSS with an area about of 79,000 ha with a spatial resolution of 30 m in the Mtsensk district (Oryol Oblast, Russia). Verification was provided by four independent field surveys (which were carried out using three various methods) in different years (2022–2023). The regression is described by a polynomial of degree 2. The coefficient of determination (R2) of the regression was 0.8. The proposed method can be widely used for mapping of SOM in the areas of transition from leached chernozems (Luvic Chernozems) to sod-podzolic (Albic Retisol) and gray forest (Luvic Phaeozems) soils or in similar nature conditions.

This article discusses the intersection of data science and sustainable development, highlighting the potential for data science tools to contribute to the achievement of the United Nations' 2030 Agenda for Sustainable Development. The article emphasizes the need for empirical evidence regarding the advantages of data science in sustainability, but also acknowledges the numerous cases where data science has been implemented to promote a transition towards a more sustainable economy. The article presents a special issue of Expert Systems that explores the application of data science methodologies to the concept of sustainability, with articles covering topics such as digital soil mapping, artificial intelligence governance, digital innovation in large corporations, sustainable supply chain design for electric vehicle battery production, quantum computing, and edge computing in professional farmer cultivation. The authors hope that this compilation of research will inspire further initiatives and collaboration between researchers, practitioners, and policymakers in the field of data science for sustainable development. [Extracted from the article]

IntroductionUnderstanding the particle size distribution (PSD) is of great importance for plant growth and soil management. In recent years, the science of soil has witnessed a significant increase in digital soil mapping (DSM) activities. In this regard, machine learning models (ML) have emerged as an alternative and tool for DSM, which are mainly used for data mining and pattern recognition purposes, and are now widely used for regression and classification tasks in all fields of science. Hence, this study was undertaken to spatially model sand, silt, and clay particles utilizing machine learning models such as Random Forest (RF), Support Vector Regression (SVR), and the Co-Kriging geostatistical model. Additionally, auxiliary variables with high spatial resolution were incorporated into the analysis. This investigation was conducted in a section of the Marvdasht plain, located in Fars province. Materials and MethodsThe present study was conducted in a part of Marvdasht plain located between 35.82´41°52' to 1.07´57°52' east longitude and 35.02´48°29' to 14.72´2°30' north latitude, and 40 km north of Shiraz with an area of about 50,000 hectares. After determining the study area boundaries, the positions of 200 sampling points were determined using the R software and the conditioned Latin hypercube sampling method. In other words, for soil feature modeling, 200 samples were taken from two depths of zero to 30 and 30 to 60 centimeters in the study area. Then, the samples were transferred to the laboratory, dried, and passed through a 2 mm sieve. Finally, the soil texture components were measured by the hydrometer method. The environmental variables used in this study are a wide range of representatives of soil-forming factors that were prepared as much as possible from sources with minimum cost and high accessibility. In total, 75 environmental variables were prepared, and the raster format related to all environmental variables, including 39 elevation and altitude variables and 36 remote sensing measurement variables, was extracted. Finally, the factor-tuning inflation variance and Boruta algorithm were used to select the optimal variables. ResultsThe minimum amount of clay was measured at 10.21% and 10.45%, respectively, and the maximum amount was 32.65% and 36.35% at the surface and subsurface depths. The average amount of clay in all samples was 37.91% and 35.61%. The average amount of sand was measured at 25.65% and 26.02% at the surface and subsurface depths, respectively. The maximum amount of sand was observed in the northern and higher parts of the study area and was equal to 54.68% and the minimum amount was predicted in the low-lying areas of the study area. Low-lying areas and sedimentary plains in the central part of the study area contained high amounts of silt. Four depth variables valley depths (VD), texture (TE), topographic wetness index (TWI), and clay index (CI) related to geomorphometric parameters and the normalized difference vegetation index (NDVI) variable related to remote sensing indices were selected as optimal variables. The RF model with R2 of 54.0% and 36.0% for predicting sand, 48.0% and 64.0% for predicting silt, and 52.0% and 49.0% for predicting clay at both surface and subsurface depths performed better than the SVR and Co-Kriging models. The most effective variable in predicting the spatial distribution of soil particles was VD with relative importance of 60% and 65% for predicting sand at the surface and subsurface depths, 70% for predicting silt at the surface depth, and 70% and 65% for predicting clay at both surface and subsurface depths, respectively. Only TE and TWI variables were more important than VD for predicting silt at subsurface depth. These results show that topographic variables are effective in the spatial variation of soil particles. Unlike clay, the highest amount of sand in both depths was observed in the northern part and the highest part of the study area, and the lowest amount was predicted in the low-lying areas of the study area. ConclusionIn general, with the aim of this research, maps of the spatial distribution of soil texture components were prepared at both surface and subsurface depths using machine learning and geostatistical approaches along with environmental covariates in a part of Marvdasht plain. Among the selected environmental covariates, topographic attributes, especially the valley depth (VD), had the highest effect in justifying the spatial prediction of soil texture components. Also, the results of comparing the performance of machine learning models supported the higher efficiency of the RF model than other models. Therefore, the approach used in this study to prepare a map of soil texture components can be useful as a guide for mapping useful soil features in areas with similar climatic and topographic conditions.

Snapshots derived from the World Soil Information Service (WoSIS) are served freely to the international community. These static datasets provide quality-assessed and standardised soil profile data that can be used to support digital soil mapping and environmental applications at broad scale levels. Since the release of the preceding snapshot in 2019, new ETL (Extract, Load, Transform) procedures for screening, ingesting and standardising disparate source data have been developed. In conjunction with this, the WoSIS data model was overhauled making it compatible with the ISO 28258 and Observations and Measurements (O&M) domain models. Additional procedures for querying, serving, and downloading the publicly available standardised data have been implemented using open software (e.g. GraphQL API). Following up on a short discussion of these methodological developments we discuss the structure and content of the "WoSIS 2023-snapshot". A range of new soil datasets was shared with us, registered in the ISRIC World Data Centre for Soils (WDC-Soils) data repository, and subsequently processed in accordance with the licences specified by the data providers. An important effort has been the processing of forest soil data collated in the framework of the EU-HoliSoils project. We paid special attention to the standardisation of soil property definitions, description of the soil analytical procedures, and standardisation of the units of measurement. The "2023 snapshot" considers the following soil chemical properties (total carbon, organic carbon, inorganic carbon (total carbonate equivalent), total nitrogen, phosphorus (extractable-P, total-P, and P-retention), soil pH, cation exchange capacity, and electrical conductivity) and physical properties (soil texture (sand, silt, and clay), bulk density, coarse fragments, and water retention), grouped according to analytical procedures that are operationally comparable. Method options are defined for each analytical procedure (e.g. pH measured in water, KCl or CaCl2 solution, molarity of the solution, and soil/solution ratio). For each profile we also provide the original soil classification (i.e. FAO, WRB and USDA system with their version) and pedological horizon designations as far as these have been specified in the source databases. Three measures for "fitness-for-intended-use" are provided to facilitate informed data use: a) positional uncertainty of the profile's site location, b) possible uncertainty associated with the operationally defined analytical procedures, and c) date of sampling. The most recent (i.e. dynamic) dataset, called wosis_latest , is freely accessible via various webservices. To permit consistent referencing and citation we also provide a static snapshot (in casu, December 2023). This snapshot comprises quality-assessed and standardised data for 228 k geo-referenced profiles. The data come from 174 countries and represent more than 900 k soil layers (or horizons) and over 6 million records. The number of measurements for each soil property vary (greatly) between profiles and with depth, this generally depending on the objectives of the initial soil sampling programmes. In the coming years, we aim to gradually fill gaps in the geographic distribution of the profiles, as well as in the soil observations themselves, this subject to the sharing of a wider selection of "public" soil data by prospective data contributors. The WoSIS 2023-snapshot is archived and freely available at https://doi.org/10.17027/isric-wdcsoils-20231130 (Calisto et al., 2023). [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]

A WebGis tool called GoProsit has been developed to support winegrowers in planning a new sustainable vineyard and in the identification of high-quality terroir in Tuscany, Central Italy, by providing various information on soils, climate, hydrological risks, and fertilization. GoProsit, hosted by the web platform GEAPP, is a free, user-friendly, and interactive Geographic Decision Support System (GDSS). Soil data behind the WebGis tool has a 1 ha resolution, achieved by processing the legacy vector-type soil database of the Tuscany Region with the DSMART (Disaggregation and Harmonization of Soil Map Units Through Resampled Classification Trees as supervised classification) algorithm, which disaggregated the map to 297,023 vineyard grid cells. Each grid cell holds climatic and pedologic information, along with physical and chemical features for each horizon of the most probable soil. GoProsit also provides soil maps in image format obtained by georeferencing about 50 historical soil maps (1969–2012). Finally, GoProsit runs and returns the outputs of six models: (a) carbon footprint, (b) potential erosion and maximum vine row length compatible with tolerable erosion, (c) potential water stress, (d) risk of runoff/waterlogging, (e) identification of suitable rootstocks, and (f) nutritional needs before planting. Statistics of the main model results for the investigated area are reported. This promising tool will soon be usable for the whole Italian territory; however, its potential makes it suitable for use in any wine-growing district. [ABSTRACT FROM AUTHOR]