This study is aimed at detailed statistical and geostatistical investigation of lead and zinc concentration in an old mining area located in the eastern part of the Upper Silesian Industrial Region. This area is rich in lead and zinc ores whose intense extraction dates back to the tenth century ad. The complexity of the area results from historical and current mining activities, as well as from a variety of different types of land management and complex geological conditions. Almost 1,000 collected soil cores were divided into two subsets: those collected at the depth of up to 20 cm and the those collected at the depth from 40 to 60 cm. Extensive analyses considered geological substrata in terms of spatial variability and spatial distributions, the type of land management, geoaccumulation indexes and enrichment factors. Lead and zinc concentration was several times higher on depths ranging from 40 to 60 cm beneath the soil surface than in the 20-cm topsoil. The results showed that clearer spatial dependence was observed for deeper soil layers then for the topmost ones, especially near mines where anthropogenic factors predominated over lithogenic ones. Weak spatial dependence was accompanied by high values of the geoaccumulation index. The lowest concentrations of pollution with geoaccumulation index below 0 and enrichment factor up to 5 observed in the forest were caused by low anthropogenic pressure and the presence of sandy soils, less capable of accumulating heavy metals.

In soil magnetometry, two types of measurements are usually performed. The first type is measurements performed on the soil surface, frequently using an MS2D sensor. The second type includes measurements of magnetic susceptibility carried out in the soil profile, usually to a depth of about 30 cm. Up to now, such measurement results were analyzed separately. However, it is possible and advantageous to integrate these two types of measurements. The goal of the study was to integrate measurements of magnetic susceptibility performed on the soil surface and in the soil profile. More specifically, the goal was to obtain 3-dimensional spatial distributions of magnetic susceptibility of the topsoil horizon. Results show that it is possible to effectively integrate measurements of magnetic susceptibility performed on the soil surface and in the soil profile. Moreover, the 3-dimensional spatial distribution that is obtained shows the magnetic susceptibility of the top 20 cm of soil, which includes the soil horizons where most of the heavy metals are accumulated. The analysis of such a spatial distribution can be very helpful in delineating areas where the heightened magnetic susceptibility is a result of the influence of anthropogenic pollution from those areas where it results from lithogenic origin. It is possible to investigate where the volumes of soil with heightened magnetic susceptibility are located in the soil profile and in this way investigate which characteristic type of soil profile it is.

Directional variograms, along the soil profile, can be useful and precise tool that can be used to increase the precision of the assessment of soil pollution. The detail analysis of spatial variability in the soil profile can be also an important part of the standardization of soil magnetometry as a screening method for an assessment of soil pollution related to the dust deposition. The goal of this study was to investigate the correlation between basic parameters of spatial correlations of magnetic susceptibility in the soil profile, such as a range of correlation and a sill, and selected magnetometric indicators of soil pollution. Magnetic indicators were an area under the curve of magnetic susceptibility versus a depth in the soil profile, values of magnetic susceptibility at depths ranging from 1 to 10 cm, and maximum and background values of magnetic susceptibility in the soil profile. These indicators were previously analyzed in the literature. The results showed that a range of correlation of magnetic susceptibility was significantly correlated with magnetic susceptibility measured at depths 1, 2, and 3 cm. It suggests that a range of correlation is a good measure of pollutants’ dispersion in the soil profile. The sill of the variogram of magnetic susceptibility was found to be significantly correlated with the area under the curve of plot of magnetic susceptibility that is related to the soil pollution. In consequence, the parameters of microscale spatial variability of magnetic susceptibility in s soil profile are important measures that take into consideration the spatial aspect of s soil pollution.