Next Article in Journal
The Influence of Hydrothermal Fluids on the REY-Rich Deep-Sea Sediments in the Yupanqui Basin, Eastern South Pacific Ocean: Constraints from Bulk Sediment Geochemistry and Mineralogical Characteristics
Next Article in Special Issue
Reducing the Negative Technogenic Impact of the Mining Enterprise on the Environment through Management of the Water Balance
Previous Article in Journal
A Robust Stochastic Approach to Mineral Hyperspectral Analysis for Geometallurgy
Previous Article in Special Issue
Heavy Metals Content in the Soils of the Tatra National Park Near Lake Morskie Oko and Kasprowy Wierch—A Case Study (Tatra Mts, Central Europe)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quality of Peri-Urban Soil Developed from Ore-Bearing Carbonates: Heavy Metal Levels and Source Apportionment Assessed Using Pollution Indices

1
Institute of Earth Sciences, University of Silesia in Katowice, Bedzinska 60, 41-200 Sosnowiec, Poland
2
Department of Agricultural and Environmental Chemistry, University of Agriculture in Krakow, Aleja Adama Mickiewicza 21, 31-120 Kraków, Poland
*
Author to whom correspondence should be addressed.
Minerals 2020, 10(12), 1140; https://doi.org/10.3390/min10121140
Submission received: 20 November 2020 / Revised: 15 December 2020 / Accepted: 16 December 2020 / Published: 19 December 2020
(This article belongs to the Special Issue Elemental Concentration and Pollution in Soil, Water, and Sediment)

Abstract

:
Pollution indices are used to assess the influence of the bedrock as a natural source of heavy-metal (HM), and anthropogenic pollution from ore mining in soils developed from ore-bearing carbonates. The research was conducted in two areas differing in geological setting and type of land use in the Upper Silesia Industrial Region, Southern Poland. Physical properties such as pH, total sulfur, total carbon and total organic carbon values, as well as total Zn, Pb, and Cd contents (ICP-OES) for 39 topsoil samples were measured. Contamination factor (Cf), degree of contamination (Cdeg), pollution load index (PLI) and geoaccumulation index (Igeo), were used to determine the deterioration of topsoil due to HM pollution. The HM content exceeded geochemical background levels by 2.5–18.1 times. Very high to moderate topsoil contamination was determined. In a shallow historical mining zone, the relative influence of particular HM was found to be in the order of Pb > Cd > Zn and, in a deep mining zone, Zn > Cd > Pb. In the topsoil developed over shallow ore bodies, the HM content was mainly (60%) due to naturally occurring HM. In the area of deeply buried ore bodies, 90% of the HM load was related to anthropogenic sources. Zn, Pb and Cd vertical distributions and the patterns of topsoil pollution differ in terms of types of mined ores, mining methods and times elapsed since mining ceased. Pollution indices are an efficient tool for distinguishing soil anthropogenic pollution and geogenic contamination.

Graphical Abstract

1. Introduction

Heavy metals (HM) in the environment originate from geological, industrial, agricultural, atmospheric, and waste sources. Soil is one of the most important environmental components at risk of HM contamination as a result of anthropogenic activities. A high concentration of HMs and metalloids (Zn, Pb, Cd, As, Tl, etc.) can be found in and around active and abandoned mines or smelting plants [1,2,3] due to the emission and dispersion of pollutants into air [4,5], water [6,7], soil [8,9,10,11,12,13,14,15,16,17], plants [4,18,19,20,21,22], and fauna [23,24]. The global average concentration of HMs in soils varies for Zn (10–300 mg·kg−1), Pb (10–150 mg·kg−1) and Cd (0.06 mg·kg−1) [25]. In general, the HM concentrations in soils are increasing over time, and the highest concentrations are observed in industrial cities, due to traffic, power plants, and other industrial activities [25]. Soil contains baseline or background concentrations of HM. The HM content is determined by the composition of the parent rock material from which the soil is derived.
The effective assessment of soil HM contamination is an important global issue [10,25,26,27,28,29]. Indicators for the geochemical assessment of the soil environment include contamination factor (Cf), degree of contamination (Cdeg), pollution load index (PLI), and geoaccumulation index (Igeo). These indicators enable the estimation of environmental risk and soil degradation due to accumulation of HMs [28,29,30]. Moreover, they facilitate differentiation between the accumulation of HMs produced by natural processes and anthropogenic activities (e.g., [31]). Essential to any assessment of the degree of soil contamination is the selection of an appropriate reference value. Although this issue has already been widely discussed [26,32,33,34,35,36], there is still a lack of unambiguous methodological findings.
The aim of the current study is to assess the quality of soil developed from Zn-Pb ore-bearing carbonates, in terms of (1) the influence of the bedrock, as a natural/geogenic source of HMs, and (2) anthropogenic pollution from ore mining and processing. We attempted a source apportionment of HMs (Zn, Pb, Cd) based on the pollution indices and using the local/on-site geochemical background. The research was carried out on the topsoil layers that are usually expected to accumulate trace elements of anthropogenic origin.
The main novelty of our work consists in taking into account data linking the indices to a broader geological context and mining history which has scarcely been published. We compare two areas located in the same Zn-Pb MVT-type deposit, which differ in terms of types of mined ores, ore-body depth, mining method and the time of mining cessation. The presented approach can be applied to the analyses of environmental risk and HM source apportionment in the abandoned mining sites worldwide.

2. Study Area

Two areas located in the north-east part of the Upper Silesia Basin (Figure 1) in the Dlugoszyn and Wilkoszyn Synclines in Jaworzno City were chosen for the study. Both synclines belong to the superior synclinorial structure with an NW-NE trending axis resulting from the Alpine orogenic movements. The inclination of the rock layers varies from 4° to 15° NW in the Dlugoszyn Syncline, and from 4° to 20° SE in the North-West limb of the Wilkoszyn Syncline. These tectonic structures comprise the Middle Triassic carbonate formation, which is composed of dolomites, limestones, and marls. The carbonate rock profile is partially altered due to the epigenetic fluid flow resulting in dolomitisation [37]. The Zn-Pb mineralisation followed the dolomitisation episode, which developed the ore-bearing dolomites with Zn-Pb ore deposits [38]. The Zn-Pb ore deposits in the area belong to the stratabound type. The ore minerals (galena, sphalerite, pyrite, marcasite, and secondary nonsulfides) form bodies of varying horizontal- and vertical extension ranging from several tens of centimetres up to several tens of meters. They consist of metasomatic-, dispersed- and cavity-filling-ores, the latter including crusted-, veined-, drusy- and breccia varieties [39]. Most of the ore bodies have a typical tabular form concordant with the rock bedding. The boundaries between the ore-bearing dolomites and the host carbonate rocks are rarely defined, but more commonly form a broad transition zone into weakly mineralized dolostones.
The areas selected for the study vary in terms of the geological setting. In the Dlugoszyn area (see point 1 in Figure 1 and Figure 2), the ore-bearing dolomites occur near the surface (0–1 m), while in the Wilkoszyn syncline (see point 2 in Figure 1 and Figure 2) the ore bodies occur at a depth of 20 to 90 m. The relatively shallow occurrence of ore mineralisation favours the development of mining. There is evidence of outcrop and underground Zn-Pb ore mining and smelting from the 12th [40] to the 20th century. The shallow ore exploitation in the Dlugoszyn area was completed at the end of the 19th century, while the deep mines and processing plants in the Wilkoszyn area were abandoned in the middle of the 20th century [41,42]. There is evidence of the historical mining and industrial activity in the landscape of the town [43]. Local ore mining left its mark on natural habitats impacting air [44], water [45,46] and soil [2,3,47,48,49] quality and imposing imprints onto flora [50,51] and fauna [52,53].

