Interpretation of the Reﬂectance Spectra of Lithium (Li) Minerals and Pegmatites: A Case Study for Mineralogical and Lithological Identiﬁcation in the Fregeneda-Almendra Area

: Reﬂectance spectroscopy has been used to identify several deposit types. However, applications concerning lithium (Li)-pegmatites are still scarce. Reﬂectance spectroscopic studies complemented by microscopic and geochemical studies were employed in the Fregeneda–Almendra (Spain–Portugal) pegmatite ﬁeld to analyze the spectral behavior of Li-minerals and ﬁeld lithologies. The spectral similarity of the target class (Li-pegmatites) with other elements was also evaluated. Lepidolite was discriminated from other white micas and the remaining Li-minerals. No diagnostic feature of petalite and spodumene was identiﬁed, since their spectral curves are dominated by clays. Their presence was corroborated (by complementary techniques) in petalite relics and completely replaced crystals, although the clay-related absorption depths decrease with Li content. This implies that clays can be used as pathﬁnders only in areas where argillic alteration is not prevalent. All sampled lithologies present similar water and/or hydroxide features. The overall mineral assemblage is very distinct, with lepidolite, cookeite, and orthoclase exclusively identiﬁed in Li-pegmatite (being these minerals crucial targets for Li-pegmatite discrimination in real-life applications), while chlorite and biotite can occur in the remaining lithologies. Satellite data can be used to discriminate Li-pegmatites due to distinct reﬂectance magnitude and mineral assemblages, higher absorptions depths, and distinct Al–OH wavelength position. The potential use of multi- and hyperspectral data was evaluated; the main limitations and advantages were discussed. These new insights on the spectral behavior of Li-minerals and pegmatites may aid in new Li-pegmatite discoveries around the world.


Introduction
Reflectance spectroscopy studies have shown a relevant economic potential to mineral exploration, either through the direct characterization and identification of the mineralogy of the deposits or through the identification of hydrothermal alteration zonation as a proxy to target mineralizations such as (i) porphyry copper deposits [1], (ii) volcanogenic massive sulfide (VMS) deposits [2,3], (iii) rare-earth element (REE) deposits [4,5], (iv) Figure 1. Location of the study area in the Iberian Peninsula and geological setting of the Fregeneda-Almendra pegmatite field (adapted from References [26,39,40]). The field stations where the samples were collected are marked with a circle for the training areas (the color is dependent on the sample type) and with a diamond for the areas misclassified as Li-pegmatite using machine learning algorithms. Li-mineral specimens were also collected in Li-pegmatite sites. Map projection: Universal Transverse Mercator zone 29N from the WGS84 datum.

Reflectance Spectroscopy Studies
The spectral studies were conducted in the laboratory employing two distinct spectroradiometers, both covering the spectral range between 350 and 2500 nm: (i) the ultrahigh-resolution SR-6500 (Spectral Evolution, Inc.) equipment (1.5 nm @ 700 nm, 3.0 nm @ 1500 nm, and 3.8 nm @ 2100 nm) and (ii) the FieldSpec 4 (ASD, Inc.) standard resolution equipment (3.0 nm @ 700 nm, 10.0 nm @ 1400 nm, and 10.0 nm @ 2100 nm). In each case, the point spectral measurements were performed using a Contact Probe with similar characteristics, namely an internal light source provided by a halogen bulb and spot size of 10 mm. For reflectance calibration, a Spectralon (Labsphere) plate with a maximum reflectance above 95% for the 250-2500-nm range and above 99% for the 400-1500-nm region. Small (less than 2 cm) or friable samples were analyzed over the white plate. When using either of the spectroradiometers, each measurement corresponded to an average of 40 scans. For the FieldSpec 4 equipment, four additional spectra were acquired in each analyzed spot, later averaged into a final spectrum, to increase the signal-to-noise ratio.
ASD ViewSpec Pro™ software [41] and SpectraGryph software [42] were used to individually analyze all spectra to exclude measurement errors or spectra with noticeable noise caused by vegetation fragments on the surface of the samples and for spectra averaging. Point measurements belonging to the same sample and with similar spectral behavior were averaged to aid visualization and allow fast interpretation of the spectra. The final Li-minerals spectral library contains reference data from 28 samples (five spodumene, 13 petalite, and ten lepidolite samples). Moreover, to complete the spectral studies for lithological classification, 45 samples from possible training areas were analyzed (five from CXG metasediments, 33 from distinct Li-pegmatite types, and seven from granitoid  [26,39,40]). The field stations where the samples were collected are marked with a circle for the training areas (the color is dependent on the sample type) and with a diamond for the areas misclassified as Li-pegmatite using machine learning algorithms. Li-mineral specimens were also collected in Li-pegmatite sites. Map projection: Universal Transverse Mercator zone 29N from the WGS84 datum.

Reflectance Spectroscopy Studies
The spectral studies were conducted in the laboratory employing two distinct spectroradiometers, both covering the spectral range between 350 and 2500 nm: (i) the ultrahigh-resolution SR-6500 (Spectral Evolution, Inc.) equipment (1.5 nm @ 700 nm, 3.0 nm @ 1500 nm, and 3.8 nm @ 2100 nm) and (ii) the FieldSpec 4 (ASD, Inc.) standard resolution equipment (3.0 nm @ 700 nm, 10.0 nm @ 1400 nm, and 10.0 nm @ 2100 nm). In each case, the point spectral measurements were performed using a Contact Probe with similar characteristics, namely an internal light source provided by a halogen bulb and spot size of 10 mm. For reflectance calibration, a Spectralon (Labsphere) plate with a maximum reflectance above 95% for the 250-2500-nm range and above 99% for the 400-1500-nm region. Small (less than 2 cm) or friable samples were analyzed over the white plate. When using either of the spectroradiometers, each measurement corresponded to an average of 40 scans. For the FieldSpec 4 equipment, four additional spectra were acquired in each analyzed spot, later averaged into a final spectrum, to increase the signal-to-noise ratio.
ASD ViewSpec Pro™ software [41] and SpectraGryph software [42] were used to individually analyze all spectra to exclude measurement errors or spectra with noticeable noise caused by vegetation fragments on the surface of the samples and for spectra averaging. Point measurements belonging to the same sample and with similar spectral behavior were averaged to aid visualization and allow fast interpretation of the spectra. The final Li-minerals spectral library contains reference data from 28 samples (five spodumene, 13 petalite, and ten lepidolite samples). Moreover, to complete the spectral studies for lithological classification, 45 samples from possible training areas were analyzed (five from CXG metasediments, 33 from distinct Li-pegmatite types, and seven from granitoid rocks), Remote Sens. 2021, 13, 3688 5 of 33 as well as 11 samples from four areas misclassified as Li-pegmatite. All samples and point measurements were photographed.
The laboratory acquired spectra were normalized through continuum removal [43] to compare the spectra obtained by both spectroradiometers and enhance the spectral absorption features. In this case, a hull quotient was calculated after automatically fitting a convex-hull to each spectrum in a Python environment using the pysptools library [44]. The extraction of the wavelength, depth, and other geometric parameters of the absorption features was made using SpectraGryph software and the pysptools library. The complete spectral library and additional information on the equipment used and spectra processing can be found in Cardoso-Fernandes et al. [45,46].

Complementary Petrographic and Geochemical Studies
Microscopic studies were conducted at the University of Porto, using a transmitted light Nikon petrographic microscope (Eclipse E 400 POL), coupled with a ZEISS Axiocam camera for both the petalite alteration products and the samples from the four areas misclassified as Li-pegmatite. Selected petalite samples with distinct degrees of alteration were crushed to quantify their Li content at Activation Laboratories Ltd. (Actlabs), ON, Canada, using the Code 8-Lithium Ore analytical package. Sample preparation and crushing were conducted at the University of Porto using manual and automated mortar agate grinders (Retsch RM 200). Li quantification was achieved through inductively coupled plasma-optical emission spectrometry (ICP-OES) after four-acid digestion (detection limit 0.001%). Afterward, Li contents were converted to Li 2 O (%) using a compound ratio of 2.1527. Mineralogical information was obtained through powder X-ray diffraction (XRD) for friable samples and through optical microscopy and Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS) for hard rock samples. The XRD analysis was performed at the National Laboratory of Energy and Geology (LNEG, Portugal) using a Panalytical Diffractometer XPERT-PRO, equipped with a PW3050/60 goniometer, with a Theta-2Theta geometry and X-ray tube with a copper (Cu) anticathode. The tube voltage and current were set to 45 kV and 40 mA, respectively. Regarding the acquisition settings, a continuous scan mode with a scan range between 3 • and 90 • of 2teta, a step size of 0.0084 • , and 4.960 s of counting time per step was used. The SEM-EDS studies were conducted at the University of Porto, in previously carbon-coated thin sections, using a Hitachi FlexSEM 1000 SEM coupled to a Bruker QUANTAX EDS system.

