Geochemistry of Central American Volcanic Arc basement samples

The Central American Volcanic Arc (CAVA) has been the subject of intensive research over the past few years, leading to a variety of distinct models for the origin of CAVA lavas with various source components. We present a new model for the NW Central American Volcanic Arc based on a comprehensive new geochemical data set (major and trace element and Sr-Nd-Pb-Hf-O isotope ratios) of mafic volcanic front (VF), behind the volcanic front (BVF) and back-arc (BA) lava and tephra samples from NW Nicaragua, Honduras, El Salvador and Guatemala. Additionally we present data on subducting Cocos Plate sediments (from DSDP Leg 67 Sites 495 and 499) and igneous oceanic crust (from DSDP Leg 67 Site 495), and Guatemalan (Chortis Block) granitic and metamorphic continental basement. We observe systematic variations in trace element and isotopic compositions both along and across the arc. The data require at least three different endmembers for the volcanism in NW Central America. (1) The NW Nicaragua VF lavas require an endmember with very high Ba/(La, Th) and U/Th, relatively radiogenic Sr, Nd and Hf but unradiogenic Pb and low d18O, reflecting a largely serpentinite-derived fluid/hydrous melt flux from the subducting slab into a depleted N-MORB type of mantle wedge. (2) The Guatemala VF and BVF mafic lavas require an enriched endmember with low Ba/(La, Th), U/Th, high d18O and radiogenic Sr and Pb but unradiogenic Nd and Hf isotope ratios. Correlations of Hf with both Nd and Pb isotopic compositions are not consistent with this endmember being subducted sediments. Granitic samples from the Chiquimula Plutonic Complex in Guatemala have the appropriate isotopic composition to serve as this endmember, but the large amounts of assimilation required to explain the isotope data are not consistent with the basaltic compositions of the volcanic rocks. In addition, mixing regressions on Nd vs. Hf and the Sr and O isotope plots do not go through the data. Therefore, we propose that this endmember could represent pyroxenites in the lithosphere (mantle and possibly lower crust), derived from parental magmas for the plutonic rocks. (3) The Honduras and Caribbean BA lavas define an isotopically depleted endmember (with unradiogenic Sr but radiogenic Nd, Hf and Pb isotope ratios), having OIB-like major and trace element compositions (e.g. low Ba/(La, Th) and U/Th, high La/Yb). This endmember is possibly derived from melting of young, recycled oceanic crust in the asthenosphere upwelling in the back-arc. Mixing between these three endmember types of magmas can explain the observed systematic geochemical variations along and across the NW Central American Arc.

Geochemistry of Pliocene-Pleistocene volcanic rocks from the Izu Arc

The igneous geochemistry of lavas and breccias from the basement of Sites 790 and 791, and pumice clasts from the Pliocene-Pleistocene sedimentary section of Sites 788, 790, 791, and 793 were studied. Arc volcanism became silicic about 1.5 m.y. before the inception of rifting in the Sumisu Rift at 2 Ma, but eruption of these silicic magmas reflects changes in stress regime, especially during the last 130,000 yr, rather than crustal anatexis. Arc magmas have had a larger proportion of slab-derived components since the inception of rifting than before, but are otherwise similar. Rift basalts and rhyolites are derived from a different source than are arc andesites to rhyolites. The rift source has less slab-derived material and is an E-MORB-like source, in contrast to an N-MORB-type source overprinted with more slab-derived material beneath the arc. Rift magma types, in the form of rare pumice and lithic clasts, preceded the rift, and the earliest magmas that erupted in the rift already differed from those of the arc. The earliest large rift eruption produced an exotic explosion breccia ("mousse") despite eruption at >1800 mbsl. Although this rock type is attributed primarily to high magmatic water content, the clasts are more MORB-like in trace element and isotopic composition than are modern Mariana Trough basalts. After rifting began, arc volcanism continued to be predominantly silicic, with individual pumice deposits containing clasts that vary in composition by about 5 wt% SiO2, or about as much as in historical eruptions of submarine Izu Arc volcanoes. The overall variations in magma composition with time during the inception of arc rifting are broadly similar in the Sumisu Rift and Lau Basin, though newly tapped OIB-type mantle seems to be present earlier during basin formation in the Sumisu than Lau case.

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.

(Figures S2a-c) Phase diagrams of the upper crust, the lower crust and the upper mantle based on averaged compositions of the oceanic crust and mantle

The age of the subducting Nazca Plate off Chile increases northwards from 0 Ma at the Chile Triple Junction (46°S) to 37 Ma at the latitude of Valparaíso (32°S). Age-related variations in the thermal state of the subducting plate impact on (a) the water influx to the subduction zone, as well as on (b) the volumes of water that are released under the continental forearc or, alternatively, carried beyond the arc. Southern Central Chile is an ideal setting to study this effect, because other factors for the subduction zone water budget appear constant. We determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the Chile Trench. The water release under forearc and arc is determined by coupling FEM thermal models of the subducting plate with stability fields of water-releasing mineral reactions for upper and lower crust and hydrated mantle. Results show that both the influx of water stored in, and the outflux of water released from upper crust, lower crust and mantle vary drastically over segment boundaries. In particular, the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments, but their release flux to the forearc is only about one fourth of the latter. This high variability over a subduction zone of < 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes. This article is protected by copyright. All rights reserved.

