3.1. Distribution of Ag, As, Bi, Cd, Cu, W and Zn in Soil Samples
Previous studies carried out in the Panasqueira area [
6,
43,
44], allowed, through the topsoil concentrations, the characterization of the anthropogenic soil contamination due to the presence of tailings and open air impoundments.
Descriptive statistical parameters (median and range) for pH, SOM, Ag, As, Bi, Cd, Cu, W and Zn are summarized in
Table 4. The individual element values are compared with the corresponding local geochemical background levels, and also with reference values from the literature.
All samples, both topsoil and subsoil, have a pH lower than 7. In particular, topsoil pH values range from 3.2 to 6.6; 35.9% of topsoil samples have pH values between 3.0 and 4.0; 59.2% are between 4.0 and 5.0; 4.2% from 5.0 to 6.0, and 1.7% from 6.0 to 7.0,
i.
e., generally very acid to acid according to the United States Department of Agriculture classification [
81]. As shown in
Figure 2 (pH (a)), it is possible to identify that the areas with higher values are located around the villages Barroca Grande, S. Francisco de Assis, Rio, Barroca, Dornelas do Zêzere and S. Martinho. Around Rio and Barroca Grande tailings, the soil pH is above the 75th percentile (<4.1). The subsoil pH ranges from 3.6 to 6.2: 23.0% of the subsoil samples have pH values between 3.0 and 4.0; 73.5% are between 4.0 and 5.0; 2.7% from 5.0 to 6.0, and 0.9% from 6.0 to 7.0,
i.
e., also generally very acid to acid. The highest pH value found in subsoil samples was located in Rio village (
Figure 2 pH (b)) on the north section of the Rio tailings. In the surroundings of S. Francisco de Assis, Barroca Grande and Dornelas do Zêzere villages, it is possible to identify some agricultural soil with pH > 4.6. Possibly these slightly higher pH values are due to land agricultural beneficiation techniques. Nevertheless, the soil of the present study area is overall classified as very acid to acid.
Figure 2.
Spatial distribution of pH and soil organic matter (SOM, %) and Ag (mg kg−1) for topsoil (a) and subsoil (b). Percentile class limits at P5, P10, P25, P50, P75, P90, P95 and P97.5. In blue are represented the main water courses, in grey the villages and in white the main tailings and dams. The green dashed lines represent the Ag background levels.
Figure 2.
Spatial distribution of pH and soil organic matter (SOM, %) and Ag (mg kg−1) for topsoil (a) and subsoil (b). Percentile class limits at P5, P10, P25, P50, P75, P90, P95 and P97.5. In blue are represented the main water courses, in grey the villages and in white the main tailings and dams. The green dashed lines represent the Ag background levels.
Table 4.
Median, minimum and maximum concentrations values of topsoil and subsoil samples from Panasqueira mine and surrounding environment compared with international data. Casegas area, considered as representative of geochemical background (Bk).
Table 4.
Median, minimum and maximum concentrations values of topsoil and subsoil samples from Panasqueira mine and surrounding environment compared with international data. Casegas area, considered as representative of geochemical background (Bk).
Units | pH | SOM | Ag | As | Bi | Cd | Cu | W | Zn |
---|
| % | mg·kg−1 | mg·kg−1 | mg·kg−1 | mg·kg−1 | mg·kg−1 | mg·kg−1 | mg·kg−1 |
---|
Topsoil (n = 122) | Median | 4.1 | 14.6 | 0.1 | 65 | 0.9 | 0.2 | 35 | 6 | 70 |
Min–Max | 3.2–6.6 | 5.7–39.0 | 0.05–1.7 | 17–1503 | 0.3–65 | 0.1–3 | 10–292 | 0.2–200 | 22–199 |
Sk | 1.59 | 1.2 | 4.72 | 6.98 | 9.21 | 3.77 | 4.80 | 5.40 | 1.38 |
Subsoil (n = 116) | Median | 4.2 | 10.2 | 0.1 | 52 | 0.6 | 0.1 | 33 | 3 | 70 |
Min–Max | 3.6–6.2 | 3.6–19.5 | 0.1–0.6 | 8–350 | 0.2–15 | 0.1–1.5 | 12–146 | 0.1–29 | 16–192 |
Sk | 2.09 | 0.23 | 1.98 | 2.05 | 5.51 | 3.65 | 2.86 | 2.51 | 1.08 |
Geochemical background | Topsoil | 4.2 | 14.5 | 0.1 | 63 | 0.8 | 0.20 | 34 | 4.4 | 68 |
Subsoil | 4.3 | 10.0 | 0.1 | 49 | 0.6 | 0.10 | 31 | 2.5 | 68 |
Bk | 3.98 | 10.5 | 0.05 | 22 | 0.3 | 0.01 | 28 | 0.4 | 58 |
Data from literature | 6.0–8.0 a | 1–20 b | 0.07 c | 11 d | 0.3 c | 0.1 d | 16 d | 1 d | 55 d |
The topsoil SOM values range from 5.7% to 39.0% (see
Table 4 and
Figure 2 SOM (a)), while the subsoil SOM values vary from 3.6% to 19.5% (see
Table 4 and
Figure 2 SOM (b)). Comparing the topsoil and subsoil pH and SOM maps (
Figure 2 and
Table 4) is possible to observe that areas with the highest SOM values are related to areas with the lowest pH values.
To evaluate the presence of possible local anomalies the Ag (
Figure 2), As, Bi, Cd (
Figure 3), Cu, W, Zn (
Figure 4) median and maximum values, observed in topsoil, were compared with the corresponding results in subsoil samples. Although the results tabulated in
Table 4 show that Ag, As, Bi, Cu, W, and Zn present highest median and maximum values for topsoil when compared with subsoil, the paired samples
t-test indicates that there are significant differences between soil samples at both depths for the elements As, Bi and W. Additionally, the results clearly show that there are higher contents in soil samples at both depths, when compared with the local geochemical background. The behavior of the topsoil element concentrations may reflect the influence of the Barroca Grande and Rio tailings and open air impoundments, possibly caused by aerial transport and deposition.
Topsoil samples, as shown in
Table 5, have many strong to very strong correlation coefficients,
i.
e., As/Ag, As/Bi, As, Cd, As/Cu, As/W, Ag/Bi, Ag/Cu, As/W, Bi/Cd, Bi/Cu, Bi/W, Cd/Cu, Cd/W, Cd/Zn, Cu/W. The strongest correlations are between Ag/Cu (0.94), As/Ag (0.90) and Bi/W (0.90). There is also a positive intermediate correlation (
r = 0.31 to 0.70) between PTEs and pH. Several authors claim that heavy metal mobility holds a positive correlation with pH [
18,
88,
89]. It is well known that pH affects heavy metal availability, since it is the major factor in controlling the ability of soil to retain heavy metals in an exchangeable form [
21]. With low SOM the pH values may become relatively more important for the partitioning of metals. Most elements exhibit a weak negative correlation with SOM. For subsoil samples the PTEs correlation also presents some strong correlation coefficients (
Table 5), but they are overall lower than in topsoil. The stronger correlations are between As/Bi (
r = 0.75) and Bi/Cu (
r = 0.69). At this depth, pH is also the most significant soil property, presenting weak to intermediate positive correlations with the PTEs, except W.
Figure 3.
Spatial distribution of As, Bi, Cd (mg·kg−1) values for topsoil (a) and subsoil (b). Percentile class limits at P5, P10, P25, P50, P75, P90, P95 and P97.5. The green dashed lines represent the elements background levels. In blue are represented the main water courses, in grey the villages and in white the main tailings and dams.
Figure 3.
Spatial distribution of As, Bi, Cd (mg·kg−1) values for topsoil (a) and subsoil (b). Percentile class limits at P5, P10, P25, P50, P75, P90, P95 and P97.5. The green dashed lines represent the elements background levels. In blue are represented the main water courses, in grey the villages and in white the main tailings and dams.
Figure 4.
Spatial distribution of Cu, W, Zn (mg·kg−1) values for topsoil (a) and subsoil (b). Percentile class limits at P5, P10, P25, P50, P75, P90, P95 and P97.5. The green dashed lines represent the elements background levels. In blue are represented the main water courses, in grey the villages and in white the main tailings and dams.
Figure 4.
Spatial distribution of Cu, W, Zn (mg·kg−1) values for topsoil (a) and subsoil (b). Percentile class limits at P5, P10, P25, P50, P75, P90, P95 and P97.5. The green dashed lines represent the elements background levels. In blue are represented the main water courses, in grey the villages and in white the main tailings and dams.
Table 5.
Pearson’s linear product-moment correlation coefficients for Ag, As, Bi, Cd, Cu, W, Zn, pH and SOM in topsoil (n = 122) and subsoil (n = 116) samples.
Table 5.
Pearson’s linear product-moment correlation coefficients for Ag, As, Bi, Cd, Cu, W, Zn, pH and SOM in topsoil (n = 122) and subsoil (n = 116) samples.
Var | | Ag | As | Bi | Cd | Cu | W | Zn | pH | SOM |
---|
Topsoil | Ag | 1.00 | 0.90 ** | 0.84 ** | 0.79 ** | 0.94 ** | 0.83 ** | 0.51 ** | 0.70 ** | −0.1 |
As | | 1.00 | 0.84 ** | 0.77 ** | 0.89 ** | 0.81 ** | 0.44 ** | 0.37 ** | −0.09 |
Bi | | | 1.00 | 0.82 ** | 0.87 ** | 0.90 ** | 0.44 ** | 0.54 ** | −0.14 |
Cd | | | | 1.00 | 0.86 ** | 0.79 ** | 0.68 ** | 0.42 ** | −0.07 |
Cu | | | | | 1.00 | 0.80 ** | 0.60 ** | 0.61 ** | −0.14 |
W | | | | | | 1.00 | 0.38 ** | 0.46 ** | 0.01 |
Zn | | | | | | | 1.00 | 0.31 ** | −0.15 |
pH | | | | | | | | 1.00 | −0.40 ** |
SOM | | | | | | | | | 1.00 |
Subsoil | Ag | 1.00 | 0.43 ** | 0.34 ** | 0.45 ** | 0.42 ** | 0.36 ** | 0.43 ** | 0.30 ** | 0.15 |
As | | 1.00 | 0.75 ** | 0.43 ** | 0.59 ** | 0.49 ** | 0.30 ** | 0.09 | −0.01 |
Bi | | | 1.00 | 0.50 ** | 0.69 ** | 0.56 ** | 0.35 ** | 0.04 | −0.01 |
Cd | | | | 1.00 | 0.53 ** | 0.43 ** | 0.57 ** | 0.17 | −0.09 |
Cu | | | | | 1.00 | 0.44 ** | 0.56 ** | 0.11 | −0.17 |
W | | | | | | 1.00 | 0.14 | −0.21 * | 0.16 |
Zn | | | | | | | 1.00 | 0.45 ** | −0.23 * |
pH | | | | | | | | 1.00 | −0.35 * |
SOM | | | | | | | | | 1.00 |
In order to identify the most important factor controlling the different PTEs spatial differences concentrations, a two-way ANOVA was performed [
90]. In this study, two models were selected to group the variables, (a) depth-SOM and (b) depth-pH. The results of the test between dependent (PTEs) and independent (depth and SOM) variables (
Table 6) showed that soil depth accounts for significant variations between the group means: Bi (
p = 0.035), Cd (
p = 0.015) and W (
p = 0.009). The independent variable SOM shows significant variations in the concentration of As (
p = 0.021), Bi (
p = 0.068), W (
p = 0.013) and Zn (
p = 0.021). There is no significant interaction between depth and SOM. While the independent variable pH (
Table 7) shows that there are significant variations with depth for most PTEs (Ag
p = 0.071; As
p = 0.027; Bi
p = 0.001; Cd
p <0.000; W
p <0.000), except for Cu (
p =0.371). pH presents a significant variation for Ag (
p = 0.031), W (
p = 0.001) and Zn (
p <0.001). The interaction between pH*depth shows no significant variations (
p >0.050), except on the concentration of Cu (
p = 0.038).
3.2. Quantitative Assessment of Soil Contamination
In this study, a simplified approach to assess soil contamination based on comparing the measured concentrations in the Panasqueira soil with the geochemical background values for Casegas was adopted.
Table 8 shows the results of contamination factors (
CF) and the modified degree of contamination (
mCd) for the selected elements in topsoil and subsoil, and also in Casegas soil.
Table 6.
Two-way ANOVA results between the dependent (PTEs) and independent (depth, SOM) variables. PTEs were subjected to log-normal transformation (As, As, Cd, Cu, W and Zn with α = 0.050; Bi with α = 0.075).
Table 6.
Two-way ANOVA results between the dependent (PTEs) and independent (depth, SOM) variables. PTEs were subjected to log-normal transformation (As, As, Cd, Cu, W and Zn with α = 0.050; Bi with α = 0.075).
Source of Variation | Dependent Variables | df | Mean Square | F | p-Value |
---|
Depth | Ag | 1 | 0.223 | 1.624 | 0.204 |
As | 1 | 0.127 | 1.652 | 0.200 |
Bi | 1 | 0.333 | 4.484 | 0.035 |
Cd | 1 | 0.926 | 5.951 | 0.015 |
Cu | 1 | 0.012 | 0.483 | 0.488 |
W | 1 | 1.486 | 6.928 | 0.009 |
Zn | 1 | 0.062 | 1.933 | 0.166 |
SOM | Ag | 1 | 0.010 | 0.076 | 0.783 |
As | 1 | 0.417 | 5.422 | 0.021 |
Bi | 1 | 0.251 | 3.371 | 0.068 |
Cd | 1 | 0.160 | 1.029 | 0.311 |
Cu | 1 | 0.007 | 0.275 | 0.601 |
W | 1 | 1.763 | 8.219 | 0.005 |
Zn | 1 | 0.173 | 5.372 | 0.021 |
Depth*SOM | Ag | 1 | 0.051 | 0.373 | 0.542 |
As | 1 | 0.007 | 0.093 | 0.760 |
Bi | 1 | 0.001 | 0.008 | 0.928 |
Cd | 1 | 0.001 | 0.006 | 0.940 |
Cu | 1 | 0.037 | 1.52 | 0.219 |
W | 1 | 0.095 | 0.44 | 0.506 |
Zn | 1 | 0.002 | 0.06 | 0.811 |
Error | Ag | 238 | 0.137 | - | - |
As | 215 | 0.077 | - | - |
Bi | 225 | 0.074 | - | - |
Cd | 233 | 0.156 | - | - |
Cu | 218 | 0.024 | - | - |
W | 215 | 0.215 | - | - |
Zn | 235 | 0.032 | - | - |
According the topsoil median values of the contamination factor, Bi, Cd and W present an extremely high degree of contamination, while Ag and As a high contamination factor. The results demonstrated that mCd values vary from the minimum 1.2 in both topsoil and subsoil to the maximum of 150 in topsoil, with median values varying from 4 in subsoil to 6 in topsoil. The cumulative frequency distribution indicates that only 7.4% of the soil samples were classified as no to low degree of pollution, with mCd values < 2.0, and the remaining soil samples (92.6%) registered moderate to ultra-high degree of pollution, with mCd values equal or greater than 2.0 (27.9% between 2 ≤ mCd < 4; 27.1% between 4 ≤ mCd < 8; 19.7% between 8 ≤ mCd < 16; 15.6% between 16 ≤ mCd < 32 and 2.5% between 32 ≤ mCd). The enrichment is more pronounced in topsoil (subsoil mCd values ranged 1.2–26.4).
Table 7.
Two-way ANOVA results between the dependent (PTEs) and independent (depth, pH) variables. PTEs were subjected to log-normal transformation (As, Bi, Cd, Cu, W, Zn with α = 0.050; Ag with α = 0.075).
Table 7.
Two-way ANOVA results between the dependent (PTEs) and independent (depth, pH) variables. PTEs were subjected to log-normal transformation (As, Bi, Cd, Cu, W, Zn with α = 0.050; Ag with α = 0.075).
Source of Variation | Dependent Variables | df | Mean Square | F | p-Value |
---|
Depth | Ag | 1 | 0.353 | 3.296 | 0.071 |
As | 1 | 0.512 | 4.930 | 0.027 |
Bi | 1 | 0.832 | 10.605 | 0.001 |
Cd | 1 | 1.810 | 13.449 | 0.000 |
Cu | 1 | 0.019 | 0.805 | 0.371 |
W | 1 | 3.055 | 14.560 | 0.000 |
Zn | 1 | 0.007 | 0.266 | 0.607 |
pH | Ag | 1 | 0.506 | 4.726 | 0.031 |
As | 1 | 0.031 | 0.301 | 0.584 |
Bi | 1 | 0.017 | 0.222 | 0.638 |
Cd | 1 | 0.020 | 0.146 | 0.703 |
Cu | 1 | 0.000 | 0.001 | 0.971 |
W | 1 | 2.315 | 11.033 | 0.001 |
Zn | 2 | 0.605 | 21.777 | 0.000 |
Depth*pH | Ag | 1 | 0.116 | 1.085 | 0.299 |
As | 1 | 0.053 | 0.510 | 0.476 |
Bi | 1 | 0.082 | 1.040 | 0.309 |
Cd | 1 | 0.179 | 1.332 | 0.250 |
Cu | 1 | 0.104 | 4.351 | 0.038 |
W | 1 | 0.321 | 1.532 | 0.217 |
Zn | 2 | 0.077 | 2.769 | 0.065 |
Error | Ag | 233 | 0.107 | - | - |
As | 230 | 0.104 | - | - |
Bi | 225 | 0.078 | - | - |
Cd | 226 | 0.135 | - | - |
Cu | 217 | 0.024 | - | - |
W | 213 | 0.210 | - | - |
Zn | 233 | 0.028 | - | - |
Table 8.
Contamination Factors (CF) and Modified Degree of Contamination (mCd) using geochemical background values. Casegas area is considered as representative of geochemical background.
Table 8.
Contamination Factors (CF) and Modified Degree of Contamination (mCd) using geochemical background values. Casegas area is considered as representative of geochemical background.
Var | AgCF | AsCF | BiCF | CdCF | CuCF | WCF | ZnCF | mCd |
---|
Bk | Min | 1.0 | 0.4 | 0.7 | 1.0 | 0.6 | 0.1 | 0.3 | 0.7 |
Med | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.1 | 1.0 | 1.4 |
Max | 8 | 5 | 9 | 30 | 6 | 26 | 2 | 13 |
Topsoil | Min | 1.0 | 0.8 | 1.0 | 5.0 | 0.4 | 0.6 | 0.4 | 1.2 |
Med | 2 | 3 | 3 | 20 | 1.2 | 17 | 1.2 | 6 |
Max | 34 | 68 | 215 | 300 | 10 | 571 | 3 | 150 |
Subsoil | Min | 1.0 | 0.4 | 0.7 | 5.0 | 0.4 | 0.1 | 0.3 | 1.2 |
Med | 2 | 2 | 2 | 10 | 1.1 | 8 | 1.2 | 4 |
Max | 12 | 16 | 51 | 150 | 5 | 84 | 3 | 26 |
Casegas is considered as an uncontaminated site, since it is outside the mine area and the influence of airborne polluted dust. The maximum mCd value of 13.1, represents a very high contamination degree, which is mostly due to the Cd and W contamination factors. In this case, however, the high mCd value is most likely due to geogenic sources, as this is a naturally enriched and mineralized zone.
The calculated
mCd values make possible the assessment of the spatial distribution of the modified degree of contamination. The first step is the determination of the spatial structure of the new variables, and the experimental variograms were used to model the data using exponential models, the extracted parameters are for topsoil: main direction = 90°; Nugget effect (
C0) = 30; Sill-Nugget effect (
C1) = 330; range of influence (length) = 1200 m; anisotropy ratio = 1.81; and for subsoil: main direction = 90°;
C0 = 0;
C1 = 19; length = 900 m; anisotropy ratio = 2.90. Estimation of the spatial distribution was then achieved by Ordinary Kriging and the respective map plotted. In
Figure 5 it is possible to observe the
mCd spatial distribution revealing areas with very high values.
Figure 5 shows that soil samples collected near the Barroca Grande tailings (A), Rio tailings (C) and the mud impoundments stand out clearly, because the soil is enriched in Ag, As, Bi, Cd, Cu, W and Zn; all soil samples from both depths taken from Barroca Grande exceed the As, Bi, Cd and W baseline values for Portugal, while at Rio only Zn in some samples (25%) has concentrations that are lower than the guide value. According to [
42], the Barroca Grande tailings and open impoundments have high As, Cd, Cu, Pb, W, and Zn concentrations (mean content in the more coarse tailings material As = 7142 mg·kg
−1 ; Cd = 56 mg·kg
−1 ; Cu = 2501 mg·kg
−1; Pb = 172 mg·kg
−1; Sn = 679 mg·kg
−1; W = 5400 mg·kg
−1 and Zn = 1689 mg·kg
−1 and mean content in the impoundment material (rejected from the mill operations) As = 44,252 mg·kg
−1 ; Cd = 491 mg·kg
−1; Cu = 4029 mg·kg
−1; Pb = 166 mg·kg
−1; Sn = 454 mg·kg
−1; W = 3380 mg·kg
−1 and Zn = 3738 mg·kg
−1). The mineralogy of these tailings consists of mainly quartz, muscovite, kaolinite, illite-montmorillonite, montmorillonite-vermiculite, and chlorite, and also arsenopyrite, wolframite, and natrojarosite.
Figure 5.
Spatial distribution of the modified degree of Contamination (mCd) for topsoil (a) and subsoil (b). Values were estimated on the basis of the concentration factors of Ag, As, Bi, Cd, Cu, W and Zn (A, Barroca Grande; B, S. Francisco de Assis; C, Rio).
Figure 5.
Spatial distribution of the modified degree of Contamination (mCd) for topsoil (a) and subsoil (b). Values were estimated on the basis of the concentration factors of Ag, As, Bi, Cd, Cu, W and Zn (A, Barroca Grande; B, S. Francisco de Assis; C, Rio).
The dams at Barroca Grande may pose a significant potential threat, due to the fine-grained nature of the materials, and their location with respect to the Casinhas stream that cross S. Francisco de Assis village. The X-Ray Diffraction (XRD) analysis of the impoundment material revealed the presence of scorodite, arsenopyrite, quartz, sphalerite, hematite, and muscovite. These tailings and impoundment material are metal-enriched at such a level, likely to be toxic to the ecosystem [
42]. The concentrations exceed the values defined for the 90th percentile of the South Portuguese Zone (As 157 mg·kg
−1; Cu 108 mg·kg
−1; Ni 62 mg·kg
−1; Pb 117 mg·kg
−1; Zn 134 mg·kg
−1), which is indicative of enrichment in trace metals. This Ag-As-Bi-Cd-Cu-W-Zn association is quite logical and is linked to the Panasqueira ore-paragenetic association. The highest
mCd values identified near Barroca Grande (A), São Francisco de Assis (B) and Rio (C) confirms that mechanical and chemical dispersion from Barroca Grande and Rio tailings and mud impoundments occurs. Most of the samples (90% or higher of total samples) collected from Barroca Grande (A), São Francisco de Assis (B) and Rio (C) villages, exhibit
mCd values >8.0, clearly indicting a very high degree of contamination.
The Panasqueira tailings impoundments have been and are affected by surface water flows (from heavy rainfall events) that have eroded the tailings from their original location and transported the materials downstream to residential areas (namely to S. Francisco de Assis). However, the superficial flat of the tailings have dried and are susceptible to wind erosion. The relative rates of water and wind erosion and transport, in Barroca Grande, São Francisco de Assis and Rio, suggest that wind processes have similar, and in many cases greater, impact on loss and local redistribution of soil in ecosystems than an eventual erosional soil enrichment.
3.3. Potential Ecological Risk Factor and Risk Index
The topsoil and subsoil samples results for individual element potential pollution factor (
EF) and potential ecological risk (
PERI) are presented in
Table 9. For both soil sample depths, As, Cu and Zn show a low potential ecological risk, with median values <40 (see
Table 2). Tungsten also exhibits a low risk in subsoil
, but a high risk in topsoil, while Bi and Ag show a moderate risk at both depths. Cadmium presents a very high ecological risk in topsoil and a high risk on subsoil. Using the median values, the topsoil risk factor is ranked as: Cd > W > Ag > Bi > As >> Cu > Zn, while for subsoil the ranking is: Cd > Ag > Bi > As > Cu = W > Zn. These results suggest a very high environmental risk, especially for Cd.
In order to estimate the global potential ecological risk in the study area, the PERI was computed. The median values classify soil samples at both depths with a very high risk (
Table 9). The cumulative analysis shows that the soil samples at both depths do not display a low risk index (<150), and only 7.4% of topsoil present a moderate risk index.
PERI classified 92.6% of topsoil samples as high to very high ecological risk. The same occurs for 61.2% of subsoil samples, which should be considered as an extensive hazard.
Figure 6 displays the
PERI spatial distribution. Topsoil has a wide area classified with a very high ecological risk, which is consistent with the wind direction, the water courses and the actual and previous exploration and beneficiation locations (
Figure 6a). Subsoil also presents a very high risk index in the same topsoil areas, but with smaller area expression. These results are consistent with those mapped by the individual elements, Ag, As, Bi, Cd, Cu, W and Zn (
Figure 2,
Figure 3 and
Figure 4), and the modified degree of contamination (
mCd;
Figure 5), with the same affected areas.
Table 9.
Statistical results of the single element potential pollution factor (EF) and potential ecological risk (PERI) for topsoil and subsoil samples.
Table 9.
Statistical results of the single element potential pollution factor (EF) and potential ecological risk (PERI) for topsoil and subsoil samples.
Var | AgEF | AsEF | BiEF | CdEF | CuEF | WEF | ZnEF | PERI |
---|
Topsoil | Min | 120 | 4 | 20 | 150 | 1 | 9 | 0 | 224 |
Mean | 120 | 23 | 127 | 893 | 3 | 573 | 1 | 1,740 |
Med | 70 | 15 | 60 | 600 | 2 | 255 | 1 | 1,020 |
Max | 1,190 | 338 | 4,300 | 9,000 | 21 | 8,571 | 3 | 23,353 |
Subsoil | Min | 30 | 4 | 10 | 125 | 2 | 2 | 0 | 173 |
Mean | 76 | 33 | 56 | 418 | 7 | 8 | 1 | 600 |
Med | 60 | 23 | 30 | 250 | 6 | 6 | 1 | 350 |
Max | 360 | 158 | 770 | 3,750 | 26 | 203 | 3 | 4,369 |
Figure 6.
Spatial distribution of Potential Ecological Risk Index (PERI) for topsoil (a) and subsoil (b) in the study area.
Figure 6.
Spatial distribution of Potential Ecological Risk Index (PERI) for topsoil (a) and subsoil (b) in the study area.
3.4. Human Health Risk Assessment
Both non-carcinogenic hazard (
HQ) and carcinogenic risk (Risk
pathway) of topsoil in the Panasqueira mine and surrounding area, through the different pathways (ingestion, dermal and inhalation), were estimated according to the human health risk assessment model [
28]. The cumulative hazard index (
HI) and total risk (Risk
total) were also characterized for multi-pathway routes in resident population.
The non-carcinogenic effects considered the most conservative exposure condition—children (1–6 years old). The potentially toxic elements defined in this study (Ag, Cd, Cu, W and Zn), apart from As, do not present a non-carcinogenic hazard for children in the Panasqueira area (maximum
HIchild-Ag,Cd,Cu,W,Zn ≤ 0.37). The As non-carcinogenic hazard median values, estimated for the different exposure routes were
HQing-As (2.75) >>
HQdrm-As (0.23) >>
HQinh-As (0.00).
HIchild-As values ranged between 0.78 and 69.50, with a median value of 2.98 ≈
HQing-As (2.70), due to the ingestion hazard quotient ranging between 0.72 and 64.10. These results (<1—safe level) indicate that there is a cause for concern for the non-cancer health effects for children living in the Panasqueira study area, mainly due to As oral ingestion, with
HQing-As showing median values above one (
Figure 7).
Figure 7.
Comparative boxplot distribution of the Non-carcinogenic Hazard Quotient for children of As for Ingestion (HQingAs), Inhalation (HQihnAs), Dermal contact (HQdrmAs) routes, the Cumulative Hazard Index (HIAs) and the sum of the Non-carcinogenic Hazard Quotient of the other defined PTEs (Ag, Cd, Cu, W and Zn) for Ingestion (HQingPTE(-As)), Inhalation (HQihnPTE(-As)), Dermal contact (HQdrmPTE(-As)) routes and the Cumulative Hazard Index (HIPTE(-As)) for topsoil samples (tailing samples were removed) in the Panasqueira area. The extremes and outliers were removed.
Figure 7.
Comparative boxplot distribution of the Non-carcinogenic Hazard Quotient for children of As for Ingestion (HQingAs), Inhalation (HQihnAs), Dermal contact (HQdrmAs) routes, the Cumulative Hazard Index (HIAs) and the sum of the Non-carcinogenic Hazard Quotient of the other defined PTEs (Ag, Cd, Cu, W and Zn) for Ingestion (HQingPTE(-As)), Inhalation (HQihnPTE(-As)), Dermal contact (HQdrmPTE(-As)) routes and the Cumulative Hazard Index (HIPTE(-As)) for topsoil samples (tailing samples were removed) in the Panasqueira area. The extremes and outliers were removed.
The carcinogenic risks adjusted both to children and adult were studied for the identified PTEs in this study. The median values of the different exposure routes and total risk were estimated. Their mean distribution were Risk
total (2.62 × 10
−4) ≈ Risk
As (2.62 × 10
−4) ≈ Risk
ing (2.39 × 10
−4) ≈ Risk
ingAs (2.39 × 10
−4) >> Risk
drm (2.26 × 10
−5) ≈ Risk
drmAs (2.26 × 10
−5) >> Risk
inh (1.32 × 10
−7) ≈ Risk
inhAs (1.32 × 10
−7) (
Figure 8). The high significance of As to cumulative carcinogenic elements risk is due to the very low risk displayed by the other elements. Furthermore, the representation of the exposure routes are distributed as: (a) Risk
ing ≈ Risk
ingAs ranged 3.97 × 10
−5–3.53 × 10
−3, with a median value of 1.52 × 10
−4; (b) Risk
drm ≈ Risk
drmAs ranged 3.76 × 10
−6–3.34 × 10
−4, with median value of 1.44 × 10
−5; and (c) Risk
inh ≈ Risk
inhAs ranged 2.21 × 10
−8–1.96 × 10
−6, with median value of 8.38 × 10
−8. The cumulative pathway Risk
total ≈ Risk
As ranged 4.35 × 10
−5–3.87 × 10
−3, with a median Risk
total of 1.66 × 10
−4 mostly due to Risk
ingAs. Intake of As may cause cancer in several human organs through ingestion, and lung and skin cancer through inhalation [
28]. Moreover, the Risk
ing represents 91.29% of cumulative risk from exposure routes, while Risk
drm (8.69%) and Risk
inh (0.05%) have less significance. Similar results were obtained in other studies [
77,
91,
92,
93,
94]. The cumulative median Risk
drm is > 1 × 10
−6 for all samples, of which 69.67% > 1 × 10
−5, and the cumulative exposure route median Risk
ing is > 1 × 10
−5 for all samples, being 72.13% > 1 × 10
−4. These results show that there is a very high As ingestion cancer risk. The samples with higher hazard are located in and around the villages of the study area, being the more representative results found nearby the large tailing piles and open air impoundments (
Figure 9). Although the median Risk
drm value is lower and between 1 × 10
−6 and 1 × 10
−4, and is considered as acceptable, it cannot, however, be negligible in the Panasqueira mining area, once all samples exceed the target value (1 × 10
−6). While the sum of exposure pathways is Risk
Cd = 1.09 × 10
−10 (<1 × 10
−6), it should be noted that it is the cumulative toxic metal and kidney index that is the main target for Cd toxicity [
26,
95].
Figure 8.
Comparative boxplot distribution of the Carcinogenic Risk adjusted to both children and adult for total PTEs (Ag, As, Cd, Cu, W and Zn) for Ingestion (Risking), Inhalation (Riskinh), Dermal contact (Riskdrm) routes, the Cumulative Risk (Risktotal) and the Carcinogenic Risk of As for Ingestion (RiskingAs), Inhalation (RiskinhAs), Dermal (RiskdrmAs) routes and the Cumulative Risk (RiskAs) for topsoil samples (tailing samples were removed) of the Panasqueira area. The extremes and outliers were removed.
Figure 8.
Comparative boxplot distribution of the Carcinogenic Risk adjusted to both children and adult for total PTEs (Ag, As, Cd, Cu, W and Zn) for Ingestion (Risking), Inhalation (Riskinh), Dermal contact (Riskdrm) routes, the Cumulative Risk (Risktotal) and the Carcinogenic Risk of As for Ingestion (RiskingAs), Inhalation (RiskinhAs), Dermal (RiskdrmAs) routes and the Cumulative Risk (RiskAs) for topsoil samples (tailing samples were removed) of the Panasqueira area. The extremes and outliers were removed.
As the risk analyses showed that ingestion is the most important pathway, further studies should be carried out in order to estimate the metal bioaccessibility in human receptors, such as a simplified
in vitro physiologically based extraction test (PBET), allowing the estimation of the percentage of the bioaccessible orally ingested fraction of each metal from soil samples [
77,
80,
96]. Additional investigations must be made on different sample media, such as vegetables, rhizosphere, irrigation and drinking water and street dust in the local villages under influence of mining activities (such as S. Francisco de Assis and Barroca) and also in non-exposed villages considered as reference areas (such as Casegas and Unhais-o-Velho). These studies will provide the necessary data and information for the determination of the influence of the Panasqueira mine activities, through more than 110 years, on the surrounding population, and should be complemented by biological studies [
97].
Figure 9.
Spatial distribution of the Carcinogenic risk adjusted for both children and adult considering (a) ingestion; (b) dermal contact; (c) inhalation routes; and (d) the cumulative Risk.
Figure 9.
Spatial distribution of the Carcinogenic risk adjusted for both children and adult considering (a) ingestion; (b) dermal contact; (c) inhalation routes; and (d) the cumulative Risk.