# Boninite-like intraplate magmas from Manihiki Plateau require ultra-depleted and enriched source components

Roman Golowin, Maxim Portnyagin, Kaj Hoernle, Folkmar Hauff, Andrey Gurenko, Dieter Garbe-Schönberg, Reinhard Werner, Simon Turner

► To cite this version:

Roman Golowin, Maxim Portnyagin, Kaj Hoernle, Folkmar Hauff, Andrey Gurenko, et al.. Boninite-like intraplate magmas from Manihiki Plateau require ultra-depleted and enriched source components. Nature Communications, 2017, 8, pp.14322. 10.1038/ncomms14322. insu-03619403

HAL Id: insu-03619403

https://insu.hal.science/insu-03619403

Submitted on 25 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

Received 18 Sep 2015 | Accepted 16 Dec 2016 | Published 9 Feb 2017

DOI: 10.1038/ncomms14322

OPEN

# Boninite-like intraplate magmas from Manihiki Plateau require ultra-depleted and enriched source components

Roman Golowin<fnref n="1" />, Maxim Portnyagin<fnref n="1" />,², Kaj Hoernle<fnref n="1" />,³, Folkmar Hauff<fnref n="1" />, Andrey Gurenko<fnref n="4" />, Dieter Garbe-Schönberg³, Reinhard Werner<fnref n="1" /> & Simon Turner⁵

The Ontong Java and Manihiki oceanic plateaus are believed to have formed through high-degree melting of a mantle plume head. Boninite-like, low-Ti basement rocks at Manihiki, however, imply a more complex magma genesis compared with Ontong Java basement lavas that can be generated by ~30% melting of a primitive mantle source. Here we show that the trace element and isotope compositions of low-Ti Manihiki rocks can best be explained by re-melting of an ultra-depleted source (possibly a common mantle component in the Ontong Java and Manihiki plume sources) re-enriched by ≤1% of an ocean-island-basalt-like melt component. Unlike boninites formed via hydrous flux melting of refractory mantle at subduction zones, these boninite-like intraplate rocks formed through adiabatic decompression melting of refractory plume material that has been metasomatized by ocean-island-basalt-like melts. Our results suggest that caution is required before assuming all Archaean boninites were formed in association with subduction processes.

<footnote n="1">¹GEOMAR Helmholtz Centre for Ocean Research Kiel, Wischhofstrasse 1-3, 24148 Kiel, Germany. ²V.I.Vernadsky Institute of Geochemistry and Analytical Chemistry, Kosigin street 19, Moscow 119991, Russia. ³Christian-Albrechts-University of Kiel, Ludewig-Meyn-Strasse 10, 24118 Kiel, Germany.</footnote>

<footnote n="4">⁴Centre de Recherches Pétrographiques et Géochimiques, UMR 7358, Université de Lorraine, 54501 Vandoeuvre-lès-Nancy, France. ⁵Department of Earth and Planetary Sciences, Macquarie University, New South Wales 2109, Australia. Correspondence and requests for materials should be addressed to R.G. (email: rgolowin@geomar.de).</footnote>
Large igneous provinces are formed over a relatively short time period (several million years) through massive volcanism related to pressure-release melting of upwelling mantle plume heads¹⁻³. Some basement lavas from the Ontong Java (OJP), Manihiki and Hikurangi Plateaus, located in the Southwest Pacific Ocean (Fig. 1), have similar chemical compositions and ages (≈ 118–125 Myr), consistent with their derivation from a common source through high degrees (≈ 30%) of mantle melting⁴⁻⁷. It has been proposed that they represent fragments of a single ‘super’ plateau⁸,⁹ and that this magmatic event, the largest in the Phanerozoic, covered ≈ 1% of Earth’s surface⁵.

The primitive Kroenke and more evolved but isotopically identical and dominant Kwaimbaita basement lavas from OJP⁶,¹⁰,¹¹ are thought to be derived from a primitive, ancient, lower-mantle reservoir¹². These lavas have flat, primitive-mantle-normalized incompatible-element patterns. The stratigraphically younger and minor Singgalo lavas¹⁰,¹¹ at OJP have geochemical characteristics similar to Manihiki high-Ti basement lavas⁴, for example, drilled at DSDP Site 317 (ref. 5). Owing to their possibly

younger stratigraphic age, we do not consider these lavas further in this paper.

In contrast to OJP, a major part of the Manihiki Plateau is composed of relatively low-Ti, high-MgO rocks, dated at 118–125 Myr (refs 4,13), and related to the basement (main) phase of plateau formation similar to OJP Kroenke/Kwaimbaita lavas. These low-Ti rocks have variable contents of the most incompatible elements, positive relative Nb-Ta anomalies, and most of them have U-shaped patterns on primitive mantle- and mid-ocean-ridge basalt (MORB-)–normalized incompatible element diagrams⁴,¹³. The boninite-like, U-shaped patterns are unique and enigmatic for plume-derived intraplate magmas, indicating depletion and re-enrichment processes involved in their genesis. In previous studies, the origin of the low-Ti rocks was explained by high-degree melting of a mixed source containing depleted mantle wedge material and small amounts of subducted volcaniclastic sediments¹³ and by extensive melting of a hybrid mantle source comprising FOZO-type peridotite and recycled oceanic crust⁴.

<figure><img src="image_2.png" /></figure>

**Figure 1** | Bathymetric map of the Manihiki Plateau. Red circles show dredge sites of glass samples used for this study from the SO193 and SO225 expeditions. Yellow diamonds show additional dredge sites where low-Ti rock samples have been recovered and reported in ref. 4, and blue boxes show dredge sites where low-Ti rocks were recovered during the SO225 expedition, illustrating the large-scale distribution of low-Ti lavas at Manihiki. The contour lines denote depth intervals of 1,000 m with the thick black contour line showing a depth of 5,000 m below sea level. The map was created based on GEBCO bathymetry data⁶⁵.
Here we present geochemical data from rare, fresh glass samples and additional major element data from bulk rock samples from the Manihiki Plateau obtained by dredging during the R/V Sonne SO193 (ref. 4) and SO225 expeditions (Fig. 1). Glass samples from the SO193 cruise yielded ages of 122.9 ± 1.6 and 124.5 ± 1.5 Myr (ref. 4). The fresh glasses provide a unique opportunity for precise geochemical characterization of the parental Manihiki melts, not possible with phenocryst-bearing whole rock samples variably affected by seawater alteration. We show that the origin of the low-Ti Manihiki rocks is related to re-melting of an isotopically distinctive, highly depleted plume source beneath Manihiki and possibly the OJP, comprising small and variable amounts of a recycled HIMU (high time-integrated $^{238}$U/$^{204}$Pb ratio)-like component.

## Results

