Wild Plants for the Phytostabilization of Phosphate Mine Waste in Semi-Arid Environments: A Field Experiment

: The management of mine waste has become an urgent issue, especially in semi-arid environments. In this context, and with an aim to inhibit the oxidation of the sulﬁde tailings of the abandoned mine of Kettara in Morocco, a store-and-release (SR) cover made of phosphate mine waste (PW) was implemented. In order to guarantee its long-term performance, phytostabilization by local wild plant species is currently the most effective and sustainable solution. This study aimed to assess the growth performance and phytostabilization efﬁciency of ﬁve local wild plant species to grow on the SR cover made of PW. A ﬁeld experiment was conducted for two growing seasons (2018 and 2019), without amendments and with the minimum of human care. PW and the aboveground and belowground parts of the studied plant species were collected and analyzed for As, Cd, Cu, Ni, and Zn. The bioconcentration factor (BCF) and translocation factor (TF) were also calculated. Despite the hostile conditions of the mining environment, the ﬁve plant species showed promising growth performances as follows: Atriplex semibaccata > Vicia sativa > Launaea arborescens > Peganum harmala > Asparagus horridus . The ﬁve plants showed high accumulation capacity of the trace elements, with the highest concentrations in belowground tissue. Principal component analysis distinguished A. semibaccata as having a high concentration of Cu and As, while Asparagus horridus had higher concentrations of Cd and Zn. In contrast, P. harmala , V. sativa , and L. arborescens demonstrated afﬁnity regarding Ni. According to the BCF (<1) and TF (<1), these plant species could be used as effective phytostabilizers of the studied trace elements. The present study showed that local wild plant species have a great potential for the phytostabilization of PW, and could ensure the long-term efﬁciency of SR cover.


Introduction
Morocco has a long-standing history and an important role in the mining trade due to the diversity of its mineral resources [1]. It exploits many minerals, mainly phosphate, lead, copper, gold, silver, and barite, among others [2]. Phosphate mining is considered a fundamental pillar of the Moroccan economy [3]. According to the U.S. Geological Survey, Mineral Commodity Summaries [2], Morocco is the world leader in terms of exports of phosphate, with 31% of the market and 50 billion tons in reserve, and the Following the preceding considerations, the question that naturally arises is the suitability or not of local wild plant species for the phytostabilization of the SR cover made of PW in a semi-arid environment. Therefore, this field experiment study aimed to evaluate and compare the performance and potential for phytostabilization of PW by five local wild plants in a semi-arid environment.

Field Experimental Design
The studied site is the abandoned mine in Kettara (470 masl, 31 • 51 36 N and 8 • 9 36 W) located 35 km north-northwest of the city of Marrakech (Figure 1a,b). The area is characterized by semi-arid bioclimatic conditions [12,36]. , tailing impoundment area (c), and the store-and-release cover design at the experimental park (d) (modified from Bossé et al. [12]).
The field experiment was performed based on a pre-designed 10 × 10 m experimental plot located in the Kettara abandoned mine experimental park. The plot was designed, constructed, and equipped in order to assess the hydrogeological behavior of the SR cover made of PW [6,12] (Figure 2). Details of substrate characteristics and different thicknesses used are given in Bossé et al. [12]. This experiment was designed in order to finalize the global plan of reclamation of the Kettara abandoned mine (Figure 1d). The experiment was conducted under local conditions without any amendments over 2 years from February 2018 to December 2019. Irrigation was only applied at the beginning of the experiment following the establishment of the seedlings, so as to better evaluate the plant performance with minimal intervention and care. Specimens of the same species were planted in agricultural soil 1 km away as controls and were treated and sampled in the same way.

Plant and Soil Sampling
In order to evaluate plant growth, two plant samplings were carried out at the end of each growing season. Shoots and leaves of each plant species were harvested, ovendried at 70 °C, and weighed to determine the average dry biomass. After the 2-year experiment, plant species were carefully collected, separated into aerial parts and roots, and stored in a plastic bag and frozen until further analysis.
PW samples were collected in the zone in contact with the belowground parts of plants. A soil sample was collected 1 km away from the control agricultural soil. Soil samples were mixed, homogenized, sieved to 2 mm, and oven-dried at 60 °C for 48 h, then stored in polyethylene bags until further analysis.