3. Material and Analytical Methods

To investigate the relationship between the ore-bearing bedrock and content of HMs (Zn, Pb, and Cd) in the soils (mostly skeletal Rendzic Leptosol), 15 soil profiles were studied (Figure 3). We tried to select homogenous soil profiles. Each soil profile was divided into separate horizons, providing 54 soil samples. Only the topsoil layers (up to 30 cm depth), characterised by the highest concentrations of HM, were considered further; these are represented by 39 soil samples in Table 1. Moreover, bedrock specimens (from each study area) were sampled and analysed to evaluate the local/natural chemical pedogenic enrichment.
Prior to the analysis, the soil samples were oven-dried at 105 °C to constant weight, sieved to 2 mm through a stainless-steel sieve, and milled into a fine powder, while the bedrock samples were crushed and ground. The soil pH was determined using a 1:5 (g:mL) ratio of soil and 1 M KCl solution, with pH meter ELMETRON CP-315 m. The concentration of total carbon (TC), total organic carbon (TOC), and total sulfur (TS) were determined using an Eltra CS-530 IR-analyser with a TIC module. In order to determine the total content of Zn, Pb, and Cd, the soil and rock samples were wet digested in a closed system in a mixture of 6 cm3 of concentrated nitric acid (65% Suprapur), and 2 cm3 of hydrochloric acid (30% Suprapur). Digestion proceeded with the use of a Perkin Elmer Multiwave 3000 Microwave Digestion, in two steps according to the program of mineralisation recommended by the manufacturer—power: 1400 W, recovery time: 5 min, hold: 25 min in the first step, 10 min in the second, fan speed: 1 in the first step, 3 in the second. After mineralisation, the samples were transferred to measuring flasks (10 cm3) with a 1% solution of Suprapur nitric acid. The HM content in the prepared solution was determined using an atomic emission spectrometer ICP-OES Optima 7300 Dual View Perkin Elmer. Each sample was analysed twice. The quality control procedure was performed using internal laboratory standard for soil material. The percentage recovery for observed elements ranged from 86.5% to 100% (Supplementary Materials Table S1). The ICP-OES analyses were performed at the Department of Agricultural and Environmental Chemistry, University of Agriculture in Krakow.
The phase composition was provided by X-ray diffraction (XRD). The analyses were performed on powdered samples using a PANalytical X’Pert Pro MPD (multipurpose diffractometer) powered by a Philips PW3040/60 X-ray generator and fitted with a 1D silicon strip detector (X’Celerator). The measurements were performed using Co Kα-radiation with a wavelength of 0.1789010 nm, an acceleration voltage of 40 kV, a current of 40 mA, and with 0.02° 2θ step sizes between the angles of 5° and 70° 2θ and a 300 s measurement time per step. The data obtained were processed using HighScore+ software and the ICSD database and PDF4+ ICDD database. All XRD analyses were performed at the Institute of Earth Sciences, Faculty of Natural Sciences, University of Silesia, Sosnowiec.

4. Soil Pollution Assessment

The concentration of HMs in the soil is related either to the natural abundance if they come from parent rocks or the anthropogenic load if they come from mining- and processing activities. The distinction between the background of the site level and the pollution load is essential to determine the nature of any anomaly, especially in industrial areas such as the selected study area.
The HM contamination in the soil samples is evaluated by applying the Cf, Cdeg, PLI, and Igeo.
The Cf is employed to determine the degree of soil pollution with a particular potentially toxic element. The calculation of this factor follows Håkonson’s [55] definition (Equation (1)):
C f = C   H M   C B ,
where CHM is the HM concentration in soil samples (mg·kg−1) and CB is reference concentration of HMs in soil in the study area (mg·kg−1). According to Håkonson [55], four contamination categories can be distinguished: Cf < 1, low; 1 < Cf < 3, moderate; 3 < Cf < 6, considerable; and Cf > 6, very high contamination.
The estimation of the degree of total contamination of surface layers in a particular core or sampling site was proposed by Håkanson [55] as a sum of the Cf for each sample, or Cdeg (Equation (2)):
C d e g =   i = 1 i = n C f ,
The classification of the Cdeg in sediments is as follows [55]: Cdeg < 6, low; 6 < Cdeg < 12, moderate; 12 < Cdeg < 24, considerable; and Cdeg > 24, high.
The status of the HM pollution in the studied soils is assessed through the PLI [56]. The PLI is calculated by obtaining the n-root from the product of n- Cf s for all the HMs, where n is the number of HMs studied (Equation (3)):
P L I = ( C F H M 1   × C F H M 2 . × C F H M n ) n ,
A PLI value < 1 indicates minimal or no metal pollution, a PLI = 1 indicates that the level of pollutants is equal to the baseline content of HM, whereas PLI > 1 means that there is pollution and the value indicates the severity of the pollution. For the latter case, Zhang et al. [57] classify soils into the following five categories: 1 < PLI < 2, moderately; 2 < PLI < 3, moderately to highly; 3 < PLI < 4, highly 4 < PLI < 5, very highly; and PLI > 5, extremely polluted.
The Igeo proposed by Müller [58] is a common criterion to evaluate the HM pollution in sediments/soil. It is defined as follows (Equation (4)):
I g e o = l o g 2 ( C H M 1.5   C B ) ,
where CHM is the HM concentration in soil samples (mg·kg−1); CB is reference concentration of HMs in the soil in the study area (mg·kg−1), and factor 1.5 is incorporated in the relationship to account for possible variation in background data due to the lithogenic effect. Müller [58] distinguished seven classes of soil quality based on HM enrichment: Class 0–Igeo < 0, uncontaminated; Class 1–0 < Igeo < 1, uncontaminated to moderately; Class 2–1 < Igeo < 2, moderately; Class 3–2 < Igeo < 3, moderately to highly; Class 4–3 < Igeo < 4, highly; Class 5–4 < Igeo < 5, highly to extremely; Class 6–Igeo > 5, extremely contaminated. Igeo Class 0 suggests the lack of contamination, while the highest Class 6 shows the extreme enrichment of the HMs relative to background values.
The above reference concentration, symbolised as CB, depends on the local natural conditions, and differs widely from one geologic unit to another. The metal background depends on the composition of the parent rock from which the soil was derived [13]. Thus, use of Clark values or off-site reference methods are obviously not appropriate, especially in assessing weakly contaminated sites [34]. The on-site reference value is more specific and sensitive in the case of trace-element soil contamination [33]. Because of the bedrock of ore-bearing dolomites in the research area, we decided to use the HM content in the carbonate bedrock (Table 2) as a local reference value. In addition, three more off-site reference values (Table 2) were applied to further define the role played by different sources in the topsoil contamination detected.

5. Results and Discussion

