Next Article in Journal
Growth and Antioxidant Responses in Iron-Biofortified Lentil under Cadmium Stress
Previous Article in Journal
Association between Pesticides in House Dust and Residential Proximity to Farmland in a Rural Region of Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Assessing the Spatial and Vertical Distributions of Potentially Toxic Elements in Soil and How the Concentrations Differ

1
Department of Soil Science and Soil Protection, Faculty of Agrobiology, Food and Natural Resources, Czech University of Life Sciences Prague, Kamýcká 129, 16500 Prague, Czech Republic
2
Department of Environmental Management, Federal University of Technology, Owerri, P.M.B. 1526, Owerri 460114, Nigeria
3
Department of Forest Protection and Entomology, Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Kamýcká 129, 16500 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 12 April 2021 / Revised: 1 June 2021 / Accepted: 3 June 2021 / Published: 31 July 2021
(This article belongs to the Section Environmental Chemistry)

Abstract

:
A healthy soil is a healthy ecosystem because humans, animals, plants, and water highly depend upon it. Soil pollution by potentially toxic elements (PTEs) is a serious concern for humankind. The study is aimed at (i) assessing the concentrations of PTEs in soils under a long-term heavily industrialized region for coal and textiles, (ii) modeling and mapping the spatial and vertical distributions of PTEs using a GIS-based ordinary kriging technique, and (iii) identifying the possible sources of these PTEs in the Jizerské Mountains (Jizera Mts.) using a positive matrix factorization (PMF) model. Four hundred and forty-two (442) soil samples were analyzed by applying the aqua regia method. To assess the PTE contents, the level of pollution, and the distribution pattern in soil, the contamination factor (CF) and the pollution load index load (PLI) were applied. ArcGIS-based ordinary kriging interpolation was used for the spatial analysis of PTEs. The results of the analysis revealed that the variation in the coefficient (CV) of PTEs in the organic soil was highest in Cr (96.36%), followed by Cu (54.94%) and Pb (49.40%). On the other hand, the mineral soil had Cu (96.88%), Cr (66.70%), and Pb (64.48%) as the highest in CV. The PTEs in both the organic soil and the mineral soil revealed a high heterogeneous variability. Though the study area lies within the “Black Triangle”, which is a historic industrial site in Central Europe, this result did not show a substantial influence of the contamination of PTEs in the area. In spite of the rate of pollution in this area being very low based on the findings, there may be a need for intermittent assessment of the soil. This helps to curtail any excessive accumulation and escalation in future. The results may serve as baseline information for pollution assessment. It might support policy-developers in sustainable farming and forestry for the health of an ecosystem towards food security, forest safety, as well as animal and human welfare.

Graphical Abstract

1. Introduction

Soil is an indispensable component of an ecosystem that directly or indirectly links and maintains the Earth’s four spheres (namely the lithosphere, biosphere, hydrosphere, and atmosphere). However, this essential potential of soil has, in recent times, been threatened by heavy metals or potentially toxic elements (PTEs). Interestingly, the chemistry of soil makes it vulnerable to high concentrations of heavy metals or PTEs. At a required concentration, most PTEs such as Cr, Cu, Fe, Mn, Zn, Ni, Mo, Co, Se, and others are essential for plants, animals, or humans [1,2]. The presence of PTEs in soil has been attracting reasonable attention because of their ecological and biological risks. Several studies in different biomes have been performed to identify the sources of PTEs in the soils [3,4,5,6,7,8,9,10,11,12,13,14,15,16].
Natural phenomena and anthropogenic activities are the two major sources that determine the concentrations of PTEs in soils [17,18,19,20,21,22,23,24]. Natural phenomena are described as the components generated from parent material, whereas anthropogenic sources primarily originate from acute human activities [17,18,25,26]. Many authors have revealed that natural sources of some PTEs (such as Pb, Cd, and Hg) have been surpassed by anthropogenic deposits into soils due to pedogenesis [17]. Industrial inputs, the combustion of fossil fuels, municipal wastewaters, and sewage sludge have been identified as anthropogenic sources of metals [11,14,27,28]. Furthermore, intensive agricultural practices have been reported to increase PTEs in soils [16,17]. In addition, automobile and vehicle emissions, road dusts, and military activities also account for increases in PTEs [4,29]. It has been estimated by some authors that agricultural practices contributed to 79.6%, 56%, and 63% of the annual concentrations of Cu, Zn, and Cd, respectively [30]. The authors further summarized that the total annual input of Pb (85%), Ni (67.5%), and Cr (43%) found in soil emanates from industrial atmospheric deposition.
The safety of plants, food, animals, and human health have been threatened by the accumulation of PTEs in soil. Toxic elements are discharged into the soil and subsequently absorbed by plants, which are consumed by animals and humans [12]. In some cases, the PTEs penetrate into surface and underground water, which are used by living organisms including humans [31,32,33].
In the Czech Republic, edible mushrooms grown in a smelting area were reported to have been contaminated by the atmospheric deposition of Pb [34]. In Germany, there has been an urgent call to address the Pb contents in plant-based foodstuffs including bread and potatoes, which are important suppliers of this metal in human meals [35]. The yearly deposition of Cr to soils in the UK was 327 tons [36]. In addition, the study reported that 126 tons out of the 327 tons were emanated from chemical fertilizers (mostly phosphate), while 83 tons originated from atmospheric deposition and 78 tons came from sewage sludge [36]. The effect of Cr is not only recorded in food crops but also in forest trees. The health of forest plants has been at risk because of exacerbated Cr content in the soil [37]. Globally, there have been reports on the effects of increased Cu, Fe, Mn, Mo, Zn, Ni, and other PTEs on soils, plants, animals, humans, and water. Therefore, the issue has become of critical concern to the government and the stakeholders, including decision makers.
The urgency of the situation demands a robust assessment with effective quantitative and qualitative analyses. An investigation of the PTEs in soil and their sources is the principal purpose for preserving and enhancing soil quality in most areas in the world. Thus, to develop reliable policies for a sustainable soil safety for an area, it is important to have good information on the soil and its contamination level. In recent years, several analyses including statistical, geostatistical, geo-accumulation index, multivariate and modeling, as well as potential ecological risk index analyses have been proposed and applied to investigate the source, degree, and spatiotemporal state of PTE pollution in the soil [3,4,5,6,12,14,24,28,29,31,32,33]. As reported by some authors, an assessment of the contents and distribution of PTEs in soil requires intensive and robust sampling to investigate the soil conditions under distinct soil types [38,39,40,41,42]. Furthermore, considering the high temporal, spatial, and vertical variability in the uppermost soil layers of a forest, a substantial number of samples need to be examined in order to thoroughly quantify the soil adequately along an extensive scale [43]. Routinely, the study of PTE content in soil has been performed following the regular laboratory chemical methods, including atomic absorption or inductively coupled plasma analysis. These methods are expensive and time-consuming and involve consecutive serial procedures with growing complications [44,45]. Thus, a systematically structured and affordable analytical method to monitor and assess the PTEs in soil on an appropriate vertical and spatial scale is necessary [46], especially when a tangible number of sampling points are considered. The flexibility and rapid accessibility of the positive matrix factorization (PMF) model in assessing soil pollutions is remarkable [6,41]. This analytical method has high functionality for the investigation of the PTEs in soil. The PMF provides a great advantage in detecting and monitoring PTEs in organic and mineral soils: it is one of the best and latest models [6,41,47,48].
The study area is located in a northern part of the Czech Republic. The area is called Jizerské hory Mountains (Jizera Mts.). The area was polluted by past accumulation of PTEs from human and natural sources. However, there have been ongoing policies and efforts by the government and the people to ameliorate the problem, yet the impact is still prevailing in the ecosystem (mainly in the soil and vegetation). This is partly because, after soil is polluted, it takes a longer period for remediation and for ecosystem recovery processes to be completed [49]. Second, a large amount of the PTEs are enriched in the acidified forest soil, and these elements are still being discharged [5,50,51]. Some authors have reported health risks from the PTEs in high-altitude mountains in Europe including the Jizera Mts. [52,53,54]. For example, in a study performed by EscartÍn and Porte [52], the authors reported that a high percentage (76%) of polycyclic aromatic hydrocarbon (PAH) metabolites were detected in trout from the Central European high Mountains lake and that this has a high health risk. There have been many studies that focused on the acidification by sulfur and nitrogen oxides in the area. Studies focusing on the spatial and vertical distributions as well as the content and hotspots of PTEs are crucial for closing the gap in sustainable pollution assessment in the area. The benefits of applying PMF to investigate PTEs in the soil is commendable [55]. This study aims (i) to assess the concentrations of PTEs in the soil under Jizera Mts. in the Liberec region of Czech Republic after long-term, heavy industrialization; (ii) to model and map the spatial and vertical distributions of the PTEs using a GIS-based ordinary kriging technique; and (iii) to identify the possible sources of these PTEs and their contamination levels in the area using a PMF model. The findings from this study may serve as a baseline for the pollution assessment of farmland and forest soil quality in the Czech Republic and in Europe. The results might support policy-developers in sustainable farming and forestry for the health of the ecosystem and for food security, forest safety, as well as animal and human welfare.

