Soil and Landscape Grid National Soil Attribute Maps - Organic Carbon (1" resolution) - Release 1

Brief description

This is Version1 of the Australian Soil Organic Carbon product of the Soil and Landscape Grid of Australia at 30m resolution.

The map gives a modelled estimate of the spatial distribution of total organic carbon in soils across Australia.

It supersedes the Release 1 product that can be found at https://doi.org/10.4225/08/547523BB0801A

The Soil and Landscape Grid of Australia has produced a range of digital soil attribute products. Each product contains six digital soil attribute maps, and their upper and lower confidence limits, representing the soil attribute at six depths: 0-5cm, 5-15cm, 15-30cm, 30-60cm, 60-100cm and 100-200cm. These depths are consistent with the specifications of the GlobalSoilMap.net project (https://esoil.io/TERNLandscapes/Public/Pages/SLGA/Resources/GlobalSoilMap_specifications_december_2015_2.pdf). The digital soil attribute maps are in raster format at a resolution of 1 arc sec (~90 x 90 m pixels).

Detailed information about the Soil and Landscape Grid of Australia can be found at - https://esoil.io/TERNLandscapes/Public/Pages/SLGA/index.html

Attribute Definition: Mass fraction of carbon by weight in the < 2 mm soil material as determined by dry combustion at 900 Celcius
Units: %;
Period (temporal coverage; approximately): 1970-2021;
Spatial resolution: 1 arc seconds (approx 30m);
Total number of gridded maps for this attribute: 18;
Number of pixels with coverage per layer: 2007M (49200 * 40800);
Data license : Creative Commons Attribution 4.0 (CC BY);
Target data standard: GlobalSoilMap specifications;
Format: Cloud Optimised GeoTIFF;
Lineage: Data on total organic carbon (TOC) concentration (%) was extracted with the Soil Data Federator (https://esoil.io/TERNLandscapes/Public/Pages/SoilDataFederator/SoilDataFederatorHelp.html) managed by CSIRO. The Soil Data Federator is a web API that compiles soil data from different institutions and government agencies throughout Australia. The laboratory methods for total organic carbon included in the study are 6A1, 6A1_UC, 6B2, 6B2b, 6B3, 6B3a. We selected TOC data from the period 1970-2020 to get a compromise between representativity of current TOC concentration and spatial coverage. The data was cleaned and processed to harmonize units, exclude duplicates and potentially wrong data entries (e.g. missing upper or lower horizon depths, extreme TOC values, unknown sampling date). Additional TOC measurements from the Biome of Australian Soil Environments (BASE) contextual data (Bisset et al., 2016) were also included in the analyses. TOC concentration for BASE samples was determined by the Walkley-Black method (method 6A1). Upper limits for TOC concentration by biome and land cover classes were set according to published literature, consistent datasets (Australian national Soil Carbon Research Program (SCaRP) and BASE, and data exploration to exclude unrealistic TOC values (e.g. maximum TOC = 30% in temperate forests, maximum TOC = 14% in temperate rainfed pasture). Since TOC concentration in Australian ecosystems has been underestimated by previous SOC maps, we did not set conservative TOC upper limits, knowing that machine learning model would likely underestimate high SOC values.

The equal-area quadratic spline function were fitted to the whole collection of pre-processed TOC data, and then values extracted for the 0-5 cm, 5-15 cm, 15-30 cm, 30-60 cm, 60-100 cm, and 100-200 cm depth intervals, following GlobalSoilMap specifications (Arrouays et al., 2014}.

Covariates: We collected a set of 57 spatially exhaustive environmental covariates covering Australia and representing proxies for factors influencing SOC formation and spatial distribution: soil properties, climate, organisms/vegetation, relief and parent material/age. The covariates were reprojected to WGS84 (EPSG:4326) projection and cropped to the same spatial extent. All covariates were resampled using billinear interpolation or aggregated to conform with a spatial resolution with grid cell of 30 m x 30 m.

Mapping: The spatial distribution of soil TOC concentration is driven by the combined influence of climate, vegetation, relief and parent materials. We thus modelled TOC concentration as a function of environmental covariates representing biotic and abiotic control of TOC. The measurement of SOC and their corresponding value of environmental covariate at same measurement locations were used to fit the mapping model.

Mapping is made with Quantile regression forest, which is similar to the popular random forest algorithm for mapping. Instead of obtaining a single statistic, that is the mean prediction from the decision trees in the random forest, we report all the target values of the leaf node of the decision trees. With QRF, the prediction is thus not a single value but a cumulative distribution of the TOC prediction at each location, which can be used to compute empirical quantile estimates.

All processing for the generation of these products was undertaken using the R programming language. R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.

Code - https://github.com/AusSoilsDSM/SLGA
Observation data - https://esoil.io/TERNLandscapes/Public/Pages/SoilDataFederator/SoilDataFederator.html
Covariate rasters - https://esoil.io/TERNLandscapes/Public/Pages/SLGA/GetData-COGSDataStore.html