# Screening of Native Plants Growing on a Pb/Zn Mining Area in Eastern Morocco: Perspectives for Phytoremediation

Said El Hasnaoui, Mouna Fahr, Catherine Keller, Clément Levard, Bernard Angeletti, Perrine Chaurand, Zine El Abidine Triqui, Abdelkarim Guedira, Laila Rhazi, Fabrice Colin, et al.

► To cite this version:

Said El Hasnaoui, Mouna Fahr, Catherine Keller, Clément Levard, Bernard Angeletti, et al. "Screening of Native Plants Growing on a Pb/Zn Mining Area in Eastern Morocco: Perspectives for Phytoremediation. Plants, 2020, 9 (11), pp.1-23. 10.3390/plants9111458. hal-02985696

HAL Id: hal-02985696

https://amu.hal.science/hal-02985696

Submitted on 2 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Screening of Native Plants Growing on a Pb/Zn Mining Area in Eastern Morocco: Perspectives for Phytoremediation

Said El Hasnaoui<sup>1,2,●</sup>, Mouna Fahr<sup>1,2,*</sup>●, Catherine Keller<sup>2,3</sup>●, Clément Levard<sup>2,3</sup>, Bernard Angeletti<sup>2,3</sup>, Perrine Chaurand<sup>2,3</sup>●, Zine El Abidine Triqui<sup>1,2</sup>, Abdelkarim Guedira<sup>1,2</sup>, Laila Rhazi<sup>4</sup>, Fabrice Colin<sup>2,3</sup> and Abdelaziz Smouni<sup>1,2,*</sup>●

Received: 1 September 2020; Accepted: 13 October 2020; Published: 29 October 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Screening of native plant species from mining sites can lead to identify suitable plants for phytoremediation approaches. In this study, we assayed heavy metals tolerance and accumulation in native and dominant plants growing on abandoned Pb/Zn mining site in eastern Morocco. Soil samples and native plants were collected and analyzed for As, Cd, Cu, Ni, Sb, Pb, and Zn concentrations. Bioconcentration factor (BCF), translocation factor (TF), and biological accumulation coefficient (BAC) were determined for each element. Our results showed that soils present low organic matter content combined with high levels of heavy metals especially Pb and Zn due to past extraction activities. Native and dominant plants sampled in these areas were classified into 14 species and eight families. Principal components analysis separated *Artemisia herba-alba* with high concentrations of As, Cd, Cu, Ni, and Pb in shoots from other species. Four plant species, namely, *Reseda alba*, *Cistus libanotis*, *Stipa tenacissima*, and *Artemisia herba-alba* showed strong capacity to tolerate and hyperaccumulate heavy metals, especially Pb, in their tissues. According to BCF, TF, and BAC, these plant species could be used as effective plants for Pb phytoextraction. *Stipa tenacissima* and *Artemisia herba-alba* are better suited for phytostabilization of Cd/Cu and Cu/Zn, respectively. Our study shows that several spontaneous and native plants growing on Pb/Zn contaminated sites have a good potential for developing heavy metals phytoremediation strategies.

**Keywords:** phytoextraction; phytostabilization; native plants; heavy metal; hyperaccumulation; tolerance; Pb/Zn mining area

## 1. Introduction

Soil contamination by heavy metals has become one of the major environmental problems around the world [1]. Anthropogenic activities such as mining, smelting, and industrial processing are the most important sources of heavy metals entering the environment [2–4]. Heavy metals persist a long time in the soils as these elements are not biodegradable [2,5]. This contamination affects the taxonomic,

<footnote n="1">1. Laboratoire de Biotechnologie et Physiologie Végétales, Centre de Biotechnologie Végétale et Microbienne Biodiversité et Environnement, Faculté des Sciences, Université Mohammed V de Rabat, 10000 Rabat, Morocco; said.elhasnaoui@um5s.net.ma (S.E.H.); z.triqui@um5s.net.ma (Z.E.A.T.); a.guedira@um5s.net.ma (A.G.)</footnote>

<footnote n="2">2. Laboratoire Mixte International Activité Minière Responsable "LMI-AMIR", IRD/UM5/INAU, 10000 Rabat, Morocco; keller@cerege.fr (C.K.); levard@cerege.fr (C.L.); angeletti@cerege.fr (B.A.); chaurand@cerege.fr (P.C.); fabrice.colin@ird.fr (F.C.)</footnote>

<footnote n="3">3. CNRS, IRD, INRAE, Collège de France, CEREGE, Aix Marseille University, 13100 Aix-en-Provence, France</footnote>

<footnote n="4">4. Centre de Biotechnologie Végétale et Microbienne Biodiversité et Environnement, Faculté des Sciences, Université Mohammed V de Rabat, 10000 Rabat, Morocco; laila.rhazi@um5.ac.ma</footnote>

<footnote>* Correspondence: mouna.fahr@um5.ac.ma or m.fahr@um5s.net.ma (M.F.); abdelaziz.smouni@um5.ac.ma (A.S.)</footnote>
functional diversity, and microbial properties of soils [6]. Moreover, due to their high concentrations in soils, metals can be transferred and accumulated not only in the plants [7], but also in the humans being, causing many serious health disorders [8]. There are several technologies to remediate heavy metals contaminated soils (e.g., excavation of contaminated material and chemical/physical treatment). However, many of these technologies are costly or do not achieve a long-term nor aesthetic solution [9,10].

Phytoremediation is a so-called "green" and "cost-effective" technique [4] whereby plant species are used to detoxify contaminated soils through immobilization or hyper/accumulation processes [11,12]. The remediation of toxic metals by plant can be divided into two most important strategies: (1) phytostabilization: plants are used to stabilize contaminated soil from the potential migration of pollutants by immobilizing contaminants in root system through absorption into the roots or precipitation in the rhizosphere [13]; and (2) phytoextraction: plants are used to remove metal pollutants from soils via root absorption and their concentration in the shoots and leaves [4,14]. Phytoextraction uses both accumulator and "hyperaccumulator" plants. These latter can accumulate unusual amounts of metals (>1% and up to 10%) in their shoots [14,15]. However, plants adaptation to environmental conditions of target regions is among factors limiting the success of the phytoremediation strategy [11]. Spontaneous and native plant species could overcome these constraints, constituting a key element in the efficiency of phytoremediation strategies. The native plants are often efficient in terms of growth, reproduction, and survival under environmental stress [16]. In this way, there is a continuing interest in searching for native plants that are tolerant to heavy metals and can promote phytoremediation [12].

The screening of spontaneous and native plant species from mining sites all over the world becomes a subject of interest. In Morocco, there are many natural metalliferous sites and metal-contaminated soils, and their metallophyte flora still poorly explored. Oued el Heimer and Touissite sites in eastern Morocco are among the most demonstrative examples of abandoned mining areas. After overexploitation during more than 70 years [17], they are a source of contamination by heavy metals because of poorly controlled waste disposal. These mining areas show an adapted native flora that has the potential to be used for different phytoremediation strategies to rehabilitate these contaminated areas. The objectives of this study are (i) screening of native plants growing on a contaminated area in the vicinity of a Pb/Zn mining area in eastern Morocco; and (ii) identification of tolerant and/or hyperaccumulator plants to suggest suitable plants for use in phytoremediation.

## 2. Materials and Methods

### 2.1. *Sample Collection*

Soil and plant sampling were collected in the "Oued el Heimer (34°26′50.5″ N 1°54′07.3″ W) and Touissite (34°28′07.7″ N 1°46′20.9″ W)" areas during May 2018. Eighty-four samples of dominant plants were collected, (about 6 species from Oued el Heimer and 8 from Touissite with 6 replicas per species). To obtain the accurate plant species information, plants were identified according to flora of Morocco [18–20]. Soil samples, from the plants rooting zone (0–20 cm depth), were also collected. Seventy-eight samples of soils were collected (about 6 samples from Oued el Heimer and 7 from Touissite with 6 replicas per sample) (Figures 1 and 2).

### 2.2. Soil and Plant Analysis

Soil samples were air-dried and sieved to <2 mm powder before use. The percentage of organic matter (OM) content were determined according to the method described by Combs et al. [21]. Soil pH was determined in 5:1 water/soil. The mixture was stirred for 1 h with a magnetic stirrer and then decanted for 30 min. The pH was determined on the supernatant directly using a pH meter "Bante210 Benchtop pH/mV Meter, Bante instruments CO., LTD. Shanghai, China" [17]. Electrical conductivity (EC) (1:5 w/v soil water suspensions) was carried out using a multi-parameter measuring instrument
"MultiLine® Multi 3510 IDS, WTW a xylem brand". For heavy metal content analysis, soil samples were dried at 72 °C to constant weight for two days and digested with a mixture following this procedure: 100–300 mg of soil material is placed in Polytetrafluoroethylene (PTFE) tubes. The soil samples are first moistened with 0.2 mL of water and then in order, 2 mL of concentrated HCl (37%, A Fisher Scientific International Company Pittsburgh, Pennsylvania, USA) and 1 mL of concentrated HNO₃ (69%, Loba Chemie PVT. LTD Colaba, Mumbai, INDIA) are added. The mixture is incubated for 16 h at room temperature. The mixture is then brought to the boil for 2 h at 150 °C by placing a glass bead above each tube for reflux. After dissolving the soil sample, the mixture is transferred quantitatively to a 15 mL tube, adjusting with distilled water [22]. Plant samples were separated into shoots and roots, washed thoroughly with running tap water to remove adhering substrate materials, rinsed twice for 15 s in cold 0.2 mM CaSO₄, and rinsed with cold distilled water. The plant samples were dried at 72 °C for 48 h and grounded to powder. Then, plant samples were digested according to the acid hydrolysis protocol described by Temminghoff and Houba [22]. The concentrations of As, Cd, Cu, Ni, Sb, Pb, and Zn in digested soil or plant samples were then analyzed in ICP-AES (inductively coupled plasma-atomic emission spectrometry; Ultima2 JY, Horiba Company, Edison, NJ, USA) according to Margui et al. [23] with the blank values being deducted from the measurement.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Plant and soil sampling area at Oued el Heimer site.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Plant and soil sampling area at Touissite site.</figcaption></figure>
## 2.3. Soil Pollution Index (PI)

The pollution index of soils is calculated by dividing the metal concentrations in soil (mg kg⁻¹) by the background metal concentration according to Chon et al. [24] and Wu et al. [25]. This background concentration corresponds to standard levels in uncontaminated soils [26]. Thus, the pollution index is calculated using the following equation:

$$PI = ([As]_{soil}/6 + [Cd]_{soil}/3 + [Cu]_{soil}/100 + [Ni]_{soil}/100 + [Pb]_{soil}/100 + [Zn]_{soil}/300 + [Sb]_{soil}/8.6)/7$$
$$[metal]_{soil}: concentration in mg kg⁻¹ \qquad (1)$$

where a PI greater than 1 corresponds to a polluted soil.

## 2.4. Soil Enrichment Factor (EF)

The enrichment factor (EF) is used to assess the impact of human activities on the natural environment. It is calculated by the following equation [27]:

$$EF = (C_i/C_{ref})_{sample}/(C_i/C_{ref})_{background} \qquad (2)$$

where (Cᵢ/Cᵢ<sub>ref</sub>) sample is the ratio of the content of the interested element and reference element in the soil sample, and (Cᵢ/Cᵢ<sub>ref</sub>) background is the ratio of the content of the interested element and reference element in natural background (Average Crust) [26] (Table S1). In this study, Scandium (Sc) was selected as a reference element [28,29].

According to Liu et al. [30] five contamination categories are generally recognized based on the enrichment factor: EF < 1, no enrichment; 1 ≤ EF < 10, minimal enrichment; 10 ≤ EF < 100, moderate enrichment; 100 ≤ EF < 1000, significant enrichment; and EF > 1000, extremely high enrichment.

## 2.5. Bioconcentration Factor, Translocation Factor, and Biological Accumulation Coefficient

To evaluate the metal enrichment characteristics of different plant species, three factors were calculated:

The biological concentration factor (BCF) is found by dividing the plant root metal concentration by the metal concentration in rhizospheric soils (Table S1) [31]. This factor shows the accumulation ability for root from the soil. It is calculated using the following formula:

$$BCF = [Metal] root/[Metal] soil \qquad (3)$$

Translocation factor (TF) is the metal ratio transfer capacity from root to shoots (Table S1) [31]. It is calculated using the following formula:

$$TF = [Metal] shoot/[Metal] root \qquad (4)$$

Biological accumulation coefficient (BAC) is defined as the concentration of metals in plant shoots divided by metal concentration in rhizospheric soils (Table S1) [32]. It is calculated using the following formula:

$$BAC = [Metal] shoot/[Metal] soil \qquad (5)$$

## 2.6. Data Analysis

Differences in metal concentrations in roots and shoots among species were studied using one-way ANOVA. This analysis was made separately for each metal in the shoots and the roots. To identify the significant pair-wise differences between species for each metal concentration, post hoc Tukey HSD tests were used.

To study the relations between metals concentrations in shoots and roots for the studied plant species, principal component analysis (PCA) was performed. The data matrix contains the average
values (6 replicates) of each metal in the shoots and in the roots of each species. The similarity between species in the ability of metal accumulation in shoots and roots was studied using ascending hierarchical classification (AHC).

Multivariate and univariate analyses were carried out using STATISTICA 10 (StatSoft Inc., Tulsa, OK, USA).

# 3. Results and Discussions

## 3.1. Site Description

