Cyrielle Bernard

Guillaume Estrade
email: guillaume.estrade@get.omp.eu

Stefano Salvi

Didier Béziat

Martin J Smith

Guillaume Estrade A⸸

Alkali pyroxenes and amphiboles: a window on rare earth elements and other high field strength elements behavior through the magmatic-hydrothermal transition of peralkaline granitic systems

Keywords: Peralkaline granite, pegmatite, pyroxene, amphibole, rare earth elements, high field strength elements

Introduction

High Field Strength Elements (HFSE, namely Nb, Ta, Zr, Hf, U, Th, Sn, and Ti), and in particular Rare Earth Elements (REE), are among raw materials considered the most critical today (European Commission 2018). They are used in many modern technologies linked to the transition to renewable energy infrastructure [START_REF] Lucas | Rare earths: science, technology, production and use[END_REF][START_REF] Goodenough | The Rare Earth Elements: Demand, Global Resources, and Challenges for Resourcing Future Generations[END_REF]), but their production is quite restricted, worldwide. REE deposits can be primary (igneous, carbonatites and alkaline systems) or secondary (placers and ion-adsorption) [START_REF] Chakhmouradian | Rare earth elements: minerals, mines, magnets (and more)[END_REF] and economic concentrations of heavy REE (HREE, Gd to Lu) are rather rare compared to those of light REE (LREE, La to Eu). The latter occur in high amounts mostly in carbonatites and [START_REF] Larsen | The Ilímaussaq intrusion-progressive crystallization and formation of layering in an agpaitic magma[END_REF][START_REF] Chakhmouradian | Rare earth elements: minerals, mines, magnets (and more)[END_REF]. Most alkaline igneous rocks are silicaundersaturated, with the silica-saturated variety (i.e., granitic) being quite rare (e.g. [START_REF] Larsen | The Ilímaussaq intrusion-progressive crystallization and formation of layering in an agpaitic magma[END_REF][START_REF] Ranløv | Compositional zoning in hydrothermal aegirine from fenites in the Proterozoic Gardar Province, South Greenland[END_REF][START_REF] Foland | Formation of cogenetic quartz and nepheline syenites[END_REF][START_REF] Kramm | Nd and Sr isotope signatures of the Khibina and Lovozero agpaitic centres, Kola Alkaline province, Russia[END_REF][START_REF] Sheard | Controls on the Concentration of Zirconium, Niobium, and the Rare Earth Elements in the Thor Lake Rare Metal Deposit, Northwest Territories, Canada[END_REF]). However, it is the latter that have high Yb/La ratios compared to other varieties (Fig 1). Pegmatites, which form from the last and most evolved part of silica saturated melts, concentrate the highest amounts of REE and other HFSE.

There is a general consensus in the literature that magmatic processes play a role for a pre-enrichment of alkaline rocks in REE and other HFSE [START_REF] Marks | A global review on agpaitic rocks[END_REF]. These include sourcing from undepleted, more or less metasomatized mantle, crustal contamination, fractional crystallization, and melt-melt immiscibility during magma cooling [START_REF] Bonin | A-type granites and related rocks: evolution of a concept, problems and prospects[END_REF][START_REF] Nardi | A-type granitic rocks in post-collisional settings in southernmost Brazil: their classification and relationship with tectonics and magmatic series[END_REF][START_REF] Chakhmouradian | Rare earth mineralization in igneous rocks: sources and processes[END_REF][START_REF] Veksler | Partitioning of elements between silicate melt and immiscible fluoride, chloride, carbonate, phosphate and sulfate melts, with implications to the origin of natrocarbonatite[END_REF]. There is also evidence that hydrothermal processes may play a key role in concentrating these elements to ore grades (e.g. [START_REF] Salvi | The role of hydrothermal processes in the granite-hosted Zr, Y, REE deposit at Strange Lake, Quebec/Labrador: evidence from fluid inclusions[END_REF][START_REF] Migdisov | Hydrothermal mobilisation of the Rare Earth Elements-a Tale of" Ceria" and[END_REF]. However, to date, the relative importance of magmatic versus hydrothermal processes is still debated, partly because the processes occurring at the magmatic-hydrothermal transition are hard to document and most studies were carried out in the laboratory (e.g. [START_REF] Migdisov | Hydrothermal transport, deposition, and fractionation of the REE: Experimental data and thermodynamic calculations[END_REF]).

To fill this gap, it is important to understand the timing of concentration of the REE in natural occurrences.

To do so, in this paper we propose for the first time a detailed study of common minerals that carry REE and other HFSE from the magmatic through the hydrothermal stages, from six silica-saturated alkaline occurrences worldwide. The goal is to make out the meaning of textures and link them to the major events that impacted the selected complexes, i.e. the magmatic-hydrothermal transition. Once it is clear which zones of the crystals grew during the magmatic phase and which grew during hydrothermal circulations, we can analyze the REE and other HFSE compositions in the different zones, and infer the timing of concentration and fractionation of the REE.

Ideal candidates to do so are the pyroxene and amphibole mineral groups, because of their ubiquitous presence in these rocks and because they have been documented to form during a large span of the evolution of alkaline rocks, from early magmatic to late hydrothermal stages (e.g. [START_REF] Nielsen | The occurrence and formation of Ti-aegirines in peralkaline syenites[END_REF][START_REF] Jones | Zirconium-bearing aegirines from Motzfeldt, south Greenland[END_REF][START_REF] Ranløv | Compositional zoning in hydrothermal aegirine from fenites in the Proterozoic Gardar Province, South Greenland[END_REF][START_REF] Marks | Trace element variations in clinopyroxene and amphibole from alkaline to peralkaline syenites and granites: implications for mineral-melt trace-element partitioning[END_REF]). In addition, they incorporate a relatively high amount of REE and other HFSE compared to other rock-forming minerals (e.g. alkali feldspar and quartz), even though they are not considered as ore as their partition coefficients for HFSE (except Eu) are not high enough [START_REF] Fedele | Trace-element partitioning between plagioclase, alkali feldspar, Ti-magnetite, biotite, apatite, and evolved potassic liquids from Campi Flegrei (Southern Italy)[END_REF]. Primary zirconosilicates can also incorporate high levels of HFSE and are common in alkaline rocks, however, they are frequently altered to secondary minerals during the hydrothermal stage. Hence, pyroxenes and amphiboles are the best minerals to evidence the behavior of HFSE throughout the entire process of alkaline granite and pegmatite crystallization.

Geological background

The six alkaline igneous complexes on which this work focuses are: Ambohimirahavavy (Estrade et al. 2014a) and Manongarivo (Donnot 1963) in Madagascar, Amis in Namibia [START_REF] Schmitt | Zr-Nb-REE mineralization in peralkaline granites from the Amis Complex, Brandberg (Namibia): evidence for magmatic pre-enrichment from melt inclusions[END_REF], Evisa in Corsica [START_REF] Bonin | Age, origin and evolution of the anorogenic complex of Evisa (Corsica): A K-Li-Rb-Sr study[END_REF], Khan Bogd in Mongolia [START_REF] Kovalenko | Geology, geochronology, and geodynamics of the Khan Bogd alkali granite pluton in southern Mongolia[END_REF], and Strange Lake in Canada [START_REF] Currie | An unusual peralkaline granite near lac Brisson, Quebec-Labrador[END_REF].

These complexes were chosen because all present evidence of hydrothermal circulation and contain high amounts of REE (from 800 ppm in a granite from Khan Bogd to 12 % in a pegmatite from Strange Lake, unpublished data). All complexes contain peralkaline rocks, i.e. their ratio of (Na 2 O+K 2 O)/Al 2 O 3 is >1. Some of them have been the subject of in-depth studies, such as the Strange Lake complex which was first studied in the late 1980's (e.g. Miller, 1986;[START_REF] Salvi | The role of hydrothermal processes in the granite-hosted Zr, Y, REE deposit at Strange Lake, Quebec/Labrador: evidence from fluid inclusions[END_REF]) and still is today [START_REF] Vasyukova | Partial melting, fractional crystallisation, liquid immiscibility and hydrothermal mobilisation-A `Recipe'for the formation of economic A-Type granite-hosted HFSE deposits[END_REF]; others are poorly known, as is the case for the Manongarivo complex, for which the only available information is a mere mention in a French compilation (Donnot, 1963).

Although they were emplaced in different geodynamic complexes and at different periods, it appears that the melt source was a metasomatized upper mantle for all six complexes. Crustal contamination and feldspars fractionation, if occurring, are two processes that can enrich alkaline granites in REE. However, the rate of partial melting at the origin of these complexes is still debated, which is an essential parameter given that a lower rate implies a greater accumulation of incompatible elements, REE included, in the newly generated melt. All complexes underwent at least one alteration episode from an orthomagmatic fluid, and some of them an additional late hydrothermal event (identified at Amis, Evisa, and Strange Lake). The composition of the fluids that circulated was measured in fluid inclusions at Ambohimirahavavy and Strange Lake, and inferred from the study of secondary mineral assemblages for the other complexes. The presence of significant amounts of Na + and Fe 3+ in an orthomagmatic fluid is reported for Evisa, the two Malagasy complexes and Strange Lake, and Ca 2+ was identified at Khan Bogd, Madagascar and Strange Lake. Locally, CO 2 and/or CH 4 were reported at Khan Bogd and Strange Lake, and CO 3 2-at Amis. The resulting conclusions are similar for all the complexes: independently of its orthomagmatic or late origin, a F-, Ca-rich fluid is systematically observed and inferred to have mobilized and fractionated the REE. The importance of this process compared to magmatic ones in concentrating and fractionating the REE is still a matter of debate. Comparing six complexes that were emplaced in different geodynamic settings and witnessed different fluid conditions (composition, temperature), is a strategy to highlight common features in order to better constrain and generalize the timing of REE and other HFSE enrichment in alkaline complexes. Below, follows a rapid summary of the main geological features of each complex; more detailed information is provided in Online Resource 1.

Amis complex, Namibia

The Amis complex and the larger Brandberg complex onto which it is juxtaposed are of the same age (132.5 to 130.5 My) and are hosted in metasediments of the Damara orogen (~ 550 My;[START_REF] Miller | The Pan-African Damara Orogen of South West Africa/Namibia[END_REF]) and volcanic rocks from the Paraná-Etendeka igneous province. Three granite and 4 pegmatite-aplite samples were used for this study. The Amis complex is mainly made of an arfvedsonite granite and aegirine-rich pegmatite-aplites. The Amis complex is highly enriched in REE, other HFSE and volatiles (H 2 O, F). Pegmatites and aplites in the north-western part of the Amis complex are particularly enriched in REE mineralizations and uranium. [START_REF] Schmitt | Zr-Nb-REE mineralization in peralkaline granites from the Amis Complex, Brandberg (Namibia): evidence for magmatic pre-enrichment from melt inclusions[END_REF] infer that hydrothermal fluids had a very limited impact on the rocks of Amis complex. They affirm that minerals are all mainly magmatic, with the exception of replacement of arfvedsonite by quartz-hematite overgrowths. Based on the study of melt inclusions, they propose that the main process of REE and other HFSE enrichment is magmatic. In contrast, [START_REF] Diehl | A general equation for estimating Fe 3+ concentrations in ferromagnesian silicates and oxides from microprobe analyses, using stoichiometric criteria[END_REF] proposed that hydrothermal fluids played a significant role, and inferred the replacement of arfvedsonite by aegirine.

Evisa complex, Corsica

The Evisa complex (290 My; [START_REF] Cocherie | Comparative use of TIMS and SHRIMP for U-Pb zircon dating of A-type granites and mafic tholeiitic layered complexes and dykes from the Corsican Batholith (France)[END_REF]) was emplaced among metaluminous granitoids. We collected 7 granite and 5 pegmatite samples from the central part of the complex. The complex mainly comprises hypersolvus (perthitic feldspar) and subsolvus (albitic and alkali feldspar) peralkaline granites [START_REF] Bonin | Les granites des complexes annulaires[END_REF]) and associated pegmatites rich in REE-bearing minerals such as monazite-(Ce), apatite, and allanite-(Ce). The hypersolvus granite probably experienced two hydrothermal events: a minor, early episode at the end of crystallization, and a more important second episode, characterized by F-rich fluids, around 200 My [START_REF] Poitrasson | Importance of late-magmatic and hydrothermal fluids on the Sm-Nd isotope mineral systematics of hypersolvus granites[END_REF][START_REF] Bonin | Eurogranites-IGCP510 2008 Joint Field-meeting-Alkaline The Permian-Triassic A-type Volcanic-Plutonic Igneous Suite of Corsica[END_REF]. The subsolvus granite is intrusive into the perthitic unit and only experienced the second hydrothermal event, which triggered the replacement of the primary zirconosilicate elpidite by a secondary assemblage [START_REF] Bonin | Les granites des complexes annulaires[END_REF]).

