Front. Earth Sci. Frontiers in Earth Science Front. Earth Sci. 2296-6463 Frontiers Media S.A. 1518194 10.3389/feart.2025.1518194 Earth Science Original Research Genesis of the Bayan Obo Fe-REE deposit: evidence from Mg isotope and geochemistry of siderite carbonatites, Inner Mongolia, China Ren et al. 10.3389/feart.2025.1518194 Ren Chenghao 1 Jia Hongxiang 1 * Zhou Qunmao 2 She Hongquan 1 Li Jinwen 1 Li Baoliang 1 1 Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing, China 2 Research Institute of Exploration and Development, PetroChina Xinjiang Oilfield Company, Karamay, Xinjiang, China

Edited by: Hongjian Zhu, Yanshan University, China

Reviewed by: Debbrota Mallick, University of Georgia, United States

Yanick Blaise Ketchaya, Hefei University of Technology, China

Jiangtao He, Hebei GEO University, China

Wenbin Jia, Chinese Academy of Geological Sciences (CAGS), China

Moise Luemba Luemba, University of Kinshasa, Democratic Republic of Congo

Anupam Banerjee, Indian Institute of Technology Kanpur, India

 *Correspondence: Hongxiang Jia, cugbjhx@126.com 27 10 2025 2025 13 1518194 28 10 2024 22 09 2025 Copyright © 2025 Ren, Jia, Zhou, She, Li and Li. 2025 Ren, Jia, Zhou, She, Li and Li

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The Bayan Obo Fe-REE deposit, a world-renowned giant polymetallic ore concentration area, is located on the northern edge of the North China Craton, which borders on the Central Asian orogenic belt. Many studies have been conducted since the primary orebody was discovered, but its genesis has always been a fiercely debated issue. In recent years, the rapid development of Mg isotope research has shown promise in tracing mantle-derived igneous rocks and carbonatites, providing a new approach to studying the formation of the Bayan Obo deposit. In this study, we conducted Mg isotope composition and whole-rock geochemical analysis on the siderite carbonatites closely associated with the formation of the Bayan Obo deposit to reveal its genesis. Results show that the siderite carbonatites have high levels of Sr (796.3–6,343.1 ppm), Ba (84.8–11593.0 ppm), and Mn (8,319.8–48680.8 ppm), as well as rare earth element abundance (2,696.8–20763.6 ppm), distinguishing them from sedimentary carbonate rocks. The δ26Mg variation range for the siderite carbonatites is −1.31 to −0.09‰, with a mean value of −0.37‰, similar to characteristics of dolomite carbonatites in the mining area. Integrating the results of this study with previous research, we propose that the Bayan Obo deposit formed through the intrusion of carbonatitic magma and subsequent metasomatic processes. The genesis of the deposit was primarily governed by mantle-derived carbonatitic magmatism, with progressive magma differentiation and evolution likely facilitating the enrichment of Fe, REE, Sr, Ba, and Y.

 siderite carbonatite Mg isotope whole-rock geochemistry igneous carbonatite genesis of ore deposit Bayan Obo section-at-acceptance Economic Geology 1 Introduction

Bayan Obo is the world’s largest rare earth element (REE) deposit, with a cumulative exploration of 1,160 Mt of iron ore, and 112 Mt of rare earth oxides (Chao et al., 1997; Xie et al., 2019). Its REE resources account for more than two-thirds of all known rare earth resources in the world. Except for iron and REE, the deposit is also rich in critical minerals such as niobium (one of the world’s largest niobium deposits), scandium, thorium, and fluorine. As a result, Bayan Obo has a significant economic and strategic position, drawing the interest of many geologists throughout the world.

The deposit was discovered in 1927 and has since undergone systematic geological investigations and, in places, mining development. Significant progress has been achieved in understanding its genesis, mineralization age, rare earth mineralogy, and other aspects, providing important contributions to REE mineralization theory (Zhang and Tao, 1996; Zhang P. S. et al., 1998; Zhang P. S. et al., 2001; Yang and Le Bas, 2004; She et al., 2023). The formation of the Bayan Obo deposit has undergone significant structural deformation and pervasive hydrothermal alteration, leading to complex geological and geochemical features. Numerous studies have been conducted on its metallogenic age, hydrodynamic history, mineralization processes, and genetic models (Smith et al., 2015; Fei et al., 2019). Previous studies have attributed the ore fluids of the Bayan Obo deposit to be source from sedimentary (Meng, 1982; Wei and Shangguan, 1983), sedimentary transformation - metasomatism (Lai, 2013; Zhang Y. X. et al., 1998; Zhang et al., 2008; Hou, 1989; Qiao et al., 1997), submarine volcanic exhalation sedimentation (Yuan et al., 1991; Yuan et al., 1995; Bai and Yuan, 1985; Bai et al., 1996), volcanic institutions (Hao et al., 2002; Xiao et al., 2006; Xiao et al., 2012), subduction fluid metasomatism (Ling et al., 2013; Ling et al., 2014), high-temperature hydrothermal and mantle fluid metasomatism (Cao et al., 1994), carbonate rock magma hydrothermal fluids (Zhou et al., 1980; Liu, 1986; Wang et al., 2002), and carbonatite aegirinization (Wang et al., 2010; Wang et al., 2018; Elliott et al., 2018; Liu et al., 2018; Yang et al., 2019). In the last two decades, the relationship between mineralization and carbonatites has received increasing attention, particularly in the Bayan Obo deposit. Earlier studies generally regarded the ore-bearing dolomite as a sedimentary dolostone, but subsequent work has reinterpreted it as an igneous dolomitic carbonatite (Yang et al., 2009; Hu et al., 2023; Li et al., 2024). In parallel, two contrasting models have been proposed to explain the evolution of ore-forming fluids: (1) one links the genesis of REE mineralization to diagenetic processes associated with the emplacement of sedimentary dolostone, whereas (2) the other attributes the mineralization to igneous carbonatite intrusions in the region. Consequently, the origin and evolution of ore-forming fluids at Bayan Obo remain debated and continue to be a key issue in understanding the deposit’s formation.

In the last two decades, increasing attention has been paid to the genetic relationship between mineralization and carbonatites, particularly regarding the reinterpretation of the Bayan Obo ore-bearing dolomite. Early studies generally considered this rock to be a sedimentary dolostone. However, subsequent research based on petrography, geochemistry, and isotopic evidence has demonstrated that it is better classified as an igneous dolomitic carbonatite related to carbonatite magmatism (Yang et al., 2009; Hu et al., 2023; Li et al., 2024). Nevertheless, debate persists regarding the diagenetic features of these ore-bearing dolomite carbonatites, which remains a key issue constraining the understanding of the formation and evolution of the Bayan Obo deposit.

Magnesium is one of the major elements in igneous carbonatites. As high-precision testing technologies for Mg stable isotopes rapidly advance (Galy et al., 2001; Galy et al., 2002; Chang et al., 2003), Mg isotope geochemistry has become widely used in geological processes, such as recreating ancient marine habitats (Gothmanna et al., 2017; Riechelmann et al., 2016), tracking material cycles (Yang et al., 2012; Wang et al., 2016), and tracing the formation of mineral deposits (Ling et al., 2013; Tang et al., 2018). The advances and applications of Mg isotopes have been reviewed (He et al., 2008; Ge and Jiang, 2008; Teng, 2017). For example, Ke et al. (2011) were the first to summarize the Mg isotope composition of various material reservoirs on Earth and planets, and summarized the fractionation mechanism of Mg isotopes under high and low temperature conditions; Chen et al. (2021) mainly summarized the fractionation behavior of Mg isotopes in different geological processes and their application in carbonate rock research, and summarized the possible genesis mechanism of low δ26Mg values in mantle derived rocks; Liu et al. (2023) systematically reviewed the application of Mg isotopes in the genesis of mineral deposits, porphyry genesis, geological thermometers, and fractionation mechanisms in endogenous geological processes. Therefore, Mg isotopic composition provides valuable information about the mineralization process and may be utilized to trace the origin of mineral deposits.

Despite that the Bayan Obo mining region is well known for siderite carbonatites associated with Fe, REE, and Nb mineralization, detailed investigations of these rocks are limited, and Mg isotopes have not previously been employed to elucidate their sources or evolutionary history. In this study, we present the first integrated investigation of Mg isotopes together with major-, trace-, and rare earth element geochemistry of siderite carbonatites, supported by field and petrographic observations. These results provide new insights into the petrogenesis of siderite carbonatites and offer important constraints on the genetic mechanisms of the Bayan Obo deposit.

 2 Geological background

The Bayan Obo Fe-REE deposit is situated around 150 km north of Baotou City in Inner Mongolia, on the northern border of the North China Craton, adjacent to the Central Asian orogenic belt. The region has undergone complex tectonic evolution since the Proterozoic era (Zhao et al., 1999; Zhao et al., 2005; Zhai et al., 2015), including multiple rift events from the late Paleoproterozoic to the Neoproterozoic (Zhai et al., 2015), multi-stage magmatic processes related to subduction in the Paleozoic era (Zhang et al., 1994; Ling et al., 2014; Gao, 2022), and metamorphosis and deformation processes (Zhang et al., 2003; Smith et al., 2015). The long-term tectonic evolution has resulted in complex geological and geochemical characteristics of mineral deposits and has also led to widespread controversy over the understanding of the mineralization background, mineralization era, and mineralization process of mineral deposits.

The Bayan Obo deposit is developed in the Middle Proterozoic Bayan Obo Group, which is divided into nine lithological sections (H1-H9), from bottom to top, H1-H7 is mainly quartz sandstone, feldspathic quartz sandstone, quartzite, with marl in the upper part of the layer, H8 is mainly dolomite, H9 is mainly potassium-rich slate. The orebodies are mainly hosted in the dolomite of the H8 lithological section, forming an east-west ore belt with a length of about 18 km and a width of about 2–3 km. Within a range of 48 km2, there are three mining areas from east to west (Figure 1). The two largest known as the Main Orebody and the East Orebody, are strongly affected by sodium and fluorine alteration, and rare earth, niobium, and iron mineralization are also very strong. In comparison, the western half of the ore-bearing zone, also known as the West Orebody, exhibits relatively weaker hydrothermal alteration and mineralization. Large areas of granite basement are exposed approximately 0.5–1 km east and south of the ore-bearing zone. The REE mineralization mainly occur in dolomite and it is hosted in fluorocarbonate series minerals and monazite (Zhang and Tao, 1996). Hundreds of carbonate dikes are distributed around the periphery of the mining region, intruding into the metamorphosed sedimentary rocks and basal gneiss layers of the Bayan Obo Group. According to the primary mineral composition, they are classified as dolomitic, calcitic, and dolomite-calcite coexisting carbonate dikes. Many scholars believe that the carbonate magmatism represented by carbonate dikes is intimately connected to Fe and REE mineralization (e.g., Liu et al., 2018; Yang et al., 2019). In the southern part of the Bayan Obo deposit, a granite intrusion which is in direct contact with the dolomite was dated to be emplaced around 281–262 Ma (Ling et al., 2014). Associated with this intrusion, a magnesium silicate skarn is developed, locally enriched in REE mineralization.

