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-5 cm, 5-15 cm, 15-30 cm, 30-60 cm, 60-100 cm and 100-200 cm. These depths are consistent with the specifications of the GlobalSoilMap.net project - GlobalSoilMaps. The digital soil attribute maps are in raster format at a resolution of 3 arc sec (~90 x 90 m pixels).
Detailed information about the Soil and Landscape Grid of Australia can be found at - SLGA- Attribute Definition: Mass fraction of carbon by weight in the < 2 mm soil material as determined by dry combustion at 900 Celsius Units: %;
- Period (temporal coverage; approximately): 1970-2021;
- Spatial resolution: 3 arc seconds (approximately 90 m);
- Total number of gridded maps for this attribute: 18;
- Number of pixels with coverage per layer: 2007M (49200 * 40800);
- Target data standard: GlobalSoilMap specifications;
- Format: Cloud Optimised GeoTIFF
Credit
We at TERN acknowledge the Traditional Owners and Custodians throughout Australia, New Zealand and all nations. We honour their profound connections to land, water, biodiversity and culture and pay our respects to their Elders past, present and emerging.
This work was jointly funded by CSIRO, Terrestrial Ecosystem Research Network (TERN) and the Australian Government through the National Collaborative Research Infrastructure Strategy (NCRIS).
We are grateful to the custodians of the soil site data in each state and territory for providing access to the soil site data, and all of the organisations listed as collaborating agencies for their significant contributions to the project and its outcomes.
Purpose
The aim is to operate an open national capability that provides access to verified, science-quality land surface dynamics data and soils information layers, plus high-end data analytics tools that integrated with other TERN observations can meet the needs of ecosystem researchers and actionable information for policy makers and natural resource managers.
Lineage
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 (see Arrouays et al., 2014). Boxplots with TOC values by biome and land cover after data cleaning and depth standardization are shown in Figure 1.
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 bilinear interpolation or aggregated to conform with a spatial resolution with grid cell of 90 m x 90 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. For the mapping we used a machine learning model called quantile regression forest.
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).
Authors
Point of Contact
Publisher
Method Documentation