Next Article in Journal
Preliminary Assessment of Radiolysis for the Cooling Water System in the Rotating Target of SORGENTINA-RF
Next Article in Special Issue
Comparative Evaluation of Technologies at a Heavy Metal Contaminated Site: The Role of Feasibility Studies
Previous Article in Journal
Allometric, Growth, and Biomass Estimation Models for Acacia dealbata Link.: A Case Study in Serra da Estrela Natural Park (Portugal)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Soil Physicochemical Characteristics and As, Cu, Pb and Zn Contamination in Non-Active Mines at the Portuguese Sector of the Iberian Pyrite Belt

1
LEAF, Linking Landscape, Environment, Agriculture and Food Research Center, Associated Laboratory TERRA, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisbon, Portugal
2
DTCA, Departamento de Tecnologias e Ciências Aplicadas, Escola Superior Agrária do Instituto Politécnico de Beja, 7801-295 Beja, Portugal
3
ICT, Instituto de Ciências da Terra, Universidade de Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
4
CESAM & Departamento de Ambiente e Ordenamento, Universidade de Aveiro, 3810-193 Aveiro, Portugal
*
Author to whom correspondence should be addressed.
Environments 2022, 9(8), 105; https://doi.org/10.3390/environments9080105
Submission received: 3 July 2022 / Revised: 29 July 2022 / Accepted: 12 August 2022 / Published: 17 August 2022

Abstract

:
This study aimed to evaluate soil physicochemical characteristics (pH, electrical conductivity, organic matter, total N, and extractable P and K), and potentially toxic elements (As, Cu, Pb, and Zn), in non-active mines located in the Portuguese sector of the Iberian Pyrite Belt (IPB). A total of 70 sampling sites were surveyed at Aljustrel and Lousal, in areas already rehabilitated, and at São Domingos, where rehabilitation was only beginning. The soils at São Domingos were very heterogeneous, with extreme values for some properties (e.g., minimum soil pH 2.0 and maximum As concentration, 4382.8 mg kg−1 dry weight basis (DW)). Aljustrel was the site that presented soils with a higher total As, Cu, Pb, and Zn concentration (median values: 441.5, 545.9, 1396.8, and 316.5 mg kg−1 DW, respectively), above the soil quality guidelines values proposed by the Portuguese Environmental Agency (18, 230, 120, and 340 mg kg−1 DW, respectively). A principal component analysis identified the most relevant soil properties to explain the data variance, which were the soil pH and Pb total concentration, followed by Cu and Zn total concentrations, allowing a separation of Aljustrel from the other mines. Pearson correlation coefficients revealed very strong associations between Pb and As, markedly found at higher concentrations in São Domingos, whereas Aljustrel had an elevated concentration of As and Pb, but also of Cu and Zn. It is evident the risk that persists in the Aljustrel mine area, which was not alleviated by the “dig, dump, and cover” techniques that were implemented to rehabilitate the area.

1. Introduction

In the recently launched European Union (EU) Soil Strategy for 2030 [1], medium- and long-term objectives are assumed (to 2030 and 2050, respectively), aiming to achieve good soil health. In the view of this strategy, protection, sustainable use, and restoration of soil will became the norm by 2050 [1], which will require all Member States to take decisive action during this decade. This type of strategic document was expected for a long time, and it is essential to leverage specific actions in some Member States where the legal void concerning soil contamination and remediation persists. In fact, the lack of dedicated soil policies at the EU level has allowed some countries to delay the publication of legislation with contaminant thresholds and procedures for the rehabilitation of contaminated sites [1]. That is the case of Portugal, which has prepared a draft proposal in 2015, with the support of the Portuguese Environmental Agency (APA), but that was not officially published to the present date. Therefore, APA produced a Technical Guide with Reference Values for the Soil [2], based on the Ontario Soil Quality Guidelines, that, despite not having a legally compulsory value, are being used to assist processes of soil quality assessment, remediation, and confirmation of the results achieved with remediation measures.
The EU Soil Strategy for 2030 will undoubtably determine that all Member States independently develop methodologies for assessing the risks of soil contamination, identify contaminated sites, and create conditions to the remediation of sites that pose a risk to human health or the environment [1]. Mining areas are an example of historical or orphan contaminated sites, where the polluter(s) can no longer be held responsible for the costs of remediation. To tackle this problem, and to deal with multiple situations of pos-closure abandonment of mine sites, in some cases intensively exploited for hundreds of years, the Portuguese government attributed to EXMIN, presently EDM (Empresa de Desenvolvimento Mineiro, S.A.), the responsibility for the rehabilitation of these historically degraded mining areas [3]. Since 2001, EDM has developed efforts to assess soil contamination at these mine areas, prioritized interventions, established pre-remediation plans, took actions to remediate some of the sites, and implemented monitoring plans for the post-remediation [4]. So far, EDM has identified 199 abandoned and contaminated areas and is now working on their rehabilitation [4]. Nevertheless, interventions mainly consisted of the application of constructive techniques, which, although essential to manage the risk at the contaminated sites, may not be sufficient to restore the soil’s healthy status. Mourinha et al. [5] have recently published a review about the trace elements’ pollution in the Portuguese sector of the Iberian Pyrite Belt (IPB) and reported some of the remediation options adopted for these IPB mines. It is important to evaluate if, in the rehabilitated areas, the soil that was left behind may still need to have its characteristics improved to be considered a healthy soil, i.e., present good chemical, biological, and physical conditions, and, therefore, be able to provide a minimum number of ecosystem services [1,6].
The mine exploitation in the Portuguese sector of the IPB was very important in pre-Roman and Roman times, and again, with intensive mine works, during part of the 19th and 20th centuries, especially in Neves-Corvo, São Domingos, Aljustrel, Lousal, and Caveira (Alentejo region) [7,8,9]. Nowadays, the exploitation in Neves-Corvo (Cu, Zn, Pb, and Ag concentrates, Lundin Mining Company) and Aljustrel (Cu concentrates, Almina Company, which reactivated this mine) is on-going, but they are following the current best available practices to avoid environmental impacts [10]. This means that they are going to be made accountable for the soil remediation after their activity has ceased. The same was not true for some representative polymetallic sulfide mines in the Portuguese sector of the IPB, which have been abandoned for decades, i.e., São Domingos, Aljustrel, Lousal, and Caveira, and were a source of continuous environmental impact to soils, water, sediments, and biota [11,12,13,14,15]. The impacts in these mine areas have been extensively assessed and documented by different authors ([5,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26] refer to some in the Portuguese sector of the IPB).
Lousal and Ajustrel mines have their rehabilitation projects finished by EDM, and São Domingos is still ongoing [5]. However, in the surroundings of these mines, there are extensive areas exhibiting degraded and potentially contaminated soils, which are important to rehabilitate, preferably improving soil health and their ecosystem functions [1,6]. It is important to consider their constraints, namely the properties that hinder soil fertility and the concentrations of the potentially toxic elements (PTEs) in those soils [27]. As a matter of fact, the great majority of the studies which assessed the abandoned mines in the Portuguese sector of the IPB, were mainly focused on their geological and geochemical aspects, and in the environmental impacts regarding PTEs concentrations [13,14,16,19,20,21,22,23,24] and less focused in the soils characteristics, which condition their recovery, namely using phytotechnologies. Nevertheless, this type of assessment has been done in some studies, with a full characterization of the soils to be used in lab-scale or pot experiments [28,29,30,31,32,33,34], but only with soil from a specific site, or considering composite samples.
This study aimed to evaluate soil physicochemical characteristics (pH, electrical conductivity (EC), organic matter (OM) content, and macronutrients content, nitrogen (N), phosphorous (P), and potassium (K)), and PTEs contamination associated with the mining activities in the Portuguese sector of the IPB (As, Cu, Pb, and Zn, which are those commonly found in excessive concentrations). The study will assess the most representative non-active mines (Lousal and São Domingos), or areas outside the concession of active mines (Aljustrel). The total As, Cu, Pb, and Zn concentrations will be further compared with the soil quality guidelines values proposed by APA [2], that, in the absence of a specific national regulatory framework, may be used to assess the likelihood of risk of soil contamination, and the need for further actions to confine and manage the risk of the contaminated site, or to take active remediation measures [1]. Using multivariate exploratory techniques, the data will be statistically analyzed to evidence the soil properties that are more relevant to explain the data variance, and to potentially discriminate between mining areas and their “signature” properties.
This field survey was performed in the scope of the project Life No_Waste (Management of biomass ash and organic waste in the recovery of degraded soils: a pilot project set in Portugal—LIFE14 ENV/PT/000369), surveying a total of 70 sites in three mining areas: Aljustrel, Lousal, and São Domingos. To our best knowledge, this is the most complete soil survey in mines of the Portuguese sector of the IPB, gathering not only contaminants concentration, but also soil properties important to achieve soil health, namely, soil pH, organic matter, and nutrients content.

