Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Earth Sci., 08 January 2026

Sec. Geochemistry

Volume 13 - 2025 | https://doi.org/10.3389/feart.2025.1709937

Metallogenic age and hydrothermal evolution of the Huangyangshan graphite deposit in eastern Xinjiang, NW China: constraint from zircon geochronology and geochemistry

Xinhao SunXinhao Sun1Yunsheng RenYunsheng Ren2Zhenjun SunZhenjun Sun2Jingmou LiJingmou Li3Mengjia Huang
Mengjia Huang4*
  • 1Chengdu Center, China Geological Survey (Geosciences Innovation Center of Southwest China), Chengdu, China
  • 2School of Earth Science, Institute of Disaster Prevention, Langfang, China
  • 3Development and Research Center, China Geological Survey, Beijing, China
  • 4Sichuan Institute of Land Science and Technology (Sichuan Center of Satellite Application Technology), Chengdu, China

The Huangyangshan graphite deposit in the East Junggar region of Xinjiang Province is the world’s first documented graphite deposit hosted within an alkaline granite intrusion. Classified as fluid-deposited graphite deposit, it exhibits features of both magmatic and hydrothermal deposits, thus holding special significance for theoretical research. In this study, we conducted LA-ICP-MS zircon U-Pb dating and trace element analyses on zircons from the Huangyangshan deposit orebodies. And comprehensive analysis has been conducted on the REE geochemical signatures of the youngest fine-grained biotite granite within the Huangyangshan pluton. Based on their textures and element compositions, zircons can be divided into two types. Type I zircons are pale, large-grained, with clear oscillatory zonation, low Th, U, and REE contents, and significant positive Ce and strong negative Eu anomalies, which can be classified as magmatic zircons. Type II zircons exhibit dark coloration, fine grain sizes, poorly developed oscillatory zonation, and elevated contents of Th, U, and REEs. They contain few inclusions and have weak LREE and HREE fractionation, consistent with typical hydrothermal zircon. Zircon U-Pb dating results yield an age of 318.8 ± 4.0 Ma for the magmatic zircons, representing the magmatic crystallization age of the graphite-bearing alkaline granite. Hydrothermal zircons have two age groups: One is 299.9 ± 4.6 Ma, which is close to the age of the fine-grained biotite granite, the latest lithoface associated with graphite mineralization. It might be derived from magmatic hydrothermal fluids, representing the graphite mineralization age. The other is 259.9 ± 3.3Ma, which is obviously later than the age of the Huangyangshan pluton. It possibly related to late regional metamorphic fluids. The fine-grained biotite granite has significantly lower total REEs than other lithofacies of the Huangyangshan pluton, with an M-type tetrad effect in REE distribution. The M-type tetrad effect is the result of intense fluid activity, possibly directly related to hydrothermal graphite mineralization. Therefore, the ∼300 Ma hydrothermal event reflected by the first-stage hydrothermal zircons has a close genetic relationship with the fine-grained biotite granite. The Huangyangshan graphite mineralization has experienced a magmatic hydrothermal process and should be classified as a hydrothermal deposit.

1 Introduction

Due to its superior physical and chemical properties—such as insulation, high-temperature resistance and electrical conductivity—graphite is widely used in both traditional industries and high-tech fields. Based on crystallization degree, natural graphite can be categorized into crystalline graphite and cryptocrystalline graphite. The crystalline graphite, characterized by large flake particle size, good washability and high industrial utilization value, has been an important strategic mineral resource in China (Yan et al., 2018). Globally, the vast majority of crystalline graphite deposits are of regional metamorphic or contact metamorphic type (Wada et al., 1995; Papineau et al., 2010; Manoel and Leite, 2018), while a small number of graphite deposits formed by precipitation from carbon-bearing hydrothermal fluids or magmas are collectively referred to as fluid-deposited type (Luque et al., 1998). Compared with metamorphic graphite deposits, fluid-deposited ones are characterized by broad range of mineralization ages, high graphite grade and crystallization degree, and wide range of carbon sources (Duke and Rumble, 1986; Luque et al., 1998; Luque et al., 2012; Luque et al., 2014; Crespo et al., 2006; Doroshkevich et al., 2007; Huizenga, 2011). However, current research on their formation mechanisms remains limited, and the understanding level is relatively low, thus it holds extremely high research value.

Since 2016, three crystalline graphite mineralization subzones have been discovered in the East Junggar of Xinjiang Province, establishing East Junggar as an important potential graphite metallogenic area in China. The crystalline graphite deposit discovered in the Huangyangshan pluton within Qitai County is the first super-large crystalline graphite deposit related to magmatic rocks in China (Zhang et al., 2017). Preliminary exploration estimates that the graphite resources exceed 70 million tons, showing great mineralization potential. Furthermore, the Huangyangshan graphite deposit also has unique characteristics such as “occurring in alkaline granite”, “the ore having unique spherical structure”, and “graphite closely associated with metal sulfides”. Previous studies have on the Huangyangshan graphite deposit have yielded insights into mineral assemblages and geochemical characteristics, diagenetic age and tectonic setting, as well as magmatic evolution processes (Bai et al., 2018; Shao, 2018; Li et al., 2019; Ai et al., 2020; Sun et al., 2021; Sun et al., 2022; Ren et al., 2022). However, due to the lack of research on hydrothermal indicator minerals, there is still a lack of attention to its magmatic-hydrothermal evolution process.

Zircon is a common accessory mineral in granitic rocks and is generally one of the first minerals to crystallize (Hoskin, 2005; Yang et al., 2014; Zeng et al., 2017). It can record the magmatic-hydrothermal evolution process and reflect the evolution of melt-fluid compositions (Belousova et al., 2006; Geisler et al., 2007; Xie et al., 2018; Yin et al., 2023). Hydrothermal zircons that crystallize directly from saturated hydrothermal fluids form across a temperature range spanning the high-temperature magmatic stage to the late medium-temperature hydrothermal stage (≈600∼300 °C; Schaltegger, 2007). These zircons are known to survive multistage hydrothermal events, even when host rocks undergo intensive modification (Zhao et al., 2017; Su et al., 2021; Gao et al., 2023). Therefore, hydrothermal zircons can be used to reliably date primary hydrothermal metasomatism or water-rock interaction events, helping to reconstruct the hydrothermal evolution process. This study focuses on zircons from the Huangyangshan graphite orebody, conducting chronological and trace element studies. Combined with previous research results, we discussed zircon genesis and age, as well as the magmatic-hydrothermal evolution process—aiming to provide a theoretical basis for the study of the mineralization process of fluid-deposited graphite deposits.

2 Regional geology

The East Junggar orogenic belt in Xinjiang Province is located in the middle segment of the Central Asian Orogenic Belt. Geographically, it is adjacent to the Junggar Basin to the west, Mongolia to the east, and the Altai and Tianshan orogenic belts to the north and south respectively. Tectonically, it comprises five units from north to south: Durat-Baytag arc, the Aermantai ophiolite belt, the Yemaquan arc, the Kalamaili ophiolite belt, and the Harlik-Dananhu arc. Geochronological data in this region indicate that the Junggar Ocean and the Altai block formed at least in the Early Ordovician and gradually closed in the Late Paleozoic (Xiao et al., 2009; Xiao et al., 2015; Zhang et al., 2013; Xu et al., 2020). The Karamaili ophiolite belt, situated between the Yemaquan arc and the Harlik–Dananhu arc, is considered to be the site where the Junggar Ocean finally disappeared. It extends in NWW direction along the Karamaili Fault (Figure 1; Liu et al., 2017).

Figure 1
Geological map of a region near Mongolia, depicting various deposits and faults. Areas are color-coded to show different volcanic-sediments and strata. Symbols indicate locations of graphite, gold (Au), copper (Cu), and tin (Sn) deposits. Notable faults include the Suijquan, Kalamailli, Kubusu, and Beitashanqian. Specific locations like Laoyaquan, Beilekuduke, and Huangyangshan are labeled. A legend and scale are included for reference.

Figure 1. Geological sketch map of the Kalamaili area, Eastern Junggar in Xinjiang Province (modified after Sun et al., 2021).

The Kalamaili area, hosting the Huangyangshan graphite deposit, is located south of the Yemaquan arc and north of the Karamaili Fault. During the Carboniferous to Permian, this area experienced a post-collisional extensional tectonic setting, leading to the emplacement of numerous post-collisional alkaline A-type granites such as Laoyaquan, Beilekuduke, and Huangyangshan plutons, which are collectively termed the Kalamaili Granite Belt. These granites all occur in the form of batholiths or stocks, including granodiorite, biotite alkaline granite, hornblende alkaline granite, and riebeckite granite, etc. Previous studies have reported their crystallization ages cluster between 280 and 330 Ma (Table 1), indicating that the alkaline magmatism in the East Junggar lasted for a relatively short period, spanning the Middle-Late Carboniferous to Early Permian. The Huangyangshan pluton is a representative result of magmatism in this granite belt.

Table 1
www.frontiersin.org

Table 1. Statistical table of isotopic ages of typical plutons in Kalamaili area.

