Remote Sensing and Mycorrhizal-Assisted Phytoremediation for the Management of Mining Waste: Opportunities and Challenges to Raw Materials Supply

: In recent times, the development of innovative processes permits the application of a circular economy approach to the management and exploitation of mining waste with respect to human health and environment, such that society is changing its fundamentally negative perception of the mining sector. This study presents the opportunities and challenges of supplying raw materials from waste using a remote sensing technique, mycorrhizal-assisted phytoremediation, and hydrometallurgical techniques to transform mining waste from a problem to a resource. Soil/mine wastes from the Sierra Pintada mine (Mendoza, Argentina) were mineralogically and chemically analyzed, and then, a mapping of the mining waste was carried out by Sentinel-2A images to identify areas with similar characteristics. The bioaccumulation of HMs by autochthonous shrubs was also determined to select accumulator plant species, and to evaluate their potential for phytoremediation of mine soils at different technological scales, when they were inoculated with arbuscular mycorrhizal fungi originated from a mining-impacted area. RMs were recovered from plant biomass by scaling in bioreactors, the depuration module, and hydrometallurgical techniques. The encouraging results highlight that this multidisciplinary approach can be applied to meet the increasing demand for RMs supply and, at the same time, to protect the environment and public health.


Introduction
Since prehistoric times, humans have exploited their territory and its resources to their own ends (progress, competition, and economic growth) with little or no attention to the environmental, health, and social consequences. The worldwide population growth, the emerging technological development, and the high demand for raw materials (RMs) have led to an ever-less-sustainable consumption of mineral resources, since they were considered as abundant and easily extractable. During industrial processes and consumption, large quantities of waste are accumulated causing notable environmental problems associated with their disposal, including ecological losses, downstream contamination, and pronounced landscape transformation (e.g., stockpiled waste rock and tailings, subsidence basins, open pits, and removal of overburden rock and topsoil) [1]. After the Aznalcollar (Spain, 1998) and Baia Borsa/Baia Mare (Romania, 2000) accidents involving mining activities (the latest in a long list), in 2006 the European Commission (EC) adopted the Mining Waste Directive on the management of waste from extractive industries [2].
Although European countries have a long mining history, they produce only 3% of the world's ore, but consume 25%-30% of the world's metal [3]. Over the last decades, due to increasing domestic production costs, competitive land use, and resource depletion, EU countries' RMs import rate has dramatically increased, representing one of the key challenges for the EU economy. Nevertheless, to address this major challenge, in 2008, the EC adopted the European RMs Initiative [4] and, in 2011, a strategy document [5] based on three pillars to secure and improve access to RMs for the EU: 1) ensuring the fair and sustainable supply of RMs from international markets; 2) fostering sustainable supply within the EU; and 3) boosting resource efficiency and promoting recycling. In 2011, the EC started publishing the lists of Critical Raw Materials (CRMs), defined as those RMs that are economically and strategically important for the European economy [5][6][7][8]. CRMs are identified in the EU according to two main dimensions-economic importance and supply risk-that change over time and thus the list is subject to regular reviews and updates [9], the latest dating back to 2020 (https://rmis.jrc.ec.europa.eu/, accessed on 15 March 2023).
Current EU goals are to: (i) reduce import dependency by improving access to RMs within the EU and from other sources; (ii) promote resource efficiency also through recycling and substitution; and (iii) to reduce environmental and social impacts [9]. In this context, mining waste plays a major role as these huge deposits represent, at the same time, a social and environmental concern as well as potential anthropogenic resources, exploitable through appropriate management strategies and innovative processing solutions. However, to achieve that, it is necessary to know accurate information about waste because, at the time of opening mining activity, the waste chemical and mineralogical composition was ignored, or this material was considered below the cut-off grade [10]. In the present, the technological progress and development of innovative processes permit the application of a circular economy approach to the management and exploitation of mining waste with respect to human health and environment, so that society changes its fundamentally negative perception of the mining sector [10]. During the pandemic situation linked to COVID-19 and the Ukraine war, the concern has transformed into a major problem creating a status of dependence on extra-EU supply for Europe and risks of market distortions. To overcome this problem, multidisciplinary approaches are being developed that consider the common interest towards the raw materials needed, acting sustainably and responsibly for the green recovery and for building a fair and climate-neutral economy.
To this end, in the last few decades, a remote sensing approach has been increasingly used in geological applications (such as mineral exploration) by exploiting the specific interactions of geological materials with electromagnetic spectrum in three wavelength regions: visible to near infrared (VNIR-from 400 nm to 1000 nm), shortwave infrared (SWIR-from 1000 nm to 2500), and thermal infrared (TIR-from 7600 nm to 14,000 nm) [11,12]. Thanks to the diagnostic spectral absorptions, unique for each material, the remote sensing represents a key tool to individuate and map the mining waste in a valuable, fast, low-cost, and sustainable way, and, at the same time, to manage wastes rich in RMs as resources for the future [11][12][13]. In particular, Sentinel-2, the European optical imaging satellite launched in 2015 as part of the European Space Agency's Copernicus Program, is used in studies to optimize the management of mining waste [14], to improve mineral estimation by band ratio [15], to map mining residues rich in RMs [16], exploiting the spatial resolution of 10-20 m, and to revisit the frequencies of 10 days and spectral resolution ranging from 490 nm to 2190 nm.
The use of biological methods to improve the quality of the environment, degraded because of anthropo-pressure, is a priority at this time [17,18]. Many studies tackle the challenge to recover RMs from mining waste by using mycorrhizal-assisted phytoremediation (MAP) [19]. The MAP is composed of a hyperaccumulating plant species associated to an Arbuscular Mycorrhizal (AM) fungus in its roots, and the addition of an enzymatic cofactor, placed in an environment constituted by volcanic ashes and soil (methodology according to the granted patent) [20]. There are many hyperaccumulating species of heavy metal(loids) (HMs) described in the literature within the families Asteraceae, Poaceae, and Brassicaceae, and most of them establish mutualistic symbiosis with AM fungi (Glomeromycota Phylum). AM fungi are capable of uptaking nutrients and HMs from soil and translocating them to the host plant. They also accumulate these elements in the intraradical structures or immobilize them in external mycelium and spores in soil [21]. Some mycorrhized plant species show high bioaccumulation of HM in their roots, and others display high translocation and bioaccumulation in the aerial parts, with exclusion of HM from the reproductive and fruiting parts [19]. These behaviors highly depend on the plant species, the nature and concentration of the HM, the environmental and physical-chemical variables, such as the presence of other metals, and pH and Eh values in soil/substrate. The substances involved in the processes, some of them being described as the results of mycorrhizal symbiosis, are fundamentally phytochelatins (involved in HM chelation and making them more bioavailable for entry into cells), metallothioneins (that favor the entry of certain HM into cells), and glomalin (glycoproteins secreted by AM fungi that sequester HM in soils). Other transporting substances are still being studied [19,22,23].
The MAP system to be applied to the territory must go through different technology maturation stages called Technology Readiness Levels (TRLs). The TRL for pilot tests or prototypes that involve environmental engineering has nine levels, from the verification of conceptual hypotheses to the replication in territory [22,24,25]. Scotti et al. [24] proposed a Vegetable Depuration Module (VDM) for scaling up the MAP at TRL 6. The VDM falls under the definition of artificial subsurface wetland, but with the difference that it is modified in its solid matrix and in its slope with an outlet to a collecting chamber. These modifications allow us to know the partition of HMs and the possibility of recycling the filtered solution by setting other modules in series or parallel order with a change in calibration parameters [22,25].
In this work, we acquired a Sentinel-2A image to analyze a simultaneous image with the same time campaigns and we built a bioreactor (BR) that simulates the VDM at a smaller scale (TRL 4). The BR is ten times smaller than the VDM longitudinally but with the same characteristics of height, slope, and fill contents. The calibration of physicalchemical-biological variables in BR permits us to extrapolate and to verify results to the VDM (TRL 6). Additionally, any variable can be adjusted for the change in scale, and a protocol emerges to be applied in territory to the extent required by the project engineering [19,24,25]. The extraction potential of the system is related to the physical-chemicalbiological variables. The hydraulic variables are related to irrigation, inflow, outflow, slope, type of vertical or horizontal flow, hydraulic retention time, and hydraulic constant. Any variation of those conditions can promote the element exclusion or bioaccumulation in plant biomass [26]. The amount of biomass obtained and the high bioextraction potential with the MAP system are essential for the recovery of secondary raw materials (SRMs) and CRMs, this is why the scaling in VDM to TRL6 is proposed.
The aims of this study were to determine the potential of Sentinel-2 for the classification and mapping of mining waste, identifying the areas in the mine with the same mineralogical, chemical, and spectral characteristics, and then to apply the MAP system followed by the hydrometallurgy technique to extract HMs of importance, in order to transform mining waste from a problem to a resource, together with the possible opportunities and challenges of exploiting wastes.

