The Saint-Honor é Carbonatite REE Zone , Qu é bec , Canada : Combined Magmatic and Hydrothermal Processes

The Saint-Honoré carbonatite complex hosts a rare earth element (REE) deposit traditionally interpreted as being produced by late-stage hydrothermal fluids that leached REE from apatite or dolomite found in the early units and concentrated the REE in the late-stage units. New evidence from deeper units suggest that the Fe-carbonatite was mineralized by a combination of both magmatic and hydrothermal crystallization of rare earth minerals. The upper Fe-carbonatite has characteristics typical of hydrothermal mineralization—polycrystalline clusters hosting bastnäsite-(Ce), which crystallized radially from carbonate or barite crystals, as well as the presence of halite and silicification within strongly brecciated units. However, bastnäsite-(Ce) inclusions in primary magmatic barite crystals have also been identified deeper in the Fe-carbonatite (below 1000 m), suggesting that primary crystallization of rare earth minerals occurred prior to hydrothermal leaching. Based on the intensity of hydrothermal brecciation, Cl depletion at depth and greater abundance of secondary fluid inclusions in carbonates in the upper levels, it is interpreted that hydrothermal activity was weaker in this deepest portion, thereby preserving the original magmatic textures. This early magmatic crystallization of rare earth minerals could be a significant factor in generating high-volume REE deposits. Crystallization of primary barite could be an important guide for REE exploration.