2. Materials and Methods

2.1. Site Description

The IPB, located in the SW of the Iberian Peninsula, forms an arch extending from Grândola (Portugal) to Seville (Spain), which is approximately 240 km long and 35 km wide [7] (Figure 1).
It was formed 350 million years ago, connected to active hydrothermal volcanism that led to the formation of a Volcano–Sedimentary Complex [7], and represents one of the most important volcanogenic massive sulfide districts in the world [7,8,9,35]. The massive sulfide ore is composed mainly by pyrite (FeS2), approximately 95%, with variable amounts of chalcopyrite (CuFeS2), sphalerite (ZnS), galena (PbS), and arsenopyrite (FeAsS), while a Cu-rich stock zone holds most of the remaining mineralization [35].
The Aljustrel mine is one of the greatest sulfide deposits of the IPB, containing six mineral masses rich in Cu and Zn (Feitais, Estação, Algares, Moinho, S. João, and Gavião) [6,20,21]. The major environmental impacts at this mine were caused by the large volume of tailings dispersed in the area, which affected, for a long period, the soils and the hydrological system in the surrounding area [15,36,37,38,39,40]. The Aljustrel mine’s historically contaminated area was rehabilitated by EDM, as reported by Mourinha et al. [5] (Figure 2).
The Lousal mine is located in the NW part of the Portuguese sector of the IPB, in a lineament of the volcano–sedimentary complex which also includes the old Caveira pyrite mine [22]. The deposit was explored, mainly for pyrite, until 1988, when it was closed due to the low grades of Cu and Zn in the mined ores [41]. The rehabilitation of the mine area begun less than 10 years after the closure of the mine, by the Fréderic Velge Foudation and the Grândola Municipality, which made efforts to preserve the mine structures to be used for science and touristic purposes [5]. Later, the EDM also rehabilitated a part of the mine area, mainly to treat the acid waters that formed in the milled ore deposited and in the open pit, avoiding their impact in the surrounding hydrological system [5,22,23,24] (Figure 3).
São Domingos is, perhaps, the most representative mine district in the Portuguese sector of the IPB, with a big pit (122 m deep), left open after the extraction works, which is now filled with acid waters [18]. São Domingos was extensively exploited from 1857 to 1966, causing considerable impacts in soils, water, and sediments, which have been assessed by different authors [12,13,18,19,25,26,42,43]. The mined raw material was crushed in a mill, located near the open cast pit, and transported to the sulfur factory at Achada do Gamo (~3 km south), where it was smelted to obtain high-level grades of Cu ore and sulfur products [16]. These works generated approximately 750,000 tons of highly heterogeneous wastes, including Roman and modern slags, smelting ashes, and pyrite-rich waste dumps [18,19], which had a considerably negative impact observed along the São Domingos stream valley and in the sulfur factory at Achada do Gamo, the major area of the waste deposits [11] (Figure 4).

2.2. Soil Sampling

Soils from the Aljustrel, Lousal, and São Domingos mines’ affected areas were surveyed in July 2016, choosing areas with greater and lesser impact on the mining activity, in a total of 70 sampling sites (14 at Aljustrel, 12 at Lousal, and 44 at São Domingos). The selection of the sampling sites had the collaboration of EDM and took into consideration the fact that Aljustrel and Lousal were already rehabilitated. Therefore, the sites that were chosen covered areas which were still potentially affected by AMD, or may need further intervention, to complement the rehabilitation already performed [5].

2.3. Soil Physicochemical Characterization

Soil samples were air-dried and sieved using a 2 mm non-metallic sieve, to collect the <2 mm fraction to further analysis. Soils’ particle-size distribution were determined using the pipet method [44] and their texture class classification was performed according to the United States Department of Agriculture (USDA) classification [45]. Physicochemical characterization of the samples was performed with parameters typically used to ascertain general soil fertility characteristics: pH (H2O), electrical conductivity (EC), total nitrogen (N), and the extractable fraction of phosphorous (P) and potassium (P). Analyses were performed according to methodologies described by Alvarenga et al. [31]: soil pH (H2O), was determined in a soil to deionized water suspension of 1:2.5 (w/v); EC was determined in a soil to deionized water suspension of 1:5 (w/v); total N was analyzed by the Kjeldahl method (NKjeldahl); total oxidizable organic carbon was determined according to the Walkley and Black wet-oxidation method [46], and converted to organic matter content (OM) by multiplying by a factor of 1.72; extractable P and K were determined using the Egner–Riehm extraction procedure, with an ammonium lactate (0.1 M) in acetic acid (0.4 M) solution, with a pH in the range 3.65–3.75, sufficient to be used in acid soils [47].

2.4. Pseudo-Total As, Cu, Pb and Z Concentrations

Pseudo-total metal concentrations (As, Cu, Pb, and Zn) were determined by inductively coupled plasma-mass spectrometry (ICP-MS), using a Thermo X Series spectrophotometer. Samples were digested with aqua regia, according to ISO 11,466 [48]. Three independent replicates were performed for each sample and blank, and control standards were measured in parallel. ICP-MS element analysis was performed following rigorous quality control proceedings, and the detection limits were: 8 μg kg−1 for As, and 1 μg kg−1 for Cu, Pb, and Zn.

2.5. Statistical Treatment of the Data

The Spearman rank correlation coefficients (R) were calculated to evaluate the associations between the soil physicochemical properties with their pseudo-total element concentrations (n = 65). Three levels of significance were used: p < 0.001, p < 0.01, and p < 0.05.
A Principal Component Analysis (PCA) was performed to identify the most meaningful parameters which describe the whole data set. This multivariate exploratory technique allows the association between variables, reducing the dimension of the data matrix, transforming the original variables into the same number of uncorrelated variables, called principal components (PCs). The eigenvalues of the PCs are a measure of their variance, the participation of the original variables in the PCs is given by the loadings, and the individual transformed observations are called scores [33].
All statistical analyses were performed with the STATISTICA 7.0 (Software™ Inc., Tulsa, OK, USA, 2004).

3. Results and Discussion

3.1. Soil Physicochemical Characteristics

