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.

Permeability measurements on sheared and unsheared samples from different IODP Holes of Leg 316 and from IODP Hole 315-C0001E and ODP Hole 190-1174B

At subduction zones, the permeability of major fault zones influences pore pressure generation, controls fluid flow pathways and rates, and affects fault slip behavior and mechanical strength by mediating effective normal stress. Therefore, there is a need for detailed and systematic permeability measurements of natural materials from fault systems, particularly measurements that allow direct comparison between the permeability of sheared and unsheared samples from the same host rock or sediment. We conducted laboratory experiments to compare the permeability of sheared and uniaxially consolidated (unsheared) marine sediments sampled during IODP Expedition 316 and ODP Leg 190 to the Nankai Trough offshore Japan. These samples were retrieved from: (1) The décollement zone and incoming trench fill offshore Shikoku Island (the Muroto transect); (2) Slope sediments sampled offshore SW Honshu (the Kumano transect) ~ 25 km landward of the trench, including material overriden by a major out-of-sequence thrust fault, termed the "megasplay"; and (3) A region of diffuse thrust faulting near the toe of the accretionary prism along the Kumano transect. Our results show that shearing reduces fault-normal permeability by up to 1 order of magnitude, and this reduction is largest for shallow (< 500 mbsf) samples. Shearing-induced permeability reduction is smaller in samples from greater depth, where pre-existing fabric from compaction and lithification may be better developed. Our results indicate that localized shearing in fault zones should result in heterogeneous permeability in the uppermost few kilometers in accretionary prisms, which favors both the trapping of fluids beneath and within major faults, and the channeling of flow parallel to fault structure. These low permeabilities promote the development of elevated pore fluid pressures during accretion and underthrusting, and will also facilitate dynamic hydrologic processes within shear zones including dilatancy hardening and thermal pressurization.

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.