Geochemical Characterization and Trace-Element Mobility Assessment for Metallic Mine Reclamation in Soils Affected by Mine Activities in the Iberian Pyrite Belt

: The geochemical characterization of the mine deposits and soils in metal mining areas is essential in order to develop an effective mine reclamation strategy. The determination of total potentially toxic element (PTE) content, together with the application of chemical extraction procedures, can give insight into the behavior of contaminants after the application of different mine reclamation solutions, as well as identify the areas where urgent action is needed. This work presents a practical application to the evaluation of the pollution potential of trace elements in soils affected by mining activities, to be used in metallic mine reclamation. The PTE behavior was assessed by single extractions in order to simulate four environmental conditions: PTE mobility under rainfall conditions, acid mine drainage, reducing conditions, and plant uptake. The spatial distribution of contaminants in the study area was evaluated by determination of PTE total content in soil samples. Trace elements with high natural mobility, such as Zn, appeared concentrated at water and sediment discharge areas, while As, Pb, and Cu contents were higher near the mine wastes. The results obtained after the extractions suggested that the highest PTE content was extracted in the complexing–reducing medium, due to the dissolution of secondary sulfates and Fe 3+ oxyhydroxides and the subsequent release of PTEs associated with those mineral phases. Reclamation strategies applied in the study area should promote efﬁcient water drainage, inﬁltration, and subsuperﬁcial water circulation in order to maintain oxidant conditions in the soil. The methodology applied in this study may constitute a valuable tool to deﬁne the geochemical constraints in metal mining areas, as well as help to develop appropriate mine reclamation solutions.


Introduction
The long story of metalliferous mining, processing, and smelting activities in the Iberian Pyrite Belt (southwest of the Iberian Peninsula) has left behind a legacy of numerous abandoned mines and considerable quantities of waste rock piles, spoil heaps, and tailing ponds [1][2][3]. Modern day extraction and treatment processes are usually very efficient, and the resulting mining wastes contain relatively low concentrations of potentially toxic elements (PTEs). However, wastes derived from historical mining operations are characterized by higher PTE concentrations and a relative abundance of ore minerals (i.e., pyrite, galena, and sphalerite), due to the use of less efficient mineral processing and extraction technologies [4,5]. When oxidation of the sulfides remaining in those mine wastes occurs, very coarse texture. In some areas, poorly developed soils over the tertiary deposits are 150 also present, characterized by a sandy texture, very limited organic matter content, and a 151 slightly acidic nature. Due to the mining operations, most of the soils present at the mine 152 site can be considered as Technosols [41], composed by mixtures of soil, rock fragments, 153 and tailing material [42,43]. Belt. This complex comprises felsic volcanic rocks, interbedded with black shales, volcano-161 clastic sandstones, and detrital sediments [30,44,45]. In the Lousal region, the volcano-162 sedimentary complex is overlain by the Tertiary-aged sediments of the Alto Sado Basin 163 [42]. The mineralization hosted is related to an exhalative deposit, presumably developed 164 in brine pools within the lower units of the volcano-sedimentary complex [44]. The ore 165 deposit is characterized by fine-grained massive sulfides with a mineralization dominated 166 by pyrite, with variable amounts of chalcopyrite, galena, sphalerite, pyrrhotine (Fe(1−x)S), 167 marcasite (FeS2), bournonite (PbCuSbS3), tetrahedrite (Cu12Sb4S13), arsenopyrite, cobaltite 168 (CoAsS), magnetite (Fe 2+ Fe 3+ 2O4), and native gold [46]. The Lousal area lies in the SW limb of an anticline affected by N-S and NE-SE late-stage faults, along a lineament of the volcano-sedimentary complex of the Iberian Pyrite Belt. This complex comprises felsic volcanic rocks, interbedded with black shales, volcanoclastic sandstones, and detrital sediments [30,44,45]. In the Lousal region, the volcano-sedimentary complex is overlain by the Tertiary-aged sediments of the Alto Sado Basin [42]. The mineralization hosted is related to an exhalative deposit, presumably developed in brine pools within the lower units of the volcano-sedimentary complex [44]. The ore deposit is characterized by fine-grained massive sulfides with a mineralization dominated by pyrite, with variable amounts of chalcopyrite, galena, sphalerite, pyrrhotine (Fe (1−x) S), marcasite (FeS 2 ), bournonite (PbCuSbS 3 ), tetrahedrite (Cu 12 Sb 4 S 13 ), arsenopyrite, cobaltite (CoAsS), magnetite (Fe 2+ Fe 3+ 2 O 4 ), and native gold [46]. The Lousal mine, located in the Grândola municipality (Alentejo Region), was exploited between 1900 and 1988, primarily for pyrite, through surface and underground works that went to depths of about 500 m [36,40,47]. During the first two decades of ex-ploitation, various gossan deposits and supergene enrichment zones were worked mainly for copper ore. Later, in the 1930s, the mining activity increased significantly and shifted toward the exploitation of deeper ore lenses, due to the mining company needs of production of high tonnages of pyrite concentrate [47]. This trend continued for the remaining of the mine life cycle, with a notable increase in the pyrite concentrates during the late 1950s and early 1960s.
The mining legacy has resulted in a large volume of wastes, ranging from barren materials to various types of waste deposits and tailing impoundments, estimated to be greater than 1 Mt [43,44]. Two main tailings and waste material impoundments can currently be identified in the Lousal mine area, one related to the mine shafts and ore mill plants, and another related to a big deposit of fine-grained pyrite close to the railway station at the NE of the mining facilities, where the ore was deposited and later transported for chemical treatment [42,44]. During the mine life cycle, most of the waste deposits were left exposed to weathering, and no effective environmental measures were taken to reduce tailings erosion and avoid the generation of AMD. Due to the fast environmental degradation of the mine wastes after the mine closure in 1988, a rehabilitation program (RELOUSAL, Rehabilitation and Integral Development of the Lousal Mine) was promoted by the mining company (SAEPEC) and the Grândola town council. As part of this program, an AMD passive treatment facility was implemented in the northeastern area of the mining complex, next to the pyrite deposits. A system of artificial wetlands and open limestone channels was constructed in order to reduce PTE content and increase the pH of surficial waters that drained into the Corona stream [48]. Additionally, a geotextile fabric was extended over a section of the pyrite deposits next to the railway station and then covered with a soil layer, in order to encapsulate the wastes and reduce the impact of weathering processes.
Currently, a European-funded mine reclamation project (LIFE RIBERMINE, LIFE ENV/ES/000181) is starting to be applied in the soils affected by the pyrite-rich mine wastes. This project will be the first worldwide to apply a combined fluvial-geomorphic and geochemical approach to the reclamation of the mining area. The topography of the affected mine soils will be redesigned and reconstructed using the GeoFluv method [49][50][51] in order to reduce accelerated water-driven erosion, and a buffering soil amendment will be extended over the reconstructed topography to prevent the chemical weathering of the mine soils and the formation of acid mine drainage. Reclamation actions are being implemented during the years 2021 and 2022.