Regional geological map of the Bayan Obo deposit, Inner Mongolia. (modified from Li et al., 2023).

Geological map displaying various rock formations and structural features of a region, labeled with codes like Jxb, Pt1γo, and Chd. It includes elements such as synclines, anticlines, ore bodies, and faults. A legend at the bottom identifies colors and patterns representing rock types and geological structures. The map covers an area near Bayan Obo Town, marked for scale with a 2 km reference line. 3 Sampling and analytical methods 3.1 Sample collection and petrography

At Bayan Obo, Mg is hosted predominantly in dolomite, with only minor amounts in silicate minerals such as aegirine and biotite. For this study, siderite carbonatites were selected for detailed investigation and compared with micritic mound dolomite, typical sedimentary carbonate rocks, and wall-rock dolomites.

Fresh siderite carbonatite samples were collected from outcrops in the Main and West Orebodies of the Bayan Obo deposit. To minimize contamination, sampling avoided weathered surfaces, hydrothermal veins, and late-stage calcite stringers. Six representative specimens were obtained at 0.5–1 m intervals. Each sample was divided into two portions: one prepared as thin sections for petrographic observation, and the other processed for geochemical and isotopic analysis. For powder preparation, weathered rinds were removed, and the fresh material was cleaned, oven-dried, and reduced to millimeter-sized chips with a jaw crusher or hand hammer. The chips were subsequently ground in an agate mortar or ball mill, sieved to 200 mesh, and reground where necessary. The resulting powders were homogenized, split into aliquots, and sealed in clean containers to ensure uniformity and minimize contamination.

The siderite carbonatite is yellow-brown in color and displays a blocky, porphyritic texture (Figure 2). Siderite phenocrysts account for ∼90% of the primary mineral assemblage, whereas secondary phases include magnetite, pyrite, monazite, and calcite (Figures 3a–d). Rare earth elements are hosted mainly in monazite. The groundmass consists largely of Mg-rich siderite formed through hydrothermal alteration and metasomatism, with partial replacement by dolomite. Hydrothermal fluids enriched in Mg, Ca, Ba, F, and REE promoted decomposition of primary siderite and subsequent formation of Mg-bearing siderite, dolomite, magnetite, and monazite. Magnetite commonly occurs along fractures in carbonate minerals, accompanied by strong alteration and metasomatism. Progressive replacement of siderite results in a mineral sequence from siderite → Mg-siderite → Mg-rich siderite → dolomite → monazite ± magnetite. Back-scattered electron (BSE) imaging reveals corresponding zoning, with decreasing Fe contents expressed as a transition from bright to darker gray tones. Hydrothermal decomposition of siderite phenocrysts also generated pyrite, sericite, and quartz (Figures 3e,f), a process analogous to alteration patterns in porphyry Cu systems.

Photographs of a representative field sample (a,b) and hand specimen (c,d) of siderite carbonatites from the Bayan Obo deposit, exhibiting a porphyritic texture. The phenocrysts are predominantly composed of siderite, while the matrix consists of dolomite, magnetite, and monazite.

Four-panel image of rocks: (a) Rock with an orange and black marker for scale; (b) Rock next to a hammer with a blue handle; (c) Rock with handwritten markings; (d) Rock with more visible markings, showing varying textures and colors.

Photomicrographs in cross-polarized transmitted light (a,c), reflected light (b,d), and scanning electron microscope backscattered-electron images (e,f) of siderite carbonatites from the Bayan Obo deposit, highlighting key geological features, ore mineral assemblages, and textures within the mineralization stages. (Abbreviations: Mag, magnetite; Mnz, monazite; Sid, siderite; Dol, dolomite; Py, pyrite; Bst, bastnäsite; Whitney and Evans, 2010).

Six microscopic images labeled a to f show mineral compositions in rock samples, with annotations indicating the presence of magnetite, monazite, siderite, dolomite, pyrite, and rare earth elements (REE). Images a, c, and d display colorful, patterned textures, while b, e, and f have darker, more muted tones. Each image is marked with scales indicating magnification and dimensions.

In this study, carbonatites are classified according to their dominant mineralogy as calcite, dolomite, and siderite carbonatites, consistent with widely used nomenclature in previous studies of the Bayan Obo deposit (e.g., Woolley and Kempe, 1989; Chao et al., 1997). These mineralogical terms correspond broadly to the geochemical categories of calcio-, magnesio-, and ferro-carbonatites (Le Bas, 1987), but we adopt the mineral-based terminology here because the classification of our samples is supported directly by petrographic observations and mineral compositions.

 3.2 Analytical methods 3.2.1 Labware cleaning

All labware and sample containers were rigorously cleaned to prevent contamination. Specifically, we now clarify that all Teflon vials were pre-cleaned by sequential soaking in 6 N HCl and Milli-Q water, while all plastic ware was soaked in 10% HCl and thoroughly rinsed before use. All acids used were either commercially ultrapure (Optima grade) or purified further by double sub-boiling distillation.

 3.2.2 Major and trace elements

Whole-rock major and trace element analyses of Bayan Obo siderite carbonatites were carried out at Yanduzhongshi Geological Analysis Laboratories Ltd., Beijing, China. Major elements were determined by X-ray fluorescence spectrometry (Zetium, PANalytical Malvern, United Kingdom) with analytical uncertainties better than 1%. Trace elements were measured by inductively coupled plasma–mass spectrometry (ICP-MS), with accuracies better than 5% relative to the GSR-2 standard. For volatile and ultra-trace elements, analytical errors did not exceed 10%.

 3.2.3 Chemical purification and MC-ICP-MS analysis of Mg isotopes

All chemical procedures were conducted in a class-1000 clean laboratory under a class-100 laminar flow hood. Mineral separates or whole-rock powders (10–50 mg, depending on Mg contents) were digested in sealed Teflon beakers or high-pressure bombs using a mixture of ultrapure HF–HNO3 (3:1, v/v), prepared either from commercial ultra-pure acids or by double sub-boiling distillation. After initial digestion, the solutions were evaporated to dryness, treated with aqua regia, dried again, and refluxed with concentrated HNO3 to eliminate residual fluorides. The final residues were dissolved in 1 N HNO3 for column chromatography. No visible precipitates were present in the final solutions, and centrifugation was applied where necessary to ensure complete dissolution without measurable Mg loss.

Magnesium purification was carried out using Savillex microcolumns loaded with 2 mL of Bio-Rad AG50W-X12 (200–400 mesh) cation-exchange resin (Chang et al., 2003). Prior to use, resins were sequentially rinsed with 8 N HNO3 and Milli-Q water (18.2 MΩ) at least three times and stored in ultrapure water. Before each separation, the resins were backwashed with Milli-Q water to remove air bubbles and minimize compaction, then pre-cleaned with alternating rinses of 4 N HNO3 + 0.5 N HF and Milli-Q water, and finally conditioned with 6 mL of 2 N HNO3. Mg was purified using a two-stage cation-exchange procedure with AG50W-X12 resin. In the first stage, Mg and Na were separated from other cations on a long resin column by elution with 2 N HCl. In the second stage, Na was removed from Mg on a short resin column by sequential elution with 0.4 N HCl (to wash away Na) and 6 N HCl (to elute Mg). The purified Mg fractions were collected, evaporated to dryness, and converted into an HNO3 medium. Prior to MC-ICP-MS analysis, the solutions were diluted to appropriate concentrations in 0.1 N HNO3. The overall Mg recovery through the purification procedure was >99.99%.

Magnesium isotope ratios were measured using a Nu Plasma high-resolution multicollector ICP-MS (MC-ICP-MS). Samples introduced through a DSN-100 desolvation nebulizer were mixed with 0.3 mol/L HNO3 prior to analysis. Sample and standard solutions were matched in concentration (∼1 μg/ml Mg), with differences maintained within 10%. Background signals were measured before each analytical block, and data were acquired as sets of ten ratios with 10 s integration times. Instrumental mass bias was corrected by sample–standard bracketing, a procedure that minimizes matrix effects and ensures accurate and reproducible Mg isotope results. The instrument control system is comparable to that of the Nu Plasma HR-MC-ICP-MS at Oxford University (Belshaw et al., 2000). Magnesium isotope compositions are reported in standard δ notation as: δ i Mg = MG i / MG 24 sample / MG i / MG 24 DSM 3 − 1 × 1000 , i = 25  or  26 , relative to the DSM3 standard, with uncertainties reported at the per mil (‰) level.

 3.2.4 Data quality

To ensure the reliability of concentration and Mg isotope compositions, procedural blanks were routinely measured and used to correct all sample measurements. Both concentration and isotopic data were corrected for blank contributions. The total blank was <10 ng for Mg, which is negligible compared to the sample Mg mass (typically 1–5 μg). The isotopic composition of the blanks was indistinguishable from the bracketing standard within analytical uncertainty. Both concentration and isotopic data were corrected for blank contributions using: R corrected = R measured × M sample   –   R blank × M blank / M sample   –   M blank , where R represents the isotopic ratio and M denotes the Mg mass of the sample or blank. The propagated uncertainties include analytical repeatability, blank corrections, and instrumental uncertainties. Long-term bracketing standards analyzed on the MC-ICP-MS demonstrated reproducibility better than ±0.07‰ (2SD) in δ26Mg over the analytical period. Procedural standards processed alongside samples confirmed >99% Mg recovery with no detectable isotope fractionation. The combined uncertainties are reported for all concentration and isotopic measurements.

 4 Analysis results 4.1 Major element oxides

