The Petyayan-Vara Carbonatite-Hosted Rare Earth Deposit (Vuoriyarvi, NW Russia): Mineralogy and Geochemistry

The Vuoriyarvi Devonian carbonatite–ijolite–pyroxenite–olivinite complex comprises several carbonatite fields: Neske Vara, Tukhta-Vara, and Petyayan-Vara. The most common carbonatites in the Tukhta-Vara and Neske-Vara fields are calciocarbonatites, which host several P, Fe, Nb, and Ta deposits. This paper focuses on the Petyayan-Vara field, in which the primary magmatic carbonatites are magnesian. The least altered magnesiocarbonatites are composed of dolomite with burbankite and are rich in REE (up to 2.0 wt. %), Sr (up to 1.2 wt. %), and Ba (up to 0.8 wt. %). These carbonatites underwent several stages of metasomatism. Each metasomatic event produced a new rock type with specific mineralization. The introduction of K, Si, Al, Fe, Ti, and Nb by a F-rich fluid (or fluid-saturated melt) resulted in the formation of high-Ti magnesiocarbonatites and silicocarbonatites, composed of dolomite, microcline, Ti-rich phlogopite, and Fe–Ti oxides. Alteration by a phosphate–fluoride fluid caused the crystallization of apatite in the carbonatites. A sulfate-rich Ba–Sr–rare-earth elements (REE) fluid (probably brine-melt) promoted the massive precipitation of ancylite and baryte and, to a lesser extent, strontianite, bastnäsite, and synchysite. Varieties of carbonatite that contain the highest concentrations of REE are ancylite-dominant. The influence of sulfate-rich Ba-Sr-REE fluid on the apatite-bearing rocks resulted in the dissolution and reprecipitation of apatite in situ. The newly formed apatite generation is rich in HREE, Sr, and S. During late-stage transformations, breccias of magnesiocarbonatites with quartz-bastnäsite matrixes were formed. Simultaneously, strontianite, quartz, calcite, monazite, HREE-rich thorite, and Fe-hydroxides were deposited. Breccias with quartz-bastnäsite matrix are poorer in REE (up to 4.5 wt. % total REE) than the ancylite-dominant rocks (up to 11 wt. % total REE).


Introduction
Rare-earth elements (REEs) are widely viewed as critical metals due to their extensive use in modern and "green" technologies [1,2]. Carbonatites, i.e., igneous rocks composed of more than 50 vol. % primary magmatic carbonates of Ca, Mg, Fe, and Mn and less than 20 wt. % of SiO 2 [3], Figure 1. The geological setting and structure of the Vuoriyarvi alkaline-ultrabasic carbonatite complex simplified from [83]. Blue circles indicate the positions of the Tukhta-Vara, Neske-Vara, and Petyayan-Vara carbonatite fields and the Kolvikskiy agpaitic nepheline-syenite massif-satellite.  [83]. Blue circles indicate the positions of the Tukhta-Vara, Neske-Vara, and Petyayan-Vara carbonatite fields and the Kolvikskiy agpaitic nepheline-syenite massif-satellite. This paper focuses on specific carbonatites from the eastern flank of the massif, which outcrops on the slope of the Petyayan-Vara hill. In this area, the country rocks are melteigites, ijolites, and pyroxenites. A considerable amount of pyroxenites was metasomatically altered, forming coarseand gigantic-grained phlogopite glimmerites, particularly in areas adjacent to carbonatite veins. Carbonatites occur as veins and lenses several hundred meters in length and up to tens of meters in width. In terms of chemical and mineral characteristics, carbonatites of Petyayan-Vara are extraordinary for the Vuoriyarvi massif. Most are magnesiocarbonatites (Figure 2), which we divided into five groups differing in mineral composition (Table 1): 1.
High-Ti magnesiocarbonatites and silicocarbonatites with microcline and areas of apatite mineralization; 3.
The listed rock varieties and their hybrids are the subjects of our study.
Minerals 2020, 10, x FOR PEER REVIEW 3 of 36 This paper focuses on specific carbonatites from the eastern flank of the massif, which outcrops on the slope of the Petyayan-Vara hill. In this area, the country rocks are melteigites, ijolites, and pyroxenites. A considerable amount of pyroxenites was metasomatically altered, forming coarse-and gigantic-grained phlogopite glimmerites, particularly in areas adjacent to carbonatite veins. Carbonatites occur as veins and lenses several hundred meters in length and up to tens of meters in width. In terms of chemical and mineral characteristics, carbonatites of Petyayan-Vara are extraordinary for the Vuoriyarvi massif. Most are magnesiocarbonatites (Figure 2), which we divided into five groups differing in mineral composition (Table 1): 1. Least altered magnesiocarbonatites; 2. High-Ti magnesiocarbonatites and silicocarbonatites with microcline and areas of apatite mineralization; 3. Magnesiocarbonatites with barium-strontium-rare earth mineralization; 4. Late calciocarbonatites; 5. Breccias of magnesiocarbonatites with rare-earth-silicate groundmass.
The listed rock varieties and their hybrids are the subjects of our study.

Sampling Procedure and Analytical Methods
For this study, 43 samples were collected from three geological cross-sections in different areas of the Petyayan-Vara field perpendicular to the strikes of late rare earth magnesiocarbonatite bodies. Five large veins and lenses, whose thicknesses range from 4 to 13 m, were discovered during fieldwork.

Petrographic Studies
Microscopic observations were made on regular and polished thin sections. For the petrographic study, an Axioplan 2 Imaging (Carl Zeiss) microscope (Carl Zeiss, Oberkochen, Germany) was used. Mineral identification was aided by Raman spectroscopy performed using a HORIBA LabRAM HR 800 spectrometer (Horiba Ltd., Kyoto, Japan) with a 514.5 nm (green) exciting laser and a nominal output power of 50 mW.

Sampling Procedure and Analytical Methods
For this study, 43 samples were collected from three geological cross-sections in different areas of the Petyayan-Vara field perpendicular to the strikes of late rare earth magnesiocarbonatite bodies. Five large veins and lenses, whose thicknesses range from 4 to 13 m, were discovered during fieldwork.

Major Elements
The contents of major and volatile components were studied with the classical methods of "wet" chemistry. When analyzing the concentrations of SiO 2 , Al 2 O 3 , Fe 2 O 3 , MgO, TiO 2 , and SrO, average portions of the powered rock were decomposed by melting with borax and soda and were then transferred into the solution. A separate sample, also melted with borax and soda, was used to determine the F content, while Cl was analyzed after sintering with KNaCO 3 . Certain samples were decomposed by the acid in the mixture (H 2 SO 4 + HNO 3 + HF) to analyze the concentrations of K 2 O, Na 2 O, and MnO; in mixture (HNO 3 + HF) to analyze P 2 O 5 ; in mixture (H 2 SO 4 + HF) to analyze FeO; and in nitric acid to calculate S content. To analyze CO 2 content, the tested powder was previously exposed to NaOH. Analyses of SiO 2 , Al 2 O 3 , Fe 2 O 3 , MgO, CaO, and SrO were carried out using the atomic absorption method; those of K 2 O, Na 2 O, and MnO with the emission method; TiO 2 and P 2 O 5 with the colourimetric method; F and Cl via direct potentiometric determination; CO 2 and FeO with titration; and S T , H 2 O + , and H 2 O − gravimetrically. The analytical errors were 1.5% for concentrations >10 wt. % and 3.5% for concentrations between 1 and 10 wt. %.

Trace Elements
The concentrations of trace elements were determined by inductively coupled plasma emission mass spectrometry (ICP-MS) on an ELAN 9000 DRC-e quadrupole mass spectrometer (PerkinElmer Inc., Shelton, CT, USA). For this procedure, each sample was completely transferred into the solution after acid decomposition in glassy carbon crucibles by gradually adding distilled HF (triple addition and evaporation to wet salts) and HNO 3 . The resulting sample solutions were diluted in such a way that the intensities of the analytical signals corresponded to the selected calibration ranges. The precision of the ICP-MS measurements was better than ±5% for most elements.

