Articles | Volume 15, issue 2
Research article
08 Feb 2024
Research article |  | 08 Feb 2024

Contribution of carbonatite and recycled oceanic crust to petit-spot lavas on the western Pacific Plate

Kazuto Mikuni, Naoto Hirano, Shiki Machida, Hirochika Sumino, Norikatsu Akizawa, Akihiro Tamura, Tomoaki Morishita, and Yasuhiro Kato

Petit-spot volcanoes, occurring due to plate flexure, have been reported globally. As the petit-spot melts ascend from the asthenosphere, they provide crucial information of the lithosphere–asthenosphere boundary. Herein, we examined the lava outcrops of six monogenetic volcanoes formed by petit-spot volcanism in the western Pacific. We then analyzed the 40Ar/39Ar ages, major and trace element compositions, and Sr, Nd, and Pb isotopic ratios of the petit-spot basalts. The 40Ar/39Ar ages of two monogenetic volcanoes were ca. 2.6 Ma (million years ago) and ca. 0 Ma. The isotopic compositions of the western Pacific petit-spot basalts suggest geochemically similar melting sources. They were likely derived from a mixture of high-μ (HIMU) mantle-like and enriched mantle (EM)-1-like components related to carbonatitic/carbonated materials and recycled crustal components. The characteristic trace element composition (i.e., Zr, Hf, and Ti depletions) of the western Pacific petit-spot magmas could be explained by the partial melting of  5 % crust bearing garnet lherzolite, with 10 % carbonatite flux to a given mass of the source, as implied by a mass-balance-based melting model. This result confirms the involvement of carbonatite melt and recycled crust in the source of petit-spot melts. It provides insights into the genesis of tectonic-induced volcanoes, including the Hawaiian North Arch and Samoan petit-spot-like rejuvenated volcanoes that have a similar trace element composition to petit-spot basalts.

1 Introduction

Among the upper-mantle-derived alkali basaltic lavas in oceanic settings, those on thicker plates away from the mid-ocean ridge, could be divided into plume-related and non-plume-related volcanoes. Plume-related North Arch and post-erosional (rejuvenated-stage) volcanoes have been reported in Hawaii and Samoa (Bianco et al., 2005; Bizimis et al., 2013; Clague and Frey, 1982; Clague and Moore, 2002; Dixon et al., 2008; Frey et al., 2000; Garcia et al., 2016; Hart et al., 2004; Konter and Jackson, 2012; Koppers et al., 2008; Reinhard et al., 2019; Yang et al., 2003). Non-plume-related intraoceanic alkali volcanoes, known as petit-spot volcanoes, probably originate where nearby plate subduction causes plate flexures and upwelling of asthenospheric magma (Hirano et al., 2006; Hirano and Machida, 2022; Machida et al., 2015, 2017; Yamamoto et al., 2014, 2018, 2020). The occurrence of petit-spot volcanism supports the presence of melt at the lithosphere–asthenosphere boundary (LAB) (below the area at least).

The occurrence of melt in the uppermost asthenosphere could be attributed to small-scale convection, the presence of hydrous or carbonatitic components, or the uplift of the lithosphere in response to plate flexure; however, the possibility of such an occurrence remains ambiguous (e.g., Bianco et al., 2005; Hua et al., 2023; Korenaga, 2020). The presence of CO2 and carbonated/carbonatitic materials is a significant factor in the formation of alkaline silica-undersaturated melt in the upper mantle (Dasgupta and Hirschmann, 2006; Dasgupta et al., 2007, 2013; Kiseeva et al., 2013; Novella et al., 2014). Experimental studies have shown that the solidus of carbonate-bearing peridotite is lower than that of CO2-free peridotite (Falloon and Green, 1989, 1990; Foley et al., 2009; Ghosh et al., 2009). Moreover, carbonatites and Si-undersaturated melts are generated through the partial melting of CO2-bearing or carbonated peridotite. The produced melts can exhibit continuous chemical variations, depending on pressure (i.e., depth). Carbonatitic melts are produced in the deep asthenosphere (300–110 km), while carbonated or alkali silicate melts are generated in the shallower upper mantle (from  110 to  75 or 60 km) (Keshav and Gudfinnsson, 2013; Massuyeau et al., 2015, 2021). Primary carbonated silicate magma and evolved alkali basalts have been simultaneously observed at the post-spreading ridge in the South China Sea (Zhang et al., 2017; Zhong et al., 2021). The occurrence of Hawaiian rejuvenated volcanoes can be attributed to a carbonatite-metasomatized source with or without silicate metasomatism (Borisova and Tilhac, 2021; Dixon et al., 2008; Zhang et al., 2022).

Submarine petit-spot volcanoes on the subducting northwestern (NW) Pacific Plate may have originated from carbonate-bearing materials and crustal components (pyroxenite and/or eclogite), based on characteristic trace elements, i.e., enriched mantle (EM)-1-like Sr, Nd, and Pb isotopic compositions, and relatively low Mg isotopic compositions (Liu et al., 2020; Machida et al., 2009, 2015). Particularly, the depletion of specific high-field-strength elements (HFSEs) (i.e., Zr, Hf, and Ti) and the abundance of CO2 in petit-spot basalts imply that their melting sources are related to carbonated materials (Hirano and Machida, 2022; Okumura and Hirano, 2013). The nature of the uppermost part of the asthenosphere beneath the oldest Pacific Plate aged 160 Ma was characterized using the eruptive ages and geochemical properties of six newly observed petit-spot volcanoes and lava outcrops. We verified the contribution of carbonatitic components and crustal materials to the melting source of petit-spot volcanoes to understand the nature of the underlying lithosphere–asthenosphere system and model the geodynamic evolution of the region.

Figure 1Geological and geophysical information of the study area. (a) Bathymetry of the western Pacific near the Mariana Trench. The red box shows the study area to the southeast of Minamitorishima (Marcus) Island (Fig. 2). The bathymetric data are adopted from ETOPO1 (NOAA National Geophysical Data Center;, last access: 10 January 2022). (b) Seafloor age map of the same area as panel (a). This study area is on a 160–170 Ma Pacific Plate, called the Jurassic quiet zone (JQZ) (Tivey et al., 2006). The present absolute motion of the Pacific Plate and the seafloor age are derived from studies by Gripp and Gordon (1990) and Müller et al. (2008), respectively. (c) The cross section P-wave tomography beneath the thick yellow line including the study area on the ETOPO1 bathymetry map (left). The bathymetric images were drawn using the Generic Mapping Tool (GMT6; Wessel et al., 2019). The tomographic image (right) was drawn using the SubMachine (Hosseini et al., 2018;, last access: 16 February 2023) on applying the data of Lu et al. (2019).

2 Background

Over the last 20 years, there has been an increase in the understanding of petit-spot volcanic settings, providing valuable insights into the nature of the lithosphere–asthenosphere system, particularly in the NW Pacific region (Hirano et al., 2006; Hirano and Machida, 2022). As with other implications, subducted petit-spot volcanic fields with geological disturbances on the seafloor play a role in controlling the hypocentral regions of megathrust earthquakes (Fujiwara et al., 2007; Fujie et al., 2020; Akizawa et al., 2022). Additionally, the vestige of hydrothermal activity due to petit-spot magmatism has recently been reported (Azami et al., 2023).

Petit-spot melts emerging from the asthenosphere, which are unrelated to the mantle plume, could play a crucial role in clarifying the nature of the LAB (Hirano and Machida, 2022). Their asthenospheric origin was supported by mid-ocean ridge basalt (MORB)-like noble-gas isotopic ratios, multiphase saturation experiment, and geochemistry (Hirano et al., 2006; Hirano and Machida, 2022; Machida et al., 2015, 2017; Yamamoto et al., 2018). The LAB is recognized as a discontinuous transition in seismic velocities at the base of the lithosphere, and its causes are attributed to hydration, melting, and mineral anisotropy, with considerations for the unique characteristics in each tectonic setting (e.g., Rychert and Shearer, 2009). The occurrence of petit-spot volcanoes confirms the existence of melt at the LAB beneath the area at least (Hirano et al., 2006). Recently, similar volcanic activities have been observed globally, including in Java (Sunda) Trench, Tonga Trench, Peru–Chile Trench, Mariana Trench, Costa Rica, North American Basin and Range, and the southern offshore of Greenland, implying the universal occurrence of petit-spot and similar magmatisms (Axen et al., 2018; Buchs et al., 2013; Falloon et al., 2022; Hirano et al., 2013, 2016, 2019; Reinhard et al., 2019; Taneja et al., 2016; Uenzelmann-Neben et al., 2012; Yamamoto et al., 2018, 2020; Zhang et al., 2019). Although the question of whether the LAB discontinuity is due to the differences in the physical properties of minerals (e.g., Hirth and Kohlstedt, 1996; Kang and Karato, 2023; Karato and Jung, 1998; Katsura and Fei, 2021; Stixrude and Lithgow-Bertelloni, 2005; Wang et al., 2006) or the presence of partial melts remains open (e.g., Audhkhasi and Singh, 2022; Chantel et al., 2016; Conrad et al., 2011; Debayle et al., 2020; Herath et al., 2022; Hua et al., 2023; Kawakatsu et al., 2009; Mierdel et al., 2007; Sakamaki et al., 2013; Yoshino et al., 2006), the occurrence of petit-spot volcanism indicates the partial melting of the asthenospheric mantle in the region because they erupted on the seafloor without hotspot and ridge activities (Hirano et al., 2006; Hirano and Machida, 2022; Machida et al., 2015, 2017; Yamamoto et al., 2014, 2018, 2020).

Figure 2Detailed bathymetry of the study area. The on-board multibeam data were surveyed during the YK10-05 and the YK18-08 cruises by the Japan Agency for Marine–Earth Science and Technology (JAMSTEC). The petit-spot knolls and outcrops were investigated during several dives as 6K#1466, 6K#1521, 6K#1522, 6K#1542, and 6K#1544. The pink-colored stars represent the sampling points. The age information was obtained in the present study and Hirano et al. (2019). The bathymetric image was drawn using the Generic Mapping Tools (GMT) software (Wessel et al., 2019).

Table 1Information of the collected western Pacific petit-spot basalts.

a The samples which have no data of palagonite and/or Mn crust thickness are due to the lack of samples or crumbled samples. b This is a reference value due to the lack of radiogenic 40Ar in this sample.

Download Print Version | Download XLSX

