Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


Recent research on the water content of large igneous provinces (LIPs) has revealed that water has a significant impact on the formation of LIPs. However, most studies focus on the water content of mafic-ultramafic rocks, while relatively little attention has been paid to the water content of continental flood basalts (CFB), which form the major part of LIPs and are characterized by huge volumes (> 1 × 105 km3) and short eruption times. Here, we determined water contents of clinopyroxene crystals from the Akesu diabase, which is co-genetic with flood basalts of the Tarim LIP in China. Based on these measurements, we obtained a water content of higher than 1.23 ± 0.49 wt.% for the parental magma to the Tarim CFB and a minimum water content of 1230 ± 490 ppm for the mantle source, thus indicating the presence of a hydrous mantle plume. Combined with previous studies, our results suggest that water plays a key role in the formation of the Tarim LIP. Additionally, the whole-rock compositions of the Akesu diabase indicate a contribution of pyroxenite in the mantle source. This is consistent with a model, in which water was brought into the Tarim mantle plume by a subducted oceanic plate that entered the deep mantle.

Free full text 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2024; 14: 9514.
Published online 2024 Apr 25. https://doi.org/10.1038/s41598-024-60213-4
PMCID: PMC11045731
PMID: 38664514

Hydrous mantle plume promoted the generation of continental flood basalts in the Tarim large igneous province

Associated Data

Supplementary Materials
Data Availability Statement

Abstract

Recent research on the water content of large igneous provinces (LIPs) has revealed that water has a significant impact on the formation of LIPs. However, most studies focus on the water content of mafic–ultramafic rocks, while relatively little attention has been paid to the water content of continental flood basalts (CFB), which form the major part of LIPs and are characterized by huge volumes (> 1 × 105 km3) and short eruption times. Here, we determined water contents of clinopyroxene crystals from the Akesu diabase, which is co-genetic with flood basalts of the Tarim LIP in China. Based on these measurements, we obtained a water content of higher than 1.23 ± 0.49 wt.% for the parental magma to the Tarim CFB and a minimum water content of 1230 ± 490 ppm for the mantle source, thus indicating the presence of a hydrous mantle plume. Combined with previous studies, our results suggest that water plays a key role in the formation of the Tarim LIP. Additionally, the whole-rock compositions of the Akesu diabase indicate a contribution of pyroxenite in the mantle source. This is consistent with a model, in which water was brought into the Tarim mantle plume by a subducted oceanic plate that entered the deep mantle.

Subject terms: Geochemistry, Volcanology

Introduction

Flood basalts are the main constituents of large igneous provinces (LIPs). Their huge magma volumes (> 1 × 105 km3) but short eruption times (usually < 1–3 Ma)1,2 imply unique geodynamic processes, and therefore the origin of LIPs has been a hot topic in Earth sciences2.

The view according to which the formation of LIPs is related to mantle plumes originated from the core-mantle boundary is widely accepted24. Previous studies have found that these mantle plumes generally have an abnormally high temperature4,5, indicating that temperature is one of the important factors to the formation of LIPs. In addition, experimental petrology has demonstrated that water lowers the solidus temperature of rocks68 and reduces the viscosity of magmas9,10, which may contribute to the formation of LIPs.

In recent years, studies of the water content in LIP magmas have gradually increased1123. Based on the small-scale mafic–ultramafic rocks (e.g., picrites, meimechites, komatiites and mafic intrusive rocks), they have constrained the water content in the LIPs, including the Emeishan LIP (2.71–3.73 wt.%)11,12, the Siberian LIP (0.25–3.88 wt.%)13,14, the Karoo LIP (1–2 wt.%)15, the Caribbean LIP (0.25–0.70 wt.%)1618, and the Tarim LIP (4.82 ± 1.00 wt.%)19. In addition, Stefano et al.20, Cabato et al.21 and Choudhary et al.22 measured the water content of melt inclusions in Snake River Plain basalts (SRPB), Columbia River basalts (CRB) and Deccan Trap basalts, respectively, and found that these basalts likely have a high water contents as well (up to 2.4 wt.%, 3.3 wt.% and 2.05 wt.%, respectively). Based on these studies above, Liu et al.11 estimated that the water content of the mantle source in the LIPs ranges from 900 to > 6000 ppm, and the H2O/Ce ratios range from 160– 400 to > 2000, which are significantly higher than those of MORBs (50–250 ppm, 150–210)24,25 or most OIBs (300–1000 ppm, < 250)26,27. Thus, they proposed that water is an important factor in the formation of LIPs1123. However, most of the water content data are from small-scale mafic–ultramafic rocks, and the data from continental flood basalts (CFB) are still relatively scares. Besides, CFB often undergo intense crystal fractionation and the host minerals of melt inclusions are not the earliest-crystallized2022 phases in previous studies, making it difficult to directly use the measured water contents to obtain information on the mantle sources. Hence, the water contents in CFB and the effects of water on CFB are not yet well understood.

The Tarim LIP is a typical LIP of the Phanerozoic Eon. It is located within the Tarim Basin in northwest China with a residual area of approximately 250,000 km2 and a maximum thickness of up to 800 m28,29 (Fig. 1a, b). Except for the small-volume kimberlites emplaced at ~ 300 Ma30, the Tarim LIP can be divided into two stages: the main stage (290–288 Ma) produced large-scale CFB28,29 and the later stage (284–278 Ma) small-scale mafic–ultramafic intrusive rocks28,29. Liu et al.11 calculated a high mantle potential temperature for picrites, indicating that the formation of the Tarim LIP is related to a mantle plume. Xia et al.19 found that late-stage intrusive rocks have a high water content. In addition, based on numerical simulation, Liu and Leng31 also proposed that the formation of the Tarim LIP is related to a water-rich mantle plume. However, the water content of the CFB in the Tarim LIP is still unclear, and therefore this work focuses on determining water contents of the Tarim CFB.

An external file that holds a picture, illustration, etc.
Object name is 41598_2024_60213_Fig1_HTML.jpg

(a) Map showing the location of the Tarim LIP (Modified after Zhou et al.32). (b) Simplified geologic map of the Tarim LIP (modified after Xu et al.29) showing the distribution of CFB and mafic–ultramafic intrusive rocks. (c) Geological sketch map of the Akesu section (modified after Turner66).

