(Table 1) Osmium concentrations and isotopic ratios for ODP Leg 125 peridotites

Mantle peridotites drilled from the Izu-Bonin-Mariana forearc have unradiogenic 187Os/188Os ratios (0.1193 to 0.1273), which give Proterozoic model ages of 820 to 1230 million years ago. If these peridotites are residues from magmatism during the initiation of subduction 40 to 48 million years ago, then the mantle that melted was much more depleted in incompatible elements than the source of mid-ocean ridge basalts (MORB). This result indicates that osmium isotopes record information about ancient melting events in the convecting upper mantle not recorded by incompatible lithophile isotope tracers. Subduction zones may be a graveyard for ancient depleted mantle material, and portions of the convecting upper mantle may be less radiogenic in osmium isotopes than previously recognized.

Geochemistry of serpentenites from DSDP Hole 84-566C

We studied a unique chrysotile-antigorite serpentinite, drilled on Deep Sea Drilling Project Leg 84 (Site 566) in the Guatemala forearc. Our in situ major and trace element data provide new constraints on possible reactions and associated trace element mobilisation during shallow serpentinite subduction. Chrysotile of the studied serpentinite, formed by the hydration of an upper mantle peridotite precursor, is partially replaced by antigorite (alone) which also occurs in 0.5 mm wide unoriented veins crosscutting the rock. Based on textural relationships and the P-T-X stability of the rock forming phases, the replacement of chrysotile by antigorite occurred at T < 300 °C, due to interaction between the chrysotile-serpentinite and an aqueous fluid. A comparison of the chemical compositions of reactant and product phases reveals that about 90% of the Cl, more than 80% of the B and about 50% of the Sr hosted originally by chrysotile was lost during fluid-assisted chrysotile-to-antigorite transformation and accompanying partial dehydration, and documents the much lower affinity of antigorite for trace element uptake than that of chrysotile. The fluid-assisted chrysotile-to-antigorite transformation and associated trace element loss documented here can occur in the shallow (< 30 km) region of subduction zones. This transformation decreases notably the Cl and B inventory of subducting serpentinites, which are regarded as one of the most important carriers of these elements into subduction zones. The evolution of serpentinites during initial subduction stages thus appears to be critical in the recycling of specific trace elements such as B or Cl from forearc to subarc depths.

Rhenium and platinum-group element concentrations and Os isotope data for Philippine Sea Plate basalts from DSDP Legs 31, 58, and 59

Platinum-group elements (PGE), rhenium and osmium isotope data are reported for basalts from Deep Sea Drilling Project cores in the Philippine Sea Plate (PSP). Lithophile trace element and isotopic characteristics indicate a range of source components including DMM, EMII and subduction-enriched mantle. MORB-like basalts possess smooth, inclined chondrite-normalised PGE patterns with high palladium-PGE/iridium-PGE ratios, consistent with previously published data for MORB, and with the inferred compatibility of PGE. In contrast, while basalts with EMII-type lithophile element chemistry possess high Pt/Ir ratios, many have much lower Pd/Ir and unusually high Ru/Ir of >10. Similarly, back-arc samples from the Shikoku and Parece-Vela basins have very high Ru/Ir ratios (>30) and Pd/Ir as low as 1.1. Such extreme Pd/Ir and Ru/Ir ratios have not been previously reported in mafic volcanic suites and cannot be easily explained by variable degrees of melting, fractional crystallisation or by a shallow-level process such as alteration or degassing. The data appear most consistent with sampling of at least two mantle components with distinct PGE compositions. Peridotites with the required PGE characteristics (i.e. low Pd, but relatively high Ru and Re) have not been documented in oceanic mantle, but have been found in sub-continental mantle lithosphere and are the result of considerable melt depletion and selective metasomatic enrichment (mainly Re). The long-term presence of subduction zones surrounding the Philippine Sea Plate makes this a prime location for metasomatic enrichment of mantle, either through fluid enrichment or infiltration by small melt fractions. The Re-Os isotope data are difficult to interpret with confidence due to low Os concentrations in most samples and the uncertainty in sample age. Data for Site 444A (Shikoku Basin) give an age of 17.7+/-1.3 Ma (MSWD = 14), consistent with the proposed age of basement at the site and thus provides the first robust radiometric age for these samples. The initial 187Os/188Os of 0.1298+/-0.0069 is consistent with global MORB, and precludes significant metasomatic enrichment of Os by radiogenic slab fluids. Re-Os data for Sites 446A (two suites, Daito Basin) and 450 (Parece-Vela Basin) indicate ages of 73, 68 and 43 Ma, which are respectively, 30, 17 and >12 Ma older than previously proposed ages. The alkalic and tholeiitic suites from Site 446A define regression lines with different 187Os/188Osinitial (0.170+/-0.033 and 0.112+/-0.024, respectively) which could perhaps be explained by preferential sampling of interstitial, metasomatic sulphides (with higher time-integrated Re/Os ratios) by smaller percentage alkalic melts. One sample, with lithophile elements indistinguishable from MORB, is Os-rich (146 pg/g) and has an initial 187Os/188Os of 0.1594, which is at the upper limit of the accepted OIB range. Given the Os-rich nature of this sample and the lack of evidence for subduction or recycled crust inputs, this osmium isotope ratio likely reflects heterogeneity in the DMM. The dataset as a whole is a striking indication of the possible PGE and Os isotope variability within a region of mantle that has experienced a complex tectonic history.