Introduction
China has practically monopolized the production of rare earth elements (REE) over the last decades [1,2].In 2009, the Chinese government established a restrictive quota on REE exportation and exploration for new sources of REE rapidly increased around the world.Future demand is expected to be strong as the market for green technologies, the main consumer of REE, is projected to increase [3][4][5][6].Several geological processes or environments can produce economic REE concentrations [1,[7][8][9].However, carbonatites are REE-hosting rocks in which extractive metallurgy is better understood.According to Chakhmouradian and Zaitsev [9], 20 of the 50 most advanced REE extraction projects outside of China are found within igneous carbonate rocks, whereas most present-day or historical production has been from carbonatites or their weathered components.
Although carbonatites are known to have high natural concentrations of REE [10][11][12], in most cases these amounts are not sufficient enough to be economical [6,9,13].Specific geological processes are required to enhance REE concentrations.Four major processes are commonly invoked.First, the precipitation of REE-bearing minerals, such as minerals of the bastnäsite group [(REE)(CO 3 )F] or the parisite group [(REE) 2 Ca(CO 3 ) 3 F 2 ], may occur from primary carbonatite magma forming idiomorphic crystals of REE minerals in the late-stage carbonatite [14 -16].Second, a carbonate melt phase can be separated from a silicate melt in order to generate a REE-rich carbonate melt and a REE-poor silicate phase [17][18][19][20][21][22][23].Third, REE are remobilized and concentrated by hydrothermal solutions.A late hydrothermal solution can leach magmatic carbonates, apatite and sulphides to enrich late stage hydrothermal minerals in Ba, F, SO 4 2− , Sr, REE and Th [24][25][26].According to Mariano [15], in most carbonatites, "hydrothermal REE minerals occur in veinlets or as interstitial fillings, and appear as fine-grained polycrystalline clusters, commonly associated with barite, fluorite, hematite, quartz, strontianite and sulphides".This third process accounts for most REE concentration in carbonatites around the world.The fourth process is the result of chemical weathering of sub-surface carbonatite.As the carbonate minerals (calcite and dolomite) are easily dissolved, Ca and Mg are removed, thereby concentrating the less mobile elements, such as REE [15,27,28].However, despite the great abundance of REE concentrated by hydrothermal processes in carbonatite, these concentrations rarely form high-tonnage economically viable deposits [15].Early magmatic REE concentrations may therefore represent an important geological process favouring the formation of world-class REE deposits in carbonatites.The Saint-Honoré carbonatite complex, exploited by Niobec Inc., hosts a niobium deposit as well as a REE zone in its central portion.Presently, only the niobium deposit is exploited.Previous studies of the Saint-Honoré REE zone [29,30] highlighted the role of a hydrothermal fluid in concentrating REE through fine-grained bastnäsite-(Ce) mineralization.In 2011-2012, this REE zone was drilled to a depth of 1250 m, providing an opportunity to characterize the deeper portion of the REE zone.The Saint-Honoré carbonatite is classified as a world-class deposit with its total proven and inferred resources estimated at 1058.6 Mt at 1.75% total rare earth elements (TREE) [31].Based on its high tonnage, the REE mineralization of the Saint-Honoré carbonatite is not the usual, expected hydrothermal REE mineralization.The main objective of this study is to demonstrate the importance of the igneous preconcentration of the REE prior to hydrothermal enrichment to create a large deposit.

Geological Setting
The Saint-Honoré carbonatite complex is located in the Saguenay region of Québec, Canada (Figure 1).This region is comprised of three petrographic domains [32,33]: (1) a gneiss complex (orthogneiss and paragneiss) deformed during the Hudsonian Orogeny (1735 Ma); (2) anorthosite and mangerite-charnockite batholiths dated to pre-to post-Grenvillian orogeny (935 Ma); and (3) younger calc-alkaline intrusions (which include the Saint-Honoré alkaline complex) related to the St. Lawrence River Rift and the associated Saguenay Graben [34].This latter extensional event provided crustal anisotropies that channeled alkaline magmatism, including the Saint-Honoré carbonatite [35].An age of 565 Ma was obtained by Doig and Barton [36] for the alkaline complex.McCauland et al. [37] report an age of 571 ± 5 Ma from phlogopite in the carbonatite while Kamenetsky et al. [38] report a Rb-Sr model age of 564 ± 8 Ma based on two points.A U-Pb age of 582.1 ± 1.8 Ma was obtained on baddeleyite from the carbonatite (Michael Higgins, UQAC, personal communication 2015).The Saint-Honoré alkaline complex is a vertical intrusion composed of a carbonatite core surrounded by alkaline silicate rocks.Fenitization is present in host rocks and is characterized by the metasomatism of pyroxenes into sodic amphiboles and carbonates, sericitization of plagioclase, and formation of many carbonatitic red to green veins [40].The alkaline silicate units are mainly composed of nepheline syenite and younger alkaline syenite (all miaskitic, (Na2O + K2O)/Al2O3 < 1: molar proportion) [40].Some xenoliths of alkaline syenite are observed in the carbonatitic part, suggesting that the carbonatite is younger than the syenites [40].The carbonatitic units have an elliptical shape oriented north-east south-west (Figure 2) where the Fe-carbonatite is in the centre, enclosed by the Mg-carbonatite which is itself enclosed by the outer Ca-carbonatite.The main carbonatite crystallized at high temperatures (1000-1150 °C, using the apatite-phlogopite geothermometer [30]), clearly suggesting a magmatic origin for the Mg-carbonatite.However, Fournier [30] obtained temperatures of 364 to 321 °C for the weathered upper zone of the Fe-carbonatite (chlorite geothermometer).Crosscutting relationships show an evolution trend from Ca-carbonatite to Mg-carbonatite to Fe-carbonatite [40].This evolution is in agreement with geochemical data: MnO, which reflects the evolution of a carbonatite complex [12], increases from the calcitic rocks to the Fe-carbonatite (Table S1).Table 1 (from Fortin-Belanger [40]) shows partial geochemical analyses of the different units.Although these older results do not have the same quality as modern ones and many of these units are not available anymore, results are in agreement with modern ones and they are fit for the purpose [41].La, Ce, and metal concentrations increase from older to younger units.This pattern is interpreted as the result of magmatic fractionation of these less compatible elements during the crystallization of carbonates [30].An enrichment in REE can also be caused by the formation of immiscible carbonatitic liquids that capture most of trace and ultra-trace elements compatible with the carbonatite melt [17][18][19][20][21][22][23].Furthermore, the REE deposit is hosted in the The Saint-Honoré alkaline complex is a vertical intrusion composed of a carbonatite core surrounded by alkaline silicate rocks.Fenitization is present in host rocks and is characterized by the metasomatism of pyroxenes into sodic amphiboles and carbonates, sericitization of plagioclase, and formation of many carbonatitic red to green veins [40].The alkaline silicate units are mainly composed of nepheline syenite and younger alkaline syenite (all miaskitic, (Na 2 O + K 2 O)/Al 2 O 3 < 1: molar proportion) [40].Some xenoliths of alkaline syenite are observed in the carbonatitic part, suggesting that the carbonatite is younger than the syenites [40].The carbonatitic units have an elliptical shape oriented north-east south-west (Figure 2) where the Fe-carbonatite is in the centre, enclosed by the Mg-carbonatite which is itself enclosed by the outer Ca-carbonatite.The main carbonatite crystallized at high temperatures (1000-1150 • C, using the apatite-phlogopite geothermometer [30]), clearly suggesting a magmatic origin for the Mg-carbonatite.However, Fournier [30] obtained temperatures of 364 to 321 • C for the weathered upper zone of the Fe-carbonatite (chlorite geothermometer).Crosscutting relationships show an evolution trend from Ca-carbonatite to Mg-carbonatite to Fe-carbonatite [40].This evolution is in agreement with geochemical data: MnO, which reflects the evolution of a carbonatite complex [12], increases from the calcitic rocks to the Fe-carbonatite (Table S1).Table 1 (from Fortin-Belanger [40]) shows partial geochemical analyses of the different units.Although these older results do not have the same quality as modern ones and many of these units are not available anymore, results are in agreement with modern ones and they are fit for the purpose [41].La, Ce, and metal concentrations increase from older to younger units.This pattern is interpreted as the result of magmatic fractionation of these less compatible elements during the crystallization of carbonates [30].An enrichment in REE can also be caused by the formation of immiscible carbonatitic liquids that capture most of trace and ultra-trace elements compatible with the carbonatite melt [17][18][19][20][21][22][23].Furthermore, the REE deposit is hosted in the youngest and most evolved unit: the Fe-carbonatite.Yttrium, which behaves similarly to heavy rare earth elements (HREE), decreases with fractionation (Table 1), suggesting an early precipitation of HREE [40].Only the light rare earth elements (LREE) are interpreted to have been concentrated by fractional crystallization to form the REE zone [29,30,40].This fractionation between HREE and LREE is common in carbonatite and can be explained by the higher compatibility of HREE in early minerals such as apatite [42].Compared to the carbonatitic units, the alkaline silicate units are significantly poorer in REE but richer in Zr (Table 1; [40]).youngest and most evolved unit: the Fe-carbonatite.Yttrium, which behaves similarly to heavy rare earth elements (HREE), decreases with fractionation (Table 1), suggesting an early precipitation of HREE [40].Only the light rare earth elements (LREE) are interpreted to have been concentrated by fractional crystallization to form the REE zone [29,30,40].This fractionation between HREE and LREE is common in carbonatite and can be explained by the higher compatibility of HREE in early minerals such as apatite [42].Compared to the carbonatitic units, the alkaline silicate units are significantly poorer in REE but richer in Zr (Table 1; [40]).Table 1.Representative analyses of different units within the Saint-Honoré complex (from [40]).The alkaline complex is covered by Ordovician limestone and the uppermost 60-120 m of the Fe-carbonatite were heavily weathered to an orange or red coloured carbonatite-possibly caused by the presence of free hematite-prior to deposition of the limestone [29,30].
According to previous studies [29], the principal economic REE minerals of the Saint-Honoré deposit are late-stage fluoro-carbonates that occur as needles of 0.5-2 μm in diameter and up to 20 μm in length.Although previous studies [29,30] suggested that mineralization was related to a late-stage orthomagmatic hydrothermal system, they were limited to the first few hundreds of metres below ground level of weathered Fe-carbonatite.Table 1.Representative analyses of different units within the Saint-Honoré complex (from [40]).The alkaline complex is covered by Ordovician limestone and the uppermost 60-120 m of the Fe-carbonatite were heavily weathered to an orange or red coloured carbonatite-possibly caused by the presence of free hematite-prior to deposition of the limestone [29,30].
According to previous studies [29], the principal economic REE minerals of the Saint-Honoré deposit are late-stage fluoro-carbonates that occur as needles of 0.5-2 µm in diameter and up to 20 µm in length.Although previous studies [29,30] suggested that mineralization was related to a late-stage orthomagmatic hydrothermal system, they were limited to the first few hundreds of metres below ground level of weathered Fe-carbonatite.

Petrographic Studies
The macroscopic observations were performed by logging 19 deep drill cores located almost entirely in the REE zone.Microscopic observations were made from 100 polished thin sections taken from these drill cores, selected to cover representative vertical and lateral distributions.The thin sections were prepared by SGS Canada Inc. (Mississauga, ON, Canada) with oil-based lubricants and ethylene glycol to preserve the highly soluble minerals, such as halite, present in these rocks.

Scanning Electron Microscope (SEM)
Imaging was obtained through a Zeiss EVOMA15 HD scanning electron microscope (SEM) coupled with an energy dispersive spectrometer-silicon drift detector (EDS-SDD) from Oxford Instruments (Model X-Max 150, Wycombe, UK) from IOS Services Géoscientifiques Inc., Chicoutimi, QC, Canada.Back-scattered electron images, chemical mapping and phase identification were obtained with the following parameters: 40 µA, 20 kV with a working distance of 11.5 µm.

Whole Rock Analysis
Whole rock analyses presented in this study are from the deepest drill hole (located in Figure 2) as it is the most representative and there is a strong lateral continuity between the different drill holes.Whole rock analyses were generated in two main batches.The first batch was lower quality routine analysis and are not presented in this study.The deepest complete vertical drill hole has been reanalyzed by SGS Canada Inc. with the highest quality available to offer a complete vertical section.Major elements were determined by XRF and trace elements by ICP-MS following a multi-acid dissolution (quality controlled by SGS Canada with ORES-101 and GRE-2 reference materials).Chlorine was determined by titration.Finally, C and S were analyzed at LabMaTer (Université du Québec à Chicoutimi) using an elemental analyzer (Horiba, EMIA 220-V, Longjumeau, France) following the method of Bédard et al. [43] and WMS-1a and NBS-88a as reference materials [43,44].All results were within the acceptable limits of uncertainty for our purposes [41].

Petrography of the Fe-Carbonatite
The upper 80 m were previously described [29,30] as a sub-lateritic carbonatite, rich in chlorite and hematite.This heavy alteration is interpreted to be the result of past tropical climate and/or sea infiltration, as suggested by the presence of the Ordovician limestone on top of the carbonatite (Figure 3).As the main objective of this paper is to present the magmatic and hydrothermal components associated with the REE mineral crystallization, the weathered uppermost section of the deposit is not considered further except for comparative purposes.infiltration, as suggested by the presence of the Ordovician limestone on top of the carbonatite (Figure 3).As the main objective of this paper is to present the magmatic and hydrothermal components associated with the REE mineral crystallization, the weathered uppermost section of the deposit is not considered further except for comparative purposes.The Fe-carbonatite is divided into two units: REE-poor Fe-carbonatite (PFeC) and REE-rich Fe-carbonatite (RFeC).Spatially, the RFeC is surrounded by the PFeC, and centimetre-to metre-scale crosscutting injections of the RFeC into the PFeC indicate that the RFeC represents a later stage than the PFeC.Considering that the REE-poor Fe-carbonatite does not hold important mineralization or critical information, it is not further discussed.

REE-Rich Fe-Carbonatite (RFeC)
This unit is the main host of REE minerals and consequently the focus of this study.At a macroscopic scale, the RFeC is primarily composed of carbonates (dolomite, dolomite with ankerite rim, ankerite, with minor amounts of calcite and siderite), barite, and purple polycrystalline clusters hosting REE mineralization (Figure 4D).Halite, quartz, pyrite, apatite, Nb-rutile, and thorianite are the minor phases with varying abundance from cluster to cluster.They do not show alteration and are interpreted to be hydrothermal.Two end-member textures are present: massive coarse-grained carbonatite (Figure 4D) and fine-grained hydraulic breccia (Figure 4A).Some of the contacts between the massive and brecciated carbonatite textures are sharp whereas others are gradual on the scale of a few centimetres or metres.The Fe-carbonatite is divided into two units: REE-poor Fe-carbonatite (PFeC) and REE-rich Fe-carbonatite (RFeC).Spatially, the RFeC is surrounded by the PFeC, and centimetre-to metre-scale crosscutting injections of the RFeC into the PFeC indicate that the RFeC represents a later stage than the PFeC.Considering that the REE-poor Fe-carbonatite does not hold important mineralization or critical information, it is not further discussed.

REE-Rich Fe-Carbonatite (RFeC)
This unit is the main host of REE minerals and consequently the focus of this study.At a macroscopic scale, the RFeC is primarily composed of carbonates (dolomite, dolomite with ankerite rim, ankerite, with minor amounts of calcite and siderite), barite, and purple polycrystalline clusters hosting REE mineralization (Figure 4D).Halite, quartz, pyrite, apatite, Nb-rutile, and thorianite are the minor phases with varying abundance from cluster to cluster.They do not show alteration and are interpreted to be hydrothermal.Two end-member textures are present: massive coarse-grained carbonatite (Figure 4D) and fine-grained hydraulic breccia (Figure 4A).Some of the contacts between the massive and brecciated carbonatite textures are sharp whereas others are gradual on the scale of a few centimetres or metres.The coarse-grained massive texture is composed of millimetre-to centimetre-sized carbonates, purple polycrystalline clusters hosting the REE mineral (bastnäsite-(Ce) and monazite-(Ce)) and barite in varying proportions.The carbonates, mainly dolomite and ankerite, are primary minerals, but exhibiting in certain places, late and locally low to medium recrystallization and/or brecciation.Some calcite crystals have an alteration rim of Fe-rich carbonate, which is siderite (Figure 5A).
The barite crystals are, in general, finely recrystallized, but minor amounts of crystals exhibit a relict core of primary barite and have the shape of euhedral primary barite (Figure 4C,F).Moreover, The coarse-grained massive texture is composed of millimetre-to centimetre-sized carbonates, purple polycrystalline clusters hosting the REE mineral (bastnäsite-(Ce) and monazite-(Ce)) and barite in varying proportions.The carbonates, mainly dolomite and ankerite, are primary minerals, but exhibiting in certain places, late and locally low to medium recrystallization and/or brecciation.Some calcite crystals have an alteration rim of Fe-rich carbonate, which is siderite (Figure 5A).such large barite crystals show homogeneous extinction suggesting they were a single euhedral large crystal (Figure 4C).Other barite crystals exhibit dissolution patterns associated with quartz and carbonate (Figure 5B,D) and halite (Figure 5C).The solubility of barite is known to increase in the presence of Cl -in fluids [55].The barite crystals are, in general, finely recrystallized, but minor amounts of crystals exhibit a relict core of primary barite and have the shape of euhedral primary barite (Figure 4C,F).Moreover, such large barite crystals show homogeneous extinction suggesting they were a single euhedral large crystal (Figure 4C).Other barite crystals exhibit dissolution patterns associated with quartz and carbonate (Figure 5B,D) and halite (Figure 5C).The solubility of barite is known to increase in the presence of Cl -in fluids [55].
The following minerals have been observed in clusters: REE fluoro-carbonates (bastnäsite-(Ce), parisite, synchysite), monazite-(Ce), hematite, quartz, halite, barite, dolomite, and ankerite.However, not all of these minerals are present in all clusters.The clusters are generally ameboidal in shape and interstitial to the carbonates (Figure 4D,B).In these clusters, electron microscopy reveals that the barite exhibits a resorbed texture (Figure 5D), associated with the presence of halite (Cl − ) and/or silicification.Silicification is common in the hydrothermal REE mineralization of carbonatites [56].Hexagonal quartz crystals are locally present with barite suggesting it may be synchronous with silicification.The REE fluoro-carbonates are present primarily as bastnäsite-(Ce) in the form of fine needles (5-20 µm in length; Figure 4A,B).The fluoro-carbonates are usually crystallized radially from the edge of a carbonate or a resorbed barite (Figure 4B).Some overgrowth of bastnäsite-(Ce) has been observed from the core (parisite-(Ce)) to the rim (bastnäsite-(Ce)).The REE fluoro-carbonates are mainly present in the form of bastnäsite-(Ce) [(Ce,La)CO 3 F] and the composition varies from parisite [Ca(Ce,La) 2 (CO 3 ) 3 F 2 ] or synchisite [Ca(Ce,La)(CO 3 ) 2 F] to bastnäsite-(Ce).The monazite-(Ce) [(La,Ce,Nd)PO 4 ], also found in the form of needles, is even smaller (less than 1 µm in length) and is interstitial to other minerals.Usually, the monazite-(Ce) is associated with µm-sized hematite.
The size of clasts within the brecciated unit is variable suggesting a variation of brecciation intensity.The breccia is also mineralized with REE.The clasts are composed of carbonate and barite, implying that they were both present prior to the brecciation.In thin sections where the brecciation intensity is low, the clasts can usually be fit together into their initial shape similar to a puzzle (jigsaw texture), suggesting hydraulic brecciation [57] (Figure 4A).Bastnäsite-(Ce) and monazite-(Ce), the most common REE minerals in the breccia, usually crystallize radially from the clasts or from the edge of the breccia as an overprint.Furthermore, no early clasts of REE minerals are observed, suggesting that the mineralization occurred post-brecciation.The transportation distance of clasts appears to be relatively short, according to the sharp angles and sizes of these clasts (mm to cm).
An infiltration by surficial water is interpreted to have occurred in the upper 250 m of the carbonatite [30].In this interval, halite, a highly soluble mineral, is completely absent.Additionally, barite is significantly less abundant and almost absent from the polycrystalline clusters.This observation could be explained by the relative solubility of barite in a Cl-rich fluid [55], as suggested by the presence of halite.The clusters of this upper section exhibit a porosity of 3-10% while there is no observable porosity in the deeper zone.Furthermore, interstitial minerals of the chlorite group in this section are observed in dissolved carbonates.
In the deepest drilled portion of the Saint-Honoré Fe-carbonatite (more than 1000 m below the surface), there are several important textural and mineralogical variations.First, grain size is significantly larger than in the upper zone.Some crystals can reach up to 10 cm in length.Furthermore, in some places, barite, which is mostly euhedral to sub-euhedral, exhibits a zoned pattern featuring a pristine core and an altered rim composed of the cluster assemblage (Figure 4F).The observed breccia has only large and angular clasts, suggesting that only low-intensity brecciation is present in the deepest zone.Finally, barite from the deeper part of the deposit show bastnäsite-(Ce) and monazite-(Ce) inclusions implying early co-crystallization of barite and fluoro-carbonates (Figures 4E  and 5E,F; Table S2).The minerals of the apatite group in the RFeC occurs in cm-scale stringers that crosscut all other minerals and units.In the PFeC, the minerals of the apatite group are reddish and their textures suggest a replacement of the carbonates.

Geochemistry of the Fe-Carbonatite
For this section, only the deepest drill core has been considered as there is a strong lateral continuity.In a Ca-Mg-Fe classification diagram [58], most samples (Table S1) from Saint-Honoré fall within the Fe-carbonatite field (Figure 6).A linear trend exists between the Ca and Fe concentrations as the Mg concentration is relatively constant as shown in Figure 6.This type of trend is usually interpreted as being the result of a magmatic evolution [15], however this evolution trend does not show any relation with depth.There is no geochemical evolution trend with depth for the major elements as suggested by SiO 2 , Al 2 O 3 , and K 2 O + Na 2 O (Figure 7) except for the weathering in the first 250 m.The strongly-weathered, uppermost portion (upper 250 m) of the Fe-carbonatite is characterized by relatively higher concentrations of SiO 2 , Al 2 O 3 , and K 2 O + Na 2 O (Figure 7), corresponding to the presence of phlogopite [29,30] and chlorite [29].
The Cl concentrations (Table S1) vary from 0.02% to 0.43% and BaO from 0.03% to 13.96%.Overall, Cl and BaO vary significantly with depth (Figure 8): Cl decrease while BaO increase with depth except for the weathered upper 250 m (Figure 8).It reflects a decrease in halite and an increase in barite proportion with depth.This relationship is in agreement with petrographic observations suggesting barite resorption was induced by a Cl-rich fluid in the upper zone [55].The RFeC of Saint-Honoré has a BaO/SrO ratio ranging from 4 to 58 (>>1).Additionally, the uppermost 250 m portion is essentially depleted in Cl and BaO (Figure 8), reflecting the preferential weathering of more soluble minerals, such as halite and barite.The presence of halite below 250 m (Figure 8) suggests minimal interaction with surficial water at depth.
The abundance of Cl and BaO is anomalously high compared to the Woolley and Kempe [10] average Cl (0.02%) and BaO (0.88%) concentrations in Fe-carbonatite.In most Fe-carbonatite samples, no primary barite is observed because Ba is substituted in calcite, dolomite, and ankerite [12].Although Woolley and Kempe [10] also documented high-BaO Fe-carbonatite with 9.6 to 20.6% BaO (interpreted as being related to barite crystallization), they excluded these anomalous data from their calculation of "ordinary" Fe-carbonatite.The Cl concentrations (Table S1) vary from 0.02% to 0.43% and BaO from 0.03% to 13.96%.Overall, Cl and BaO vary significantly with depth (Figure 8): Cl decrease while BaO increase with depth except for the weathered upper 250 m (Figure 8).It reflects a decrease in halite and an increase in barite proportion with depth.This relationship is in agreement with petrographic observations suggesting barite resorption was induced by a Cl-rich fluid in the upper zone [55].The RFeC of Saint-Honoré has a BaO/SrO ratio ranging from 4 to 58 (>>1).Additionally, the uppermost 250 m portion is essentially depleted in Cl and BaO (Figure 8), reflecting the preferential weathering of more soluble minerals, such as halite and barite.The presence of halite below 250 m (Figure 8) suggests minimal interaction with surficial water at depth.
The abundance of Cl and BaO is anomalously high compared to the Woolley and Kempe [10] average Cl (0.02%) and BaO (0.88%) concentrations in Fe-carbonatite.In most Fe-carbonatite samples, no primary barite is observed because Ba is substituted in calcite, dolomite, and ankerite [12].Although Woolley and Kempe [10] also documented high-BaO Fe-carbonatite with 9.6 to 20.6% BaO (interpreted as being related to barite crystallization), they excluded these anomalous data from their calculation of "ordinary" Fe-carbonatite.The Saint-Honoré RFeC has strong REE concentrations where the sum of LREE are between 30,000 and 100,000 times the chondrite values and HREE are between 500 and 1000 chondrites (Table 2, Figure 9, Table S1).The REE concentrations do not follow the magmatic evolution trend.Furthermore, no spatial relationship between REE concentration and depth or lateral distance is observed within the  The Saint-Honoré RFeC has strong REE concentrations where the sum of LREE are between 30,000 and 100,000 times the chondrite values and HREE are between 500 and 1000 chondrites (Table 2, Figure 9, Table S1).The REE concentrations do not follow the magmatic evolution trend.Furthermore, no spatial relationship between REE concentration and depth or lateral distance is observed within the The Saint-Honoré RFeC has strong REE concentrations where the sum of LREE are between 30,000 and 100,000 times the chondrite values and HREE are between 500 and 1000 chondrites (Table 2, Figure 9, Table S1).The REE concentrations do not follow the magmatic evolution trend.Furthermore, no spatial relationship between REE concentration and depth or lateral distance is observed within the RFeC.The upper and deeper zone REE patterns have different La/Lu ratio: the upper zone (250-1000 m) has (La/Lu) cn ≈ 460 and the deeper zone (below 1000 m), has (La/Lu) cn ≈ 3100 (Table S1).
The LREE to HREE ratio is more important at depth.In general, the HREE concentrations of the Saint-Honoré Fe-carbonatite are relatively similar to non-mineralized carbonatites [10], but the LREE are one to two orders of magnitude higher (Figure 9).The REE patterns are strongly fractionated which is typical of REE deposits in carbonatites [8,9,15,59,60].In general, based on >15,000 REE analyses performed in the RFeC (IAMGold unreleased data), no significant higher-grade zones or vertical variations exist.The mineralization is broadly and homogeneously distributed throughout the Fe-carbonatite.RFeC.The upper and deeper zone REE patterns have different La/Lu ratio: the upper zone (250-1000 m) has (La/Lu)cn ≈ 460 and the deeper zone (below 1000 m), has (La/Lu)cn ≈ 3100 (Table S1).The LREE to HREE ratio is more important at depth.In general, the HREE concentrations of the Saint-Honoré Fe-carbonatite are relatively similar to non-mineralized carbonatites [10], but the LREE are one to two orders of magnitude higher (Figure 9).The REE patterns are strongly fractionated which is typical of REE deposits in carbonatites [8,9,15,59,60].In general, based on >15,000 REE analyses performed in the RFeC (IAMGold unreleased data), no significant higher-grade zones or vertical variations exist.The mineralization is broadly and homogeneously distributed throughout the Fe-carbonatite.

Discussion
The magmatic component of carbonatite-hosted REE deposits is, in general, less often described either as it is overprinted by hydrothermal activity or because access to the magmatic evidence is physically difficult.Via numerous drill cores down to 1200 m below the surface, the Saint-Honoré deposit provides an exceptional access to deeper portions of a REE-bearing Fe-carbonatite.Consequently, observations of the more primitive/pristine units are possible.
The proposed hypothesis that late-stage hydrothermal REE remobilization occurred in the RFeC of Saint-Honoré is supported by a number of observations: (1) hydraulic brecciation; (2) in the breccia, REE minerals are precipitated post-brecciation (for example, the minerals crystalized radially (Figure 4B) from clast edges); (3) the mineralization is associated with minerals typical of hydrothermal REE deposition in carbonatites, such as quartz, hematite, barite and late calcite (Figure 5D); (4) in the late-stage polycrystalline clusters, REE mineralization overprints primary barite and crystallization occurs radially from the barite edges (Figure 4B); (5) mineralization occurred in the form of tiny needles (Figures 4B and 5B,D), suggesting a rapid lower-temperature crystallization (Figures 4B and 5B,D); and (6) La/Yb ratio are similar to hydrothermal REE deposits such as Bear Lodge (Figure 10) [59,62].
Halite is commonly observed in the dissolution shadow of barite.Halite is also present as deep late-stage veins in Ca-and Mg-carbonatite mined for niobium [63].Most of the halite observed is late-stage (in polycrystalline clusters, associated with dissolution textures, vein and breccia filling), and interpreted to be related to the hydrothermal activity although halite could have a magmatic origin in the Saint-Honoré Mg-carbonatite as suggested by Kamenetsky et al. [38].However, subsequent work by Tremblay et al. [63] proposed that halite minute crystals within pyrochlore grains are hydrothermal and Na comes from pyrochlore alteration.Given that the deposit is overlain by Ordovician limestone, it also could be argued that the halite is of marine origin.However, neither sylvite nor gypsum, typical of seawater infiltration, are present.Furthermore, chlorine is rarely observed in a carbonatite-related fluid.The best example is Ol Doinyo Lengai, which features salt fumaroles [64].However, fenitization is a clear indication that Na-and/or K-rich fluids were associated with the carbonatites.Furthermore, the existence of halides within the upper part of the RFeC, is probably hidden or lost because of its extreme solubility (e.g.dissolved by meteoric water or during drilling, sampling or thin-section preparation), and most studies have access only to near-surface samples.The experimental data of Williams-Jones et al. [28] suggest that Cl − is an effective ligand for transporting REE in solution.Migdisov and Williams-Jones [65] demonstrated that Cl − as well as sulphate complexes are potent ligands for mobilizing REE as observed in Wicheeta carbonatite [66].
Different sources of REE in carbonatites have been proposed: (1) magmatic REE concentrations in the orthomagmatic fluid generated by fractional crystallization [30]; (2) remobilization of REE leached from primary magmatic minerals, such as carbonates or apatite [6,15,59,67,68] where REE are not an essential constituent (Table S2); and (3) remobilization of REE from early magmatic REE minerals, where REE are essential elements of the minerals from an initially REE-enriched melt [7][8][9].To identify the best source model for the Saint-Honoré REE deposit, the most pristine magmatic Fe-carbonatites were identified.
Barite in REE deposits is usually a late-stage mineral precipitated from an orthomagmatic hydrothermal solution [69].However, the barite of the RFeC presents many lines of evidence for early magmatic crystallization, such as: (1) pseudomorph of early barite by polycrystalline recrystallized barite (Figure 4C,F); (2) clasts of barite in the breccia; and (3) resorption textures in the barite of polycrystalline clusters (Figure 5D).The magmatic crystallization of barite from carbonatite magma is not common, but experimental data demonstrated that this is possible when Ba abundance is high [70,71].For example, Jones and Wyllie [72] made a glass out of the Mountain Pass carbonatite with over 6% Ba, thereby demonstrating that such a liquid (carbonate, barite, and La-bearing mineral) is clearly possible.The Mountain Pass deposit, having up to 25% BaO [73,74], is a prime example of a Ba-rich carbonatite.Usually, the BaO/SrO ratio from Ca-carbonatite to Fe-carbonatite increases from 0.4 to 1, reflecting an increasing BaO concentration [10,12].For the Ba-rich Fe-carbonatites of the Saint-Honoré carbonatites, the BaO/SrO ratio is very high (>> 1) because Sr has not been enriched to the same degree as Ba.Moreover, the Mountain Pass REE deposit is a well-known example of primary euhedral barite crystals in a carbonatite.In the RFeC, the vertical geochemical variation in Cl and Ba (Figure 8) appears to reflect the intensity of hydrothermal activity.As mentioned, barite is more soluble in a Cl-rich fluid [55].The presence of Cl can be related to the hydrothermal fluids or can be sourced from magmatic activity [38,63].Therefore, the upper part of the deposit, which is rich in Cl and poor in Ba, could have been exposed to more intense hydrothermal activity, resulting in extensive resorption of primary barite.Inversely, the lower part, which is richer in Ba and poorer in Cl, is interpreted as a better preserved portion that was exposed to lesser hydrothermal activity.Petrographic observations also support this lower hydrothermal activity in the lower portion of the RFeC: (1) the brecciation is significantly less intense; (2) the crystal sizes are larger; and (3) the barite crystals feature a pristine core and an altered rim (Figure 4C,F).Furthermore, the preliminary results from fluid inclusion gas analyses by solid probe mass spectrometry [75] revealed a progressive decrease in secondary fluid inclusions trapped in carbonates with depth [76].Moreover, the REE pattern of the deeper zone are typical of magmatic REE deposits such as Mountain Pass (Figure 10) [73,74].The upper part of the RFeC has a chondrite-normalized REE pattern similar to the Bear Lodge hydrothermal-influenced carbonatite ((La/Lu) cn ≈ 540 [59]; Figure 10) while the deeper zone of the RFeC has a chondrite-normalized REE pattern more similar to the magmatic Mountain Pass deposit ((La/Lu) cn ≈ 3900; Figure 10) [73,74].as Ba.Moreover, the Mountain Pass REE deposit is a well-known example of primary euhedral barite crystals in a carbonatite.In the RFeC, the vertical geochemical variation in Cl and Ba (Figure 8) appears to reflect the intensity of hydrothermal activity.As mentioned, barite is more soluble in a Cl-rich fluid [55].The presence of Cl can be related to the hydrothermal fluids or can be sourced from magmatic activity [38,63].Therefore, the upper part of the deposit, which is rich in Cl and poor in Ba, could have been exposed to more intense hydrothermal activity, resulting in extensive resorption of primary barite.Inversely, the lower part, which is richer in Ba and poorer in Cl, is interpreted as a better preserved portion that was exposed to lesser hydrothermal activity.Petrographic observations also support this lower hydrothermal activity in the lower portion of the RFeC:  As mentioned previously, bastnäsite-(Ce) crystals hosted in primary barite (Figures 4E and 5F; Table S2) are observed in the lower part of the RFeC.This occurrence is clearly related to inclusions and not exsolution as: (1) barite and bastnäsite-(Ce) have different crystallographic systems and (2) the orientation of bastnäsite-(Ce) crystals in barite are not parallel to the barite crystal cleavages (Figure 4E).This relationship clearly indicates an early crystallization of magmatic bastnäsite-(Ce).The co-crystallization of REE minerals and barite from a carbonatite magma has already been As mentioned previously, bastnäsite-(Ce) crystals hosted in primary barite (Figures 4E and 5F; Table S2) are observed in the lower part of the RFeC.This occurrence is clearly related to inclusions and not exsolution as: (1) barite and bastnäsite-(Ce) have different crystallographic systems and (2) the orientation of bastnäsite-(Ce) crystals in barite are not parallel to the barite crystal cleavages (Figure 4E).This relationship clearly indicates an early crystallization of magmatic bastnäsite-(Ce).The co-crystallization of REE minerals and barite from a carbonatite magma has already been demonstrated; Mountain Pass is probably the best natural example of primary bastnäsite-(Ce) crystals co-occurring with barite [73,74,78].Experimental studies [14, 71,72] have demonstrated the possibility of coeval magmatic barite and REE mineral crystallization with relatively high concentrations of Ba and REE in magma.The principal difference between the RFeC and the PFeC of Saint-Honoré is the presence of both barite and REE minerals.More experimental studies need to be performed to understand the role of barite and the related formation conditions.Despite the dominant interpretation of hydrothermal REE mineralization in carbonatites, early magmatic crystallization of REE minerals could be an important process in generating world-class carbonatite-hosted REE deposits.Song et al. [79] show that REE, Ba, Sr, Mo, and W will partition strongly in the carbonatitic within magma rather than within fluid, hence suggesting that a magmatic stage is probably more important than previously recognized.A magmatic stage could even be the precursor to many hydrothermal REE deposits, given that their deeper portions are often not accessible for documenting such processes.
From a petrogenetic perspective, Girnis et al. [80] studied the experimental partitioning of selected trace elements between mantle minerals and silico-carbonate melts at 6-12 GPa and 1300-1700 • C.They show that high La/Lu and BaO/SrO ratios can be produced by a low primary mantle melting ratio as LREE and BaO are more incompatible than HREE and SrO in mantle minerals.Such processes could have led to the primary carbonatite magmas that afterward differentiated to produce the LREE and BaO rich Fe-carbonatite as it could have been the case in the Saint-Honoré deposit.
The observations presented in this study significantly improve the ore genesis model of the Saint-Honoré REE deposit.First, the carbonatite complex evolves from Ca-carbonatite to Fe-carbonatite by crystal fractionation, concentrating incompatible elements, such as REE, in the late-stage of Fe-carbonatite formation.During the crystallization of the Fe-carbonatite, the concentrations of both REE and barite continuously increase in the residual liquid until the point of oversaturation, at which point bastnäsite-(Ce) and barite begin to crystallize.Depending on the activity of PO 4 , monazite-(Ce) may also crystalize at this stage, however no primary monazite-(Ce) was observed.In the Saint-Honoré complex, the combination of initial low partial melting of the mantle, crystal fractionation and partitioning of REE in the magma instead of the fluid could explain this high concentration of REE and Ba during the magmatic stage.A late-stage fluid, probably orthomagmatic and rich in SiO 2 , FeO, H 2 O, CO 2 , and Cl, dissolves the primary REE minerals and enriches the REE in solution.The initial Fe-carbonatite exhibits evidence of hydrothermal fluid-induced cracking of the primary carbonatite up to the total breakdown in the form of a hydraulic breccia.This hydrothermal brecciation induces spontaneous decreases in pressure and causes the crystallization of tiny needles of REE minerals around clasts and partially dissolved minerals such as barite.Finally, the upper part of the Fe-carbonatite, post-magmatic, is partially leached by surficial (meteoric or marine) water, evidenced by the absence of halite and barite in the uppermost zone and, more extensively, by the strong near-surface weathering.The precipitation of barite, the interactions with carbonates and the crystallization of apatite (Table S2) are all important mechanisms that precipitate REE minerals, such as fluoro-carbonates and monazite-(Ce).These three mechanisms are present in the Saint-Honoré REE deposit and probably played key roles in the hydrothermal precipitation of REE mineralization.

Conclusions
In earlier studies, the REE mineralization at Saint-Honoré has been interpreted to be purely hydrothermal.In this study, access to deeper portions of the carbonatite provided crucial evidence showing magmatic crystallization of bastnäsite-(Ce).Considering the REE distribution is associated with hydrothermal mineralization, the early magmatic crystallization appears critical to generate larger deposits and larger REE crystals which are more profitable to exploit.
Furthermore, the crystallization of primary barite appears to be an important indication of magmatic crystallization of REE minerals.This study and other experiments indicate that carbonatitic magma rich in Ba and REE will crystallize barite and REE minerals.More specific experiments should be conducted to understand whether the concentration of Ba positively affects the solidus of REE minerals or whether the concentration of Ba is a proxy for the concentration of REE in the magma or a combination of both.For future exploration, the presence of primary barite in a carbonatite can be considered to be a good indicator of the carbonatite potential for hosting significant REE deposits.

Figure 1 .
Figure 1.Study site location and local geology of the Saint-Honoré alkaline complex.Map modified from [39].Note that the Saint-Honoré Complex is covered by the limestone and its location is interpreted from diamond drill hole surveys.

Figure 1 .
Figure 1.Study site location and local geology of the Saint-Honoré alkaline complex.Map modified from [39].Note that the Saint-Honoré Complex is covered by the limestone and its location is interpreted from diamond drill hole surveys.

Figure 2 .
Figure 2. Geological map of the Saint-Honoré alkaline complex (below the Ordovician limestone).The location of deepest drill core used for geochemistry is localized by the black dot in the center of the ferro-carbonatite.The figure is modified from Niobec Inc. (Guillaume Matton, personal communication 2018) (UTM NAD 83 projection, Zone 19N).

Figure 2 .
Figure 2. Geological map of the Saint-Honoré alkaline complex (below the Ordovician limestone).The location of deepest drill core used for geochemistry is localized by the black dot in the center of the ferro-carbonatite.The figure is modified from Niobec Inc. (Guillaume Matton, personal communication 2018) (UTM NAD 83 projection, Zone 19N).

19 Figure 4 .
Figure 4. Selected images of the REE mineralization ordered by increasing depth: (A) Scanning electron microscope (SEM) image of low-intensity hydraulic breccia as suggested by the carbonate puzzle-shaped clasts (approximate depth: 235 m); (B) SEM image of a polycrystalline cluster showing an ameboidal shape, radial bastnäsite-(Ce), monazite-(Ce), barite, carbonate, and hematite (approximate depth: 381 m); (C) Photomicrograph (polarized natural light) of a large barite crystal showing a sharp pseudo-hexagonal shape (euhedral) with a pristine core, the external dashed line highlight the barite edge while the internal dashed line separate the recrystallized outer rim from the pristine core (approximate depth: 500 m); (D) Photograph from the recovered drill core showing REE-rich Fe-carbonatite with purple polycrystalline clusters hosting the REE minerals; the buff coloration in the lower right-hand side of the photograph appears on carbonates after a few days of exposure to air and humidity (approximate depth: 500 m); (E) SEM image of bastnäsite-(Ce) inclusions in a barite crystal suggesting early crystallization of REE minerals or co-crystallization of barite and bastnäsite-(Ce) (and parisite-(Ce)) (approximate depth: 910 m); (F) Photograph from the recovered drill core showing barite with a pristine core and an altered rim (approximate depth: 1100 m).Cb = carbonate, Brt = barite, Bsn = bastänasite-(Ce), Hem = hematite.These images are representative of all drill cores and not only the deepest one used for geochemistry.

Figure 4 .
Figure 4. Selected images of the REE mineralization ordered by increasing depth: (A) Scanning electron microscope (SEM) image of low-intensity hydraulic breccia as suggested by the carbonate puzzle-shaped clasts (approximate depth: 235 m); (B) SEM image of a polycrystalline cluster showing an ameboidal shape, radial bastnäsite-(Ce), monazite-(Ce), barite, carbonate, and hematite (approximate depth: 381 m); (C) Photomicrograph (polarized natural light) of a large barite crystal showing a sharp pseudo-hexagonal shape (euhedral) with a pristine core, the external dashed line highlight the barite edge while the internal dashed line separate the recrystallized outer rim from the pristine core (approximate depth: 500 m); (D) Photograph from the recovered drill core showing REE-rich Fe-carbonatite with purple polycrystalline clusters hosting the REE minerals; the buff coloration in the lower right-hand side of the photograph appears on carbonates after a few days of exposure to air and humidity (approximate depth: 500 m); (E) SEM image of bastnäsite-(Ce) inclusions in a barite crystal suggesting early crystallization of REE minerals or co-crystallization of barite and bastnäsite-(Ce) (and parisite-(Ce)) (approximate depth: 910 m); (F) Photograph from the recovered drill core showing barite with a pristine core and an altered rim (approximate depth: 1100 m).Cb = carbonate, Brt = barite, Bsn = bastänasite-(Ce), Hem = hematite.These images are representative of all drill cores and not only the deepest one used for geochemistry.

Figure 5 .Figure 5 .
Figure 5. Selected SEM images of barite, carbonates, and halite by increasing depth: (A) Dolomite (Mg-Cb) with ankerite (Fe-Cb) rims from a depth of 350 m; (B) Relic barite crystals showing dissolution in quartz and calcite in a REE-fluorocarbonate matrix from a depth of about 400 m; (C) A small halite (Hl) crystal with carbonates from a depth of about 400 m; (D) similar texture than b) from a depth of about 650 m; (E) Barite from a depth of about 800 m with elongated bastnäsite-(Ce) crystal; (F) Barite from deeper than (B,D) (about 1100 m) with 50 μm bastnäsite-(Ce) inclusions, some of which are subidiomorph.Cb = carbonate; Mg-Cb = dolomite; Fe-Cb = ankerite; Brt = barite; Bsn = REE-bastänasite-(Ce); Hem = hematite; HI = halite; Qtz = quartz.These images are representative of all drill cores and not only the deepest one used for geochemistry.