Soil and Plant Tissue Analyses
Physico-chemical soil properties were determined according to the standard procedures published in the Official Methods of Soil Analysis [37]. Fine particles were analyzed for the determination of electrical conductivity (1:5 w/v soil:water) and pH (1:2.5 w/v soil:water) using an LF 92 WTW multiparameter probe (Cond-pH 1970i; WTW GmbH, Weilheim, Germany).
Trace element concentrations were determined by adding 7 mL HNO3 (65% w/w), 2 mL of HF (40% w/w), and 1 mL of HClO4 (60% w/w) to 500 mg of the PW and control soil samples, previously calcinated in a muffle furnace at 500 °C. Then each sample solution was diluted to 50 mL. The samples were analyzed in triplicate for trace elements (As, Cd,

Plant Species
The five plant species used in this field experiment were selected according to their performance from preliminary greenhouse results [36], and to their botanical and ecological characteristics: namely, Vicia sativa L., Asparagus horridus L., Peganum harmala L., Launaea arborescens (Batt.) Murb., and Atriplex semibaccata R. Br. Fourteen plants of each species of similar sizes, growth (2 months), and general aspects were randomly transplanted into the plot in February 2018 ( Figure 2).
The experiment was conducted under local conditions without any amendments over 2 years from February 2018 to December 2019. Irrigation was only applied at the beginning of the experiment following the establishment of the seedlings, so as to better evaluate the plant performance with minimal intervention and care. Specimens of the same species were planted in agricultural soil 1 km away as controls and were treated and sampled in the same way.

Plant and Soil Sampling
In order to evaluate plant growth, two plant samplings were carried out at the end of each growing season. Shoots and leaves of each plant species were harvested, oven-dried at 70 • C, and weighed to determine the average dry biomass. After the 2-year experiment, plant species were carefully collected, separated into aerial parts and roots, and stored in a plastic bag and frozen until further analysis.
PW samples were collected in the zone in contact with the belowground parts of plants. A soil sample was collected 1 km away from the control agricultural soil. Soil samples were mixed, homogenized, sieved to 2 mm, and oven-dried at 60 • C for 48 h, then stored in polyethylene bags until further analysis.

Soil and Plant Tissue Analyses
Physico-chemical soil properties were determined according to the standard procedures published in the Official Methods of Soil Analysis [37]. Fine particles were analyzed for the determination of electrical conductivity (1:5 w/v soil:water) and pH (1:2.5 w/v soil:water) using an LF 92 WTW multiparameter probe (Cond-pH 1970i; WTW GmbH, Weilheim, Germany).
Trace element concentrations were determined by adding 7 mL HNO 3 (65% w/w), 2 mL of HF (40% w/w), and 1 mL of HClO 4 (60% w/w) to 500 mg of the PW and control soil samples, previously calcinated in a muffle furnace at 500 • C. Then each sample solution was diluted to 50 mL. The samples were analyzed in triplicate for trace elements (As, Cd, Cu, Ni, and Zn) using inductively coupled plasma-atomic emission spectrometry (ICP-AES) according to ISO 22,036 (2008) [38].
Plant samples were carefully washed with tap water, followed by deionized water, oven-dried at 80 • C for 72 h, and then ground into fine powder using a pestle and mortar. The concentrations of trace elements in plant tissues were determined according to the method of Temminghoff and Houba [39]. Plant tissue samples (200 mg) were fully mineralized with HNO 3 (65%), HF (40%), and H 2 O 2 (30%) in a microwave oven and the concentrations of trace elements (As, Cd, Cu, Ni, and Zn) were determined using ICP-AES (Optima 3100 RL; Perkin Elmer Waltham, MA, USA).

Phytoremediation Factors and Phytostabilization Efficiency
The phytostabilization efficiency of the selected plants and their accumulation capacity for trace elements in their belowground and aboveground fractions were assessed by the calculation of two key factors: the translocation factor (TF) and the bioconcentration factor (BCF).
The TF indicates a plant's efficiency in translocating accumulated trace elements from its belowground parts to its aboveground parts. It is calculated according to Ali et al. [40], as in Equation (1). The BCF indicates the efficiency of a plant species to accumulate trace elements in its tissue from the surrounding substrate and is calculated according to Ali et al. [40], using Equation (2

Statistical Analysis
Statistical differences in Zn, Cu, As, Ni, and Cd concentrations in different plant species and dry biomass means were assessed using one-way analysis of variance (ANOVA) followed by Tukey's post hoc test to determine the significant difference between years and plant species. All statistical analyses were performed using SPSS Statistics version 21 (IBM Corp., Armonk, NY, USA), with p-values ≤ 0.05 considered statistically significant. A principal component analysis (PCA) was performed using R 4.0.3 software (R Core Team, 2020) to distinguish the tolerance preference of plant species regarding the accumulation of trace elements.

Trace Element Content in Substrate
The pH, electrical conductivity (EC), and trace elements in PW and control samples are summarized in Table 1. The results showed that PW presented high concentra- tions of Cu (118.14 ± 33.71 mg kg −1 ) and Zn (96.09 ± 11.23 mg kg −1 ), followed by Ni (21.10 ± 2.48 mg kg −1 ) and Cd (11.13 ± 1.96 mg kg −1 ), whereas As presented the lowest value at 9.61 ± 1.22 mg kg −1 . Moreover, the results showed a significant difference (p < 0.05) in the concentrations of Cu, Zn, and Cd in the PW compared with those in the control soil. According to Kabata-Pendias and Pendias [41], Kabata-Pendias [42], and Kloke et al. [43], the threshold limit concentrations values (TLVs) range from 100 to 150 mg kg −1 for Cu, 300 mg kg −1 for Zn, 50 mg kg −1 for Ni, and 20 mg kg −1 for As. Therefore, the concentrations of these trace elements (Cu, Zn, Ni, and As) present in PW do not exceed the level of the environmental quality standards, except for the concentration of Cd, which exceeded the TLV (3 mg kg −1 ), which can be considered as a potential source of surrounding environment contamination [43]. These results are consistent with those of Hakkou et al. [11]; the same authors showed, using kinetic testing, that PW does not exhibit a significant generation of contaminants despite the presence of the abovementioned contaminants. The PW had a pH near neutrality (7.42 ± 0.4) and a relatively high EC (285.2 ± 9.4 µs/cm). The results of Zine et al. [36] showed similar values for the same substrates. Values are presented as mean ± SD (n = 3). Different letters in the same row indicate significant differences (p < 0.05).

Plant Biomass
The plant biomass values of the studied plants are presented in Figure 3. Comparisons of the biomass produced by each species revealed significant differences between plant species in both 2018 and 2019. The growth of the seedlings of Atriplex semibaccata R. Br. and Vicia sativa L., followed by Launaea arborescens (Batt.) and Peganum harmala L., Murb., demonstrated a high adaptation to the harsh climatic conditions and the plants did not seem be experiencing stress due to the hostile characteristics of the PW. In contrast, the plants of Asparagus horridus L. showed low biomass values, which reflects the species' relatively limited capacity to adapt to environmental stress, since according to Zine et al. [36], it showed a higher growth speed in a greenhouse experiment, in contrast to other species such as Launaea arborescens (Batt.) and Peganum harmala L., Murb.
Although we observed a contrasting performance between 2018 and 2019, notably for Atriplex semibaccata and Vicia sativa (Figure 3), these differences were not considered statistically significant (p < 0.05), probably due to the contrast in terms of precipitation between 2018 and 2019. In general, chamaephyte growth is known to be highly modified by the precipitation rate [44,45]. A. semibaccata reached impressive values after 2 years of cultivation, with an average annual productivity per plant of almost 135.35 ± 6.90 g m −2 . Pioneer plant species, such as Vicia sativa and A. semibaccata, have been reported to be characterized by a high tolerance to mining waste, with a rapid growth rate and excellent soil coverage [46,47]. These properties confer advantages to these kinds of plant species as promising candidates for the phytomanagement of mining waste. According to Conesa et al. [48], Moreno-Jiméenez et al. [49], and Yuanyuan et al. [50], in general, these local plant species are often genetically better equipped in terms of survival, growth, and the ability to reproduce under environmental stress than other exotic plant species.

Trace Element Concentrations in Plant Species
The concentrations of metal and metalloid trace elements in aboveground and belowground plant fractions were evaluated for the five plant species, and the corresponding results are presented in Figure 4. The results show that As, Cd, Cu, Ni, and Zn concentrations in aboveground plant parts (shoots) were very low compared with those in belowground tissues (roots). Concerning the accumulation of the studied trace elements, the results reveal a common observation for Vicia sativa L., Peganum harmala L., Launaea arborescens (Batt.) Murb., and Atriplex semibaccata R. Br, as the concentrations of the different trace elements in the aboveground and belowground parts of these five plant species all obeyed the following order: Zn > Cu > Ni > Cd > As. In contrast, the results for Asparagus horridus L. followed a different order: Zn > Cu > Cd > Ni > As.
For each metal considered, there was a significant difference (p < 0.05) between plant species, except for Ni, which shows similar values, especially in the aboveground parts of Vicia sativa L., Asparagus horridus L., and Peganum harmala L.
Trace element concentrations varied in the aboveground and belowground parts according to the species (Figure 4). Among the five species, Atriplex semibaccata R. Br. and Vicia sativa L. had the highest concentration of As in shoots at 4.27 ± 0.08 mg kg⁻ 1 dw and 2.14 ± 0.07 mg kg⁻ 1 dw, and 2.23 ± 0.08 mg kg⁻ 1 dw and 1.91 ± 0.31 mg kg⁻ 1 dw, in roots and shoots, respectively. Asparagus horridus L. had the highest concentration of Cd, with 7.82 ± 0.70 and 5.75 ± 0.75 mg kg⁻ 1 dw, followed by Vicia sativa L., with 3.59 ± 0.45 and 3.60 ± 0.63 mg kg⁻ 1 dw in roots and shoots, respectively. Concerning the concentration of Cu, Atriplex semibaccata R. Br. had the highest concentrations in roots and shoots, with values at 29.03 ± 2.01 and 17.39 ± 2.23 mg kg⁻ 1 dw in roots and shoots, respectively, followed by Peganum harmala L. with 14.11 ± 0.19 mg kg⁻ 1 dw in roots, and 11.80 ± 2.90 mg kg⁻ 1 dw in shoots.
The concentrations of As, Cd, and Cu in the roots and shoots of Atriplex semibaccata R. Br. and those of Asparagus horridus L. were two or three times higher in comparison to those measured in the roots and shoots of the other species. The Zn concentration was higher in Asparagus horridus L. and Atriplex semibaccata R. Br. Regarding the concentrations of Ni, the plant species can be arranged in two groups, each with similar results, the first

Trace Element Concentrations in Plant Species
The concentrations of metal and metalloid trace elements in aboveground and belowground plant fractions were evaluated for the five plant species, and the corresponding results are presented in Figure 4. The results show that As, Cd, Cu, Ni, and Zn concentrations in aboveground plant parts (shoots) were very low compared with those in belowground tissues (roots). Concerning the accumulation of the studied trace elements, the results reveal a common observation for Vicia sativa L., Peganum harmala L., Launaea arborescens (Batt.) Murb., and Atriplex semibaccata R. Br, as the concentrations of the different trace elements in the aboveground and belowground parts of these five plant species all obeyed the following order: Zn > Cu > Ni > Cd > As. In contrast, the results for Asparagus horridus L. followed a different order: Zn > Cu > Cd > Ni > As.
For each metal considered, there was a significant difference (p < 0.05) between plant species, except for Ni, which shows similar values, especially in the aboveground parts of Vicia sativa L., Asparagus horridus L., and Peganum harmala L.
Trace element concentrations varied in the aboveground and belowground parts according to the species (Figure 4). Among the five species, Atriplex semibaccata R. Br. and Vicia sativa L. had the highest concentration of As in shoots at 4.27 ± 0.08 mg kg −1 dw and 2.14 ± 0.07 mg kg −1 dw, and 2.23 ± 0.08 mg kg −1 dw and 1.91 ± 0.31 mg kg −1 dw, in roots and shoots, respectively. Asparagus horridus L. had the highest concentration of Cd, with 7.82 ± 0.70 and 5.75 ± 0.75 mg kg −1 dw, followed by Vicia sativa L., with 3.59 ± 0.45 and 3.60 ± 0.63 mg kg −1 dw in roots and shoots, respectively. Concerning the concentration of Cu, Atriplex semibaccata R. Br. had the highest concentrations in roots and shoots, with values at 29.03 ± 2.01 and 17.39 ± 2.23 mg kg −1 dw in roots and shoots, respectively, followed by Peganum harmala L. with 14.11 ± 0.19 mg kg −1 dw in roots, and 11.80 ± 2.90 mg kg −1 dw in shoots. tion threat to animal health [14]. However, according to the National Research Council [31], the concentration of Cd, ranging from 3 to 100 mg kg⁻ 1 , found in phosphate sedimentary rocks, could be a source of contamination of the trophic chain, even if its mobility and bioavailability to plants is highest in acid soil than in alkaline one [31]. Therefore, the aboveground biomass of these plants can be valorized and used in the production of nonfood products [51,52].

Bioconcentration and Translocation Factors
The BCF and TF values of the five local plant species for the studied trace elements are presented in Table 2. The results showed that all BCF values were less than one for the five studied trace elements (As, Cd, Cu, Ni, and Zn). Similarly, all plants showed a TF less than 1. However, TF values were higher than BCF values for the same trace element and for the same plant species. In addition, the values of BCF reported in Table 2 are remarkably low and indicate that the concentrations of trace elements found in the aboveground The concentrations of As, Cd, and Cu in the roots and shoots of Atriplex semibaccata R. Br. and those of Asparagus horridus L. were two or three times higher in comparison to those measured in the roots and shoots of the other species. The Zn concentration was higher in Asparagus horridus L. and Atriplex semibaccata R. Br. Regarding the concentrations of Ni, the plant species can be arranged in two groups, each with similar results, the first comprising Vicia sativa L., Asparagus horridus L., and Peganum harmala L., and the second comprising Launaea arborescens (Batt.) Murb. and Atriplex semibaccata R. Br.
All concentrations of trace elements in the aboveground parts of the studied plants were below the domestic animal toxicity limits for cattle grazing according to the Committee on Minerals and Toxic Substances in Diets and Water for Animals, National Research Council [31] (As < 100 mg kg −1 , Cd < 10 mg kg −1 , Cu < 500 mg kg −1 , Ni < 30 mg kg −1 , and Zn < 500 mg kg −1 ). Thus, the harvested parts of the plants do not present a contamination threat to animal health [14]. However, according to the National Research Council [31], the concentration of Cd, ranging from 3 to 100 mg kg −1 , found in phosphate sedimentary rocks, could be a source of contamination of the trophic chain, even if its mobility and bioavailability to plants is highest in acid soil than in alkaline one [31]. Therefore, the aboveground biomass of these plants can be valorized and used in the production of non-food products [51,52].

Bioconcentration and Translocation Factors
The BCF and TF values of the five local plant species for the studied trace elements are presented in Table 2. The results showed that all BCF values were less than one for the five studied trace elements (As, Cd, Cu, Ni, and Zn). Similarly, all plants showed a TF less than 1. However, TF values were higher than BCF values for the same trace element and for the same plant species. In addition, the values of BCF reported in Table 2 are remarkably low and indicate that the concentrations of trace elements found in the aboveground parts of the five local plant species were much lower than the concentrations found in the PW.
Statistical analysis demonstrated different responses for the five plant species regarding trace elements. Lee et al. [53] reported that BCF and TF values are particularly dependent on the plant species used, on the degree of substrate contamination, its chemical properties, and on the time of year. This is due to the fact that the absorption and reaction of trace elements differs from one element to another, as well as from one species to another [54].
The assessment and selection of plants for phytomanagement for phytoremediation purposes depend on BCF and TF values [40,55]. According to Nouri et al. [56] and Ali et al. [40], both these factors are crucial in screening plants for the phytostabilization of trace elements in mining waste. Indeed, plant species with a BCF < 1 and a TF < 1 are suitable for use in phytoremediation through phytostabilization, while plant species with both BCF and TF values greater than one are suitable for use in phytoextraction [57].
The BCF and TF values show good agreement with a previous study conducted under greenhouse conditions [36], wherein the same plant species gave similar results. Therefore, our results cast new light on and confirm the capacity of the selected local plant species for the management of PW through their phytostabilization in semi-arid environments. According to Mendez and Maier [19], the selection of suitable plant species candidates is a crucial issue, and requires specific attention. The plant species destined to be used in the phytomanagement of mining waste by phytostabilization should produce significant biomass, have a dense root system, and the translocation of trace elements from belowground to aboveground parts should be as minimal as possible; they should also be well adapted to local environmental stress. The selection of local plant species for the phytostabilization of SR made of PW requires another specific condition regarding the length of the roots of the plants. The length of the roots should not exceed 1 m, since this is the thickness of SR required to ensure that the oxidation of the sulfide tailings of Kettara is prevented [11,12]. Therefore, this research exploited the use of chamaephytes known for their limited root system [58]. Values are expressed as mean ± SD (n = 3). Different letters in the same column designate significant differences (p < 0.05).

Candidate Local Plant Species for Phytostabilization
The biplot of the principal component analysis (PCA) carried out on the concentrations of trace elements in the aboveground parts (shoots) of the five plant species (Figure 5a) accounted for 72.9% of the total variance. Dim1 (43.6% of the total inertia) represents the capacity of the local plant species to accumulate diverse trace elements synchronously. It distinguished Atriplex semibaccata with high concentrations of Cu in its aboveground parts, and it separated Asparagus horridus with significant concentrations of Cd and Zn in its aboveground parts from the other plant species, namely, Vicia sativa, Peganum harmala, and Launaea arborescens, which contained lower concentrations of these trace elements. Dim2 (29.3% of the total inertia) represents the concentration gradients of As, Cd, Cu, Ni, and Zn ( Figure 5a). This dimension highlights a relevant affinity of Atriplex semibaccata for Cu, and of Asparagus horridus for Cd and Zn, which are separated from the other plant species that are only weakly correlated with Cu, Cd, and Zn. It distinguished Atriplex semibaccata with high concentrations of Cu in its aboveground parts, and it separated Asparagus horridus with significant concentrations of Cd and Zn in its aboveground parts from the other plant species, namely, Vicia sativa, Peganum harmala, and Launaea arborescens, which contained lower concentrations of these trace elements. Dim2 (29.3% of the total inertia) represents the concentration gradients of As, Cd, Cu, Ni, and Zn (Figure 5a). This dimension highlights a relevant affinity of Atriplex semibaccata for Cu, and of Asparagus horridus for Cd and Zn, which are separated from the other plant species that are only weakly correlated with Cu, Cd, and Zn. The biplot of the PCA carried out on trace elements in belowground parts (roots) of the plant species (Figure 5b) explains 85.8% of the total variance. Dim1 (49.1% of the total inertia) explains the main information concerning the plant species' capability to accumulate different trace elements simultaneously (Figure 5b). Atriplex semibaccata presented a high concentration of Cu and As, and Aspragaus horridus had higher concentrations of Cd and Zn; this indicated that the belowground parts of these plant species have accumulative capacity for these trace elements. The other plant species, namely Peganum harmala, Vicia sativa, and Launaea arborescens, had a preference for accumulating Ni.  The biplot of the PCA carried out on trace elements in belowground parts (roots) of the plant species (Figure 5b) explains 85.8% of the total variance. Dim1 (49.1% of the total inertia) explains the main information concerning the plant species' capability to accumulate different trace elements simultaneously (Figure 5b). Atriplex semibaccata presented a high concentration of Cu and As, and Aspragaus horridus had higher concentrations of Cd and Zn; this indicated that the belowground parts of these plant species have accumulative capacity for these trace elements. The other plant species, namely Peganum harmala, Vicia sativa, and Launaea arborescens, had a preference for accumulating Ni.
According to Krzesłowska [59] and Gajic et al. [60], the uptake of trace elements and their flow and translocation depends on their concentrations, availability in the substrate, local climatic conditions, and their affinity to plant species, as well as the plant species itself. Indeed, the reaction of metallic and metalloid trace elements differs from one element to another as well as from one plant species to another [54]. The intrinsic properties of each plant species, such as the adaptive mechanisms involved in the uptake and the translocation of trace elements, are considered the decisive factors of the affinity of plant species for trace elements [60]. Indeed, plant species differ from one another in the genetic expression of transporter genes controlling the absorption and efflux of trace elements through their tissues [61,62]. Therefore, the use of local wild plant species is highly recommended given that they are often genetically better equipped in terms of growth and the ability to reproduce under the prevailing environmental conditions [49,50].
The plant species used in this experiment are predominantly perennial, namely, Asparagus horridus L., Atriplex semibaccata R. Br., Launaea arborescens (Batt.) Murb., and Peganum harmala L. According to Grime [63], these plant species have a long life cycle with high vegetative growth, which characterizes their K-selection strategy. Vicia sativa L. belongs to the legume plant family, which is considered a keystone of any plant community because of its ability for the fixation of atmospheric nitrogen by using bacteria (Rhizobium Frank 1889) [64]. This promotes the natural processes of plant establishment, increases their rate of growth, and improves soil recovery [60]. Asparagus horridus L., Atriplex semibaccata R. Br., and Launaea arborescens (Batt.) Murb. are known as plant species of semi-arid environments, and they are not very demanding in terms of water and nutrients; therefore, they also offer an interesting choice for the phytomanagement of a mining site in this environment [40]. These plant species are also important in the colonization of waste disposal sites because they adopt R-selection as a strategy of propagation by producing a large number of seeds [50,60]. According to Skousen et al. [65], these kinds of plant species can be competitive and their development delays the appearance of woody species which represent a threat to the long-term performance of SR cover made of PW.

Conclusions
The present field experiment demonstrated the high capacity of the five tested plant species, namely, Asparagus horridus L., Atriplex semibaccata R. Br., Launaea arborescens (Batt.) Murb., Peganum harmala L., and Vicia sativa L., to accumulate several trace elements, e.g., As, Cd, Cu, Ni, and Zn, with the highest concentrations in the belowground parts compared to those in aboveground part tissues. The concentrations of the studied element traces in the plant were below the domestic animal toxicity limits values, with careful consideration of the potential contamination of the trophic chain by cadmium. However, regarding the BCF and TF values, these local plants can be considered excellent candidates for the phytostabilization of PW. Atriplex semibaccata R. Br. presents a high affinity for Cu and As, Aspragaus horridus L. for Cd and Zn, whilst Peganum harmala L., Vicia sativa L., and Launaea arborescens (Batt.) Murb. possess a good tolerance to Ni. Overall, we conclude that the best candidates for establishing a plant community for the phytostabilization of PW and ensure the long-term performance of the SR cover are Atriplex semibaccata R. Br. and Vicia sativa L. as keystone species of the community, accompanied by Launaea arborescens (Batt.) Murb. and Peganum harmala L.