Shear experiments from IODP Expedition 322 and 333 samples

The mechanical behavior of the plate boundary fault zone is of paramount importance in subduction zones, because it controls megathrust earthquake nucleation and propagation as well as the structural style of the forearc. In the Nankai area along the NanTroSEIZE (Kumano) drilling transect offshore SW Japan, a heterogeneous sedimentary sequence overlying the oceanic crust enters the subduction zone. In order to predict how variations in lithology, and thus mechanical properties, affect the formation and evolution of the plate boundary fault, we conducted laboratory tests measuring the shear strengths of sediments approaching the trench covering each major lithological sedimentary unit. We observe that shear strength increases nonlinearly with depth, such that the (apparent) coefficient of friction decreases. In combination with a critical taper analysis, the results imply that the plate boundary position is located on the main frontal thrust. Further landward, the plate boundary is expected to step down into progressively lower stratigraphic units, assisted by moderately elevated pore pressures. As seismogenic depths are approached, the décollement may further step down to lower volcaniclastic or pelagic strata but this requires specific overpressure conditions. High-taper angle and elevated strengths in the toe region may be local features restricted to the Kumano transect.

Composition of matrix glasses and melt inclusions of fallout tephra from ODP Hole 125-782A

We studied the systematics of Cl, F and H2O in Izu arc front volcanic rocks using basaltic through rhyolitic glass shards and melt inclusions (Izu glasses) from Oligocene to Quaternary distal fallout tephra. These glasses are low-K basalts to rhyolites that are equivalent to the Quaternary lavas of the Izu arc front (Izu VF). Most of the Izu glasses have Cl ~400-4000 ppm and F ~70-400 ppm (normal-group glasses). Rare andesitic melt inclusions (halogen-rich andesites; HRA) have very high abundances of Cl (~6600-8600 ppm) and F (~780-910 ppm), but their contents of incompatible large ion lithophile elements (LILE) are similar to the normal-group glasses. The preeruptive H2O of basalt to andesite melt inclusions in plagioclase is estimated to range from ~2 to ~10 wt% H2O. The Izu magmas should be undersaturated in H2O and the halogens at their preferred levels of crystallization in the middle to lower crust (~3 to ~11 kbar, ~820° to ~1200°C). A substantial portion of the original H2O is lost due to degassing during the final ascent to surface. By contrast, halogen loss is minor, except for loss of Cl from siliceous dacitic and rhyolitic compositions. The behavior of Cl, F and H2O in undegassed melts resembles the fluid mobile LILE (e.g.; K, Rb, Cs, Ba, U, Pb, Li). Most of the Cl (>99%), H2O (>95%) and F (>53%) in the Izu VF melts appear to originate from the subducting slab. At arc front depths, the slab fluid contains Cl = 0.94+/-0.25 wt%, F = 990+/-270 ppm and H2O = 25+/-7 wt%. If the subducting sediment and the altered basaltic crust were the only slab sources, then the subducted Cl appears to be almost entirely recycled at the Izu arc (~77-129%). Conversely, H2O (~13-22% recycled at arc) and F (~4-6% recycled) must be either lost during shallow subduction or retained in the slab to greater depths. If a seawater-impregnated serpentinite layer below the basaltic crust were an additional source of Cl and H2O, the calculated percentage of Cl and H2O recycled at arc would be lower. Extrapolating the Izu data to the total length of global arcs (~37000 km), the global arc outflux of fluid-recycled Cl and H2O at subduction zones amounts to Cl ~2.9-3.8 mln ton/yr and H2O ~70-100 mln ton/yr, respectively - comparable to previous estimates. Further, we obtain a first estimate of global arc outflux of fluid-recycled F of ~0.3-0.4 mln ton/yr. Despite the inherent uncertainties, our results support models suggesting that the slab becomes strongly depleted in Cl and H2O in subduction zones. In contrast, much of the subducted F appears to be returned to the deep mantle, implying efficient fractionation of Cl and H2O from F during the subduction process. However, if slab devolatilization produces slab fluids with high Cl/F (~9.5), slab melting will still produce components with low Cl/F ratios (~0.9), similar to those characteristic of the upper continental crust (Cl/F ~0.3-0.9).

