Natural and anthropogenic activities affect soil heavy metal pollution at different spatial scales. Quantifying the spatial variability of soil pollution and its driving forces at different scales is essential for pollution mitigation opportunities. This study applied a multivariate factorial kriging technique to investigate the spatial variability of soil heavy metal pollution and its relationship with environmental factors at multiple scales in a highly urbanized area of Guangzhou, South China. We collected 318 topsoil samples and used five types of environmental factors for the attribution analysis. By factorial kriging, we decomposed the total variance of soil pollution into a nugget effect, a short-range (3 km) variance and a long-range (12 km) variance. The distribution of patches with a high soil pollution level was scattered in the eastern and northwestern parts of the study domain at a short-range scale, while they were more clustered at a long-range scale. The correlations between the soil pollution and environmental factors were either enhanced or counteracted across the three distinct scales. The predictors of soil heavy metal pollution changed from the soil physiochemical properties to anthropogenic dominated factors with the studied scale increase. Our study results suggest that the soil physiochemical properties were a good proxy to soil pollution across the scales. Improving the soil physiochemical properties such as increasing the soil organic matter is essentially effective across scales while restoring vegetation around pollutant sources as a nature-based solution at a large scale would be beneficial for alleviating local soil pollution.
Keywords: factorial kriging; heavy metal; multiple scales; soil pollution.