**Composition of low-Ti glasses.** Relics of fresh glasses were found in strongly palagonitized hyaloclastites and pillow-lava margins. The glasses contain very rare vesicles and fresh olivine microphenocrysts with inclusions of chromium-rich spinel. The glasses are tholeiitic ($SiO_2$ = 51.1–52.4 wt%; $Na_2O + K_2O$ = 1.5–1.8 wt%), characterized by very low $TiO_2$ (0.27–0.44 wt%) and relatively high MgO (9.0–9.8 wt%) contents (Fig. 2a,b, Supplementary Table 1). Some glass and whole-rock samples have $SiO_2$ contents above 52 and as high as 56 wt% at MgO > 8wt%, which is characteristic for boninites¹⁴. High CaO contents (12.9–13.2 wt%) reflect insignificant clinopyroxene fractionation consistent with its absence as a phenocryst phase. The low-Ti glasses also have low but variable $K_2O$ contents (0.02–0.17 wt%) at relatively constant and high Mg# (molar Mg/(Mg + Fe) = 0.65–0.67). Compared with primitive Kroenke lavas from the OJP basement with similar MgO content (9–10 wt% MgO; ref. 15), the low-Ti Manihiki glasses have higher $SiO_2$ (51.4–52.4 versus ~48.9–50.4 wt%) and lower $TiO_2$ (0.27–0.44 versus ~0.68–0.78 wt%) contents. The low-Ti Manihiki glasses contain 0.13–0.25 wt% $H_2O$, 350–546 p.p.m. S and 126–242 p.p.m. Cl (Supplementary Table 1 and ref. 4), which are similar to those in quenched glasses and olivine-hosted melt inclusions from OJP¹⁵,¹⁶. The most $K_2O$-rich glasses (0.09–0.18 wt% $K_2O$) have Cl/$K_2O$ = 0.10–0.14, $H_2O/K_2O$ = 1.4–1.5 and $H_2O/Ce$ ~200, which are within the range of typical MORB and OIB melts¹⁷⁻¹⁹. The $K_2O$-poor glasses (~0.02 wt% $K_2O$) have high Cl/$K_2O$ = 0.61–0.71, $H_2O/K_2O$ = 6.3–7.0 and $H_2O/Ce$ = 1,191–1,258. The relative enrichment of the glasses in Cl and $H_2O$ may indicate magma assimilation of a sea-water-derived component in the crust¹⁶,¹⁹,²⁰. However, Cl and $H_2O$ in the Manihiki glasses exhibit strong positive correlations with incompatible trace elements, whose concentrations are very low in seawater (for example, Th, Nb, rare earth elements (REE)). Such correlations preclude significant $H_2O$ and Cl degassing from magmas before and during eruption and are unexpected if magmas assimilated variable amounts of seawater-derived components, typically brine²¹. Instead, the relative Cl and $H_2O$ enrichment in low-Ti glasses may reflect their source composition, which might be similar to the incompatible-element depleted sources of komatiites from the Abitibi Belt, Canada and Gorgona Island, eastern Pacific²²,²³.

The low-Ti Manihiki glass samples exhibit variable contents of highly incompatible elements and display U-shaped patterns on incompatible-element diagrams with concentrations normalized to N-MORB (Fig. 2c). U-shaped patterns are also particularly well expressed by the REE normalized to chondrites (Fig. 2d). The U-shaped patterns in Manihiki glasses are similar to those in high-Ca boninites from the Troodos ophiolite on Cyprus and the Tonga Arc, whereas Kroenke-type OJP rocks have relatively flat trace element patterns (Fig. 2c). The low-Ti Manihiki glasses are

more depleted in moderately incompatible elements (Lu to Zr) than primitive Kroenke lavas, whereas more incompatible elements (for example, light rare earth elements (LREE), Nb, Ta, Th, U, Ba, Rb) in low-Ti Manihiki glasses range from being more depleted to more enriched (Fig. 2c). Both, low-Ti Manihiki and Kroenke/Kwaimbaita lavas show strong relative depletion in K relative to LREE and relative enrichment in Nb-Ta in the normalized incompatible-element patterns that distinguishes them from subduction-zone-related boninites from Troodos and Tonga, as does the lack of strong relative Pb enrichment of the low-Ti Manihiki glasses (Fig. 2c).

Isotope compositions of the low-Ti glass samples fall into two clusters forming positive correlations on the $^{87}$Sr/$^{86}$Sr versus $^{143}$Nd/$^{144}$Nd and the Pb isotope diagrams, and negative correlations on $^{206}$Pb/$^{204}$Pb versus $^{143}$Nd/$^{144}$Nd and $^{87}$Sr/$^{86}$Sr isotope diagrams (Fig. 3, Supplementary Fig. 1). The most incompatible-element-depleted low-Ti tholeiites have isotopic compositions similar to enriched mid-ocean-ridge basalts (E-MORB; Fig. 3, Supplementary Fig. 1). Glasses with the strongest enrichment in the most incompatible elements have the most radiogenic Pb (Supplementary Table 1) and the least radiogenic Sr and Nd isotope signatures, and trend towards the compositions of some HIMU-like ocean-island basalts (OIBs).

## Discussion

Manihiki low-Ti tholeiites share a number of characteristics with boninites from subduction-related settings, which are defined as rocks with >52 wt% $SiO_2$, >8 wt% MgO and <0.5 wt% $TiO_2$ (ref. 14). Two glass samples, as well as a significant number of whole rock samples, fulfil the major element criteria for boninites in terms of $SiO_2$, $TiO_2$ and MgO content (Fig. 2a,b, Supplementary Table 1, 2). In addition, both the low-Ti Manihiki rocks and boninites have U-shaped incompatible element patterns with low Sm/Yb but high La/Sm (Fig. 2c,d), and low Ti/V (<12) ratios. Based on high CaO/$Al_2O_3$ ratios (0.83–0.88), the mafic low-Ti Manihiki tholeiites are similar to high-Ca-boninites²⁴, such as, for example, those from the Tonga Arc²⁵ and the Troodos upper pillow lavas on Cyprus²⁶, indicating their origin from a clinopyroxene-bearing, depleted lherzolite or harzburgite source²⁷.

A major process for generating the characteristic compositional features of boninites is melting of highly depleted (through previous melt extraction) peridotitic mantle, triggered through addition of a hydrous slab-derived fluid or melt, enriched in incompatible elements²⁵. The hydrous fluid/melt lowers the solidus temperature of the depleted peridotite, causing it to melt²⁴,²⁸. Several geochemical characteristics of the low-Ti Manihiki magmas, however, are not consistent with a subduction origin. There is no evidence of strong, selective enrichment in highly fluid-mobile elements (for example, Rb, Sr, U, Pb). Further, the strong enrichment in Nb and Ta compared with other incompatible elements is unusual for subduction-related magmas (Fig. 2c). In addition, the Manihiki glasses have low $H_2O$ contents (0.13–0.25 wt%, see Supplementary Table 1 and ref. 4). Finally, palaeo-reconstruction studies indicate that the plateau basement formed in an oceanic intraplate setting, possibly near a mid-ocean-ridge spreading centre, but far away from a subduction zone⁸,²⁹.

In contrast to the enriched fluid or melt component in boninites, which is derived from the subducting slab, we show in the following discussion that the re-enrichment of the ultra-depleted mantle source at Manihiki most likely occurred through addition of an enriched HIMU-like mantle plume component to the source. Unlike hydrous fluids involved in boninite petrogenesis, the small amounts (≤1%, see below) of this HIMU
<figure><img src="image_3.png" /><figcaption>**Figure 2** | Selected major oxide and incompatible element compositions of low-Ti Manihiki rocks and glasses. (a) MgO versus SiO₂ and (b) MgO versus TiO₂ plots for low-Ti Manihiki glasses (red circles) and whole rock samples (white circles: this study, Supplementary Table 2; grey circles: data from refs 4,13). All major element oxide data were normalized on a volatile-free basis to 100%. Compositional fields for Tonga and Troodos boninites²⁶,²⁷, Kroenke/Kwaimbaita OJP rocks and glasses⁶,¹⁵, melt inclusions from Ontong Java¹⁶ and Pacific MORB⁶⁶ are shown for reference. Compositional range of boninites (≥52 wt% SiO₂, >8 wt% MgO and <0.5 wt% TiO₂) is shown by grey field with dashed outer line¹⁴. (c) Incompatible element contents in low-Ti Manihiki glasses (black patterns), high-Ca boninites from Tonga (green patterns) and Troodos (blue patterns) normalized to N-MORB. (d) Rare earth elements in the same samples normalized to chondritic abundances. Most low-Ti glass samples from Manihiki have U-shaped patterns similar to boninites, reflecting in part re-enrichment of the most incompatible elements in a depleted source. Direct comparison of incompatible elements shows that boninite-like low-Ti Manihiki glasses do not show enrichment in fluid-mobile elements or negative Nb-Ta anomalies characteristic of boninites. Representative compositions for Troodos and Tonga boninites⁶⁷,⁶⁸ and OJP Kroenke-type rocks (OJP Kroenke; yellow field)⁶ are also shown. OJP inclusions¹⁶ in d have similar labelling as in a,b. N-MORB and chondrite compositions for normalization are from ref. 32.</figcaption></figure>