Study Area and Sampling Sites
The Sierra Pintada Uranium mine, with an area of 2007 Ha, is located approximately 12 km to the Southwest of Villa 25 de Mayo in Mendoza province ( Figure 1). It was discovered in 1968 through aerial radiometric prospecting, identifying a series of anomalies with subsequent corroboration of Uranium presence within the area. The area belongs to the so-called San Rafael Block, for the Sierra Pintada area. It is constituted by the basement of ancient lands (Devonian and Carboniferous), affected by intense tectonism and a profuse volcanic activity developed in the Lower Permian that led to the accumulation of an important sequence of volcanic clastic rocks [27]. The large volumes of intermediate and acid ignimbrites, as consequences of permo-triassic volcanism, generated favorable conditions for the concentration of Uranium deposits [28]. From 1975, the Sierra Pintada Manufacturing Mining Complex (SPMMC) was put into operation for obtaining Ammonium Diuranate as a final product, thus converting it as the largest Argentine uranium deposit to date. However, in 1997, the mining complex was forced to stop its activities due to the drop in the international price of commodities.
In  Figure 1). About 1 Kg of soil sample was collected for further mineralogical and chemical characterization. Plant material was chemically determined to determine potential hyperaccumulating plant species. The geographical coordinates, pictures, and a brief description of each sampling site were registered.

