Mafic Intrusions in the Cuona Area, Eastern Tethyan Himalaya: Early Kerguelen Mantle Plume Activity and East Gondwana Rifting

Abstract

The widespread occurrence of Mesozoic ocean island basalt (OIB)-like igneous rocks in the Southern Tibetan Himalayan Belt provides important constraints on the rifting of East Gondwana. This study undertook a petrological, geochronological, and geochemical investigation of mafic intrusive rocks in the Cuona area of the eastern Tethyan Himalayan Belt. The mafic intrusions have OIB-type geochemical signatures, including diabase porphyrite, gabbro, and diabase. Zircon U–Pb dating indicates that the diabase porphyrite formed at 135.0 ± 1.6 Ma. The diabase porphyrite and gabbro are enriched in high-field-strength elements (Nb and Ti) and large-ion lithophile elements (Sr and Pb) and experienced negligible lithospheric mantle or crustal contamination. The diabase is enriched in large-ion lithophile elements (LILEs, e.g., La and Ce) and depleted in high-field-strength elements (HFSEs, e.g., Ru, Zr and Ti). In general, the mafic intrusions exhibit significant light REE enrichment and heavy REE depletion and have no Eu anomalies. Whole-rock neodymium (ε_Nd(t) = 1.55) and zircon Hf (ε_Hf(t) = 0.60–3.73) isotopic compositions indicate derivation of the magma from enriched type I mantle. We propose that the diabase porphyrite and diabase formed in a continental margin rift setting, influenced by the Kerguelen mantle plume, and represent magmatism related to the breakup of East Gondwana. However, the gabbro formed in a relatively stable continental intraplate environment, likely derived from deep magmatic processes associated with the Kerguelen mantle plume. Our results provide new constraints on the early activity of the Kerguelen mantle plume and offer insights into the breakup and tectonic evolution of East Gondwana.

1. Introduction

The Tethyan Himalayan Belt is a key geological unit of the Tibetan Plateau and records continental rifting, lithospheric thinning, and mantle plume activity [1,2,3,4,5,6,7,8,9,10,11]. Mantle plume activity links large igneous provinces (LIPs) and the rifting of (super)continents [12]. The Kerguelen mantle plume was one of the largest mantle plumes in Earth history and had a significant role in volcanic activity in the central Indian Ocean and rifting of East Gondwana [13,14,15,16,17,18] (Figure 1a). As such, the Kerguelen mantle plume has been a focus of research. Some studies have suggested that during the Early Cretaceous, the Kerguelen mantle plume may have accelerated the rifting of East Gondwana, leading to the subsequent separation of the Indian, Australian, and Antarctic plates [13,16].

The Tethyan Himalayan Belt contains extensive Late Jurassic to Early Cretaceous mafic igneous rocks, which can be broadly classified into three types: ocean island basalts (OIBs), normal-type mid-ocean ridge basalts (N-MORBs), and transitional rocks. These rocks generally formed in a tectonic setting characterized by lithospheric extension and thinning along the northern margin of Gondwana [19,20,21,22,23,24]. OIB-type mafic rocks are considered to be the products of rifting of the East Gondwana continent and provide key insights into the regional magmatism and tectonic evolution [25,26,27,28,29]. Previous studies have suggested that an extensional tectonic setting existed during the Early Cretaceous in both the eastern and western segments of the Tethyan Himalayan tectonic belt [30]. However, the petrogenesis and tectonic setting of igneous rocks in this region are debated. Some studies have proposed that Early Cretaceous mafic magmatism in the eastern Tethyan Himalayan Belt was closely associated with the Kerguelen mantle plume, suggesting the plume was the primary driver of rifting in East Gondwana [15,25,27,30,31,32,33]. In contrast, other studies have proposed that the Early Cretaceous magmatism resulted from partial melting of enriched asthenospheric mantle, with no direct connection to the Kerguelen mantle plume. In this case, rifting of East Gondwana was driven primarily by passive rifting between the Indian and Australian–Antarctic plates [7,26,27,28,29,30,31,32,33,34].

To better understand the formation of mafic intrusions in this region and the relationship between the Kerguelen mantle plume and rifting of East Gondwana, we undertook a petrological, geochronological, and geochemical study of mafic intrusions in the Cuona area of the eastern Tethyan Himalayan Belt. By determining the ages, sources, and tectonic settings of the intrusions, we provide new constraints on the early activity of the Kerguelen mantle plume and the rifting and tectonic evolution of East Gondwana.

Figure 1

(a) Distribution of igneous rocks related to the Kerguelen mantle plume in the Indian Ocean and adjacent continental margins [2,4,31]. (b) Geological map of the Tethyan Himalaya and adjacent areas [35,36]. (c) Geological map of the Cuomei large igneous province showing the distribution of igneous rocks with ages of >130 Ma [10,16].

[Figure omitted. See PDF]

2. Geological Setting and Sample Descriptions

The Himalayan Orogenic Belt was part of the Rodinia and Gondwana supercontinents and has experienced both convergent and passive continental margin tectonism [37,38]. The Himalayan Orogenic Belt is divided into the Tethyan Himalaya, Great Himalaya, and Lesser Himalaya [35]. The Tethyan Himalayan Belt is the northernmost tectonic unit of the Himalayan Orogenic Belt and is located south of the Yarlung Zangbo Suture Zone and north of the High Himalayan Belt. The Tethyan Himalayan Belt is E–W-trending (Figure 1b). The Tethyan Himalaya consists primarily of Paleozoic–Mesozoic marine sedimentary rocks, within which mafic volcanic rocks are widely developed. The magmatism in this region occurred mainly from the Permian to Cretaceous [39,40]. The Tethyan Himalayan Belt can be divided into southern and northern zones. The southern zone consists of weakly metamorphosed sedimentary rocks from the Indian continental margin, while the northern zone is characterized by deep-water sedimentary rocks, including flysch and ophiolitic fragments.

Igneous rocks are widely distributed in the Tethyan Himalayan Belt and mainly intrude Triassic, Jurassic, and Lower Cretaceous sedimentary rocks. The igneous rocks consist mainly of basaltic volcanic rocks, diabase dikes/sills, and gabbroic intrusions, which are typically exposed in an E–W-trending belt. The study area is located in the Cuona region of the eastern Tethyan Himalaya (Figure 1c). Jurassic and Cretaceous sedimentary rocks interlayered with volcanic rocks are widely distributed in this region, occurring in the Middle to Upper Jurassic Zhela, Lure, and Sangxiu formations. The Zhela Formation consists mainly of siltstone, black shale, and graywacke interlayered with chert, limestone, and volcanic rocks. The Lure Formation is dominated by dark gray muddy (micritic) limestone, which is locally interlayered with basalt, as well as dark gray calcareous silty slate interbedded with dark gray, muddy, fine-grained limestone.

Seven samples were collected from near Juela Township in Cuona County, within the eastern Tethyan Himalayan Belt; the main rock types are diabase porphyrite (21CN-01, 21CN-02, 21CN-04, 21CN-08), gabbro (21CN-06), and diabase (21CN-07-1, 21CN-07-2) (Figure 2). The samples were collected from mafic dikes in the Lure and Zhela formations. Field investigations revealed that these intrusions have clear intrusive contacts and are extensively exposed at the surface (Figure 3a–c).

The diabase porphyrite has gray–green fresh surfaces and weathers to a light brown color (Figure 3d). It is characterized by a diabasic texture and massive structure. The primary minerals are plagioclase (~50 vol.%) and clinopyroxene (~35 vol.%). Plagioclase occurs mainly as subhedral–euhedral lath-like and tabular crystals that exhibit significant alteration, with some grains being opaque due to sericitization and chloritization. These crystals display oscillatory zoning visible under polarized light, suggesting rapid crystallization from a magmatic melt. Some plagioclase grains show fine-grained sericite filling the interstitial spaces between crystals, indicating low-temperature alteration. Clinopyroxene forms anhedral granular to short prismatic crystals that exhibit chloritization. The groundmass consists mainly of microcrysts of plagioclase, pyroxene, and biotite. Minor minerals include quartz and amphibole, with quartz appearing as small, subhedral grains often found in the interstitial spaces. Accessory minerals include magnetite, which forms opaque grains, and apatite, which appears as fine, needle-like crystals (Figure 3g).

The gabbro has dark gray fresh surfaces and weathers to a yellow–brown color (Figure 3e). It is characterized by a gabbroic texture and massive structure. The primary minerals are plagioclase (~45 vol.%) and clinopyroxene (~40 vol.%). Plagioclase is commonly subhedral, altered to clay minerals, and has a patchy texture indicative of replacement by K-feldspar. Clinopyroxene typically occurs as subhedral–anhedral, short–prismatic crystals that are interstitial to plagioclase and exhibit partial chloritization. The accessory minerals are magnetite, titanite, and apatite, which occur as irregular or acicular crystals that are randomly distributed. Micro-cracks are well developed and infilled with crystalline quartz and carbonate veinlets (Figure 3h). Additionally, some clinopyroxene crystals display distinct edge alteration features, which may be related to compositional exchange during the crystal growth process. The plagioclase shows some degree of fractures and altered textures, which may suggest that the rock has undergone an alteration process, such as weathering or hydrothermal alteration. Some titanite crystals exhibit a subtle metallic luster, which may reflect the high-temperature environment in which they formed.