mantle component probably did not have a major affect on the degree of mantle melting, and instead were mainly responsible for the variability in the concentrations of highly incompatible elements and Sr-Nd-Pb isotope ratios in low-Ti Manihiki magmas.

Variations of radiogenic isotope ratios (Fig. 3; Supplementary Fig. 1; Supplementary Table 1) in the low-Ti Manihiki glasses indicate that they do not simply represent different melt fractions from a single homogeneous mantle source. Following the criteria of ref. 30, linear correlations between incompatible element and isotope ratios with the same denominator (Supplementary Fig. 2) allow us to narrow down the range of possible scenarios for the origin of Manihiki magmas, indicating either mixing of two magmas from different sources ('magma mixing') or fixed degree

of melting of a heterogeneous mantle source, comprising different proportions of depleted and enriched components ('source mixing'), for example, metasomatism of a depleted source with variable amounts of an enriched component. Although these two petrogenetic scenarios cannot be distinguished solely on the basis of incompatible-element (Supplementary Fig. 3) and isotope systematics³⁰, they imply different mass fractions of the two components and consequently different effects on the major element composition of the magmas.

Before we can determine the compositions of the depleted and enriched end members, we need to calculate the compositions of primary Manihiki magmas from the composition of the low-Ti glasses. We calculated the primary magma compositions using the PRIMELT2 model³¹. The glasses are saturated with olivine
<figure><img src="image_4.png" /><figcaption>Figure 3 | Radiogenic isotopic composition of low-Ti Manihiki rocks.</figcaption></figure>

Initial radiogenic isotope ratio plots (with correction for radiogenic ingrowth to 120 Myr ago) of (a) $^{206}_{238}\text{Pb}/^{204}_{232}\text{Pb}$ versus $^{208}_{204}\text{Pb}/^{204}\text{Pb}$ and (b) $^{206}_{204}\text{Pb}/^{204}\text{Pb}$ versus $^{143}_{144}\text{Nd}/^{144}\text{Nd}$ for low-Ti glasses. The stars represent the estimated UDM (red) and enriched (blue) end members for the Manihiki low-Ti tholeiites. The enriched melt end member for the Manihiki low-Ti tholeiites has characteristics similar to HIMU-like OIB lavas from Rurutu and Tubuai⁶⁹ of the Cook-Austral Chain (see discussion, Supplementary Fig. 2 for calculation of the Manihiki-enriched end member). The black lines show mixing of UDM with variable amounts (indicated at tick marks in wt%) of the enriched OIB-like component represented by a mafic Rurutu lava sample from ref. 34. Note: $^{143}_{144}\text{Nd}/^{144}\text{Nd}$ ratio for intermediate sample SO225DR12-2 was not included for end member calculations (see 'Sample material' in the 'Methods' section). Isotope data for the Kroenke/Kwaimbaita field were taken from refs 10,11. The MORB-field data⁷⁰ were projected to 120 Myr using the following parent/daughter ratios from ref. 4: $^{147}_{145}\text{Sm}/^{144}\text{Nd} = 0.25$, $^{87}_{86}\text{Rb}/^{86}\text{Sr} = 0.005$, $^{238}_{204}\text{U}/^{204}\text{Pb} = 10$, $^{235}_{204}\text{U}/^{204}\text{Pb} = 0.073$, $^{232}_{204}\text{Th}/^{204}\text{Pb} = 40$. The Indian MORB⁷⁰ field was projected to 120 Myr using calculated parent/daughter ratios based on the incompatible element composition of E-DMM from ref. 33. The fields for the Tubuai and young Rurutu lavas were also projected to 120 Myr ago and plotted based on compiled literature data and parent/daughter ratios from the calculated OIB source from ref. 36 and references therein.

and spinel, which are the only liquidus phases, and thus the application of the model, involving incremental addition of olivine to the evolved melt, is straightforward. The calculated primary magmas are in equilibrium with olivine Fo$_{91.1-91.3}$, have high MgO (15.7–17.0 wt%) and CaO (10.4–11.0 wt%), and are peridotite-derived compositions (Supplementary Table 3). Concentrations of incompatible elements in the primary magmas

were then calculated by dilution of the trace element concentrations in glasses proportional to the amount of olivine added to the glass compositions.

To estimate the composition of the depleted source component, we first carried out forward modelling of pooled fractional melting, using a primitive mantle composition³² and bulk partition coefficients for peridotite melting³³ to derive melts with similar concentrations of moderately incompatible elements, for example, heavy rare earth elements (HREE), to those in the calculated low-Ti primary melts (Fig. 4, Supplementary Fig. 4). The results of the forward modelling show that up to 15% fractional melting of fertile or moderately depleted mantle peridotite cannot explain the strong depletion of the Manihiki melts in moderately incompatible elements (Fig. 4a). The required depletion can only be achieved with multi-stage melting, when a significantly depleted source (through previous melting) is re-melted, and the magmas produced during the second stage of melting do not mix with magmas of the initial stage of melting, implying a time gap between these two melting events. The best-fit parameters of the two-stage melting scenario obtained by the least-squares fitting of the HREE concentrations in the modelled melts to the Manihiki primary magmas are ~11% melting of fertile peridotite during the first stage and ~9% melting of the residual mantle during the second stage. Modelling with more depleted sources than primitive mantle (for example, Depleted MORB Mantle (DMM) type source³³) results in decreasing the estimate for the total extent of melting. Taking into account the uncertainties in the physical parameters of melting and the initial source composition, the estimates in Fig. 4a only provide a possible (but not unique) origin for the depleted end member of the low-Ti Manihiki primary melts. Although the degrees of melting of the individual melting stages cannot be fully constrained, the conclusion about the necessity of two (or multi)-stage melting to explain the concentrations of moderately incompatible elements in primary Manihiki magmas is robust. In conclusion, the modelling results show that the depleted end member involved in the origin of the low-Ti Manihiki magmas could be highly depleted peridotite (that is, ultra-depleted mantle, UDM; Supplementary table 4) residual after ~11% fractional melting and melt removal.

Re-melting of residual mantle peridotite (and subsequent re-enrichment in incompatible elements) points to a similar type of origin for the low-Ti intraplate rocks and subduction-related boninites²⁴. Since the low-Ti rocks do not show any evidence for a subduction-related origin as discussed above, elevated temperature and further decompression melting related to the emplacement of a hot mantle plume head are likely to have driven the melting under the Manihiki Plateau rather than a hydrous fluid/melt. The overall lower SiO₂ content of the low-Ti tholeiites compared with boninites can also be explained through dry refractory lherzolite melting, since hydrous melting of refractory lherzolite would increase the contribution of orthopyroxene to the melt resulting in a higher SiO₂ melt content²⁸.

A recent study described ultra-depleted melt inclusions in olivines from a Kroenke-type lava (Fig. 2a,b,d), which were proposed to originate from a previously unrecognized ultra-depleted component in the OJP mantle source¹⁶. Since OJP, Hikurangi and Manihiki Plateaus are likely to have formed as a single plateau originally⁸,⁹, they presumably shared some similar source material. The recognition of UDM in the source of the low-Ti Manihiki magmas suggests that this mantle component may be involved in the origin of both OJP and Manihiki magmas.