2. Materials and Methods

2.1. Study Area

The study covered about 110 square kilometers in the Jizera Mts. The height above sea level of the area ranges from 600 to 1100 m. The average yearly temperature falls between 3 and 6 °C, which is contingent upon the altitude. The annual precipitation reaches about 1500 mm at the top of the mountains. Most areas are covered by forests (Figure 1), though in some areas, the regeneration of trees has been slow after intensive forest decline in the 1980s and 1990s [5,13]. Coniferous species, namely Norway spruce (Picea abies) and the European beech (Fagus sylvatica), are key forest trees. There are also areas with pockets of peatbogs. PTE pollution in the area is considered to have been emanated from atmospheric deposition released from the coal, textile, and steel industries and from agricultural activities. Geologically, the area is characterized by principal acidic bedrocks such as granite (granodiorite) and gneiss. Haplic/Entic Podzols, Stagnosols, and Cambisols are the predominant soils [56,57,58]. In most of the area, especially in the higher altitudes, the mor form of humus dominates while the moder humus type is observable only at lower altitudes [59]. The value of the soil pH was relatively low (Table 1), thus contributing to the high acidic condition of the area.

2.2. Soil Sampling and Laboratory Analysis

At every 3 km, soil samples were collected for both organic soil and mineral soil (to the depth of 30 cm). The samples were collected in 3 replicates for each sampling point, and the average value of the sampled points was used for the analysis. The sampling points were located using a handheld GPS system, while samples were collected using either a push probe or bucket auger depending on the terrain. A total of 221 samples each were collected from organic soil (org) and mineral soil (A) across the study area. The collected soil samples were stored in well-labelled plastic bags and taken to the laboratory. The collected samples were air-dried, ground, and sieved with a mesh of size 2.0 mm to obtain a pulverized sample.

Chemical Analysis and Instrument

The presence of elements such as Cr, Cu, Pb, Mn, and Fe in the soil were extracted using the aqua regia standard method (ISO 11466:1995, 1995) to determine their pseudo-total content [60]. For quality control (QC) of the method, the standard addition technique was adopted. For example, the QC of the concentration determination was guaranteed using the SRM 2711 (Montana II soil) reference material (National Institute of Standards and Technology, Gaithersburg, MD, USA). The values achieved were consistent with the reference data. The recovery differences were generally < 10% (n = 3). The detection limits for the elements based on the applied method were as follows: Cr (0.03 mg L−1), Cu (0.015 mg L−1), Pb (0.05 mg L−1), Mn (0.05 mg L−1), and Fe (0.15 mg L−1). The presence of Fe and Mn in the soils were also investigated; their concentrations posed no threat in the area because their concentrations were far below the EU and world recommended limits.

2.3. Contamination Level Analysis for PTEs

The PTE pollution status of the study area was assessed through various contamination assessment indices, namely the contamination factor (CF) and the pollution index (PLI).

2.3.1. Contamination Factor

CF is defined as the ratio of metal content in the sample to the background value of the same metal. It is given by the following:
CF = C (metal) Sample/C (metal) background value
where C (metal) is the concentration of metal analyzed from sampled soil and where C (metal) background value is the geochemical background value (or concentration) of that metal.
It is important to state here that the baseline values used were the world values [10].

2.3.2. Pollution Load Index (PLI)

The PLI is an estimation and was first proposed by [60]. The pollution load index has been in use for the detection of pollution. It is robust and effective in the comparison of pollution levels in space and time. The PLI was calculated based on the concentration factor of each PTE by focusing on the background value in the soil, where CF is the contamination factor earlier stated (Equation (1)) and the letter ‘n’ signifies the number of metals studied. A pollution load index less than 1 indicates the optimal soil quality, and a PLI that is equal to 1 proves that only the baseline levels of contaminants are present, while a PLI greater than 1 infers the degradation of the quality of the site by [61].
The pollution load index (PLI) equation is given by the following:
PLI = n√ (CF1 × CF2 × CF3 × … × CFn)
where CF is the contamination factor derived for each metal and where n is the number of metals.

2.4. Source Apportionment Using a Positive Matrix Factorization (PMF) Model

The positive matrix factorization (PMF) model is an effective method acquired from the software EPA-PMF v 5.0, Washington DC, USA [55]. It was applied to determine the contribution of PTE sources to contamination in the study area. The mathematical method is a receptor model used in calculating the contribution of the sources to samples built on the composition or fingerprints of the sources. The PMF model apportions the collaborations of elements in soil composition by solving chemical mass balance:
C i j =   K = 1 p G i k + F k j + E i j
where Cij signifies the content of PTEs j in soil sample I, p represents the number of factors (i.e., pollution sources), Gik shows contribution of factor k to soil sample I, Fkj denotes the content of PTEs j in factor k, and Eij stands for the residual.
Additional information on the procedures, methods, and formulas used in this study for determining the soil or site contamination level through the PMF model was followed as specified by [55,62] and as applied by [3,6].

2.5. Statistical Analysis and Spatial Modeling

Basic statistical parameters (such as mean, median, minimum, maximum, standard deviation, and coefficient of variance) were first calculated for each soil property based on horizon. Positive matrix factorization (PMF, EPA version 5.0, Washington, DC, USA) was used for the estimation of source apportionment and contamination level of the PTEs. To determine the relationship between the PTEs in organic and mineral soils, an ANOVA and correlation analysis were used. Ordinary kriging interpolation was used in determining the differences and/or similarities among sites with a proportional distance among them. The interpolation technique enhanced the creation of the spatial distribution maps of the PTEs of the study area. ArcGIS, version 10.7.1, CA, USA [63], was used for processing and visualizing of the spatial data. By applying the ordinary kriging technique, maps of the spatial distribution of these soil properties were generated [64]. The result was validated using the mean error [65,66]. In other words, to determine the accuracy of the produced maps, the mean error (ME) was used for the validation. The formula is shown below in Equation (4):
ME = i = 1 n x 1 . i x 2 . j / n
where x1 is prediction of the variable x, x2 is measure of that variable, and n is number of records.

