Introduction
Micas are one of the most widespread and probably the most compositionally diverse group of the common rock-forming minerals. As such, they are commonly the subject of petrological and mineralogical studies (see reviews in Bailey Reference Bailey1984 and Mottana et al., Reference Mottana, Sassi, Thompson and Guggenheim2002). The crystal structure of mica is favourable for incorporation of a broad spectrum of minor and trace elements including lithium. Contents of Li vary from dozens of ppm in phlogopite from melanocratic rocks (Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a) up to ca. 6.5 wt.% Li2O in polylithionite from fractionated granites and pegmatites (Černý et al., Reference Černý, Staněk, Novák, Baadsgaard, Rieder, Ottolini, Kavalová and Chapman1995). In addition, the Li-rich micas, zinnwaldite and lepidolite, are objects of economic interest as potential sources of Li metal for battery technologies (e.g. the Cínovec/Zinnwald Li deposit, https://www.europeanmet.com). The principles of the classification of Li-bearing micas were defined by Foster (Reference Foster1960), and an updated classification approved by the International Mineralogical Association (IMA) was published by Rieder et al. (Reference Rieder, Cavazzini, Dyakonov, Frank-Kamenetskii, Gottardi, Guggenheim, Koval, Muller, Neiva, Radoslovich, Robert, Sassi, Takeda, Weis and Wones1999).
The correct study of the composition of Li-bearing micas is not trivial. Before the end of the 1980s (Rieder Reference Rieder1970; Monier et al., Reference Monier, Charoy, Cuney, Ohnenstetter and Robert1987; Stone et al., Reference Stone, Exley and George1988, etc.) and rarely later (Du Bray Reference Du Bray1994; Foord et al., Reference Foord, Černý, Jackson, Sherman and Eby1995; Černý et al., Reference Černý, Staněk, Novák, Baadsgaard, Rieder, Ottolini, Kavalová and Chapman1995, Roda et al., Reference Roda, Pesquera and Velasco1995, Brigatti et al., Reference Brigatti, Lugli, Poppi, Ford and Kile2000), in common with other minerals, micas were analysed using classical methods of wet chemistry applied to monomineralic concentrates of the best achievable purity. This method permitted the determination of major and trace element concentrations, including Li, although only as average values of the entire mica sample.
Since the 1980s, electron-probe microanalysis (EPMA) has become a standard method of mineral analysis enabling the study of compositional zoning and heterogeneity of individual mineral grains. Unfortunately, Li cannot be analysed by this method. As a consequence, several researchers have attempted to evaluate the Li contents by a calculated estimation based on EPMA data.
Initially, based on experiments, Monier and Robert (Reference Monier and Robert1986) proposed that Li=F in atomic values. Subsequently, on the basis of analyses of Li-mica concentrates from Cornubian granites, Stone et al. (Reference Stone, Exley and George1988) proposed the equation Li2O = 0.236×SiO2 – 7.56 (wt.%) for trioctahedral micas. Tindle and Webb (Reference Tindle and Webb1990) modified this equation to Li2O = (0.287×SiO2) – 9.552 (in wt.%). Tischendorf et al. (Reference Tischendorf, Gottesmann, Förster and Trumbull1997), on the basis of their evaluation of more than 1200 published mica compositions worldwide, proposed a very similar equation Li2O = (0.289×SiO2) – 9.658 for trioctahedral Mg-poor micas, and equations Li2O = 0.3935×F1.326 and Li2O = 1.579 × Rb2O1.45 (in wt.%) for dioctahedral micas. The equation, valid for Mg-rich (>6 wt.% MgO) trioctahedral micas Li2O = [2.7/(0.35+MgO)] – 0.13 (in wt.%) (Tischendorf et al., Reference Tischendorf, Gottesmann, Förster and Trumbull1997), was subsequently revised to Li2O = [2.1/(0.356+MgO)] – 0.088 (in wt.%) (Tischendorf et al., Reference Tischendorf, Gottesmann and Förster1999). These Tischendorf equations quickly became generally accepted and have been applied to major element data obtained by EPMA in most mica-related papers in the past decades (Roda-Robles et al., Reference Roda-Robles, Pesquera, Gil-Grespo, Torres-Ruiz and De Parseval2006, Wang et al., Reference Wang, Hu, Zhang, Fontan, De Parseval and Jiang2007, Legros et al., Reference Legros, Marignac, Mercadier, Cuney, Richard, Wang, Charles and Lespinasse2016, Zhu et al., Reference Zhu, Wang, Marignac, Cuney, Mercadier, Che and Lespinasse2018, Bouguebrine et al., Reference Bouguebrine, Bouabsa and Marignac2023, etc.).
Neiva (Reference Neiva1987) combined arithmetic means of several electron microprobe spot analyses with chemical analyses of mica concentrates for trace elements (Li analysed by atomic absorption spectroscopy) to obtain average mica compositions. A rather specific approach by Roda et al. (Reference Roda, Keller, Pesquera and Fontan2007) was the application of EPMA + chemical analyses to monomineral concentrates, using Li estimation by another, slightly modified Tischendorf´s equation Li2O = 0.72×F – 0.612 (wt.%) for the construction of structural formulae, although Li was analysed by AAS. In similar settings, Vieira et al. (Reference Vieira, Roda-Robles, Pesquera and Lima2011) applied the equation Li2O = 0.5387×F – 0.1205 (wt.%), and similarly Marignac et al. (Reference Marignac, Cuney, Cathelineau, Lecomte, Carocci and Pinto2020) used another of the Tischendorf-proposed equations: Li2O = 0.697×F + 1.026, for micas from the Panasqueira W deposit.
Several authors have combined major element EPMA compositions with LA–ICP–MS data for trace elements. However, the apparent contents of Li were calculated according to Tischendorf et al. (Reference Tischendorf, Gottesmann, Förster and Trumbull1997) or by similar equations from Si or F values (Martins et al., Reference Martins, Roda-Robles, Lima and De Parseval2012, Xie et al., Reference Xie, Wang, Groat, Zhu, Huang and Cempírek2015, Li et al., Reference Li, Huang, He, Li, Yu and Chen2015, Legros et al., Reference Legros, Marignac, Tabary, Mercadier, Richar, Cuney, Wang, Charles and Lespinasse2018, Launay et al., Reference Launay, Sizaret, Lach, Melleton, Gloaguen and Poujol2021). This approach was probably justified by the difficulty in achieving identical analysis locations for the two analytical methods. For example, Van Lichtervelde et al. (Reference Van Lichtervelde, Grégoire, Linnen, Béziat and Salvi2008) applied EPMA (major elements) and LA–ICP–MS (minor elements including Li) to Tanco micas but used the equation Li2O = 0.782×F + 0.013 (wt.%, F from EPMA) for structural formulae.
Grew et al. (Reference Grew, Bosi, Ros, Kristiansson, Gunter, Halenius, Trumbull and Yates2018) using EMPA interpreted the Li–Al micas from the Sinceni pegmatite as a solid solution of muscovite, polylithionite and trilithionite components, resulting in the equation Li2O = F/9.34×5.92 (wt.%) which was considered as valid for two di- and trioctahedral mica species. Nevertheless, the existence of the muscovite–lepidolite solid solution was questioned recently by Sulcek et al. (Reference Sulcek, Marler and Fechtelkord2023) who interpret the transitional compositions as a mechanical mixture of muscovite with a polylithionite–trilithionite solid solution.
Michaud and Pichavant (Reference Michaud and Pichavant2020) and Monnier et al. (Reference Monnier, Salvi, Melleton, Lach, Pochon, Baily, Beziat and De Parseval2022) have published extensive EPMA and LA–ICP–MS data for micas from Argemela, Portugal and Beauvoir, France, respectively, but provided no correlation of the two analytical methods.
Published analyses of micas with comprehensively analysed and evaluated Li contents are still surprisingly scarce. A combination of EPMA with ion microprobe Li analyses was used by Henderson et al. (Reference Henderson, Martin and Mason1989), Černý et al. (Reference Černý, Staněk, Novák, Baadsgaard, Rieder, Ottolini, Kavalová and Chapman1995) and Charoy et al. (Reference Charoy, Chaussidon and Noronha1995). A combination of EPMA with LA–ICP–MS has gradually become the standard operating procedure for in situ complex Li-mica analyses (Roda-Robles et al., Reference Roda-Robles, Pesquera, Gil-Grespo and Torres-Ruiz2012, Petrík et al., Reference Petrík, Čík, Miglierini, Vaculovič, Dianiška and Ozdín2014, Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a, Reference Breiter, Hložková, Korbelová and Vašinová Galiová2019, Reference Breiter, Ďurišová, Korbelová, Lima, Vašinová Galiová, Hložková and Dosbaba2022, Reference Breiter, Ďurišová, Korbelová, Vašinová Galiová and Hložková2023a, Reference Breiter, Vašinová Galiová, Hložková, Korbelová, Kynický and Costi2023b).
Most of the aforementioned methods of Li estimation are based on a statistical treatment of datasets, i.e. on the correlation between the contents of Li and other elements, mostly SiO2, F, and Rb2O, and are only valid for the tri- or dioctahedral micas. Although the general error of estimation can be minimised by processing large data sets, the error of the estimates for individual samples/spots remains large (Thiergärtner, Reference Thiergärtner2010). Methods based on structural considerations (Monier and Robert, Reference Monier and Robert1986, Grew et al., Reference Grew, Bosi, Ros, Kristiansson, Gunter, Halenius, Trumbull and Yates2018) have been treated as independent of octahedral site occupancy; their use in petrological practice are evaluated below in the Discussion.
This paper is based on an extensive set of complex local analyses obtained in our labs in the past years (Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a, Reference Breiter, Hložková, Korbelová and Vašinová Galiová2019, Reference Breiter, Ďurišová, Korbelová, Lima, Vašinová Galiová, Hložková and Dosbaba2022, Reference Breiter, Ďurišová, Korbelová, Vašinová Galiová and Hložková2023a, Reference Breiter, Vašinová Galiová, Hložková, Korbelová, Kynický and Costi2023b). The objectives are to show (1) the bias in Li content estimations according to the Tischendorf´s proposals; and (2) the limits of using estimated values in petrological practice.
Samples
The datasets evaluated comprise micas from: representative granitoids of the Bohemian Massif, Czech Republic, ranging from phlogopites from ultramafic dykes to zinnwaldites from stanniferous rare-metal granites (RMG) of the Nejdek and Cínovec plutons (Breiter Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a, Reference Breiter, Hložková, Korbelová and Vašinová Galiová2019); muscovites to lepidolites from the Argemela rare-metal granite, Portugal (Breiter et al., Reference Breiter, Ďurišová, Korbelová, Lima, Vašinová Galiová, Hložková and Dosbaba2022); muscovites from the Panasqueira tungsten deposit, Portugal (Breiter et al., Reference Breiter, Ďurišová, Korbelová, Vašinová Galiová and Hložková2023a); zinnwaldites to lepidolites from the Beauvoir rare-metal granite, France; biotite, muscovite, phengite, zinnwaldite and lepidolite from the Orlovka Ta deposit, Siberia; biotite to zinnwaldite from the Wiborg batholith, Finland; and biotite to lepidolite from the Madeira pluton, Brazil (all in Breiter et al., Reference Breiter, Vašinová Galiová, Hložková, Korbelová, Kynický and Costi2023b). Though the Nejdek, Argemela, Panasqueira and Beauvoir plutons represent strongly peraluminous P-rich granites, the Cínovec, Wiborg and Orlovka plutons represent only slightly peraluminous P-poor post-orogenic granites. The Madeira pluton is an example of a transition from meta-aluminous to peralkaline anorogenic granites.
For additional geological information on the plutons investigated see Badanina et al. (Reference Badanina, Veksler, Thomas, Syritso and Trumbull2004) for Orlovka, Costi et al. (Reference Costi, Dall´Agnol, Pichavant and Rämö2009) and Bastos Neto et al. (Reference Bastos Neto, Pereira, Ronchi, Lima and Frantz2009) for Madeira, Raimbault et al. (Reference Raimbault, Cuney, Azencott, Duthou and Joron1995) and Monnier et al. (Reference Monnier, Salvi, Melleton, Lach, Pochon, Baily, Beziat and De Parseval2022) for Beauvoir, Marignac et al. (Reference Marignac, Cuney, Cathelineau, Lecomte, Carocci and Pinto2020) and Launay et al. (Reference Launay, Sizaret, Lach, Melleton, Gloaguen and Poujol2021) for Panasqueira, Michaud and Pichavant (Reference Michaud and Pichavant2020) for Argemela, Lukkari et al. (Reference Lukkari, Thomas and Haapala2009) for Wiborg, and Breiter et al. (Reference Breiter, Ďurišová, Hrstka, Korbelová, Hložková Vaňková, Vašinová Galiová, Kanický, Rambousek, Knésl, Dobeš and Dosbaba2017b) for Cínovec.
Data from our laboratory are supplemented with EPMA + SIMS data from Cornwall, UK (Henderson et al., Reference Henderson, Martin and Mason1989) in order to: (1) show data obtained using methods other than LA–ICP–MS; and (2) to discuss Cornwall as one of the classical provinces of rare-metal granites. In addition, a limited published dataset from the Tanco pegmatite (Van Lichtervelde et al., Reference Van Lichtervelde, Grégoire, Linnen, Béziat and Salvi2008) allows comparison with analysed (LA–ICP–MS) and estimated (by the authors on the basis of F contents) Li contents in micas from a rather different environment of extremely fractionated Li–Cs–Ta-enriched (LCT) pegmatite.
In the following text we used the IMA terminology (Rieder et al., Reference Rieder, Cavazzini, Dyakonov, Frank-Kamenetskii, Gottardi, Guggenheim, Koval, Muller, Neiva, Radoslovich, Robert, Sassi, Takeda, Weis and Wones1999) supplemented with group-names biotite (solid solution of annite, phlogopite and siderophyllite), lepidolite (solid solution of polylithionite and trilithionite), phengite (dioctahedral mica close to muscovite–celadonite join) and zinnwaldite (trioctahedral mica close to the siderophyllite–polylithionite join).
Methods
The contents of major elements in micas were analysed using a CAMECA SX100 electron microprobe housed at the Institute of Geology of the Czech Academy of Sciences, Praha. An accelerating voltage of 15 kV, a beam current of 10 nA and a beam diameter of 2 μm were applied. The following standards were used: Na, Al – jadeite; Mg, Si, Ca – diopside; K – leucite; Ti – rutile; Mn – MnCr2O4; Fe – magnetite; F – fluorite; Rb – RbCl; Cs – pollucite; and Zn – willemite. Counting times on each peak were optimised for individual elements according to their expected concentrations (10–60 s), and half that time was used to obtain background counts. X-ray lines and background offsets were selected to minimise interference. The X-Phi correction procedure (Merlet, Reference Merlet1994) was applied. An in-house standard of Li–Fe-mica was analysed in every analytical session to monitor stability of results, namely F and Rb. Analytical data of micas were recalculated to the proposed structural formulae based on 44 negative charges. The average detection limits (3σ) under the operating conditions were as follows: 0.01 wt.% for Mg; 0.02 wt.% for Ca and Na; 0.05 wt.% for Al, Mn, Rb, Cs and K; 0.06 wt.% for Si and Ti; 0.07 wt.% for Zn; 0.08 wt.% for F; and 0.19 wt.% for Fe.
The contents of trace elements in mica samples were analysed using two different LA–ICP–MS instrumentation configurations. The solid-state Nd:YAG laser (UP 213) working at a wavelength of 213 nm (New Wave Research, Inc., Fremont, California, USA) coupled to quadrupole-based ICP mass spectrometer Agilent 7500ce and installed at the Department of Chemistry, Masaryk University Brno, with an average detection limit of 1.8 ppm Li was used to analyse all the samples from the Bohemian Massif (Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017, Reference Breiter, Hložková, Korbelová and Vašinová Galiová2019). The second LA–ICP–MS instrumentation, housed at Faculty of Chemistry, Brno University of Technology and BIC Brno, consists of the ArF* excimer laser ablation system Analyte Excite+ (Teledyne CETAC Technologies, Omaha, Nebraska, USA), which emitted the laser beam at a wavelength of 193 nm, connected to quadrupole ICP mass spectrometer Agilent 7900. The average detection limit was 1.7 ppm Li. This instrumentation was used to analyse all other samples (Breiter et al., Reference Breiter, Ďurišová, Korbelová, Lima, Vašinová Galiová, Hložková and Dosbaba2022, Reference Breiter, Ďurišová, Korbelová, Vašinová Galiová and Hložková2023a, Reference Breiter, Vašinová Galiová, Hložková, Korbelová, Kynický and Costi2023b). In all cases, Li contents were quantified using standards SRM NIST 610 and 612, and Si and Al as internal reference elements. An in-house mica standard was included in each analytical session to monitor the stability of results. For details on all analytical settings see the referred papers.
This contribution is based on 3000 individual spot analyses for which lithium contents were obtained using LA–ICP–MS and combined with ca. 2000 individual spots where EPMA contents of major elements were obtained. The purpose of the older studies from the Bohemian Massif was to evaluate Li contents in micas to assess the potential of micas as Li resources (Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a, Reference Breiter, Ďurišová, Hrstka, Korbelová, Hložková Vaňková, Vašinová Galiová, Kanický, Rambousek, Knésl, Dobeš and Dosbaba2017b). In such cases, 15–20 EPMA spots were combined with 20–25 laser-ablation spots spread over one thin section and resulting means were presented as one sample. Later, with the intention to express better local compositional variability, typically five mica grains within each thin section were evaluated (Breiter et al., Reference Breiter, Ďurišová, Korbelová, Vašinová Galiová and Hložková2023a, Reference Breiter, Vašinová Galiová, Hložková, Korbelová, Kynický and Costi2023b); i.e. 3–5 laser ablation spots were combined with 2–3 EPMA spots in each mica grain according to their size. The means of both methods were coupled and are presented here as one analysis, i.e. usually 5 analyses per one sample. In petrological studies from Argemela (Breiter et al., Reference Breiter, Ďurišová, Korbelová, Lima, Vašinová Galiová, Hložková and Dosbaba2022) and the Bohemian Massif (unpublished), one EPMA and one laser ablation spot, successively in exactly the same place were realised, and coupled data from each spot are presented here, in total giving 560 analyses plotted in the figures.
Compositions of the micas were recalculated to the proposed structural formulae on the basis of 44 negative charges. The water content (OH−) and trivalent iron content were not calculated due to the large uncertainties of such calculations for Li-rich micas. All results are available in the Supplementary Tables with an indication of how many LA-analyses of Li were coupled to obtain the value presented. In addition, the totals of octahedral occupation, crucial for choosing the correct equation, are shown. The same data are presented in Figs 1–3.
Results
In the text below, estimations of Li contents according to Tischendorf et al. (Reference Tischendorf, Gottesmann, Förster and Trumbull1997, Reference Tischendorf, Gottesmann and Förster1999) are compared with the real Li contents analysed by LA–ICP–MS. Equation Li2O = (0.289×SiO2) – 9.658 was used for trioctahedral Mg-poor micas, equation Li2O = [2.1/(0.356+MgO)] – 0.088 for Mg-rich trioctahedral micas, and equations Li2O = 0.3935×F1.326 and Li2O = 1.579×Rb2O1.45 were applied to dioctahedral micas.
Data for trioctahedral micas are illustrated in Fig. 1. The Li-enriched Mg-poor micas from the Bohemian Massif, (Nejdek and Cínovec plutons from the Erzgebirge) show a good correlation between the measured and estimated values, namely in the interval between 1–3 wt.% Li2O (Fig. 1a). In the case of Li-poor biotite from the Central Bohemian Pluton and South Bohemian Pluton (Li2O<0.5 wt.%, Fig. 1b) and Li-rich zinnwaldite (Li2O>3 wt.%, Fig. 3a), the dispersion is rather larger. It is worth mentioning that no difference was found in reproducibility of estimated values between magmatic and hydrothermal (greisen) micas.
Data for Mg-rich trioctahedral micas, phlogopite to Mg-rich annite (Fe/(Fe+Mg) = 0.19–0.57) from the Bohemian Massif are shown on Fig. 1b. All samples represent ultramafic dykes, durbachitic rocks and geochemically less evolved granitoids of the Central Bohemian Pluton (details in Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a). As is evident from the figure, the estimations give a generally valid assumption for low Li contents, but the dispersion of individual analysed vs. estimated values is rather wide, giving both over- and underestimated values.
The Madeira pluton in central Brazil comprises three mica-bearing rock types. Although the estimations of Li in Li-poor biotite (<0.5 wt.% Li2O) from the amphibole–biotite granite and slightly Li-enriched biotite from the biotite granite give acceptable results, Li contents in all samples of Li-enriched annite (0.5–1.5 wt.% Li2O) from the peralkaline granites are strongly overestimated (Fig. 1c). The explanation is simple: micas from Al-deficient peralkaline granites contain unusually low IVAl, i.e. they are relatively Si-enriched compared to micas from peraluminous rocks. Higher than ideal Si occupation in the tetrahedra is compensated by vacancies in the octahedral layer. As a consequence, calculations based on SiO2 contents give unrealistically high Li values. Lepidolite from the cryolite granite (5.6–6.2 wt.% Li2O) gives a good correlation.
The Orlovka layered granite pluton comprises several facies bearing dioctahedral micas and three facies with trioctahedral mica species: biotite; zinnwaldite; and lepidolite. Whereas Li estimations for zinnwaldite and lepidolite give very acceptable results, all biotite values are slightly overestimated due to high IVSi occupancy (Fig. 1d).
The rather monotonous Wiborg batholith contains Li-poor biotite while the adjacent Kimi stock contains Li-rich biotite in the matrix and large flakes of macroscopically black mica in the border pegmatite facies. The Li estimations are generally acceptable for Li-poor micas, although a closer look at Wiborg biotite reveals a relative large error for grains, probably affected by mild alteration (Fig. 3b). The values for pegmatitic Li–Fe-mica are strongly overestimated (Fig. 1e): the unusually high IVSi-occupation here decreases the hexahedral occupation to ca. 4.9, i.e. this mica, regardless of its biotite appearance, is already dioctahedral and the use of this calculation method is not entirely correct. (Note that all equations proposed for dioctahedral micas in this case give similarly overrated values.)
Micas in the vertically zoned Beauvoir granite stock evolved upwards from zinnwaldite to lepidolite. Analysed and calculated values correlate well, albeit the majority of the estimated values are underestimated by ca. 0.5–1 wt.% Li2O (Fig. 1f). This reflects the strongly peraluminous character of Beauvoir rocks as is evident in the high IVAl and low Si contents (siderophyllite component) of the micas. The equation best fitting this set of Li-rich trioctahedral micas from Beauvoir is Li2O = 0.257 × SiO2 – 7.388 (wt.%).
Leucogranites from Argemela contain mica species of the muscovite–lepidolite series, i.e. micas crossing the di/tri-octahedral occupation border; the trioctahedral mica (lepidolite) data are presented in Fig. 1g. Because the exact Li content is not known one cannot reliably calculate the layer occupation and choose the appropriate equation for Li estimation, in addition compositionally named dioctahedral micas with occupation approaching the di/tri-octahedral border are shown here. Hence the general overestimation of lithium is noticeable, still increasing with decreasing value of octahedral occupation. The reason is the low Fe-content combined with the high content of Si in this specific mica.
Mica varieties from Li-rich biotite to lepidolite from Cornwall (Henderson et al., Reference Henderson, Martin and Mason1989) are shown in Fig. 1h. With the exception of one relatively Li-poor sample, the estimated values are generally valid. However, on closer inspection nearly identical estimated values of 4.0–4.5 wt.% were obtained for samples with real contents between 3.5 and 5.2 wt.% Li2O. This implies relatively significant differences between the real and estimated values in many situations.
Data for dioctahedral micas of the muscovite–Li-muscovite and muscovite–Li-phengite series are presented in Fig. 2, with values estimated using F and Rb2O according to Tischendorf et al. (Reference Tischendorf, Gottesmann, Förster and Trumbull1997) shown for comparison (Fig. 2c,d). Samples from all presented RMG plutons show large differences between Li values estimated from the F and Rb contents. Values calculated from F contents are generally better correlated with the measured values, but strongly overestimated, whereas calculations based on Rb contents tend to be underestimated. Note however that the muscovites are generally Li-poor (mostly below 0.8 wt.% Li2O in Argemela, and 0.5 wt.% Li2O in Panasqueira), so relative errors of estimation for such low Li are inevitably higher.
A notable difference between magmatic and late hydrothermal muscovite was found at Argemela. The F-based Li contents are overestimated but show a good positive correlation for the magmatic muscovite cores, whereas an apparent negative correlation was found in late hydrothermal muscovite overgrowth samples. Values calculated from Rb show a rather poorer correlation with analytical data (Fig. 2b). Both of the F- and Rb-based estimations depend on the octahedral occupation and spots with higher occupation approaching the di/trioctahedral border tend to show extreme dispersion (F-based estimation, Fig. 2a), or to be underestimated (Rb-based estimation, Fig. 2b). Unacceptable dispersion in the case of F-based estimation of micas with octahedral occupation between 4.76–5.00 (Fig. 2a) is caused by a large variation in the F content of this population.
Muscovite from the Panasqueira pluton is poor in Li. The F-based estimation gives rather high values (Fig. 2c), whereas the Rb-based estimation is relatively better in this specific case (Fig. 2d). Muscovite and phengite from Orlovka show overestimated, but well correlated F-based Li values (Fig. 2e) and strongly underestimated values with large dispersion in the case of Rb-based estimation (Fig. 2f). Strong underestimation of Rb-based values at Orlovka and Argemela can be explained by preferential incorporation of Rb in coexisting K-feldspar.
Discussion
The first attempt to estimate the Li contents of mica from microprobe data was made on the basis of results from an experimental study of the muscovite–biotite miscibility gap by Monier and Robert (Reference Monier and Robert1986) who suggested Li=F in atomic proportions. This proposal found only limited response (i.e. Černý et al., Reference Černý, Chapman, Teerstra and Novák2003), although recent experiments with micas of the Li-muscovite–lepidolite series (Sulcek et al., Reference Sulcek, Marler and Fechtelkord2023) confirmed the strong Li–F affinity in this mica species. The application of the Monier and Robert (Reference Monier and Robert1986) equation to Argemela samples shows a good positive correlation but a strong overestimation for low Li contents combined with a wide data dispersion in the case of higher Li contents (Fig. 3c). We can only speculate that the promising experimental results were based on stable conditions during the experiments, whereas natural micas have crystallised under highly variable Si–Al–Li–F(+Fe) ratios.
The majority of published, and used in petrological practice, attempts of Li estimation (Stone et al., Reference Stone, Exley and George1988, Tindle and Webb Reference Tindle and Webb1990, Tischendorf et al., Reference Tischendorf, Gottesmann, Förster and Trumbull1997, Reference Tischendorf, Gottesmann and Förster1999) are based on statistical processing of compositional data of various sizes. If the number of evaluated analyses exceeds 1000, as in the paper of Tischendorf et al. (Reference Tischendorf, Gottesmann, Förster and Trumbull1997), the calculated error of the proposed equation might appear to be acceptable. Nonetheless, even in these cases the difference between the estimated Li and measured contents in individual samples remains unacceptably wide, as shown in Fig. 3a,b. This is because lithium is incorporated in the crystal lattice of trioctahedral micas in two possible ways: (1) Fe+Fe ⇔ Al+Li and (2) Al+Fe ⇔ Si+Li. The first substitution operates only in the octahedral layer, whereas the second combines changes in the composition of the octahedral layer and in the tetrahedra. In natural samples, the two types of reactions combine in different proportions depending on the actual activity of Si, Al, and Fe (Breiter et al., Reference Breiter, Vaňková, Vašinová Galiová, Korbelová and Kanický2017a). In dioctahedral micas of the muscovite–Li-muscovite–lepidolite series, theoretical exchange vectors: (1) VIAl+2□ ⇔ 3Li; and (2) IVAl+VIAl+□ ⇔ Si+2Li are combined, although actual miscibility between Li-muscovite and lepidolite has been recently questioned by Sulcek et al. (Reference Sulcek, Marler and Fechtelkord2023). Moreover, the sum of cations in the octahedral layer is not constant, changing continuously from 6 to 4 without any discontinuities between tri- and dioctahedral micas. Furthermore, part of the iron present in the samples is typically oxidised to Fe3+ (up to 10% according to unpublished author's data from the Erzgebirge). This suggests that a correct calculation of Li contents from the contents of the major elements (EPMA) is not practical.
Another important fact to be stressed is that the Li–Si relationships in trioctahedral micas and the Li–F relationship in the muscovite–lepidolite series are at least partly structurally conditioned, albeit being disturbed by changes in local PTX conditions. No structural relationship exists between Li and Rb in dioctahedral micas. Consequently, equations Li2O = 1.579×Rb2O1.45 proposed for dioctahedral micas only express the general trend of enrichment of relatively incompatible elements during magma fractionation. In this process, the relative behaviour of Li versus Rb is controlled by several factors, such as the regional specialisation and timing of mica versus K-feldspar crystallisation and water saturation/separation (Fig. 3d). Thus, the Rb-based estimation of Li in dioctahedral micas is even less correct and a riskier procedure than using other proposed estimations.
All proposed equations are sensitive to octahedral site occupation, i.e. are valid for di- or trioctahedral species only, which complicates their correct application to mica of the muscovite–lepidolite series. Especially in LCT pegmatites, micas evolve from muscovite to lepidolite, crossing the formal di- and trioctahedral border. In such cases, some authors applied an estimation based on the F content for the whole series. Equation Li2O = 0.782×F + 0.013 was inferred for the Tanco pegmatite by Van Lichtervelde et al. (Reference Van Lichtervelde, Grégoire, Linnen, Béziat and Salvi2008) based on a limited number of samples analysed using LA–ICP–MS (Fig. 3e). This “locally adjusted” equation was then used for Li addition to EPMA data of other samples. Such a method of estimation seems to be the best possible option for Li estimation when actual local Li analysis is not possible. Grew et al. (Reference Grew, Bosi, Ros, Kristiansson, Gunter, Halenius, Trumbull and Yates2018) in study of Li–Al micas from LCT Sinceni pegmatite, Swaziland, concluded that these micas represent a fine-grained admixture of Li,F-free muscovite and fully Li,F-saturated masutomilite, polylithionite and trilithionite components; this resulted in an empirical equation Li2O = F/9.34×5.92 wt.%, which can be simplified to Li2O = 0.644×F (wt.%). Notice the relatively small difference to the equation by Van Lichtervelde et al. (Reference Van Lichtervelde, Grégoire, Linnen, Béziat and Salvi2008). We applied this equation to Argemela samples (Fig. 3f) and also found in this case, a marked dependence on the octahedral occupation: a decrease in occupation (<4.75) causes a relative rise in estimated Li values. The reason is probably that the Sinceni micas originate at stable conditions, whereas Argemela samples represent several successive stages of magmatic/hydrothermal evolution.
Monier and Robert (Reference Monier and Robert1986) conclusions were that the gap between di- and trioctahedral Li-free micas is large, and its width decreases with increasing Li content and that for a sufficient Li content (Li > 0.6 apfu based on 11 oxygens) a single Li-mica crystallised. In contrast, Sulcek et al. (Reference Sulcek, Marler and Fechtelkord2023) interpreted the Li–Al micas as mixture of dioctahedral Li,F-free muscovite with trioctahedral Li,F-saturated solid solution of trilithionite and polylithionite, i.e. an immiscibility gap divides all di- and trioctahedral micas and not only Li-poor species. Although the gap between muscovite and biotite was intuitively expected (muscovite–biotite granites are common, mixed mica compositions are not known), the gap between muscovite and lepidolite is seemingly at odds with thousands of electron microprobe analyses, covering the entire range between muscovite and trilithionite. This observation undoubtedly calls for further research. We can only note that the Li/F ratios in ideal muscovite, trilithionite and polylithionite are 0/0, 1.5/2 and 2/2, respectively (Rieder et al., Reference Rieder, Cavazzini, Dyakonov, Frank-Kamenetskii, Gottardi, Guggenheim, Koval, Muller, Neiva, Radoslovich, Robert, Sassi, Takeda, Weis and Wones1999). If Sulcek et al., are correct, the Li/F ratio in mixed micas should be within the interval from 1.5/2 to 2/2, i.e. 0.75–1. In actual Li–Al micas we found the Li/F ratio to be ca. 1 at Beauvoir, 0.5 at Argemela and only 0.2 at Panasqueira. Thus, the Li/F ratio in magmatic/hydrothermal micas is influenced not only structurally, but also, and most importantly, by the composition of the parental magma/fluid.
Conclusive recommendation
As demonstrated above none of the lithium estimation methods can fully compensate for the actual detemination of Li by LA–ICP–MS or SIMS at the same analytical site as used to obtain EMPA data. Generally, SiO2-based estimations of Li contents for trioctahedral micas provide a better match to the analysed values than F-based estimations for dioctahedral micas. The Rb-based estimations for dioctahedral micas do not provide applicable results. The use of Si- and F-based estimations can be accepted for general petrological studies where arithmetic means of a significant number of analysed spots, at least in part eliminating the errors of individual estimations, give realistic values for particular mica species. In the case of zoned micas or a superposition of two or more mica populations, it has to be considered that the substitutional mechanisms and hence also the Li–Si(F) relations might differ markedly in the individual mica zones or populations. This has hampered the use of the estimation method for a detailed interpretation of mica chemistry.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/mgm.2023.72.
Acknowledgements
This study was supported by the Czech Science Foundation Project No. EXPRO 19-29124X at the BIC Brno and by RVO 67985831 at the Institute of Geology of the Czech Academy of Sciences. Louis Raimbault, Hilton Tulio Costi and Elena Badanina are thanked for providing samples from Beauvoir, Madeira and Orlovka, respectively. Associate Editor E.S. Grew, Igor Petrík and one anonymous reviewer are thanked for a detail constructive review which helped us to improve the original version of the manuscript significantly.
Competing interests
The authors declare none.