The major element compositions of the samples are shown in Table 1. Siderite carbonatites are dominated by TFe2O3, MgO, CaO, and MnO, with TFe2O3 content between 8.83% and 53.35%, MgO content between 6.60% and 17.86%, CaO content between 4.54% and 26.58%, and MnO content between 1.03% and 4.94%, while SiO2, Al2O3, Na2O, K2O, TiO2 and P2O5 content are low. Compared with global average carbonatites, Bayan Obo siderite carbonatites are notably Mg-enriched and Ca-deficient. Their chemical composition is broadly similar to dolomitic carbonatites and wall-rock dolomites, but differs significantly from sedimentary dolomites and micritic mound dolomites in the deposit. Previous studies classified carbonatites as calcareous, magnesian, or ferritic carbonatites based on CaO, MgO, and TFe2O3+MnO contents (Woolley and Kempe, 1989). As shown in Figure 4, Bayan Obo siderite carbonatites predominantly fall into the ferro (iron-rich) carbonatite category.

The major element contents (wt.%) of siderite carbonatite, carbonatite dyke, dolomite carbonatite, sedimentary carbonate rock and microcrystalline mound dolomite from the Bayan Obo deposit.

Sample types Samples SiO2 TiO2 Al2O3 TFe2O3 MnO MgO CaO Na2O K2O P2O5 LOI F S Total Data sources
Siderite carbonatite BY23207B 1.91 0.04 0.71 53.35 4.94 6.6 4.54 0.3 0.22 0.06 26.72 - - 99.38 This Paper
B2133D 3.01 0.37 0.92 15.11 2.83 13.89 25.77 0.3 0.07 1.24 36.43 - - 99.93
BY21801D 3.1 0.07 0.94 30.44 3.16 13.39 16.24 0.3 0.28 0.45 31.61 - - 99.97
BY23116 2.02 0.01 0.95 8.83 1.03 17.86 26.58 0.3 0.15 0.26 41.26 - - 99.25
BY23005D 3.7 0.21 0.87 27.08 3.13 14.49 16.61 0.3 0.25 0.2 33.1 - - 99.94
B5334-1 0.55 0.08 0.69 35.32 3.82 11.22 16.4 0.3 0.04 0.37 30.1 - - 98.88
B5335-2 2.04 0.09 0.79 36.93 3.76 13.29 14.63 0.3 0.21 0.58 25.42 - - 98.03
Carbonatite dyke BN-09-13 1.61 0.02 0.14 1.59 0.36 8.82 15.61 0.26 0.07 7.6 25.68 0.12 1.03 61.8 Sun et al. (2012)
BN-09-35 6.59 0.12 0.98 8.24 1.46 8.45 30.56 0.13 0.73 3.99 33.37 0.24 0.01 94.6
Dolomite carbonatite BN-09-44 0.71 0.01 0.05 3.58 0.66 18.64 29.69 0.07 0.03 0.54 44.28 0.04 0.01 98.3
BE-09-16 7.7 0 0.15 8.56 2.86 10.9 22.3 0.97 0.34 0.44 31.04 0.89 1.03 85.3
BE-09-22 1.93 0.13 0.16 13.72 1.69 12.86 23.96 0.25 0.12 0.38 34.06 1.82 0.53 89.2
BE-09-30 0.22 0.18 0.06 26.25 2.38 5.08 29.75 0.18 0.02 0.61 20.88 9.36 0.37 85.6
Sedimentary carbonate rock BN-09-64 1.77 0.05 0.08 7.4 0.96 12.02 30.3 0.17 0.09 4.75 33.39 1.82 0.58 91
BN-09-68 26.07 0.07 1.33 0.96 0.07 14.94 21.64 0.02 0.44 0.04 33.75 - - 99.3
Microcrystalline mound dolomite BH-09-14 2.18 0.01 0.14 0.59 0.09 20.78 29.75 0.04 0.04 0.01 45.04 - - 98.7
BH-09-17 2.76 0.01 0.07 0.63 0.08 20.64 29.74 0.05 0.02 0.01 45.03 - - 99

Carbonatite classification diagram of CaO-MgO-Fe2O3 T + MnO (Woolley and Kempe, 1989; Le Maitre, 2002) for siderite carbonatites from the Bayan Obo deposit.

Ternary diagram illustrating the classification of carbonatites by their chemical composition, with axes for CaO, MgO, and TFe₂O₃ + MnO. The diagram is divided into fields for magnesiocarbonatite, ferrocarbonatite, and calciocarbonatite. Blue squares represent siderite carbonatite samples distributed across these fields. 4.2 Trace and rare earth elements

The trace and rare earth element compositions of the Bayan Obo siderite carbonatites differ markedly from those of sedimentary carbonate rocks and are broadly comparable with dolomite carbonatites (Table 2). The total REE (ΣREE) ranges from 2,696.75 to 20,763.57 ppm. The REE patterns are right-skewed, indicating pronounced light rare earth elements (LREE) enrichment (La/YbN = 129.15–1,681.51) and heavy rare earth elements (HREE) depletion. Ce shows slight positive anomalies, while Eu exhibits moderate negative anomalies, suggesting variable oxidation states during mineral formation (Figure 5b). Trace elements such as Ba, Th, and Nb show significant positive anomalies, likely reflecting contributions from accessory minerals (e.g., monazite, bastnäsite), whereas Rb, U, Ta, and Zr display strong negative anomalies (Figure 5a). Compared with sedimentary carbonate rocks, the siderite carbonatites are notably enriched in Mg, Fe, and LREE, consistent with their magmatic-hydrothermal origin. The geochemical characteristics provide insights into the fractionation processes and metallogenic environment of the deposit.

Analysis results of rare earth and trace elements (10–6) in siderite carbonatite from the Bayan Obo deposit.

Smples La Ce Pr Nd Sm Eu Ga Gd Tb Dy Ho Er Tm Yb Lu Y ∑REE LREE
BY23207B 1,088.03 4,548.17 917.47 3,236.05 157.37 19.93 32.36 216.66 16.13 19.77 3.19 18.04 0.54 2.62 0.3 76.54 10,244.27 9,967.02
B2133D 4,131.79 9,701.5 1,340.81 4,723.83 340.61 44.57 59.98 375.46 29.2 38.17 5.25 27.19 0.75 4.03 0.41 124.12 20,763.57 20,283.11
BY21801D 1951.33 4,298.62 704.79 2,333.44 133.54 20.13 30.84 182.1 11.4 17.06 2.31 13.92 0.37 1.9 0.21 66.47 9,671.12 9,441.85
BY23116 1,351.76 1,928.65 265.17 743.66 59.1 12.76 14.68 70.18 4.93 8.05 0.83 4.51 0.11 1.01 0.11 15.44 4,450.83 4,361.1
BY23005D 607.43 1,102.89 195.19 667.12 49.88 8.48 9.79 46.19 4.1 7.86 1.06 4.98 0.22 1.22 0.13 26.36 2,696.75 2,630.99
B5334-1 426.71 2,895.3 866.89 3,309.24 174.88 20.22 29.16 196.5 13.46 14.49 2.38 17.52 0.47 2.37 0.26 64.73 7,940.69 7,693.24
B5335-2 3,821.09 6,582.51 857.23 2,629.85 128.15 19.74 38.46 217.34 12.3 14.68 1.99 14.81 0.29 1.63 0.2 38.65 14,301.81 14,038.57
HREE LREE/HREE (La/Yb) N δEu δCe Rb K Ba Th U Nb Ta Sr Nd P Zr Hf Ti
277.25 35.95 297.88 0.33 1.12 5.61 1783.15 4,050.82 234.1 0.5 996.96 0.11 796.3 3,236.05 266.953 1.35 0.29 120.35
480.46 42.22 735.42 0.38 1.01 1.15 585.031 117.78 184.23 0.85 3,370.34 2.22 956.32 4,723.83 5,424.074 12.02 1.83 1732.43
229.27 41.18 736.68 0.39 0.9 9.91 2,301.9 84.76 78.17 0.04 229.05 0.13 852.12 2,333.44 1941.472 7.11 1.17 252.79
89.73 48.6 960.02 0.61 0.79 4.55 1,205.94 2,465.53 123.65 0.32 39.64 0.05 6,343.1 743.66 1,147.368 2.61 0.2 26.83
65.76 40.01 357.14 0.54 0.79 8.87 2,110.83 137.35 25.14 0.07 495.38 0.68 913.64 667.12 866.579 7.62 1.23 810.9
247.45 31.09 129.15 0.33 1.17 1.25 344.135 7,161.79 173.87 0.18 297.57 0.03 1,412.13 3,309.24 1,604.601 4.09 0.62 278.45
263.24 53.33 1,681.51 0.36 0.89 9.41 1765.6 11,593.01 168.84 0.03 95.61 0.02 961.63 2,629.85 2,518.053 2.82 0.65 258.33

Linear diagram of δ25MgDSM3 versus δ26MgDSM3 for siderite carbonatites from the Bayan Obo deposit.

Scatter plot showing a positive linear relationship between δ²⁶Mgₒₛₘ₃ (x-axis) and δ²⁵Mgₒₛₘ₃ (y-axis). The line of best fit has the equation y = 0.493x + 0.0046 with R² = 0.9973, indicating a strong correlation. Data points are closely aligned with the line. 4.3 Mg isotopes

Test data correctness is often checked by comparing numerical values to established reference materials. For Mg isotopes, however, no internationally recognized standard currently exists due to limited data. All δ26Mg and δ25Mg values in this study fall along the mass-dependent fractionation line (Figure 6), indicating negligible isobaric interferences during MC-ICP-MS analysis. Mg isotope results for the samples are listed in Table 3, with δ26Mg ranging from −1.31 to −0.09‰ and δ25Mg from −0.64 to −0.08‰. Figure 7 places these results in context, comparing them with previously reported ranges for carbonate wall rocks (δ26Mg = −1.67–−1.50‰), dolomites (−2.31–−1.05‰), dolomite carbonatites (−1.12–−0.31‰), sedimentary carbonate rocks (−1.53–−1.51‰), and mantle rocks (−0.52–0.23‰) (Galy et al., 2002; Chang et al., 2003; Tipper et al., 2006b; Wiechert and Halliday, 2007; Handler et al., 2009; Bourdon et al., 2010; Teng et al., 2010a; Yang et al., 2012; Sun et al., 2012; Sun, 2013).