The studied area “Oued el Heimer—Touissite” is located in the North-West of Morocco, at the mountainous band (Horst chain), which stretches over a length of 100 km and a width of 20–30 km. This chain is North limited by the corridor Oujda-Taourirt, on the West by the plain of Tafrata, in the South by the vast domain of high plateaus, and in the East by the Morocco–Algerian border [17]. Semi-arid Mediterranean climate conditions with a dry and hot summer of variable duration and wet winter prevail in this area. Annual rainfall rates range from 90 to 100 mm/year [33].

Oued el Heimer foundry is close to the village of Oued el Heimer, which is located approximately 33 km southeast of Oujda town. This smelter, closed in 2013, specialized in the treatment, fusion, and refining of Pb and silver (Ag). Oued el Heimer area is subjected to intense pollution by a massive accumulation of hazardous wastes in the neighborhood of the smelter. These wastes form slag particle mounds at the banks of the Oued el Heimer creek. In the period of floods, waters carry important quantities of wastes. This natural leaching constitutes a significant source of contamination of grounds and water sources (Figure 1) [17].

The Touissite mining site, located thirty kilometers south the Oujda city, was closed in 2002, following 75 years' period of exploitation. Source of pollution in these areas originates mainly from mining waste which forms sandy dams on large surfaces. These sands are easily transported by winds and rainfall. Some of these dams were covered by residual materials directly stemming from the ore extracting process (Figure 2) [17].

The Pb-Zn ores of the Touissite area substitute recrystallized Liassic dolomites at the base of a block-faulted Jurassic section, which rests on Paleozoic basement rocks. Galena, sphalerite with pyrite and marcasite dominate the mineral assemblages with lesser chalcopyrite, bornite, and tetrahedrite. Pb-Sb sulphosalts (bournonite and stephanite) and Ag-bearing species (argentite and native silver) are rather scarce [34–36] (Figure 3).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Geology of the Touissite area (Ministry of energy and mines geology department, 1985).</figcaption></figure>
## 3.2. Soil Properties and Metals Concentrations

The soils were slightly acidic to neutral with a pH ranging from 6.45 to 7.51 at the Oued el Heimer site and from 6.91 to 7.76 at Touissite (Table 1). These results confirm the data reported by Smouni et al. [17] related to soils of the same area of Oued el Heimer. Midhat et al. [8] report pH values of 7.4 in soils of Sidi Bou-Othmane mining site located in Marrakech (Morocco). Navarro et al. [37], however, found pH values ranging from 3.1 to 8.4 in the Pb-Zn mined areas of Cabezo Rajao in Spain. Many authors have shown that soil pH plays a crucial role in the solubility and bioavailability of metals [38]. Bliefert and Perraud [39] reported that solubility of Pb decreases with increasing pH. Compared to pH, EC values showed greater variability ranging from 0.40 to 2.28 and from 0.23 to 0.38 mS/cm respectively at Oued el Heimer and Touissite. These values indicate low salinities of both studied site soils [40]. Soil samples from Oued el Heimer presented a low organic matter (OM) content that varies respectively between 0.41% and 2.80%. The same trend was observed for soil samples from the Touissite site with 0.28% to 2.07% for organic matter (Table 1). Low OM contents reflect low biological activities and plant developments. These may be due to specific environmental conditions combined to high amount of heavy metals in soils. Indeed, some authors suggest that high levels of Zn in soils, negatively influence the activity of micro-organisms and earthworms, thus inhibiting the breakdown of organic matter [41].

Heavy metals concentrations in soils sampled from Oued el Heimer and Touissite sites are presented in Table 1. The soils from Touissite site have As concentrations varying between 43 and 82.90 mg kg⁻¹ dry weight (DW), between 15 and 36 mg kg⁻¹ DW for Cd, 328 and 1405 mg kg⁻¹ DW for Cu, 9 and 16 mg kg⁻¹ DW of Ni, 6445 and 18,324 mg kg⁻¹ DW of Pb, 96.71 and 242.80 mg kg⁻¹ DW of Sb, and between 2096 and 5385 mg kg⁻¹ for Zn. The area of Oued el Heimer, where the smelter is located, has a content of As ranging between 18.7 and 466 mg kg⁻¹ DW, Cd from 32 to 280 mg kg⁻¹ DW, Cu from 35 to 592 mg kg⁻¹ DW, Ni between 12 and 44.4 mg kg⁻¹ DW, Pb between 611 and 12,461 mg kg⁻¹ DW, Sb from 33.90 to 247.10 mg kg⁻¹ DW, and Zn between 318 and 43,540 mg kg⁻¹ DW. This area is characterized by high Zn and Pb levels. Therefore, soils from these areas can be considered as highly contaminated [26]. Pollution index (PI) recorded at both sites ranged from 3.62 to 67.90 at Oued el Heimer and from 14.00 to 38.51 at Touissite (Table 1). This expresses the extremely deleterious effect of these mining activities and abandoned wastes on the environment.

Enrichment factor (EF) defines the concentrations of a given element at a specific location as compared to the average natural occurrences [42]. In Oued el Heimer site, the highest EF values for As, Cd, Cu, Ni, Pb, Zn, and Sb were 364.8, 3656.4, 98.8, 8.2, 2995, 2749.6, and 2393.3, respectively (Table 1). At Touissite mine, the highest enrichment factor (EF) for As, Cd, Cu, Ni, Pb, Zn, and Sb were respectively 52.5, 410.4, 166.5, 2.3, 2758.3, 1092.2, and 1877.8, (Table 1). According to EF categories [30], all metals in both areas were affected by human activities. Among them, the EF values of Ni at both sites are less than 10 (1 ≤ EF < 10), which was moderately enriched and affected both by a natural source and human activities [30]. The enrichment factor values for Cu (Oued el Heimer site) and As (Touissite site) belonging to category (10 ≤ EF < 100), indicate moderate enrichment. The EF values for As (Oued el Heimer site), Cd, and Cu (Touissite site) belonging to category (100 ≤ EF <1000), indicate that As, Cd, and Cu were significantly enriched.

The EF values for Cd at Oued el Heimer site and Pb, Zn, and Sb at both sites attached to category (EF > 1000), indicate extremely high enrichment, and these values revealed that heavy metals contents are affected mainly by anthropogenic activities (Human factors) [30].
<table><caption>**Table 1.** Soil properties (pH, EC: Electrical conductivity, OM: Organic Matter) and range concentrations (mg kg⁻¹ DW) of As, Cd, Cu, Ni, Pb, Sb, and Zn in soils sampled from studied areas.</caption><thead><tr><th>Site</th><th></th><th>As</th><th>Cd</th><th>Cu</th><th>Ni</th><th>Pb</th><th>Zn</th><th>Sb</th><th>PI</th><th>pH</th><th>EC (mS/cm)</th><th>% OM</th></tr></thead><tbody><tr><td rowspan="2">Oued el Heimer</td><td>Concentration range in soils (mg kg<sup>-1</sup>)</td><td>18.7–466</td><td>32–280</td><td>35–592</td><td>12–44.4</td><td>611–12,461</td><td>318–43,540</td><td>33.9–247.1</td><td>3.62–67.9</td><td>6.4–7.51</td><td>0.40–2.28</td><td>0.41–80%</td></tr><tr><td>Enrichment factor (EF)</td><td>15–364.8</td><td>233–3656</td><td>5.3–98.8</td><td>2.5–8.2</td><td>131.7–2995</td><td>26.5–2749</td><td>285.7–2393</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Touissite</td><td>Concentration range in soils (mg kg<sup>-1</sup>)</td><td>43–82.9</td><td>15–36</td><td>328–1405</td><td>9–16</td><td>6445–18,324</td><td>2096–5387</td><td>96.7–242.8</td><td>14.0–38.5</td><td>6.9–7.7</td><td>0.23–0.38</td><td>0.28–2.07%</td></tr><tr><td>Enrichment factor (EF)</td><td>22.2–52.5</td><td>133–410.4</td><td>29.5–166.5</td><td>1.7–2.3</td><td>833.5–2758</td><td>115.2–1092</td><td>550.3–1877</td><td></td><td></td><td></td><td></td></tr></tbody></table>
Many reports recorded different degrees of soil metal concentrations in the vicinity of mining sites and smelters. For example, Navarro et al. reported 19 and 53.1 mg kg⁻¹ DW of Pb and Zn, respectively, in Cabezo Rajao abandoned mine (Pb–Zn) in SE Spain [37]. El Hachimi et al. [43] reported the same trend of heavy metal contents in some Pb and Zn mining site in High Moulouya from Morocco, 5547 mg Pb kg⁻¹ and 7500 mg Zn kg⁻¹ at Zaida mine, 10,520 and 9075 mg kg⁻¹ of Pb and Zn, respectively, at Mibladen mine, and 2101 mg kg⁻¹ of Pb and 3125 mg kg⁻¹ of Zn at Aouli. Furthermore, Yassir et al. [40] reported 321.7 and 723.4 mg kg⁻¹ of Pb and Zn, respectively, at Sidi Bou Othmane abandoned mine (Marrakech, Morocco). The concentrations of Pb and Zn in the Ahangaran mine at Malayer City, South of Hamedan in Iran, are respectively about 8955.4 and 12,963.1 mg kg⁻¹ DW [44]. The Pb and Zn concentrations found in Oued el Heimer and Touissite were much lower than those reported in two mining districts (Maline and Les Avinières) located in Southern France (more than 84,130 mg Pb kg⁻¹ and 91,454 mg Zn kg⁻¹) by Escarré et al. [45]. Other studies reported approximately the same range of As, Cd, and Cu recorded in soils from our study [37,40]. For Antimony (Sb), many reports recorded 100.6–5045 mg kg⁻¹ DW [46], 527–11,798 mg kg⁻¹ DW [47], and between 74.2 and 16,389 mg kg⁻¹ of Sb [48] in different mining sites in China. Moreover, about 2900 mg kg⁻¹ of Sb were measured at Sb smelting sites in Japan [49]. According to previous studies, the majority of abandoned mines present high levels of soil metals contamination due to anthropogenic activities. We report similar results. In this context, identification of native plants and characterization of their heavy metals tolerance and accumulation potential may offer efficient solutions for contaminated area restoration.

### 3.3. Screening of Native Plant Species in Oued el Heimer and Touissite Areas

A total of 14 plant species were identified and classified to 8 families according to the nomenclature of "Flore Pratique du Maroc" [18-20] (Table 2) with a dominance of the Brassicaceae, Fabaceae, and Asteraceae families. At Oued el Heimer site, six samples were collected downstream of the el Heimer Oued in the area located 500 m far from the smelter and surrounded by pin forest and identified as *Cistus libanotis*, *Artemisia herba-alba*, and *Capsella bursa-pastoris* belonging to Cistaceae, Asteraceae, and Brassicaceae families, respectively. Other samples taken directly from the slag were identified as *Hirschfeldia incana*, *Stipa tenacissima* and *Agathophora alopecuroides* belonging respectively to Brassicaceae, Poaceae, and Amaranthaceae families. In the Touissite site covered by mining waste, eight samples were identified belonging to eight species: *Reseda alba*, *Convolvulus althaeoides*, *Hedysarum spinosissimum*, *Phragmites communis*, *Lotus corniculatus*, *Capsella bursa-pastoris*, *Scolymus hispanicus* L., and *Rapistrum rigosum*. The sampled plant species were divided into two life-forms, perennial and annual according to the flora of Morocco [18-20]. All plant species belong to the most important groups of worldwide spontaneously distributed in the Mediterranean region [50-57]. Many authors reported their occurrence in both uncontaminated and contaminated soils, including mining sites ecosystems. For example, several species of Brassicaceae, Cistaceae, Fabaceae, and Resedaceae families grow naturally in Pb-Zn and Sb contaminated soils in Spain [58,59]. Furthermore, Asteraceae, Poaceae, and Amaranthaceae are found in the vicinity of different mining sites in southern central Morocco [8]. The majority of these families could also be found at other mining sites located outside of the Mediterranean region such as in 18 metalliferous areas of the Austrian Alps [60], in the lead-zinc mine in Iran [44,61,62], in the gold mining site, in Colombia [63], in Zn rich mining soils in China [41], in Pb-Zn mining areas in Vietnam [64], in a metal-contaminated site in Florida [16] and in As-contaminated soils in the Czech Republic [65].

### 3.4. Heavy Metals Concentration in Plants

Metal roots and shoots contents are given in Table 3 that also provide phytotoxic levels and hyperaccumulation threshold issued from Kabata-Pendias and pendias [66]; Kramer [67] and Bioconcentration factors (BCF), translocation factors (TF), and biological accumulation coefficient (BAC) in Table 4.
<table><caption>**Table 2.** List of sampled plants in the vicinity of the studied sites.</caption><thead><tr><th>Site</th><th>Plants</th><th>Family</th><th>Life Span</th></tr></thead><tbody><tr><td rowspan="8">Touissite</td><td>Reseda alba</td><td>Resedaceae</td><td>Annual</td></tr><tr><td>Convolvulus althaeoides</td><td>Convolvulaceae</td><td>Perennial</td></tr><tr><td>Hedysarum spinosissimum</td><td>Fabaceae</td><td>Annual</td></tr><tr><td>Phragmites communis</td><td>Poaceae</td><td>Perennial</td></tr><tr><td>Lotus corniculatus</td><td>Fabaceae</td><td>Perennial</td></tr><tr><td>Capsella bursa-pastoris</td><td>Brassicaceae</td><td>Annual</td></tr><tr><td>Scolymus hispanicus</td><td>Asteraceae</td><td>Perennial</td></tr><tr><td>Rapistrum rigosum</td><td>Brassicaceae</td><td>Annual</td></tr><tr><td rowspan="5">Oued el Heimer</td><td>Cistus libanotis</td><td>Cistaceae</td><td>Perennial</td></tr><tr><td>Agathophora alopecuroides</td><td>Amaranthaceae</td><td>Perennial</td></tr><tr><td>Hirschfeldia incana</td><td>Brassicaceae</td><td>Perennial</td></tr><tr><td>Stipa tenacissima</td><td>Poaceae</td><td>Perennial</td></tr><tr><td>Artemisia herba-alba</td><td>Asteraceae</td><td>Perennial</td></tr><tr><td></td><td>Capsella bursa-pastoris</td><td>Brassicaceae</td><td>Annual</td></tr></tbody></table>