Sampling Design
The sampling campaign was conducted in January 2019. Surficial and deep soil samples were collected from the mine soils. Deep samples were collected following two transects (Transect W and Transect E), parallel to the direction of maximum slope ( Figure 2). In each transect, several trenches of 1 m depth were excavated with a backhoe, and soil samples were collected from the trench walls at different depths. Samples collected down to a depth of 20 cm were considered as surficial samples, while samples collected below the first 20 cm of the soil profile were considered as deep samples. In total, 28 samples were collected using this procedure, 19 in Transect W and nine in Transect E. Alongside the samples collected in the two transects, 16 surficial samples were collected following the flow trajectory of two dry ephemeral channels developed in the soils affected by mining activities, and one sample (LOS17) was collected from sediments accumulated in a catch basin, a part of a drainage channel built to transport ADM leakages to the treatment system. The location of these samples corresponds to the areas where the topographic reconstruction of the mine wastes is expected to be implemented. Sampling was carried out by digging down below the first 10 cm of the soil and taking the sample directly with a hand shovel. All the samples were air-dried, passed through a 2 mm sieve, homogenized, and stored in plastic bags at room temperature until laboratory analysis. a hand shovel. All the samples were air-dried, passed through a 2 mm sieve, homoge-223 nized, and stored in plastic bags at room temperature until laboratory analysis.