The diabase has grayish blue fresh surfaces and weathers to a yellow–brown color (Figure 3f). It has a diabasic texture and massive structure. The primary mineral is plagioclase (~55 vol.%), which is subhedral–euhedral and exhibits significant alteration, including chloritization, carbonatization, and silicification. Some plagioclase microcrysts are irregularly arranged, indicative of variable degrees of alteration. These crystals dis- play oscillatory zoning, visible under polarized light, suggesting rapid crystallization from a magmatic melt. Secondary minerals include chlorite and calcite, and accessory minerals are magnetite and apatite (Figure 3i). Magnetite occurs as opaque, subhedral to anhedral grains, while apatite appears as fine, needle-like crystals interspersed in the groundmass. The rock also contains minor quartz, which appears as small, subhedral grains, commonly located in the interstitial spaces between plagioclase crystals.

3. Methods

3.1. Zircon U–Pb Dating

Two samples (21CN-02) was selected for zircon U–Pb dating by laser ablation–inductively coupled plasma–mass spectrometry (LA–ICP–MS). Zircon separation was conducted by the Shangyi Rock and Mineral Testing Technology Service Company Limited, Langfang, Hebei Province, China. Zircon mounting and cathodoluminescence (CL) imaging were undertaken by the Yujin Technology Company Limited in Chongqing, China. The analyses were carried out in the Guangxi Key Laboratory of Hidden Metallic Ore Exploration, Guilin University of Technology, Guilin, Guangxi, China. The LA–ICP–MS instrument comprised a New Wave 193 nm LA system coupled to an Agilent 7500 ICP–MS instrument. The laser was operated at a frequency of 5 Hz, and high-purity He was used as the carrier gas. A laser spot size of 32 μm was used for the analyses. Zircon ages were calculated by use of the standard zircon Plesovice (PL) as an external standard, and elemental concentrations were calibrated by use of the synthetic silicate glass standard NIST610 from the United States National Institute of Standards and Technology (NIST). Data processing was performed by use of ICP–MS DataCal 7.2 software [41], with common Pb corrected following [42]. Zircon U–Pb ages were calculated by use of concordia diagrams and weighted-mean age calculations, both using Isoplot 3.2 software [43]. The detailed instrumental operating conditions and data-processing methods were described by Wang et al., (2006) [44] and Yuan et al., (2004) [45].

3.2. Major and Trace Elements

The contents of major and trace elements in seven samples of diabase porphyry (21CN-01, 21CN-02, 21CN-04, 21CN-08), gabbro (21CN-06), and diabase (21CN-07-1, 21CN-07-2) were determined in the Key Laboratory of Guangxi Concealed Metal Ore Exploration of Guilin University of Technology. The rock samples were washed and dried, and weathered surfaces were removed. The samples were then crushed to <5 mm and ground into powder finer than 200 mesh. Major elements were analyzed by use of an alkali fusion glass bead method, with a ZSX Primus II X-ray fluorescence spectrometer. The precision of the major element analyses was better than ±10%. We determined trace elements with an Agilent 7500cx ICP–MS instrument after acid digestion, achieving a precision better than ±5%. Our analytical procedures followed Li et al., (2016) [46].

3.3. Sr–Nd Isotopes

The Sr–Nd isotopic compositions of sample (21CN-02) was determined at the Guangxi Key Laboratory of Hidden Metallic Ore Exploration, Guilin University of Technology. The samples were fully dissolved by acid digestion. Rubidium, Sr, and the rare earth elements (REEs) were separated by use of cation exchange resin, and Sm and Nd were further separated by use of HDEHP resin. The Sr and Nd fractions were analyzed with a Neptune Plus multi-collector–ICP–MS instrument. The ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd ratios were corrected for instrumental mass bias by normalization to ⁸⁶Sr/⁸⁸Sr = 0.1194 and ¹⁴⁶Nd/¹⁴⁴Nd = 0.7219, respectively. The details of the Sr–Nd isotopic analyses were described by Liang et al., (2003) [47].

3.4. Zircon Hf Isotopes

The Zircon Hf isotopic analysis was conducted on one sample (21CN-02) at the Guangxi Key Laboratory of Hidden Metallic Ore Exploration, Guilin University of Technology. The Hf isotopic analyses were undertaken near the zircon U–Pb dating sites. The Hf isotopic analyses were undertaken with a Neptune Plus multi-collector–ICP–MS instrument coupled to a GeoLas HD excimer ArF LA system. The laser spot diameter was set to 32 μm, and the laser frequency was 8 Hz. Crustal-model Hf isotopic ages were calculated by assuming the parental magma was derived from average continental crust. The ¹⁷⁶Lu/¹⁷⁷Hf ratio of 0.015 represents that of a depleted mantle source [48].

4. Results

4.1. Zircon U–Pb Ages

The zircon U–Pb dating results are presented in Table 1. The zircon grains from both samples are primarily euhedral–subhedral, with lengths of 50–100 μm, aspect ratios of 1:1 to 2:1, and prismatic shapes (Figure 4a). The U contents in the zircons from the diabase porphyrite sample (21CN-02) are 494–1050 ppm, and the Th contents are 462–1432 ppm. The Th/U ratios are 0.94–1.59, with an average of 1.17. These values are consistent with those of typical magmatic zircon.

After excluding discordant ages and those with large errors, the analyses were used to construct concordia diagrams. The weighted-mean ²⁰⁶Pb/²³⁸U age for sample 21CN-02 is 135.0 ± 1.6 Ma (MSWD = 0.68; n = 13) (Figure 4b).

4.2. Major and Trace Elements

The major element data are presented in Table 2. The SiO₂ contents of the mafic intrusive rock samples are 49.24–51.95 wt.%. The Al₂O₃ contents are 10.51–15.58 wt.%, and the Na₂O contents are 2.18–5.38 wt.%. The average K₂O content is 0.61 wt.%, and the total alkali (K₂O + Na₂O) contents are 2.20–5.40 wt.%. Mg^# values vary from 44.2 to 66.2. The loss-on-ignition (LOI) values are 1.26–5.42 wt.%, indicating these mafic intrusive rocks have undergone variable degrees of alteration. In a Nb/Y–Zr/TiO₂ diagram, the samples plot in the subalkaline and alkaline basalt fields (Figure 5a), indicating transitional geochemical signatures between subalkaline and alkaline affinities. In an A/CNK–A/NK diagram, the samples plot in the metaluminous field (Figure 5b). Due to their relatively high alumina saturation index (A/CNK; 0.52–0.92) and Laitman index (σ; 0.76–3.77) values, the samples plot in the calc-alkaline to metaluminous fields on an AR–SiO₂ classification diagram (Figure 5c), showing a trend in transition from sub-alkaline to alkaline components. In a SiO₂–K₂O diagram, the samples plot in the calc-alkaline to low-K (tholeiitic) fields (Figure 5d).

The trace element data are listed in Table 2. The total REE contents of the mafic intrusive rock samples range from 98.2 to 308.2 ppm. The REE patterns are light-REE-enriched (Figure 6a). The light/heavy REE ratios vary from 6.55 to 13.84, and (La/Yb)_N = 4.25–13.09. δEu values range from 0.93 to 1.09, with an average of 1.02, indicating there are no significant Eu anomalies. In a primitive-mantle-normalized trace element diagram (Figure 6b), the diabase porphyrite and gabbro have similar trace element characteristics, with enrichments in high-field-strength elements (HFSEs, e.g., Nb and Ti) and large-ion lithophile elements (LILEs, e.g., Sr and Pb), similar to OIBs. The diabase is enriched in large-ion lithophile elements (LILEs, e.g., La and Ce) and depleted in high-field-strength elements (HFSEs, e.g., Ru, Zr, and Ti).

4.3. Sr–Nd Isotopes

The Sr–Nd isotopic data are listed in Table 3. The diabase has (⁸⁷Sr/⁸⁶Sr)_i = 0.705583, (¹⁴³Nd/¹⁴⁴Nd)_i = 0.512541, and ε_Nd(t) = 1.51. The diabase porphyrite has relatively high (⁸⁷Sr/⁸⁶Sr)_i and low ε_Nd(t) values, with T_DM ages of 1168 Ma (Figure 7).

4.4. Zircon Hf Isotopes

The in situ zircon Hf isotopic data are listed in Table 4. The diabase porphyrite has an average zircon ¹⁷⁶Lu/¹⁷⁷Hf value of 0.000803 (1σ; n = 14) and an average zircon ¹⁷⁶Hf/¹⁷⁷Hf value of 0.2828746 (1σ; n = 14), which represents the initial ¹⁷⁶Hf/¹⁷⁷Hf ratio, given the low Lu/Hf ratio. The ε_Hf(t) values range from 0.60 to 3.73 (Figure 8), with T_DM Hf model ages of 641–776 Ma.

5. Discussion

5.1. Ages

Zircon U–Pb dating was undertaken on two samples of mafic intrusive rock from the Cuona area in the eastern Tethyan Himalayan Belt of southern Tibet. The diabase porphyrite sample (21CN-02) has a crystallization age of 135.0 ± 1.6 Ma. The age of diabase porphyrite indicates that there is Early Cretaceous mafic magmatism in this area.

We compiled crystallization ages for Late Jurassic to Early Cretaceous mafic rocks in the Tethyan Himalayan Belt and constructed an age histogram (Figure 9). On the basis of previous studies, Late Jurassic to Early Cretaceous mafic magmatism in the Tethyan Himalayan Belt occurred mainly during 147–129 Ma. The crystallization age of the diabase porphyrite of this study (135.0 ± 1.6 Ma) is generally consistent with diabase crystallization ages (137–131 Ma) previously reported from the Cuona area [2,4]. It is also consistent with the Late Jurassic to Early Cretaceous magmatism in the Tethyan Himalayan Belt (Figure 9). Therefore, the studied mafic intrusive rocks may be the products of two distinct magmatic events. The diabase porphyrite corresponds to Late Jurassic to Early Cretaceous magmatism in the Tethyan Himalayan Belt (147–129 Ma), whereas the gabbro may represent the initial activity of the Kerguelen mantle plume during the Late Jurassic to Early Cretaceous.