Results
In this section, selected spectra from the spectral library are presented, taking into account the previously designated goals. The main mineralogy inferred from the spectra is presented, but a detailed analysis spectrum per spectrum can be found in the Supplementary Materials (Tables S1-S3).

Li-Minerals Spectral Database
The spectral database was built for three distinct Li-minerals: lepidolite ( Figure 2), spodumene (Figure 3), and petalite ( Figure 4). For this study, monomineralic samples were preferred to rock samples to avoid the interference of other minerals. As can be seen in Figure 2, all lepidolite spectra show the white mica typical features due to the hydroxide ion (OH − ) absorption around~1409 nm associated with the principal, sharp Al-OH absorption feature (between 2195 nm and 2205 nm) and the Al-OH secondaries (around~2348 nm and~2435 nm, respectively) [47][48][49]. Other white mica-related, weak absorption features located~1537 nm,~2002 nm, and~2114 nm are due to OH − and OH − combined with lattice modes [49,50]. The depth of the water feature (~1909 nm) decreases from Figure 2 (a-c), and the slight inflection around~1460 nm indicates that montmorillonite is also present, mixed with white mica features [1,49]. The remaining spectra (d-k) have higher purity, allowing to characterize the lepidolite due to: (i) a peak in the red region, caused by lepidolite purple color and (ii) multiple absorption features in the visible and near-infrared (VNIR) region, together with the drop off to extreme blue observed in the raw reflectance curves (left panel of Figure 2). These features are interpreted to be related to ferric iron (Fe 3+ ) and manganese (Mn 3+ ) incorporated into the structure of lepidolite [50]. Ferrous iron (Fe 2+ ) may cause the absorption at 720 nm ( Figure 2 (e)) [16,51].  Figure 2). These features are interpreted to be related to ferric iron (Fe 3+ ) and manganese (Mn 3+ ) incorporated into the structure of lepidolite [50]. Ferrous iron (Fe 2+ ) may cause the absorption at 720 nm ( Figure 2e) [16,51].

Figure 2.
Spectral signature of lepidolite (Li-bearing mica): reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Each spectrum has a database code [45] providing information about the type of sample (L-Li-bearing minerals), mineralogy (Lep-lepidolite), sample provenance (FE-Feli mine and AL-Almendra), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4). Dotted curves: fresh samples; full curves: samples with visible signs of weathering.
The reference spectral data of spodumene was acquired considering pure or spodumene intergrowth with quartz. Despite this, it is not expected that the associated quartz will obscure the spectral behavior of spodumene because the main diagnostic features of the former do not occur in the VNIR and shortwave infrared (SWIR) regions [47,52]. Figure 3 shows the main absorption features located around ~1414 nm, ~1910 nm, and ~2209 nm for spodumene + quartz samples and around ~1412 nm, ~1913 nm, and ~2205 nm for spodumene samples. The presence of molecular water is evidenced by the broader, asymmetric features approximately at ~1410 nm (OH − + H2O) and ~1910 nm (H2O) [48,52]. Molecular water associated with (i) higher depth of the water feature in comparison with the Al-OH fundamental absorption (around ~2200 nm) and (ii) a slight inflection around ~1467 nm are diagnostic of smectite minerals-in particular, of montmorillonite [1,48,49]. The shallow Al-OH secondaries at ~2350 nm and ~2437 nm for spodumene + quartz, and at ~2352 nm and ~2440 nm for spodumene, indicate the existence of illite mixed with montmorillonite and possibly kaolinite [48,49]. Iron incorporation (Fe 3+ and possibly Fe 2+ ) into both spodumene and spodumene + quartz samples could be identified due to: (i) the drop Figure 2. Spectral signature of lepidolite (Li-bearing mica): reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Each spectrum has a database code [45] providing information about the type of sample (L-Li-bearing minerals), mineralogy (Leplepidolite), sample provenance (FE-Feli mine and AL-Almendra), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4). Dotted curves: fresh samples; full curves: samples with visible signs of weathering.
The reference spectral data of spodumene was acquired considering pure or spodumene intergrowth with quartz. Despite this, it is not expected that the associated quartz will obscure the spectral behavior of spodumene because the main diagnostic features of the former do not occur in the VNIR and shortwave infrared (SWIR) regions [47,52]. Figure 3 shows the main absorption features located around~1414 nm,~1910 nm, and 2209 nm for spodumene + quartz samples and around~1412 nm,~1913 nm, and~2205 nm for spodumene samples. The presence of molecular water is evidenced by the broader, asymmetric features approximately at~1410 nm (OH − + H 2 O) and~1910 nm (H 2 O) [48,52]. Molecular water associated with (i) higher depth of the water feature in comparison with the Al-OH fundamental absorption (around~2200 nm) and (ii) a slight inflection around 1467 nm are diagnostic of smectite minerals-in particular, of montmorillonite [1,48,49]. The shallow Al-OH secondaries at~2350 nm and~2437 nm for spodumene + quartz, and at~2352 nm and~2440 nm for spodumene, indicate the existence of illite mixed with montmorillonite and possibly kaolinite [48,49]. Iron incorporation (Fe 3+ and possibly Fe 2+ ) into both spodumene and spodumene + quartz samples could be identified due to: (i) the drop off to the blue observed in the reflectance curves (left panel of Figure 3) [50] and Remote Sens. 2021, 13, 3688 7 of 33 (ii) a slight inflection around~2240 nm related with the Fe-OH bond [49] (Figure 3 (a-h)). Spectra h-k show additional features diagnostic of kaolinite, such as (i) a small inflection around~2165 nm and (ii) the presence of an Al-OH triplet at~2312,~2350, and 2380 nm [8,48,49,53]. off to the blue observed in the reflectance curves (left panel of Figure 3) [50] and (ii) a slight inflection around ~2240 nm related with the Fe-OH bond [49] (Figure 3a-h). Spectra h-k show additional features diagnostic of kaolinite, such as (i) a small inflection around ~2165 nm and (ii) the presence of an Al-OH triplet at ~2312, ~2350, and ~2380 nm [8,48,49,53]. Each spectrum has a database code [45] providing information about the type of sample (L-Li-bearing minerals), mineralogy (Spd-spodumene), sample provenance (ALB-Alberto mine and HD-Hinojosa de Duero), a three-digit numerical part representing the GPS code, a one-or twodigit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
Most petalite spectra ( Figure 4) are very similar to spodumene, showing the same main absorptions (Figure 4a-q). Montmorillonite is identified in all samples by the broad, asymmetric, and acute absorptions around ~1413 nm and ~1911 nm, associated with the Al-OH feature at ~2206 nm [48,49], shallower than the ~1911 nm one, and the slight inflection around ~1467 nm [1,49]. It is the sole mineral identified in spectra a-c, e, and o-q ( Figure 4). Others (d and g-n) show additional shallow Al-OH secondary absorptions at ~2356 nm and ~2445 nm, indicating the mixture with illite [48,49]. The presence of the Al-OH triplet and the slight inflection around ~2165 nm allowed the identification of kaolinite in the spectra d and j ( Figure 3) [8,48,49,53]. Figure 4l shows a sharp absorption at ~2312 nm and a shallow feature at ~2390 nm, indicating that another smectite group mineral may be present, probably saponite, with diagnostic Mg-OH features at ~2310 nm and ~2380 nm [48,49]. The last three spectra of Figure 4 show distinct mineral assemblages. In spectra r and s, besides the montmorillonite features (broad, asymmetric water absorptions), it was possible to identify: (i) a deep symmetric but broad main Al-OH feature, (ii) Figure 3. Spectral signature of spodumene (a Li-bearing pyroxene): reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Dotted curves: fresh samples; full curves: samples with visible signs of weathering. Each spectrum has a database code [45] providing information about the type of sample (L-Li-bearing minerals), mineralogy (Spd-spodumene), sample provenance (ALB-Alberto mine and HD-Hinojosa de Duero), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
Most petalite spectra ( Figure 4) are very similar to spodumene, showing the same main absorptions (Figure 4 (a-q)). Montmorillonite is identified in all samples by the broad, asymmetric, and acute absorptions around~1413 nm and~1911 nm, associated with the Al-OH feature at~2206 nm [48,49], shallower than the~1911 nm one, and the slight inflection around~1467 nm [1,49]. It is the sole mineral identified in spectra a-c, e, and o-q ( Figure 4). Others (d and g-n) show additional shallow Al-OH secondary absorptions at~2356 nm and~2445 nm, indicating the mixture with illite [48,49]. The presence of the Al-OH triplet and the slight inflection around~2165 nm allowed the identification of kaolinite in the spectra d and j ( Figure 3) [8,48,49,53]. Figure 4 (l) shows a sharp absorption at~2312 nm and a shallow feature at~2390 nm, indicating that another smectite group mineral may be present, probably saponite, with diagnostic Mg-OH features at~2310 nm and~2380 nm [48,49]. The last three spectra of Figure 4 show distinct mineral assemblages. In spectra r and s, besides the montmorillonite features (broad, asymmetric water absorptions), it was possible to identify: (i) a deep symmetric but broad main Al-OH feature, (ii) two OH − absorptions located at~1574 nm and~1834 nm, and (iii) a deep, broad, asymmetric Mg-OH feature at~2355 nm that allowed the recognition of cookeite (Libearing chlorite) [19,25,54,55]. Spectrum t ( Figure 4) shows illite/white mica features: (i) the sharp and symmetric OH − and Al-OH features at~1413 nm and 2198 nm, respectively and (ii) the well-pronounced Al-OH secondaries at 2345 nm and 2434 nm [47][48][49]. The water feature is deep and asymmetrical but shorter than the Al-OH, diagnosing the existence of illite [49]. Looking at the VNIR portion of the spectra, it was possible to notice a fall off to extreme blue in the reflectance curves (left panel of Figure 4), indicating that iron could be present in several samples [53]. Both ferric and ferrous iron were identified, the latter mainly through its characteristic ramp-like feature around~1150 nm [56]. Notwithstanding, similar features around~1150 nm were interpreted as water features in montmorillonite spectra by Hunt and Salisbury [47]. two OH − absorptions located at ~1574 nm and ~1834 nm, and (iii) a deep, broad, asymmetric Mg-OH feature at ~2355 nm that allowed the recognition of cookeite (Li-bearing chlorite) [19,25,54,55]. Spectrum t ( Figure 4) shows illite/white mica features: (i) the sharp and symmetric OH − and Al-OH features at ~1413 nm and 2198 nm, respectively and (ii) the well-pronounced Al-OH secondaries at 2345 nm and 2434 nm [47][48][49]. The water feature is deep and asymmetrical but shorter than the Al-OH, diagnosing the existence of illite [49]. Looking at the VNIR portion of the spectra, it was possible to notice a fall off to extreme blue in the reflectance curves (left panel of Figure 4), indicating that iron could be present in several samples [53]. Both ferric and ferrous iron were identified, the latter mainly through its characteristic ramp-like feature around ~1150 nm [56]. Notwithstanding, similar features around ~1150 nm were interpreted as water features in montmorillonite spectra by Hunt and Salisbury [47]. Each spectrum has a database code [45] providing information about the type of sample (L-Li-bearing minerals); mineralogy (Pet-petalite); sample provenance (ALB-Alberto mine and BJ-Bajoca mine, FG-Fregeneda, and HD-Hinojosa de Duero); a three-digit numerical part representing the GPS code; a one-or two-digit code representing the spectra number or the averaged spectra numbers; and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4). Each spectrum has a database code [45] providing information about the type of sample (L-Li-bearing minerals); mineralogy (Pet-petalite); sample provenance (ALB-Alberto mine and BJ-Bajoca mine, FG-Fregeneda, and HD-Hinojosa de Duero); a three-digit numerical part representing the GPS code; a one-or two-digit code representing the spectra number or the averaged spectra numbers; and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
Lastly, the position of the~2200 nm Al-OH feature was evaluated for all Li-minerals ( Figure A1, Appendix A) as a possible factor for mineral discrimination.