Structures in the Kalamaili area are dominated by faults, with overall structural lines oriented NW-SE direction and characterized by multi-stage activities. From north to south, the main faults include the southern segment of the Zaheba-Almantai Fault, the Kupu Fault, the Qingshui-Sujiquan Fault, and the Kalamaili Fault. These faults dissect stratigraphic units and control the formation of most intrusive rocks and ore deposits in the area (Figure 1). Polymetallic mineralization dominated by Au, Cu, and Cr is mainly distributed along the Kalamaili fault zone, with faults exhibiting a clear ore-controlling role (Xu, 2010; Sun, 2018). In contrast, graphite and tin mineralization are closely associated with the alkaline granite belt distributed along the northern margin of the Kalamaili Fault, and have close relationship with magmatism (Tang et al., 2007; Chen W. et al., 2018; Ai et al., 2020). The study object of this paper, the Huangyangshan super-large graphite deposit, occurs within the Huangyangshan alkaline intrusion.

3 Deposit geology

The Huangyangshan pluton, which hosts the graphite deposit, is situated on the southeastern side of the Kalamaili granite belt and the northeastern side of the Qingshui-Sujiquan Fault. It exhibits an overall NW-SE trending elliptical shape and intrudes into tuffaceous sandstones and siltstones of the Middle Devonian Beitaishan Formation and Lower Carboniferous Jiangbasitao Formation, characterized by a concentric ring structure (Figure 2a). Based on mineral grain size and types of dark-colored minerals, the pluton is divided into five lithofacies. All lithofacies share characterized of high silica, low aluminum, depleted calcium and magnesium, enriched alkalis, and high differentiation. They all belong to the high-K calc-alkaline series and representing typical alkaline A-type granite (Yang et al., 2009; Guo et al., 2010; Ai et al., 2020). Numerous dark-colored fine-grained intermediate-mafic enclaves are distributed within the pluton, predominantly in riebeckite and amphibole alkaline granites. Petrographic characteristics indicate a magmatic crystallization origin for these enclaves, while mineral geochemical studies reveal that the Huangyangshan pluton experienced significant mixing of mafic mantle-derived magma and felsic crust-derived magma (Yang et al., 2011; Sun et al., 2021). At present, due to the geological characteristics of the Huangyangshan deposit overlapping with two genetic types, its genesis remains debated: One view classifies it as magmatic type (Bai et al., 2018; Ren et al., 2022), and the other as hydrothermal type (Li et al., 2015; Ai et al., 2020).

Figure 2
Geological map detailing various rock formations and structures in an area. It includes color-coded geological units, such as medium to fine-grained biotite granite, diorite veins, and graphite-bearing amphibole granite. The map features faults, ore veins, exploration lines, and sample locations, with an inset showing a close-up of specific exploration sites. A cross-section below illustrates subsurface rock layers and drill hole positions like ZK704 and ZK701. A legend clarifies symbols and colors, indicating formations and features. The map includes scale bars and directional arrows for orientation.

Figure 2. (a) Geological sketch map of the Huangyangshan pluton (modified after Bai et al., 2018). (b) Geological map of the No.1 and 2 ore bodies of the Huangyangshan graphite deposit. (c) Cross-section along line 7–7’ in (b) (modified after Sun et al., 2020).

To date, a total of nine graphite orebodies have been discovered within the Huangyangshan pluton, containing three large or larger scale ones. Orebody ① and ② are hosted at the boundary between medium-to-fine-grained amphibole granite and fine-grained biotite granite, and Orebody ④ occurs within fine-grained biotite granite. All orebodies are lithologically characterized as medium-to-fine-grained graphite-bearing alkaline granite (Figure 3a). Orebody ① is a concealed orebody with weak surface-exposed graphite mineralization. It exhibits a lenticular shape, wide in the middle and narrow in the N-S direction (Figure 2b), and occurs with a gentle dip (Figure 2c). Orebody ② is an outcropping ore body with whole-rock mineralization, extending approximately E-W. It presents a saddle-like shape at the surface (Figure 2b) and transitions to a steep dip at depth (Figure 2c).

Figure 3
Geological image set showing various mineral formations and rock textures in different panels. Panel a shows a boundary between fine-grained biotite granite and a graphite orebody. Panel b displays a rock with graphite formations labeled

Figure 3. (a) Contact boundary between the orebody and country rocks at surface in the Huangyangshan deposit; (b) Spherical and irregular aggregates of graphite; (c) Primary spherical graphite and its internal silicate mineral cores in drill cores; (d) Subhedral-anhedral amphibole in graphite ore (−); (e) Flaky biotite replaced by graphite (−); (f) Sulfides enclosed in spherical graphite ore from drill cores; (g) Anhedral pyrrhotite containing chalcopyrite and pentlandite in graphite ore (−); (h) Flaky graphite associated with albitization (+); (i) Veined graphite and altered biotite within it (−). Gr - Graphite; Amp - Amphibole; Bt - Biotite; Kfs - K-feldspar; Qz - Quartz; Pl - Plagioclase; Ab - Albite; Po - Pyrrhotite; Ccp - Chalcopyrite; Pn - Pentlandite.

The main ore mineral in the orebodies is graphite, mostly occurring as flaky or scaly aggregates forming spherical or irregular masses (Figure 3b). Associated silicate minerals—predominantly feldspar, quartz, amphibole, biotite, and pyroxene—are commonly found within these graphite spheres (Figure 3c). They are closely associated with graphite and cut by it. Amphibole occurs as euhedral-subhedral prismatic grains, with partial chloritization (Figure 3d). Biotite appears as brown flakes (Figure 3e) and is replaced and cut by graphite, suggesting that graphite formed at post-magmatic stage. Symbiotic metal sulfides are predominantly developed in deep ores (Figure 3f), mainly including pyrrhotite, chalcopyrite, pentlandite, and ilmenite. Their distribution shows a positive correlation with graphite (higher graphite content corresponds to higher pyrrhotite content), suggesting a close genetic relationship and graphite forming after pyrrhotite crystallization (Figure 3g). Based on petrographic observation, most flake graphite is closely associated with albitization (Figure 3h), indicating a high-temperature hydrothermal crystallization origin. Furthermore, within spheroidal graphite, veined graphite is observed filling interstices between quartz and feldspar in an anhedral interstitial texture (Figure 3i). Some low-temperature hydrothermal alterations like chloritization and argillization are developed within these veins, suggesting formation during the late hydrothermal evolution stage.

4 Samples and analytical methods

The graphite-bearing alkaline granite sample (ZK701) used for zircon U-Pb dating in this study was collected from the drill core of Orebody ② (Figure 2c), with weak weathering. The sample is mainly composed of graphite (∼10%), quartz (∼40–50%), K-feldspar (∼30–40%), amphibole (∼4–8%), and biotite (<3%), with partial albitization, chloritization and carbonatization alteration.

Zircon separation, mount preparation, and cathodoluminescence (CL) imaging were conducted by Beijing Zircon Age Leading Technology Co., Ltd. Initially, zircon grains with good transparency, high crystallinity, and relatively large grain sizes were handpicked under microscope. These grains were then embedded in epoxy resin, mounted on a circular stub, and polished until their cores were exposed, followed by CL imaging. Based on comparisons of zircon grains in the CL images, those with well-developed crystal forms, high transparency, minimal fractures and inclusions, and distinct zoning were selected for subsequent analysis. The U-Th-Pb isotopies and trace element analysis were performed simultaneously on zircon, using an Agilent 7900 type inductively coupled plasma-mass spectrometry (ICP-MS) equipped with a GeolasPro 193 nm ArF excimer laser ablation (LA) system at the Key Laboratory of Mineral Resources Evaluation in Northeast Asia, Ministry of Natural Resources. Helium (He) was used as the carrier gas for ablated materials, and the synthetic silicate glass standard reference material NIST 610, developed by the National Institute of Standards and Technology (USA), was employed as the reference for element concentrations. The analyses were conducted using a laser repetition rate of 7 Hz, an energy density of 3 J/cm2, and a spot size of 32 μm. For dating purposes, the zircon standard 91500—with its U-Th-Pb isotopic ratios referenced from Wiedenbeck et al. (1995)—served as the external isotopic calibrant, analyzed twice every 8–10 sample analyses. For each analysis, time-resolved signals were examined to detect the presence of inclusions or common Pb. Zircon trace element concentrations were calibrated using the glass standard NIST 610 as the external standard, with the theoretical Si content of zircon employed as the internal standard. Data processing was performed using the software ICPMSDataCal (Liu et al., 2008; Liu et al., 2010), and age calculations were carried out using the program Isoplot/Ex v. 3.0 (Ludwig, 2003).

5 Results

5.1 Zircon petrography and trace element compositions

