Comparison of point estimation techniques in the spatial analysis of radium-226, radium-228 and potassium-40 in soil

Many environmental surveys require the implementation of estimation techniques to determine the spatial distribution of the variable being investigated. Traditional methods of interpolation and estimation, for example, inverse distance squared and triangulation often ignore features of the data set such as anisotropy which may have a significant impact on the quality of the estimates produced. Geostatistical techniques may offer an improved method of estimation by modelling the spatial continuity of the variable using semi-variogram analysis. The theoretical model fitted to the semi-variogram is then used in the assignation of weighting factors to the samples surrounding the location to be estimated. This paper outlines the results of a comparison between three common estimation methods, polygonal, triangulation and inverse distance squared and a geostatistical method, in the estimation of soil radionuclide activities. The geostatistical estimation method known as kriging performed best over a range of parameters used to test the performance of the methods. Kriging exhibited the best correlation between actual and estimated values, the narrowest error distribution and the lowest overall estimation error. Polygonal estimation was best at reproducing the data set distribution. Conditional bias was evident in all the methods, low values being over-estimated and high values being under-estimated.