Khan Bogd complex, Mongolia

The Khan Bogd complex is located in the southern Gobi Desert, at the transition between island-arc calc-alkaline differentiated volcanics (329±5 My) and rift-related bimodal basalt-comendite-alkali granite association (318-290 My), and is dated at approximately 290 Ma [START_REF] Kovalenko | Geology, geochronology, and geodynamics of the Khan Bogd alkali granite pluton in southern Mongolia[END_REF]. It is among the largest alkaline granite plutons in the world, with a surface area of 1500 km². We collected 3 granite and 7 pegmatite samples, mostly on the western side of the complex. The pluton consists of two ring bodies: a western peralkaline arfvedsonitebearing granite and a later emplaced, eastern aegirine granite. They are both associated with pegmatites. The pluton has been investigated for mining purposes and estimates indicate grades between 0.3 and 4.5 % REE, the highest grade being located at the top of the granite plutons (estimation of the tonnage is not available, [START_REF] Kovalenko | Endogenous rare metal ore formations and rare metal metallogeny of Mongolia[END_REF]. The origin of this pluton is linked to a mantle plume and the continental subduction of the South Mongolian Hercynides, but the exact magmatic process is still debated [START_REF] Kovalenko | Geology, geochronology, and geodynamics of the Khan Bogd alkali granite pluton in southern Mongolia[END_REF][START_REF] Kynicky | Distribution and evolution of zirconium mineralization in peralkaline granites and associated pegmatites of the Khan Bogd complex, southern Mongolia[END_REF]. [START_REF] Kynicky | Distribution and evolution of zirconium mineralization in peralkaline granites and associated pegmatites of the Khan Bogd complex, southern Mongolia[END_REF] propose the release of a silica-saturated orthomagmatic fluid that, similarly to Evisa, triggered the replacement of primary elpidite.

Manongarivo and Ambohimirahavavy complexes, Madagascar

Ambohimirahavavy and Manongarivo are the biggest complexes of the Ampasindava province in northwestern Madagascar, yet are not well documented. They were emplaced at circa 24 Ma into marine-shelf carbonates and marine-fluvial siliciclastic sediments of the Isalo Group [START_REF] Thomas | Geological evolution of the Neoproterozoic Bemarivo Belt, northern Madagascar[END_REF][START_REF] Cucciniello | The age and petrogenesis of alkaline magmatism in the Ampasindava Peninsula and Nosy Be archipelago, northern Madagascar | SpringerLink[END_REF]. We collected 3 granite and 4 pegmatite samples from the southern part of the Ambohimirahavavy complex. Due to poor outcropping conditions and low occurrence of pegmatites, we could only obtain 3 pegmatite samples from Manongarivo. The samples from Ambohimirahavavy were selected according to previous studies from Estrade et al. (2014a, b). Both complexes were likely emplaced in a rifting context and are made of two ring-shaped intrusions side by side, composed of nepheline syenite, alkali feldspar syenite, biotite granite, peralkaline pegmatites, and various volcanic rocks. Due to the recent discovery of rare metals in economic concentrations in ion-adsorption clays, the province drew the attention of a mining exploration company (Estrade et al. 2014b). Granite and pegmatite dykes are the most enriched lithologies in REE and other

HFSE. An orthomagmatic fluid metasomatized the complexes in the very last stages of magmatic evolution (Estrade et al. 2014b) and caused the replacement of primary zirconosilicates (eudialyte group minerals (EGM)

at Ambohimirahavavy, unknown at Manongarivo) by an assemblage of secondary HFSE-bearing minerals [START_REF] Lacroix | Extended Voigt-based analytic lineshape method for determining Ndimensional correlated hyperfine parameter distributions in Mössbauer spectroscopy[END_REF][START_REF] Rakotovao | Pétrologie du complexe alcalin cénozoïque d'Ambohimirahavavy, presqu'île d'Ampasindava, nord-ouest de Madagascar[END_REF][START_REF] Estrade | Crystallization and destabilization of eudialyte-group minerals in peralkaline granite and pegmatite: a case study from the Ambohimirahavavy complex, Madagascar[END_REF].

Strange Lake complex, Canada

The Strange Lake complex was emplaced 1240 My ago into Paleoproterozoic gneisses and quartz monzonite [START_REF] Miller | Structural and textural evolution of the Strange Lake peralkaline rare-element (NYF) granitic pegmatite, Quebec-Labrador[END_REF] and is considered to represent an extension of the Gardar rift in Greenland [START_REF] Pillet | Chronologie Rb/Sr et K/Ar du granite peralcalin du lac Brisson, Labrador central, Nouveau-Québec[END_REF][START_REF] Boily | The role of magmatic and hydrothermal processes in the chemical evolution of the Strange Lake plutonic complex, Quebec-Labrador[END_REF]Siegel et al. 2017a). We chose 3 granite and 4 pegmatite samples which were selected according to previous studies from Salvi and Williams- [START_REF] Salvi | The role of hydrothermal processes in the granite-hosted Zr, Y, REE deposit at Strange Lake, Quebec/Labrador: evidence from fluid inclusions[END_REF][START_REF] Salvi | Alteration, HFSE mineralisation and hydrocarbon formation in peralkaline igneous systems: Insights from the Strange Lake Pluton, Canada[END_REF]1995). The Strange Lake complex consists of a hypersolvus granite, while a younger, largely metasomatized transsolvus granite (rare perthite plus two distinct feldspars, previously classified as subsolvus, e.g. [START_REF] Boily | The role of magmatic and hydrothermal processes in the chemical evolution of the Strange Lake plutonic complex, Quebec-Labrador[END_REF] surrounds it, forming the majority of the complex [START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF]Siegel et al. 2017b). These plutons are partially bounded by an outwardly-dipping fracture associated with fluorite and hematite breccia. Two zones of pegmatites have been reported: the Main zone in the center of the complex, and the B zone on the northwestern edge. These magmatic events correlate with a progressive enrichment in REE and other HFSE, from HFSE-poor hypersolvus granite to HFSE-rich transsolvus granite and pegmatites [START_REF] Miller | Structural and textural evolution of the Strange Lake peralkaline rare-element (NYF) granitic pegmatite, Quebec-Labrador[END_REF] [START_REF] Gowans | Quest Rare Minerals Ltd.: Strange Lake Resource Estimation Grigor'eva AA, Zubkova NV, Pekov IV, et al ([END_REF]. The complex also contains significant amounts of ZrO 2 , Nb 2 O 5 , and BeO. Two main extensive hydrothermal events affected the complex. The first event is attributed to circulation of a hot (≥ 300 °C) orthomagmatic brine consisting of K-and Na-enriched aqueous and carbonic phases [START_REF] Salvi | Alteration, HFSE mineralisation and hydrocarbon formation in peralkaline igneous systems: Insights from the Strange Lake Pluton, Canada[END_REF][START_REF] Vasyukova | Closed system fluid-mineral-mediated trace element behaviour in peralkaline rare metal pegmatites: Evidence from Strange Lake[END_REF]. The second hydrothermal event consists of a cooler fluid (100 -200 °C) resulting from mixing of a meteoric fluid and a fluid originating from the granites [START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF].

Features common to pegmatites and granites from all six complexes

Pegmatites in all complexes are heterogeneous in texture, i.e. they show mineralogical layering visible from the scale of the outcrop to that of a thin section (Fig 3), whereas granites are generally homogeneous. In pegmatites, the layers have grain sizes ranging from a few µm to about 40 cm and locally more (Fig 3).

