Spatial Distribution of Lead (Pb) in Soil: A Case Study in a Contaminated Area of the Czech Republic
For decades, the Příbram district in the Czech Republic has been affected by industrial and mining activities. These activities are important sources of heavy metal pollutants that are detrimental to soil quality. A recent study examined visible–near-infrared (VNIR), shortwave-infrared (SWIR) and X-ray fluorescence (XRF) spectroscopy to model soil lead (Pb) content in a selected area located in the Příbram district. Following that study, and using the same chemical and geographical data, we examined the spatial distribution of Pb content in the soil, with a combination of different traditional spatial analyses (Moran’s I, hotspot analysis, and Kriging) that were significantly validated. One of the novel points of this work is the use of the Getis–Ord hotspot analysis before the execution of a Kriging interpolation model to better emphasize clustering patterns. The results indicated that Pb was a spatially dependent soil property and through extensive in situ sampling, it was possible to generate a very accurate Kriging interpolation model. The high-Pb hotspots coincided with topographic obstacles that were modeled using topographic profiles extracted from the open-source Google Earth platform, indicating that Pb content does not always exhibit a direct relationship with topographic height as a result of runoff, due to the contribution of topographic steps. This observation provides a new perspective on the relationship between Pb content and topographic patterns.