Kriging method evaluation for assessing the spatial distribution of urban soil lead contamination

Describing contaminant spatial distribution is an integral component of risk assessment. Application of geostatistical techniques for this purpose has been demonstrated previously. These techniques may provide both an estimate of the concentration at a given unsampled location, as well as the probability that the concentration at that location will exceed a critical threshold concentration. This research is a comparative study between multiple indicator kriging and kriging with the cumulative distribution function of order statistics, with both local and global variograms. The aim was to determine which of the four methods is best able to delineate between "contaminated" and "clean" soil. The four methods were validated with a subset of data values that were not used in the prediction. Method performance was assessed by calculating the root mean square error (RMSE), analysis of variance, the proportion of sites misclassified by each method as either "clean" when they were actually "contaminated" or vice versa, and the expected loss for each misclassification type. The data used for the comparison were 807 topsoil Pb concentrations from the inner-Sydney suburbs of Glebe and Camperdown, Australia. While there was very little difference between the four methods, multiple indicator kriging was found to produce the most accurate predictions for delineating "clean" from "contaminated" soil.