Digital soil property mapping and uncertainty estimation using soil class probability rasters

article i nfo Digital soil property mapping Weighted means Prediction interval Uncertainty estimation Probability rasters The objective of theworkpresentedinthispaperwastoinvestigatehowrastersof theprobabilitiesofoccurrence of soil classesmay beused tocreate digital soil propertymaps and mapsof theirassociateduncertainties.The ap- proach we present is formalised in an algorithm we developed called "Digital Soil Property Mapping Using Soil Class Probability Rasters" (PROPR). Thesoilclassprobabilityrasterswerederivedpreviouslyfromaspatialdisaggregationofthe1:250,000-scaleDal- rymple Shire legacy soil polygon map from central Queensland, Australia. We created digital soil property maps of soil pH 1:5 H2O and their uncertainties (as indicated by estimates of the limits of the 90% prediction interval) at six depth increments down the soil profile (0- 5c m, 5-15 cm, 15-30 cm, 30-60 cm, 60-100 cm, 100-200 cm). The calculation of the weighted mean soil pH value for each depth incre- ment at each grid cell was based on reference pH values for each soil class and used the probabilities of occur- rence at each grid cell as weights. The calculation of the prediction interval limits for each depth increment involved sampling from the triangular distribution of the soil pH of each soil class using the soil class probabilitiesat each grid cell as weights in order to identify the number of samples to draw from each distribution. The 90% prediction interval limits were then es- timated as the 5th and 95th percentiles of the distribution of samples drawn from the soil classes' triangular distributions. The maps of soil pH displayed strong spatial patterns. Soil pH generally increased with depth. Uncertainty gener- ally increased with depth. Validation on 300 randomly-selected soil profiles returned a Lin's concordance corre- lation coefficient of 0.193 at the surface increasing to 0.266 at depth. RMSE increased with depth from about 0.75 pH units at the surface to 1.15 at depth. Soil class probability rasters are useful for generating digital soil property maps and maps of the associated uncertainties. Validation left room for improvement but the quality of the results is probably strongly affected by the quality of the spatial disaggregation that produced the soil class probability rasters. The PROPR approach may be useful in situations where profile observations are limiting but where legacy soil maps are available. Generation of the soil class probability rasters to use in PROPR is a predictive exercise in itself and so is also sub- ject to uncertainty. The probability rasters likely can be derived by several methods including logistic regression and data mining; the probability rasters we used were derived via spatial disaggregation of a legacy soil polygon map. PROPR may be useful in situations where profile observations are limiting but where legacy soil maps are avail- able. The soil class probability rasters need to be produced separately. Information on the within-soil-class vari- ability of the target soil property at each depth increment must be known in order to establish the triangular distributions for the uncertainty estimation. PROPR may reduce reliance on having sufficient soil profile observations in areas wheresuch data is limiting. We usedavailableprofileobservationstoestimatethewithin-soil-classvariabilityinordertoestablishtriangulardis-