Therefore, to obtain meaningful data, samples for this study were selected from zones of relatively fine grain size (less than a few cm).

All pegmatites and granites are dominated by quartz, alkali feldspar (perthites, albite, orthoclase, and/or microcline) that are commonly strongly albitized, plus alkali amphiboles and, mostly in pegmatites, aegirine Ca and Na zirconosilicates, which are characteristic minerals of agpaitic rocks, are mostly eudialyte-group minerals (EGM) at Ambohimirahavavy [START_REF] Lacroix | Extended Voigt-based analytic lineshape method for determining Ndimensional correlated hyperfine parameter distributions in Mössbauer spectroscopy[END_REF][START_REF] Estrade | Crystallization and destabilization of eudialyte-group minerals in peralkaline granite and pegmatite: a case study from the Ambohimirahavavy complex, Madagascar[END_REF], and elpidite at Evisa, Khan Bogd and Strange Lake [START_REF] Bonin | Peralkaline granites in Corsica: some petrological and geochemical constraints[END_REF][START_REF] Salvi | The role of hydrothermal processes in the granite-hosted Zr, Y, REE deposit at Strange Lake, Quebec/Labrador: evidence from fluid inclusions[END_REF]Grigor'eva et al. 2011). These zirconosilicates were partially to totally replaced by secondary miaskitic mineral assemblages, which commonly form pseudomorphs [START_REF] Poitrasson | Importance of late-magmatic and hydrothermal fluids on the Sm-Nd isotope mineral systematics of hypersolvus granites[END_REF]Estrade et al. 2014b;[START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF]. The mineralogy of the pseudomorphs varies in the different complexes and consists either of different Zr-and/or REE-bearing minerals (Fig. 4b), or exclusively zircon plus quartz, i.e. typical miaskitic mineral assemblage (Fig 4c 4d). Pseudomorphs from Evisa consist of Ca-rich mineral assemblages such as fluorite and parisite-(Ce), the latter growing mostly at the expense of bastnäsite-(Ce). At Manongarivo and Amis, the precursor mineral has not been identified because partial replacement has not been observed in these complexes, so far.

Analytical methods

Polished thin sections were prepared from all samples from the six complexes and were studied using transmitted light microscopy, cathodoluminescence, Scanning Electron Microscopy (SEM), microprobe, and Laser-Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICPMS) in order to determine mineral textural relationships and amphiboles and pyroxene zoning.

SEM images in back-scattered electron (BSE) mode were collected with a Jeol JSM6360LV SEM coupled with a Bruker silicon drift detector and interpreted with the Esprit software at the Géosciences Environnement Toulouse laboratory, using an accelerating voltage of 20 kV, and a detection time of 40 s.

Mineral compositions were determined at the Centre Raymond Castaing in Toulouse, using a CAMECA SX Five electron microprobe (EPMA), with an accelerating voltage of 15 kV, a beam current of 20 nA for pyroxene, 10 nA for amphibole, and a beam diameter of about 2 μm. Quantitative elemental maps were performed on the same instrument, using an accelerating voltage of 15 kV, a beam current of 20 nA and 100 nA for major and minor elements respectively. Counting time was 1 s on each pixel with a step of 0.5 to 4 µm depending on crystal size. Uncertainties were calculated using the method of [START_REF] Ancey | Application des méthodes statistiques en microanalyse[END_REF]. Standards used for calibration as well as detection limits based on repeated measurements are provided in Online Ressource 2.

To resolve the oxidation state of iron measured by microprobe, calculations of the valence were made following Droop (1987). In order to verify these measurements, a 57 Fe Mössbauer spectrum was collected at the Laboratoire de Chimie de Coordination -University of Toulouse on a constant-acceleration conventional spectrometer, at 80 K and with a 1.85 GBq source of 57Co, over a ±12 mm/s range with 512 channels. The spectra were acquired via a compact detector c-system and recorded with a Canberra multichannel analyzer, coupled to a computer with Recoil Mossbauer Analysis Software 1.05 from Lagarec and Rancourt (1997).

Uncertainties were calculated using the covariance matrix. Pyroxene end-members were then calculated following [START_REF] Mann | Influence of oxygen fugacity on mineral compositions in peralkaline melts: The Katzenbuckel volcano, Southwest Germany[END_REF]. Amphiboles end-members were calculated using the spreadsheet available from [START_REF] Locock | An Excel spreadsheet to classify chemical analyses of amphiboles following the IMA 2012 recommendations[END_REF]. only 5 elements in order to obtain a better signal. LA-ICP-MS maps provide detailed spatial information on zonation in crystals that would not be available otherwise. Analysis were semi-quantitative and compared to microprobe data to obtain the real concentration.

