Lithium isotopic composition of pore fluids and sediments in the Costa Rica subduction zone

Pore fluid and sediment Li concentrations and isotopic ratios provide important insights on the hydrology, sediment contribution to the arc volcanoes and fluid-sediment reactions at the dominantly non-accretionary Costa Rica subduction zone. Ocean Drilling Program Site 1039 in the trench axis provides a reference section of 400 m of the incoming sediments, and Site 1040, situated arcward from the trench, consists of a deformed sedimentary wedge and apron sediments, the décollement, and the partially dewatered underthrust sediment section. At the reference site, pore fluids show important isotopic variations (delta6Li=-21.7 to -37.8 per mil), reflecting the interplay of in situ alteration of volcanic material and ion exchange with clay minerals. In the basal section, a reversal of Li concentration and delta6Li toward seawater values is observed, providing supporting evidence for a lateral seawater flow system in the upper oceanic basement underlying this sediment section. At Site 1040, pore fluid of the lower deformed wedge sediments and within the décollement is enriched in Li and the isotopic compositions are relatively light, suggesting infiltration of a deep-seated fluid. The delta6Li value of -22 per mil of this Li-enriched fluid (261 µM), when compared with the delta6Li value of the subducted sediment section (-11 per mil), suggests that the deep source fluid originates from mineral fluid dehydration and transformation reactions at temperatures of 100 to 150°C, consistent with the temperature range of the up-dip seismogenic zone and of transformation of smectite to illite. The distribution of Li and its isotopes in the underthrust section are similar to those at the reference site, indicating near complete subduction of the incoming sediments and that early dewatering of the underthrust sediments occurs predominantly by lateral flow into the ocean. The hemipelagic clay-rich sediment section of the subducting plate carries most of the Li into this subduction zone, and the pelagic diatomaceous and nannofossil calcareous oozes contain little Li. The Li isotopes of both the clay-rich hemipelagic sediments and of the pelagic oozes are, however, similar, with delta6Li values of -9 to -12 per mil. The observations that (1) the delta6Li values of the underthrust sediments are distinctly lower than that of the mantle, and (2) the lavas of the Costa Rican volcanoes are enriched in Li and 7Li, provide an approximation of the contribution of the subducted sediments to the arc volcanoes. A first order mass balance calculation suggests that approximately half of the Li flux delivered by subducted sediments and altered oceanic crust into the Middle American Trench is recycled to the Costa Rican arc and at most a quarter of sedimentary Li is returned into the ocean through thrust faults, primarily the décollement thrust.

Geochemistry of marine sediments of the world's major subduction zones

This paper presents new major and trace-element data and Lu-Hf and Sm-Nd isotopic compositions for representative suites of marine sediment samples from 14 drill sites outboard of the world's major subduction zones. These suites and samples were chosen to represent the global range in lithology, Lu/Hf ratios, and sediment flux in subducting sediments worldwide. The data reported here represent the most comprehensive data set on subducting sediments and define the Hf-Nd isotopic variations that occur in oceanic sediments and constrain the processes that caused them. Using new marine sediment data presented here, in conjunction with published data, we derive a new Terrestrial Array given by the equation, epsilon-Hf = 1.55 * epsiolon-Nd + 1.21. This array was calculated using >3400 present-day Hf and Nd isotope values. The steeper slope and smaller y-intercept of this array, compared to the original expression (epsilon-Hf = 1.36 * epsilonNd + 2.89; Vervoort et al., 1999, doi:10.1016/S0012-821X(99)00047-3) reflects the use of present day values and the unradiogenic Hf of old continental samples included in the array. In order to examine the Hf-Nd isotopic variations in marine sediments, we have classified our samples into 5 groups based on lithology and major and trace-element geochemical compositions: turbidites, terrigenous clays, and volcaniclastic, hydrothermal and hydrogenetic sediments. Compositions along the Terrestrial Array are largely controlled by terrigenous material derived from the continents and delivered to the ocean basins via turbidites, volcaniclastic sediments, and volcanic inputs from magmatic arcs. Compositions below the Terrestrial Array derive from unradiogenic Hf in zircon-rich turbidites. The anomalous compositions above the Terrestrial Array largely reflect the decoupled behavior of Hf and Nd during continental weathering and delivery to the ocean. Both terrigenous and hydrogenetic clays possess anomalously radiogenic Hf, reflecting terrestrial sedimentary and weathering processes on the one hand and marine inheritance on the other. This probably occurs during complementary processes involving preferential retention of unradiogenic Hf on the continents in the form of zircon and release of radiogenic Hf from the breakdown of easily weathered, high Lu-Hf phases such as apatite.