Distribution of Mg isotope compositions of samples from the Bayan Obo deposit. (Data sources, a - Galy et al., 2002; Chang et al., 2003; Tipper et al., 2006b; b - Wiechert and Halliday, 2007; Handler et al., 2009; Bourdon et al., 2010; Teng et al., 2010a; Yang et al., 2012; c - Sun et al., 2012; Sun, 2013).

A scatter plot showing isotopic data of various carbonate types. Dolomite (Previous data-a) is marked with blue diamonds, Peridotite (Previous data-b) with yellow diamonds, and other carbonates like Microcrystalline dolomite, Mesoproterozoic sedimentary dolomite, and Carbonatite dykes are represented in different colors. The x-axis represents δ²⁶Mg₍DSM3₎ values ranging from -2.5 to 0. "This study" indicates data with light blue diamonds.

Mg isotopic compositions of siderite carbonatite, carbonatite dyke, dolomite carbonatite, sedimentary carbonate rock and microcrystalline mound dolomite from the Bayan Obo deposit.

Sample types Samples δ26/24MgDSM3‰ 2SD δ25/24MgDSM3‰ 2SD n Data sources
Siderite carbonatite BY23207B −0.26 0.05 −0.14 0.04 3 This Paper
B2133D −1.31 0.05 −0.64 0.04 3
BY21801D −0.19 0.01 −0.09 0.02 3
BY23116 −0.09 0.05 −0.04 0.03 3
BY23005D −0.2 0.04 −0.09 0.03 3
B5334-1 −0.21 0.03 −0.08 0.02 3
B5335-2 −0.3 0.02 −0.15 0.01 3
Carbonatite dyke BN-09-13 −0.14 0.18 −0.05 0.08 4 Sun et al. (2012)
BN-09-35 −0.34 0.11 −0.16 0.05 7
Dolomite carbonatite BN-09-44 −0.81 0.06 −0.2 0.03 12
BE-09-16 −1.12 0.2 −0.56 0.11 2
BE-09-22 −0.31 0.09 −0.14 0.04 10
BE-09-30 −1.11 0.09 −0.55 0.06 4
Sedimentary carbonate rock BN-09-64 −1.5 0.11 −0.76 0.08 5
BN-09-68 −1.67 0.03 −0.85 0 2
Microcrystalline mound dolomite BH-09-14 −1.99 0.04 −1.01 0 2
BH-09-17 −1.93 0.13 −0.96 0.05 2

Sr + Ba vs. REE scatter diagram for the siderite carbonatites from the Bayan Obo deposit.

Graph illustrating the relationship between (Sr+Ba) and REE in parts per million. Three zones are outlined: endogenetic rocks, carbonate rocks including siderite carbonatite indicated by blue squares, and sedimentary and metamorphic rocks. The x-axis represents REE/ppm, and the y-axis represents (Sr+Ba)/ppm, both on logarithmic scales from 1 to 10^5. 5 Discussion 5.1 Features of the siderite carbonatite source area

The origin of carbonatites has been a subject of debate, reflecting the complexity of the source region’s composition. Some scholars suggest that certain carbonatites may not directly originate from mantle-derived melts, but instead form through a contact metasomatic process within a permeable magma-hydrothermal system (Lentz, 1999). However, the mainstream view generally holds that carbonatites primarily result from low-degree partial melting of the mantle, and differentiation processes such as fractional crystallization or hydrous melt immiscibility (Le Maitre, 2002; Jones et al., 2013; Lentz, 2014; Vladykin and Pirajno, 2021). The siderite carbonatites in this study have similar geochemical characteristics to the dolomite carbonatites (Figure 7), suggesting that they may have originated from the same magma source. Dolomite carbonatites have an ε Nd(t) value ranging from −2.72 to −0.51, and a low (87Sr/86Sr)i value (0.70341–0.70590) (Zhang Z. Q. et al., 2001), indicating little contamination from lower crust minerals (Simonetti et al., 1995; Hou et al., 2006). Research suggests that the genesis of carbonatites with εNd<0 may be related to the mixing of High U/Pb mantle unit - Enriched mantle I (HIMU-EM I) endmembers (Simonetti et al., 1995; Tilton et al., 1998). Therefore, the source area of the Bayan Obo siderite carbonatites may be a complicated mantle source generated by the mixing of HIMU and EM II mantle end components.

Fluid metasomatism in subducting marine sediments is thought to be the cause of rare earth enrichment in mantle source area associated with carbonate rocks (Hou et al., 2015). The Bayan Obo siderite carbonatites exhibit high Sr (796.3–6,343.1 ppm), Ba (84.8–11,593.0 ppm), and LREE (2,631.0–20,283.1 ppm) (Figure 8), as well as a high Ba/Th ratio (0.64–68.66) characteristics, which supports the fluid metasomatism model in sediments. Diving marine sediments are rich in REE and CO2 fluids, which typically have high concentrations of large ion lithophile elements (LILE) and relatively low quantities of high field strength elements (HFSE), resulting in high LILE/HFSE ratios (Turner et al., 1997; Elburg et al., 2022). This type of CO2 rich fluid metasomatism produces carbonate rocks enriched in LILE and REE. During partial melting, HFSE, HREE, and others are retained in the titanium iron oxide, whereas LILE and LREE preferentially penetrate the melt (Foley et al., 2000; Gaetani et al., 2008; Hammouda et al., 2009). This explains the reason for the enrichment of rare earth elements in the mantle source region.

(a) Chondrite-normalized REE pattern diagram and (b) primitive mantle-normalized trace element spider diagram for the siderite carbonatites from the Bayan Obo deposit (Data sources: Ling et al., 2011; Sun et al., 2012).

Two logarithmic scale graphs compare geochemical data. Graph (a) shows Sample/Chondrite ratios for various elements with siderite carbonatite data in blue squares. Graph (b) presents Sample/Primitive Mantle ratios, highlighting dolomite carbonatite in yellow and sedimentary carbonate rock in grey, with element values labeled below. 5.2 Mg isotope system in nature

Magnesium is a primary rock-forming element widely distributed in the mantle, crust, and hydrosphere. In the mantle and crust, Mg mainly occurs in silicate minerals such as olivine, pyroxene, and mica, as well as in carbonates such as dolomite. In the hydrosphere, Mg is the most abundant metallic element in both river water and seawater. Previous studies have documented systematic isotopic differences among these reservoirs (Figure 9).

Mg isotope composition of natural materials. (Data sources: Galy et al., 2002; Chang et al., 2003; Chang et al., 2004; Tipper et al., 2006a; Tipper et al., 2006b; Tipper et al., 2008a; Tipper et al., 2008b; Buhl et al., 2007; Wiechert and Halliday, 2007; Brenot et al., 2008; Pogge von Strandmann et al., 2008; Handler et al., 2009; Yang et al., 2012; Bourdon et al., 2010; Teng et al., 2010a; Ling et al., 2011).

Scatter plot showing δ²⁶Mg values for different materials. Blue diamonds represent foraminifera, purple squares denote stalactite, grey diamonds indicate limestone, yellow squares for dolomite, blue bars for seawater, green bar for river water, and red circles for soil. A shaded gray area indicates mantle peridotite and basalt, centered around zero δ²⁶Mg value. The x-axis ranges from -6 to 1.

The upper mantle exhibits a relatively homogeneous Mg isotopic composition. Chondrites, mantle peridotites, and basalts yield δ26Mg values of −0.49 to +0.06‰, averaging around −0.24‰ (Ke et al., 2011), indicating a narrow mantle range. In contrast, the continental crust, especially the upper crust, shows significant isotopic heterogeneity. Soils are commonly enriched in heavy Mg isotopes relative to basalt, whereas sedimentary carbonates are generally isotopically lighter. Carbonates display large variations (δ26Mg = −4.84 to −1.00‰; Galy et al., 2002; Tipper et al., 2006a; Tipper et al., 2006b; Pogge von Strandmann et al., 2008). Within this group, dolomite tends to have relatively heavier values (−2.29 to −1.09‰), whereas limestones are lighter (−4.47 to −2.43‰). Such differences reflect equilibrium fractionation among dolomite, calcite, and water (Galy et al., 2002). Stalactites and foraminifera yield even lighter δ26Mg values, indicating significant biological or low-temperature fractionation (Chang et al., 2004; Pogge von Strandmann et al., 2008).

In the hydrosphere, seawater is isotopically uniform (δ26Mg ≈ −0.83‰; Chang et al., 2003; Young and Galy, 2004; Ling et al., 2011), consistent with its long residence time (>10 Ma; Li, 1982). By contrast, continental waters are highly variable: river waters range from −2.08 to −0.52‰ (Young and Galy, 2004; Tipper et al., 2006a; Tipper et al., 2006b; Brenot et al., 2008).

The mechanisms controlling these variations are well established. High-temperature magmatic processes produce only minor Mg isotope fractionation (Teng et al., 2007; Teng et al., 2010a; Liu et al., 2010), whereas low-temperature water–rock interaction induces large isotopic shifts (Galy et al., 2002; Chang et al., 2004; Young and Galy, 2004; Pogge von Strandmann et al., 2008). Weathering of silicate rocks preferentially removes light Mg isotopes into solution, leaving soils isotopically heavy and rivers/seawater relatively light (Tipper et al., 2006a; Tipper et al., 2006b; Brenot et al., 2008). During carbonate precipitation, sediments preferentially incorporate light isotopes (Galy et al., 2002; Chang et al., 2004).

Against this global framework, our new data show that the Bayan Obo siderite and dolomite carbonatites have distinctly heavier δ26Mg values than Mesoproterozoic sedimentary dolomites from the same region (Figure 9). This isotopic distinction demonstrates that the Bayan Obo ore-bearing carbonatites cannot be explained by a simple sedimentary origin, and instead show closer affinity to mantle values with only limited modification from wall-rock interaction. These results provide a critical baseline for constraining the genetic models of the deposit.

 5.3 Constraints on genesis of the Bayan Obo deposit