Geochemistry of ocean floor and forearc serpentinites

We provide new insights into the geochemistry of serpentinites from mid-ocean ridges (Mid-Atlantic Ridge and Hess Deep), passive margins (Iberia Abyssal Plain and Newfoundland) and fore-arcs (Mariana and Guatemala) based on bulk-rock and in situ mineral major and trace element compositional data collected on drill cores from the Deep Sea Drilling Project and Ocean Drilling Program. These data are important for constraining the serpentinite-hosted trace element inventory of subduction zones. Bulk serpentinites show up to several orders of magnitude enrichments in Cl, B, Sr, U, Sb, Pb, Rb, Cs and Li relative to elements of similar compatibility during mantle melting, which correspond to the highest primitive mantle-normalized B/Nb, B/Th, U/Th, Sb/Ce, Sr/Nd and Li/Y among subducted lithologies of the oceanic lithosphere (serpentinites, sediments and altered igneous oceanic crust). Among the elements showing relative enrichment, Cl and B are by far the most abundant with bulk concentrations mostly above 1000 µg/g and 30 µg/g, respectively. All other trace elements showing relative enrichments are generally present in low concentrations (µg/g level), except Sr in carbonate-bearing serpentinites (thousands of µg/g). In situ data indicate that concentrations of Cl, B, Sr, U, Sb, Rb and Cs are, and that of Li can be, increased by serpentinization. These elements are largely hosted in serpentine (lizardite and chrysotile, but not antigorite). Aragonite precipitation leads to significant enrichments in Sr, U and B, whereas calcite is important only as an Sr host. Commonly observed brucite is trace element-poor. The overall enrichment patterns are comparable among serpentinites from mid-ocean ridges, passive margins and fore-arcs, whereas the extents of enrichments are often specific to the geodynamic setting. Variability in relative trace element enrichments within a specific setting (and locality) can be several orders of magnitude. Mid-ocean ridge serpentinites often show pronounced bulk-rock U enrichment in addition to ubiquitous Cl, B and Sr enrichment. They also exhibit positive Eu anomalies on chondrite-normalized rare earth element plots. Passive margin serpentinites tend to have higher overall incompatible trace element contents than mid-ocean ridge and fore-arc serpentinites and show the highest B enrichment among all the studied serpentinites. Fore-arc serpentinites are characterized by low overall trace element contents and show the lowest Cl, but the highest Rb, Cs and Sr enrichments. Based on our data, subducted dehydrating serpentinites are likely to release fluids with high B/Nb, B/Th, U/Th, Sb/Ce and Sr/Nd, rendering them one of the potential sources of some of the characteristic trace element fingerprints of arc magmas (e.g. high B/Nb, high Sr/Nd, high Sb/Ce). However, although serpentinites are a substantial part of global subduction zone chemical cycling, owing to their low overall trace element contents (except for B and Cl) their geochemical imprint on arc magma sources (apart from addition of H2O, B and Cl) can be masked considerably by the trace element signal from subducted crustal components.

Element and isotope composition of selected volcanics from the Bicol arc, Philippines

Pliocene to recent volcanic rocks from the Bulusan volcanic complex in the southern part of the Bicol arc (Philippines) exhibit a wide compositional range (medium- to high-K basaltic-andesites, andesites and a dacite/rhyolite suite), but are characterised by large ion lithophile element enrichments and HFS element depletions typical of subduction-related rocks. Field, petrographic and geochemical data indicate that the more silicic syn- and post-caldera magmas have been influenced by intracrustal processes such as magma mixing and fractional crystallisation. However, the available data indicate that the Bicol rocks as a group exhibit relatively lower and less variable 87Sr/86Sr ratios (0.7036-0.7039) compared with many of the other subduction-related volcanics from the Philippine archipelago. The Pb isotope ratios of the Bicol volcanics appear to be unlike those of other Philippine arc segments. They typically plot within and below the data field for the Philippine Sea Basin on 207Pb/204Pb versus 206Pb/204Pb and 208Pb/204Pb versus 206Pb/204Pb diagrams, implying a pre-subduction mantle wedge similar to that sampled by the Palau Kyushu Ridge, east of the Philippine Trench. 143Nd/144Nd ratios are moderately variable (0.51285-0.51300). Low silica (<55 wt%) samples that have lower 143Nd/144Nd tend to have high Th/Nd, high Th/Nb, and moderately low Ce/Ce* ratios. Unlike some other arc segments in the Philippines (e.g. the Babuyan-Taiwan segment), there is little evidence for the involvement of subducted terrigenous sediment. Instead, the moderately low 143Nd/144Nd ratios in some of the Bicol volcanics may result from subduction of pelagic sediment (low Ce/Ce*, high Th/Nd, and high Th/Nb) and its incorporation into the mantle wedge via a slabderived partial melt.

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.