227
Elemental concentration data were achieved by means of energy-dispersive X-ray 228 fluorescence spectrometry (XRF), using a Malvern PANalytical, mod. Epsilon1 spectrom-229 eter at the Instituto de Geología Aplicada, IGEA (Castilla-La Mancha University). All the 230 collected soil samples were ground to fine powder (approximate grain size of 100 µ m) 231 using an agate ball mill prior to analysis. Sieved samples were placed into cylindrical sam-232 ple holders with a 3.6 µ m Mylar filter and then introduced in the XRF spectrometer. An 233 extended analysis time (23 min) was used to achieve the minimum standard deviation of 234 each data. The fluorescence spectra of each sample were corrected and adapted before 235 quantification to attenuate the effects of interelement interferences. From the total list of 236 major, minor, and trace elements analyzed, As, Cd, Cu, Fe (Fe2O3), Pb, S (SO3), and Zn 237 were specially chosen for this study because of their abundance and variability in these 238 types of mine wastes, and because most of them are included in the priority contaminant 239 list of environmental protection agencies [52]. To ensure the quality of the analyses, du-240 plicate samples were analyzed to check precision, and a certified reference material (SRM 241 2711) was analyzed for accuracy. The precision results of the studied elements were as 242 follows: Zn (87.3%), As (82.9%), Pb (93.1%), SO3 (93.5%), Fe2O3 (88.4%), and Cu (85.0%). 243 The recovery percentages for the certified elements in SRM 2711 were as follows (in %): S 244 (178.7-182.

Geochemical Analysis
Elemental concentration data were achieved by means of energy-dispersive X-ray fluorescence spectrometry (XRF), using a Malvern PANalytical, mod. Epsilon1 spectrometer at the Instituto de Geología Aplicada, IGEA (Castilla-La Mancha University). All the collected soil samples were ground to fine powder (approximate grain size of 100 µm) using an agate ball mill prior to analysis. Sieved samples were placed into cylindrical sample holders with a 3.6 µm Mylar filter and then introduced in the XRF spectrometer. An extended analysis time (23 min) was used to achieve the minimum standard deviation of each data. The fluorescence spectra of each sample were corrected and adapted before quantification to attenuate the effects of interelement interferences. From the total list of major, minor, and trace elements analyzed, As, Cd, Cu, Fe (Fe 2 O 3 ), Pb, S (SO 3 ), and Zn were specially chosen for this study because of their abundance and variability in these types of mine wastes, and because most of them are included in the priority contaminant list of environmental protection agencies [52]. To ensure the quality of the analyses, duplicate samples were analyzed to check precision, and a certified reference material (SRM 2711) was analyzed for accuracy. The precision results of the studied elements were as follows: Zn (87.3%), As (82.9%), Pb (93.1%), SO 3

Chemical Simple Extractions
In order to assess the mobility of PTEs (Pb, Zn, Cd, Cu, and As) present in the soils affected by mining activities in the Lousal mine, five selective extraction procedures were applied to the surficial samples: water extraction, acidic medium extraction, complexingreducing medium extraction, EDTA extraction, and DTPA extraction. These selective extractions allowed for the simulation of the potential mobility of PTEs under different environmental situations. Due to the strong association and affinity of heavy metals with fine-grained soil components, all samples were ground to fine powder (grain size < 63 µm) using an agate ball mill prior to the chemical extractions.
Water extraction was used to assess the natural mobility of PTEs soluble under actual conditions, which represent the most mobile fraction of trace elements in the sediment samples. Dilute aqueous salt solutions (CaCl 2 ) are often used to determine the watersoluble PTE fraction in soils [53,54]. In this study, however, we used distilled water, as its composition is more similar to that of natural rainfall water. PTE content was determined in the extracts of a 1:2 (w/v) soil/water suspension [55]. The soil/water mixture was then agitated for 24 h using a Multi-Purpose Tube Rotator (Fisherbrand) and centrifuged at 3000 rpm for 5 min in order to separate the supernatant solution from the solid residue. The solution was then filtered with 0.45 µm pore spacing and kept refrigerated prior to analysis.
Acidic medium extraction was applied to simulate the effect of AMD in the selected materials and assess the potential mobilization of PTEs under acidic conditions. PTE content was determined from the extracts of a 1:50 (w/v) soil/acid solution suspension. For this extraction, 1 g of sample was mixed with 50 mL of 0.1 mol −1 HNO 3 solution [56]. The suspension was agitated for 1 h and then centrifuged at 3000 rpm for 20 min. The liquid fraction was finally filtered with 0.45 µm pore spacing and kept refrigerated before laboratory analysis.
A complexing-reducing medium extraction was used to assess the PTE mobilization under anoxic conditions, using the dithionite-citrate extraction method [57]. The dithionitecitrate extraction partially mobilizes goethite, hematite, other non-crystalline Fe-oxides, and Fe-Al organic complexes present in the soil sample, thus allowing for the assessment of the PTE fraction adsorbed to those non-crystalline Fe compounds [57]. For each sample, 1 g was weighed and mixed with 40 mL of sodium citrate (chelating agent) and 5 mL of a 1 N NaHCO 3 solution (buffering solution). The obtained soil/reagent suspension was heated up to 80 • C and 1 g of sodium dithionite was added. While heated, it was continually mixed for 1 min and agitated intermittently for another 15 min. The final suspension was centrifuged at 3000 rpm for 10 min, and the liquid fraction was filtered with 0.45 µm pore spacing and kept refrigerated until laboratory analysis.
EDTA extraction was used to assess the PTE fraction assimilable by plants. EDTA is a widely used complexing agent that is assumed to extract metals on exchange sites of both organic and inorganic complexes, and the PTE fraction mobilized by this selective extraction has been proven to show a good correlation with plant uptake [58,59]. PTE content was determined in the extracts of a 1:5 (w/v) soil/EDTA suspension. For each collected sample, 5 g was weighted and placed in a volumetric flask, and 25 mL of 0.05 M Na 2 -EDTA solution was added. The soil/EDTA suspension was agitated for 2 h and then centrifuged at 2000 rpm for 20 min. The supernatant solution was separated from the solid residue, filtered with 0.45 µm pore spacing, and kept refrigerated prior to analysis. DTPA extraction has a similar application as the EDTA extraction, as it is commonly used to estimate plant-available PTE content in soils. EDTA extraction has been used successfully under many different soil conditions, while DTPA is often used as an extractant for micronutrients in moderate to strongly acidic soils [60]. PTE content was determined following the procedure developed by Lindsay and Norvell [61]. For each sample, 5 g was weighted and mixed with 10 mL of extractant solution (0.05 N DTPA, 0.01 N CaCl 2 , and 0.1 M triethanolamine). The soil/solution suspension was agitated for 2 h and then centrifuged at 3000 rpm for 20 min. The liquid fraction was filtered with 0.45 µm pore spacing and kept refrigerated before laboratory analysis.
The Pb, Zn, Cd, and Cu content in the solutions obtained by the application of the five selective extraction procedures was measured by atomic absorption spectrometry (AAS) using a ContrAA-800 high-intensity continuous spectrum atomic absorption spectrometer (Analytik Jena) at the Castilla-La Mancha University. To ensure the quality of the analyses, duplicate samples were analyzed to check precision. Each sample was analyzed thrice to verify the obtained results. In order to avoid wrong measurements due to the effect of matrix interference, a reagent blank was prepared for each of the extraction procedures, and calibration standards were prepared using the same reagent matrix used in each extraction. A calibration curve was made for each of the analyzed trace elements, and a certified standard was analyzed for accuracy. The coefficient of determination (R 2 ) for the studied elements was as follows: Pb (0.998-0.999), Zn (0.970-0.999), Cd (0.987-0.999), Cu (0.997-0.999). The recovery percentages for the analyzed elements in the certified standard were as follows (in %): Pb (95-101), Zn (97-98), Cd (98-102), Cu (94-98). The As content was determined by atomic fluorescence spectrometry using a coupled HPCL (highperformance liquid chromatography) short column and hydride generation spectrometer (PSA Millenium Excalibur 10.055) at the Castilla-La Mancha University. A fraction of each sample was submitted to a chemical reduction prior to the laboratory analysis, which reduces As(V) to As(III) and facilitates As detection when using the hydride generation spectrometer. In order to do that, a solution of Kl/ascorbic acid was applied to the soil samples and left at room temperature for 30 min. The quality control of the analyses was performed following the same procedures as those previously described for the analysis of Pb, Zn, Cd, and Cu. The coefficient of determination (R 2 ) for the calibration curves of As ranged between 0.998 and 0.999. The recovery percentages of As for each extraction were as follows (in %): water extraction (104), acidic medium (104), complexing-reducing medium (107), EDTA (108), and DTPA (107).

Mineralogical Analysis
The mineralogy was defined through X-ray diffraction (XRD). Original soil samples and residues after each extraction were air-dried and ground in an agate mortar to fine powder (grain size < 63 µm). Mineral phases were identified by powder XRD using a Bruker D8 Advance diffractometer equipped with a Cu anticathode. Diffractograms were recorded continuously at 2θ angles from 2 • to 68 • , with 0.02 stepping intervals and 1 s per step. The semiquantitative analysis was carried out according to the Chung method [62][63][64], using XPowder software [65]. The Chung method is based on the determination of the reference intensity ratios (RIR) of the existing phases, which allows the intensity calculations to be normalized according to the assumption that the sum of all phases in the sample is equal to 100%. The XRD analyses were performed at the Complutense University of Madrid (Laboratory of Geological Techniques, CAI, Complutense University of Madrid, Madrid, Spain).

Multivariate Analysis
A multivariate analysis was applied for the geochemical and mineralogical variables, including a correlation matrix and principal component and factor analysis, using Statgraphics 19 software [66]. The multivariate analysis was used to determine potential relationships between the PTE total content and the mineralogical composition of the original soil samples. The correlation matrix was applied to identify the relationship between trace element and mineralogy in the original soil samples, using Pearson's product moment correlation coefficient [67]. Factorial analysis permits a statistical approximation for analyzing interrelations between a large number of variables. The factorial analysis was carried out by the principal component extraction method, using a varimax normalized rotation of the factors [68].

Assessment of Potential Environmental Risk
A significant number of indicators designed to approximate the quality of soils can readily be found in the literature [69][70][71]. In this study, assessment of soil contamination level was performed by the quantification of the contamination factor (CF) and the pollution load index (PLI). These indicators are useful tools to assess the degree of trace-element contamination in soils, and they allow for the comparison of results with other potentially polluted areas.

Contaminant Factor
The contaminant factor (CF) was calculated as the ratio between the PTE concentrations in soil samples using the Portuguese Background values [72].
CF < 1 indicates low contamination, 1 < CF < 3 indicates moderate contamination, 3 < CF < 6 indicates considerable contamination, and CF > 6 indicates very high contamination [73]. CF values were calculated for PTEs with available geochemical background values: Pb, Zn, Cd, As, and Cu. In the case of Cd, only some of the surficial soil samples had an associated CF value, due to most of the samples showing Cd concentrations below the analytical detection limit.

Pollution Load Index
For the entire soil sampling site, the pollution load index (PLI) was determined as the n-th root of the product of the n CF. In samples with quantifiable Pb, Zn, Cd, As, and Cu, the n value used was 6. In samples with one or more PTEs below the analytical detection limit, the n value was reduced by the amount of PTE not quantified. In general, PLI < 1 indicates low contamination, 1 < PLI < 3 indicates moderate contamination, 3 < PLI < 6 indicates considerable contamination, and PLI > 6 indicates very high contamination.

Results and Discussion
The results of the geochemical and mineralogical characterization of the soil samples affected by mining activities, as well as those obtained from the assessment of PTE mobility under different environmental conditions, are presented here.

Geochemical Characterization
Total SiO 2 , SO 3 , Fe 2 O 3 , and trace-element concentrations obtained from surficial and deep samples are presented in Tables 1 and 2. The general composition of the two transects and surficial samples was similar, characterized by notably high SiO 2 (41.25 to 62.58 wt.%) and Fe 2 O 3 content (8.73 to 57.86 wt.%), as well as moderate to high SO 3 content (0.06 to 18.59 wt.%). Table 1. Geochemical characterization of surficial soil samples.  Table 2. Geochemical characterization of soil samples collected from the two transects (surficial and deep soil samples). Trace elements such as As, Cu, Zn, and Pb were not detected in very high concentrations, although their variations at different sampling depths and locations showed an anomalous distribution. The As content ranged between 18 and 1230 mg·kg −1 , with maximum values in samples LOS17 and LO2.1-01 (Tables 1 and 2), and Pb ranged from 25 to 5730 mg±kg −1 (maximum value in sample LO2.1-02, at the beginning of transect E). Other PTEs that showed noteworthy concentration values were Cu (up to 1320 mg·kg −1 , sample LOS6) and Zn (up to 546 mg·kg −1 , in sample LOS6). In the case of Cd, the soil samples presented a very low content, ranging from a minimum value of 0.2 mg·kg −1 to a maximum value of only 2.6 mg·kg −1 in sample LOS7. Out of the 45 soil samples collected, 31 showed Cd content below the analytical detection limit.

Fe 2 O 3 (%) SiO 2 (%) SO 3 (%) Pb
The distribution of total concentration for the elements of interest is discussed next. For Fe 2 O 3 , the general tendency showed higher values near the surface and a progressive decrease with increasing sampling depth ( Figure 3). In transect W, it is possible to observe that samples collected from a higher topographic position had lower concentrations than those in the low section of the transect. In this transect, Fe 2 O 3 content seemed to be fairly consistent within the deep samples, while a progressive increase in concentration was found for the superficial samples from the beginning of the transect (at higher topographic location) to the end of the transect. Taking into account that transect W follows the trajectory of an ephemeral channel, this increase in concentration could surely be related to the mobilization of surface iron oxides by water runoff. The highest value of Fe 2 O 3 (58 wt.%) was registered in sample LOS17, collected from the sediments accumulated in a catch basin.  Similarly to the aforementioned trace elements, the Cu contents also showed a general increase toward the lower section of transect W, with generally higher values in the surficial samples ( Figure 3). This tendency appeared disrupted in the beginning of the transect and the western section of the studied area, with anomalously high values of 1320 mg·kg −1 and 719 mg·kg −1 at surface level and at a depth of 50 cm, respectively. On the other hand, transect E showed a contrasting distribution of Cu content, with higher values related to the samples located at the beginning of the transect, near the encapsulated mine wastes. The lowest Cu values appeared at the end of transect E, in samples collected in a ridge between two small depressions or swales produced by water runoff.
The Zn concentration values were fairly consistent throughout the sampling area, with no anomalously high or low values ( Figure 3). Maximum Zn values (546 mg·kg −1 ) appeared at upper positions in the slope, near the encapsulated materials in the western section of the study area. Aside from those values, the general tendency showed a slight increase in Zn content toward the lower end of the slope, at lower topographic positions, for both surficial and deep samples. Similarly to Fe2O3 and SO3, this may indicate that Zn content was mobilized by water runoff from the materials at higher topographic positions, with relatively high Zn content, to the lower sections of the ephemeral channels. Samples collected from ridges on the slope (LO1.1 and LO1.2 in transect W, and LO2.3 in transect E) exhibited the lowest Zn values, confirming that Zn content in soils may be related to mobilization processes by water runoff.
The concentration of Pb was strongly anomalous in the lowest section of transect W and in the upper samples in transect E, showing strongly fluctuating values ( Figure 3). This element appeared to be concentrated near the surface, although the highest obtained value (5730 mg·kg −1 ) corresponded to a depth of approximately 20 cm in the samples located at the upper section of transect E. Although Pb content was generally low for most of the collected samples, it is worth noting that, similarly to the Zn content, Pb concentration appeared to be lower in samples collected at ridge positions in the slope (LO1.1 and LO1.2 in transect W, and LO2.3 in transect E), while generally higher values of Pb appeared in sediment samples collected from the ephemeral channels. In transect E, however, this tendency was reversed, with higher SO 3 values near the encapsulated materials at the beginning of the transect, most probably due to the proximity to the pollution source. The absolute highest SO 3 value was measured in the surficial soil sample LOS17 (18.59 wt.%).
The As concentrations showed a higher variability than that of Fe 2 O 3 , with high values specifically located at a depth of 50 cm in transect W and no clear distribution in transect E (Figure 3). In transect W, the As concentration progressively increased toward its lower section, near the AMD treatment pond, in both surficial and deep samples. The maximum values, however, corresponded to the samples located at the beginning of transect E, next to encapsulated pyrite wastes. These samples showed a significantly anomalous As concentration with a drastic reduction with depth, going from 1230 mg·kg −1 to 276 mg·kg −1 in less than 30 cm. Similarly to Fe 2 O 3 and SO 3 , the highest As content was found in the sample LOS17, probably related to As adsorption processes to Fe 3+ oxyhydroxides such as goethite.
Similarly to the aforementioned trace elements, the Cu contents also showed a general increase toward the lower section of transect W, with generally higher values in the surficial samples ( Figure 3). This tendency appeared disrupted in the beginning of the transect and the western section of the studied area, with anomalously high values of 1320 mg·kg −1 and 719 mg·kg −1 at surface level and at a depth of 50 cm, respectively. On the other hand, transect E showed a contrasting distribution of Cu content, with higher values related to the samples located at the beginning of the transect, near the encapsulated mine wastes. The lowest Cu values appeared at the end of transect E, in samples collected in a ridge between two small depressions or swales produced by water runoff.
The Zn concentration values were fairly consistent throughout the sampling area, with no anomalously high or low values ( Figure 3). Maximum Zn values (546 mg·kg −1 ) appeared at upper positions in the slope, near the encapsulated materials in the western section of the study area. Aside from those values, the general tendency showed a slight increase in Zn content toward the lower end of the slope, at lower topographic positions, for both surficial and deep samples. Similarly to Fe 2 O 3 and SO 3 , this may indicate that Zn content was mobilized by water runoff from the materials at higher topographic positions, with relatively high Zn content, to the lower sections of the ephemeral channels. Samples collected from ridges on the slope (LO1.1 and LO1.2 in transect W, and LO2.3 in transect E) exhibited the lowest Zn values, confirming that Zn content in soils may be related to mobilization processes by water runoff.
The concentration of Pb was strongly anomalous in the lowest section of transect W and in the upper samples in transect E, showing strongly fluctuating values ( Figure 3). This element appeared to be concentrated near the surface, although the highest obtained value (5730 mg·kg −1 ) corresponded to a depth of approximately 20 cm in the samples located at the upper section of transect E. Although Pb content was generally low for most of the collected samples, it is worth noting that, similarly to the Zn content, Pb concentration appeared to be lower in samples collected at ridge positions in the slope (LO1.1 and LO1.2 in transect W, and LO2.3 in transect E), while generally higher values of Pb appeared in sediment samples collected from the ephemeral channels.
The heterogeneity seen for the PTE content in surface and profound soil samples is common in small and highly anthropized mining areas. In those areas, the mine waste layout in heterogeneous mine dumps, usually mixing ore and gangue minerals, together with the PTE and particle dispersion to the environment through operational activities, leaching, and meteorological events, usually produces a distribution of PTEs with no apparent structure or clear tendency.
The surficial distribution of elements of interest is also shown in Figure 4. It includes information obtained from the surficial soil samples (LOS1 to LOS17) and samples collected from the first 20 cm of the soil profile in transect W and E.
Isocontent maps for Fe 2 O 3 exhibit two clearly distinct areas of surficial concentration of this PTE: (1) an area of sediment accumulation at the western section of the study area, corresponding to the lower section of one of the ephemeral channels, and (2) the upper section of transect E, next to the encapsulated mine waste deposit. In the first case, the distribution of high concentrations of PTEs could be related to the mobilization by runoff of sediments and waste materials from the upper section of the study area and consequent accumulation in the lower section of the channel. In the second case, the high content of trace elements in the surficial samples could be due to the proximity to the encapsulated mine wastes. The geotextile fabric and soil layer that covers this area reduce the amount of infiltrating water that makes contact with the mine wastes, decreasing the impact of weathering processes and the mobilization of PTEs to lower sections of the slope.
As and Pb showed a similar surficial distribution, with higher values at the upper section of transect W and transect E. These two trace elements had notably distinct mobilization conditions. Therefore, it is reasonable to think that the coincidence in their distribution pattern correlated to the proximity to the source of PTE contamination (in this case, the encapsulated mine wastes at the top of the deposit) and not to the mobilization processes. Cu concentration in soil samples followed a similar pattern, with higher values on upper topographic positions. In this case, however, maximum concentration values corresponded to transect W and the western section of the study area, while the eastern section exhibited low Cu content.
The distribution of Zn presented a different pattern compared to the other trace elements, with higher concentrations at the central and western section of the study zone. The maximum surficial concentrations were observed at the center of the study area, next to the AMD treatment system. This location was considered the point of discharge for all the superficial runoff coming from this area of the study zone, and it showed signs of being periodically flooded. The high PTE content could, therefore, be understood as the consequence of lixiviation and mobilization of the highly mobile element Zn from upper topographic positions. Other areas of high Zn content were scattered through the western section of the study area, with no clear correlation between them.
Although not considered as a PTE, the S (SO 3 ) content also showed a similar spatial distribution as Zn, with higher values located at runoff discharge points at the center of the study area and near the pollution source at the upper sections of the ephemeral channels. In this case, the S content was indicative of the presence of sulfates in soils, which is typical in areas affected by AMD processes.  Figure 4. Isocontent maps of total Fe2O3, SO3, and trace-element concentrations obtained from surficial soil samples. 487 As and Pb showed a similar surficial distribution, with higher values at the upper 488 section of transect W and transect E. These two trace elements had notably distinct mobi-489 lization conditions. Therefore, it is reasonable to think that the coincidence in their distri-490 bution pattern correlated to the proximity to the source of PTE contamination (in this case, 491 the encapsulated mine wastes at the top of the deposit) and not to the mobilization pro-492 cesses. Cu concentration in soil samples followed a similar pattern, with higher values on 493 upper topographic positions. In this case, however, maximum concentration values cor-494 responded to transect W and the western section of the study area, while the eastern sec-495 tion exhibited low Cu content. 496
No significant differences in mineralogical composition were observed between samples collected in upper, mid, and lower sections of the ephemeral channels, with the exception of high secondary gypsum content in the LOS13 sample and a slightly higher jarosite/plumbojarosite and alunite content in some of the samples collected in upper topographic positions. Nevertheless, the total content of soluble secondary sulfates (jarosite, plumbojarosite, and alunite) in the surficial soil samples was still considerably low (up to 9%), which may indicate the presence of active dissolution and mobilization process of sulfates in the study area.

Multivariate Analysis
A geostatistical multivariate analysis was carried in order to confirm possible relationships between the total PTE and mineralogical content of the surficial soil samples (LOS1 to LOS17). The Pearson correlation analysis demonstrated a positive correlation between the jarosite/plumbojarosite content and Fe 2 O 3 content (Pearson correlation coefficient, R = 0.9451), between SO 3 and jarosite/plumbojarosite content (R = 09344), and between Pb and jarosite/plumbojarosite content (R = 0.9469). No other relevant positive or negative correlations were found in the surficial soil samples.
A factor analysis for the PTE content and soil mineralogy was also carried out to help diagnose relationships between variables ( Table 4). The first four extracted factors were the only ones to display eigenvalues higher than 1. Using the four-factor model, the final factor solution represented 79.97% of the variance in the analyzed data. Factor 1 (F1) contained Fe 2 O 3 , SO 3 , Pb, As, and jarosite/plumbojarosite, all with a positive varimax rotated factor, as well as SiO 2 and feldspars with a negative factor. This factor grouped the elements and mineralogy commonly related to AMD processes (Fe 2 O 3 , SO 3 , Pb, As, and jarosite/plumbojarosite) and the components of the soil samples related to the host rock (SiO 2 and feldspars). Factor 2 included alteration minerals such as phyllosilicates and goethite (positive factor), as well as quartz and Cd (negative factor). The presence of Cd in F2 may not be accurate due to the very low Cd content detected in all soil samples. Factor 3 comprised alunite and gypsum (positive factor), which are both soluble sulfate phases. Lastly, F4 contained Zn and Cu (positive factor), with no clear correlation between them.  Figure 5 shows the score plots for F2-F1 and F3-F1. In the first case, samples with relative low quartz content plotted above the Y-axis (F2). Only a few of the analyzed samples (LOS1, LOS2, LOS3, and LOS17) showed relative high values of trace elements and mineral phases associated with AMD (Fe 2 O 3 , SO 3 , Pb, As, and jarosite/plumbojarosite) in comparison to the host rock phases (feldspars, SiO 2 -quartz). In the score plot for F3-F1, most of the samples plotted in the lower left margin of the graphic, which exemplified a low soluble sulfate content (alunite, gypsum, and jarosite/plumbojarosite), in agreement with the results obtained from the mineralogical characterization of the soil samples.   Figure 5 shows the score plots for F2-F1 and F3-F1. In the first case, samples w relative low quartz content plotted above the Y-axis (F2). Only a few of the analyzed sam ples (LOS1, LOS2, LOS3, and LOS17) showed relative high values of trace elements a mineral phases associated with AMD (Fe2O3, SO3, Pb, As, and jarosite/plumbojarosite) comparison to the host rock phases (feldspars, SiO2-quartz). In the score plot for F3-F most of the samples plotted in the lower left margin of the graphic, which exemplified low soluble sulfate content (alunite, gypsum, and jarosite/plumbojarosite), in agreeme with the results obtained from the mineralogical characterization of the soil samples.

Water Extraction
The mobility of PTEs under natural conditions varied greatly between each of t analyzed trace elements (

Water Extraction
The mobility of PTEs under natural conditions varied greatly between each of the analyzed trace elements ( Table 5). The obtained results show that the highest concentration values in the soluble extraction corresponded to Zn, with a mean value of 69.8 mg·kg −1 , while Cu presented intermediate values around 23.4 mg·kg −1 , with highly contrasting maximum and minimum values of 99.6 mg·kg −1 and 0.1 mg·kg −1 , respectively. The As extracted exhibited low mobility rates under natural conditions, with an average value of 0.04 mg·kg −1 of extractable As. Both Pb and Cd concentrations appeared below the analytical detection limit in all the measured samples; in the case of Pb, this was an indication of the low mobility of the element, whereas, for Cd, this fact was related to the very low total concentrations found.
The water extraction percentage results show that release of Zn was considerably high, with an average value of 23.4% total Zn mobilized from the original soil samples. The extraction percentages for the other PTEs were notably lower, however. Both Cu and As presented the lowest quantified absolute values of extraction percentage after water extraction, with minimum values of 0.03‰ and 0.002‰ of total Cu and As, respectively. In the cases of Pb and Cd, the extraction percentages could not be accurately estimated, due to the concentration values appearing below the analytical detection limit.
In general, there were no significant variations in the mineralogical composition of the samples after the water extraction ( Figure 6), with the exception of a slight increase in some mineral phases, such as quartz and phyllosilicates, as a result of the disappearance of gypsum in some samples, i.e., samples LOS8 and LOS11. On the other hand, Cu content in the soluble extract appeared lower in comparison 639 with the other extraction procedures, with an average value of 11.5 mg·kg −1 . As with the 640 other extractions, Cd content could not be quantified and appeared below the analytical 641 detection limit. 642 With regard to the extraction percentages results, PTE mobilization was considerably 643 high for Zn (30% of total Zn content) and As (mean value of 22.8%), while the Cu extrac-644 tion percentage was relatively low, with an average value of 3%. In contrast to the previ-645

Acidic Medium Extraction
Mobilization of PTEs under acidic conditions showed relatively similar results to the water extraction values (Table 5). In conditions comparable to AMD, Zn appeared as the most mobile PTE, with an average value of 72.6 mg·kg −1 in the extraction solution. When submitted to the acidic medium extraction, Cu showed the highest concentration value in comparison to other extracting conditions, with values ranging from 104.3 mg·kg −1 to 10.4 mg·kg −1 , and an average value of 42.8 mg·kg −1 . The As content in the soluble fraction after the acidic medium extraction was the lowest of the quantified PTEs, with its maximum value reaching only 4.048 mg·kg −1 (mean value of 1.3 mg·kg −1 ). Similarly to the water extraction, Pb and Cd appeared below the analytical detection limit in all samples.
Regarding the extraction percentage results, Zn showed the highest mobilization percentage, with 24.7% of total Zn content, followed by Cu, with a mean value of 10.9%. The results showed that, under acidic conditions, As was expected to be almost immobile, with an average extraction percentage of 0.29%. The extraction percentages of Pb and Cd could not be estimated, due to the lack of precisely quantified values.
The mineralogy of residues after the extraction showed the complete disappearance of water-soluble phases such as gypsum in samples LOS8, LOS11, and LOS13. The rest of the mineral phases remained unaffected, with a small relative increase in concentration due to the decrease of the aforementioned soluble phases.

Complexing-Reducing Medium
Amongst the five extraction procedures, soil samples subjected to the complexingreducing medium extraction showed the highest PTE mobility regarding As, Zn, and Pb content, with average concentration values of 111.1 mg·kg −1 , 84.5 mg·kg −1 , and 55.8 mg·kg −1 , respectively (Table 5). Regarding the Pb content, due to the fact that some samples showed Pb content below analytical detection limit, we considered the Pb content of those samples as half the detection limit value (DL/2 = 0.75 mg·kg −1 ), in order to estimate its mean value. By doing this, samples with low extractable Pb content could also be taken into consideration when estimating the mean value.
On the other hand, Cu content in the soluble extract appeared lower in comparison with the other extraction procedures, with an average value of 11.5 mg·kg −1 . As with the other extractions, Cd content could not be quantified and appeared below the analytical detection limit.
With regard to the extraction percentages results, PTE mobilization was considerably high for Zn (30% of total Zn content) and As (mean value of 22.8%), while the Cu extraction percentage was relatively low, with an average value of 3%. In contrast to the previous extraction procedures, quantifiable lead content was released under complexing-reducing conditions in some samples. The extractable Pb released in this medium reached notably high values, with some samples exhibiting extraction percentages above 50%, reaching a maximum value of 95.27% total Pb content.
X-ray diffraction results showed a complete depletion of hydrated and secondary sulfates, such as gypsum, jarosite, and alunite, in all samples except for sample LOS17, where the jarosite/plumbojarosite content was reduced from 32 wt.% to 8 wt.%. Fe 3+ oxyhydroxides such as goethite were also absent in the residual sediment samples. The disappearance of jarosite/plumbojarosite and goethite could be interpreted as a consequence of the reduction of ion Fe 3+ to Fe 2+ , leading to the dissolution of those mineral phases. Quartz, phyllosilicate, and feldspar contents were not affected or showed a relative increase in response to the disappearance of the previously mentioned mineral phases.

EDTA
Mobilization of PTEs in soils after the EDTA extraction was slightly lower than that for other extraction methods, especially for Zn, with an average concentration of 50.5 mg·kg −1 ( Table 5). It is worth noting that, in contrast to the other extraction procedures, Zn content appeared below the analytical detection limit in some of the analyzed samples. Cu showed moderate to low concentration values, with an average of 37.2 mg·kg −1 . Extractable As content was slightly lower than that detected in acidic medium extraction, with values between 1.308 mg·kg −1 and 0.034 mg·kg −1 . In all considered samples, extractable Pb and Cd appeared below the analytical detection limit.
Extraction percentages were generally low for all studied PTEs, with Zn and Cu as the trace elements that presented higher mobilization percentages, with 17% and 8.7% of total PTE content, respectively. Arsenic exhibited an almost immobile behavior in conditions comparable to plant uptake, with extraction percentages that did not exceed 0.341% of total As. The extraction percentages of Pb and Cd could not be precisely estimated; therefore, we chose not to include them, in order to avoid unnecessary error.
The mineralogical content in the residue after the EDTA extraction was almost identical to that of the original soil samples. The only mineral phase that showed a noticeable variation was gypsum, with a decrease in concentration. A relative increase in some phases, such as goethite and feldspars, could take place, as a result of the decrease of the hydrated sulfates.

DTPA
Results obtained from soils submitted to the DTPA extraction were similar to the EDTA extraction (Table 5). Cu content in the extraction solution was almost identical to the EDTA results, with an average concentration of 37.5 mg·kg −1 . Zn and As showed slightly higher values, with 65.1 mg·kg −1 and 0.5 mg·kg −1 , respectively, suggesting that the DTPA extraction procedure was able to release a higher amount of PTEs when used in strongly acidic soils. Moreover, the DTPA extraction procedure was the only methodology used that was able to extract enough Pb in all soil samples to be quantified above the analytical detection limit, with values ranging from 9.1 mg·kg −1 to 0.56 mg·kg −1 .
Once the extractable PTE content was converted into percentages, Zn showed the highest values (20.7% of total Zn), followed by Cu (9%). For Pb and As, the percentages were notably low, ranging from 0.03% to 5.16% (average value of 1.8%) and from 0.294% to 0.010%, respectively.
The obtained mineralogical data were very similar to the original soil samples, with a small decrease or disappearance of gypsum in samples LOS13, LOS8, and LOS11, and with a consequent relative increase in the other mineral phases. Goethite and jarosite/plumbojarosite still appeared in most of the samples, but with a slight reduction in content.
Results obtained from the application of the five selective extraction procedures are summarized in Figures 7 and 8. Cu showed moderate to low concentration values, with an average of 37.2 mg·kg −1 . Ex-663 tractable As content was slightly lower than that detected in acidic medium extraction, 664 with values between 1.308 mg·kg −1 and 0.034 mg·kg −1 . In all considered samples, extracta-665 ble Pb and Cd appeared below the analytical detection limit. 666 Extraction percentages were generally low for all studied PTEs, with Zn and Cu as 667 the trace elements that presented higher mobilization percentages, with 17% and 8.7% of 668 total PTE content, respectively. Arsenic exhibited an almost immobile behavior in condi-669 tions comparable to plant uptake, with extraction percentages that did not exceed 0.341% 670 of total As. The extraction percentages of Pb and Cd could not be precisely estimated; 671 therefore, we chose not to include them, in order to avoid unnecessary error. 672 The mineralogical content in the residue after the EDTA extraction was almost iden-673 tical to that of the original soil samples. The only mineral phase that showed a noticeable 674 variation was gypsum, with a decrease in concentration. A relative increase in some 675 phases, such as goethite and feldspars, could take place, as a result of the decrease of the 676 hydrated sulfates. Results obtained from soils submitted to the DTPA extraction were similar to the 679 EDTA extraction (Table 5). Cu content in the extraction solution was almost identical to 680 the EDTA results, with an average concentration of 37.5 mg·kg −1 . Zn and As showed 681 slightly higher values, with 65.1 mg·kg −1 and 0.5 mg·kg −1 , respectively, suggesting that the 682 DTPA extraction procedure was able to release a higher amount of PTEs when used in 683 strongly acidic soils. Moreover, the DTPA extraction procedure was the only methodol-684 ogy used that was able to extract enough Pb in all soil samples to be quantified above the 685 analytical detection limit, with values ranging from 9.1 mg·kg −1 to 0.56 mg·kg −1 . 686 Once the extractable PTE content was converted into percentages, Zn showed the 687 highest values (20.7% of total Zn), followed by Cu (9%). For Pb and As, the percentages 688 were notably low, ranging from 0.03% to 5.16% (average value of 1.8%) and from 0.294% 689 to 0.010%, respectively. 690 The obtained mineralogical data were very similar to the original soil samples, with 691 a small decrease or disappearance of gypsum in samples LOS13, LOS8, and LOS11, and 692 with a consequent relative increase in the other mineral phases. Goethite and jaro-693 site/plumbojarosite still appeared in most of the samples, but with a slight reduction in 694 content. 695 Results obtained from the application of the five selective extraction procedures are 696 summarized in Figures 7 and 8. In conclusion, the PTE release after the water extraction depended on the natural 701 mobility of the trace elements and mineralogical composition of samples, affecting only 702 soluble phases, mainly hydrated sulfates (Figure 8). 703 704 Figure 8. Summary of affected mineral phases and PTEs extracted in the studied soil samples. 705 PTEs with high mobility, such as Zn, showed higher extraction percentages, while 706 As and Pb were minimally affected. PTEs extracted after the acidic medium procedure 707 showed slightly higher mobilization percentages as compared to the water extraction. The 708 lower pH of this medium favored metal mobilization and attacked soluble sulfates and 709 carbonates. However, the release of PTEs under those conditions was not especially high, 710 due to the lack of a carbonate-associated PTE fraction in the sediment samples. In the 711 complexing-reducing medium, hydrated sulfates, jarosite/plumbojarosite, alunite, and 712 goethite were attacked. In this case, PTE release was especially high for As and Pb, which 713 means that those elements were probably related to the aforementioned mineralogical 714 phases. Specifically, most of the As content in soils in the Lousal mine samples appeared 715 to be adsorbed to goethite, while Pb contents were probably related to secondary sulfates 716 such as plumbojarosite. 717 Lastly, the bioavailable fraction assessed with the EDTA and DTPA procedures was 718 low for all of the PTEs analyzed. Zn and Cu showed medium to low extraction percent-719 ages, and Pb and As exhibited an almost immobile behavior under conditions similar to 720 plant uptake. The low Cd content obtained in all the extraction procedures (below the 721 analytical detection limit) was related to the low total Cd content in the original samples. 722

723
The contaminant factor (CF) values for As and Cu were very high in almost all the 724 soil samples analyzed (Tables 6 and 7). In general, both surficial and deep soil samples 725 showed similar CF values for those PTEs, with only a slightly higher CF for As in surficial 726 samples. Regarding Cu, values of CF varied between considerably contaminated and very 727 contaminated for almost the entirety of the samples and showed a slightly lower standard 728 deviation than obtained for As. In the case of Pb and Zn, CF values were lower than those 729 determined for As and Cu and varied from low contamination to very high contamina-730 tion. CF values obtained for Zn were fairly constant throughout all analyzed samples, 731 without surpassing values of CF = 6, suggesting very low to moderate contamination. In conclusion, the PTE release after the water extraction depended on the natural mobility of the trace elements and mineralogical composition of samples, affecting only soluble phases, mainly hydrated sulfates ( Figure 8).
PTEs with high mobility, such as Zn, showed higher extraction percentages, while As and Pb were minimally affected. PTEs extracted after the acidic medium procedure showed slightly higher mobilization percentages as compared to the water extraction. The lower pH of this medium favored metal mobilization and attacked soluble sulfates and carbonates. However, the release of PTEs under those conditions was not especially high, due to the lack of a carbonate-associated PTE fraction in the sediment samples. In the complexing-reducing medium, hydrated sulfates, jarosite/plumbojarosite, alunite, and goethite were attacked. In this case, PTE release was especially high for As and Pb, which means that those elements were probably related to the aforementioned mineralogical phases. Specifically, most of the As content in soils in the Lousal mine samples appeared to be adsorbed to goethite, while Pb contents were probably related to secondary sulfates such as plumbojarosite.
Lastly, the bioavailable fraction assessed with the EDTA and DTPA procedures was low for all of the PTEs analyzed. Zn and Cu showed medium to low extraction percentages, and Pb and As exhibited an almost immobile behavior under conditions similar to plant uptake. The low Cd content obtained in all the extraction procedures (below the analytical detection limit) was related to the low total Cd content in the original samples.

Assessment of Potential Environmental Risk
The contaminant factor (CF) values for As and Cu were very high in almost all the soil samples analyzed (Tables 6 and 7). In general, both surficial and deep soil samples showed similar CF values for those PTEs, with only a slightly higher CF for As in surficial samples. Regarding Cu, values of CF varied between considerably contaminated and very contaminated for almost the entirety of the samples and showed a slightly lower standard deviation than obtained for As. In the case of Pb and Zn, CF values were lower than those determined for As and Cu and varied from low contamination to very high contamination. CF values obtained for Zn were fairly constant throughout all analyzed samples, without surpassing values of CF = 6, suggesting very low to moderate contamination. In the case of Pb, the highest values were quantified in samples located at upper positions in the slope, near the encapsulated materials (LO2.1 and LOS6), and in samples collected from areas of sediment accumulation (LO1.4 and LOS17), showing a very high degree of contamination for this element in those areas. The contaminant factor for Cd varied from low to considerable contamination in some of the surficial samples, whereas it could not be estimated in the remaining samples, due to the total concentration values being below the analytical detection limit. Table 6. Contamination factor (CF) and pollution load index (PLI) for Pb, Zn, Cd, As, and Cu in surficial soil samples.
According to the pollution load index results (PLI) for the selected elements, the soils affected by mine activity in the Lousal mine presented contamination varying from a considerable to a very high degree. Similarly to the CF results, the highest degree of contamination was found in areas near the encapsulated materials and at the discharge points of superficial runoff coming from the slope, which showed signs of being periodically flooded.
The results obtained after applying the five proposed extraction procedures indicate that the highest PTE mobilization took place in a complexing-reducing medium. These results suggest that, if the overall environmental conditions of the Lousal mine wastes and soils changed from an oxidizing medium to a reducing medium, mobilization of PTEs such as As and Pb could be an important issue. As seen from the geochemical characterization of the original sediment samples, the spatial distribution of As and Pb content in soils did not seem to be influenced by active mobilization processes. Under the current oxidant conditions, those PTEs exhibited low mobility, indicating that they may be associated with mineral phases that are stable in an acidic and oxidizing medium, such as goethite and plumbojarosite [7,74]. However, a change to reducing conditions, which can be caused by the partial encapsulation of the mine wastes and soils or by bad drainage, could lead to the dissolution of the Fe 3+ oxyhydroxides and secondary sulfates, thus favoring the release of As and Pb into the environment and nearby water courses [6,75].
PTE extraction percentages obtained after acidic medium extraction were similar to those obtained by water extraction, suggesting that the trace-element mobilization promoted by AMD did not seem to have a significant impact in the study area. However, although the acidity and chemical characteristics of runoff water did not play a major role in the release, transport, and accumulation of PTEs, the preferential surface runoff flow paths and areas of water and sediment accumulation seemed to condition the spatial distribution of certain trace elements. For example, PTEs with higher natural mobility such as Zn and ions such as SO 3 associated with water-soluble phases appeared in especially high concentrations at discharge areas of surficial water runoff and sediment deposition. Remediation strategies focusing on the management of those discharge areas and preferential runoff flow paths coupled with adequate chemical treatment actions would allow for an easier control of PTEs with high natural mobility. Regarding the pollution impact on the vegetal communities developed in the analyzed soils of the Lousal mining area, the data obtained for the PTEs content assimilable by vegetation (EDTA and DTPA extractions) show that PTE plant uptake does not seem to be a huge concern in the study area. Field observations agree with those results (Figure 9), as PTE content in soils did not seem to heavily affect plant growth in the mine wastes, aside from the limitation imposed by the acidic pH on the plant community, restricting the diversity of vegetation to species tolerable to acidic soils, such as the genera Cistus and Acacia. Moreover, the existence of a well-developed plant community favored the presence of oxidant conditions in the superficial soil layer, limiting the mobility of contaminants such as As and Pb.

787
A combined fluvial-geomorphic and geochemical reclamation approach will be ap-788 plied in the Lousal mining area throughout 2021-2022, within the framework of the Euro-789 pean-funded mine reclamation project LIFE RIBERMINE (LIFE ENV/ES/000181). The to-790 pography of the affected mine soils will be redesigned and reconstructed, and a soil 791 amendment will be extended over the new topography. In light of the results obtained 792 from the geochemical characterization and PTE mobility assessment of the soil samples, 793 Figure 9. View of the study area, with shrub vegetation (Cistus ladanifer), Acacia, and Eucalyptus growing in soils affected by mining activity.

Application to Metallic Mine Reclamation
A combined fluvial-geomorphic and geochemical reclamation approach will be applied in the Lousal mining area throughout 2021-2022, within the framework of the European-funded mine reclamation project LIFE RIBERMINE (LIFE ENV/ES/000181). The topography of the affected mine soils will be redesigned and reconstructed, and a soil amendment will be extended over the new topography. In light of the results obtained from the geochemical characterization and PTE mobility assessment of the soil samples, if these reclamation strategies are bound to succeed, the risk of introducing reducing conditions with those reclamation actions must be taken into account and prevented.
The GeoFluv method will be applied in order to redesign the original topography. The GeoFluv method is a fluvial-geomorphic reclamation method that aims to reconstruct mature landforms and drainage networks that would develop within a natural catchment, after thousands of years of work performed by geomorphic processes [49][50][51]. This reclamation method essentially "compresses time", building geomorphically stable and functional drainage networks that are able to efficiently manage the surficial runoff of the reclaimed area. The resulting landscape is composed of a system of ridges and swales that drain the runoff water into meandering channels, similarly to the natural landscapes in the surrounding areas.
In the Lousal mining area, this will be implemented with the reconstruction of two small drainage networks (Figure 10), which will convey surficial water flow into a series of constructed channels. GeoFluv landscape, featuring two small drainage networks and a system of ridges, swales, and me-810 andering channels to convey surficial runoff. 811 The new reclamation landscape would favor surficial water flow within preferential 812 paths (swales and channels), thus preventing the risk of formation of periodically flooded 813 areas that would introduce reducing conditions and accelerated mobilization of As and 814 Pb. This reclamation approach would also facilitate the control and treatment of other 815 contaminants with high natural mobility, such as Zn or Cd. The management of the surf-816 The new reclamation landscape would favor surficial water flow within preferential paths (swales and channels), thus preventing the risk of formation of periodically flooded areas that would introduce reducing conditions and accelerated mobilization of As and Pb. This reclamation approach would also facilitate the control and treatment of other contaminants with high natural mobility, such as Zn or Cd. The management of the surficial drainage though designed flow paths allows focalizing the chemical treatment actions aimed at neutralizing surficial waters and removing those contaminants that are more easily mobilized by water and particle transport processes. In this regard, limestone blocks will be extended over the riverbed of the designed channels in order to increase the pH of surficial waters and favor the precipitation of heavy metals [76] (Figure 11).
Geosciences 2021, 11, x FOR PEER REVIEW 2 tration and subsuperficial water circulation should prevent water stagnation and, t fore, maintain oxidant conditions throughout the soil profile, reducing the risk of A Pb mobilization. Compaction of the surficial soil horizon should be avoided in ord guarantee a good water flow through the buffering soil amendment. Figure 11. Schematic view of the limestone channels and soil amendment designed with framework of the LIFE RIBERMINE reclamation project.

Conclusions
The Lousal mining area was exploited from 1900 to 1988, resulting in the produ of a large volume of pyrite-rich mine wastes, ranging from barren materials to va types of waste deposits and tailing impoundments. The geochemical and mineralo characterization of the studied area shows that the mine wastes underwent intense w ering processes, producing an important contamination of the adjacent soils, which led to the release and mobilization of PTEs into nearby water courses. Total PTE re indicate that the soils affected by mining activities were highly contaminated with A Cu, while Zn and Pb content ranged from low to very high, depending on the anal samples. Cadmium levels were found to be very low in most of the soil samples.
The analysis of the spatial distribution of PTE content, together with the sele extraction procedures, proved to be powerful tools for the general study of PTE mo and pollution potential in areas affected by mining activities. The results show that with high natural mobility, such as Zn, and elements related to soluble phases (gyp and secondary sulfates) such as S appear concentrated at water and sediment disch areas, while As, Pb, and Cu contents were generally higher in high topographic posi of the slope, near the contamination source.
The results obtained after applying the proposed extraction methodologies ind that the highest mobilization took place in the complexing-reducing medium. This pecially evident in the case of otherwise near-immobile elements such as As and Pb, ciated with mineralogical phases that are stable in oxidizing and acidic conditions, as goethite and jarosite/plumbojarosite. These results allow suggesting that, if a ch from oxidant to reducing conditions takes place, PTE mobilization could be an impo A buffering soil amendment will be extended over the new reclamation topography. The amendment will be composed of two distinct soil horizons with specific purposes: to provide a good foundation for plant growth and to favor the chemical stabilization of contaminants ( Figure 11). The most surficial horizon (horizon A, 10 cm of thickness) will be composed of poultry manure, topsoil, and non-expansive clays, while the lower horizon (horizon B, 15 cm of thickness) will be made of a mixture of clays and limestone gravel. This horizon will help neutralize infiltrating water and, most importantly, facilitate subsuperficial water flow, due to the high limestone gravel content. Enhanced infiltration and subsuperficial water circulation should prevent water stagnation and, therefore, maintain oxidant conditions throughout the soil profile, reducing the risk of As and Pb mobilization. Compaction of the surficial soil horizon should be avoided in order to guarantee a good water flow through the buffering soil amendment.

Conclusions
The Lousal mining area was exploited from 1900 to 1988, resulting in the production of a large volume of pyrite-rich mine wastes, ranging from barren materials to various types of waste deposits and tailing impoundments. The geochemical and mineralogical characterization of the studied area shows that the mine wastes underwent intense weathering processes, producing an important contamination of the adjacent soils, which also led to the release and mobilization of PTEs into nearby water courses. Total PTE results indicate that the soils affected by mining activities were highly contaminated with As and Cu, while Zn and Pb content ranged from low to very high, depending on the analyzed samples. Cadmium levels were found to be very low in most of the soil samples.
The analysis of the spatial distribution of PTE content, together with the selective extraction procedures, proved to be powerful tools for the general study of PTE mobility and pollution potential in areas affected by mining activities. The results show that PTEs with high natural mobility, such as Zn, and elements related to soluble phases (gypsum and secondary sulfates) such as S appear concentrated at water and sediment discharge areas, while As, Pb, and Cu contents were generally higher in high topographic positions of the slope, near the contamination source.
The results obtained after applying the proposed extraction methodologies indicate that the highest mobilization took place in the complexing-reducing medium. This is especially evident in the case of otherwise near-immobile elements such as As and Pb, associated with mineralogical phases that are stable in oxidizing and acidic conditions, such as goethite and jarosite/plumbojarosite. These results allow suggesting that, if a change from oxidant to reducing conditions takes place, PTE mobilization could be an important environmental concern. Reclamation actions sought to be applied to these soils, including fluvial-geomorphic reclamation approaches and the extension of a buffering soil amendment, will favor an efficient water drainage, infiltration, and subsuperficial water circulation, thus maintaining the overall oxidant conditions in the soil.
Geo-environmental characterization, reclamation, and adaptation of closed and abandoned mines are fundamental in order to revalorize the mining heritage, particularly in regions with a long mining history, which in turn brings significant educational, environmental, socioeconomic, and eco-touristic benefits [77]. With this study, we provided an initial geochemical framework that can be directly applied to the planification and development of an effective reclamation strategy of the Lousal mine wastes. This methodological procedure is not restricted to this specific area, but can be implemented elsewhere in order to provide the necessary information about the geochemical characteristics and constraints of the mine wastes and soils affected by metallic mining, before developing an environmentally sound reclamation strategy.
Author Contributions: R.S.-D. was involved with the conceptualization of the study, contributed to the acquisition, treatment and analysis of data, wrote the original manuscript and the ongoing revised versions; M.L.G.L. contributed with the conceptualization of the study, was involved in the acquisition, treatment and analysis of data and provided ongoing reviews to paper drafts; J.M.E. participated with the conceptualization of the study, was involved in the acquisition, treatment and analysis of data and provided ongoing reviews to paper drafts; E.M.G.-N. was involved in acquisition and interpretation of laboratory data; P.H. contributed with the conceptualization of the study, was involved in the acquisition, treatment and analysis of data and provided ongoing reviews to paper drafts; E.C. was involved with the conceptualization of the study, was involved in the acquisition, treatment and analysis of data and provided ongoing reviews to paper drafts. All authors have read and agreed to the published version of the manuscript.