The genesis of the Bayan Obo ore-bearing dolomite carbonatites has been widely debated, with models including: (1) microbial or diagenetic micritic mounds genesis, which emphasizes microbial or early diagenetic processes (Qiao et al., 1997; Zhang P. S. et al., 1998); (2) normal sedimentary genesis, proposing a typical marine carbonate sedimentation origin (Meng, 1982; Meng and Drew, 1992; Wei and Shangguan, 1983); and (3) primary carbonatite genesis, which interprets the dolomite carbonatites as direct products of mantle-derived carbonatitic magmatism (Bai et al., 1996; Le Bas et al., 1997; Mitchell, 2005; Yang and Le Bas, 2004). Mineralogical and geochemical features suggest classification of the Bayan Obo siderite carbonatites as magmatic carbonatites (Mitchell, 2005), and the Mg isotope data presented here provide new constraints on these interpretations.

5.3.1 Possibility of the formation of micrite mound

Micrite mound, also known as carbonate mud mound, is mainly composed of plaster, with only a small number of organisms and biological debris (Fan and Zhang, 1985). Macroscopically, it is generally a discus-shaped body having a convex top and a flat bottom, with the thickness ranging from a few meters to many tens of meters. It appears in a strip parallel to the ancient coastline in areas with deeper water slopes (Qiao et al., 1997). Typical micrite mound in China include: the top micrite mound of the Cambrian system in Xishan, Beijing, and the top micrite mound of the Sailinhudong Group in the Heinaobao area of Inner Mongolia (about 25 km southeast of the Bayan Obo mining area) (Qiao et al., 1997). Qiao et al. (1997) compared the ore-bearing dolomite carbonatite of the Bayan Obo deposit with the micrite mound of the Sailinhudong area, and suggested that the two may belong to the same stratigraphic horizon based on their macroscopic geological features. Subsequently, Zhang Y. X. et al. (1998), Zhang et al. (2008) argued that the micritic carbonates in the Bayan Obo deposit formed as products of seafloor hydrothermal activity (with CO2 involvement), primarily by chemical sedimentation, with REE, Nb, Na, and F sourced from deep hydrothermal fluids, and Ca, Mg, and Fe mainly derived from seawater.

Our Mg isotope results provide new geochemical constraints on this issue. The Bayan Obo siderite carbonatites (δ26Mg = −1.31–−0.09‰) and dolomite carbonatites (δ26Mg = −1.12–−0.31‰) display significantly higher values than Mesoproterozoic sedimentary dolomites (δ26Mg = −2.50–−2.00‰), but are distinguishable from the Sailinhudong micrite mound dolomites (δ26Mg = −1.99–−1.93‰) (Figure 9). These observations indicate that, although the Mg isotopic composition of Bayan Obo carbonatites partly overlaps with sedimentary fields, they are clearly shifted relative to the Sailinhudong micrite mound dolomites. This offset suggests that the Bayan Obo carbonatites cannot be simply interpreted as equivalents of the Sailinhudong micrite mounds, and therefore Mg isotopes highlight the difficulty of establishing a direct genetic relationship between the two.

 5.3.2 Possibility of normal sedimentary genesis

Previous study indicates that the Mg isotope composition of sedimentary dolomites varies significantly (δ26Mg is −2.29–−1.09‰; Galy et al., 2002; Chang et al., 2004; Tipper et al., 2006b; Brenot et al., 2008). This might be related to factors like varying ages or geological backgrounds of dolomite. To compare with the dolomite carbonatites from the Bayan Obo deposit, it is necessary to obtain the Mg isotope composition of the sedimentary dolomites in the H8 rock section of the Bayan Obo Group. This can be obtained from the Kuangou north sedimentary dolomite, whose δ26Mg is −1.67 –−1.54‰. This should represent the Mg isotope composition characteristics of the dolomites in the H8 rock section of the Bayan Obo Group. Moreover, the δ26Mg of the Middle Proterozoic dolomite in Pingquan area is −1.81–−1.53‰, consistent with the Mg isotope of the sedimentary dolomites in Kuangou north. This suggests that the sedimentary dolomites in Bayan Obo Pingquan area during the Middle Proterozoic period has a consistent Mg isotope composition. The δ26Mg of sedimentary dolomites from the Middle Proterozoic in the Bayan Obo area should be −1.81–−1.53‰.

The δ26Mg isotope compositions of the dolomite carbonatites and extracted siderite carbonatite samples from Bayan Obo are significantly heavier than those of Mesoproterozoic sedimentary dolomites, with none of the samples falling within the δ26Mg range of the sedimentary dolomites (Figure 9). This suggests that the ore-bearing dolomite carbonatites are not normally deposited dolomites. The elevated δ26Mg values may result from magmatic differentiation and hydrothermal processes, which preferentially retain heavy Mg isotopes in the residual carbonatites. In addition, fluid–rock interaction, partial recrystallization, or carbonate precipitation from Mg-rich deep fluids could further modify the Mg isotope composition, distinguishing these carbonatites from sedimentary dolomites formed under surface depositional conditions.

 5.3.3 Possibility of magmatic genesis

The Mg isotopic ratios of the carbonatites exhibit a broader range than those of the Mesoproterozoic sedimentary dolomites (Figure 6). This variability likely reflects the superposition of several processes: (1) magmatic differentiation and fractional crystallization, which generate isotopic heterogeneity in the carbonatite system; (2) hydrothermal alteration and fluid–rock interaction, which preferentially mobilize light Mg isotopes and shift isotopic compositions; and (3) post-depositional recrystallization of carbonates, which may further modify primary signatures. Together, these processes account for the wider spread of δ26Mg values in the carbonatites compared with the relatively homogeneous dolomites.

Magnesium isotope fractionation is strongly temperature dependent. At magmatic temperatures, equilibrium fractionation is generally limited (Pogge von Strandmann et al., 2011; Ling et al., 2011), whereas at low temperatures water–rock interaction can produce significant fractionation (Galy et al., 2002; Chang et al., 2004; Young and Galy, 2004; Pogge von Strandmann et al., 2008; Teng et al., 2010a; Teng et al., 2010b). During silicate magmatic differentiation, light Mg isotopes partition preferentially into fluids, enriching residual silicates in heavy isotopes, whereas rivers and seawater are relatively enriched in light isotopes (Tipper et al., 2006a; Tipper et al., 2006b; Tipper et al., 2008a; Tipper et al., 2008b; Brenot et al., 2008). In contrast, carbonate precipitation favors incorporation of light isotopes into the solid phase (Galy et al., 2002; Chang et al., 2004).

In the Bayan Obo siderite carbonatites, some samples yield exceptionally light δ26Mg values (−1.31‰), far lower than mantle compositions. Possible causes include magmatic differentiation, weathering, hydrothermal alteration, wall-rock assimilation, and fluid exsolution during magma ascent (Teng et al., 2011; Chen et al., 2020). Experimental studies show that silicate–carbonate isotope exchange at high temperatures favors retention of heavy Mg isotopes in silicate melts, leaving carbonates isotopically lighter (Schauble, 2011; Macris et al., 2013). Liquid immiscibility and fractional crystallization of carbonatite melts may further produce ∼0.2‰ fractionation (Li et al., 2016), though this alone cannot account for the extremely light δ26Mg values observed.

Crustal contamination is unlikely to be the sole explanation, as crustal rocks generally overlap mantle values (−0.40–+0.12‰; Rudnick and Gao, 2014; Huang et al., 2013; Cheng et al., 2017). Instead, assimilation of Mg-rich carbonate wall rocks is more consistent with the data. Local dolostones exhibit lighter δ26Mg values (−1.81–−1.53‰), and petrographic evidence reveals marble xenoliths within the carbonatites. Thus, wall-rock assimilation, likely enhanced by fluid exsolution during magma ascent, offers a plausible mechanism for the isotopically light signatures in Bayan Obo siderite carbonatites.

Overall, the Mg source of the Bayan Obo carbonatites is dominantly mantle-derived, but the observed isotopic heterogeneity reflects overprinting by fractional crystallization, liquid immiscibility, hydrothermal alteration, and particularly wall-rock assimilation.

Independent isotopic systems corroborate this interpretation. (1) S isotopes: Sulfides (pyrite, galena) display mantle-like sulfur isotope compositions, whereas whole rocks and barites fall between mantle and seafloor sedimentary sulfate endmembers (Ding et al., 2003). (2) Sr–Nd isotopes: The deposit shows clear mantle signatures (Zhang et al., 2003). (3) C–O isotopes: δ13CV-PDB and δ18OV-SMOW values of siderite carbonatites (−3.66–−0.13‰ and 10.89–13.66‰, respectively) resemble dolomite carbonatites and fall within the range between primary igneous carbonatites and sedimentary rocks (Liu, 1986; Cao et al., 1994; Wei et al., 2022). In contrast, sedimentary dolostones and other carbonates typically yield δ13C > −4‰ and δ18O > 18‰ (Tang, 2022). These combined data indicate that Bayan Obo mineralization was driven by mantle-derived carbonatitic magmatism with localized carbonate assimilation during magma degassing.

An alternative model proposes that sedimentary carbonate rocks were directly replaced by mantle-derived carbonatites or fluids. This scenario is inconsistent with current observations. First, the majority of siderite carbonatite δ26Mg values fall within the mantle range, with no samples overlapping sedimentary endmembers, indicating a dominant mantle Mg source. Second, if sedimentary carbonate rocks had been extensively replaced, isotopic compositions would more strongly retain sedimentary characteristics at lower δ26Mg values, reflecting incomplete metasomatism. Instead, Bayan Obo siderite carbonatites are characterized by low SiO2 and CaO, but high TFe2O3, MgO, MnO, and REE (Table 1), with pronounced enrichment in LREE and incompatible elements such as Sr and Ba. These geochemical traits match mantle-derived igneous carbonatites and sharply contrast with sedimentary carbonates, which typically contain <200 ppm Sr + Ba and <25 ppm REE (Kang et al., 2024).

Taken together, isotopic and geochemical evidence demonstrates that the Bayan Obo siderite carbonatites are best explained as products of mantle-derived carbonatitic magmatism that experienced substantial modification by fractional crystallization, immiscibility, hydrothermal alteration, and wall-rock assimilation, rather than by wholesale replacement of sedimentary carbonate rocks.

 6 Conclusion