The petit-spot volcanic province on the abyssal plain of the western Pacific is surrounded by Cretaceous seamounts and oceanic islands of the western Pacific Seamount Province (Koppers et al., 2003) and is located  100 km southeast of the Minamitorishima (Marcus) Island (Fig. 1a). The study area corresponds to the oldest portion of the Pacific Plate, aged at 160 Ma, and the foot of the outer-rise bulge related to the Mariana subduction system (Hirano et al., 2019; Fig. 1b). Despite several seamounts crosscutting, a subduction-related forebulge in front of the Mariana Trench was detected in satellite gravity maps and has been numerically modeled (Bellas et al., 2022; Hirano et al., 2019; Zhang et al., 2014, 2020). Petrography, geochemistry, and geochronology of petit-spot basalts and zircons in peperites collected from a knoll suggest that petit-spot magmas in this region ascend from the asthenosphere along the concavely flexed plate in response to subduction into the Mariana Trench at younger than ∼3 Ma (Yamamoto et al., 2018; Hirano et al., 2019). Below the study area, a low seismic velocity zone is observed under the lithosphere (Lu et al., 2019; Fig. 1c). Notwithstanding the low-velocity anomalies crosscutting the lower mantle (Fig. 1c), no active hotspots (i.e., heat supplies) have been reported around the western Pacific petit-spot province, which is surrounded by Marcus–Wake Seamount chains, including Minamitorishima Island and Paleogene intraplate volcanoes (Koppers et al., 2003; Aftabuzzaman et al., 2021; Hirano et al., 2021). Other petit-spot lava outcrops were observed in a volcanic cluster during three research cruises, using the research vessel (R/V) Yokosuka (YK16-01, YK18-08, and YK19-05S), with five dives using the submersible, Shinkai 6500 (6K#1466, 6K#1521, 6K#1522, 6K#1542, and 6K#1544; Fig. 2); here, fresh basalts were collected. Information related to the sampling point, depth, and thickness of palagonite rind and manganese crust, as well as the age of the western Pacific petit-spot basalts are provided in Table 1.

3 Field observations, sample locations, and petrography

Here, the eruption sites of monogenetic volcanoes or lava outcrops are approximately aligned with each dive site numbered 6K#1466, #1521, #1522, #1542, and #1544 and conducted using Shinkai 6500. The 6K#1466 dive was conducted at two types of monogenetic volcanoes, categorized as glassy type (R3) and crystalline and vesicular type (R6 and R7), based on the geochemical and petrographic descriptions and occurrence of basaltic samples.

Figure 3Bathymetric map with photos of the outcrop, the collected samples, and their photomicrographs with detailed bathymetry of the sampling points. (a) The 6K#1466, (b) 6K#1521, (c) 6K#1522, and (d) 6K#1542 and 6K#1544 dives, using Shinkai 6500 by JAMSTEC. The 1 m gridded bathymetry of the 6K#1544 dive is shown in panel (d), obtained using an MBES equipped with Shinkai 6500 over a 100 m resolution map obtained using the surface ship, R/V Yokosuka (Kaneko et al., 2022). The photomicrographs of representative samples are shown for plane-polarized light (PPL), cross-polarized light (XPL), and backscatter electron (BSE). Ol is olivine, Cpx is clinopyroxene, Mgt is magnetite, and Spl is spinel. The bathymetric images were drawn using the GMT (Wessel et al., 2019). Photos of seafloor lava outcrops were provided by the cruise report of YK16-01, YK18-08, and YK19-05S cruises in the Data and Sample Research System for Whole Cruise Information (DARWIN) by JAMSTEC ( last access: 18 December 2023).

3.1 YK16-01 cruise and 6K#1466 dive

During the YK16-01 cruise, a small conical knoll (ca. 0.04 km3) was investigated by a submersible dive, namely 6K#1466 (Figs. 2 and 3a). The lava flows, which were observed in a hollow lava tube resulting in a sediment-rolling/disturbing eruption, were located  600 m south of the top of the knoll, featuring extremely fresh and glassy samples (6K#1466R3-001 and R3-004 basalts) (Fig. 3a). Vesicular pillow basalts were collected on the western slope of the knoll (samples 6K#1466R6-001, R7-001, and R7-003; Fig. 3a). While the strong acoustic reflection could not entirely distinguish the petit-spot lava fields in ferromanganese nodule fields, the 6K#1466 dive revealed lava outcrops using a sub-bottom profiler (SBP) and a multinarrow beam echo sounder (MBES). Specifically, the petit-spot lava field, as an acoustically opaque layer, exhibited a vigorous backscattering intensity in the MBES, along with the distributions of the basement and sediment layers in the SBP.

The 6K#1466R3-001 and R3-004 samples were extremely fresh glassy basalts. The samples exhibited similar petrographic features (Fig. 3a). These samples were enveloped by a 3.0–4.5 mm thick palagonite layer (hydrated quenched glass), with their outermost parts being surrounded by a 5.8–7.2 mm thick ferromanganese crust (Fig. 3a). They were less vesicular (< 3 vol. %) and were dominantly composed of basaltic glass, euhedral–subhedral olivine microphenocrysts ( 100–500 µm in size), ferrotitanium oxide (< 50 µm in size), and minor plagioclase ( 500 µm in size) (Fig. 3a). No secondary phases such as clay minerals were observed.

The 6K#1466R6-001, R7-001, and R7-003 basalts, which were covered with a 4.3–5.2 mm thick ferromanganese crust over 5.5–6.6 mm thick palagonite rinds, exhibited high vesicularity (20 vol. %–40 vol. %) (Fig. 3a). Mikuni et al. (2022) reported certain pyroxene-dominated xenocrysts and peridotite xenoliths. The basaltic groundmass was characterized by needle-shaped clinopyroxene (50–400 µm in size), subhedral olivine that was partly with aureoles of iddingsite (up to 100 µm in size), ferrotitanium oxide, minor spinel (up to 10 µm in size), glass, and crystallite, notably without remarkable phenocrysts (Fig. 3a). The photomicrograph of R6-001 is shown in Fig. 3a.

3.2 YK18-08 cruise and 6K#1521 and 6K#1522 dives

Two submersible dives (6K#1521 and #1522) were conducted during the YK18-08 cruise to investigate petit-spot volcanoes. During the 6K#1521 dive, a small lava outcrop was identified in the abyssal plain by tracing a strong acoustic reflection, which was expected to originate from intrusive rock bodies, in the sedimentary layer detected by deep-sea SBP equipped on Shinkai 6500. The strong reflective surface gradually became shallow during the navigation, revealing the small lava outcrop (Figs. 2 and 3b). Fresh and massive (nonvesicular) basalts were collected from this outcrop (samples 6K#1521R04 and R05; Fig. 3b). The samples obtained from the 6K#1522 dive at a seamount exhibited highly irregular shapes, and massive lava flows, pillows, and lava breccia were observed (Fig. 3c). All the samples were fresh vesicular basalts (6K#1522R01, R02, R05, R12, R13, R16, and R17; Fig. 3c).

The fresh, massive, and nonvesicular basalts were collected during the 6K#1521 dive (R04 and R05) and comprised euhedral olivine microphenocrysts (150–400 µm in size), two types of ferrotitanium oxide (50–150 µm in size), and crystallite (Fig. 2b). Secondary phases were not observed. They were covered with a 5.6–5.9 mm thick ferromanganese crust and a  1.0 mm thick palagonite rind (Fig. 3b); however, R05 did not have palagonite rinds. The photomicrograph of R04 is shown in Fig. 3b.

The seven fresh basalts collected during the 6K#1522 dive (6K#1522R01, R02, R05, R12, R13, R16, and R17), exhibited high vesicularity (20 vol. %–40 vol. %) with 2.9–6.0 mm thick palagonite rinds covered with 2.7–5.9 mm thick ferromanganese crusts (Fig. 3c). Euhedral–subhedral olivine microphenocrysts (glomeroporphyritic; 30–200 µm in size), radial–needle-shaped clinopyroxene, iddingsite (< 200 µm in size), spinel, and glass with minor xenocrystic olivines were observed (Fig. 3c). The photomicrograph of R01 is shown in Fig. 3c.

3.3 YK19-05S cruise and 6K#1542 and 6K#1544 dives

A petit-spot knoll and associated lava flows were investigated by the 6K#1542 and #1544 dives during the YK19-05S cruise (Fig. 2). During the 6K#1542 dive, geological survey and rock sampling were conducted from two points on the eastern slope of the knoll (Figs. 2 and 3d). The 6K#1542R03 and R05 basalts were collected from the lava–breccia field covered with a thin ferromanganese crust (Fig. 3d). Additionally, samples R06 and R09 were obtained from the lobate-surface lava between tubular lavas closer to the summit than R03 and R05 (Fig. 3d).

High-resolution (1 m scale) bathymetric mapping was successfully conducted during the 6K#1544 dive, which can contribute to future oceanographic investigations using a human-occupied vehicle (Kaneko et al., 2022). Several mounds, 10–20 m in height and a few hundred meters in diameter, were recognized during this acoustic survey (Fig. 3d). We observed these mounds and collected samples from outcrops during the second half of the dive. Furthermore, pillow lavas, tumuli, and lava–breccia were observed, and basaltic samples (6K#1544R04, R05, and R06) were collected (Fig. 3d).

Four vesicular basalts (10 vol. %–30 vol. % vesicularity; 6K#1542R03, R05, R06, and R09) were covered with 4.3–4.4 mm thick ferromanganese crusts. The outer palagonitic rinds were 3.2–3.4 mm thick (Fig. 3d). Euhedral–subhedral olivine microlites (up to sizes of 300,µm) and microphenocrysts were glomeroporphyritic (Fig. 3d). The groundmass was dominated by needled dendritic clinopyroxenes ( 100 µm in size), along with olivine, spinel, glass, and xenocrystic olivine megacrysts. The photomicrograph of R06 is shown in Fig. 3d.

Basaltic samples from the 6K#1544 dive (6K#1544R04, R05, and R06) were covered with ferromanganese crust (5.0–5.8 mm thick) over palagonitic rinds (3.4–4.4 mm thick). All the samples exhibited high vesicularity in the range of 20 vol. %–35 vol. % (Fig. 3d). They comprised olivine microphenocrysts (30–250 µm in size; euhedral–subhedral or columnar), clinopyroxene (< 100 µm; needled, columnar, radial or dendritic shape), spinel, and glass without secondary phases (Fig. 3d).

The photomicrograph of R04 is shown in Fig. 3d. During macroscopic observations, practically all the basalts from the 6K#1542 and 6K#1544 dives exhibited similar vesicularity and freshness. Their geochemical features were also similar to each other and are described in Sects. 5.1 and 5.2.

4 Analytical methods

4.1 Major and trace element analysis of volcanic glass, mineral, and whole rock

Major element compositions of glasses and minerals were determined using an electron probe micro-analyzer (EPMA). JXA-8900R at the Atmosphere and Ocean Research Institute (AORI), the University of Tokyo, was used for glass analysis and JXA-iHP200F at the Geological Survey of Japan (GSJ), National Institute of Advanced Industrial Science and Technology (AIST), was used for mineral analysis. The analyses were performed using an accelerating voltage of 15 kV, a beam current of 12 nA, and a beam diameter of 10 µm for glass and 2 µm for mineral. A peak counting time of 20 s and a background counting time of 10 s were used, except for Ni, for which a peak counting time of 30 s and a background counting time of 15 s were used. For Na analysis of glass, the peak counting time was 5 s, and the background counting time was 2 s. Natural and synthetic minerals were used as standards, and data were corrected using a ZAF online correction program (Akizawa et al., 2021). Major element composition of glass was determined by the mean value of 10 analytical points.

Trace element compositions of minerals were determined using a laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS; New Wave Research UP213 and Agilent 7500s) at Kanazawa University. The Nd:YAG deep UV (ultraviolet) laser's wavelength is 213 nm. The analyses were conducted with 100,µm spot size. A repetition frequency of 6 Hz and a laser energy density of 8 J cm−2 were used. NIST612 glass (distributed by National Institute of Standards and Technology) was employed for calibration, using the preferred values of Pearce et al. (1997). Data reduction was undertaken with 29Si as the initial standard, and SiO2 concentrations were obtained by an electron microprobe analysis (Longerich et al., 1996). BCR-2G (distributed by the United States Geological Survey) was used as a secondary standard to assess the precision of each analytical session (Jochum and Nohl, 2008).

Whole-rock major and trace element compositions of rock samples were analyzed by Activation Laboratories Ltd., Canada, using the Code 4Lithoresearch Lithogeochemistry and ultratrace5 Exploration Geochemistry Package. The former package uses lithium metaborate/tetraborate fusion with inductively coupled plasma–optical emission spectrometry (FUS-ICP-OES) and inductively coupled plasma–mass spectroscopy (FUS-ICP-MS) for the major and trace element analyses, respectively. The latter package uses inductively coupled plasma–optical emission spectrometry (ICP-OES) and inductively coupled plasma–mass spectroscopy (ICP-MS) for the major and trace element analyses, respectively.

4.2Sr, Nd, and Pb isotope analysis

4.2.1 Acid leaching

Acid leaching was conducted for the selected basaltic samples on the basis of the procedure of Weis and Frey (1991, 1996) as follows: (1) about 0.3–0.4 or 0.6 g of rock powder is weighed into an acid-washed 15 mL Teflon vial (Savillex®). (2) Then 10 or 12 mL of 6 N (N is for normality) HCl were added and heated at 80 C for 20–30 min. (3) After heating, the suspension is ultra-sonicated in 60 C water for 20 min. (4) The supernatant is decanted. Steps (2) to (4) were repeated more than four times (up to six times) until the supernatant become clear or pale yellow to colorless. (5) TAMAPURE-AA ultrapure water (Tama Chemicals Co., Ltd.), which includes a lower Pb blank than Milli-Q H2O, was added instead of 6 N HCl, and the suspension is ultra-sonicated for 20 min. This step is conducted twice. (6) The leached rock powder is dried on a hot plate at 120 C. (7) After cooling, the powder is weighed.

4.2.2 Extraction of Pb, Sr, and Nd

The extraction of Pb, Sr, and Nd was performed, following the procedures of Tanimizu and Ishikawa (2006) and Machida et al. (2009). First,  50 to  100 mg of rock powder was weighted in a 7 mL Teflon vial (designated as vial A) and digested using mixed acid composed of HF and HBr. The separation was conducted by cation exchange resin (AG-1X8; Bio-Rad Laboratories, Inc.) on the basis of procedures described in Tanimizu and Ishikawa (2006). All fractions from the first and second supernatant loading (0.5 M HBr) to the elution of other elements (mixed acid composed of 0.25 M HBr and 0.5 M HNO3) were collected in another 7 mL Teflon vial (designated as vial B) for Sr and Nd separation. Finally, Pb was extracted by 1 mL of 1 M HNO3 in another 7 mL Teflon vial (designated as vial C). The procedural blanks for Pb totaled less than 23 pg.

The Sr- and Nd-bearing solution in vial B was transferred into vial A and contained residues of digested samples. A total of 2 mL of HClO4 and 2 mL HNO3 was further added to vial A, and the residue was dissolved at 110 C. Both Sr and Nd were separated by column with a cation exchange resin (AG50W-8X; Bio-Rad Laboratories, Inc.) and a Ln resin (Eichrom Technologies LLC) on the basis of the procedures described in Machida et al. (2009). The separated Sr and Nd were further purified by column separation with a cation exchange resin. The total procedural blanks for Sr and Nd were less than 100 pg.

4.2.3 Analytical procedure

Pb isotopic ratios were obtained using the multicollector ICP-MS (MC-ICP-MS; NEPTUNE Plus, Thermo Fisher Scientific), with nine Faraday collectors, at the Chiba Institute of Technology (CIT), Japan. The NIST SRM-981 Pb standard was also analyzed and yielded average values of 206Pb/204Pb= 16.9303 ± 0.0005, 207Pb/204Pb= 15.4828 ± 0.0006, and 208Pb/204Pb= 36.6710 ± 0.0016. These correspond to previous values determined using MC-ICP-MS with Tl normalization, but they were slightly lower than values determined by thermal ionization mass spectrometry (TIMS) in Tanimizu and Ishikawa (2006) from the 207Pb204Pb double spike. Reproducibility was monitored by an analyses of the JB-2 GSJ standard, and the obtained values were 206Pb/204Pb= 18.3326 ± 0.0005, 207Pb/204Pb= 15.5453 ± 0.0006, and 208Pb/204Pb= 38.2240 ± 0.0017.

Sr and Nd isotopic analyses for powdered rocks and glasses were conducted using TIMS (Triton XT, Thermo Fisher Scientific), with nine Faraday collectors, at CIT. A total of 1.5 µL of 2.5 M HCl and 0.5 M HNO3 was used for loading of separated Sr and Nd of sample on the single- and double-Re filament, respectively. The measured isotopic ratios were corrected for instrumental fractionation by adopting the 86Sr/85Sr value to be 0.1194 and that of 146Nd/144Nd to be 0.7219. The average value for the NIST SRM-987 Sr standard was 0.710239 ± 0.000005 (2σ, n=2) and that for the GSJ JNdi-1 Nd standard was 0.512103 ± 0.000005 (2σ, n=2). They agree well with values from the literature for the NIST SRM-987 (87Sr/86Sr= 0.710252–0.710256; Weis et al., 2006) and JNdi-1 (143Nd/144Nd= 0.512101; Wakaki et al., 2007). Consequently, we did not correct the values of the unknowns for offsets between the measurements and the values for the Sr and Nd standards.

4.340Ar/39Ar dating

Samples for 40Ar/39Ar dating were prepared by separating crystalline groundmass after crushing them to sizes between 100 and 500 µm. The separated groundmass samples were leached by HNO3 (1 mol L−1) for 1 h to remove clays and altered materials. All samples were wrapped in aluminum foil, along with JG-1 biotite (Iwata, 1998), K2SO4, and CaF2q flux monitors. Anything amorphous (e.g., quenched glass) was removed because 39Ar may move from one phase to another in a process known as recoil. This can create a disturbed age spectrum when 39Ar is produced from 39K in amorphous material through interaction with fast neutrons during the irradiation of the sample. Samples were irradiated for 6.6 d in the Kyoto University Research Reactor (KUR), Kyoto University. Argon extraction and isotopic analyses were undertaken at the Graduate School of Arts and Sciences, the University of Tokyo. The sample gases were extracted by incremental heating of 10 or 11 steps between 600 and 1500 C. The analytical methods used are the same as those used by Ebisawa et al. (2004) and Kobayashi et al. (2021).

5 Results

To describe the geochemical and chronological results, each sample group was denoted by its dive number; e.g., the sample group obtained from the 6K#1521 dive was labeled 1521 samples or basalts. The basalts from the 6K#1466 dive were divided into two groups for R3 (collected from the seafloor south of the knoll) and R6–R7 (sampled on the knoll), based on their geographical, petrological, and compositional differences. The mineral compositions of each petit-spot basalt are shown in Fig. S1 in the Supplement and Tables S1–S3 in the Supplement.

Table 2Major element compositions of western Pacific petit-spot basalts.

FeOT as total values. Mg#= 100 ×Mg/[Mg+Fe2+]molar. “–” not detected.  Analyzed by Activation Laboratories Ltd. (Actlabs).

Download Print Version | Download XLSX

Table 3Trace element compositions of western Pacific petit-spot basalts.

“–” not detected.  Analyzed by Activation Laboratories Ltd. (Actlabs).

Download Print Version | Download XLSX

Figure 4Relationships between the SiO2 and alkali contents. (a) Total alkali vs. silica diagram, using the platform of Le Bas et al. (1986). The dividing line of alkaline and sub-alkaline is from Irvine and Baragar (1971). The data are plotted as the total 100 wt %. The triangles and circles show the whole-rock and quenched-glass compositions, respectively. The compositions of the NW Pacific petit-spots are represented by gray triangles (Hirano and Machida, 2022). The data of the 1203 and 1206 basalts are from Hirano et al. (2019), and those of the 1466R7 basalts are from Mikuni et al. (2022). (b) K2O vs. Na2O diagram. The maximum K2O/Na2O value of kimberlite is from PetDB database (, last access: 15 March 2022). The data of ocean island basalt (OIB) and MORB are compiled from Stracke et al. (2022) as “Expert datasets” in the GEOROC database (, last access: 31 October 2023).

5.1 Major and trace element compositions

The major and trace element compositions for the whole rock and glass of the petit-spot basalts are listed in Table 2 and 3, respectively. The basalt compositions for a petit-spot knoll were reported by Hirano et al. (2019) (expressed as 1203 and 1206 in each figure). The data are discussed, along with the reported NW Pacific petit-spots (Hirano and Machida, 2022). Using a total alkali vs. silica (TAS) diagram, virtually all the samples were classified as alkalic rocks, but the 1542 and 1544 basalts were plotted near the boundary between alkalic and non-alkalic (Fig. 4a). Two petit-spot basalts (1466R7-001 and R7-003) from the petit-spot knoll were notably silica-undersaturated (i.e., SiO2= 39.3–39.4 wt %) and classified as foidite (Mikuni et al., 2022). All the western Pacific petit-spot basalts, except for the 1466R7 basalts, were sodic (K2O/Na2O= 0.24–0.58) and were notably distinguished from the potassic NW Pacific petit-spots (Fig. 4b).

Figure 5Selected major-element oxides against MgO. The symbols and compiled data correspond to those in Fig. 3.


Figure 6Selected trace element ratios against MgO. The symbols and compiled data correspond to those in Fig. 3.


Selected major element oxides and trace element ratios vs. MgO plots for the petit-spot basalts are shown in Figs. 5 and 6, respectively. The MgO concentrations of the 1466R3 and 1521 samples each exhibiting similar petrographic features (i.e., nonvesicular and glassy) were characterized by values (4.0–4.4 wt %) lower than those of other vesicular samples (6.6–9.3 wt %). The K2O, Na2O, Al2O3, and SiO2 contents negatively correlated with MgO (Fig. 5a–d). The CaO, FeOT, and CaO/Al2O3 abundances exhibited positive correlations with MgO (Fig. 5e–g). The TiO2 concentrations exhibited no correlations with MgO (Fig. 5h), as well as the selected trace element ratios (Fig. 6a–g), except for the Sm/Hf ratio with positive correlations (Fig. 6h). The Sm/Hf ratio also negatively correlated with SiO2 (Fig. S2). The study samples exhibited whole-rock loss on ignition (LOI) in the range of 0.67–1.72 wt %, excluding two relatively altered samples, 1466R7-001 (LOI = 2.68 wt %) and R7-003 basalts (LOI = 6.29 wt %).

Figure 7Primitive mantle (PM)-normalized (Sun and McDonough, 1989) trace element patterns (a–g) and element ratios (h). (g) The compositional range of the study samples and NW Pacific petit-spots (Hirano and Machida, 2022). (h) The Ba/Nb and Sm/Hf ratios of the petit-spot basalts to distinguish the three groups, after Machida et al. (2015). The data of 1203 and 1206 basalts and 1466R7 basalts are from Hirano et al. (2019) and Mikuni et al. (2022), respectively. The symbols and compiled data in panel (h) correspond to those in Fig. 3.

Figure 8Alteration sensitive elements (Ba and U) vs. insensitive elements (Nb and Th). The symbols and compiled data correspond to those in Fig. 3.


The primitive mantle (PM)-normalized (Sun and McDonough, 1989) trace element patterns for the petit-spot basalts, including those reported by a previous study (Hirano et al., 2019), were shown for each dive compared to the representative ocean island basalt (OIB) in Fig. 7a–f. The petit-spot basalts generally showed high light rare earth element (LREE) or heavy rare earth element (HREE) ratios. Negative Zr, Hf, Ti, and Y anomalies were commonly observed in these western Pacific petit-spots, as well as those of the NW Pacific petit-spots (Fig. 7g). The 1466 basalts collected on the seafloor south of the knoll (1466R3-001 and 1466R3-004 basalts) were compositionally different from those obtained on the knoll (1466R7-001 and 1466R7-003 samples). The basalts from the 6K#1542 and #1544 dives, collected from nearby locations, had the same compositions in major and trace element ratios in both whole rock and glass, respectively (Figs. 4–6, and 7e and f). These samples in the Ba/Nb and Sm/Hf diagrams were plotted in the range of Group 3 in the discrimination of the NW Pacific petit-spot basalts (Machida et al., 2015), indicating their negative Zr and Hf anomalies without notable U, Th, Nb, and Ta anomalies in the PM-normalized trace element patterns (Fig. 7h). The Sm/Hf ratio of the differentiated 1466R3 samples was lower than that of other samples. A positive correlation between fluid mobile and immobile elements, Ba vs. Nb (Fig. 8a) and U vs. Th (Fig. 8b), respectively, was observed, excluding the Ba of the 1466R7 samples (Fig. 8a).

Figure 9SrNdPb isotopic variations in the petit-spot basalts. The mantle endmembers are derived from a study by Zindler and Hart (1986). The open triangles (a–d) represent the acid-leached samples. Carbonatite data were compiled from GEOROC (, last access: 1 November 2023) with Bizimis et al. (2003). Carbonatite data with 87Sr/86Sr > 0.706 by GEOROC were eliminated. The northwestern (NW) Pacific petit-spots and petit-spots off the Tonga Trench are from Hirano and Machida (2022) and Reinhard et al. (2019), respectively. The petit-spots off the Java Trench are from Taneja et al. (2016) and Falloon et al. (2022). The data of 1203 and 1206 basalts are from Hirano et al. (2019). The data of the Marcus–Wake Seamount chains are from studies by Konovalov and Martynov (1992), Koppers et al. (2003), Konter et al. (2008), Natland (1976), Smith et al. (1989), and Staudigel et al. (1991). The Northern Hemisphere reference line (NHRL) and low-Nd (LoNd) arrays are from studies by Hart (1984) and Hart et al. (1986), respectively. (e) The three-dimensional (3D) plot of the SrNdPb isotopic compositions. The compilation and mantle endmembers correspond to panels (a)(d). The color usages of the plots were the same as panels (a)(d).

Table 4Sr, Nd, and Pb isotopic compositions of western Pacific petit-spot basalts and measured standards.

Errors shown in parentheses represent 2σ and apply to the last two digits.

Download Print Version | Download XLSX

5.2SrNdPb isotopic composition

The Sr, Nd, and Pb isotopic compositions of the leached and unleached whole rock and fresh glasses in this study (presented in Table 4) were in practically identical ranges of 87Sr/86Sr (0.703412–0.704424), 143Nd/144Nd (0.512694–0.512890), 206Pb/204Pb (18.6582–18.7778), 207Pb/204Pb (15.5086–15.5749), and 208Pb/204Pb (38.6506–38.8041) despite their different locations (Fig. 9a–d; Table 4). The isotopic compositions of the quenched glass and whole rock were identical, indicating that the characteristics of the melting source could be obtained through the geochemistry of the young and fresh volcanic quenched glass. The leached and unleached materials of the same sample also had similar isotopic ratios, except for the 1466R7-003 basalt, which had a relatively high LOI (6.29 wt %) (Fig. 9a–d). The SrNdPb isotopic three-dimensional (3D) plot is shown in Fig. 9e.

Figure 10Geochronological data. (a) The 40Ar/39Ar ages of the 6K#1466R6-001 and 6K#1522R01 basalts. The errors show a 2σ confidence level. (b) Estimated relative ages using the thickness of ferromanganese crust (green bands) and palagonite (hydrated quenched-glass rind; red bands) covered with petit-spot basalts. These values were estimated using the average for each sample (n= 20). The UPb age of zircon in the 6K#1203 and 1206 peperites are from Hirano et al. (2019).

5.3 Age determination and estimation

The 40Ar/39Ar ages were determined for two samples (1466R6-001 and 1522R01) (Fig. 10a; Table S4). The secondary material (e.g., alteration products) plausibly causes the recoil loss and redistribution of Ar during irradiation of samples, particularly the fine-grained groundmass separates of submarine basalt (Koppers et al., 2000). This effect is negligible for 40Ar/39Ar dating samples in this study because the total K/Ca ratios estimated using the irradiated 39ArK/37ArCa ratio (0.089 for 1466R6, 0.080 for 1522R01; Table S4) are mostly correspond to the bulk K/Ca ratios calculated using the major element compositions of Table 2 (0.088 for 1466R6-001; 0.076 for 1522R01). This is supported by the rock descriptions recognizing no secondary materials of crystalline 40Ar/39Ar specimens. The 1466R6-001 sample had a plateau age of 3.03 ± 0.18 Ma in seven fractions comprising 94.1 % released 39Ar. However, the plateau age was recognized as apparently old, owing to excess 40Ar, as indicated by the initial 40Ar/36Ar ratio of 325 ± 15, which exceeded the atmospheric ratio (296.0; Nier, 1950) in the inverse isochron. The inverse isochron age of 2.56 ± 0.34 Ma showed the best age estimate for the 1466R6-001 basalt (Fig. 10a). The 1522R01 sample released almost no radiogenic daughter nuclide of 40Ar in the KAr age system (Fig. 10a).

The ranges of eruption age were estimated for all the samples, using the average thickness (n=20) of ferromanganese crust and palagonite rind (hydrated quenched glass) with their deposition/formation rates on the seafloor (ferromanganese crust; 1–10 mm Myr−1; Hein et al., 1999; palagonite, 0.03–0.3 mm Myr−1; Moore et al., 1985) (Fig. 10b). Using this approach, the western Pacific petit-spots were expected to have erupted later than ca. 9 Ma. The ranges of eruption age estimated from palagonite rind did not overlap with those from ferromanganese crust showing older durations, although they had general correlations (Fig. 10b). The 40Ar/39Ar ages of two samples and the UPb age of zircon in the 1203 and 1206 peperites (Hirano et al., 2019) were overlaid within these ranges.

6 Discussion

6.1 Eruptive setting of western Pacific petit-spots

In this study, two crystalline petit-spot basalts were subjected to 40Ar/39Ar dating. A previously investigated petit-spot knoll in this region (examined during the 6K#1203 and #1206 dives) was dated at younger than 3 Ma through the UPb dating of eight zircons in peperites (Fig. 10b) (Hirano et al., 2019). The results revealed that the silica-undersaturated vesicular basalt of 1466R6-001, hosting ultramafic xenoliths (Mikuni et al., 2022), exhibited a 40Ar/39Ar age of 2.56 ± 0.34 Ma (Fig. 10). On the contrary, the fresh vesicular basalt of 1522R01, which erupted at the foot of the 100 Ma Takuyo–Daigo Seamount (Fig. 2) (Nozaki et al., 2016), did not exhibit radiogenic 40Ar to indicate its young age ( 0 Ma) (Fig. 10). The ranges of eruption ages were estimated using the average thickness of ferromanganese crust and palagonite rind (seawater-hydrated quenched glass) with their deposition/formation rates on the seafloor. The 40Ar/39Ar and zircon UPb ages were within these ranges (Fig. 10). The petit-spot volcanic field is surrounded by Cretaceous seamounts (Koppers et al., 2003) and irregular Paleogene volcanoes (Aftabuzzaman et al., 2021; Hirano et al., 2021). However, no zero-aged hotspots were observed in this region, and the P-wave tomographic image of the surface to the core–mantle boundary of the study area did not exhibit a plume-like low-velocity zone (Fig. 1c; Lu et al., 2019). Furthermore, the MORB-like to more depleted noble-gas isotopic compositions of the petit-spot knoll (investigated by 6K#1203 and #1206 dives) suggested its upper-mantle origin (Yamamoto et al., 2018). Along with the outer-rise bulge in front of the Mariana Trench detected through a positive gravitational anomaly (Hirano et al., 2019), these data suggest that the western Pacific petit-spot volcanoes could have erupted at  0–3 Ma, owing to the flexure of the subducting Pacific Plate into the Mariana and Izu–Ogasawara trenches.

The petit-spot basalts from the 6K#1542 and #1544 dives could have originated from the same eruptive source, based on their similar petrographic and geochemical features, despite a distance of  6.8 km between both (Figs. 3d and 4–9). Contrarily, in terms of their petrography and geochemistry, the basalts from the 6K#1466 dive are distinguished between the samples from the lava flows on the abyssal plain (1466R3-001 and 1466R3-004 samples) and the samples from the knoll site (1466R6-001, 1466R7-001, and 1466R7-003 samples). The 1466R3 basalts were collected at a lava outcrop 600 m south of the knoll, and the 1466R6 and 1466R7 samples were collected on the western slope of the knoll (Fig. 3a). The 1466R3 series are glassy, with a high SiO2 content (50.6–51.6 wt %), including minor plagioclase and fewer vesicles (Figs. 3a and 4a). However, the 1466R6–R7 series exhibited silica-undersaturated compositions (SiO2= 39.3–39.4 wt %) and high vesicularities (20 vol. %–40 vol. %) (Figs. 3b and 4a). Combining these observations with the differences in MgO contents and trace element compositions, the 1466R3 and 1466R6–R7 basalts are implied to have different parental magmas (Figs. 6 and 7b). Generally, vesicular samples (1203, 1206, 1466R7, 1522, 1542, and 1544 basalts) are relatively primary (i.e., MgO > 6.63 wt %), whereas nonvesicular samples (1466R3 and 1521 basalts) are evolved (i.e., MgO < 4.43 wt %). This correlates with the compositions of olivine microphenocrysts in the low-forsterite content (Fo# = 100 ×Mg/[Mg+Fe2+]cation) of olivine in evolved basalts and the high Fo# of olivine in the relatively primary basalts (Fig. S1a–c).

The CI chondrite-normalized REE ratios of these samples are within those of OIBs, and the REE patterns exhibit HREE-depleted patterns (Fig. S3). However, among the western Pacific petit-spots, each volcano shows distinct REE and trace element ratios (i.e., parental magmas) (Figs. 6 and S3). Considering the absence of correlation between MgO and the trace element ratios, it is suggested that each volcano could have originated from isolated sources (i.e., melt ponds) with varying chemical compositions and degrees of melting (Fig. 6). On the contrary, the radiogenic Sr, Nd, and Pb isotopic ratios of the samples are nearly identical, indicating equivalent components in the source (Fig. 9).

In summary, (1) the western Pacific petit-spot volcanoes erupted at  0–3 Ma, owing to the plate flexure related to the subduction of the Pacific Plate into the Mariana Trench (Figs. 1 and 2). (2) The 1542 and 1544 samples originated during the same magmatic event (Fig. 3d). However, the basalts from the 6K#1466 dive were divided into two parental magmas (1466R3 and 1466R6–R7 basalts) (Fig. 3a). (3) Each volcano originated from an isolated source and/or ascending processes, as indicated by independent trace element ratios. Despite this, the geochemical components involved in the source were similar among the western Pacific petit-spot volcanoes due to the nearly identical Sr, Nd, and Pb isotopic compositions (Figs. 6 and 9). The variation in trace element compositions among the volcanoes is plausibly attributed to the degree of contribution of carbonatite flux and/or the recycled crustal component to the source, as discussed below.

6.2 Petit-spot magma composition and its evaluation

Post-eruption alteration in seawater may have affected the chemical composition of oceanic basalts. Thus, various approaches, including petrographic observation, geochemical investigation, and acid leaching, have been employed to evaluate the primary features and the removal of this effect for isotopic analysis (Hanano et al., 2009; Melson et al., 1968; Miyashiro et al., 1971; Nobre Silva et al., 2009; Resing and Sansone, 1999; Staudigel and Hart, 1983; Zakharov et al., 2021). The study samples exhibit whole-rock LOI of < 1.72 wt %, except for two relatively altered samples, 1466R7-001 (LOI = 2.68 wt %) and R7-003 basalts (LOI = 6.29 wt %). Pristine quenched glasses are preserved in most of the samples, excluding three exceptional samples (1466R6-001, R7-001, and R7-003 basalts). Positive correlations exist between the alteration-insensitive (e.g., Nb and Th) and alteration-sensitive (e.g., Ba and U) incompatible elements, indicating that the effect of seawater alteration was not extensive, except for the 1466R7-001 and R7-003 basalts (Fig. 8). Despite originating from different volcanic edifices, the positive correlation of all the study samples is attributed to the chemical similarity of source compositions for certain elements (i.e., the Ba/Nb and U/Th ratios are nearly constant among the samples), as well as the Sr, Nd, and Pb isotopic compositions (Fig. 9). These findings demonstrate that most of the petit-spot basalts were largely unaffected by seawater alteration, with a few exceptions, i.e., 1466R7-001 and R7-003 basalts.

The MgO (4–9 wt %), Ni (< 263 ppm), and Cr (< 350 ppm) contents in the samples are lower than the expected values of primary mantle-derived melt (MgO > 10 wt %, Ni > 400 ppm, Cr > 1000 ppm; Frey et al., 1978). Similarly, the Mg# (100 ×Mg/[Fe2++Mg]molar) values range from 41 to 57 (Table 2) against the primary basaltic melt, which is equilibrated with the upper mantle (Mg#= 66–75; Irving and Green, 1976). No phenocrysts were observed (only microphenocrysts), despite such differentiated compositions, as well as most of the NW Pacific petit-spot basalts. This suggests that the western Pacific petit-spots experienced crystal fractionation in the lithosphere, as well as the case in the NW Pacific petit-spot (Machida et al., 2017; Valentine and Hirano, 2010; Hirano, 2011; Yamamoto et al., 2014). Consequently, calculating the primary composition of the petit-spot basalts using the mineral modal composition on the thin section was not possible. However, the major element trends of the samples indicate the crystal fractionation of the same phases. Negative trends of the Al2O3 content and the positive trends in CaO and CaO/Al2O3 content with decreasing MgO indicate the occurrence of olivine, spinel, and clinopyroxene fractionation (Fig. 5c, e, and g). The absence of visible correlations of K2O, Na2O, SiO2, and TiO2 contents against MgO suggests insignificant fractionation of plagioclase and FeTi oxides. The FeTi oxides as minor phases in the groundmasses and plagioclases were only observed in the most differentiated 1466R3-001 and R3-004 basalts (Figs. 3 and 5a, b, d, and h). However, these major elemental trends should be interpreted as apparent because each petit-spot volcano originated from an isolated parental magma with a different chemical composition or degree of partial melting, as discussed above.

The melting source of alkali basalts can be determined more effectively by examining their trace element composition rather than major elements (Hofmann, 2003; Machida et al., 2014, 2015). Trace element composition of magma, however, could be modified by crustal and/or mantle assimilation and fractionation of specific minerals. The relatively primitive basalts (1203, 1206, 1466R6, R7, 1522, 1542, and 1544 samples) contained xenocrystic olivines and partly ultramafic xenoliths, suggesting a rapid magma ascent (Hirano et al., 2019; Mikuni et al., 2022; Fig. S4). However, since the stagnation of ascending petit-spot magma could lead to the formation of fertile peridotite and pyroxene-rich veins in the middle to lower depths of the lithosphere (Mikuni et al., 2022; Pilet et al., 2016), the chemical composition of the petit-spot magma could be modified through assimilation with ambient lithospheric peridotite. According to Hirano and Machida (2022), ascending silica-undersaturated melt would predominantly consume orthopyroxene (± spinel) and result in a more silicic composition with Zr and Hf depletion. This is due to the relatively higher ZrHf partition of orthopyroxene compared to other trace elements (Pilet et al., 2008; Shaw, 1999; Tamura et al., 2019). The orthopyroxenes of fertile pyroxenites and lherzolite xenoliths metasomatized by petit-spot melts exhibit Zr and Hf enrichment (Mikuni et al., 2022; Fig. S5). If this silica-enrichment (i.e., melt–rock interaction) were significant, then a positive correlation between SiO2 and Sm/Hf would be expected as a mantle assimilation trend. However, the samples exhibited a negative correlation, similar to those of the NW Pacific petit-spots (Hirano and Machida, 2022) (Fig. S2). Considering the relation between the Sm and Hf partition coefficients of clinopyroxene (i.e., DHf <DSm; McKenzie and O'Nions, 1991; Kelemen et al., 2003), we suggest that the negative correlation between the Sm/Hf and SiO2 in the petit-spot basalts probably reflects the crystal fractionation of clinopyroxene rather than mantle assimilation. The Ba/Nb ratios of the samples are nearly constant and do not correlate with the MgO and SiO2 contents (Figs. 6g and S2g). The lack of correlation between other trace element ratios, excluding Sm/Hf and Ba/Nb (i.e., La/Y, La/Lu, Sm/Yb, La/Sm, Nb/Ta, Zr/Hf), and the MgO concentration suggests that crystal fractionation may not have been involved in those of the incipient melt (Fig. 6). However, independently tracking the evolution of the trace element composition for each volcano is challenging, given that each volcano originated from isolated sources. Thus, considering the observations above, the fresh and zero-aged 1522 basalts (having the highest Sm/Hf ratios and lowest SiO2 contents among the fresh samples and higher MgO contents) were selected for further analysis with geochemical modeling. Given that the 1522 samples had MgO in the range of 6.63–7.36 wt %, olivine was expected to be the dominant phase of crystal fractionation (Asimow and Langmuir, 2003; Helz and Thornber, 1987; Herzberg, 2006). By applying the olivine maximum fractionation model (Takahashi, 1986; Tatsumi et al., 1983) to test two samples, it was noted that 7 %–9 % olivine addition was required to achieve the olivine composition, corresponding to the mantle olivine array in the NiO and Fo# spaces (Fig. S6a and b). The calculated primary trace element contents did not considerably differ from those of the analytical compositions (Table S5 and Fig. S6). Thus, the 1522 basalts were assumed to be the most primary petit-spot basalt samples and were used to evaluate the geochemical modeling results.

6.3 Melting source of western Pacific petit-spots

The depletions observed in specific elements (e.g., Ta, Zr, Hf, and Ti) in the petit-spot basalts potentially demonstrate the involvement of carbonatitic materials in conjunction with a large amount of CO2 and lower Mg isotopic ratio than that of the normal mantle (Bizimis et al., 2003; Dasgupta et al., 2009; Hirano and Machida, 2022; Hoernle et al., 2002; Liu et al., 2020; Okumura and Hirano, 2013). Other oceanic lavas originating from the asthenosphere (e.g., Hawaiian rejuvenated lavas and North Arch volcanoes) exhibited characteristic trace element signatures (i.e., Zr and Hf depletion), similar to those of petit-spot lavas. This implies that their melting sources were involved with carbonatitic materials with or without plume-derived components (Fig. S7; Borisova and Tilhac, 2021; Clague and Frey, 1982; Clague et al., 1990; Dixon et al., 2008; Yang et al., 2003). Additionally, the involvement of recycled crustal components was inferred from the geochemical features of the petit-spot basalts, and the upper mantle was revealed to be heterogeneous (Liu et al., 2020; Machida et al., 2009, 2015). Such a scenario of the source for petit-spot magma aligns with the previously suggested petrogenesis of alkaline rocks explained by the addition of CO2-rich components and/or recycled crustal materials with or without sediment to the mantle (e.g., Dasgupta et al., 2007; Hofmann, 1997). Conversely, the melting of an amphibole-rich metasomatic vein explains the major and trace element composition of alkali basalts (Pilet et al., 2008; Pilet, 2015). However, the experimentally produced melts exhibit Pb depletion and a positive NbTi anomaly in the PM-normalized trace element patterns (Fig. S8), which is inconsistent with the petit-spot basalts (Fig. 7). Moreover, Juriček and Keppler (2023) demonstrated that amphibole dehydration is not the cause for the oceanic LAB through high-pressure experiments under the realistic conditions. The fertile pyroxenitic xenoliths and pyroxene xenocrysts in the 1466R6 and R7 basalts, originating from the metasomatic vein related to prior petit-spot magmatism, had neither amphiboles nor other hydrous minerals (Mikuni et al., 2022).

To explore the involvement of carbonatitic and crustal components in petit-spot melts, a partial melting model of the heterogeneous mantle is presented. The involvement of carbonatitic fluids and recycled materials in the genesis of petit-spot melts has been suggested, and the open-system model with carbonatite influx from the outer system was employed, using OSM-4 by Ozawa (2001) and by referring the parameters of Borisova and Tilhac (2021). This model is based on the mass conservation equations of one-dimensional steady-state melting. In this study, the model asset the critical melt fraction (αc; mass fraction of melt when melt separation begins, leading to the melt connectivity threshold) at 0.005 or 0.01. The system opens to fluxing at a constant melt separation rate (γ) when the system reaches the αc. The final trapped melt fraction (αf; mass fraction of melt trapped in the residue) was fixed at  0 (it was calculated as 10−6, owing to mass balance). We calculated the trace element composition of partial melts at various degrees of melting (F), as well as a few rates of influx (β) and melt separation (γ). We assumed a primitive mantle (PM) source as the lherzolite with or without a normal (N)-MORB source as the recycled oceanic crust (Sun and McDonough, 1989), such as pyroxenite and eclogite. The recycled crust (N-MORB component) was mixed in the source as compositional heterogeneity and calculated as 0.05N-MORB + 0.95 PM for trace element concentration. The mineral phases and their proportions considered were derived only from garnet lherzolite (i.e., olivine, orthopyroxene, clinopyroxene, and garnet). The mineral mode of garnet lherzolite (olivine 55 %, orthopyroxene 20 %, clinopyroxene 15 %, and garnet 10 %) and the melting reaction mode (olivine 8 %, orthopyroxene 19 %, clinopyroxene 81 %, and garnet 30 %) are based on studies by Johnson et al. (1990) and Walter (1998), respectively. The proportion of olivine and garnet was also changed to assess the effect of the garnet modal ratio on the produced melt composition. In this situation, the clinopyroxene is consumed at a degree of partial melting of  19 %; hence, the system was calculated up to 18 % partial melting. The carbonatite melt used in this model as a influx is the average carbonatite from a study by Bizimis et al. (2003). The partition coefficient of trace elements is generally based on a study by McKenzie and O'Nions (1991, 1995), excluding Ti for clinopyroxene and garnet (Kelemen et al., 2003). The variables of β (influx rate) and γ (melt separation rate) were changed during the modeling within the mass balance (γβ+1). The modeled melts were outputted as total melt, considering the instantaneous and accumulated melts. For the carbonatite composition, the value of the average carbonatite from Bizimis et al. (2003) is applied because the chemical composition of carbonatite is largely diverse, and this value is recommended for geochemical modeling (Bizimis et al., 2003). The parameters are detained in Table S6. Consequently, partial melting of garnet lherzolite with a 10 % carbonatite influx to a given mass of source (i.e., garnet lherzolite) can provide a rough explanation of the trace element pattern of petit-spot basalts (Fig. 11a–e). The most plausible for petit-spot magma generation involves the presence of a 5 % crustal component in the source (Fig. 11b and d). In addition, having slightly less garnet in the lherzolite source than the modal ratio of Johnson et al. (1990) offers a better fit for petit-spot characteristics (Fig. 11b). In both scenarios, incorporating a crustal component in the source produces more plausible outcomes (Fig. 11a–d). The higher carbonatite influx (β=1.0) could not explain the trace element composition of the petit-spot basalts (Fig. 11f). A melt connectivity threshold (αc) of 0.01 is considered plausible, as a higher connectivity of melt (i.e., lower αc value) leads to enrichment of large-ion lithophile elements (LILEs) and LREEs (Fig. 11g). The results also indicate that the melt separation ratio has no significant impact on the trace element composition of the calculated melts (Fig. 11d and e). Thereafter, we concluded that the partial melting of  5 % crustal component bearing garnet lherzolite with  10 % carbonatite flux to a given mass of the source plausibly explains the melting source of petit-spot volcanoes (Fig. 11b and d). Assuming that the trace element composition of 1203, 1206, 1542, and 1544 basalts is also primitive, they could be explained by the partial melting of garnet lherzolite with 5 % crustal component and lower carbonatite influx rate (β=0.03) (Fig. S9). Actually, the 1203, 1206, 1542, and 1544 basalts exhibited similar MgO contents and Mg# to those of the 1522 basalts (Fig. 4 and Table 2). These results provide quantitative evidence regarding the petrogenesis of petit-spots, i.e., the contribution of carbonatite melt and recycled oceanic crust.

Figure 11Geochemical modeling for the PM-normalized trace element pattern. The calculated hypothetical melts are a production of carbonatite influx melting of garnet lherzolite with or without 5 % crustal component. Detailed information of the parameters is described in Sect. 6.3 and Table S6. F is the degree of melting (%). The trace element composition of the western Pacific petit-spot basalts from the 6K#1522 dive is shown as black lines for comparison. The PM composition of lherzolite and the N-MORB composition of the recycled crust were based on a study by Sun and McDonough (1989). The influx carbonatite is the average carbonatite from a study by Bizimis et al. (2003). The parameters used in the open-system melting models were as follows: ac is a critical melt fraction, af is a final trapped melt fraction, β is a melt influx rate, and γ is a melt separation rate. Model results are compared by varying each parameter, i.e., garnet modal ratio and presence of crustal material (a–d), melt separation rate (d, e), carbonatite influx rate (d, f), and critical melt fraction (d, g). Each figure is expressed based on the difference from the condition in panel (d).


Figure 12Schematic illustration of the magmatic processes of the western Pacific petit-spot volcanoes. Carbonatitic melt and recycled oceanic crust potentially induce the partial melting of the asthenospheric mantle beneath the western Pacific region. Carbonatitic melt might have originated from a carbon-rich component that horizontally migrated from a subduction zone (Hammouda et al., 2020) or a redox melting of reduced carbon in the deep mantle (Chen et al., 2022; Grassi and Schmidt, 2011; Rohrbach et al., 2007). Petit-spot magma stagnated in the lithosphere with fractional crystallization and melt–rock interaction (Mikuni et al., 2022), and they have erupted at  0–3 Ma.


Although the melting source included small proportions of carbonatite melt and crustal components, these components could have contributed to the isotopic composition, owing to their abundant incompatible elements, as opposed to the ambient mantle. Determination of the Sr, Nd, and Pb isotopic compositions indicated that they had geochemically identical prevalent mantle (PREMA)-like sources (Fig. 9). Contrary to those of the NW Pacific petit-spots, which exhibit EM-1 isotopic composition (Machida et al., 2009; Liu et al., 2020), the samples herein did not align with any mantle isotopic endmembers (i.e., depleted MORB mantle (DMM); EM-1 and EM-2; and high-μ (HIMU); Fig. 9). In the Pb isotopic space, the present samples did not correlate with those of the neighboring HIMU-like Cretaceous seamounts (Fig. 9a) (N Wake and S Wake seamounts; Konter et al., 2008; Koppers et al., 2003; Natland, 1976; Smith et al., 1989; Staudigel et al., 1991). For the melting source of the NW Pacific petit-spot basalts, the involvement of the eclogite/pyroxenite endmember as recycled oceanic crust and the carbonated endmember was suggested. This suggestion was based on the major and trace elements and the Mg, Sr, Nd, and Pb isotopic compositions with Mg diffusion modeling (Liu et al., 2020). The higher FeO / MnO ratios observed in the present melts (65.9–78.0), compared to those of partial melts originating from peridotite (50–60), are attributed to the presence of recycled pyroxenite (Herzberg, 2011), potentially contributing to crustal components in the melting source. However, the western Pacific petit-spots in this study uniformly displayed a PREMA-like isotopic signature without extreme endmember contributions, as described previously (Fig. 9). Such isotopic compositions with the world's petit-spots can be possibly explained by the diverse mixing proportion of HIMU and EM-1 components (Fig. 9e). The isotopic compositions of the NW Pacific petit-spots (off the Japan Trench), Samoan petit-spots (off the Tonga Trench), petit-spot dikes in Christmas Island (off the Java Trench), and western Pacific petit-spots (off the Mariana Trench in this study) are roughly along the HIMU–EM-1 mixing line (Fig. 9e). Furthermore, the isotopic compositions of global carbonatites can generally be explained by the mixing of HIMU and EM-1 (Bell and Tilton, 2002; Hoernle et al., 2002; Hulett et al., 2016). The contributions of the carbonated material/carbonatite and crustal components to the melting source were suggested in relation to the origin of HIMU and EM-1 (Collerson et al., 2010; Hanyu et al., 2011; Wang et al., 2018; Weiss et al., 2016; Workman et al., 2004; Zindler and Hart, 1986). However, the determination of EM-1 and HIMU components as carbonated components and recycled crust, respectively, is challenging due to the varied perspectives on each tectonic setting for the mantle endmember. The variability in the global carbonatite isotopic compositions poses challenges in determining their representative isotope ratios (Fig. 9). Despite these challenges hindering a quantitative isotopic mixing model, the HIMU-EM-1-like trend observed in global petit-spot volcanoes suggests the involvement of carbonatitic and recycled crustal materials. In conclusion, the mass balance models applied to trace elements and the isotopic variations in the petit-spot volcanoes confirmed the contribution of carbonatite melt and the recycled oceanic crust to the melting source of the western Pacific petit-spots (Fig. 12). Experimental studies have revealed the diverse petrogenesis scenarios of carbonatite and carbonatitic alkali-rich magma under high pressures (Dasgupta et al., 2006; Ghosh et al., 2009). The geochemistry of petit-spot basalts, including Mg isotopes, suggested that the conceivable origin of carbonatite related to the petit-spot melt is subducted carbonated pelite, pyroxenite and/or eclogite, or peridotite stored as diamond or metal carbide in the reduced lower portion of the upper mantle (Liu et al., 2020; Rohrbach et al., 2007). For instance, subducted carbonated pelite would melt under high pressure (> 8 GPa) through oxidation at the redox boundary where the iron–wüstite (IW) buffer changes to the quartz–fayalite–magnetite (QFM) buffer (i.e., redox melting; Grassi and Schmidt, 2011). Chen et al. (2022) demonstrated that the alkali-rich carbonatite melt could occur at a pressure exceeding 6 GPa, particularly exhibiting K-rich and Na-rich carbonatites under 6–12 and > 12 GPa, respectively. This pressure-dependent alkalinity of the resulting carbonatite melts could potentially account for the differences between potassic NW Pacific petit-spot lavas and present sodic petit-spot lavas (Fig. 4b). On the other hand, an experimental study highlighted the presence of a carbonate-rich layer in the LAB owing to the horizontally spread carbonate from around the wedge mantle rather than upwelling from the deep mantle (Hammouda et al., 2020). Several high-pressure–temperature experiments and modeling revealed that the chemical composition of intraplate magmas originating from the upper mantle depends on their original depth. Specifically, the carbonatitic melt can be generated beneath thick cratonic lithosphere ( 250–200 km), kimberlitic melt could be produced at > 120 km in depth, and alkali basalt could occur at 100–60 km depth by the partial melting of the original CO2- and H2O-bearing mantle (Massuyeau et al., 2021). This depth-dependent variation in composition, i.e., K-rich kimberlite to alkali basalt, may provide an explanation for the geochemical gap between K-rich NW Pacific petit-spots and K-poor western Pacific petit-spots (Fig. 4b). Although the multiple origins of carbonatite are merely suggested and remain unclear, carbon-rich components play a key role in the partial melting of mantle at the LAB (Sifré et al., 2014), constituting the source of petit-spot magma.

7 Conclusion

The occurrence of petit-spot volcanism supports partial melting at the LAB, carrying significant implications for the characteristics of this geophysical discontinuity. Numerous instances of petit-spot magmatism occurred on the western Pacific Plate at  0–3 Ma, originating from similar PREMA-like melting sources, based on 40Ar/39Ar dating and the Sr, Nd, and Pb isotopic compositions. The mass-balance-based open-system modeling for trace elements revealed that the western Pacific petit-spot magma was generated by the partial melting of a small amount (5 %) of oceanic crust bearing garnet lherzolite with 3 %–10 % carbonatite influx to a given mass of the source. The isotopic compositions of Sr, Nd, and Pb of the study samples, in conjunction with those of the NW Pacific petit-spots, i.e., petit-spots off the Tonga and Java trenches, could be explained by mixing the EM-1-like and HIMU-like components, contributing to the subducted carbonated/crustal materials. The tectonic-induced magmatism, such as a petit-spot, may follow a similar melting mechanism.

Data availability

The data newly analyzed in this study and results of geochemical modeling are included in digital format in the online data repository (Mikuni et al., 2024) of this paper (see Tables 1 to 4 and S1 to S6).


The supplement related to this article is available online at:

Author contributions

KM and NH conceived the project and performed all experiments. SM and YK contributed the Sr, Nd, and Pb isotopic analysis using TIMS and MC-ICP-MS. HS contributed the 40Ar/39Ar dating. NA, AT, and TM helped and performed EPMA and LA-ICP-MS analyses. SM and NH conducted the research cruises to gain the rock samples. All authors interpreted the data and wrote the article, providing comments and improvements.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.


We would like to thank the captains, crews, and shipboard scientific parties of the R/V Yokosuka and the operating team of the submersible Shinkai 6500 for their great work during the YK16-01, YK18-08, and YK19-05S cruises. We used the submersible photos, rock samples, and survey information for these cruises provided in the Data and Sample Research System for Whole Cruise Information by JAMSTEC (, last access: 18 December 2023). The Kyoto University Research Reactor Institute is gratefully acknowledged for their assistance of undertaking the radiometric dating. We would like to express our great appreciation to Tatsuki Tsujimori (ORCID: 0000-0001-9202-7312) for his effort in the management of the laboratory at Tohoku University. We also thank Ryo Fukushima (ORCID: 0000-0003-2683-6757) for improving the wording in the paper. We are really grateful to Yuki Matamura, Yosuke Shimbo, and Yuya Jindo for their help and discussion on scientific matters. The authors would like to thank Enago (, last access: 22 December 2023) for the English language review. This paper was reviewed and improved by two anonymous reviewers, topic editor Massimo Coltorti, executive editor Andrea Di Muro, and the editorial board of Solid Earth. This research has been supported by the Cooperative Program (grant nos. 106 and 202) of Atmosphere and Ocean Research Institute, the University of Tokyo.

Financial support

This research has been supported by the Japan Society for the Promotion of Science (grant nos. 17K05715, 18H03733, and 20K04098).

Review statement

This paper was edited by Massimo Coltorti and reviewed by two anonymous referees.


Aftabuzzaman, M. R., Yomogoda, K., Suzuki, S., Takayanagi, H., Ishigaki, A., Machida, S., Asahara, Y., Yamamoto, K., Hirano, N., Sano, S.-I., Chiyonobu, S., Bassi, D., and Iryu, Y.: Multi-approach characterization of shallow-water carbonates off Minamitorishima and their depositional settings/history, Isl. Arc, 30, e12400,, 2021. 

Akizawa, N., Ohara, Y., Okino, K., Ishizuka, O., Yamashita, H., Machida, S., Sanfilippo, A., Basch, V., Snow, J. E., Sen, A., Hirauchi, K.-I., Michibayashi, K., Harigane, Y., Fujii, M., Asanuma, H., and Hirata, T.: Geochemical characteristics of back-arc basin lower crust and upper mantle at final spreading stage of Shikoku Basin: an example of Mado Megamullion, Prog. Earth Planet. Sci., 8, 65,, 2021. 

Akizawa, N., Hirano, N., Matsuzaki, K. M., Machida, S., Tamura, C., Kaneko, J., Iwano, H., Danhara, T., and Hirata, T.: A direct evidence for disturbance of whole sediment layer in the subducting Pacific plate by petit-spot magma–water/sediment interaction, Mar. Geol., 444, 106712,, 2022. 

Asimow, P. D. and Langmuir, C. H.: The importance of water to oceanic mantle melting regimes, Nature, 421, 815–820,, 2003. 

Audhkhasi, P. and Singh, S. C.: Discovery of distinct lithosphere-asthenosphere boundary and the Gutenberg discontinuity in the Atlantic Ocean, Sci. Adv., 8, eabn5404,, 2022. 

Axen G. J., van Wijk, J. W., and Currie, C. A.: Basal continental mantle lithosphere displaced by flat-slab subduction, Nat. Geosci., 11, 961–964,, 2018. 

Azami, K., Machida, S., Hirano, N., Nakamura, K., Yasukawa, K., Kogiso, T., Nakanishi, M., and Kato, Y.: Hydrothermal ferromanganese oxides around a petit-spot volcano on old and cold oceanic crust, Commun. Earth Environ., 4, 191,, 2023. 

Bell, K. and Tilton, G. R.: Probing the mantle: the story from carbonatites, Eos, 83, 273–277,, 2002.  

Bellas, A., Zhong, S., and Watts, A. B.: Reconciling lithospheric rheology between laboratory experiments, field observations and different tectonic settings, Geophys. J. Int., 228, 857–875,, 2022. 

Bianco, T. A., Ito, G., Becker, J. M., and Garcia, M. O.: Secondary Hawaiian volcanism formed by flexural arch decompression, Geochem. Geophy. Geosy., 6, Q08009,, 2005. 

Bizimis, M., Salters, V. J. M., and Dawson, J. B.: The brevity of carbonatite sources in the mantle: evidence from Hf isotopes, Miner. Petrol., 145, 281–300,, 2003. 

Bizimis, M., Salters, V. J. M., Garcia, M. O., and Norman, M. D.: The composition and distribution of the rejuvenated component across the Hawaiian plume: Hf-Nd-Sr-Pb isotope systematics of Kaula lavas and pyroxenite xenoliths, Geochem. Geophy. Geosy., 14, 4458–4478,, 2013. 

Borisova, A. Y. and Tilhac, R.: Derivation of Hawaiian rejuvenated magmas from deep carbonated mantle sources: A review of experimental and natural constraints, Earth. Sci. Rev., 222, 103819,, 2021. 

Buchs, D. M., Pilet, S., Cosca, M., Flores, K. E., Bandini, A. N., and Baumgartner, P. O.: Low-volume intraplate volcanism in the Early/Middle Jurassic Pacific basin documented by accreted sequences in Costa Rica, Geochem. Geophy. Geosy.,14, 1552–1568,, 2013. 

Chantel, J., Manthilake, G., Andrault, D., Novella, D., yu, T., and Wang, Y.: Experimental evidence supports mantle partial melting in the asthenosphere, Sci. Adv., 2, e1600246,, 2016. 

Chen, X., Wang, M., Inoue, T., Liu, Q., Zhang, L., and Bader, T.: Melting of carbonated pelite at 5.5–15.5 GPa: implications for the origin of alkali-rich carbonatites and the deep water and carbon cycles, Miner. Petrol., 177, 2,, 2022. 

Clague, D. A. and Frey, F. A.: Petrology and Trace element Geochemistry of the Honolulu Volcanics, Oahu: Implications for the Oceanic Mantle below Hawaii, J. Petrol., 23, 447–504,, 1982. 

Clague, D. A. and Moore, J. G.: The proximal part of the giant submarine Wailau landslide, Molokai, Hawaii, J. Volcanol. Geoth. Res., 113, 259–287,, 2002. 

Clague, D. A., Holcomb, R. T., Sinton, J. M., Detrick, R. S., and Torresan, M. E.: Pliocene and Pleistocene alkali flood basalts on the seafloor north of the Hawaiian island, Earth Planet. Sc. Lett., 98, 175–191,, 1990. 

Collerson, K. D., Williams, Q., Ewart, A. E., and Murphy, D. T.: Origin of HIMU and EM-1 domains sampled by ocean island basalts, kimberlites and carbonatites: The role of CO2-fluxed lower mantle melting in thermochemical upwellings, Phys. Earth Planet. In., 181, 112–131,, 2010. 

Conrad, C. P., Bianco, T. A., Smith, E. I., and Wessel, P.: Patterns of intraplate volcanism controlled by asthenospheric shear, Nat. Geosci., 4, 317–321,, 2011. 

Dasgupta, R. and Hirschmann, M. M.: Melting in the Earth's deep upper mantle caused by carbon dioxide, Nature, 440, 659–662,, 2006. 

Dasgupta, R., Hirschmann, M. M., and Stalker, K.: Immiscible Transition from Carbonate-rich to Silicate-rich Melts in the 3 GPa Melting Interval of Eclogite + CO2 and Genesis of Silica-undersaturated Ocean Island Lavas, J. Petrol., 47, 647–671,, 2006. 

Dasgupta, R., Hirschmann, M. M., and Smith, N. D.: Partial Melting Experiments of Peridotite + CO2 at 3 GPa and Genesis of Alkalic Ocean Island Basalts, J. Petrol., 48, 2093–2124,, 2007. 

Dasgupta, R., Hirschmann, M. M., McDonough, W. F., Spiegelman, M., and Withers, A.: Trace element partitioning between garnet lherzolite and carbonatite at 6.6 and 8.6 GPa with applications to the geochemistry of the mantle and of mantle-derived melts, Chem. Geol., 262, 57–77,, 2009. 

Dasgupta, R., Mallik, A., Tsuno, K., Withers, A. C., Hirth, G., and Hirschmann, M. M.: Carbon-dioxide-rich silicate melt in the Earth's upper mantle, Nature, 493, 211–215,, 2013. 

Debayle, E., Bodin, T., Durand, S., and Ricard, Y.: Seismic evidence for partial melt below tectonic plates, Nature, 586, 555–559,, 2020. 

Dixon, J., Clague, D. A., Cousens, B., Monsalve, M. L., and Uhl, J.: Carbonatite and silicate melt metasomatism of the mantle surrounding the Hawaiian plume: evidence from volatiles, trace elements, and radiogenic isotopes in rejuvenated-stage lavas from Niihau, Hawaii, Geochem. Geophy. Geosy., 9, Q09005,, 2008. 

Ebisawa, N., Sumino, H., Okazaki, R., Takigami, Y., Hirano, N., Nagao, K., and Kaneoka, I.: Construction of I-Xe and 40Ar39Ar dating system using a modified VG3600 noble gas mass spectrometer and the first I-Xe data obtained in Japan, J. Mass Spectrom. Soc. Jpn., 52, 219–229,, 2004. 

Falloon, T. J. and Green, D. H.: The solidus of carbonated, fertile peridotite, Earth Planet. Sc. Lett., 94, 364–370,, 1989. 

Falloon, T. J. and Green, D. H.: Solidus of carbonated fertile peridotite under fluid-saturated conditions, Geology, 18, 195–199,<0195:SOCFPU>2.3.CO;2, 1990. 

Falloon, T. J. Hoernle, K., Schaefer, B. F., Bindeman, I. N., Hart, S. R., Garbe-Schonberg, D., and Duncan, R. A.: Petrogenesis of Lava from Christmas Island, Northeast Indian Ocean: Implications for the Nature of Recycled Components in Non-Plume Intraplate Settings, Geosciences, 12, 118,, 2022. 

Foley, S. F., Yaxley, G. M., Rosenthal, A., Buhre, S., Kiseeva, E. S., Rapp, R. P., and Jacob, D. E.: The composition of near-solidus melts of peridotite in the presence of CO2 and H2O between 40 and 60 kbar, Lithos, 112, 274–283,, 2009. 

Frey, F. A., Green, D. H., and Roy, S. D.: Integrated Models of Basalt Petrogenesis: A Study of Quartz Tholeiites to Olivine Melilitites from South Eastern Australia Utilizing Geochemical and Experimental Petrological Data, J. Petrol., 19, 463–513,, 1978. 

Frey, F. A., Clague, D., Mahoney, J. J., and Sinton, J. M.: Volcanism at the edge of the Hawaiian plume: Petrogenesis of submarine alkali lavas from the North Arch volcanic field, J. Petrol., 41, 667–691,, 2000. 

Fujie, G., Kodaira, S., Nakamura, Y., Morgan, J. P. Dannowski, A., Thorwart, M., Grevemeyer, I., and Miura, S.: Spatial variations of incoming sediments at the northeastern Japan arc and their implications for megathrust earthquakes, Geology, 48, 614–619,, 2020. 

Fujiwara, T., Hirano, N. Abe, N., and Takizawa, K.: Subsurface structure of the “petit-spot” volcanoes on the northwestern Pacific Plate, Geophys. Res. Lett., 34, L13305,, 2007. 

Garcia, M. O., Weis, D., Jicha, B. R., Ito, G., and Hanano, D.: Petrology and geochronology of lavas from Ka`ula Volcano: Implications for rejuvenated volcanism of the Hawaiian mantle plume, Geochim. Cosmochim. Ac., 185, 278–301,, 2016. 

Ghosh, S., Ohtani, E., Litasov, K. K., and Terasaki, H.: Solidus of carbonated peridotite from 10 to 20 GPa and origin of magnesiocarbonatite melt in the Earth's deep mantle, Chem. Geol., 262, 17–28,, 2009. 

Grassi, D. and Schmidt, M. W.: The Melting of Carbonated Pelites from 70 to 700 km Depth, J. Petrol., 52, 765–789,, 2011. 

Gripp, A. E. and Gordon, R. G.: Current plate velocities relative to the hotspots incorporating the NUVEL-1 global plate motion model, Geophys. Res. Lett., 17, 1109–1112,, 1990. 

Hammouda, T., Manthilake, G., Goncalves, P., Chantel, J., Guignard, J., Crichton, W., and Gaillard, F.: Is There a Global Carbonate Layer in the Oceanic Mantle?, Geophys. Res. Lett., 48, e2020GL089752,, 2020. 

Hanano, D., Scoates, J. S., and Weis, D.: Alteration mineralogy and the effect of acid-leaching on the Pb-isotope systematics of ocean-island basalts, Am. Mineral., 94, 17–26,, 2009. 

Hanyu, T., Tatsumi, Y., Senda, R., Miyazaki, T., Chang, Q., Hirahara, Y., Takahashi, T., Kawabata, H., Suzuki, K., Kimura, J-I., and Nakai, S.: Geochemical characteristics and origin of the HIMU reservoir: A possible mantle plume source in the lower mantle, Geochem. Geophy. Geosy., 12, Q0AC09,, 2011. 

Hart, S. R.: A large-scale isotope anomaly in the Southern Hemisphere mantle, Nature, 309, 753–757,, 1984. 

Hart, S. R., Gerlach, D. C., and White, W. M.: A Possible new Sr-Nd-Pb mantle array and consequences for mantle mixing, Geochim. Cosmochim. Ac., 50, 1551–1557,, 1986. 

Hart, S. R., Coetzee, M., Workman, R. K., Blusztajn, L., Johnson, K. T. M., Sinton, J. M., Steinberger, B., and Hawkins, J. W.: Genesis of the Western Samoa seamount province: Age, geochemical fingerprint and tectonics, Earth Planet. Sc. Lett., 227, 37–56,, 2004. 

Hein, J. R., Koschinsky, A., Bau, M., Manheim, F. T., Kang, J. K., and Roberts, L.: Cobalt-rich ferromanganese crusts in the Pacific, in: Handbook of Marine Mineral Deposits, edited by: Cronan D. S., CRC Press, Boca Raton, Florida, 239–279, (last access: 6 February 2024), 1999. 

Helz, R. T. and Thornber, C. R.: Geochemistry of Kilauea Iki lava lake, Hawaii, B. Volcanol., 49, 651–658,, 1987. 

Herath, P., Stern, T. A., Savage, M. K., Bassett, D., and Henrys, S.: Wide-angle seismic reflections reveal a lithosphere-asthenosphere boundary zone in the subducting Pacific Plate, New Zealand, Sci. Adv., 8, eabn5697,, 2022. 

Herzberg, C.: Petrology and thermal structure of the Hawaiian plume from Mauna Kea volcano, Nature, 444, 605–609,, 2006. 

Herzberg, C.: Identification of Source Lithology in the Hawaiian and Canary Islands: Implications for Origins, J. Petrol., 52, 113–146,, 2011. 

Hirano, N.: Petit-spot volcanism: a new type of volcanic zone discovered near a trench, Geochem. J., 45, 157–167,, 2011. 

Hirano, N. and Machida, S.: The mantle structure below petit-spot volcanoes, Commun. Earth Environ., 3, 110,, 2022. 

Hirano, N., Takahashi, E., Yamamoto, J., Abe, N., Ingle, S. P., Kaneoka, I., Hirata, T., Kimura, J.-I., Ishii, T., Ogawa, Y., Machida, S., and Suyehiro, K.: Volcanism in response to plate flexure, Science, 313, 1426–1428., 2006. 

Hirano, N., Machida, S., Abe, N., Morishita, T., Tamura, A., and Arai, S.: Petit-spot lava fields off the central Chile trench induced by plate flexure, Geochem. J., 47, 249–257,, 2013. 

Hirano, N., Nakanishi, M., Abe, N., and Machida, S.: Submarine lava fields in French Polynesia, Mar. Geol., 373, 39–48,, 2016. 

Hirano, N., Machida, S., Sumino, H., Shimizu, K., Tamura, A., Morishita, T., Iwano, H., Sakata, S., Ishii, T., Arai, S., Yoneda, S., Danhara, T., and Hirata, T.: Petit-spot volcanoes on the oldest portion of the Pacific Plate, Deep-Sea Res. Pt. I, 154, 103142,, 2019. 

Hirano, N., Sumino, H., Morishita, T., Machida, S., Kawano, T., Yasukawa, K., Hirata, T., Kato, Y., and Ishii, T.: A Paleogene magmatic overprint on Cretaceous seamounts of the western Pacific, Isl. Arc, 30, e12386,, 2021. 

Hirth, G. and Kohlstedt, D. L.: Water in the oceanic upper mantle: implications for rheology, melt extraction and the evolution of the lithosphere, Earth Planet. Sc. Lett., 144, 93–108,, 1996. 

Hoernle, K., Tilton, G., Le Bas, M. J., Duggem, S., and Garbe-Schönberg, D.: Geochemistry of oceanic carbonatites compared with continental carbonatites: mantle recycling of oceanic crustal carbonate, Miner. Petrol., 142, 520–542,, 2002. 

Hofmann, A. W.: Mantle geochemistry: the message from oceanic volcanism, Nature, 385, 219–229,, 1997. 

Hofmann, A. W.: Sampling mantle heterogeneity through oceanic basalts: isotopes and trace elements, in: Treatise on Geochemistry, 2, The Mantle and Core, edited by: Carson, R. W., Elsevier, 61–101,, 2003. 

Hosseini, K., Matthews, K. J., Sigloch, K., Shephard, G. E., Domeier, M., and Tsekhmistrenko, M.: SubMachine: Web-Based tools for exploring seismic tomography and other models of Earth's deep interior, Geochem. Geophy. Geosy., 19, 1464–1483,, 2018. 

Hua, J., Fisher, K. M., Becker, T. W., Gazel, E., and Hirth, G.: Asthenospheric low-velocity zone consistent with globally prevalent partial melting, Nat. Geosci., 16, 175–181,, 2023. 

Hulett, S. R., Simonetti, A., Rasbury, E. T., and Hemming, N. G.: Recycling of subducted crustal components into carbonatite melts revealed by boron isotopes, Nat. Geosci., 9, 904–908,, 2016. 

Irvine, T. N. and Baragar, W. R. A.: A Guide to the Chemical Classification of the Common Volcanic Rocks, Can. J. Earth Sci., 8, 523–548,, 1971. 

Irving, A. J. and Green, D. H.: Geochemistry and petrogenesis of the newer basalts of Victoria and South Australia, J. Geol. Soc. Aust., 23, 45–66,, 1976. 

Iwata, N.: Geochronological study of the Deccan volcanism by the 40Ar39Ar method, Doctor Thesis, University of Tokyo, 168 pp.,, 1998. 

Jochum, K. P. and Nohl, U.: Reference materials in geochemistry and environmental research and the GeoReM database, Chem. Geol., 253, 50–53,, 2008. 

Johnson, K. T. M., Dick, H. J. B., and Shimizu, N.: Melting in the oceanic upper mantle: An ion microprobe study of diopsides in abyssal peridotites, J. Geophys. Res., 95, 2661–2678,, 1990. 

Juriček, M. P and Keppler, H.: Amphibole stability, water storage in the mantle, and the nature of the lithosphere–asthenosphere boundary, Earth Planet. Sc. Lett., 608, 118082,, 2023. 

Kaneko, J., Machida, S., Hirano, N., Kasaya, T., and Kumagai, H.: Near bottom MBES survey mounted on a HOV at 5500 m depth, in: Oceans Conference Record (IEEE) 2022, 1–24 February 2022, Chennai, India, 1–5,, 2022. 

Kang, L. and Karato, S. -I.: Hydrogen Partitioning Between Olivine and Orthopyroxene: Implications for the Lithosphere-Asthenosphere Structure, J. Geophys. Res., 128, e2022JB025259,, 2023. 

Karato, S.-I. and Jung, H.: Water, partial melting and the origin of the seismic low velocity and high attenuation zone in the upper mantle, Earth Planet. Sc. Lett., 157, 193–207,, 1998. 

Katsura, T. and Fei, H.: Asthenosphere dynamics based on the H2O dependence of element diffusivity in olivine, Natl. Sci. Rev., 8, nwaa278., 2021. 

Kawakatsu, H., Kumar, P., Takei, Y., Shinohara, M., Kanazawa, T., Araki, E., and Suyehiro, K.: Seismic Evidence for Sharp Lithosphere–Asthenosphere Boundaries of Oceanic Plates, Science, 324, 499–502,, 2009. 

Kelemen, P. B., Yogodzinskim G. M., and Scholl, D. W.: Along-strike variation in the Aleutian Island Arc: genesis of high Mg# andesite and implications for continental crust, in: Inside the subduction Factory, edited by: Eiler, J., American Geophysical Union, Geophysical Monograph, 138, 223–276,, 2003. 

Keshav, S. and Gudfinnsson, G. H.: Silicate liquid-carbonatite liquid transition along the melting curve of model, vapor-saturated peridotite in the system CaO-MgO-Al2O3-SiO2-CO2 from 1.1 to 2 GPa, J. Geophys. Res., 118, 3341–3353,, 2013. 

Kiseeva, E. S., Litasov, K. D., Yaxley, G. M., Ohtani, E., and Kamenetsky, V. S.: Melting and Phase Relations of Carbonated Eclogite at 9–21 GPa and the Petrogenesis of Alkali-Rich Melts in the Deep Mantle, J. Petrol., 54, 1555–1583,, 2013. 

Kobayashi, M., Sumino, H., Saito, T., Nagao, K.: Determination of halogens in geological reference materials using neutron irradiation noble gas mass spectrometry, Chem. Geol., 582, 120420,, 2021. 

Konovalov, Y. I. and Martynov, Y. A.: Volcanic complex of the La Mont Guyot; Marcus-Wake Uplift, Pacific Ocean, Pacific Geology, 5, 40–47, 1992. 

Konter, J. G. and Jackson, M. G.: Large volumes of rejuvenated volcanism in Samoa: Evidence supporting a tectonic influence on late-stage volcanism, Geochem. Geophy. Geosy., 13, Q0AM04,, 2012. 

Konter, J. G., Hanan, B. B., Blicher-Toft, J., Koppers, A. A. P., Plank, T., and Staudigel, H.: One hundred million years of mantle geochemical history suggest the retiring of mantle plumes is premature, Earth Planet. Sc. Lett., 275, 285–295,, 2008. 

Koppers, A. A. P., Staudigel, H., and Wijbrans, J. R.: Dating crystalline groundmass separates of altered Cretaceous seamount basalts by the Ar40/Ar39 incremental heating technique, Chem. Geol., 166, 139–158., 2000. 

Koppers, A. A. P., Staudigel, H., Pringle, M. S., and Wijbrans, J. R.: Short-lived and discontinuous intra-plate volcanism in the South Pacific: hotspots or extensional volcanism?, Geochem. Geophy. Geosy., 4, 1089,, 2003. 

Koppers, A. A. P., Russell, J. A., Jackson, M. G., Konter, J., Staudigel, H., and Hart, S. R.: Samoa reinstated as a primary hotspot trail, Geology, 36, 435–438,, 2008. 

Korenaga, J.: Plate tectonics and surface environment: Role of the oceanic upper mantle, Earth Sci. Rev., 205, 103185,, 2020. 

Le Bas, M. J., Le Maitre, R., Strackeisen, A., and Zanettin, B.: A chemical classification of volcanic rocks based on the total alkali–silica diagram, J. Petrol., 27, 745–750,, 1986. 

Liu, J., Hirano, N., Machida, S., Xia, Q., Tao, C., Liao, S., Liang, J., Li W., Yang, W. Zhang, G., and Ding, T.: Melting of recycled ancient crust responsible for the Gutenberg discontinuity, Nat. Commun., 11, 172,, 2020. 

Longerich, H. P., Jackson, S. E., and Gunther, D.: Laser ablation inductively coupled plasma mass spectrometric transient signal data acquisition and analyte concentration calculation, J. Anal. Atom. Spectrom., 11, 899–904,, 1996. 

Lu, C., Grand, S. P., Lai, H., and Garnero, E. J.: TX2019slab: A New P and S Tomography Model Incorporating Subducting Slabs, J. Geophys. Res., 124, 11549–11567,, 2019. 

Machida, S., Hirano, N., and Kimura, J.-I.: Evidence for recycled material in Pacific upper mantle unrelated to plumes, Geochim. Cosmochim. Ac., 73, 3028–3037,, 2009. 

Machida, S., Orihashi, Y., Magnani, M., Neo, N., Wilson, S., Tanimizu, M., Yoneda, S., Yasuda, A., and Tamaki, K.: Regional mantle heterogeneity regulates melt production along the Réunion hotspot-influenced Central Indian Ridge, Geochem. J., 48, 433–449,, 2014. 

Machida, S., Hirano, N., Sumino, H., Hirata, T., Yoneda, S., and Kato, Y.: Petit-spot geology reveals melts in upper-most asthenosphere dragged by lithosphere, Earth Planet. Sc. Lett., 426, 267–279,, 2015. 

Machida, S., Kogiso, T., and Hirano, N.: Petit-spot as definitive evidence for partial melting in the asthenosphere caused by CO2, Nat. Commun., 8, 14302,, 2017. 

Massuyeau, M., Gardés, E., Morizet, Y., and Gaillard, F.: A model for the activity of silica along the carbonatite–kimberlite–mellilitite–basanite melt compositional joint, Chem. Geol., 418, 206–216,, 2015. 

Massuyeau, M., Gardés, E., Rogerie, G., Aulbach, S., Tappe, S., Le Trong, E., Sifré, D., and Gaillaer, F.: MAGLAB: A computing platform connecting geophysical signatures to melting processes in Earth's mantle, Phys. Earth Planet. In., 314, 106638,, 2021. 

McKenzie, D. and O'Nions, R. K.: Partial melt distributions from inversion of rare Earth element concentrations, J. Petrol., 32, 1021–1091,, 1991. 

McKenzie, D. and O'Nions, R. K.: The Source Regions of Ocean Island Basalts, J. Petrol., 36, 133–159,, 1995. 

Melson, W. G., Thompson, G., and van Andel, T. H.: Volcanism and metamorphism in the Mid-Atlantic Ridge, 22 N latitude, J. Geophys. Res., 73, 5925–5941,, 1968. 

Mierdel, K., Keppler, H., Smyth, J. R., and Langenhorst, F.: Water solubility in aluminous orthopyroxene and the origin of Earth's Asthenosphere, Science, 315, 364–368,, 2007. 

Mikuni, K., Hirano, N., Akizawa, N., Yamamoto, J., Machida, S., Tamura, A., Hagiwara, Y., and Morishita, T.: Lithological structure of western Pacific lithosphere reconstructed from mantle xenoliths in a petit-spot volcano, Prog. Earth Planet. Sci., 9, 62,, 2022. 

Mikuni, K., Hirano, N., Machida, S., Sumino, H., Akizawa, N., Tamura, A., Morishita, T., Kato, Y.: Results of geochemical analysis and modeling of Mikuni et al. “Contribution of carbonatite and recycled oceanic crust to petit-spot lavas on the western Pacific Plate”, Version 1.0, Interdisciplinary Earth Data Alliance (IEDA) [data set],, 2024. 

Miyashiro, A., Shido, F., and Ewing, M.: Metamorphism on the Mid-Atlantic Ridge near 24 and 30 N, Philos. T. R. Soc. Lond., 268, 589–603,, 1971. 

Moore, J. G., Fornari, D. J., and Clague, D. A.: Basalts from the 1877 Submarine Eruption of Mauna Loa, Hawaii; New Data on the Variation of Palagonitization Rate with Temperature, United States Geol. Surv. Bull. 1663, US Geological Survye, 1–11,, 1985. 

Müller, R. D., Sdrolias, M., Gaina, C., and Roest, W. R.: Age, spreading rates, and spreading asymmetry of the world's ocean crust, Geochem. Geophy. Geosy., 9, Q04006,, 2008. 

Natland, J.: Petrology of Volcanic Rocks Dredged from Seamounts in the Line Islands, Initial Rep. Deep Sea, 33, 749–777., 1976. 

Nier, A.: A redetermination of the relative abundances of the isotopes of carbon, nitrogen, oxygen, argon, and potassium, Phys. Rev., 77, 789–793,, 1950. 

Nobre Silva, I. G., Weis, D., Barling, J., and Scoates, J. S.: Leaching systematics and matrix elimination for the determination of high-precision Pb isotope compositions of ocean island basalts, Geochem. Geophy. Geosy., 10, Q08012,, 2009. 

Novella, D., Keshav, S., Gudfinnsson, G. H., and Ghosh, S.: Melting phase relations of model carbonated peridotite from 2 to 3 GPa in the system CaO-MgO-Al2O3-SiO2-CO2 and further indication of possible unmixing between carbonatite and silicate liquids, J. Geophys. Res., 119, 2780–2800,, 2014. 

Nozaki, T., Tokumaru, A., Takaya, Y., Kato, Y., Suzuki, K., and Urabe, T.: Major and trace element compositions and resource potential of ferromanganese crust at Takuyo Daigo Seamount, northwestern Pacific Ocean, Geochem. J., 50, 527–537,, 2016. 

Okumura, S. and Hirano, N.: Carbon dioxide emission to earth's surface by deep-sea volcanism, Geology, 41, 1167–1170,, 2013. 

Orihashi, Y., Maeda, J., Tanaka, R., Zeniya, R., and Niida, K.: Sr and Nd isotopic data for the seven GSJ rock reference samples; JA-1, JB-1a, JB-2, JB-3, JG-1a, JGb-1 and JR-1, Geochem. J., 32, 205–211,, 1998. 

Ozawa, K.: Mass balance equations for open magmatic systems: Trace element behavior and its application to open system melting in the upper mantle. J. Geophys. Res., 106, 13407–13434,, 2001. 

Pearce, N. J. G., Perkins, W. T., Westgate, J. A., Gorton, M. P., Jackson, S. E., Neal, C. R., and Chenery, S. P.: A compilation of new and published major and trace element data for NIST SRM 610 and NIST SRM 612 glass reference materials, Geostanddard. Newslett., 21, 115–144,, 1997. 

Pilet, S.: Generation of low-silica alkaline lavas: Petrological constrains, models, and thermal implications, in: The Interdisciplinary Earth: A Volume in Honor of Don L. Anderson, edited by: Foulger, G. R., Lustrino, M., and King, S. D., The Geological Society of America,, 2015. 

Pilet, S., Baker, M. B., and Stolper, E. M.: Metasomatized Lithosphere and the Origin of Alkaline Lavas, Science, 320, 916–919,, 2008. 

Pilet, S., Abe, N., Rochat, L., Kaczmarek, M.-A., Hirano, N., Machida, S., Buchs, D. M., Baumgarther, P. O., and Müntener, O.: Pre-subduction metasomatic enrichment of the oceanic lithosphere induced by plate flexure, Nat. Geosci., 9, 898–903,, 2016. 

Reinhard, A. A., Jackson, M. G., Blusztajn, J., Koppers, A. A. P., Simms, A. R., and Konter, J. G.: “Petit Spot” Rejuvenated Volcanism Superimposed on Plume-Derived Samoan Shield Volcanoes: Evidence From a 645 m Drill Core From Tutuila Island, American Samoa, Geochem. Geophy. Geosy., 20, 1485–1507,, 2019. 

Resing, J. A. and Sansone, F. J.: The chemistry of lava–seawater interactions: the generation of acidity, Geochim. Cosmochim. Ac., 63, 2183–2198,, 1999. 

Rohrbach, A., Ballhaus, C., Golla-Schindler, U., Ulmer, P., Kamenetsky, V. S., and Kuzmin, D. V.: Metal saturation in the upper mantle, Nature, 449, 456–458,, 2007. 

Rychert, C. A. and Shearer, P. M.: A global view of the lithosphere–asthenosphere boundary, Science, 324, 495–498,, 2009. 

Sakamaki, T., Suzuki, A., Ohtani, E., Terasaki, H., urakawa, S., Katayama, Y., Funakoshi, K.-I., Wang, Y. Hernlund, J. H., and Ballmer, M. D.: Ponded melt at the boundary between the lithosphere and asthenosphere, Nat. Geosci., 6, 1041–1044,, 2013. 

Shaw, C. S. J.: Dissolution of orthopyroxene in basanitic magma between 0.4 and 2 GPa: Further implications for the origin of Si-rich alkaline glass inclusions in mantle xenoliths, Miner. Petrol., 135, 114–132,, 1999. 

Sifré, D., Gardés, E., Massuyeau, M., Hashim, L., Hier-Majumder, S., and Gaillard, F.: Electrical conductivity during incipient melting in the oceanic low-velocity zone, Nature, 509, 81–85,, 2014. 

Smith, W. H. F., Staudigel, H., Watts, A. B., and Pringle, M. S.: The Magellan seamounts: early Cretaceous record of the South Pacific isotopic and thermal anomaly, J. Geophys. Res., 94, 10501–10523,, 1989. 

Staudigel, H. and Hart, S. R.: Alteration of basaltic glass: processes and significance for the oceanic crust-sewater budget, Geochim. Cosmochim. Ac., 47, 337–350,, 1983. 

Staudigel, H., Park, K. H., Pringle, M., Rubenstone, J. L., Smith, W. H. F., and Zindler, A.: The longevity of the South-Pacific isotopic and thermal anomaly, Earth Planet. Sc. Lett., 102, 24–44,, 1991. 

Stixrude, L. and Lithgow-Bertelloni, C.: Thermodynamics of mantle minerals — I. Physical properties, Geophys. J. Int., 162, 610–632,, 2005. 

Stracke, A., Michael, W., Felix, G., Paul, B., and Erin, T.: Major and trace element concentrations and Sr, Nd, Hf, Pb isotope ratios of global mid ocean ridge and ocean island basalts, GRO data, V1, DIGIS [data set],, 2022. 

Sun, S.-S. and McDonough, W. F.: Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes, Geol. Soc. Spec. Publ., 42, 313–345,, 1989. 

Takahashi, E.: Origin of basaltic magmas: Implications from peridotite melting experiments and an olivine fractionation model, B. Volcanol. Soc. Jpn., 2nd Ser, 30, S17–S40,, in Japanese with English abstract, 1986. 

Tamura, A., Arai, S., Takeuchi, M., Miura, M., and Pirnia, T.: Compositional heterogeneity of a websterite xenolith from Kurose, southwest Japan: insights into the evolution of lower crust beneath the Japan Arc, Eur. J. Mineral., 31, 35–47,, 2019. 

Taneja, R., Rushmer, T., Blichert-Toft, J., Turner, S., and O'Neill, C.: Mantle heterogeneities beneath the Northeast Indian Ocean as sampled by intra-plate volcanism at Christmas Island, Lithos, 262, 561–575,, 2016. 

Tanimizu, M. and Ishikawa, T.: Development of rapid and precise Pb isotope analytical techniques using MC-ICPMS and new results for GSJ rock reference samples, Geochem. J., 40, 121–133., 2006. 

Tatsumi, Y., Sakuyama, M., Fukuyama, H., and Kushiro, I.: Generation of arc basalt magmas and thermal structure of the mantle wedge in subduction zones, J. Geophys. Res., 88, 5815–5825,, 1983. 

Tivey, M. A., Sager, W. W., Lee, S.-M., and Tominaga, M.: Origin of the Pacific Jurassic quiet zone, Geology, 34, 789–792,, 2006. 

Uenzelmann-Neben, G., Schmidt, D. N., Niessen, F., and Stein, R.: Intraplate volcanism off South Greenland: caused by glacial rebound?, Geophys. J. Int., 190, 1–7,, 2012. 

Valentine, G. A. and Hirano, N.: Mechanisms of low-flux intraplate volcanic fields—Basin and Range (North America) and northwest Pacific Ocean, Geology, 38, 55–58,, 2010. 

Wakaki, S., Shibata, S.-N., and Tanaka, T.: Isotope ratio measurements of trace Nd by the total evaporation normalization (TEN) method in thermal ionization mass spectrometry, Int. J. Mass Spectrom., 264, 157–163,, 2007. 

Walter, M. J.: Melting of garnet peridotite and the origin of komatiite and depleted lithosphere, J. Petrol., 39, 29–60,, 1998. 

Wang, D., Mookherjee, M., Xu Y., and Karato, S.-I.: The effect of water on the electrical conductivity of olivine, Nature, 443, 977–980,, 2006. 

Wang, X.-J., Chen, L.-H., Hofmann, A. W., Hanyu, T., Kawabata, H., Zhong, Y., Xie, L.-W., Shi, J.-H., Miyazaki, T., Hirata, Y., Takahashi, T., Senda, R., Chang, O., Vaglarov, B. S., and Kimura, J.-I.: Recycled ancient ghost carbonate in the Pitcairn mantle plume, P. Natl. Acad. Sci. USA, 115, 8682–8687,, 2018. 

Weis, D. and Frey, F. A.: Isotope geochemistry of the Ninetyeast Ridge basement basalts: Sr, Nd, and Pb evidence for involvement of the Kerguelen hot spot, Proc. Ocean Drill. Program Sci. Results, 121, 591–610, 1991. 

Weis, D. and Frey, F. A.: Role of the Kerguelen Plume in generating the eastern Indian Ocean seafloor. J. Geophys. Res., 101, 13381–13849,, 1996. 

Weis, D., Kieffer, B., Maerschalk, C., Barling, J., de Jong, J., Williams, G. A., Hanano, D., Pretorius, W., Mattielli, N., Scoates, J. S., Goolaerts, A., Friedman, R. M., and Mahoney, J. B.: High-precision isotopic characterization of USGS reference materials by TIMS and MC-ICP-MS, Geochem. Geophy. Geosy., 7, Q08006,, 2006. 

Weiss, Y., Class, C., Goldstein, S. L., and Hanyu, T.: Key new pieces of the HIMU puzzle from olivines and diamond inclusions, Nature, 537, 666–670,, 2016. 

Wessel, P., Luis, J. F., Uieda, L., Scharroo, R., Wobbe, F., Smith, W. H. F., and Tian, D.: The Generic Mapping Tools version 6, Geochem Geophys Geosyst., 20, 5556–5564,, 2019. 

Workman, R. K., Hart, S. R., Jackson, M., Regelous, M., Farley, K. A., Blusztajn, J., Kurz, M., and Staudigel, H.: Recycled metasomatized lithosphere as the origin of the Enriched Mantle II (EM2) end-member: Evidence from the Samoan Volcanic Chain, Geochem. Geophy. Geosy., 5, Q04008,, 2004. 

Yamamoto, J., Korenaga, J., Hirano, N., and Kagi, H.: Melt-rich lithosphere-asthenosphere boundary inferred from petit-spot volcanoes, Geology, 42, 967–970,, 2014. 

Yamamoto, J., Kawano, T., Takahata, N., and Sano, Y.: Noble gas and carbon isotopic compositions of petit-spot lavas from southeast of Marcus Island, Earth Planet. Sc. Lett., 497, 139–148,, 2018. 

Yamamoto, J., Hirano, N., and Kurz, M. D.: Noble gas isotopic compositions of seamount lavas from the central Chile trench: Implications for petit-spot volcanism and the lithosphere asthenosphere boundary, Earth Planet. Sc. Lett., 552, 116611,, 2020. 

Yang, H.-J., Frey, F. A., and Clague, D. A.: Constraints on the Source Components of Lavas Forming the Hawaiian North Arch and Honolulu Volcanics, J. Petrol., 44, 603–627,, 2003.  

Yoshino, T., Matsuzaki, T., Yamashita, S., and Katsura T.: Hydrous olivine unable to account for conductivity anomaly at the top of the asthenosphere, Nature, 443, 973–976,, 2006. 

Zakharov, D. O., Tanaka, R., Butterfield, D. A., and Nakamura, E.: A New Insight Into Seawater–Basalt Exchange Reactions Based on Combined δ18OΔ17O87Sr/86Sr Values of Hydrothermal Fluids From the Axial Seamount Volcano, Pacific Ocean, Front. Earth Sci., 9, 691699,, 2021. 

Zhang, F., Lin, J., and Zhan, W.: Variations in oceanic plate bending along the Mariana trench, Earth Planet. Sc. Lett., 401, 206–214,, 2014. 

Zhang, G., Wang, S., Huang, S., Zhan, M., and Yao, J.: CO2-rich rejuvenated stage lavas on Hawaiian Islands, Geochem. Geophy. Geosy., 23, e2022GC010525,, 2022. 

Zhang, G. L., Chen, L. H., Jackson, M., and Hofmann, A. W.: Evolution of carbonated melt to alkali basalt in the South China Sea, Nat. Geosci., 10, 229–235,, 2017. 

Zhang, J., Xu, M., and Sun, Z.: Lithospheric flexural modelling of the seaward and trenchward of the subducting oceanic plates, Int. Geol. Rev., 62, 908–923,, 2020. 

Zhang, W., Johnston, S., and Currie, C. A.: Kimberlite magmatism induced by west-dipping subduction of the North American plate, Geology, 47, 395–398,, 2019. 

Zhong, Y., Zhang, G.-L., Zhong, L.-F., Chen, L.-H., and Wang, X.-J.: Post-spreading volcanism triggered by CO2 along the South China Sea fossil spreading axis, Lithos, 404–405, 106478,, 2021. 

Zindler, A. and Hart, S.: Chemical geodynamics, Annu. Rev. Earth Pl. Sci., 14, 493–571,, 1986. 

Short summary
Plate tectonics theory is the motion of rocky plates (lithosphere) over ductile zones (asthenosphere). The causes of the lithosphere–asthenosphere boundary (LAB) are controversial; however, petit-spot volcanism supports the presence of melt at the LAB. We conducted geochemistry, geochronology, and geochemical modeling of petit-spot volcanoes on the western Pacific Plate, and the results suggested that carbonatite melt and recycled oceanic crust induced the partial melting at the LAB.