A full description of the methods used to generate this product can be found at -
AusSoilDSM.
We employed standard Digital Soil Modelling (DSM) (McBratney et. al., 2002) methods utilizing publicly available soil observation data and publicly available environmental covariate data in an environmental correlation approach using machine learning to map the soil properties of volumetric (mm/mm) Drained Upper Limit (DUL) and Soil Lower Limit (L15) across the entire continent at 6 standard depths at 90 m pixel resolution.
We used pedotransfer functions for estimating Drained Upper Limit - 1/3 bar (DUL) and Lower Limit - 15 bar (L15) from readily available soil attribute data using data from the National Soil Site Collation (NSSC) (Searle, 2014). Soil property data was obtained using the TERN SoilDataFederator (SDF).
The spatial modelling of DUL and L15 is done at six standard depth intervals conforming to the GlobalSoilMap Specifications. (GlobalSoilMap Science Committee, 2015) ie 0-5 cm, 5-15 cm, 15-30 cm, 30-60 cm, 60-100 cm and 100-200 cm. To facilitate modelling at these standard depths the observed data set depths were harmonized to these depths using a mass preserving spline method as described Bishop (1999). A total of 20545 soil profiles were splined in this way and used as inputs to the spatial modelling.
We utilised the publicly available Terrestrial Ecosystem Research Network (TERN) raster covariate stack. It is comprised of 154 individual raster data layers.
The covariate stack was used as the independent variable data for the predictions across all grid cells and at each depth.
Fifty bootstrapped model realisations using the Cubist machine learning algorithm (Quinlan, 1992) were generated and were used to predicted DUL and L15 values (mean of the bootstrap realisations) and estimate upper and lower confidence intervals (5% and 95%) across the entire continent.
The Available Water Capacity values were calculated by subtracting the L15 values of each layer from the DUL values of each layer and the upper and lower confidence intervals were estimated by combining the variances of the upper and lower confidence intervals of L15 and DUL.
To estimate the Total Available Volumetric Water Capacity (mm) to 1 and 2 meters we summed all the AWC layer values converted to mm of water to the estimated soil depth (Australia-wide 3D digital soil property maps - Depth of Soil (3 arc second resolution) Version 2) or the designated depth of the product - which ever was shallowest.
All processing for the generation of these products was undertaken using the R programming language (R Core Team, 2020).