### 3.4.1. Arsenic (As)

Arsenic (As) total concentrations in the roots and shoots was significantly different between species (Tables 3 and 4). As contents in the roots of *Stipa tenacissima* and *Artemisia herba-alba* were significantly (*p* < 0.05) higher as compared to those in the other plant species. While *Rapistrum rigosum* accumulate the low concentration of this metal (1.3 mg kg⁻¹ DW). In shoots, As amount in shoots of *Artemisia herba-alba* followed by *Stipa tenacissima* and *Cistus libanotis* were significantly higher compared to other species. *Artemisia herba-alba* present significantly the highest value up to 50.7 mg kg⁻¹ DW while the lowest accumulation of As was found in *Rapistrum rigosum* (0.3 mg kg⁻¹ DW). Two plant species from the Oued el Heimer site reach the maximum values of As in their roots and shoots, including *Stipa tenacissima* and *Artemisia herba-alba*. However, all plant species collected around the investigated sites accumulated As concentrations less than 80 mg kg⁻¹ in their roots and shoots, which is a phytotoxic As concentration in plants [67]. Arsenic hyperaccumulators are defined as plants that accumulate more than 1000 mg kg⁻¹ DW in their shoots [68]. In this study, none of the 14 plant species showed Arsenic (As) concentrations in their tissues (Table 3) that exceed these thresholds; thus, none of them is As hyperaccumulators.

Several authors used BCF, TF, and BAC to evaluate the metal accumulation efficiency in plants and to estimate the plants potential for phytostabilization and/or phytoextraction [69,70]. BCF, TF, and BAC with values higher than 1 are used as criteria to select the plants used for phytoextraction, while BCF > 1 and TF < 1 had been used to evaluate the potential of plants for phytostabilization [71]. None of the collected plant species has BCF > 1 (Table 5). However, nine of them have a TF slightly >1 such as *C. libanotis* (2.87), *R. alba* (2.37), *C. althaeoides* (2.12), *C. bursa-pastoris* (1.53) from Touissite and *S. hispanicus* L. (1.32), *P. communis* (1.31), *C. bursa-pastoris* from Oued el Heimer (1.13), *H. incana* (1.10), and 1.04 in *A. herba-alba*. All plant species have a BCF and BAC < 1 (Table 5) indicating that none of the fourteen plant species has the potential to be used in As phytostabilization (BCF > 1 and TF < 1) or for As phytoextraction (BCF, TF, and BAC > 1).
<table><caption>**Table 3.** Pair-wise differences of heavy metals concentrations (mg kg⁻¹ DW) in the roots/shoots of plants collected from the studied sites in the oriental region of Morocco.</caption><thead><tr><th rowspan="2">Site</th><th rowspan="2">Plants</th><th colspan="2">As</th><th colspan="2">Cd</th><th colspan="2">Cu</th><th colspan="2">Ni</th><th colspan="2">Pb</th><th colspan="2">Zn</th><th colspan="2">Sb</th></tr><tr><th>Root</th><th>Shoot</th><th>Root</th><th>Shoot</th><th>Root</th><th>Shoot</th><th>Root</th><th>Shoot</th><th>Root</th><th>Shoot</th><th>Root</th><th>Shoot</th><th>Root</th><th>Shoot</th></tr></thead><tbody><tr><td rowspan="10">Touissite</td><td><i>R. alba</i></td><td>(2.6) <sup>a</sup></td><td>(6.3) <sup>bc</sup></td><td>(2.1) <sup>a</sup></td><td>(4.1) <sup>ab</sup></td><td>(45.1) <sup>ab</sup></td><td>(46.1) <sup>de</sup></td><td>(2.9) <sup>abc</sup></td><td>(2.0) <sup>ab</sup></td><td>(322.7) <sup>a</sup></td><td>(1607.5) <sup>f</sup></td><td>(254.3) <sup>a</sup></td><td>(465.5) <sup>c</sup></td><td>(6.1) <sup>ab</sup></td><td>(21.7) <sup>abcd</sup></td></tr><tr><td><i>C. althaeoides</i></td><td>(2.4) <sup>a</sup></td><td>(5.1) <sup>abc</sup></td><td>(2.9) <sup>a</sup></td><td>(1.3) <sup>a</sup></td><td>(59.7) <sup>ab</sup></td><td>(35.4) <sup>bcd</sup></td><td>(1.89) <sup>ab</sup></td><td>(2.2) <sup>ab</sup></td><td>(714.3) <sup>abc</sup></td><td>(426.3) <sup>abcd</sup></td><td>(286.1) <sup>a</sup></td><td>(207.4) <sup>ab</sup></td><td>(16.0) <sup>bcd</sup></td><td>(13.6) <sup>abc</sup></td></tr><tr><td><i>H. spinosissimum</i></td><td>(6.0) <sup>a</sup></td><td>(3.6) <sup>ab</sup></td><td>(8.7) <sup>a</sup></td><td>(3.9) <sup>ab</sup></td><td>(58.5) <sup>ab</sup></td><td>(32.6) <sup>bcd</sup></td><td>(4.1) <sup>abc</sup></td><td>(2.3) <sup>ab</sup></td><td>(983.6) <sup>abc</sup></td><td>(253.9) <sup>ab</sup></td><td>(347.0) <sup>ab</sup></td><td>(314.0) <sup>abc</sup></td><td>(35.1) <sup>e</sup></td><td>(33.1) <sup>cde</sup></td></tr><tr><td><i>P. communis</i></td><td>(3.8) <sup>a</sup></td><td>(5.0) <sup>abc</sup></td><td>(2.7) <sup>a</sup></td><td>(2.0) <sup>ab</sup></td><td>(66.7) <sup>ab</sup></td><td>(65.8) <sup>ef</sup></td><td>(2.2) <sup>ab</sup></td><td>(2.6) <sup>ab</sup></td><td>(2305.7) <sup>d</sup></td><td>(720.4) <sup>bcd</sup></td><td>(198.9) <sup>a</sup></td><td>(432.5) <sup>c</sup></td><td>(22.1) <sup>cd</sup></td><td>(25.5) <sup>abcde</sup></td></tr><tr><td><i>L. corniculatus</i></td><td>(5.3) <sup>a</sup></td><td>(3.0) <sup>ab</sup></td><td>(3.0) <sup>a</sup></td><td>(1.8) <sup>a</sup></td><td>(117.3) <sup>bc</sup></td><td>(39.6) <sup>bcd</sup></td><td>(4.1) <sup>abc</sup></td><td>(3.3) <sup>ab</sup></td><td>(1493.0) <sup>bcd</sup></td><td>(832.4) <sup>bcd</sup></td><td>(276.6) <sup>a</sup></td><td>(170.5) <sup>a</sup></td><td>(27.0) <sup>de</sup></td><td>(14.8) <sup>abcd</sup></td></tr><tr><td><i>C. bursa-pastoris</i></td><td>(3.8) <sup>a</sup></td><td>(5.9) <sup>bc</sup></td><td>(4.4) <sup>a</sup></td><td>(4.2) <sup>ab</sup></td><td>(65.7) <sup>ab</sup></td><td>(49.0) <sup>def</sup></td><td>(4.1) <sup>abc</sup></td><td>(4.4) <sup>b</sup></td><td>(849.6) <sup>abc</sup></td><td>(733.8) <sup>bcd</sup></td><td>(587.0) <sup>b</sup></td><td>(666.9) <sup>d</sup></td><td>(27.0) <sup>de</sup></td><td>(55.5) <sup>e</sup></td></tr><tr><td><i>S. hispanicus L.</i></td><td>(3.6) <sup>a</sup></td><td>(4.7) <sup>abc</sup></td><td>(5.2) <sup>a</sup></td><td>(5.4) <sup>ab</sup></td><td>(48.8) <sup>ab</sup></td><td>(44.3) <sup>cde</sup></td><td>(3.6) <sup>abc</sup></td><td>(4.2) <sup>b</sup></td><td>(798.7) <sup>abc</sup></td><td>(972.7) <sup>de</sup></td><td>(307.9) <sup>a</sup></td><td>(452.0) <sup>c</sup></td><td>(12.3) <sup>bc</sup></td><td>(17.2) <sup>abcd</sup></td></tr><tr><td><i>R. rigosum</i></td><td>(1.3) <sup>a</sup></td><td>(0.3) <sup>a</sup></td><td>(1.7) <sup>a</sup></td><td>(1.8) <sup>a</sup></td><td>(13.6) <sup>a</sup></td><td>(10.0) <sup>a</sup></td><td>(1.3) <sup>a</sup></td><td>(0.5) <sup>a</sup></td><td>(454.8) <sup>ab</sup></td><td>(47.1) <sup>a</sup></td><td>(239.9) <sup>a</sup></td><td>(326.8) <sup>abc</sup></td><td>(1.6) <sup>a</sup></td><td>(0.5) <sup>a</sup></td></tr><tr><td rowspan="6">Oued el Heimer</td><td><i>C. libanotis</i></td><td>(3.4) <sup>a</sup></td><td>(9.7) <sup>c</sup></td><td>(22.7) <sup>a</sup></td><td>(16.2) <sup>bc</sup></td><td>(23.6) <sup>ab</sup></td><td>(25.5) <sup>abcd</sup></td><td>(2.4) <sup>ab</sup></td><td>(3.1) <sup>ab</sup></td><td>(1219.2) <sup>abcd</sup></td><td>(1261.8) <sup>ef</sup></td><td>(135.6) <sup>a</sup></td><td>(161.7) <sup>a</sup></td><td>(8.3) <sup>ab</sup></td><td>(9.4) <sup>abc</sup></td></tr><tr><td><i>A. alopecuroides</i></td><td>(9.1) <sup>a</sup></td><td>(4.6) <sup>abc</sup></td><td>(4.9) <sup>a</sup></td><td>(18.0) <sup>bcd</sup></td><td>(26.0) <sup>ab</sup></td><td>(16.3) <sup>ab</sup></td><td>(4.0) <sup>abc</sup></td><td>(4.2) <sup>abc</sup></td><td>(235.1) <sup>a</sup></td><td>(293.8) <sup>ab</sup></td><td>(142.5) <sup>a</sup></td><td>(168.2) <sup>a</sup></td><td>(5.2) <sup>ab</sup></td><td>(6.2) <sup>abc</sup></td></tr><tr><td><i>H. incana</i></td><td>(6.0) <sup>a</sup></td><td>(6.6) <sup>bc</sup></td><td>(25.8) <sup>a</sup></td><td>(32.9) <sup>d</sup></td><td>(37.3) <sup>ab</sup></td><td>(19.5) <sup>ab</sup></td><td>(2.9) <sup>abc</sup></td><td>(2.3) <sup>ab</sup></td><td>(441.3) <sup>ab</sup></td><td>(343.7) <sup>abc</sup></td><td>(273.4) <sup>a</sup></td><td>(309.7) <sup>abc</sup></td><td>(5.3) <sup>ab</sup></td><td>(6.4) <sup>ab</sup></td></tr><tr><td><i>S. tenacissima</i></td><td>(59.3) <sup>b</sup></td><td>(19.6) <sup>d</sup></td><td>(241.2) <sup>b</sup></td><td>(28.9) <sup>cd</sup></td><td>(237.9) <sup>dc</sup></td><td>(29.7) <sup>abcd</sup></td><td>(6.4) <sup>bc</sup></td><td>(3.8) <sup>ab</sup></td><td>(3785.7) <sup>e</sup></td><td>(1146.3) <sup>ef</sup></td><td>(637.0) <sup>b</sup></td><td>(322.9) <sup>abc</sup></td><td>(167.7) <sup>f</sup></td><td>(32.1) <sup>bcd</sup></td></tr><tr><td><i>A. herba-alba</i></td><td>(48.7) <sup>b</sup></td><td>(50.7) <sup>e</sup></td><td>(50.2) <sup>a</sup></td><td>(56.2) <sup>e</sup></td><td>(203.1) <sup>cd</sup></td><td>(72.3) <sup>f</sup></td><td>(6.5) <sup>c</sup></td><td>(8.1) <sup>c</sup></td><td>(1748.7) <sup>cd</sup></td><td>(4672.2) <sup>g</sup></td><td>(577.0) <sup>b</sup></td><td>(357.0) <sup>bc</sup></td><td>(13.0) <sup>bc</sup></td><td>(41.6) <sup>de</sup></td></tr><tr><td><i>C. bursa-pastoris</i></td><td>(4.5) <sup>a</sup></td><td>(5.1) <sup>abc</sup></td><td>(41.4) <sup>a</sup></td><td>(30.0) <sup>cd</sup></td><td>(28.7) <sup>ab</sup></td><td>(18.2) <sup>abc</sup></td><td>(1.7) <sup>ab</sup></td><td>(2.1) <sup>ab</sup></td><td>(629.7) <sup>abc</sup></td><td>(920.9) <sup>cde</sup></td><td>(395.3) <sup>ab</sup></td><td>(396.1) <sup>bc</sup></td><td>(8.7) <sup>ab</sup></td><td>(12.9) <sup>abcd</sup></td></tr><tr><td colspan="2">Phytotoxic concentrations of metals [66,67]</td><td>2–80</td><td>0.1–3</td><td>20–30</td><td>10–50</td><td>0.6–28</td><td>100–300</td><td>5–10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Hyperaccumulation threshold [67]</td><td>&gt;1000</td><td>&gt;100</td><td>&gt;1000</td><td>&gt;1000</td><td>&gt;1000</td><td>&gt;10,000</td><td>&gt;1000</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Explanation:** The values in the table represent the average of heavy metals contents in roots & shoots for each plant species. The averages followed by the same letters do not differ significantly.
<table><caption>**Table 4.** Results of variance analysis (One-way ANOVA) comparing the difference in metal concentration among species (dF = 13).</caption><thead><tr><th rowspan="2"></th><th colspan="2">Shoots</th><th colspan="2">Roots</th></tr><tr><th>F</th><th>p</th><th>F</th><th>p</th></tr></thead><tbody><tr><td>As</td><td>111.03</td><td>&lt;0.0001</td><td>46.19</td><td>&lt;0.0001</td></tr><tr><td>Cd</td><td>29.74</td><td>&lt;0.0001</td><td>19.95</td><td>&lt;0.0001</td></tr><tr><td>Cu</td><td>14.82</td><td>&lt;0.0001</td><td>10.53</td><td>&lt;0.0001</td></tr><tr><td>Ni</td><td>6.91</td><td>&lt;0.0001</td><td>3.06</td><td>0.0180</td></tr><tr><td>Pb</td><td>101.65</td><td>&lt;0.0001</td><td>14.69</td><td>&lt;0.0001</td></tr><tr><td>Zn</td><td>18.43</td><td>&lt;0.0001</td><td>8.16</td><td>&lt;0.0001</td></tr><tr><td>Sb</td><td>7.54</td><td>&lt;0.0001</td><td>250.93</td><td>&lt;0.0001</td></tr></tbody></table>