The pH values of the soil in the Dlugoszyn area (Table 1) fluctuate between 4.26 and 7.24 (median 6.57), while in the Wilkoszyn area, the pH varies from 6.01 to 7.6 (median 6.77), indicating that the examined soil samples are extremely acidic to slightly alkaline, but mostly neutral. Generally, the pH value decreases with depth. The TOC values obtained for the soil range from 0.2% to 23.12% (median 1.72%) in the Dlugoszyn area and from 0.68% to 13.84% (median 2.04%) in the second study area, while TS varies from 0 to 0.14% (median 0.01%) and 0.01 to 0.12% (median 0.02%), in the two study areas, respectively (Table 1). Both parameters diminish with depth. The TOC levels observed are characteristic for rendzina skeletal soils.
The concentrations of Zn, Pb, and Cd in the soil differ significantly between the elevated (Dlugoszyn) and buried (Wilkoszyn) ore-bearing dolomite areas (Table 2). The Dlugoszyn soil samples contain 68–3329 mg·kg−1 d.m. of Zn, 18–1877 mg·kg−1 d.m. of Pb, and 0.31–37.58 mg·kg−1 d.m. of Cd, with weighted arithmetic means 1351 mg·kg−1 d.m., 550 mg kg−1 d.m., and 13.01 mg·kg−1 d.m., respectively. The pollution load of the same elements in the Wilkoszyn area ranges as follows: 480–3077 mg.kg−1 d.m. for Zn, 118–1035 mg·kg−1 d.m. for Pb, and 2.34–27.66 mg·kg−1 d.m. for Cd, with weighted arithmetic means 1549 mg.kg−1 d.m., 410 mg·kg−1 d.m., and 10.43 mg·kg−1 d.m., respectively. For the soil derived from ore-bearing dolomites occurring close to the surface in the Dlugoszyn area, the mean concentrations of the studied metals exceed 2.5-, 9- and 3.8-fold of the Zn, Pb, and Cd content in bedrock, respectively. For the soil developed in the Wilkoszyn area from deeper ore-bearing dolomites, the mean values calculated are 18.1- (for Zn), 5.1- (for Pb), and 8.6-fold (for Cd) higher than the baseline level (compare Table 2), indicating intensive soil contamination in both study areas.
The main crystalline phases in the topsoil are in the 79.5–84 wt.% (Table 3) represented by quartz (Figure 4). Goethite is abundant, especially in the Wilkoszyn area (11.5 wt.%). The large concentration of dolomite (3–5%), and some calcite (to 0.5 wt.%) were also detected. The K-feldspars (microcline and orthoclase), albite, plagioclases, and kaolinite occur there as the accessory minerals. The low-crystalline goethite (7.16Å in Figure 4), is accompanied by amorphous iron oxides or sulfides (marcasite, hematite—2.695 Å) and/or iron oxide/hydroxide. Additionally, the pyroxene (augite-diopside) or spinel Zn-Al-Fe family phases were documented by diffraction peak 4.675 Å.
On the basis of the trace-element contents of the sampled soil, the natural concentration and the contamination from anthropogenic sources cannot be distinguished. To assess the degree of anthropogenic influence on the soils, the selected pollution indices (Cf, Cdeg, PLI, and Igeo) were calculated.
The Cf was calculated for each soil sample relative to the bedrock value for both study areas (Table 2). The Cf of Zn, Pb, and Cd in the Dlugoszyn area is in the range of 1.0–6.0, 2.8–30.6, and 2.1–9.6, with a mean value of 2.4, 9.1, and 3.8, respectively (Table 4). This indicates considerable to very high (80% of results for Pb) and moderate to considerable (100% for Zn and 80% for Cd) topsoil contamination level (Table 5). In the Wilkoszyn area, Cf has values of 6.9–35.8, 1.1–6.1 and 1.9–22.6 for Zn, Pb, and Cd, respectively, with a mean value of 14.8, 4.1, and 7.0, respectively (Table 4). Such Cf levels point to considerable to very high contamination degree with Zn (100% of results) and Cd (90%), and moderate to considerable pollution with Pb (80%) (Table 5).
The Cdeg index determines the degree of overall contamination of a particular sample in the study areas. Based on the Cdeg (see Table 4 and Table 5), 80% of the calculated values for the Dlugoszyn samples are recognised as moderately contaminated, while 20% are very highly contaminated. In the Wilkoszyn samples, 60% indicate very high (30%) and moderate (10%) pollution of the topsoil (Table 5).
The PLI indicates deterioration of the soil quality due to metal pollution. In the current study, the PLI obtained for the Dlugoszyn samples indicates moderate (20% of results), moderate to high (40%), high (20%) and extreme (20%) topsoil pollution (Table 5), with values ranging from 1.8 to 12.1, and mean equal to 4.3 (Table 4). In the Wilkoszyn samples, the PLI values range from 1.9 to 16.5, with a mean of 7.2 (Table 4) pointing to extremely (70% of results), very high (20%), and moderately (10%) polluted topsoil (Table 5).
The pollution indices indicate the moderately to extremely high level of topsoil contamination and reveal a difference in the level of topsoil contamination between the study areas (Table 5).
The Igeo ranges from −0.6 to 2.0 for Zn, with a mean value of 0.4; 0.9 to 4.3 for Pb, with a mean value of 1.9, and 0.5 to 2.7 for Cd, with a mean value of 1.1 in the Dlugoszyn area (Table 4). The values of Igeo for HMs decrease in the sequence Pb > Cd > Zn. The results classify the area as highly polluted with Pb and moderately polluted with Zn and Cd (Table 5). In contrast, in the Wilkoszyn area, the Igeo varies widely from 1.2–4.6 for Zn, with mean value 3.0; −0.5–2.0 for Pb, with mean value 1.3 and 0.3–3.9 for Cd, with mean value 1.8 (Table 4). The HM pollution follows the order Zn > Cd > Pb. Using Igeo, we can determine the potential source of pollution. The Igeo values obtained for the elevated ore-bearing dolomites in the Dlugoszyn area indicate mainly baseline levels of pollution for Zn and Cd, as 80% Zn and 80% Cd belong to Classes 0 and 1. The only exception is pollution with Pb (see Table 5). Most of the Zn and Cd are derived from geogenic sources, and most of the Pb (80%) is derived from anthropogenic sources. As such, in the Dlugoszyn area, the deterioration in soil quality is due to Pb. In the Wilkoszyn area, with buried ore-bearing dolomites, moderate to heavy pollution by Pb and Cd is observed, as 80% Pb and 90% Cd belong to Class 2 and 3 (Table 5). Heavy to extreme pollution by Zn is identified, as 70% Zn belongs to Class 4 and 5 (Table 5). Based on these results and the fact that concentrations of HMs decrease with depth, the HMs in the topsoil of the Wilkoszyn location is of anthropogenic origin. The graphical model of HM immissions into the topsoil of the Dlugoszyn and Wilkoszyn areas is presented in the Figure 5, while the complete set of soil profiles studied is put into the Supplementary Materials (Figure S1).
When calculating pollution indices, the selection of an appropriate reference value is an important and complicated issue. Opinions about the best value to use are divided. Some researchers favour the trace-element contents in the deeper soil horizon [32,36,62] or the average values in soil from the study area [30,63] as correct. Others consider the trace-element contents in the bedrock [64,65] as best. The accuracy and sensitivity of the on-site reference value in identifying trace-element soil contamination was evaluated by Desaules [33]. The specific geological setting of both study areas forced us to calculate the selected pollution indices with reference to the bedrock.
We also take off-site reference values defined as the median for topsoil in (1) the Cracow-Silesia region [60], (2) S Poland [59] and (3) Europe [61] into consideration. Outcomes of these examination are listed in Table 6, Table 7 and Table 8. In applying the off-site reference values to topsoil in the area where shallow ore-bearing dolomites occur (the Dlugoszyn area), a significantly higher degree of contamination is observed. Additionally, anthropogenic sources are mainly responsible for the HM presence in the soil. The percentage of class distribution shows (Table 6, Table 7 and Table 8) that 80–100% of Igeo values belong to classes 2–6. The only exception refers to Pb when the S Poland median value is used. In the case of deeper situated Zn-Pb mineralization (the Wilkoszyn area), comparable results were obtained using the off-site and on-site reference values. Significant differences were noted when the reference value for European topsoil was considered. Using the latter reference value, the extremely high soil contamination (Cf, Cdeg) here is attributed mainly to human activity (Igeo; see Table 8). In our opinion, and especially in the case where the ore-bearing dolomite is near surface, use of the local reference value to assess topsoil contamination provides more reliable results than use of off-site reference values.
Diatta et al. [9] observed similar results, with significantly high HM levels in soils impacted by the Miasteczko Slaskie Zn smelter, in part of the Upper Silesia Industrial Region. The mean Zn, Pb, and Cd contents in mg.kg−1 reported are 1062.97, 781.91, and 12.32, respectively. Using the off-site reference value to Igeo, Cf, and Cdeg calculations, extremely high contamination was observed for Zn, Cd, and Pb, in ascending order [9]. Our findings differ slightly in terms of the influence of particular HMs. In our opinion, the difference between the orders of pollutant intensity acquired for individual sites resulted primarily from the local industrial land use and the kind of reference value.
The historical mining of Zn-Pb ores in the Dlugoszyn area was based on the dominant mineral, galena (Pb sulfide), which occurs in shallow ore bodies which are easily recognised. Moreover, Pb has a relatively low melting point, allowing for ease of processing [41]. As such, ore mining in the area started in the early medieval period and resulted in the most severe Pb pollution in the Dlugoszyn topsoil. The levels of Zn and Cd are derived from the metal-rich bedrock.
The very high Zn and Pb contamination observed in the Wilkoszyn area probably results from the 20th century underground Zn-Pb ore mining and processing, when galena, Zn sulfides, and nonsulfides were mined and processed in the area. As such, it seems probable that recent industrial operations are a principal cause of the high levels of Zn and Pb accumulated in the Wilkoszyn topsoil.
Generally, soil samples from arable fields and home gardens across the Upper Silesia Industrial Region are characterised by high levels of toxic HMs. This is reflected in the HM concentrations in locally cultivated vegetables, which are well above the permissible levels [11]. Both study areas are used as peri-urban agricultural lands. The very high topsoil HM levels reported in this paper are not as harmful due to the pH value (Table 1) observed in studied soil. The HMs are regarded as highly soluble and more toxic in an acidic soil environment [66,67].

