Geochemistry of serpentinized peridotites from ODP Holes 207-1272A and 207-1274A

Despite the key importance of altered oceanic mantle as a repository and carrier of light elements (B, Li, and Be) to depth, its inventory of these elements has hardly been explored and quantified. In order to constrain the systematics and budget of these elements we have studied samples of highly serpentinized (>50%) spinel harzburgite drilled at the Mid-Atlantic Ridge (Fifteen-Twenty Fracture zone, ODP Leg 209, Sites 1272A and 1274A). In-situ analysis by secondary ion mass spectrometry reveals that the B, Li and Be contents of mantle minerals (olivine, orthopyroxene, and clinopyroxene) remain unchanged during serpentinization. B and Li abundances largely correspond to those of unaltered mantle minerals whereas Be is close to the detection limit. The Li contents of clinopyroxene are slightly higher (0.44-2.8 µg/g) compared to unaltered mantle clinopyroxene, and olivine and clinopyroxene show an inverse Li partitioning compared to literature data. These findings along with textural observations and major element composition obtained from microprobe analysis suggest reaction of the peridotites with a mafic silicate melt before serpentinization. Serpentine minerals are enriched in B (most values between 10 and 100 µg/g), depleted in Li (most values below 1 µg/g) compared to the primary phases, with considerable variation within and between samples. Be is at the detection limit. Analysis of whole rock samples by prompt gamma activation shows that serpentinization tends to increase B (10.4-65.0 µg/g), H2O and Cl contents and to lower Li contents (0.07-3.37 µg/g) of peridotites, implying that-contrary to alteration of oceanic crust-B is fractionated from Li and that the B and Li inventory should depend essentially on rock-water ratios. Based on our results and on literature data, we calculate the inventory of B and Li contained in the oceanic lithosphere, and its partitioning between crust and mantle as a function of plate characteristics. We model four cases, an ODP Leg 209-type lithosphere with almost no igneous crust, and a Semail-type lithosphere with a thick igneous crust, both at 1 and 75 Ma, respectively. The results show that the Li contents of the oceanic lithosphere are highly variable (17-307 kg in a column of 1 m * 1 m * thickness of the lithosphere (kg/col)). They are controlled by the primary mantle phases and by altered crust, whereas the B contents (25-904 kg/col) depend entirely on serpentinization. In all cases, large quantities of B reside in the uppermost part of the plate and could hence be easily liberated during slab dehydration. The most prominent input of Li into subduction zones is to be expected from Semail-type lithosphere because most of the Li is stored at shallow levels in the plate. Subducting an ODP Leg 209-type lithosphere would mean only very little Li contribution from the slab. Serpentinized mantle thus plays an important role in B recycling in subduction zones, but it is of lesser importance for Li.

Stable chlorine isotopes in pore waters from DSDP Hole 79-546 and ODP Hole 131-808C

Stable Cl isotope ratios, measured in marine pore waters associated with the Barbados and Nankai subduction zones, extend significantly (to ~-8 per mil) the range of d37Cl values reported for natural waters. These relatively large negative values, together with geologic and chemical evidence from Barbados and Nankai and recent laboratory data showing that hydrous silicate minerals (i.e., those with structural OH sites) are enriched up to 7.5 per mil in 37Cl relative to seawater, strongly suggest that the isotopic composition of Cl in pore waters from subduction zones reflects diagenetic and metamorphic dehydration and transformation reactions. These reactions involve clays and/or other hydrous silicate phases at depth in the fluid source regions. Chlorine therefore cannot be considered geochemically conservative in these systems. The uptake of Cl by hydrous phases provides a mechanism by which Cl can be cycled into the mantle through subduction zones. Thus, stable Cl isotopes should help in determining the extent to which Cl and companion excess volatiles like H2O and CO2 cycle between the crust and mantle.

(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.