5.2. Fractional Crystallization and Crustal Contamination

Alteration can complicate the use of geochemical data to constrain the petrogenesis of igneous rocks. Some trace elements in mafic rocks, such as REEs, HFSEs, and transition metals (e.g., V, Cr, and Ni), are barely affected by low-temperature alteration or weathering [55,56,57]. The mafic rock samples from the Cuona area have LOI values of 1.26–5.42 wt.%, indicating they have undergone variable degrees of alteration. Therefore, elements that are relatively immobile during alteration are used in the following discussion.

Primary magmas from a mantle source typically have high Mg^# values (73–81) and elevated Cr (>1000 ppm) and Ni (>400 ppm) concentrations [58]. The studied samples have Mg^# = 43.80–66.18, Cr = 14.8–451 ppm, and Ni = 6.5–259 ppm, which are lower than the values for primary mantle-derived magmas. This indicates that the samples have undergone fractional crystallization. Both Cr and Ni exhibit positive correlations with MgO (Figure 10d–e), indicative of the early fractional crystallization of olivine or clinopyroxene. δEu values of the mafic rock samples are 0.93–1.09, with an average of 1.02, indicative of no significant plagioclase fractionation. Al₂O₃ exhibits a negative correlation with MgO (Figure 10a), which is also indicative of little plagioclase fractionation. These results suggest the magma chemistry was affected mainly by clinopyroxene fractionation. Moreover, CaO exhibits a positive correlation with MgO and a negative correlation with SiO₂ (Figure 10b,f), which is consistent with clinopyroxene fractionation being dominant. TiO₂ exhibits a negative correlation with MgO (Figure 10c) because TiO₂ was an incompatible element that was progressively enriched in the magma during early fractionation. As the magma evolved further, late fractional crystallization of magnetite occurred, which caused a decrease in TiO₂ concentrations.

The diabase porphyrite samples have La/Ta = 18.2–24.4 and La/Sm = 2.33–3.12, while the gabbro samples have La/Ta = 17.0 and La/Sm = 2.50, indicating these rocks were not significantly contaminated by the crust or lithospheric mantle [59]. Moreover, the diabase porphyrite samples have Zr/Hf = 38.08–38.62 and Th/Ta = 1.5–2.40, while the gabbro samples have Zr/Hf = 38.13 and Th/Ta = 1.45, which are both different from typical crustal values [53,60]. The diabase samples have La/Ta = 26.1–27.6 and La/Sm = 2.80–3.06, indicative of no significant crustal contamination but some degree of lithospheric mantle contamination. In addition, the diabase samples have Zr/Hf = 38.91–40.51 and Th/Ta = 1.94–2.12, which are also markedly different from typical crustal values. In a La/Sm–La/Nb diagram (Figure 11a), all samples plot far from the crustal contamination trend line, indicative of insignificant crustal contamination. In a (Th/Ta)_PM–(La/Nb)_PM diagram (Figure 11b), the diabase porphyrite and gabbro samples plot near the primitive mantle field and away from the crustal field. Some diabase porphyrite and gabbro samples plot in the OIB-type igneous rock field associated with the Kerguelen mantle plume. The diabase samples plot farther away from the middle and upper crustal fields and closer to the lower crustal field. In a Nb/Yb–Th/Yb diagram (Figure 11c), the diabase porphyrite, gabbro, and diabase samples all exhibit characteristics similar to those of rocks from the Kerguelen Islands. They exhibit an overall trend towards the intraplate enrichment field, with no involvement of subduction-related or crustal materials. In a (La/Nb)_PM–(Th/Nb)_PM diagram (Figure 11d), the diabase porphyrite and gabbro samples plot in the Cuomei LIP field, close to the primitive mantle values. In contrast, the diabase samples plot closer to the lower crustal field, indicative of possible interactions with the lower crust.

In summary, the diabase porphyrite and gabbro show no features of crustal or lithospheric mantle contamination, and the diabase is contaminated by a certain degree of lithospheric mantle.

5.3. Magma Sources and Tectonic Setting

In a (⁸⁷Sr/⁸⁶Sr)_i–(¹⁴³Nd/¹⁴⁴Nd)_i diagram (Figure 7a), the samples plot near the field for the Kerguelen Islands and Kerguelen mantle plume. In a (⁸⁷Sr/⁸⁶Sr)_i–ε_Nd(t) diagram (Figure 7b), the samples plot within the fields for the Cuomei LIP and Kerguelen mantle plume and have similar isotopic compositions as the gabbro of the Lakang Formation in the southeastern Tethyan Himalayan Belt (Wang et al., 2022) [18]. In a (¹⁴³Nd/¹⁴⁴Nd)_i–(¹⁷⁶Hf/¹⁷⁷Hf)_i diagram (Figure 8b), the diabase porphyrite plots within the field for the Kerguelen Islands, and its isotopic composition is similar to that of the Cuomei LIP (Xia et al., 2012) [22], with (⁸⁷Sr/⁸⁶Sr)_i = 0.7047 and ε_Nd(t) = 1.5. These characteristics suggest the diabase porphyrite has similar isotopic compositions as other Early Cretaceous OIB-type mafic rocks in the Tethyan Himalayan Belt, and were associated with the Kerguelen mantle plume.

The average Zr/Nb and La/Nb ratios of the studied mafic intrusive rocks are 9.42 and 1.28, respectively, which are close to the values for EM I (Zr/Nb = 4.2–11.5 and La/Nb = 0.86–1.19). The geochemical characteristics of the samples are relatively uniform. In La–La/Nb and Nb–Nb/Th diagrams (Figure 12a,b), all samples plot in the OIB field. In a Zr/Nb–Y/Nb diagram (Figure 12c), all samples plot in the enriched mantle region. In a Ta/Yb–Th/Yb diagram (Figure 12d), the samples plot in the enriched mantle field and close to that for active continental margins. The OIB-type mafic intrusions in this region have EM I geochemical features, similar to other such rocks in several locations within the Tethyan Himalayan Belt in southern Tibet [19,24,29].

In Zr–Zr/Y and Zr–TiO₂ diagrams, all the samples plot in the intraplate basalt field (Figure 13a,b). In a Ta/Hf–Th/Hf tectonic discrimination diagram, the diabase porphyrite samples plot in the initial rift basalt field, the gabbro samples plot in the oceanic intraplate basalt field (including OIBs, seamounts, transitional-type MORBs, and enriched-type MORBs), and the diabase samples plot in the boundary region between the initial and intracontinental rift fields (Figure 13c). In a 2Nb–Zr/4–Y tectonic discrimination diagram, the diabase porphyrite and diabase samples plot in the intraplate alkaline basalt + intraplate tholeiitic basalt field, while the gabbro samples plot in the intraplate alkaline basalt field (Figure 13d). Compared with previous studies, the diabase porphyrite and diabase have similar geochemical characteristics to volcanic rocks in the Sangxiu, Weimei, and Zhela formations in the eastern Tethyan Himalayan Belt [1,8]. This similarity suggests that the diabase porphyrite and diabase may have formed in a continental rift environment, while the gabbro, with its distinct composition, could potentially represent an continental intraplate environment. However, the age differences between these rock types have not been clearly established, and further investigations are needed to support these interpretations.

5.4. Geological Significance

Although most Cretaceous magmatism in East Gondwana was unrelated to the Kerguelen mantle plume, the Early Cretaceous magmatism in the Tethyan Himalayan Belt was closely associated with the plume [26,34]. Previous studies have shown that the OIB-type mafic magmatism in this tectonic belt began at ca. 147 Ma and ended at 125 Ma [2,4,5,24,29,69]. According to an Early Cretaceous paleogeographic reconstruction of Gondwana, the Tethyan Himalayan Belt was located along a passive continental margin at the time. There is a lack of arc-related igneous rocks that coexist with the Early Cretaceous OIB-type mafic rocks, which rules out the possibility of their formation in a subduction setting [70,71]. During the Late Jurassic to Early Cretaceous, the Cuona area was located along the passive continental margin of eastern Gondwana, adjacent to the Neo-Tethys Ocean. The rapid expansion of the Neo-Tethys Ocean caused intense extensional tectonism south of the mid-ocean ridge, leading to the formation of rift structures [8,27,72]. OIB-type igneous rocks are generally considered to be the products of a mantle plume or hotspot activity [73]. The Early Cretaceous OIB-type mafic rocks in the Tethyan Himalayan Belt were closely related to the Kerguelen mantle plume, which had a key role in the rifting of eastern Gondwana [16,25,28,30].

The diabase porphyrite formed at 135.0 ± 1.6 Ma, which is slightly earlier than the peak of Kerguelen mantle plume activity (ca. 132 Ma; Zhu et al., 2008; Zhu et al., 2009; Ma et al., 2016) [2,16,74], and was coeval with the Bunbury basalts (137–136 Ma) and the oceanic magnetic anomaly (136 Ma) that record the rifting of eastern Gondwana [34,75]. Furthermore, the diabase porphyrite has distinct OIB-type geochemical characteristics, with enrichments in light REEs and HFSEs, indicative of a mantle plume source. Its age coincides with the peak of mafic magmatism during the Late Jurassic–Early Cretaceous (147–129 Ma) in the Tethyan Himalayan Belt (Figure 9). Therefore, the diabase porphyrite was formed in a continental rift setting associated with the Kerguelen mantle plume, and represents a product of the rifting of eastern Gondwana.