### 3.4.2. Cadmium (Cd)

The Cd content in shoots and roots were significantly different between species (Tables 3 and 4). Cd content in roots of *Stipa tenacissima* up to 241.2 mg kg⁻¹ DW was significantly much higher than those observed in other plants (Tables 3 and 4). While the lowest accumulation of Cd was recorded in *R. rigosum* (1.7 mg kg⁻¹ DW) in roots. Midhat et al. [8] indicate that the plants from *Stipa* genus, such as *Stipa capensis Thunb* accumulates 7.82 mg kg⁻¹ DW in their root, which means that *Stipa tenacissima* is characterized by a strong capacity to take up Cd from the soil compared to other species from the same genus. In shoots, the highest concentrations of Cd are found significantly (Table 4) in *Artemisia herba-alba* followed by *H. incana*, *C. bursa-pastoris*, *S. tenacissima*, *A. alopecuroides*, and *C. libanotis* harvested from Oued el Heimer site compared to those in the other plant species and reach respectively 56.1, 32.9, 30.0, 28.9, 18.0, and 16.2 mg kg⁻¹ DW (Table 3). Many reports recorded various amounts of Cd in shoot parts such as 0.62 mg kg⁻¹ DW in *Artemisia vulgaris* growing near a mining site in Northern Vietnam [72], 3.28 and 4 mg kg⁻¹ DW in *Stipa capensis Thunb* and *Stipa barbata Desf.*, respectively, growing around lead-zinc mine [8,73]. Plants generally exhibit measurable Cd amounts, especially in roots, but also in shoots resulting from uptake and translocation [8]. According to Kramer [67], *R. alba*, *H. spinosissimum*, *C. bursa-pastoris* (Touissite), *S. hispanicus* L., *C. libanotis*, *A. alopecuroides*, *H. incana*, *S. tenacissima*, *A. herba-alba*, and *C. bursa-pastoris* (Oued el Heimer) accumulate more than the Cd phytotoxic concentrations (0.1 to 3 mg kg⁻¹) in plant and are considered Cd tolerant species. A shoots concentration higher than 100 mg kg⁻¹ is used as a criterion for Cd hyperaccumulation [68]. On this basis, none of the 14 plant species can be considered as Cd hyper accumulator.

The BCF for Cd varied from 0.07 in *P. communis* to 2.72 in *S. tenacissima* (Table 5). Five plant species have a TF > 1, namely *A. alopecuroides* (3.67), *R. alba* (1.91), *H. incana* (1.27), *A. herba-alba* (1.11), and *S. hispanicus* L. (1.03) (Table 5). The values of the translocation factor are higher than the bioconcentration factor. This indicates that plants translocate Cd from roots to shoots [16]. The value of biological accumulation coefficient (BAC) higher than one is recorded only in *A. herba-alba* (1.10) (Table 5). Based BCF, TF, and BAC values, none of the plant species has the potential to be used in Cd phytoextraction [14]. *S. tenacissima* has BCF = 2.72 and TF = 0.11, also showed a high accumulation of Cd in its roots, which suggests that it is a suitable plant for Cd phytostabilization. This native plant species with both the capacity to accumulate high amounts of Cd in its roots and high biomass production could be used to minimize the migration of Cd in soil. This process reduces metal mobility and leaching towards groundwater's [74].
<table><caption>**Table 5.** Bioconcentration factor (BCF), Translocation factor (TF) and Biological accumulation coefficient (BAC) for the studied plants species.</caption><thead><tr><th rowspan="2">Site</th><th rowspan="2">Plants</th><th colspan="15">Trace Metal Elements</th></tr><tr><th colspan="2">As</th><th colspan="2">Cd</th><th colspan="2">Cu</th><th colspan="2">Ni</th><th colspan="2">Pb</th><th colspan="2">Zn</th><th colspan="2">Sb</th></tr><tr><th></th><th></th><th>BCF</th><th>TF</th><th>BAC</th><th>BCF</th><th>TF</th><th>BAC</th><th>BCF</th><th>TF</th><th>BAC</th><th>BCF</th><th>TF</th><th>BAC</th><th>BCF</th><th>TF</th><th>BAC</th></tr></thead><tbody><tr><td rowspan="10">Touissite</td><td><i>R. alba</i></td><td>0.04</td><td><b>2.37</b></td><td>0.10</td><td>0.08</td><td><b>1.91</b></td><td>0.15</td><td>0.11</td><td><b>1.02</b></td><td>0.11</td><td>0.32</td><td>0.66</td><td>0.21</td><td>0.02</td><td><b>4.98</b></td><td>0.12</td><td>0.04</td><td><b>1.83</b></td><td>0.08</td><td>0.04</td><td><b>3.55</b></td><td>0.15</td></tr><tr><td><i>C. althaeoides</i></td><td>0.03</td><td><b>2.12</b></td><td>0.07</td><td>0.08</td><td>0.62</td><td>0.04</td><td>0.12</td><td>0.59</td><td>0.07</td><td>0.13</td><td><b>1.17</b></td><td>0.15</td><td>0.09</td><td>0.59</td><td>0.05</td><td>0.05</td><td>0.72</td><td>0.03</td><td>0.10</td><td>0.85</td><td>0.08</td></tr><tr><td><i>H. spinosissimum</i></td><td>0.13</td><td>0.59</td><td>0.08</td><td>0.57</td><td>0.44</td><td>0.24</td><td>0.17</td><td>0.55</td><td>0.09</td><td>0.25</td><td>0.57</td><td>0.14</td><td>0.15</td><td>0.25</td><td>0.03</td><td>0.10</td><td>0.90</td><td>0.10</td><td>0.36</td><td>0.94</td><td>0.34</td></tr><tr><td><i>P. communis</i></td><td>0.06</td><td><b>1.31</b></td><td>0.08</td><td>0.07</td><td>0.73</td><td>0.05</td><td>0.04</td><td>0.99</td><td>0.04</td><td>0.18</td><td><b>1.17</b></td><td>0.20</td><td>0.12</td><td>0.31</td><td>0.03</td><td>0.04</td><td><b>2.17</b></td><td>0.08</td><td>0.09</td><td><b>1.15</b></td><td>0.10</td></tr><tr><td><i>L. corniculatus</i></td><td>0.08</td><td>0.57</td><td>0.04</td><td>0.19</td><td>0.55</td><td>0.1</td><td>0.14</td><td>0.33</td><td>0.05</td><td>0.25</td><td>0.80</td><td>0.2</td><td>0.08</td><td>0.55</td><td>0.05</td><td>0.13</td><td>0.60</td><td>0.08</td><td>0.12</td><td>0.54</td><td>0.06</td></tr><tr><td><i>C. bursa-pastoris</i></td><td>0.04</td><td><b>1.53</b></td><td>0.07</td><td>0.12</td><td>0.95</td><td>0.11</td><td>0.09</td><td>0.74</td><td>0.06</td><td>0.27</td><td><b>1.08</b></td><td>0.29</td><td>0.09</td><td>0.88</td><td>0.08</td><td>0.13</td><td><b>1.13</b></td><td>0.15</td><td>0.11</td><td><b>2.05</b></td><td>0.22</td></tr><tr><td><i>S. hispanicus</i> L.</td><td>0.06</td><td><b>1.32</b></td><td>0.07</td><td>0.20</td><td><b>1.03</b></td><td>0.21</td><td>0.14</td><td>0.90</td><td>0.12</td><td>0.27</td><td><b>1.15</b></td><td>0.32</td><td>0.10</td><td><b>1.21</b></td><td>0.12</td><td>0.1</td><td><b>1.46</b></td><td>0.14</td><td>0.07</td><td><b>1.39</b></td><td>0.1</td></tr><tr><td><i>R. rigosum</i></td><td>-</td><td>0.2</td><td>-</td><td>-</td><td>0.8</td><td>-</td><td>-</td><td>0.7</td><td>-</td><td>-</td><td>0.4</td><td>-</td><td>-</td><td>0.1</td><td>-</td><td>-</td><td><b>1.36</b></td><td>-</td><td>-</td><td>0.3</td><td>-</td></tr><tr><td><i>C. libanotis</i></td><td>0.18</td><td><b>2.87</b></td><td>0.51</td><td>0.60</td><td>0.71</td><td>0.42</td><td>0.67</td><td><b>1.08</b></td><td>0.71</td><td>0.17</td><td><b>1.27</b></td><td>0.22</td><td><b>1.99</b></td><td><b>1.03</b></td><td><b>2.06</b></td><td>0.42</td><td><b>1.19</b></td><td>0.50</td><td>0.24</td><td><b>1.13</b></td><td>0.27</td></tr><tr><td><i>A. alopecuroides</i></td><td>0.48</td><td>0.50</td><td>0.24</td><td>0.15</td><td><b>3.67</b></td><td>0.56</td><td>0.74</td><td>0.62</td><td>0.05</td><td>0.23</td><td><b>1.05</b></td><td>0.24</td><td>0.38</td><td><b>1.24</b></td><td>0.48</td><td>0.07</td><td><b>1.18</b></td><td>0.08</td><td>0.03</td><td><b>1.19</b></td><td>0.04</td></tr><tr><td><i>H. incana</i></td><td>0.01</td><td><b>1.10</b></td><td>0.01</td><td>0.09</td><td><b>1.27</b></td><td>0.11</td><td>0.06</td><td>0.52</td><td>0.03</td><td>0.06</td><td>0.79</td><td>0.05</td><td>0.03</td><td>0.77</td><td>0.03</td><td>0.01</td><td><b>1.13</b></td><td>0.01</td><td>0.02</td><td><b>1.2</b></td><td>0.03</td></tr><tr><td><i>S. tenacissima</i></td><td>0.14</td><td>0.32</td><td>0.04</td><td><b>2.72</b></td><td>0.11</td><td>0.32</td><td><b>1.55</b></td><td>0.12</td><td>0.19</td><td>0.40</td><td>0.59</td><td>0.23</td><td>0.51</td><td>0.30</td><td>0.09</td><td>0.01</td><td>0.50</td><td>0.00</td><td>0.89</td><td>0.19</td><td>0.17</td></tr><tr><td><i>A. herba-alba</i></td><td>0.86</td><td><b>1.04</b></td><td>0.90</td><td>0.80</td><td><b>1.11</b></td><td><b>1.10</b></td><td><b>3.63</b></td><td>0.35</td><td><b>1.29</b></td><td>0.54</td><td><b>1.24</b></td><td>0.62</td><td>0.48</td><td><b>2.67</b></td><td><b>1.30</b></td><td><b>1.69</b></td><td>0.61</td><td>1.04</td><td>0.38</td><td><b>3.2</b></td><td><b>1.22</b></td></tr><tr><td><i>C. bursa-pastoris</i></td><td>0.01</td><td><b>1.13</b></td><td>0.01</td><td>0.86</td><td>0.72</td><td>0.16</td><td>0.05</td><td>0.63</td><td>0.03</td><td>0.07</td><td><b>1.24</b></td><td>0.09</td><td>0.05</td><td><b>1.46</b></td><td>0.07</td><td>0.03</td><td><b>1.00</b></td><td>0.03</td><td>0.03</td><td><b>1.48</b></td><td>0.05</td></tr></tbody></table>

**Bold values** indicate BCF, TF and BAC values greater than 1.0.
### 3.4.3. Copper (Cu)

