Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Nature Communications volume 15, Article number: 4088 (2024 ) Cite this article Arc Smelter
Compositions of island arc and back-arc basin basalts are often used to trace the recycling of subducted materials. However, the contribution of subducted components to the mantle source during initial arc rifting before back-arc basin spreading is not yet well constrained. The northernmost Mariana arc is ideal for studying this because the transition from rifting to back-arc spreading is happening here. Here we report major and trace element and Pb isotopic compositions of olivine-hosted melt inclusions from lavas erupted during initial rifting at 24°N (NSP-24) and compare them with those in active arc front at 21°N and mature back-arc basin at 18°N. NSP-24 high-K melt inclusions have highly radiogenic Pb compositions and are close to those of the HIMU end-member, suggesting the presence of this component in the magma source. The HIMU-like component may be stored in the over-riding plate and released into arc magma with rifting. HIMU-type seamounts may be subducted elsewhere beneath the Mariana arc, but obvious HIMU-type components appear only in the initial stages of arc rifting due to the low melting degree and being consumed during the process of back-arc spreading.
Subduction zones are important places for chemical recycling in the Earth1. Oceanic island arcs, where continental crust is absent, are key locations for identifying various subduction components and recycling processes2. Arc magmas are characterized by the enrichment of large-ion lithophile elements (LILE) and light rare-earth elements (LREE), relative to mid-ocean ridge basalts (MORB), which is mainly attributed to the dehydration and melting of subducted components, including sediments, altered oceanic crust (AOC)3,4,5 and seamounts6,7,8. The Izu-Bonin-Mariana (IBM) arc is a typical oceanic island arc, where sediment, AOC, and seamounts of the Pacific Plate subduct, making it an ideal system to study subducted component influences on magma genesis4,9,10,11,12. Geochemical studies of IBM lavas show that the contribution of subduction components varies along the length of the arc12. In the Izu-Bonin arc, hydrous fluid released from pelagic sediment and AOC dominate contributions to the mantle wedge, while in the northern Mariana arc, there is a greater contribution of sediment melt4,10,12. The mantle source beneath the Mariana back-arc basin gradually changes from subduction-modified to MORB-like mantle with decreasing latitude and basin widening13,14,15,16. In addition, a few arc volcanoes (Iwo Jima and Kita Iwo Jima) in the southern Izu-Bonin arc have lavas with highly radiogenic Pb compositions (206Pb/204Pb > 19.4), reflecting the result of subducted seamounts with HIMU (high-μ or high time-integrated 238U/204Pb) characteristics10,12,17.
The North Northern Seamount Province (N-NSP) of the Mariana arc between 23°N and 24°N is immediately north of where the Mariana Trough terminates (Fig. 1) and is undergoing northward-propagating back-arc rifting, but seafloor spreading to form oceanic crust has not yet begun. Stern et al.18 suggested that magmas that erupted in this location just prior to initial arc rifting have different mantle source compositions than normal arc magmas. High-precision Pb isotopic studies of whole-rocks show that the compositions of the mantle source beneath this area are complex, involving subducted pelagic sediment, AOC and HIMU-like Cretaceous seamounts, and volcaniclastics10,19. Compared to whole-rocks, olivine-hosted melt inclusions are less affected by later magmatic processes such as magma mixing and fractional crystallization, because they are protected by their host minerals, thus retaining more complete compositional information of the primitive magma and the mantle source20,21,22,23. In addition, melt inclusions also record information about magmatic evolution prior to eruption20, which can help us better understand how subduction components migrate through the mantle wedge. Therefore, in tandem with studying magmatic products during arc rifting, identifying the composition of mantle sources before rifting begins allows us to better understand lithospheric processes at the initiation of rifting, and sheds light on the various subduction components19. In this paper, we report major and trace elements and Pb isotopic analyses of olivine-hosted melt inclusions from basalts at three sites in the northern Mariana arc. These include samples collected from the eastern scarp of the remnant arc West Mariana Ridge (WMR) which is the conjugate pair of the northernmost North Mariana Seamount Province around 24°N (NSP-24; Fig. 1), the Mariana arc at 21°N (NSP-21) and the back-arc basin at 18°N (MT-18). We use these data to identify trace element and Pb isotopic compositions of mantle sources, compare and analyze their spatial and temporal variabilities at different tectonic evolutionary stages, and illuminate the causes controlling these variations.
Map of the Mariana subduction zone showing the locations of the studied samples (a). The base map data is from SRTM15 + V2.5.593. Previously studied whole-rock samples of the corresponding locations and the detailed sample information are shown in Supplementary Data 1. Insets show detailed locations of samples from the NSP-24 (b), NSP-21 (c) and MT-18 (d). MOR - Mid-Ocean Ridge; NSP - Northern Seamount Province; CIP - Central Islands Province; SSP - Southern Seamount Province; NMT - Northern Mariana Trough; CMT - Central Mariana Trough; SMT - Southern Mariana Trough.
The IBM convergent margin marks the eastern edge of the Philippine Sea Plate and formed in response to the westward subduction of the Pacific Plate (Fig. 1). The Mariana Arc and Mariana Trough back-arc basin are different but related magma-tectonic systems. The magmatic arc is subdivided into four segments, from north to south: the southern Volcano Arc (VA, 24–28°N), the Northern Seamount Province (NSP, 20.7–24°N), the Central Islands Province (CIP, 16–20.7°N), and the Southern Seamount Province (SSP, south of 16°N)11,24. The Mariana Trough is subdivided into the Northern Mariana Trough (NMT, 21–24°N), the Central Mariana Trough (CMT, 17.6–21°N), and the Southern Mariana Trough (SMT, 12.5–17.5°N)25,26. The Mariana Trough south of 19°15′N is dominated by seafloor spreading which is propagating northward into the NMT, which is characterized by rifting14,27,28. We focus our study on the CMT and NSP (Fig. 1), which have different mantle sources affected by different proportions and types of subduction components. The CMT includes a mature back-arc spreading ridge with spreading rates ranging from 3.0 to 4.4 cm/a29,30; volcanism occurs near the main spreading axis27,31. The CMT has a maximum width of ~250 km at 18° N and back-arc spreading started at about 6 Ma32,33. CMT crustal structure is similar to that of a slow-spreading mid-oceanic ridge; seismic profiles show that this crust is about 6 km thick11. CMT lavas have geochemical and isotopic characteristics of MORB with minor involvement of subduction components14,28,34,35,36,37. The depth to the subducted slab beneath Mariana arc volcanoes along the magmatic front is ~130 km38. The NSP can be further subdivided into north NSP (N-NSP, 23–24°N) and south NSP (S-NSP, 20.7–23°N) reflecting differences in enriched components in the magma source25,26,39,40 (Fig. 1). The N-NSP is where the propagating Mariana Trough rift intersects the Mariana arc, and the WMR (remnant arc) terminates here19,41. The NMT extensional zone terminates at Nikko volcano (23.1°N)14. The northward termination of the Mariana Trough may also have been caused by the ongoing collision of the Ogasawara Plateau with the IBM arc at 26 °N42 (Fig. 1).
N-NSP and southern VA arc lavas are shoshonitic with distinct enrichments in incompatible major and trace elements and radiogenic isotopes; this arc segment is also called the Alkalic Volcano Province (AVP)26. In contrast, S-NSP lavas are not as enriched in incompatible trace elements as those of the N-NSP10,26,40,43. This study focuses on samples from around 18°N in the CMT, 21°N in the S-NSP, and 24°N in the WMR near the N-NSP (hereafter MT-18, NSP-21, and NSP-24, respectively) (Fig. 1). These correspond to magmatism associated with back-arc spreading, active arc, and arc rifting, respectively14. All studied samples are basalts, with olivine, clinopyroxene, and plagioclase as the main phenocrysts (Supplementary Fig. 1 and Supplementary Data 1, 2). Olivine phenocrysts in all samples are euhedral to subhedral and show no obvious compositional zonation (Supplementary Fig. 1). Round or elliptical melt inclusions 20–60 μm across are randomly distributed in some olivine phenocrysts (Supplementary Fig. 1), indicating a primary origin44. Compositions of all studied samples plot within the range of the previous whole-rock data (Supplementary Figs. 2–4), suggesting that they are representative.
The chemical compositions of host olivines and melt inclusions are presented in Supplementary Data 3. Olivine compositions vary significantly in different sampling areas, with Forsterite contents [Fo, atomic 100×Mg/(Mg + Fe)] decreasing from NSP-24 (83.8–90.7) to MT-18 (76.6–87.8) to NSP-21 (70.8–78.5). Almost all host olivines have high CaO contents (average = 0.18 wt.%, n = 300), indicating that they crystallized from magma, rather than being mantle-derived xenocrysts (CaO <0.1 wt.%45).
Melt inclusion compositions may be modified by crystallization and/or re-equilibration after being entrapped by their host minerals, i.e., post-entrapment crystallization (PEC)20,46,47. The measured major element compositions of melt inclusions were corrected for PEC Fe-loss using the method described by ref. 46 and recalculated with PETROLOG 3 software47. Olivine-melt equilibrium was set by applying the model of ref. 48. The initial FeOt content of the trapped melt can be obtained by using the FeOt fractionation trend of whole-rock samples20,46, which were assumed to be 8.0 wt.% in MT-18, 9.0 wt.% in NSP-21 and 9.0 wt.% in NSP-2426,49. These corrections mainly affect MgO and FeO concentrations but have little effect on other elements, thus the overall geochemical trends and composition changes remain unaffected. We focus on the ratios of incompatible trace elements that are not significantly affected by PEC Fe-loss.
Compositions of melt inclusions after correcting for PEC Fe-loss are listed in Supplementary Data 3 and displayed on the total alkali-silica (TAS) diagram. On the TAS diagram, melt inclusions show larger compositional variation than their host whole-rocks (Supplementary Fig. 2). Most MT-18 and NSP-21 melt inclusions have compositions similar to those of their host whole-rocks, whereas some NSP-24 melt inclusion compositions are more primitive (Supplementary Fig. 2). MT-18 melt inclusions are subalkaline, low-K tholeiites. Some NSP-21 melt inclusions have high enough alkali contents to be classified as alkaline series; these are the first alkaline compositions reported from the S-NSP. According to their alkali contents, we identified two groups of both NSP-21 and NSP-24 melt inclusions. NSP-21 and NSP-24 group 1 melt inclusions plot in the medium-K calc-alkaline series field (Supplementary Fig. 2b). NSP-24 group 2 melt inclusions plot in the high-K calc-alkaline series field, whereas NSP-21 group 2 melt inclusions plot in high-K calc-alkaline to shoshonitic fields (Supplementary Fig. 2b). The major oxides versus MgO contents of different group melt inclusions do not covary systematically (Supplementary Fig. 3), indicating that fractional crystallization was not the only cause of magmatic variability but that different melt sources were also important.
In primitive mantle-normalized trace element “spider” diagrams, MT-18 melt inclusions are similar to E-MORB, except for slightly depleted Nb and Ta and enriched Pb contents (Fig. 2a). In contrast, NSP-21 and NSP-24 melt inclusions are markedly enriched in LILEs and REEs and depleted in high-field-strength elements (HFSEs), exhibiting arc-like trace element distribution patterns (Fig. 2c, e). NSP-21 group 1 melt inclusions are chemically similar to their host whole-rocks (Fig. 2 and Supplementary Fig. 5), whereas group 2 melt inclusions are strongly enriched in Li, Rb, Ba, and K but depleted in Cu, Zn, and Pb (Fig. 2 and Supplementary Fig. 4). NSP-24 group 1 melt inclusions also have similar trace element distribution patterns to their host whole-rocks and other AVP basaltic lavas (e.g., North Hiyosh), but incompatible trace element contents (e.g. Rb, K, P) of group 2 melt inclusions are higher (Fig. 2 and Supplementary Fig. 4, 5). Chondrite-normalized REE patterns of MT-18 melt inclusions are similar to those of E-MORB, showing little fractionation of LREEs and HREEs (Fig. 2b, d, f). In contrast, NSP-21 and NSP-24 melt inclusions show similar strong fractionation of LREEs from HREEs (Fig. 2d, f). NSP-24 group 2 melt inclusions have higher LREE contents than NSP-24 group 1 melt inclusions (Fig. 2f).
Trace element patterns normalized to primitive mantle concentrations (a, c, e), and rare earth element patterns normalized to chondritic values (b, d, f) for melt inclusions and their host whole rocks. Normalized-primitive mantle and -chondrite values are from ref. 94. There are two different groups of melt inclusions in NSP-21 and NSP-24, respectively.
The Pb isotopic compositions of olivine-hosted melt inclusions are listed in Supplementary Data 3. Each of the five types of melt inclusions defined above has a narrow range of Pb isotope compositions. The total range that encompasses most of the global variation seen in oceanic basalts and that is far beyond the host whole-rocks (Fig. 3a). There is no significant relationship between Pb isotope variations of different group melt inclusions and the extent of magmatic fractionation as captured by MgO contents or by the Fo content of host olivines (Supplementary Fig. 6).
The insets show the range of all melt inclusions. Pb isotope characteristics of olivine-hosted melt inclusions compared with modeling results (a). Locations of EM I and FOZO are from ref. 95. Fields for Pacific Ocean-type mantle, Indian Ocean-type mantle, EM I OIB, EM II OIB, and HIMU are from the GEOROC database (http://georoc.mpch-mainz.gwdg.de/georoc). Basaltic whole-rocks fields of NSP, CIP, MT (16–22°N), and AVP (except for Iwo Jima and Kita-Iwo Jima) are compiled from literatures4,10,15,19,35,36,37. The fields of Iwo Jima and Kita-Iwo Jima from southern Izu-Bonin arc with HIMU characteristics are from ref. 10,12,17. The field of Magellan seamount with HIMU characteristics is from ref. 71. Error bars with the 207Pb/206Pb and 208Pb/206Pb data are 2 standard errors (2SE). Modeling calculations indicate that mixing DMM with subducted sediments (SED) and AOC, and then mixing with ~1.5% HIMU component could produce the Pb isotope compositions in NSP-24 group 2 melt inclusions. Modeling parameters can be found in the Supplementary Data 6. Primitive mantle-normalized spidergrams of trace element modeling results of NSP-24 calculated composite mantle component (b). The DMM is mixed with a bulk slab component comprising a mixture of 70% subducted sediment and 30% AOC, and then mixed with 1.5% HIMU lava to yield a composite component. This composite source is assumed to be contained in garnet-lherzolite consisting of 60% olivine, 20% orthopyroxene, 10% clinopyroxene, and 10% garnet96. About 1–3% partial melting of this composite component can produce trace element compositions that are similar to those of NSP-24 group 2 melt inclusions. The calculations use partition coefficients for silicate/basaltic melt from the GERM-KdD database (https://kdd.earthref.org/KdD/).
The 207Pb/206Pb and 208Pb/206Pb values of MT-18 melt inclusions range from 0.8344 to 0.8612 and 2.0607 to 2.1145, with mean values of 0.8467 ± 0.0117 (2 SD, n = 34) and 2.0809 ± 0.0268 (2 SD, n = 34), respectively. These values all plot within the Indian Ocean-type mantle field (Fig. 3a). Actually, the radiogenic isotope composition of the mantle wedge is very commonly observed to be similar to Indian Ocean-type mantle for almost all Western Pacific Subduction zones (e.g., Kamchatka, Kuriles, IBM, Ryukyu)50,51,52,53,54. The 207Pb/206Pb and 208Pb/206Pb ratios of NSP-21 melt inclusions range from 0.8223 to 0.8780 and 2.0380 to 2.1503. NSP-21 group 2 melt inclusions have more old uranogenic and more thorogenic Pb (average 207Pb/206Pb = 0.8653 ± 0.0215, 208Pb/206Pb = 2.1243 ± 0.0494, 2 SD, n = 11) than group 1 melt inclusions (average 207Pb/206Pb = 0.8253 ± 0.0041, 208Pb/206Pb = 2.0430 ± 0.0069, 2 SD, n = 31) and tend towards the enriched mantle type I (EM I) field (Fig. 3a). The two groups of NSP-24 melt inclusions also have very different Pb isotopic compositions. Group 1 melt inclusions have a range of 207Pb/206Pb and 208Pb/206Pb ratios from 0.8202 to 0.8300 and 2.0343 to 2.0508 with mean values of 0.8246 ± 0.0054 (2 SD, n = 24) and 2.0433 ± 0.0096 (2 SD, n = 24), respectively. These values fall within the range of previously studied NSP-24 whole-rocks (Fig. 3a). NSP-24 group 1 melt inclusion Pb isotopes are remarkably similar to NSP-21 group 1 melt inclusions. The 207Pb/206Pb and 208Pb/206Pb ratios of NSP-24 group 2 melt inclusions are lower than those of NSP-24 group 1 melt inclusions, which range from 0.7979 to 0.8078 and 2.0062 to 2.0253, with mean values of 0.8028 ± 0.0043 (2 SD, n = 33) and 2.0119 ± 0.0087 (2 SD, n = 33), respectively. Pb isotopic compositions of NSP-24 group 2 melt inclusions fall outside the Pacific- and Indian Ocean-type mantle range and trend towards the HIMU end-member (Fig. 3a).
Trace element variations between the different groups of melt inclusions in NSP-21 and NSP-24 might result from the fractional crystallization of minerals and partial melting of the mantle source. Although the MgO contents of different group melt inclusions vary, they do not differ significantly in the major minerals likely to be involved in fractional crystallization (such as olivine, pyroxene, and plagioclase) (Supplementary Fig. 3). Nor can fractional crystallization explain the large and systematic variations in Pb isotopic compositions seen between group 1, NSP-21 group 2 and NSP-24 group 2 melt inclusions.
Partial melting model calculations attributed geochemical variations along the Mariana magmatic front to varying degrees of melting of a homogeneous mantle43. Geochemical studies of NMT lavas also show that the degree of mantle partial melting gradually increases northward as the rift propogated in that direction14. In contrast, Bloomer et al.26 argued that melts of the Mariana arc magmatic front were produced by similar degrees of mantle partial melting. In this study, we estimate the partial melting degree of different groups of melt inclusions from ratios of incompatible trace elements (Supplementary Fig. 7). Although the results may vary as the simulation parameters change, the overall modeling trend shows that Mariana back-arc basin basalt is mainly derived from partial melting of spinel-bearing mantle peridotite, which is consistent with the results of ref. 43. In contrast, the Mariana arc magma source contains residual garnet, and the role of garnet increases northward (i.e., NSP-24 > NSP-21) (Supplementary Fig. 7). Different groups of NSP-21 melt inclusions show similar degrees of partial melting, whereas NSP-24 high-K melt inclusions appear to represent lower degrees of melting than do medium-K melt inclusions (Supplementary Fig. 7). However, the difference in melting degree is not enough to control observed incompatible element variations (such as K, Li, Rb, and Ba). Accordingly, the distinctive element compositions of these different groups of melt inclusions may reflect sources affected by variable subduction components rather than different degrees of partial melting of a homogeneous mantle source. This conclusion is also consistent with Pb isotopic variations.
The source of most arc and back-arc magmas is the mantle wedge, which is affected by the addition of subduction components and exhibits enrichment of LILEs relative to HFSEs and REEs1. On the basis of trace element data, Peate and Pearce40 proposed that the NSP subduction component is dominated by silicate melts derived from subducted pelagic sediments and AOC. NSP-21 and NSP-24 melt inclusions are enriched in LILEs and depleted in HFSEs, including negative Nb, Ta, and Ti anomalies (Fig. 2). These can be explained by residual or fractionating minerals containing these elements or by a depleted mantle wedge. Some elements, owing to their differential compatibility in magma, can be used to further constrain the types of subduction components. We use incompatible trace element ratios, which change little during fractional crystallization and partial melting or during post-entrapment crystallization, to constrain the contributions of subduction components to the mantle source. For example, Th is less soluble than REEs and HFSEs in slab-derived fluids, whereas Ba is strongly partitioned into fluids55,56. In addition, Th and Ba abundances in subducted sediments are much higher than their abundances in the mantle5,57.
As is illustrated in Fig. 4, NSP-21 and NSP-24 melt inclusions have highly variable Th/Yb and Ba/Yb ratios, suggesting that subducted sediment contributions are more variable than those of subducted AOC (Fig. 4). Although MT-18 melt inclusions have E-MORB-like trace element distribution patterns, they still show some influence of subduction components, such as negative Nb, Ta and positive Pb anomalies (Fig. 2), slightly elevated Th/Yb and Ba/Yb ratios (Fig. 4) and high Ba/Nb ratios (>7; Supplementary Data 3) as discussed by ref. 58. In addition, the Pb isotopic compositions of MT-18 melt inclusions covers almost the full range of 16-22°N Mariana Trough basalt compositions (Fig. 3a). This suggests that a wide range of subduction components and small-scale mantle heterogeneity interact beneath the Mariana Trough, which also explains the large variation of whole-rock geochemical compositions in the study area (Supplementary Figs. 2–4). Nevertheless, subduction components have less influence on the MT-18 magma source compared to the NSP magma source (Figs. 2, 4).
Plot of Th/Yb versus Nb/Yb (a) and Ba/Yb versus Th/Yb (b) for studied melt inclusions and their host whole-rocks. The red arrow represents mantle source enrichment and the blue arrow shows the addition of the subduction component. Subducted sediment is an important factor in changing the trace element ratios of samples. End-member data are from ref. 94 and Supplementary Data 6.
NSP-21 K-rich melt inclusions have especially high alkali contents and can be classified as alkaline series (Supplementary Fig. 2a), similar to the high-K melts of Kasuga 2 and 3 seamounts from the Kasuga Cross-Chain in the northern Mariana back-arc basin2,49, suggesting that a high-K magmatic component is present beneath the volcanic front, as well as behind it in the Kasuga Cross-Chain (Fig. 1). In addition, NSP-21 K-rich melt inclusions have higher 208Pb/206Pb and 207Pb/206Pb ratios and Li, Rb, and Ba contents, but lower Cu, Zn, and Pb contents, than do coexisting NSP-21 medium-K melt inclusions (Figs. 2, 3, and Supplementary Fig. 4). Slab-derived fluids with high Li, Rb, and Ba contents can metasomatize the subarc mantle and contribute these elements to arc magmas. Potassium abundances in some arc lavas are consistently high59,60, suggesting that a large amount of potassium is retained in subducted sediments by mica (biotite and/or phengite) to considerable depth, where it may be released by partial melting61,62. When sulfides (e.g., chalcopyrite, sphalerite, and galena) crystallize and/or separate from magma, Cu, Zn, and Pb partition into them63, resulting in residual magmas that are highly depleted in chalcophile elements64,65. It should be noted that NSP-21 high-K melt inclusions have low MgO content (<3.07-3.58 wt.%) and a slight correlation between Pb isotopes and MgO (Supplementary Figs. 3, 6), indicating that these melts were products of fractional crystallization. Therefore, it is reasonable to speculate that fractional crystallization and/or sulphide immiscibility control the element anomalies coupled with subducted K-rich slab components that also affect the Pb isotopic compositions of the NSP-21 K-rich melt inclusions.
Except for a few elements (e.g., Rb, K, P), the two groups of NSP-24 melt inclusions have similar trace element distribution patterns, while their Pb isotopes differ significantly (Figs. 2, 3a). Medium-K melt inclusions in NSP-24 and NSP-21 have very similar trace element contents and Pb isotopic compositions, which are consistent with NSP whole-rock compositions (Figs. 2, 3a). In contrast, NSP-24 high-K melt inclusions show unusually low 207Pb/206Pb and 208Pb/206Pb ratios (note that the lowest 207Pb/206Pb and 208Pb/206Pb ratios correspond to the most radiogenic Pb), which is not reflected in NSP basaltic whole-rock compositions but is close to andesitic-dacitic lavas from Iwo Jima (24°45′N) and Kita Iwo Jima (25° 30′N) in the nearby VA10,12,17 (Figs. 1, 3a). However, NSP-24 high-K melt inclusions are all basaltic and have different trace elements distribution patterns than Iwo Jima and Kita Iwo Jima lavas (Supplementary Figs. 2, 5), indicating different mantle source components and/or magmatic evolution processes. This suggests that medium-K melt inclusions represent partial melts of ambient NSP mantle, while the high-K melt inclusions come from different mantle source regions. Figure 3a shows that these low 207Pb/206Pb and 208Pb/206Pb ratios cannot be explained by mixing of depleted MORB mantle (DMM), subduction sediment, and AOC, but rather point to a HIMU end-member component. This is a key observation from our melt inclusion data.
As one of the mantle end-members, the HIMU reservoir is defined by the composition of lavas from Mangaia and Tubuaii islands in the southwest Pacific Ocean and St. Helena in the south Atlantic Ocean, which have radiogenic Pb ratios (206Pb/204Pb > 20.5)66,67. To determine how the HIMU component was added to the magma source of NSP-24 high-K melt inclusions, we carried out a simulation calculation using a bulk components mixing model68,69. The results show that the injection of less than 1.5% HIMU component into a composite of 99% depleted mantle and 1% mixture of subducted sediment and AOC (70:30) can explain the Pb isotopic compositions of NSP-24 high-K melt inclusions (Fig. 3a). To further test the possible petrogenetic link between the HIMU mantle component and the NSP magma source, we calculated trace element distribution patterns (Fig. 3b). Using results from the Pb isotope simulation (Fig. 3a), the calculated melt trace element patterns show that about 1–3% partial melting of this composite produces trace element patterns that are strikingly similar to those of NSP-24 melt inclusions (Fig. 3b).
The formation of the HIMU end-member with a high 206Pb/204Pb ratio requires long-term isolation (~ 2 to 3 Ga), and time-integrated enrichment of Th and U relative to Pb70,71. Two regions in the Earth where this can happen are the mantle transition zone (MTZ)72,73 and the base of the lower mantle near the core-mantle boundary74,75. Geophysical evidence does not reveal any deep-rooted mantle plumes beneath the Mariana arc, which means that a direct contribution of HIMU from the deep mantle to the NSP magmatic source is unrealistic. In contrast, there are many ocean island basalt (OIB)-type seamounts with HIMU-like Pb isotopic characteristics outboard of the Mariana Trench (e.g., Wake and Magellan seamount chains; Figs. 1, 3)76,77,78,79 which are subducted and provide more likely sources of the HIMU component (Fig. 3a). Therefore, we propose that the subduction of these HIMU seamounts is responsible for generating NSP enriched magmas (Fig. 5).
The base map data is from SRTM15 + V2.5.593. Subducted seamounts together with associated volcanogenic sediments, AOC, and mantle form a HIMU-like component, which is diapiric and stored in the lithospheric mantle. The HIMU-bearing component is remobilized and undergoes a small degree of decompression melting when rifting happens with the crust thinned and melt in this initial stage of the rift has significant HIMU characteristics (a); Back-arc basin spreading center is located far away from the trench and is mainly fed by a large amount of MORB-like magma, no HIMU component is injected or this component is consumed or swamped with a large partial melting degree during the back-arc spreading stage (b). WMR – West Mariana Ridge; MA – Mariana Arc; PSP – Philippine Sea Plate.
For the IBM system, there are abundant seamounts in the western Pacific outboard of the Mariana arc (Fig. 1), many of which have HIMU characteristics and are subducting with the Pacific Plate76,80,81. In addition, the HIMU-type seamount subduction signal has also been recognized in the central Mariana forearc metabasalts (16–18°N)8,82. Pb isotopic compositions of NSP-21 and NSP-24 medium-K melt inclusions are close to those of the FOZO end-member (Fig. 3a). FOZO (FOcal Zone83) is usually interpreted to have formed in a manner similar to that of HIMU, but its overall formation age is younger than that of HIMU79,84. As shown in Fig. 3a, the addition of a small amount of HIMU components (~ 0.5%) to the mantle source can explain the Pb isotopic variation of NSP-24 and NSP-21 medium-K melt inclusions and most CIP, NSP, and AVP whole-rock compositions. Therefore, we speculate that a slight HIMU-like component exists widely in the Mariana subduction system, as previously suggested10,19,85.
In contrast to arc magmas such as NSP-21 melt, obvious HIMU characteristics only appear in NSP-24 melt inclusions associated with the earliest stage of rifting (Fig. 1). This suggests that the HIMU component is only mobilized into arc magma by rifting. HIMU components from subducted seamounts may migrate and be stored in the lithospheric mantle. Rifting thins, heats and decompresses the lithosphere, remobilizing some components in it, including rifting HIMU-bearing components (Fig. 5a). This interpretation is consistent with the high Nb/Ta ratios in NSP25 and the anomalous low-velocity zone in the upper mantle beneath the region of rift propagation86. As rifting progresses, asthenospheric magma sources dominate over lithospheric sources, more MORB-type components are involved with the HIMU component gradually becoming less important. Once spreading begins, the HIMU component is swamped by upwelling asthenospheric MORB-type magma and the HIMU signature disappears (Fig. 5b). The element and isotope compositions of MT-18 melt inclusions show MORB-like characteristics (Figs. 2, 3). This may be because the back-arc basin spreading center is located far away from the trench, and is mainly fed by a large amount of MORB-like magma. If so, then no lithospheric HIMU component is present, or it has been diluted. Our melt inclusion study demonstrates that early rift magma retains more complete information about subducted components and that a significant HIMU-like component in the Mariana subduction system may be more widespread than previously thought, at least to the north of 24°N, which may be the key to producing arc potassic/shoshonitic magmas.
Selection and preparation of melt inclusions were carried out at the Guangzhou Institute of Geochemistry, Chinese Academy of Sciences. Melt inclusions were re-homogenized and quenched using a high-temperature furnace following procedures in ref. 87. A platinum packet containing olivine crystals was heated at 1250 °C for 10 min to melt the crystal-glass mixture. Then, the packet was quickly raised to the top of the furnace tube to achieve quenching. Melt inclusion-bearing olivines were selected, mounted on epoxy resin disks, and polished until the melt inclusions were exposed for analysis. The major and trace element analysis of host whole-rocks are presented in the Supplementary Information.
Major elemental compositions of melt inclusions and host olivine phenocrysts were analyzed with a JEOL JXA-8230 electron probe microanalyzer (EPMA) at the Key Laboratory of Submarine Geosciences and Prospecting Techniques, Ministry of Education, Ocean University of China, Qingdao. The operating conditions were 15 kV accelerating voltage, 20 nA beam current, and 2 μm beam diameter for melt inclusions and host olivines. Na and K were measured first with 10 s for the peak and 5 s for the background to minimize volatilization loss. Analytical precisions for major and trace elements were better than 1% and 5%, respectively.
After EPMA analysis of melt inclusions, their trace elements were measured using inductively coupled plasma mass spectrometry (ICP-MS) equipped with a 193-nm (ArF) laser ablation system in the Institute of Oceanology, Chinese Academy of Sciences. The carbon coating was completely removed from each epoxy resin disk before laser ablation. Laser conditions were set to 25 μm beam size, ~4 J cm−2 energy density, and 5 Hz repetition rate. Each analysis consisted of a 20 s gas blank with the laser off and a 20 s sample ablation signal with the laser on. Every eight sample analyses were followed by an analysis of the external reference glasses of BCR-2G, BHVO-2G, and BIR-1G to correct for sensitivity drift and mass discrimination. Si determined by EPMA was used as the internal standard element for correction. All analyzed data were processed using the offline software ICPMSDataCal 11.888. Results show that the measured values of reference standards are in good agreement with the recommended values (Supplementary Data 4). The accuracy and precision of most trace element analyses were better than 10% and 20%, respectively (Supplementary Data 4).
After EPMA and LA-ICP-MS measurements, melt inclusions with diameters >40 μm were chosen for in situ Pb isotope analysis using a Neptune Plus multiple collector inductively coupled plasma mass spectrometer (MC-ICP-MS, Thermo Fisher Scientific, Bremen, Germany) coupled to a 193 nm (ArF) COMPexPro 102 laser ablation system in the Guangxi Key Laboratory of Hidden Metallic Ore Deposits Exploration, Guilin University of Technology, Guilin, China. Prior to this procedure, each disk was washed with purified 0.1% nitric acid and ultrasonically cleaned. Detailed procedures for Pb isotope analysis followed those of ref. 89. The laser parameters were set with a repetition rate of 6 Hz, laser beam diameter of 32 μm, and energy of 5 J cm−2. Helium was the carrier gas (650 ml min−1). Before samples were lasered, 30 s were used to detect the gas blank, and then 30 s were used to laser ablate the samples and collect the ablation signals. In order to improve instrumental sensitivity, an X-skimmer cone and a Jet sample cone were used. All isotopic signals were detected with ion counters under static mode. The NKT-1G international basaltic glass standard90 was chosen as the external standard. Basaltic glass standards NIST 614 and BHVO-2G were analyzed before and after each batch of eight melt inclusion analyses to monitor instrumental drift. Here, we present the ratios of 20xPb/206Pb (x = 208, 207) and 20yPb/204Pb (y = 208, 207, 206) of the melt inclusions. Data processing followed ref. 91. Instrumental mass bias was corrected using the standard–sample–standard bracketing method and any outliers (>2 SD) were excluded. The Pb isotope measurement of every unknown sample was carried out before and after the measurement of standard glass for mass bias correction. Analyses of NIST 614 and BHVO-2G during this study yielded averaged 20yPb/204Pb and 20xPb/206Pb consistent with the recommended values. For most analyses, precisions were better than 0.39% for 20xPb/206Pb and 1.22% for 20yPb/204Pb (2RSD, n = 20), whereas accuracies were better than 0.25% for 20xPb/206Pb and 0.28% for 20yPb/204Pb (Supplementary Data 5). Due to the small melt inclusion lasered size (32μm) and low relative abundance of 204Pb, precisions of 20yPb/204Pb are lower than the 20xPb/206Pb (Supplementary Data 5). Therefore, we emphasize ratios of 20xPb/206Pb for melt inclusions in this study.
The authors declare that the data generated or analyzed during this study are available within its supplementary data in a Source Data file deposited at EarthChem Library (https://doi.org/10.60520/IEDA/113107)92. Source data are provided with this paper.
Stern, R. J. Subduction zones. Rev. Geophys. 40, 1012 (2002).
Stern, R. J., Jackson, M. C., Fryer, P. & Ito, E. O. Sr, Nd and Pb isotopic composition of the Kasuga cross-chain in the Mariana Arc: A new perspective on the K-h relationship. Earth Planet. Sci. Lett. 119, 459–475 (1993).
Article ADS CAS Google Scholar
Tatsumi, Y., Hamilton, D. L. & Nesbitt, R. W. Chemical characteristics of fluid phase released from a subducted lithosphere and origin of arc magmas: Evidence from high-pressure experiments and natural rocks. J. Volcanol. Geotherm. Res. 29, 293–309 (1986).
Article ADS CAS Google Scholar
Elliott, T., Plank, T., Zindler, A., White, W. & Bourdon, B. Element transport from slab to volcanic front at the Mariana arc. J. Geophys. Res. Solid Earth 102, 14991–15019 (1997).
Hermann, J. & Rubatto, D. Accessory phase control on the trace element signature of sediment melts in subduction zones. Chem. Geol. 265, 512–526 (2009).
Article ADS CAS Google Scholar
Staudigel, H., Koppers, A., Plank, T. & Hanan, B. Seamounts in the Subduction Factory. Oceanography 23, 176–181 (2010).
Timm, C. et al. Louisville seamount subduction and its implication on mantle flow beneath the central Tonga–Kermadec arc. Nat. Commun. 4, 1720 (2013).
Article ADS PubMed Google Scholar
Deng, J. et al. Geochemistry of subducted metabasites exhumed from the Mariana forearc: Implications for Pacific seamount subduction. Geosci. Front. 12, 101117 (2021).
Hochstaedter, A. et al. Across-arc geochemical trends in the Izu-Bonin arc: Contributions from the subducting slab. Geochem. Geophys. Geosyst. 2, 1019 (2001).
Sun, C. H. & Stern, R. J. Genesis of Mariana shoshonites: Contribution of the subduction component. J. Geophys. Res. Solid Earth 106, 589–608 (2001).
SternR. J., Fouch, M. J. & Klemperer, S. L. An overview of the Izu-Bonin-Mariana subduction factory, in Inside the Subduction Factory, Geophys. Monogr. Ser., 138, edited by J. Eiler, 175–222, AGU, Washington, D. C (2003).
Ishizuka, O. et al. Processes controlling along‐arc isotopic variation of the southern Izu-Bonin arc. Geochem. Geophys. Geosyst., 8, Q06008 (2007).
Stern, R. J. et al. Enriched back-arc basin basalts from the northern Mariana Trough: implications for the magmatic evolution of back-arc basins. Earth Planet. Sci. Lett. 100, 210–225 (1990).
Article ADS CAS Google Scholar
Gribble, R. F., Stern, R. J., Newman, S., Bloomer, S. H. & O’Hearn, T. Chemical and isotopic composition of lavas from the northern Mariana trough: implications for magmagenesis in back-arc basins. J. Petrol. 39, 125–154 (1998).
Article ADS CAS Google Scholar
Chen, Y., Niu, Y., Xue, Q., Gao, Y. & Castillo, P. An iron isotope perspective on back-arc basin development: Messages from Mariana Trough basalts. Earth Planet. Sci. Lett. 572, 117133 (2021).
Lin, P. N., Stern, R. J., Morris, J. & Bloomer, S. H. Nd- and Sr-isotopic compositions of lavas from the northern Mariana and southern Volcano arcs: Implications for the origin of island arc melts. Contrib. Mineral. Petrol. 105, 381–392 (1990).
Article ADS CAS Google Scholar
Tatsumoto, M. Isotopic composition of lead in volcanic rocks from Hawaii, Iwo Jima, and Japan. J. Geophys. Res. 71, 1721–1733 (1966).
Article ADS CAS Google Scholar
Stern, R. J., Smoot, N. C. & Rubin, M. Unzipping of the Volcano Arc, Japan. Tectonophysics 102, 153–174 (1984).
Article ADS CAS Google Scholar
Ishizuka, O. et al. Migrating shoshonitic magmatism tracks Izu-Bonin-Mariana intra-oceanic arc rift propagation. Earth Planet. Sci. Lett. 294, 111–122 (2010).
Article ADS CAS Google Scholar
Danyushevsky, L. V., McNeill, A. W. & Sobolev, A. V. Experimental and petrological studies of melt inclusions in phenocrysts from mantle-derived magmas: an overview of techniques, advantages and complications. Chem. Geol. 183, 5–24 (2002).
Article ADS CAS Google Scholar
Ren, Z. Y., Ingle, S., Takahashi, E., Hirano, N. & Hirata, T. The chemical structure of the Hawaiian mantle plume. Nature 436, 837–840 (2005).
Article ADS CAS PubMed Google Scholar
Sobolev, A. V., Hofmann, A. W. & Nikogosian, I. K. Recycled oceanic crust observed in ‘ghost plagioclase’ within the source of Mauna Loa lavas. Nature 404, 986–990 (2000).
Article ADS CAS PubMed Google Scholar
Zajacz, Z. & Halter, W. LA-ICPMS analyses of silicate meltinclusions in co-precipitated minerals: Quantification, data analysis and mineral/melt partitioning. Geochim. Cosmochim. Acta 71, 1021–1040 (2007).
Article ADS CAS Google Scholar
Stern, R. J., Bloomer, S. H., Lin, P.-N., Ito, E. & Morris, J. Shoshonitic magmas in nascent arcs: New evidence from submarine volcanoes in the northern Marianas. Geology 16, 426–430 (1988).
Article ADS CAS Google Scholar
Pearce, J. A., Stern, R. J., Bloomer, S. H. & Fryer, P. Geochemical mapping of the Mariana arc-basin system: implications for the nature and distribution of subduction components. Geochem. Geophys. Geosyst. 6, Q07006 (2005).
Bloomer, S. H., Stern, R. J., Fisk, E. & Geschwind, C. H. Shoshonitic volcanism in the northern Mariana Arc: 1. petrographic and major and trace element characteristics. J. Geophys. Res. 94, 4469–4496 (1989).
Article ADS CAS Google Scholar
Fryer, P. Geology of the Mariana Trough, in Back-Arc Basins: Tectonics and Magmatism (Springer, New York 1995).
Martinez, F., Fryer, P., Baker, N. A. & Yamazaki, T. Evolution of backarc rifting: Mariana Trough, 20°-24°N. J. Geophys. Res. Solid Earth 100, 3807–3827 (1995).
Stern, R. J., Bloomer, S. H., Martinez, F., Yamazaki, T. & Harrison, T. M. The composition of back-arc basin lower crust and upper mantle in the Mariana Trough: a first report. Isl. Arc. 5, 354–372 (1996).
Takahashi, N., Kodaira, S., Tatsumi, Y., Kaneda, Y. & Suyehiro, K. Structure and growth of the Izu-Bonin-Mariana arc crust: 1. Seismic constraint on crust and mantle structure of the Mariana arc-back-arc system. J. Geophys. Res. Solid Earth 113, B01104 (2008).
Hussong, D. M. & Fryer, P. Back-arc seamounts and the Sea-MARC II seafloor mapping system. Eos. Trans. AGU 64, 627–632 (1983).
Seama, N. & Yamazaki, T. Asymmetric sea-floor spreading in the central Mariana Trough. Eos Trans. AGU 79, F908 (1998). Fall Meet. Suppl.
Yamazaki, T. & Stern, R. J. Topography and magnetic vector anomalies in the Mariana Trough. JAMSTEC J. Deep-Sea Res. 13, 31–45 (1997).
Gribble, R. F., Stern, R. J. & Bloomer, S. H. MORB mantle and subduction components interact to generate basalts in the southern Mariana Trough back-arc basin. Geochim. Cosmochim. Acta 60, 2153–2166 (1996).
Article ADS CAS Google Scholar
Yan, Q. et al. Geochemistry of axial lavas from the mid-and southern Mariana Trough, and implications for back-arc magmatic processes. Mineral. Petrol. 113, 803–820 (2019).
Article ADS CAS Google Scholar
Li, X. et al. Across-arc variations in Mo isotopes and implications for subducted oceanic crust in the source of back-arc basin volcanic rocks. Geology 49, 1165–1170 (2021).
Article ADS CAS Google Scholar
Zhao, G., Luo, W., Lai, Z., Tian, L. & Xu, C. Influence of subduction components on magma composition in back-arc basins: a comparison between the Mariana and Okinawa troughs. Geol. J. 51, 357–367 (2016).
Hayes, G. P. et al. Slab 2, a comprehensive subduction zone geometry model. Science 362, 58–61 (2018).
Article ADS CAS PubMed Google Scholar
Meen, J. K., Stern, R. J. & Bloomer, S. H. Evidence for magma mixing in the Mariana arc system. Isl. Arc. 7, 443–459 (1998).
Peate, D. W. & Pearce, J. A. Causes of spatial compositional variations in Mariana arc lavas: Trace element evidence. Isl. Arc. 7, 479–495 (1998).
Karig, D. E. Remnant arcs. Geol. Soc. Am. Bull. 83, 1057–1068 (1972).
Hsui, A. T. & Youngquist, S. A dynamic model of the curvature of the Mariana Trench. Nature 318, 455–457 (1985).
Lin, P. N., Stern, R. J. & Bloomer, S. H. Shoshonitic volcanism in the Northern Mariana Arc: 2. Large‐ion lithophile and rare earth element abundances: Evidence for the source of incompatible element enrichments in intraoceanic arcs. J. Geophys. Res. Solid Earth 94, 4497–4514 (1989).
Schiano, P. Primitive mantle magmas recorded as silicate melt inclusions in igneous minerals. Earth Sci. Rev. 63, 121–144 (2003).
Article ADS CAS Google Scholar
Gavrilenko, M. et al. A calcium-in-olivine geohygrometer and its application to subduction zone magmatism. J. Petrol. 57, 1811–1832 (2016).
Article ADS CAS Google Scholar
Danyushevsky, L. V., Della-Pasqua, F. N. & Sokolov, S. Re-equilibration of melt inclusions trapped by magnesian olivine phenocrysts from subduction-related magmas: petrological implications. Contrib. Mineral. Petrol. 138, 68–83 (2000).
Article ADS CAS Google Scholar
Danyushevsky, L. V. & Plechov, P. Petrolog 3: integrated software for modeling crystallization processes. Geochem. Geophys. Geosyst. 12, Q07021 (2011).
Ford, C. E., Russel, D. G., Craven, J. A. & Fisk, M. R. Olivine-liquid equilibria: temperature, pressure and composition dependence of the crystal/liquid cation partition coefficients for Mg, Fe2+, Ca and Mn. J. Petrol. 24, 256–265 (1983).
Article ADS CAS Google Scholar
Fryer, P., Gill, J. B. & Jackson, M. C. Volcanologic and tectonic evolution of the Kasuga seamounts, northern mariana trough: Alvin submersible investigations. J. Volcanol. Geotherm. Res. 79, 277–311 (1997).
Article ADS CAS Google Scholar
Pearce, J. A., Kempton, P. D., Nowell, G. M. & Noble, S. R. Hf-Nd element and isotope perspective on the nature and provenance of mantle and subduction components in western pacific arc-basin systems. J. Petrol. 40, 1579–1611 (1999).
McCarthy, A. et al. Volcaniclastic sandstones record the influence of subducted Pacific MORB on magmatism at the early Izu-Bonin arc. Geochim. et. Cosmochim. Acta 296, 170–188 (2021).
Article ADS CAS Google Scholar
Martynov, A. Y., Kimura, J. I., Martynov, Y. A. & Rybin, A. V. Geochemistry of late Cenozoic lavas on Kunashir Island, Kurile arc. Isl. Arc 19, 86–104 (2010).
Shu, Y. et al. Tracing subducted sediment inputs to the Ryukyu arc-Okinawa Trough system: Evidence from thallium isotopes. Geochim. et. Cosmochim. Acta 217, 462–491 (2017).
Article ADS CAS Google Scholar
Shu, Y. et al. Sources of dehydration fluids underneath the Kamchatka arc. Nat. Commun. 13, 4467 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Brenan, J. M., Shaw, H. F., Ryerson, F. J. & Phinney, D. L. Mineral-aqueous fluid partitioning of trace elements at 900 °C and 2.0 gpa: constraints on the trace element chemistry of mantle and deep crustal fluids. Geochim. Cosmochim. Acta 59, 3331–3350 (1995).
Article ADS CAS Google Scholar
Kessel, R., Schmidt, M. W., Ulmer, P. & Pettke, T. Trace element signature of subduction-zone fluids, melts and supercritical liquids at 120-180 km depth. Nature 437, 724–727 (2005).
Article ADS CAS PubMed Google Scholar
Plank, T. Constraints from thorium/lanthanum on sediment recycling at subduction zones and the evolution of the continents. J. Petrol. 46, 921–944 (2005).
Article ADS CAS Google Scholar
Woodhead, J., Stern, R. J., Pearce, J., Hergt, J. & Vervoort, J. Hf-Nd isotope variation in Mariana Trough basalts: The importance of “ambient mantle” in the interpretation of subduction zone magmas. Geology 40, 539–542 (2012).
Article ADS CAS Google Scholar
Dickinson, W. R. Potash-depth (K-h) relations in continental margin and intra-oceanic arcs. Geology 3, 53–56 (1975).
Price, R. C., Stewart, R. B., Woodhead, J. D. & Smith, I. E. M. Petrogenesis of high-K arc magmas: Evidence from Egmont Volcano, North Island, New Zealand. J. Petrol. 40, 167–197 (1999).
Article ADS CAS Google Scholar
Johnson, M. C. & Plank, T. Dehydration and melting experiments constrain the fate of subducted sediments. Geochem. Geophys. Geosyst. 1, 1007 (1999).
Schmidt, M. W., Vielzeuf, D. & Auzanneau, E. Melting and dissolution of subducting crust at high pressures: The key role of white mica. Earth Planet. Sci. Lett. 228, 65–84 (2004).
Article ADS CAS Google Scholar
Pitcairn, I. K. et al. Sources of metals and fluids in orogenic gold deposits: insights from the Otago and Alpine Schists. N.Z. Econ. Geol. 101, 1525–1546 (2006).
Kamenetsky, V. S. et al. Parental basaltic melts and fluids in eastern Manus backarc Basin: implications for hydrothermal mineralisation. Earth. Planet. Sci. Lett. 184, 685–702 (2001).
Article ADS CAS Google Scholar
Lee, C.-T. A. et al. Copper systematics in arc magmas and implications for crust-mantle differentiation. Science 336, 64–68 (2012).
Article ADS CAS PubMed Google Scholar
Zindler, A. & Hart, S. R. Chemical geodynamics. Annu. Rev. Earth Planet. Sci. 14, 493–571 (1986).
Article ADS CAS Google Scholar
Homrighausen, S. et al. Global distribution of the HIMU end member: Formation through Archean plume-lid tectonics. Earth-Sci. Rev. 182, 85–101 (2018).
Article ADS CAS Google Scholar
Marschall, H. R. & Schumacher, J. C. Arc magmas sourced from mélange diapirs in subduction zones. Nat. Geosci. 5, 862–867 (2012).
Article ADS CAS Google Scholar
Nielsen, S. G. & Marschall, H. R. Geochemical evidence for mélange melting in global arcs. Sci. Adv. 3, e1602402 (2017).
Article ADS PubMed PubMed Central Google Scholar
Stracke, A., Hofmann, A. W. & Hart, S. R. FOZO, HIMU, and the rest of the mantle zoo. Geochem. Geophys. Geosyst. 6, Q05007 (2005).
Cabral, R. A. et al. Anomalous sulphur isotopes in plume lavas reveal deep mantle storage of Archaean crust. Nature 496, 490–493 (2013).
Article ADS CAS PubMed Google Scholar
Kuritani, T., Ohtani, E. & Kimura, J.-I. Intensive hydration of the mantle transition zone beneath China caused by ancient slab stagnation. Nat. Geosci. 4, 713–716 (2011).
Article ADS CAS Google Scholar
Mazza, S. E. et al. Sampling the volatile-rich transition zone beneath Bermuda. Nature 569, 398 (2019).
Article ADS CAS PubMed Google Scholar
Weiss, Y., Class, C., Goldstein, S. L. & Hanyu, T. Key new pieces of the HIMU puzzle from olivines and diamond inclusions. Nature 537, 666–670 (2016).
Article ADS CAS PubMed Google Scholar
Zhang, X. Y. et al. Zinc isotopic evidence for recycled carbonate in the deep mantle. Nat. Commun. 13, 6085 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Wei, X. et al. Co-Occurrence of HIMU and EM1 Components in a Single Magellan Seamount: Implications for the Formation of West Pacific Seamount Province. J. Petrol. 63, egac022 (2022).
Koppers, A. A. P., Staudigel, H., Pringle, M. S. & Wijbrans, J. R. Short-lived and discontinuous intraplate volcanism in the South Pacific: hot spots or extensional volcanism? Geochem. Geophys. Geosyst. 4, 1089 (2003).
Chauvel, C., Marini, J. C., Plank, T. & Ludden, J. N. Hf-Nd input flux in the Izu‐Mariana subduction zone and recycling of subducted material in the mantle. Geochem. Geophys. Geosyst. 10, Q01001 (2009).
Wei, X. et al. Mid-Cretaceous Wake seamounts in NW Pacific originate from secondary mantle plumes with Arago hotspot composition. Chem. Geol. 587, 120632 (2022a).
Staudigel, H., Park, K. H., Pringle, M. S., Rubenstone, J. L. & Zindler, A. The longevity of the South Pacific isotopic and thermal anomaly. Earth Planet. Sci. Lett. 102, 24–44 (1991).
Watts, A. B., Koppers, A. A. P. & Robinson, D. P. Seamount subduction and earthquakes. Oceanography 23, 166–173 (2010).
Liu, H. et al. Slab dehydration and potassium-lithium recycling in the forearc indicated by potassium and lithium isotope compositions of exhumed metabasites. Geochim. Cosmochim. Acta 360, 16–35 (2023).
Article ADS CAS Google Scholar
Hart, S. R., Hauri, E. H., Oschmann, L. A. & Whitehead, J. A. Mantle plumes and entrainment: isotopic evidence. Science 256, 517–520 (1992).
Article ADS CAS PubMed Google Scholar
Thirlwall, M. F. Pb isotopic and elemental evidence for OIB derivation from young HIMU mantle. Chem. Geol. 139, 51–74 (1997).
Article ADS CAS Google Scholar
Hickey-Vargas, R. A refractory HIMU component in the sources of island arc magma. Nature 360, 57–59 (1992).
Article ADS CAS Google Scholar
Isse, T. et al. Seismic structure of the upper mantle beneath the Philippine Sea from seafloor and land observation: implications for mantle convection and magma genesis in the Izu-Bonin-Mariana subduction zone. Earth Planet. Sci. Lett. 278, 107–119 (2009).
Article ADS CAS Google Scholar
Ren, Z. Y. et al. Primary magmas and mantle sources of Emeishan basalts constrained from major element, trace element and Pb isotope compositions of olivine hosted melt inclusions. Geochim. Cosmochim. Acta 208, 63–85 (2017).
Article ADS CAS Google Scholar
Liu, Y. S. et al. Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS. Chin. Sci. Bull. 55, 1535–1546 (2010).
Yu, H. et al. Improved in situ analysis of lead isotopes in low-Pb melt inclusions using laser ablation-multi-collector-inductively coupled plasma-mass spectrometry. Rapid. Commun. Mass. Spectrom. 36, e9383 (2022).
Article CAS PubMed Google Scholar
Elburg, M., Vroon, P., Wagt, B. V. D. & Tchalikian, A. Sr and Pb isotopic composition of five USGS glasses (BHVO-2G, BIR-1G, BCR-2G, TB-1G, NKT-1G). Chem. Geol. 223, 196–207 (2005).
Article ADS CAS Google Scholar
Zhang, L. et al. Lead isotope analysis of melt inclusions by LA-MC-ICP-MS. J. Anal. Spectrom. 29, 1393–1405 (2014).
Li, X. 2024. Melt inclusion data from the Mariana subduction zone, Version 1.0. Interdisciplinary Earth Data Alliance (IEDA). https://doi.org/10.60520/IEDA/113107 (2024).
Tozer, B. et al. Global bathymetry and topography at 15 arc sec: SRTM15. Earth Space Sci. 6, 1847–1864 (2019).
Sun, S. S. & McDonough, W. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geol. Soc. Lond. Spec. Publ. 42, 313–345 (1989).
Saal, A. E. et al. Pb isotopic variability in melt inclusions from the EMI - EMII - HIMU mantle end-members and the role of the oceanic lithosphere. Earth Planet. Sci. Lett. 240, 605–620 (2005).
Article ADS CAS Google Scholar
Walter, M. J. Melting of garnet peridotite and the origin of komatiite and depleted lithosphere. J. Petrol. 39, 29–60 (1998).
Article ADS CAS Google Scholar
We thank Le Zhang and Xiaohong Wang for assisting with the melt inclusion preparation and trace element analysis, and thank to Zhigang Zeng for providing part of the samples in our preliminary research, which provided us with inspiration for this study. Funding for this research was provided by the National Natural Science Foundation of China (NSFC) (42006052) to X.H.L., NSFC (42121005, 42488201) and Scientific and Technological Innovation Project of Laoshan Laboratory (LSKJ202204401) to S.Z.L., the Fundamental Research Funds for Central Universities (2023000-842172002、842172003) to X.H.L., and NSFC (42072061) and the Shandong Excellent Young Scientist Grant (ZR2022YQ32) to L.C. This is UTD Geosciences contribution #1714.
Frontiers Science Center for Deep Ocean Multispheres and Earth System, Key Lab of Submarine Geosciences and Prospecting Techniques, MOE and College of Marine Geosciences, Ocean University of China, Qingdao, 266100, China
Xiaohui Li, Sanzhong Li, Zhiqing Lai, Yanhui Suo & Long Chen
Laboratory for Marine Mineral Resources, Qingdao Marine Science and Technology Center, Qingdao, 266237, China
Xiaohui Li, Sanzhong Li, Yanhui Suo & Long Chen
Key Laboratory of Marine Geology and Environment, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, 266071, China
Institute of Geology and Geoinformation, Geological Survey of Japan/AIST, Central 7, 1-1-1, Higashi, Tsukuba, Ibaraki, 305-8567, Japan
Japan Agency for Marine-Earth Science and Technology, 2-15 Natsushima, Yokosuka, Kanagawa, 237-0061, Japan
Department of Sustainable Earth Systems Science, University of Texas at Dallas, Richardson, TX, 75080, USA
UCD School of Earth Sciences, University College Dublin, Belfield, Dublin, 4, Ireland
Guangxi Key Laboratory of Hidden Metallic Ore Deposits Exploration, Guilin University of Technology, Guilin, 541006, China
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
X.H.L. and S.Z.L. conceived and designed the study. X.H.L., O.I., R.J.S., S.Z.L., I.S. and Y.H.S. wrote the manuscript. X.H.L., Z.Q.L. and H.X.Y. analyzed the data. X.H.L., O.I., R.J.S., S.Z.L. and L.C. interpreted the data. All authors contributed to discussions.
The authors declare no competing interests.
Nature Communications thanks James Gill and Sune Nielsen for their contribution to the peer review of this work. A peer review file is available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Li, X., Ishizuka, O., Stern, R.J. et al. A HIMU-like component in Mariana Convergent Margin magma sources during initial arc rifting revealed by melt inclusions. Nat Commun 15, 4088 (2024). https://doi.org/10.1038/s41467-024-48308-y
DOI: https://doi.org/10.1038/s41467-024-48308-y
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Nature Communications (Nat Commun) ISSN 2041-1723 (online)
Metallurgical Engineer Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.