The soils presented a high variability in their texture classification, given that they cover different areas, but they were mainly characterized by a coarse granulometry (e.g., sand, loamy sand, sandy loam, 64% of the samples), with a lower number of samples with medium to fine granulometry (e.g., sandy clay loam, loam, clay loam, and clay, USDA texture classes) [45], with a median fraction of clay, silt, and sand, in the different sampling areas, of 15, 27, and 55%, respectively.
Considering the soils’ physicochemical properties that may hinder their fertility (Figure 5), it is possible to have an idea about their range of variation. In the São Domingos mine, the area covered was larger and, as the rehabilitation was in its initial phase, the sampled soils evidenced extreme characteristics, a higher number of outliers, and, for some properties, higher or lower extreme values than in the other mines. That was the case of the soil pH (Figure 5a), with soils’ pH values < 2 in some locations of São Domingos, despite its median value being higher that that found for Lousal and Aljustrel soils, and with pH < 4 in both mines, considered as very acid soils [49]. However, that higher dispersion of values in the soils sampled at São Domingos were also felt for the soil EC, soil OM content, total N, and extractable P (Figure 5b–e), being a consequence not only of the heterogeneity of the area, but also of the fact that the rehabilitation program was only beginning. The rehabilitation operations that were performed by EDM at Lousal and Aljustrel [5], which have removed and confined the most contaminated soils/wastes, had the benefit of making the whole area more homogeneous, which would facilitate further interventions (e.g., the implementation of phytotechnologies). In addition, some of the soils samples at São Domingos could be considered “natural soils”, not submitted to the effect of AMD, which could have influenced this higher dispersion of values.
The soils were also characterized by very low OM content, with median values < 1% OM (g/100 g DW) for the three mines (Figure 5c) [49], and low macronutrients content, especially with total N and extractable P (Figure 5d,e). These values of extractable P in soils were considered as very low for this macronutrient in the three mines [49]. For extractable K, the values were different, depending on the mine, with concentrations typical of soils with low extractable K concentrations, at Aljustrel and Lousal, and with higher concentrations at São Domingos (Figure 5f).
Previous studies have alerted to the low quality of soils in some of these sites. Santos et al. [29] evaluated some of these properties in the São Domingos mine area and referred that, in some cases, the soils were developed on mining wastes composed by gossanous materials and host rocks, with very low pH values (4.53), low OM content (1.687%, w/w), as well as residual concentration for some macronutrients, especially N and P plant (0.004% N (w/w) and 2.23 mg kg−1 for extractable P, dry weight, DW) [29].
Considering the abovementioned soil characteristics, and if the chemical stabilization of the soils is intended, the main characteristic to be corrected would be the soil pH. The use of lime-based amendments (e.g., calcium carbonate, calcium oxide, calcium hydroxide) [31,50], or waste-derived materials with alkaline pH (e.g., biomass ashes), would, therefore, be crucial to increase the soils’ pH and, potentially, reduce the PTEs extractability [5,51]. Taking this into consideration, the agronomical recommendations to increase the soil pH(H2O) to 6.5, for soils with low OM content (between 0.5–1%), and soil pH(H2O) < 4.5, the lime application ratio should be 4–5 t CaCO3 ha−1, depending on the soil texture. That target pH can be considered ideal to have lower PTEs solubility and mobility and, very important, to allow seed germination and plant survival [49,52].
However, the lime-based amendments are, usually, poor in OM and nutrients [51]. Therefore, for a successful natural attenuation, or eventually, if a phytomanagement strategy is intended, organic amendments, rich in OM, N, and P, should also be considered [27,31,32,53,54].

3.2. Potentially Toxic Elements Concentrations

The upper layers of soil affected by mining activities can contain elevated concentrations of PTEs, which depend on the ore that was mined. In the case of the mines located at the Portuguese sector of the IPB, the PTEs that are often found in elevated concentrations are As, Cu, Pb, and Zn, and these were the ones which were assessed in this study.
As for the pseudo-total concentrations for As (Figure 6a), and Pb (Figure 6c) (aqua-regia extractable), their values surpass the recommended limit values for soils intended for industrial use, considering the APA’s Technical Guide [2]. Median concentration for As in the three mines were, in fact, several times higher than that limit value (18 mg As kg−1 DM), with a higher dispersion of values for the São Domingos mine, higher upper limit concentrations, and a higher median (Figure 6a). High concentration for total As in the São Domingos mine were also reported by Freitas et al. [18], 1291.0 mg kg−1 DW; Santos et al. [29], 2600 mg kg−1, and Alvarenga et al. [26], 1956 mg kg−1. However, some of these authors have alerted to the fact that, despite the high total concentrations for As, the extractable fractions (e.g., with CaCl2 0.01 M, which can be considered surrogate measures of their bioavailable fraction), are, usually, very low, <1% [26].
As for Pb, the limit value proposed in the Technical Guide of APA [2] (120 mg Pb kg−1 DM) was surpassed by more samples in Aljustrel, followed by São Domingos, than in the Lousal mine. In fact, the samples at Lousal presented low median Pb content, lower than the proposed limit value. The samples from Aljustrel can be considered the most problematic regarding Pb contamination (Figure 6c). Other authors have also reported high concentration for total Pb in Aljustrel, such as Alvarenga et al. [28], 3500 for Pb mg kg−1, and Candeias et al. [20,21], which reported values as high as 20,000 mg Pb kg−1 of soil. Extremely high Pb concentrations were also reported by other authors in the São Domingos mine area, such as Freitas et al. [18] 12,217.5 mg Pb kg−1 DW, Santos et al. [29], 7300 mg kg−1, and Alvarenga et al. [25,26], up to 10,795 mg kg−1.
Copper is also more problematic in the Aljustrel’s soil samples (Figure 6b), with a median concentration (546 mg Cu kg−1 DW) above the threshold value recommended (230 mg Cu kg−1 DW), more than doubling its value, and with samples that can go as high as 4383 mg Cu kg−1 DW (Figure 6b). These values were in accordance with others reported in previous studies by Candeias et al. [20,21], that reported concentrations for Cu up to 5414 mg kg−1, or by Alvarenga et al. [28], reporting 1800 mg Cu kg−1 of soil. The same was not true for the samples collected at Lousal and São Domingos, which were characterized by, generally, lower Cu pseudo-total concentrations (Figure 6b), below the threshold of the recommended value in most of the samples. That general high concentration for Cu can, therefore, be considered a distinctive characteristic of the soils at Aljustrel. This is a fact that is important to be taken into consideration, since, in Aljustrel, most of the samples were collected in soils which were already submitted to a rehabilitation project [5], which means that the soils that were considered rehabilitated still have a high Cu total concentration, which means that the mobility/bioavailability should be thoroughly monitored and controlled.
That same pattern was followed for Zn (Figure 6d): higher concentration in Aljustrel than in Lousal and São Domingos. However, in this case, 75% of the samples at Lousal and São Domingos had pseudo-total Zn concentration below the recommended limit value (340 mg Zn kg−1 DW) [2], and even the median value for Zn in Aljustrel (316 mg Zn kg−1 DW) was below that threshold limit value. Therefore, it is fair to say that, even though some samples that were found had an elevated concentration of Zn (979 and 840 mg Zn kg−1 DW as upper limits at Aljustrel and São Domingos, respectively), Zn was not the most problematic PTE found in these mine areas.
The concentration of PTEs at these abandoned mines in the IPB have been thoroughly studied by different authors, and Mourinha et al. [5] have recently published a review with indicative ranges for the concentrations of As, Cu, Pb, and Zn found in different mines in the IPB, both in Portugal and in Spain.

3.3. Evaluation of the Association between the Variables—Pearson Correlations