The Tarim CFB generally have low MgO contents28,29,32,33 (Supplementary Fig. S1), indicating that they have undergone intense fractional crystallization, which means that their water contents do not represent those of the primary melts. Recently, Wu et al.34 reported the presence of sheet-like mafic intrusive rocks (diabases) with a thickness of 6–8 m in the Akesu area (Fig. 1c). These rocks intruded into the Lower Cambrian Sugetbark Formation, forming a discordant body. Geochronological studies indicate that the formation age of these rocks is between 292 and 290 Ma34,35, which is consistent with the age of the Tarim CFB rather than that of the mafic–ultramafic intrusive rocks belonging to the later magmatic stage28,29,32,33,36. In addition, the trace element characteristics and Sr–Nd isotope compositions of the Akesu diabase are similar to those of the Tarim CFB32. These results confirm that the Akesu diabase is co-genetic with the Tarim CFB. It is noteworthy that the MgO content of the Akesu diabase (up to 8.71 wt.%)35 is generally higher than that of other Tarim CFB (2.05–6.99 wt.%)28,29,32,33 (Supplementary Fig. S1), indicating that the diabase is less evolved and can better reflect the information on the mantle source of the Tarim CFB.

In this study, six samples of the Akesu diabase with large clinopyroxene phenocrysts (18SK3-2, 18SK3-4, 13-AKS-19a, 14-AKS-19a, 11-AKS-19a and 11-AKS-19b) from previous studies34,35 were selected (Supplementary Fig. S1 and S2). The diabase generally exhibits a low porosity, appearing as dense blocks with localized spherical weathering. Petrographic microscope shows that these diabases are ophitic texture, but commonly contain large clinopyroxene phenocrysts (Supplementary Fig. S3). The matrix is primarily composed of plagioclase and pyroxene microlites, and volcanic glass34,35. Water contents of clinopyroxene phenocrysts were measured using unpolarized infrared spectra and major and trace element contents of clinopyroxene phenocrysts of these samples were determined with an electron microprobe and LA-ICP-MS instrument, respectively. Based on mass balance calculations, the water content and H2O/Ce ratio in the mantle source were estimated and the impact of water on the formation of LIPs is discussed.

Results

Water content and chemical composition of clinopyroxene in the Akesu diabase

Major and trace element compositions of clinopyroxene phenocrysts are shown in Supplementary Table S1. The Mg# (Mg# = molar ratio of 100 × Mg2+/(Mg2+ + Fe2+)) of clinopyroxene phenocrysts ranges from 65.6 to 74.0, which are too low to be in equilibrium with the melt of the Akesu diabase (Supplementary Fig. S4). Based on the partition coefficients of trace elements between clinopyroxene and melt37, the trace element composition of the equilibrated melt was calculated. As shown in Supplementary Fig. S5, the trace elements pattern of the equilibrated melt is consistent with that of Akesu diabase, which features indicate that the clinopyroxene phenocrysts were crystallized from the magma of their host rocks37.

In the OH stretching vibration region (3000–3800 cm−1), the clinopyroxene in the Akesu diabase exhibits three main bands at ~ 3640 cm−1, ~ 3530 cm−1, and ~ 3460 cm−1 (Fig. 2a), of which the band at 3640 cm−1 is strongest. These results are consistent with data from clinopyroxene phenocrysts in basalts reported previously38,39. The maximum linear absorbance band for all clinopyroxenes is less than 0.15, meaning that the unpolarized infrared spectra measurement is reliable with an error < 30%40. Based on the Beer-Lambert law, the calculated water content of clinopyroxene grains in the Akesu diabase ranges from 20.7 to 253.0 ppm (see Supplementary Table S1 and Supplementary Fig. S6. The detailed calculation process is reported in the Supplementary Text).

An external file that holds a picture, illustration, etc.
Object name is 41598_2024_60213_Fig2_HTML.jpg

(a) Representative IR spectra of OH in clinopyroxene phenocrysts from the Akesu diabase. (b) Representative IR of profiles for clinopyroxene phenocrysts from the Akesu diabase. The vertical axis is normalized absorbance to 1 cm and the baseline has been corrected. The dashed lines mark the position of individual OH bands. The insert shows spectra without baseline correction.

The water content of the melt in equilibrium with the clinopyroxene was calculated using mass balance and the coefficient DH2Ocpx-melt using the following equation 41:

DH2Ocpx-melt=Exp (-4.2+6.5×XAlIVcpx-XCacpx)
1

where XAlIVcpx and XCacpx are the stoichiometric numbers of four-coordinated Al and Ca in the chemical equation of clinopyroxene phenocrysts (under 6 oxygen atoms). The calculated water content of the melt equilibrated with clinopyroxene phenocrysts ranges from 0.17 to 1.87 wt.%.

For clinopyroxene phenocrysts with the maximum linear absorbance band less than 0.3, the maximum deviation in measuring of a single clinopyroxene phenocryst is less than 19%40,42. Taking into account an error of 10% for the I, the error of calculated the water content of a single clinopyroxene phenocryst by unpolarized light is < 30%,40,42,43. Considering the typical uncertainty of DH2Ocpx-melt 41(~ 10%), the total uncertainty of the calculated water content in the melt is less than 40%40,42,43.

Discussion

The effects of alteration on the Akesu diabase

In surface environments, rocks may undergo alteration, leading to the changes in the water content. So, it is necessary to assess the effects of alteration, on the Akesu diabase. before evaluating its water content.

The Akesu diabase exhibits high levels of loss on ignition (LOI 3.2–5.3 wt.%)35, indicating a relatively strong effect of alteration. In particular, some samples with high LOI values have anomalous contents of fluid-mobile elements (e.g., K, Rb; Supplementary Fig. S7 and S8), indicating that their chemical composition has been strongly affected by alteration. However, the K and Rb contents of the samples selected for this study are relatively consistently low (Supplementary Fig. S7), suggesting that the effect of alteration on this samples is limited. Furthermore, the infrared spectroscopy is highly sensitive to the hydroxy, can effectively identify different types of structural water and molecular water, allowing for the check whether minerals have undergone alteration or not. Therefore, we were able to select fresh clinopyroxene grains for analysis of the water content and chemical composition to avoid the influence of alteration.

Crustal contamination and fractional crystallization

During magma evolution, the water content of melt may be influenced by crustal contamination and fractional crystallization. Compared to the mantle, the continental crust is generally enriched in U and Pb, but depleted in Ce and Nb44. Thus, when magma is contaminated by the material from continental crust, a positive correlation in the ratios of Nb/U to Ce/Pb or Nb/U to Nb can be identified. However, for samples with lesser influence from alteration, the lack of a relationship between Nb/U and Ce/Pb or Nb/U and Nb (Supplementary Fig. S9) suggests that crustal contamination did not significantly affect the Akesu diabase.