To check for the absence of REE mineral inclusions in the pyroxene, we analyzed one sample from the Amis complex by Scanning Transmission Electron Microscopy (STEM) at the Centre Raymond Castaing in Toulouse. Unlike TEM, with STEM the electron beam is focused on a specific point of the sample and scans its surface instead of being focused on the entire surface. This allows a better resolution as well as local chemical identification (e.g. [START_REF] Pennycook | Scanning Transmission Electron Microscopy: Imaging and Analysis[END_REF]. The Focused Ion Beam (FIB) was used to prepare the sample for STEM. The FIB is a Hélios 600i operated at 5 kV, and the TEM a JEOL cold-FEG JEM-ARM200F operated at 200 kV equipped with a Cs probe corrector reaching a spatial resolution of 0.078 nm. Images were acquired in high-angle dark field (HAADF) as well as bright field (BF) modes. EDX spectra were recorded on a JEOL CENTURIO SDD detector.

Amphibole and pyroxene textural and chemical description

Textural features 4.1.1. Amphiboles

Amphiboles occur in granites and pegmatites as well-developed euhedral crystals (Fig 4e). No significant textural difference was found between amphiboles of the six complexes. They all appear pleochroic, dark green or dark blue to black under plain polarized light (PPL). Their size is highly variable from one complex to another as well as inside pegmatites (typically ranging from 20 µm to 10 cm), but is always comparable with the size of

Pyroxene

Pyroxene is ubiquitous in pegmatites and granites of all six complexes but present three different (3) In granites as well as in some pegmatites, anhedral pyroxene is found replacing amphibole (Fig 4e). This pyroxene is also weakly oscillatory zoned. As it replaces amphibole only partially, it's size is smaller than the previously described pyroxene crystals and depends on the original amphibole size. This pyroxene is associated with Fe-and/or Ti-oxides and can sometimes keep the 120° cleavages from the original amphibole.

Based on occurrence, size, color, zoning and mineralogical associations, we distinguish two main types of pyroxene in the six complexes we focus on in this study. Type-I matches euhedral zoned pyroxene, and our study focuses mainly on this type. Based on the observation that they present anhedral shape, oscillatory zoning, and are associated with Fe-and/or Ti-oxides, we regroup anhedral single pyroxene crystals and pyroxene replacing amphibole under the denomination type-II.

Major element composition

Amphiboles

Since they present the same textural features, all amphiboles were identified based solely on their chemistry, following the classification of [START_REF] Hawthorne | Nomenclature of the amphibole supergroup[END_REF] and in accordance with [START_REF] Leake | Nomenclature of amphiboles; report of the Subcommittee on Amphiboles of the International Mineralogical Association Commission on new minerals and mineral names[END_REF].

The general formula of amphibole is AB 2 C 5 T 8 O 22 W 2 where, in this study, A = Na; B = Na, Ca; C = Fe 2+ , Fe 3+ , Li; T = Ca, Al; and W = OH, F. The rootname is given to these amphiboles based on the composition of the sites A, B and C, and the prefix is added according to the composition of the sites C and W. By strictly following these rules, it appears that some complexes contain only one amphibole type, common to granites and pegmatites, and some contain different amphiboles (Table 1). Ambohimirahavavy (ferro-ferri-fluoro-leakeite; However, it is important to keep in mind that amphibole is a solid solution, thereby two amphiboles can have different names but be close in composition, and inversely can have the same name but present a compositional range. To avoid these artificial hiatuses and subsets, we focused on the concentration of elements involved, when attributing amphibole names in this study. Hence, the distinction between katophorite and other amphiboles is based on the Ca/(Ca+Na) ratio; being >0.25 (close to 0.4), this makes it a Na-Ca amphibole, compared to values <0.1 for other amphiboles, which are then considered Na-amphiboles. Thereby, ferro-ferrikatophorite at Evisa exhibits a significant compositional gap with the other amphiboles from this study. The distinction between Na-amphiboles is more subtle, and is based on the amount of Li and Fe 2+ /Fe 3+ ratio in the C site. Ideally, leakeite has 1 Li apfu; this value is on average 0.6 apfu at Ambohimirahavavy and Evisa, and 0.7 apfu at Amis. Arfvedsonites from the other complexes have average Li contents of 0.1-0.3 apfu, validating the distinction proposed here, although these values are intermediate between the two categories. The prefixes ferro-and ferri-refer to the presence of Fe 2+ and Fe 3+ , respectively, in the C site.The prefix fluoro-indicates that the F concentration is systematically higher than the concentration of other elements in the W site. In this study, the amounts of F and OH, respectively, are close to 1.1 and 0.9 apfu in arfvedsonite from all complexes but Khan Bogd, where it is closer to 0.9 apfu F and 1.1 apfu OH. Strictly complying with the definition would result in two different names (F and OH varieties), which would vary with each complex. However, this would lead to add unnecessary confusion to the nomenclature, therefore, because the compositional differences are only minor, we chose to ignore this parameter.

Within a given category, amphiboles are globally similar in composition in all complexes, including Fe 2+ /Fe 3+ (around 3 for fluoro-arfvedsonite, 2 for leakeite, and 3.5 for katophorite), with only local variations (Table 1). Chlorine was not detected by microprobe analyses in any of the amphiboles. The main distinction regarding major element composition is found in (fluoro)-arfvedsonite from Manongarivo and Khan Bogd, which both have high contents of Mn (respectively 1.92 and 1.66 % vs 0.50 % in the other complexes, Fig 5d)

and Fe 2+ /Fe 3+ (4.00 and 5.90 vs 2.40). Amphibole from Ambohimirahavavy is also rich in Mn (1.99 %) but does not show a high Fe 2+ /Fe 3+ value. Fluoro-arfvedsonite from Amis is as rich in F as ferro-ferri-fluoro-leakeite from all other complexes (around 1.6 apfu), which is unusual (Deer et al. 1997b). Amphiboles from pegmatites are generally slightly richer in HFSE, Li and Fe 3+ than amphiboles from granites.

Pyroxene

Based on composition, pyroxene in granites and pegmatites consists of 85 to 97 % aegirine endmember. Therefore, according to the classification of [START_REF] Morimoto | Nomenclature of Pyroxenes[END_REF], we will refer to it as aegirine in the following sections of this paper. Pyroxene from Ambohimirahavavy and Khan Bogd is more Ca-rich than the other complexes, but it is still classified as aegirine.

Type-I aegirine crystals show a characteristic zoning pattern in pegmatites from all locations: a core-to- Mg, Mn and Zr compared to rims, which in turn contain higher Al, Ti, Na and Fe 3+ . [START_REF] Shearer | Sector-zoned aegirine from the Ilimaussaq alkaline intrusion, South Greenland: Implications for trace-element behavior in pyroxene[END_REF] supplemented the description of aegirine from Ilímaussaq, by reporting Sr and REE enrichment in cores compared to rims. In our analyses, we did not observe a systematic difference in Mn and Sr contents between core and rims.

Type-II aegirine is not sector-zoned but shows oscillatory zoning and we observed a wide variety of compositions depending on which layer was analyzed (Table 2). This zoning is irregular, with no core, made of compositionally varying layers with variable thickness. No common pattern was found among any type-II aegirine crystals, beside their heterogeneity. Many elements are involved, including Fe, Mn, Ca, and Ti (see Online Resource 8). No significant chemical differences were found between type-II aegirine from granites and pegmatites. Type-II aegirine is commonly depleted in Ti and Sn compared to type-I aegirine, except at Khan Bogd. Mg is commonly found in higher concentration in type-II aegirine (Table 2). No significant compositional difference was found between type-II aegirine replacing amphibole in granite and those in pegmatites. At

Strange Lake, all observed aegirine crystals are of type-II (Table 2; Salvi and Williams-Jones 1990; Roelofsen 1997).

Trace element composition

Amphiboles

Amphiboles from all six complexes have very similar trace element patterns (Table 1). Fluoroarfvedsonite and ferri-fluoro-leakeite have practically the same composition, with Li contents alone being responsible for the transition from one kind to the other. This cation enters the amphibole structure according to the substitution: C Li + + C Fe 3+ = 2 C Fe 2+ , described by [START_REF] Hawthorne | Fluor-ferro-leakeite, NaNa2 (Fe2+ 2Fe3+ 2Li) Si8O22F2, a new alkali amphibole from the Canada Pinabete pluton, Questa, New Mexico, USA[END_REF]. Ferro-ferri-katophorite, which is more calcic than the two other amphiboles, relates to fluoro-arfvedsonite through the substitution: B Ca 2+ + T Al 3+ = B Na + + T Si 4+ (Deer et al. 1997b). This last substitution increases the size of the C site, leaving more space to incorporate LREE. Similarly, the coupled substitution (Mg, Fe 2+ )  Fe 3+ occurring in the C site of amphiboles enhances HREE incorporation [START_REF] Bottazzi | Distinct site preferences for heavy and light REE in amphibole and the prediction of Amph/L D REE[END_REF]Siegel et al. 2017b) [START_REF] Hawthorne | Li-bearing arfvedsonitic amphiboles from the Strange Lake peralkaline granite, Quebec[END_REF]. Mn and Zn covary in amphiboles from Ambohimirahavavy and Amis as well, but in a different way (respectively 4:1 and 1:3), while at Manongarivo and Khan Bogd they do not appear to be related.

Pyroxene

Type-I aegirine, similarly to amphiboles, presents comparable trace elements patterns but different absolute concentrations, from one complex to the other (Table 2). Cores of type-I aegirine from the two Malagasy complexes have high Zr concentrations compared to other complexes (> 8000 vs 2700 ppm on average for the other complexes). Khan Bogd aegirine is poorer in Li (17 vs 115 ppm), LREE (12 vs 80 ppm on average) and HREE (21 vs 125 ppm on average). Type-I aegirine from Amis, similarly to ferro-ferri-fluoro-leakeite in this complex, shows quite a different trace chemistry than aegirine from other complexes. Average Li is low (50 vs 115 ppm), while other metals are high, e.g. Zn (1780 vs 435 ppm), Sn (3740 vs 700 ppm), Pb (39 vs 7 ppm) and U (29 vs 0.5 ppm) (Table 2 andFig To verify that the concentrations of REE measured in aegirine truly result from their incorporation in the crystal structure and not from small mineral inclusions, we used a STEM approach to investigate the core of a type-I aegirine, where REE are most concentrated. We did not detect any REE mineral inclusion, validating the premise that the REE are indeed incorporated in the pyroxene's structure. The main substitution involved in this process is common to type-I aegirine from all pegmatites. The classic aegirine substitution is M2 Ca 2+ + M1 (Mg 2+ , Fe 2+ ) ↔ M2 Na + + M1 Fe 3+ . (Eq. 1)

However, in aegirine from this study there is not enough Fe 3+ to equilibrate this reaction, hence Al, Zr and Ti, present in aegirine, also contribute to compensate this deficiency, resulting in a new substitution mechanism:

M2 Ca 2+ + M1 (Mg 2+ , Fe 2+ ) ↔ M2 Na + + M1 (Fe 3+ , Al 3+ , Zr 4+ , Ti 4+ ). (Eq. 2)

However, REE can enter both M1 and M2 sites along with Fe 3+ and Na + [START_REF] Beard | Clinopyroxene/melt trace element partitioning in sodic alkaline magmas[END_REF]. In clinopyroxenes, incorporation of REE, as well as other HFSE, is controlled by five crystal-chemistry driven mechanisms, with site parameters being determined by the size of the major elements. The first mechanism is decreasing temperature and increasing melt alkalinity, which usually lowers the REE compatibility. The second mechanism is similar to that for amphiboles and is defined by the substitution Ca → Na, which increases the size of the M2 site leaving more space to incorporate LREE. It also decreases the size of the M1 site, making it more adequate to incorporate HREE. The third mechanism takes place in the M1 crystallographic site, where the combined substitution (Mg, Fe 2+ ) → Fe 3+ leads to a charge effects that allows more HREE to enter this site along with Fe 3+ [START_REF] Marks | Trace element variations in clinopyroxene and amphibole from alkaline to peralkaline syenites and granites: implications for mineral-melt trace-element partitioning[END_REF][START_REF] Beard | Clinopyroxene/melt trace element partitioning in sodic alkaline magmas[END_REF]). The fourth mechanism, described by [START_REF] Mollo | Clinopyroxene-melt element partitioning during interaction between trachybasaltic magma and siliceous crust: Clues from quartzite enclaves at Mt. Etna volcano[END_REF], is also due to a charge effect. These authors propose the coupled substitution M1 (Al, Fe 3+ ) + T Al ↔ M2 (Mg, Fe 2+ ) + T Si, (Eq. 3) whereby the increase in Al at the expense of Si creates a charge deficiency compensated by the incorporation of REE in the M2 site. In addition, the growth from an environment rich in REE can result in incorporation of high levels of REE in aegirine. Indeed, the adsorption of REE onto aegirine surface during supersaturation may be too fast to equilibrate as crystal growth rate exceeds the internal diffusion rate [START_REF] Smith | The rare earth elements and uranium in garnets from the Beinn an Dubhaich Aureole, Skye, Scotland, UK: constraints on processes in a dynamic hydrothermal system[END_REF].

All of the above mechanisms have the potential of favoring the incorporation of REE into pyroxene, however, aegirine in our study having low T Al and practically no Mg or Fe 2+ in the M2 site, the mechanism of Mollo et al. is probably not significant in this case. The high REE concentration in aegirine of this study is therefore probably the result of the four other mechanisms that rely on size, charge, and kinetics effect at M sites of aegirine.

Ce anomalies

Some amphibole and aegirine chondrite-normalized spectra, mostly from one pegmatite sample from Amis, display a peak in Ce concentration (Fig 8 and9). The occurrence of this peak being restricted to a few samples, mass interference is excluded. In addition, ICP-MS signals showing peaks during ablation were all excluded so the presence of this Ce peak is not due to mineral inclusions. This peak may be due to the presence of Ce 4+ instead of the usual Ce 3+ . Indeed, Ce 4+ has a smaller ionic radius (0.97 Å) than Ce 3+ (1.143 Å) for an eight-fold coordination, thereby Ce 4+ substitutes more easily with Fe 2+ (0.92 Å) and Mg 2+ (0.89 Å; Shannon 1976). The change in oxidation state is commonly interpreted to indicate an increase in oxygen fugacity in the environment that can be either low temperature magmatic (<600 °C, [START_REF] Trail | Ce and Eu anomalies in zircon as proxies for the oxidation state of magmas[END_REF] or hydrothermal [START_REF] Xu | Origin of heavy rare earth mineralization in South China[END_REF].

Discussion

Origin of amphiboles

In the previous sections we have seen that three different amphiboles can be distinguished based on Na rises in the melt as differentiation occurs, until Li-and Na-rich amphiboles crystallize instead of their Lipoor, Mg-, Ca-rich equivalents. Based on various pieces of evidence, such as having a similar age to that of granite emplacement [START_REF] Kovalenko | Geology, geochronology, and geodynamics of the Khan Bogd alkali granite pluton in southern Mongolia[END_REF], presence of melt inclusions [START_REF] Schmitt | Zr-Nb-REE mineralization in peralkaline granites from the Amis Complex, Brandberg (Namibia): evidence for magmatic pre-enrichment from melt inclusions[END_REF], occurrence as well-developed euhedral crystals, and dissolution/corrosion textures (Estrade et al. 2014a;[START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF], it is accepted in the literature that amphiboles are magmatic phases in peralkaline granites. Hence, in accordance with our observations and those in the literature about the six studied complexes, we conclude that amphiboles of all types in our samples are magmatic. Their variation in trace element composition is directly linked to the environment they grew in, a highly evolved melt, enriched in incompatible elements. Among incompatible elements are the REE, which are known to be preferentially incorporated in alkaline amphiboles along with Na for LREE and Fe 3+ for HREE [START_REF] Bottazzi | Distinct site preferences for heavy and light REE in amphibole and the prediction of Amph/L D REE[END_REF]Siegel et al. 2017b). Due to this crystallo-chemical control, amphiboles can incorporate a lot of REE at the magmatic stage, with a preference for HREE. Slight compositional differences in REE, Zn, Sn, Nb, Zr and Pb between amphiboles of the 6 complexes (Fig 5d) likely reflect parental melt compositions. For example, the presence of ferro-ferri-katophorite in Evisa probably indicates an effect of melt contamination in Ca by the surrounding aluminous granites.

As previously described, amphiboles are commonly altered to Fe-and Ti-oxides and to type-II aegirine.

Based on textural evidence that at Strange Lake this replacement occurs preferentially along cleavages in altered granites and pegmatites, as well as on mass-balance calculations, Salvi and Williams [START_REF] Salvi | The role of hydrothermal processes in the granite-hosted Zr, Y, REE deposit at Strange Lake, Quebec/Labrador: evidence from fluid inclusions[END_REF] identified this pyroxene occurrence as hydrothermal. Because we observed this replacement in all of our samples, it follows that a late hydrothermal event is a common feature in the complexes studied and, by extrapolation, most likely in alkaline rocks in general; indeed, a hydrothermal event was also documented in silica-undersaturated alkaline rocks such as at Ilímaussaq (e.g. [START_REF] Borst | Zirconosilicates in the kakortokites of the Ilímaussaq complex, South Greenland: Implications for fluid evolution and high-field-strength and rare-earth element mineralization in agpaitic systems[END_REF] or at Khibina and Lovozero (e.g. [START_REF] Arzamastsev | Contact interaction of agpaitic magmas with basement gneisses: an example of the Khibina and Lovozero massifs[END_REF]. [START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF] and [START_REF] Vasyukova | Closed system fluid-mineral-mediated trace element behaviour in peralkaline rare metal pegmatites: Evidence from Strange Lake[END_REF] show that at Strange Lake this event is due to a relatively high temperature, Ca-rich fluid. Based on the similarity of amphiboles alteration in all six complexes, we infer that the circulation of an orthomagmatic fluid is a process common to all complexes.

Origin of pyroxene

Core-to-rim zonation

Even though both types of aegirine identified in all complexes show some kind of zonation, only type-I

shows a core-to-rim zonation, with a core enriched in Ca, Zr, Sn, Hf, Sc, REE and depleted in Na, Al and Fe compared to the rims (Fig 6). The rims also have a higher Fe 3+ / Fe 2+ ratio than the core. This pattern is systematic for type-I aegirine crystals from all complexes. Type-II aegirine is oscillatory zoned and can replace amphibole.

The sharp compositional change between the core and the rims in type-I aegirine ( ). This process increases the amount of aegirine endmember compared to hedenbergite endmember which, along with the hydrothermal environment being more oxidizing than the silicate melt, increases the Fe 3+ / Fe 2+ ratio in the rims.

The absolute amount of Al is very low in aegirine of this study, so the observed increase in this element in aegirine rims can be linked only to charge effects, with Al 3+ being incorporated along with Fe 3+ .

Because of their similar behavior, Zr and Hf both occur in the same growth zones. Based on the many arguments in favor of a magmatic core and hydrothermal rims, we propose two hypotheses as to the crystallization timing of type-I aegirine. [START_REF] Jones | Zirconium-bearing aegirines from Motzfeldt, south Greenland[END_REF] asserted that formation of Zr-enriched aegirine is enhanced by low oxygen fugacity, and can only occur if no other Zr-bearing phase is growing. Pseudomorphs after primary zirconosilicates were found around aegirine crystals (Fig 4b), hence it is impossible that aegirine grew after the magmatic zirconosilicates. Our first hypothesis is that, in accordance with the observations made by [START_REF] Jones | Zirconium-bearing aegirines from Motzfeldt, south Greenland[END_REF], type-I aegirine cores grew before magmatic zirconosilicates, in the magma chamber. In the absence of any mineral Zr has a strong affinity for, Zr partitioned mostly into aegirine. These primary phenocrysts were brought into pegmatites, where primary zirconosilicates began to crystallize, incorporating most of the available Zr.. Moving on to the hydrothermal stage, primary zirconosilicates were destabilized and replaced by minerals richer in Zr, including zircon, which form the pseudomorphs. Aegirine rims were then in competition with zircon, which has a high Zr uptake and therefore accounts for the decrease in Zr concentration in aegirine rims. In this hypothesis, aegirine necessarily stopped growing during their migration from the magma chamber to the pegmatites, or it would have resulted in another zonation in the magmatic core.

Our second hypothesis is that aegirine core and primary zirconosilicates grew at the same time. In this case the conclusions of [START_REF] Jones | Zirconium-bearing aegirines from Motzfeldt, south Greenland[END_REF] do not apply. This may be because their study was based on rocks from the magmatic chamber of Ilímaussaq which had plenty of time to equilibrate during crystallization, whereas our study is based on pegmatites, which crystallize quickly. Anyhow, in this hypothesis the competition between the two minerals resulted in the preference of Zr to go into primary zirconosilicates, but because of the high content of Zr in the melt and its low mobility [START_REF] Duggan | Zirconium-rich sodic pyroxenes in felsic volcanics from the Warrumbungle Volcano, Central New South Wales, Australia[END_REF], affinity for aegirine was still high enough to enrich it up to the 2 wt% that we have measured. The subsequent hydrothermal stage is similar to the one in our previous hypothesis, with higher Zr uptake by zircon rather than by zirconosilicates accounting for the decrease in Zr concentration in aegirine rims. A similar behavior was observed for Sn and is most likely due to a similar phenomenon. No Sn-based mineral was found in pseudomorphs, but Sn 4+ can easily enter in the composition of secondary Ti-oxides because its radius is similar to the Ti 4+ ion in eight-fold coordination [START_REF] Shannon | Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides[END_REF]). In

Sector zoning

Sector zoning occurs in the most outer rims in type-I aegirine, involves variations in Ti, Ca, Fe and REE, and was observed in all complexes. Many explanations have been invoked as to the origin of sector zoning, and the latest agree on the fact that it is linked to a crystal growth faster than element diffusion rates [START_REF] Ubide | Sector-zoned clinopyroxene as a recorder of magma history, eruption triggers, and ascent rates[END_REF]. If this condition is respected, it will result a disequilibrium between the different crystal faces that will not all incorporate the same elements as it normally would [START_REF] Strong | Formation of the hour-glass structure in augite[END_REF][START_REF] Ferguson | On hour-glass sector zoning in clinopyroxene[END_REF][START_REF] Larsen | Sector zoned aegirine from the Ilímaussaq alkaline intrusion, South Greenland[END_REF]. Elements are however incorporated on the different crystal faces in a matter so that charge balance is respected [START_REF] Ubide | Sector-zoned clinopyroxene as a recorder of magma history, eruption triggers, and ascent rates[END_REF]. The crystallization rate could also impact the size of protosites (partially formed site on the surface of a growing crystal), which are in equilibrium with the crystallizing environment, and preferentially incorporate elements depending on their charge/radius ratio. If the growth rate is higher than ionic diffusion in the crystal, the protosites will not have enough time to re-equilibrate in size and incorporate different elements on the different faces of the crystal [START_REF] Nakamura | Origin of sector-zoning of igneous clinopyroxenes[END_REF]). This higher growth rate is likely related to a lower temperature due to the change from magmatic to hydrothermal environment [START_REF] Barkov | Anomalous Cr-rich zones in sector-zoned clinopyroxene macrocrysts in gabbro[END_REF]. However, these considerations must be taken carefully, as most of the research is based on augite and not aegirine.

Incorporation of Ti in aegirine is uncommon [START_REF] Rønsbo | Titan-aegirine from early Tertiary ash layers in northern Denmark[END_REF][START_REF] Dyulgerov | Unusual Ti and Zr aegirine-augite and potassic magnesio-arfvedsonite in the peralkaline potassic oversaturated Buhovo-Seslavtzi complex, Bulgaria[END_REF]. Many explanations have been evoked: high temperature, low fO 2 , high fO 2 , low pressure, high activities of TiO 2 and SiO 2 , high activity of alkalis, or a peculiar chemical composition of the environment such as low Zr and Na and high Ca [START_REF] Ferguson | On hour-glass sector zoning in clinopyroxene[END_REF][START_REF] Flower | Phase relations of titan-acmite in the system Na2O-Fe2O3-Al2O3-TiO2-SiO2 at 1000 bars total water pressure[END_REF][START_REF] Larsen | Clinopyroxenes and coexisting mafic minerals from the alkaline Ilimaussaq intrusion, South Greenland[END_REF][START_REF] Rønsbo | Titan-aegirine from early Tertiary ash layers in northern Denmark[END_REF][START_REF] Nielsen | The occurrence and formation of Ti-aegirines in peralkaline syenites[END_REF][START_REF] Brousse | Crystallization trends of pyroxenes from agpaitic phonolites (Cantal, France)[END_REF][START_REF] Dyulgerov | Unusual Ti and Zr aegirine-augite and potassic magnesio-arfvedsonite in the peralkaline potassic oversaturated Buhovo-Seslavtzi complex, Bulgaria[END_REF]. In our samples, a high temperature is ruled out as Ti incorporation occurs during the hydrothermal stage. As described in the previous section, Zr and Ca concentration are low in the crystallizing environment, and Na is high. Hence, and in accordance with [START_REF] Dyulgerov | Unusual Ti and Zr aegirine-augite and potassic magnesio-arfvedsonite in the peralkaline potassic oversaturated Buhovo-Seslavtzi complex, Bulgaria[END_REF], [START_REF] Larsen | Clinopyroxenes and coexisting mafic minerals from the alkaline Ilimaussaq intrusion, South Greenland[END_REF][START_REF] Flower | Phase relations of titan-acmite in the system Na2O-Fe2O3-Al2O3-TiO2-SiO2 at 1000 bars total water pressure[END_REF], we suggest that in our rocks the incorporation of Ti into hydrothermal aegirine is favored by a high activity of Na and the low availability of Zr. Titanium may have been provided by the alteration of magmatic Ti-bearing minerals, such as astrophyllite (found at Manongarivo, Strange Lake, Khan Bogd and Amis, Fig 4a) chevkinite-(Ce) (found at Ambohimirahavavy, Manongarivo and Evisa), or others such as narsarsukite, låvenite, and aenigmatite.

Oscillatory zoning

Oscillatory zoning is present in both types of aegirine. In type-I, it is found in the sector zoned rims C (crystallographic sector (110)) and D (crystallographic sector (100)) as well as in zone B, intermediate between the core and the rims. Considering that zone B in type-I aegirine is only present in our samples in Khan Bogd, and was locally described in aegirine from the South Gardar Province, South Greenland by [START_REF] Ranløv | Compositional zoning in hydrothermal aegirine from fenites in the Proterozoic Gardar Province, South Greenland[END_REF], a general process to explain its formation is difficult to establish. Like zones C and D, zone B is poorer in Ca, Zr, Hf, Sn and REE than zone A. It is not specifically enriched in any element, but shows oscillatory zoning in Ti and Mn. We believe that zone B development is related to a specificity in the history of the Khan Bogd complex [START_REF] Kynicky | Distribution and evolution of zirconium mineralization in peralkaline granites and associated pegmatites of the Khan Bogd complex, southern Mongolia[END_REF]. In this complex, magmatic zone A aegirine would crystallize before or in competition with primary elpidite. At the hydrothermal stage, primary elpidite altered mostly into a secondary Type-II aegirine either replaces amphibole or forms isolated crystals. Both kinds are found associated with Fe-and/or Ti-oxides, and are randomly zoned in many elements (Ti, Ca, Mn, Fe, Hf, Sn, Zr at least). The formation of aegirine from amphibole is known to be hydrothermal (e.g. [START_REF] Salvi | The role of hydrothermal processes in the granite-hosted Zr, Y, REE deposit at Strange Lake, Quebec/Labrador: evidence from fluid inclusions[END_REF][START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF]. For arfvedsonite, it follows the reaction provided by Salvi and Williams-Jones (1997):

3Na 3 Fe 5 Si 8 O 22 (F,OH) 2 + 2H 2 O = 9NaFeSi 2 O 6 + 2Fe 3 O 4 + 6SiO 2 + 5H 2 .
(Eq. 4)

The released Fe generally forms hematite or magnetite. Considering that zoning patterns are the same for all type-II aegirine crystals and that they are systematically found associated with Fe-and/or Ti-oxides, we consider that they all grew replacing amphibole, totally or partially. Hence, we consider type-II aegirine as hydrothermal.

REE behavior in amphiboles and aegirine

As mentioned in the introduction, amphiboles and pyroxenes, in agreement with their crystal-chemical parameters, can incorporate up to several hundreds of ppm of REE. In addition, their alkaline nature enhances the fractionation of light versus heavy REE through the preferential incorporation of HREE in their structure, resulting in the globally HREE-enriched patterns observed (Fig. 9; [START_REF] Beard | Clinopyroxene/melt trace element partitioning in sodic alkaline magmas[END_REF]. From the amphibole compositions, we know that the melt was rich in F (Table 1). [START_REF] Beard | The effect of fluorine on clinopyroxene / melt trace-element partitioning[END_REF] showed that the more F in the melt, the lower the pyroxene-melt partition coefficient, hence the less REE are incorporated in pyroxene; this affects the LREE more than the HREE and hence contributes to the observed globally HREE-enriched patterns.

The REE compositional patterns observed in type-I aegirine are similar to the compositional trends for Zr, Hf and Sn, i.e. they are globally enriched in the core compared to the rims, except for Khan Bogd where the opposite is true (Fig 10, Table 3). Given that the REE concentrations in aegirine at Khan Bogd are very low (2-7 ppm for Ce, 4-15 ppm for Y), we believe that this reverse trend is actually not significant (Fig 10) site. However, this mechanism is well defined only for pyroxenes whose aegirine component is not higher than 50 %. In our samples, the substitution Ca → Na simply cannot be responsible for the zonation observed, as more REE are found in the C sector, where Ca concentration is also higher. In addition, in Ca-rich pyroxenes such as diopside and augite, Ca was not reported to fractionate between sectors (e.g. [START_REF] Ubide | Sector-zoned clinopyroxene as a recorder of magma history, eruption triggers, and ascent rates[END_REF], which is another argument for the occurrence of a different exchange mechanism in our samples.

In addition to the above variations in bulk abundance, cores and rims also preferentially incorporate REE with lower or higher atomic numbers, creating a fractionation between the heavy and light REE. To evaluate this process, we calculated the percentages of REE enrichment or depletion from core to rim in type-I pyroxene. To do so, data obtained by LA-ICPMS were assigned either to the core or to the rim, based on their Zr content and zoning SEM observations. By comparing these two parameters, we were able to determine a threshold value, unique to each complex, below which a given aegirine analysis should be considered as rim, and above which core. These data were averaged out to obtain a single value each for core and rim (Fig 10). Finally, we calculated the extent of loss or enrichment between core and rim for each REE, according to: From the calculation of the percentage of depletion of REE from core to rim (Fig 11) we see that Amis, Manongarivo and Evisa have a trend with a more important depletion of LREE than HREE: from 40-80 % to 20-40 % depletion, respectively. Ambohimirahavavy shows the opposite trend: around 30 % of LREE depletion for 45% of HREE depletion. The latter trend is not as obvious as for other complexes, but it is still significant.

General model

We propose a global model of formation for amphiboles and aegirine in pegmatites (Fig 12). During the magmatic stage, two hypotheses remain. In the first one, amphibole and zone A of type-I aegirine grew at early magmatic stage, concentrating Zr and REE. In the second hypothesis, type-I aegirine core, amphibole and zirconosilicates grew in competition. Most Zr was incorporated in the zirconosilicates, as well as REE if the primary zirconosilicates was EGM. However, highly differentiated alkaline melts being particularly enriched in Zr and REE, there is an extra supply in these metals which, combined with a fast crystallization, could be incorporated by the structures of amphiboles and aegirine (Fig 12a).

During the orthomagmatic hydrothermal stage, amphiboles did not grow anymore, and type-I aegirine was in competition with secondary Zr-and REE-bearing minerals (generally zircon and bastnäsite-(Ce)) that formed pseudomorphs after primary zirconosilicates. Since HFSE have much higher affinity for these secondary minerals, their concentrations globally decreased in the co-crystallizing aegirine rims. It is also during this stage that type-II aegirine crystals nucleated. Nevertheless, competition does not explain the observed LREE-HREE fractionation. Residual melt, depleted in incompatible elements, likely crystallized late quartz and feldspar.

Given the overwhelming evidence for circulation of an orthomagmatic fluid in all complexes, we propose that fractionation results from variations in fluid composition. From the presence of secondary minerals such as bastnäsite-(Ce) and fluorite, and the study of fluid inclusions at Ambohimirahavavy and Strange Lake [START_REF] Estrade | Le complexe cénozoïque alcalin d'Ambohimirahavavy à Madagascar : origine, évolution et minéralisations en métaux rares[END_REF][START_REF] Vasyukova | Fluid evolution in the Strange Lake granitic pluton, Canada: Implications for HFSE mobilisation[END_REF], we know that the fluids contained at least the anions Cl -, F -, and CO 3 2-. These species have the ability to form stable complexes with the REE, to different extents. F -and CO 3 2-together are mostly depositional ligands [START_REF] Migdisov | Hydrothermal transport, deposition, and fractionation of the REE: Experimental data and thermodynamic calculations[END_REF]), thus precipitate with REE typically forming the LREEbearing bastnäsite-(Ce). On the other hand, Cl -is a weak but abundant ligand that can transport LREE more easily than the heavy ones [START_REF] Migdisov | Hydrothermal transport, deposition, and fractionation of the REE: Experimental data and thermodynamic calculations[END_REF]. If present as a transporting ligand, F -also transports more LREE than HREE [START_REF] Beard | The effect of fluorine on clinopyroxene / melt trace-element partitioning[END_REF]. We can thus infer that, depending on the relative concentrations of each ligand, fractionation of light vs heavy REE would be specific to each complex (

Conclusions

By comparing amphibole and pyroxene in six alkaline complexes worldwide, our study highlights that despite different geodynamic contexts, mineralogy, and REE enrichment, there is a general crystallization process common to all studied alkaline pegmatites. We believe this process can be generalized to all alkaline SiO 2 -saturated pegmatites worldwide, with local variations in mineralogy and REE-enrichment and fractionation rates. This crystallization process includes at least two main stages, respectively magmatic and hydrothermal.

From the study of the composition of amphibole and type-I aegirine core, we document the magmatic stage, during which a systematic primary enrichment in REE, HFSE and a LREE-HREE fractionation occurs. This enrichment is common to all complexes and linked to the properties of the extremely differentiated alkaline magma. The enrichment rate, however, differs from one complex to another. The REE fractionation is most likely due to crystallographic properties of aegirine and alkali amphiboles that preferentially incorporate HREE and do not represent an inherent property of the melt.

From the study of the texture and composition of type-I aegirine rims and type-II aegirine, we document similarities between complexes during the hydrothermal stage, at which amphiboles stop growing, type-I aegirine grows sector-zoned rims, and type-II aegirine begins replacing amphibole along with Fe-and Ti-oxides.

The drop in REE concentration in hydrothermal aegirine of all complexes shows that REE hydrothermal enrichment is globally less effective than magmatic REE enrichment, but still effective. The observed REE fractionation in all complexes show that hydrothermal fluids also have the ability to fractionate LREE-HREE through their transportation associated with ligands such as F -, Cl -. The rate of fractionation as well as the amount of mobilized REE depends on the local properties of the circulating fluid(s). These observations show that hydrothermal flow is mandatory in order to enhance the concentration of the REE, and specially HREE, in alkaline pegmatites to ore levels.

Table 1 Mean values for the major and trace element composition, obtained by microprobe and LA-ICPMS respectively, of amphiboles from pegmatites and granites from the six complexes. Abbreviations: no.: number of replicates; F-arf: fluoro-arfvedsonite; F-f-f-lea: ferro-ferri-fluoro-leakeite; F-f-kato: ferro-ferri-katophorite; b.d.l: below detection limit; Ambo: Ambohimirahavavy; KB: Khan bogd; M: Manongarivo; SL: Strange Lake

  placer deposits and are commonly hosted in fluorcarbonates such as bastnäsite-(Ce), parisite-(Ce) and synchysite-(Ce), and phosphate such as monazite-(Ce). Ion-adsorption deposits provide most of the global HREE production (Chakhmouradian and Wall 2012; U. S. Geological Survey 2019). In alkaline igneous rocks, REE are hosted mainly in monazite-(Ce), xenotime-(Y), fergusonite-(Y), loparite-(Ce) and eudialyte-group minerals (EGM)

(

  Fig 4). Common minerals in both rock types include zircon, pyrochlore group minerals (PGM) and Fe-and Tioxides. Other minerals present mostly in pegmatites include Ca and Na zirconosilicates, fluorite, astrophyllite and REE-bearing minerals (typically bastnäsite-(Ce), monazite-(Ce), xenotime-(Y), chevkinite-(Ce), allanite-(Ce), aeschynite-(Y), fergusonite-(Y), britholite-(Y), synchysite-(Ce)). Ti-bearing minerals such as chevkinite-(Ce) and astrophyllite show signs of alteration such as dissolved edges (Fig 4a).

  Trace element concentrations were determined in situ by LA-ICPMS at the Géosciences Environnement Toulouse laboratory, University of Toulouse, on 30-μm polished thin-sections, using a New Wave Research ESI 213 laser coupled to a Thermo-Fisher Element-XR high-resolution ICP-MS. Laser beam diameter varied from 30 to 50 µm depending on mineral size, and data were acquired by ablating lines instead of spots in order to prevent intersecting solid inclusion or internal zoning. NIST synthetic glass certified reference materials SRM 610 and SRM 612 were used as external and quality control standards, respectively. Each analysis was normalized using SiO 2 values previously determined by EMPA. The relative precision ranged from 5 to 10% and detection limits, based on repeated measurements, are provided in Online Ressource 3. Trace elements maps were performed at the Trinity College in Dublin according to the method provided in[START_REF] Ubide | High-resolution LA-ICP-MS trace element mapping of igneous minerals: In search of magma histories[END_REF] using photon machines G2 193 nm UV laser with a Helex two-volume cell coupled to a Thermo iCAPQ ICP-MS. Laser beam spot for maps was 12 µm square, scan speed 8 µm/s and repetition rate 41 Hz. NIST SRM 610 glass was used as standard. A beam size of 12 µm square typically allows a spatial resolution of 7-10 µm in x and y. This relatively new ablation technique removes only 3-4 µm of material, which allows several maps to be performed on the same area. A first map with all REE and other HFSE was first established, before acquiring a second map with

  quartz and alkali feldspar crystals in the same rock. Except at Evisa(Fig 4j), amphiboles are not zoned, but they display dissolution/corrosion textures on their rims. They commonly have small quartz and feldspar inclusions, and are in equilibrium with these minerals. Replacement of amphiboles by aegirine and Ti-and/or Fe-oxides is a common feature of peralkaline granitoids[START_REF] Marks | A global review on agpaitic rocks[END_REF] and occurs in all pegmatites and some granites (Fig4e).

  textural features. (1) In pegmatites which contain only minor amphibole, it occurs as light to dark green euhedral crystals under PPL ranging from 50 µm to 10 cm. These pyroxene crystals present a core-to-rim zonation pattern, with additional sector zoning in the rims, which are visible under PPL(Fig 4f). Similarly to amphiboles, these crystals are in equilibrium with quartz and alkali feldspar grains. In the rims of this pyroxene, we also found fluid inclusions(Fig 4g) and mineral inclusions that match the composition of those met in pseudomorphs (Fig4b and h). These pyroxene crystals were found in all studied complexes but Strange Lake. At Khan Bogd, these crystals present an additional zoning area between the core and the sector zoned rims(Fig 4k). (2) In some pegmatites, pyroxene occurs as anhedral crystals of size similar to the one we just described. These crystals are oscillatory zoned, from yellowish to green under PPL(Fig 4i). They commonly include Fe-and/or Ti-oxides.

Na 3

 3 Fe 4 LiSi 8 O 22 F 2 ), Manongarivo and Strange Lake (fluoro-arfvedsonite; Na 3 Fe 5 Si 8 O 22 F 2 ), and Khan Bogd (arfvedsonite; Na 3 Fe 5 Si 8 O 22 (F,OH) 2 ) contain only one amphibole. At Amis, amphibole is ferro-ferri-fluoroleakeite in pegmatites and fluoro-arfvedsonite in the granite. At Evisa, granites and pegmatites contain a zoned amphibole with a core of fluoro-arfvedsonite and a rim of ferro-ferri-fluoro-leakeite (Fig 4L). In addition, (Fe 2+ ,Mg) 4 Fe 3+ (Si 7 Al)O 22 (OH) 2 ) in their core and fluoro-arfvedsonite in their rims (Fig 4j).

  rim zoning visible in BSE images (Fig 6n) and a sector zoning visible even by optical microscopy (Fig 6o). This pattern defines 3 main zones, where zone A is the core and zones C and D form sector-zoned rims. Zone C matches the crystallographic sector (110), and zone D sector (100) (Ubide et al. 2019). In a few cases, an additional zone matching the crystallographic sector (010) is also present (e.g. ESM 6). All zones are optically distinguished by different tones of green, attributed to different Ti contents (Fig 6g) (e.g. Strong 1969; Ferguson 1973; Nielsen 1979). The contacts between different zoning are sharp in all occurrences (Fig 6). Zone A, the core, is systematically enriched in Ca, Sn, Hf, Zr (Fig 7a, b, c, g) and depleted in Na, Al and Fe 3+ compared to the other zones (

  7g). REE and Sc contents in cores of type-I aegirine (zone A) are higher than in the rims (zones B, C and D) for all complexes(Fig 6, 9 and Table3; similar maps of aegirine crystals from each complex are provided in Online Resource 4 to 8). Inside the rims, sector zone C is systematically enriched in REE compared to sector D (Fig6k, l, m). Nb and Ta do not have a constant behavior: their contents can be higher in sector zone C, be present only in fractures, or be enriched in an independent zonation, with all three cases being possible within the same sample. Similarly to amphiboles, the chondrite-normalized REE patterns of aegirine show high values for HREE relative to LREE, with lower medium REE. This particular pattern is known to be specific to Na-and Fe-enriched aegirine[START_REF] Mahood | Trace-element partitioning in pantellerites and trachytes[END_REF][START_REF] Shearer | Sector-zoned aegirine from the Ilimaussaq alkaline intrusion, South Greenland: Implications for trace-element behavior in pyroxene[END_REF]. The absolute REE concentrations, are quite similar from one complex to another, although small variations do exist(Fig 9).Trace element concentrations in type-II aegirine are highly dependent on the zoning pattern. This is notably true for Sn, Hf and Zr. Compared to type-I aegirine, Sn covaries with Ti instead of Zr. The shape of type-II aegirine REE patterns is similar to that of type-I, i.e. high HREE relative to LREE, with a depression in medium REE and a negative Eu anomaly. There is no common pattern in the REE contents of type-II aegirine: REE contents are higher in type-II than in type-I aegirine at Khan Bogd, generally lower at Evisa, and have similar values in Madagascar and at Amis (Fig7g).

  their major and trace chemistry: fluoro-arfvedsonite (Na 3 Fe 5 Si 8 O 22 (F,OH) 2 ), ferri-fluoro-leakeite (Na 3 Fe 4 LiSi 8 O 22 F 2 ), and ferro-ferri-katophorite (Na 2 Ca(Fe 2+ ,Mg) 4 Fe 3+ (Si 7 Al)O 22 (OH) 2 ). All of them occur in granites and pegmatites as euhedral crystals and contain inclusions of quartz and feldspar grains (Fig 4e), suggesting co-crystallization. Only amphiboles from Evisa have a variable composition, from fluoroarfvedsonite or ferro-ferri-katophorite in cores, respectively to ferri-fluoro-leakeite or fluoro-arfvedsonite in the rims. In the other complexes, only one amphibole type is present in granites and in pegmatites. Based on textural observations, we believe these changes are the sign of magmatic evolution. Indeed, the concentration of Li and

  Fig 6) suggests a sudden change in the crystallizing environment. Based on the scarce presence of type-I aegirine where amphiboles are present, we infer that they grew in competition. Only the rims C and D contain fine oscillatory zoning (e.g. TiO 2 on Fig 6g), fluid inclusions (Fig 4g), and, at Ambohimirahavavy, mineral inclusions that match the composition of those measured in hydrothermal pseudomorphs (e.g. zircon, bastnäsite-(Ce), Fig 4b and h). Therefore, we infer that type-I aegirine has a magmatic core (zone A on Fig 6p) while the rims are hydrothermal (zones C and D in Fig 6p). Zone B being only present at Khan Bogd, it will be discussed lower in this section.Chemical zoning in aegirine of our rocks matches that observed by[START_REF] Piilonen | The crystal chemistry of aegirine from Mont Saint-Hilaire, Quebec[END_REF] in aegirine from Mont 20 Saint-Hilaire, which they describe as a magmatic pyroxene overgrown by hydrothermal fibrous aegirine. More details about each element in the different zones is provided in the following paragraphs.The variations in Na, Ca, Al, and Fe can be explained by a common process. Type-I aegirine magmatic cores are enriched in Ca, and its hydrothermal rims are richer in Na, Al and Fe 3+ . Evidence such as widespread feldspar albitization, fluid inclusions composition, and the hydrothermal growth of aegirine without amphibole indicate that circulating hydrothermal fluids are oxidizing (above the magnetite-hematite buffer), of relatively high pH and rich in NaCl among other phases[START_REF] Salvi | The role of hydrothermal processes in concentrating high-field strength elements in the Strange Lake peralkaline complex, northeastern Canada[END_REF][START_REF] Smith | Metasomatic silicate chemistry at the Bayan Obo Fe-REE-Nb deposit, Inner Mongolia, China: contrasting chemistry and evolution of fenitising and mineralising fluids[END_REF][START_REF] Estrade | Le complexe cénozoïque alcalin d'Ambohimirahavavy à Madagascar : origine, évolution et minéralisations en métaux rares[END_REF][START_REF] Gysi | Lithogeochemical vectors for hydrothermal processes in the Strange Lake peralkaline granitic REE-Zr-Nb deposit[END_REF]. Depletion in Ca and enrichment in Na and Fe while moving on from aegirine magmatic core to hydrothermal rims can then be explained by intake of Na from the hydrothermal fluids and simultaneous crystallization of Ca-rich secondary minerals in the pseudomorphs (fluorite, bastnäsite-(Ce)

  addition, Sn is known to have a high solubility in Cl-rich fluids (e.g.[START_REF] Keppler | Partitioning of Cu, Sn, Mo, W, U, and Th between melt and aqueous fluid in the systems haplogranite-H 2 O-HCl and haplogranite-H 2 O-HF[END_REF], making it easily mobilized. This results in a much higher Sn content in the magmatic zone A than in the rims of aegirine. At Amis, amphibole in pegmatite as well as aegirine are both enriched in Sn, Zn and Pb in comparison with the other complexes of this study (Fig 7g). Ambohimirahavavy and Manongarivo have similar trace elements concentration (Fig 7g) and formed in the same geological context, a few km away from each other. Both amphibole and pyroxene from Khan Bogd have low trace elements concentrations. All these observations tend to indicate that the variations of absolute concentration in REE, Zn, Sn, Zr, U and Pb between aegirine crystals of the 6 complexes (Fig 7g), similarly to amphiboles, likely reflect parent melt compositions.

  hydrated Ca-and REE-rich elpidite (Fig 4i), and zone B grew simultaneously. The uptake of Ca and Zr by secondary elpidite from the fluid led to the lower concentration of these elements observed in zone B. The final stage, similarly to other complexes, was a more or less extensive pseudomorphism of elpidite and the formation of zones C and D on aegirine.Oscillating concentrations in zones C and D are most likely produced by autocatalytic surface attachment, or boundary layer effect[START_REF] Ortoleva | Geochemical self-organization I; reaction-transport feedbacks and modeling approach[END_REF] London 2008), which occurs when the crystal growth rate exceeds the diffusion rate of components in the fluid. Consider an initial solution saturated in component A, whose adsorption on the crystal surface is fast. As the adsorption progresses, the interface between the crystal and the solution becomes depleted in component A, and switches to saturation in component B. Component B is then adsorbed at the surface of the crystal, until the interface crystal-solution becomes saturated in component A again. This process results in a rhythmic but irregular zoning of aegirine crystals in zones C and D, consistent with a growth in a rapidly changing environment, which is a hydrothermal fluid.

  . The difference in REE concentration between aegirine zone A and zones B, C and D is due to the competition between aegirine and zirconosilicates for these elements, during the magmatic stage for zone A, and during the hydrothermal stage for the zones B, C and D. Scandium is sometimes grouped with REE, hence, it is not surprising to observe that, similarly to REE, it is enriched in the core of type-I aegirine.During the hydrothermal stage, aegirine is in competition with an altered, REE-enriched elpidite for zone B, and with REE-bearing minerals in pseudomorphs (e.g. bastnäsite-(Ce)) for zones C -(110) and D -(100). We therefore suggest that it is a competition for REE between aegirine and coprecipitating minerals having each a different REE affinity which explains the variations in REE concentrations among the aegirine core-to-rim zones. REE are also affected by sector zoning and are more concentrated in zone C than D, along with Ti and Ca.[START_REF] Beard | Clinopyroxene/melt trace element partitioning in sodic alkaline magmas[END_REF] studied experimentally the fractionation of REE in clinopyroxene in peralkaline melts and showed that the more Na and less Ca in M2 site of pyroxene, the smaller the M1 site, and the more HREE are incorporated. In addition, the authors mention that physico-chemical parameters not recorded in the composition of clinopyroxene, such as pressure, also play a significant role in the size of the M1

  the percentage difference between core and rim, C and R the mean values obtained for core and rim, respectively. Error bars shown in Fig 11 were obtained using common error calculation formulae, namely (Eq. 6) for the mean and (Eq. 7) for the percentage. error, n the number of measurements, d err the error on the percentage variation from core to rim, err C the mean error on core data, and err R the mean error on rim data. Resulting values indicate an enrichment in that element in the rim compared to the core when positive, while negative ones indicate a depletion (Fig 11).

  Fig 12b). As long as fluid flow continued, aegirine rims kept growing, although REE incorporation and LREE-HREE fractionation kept diminishing. By the time the hydrothermal stage waned, the two aegirine types were in equilibrium with the pseudomorphs and the partially replaced amphiboles (Fig 12c).

Figures caption Fig. 1 AFig. 2 Fig. 3 Fig. 4 aFig. 5 Fig. 6 Fig. 7 Fig. 8 Fig. 9

 caption123456789 Figures caption

  Table3). In pegmatites, we measured up to 2 wt% ZrO 2 in aegirine cores; although uncommon, Zr-rich aegirine was also described in other localities such as in nepheline syenites from the Motzfeldt Centre, South Greenland[START_REF] Jones | Zirconium-bearing aegirines from Motzfeldt, south Greenland[END_REF]; up to 7 wt% ZrO 2 ) and in metaluminous trachytes from the Warrumbungle Volcano, Australia[START_REF] Duggan | Zirconium-rich sodic pyroxenes in felsic volcanics from the Warrumbungle Volcano, Central New South Wales, Australia[END_REF] up to 14.5 wt% ZrO 2 ). Sector zone C is particularly enriched in Ti and Ca (Fig 6e, g), while D is in Fe (Fig 6c). Two major growth layers are visible in the rims with, in sector C, Ti and Ca decreasing towards the rims of each layer (Fig 6e, g). Al content is low in aegirine from all sixcomplexes (about 0.3 % Al 2 O 3 compared to standard average values of 1.2 % reported in aegirine; e.g.Deer et al. 1997a). At Khan Bogd, an additional zone, referred to as zone B, is found between the core and the sectorzoned rims of the crystals (Fig 4k); it shows random oscillations due to variations in many elements, including Ti