Characterization of Petalite Alteration Products
The knowledge on the alteration associations can help to refine satellite image processing algorithms and to better understand the spectral behavior of fresh samples. During the field campaigns, it was possible to identify in the Bajoca mine front white to greenish clays preserving petalite cleavage (Figure 5a,b) [57] and yellow to greenish petalite crystals in Hinojosa de Duero. Additionally, in the Bajoca mine front, some petalite relics were found within the heavily altered portions, mainly consisted of friable clay-like minerals. Petalite relics and clayey alteration products were characterized through petrography, XRD, and reflectance spectroscopy. Table 1 summarizes these results and the ICP-OES Li quantification. Appendix A ( Figure A6) provides the XRD diffractograms. Regarding the spectral behavior of the alteration assemblages ( Figure 6), montmorillonite features such as the broad, asymmetric OH − (~1413 nm) and H2O (~1912 nm) absorptions with sharp minima, the Al-OH (~2200 nm) feature, and the small inflection around ~1467 nm [1,48,49] suggests its identification in spectra from a to j. The X-ray diffraction technique identified montmorillonite only in two of the four analyzed samples (Table 1 and Appendix A, Figure A6). In spectra from b to h, illite and/or white mica are also detected by (i) the Al-OH secondaries at ~2350 nm and ~2440 nm, (ii) the slight OH − feature at around ~1790 nm, and (iii) the proportion of Al-OH (~2200 nm) absorption depth to water absorption depth ranging around 1/2 [47][48][49]. The coexistence of illite and white mica was confirmed three samples analyzed through XRD (Table 1 and Appendix A, Figure A6). In Figure 6, the spectra i-k represent more complex mineral assemblages. Kaolinite can be present in spectra i and j because of the Al-OH absorption triplet (~2310 nm-~2350 nm-~2387 nm) and the slight inflection around ~2165 nm [8,48,49,53]. Nonetheless, the deeper 2311-nm absorption in spectra j indicates that a magnesium-bearing mineral such as saponite or sepiolite is also present [48,49]. Oppositely, spectrum k of Figure  6 is dominated by illite and/or white mica features: (i) deep, sharp, and symmetric main Al-OH feature (2196 nm) associated with a sharp, symmetric OH − features (1409 nm); (ii) very pronounced Al-OH secondaries (at 2350 nm and 2433 nm); and (iii) minor OH − feature at 2105 nm [47][48][49]. Most of the spectra of Figure 6 show iron-related absorption features, as observed in petalite fresh samples ( Figure 4). The wavelength position allowed the discrimination of both ferrous and ferric iron [56].  Table 1. Characterization and mineral identification of petalite's alteration products through petrography, XRD, and reflectance spectroscopy. Li quantification through ICP-OES is also shown. Alb-albite; Cok-cookeite (Li-chlorite); Il-illite; Kao-kaolinite; Mnt-montmorillonite; Msc-muscovite; Non-nontronite; Or-orthoclase; Pet-petalite; Plplagioclase; Qtz-quartz; Sap-saponite; Sep-sepiolite; Ser-sericite; WM-white mica. Sample provenance: BJ-Bajoca mine; HD-Hinojosa de Duero. Figure 5c-h shows representative photomicrographs of petalite's alteration minerals. Thin sections of the less altered samples allowed identifying petalite in two out of the three. In the third, petalite is completely replaced despite preserving its cleavage, as previously described by Cardoso-Fernandes et al. [57]. Petalite alteration concentrates along the cleavage (Figure 5c,d), filled by fine white mica and pale brown clays (under parallel polarized light-ppl; Figure 5c), identified as montmorillonite through SEM-EDS [57]. Montmorillonite also forms spherical aggregates more or less aligned (Figure 5c,d). Local sericitization was observed on the crystal borders (Figure 5e,f). Figure 5g,h shows the formation of cookeite (associated with white mica) and clay minerals with pseudospherulitic habit, often interpreted as crystals growing filling pre-existing voids [58].

