4.1. Metal Partition from AMD to Solid Phases
During the AMD treatment, the progressive increase of pH along the DAS-Ca column promoted the formation of an “ochre horizon” in the first centimeters of the solid residue followed by a “white horizon” (
Figure 1a). Samples Ca1, Ca2, and Ca3 were taken from the first centimeters of the profile (3.5, 7.5, 8.5 cm, respectively), while samples Ca4 and Ca5 were collected in the white horizon (10 and 15 cm, respectively). The chemical conditions developed along the column, triggered Fe depletion from solution (Delgado et al., 2019), probably forming low-crystalline oxyhydroxides or oxy-hydroxysulfates. such as H-jarosite [(H3O)Fe3(SO4)2(OH)6], K-jarosite [KFe3(SO4)2(OH)6], schwertmannite [Fe8O8(OH)6(SO4)], ferrihydrite [Fe(OH)3], hematite [α-Fe2O3], goethite [α-FeOOH], lepidocrocite [γ-FeOOH], and maghemite [γ-Fe2O3].
Indeed, the “ochre horizon” (
Figure 2a) exhibited the highest concentrations of Fe, with retention rates exceeding 90% (Ca1; 336 g/kg) and 80% (Ca3; 149 g/kg) of the total metal content (
Figure 2a;
Table 2), as confirmed by the corresponding visual observations. Lower Fe concentrations were measured in sample Ca2 (30 g/kg), accounting for 10% of the total element proportion. This can be attributed to significant dissolution of the calcium substrate and unexpected local alkaline conditions in this section of the DAS-Ca (
Figure 1c). This is further supported by the high concentrations of Cu (42 g/kg) observed in this horizon profile, indicated by the presence of green precipitates (
Figure 1c). These AMDp are commonly associated with the precipitation of divalent metal hydroxides, as described by [
32], when the pH reaches neutral to alkaline conditions. This supports the high concentrations of Al (43 g/kg), Mg (1300 mg/kg), Mn (1223 mg/kg), and Zn (3860 mg/kg) in this layer of the profile.
The deeper ‘white horizon’ likely corresponds to Al precipitates formed when pH increases, as described by [
14]. This is confirmed by the higher Al concentration through the column, with the highest value (56 g/kg) present in sample Ca4 (
Table 2).
SEP data of the surficial sample, Ca1, indicates the presence of high concentrations of Ca (760 mg/kg). Besides, Ca increases at greater depths, so Ca3 shows 17 g/kg associated with the labile fraction of the residues (
Figure 2b,d, respectively). Furthermore, in the exchangeable fraction (F2), the release of low crystalline iron hydroxide seems evident, along with weakly adsorbed trace elements onto these compounds, as previously described other authors [
40]. The presence of iron and trace metals in the F2 of sample Ca1 (
Figure 2b) endorse this interpretation.
The movement of trace elements from solution to solid Fe-Al phases involves adsorption and/or coprecipitation mechanisms. Research on ochre-like deposits rich in iron in acidic conditions indicates that jarosite and schwertmannite may have the greatest amounts of As and Cr [
41], as well as other metals like Cu, Pb, and Zn [
42]. Furthermore, ferrihydrite appears to be related to high concentrations of Pb and relatively high content of Zn and Ni [
43]. In this regard, the main concentration of As associated with the higher concentration of Fe and K in samples Ca1 and Ca3 (
Figure 2b,d) was determined, presumably associated to schwertmannite and jarosite. Contrary to expectations, although significant concentrations of As, Cu, and Pb associated with fractions F3 and F4 (involving Fe(III) oxyhydroxides and oxyhydroxysulfates) were leached, concentrations of up to 2536 mg/kg of As (e.g. Ca1 sample) were extracted in F2, the fraction associated with carbonates (
Figure 2b). This could be attributed to the dissolution of amorphous or not ordered Fe(III) precipitates, along with the release of adsorbed elements such as As, Cu, Mn, and Zn, induced by the exposure to acetic acid at pH 4.5 within fraction F2. Additionally, excluding the samples from the DAS-Mg, trace metals Cu, Cd, Pb and Zn showed considerable concentrations associated with this section of the iron profile (up to 105, 38.1, 27.7 and 42.2 mg/kg, respectively), confirming the link between Fe-rich mineral phases and trace elements retention.
According to [
15], aluminum may be trapped as hydroxides or hydroxysulfates in the white precipitates, which were mainly concentrated in sample Ca4, as shown in
Figure 1a. The high Al concentrations extracted in F3 and F4 confirmed the dissolution of Al-rich sulphate phases such as alunite [KAl3(SO4)2(OH)6], basaluminite [Al4(SO4)(OH)10·5(H2O)], jurbanite [Al(SO4)(OH)·5(H2O)] or Al-hydroxides [Al(OH)3(am), gibbsite AlOH3, boehmite γ-AlO(OH), and diaspore AlO(OH)] that would be important as the pH rises above 6 [44-45; 25]. Furthermore, the presence of alunite could be elucidated by considering the relative concentration of potassium in this profile section. Also, elevated concentrations of Cu, Ca, and Al detected in sample Ca2 seem to support this statement. These precipitates frequently contain gypsum (confirmed by the high concentration of Ca associated to F1 fraction,
Figure 2c,e). Furthermore, according to the hydrochemical model developed by [
2] in these sections of the profile, the dominant aqueous species (AlSO4+, Al(OH)2+ and AlOH+2) and the pH range between 5.8 and 6.4 are consistent with the precipitation of Al-phases.
It is noted that the neutral conditions obtained in the final section of the DAS-Ca are not usually enough for the formation of divalent metal hydroxides [
32]. According to [
46], the most effective retention of divalent metals is achieved when the pH is between 8 and 10 because important desorption processes can occur at higher pH values. Based on these premises, the presence of high concentration of Ca at the final section of the DAS-Ca (F2 of samples Ca4 and Ca5,
Figure 2) added to the addition of alkalinity by the magnesia dissolved in the DAS-Mg (supported by the presence of Mg in the F1 of Ca5) appear to provide an ideal scenario for divalent metal retention. However, the control over the fractionation of divalent metals appears to be mainly exerted by carbonate phases [
2]. These conditions were obtained towards the end of DAS-Ca (sample Ca5), were accompanied by significant concentrations of Ca, Mg, Cu, or Fe (until 231, 2.3, 9.4, and 1.7 g/kg, respectively) associated with the exchangeable phase (
Figure 2f). Additionally, the fractionation of divalent elements in this zone of the profile seems to be somehow influenced by the availability of aluminum, with these phases extracted in F3 of the SEP, as previously suggested by [
35].
The samples at 1 and 12cm in the DAS-Mg (Mg1 and Mg2, respectively) exhibited the highest concentrations of Cd, Co, Cu, Ni, Pb and Zn (
Table 2 and
Figure 2g-h). The DAS-Mg hydrochemical model indicated that the primary carbonate species in water quickly vanish upon entering the DAS-Mg column because of the formation of carbonate minerals [
2], such as malachite [Cu3(CO3)2(OH)2], azurite and to a lesser extent smithsonite [ZnCO3:1H2O]. These premises are congruent with the results of the partitioning process observed in the experimental columns, where the main retention of Cu and Zn (900 and 516 g/kg, respectively) were associated with the sample Mg1 (exchangeable fraction,
Figure 2g), which was supported by a noticeable change in color (bluish green) of the precipitates at the DAS-Mg surface. It is also important to note that most partition of Zn from solution to the solid phase occurred in the final section of the DAS-Mg, with Mg1 (around 40%) and Mg2 (80%) being the most significant (
Figure 2a). In calcareous settings, when there are relatively high levels of divalent metals present (such as in the DAS- Mg column), the amount of these metals in water (including Cd and Mn) is typically controlled by the formation of solid metal-carbonate compounds through precipitation [
47,
48,
49]. In fact, Zn and Mn showed a clear association with the carbonate fraction of the AMDp. Also, while Cd exhibits some association with water-soluble fraction along with Co (894 and 1360 mg/kg, respectively), its higher concentrations (3558 and 3202 mg/kg) appear to be associated with the F2 fraction (
Figure 2h). Other divalent metal as Ni and Cr seems plot similar behavior associated with these fractions, showing their higher concentrations (e.g. Ni, 1878 mg/kg in F2,
Figure 2h).
4.2. DAS Mineralogical Assemblages
In the diffraction pattern of the surficial sample of DAS-Ca (Ca1), the presence of amorphous and low crystalline phases required background subtraction, so the peaks of some minerals become evident (
Figure 3a).
In the “iron profile”, the precipitation of the first oxyhydroxysulfates, such as schwertmannite and jarosite, could trap toxic elements like As, Cr, Cu, Pb, and Zn from solution [
41,
42]. According to [
50], these are the controlling phases at low pH values. Observations by SEM-EDS showed a clear association between S and Al-Fe, and dissemination of As and Pb trapped in this column section (
Figure 4a). Other sulphate species could play an important role here, depending on the availability of certain elements in the solution. The high proportion of Fe3+ in solution could led to coquimbite formation (Fe3+2(SO4)3·9H2O), and if the water contains more than 1500 mg/L sulfate, the presence of gypsum is expected due to the dissolution of the reactive calcium across the substrate. Additionally, brushite (CaHPO4·2H2O) seems to crystallize under these conditions [
51] in an AMD wetland restoration treatment. In this line, the hydrochemical data reported by [2; 26] suggest that during the initial acidity neutralization process at low pH, several phosphates could crystallize, including brushite (
Figure 3a) and other not detected like monazite (La, Ce) which could retain REE from the solution.
Sample Ca2 (
Figure S1, Supplementary Data), was dominated by amorphous phases, but also presented calcite and gypsum, which accounted 54% and 45% of the crystalline phases, respectively (
Table 3). This sample corresponded to a zone of the profile where the nature of the alkaline reagent could be observed. However, the relatively high concentration of Cu, Zn (
Figure 2c), and gypsum (rhombohedral) in this section seems to be a consequence of the beginning neutralization processes, as evidenced by the green-white color of the precipitates. In fact, the SEM images confirmed the presence of some Al-hydroxide amorphous precipitates and discrete Cu carbonated phases. Low presence of Fe associated to S and Al could explain the commonly described substitution Fe to Al in the Jarosite-Alunite group [
52,
53] as the pH increases (
Figure 4b-c).
Likewise, in Ca1 sample the background of Ca3 (8.5 cm deep profile) was subtracted to improve the identification of mineral phases. The results, consistent with the modelled hydrochemical data [
2], pointed the presence of schwertmannite (5.8%) and jarosite (6.9%), with presence of plumbojarosite (Pb0.5Fe3+3(SO4)2(OH)6) (
Figure S2). Although distinguishing among these phases is challenging [
41], these results were also checked with the chemical data (
Table 2). In this case, a preferential retention of Pb can be observed, supporting DR-X interpretations. Regarding trace elements, SEP results (F3 and F4 SEP extractant) and potassium distribution among the fractions could indicate a relative affinity between both As-schwertmannite and Cu-jarosite. Additionally, other sulphate phases such as melanterite (FeSO4·7H2O; 4%), gypsum (74%), or phosphates like brushite (5.9%) are linked to Fe and Ca reactivity in this section. The hydrochemistry of the water column in Ca3 and deeper strata suggested that Fe di-sulfate species were leading [FeSO4+ (73%) and Fe (SO4)2− (14%)], supporting the formation of melanterite, which would increase the proportion of hydrolyzed species, promoting the precipitation of saturated Fe phases like schwertmannite and jarosite [
54].
On the other hand, at this profile position, a low 2θ crest seems to indicate the presence of Al-phases like boehmite (3.4%,
Table 3). According to the hydrochemical model, at pH values greater than 4 (10cm depth), and in line with the data described by [
41], the presence of Al oxyhydroxides (alunite and basaluminite, among others) can be expected. As previously outlined, the substitution of Fe by Al becomes more noticeable in this section of the profile (
Figure 4c). In fact, SEM data allowed the distinction of a common pattern of distribution of Al, S, and Mg (
Figure 4d), which could imply the presence of other sulfate phases such as Epsomite (MgSO4) or pickeringite (MgAl2(SO4)4 2⸳2H2O). Although their presence cannot be confirmed, the SEP data also support this argument, as significant concentrations of Mg (280 mg/kg) associated with the labile fraction (F1) can be observed.
The samples at depths Ca4 and Ca5 did not show more Al-rich phases, except for the presence of diaspore, but in a low proportion (1.15% of the crystalline phases in sample Ca4). Both samples were characterized by the high presence of calcite from the reactive substrate, higher than 50% (
Table 3,
Figures S3, S4). Additionally, the Ca4 sample showed representative amounts of gypsum (44%) and dolomite (3%), probably derived from pH increase that could also promote the precipitation of this carbonate-type divalent metal phases in presence of high Mg concentration [
2]. Calcite and gypsum are the most frequent phases in the final section of the DAS-Ca and in the DAS-Mg column (
Table 3), due to the large amount of calcium-reactive available (evidenced by the high proportion in sample Ca5,
Figure 4e). In fact, SEM study of the surface precipitates in the DAS-Mg column suggested that the white aggregates (Mg-SB) were essentially composed of gypsum (
Figure 5a), while the green-blue ones (Mg-SV), also exhibited a nanocrystalline lattice composed of Fe, Cu and C (
Figure 5b). XRD suggested this grid was composed by metal(II) hydroxide carbonate minerals (rosasite group, including malachite-Cu) and siderite.
Under the conditions of the DAS-Mg column with pH > 6, newly formed carbonate and bicarbonate phases could play an important role in removing divalent metals from solution [
2]. The presence of dolomite in this section was explained by the dissolution of the reactive magnesia in the DAS-Mg. Thus, monohydrocalcite [CaCO3⸳(H2O)], malachite [Cu2CO3(OH)2], azurite [Cu3(CO3)2(OH)2], siderite and rhodochrosite [MnCO3] were identified (
Table 3). The efficient removal of Cu and Mn within the initial centimeters of the DAS-Mg (Mg1,
Figure 3b) was evident so subspherical Cu-C aggregates could be recognized by SEM (
Figure 5c).
Although iron chemistry is typically broad in more acidic environments (e.g., DAS-Ca), small amounts of available Fe associated with carbonate group as a divalent metal and precipitated as siderite [Fe2+(CO3)] (
Figure 3b and
Figure S5;
Table 3), representing 6.5% and 3.1% of the crystalline phases in samples Mg1 and Mg2, respectively. SEP results clearly revealed high concentrations of Cu, Fe, Mn, and Zn in the DAS-Mg (
Figure 2g and
Figure 2h) associated with carbonate system (F2). Additionally, based on the hydrochemical data of these environments (high Mn concentrations) and the mineralogical interpretations of the Das-Mg samples, the presence of rhodochrosite can be confirmed (
Figure 3b and 5S; 9% of the Mg2 sample).
According to [
32], metal hydroxides rapidly precipitate in these experiments, while sulphates precipitate more slowly and progressively. Although the literature suggests the presence of brucite [Mg(OH)2], it only appeared in a negligible proportion (0.6% over the total crystalline phases of the sample Mg2,
Figure S5) at the deepest part of the Mg-column. However, SEM data demonstrated the presence of Na and Mg associated, forming small-sized spherical aggregates (
Figure 5d). This was also supported by the concentration of Mg and Na determined in F1 of DAS-Mg samples (
Figure 2g-h), which evidenced carbonate-type divalent metal phases play an important role in the partition processes between AMD and solids in the final section of the treatment [
2]. Moreover, the alkaline conditions facilitate the saturation of carbonated phases, such as calcite or dolomite, which could be recognized as subspherical aggregates of Ca-Mg, displaying the same compositional pattern as C (
Figure 5d). These phases could play a role in coprecipitation-sorption processes, influencing the solubility of Zn, Mn, and other trace elements [
56] (Jensen et al., 2002). In fact, the higher Zn and Mn concentrations were associated with the F2 of the Mg1 and Mg2 samples (516 and 465 g/kg of Zn and 13 and 75 g/kg of Mn, respectively). The data revealed that Zn is the divalent metal whose partition towards the solid phase undergoes in the late stages of the treatment (DAS-Mg) but it is depleted by carbonate mineral phases such as aurichalcite [(Zn, Cu2+)5 (CO3)2⸳(OH)6] (7% in the Mg2 samples). SEP data showed that 516 g/kg of Zn, and 906 g/Kg of Cu were trapped in the carbonated phases (F2) of Mg1 sample, supporting this assumption.
On the other hand, due to the high availability of sulphate hydroxides, newly formed phases linked to divalent metals such as brochantite [Cu4SO4(OH)6], antlerite [(Cu2+)3SO4(OH)4], goslarite [ZnSO4·7H2O], and anglesite [PbSO4], among others, are frequently described [
47]. However, the mineralogical characterization of the samples is consistent with the hydrochemical model (
Figure 3b in [
2]) that shows low depletion of aqueous sulphated phases in the final section of the treatment. While hydrated sulphates are frequently described as the main controlling factor in the partitioning process of divalent metals such as Cu, Zn, or even Fe [
57], no sulphate phases involving these elements were identified in the samples of the present study.
Finally, the bottom sample of the treatment (“Arenas del Pindo” – sandy drainage;
Table 3 and
Figure S6) is characterized by the presence of quartz (72%), albite (8.7%), K-feldspar (microcline-orthoclase 2.5%), and phyllosilicates (mainly biotite 12% and chlorite, 4% of the total). This mineralogy explains the presence of crystalline quartz in some samples such as Ca1 or Mg2, probably accidentally transferred along the samples during the columns opening process.
4.3. Possible Environmental Implication of Residues Management
In the literature, inconsistent predictions regarding the stability of PTE (Potentially Toxic Elements) content in AMDp have been reported, depending on various factors such as the type of leaching test employed, the extracting agent used, the duration of the test, the presence of oxygen, and the regulatory criteria applied [
58]. In order to make more accurate predictions about how these residues will behave in the environment and to carry out an appropriate sustainable management option, metal stability in multi-step passive treatment systems require a comprehensive study [
17], combining chemical and mineralogical characterization [
59,
60]. Additionally, this approach could indicate the potential economic benefits of the generated residues [
61].
SEP is commonly employed to assess the possible movement of pollutants in different simulated environmental scenarios. However, it is crucial to carefully select and test reagents to prevent incomplete extraction and potential alterations in the outcomes of the experiment [62; 58]. In this sense, studies have been carried out on samples containing high levels of iron oxides and iron oxyhy-droxysulfates [
21], and latter improved by [
35] for phases like hydrobasaluminite when Al content is associated to Fe(III) minerals. This procedure delineates five distinct fractions: water-soluble, exchangeable, poorly crystallized oxyhydroxides, well-crystallized oxyhydroxides, and residual. The distinction among these fractions depends on their crystallinity, which plays a key role in determining the release of contaminants.
Based on this improved SEP, this study has evaluated the potential leaching of PTEs from AMDp in the different weathering scenarios (
Figure 6) by comparing with limit values established by the criterion continuous concentration (CCC) from the recommended water quality criteria adopted by the European Community [
63]. This value represents a potential negative effect for the aquatic life when uncontrolled leachates are generated from an inappropriate disposal of the AMDp. The proposed scenarios are Scenario 1: Contact of the residue with rainwater (emulated by F1 SEP fraction) which may occur in uncovered disposal facilities and are considered highly dangerous to the environment according to regulations [
64]; Scenario 2: Interaction between the residue and weak acidic leachates under reducing conditions (F1+F2 fraction). These conditions are relevant to scenarios such as underground disposal facilities; and finally, Scenario 3: Exposure of the residues to weak acidic leachates in oxidizing environments (F1+F3+F4 fraction). These conditions aim to simulate disposal scenarios in surface impoundments or waste piles, as previously has been described by [
12].
The data suggests that, although some samples may occasionally fall below the CCC limit, when the resulting whole-body residue from any DAS system is disposed of, it could pose danger to aquatic ecosystems and the life associated with them. In scenario 1 (exposure to rainfall), Pb, Zn, Al, and Fe exceed the proposed values for DAS-Ca, while Cd, Cr, Zn, and Fe behave similarly for DAS-Mg. For the rest of the more aggressive conditions, scenario 2 (underground disposal) and scenario 3 (surface impoundments or waste piles), the AMD residues generate leachates that surpass the threshold for most of the elements studied, making them incompatible with aquatic life. In this regard, these metal-rich residues from passive treatments (both laboratory scale and field scale), when stored in an oxidizing environment, demonstrated significant leaching of metals and sulfates [58; 65]. Several studies have shown that residues from Fe-rich AMD are unstable and have a high risk of leaching. They should be stored underwater or in a neutral pH environment to prevent contamination release and ensure safe disposal [5; 6]. Meanwhile, specific studies on MgO-DAS also recommend storing them in a dry (not mix acidic or alkaline waters) environment [
15]. More recent studies suggest that in addition to being stored under anoxic conditions in the presence of water to prevent the formation of new AMD, these residues should undergo stabilization/solidification before any storage [
8].
In addition,
Table 4 presents the regulatory limits for waste acceptance at landfills in the EU [
64] in comparison to the concentrations of elements detected in AMDp (DAS-Ca and DAS-Mg) following the F1 leaching fraction. This fraction denotes the metals and metalloids capable of mobilization upon contact with water and are considered as the leaching limit values to be applied to the waste for storage in the three types of landfills proposed by the Royal Decree 646/2020, under the auspices of the [
64].
According to the results obtained, the DAS-Ca residues would be classified as hazardous since the concentration of Zn (104 mg/kg) exceeds the reference value, while Cu (50 mg/kg) and Pb (10 mg/kg) are at the limit value. Meanwhile, the precipitates of the DAS-Mg would pose storage issues in landfills due to the excessive leached concentrations of Cd (901), Cu (164), and Zn (206 mg/kg). Under these circumstances, the EU regulation [
64] recommends additional treatment (such as further neutralization, encapsulation, etc.) before the storage of these materials in any type of regulated landfill. Also, recent directives established and emphasized, in general terms, the necessity of subjecting waste destined for landfills to appropriate pretreatment [
66], considering as an essential requirement to ensure that landfill operations are conducted without endangering human health and the environment. Additionally, in accordance with the circularity principle, it is important to note that these treatments generate metal-rich solid residues considered generally hazardous wastes, but that can also present an opportunity for the potential recovery of metals and other critical elements [58; 67-68]. In this regard, the European Regional Development Fund (FEDER, 2021-2027) aims to enhance these efforts by fostering the transition to a circular and resource-efficient economy (Raw Materials Initiative) through the generation of added value from waste. This initiative contributes to reducing the amount of waste sent to landfills and conserving natural resources by promoting reuse and recycling.