The isotopic composition of the UDM end member can place important constraints on its origin. Since the modelled UDM incompatible element pattern shows systematically decreasing
<figure><img src="image_5.png" /><figcaption>**Figure 4** | Two-stage mantle melting and source mixing to generate low-Ti primary Manihiki magmas. (a) Primitive-mantle-normalized multi-element plot of incompatible elements showing the compositions of primary low-Ti melts calculated from the compositions of Manihiki glasses (Supplementary Table 3; red circles); the ranges of melt compositions produced by pooled (aggregate) fractional melting from 1% to 15% of primitive mantle (PM) and of residual primitive mantle after 5% (PM-5%) and 10% (PM-10%) fractional melting; and residual mantle compositions after 5% and 10% previous melt extraction (yellow and green triangles, respectively). (b) The compositions of second stage melts produced by 8.7% pooled fractional melting of an UDM source. The UDM source was previously generated by 10.6% melt extraction from PM and subsequent enrichment/metasomatism with variable amounts (indicated in wt%) of melt corresponding in composition to the mafic Rurutu sample #RRT-B-30 from ref. 34. Details of the melt modelling are given in the discussion and input parameters in Supplementary Table 4.</figcaption></figure>

abundances of elements as they become more incompatible with the most incompatible elements (Rb through Ta) showing near zero abundances (<0.001; Fig. 4a, Supplementary Table 4), the most-depleted samples cannot serve as the depleted end member, due to their U-shaped patterns. Following the approach of ref. 30, the isotope composition of the UDM was estimated from linear correlations of the Manihiki low-Ti lavas on plots of x/Pb versus Pb isotope ratios, x/Nd versus $^{143}$Nd/$^{144}$Nd and x/Sr versus $^{87}$Sr/$^{86}$Sr (where x is a highly incompatible element, for example, Ba, which has a negligible—near zero—concentration in the residual mantle after fractional melting). Since the Ba concentration in the UDM is near to zero, ratios of Ba/Pb, Ba/Nd and Ba/Sr will also approach zero and thus the UDM will have the following isotopic composition: $^{143}$Nd/$^{144}$Nd = 0.51285 ± 0.00001, $^{206}$Pb/$^{204}$Pb = 18.00 ± 0.16, $^{207}$Pb/$^{204}$Pb = 15.48 ± 0.03, $^{208}$Pb/$^{204}$Pb = 38.07 ± 0.11, $^{87}$Sr/$^{86}$Sr = 0.70373 ± 0.00007 (Supplementary Fig. 2, Supplementary Table 5; see also Fig. 3, Supplementary Fig. 1). Higher $^{208}$Pb/$^{204}$Pb and more radiogenic $^{87}$Sr/$^{86}$Sr compared with Pacific MORB (Fig. 3, Supplementary Fig. 1) are not consistent with an origin of the UDM end member through melting of depleted Pacific MORB mantle in the lithosphere or asthenosphere above the plume at the time of formation of the Ontong Java and Manihiki plateaus. Furthermore, the more radiogenic $^{143}$Nd/$^{144}$Nd and less radiogenic Pb isotope ratios exclude any connection of the UDM to the presumably dominant Kroenke/Kwaimbaita lavas at OJP. The very high $^{147}$Sm/$^{144}$Nd parent/daughter ratio of ~0.57 in the UDM end member (assuming 10% previous melt extraction from its source, see Supplementary Table 4) places strong constraints on the timing between the depletion and re-enrichment (metasomatic) events. If the interval was ≥18 Myr, the εNd(t) for the UDM would have been <0. Therefore, the time interval between the events is likely to have been within a few million years unless the original UDM source was isotopically very enriched having very unradiogenic Nd isotopic composition (reflecting very low time-integrated

$^{147}$Sm/$^{144}$Nd). Even the $^{147}$Sm/$^{144}$Nd of the two depleted samples (derived from the re-enriched UDM source) is high (up to 0.324) and thus the depletion event of the UDM source is unlikely to have taken place more than several hundred million years ago or the original UDM source was very enriched (εNd(t)<0 if more than 400 Myr ago). The elevated $^{208}$Pb/$^{204}$Pb (at its $^{206}$Pb/$^{204}$Pb, plotting above the Pacific MORB field) and relatively radiogenic $^{87}$Sr/$^{86}$Sr isotopic composition of the UDM end member overlaps almost completely with the range of Indian MORB, suggesting that it was ultimately derived from a source with a composition similar to the source of Indian MORB.

The primitive-mantle-normalized incompatible element patterns of the most enriched low-Ti glasses show enrichment in the moderately to highly incompatible elements, relative enrichments in Nb and Ta and relative depletions in K, Pb and Zr-Hf, general characteristics of HIMU-type lavas (Fig. 4). Below we will test the source and melt mixing models to see which can best explain the compositions of the low-Ti melts.

First, we test the source-mixing model (Fig. 4b) using the composition of primitive (MgO > 10 wt%) HIMU-like lavas from the Rurutu Island, Austral archipelago, Central Pacific³⁴, which have isotope compositions similar to the most LREE-enriched low-Ti rocks, to serve as the enriched end member. The model reproduces the general shape of the incompatible-element patterns and abundances of U, Nb, Ta and LREE in the primitive low-Ti Manihiki magmas well (Fig. 4b). Small discrepancies are however evident, primarily for fluid-mobile incompatible elements (Rb, Ba, U and Sr). A slightly better fit for the incompatible elements in the low-Ti Manihiki primary magmas can be achieved if the most enriched Macquarie Island ophiolite glasses³⁵ (Supplementary Fig. 4a), which also have HIMU-like incompatible element compositions but distinct (Pacific-MORB-like) isotopic compositions, are assumed for the enriched end member. The addition of ~0.06 to 1% of the Rurutu or Macquarie enriched end members to the UDM source can reproduce the full range of low-Ti Manihiki incompatible-
element patterns. A similarly good fit to the Manihiki magmas can be obtained by modelling direct mixing of ultra-depleted melt from the UDM source and an enriched Rurutu or Macquarie melt (Supplementary Fig. 4b,c). The melt-mixing model, however, requires an order of magnitude more enriched melt (1–10%) to contribute to the Manihiki primary magmas than the source-mixing model. The melt-mixing model was therefore rejected as less probable due to the absence of correlations between major and trace elements in Manihiki glasses, expected from the relatively large mass fraction of the enriched melt required to fit the trace-element patterns. The origin of the magmas can thus be qualitatively attributed to melting ultra-depleted mantle metasomatized by the enriched melt component, present either cryptically or as a different lithology.

The possible isotopic composition of the enriched component was calculated by solving the general mixing equations for isotope ratios³⁰ to fit the estimated mass fraction of the enriched component in the UDM obtained from trace element modelling and assuming Sr, Nd and Pb concentrations in the enriched component equal to the Rurutu melt used for the mixing model (see Fig. 4b). The estimated initial isotope composition of the enriched component is ¹⁴³Nd/¹⁴⁴Nd = 0.51279 ± 0.00001, ²⁰⁶Pb/²⁰⁴Pb = 19.87 ± 0.16, ²⁰⁷Pb/²⁰⁴Pb = 15.66 ± 0.03, ²⁰⁸Pb/²⁰⁴Pb = 39.17 ± 0.11, ⁸⁷Sr/⁸⁶Sr = 0.70280 ± 0.00007 (see Fig. 3, Supplementary Fig. 1, Supplementary Table 5). The end member ratios are very similar to those in the most-enriched low-Ti glass sample SO225DR12-3 with ²⁰⁶Pb/²⁰⁴Pb(t) = 19.78, ²⁰⁷Pb/²⁰⁴Pb(t) = 15.65, ²⁰⁸Pb/²⁰⁴Pb(t) = 39.12, ⁸⁷Sr/⁸⁶Sr(t) = 0.70282 and ¹⁴³Nd/¹⁴⁴Nd(t) = 0.51279.

The enriched Manihiki end member has Sr-Nd-Pb isotopic ratios similar to lavas from sources with HIMU-type Ocean Island Basalt lavas from the Cook-Austral chain, such as those from Rurutu (young lavas) and Tubuai³⁴,³⁶ at 120 Myr ago (Fig. 3; Supplementary Fig. 1). Following the prevailing point of view, the HIMU-like enriched component may originate from recycled oceanic crust³⁷.