Based on morphological, color, and cathodoluminescence (CL) texture characteristics, zircons from the graphite-bearing alkaline granite can be divided into two distinct types (Figure 4). Type I zircons are euhedral to subhedral prismatic grains with a relatively light color (white or light gray) and larger grain sizes, with most having a long axis exceeding 100 μm. They display distinct oscillatory zoning (Figure 4a). Some zircons display core-rim textures with dark rims, and certain zircons contain inherited cores and mineral inclusions (Figure 4a), indicating a magmatic origin with subsequent intense hydrothermal alteration. CL image characteristics of Type II zircons significantly differ from those of magmatic zircons. These zircons are dark black in color, mostly occurring as euhedral short prisms or granular grains with relatively small sizes (long axis generally ranging from 70 to 130 μm). They either lack oscillatory zoning or have indistinct zoning, with a small number of inclusions visible internally (Figure 4b).

Figure 4
Ten zircon crystals with ages labeled range from 256 to 327 million years. Crystals are split into two sections: (a) with five crystals aged 315 to 327 million years, and (b) with five crystals aged 256 to 300 million years. Each crystal is marked with a yellow circle, and the scale indicates one hundred micrometers.

Figure 4. CL images and analytical spot locations of representative zircons from the graphite-bearing granite in drillhole ZK701. (a) CL images and analytical spot locations of pale magmatic zircons. (b) CL images and analytical spot locations of dark hydrothermal zircons.

The REE compositions of zircons (analyzed by LA-ICPMS) are presented in Table 2. Type I zircons have relatively low Th and U contents, with Th ranging from 34 to 137 ppm (mostly <100 ppm) and U from 101 to 331 ppm (mostly <200 ppm), yielding Th/U ratios of 0.30–0.48. They contain total rare earth elements (∑REE) of 652–1784 ppm, with an average of 963 ppm, characterized by depletion in light rare earth elements (LREE) and enrichment in heavy rare earth elements (HREE). The LREE contents vary from 6 to 75 ppm, whereas the HREE contents range from 611 to 1735 ppm, resulting in LREE/HREE ratios of 0.01–0.1 (Table 2). These zircons exhibit similar chondrite-normalized REE patterns with strong positive Ce anomalies (Ce/Ce* = 7.61–316.24). Additionally, they have extremely low La and Pr contents, showing significant positive Ce anomalies and relatively strong negative Eu anomalies in the chondrite-normalized REE distribution diagram (Figure 5a).

Table 2
www.frontiersin.org

Table 2. Zircon rare earth trace element composition of the graphite-bearing granite in drillhole ZK701.

Figure 5
Two line graphs compare rare earth element patterns in zircon samples. Graph (a) shows magmatic zircon with upward trends and graph (b) shows hydrothermal zircon with distinct patterns. Both use a logarithmic scale for sample-to-chondrite ratios, with elements including Lanthanum, Cerium, and Lutetium on the x-axis.

Figure 5. Chondrite-normalized REE patterns for the zircons from the graphite-bearing granite in drillhole ZK701 normalized to the chondrite composition of Sun and McDonough (1989). (a) Chondrite-normalized REE patterns of magmatic zircons. (b) Chondrite-normalized REE patterns of hydrothermal zircons.

Type II zircons exhibit extremely high Th and U contents, with Th ranging from 1149 to 8724 ppm (mostly >3000 ppm) and U from 2124 to 11693 ppm (mostly >4000 ppm), yielding Th/U ratios between 0.53 and 1.01—significantly higher than those of magmatic zircons. Compared to Type I zircons, these zircons generally have higher rare earth element (REE) contents (∑REE = 1192–24533 ppm, averaging 6168 ppm) (Table 2). They display less pronounced fractionation between light and heavy REEs, with LREE contents varying from 76 to 13385 ppm and HREE contents ranging from 995 to 11148 ppm, resulting in LREE/HREE ratios of 0.02–1.2. In the chondrite-normalized REE distribution diagram (Figure 5b), their patterns are relatively flat. Y contents (mostly >1%) are significantly higher than those in Type I zircons (mostly 0.1%–0.2%). Additionally, Type II zircons have high La contents (2.56–1609.66 ppm) and low (Sm/La)N (0.59–8.19) and Ce/Ce* (1.8–3.96) values, indicating weaker positive Ce anomalies.

5.2 Zircon U-Pb dating

Due to the presence of two zircon types, a total of 45 zircon grains were selected for analysis from sample ZK701, including 18 magmatic zircons and 27 hydrothermal zircons, yielding 43 valid zircon spot ages (Table 3). These spot ages range from 256 ± 6 Ma to 327 ± 9 Ma, all lying on or near the U-Pb concordia line (Figure 6). Based on their zircon types and age distribution characteristics, they can be clearly divided into three groups:

Table 3
www.frontiersin.org

Table 3. LA-ICP-MS zircon U-Pb dating data of graphite-bearing granite in drillhole ZK701 (1σ values represent the relative uncertainties).

Figure 6
Scatter plot with ellipses illustrating ^{206}Pb/^{238}U versus ^{207}Pb/^{235}U data and a regression line. Inset bar charts show ages with means and uncertainties: 318.8±4.0 (n=18), 299.9±4.6 (n=12), and 259.9±3.3 (n=13).

Figure 6. U–Pb concordia diagrams and weighted mean ages for the magmatic and hydrothermal zircons from the graphite-bearing granite in drillhole ZK701.

Group 1 consists of 18 magmatic zircons with relatively older ages ranging from 315 ± 9 Ma to 327 ± 9 Ma. The weighted mean age is 318.8 ± 4.0 Ma (MSWD = 0.10), representing the magmatic crystallization age of the graphite-bearing medium-fine-grained alkaline granite. Group 2 comprises 12 hydrothermal zircons with ages concentrated between 297 ± 6 Ma and 303 ± 8 Ma. The weighted mean age is 299.9 ± 4.6 Ma (MSWD = 0.052), slightly younger than the magmatic crystallization age of the graphite-bearing granite and consistent with the crystallization age of the fine-grained biotite alkali-feldspar granite phase (Sun et al., 2021), indicating that the hydrothermal event forming these zircons occurred in the Late Carboniferous. Group 3 includes 13 hydrothermal zircons with ages ranging from 256 ± 6 Ma to 267 ± 7 Ma. The weighted mean age is 259.9 ± 3.3 Ma (MSWD = 0.31), significantly younger than the pluton crystallization age, suggesting a later overprinting hydrothermal event.

6 Discussion

6.1 Genesis of the two types of zircons

During the crystallization differentiation of highly fractionated granites, the crystallization of anhydrous or low-water minerals gradually increases the water content of the system, leading to a magmatic-hydrothermal transitional stage. Zircons can crystallize throughout this process, forming both magmatic and hydrothermal zircons. Notably, magmatic zircons formed in the early magmatic stage may be altered by later exsolved fluids to form hydrothermal zircons, ultimately resulting in the high abundance of hydrothermal zircons in highly evolved granites (Reed et al., 2000; Meng et al., 2024). This phenomenon serves as a direct indicator that the magmatic system has entered the magmatic-hydrothermal evolutionary stage.

According to previous studies, magmatic zircons are characterized by colorless transparency, distinct oscillatory zoning, low Th and U contents, extremely low LREE contents, significant LREE-HREE fractionation, and pronounced positive Ce anomalies and negative Eu anomalies (Hoskin, 2005; Fu et al., 2009; Li et al., 2014). In contrast, hydrothermal zircons typically appear as dark brown or dark translucent crystals with spotty or spongy textures, contain diverse mineral inclusions, exhibit vague or absent oscillatory zoning, have relatively high LREE contents, weak LREE-HREE fractionation, and absent or weak Ce and Eu anomalies (Geisler et al., 2003; Kozlik et al., 2016; Takehara et al., 2018; Jiang et al., 2019). Based on these criteria, Type I zircons in this study can be inferred as magmatic and Type II zircons as hydrothermal (Table 2). This classification is further supported by geochemical diagrams. In the ΣREE vs. ΣLREE and U vs. Th diagrams (Figures 7a,b), Type I zircons (with low element contents) and Type II zircons (with high element contents) plot in the magmatic zircon and hydrothermal zircon fields, respectively. In the La versus (Sm/La)N and (Sm/La)N versus Ce/Ce* diagrams (Figures 7c,d), most data points of Type I and Type II zircons fall within or near the fields of magmatic and hydrothermal zircons, respectively. These observations further confirm the presence of two zircon types in the Huangyangshan graphite deposit, indicating intense hydrothermal evolution during the mineralization process.

Figure 7
Four scatter plots comparing different chemical properties of Type I and Type II zircons. (a) Shows LREE versus REE with a magmatic to hydrothermal trend. (b) Displays Th versus U, indicating Th/U ratios. (c) Illustrates Sm/La versus La with shaded magmatic and hydrothermal zircons. (d) Presents Ce/Ce* versus Sm/La, similarly highlighting magmatic and hydrothermal zones. Green and black dots differentiate the zircon types.

Figure 7. Discriminant diagrams for the zircons from the graphite-bearing granite in drillhole ZK701. (a) ΣREE vs. ΣLREE diagram; (b) U vs. Th diagram; (c) La vs. (Sm/La)N and (d) (Sm/La)N vs. Ce/Ce* diagram (base discrimination diagrams after Hoskin, 2005).