Among the collected plants, Cu content in plants range from 10 mg kg⁻¹ DW (*R. rigosum* shoot) to 237.8 mg kg⁻¹ DW (*S. tenacissima* root). Cu contents in roots of *S. tenacissima* and *A. herba-alba* were significantly higher as compared to those in the other species, whereas the concentrations in shoots were found to be significantly (*p* < 0.05) higher in *A. herba-alba*, *P. communis*, and *C. bursa-pastoris* (from Touissit area) compared to other plant species (Tables 3 and 4). Other studies reported 23.42 mg kg⁻¹ DW in *Scolymus hispanicus* L. and 0.07 mg kg⁻¹ DW in *Convolvulus althaeoides* [8]. This variation of accumulation of Cu depends in part upon the soils physicochemical proprieties [75]. Nine plant species, namely *R. alba*, *C. althaeoides*, *H. spinosissimum*, *P. communis*, *L. corniculatus*, *C. bursa-pastoris* (Touissite), *S. hispanicus* L., *S. tenacissima*, and *A. herba-alba*, accumulate Cu above phytotoxic concentrations (20–30 mg kg⁻¹) in the plant tissues, which suggested that they are tolerant to Cu [67]. Copper hyperaccumulator plants are defined as plants that accumulate more than 1000 mg Cu/kg DW in aboveground tissues [76]. On this basis, none of the collected plant species can be considered as Cu hyperaccumulators.

The lowest BCF value was observed in *P. communis* (0.04) and the maximum BCF values were measured in the tissues of *A. herba-alba* (3.63) and *S. tenacissima* (1.55) (Table 5). The TF ranges from 0.12 in *S. tenacissima* to 1.08 in *C. libanotis*. For the biological accumulation coefficient, only *A. herba-alba* presented a BAC equal to 1.29. According to these values, none of the 14 plant species has the potential for Cu phytoextraction from the soil. However, the BCF was >1 and TF < 1 in both *A. herba-alba* and *S. tenacissima*. These plant species showed a high ability to tolerate and accumulate Cu in roots and are then suggested suitable for Cu phytostabilization.

### 3.4.4. Nickel (Ni)

In the present study, all plant species showed a low Ni content in their roots ranged between 1.3 mg kg⁻¹ DW in *R. rigosum* to 6.5 in *A. herba-alba*. The amount observed in *A. herba-alba* was significantly higher compared to the amount found in roots of other plant species (Tables 3 and 4). For the shoot parts, the lowest accumulation of Ni up to 0.5 mg kg⁻¹ DW is found in *R. rigosum* while the highest value up to 8.1 mg kg⁻¹ DW was reached significantly in *A. herba-alba* followed by *C. bursa-pastoris* (4.4 mg kg⁻¹ DW) and *S. hispanicus* (4.2 mg kg⁻¹ DW) compared to other species (Tables 3 and 4). According to Krämer [67], none of the plants sampled from both sites exceed phytotoxic Ni concentrations. This result from low concentration of Ni in soils. According to Reeves [77], none of plants species showed a Ni hyperaccumulation characteristic.

None of the plants harvested from Oued el Heimer and Touissite has the BCF > 1 (Table 5). However, eight plant species have translocation factors superior to 1 namely *C. libanotis* (1.27), *A. herba-alba* (1.24), *C. bursa-pastoris* (Oued el Heimer site) (1.24), *P. communis* (1.17), *C. althaeoides* (1.17), *S. hispanicus* L. (1.15), *C. bursa-pastoris* (Touissite site) (1.08), and *A. alopecuroides* (1.05) (Table 5). These plants have a good capacity to translocate Ni from roots to the shoots [14]. BAC values range from 0.05 in *H. incana* to 0.62 in *A. herba-alba*. According to Yoon et al. [16], none of these plant species has the potential to be used for phytostabilization (BCF > 1 and TF < 1) or phytoextraction (BCF, TF, and BAC > 1) of Ni.

### 3.4.5. Lead (Pb)

Pb concentrations in plant roots ranged from 235.1 to 3785.7 mg kg⁻¹ DW in *A. alopecuroides* and *S. tenacissima*, respectively (Table 3). This level of Pb in *S. tenacissima* was significantly higher (*p* < 0.05) to other plant species, followed by *P. communis*, and *A. herba-alba*. No considerable difference of Pb content was observed in the other plant species. Moreover, in eight plant species, the Pb amounts in roots were higher than those found in the shoots, indicating a low translocation rate from roots to shoots. Many studies reported Pb concentrations in different plant species growing from contaminated soils. Midhat et al. [8] and Brunetti et al. [78] recorded about 0.7 and 282.33 mg kg⁻¹ DW of Pb in
roots respectively of *Stipa austroitalica* and *Stipa capensis* Thunb. However, Nouri et al. [61] reported up to 1743 mg kg⁻¹ DW of Pb in the root of *Reseda alba* L growing in Pb-Zn mine in Iran. These values are higher than those observed in the root of *R.a alba* (322.7 mg kg⁻¹ DW) collected from Touissite site and lower than those found in *S. tenacissima* collected from Oued el Heimer site. The variation of Pb accumulation in the roots of *R. alba*, and the genus of *Stipa* can significantly be related to the physicochemical properties of soils.

In shoot, the concentrations of Pb vary from 47.1 mg kg⁻¹ DW in *R. rigosum* to 4672. 2 mg kg⁻¹ DW in *A. herba alba* (belonging to Asteraceae family) (Table 3). Pb content found in *A. herba alba* was significantly higher than those observed in *R. alba*, *C. libanotis*, *S. tenacissima*, and other plant species (Tables 3 and 4). Many authors report high heavy metals hyperaccumulation ability in several plants belonging to the Asteraceae family. Nouri et al. reported high level of Pb up to 9017 mg kg⁻¹ DW in shoots of *Scariola orientalis* growing in the vicinity of the Ahangaran lead-zinc mine (Hamedan, Iran) [61]. *Aster prorerus Hemsl*, a plant belonging to the Asteraceae family, collected near the smelter of the Baoshan mining project (China), showed a Pb concentration up to 3677 mg kg⁻¹ DW in shoots [79]. In our study, four plant species namely *A. herba-alba* (4672.2 mg kg⁻¹ DW), *R. alba* (1607.51 mg kg⁻¹ DW), *C. libanotis* (1261.8 mg kg⁻¹ DW), and *S. tenacissima* (1146.3 mg kg⁻¹ DW) had a Pb concentration higher than 1000 mg kg⁻¹ DW in their shoots, which exceeded the threshold for Pb hyperaccumulation in plants [14,75]. All studied plants sampled from both sites showed Pb concentrations higher than the phytotoxic one (0.6–28 mg kg⁻¹) which means that they are tolerant to this metal.

The BCF value varied from 0.02 (*R. alba*) to 1.99 in *C. libanotis*. While the TF > 1 was recorded in six plant species namely *R. alba* (4.98), *A. herba alba* (2.67), *C. bursa-pastoris* (Oued el Heimer) (1.46), *A. alopecuroides* (1.24), *S. hispanicus* L. (1.21), and *C. libanotis* (1.03). The minimum and maximum values of BAC found in *R. rigosum* (0.001) and *C. libanotis* (2.06), respectively. According to Yoon et al. [16], three plant species from Oued el Heimer site have the potential for Pb phytoextraction. *C. libanotis* with BCF, TF, and BAC superior to 1, followed by *S. tenacissima*, which although it accumulates more than 1000 mg kg⁻¹ of Pb in its shoots, does not meet the second criterion for hyperaccumulation (TF > 1). This plant remains a good candidate for Pb phytoextraction. The third plant is *A. herb-alba* which accumulates up to 1000 mg kg⁻¹ in its shoots with a high translocation factor (2.67). *R. alba* from Touissite also can be used in Pb phytoextraction because it accumulates up to >1000 of Pb mg kg⁻¹ in shoot parts with TF > 1. This result indicates that this plant species exhibited a good capacity to translocate Pb from roots to shoots. According to the criteria of phytostabilization (BCF > 1 and TF < 1), none of the collected plants has the potential for Pb phytostabilization [13].

### 3.4.6. Zinc (Zn)

Zn content in shoots and roots were significantly different between species (Tables 3 and 4). Zn concentrations in roots range between 135.6 and 637 mg kg⁻¹ DW in *C. libanotis* and *S. tenacissima*, respectively (Table 3). However, no significant difference in Zn level in the roots of *S. tenacissima*, *A. herba-alba*, and *C. bursa-pastoris* (Touissite). The maximum Zn accumulated in shoot was 666.9 mg kg⁻¹ DW in *C. bursa-pastoris* collected from Touissite. Zn content in *C. bursa-pastoris* shoots (from Touissite) was significantly higher compared to other 13 plant species. Nine plant species could accumulate concentrations of Zn higher than the phytotoxic levels (100–300 mg kg⁻¹) in their tissues and were considered as Zn tolerant species, i.e., *R. alba*, *S. hispanicus*, *P. communis*, *C. bursa-pastoris* (from both site), *H. spinosissimum*, *R. rigosum*, and *S. tenacissima*, *A. herba-alba* (Table 3). Many authors recorded several plant species growing on contaminated soils that could accumulate a high concentration of Zn in different plant tissues, especially in roots. Nouri et al. reported 73 and 2938 mg Zn kg⁻¹ DW in roots of *Reseda alba* L and *Reseda lutea* L, collected in the vicinity of a Pb-Zn mine in Iran [61]. Midhat et al. [8] reported 856.19 mg Zn kg⁻¹ DW in *Stipa capensis* Thunb roots sampled from Sidi Bou-Othmane mining site (Marrakech city, Morocco). Regarding hyperaccumulation criteria, none of the plant species collected from the two areas exceeds the threshold of Zn hyperaccumulation > 10,000 mg kg⁻¹ in their shoots and then cannot be considered as Zn hyper accumulator [14].
The bioconcentration factor (BCF) varied from 0.01 in *H. incana* and *S. tenacissima* to 1.69 in *A. herba-alba*. Nine plant species have the TF > 1 namely *P. communis* (2.17), *R. alba* (1.83), *S. hispanicus* L. (1.46), *R. rigosum* (1.36), *C. libanotis* (1.19), *A. alopecuroides* (1.18), *C. bursa-pastoris* (Touissite) (1.13), *H. incana* (1.13), and *C. bursa-pastoris* (Oued el Heimer) (1.00) (Table 5). The highest BAC value was observed in *A. herba-alba* (1.04), and the lowest value was reported in *S. tenacissima* (0.001). According to Yoon et al. [16], none of the harvested plant species from both areas has the potential for phytoextraction of Zn (BCF, TF, and BAC > 1). *A. herba-alba*, with BCF > 1 and TF < 1, could be a suitable candidate for Zn phytostabilization to minimize Zn migration in soils [16,62].

### 3.4.7. Antimony (Sb)

Sb content in shoots and roots were significantly different between species. The maximum Sb accumulated in roots was 167.7 mg kg⁻¹ DW in *S. tenacissima* and the minimum was 1.6 mg kg⁻¹ DW in *R. rigosum*. Sb content in roots of *S. tenacissima* was significantly higher than those observed in roots of other plants species (Tables 3 and 4). *C. bursa-pastoris* (Touissite) and *A. herba-alba* accumulated the highest levels of Sb in their shoots (55.5 and 41.6 mg kg⁻¹ DW, respectively), there was no significant between those two species, but a significant variation of Sb was observed in shoots between other plant species (Table 3). However, the lowest concentration of Sb 0.5 mg kg⁻¹ DW in shoots was found in *R. rigosum* (Table 3). These plant species were grown in soils with maximum Sb concentrations of 247.1 and 242.8 mg kg⁻¹ in Oued el Heimer and Touissite, respectively. Many authors indicate that the high levels of Sb in soils do not necessarily result in significant accumulation of Sb in plant tissues [80]. The concentration of this metal in plants depends on three principal factors, i.e., 1: the Sb bioavailability in soils; 2: Sb speciation; and 3: variations of coexisting ions in soils such as calcium (Ca) and phosphorus (P) [80]. In the present study, ten plant species accumulated more than phytotoxic concentration of Sb (5–10 mg kg⁻¹) in their tissues and are expected as Sb tolerant plants, i.e., *R. alba*, *C. althaeoides*, *H. spinosissimum*, *P. communis*, *L. corniculatus*, *C. bursa-pastoris* (Touissite), *S. hispanicus* L., *S. tenacissima*, *A. herba-alba*, and *C. bursa-pastoris* (Oued el Heimer). All plant species did not accumulate >1000 mg kg⁻¹ in their shoots, meaning that none of the harvested plants could be considered as Sb hyperaccumulator.

All plant species sampled from both sites exhibited a BCF less than 1 (Table 5). However, nine plant species showed a TF > 1 namely *R. alba* (3.55), *A. herba-alba* (3.2), *C. bursapastoris* (Touissite) (2.05), *S. hispanicus* L. (1.39), *A. alopecuroides* (1.19), *P. communis* (1.15), *H. incana* (1.2), and *C. libanotis* (1.13) and only one plant species, i.e., *A. herba-alba* (1.22) has recorded BAC > 1 (Table 5). According to Yoon et al. [16] none of the fourteen plant species is suitable for phytostabilization or phytoextraction of Sb.

In general, the majority of collected plant species from both prospected sites demonstrate a high ability to survive in this environment. However, all of them presented high heavy metal concentrations in different plant tissues that exceeded the phytotoxic levels without sustaining toxicity. BCF, TF, and BAC values indicate that only four plant species exhibited strong heavy metals tolerance and hyperaccumulation and then have a good potential to be used in phytoremediation strategies.

### 3.5. *Plants Polymetallic Accumulation Ability*