Sample
Regarding the spectral behavior of the alteration assemblages ( Figure 6), montmorillonite features such as the broad, asymmetric OH − (~1413 nm) and H 2 O (~1912 nm) absorptions with sharp minima, the Al-OH (~2200 nm) feature, and the small inflection around~1467 nm [1,48,49] suggests its identification in spectra from a to j. The X-ray diffraction technique identified montmorillonite only in two of the four analyzed samples (Table 1 and Appendix A, Figure A6). In spectra from b to h, illite and/or white mica are also detected by (i) the Al-OH secondaries at~2350 nm and~2440 nm, (ii) the slight OH − feature at around~1790 nm, and (iii) the proportion of Al-OH (~2200 nm) absorption depth to water absorption depth ranging around 1/2 [47][48][49]. The coexistence of illite and white mica was confirmed three samples analyzed through XRD (Table 1 and Appendix A, Figure A6). In Figure 6, the spectra i-k represent more complex mineral assemblages. Kaolinite can be present in spectra i and j because of the Al-OH absorption triplet (~2310 nm-~2350 nm-~2387 nm) and the slight inflection around 2165 nm [8,48,49,53]. Nonetheless, the deeper 2311-nm absorption in spectra j indicates that a magnesium-bearing mineral such as saponite or sepiolite is also present [48,49]. Oppositely, spectrum k of Figure 6 is dominated by illite and/or white mica features: (i) deep, sharp, and symmetric main Al-OH feature (2196 nm) associated with a sharp, symmetric OH − features (1409 nm); (ii) very pronounced Al-OH secondaries (at 2350 nm and 2433 nm); and (iii) minor OH − feature at 2105 nm [47][48][49]. Most of the spectra of Figure 6 show iron-related absorption features, as observed in petalite fresh samples ( Figure 4). The wavelength position allowed the discrimination of both ferrous and ferric iron [56].
Remote Sens. 2021, 13, x FOR PEER REVIEW 12 of 35 Figure 6. Spectral signature of the petalite's alteration products: reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Each spectrum name is composed by the sample provenance (BJ-Bajoca mine and HD-Hinojosa de Duero), a field code, a one-digit number representing the spectra, and the spectroradiometer employed (asd-FieldSpec 4).
As done for the spectra of Li-minerals, the position of the main Al-OH feature of the petalite alteration products was also assessed in Figure A1 (Appendix A). Additionally, any possible correlation between the Li content and the absorptions features in terms of wavelength, depth, area, and full-width at half maximum (FWHM) was evaluated ( Figure  7). . Spectral signature of the petalite's alteration products: reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Each spectrum name is composed by the sample provenance (BJ-Bajoca mine and HD-Hinojosa de Duero), a field code, a one-digit number representing the spectra, and the spectroradiometer employed (asd-FieldSpec 4).
As done for the spectra of Li-minerals, the position of the main Al-OH feature of the petalite alteration products was also assessed in Figure A1

Spectral Studies for Lithological Classification
This section concerns the comparison of the spectral behavior of the main field lithologies. Additionally, a description of the areas misclassified as Li-pegmatite is presented, as well as the spectral behavior of the samples collected in each of the four areas.
Stacked reflectance and continuum removed curves of the metasedimentary formations from the CXG are presented in Figure 8a-g. White mica diagnostic features such as (i) deepest, sharp Al-OH feature around ~2214 nm, (ii) sharp and symmetric OH − feature (~1413 nm) associated with the absence of a water feature (expected around ~1900 nm), and (iii) weak to pronounced Al-OH secondaries at ~2346 nm and ~2451 nm dominate the spectra of Figure 8a-c [47][48][49]. The spectra of fresh samples (a-c) are also characterized by a ramp-like absorption feature due to the presence of ferrous iron (Fe 2+ ) [56], together with the ~750 nm, ~920 nm, and ~1150 nm Fe 3+ and Fe 2+ features [47,50,54]. These VNIR absorptions, associated with a slight inflection around ~2250 nm, related to the Fe-OH bonds [6,49] or Al(Mg, Fe)-OH bonds [59,60], and a slightly deeper 2350-nm absorption when comparing with the ~2451 nm Al-OH secondary, indicating the contribution of the Mg-OH bonds [48,49], are indicative of the presence of chlorite in fresh samples [47,49,54,59]. Despite the absorption features of spectra d-f of Figure 8 not being evident, there is (i) a weak, broad, rounded water absorption at ~1926 nm that may be associated with fluid inclusions in quartz [49] and (ii) ferric iron-related (Fe 3+ ) absorptions around ~942 nm evidencing sample oxidation ( Figure 8) [49,56]. Finally, spectra g shows the calcite diagnostic features-namely, the carbonate ion (CO3 2− ) feature at 2339 nm and weaker carbonate features around ~1900-2000 nm [49].