1. The Bayan Obo siderite carbonatites are characterized by high concentrations of TFe2O3, MgO, and MnO, and low contents of SiO2, K2O, Na2O, and P2O5—geochemical features typical of ferro carbonatites. In terms of trace element composition, these rocks are enriched in Ba, Sr, Th, and REE, exhibiting geochemical signatures similar to those of the ore-bearing dolomite carbonatites, but distinct from sedimentary carbonate rocks. When integrated with previous studies, these characteristics support the interpretation that the Bayan Obo siderite carbonatites may have originated from a mantle source influenced by mixing between HIMU and EM I components.

2. The δ26Mg values of the Bayan Obo siderite carbonatites range from −1.31 to −0.09‰, clustering near the mantle endmember. Their Mg isotope compositions closely resemble those of the ore-bearing dolomite carbonatites, supporting a mantle-derived magmatic origin for both rock types.

3. Combining the results of this study with previous research, we propose that the Bayan Obo deposit formed through carbonatitic magma intrusion and associated metasomatism. The genesis of the deposit was primarily controlled by mantle-derived carbonatitic magmatism, with subsequent magma differentiation and evolution likely driving the enrichment of Fe, REE, Sr, Ba, and Y.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

 Author contributions

CR: Writing – original draft, Conceptualization, Methodology. HJ: Writing – review and editing, Supervision. QZ: Writing – review and editing. HS: Writing – review and editing, Supervision, Formal analysis, Funding acquisition. JL: Writing – review and editing. BL: Writing – review and editing.

 Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The work was supported by the National Science and Technology Major Project (2024ZD1001001), Basic Research Fund of the Chinese Academy of Geological Sciences (KK2108), National Natural Science Foundation of China (42072114), National Key Research and Development Program of China (2022YFC2905301), and Project for the Application and Transformation of Research Outcomes (HE2513).

The authors would like to thank the editor and six reviewers for their helpful and constructive comments that greatly contributed to the improvement of the manuscript. Thanks to Dongsheng Wang from the Institute of Mineral Resources, Chinese Academy of Geological Sciences for his assistance during the field work and the article revision process.

 Conflict of interest