Chemical and Mineralogical Analysis of Mining Waste
In the Soil Laboratory of the Uranium Mining Environmental Remediation Program (PRAMU) of the Nuclear Chemistry Department (CNEA), soil subsamples (about 50 g) were used for the determination of pH (1:3) with a Pocket Hanna Combo digital pHmeter, and the heavy metal(loids) (HMs) concentration by Neutron Activation Analysis (NAA). Previously, the samples were dried in an oven for 24 h and packed in quartz ampoules for irradiation together with the certified reference materials (NIST SRM 1633b Coal Fly Ash and Chinese Soil GBW 07405). Subsequently, the samples and the reference materials were irradiated for 4.5 h in the RA-3 reactor. Two measurements were carried out with decay times of approximately 7 and 30 days, and the following HMs were determined: As, Ce, Co, Cr, Cs, Eu, Fe, Hf, La, Nd, Rb, Sb, Sc, Sm, Sr, Th, and U. The Ortec GeHP detectors (30% efficiency, 1.9 keV resolution for the 60Co 1332.5 keV peak), the Ortec 919E multichannel buffer module and Gamma Vision data acquisition software was used. The values of concentrations were estimated through programs developed in the laboratory. Determinations were performed using neutron activation analysis (NAA) by the Nuclear Analytical Techniques Laboratory (TAN, CNEA).
The mineralogy determination was carried out by the Institute of Environmental Geology and Geoengineering (National Research Council, Montelibretti, Italy). For that, pulverized mining wastes were analyzed with the X-Ray Powder Diffraction technique (XRPD) by using a Bruker D2 Phaser XRD, equipped with Lynxeye detector, Cu-Ka radiation and working at 30 kV and 10 mA, a fast instrument for determining the phase content of polycrystalline material. Each sample exhibited a distinctive XRPD pattern and semi-quantitative phase composition was obtained by using Crystallography Open Database of the DIFFRAC.EVA software.

Chemical Analysis of Plant Samples
Root samples and aerial parts of individual plants of Senecio subulatus, Bacharis salicifolia, and Sena aphylla were collected in each sampling site and the content of HMs was determined through NAA as previously described. Firstly, plant samples were ground using a cryogenic mill and lyophilized for 24 h, then, they were pelletized and packaged for irradiation, and subsequently measured. Two measurements were made with decay times of approximately 7 and 30 days for the determination of: As, Ce, Co, Cr, Cs, Eu, Fe, Hf, La, Rb, Sc, Sm, Sr, Th, U, and Zn.

Analysis of Sentinel-2A Image in the Sierra Pintada Mine Area
The spatial, spectral, and temporal resolution of satellite data together with data availability, data storage, and data handling limitations, have prevented the use of remote sensing in many potential applications. To try and overcome these limitations, the European Space Agency (ESA, Paris, France) and European Commission have developed a series of next-generation Earth Observation missions, called Sentinels, for land, ocean, and atmospheric monitoring, and each Sentinel mission is based on a constellation of satellites to provide accurate, timely, and free datasets. In this study, the multispectral Sentinel-2A image was downloaded from the ESA Copernicus Open Access hub and acquired in the Sierra Pintada mine area on 12 December 2018. The spatial resolution of the Sentinel-2A image is of 10 m in the visible and near-infrared range and 20 m in the short-wave infrared range, with spectral resolution ranging from 490 nm to 2190 nm. The 13 georeferenced samples, collected during the in situ sampling campaign, were overlaid on the satellite image and, for each sample, the spectral signature was acquired. Then, all the spectral signatures and mineralogical data were used as input for cluster analysis to obtain groups of spectrally and mineralogically similar samples. In particular, the average linkage between groups criterion (UPGMA algorithm) and Euclidean similarity index have been used. Afterwards, the supervised pixel-based classification was applied using the Spectral Angle Mapper (SAM) algorithm (see Guglietta et al. (2020) [14]).

