Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

Antarctic temperature changes during the last millennium: evaluation of simulations and reconstructions

Temperature changes in Antarctica over the last millennium are investigated using proxy records, a set of simulations driven by natural and anthropogenic forcings and one simulation with data assimilation. Over Antarctica, a long term cooling trend in annual mean is simulated during the period 1000–1850. The main contributor to this cooling trend is the volcanic forcing, astronomical forcing playing a dominant role at seasonal timescale. Since 1850, all the models produce an Antarctic warming in response to the increase in greenhouse gas concentrations. We present a composite of Antarctic temperature, calculated by averaging seven temperature records derived from isotope measurements in ice cores. This simple approach is supported by the coherency displayed between model results at these data grid points and Antarctic mean temperature. The composite shows a weak multi-centennial cooling trend during the pre-industrial period and a warming after 1850 that is broadly consistent with model results. In both data and simulations, large regional variations are superimposed on this common signal, at decadal to centennial timescales. The model results appear spatially more consistent than ice core records. We conclude that more records are needed to resolve the complex spatial distribution of Antarctic temperature variations during the last millennium.

Qualitative Distinction of Autotrophic and Heterotrophic Processes at the Leaf Level by Means of Triple Stable Isotope (C–O–H) Patterns

Foliar samples were harvested from two oaks, a beech, and a yew at the same site in order to trace the development of the leaves over an entire vegetation season. Cellulose yield and stable isotopic compositions (δ13C, δ18O, and δD) were analyzed on leaf cellulose. All parameters unequivocally define a juvenile and a mature period in the foliar expansion of each species. The accompanying shifts of the δ13C-values are in agreement with the transition from remobilized carbohydrates (juvenile period), to current photosynthates (mature phase). While the opponent seasonal trends of δ18O of blade and vein cellulose are in perfect agreement with the state-of-art mechanistic understanding, the lack of this discrepancy for δD, documented for the first time, is unexpected. For example, the offset range of 18 permil (oak veins) to 57 permil (oak blades) in δD may represent a process driven shift from autotrophic to heterotrophic processes. The shared pattern between blade and vein found for both oak and beech suggests an overwhelming metabolic isotope effect on δD that might be accompanied by proton transfer linked to the Calvin-cycle. These results provide strong evidence that hydrogen and oxygen are under different biochemical controls even at the leaf level.

39Ar-40Ar geochronology of mono- and polymetamorphic basement rocks

The 39Ar-40Ar technique is often used to date the metamorphic evolution of basement rocks. The present review article examines systematic aspects of the K-Ar decay system in different mineral chronometers frequently found in mono- and polymetamorphic basements (amphibole, biotite, muscovite/phengite, K-feldspar). A key observation is that the measured dissolution rate of silicates in aqueous fluids is many orders of magnitude faster, and has a much lower activation energy, than the rate of Fickian diffusion of Ar. The effects of this inequality are patchy age zonations, very much like those observed in many U-Pb chronometers, unaccompanied by intra-crystalline bell¬shaped Ar loss profiles. Recognizing the importance of the respective rate constants in field situations leads to re-evaluating the ages and the interpretive paradigms in classic examples such as the Central Alpine "Lepontine" amphibolite event and the Western Alpine eclogitic event.

Apatite as an indicator of fluid salinity: An experimental study of chlorine and fluorine partitioning in subducted sediments

In order to constrain the salinity of subduction zone fluids, piston-cylinder experiments have been conducted to investigate the partitioning behaviour of Cl and F in subducted sediments. These experiments were performed at H2O-undersaturated conditions with a synthetic pelite starting composition containing 800 ppm Cl, over a pressure and temperature range of 2.5–4.5 GPa and 630–900 °C. Repetitive experiments were conducted with 1900 ppm Cl + 1000 ppm F, and 2100 ppm Cl. Apatite represents the most Cl-abundant mineral phase, with Cl concentration varying in the range 0.1–2.82 wt%. Affinity for Cl decreases over the following sequence: aqueous fluid > apatite ⩾ melt > other hydrous minerals (phengite, biotite and amphibole). It was found that addition of F to the Cl-bearing starting composition significantly lowers the Cl partition coefficients between apatite and melt (DClAp–melt) and apatite and aqueous fluid (DClAp–aq). Cl–OH exchange coefficients between apatite and melt (KdCl–OHAp–melt) and apatite and aqueous fluid (KdCl–OHAp–aq) were subsequently calculated. KdCl–OHAp–melt was found to vary from 1 to 58, showing an increase with temperature and a decrease with pressure and displaying a regular decrease with increasing H2O content in melt. Mole fractions of Cl and OH in melt were calculated based on an ideal mixing model for H2O, OH, O, Cl and F. The Cl contents of other hydrous minerals (phengite, biotite and amphibole) fall between 200 and 800 ppm, with resultant Cl partition coefficients from 0.02 to 0.49, appearing independent of the bulk Cl and F content. Preliminary data from this study show that the partitioning behaviour of F is strongly in favour of apatite relative to melt and phengite, with DFAp–melt = 15–51. Apatites from representative eclogite facies metasediments were examined and found to have low Cl contents close to ∼100 ppm. Calculations using our experimentally determined KdCl–OHAp–aq of 0.004 at 2.5 GPa, 630 °C indicate a low salinity character (0.5–2 wt% NaCleq) for the fluid formed during dehydration of subducted oceanic sediment at ∼80 km depth.

On the theory of high and low pressure areas: The significance of divergence in pressure areas

This is the edited and translated version of the article by Richard Scherhag “Zur Theorie der Hoch- und Tiefdruckgebiete. Die Bedeutung der Divergenz in Druckfeldern” (On the theory of high and low pressure areas: The significance of divergence in pressure areas), which was published in Meteorologische Zeitschrift 51, 129–138.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.