6. Conclusions

The application of various pollution indices (Cf, Cdeg, PLI, and Igeo) and the local background value CB enable the detection of Zn, Pb, and Cd pollution in soils. Depending on the index used, we determined very high to considerable (according to the classification proposed by [55]) topsoil contamination. The influence intensity of particular metals differs between studied areas. In the case of shallow historical mining (the Dlugoszyn area), it follows the order Pb > Cd > Zn, while in the place of deep mining and processing (the Wilkoszyn area) it decreases in the sequence Zn > Cd > Pb. With the geo-accumulation index, it is possible to discriminate between natural and anthropogenic HMs in the soils developed over the ore-bearing formation in the two areas. Based on Igeo values, we consider that the HM content in the topsoil developed over the shallowly occurring ore bodies in dolomites (the Dlugoszyn area) is mainly (60%) connected with the natural presence of metals. In the Wilkoszyn area where ore bodies are more deeply buried, 90% of the HM load is related to anthropogenic sources. An accurate assessment of soil quality, based on HM content, is possible only with the combined use of various pollution indices.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-163X/10/12/1140/s1, Table S1. Zn, Pb, Cd content mg. kg−1 d. m. determined in the internal laboratory standard for soil material (n = 20), Figure S1. The complete set of soil profiles for D-Dlugoszyn and W-Wilkoszyn areas

Author Contributions

Conceptualization: K.S., L.T.; Methodology: K.S.; Investigation: K.S., T.C., T.H., M.O., J.Z.; Formal analysis: K.S.; Resources: T.C.; Writing—original draft: K.S.; Writing—review & editing: L.T.; Funding acquisition: L.T.; Supervision: L.T.; Data curation and Visualization: K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the University of Silesia, Institute of Earth Sciences (WNP/INOZ/2020_ZB32).

Acknowledgments