Potential Phytoremediation of Autochthonous Plants
Plants of Senecio subulatus, Bacharis salicifolia, and Sena aphylla species were found in all sampling sites. These plant species are autochthonous from the study area, and they are characterized by their rapid growth and their great capacity to colonize disturbed areas by mining activity. Furthermore, considering the impact level and the possibility of quantifying the concentration of some HMs in plant biomass of these species, it was decided to use them as the MAP system.
The bioconcentration factor (BCF) and translocation factor (TF) were used to evaluate the potential phytoremediation of plants. The BCF is the concentration of HMs in total biomass (ppm)/concentration HMs in soil (ppm); TF is the concentration HMs in aerial part (ppm)/concentration HMs in roots (ppm). A BCF value higher than 1 indicated that a plant was an accumulator. TF value determines plant efficiency in HMs translocation from the root to the shoot. A plant was considered efficient in metal translocation from root to shoot when TF was higher than 1 due to an efficient metal transport system. TF values less than 1 indicated ineffective metal transfer and suggested that these plants accumulate metals in the roots and rhizomes more than in shoots or leaves. The BCF and TF were determined for the vegetation that naturally grew in the SPMMC and the BCF of plant biomass in the BRs assay was also determined.

Mycorrhizal-Assisted Phytoremediation (MAP) in Bioreactors (BR) and Vegetable Depuration Module (VDM)
For the MAP system, the following plant species were considered: Senecio subulatus and Bacharis salicifolia plant species. They were inoculated with a mixed inoculum of mycelium, spores, and colonized root fragments composed by several AM fungi, previously isolated from Paramillos de Uspallata mine (Mendoza, Argentina). The AM fungal inoculum was obtained from a trap culture using Medicago sativa (L.) as host plant, isolated and propagated by Banco de Glomeromycota in vitro (FCEN UBA, https://bgiv.com.ar/, accessed on 12 May 2010) as described in [29].
Four BRs, placed in the Bioenvironmental Laboratory (CNEA-FRSR, Mendoza), were made with polycarbonate containers (30 cm × 60 cm and 80 cm in height and with 10 cm for the collector chamber to measure the volume/flow of the leachate generated by the substrate irrigation (Figure 2a)). Each BR was filled with three layers of stones of different granulometry. The first deeper layer (thickness of 10 cm) is represented by coarse gravel (average diameter = 10 cm); the second one (thickness of 15 cm) by medium gravel (average diameter = 5 cm); the third layer (thickness of 20 cm) by fine gravel (average diameter = 1 cm); and the upper layer (thickness of 15 cm) by mine soil samples, volcanic ashes (1:1, v/v), and 350 ppm ZnSO4 (contaminated soil, CS) (Figure 2a). The control treatment was performed with commercial soil (BS). The BR and test chamber were connected to a collection chamber (10 cm × 30 cm and 80 cm deep) to collect percolated water. These BR had a relation of 1:10 with the vegetable depuration module (VDM, TRL 6), and conserved the dimensions, slope, filter layers, physicochemical (pH, Eh), and hydraulic parameters (hydraulic constant, hydraulic retention time, and vertical income flow) [19]. The depth and slope of the VDM and the BR are 80 cm and 6%, respectively. The VDM is located at the Bioenvironmental Laboratory (CNEA-FRSR, Mendoza) and a description of its performance is detailed in Scotti et al. [19] (Figure 2c). The size of VDM is 3 m in width and 5 m in length with a collection chamber of 3 m in width and 1 m in length (Figure 2b). According to Scotti et al. [25], the fill is made up of layers of coarse gravel, fine gravel, pellets, and the substrate to be bioremediated.

Hydraulic Parameters in BR and VDM
To model the hydraulic flow in the BR and the VDM, Darcy's law was used. This law considers the inlet and outlet flow when field capacity is reached by this Equation (1): where Ks: hydraulic constant (m/days); Q: average flow rate at which it enters and which it exits (m 3 /days); Ac: perpendicular area to the flow (m 2 ) = 0.3 m × 0.5 m; s: slope (m/m) = 0.03 m/0.5 m. Furthermore, the hydraulic retention time (th) was measured by registering the time it took the flow to cross down the different layers and to exit towards the collecting chamber.
Additionally, the water-holding capacity (WHC) of the BR was determined by Equation (2): where MWHC: mass of the water in grams; Mt: total mass of the container and wet soil in grams; Ms: total mass of the container and dry soil in grams.

Experimental Design and Measurement of Parameters
Four BR were performed, as described above, under two treatments: control soil (BS) and contaminated mine soil substrate (CS), each one inoculated by AM fungi and uninoculated. In each BR, five plants of Senecio subulatus and eight plants of Bacharis salicifolia were grown. One of the BR was inoculated with the AM fungal inocula and the other one was not inoculated. After 4 months, plants were harvested from each BR. The shoots were separated into leaves and flowers, and the roots were carefully rinsed with distilled water to eliminate substrate particles. Fresh weights of leaves, flowers, and roots were registered, then dried in an oven at 80 °C for 72 h to obtain dry biomass.
The BCF and TF were determined as previously described for the following heavy metal(loids): Cr, Zn, P, Ni, Rb, Sr, Cu, and Mn [17]. Furthermore, the BioPotential Extraction (BPBR/VDM (g)) of the system, defined as the HMs concentration in total biomass (ppm) * density of biomass plants in Br and VDM (g), was calculated.

Map of Characterized Mining Waste
The dendrogram of cluster analysis has identified three main groups (Cluster 1, Cluster 2, and Cluster 3) of mining waste and three samples (S3P2, S5P1, and S5P2) that do not belong to any cluster with different spectral characteristics and mineralogical phases (Figures 3 and 4). Each cluster represents a distinct class with the samples used, respectively, as training set to classify the Sentinel-2A image and as ground control point to validate the results of SAM ( Figure 5). In particular, the results of accuracy include the user's accuracy (84.13%), producer's accuracy (89.6%), overall accuracy (90.73%), and kappa coefficient (0.74).  The results of SAM classification have been used to create a map of mining waste that highlights the location of the classes in the Sierra Pintada mine with the major mineralogical phases ( Figure 5). Class A is represented by samples from Cluster 1 (S4P1, S4P2, and S6P2) with 32.89% of andesine, 39.15% of quartz, 9.89% of muscovite, 3.96 % calcite, and 2.47% chlorite. In particular, the sample S4P1 belonging to this class derives from the processing activities of the mine, a stand with more than 20 years, with the natural native vegetation that is growing all around and the wind and water erosion processes that are accumulating material above the processing mine waste. Class B is characterized by samples from Cluster 2 (S2P3, S2P2, and S1P2) with 62% of andesine, 24.56% of quartz, 6.47% of muscovite, and 5.65% of hornblende, and this class includes the samples collected in the natural area along the Tigre River. The class C is characterized by samples from Cluster 3 (S3P1, S6P1, S2P1, and S1P1) with natural recolonization of the area by native vegetation. The main minerals are represented by quartz (70.24%), chlorite (11.97%), muscovite (10%), dolomite (4.29%), and calcite (3.52%). Furthermore, S3P2, S5P1, and S5P2, not belonging to the other three clusters, are thus different classes characterized by different mineralogical composition. S3P2 contains waste from the extraction process with 77.48% of quartz, 10.74% of chlorite, and 12.79% of muscovite. S5P1 and S5P2 have been collected in the mine area with waste from the processing mine activity accumulated over the time and with different main mineralogical phases: S5P1 with exclusively of 100% gypsum and S5P2 with quartz (84.14%) and albite (8.72%). The map has been used to identify the different kinds of waste that are selected for further phytoremediation analysis: S4P1 from Class A (in yellow), S1P2 from Class B (in red), S2P1, S3P1, and S6P1 from Class C (in green), S5P1 (in cyan), and S5P2 (in magenta).

Potential Natural Phytoremediation
The results in Tables 1 and 2 represent the BCF and TF values for sampled vegetation grown under natural conditions at each site from the SPMMC. A value of BCF higher than 1 or close to 1 (in bold) indicates that the plant accumulated HMs in their biomass. These results support previous reports about members of the Asteraceae family being capable of tolerating and accumulating large amounts of HMs and radionuclides in their tissues, and usually establishing the AM symbiosis in their roots [30].  Regarding the TF, the results in Table 2 indicate that the HMs passed from the roots to the aerial part when the value of TF is greater than 1, whilst a TF < 1 indicates low translocation to the aerial part, and possibly the stabilization of HMs in the roots or soil. The highest value of TF for Co, Eu, and Sr were found in the biomass of Senecio subulatus grown at the S4P1 site (Class A) with 14.35, 3.78, and 4.28, respectively. The TF values for Bacharis salicifolia plants found at S1P2 site (Class B) were 3.78 for Fe and 4.05 for Sm. With regards to the ability of Sena aphylla to translocate HMs from roots to aerial parts (Table 2), it was observed that the highest values of TF were obtained for La (12.27), Sm (47.14), Cr (3.5), Eu (1.74), and Sr (1.17) in plants grown at the S2P1 and S6P1 sites (Class C), Zn (4.0) from S5P2, and Fe (1.25) from the S4P1 (Class A) sites. For the rest of the elements, low values of TF were obtained.
In Sena aphyla plants from the S5P1 site, the TF values for As and Th were close to 1, and these elements could be considered as transferable in Sena aphylla. HMs with very low values of TF, such as U and Th (S5P2 site), Co (S5P1 site), and Zn (S2P1 site), could be stabilized in the root tissues or soil. Although the TF indicates translocation to the aerial part of Sena aphylla for some elements, the low BCF values indicated a weak translocation of these elements from the soil to roots, or probably a strong sequestration of HMs into the soil matrix. The efficiency of Sena aphylla to uptake the elements could depend on root exudates, soil pH, texture, and content of organic matter. Otherwise, soil microorganisms and AM fungi could greatly influence the sequestration and/or assimilation processes of Sena aphylla plants [17].
A positive correlation between the concentrations of HMs in soil with the concentrations in roots and aerial parts was detected in Sena aphylla plants grown at the S5P1, S5P2, and S1P2 sites (Cluster 2) (Figures 3 and 5). In the remaining sites, there was a positive correlation between the concentration of the elements in the root and in the aerial part.
The HMs concentrations in the roots of Sena aphylla were higher than the values in the aerial parts for Zn, Cs, La, Sc, Rb, Co, Cr, Eu, Fe, Sr, and Sm, which are transferred from the root to the aerial part while others (such as U, Th, and Co) did not, and this depended on the low value of TF (TF < 1). The results indicate that these elements are probably stabilized in the root system. Senecio subulatus and Bacharis salicifolia plants were able to translocate Sr from roots to the aerial part. This is of great importance given that Sr is considered as an RM.
According to the Total BCF (Table 2) of the studied vegetation developed under natural conditions in the SPMMC, Senecio subulatus, Sena aphylla, and Bacharis salicifolia had the capacity to accumulate great amounts of Zn and Sr. Figure 6 shows that the presence of Zn is only for the S4P1 (Class A), S5P1, S5P2, and S6P1 (Class C) sites. In particular, the BCF for Zn reached 47.5 in Bacharis salicifolia plants at the S4P1 site (Class A), and the lowest value in Senecio subulatus at the S5P1 site (BCF = 11.3), while Sena aphylla at S6P1 (Class C) resulted in the highest value (BCF = 140).
It is noteworthy that the S4P1 (Class A) and S6P1 (Class C) sites with plants with high total BCF values of Zn (47.5 and 140, respectively), had slightly alkaline soil pH values   Figure 7 showed the total BCF for Sr in plant biomass of the three sampled plant species. It is important to note that the value of the total BCF for Sr in Sena aphylla biomass was higher than the other plant species. In particular, the BCF (1.35) at the S6P1 site (Class C) was higher compared to the other plant species within the same site, and the other mine sites. Furthermore, the BCF values of Sr at the S4P1 site (Class A) were equal to 0.41 in Sena aphylla, 0.31 in Bacharis salicifolia, and 0.29 in Senecio subulatus. At the S1P2 site (Class B), they were similar to the previous class (Class A), even if the samples were from different sites within the mining area. Class A represents mining waste from processing activities but with accumulating materials (i.e., natural rocks) on its surface by erosion processes, while Class B is characterized by samples from natural areas. The lowest values of BCF (≤0.22) for Sr were highlighted in Senecio subulatus and Bacharis salicifolia plants grown at the S1P2 (Class B), S5P1, and S6P1 (Class C) mine sites.
The O.M percentages were lower at S4P1 (Class A) site (O.M = 0.3%) and higher at the S5P1 site (O.M = 1.5%). In particular, the S4P1 (Class A) and S5P1 represent areas impacted by mining activity but with different re-vegetation processes due to site management and topographic conditions. The S4P1 site is characterized by a flat topography with a high evaporation rate, while the S5P1 is a terrain that is more uneven, with a high accumulation of soil moisture. The natural conditions of the other sites showed the typical percentages of O.M of these semi-arid zones [27]. The cluster analysis of the total BCF of HMs, the pH values, and O.M contents in soil determined clusters according to the degree of disturbance or activity related to the substrate, but also to the presence of the different plant species (Figure 8). Despite the S4P1 site being characterized by a high level of disturbance in the substrate compared to the S5P1, S5P2, and S6P1 sites, they were grouped by the HMs assimilation capacity in Sena aphylla and Senecio subulatus plants. Similar results were found with Senecio subulatus and Sena aphylla in the S5P1 site, although there was a distance between S. aphylla in comparison to Senecio subulatus in the S6P1 site (low level of disturbance). Furthermore, B. salicifolia in the S4P1 site was distant compared to in the S1P2, S2P1, and S3P1 sites (low level of disturbance). These results could be explained by the typical ecological behaviors of the colonizing plant species since they usually depend on adapted mechanisms (such as competition) to establish in areas with few resources.

MAP Application in BR and Estimation of the Bioextraction Potential (BP) in VDM
Bacharis salicifolia and Senecio subulatus plant species were chosen for the BR experiment to recover Sr and Zn from plant biomass. Then, the results obtained in the BR were projected in the VDM. Sena aphylla could not be scaled in BR since it did not survive the nursery stage, even though the percentage of germinated seeds was good (300 seeds 90% germination). Bacharis salicifolia and Senecio subulatus were placed in the BR containing mining waste from the S4P1 (Class A) site (contaminated soil, CS treatment) and the other plants in BR with the control soil (blank soil, BS), under the protocol previously described.  Table 3. The inoculated plants of Bacharis salicifolia had the highest aerial biomass. Biological parameters such as the type of fungal inoculum MA, the plant biomass, and the physical-chemical and hydraulic parameters, such as hydraulic constant (Kh), water holding capacity (WHC), inlet flow, pH, Eh, and enzymatic cofactors, are important in carrying out scaling in the VDM and in retrieving items of interest, such as RMs. It is necessary to establish the experimental conditions of the MAP system in the BR to later scale it in VDM. In this sense, the protocol to be applied was: inflow every 2 days at field capacity, with no outlet to the collection chamber, vertical flow, mycorrhizal tension in the impacted area of the Paramillos mine, and density of 9 plants per BR. The Kh value was 165 ± 36 m 3 /d/m 2 (mean ± SD), and the WHC resulted in 5.4 kg. These results will be considered for scaling-up in the VDM.
In previous experiments, we observed that the mentioned physical-chemical-biological variables were relevant to predict possible competitions between HMs at the membrane-cell level that could modify the patterns of entry of elements, or the presence of fungal-bacterial consortia that could exude chemical substances that alter the solubility of the different HMs present in soils [17][18][19]. To verify the statistical significance of the biomass obtained from the plant species in the BRs, a Tukey Test was performed (Table  S1). The post hoc range test identified homogeneous subsets of means that did not differ from each other. No significant differences in the root biomass of B. salicifolia in the BS and inoculated (M+) and CS treatments were observed. On the contrary, the aerial parts of the plants inoculated with AM fungi (M+) of the BS and CS treatments presented more biomass than the non-inoculated one. On the other hand, the inoculated plants (M+) of S. subulatus had a higher biomass than the uninoculated plants of the BS treatment whilst, in contaminated soil, the plants did not survive under these conditions in the BRs. On the other hand, S. aphylla did not survive in the BRs.
The bioaccumulation of Sr and Zn of the studied plant species are shown in Supplemental Table 1; Sena aphylla was a good accumulator, followed by Bacharis salicifolia, and finally Senecio subulatus for Sr, whilst Bacharis salicifolia was the best bioaccumulator for Zn.
The bioextraction potential (BP) of the inoculated Bacharis salicifolia, previously analyzed under BR conditions, was scaled to the VDM with the following parameters: 12 plants per 0.15 m 2 , 2 m 3 of substrate (soil: volcanic ash) at field capacity conditions. The projection of BP in the VDM are shown in Table 4. Although the three plant species resulted in a high BCF for Zn and Sr elements, only Bacharis salicifolia displayed adequate characteristics (resistance to manipulation, biomass) for scaling-up in the mycorrhizal-assisted phytoremediation. Sr is considered a CRM by the fifth technical assessment in 2023 for the EU, and Zn, even if not considered critical, is important for the EU economy [3].
The MAP system was established with native plant and AM fungal species, both adapted to mine areas to efficiently extract HMs, including RMs. The results showed that the studied native species were much more efficient than non-native species associated with AM fungi from non-impacted areas. In Table 5, the BP values of Zn and Sr obtained in this experiment (678 g and 1576 g, respectively) far exceed those obtained in previous report (159 g and 30.2 g, respectively) [31]. These results may be due to the application of the AM fungal inoculum adapted to grow in metalliferous or contaminated soils that possibly show specific mechanisms for the uptake and tolerance of toxic elements [32]. Several authors have confirmed that the inoculation with AM fungi adapted to these soils increases the efficiency of bioremediation processes through greater absorption by the external mycelium and translocation to the plant, providing better protection against stress [32,33]. This result confirms the importance of AM fungal strains adapted to these conditions, for their subsequent application in degraded soils.

Application of Hydrometallurgical Process
There is an increasing need to develop innovative solutions, to implement environmentally sustainable practices in the recovery of SRMs and CRMs accumulated, and to concentrate on plant biomass by application of phytoremediation techniques [19,31]. Hydrometallurgical processes have a wide range of useful applications, in the mining and environmental sector for the remediation of polluted sites, for the recovery and purification of SRMs and CRMs from wastes, and for their commercial reuse as valuable metals [19,[34][35][36][37][38][39].
The results achieved in this experimental investigation encourage the subsequent recovery of SRMs and CRMs by hydrometallurgical techniques, with final purification by selective electrodeposition, as a viable environmentally friendly and cost-effective option. They allow, using specific innovative reagents and the application of electrochemical techniques, the selective recovery of precious valuable metals at high degree of purity. Therefore, following this investigation, a first phase by laboratory experiments will be carried out, assuming an integrated process circuit consisting of a leaching phase of accumulated and concentrated SRMs and CRMs, and subsequent purification by electrochemical methods. The main aim of the leaching phase will be to reach a recovery of 95% of purified metals, accumulated and concentrated on plant biomass, such as demonstrated in previous studies [19,31,37,38]. Regarding the purification phase by electrochemical methods, the scope will be to reach a purity level over 99% [37][38][39] that will permit an industrial reuse of SRMs and CRMs in the context of a circular economy approach. In a subsequent step, an integrated flowsheet will be drawn, with the aim of outlining a first scale-up of the process for industrial application.

Conclusions
Considering the actual worldwide economic and geopolitical situation, the technological advances in the field of remote sensing techniques (i.e., new sensors with high spatial and spectral resolution) and processing software, represent an opportunity to map the potential deposits of RMs, and then to supply these valuables resources from mining waste in a direct, cost-effective, and rapid way.
In this study, we have investigated the potential of Sentinel-2 to identify and classify the different types of mining waste considering similar mineralogical phases and spectral signatures, and then to create a map to recognize subjected sites for further phytoremediation and extraction of RMs from mine wastes. In previous works, we applied this methodology on mining wastes of an open mine, obtaining encouraging results. In this study, we applied this methodology to the Sierra Pintada mine, a stand with more than 20 years, with different mineralogical and spectral characteristics, in comparison with mines of previous case studies.
The inoculation of plants with AM fungal species, isolated from the Paramillos de Uspallata mine, significantly improved the efficiency of the phytoextraction and bioaccumulation of autochthonous plant species. The AM fungi were found in soils with high content of Zn and Sr within the Paramillos de Uspallata mine, thus this could be related to a positive effect on these RMs by improving their bioaccumulation in the biomass of mycorrhizal autochthonous plant species. Three autochthonous plant species capable of accumulating RMs under natural conditions were chosen to optimize the bioextraction potential in the BR, through the adjustment of physical-chemical and biological parameters, prior to scaling-up to the VDM for obtaining an abundant extractive biomass of Sr and Zn. Only two plants species (Baccharis salicifolia and Senecio subulatus) were tolerant of the contamination condition of the soils and survived, but only Bacharis salicifolia had a good development of biomass with high extraction of RMs elements. For further studies, conditions in the BR will be optimized so that all plant species can be tested.
Finally, the remote sensing technique will be applied by multitemporal analysis from different sensors (i.e., Landsat, Sentinel-2), to highlight the process of recolonization by native vegetation in the mine and to identify the mining area for the phytoremediation process.