While the MgO content of the Akesu diabase is higher compared to other Tarim CFB, it is still relatively low in general, which means that the Akesu diabase has undergone intense fractional crystallization. As shown in Supplementary Fig. S10, the Akesu diabase has a low concentration of Ni and Cr, and the Ni and Cr display a positive correlation with MgO, respectively. These phenomenon in Supplementary Fig. S10 suggests that Akesu diabase primarily underwent fractional crystallization of the olivine and clinopyroxene, as indicated by the enrichment of Ni in olivine and Cr in clinopyroxene. Moreover, if fractional crystallization of plagioclase takes place, the magma would exhibit negative anomalies in Eu due to the compatibility of Eu with plagioclase. However, the absence of an obviously Eu anomaly in the rare earth element pattern diagram of the Akesu diabase (Supplementary Fig. S8) indicates that the crystal differentiation of plagioclase is insignificant.

Estimation of the melt water content of the Akesu diabase

The calculated trace element composition of the melt in equilibrium with the clinopyroxene phenocrysts displays a similar pattern to the Akesu diabase (Supplementary Fig. S5c, d), indicating that these clinopyroxene phenocrysts crystallized from the Akesu diabase magma. This means that clinopyroxene grains can be utilized to infer the water content of the magma.

During crystallization, H is an incompatible element and consequently the water content of the residual melt would increase. To eliminate the effects of crystallization, the earliest-crystallized phenocrysts should be used to calculate the water content of the primitive melt. However, all clinopyroxene phenocrysts in the diabase have relatively low Mg# (65.6–74.0) as a result of strong crystal fractionation and hence the water content in the melt was modified by crystal fractionation.

As mentioned earlier, the water content of the Akesu diabase is influenced by fractional crystallization of olivine and clinopyroxene. In order to constrain the water content of the primary melt, we incrementally added olivine and clinopyroxene to the melts in certain proportions, with 1% incremental steps. The calculation was stopped when the melt Mg# reaches 72. The simulation process involved the design of three different scenarios. In the first scenario, olivine and clinopyroxene were introduced into the melt in a 1:1 ratio. In the second scenario, the ratio was adjusted to 2:1, with a higher proportion of olivine to clinopyroxene. The third scenario involved solely adding olivine to the melt. For the calculation of H2O content in melt, we employed specific H partition coefficients and Fe–Mg exchange coefficients. The H partition coefficient for olivine was set at 0.000245, while the Kd for olivine was determined as 0.3046. Additionally, we utilized a H partition coefficient of 0.015 (calculated from clinopyroxene in Akesu diabase) and Kd of 0.3346 for clinopyroxene.

The simulation results, as depicted in Supplementary Fig. S11, indicate an overestimation of the water content in the melt by 52.0% (1:1), 43.2% (2:1), and 30.8% (only olivine) in the three scenarios, respectively.

Considering that olivine and clinopyroxene are the primary crystallized phases during magma evolution in Akesu diabase, we utilized a 1:1 ratio of olivine to clinopyroxene to estimate the minimum water content. Then the water content of the equilibrated melts was determined to be 0.11 to 1.23 wt.%.

On the other hand, H profile analyses across clinopyroxene phenocrysts show that the water content at the rim is generally lower than that in the core (Fig. 2b). This indicates outward diffusion of H in clinopyroxene phenocrysts.

As shown in the Fig. 3a, calculated water content of the equilibrated melt varies widely and demonstrate a slight decrease with the evolution of magma, (especially for 18SK3-2). In theory, due to the strong incompatibility of hydrogen (H), the water content in magma tends to continuously increase during the process of fractional crystallization47. However, this is different from features observed in Fig. 3, suggesting the degassing of magma. Furthermore, the Akesu diabase exhibits a significant variation in water content, as depicted in Fig. 3a. This wide range of water content suggests the occurrence of degassing processes or the lack of homogenization in clinopyroxene phenocrysts.

An external file that holds a picture, illustration, etc.
Object name is 41598_2024_60213_Fig3_HTML.jpg

(a) Water contents of equilibrated melts and (b) calculated H2O/Ce ratios of equilibrated melts versus Mg# values of clinopyroxene phenocrysts in the Akesu diabase. The inserted schematic plot illustrates the variation in H2O-Mg# trends during the processes of crystallization and degassing (modified after Wade et al.47). The partition coefficient for water between clinopyroxene and equilibrated melt is from Eq. 10 in O’Leary et al.41. The partition coefficient for Ce between clinopyroxene and equilibrated melt is from Eqs. 304, 306, and 307 in Bédard48. The arrows show the effects of degassing47.

In order to minimize the influence of whole-rock composition, the H2O/Ce ratio of the equilibrated melt was calculated based on the major and trace element compositions of clinopyroxene phenocrysts. The partition coefficient of Ce was estimated by averaging values obtained from Eqs. 304, 306, and 307 in Bédard48, while the partition coefficient of H2O was determined using Eq. 10 in O'Leary et al.41. Figure 3b demonstrates a notable positive correlation between the H2O/Ce ratio and the Mg# of the samples. Previous studies have shown that H2O and Ce exhibit similar incompatibility in basaltic melts25,26. This suggests that there should be minimal changes in the H2O/Ce ratio during partial melting and fractional crystallization. When comparing the variations of 1/Ce and H2O/Ce, it is evident that the trend in 1/Ce shows much smaller changes compared to the trend in H2O/Ce (Supplementary Fig. S12). This indicates that the observed trend in Fig. 3b can be attributed to degassing rather than an enrichment of Ce in the melt.

According to the research conducted by Patkó et al.49, clinopyroxene found in peridotite samples from the northern Pannonian Basin exhibited infrared spectra where the dominant OH absorption band at 3630 cm−1 was not prominent (the ratio of absorption intensities at approximately 3630 cm−1 and ~ 3525 cm−1 wavenumbers was less than 1.2). They concluded that these clinopyroxene phenocrysts experienced equilibration under conditions of lower water activity, resulting in the loss of water content. However, in the case of Akesu diabase, the clinopyroxene phenocrysts displayed consistent features with other basalts40,42,43, showing a dominant OH absorption band at 3630 cm−1. This discrepancy indicates that the clinopyroxene phenocrysts in Akesu diabase did not undergo a similar process as observed in the peridotite samples from the northern Pannonian Basin49.

Consequently, we can conclude that the water content of the primitive melt is underestimated. Considering both the impact of crystal fractionation and water diffusion, it is likely that the water content of the primitive melt in the Tarim CFB was higher than 1.23 ± 0.49 wt.%.

The estimated water content of the Tarim CFB is similar to those of the other CFB (SRPB, 0.2–3.3 wt.%20; CRB, 0.6–4.24 wt.%21; Deccan flood basalt, 0.94–2.05 wt.%22 and Emeishan flood basalt, 2.71 ± 0.95 wt.%)23 and the small-scale mafic–ultramafic rocks in LIPs (0.25–4.83 wt.%)1118 and is significantly higher than that of OIBs (0.1–0.3 wt.%)50 and MORBs (0.3–1.0 wt.%)10.