(Tables 1-3) Pb and Sr isotopic data for sediments and rocks from DSDP holes in the Mariana island-arc system

Analyses of the isotopic composition of Pb in (1) western Pacific Ocean sediments [Jurassic(?) to Pleistocene in age, including clays and biogenic oozes], (2) Pacific Ocean basaltic rocks, (3) Mariana frontal arc volcanic rocks (Eocene to Miocene), and (4) Mariana active arc volcanic rocks [Pliocene (?) to Holocene] indicate that Pacific Ocean sediments could not have been a significant component of the source material for the Mariana arc volcanic rocks. Calculations involving the average concentrations and isotopic compositions of Pb in oceanic sediments, sea-floor basaltic rocks, and the Mariana arc volcanic rocks suggest that the sediment component must have been less than 1 percent of this source material. The Pb isotopic compositions of the Mariana arc volcanic rocks lie, within experimental error, along the trend of available Pacific Ocean basalt analyses in versus 207Pb/204Pb versus 206Pb/204Pb and 208Pb/204Pb versus 206Pb/204Pb diagrams. Isotopic analyses of Pb in Pacific Ocean sediments do not lie along this trend; they have higher 207Pb/204Pb and 208Pb/204Pb values for comparable 206Pb/204Pb ratios. Clayey sediments generally have higher 208Pb/204Pb and 207Pb/204Pb ratios than biogenic oozes regardless of the age of the sediment. Comparison of combined Sr and Pb isotopic analyses for (1) mantle-derived materials erupted through oceanic crust, (2) altered ocean-floor basaltic rocks, and (3) volcanic rocks from oceanic island arcs suggests that the Mariana arc volcanic rocks were derived, at least in part, from altered Pacific lithosphere subducted beneath the Mariana arc. Unaltered basalts from the Mariana inter-arc basin (Mariana Trough) have Pb and Sr isotopic compositions that are very similar to those reported for some Hawaiian volcanic rocks but distinct from Mariana active and frontal arc compositions. These observations, in addition to existing major-and trace-element data, support a mantle origin for the interarc basin volcanic rocks. Dacites dredged from the Mariana remnant arc (South Honshu Ridge) have Pb isotopic compositions that are within experimental error of the active-arc analyses, consistent with a genetic relation.

Chemistry and boron isotope ratios of Izu arc front volcanic rocks

The process of fluid release from the subducting slab beneath the Izu arc volcanic front (Izu VF) was examined by measuring B concentrations and B isotope ratios in the Neogene fallout tephra (ODP Site 782A). Both were measured by secondary ion mass spectrometry, in a subset of matrix glasses and glassy plagioclase-hosted melt inclusions selected from material previously analyzed for major and trace elements (glasses) and radiogenic isotopes (Sr, Nd, Pb; bulk tephra). These tephra glasses have high B abundances (~10-60 ppm) and heavy delta11B values (+4.5? to +12.0?), extending the previously reported range for Izu VF rocks (delta11B, +7.0? to +7.3?). The glasses show striking negative correlations of delta11B with large ion lithophile element (LILE)/Nb ratios. These correlations cannot be explained by mixing two separate slab fluids, originating from the subducting sediment and the subducting basaltic crust, respectively (model A). Two alternative models (models B and C) are proposed. Model B proposes that the inverse correlations are inherited from altered oceanic crust (AOC), which shows a systematic decrease of B and LILE with increasing depth (from basaltic layer 2A to layer 3), paralleled by an increase in delta11B (from ~ +1? to +10? to +24?). In this model, the contribution of sedimentary B is insignificant (<4% of B in the Izu VF rocks). Model C explains the correlation as a mixture of a low-delta11B (~ +1?) 'composite' slab fluid (a mixture of metasediment- and metabasalt-derived fluids) with a metasomatized mantle wedge containing elevated B (~1-2 ppm) and heavy delta11B (~ +14?). The mantle wedge was likely metasomatized by 11B-rich fluids beneath the outer forearc, and subsequently down dragged to arc front depths by the descending slab. Pb-B isotope systematics indicate that, at arc front depths, ~ 53% of the B in the Izu VF is derived from the wedge. This implies that the heavy delta11B values of Izu VF rocks are largely a result of fluid fractionation, and do not reflect variations in slab source provenance (i.e. subducting sediment vs. basaltic crust). Since the B content of the peridotite at the outer forearc (7-58 ppm B, mean 24 +/- 16 ppm) is much higher than beneath the arc front (~1-2 ppm B), the hydrated mantle wedge must have released a B-rich fluid on its downward path. This 'wedge flux' can explain (1) the across-arc decrease in B and delta11B (e.g. Izu, Kuriles), without requiring a progressive decrease in fluid flux from the subducting slab, and (2) the thermal structure of volcanic arcs, as reflected in the B and delta11B variations of volcanic arc rocks.

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.