Updating the Australian digital soil texture mapping (Part 2*): spatial modelling of merged field and lab measurements
Brendan Malone
A C and Ross Searle B
+ Author Affiliations
- Author Affiliations
A CSIRO Agriculture and Food, Clunies Ross Street, Black Mountain, ACT 2601, Australia.
B CSIRO Agriculture and Food, 306 Carmody Road, St Lucia, Qld 4067, Australia.
C Corresponding author. Email: brendan.malone@csiro.au
Soil Research 59(5) 435-451 https://doi.org/10.1071/SR20284
Submitted: 30 September 2020 Accepted: 11 February 2021 Published: 27 April 2021
Journal Compilation © CSIRO 2021 Open Access CC BY
Abstract
Malone and Searle (2021) described a new approach to convert field measured soil texture categories into quantitative estimates of the proportion of clay, silt and sand fractions. Converted data can seamlessly integrate with laboratory measured data into digital soil mapping workflow. Here, we describe updating the Australian national coverages of clay, sand and silt content. The approach, based on machine learning, predicts each soil texture fraction at 90 m grid cell resolution, at depths 0–5 cm, 5–15 cm, 15–30 cm, 30–60 cm, 60–100 cm and 100–200 cm. The approach accommodates uncertainty in converting field measurements to quantitative estimates of texture fractions. Existing methods of bootstrap resampling were exploited to predict uncertainties, which are expressed as 90% prediction intervals about the mean prediction at each grid cell. The models and the prediction uncertainties were assessed by an external validation dataset. Results were compared with Version 1 Soil and Landscape Grid of Australia (v1.SLGA) (Viscarra Rossel et al. 2015). All predictive and functional accuracy diagnostics demonstrate improvements compared with v1.SLGA. Improvements were noted for the sand and clay fraction mapping with average improvement of 3% and 2%, respectively, in the RMSE estimates. Marginal improvements were made for the silt fraction mapping, which was relatively difficult to predict. We also made comparisons with recently released World Soil Grid products (v2.WSG) and made similar conclusions. This work demonstrates the need to continually revisit and if necessary, update existing versions of digital soils maps when new methods and efficiencies evolve. This agility is a key feature of digital soil mapping. However, without a companion program of new data acquisition through strategic field campaigns, continued re-modelling of existing data does have its limits and an eventual model skill ceiling will be reached which may not meet expectations for delivery of accurate national scale digital soils information.
Keywords: soil texture, digital soil mapping, spatial resolution, compositional data analysis, isometric log-ratio, uncertainty analysis, functional accuracy, pedometric
References
Adhikari K, Mishra U, Owens PR, Libohova Z, Wills SA, Riley WJ, Hoffman FM, Smith DR (2020) Importance and strength of environmental controllers of soil organic carbon changes with scale. Geoderma 375, 114472| Importance and strength of environmental controllers of soil organic carbon changes with scale.Crossref | GoogleScholarGoogle Scholar |
Aitchison J (1986) 'The Statistical Analysis of Compositional Data.' (Chapman & Hall: London)
Arrouays D, McBratney A, Minasny B, Hempel J, Heuvelink G, Macmillan RA, Hartemink A, Lagacherie P, McKenzie N (2014) The GlobalSoilMap project specifications. In 'GlobalSoilMap: Basis of the Global Spatial Soil Information System. Proceedings of the 1st GlobalSoilMap Conference'. Eds D Arrouays, N McKenzie, J Hempel, A Richer-de-Forges, A McBratney.) pp. 9–12. (CRS Press/Balkema: The Netherlands)
Bishop TFA, McBratney AB, Laslett GM (1999) Modelling soil attribute depth functions with equal-area quadratic smoothing splines. Geoderma 91, 27–45.
| Modelling soil attribute depth functions with equal-area quadratic smoothing splines.Crossref | GoogleScholarGoogle Scholar |
Bivand RS, Pebesma E, Gomez-Rubio V (2013) ‘Applied spatial data analysis with R.’ (Springer: New York)
Bowman G, Hutka J (2002) Particle size analysis. In ‘Soil physical measurement and interpretation for land evaluation’. (Eds NJ McKenzie, K Coughlan, HP Cresswell.) pp 224–239. (CSIRO Publishing: Melbourne, Vic)
Breiman L (2001) Random Forests. Machine Learning 45, 5–32.
| Random Forests.Crossref | GoogleScholarGoogle Scholar |
Carlile P, Bui E, Moran C, Minasny B, McBratney AB (2001) ‘Estimating soil particle size distributions and percent sand, silt and clay for six texture classes using the Australian Soil Resource Information System point database.’ (CSIRO Land and Water Technical Report 29/01: Canberra)
Chen S, Mulder VL, Martin MP, Walter C, Lacoste M, Richer-de-Forges AC, Saby NPA, Loiseau T, Hu B, Arrouays D (2019) Probability mapping of soil thickness by random survival forest at a national scale. Geoderma 344, 184–194.
| Probability mapping of soil thickness by random survival forest at a national scale.Crossref | GoogleScholarGoogle Scholar |
Chen S, Mulder VL, Heuvelink GBM, Poggio L, Caubet M, Román Dobarco M, Walter C, Arrouays D (2020) Model averaging for mapping topsoil organic carbon in France. Geoderma 366, 114237
| Model averaging for mapping topsoil organic carbon in France.Crossref | GoogleScholarGoogle Scholar |
Christensen WF (2011) Filtered Kriging for Spatial Data with Heterogeneous Measurement Error Variances. Biometrics 67, 947–957.
| Filtered Kriging for Spatial Data with Heterogeneous Measurement Error Variances.Crossref | GoogleScholarGoogle Scholar | 21361891PubMed |
Coventry RJ, Fett DER (1979) A pipette and sieve method of particle‐size analysis and some observations on its efficacy. CSIRO Australia, Division of Soils, Divisional Report No. 38: Australia.
Cressie NAC (1991) ‘Statistics for Spatial Data.’ (Wiley: New York)
Czarnecki WM, Podolak ITK (2013) Machine Learning with Known Input Data Uncertainty Measure. In: Computer Information Systems and Industrial Management. CISIM 2013. Lecture Notes in Computer Science, vol 8104. (Eds K Saeed, R Chaki, A Cortesi, S Wierzchoń) pp 379–388. (Springer, Berlin, Heidelberg)
Delhomme JP (1978) Kriging in the hydrosciences. Advances in Water Resources 1, 251–266.
| Kriging in the hydrosciences.Crossref | GoogleScholarGoogle Scholar |
Efron B, Tibshirani R (1997) Improvements on Cross-Validation: The. 632+ Bootstrap Method. Journal of the American Statistical Association 92, 548–560.
| Improvements on Cross-Validation: The. 632+ Bootstrap Method.Crossref | GoogleScholarGoogle Scholar |
Egozcue JJ, Pawlowsky-Glahn V, Mateu-Figueras G, Barceló-Vidal C (2003) Isometric Logratio Transformations for Compositional Data Analysis. Mathematical Geology 35, 279–300.
| Isometric Logratio Transformations for Compositional Data Analysis.Crossref | GoogleScholarGoogle Scholar |
Gallant J, Read A, Dowling T (2012) Building the national one-second digital elevation model for Australia. In ‘Water Information Research and Development Alliance: Science Symposium Proceedings’, Melbourne, 1–5 August 2011. (Eds J Sims, L Merrin, R Ackland, N Herron) pp. 250–245. (CSIROL Melbourne)
Gee GW, Bauder JW (1986) Particle-size analysis. In ‘Methods of soil analysis. Part 1. 2nd ed’. (Ed. A Klute.) pp. 382–411. (ASA and SSSA: Madison, WI)
Grundy MJ, Rossel RAV, Searle RD, Wilson PL, Chen C, Gregory LJ (2015) Soil and Landscape Grid of Australia. Soil Research 53, 835–844.
| Soil and Landscape Grid of Australia.Crossref | GoogleScholarGoogle Scholar |
Harwood T, Ferrier S, Harman I, Ota N, Perry J, Williams K (2014) ‘Gridded continental climate variables for Australia.’ (CSIRO Land and Water: Canberra)
Hijmans RJ (2019) ‘raster: Geographic Data Analysis and Modeling.’ (R package version 2.9–5: https://CRAN.R-project.org/package=raster)
Johnston RM, Barry SJ, Bleys E, Bui EN, Moran CJ, Simon DAP, Carlile P, McKenzie NJ, Henderson BL, Chapman G, Imhoff M, Maschmedt D, Howe D, Grose C, Schoknecht N, Powell B, Grundy M (2003) ASRIS: the database. Soil Research 41, 1021–1036.
| ASRIS: the database.Crossref | GoogleScholarGoogle Scholar |
Kidd D, Searle R, Grundy M, McBratney A, Robinson N, O’Brien L, Zund P, Arrouays D, Thomas M, Padarian J, Jones E, Bennett JM, Minasny B, Holmes K, Malone BP, Liddicoat C, Meier E, Stockmann U, Wilson P, Wilford J, Payne J, Ringrose-Voase A, Slater B, Odgers N, Gray J, van Gool D, Andrews K, Harms B, Stower L, Triantafilis J (2020) Operationalising digital soil mapping – Lessons from Australia. Geoderma Regional 23, e00335
| Operationalising digital soil mapping – Lessons from Australia.Crossref | GoogleScholarGoogle Scholar |
Kuhn M, Wing J, Weston S, Williams A, Keefer C, Engelhardt A, Cooper T, Mayer Z, Kenkel B, Benesty M, Lescarbeau R, Ziem A, Scrucca L, Tang Y, Candan C, Hunt T (2019) ‘caret: Classification and Regression Training. R package version 6.0–84.’ https://CRAN.R-project.org/package=caret
Lark RM, Bishop TFA (2007) Cokriging particle size fractions of the soil. European Journal of Soil Science 58, 763–774.
| Cokriging particle size fractions of the soil.Crossref | GoogleScholarGoogle Scholar |
Malone BP, Searle R (2020) Improvements to the Australian national soil thickness map using an integrated data mining approach. Geoderma 377, 114579
| Improvements to the Australian national soil thickness map using an integrated data mining approach.Crossref | GoogleScholarGoogle Scholar |
Malone BP, Searle R (2021) Updating the Australian digital soil texture mapping: Part 1. Re-calibration of field soil texture class centroids. Soil Research
| Updating the Australian digital soil texture mapping: Part 1. Re-calibration of field soil texture class centroids.Crossref | GoogleScholarGoogle Scholar |
Malone BP, McBratney AB, Minasny B, Laslett GM (2009) Mapping continuous depth functions of soil carbon storage and available water capacity. Geoderma 154, 138–152.
| Mapping continuous depth functions of soil carbon storage and available water capacity.Crossref | GoogleScholarGoogle Scholar |
Malone BP, Minasny B, Odgers NP, McBratney AB (2014) Using model averaging to combine soil property rasters from legacy soil maps and from point data. Geoderma 232–234, 34–44.
| Using model averaging to combine soil property rasters from legacy soil maps and from point data.Crossref | GoogleScholarGoogle Scholar |
Malone BP, Minasny B, McBratney A (2017) ‘Using R for Digital Soil Mapping.’ (Springer: The Netherlands)
Malone B, Luo Z, He D, Viscarra Rossel R, Wang E (2020) Bioclimatic variables as important spatial predictors of soil hydraulic properties across Australia’s agricultural region. Geoderma Regional 23, e00344
| Bioclimatic variables as important spatial predictors of soil hydraulic properties across Australia’s agricultural region.Crossref | GoogleScholarGoogle Scholar |
Marshall TJ (1956) A plummet balance for measuring the size distribution of soil particles. The Journal of Applied Science 7, 142–147.
McBratney AB, Mendonça Santos ML, Minasny B (2003) On digital soil mapping. Geoderma 117, 3–52.
| On digital soil mapping.Crossref | GoogleScholarGoogle Scholar |
McKenzie N, Jacquier D, Ashton L, Cresswell H (2000) ‘Estimation of soil properties using the atlas of Australian soils. Technical Report 11/00.’ (CSIRO Land and Water: Canberra, ACT)
Minasny B, McBratney AB (2001) The Australian soil texture boomerang: a comparison of the Australian and USDA/FAO soil particle-size classification systems. Soil Research 39, 1443–1451.
| The Australian soil texture boomerang: a comparison of the Australian and USDA/FAO soil particle-size classification systems.Crossref | GoogleScholarGoogle Scholar |
Muzzamal M, Huang J, Nielson R, Sefton M, Triantafilis J (2018) Mapping Soil Particle-Size Fractions Using Additive Log-Ratio (ALR) and Isometric Log-Ratio (ILR) Transformations and Proximally Sensed Ancillary Data. Clays and Clay Minerals 66, 9–27.
| Mapping Soil Particle-Size Fractions Using Additive Log-Ratio (ALR) and Isometric Log-Ratio (ILR) Transformations and Proximally Sensed Ancillary Data.Crossref | GoogleScholarGoogle Scholar |
Northcote, K, Beckmann, G, Bettenay, E, Churchward, H, Van Dijk, D, Dimmock, G, Hubble, G, Isbell, R, McArthur, W, Murtha, G, Nicolls, K, Paton, T, Thompson, C, Webb, A, Wright, M (1960–1968) ‘Atlas of Australian soils, sheets 1 to 10.’ (Melbourne University Press: Melbourne)
Odeh I, Todd A, Triantafilis J (2003) Spatial Prediction of Soil Particle-Size Fractions As Compositional Data. Soil Science 168, 501–515.
| Spatial Prediction of Soil Particle-Size Fractions As Compositional Data.Crossref | GoogleScholarGoogle Scholar |
Odgers NP, McBratney AB, Minasny B (2011) Bottom-up digital soil mapping. I. Soil layer classes. Geoderma 163, 38–44.
| Bottom-up digital soil mapping. I. Soil layer classes.Crossref | GoogleScholarGoogle Scholar |
R Core Team (2018) 'R: A Language and Environment for Statistical Computing.' (R Foundation for Statistical Computing: Vienna, Austria)
Ruddell BL, Drewry DT, Nearing GS (2019) Information Theory for Model Diagnostics: Structural Error is Indicated by Trade-Off Between Functional and Predictive Performance. Water Resources Research 55, 6534–6554.
| Information Theory for Model Diagnostics: Structural Error is Indicated by Trade-Off Between Functional and Predictive Performance.Crossref | GoogleScholarGoogle Scholar |
Searle R (2014) The Australian site data collation to support the GlobalSoilMap. In ‘GlobalSoilMap: Basis of the global spatial soil information system. CRC Press’. Eds D Arrouays, AB McBratney, J Hempel, AC Richer-de-Forges.) pp. 127–133. (CRC Press: London, UK)
Searle R, Grundy MJ, McBratney AB, Gregory LJ, Wilson PL, Malone BP, Stenson M (2019) ‘Phased development of digital soil infrastructure for Australia and its contribution to global initiatives, 2019 Joint workshop for Digital Soil Mapping and GlobalSoilMap IUSS Working groups’. Santiago, Chile.
Searle R, Kidd D, Grundy M, McBratney A, Robinson N, O’Brien L, Zund P, Arrouays D, Thomas M, Padarian J, Jones E, Bennett JM, Minasny B, Holmes K, Malone BP, Liddicoat C, Meier E, Stockmann U, Wilson P, Wilford J, Payne J, Ringrose-Voase A, Slater B, Odgers N, Gray J, van Gool D, Andrews K, Harms B, Stower L, Triantafilis J (2021) Digital Soil Mapping and Assessment for Australia and Beyond: A Propitious Future. Geoderma Regional 24, e00359
| Digital Soil Mapping and Assessment for Australia and Beyond: A Propitious Future.Crossref | GoogleScholarGoogle Scholar |
Searle R, Stenson M, Wilson PL, Gregory LJ, Singh R, Malone BP (2021b) ‘Soil data, united, will never be defeated – The SoilDataFederator, Joint Australian and New Zealand Soil Science Societies Conference’. Cairns, QLD.
Shrestha DL, Solomatine DP (2006) Machine learning approaches for estimation of prediction interval for the model output. Neural Networks 19, 225–235.
| Machine learning approaches for estimation of prediction interval for the model output.Crossref | GoogleScholarGoogle Scholar | 16530384PubMed |
Solomatine DP, Shrestha DL (2009) A novel method to estimate model uncertainty using machine learning techniques. Water Resources Research 45,
| A novel method to estimate model uncertainty using machine learning techniques.Crossref | GoogleScholarGoogle Scholar |
Somarathna PDSN, Minasny B, Malone BP, Stockmann U, McBratney AB (2018) Accounting for the measurement error of spectroscopically inferred soil carbon data for improved precision of spatial predictions. The Science of the Total Environment 631–632, 377–389.
| Accounting for the measurement error of spectroscopically inferred soil carbon data for improved precision of spatial predictions.Crossref | GoogleScholarGoogle Scholar | 29525716PubMed |
Stockmann U, Austin J, Gallant J, Cocks B, Glover M, Thomas M, Verburg K (2020) ‘Macquarie-Bogan floodplain Plant Available Water Capacity prediction case study.’ (CSIRO Technical Report: Canberra, Australia)
Taylor J, Minasny B (2006) A protocol for converting qualitative point soil pit survey data into continuous soil property maps. Australian Journal of Soil Research 44, 543–550.
| A protocol for converting qualitative point soil pit survey data into continuous soil property maps.Crossref | GoogleScholarGoogle Scholar |
van den Boogaart KG, Tolosana-Delgado R, Bren M (2018) ‘compositions: Compositional Data Analysis.’ (R package version 1.40–2: https://CRAN.R-project.org/package=compositions)
Viscarra Rossel RA (2011) Fine-resolution multiscale mapping of clay minerals in Australian soils measured with near infrared spectra. Journal of Geophysical Research. Earth Surface 116,
| Fine-resolution multiscale mapping of clay minerals in Australian soils measured with near infrared spectra.Crossref | GoogleScholarGoogle Scholar |
Viscarra Rossel RA, Chen C, Grundy MJ, Searle R, Clifford D, Campbell PH (2015) The Australian three-dimensional soil grid: Australia’s contribution to the GlobalSoilMap project. Soil Research 53, 845–864.
| The Australian three-dimensional soil grid: Australia’s contribution to the GlobalSoilMap project.Crossref | GoogleScholarGoogle Scholar |
Wadoux AMJC, Samuel-Rosa A, Poggio L, Mulder VL (2020) A note on knowledge discovery and machine learning in digital soil mapping. European Journal of Soil Science 71, 133–136.
| A note on knowledge discovery and machine learning in digital soil mapping.Crossref | GoogleScholarGoogle Scholar |
Wiesmeier M, Urbanski L, Hobley E, Lang B, von Lützow M, Marin-Spiotta E, van Wesemael B, Rabot E, Ließ M, Garcia-Franco N, Wollschläger U, Vogel H-J, Kögel-Knabner I (2019) Soil organic carbon storage as a key function of soils - A review of drivers and indicators at various scales. Geoderma 333, 149–162.
| Soil organic carbon storage as a key function of soils - A review of drivers and indicators at various scales.Crossref | GoogleScholarGoogle Scholar |
Wilford J (2012) A weathering intensity index for the Australian continent using airborne gamma-ray spectrometry and digital terrain analysis. Geoderma 183–184, 124–142.
| A weathering intensity index for the Australian continent using airborne gamma-ray spectrometry and digital terrain analysis.Crossref | GoogleScholarGoogle Scholar |
Wright MN, Ziegler A (2017) ranger: A Fast Implementation of Random Forests for High Dimensional Data in C++ and R. Journal of Statistical Software 77, 17
| ranger: A Fast Implementation of Random Forests for High Dimensional Data in C++ and R.Crossref | GoogleScholarGoogle Scholar |