High water content and H2O/Ce ratio in the mantle source of the Tarim CFB

The water content in the mantle source of the Tarim CFB can be calculated using batch melting and fractional melting models and the following equations:

Batch melting

C0=CL[F+1-F×D]
2

Fractional melting

C0=CLF[1-(1-F)1D]
3

where C0 is the water content in the mantle source, CL is the water content of the melt, F is the degree of partial melting, and D is the bulk partition coefficient of water between melt and the mantle minerals (D = 0.0088–0.013)7.

Based on the fractionation of rare earth elements, the degree of partial melting for the Akesu diabase magma was estimated to be approximately 10% (Fig. 4). Our calculation indicates that the minimum water content in the mantle source of the Tarim CFB is 1330 ± 530 ppm for the batch melting model and 1230 ± 490 for the fractional melting model (Supplementary Table S2).

An external file that holds a picture, illustration, etc.
Object name is 41598_2024_60213_Fig4_HTML.jpg

Sm/Yb versus La/Yb diagram showing partial melting degrees for the Akesu diabase magma and other Tarim CFB (modified after Cheng et al.54). The mantle source is assumed to be composed of carbonated peridotite and pyroxenite, with peridotite containing 55% olivine, 13% orthopyroxene, 15% clinopyroxene 15% garnet and 2% calcite. Pyroxenite is composed of clinopyroxene. The dashed line represents mixing curves between melts derived from peridotite and pyroxenite. The dashed line represents mixing curves between melts derived from peridotite and pyroxenite. The data symbols with dashed borders represent samples that have been strongly affected by alteration. The chemical compositions of the Akesu diabase samples selected for water content measurements and the other Akesu diabase samples are all from Cheng et al.35. The contents of La, Sm and Yb are from McDonough and Sun67, and mineral/melt partition coefficients for La, Sm and Yb are from Rollison68. The data of other Tarim CFB are from Keping basalts28,29,32,33,54.

As shown in Fig. 5a, the mantle source of other CFB, including SRPB (940–2257 ppm)20, CRB (2257–3869 ppm)21, and Emeishan basalt (1357–1579 ppm)23, exhibit similar water contents to that of the Akesu diabase. In the case of small-scale mafic–ultramafic rocks found in LIPs, the water content in their mantle source ranges from 750 to higher than 6000 ppm1119 (Fig. 5a). It is evident that the mantle source of both CFB and mafic–ultramafic rocks in LIPs contain much more water than those of MORBs (50–250 ppm)24,25 or OIB (300–1000 ppm)26,27.

An external file that holds a picture, illustration, etc.
Object name is 41598_2024_60213_Fig5_HTML.jpg

Comparison of (a) the water contents and (b) H2O/Ce ratios in the mantle source of LIPs and OIBs (modified after Liu et al.11). The purple boxes represent small-scale mafic–ultramafic rocks in LIPs, and the green boxes represent CFB. The ranges of Hawaii and other OIBs are from Bizimis and Peslier26. Data of LIPs are from the Siberian Trap13,14, Emeishan LIP11,12,23, Caribbean LIP1618, Karoo LIP15, Columbia River basalts (CRB)21, Snake River Plain basalts (SRPB)20 and Tarim LIP19,52.

As mentioned earlier, the H2O/Ce ratio in magma remains stable during processes of partial melting and fractional crystallization51. This indicates that H2O/Ce in the melt can provide information on the ratio in the mantle source. Our calculations reveal that the H2O/Ce ratios of the melt ranges from 30 to 500 (as illustrated in Fig. 3b). However, considering the effects of hydrogen diffusion, the H2O/Ce ratio of the mantle source likely exceeds 500.

According to Fig. 5b, the H2O/Ce ratios of the Akesu diabase are similar to those in other LIPs (150 to  > 2000)1121,23, and the majority of LIPs have higher H2O/Ce ratios than MORBs (150–210)24,25. The results indicate that the water content in the mantle source of the Tarim CFB is higher than 1230 ± 490 ppm and H2O/Ce ratio exceeds 500. When combined with previous data from late-stage mafic intrusive rocks (4820 ppm, 1043)19,52, it can be concluded that the water content and H2O/Ce ratio in the mantle source of the Tarim LIP, falling within the range of LIPs (Fig. 5), are significantly higher than those of MORBs and OIBs24,25.

The presence of water in the mantle source can lower the solidus of rocks, cause melting to occur at greater depths, and increase the volume of the generated magma68. This likely contributed to the formation of the large-scale CFB observed in the Tarim LIP. Additionally, water can decrease the density and viscosity of magma, which in turn can increase the rate of magma ascent, promoting the formation and upward movement of mantle plumes. These observations suggest that water plays a key role in the formation of the Tarim LIP.

The origin of the water in the Tarim CFB

As a part of the Tarim CFB, the Akesu diabase exhibits lower CaO and higher Zn/Fe ratios than melts derived from peridotite (Fig. 6), which suggests the existence of pyroxenite components in the mantle source. Previous studies5355 have shown that the Tarim continental flood basalts have lower Mg and higher Zn isotope ratios compared to the primitive mantle, suggesting the involvement of recycled oceanic plate in the mantle source. During the subduction of plates, basaltic oceanic plates, in the form of eclogite, can be subjected to partial melting and react with the peridotite mantle. This process can result in the formation of pyroxenite56. Taken together, these findings suggest that the mantle source of the Tarim CFB involved a subducted oceanic plate component.

An external file that holds a picture, illustration, etc.
Object name is 41598_2024_60213_Fig6_HTML.jpg

(a) CaO and (b) 10,000*Zn/Fe versus MgO diagram for bulk rock samples from the Akesu diabase and other Tarim CFB. The line separating pyroxenite and peridotite melts in (a) is from Herzberg and Asimow69. The ranges of 10,000*Zn/Fe in peridotite and pyroxenite melts in (b) are from Le Roux et al.70. The data sources for the Akesu diabase and other Tarim CFB are same as in Fig. 4.

Subducted oceanic plates play a crucial role in transporting water into the Earth's interior. During subduction of oceanic plates, fluids are released and metasomatize the mantle wedge, leading to the generation of water-rich island arc basalts (e.g., Ryukyu)57. Due to the depletion of high field strength elements (e.g., Nb, Ta) in the released fluids, island arc basalts typically exhibit significant negative anomalies in Nb and Ta57, However, this is inconsistent with the trace element features of the Akesu diabase (Supplementary Fig. S8), which indicates that the water in the mantle source of Akesu diabase does not originate from the metasomatized mantle wedge.