It should be noted that typical hydrothermal zircons usually exhibit irregular crystal forms, well-developed internal pores, and dissolution features—traits that do not fully align with the characteristics of Type II zircons. Instead, Type II zircons have relatively well-formed euhedral forms and few inclusions. Additionally, some Type II zircons contain residual light-colored magmatic zircon cores (Figure 4b), indicating that these hydrothermal zircons did not crystallize directly from Zr-saturated hydrothermal fluids but were likely altered or modified by hydrothermal fluids. This is consistent with the core-rim textures observed in some magmatic zircons (Figure 4a), suggesting that the dark rims formed due to hydrothermal alteration after the crystallization of magmatic zircons. After hydrothermal metasomatism, the Th, U, and LREE contents of the original light-colored magmatic zircons increased, which is the main reason for the darkening of hydrothermal zircon. Type II zircons represent the hydrothermal disturbance events experienced during the late evolutionary stage of magmatic zircon evolution.

6.2 Graphite mineralization age

Previous studies have shown that hydrothermal zircons can either crystallize directly from medium-low temperature hydrothermal fluids or form via alteration along the margins or fractures of primary zircons. This process may reset the U-Pb chronometric system of altered zircons and effectively recording the timing of hydrothermal activity (Rubin et al., 1993; Van Lichtervelde et al., 2009; Xia et al., 2021). Petrographic studies indicate that hydrothermal alterations in the Huangyangshan graphite ores are dominated by high-temperature albitization, greisenization, and low-temperature chloritization and sericitization. Among these, high-temperature hydrothermal alterations such as albitization are paragenetic with graphite ores of the main mineralization stage, whereas low-temperature hydrothermal alterations are mostly distributed in granitic matrices and some light-colored felsic cores, showing a weak genetic relationship with spherical graphite (Sun, 2022). This suggests that graphite mineralization is associated with high-temperature hydrothermal processes.

All samples in this study were collected from drill cores rich in graphite of the main mineralization stage, and the hydrothermal fluid alteration features of Type II hydrothermal zircons are consistent with the aforementioned alteration scenario. Therefore, hydrothermal activity is closely related to graphite mineralization, and the ages of hydrothermal zircons can reflect the crystallization age of graphite. Additionally, hydrothermal zircons mostly form under medium-low temperature conditions (<450 °C) and only zircons subjected to thorough hydrothermal alteration can be used for micro-area dating (Zhu and Song, 2006; Li, 2009). However, graphite in the Huangyangshan deposit formed during an earlier high-temperature hydrothermal stage, so the measured ages of hydrothermal zircons might be slightly younger than the actual crystallization age of graphite.

Due to the lack of suitable dating minerals, previous studies on the Huangyangshan graphite deposit primarily relied on zircon U-Pb dating of the ore-forming pluton to indirectly constrain the mineralization age (Ai et al., 2020; Bai, 2021). The authors previously conducted Re-Os isotopic dating on Huangyangshan graphite, yielding an isochron age of 332 ± 53 Ma with a large error range (Sun et al., 2021)—this can only serve as a reference for corroborating the mineralization age. In this study, three groups of zircon ages were obtained, each with distinct geological chronological significance: (1) Magmatic zircon age (318.8 ± 4.0 Ma). As the only magmatic zircon age group in the graphite ores, it represents the magmatic crystallization age of the graphite-bearing alkaline granite. This age is consistent with that of the medium-fine-grained amphibole granite surrounding Orebodies ① and ② (318.3 ± 4.0 Ma; Sun, 2022), indicating that the graphite-bearing alkaline granite and the country rocks formed contemporaneously as products of the same magmatic event, both emplaced during the early Late Carboniferous. (2) Earlier hydrothermal zircon ages (299.9 ± 4.6 Ma). Slightly younger than the magmatic crystallization age but close to the crystallization age of the fine-grained biotite granite (301.1 ± 3.6 Ma)—a lithofacies closely associated with graphite orebodies—suggesting that the formation of this lithofacies may relate to this hydrothermal event. Flaky graphite aggregates cut through silicate minerals and are closely paragenetic with high-temperature hydrothermal alterations such as albitization (Figures 3d,e). Thus, the graphite formed during the high-temperature hydrothermal stage (not long after magmatic crystallization). In addition, the temperature of fluid inclusions developed in the main graphite metallogenic period is above 450 °C, and their H-O isotope results are projected near the magmatic water area, indicating a magmatic hydrothermal origin (Ai et al., 2020). Thus, this hydrothermal event most likely originated from magmatic fluids exsolved during late magmatic evolution stage, sharing a genetic link with graphite formation, and effectively constrains the lower limit of the main mineralization stage for graphite. (3) Later hydrothermal zircon ages (259.9 ± 3.3 Ma). Significantly different from the crystallization age of the Huangyangshan pluton but consistent with the mineralization ages of regional hydrothermal deposits such as the Shuangquan and Kubusu gold deposits (Xu, 2010; Han, 2014)—indicating a regional hydrothermal event of sufficient scale to affect the Huangyangshan pluton. Previous studies suggest these gold deposits along the Kalamaili Fault are typical orogenic gold deposits, with fluid sources dominated by metamorphic water mixed with minor magmatic or meteoric water (Xu et al., 2010). Therefore, this hydrothermal event in the Huangyangshan graphite orebodies may also relate to metamorphic fluids transported by the Kalamaili Fault and its secondary fractures, though its genetic relationship with graphite requires further investigation.

The Re-Os isotopic dating result for graphite is approximately 330 Ma, but no ∼330 Ma geological bodies have been found in the Huangyangshan pluton to date—the oldest crystallization age of the medium-fine grained riebeckite granite is only 322.7 ± 4.5 Ma (Sun, 2022). Although magmatic metallic sulfides and phenocrysts of silicate minerals such as pyroxene and plagioclase in graphite ores formed during early magmatic evolution, mineral paragenetic relationships indicate that graphite mineralization was significantly later than these early-stage minerals, suggesting graphite’s crystallization age is much younger than 330 Ma.

In this study, zircons extracted from graphite ores exhibit a dichotomous characteristic in both age and genesis—distinctly different from previous studies, where only pale-colored magmatic zircons were found in graphite-free plutons (Ai et al., 2020; Sun et al., 2021). Magmatic zircons in the pluton show consistent morphology, age, and trace element characteristics, with rare captured zircons. This indicates low interference from ancient zircons. For a super-large deposit like the Huangyangshan deposit, the large scale of the ore-forming pluton and the relatively long duration of magmatic-hydrothermal processes are also key factors contributing to this age dichotomy. Based on the ages of magmatic and hydrothermal zircons in graphite ores, as well as the crystallization age of the fine-grained biotite granite that is closely related to graphite, the orebody’s magmatic crystallization age is 318.8 ± 4.0 Ma, the age of the hydrothermal event associated with graphite mineralization is 299.9 ± 4.6 Ma, and the fine-grained biotite granite’s age is 301.1 ± 3.6 Ma (Sun, 2022). These ages constrain that the graphite crystallized in the early high-temperature hydrothermal stage is no older than the magmatic crystallization age and no younger than the crystallization age of hydrothermal zircons.

Therefore, this study proposes that the Huangyangshan graphite deposit formed ∼ 300 Ma, contemporaneous with the late-stage biotite granite in the Huangyangshan pluton. This age is consistent with the Rb-Sr age of the Ganliangzi tin deposit (305 ± 25 Ma, Chen et al., 1999) and the molybdenite Re-Os age of the Sareshike tin deposit (307 ± 11 Ma, Tang et al., 2007)—both related to alkaline granites in this region. Therfore, the Huangyangshan graphite deposit formed in the late Late Carboniferous.

6.3 Hydrothermal processes in the Huangyangshan graphite deposit

Based on the characteristics of zircons in graphite-bearing granite, the Huangyangshan graphite deposit underwent an evolutionary process from magmatism to hydrothermal activity—findings that hold great significance for resolving the dispute over the genesis of the Huangyangshan deposit. In terms of whole-rock compositions (Sun et al., 2021), the four earlier-crystallized lithofacies of the Huangyangshan pluton are characterized by significant LREE enrichment and relatively high positive εHf(T) and εNd(T) values. This suggests they were primarily derived from juvenile crustal materials at the base of the lower crust and experienced extensive magma mixing and fractional crystallization. However, compared with other lithofacies, the youngest fine-grained biotite granite—whose age is close to that of the first group of hydrothermal zircons—exhibits similar major element characteristics but distinctly different REE distribution patterns. This indicates that while sharing a common magmatic source, these lithofacies underwent distinct late-stage evolutionary processes.