Figure 7 .
Figure 7. Abundance of silica, aluminium and potassium (wt %) as a function of depth within the Saint-Honoré Fe-carbonatite.Note that the uppermost 250 m are strongly weathered.

Figure 8 .
Figure 8. BaO and Cl (in wt %) as a function of depth within the Saint-Honoré Fe-carbonatite.Note that the uppermost 250 m are strongly weathered.

Figure 7 . 19 Figure 7 .
Figure 7. Abundance of silica, aluminium and potassium (wt %) as a function of depth within the Saint-Honoré Fe-carbonatite.Note that the uppermost 250 m are strongly weathered.

Figure 8 .
Figure 8. BaO and Cl (in wt %) as a function of depth within the Saint-Honoré Fe-carbonatite.Note that the uppermost 250 m are strongly weathered.

Figure 8 .
Figure 8. BaO and Cl (in wt %) as a function of depth within the Saint-Honoré Fe-carbonatite.Note that the uppermost 250 m are strongly weathered.

Figure 9 .
Figure 9. REE patterns (full range) of the upper REE-rich Fe-carbonatite (Upper RFeC-pink) and the lower portion of the REE-rich Fe-carbonatite (Lower RFeC-blue).Chondrite data from McDonough and Sun [61].

Figure 9 .
Figure 9. REE patterns (full range) of the upper REE-rich Fe-carbonatite (Upper RFeC-pink) and the lower portion of the REE-rich Fe-carbonatite (Lower RFeC-blue).Chondrite data from McDonough and Sun [61].

( 1 )
the brecciation is significantly less intense; (2) the crystal sizes are larger; and (3) the barite crystals feature a pristine core and an altered rim (Figure4C,F).Furthermore, the preliminary results from fluid inclusion gas analyses by solid probe mass spectrometry[75] revealed a progressive decrease in secondary fluid inclusions trapped in carbonates with depth[76].Moreover, the REE pattern of the deeper zone are typical of magmatic REE deposits such as Mountain Pass (Figure10)[73,74].The upper part of the RFeC has a chondrite-normalized REE pattern similar to the Bear Lodge hydrothermal-influenced carbonatite ((La/Lu)cn ≈ 540 [59]; Figure10) while the deeper zone of the RFeC has a chondritenormalized REE pattern more similar to the magmatic Mountain Pass deposit ((La/Lu)cn ≈ 3900; Figure10)[73,74].

Table 2 .
Average REE concentrations and standard deviation of the average for the RFeC.
Note: Std Dev = standard deviation of the mean.

Table 2 .
Average REE concentrations and standard deviation of the average for the RFeC.
Note: Std Dev = standard deviation of the mean.