5.5. Kerguelen Mantle Plume and Evolution of Eastern Gondwana

The initial stages of continental rifting are typically associated with LIPs [76,77]. There are differing views on the relationship between the rifting of eastern Gondwana and the Kerguelen mantle plume. Some studies have suggested there is no direct relationship between the two [26]. However, there is some evidence that the Kerguelen mantle plume had a key role in the rifting of East Gondwana. For example, the bimodal magmatism during 147–130 Ma, the formation of A-type granites [6,25], and the development of the Indian Ocean Ridge at 136 Ma were all associated with mantle plume activity. Some studies have suggested that the small magmatic provinces on East Gondwana were associated with the activity of the Kerguelen mantle plume [16,18,26,34]. The magmatism associated with the plume lasted for an extended period, with intermittent eruptions for potentially >100 Myr on both sides of the Southeast Indian Ridge [78]. These lines of evidence suggest that the prolonged magmatic activity of the Kerguelen mantle plume played a key role in the rifting of East Gondwana and potentially caused the rifting.

As East Gondwana began to rift, the Indian Plate gradually moved northward, and the Kerguelen mantle plume activity became increasingly focused in the Tethyan Himalayan region. At ca. 147 Ma, East Gondwana underwent initial rifting and the Cuona area in the eastern Tethyan Himalayan Belt was located in the central part of the Kerguelen mantle plume head (Figure 14a; Zhang et al., 2023) [10]. As the rifting of East Gondwana proceeded, the Cuona area began to drift towards the edge of the Kerguelen mantle plume, and the mantle plume magmatism peaked, which formed a large volume of OIB-type basalts at ca. 135 Ma (Figure 14b). With the further expansion of the mantle plume head, younger OIB-type volcanic rocks (125 Ma) formed (Dong et al., 2019) [70]. The bimodal magmatism (118–115 Ma) and N-MORB-type basalts (ca. 120 Ma) in the eastern Tethyan Himalaya are the products of extension rather than Kerguelen mantle plume activity [2,3,5,16]. Furthermore, the Australia–Antarctic plate had a relatively stable paleo-latitude during 140–120 Ma (Torsvik et al., 2012) [79], indicating the Indian Plate had completely separated from East Gondwana by ca. 125 Ma (Figure 14c).

Our results indicate that early activity of the Kerguelen mantle plume not only facilitated the rifting of East Gondwana, but also affected the tectonic and magmatic evolution of the Tethyan Himalaya. Our study has provided new constraints on the early activity of the Kerguelen mantle plume, which had a key role in the geological evolution of the Tethyan Himalaya, and the origins of basaltic magmatism in continental margin rift and intraplate settings during the rifting of East Gondwana.

6. Conclusions

We investigated the petrology, geochronology, and geochemistry of mafic intrusive rocks in the Cuona region of the eastern Tethyan Himalaya. Our main conclusions are as follows:

• (1). The diabase porphyrite formed in the Early Cretaceous (135.0 ± 1.6 Ma). The diabase porphyrite and gabbro have a similar petrogenesis to the numerous OIB-type mafic igneous rocks in the Cuomei LIP, displaying typical OIB-type geochemical signatures without contamination by lithospheric mantle or crust. These rocks were likely derived from EMI-like enriched mantle. The diabase porphyrite formed in a continental rift setting related to the Kerguelen mantle plume, whereas the gabbro should be the product of the deep magma of the Kerguelen mantle plume in a stable continental intraplate environment.

• (2). The mafic intrusive rocks in the Cuona region are subalkaline to alkaline. There is significant fractionation between the light (enriched) and heavy (depleted) REEs and no significant Eu anomalies. The trace element characteristics of the diabase porphyrite and gabbro are similar, with both being enriched in HFSEs and LILEs. The diabase is enriched in large-ion lithophile elements (LILEs, e.g., La and Ce) and depleted in high-field-strength elements (HFSEs, e.g., Ru, Zr and Ti).

• (3). The age of the diabase porphyrite (135 Ma) is slightly older than the peak of Kerguelen mantle plume activity (132 Ma), indicating it was the product of East Gondwana continental rifting associated with the effects of the Kerguelen mantle plume.

Author Contributions

Conceptualization, F.Y. and C.M.; methodology, Z.Q.; software, Z.L.; validation, F.Y.; formal analysis, J.M.; investigation, Z.L.; resources, L.W.; data curation, Z.Q.; writing—original draft preparation, C.S.; writing—review and editing, Z.K.; visualization, L.Z.; supervision, F.Y.; project administration, C.S.; funding acquisition, Z.K. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 2. Simplified geological map of the Cuona area in the eastern segment of the Tethyan Himalayan Belt (modified from the 1:50,000 Juela Geological Map).

Figure 3. (a–c) Field photographs, (d–f) hand specimens, and (g–i) microscopic characteristics of mafic intrusions in the Cuona area. Pl = plagioclase; Cpx = clinopyroxene; Mt = magnetite; Bt = biotite; Kfs = K-feldspar; Cal = calcite; Chl = chlorite.

Figure 4. (a) Cathodoluminescence (CL) images from mafic intrusions in the Cuona area, (b) U–Pb concordia diagrams of zircon grains from mafic intrusions in the Cuona area.

Figure 5. Geochemical classification diagrams for mafic intrusive rocks in the Cuona area and related rocks [49,50,51,52]. (a) Nb/Y–Zr/TiO2 diagram; (b) A/CNK–A/NK diagram; (c) AR–SiO2 classification diagram; (d) SiO2–K2O diagram. Data for regional OIB-type mafic rocks are from Zhu et al., (2008), Chen et al., (2021), and Yang et al., (2022) [2,4,29].

Figure 6. (a) Chondrite-normalized rare earth element (REE) patterns and (b) primitive-mantle-normalized trace element diagram for mafic intrusive rocks in the Cuona area and related rocks (chondrite, primitive mantle, OIB, N-MORB, and E-MORB data are from Sun and McDonough, 1989 [53]; OIB-type mafic rock data associated with the Kerguelen mantle plume are from Zhu et al., 2008; Chen et al., 2021; Yang et al., 2022) [2,4,29].

Figure 7. (a) (87Sr/86Sr)i versus (143Nd/144Nd)i and (b) (87Sr/86Sr)i versus εNd(t) diagrams for mafic intrusive rocks in the Cuona area and related rocks (Kerguelen mantle plume, Kerguelen Islands, Kerguelen Plateau, Ninetyeast Ridge, and Rajmahal–Bengal–Sylhet Traps data are from Olierook et al., 2019; Peng et al., 2022; gabbro data are from Wang et al., 2022) [8,18,26].

Figure 8. (a) Zircon εHf(t) and (b) (143Nd/144Nd)i–(176Hf/177Hf)i diagrams for mafic intrusive rocks in the Cuona area and related rocks (Kerguelen Islands and Tethyan Himalayan Igneous Province data are from Chen et al., 2018 and Chen et al., 2021) [3,4].

Figure 9. Histogram of ages for Late Jurassic to Early Cretaceous mafic rocks in the Tethyan Himalaya (age data are from Zhu et al., 2009; Hou Chenyang, 2017; Shi et al., 2017; Chen et al., 2021; Deng Jin, 2021; Liang et al., 2023) [4,16,36,54].

Figure 10. Major elements variation diagram for mafic rocks from the Cuona area. (a) MgO-Al2O3 diagram; (b) MgO-CaO diagram; (c) MgO-TiO2 diagram; (d) MgO-Cr diagram; (e) MgO-Ni diagram; (f) SiO2-CaO diagram.

Figure 11. Plots used to identify crustal contamination of mafic intrusions in the Cuona area: (a) La/Sm versus La/Nb; (b) Th/Ta versus La/Nb; (c) Nb/Yb versus Th/Yb; (d) (La/Nb)PM versus (Th/Nb)PM [1,2,34]. Data sources: upper, middle, and lower crust [61]; primitive mantle [62]; lithospheric mantle [63]; OIB-type igneous rocks related to the Kerguelen mantle plume [1]. OIB = ocean island basalt; E-MORB = enriched mid-ocean ridge basalt; N-MORB = normal mid-ocean ridge basalt; UC = upper crust; MC = middle crust; LC = lower crust; PM = primitive mantle; CLM = continental lithospheric mantle.

Figure 12. (a) La versus La/Nb ((a) is from Shuguang Li, 1993) [64], (b) Nb versus Nb/Th ((b) is from Shuguang Li, 1993) [64], (c) Zr/Nb versus Y/Nb ((c) is from Wilson, 1997) [58], and (d) Ta/Yb versus Th/Yb diagrams ((d) is from Meschede, 1986) [65] for mafic intrusive rocks in the Cuona area. OIB = ocean island basalt; MORB = mid-ocean ridge basalt; IAB = island arc basalt; N-MORB = normal mid-ocean ridge basalt; CLM = continental lithospheric mantle.