This study was undertaken in the framework of the activities of the University of Silesia in Katowice and was funded by the University of Silesia, Institute of Earth Sciences (WNP/INOZ/2020_ZB32). We wish to thank Tomasz Krzykawski from the University of Silesia in Katowice for his assistance with XRD analysis. We are grateful to Pádhraig S. Kennan from University College Dublin for valuable remarks and suggestions. Two anonymous reviewers are thanked for improving the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Norgate, T.E.; Jahanshahi, S.; Rankin, W.J. Assessing the environmental impact of metal production processes of metal production processes. J. Clean. Prod. 2007, 15, 838–848. [Google Scholar] [CrossRef]
  2. Chrastný, V.; Vaněk, A.; Teper, L.; Cabala, J.; Procházka, J.; Pechar, L.; Drahota, P.; Penížek, V.; Komárek, M.; Novák, M. Geochemical position of Pb, Zn and Cd in soils near the Olkusz mine/smelter, South Poland: Effects of land use, type of contamination and distance from pollution source. Environ. Monit. Assess. 2012, 184, 2517–2536. [Google Scholar] [CrossRef]
  3. Vaněk, A.; Chrastný, V.; Komárek, M.; Penížek, V.; Teper, L.; Cabala, J.; Drábek, O. Geochemical position of thallium in soils from a smelter-impacted area. J. Geochem. Explor. 2013, 124, 176–182. [Google Scholar] [CrossRef]
  4. Kalinovic, T.S.; Serbula, S.M.; Radojevic, A.A.; Kalinovic, J.V.; Steharnik, M.M.; Petrovic, J.V. Elder, linden and pine biomonitoring ability of pollution emitted from the copper smelter and the tailings ponds. Geoderma 2016, 262, 266–275. [Google Scholar] [CrossRef]
  5. Ghayoraneh, M.; Qishlaqi, A. Concentration, distribution and speciation of toxic metals in soils along a transect around a Zn/Pb smelter in the northwest of Iran. J. Geochem. Explor. 2017, 180, 1–14. [Google Scholar] [CrossRef]
  6. Varol, M. Assessment of heavy metal contamination in sediments of the Tigris River (Turkey) using pollution indices and multivariate statistical techniques. J. Hazard. Mater. 2011, 195, 355–364. [Google Scholar] [CrossRef] [PubMed]
  7. Ogbeibu, A.E.; Omoigberale, M.O.; Ezenwa, I.M.; Eziza, J.O.; Igwe, J.O. Using pollution load index and geoaccumulation index for the assessment of heavy metal pollution and sediment quality of the Benin River, Nigeria. Nat. Environ. 2014, 2, 1–9. [Google Scholar] [CrossRef]
  8. Loska, K.; Wiechuła, D.; Korus, I. Metal contamination of farming soils affected by industry. Environ. Int. 2004, 30, 159–165. [Google Scholar] [CrossRef]
  9. Diatta, J.B.; Chudzinska, E.; Wirth, S. Assessment of heavy metal contamination of soils impacted by a zinc smelter activity. J. Elementol. 2008, 13, 5–16. [Google Scholar]
  10. Li, Z.; Ma, Z.; van der Kuijp, T.J.; Yuan, Z.; Huang, L. A review of soil heavy metal pollution from mines in China: Pollution and health risk assessment. Sci. Total Environ. 2014, 468–469, 843–853. [Google Scholar] [CrossRef]
  11. Dziubanek, G.; Piekut, A.; Rusin, M.; Baranowska, R.; Hajok, I. Contamination of food crops grown on soils with elevated heavy metals content. Ecotoxicol. Environ. Saf. 2015, 118, 183–189. [Google Scholar] [CrossRef] [PubMed]
  12. Ettler, V. Soil contamination near non-ferrous metal smelters: A review. Appl. Geochem. 2016, 64, 56–74. [Google Scholar] [CrossRef] [Green Version]
  13. Nouri, M.; Haddioui, A. Assessment of metals contamination and ecological risk in ait Ammar abandoned iron mine soil, Morocco. Ekológia 2016, 35, 32–49. [Google Scholar] [CrossRef] [Green Version]
  14. Obiora, S.C.; Chukwu, A.; Davies, T.C. Heavy metals and health risk assessment of arable soils and food crops around Pb-Zn mining localities in Enyigba, southeastern Nigeria. J. Afr. Earth Sci. 2016, 116, 182–189. [Google Scholar] [CrossRef]
  15. Xue, S.; Shi, L.; Wu, C.; Wu, H.; Qin, Y.; Pan, W.; Hartley, W.; Cui, M. Cadmium, lead, and arsenic contamination in paddy soils of a mining area and their exposure effects on human HEPG2 and keratinocyte cell-lines. Environ. Res. 2017, 156, 23–30. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, J.; Wang, J.; Xiao, T.; Bao, Z.; Lippold, H.; Luo, X.; Yin, M.; Ren, J.; Chen, Y.; Linghu, W. Geochemical dispersal of thallium and accompanying metals in sediment profiles from a smelter-impacted area in South China. Appl. Geochem. 2018, 88, 239–246. [Google Scholar] [CrossRef]
  17. Wu, J.; Long, J.; Liu, L.; Li, J.; Liao, H.; Zhang, M.; Zhao, C.; Wu, Q. Risk assessment and source identification of toxic metals in the agricultural soil around a Pb/Zn mining and smelting area in Southwest China. Int. J. Environ. Res. Public Health 2018, 15, 1838. [Google Scholar] [CrossRef] [Green Version]
  18. Wiggenhauser, M.; Bigalke, M.; Imseng, M.; Müller, M.; Keller, A.; Murphy, K.; Kreissig, K.; Rehkämper, M.; Wilcke, W.; Frossard, E. Cadmium isotope fractionation in soil-wheat systems. Environ. Sci. Technol. 2016, 50, 9223–9231. [Google Scholar] [CrossRef] [Green Version]
  19. Dinis, L.; Savard, M.M.; Gammon, P.; Bégin, C.; Vaive, J. Influence of climatic conditions and industrial emissions on spruce tree-ring Pb isotopes analyzed at ppb concentrations in the Athabasca oil sands region. Dendrochronologia 2016, 37, 96–106. [Google Scholar] [CrossRef]
  20. Wei, X.; Lyu, S.; Yu, Y.; Wang, Z.; Liu, H.; Pan, D.; Chen, J. Phylloremediation of air pollutants: Exploiting the potential of plant leaves and leaf-associated microbes. Front. Plant Sci. 2017, 28, 1318. [Google Scholar] [CrossRef]
  21. Liu, J.; Wei, X.; Zhou, Y.; Tsang, D.C.W.; Bao, Z.; Yin, M.; Lippold, H.; Yuan, W.; Wang, J.; Feng, Y.; et al. Thallium contamination, health risk assessment and source apportionment in common vegetables. Sci. Total Environ. 2020, 703, 135547. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, J.; Jiang, Y.; Sun, J.; She, J.; Yin, M.; Fang, F.; Xiao, T.; Song, G.; Liu, J. Geochemical transfer of cadmium in river sediments near a lead-zinc smelter. Ecotoxicol. Environ. Saf. 2020, 196, 110529. [Google Scholar] [CrossRef] [PubMed]
  23. Durkalec, M.; Szkoda, J.; Kolacz, R.; Opalinski, S.; Nawrocka, A.; Zmudzki, J. Bioaccumulation of Lead, Cadmium and Mercury in Roe Deer and Wild Boars from areas with different levels of toxic metal pollution. Int. J. Environ. Res. 2015, 9, 205–212. [Google Scholar] [CrossRef]
  24. Kapusta, P.; Sobczyk, Ł. Effects of heavy metal pollution from mining and smelting on enchytraeid communities under different land management and soil conditions. Sci. Total Environ. 2015, 536, 517–526. [Google Scholar] [CrossRef]
  25. Weissmannová, H.D.; Pavlovský, J. Indices of soil contamination by heavy metals—Methodology of calculation for pollution assessment (minireview). Environ. Monit. Assess. 2017, 189, 616. [Google Scholar] [CrossRef]
  26. Barbieri, M. The importance of enrichment factor (EF) and geoaccumulation index (Igeo) to evaluate the soil contamination. J. Geol. Geophys. 2016, 5, 237. [Google Scholar] [CrossRef]
  27. Baran, A.; Wieczorek, J. Application of geochemical and ecotoxicity indices for assessment of heavy metals content in soils. Arch. Environ. Prot. 2015, 41, 54–63. [Google Scholar] [CrossRef] [Green Version]
  28. Baran, A.; Wieczorek, J.; Mazurek, R.; Urbański, K.; Klimkowicz-Pawlas, A. Potential ecological risk assessment and predicting zinc accumulation in soils. Environ. Geochem. Health 2018, 40, 435–450. [Google Scholar] [CrossRef] [Green Version]
  29. Wieczorek, J.; Baran, A.; Urbański, K.; Mazurek, R.; Klimowicz-Pawlas, A. Assessment of the pollution and ecological risk of lead and cadmium in soils. Environ. Geochem. Health 2018, 40, 2325–2342. [Google Scholar] [CrossRef]
  30. Zang, Z.; Li, Y.; Li, H.; Guo, Z.; Zhang, R. Spatiotemporal variation and pollution assessment of Pb/Zn from smelting activities in China. Int. J. Environ. Res. Public Health 2020, 17, 1968. [Google Scholar] [CrossRef] [Green Version]
  31. Mazurek, R.; Kowalska, J.B.; Gąsiorek, M.; Zadrożny, P.; Wieczorek, J. Pollution indices as comprehensive tools for evaluation of the accumulation and provenance of potentially toxic elements in soils in Ojców National Park. J. Geochem. Explor. 2019, 201, 13–30. [Google Scholar] [CrossRef]
  32. Blaser, P.; Zimmermann, S.; Luster, J.; Shotyk, W. Critical examination of trace element enrichments and depletions in soils: As, Cr, Cu, Ni, Pb, and Zn in Swiss forest soils. Sci. Total Environ. 2000, 249, 257–280. [Google Scholar] [CrossRef]
  33. Desaules, A. Critical evaluation of soil contamination assessment methods for trace metals. Sci. Total Environ. 2012, 426, 120–131. [Google Scholar] [CrossRef]
  34. Wu, J.; Teng, Y.; Lu, S.; Wang, Y.; Jiao, X. Evaluation of soil contamination indices in a mining area of Jiangxi, China. PLoS ONE 2014, 9, e112917. [Google Scholar] [CrossRef]
  35. Maanan, M.; Saddik, M.; Maanan, M.; Chaibi, M.; Assobhei, O.; Zourarah, B. Environmental and ecological risk assessment of heavy metals in sediments of Nador lagoon, Morocco. Ecol. Indic. 2015, 48, 616–626. [Google Scholar] [CrossRef]
  36. Rachwał, M.; Kardel, K.; Magiera, T.; Bens, O. Application of magnetic susceptibility in assessment of heavy metal contamination of Saxonian soil (Germany) caused by industrial dust deposition. Geoderma 2017, 295, 10–21. [Google Scholar] [CrossRef]
  37. Leach, D.L.; Viets, J.G.; Powell, J.W. Textures of ores from the Silesian-Cracow zinc-lead deposits, Poland: Clues to the ore-forming environment. Pol. Geol. Inst. Pap. 1996, 154, 37–50. [Google Scholar]
  38. Szuwarzynski, M. Ore bodies in the Silesian-Cracow Zn-Pb ore district, Poland. Pol. Geol. Inst. Pap. 1996, 154, 9–24. [Google Scholar]
  39. Górecka, E. Mineral sequence development in the Zn-Pb deposits of the Silesian-Cracow area, Poland. Pol. Geol. Inst. Pap. 1996, 154, 37–50. [Google Scholar]
  40. Molenda, D. Mining of Lead and Silver in the Silesia-Cracow Region Up to the Middle of the 16th Century; Ossolineum: Wrocław, Poland, 1963. [Google Scholar]
  41. Panek, S. The History of Lead and Zinc Ore Mining in Jaworzno-Trzebinia-Chrzanow. Unpublished work, Jaworzno. 1995. [Google Scholar]
  42. Cabala, J.; Sutkowska, K. Past exploitation and processing of Zn-Pb ore influence on the industrial soil minerals composition. Olkusz and Jaworzno district. Pr. Nauk. Inst Górnictwa Polit Wrocł Studia Mater. 2006, 117, 13–22. [Google Scholar]
  43. Sutkowska, K.; Teper, L. Anthropogenic forms of urban post-mining landscape by example of Jaworzno. Prz. Górn. 2012, 68, 70–77. [Google Scholar]
  44. Teper, E. Dust-particle migration around flotation tailings ponds: Pine needles as passive samplers. Environ. Monit. Assess. 2009, 154, 383–391. [Google Scholar] [CrossRef] [PubMed]
  45. Aleksander-Kwaterczak, U.; Helios-Rybicka, E. Contaminated sediments as a potential source of Zn, Pb, and Cd for a river system in the historical metalliferous ore mining and smelting industry area of South Poland. J. Soils Sediments 2009, 9, 13. [Google Scholar] [CrossRef]
  46. Aleksander-Kwaterczak, U.; Ciszewski, D.; Szarek-Gwiazda, E.; Kwandrans, J.; Wilk-Woźniak, E.; Waloszek, A. The influence of historical activity of the Zn-Pb ore mine in Chrzanów on the aquatic environment quality of the Matylda valley. Górnictwo Geol. 2010, 5, 21–30. [Google Scholar]
  47. Chlopecka, A.; Bacon, J.R.; Wilson, M.J.; Kay, J. Forms of Cadmium, Lead and Zinc in contaminated soils from Southwest Poland. J. Environ. Qual. 1996, 25, 69–79. [Google Scholar] [CrossRef]
  48. Cabala, J.; Teper, L. Metalliferous constituents of rhizosphere soils contaminated by Zn-Pb mining in southern Poland. Water Air Soil Pollut. 2007, 178, 351–362. [Google Scholar] [CrossRef]
  49. Sutkowska, K.; Teper, L.; Vaněk, A.; Czech, T.; Baran, A. Effect of historical zinc processing on soil: A case study in Southern Poland. In Proceedings of the 3rd World Congress on New Technologies (NewTech’17), Rome, Italy, 6–8 June 2017; p. 110. [Google Scholar] [CrossRef]
  50. Vaněk, A.; Chrastný, V.; Teper, L.; Cabala, J.; Penížek, V.; Komárek, M. Distribution of thallium and accompanying metals in tree rings of Scots pine (Pinus sylvestris L.) from a smelter-affected area. J. Geochem. Explor. 2011, 108, 73–80. [Google Scholar] [CrossRef]
  51. Pająk, M.; Halecki, W.; Gąsiorek, M. Accumulative response of Scots pine (Pinus sylvestris L.) and silver birch (Betula pendula Roth) to heavy metals enhanced by Pb-Zn ore mining and processing plants: Explicitly spatial considerations of ordinary kriging based on a GIS approach. Chemosphere 2017, 168, 851–859. [Google Scholar] [CrossRef]
  52. Augustyniak, M.; Orzechowska, H.; Kędziorski, A.; Sawczyn, T.; Doleżych, B. DNA damage in grasshoppers’ larvae—Comet assay in environmental approach. Chemosphere 2014, 96, 180–187. [Google Scholar] [CrossRef]
  53. Dmowski, K.; Rossa, M.; Kowalska, J.; Krasnodębska-Ostręga, B. Thallium in spawn, juveniles, and adult common toads (Bufo bufo) living in the vicinity of a zinc-mining complex, Poland. Environ. Monit. Assess. 2015, 187, 4141. [Google Scholar] [CrossRef] [Green Version]
  54. Kurek, S.; Sobczyński, P.; Szuwarzyński, M.; Wojnar, E. The characteristics of the Chrzanow region deposits. Pol. Geol. Inst. Pap. 1977, 45, 45–71. [Google Scholar]
  55. Håkanson, L. An ecological risk index for aquatic pollution controls sediment logical approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  56. Tomlinson, D.L.; Wilson, J.G.; Harris, C.R.; Jeffrey, D.W. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgol. Wiss. Meeresunters. 1980, 33, 566–575. [Google Scholar] [CrossRef] [Green Version]
  57. Zhang, C.; Qiao, Q.; Piper, D.; Huang, B. Assessment of heavy metal pollution from a Fe-smelting plant in urban river sediments using environmental magnetic and geochemical methods. Environ. Pollut. 2011, 159, 3057–3070. [Google Scholar] [CrossRef] [PubMed]
  58. Müller, G. Index of geo-accumulation in sediments of the Rhine River. Geo. J. 1969, 2, 108–118. [Google Scholar]
  59. Detailed Geochemical Map of Upper Silesia. Sheet: Sławków, Olkusz, Nowa Góra, Myślachowice, Chrzanów, Dabrowa Górnicza, Strzemieszyce, Jaworzno and Libiąż; Pasieczna, A., Ed.; Polish Geological Institute: Warszawa, Poland, 2008. [Google Scholar]
  60. Lis, J.; Pasieczna, A. Geochemical Atlas of Upper Silesia 1:200,000; Publications of the Polish Geological Institute: Warszawa, Poland, 1995. [Google Scholar]
  61. Salminen, R. (Ed.) Geochemical Atlas of Europe 2005 Part 1 and 2; Geol. Survey of Finland: Espoo, Finland, 2005. [Google Scholar]
  62. Sutherland, R.A.; Tolosa, C.A.; Tack, F.M.G.; Verloo, M.G. Characterization of Selected Element Concentrations and Enrichment Ratios in Background and Anthropogenically Impacted Roadside Areas. Arch. Environ. Contam. Toxicol. 2000, 38, 428–438. [Google Scholar] [CrossRef] [PubMed]
  63. Gupta, S.; Jena, V.; Matic, N.; Kapralova, V.; Solanki, J.S. Assessment of geo-accumulation index of heavy metal and source of contamination by multivariate factor analysis. Int. J. Hazard. Mater. 2014, 2, 18–22. [Google Scholar]
  64. Szuszkiewicz, M.; Łukasik, A.; Magiera, T.; Mendakiewicz, M. Combination of geo- pedo- and technogenic magnetic and geochemical signals in soil profiles Diversification and its interpretation: A new approach. Environ. Pollut. 2016, 214, 464–477. [Google Scholar] [CrossRef]
  65. Szuszkiewicz, M.; Petrovsky, E.; Łukasik, A.; Gruba, P.; Grison, H.; Szuszkiewicz, M.M. Technogenic contamination or geogenic enrichment in Regosols and Leptosols? Magnetic and geochemical imprints on topsoil horizons. Geoderma 2021, 381, 114685. [Google Scholar] [CrossRef]
  66. Kabata-Pendias, A.; Pendias, H. Biogeochemistry of Trace Elements; PWN: Warszawa, Poland, 1999. [Google Scholar]
  67. Romero-Baena, A.J.; González, I.; Galán, E. Soil pollution by mining activities in Andalusia (South Spain)—The role of mineralogy and geochemistry in three case studies. J. Soils Sediments 2018, 18, 2231–2247. [Google Scholar] [CrossRef]