To have an idea about the associations between the soil physicochemical characteristics which were analyzed and the As, Cu, Pb, and Zn total concentration, Pearson correlation coefficients (r) were calculated (Table 1). Three levels of significance were used to classify the correlation as: moderate (p < 0.05), high (p < 0.01), or very high (p < 0.001), evidencing the strength of the associations.
The organic matter content was highly positively correlated with the extractable P and K concentrations (p < 0.001), which was not surprising, given the importance of soil OM as an indicator of soil quality and it was used to assess its fertility. Soil OM also had a high positive correlation with soil pH (p < 0.01), meaning that less acidic soils, with higher pH values, were also those with the higher OM content.
Another interesting correlation which was found, between general physicochemical properties, was the very high negative correlation found between soil pH and its EC (p < 0.001) (Table 1); it is a fact, for a considerable number of metal salts, their solubility increases as the soil pH decreases [55]. This case was not an exception: considering all the results (n = 65), higher soluble salt concentrations (higher EC values) were found in the most acidic soils.
As for the PTEs correlations, very strong associations were found between total Pb and total As, and between total Cu and total Zn (p < 0.001), meaning that they happen to appear together as soils’ contaminants in these mines. These associations corroborate the results discussed in Section 3.2: As and Pb were markedly found at higher concentrations in São Domingos, which, on the other hand, had soils not so markedly contaminated by Cu and Zn. Regarding Aljustrel, not only As and Pb were found in elevated concentrations, but also Cu and Zn, strengthening these associations regarding PTEs affecting soils in the Portuguese sector of the IPB.

3.4. Principal Component Analysis

PCA was carried out on 10 variables (pH, EC, OM content, NKjeldahl, extractable P, extractable K, and As, Cu, Pb, and Zn pseudo-total concentrations in the soils), with the results obtained for the 65 soil samples collected in the three mines, to identify those which can be considered key parameters to explain data variability.
The original 10 variables could be reduced to four principal components, explaining about 70% of the total variance (Table 2). The first component (PC1) explained 24.6% of the variance of the original variables, and pH was the only variable with a high positive factor loading on PC1 (>0.7), which means that, by itself, pH was mainly responsible for the samples’ variability. The second component (PC2) explained 17.6% of the variance, whereas the third component explained 15.6%. Lead, Cu, and Zn pseudo-total concentrations had significant factor loadings on these principal components, Pb with PC2, with a negative correlation coefficient, and Cu and Zn on PC3, both with positive correlations. A fourth component (PC4) was necessary to explain a minor part of the variance (12.0%), and only one soil property correlated significantly with this component, with a negative value, which was soil NKjeldahl content, the only characteristic related with soil fertility which was ascribed to one PC. Some of the properties were not significantly correlated to the principal components, meaning that they are not responsible for the variability observed in the data set: EC, OM, extractable P and K, and As pseudo-total concentration.
Figure 7 shows the combined plot of scores and loadings on the first three main principal components: PC1 versus PC2 (Figure 7a), PC1 versus PC3 (Figure 7b), and PC2 versus PC3 (Figure 8). The most important parameters for the definition of each component are shown at the edge of each axis, indicating the direction in which the value of the parameter increases.
Results from the PCA emphasizes the importance of soil pH, and, to a lesser extent, Cu, Pb, and Zn pseudo-total concentration, to explain the variability of the soil properties in these mine areas. These properties were just associated between themselves, and not with the other characterized soil properties, except for pH, since highly acidic soils, subject to AMD, are also those with the highest contamination level. The PCA did not allow a clear clustering of samples from the three mines, which was somehow expectable, given the similarity of these mining areas, all located in the IPB. In fact, it was possible to confirm some of the specificities of each mine site.
Considering the PC1 versus PC2 plan (Figure 7a), the samples located in the left quadrants were those with lower pH values, evidencing the more acidic soils (64% of Aljustrel’s samples). Distant from these samples, in the positive part of PC1, it is possible to find two samples from São Domingos with unusually high pH values: SD8 and SD5, with a soil pH of 7.02 and 6.22, respectively, but also, a great part of the São Domingos samples (59%), corresponding to acid soils, but not as acidic as in Aljustrel (see also Figure 5a). On the other hand, the samples located downward from the origin axis of PC2, on its negative part, were those which presented higher Pb pseudo-total concentrations. There, it is possible to find the majority of the Aljustrel’s samples (57%), and a smaller number of samples from São Domingos (38%) and Lousal (17%).
The PC1 vs. PC3 map (Figure 7b) evidenced that the majority of the Aljustrel’s samples, located in the positive part of PC3 with a positive correlation with Cu and Zn pseudo-total concentrations, corresponded to the locations with higher concentration on those PTEs. PC3 was not able to discriminate the samples from São Domingos, which were almost all concentrated near the origin of that axis, except for SD17, in the positive part of that chart, with an exceptionally high total concentration for both elements (3383.1 mg Cu kg−1 and 8045 mg Zn kg−1).
The plan of PC2 versus PC3 plan (Figure 8) emphasized the samples with, simultaneously, the highest Pb, Cu, and Zn pseudo-total concentration, which were those located in the upper-left quadrant, where, again, it is possible to find some samples from Aljustrel (A1, A3, A4, A5, A8, A9, and A13). Since some of these samples were collected in an area that was already rehabilitated by EDM using constructive techniques [5], it is fair to say that further intervention actions should be considered, to control not only the low pH values found, but also the high pseudo-total PTEs concentrations at the site.

4. Conclusions

The PCA identified the most relevant soil properties in order to explain data variance in the mine areas located in the Portuguese sector of the IPB, which were soil pH (24.6% of the variance, with a high positive factor loading on PC1), Pb total concentration (17.6% of the variance, with a high negative factor loading with PC2), followed by Cu and Zn total concentrations (15.6% of the variance, high positive factor loadings on PC3), allowing a separation of the Aljustrel’s samples from those collected in the other mines. The Pearson correlation coefficients revealed very strong associations between Pb and As, a consequence of the co-existence of galena and arsenopyrite, and between Cu and Zn (p < 0.001), meaning that they happen to appear together as soils’ contaminants in these mines. In fact, As and Pb were markedly found at a higher concentration in São Domingos, whereas Aljustrel had an elevated concentration for As and Pb, and for Cu and Zn. Therefore, it is evident the risk that persists in the Aljustrel mine area, with high concentration of As, Cu, and Pb, low pH, and degraded nutritional status, which was not alleviated by the techniques that were implemented to rehabilitate the area, which consisted mainly in “digging, dumping, and covering”.
Further remediation measures in these mining areas are, therefore, recommended. That is true even for the mines which were already rehabilitated, especially Aljustrel. Nevertheless, it is important to emphasize that these values should also be compared with background concentrations, because these soils are naturally enriched in some of these PTEs. Comparison with the APA standards can provide some reference but should not be the only criterion to assess the need to decontaminate, nor the level of concentration that should be targeted, and the risk should also consider PTEs availability to plants and soil organisms.
One possibility to control the risk, once the complete and rapid removal of the contaminated soils is expensive, disruptive, and with high environmental and social impacts, is by rendering the PTEs less mobile and bioavailable, by the application of assisted phytoremediation, selecting proper soil amendments and native plant species, strategies that are being recommended to these types of derelict environments.

Author Contributions

Conceptualization, P.A.; methodology and validation, N.C. and S.M.R.; formal analysis, P.A., P.P. and S.M.R.; investigation, C.M., N.C., S.M.R. and P.A.; resources, P.A., P.P. and S.M.R.; writing—original draft preparation, P.A.; writing—review and editing, P.A., P.P. and S.M.R.; visualization, P.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the project Life No_Waste-LIFE14 ENV/PT/000369—“Management of biomass ash and organic waste in the recovery of degraded soils: a pilot project set in Portugal”, and by national funds provided by FCT—Fundação para a Ciência e a Tecnologia, I.P., project UID/AGR/04129/2020 (LEAF), and ICT project (UIDB/04683/2020) with the reference POCI-01-0145-FEDER-007690.

Acknowledgments