Figure 13. (a) Zr/Y versus Zr ((a) is from Pearce and Norry, 1979) [66], (b) Zr versus Ti ((b) is from Pearce, 1982) [67], (c) Th/Hf versus Ta/Hf ((c) is from Wang et al., 2001) [68], and (d) 2Nb–Zr/4–Y tectonic discrimination diagrams ((d) are from Meschede, 1986) [65] for mafic intrusive rocks in the Cuona area. WPB = within-plate basalt; IAB = island arc basalt; MORB = mid-ocean ridge basalt; A1 = within-plate alkaline basalts; A2 = within-plate alkaline basalts + within-plate tholeiites; B = E-MORBs; C = within-plate tholeiites + volcanic arc basalts; D = volcanic arc basalts + E-MORBs.

Figure 14. Model of the evolution of the Kerguelen mantle plume and East Gondwana (modified from Gibbons et al., 2013; Olierook et al., 2019; Wang et al., 2022; Liang et al., 2023) [9,18,26,77]. (a) At ca. 147 Ma, East Gondwana underwent initial rifting and the Cuona area in the eastern Tethyan Himalayan Belt was located in the central part of the Kerguelen mantle plume head. (b) At ca. 135 Ma, as the rifting of East Gondwana proceeded, the Cuona area began to drift towards the edge of the Kerguelen mantle plume, and the mantle plume magmatism peaked, which formed a large volume of OIB-type basalts. (c) By ca. 125 Ma, the Indian Plate had completely separated from East Gondwana, as indicated by the paleo-latitude stability of the Australia–Antarctic Plate.

Table 1

Zircon U-Th-Pb isotopic dating results of mafic intrusive rocks in the Cuona region.

Sample Point	Content/10−6			Th/U	Isotopic Ratio				Isotopic Age/Ma
	Pb	Th	U		207Pb/235U	1σ	206Pb/238U	1σ	207Pb/235U	1σ	206Pb/238U	1σ	Con*
21CN-02-06	21.77792	755.0443	739.7702	1.02	0.14585	0.00607	0.02159	0.00023	138.2	5.4	137.7	1.4	99%
21CN-02-11	32.31679	1302.509	1049.637	1.24	0.14207	0.00551	0.02125	0.00025	134.9	4.9	135.5	1.6	99%
21CN-02-14	14.03134	462.2057	493.6375	0.94	0.14318	0.00859	0.02134	0.00030	135.9	7.6	136.1	1.9	99%
21CN-02-16	18.90271	701.7517	629.2832	1.12	0.14162	0.00695	0.02179	0.00029	134.5	6.2	138.9	1.8	96%
21CN-02-17	16.85502	588.0218	581.0979	1.01	0.13800	0.00694	0.02094	0.00030	131.3	6.2	133.6	1.9	98%
21CN-02-30	15.79062	511.8948	516.887	0.99	0.15149	0.01039	0.02152	0.00037	143.2	9.2	137.2	2.3	95%
21CN-02-31	30.23075	1272.806	971.9639	1.31	0.14006	0.00910	0.02047	0.00025	133.1	8.1	130.6	1.6	98%
21CN-02-33	30.33958	1412.224	916.9356	1.54	0.13682	0.00604	0.02101	0.00024	130.2	5.4	134.1	1.5	97%
21CN-02-34	24.13773	930.2534	796.9868	1.17	0.13195	0.00561	0.02113	0.00025	125.9	5.0	134.8	1.6	93%
21CN-02-37	23.66793	838.1252	767.2604	1.09	0.13738	0.00566	0.02177	0.00025	130.7	5.1	138.8	1.6	93%
21CN-02-38	29.24066	1148.345	976.7321	1.18	0.13456	0.00550	0.02063	0.00026	128.2	4.9	131.6	1.6	97%
21CN-02-39	29.71245	1432.31	901.637	1.59	0.13925	0.00604	0.02096	0.00024	132.4	5.4	133.7	1.5	98%
21CN-02-40	16.04238	557.4163	553.4698	1.01	0.13057	0.00650	0.02103	0.00029	124.6	5.8	134.1	1.8	92%

Con*: Concordance.

Table 2

Major (wt.%), trace (ppm), and rare earth element (ppm) analysis results of mafic intrusive rocks from the Cuona region.

Sample	21CN-01	21CN-02	21CN-04	21CN-06	21CN-07-1	21CN-07-2	21CN-08
Rock type	Diabase porphyrite	Diabase porphyrite	Diabase porphyrite	Gabbro	Diabase	Diabase	Diabase porphyrite
Al2O3	13.06	13.31	13.70	12.82	13.76	14.68	9.99
CaO	7.28	9.67	10.77	6.67	5.79	4.20	8.68
Fe2O3T	11.52	10.04	11.00	12.80	13.10	12.05	13.17
K2O	0.80	0.98	0.87	1.40	0.03	0.02	0.02
MgO	6.12	6.62	6.49	4.28	4.60	4.09	11.06
MnO	0.16	0.15	0.16	0.20	0.19	0.18	0.17
Na2O	3.98	3.20	2.23	3.15	4.54	5.07	2.07
P2O5	0.40	0.28	0.30	0.45	1.56	1.43	0.29
SiO2	49.05	48.95	49.55	50.89	46.38	48.94	46.94
TiO2	3.99	2.74	2.85	4.58	4.25	3.55	2.68
LOI	2.74	3.17	1.26	2.62	5.42	4.89	4.90
Total	99.11	99.12	99.19	99.87	99.62	99.10	99.93
Li	18.69	9.94	14.80	13.74	25.97	21.01	31.60
Sc	35.47	33.39	38.37	14.39	18.57	13.82	19.99
Ti	24,816.87	17,116.76	17,443.96	28,228.77	27,042.95	22,586.42	16,904.02
V	272.22	231.41	247.41	251.02	243.71	171.81	242.08
Cr	298.87	440.10	179.86	14.81	98.47	14.96	451.08
Mn	1104.18	1066.03	1110.47	1429.79	1355.65	1174.42	1171.26
Co	40.97	39.17	42.95	34.39	29.58	23.20	63.07
Ni	124.74	154.79	67.24	17.70	42.62	6.51	259.45
Cu	36.88	38.13	26.41	13.99	19.85	10.83	32.69
Zn	128.99	109.63	113.31	127.48	135.54	150.72	124.56
Ga	25.89	14.11	14.08	16.11	16.94	14.28	9.41
Rb	17.93	29.81	29.64	18.02	1.67	1.21	0.38
Sr	548.45	437.22	578.08	672.44	1046.50	587.62	79.02
Y	30.22	23.25	25.58	31.12	38.76	40.01	25.63
Zr	240.84	165.96	177.32	275.20	215.56	239.69	162.32
Nb	20.79	15.30	15.40	23.31	29.46	33.46	12.94
Mo	1.78	2.24	0.65	1.53	1.75	1.58	0.72
Cs	0.18	0.32	0.60	0.31	1.12	0.48	0.72
Ba	538.48	195.03	166.93	236.48	214.15	100.33	22.52
La	24.41	16.89	18.18	21.48	44.31	48.41	13.21
Ce	52.31	36.36	38.76	52.73	109.73	116.73	32.96
Pr	6.96	4.87	5.20	6.55	13.67	14.22	4.58
Nd	31.51	22.21	23.54	33.38	62.53	72.87	21.13
Sm	7.84	5.57	6.09	8.59	15.84	15.83	5.67
Eu	2.51	1.87	1.99	2.57	5.35	5.15	1.93
Gd	7.62	5.59	6.08	8.33	14.31	14.20	5.74
Tb	1.19	0.89	0.95	1.33	1.94	1.94	0.94
Dy	6.64	5.00	5.46	7.45	9.71	9.76	5.47
Ho	1.25	0.96	1.02	1.40	1.66	1.68	1.03
Er	3.15	2.40	2.58	3.49	3.81	3.91	2.66
Tm	0.41	0.32	0.35	0.47	0.45	0.47	0.37
Yb	2.52	1.93	2.08	2.83	2.50	2.65	2.23
Lu	0.35	0.27	0.29	0.39	0.33	0.35	0.31
Hf	6.32	4.34	4.59	7.22	5.54	5.92	4.24
Ta	1.00	0.93	0.79	1.27	1.70	1.76	0.61
Tl	0.07	0.05	0.07	0.11	0.01	0.01	0.14
Pb	2.41	1.53	2.25	3.01	3.19	2.75	2.15
Th	1.88	1.40	1.89	1.83	3.29	3.72	1.34
U	0.48	0.33	0.42	0.60	0.74	0.84	0.33
A/CNK	0.63	0.56	0.57	0.68	0.76	0.92	0.52
A/NK	1.76	2.10	2.97	1.91	1.83	1.76	2.92
AR	1.61	1.44	1.29	1.25	1.61	1.61	1.74
Mg#	55.32	60.58	57.89	43.80	45.00	44.17	66.18
FeOT	10.76	9.42	10.11	11.84	12.51	11.51	12.47
∑REE	148.67	105.14	112.56	150.98	286.14	308.16	98.23
LREE	133.16	93.37	99.84	133.62	265.73	287.40	85.22
HREE	15.52	11.77	12.73	17.36	20.40	20.76	13.01
LREE/HREE	8.58	7.93	7.85	7.70	13.02	13.84	6.55
Eu*	0.99	1.03	1.00	0.93	1.09	1.05	1.04
(La/Yb)N	6.96	6.26	6.28	5.45	12.70	13.09	4.25
(Sm/Yb)N	3.46	3.20	3.25	3.37	7.03	6.62	2.83
(La/Sm)N	2.01	1.96	1.93	1.62	1.81	1.98	1.51
(Tb/Yb)N	2.16	2.10	2.10	2.14	3.54	3.34	1.92

Notes: AR = (Al₂O₃ + CaO + Na₂O + K₂O)/(Al₂O₃ + CaO-Na₂O-K₂O).