3. Results and Discussion

3.1. General Description of PTEs Concentrations and Their Spatial Distribution in the Soil

The basic statistical characteristics of the studied PTEs including Cr, Cu, Fe, Mn, and Pb for the organic soil and mineral soil have been described in Table 2. The coefficient of variation (CV) defines the degree of variations within PTE concentrations [67]. A coefficient of variation value less than 20% represents low variability, and a CV that falls between 21–50% is regarded as moderate variability. On the other hand, a CV ranging from 50–100% signifies high variability, while a CV greater than 1 (that is >100%) is described as having extremely high heterogeneity. In this study, the CV of the PTEs in the organic soil increased in the following order: Fe < Pb < Cu < Cr < Mn, accounting for 46.31%, 49.40%, 54.94%, 96.36%, and 97.06%, respectively.
The CV of the PTEs for mineral soil was also in ascending order: Fe (60.12%) < Pb (64.48%) < Cr (66.70%) < Mn (81.9%) < Cu (96.88%). The results derived from the CV revealed a high variability between the PTEs in the mineral soil. Similarly, in the organic soil, the CV for Cr, Mn, and Cu indicated a high variability (Table 2). In general, both the organic soil and the mineral soil revealed high heterogeneity (or variability). All of the PTEs showed relatively high variability in both soil horizons except for the Fe (46.31%) in the organic soil. The spatial distribution of the heterogeneity of the PTEs suggest that the metals are enriched by intensive sources of from the industrial, commercial, domestic and agricultural sectors [3,5,13]. However, the content of the PTEs varies between the soil horizons, yet the organic soil had higher mean values across the metals, excluding Mn. The content of manganese was 18.4 mg kg−1 higher in the mineral soil compared with its content in the organic soil. This finding agrees with a report by other authors in the same region [21]. Studies have shown that, in addition to human activities and their associated soil acidifications [13], the geological bedrock of the area also contributes to accruing PTE concentrations [15]. The study area has Podzols and Dystric Cambisols as the prevailing soils [5], and this might have contributed to the high contents of Pb, Fe, and Mn. The concentration of Fe in the study area is remarkable when compared with other metals. This could be attributed to the high acidic soil status of the area (Table 1). As has been earlier reported, the concentration of Fe in the soil solution at optimal soil pH falls between 30 and 550 μg L−1, but in high acid soils, it may exceed 2000 μg L−1 [68]. Higher concentrations of the PTEs were found in this study area relative to the neighboring regions in the country [3]. In comparison with the European value [10], the world value [10], and the Crati Basin value [69], the contents of Pb in both the organic soil and the mineral soil were higher. The exceptional content of Pb in the study area might be attributed to past intensive anthropogenic activities and the prevailing geological formation of Cambisols. Lead has been reported to exhibit the highest content in a Cambisols soil group [10].

3.2. Relationships among the PTE Concentrations in the Organic and Mineral Soils

The correlation analysis for the PTE concentrations in the organic and mineral soils revealed that the content of Cr in the mineral soil (Cr_tot_A) showed a significant and strong positive correlation with Cu, Fe, Mn, and Pb in the mineral soil (Table 3). Chromium has a strong relationship with other elements because it is easily mobilized in acidic soils and our study area is highly acidic [5,10]. Furthermore, in the mineral soil, Pb is significantly correlated with Cu and Mn. In the organic soil, Pb has a significant and strong positive relation with Cu and Fe. The strong relation between Pb, Fe, and Mn was documented earlier [10]. The concentrations of Pb in Fe–Mn nodules can be as high as 20,000 mg kg−1 [70]. Most of the negative correlations between the PTEs occurred in the inter-horizon and not within the same horizon. The correlation between the elements in the same soil horizon showed more positive relationships. This could be described by the likelihood that they shared the same origin. Furthermore, the correlation of the PTEs in the mineral soil revealed stronger relationships when compared with the correlation in the organic soil. This probably proved that these PTEs are more closely associated in the mineral soil relative to the organic soil. This finding was consistent with a recent report by other authors on the same issue [3].
The ANOVA in Table 4 was used to analyze the distribution of PTE contents in relation to the organic soil and the mineral soil horizons. It was revealed that all of the elements, with the exception Mn, showed significant relationships in both the organic soil and the mineral soil. There have been few studies within and outside the study area that focused on the relationship between soil horizon [8,21]. Consistent with our study, many authors have reported a significant relationship between soil horizon, elevation, and the concentrations of metals. For example, in the Suxian district of Chenzhou City in Hunan Province of China, it was revealed that heavy metal concentrations decreased at low elevation but increased considerably with increasing elevations [8]. Other studies have affirmed that fine-particle metals including Cr and Cu accumulate more at lower elevations [11].
In the organic soil, the highest contents of Cr, Pb, Fe, and Mn were found in the northern and central parts of the area (Figure 2). Studies have revealed a close association between Mn and Fe. Manganese is described as a member of the iron family, and both elements are closely linked in geochemical processes [10]. The author further stressed that Mn cycles follow Fe cycles in various terrestrial environments. Copper on the other hand had a concentration hotspot that spread from the northeast to the northwest.
The kriged map of Cu and Pb distributions showed almost the same pattern in the mineral soil (Figure 3). They tend to have higher concentrations towards the east in this area. In this mineral soil, Fe and Cr showed extensive spatial distribution patterns that spread from the northeast, through the central region, and to the northwestern part of the mapped area. On the other hand, the kriged map of Mn distribution showed higher concentrations within the northwest and north-central parts of the area. In general, for both organic and mineral soils, the northern and central parts of the kriged maps revealed more distribution of the elements when compared with the southern part. This could be explained by the historical distribution pattern of industrial and agricultural activities in the study area, which were mainly located in the northern and central parts [71].

3.3. Source Apportionment by the Positive Matrix Factorization (PMF) Model