The biplot one-half of the principal components analysis (PCA) performed on metals concentrations in the shoots of the studied plant species (Figure 4A) explain 83.33% of the total variance. The F1 axis (58.78% of the total inertia) represent the species ability to accumulate different heavy metals simultaneously. It separated *Artemisia herba-alba* with high concentrations of As, Cd, Cu, Ni, and Pb in shoots from the other species (*Rapistrum rigosum*, *Convolvulus altheoides*, *Lotus corniculatus*, *Hirschfeldia incana*, etc.) with low concentrations of these metals in the shoots (Figure 4A). The F2 axis (24.55% of the total inertia) represents Zn and Sb concentration gradients (Figure 4A). It separates *Capsella bursa-pastoris* (from Touissite site), *Phragmites communis*, *Reseda alba*, *Scolymus hispanica*, and *Hedysarum spinosissimum* which are highly correlated with Zn and Sb concentrations in shoots, from the other species weakly correlated with Zn and Sb concentrations (Figure 4A). The clusters of plants species
identified with PCA was confirmed by dendrogram obtained by ascending hierarchical classification (AHC) (Figure S1A).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Plot of axes 1 and 2 of the principal components analysis (PCA) on the means of the metals concentrations separately in the shoots (A) and the roots (B) for the studied species. Abbreviation of species names: Ag.Alo: Agathophora alopecuroides; Ci.lib: Cistus libanotus; Ca.bur_T: Capsella bursa-pastoris Touissite site; Ca.bur_O: Capsella bursa-pastoris Oued El Heimer site; Lo.cor: Lotus corniculatus; Ph.com: Phragmites communis; He spi: Hedysarum spinosissimum; Co.alt: Convolvulus althaeoides; Re.alb: Reseda alba; Sc.his: Scolymus hispanicus; Ra.rig: Rapistrum rigosum; Hi inc: Hirschfeldia incana; St.ten: Stipa tenacissima; Ar.her: Artemisia herba-alba.</figcaption></figure>

The biplot one-half of the PCA performed on metals concentrations in roots of studied species (Figure 4B) explain 81.03% of the total variance. The F1 axis (77.54% for the total inertia) explains the major part of the information which represents the species ability to accumulate different heavy metals simultaneously (Figure 4B). According to this axis, *Stipa tenacissima* followed by *Artemisia herba-alba* presented a high correlation with As, Cd, Cu, Ni, Pb, Sb, and Zn indicating that these species have strong accumulation ability on roots of all studied metals (Figure 4B). The other group (*Capsella bursa-pastoris* "Touissite", *Reseda alba*, and *Hedysarum spinosissimum*) form another cluster (Figure S1B). This shows that these species have a preferential accumulation capacity for specific metals (e.g., *Cistus libanotis* for Pb (1607.5 mg kg⁻¹ DW); *Stipa tenacissima* for Pb (1146.3 mg kg⁻¹ DW))

In our knowledge, the uptake and translocation of metal ions by the plants depends not only on the concentration of these elements in the soil, but also on metal plant affinity, plant species and metal speciation [81]. Indeed, biotic (microorganisms)and abiotic soil characteristics (physical and chemical properties), have a significant influence on heavy metals content in plants [82]. As well, there is often an antagonism between ions in plants. As an example, antagonistic decrease in the uptake and accumulation of heavy metals in plants, e.g., ions of Pb or Cd was observed with the increase of the concentration of trace elements in the soil, e.g., Zn [83]. The potential antagonism between Pb and Zn occurs only in the case of some plants. According to our results, native wild plants *Stipa tenacissima* and *Artemisia herba-alba* showed cooperative accumulation of all metals in roots. Interestingly, *Artemisia herba-alba* demonstrated strong translocation capacity in shoots for As, Cd, Cu, Ni, and Pb.
## 3.6. Suitable Plant Species for Use in Phytoremediation Strategies

### 3.6.1. *Stipa tenacissima* Suitable Plant for Cd/Cu Phytostabilization and Pb Phytoextraction

Among all plant species, only *Stipa tenacissima* from Oued el Heimer presented a high amount of Cd and Cu up to 241.2 mg kg⁻¹ DW and 237.9 mg kg⁻¹ DW, in roots respectively, with TF < 1. *Stipa tenacissima* is a perennial plant belonging to Poaceae family, commonly known as "Halfa" [84]. This species has already identified in other mining sites. Some authors reported that *Stipa tenacissima* species from Zaida mining site (Morocco) could accumulate up to 58.4 and 41.8 mg kg⁻¹ DW of Cu in their roots and shoots, respectively [85]. Other species of genus *Stipa* from Sidi Bou-Othmane mining site (Marrakech, Morocco) accumulate 50.88 mg kg⁻¹ DW of Cu and 0.87 mg kg⁻¹ DW of Cd in their roots [8], with roots without being affected. Regarding the BCF and TF values, *Stipa tenacissima* presents a hyper tolerance characteristic with an excellent capacity to accumulate Cd and Cu in their roots. This native plant could be a suitable candidate for Cd/Cu phytostabilization.

*Stipa tenacissima* also accumulates up to 1226.9 mg kg⁻¹ of Pb in their shoots. This concentration is higher than 193.2 mg kg⁻¹ of Pb recorded in *Stipa tenacissima* plants harvested at Zaida Pb mining site (high Moulouya, Morocco) [86]. Our result showed that Pb accumulated in *Stipa tenacissima* shoots exceeds the hyperaccumulation thresholds, more than 1000 mg Pb/kg DW. Since the roots accumulated more Pb up to 3785.7 mg kg⁻¹ than the shoots (TF < 1), the term hyper accumulator is not fully appropriate. However, it appears that *Stipa tenacissima* exhibits a good Pb accumulation capacity and can be a suitable candidate for Pb phytoextraction.

### 3.6.2. *A. herba-alba* Suitable Plant for Cu/Zn Phytostabilization and Pb Phytoextraction

The genus *Artemisia* belongs to Asteraceae family, that include many species known as accumulators of several metals [86–88]. In the present study, high amounts of Cu up to 203.1 mg kg⁻¹ DW and Zn up to 576.9 mg kg⁻¹ DW have been observed in the roots of *A. herba-alba*. The Cu and Zn concentration accumulated by *A. herba-alba* roots were higher than those reported by Baghdad et al. [85] at Zaida mining site (56.2 mg kg⁻¹ DW of Cu and 133.9 mg kg⁻¹ DW of Zn). Some *Artemisia* species are widespread on contaminated sites and have a large biomass production and high ability to accumulate heavy metals, especially Zn [88]. Alirzayeva et al. [88] reported up to 200.71 mg kg⁻¹ DW of Zn in roots of *A. scoparia*, and up to 112.24 mg kg⁻¹ DW in roots of *A. fragrans* was reported by Ashraf et al., [89]. In our study, *A. herba-alba* characterized by BCF > 1 and a high amount of Zn and Cu in its roots and is then a good candidate for Cu/Zn phytostabilization.

*Artemisia herba-alba* accumulates also a large amount of Pb in roots (1748.7 mg kg⁻¹ DW) and shoots (4672.2 mg kg⁻¹ DW). These amounts are higher than those recorded in the same plant species growing at different Pb mining sites in Morocco (704.0 and 660.0 mg kg⁻¹ DW of Pb in its roots and shoots, respectively) [85]. *Artemisia lancangensis*, another plant from *Artemisia genus* growing at the lead-zinc mine area in China, accumulates 12.83 mg kg⁻¹ DW of Pb in their shoots [90]. Other authors recorded 892 mg kg⁻¹ DW in root and 834 mg kg⁻¹ DW Pb in the shoot of *Artemisia lactiflora* Wall from Beiya Pb mine area, China [91], 11.58 mg kg⁻¹ DW in the shoot of *Artemisia vulgaris* growing near mining and industry in upper Silesia (southern Poland) [92] and 375, 446, 539 and 3677 mg kg⁻¹ DW of Pb in *Artemisia taurica*, *Artemisia argyi*, *Artemisia japonica*, and *Aster prorerus*, respectively growing on polluted soils, with Pb concentrations ranging from 568 to 49,294 mg kg⁻¹ at the Baoshan mining area (China) [79]. The higher amount of Pb in the shoot of *Artemisia herba-alba* can be due to the strong capacity of metal uptake and an enhanced xylem loading capacity for metals [79]. According to previous studies, this is the first time that *A. herba-alba* is identified as hyper-tolerant and hyper accumulator of Pb. This is not a general characteristic of Asteraceae family. Our result showed that Pb accumulated in *Artemisia herba-alba* exceeds the hyperaccumulation threshold of 1000 mg Pb kg⁻¹ DW. Since the shoots accumulated more Pb than roots (TF = 2.67), the term hyper accumulator is fully appropriate, which means that *A. herba-alba* evolved high Pb transport ability from the roots to the shoots, and can be suggested as a good candidate for Pb phytoextraction.
### 3.6.3. *Reseda alba* a Powerful Plant Species for Pb Phytoextraction

R. alba belonging to the Resedaceae family, was known for its potential to accumulate Pb, Zn, and Cd in roots and shoots [13,93]. In our study, this species accumulates 1607.5 mg kg⁻¹ DW of Pb in their shoots, with high TF = 4.98. However, Nouri et al. [61] indicate that *Reseda alba* from the Ahangaran lead-zinc mine in Hamedan, Iran (soil: 9535 mg kg⁻¹ of Pb), can accumulate up to 703 mg kg⁻¹ DW of Pb in aboveground parts. The change in soil conditions can have many impacts on the speciation and distribution of metals [41]. For example, pH plays an important role in the solubility and mobility of heavy metals. Bliefert et al. [39] reported that pH more than 7.5 decreases significantly the mobility of Pb while the mobility of this metal is greater in acidic soils. To our knowledge, this is the first time that *Reseda alba* is described as high Pb accumulator since its Pb content exceeded the hyperaccumulator threshold of 1000 mg kg⁻¹ DW and TF > 1. R. alba then could be an excellent candidate for the Pb phytoextraction.

### 3.6.4. *Cistus libanotis* a Good Candidate for Pb Phytoextraction

*Cistus libanotis* is a fast-growing perennial plant with high biomass [94]. In the present study, *Cistus libanotis* can accumulate up to 1261.8 mg kg⁻¹ DW of Pb in shoots, with TF > 1. These results confirmed that *Cistus libanotis* is a good candidate for use in phytoremediation, especially in phytoextraction. *Cistus libanotis* was reported previously by our team as an excellent plant for Pb phytoextraction [95]. Plants belonging to *Cistus* genus show high plasticity. They can grow both in contaminated and non-contaminated soils being tolerant to high concentrations of trace elements, like As and Pb [96,97]. *C. albidus* and *C. salviifolius* respectively accumulate higher amounts of Pb up to 4831 and 2896 mg kg⁻¹ of Pb in their roots than in their shoots and were considered suitable for phytostabilization of Pb contaminated soils [94]. The same study indicates that *Cistus libanotis* can accumulate a large amount of Pb in its roots and shoots (10,340 and 2055 mg Pb kg⁻¹ DW, respectively) when grown in hydroponic conditions [94]. Santos et al. [98] suggest that *C. ladanifer* is suitable for phytostabilization of mine soils rich in Pb, As, Cu, and Zn. However, metal accumulation is not a general pattern of this genus. In the present study, *Cistus libanotis* collected from Oued el Heimer site can be a good candidate for Pb hyperaccumulation.

## 4. Conclusions

Long term mining activities in Oued el Heimer and Touissite areas in eastern Morocco has caused strong heavy metal contaminations related to massive hazardous and poorly managed waste deposits. Soils are characterized by high levels of contamination especially Pb, Cd, Zn, and the critical metal Sb. Native plant species identified in these areas present a high ability to tolerate and accumulate several metals in their tissues. Among the 14 plant species collected from these areas, only four plants are identified as Pb hyperaccumulators, namely *Reseda alba*, *Stipa tenacissima*, *Artemisia herba-alba*, and *Cistus libanotis*. Regarding the high amounts of metals in their roots and shoots and BCF, TF, and BAC values, these plant species, classified as hyper-tolerant with a strong capacity to accumulate several metals, are promising native plants for Pb phytoextraction. *Stipa tenacissima* and *Artemisia herba-alba* have the potential to be used in phytostabilization for Cd/Cu and Cu/Zn, respectively. Further investigations will be led to increase the value of metal-phytomining potential from Oued el Heimer and Touissite soils using suitable selected plants species.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2223-7747/9/11/1458/s1, Table S1: Heavy metal concentrations and reference element (mgKg⁻¹ DW) in the rhizospheric soils samples collected from studied area in the oriental region of Morocco., Figure S1. Dendrogram obtained by Ascending Hierarchical Classification (AHC) carried out on the metals concentrations in the shoots (A) and the roots (B) of the studied plants.

**Author Contributions:** Conceptualization, M.F. and A.S.; Investigation, S.E.H., B.A., and P.C.; Formal analysis, L.R.; resources, M.F., C.K., F.C., and A.S.; writing—original draft preparation, S.E.H., M.F., and A.S.; writing—review and editing, M.F., C.K., C.L., P.C., Z.E.A.T., A.G., L.R., F.C., and A.S.; supervision, M.F., and A.S.; project
administration, M.F., F.C., and A.S.; funding acquisition, M.F., C.K., F.C., and A.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** The APC was funded by "Institute of Research for Development—France IRD: JEAI—PhytoMet 2018–2020".

**Acknowledgments:** Published with the support of "IRD: JEAI—PhytoMet 2018–2020; LMI ‘AMIR’ 2020–2025", APIC 2019 CEREGE—France and Programme Prioritaire de recherche PPR- Ministère de l’Enseignement Supérieur, de la Recherche Scientifique et de la Formation des Cadres. We thank the four anonymous reviewers and editors for their time and fruitful comments.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**References**

1. Baize, D.; Tercé, M. *Les Éléments Traces Métalliques Dans Les Sols: Approches Fonctionnelles et Spatiales*; Editions Quae: Versailles, France, 2002.

2. Adriano, D.C. *Trace Elements in Terrestrial Environments*; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2001; pp. 219–261.

3. Singh, A.; Sharma, R.K.; Agrawal, S.B.; Marshall, F.M. Health risk assessment of heavy metals via dietary intake of foodstuffs from the wastewater irrigated site of a dry tropical area of India. *Food Chem. Toxicol.* **2010**, *48*, 611–619. [CrossRef]