As the oceanic plates continue to subduct, released fluids can further metasomatize the continental lithospheric mantle, the released fluids can further metasomatize the lithospheric mantle, forming a hydrated metasomatized lithospheric mantle similar to the Mesozoic North China Craton (NCC)40,58. According to previous studies, the Tarim CFB can be divided into Group-1 and Group-2 based on the Sr–Nd isotope characteristics59. Compared to Group-2, Group-1 displays an enriched Sr–Nd isotope composition, which is supposed to originate from the lithospheric mantle. As shown in Supplementary Fig. S13, The Sr–Nd isotope composition of the Akesu diabase falls within the range of Group-1. Additionally, using machine learning, Zhang et al.60 suggested that the mantle source of the Tarim LIP is hydrous, and the fluid activity in the Tarim LIP decreases from northeast to southwest, which may be correlated with fluids released from earlier southward subduction of oceanic plate. This propose that the Akesu diabase may originate from the metasomatized lithospheric mantle.

The analogy can be drawn to the Mesozoic NCC, which serves as a prime example of a lithospheric mantle that has undergone metasomatism and is enriched in water40,58. The NCC experienced the continuous subduction of the Pacific Plate during the Mesozoic. As a result, fluids derived from the subducting Pacific Plate infiltrated the lithospheric mantle, leading to the development of a mantle source with significantly high water content. Notable examples of water-enriched basalts in the eastern region of the NCC include Feixian (> 1000 ppm in mantle source)40, Yixian (> 1500 ppm in mantle source)58, and Sihetun (> 1200 ppm in mantle source)58. These Mesozoic basalts consistently exhibit trace element characteristics of Nb–Ta depletion, as illustrated in Supplementary Fig. S8.

While the Akesu diabase also displays weak negative anomalies in Nb–Ta, they are distinctly different from those observed in the Mesozoic NCC basalts (as shown in Supplementary Fig. S8). This disparity suggests that the water found in the mantle source of the Akesu diabase may not primarily originate from a metasomatized lithospheric mantle in a manner similar to the Mesozoic NCC. Furthermore, a high water content (> 1000 ppm) in the lithospheric mantle can potentially induce destabilization and thinning processes40. However, the Tarim Craton has been recognized as a region characterized by geological stability throughout its history61,62.

It is worth noting that the oceanic plates can further subduct into the deep mantle. Seismological observations have shown that oceanic plates can subduct into the lower mantle63, and experimental petrology has demonstrated the existence of hydrous mineral phases in subducted oceanic plates that can remain stable under the high pressure and temperature conditions of the lower mantle64. This suggests that the water from subducted plates can be transported into the deep mantle. When subducted material enters into the mantle plume, it can give rise to the formation of water-rich basalts, which do not exhibit negative anomalies in Nb and Ta, such as Emeishan LIP11,12,23 and Columbia River Basalts21. Moreover, Zhang et al.60 suggested the mantle plume in the Tarim LIP might be hydrous. Additionally, based on numerical simulations, Liu and Leng31 concluded that the mantle plume of Tarim LIP exhibit characteristics of a water-rich composition. All of these studies indicate that the water in the mantle source of Tarim may originate from a mantle plume and be associated with subducted oceanic plates.

On the other hand, the finding of hydrous ringwoodite in ultra-deep diamonds suggests that the mantle transition zone could serve as a significant reservoir for storing water in the Earth's interior65. As the mantle plume passes through the mantle transition zone, melting can occur and transport water from the mantle transition zone into the mantle plume. This could also be a contributing factor to the high water content observed in the Tarim mantle plume.

Conclusions

Considering the effects of fractional crystallization and degassing, it was estimated that the water content of the primitive Tarim CFB magma is higher than 1.23 ± 0.49 wt.%, and the minimum water content and H2O/Ce ratio in the mantle source is 1230 ± 490 ppm and 500, respectively, indicating the potential presence of a hydrous mantle plume under the Tarim Basin.

Because of the existence of a recycled oceanic plate component in the mantle source of the Tarim CFB, the water in the Tarim mantle plume may have originated from the hydrated subducted oceanic plates or the mantle transition zone.

Methods

Analysis of the water content in clinopyroxene phenocrysts

The analysis of the water content in clinopyroxene phenocrysts was carried out in the FTIR (Fourier-transform infrared spectroscopy) laboratory at the School of Earth Sciences, Zhejiang University. The water content of clinopyroxene crystals was analyzed using a Nicolet iS50 FTIR spectrometer equipped with a Continuμm microscope under unpolarized light, with a liquid nitrogen-cooled MCT-A detector. Dry air was used to sweep the infrared equipment and optical path throughout the process. During measurement, the background and number of scans were set to 128 times, the resolution was 4 cm−1, and the wavelength range was 1000–4500 cm−1. The measurements were focused on the center of crystals avoiding cracks or inclusions. The spot radius varied from 30 to 50 μm depending on the sample. For large crystals, multi-point measurements of the core–edge were performed. Finally, the thickness of each infrared analysis area was measured using a Digimatic Indicator ID-F150.

In-situ analysis of major elements in clinopyroxene phenocrysts

The in-situ major element analysis of clinopyroxene phenocrysts was conducted using a Shimadzu EPMA 1720 electron probe at the Electron Probe Microanalysis Laboratory of Zhejiang University's School of Earth Sciences.

The measurement conditions were as follows: acceleration voltage of 15 kV, current of 20 nA, beam spot diameter of 1 μm, and counting time of 30 s for each element (10 s for characteristic bands and 10 s for front and rear background each). Natural minerals and artificially synthesized oxides were used as standard samples, and the ZAF correction process was used to obtain the major elements content of clinopyroxene phenocrysts. For the obtained data, see Supplementary Table S1.

In-situ analysis of trace elements in clinopyroxene phenocrysts

In-situ trace element analyses of clinopyroxene phenocrysts were conducted in the LA-ICP-MS (laser ablation inductively coupled plasma mass spectrometry) laboratory of the School of Earth Sciences, Zhejiang University. The parameters of Analyte HE laser ablation system were as follows: laser wavelength of 193 nm, maximum energy density of 45 J/cm2, frequency of 6 Hz, pulse energy of 4 J/cm2, and ablation beam spot diameter of 60 μm. During measurement, the blank background acquisition time was 10 s, the sample ablation time was 40 s, and there was a 30 s purge between two ablations to clean the injection system. The total time for single-point analysis was 65 s. High-purity helium gas was used as the carrier gas, and the flow rates of this gas in the sample chamber and channel were set to 0.6 and 0.3 L/min, respectively. The mass spectrometer model was iCAPRQ (Thermofisher), and the parameters were set to: 14 L/min plasma cooling gas flow, 0.8 L/min auxiliary gas flow, 0.9 L/min sample gas flow, 1500 W radio frequency power, and a dwell time of 10 ms for each element.