Despite having a relatively low total REE content, the fine-grained biotite granite displays an M-type tetrad effect in its REE distribution patterns (Figure 8). The “tetrad effect” of REE is a unique REE pattern in granites, where REEs are divided into four groups (La-Ce-Pr-Nd, Pm-Sm-Eu-Gd, Gd-Tb-Dy-Ho, and Er-Tm-Yb-Lu) using Nd/Pm, Gd, and Ho/Er as boundaries. These groups form four convex or concave curves, termed M-type and W-type REE tetrad effects, respectively (Masuda et al., 1987). This phenomenon is commonly observed in REE-rich alkaline granites and is generally associated with highly evolved magmatic systems, magma interactions with volatile-rich fluids, and fluid (including hydrothermal) activities (Zhao et al., 1999; Zhao et al., 2010). Calculations of REE tetrad effect parameters for the fine-grained biotite granite (Table 4) show that the central element parameters of the first subgroup (CeN/Ce*, PrN/Pr*) and third subgroup (TbN/Tb*, DyN/Dy*) are all greater than 1, exhibiting a distinct M-type. The central element parameters of the fourth subgroup (TmN/Tm*, YbN/Yb*) are approximately 1, indicating an insignificant tetrad effect. The second subgroup deviates from the convex curve due to a strong negative Eu anomaly (Figure 8). Compared with other highly fractionated granites with REE tetrad effects in China, the Huangyangshan fine-grained biotite granite has smaller parameters such as CeN/Ce* and PrN/Pr* but a more significant negative Eu anomaly, with REE distribution patterns similar to those of the Qianlishan and Baerzhe plutons (Zhao et al., 1999).

Figure 8
Logarithmic scale graph showing Sample/Chondrite ratios for various elements. Pink circles represent fine-grained biotite granite, blue and green circles indicate other lithofacies and the Huangyangshan pluton. Notable variation occurs at Europium (Eu) with a dip in values. Elements La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu are plotted on the x-axis.

Figure 8. Chondrite-normalized REE patterns for different granite phases of the Huangyangshan pluton (data from Sun et al., 2021).

Table 4
www.frontiersin.org

Table 4. REE tetrad effect parameters of the fine-grained biotite granite (data from Sun et al., 2021).

Currently, two main viewpoints explain the origin of the REE tetrad effect in granites: (1) Highly fractional crystallization of granitic magma and the crystallization of REE-rich accessory minerals such as monazite, xenotime, and apatite can cause residual melts to exhibit the tetrad effect (Yurimoto et al., 1990; Seward and Barnes, 1997); (2) Intense melts-fluids interactions can induce the tetrad effect in rocks (Bau, 1996; Kawabe, 1999). As mentioned above, while the Huangyangshan pluton underwent significant fractional crystallization during its overall evolution, other lithofacies do not exhibit the REE tetrad effect—suggesting fractional crystallization is not the primary driver of this phenomenon. Additionally, the fine-grained biotite granite lacks typical REE minerals. In contrast, hydrothermal Sn deposits occur in the riebeckite granite phase in the northern part of the pluton, further indicating that REE minerals is not the main cause of the REE tetrad effect.

Rare earth elements are generally considered stable, immobile elements in magmas. However, existing studies have shown that REE mobility increases significantly with the risng proportion of hydrothermal fluids in melts or rocks (Alderton et al., 1980; Veksler et al., 2005). Compared with other lithofacies, the fine-grained biotite granite has a significantly lower total REE content (ΣREE = 40.54–43.16 ppm), with a particularly marked decrease in LREEs (LaN/YbN = 0.64–0.85) (Sun et al., 2021), indicating extensive REE precipitation during this magmatic stage. The intensity of the REE tetrad effect can be expressed by the parameter TE1 [=(Ce/Ce* × Pr/Pr*)0.5]. The K/Rb and La/Yb ratios of melts reflect the degree of magmatic fractional crystallization, while the Zr/Hf and Y/Ho ratios of rocks indicate melt-fluid interactions. Lower Zr/Hf ratios and higher Y/Ho ratios imply stronger non-CHARAC trace element behavior in melts and more significant fluid activity (Bau, 1996; Chen Y. F. et al., 2018). Plots of TE1 against these elemental ratios (Figures 9a,b) reveal a clear correlation between TE1 and Zr/Hf and Y/Ho ratios, indicating a higher degree of REE tetrad effect corresponds to stronger non-CHARAC behavior in melts and more significant fluid influence. In contrast, no obvious correlation exists between TE1 and K/Rb or La/Yb ratios (Figures 9c,d), suggesting no direct link between the degree of magmatic fractional crystallization and the intensity of the REE tetrad effect. Melt-fluid interactions induced REE fractionation (manifested as the tetrad effect) and disrupted the stability of the magmatic environment, causing significant system decompression and cooling. This conclusion is further corroborated by the low-temperature crystallization stage experienced by Mg-biotite and ferrohornblende formed via hydrothermal alteration in the orebody (Sun et al., 2022). Cooling is the most effective mechanism for graphite precipitation (Luque and Rodas, 1999; Luque et al., 2012) and the dominate factor driving graphite precipitation in the Huangyangshan deposit (Sun, 2022). Consequently, carbonaceous materials precipitated and crystallized to form graphite during the fluid-melt interaction process. This indicates that the melt-fluid metasomatism responsible for the REE tetrad effect in the Huangyangshan fine-grained biotite granite may be directly related to the hydrothermal activity within the graphite orebody.

Figure 9
Four scatter plots labeled (a) to (d) analyze fine-grained biotite granite. Plot (a) shows Zr/Hf values decreasing with TE1. Plot (b) shows Y/Ho values increasing with TE1. Plot (c) depicts K/Rb with minimal variation across TE1. Plot (d) presents La/Yb with little variation in TE1 values. Pink circles represent data points, and gray arrows indicate trends.

Figure 9. TE1 vs. Zr/Hf (a), Y/Ho (b), K/Rb (c) and La/Yb (d) diagrams of the Huangyangshan fine-grained biotite granite (data from Sun et al., 2021). Note that the K/Rb and La/Yb ratios reflect the degree of magmatic fractional crystallization, while the Zr/Hf and Y/Ho ratios indicate melt-fluid interactions. The correlation between TE1 and those ratios can reflect the connection of REE tetrad effect and above factors.

The tetrad effect of fine-grained biotite granite indicates its evolution is related to fluid-melt interaction. Meanwhile, hydrothermal zircons in the graphite-bearing alkaline granite confirm the pluton underwent a significant post-magmatic hydrothermal evolution process. Moreover, the zircon U-Pb ages are consistent within error. Therefore, the formation of the fine-grained biotite granite is closely related to the hydrothermal event around 300 Ma recorded by the first stage of hydrothermal zircons and the formation of graphite orebodies. This confirms that hydrothermal events related to graphite mineralization existed during the late magmatic evolution stage of the Huangyangshan graphite deposit.

7 Conclusion

1. Both magmatic zircons and hydrothermal zircons are present in the ores of the Huangyangshan graphite deposit. These hydrothermal zircons have been altered or modified by hydrothermal fluids, indicating that the Huangyangshan graphite deposit has undergone intense hydrothermal evolution, and the mineralization process may be related to hydrothermal activity.

2. The zircon U-Pb ages can be divided into three groups: magmatic zircon age (318.8 ± 4.0 Ma), early hydrothermal zircon ages (299.9 ± 4.6 Ma), and late hydrothermal zircon ages (259.9 ± 3.3 Ma). They represent the magmatic crystallization age of the orebody, the graphite mineralization age, and a later hydrothermal event related to regional metamorphic fluids, respectively.

3. The M-type REE tetrad effect in the fine-grained biotite granite, driven by strong fluid-melt interaction, lead to significant decompression and cooling, providing a necessary condition for graphite precipitation. This hydrothermal activity around 300 Ma might be closely related to the graphite mineralization.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

XS: Conceptualization, Data curation, Formal Analysis, Investigation, Writing – original draft, Writing – review and editing. YR: Funding acquisition, Investigation, Project administration, Supervision, Writing – review and editing. ZS: Formal Analysis, Investigation, Methodology, Writing – original draft. JL: Data curation, Investigation, Methodology, Writing – original draft. MH: Data curation, Funding acquisition, Visualization, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was financially supported by the National Natural Science Foundation of China (42472138).

Acknowledgments

We sincerely appreciate the geologists from the Xinjiang Branch of China National Geological Exploration Center of Building Materials Industry for their support of our fieldwork.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not 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

Ai, J., Lu, X. B., Li, Z. W., and Wu, Y. L. (2020). Genesis of the graphite orbicules in the Huangyangshan graphite deposit, Xinjiang, China: evidence from geochemical, isotopic and fluid inclusion data. Ore Geol. Rev. 122, 103505. doi:10.1016/j.oregeorev.2020.103505

CrossRef Full Text | Google Scholar

Alderton, D. H. M., Pearce, J. A., and Potts, P. J. (1980). Rare earth element mobility during granite alteration: evidence from southwest England. Earth Planet. Sci. Lett. 49, 149–165. doi:10.1016/0012-821x(80)90157-0

CrossRef Full Text | Google Scholar

Bai, J. K. (2021). Metallogeny of the graphite deposits in the Eastern Junggar area. Xi’an, China: Northwest University.

Google Scholar

Bai, J. K., Chen, J. L., and Peng, S. X. (2018). Geochronology and geochemistry of ore-bearing intrusions from Huangyangshan magmatic hydrothermal graphite deposit in Qitai county, Xinjiang. Acta Petrol. Sin. 34 (8), 2327–2340.

Google Scholar

Bau, M. (1996). Controls on the fractionation of isovalent trace elements in magmatic and aqueous systems: evidence from Y/Ho, Zr/Hf and lanthanide tetrad effect. Contrib. Mineral. Petrol. 123, 323–333. doi:10.1007/s004100050159