There have been many reports and studies that affirmed that our study area is located within the vicinity of various anthropogenic activities including mining, intensive agriculture, automobile gas emissions, and acute biological sludge, which might affect the soil [4,52,53,72,73]. It is important to examine the sources of PTEs in the study area. Therefore, the PMF model was adopted as one of the best and latest models with high functionalities for effective PTE source identification [6]. The validity and reliability of the analysis are based on minimum Q to model the residual matrix that influences a substantial number of variables. To derive the best result, the PMF model was run 20 times, while the best outputs (which were Run 8 and Run 20) were selected following the software developer’s guide [55]. The PMF analysis produced six factors (see Figure 4 and Figure 5) and disclosed the origin of contributions based on each PTE [62].
In the organic soil, factor 1 was dominated by Pb and Cr, with a factor loading of 73.5% and 23.7%, while in the mineral soil, Cu and Pb had factor loads of 92.2% and 18.3%, respectively (Figure 4 and Figure 5). Factor 2 was characterized by Cr and Mn, and these accrued factor loadings of 45.1% and 15.1%, respectively, in the organic soil layer. On the other hand, in the mineral soil, Mn and Cu accumulated factor loadings of 28.3% and 6.7%, respectively. The factor 3 load consisted mainly of Mn (72.6%) and Cu (4.9%) in the organic soil, while Fe (64.8%) and Cr (11.9%) accrued in the mineral soil. Furthermore, factor 4 load was characterized by Cu (51.7%) and Fe (18.3%) in the organic soil while Mn (71.1) and Pb (26.4%) accumulated in the mineral soil. The factor 5 had Fe (42.5%), Cr (19.6%), and Pb (14.2%), accounting for the highest elements in the organic soil, while Cr (72.9%) and Fe (20.1%) accrued in the mineral soil. The 6th factor (which is the last factor) revealed that Cu (15.7%) had the highest factor load in the organic and mineral soils relative to all of the studied elements. Cr, Cu, and Pb accumulations in the study probably confirmed intensive pollution from many sources such as industrial, agricultural, commercial, and municipal activities and wastes [28]. These PTEs might have also been deposited during weathering because, in mineral forms, most of the elements are oxidized, released, and reprecipitated in the soil [74]. However, the sources of Mn and Fe in the soil are natural sources, posing no threats if the concentrations do not exceed the maximum allowable limits. However, the sources of Mn and Fe in the soil are natural sources, posing no threats if the concentrations do not exceed the maximum allowable limits.

3.4. Contamination Factor and Pollution Load Index of PTEs for Organic Soil and Mineral Soil

The research revealed that, in organic soil, there was no obvious pollution recorded in any of the soil sample points except in samples 33, 36, 65, 81, 103, 113, 147, 157, and 193. On the other hand, in the mineral horizon, the soil at some sites were polluted, observed only in samples 111, 126, 213, and 221 horizon (Table S1 and Table S2). This finding, however, is in contrast with the reports from previous studies in the same area. For example, there are studies showing that, in the past few decades, the area was among the major pollution zones caused by industrial and agricultural activities [5,75]. Furthermore, as one of the regions located in northern Bohemia, the study area has been documented as a region characterized by power and coal production from the 1950s to the 1980s [75]. In addition to intensive agriculture, industrial activities, and geological processes, the study area has some peatbogs, and these increased the pollution of the area by PTEs to a large extent [76,77]. This study area lies in the “Black Triangle”, which is affected by industrial activities linked to the extraction and exploitation of coal and other natural resources on the sides of the Jizera Mts areas [20,21,72,73]. The study area is historically susceptible to pollution, but the current results revealed otherwise. This could be due to measures put in place by authorities to ameliorate the impact of industrial activities in the area.

4. Conclusions

In both organic and the mineral soils, a high variability in the PTEs was observed. The spatial distribution and the heterogeneity of the PTEs suggested that there were widely distributed sources of the metals’ enrichment from the industrial, commercial, domestic, and agricultural sectors. Lead revealed a high concentration level. Chromium showed a strong relationship to other elements investigated, while Pb has a significant and strong positive relation with Cu. This was probably because Cr is easily mobilized in acidic soils and our study area is highly acidic.
The correlation between the elements in the same soil layer showed more positive relationships, while the correlation of the PTEs in the mineral soil revealed stronger relationships when compared with the correlation in the organic soil. The findings revealed that all the elements in exemption of Mn indicated significant relationships in both the organic soil and the mineral soil. Meanwhile, the concentrations of Mn and Fe were not harmful in the study area. On the other hand, Mn and Fe we below the EU and World limits. In both the organic soil and the mineral soil, the northern part of the kriged maps revealed more distribution of the elements when compared with either the southern part. This implied that the concentrations of the elements were oriented towards the northern part of the study region.
The applications of the Positive Matrix Factorization (PMF) model, ArcGIS-based ordinary kriging, and contamination level analysis were effective for source identification, hotspot location, and assessment of the contamination level of the PTEs. In the organic soil, there was no obvious pollution recorded in any of the soil sample points except in samples 33, 36, 65, 81, 103, 113, 147, 157, and 193. On the other hand, in the mineral horizon, some deteriorated site quality was observed only in samples 111, 126, 213, and 221. The method and results presented might be applicable in coniferous and broad-leaf tree-dominated highlands with a history of industrial activities and atmospheric acidifications. The methods are suitable for measuring the distribution and concentration of elements.
The current result revealed that there is no evidence of pollution by PTEs in the Jizera Mts area. In contrast, this area lies within the “Black Triangle”, which was affected by industrial activities linked with the extraction and exploitation of coal and other natural resources in Central Europe. However, the rate of pollution in the area is very low based on the findings of this study. There may be a need for intermittent assessment of the soil. This regular assessment will help to curtail the possibility of any excessive accumulation and escalation in the future. The findings from this study may serve as a baseline for pollution assessments of farmland and forest soil quality in the Czech Republic and in Europe. The results might support policy-developers in sustainable farming and forestry for the health of the ecosystem to achieve food security, forest safety, as well as animal and human welfare.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/toxics9080181/s1, Table S1: Contamination Factor (CF) and Pollution Loading Index (PLI) for PTEs in organic soil (n = 221), Table S2: Contamination Factor (CF) and Pollution Loading Index (PLI) for PTEs in mineral soil (n = 221).

Author Contributions

Conceptualization, S.K.A. and L.B.; data curation, C.N.; formal analysis, C.N.; funding acquisition, L.B. and S.K.A.; investigation, V.Y.O.S.; methodology, C.N.; resources, all authors; software, C.N.; visualization, all authors; writing—original draft, C.N. and S.K.A.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded (1) by the internal PhD Grant No. SV20-5-21130 of the Faculty of Agrobiology, Food and Natural Resources of the Czech University of Life Sciences Prague (CZU) and (2) with support from the Centre of Excellence (a center for the investigation of synthesis and transformation of nutritional substances in the food chain in interactions with potential risk implications of anthropogenic origin: a comprehensive assessment of the soil contamination risks for the quality of agricultural products, NutRisk Centre, supported by the European Union and Ministry of Education, Youth and Sports of the Czech Republic (project No. CZ.02.1.01/0.0/0.0/16_019/0000845)).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are contained within the article and Supplementary Materials.

Acknowledgments