Partitioning of deformation within a subduction channel during exhumation of high-pressure rocks: a case study from the Western Alps

The metamorphic belt of the Western Alps was subjected to widespread extensional tectonism at the end of the Eocene (ca. 45-35 Ma). Extension was accommodated by hinterland-directed movements along gently inclined extensional shear zones, which facilitated rapid exhumation of high-pressure and ultra-high-pressure rocks. This deformation resulted in a normal metamorphic sequence. Extension in the inner parts of the Western Alps was coeval with shortening at the front of the belt (foreland-directed thrusts), which took place during decompression, and emplaced higher grade metamorphic units over lower grade metamorphic rocks, thus forming an inverse metamorphic sequence. Two mechanisms for this extensional episode are discussed: (1) collapse of an overthickened lithosphere, and (2) internal readjustments within the orogenic wedge due to subduction channel dynamics. We favour the latter mechanism because it can account for the development of the observed inverse and normal metamorphic sequences along foreland-directed thrusts and hinterland-directed detachments, respectively. This hypothesis is supported by published structural, metamorphic and geochronological data from four geological transects through the Western Alps. This study also emphasizes the importance of post-shearing deformation (e.g. horizontal buckling versus vertical flattening), which can modify the distribution of hinterland- and foreland-directed shear zones in orogenic belts. (c) 2006 Elsevier Ltd. All rights reserved.

Geochemical variation among small eruptive centers in the central SVZ of the Andes: An evaluation of subduction, mantle and crustal influences

Subduction zone magmatism is an important and extensively studied topic in igneous geochemistry. Recent studies focus on from where arc magmas are generated, how subduction components (fluids or melts) are fluxed into the source of the magmas, and whether or how the subduction components affect partial melting processes beneath volcanic arcs at convergent boundaries. ^ At 39.5°S in the Central Southern Volcanic Zone of the Andes, Volcano Villarrica is surrounded by a suite of Small Eruptive Centers (SEC). The SECs are located mostly to the east and northeast of the stratovolcano and aligned along the Liquine-Ofqui Fault Zone, the major fracture system in this area. Former studies observed the geochemical patterns of the SECs differ distinctively from those of V. Villarrica and suggested there may be a relationship between the compositions of the volcanic units and their edifice sizes. This work is a comprehensive geochemical study on the SECs near V. Villarrica, using a variety of geochemical tracers and tools including major, trace and REE elements, Li-Be-B elements, Sr-Nd-Pb isotopes and short-lived isotopes such as U-series and 10Be. In this work, systematic differences between the elemental and isotopic compositions of the SECs and those of V. Villarrica are revealed and more importantly, modeled in terms of magmatic processes occurring at continental arc margins. Detailed modeling calculations in this work reconstruct chemical compositions of the primary magmas, source compositions, compositions and percentages of different subduction endmembers mixed into the source, degrees of partial melting and different time scales of the SECs and V. Villarrica, respectively. Geochemical characteristics and possible origins of the two special SECs—andesitic Llizan, with crustal signatures, and Rucapillan, to the northwest toward the trench, are also discussed in this work. ^

Benthic foraminiferal record of the Pleistocene uplift of the sedimentary deposits of the Burica Peninsula (Costa Rica-Panama) as a result of Cocos Ridge subduction beneath the Central American Arc

This study used paleobathymetric interpretations from benthic foraminiferal assemblages to determine the timing of the initial subduction of the Cocos Ridge beneath the Costa Rica-Panama Arc and the associated amount of uplift. Forty-seven samples from four stratigraphic sections of the Burica and the Armuelles formations yielded 217 identified species. Foraminiferal paleoecology and cluster analyses indicated paleobathymetric differences within and between the sections. Maximum and minimum uplift rates were calculated from the biochronology, elevation, paleobathymetry and estimated eustatic changes. Large decreases in paleobathymetry were mainly due to the initial middle Pleistocene subduction of the Cocos Ridge in less than 0.5 Ma. Uplift rates were uneven across the Burica Peninsula, as follows: La Vaca 4.5–0.8 m/ky, Rabo de Puerco 2.8–0.7 m/ky, San Bartolo–Chiquito 2.7–0.8 m/ky, and eastern coast 8.0–1.5 m/ky. These differences probably resulted from Cocos Ridge asymmetry and differential uplift of tectonic blocks.