Table 3

Sr-Nd isotope analysis results of mafic intrusive rocks in the Cuona region.

Sample	Sm(ppm)	Nd(ppm)	Rb(ppm)	Sr(ppm)	147Sm/144Nd	143Nd/144Nd	(143Nd/144Nd)i	εNd(t)	87Rb/86Sr	87Sr/86Sr	(87Sr/86Sr)i	TDM(Ga)
21CN-02	5.57	22.21	29.81	437.22	0.15171	0.51268	0.51254	1.55	0.19725	0.70597	0.70558	1.1668

Table 4

Lu-Hf isotope analysis results of zircons from mafic intrusive rocks in the Cuona region.

Sample Point	Age (Ma)	176Yb/177Hf	176Lu/177Hf	176Hf/177Hf	±2σ	εHf(0)	εHf(t)	(176Hf/177Hf)i	TDM(Hf)	TDM (Ma)	fLu/Hf
21CN-02-02	143	0.02742	0.00063	0.28272	0.000016	−1.83040	1.30543	0.28272	745	1113	−0.981
21CN-02-03	142	0.03072	0.00073	0.28273	0.000017	−1.65096	1.46297	0.28273	739	1102	−0.978
21CN-02-04	149	0.03014	0.00073	0.28271	0.000016	−2.25922	1.00827	0.28271	763	1137	−0.978
21CN-02-11	136	0.02601	0.00061	0.28278	0.000018	0.20608	3.18877	0.28278	664	988	−0.982
21CN-02-12	140	0.02556	0.00059	0.28277	0.000016	−0.16955	2.90089	0.28277	678	1009	−0.982
21CN-02-14	136	0.02367	0.00055	0.28274	0.000018	−1.20716	1.77511	0.28274	719	1078	−0.983
21CN-02-17	134	0.02746	0.00065	0.28280	0.000016	0.79427	3.73320	0.28279	641	951	−0.980
21CN-02-19	148	0.02574	0.00062	0.28274	0.000018	−1.10972	2.13618	0.28274	716	1064	−0.981
21CN-02-20	141	0.03373	0.00081	0.28271	0.000017	−2.39056	0.70116	0.28270	770	1150	−0.976
21CN-02-31	131	0.05426	0.00139	0.28276	0.000035	−0.42895	2.44373	0.28276	701	1031	−0.958
21CN-02-36	139	0.04103	0.00097	0.28276	0.000017	−0.61020	2.43815	0.28275	701	1038	−0.971
21CN-02-38	132	0.03465	0.00084	0.28276	0.000016	−0.41738	2.47727	0.28276	692	1030	−0.975
21CN-02-39	134	0.04543	0.00110	0.28277	0.000017	−0.29332	2.64530	0.28276	691	1021	−0.967
21CN-02-40	140	0.03928	0.00102	0.28270	0.000020	−2.47044	0.59930	0.28270	776	1156	−0.969

References

