Production of high quality interpolation maps of heavy metals is important for risk assessment of environmental pollution. In this paper, the spatial correlation characteristics information obtained from Moran’s I analysis was used to supplement the traditional geostatistics. According to Moran’s I analysis, four characteristics distances were obtained and used as the active lag distance to calculate the semivariance. Validation of the optimality of semivariance demonstrated that using the two distances where the Moran’s I and the standardized Moran’s I, Z(I) reached a maximum as the active lag distance can improve the fitting accuracy of semivariance. Then, spatial interpolation was produced based on the two distances and their nested model. The comparative analysis of estimation accuracy and the measured and predicted pollution status showed that the method combining geostatistics with Moran’s I analysis was better than traditional geostatistics. Thus, Moran’s I analysis is a useful complement for geostatistics to improve the spatial interpolation accuracy of heavy metals.

This study explored the spatial pattern of heavy metals in Beijing agricultural soils using Moran’s I statistic of spatial autocorrelation. The global Moran’s I result showed that the spatial dependence of Cr, Ni, Zn, and Hg changed with different spatial weight matrixes, and they had significant and positive global spatial correlations based on distance weight. The spatial dependence of the four metals was scale-dependent on distance, but these scale effects existed within a threshold distance of 13 km, 32 km, 50 km, and 29 km, respectively for Cr, Ni, Zn, and Hg. The maximal spatial positive correlation range was 57 km, 70 km, 57 km, and 55 km for Cr, Ni, Zn, and Hg, respectively and these were not affected by sampling density. Local spatial autocorrelation analysis detected the locations of spatial clusters and spatial outliers and revealed that the pollution of these four metals occurred in significant High-high spatial clusters, Low-high, or even High-low spatial outliers. Thus, three major areas were identified and should be receiving more attention: the first was the northeast region of Beijing, where Cr, Zn, Ni, and Hg had significant increases. The second was the southeast region of Beijing where wastewater irrigation had strongly changed the content of metals, particularly of Cr and Zn, in soils. The third area was the urban fringe around city, where Hg showed a significant increase.