High-Ti Carbonatites with Areas of Apatite Mineralization
High-Ti carbonatite (with over 1 wt. % of TiO2) is dense red and brown rock with a fine crystalline carbonate matrix. Macroscopically, this rock resembles the classical rødbergite of the Fen massif (Norway) [58,88]. The high-Ti carbonatites of the Petyayan-Vara are localized within selvages of carbonatite bodies and constitute veins several meters in thickness. The morphology of the contacts between these veins and magnesiocarbonatites is irregular but sharp (Figure 4a). Several thin (of about 1 cm and less in thickness) veinlets of high-Ti carbonatite cut least altered magnesiocarbonatites.
In all high-Ti carbonatites of the Petyayan-Vara, the principal minerals are dolomite, sodiumand barium-free microcline, and and Fe-Ti oxides. The mineral abundance of these carbonatites is variable, and quartz, calcite, Ti-rich phlogopite, apatite, Ti-rich aegirine, and albite are present in lesser amounts, or absent. Sulfides, magnetite, zircon, pyrochlore, strontianite, baryte, ancylite-(Ce), and monazite were found as accessory phases. Feldspars and sulfides (mainly pyrite) are widespread only in high-Ti carbonatites (in later metasomatic rocks, sulfur is accommodated almost exclusively

High-Ti Carbonatites with Areas of Apatite Mineralization
High-Ti carbonatite (with over 1 wt. % of TiO 2 ) is dense red and brown rock with a fine crystalline carbonate matrix. Macroscopically, this rock resembles the classical rødbergite of the Fen massif (Norway) [58,88]. The high-Ti carbonatites of the Petyayan-Vara are localized within selvages of carbonatite bodies and constitute veins several meters in thickness. The morphology of the contacts between these veins and magnesiocarbonatites is irregular but sharp (Figure 4a). Several thin (of about 1 cm and less in thickness) veinlets of high-Ti carbonatite cut least altered magnesiocarbonatites.
Potassium feldspar is predominantly a microcline of an unusual goldish-yellow color with no visible twinning in the thin section. Microcline occurs as xenomorphic spongy grains like those of Dol-1 (see Figure 4b). Chemically, this phase is pure (Na, Ba, and Sr are below detection limit of WDS). However, this phase contains abundant fluid and solid-phase microinclusions, giving it a turbid appearance. The microcline penetrates both the groundmass and late mineralized fractures. Dolomite in these rocks is represented by two generations. The earliest generation (Dol-1), comprising the bulk of high-Ti carbonatites, is turbid, spongy, and xenomorphic ( Figure 4b). This generation is chemically heterogeneous due to variability in its Mg/Fe ratio (19 ± 8 mol. % Ank). Its manganese content is low (3 ± 4 mol. % Kut). The spongy dolomite is chemically similar to dolomite Dol-1 from the unaltered magnesiocarbonatites (Fe-rich, without Sr). Dol-1 from high-Ti carbonatites hosts abundant inclusions of Fe oxides and hydroxides (mostly goethite and hematite). The late generation dolomite Dol-2 is adjacent to the selvages of mineralized fractures. Dol-2 differs from Dol-1 in its compact morphology, optical clarity, idiomorphic appearance, and fine oscillatory zoning (Figure 4b). Compared with Dol-1, Dol-2 contains significantly less Mg and Mn, and its content is much less variable (4 ± 3 mol. % Ank, 3 ± 2 mol. % Kut, without Sr). Dol-2 occurs sporadically in the groundmass composed of Dol-1. There, Dol-2 occurs as grains isolated from spongy Dol-1 by rims of titanium oxides (Figure 4c).
Potassium feldspar is predominantly a microcline of an unusual goldish-yellow color with no visible twinning in the thin section. Microcline occurs as xenomorphic spongy grains like those of Dol-1 (see Figure 4b). Chemically, this phase is pure (Na, Ba, and Sr are below detection limit of WDS). However, this phase contains abundant fluid and solid-phase microinclusions, giving it a turbid appearance. The microcline penetrates both the groundmass and late mineralized fractures. Where coexistent with phlogopite, microcline typically replaces the mica pseudomorphically. In high-Ti carbonatites, the bulk of the silica is concentrated in microcline. Moreover, in many high-Ti carbonatites samples, the SiO 2 content exceeds 20 wt. % and are therefore, nominally classified as silicocarbonatites.
High-Ti carbonatites host 5 vol. % and greater titanium oxides (Ti-oxides hereafter), which are intergrown with dolomite (mainly Dol-1), microcline, and phlogopite ( Figure 4d). Ti-oxides occur as three polymorphic varieties: anatase, rutile, and brookite. In high-Ti carbonatites, their total abundance is high (from 0.9 to 4.4 wt. % TiO 2 in the rock, with a median of 2.3 wt. % for 11 samples). The most common polymorph is brookite, while anatase is less abundant. Sparse rutile occurs paramorphically on both minerals. In association with the late generations of Nb-rich Ti-oxides and late dolomite, well-preserved idiomorphic grains of pyrochlore were detected. Typically, these grains' compositions are Ca 1.1 Na 0.9 Fe 0.1 Nb 1.8 Ti 0.2 O 6 F, which corresponds to fluorcalciopyrochlore [89]. Hydroxycalciopyrochlore and kenoplumbopyrochlore are less abundant [90]. Magnesiocarbonatites were overprinted by Ti-rich carbonatitic fluids as evidenced by fine veinlets. Only increased TiO 2 content (up to 1.1 wt. %) differentiates magnesiocarbonatites with veinlets from the other primary dolomite carbonatites, as described in Section 4.1.1. Occasionally, pyrochlore clusters (+ magnetite and Ti-depleted phlogopite) occur without Ti-oxides in magnesiocarbonatites, bordering the high-Ti rocks. In these rocks, the concentration of Nb reaches 1500 ppm, with an average Nb concentration of 400-800 ppm in high-Ti varieties and 10-100 ppm in other carbonatites of the Petyayan-Vara. Niobium deportment in these rocks is shown in the Nb-Nb/Ti diagram ( Figure 5). Where coexistent with phlogopite, microcline typically replaces the mica pseudomorphically. In high-Ti carbonatites, the bulk of the silica is concentrated in microcline. Moreover, in many high-Ti carbonatites samples, the SiO2 content exceeds 20 wt. % and are therefore, nominally classified as silicocarbonatites.
High-Ti carbonatites host 5 vol. % and greater titanium oxides (Ti-oxides hereafter), which are intergrown with dolomite (mainly Dol-1), microcline, and phlogopite ( Figure 4d). Ti-oxides occur as three polymorphic varieties: anatase, rutile, and brookite. In high-Ti carbonatites, their total abundance is high (from 0.9 to 4.4 wt. % TiO2 in the rock, with a median of 2.3 wt. % for 11 samples). The most common polymorph is brookite, while anatase is less abundant. Sparse rutile occurs paramorphically on both minerals. In association with the late generations of Nb-rich Ti-oxides and late dolomite, well-preserved idiomorphic grains of pyrochlore were detected. Typically, these grains' compositions are Ca1.1Na0.9Fe0.1Nb1.8Ti0.2O6F, which corresponds to fluorcalciopyrochlore [89]. Hydroxycalciopyrochlore and kenoplumbopyrochlore are less abundant [90]. Magnesiocarbonatites were overprinted by Ti-rich carbonatitic fluids as evidenced by fine veinlets. Only increased TiO2 content (up to 1.1 wt. %) differentiates magnesiocarbonatites with veinlets from the other primary dolomite carbonatites, as described in Section 4.1.1. Occasionally, pyrochlore clusters (+ magnetite and Ti-depleted phlogopite) occur without Ti-oxides in magnesiocarbonatites, bordering the high-Ti rocks. In these rocks, the concentration of Nb reaches 1500 ppm, with an average Nb concentration of 400-800 ppm in high-Ti varieties and 10-100 ppm in other carbonatites of the Petyayan-Vara. Niobium deportment in these rocks is shown in the Nb-Nb/Ti diagram ( Figure 5).  Figure 5 are grouped as follows: (1) rocks that contain pyrochlore mineralization and are free of Ti-oxides (with high Nb content and Nb/Ti ratio of 0.4-1.2); (2) high-Ti carbonatites and least altered magnesiocarbonatites hosting TiO2-carbonate veinlets (with high Nb content and a Nb/Ti ratio of 0.01-0.03); and (3) rocks free of both pyrochlore and Ti-oxides (with low Nb content and a variable Nb/Ti ratio). The high Nb contents in combination with the low Nb/Ti ratios suggest that Tioxides are the principal niobium concentrators in high-Ti carbonatites (as well as the Petyayan-Vara carbonatites in general), which is typical of the late-stage carbonatites of many other complexes worldwide (e.g., [91][92][93]).

Samples in
In most high-Ti rocks, fluorapatite is scarce and P2O5 content is about 0.2-0.5 wt. %, which is similar to the P2O5 content in other varieties of the Petyayan-Vara carbonatites. In contrast, some areas of high-Ti rocks contain up to 7.5 wt. % P2O5. In these rocks composed of apatite + calcite ± Ti-rich aegirine assemblages, apatite occurs as a vein network intersecting the groundmass. Like the veinlets with Ti-oxides, those with apatite were found in several samples of the least altered magnesiocarbonatites surrounding the high-Ti carbonatites. Apatite forms clusters of anhedral grains (Figure 6a), rims on various minerals (e.g., sulfides, see Figure 6b), and radial aggregates of small euhedral crystals (Figure 6c). Anhedral grains represent the earliest apatite generation (Ap-1) and mainly occur in the carbonate groundmass as apatite-rich spots, while the rims and idiomorphic  Figure 5 are grouped as follows: (1) rocks that contain pyrochlore mineralization and are free of Ti-oxides (with high Nb content and Nb/Ti ratio of 0.4-1.2); (2) high-Ti carbonatites and least altered magnesiocarbonatites hosting TiO 2 -carbonate veinlets (with high Nb content and a Nb/Ti ratio of 0.01-0.03); and (3) rocks free of both pyrochlore and Ti-oxides (with low Nb content and a variable Nb/Ti ratio). The high Nb contents in combination with the low Nb/Ti ratios suggest that Ti-oxides are the principal niobium concentrators in high-Ti carbonatites (as well as the Petyayan-Vara carbonatites in general), which is typical of the late-stage carbonatites of many other complexes worldwide (e.g., [91][92][93]).

Samples in
In most high-Ti rocks, fluorapatite is scarce and P 2 O 5 content is about 0.2-0.5 wt. %, which is similar to the P 2 O 5 content in other varieties of the Petyayan-Vara carbonatites. In contrast, some areas of high-Ti rocks contain up to 7.5 wt. % P 2 O 5 . In these rocks composed of apatite + calcite ± Ti-rich aegirine assemblages, apatite occurs as a vein network intersecting the groundmass. Like the veinlets with Ti-oxides, those with apatite were found in several samples of the least altered magnesiocarbonatites surrounding the high-Ti carbonatites. Apatite forms clusters of anhedral grains (Figure 6a), rims on various minerals (e.g., sulfides, see Figure 6b), and radial aggregates of small euhedral crystals (Figure 6c). Anhedral grains represent the earliest apatite generation (Ap-1) and mainly occur in the carbonate groundmass as apatite-rich spots, while the rims and idiomorphic crystals of apatite (Ap-2) are distributed along the fractures. Anhedral grains of Ap-1 are almost free of impurities (up to a maximum of 0.5 wt. % Na 2 O, 0.8 wt. % SiO 2 , and 0.9 wt. % SrO; the content of other impurities is below the detection limit).
Minerals 2020, 10, x FOR PEER REVIEW 9 of 36 crystals of apatite (Ap-2) are distributed along the fractures. Anhedral grains of Ap-1 are almost free of impurities (up to a maximum of 0.5 wt. % Na2O, 0.8 wt. % SiO2, and 0.9 wt. % SrO; the content of other impurities is below the detection limit). In contrast, the composition of both the rims and idiomorphic crystals of late apatite Ap-2 is highly variable, even within a single grain. On the BSE images, this compositional variability is expressed as fine oscillatory zoning ( Figure 6d). In Ap-2, SO3 content varies in the range of 0.4-1.3 wt. %, SrO concentration reaches 3.0 wt. %, Na2O reaches 3.0 wt. %, and ThO reaches 0.9 wt. %. Total REE2O3 content is 8.8 wt. %, with increased content of heavy REE (up to 0.6 wt. % Gd2O3 and 1.9 wt. % Y2O3). A positive correlation between Sr and S indicates the joint occurrence of these elements in the structure of Ap-2 ( Figure 7a). A similar positive correlation is observed between Na and REE ( Figure 7b). However, the correlation between the contents of (Sr + S) and (Na + REE) is inverse ( Figure 7c). The zoning in Ap-2 can thus be explained by the successive alternation of apatite sites enriched with either Sr and S or Na and REE (the dark and light apatite zones in Figure 6d, respectively). Nevertheless, Sr and S mainly accumulated during the initial stage of Ap-2 growth, and Na and REE mostly accumulated during the later stage (see Figure 6d). In contrast, the composition of both the rims and idiomorphic crystals of late apatite Ap-2 is highly variable, even within a single grain. On the BSE images, this compositional variability is expressed as fine oscillatory zoning ( Figure 6d). In Ap-2, SO 3 content varies in the range of 0.4-1.3 wt. %, SrO concentration reaches 3.0 wt. %, Na 2 O reaches 3.0 wt. %, and ThO reaches 0.9 wt. %. Total REE 2 O 3 content is 8.8 wt. %, with increased content of heavy REE (up to 0.6 wt. % Gd 2 O 3 and 1.9 wt. % Y 2 O 3 ). A positive correlation between Sr and S indicates the joint occurrence of these elements in the structure of Ap-2 ( Figure 7a). A similar positive correlation is observed between Na and REE ( Figure 7b). However, the correlation between the contents of (Sr + S) and (Na + REE) is inverse ( Figure 7c). The zoning in Ap-2 can thus be explained by the successive alternation of apatite sites enriched with either Sr and S or Na and REE (the dark and light apatite zones in Figure 6d, respectively). Nevertheless, Sr and S mainly accumulated during the initial stage of Ap-2 growth, and Na and REE mostly accumulated during the later stage (see Figure 6d). Calcite is associated with apatite in veins and pseudomorphically replaces dolomite near the late idiomorphic apatite crystals (see Figure 4c). In pseudomorphs, calcite commonly hosts inclusions of Fe-hydroxides.
Zircon occurs as microinclusions, anhedral segregations in apatite, and overgrowths on apatite edges. Similar relationships between zircon and sulfate-rich fluorapatite were discovered in other complexes-for example, in the Khaluta carbonatite massif, Russia [79].
As noted above, strontianite, baryte, and ancylite-(Ce) are accessory phases in high-Ti carbonatites. The individual grains and small clusters of these minerals were detected in cavities and fractures of high-Ti carbonatite (Figure 8a), often in combination with late HREE-rich apatite.
Monazite-(Ce) and bastnäsite-(Ce) crystallized later than the Ba-Sr-REE assemblages. In the intersections between the monazite veinlets and cavity-filling Ba-Sr-REE assemblage, ancylite is pseudomorphically replaced by monazite ( Figure 8b). In most cases, monazite is spatially associated with apatite ( Figure 8c). In some samples, both Dol-1 and Dol-2 are partially replaced by monazite ( Figure 8d). Monazite is enriched in CaO, SrO, and BaO (by several percent) and contains up to 5.4 wt. % SO3. This, however, is typical for carbonatites of the Vuoriyarvi massif [94]. Th-bearing monazite (up to 2.3 wt. % ThO) often accompanies Th-bearing apatite. Most analyses of monazite showed a mass deficit. In addition to monazite, rhabdophane is also common in carbonatites [Ce(PO4)•H2O] (e.g., [71,80,[95][96][97]). The proportions of REE:P in this mineral are the same as those in monazite, but rhabdophane analyses show a mass deficit due to the presence of water in its structure. These minerals also differ in structure: Monazite crystallizes in a monoclinic system, space group P21/n, Z = 4; and rhabdophane crystallizes in a hexagonal system, space group P6222, Z = 2. To verify the diagnosis of monazite, we examined its structure by Raman spectroscopy. The Raman spectra of rhabdophane and monazite in the region of 400-1200 cm −1 show the same set of bands with similar positions ( Table 2). Because of the water in the structure of rhabdophane, its spectrum contains additional bands near 3500 cm −1 [98][99][100][101]. The monazite spectra from the Petyayan-Vara carbonatites show no bands in this region ( Figure 3a); therefore, these monazite grains are anhydrous. Spectra of monazite show all bands characteristic of this mineral and most show an additional band of variable intensity at 1088 cm −1 , the main band of calcite [102,103]. Thus, an excess of Ca is most likely a consequence of the monazite-calcite intergrowth, and the mass deficit is caused by the presence of CO3 bound in calcite, which was not considered in the analyses. Calcite is associated with apatite in veins and pseudomorphically replaces dolomite near the late idiomorphic apatite crystals (see Figure 4c). In pseudomorphs, calcite commonly hosts inclusions of Fe-hydroxides.
Zircon occurs as microinclusions, anhedral segregations in apatite, and overgrowths on apatite edges. Similar relationships between zircon and sulfate-rich fluorapatite were discovered in other complexes-for example, in the Khaluta carbonatite massif, Russia [79].
As noted above, strontianite, baryte, and ancylite-(Ce) are accessory phases in high-Ti carbonatites. The individual grains and small clusters of these minerals were detected in cavities and fractures of high-Ti carbonatite (Figure 8a), often in combination with late HREE-rich apatite.
Monazite-(Ce) and bastnäsite-(Ce) crystallized later than the Ba-Sr-REE assemblages. In the intersections between the monazite veinlets and cavity-filling Ba-Sr-REE assemblage, ancylite is pseudomorphically replaced by monazite ( Figure 8b). In most cases, monazite is spatially associated with apatite ( Figure 8c). In some samples, both Dol-1 and Dol-2 are partially replaced by monazite ( Figure 8d). Monazite is enriched in CaO, SrO, and BaO (by several percent) and contains up to 5.4 wt. % SO 3 . This, however, is typical for carbonatites of the Vuoriyarvi massif [94]. Th-bearing monazite (up to 2.3 wt. % ThO) often accompanies Th-bearing apatite. Most analyses of monazite showed a mass deficit. In addition to monazite, rhabdophane is also common in carbonatites [Ce(PO 4 )·H 2 O] (e.g., [71,80,[95][96][97]). The proportions of REE:P in this mineral are the same as those in monazite, but rhabdophane analyses show a mass deficit due to the presence of water in its structure. These minerals also differ in structure: Monazite crystallizes in a monoclinic system, space group P2 1 /n, Z = 4; and rhabdophane crystallizes in a hexagonal system, space group P6 2 22, Z = 2. To verify the diagnosis of monazite, we examined its structure by Raman spectroscopy. The Raman spectra of rhabdophane and monazite in the region of 400-1200 cm −1 show the same set of bands with similar positions ( Table 2). Because of the water in the structure of rhabdophane, its spectrum contains additional bands near 3500 cm −1 [98][99][100][101]. The monazite spectra from the Petyayan-Vara carbonatites show no bands in this region ( Figure 3a); therefore, these monazite grains are anhydrous. Spectra of monazite show all bands characteristic of this mineral and most show an additional band of variable intensity at 1088 cm −1 , the main band of calcite [102,103]. Thus, an excess of Ca is most likely a consequence of the monazite-calcite intergrowth, and the mass deficit is caused by the presence of CO 3 bound in calcite, which was not considered in the analyses.  Bastnäsite occurs as dense, cavity-filling, xenomorphic masses. Sporadically, porous bastnäsite ( Figure 8e) was observed in association with collomorphic quartz filling the latest cracks. Rarely, quartz crystals host bastnäsite microinclusions (Figure 8f). The F content in bastnäsite varies from 0 to 5 wt. %. Although hydroxylbastnäsite-(Ce) is a rare mineral, the Petyayan-Vara field is its type locality [106]. Raman spectroscopy was used to check whether variations in the F content in bastnäsite are related to its substitution with (OH). The spectra of bastnäsite and hydroxylbastnäsite are fairly similar (see Table 2). However, bastnäsite shows one peak in the 1080-1090 cm −1 region and no peaks in the 3500-3600 cm −1 region, while hydroxylbastnäsite shows three peaks in the 1080-1090 cm −1 region and several peaks in the region of 3500-3600 cm −1 due to the presence of an (OH)-group in the latter [104]. Raman spectra of both bastnäsite sensu stricto (with F) and hydroxylbastnäsite (with OH) were observed in high-Ti carbonatites (Figure 9b and Figure 9c, respectively). Further research showed that these two phases coexist in all other bastnäsite-containing varieties of the Petyayan-Vara  Bastnäsite occurs as dense, cavity-filling, xenomorphic masses. Sporadically, porous bastnäsite ( Figure 8e) was observed in association with collomorphic quartz filling the latest cracks. Rarely, quartz crystals host bastnäsite microinclusions (Figure 8f). The F content in bastnäsite varies from 0 to 5 wt. %. Although hydroxylbastnäsite-(Ce) is a rare mineral, the Petyayan-Vara field is its type locality [106]. Raman spectroscopy was used to check whether variations in the F content in bastnäsite are related to its substitution with (OH). The spectra of bastnäsite and hydroxylbastnäsite are fairly similar (see Table 2). However, bastnäsite shows one peak in the 1080-1090 cm −1 region and no peaks in the 3500-3600 cm −1 region, while hydroxylbastnäsite shows three peaks in the 1080-1090 cm −1 region and several peaks in the region of 3500-3600 cm −1 due to the presence of an (OH)-group in the latter [104]. Raman spectra of both bastnäsite sensu stricto (with F) and hydroxylbastnäsite (with OH) were observed in high-Ti carbonatites (Figures 9b and 9c, respectively). Further research showed that these two phases coexist in all other bastnäsite-containing varieties of the Petyayan-Vara carbonatites (see Sections 4.1.3 and 4.1.5). For simplicity, both minerals are hereafter referred to as "bastnäsite".
Both monazite and bastnäsite are often surrounded by or intergrown with oxides and hydroxides of Fe (hereafter, Fe-oxides).  from high-Ti carbonatites of the Petyayan-Vara (red and blue lines). Spectra (black lines) and ID references of minerals are from the RRUFF database [107]. The positions of the main bands (see Table  2) of monazite (m), rhabdophane (r), calcite (c), hydroxylbästnasite (h), and bästnasite (b) are plotted.

Barium-Strontium-Rare Earth Carbonatites
Barium-strontium-rare earth carbonatites are the most abundant rocks among carbonatite veins. Several variations were observed petrographically. Their common feature is the abundance of baryte, ancylite, and/or strontianite. The observed mineral assemblages are named according to the dominant mineral (in the sequence of formation): (1) Baryte, (2) ancylite, and (3) strontianite. The first and second assemblages are the most abundant, and the third was detected at only a few sites.
Minerals of the baryte assemblage reside in fine-grained, brown colored mineralized dissolution cavities (Figure 6a) of the magnesiocarbonatites. There, baryte and late-stage dolomite constitute collomorphic clusters (Figure 6b). Compositionally, late-stage dolomite is similar to the Dol-2 from high-Ti carbonatites (Fe-poor and Sr-free, see Section 4.1.2). A significant part of this assemblage is composed of oxides and hydroxides of Fe and Mn, which commonly occur as zoned spherical aggregates with a collomorphic appearance (Figure 6c). The accessory mineral phases of the baryte assemblage are represented by sulfur-bearing monazite, which is often overgrown by bastnäsite (Figure 6d), as well as hollandite and norsethite. Most bastnäsite is hydrated (hydroxylbastnäsite is Figure 9. Raman spectra of (a) monazite-(Ce), (b) hydroxylbästnasite-(Ce), and (c) bästnasite-(Ce) from high-Ti carbonatites of the Petyayan-Vara (red and blue lines). Spectra (black lines) and ID references of minerals are from the RRUFF database [107]. The positions of the main bands (see Table 2) of monazite (m), rhabdophane (r), calcite (c), hydroxylbästnasite (h), and bästnasite (b) are plotted.

Barium-Strontium-Rare Earth Carbonatites
Barium-strontium-rare earth carbonatites are the most abundant rocks among carbonatite veins. Several variations were observed petrographically. Their common feature is the abundance of baryte, ancylite, and/or strontianite. The observed mineral assemblages are named according to the dominant mineral (in the sequence of formation): (1) Baryte, (2) ancylite, and (3) strontianite. The first and second assemblages are the most abundant, and the third was detected at only a few sites.
Minerals of the baryte assemblage reside in fine-grained, brown colored mineralized dissolution cavities (Figure 6a) of the magnesiocarbonatites. There, baryte and late-stage dolomite constitute collomorphic clusters (Figure 6b). Compositionally, late-stage dolomite is similar to the Dol-2 from high-Ti carbonatites (Fe-poor and Sr-free, see Section 4.1.2). A significant part of this assemblage is composed of oxides and hydroxides of Fe and Mn, which commonly occur as zoned spherical aggregates with a collomorphic appearance (Figure 6c). The accessory mineral phases of the baryte assemblage are represented by sulfur-bearing monazite, which is often overgrown by bastnäsite (Figure 6d), as well as hollandite and norsethite. Most bastnäsite is hydrated (hydroxylbastnäsite is prevalent). Calcite and ancylite occur locally. In samples composed predominantly of this assemblage, the BaO content can reach 17 wt. %.
The ancylite assemblage consists primarily of ancylite (main phase), baryte, strontianite, calcite, and, occasionally, quartz. Rocks with this assemblage have a reddish color due to the abundance of Fe oxides and hydroxides inside the ancylite grains and interstices. Commonly, this mineralization fills either submillimeter veinlets or small cavities. Such cavities were observed in magnesiocarbonatite, baryte-rich carbonatite, and rarely in high-Ti carbonatite. This suggests a stage during which dissolution occurred between the formation of baryte and ancylite assemblages (see inset on Figure 10a). Magnesiocarbonatite samples with filled cavities and veinlets have a total REE oxide content of 1-3 wt. %. In many Petyayan-Vara carbonatite veins, the total REE oxide contents are much higher (≥10 wt. % total REE 2 O 3 ). This enrichment is explained by the presence of (Figure 11a The ancylite assemblage consists primarily of ancylite (main phase), baryte, strontianite, calcite, and, occasionally, quartz. Rocks with this assemblage have a reddish color due to the abundance of Fe oxides and hydroxides inside the ancylite grains and interstices. Commonly, this mineralization fills either submillimeter veinlets or small cavities. Such cavities were observed in magnesiocarbonatite, baryte-rich carbonatite, and rarely in high-Ti carbonatite. This suggests a stage during which dissolution occurred between the formation of baryte and ancylite assemblages (see inset on Figure 10a). Magnesiocarbonatite samples with filled cavities and veinlets have a total REE oxide content of 1-3 wt. %. In many Petyayan-Vara carbonatite veins, the total REE oxide contents are much higher (≥10 wt. % total REE2O3). This enrichment is explained by the presence of (Figure 11a    The accessory phases in this assemblage are bastnäsite, synchysite, monazite, and thorite. Several bastnäsite morphologies were detected. First, peculiar randomly orientated bastnäsite laths and/or tabular bastnäsite syntaxial intergrowths with synchysite occur in late pore-filling quartz and calcite surrounding ancylite-baryte-strontianite aggregates (Figure 12a). Similar syntaxial intergrowths of REE fluorocarbonates occur in many other carbonatite complexes (e.g., [12,14,48,78,80,108]). These intergrowths are commonly interpreted as the results of Ca 2+ and CO3 2− activity changes during mineral formation. Second, we observed the local pseudomorphic replacement of ancylite with anhedral aggregates of bastnäsite (Figure 12b). A similar replacement was found, for example, in carbonatites of Wicheeda, Canada [14] (Figure 7b). Third, small cavities often contain fibroradial bastnäsite aggregates (Figure 12c), which are also common for REE carbonatites worldwide (e.g., [45,68,71,108,109]). Along with monazite and thorite, bastnäsite is one of the latest minerals of the ancylite assemblage. Monazite is chemically similar to the monazite in high-Ti carbonatites and associated with trace thorite (Figure 12d). In the Petyayan-Vara field, thorite is rich in REE (1.5-1.7 wt. % Ce2O3, 1.5-2.0 wt. % Nd2O3, 1.4-1.7 wt. % Sm2O3, 1.8-2.8 wt. % Gd2O3, 1.2-2.8 wt. % Y2O3), and P2O5 (1.3-3.6 wt. %). Fe-oxides have overgrown many primary minerals or fully replaced them. The most pronounced replacement by Fe-oxides is manifested in ancylite grains (see Figure 12c,d). Occasionally, HREE-rich thorite cores are found inside Fe-oxide aggregates ( Figure  12e). The accessory phases in this assemblage are bastnäsite, synchysite, monazite, and thorite. Several bastnäsite morphologies were detected. First, peculiar randomly orientated bastnäsite laths and/or tabular bastnäsite syntaxial intergrowths with synchysite occur in late pore-filling quartz and calcite surrounding ancylite-baryte-strontianite aggregates (Figure 12a). Similar syntaxial intergrowths of REE fluorocarbonates occur in many other carbonatite complexes (e.g., [12,14,48,78,80,108]). These intergrowths are commonly interpreted as the results of Ca 2+ and CO 3 2− activity changes during mineral formation. Second, we observed the local pseudomorphic replacement of ancylite with anhedral aggregates of bastnäsite (Figure 12b). A similar replacement was found, for example, in carbonatites of Wicheeda, Canada [14] (Figure 7b). Third, small cavities often contain fibroradial bastnäsite aggregates (Figure 12c), which are also common for REE carbonatites worldwide (e.g., [45,68,71,108,109]). Along with monazite and thorite, bastnäsite is one of the latest minerals of the ancylite assemblage. Monazite is chemically similar to the monazite in high-Ti carbonatites and associated with trace thorite (Figure 12d). In the Petyayan-Vara field, thorite is rich in REE (1.  Figure 12c,d). Occasionally, HREE-rich thorite cores are found inside Fe-oxide aggregates (Figure 12e). The main minerals of the strontianite assemblage are milky-white strontianite and calcite. These minerals occupy centimeter-scale veins and lenses and lenses that intersect the magnesiocarbonatites with a red ancylite mixture (Figure 13a). Near these veins, the polymineral pseudomorphs of the least altered magnesiocarbonatites have reduced baryte content and a higher abundance of quartz and bastnäsite (Figure 13b). In bastnäsite, the amount of F gradually decreases from core to rim. The rims usually consist of hydroxylbastnäsite (see the inset in Figure 13b).
In proximity to the baryte and ancylite assemblages and strontianite veinlets, dolomite has elevated contents of Fe and Mn, including the occurrences of ankerite (Ank91Dol8Kut1) and kutnohorite (Kut48Dol46Ank6). A similar increase of the Fe and Mn contents in the dolomite of late carbonatites was described for the Sallanlatva massif [110], located 25 km to the west of the Vuoriyarvi massif. The main minerals of the strontianite assemblage are milky-white strontianite and calcite. These minerals occupy centimeter-scale veins and lenses and lenses that intersect the magnesiocarbonatites with a red ancylite mixture (Figure 13a). Near these veins, the polymineral pseudomorphs of the least altered magnesiocarbonatites have reduced baryte content and a higher abundance of quartz and bastnäsite (Figure 13b). In bastnäsite, the amount of F gradually decreases from core to rim. The rims usually consist of hydroxylbastnäsite (see the inset in Figure 13b).
In proximity to the baryte and ancylite assemblages and strontianite veinlets, dolomite has elevated contents of Fe and Mn, including the occurrences of ankerite (Ank 91 Dol 8 Kut 1 ) and kutnohorite (Kut 48 Dol 46 Ank 6 ). A similar increase of the Fe and Mn contents in the dolomite of late carbonatites was described for the Sallanlatva massif [110], located 25 km to the west of the Vuoriyarvi massif.

Late Calciocarbonatites
At sites proximal to the magnesiocarbonatites with abundant ancylite mineralization, giant-grained monomineralic calcite veins are typical. Macroscopically, these veins are distinguishable by their beige-pale brown colours (Figure 13c). In the exocontacts of these veins, magnesiocarbonatites are brecciated and ferruginated. Calcite veins have rhombic dolomite inclusions coated with hydroxide films and fragments of the host rocks (Figure 13d
important REE ore types, while ancylite rocks are less commonly ore grade [57,68], and the ancylitequartz association is extremely rare. The dolomite fragments in the breccias were repeatedly overgrown by later dolomite generations (Figure 14e). The early dolomite is Dol-1 from high-Ti rocks both morphologically (porous) and chemically (9 ± 2 mol. % Ank, 2 ± 1 mol. % Kut, without Sr). The surrounding rims of dolomite Dol-2, optically turbid from the many inclusions of Fe-oxides, are chemically (4 ± 2 mol. % Ank, 2 ± 1 mol. % Kut) similar to dolomite Dol-2 from the other carbonatite varieties of the Petyayan-Vara. The needles and films of Fe oxides/hydroxides were formed on the surfaces of Dol-2 (see Figure 14e). The latest dolomite overgrowths (Dol-3) are oscillatory zoned in the BSE-images due to their Fe/Mg contents and non-porous and optically clear properties. Dol-3 is rich in Fe (14 ± 5 mol. % Ank), very poor in Mn (<1 mol. % Kut), and free of Sr.

"Hybrid Rocks" and the Mineral Paragenesis
In addition to the rocks discussed above, "hybrid" carbonatites occur in the Petyaian-Vara field. In these rocks, the mineral associations described in

"Hybrid Rocks" and the Mineral Paragenesis
In addition to the rocks discussed above, "hybrid" carbonatites occur in the Petyaian-Vara field. In these rocks, the mineral associations described in Sections 4.1.1-4.1.5 are superimposed on each other in various combinations. Therefore, hybrid carbonatites record mineral paragenesis. For example, sample 15K-12.0 (Figure 15a,b) contains minerals and textures representative of most assemblages of the Petyayan-Vara carbonatites (in the formation sequence):
In the space between the dolomite fragments separated by brecciation, minerals were deposited in several stages. The earliest mineral is baryte, representing the assemblage of baryte-dominant carbonatites. Baryte is not as idiomorphic as dolomite (Figure 15c) but is more euhedral than all other minerals; 3.
After baryte and dolomite, euhedral crystals of irregular polygonal and rhombic shapes were crystallized (Figure 15e). This morphology is characteristic of ancylite from the studied carbonatites. Ancylite crystals correspond to the formation stages of ancylite-dominant carbonatites. A detailed study of these crystals showed that they represent pseudomorphs consisting of anhedral bastnäsite fragments inside a matrix of intergrown ancylite and quartz (Figure 15f). Inclusions of ancylite and strontianite also occur in bastnäsite; 4.
Fine crystalline quartz occupies the most remaining space. It contains spongy grains (Figure 15g) and clusters of bastnäsite needles (Figure 15h). These minerals correspond to the mineral assemblage of breccias with a REE and silicate mineral matrix;

5.
On the walls of the cavities that remained after the quartz deposition, quartz crystals with well-crystallized tips were formed (Figure 15i). In the final stage, all remaining empty space was filled with strontianite assemblage.
These observations combined with our data obtained for other carbonatite varieties of the Petyayan-Vara suggest a paragenetic sequence similar to that for the Bear Lodge carbonatites ( [72]; cf. their Table 1 and our Figure 16) and others of the Kola Peninsula ( [78] and references therein).
3. After baryte and dolomite, euhedral crystals of irregular polygonal and rhombic shapes were crystallized (Figure 15e). This morphology is characteristic of ancylite from the studied carbonatites. Ancylite crystals correspond to the formation stages of ancylite-dominant carbonatites. A detailed study of these crystals showed that they represent pseudomorphs consisting of anhedral bastnäsite fragments inside a matrix of intergrown ancylite and quartz (Figure 15f). Inclusions of ancylite and strontianite also occur in bastnäsite; 4. Fine crystalline quartz occupies the most remaining space. It contains spongy grains (Figure 15g) and clusters of bastnäsite needles (Figure 15h). These minerals correspond to the mineral assemblage of breccias with a REE and silicate mineral matrix; 5. On the walls of the cavities that remained after the quartz deposition, quartz crystals with wellcrystallized tips were formed (Figure 15i). In the final stage, all remaining empty space was filled with strontianite assemblage.
These observations combined with our data obtained for other carbonatite varieties of the Petyayan-Vara suggest a paragenetic sequence similar to that for the Bear Lodge carbonatites ( [72]; cf. their Table 1 and our Figure 16) and others of the Kola Peninsula ( [78] and references therein).

Geochemistry
The chemical differences among the Petyayan-Vara carbonatites are illustrated by the primitive mantle-and chondrite-normalized plots (Figure 17a,b).

Geochemistry
The chemical differences among the Petyayan-Vara carbonatites are illustrated by the primitive mantle-and chondrite-normalized plots (Figure 17a,b). All carbonatites of the Petyayan-Vara are rich in Ba and Sr. High-Ti carbonatites are distinguished from other varieties by higher contents of K, Ti, Nb, and Ta (as well as Al, Si, and Fe) due to the abundance of potassium feldspar, phlogopite, and titanium oxides. Apatite-bearing high-Ti carbonatites are additionally enriched with P, Zr, Hf, U, and HREE. The least altered magnesiocarbonatites have LaCN/YbCN of 670 to 3170 and total REE contents of 0.32-2.03 wt. %. High-Ti carbonatites and their apatite-rich areas have low LaCN/YbCN ratios (60-330 and 20-160, respectively) and low REE contents (0.15-0.74 wt. % and 0.18-1.14 wt. % total REE, respectively). Ba-Sr-REE carbonatites (including brecciated magnesiocarbonatites with quartz-bastnäsite matrix) are characterized by highly fractionated REE distributions (LaCN/YbCN = 1200-7600) with a maximum REE content of 5-11 wt. %. For ancylite-dominant rocks, the LaCN/YbCN ratios are the highest (more than 6000), and the total REE content is more than 10 wt. %. Brecciated magnesiocarbonatites with quartzbastnäsite matrix have LaCN/YbCN of 1500-2700. Calciocarbonatites are characterized by low REE contents (0.60-1.14 wt. %) and average LaCN/YbCN ratios (240-500). There is a positive correlation between the LaCN/YbCN ratios and total REE content. The LaCN/YbCN ratio is closely related to the mineral composition. Low LaCN/YbCN values are characteristic of rocks in which HREE-apatite Ap-2 is present. High LaCN/YbCN values were calculated for rocks with bastnäsite and (higher) with ancylite-i.e., with minerals containing mainly LREE. The only independent REE host mineral in calciocarbonatites is ancylite, trapped with wall rock fragments. Studies of the Aley complex (British All carbonatites of the Petyayan-Vara are rich in Ba and Sr. High-Ti carbonatites are distinguished from other varieties by higher contents of K, Ti, Nb, and Ta (as well as Al, Si, and Fe) due to the abundance of potassium feldspar, phlogopite, and titanium oxides. Apatite-bearing high-Ti carbonatites are additionally enriched with P, Zr, Hf, U, and HREE. The least altered magnesiocarbonatites have La CN /Yb CN of 670 to 3170 and total REE contents of 0.32-2.03 wt. %. High-Ti carbonatites and their apatite-rich areas have low La CN /Yb CN ratios (60-330 and 20-160, respectively) and low REE contents (0.15-0.74 wt. % and 0.18-1.14 wt. % total REE, respectively). Ba-Sr-REE carbonatites (including brecciated magnesiocarbonatites with quartz-bastnäsite matrix) are characterized by highly fractionated REE distributions (La CN /Yb CN = 1200-7600) with a maximum REE content of 5-11 wt. %. For ancylite-dominant rocks, the La CN /Yb CN ratios are the highest (more than 6000), and the total REE content is more than 10 wt. %. Brecciated magnesiocarbonatites with quartz-bastnäsite matrix have La CN /Yb CN of 1500-2700. Calciocarbonatites are characterized by low REE contents (0.60-1.14 wt. %) and average La CN /Yb CN ratios (240-500). There is a positive correlation between the La CN /Yb CN ratios and total REE content. The La CN /Yb CN ratio is closely related to the mineral composition. Low La CN /Yb CN values are characteristic of rocks in which HREE-apatite Ap-2 is present. High La CN /Yb CN values were calculated for rocks with bastnäsite and (higher) with ancylite-i.e., with minerals containing mainly LREE. The only independent REE host mineral in calciocarbonatites is ancylite, trapped with wall rock fragments. Studies of the Aley complex (British Columbia, Canada) and the Bear Lodge complex (USA) have shown that calcite in hydrothermal carbonatites usually has a flat REE distribution spectrum [115]. The contribution of REE from calcite is likely responsible for the decreased La CN /Yb CN values for calciocarbonatites.
The interrelationships of element concentrations in the studied rocks were analyzed using factor analysis. The basic principles of this method are available in [116,117]. The analytical dataset of all carbonatite varieties from the Petyayan-Vara field (43 samples), comprising the contents of major and trace elements (Table S1), was subjected to a factor analysis in its R-modification ("by variables"). Factors were extracted by the method of principal components without subsequent rotation. The two major (in terms of weight) factors are: where the subscript at factor (F) is the fraction of the total variance explained by the factor (in %), and the subscript at the element is its loading on the factor multiplied by 100. The numerator contains elements that directly correlate with the factor, and the denominator includes elements with inverse correlation. Only elements with significant factor loadings are shown (the correlation is insignificant at r < |r crit |; for n = 43, the critical value module r crit is 0.30 at a significance level of p = 0.05). On the factor plane determined by these factors, the figurative points of the elements shape the fields that reflect the chemical specificity of the Petyayan-Vara carbonatites ( Figure 18). In order to provide a mineralogical meaning to the geochemical variables, the technique outlined by [118] was used. As demonstrated by this study, in most cases, each component is concentrated primarily in one or two minerals. For example, almost all Ti, Nb, and Ta are in titanium oxides; apatite hosts P, HREE, and, along with phlogopite, F; zircon contains all Zr and Hf; apatite and bastnäsite take a leading role in the distribution of Th; ancylite and strontianite incorporate most Sr; almost the entire volume of LREE is concentrated in ancylite and, to a lesser extent, in bastnäsite.
The considered factors distinguish the main mineral assemblages in Figure 18. The fields of the mineral assemblages overlap in only two cases: (1) When a common mineral is present in two assemblages, e.g., baryte found in both unaltered and Ba-Sr-REE carbonatites; and (2) when an element is present in two or more minerals from different assemblages, e.g., Th contained in both apatite and bastnäsite (the apatite and Ba-Sr-REE assemblages, respectively). The figurative points of such "common" components are shifted towards one of the axes. That is, by virtue of this "commonality", the factor corresponding to the orthogonal axis cannot indicate which of the assemblages to attribute to the element. This is especially important for calcium points, which are displaced to the intersection of the axes due to calcite's involvement in virtually all the identified mineral assemblages.
The first factor separated high-Ti carbonatites from all others. The values of the factor loading indicate that the formation process of high-Ti rocks was accompanied by an increase in K, Na, Al, Si, and the entire spectrum of HFSE (specifically Ti, Nb, Ta, Zr, and Hf). The second factor reflects the processes related to the accumulation of rare-earth elements (REEs). These processes are (1) phosphorus addition to the high-Ti rocks, accompanied by abundant apatite crystallization (hereafter, "apatitization") and (2) Ba-Sr-REE mineralization of the magnesiocarbonatites. Heavy rare earth elements (Y, Dy-Lu) are precipitated mainly during the first process, and light rare earth elements (La-Gd) are precipitated during the second process. Apatitization was also accompanied by a separation of HFSE: Ti, Nb, and Ta are concentrated in high-Ti carbonatites, while Zr and Hf are concentrated only in the apatite-rich sites of these rocks. Notably, the correlation coefficients between the concentration of P 2 O 5 and the concentrations of Zr, Hf, and HREE are high (from 0.85 to 0.93 at n = 43).
Minerals 2020, 10, x FOR PEER REVIEW 22 of 36 Figure 18. Points of the studied major and trace elements on the factor plane formed by the two factors with maximum eigenvalues. The colored fields distinguish the groups of elements associated with certain minerals from one or other mineral assemblages. Illustration is based on the work of [118].
We note that factor analysis was carried out at the early stage of the study (before petrographic and mineralogical investigations). Two data sets were subjected to statistical processing: (1) chemical compositions of rock samples (major and trace elements) and (2) a combined selection of chemical compositions and X-ray diffraction patterns. Factor analysis of a set of chemical compositions revealed groups of elements similarly distributed in the rocks ( Figure 18). The projection of samples onto this factor plane yielded the first classification of carbonatites. Thus, the statistical data processing of a combination of chemical compositions and X-ray diffraction patterns provides the following advantages: (1) simplified determination of minerals, (2) revealing the effect of an identified mineral on the distribution of each element, and (3) tracing the connection between results of the statistical processing of geochemical data and minerals/mineral assemblages ( Figure 18, the relationship between elements and minerals) [118]. The detailed petrographic study presented in this work fully confirmed the results obtained using statistical methods. This shows the high efficiency of the described statistical approaches for the study of such complex geological objects as late carbonatites.

The Protolith of Carbonatites of the Petyayan-Vara Field
Except for the Petyayan-Vara field, in all parts of the Vuoriyarvi massif, the most common primary carbonatites are calciocarbonatites. It is believed that magnesiocarbonatites, which dominate in the Petyayan-Vara field, are more mature products of magmatic differentiation than Figure 18. Points of the studied major and trace elements on the factor plane formed by the two factors with maximum eigenvalues. The colored fields distinguish the groups of elements associated with certain minerals from one or other mineral assemblages. Illustration is based on the work of [118].
We note that factor analysis was carried out at the early stage of the study (before petrographic and mineralogical investigations). Two data sets were subjected to statistical processing: (1) chemical compositions of rock samples (major and trace elements) and (2) a combined selection of chemical compositions and X-ray diffraction patterns. Factor analysis of a set of chemical compositions revealed groups of elements similarly distributed in the rocks ( Figure 18). The projection of samples onto this factor plane yielded the first classification of carbonatites. Thus, the statistical data processing of a combination of chemical compositions and X-ray diffraction patterns provides the following advantages: (1) simplified determination of minerals, (2) revealing the effect of an identified mineral on the distribution of each element, and (3) tracing the connection between results of the statistical processing of geochemical data and minerals/mineral assemblages ( Figure 18, the relationship between elements and minerals) [118]. The detailed petrographic study presented in this work fully confirmed the results obtained using statistical methods. This shows the high efficiency of the described statistical approaches for the study of such complex geological objects as late carbonatites.

The Protolith of Carbonatites of the Petyayan-Vara Field
Except for the Petyayan-Vara field, in all parts of the Vuoriyarvi massif, the most common primary carbonatites are calciocarbonatites. It is believed that magnesiocarbonatites, which dominate in the Petyayan-Vara field, are more mature products of magmatic differentiation than calciocarbonatites [119]. Consequently, the Petyayan-Vara magnesiocarbonatites were localized in the most highly evolved part of the Vuoriyarvi complex.
The magnesiocarbonatites with minor baryte and hexagonal polymineral pseudomorphs (consisting of REE-carbonates, baryte, and strontianite; see Section 4.1.1) are only slightly metasomatically altered and are probably most representative of the least altered carbonatites in the Petyayan-Vara field. The presence of carbocernaite, burbankite, and calcite inclusions in poikilitic dolomite are typical of primary magmatic dolomite ( [120], cf. their Figure 5e and our Figure 3d). The polymineral pseudomorphic replacement textures are evidence of alteration. As shown for the various alkaline complexes, such pseudomorphs are usually formed after burbankite [71,78,121], an early syn-magmatic mineral [97]. Thus, the presence of altered burbankite is of interest, because the "early magmatic crystallization of rare earth minerals could be a significant factor in generating high-volume REE deposits" [16].
Magmatic magnesiocarbonatites were protoliths for all other (carbo)-hydrothermal-metasomatic carbonatite varieties. The specific mineralogical and geochemical features (including high concentrations of REE) were produced in these carbonatites during a series of metasomatic events.

Formation of HFSE-Rich (High-Ti) Carbonatites
High-Ti carbonatites are of particular interest due to the discovery of the unique Morro dos Seis Lagos Ti-Nb-REE deposit (Brazil) [11], wherein Ti-oxides are the main hosts of Nb. High-Ti carbonatites are reported at Magnet Cove (USA) [91], Gross Brukkaros (Namibia) [92], and Salpeterkop (South Africa) [93]. In terms of mineral associations, the high-Ti carbonatites of the Petyayan-Vara are akin to the "classical" rødbergites of the Fen massif [58,88]. However, there are some minor differences. First, in contrast to rødbergites, feldspars and mica from high-Ti carbonatites of the Petyayan-Vara do not contain barium, and baryte is minor. Second, unlike high-Ti carbonatites, rødbergites of the Fen massif lack titanium oxides. High-Ti carbonatites may represent specific varieties of rødbergites, but the mechanism of their formation and the source of their enrichment in the HFSE are not completely clear.
At present, we cannot state that magnesiocarbonatites were the only protolith of high-Ti carbonatites of the Petyayan-Vara field. Titaniferous carbonatites could possibly have formed by the (auto) metasomatism of earlier carbonatites of a different variety (not magnesiocarbonatite). Nonetheless, the formation of high-Ti carbonatites was obviously controlled by the influence from an additional portion of a fluid or fluid-saturated melt.
Geochemically, high-Ti carbonatites of the Petyayan-Vara are unusual for carbonatites. Commonly, carbonatites contain minor amounts of K, Al, Si, and Ti. According to the results of the study of natural materials [122] and the experimental data [123][124][125][126], the low contents of the indicated elements are caused by their partitioning into the silicate part during immiscible separation of carbonatite and silicate melt. Moreover, the crystallization of alkaline feldspars, which are widespread in high-Ti carbonatites, requires abnormally high (for carbonatite melts) SiO 2 activity [9,127]. This high SiO 2 activity is explained by assimilation and/or by superimposed hydrothermal alteration [9,127,128]. This is consistent with the data on Sr and Nd isotope systems in silicocarbonatites [129]. General ideas on the nature of high-Ti carbonatites correspond to the above-mentioned hypotheses.
The source of titanium (as well as K, Al, and Si) could be a specific residual silicate-carbonatite melt and/or fluid. The presence of Ti-rich fluorophlogopite and fluorocalciopyrochlore in high-Ti carbonatites of the Petyayan-Vara field indicates that the source was rich in fluorine. According to reports on the partitioning of elements between silicate melts and immiscible fluoride, chloride, carbonate, phosphate, and sulfate melts [125], this type of residual melt most likely accumulated HFSE, including titanium. An alternative source of Ti, K, Al, and Si could be melt-assimilated or fluid-remobilized silicate rock (s). For example, Kapustin [130] established a compositional relationship between carbonatites and associated aluminosilicate rocks, concluding that Ti and several other elements could be assimilated from host rocks. However, in the Petyayan-Vara field, the assimilation of wall rock material seems to be improbable. The pyroxenites hosting carbonatites are almost free of Al and K. Na, which is abundant in pyroxenites, is absent in high-Ti carbonatites. Hypothetically, the neighbouring gneisses could also have been assimilated by carbonatite melts deeper in the complex, increasing Al and Si concentrations. As shown in [131], contamination of the kimberlite melt with crustal components results in the appearance of a Ti-rich potassium-aluminosilicate melt, captured as the inclusions in diamonds. Ti-rich melts of similar compositions were obtained in experiments modelling the melting of carbonatized pelites [132].

Connection between the Local HREE Enrichment and Apatitization
In general, the REE content in high-Ti carbonatites is relatively low (up to 1.1 wt. % total REE). However, high-Ti carbonatites and especially their apatite-rich areas are noticeably enriched in HREE compared to all other rocks of the Petyayan-Vara field, including Ba-Sr-REE carbonatites. Late carbonatites with HREE mineralization associated with late-stage hydrothermal phenomena are found in many complexes (e.g., Tundulu, Malawi [80]; Songwe-Hill, Malawi [47]; Bear Lodge, USA [68]; Lesser Qinling, China [34]; Fen, Norway [58]). Researchers have presented several hypotheses for the HREE enrichment of late carbonatites. Possible explanations include HREE enrichment of the primary source, the passive enrichment of minerals due to the selective removal of LREE, and selective LREE or HREE transport by fluids rich in sulfur, phosphorus, chlorine, and/or fluorine. More than one mechanism is possible for the enrichment of late carbonatites with heavy rare earth elements. As a result of these mechanisms, xenotime-(Y) [34,69,133], brockite-(Y) [69], mckelveyite-(Y), ewaldite, donnayite-(Y) [78], and other HREE mineral phases were formed in these rocks. In carbonatites of the Petyayan-Vara and a number of other complexes (e.g., [47,80,81]), the accumulation of HREE is related to fluorapatite (see Section 4.2).
In the Petyayan-Vara carbonatite field, apatitization occurs in high-Ti carbonatites and the surrounding rocks. However, the early apatite Ap-1, which precipitated either simultaneously or immediately after the primary mineral assemblage of the high-Ti carbonatites, is free of impurities. The activity of both LREE and HREE was thus relatively low during the formation of high-Ti carbonatites. The HREE enrichment of the Petyayan-Vara carbonatites is related to alteration and late-stage apatite crystallization. The morphology of late Ap-2 rich in REE is typical of hydrothermal apatite ( [134]; cf. Figure 10h therein and Figure 5c in this contribution). Nevertheless, the solubility of REE in phosphate solutions is low ( [135,136] and the references therein). Therefore, phosphate solutions are unlikely transport agents of REE. According to the experimental data, "the main REE transporting ligands are chloride and sulfate, whereas fluoride, carbonate, and phosphate likely play an important role as depositional ligands" [136].
HREE-rich Ap-2 appears to have formed by in situ dissolution of Ap-1 and recrystallization. The chemical composition of Ap-2 (high Sr, REE, and S concentrations) relates it to the minerals of Ba-Sr-REE carbonatites (baryte, strontianite, and REE-carbonates). The high sulfur content in Ap-2 suggests a higher activity of this component during crystallization ( [137] and the references therein). Moreover, the fluids controlling the formation of late apatite influenced the dolomite of the protolith in the same way as the fluids, which caused ancylite formation. In both cases, dolomite was pseudomorphically replaced by the calcite hosting many inclusions of Fe oxides (cf. Figures 6c and 11d). This illustrates the relationship between the HREE-enrichment of high-Ti carbonatites and the stage of alteration which led to the formation of baryte-and ancylite-dominant carbonatites in the Petyayan-Vara field.
We suggest that the mechanism for HREE precipitation was the interaction of orthomagmatic carbonatitic Ba-Sr-REE-S-rich fluid with phosphate-bearing rocks formed simultaneously with high-Ti carbonatites or immediately after them. As a result of exposure to this fluid, dolomite was locally replaced with calcite, the original apatite Ap-1 was partially dissolved and reprecipitated in situ, and Ba-Sr-REE minerals crystallized in cracks and cavities. Most Ba, Sr, and LREEs were removed through the cracks, whereas HREEs were absorbed by the growing Ap-2. A similar enrichment mechanism was described in [138].
The same HREE enrichment mechanism may have occurred in carbonatites of other complexes geochemically similar to those of the Petyayan-Vara. For instance, in the Montviel Carbonatite Complex (Abitibi, Canada), the Ba-Sr metasomatic zone with LREE mineralization is separated from the P-F metasomatic zone containing HREE minerals [13].

The Main Stage of REE Mineralization (the Formation of Ancylite Ores)
From an economic standpoint, carbonatites of this field are unique because ancylite is the main REE ore mineral. In most known carbonatite complexes, the principal REE host is bastnäsite (REE-carbonate). Frequently, bastnäsite occurs in association with strontianite (Sr-carbonate). However, despite the abundance of both Sr and REE in late carbonatites, ancylite (REE-Sr-carbonate) forms large accumulations on rare occasions. Unlike bastnäsite [139][140][141], the physicochemical stability limits of ancylite are not clear. Factors contributing to the accumulation and preservation of economically significant ancylite contents include the following: 1.
The absence of supergene processes, which often destabilize ancylite. For example, in the Bear Lodge alkaline complex, a clear vertical zonation was discovered [68,71]. Ancylite-strontianitebaryte mineralization occurs in deeper parts of the section and is replaced by REE fluorocarbonates in the shallower oxide zone (100 m from the surface). Carbonatite samples of the Petyayan-Vara field were collected from the surface, where most of their mineral assemblages are similar to those found only in the deeper parts of the Bear Lodge alkaline complex. The weak influence of supergene alteration on the rocks of the Vuoriyarvi massif is explained by specific aspects of the recent geological history of the Kola region. The Svecofennian shield has undergone several periods of glaciation. Any weathering crusts that may have formed before the latest glaciations were eroded by glacial abrasion. Subsequently, the exposed rocks were protected from aggressive chemical weathering by the subarctic climate. This facilitated the preservation of ancylite ores in the near-surface parts of the Vuoriyarvi massif.
As mentioned above, chloride and sulfate are considered to be the main ligands for REE transport in hydrothermal fluids (e.g., [14,62,[144][145][146]). We suggest that the Ba-Sr-REE-rich sulfate (+Cl?) fluid played an important role in the transport of REE within the Petyayan-Vara field.
Baryte-ancylite assemblages seem to offer reliable evidence for the correlation between REE enrichment and the involvement of sulfur in the system. Probably carbonatites of Petyayan-Vara derived from a specific sulfur-rich source, similar to that described in the West Transbaikalia carbonatite province [79]. A large baryte deposit (244.5 million tons of ore with 11-14 wt. % of BaSO 4 and 0.4 wt. % of REE 2 O 3 , [83]) of the Sallanlatva carbonatite complex is located approximately 20 km to the west of the Vuoriyarvi massif. Thus, some sulfur-rich sources seem to have had a regional impact on carbonatite complexes. The lithological control of sulfur mineralization is observed in the studied area (i.e., sulfides occur in high-Ti rocks, and sulfates are distributed in REE-Sr-Ba rocks). The transition from sulfide mineralization to sulfate mineralization presumably was caused by a change in oxidizing conditions in the fluid during the waning stages of metasomatism. As shown for the carbonatites from the Bear Lodge, the reason for this phenomenon could be an influx of external (meteoric) fluid [70,72].
The mineral composition of carbonatites at Petyayan-Vara does not reflect the influence of chlorine. However, barium is poorly transported by sulfur-rich fluid due to precipitation of baryte by the reaction of Ba 2+ with sulfate [147]. In the empirical modelling study of [147], baryte solubility was conditioned by Cl content (the main factor) in the fluid, as well as temperature and pressure. Thus, the fluid most likely contained both sulfur and chlorine.
The similar mineral composition of burbankite pseudomorphs and ancylite assemblage indicates their formation during a single event. We suggest that S (+Cl?)-rich fluid caused decomposition of burbankite from primary magmatic carbonatites, thereby accumulating Sr, Ba, and REE, and then re-deposited these elements into cracks and crushing zones. This explains the specific appearance of ancylite-dominant rocks (see Figure 11a-c). As showed Andersen et al. [72], fluids responsible for the burbankite-to-ancylite conversion are predominantly magmatic (carbohydrothermal).
The formation of the baryte mineral assemblage instead of the ancylite assemblage is explained by the mechanism proposed in [36]. Initially, fluids had a high capacity to transport REE, and, consequently, no REE minerals were deposited. Cooling caused the crystallization of a large amount of baryte [147], which captured the bulk of the S from the fluid. As a result, the capacity of the latter to transport REE decreased, triggering ancylite precipitation in the waning stage of (carbo)-hydrothermal activity. Evidence of a decrease in sulfur content in the fluid is the absence of celestine in the ancylite assemblage. Additionally, both the drop in pressure when the fluid infiltrated the fractured and brecciated zones and the pH increase under fluid-rock interactions could contribute to the deposition of ancylite [136].
The origin of coarse-grained calcite carbonatites, which are associated with ancylite ores, is enigmatic. They may have been derived from the Ba-Sr-REE fluid phase depleted in REE, Sr, and Ba due to the precipitation of strontianite, ancylite, and baryte. Alternatively, these rocks could have originated from calciocarbonate melt residua injected after the cessation of the ancylite crystallization. The first scenario is supported by the high aCa 2+ during Ba-Sr-REE mineralization, in which dolomite is replaced by calcite near ancylite-rich areas. Syntaxial intergrowths of bastnäsite and synchysite also indicate an increase in aCa 2+ [80,148]. However, the juxtaposition of Ba-Sr-REE carbonatites and late calciocarbonatites may be incidental. In this case, the latter are products of barren Ca-rich fluid circulated during the final stage of the Vuoriyarvi complex's formation.

Late-Stage Processes
Late-stage processes yielded breccias of magnesiocarbonatite with a matrix of quartz-bastnäsite, monazite, Fe hydroxides, quartz, and late dolomite and calcite generations. Of particular interest are the breccias of magnesiocarbonatite with a quartz-bastnäsite matrix. In addition to alkaline-ultrabasic carbonatite complexes, similar rocks, dominated by quartz and REE-fluorocarbonates occur in other types of alkaline massifs (for example, in the Jbel Boho alkaline complex, Morocco [148]; Fazenda Varela carbonatites of the Lages alkaline complex, Brazil [149]). The origin of such rocks is attributed to the remobilization of rare-earth elements from alkaline rocks by the late-stage fluids [148]. At the Palabora complex, it was shown that thorium can be used as an indicator of remobilization [65]. In the Palabora carbonatites, Th is abundant in early minerals crystallized during the main stage of REE mineralization. During remobilization, Th was incorporated into the secondary phases. In the Petyayan-Vara carbonatites, the early Th-bearing phase is apatite of the second generation crystallized during the early stage of the Ba-Sr-REE metasomatism. In minerals of the main Ba-Sr-REE deposition stage, Th content is below the detection limit. Bastnäsite contains up to several percent of Th. Apparently, the quartz-bastnäsite matrix in the dolomite carbonatite breccias also formed as a result of remobilizing the REE from preexisting minerals (primarily ancylite).
As shown above, ancylite is often pseudomorphically replaced by bastnäsite (Figures 12b, 14d and 15f). Similar pseudomorphs associated with late hydrothermal processes were found in the Bear Lodge carbonatites [68,71]. We detected the replacement of ancylite by bastnäsite exclusively in quartz-bearing carbonatites of Petyayan-Vara. This mineral replacement was obviously associated with an increase in silica activity. A likely source of silica is the host silicate rock.
In Bear Lodge carbonatites, it was shown that strontium is lost during the replacement of ancylite by bastnäsite [71]. A similar process could explain the strontianite assemblage in Ba-Sr-REE carbonatites and the bulk of strontianite in hybrid rocks (see Section 4.1.6). In other complexes, strontianite-dominant rocks are rather rare (e.g., [19,150]). From our observations of the paragenetic relationships, strontainite-dominant rocks are some of the latest rocks to form (e.g., Figure 7 in [150]).
Ancylite was not abundant in magnesiocarbonatite which was protolithic to breccias with a quartz-bastnäsite matrix. The unique matrix mineralogy is a result of remobilized REEs introduced with the silica-rich fluid. The absence of sulfur in the breccias indicates the negligible role of S in REE mobility during the later hydrothermal stages. A specific feature of the Petyayan-Vara carbonatites is the predominance of hydroxylbastnäsite. The listed characteristics point to low fluorine activity and high (OH) − activity in the fluid. These conditions favor thorium transfer [151]. As shown in the inclusion study of quartz-synchysite rocks in the Jbel Boho alkaline complex, the most probable REE transport ligand in such fluids is Cl [148].
Near early-crystallized apatite (Ap-1) and in apatite-free rocks, monazite is poor in Th but rich in Ba and S. Monazite in proximity to Ap-2 has higher thorium content. This substitution and the spatial association of apatite and monazite indicate the growth of the latter at the expense of apatite. A similar mechanism of monazite formation was described for the late-stage carbonatites of the Cummins Range Carbonatite Complex, Australia [7], the Fen complex, Norway [58], the Sokli complex, Finland [152], the Palabora complex, South Africa [65], and others worldwide. Presumably, thorite or huttonite (complementary to monazite) is the product of the selective concentration of elements (Th and HREE) released during apatite decomposition. A similar selective redistribution of LREE and Th-HREE has been described in carbonatites of Nolans Bore (Australia) [8].
Collomorphic monazite (see Figure 10d) was formed by the decomposition of apatite during a late hydrothermal stage. In addition to collomorphic textures, needle and fibrous monazite forms, characteristic of supergene REE mineralization (e.g., [71,153]), are present (e.g., Figure 12c). The grain morphologies of Fe oxides/hydroxides are analogous to those of monazite (collomorphic, needle, and fibrous). Some of the Fe oxides and hydroxides were likely formed by the classical decomposition scheme of Fe-containing carbonates [88,154] and trapped as inclusions within newly formed carbonates. However, a substantial amount of Fe migrated during late hydrothermal and, possibly, supergene stages. In our study area, ferruginization accompanied all late-stage processes.

Conclusions and Further Perspectives
The main conclusions of this study are as follows: 1.
The primary rocks of the Petyayan-Vara field were magmatic magnesiocarbonatites, consisting of dolomite with burbankite, minor magnetite, and possibly sulfides. These rocks were the protolith of all other carbonatites; 2.
Mineralogical and geochemical features, including high REE contents that are characteristic of the rocks of the Petyayan-Vara field, are the result of several metasomatic events; 3.
The introduction of a carbonatitic fluid (or fluid-saturated melt) rich in K, Al, Si, Fe, and Ti has resulted in the formation of high-Ti carbonatites containing Nb-rich Ti-oxides. Local apatitization is superimposed on high-Ti carbonatites and the surrounding rocks; 4.
The most pervasive alteration of magnesiocarbonatites was caused by a magmatic S (+Cl?)-rich fluid with high concentrations of Ba, Sr, and REE. At this stage, baryte-strontianite-ancylite ores with high REE contents (up to 11 wt. % total REE) were formed; 5.
Initially, Ba-Sr-REE-S fluid had a high capacity to transport REE, and, consequently, no REE minerals were deposited. Cooling caused the crystallization of a large amount of baryte, which captured the bulk of the S from the fluid. As a result, the capacity of the latter to transport REE decreased, triggering ancylite precipitation in the waning stage of (carbo)-hydrothermal activity. This mechanism explains the formation of the baryte mineral assemblage instead of the ancylite assemblage;

6.
Interaction with the Ba-Sr-REE-S fluid caused replacement of early-formed apatite with a secondary HREE-rich generation of apatite; 7.
Late-stage processes yielded breccias of magnesiocarbonatites with a quartz-bastnäsite matrix, minerals of the strontianite assemblage, monazite, Fe hydroxides, quartz, as well as late dolomite and calcite generations.
The similarity of the mineralogical and geochemical features of the Petyayan-Vara carbonatites and late carbonatites of various complexes worldwide suggests common mechanisms of formation. A number of issues are yet to be clarified: • The universality of the proposed HREE-enrichment mechanisms for late carbonatites; • Ancylite stability field limits; • The mechanism(s) of joint fluid migration of REE and Ba.
A clear correspondence between the results of petrographic and statistical studies points out the high efficiency of the described statistical approaches for the study of such complex geological objects as late carbonatites.