This paper describes the technical development and accuracy assessment of the most recent and improved version of the SoilGrids system at 250m resolution (June 2016 update). SoilGrids provides global predictions for standard numeric soil properties (organic carbon, bulk density, Cation Exchange Capacity (CEC), pH, soil texture fractions and coarse fragments) at seven standard depths (0, 5, 15, 30, 60, 100 and 200 cm), in addition to predictions of depth to bedrock and distribution of soil classes based on the World Reference Base (WRB) and USDA classification systems (ca. 280 raster layers in total). Predictions were based on ca. 150,000 soil profiles used for training and a stack of 158 remote sensing-based soil covariates (primarily derived from MODIS land products, SRTM DEM derivatives, climatic images and global landform and lithology maps), which were used to fit an ensemble of machine learning methods—random forest and gradient boosting and/or multinomial logistic regression—as implemented in the packages , , and . The results of 10–fold cross-validation show that the ensemble models explain between 56% (coarse fragments) and 83% (pH) of variation with an overall average of 61%. Improvements in the relative accuracy considering the amount of variation explained, in comparison to the previous version of SoilGrids at 1 km spatial resolution, range from 60 to 230%. Improvements can be attributed to: (1) the use of machine learning instead of linear regression, (2) to considerable investments in preparing finer resolution covariate layers and (3) to insertion of additional soil profiles. Further development of SoilGrids could include refinement of methods to incorporate input uncertainties and derivation of posterior probability distributions (per pixel), and further automation of spatial modeling so that soil maps can be generated for potentially hundreds of soil variables. Another area of future research is the development of methods for multiscale merging of SoilGrids predictions with local and/or national gridded soil products (e.g. up to 50 m spatial resolution) so that increasingly more accurate, complete and consistent global soil information can be produced. SoilGrids are available under the Open Data Base License.

BackgroundSoils are widely recognized as a non-renewable natural resource and as biophysical carbon sinks. As such, there is a growing requirement for global soil information. Although several global soil information systems already exist, these tend to suffer from inconsistencies and limited spatial detail.Methodology/Principal FindingsWe present SoilGrids1km — a global 3D soil information system at 1 km resolution — containing spatial predictions for a selection of soil properties (at six standard depths): soil organic carbon (g kg−1), soil pH, sand, silt and clay fractions (%), bulk density (kg m−3), cation-exchange capacity (cmol+/kg), coarse fragments (%), soil organic carbon stock (t ha−1), depth to bedrock (cm), World Reference Base soil groups, and USDA Soil Taxonomy suborders. Our predictions are based on global spatial prediction models which we fitted, per soil variable, using a compilation of major international soil profile databases (ca. 110,000 soil profiles), and a selection of ca. 75 global environmental covariates representing soil forming factors. Results of regression modeling indicate that the most useful covariates for modeling soils at the global scale are climatic and biomass indices (based on MODIS images), lithology, and taxonomic mapping units derived from conventional soil survey (Harmonized World Soil Database). Prediction accuracies assessed using 5–fold cross-validation were between 23–51%.Conclusions/SignificanceSoilGrids1km provide an initial set of examples of soil spatial data for input into global models at a resolution and consistency not previously available. Some of the main limitations of the current version of SoilGrids1km are: (1) weak relationships between soil properties/classes and explanatory variables due to scale mismatches, (2) difficulty to obtain covariates that capture soil forming factors, (3) low sampling density and spatial clustering of soil profile locations. However, as the SoilGrids system is highly automated and flexible, increasingly accurate predictions can be generated as new input data become available. SoilGrids1km are available for download via http://soilgrids.org under a Creative Commons Non Commercial license.

Abstract. SoilGrids produces maps of soil properties for the entire globe at medium spatial resolution (250 m cell size) using state-of-the-art machine learning methods to generate the necessary models. It takes as inputs soil observations from about 240 000 locations worldwide and over 400 global environmental covariates describing vegetation, terrain morphology, climate, geology and hydrology. The aim of this work was the production of global maps of soil properties, with cross-validation, hyper-parameter selection and quantification of spatially explicit uncertainty, as implemented in the SoilGrids version 2.0 product incorporating state-of-the-art practices and adapting them for global digital soil mapping with legacy data. The paper presents the evaluation of the global predictions produced for soil organic carbon content, total nitrogen, coarse fragments, pH (water), cation exchange capacity, bulk density and texture fractions at six standard depths (up to 200 cm). The quantitative evaluation showed metrics in line with previous global, continental and large-region studies. The qualitative evaluation showed that coarse-scale patterns are well reproduced. The spatial uncertainty at global scale highlighted the need for more soil observations, especially in high-latitude regions.