and Mn. The same zoning pattern, including zone B, was also observed by

[START_REF] Ranløv | Compositional zoning in hydrothermal aegirine from fenites in the Proterozoic Gardar Province, South Greenland[END_REF]

in nepheline syenite from the Narssaq Peninsula. However, they only documented variations in Zr, Ti, Al, Na, Ca and Fe.

[START_REF] Larsen | Clinopyroxenes and coexisting mafic minerals from the alkaline Ilimaussaq intrusion, South Greenland[END_REF]

and

[START_REF] Piilonen | The crystal chemistry of aegirine from Mont Saint-Hilaire, Quebec[END_REF]

respectively described a core-to-rim zoning in aegirine from nepheline syenite in the Ilímaussaq and Mont Saint-Hilaire complexes. They report cores enriched in Ca, Fe 2+ ,

  is richer than ferro-ferri-fluoro-leakeite from other complexes in Ti (0.18 vs 0.07 apfu), Sn (360 vs 32 ppm), Zn (0.24 vs 0.05 apfu), Pb (97 vs 7 ppm), LREE (168 vs 37 ppm on average) and HREE (140 vs 83 ppm on average) (Table1, Fig 5d). Mn and Zn vary together on a 1:1 relation at Strange Lake, as mentioned in the paper by

. Amphibole chondrite-normalized REE patterns are systematically enriched in HREE relative to LREE, with a depression in medium REE (MREE, from Eu to Dy). This pattern is unusual compared to more Ca-rich amphibole compositions in which REE patterns are usually enriched in LREE or in MREE (e.g.

[START_REF] Marks | Trace element variations in clinopyroxene and amphibole from alkaline to peralkaline syenites and granites: implications for mineral-melt trace-element partitioning[END_REF][START_REF] Coint | Use of trace element abundances in augite and hornblende to determine the size, connectivity, timing, and evolution of magma batches in a tilted batholith[END_REF]

. All spectra display a negative Eu anomaly suggesting early plagioclase fractionation (Fig 8). The absolute REE concentrations, though somewhat varying from one complex to another, are quite similar (Fig 8). Despite their similar patterns, the absolute concentrations of REE and other HFSE can vary in amphiboles from one complex to the other. Fluoro-arfvedsonite from Manongarivo and Khan Bogd both have low quantities of Sn (4 and 10 ppm) and HREE (35 and 50 ppm on average). Fluoro-arfvedsonite from Manongarivo also has less LREE than fluoro-arfvedsonite from other complexes (11 vs 60 ppm on average). Ferro-ferri-katophorite from Evisa displays high enrichments in LREE and HREE (respectively 224 and 290 ppm) (Fig 8). Ferro-ferri-fluoro-leakeite from Amis is notably different from other complexes in terms of trace elements: it

Table 2

 2 Mean values (number of analysis given in bold) for the major and trace element composition, obtained by microprobe and LA-ICPMS respectively, of aegirine from pegmatites and granites from six complexes worldwide in both pegmatites and granites. Granites only comprise type-II aegirine. T I, T II: aegirine type-I, II.

	Abbreviations: no.: number of replicates; b.d.l: below detection limit; Ambo: Ambohimirahavavy; KB: Khan
	bogd; M: Manongarivo; SL: Strange Lake

Table 3

 3 Mean composition of cores and rims of type-I aegirine from 5 complexes. Abbreviations: no.: number of

replicates; b.d.l: below detection limit; Ambo: Ambohimirahavavy; KB: Khan bogd; M: Manongarivo

This work was supported by an INSU/TelluS grant from CNRS (French National Center for Scientific Research). We thank Anthony Williams-Jones, Olga Vasyukova and Sam Broom-Fendley for providing some of the rock samples from the Strange Lake and Amis complexes and David Chew who helped with LA-ICP-MS mapping. Kathryn Goodenough provided helpful discussion. We also wish to thank CPM reviewers Charles