1. Zhu, D.; Pan, G.; Mo, X.; Liao, Z.; Jiang, X.; Wang, L.; Zhao, Z. Petrogenesis of volcanic rocks in the Sangxiu Formation, central segment of Tethyan Himalaya: A probable example of plume–lithosphere interaction. J. Asian Earth Sci.; 2007; 29, pp. 320-335. [DOI: https://dx.doi.org/10.1016/j.jseaes.2005.12.004]

2. Zhu, D.; Mo, X.; Pan, G.; Zhao, Z.; Dong, G.; Shi, Y.; Liao, Z.; Wang, L.; Zhou, C. 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?. Lithos; 2008; 100, pp. 147-173. [DOI: https://dx.doi.org/10.1016/j.lithos.2007.06.024]

3. Chen, S.; Fan, W.; Shi, R.; Liu, X.; Zhou, X. 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.; 2018; 491, pp. 21-33. [DOI: https://dx.doi.org/10.1016/j.epsl.2018.03.034]

4. Chen, S.; Fan, W.; Shi, R.; Xu, J.; Liu, Y. The Tethyan Himalaya Igneous Province: Early Melting Products of the Kerguelen Mantle Plume. J. Petrol.; 2021; 62, pp. 301-320. [DOI: https://dx.doi.org/10.1093/petrology/egab069]

5. Shi, Y.; Hou, C.; Anderson, J.L.; Yang, T.; Bian, W.; Jin, J. Zircon SHRIMP U–Pb age of Late Jurassic OIB-type volcanic rocks from the Tethyan Himalaya: Constraints on the initial activity time of the Kerguelen mantle plume. Acta Geochim.; 2018; 37, pp. 441-455. [DOI: https://dx.doi.org/10.1007/s11631-017-0239-2]

6. Tian, Y.; Gong, J.; Chen, H.; Guo, L.; Xu, Q.; Chen, L.; Lin, X.; Cheng, X.; Yang, R.; Zhao, L. et al. Early Cretaceous bimodal magmatism in the eastern Tethyan Himalayas, Tibet: Indicative of records on precursory continental rifting and initial breakup of eastern Gondwana. Lithos; 2018; 324, pp. 699-715. [DOI: https://dx.doi.org/10.1016/j.lithos.2018.12.001]

7. Zeng, Y.; Xu, J.; Chen, J.; Wang, B.; Huang, F.; Yu, H.; Chen, X.; Zhao, P. Breakup of Eastern Gondwana as Inferred from the Lower Cretaceous Charong Dolerites in the Central Tethyan Himalaya, Southern Tibet. Palaeogeogr. Palaeoclimatol. Palaeoecol.; 2019; 515, pp. 70-82. [DOI: https://dx.doi.org/10.1016/j.palaeo.2017.10.010]

8. Peng, W.; Yang, T.; Shi, Y.; Bian, W.; Ma, Y.; Wang, S.; Jiao, X.; Ma, J.; Kang, Y.; Zhang, S. et al. Role of the Kerguelen Mantle Plume in Breakup of Eastern Gondwana: Evidence from Early Cretaceous Volcanic Rocks in the Eastern Tethyan Himalaya. Palaeogeogr. Palaeoclimatol. Palaeoecol.; 2022; 588, 110823. [DOI: https://dx.doi.org/10.1016/j.palaeo.2021.110823]

9. Liang, J.; Bian, W.; Jiao, X.; Peng, W.; Ma, J.; Wang, S.; Ma, Y.; Zhang, S.; Wu, H.; Li, H. et al. Geochronological Results from the Zhela Formation Volcanics of the Tethyan Himalaya and Their Implications for the Breakup of Eastern Gondwana. Sci. Rep.; 2023; 13, 20035. [DOI: https://dx.doi.org/10.1038/s41598-023-47268-5]

10. Zhang, Z.; Li, G.; He, X.; Zhang, L.; Dong, S.; Liang, W.; Han, S. The evolution of Kerguelen mantle plume and breakup of eastern Gondwana: New insights from multistage Cretaceous magmatism in the Tethyan Himalaya. Gondwana Res.; 2023; 119, pp. 68-85. [DOI: https://dx.doi.org/10.1016/j.gr.2023.03.009]

11. Meng, F.; Meng, Y.; Wei, Y. Identification of the EMII-type Mantle End-member of the Kerguelen Mantle Plume and Its Implications for the Initial Breakup of Eastern Gondwana: New Constraints from the Cretaceous Diabase in the Tethyan Himalaya. J. Asian Earth Sci.; 2024; 272, 106250. [DOI: https://dx.doi.org/10.1016/j.jseaes.2024.106250]

12. Duncan, R.A.; Hooper, P.R.; Rehacek, J.; Marsh, J.S.; Duncan, A.R. The timing and duration of the Karoo igneous event, southern Gondwana. J. Geophys. Res. Solid Earth; 1997; 102, pp. 18127-18138. [DOI: https://dx.doi.org/10.1029/97JB00972]

13. Pirajno, F.; Santosh, M. Mantle plumes, supercontinents, intracontinental rifting and mineral systems. Precambrian Res.; 2015; 259, pp. 243-261. [DOI: https://dx.doi.org/10.1016/j.precamres.2014.12.016]

14. Müller, R.D.; Royer, J.Y.; Lawver, L.A. Revised plate motions relative to the hotspots from combined Atlantic and Indian Ocean hotspot tracks. Geology; 1993; 21, pp. 275-278. [DOI: https://dx.doi.org/10.1130/0091-7613(1993)021<0275:RPMRTT>2.3.CO;2]

15. Bryan, S.E.; Constantine, A.E.; Stephens, C.J.; Ewart, A.; Schön, R.W.; Parianos, J. Early Cretaceous volcano-sedimentary successions along the eastern Australian continental margin: Implications for the break-up of eastern Gondwana. Earth Planet. Sci. Lett.; 1997; 153, pp. 85-102. [DOI: https://dx.doi.org/10.1016/S0012-821X(97)00124-6]

16. Zhu, D.; Chung, S.; Mo, X.; Zhao, Z.; Niu, Y.; Song, B.; Yang, Y. The 132 Ma Comei-Bunbury Large Igneous Province: Remnants Identified in Present-day Southeastern Tibet and Southwestern Australia. Geology; 2009; 37, pp. 583-586. [DOI: https://dx.doi.org/10.1130/G30001A.1]

17. Wei, Y.; Liang, W.; Shang, Y.; Zhang, B.; Pan, W. Petrogenesis and Tectonic Implications of ~130 Ma Diabase Dikes in the Western Tethyan Himalaya (Western Tibet). J. Asian Earth Sci.; 2017; 143, pp. 236-248. [DOI: https://dx.doi.org/10.1016/j.jseaes.2017.04.008]

18. Wang, Z.; Shi, Y.; Yang, T.; Lawford Anderson, j.; Zhao, J.; Peng, X.; Hou, C.; Sun, H.; Bian, W.; Ma, Y. Early Activity of the Kerguelen Mantle Plume: Geochronology, Geochemistry and Sr-Nd-Pb Isotopes of Mafic Dykes and Sills from the Tethyan Himalaya. Int. Geol. Rev.; 2022; 65, pp. 512-526. [DOI: https://dx.doi.org/10.1080/00206814.2022.2053888]

19. Jiang, S.; Nie, F.; Hu, P.; Liu, Y.; Lai, X. Geochemical Characteristics of Mafic Dykes in Southern Tibet. Acta Geol. Sin.; 2007; 81, pp. 60-71. (In Chinese)

20. Tong, J.; Liu, J.; Zhong, H.; Xia, J.; Lu, R.; Li, Y. Zircon U-Pb Dating, Geochemistry, and Tectonic Significance of Mafic Dykes in the Luozha Area of Southern Tibet. Geol. Bull. China; 2007; 26, pp. 1654-1664. (In Chinese)

21. Qiu, B.; Zhu, D.; Zhao, Z.; Wang, L. Western Extension and Significance of the Cuomei Large Igneous Province in Southern Tibet. Acta Petrol. Sin.; 2010; 26, pp. 2207-2216. (In Chinese)

22. Xia, Y.; Zhu, D.; Zhao, Z.; Wang, Q.; Yuan, S.; Chen, Y.; Mo, X. Whole-Rock Geochemistry and Zircon Hf Isotopes of OIB-Type Mafic Rocks in the Cuomei Large Igneous Province, Southeastern Tibet. Acta Petrol. Sin.; 2012; 28, pp. 1588-1602. (In Chinese)

23. Zhu, D.; Xia, Y.; Qiu, B.; Wang, Q.; Zhao, Z. Why Propose the Early Cretaceous Cuomei Large Igneous Province in Southeastern Tibet?. Acta Petrol. Sin.; 2013; 29, pp. 3659-3670. (In Chinese)

24. Ding, F.; Gao, J.; Xu, K. Geochemical and Geochronological Characteristics of Mafic Intrusions in the Rongbu Area of Southern Tibet and Their Geological Significance. Acta Petrol. Sin.; 2020; 36, pp. 391-408. (In Chinese)

25. Zhou, Q.; Liu, Z.; Lai, Y.; Wang, G.; Liao, Z. Petrogenesis of Mafic and Felsic Rocks from the Comei Large Igneous Province, South Tibet: Implications for the Initial Activity of the Kerguelen Plume. Geol. Soc. Am. Bull.; 2018; 130, pp. 811-824. [DOI: https://dx.doi.org/10.1130/B31653.1]

26. Olierook, H.K.; Jiang, Q.; Jourdan, F.; Chiaradia, M. Greater Kerguelen Large Igneous Province Reveals No Role for Kerguelen Mantle Plume in the Continental Breakup of Eastern Gondwana. Earth Planet. Sci. Lett.; 2019; 511, pp. 244-255. [DOI: https://dx.doi.org/10.1016/j.epsl.2019.01.037]

27. Tian, Y. Early Cretaceous Bimodal Magmatic Activity in the Eastern Tethyan Himalayas: A Marker for Continental Rifting and the Initial Breakup of East Gondwana. Lithos; 2019; 324–325, pp. 699-715. (In Chinese) [DOI: https://dx.doi.org/10.1016/j.lithos.2018.12.001]

28. Lian, D.; Yang, J.; Liu, F.; Wu, Q.; Xu, Z. Fingerprints of the Kerguelen Mantle Plume in Southern Tibet: Evidence from Early Cretaceous Magmatism in the Tethyan Himalaya. J. Geol.; 2021; 129, pp. 207-231. [DOI: https://dx.doi.org/10.1086/714173]

29. Yang, C.; Li, Y.; Zhang, J.; Zhang, G. Age, Genesis, and Tectonic Significance of the Karlong Gabbro in the Langkazi Area of Southern Tibet. Geol. Sci.; 2022; 57, pp. 399-414. (In Chinese)

30. Liu, Z.; Zhou, Q.; Lai, Y.; Qing, C.; Li, X.; Wu, J.; Xia, X. Petrogenesis of the Early Cretaceous Laguila Bimodal Intrusive Rocks from the Tethyan Himalaya: Implications for the Break-up of Eastern Gondwana. Lithos; 2015; 236–237, pp. 190-202. [DOI: https://dx.doi.org/10.1016/j.lithos.2015.09.006]

31. Coffin, M.F.; Pringle, M.S.; Duncan, R.A.; Gladczenko, T.P.; Storey, M.; Müller, R.D. Kerguelen Hotspot Magma Output since 130 Ma. J. Petrol.; 2002; 43, pp. 1121-1139. [DOI: https://dx.doi.org/10.1093/petrology/43.7.1121]

32. Ingle, S.; Weis, D.; Doucet, S.; Mattielli, N. Hf isotope constraints on mantle sources and shallow-level contaminants during Kerguelen hot spot activity since ~120 Ma. Geochem. Geophys. Geosystems; 2003; 4, pp. 767-782. [DOI: https://dx.doi.org/10.1029/2002GC000482]

33. Liu, Z.; Kang, Z.; Yang, F.; Qiao, Z.; Ma, C.; Shen, Y. Geochronology, Geochemistry, and Petrogenesis of Volcanic Rocks from the Eastern Tethyan Himalayas. Geol. Explor. J.; 2024; 39, pp. 219-231. (In Chinese)

34. Olierook, H.K.; Jourdan, F.; Merle, R.E.; Timms, N.E.; Kusznir, N.; Muhling, J.R. Bunbury Basalt: Gondwana breakup products or earliest vestiges of the Kerguelen mantle plume?. Earth Planet. Sci. Lett.; 2016; 440, pp. 20-32. [DOI: https://dx.doi.org/10.1016/j.epsl.2016.02.008]

35. Yin, A.; Harrison, T.M. Geologic Evolution of the Himalayan-Tibetan Orogen. Annu. Rev. Earth Planet. Sci.; 2000; 28, pp. 211-280. [DOI: https://dx.doi.org/10.1146/annurev.earth.28.1.211]

36. Hou, C. Zircon Geochronology, Geochemistry, and Geological Significance of the Lacon Formation Volcanic Rocks in the Tethyan Himalayas. Master’s Thesis.; China University of Geosciences: Beijing, China, 2017; (In Chinese)

37. Boger, S.D. Antarctic Geology: Extending a Gondwanan Perspective. Gondwana Res.; 2011; 19, pp. 548-566.

38. Merdith, A.S.; Collins, A.S.; Williams, S.E.; Pisarevsky, S.; Foden, J.D.; Archibald, D.B.; Blades, M.L.; Alessio, B.L.; Armistead, S.; Plavsa, D. et al. A Full-Plate Global Reconstruction of the Neoproterozoic. Gondwana Res.; 2017; 50, pp. 84-134. [DOI: https://dx.doi.org/10.1016/j.gr.2017.04.001]

39. Huang, Y.; Liang, W.; Zhang, L.; Li, G.; Huang, C.; Xia, X.; Dong, S.; Wu, J. Initial Rifting of the Tethyan Himalayas-Indo-Asian Crust: Evidence from Late Cretaceous OIB-Type Basalts in Southern Tibet. Earth Sci.; 2018; 43, pp. 2651-2663. (In Chinese)

40. Pan, G.; Wang, L.; Yin, F.; Geng, Q.; Li, G.; Zhu, D. Review, Progress, and Prospects of the Formation and Evolution of the Tibetan Plateau. Sediment. Tethys Geol.; 2022; 42, pp. 151-175. (In Chinese)

41. Liu, Y.; Gao, S.; Hu, Z.; Gao, C.; Zong, K.; Wang, D. Continental and Oceanic Crust Recycling-induced Melt-Peridotite Interactions in the Trans-North China Orogen: U-Pb Dating, Hf Isotopes and Trace Elements in Zircons from Mantle Xenoliths. J. Petrol.; 2009; 51, pp. 392-399. [DOI: https://dx.doi.org/10.1093/petrology/egp082]

42. Andersen, T. Correction of common lead in U-Pb analyses that do not report ²⁰⁴Pb. Chem. Geol.; 2002; 192, pp. 59-79. [DOI: https://dx.doi.org/10.1016/S0009-2541(02)00195-X]

43. Ludwig, K. A Plotting and Regression Program for Radiogenic Isotope Data; Open-File Report; 91-445 U.S. Geological Survey: Sunrise Valley Drive Reston, VA, USA, 1991.

44. Wang, X.; Zhou, J.; Qiu, J.; Zhang, W.; LIU, X.; Zhang, G. LA-ICP-MS U-Pb zircon geochronology of the Neoproterozoic igneous rocks from northern Guangxi, south China: Implications for tectonic evolution. Precambrian Res.; 2006; 145, pp. 111-130. [DOI: https://dx.doi.org/10.1016/j.precamres.2005.11.014]

45. Yuan, H.; Gao, S.; Liu, X.; Li, H.; Günther, D.; Wu, F. Accurate U-Pb Age and Trace Element Determinations of Zircon by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry. Geostand. Geoanal. Res.; 2004; 28, pp. 353-370. [DOI: https://dx.doi.org/10.1111/j.1751-908X.2004.tb00755.x]

46. Li, X.Y.; Wang, H.Y.; Fan, G.Q. Accurate determination of more than 40 trace elements in rock samples by ICP-MS. Priv. Sci. Technol.; 2016; 18.(In Chinese)

47. Liang, X.; Wei, G.; Li, X.; Liu, Y. Precise Determination of ¹⁴³Nd/¹⁴⁴Nd and Sm/Nd Ratios Using MC-ICPMS. Geochimica; 2003; 32, pp. 91-96. (In Chinese)

48. Griffin, W.L.; Pearson, N.J.; Belousova, E.; Jackson, S.E.; Achterbergh, V.E.; O’Reilly, S.Y.; Shee, S.R. The Hf isotope composition of cratonic mantle: LAM-MC-ICP-MS analysis of zircon megacrysts in kimberlites. Geochim. Cosmochim. Acta; 2000; 64, pp. 133-147. [DOI: https://dx.doi.org/10.1016/S0016-7037(99)00343-9]

49. Pearce, J. A User’s Guide to Basalt Discrimination Diagrams. Trace Element Geochemistry of Volcanic Rocks: Applications for Massive Sulphide Exploration; Geological Association of Canada Short Course Notes; Geological Association of Canada: Toronto, ON, Canada, 1996; Volume 12, pp. 79-113.

50. Maniar, P.D.; Piccoli, P.M. Tectonic Discrimination of Granitoids. Geol. Soc. Am. Bull.; 1989; 101, pp. 635-643. [DOI: https://dx.doi.org/10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2]

51. Wright, J.B. A Simple Alkalinity Ratio and Its Application to Questions of Non-Orogenic Granite Genesis. Geol. Mag.; 1969; 106, pp. 370-384. [DOI: https://dx.doi.org/10.1017/S0016756800058222]

52. Winchester, J.A.; Floyd, P.A. Geochemical discrimination of different magma series and their differentiation products using immobile elements. Chem. Geol.; 1977; 20, pp. 325-343. [DOI: https://dx.doi.org/10.1016/0009-2541(77)90057-2]

53. Sun, S.S.; McDonough, W.F. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. Geol. Soc. Lond. Spec. Publ.; 1989; 42, pp. 313-345. [DOI: https://dx.doi.org/10.1144/GSL.SP.1989.042.01.19]

54. Deng, J. Paleomagnetic and Zircon U-Pb Geochronology of the Zhela Formation Volcanic Rocks in the Tethyan Himalayas. Master’s Thesis; China University of Geosciences: Beijing, China, 2021; (In Chinese)

55. Ludden, J.N.; Thompson, G. Behaviour of Rare Earth Elements During Submarine Weathering of Tholeiitic Basalt. Nature; 1978; 274, pp. 147-149. [DOI: https://dx.doi.org/10.1038/274147a0]

56. Bienvenu, P.; Bougault, H.; Joron, J.L.; Treuil, M.; Dmitriev, L. MORB Alteration: Rare-Earth Element/Non-Rare-Earth Hygromagmaphile Element Fractionation. Chem. Geol.; 1990; 82, pp. 1-14. [DOI: https://dx.doi.org/10.1016/0009-2541(90)90070-N]

57. Staudigel, H.; Plank, T.; White, B.; Schmincke, H. Geochemical Fluxes During Seafloor Alteration of the Basaltic Upper Oceanic Crust: DSDP Sites 417 and 418. Subduction: Top to Bottom; Geophysical Monograph Series; Advancing Earth and Space Sciences: Washington, DC, USA, 2013; Volume 96, pp. 19-38.

58. Wilson, M.; Hyman, U. Igneous Petrogenesis; Springer: Dordrecht, The Netherlands, 1997.

59. Lassiter, J.C.; Depaolo, D.J. Plume/Lithosphere Interaction in the Generation of Continental and Oceanic Flood Basalts: Chemical and Isotopic Constraints; American Geophysical Union (AGU): Washington, DC, USA, 2013.

60. Weaver, B.L. The origin of ocean island basalt end-member compositions: Trace element and isotopic constraints. Earth Planet. Sci. Lett.; 1991; 104, pp. 381-397. [DOI: https://dx.doi.org/10.1016/0012-821X(91)90217-6]

61. Rudnick, R.L.; Gao, S. Composition of the Continental Crust. Treatise Geochem.; 2003; 3, pp. 1-64.

62. Taylor, S.R.; Mclennan, S.M. The continental crust: Its composition and evolution. J. Geol.; 1985; 94, pp. 57-72. [DOI: https://dx.doi.org/10.1086/629067]

63. McDonough, W.F. Contraints on the composition of the continental lithospheric mantle. Earth Planet. Sci. Lett.; 1990; 101, pp. 1-18. [DOI: https://dx.doi.org/10.1016/0012-821X(90)90119-I]

64. Shuguang, L. Ba-Nb-Th-La diagrams used to identify tectonic environments of ophiolite. Acta Petrol. Sin.; 1993; 9, 146.

65. Meschede, M. A method of discriminating between different type of mid-ocean ridge basalts and continental tholeiites with the Nb-Zr-Y diagram. Chem. Geol.; 1998; 56, pp. 207-218. [DOI: https://dx.doi.org/10.1016/0009-2541(86)90004-5]

66. Pearce, J.A.; Norry, M.J. Petrogenetic Implications of Ti, Zr, Y, and Nb Variations in Volcanic Rocks. Contrib. Miner. Petrol.; 1979; 69, pp. 33-47. [DOI: https://dx.doi.org/10.1007/BF00375192]

67. Pearce, J. Trace element characteristics of lavas from destructive plate boundaries. Andesites: Orogenic Andesites and Related Rocks; John and Wiley and Sons: Hoboken, NJ, USA, 1982.

68. Wang, Y.; Zhang, C.; Xiu, S. Th/Hf-Ta/Hf Discrimination Diagram for Tectonic Environments of Basalt Formation. Acta Petrol. Sin.; 2001; 17, pp. 413-421. (In Chinese)

69. Dong, L.; Li, G.; Zhang, Z.; Huang, Y.; Cao, H. Zircon U-Pb Ages, Geochemistry, and Sr-Nd-Pb Isotopes of Diabase Dykes from the Xuesha Area, Southern Tibet. J. Chengdu Univ. Technol. (Sci. Technol. Ed.); 2019; 46, pp. 734-745. (In Chinese)

70. Yang, G.; Li, Y.; Xiao, W.; Tong, L. OIB-type rocks within West Junggar ophiolitic melanges: Evidence for the accretion of seamounts. Earth-Sci. Rev.; 2015; 150, pp. 477-496. [DOI: https://dx.doi.org/10.1016/j.earscirev.2015.09.002]

71. Liu, H.; Wang, Y.; Cawood, P.A.; Guo, X. Episodic slab rollback and back-arc extension in the Yunnan-Burma region: Insights from Cretaceous Nb-enriched and oceanic-island basalt-like mafic rocks. Geol. Soc. Am. Bull.; 2017; 129, B31604. [DOI: https://dx.doi.org/10.1130/B31604.1]

72. Ren, C. Geological Geochronology, Geochemistry, and Geological Significance of Early Cretaceous Magmatic Rocks in Cuomei, Tibet. Bachelor’s Thesis; China University of Geosciences: Beijing, China, 2015; (In Chinese)

73. Zhu, D.; Pan, G.; Mo, X.; Liao, Z.; Jiang, X.; Wang, L.; Zhao, Z. Geochemistry and Petrogenesis of the Basalts from the Sanxiu Formation in the Central Tethyan Himalayas. Geochimica; 2005; pp. 7-19. (In Chinese)

74. Ma, Y.; Yang, T.; Bian, W.; Zhang, S.; Wu, H.; Li, H. Early Cretaceous paleomagnetic and geochronologic results from the Tethyan Himalaya: Insights into the Neotethyan paleogeography and the India-Asia collision. Sci. Rep.; 2016; 6, 21605. [DOI: https://dx.doi.org/10.1038/srep21605]

75. Davis, J.K.; Lawver, L.A.; Norton, I.O.; Gahagan, L.M. New Somali Basin magnetic anomalies and a plate model for the early Indian Ocean. Gondwana Res.; 2016; 34, pp. 16-28. [DOI: https://dx.doi.org/10.1016/j.gr.2016.02.010]

76. Storey, B.C. The role of mantle plumes in continental breakup: Case histories from Gondwanaland. Nature; 1995; 377, pp. 301-308. [DOI: https://dx.doi.org/10.1038/377301a0]

77. Gibbons, A.D. The breakup of East Gondwana: Assimilating constraints from Cretaceous ocean basins around India into a best-fit tectonic model. J. Geophys. Res. Solid Earth; 2013; 118, pp. 808-822. [DOI: https://dx.doi.org/10.1002/jgrb.50079]

78. Jiang, Q.; Jourdan, F.; Olierook, H.K. Longest continuously erupting large igneous province driven by plume-ridge interaction. Geology; 2020; 49, pp. 206-210. [DOI: https://dx.doi.org/10.1130/G47850.1]

79. Torsvik, T.H.; Van der Voo, R.; Preeden, U.; Niocaill, C.M.; Steinberger, B.; Doubrovine, P.V.; Hinsbergen, D.; Domeier, M.; Gaina, C.; Tohver, E. Phanerozoic polar wander, palaeogeography and dynamics. Earth-Sci. Rev.; 2012; 114, pp. 325-368. [DOI: https://dx.doi.org/10.1016/j.earscirev.2012.06.007]