Figure 1. Geological context of the study area, without Quaternary sediments (after [54]).
Figure 1. Geological context of the study area, without Quaternary sediments (after [54]).
Minerals 10 01140 g001
Figure 2. Schematic geological cross-section of the study area showing location of ore bodies, without Quaternary sediments (after [54]).
Figure 2. Schematic geological cross-section of the study area showing location of ore bodies, without Quaternary sediments (after [54]).
Minerals 10 01140 g002
Figure 3. Characteristic examples of studied soil profiles in elevated- (profile D5) and buried (profile W3) ore-bearing dolomitic areas (photographs: T. Hulok).
Figure 3. Characteristic examples of studied soil profiles in elevated- (profile D5) and buried (profile W3) ore-bearing dolomitic areas (photographs: T. Hulok).
Minerals 10 01140 g003
Figure 4. XRD patterns of topsoil from D-Dlugoszyn (red line) and W-Wilkoszyn (blue line) area. Symbol explanations: C—calcite, D—dolomite, F—feldspars, G—goethite, Q—quartz.
Figure 4. XRD patterns of topsoil from D-Dlugoszyn (red line) and W-Wilkoszyn (blue line) area. Symbol explanations: C—calcite, D—dolomite, F—feldspars, G—goethite, Q—quartz.
Minerals 10 01140 g004
Figure 5. Graphical model of anthropogenic (red arrows) and geogenic (green arrows) HM immissions into the topsoil of the D-Dlugoszyn and W-Wilkoszyn areas.
Figure 5. Graphical model of anthropogenic (red arrows) and geogenic (green arrows) HM immissions into the topsoil of the D-Dlugoszyn and W-Wilkoszyn areas.
Minerals 10 01140 g005
Table 1. Geochemical properties of the topsoil samples. D—Dlugoszyn area, W—Wilkoszyn area.
Table 1. Geochemical properties of the topsoil samples. D—Dlugoszyn area, W—Wilkoszyn area.
Sample PointCoordinatesSample
Number
Depth
[cm]
pHTC
[%]
TOC
[%]
TS
[%]
HM Content [mg·kg−1]
ZnPbCd
D150°14′19.2″ N
19°14′30″ E
a0–26.4010.089.830.113329 ± 841674 ± 3037.58 ± 3.27
b2–257.043.292.780.033319 ± 791877 ± 3932.40 ± 1.63
D250°14′58.5″ N
19°14′18.7″ E
a0–24.2620.9220.920.14701 ± 17341 ± 611.75 ± 0.56
b2–95.269.889.880.08907 ± 8353 ± 314.72 ± 0.36
c9–196.231.351.350.01482 ± 4897 ± 104.54 ± 0.46
d19–267.240.210.200.00277 ± 3055 ± 82.14 ± 0.27
D350°15′7.2″ N
19°14′25″ E
a0–36.0422.0121.220.131438 ± 66234 ± 127.46 ± 0.51
b3–116.659.897.380.041291 ± 14329 ± 0.211.29 ± 0.05
c11–227.032.061.720.011170 ± 21248 ± 0.36.99 ± 0.07
D450°15′7.1″ N
19°14′33″ E
a0–26.418.558.500.071113 ± 12275 ± 39.59 ± 0.45
b2–146.717.146.280.071037 ± 16269 ± 89.50 ± 0.13
c14–346.911.411.390.01937 ± 23227 ± 68.11 ± 0.3
D550°15′18.2″ N
19°14′25.6″ E
a0–25.5523.1223.120.17926 ± 32283 ± 1110.46 ± 0.47
b2–146.067.167.160.07797 ± 18267 ± 210.45 ± 0.04
c14–195.960.200.200.0068 ± 818 ± 20.31 ± 0.01
W150°13’36.1” N
19°16’48.3” E
a0–26.745.765.760.062975 ± 45440 ± 226.83 ± 0.23
b2–186.973.403.400.043077 ± 54456 ± 127.66 ± 0.13
W250°13’36.2” N
19°16’58.1” E
a0–26.9311.5711.080.101861 ± 4673 ± 814.45 ± 0.25
b2–177.366.245.070.051831 ± 9711 ± 611.86 ± 0.19
W350°13’28.6” N
19°16’58.7” E
a0–26.349.229.220.072455 ± 34885 ± 1317.93 ± 0.14
b2–126.677.477.470.072761 ± 341035 ± 1821.44 ± 0.34
c12–246.690.680.680.011388 ± 2658 ± 143.30 ± 0.01
W450°13’27.7” N
19°16’44.5” E
a0–26.7413.8913.840.121956 ± 3329 ± 1913.70 ± 1.06
b2–207.087.117.100.072547 ± 9502 ± 620.57 ± 0.30
c20–307.102.032.030.022771 ± 27328 ± 312.59 ± 0.31
W550°13’27.8” N
19°16’35.4” E
a0–27.053.493.490.041063 ± 9403 ± 27.42 ± 0.05
b2–187.243.653.650.041052 ± 6400 ± 147.32 ± 0.30
c18–287.600.830.820.011124 ± 48207 ± 84.39 ± 0.21
W650°13’9” N
19°16’11” E
a0–36.212.042.040.02630 ± 15 290 ± 46.28 ± 0.29
b3–116.401.651.650.02583 ± 17275 ± 54.34 ± 0.01
W750°13’8” N
19°15’56” E
a0–36.996.336.110.071063 ± 11289 ± 15.81 ± 0.14
b3–197.003.793.700.051169 ± 1307 ± 27.09 ± 0.02
c19–267.201.451.170.01639 ± 1118 ± 12.34 ± 0.01
W850°13’15” N
19°15’56” E
a0–46.473.443.440.041003 ± 3242 ± 14.38 ± 0.03
b4–206.502.782.780.031080 ± 7298 ± 16.05 ± 0.02
W950°13’15” N
19°16’10” E
a0–36.011.591.590.02480 ± 9146 ± 53.90 ± 0.04
b3–136.121.601.600.02497 ± 2169 ± 13.95 ± 0.01
W1050°13’25” N
19°15’59” E
a0–36.532.342.340.02652 ± 7207 ± 23.71 ± 0.01
b3–176.562.312.310.03727 ± 3219 ± 14.20 ± 0.02
Table 2. HM concentrations in topsoil samples and baseline CB values (mg·kg−1 dry weight).
Table 2. HM concentrations in topsoil samples and baseline CB values (mg·kg−1 dry weight).
Sample PointWeighted Arithmetic Mean
(mg·kg−1 Dry Weight)
Statistical Parameters
ZnPbCd ZnPbCd
D13320186132.82contentmin. 68180.31
D25581737.19max.3329187737.58
D312512758.62weighted arithmetic mean 135155013.01
D49822448.69median9372699.59
D56192037.73standard deviation9105359.83
coefficient of variation [%]7712383.29
W1306645427.57
W211594927.04
W310323614.14
W4258243217.45 ZnPbCd
W510543407.34contentmin.4801182.34
W65952794.87max.3077103527.66
W710142545.66weighted arithmetic mean154941010.43
W811493314.43median11023176.68
W9287872.29standard deviation8492327.53
W107142174.11coefficient of variation [%]11023176.68
baseline CB values
ZnPbCd
1552613.41carbonate bedrock for Dlugoszyn area
286811.22carbonate bedrock for Wilkoszyn area
3104441.30median for topsoil for S Poland [59]
4200842median for topsoil in Cracow-Silesia region [60]
548150.15median for topsoil in Europe [61]
Table 3. Semiquantitative mineralogical analysis of XRD results.
Table 3. Semiquantitative mineralogical analysis of XRD results.
Mineral NameD AreaW Area
[%]
Quartz84.079.5
Dolomite4.53.0
Goethite0.511.5
Microcline3.02.5
Orthoclase3.01.0
Albite2.52.5
Kaolinite2.0-
Calcite0.5-
Table 4. Assessment of HM pollution indices in studied topsoil in relation to the local carbonate bedrock in the D-Dlugoszyn and W-Wilkoszyn areas as a background value.
Table 4. Assessment of HM pollution indices in studied topsoil in relation to the local carbonate bedrock in the D-Dlugoszyn and W-Wilkoszyn areas as a background value.
Sample PointCfCdegPLIIgeo
ZnPbCd ZnPbCd
D16.030.69.646.212.12.04.32.7
D21.02.82.16.01.8−0.60.90.5
D32.34.52.59.33.00.61.60.8
D41.84.02.58.32.60.21.40.8
D51.13.32.36.72.0−0.41.20.6
W135.85.622.664.016.54.61.93.9
W213.56.15.825.47.83.22.01.9
W312.04.53.419.95.73.01.61.2
W430.15.314.349.813.24.31.83.3
W512.34.96.023.27.13.01.72.0
W66.93.44.014.44.62.21.21.4
W711.83.14.619.65.63.01.11.6
W813.44.13.621.15.83.21.41.3
W93.31.11.96.31.91.2−0.50.3
W108.32.73.414.44.22.50.81.2
In bold—strongly contaminated.
Table 5. Percentage of class distribution of topsoil pollution for indices considered. D—Dlugoszyn area, W—Wilkoszyn area.
Table 5. Percentage of class distribution of topsoil pollution for indices considered. D—Dlugoszyn area, W—Wilkoszyn area.
Pollution IndexClasses of Soil QualityStudied Area
DW
ZnPbCdZnPbCd
Cf [%]very high contamination02020901020
considerable contamination20600107070
moderate contamination80208002010
low contamination000000
Cdeg [%]very high contamination2030
considerable contamination060
moderate contamination8010
low contamination00
PLI [%]extremely high pollution2070
very high pollution020
high pollution200
moderate to high pollution400
moderate pollution2010
unpolluted00
Igeo [%]Class 6 extremely polluted000000
Class 5 highly to extremely polluted02002000
Class 4 heavy polluted0005000
Class 3 moderately to heavily polluted20020201030
Class 2 moderately polluted0600107060
Class 1 unpolluted to moderately polluted40208001010
Class 0 unpolluted40000100
Table 6. Assessment of HM pollution indices in studied topsoil in relation to the median for topsoil in the Cracow-Silesia region as a background value [60].
Table 6. Assessment of HM pollution indices in studied topsoil in relation to the median for topsoil in the Cracow-Silesia region as a background value [60].
Sample PointCfCdegPLIIgeo
ZnPbCd ZnPbCd
D131.942.325.299.532.44.44.84.1
D25.43.95.514.84.91.81.41.9
D312.06.36.624.97.93.02.12.1
D49.45.66.721.77.12.71.92.2
D55.94.65.916.55.52.01.62.0
W129.510.321.261.018.64.32.83.8
W211.111.25.427.78.82.92.91.9
W39.98.23.221.36.42.72.51.1
W424.89.813.448.114.84.02.73.2
W510.19.15.624.98.02.82.61.9
W65.76.33.715.85.11.92.11.3
W79.75.84.419.96.32.71.91.5
W811.07.53.422.06.62.92.31.2
W92.82.01.86.52.10.90.40.2
W106.94.93.215.04.72.21.71.1
In bold—strongly contaminated.
Table 7. Assessment of HM pollution indices in studied topsoil in relation to the median for topsoil in S Poland as a background value [59].
Table 7. Assessment of HM pollution indices in studied topsoil in relation to the median for topsoil in S Poland as a background value [59].
Sample PointCfCdegPLIIgeo
ZnPbCd ZnPbCd
D116.622.116.455.218.23.53.93.5
D22.82.13.68.42.70.90.51.3
D36.33.34.313.84.52.11.11.5
D44.92.94.312.24.01.71.01.5
D53.12.43.99.43.11.00.71.4
W15.413.815.334.510.53.41.93.2
W25.93.55.815.24.92.02.01.2
W34.32.15.211.53.61.81.50.5
W45.18.712.926.88.33.11.82.5
W54.83.75.313.74.51.81.71.3
W63.03.32.48.72.91.01.10.7
W75.13.02.810.93.51.81.00.9
W85.73.92.211.93.71.91.40.6
W91.41.01.13.61.2−0.1−0.5−0.4
W103.62.62.18.22.71.30.80.5
In bold—strongly contaminated.
Table 8. Assessment of HM pollution indices in studied topsoil in relation to the median for European topsoil as a background value [61].
Table 8. Assessment of HM pollution indices in studied topsoil in relation to the median for European topsoil as a background value [61].
Sample PointCfCdegPLIIgeo
ZnPbCd ZnPbCd
D169.2124.0218.8412.0123.45.56.47.2
D211.611.547.971.118.63.02.95.0
D326.118.457.4101.930.24.13.65.3
D420.516.357.994.726.83.83.45.3
D512.913.551.578.020.83.13.25.1
W163.930.3183.8278.070.95.44.36.9
W224.232.846.9103.933.44.04.55.0
W321.524.127.673.224.33.84.04.2
W453.828.8116.4199.056.55.24.36.3
W521.926.748.997.530.63.94.25.0
W612.418.632.563.519.63.03.64.4
W721.116.937.875.823.83.83.54.7
W823.922.029.675.525.04.03.94.3
W96.05.815.327.08.12.01.93.3
W1014.914.527.456.718.13.33.34.2
In bold—strongly contaminated.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sutkowska, K.; Teper, L.; Czech, T.; Hulok, T.; Olszak, M.; Zogala, J. Quality of Peri-Urban Soil Developed from Ore-Bearing Carbonates: Heavy Metal Levels and Source Apportionment Assessed Using Pollution Indices. Minerals 2020, 10, 1140. https://doi.org/10.3390/min10121140

AMA Style

Sutkowska K, Teper L, Czech T, Hulok T, Olszak M, Zogala J. Quality of Peri-Urban Soil Developed from Ore-Bearing Carbonates: Heavy Metal Levels and Source Apportionment Assessed Using Pollution Indices. Minerals. 2020; 10(12):1140. https://doi.org/10.3390/min10121140

Chicago/Turabian Style

Sutkowska, Katarzyna, Leslaw Teper, Tomasz Czech, Tomasz Hulok, Michał Olszak, and Jan Zogala. 2020. "Quality of Peri-Urban Soil Developed from Ore-Bearing Carbonates: Heavy Metal Levels and Source Apportionment Assessed Using Pollution Indices" Minerals 10, no. 12: 1140. https://doi.org/10.3390/min10121140

APA Style

Sutkowska, K., Teper, L., Czech, T., Hulok, T., Olszak, M., & Zogala, J. (2020). Quality of Peri-Urban Soil Developed from Ore-Bearing Carbonates: Heavy Metal Levels and Source Apportionment Assessed Using Pollution Indices. Minerals, 10(12), 1140. https://doi.org/10.3390/min10121140

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop