Implications of subduction rehydration for earth's deep water cycle

The “standard model” for the genesis of the oceans is that they are exhalations from Earth’s deep interior continually rinsed through surface rocks by the global hydrologic cycle. No general consensus exists, however, on the water distribution within the deeper mantle of the Earth. Recently Dixon et al. [2002] estimated water concentrations for some of the major mantle components and concluded that the most primitive (FOZO) are significantly wetter than the recycling associated EM or HIMU mantle components and the even drier depleted mantle source that melts to form MORB. These findings are in striking agreement with the results of numerical modeling of the global water cycle that are presented here. We find that the Dixon et al. [2002] results are consistent with a global water cycle model in which the oceans have formed by efficient outgassing of the mantle. Present-day depleted mantle will contain a small volume fraction of more primitive wet mantle in addition to drier recycling related enriched components. This scenario is consis-tent with the observation that hotspots with a FOZO-component in their source will make wetter basalts than hotspots whose mantle sources contain a larger fraction of EM and HIMU components

Across-arc geochemical variations in the Southern Volcanic Zone, Chile (34.5- 38.0°S): Constraints on Mantle Wedge and Input Compositions

Crustal assimilation (e.g. Hildreth and Moorbath, 1988) and/or subduction erosion (e.g. Stern, 1991; Kay et al., 2005) are believed to control the geochemical variations along the northern portion of the Chilean Southern Volcanic Zone. In order to evaluate these hypotheses, we present a comprehensive geochemical data set (major and trace elements and O-Sr-Nd-Hf-Pb isotopes) from Holocene primarily olivine-bearing volcanic rocks across the arc between 34.5-38.0°S, including volcanic front centers from Tinguiririca to Callaqui, the rear arc centers of Infernillo Volcanic Field, Laguna del Maule and Copahue, and extending 300 km into the backarc. We also present an equivalent data set for Chile Trench sediments outboard of this profile. The volcanic arc (including volcanic front and rear arc) samples primarily range from basalt to andesite/trachyandesite, whereas the backarc rocks are low-silica alkali basalts and trachybasalts. All samples show some characteristic subduction zone trace element enrichments and depletions, but the backarc samples show the least. Backarc basalts have higher Ce/Pb, Nb/U, Nb/Zr, and Ta/Hf, and lower Ba/Nb and Ba/La, consistent with less of a slab-derived component in the backarc and, consequently, lower degrees of mantle melting. The mantle-like δ18O in olivine and plagioclase phenocrysts (volcanic arc = 4.9-5.6 and backarc = 5.0-5.4 per mil) and lack of correlation between δ18O and indices of differentiation and other isotope ratios, argue against significant crustal assimilation. Volcanic arc and backarc samples almost completely overlap in Sr and Nd isotopic composition. High precision (double-spike) Pb isotope ratios are tightly correlated, precluding significant assimilation of older sialic crust but indicating mixing between a South Atlantic Mid Ocean-Ridge Basalt (MORB) source and a slab component derived from subducted sediments and altered oceanic crust. Hf-Nd isotope ratios define separate linear arrays for the volcanic arc and backarc, neither of which trend toward subducting sediment, possibly reflecting a primarily asthenospheric mantle array for the volcanic arc and involvement of enriched Proterozoic lithospheric mantle in the backarc. We propose a quantitative mixing model between a mixed-source, slab-derived melt and a heterogeneous mantle beneath the volcanic arc. The model is consistent with local geodynamic parameters, assuming water-saturated conditions within the slab

Fluid release from the subducted Cocos plate and partial melting of the crust deduced from magnetotelluric studies in southern Mexico: implications for the generation of volcanism and subduction dynamics

In order to study electrical conductivity phenomena that are associated with subduction related fluid release and melt production, magnetotelluric (MT) measurements were carried out in southern Mexico along two coast to coast profiles. The conductivity-depth distribution was obtained by simultaneous two-dimensional inversion of the transverse magnetic and transverse electric modes of the magnetotelluric transfer functions. The MT models demonstrate that the plate southern profile shows enhanced conductivity in the deep crust. The northern profile is dominated by an elongated conductive zone extending >250 km below the Trans-Mexican Volcanic Belt (TMVB). The isolated conductivity anomalies in the southern profile are interpreted as slab fluids stored in the overlying deep continental crust. These fluids were released by progressive metamorphic dehydration of the basaltic oceanic crust. The conductivity anomalies may be related to the main dehydration reactions at the zeolite → blueschist → eclogite facies transitions and the breakdown of chlorite. This relation allows the estimation of a geothermal gradient of ∼8.5°C/km for the top of the subducting plate. The same dehydration reactions may be recognized along the northern profile at the same position relative to the depth of the plate, but more inland due to a shallower dip, and merge near the volcanic front due to steep downbending of the plate. When the oceanic crust reaches a depth of 80–90 km, ascending fluids produce basaltic melts in the intervening hot subcontinental mantle wedge that give rise to the volcanic belt. Water-rich basalts may intrude into the lower continental crust leading to partial melting. The elongated highly conductive zone below the TMVB may therefore be caused by partial melts and fluids of various origins, ongoing migmatization, ascending basaltic and granitic melts, growing plutons as well as residual metamorphic fluids. Zones of extremely high conductance (>8000 S) in the continental crust on either MT profile might indicate extinct magmatism

Intraplate seismicity and releated mantle hydration at the Nicaraguan trench outer rise

We examine micro-earthquake records from a dense temporary array of ocean bottom seismometers (OBS) and hydrophones that has been installed from September to November 2005 at the trench outer rise offshore Nicaragua. Approximately 1.5 locatable earthquakes per day within the array of 110 × 120 km show the high seismic activity in this region. Seismicity is restricted to the upper ∼15 km of the mantle and hence where temperatures reach 350–400 °C, which is smaller than values observed for large mantle intraplate events (650 °C). Determination of moment tensor solutions suggest a change of the stress region from tensional in the upper layers of the oceanic plate to compressional beneath. The neutral plane between both regimes is located at ∼6–9 km beneath Moho and thus very shallow. Fluids, which are thought to travel through the tensional fault system into the upper mantle, may not be able to penetrate any deeper. The earthquake catalogue, which seems to be complete for magnitudes above Mw = 1.6–1.8, suggests a strong change of the lithospheric rheology when approaching the trench. And b-factors, that is the ratio between small and large earthquakes increase significantly in the closest 20 km to the trench axis, implying that the crust and upper mantle is massively weakened and hence ruptures more frequently but under less release of stress. We explain this with a partly serpentinized upper mantle

Geochemical variations in the Central Southern Volcanic Zone, Chile (38-43°S): The role of fluids in generating arc magmas

We present new Sr-Nd-Pb-Hf-O isotope data from the volcanic arc (VA, volcanic front and rear arc) in Chile and the backarc (BA) in Argentina of the Central Southern Volcanic Zone in Chile (CSVZ; 38-43°S). Compared to the Transitional (T) SVZ (34.5-38°S; Jacques et al., 2013), the CSVZ VA has erupted greater volumes over shorter time intervals (Völker et al., 2011) and produced more tholeiitic melts. Although the CSVZ VA monogenetic cones are similar to the TSVZ VA samples, the CSVZ VA stratovolcanoes have higher ratios of highly fluid-mobile to less fluid-mobile trace elements (e.g. U/Th, Pb/Ce, Ba/Nb) and lower more- to less-incompatible fluid-immobile element ratios (e.g. La/Yb, La/Sm, Th/Yb, Nb/Yb), consistent with an overall higher fluid flux and greater degree of flux melting beneath the CSVZ stratovolcanoes compared to the CSVZ monogenetic centers and the TSVZ VA. The CSVZ monogenetic centers overlap the TSVZ in Sr and Nd isotopes, but the stratovolcanoes are shifted to higher Sr and/or Nd isotope ratios. The Pb isotopic composition of the CSVZ overlaps the TSVZ, which is clearly dominated by the composition of the trench sediments, but the CSVZ monogenetic samples extend to less radiogenic Pb isotope ratios. δ18Omelt from the CSVZ stratovolcano samples are below the MORB range, whereas the CSVZ monogenetic and the TSVZ samples fall within and slightly above the MORB range. The Nd and Hf isotopic ratios of the CSVZ VA extend to more radiogenic compositions than found in the TSVZ VA, indicating a greater contribution from a more depleted source. These correlations are interpreted to reflect derivation of fluids from hydrothermally altered oceanic crust and/or serpentinized upper mantle of the subducting plate. CSVZ BA basalts largely overlap TSVZ BA basalts, displaying less or no subduction influence compared to the VA, but some CSVZ BA basalts tap more enriched mantle, possibly subcontinental lithosphere, with distinctively lower Nd and Hf and elevated 207Pb/204Pb and 208Pb/204Pb isotope ratios

Crustal structure of the Peruvian continental margin from wide-angle seismic studies

Active seismic investigations along the Pacific margin off Peru were carried out using ocean bottom hydrophones and seismometers. The structure and the P-wave velocities of the obliquely subducting oceanic Nazca Plate and overriding South American Plate from 8°S to 15°S were determined by modelling the wide-angle seismic data combined with the analysis of reflection seismic data. Three detailed cross-sections of the subduction zone of the Peruvian margin and one strike-line across the Lima Basin are presented here. The oceanic crust of the Nazca Plate, with a thin pelagic sediment cover, ranging from 0–200 m, has an average thickness of 6.4 km. At 8°S it thins to 4 km in the area of Trujillo Trough, a graben-like structure. Across the margin, the plate boundary can be traced to 25 km depth. As inferred from the velocity models, a frontal prism exists adjacent to the trench axis and is associated with the steep lower slope. Terrigeneous sediments are proposed to be transported downslope due to gravitational forces and comprise the frontal prism, characterized by low seismic P-wave velocities. The lower slope material accretes against a backstop structure, which is defined by higher seismic P-wave velocities, 3.5–6.0 km s−1. The large variations in surface slope along one transect may reflect basal removal of upper plate material, thus steepening the slope surface. Subduction processes along the Peruvian margin are dominated by tectonic erosion indicated by the large margin taper, the shape and bending of the subducting slab, laterally varying slope angles and the material properties of the overriding continental plate. The erosional mechanisms, frontal and basal erosion, result in the steepening of the slope and consequent slope failure

Structure and serpentinization of the subducting Cocos plate offshore Nicaragua and Costa Rica

Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 12 (2011): Q06009, doi:10.1029/2011GC003592.The Cocos plate experiences extensional faulting as it bends into the Middle American Trench (MAT) west of Nicaragua, which may lead to hydration of the subducting mantle. To estimate the along strike variations of volatile input from the Cocos plate into the subduction zone, we gathered marine seismic refraction data with the R/V Marcus Langseth along a 396 km long trench parallel transect offshore of Nicaragua and Costa Rica. Our inversion of crustal and mantle seismic phases shows two notable features in the deep structure of the Cocos plate: (1) Normal oceanic crust of 6 km thickness from the East Pacific Rise (EPR) lies offshore Nicaragua, but offshore central Costa Rica we find oceanic crust from the northern flank of the Cocos Nazca (CN) spreading center with more complex seismic velocity structure and a thickness of 10 km. We attribute the unusual seismic structure offshore Costa Rica to the midplate volcanism in the vicinity of the Galápagos hot spot. (2) A decrease in Cocos plate mantle seismic velocities from ∼7.9 km/s offshore Nicoya Peninsula to ∼6.9 km/s offshore central Nicaragua correlates well with the northward increase in the degree of crustal faulting outboard of the MAT. The negative seismic velocity anomaly reaches a depth of ∼12 km beneath the Moho offshore Nicaragua, which suggests that larger amounts of water are stored deep in the subducting mantle lithosphere than previously thought. If most of the mantle low velocity zone can be interpreted as serpentinization, the amount of water stored in the Cocos plate offshore central Nicaragua may be about 2.5 times larger than offshore Nicoya Peninsula. Hydration of oceanic lithosphere at deep sea trenches may be the most important mechanism for the transfer of aqueous fluids to volcanic arcs and the deeper mantle.This work was funded by the U.S. National Science Foundation MARGINS program under grants OCE0405556, OCE 0405654, and OCE 0625178

Seismic evidence for fluids in fault zones on top of the subducting Cocos Plate beneath Costa Rica

Author Posting. © The Authors, 2010. This article is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 181 (2010): 997-1016, doi:10.1111/j.1365-246X.2010.04552.x.In the 2005 TICOCAVA explosion seismology study in Costa Rica we observed crustal turning waves with a dominant frequency of ~10 Hz on a linear array of short-period seismometers from the Pacific Ocean to the Caribbean Sea. On one of the shot records, from Shot 21 in the backarc of the Cordillera Central, we also observed two seismic phases with an unusually high dominant frequency (~20 Hz). These two phases were recorded in the forearc region of central Costa Rica and arrived ~7 s apart and 30 to 40 s after the detonation of Shot 21. We considered the possibility that these secondary arrivals were produced by a local earthquake that may have happened during the active-source seismic experiment. Such high-frequency phases following Shot 21 were not recorded after Shots 22, 23, and 24, all in the backarc of Costa Rica, which might suggest that they were produced by some other source. However, earthquake dislocation models cannot produce seismic waves of such high frequency with significant amplitude. In addition, we would have expected to see more arrivals from such an earthquake on other seismic stations in central Costa Rica. We therefore investigate whether the high-frequency arrivals may be the result of a deep seismic reflection from the subducting Cocos plate. The timing of these phases is consistent with a shear wave from Shot 21 that was reflected as a compressional (SxP) and a shear (SxS) wave at the top of the subducting Cocos slab between 35 and 55 km depth. The shift in dominant frequency from ~10 Hz in the downgoing seismic wave to ~20 Hz in the reflected waves requires a particular seismic structure at the interface between the subducting slab and the forearc mantle in order to produce a substantial increase in reflection coefficients with frequency. The spectral amplitude characteristics of the SxP and SxS phases from Shot 21 are consistent with a very high Vp/Vs ratio of 6 in ~5 m thick, slab-parallel layers. This result suggests that a system of thin shear zones near the plate interface beneath the forearc is occupied by hydrous fluids under near-lithostatic conditions. The overpressured shear zone probably takes up fluids from the downgoing slab, and it may control the lower limit of the seismogenic zone.This work was funded by the US National Science Foundation MARGINS programme