The authors acknowledge the important support given by EDM (Empresa de Desenvolvimento Mineiro, S.A.), and their collaborators Edgar Carvalho and Catarina Diamantino, which collaborated in the sampling campaign. The authors are also thankful to Marta Baioneta Martins for the map preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Commission. EU Soil Strategy for 2030. Communication from the Commission to the European Parliament, the Council, the European Economic and Social Committee and the Committee of the Regions. Brussels, 17 November 2021. COM(2021) 699 Final. 2021. Available online: https://ec.europa.eu/info/law/better-regulation/have-your-say/initiatives/12634-Healthy-soils-new-EU-soil-strategy_en (accessed on 13 January 2022).
  2. Agência Portuguesa do Ambiente. Solos Contaminados–Guia Técnico. Valores de Referência Para o Solo. 2019. Available online: https://apambiente.pt/avaliacao-e-gestao-ambiental/guias-tecnicos-0 (accessed on 29 June 2021).
  3. Decree-Law No. 198-A/2001, 6 July 2001. Estabelece o Regime Jurídico de Concessão do Exercício da Actividade de Recuperação Ambiental das Áreas Mineiras Degradadas. Available online: https://data.dre.pt/eli/dec-lei/198-a/2001/07/06/p/dre/pt/html (accessed on 30 June 2021).
  4. EDM (Empresa de Desenvolvimento Mieiro, S.A.). Inventariação de Áreas Mineiras. 2017. Available online: https://edm.pt/area-ambiental/inventariacao-de-areas-mineiras/ (accessed on 28 June 2021).
  5. Mourinha, C.; Palma, P.; Alexandre, C.; Cruz, N.; Rodrigues, S.M.; Alvarenga, P. Potentially Toxic Elements’ Contamination of Soils Affected by Mining Activities in the Portuguese Sector of the Iberian Pyrite Belt and Optional Remediation Actions: A Review. Environments 2022, 9, 11. [Google Scholar] [CrossRef]
  6. Adhikari, K.; Hartemink, A.E. Linking soils to ecosystem services—A global review. Geoderma 2016, 262, 101–111. [Google Scholar] [CrossRef]
  7. Leistel, J.M.; Marcoux, E.; Thiéblemont, D.; Quesada, C.; Sánchez, A.; Almodóvar, G.R.; Pascual, E.; Sáez, R. The Volcanic-Hosted Massive Sulphide Deposits of the Iberian Pyrite Belt. Miner. Depos. 1997, 33, 2–30. [Google Scholar] [CrossRef]
  8. Tornos, F. Environment of Formation and Styles of Volcanogenic Massive Sulfides: The Iberian Pyrite Belt. Ore Geol. Rev. 2006, 28, 259–307. [Google Scholar] [CrossRef]
  9. Tornos, F.; Pamo, E.; España, J. The Iberian Pyrite Belt. In Contextos Geológicos Españoles: Una Aproximación al Patrimonio Geológico de Relevancia Internacional; Instituto Geológico y Minero: Madrid, Spain, 2008; pp. 56–64. ISBN 978-84-7840-754-5. [Google Scholar]
  10. Law, N. 54/2015, 22 June 2015. Bases do Regime Jurídico da Revelação e do Aproveitamento dos Recursos Geológicos Existentes no Território Nacional, Incluindo os Localizados no Espaço Marítimo Nacional. Available online: https://data.dre.pt/eli/lei/54/2015/06/22/p/dre/pt/html. (accessed on 30 June 2021).
  11. Oliveira, J.M.S.; Farinha, J.; Matos, J.; Paula, Á.; Rosa, C.; Machado, M.J.C.; Daniel, F.S.; Martins, L.; Leite, M.R.M. Diagnóstico Ambiental das Principais Áreas Mineiras Degradadas do País. Bol. Minas 2002, 39, 67–85. [Google Scholar]
  12. Matos, J.X.; Martins, L.P. Reabilitaçao ambiental de áreas mineiras do sector português da Faixa Piritosa Ibérica: Estado da arte e perspectivas futuras. Bol. Geol. Y Min. 2006, 117, 289–304. [Google Scholar]
  13. Álvarez-Valero, A.M.; Pérez-López, R.; Matos, J.; Capitán, M.A.; Nieto, J.M.; Sáez, R.; Delgado, J.; Caraballo, M. Potential Environmental Impact at São Domingos Mining District (Iberian Pyrite Belt, SW Iberian Peninsula): Evidence from a Chemical and Mineralogical Characterization. Environ. Geol. 2008, 55, 1797–1809. [Google Scholar] [CrossRef]
  14. Pérez-López, R.; Delgado, J.; Nieto, J.M.; Márquez-García, B. Rare Earth Element Geochemistry of Sulphide Weathering in the São Domingos Mine Area (Iberian Pyrite Belt): A Proxy for Fluid–Rock Interaction and Ancient Mining Pollution. Chem. Geol. 2010, 276, 29–40. [Google Scholar] [CrossRef]
  15. Alvarenga, P.; Guerreiro, N.; Simões, I.; Imaginário, M.J.; Palma, P. Assessment of the Environmental Impact of Acid Mine Drainage on Surface Water, Stream Sediments, and Macrophytes Using a Battery of Chemical and Ecotoxicological Indicators. Water 2021, 13, 1436. [Google Scholar] [CrossRef]
  16. Quental, L.; Brito, M.G.; Sousa, A.J.; Abreu, M.M.; Batista, M.J.; Oliveira, V.; Vairinho, M.; Tavares, T. Utilização de imagens hiperespectrais na avaliação da contaminação mineira em S. Domingos, Faixa Piritosa, Alentejo. In Proceedings of the VI Congresso Nacional de Geologia, Monte de Caparica, Portugal, 4–6 June 2003. [Google Scholar]
  17. Abreu, M.M.; Batista, M.J.; Magalhães, M.C.F.; Matos, J.X. Acid mine drainage in the Portuguese Iberian Pyrite Belt. In Mine Drainage and Related Problems; Brock, C.R., Ed.; Nova Science Pub Inc.: New York, NY, USA, 2010; pp. 71–118. [Google Scholar]
  18. Freitas, H.; Prasad, M.N.V.; Pratas, J. Plant Community Tolerant to Trace Elements Growing on the Degraded Soils of São Domingos Mine in the South East of Portugal: Environmental Implications. Environ. Int. 2004, 30, 65–72. [Google Scholar] [CrossRef]
  19. Pérez-López, R.; Álvarez-Valero, A.M.; Nieto, J.M.; Sáez, R.; Matos, J.X. Use of Sequential Extraction Procedure for Assessing the Environmental Impact at Regional Scale of the São Domingos Mine (Iberian Pyrite Belt). Appl. Geochem. 2008, 23, 3452–3463. [Google Scholar] [CrossRef]
  20. Candeias, C.; Ferreira da Silva, E.; Salgueiro, A.R.; Pereira, H.G.; Reis, A.P.; Patinha, C.; Matos, J.X.; Ávila, P.H. Assessment of Soil Contamination by Potentially Toxic Elements in the Aljustrel Mining Area in Order to Implement Soil Reclamation Strategies. Land Degrad. Dev. 2011, 22, 565–585. [Google Scholar] [CrossRef]
  21. Candeias, C.; Ferreira da Silva, E.; Salgueiro, A.R.; Pereira, H.G.; Reis, A.P.; Patinha, C.; Matos, J.X.; Ávila, P.H. The Use of Multivariate Statistical Analysis of Geochemical Data for Assessing the Spatial Distribution of Soil Contamination by Potentially Toxic Elements in the Aljustrel Mining Area (Iberian Pyrite Belt, Portugal). Environ. Earth Sci. 2011, 62, 1461–1479. [Google Scholar] [CrossRef]
  22. Ferreira da Silva, E.; Patinha, C.; Reis, P.; Fonseca, E.C.; Matos, J.X.; Barrosinho, J.; Oliveira, J.M.S. Interaction of Acid Mine Drainage with Waters and Sediments at the Corona Stream, Lousal Mine (Iberian Pyrite Belt, Southern Portugal). Environ. Geol. 2006, 50, 1001–1013. [Google Scholar] [CrossRef]
  23. Ferreira da Silva, E.; Fonseca, E.C.; Matos, J.X.; Patinha, C.; Reis, P.; Oliveira, J.M.S. The Effect of Unconfined Mine Tailings on the Geochemistry of Soils, Sediments and Surface Waters of the Lousal Area (Iberian Pyrite Belt, Southern Portugal). Land Degrad. Dev. 2005, 16, 213–228. [Google Scholar] [CrossRef]
  24. Luís, A.T.; Teixeira, P.; Almeida, S.F.P.; Matos, J.X.; da Silva, E.F. Environmental Impact of Mining Activities in the Lousal Area (Portugal): Chemical and Diatom Characterization of Metal-Contaminated Stream Sediments and Surface Water of Corona Stream. Sci. Total Environ. 2011, 409, 4312–4325. [Google Scholar] [CrossRef]
  25. Alvarenga, P.; Palma, P.; de Varennes, A.; Cunha-Queda, A.C. A Contribution towards the Risk Assessment of Soils from the São Domingos Mine (Portugal): Chemical, Microbial and Ecotoxicological Indicators. Environ. Pollut. 2012, 161, 50–56. [Google Scholar] [CrossRef]
  26. Alvarenga, P.; Laneiro, C.; Palma, P.; de Varennes, A.; Cunha-Queda, C. A Study on As, Cu, Pb and Zn (Bio)Availability in an Abandoned Mine Area (São Domingos, Portugal) Using Chemical and Ecotoxicological Tools. Environ. Sci. Pollut. Res. 2013, 20, 6539–6550. [Google Scholar] [CrossRef]
  27. Kidd, P.; Mench, M.; Álvarez-López, V.; Bert, V.; Dimitriou, I.; Friesl-Hanl, W.; Herzig, R.; Olga Janssen, J.; Kolbas, A.; Müller, I.; et al. Agronomic Practices for Improving Gentle Remediation of Trace Element-Contaminated Soils. Int. J. Phytoremediat. 2015, 17, 1005–1037. [Google Scholar] [CrossRef]
  28. Alvarenga, P.M.; Araújo, M.F.; Silva, J.A.L. Elemental Uptake and Root-Leaves Transfer in Cistus Ladanifer L. Growing in a Contaminated Pyrite Mining Area (Aljustrel-Portugal). Water Air Soil Pollut. 2004, 152, 81–96. [Google Scholar] [CrossRef]
  29. Santos, E.S.; Abreu, M.M.; Nabais, C.; Saraiva, J.A. Trace Elements and Activity of Antioxidative Enzymes in Cistus ladanifer L. Growing on an Abandoned Mine Area. Ecotoxicology 2009, 18, 860–868. [Google Scholar] [CrossRef]
  30. Santos, E.S.; Magalhães, M.C.F.; Abreu, M.M.; Macías, F. Effects of Organic/Inorganic Amendments on Trace Elements Dispersion by Leachates from Sulfide-Containing Tailings of the São Domingos Mine, Portugal. Time Evaluation. Geoderma 2014, 226–227, 188–203. [Google Scholar] [CrossRef]
  31. Alvarenga, P.; Gonçalves, A.P.; Fernandes, R.M.; de Varennes, A.; Vallini, G.; Duarte, E.; Cunha-Queda, A.C. Evaluation of Composts and Liming Materials in the Phytostabilization of a Mine Soil Using Perennial Ryegrass. Sci. Total Environ. 2008, 406, 43–56. [Google Scholar] [CrossRef] [PubMed]
  32. Alvarenga, P.; Gonçalves, A.P.; Fernandes, R.M.; de Varennes, A.; Duarte, E.; Cunha-Queda, C.A.A.; Vallini, G. Reclamation of a Mine Contaminated Soil Using Biologically Reactive Organic Matrices. Waste Manag. Res. 2009, 27, 101–111. [Google Scholar] [CrossRef] [PubMed]
  33. Alvarenga, P.; Gonçalves, A.P.; Fernandes, R.M.; de Varennes, A.; Vallini, G.; Duarte, E.; Cunha-Queda, A.C. Organic residues as immobilizing agents in aided phytostabilization: (I) Effects on soil chemical characteristics. Chemosphere 2009, 74, 1292–1300. Available online: http://www.sciencedirect.com/science/article/pii/S0045653508014793 (accessed on 21 March 2021). [CrossRef] [PubMed]
  34. Alvarenga, P.; Palma, P.; Gonçalves, A.P.; Fernandes, R.M.; de Varennes, A.; Vallini, G.; Duarte, E.; Cunha-Queda, A.C. Organic Residues as Immobilizing Agents in Aided Phytostabilization: (II) Effects on Soil Biochemical and Ecotoxicological Characteristics. Chemosphere 2009, 74, 1301–1308. [Google Scholar] [CrossRef]
  35. Solomon, M.; Tornos, F.; Large, R.R.; Badham, J.N.P.; Both, R.A.; Zaw, K. Zn–Pb–Cu Volcanic-Hosted Massive Sulphide Deposits: Criteria for Distinguishing Brine Pool-Type from Black Smoker-Type Sulphide Deposition. Ore Geol. Rev. 2004, 25, 259–283. [Google Scholar] [CrossRef]
  36. Durães, N.; Bobos, I.; da Silva, E.F. Speciation and Precipitation of Heavy Metals in High-Metal and High-Acid Mine Waters from the Iberian Pyrite Belt (Portugal). Environ. Sci. Pollut. Res. 2017, 24, 4562–4576. [Google Scholar] [CrossRef]
  37. Luís, A.; Grande, J.; Davila, J.M.; Aroba, J.; Durães, N.; Almeida, S.; de la Torre, M.; Sarmiento, A.M.; Fortes, J.C.; Ferreira da Silva, E.F.; et al. Application of Fuzzy Logic Tools for the Biogeochemical Characterisation of (Un)Contaminated Waters from Aljustrel Mining Area (South Portugal). Chemosphere 2018, 211, 736–744. [Google Scholar] [CrossRef]
  38. Luís, A.T.; Durães, N.; de Almeida, S.F.P.; da Silva, E.F. Integrating Geochemical (Surface Waters, Stream Sediments) and Biological (Diatoms) Approaches to Assess AMD Environmental Impact in a Pyritic Mining Area: Aljustrel (Alentejo, Portugal). J. Environ. Sci. 2016, 42, 215–226. [Google Scholar] [CrossRef]
  39. Luís, A.T.; Teixeira, P.; Almeida, S.F.P.; Ector, L.; Matos, J.X.; Ferreira da Silva, E.A. Impact of Acid Mine Drainage (AMD) on Water Quality, Stream Sediments and Periphytic Diatom Communities in the Surrounding Streams of Aljustrel Mining Area (Portugal). Water Air Soil Pollut. 2009, 200, 147–167. [Google Scholar] [CrossRef]
  40. Maia, F.; Pinto, C.; Waerenborgh, J.C.; Gonçalves, M.A.; Prazeres, C.; Carreira, O.; Sério, S. Metal Partitioning in Sediments and Mineralogical Controls on the Acid Mine Drainage in Ribeira Da Água Forte (Aljustrel, Iberian Pyrite Belt, Southern Portugal). Appl. Geochem. 2012, 27, 1063–1080. [Google Scholar] [CrossRef]
  41. Relvas, J.; Pinto, A.; Matos, J. Lousal, Portugal: A Successful Example of Rehabilitation of a Closed Mine in the Iberian Pyrite Belt. Soc. Geol. Appl. Miner. Depos. SGA News 2012, 31, 1–16. [Google Scholar]
  42. Andráš, P.; Matos, J.X.; Turisová, I.; Batista, M.J.; Kanianska, R.; Kharbish, S. The Interaction of Heavy Metals and Metalloids in the Soil–Plant System in the São Domingos Mining Area (Iberian Pyrite Belt, Portugal). Environ. Sci. Pollut. Res. 2018, 25, 20615–20630. [Google Scholar] [CrossRef]
  43. Pereira, R.; Sousa, J.P.; Ribeiro, R.; Gonçalves, F. Microbial Indicators in Mine Soils (S. Domingos Mine, Portugal). Soil Sediment Contam. Int. J. 2006, 15, 147–167. [Google Scholar] [CrossRef]
  44. Gee, G.W.; Bauder, J.W. Particle-size analysis. In Methods of Soil Analysis. Part 1. Physical and Mineralogical Methods; Klute, A., Ed.; Soil Science Society of America: Madison, WI, USA, 1986; pp. 383–412. [Google Scholar]
  45. USDA. Soil Science Division Staff. Soil Survey Manual. In USDA Handbook 18; Ditzler, C., Scheffe, K., Monger, H.C., Eds.; Government Printing Office: Washington, DC, USA, 2017. Available online: https://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/ref/?cid=nrcs142p2_054253#soil_texture (accessed on 3 September 2021).
  46. Walkley, A.; Black, J.A. An examination of the Degtjareff method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Sci. 1934, 37, 29–38. [Google Scholar] [CrossRef]
  47. Riehm, H. Die ammoniumlaktatessigsaure-methode zur bestimmung der leichtoeslichen phosphosaure in karbonathaltigen boden. Agrochimica 1958, 3, 49–65. [Google Scholar]
  48. ISO 11466; Soil Quality-Extraction of Trace Elements Soluble in Aqua Regia. International Organisation for Standardisation: Genève, Switzerland, 1995.
  49. LQARS. Manual de Fertilização de Culturas; Laboratório Químico Agrícola Rebelo da Silva, Instituto Nacional de Investigação Agrária e Veterinária–INIAV: Oeiras, Portugal, 2006. [Google Scholar]
  50. Alvarenga, P.; Gonçalves, A.P.; Fernandes, R.M.; de Varennes, A.; Duarte, E.; Vallini, G.; Cunha-Queda, A.C. Effect of organic residues and liming materials on metal extraction from a mining-contaminated soil. Bioremediat. J. 2008, 12, 58–69. Available online: http://www.tandfonline.com/doi/abs/10.1080/10889860802059909?journalCode=bbrm20 (accessed on 18 May 2018). [CrossRef]
  51. Alvarenga, P.; Rodrigues, D.; Mourinha, C.; Palma, P.; de Varennes, A.; Cruz, N.; Tarelho, L.A.C.; Rodrigues, S. Use of wastes from the pulp and paper industry for the remediation of soils degraded by mining activities: Chemical, biochemical and ecotoxicological effects. Sci. Total Environ. 2019, 686, 1152–1163. [Google Scholar] [CrossRef]
  52. Santos, J.Q. Fertilização, Fundamentos da Utilização dos Adubos e Corretivo; Euroagro, C., Ed.; Publicações Europa-América: Lisboa, Portugal, 1996. [Google Scholar]
  53. Alvarenga, P.; Varennes, A.; Cunha-Queda, A.C. The Effect of Compost Treatments and A Plant Cover with Agrostis tenuis on the Immobilization/Mobilization of Trace Elements in a Mine-Contaminated Soil. Int. J. Phytoremediat. 2014, 16, 138–154. [Google Scholar] [CrossRef]
  54. Touceda-González, M.; Álvarez-López, V.; Prieto-Fernández, Á.; Rodríguez-Garrido, B.; Trasar-Cepeda, C.; Mench, M.; Puschenreiter, M.; Quintela-Sabarís, C.; Macías-García, F.; Kidd, P.S. Aided Phytostabilisation Reduces Metal Toxicity, Improves Soil Fertility and Enhances Microbial Activity in Cu-Rich Mine Tailings. J. Environ. Manag. 2017, 186, 301–313. [Google Scholar] [CrossRef] [PubMed]
  55. Kabata-Pendias, A. Trace Elements in Soils and in Plants, 4th ed.; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