Aim It has been suggested that on a global scale, fire activity changes along the productivity/aridity gradient following a humped relationship, i.e. the intermediate fire–productivity hypothesis. This relation should be driven by differing relative roles of the main fire drivers (weather and fuel) along the productivity gradient. However, the full intermediate fire–productivity model across all world ecosystems remains to be validated. Location The entire globe, excluding Antarctica. Methods To test the intermediate fire–productivity hypothesis, we use the world ecoregions as a spatial unit and, for each ecoregion, we compiled remotely sensed fire activity, climate, biomass and productivity information. The regression coefficient between monthly MODIS fire activity and monthly maximum temperature in each ecoregion was considered an indicator of the sensitivity of fire to high temperatures in the ecoregion. We used linear and generalized additive models to test for the linear and humped relationships. Results Fire occurs in most ecoregions. Fire activity peaked in tropical grasslands and savannas, and significantly decreased towards the extremes of the productivity gradient. Both the sensitivity of fire to high temperatures and above‐ground biomass increased monotonically with productivity. In other words, fire activity in low‐productivity ecosystems is not driven by warm periods and is limited by low biomass; in contrast, in high‐productivity ecosystems fire is more sensitive to high temperatures, and in these ecosystems, the available biomass for fires is high. Main conclusion The results support the intermediate fire–productivity model on a global scale and suggest that climatic warming may affect fire activity differently depending on the productivity of the region. Fire regimes in productive regions are vulnerable to warming (drought‐driven fire regime changes), while in low‐productivity regions fire activity is more vulnerable to fuel changes (fuel‐driven fire regime changes).

Abstract. The World Soil Information Service (WoSIS) provides quality-assessed and standardised soil profile data to support digital soil mapping and environmental applications at broadscale levels. Since the release of the first “WoSIS snapshot”, in July 2016, many new soil data were shared with us, registered in the ISRIC data repository and subsequently standardised in accordance with the licences specified by the data providers. Soil profile data managed in WoSIS were contributed by a wide range of data providers; therefore, special attention was paid to measures for soil data quality and the standardisation of soil property definitions, soil property values (and units of measurement) and soil analytical method descriptions. We presently consider the following soil chemical properties: organic carbon, total carbon, total carbonate equivalent, total nitrogen, phosphorus (extractable P, total P and P retention), soil pH, cation exchange capacity and electrical conductivity. We also consider the following physical properties: soil texture (sand, silt, and clay), bulk density, coarse fragments and water retention. Both of these sets of properties are grouped according to analytical procedures that are operationally comparable. Further, for each profile we provide the original soil classification (FAO, WRB, USDA), version and horizon designations, insofar as these have been specified in the source databases. Measures for geographical accuracy (i.e. location) of the point data, as well as a first approximation for the uncertainty associated with the operationally defined analytical methods, are presented for possible consideration in digital soil mapping and subsequent earth system modelling. The latest (dynamic) set of quality-assessed and standardised data, called “wosis_latest”, is freely accessible via an OGC-compliant WFS (web feature service). For consistent referencing, we also provide time-specific static “snapshots”. The present snapshot (September 2019) is comprised of 196 498 geo-referenced profiles originating from 173 countries. They represent over 832 000 soil layers (or horizons) and over 5.8 million records. The actual number of observations for each property varies (greatly) between profiles and with depth, generally depending on the objectives of the initial soil sampling programmes. In the coming years, we aim to fill gradually gaps in the geographic distribution and soil property data themselves, this subject to the sharing of a wider selection of soil profile data for so far under-represented areas and properties by our existing and prospective partners. Part of this work is foreseen in conjunction within the Global Soil Information System (GloSIS) being developed by the Global Soil Partnership (GSP). The “WoSIS snapshot – September 2019” is archived and freely accessible at https://doi.org/10.17027/isric-wdcsoils.20190901 (Batjes et al., 2019).