CrossRef Full Text | Google Scholar

Belousova, E. A., Griffin, W. L., and O'Reilly, A. S. Y. (2006). Zircon crystal morphology, trace element signatures and Hf isotope composition as a tool for petrogenetic modelling: examples from Eastern Australian granitoids. J. Petrol. 47 (2), 329–353. doi:10.1093/petrology/egi077

CrossRef Full Text | Google Scholar

Chen, F. W., Li, H. Q., Cai, H., Liu, Q., and Chang, H. L. (1999). Chronology and origin of the Ganliangzi tin ore field, Xinjiang. Min. Depos. 18 (1), 91–97.

Google Scholar

Chen, W., Chen, B., and Sun, K. K. (2018). Petrogenesis of the maogongdong highly differentiated granite in the dahutang tungsten ore field,Jiangxi province. Acta Petrol. Sin. 34 (6), 1704–1724.

Google Scholar

Chen, Y. F., Wang, Y. W., Wang, J. B., Wang, L. J., Tang, P. Z., Shi, Y., et al. (2018). Greisenized alteration-mineralization geochemistry of the tin deposit related to A-type granite: case study on the kamusite and ganliangzi deposits, Xinjiang. Earth Sci. 43 (09), 3154–3168.

Google Scholar

Crespo, E., Luque, F. J., Rodas, M., Wada, H., and Gervilla, F. (2006). Graphite-sulfide deposits in Ronda and Beni Bousera peridotites (Spain and Morocco) and the origin of carbon in mantle-derived rocks. Gondwana Res. 9, 279–290. doi:10.1016/j.gr.2005.10.003

CrossRef Full Text | Google Scholar

Doroshkevich, A. G., Wall, F., and Ripp, G. S. (2007). Magmatic graphite in dolomite carbonatite at Pogranichnoe, North Transbaikalia, Russia. Contrib. Mineral. Petrol. 153, 339–353. doi:10.1007/s00410-006-0150-z

CrossRef Full Text | Google Scholar

Duke, E. F., and Rumble, I. I. I. D. (1986). Textural and isotopic variations in graphite from plutonic rocks, south central New Hampshire. Contrib. Mineral. Petrol. 93, 409–419. doi:10.1007/bf00371711

CrossRef Full Text | Google Scholar

Feng, H. Q., Zhang, B. W., Gao, D., Zhan, X. Z., and Zheng, Y. Y. (2019). Geochronology, geochemical characteristics and tectonic setting of the kamusite A2-type granite in the Eastern Junggar. Bull. Mineral. Petrol. Geochem. 38 (06), 1154–1169.

Google Scholar

Fu, B., Mernagh, T. P., Kita, N. T., Kemp, A. I. S., and Valley, J. W. (2009). Distinguishing magmatic zircon from hydrothermal zircon: a case study from the gidginbung high-sulphidation Au-Ag-(Cu) deposit, SE Australia. Chem. Geol. 259, 131–142. doi:10.1016/j.chemgeo.2008.10.035

CrossRef Full Text | Google Scholar

Gan, L., Tang, H. F., and Han, Y. J. (2010). Geochronology and geochemical characteristics of the yemaquan granitic pluton in East Junggar, Xinjiang. Acta Petrol. Sin. 26 (8), 2374–2388.

Google Scholar

Gao, J. G., Liang, T., Feng, Y. G., Huang, F., Wang, D. H., San, J. Z., et al. (2023). Genetic relationship between the complex massif and pegmatites in halong rare metal ore⁃concentrated area, Altay. Acta Petrol. Sin. 39 (7), 1908–1930. doi:10.18654/1000-0569/2023.07.03

CrossRef Full Text | Google Scholar

Gao, F., Pei, X. Z., Cheng, Y. F., Guo, R. Q., Sun, C. Y., Liu, X. Q., et al. (2024). Petrogenesis of the early Permian A-type granites in eastern North Tianshan, Central Asian orogenic belt: implications for the final geodynamic evolution of the Kalamaili suture zone. J. Asian Earth Sci. 264, 106057. doi:10.1016/j.jseaes.2024.106057

CrossRef Full Text | Google Scholar

Geisler, T., Rashwan, A. A., Rahn, M. K. W., Poller, U., Zwingmann, H., Pidgeon, R. T., et al. (2003). Low-temperature hydrothermal alteration of natural metamict zircons from the eastern desert, Egypt. Mineral. Mag. 67, 485–508. doi:10.1180/0026461036730112

CrossRef Full Text | Google Scholar

Geisler, T., Schaltegger, U., and Tomaschek, F. (2007). Re-equilibration of zircon in aqueous fluids and melts. Elements 3 (1), 43–50. doi:10.2113/gselements.3.1.43

CrossRef Full Text | Google Scholar

Guo, F. F., Jiang, C. Y., Lu, R. H., Xia, Z. D., Ling, J. L., and Guo, N. X. (2010). Petrogenesis of the Huangyangshan alkali granites in Kalamaili area, northern Xinjiang. Acta Petrol. Sin. 26 (8), 2357–2373.

Google Scholar

Han, Q. (2014). Study of metallogenic regularity and mineralization model of Kalamaili mineralization belt in East Junggar. Urumqi: Xinjiang University.

Google Scholar

Han, B. F., Ji, J. Q., Song, B., Chen, L. H., and Zhang, L. (2006). Late Paleozoic vertical growth of continental crust around the Junggar Basin, Xinjiang, China (part I): timing of post-collisional plutonism. Acta Petrol. Sin. 22 (5), 1077–1086.

Google Scholar

Han, Y. J., Tang, H. F., and Gan, L. (2012). Zircon U-Pb ages and geochemical characteristics of the laoyaquan A-type granites in East Junggar,North Xinjiang,China. Acta Min. Sin. 32 (2), 193–199.

Google Scholar

Hoskin, P. W. O. (2005). Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the jack hills, Australia. Geochim. Cosmochim. Acta 69, 637–648. doi:10.1016/j.gca.2004.07.006

CrossRef Full Text | Google Scholar

Hu, W. L. (2016). Geochemical characteristics and geological significance of late carboniferous granites from the laoyaquan in Eastern Junggar, Xinjiang. Lanzhou: Lanzhou University.

Google Scholar

Huizenga, J. M. (2011). Thermodynamic modeling of a cooling C-O-H fluid-graphite system: implications for hydrothermal graphite precipitation. Min. Deposita 46, 23–33. doi:10.1007/s00126-010-0310-y

CrossRef Full Text | Google Scholar

Jiang, W. C., Li, H., Evans, N. J., and Wu, J. H. (2019). Zircon records multiple magmatic-hydrothermal processes at the giant Shizhuyuan W-Sn-Mo-Bi polymetallic deposit, South China. Ore Geol. Rev. 115, 103160. doi:10.1016/j.oregeorev.2019.103160

CrossRef Full Text | Google Scholar

Kawabe, I. (1999). Thermochemical parameters for solution of lanthanide (III) ethylsulphate and trichloride hydrate series: tetrad effects and hydration change in aqua Ln3+ ion series. Geochem. J. 33 (4), 249–265. doi:10.2343/geochemj.33.249

CrossRef Full Text | Google Scholar

Kozlik, M., Raith, J. G., and Gerdes, A. (2016). U-Pb, Lu–Hf and trace element characteristics of zircon from the felbertal scheelite deposit (Austria): new constraints on timing and source of W mineralization. Chem. Geol. 421, 112–126. doi:10.1016/j.chemgeo.2015.11.018

CrossRef Full Text | Google Scholar

Li, C. M. (2009). A review on the minerageny and situ microanalytical dating techniques of zircons. Geol. Surv. Res. 32 (03), 161–174.

Google Scholar

Li, Y. C., Yang, F. Q., Zhao, C. S., Zhang, Y., Yan, S. H., Dai, J. Z., et al. (2007). SHRIMP U-Pb zircon dating of the Beilekudulke pluton in Xinjiang and its geological implications. Acta Petrol. Sin. 23 (10), 2483–2492.

Google Scholar

Li, H., Watanabe, K., and Yonezu, K. (2014). Zircon morphology, geochronology and trace element geochemistry of the granites from the Huangshaping polymetallic deposit, South China: implications for the magmatic evolution and mineralization processes. Ore Geol. Rev. 60, 14–35. doi:10.1016/j.oregeorev.2013.12.009

CrossRef Full Text | Google Scholar

Li, C., Wang, D. H., Zhao, H., Pei, H. X., Li, X. W., Zhou, L. M., et al. (2015). Minerogenetic regularity of graphite deposits in China. Min. Depo. 34 (6), 1223–1236.

Google Scholar

Li, Z. W., Zhang, X. L., Zhang, L., Chen, Z. G., Chen, J. Y., Wu, Y. L., et al. (2019). Geological characteristics and resource potential of Huangyangshan crystalline graphite deposit in Qitai county, Xinjiang. China Non-metallic Min. Indus. S1, 25–34.

Google Scholar