Artificially synthesized glass NIST SRM 610 was used as the external standard, and the Ca content of clinopyroxene analyzed by electron probe was used as the internal standard for result calibration. The standard samples used during analysis were: NIST SRM 610 (artificially synthesized glass), SRM 612 (artificially synthesized glass), BCR-2G (natural glass) and BHVO-2G (natural glass from Hawaii volcanoes). The analysis accuracy of all trace elements was better than 20%. The results are shown in Supplementary Table S1.

Supplementary Information

Acknowledgements

This study was supported by the National Natural Science Foundation of China (grant nos. 41630205 and 41872055), the Fundamental Research Funds for the Central Universities (grant nos. B240201021) and the National Postdoctoral Program for Innovative Talents (grant no. BX201700206).

Author contributions

Q.-K.X and H.C. conceived this study. Y.B. conducted the analyses. H.C. and Y.B. wrote the manuscript with inputs from E.H., T.K., H.-X.W, F.-Q.Z., J.L., X.-Y.G. and Q.-K.X. All authors contributed the interpretation of the results.

Data availability

All the data used in this study are listed in Supplementary Data.

Footnotes

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-60213-4.

References

1. Ernst RE, Buchan KL, Campbell IH. Frontiers in large igneous province research. Lithos. 2005;79:271–297. 10.1016/j.lithos.2004.09.004. [CrossRef] [Google Scholar]
2. Ernst, R. E. Large Igneous Provinces. 1–653 (Cambridge University Press, 2014).
3. Campbell, I. H. Identification of ancient mantle plumes in Mantle Plumes: Their identification Through Time (eds. Ernst, R. E. & Buchan, K. L.) 5–21 (Geological Society of America, 2001).
4. Herzberg C, Gazel E. Petrological evidence for secular cooling in mantle plumes. Nature. 2009;458:619–622. 10.1038/nature07857. [Abstract] [CrossRef] [Google Scholar]
5. Simmons NA, Forte AM, Boschi L, Grand SP. GyPSuM: A joint tomographic model of mantle density and seismic wave speeds. J. Geophys. Res. Solid Earth. 2010;115:B12310. 10.1029/2010JB007631. [CrossRef] [Google Scholar]
6. Asimow PD, Langmuir CH. The importance of water to oceanic mantle melting regimes. Nature. 2003;421:815–820. 10.1038/nature01429. [Abstract] [CrossRef] [Google Scholar]
7. Hirschmann MM, Tenner T, Aubaud C, Withers AC. Dehydration melting of nominally anhydrous mantle: The primacy of partitioning. Phys. Earth Planet. Interiors. 2009;176:54–68. 10.1016/j.pepi.2009.04.001. [CrossRef] [Google Scholar]
8. Green DH, Hibberson WO, Kovács I, Rosenthal A. Water and its influence on the lithosphere–asthenosphere boundary. Nature. 2010;467:448–451. 10.1038/nature09369. [Abstract] [CrossRef] [Google Scholar]
9. Hirth G, Kohlstedt DL. Water in the oceanic upper mantle: implications for rheology, melt extraction and the evolution of the lithosphere. Earth Planet. Sci. Lett. 1996;144:93–108. 10.1016/0012-821X(96)00154-9. [CrossRef] [Google Scholar]
10. Dixon JE, Dixon TH, Bell D, Malservisi R. Lateral variation in upper mantle viscosity: Role of water. Earth Planet. Sci. Lett. 2004;222:451–467. 10.1016/j.epsl.2004.03.022. [CrossRef] [Google Scholar]
11. Liu J, Xia QK, Kuritani T, Hanski E, Yu HR. Mantle hydration and the role of water in the generation of large igneous provinces. Nat. Commun. 2017;8:1824. 10.1038/s41467-017-01940-3. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
12. Liu, J. et al. Compositional variation of Picrites in the Emeishan large igneous province modulated by water in the mantle plume. J. Geophys. Res. Solid Earth127, e2021JB023584 10.1029/2021JB023584 (2022).
13. Sobolev AV, Krivolutskaya NA, Kuzmin DV. Petrology of the parental melts and mantle sources of Siberian trap magmatism. Petrology. 2009;17:253–286. 10.1134/S0869591109030047. [CrossRef] [Google Scholar]
14. Ivanov AV, et al. Volatile concentrations in olivine-hosted melt inclusions from meimechite and melanephelinite lavas of the Siberian traps large igneous province: Evidence for flux-related high-Ti, high-Mg magmatism. Chem. Geol. 2018;483:442–462. 10.1016/j.chemgeo.2018.03.011. [CrossRef] [Google Scholar]
15. Heinonen JS, Luttinen AV. Mineral chemical evidence for extremely magnesian subalkaline melts from the Antarctic extension of the Karoo large igneous province. Mineral. Petrol. 2010;99:201–217. 10.1007/s00710-010-0115-9. [CrossRef] [Google Scholar]
16. Shimizu K, et al. CO2-rich komatiitic melt inclusions in Cr-spinels within beach sand from Gorgona Island, Colombia. Earth Planet. Sci. Lett. 2009;288:33–43. 10.1016/j.epsl.2009.09.005. [CrossRef] [Google Scholar]
17. Gurenko AA, Kamenetsky VS, Kerr AC. Oxygen isotopes and volatile contents of the Gorgona komatiites, Colombia: A confirmation of the deep mantle origin of H2O. Earth Planet. Sci. Lett. 2016;454:154–165. 10.1016/j.epsl.2016.08.035. [CrossRef] [Google Scholar]
18. Trela J, et al. The hottest lavas of the Phanerozoic and the survival of deep Archaean reservoirs. Nat. Geosci. 2017;10:451–456. 10.1038/ngeo2954. [CrossRef] [Google Scholar]
19. Xia QK, et al. High water content in primitive continental flood basalts. Sci. Rep. 2016;6:25416. 10.1038/srep25416. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
20. Stefano CJ, Mukasa SB, Andronikov A, Leeman WP. Water and other volatile systematics of olivine-hosted melt inclusions from the Yellowstone hotspot track. Contrib. Mineral. Petrol. 2011;161:615–633. 10.1007/s00410-010-0553-8. [CrossRef] [Google Scholar]
21. Cabato JA, Stefano CJ, Mukasa SB. Volatile concentrations in olivine-hosted melt inclusions from the Columbia River flood basalts and associated lavas of the Oregon Plateau: Implications for magma genesis. Chem. Geol. 2015;392:59–73. 10.1016/j.chemgeo.2014.11.015. [CrossRef] [Google Scholar]
22. Choudhary BR, Santosh M, De Vivo B, Jadhav G, Babu E. Melt inclusion evidence for mantle heterogeneity and magma degassing in the Deccan large Igneous Province, India. Lithos. 2019;346:105135. 10.1016/j.lithos.2019.07.002. [CrossRef] [Google Scholar]
23. Sun H, Xia Q, Liu J, Wang Z, Bi Y. Water content in Ertan basalts: Implications for the regional distribution of water in the mantle sources of Emeishan large igneous province. Geol. J. China Univ. 2021;27:121–132. 10.16108/j.issn1006-7493.2019108. [CrossRef] [Google Scholar]
24. Workman RK, Hart SR. Major and trace element composition of the depleted MORB mantle (DMM) Earth Planet. Sci. Lett. 2005;231:53–72. 10.1016/j.epsl.2004.12.005. [CrossRef] [Google Scholar]
25. Simons K, Dixon J, Schilling JG, Kingsley R, Poreda R. Volatiles in basaltic glasses from the Easter-Salas y Gomez seamount chain and easter microplate: Implications for geochemical cycling of volatile elements. Geochem. Geophys. Geosyst. 2002;3:1–29. 10.1029/2001GC000173. [CrossRef] [Google Scholar]
26. Bizimis M, Peslier AH. Water in Hawaiian garnet pyroxenites: Implications for water heterogeneity in the mantle. Chem. Geol. 2015;397:61–75. 10.1016/j.chemgeo.2015.01.008. [CrossRef] [Google Scholar]
27. Tenner TJ, Hirschmann MM, Withers AC, Ardia P. H2O storage capacity of olivine and low-Ca pyroxene from 10 to 13 GPa: Consequences for dehydration melting above the transition zone. Contrib. Mineral. Petrol. 2012;163:297–316. 10.1007/s00410-011-0675-7. [CrossRef] [Google Scholar]
28. Yang S, et al. Early Permian Tarim large igneous province in northwest China. Sci. China Earth Sci. 2013;56:2015–2026. 10.1007/s11430-013-4653-y. [CrossRef] [Google Scholar]
29. Xu YG, Wei X, Luo ZY, Liu HQ, Cao J. The early Permian Tarim large igneous province: Main characteristics and a plume incubation model. Lithos. 2014;204:20–35. 10.1016/j.lithos.2014.02.015. [CrossRef] [Google Scholar]
30. Zhang DY, et al. Perovskite and baddeleyite from kimberlitic intrusions in the Tarim large igneous province signal the onset of an end-Carboniferous mantle plume. Earth Planet. Sci. Lett. 2013;361:238–248. 10.1016/j.epsl.2012.10.034. [CrossRef] [Google Scholar]
31. Liu, H. & Leng, W. Tarim large igneous province caused by a wide and wet mantle plume. J. Geophys. Res. Solid Earth125, e2019JB019001 10.1029/2019JB019001 (2020).
32. Zhou MF, et al. OIB-like, heterogeneous mantle sources of Permian basaltic magmatism in the western Tarim Basin, NW China: Implications for a possible Permian large igneous province. Lithos. 2009;113:583–594. 10.1016/j.lithos.2009.06.027. [CrossRef] [Google Scholar]
33. Wei X, Xu YG, Feng YX, Zhao JX. Plume-lithosphere interaction in the generation of the Tarim large igneous province, NW China: Geochronological and geochemical constraints. Am. J. Sci. 2014;314:314–356. 10.2475/01.2014.09. [CrossRef] [Google Scholar]
34. Wu H, et al. Discovery of Permian mafic sills intrusion event in the Sinian system, Northwest Tarim block. Acta Geol. Sin. 2020;94:1869–1882. 10.19762/j.cnki.dizhixuebao.2020079. [CrossRef] [Google Scholar]
35. Cheng XX, et al. The Permian mafic intrusive events in the northwestern margin of the Tarim Basin and their tectonic significance. Acta Petrol. Sin. 2022;38:743–764. 10.18654/1000-0569/2022.03.09. [CrossRef] [Google Scholar]
36. Zhong Y-T, et al. Constraining the duration of the Tarim flood basalts (northwestern China): CA-TIMS zircon U-Pb dating of tuffs. Geol. Soc. Am. Bull. 2021;134:325–334. 10.1130/b36053.1. [CrossRef] [Google Scholar]
37. Hauri EH, Wagner TP, Grove TL. Experimental and natural partitioning of Th, U, Pb and other trace elements between garnet, clinopyroxene and basaltic melts. Chem. Geol. 1994;117:149–166. 10.1016/0009-2541(94)90126-0. [CrossRef] [Google Scholar]
38. Bell DR, Rossman GR. Water in Earth's mantle: The role of nominally anhydrous minerals. Science. 1992;255:1391–1397. 10.1126/science.255.5050.1391. [Abstract] [CrossRef] [Google Scholar]
39. Ingrin J, Skogby H. Hydrogen in nominally anhydrous upper-mantle minerals: Concentration levels and implications. Eur. J. Mineral. 2000;12:543–570. 10.1127/0935-1221/2000/0012-0543. [CrossRef] [Google Scholar]
40. Xia QK, et al. High water content in Mesozoic primitive basalts of the North China Craton and implications on the destruction of cratonic mantle lithosphere. Earth Planet. Sci. Lett. 2013;361:85–97. 10.1016/j.epsl.2012.11.024. [CrossRef] [Google Scholar]
41. O'Leary JA, Gaetani GA, Hauri EH. The effect of tetrahedral Al3+ on the partitioning of water between clinopyroxene and silicate melt. Earth Planet. Sci. Lett. 2010;297:111–120. 10.1016/j.epsl.2010.06.011. [CrossRef] [Google Scholar]
42. Liu J, Xia QK, Deloule E, Chen H, Feng M. Recycled oceanic crust and marine sediment in the source of alkali basalts in Shandong, eastern China: Evidence from magma water content and oxygen isotopes. J. Geophys. Res. Solid Earth. 2015;120:8281–8303. 10.1002/2015JB012476. [CrossRef] [Google Scholar]
43. Chen H, Xia QK, Ingrin J, Deloule E, Bi Y. Heterogeneous source components of intraplate basalts from NE China induced by the ongoing Pacific slab subduction. Earth Planet. Sci. Lett. 2017;459:208–220. 10.1016/j.epsl.2016.11.030. [CrossRef] [Google Scholar]
44. Rudnick, R., Gao, S., Holland, H. & Turekian, K. Composition of the Continental Crust in the Crust (ed. Rudnick, R. L.). 1–64 (Elsevier, 2003).
45. Hauri EH, Gaetani GA, Green TH. Partitioning of water during melting of the Earth's upper mantle at H2O-undersaturated conditions. Earth Planet. Sci. Lett. 2006;248:715–734. 10.1016/j.epsl.2006.06.014. [CrossRef] [Google Scholar]
46. Roeder P, Emslie R. Olivine-liquid equilibrium. Contrib. Mineral. Petrol. 1970;29:275–289. 10.1007/BF00371276. [CrossRef] [Google Scholar]
47. Wade JA, et al. Prediction of magmatic water contents via measurement of H2O in clinopyroxene phenocrysts. Geology. 2008;36:799–802. 10.1130/G24964A.1. [CrossRef] [Google Scholar]
48. Bédard JH. Parameterizations of calcic clinopyroxene—Melt trace element partition coefficients. Geochem. Geophys. Geosyst. 2014;15:303–336. 10.1002/2013GC005112. [CrossRef] [Google Scholar]
49. Patkó L, et al. Extremely low structural hydroxyl contents in upper mantle xenoliths from the Nógrád-Gömör Volcanic Field (northern Pannonian Basin): Geodynamic implications and the role of post-eruptive re-equilibration. Chem. Geol. 2019;507:23–41. 10.1016/j.chemgeo.2018.12.017. [CrossRef] [Google Scholar]
50. Danyushevsky L, et al. The H2O content of basalt glasses from southwest Pacific back-arc basins. Earth Planet. Sci. Lett. 1993;117:347–362. 10.1016/0012-821X(93)90089-R. [CrossRef] [Google Scholar]
51. Dixon JE, Leist L, Langmuir C, Schilling JG. Recycled dehydrated lithosphere observed in plume-influenced mid-ocean-ridge basalt. Nature. 2002;420:385–389. 10.1038/nature01215. [Abstract] [CrossRef] [Google Scholar]
52. Wei X, Xu YG, Luo ZY, Zhao JX, Feng YX. Composition of the Tarim mantle plume: Constraints from clinopyroxene antecrysts in the early Permian Xiaohaizi dykes, NW China. Lithos. 2015;230:69–81. 10.1016/j.lithos.2015.05.010. [CrossRef] [Google Scholar]
53. Cheng Z, et al. Petrogenesis of nephelinites from the Tarim large igneous province, NW China: Implications for mantle source characteristics and plume-lithosphere interaction. Lithos. 2015;220:164–178. 10.1016/j.lithos.2015.02.002. [CrossRef] [Google Scholar]
54. Cheng Z, Zhang Z, Xie Q, Hou T, Ke S. Subducted slab-plume interaction traced by magnesium isotopes in the northern margin of the Tarim large igneous province. Earth Planet. Sci. Lett. 2018;489:100–110. 10.1016/j.epsl.2018.02.039. [CrossRef] [Google Scholar]
55. Wang C, et al. Geochemical and O–C–Sr–Nd isotopic constraints on the petrogenetic link between aillikites and carbonatites in the Tarim large igneous province. J. Petrol. 2021;62:017. 10.1093/petrology/egab017. [CrossRef] [Google Scholar]
56. Sobolev AV, et al. The amount of recycled crust in sources of mantle-derived melts. Science. 2007;316:412–417. 10.1126/science.1138113. [Abstract] [CrossRef] [Google Scholar]
57. Chen, Z. et al. Zinc isotopes of the Mariana and Ryukyu arc‐related lavas reveal recycling of forearc serpentinites into the subarc mantle. J. Geophys. Res. Solid Earth126, e2021JB022261 10.1029/2021JB022261 (2021).
58. Geng X, et al. Thermal-chemical conditions of the North China Mesozoic lithospheric mantle and implication for the lithospheric thinning of cratons. Earth Planet. Sci. Lett. 2019;516:1–11. 10.1016/j.epsl.2019.03.012. [CrossRef] [Google Scholar]
59. Yu X, et al. Petrogenetic model of the Permian Tarim large igneous province. Sci. China Earth Sci. 2017;60:1805–1816. 10.1007/s11430-016-9098-7. [CrossRef] [Google Scholar]
60. Zhang, R. et al. Formation of Tarim large igneous province and strengthened lithosphere revealed through machine learning. J. Geophys. Res. Solid Earth128, e2022JB025772 10.1029/2022JB025772 (2023).
61. Cao HZ, He LJ, Zhang LY. Inversion of background thermal history since the formation of the Tarim Craton. Chin. J. Geophys. (in Chinese) 2019;62:236–247. 10.6038/cjg2019L0802. [CrossRef] [Google Scholar]
62. Zhang C-L, et al. Tectonic framework and evolution of the Tarim Block in NW China. Gondwana Res. 2013;23:1306–1315. 10.1016/j.gr.2012.05.009. [CrossRef] [Google Scholar]
63. Rodríguez EE, et al. Mantle dynamics of the Andean subduction zone from continent-scale teleseismic S-wave tomography. Geophys. J. Int. 2021;224:1553–1571. 10.1093/gji/ggaa536. [CrossRef] [Google Scholar]
64. Ohtani E. Hydrous minerals and the storage of water in the deep mantle. Chem. Geol. 2015;418:6–15. 10.1016/j.chemgeo.2015.05.005. [CrossRef] [Google Scholar]
65. Pearson D, et al. Hydrous mantle transition zone indicated by ringwoodite included within diamond. Nature. 2014;507:221–224. 10.1038/nature13080. [Abstract] [CrossRef] [Google Scholar]
66. Turner SA. Sedimentary record of late Neoproterozoic rifting in the NW Tarim Basin, China. Precambrian Res. 2010;181:85–96. 10.1016/j.precamres.2010.05.015. [CrossRef] [Google Scholar]
67. McDonough WF, Sun SS. The composition of the Earth. Chem. Geol. 1995;120:223–253. 10.1016/0009-2541(94)00140-4. [CrossRef] [Google Scholar]
68. Rollinson, H. R. Using Geochemical Data: Evaluation, Presentation, Interpretation. 1–384 (Routledge, 1993).
69. Herzberg, C. & Asimow, P. D. Petrology of some oceanic island basalts: PRIMELT2. XLS software for primary magma calculation. Geochem. Geophys. Geosyst.9, Q09001 10.1029/2008GC002057 (2008).
70. Le Roux V, Lee CT, Turner S. Zn/Fe systematics in mafic and ultramafic systems: Implications for detecting major element heterogeneities in the Earth’s mantle. Geochim. Cosmochim. Acta. 2010;74:2779–2796. 10.1016/j.gca.2010.02.004. [CrossRef] [Google Scholar]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

Similar Articles 


To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.


Funding 


Funders who supported this work.

National Natural Science Foundation of China (2)

National Natural Science Foundation of China (National Science Foundation of China) (2)

National Postdoctoral Program for Innovative Talents (1)