We appreciate Karel Němeček for providing the study area shapefile.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alloway, B.J. Micronutrient Deficiencies in Global Food Production; Springer: Heidelberg, Germany, 2008. [Google Scholar]
  2. Eriksson, J.; Dahlin, S.A.; Sohlenius, G.; Söderström, M.; Öborn, I. Spatial patterns of essential trace element concentrations in Swedish soils and crops. Geoderma Reg. 2017, 10, 163–174. [Google Scholar] [CrossRef]
  3. Agyeman, P.C.; Ahado, S.K.; Kingsley, J.; Kebonye, N.M.; Biney, J.K.M.; Borůvka, L.; Vasat, R.; Kocarek, M. Source apportionment, contamination levels, and spatial prediction of potentially toxic elements in selected soils of the Czech Republic. Environ. Geochem. Health 2020, 43, 601–620. [Google Scholar] [CrossRef] [PubMed]
  4. Barker, B.J.; Clausen, J.L.; Douglas, T.L.; Bednar, A.J.; Griggs, C.S.; Martin, W.A. Environmental impact of metals resulting from military training activities: A review. Chemosphere 2021, 265, 129110. [Google Scholar] [CrossRef]
  5. Borůvka, L.; Vašát, R.; Němeček, K.; Novotný, R.; Šrámek, V.; Vacek, O.; Pavlů, L.; Fadrhonsová, V.; Drábek, O. Application of regression-kriging and sequential Gaussian simulation for the delineation of forest areas potentially suitable for liming in the Jizera Mountains region, Czech Republic. Geoderma Reg. 2020, 21, e00286. [Google Scholar] [CrossRef]
  6. Chen, L.; Lowenthal, D.H.; Watson, J.G.; Koracin, D.; Kumar, N.; Knipping, E.M.; Wheeler, N.; Craig, K.; Reid, S. Toward effective source apportionment using positive matrix factorization: Experiments with simulated PM2.5 data. J. Air Waste Manag. Assoc. 2010, 60, 43–54. [Google Scholar] [CrossRef]
  7. Chien, C.S.; Wang, H.; Chen, Y.; Wang, M.; Liu, C. Removal of heavy metals from contaminated paddy soils using chemical reductants coupled with dissolved organic carbon solutions. J. Hazard. Mater. 2021, 403, 123549. [Google Scholar] [CrossRef] [PubMed]
  8. Ding, Q.; Cheng, G.; Wang, Y.; Zhuang, D. Effects of natural factors on the spatial distribution of heavy metals in soils surrounding mining regions. Sci. Total Environ. 2017, 578, 577–585. [Google Scholar] [CrossRef]
  9. Gerdol, R.; Bragazza, L. Effects of altitude on element accumulation in alpine moss. Chemosphere 2020, 64, 810–816. [Google Scholar] [CrossRef] [PubMed]
  10. Kabata-Pendias, A. Trace Elements in Soils and Plants, 4th ed.; CRC Press Taylor & Francis Group: New York, NY, USA, 2011. [Google Scholar]
  11. Liu, X.; Shi, H.; Bai, Z.; Zhou, W.; Liu, K.; Wang, M.; He, Y. Heavy metal concentrations of soils near the large opencast coal mine pits in China. Chemosphere 2020, 244, 125360. [Google Scholar] [CrossRef]
  12. Marcantonio, R.A.; Field, S.P.; Sesay, P.B.; Lamberti, G.A. Identifying human health risks from precious metal mining in Sierra Leone. Reg. Environ. Chang. 2021, 21, 1–12. [Google Scholar] [CrossRef]
  13. Pavlů, L.; Drábek, O.; Stejskalová, Š.; Tejnecký, V.; Hradilová, M.; Nikodem, A.; Borůvka, L. Distribution of aluminium fractions in acid forest soils: Influence of vegetation changes. iForest Biogeosci. For. 2018, 11, 721–727. [Google Scholar] [CrossRef]
  14. Peli, M.; Bostick, B.C.; Barontini, S.; Lucchini, R.G.; Ranzi, R. Profiles and species of Mn, Fe and trace metals in soils near a ferromanganese plant in Bagnolo Mella (Brescia, IT). Sci. Total Environ. 2021, 755, 143123. [Google Scholar] [CrossRef]
  15. Salvador-Blanes, S.; Cornu, S.; Bourennane, H.; King, D. Controls of the spatial variability of Cr concentration in topsoils of a central French landscape. Geoderma 2006, 132, 143–157. [Google Scholar] [CrossRef]
  16. Sungur, A.; Kavdir, Y.; Özcan, H.; İlay, R.; Soylak, M. Geochemical fractions of trace metals in surface and core sections of aggregates in agricultural soils. Catena 2021, 197, 104995. [Google Scholar] [CrossRef]
  17. Kalkhajeh, Y.K.; Huang, B.; Hu, W.; Ma, C.; Gao, H.; Thompson, M.L.; Hansen, H.C.B. Environmental soil quality and vegetable safety under current greenhouse vegetable production management in China. Agric. Ecosyst. Environ. 2021, 307, 107230. [Google Scholar] [CrossRef]
  18. Tian, H.; Zhang, C.; Qi, S.; Kong, X.; Yue, X. Concentration and Spatial Distribution of Potentially Toxic Elements in Surface Soil of a Peak-Cluster Depression, Babao Town, Yunnan Province, China. Int. J. Environ. Res. Public Health 2021, 18, 3122. [Google Scholar] [CrossRef]
  19. Shar, S.; Reith, F.; Ball, A.S.; Shahsavar, E. Long-term Impact of Gold and Platinum on Microbial Diversity in Australian Soils. Microb. Ecol. 2021, 81, 977–989. [Google Scholar] [CrossRef]
  20. Matys Grygar, T.; Nováková, T.; Bábek, O.; Elznicová, J.; Vadinová, N. Robust assessment of moderate heavy metal contamination levels in floodplain sediments: A case study on the Jizera River, Czech Republic. Sci. Total Environ. 2013, 452–453, 233–245. [Google Scholar] [CrossRef]
  21. Kváčová, M.; Ash, C.; Borůvka, L.; Pavlů, L.; Nikodem, A.; Němeček, K.; Tejnecký, V.; Drábek, O. Contents of Potentially Toxic Elements in Forest Soils of the Jizera Mountains Region. Environ. Model. Assess. 2015, 20, 183–195. [Google Scholar] [CrossRef]
  22. Demková, L.; Jezný, T.; Bobuľská, L. Assessment of soil heavy metal pollution in a former mining area before and after the end of mining activities. Soil Water Res. 2017, 12, 229–236. [Google Scholar] [CrossRef] [Green Version]
  23. Kabala, C.; Galka, B.; Jezierski, P. Assessment and monitoring of soil and plant contamination with trace elements around Europe’s largest copper ore tailings impoundment. Sci. Total Environ. 2020, 738, 139918. [Google Scholar] [CrossRef]
  24. Juhos, K.; Czigány, S.; Madarász, B.; Ladányi, M. Interpretation of soil quality indicators for land suitability assessment using multivariate approach for Central European arable soils. Ecol. Indic. 2019, 99, 261–272. [Google Scholar] [CrossRef]
  25. Ahmadi, M.; Akhbarizadeh, R.; Haghighifard, N.J.; Barzegar, G.; Jorfi, S. Geochemical determination and pollution assessment of heavy metals in agricultural soils of south western of Iran. J. Environ. Health Sci. Eng. 2019, 17, 657–669. [Google Scholar] [CrossRef]
  26. Hoaghia, M.A.; Levei, E.A.; Cadar, O.; Senila, M.; Hognogi, G.G. Assessment of metal contamination and ecological risk in urban soils situated near a metallurgical complex. Environ. Eng. Manag. J. 2017, 16, 1623–1630. [Google Scholar] [CrossRef]
  27. Alloway, B. Heavy Metals in Soils; Chapman and Hall: London, UK, 1995. [Google Scholar]
  28. Nwaogu, C.; Ogbuagu, D.H.; Abrakasa, S.; Olawoyin, M.A.; Pavlů, V. Assessment of the impacts of municipal solid waste dumps on soils and plants. Chem. Ecol. 2017, 33, 589–606. [Google Scholar] [CrossRef]
  29. Mondal, S.; Singh, G. Pollution evaluation, human health effect and tracing source of trace elements on road dust of Dhanbad, a highly polluted industrial coal belt of India. Environ. Geochem. Health 2021, 43, 2081–2103. [Google Scholar] [CrossRef] [PubMed]
  30. Luo, X.S.; Yu, S.; Zhu, Y.G.; Li, X.D. Trace metal contamination in urban soils of China. Sci. Total Environ. 2012, 421–422, 17–30. [Google Scholar] [CrossRef] [PubMed]
  31. Bhuiyan, M.A.H.; Bodrud-Doza, M.; Rakib, M.A.; Saha, B.B.; Didar-Ul Islam, S.D. Appraisal of pollution scenario, sources and public health risk of harmful metals in mine water of Barapukuria coal mine industry in Bangladesh. Environ. Sci. Pollut. Res. 2021, 28, 22105–22122. [Google Scholar] [CrossRef]
  32. Simionov, I.A.; Cristea, D.S.; Petrea, S.M.; Mogodan, A.; Nicoara, M.; Plavan, G.; Baltag, E.S.; Jijie, R.; Strungaru, S.A. Preliminary investigation of lower Danube pollution caused by potentially toxic metals. Chemosphere 2021, 264, 128496. [Google Scholar] [CrossRef]
  33. Tamilmani, A.; Venkatesan, G. Assessment of trace metals and its pollution load indicators in water and sediments between Upper and Grand Anicuts in the Cauvery. Int. J. Environ. Sci. Technol. 2021. [Google Scholar] [CrossRef]
  34. Komárek, M.; Ettler, V.; Chrastný, V.; Mihaljevic, M. Lead isotopes in environmental sciences: A review. Environ. Int. 2008, 34, 562–577. [Google Scholar] [CrossRef]
  35. Brüggemann, J.; Dörfner, H.H.; Hecht, H.; Kumpulainen, J.T.; Westermair, T. Status of Trace Elements in Staple Foods from Germany 1990–1994. In Trace Elements, Natural Antioxidants and Contaminants in European Foods and Diets; Kumpulainen, J.T., Ed.; FAO: Rome, Italy, 1996; Volume 49, p. 5. [Google Scholar]
  36. Nicholson, F.A.; Smith, S.R.; Alloway, B.J.; Carlton-Smith, C.; Chambers, B.J. An inventory of heavy metals inputs to agricultural soils in England and Wales. Sci. Total Environ. 2003, 311, 205–219. [Google Scholar] [CrossRef]
  37. Staniszewski, P.; Bilek, M.; Szwerc, W.; Tomusiak, R.; Osiak, P.; Kocjan, R.; Moskalik, T. The effect of tree age, daily sap volume and date of sap collection on the content of minerals and heavy metals in silver birch (Betula pendula Roth) tree sap. PLoS ONE 2020, 15, e0244435. [Google Scholar] [CrossRef]
  38. Ren, H.Y.; Zhuang, D.F.; Singh, A.; Pan, J.J.; Qiu, D.S.; Shi, R.H. Estimation of as and Cu contamination in agricultural soils around a mining area by reflectance spectroscopy: A case study. Pedosphere 2009, 19, 719–726. [Google Scholar] [CrossRef]
  39. Gholizadeh, A.; Boruvka, L.; Saberioon, M.; Vasat, R. Visible, near-infrared, and mid-infrared spectroscopy applications for soil assessment with emphasis on soil organic matter content and quality: State-of-the-art and key issues. Appl. Spectrosc. 2013, 67, 1349–1362. [Google Scholar] [CrossRef]
  40. Gholizadeh, A.; Boruvka, L.; Vasat, R.; Saberioon, M.; Klement, A.; Kratina, J.; Tejnecky, V.; Drabek, O. Estimation of potentially toxic elements contamination in anthropogenic soils on a Brown coal mining dumpsite by reflectance spectroscopy: A case study. PLoS ONE 2015, 10, e0117457. [Google Scholar] [CrossRef] [Green Version]
  41. Gholizadeh, A.; Saberioon, M.; Ben-Dor, E.; Rossel, R.A.V.; Boruvka, L. Modelling potentially toxic elements in forest soils with viseNIR spectra and learning algorithms. Environ. Pollut. 2020, 267, 115574. [Google Scholar] [CrossRef]
  42. Sun, X.; Zhang, L.; Lv, J. Spatial assessment models to evaluate human health risk associated to soil potentially toxic elements. Environ. Pollut. 2021, 268, 115699. [Google Scholar] [CrossRef] [PubMed]
  43. Chodak, M.; Khanna, P.; Horvath, B.; Beese, F. Near infrared spectroscopy for determination of total and exchangeable cations in geologically heterogeneous forest soils. J. Near Infrared Spectrosc. 2004, 12, 315–324. [Google Scholar] [CrossRef]
  44. Hang, X.; Wang, H.; Zhou, J.; Ma, C.; Du, C.; Chen, X. Risk assessment of potentially toxic element pollution in soils and rice (Oryza sativa) in a typical area of the Yangtze River Delta. Environ. Pollut. 2009, 157, 2542–2549. [Google Scholar] [CrossRef]
  45. Liu, M.; Wang, T.; Skidmore, A.K.; Liu, X. Heavy metal-induced stress in rice crops detected using multi-temporal Sentinel-2 satellite images. Sci. Total Environ. 2018, 637–638, 18–29. [Google Scholar] [CrossRef] [PubMed]
  46. Arellano, P.; Tansey, K.; Balzter, H.; Boyd, D.S. Detecting the effects of hydrocarbon pollution in the Amazon forest using hyperspectral satellite images. Environ. Pollut. 2015, 205, 225–239. [Google Scholar] [CrossRef] [PubMed]
  47. Kong, F.; Chen, Y.; Huang, L.; Yang, Z.; Zhu, K. Human health risk visualization of potentially toxic elements in farmland soil: A combined method of source and probability. Ecotoxicol. Environ. Saf. 2021, 211, 111922. [Google Scholar] [CrossRef]
  48. Wang, Z.; Chen, X.; Yu, D.; Zhang, L.; Wang, J.; Lv, J. Source apportionment and spatial distribution of potentially toxic elements in soils: A new exploration on receptor and geostatistical models. Sci. Total Environ. 2021, 759, 143428. [Google Scholar] [CrossRef]
  49. Skeffington, R.A.; Cosby, B.J.; Whitehead, P.G. Long-term predictions of ecosystem acidification and recovery. Sci. Total Environ. 2016, 568, 381–390. [Google Scholar] [CrossRef] [Green Version]
  50. Akselsson, C.; Hultberg, H.; Karlsson, P.E.; Pihl Karlsson, G.; Hellsten, S. Acidification trends in south Swedish forest soils 1986–2008-Slow recovery and high sensitivity to sea-salt episodes. Sci. Total Environ. 2013, 444, 271–287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Marx, A.; Hintze, S.; Sanda, M.; Jankovec, J.; Oulehle, F.; Dusek, J.; Vitvar, T.; Vogel, T.; van Geldern, R.; Barth, J.A.C. Acid rain footprint three decades after peak deposition: Long-term recovery from pollutant sulphate in the Uhlirska catchment (Czech Republic). Sci. Total Environ. 2017, 598, 1037–1049. [Google Scholar] [CrossRef]
  52. EscartÍn, E.; Porte, C. Biomonitoring of PAH Pollution in High-Altitude Mountain Lakes through the Analysis of Fish Bile. Environ. Sci. Technol. 1999, 33, 406–409. [Google Scholar] [CrossRef]
  53. Kopáček, J.; Hejzlar, J.; Krám, P.; Oulehle, F.; Posch, M. Effect of industrial dust on precipitation chemistry in the Czech Republic (Central Europe) from 1850 to 2013. Water Res. 2016, 103, 30–37. [Google Scholar] [CrossRef] [PubMed]
  54. Křeček, J.; Palán, L.; Stuchlík, E. Impacts of land use policy on the recovery of mountain catchments from acidification. Land Use Policy 2019, 80, 439–448. [Google Scholar] [CrossRef]
  55. U.S. EPA. EPA Positive Matrix Factorization (PMF) 3.0 Model; U.S. Environmental Protection Agency: Research Triangle Park, NC, USA, 2010. Available online: http://www.epa.gov/heasd/products/pmf/pmf.html (accessed on 20 March 2021).
  56. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015. In International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports No. 106; FAO: Rome, Italy, 2015. Available online: http://www.fao.org/publications/card/en/c/942e424c-85a9-411d-a739-22d5f8b6cc41 (accessed on 2 December 2019).
  57. Vacek, O.; Vašát, R.; Borůvka, L. Quantifying the pedodiversity-elevation relations. Geoderma 2020, 373, 114441. [Google Scholar] [CrossRef]
  58. Tessier, A.; Campbell, P.G.C.; Bisson, M. Sequential extraction procedure for the speciation of particulate trace metals. Anal. Chem. 1979, 51, 844–851. [Google Scholar] [CrossRef]
  59. Pavlů, L.; Borůvka, L.; Drábek, O.; Nikodem, A. Effect of natural and anthropogenic acidification on aluminium distribution in forest soils of two regions in the Czech Republic. J. For. Res. 2021, 32, 363–370. [Google Scholar] [CrossRef] [Green Version]
  60. Melo, B.A.G.; Motta, F.L.; Santana, M.H.A. Humanic acids: Structural properties and multiple functionalities for novel technical development. Mater. Sci. Eng. 2016, 62, 967–974. [Google Scholar] [CrossRef] [PubMed]
  61. 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änder Meeresuntersuchungen 1980, 33, 566–575. [Google Scholar] [CrossRef] [Green Version]
  62. Norris, G.; Duvall, R.; Brown, S.; Bai, S. EPA Positive Matrix Factorization (PMF) 5.0 Fundamentals and User Guide; U.S. Environmental Protection Agency, Office of Research and Development: Washington, DC, USA, 2014.
  63. ESRI. ArcGIS Desktop: Release 10; Environmental Systems Research Institute: Redlands, CA, USA, 2019. [Google Scholar]
  64. Webster, R.; Oliver, M.A. Geostatistics for Environmental Scientists, 2nd ed.; John Wiley & Sons Ltd.: Chichester, UK, 2007. [Google Scholar]
  65. Brus, D.J.; Kempen, B.; Heuvelink, G.B.M. Sampling for validation of digital soil maps. Eur. J. Soil Sci. 2011, 62, 394–407. [Google Scholar] [CrossRef]
  66. Vašát, R.; Pavlů, L.; Borůvka, L.; Drábek, O.; Nikodem, A. Mapping the topsoil pH and humus quality of forest soils in the north bohemian Jizerské hory Mts. Region with ordinary, universal and regression kriging: Cross-validation comparison. Soil Water Res. 2013, 8, 97–104. [Google Scholar] [CrossRef] [Green Version]
  67. Karami, M.; Afyuni, M.; Khoshgoftarmanesh, A.H.; Papritz, A.; Schulin, R. Grain zinc, iron, and copper concentrations of wheat grown in central Iran and their relationships with soil and climate variables. J. Agric. Food Chem. 2009, 57, 10876–10882. [Google Scholar] [CrossRef]
  68. Kabata-Pendias, A.; Wiacek, K. Excessive uptake of heavy metals by plants from contaminated soil. Soil Sci. Ann. 1985, 36, 33. [Google Scholar]
  69. Guagliardi, I.; Cicchella, D.; De Rosa, R. A geostatistical approach to assess concentration and spatial distribution of heavy metals in urban soils. Water Air Soil Pollut. 2012, 223, 5983–5998. [Google Scholar] [CrossRef]
  70. Kabata-Pendias, A.; Pendias, H. Biogeochemistry of Trace Elements, 2nd ed.; Wyd. Nauk PWN: Warsaw, Poland, 1999. [Google Scholar]
  71. Černik, J.; Kunc, J.; Martinat, S. Territorial-technical and socio-economic aspects of successful brownfield regeneration: A case study of the liberec region (Czech Republic). Geogr. Tech. 2016, 11, 22–38. [Google Scholar] [CrossRef] [Green Version]
  72. Novák, M.; Erel, Y.; Zemanová, L.; Bottrel, S.H.; Adamová, M. A comparison of lead pollution record in Sphagnum peat with known historical Pb emission rates in the British Isles and the Czech Republic. Atmos. Environ. 2008, 42, 8997–9006. [Google Scholar] [CrossRef]
  73. Oulehle, F.; Hruska, J. Tree species (Picea abies and Fagus sylvatica) effects on soil water acidification and aluminium chemistry at sites subjected to long term acidification in the Ore Mts., Czech Republic. J. Inorg. Biochem. 2005, 99, 1822–1829. [Google Scholar] [CrossRef]
  74. Wei, W.; Ma, R.; Sun, Z.; Zhou, A.; Bu, J.; Long, X.; Liu, Y. Effects of Mining Activities on the Release of Heavy Metals (HMs) in a Typical Mountain Headwater Region, the Qinghai-Tibet Plateau in China. Int. J. Environ. Res. Public Health 2018, 15, 1987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Maján, G.; Kozak, M.; Püspöki, Z.; McIntosh, R.; Miko, L. Environmental geological examination of chromiumcontamination in Eastern Hungary. Environ. Geochem. Health 2001, 23, 229–233. [Google Scholar] [CrossRef]
  76. Glasshiem, E. Most, the Town that Moved: Coal, Communists and the ‘Gypsy Question’ in Post-War Czechoslovakia. Environ. Hist. 2007, 13, 447–476. [Google Scholar] [CrossRef]
  77. Zuna, M.; Mihaljevic, M.; Sebek, O.; Ettler, V.; Handley, M.; Navrátil, T.; Goliás, V. Recent lead deposition trends in the Czech Republic as recorded by peat bogs and tree rings. Atmos. Environ. 2011, 45, 4950–4958. [Google Scholar] [CrossRef]