Liu, Y. S., Hu, Z. C., Gao, S., Günther, D., Xu, J., Gao, C. G., et al. (2008). In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard. Chem. Geol. 257, 34–43. doi:10.1016/j.chemgeo.2008.08.004

CrossRef Full Text | Google Scholar

Liu, Y. S., Hu, Z. C., Zong, K. Q., Gao, C. C., Chen, H. H., Xu, J., et al. (2010). Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS. Chin. Sci. Bull. 55, 1535–1546. doi:10.1007/s11434-010-3052-4

CrossRef Full Text | Google Scholar

Liu, W., Liu, X. J., and Liu, L. J. (2013). Underplating generated A- and I-type granitoids of the East Junggar from the lower and the upper oceanic crust with mixing of mafic magma: insights from integrated zircon U–Pb ages, petrography, geochemistry and Nd–Sr–Hf isotopes. Lithos 179, 293–319. doi:10.1016/j.lithos.2013.08.009

CrossRef Full Text | Google Scholar

Liu, X. J., Xiao, W. J., Xu, J. F., Castillo, P. R., and Shi, Y. (2017). Geochemical signature and rock associations of ocean ridge-subduction: evidence from the Karamaili Paleo-Asian ophiolite in east Junggar, NW China. Gondwana Res. 48, 34–49. doi:10.1016/j.gr.2017.03.010

CrossRef Full Text | Google Scholar

Ludwig, K. R. (2003). User's manual for isoplot 3.00: a geochronological toolkit for microsoft excel 4. Berkeley: Berkeley Geochronology Centre Special Publication, 74.

Google Scholar

Luque, F. J., Pasteris, J. D., Wopenka, B., Rodas, M., and Barrenechea, J. F. (1998). Natural fluid-deposited graphite: mineralogical characteristics and mechanisms of formation. Am. J. Sci. 298, 471–498. doi:10.2475/ajs.298.6.471

CrossRef Full Text | Google Scholar

Luque, F. J., Crespo, E., Barrenechea, J. F., and Ortega, L. (2012). Carbon isotopes of graphite: implications on fluid history. Geosci. Front. 3 (2), 197–207. doi:10.1016/j.gsf.2011.11.006

CrossRef Full Text | Google Scholar

Luque, F. J., Huizenga, J. M., Crespo, E., and Wada, H. (2014). Vein graphite deposits: geological settings, origin, and economic significance. Min. Depo. 49, 261–277. doi:10.1007/s00126-013-0489-9

CrossRef Full Text | Google Scholar

Luque, F. J., and Rodas, M. (1999). Constraints on graphite crystallinity in some Spanish fluid-deposited occurrences from different geologic settings. Mine. Depo. 34, 215–219.

CrossRef Full Text | Google Scholar

Manoel, T. N., and Leite, J. A. D. (2018). On the origin of the Neoproterozoic peresopolis graphite deposit, Paraguay belt, Brazil. J. South Amer. Earth Sci. 84, 104–112. doi:10.1016/j.jsames.2018.03.007

CrossRef Full Text | Google Scholar

Masuda, A., Kawakami, O., Dohmoto, Y., and Takenaka, T. (1987). Lanthanide tetrad effects in nature: two mutually opposite types, W and M. Geochem. J. 21 (3), 119–124. doi:10.2343/geochemj.21.119

CrossRef Full Text | Google Scholar

Meng, A. P., Dai, G. Q., Zhang, N., and Liu, Z. C. (2024). Internal morphology and compositional characteristics of zircons from the two-mica granites in the Luozha area, Southern Xizang: insights into the magmatic-hydrothermal evolution of the Himalayan leucogranites. Geotect. Metallogenia 48 (3), 529–546.

Google Scholar

Papineau, D., De Gregorio, B. T., Cody, G. D., Fries, M. D., Mojzsis, S. J., Steele, A., et al. (2010). Ancient graphite in the Eoarchean quartz-pyroxene rocks from akilia in southern West Greenland I: petrographic and spectroscopic characterization. Geochim. Cosmochim. Acta 74 (20), 5862–5883. doi:10.1016/j.gca.2010.05.025

CrossRef Full Text | Google Scholar

Reed, M. J., Candela, P. A., and Piccoli, P. M. (2000). The distribution of rare earth elements between monzogranitic melt and the aqueous volatile phase in experimental investigations at 800 °C and 200 MPa. Contrib. Mineral. Petrol. 140 (2), 251–262. doi:10.1007/s004100000182

CrossRef Full Text | Google Scholar

Ren, Y. S., Li, J. M., Sun, X. H., Li, Z. W., and Sun, Z. J. (2022). Genesis of metal sulfides and its significance on graphite mineralization in the Huangyangshan graphite deposit, East Junggar, Xinjiang Province, China. Minerals 12, 1450. doi:10.3390/min12111450

CrossRef Full Text | Google Scholar

Rubin, J. N., Henry, C. D., and Price, J. G. (1993). The mobility of zirconium and other “immobile” elements during hydrothermal alteration. Chem. Geol. 110 (1), 29–47. doi:10.1016/0009-2541(93)90246-f

CrossRef Full Text | Google Scholar

Schaltegger, U. (2007). Zircon tiny but timely: hydrothermal zircon. Elements 3, 51–79. doi:10.2113/gselements.3.1.51

CrossRef Full Text | Google Scholar

Seward, T. M., and Barnes, H. L. (1997). “Metal transport by hydrothermal ore fluids,” in Geochemistry of hydrothermal ore deposits. 3rd Edition (New York: John Wiley and Sons), 435–486.

Google Scholar

Shao, B. Q. (2018). The geochemical characeristics and deposit genesis analysis of the Huangyangshan graphite deposit in Qitai county of Xinjiang autonomous region. Xi’an, China: Chang’an University.

Google Scholar

Su, H. M., Jiang, S. Y., Shao, J. B., Zhang, D. Y., Wu, X. K., and Huang, X. Q. (2021). New identification and significance of Early Cretaceous mafic rocks in the interior South China block. Sci. Rep. 11, 11396. doi:10.1038/s41598-021-91045-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, W. L. (2018). A study on the mineralization fluid, ore-bearing geochronology and geochemistry of wulunbulake copper deposit in Xinjiang. China: Xi’an: Chang’an University.

Google Scholar

Sun, X. H. (2022). Magmatism, mineralization and ore-prospecting direction of the Huangyangshan super-large graphite deposit in Qitai county, Xinjiang province. Changchun: Jilin University.

Google Scholar

Sun, S. S., and McDonough, W. F. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geol. Soc. Spec. Publ. 42 (1), 313–345. doi:10.1144/gsl.sp.1989.042.01.19

CrossRef Full Text | Google Scholar

Sun, X. H., Ren, Y. S., Sun, Z. J., Li, Z. W., and Wang, C. Y. (2020). Characteristics, ore sources and genesis of Huangyangshan graphite deposit in Qitai county, Xinjiang. Acta Petrol. Sin. 37 (6), 1867–1882.

Google Scholar

Sun, X. H., Ren, Y. S., Sun, Z. J., Wang, C. Y., and Li, Z. W. (2021). Geochronology and geochemical properties of the large-scale graphite mineralization associated with the Huangyangshan alkaline pluton, Eastern Junggar, Xinjiang, NW China. Geochemistry 81 (4), 125820. doi:10.1016/j.chemer.2021.125820

CrossRef Full Text | Google Scholar

Sun, X. H., Ren, Y. S., Li, J. M., Huang, M. J., Sun, Z. J., and Li, Z. W. (2022). Magma evolution and constraints on the graphite mineralization hosted by the Huangyangshan alkaline granite suite in the East Junggar of Xinjiang province: evidence from in situ analyses of silicate minerals. Minerals 12, 1458. doi:10.3390/min12111458

CrossRef Full Text | Google Scholar

Takehara, M., Horie, K., Hokada, T., and Kiyokawa, S. (2018). New insight into disturbance of U-Pb and trace-element systems in hydrothermally altered zircon via SHRIMP analyses of zircon from the duluth gabbro. Chem. Geol. 484, 168–178. doi:10.1016/j.chemgeo.2018.01.028

CrossRef Full Text | Google Scholar

Tang, H. F., Qu, W. J., Su, Y. P., Hou, G. S., Du, A. D., and Cong, F. (2007). Genetic connection of Sareshike tin deposit with the alkaline A-type granites of Sabei body in Xinjiang: constraint from isotopic ages. Acta Petrol. Sin. 8, 1989–1997.

Google Scholar

Van Lichtervelde, M., Melcher, F., and Wirth, R. (2009). Magmatic vs. hydrothermal origins for zircon associated with tantalum mineralization in the tanco pegmatite, Manitoba, Canada. Am. Mineral. 94, 439–450. doi:10.2138/am.2009.2952

CrossRef Full Text | Google Scholar

Veksler, I. V., Dorfman, A. M., Kamenetsky, M., Dulski, P., and Dingwell, D. B. (2005). Partitioning of lanthanides and Y between immiscible silicate and fluoride melts, fluorite and cryolite and the origin of the lanthanide tetrad effect in igneous rocks. Geochim. Cosmochim. Acta 69, 2847–2860. doi:10.1016/j.gca.2004.08.007

