Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Geochemistry of lavas from the Lau-Tonga arc and back-arc system

Detailed comparison of mineralogy, and major and trace geochemistry are presented for the modern Lau Basin spreading centers, the Sites 834-839 lavas, the modern Tonga-Kermadec arc volcanics, the northern Tongan boninites, and the Lau Ridge volcanics. The data clearly confirm the variations from near normal mid-ocean-ridge basalt (N-MORB) chemistries (e.g., Site 834, Central Lau Spreading Center) to strongly arc-like (e.g., Site 839, Valu Fa), the latter closely comparable to the modern arc volcanoes. Sites 835 and 836 and the East Lau Spreading Center represent transitional chemistries. Bulk compositions range from andesitic to basaltic, but lavas from Sites 834 and 836 and the Central Lau Spreading Center extend toward more silica-undersaturated compositions. The Valu Fa and modern Tonga-Kermadec arc lavas, in contrast, are dominated by basaltic andesites. The phenocryst and groundmass mineralogies show the strong arc-like affinities of the Site 839 lavas, which are also characterized by the existence of very magnesian olivines (up to Fo90-92) and Cr-rich spinels in Units 3 and 6, and highly anorthitic plagioclases in Units 2 and 9. The regional patterns of mineralogical and geochemical variations are interpreted in terms of two competing processes affecting the inferred magma sources: (1) mantle depletion processes, caused by previous melt extractions linked to backarc magmatism, and (2) enrichment in large-ion-lithophile elements, caused by a subduction contribution. A general trend of increasing depletion is inferred both eastward across the Lau Basin toward the modern arc, and northward along the Tongan (and Kermadec) Arc. Numerical modeling suggests that multistage magma extraction can explain the low abundances (relative to N-MORB) of elements such as Nb, Ta, and Ti, known to be characteristic of island arc magmas. It is further suggested that a subduction jump following prolonged slab rollback could account for the initiation of the Lau Basin opening, plausibly allowing a later influx of new mantle, as required by the recognition of a two-stage opening of the Lau Basin.

Rapid magmatic processes accompany arc-continent collision: the Western Bismarck arc, Papua New Guinea

New U-Th-Ra, major and trace element, and Sr-Nd-Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc-continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc-continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated Pb-208/Pb-204. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. U-238 and Ra-226 excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.

Geochemistry of sediments beneath the Tonga-Kermadec arc

