Geostatistics and spatial autocorrelation analysis, in combination
with geographic information system (GIS) techniques are
widely employed to quantify the spatial distribution of soil
properties to reduce uncertainties, minimizing investigation costs and to identify pollution sources (Zhang and McGrath, 2004; Liu
et al., 2006a,b; Zhang et al., 2009).