Spectral Studies for Lithological Classification
This section concerns the comparison of the spectral behavior of the main field lithologies. Additionally, a description of the areas misclassified as Li-pegmatite is presented, as well as the spectral behavior of the samples collected in each of the four areas.
Stacked reflectance and continuum removed curves of the metasedimentary formations from the CXG are presented in Figure 8 (a-g). White mica diagnostic features such as (i) deepest, sharp Al-OH feature around~2214 nm, (ii) sharp and symmetric OH − feature (~1413 nm) associated with the absence of a water feature (expected around~1900 nm), and (iii) weak to pronounced Al-OH secondaries at~2346 nm and~2451 nm dominate the spectra of Figure 8 (a-c) [47][48][49]. The spectra of fresh samples (a-c) are also characterized by a ramp-like absorption feature due to the presence of ferrous iron (Fe 2+ ) [56], together with the~750 nm,~920 nm, and~1150 nm Fe 3+ and Fe 2+ features [47,50,54]. These VNIR absorptions, associated with a slight inflection around~2250 nm, related to the Fe-OH bonds [6,49] or Al(Mg, Fe)-OH bonds [59,60], and a slightly deeper 2350-nm absorption when comparing with the~2451 nm Al-OH secondary, indicating the contribution of the Mg-OH bonds [48,49], are indicative of the presence of chlorite in fresh samples [47,49,54,59]. Despite the absorption features of spectra d-f of Figure 8 not being evident, there is (i) a weak, broad, rounded water absorption at~1926 nm that may be associated with fluid inclusions in quartz [49] and (ii) ferric iron-related (Fe 3+ ) absorptions around~942 nm evidencing sample oxidation ( Figure 8) [49,56]. Finally, spectra g shows the calcite diagnostic features-namely, the carbonate ion (CO 3 2− ) feature at 2339 nm and weaker carbonate features around~1900-2000 nm [49].
Regarding the granite samples, continuum removed curves (Figure 8 (h-k)), evidence a mixture of montmorillonite/illite and white mica features: (i) asymmetric OH-(~1413 nm) and water (~1913 nm) features, (ii) broad water feature with sharp minima and similar depth to the Al-OH one (~2206 nm), (iii) a small inflection around~1467 nm [1,48,49], (iv) Al-OH secondaries at~2351 nm and~2443 nm, and (v) weaker OH-absorption at 2110 nm [47][48][49]. Fresher samples (spectra h-i) show distinctive VNIR features at~960 nm (Fe 3+ ) and~1150 nm (Fe 2+ ), inferring the presence of biotite. Despite the dominance of montmorillonite and white mica feature, which overlap with the main kandite group absorption features, altered samples (spectra j-k) show slight inflections at~2165 nm and 2382 nm that indicate that a kandite group mineral, like kaolinite or halloysite, is present in the mineral mixture [8,49]. Additionally, goethite was identified in altered samples through its diagnostic ferric iron features at~670 nm and~950 nm and peak in the red region (~750 nm) [49,56]. Each spectrum has a database code [45] providing information about the type of sample (TA-training areas), type of lithology (M-CXG metasediments, G-granite and granodiorite), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
Regarding the granite samples, continuum removed curves (Figure 8h-k), evidence a mixture of montmorillonite/illite and white mica features: (i) asymmetric OH-(~1413 nm) and water (~1913 nm) features, (ii) broad water feature with sharp minima and similar depth to the Al-OH one (~2206 nm), (iii) a small inflection around ~1467 nm [1,48,49], (iv) Al-OH secondaries at ~2351 nm and ~2443 nm, and (v) weaker OH-absorption at ~2110 nm [47][48][49]. Fresher samples (spectra h-i) show distinctive VNIR features at ~960 nm (Fe 3+ ) and ~1150 nm (Fe 2+ ), inferring the presence of biotite. Despite the dominance of montmorillonite and white mica feature, which overlap with the main kandite group absorption features, altered samples (spectra j-k) show slight inflections at ~2165 nm and ~2382 nm that indicate that a kandite group mineral, like kaolinite or halloysite, is present  [45] providing information about the type of sample (TA-training areas), type of lithology (M-CXG metasediments, G-granite and granodiorite), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
The granodiorite samples stacked reflectance and continuum removed curves are displayed in Figure 8 (l-q). White mica features dominate the spectra n to q, as evidenced by (i) the deep, narrow Al-OH absorption at~2207 nm, (ii) the associated Al-OH secondaries (at~2349 nm and~2438 nm), and (iii) minor OH − combined with lattice modes absorption features around~2110 nm [47][48][49]. White mica is also present in Figure 8 (m). The ubiquitous biotite presence in all spectra (l to q) was confirmed by its diagnostic (i) 2253-nm absorption feature, related to the either Fe-OH bonds [6,49] or Al(Mg, Fe)-OH bonds [59,60] and (ii) Mg-OH absorptions at~2330-2350 nm and~2380 nm [48,49]. The higher relative depth of the Mg-OH feature compared to the~2250-nm absorption and the slight inflection/weak absorption around~2430 nm indicate that biotite is Mg-rich (possibly mixed with phlogopite-for example, spectra q) [49,53]. Biotite also occurs mixed with chlorite in spectrum l, being chlorite identified through its diagnostic VNIR features (~750-nm,~915-nm, and~1150-nm ferric and ferrous iron absorptions) [47,50,54,56]. The existence of a doublet OH − absorption in spectra o to q, associated with the slight inflection at~2166 nm and the Al-OH absorptions at~2349 nm and~2383 nm is characteristic of the kandite group [8,48,49,53]. The prominent and deep water feature (~1915 nm) is diagnostic of halloysite [48,49]. Lastly, spectra q of Figure 8 shows a broad, asymmetric water feature that coupled with a unique OH − feature (at~1413 nm), and an Al-OH feature (at 2207 nm) indicates that montmorillonite and/or illite may be present [48,49].
Finally, the Li-bearing pegmatites were separated according to the dominant mineralogy into (i) lepidolite-bearing and (ii) petalite-bearing. Stacked reflectance and continuum removed curves of both lepidolite-and petalite-bearing pegmatites are presented in Figure 9.
As shown in Figure 9 (a-f), white mica related absorptions dominate all spectra from lepidolite-bearing pegmatites (sharp and symmetric main OH − and Al-OH, coupled with the associated secondaries at~2347 nm and~2434 nm and weaker OH − combined with lattice modes absorptions at 1835 nm,~1997 nm, and 2116 nm [47][48][49]). However, features in the VNIR region such as (i) the~551-nm absorption, associated with the incorporation of Fe 3+ and/or Mn 3+ , and (ii) the peak in the red region (left panel of Figure 9), allowed discrimination of the white mica as lepidolite in spectra c-f [50]. In the VNIR region, another iron-related feature was observed at~960 nm (Fe 3+ ) [56]. Some clay minerals, such as montmorillonite and illite, may be present in spectra a, d, and f (Figure 9), identified by the relatively deep and/or asymmetric of the water feature at~1912 nm and the slight inflection around~1467 nm [1,49].
Like the petalite samples, mixtures of montmorillonite and illite dominate most of the spectra of the petalite-bearing pegmatites (Figure 9 (g,h,j-q)). The montmorilloniterelated inflection at~1467 nm is present in all the mentioned spectra [1]. In spectra g, i and j, montmorillonite features prevail over illite features: (i) higher asymmetry of the main OH − feature (~1410 nm), (ii) higher depth of the water feature (~1912 nm) in comparison with the main Al-OH absorption (2192-2208 nm), and (iii) Al-OH secondaries (~2349 nm and~2439 nm) absent or poorly developed [48,49]. In spectra o-q, the opposite happens, with illite being predominant due to: (i) narrower and/or more symmetric main OH − feature, (ii) weaker OH − absorptions 1620 nm and~2108 nm, (iii) higher depth of the main Al-OH feature in relation to the water one, and (iv) well-developed Al-OH secondaries [48,49]. Spectra k-n represent a mixture of montmorillonite and illite. The rounded water absorption in spectra h indicates that free water is present, which could be due to fluid inclusions in minerals such as quartz [49]. White mica, which is discriminated from illite by the absence of a water feature [49], predominates spectra r-t, although the white mica present in spectra r is a hydrated white mica, probably sericite because the depth of the water feature is small than what would be expected for illite [49]. Figure 9 (w) shows diagnostic cookeite features, namely: (i) a broad, symmetric main Al-OH feature, (ii) OH − absorptions at~1581 nm and~1835 nm, and (iii) a deep and broad Mg-OH feature at 2358 nm [19,25,54,55]. Cookeite also occurs mixed with white mica features in spectra t and mixed with montmorillonite/illite in spectra q. Kaolinite was identified in Figure 9 (g,j) through its characteristic Al-OH triplet and slight inflection at 2166 nm [8,48,49,53]. Figure 9 (u,v) shows spectral behavior correspondent to orthoclase, where the absorptions at~1410 nm and~1912 nm may be the result of aqueous fluid inclusions [47,55]. Additionally, the~2110-nm feature is diagnostic of ammonium (NH 4 + ), but its relatively shallow depth can indicate only small amounts of NH 4 + substituting K + and not the NH 4 + feldspar (buddingtonite) [49,61,62]. When looking at the VNIR range, especially in the raw reflectance curves (left panel of Figure 9), goethite was identified in spectra g, h, l, and r, through its diagnostic ferric iron absorptions (Fe 3+ ) at~960 nm and 660 nm and peak in the red region around~760 [49,56]. Other iron-related features were also observed.
Remote Sens. 2021, 13, x FOR PEER REVIEW 17 of 35 Figure 9. Spectral signature of the outcropping lithologies, namely lepidolite-bearing pegmatite and petalite-bearing pegmatite: reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Dotted curves: fresh samples; full curves: samples with visible signs of weathering. Each spectrum has a database code [45] providing information about the type of sample (TA-training areas); type of lithology (LP-Li-pegmatite); sample provenance (AL-Almendra, BJ-Bajoca mine, and HD-Hinojosa de Duero); a three-digit numerical part representing the GPS code; a one-or two-digit code representing the spectra number or the averaged spectra numbers; and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
In the case of the false-positive area no. 2 ( Figure A2f,g), it was possible to recognize hydrated kaolinite or kaosmectite due to the (i) slight inflections ~1407 nm and ~2197 nm indicative of the OH − and Al-OH doublets; (ii) very shallow Al-OH triplet at (~2313 nm, Figure 9. Spectral signature of the outcropping lithologies, namely lepidolite-bearing pegmatite and petalite-bearing pegmatite: reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Dotted curves: fresh samples; full curves: samples with visible signs of weathering. Each spectrum has a database code [45] providing information about the type of sample (TA-training areas); type of lithology (LP-Li-pegmatite); sample provenance (AL-Almendra, BJ-Bajoca mine, and HD-Hinojosa de Duero); a three-digit numerical part representing the GPS code; a one-or two-digit code representing the spectra number or the averaged spectra numbers; and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).
Overall, the misclassified regions correspond to (i) wine or olive plantations, where the agricultural activities exposed metasedimentary boulders, or (ii) hillside and slope areas, where many rock boulders are naturally exposed. Petrography and reflectance spectroscopy studies were conducted to fully characterize the collected samples and understand if there can be any spectral confusion with the target class (Li-pegmatite) and minerals. The results are summarized in Table A1 (Appendix A), while the spectral curves and photomicrographs are presented in Figures A2 and A3, respectively.
As seen in Table A1, most of the samples correspond to CXG metasediments, although some quartzite and pegmatite samples were in some false-positive areas. The microscopic studies of CXG samples allowed characterization as interbed phyllites with metagraywacke layers, with well-preserved sedimentary features such as bedding and load casts ( Figure A3a,b). The mineralogical assemblage identified, composed of quartz, white mica (detrital and metamorphic origin), biotite, chlorite, plagioclase, K-feldspar, detrital zircon, tourmaline, apatite, iron oxides, and other opaque minerals ( Figure A3), is in line with the described in the literature [63]. Quartzite samples have tiny mica crystals in the matrix and exhibit an inherited clastic texture, compatible with the metapsamites described for the São Grabiel's Ordovician formation [37]. Lastly, no signs of Li mineralization were found in the pegmatite samples, which characterizes it as a barren pegmatite. Figure A2 (Appendix A) shows the spectral signature of samples collected in falsepositive sites. The continuum removed curves related to the area no. 1 ( Figure A2 (a-e)) reveal white mica in spectra c (sharp, symmetric main OH − and Al-OH absorptions at 1414 nm and~2210 nm, coupled with the Al-OH secondaries at~2350 nm and~2438 nm and the absence of water absorption). Illite was recognized in spectra b through similar features but with its typical more profound water feature (~1925 nm) [47][48][49]. In Figure A2 (e), montmorillonite deep water absorption is spectrally mixed with Mg-OH-related features at~2350 nm and~2383 nm (slight inflection), possibly due to the presence of chlorite or biotite (the VNIR features are not clear enough to discriminate the magnesium-bearing mineral) [48,49,54]. The wavelength position and rounded aspect of the water feature indicate the presence of free water, possibly related to fluid inclusions in quartz or the presence of iron oxides/hydroxides [49]. The ferric iron features at~680 nm and~960 nm and peak in the red region (~760 nm) allowed identifying the goethite in spectra a, while the remaining VNIR absorptions indicate that other iron oxides/hydroxides are present in spectra b, d, and e of Figure A2 [49,56].
In the case of the false-positive area no. 2 ( Figure A2 (f,g)), it was possible to recognize hydrated kaolinite or kaosmectite due to the (i) slight inflections~1407 nm and~2197 nm indicative of the OH − and Al-OH doublets; (ii) very shallow Al-OH triplet at (~2313 nm, 2349 nm, and~2384 nm); and (iii) the deep water feature at~1921 nm [8,48,49,53]. Spectra g of Figure A2 is dominated by illite features, although there is a slight inflection at a wavelength characteristic of the Mg-OH bond (2379 nm) [47][48][49]. Identical OH − and Al-OH features but with less profound water feature allowed to identify white mica in spectra h [47][48][49], while the rounded water feature has indicated the presence of quartz with fluid inclusions [49]. The remaining spectra (f and g) also show a distinctive wavelength position and roundness of the water feature, indicating free water. Based on the VNIR features, goethite was identified in Figure A2 (f,g), while spectra h also shows some ferric iron absorptions [49,56].
Regarding the spectra of the false-positive area no. 3 ( Figure A2 (i-m)), the similar depth of the water (~1914 nm) and main Al-OH (~2208 nm) features indicate the spectral mixture montmorillonite (~1465-nm inflection and broad, deep, asymmetric water feature with sharp minima) and illite (typical Al-OH secondaries) in spectra i and k [1,48,49]. Montmorillonite also occurs mixed with chlorite/biotite in spectra j (well-pronounced Fe-OH or Al(Mg, Fe)-OH absorption at 2254 nm, associated with the deep Mg-OH absorption at~2348 nm), but the VNIR features are not clear enough to discriminate the two minerals [6,48,49,59,60]. Figure A2 (l) shows a predominance of montmorillonite features (absence of Al-OH secondaries, inflection at 1467 nm, and deep water feature), but the shape of the 1914-nm (water) feature is not typical of the smectite group, with the rounded minima indicating the presence of free water [1,48,49]. Spectra m is dominated by chlorite features (diagnostic VNIR absorptions, inflection at~2000 nm, and Mg-OHrelated absorption at 2348 nm) [47][48][49]54,59] mixed with either white mica or illite [47][48][49]. Looking at the VNIR features in both reflectance and continuum removed curves, goethite can be present in spectra i, while other iron oxides and/or hydroxides (that not goethite or hematite) are present in spectra j-l [49,56].
Finally, in the false-positive area no. 4 ( Figure A2 (n-q)), white mica (sharp, symmetric OH − and Al-OH features, associated with the Al-OH secondaries and the absence of a prominent water feature), is mixed with chlorite (diagnostic VNIR absorptions and 2000-nm inflection) in spectra n and q [47][48][49]54,59]. Spectra o and p show white mica features but a more pronounced, asymmetric water feature at~1914 nm, indicating probable sericite (water absorption deeper than typical muscovite but not deep enough as illite) [49]. Ferric iron-related absorptions occur in spectra o and p of Figure A2.
The reflectance magnitude of the different classes is compared in Figure 10. Additionally, the minimum, maximum, and mean reflectance of each class are also compared in Figure A4 (Appendix A). tion at ~2348 nm), but the VNIR features are not clear enough to discriminate th minerals [6,48,49,59,60]. Figure A2l shows a predominance of montmorillonite fe (absence of Al-OH secondaries, inflection at 1467 nm, and deep water feature), b shape of the 1914-nm (water) feature is not typical of the smectite group, with the ro minima indicating the presence of free water [1,48,49]. Spectra m is dominated by ch features (diagnostic VNIR absorptions, inflection at ~2000 nm, and Mg-OH-relat sorption at 2348 nm) [47][48][49]54,59] mixed with either white mica or illite [47][48][49]. Lo at the VNIR features in both reflectance and continuum removed curves, goethite present in spectra i, while other iron oxides and/or hydroxides (that not goethite or atite) are present in spectra j-l [49,56].
Finally, in the false-positive area no. 4 ( Figure A2n-q), white mica (sharp, sym OH − and Al-OH features, associated with the Al-OH secondaries and the absenc prominent water feature), is mixed with chlorite (diagnostic VNIR absorptions and nm inflection) in spectra n and q [47][48][49]54,59]. Spectra o and p show white mica fe but a more pronounced, asymmetric water feature at ~1914 nm, indicating probable ite (water absorption deeper than typical muscovite but not deep enough as illite Ferric iron-related absorptions occur in spectra o and p of Figure A2. The reflectance magnitude of the different classes is compared in Figure 10. tionally, the minimum, maximum, and mean reflectance of each class are also com in Figure A4 (Appendix A).  Additionally, the wavelength position and depth of absorption features are among the most common criteria for lithological discrimination. Figure 11a shows the distribution of the wavelength of the main Al-OH feature according to the distinct types of samples containing white mica, and Figure 11b relates the position and depth of their Al-OH feature. To better discriminate the CXG metasediments from the granitoids based on the spectral mineralogy, the chlorite and/or biotite composition was evaluated through the Mg-OH wavelength position ( Figure A5, Appendix A). Only the Mg-OH absorption was used, because not all spectra presented the~2250-nm feature well-defined.
Lastly, the spectral resolutions of different satellite or drone-borne sensors (such as Sentinel-2, Landsat-8 Operational Land Imager (OLI), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Prototype Research Instruments and Space Mission technology Advancement (PRISMA), and HySpex Mjolnir) were considered to evaluate the potential use of remote sensing data in Li-pegmatite exploration ( Figure 12). ples containing white mica, and Figure 11b relates the position and depth of their Al-OH feature. To better discriminate the CXG metasediments from the granitoids based on the spectral mineralogy, the chlorite and/or biotite composition was evaluated through the Mg-OH wavelength position ( Figure A5, Appendix A). Only the Mg-OH absorption was used, because not all spectra presented the ~2250-nm feature well-defined.   [49,64]) and (b) relation between depth and position of the Al-OH feature. Lepidolite samples (squares) were plotted for comparison. Circles correspond to samples from training areas, while diamonds represent false-positive areas. Full symbols correspond to samples where white mica was exclusively identified; half-full symbols correspond to mixtures of white mica with minerals that do not overlap the main Al-OH absorption; symbols without fill represent spectra where white mica occurs mixed with minerals whose feature overlaps the Al-OH absorption. Al-OH absorption; symbols without fill represent spectra where white mica occurs mixed with minerals whose feature overlaps the Al-OH absorption.
Lastly, the spectral resolutions of different satellite or drone-borne sensors (such as Sentinel-2, Landsat-8 Operational Land Imager (OLI), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Prototype Research Instruments and Space Mission technology Advancement (PRISMA), and HySpex Mjolnir) were considered to evaluate the potential use of remote sensing data in Li-pegmatite exploration (Figure 12).