On the basis of our new and previously published data, OJP and Manihiki Plateau, as well as possibly Hikurangi Plateau, were likely formed by melting of a compositionally heterogeneous plume source. The dominant plume component beneath the OJP and possibly the Hikurangi part of the original super-plateau (before breakup into individual plateau fragments) was the source of Kroenke/Kwaimbaita basement lavas⁵,¹⁰,¹¹ (Fig. 5), which possibly represented a nearly primitive mantle source from the lower mantle¹². The UDM component might also be present in the OJP source, but thus far has never been recognized in whole rock samples, but only as melt inclusions in olivine from one OJP Kroenke-type rock¹⁶. The dominant plume component beneath Manihiki appears to be the UDM metasomatized by HIMU-like melts. Together the UDM (mantle from which melt has been previously extracted) and enriched HIMU-like (recycled ocean crust) end members may represent relatively young, subducted Indian-MORB-like oceanic lithosphere, which stalled in the transition zone. Such oceanic lithosphere could have been entrained and recycled by the OJP/Manihiki/Hikurangi plume head as it upwelled from the lower mantle. Incorporation of the dense subducted ocean lithosphere could also have served to reduce the buoyancy of the upwelling plume head, such that it did not cause major uplift of the seafloor³⁸. This could explain why the OJP, Manihiki and Hikurangi Plateaus remained largely submarine throughout their history on the seafloor. A piece of this entrained oceanic crust (eclogite) may still be located beneath OJP as proposed in a recent study³⁹ to explain the fast shear wave velocity of 4.75 km s⁻¹ observed beneath OJP. Regardless of the origin of the UDM and HIMU-like components, the composition of the mantle source below OJP (Kroenke/Kwaimbaita with

<figure><img src="image_6.png" /><figcaption>**Figure 5 |** A schematic model for the formation of the Ontong Java and Manihiki Plateaus from a geochemically heterogeneous mantle-plume head. The presence of ultra-depleted melts at OJP and Manihiki (UDM) can be explained with a geochemically heterogeneous plume head containing different components beneath the two plateau segments. Beneath OJP, the dominant Kroenke/Kwaimbaita source component upwells and generates the dominant Kroenke/Kwaimbaita basement lavas with flat primitive-mantle-normalized incompatible element patterns. The presumed dominant source component beneath Manihiki is the UDM end member metasomatized by ≤1% HIMU-type melts. We propose that both components beneath Manihiki are derived through recycling of Indian-mantle-like oceanic lithosphere. Ultra-depleted melt inclusions in olivines in a Kroenke sample from OJP¹⁶ indicate that smaller proportions of a possibly common UDM source may also be present beneath the OJP.</figcaption></figure>

minor UDM) and Manihiki (mainly UDM with some HIMU-like component) was largely different (Fig. 5). In conclusion, our results indicate that melting of incompatible-element ultra-depleted and subsequently re-enriched mantle is possible on a large-scale in intraplate settings during the formation of large igneous provinces, especially under young, thin oceanic lithosphere.

The presence of boninites in Archaean rock formations has been cited as evidence that subduction, and thus plate tectonics, initiated in the Archaean⁴⁰. Recently, comparison of rock sequences including boninites from a modern arc system (Izu-Bonin-Mariana) and an Archaean greenstone belt in Canada suggested that subduction of oceanic crust has been operating since the Hadean⁴⁰. High degrees of alteration, however, are common in Archaean rocks, affecting in particular major and mobile trace elements. The results of our study show that rocks with affinities to boninites in terms of major and trace elements (Fig. 2) can also originate in intraplate settings. Such melts may have been more abundant in earlier Earth history due to a hotter geotherm or enhanced plume-related magmatism. Thus, it is important to distinguish Archaean boninites that formed in subduction versus intraplate settings. The principal difference between intraplate boninite-like and most subduction-related boninites may be whether the Nb-Ta anomaly is positive or negative and if Nb/U normalized to N-MORB is greater or less than one. For example, in the Abitibi greenstone belt of Canada, some low-Ti rocks⁴¹ have trace element patterns very similar to those of the low-Ti Manihiki tholeiites (Supplementary Fig. 5). Thus, these Archaean sequences could possibly be plume-related unlike greenstones from the Nuvvuagittuq supracrustal belt in northern Quebec⁴² or Isua, Greenland⁴³, which have arc-like boninite signatures with low Nb/La⁴⁰,⁴². As a consequence, subduction-related Archaean boninites may not be as abundant as previously thought, placing some questions as to the exact initiation and on the extent of subduction, and thus of large-scale plate tectonics, in the early Earth. Future appraisal of Archaean mafic rocks may be able to exploit this distinction.
## Methods

**Sample material.** Analysed samples were dredged during the R/V Sonne 225 expedition at different locations along the flanks of the Danger Island and Suvorov Troughs (Fig. 1). Detailed information about sampling locations and sample descriptions can be found in ref. 44. All the glass samples contain fresh and partly altered olivine crystals, up to 0.5 mm in size. One relatively fresh, vesicular, olivine-phyric basalt sample with fresh glass (sample SO225DR24-10) was recovered at the eastern flank of the Suvorov Trough, containing fresh olivine phenocrysts (0.25–1.00 mm size). Concentrations of mobile elements and oxides (e.g. Rb, Ba, U, Pb, K₂O, P₂O₅) correlate with immobile elements such as Zr, Nb, La and Ce and thus do not appear to have been affected by seawater alteration. The glass samples from the SO193 expedition (ref. 4 and references therein) were re-analysed with electron micro-probe (EMP) and laser ablation inductively coupled plasma source spectrometry (LA-ICP-MS) using the same analysing conditions and standards for better data comparison. The sample SO225DR12-2 contained only small amounts of fresh glass material and thus it was not possible to replicate the radiogenic isotope composition. This sample shows slightly higher $^{143}$Nd/$^{144}$Nd and $^{87}$Sr/$^{86}$Sr ratios compared with the more-enriched glass sample compositions, plotting slightly off the estimated correlation line and therefore data for this sample have not been used to estimate the end member compositions.

Whole-rock samples from the SO225 and SO193 cruise are fine-grained and show different degrees of seawater alteration from fresh to strongly altered samples (loss of ignition (LOI): 1–13 wt%, see Supplementary Table 2).

**Electron microprobe analysis.** Major elements, S and Cl in Manihiki glasses were analysed with a JEOL JXA 8200 wavelength dispersive microprobe at GEOMAR (Kiel, Germany) using a defocused to 5 µm electron beam, 15 kV acceleration voltage, 6 nA beam current and other conditions as in ref. 45. All elements except Cl, S and F were calibrated using basaltic reference glasses USNM 111240/52 VG-2 and USNM 113498/1 VG-A99 (ref. 46). S and Cl were calibrated using scapolite USNM R6600-1 (ref. 46) and F on rhyolite glass KN-18 (ref. 47). The data accuracy was monitored using basaltic glasses USNM 111240/52 VG-2 and USNM 113498/1 VG-A99 (Supplementary Table 6a) and rhyolite glass USNM 72854 VG568 (ref. 46). The data reported in Supplementary Table 1 represent average compositions obtained from 20 to 25 analyses performed on 10–12 glass shards for every sample. All glasses were homogeneous within the precision of analyses estimated from replicate measurements of reference glasses.

