Soil provides ecosystem services, supports human health and habitation, stores carbon and regulates emissions of greenhouse gases. Unprecedented pressures on soil from degradation and urbanization are threatening agro-ecological balances and food security. It is important that we learn more about soil to sustainably manage and preserve it for future generations. To this end, we developed and analyzed a global soil visible–near infrared (vis–NIR) spectral library. It is currently the largest and most diverse database of its kind. We show that the information encoded in the spectra can describe soil composition and be associated to land cover and its global geographic distribution, which acts as a surrogate for global climate variability. We also show the usefulness of the global spectra for predicting soil attributes such as soil organic and inorganic carbon, clay, silt, sand and iron contents, cation exchange capacity, and pH. Using wavelets to treat the spectra, which were recorded in different laboratories using different spectrometers and methods, helped to improve the spectroscopic modelling. We found that modelling a diverse set of spectra with a machine learning algorithm can find the local relationships in the data to produce accurate predictions. The spectroscopic models that we derived are parsimonious and robust, and using them we derived a harmonized global soil attribute dataset, which might serve to facilitate research on soil at the global scale. This spectroscopic approach should help to deal with the shortage of data on soil to better understand it and to meet the growing demand for information to assess and monitor soil at scales ranging from regional to global. We hope that this work might reinvigorate our community’s discussion towards larger, more coordinated collaborations and encourage other contributions. We also hope that use of the database will deepen our understanding of soil so that we might sustainably manage it and push the research outcomes of the soil, earth and environmental sciences towards applications that we have not yet dreamed of

Soil spatial information has traditionally been presented as polygon maps at coarse scales. Solving global and local issues, including food security, water regulation, land degradation, and climate change requires higher quality, more consistent and detailed soil information. Accurate prediction of soil variation over large and complex areas with limited samples remains a challenge, which is especially significant for China due to its vast land area which contains the most diverse soil landscapes in the world. Here, we integrated predictive soil mapping paradigm with adaptive depth function fitting, state-of-the-art ensemble machine learning and high-resolution soil-forming environment characterization in a high-performance parallel computing environment to generate 90-m resolution national gridded maps of nine soil properties (pH, organic carbon, nitrogen, phosphorus, potassium, cation exchange capacity, bulk density, coarse fragments, and thickness) at multiple depths across China. This was based on approximately 5000 representative soil profiles collected in a recent national soil survey and a suite of detailed covariates to characterize soil-forming environments. The predictive accuracy ranged from very good to moderate (Model Efficiency Coefficients from 0.71 to 0.36) at 0-5 cm. The predictive accuracy for most soil properties declined with depth. Compared with previous soil maps, we achieved significantly more detailed and accurate predictions which could well represent soil variations across the territory and are a significant contribution to the GlobalSoilMap.net project. The relative importance of soil-forming factors in the predictions varied by specific soil property and depth, suggesting the complexity and non-stationarity of comprehensive multi-factor interactions in the process of soil development.

Soil organic carbon (SOC) changes under future climate warming are difficult to quantify in situ. Here we apply an innovative approach combining space-for-time substitution with meta-analysis to SOC measurements in 113,013 soil profiles across the globe to estimate the effect of future climate warming on steady-state SOC stocks. We find that SOC stock will reduce by 6.0 ± 1.6% (mean±95% confidence interval), 4.8 ± 2.3% and 1.3 ± 4.0% at 0–0.3, 0.3–1 and 1–2 m soil depths, respectively, under 1 °C air warming, with additional 4.2%, 2.2% and 1.4% losses per every additional 1 °C warming, respectively. The largest proportional SOC losses occur in boreal forests. Existing SOC level is the predominant determinant of the spatial variability of SOC changes with higher percentage losses in SOC-rich soils. Our work demonstrates that warming induces more proportional SOC losses in topsoil than in subsoil, particularly from high-latitudinal SOC-rich systems.