Spectral Behavior of Li-Bearing Minerals
The collected lepidolite (Figure 2), spodumene (Figure 3), and petalite ( Figure 4) spectra were compared, and, despite all three minerals showing significant absorption features in the same wavelength regions, the shape, symmetry, and relative depth of the features allowed differentiation of lepidolite from the remaining two Li-minerals. This is because lepidolite has diagnostic absorption features at around 550 nm and 960 nm, as well as a drop-off of the extreme blue (which reflect the purple color of lepidolite) and the white mica typical features due to the hydroxide ion (OH − ) absorption around~1409 nm associated with the principal, sharp Al-OH absorption feature (between 2195 nm and Remote Sens. 2021, 13, 3688 20 of 33 2205 nm) and the Al-OH secondaries (around~2348 nm and~2435 nm, respectively). These VNIR absorption features allow discriminating lepidolite from other white micas [50]. Since hyperspectral and multispectral air-or satellite-borne sensors cover the VNIR wavelength region, there is a high potential for lepidolite identification using remote sensing data. Previous studies where lepidolite was recognized and discriminated from other white micas using airborne HyMap hyperspectral data [65,66] corroborate this statement, although the position of the Al-OH features also played a fundamental part in that discrimination [66] (this topic will be further discussed below). All three Li-minerals show diagnostic ironrelated features in the VNIR region, probably due, in part, to the substitution of aluminum in the minerals' structures (e.g., References [50,67]). In spodumene, iron substitution is observed, as well in the SWIR region (~2240 nm).
Concerning the petalite and spodumene, it was not possible to identify any diagnostic feature because the alteration minerals (montmorillonite with minor illite and kaolinite) dominate their spectral curves. Though, spodumene presents iron-related absorptions at 2240 nm not visible in petalite. Conversely, some petalite samples showed additional cookeite (Li-chlorite) and saponite (Mg-smectite) features not recognized in spodumene samples. However, this is not enough to distinguish between the two Li-aluminosilicates, because cookeite and other clay minerals (such as kaolinite or smectites) can form from either petalite or spodumene alteration (e.g., [68][69][70] and references therein). These results also indicate that it would not be possible to spectrally discriminate spodumene or petalite from montmorillonite/illite, based on the absorption features.
Moreover, several authors have noticed a shift in the wavelength position of the main Al-OH absorption due to different Al contents in white micas [7,49,52,64]. As it can be seen in Figure A1 (Appendix A), the lepidolite samples show the Al-OH feature at shorter wavelengths (≤2005 nm), which indicates a higher Al content [7,49,64]. Although the Al-OH feature position of spodumene and petalite overlaps the lepidolite range, there is a tendency for the Al-OH feature to be located at longer wavelengths ( Figure A1). These wavelength differences could be caused by the spodumene and petalite substitution by montmorillonite and kaolinite instead of a higher Al content in white micas, especially because they are not present in all spectra.
Such subtle changes in wavelength positions are not detectable through multispectral sensors but are noticeable using hyperspectral data. For example, in the previous work conducted by Oshigami et al. [66], the wavelength position of the Al-OH absorptions allowed the distinction of lepidolite from Al-poor and Al-rich muscovite by applying different thresholds to the employed mineral indices. Despite some overlap in the position of the Al-OH feature between lepidolite and the other Li-minerals, a similar approach could be employed to hyperspectral data to map lepidolite-and spodumene/petalitebearing areas.