The fate of subducted sediment and the extent to which it is dehydrated and/or melted before incorporation into arc lavas has profound implications for the thermo-mechanical nature of the mantle wedge and models for crustal evolution. In order to address these issues, we have undertaken the first measurements of 10Be and light elements in lavas from the Tonga-Kermadec arc and the sediment profile at DSDP site 204 outboard of the trench. The 10Be/9Be ratios in the Tonga lavas are lower than predicted from flux models but can be explained if (a) previously estimated sediment contributions are too high by a factor of 2-10, (b) the top 1-22 m of the incoming sediment is accreted, (c) large amounts of sediment erosion are proposed, or (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Tonga. The lavas form negative Th/Be-Li/Be arrays that extend from a depleted mantle source composition to lower Th/Be and Li/Be ratios than that of the bulk sediment. Thus, these arrays are not easily explained by bulk sediment addition and, using partition coefficients derived from experiments on the in-coming sediment, we show that they are also unlikely to result from fluid released during dehydration of the sediment (or altered oceanic crust). However, partial melts of the dehydrated sediment residue formed at ~800 °C during the breakdown of amphibole +/- plagioclase and in the absence of cordierite have significantly lowered Th/Be ratios. The lava arrays can be successfully modelled as 10-15% partial melts of depleted mantle after it has been enriched by the addition of 0.2-2% of these partial melts. Phase relations suggest that this requires that the top of the subducting crust reaches temperatures of ~800 °C by the time it attains ~ 80 km depth which is in excellent agreement with the results of recent numerical models incorporating a temperature-dependent mantle viscosity. Under these conditions the wet basalt solidus is also crossed yet there is no recognisable eclogitic signal in the lavas suggesting that on-going dehydration or strong thermal gradients in the upper part of the subducting plate inhibit partialmelting of the altered oceanic crust.

Utility of thallium and geochemical analysis of Mariana arc sediments from ODP Hole 129-801A and 129-801B

We provide the first exploration of thallium (Tl) abundances and stable isotope compositions as potential tracers during arc lava genesis. We present a case study of lavas from the Central Island Province (CIP) of the Mariana arc, supplemented by representative sedimentary and altered oceanic crust (AOC) inputs from ODP Leg 129 Hole 801 outboard of the Mariana trench. Given the large Tl concentration contrast between the mantle and subduction inputs coupled with previously published distinctive Tl isotope signatures of sediment and AOC, the Tl isotope system has great potential to distinguish different inputs to arc lavas. Furthermore, CIP lavas have well-established inter island variability, providing excellent context for the examination of Tl as a new stable isotope tracer. In contrast to previous work (Nielsen et al., 2006b), we do not observe Tl enrichment or light epsilon 205Tl (where epsilon 205Tl is the deviation in parts per 10,000 of a sample 205Tl/203Tl ratio compared to NIST SRM 997 Tl standard) in the Jurassic-aged altered mafic ocean crust subducting outboard of the Marianas (epsilon 205Tl = - 4.4 to 0). The lack of a distinctive epsilon 205Tl signature may be related to secular changes in ocean chemistry. Sediments representative of the major lithologies from ODP Hole Leg 129 801 have 1-2 orders of magnitude of Tl enrichment compared to the CIP lavas, but do not record heavy signatures (epsilon 205Tl = - 3.0 to + 0.4), as previously found in similar sediment types (epsilon 205Tl > + 2.5; Rehkämper et al., 2004). We find a restricted range of epsilon 205Tl = - 1.8 to - 0.4 in CIP lavas, which overlaps with MORB. One lava from Guguan falls outside this range with epsilon 205Tl = + 1.2. Coupled Cs, Tl and Pb systematics of Guguan lavas suggests that this heavy Tl isotope composition may be due to preferential degassing of isotopically light Tl. In general, the low Tl concentrations and limited isotopic range in the CIP lavas is likely due to the unexpectedly narrow range of epsilon 205Tl found in Mariana subduction inputs, coupled with volcaniclastic, rather than pelagic sediment as the dominant source of Tl. Much work remains to better understand the controls on Tl processing through a subduction zone. For example, Tl could be retained in residual phengite, offering the potential exploration of Cs/Tl ratios as a slab thermometer. However, data for Tl partitioning in phengite (and other micas) is required before developing this application further. Establishing a database of Tl concentrations and stable isotopes in subduction zone lavas with different thermal parameters and sedimentary inputs is required for the future use of Tl as a subduction zone tracer.

Geochemistry and radiogenic isotope ratios of the volcanic sequences from Izu-Bonin Arc, Japan

New radiogenic isotope and trace element data are presented for the volcanic sequences along 600 km of the active Izu-Bonin arc, the Oligocene Izu arc, and their associated rift basins. As with many intra-oceanic island arcs, the Pliocene-Recent Izu-Bonin frontal-arc lavas are highly depleted in Zr, Nb and the rare-earth elements relative to typical mid-ocean ridge basalt (MORB), indicating that the mantle wedge source has undergone a previous episode of melting. Ratios between these elements (such as Nb/Zr and La/Sm), as well as 143Nd/144Nd, do not vary significantly along the length of the frontal-arc. These parameters suggest that each of the arc volcanoes is derived from similar melt fractions of the mantle wedge. However, Ba/Zr, Ba/Rb and 87Sr/86Sr increase along the frontal-arc to the north. This leads us to propose that a variable enrichment in Ba and radiogenic Sr is superimposed on the mantle wedge. Sr-Nd and Pb-Nd isotope variation indicate that both Sr and Pb become more radiogenic after fluid addition. However, Pb isotope ratios do not correlate with increases in Pb concentration or ratios such as Ba/Zr and Nb/Pb. In other words, the Pb isotopic composition of the arc lavas appears to be independent of the amount of Pb introduced by subduction fluids into the mantle source. This buffering of Pb isotopes along the frontal-arc means that the isotopic composition of the lavas is indistinguishable from that of the fluid. Isotopic mixing models presented for the arc are only illustrative of the many plausible combinations of components and quantities. Despite this, we are able to determine that the mantle wedge has isotopic characteristics similar to Indian Ocean MORB, and that the subduction-fluid solute is primarily derived from subducted oceanic basalt with a <2% contribution from subducted sediment. Lavas in the Oligocene Izu arc and fore-arc basin were derived from a mantle wedge of similar composition to the active arc. Despite levels of Pb enrichment comparable to those of the modern arc, the Pb isotopes of the Oligocene volcanics indicate a lower sediment input into the melting region.

Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

Focused fluid transfer through the mantle above subduction zones

Volcanic arcs above subduction zones are enriched in volatiles and fluid-mobile elements with respect to mid-oceanic ridge basalts. There is general consensus that this particular subduction zone signature is generated by fluid-induced extraction of these elements from subducted oceanic crust and its sedimentary cover. However, how these fluids are transferred through the mantle wedge to the locus of partial melting and what modification the fluids will experience is unresolved. Here we investigate the interaction of slab fluids with the mantle wedge through a series of high-pressure experiments. We explore two end-member processes of focused and porous reactive flow of hydrous slab melts through the mantle. Transfer by porous flow leads to the formation of hydrous minerals that sequester fluid-mobile elements and residual fluids characterized by trace element patterns inconsistent with typical arc lavas. In contrast, no hydrous minerals are formed in the reaction zone of experiments mimicking focused flow, and the typical trace element signature acquired during fluid extraction from the slab is preserved, indicating that this is an efficient process for element transfer through the mantle wedge.

Chemical and isotopic compositions of sediments from DSDP Hole 21-204

The highly depleted intra-oceanic Tonga-Kermadec island arc forms an endmember of arc systems and a unique location in which to isolate the effects of the slab flux. High precision TIMS uranium, thorium, strontium, neodymium, and lead isotopes, along with complete major and trace element data, have been obtained on an extensive sample set comprising fifty-eight lavas along the arc as well as nineteen samples of the subducting sediments at DSDP site 204 just to the east of the Tonga-Kermadec trench. Ca/Ti and Al/Ti ratios extend from values appropriate to an N-MORB source in the southern Kermadecs to very high ratios in Tonga interpreted to reflect increasing degrees of depletion of the mantle wedge due to backarc basalt extraction. The isotope data emphasize the need for four components in the petrogenesis of the lavas: (1) the mantle wedge; (2) a component with elevated 207Pb/204Pb towards which the Kermadec and southern Tongan lavas extend; (3) a component characterised by high 206Pb/204Pb, Ta/Nd, and low 143Nd/144Nd observed only in the northernmost Tongan islands of Tafahi and Niuatoputapu; (4) a fluid component characterised by strong enrichments of Rb, Ba, U, K, Ph, and Sr, relative to Th, Zr, and the REE and producing large 238U excesses ((230Th/238U) = 0.8-0.5) in the more depleted lavas. The mantle wedge (Component 1) is isotopically similar to the source of the Lau BABB. Component 2 is average pelagic sediment on the downgoing Pacific plate as observed at DSDP sites 595/596 and in the upper sections of the sediment pile at DSDP site 204. Mass balance calculations indicate that less than 0.5% is recycled into the arc lavas; essentially all the subducted sediment is returned to the upper mantle (~0.03 km**3/yr). Exceptionally low concentrations of Ta and Nb relative to Th and the LREE requires that this sediment component is added as a partial melt which was in equilibrium with residual rutile or ilmenite. Component 3 is identified as volcaniclastics from the Louisville Ridge which comprise the lower 44 m of the sediment section intersected at DSDP site 204. These volcaniclastics are spatially restricted to the vicinity of the Louisville Ridge and provide a unique sediment tracer which can be used to show that it takes 4 Myr from the time of subduction to its first appearance in the arc lava signature. Component 4, the fluid contribution to the lava source is inferred to contribute ~1 ppm Rb, 10 ppm Ba, 0.02 ppm U, 600 ppm K, 0.2 ppm Ph, and 30 ppm Sr. It has 87Sr/86Sr = 0.7035 and 206Pb/204Pb = 18.5 and thus it is inferred to have been derived from dehydration of the subducting altered oceanic crust. U-Th isotope disequilibria reflect the time since fluid release from the subducting slab and a reference line through the lowest (230Th/232Th) lavas constrains this to be 30000-50000 yr. The U-Th and Th-Ra isotope systematics are decoupled, and it is suggested that Th-Ra isotope disequilibria record the time since partial melting and thus indicate rapid channelled magma ascent. Olivine gabbro xenoliths from Raoul are interpreted as cumulates to their host lavas with which they form zero age U-Th isochrons indicating that minimal time was spent in magma chambers. The subduction signature is not observed in lavas from the backarc island of Niuafo'ou. These were derived from partial melting of fertile peridotite at 130-160 km depth with melt rates around 0.0002 kg/m**3/yr.

Chemical and boron isotopic compositions of volcanics from ODP Sites 135-834 to 135-840

The influence of fluid flux on petrogenesis in the Tonga-Kermadec Arc was investigated using ion microprobe measurements of B/Be and boron isotope ratios (11B/10B) to document the source and relative volumes of the fluids released from the subducting oceanic plate. We analyzed young lavas from eight different islands along the Tonga-Kermadec Arc, as well as glass shards in volcanic sediments from Ocean Drilling Program (ODP) Site 840, which record the variations in the chemistry of Tonga magmatism since 7 Ma. B/Be is variable (5.8-122), in young Tonga-Kermadec Arc lavas. In contrast, glass shards from around 3 to 4 Ma old volcanic sediments at Site 840 have the highest B/Be values yet reported for arc lavas (18-607). These values are too high to be related simply to a sediment influence on petrogenesis. Together with very high d11B values (-11.6 to +37.5) for the same shards and lavas these data indicate that most of the B is derived from fluid escaped from the subducting altered Pacific oceanic crust, rather than from sediment. High d11B values also reflect large degrees of isotopic fractionation in this cold fast subduction zone. Lower d11B values noted in the Kermadec Arc (17 to -4.4) are related to the influence of sediment eroded from New Zealand and slower convergence. High fluid flux (B/Be) is synchronous in Tonga and the Marianas at 3 to 4 Ma and may be related to acceleration of the Pacific Plate just prior to this time. The timing of maximum B/Be at 3 to 4 Ma correlates with maximum light rare earth (LREE) and high field strength element depletion. This suggests maximum degrees of partial melting at this time. Although thinning of the arc lithosphere during rifting to form the Lau Basin is expected to influence the arc geochemistry, variable aqueous fluid flux from the subducting plate alone appears capable of explaining boron and other trace element systematics in the Tonga-Kermadec Arc with no indication of slab melting.

Geochemistry of Lesser Antilles forearc sediments

Oceanic sediments contain the products of erosion of continental crust, biologic activity and chemical precipitation. These processes create a large diversity of their chemical and isotopic compositions. Here we focus on the influence of the distance from a continental platform on the trace element and isotopic compositions of sediments deposited on the ocean floor and highlight the role of zircons in decoupling high-field strength elements and Hf isotopic compositions from other trace elements and Nd isotopic compositions. We report major and trace element concentrations as well as Sr and Hf isotopic data for 80 sediments from the Lesser Antilles forearc region. The trace-element characteristics and the Sr and Hf isotopic compositions are generally dominated by detrital material from the continental crust but are also variably influenced by chemical or biogenic carbonate and pure biogenic silica. Next to the South American continent, at DSDP Site 144 and on Barbados Island, sediments, coarse quartz arenites, exhibit marked Zr and Hf excesses that we attribute to the presence of zircon. In contrast, the sediments from DSDP Site 543, which were deposited farther away from the continental platform, consist of fine clay and they show strong deficiencies in Zr and Hf. The enrichment or depletion of Zr-Hf is coupled to large changes in Hf isotopic compositions (-30 < epsilon-Hf < +4) that vary independently from the Nd isotopes. We interpret this feature as a clear expression of the "zircon effect" suggested by Patchett and coauthors in 1984. Zircon-rich sediments deposited next to the South American continent have very low epsilon-Hf values inherited from old zircons. In contrast, in detrital clay-rich sediments deposited a few hundred kilometers farther north, the mineral fraction is devoid of zircon and they have drastically higher epsilon-Hf values inherited from finer, clay-rich continental material. In the two DSDP sites, average Hf isotopes are very unradiogenic relative to other oceanic sediments worldwide (epsilon-Hf = -14.4 and -7.4) and they define the low Hf end member of the sedimentary field in Hf-Nd space. Their compositions correspond to end members that, when mixed with mantle, are able to reproduce the pattern of volcanic rocks from the Lesser Antilles. More generally, we find a relationship between Nb/Zr ratios and the vertical deviation of Hf isotope ratios from the Nd-Hf terrestrial array and we suggest that this relationship can be used as a tool to distinguish sediment input from fractionation during melting during the formation of arc lavas.

Trace element abundances and Sr and Nd isotope compositions in bulk sediments from DSDP Legs 6 and 20

Rare earth element and Nd isotopic data for ten representative samples of Lower Cretaceous to Miocene pelagic sediments from the western Pacific indicate a wide range of compositions for sediments being subducted beneath the Mariana and Volcano arcs. All samples are enriched in light rare earth elements and show negative Eu and Ce anomalies. The values of e-Nd range from +0.6 to -7.3. These data are used to calculate the Bulk Western Pacific Sediment (BWPS), which is characterized by low Sr/Nd (10), Ba/La (13), and e-Nd (-5.2) and high 87Sr/86Sr (0.7078) compared to that of Mariana and Volcano arc lavas. This composite sediment is used to refine a mixing model for the origin of Mariana and Volcano arc melts. Some lavas from the northern Mariana Arc have Ba/La higher than that of BWPS, which indicates that a third component is required. The high Ba/La in the mantle source for these lavas is interpreted to result from multiple episodes of fluid fractionation. The mixing model indicates that a minor amount of sediment and a low proportion of metasomatic fluid fluxes the mantle source at a late stage when the subarc mantle is already highly metasomatized. This model also suggests that the mantle source for arc melts is affected more by metasomatic fluids than by melting or bulk mixing of sediments.

Element and isotope compostion of basalts and sediments from the Japan Sea

Ocean Drilling Program Legs 127 and 128 in the Yamato Basin of the Japan Sea, a Miocene-age back-arc basin in the western Pacific Ocean, recovered incompatible-element-depleted and enriched tholeiitic dolerites and basalts from the basin floor, which provide evidence of a significant sedimentary component in their mantle source. Isotopically, the volcanic rocks cover a wide range of compositions (e.g., 87Sr/86Sr = 0.70369 - 0.70503, 206Pb/204Pb = 17.65 - 18.36) and define a mixing trend between a depleted mantle (DM) component and an enriched component with the composition of EM II. At Site 797, the combined isotope and trace element systematics support a model of two component mixing between depleted, MORB-like mantle and Pacific pelagic sediments. A best estimate of the composition of the sedimentary component has been determined by analyzing samples of differing lithology from DSDP Sites 579 and 581 in the western Pacific, east of the Japan arc. The sediments have large depletions in the high field strength elements and are relatively enriched in the large-ion-lithophile elements, including Pb. These characteristics are mirrored, with reduced amplitudes, in Japan Sea enriched tholeiites and northeast Japan arc lavas, which strengthens the link between source enrichment and subducted sediments. However, Site 579/581 sediments have higher LILE/REE and lower HFSE/REE than the enriched component inferred fiom mixing trends at Site 797. Sub-arc devolatilization of the sediments is a process that will lower LILE/REE and raise HFSE/REE in the residual sediment, and thus this residual sediment may serve as the enriched component in the back-arc basalt source. Samples from other potential sources of an enriched, EM II-like component beneath Japan, such as the subcontinental lithosphere or crust, have isotopic compositions which overlap those of the Japan Sea tholeiites and are not "enriched" enough to be the EM II end-member.

Sr, Nd and Pb isotopic ratios of igneous rocks from DSDP holes (Table 1)

Igneous rocks from the Philippine tectonic plate recovered on Deep Sea Drilling Project Legs 31, 58 and 59 have been analyzed for Sr, Nd and Pb isotope ratios. Samples include rocks from the West Philippine Basin, Daito Basin and Benham Rise (40-60 m.y.), the Palau-Kyushu Ridge (29-44 m.y.) and the Parece Vela and Shikoku basins (17-30 m.y.). Samples from the West Philippine, Parece Vela and Shikoku basins are MORB (mid-ocean ridge basalt)-like with 87Sr/86Sr = 0.7026 - 0.7032, 143Nd/144Nd = 0.51300 - 0.51315, and 206Pb/204Pb = 17.8 - 18.1. Samples from the Daito Basin and Benham Rise are OIB (oceanic island basalt)-like with 87Sr/86Sr = 0.7038 - 0.7040, 143Nd/144Nd = 0.51285 - 0.51291 and 206Pb/204Pb = 18.8 - 19.2. All of these rocks have elevated 207Pb/204Pb and 208Pb/204Pb compared to the Northern Hemisphere Regression Line (NHRL) and have delta207Pb values of 0 to +6 and delta208Pb values of +32 to +65. Lavas from the Palau-Kyushu Ridge, a remnant island arc, have 87Sr/86Sr = 7032 - 0.7035, 143Nd/144Nd = 0.51308 - 0.51310 and 206Pb/204Pb = 18.4 - 18.5. Unlike the basin magmas erupted before and after them, these lavas plot along the NHRL and have Pb-isotope ratios similar to modern Pacific plate MORB's. This characteristic is shared by other Palau-Kyushu Arc volcanic rocks that have been sampled from submerged and subaerial portions of the Mariana fore-arc. At least four geochemically distinct magma sources are required for these Philippine plate magmas. The basin magmas tap Source 1, a MORB-mantle source that was contaminated by EMI (enriched mantle component 1 (Hart, 1988, doi:10.1016/0012-821X(88)90131-8)) and Source 2, an OIB-like mantle source with some characteristics of EMII (enriched mantle component 2 (Hart, 1988)). The arc lavas are derived from Source 3, a MORB-source or residue mantle including Sr and Pb from the subducted oceanic crust, and Source 4, MORB-source or residue mantle including a component with characteristics of HIMU (mantle component with high U/Pb (Hart, 1988)). These same sources can account for many of the isotopic characteristics of recent Philippine plate arc and basin lavas. The enriched components in these sources which are associated with the DUPAL anomaly were probably introduced into the asthenosphere from the deep mantle when the Philippine plate was located in the Southern Hemisphere 60 m.y.b.p.