Soil and Landscape Grid National Soil Attribute Maps - Total Phosphorus (3" resolution) - Release 2 

This is Version 2 of the Australian Total Soil Phosphorus product of the Soil and Landscape Grid of Australia.

It supersedes the Release 1 product that can be found at https://doi.org/10.25919/7j78-md43

The map gives a modelled estimate of the spatial distribution of total phosphorus 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 (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 3 arc sec (~90 x 90 m pixels). An additional measure of model reliability is through assessment of model extrapolation risk. This measure provides users a spatial depiction where model estimates are made within the domain of the observed data or not.

Detailed information about the Soil and Landscape Grid of Australia can be found at - SLGA

• Attribute Definition: Total soil phosphorus;
• Units: % (percentage of fine soil mass);
• Period (temporal coverage; approximately): 1950-2021;
• Spatial resolution: 3 arc seconds (approx 90m);
• Total number of gridded maps for this attribute: 24;
• 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;

 

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. 

The map gives a modelled estimate of the spatial distribution of total phosphorus 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. 

A full description of the methods used to generate this product can be found at - https://aussoilsdsm.esoil.io/slga-version-2-products/total-soil-phosphorus

The first effort to derive national digital soil mapping of total soil phosphorus (expressed as a percentage of fine soil mass) is published and available on the CSIRO Data Access Portal among other places. The present work sort to update this mapping as part of ongoing efforts to expand and improve Australia’s national mapping and characterisation of its soil resources. Collectively these national soil mapping efforts constitute the Soil and Landscape Grid of Australia. The original work has been deemed as Version 1 (completed 2015), while the new work logically is Version 2 (completed 2023). This work has been made possible through support and funding from Australia’s National Collaborative Research Infrastructure Strategy (NCRIS) via the Terrestrial and Ecosystem Research Network.

As with the first effort, digital soil mapping is the underpinning framework for the ultimate creation of soil maps in this instance.

As with the other more recent national digital soil mapping efforts, the SoilDataFederator (Searle 2020) has been instrumental in the dynamic collation of disparate soil observational datasets from across the country. These data have been sourced mainly from each State and Territory Government departments tasked with soil survey and collection. Plus there are other data contributions from Universities and to a lessor extent individual research groups. The SoilDataFederator also taps into the larger CSIRO developed Natsoil database (CSIRO 2020) which holds the data related to research projects and field stations that CSIRO has managed.

The improvement in digital soil mapping has come about via several mechanism.

1. A huge expansion of the available library of data corresponding to each of the main soil state factors has been made possible (Searle et al. 2022). This is through acquisition of new data sets and improvement of others compared with those used for version 1.

2. Adoption of machine learning to derive empirical relationships between target variable (total soil nitrogen content) and various data related to the state factors that help determine and control soil variability across landscapes, here the Australian continent and very nearshore islands. While the adoption of ML is not an entirely new advancement, the coupling of it with additional data, and integration of it within a psedo-3D predictive framework permit an improved ability to spatially and vertically characterise soils than Version 1 did.

3. Together with a more powerful and streamlined predictive modelling approach, the quantification of uncertainties draws on the use of the UNEEC (Uncertainty Estimation based on Empirical Errors and Clustering; Shrestha and Solomatine 2006) approach instead of bootstrapping approach so that prediction interval bounds are more custom to the variations in state factor information. Bootstrapping tends to create uniform prediction interval ranges, whereas UNEEC can distinguish areas of relatively lower and higher uncertainties based on differences in soil and landscape characteristics. Therefore, for Version 2, the uncertainties are more custom and tightly defined to the environment they are quantified in.

4. An approach to understand and characterise issues of model extrapolation has been developed. This seeks to highlight areas where there is high confidence that models are going be unreliable, because these areas are outside the range of the underpinning data used in modelling. This issue is addressed via combination of data geometric and distance-based techniques.

The sequence of steps below were carried out to develop the Version 2 products:

- Prepared point and covariate data, including filtering, cleansing, and harmonisation.

- Point data intersection with covariates.

- Creation of model and test data sets.

- Ranger model hyperparameter value optimisation.

- Ranger model fitting with best hyperparameters.

- Spatialisation of ranger models.

- Uncertainty analysis with UNEEC method including rudimentary optimisation of class number size.

- Spatialisation of model uncertainties.

- Model extrapolation work with count of observation and boundary method (point data).

- Ranger model fitting of extrapolation outcomes.

- Spatialisation of model extrapolation outcomes.

- Model evaluations with both test data and against SLGA Version 1 products.

- Delivery of digital soil mapping outputs and computer code to repository.