Figure 1. Sampling points and land use–land cover of the Jizera Mts. area derived from the CORINE database (central Map), the location of the Jizera Mts. in Liberec region in the northern part of the Czech Republic (top left and down left maps), and the location of the Czech Republic in Europe (top right map).
Figure 1. Sampling points and land use–land cover of the Jizera Mts. area derived from the CORINE database (central Map), the location of the Jizera Mts. in Liberec region in the northern part of the Czech Republic (top left and down left maps), and the location of the Czech Republic in Europe (top right map).
Toxics 09 00181 g001
Figure 2. Spatial and vertical distributions of soil characteristics in the organic soil (org) assessed using ordinary kriging (all of the elements are reported in mg kg−1). Though Fe and Mn posed no threat in the region, they were mapped/modeled to draw inferences on Cr, Cu, and Pb. In other words, the study attempted to assess if the presence of Fe and Mn in the soil influenced the vertical and spatial distributions of the three other PTEs (namely Cr, Cu, and Pb) in the different soil horizons.
Figure 2. Spatial and vertical distributions of soil characteristics in the organic soil (org) assessed using ordinary kriging (all of the elements are reported in mg kg−1). Though Fe and Mn posed no threat in the region, they were mapped/modeled to draw inferences on Cr, Cu, and Pb. In other words, the study attempted to assess if the presence of Fe and Mn in the soil influenced the vertical and spatial distributions of the three other PTEs (namely Cr, Cu, and Pb) in the different soil horizons.
Toxics 09 00181 g002aToxics 09 00181 g002b
Figure 3. Spatial and vertical distributions of soil characteristics in the surface mineral soil (A) made by ordinary kriging. (all of the elements are in mg kg−1). Though Fe and Mn posed no threat in the region, they were mapped/modeled to draw inferences on Cr, Cu, and Pb. In other words, the study attempted to assess if the presence of Fe and Mn in the soil influenced the vertical and spatial distribution of the other three PTEs (namely Cr, Cu, and Pb) in the different soil horizons.
Figure 3. Spatial and vertical distributions of soil characteristics in the surface mineral soil (A) made by ordinary kriging. (all of the elements are in mg kg−1). Though Fe and Mn posed no threat in the region, they were mapped/modeled to draw inferences on Cr, Cu, and Pb. In other words, the study attempted to assess if the presence of Fe and Mn in the soil influenced the vertical and spatial distribution of the other three PTEs (namely Cr, Cu, and Pb) in the different soil horizons.
Toxics 09 00181 g003aToxics 09 00181 g003b
Figure 4. Source fingerprint of the total organic soil (tot-org) from the PMF model analysis showing the percentage contribution of PTEs. Note that the source and availability of Fe and Mn posed no risk to the soils.
Figure 4. Source fingerprint of the total organic soil (tot-org) from the PMF model analysis showing the percentage contribution of PTEs. Note that the source and availability of Fe and Mn posed no risk to the soils.
Toxics 09 00181 g004
Figure 5. Source fingerprint of the total mineral soil (tot-A) from the PMF model analysis showing the percentage contribution of PTEs. Note that the source and availability of Fe and Mn posed no risk to the soils.
Figure 5. Source fingerprint of the total mineral soil (tot-A) from the PMF model analysis showing the percentage contribution of PTEs. Note that the source and availability of Fe and Mn posed no risk to the soils.
Toxics 09 00181 g005
Table 1. Mean values of the physiochemical characteristics of the soil in the study area.
Table 1. Mean values of the physiochemical characteristics of the soil in the study area.
Properties (Unit)Organic SoilMineral Soil
Sand (%)29.728.2
Silt (%)44.225.3
Clay (%)26.146.5
TextureSandy clay-loamClay-loam
N (%)1.60.5
C (%)30.97.5
S (%)0.340.26
P (mg kg−1)946.9386.2
K (mg kg−1)811.6935.3
Ca (mg kg−1)915.2327.9
Mg (mg kg−1)839.51078.1
Al (mg kg−1)9473.58614.4
pH3.63.8
Table 2. Basic statistical characteristics of the soil PTE concentrations in the study area.
Table 2. Basic statistical characteristics of the soil PTE concentrations in the study area.
Soil HorizonsParameters CrCuFe Mn Pb
Organic soilCount221221221221221
Mean11.016.27357.8149.699.2
Median9.115.570107392.9
Mode7.218.510,20032104
Minimum3.12.310041.07.1
Maximum85.281.921,0001650339
Std dev10.68.93407.5145.249
Coef of Var. (CV)96.3654.9446.3197.0649.4
Mineral soilCount221221221221221
Mean4.56.46744.3168.065.6
Median3.83.86194.468.458.8
Mode3.91.03610248111
Minimum0.40.2159.30.56.7
Maximum26.538.324,274.01940.0281.0
Std dev3.06.24054.5137.642.3
Coef of Var. (CV) *
Czech Republic
66.7
<11.0
96.88
<16.0
60.12
>8000
81.9
<150.0
64.48
<60.0
** European mean value94.817.338,00052432
** World mean value59.538.935,00048827
** Crati Basin90.5444.3654,700130063.67
* Authors’ estimates from most publications in the Czech Republic on the issue as there was no official existing baseline; ** Kabata-Pendias [10]. All parameters and numbers are reported in mg kg−1, while CV is reported in %. Fe and Mn also showed reasonable variability; they posed no threat because their concentrations were far below the EU and world limits.
Table 3. Summary of correlation analyses between the PTE concentrations in the organic and mineral soils.
Table 3. Summary of correlation analyses between the PTE concentrations in the organic and mineral soils.
ParametersCr_tot_orgCu_tot_orgFe_tot_orgMn_tot_orgPb_tot_orgCr_tot_ACu_tot_AFe_tot_AMn_tot_APb_tot_A
Cr_tot_org1.00
Cu_tot_org0.031.00
Fe_tot_org0.77 *0.43 *1.00
Mn_tot_org0.52 **−0.400.191.00
Pb_tot_org0.000.71 **0.64 *−0.56 *1.00
Cr_tot_A0.53 *−0.140.100.58−0.211.00
Cu_tot_A0.10−0.13−0.100.20−0.26 *0.78 *1.00
Fe_tot_A0.030.160.50 *−0.060.430.56 **0.001.00
Mn_tot_A0.10−0.22 *−0.050.73 *−0.59 **0.60 *0.58 *0.031.00
Pb_tot_A0.05−0.12−0.080.10−0.100.76 **0.81 **0.190.54 *1.00
* = Correlation is significant at the 0.01 p-value, at <0.05; ** = correlation is significant at the 0.05 p-value; tot_org = total concentration in organic soil; tot_A = total concentration in the mineral soil.
Table 4. Summary of ANOVA for PTE concentrations for the tot-org and tot-A horizons.
Table 4. Summary of ANOVA for PTE concentrations for the tot-org and tot-A horizons.
Soil ParametersdfF-Statisticsp-Value *
Cr_tot_org2202.120.019
Cu_tot_org220−3.73<0.001
Fe_tot_org220−1.40<0.001
Mn_tot_org2203.310.685
Pb_tot_org2200.63<0.001
Cr_tot_A2202.060.016
Cu_tot_A220−1.940.021
Fe_tot_A220−1.13<0.001
Mn_tot_A2204.910.283
Pb_tot_A2200.820.041
* Figures (or values) in bold are significant at the 0.05 confidence level. tot_org = total concentration in organic soil; tot_A = total concentration in the mineral soil.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ahado, S.K.; Nwaogu, C.; Sarkodie, V.Y.O.; Borůvka, L. Modeling and Assessing the Spatial and Vertical Distributions of Potentially Toxic Elements in Soil and How the Concentrations Differ. Toxics 2021, 9, 181. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9080181

AMA Style

Ahado SK, Nwaogu C, Sarkodie VYO, Borůvka L. Modeling and Assessing the Spatial and Vertical Distributions of Potentially Toxic Elements in Soil and How the Concentrations Differ. Toxics. 2021; 9(8):181. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9080181

Chicago/Turabian Style

Ahado, Samuel Kudjo, Chukwudi Nwaogu, Vincent Yaw Oppong Sarkodie, and Luboš Borůvka. 2021. "Modeling and Assessing the Spatial and Vertical Distributions of Potentially Toxic Elements in Soil and How the Concentrations Differ" Toxics 9, no. 8: 181. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9080181

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