Author QZ was employed by PetroChina Xinjiang Oilfield Company.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer WJ declared a shared affiliation with the authors CR, HJ, HS, JL, BL to the handling editor at time of review.

 Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Bai G. Yuan Z. X. (1985). Geology of carbonatite and related ore deposits. Beijing: Bulletin of Institute of Mineral Deposit. Bai G. Wu C. Y. Yuan Z. X. Zhang Z. Q. Zheng L. (1996). Demonstration on the geological features and genesis of the Bayan Obo ore deposit. Beijing: Geological Publishing House. Belshaw N. Zhu X. Guo Y. O'Nions R. (2000). High precision measurement of iron isotopes by plasma source mass spectrometry. Int. J. Mass Spectrom. 197, 191–195. 10.1016/s1387-3806(99)00245-6 Bourdon B. Tipper E. T. Fitoussi C. Stracke A. (2010). Chondritic mg isotope composition of the earth. Geochimica Cosmochimica Acta 74 (17), 5069–5083. 10.1016/j.gca.2010.06.008 Brenot A. Cloquet C. Vigier N. Carignan J. France-Lanord C. (2008). Magnesium isotope systematics of the lithologically varied Moselle River basin, France. Geochimica Cosmochimica Acta 72, 5070–5089. 10.1016/j.gca.2008.07.027 Buhl D. Immenhauser A. Smeulders G. Kabiri L. Richter D. K. (2007). Time series δ26Mg analysis in speleothem calcite: kinetic versus equilibrium fractionation, comparison with other proxies and implications for palaeoclimate research. Chem. Geol. 244 (3-4), 715–729. 10.1016/j.chemgeo.2007.07.019 Cao R. L. Zhu S. H. Wang J. W. (1994). Material source and genesis of Bayan Obo Fe-REE deposit. Chin. Sci. Geoscience 24 (12), 1298–1307. Chang V. T. C. Makishima A. Belshaw N. S. O'Nions R. K. (2003). Purification of Mg from low-Mg biogenic carbonates for isotope ratio determination using multiple collector ICP-MS. J. Anal. Atomic Spectrom. 18 (4), 296–301. 10.1039/b210977h Chang V. T. C. Williams R. Makishima A. Belshawl N. S. O’Nions R. K. (2004). Mg and Ca isotope fractionation during CaCO3 biomineralisation. Biochem. Biophysical Res. Commun. 323, 79–85. 10.1016/j.bbrc.2004.08.053 15351704 Chao E. C. T. Back J. M. Minkin J. A. Tatsumoto M. Wang J. W. Conrad J. E. (1997). The sedimentary carbonate-hosted giant Bayan Obo REE-Fe-Nb ore deposit of inner Mongolia, China: a cornerstone example for giant polymetallic ore deposits of hydrothermal origin. Cent. Integr. Data Anal. Wis. Sci. Cent. Chen W. Liu H. Y. Lu J. Jiang S. Y. Simonetti A. Xu C. (2020). The formation of the ore-bearing dolomite marble from the giant Bayan Obo REE-Nb-Fe deposit, Inner Mongolia: insights from micron-scale geochemical data. Miner. Deposita 55, 131–146. 10.1007/s00126-019-00886-4 Chen J. Gong Y. L. Chen L. Xiang M. Tian S. H. (2021). New advances in magnesium isotope geochemistry and its application to carbonatite rocks. Earth Sci. 46 (12), 4366–4389. 10.3799/dqkx.2021.140 Cheng Z. G. Zhang Z. C. Hou T. Santosh M. Chen L. Ke S. (2017). Decoupling of Mg-C and Sr-Nd-O isotopes traces the role of recycled carbon in magnesiocarbonatites from the tarim large igneous Province. Geochimica Cosmochimica Acta 202, 159–178. 10.1016/j.gca.2016.12.036 Ding T. P. Jiang S. Y. Tian S. H. (2003). The genesis of the ore hosting “Dolomitic Marble” in the Bayun Obo deposit, Inner Mongolia, China: constrained by isotopic results. Acta Geosci. Sin. 24 (6), 535–542. Elburg M. A. Van Bergen M. Hoogewerff J. Foden J. Vroon P. Zulkarnain I. (2022). Geochemical trends across an arc-continent collision zone: magma sources and slab-wedge transfer processes below the Pantar Strait volcanoes, Indonesia. Geochimica Cosmochimica Acta 66 (15), 2771–2789. 10.1016/s0016-7037(02)00868-2 Elliott H. Wall F. Chakhmouradian A. R. Siegfried P. R. Dahlgren S. Weatherley S. (2018). Fenites associated with carbonatite complexes: a review. Ore Geol. Rev. 93, 38–59. 10.1016/j.oregeorev.2017.12.003 Fan J. S. Zhang W. (1985). On the basic concept and classification of organic reefs and their main identifying criteria. Acta Petrol. Sin. 3, 45–59. Fei X. J. Zeng P. S. Yan C. J. Wen L. G. (2019). The regional geological evolution history of the Bayan Obo REE-Nb-Fe deposit derived from isotopic ages. Geol. Explor. 55 (2), 04610471. Foley S. F. Barth M. G. Jenner G. A. (2000). Rutile/melt partition coefficients for trace elements and an assessment of the influence of rutile on the trace element characteristics of subduction zone magmas. Geochimica Cosmochimica Acta 64 (5), 933–938. 10.1016/s0016-7037(99)00355-5 Gaetani G. A. Asimow P. D. Stolper E. M. (2008). A model for rutile saturation in silicate melts with applications to eclogite partial melting in subduction zones and mantle plumes. Earth Planet. Sci. Lett. 272 (3-4), 720–729. 10.1016/j.epsl.2008.06.002 Galy A. Belshaw N. S. Halicz L. O'Nions R. K. (2001). High-precision measurement of magnesium isotopes by multiple-collector inductively coupled plasma mass spectrometry. Int. J. Mass Spectrom. 208, 89–98. 10.1016/s1387-3806(01)00380-3 Galy A. Bar-Matthews M. Halicz L. O'Nions R. K. (2002). Mg isotopic composition of carbonate: insight from speleothem formation. Earth Planet. Sci. Lett. 201, 105–115. 10.1016/s0012-821x(02)00675-1 Gao Z. Y. (2022). PhD dissertation. Study on Mesoproterozoic and late Paleozoic tectonic setting, magmatic source characteristics and rare earth metallogenic mechanism in Bayan Obo. Xian: Chang’an University. Ge L. Jiang S. Y. (2008). Recent advances in research on magnesium isotope geochemistry. Acta Petrologica Mineralogica 27 (4), 367–374. Gothmanna M. Stolarski J. Adkinsj F. Higgins J. A. (2017). A Cenozoic record of seawater Mg isotopes in well-preserved fossil corals. Geology 45 (11), 1039–1042. 10.1130/g39418.1 Hammouda T. Moine B. Devidal J. L. Vincent C. (2009). Trace element partitioning during partial melting of carbonated eclogites. Phys. Earth Planet. Interiors 174 (1-4), 60–69. 10.1016/j.pepi.2008.06.009 Handler M. R. Baker J. A. Schiller M. Bennett V. C. Yaxley G. M. (2009). Magnesium stable isotope composition of Earth’s upper mantle. Earth Planet. Sci. Lett. 282, 306–313. 10.1016/j.epsl.2009.03.031 Hao Z. G. Wang X. B. Li Z. Xiao G. W. Zhang T. R. (2002). Petrological study of alkaline basic dyke and carbonatite dyke in Bayan Obo, Inner Mongolia. Acta Petrologica Mineralogica 21 (4), 429–444. He X. X. Li S. Z. Tang S. H. (2008). Advances in the study of Mg isotopes application. Acta Petrologica Mineralogica 27 (5), 472–476. Hou Z. L. (1989). The Bayan Obo Fe-Nb-REE-deposit: its basic geological features, metallogenesis and genetic model. Geol. Prospect. 7, 1–5+22. Hou Z. Q. Tian S. H. Yuan Z. X. Xie Y. L. Yin S. P. Yi L. S. (2006). The Himalayan collision zone carbonatites in Western Sichuan, SW China: petrogenesis, mantle source and tectonic implication. Earth Planet. Sci. Lett. 244 (1-2), 234–250. 10.1016/j.epsl.2006.01.052 Hou Z. Q. Liu Y. Tian S. H. Yang Z. Xie Y. (2015). Formation of carbonatite-related giant rare-earth-element deposits by the recycling of marine sediments. Sci. Rep. 5 (1), 10231–10. 10.1038/srep10231 26035414 Hu L. Li Y. K. Sun S. Li R. P. Ke C. H. Wang A. J. (2023). Identification of new igneous carbonatites in the Bayan Obo area, Inner Mongolia. Geol. China 50 (6), 1788–1803. Huang K. J. Teng F. Z. Amira E. Li W. Y. Bao Z. Y. (2013). Magnesium isotopic variations in loess: origins and implications. Earth Planet. Sci. Lett. 374, 60–70. 10.1016/j.epsl.2013.05.010 Jones A. P. Genge M. Carmody L. (2013). Carbonate melts and carbonatites. Rev. Mineralogy Geochem. 75, 289–322. 10.2138/rmg.2013.75.10 Kang Q. Q. Chen Z. L. Pan J. Y. Li P. Liu Y. L. Li L. (2024). Spatial and temporal distribution, geochemical characteristics of carbonatites and their relationship with U-REE mineralization in the Xiaoqinling area, Shaanxi Province. J. Geomechanics 30 (1), 147–167. Ke S. Liu S. A. Li W. Y. (2011). Advances and application in magnesium isotope geochemistry. Acta Petrol. Sin. 27 (2), 383–397. Lai X. D. (2013). Dissertation on Ph.D. study on the genesis of Bayan Obo REE-Fe ore deposit, Inner Mongolia. Hefei: University of Science and Technology of China. Le Bas M. J. (1987). Nephelinites and carbonatites. Geol. Soc. 30 (1), 53–83. 10.1144/gsl.sp.1987.030.01.05 Le Bas M. J. Spiro B. Yang X. M. (1997). Oxygen, carbon and strontium isotope study of the carbonatitic complex at Bayan Obo Fe-Nb-REE deposit, Inner Mongolia, N China. Mineral. Mag. 61, 531–541. 10.1180/minmag.1997.061.407.05 Le Maitre R. W. (2002). Igneous rocks: a classification and glossary of terms. Cambridge: Cambridge University Press, 1–236. Lentz D. R. (1999). Carbonatite genesis: a reexamination of the role of intrusion-related pneumatolytic skarn processes in limestone melting. Geology 27, 335–338. 10.1130/0091-7613(1999)027<0335:cgarot>2.3.co;2 Lentz D. R. (2014). Reexamination of the genesis of the Bayan Obo Fe-REE-Nb deposit: autometasomatic oxidation of an extremely fractionated ferrocarbonatite. Acta Geol. Sinica-English Ed. 88 (s2), 361–363. 10.1111/1755-6724.12372_7 Li Y. H. (1982). A brief discussion on the mean oceanic residence time of elements. Earth Planet. Sci. Lett. 46, 2671–2675. 10.1016/0016-7037(82)90386-6 Li W. Y. Teng F. Z. Halama R. Keller J. Klaudius J. (2016). Magnesium isotope fractionation during carbonatite magmatism at Oldoinyo Lengai, Tanzania. Earth Planet. Sci. Lett. 444, 26–33. 10.1016/j.epsl.2016.03.034 Li Y. K. Ke C. H. She H. Q. (2023). Geology and mineralization of the Bayan Obo supergiant carbonatite-type REE-Nb-Fe deposit in Inner Mongolia, China: a review. China Geol. 6, 716–750. 10.31035/cg2023082 Li X. C. Fan H. R. Su J. H. Groves D. I. Yang K. F. Zhao X. F. (2024). Giant rare earth element accumulation related to voluminous, highly evolved carbonatite: a microanalytical study of carbonate minerals from the Bayan Obo deposit, China. Econ. Geol. 119 (2), 373–393. 10.5382/econgeo.5060 Ling M. X. Sedaghatpour F. Teng F. Z. Hays P. D. Strauss J. Sun W. (2011). Homogeneous magnesium isotopic composition of seawater: an excellent geostandard for Mg isotope analysis. Rapid Commun. Mass Spectrom. 25 (19), 2828–2836. 10.1002/rcm.5172 21913261 Ling M. X. Liu Y. L. Williams I. S. Teng F. Z. Yang X. Y. Ding X. (2013). Formation of the world’s largest REE deposit through protracted fluxing of carbonatite by subduction-derived fluids. Nat. Sci. Rep. 3 (3), 1776. 10.1038/srep01776 Ling M. X. Zhang H. Li H. Liu Y. L. Liu J. Li L. Q. (2014). The Permian-Triassic granitoids in Bayan Obo, North China Craton: a geochemical and geochronological study. Lithos 190, 430–439. 10.1016/j.lithos.2014.01.002 Liu T. G. (1986). A discussion on the genesis of dolomite in Bayan Obo, inner Mongolia—with emphasis on the composition of oxygen and carbon isotopes. Geol. Rev. 32 (2), 150–159. Liu S. A. Teng F. Z. He Y. Ke S. Li S. (2010). Investigation of magnesium isotope fractionation during granite differentiation: implication for Mg isotopic composition of the continental crust. Earth Planet. Sci. Lett. 297, 646–654. 10.1016/j.epsl.2010.07.019 Liu S. Fan H. R. Yang K. F. Hu F. F. Rusk B. Liu X. (2018). Fenitization in the giant Bayan Obo REE-Nb-Fe deposit: implication for REE mineralization. Ore Geol. Rev. 94, 290–309. 10.1016/j.oregeorev.2018.02.006 Liu J. W. Tian S. H. Wang L. (2023). Application of magnesium stable isotopes for studying important geological processes-a review. Earth Sci. Front. 30 (3), 399–424. Macris C. A. Young E. D. Manning C. E. (2013). Experimental determination of equilibrium magnesium isotope fractionation between spinel, forsterite, and magnesite from 600 to 800 °C. Geochimica Cosmochimica Acta 118, 18–32. 10.1016/j.gca.2013.05.008 Meng Q. R. (1982). The genesis of the host rock-dolomite of the Bayan Obo iron ore deposits and the analysis of its sedimentary environment. Geol. Rev. (05), 481–489+508. Meng Q. R. Drew L. J. (1992). Study on oxygen and carbon isotope and the implication for genesis of Bayan Obo ore-bearing H8 dolomite. Contributions Geol. Mineral Resour. Res. 7 (2), 46–54. Mitchell R. H. (2005). Carbonatites and carbonatites and carbonatites. Can. Mineralogist 43 (6), 2049–2068. 10.2113/gscanmin.43.6.2049 Pogge von Strandmann P. A. E. Burton K. W. James R. H. van Calsteren P. Gislason S. R. Sigmarsson O. (2008). The influence of weathering processes on riverine magnesium isotopes in a basaltic terrain. Earth Planet. Sci. Lett. 276, 187–197. 10.1016/j.epsl.2008.09.020 Pogge von Strandmann P. A. E. Elliott T. Marschall H. R. Coath C. Lai Y. J. Jeffcoate A. B. (2011). Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths. Geochimica Cosmochimica Acta 75, 5247–5268. 10.1016/j.gca.2011.06.026 Qiao X. F. Gao L. Z. Peng Y. Zhang Y. X. (1997). Composite stratigraphy of the Sailinhuodong Group and ore-bearing micrite mound in the Bayan Obo deposit, Inner Mongolia, China. Acta Geol. Sin. 71 (3), 202–211. Riechelmann S. Mavromatis V. Buhl D. Dietzel M. Eisenhauer A. Immenhauser A. (2016). Impact of diagenetic alteration on brachiopod shell magnesium isotope (δ26Mg) signatures: experimental versus field data. Chem. Geol. 440, 191–206. 10.1016/j.chemgeo.2016.07.020 Rudnick R. L. Gao S. (2014). “Composition of the continental crust,” in Treatise on geochemistry (Amsterdam: Elsevier), 4, 1–64. 10.1016/b0-08-043751-6/03016-4 Schauble E. A. F. (2011). First-principles estimates of equilibrium magnesium isotope fractionation in silicate, oxide, carbonate and hexaaquamagnesium (2+) crystals. Geochimica Cosmochimica Acta. 75 (3), 844–869. 10.1016/j.gca.2010.09.044 She H. D. Fan H. R. Santosh M. Li X. C. Yang K. F. Wang Q. W. (2023). Paleozoic remelting of carbonatite in Bayan Obo (China): further insights into the formation of a giant REE deposit. Gondwana Res. 119, 172–185. 10.1016/j.gr.2023.03.018 Simonetti A. Bell K. Viladkar S. G. (1995). Isotopic data from the amba dongar carbonatite complex, west-central India: evidence for an enriched mantle source. Chem. Geol. 122 (1-4), 185–198. 10.1016/0009-2541(95)00004-6 Smith M. P. Campbell L. S. Kynicky J. (2015). A review of the genesis of the world class Bayan Obo Fe-REE-Nb deposits, Inner Mongolia, China: multistage processes and outstanding questions. Ore Geol. Rev. 64, 459–476. 10.1016/j.oregeorev.2014.03.007 Sun J. (2013). “The origin of the Bayan Obo deposit, Inner Mongolia, China: the iron and magnesium isotope constraints,”. PhD dissertation (Beijing: China University of Geosciences). Sun J. Fang N. Li S. Z. Chen Y. L. Zhu X. K. (2012). Magnesium isotopic constraints on the genesis of Bayan Obo ore deposit. Acta Petrol. Sin. 28 (9), 2890–2902. Tang H. Y. (2022). “Genesis and resource potential analysis of the Bayan Obo Fe-REE deposit,”. PhD dissertation (Jilin: Jilin University). Tang Q. Y. Bao J. Dang Y. X. Ke S. Zhao Y. (2018). Mg-Sr-Nd isotopic constraints on the genesis of the giant Jinchuan Ni-Cu-(PGE) sulfide deposit, NW China. Earth Planet. Sci. Lett. 502, 221–230. 10.1016/j.epsl.2018.09.008 Teng F. Z. (2017). Magnesium isotope geochemistry. Rev. Mineralogy Geochem. 82 (1), 219–287. 10.2138/rmg.2017.82.7 Teng F. Z. Wadhwa M. Helz R. T. (2007). Investigation of magnesium isotope fractionation during basalt differentiation: implications for a chondritic composition of the terrestrial mantle. Earth Planet. Sci. Lett. 261, 84–92. 10.1016/j.epsl.2007.06.004 Teng F. Z. Li W. Y. Ke S. Marty B. Dauphas N. Huang S. (2010a). Magnesium isotopic composition of the Earth and chondrites. Geochimica Cosmochimica Acta 74, 4150–4166. 10.1016/j.gca.2010.04.019 Teng F. Z. Li W. Y. Rudnick R. L. Gardner L. R. (2010b). Contrasting lithium and magnesium isotope fractionation during continental weathering. Earth Planet. Sci. Lett. 300, 63–71. 10.1016/j.epsl.2010.09.036 Teng F. Z. Dauphas N. Helz R. T. Gao S. Huang S. (2011). Diffusion-driven magnesium and iron isotope fractionation in Hawaiian olivine. Earth Planet. Sci. Lett. 308, 317–324. 10.1016/j.epsl.2011.06.003 Tilton G. R. Bryce J. G. Mateen A. (1998). Pb–Sr–Nd isotope data from 30 and 300 Ma collision zone carbonatites in Northwest Pakistan. J. Petrology, 11–12. 10.1093/petroj/39.11-12.1865 Tipper E. T. Galy A. Bickle M. J. (2006a). Riverine evidence for a fractionated reservoir of Ca and Mg on the continents: implications for the oceanic Ca cycle. Earth Planet. Sci. Lett. 247, 267–279. 10.1016/j.epsl.2006.04.033 Tipper E. T. Galy A. Gaillardet J. Bickle M. J. Elderfield H. Carder E. A. (2006b). The magnesium isotope budget of the modern ocean: constraints from riverine magnesium isotope ratios. Earth Planet. Sci. Lett. 250, 241–253. 10.1016/j.epsl.2006.07.037 Tipper E. T. Galy A. Bickle M. J. (2008a). Calcium and magnesium isotope systematics in rivers draining the Himalaya-Tibetan-Plateau region: lithological or fractionation control? Geochimica Cosmochimica Acta 72, 1057–1075. 10.1016/j.gca.2007.11.029 Tipper E. T. Louvat P. Capmas F. Galy A. Gaillardet J. (2008b). Accuracy of stable Mg and Ca isotope data obtained by MC-ICP-MS using the standard addition method. Chem. Geol. 257, 65–75. 10.1016/j.chemgeo.2008.08.016 Turner S. Hawkesworth C. Rogers N. Bartlett J. Worthington T. Hergt J. (1997). 238U-230Th disequilibria, magma petrogenesis, and flux rates beneath the depleted Tonga-Kermadec island arc. Geochimica Cosmochimica Acta 61 (22), 4855–4884. 10.1016/s0016-7037(97)00281-0 Vladykin N. V. Pirajno F. (2021). Types of carbonatites: geochemistry, genesis and mantle sources. Lithos 386-387, 105982. 10.1016/j.lithos.2021.105982 Wang X. B. Hao Z. G. Li Z. Xiao G. W. Zhang T. R. (2002). A typical alkaline rock-carbonatite complex in Bayan Obo, Inner Mongolia. Acta Geol. Sin. (04), 501–582. Wang K. Y. Fan H. R. Yang K. F. Hu F. F. Wu C. M. Hu F. Y. (2010). Calcite-dolomite geothermometry of Bayan Obo carbonatites. Acta Petrol. Sin. 26 (04), 1141–1149. Wang S. J. Teng F. Z. Scott J. M. (2016). Tracing the origin of continental HIMU-like intraplate volcanism using magnesium isotope systematics. Geochimica Cosmochimica Acta 185, 78–87. 10.1016/j.gca.2016.01.007 Wang K. Y. Zhang J. E. Fang A. M. Dong C. Hu F. Y. (2018). Genesis of the Bayan Obo deposit, Inner Mongolia: the fenitized mineralization in the ore bodies and its relation to the ore-bearing dolomite. Acta Petrol. Sin. 34 (03), 785–798. Wei J. Y. Shangguan Z. G. (1983). Oxygen isotope composition of magnetite and hematite in Bayan Ebo iron deposit, Inner Mongolia. Chin. J. Geol. (03), 217–224. Wei C. W. Deng M. Xu C. Chakhmouradian A. R. Smith M. P. Kynicky J. (2022). Mineralization of the Bayan Obo rare earth element deposit by recrystallization and decarbonation. Econ. Geol. 117 (6), 1327–1338. 10.5382/econgeo.4926 Whitney D. L. Evans B. W. (2010). Abbreviations for names of rock-forming minerals. Am. mineralogist 95 (1), 185–187. 10.2138/am.2010.3371 Wiechert U. Halliday A. N. (2007). Non-chondritic magnesium and the origins of the inner terrestrial planets. Earth Planet. Sci. Lett. 256, 360–371. 10.1016/j.epsl.2007.01.007 Woolley A. R. Kempe D. R. C. (1989). “Carbonatites: nomenclature, average chemical compositions, and element distribution,” in Carbonatites: genesis and evolution. Editor Bell K. (London, UK: Unwin Hyman), 1–14. Xiao R. G. Liu J. D. Fei H. C. Wang C. Z. Huang X. L. (2006). Types and genesis of Na-Rich rocks in the Bayan Obo REE-Nb-Fe deposit, Inner Mongolia, China. Geoscience 01, 151–164. Xiao R. G. Fei H. C. Wang A. J. Yan K. (2012). A research review of the formation mechanism of the ore-bearing rocks in the Bayan Obo REE-Nb-Fe deposit, Inner Mongolia. Acta Geol. Sin. 86 (05), 735–752. Xie Y. L. Qu Y. W. Yang Z. F. Liang P. Zhong R. C. Wang Q. W. (2019). Giant Bayan Obo Fe-Nb-REE deposit: progresses, controversaries and new understandings. Mineral. deposits 38 (05), 983–1003. Yang X. M. Le Bas M. J. (2004). Chemical compositions of carbonate minerals from Bayan Obo, Inner Mongolia, China: implications for petrogenesis. Lithos 72 (1-2), 97–116. 10.1016/j.lithos.2003.09.002 Yang X. Y. Sun W. D. Zhang Y. X. Zheng Y. F. (2009). Geochemical constraints on the genesis of the Bayan Obo Fe-Nb-REE deposit in Inner Mongolia, China. Earth Planet. Sci. Lett. 73, 1417–1435. 10.1016/j.gca.2008.12.003 Yang W. Teng F. Z. Zhang H. F. Li S. G. (2012). Magnesium isotopic systematics of continental basalts from the North China Craton: implications for tracing subducted carbonate in the mantle. Chem. Geol. 328, 185–194. 10.1016/j.chemgeo.2012.05.018 Yang K. Fan H. Pirajno F. Li X. (2019). The Bayan Obo (China) giant REE accumulation conundrum elucidated by intense magmatic differentiation of carbonatite. Geology 47, 1198–1202. 10.1130/g46674.1 Young E. D. Galy A. (2004). The isotope geochemistry and cosmochemistry of magnesium. Rev. Mineral. Geochem. 55, 197–230. 10.2138/gsrmg.55.1.197 Yuan Z. X. Bai G. Wu C. Y. Zhang Z. Q. Ye X. J. (1991). Metalloge epoch and genesis of the Bayan Obo Niobium-iron deposit, Inner Mongolia. Mineral. Deposits 10 (1), 59–60. Yuan Z. X. Bai G. Wu C. Y. Ding X. S. Zhang Z. Q. (1995). Petrological features of volcanic rocks in H9 formation of the Bayan Obo ore district, Inner Mongolia, and their significance. Mineral. Deposits 14 (3), 197–205. Zhai M. G. Hu B. Zhao T. P. Peng P. Meng Q. R. (2015). Late Paleoproterozoic-Neoproterozoic multi-rifting events in the North China Craton and their geological significance: a study advance and review. Tectonophysics 662, 153–166. 10.1016/j.tecto.2015.01.019 Zhang P. S. Tao K. J. (1996). Mineralogy of Bayan Obo. Beijing: Geology Press, 1–182. Zhang Z. Q. Yuan Z. X. Tang S. H. Wang J. H. Bai G. (1994). New data for ore-forming age of the Bayan Obo REE deposit, Inner Mongolia. Acta Geosci. Sin. 29 (1-2), 85–94. Zhang P. S. Tao K. J. Yang Z. M. (1998a). Rare Earth mineralogy in China. Beijing: Geology Press, 1–234. Zhang Y. X. Peng Y. Qiao X. F. Gao L. Z. Yang X. Y. (1998b). New understanding of the genesis of ore bearing dolomite in the Bayan Obo deposit. Geol. Rev. (01), 70–76. Zhang P. S. Tao K. J. Yang Z. M. Yang X. M. Song R. K. (2001a). Genesis of rare earths, niobium and tantalum minerals in Bayan Obo ore deposit of China. J. Chin. Rare Earth Soc. (02), 97–102. Zhang Z. Q. Tang S. H. Yuan Z. X. Bai G. Wang J. H. (2001b). The Sm-Nd and Rb-Sr isotopic systems of dolomites in the Bayan Obo ore deposit, Inner Mongolia, China. Acta Petrol. Sin. 17 (4), 637–642. Zhang Z. Q. Tang S. H. Wang J. H. Yuan Z. X. Bai G. (2003). Information about ore deposit formation in different epochs: age of the west orebodies of the Bayan Obo deposit with a discussion. Geol. China 30 (2), 130–137. Zhang Y. X. Jiang S. Q. Zhang Q. L. Lai X. D. Peng Y. Yang X. Y. (2008). A discussion on forming time of the Bayan Obo group and ore-forming time of the Bayan Obo giant REE-Nb-Fe deposit, Inner Mongolia. Geol. China. 35 (06), 1129–1137. Zhao G. C. Wilde S. A. Cawood P. A. Li L. Z. (1999). Tectonothermal history of the basement rocks in the Western zone of the North China Craton and its tectonic implications. Tectonophysics 310 (310), 37–53. 10.1016/s0040-1951(99)00152-3 Zhao G. C. Sun M. Wilde S. A. Li S. Z. (2005). Late Archean to Paleoproterozoic evolution of the North China Craton: key issues revisited. Precambrian Res. 136, 177–202. 10.1016/j.precamres.2004.10.002 Zhou Z. L. Li G. Y. Song T. Y. Liu Y. G. (1980). On the geological characteristics and genesis of the dolomites at Bayan Obo, Nei Mongal (Inner Mongolia). Geol. Rev. 26 (1), 35–42.