**X-ray fluorescence.** Whole-rock samples (powders) were analysed for major elements by X-ray fluorescence spectrometry (XRF) on fused pellets using a Magix Pro PW 2540 XRF at the Institute of Mineralogy and Petrography at the University of Hamburg, Germany. International reference samples GB1-1, JB-3, JB-2, JA-3 and JB-3 provided by the Geological Survey of Japan (https://gbank.gsj.jp/geostandards/gsj1mainj.html)⁴⁸ were measured along with the samples (Supplementary Table 6b).

**Laser ablation inductively coupled plasma mass spectrometry.** Trace element concentrations were determined by LA-ICP-MS using an Agilent 7500s quadrupole instrument coupled with a GeoLas Pro 193 nm excimer laser system at the Institute of Geosciences (IfG), Kiel University. Sample mounts were mounted in a two-volume ablation cell⁴⁵, and the ablated sample aerosol was transported with helium as the carrier gas and mixed with argon prior to introduction into the spectrometer. Glass shards were analysed with 60–80 µm crater size, 10 J cm⁻² laser energy and 10 Hz laser frequency. The total analysis time for each single analysis was set to 80 s (20 s background and 60 s counting time during laser ablation). The measured intensities were converted to concentrations in the Glitter software by using $^{43}$Ca as internal standard and Ca concentrations from EMPA. The SRM NIST612 standard⁵⁰ was used for initial calibration. Matrix correction was applied by using data for KL2-G reference glass⁵¹. Reference glasses GOR128-G and BCR-2G were analysed as unknown, and the data are listed in Supplementary Table 6c. Consistency between EMPA and LA-ICP-MS data was checked using Ti contents, which agreed within 10% rel. between the two techniques. Analytical precision for multiple analyses was typically better than 3% rel. for most elements. Further details of the setup and calibration strategy can be found in ref. 52. The LA-ICP-MS data in Supplementary Table 1 are average trace element concentrations obtained from three to five spots for every sample. In addition, whole-rock powder of the sample SO225DR24-10 was analysed for trace elements using an Agilent 7500cs ICP-MS instrument following acid digestion as described in ref. 53. The trace element data generated by laser ablation ICP-MS on glasses and solution ICP-MS on the bulk sample show excellent agreement on p.p.m. to sub-p.p.m. concentrations, except for elements easily affected in bulk rocks by seawater alteration (for example, Cs, Sb, Rb, U). Elements concentrated in olivine and chromium spinel (Ni, Cr) have systematically higher concentrations in the whole-rock solution ICP-MS data.

**Secondary ion mass spectrometry.** Concentrations of H₂O were determined using the CAMECA IMS 1280 HR2 instrument at CRPG (Nancy, France), following the approach taken by ref. 54. The samples (previously analysed for major element contents by EMPA) were carefully re-polished to remove any

residual carbon from the carbon-coating process (finishing with the 0.25 µm grain size Al₂O₃ suspension), ultrasonically cleaned, removed from the epoxy mounts and then remounted by pressing them into indium metal. The samples were sputtered with a 10 kV $^{133}$Cs⁺ primary beam with a beam current of 0.8–1.5 nA focused to 5–10 µm spots and rastered to 20 × 20 µm. A normal-incidence electron flood gun was used to compensate for sample charge. A field aperture of ~1,000 µm was used to eliminate secondary ion signal that may come from spot margins due to residual carbon surviving the cleaning process. A mass resolving power (MRP) of ~5,000, enough to resolve $^{17}$O from $^{16}$OH and $^{29}$SiH from $^{30}$Si peaks, was applied. The $^{12}$C⁻ (6s), $^{16}$OH⁻ (4), $^{27}$Al⁻ (2s) and $^{30}$Si⁻ (2s) ions were counted after 300 s pre-sputtering during 12 cycles. To establish the calibration curves, reference glasses (experimental glasses M9, M15 and M47 from the Institut für Mineralogie, Leibniz Universität Hannover⁵⁵, as well as natural standard reference glasses ETNA II-6, ETNA II-7 and ALV 981 R23 (ref. 22 and references therein)) with CO₂ = 0–365 p.p.m., H₂O = 0.0–1.5 wt.% and SiO₂ = 48–50 wt.% were analysed. Olivine was also analysed together with the glasses to account for H₂O and CO₂ background.

**Thermal ionization mass spectrometry.** Sr-Nd-Pb isotope analyses were carried out at GEOMAR Helmholtz Centre for Ocean Research, Kiel, Germany by thermal ionization mass spectrometry (TIMS) following the methods outlined in ref. 56. In brief, before sample dissolution in hot HF-HNO₃, approximately 100–250 mg hand-picked glass (125–250 µm) and whole rock (250–500 µm) chips were leached in warm 2 N HCl at 70 °C for 1 h and subsequently triple-rinsed in 18 MΩ water to minimize the effects of alteration and sample handling. Ion chromatography followed established standard procedures⁵⁷,⁵⁸. Pb isotope analyses were performed on a Finnigan MAT 262 RPQ²⁺ using the Pb double-spike (DS) technique of ref. 59. DS-corrected NBS981 values obtained over the course of the study are $^{206}$Pb/$^{204}$Pb = 16.9416 ± 0.0027, $^{207}$Pb/$^{204}$Pb = 15.4990 ± 0.0027, $^{208}$Pb/$^{204}$Pb = 36.7247 ± 0.0069, $^{207}$Pb/$^{206}$Pb = 0.914849 ± 0.000043 and $^{208}$Pb/$^{206}$Pb = 2.16772 ± 0.00013 (n = 84; 2σ external reproducibility) and they compare well with published double- and triple-spike data⁶⁰⁻⁶⁴. Nd and Sr isotope ratios were measured on a Thermo Fisher TRITON TIMS. Within-run bias correction used $^{146}$Nd/$^{144}$Nd = 0.7219 and $^{86}$Sr/$^{88}$Sr = 0.1194. NBS987 and La Jolla reference material was measured five to six times each turret to obtain a normalization value relative to our preferred values of $^{87}$Sr/$^{86}$Sr = 0.710250 and $^{143}$Nd/$^{144}$Nd = 0.511850, respectively. The normalization value was then applied to the sample data of each turret. Over the course of the study, the long-term machine drift compensated external 2σ uncertainties (2SD) are ± 0.000011 (n = 25) for NBS987 and ± 0.000007 (n = 33) for La Jolla. Within-run 2 s.e. errors are equal or smaller than 2 s.d. of the reference materials. Total chemistry blanks were 10 to 30 pg for Pb and below 150 pg for Sr and Nd and thus negligible. Replicate analyses by means of a second sample digest were carried out for SO225 ROV-3-6 whole-rock chips. While Sr and Nd isotope ratios are reproduced within 2 s.d. of the standards, only $^{207}$Pb/$^{204}$Pb of the Pb isotope ratios matches this criteria. The limited reproducibility of Pb isotope ratios involving $^{206}$Pb and $^{208}$Pb is ascribed to the combined effects of (i) seawater interaction during low temperature alteration that can lead to a heterogeneous enrichment of uranium causing variable degrees of $^{206}$Pb ingrowth over time and (ii) Pb removal at high temperatures leading to variable U/Pb and Th/Pb and thus $^{206}$Pb and $^{208}$Pb ingrowth. The slight variations caused by alteration, however, do not affect the overall scientific interpretations and conclusions derived from the Pb isotope data.

**Data availability.** The authors declare that all data supporting the findings of this study are available within the paper and its Supplementary Information files.

## References

1. Richards, M. A., Duncan, R. A. & Courtillot, V. E. Flood basalts and hot-spot tracks: plume heads and tails. *Science* **246**, 103–107 (1989).

2. Campbell, I. H. & Griffiths, R. W. Implications of mantle plume structure for the evolution of flood basalts. *Earth Planet. Sci. Lett.* **99**, 79–93 (1990).

3. Griffiths, R. W. & Campbell, I. H. Stirring and structure in mantle starting plumes. *Earth Planet. Sci. Lett.* **99**, 66–78 (1990).

4. Timm, C. et al. Age and geochemistry of the oceanic Manihiki Plateau, SW Pacific: new evidence for a plume origin. *Earth Planet. Sci. Lett.* **304**, 135–146 (2011).

5. Hoernle, K. et al. Age and geochemistry of volcanic rocks from the Hikurangi and Manihiki oceanic Plateaus. *Geochim. Cosmochim. Acta* **74**, 7196–7219 (2010).

6. Fitton, J. G. & Godard, M. in *Origin and Evolution of the Ontong Java Plateau* Vol. 229 (eds Fitton, J. G., Mahoney, J. J., Wallace, P. J. & Saunders, A. D.) 151–178 (Geological Society, Special Publications, 2004).

7. Herzberg, C. in *Origin and Evolution of the Ontong Java Plateau* Vol. 229 (eds Fitton, J. G., Mahoney, J. J., Wallace, P. J. & Saunders, A. D.) 179–183 (Geological Society, Special Publications, 2004).

8. Taylor, B. The single largest oceanic plateau: Ontong Java–Manihiki–Hikurangi. *Earth Planet. Sci. Lett.* **241**, 372–380 (2006).
9. Davy, B., Hoernle, K. & Werner, R. Hikurangi Plateau: crustal structure, rifted formation, and Gondwana subduction history. *Geochem. Geophys. Geosyst.* **9**, Q07004 (2008).

10. Tejada, M. L. G. et al. in *Origin and Evolution of the Ontong Java Plateau* Vol. 229 (eds Fitton, J. G., Mahoney, J. J., Wallace, P. J. & Saunders, A. D.) 133–150 (Geological Society, Special Publications, 2004).

11. Tejada, M. L. G., Mahoney, J. J., Neal, C. R., Duncan, R. A. & Petterson, M. G. Basement geochemistry and geochronology of Central Malaita, Solomon Islands, with implications for the origin and evolution of the Ontong Java Plateau. *J. Petrol.* **43**, 449–484 (2002).

12. Jackson, M. G. & Carlson, R. W. An ancient recipe for flood-basalt genesis. *Nature* **476**, 316–319 (2011).

13. Ingle, S. et al. Depleted mantle wedge and sediment fingerprint in unusual basalts from the Manihiki Plateau, central Pacific Ocean. *Geology* **35**, 595–598 (2007).

14. Le Bas, M. J. IUGS Reclassification of the high-Mg and Picritic Volcanic Rocks. *J. Petrol.* **41**, 1467–1470 (2000).

15. Roberge, J., White, R. V. & Wallace, P. J. in 2004, *Origin and Evolution of the Ontong Java Plateau* Vol. 229 (eds Fitton, J. G., Mahoney, J. J., Wallace, P. J. & Saunders, A. D.) 239–257 (Geological Society, London, Special Publications, 2004).

16. Jackson, M. G. et al. Ultra-depleted melts in olivine-hosted melt inclusions from the Ontong Java Plateau. *Chem. Geol.* **414**, 137 (2015).

17. Dixon, J. E., Leist, L., Langmuir, C. & Schilling, J. G. Recycled dehydrated lithosphere observed in plume-influenced mid-ocean-ridge basalt. *Nature* **420**, 385–389 (2002).

18. Michael, P. Regionally distinctive sources of depleted morb - evidence from trace-elements and H₂O. *Earth Planet. Sci. Lett.* **131**, 301–320 (1995).

19. Jackson, M. G. et al. Deeply dredged submarine HIMU glasses from the Tuvalu Islands, Polynesia: implications for volatile budgets of recycled oceanic crust. *Geochem. Geophys. Geosyst.* **16**, 3210–3234 (2015).

20. Michael, P. J. Implications for magmatic processes at Ontong Java Plateau from volatile and major element contents of Cretaceous basalt glasses. *Geochem. Geophys. Geosyst.* **1**, 1–17 (2000).

21. Kendrick, M. A., Arculus, R., Burnard, P. & Honda, M. Quantifying brine assimilation by submarine magmas: examples from the Galápagos Spreading Centre and Lau Basin. *Geochim. Cosmochim. Acta* **123**, 150–165 (2013).

22. Sobolev, A. V. et al. Komatiites reveal an Archean hydrous deep-mantle reservoir. *Nature* **531**, 628–632 (2016).

23. Kamenetsky, V. S., Gurenko, A. A. & Kerr, A. C. Composition and temperature of komatiite melts from Gorgona Island, Colombia, constrained from olivine-hosted melt inclusions. *Geology* **38**, 1003–1006 (2010).

24. Crawford, A. J., Falloon, T. & Green, D. H. in *Boninites: And related rocks* (ed. Crawford, A. J.) 1–49 (Unwin Hyman, 1989).

25. Sobolev, A. V. & Danyushevsky, L. V. Petrology and geochemistry of Boninites from the North Termination of the Tonga Trench - constraints on the generation conditions of primary high-Ca Boninite magmas. *J. Petrol.* **35**, 1183–1211 (1994).

26. Cameron, W. E. Petrology and origin of primitive lavas from the Troodos Ophiolite, Cyprus. *Contrib. Mineral. Petrol.* **89**, 239–255 (1985).

27. Falloon, T. J., Green, D. H. & McCulloch, M. T. in *Boninites and Related Rocks* (ed. Crawford, A. J.) 357–395 (Unwin Hyman, 1989).

28. Falloon, T. & Danyushevsky, L. V. Melting of refractory mantle at 1.5, 2 and 2.5 GPa under anhydrous and H₂O-undersaturated conditions: implications for the petrogenesis of high-Ca Boninites and the influence of subduction components on mantle melting. *J. Petrol.* **41**, 257–283 (2000).

29. Viso, R., Larson, R. & Pockalny, R. Tectonic evolution of the Pacific–Phoenix– Farallon triple junction in the South Pacific Ocean. *Earth Planet. Sci. Lett.* **233**, 179–194 (2005).

30. Langmuir, C. H., Vocke, R. V., Hanson, G. N. & Hart, S. R. A general mixing equation with applications to Icelandic basalts. *Earth Planet. Sci. Lett.* **37**, 380–392 (1978).

31. Herzberg, C. & Asimow, P. D. Petrology of some oceanic island basalts: PRIMELT2.XLS software for primary magma calculation. *Geochem. Geophys. Geosyst.* **9**, Q09001 (2008).

32. Sun, S. S. & McDonough, W. F. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. *Geol. Soc. Lond.* **42**, 313–345 (1989).

33. Workman, R. K. & Hart, S. R. Major and trace element composition of the depleted MORB mantle (DMM). *Earth Planet. Sci. Lett.* **231**, 53–72 (2005).

34. Willbold, M. & Stracke, A. Trace element composition of mantle end-members: Implications for recycling of oceanic and upper and lower continental crust. *Geochem. Geophys. Geosyst.* **7**, Q04400 (2006).

35. Kamenetsky, V. S. et al. Enriched end-member of primitive MORB melts: petrology and geochemistry of glasses from Macquarie Island (SW Pacific). *J. Petrol.* **41**, 411–430 (2000).

36. Stracke, A., Bizimis, M. & Salters, V. J. M. Recycling oceanic crust: Quantitative constraints. *Geochem. Geophys. Geosyst.* **4**, 1–33 (2003).

37. Hofmann, A. W. & White, W. M. Mantle plumes from ancient oceanic-crust. *Earth Planet. Sci. Lett.* **57**, 421–436 (1982).

38. Dannberg, J. & Sobolev, S. V. Low-buoyancy thermochemical plumes resolve controversy of classical mantle plume concept. *Nat. Commun.* **6**, 6960 (2015).

39. Covellone, B. M., Savage, B. & Shen, Y. Seismic wave speed structure of the Ontong Java Plateau. *Earth Planet. Sci. Lett.* **420**, 140–150 (2015).

40. Turner, S., Rushmer, T., Reagan, M. & Moyen, J. F. Heading down early on? Start of subduction on Earth. *Geology* **42**, 139–142 (2014).

41. Kerrich, R. & Xie, Q. L. Compositional recycling structure of an Archean super-plume: Nb-Th-U-LREE systematics of Archean komatiites and basalts revisited. *Contrib. Mineral. Petrol.* **142**, 476–484 (2002).

42. O'Neil, J., Francis, D. & Carlson, R. W. Implications of the Nuvvuagittuq Greenstone Belt for the formation of Earth's early crust. *J. Petrol.* **52**, 985–1009 (2011).

43. Polat, A., Hofmann, A. W. & Rosing, M. T. Boninite-like volcanic rocks in the 3.7-3.8 Ga Isua greenstone belt, West Greenland: geochemical evidence for intra-oceanic subduction zone processes in the early Earth. *Chem. Geol.* **184**, 231–254 (2002).

44. Werner, R., Nürnberg, D., Hauff, F. & Party, S. S. S. RV SONNE Fahrtbericht/Cruise Report SO225, Manihiki II Leg 2, The Manihiki Plateau - Origin, structure and effects of oceanic plateaus and Pleistocene dynamic of the West Pacific warm water pool ; 19.11.2012 - 06.01.2013, Suva/Fiji - Auckland/New Zealand. Report No (Kiel, Germany, 2013).

45. Portnyagin, M. et al. H₂O-rich melt inclusions in fayalitic olivine from Hekla volcano: implications for phase relationships in silicic systems and driving forces of explosive volcanism on Iceland. *Earth Planet. Sci. Lett.* **357**, 337–346 (2012).

46. Jarosewich, E., Nelen, J. A. & Norberg, J. A. Reference samples for electron microprobe analysis. *Geostand. Newslett.* **4**, 43–47 (1980).

47. Mosbah, M., Metrich, N. & Massiot, P. Pigme fluorine determination using a nuclear microprobe with application to glass inclusions. *Nucl. Instrum. Methods Phys. Res. B Beam Interact. Mater. Atoms* **58**, 227–231 (1991).

48. Imai, N., Terashima, S., Itoh, S. & Ando, A. 1994 Compilation of analytical data for minor and trace elements in seventeen GSJ geochemical reference samples, 'Igneous Rock Series'. *Geostand. Newslett.* **19**, 135–213 (1995).

49. Fricker, M. B. et al. High spatial resolution trace element analysis by LA-ICP-MS using a novel ablation cell for multiple or large samples. *Int. J. Mass Spectrom.* **307**, 39–45 (2011).

50. Jochum, K. P. et al. Determination of reference values for NIST SRM 610-617 glasses following ISO guidelines. *Geostand. Geoanal. Res.* **35**, 397–429 (2011).

51. Jochum, K. P. et al. MPI-DING reference glasses for *in situ* microanalysis: new reference values for element concentrations and isotope ratios. *Geochem. Geophys. Geosyst.* **7**, Q02008 (2006).

52. Garbe-Schönberg, D. & Müller, S. Nano-particulate pressed powder tablets for LA-ICP-MS. *J. Anal. Atom. Spectrom.* **29**, 990 (2014).

53. Garbe-Schönberg, C.-D. Simultaneous determination of thirty-seven trace elements in twenty-eight international rock standards by ICP-MS. *Geostand. Newslett.* **17**, 81–97 (1993).

54. Hauri, E. H. et al. Matrix effects in hydrogen isotope analysis of silicate glasses by SIMS. *Chem. Geol.* **235**, 352–365 (2006).

55. Shishkina, T. A., Botcharnikov, R. E., Holtz, F., Almeev, R. R. & Portnyagin, M. V. Solubility of H₂O- and CO₂-bearing fluids in tholeiitic basalts at pressures up to 500MPa. *Chem. Geol.* **277**, 115–125 (2010).

56. Jacques, G. et al. Across-arc geochemical variations in the Southern Volcanic Zone, Chile (34.5–38.0°S): Constraints on mantle wedge and slab input compositions. *Geochim. Cosmochim. Acta* **123**, 218–243 (2013).

57. Hoernle, K. et al. Arc-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. *Nature* **451**, 1094–1098 (2008).

58. Hoernle, K. & Tilton, G. R. Sr-Nd-Pb isotope data for Fuerteventura (Canary Islands) basal complex and subaerial volcanics: applications to magma genesis and evolution. *Schweiz. Mineral. Petrogr. Mitt.* **71**, 3–18 (1991).

59. Hoernle, K. et al. On- and off-axis chemical heterogeneities along the South Atlantic Mid-Ocean-Ridge (5-11°S): Shallow or deep recycling of ocean crust and/or intraplate volcanism? *Earth Planet. Sci. Lett.* **306**, 86–97 (2011).

60. Baker, J., Peate, D., Waight, T. & Meyzen, C. Pb isotopic analysis of standards and samples using a 207Pb-204Pb double spike and thallium to correct for mass bias with a double-focusing MC-ICP-MS. *Chem. Geol.* **211**, 275–303 (2004).

61. Baker, J. A., Peate, D. W., Waight, T. E. & Thirlwall, M. F. Reply to the: Comment on 'Pb isotopic analysis of standards and samples using a 207Pb-204Pb double spike and thallium to correct for mass bias with a double focusing MC-ICP-MS' by Baker *et al*. *Chem. Geol.* **217**, 175–179 (2005).

62. Galer, S. J. G. & Abouchami, W. Practical application of lead triple spiking for correction of instrumental mass discrimination. *Mineral. Mag.* **62A**, 491–492 (1998).

63. Thirlwall, M. F. Inter-laboratory and other errors in Pb isotope analyses investigated using a ²⁰⁰Pb-²⁰⁴Pb double spike. *Chem. Geol.* **163**, 299–322 (2000).
64. Thirlwall, M. F. Multicollector ICP-MS analysis of Pb isotopes using a ²⁰⁷Pb-²⁰⁴Pb double spike demonstrates up to 400 ppm/amu systematic errors in Tl-normalization. *Chem. Geol.* **184**, 255–279 (2002).

65. Smith, W. H. F. & Sandwell, D. T. Global sea floor topography from satellite altimetry and ship depth soundings. *Science* **277**, 1956–1962 (1997).

66. Jenner, F. E. & O'Neill, H. S. Major and trace analysis of basaltic glasses by laser-ablation ICP-MS. *Geochem. Geophys. Geosyst.* **13**, Q03003 (2012).

67. Falloon, T. J. *et al*. Boninites and Adakites from the Northern Termination of the Tonga Trench: implications for Adakite Petrogenesis. *J. Petrol.* **49**, 697–715 (2007).

68. König, S., Münker, C., Schuth, S. & Garbe-Schönberg, D. Mobility of tungsten in subduction zones. *Earth Planet. Sci. Lett.* **274**, 82–92 (2008).

69. Stracke, A., Hofmann, A. W. & Hart, S. R. Fozo, Himu, and the rest of the mantle zoo. *Geochem. Geophys. Geosyst.* **6**, 1–20 (2005).

70. Meyzen, C. M. *et al*. Isotopic portrayal of the Earth's upper mantle flow field. *Nature* **447**, 1069–1074 (2007).

## Acknowledgements

We thank Captain Mallon, the R/V Sonne crew and the scientific team of the SO225 cruise for excellent cooperation. We are grateful to M. Thöner, U. Westernströer and S. Hauff for support with EMP, LA-ICP-MS and TIMS analyses respectively. Insightful comments of R. Arculus, M. Jackson and M. L. Tejada helped to substantially improve the manuscript. This study was funded by the German Ministry of Education and Research (BMBF, Manihiki II project; # 03G0225A). This is a CRPG contribution N° 2444.

## Author contributions

K.H., F.H. and R.W. conceived and secured funding for the project; R.W. organized and lead the sampling cruise with help from M.P., F.H. and R.G.; F.H. generated the isotope

data; A.G. and R.G. the H₂O analyses; and R.G., M.P. and D.G.-S. the major and trace element data; R.G. interpreted and modelled the data and wrote the paper together with K.H., M.P. and S.T.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Golowin, R. *et al*. Boninite-like intraplate magmas from Manihiki Plateau require ultra-depleted and enriched source components. *Nat. Commun.* **8**, 14322 doi: 10.1038/ncomms14322 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this

article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material.

To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017