Focused fluid transfer through the mantle above subduction zones

Volcanic arcs above subduction zones are enriched in volatiles and fluid-mobile elements with respect to mid-oceanic ridge basalts. There is general consensus that this particular subduction zone signature is generated by fluid-induced extraction of these elements from subducted oceanic crust and its sedimentary cover. However, how these fluids are transferred through the mantle wedge to the locus of partial melting and what modification the fluids will experience is unresolved. Here we investigate the interaction of slab fluids with the mantle wedge through a series of high-pressure experiments. We explore two end-member processes of focused and porous reactive flow of hydrous slab melts through the mantle. Transfer by porous flow leads to the formation of hydrous minerals that sequester fluid-mobile elements and residual fluids characterized by trace element patterns inconsistent with typical arc lavas. In contrast, no hydrous minerals are formed in the reaction zone of experiments mimicking focused flow, and the typical trace element signature acquired during fluid extraction from the slab is preserved, indicating that this is an efficient process for element transfer through the mantle wedge.

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.

Extension and Partitioning in an Oblique Subduction Zone, New Zealand: Constraints from Three-Dimensional Numerical Modeling

Contraction, strike slip, and extension displacements along the Hikurangi margin northeast of the North Island of New Zealand coincide with large lateral gradients in material properties. We use a finite- difference code utilizing elastic and elastic-plastic rheologies to build large- scale, three-dimensional numerical models which investigate the influence of material properties on velocity partitioning within oblique subduction zones. Rheological variation in the oblique models is constrained by seismic velocity and attenuation information available for the Hikurangi margin. We compare the effect of weakly versus strongly coupled subduction interfaces on the development of extension and the partitioning of velocity components for orthogonal and oblique convergence and include the effect of ponded sediments beneath the Raukumara Peninsula. Extension and velocity partitioning occur if the subduction interface is weak, but neither develops if the subduction interface is strong. The simple mechanical model incorporating rheological variation based on seismic observations produces kinematics that closely match those published from the Hikurangi margin. These include extension within the Taupo Volcanic Zone, uplift over ponded sediments, and dextral contraction to the south.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

Heat-flow studies in the Peru Trench subduction zone

New heat-flow values were obtained in the central Peru Trench area during site surveys and drilling of Ocean Drilling Program (ODP) Leg 112 by measuring temperatures with ordinary surface heat-flow probes and in the drill holes and by estimating from bottom-simulating reflectors resulting from gas hydrates. The values determined by these methods are consistent with each other within the limits of error. When combined with existing data, heat-flow distribution from the trench to the coast was delineated. Heat flow is lower than 40 mW/m**2 at the bottom of the trench and 40 to 50 mW/m**2 on the landward slope. The low heat flow at the trench bottom can be explained partly by a high sedimentation rate. Heat flow is variable about where the Mendana Fracture Zone meets the trench. This low heat flow might result from hydrothermal circulation in the fracture zone, which some scientists believe is a new propagating rift. On the landward slope, no significant difference in heat flow is recognized between the northern side and the southern side of the fracture zone, in spite of differences in the age of the subducting plate and the tectonic history. Heat flow on the landward slope may be slightly higher than that in most other subduction zones.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Geodynamic reconstructions of the South America-Antarctica plate system

The South America-Antarctica plate system shows many oceanic accretionary systems and subduction zones that initiated and then stopped. To better apprehend the evolution of the system, geodynamic reconstructions (global) have been created from Jurassic (165 Ma) to present, following the techniques used at the University of Lausanne. However, additional synthetic magnetic anomalies were used to refine the geodynamics between 33 Ma and present. The reconstructions show the break up of Gondwana with oceanisation between South America (SAM) and Antarctica (ANT), together with the break off of `Andean' geodynamical units (GDUs). We propose that oceanisation occurs also east and south of the Scotian GDUs. Andean GDUs collide with other GDUs crossing the Pacific. The west coast of SAM and ANT undergo a subsequent collision with all those GDUs between 103 Ma and 84 Ma, and the Antarctic Peninsula also collides with Tierra del Fuego. The SAM-ANT plate boundary experienced a series of extension and shortening with large strike-slip component, culminating with intra-oceanic subduction leading to the presence of the `V-' and anomalies in the Weddell Sea. From 84 Ma, a transpressive collision takes place in the Scotia region, with active margin to the east. As subduction propagates northwards into an old and dense oceanic crust, slab roll-back initiates, giving rise to the western Scotia Sea and the Powell Basin opening. The Drake Passage opens. As the Scotian GDUs migrate eastwards, there is enough space for them to spread and allow a north-south divergence with a spreading axis acting simultaneously with the western Scotia ridge. Discovery Bank stops the migration of South Orkney and `collides with' the SAM-ANT spreading axis, while the northern Scotian GDUs are blocked against the Falkland Plateau and the North-East Georgia Rise. The western and central Scotia and the Powell Basin spreading axes must cease, and the ridge jumps to create the South Sandwich Islands Sea. The Tierra del Fuego-Patagonia region has always experienced mid-oceanic ridge subduction since 84 Ma. Slab window location is also presented (57-0 Ma), because of its important implication for heat flux and magmatism. (C) 2011 Elsevier Ltd. All rights reserved.