A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Geochemistry of oceanic crust beneath Costa Rica

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge (Feigenson et al., 2004, doi:10.1029/2003GC000621; Herrstom et al., 1995, doi:10.1130/0091-7613(1995)023<0617:VILCAW>2.3.CO;2; Abratis and Woerner, 2001) or eroded fore-arc complexes (Goss and Kay, 2006, doi:10.1029/2005GC001163) but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (~63-190 mm/yr) and are comparable to the magnitude of subducting Cocos plate motion (approx85 mm/yr). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.

Tephrostratigraphy and radiometric dating of authigenic carbonates

At convergent margins, fluids rise through the forearc in response to consolidation of the upper plate and dewatering of the subducting plate, and produce various cold-seep-related features on the seafloor (mud diapirs, mud mounds). At the Central American forearc, authigenic carbonates precipitated from rising fluids within such structures during active venting while typical mixed-mud sediments were ejected onto the surrounding seafloor where they became intercalated with normal pelagic background sediments, indicating that mud mounds evolved unsteadily through alternating active and inactive phases. Intercalated regional ash layers from Plinian eruptions at the Central American volcanic arc provide time marks that constrain the ages of mud ejection activity. U/Th dating of drill core samples of authigenic carbonate caps of mud mounds yields ages agreeing well with those constrained by ash layers and showing that carbonate caps grow inward rather than outward during active venting. Both dating approaches show that offshore Nicaragua and Costa Rica (1) active and inactive phases can occur simultaneously at neighboring mounds, (2) mounds along the forearc have individual histories of activity, but there are distinct time intervals when nearly all mounds have been active or inactive, (3) lifetimes of mounds reach several hundred thousand years, and (4) highly active periods last 10-50 k.y. with intervening periods of >10 k.y. of relative quiescence.

(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 pore water and sediments offshore Nicaragua and Costa Rica

Four volcanic ash-bearing marine sediment cores and one ash-free reference core were examined during research cruise RV Meteor 54/2 offshore Nicaragua and Costa Rica to investigate the chemical composition of pore waters related to volcanic ash alteration. Sediments were composed of terrigenous matter derived from the adjacent continent and contained several distinct ash layers. Biogenic opal and carbonate were only minor components. The terrigenous fraction was mainly composed of smectite and other clay minerals while the pore water composition was strongly affected by the anaerobic degradation of particulate organic matter via microbial sulphate reduction. The alteration of volcanic matter showed only a minor effect on major element concentrations in pore waters. This is in contrast to prior studies based on long sediment cores taken during the DSDP, where deep sediments always showed distinct signs of volcanic ash alteration. The missing signal of ash alteration is probably caused by low reaction rates and the high background concentration of major dissolved ions in the seawater-derived pore fluids. Dissolved silica concentrations were, however, significantly enriched in ash-bearing cores and showed no relation to the low but variable contents of biogenic opal. Hence, the data suggest that silica concentrations were enhanced by ash dissolution. Thus, the dissolved silica profile measured in one of the sediment cores was used to derive the in-situ dissolution rate of volcanic glass particles in marine sediments. A non-steady state model was run over a period of 43 kyr applying a constant pH of 7.30 and a dissolved Al concentration of 0.05 ?M. The kinetic constant (AA) was varied systematically to fit the model to the measured dissolved silica-depth profile. The best fit to the data was obtained applying AA = 1.3 * 10**-U9 mol of Si/cm**2/ s. This in-situ rate of ash dissolution at the seafloor is three orders of magnitude smaller than the rate of ash dissolution determined in previous laboratory experiments. Our results therefore imply that field investigations are necessary to accurately predict natural dissolution rates of volcanic glasses in marine sediments.