Figure 1. Map evidencing the upper limit of the Iberian Pyrite Belt, with the most important mining districts that were responsible for extensive environmental impacts: Lousal (and Caveira), Aljustrel, and São Domingos.
Figure 1. Map evidencing the upper limit of the Iberian Pyrite Belt, with the most important mining districts that were responsible for extensive environmental impacts: Lousal (and Caveira), Aljustrel, and São Domingos.
Environments 09 00105 g001
Figure 2. Google Earth image (24 May 2022) evidencing the sampling sites at the Aljustrel mine area, covering an area already rehabilitated by EDM [5]. A1 and A4: sites near the confined waste materials, Algares area; A2 and A3: sites near the perimetral channels that collect drainage waters; A5–A8: area near recent tailing deposits; A9–A10: soils with scarce vegetation, but eventually representative of natural soils; and A11–A14: soils near old tailing deposits.
Figure 2. Google Earth image (24 May 2022) evidencing the sampling sites at the Aljustrel mine area, covering an area already rehabilitated by EDM [5]. A1 and A4: sites near the confined waste materials, Algares area; A2 and A3: sites near the perimetral channels that collect drainage waters; A5–A8: area near recent tailing deposits; A9–A10: soils with scarce vegetation, but eventually representative of natural soils; and A11–A14: soils near old tailing deposits.
Environments 09 00105 g002
Figure 3. Google Earth image (24 May 2022) evidencing the sampling sites at the Lousal mine area. L1–L5: area potentially affected by acid mine drainage from old tailing deposits; L6–L10; area potentially affected from the pyrite road and railway transport network; and L11–L12: sampling points near the exit of an old mine gallery, with tailing deposits.
Figure 3. Google Earth image (24 May 2022) evidencing the sampling sites at the Lousal mine area. L1–L5: area potentially affected by acid mine drainage from old tailing deposits; L6–L10; area potentially affected from the pyrite road and railway transport network; and L11–L12: sampling points near the exit of an old mine gallery, with tailing deposits.
Environments 09 00105 g003
Figure 4. Google Earth image (24 May 2022) evidencing the sampling sites at the São Domingos mine area, before the beginning of the EDM rehabilitation project. SD1–SD22: area affected by the original mining works and by the acid mine drainage with origin in the different wastes which were deposited along the São Domingos stream valley; SD23–SD30: soils potentially affected by acid mine drainage, located downstream from Achada do Gama sulfur factory; and SD31–SD44: soils with diverse characteristics, scarce vegetation, but eventually representative of soils less affected by the mine works.
Figure 4. Google Earth image (24 May 2022) evidencing the sampling sites at the São Domingos mine area, before the beginning of the EDM rehabilitation project. SD1–SD22: area affected by the original mining works and by the acid mine drainage with origin in the different wastes which were deposited along the São Domingos stream valley; SD23–SD30: soils potentially affected by acid mine drainage, located downstream from Achada do Gama sulfur factory; and SD31–SD44: soils with diverse characteristics, scarce vegetation, but eventually representative of soils less affected by the mine works.
Environments 09 00105 g004
Figure 5. Boxplot representation of the most representative physicochemical properties to evaluate the fertility status of the soils in the different mining areas (sampling campaign July 2016): Aljustrel (n = 14), Lousal (n = 12), São Domingos (n = 44). (a) pH (H2O, 1:5 v/v); (b) Electrical conductivity (EC, mS/cm, 1:5 v/v); (c) Organic matter content (OM, % w/w); (d) Kjeldahl N (%, w/w); (e) extractable P (mg kg-1); and (f) extractable K (mg kg-1). The mark across the box represents the median and the bottom, and the top of the box represents the 25–75% interval. All results report to a dry matter basis.
Figure 5. Boxplot representation of the most representative physicochemical properties to evaluate the fertility status of the soils in the different mining areas (sampling campaign July 2016): Aljustrel (n = 14), Lousal (n = 12), São Domingos (n = 44). (a) pH (H2O, 1:5 v/v); (b) Electrical conductivity (EC, mS/cm, 1:5 v/v); (c) Organic matter content (OM, % w/w); (d) Kjeldahl N (%, w/w); (e) extractable P (mg kg-1); and (f) extractable K (mg kg-1). The mark across the box represents the median and the bottom, and the top of the box represents the 25–75% interval. All results report to a dry matter basis.
Environments 09 00105 g005
Figure 6. Boxplot representation of the pseudo-total concentrations (Aqua Regia extraction) of the most representative potentially toxic elements which affect pyrite mines in the Portuguese sector of the IPB (sampling campaign July 2016): Aljustrel (n = 14), Lousal (n = 12), São Domingos (n = 39). (a) As (mg/kg); (b) Cu (mg/kg); (c) Pb (mg/kg); and (d) Zn (mg/kg), all concentrations in a dry weight (DW) basis. The mark across the box represents the median and the bottom, and the top of the box represents the 25–75% interval. All results report to a dry matter basis. (*) Recommended limit values for the remediation of soil intended for industrial use from the APA Technical Guide are presented in red and also as a horizontal line [2].
Figure 6. Boxplot representation of the pseudo-total concentrations (Aqua Regia extraction) of the most representative potentially toxic elements which affect pyrite mines in the Portuguese sector of the IPB (sampling campaign July 2016): Aljustrel (n = 14), Lousal (n = 12), São Domingos (n = 39). (a) As (mg/kg); (b) Cu (mg/kg); (c) Pb (mg/kg); and (d) Zn (mg/kg), all concentrations in a dry weight (DW) basis. The mark across the box represents the median and the bottom, and the top of the box represents the 25–75% interval. All results report to a dry matter basis. (*) Recommended limit values for the remediation of soil intended for industrial use from the APA Technical Guide are presented in red and also as a horizontal line [2].
Environments 09 00105 g006
Figure 7. Scores of each sample on the three main principal components: (a): PC1 versus PC2; (b) PC1 versus PC3. The most important parameters for the definition of the two components are shown on the edge of each axis, indicating the direction in which the value of the parameter increases (pH: soil pH(H2O); and Cutotal, Pbtotal and Zntotal: pseudo-total elements concentration aqua-regia digestion). The color of the label indicates the mine: red for Aljustrel; blue for Lousal; and green for São Domingos.
Figure 7. Scores of each sample on the three main principal components: (a): PC1 versus PC2; (b) PC1 versus PC3. The most important parameters for the definition of the two components are shown on the edge of each axis, indicating the direction in which the value of the parameter increases (pH: soil pH(H2O); and Cutotal, Pbtotal and Zntotal: pseudo-total elements concentration aqua-regia digestion). The color of the label indicates the mine: red for Aljustrel; blue for Lousal; and green for São Domingos.
Environments 09 00105 g007
Figure 8. Scores of each sample on the principal components, PC2 versus PC3. The most important parameters for the definition of the two components are shown on the edge of each axis, indicating the direction in which the value of the parameter increases (Cutotal, Pbtotal, and Zntotal: pseudo-total element concentrations aqua-regia digestion). The color of the label indicates the mine: red for Aljustrel; blue for Lousal; and green for São Domingos.
Figure 8. Scores of each sample on the principal components, PC2 versus PC3. The most important parameters for the definition of the two components are shown on the edge of each axis, indicating the direction in which the value of the parameter increases (Cutotal, Pbtotal, and Zntotal: pseudo-total element concentrations aqua-regia digestion). The color of the label indicates the mine: red for Aljustrel; blue for Lousal; and green for São Domingos.
Environments 09 00105 g008
Table 1. Pearson’s correlation coefficients between soil physicochemical properties and potentially toxic elements pseudo-total concentrations (n = 65). Marked correlations are significant at (*) p < 0.05; (**) p < 0.01; and (***) p < 0.001, evidencing the strength of the associations: moderate, high, and very high, respectively.
Table 1. Pearson’s correlation coefficients between soil physicochemical properties and potentially toxic elements pseudo-total concentrations (n = 65). Marked correlations are significant at (*) p < 0.05; (**) p < 0.01; and (***) p < 0.001, evidencing the strength of the associations: moderate, high, and very high, respectively.
pHECOMNKjeldalhPextractableKextractableAstotalCutotalPbtotal
EC−0.54 ***--------
OM0.33 **−0.15-------
NKjeldalh−0.03−0.05−0.05------
Pextractable0.31 *−0.030.68 ***0.02-----
Kextractable0.28 *−0.220.51 ***0.060.09----
Astotal−0.36 **0.19−0.050.39 **0.04−0.20---
Cutotal0.06−0.070.11−0.06−0.080.27 *−0.04--
Pbtotal−0.230.200.12−0.040.13−0.160.42 ***0.30 *-
Zntotal0.24−0.03−0.01−0.040.05−0.13−0.120.43 ***0.15
Table 2. Factor loadings of each soil parameter along the first four principal components, PC1, PC2, PC3, and PC4, which resulted from the PCA. Marked correlations, in red, are significant (correlation coefficient > 0.7).
Table 2. Factor loadings of each soil parameter along the first four principal components, PC1, PC2, PC3, and PC4, which resulted from the PCA. Marked correlations, in red, are significant (correlation coefficient > 0.7).
Soil ParameterPC1PC2PC3PC4
pH0.7900.1550.065−0.150
EC−0.573−0.261−0.0150.442
OM0.697−0.501−0.3030.218
NKjeldalh−0.127−0.193−0.364−0.781
Pextractable0.501−0.515−0.3730.311
Kextractable0.620−0.059−0.129−0.188
Astotal−0.486−0.587−0.292−0.351
Cutotal0.220−0.3780.702−0.224
Pbtotal−0.240−0.7640.2400.088
Zntotal0.182−0.2150.734−0.105
Explained variance (%)24.617.615.612.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alvarenga, P.; Mourinha, C.; Palma, P.; Cruz, N.; Rodrigues, S.M. Assessment of Soil Physicochemical Characteristics and As, Cu, Pb and Zn Contamination in Non-Active Mines at the Portuguese Sector of the Iberian Pyrite Belt. Environments 2022, 9, 105. https://doi.org/10.3390/environments9080105

AMA Style

Alvarenga P, Mourinha C, Palma P, Cruz N, Rodrigues SM. Assessment of Soil Physicochemical Characteristics and As, Cu, Pb and Zn Contamination in Non-Active Mines at the Portuguese Sector of the Iberian Pyrite Belt. Environments. 2022; 9(8):105. https://doi.org/10.3390/environments9080105

Chicago/Turabian Style

Alvarenga, Paula, Clarisse Mourinha, Patrícia Palma, Nuno Cruz, and Sónia Morais Rodrigues. 2022. "Assessment of Soil Physicochemical Characteristics and As, Cu, Pb and Zn Contamination in Non-Active Mines at the Portuguese Sector of the Iberian Pyrite Belt" Environments 9, no. 8: 105. https://doi.org/10.3390/environments9080105

APA Style

Alvarenga, P., Mourinha, C., Palma, P., Cruz, N., & Rodrigues, S. M. (2022). Assessment of Soil Physicochemical Characteristics and As, Cu, Pb and Zn Contamination in Non-Active Mines at the Portuguese Sector of the Iberian Pyrite Belt. Environments, 9(8), 105. https://doi.org/10.3390/environments9080105

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