CrossRef Full Text | Google Scholar

Wada, H., Tomita, T., Matsuura, K., Tuchi, K., Ito, M., and Morikiyo, T. (1995). Graphitization of carbonaceous matter during metamorphism with references to carbonate and pelitic rocks of contact and regional metamorphisms, Japan. Contrib. Mineral. Petrol. 118, 217–228. doi:10.1007/BF00306643

CrossRef Full Text | Google Scholar

Wiedenbeck, M., Alle, P., Corfu, F., Griffin, W. L., Meier, M., Oberli, F., et al. (1995). Three natural zircon standards for U-Th-Pb, Lu-Hf, trace-element and REE analyses. Geostand. Newsl. 19 (1), 1–23. doi:10.1111/j.1751-908x.1995.tb00147.x

CrossRef Full Text | Google Scholar

Xia, Q., Li, Y. S., Li, X. F., and Kou, Y. C. (2021). Zircon U-Pb ages, trace elements and Hf isotopes of the mineralized granodiorite dikes in the anfangba gold deposit, West Qinling. Acta Petrol. Sin. 37 (6), 1713–1730. doi:10.18654/1000-0569/2021.06.05

CrossRef Full Text | Google Scholar

Xiao, W. J., Windley, B. F., Yuan, C., Sun, M., Han, C. M., Lin, S. F., et al. (2009). Paleozoic multiple subduction-accretion processes of the southern altaids. Am. J. Sci. 309, 221–270. doi:10.2475/03.2009.02

CrossRef Full Text | Google Scholar

Xiao, W. J., Windley, B., Sun, S., Li, J. L., Huang, B. C., Han, C. M., et al. (2015). A tale of amalgamation of three permo-triassic collage systems in Central Asia: oroclines, sutures, and terminal accretion. Annu. Rev. Earth Planet. Sci. 43, 477–507. doi:10.1146/annurev-earth-060614-105254

CrossRef Full Text | Google Scholar

Xie, L., Wang, Z. J., Wang, R. C., Zhu, J. C., Che, X. D., Gao, J. F., et al. (2018). Mineralogical constraints on the genesis of W–Nb–Ta mineralization in the laiziling granite (Xianghualing district, south China). Ore Geol. Rev. 95, 695–712. doi:10.1016/j.oregeorev.2018.03.021

CrossRef Full Text | Google Scholar

Xu, B. (2010). Research on tectonic-magma ore-controlling regularity and prospecting forecasting of gold metallogenic belt in Kalamaili area, Xinjiang. Beijing. Beijing: China University of Geosciences.

Google Scholar

Xu, B., Lu, Y. M., Gu, X. X., and Zhang, W. Z. (2010). Characteristics of ore-forming fluids and O, H, S isotopic geochemistry of the gold deposits in the Kalamaili area, Xinjiang, China. Earth Sci. Front. 17 (4), 227–240.

Google Scholar

Xu, Q. Q., Zhao, L., Niu, B. G., Zheng, R. G., Yang, Y. Q., and Liu, J. H. (2020). Early Paleozoic arc magmatism in the Kalamaili orogenic belt, Northern Xinjiang, NW China: implications for the tectonic evolution of the East Junggar terrane. J. Asian Earth Sci. 194, 104072. doi:10.1016/j.jseaes.2019.104072

CrossRef Full Text | Google Scholar

Yan, L. Y., Gao, S. X., Chen, Z. G., Jiao, L. X., Sun, L., Liu, Y. F., et al. (2018). Metallogenic characteristics and metallogenic zoning of graphite deposits in China. Geol. China 45 (3), 421–440.

Google Scholar

Yang, G. X., Li, Y. J., Si, G. H., Wu, H. E., Jin, C., and Zhang, Y. Z. (2010). LA-ICP MS U-Pb zircon dating of the Beilekuduke granite in Kalamaili area, East Junggar, Xinjiang, China and its geological implication. Geotecto. Metallog. 34 (1), 133–138.

Google Scholar

Yang, G. X., Li, Y. J., Wu, H. G., Zhong, X., Yang, B. K., Yang, C. X., et al. (2011). Geochronological and geochemical constrains on petrogenesis of the Huangyangshan A-type granite from the East Junggar, Xinjiang, NW China. J. Asian Earth Sci. 40, 722–736. doi:10.1016/j.jseaes.2010.11.008

CrossRef Full Text | Google Scholar

Yang, G. X., Li, Y. J., Wu, H. E., Si, G. H., Jin, Z., and Zhang, Y. Z. (2009). LA-ICP-MS zircon U-Pb dating of the Huangyangshan pluton and its enclaves from Kalamaili area eastern Junggar, Xinjiang, and geological implications. Acta. Petrol. Sin. 25 (12), 3197–3207.

Google Scholar

Yang, W. B., Niu, H. C., Shan, Q., Sun, W. D., Zhang, H., Li, N. B., et al. (2014). Geochemistry of magmatic and hydrothermal zircon from the highly evolved baerzhe alkaline granite: implications for Zr-REE-Nb mineralization. Min. Depo. 49, 451–470. doi:10.1007/s00126-013-0504-1

CrossRef Full Text | Google Scholar

Yin, R., Sun, X. M., Wang, S. W., Wang, R. C., Ran, M. L., Wu, B., et al. (2023). Zircons in NYF-type pegmatites in the emeishan large igneous province, SW China: a record of Nb and REE mineralization processes. Ore Geol. Rev. 162, 105700. doi:10.1016/j.oregeorev.2023.105700

CrossRef Full Text | Google Scholar

Yurimoto, H., Duke, E. F., Opapike, J. J., and Shearer, C. K. (1990). Are discontinuous chondrite-normalized REE patterns in pegmatic granite systems the results of monazite fractionation? Geochim. Cosmochim. Acta 54, 2141–2145. doi:10.1016/0016-7037(90)90277-r

CrossRef Full Text | Google Scholar

Zeng, L. J., Niu, H. C., Bao, Z. W., and Yang, W. B. (2017). Chemical lattice expansion of natural zircon during the magmatic-hydrothermal evolution of A-type granite. Am. Mineral. 102, 655–665. doi:10.2138/am-2017-5840

CrossRef Full Text | Google Scholar

Zhang, Y. Y., Pe-Piper, G., Piper, D. J. W., and Guo, Z. (2013). Early Carboniferous collision of the Kalamaili orogenic belt, North Xinjiang, and its implications: evidence from molasse deposits. Geol. Soc. Am. Bull. 125, 932–944. doi:10.1130/b30779.1

CrossRef Full Text | Google Scholar

Zhang, X. L., Fan, W. J., Li, Z. W., Chen, Z. G., and Chen, J. Y. (2017). The discovery of a superlarge magmatic graphite deposit in Huangyangshan area, Qitai county, Xinjiang. China Geol. 44 (5), 1033–1034.

Google Scholar

Zhao, Z. H., Xiong, X. L., and Han, X. D. (1999). Discussion on the formation mechanism of the tetrad effect of rare earth elements in granites -- a case study of qianlishan and balzhe granite. Sci. Sin. (Terrae) (04), 331–338.

Google Scholar

Zhao, Z. H., Bao, Z. W., and Qiao, Y. L. (2010). A special “M” and “W” compound rare earth element four group effect: taking Shuiquangou alkaline syenite as an example. Chin. Sci. Bull. 55 (15), 1474–1488.

Google Scholar

Zhao, X. F., Zhou, M. F., Su, Z. K., Li, X. C., Chen, W. T., and Li, J. W. (2017). Geology, geochronology, and geochemistry of the dahongshan Fe-Cu-(Au-Ag) deposit, southwest China: implications for the formation of iron oxide copper-gold deposits in intracratonic rift settings. Econ. Geol. 112, 603–628. doi:10.2113/econgeo.112.3.603

CrossRef Full Text | Google Scholar

Zhu, Y. F., and Song, B. (2006). Petrology and SHRIMP chronology of mylonitized tianger granite, Xinjiang: also about the dating on hydrothermal zircon rim in granite. Acta Petrol. Sin. (01), 135–144.

Google Scholar

Keywords: Huangyangshan graphite deposit, zircon, trace elements, U-Pb dating, hydrothermal mineralization, Eastern Junggar

Citation: Sun X, Ren Y, Sun Z, Li J and Huang M (2026) Metallogenic age and hydrothermal evolution of the Huangyangshan graphite deposit in eastern Xinjiang, NW China: constraint from zircon geochronology and geochemistry. Front. Earth Sci. 13:1709937. doi: 10.3389/feart.2025.1709937

Received: 21 September 2025; Accepted: 16 December 2025;
Published: 08 January 2026.

Edited by:

Xiao-Ping Xia, Chinese Academy of Sciences (CAS), China

Reviewed by:

Tao Hong, Sun Yat-sen University, China
Qinglin Xu, Shandong University of Science and Technology, China

Copyright © 2026 Sun, Ren, Sun, Li and Huang. 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.

*Correspondence: Mengjia Huang, aG1qdXV1cEAxNjMuY29t

Disclaimer: 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.