Abstract
The Zhegu region, located in southern Tibet, is positioned within the central and eastern segments of the Tethys Himalayan tectonic belt. In this area, mafic igneous rocks are predominantly intrusion into Jurassic strata, occurring as vein-like bodies. This study presents zircon U-Pb age determinations and whole-rock geochemical analyses of diabase and gabbro samples from the region, aimed at elucidating their petrogenesis and geodynamic background. The zircon U-Pb ages yield crystallization ages of 130.7 ± 1.5 Ma for diabase and 131.6 ± 2.5 Ma for gabbro, both of which are consistent with the crystallization ages of ocean island basalt (OIB)-type mafic rocks in the Comei Large Igneous Province (130–136 Ma). Geochemical data reveal that these mafic rocks are characterized by elevated TiO2, FeOT, and P2O5 contents, alongside relatively low MgO content, indicative of tholeiitic affinities. They exhibit enrichment in light rare earth elements (LREEs) and high field strength elements (HFSEs) such as Nb, Hf and Y, resembling OIB signatures. These rocks show evidence of fractional crystallization without significant crustal contamination. Melting models suggest that the magmas originated from partial melting of a garnet-lherzolite mantle source. The petrogenetic characteristics of these rocks reflect interactions between the Kerguelen mantle plume and the overlying lithospheric mantle.
1 Introduction
The Qinghai-Tibet Plateau is the youngest, highest, and largest plateau on Earth. Its formation resulted from the closure of the New Tethys Ocean, following the collision between the Indian Plate and the Eurasian Plate. This tectonic event led to the development of the Himalayan orogenic belt, which serves as a natural laboratory for studying plate tectonics and the geological evolution of the Tethys region. Situated in the northern part of the Himalayan Mountain system, the Tethys Himalayan Belt extends along the northern edge of the Indian subcontinent, marking the forefront of the Cenozoic collision orogeny between the Indian and Eurasian plates. The magmatic activity and tectonic evolution preceding this orogenic event, particularly during the Mesozoic, have garnered significant interest in recent years as they are key to understanding the boundary conditions and material processes involved in the Cenozoic orogeny of the Himalayan Mountain system (Yin, 2006; Zhu et al., 2009a; 2013; ; Yang et al., 2022; Yu and Liu, 2023).
Previous studies have identified extensive Early Cretaceous mafic igneous rock formations in the eastern region of the Tethys Himalayan belt, particularly in areas such as Comei, Cona, and Longzi. These formations primarily span two time intervals: 140–150 Ma and 130–136 Ma (Zhu et al., 2007, 2008a, b; Xia et al., 2012; Wang et al., 2016; ; ). Based on their geochemical characteristics, these mafic rocks could be classified into three types: Ocean Island Basalt (OIB)-like, Normal Mid-Ocean Ridge Basalt (N-MORB)-like, and Enriched Mid-Ocean Ridge Basalt (E-MORB)-like. These rocks are believed to have formed in an extensional environment characterized by significant lithospheric stretching and thinning along the northern margin of the Gondwana continent (Zhu et al., 2008a; b; Xia et al., 2012; Wang et al., 2016; ; ). The onset of this extensional environment is thought to be related to continental breakup driven by the Kerguelen mantle plume (e.g., Zhu et al., 2009b). These OIB-like, E-MORB-like and N-MORB like igneous rocks are collectively referred to as the Comei Large Igneous Province (LIP), which is hypothesized to represent the early-stage magmatism associated with the Kerguelen mantle plume (e.g., Zhu et al., 2013). There are at least three distinct episodes of Cretaceous mafic magmatism within the Comei LIP area, and only the magmatic event earlier than 120 Ma is considered to be directly linked to Kerguelen mantle plume activity (Wang et al., 2016, 2024). Consequently, distinguishing the igneous rocks associated with the Comei LIP is critical for advancing our understanding of the Kerguelen mantle plume.
Mafic to intermediate dykes are extensivein the Zhegu region, yet detailed investigations into their petrogenesis and potential links to Comei LIP remain lacking. In this study, we present zircon U-Pb dating and whole rock major and trace elements of mafic rocks in the Zhegu area to elucidate their formation age, petrogenesis and tectonic settings. Additionally, we compile data on Early Cretaceous magmatic rocks exhibiting OIB and N/E-MORB like affinities within the Comei LIP. These data provide critical constraints for understanding the relationship between mafic rocks in the Zhegu area and the Comei LIP.
2 Geological background and samples
The Himalayan tectonic belt is subdivided into four tectonic units, arranged from north to south: the Tethys Himalayan, the High Himalayan, the Low Himalayan, and the Sub-Himalayan (Figure 1b). The Tethys Himalayan lies between the Yarlung Zangbo Suture Zone (IYS) and the Southern Tibetan Detachment System (STDS) and is part of the northern section of Greater India in paleogeographic terms. It is regarded as the typical passive continental margin since the Late Triassic period along the northern margin of Greater India (Yu and Wang, 1990).
FIGURE 1
The study area is located in the eastern part of the Tethys Himalayan region and is classified within the Kangmar Lhunze stratigraphic zone of the Gangdise Himalayan stratigraphic region. These strata are primarily controlled by regional faults with a northwest-southeast orientation (Figure 1b). Exposed strata in study area predominantly comprise Mesozoic to Cenozoic sedimentary rocks, particularly those from the Triassic and Jurassic periods. This stratigraphic sequence begins with the Upper Triassic Nieru Formation (T3n), characterized by gray-black, thin to medium-layered silty slate interbedded with gray, medium-layered feldspar quartz sandstone and fine sandstone. This is followed by the Lower to Middle Jurassic Ridang Formation (J1r), which comprises gray, dark gray, and gray-black mudstone, siltstone, sandstone, and shale. The Lower to Middle Jurassic Lure Formation (J1-2l) presents a lithological assemblage of gray mudstone and sandstone, mudstone and limestone, as well as interbedded mudstone and limestone. The sequence continues with the Middle Jurassic Zhela Formation (J2z), primarily composed of mudstone sandstone and silty mudstone. The region is characterized by numerous thrust faults oriented nearly east-west and northwest-southeast, as well as extensional faults oriented nearly north-south. Magmatic activity in this area was intensive, with volcanic rocks, mainly basalts, occurring predominantly within the Jurassic strata. Intrusive rocks are also widespread, often as veins in strata of the Jurassic Lure and Zhela Formation, with lithologies including pyroxenite, gabbro, diabase, and diorite. Notably, gabbro and diabase are exposed within the Zhegu area of the Comei LIP (Figure 1c), and extending in an east-west direction for approximately 20 km in length and 0.5–2 km in width.
The gabbros in this study are gray-black, massive, and predominantly composed of plagioclase (∼45%), clinopyroxene (∼40%), biotite (∼6%), amphibole (∼4%) and quartz (∼4%). Plagioclase occurs as subhedral columnar grains ranging from 0.2 to 1 mm in diameter, with some grains altered to sericite and clay minerals. Clinopyroxene exhibits a subhedral granular morphology, with grain sizes of 0.3–1 mm, and is partially altered to hornblende and chlorite. Biotite is platy, with lengths ranging from 0.3 to 1 mm. Hornblende is subhedral granular and irregular, with 0.3–1 mm in diameter. Quartz, is anhedral and ranges from 0.1 to 0.2 mm in size (Figures 2a,c,e).
FIGURE 2
The diabases are grayish-green, massive, and display an ophitic texture. They are primarily composed of plagioclase (∼45%), clinopyroxene (∼30%), biotite (∼5%), amphibole (∼3%), Fe-Ti oxide (∼4%), and quartz (1%). Plagioclase occurs as euhedral columnar grains, 0.5–1.5 mm in diameter, and is significantly altered to sericite and clay minerals. Clinopyroxene appears as subhedral granular shape, with diameters of 0.3–1 mm. Biotite is euhedral and platy, measuring 0.5–1.1 mm in size. Amphibole is diamond-shaped morphology and is altered to chlorite. Quartz, which, fills voids within clinopyroxene and plagioclase, is anhedral and granular. Additionally, Fe-Ti oxides occur as anhedral granular grains (Figures 2b,d,f).
3 Analytical methods
The gabbro (sample numbers XML001, XML002, XML003, XML004) and diabase (sample numbers ZGL001, ZGL002) were collected for whole-rock major and trace elements analyses, and representative samples (ZGL001, XML001) were selected for zircon U-Pb dating.
Zircons were collected from whole-rock samples using the standard crushing, sieving, heavy liquid, and magnetic separation techniques, and selected under a binocular microscope. These zircons were mounted in an epoxy disk with the zircon standard SL13 from the Australian National University Research School of Earth Sciences and several TEM zircons for reference from the Australian Geological Survey. Subsequently, they were polished down to expose the interior texture. Cathodoluminescence (CL) imaging was obtained at the Institute of Mineral Resources, Chinese Academy of Geological Sciences, using an Oxford MINICL detector.
SHRIMP zircon U-Pb analyses were conducted at the Beijing SHRIMP Center, using the standard zircon TEM (∼417 Ma) for the calibration of instrumental mass fractionation. The U, Th, and Pb contents of the zircons were calibrated using the SL13 zircon, which has an age of ∼572 Ma and a U content of 238 ppm. The instrument and its comprehensive operating principles were the same as previously described (; Sun et al., 2018).
The compositions of major and trace elements were conducted at the Hunan Institute of Mineral Testing and Utilization. Major elements were analyzed by an atomic fluorescence spectrometer (AFS-830A, Jitian Instruments, Beijing, China) and atomic absorption spectrometer (Z-2300, Hitachi Limited, Tokyo, Japan). Trace elements were analyzed by ICP-OES (ICAP6300, Thermo Fisher Scientific, Waltham, MA, USA) and rare earth elements were analyzed by ICP-MS (Thermo X2, Thermo Fisher Scientific, Waltham, MA, USA). The analytical precision for major elements was better than 5%, and that of the rare earth and trace elements was better than 10% ().
4 Results
4.1 Zircon U-Pb ages
Supplementary Table 2 summarizes the U-Pb age data for mafic igneous rocks from the Zhegu area. The zircons are euhedral to subhedral, with lengths ranging from 110 to 180 μm and length/width ratios of 1:1 to 3:1 (Figure 3a). Cathodoluminescence (CL) imagings reveal that many zircons display characteristic oscillatory zoning (Figure 3a). The thorium (Th) and uranium (U) contents of the zircon samples vary significantly between different rock types. For the ZGL001 sample (diabase), Th and U contents range from 157 to 943 ppm and from 140 to 690 ppm, respectively. In contrast, the XML001 sample (gabbro) shows Th and U contents ranging from 578 to 5918 ppm and from 446 to 2052 ppm, respectively. The Th/U ratios of zircons, ranging from 1.16 to 2.19 in the diabase and from 1.30 to 2.98 in the gabbro, are consistently greater than 0.5 which is the characteristics of magmatic zircon (; ). Zircon U-Pb dating results for both diabase and gabbro indicate that all data points lie on or near the Concordia line (Figure 3). The calculated Concordia ages are 130.7 ± 1.5 Ma (MSDW = 1.3, Figure 3b) for the diabase and 131.6 ± 2.5 Ma (MSDW = 3.3, Figure 3c) for the gabbro.
FIGURE 3
4.2 Whole-rock major and trace elements
The contents of major and trace elements of mafic rocks from the Zhegu area of southern Tibet are presented in Supplementary Table 3.
4.2.1 Major elements
The mafic igneous rocks in the Zhegu area exhibit SiO2 contents ranging from 50.2 to 52.6 wt.%, plotting within the basalt region on the TAS diagram (Figure 4a). These rocks are characterized by relatively high TiO2 (2.83–4.98 wt.% with an average of 3.74 wt.%), FeOT (9.07–11.6 wt.% with an average of 9.93 wt.%), and P2O5 content (0.47–1.22 wt.% with an average of 0.63 wt.%). The MgO content is relatively low, ranging from 3.32 to 5.39 wt.%. The Mg# values, which range from 37.5 to 50.2, are significantly lower than those of primary basaltic magma (Mg# = 68–75; Wilson, 1989). These rocks were plotted in the area of the tholeiitic series (Figure 4b). The aluminum saturation indexs (A/CNK) range from 0.63 to 0.77. Additional geochemical parameters further characterize their composition: the Ritter index (σ) ranges from 1.92 to 2.96, the alkalinity rate (AR) ranges from 1.53 to 1.76, and the total alkali content (ALK) ranges from 4.38 to 5.64 wt.%. The consolidation index (SI) varies between 17.2 and 26.8, while the differentiation index (DI) spans from 41.8 to 53.8.
FIGURE 4
4.2.2 Trace elements
The mafic rocks from the Zhegu area in southern Tibet exhibit total rare earth element (REE) contents ranging from 225 to 296 ppm, with a notable enrichment in light rare earth elements (LREE) (LREE/HREE = 2.37–2.73). These rocks also show significant fractionation between LREE and heavy rare earth elements (HREE), with La/Yb ratios ranging from 6.09 to 8.20. The resulting right-skewed REE patterns closely resemble those of OIB (Figure 5a) (Sun and McDonough, 1989). Notably, there are no discernible anomalies in Eu and Ce. The trace element spider diagram reveals enrichment in Nb and Hf, coupled with a depletion of Y, and the distribution patterns closely resemble those of OIB (Figure 5b). The ratios of Ti/Y, Th/Ta, Th/Yb, Ce/Zr, and Zr/Y are 377–642, 0.64–1.75, 0.79–1.55, 0.21–0.31, and 5.53–9.30, respectively, aligning closely with those observed in OIB basalts (Supplementary Table 3; ; Zhu et al., 2005).
FIGURE 5
5 Discussion
5.1 Formation ages
The crystallization ages of diabase and gabbro in this study are 130.7 ± 1.5 Ma and 131.6 ± 2.5 Ma, respectively. These crystallization ages are comparable to those of OIB-like igneous rocks in the southern and western regions of the Comei LIP, where the peak magmatic activity occurred ∼132 Ma. Furthermore, these ages align with the early-stage magma activity associated with the Kerguelen mantle plume event (Zhu et al., 2009a; ). The close temporal coincidence suggests a genetic relationship between the mafic igneous rocks in the Zhegu area and early-stage magma activity of the Kerguelen mantle plume.
5.2 Petrogenesis
5.2.1 Crustal contamination
We propose that the mafic rocks in the Zhegu area have undergone limited crustal contamination. The La/Sm ratio is generally considered stable, but it could increase significantly (>5.00) during crustal contamination (; Zhang et al., 2004). The La/Sm ratio of the Zhegu mafic rocks ranges from 2.94 to 3.38, and in the La/Sm vs. La/Nb diagram (Figure 6a), the samples do not exhibit a trend consistent with crustal contamination. Similarly, the Ce vs. Nb/Th diagram (Figure 6b) shows no linear trend, further indicating minimal crustal influence (; Rudnick and Gao, 2003). In the (Th/Ta)PM vs. (La/Nb)PM diagram (Figure 6c), the Zhegu mafic rocks are distinctly separated from the upper, middle, and lower crust, resembling the characteristics of Hawaiian Island basalt, 90° E ridge basalt, and Kerguelen OIB. The (La/Nb)PM vs. (Th/Nb)PM diagram (shows that the samples in this study plot in the region of basalts unaffected by crustal contamination (Figure 6d; ). Additionally, crustally contaminated magma typically shows negative anomalies in Nb, Ta, and Ti (Thompson et al., 1984; Wilson, 1989), which are not apparent in the Zhegu mafic rocks. The Nb content in the upper, middle, and lower crust (12 ppm, 10 ppm, 5 ppm, respectively) is significantly lower than that in these mafic rocks (30.4–42.8 ppm), further suggesting minimal crustal contamination.
FIGURE 6
Overall, the geochemical characters of the Zhegu mafic rocks strongly indicate that they have not been significantly contaminated by crustal material.
5.2.2 Magma mechanism
Partial melting and fractional crystallization are two primary mechanisms driving magmatism and its evolution. These processes can be differentiated using the La-La/Sm discriminant diagram (Figure 7a), where the samples exhibit a horizontal linear distribution, suggesting fractional crystallization as the dominant process.
FIGURE 7
Additionally, the Mg# of the samples ranges from 37.5 to 50.2, which is significantly lower than the Mg# of primary basaltic magma (68–75) (Wilson, 1989), indicating a high degree of magma differentiation. The concentrations of compatible elements, such as Ni and Cr, are notably low, ranging from 5.6 to 63.6 ppm and 8.0–100 ppm, respectively. The Harker diagrams show a positive correlation between Mg# and Ni, Cr (Figures 8a,b), indicating the significant fractional crystallization of olivine and clinopyroxene (; ). As Mg# decreases, the CaO and Al2O3 contents either increase or remain stable (Figures 8f,h), with no evidence of a negative Eu anomaly in the REE patterns (Figure 5). It is worth noting that these mafic rocks contain significant amounts of plagioclase (Figure 2), suggesting that the plagioclase did not separate from the magma following its crystallization. Furthermore, the correlations between Mg# and Zr and P2O5 are negative and not obvious respectively (Figures 8c,g), implies that the crystallization of minerals such as apatite and zircon was not significant (Wang et al., 2018).
FIGURE 8
Beyond fractional crystallization, partial melting in the mantle source also plays a crucial role in shaping the trace element characteristics of these samples. The La/Sm and Sm/Yb ratios are particularly informative for inferring the source characteristics and the degree of partial melting of mantle-derived samples (). Melts derived from garnet- or spinel-bearing peridotite are expected to be enriched in LREE, with relatively depleted MREE. Heavy rare earth elements (HREE), such as Yb, are compatible with garnet but incompatible with other mantle phases (; ). As a result, HREE enrichment would occur if garnet were present in the source. In the spinel-bearing peridotite stable field, an increase in the degree of partial melting leads to a decrease in the La/Sm ratio of the melt, while the Sm/Yb ratio remains relatively constant. In contrast, within the garnet-bearing peridotite stability field, increasing partial melting results in a significant decrease in both the La/Sm and Sm/Yb ratios. Figure 7b illustrates this concept, comparing potential mantle sources, such as depleted MORB mantle (DMM), primary mantle (PM), and Kerguelen plume-derived magmas (sample s749C-16R-16 from Site 749, ). The simulation calculation of variations in REE ratios, as depicted in the La/Sm-Sm/Yb diagram (Figure 7b), suggests that samples display characteristics typical of ocean island basalts (OIB). These samples are plotted within the 5%–10% partial melting range of garnet peridotite in the primitive mantle (PM), implying that the mafic rocks in the Zhegu area primarily formed through low-degree partial melting of garnet-bearing peridotite.
5.2.3 Magmatic source
This study investigates the ratios of incompatible elements (Supplementary Table 3; Figure 9c) to characterize the trace element composition of the source region. The results suggest that the incompatible element ratios of the mafic rocks in the Zhegu area closely resemble those of the Sangxiu Formation basalt, Hawaiian basalt, Kerguelen OIB, Emeishan high-Ti basalt, and 90ºE ridge basalt. This resemblance suggests that the source of mafic rocks in the Zhegu region is analogous to that of hotspot or mantle plume magmas. In the OIB source discrimination diagram, all samples fall within or near the OIB field (Figure 9). Coupled with the similar geochemical, REE and trace element distribution patterns (Figure 5), this indicates that the mafic rocks in the Zhegu area likely originate from a mantle source akin to OIB, potentially the product of hotspot or mantle plume activity (Zhu et al., 2005; Yang et al., 2015; ). Previous research indicates that mafic dykes originating from the depleted mantle should have a Zr/Nb ratio greater than 18 (Zheng et al., 2020). However, the Zr/Nb ratios of the Zhegu mafic rocks range from 6.95 to 13.5, all below 18. Additionally, in the Ta/Yb-Th/Yb diagram (Figure 9c), samples in the Zhegu area show an enriched mantle characteristic. In addition, the mafic rocks in the Zhegu area also have been influenced by lithospheric mantle materials. Firstly, the enrichment of LILE (such as Rb and Ba) in the mafic rocks of the Zhegu area suggests the involvement of the lithospheric mantle (). Secondly, the incompatible element ratios of the Zhegu mafic rocks are similar to those of the Sangxiu Formation basalts in Supplementary Table 3, which previous studies have attributed to the influence of the lithospheric mantle (Zhu et al., 2005).
FIGURE 9
In summary, we propose that the source of the Zhegu mafic rocks is similar to that of OIB-type mafic dyke swarms within the Tethys Himalayan tectonic belt of southern Tibet (; ; ; ). These rocks represent products of the interaction between enriched hotspots or mantle plume sources and the lithospheric mantle.
5.3 Tectonic setting
Mafic dyke swarms typically form as a result of the rapid intrusion of mafic magma from deep sources into the near-surface crust under conditions of crustal extension. These swarms can provide valuable insights into the geodynamic settings of the region.
The OIB are a prominent type of igneous rocks found within ocean basins and are considered a representative product of magmatic processes occurring in these regions. Similarly, in continental rifts, igneous rocks with geochemical characteristics similar to OIB can also develop. In addition to intraplate tectonic settings, OIB-like magmatic rocks can form in back-arc extensional rifts associated with subduction-related arc magmatism (; , ). However, since the Early Cretaceous OIB-like basalts in the Tethys Himalayan Belt are not contemporaneous with arc magmatism (Yang et al., 2022), we tend to exclude this possibility.
In the Zr-Zr/Y and Ta/Hf-Th/Hf diagrams (Figures 10a,b), all samples fall within the intraplate basalt region, aligning with previous studies (Zhu et al., 2004, 2005; ; Ren et al., 2015), suggesting that these rocks formed in intraplate tectonic environments. Zhu et al. (2009a), Zhu et al. (2013) proposed that extensive magmatic activity during the Early Cretaceous in the Comei LIP contributed to the formation of OIB-type mafic igneous rocks, which are distributed across regions such as Nagarze, Lhozhag, Cona, and Comei. Further research has indicated that the mantle plume head of the Comei LIP is compositionally similar to that of the Kerguelen mantle plume, suggesting a common magma source (). Recent chronology and geochemical studies also support the hypothesis that the Early Cretaceous OIB-type basalts in the Tethys Himalayan Belt are closely related to early-stage activity of the Kerguelen Mantle plume (Zhu et al., 2005, 2008b; 2013; ; Xia et al., 2012; ; Wang et al., 2016; Wei et al., 2017; Zhou et al., 2017; Tian et al., 2019). The OIB-type mafic rocks of the Comei LIP are characterized by high TiO2 (average of ∼3.0 wt.%) and P2O5 (average of ∼0.4 wt.%) with no Nb-Ta negative anomalies (Xia et al., 2012). This study finds that the mafic rocks in the Zhegu area also exhibit TiO2 (average of ∼3.74 wt.%) and P2O5 (average of ∼0.63 wt.%), and similarly lack negative Nb-Ta anomalies, resembling the characteristics of the mafic rocks in the Comei LIP. Additionally, the formation ages of diabase and gabbro (130.7 ± 1.5 Ma and 131.6 ± 2.5 Ma, respectively) in this study coincide with the age range of OIB-type rocks in the Comei LIP (130–136 Ma) (Zhu et al., 2009b), which are products of the early-stage activity of the Kerguelen mantle plume.
FIGURE 10
In conclusion, the mafic rocks in the Zhegu area are part of the Comei Large Igneous Province and likely represent products of early Kerguelen mantle plume activity.
6 Conclusion
(1) Zircon U-Pb dating indicates that diabase and gabbro crystallized during the Early Cretaceous period, with ages of 130.7 ± 1.5 Ma and 131.6 ± 2.5 Ma, respectively, which are comparable to the formation age of OIB-type mafic rocks found in the Comei Large Igneous Province. The trace elements distribution patterns of mafic igneous rocks in Zhegu area are similar with that of OIB.
(2) The mafic rocks in the Zhegu area experiences fractional crystallization. And the assimilation of crust materials was insignificant. They are products of the interaction between enriched hotspots/mantle plumes and the lithospheric mantle, which was associated with the early activity of the Kerguelen mantle plume.
Statements
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 authors.
Author contributions
MC: Data curation, Funding acquisition, Investigation, Project administration, Writing – original draft, Writing – review and editing. SS: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review and editing. YL: Formal Analysis, Funding acquisition, Methodology, Resources, Software, Writing – original draft. YM: Investigation, Methodology, Validation, Visualization, Writing – original draft. YT: Data curation, Investigation, Software, Validation, Writing – original draft. XL: Formal Analysis, Methodology, Resources, Visualization, Writing – original draft. MZ: Methodology, Resources, Software, Writing – original draft. XH: Formal Analysis, Funding acquisition, Project administration, Validation, Writing – original draft. TW: Investigation, Project administration, Software, Validation, Writing – original draft. HZ: Formal Analysis, Methodology, Validation, Visualization, Writing – original draft. KX: Formal Analysis, Resources, Validation, Visualization, Writing – original draft. CC: Formal Analysis, Investigation, Methodology, Software, Writing – original draft. JZ: Data curation, Methodology, Supervision, Writing – review and editing. WG: Funding acquisition, Project administration, Supervision, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the National Key R&D Program of China (2023YFF0807103), the National Natural Science Foundation of China (42373049), the China Geological Survey Project (1212011220659, DD20230527, ZD20220309, DD20243079), the Taishan Scholar Program of Shandong (tspd20230609) and Laoshan Laboratory (LSKJ202204100).
Conflict of interest
The 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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2025.1513583/full#supplementary-material
References
1
AldanmazE.PearceJ. A.ThirlwallM. F.MitchellJ. G. (2000). Petrogenetic evolution of late Cenozoic, post-collision volcanism in western Anatolia, Turkey. J. Volcanol. Geotherm. Res.102 (1), 67–95. 10.1016/S0377-0273(00)00182-7
2
BaksiA. K. (1995). Petrogenesis and timing of volcanism in the Rajmahal flood basalt province, northeastern India. Chem. Geol.121 (1), 73–90. 10.1016/0009-2541(94)00124-Q
3
ChenS. S.FanW. M.ShiR. D.LiuX. H.ZhouX. J. (2018). 118–115 Ma magmatism in the Tethyan Himalaya igneous province: Constraints on early cretaceous rifting of the northern margin of greater India. Earth Planet. Sci. Lett.491, 21–33. 10.1016/j.epsl.2018.03.034
4
ChengM.HuX.TangY.DengZ.MinY.ChenS.et al (2024). Pan-African and early paleozoic orogenic events in southern Tibet: Evidence from geochronology and geochemistry of the Kangbuzhenri Gneissic granite in the Zhegu area. Minerals14 (8), 845. 10.3390/min14080845
5
ChengM.LouY. L.TangR.LiaoJ.ChenS. M.LiY.et al (2022). Zircon U-Pb dating, geochemistry characteristics and tectonic implications of mafic dykes in the Chaiwa area, Southern Tibet. Acta Petrologica Mineralogica41 (3), 504–518. 10.3969/j.issn.1000-6524.2022.03.002
6
CoffinM.EldholmO. (1992). Volcanism and continental break-up: A global compilation of large igneous provinces. Geol. Soc.68. 17–30. 10.1144/GSL.SP.1992.068.01.02
7
CoffinM. F.PringleM. S.DuncanR. A.GladczenkoT. P.StoreyM.MüllerR. D.et al (2002). Kerguelen hotspot magma output since 130 Ma. J. Petrol.43 (7), 1121–1137. 10.1093/petrology/43.7.1121
8
CondieK. C.FreyB. A.KerrichR. (2002). The 1.75-Ga Iron King Volcanics in west-central Arizona: a remnant of an accreted oceanic plateau derived from a mantle plume with a deep depleted component. Lithos64 (1), 49–62. 10.1016/S0024-4937(02)00158-5
9
DingF.GaoJ. G.XuK. Z. (2020). Geochemistry, geochronology and geological significances of the basic dykes in Rongbu area, southern Tibet. Acta Petrol. Sin.36 (2), 391–408. 10.18654/1000-0569/2020.02.04
10
FreyF. A.CoffinM. F.WallaceP. J.WeisD.ZhaoX.WiseS. W.et al (2000). Origin and evolution of a submarine large igneous province: the Kerguelen Plateau and Broken Ridge, southern Indian Ocean. Earth Planet. Sci. Lett.176 (1), 73–89. 10.1016/S0012-821X(99)00315-5
11
FreyF. A.GreenD. H.RoyS. D. (1978). Integrated models of basalt petrogenesis: A study of quartz Tholeiites to olivine melilitites from south Eastern Australia utilizing geochemical and experimental petrological data. J. Petrology19 (3), 463–513. 10.1093/petrology/19.3.463
12
FreyF. A.McNaughtonN. J.NelsonD. R.deLaeterJ. R.DuncanR. A. (1996). Petrogenesis of the Bunbury Basalt, western Australia: interaction between the Kerguelen plume and Gondwana lithosphere?Earth Planet. Sci. Lett.144 (1), 163–183. 10.1016/0012-821X(96)00150-1
13
FreyF. A.WeisD.BorisovaA. Y.XuG. (2002). Involvement of continental crust in the formation of the cretaceous Kerguelen plateau: New perspectives from ODP leg 120 sites. J. Petrology43 (7), 1207–1239. 10.1093/petrology/43.7.1207
14
GaoX.YuS.PengY.LvP.WangM.LiuY.et al (2021). Insights into OIB-like magmatism contemporaneous with oceanic subduction: Petrogenetic constraints on the Kendelong metagabbro in the North Qaidam. Lithos392-393, 106130. 10.1016/j.lithos.2021.106130
15
GaoZ.ZhangH. F.YangH.PanF. B.LuoB. J.GuoL.et al (2018). Back-arc basin development: Constraints on geochronology and geochemistry of arc-like and OIB-like basalts in the Central Qilian block (Northwest China). Lithos310-311, 255–268. 10.1016/j.lithos.2018.04.002
16
HessP. C. (1992). Phase equilibria constraints on the origin of ocean floor basalts, 71. Geophysical Monograph.
17
HoskinP. W. O.BlackL. P. (2000). Metamorphic zircon formation by solid-state recrystallization of protolith igneous zircon. J. Metamorph. Geol.18 (4), 423–439. 10.1046/j.1525-1314.2000.00266.x
18
HuangY.LiangW.ZhangL. K.LiG. M.HuangC. M.XiaoX. B.et al (2018). The initial break-up between Tethyan-Himalaya and Indian Terrane: evidences from Late Cretaceous OIB-type basalt in southern Tibet. Earth Sci.43 (08), 2651–2663.
19
IrvingA. J.FreyF. A. (1984). Trace element abundances in megacrysts and their host basalts: Constraints on partition coefficients and megacryst genesis. Geochimica Cosmochimica Acta48 (6), 1201–1221. 10.1016/0016-7037(84)90056-5
20
JiangS. H.NieF. J.HuP.LiuY. (2010). An important spreading event of the neo-tethys ocean during the late Jurassic and early cretaceous: Evidence from Zircon U-Pb SHRIMP dating on diabase in Nagarzê, southern Tibet. Acta Geol.80, 522–527. 10.1111/j.1755-6724.2006.tb00272.x
21
JiangS. H.NieF. J.HuP.LiuY.LaiX. R. (2007). Geochemical characteristics of the mafic dyke swarms in South Tibet. Acta Geol. Sin.81 (01), 60–71. 10.3321/j.issn:0001-5717.2007.01.008
22
KentW.SaundersA. D.KemptonP. D.GhoseN. C. (1997). Rajmahal basalts, eastern India: Mantle sources and melt distribution at a volcanic rifted margin, 100. Washington DC American: Geophysical Union Geophysical Monograph.
23
KinzlerR. J. (1997). Melting of mantle peridotite at pressures approaching the spinel to garnet transition: Application to mid-ocean ridge basalt petrogenesis. J. Geophys. Res. Solid Earth102 (B1), 853–874. 10.1029/96JB00988
24
LassiterJ. C.DePaoloD. J. (1997). Plume/Lithosphere Interaction in the generation of continental and oceanic flood basalts: chemical and isotopic constraints. Large Igneous.
25
LiS. G. (1993). Ba-Nb-Th-La diagrams used to identify tectonic enviroments of ophiolite. Acta Petrol. Sin. (02), 146–157.
26
LiY.DingF.YangL.XieX.XuK.LiQ.et al (2019). Geochemical characteristics and geological significance of basalts in the Sangxiu formation, Rongbu area, southern Tibet. Arabian J. Geosciences12 (24), 757. 10.1007/s12517-019-4880-4
27
LiangF. H.XuZ. Q.BaD. Z.XuX. Z.LiuF.XiongF. H. (2011). Tectonic occurrence and emplacement mechanism of ophiolites from Luobusha-Zedang, Tibet. Acta Petrol. Sin.27 (11), 3255–3268.
28
LinB.TangJ.ChenY.BakerM.SongY.YangH.et al (2019). Geology and geochronology of Naruo large porphyry-breccia Cu deposit in the Duolong district, Tibet. Gondwana Res.66, 168–182. 10.1016/j.gr.2018.07.009
29
LiuD. Y.JianP.ZhangQ.ZhangF. Q.ShiY. R.ShiG. H.et al (2003). SHRIMP dating of adakites in the tulingkai ophiolite, Innr Mongolia: evidence for the early paleozoic subduction. Acta Geol. Sin. (03), 317–327. 10.3321/j.issn:0001-5717.2003.03.004
30
LiuF.YangJ. S.ChenS. Y.LiZ. L.LianD. Y.ZhouW. D.et al (2015). Geochemistry and Sr-Nd-Pb isotopic composition of mafic rocks in the western part of Yarlung Zangbo suture zone: Evidence for intra-oceanic supra-subduction within the Neo-Tethys. Geol. China36 (04), 441–454. 10.3969/j.issn.1000-3657.2013.03.007
31
LiuH.WangY.CawoodP. A.GuoX. (2017). Episodic slab rollback and back-arc extension in the Yunnan-Burma region: Insights from Cretaceous Nb-enriched and oceanic-island basalt–like mafic rocks. GSA Bull.129 (5-6), 698–714. 10.1130/B31604.1
32
LiuZ.ZhouQ.LaiY.QingC.LiY.WuJ.et al (2015). Petrogenesis of the early cretaceous Laguila bimodal intrusive rocks from the Tethyan Himalaya: Implications for the break-up of eastern Gondwana. Lithos236-237, 190–202. 10.1016/j.lithos.2015.09.006
33
MahoneyJ. J.JonesW. B.FreyF. A.SaltersV. J. M.PyleD. G.DaviesH. L. (1995). Geochemical characteristics of lavas from Broken ridge, the naturaliste plateau and southernmost Kerguelen plateau: Cretaceous plateau volcanism in the southeast Indian ocean. Chem. Geol.120 (3), 315–345. 10.1016/0009-2541(94)00144-W
34
McDonoughW. F. (1990). Constraints on the composition of the continental lithospheric mantle. Earth Planet. Sci. Lett.101 (1), 1–18. 10.1016/0012-821X(90)90119-I
35
McKenzieD.O'NionsR. K. (1991). Partial melt distributions from Inversion of rare earth element concentrations. J. Petrology32 (5), 1021–1091. 10.1093/petrology/32.5.1021
36
MiddlemostE. A. K. (1994). Naming materials in the magma/igneous rock system. Earth-Science Rev.37 (3), 215–224. 10.1016/0012-8252(94)90029-9
37
MiyashiroA. (1974). Volcanic rock series in island arcs and active continental margins. Am. J. Sci.274 (4), 321–355. 10.2475/ajs.274.4.321
38
MysenB. O. (1979). Trace-element partitioning between garnet peridotite minerals and water-rich vapor; experimental data from 5 to 30 kbar. Am. Mineralogist64 (3-4), 274–287.
39
NealC. R.MahoneyJ. J.ChazeyW. J. (2002). Mantle sources and the highly variable role of continental lithosphere in basalt petrogenesis of the Kerguelen plateau and Broken ridge LIP: Results from ODP leg 183. J. Petrology43 (7), 1177–1205. 10.1093/petrology/43.7.1177
40
PanG. T.WangL. Q.ZhuD. C. (2004). Thoughts on some important scientific problems in regional geological survey of the Qinghai-Tibet Plateau. Geol. Bull. China23 (01), 12–19.
41
PearceJ. A. (1982). Trace element characteristics of lavas from destructive plate boundaries.
42
PearceJ. A. (2008). Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust. Lithos100 (1), 14–48. 10.1016/j.lithos.2007.06.016
43
PearceJ. A.NorryM. J. (1979). Petrogenetic implications of Ti, Zr, Y, and Nb variations in volcanic rocks. Contributions Mineralogy Petrology69 (1), 33–47. 10.1007/BF00375192
44
QiuB. B. (2011). The westward extension of Comei fragmented large igneous province in southem Tibet and its implications. Acta Petrol. Sin.26 (7), 2207–2216.
45
QiuB. B.ZhuD. C.ZhaoZ. D.WangL. Q. (2010). The westward extension of Comei fragmented large igneous province in southern Tibet and its implications. Acta petrol. Sin.26 (07), 2207–2216.
46
RenC.LiuS.ZhuL. D.PanJ. T.GaoW. X. (2014). Geochemistry and Zircon SHRIMP U-Pb dating and their tectonic significance for intermediate-basic dyke in the Gudui region, South Tibet. J. Sichuan Geol.34 (04), 496–500.
47
RenC.MaF. Z.ZhuZ. H.ZhangJ. G.HuangR. C.ZhangY. (2015). U-Pb SHRIMP zircon ages of the mafic-ultramafic rocks from Chigu Co of South Tibet and their geological significances. Geol. China42 (4), 881–890.
48
RudnickR.GaoS. (2003). Composition of the continental crust. Treatise geochem.3, 1–64. 10.1016/B0-08-043751-6/03016-4
49
ShawD. M. (1970). Trace element fractionation during anatexis. Geochimica Cosmochimica Acta34 (2), 237–243. 10.1016/0016-7037(70)90009-8
50
SunS. J.IrelandT. R.ZhangL.ZhangR.ZhangC.SunW. (2018). Palaeoarchaean materials in the Tibetan Plateau indicated by zircon. Int. Geol. Rev.60 (8), 1061–1072. 10.1080/00206814.2017.1367967
51
SunS. S.McDonoughW. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes, 42. London: Special Publications, 313–345.
52
TaylorS. R.McLennanS. M. (1985). The continental crust: Its composition and evolution. Oxford: Blackwell Scientific Publications, 312pp.
53
ThompsonR. N.MorrisonM. A.DickinA. P.HendryG. L. (1984). An assessment of the relative roles of crust and mantle in magma genesis: an elemental approach, 310, 549–590.
54
TianY.GongJ.ChenH.GuoL.XuQ.ChenL.et al (2019). Early Cretaceous bimodal magmatism in the eastern Tethyan Himalayas, Tibet: Indicative of records on precursory continental rifting and initial breakup of eastern Gondwana. Lithos324-325, 699–715. 10.1016/j.lithos.2018.12.001
55
TongJ. S.LiuJ.ZhongH. M.XiaJ.LiY. H. (2007). Zircon U- Pb dating and geochemistry of mafic dike swarms in the Lhozag area, southern Tibet, China, and their tectonic implications. Geol. Bull. China26 (12), 1654–1664.
56
WalterM. J. (1998). Melting of garnet peridotite and the origin of komatiite and depleted lithosphere. J. Petrology39, 29–60. 10.1093/petrology/39.1.29
57
WangX. C.XiaB.LiuW. L.ZhongY.HuX. C.GuanY.et al (2018). Geochronology, geochemistry and petrogenesis of the Pungco Ophiolite, Tibet. Geotect. Metallogenia42 (03), 550–569. 10.16539/j.ddgzyckx.2018.03.010
58
WangY.ZengL.HouK.GaoL. E.WangQ.ZhaoL.et al (2022). Mantle source components and magmatic evolution for the Comei large igneous province: Evidence from the early cretaceous Niangzhong mafic magmatism in Tethyan Himalaya. J. Earth Sci.33 (1), 133–149. 10.1007/s12583-021-1464-5
59
WangY. L.ZhangC. J.XiuS. Z. (2001). Th/Hf-Ta/Hf idenfication of tectonic setting of basalts. Acta Petrol. Sin.17, 413–421.
60
WangY. Y.GaoL. E.ZengL. S.ChenF. K.HouK. J.WangQ.et al (2016). Multiple phases of cretaceous mafic magmatism in the Gyangze-Kangma area, Tethyan Himalaya, southern Tibet. Acta Petrol. Sin.32 (12), 3572–3596.
61
WangY. Y.ZengL. S.GaoL. E.ZhaoL. H.YanL. L. (2024). Early arrival of the Reunion plume at the base of the Himalaya?Earth Planet. Sci. Lett.639, 118756. 10.1016/j.epsl.2024.118756
62
WeiY.LiangW.ShangY.ZhangB.PanW. (2017). Petrogenesis and tectonic implications of ∼130Ma diabase dikes in the western Tethyan Himalaya (western Tibet). J. Asian Earth Sci.143, 236–248. 10.1016/j.jseaes.2017.04.008
63
WilsonM. (1989). Igneous petrogenesis. London: Unwin Hyman.
64
WorkmanR. K.HartS. R. (2005). Major and trace element composition of the depleted MORB mantle (DMM). Earth Planet. Sci. Lett.231 (1), 53–72. 10.1016/j.epsl.2004.12.005
65
XiaY.ZhuD. C.ZhaoZ. D.WangQ.YuanS. H.ChenY.et al (2012). Whole-rock geochemistry and zircon Hf isotope of the OIB-type mafic rocks from the Comei Large Igneous Province in southeastern Tibet. Acta Petrol. Sin.28 (5), 1588–1602.
66
YangC. Y.LiY. B.ZhangJ. S.ZhangG. (2022). Geochronology, petrogenesis and tectonic significance of the Kalong gabbro from Nagarze, South Tibet. Chin. J. Geol.57 (02), 399–414.
67
YangG.LiY.XiaoW.TongL. (2015). OIB-type rocks within West Junggar ophiolitic mélanges: Evidence for the accretion of seamounts. Earth-Science Rev.150, 477–496. 10.1016/j.earscirev.2015.09.002
68
YinA. (2006). Cenozoic tectonic evolution of the Himalayan orogen as constrained by along-strike variation of structural geometry, exhumation history, and foreland sedimentation. Earth-Science Rev.76 (1), 1–131. 10.1016/j.earscirev.2005.05.004
69
YuG. M.WangC. S. (1990). Tibet titis Deposition geology. Beijing: Beijing Geological Publishing House.
70
YuY. Z.LiuZ. C. (2023). Origin and genetic mechanism of the early cretaceous metabasite from the Ramba area, Tethyan Himalaya geological. J. China Univ.29 (04), 527–542.
71
ZhangY. X.XieC. M.YuY. P.ZhangH. Y.DongY. C. (2018). The Early Jurassic subduction of Neo-Tethyan oceanic slab: Constraints from zircon U-Pb age and Hf isotopic compositions of Sumdo high-Mg diorite. Geol. Bull. China37 (08), 1387–1399.
72
ZhangZ. C.WangF. S.HaoY. L.MahoneyJ. J. (2004). Geochemistry of the Picrites and associated basalts from the Emeishan large igneous basalt province and constraints on their source region. Acta Geol. Sin. (02), 171–180. 10.3321/j.issn:0001-5717.2004.02.005
73
ZhengJ. L.LiuT.XuL. M.LiangZ. K.GuoX. Y.WangD. K.et al (2020). Geochemical and chronologic characteristics of ultramafic rocks in Gaxian ophiolitic mélange of the Da Hinggan Mountains and their geological significance. Geol. Bull. China39 (4), 480–490.
74
ZhongH. M.TongJ. S.XiaJ.LuR. K.QiuJ. Q. (2005). Characteristics and tectonic setting of volcanic rocks of the Sangxiu Formation in the southern part of Yamzho Yumco, southern Tibet. Geol. Bull. China24 (1), 72–79.
75
ZhouQ.LiuZ.LaiY.WangG. C.LiaoZ. W.LiY. X.et al (2017). Petrogenesis of mafic and felsic rocks from the Comei large igneous province, South Tibet: Implications for the initial activity of the Kerguelen plume. GSA Bull.130 (5-6), 811–824. 10.1130/B31653.1
76
ZhuD. C.ChungS. L.MoX. X.ZhaoZ. D.NiuY.SongB.et al (2009a). The 132 Ma Comei-Bunbury large igneous province: Remnants identified in present-day southeastern Tibet and southwestern Australia. Geology37 (7), 583–586. 10.1130/G30001A.1
77
ZhuD. C.MoX. X.PanG. T.ZhaoZ. D.DongG. C.ShiY. R.et al (2008b). Petrogenesis of the earliest Early Cretaceous mafic rocks from the Cona area of the eastern Tethyan Himalaya in south Tibet: Interaction between the incubating Kerguelen plume and the eastern Greater India lithosphere?Lithos100 (1), 147–173. 10.1016/j.lithos.2007.06.024
78
ZhuD. C.MoX. X.ZhaoZ. D.NiuY. L.ChungS. L. (2008a). Whole-rock elemental and zircon Hf isotopic geochemistry of mafic and ultramafic rocks from the Early Cretaceous Comei large igneous province in SE Tibet: constraints on mantle source characteristics and petrogenesis. Himal. J. Sci.5, 178–180. 10.3126/hjs.v5i7.1350
79
ZhuD. C.MoX. X.ZhaoZ. D.NiuY. L.PanG. T.WangL. T.et al (2009b). Permian and early cretaceous tectonomagmatism in southern Tibet and Tethyan evolution: New perspective. Earth Sci. Front.16 (2), 001–020.
80
ZhuD. C.PanG. T.MoX. X.LiaoZ. L.JiangX. S.WangL. Q.et al (2005). Geochemistry and petrogenesis of the Sangxiu Formation basalts in the central segment of Tethyan Himalaya. Geochim. (01), 7–19.
81
ZhuD. C.PanG. T.MoX. X.LiaoZ. L.JiangX. S.WangL. Q.et al (2007). Petrogenesis of volcanic rocks in the Sangxiu Formation, central segment of Tethyan Himalaya: A probable example of plume–lithosphere interaction. J. Asian Earth Sci.29 (2), 320–335. 10.1016/j.jseaes.2005.12.004
82
ZhuD. C.WangL. Q.PanG. T.MoX. X.LiaoZ. L.JiangX. S.et al (2004). Discrimination of OIB-type magma and its significances of basalts from middle Jurassic Zhela Formation in the central belt of Tethyan Himalayas, South Tibet. Geol. Sci. Technol. Intell.23 (03), 15–24. 10.3969/j.issn.1000-7849.2004.03.004
83
ZhuD. C.XiaY.QiuB. B.WangQ.ZhaoZ. D. (2013). Why do we need to propose the Early Cretaceous Comei large igneous province in southeastern Tibet?Acta Petrol. Sin.29 (11), 3659–3670.
Summary
Keywords
zircon U-Pb age, whole-rock geochemistry, OIB-like type mafic rocks, Kerguelen mantle plume, southern Tibet
Citation
Cheng M, Sun S, Lou Y, Min Y, Tang Y, Li X, Zhang M, Hu X, Wan T, Zou H, Xu K, Chen C, Zhang J and Guo W (2025) Geochronology and geochemistry of mafic igneous rocks in the Zhegu area of southern Tibet. Front. Earth Sci. 13:1513583. doi: 10.3389/feart.2025.1513583
Received
18 October 2024
Accepted
24 June 2025
Published
11 July 2025
Volume
13 - 2025
Edited by
Li Tian, China University of Geosciences Wuhan, China
Reviewed by
Magdalena Matusiak-Malek, University of Wrocław, Poland
Debajyoti Paul, Indian Institute of Technology Kanpur, India
Updates
Copyright
© 2025 Cheng, Sun, Lou, Min, Tang, Li, Zhang, Hu, Wan, Zou, Xu, Chen, Zhang and Guo.
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: Saijun Sun, sunsaijun06@163.com; Junjie Zhang, zhangjunjie@qdio.ac.cn; Wei Guo, gw02108106@126.com
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.