4. Saxena, G.; Purchase, D.; Mulla, S.I.; Saratale, G.D.; Bharagava, R.N. Phytoremediation of Heavy Metal-Contaminated Sites: Eco-environmental Concerns, Field Studies, Sustainability Issues, and Future Prospects. *Rev. Environ. Contam. and Toxicol. Vol. 201* **2019**, *249*, 71–131. [CrossRef]

5. Sanyaolu, V.; Adeniran, A. Determination of Heavy Metal Fallout on the Surrounding Flora and Aquifer: Case Study of a Scrap Metal Smelting Factory in Odogunyan Area, Ikorodu, Lagos-State, Nigeria. *Int. Res. J. Environ. Sci.* **2014**, *3*, 93–100.

6. Angelovičová, L.; Fazekašová, D. Contamination of the soil and water environment by heavy metals in the former mining area of Rudňany (Slovakia). *Soil Water Res.* **2014**, *9*, 18–24. [CrossRef]

7. Lenart-Boroń, A.; Wolny-Koładka, K. Heavy metal concentration and the occurrence of selected microorganisms in soils of a steelworks area in Poland. *Plant, Soil Environ.* **2016**, *61*, 273–278. [CrossRef]

8. Midhat, L.; Ouazzani, N.; Hejjaj, A.; Ouhammou, A.; Mandi, L. Accumulation of heavy metals in metallophytes from three mining sites (Southern Centre Morocco) and evaluation of their phytoremediation potential. *Ecotoxicol. Environ. Saf.* **2019**, *169*, 150–160. [CrossRef]

9. Cao, R.X.; Ma, L.Q.; Chen, M.; Singh, S.P.; Harris, W.G. Phosphate-induced metal immobilization in a contaminated site. *Environ. Pollut.* **2003**, *122*, 19–28. [CrossRef]

10. Mulligan, C.; Yong, R.; Gibbs, B. Remediation technologies for metal-contaminated soils and groundwater: An evaluation. *Eng. Geol.* **2001**, *60*, 193–207. [CrossRef]

11. Mahar, A.; Wang, P.; Ali, A.; Awasthi, M.K.; Lahori, A.H.; Wang, Q.; Li, R.; Zhang, Z. Challenges and opportunities in the phytoremediation of heavy metals contaminated soils: A review. *Ecotoxicol. Environ. Saf.* **2016**, *126*, 111–121. [CrossRef] [PubMed]

12. McGrath, S.; Zhao, F.-J. Phytoextraction of metals and metalloids from contaminated soils. *Curr. Opin. Biotechnol.* **2003**, *14*, 277–282. [CrossRef]

13. Lorestani, B.; Cheraghi, M.; Yousefi, N. Phytoremediation potential of native plants growing on a heavy metals contaminated soil of copper mine in Iran. *World Acad. Sci. Eng. Technol.* **2011**, *77*, 377–382.

14. Baker, A.J.; Brooks, R. Terrestrial higher plants which hyperaccumulate metallic elements. A review of their distribution, ecology and phytochemistry. *Biorecovery* **1989**, *1*, 81–126.

15. Terry, N. Phytoremediation of selenium-contaminated soil and water. *Glob. Adv. Selenium Res. Theory Appl.* **2015**, 197–198. [CrossRef]

16. Yoon, J.; Cao, X.; Zhou, Q.; Ma, L.Q. Accumulation of Pb, Cu, and Zn in native plants growing on a contaminated Florida site. *Sci. Total Environ.* **2006**, *368*, 456–464. [CrossRef]

17. Smouni, A.; Ater, M.; Auguy, F.; Laplaze, L.; El Mzibri, M.; Berhada, F.; Maltouf, A.F.; Doumas, P. Évaluation de la contamination par les éléments-traces métalliques dans une zone minière du Maroc oriental. *Cah. Agric.* **2010**, *19*, 273–279. [CrossRef]
18. Fennane, M.; Ibn Tattou, M.; El Oulaidi, J. Flore pratique du Maroc. In *Travaux de l'Institut Scientifique*; Série Botanique; Travaux de l'Institut Scientifique: Rabat, Morocco, 2014; Volume 3, p. 40.

19. Fennane, M.; Ibn Tattou, M.; Mathez, J.; Ouyahya, A.; El Oualidi, J. Flore Prafique du Maroc. Manuel de determinafion des plantes vasculaires. In *Travaux de L'institut Scientifique*; Série Botanique; Travaux de l'Institut Scientifique: Rabat, Morocco, 1999; Volume 1, p. 36.

20. Fennane, M.; Ibn Tattou, M.; Mathez, J.; Ouyahya, A.; El Oualidi, J. Flore pratique du Maroc. In *Travaux Institut Scientifique, Vol. 2*; Série Botanique N°38; Travaux de l'Institut Scientifique: Rabat, Morocco, 2007; Volume 38.

21. Combs, S.; Nathan, M. Soil organic matter. Recommended chemical soil test procedures for the north central region. *North Cent. Reg. Res. Pub.* **1998**, *221*, 53–58.

22. Temminghoff, E.E.; Houba, V.J. *Plant Analysis Procedures*; Springer: Berlin/Heidelberg, Germany, 2004; Volume 179.

23. Marguí, E.; Queralt, I.; Carvalho, M.; Hidalgo, M. Comparison of EDXRF and ICP-OES after microwave digestion for element determination in plant specimens from an abandoned mining area. *Anal. Chim. Acta* **2005**, *549*, 197–204. [CrossRef]

24. Chon, H.-T.; Ahn, J.-S.; Jung, M.C. Seasonal Variations and Chemical Forms of Heavy Metals in Soils and Dusts from the Satellite Cities of Seoul, Korea. *Environ. Geochem. Health* **1998**, *20*, 77–86. [CrossRef]

25. Wu, J.; Teng, Y.; Lu, S.; Wang, Y.; Jiao, X. Evaluation of Soil Contamination Indices in a Mining Area of Jiangxi, China. *PLoS ONE* **2014**, *9*, e112917. [CrossRef]

26. Kabata-Pendias, A.; Pendias, H. *Trace Elements in Plants and Soils*; CRC Press: Boca Raton, FL, USA, 1984; pp. 233–237.

27. Bern, C.R.; Walton-Day, K.; Naftz, D.L. Improved enrichment factor calculations through principal component analysis: Examples from soils near breccia pipe uranium mines, Arizona, USA. *Environ. Pollut.* **2019**, *248*, 90–100. [CrossRef]

28. Blaser, P.; Zimmermann, S.; Luster, J.; Shotyk, W. Critical examination of trace element enrichments and depletions in soils: As, Cr, Cu, Ni, Pb, and Zn in Swiss forest soils. *Sci. Total Environ.* **2000**, *249*, 257–280. [CrossRef]

29. Zhang, L.; Ye, X.; Feng, H.; Jing, Y.; Ouyang, T.; Yu, X.; Liang, R.; Gao, C.; Chen, W. Heavy metal contamination in western Xiamen Bay sediments and its vicinity, China. *Mar. Pollut. Bull.* **2007**, *54*, 974–982. [CrossRef]

30. Liu, Y.; Li, S.; Sun, C.; Qi, M.; Yu, X.; Zhao, W.; Li, X. Pollution Level and Health Risk Assessment of PM2.5-Bound Metals in Baoding City Before and After the Heating Period. *Int. J. Environ. Res. Public Health* **2018**, *15*, 2286. [CrossRef]

31. Ahmadpour, P.; Ahmadpour, F.; Mahmud, T.M.M.; Abdu, A.; Soleimani, M.; Tayefeh, F.H. Phytoremediation of heavy metals: A green technology. *Afr. J. Biotechnol.* **2012**, *11*, 14036–14043. [CrossRef]

32. Nazir, A.; Malik, R.N.; Ajaib, M.; Khan, N.; Siddiqui, M.F. Hyperaccumulators of heavy metals of industrial areas of Islamabad and Rawalpindi. *Pak. J. Bot.* **2011**, *43*, 1925–1933.

33. Lamin, H.; Alami, S.; Bouhnik, O.; Elfaik, S.; Abdelmoumen, H.; Bedmar, E.J.; Idrissi, M.M.-E. Nodulation of Retama monosperma by Ensifer aridi in an Abandonned Lead Mine Soils in Eastern Morocco. *Front. Microbiol.* **2019**, *10*, 1456. [CrossRef]

34. Claveau, J.; Paulhac, J.; Pellerin, J. The lead and zinc deposits of the Bou Beker-Touissit area, eastern French Morocco. *Econ. Geol.* **1952**, *47*, 481–493. [CrossRef]

35. Dupuy, J.-J.; Touray, J.-C. Multistage ore deposition at the Oued Mekta strata-bound lead deposit, Touissit-Bou Beker District, eastern Morocco. *Econ. Geol.* **1986**, *81*, 1558–1561. [CrossRef]

36. Makhoukhi, S.; Marignac, C.; Pironon, J.; Schmitt, J.; Marrakchi, C.; Bouabdelli, M.; Bastoul, A. Aqueous and hydrocarbon inclusions in dolomite from Touissit-Bou Beker district, Eastern Morocco: A Jurassic carbonate hosted PbZn (Cu) deposit. *J. Geochem. Explor.* **2003**, *78*, 545–551. [CrossRef]

37. Navarro, M.; Pérez-Sirvent, C.; Martínez-Sánchez, M.; Vidal, J.; Tovar, P.; Bech, J. Abandoned mine sites as a source of contamination by heavy metals: A case study in a semi-arid zone. *J. Geochem. Explor.* **2008**, *96*, 183–193. [CrossRef]

38. Ma, X.; Zuo, H.; Tian, M.; Zhang, L.; Meng, J.; Zhou, X.; Min, N.; Chang, X.; Liu, Y. Assessment of heavy metals contamination in sediments from three adjacent regions of the Yellow River using metal chemical fractions and multivariate analysis techniques. *Chemosphere* **2016**, *144*, 264–272. [CrossRef] [PubMed]
39. Bliefert, C.; Perraud, R. *Chimie de l'environnement: Air, Eau, Sols, Déchets*, 1st ed.; De Boeck Université: Bruxelles, Belgique, 2001.

40. Yassir, B.; Sana, E.F.; Alain, P. Contamination by trace elements of agricultural soils around Sidi Bou Othmane in abandoned mine tailings in Marrakech, Morocco. *Pollution Winter.* **2016**, *2*, 93–101.

41. Zhu, G.; Xiao, H.; Guo, Q.; Song, B.; Zheng, G.; Zhang, Z.; Zhao, J.; Okoli, C.P. Heavy metal contents and enrichment characteristics of dominant plants in wasteland of the downstream of a lead-zinc mining area in Guangxi, Southwest China. *Ecotoxicol. Environ. Saf.* **2018**, *151*, 266–271. [CrossRef] [PubMed]

42. Barbieri, M. The Importance of Enrichment Factor (EF) and Geoaccumulation Index (Igeo) to Evaluate the Soil Contamination. *J. Geol. Geophys.* **2016**, *5*, 1–4. [CrossRef]

43. El Hachimi, M.L.; Fekhaoui, M.; El Abidi, A.; Rhoujatti, A. Contamination des sols par les métaux lourds à partir de mines abandonnées: Le cas des mines Aouli-Mibladen-Zeïda au Maroc. *Cah. Agric.* **2014**, *23*, 213–219.

44. Lorestani, B.; Cheraghi, M.; Yousefi, N. The Potential of Phytoremediation Using Hyperaccumulator Plants: A Case Study at a Lead-Zinc Mine Site. *Int. J. Phytoremediation* **2012**, *14*, 786–795. [CrossRef]

45. Escarré, J.; Lefèbvre, C.; Raboyeau, S.; DosSantos, A.; Gruber, W.; Marel, J.C.C.; Frérot, H.; Noret, N.; Mahieu, S.; Collin, C.; et al. Heavy Metal Concentration Survey in Soils and Plants of the Les Malines Mining District (Southern France): Implications for Soil Restoration. *Water Air Soil Pollut.* **2010**, *216*, 485–504. [CrossRef]

46. He, M. Distribution and phytoavailability of antimony at an antimony mining and smelting area, Hunan, China. *Environ. Geochem. Health* **2007**, *29*, 209–219. [CrossRef]

47. Okkenhaug, G.; Zhu, Y.-G.; Luo, L.; Lei, M.; Li, X.; Mulder, J. Distribution, speciation and availability of antimony (Sb) in soils and terrestrial plants from an active Sb mining area. *Environ. Pollut.* **2011**, *159*, 2427–2434. [CrossRef]

48. Li, J.; Wei, Y.; Zhao, L.; Zhang, J.; Shangguan, Y.; Li, F.; Hou, H. Bioaccessibility of antimony and arsenic in highly polluted soils of the mine area and health risk assessment associated with oral ingestion exposure. *Ecotoxicol. Environ. Saf.* **2014**, *110*, 308–315. [CrossRef]

49. Takaoka, M.; Fukutani, S.; Yamamoto, T.; Horiuchi, M.; Satta, N.; Takeda, N.; Oshita, K.; Yoneda, M.; Morisawa, S.; Tanaka, T. Determination of Chemical Form of Antimony in Contaminated Soil around a Smelter Using X-ray Absorption Fine Structure. *Anal. Sci.* **2005**, *21*, 769–773. [CrossRef] [PubMed]

50. Al-Shehbaz, I.A. A generic and tribal synopsis of the Brassicaceae (Cruciferae). *TAXON* **2012**, *61*, 931–954. [CrossRef]

51. Thanos, C.A.; Georghiou, K.; Kadis, C.; Pantazi, C. Cistaceae: A plant family with hard seeds. *Isr. J. Plant Sci.* **1992**, *41*, 251–263.

