It supersedes the Release 1 product that can be found at https://doi.org/10.4225/08/546EE212B0048
The map gives a modelled estimate of the spatial distribution of Bulk Density in soils across Australia.
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: Bulk Density of the whole soil (including coarse fragments) in mass per unit volume by a method equivalent to the core method;
- Units: g/cm3;
- Period (temporal coverage; approximately): 1950-2021;
- Spatial resolution: 3 arc seconds (approx 90Â m);
- 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.
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 map gives a modelled estimate of the spatial distribution of Bulk Density in soils across Australia.
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
Since the first version, more measured site data has been made available and retrievable via the Australian SoilDataFederator. In 2014 there were 3776 sites with measured whole soil bulk density. For the new update, 6116 sites had measured data. Because of usually strong empirical relationships between bulk density, soil texture and soil carbon, the use of pedotransfer functions (to predict bulk density from soil texture and soil carbon) was performed with the intention of increasing data density and spatial coverage of data that would ultimately improve digital soil mapping prediction skill. This added a further 15735 sites after building a spatial pedotransfer function using a dataset of 12308 cases (3939 sites with bulk density, soil carbon and soil texture data).
The basic steps of the work entailed:
- Use soil data federator to get pertinent soils observation data.
- Develop spatial pedotransfer function prediction whole soil bulk density using soil carbon and texture data.
- Compile measured and inferred whole soil bulk density data (86306 cases), then setting aside a dataset of 7500 cases for external model evaluation.
- Predictive models using random forest algorithm with 78806 data cases fitted. To account for uncertainties in pedotransfer function inferred data, Monte Carlo simulations were performed from the pedotransfer function model. Simulation was repeated 100 times.
- Predictive model uncertainties quantified using UNEEC approach (Uncertainty Estimation based on local errors and Clustering).
- Quantification of model extension limits derived using hybrid method involving multivariate convex hull analysis and count of observations.
- Digital soil maps with quantified uncertainties (5th and 95th prediction interval limits) and assessment of model extrapolation risk were produced at 90Â m resolution for the following depths: 0-5Â cm, 5-15Â cm, 15-30Â cm, 30-60Â cm, 60-100Â cm, 100-200Â cm.
All processing for the generation of these products was undertaken using the R programming language (R Core Team, 2020).