Effects of Alteration in the Spectral Signature of Li-Bearing Minerals
Comparing the reflectance spectroscopy of altered petalite results ( Figure 6) with the complementary petrographic/geochemical techniques (Table 1), it was possible to corroborate the presence of montmorillonite and illite in the same samples. Overall, the combination of methods was complementary, despite some contradictory mineral identification between XRD and reflectance spectroscopy. The XRD identified nontronite and Fe-bearing smectite, while the reflectance spectroscopy indicated montmorillonite and possibly saponite. Nonetheless, previous studies have shown that reflectance spectroscopy has a greater capacity to distinguish among different clay minerals and between clay minerals and illite and muscovite [48,52].
In general, the petalite alteration products' spectral behavior ( Figure 6) is very similar to the one observed for petalite (Figure 4), with a dominance of montmorillonite and illite features. Kaolinite and saponite (Mg-smectite) were also observed. Regarding the main Al-OH feature ( Figure A1, Appendix A), petalite alteration products show the widest range, reaching higher wavelengths than the ones observed for the three Li-bearing minerals.
There is also a spectral similarity between petalite relics (HD.479 and BJ4-01) and entirely replaced crystals. The petalite relics found within highly altered areas show Li contents compatible with unaltered petalite [13], which may be explained by superficial alteration cover and the concentration of alteration on fractures, cleavage, and crystal border, as observed in Figure 5c-f and Cardoso-Fernandes et al. [57]. Figure 7 shows an inverse correlation between the Li content and the hydroxyl and water features depth (Pearson's r = 0.95 and 0.88, respectively). However, it is not possible to estimate the Li content based on these parameters, because, at low Li concentrations, no trend was observed. From samples with preserved petalite relics to the wholly altered ones, there is an increase in the OH − and H 2 O absorption depth explained by higher amounts of alteration minerals since mineral abundance is related to the absorption intensity [56,64,71]. On its turn, a higher petalite substitution for alteration minerals such as clays leads to a decrease in Li content (Table 1 and Figure 7). The results indicate that the absorption intensity could be used to discriminate fresher petalite samples.
Despite the different climate, the resultant mineral alteration assemblage (montmorillonite and other smectite group minerals, illite, kaolinite) is very similar to the spodumene and petalite's supergene alteration paragenesis found in Brazil under tropical to subtropical climate [25]. Previous studies have shown that these alteration minerals can result from either supergene (weathering) or hydrothermal metasomatic alteration in the late stages of pegmatite crystallization [57,68,72,73], which would explain the occurrence of the same mineral paragenesis under distinct weathering conditions. At the same time, the presence of clay minerals in microfractures in spodumene was already described by Clark et al. [50], causing the spectral behavior of spodumene to be dominated by the clay's absorption features. Therefore, it is not surprising that spodumene and petalite samples show similar spectral curves because the two minerals seem to undergo the same alteration phenomena.
Taking this into account, it is indisputable that the alteration deeply affects the spectral behavior of Li-minerals even when clay minerals are present in minimal amounts. Therefore, any possible spectral signature of petalite or spodumene in the VNIR-SWIR spectral range is obscured by the very pronounced clay-related features. This poses advantages and limitations in satellite-based Li-exploration. Since spodumene and petalite's spectral features are similar and because different alteration degrees produce the same spectral signal, the ability to spectrally identify spodumene or petalite using the same methodology is not compromised. This is in line with the results obtained with classical and machine learning techniques in the study area [21,23]. Additionally, the results indicate no need for different methodological approaches worldwide since the same mineral assemblages could be found in distinct climates. The work of Santos et al. [22] corroborated this interpretation, because the methods developed for Li-pegmatite identification in the Fregeneda-Almendra area also showed potential for Li-exploration in the Minas Gerais region (Brazil). The main limitation is the impossibility of the spectral features allow discrimination between the areas where spodumene and/or petalite occur from the areas where the weathering of other materials produced the clay minerals. This means that as in other deposits, clay minerals can be used as pathfinders for petalite or spodumene. Montmorillonite is identified through the diagnostic, broad, asymmetric, and acute absorptions around~1413 nm and~1911 nm; the Al-OH feature at~2206 nm [48,49] shallower than the~1911 nm one; and the slight inflection around~1467 nm [1,49]. Kaolinite has also diagnostic features, such as (i) a small inflection around~2165 nm; (ii) the presence of an Al-OH triplet at 2312,~2350, and~2380 nm; and (iii) double OH − and Al-OH absorptions (not visible in poorly crystalline kaolinite) [8,48,49,53]. However, this pathfinder approach would only work in areas where argillic alteration is not predominant. In future attempts, the variation of absorption intensity related to proportions of Li content could be used to try to discriminate Li-bearing minerals from barren, clay-bearing areas as observed in Figure 7. The co-occurrence of cookeite (formed after petalite or spodumene alteration) could also be used to highlight Li-bearing areas.

Discrimination of Li-Pegmatites from the Remaining Lithologies
The last objective of this work was to understand if Li-pegmatites can be discriminated from the remaining field lithologies. In general, it was possible to ascertain that all sampled lithologies present similar absorption features caused by water and/or hydroxyl present in the structure of rock-forming minerals (Figures 8 and 9 and Figure A2). This is due to the presence of white mica/illite in samples from all field lithologies and false-positive areas or due to the presence of montmorillonite in some samples.
Lepidolite, cookeite, and orthoclase were exclusively identified in Li-pegmatite samples. These minerals can be used, therefore, as targets or pathfinders for Li-bearing pegmatites. Lepidolite diagnostic features were already described in Section 5.1, while cookeite shows several diagnostic features, namely: (i) a broad, symmetric main Al-OH feature, (ii) two OH − absorptions at~1580 nm and~1835 nm, and (iii) a deep and broad Mg-OH feature at~2355 nm [19,25,54,55]. Orthoclase can be indirectly identified through (i) aqueous fluid inclusions with absorptions at~1410 nm and~1912 nm [47,55] and (ii) NH 4 + substituting K + , with a diagnostic absorption at~2110 nm [49,61,62]. Chlorite takes place in CXG and granodiorite samples. Biotite occurs in both granite and granodiorite samples, but the exclusive occurrence of phlogopite in the granodiorite allows discrimination between the two. Chlorite and/or biotite were also observed in false-positives areas 1, 3, and 4. Montmorillonite and minor kaolinite compound the alteration assemblage of Li-pegmatites, granite, granodiorite, and false-positive areas 1-3. Altered samples of Li-pegmatite, granite, CXG, and false-positive areas 1-3 present iron oxides, predominantly goethite.
As verified, except for certain exclusive minerals or distinct mineral assemblages, other factors may need to be considered to achieve lithological discrimination fully. When considering the reflectance magnitude of the different classes (Figures 10 and A4), in general, altered samples show higher reflectance values when compared to the fresh ones, which is due to alteration minerals (such as clays and iron-oxides). Additionally, petaliteand lepidolite-bearing pegmatites present similar mean reflectance values, although the former can reach higher maximum values. The granite samples demonstrate the highest mean reflectance ( Figure A4). These relative differences in magnitude help to discriminate between lithological types in data from multi-or hyperspectral sensors.
The wavelength position of given absorption features can also be helpful for lithological discrimination (Figures 11 and A5). The metasediments samples of the CXG show the Al-OH absorption position at longer wavelengths which indicates that the white mica is Al-poor ( Figure 11); Li-pegmatite samples present the main Al-OH absorption at shorter wavelengths, similar to lepidolite; granitoid samples show an intermediate wavelength position and, therefore, Al content [7,49,52,64]. Despite the correlation between white mica composition and the wavelength position of the Al-OH absorption, more recent studies have shown that Al-rich phengites can overlap the typical muscovite wavelength range [7,8]. Moreover, lepidolite and Li-rich micas are characterized by a lower Al content than muscovite due to Li substituting Al in octahedral sites [74]. This indicates that another factor is influencing the smaller wavelength of the~2200 nm feature. Additional research, namely detailed mineral chemistry studies, is needed to fully understand this issue. Nonetheless, there is a clear 2205 nm threshold that is useful for remote sensing applications, at least in this study area. Figure 11b indicates that the Al-OH position varies within the different sampled lithologies as the depth of the Al-OH absorption, which is more profound in Li-pegmatite and shallower in CXG samples (either from the training or false-positive areas). Considering this, it is possible to differentiate the CXG metasediments (host rocks) from the target class (Li-pegmatites). The difference between the two classes is even more pronounced when only the most spectrally pure samples are considered (full and half-full symbols of Figure 11). Such wavelength and depth variations could be evaluated using a hyperspectral sensor.
In what concerns the Mg-OH wavelength position ( Figure A5; Appendix A), it is not possible to discriminate between the metasedimentary rocks and granitoids due to the overlap of wavelengths. However, chlorite presents a similar wavelength position and possibly composition in both training and false-positive areas. Except for one Mg-rich chlorite/biotite in the granodiorite samples, all remaining features are compatible with intermediate Mg/Fe compositions [49,59]. In the case of chlorite, the iron-related features in the VNIR range, the shape, position, and relative depth of the~2250 nm and~2350 nm absorptions are compatible with intermediate compositions with both Fe and Mg present. Nonetheless, it should be noted that the overlap with the white mica Al-OH secondaries can lead to shifts in the Mg-OH position.

Potential of Li-Exploration Using Satellite Data
The previous results indicate that the potential use of remote sensing data is highly dependent on the sensors. As seen in Figure 12, the spectral behavior of the CXG metasediments is very distinct from the remaining lithologies, even for the Sentinel-2 and Landsat-8 OLI. However, the Al-OH absorption around~2200 nm is only recognized in higher resolution sensors (ASTER, PRISMA, and Mjolnir). The existence of a band in this region increases the separability of the CXG metasedimentary class. The granodiorite samples are also spectrally distinct from Li-bearing pegmatites. The most significant similarity is observed between granite and Li-pegmatite samples. The hyperspectral sensors can identify the same wavelength shifts observed in the laboratory ( Figure 11) that can help to discriminate them. Nonetheless, PRISMA has advantages over the Mjolnir sensor since the VNIR features can increase the separability of the two classes. On the other hand, when on board aircraft and UAVs, the Mjolnir sensor has a better spatial resolution, reducing the components in spectral mixtures.
Considering the results, it is possible to infer that although the distinct lithologies present similar absorption features, it is possible to spectrally discriminate Li-pegmatites from other lithologies using remote sensing data due to (i) the higher reflectance magnitude, (ii) more profound Al-OH absorption depths, and (iii) shifts in the wavelength position of the Al-OH features. Furthermore, distinct mineral assemblages result in particular spectral behavior (e.g., Li-pegmatites lack chlorite/biotite features in the VNIR region while presenting important and exclusive target minerals). The false-positive areas, for example, despite showing similar clay-related absorptions, also display distinct VNIR features. Of course, the ability to discriminate the target class from the remaining outcropping lithologies increases with an increasing spectral resolution (Figure 12), showing the importance of hyperspectral sensors such as PRISMA or Mjolnir. The capability to spectrally discriminate Li-pegmatites using hyperspectral data needs to be validated in further studies in a real-life application and not only using simulated hyperspectral data (e.g., [8,65,66]).

Conclusions
Analyzing the wavelength, shape, and symmetry of diagnostic absorption features, as well as the overall reflectance magnitude and absorption depths of a Li-dedicated spectral library, allowed the identification of the spectral mineralogy and characterization of samples of (i) different Li-bearing minerals, (ii) several field lithologies of the study area, and (iii) samples from areas previously misclassified as Li-pegmatite using a machine learning algorithm. Ultimately, the results obtained allowed evaluation of the potential of remote sensing and spectral data in Li-exploration, thus fulfilling the proposed objectives and addressing the main research questions:

1.
Can Li-bearing minerals be distinguished from each other and other minerals? The spectra obtained allowed the discrimination of lepidolite from spodumene, petalite, and other white micas. However, it was not possible to distinguish petalite from spodumene because clay-related absorption features dominate the spectra of both Li-minerals, thus obscuring any likely diagnostic spectral behavior. This implies that spodumene and petalite can only be indirectly identified. Several factors, such as the VNIR features or the position of the main Al-OH absorption, enable direct lepidolite identification through multi-or hyperspectral satellite data.

2.
How can the alteration products affect the spectral signature of Li-bearing minerals? The complementary petrography and geochemical studies proved that the same alteration minerals observed in previous studies occur in either well-preserved or completely altered samples. However, the clay-related absorption depths are much smaller for petalite relics. Further studies are needed to understand if the clay minerals formed from petalite alteration in the Fregeneda-Almendra region are related to weathering or metasomatic processes. Therefore, the effect of alteration paragenesis is paramount in the spectral behavior of spodumene and petalite, thus limiting the possibility to discriminate areas where these Li-minerals occur from heavily altered, barren areas through multispectral data. On the other hand, using hyperspectral data could overcome this limitation, either through the possibility of recognizing variable absorption depths (distinct amounts of clay minerals) or by the identification of other alteration minerals exclusively related to Li-minerals such as cookeite. Notwithstanding, in areas where clay-related alteration is not prevalent, the same methodology can be employed, using montmorillonite and kaolinite as pathfinders to identify spodumene and petalite worldwide since the alteration paragenesis is similar to the described in tropical conditions.

3.
Can Li-pegmatites be spectrally discriminated from other lithologies in a supervised classification attempt? Despite the similar absorption features, we conclude that the satellite data can map Li-pegmatites through several criteria such as (i) the distinctive reflectance magnitude, (ii) higher absorptions depths, (iii) distinct wavelength position of key absorption features, and (iv) particular overall spectral behavior reflecting different mineral assemblages. Minerals such as lepidolite, cookeite, and orthoclase can work as targets for Li-pegmatites. Further studies are needed to evaluate the performance of hyperspectral data in real-life conditions and compare it with the laboratory-simulated hyperspectral data.
The results and conclusions withdrawn in this work provide valuable insights on the spectral behavior of pegmatite-hosted Li mineralization and minerals, allowing for more reasoned and directed satellite-based approaches in the future, based on the proposed strategies and pathfinders. Additional studies are needed to continue improving the knowledge on the alteration paragenesis and to understand the smaller wavelength of the~2200 nm feature in lepidolite and its relationship with the Li and Al contents in the octahedral layer.

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 the data; in the writing of the manuscript; or in the decision to publish the results.     Figure A2. Spectral signature of samples collected in areas misclassified as Li-pegmatite in a previous classification attempt: reflectance curves (to the left) and continuum removed curves (to the right). The main absorption features are highlighted, and sample photographs are provided. Dotted curves: fresh samples; full curves: samples with visible signs of weathering. Each spectrum has a database code [45] providing information about the type of sample (SC-spectral confusion), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4). Each spectrum has a database code [45] providing information about the type of sample (SC-spectral confusion), a three-digit numerical part representing the GPS code, a one-or two-digit code representing the spectra number or the averaged spectra numbers, and a reference to the spectroradiometer employed (sr-SR-6500; asd-FieldSpec 4).       Figure A4. Minimum, maximum, and mean reflectance magnitudes of the different training a lithologies. Figure A5. Comparison of the Mg-OH wavelength position of biotites and chlorites from the CXG metasediments and granitoid rocks. Full symbols correspond to biotite; empty symbols correspond to chlorite; half-full symbols represent chlorite and biotite mixtures or samples where it was not possible to discriminate between the two minerals.