52. Al-Sghair, F.G.; Mahklouf, M.H.; Abudaya, E.A. Species Diversity and Floristic Analysis of the Family Poaceae in Libya Depending on the Flora of Libya. *Adv. Biosci. Bioeng.* **2019**, *7*, 13. [CrossRef]

53. Martín-Bravo, S.; Escudero, M. Biogeography of flowering plants: A case study in mignonettes (Resedaceae) and sedges (Carex, Cyperaceae). In *Global Advances in Biogeography*; Intech: Rijeka, Croatia, 2012; pp. 257–290.

54. Hassan, N.S.; Hartmann, H.E.K.; Liede-Schumann, S. Conspectus of Aizoaceae, Gisekiaceae and Molluginaceae of Egypt and the Sudan. *Feddes Repert.* **2005**, *116*, 1–42. [CrossRef]

55. Aykurt, C.; Suembuel, H. Varieties and chorology of Convolvulus oleifolius Desr. (Convolvulaceae) in Turkey. *Biol. Divers. Conserv.* **2010**, *3*, 155–162.

56. Brullo, S.; Minissale, P.; Salmeri, C.; Del Galdo, G.G. Taxonomic investigations on Psoralea palaestina (Fabaceae), a critical species of Mediterranean flora. *Phytotaxa* **2016**, *266*, 61. [CrossRef]

57. Bengtson, A.; Anderberg, A.A. Species diversification in the Mediterranean genus Chiliadenus (Inuleae-Asteraceae). *Plant. Syst. Evol.* **2018**, *304*, 853–860. [CrossRef]

58. Gisbert, C.; Clemente, R.; Navarro-Aviñó, J.; Baixauli, C.; Ginér, A.; Serrano, R.; Walker, D.J.; Bernal, M.P. Tolerance and accumulation of heavy metals by Brassicaceae species grown in contaminated soils from Mediterranean regions of Spain. *Environ. Exp. Bot.* **2006**, *56*, 19–27. [CrossRef]

59. Casado, M.; Anawar, H.M.; Garcia-Sanchez, A.; Santa-Regina, I. Antimony and Arsenic Uptake by Plants in an Abandoned Mining Area. *Commun. Soil Sci. Plant. Anal.* **2007**, *38*, 1255–1275. [CrossRef]

60. Wenzel, W.; Jockwer, F. Accumulation of heavy metals in plants grown on mineralised soils of the Austrian Alps. *Environ. Pollut.* **1999**, *104*, 145–155. [CrossRef]
61. Nouri, J.; Lorestani, B.; Yousefi, N.; Khorasani, N.; Hasani, A.H.; Seif, F.; Cheraghi, M. Phytoremediation potential of native plants grown in the vicinity of Ahangaran lead-zinc mine (Hamedan, Iran). *Environ. Earth Sci.* **2010**, *62*, 639-644. [CrossRef]

62. Mahdavian, K.; Ghaderian, S.M.; Torkzadeh-Mahani, M. Accumulation and phytoremediation of Pb, Zn, and Ag by plants growing on Koshk lead-zinc mining area, Iran. *J. Soils Sediments* **2015**, *17*, 1310-1320. [CrossRef]

63. Marrugo-Negrete, J.; Marrugo-Madrid, S.; Pinedo-Hernández, J.; Durango-Hernández, J.; Díez, S. Screening of native plant species for phytoremediation potential at a Hg-contaminated mining site. *Sci. Total Environ.* **2016**, *542*, 809-816. [CrossRef]

64. Ha, N.T.H.; Sakakibara, M.; Sano, S.; Nhuan, M.T. Uptake of metals and metalloids by plants growing in a lead-zinc mine area, Northern Vietnam. *J. Hazard. Mater.* **2011**, *186*, 1384-1391. [CrossRef]

65. Tremlová, J.; Vašíčková, I.; Száková, J.; Goessler, W.; Steiner, O.; Najmanová, J.; Horáková, T.; Tlustos, P. Arsenic compounds occurring in ruderal plant communities growing in arsenic contaminated soils. *Environ. Exp. Bot.* **2016**, *123*, 108-115. [CrossRef]

66. Kabata-Pendias, A.; Pendias, H. Trace Elements in Soils and Plants, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2001.

67. Krämer, U. Metal Hyperaccumulation in Plants. *Annu. Rev. Plant. Biol.* **2010**, *61*, 517-534. [CrossRef]

68. Ghaderian, Y.S.M.; Lyon, A.J.E.; Baker, A.J.M. Seedling mortality of metal hyperaccumulator plants resulting from damping off by Pythium spp. *New Phytol.* **2000**, *146*, 219-224. [CrossRef]

69. Mattina, M.I.; Lannucci-Berger, W.; Musante, C.; White, J.C. Concurrent plant uptake of heavy metals and persistent organic pollutants from soil. *Environ. Pollut.* **2003**, *124*, 375-378. [CrossRef]

70. Ali, H.; Khan, E.; Sajad, M.A. Phytoremediation of heavy metals—Concepts and applications. *Chemosphere* **2013**, *91*, 869-881. [CrossRef]

71. Sudmoon, R.; Neeratanaphan, L.; Thamsenanupap, P.; Tanee, T. Hyperaccumulation of Cadmium and DNA Changes in Popular Vegetable, Brassica chinensis L. *Int. J. Environ. Res.* **2015**, *9*, 433-438.

72. Bui, A.T.K.; Nguyen, H.T.H.; Nguyen, M.N.; Tran, T.-H.T.; Vu, T.V.; Nguyen, C.H.; Reynolds, H.L. Accumulation and potential health risks of cadmium, lead and arsenic in vegetables grown near mining sites in Northern Vietnam. *Environ. Monit. Assess.* **2016**, *188*, 525. [CrossRef] [PubMed]

73. Hesami, R.; Salimi, A.; Ghaderian, S.M. Lead, zinc, and cadmium uptake, accumulation, and phytoremediation by plants growing around Tang-e Douzan lead-zinc mine, Iran. *Environ. Sci. Pollut. Res.* **2018**, *25*, 8701-8714. [CrossRef] [PubMed]

74. Fait, S.; Fakhi, S.; Elmzibri, M.; Malek, O.A.; Rachdi, B.; Faiz, Z.; Fougrach, H.; Badri, W.; Smouni, A.; Efahr, M. Behavior of As, Cd, Co, Cr, Cu, Pb, Ni, and Zn at the soil/plant interface around an uncontrolled landfill (Casablanca, Morocco). *Remediat. J.* **2018**, *28*, 65-72. [CrossRef]

75. Van Der Ent, A.; Baker, A.J.M.; Reeves, R.D.; Pollard, A.J.; Schat, H. Hyperaccumulators of metal and metalloid trace elements: Facts and fiction. *Plant. Soil* **2012**, *362*, 319-334. [CrossRef]

76. Kabata-Pendias, A. Trace Elements in Soils and Plants; CRC Press: Boca Raton, FL, USA, 2000.

77. Reeves, R. The hyperaccumulation of nickel by serpentine plants. In *The Vegetation of Ultramafic (Serpentine) Soils*; Baker, A.J.M., Proctor, J., Reeves, R.D., Eds.; Intercept Ltd.: Andover, UK, 1992; pp. 253-277.

78. Brunetti, G.; Soler-Rovira, P.; Farrag, K.; Senesi, N. Tolerance and accumulation of heavy metals by wild plant species grown in contaminated soils in Apulia region, Southern Italy. *Plant. Soil* **2008**, *318*, 285-298. [CrossRef]

79. Pan, P.; Lei, M.; Qiao, P.; Zhou, G.; Wan, X.; Chen, T. Potential of indigenous plant species for phytoremediation of metal(loid)-contaminated soil in the Baoshan mining area, China. *Environ. Sci. Pollut. Res.* **2019**, *26*, 23583-23592. [CrossRef]

80. Feng, R.; Wei, C.; Tu, S.; Ding, Y.; Wang, R.; Guo, J. The uptake and detoxification of antimony by plants: A review. *Environ. Exp. Bot.* **2013**, *96*, 28-34. [CrossRef]

81. Krzesłowska, M.; Samardakiewicz, S.; Woźny, A. Trace Metals in: *The Reactions of Plant Cells to Stress Factors*; Woźny, A., Goździcka-Józefiak, A., Eds.; Wyd. Nauk. Uniwersytetu w Poznaniu, Poznań: Poznań, Poland, 2010; Volume 2, pp. 90-146.

82. Fahr, M.; Laplaze, L.; Bendaou, N.; Hocher, V.; El Mzibri, M.; Bogusz, D.; Smouni, A. Effect of lead on root growth. *Front. Plant. Sci.* **2013**, *4*, 175. [CrossRef]
83. Aslan, M.; Yılmaz, Y.Z.; Yılmaz, Y.Z.; Ünlü, M.Y.; Türkmen, N. Sorption of Cadmium and Effects on Growth, Protein Content, and Photosynthetic Pigment Composition of Nasturtium officinale R. Br. and Mentha aquatica L. *Bull. Environ. Contam. Toxicol.* **2003**, *71*, 323–329. [CrossRef]

84. Nedjimi, B. Seasonal growth and translocation of some major and trace elements in two Mediterranean grasses (Stipa tenacissima Loefl. ex L. and Lygeum spartum Loefl. ex L.). *Acta Oecologica* **2018**, *89*, 43–50. [CrossRef]

85. Baghdad, B.; Naimi, M.; Bouabdli, A.; Sonnet, P.; Lutts, S. Heavy Metals in Tailings, Soils and Vegetation of an Abandoned Lead Mine Land in Morocco. *Ecol. Mediterr.* **2006**, *32*, 85–91. [CrossRef]

86. Samkaeva, L.; Revin, V.; Revin, Y.; Kulagin, A.; Novikov, O.; Pugayev, S. Studying of heavy metal accumulation by plant. *Biotech* **2001**, *1*, 54–59.

87. Takeda, R.; Noriyoshi; Matsumoto, S.; Komemushi, S.; Sawabe, A. Accumulation of Heavy Metals by Japanese Weeds and Their Seasonal Movement. *Contam. Soils, Sediments Water* **2006**, *349–359*. [CrossRef]

88. Alirzayeva, E.G.; Shirvani, T.S.; Alverdiyeva, S.M.; Shukurov, E.S.; Öztürk, L.; Ali-zade, V.M.; Çakmak, İ. Heavy metal accumulation in Artemisia and foliaceous lichen species from the Azerbaijan flora. *For. Snow Landsc. Res.* **2006**, *80*, 339–348.

89. Ashraf, M.; Hayat, M.Q.; Mumtaz, A.S. A study on elemental contents of medicinally important species of Artemisia L. (Asteraceae) found in Pakistan. *J. Med. Plants Res.* **2010**, *4*, 2256–2263.

90. Yanqun, Z.; Li, Y.; Schvartz, C.; Langlade, L.; Fan, L. Accumulation of Pb, Cd, Cu and Zn in plants and hyperaccumulator choice in Lanping lead-zinc mine area, China. *Environ. Int.* **2004**, *30*, 567–576. [CrossRef]

91. Liu, X.; Gao, Y.; Khan, S.; Duan, G.; Chen, A.; Ling, L.; Zhao, L.; Liu, Z.; Wu, X. Accumulation of Pb, Cu, and Zn in native plants growing on contaminated sites and their potential accumulation capacity in Heqing, Yunnan. *J. Environ. Sci.* **2008**, *20*, 1469–1474. [CrossRef]

92. Wójcik, M.; Sugier, P.; Siebielec, G. Metal accumulation strategies in plants spontaneously inhabiting Zn-Pb waste deposits. *Sci. Total Environ.* **2014**, *487*, 313–322. [CrossRef]

93. Chaabani, S.; Abdelmalek-Babbou, C.; Ben Ahmed, H.; Chaabani, A.; Sebei, A. Phytoremediation assessment of native plants growing on Pb-Zn mine site in Northern Tunisia. *Environ. Earth Sci.* **2017**, *76*, 585. [CrossRef]

94. El Mamoun, I.; Fahr, M.; Mohammed, A.; Najib, B.; Abidine, T.Z.-E.; Abdelkarim, G.; Didier, B.; Laurent, L.; Abdelaziz, S. Zinc, lead, and cadmium tolerance and accumulation in Cistus libanotis, Cistus albidus, and Cistus salviifolius: Perspectives on phytoremediation. *Remediat. J.* **2020**, *30*, 73–80. [CrossRef]

95. Laplaze, L.; Doumas, P.; Brhada, F.; Ater, M.; Smouni, A. Use of Cistus Libanotis to Clean Heavy Metals Containing Soils. PATENTSCOPE International Application No. PCT/EP2010/056449, 18 November 2010.

96. De La Fuente, V.; Rufo, L.; Rodríguez, N.; Amils, R.; Zuluaga, J. Metal Accumulation Screening of the Río Tinto Flora (Huelva, Spain). *Biol. Trace Element Res.* **2009**, *134*, 318–341. [CrossRef] [PubMed]

97. Abreu, M.M.; Santos, E.S.; Magalhães, M.C.F.; Fernandes, E. Trace elements tolerance, accumulation and translocation in Cistus populifolius, Cistus salviifolius and their hybrid growing in polymetallic contaminated mine areas. *J. Geochem. Explor. 2012*, *123*, 52–60. [CrossRef]

98. Santos, E.S.; Abreu, M.M.; Nabais, C.; Magalhães, M.C.F. Trace element distribution in soils developed on gossan mine wastes and Cistus ladanifer L. tolerance and bioaccumulation. *J. Geochem. Explor. 2012*, *123*, 45–51. [CrossRef]

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).