Geochronology of metasomatic events

In order to date any geological event, suitable mineral geochronometers that record that and only that event must be identified and analyzed. In the case of metasomatism, recrystallisation is a key process that controls both the petrology and the isotopic record of minerals. It can occur both in the form of complete neocrystallisation (e.g. in a vein) and in the form of pseudomorphism, whereby dissolution/reprecipitation at the submicroscopic scale plays a central role. Recrystallisation may be complete or not, raising the possibility that relicts of a pre-metasomatic assemblage may be preserved. Because recrystallisation is energetically less costly at almost any temperature than diffusion, and because radiogenic isotopes (except 4He) never diffuse faster than major elements forming the mineral structure, there is a strong causal link between petrographic relicts and isotopic inheritance (as demonstrated for zircon, monazite, titanite, amphibole, K-feldspar, biotite, and muscovite). Metasomatic assemblages commonly contain such mixtures between relicts and newly formed phases, whose geochronology is slightly more complex than that of simple, ideal systems, but can be managed by techniques that have become routine in the last decade and which are described in this chapter. Because of its crucial role in controlling the isotope systematics, the petrogenesis of a mineral needs to be understood in extreme detail, especially using microchemical analyses and micro-imaging techniques, before mineral ages can be correctly interpreted. As the occurrence of recrystallization is limited by the availability of water, minerals act as “geohygrometers” that allow constraints to be placed on the nature and age of fluid circulation episodes, especially metasomatic events.

Ar diffusion in K-feldspar: present and absent

The mechanisms of Ar release from K-feldspar samples in laboratory experiments and during their geological history are assessed here. Modern petrology clearly established that the chemical and isotopic record of minerals is normally dominated by aqueous recrystallization. The laboratory critique is trickier, which explains why so many conflicting approaches have been able to survive long past their expiration date. Current models are evaluated for self-consistency; especially Arrhenian non-linearity leads to paradoxes. The models’ testable geological predictions suggest that temperature-based downslope extrapolations often overestimate observed geological Ar mobility substantially. An updated interpretation is based on the unrelatedness of geological behaviour to laboratory experiments. The isotopic record of K-feldspar in geological samples is not a unique function of temperature, as recrystallisation promoted by aqueous fluids is the predominant mechanism controlling isotope transport. K-feldspar should therefore be viewed as a hygrochronometer. Laboratory degassing proceeds from structural rearrangements and phase transitions such as are observed in situ at high temperature in Na and Pb feldspars. These effects violate the mathematics of an inert Fick’s Law matrix and preclude downslope extrapolation. The similar upward-concave, non-linear shapes of Arrhenius trajectories of many silicates, hydrous and anhydrous, are likely common manifestations of structural rearrangements in silicate structures.

K-feldspar hygrochronology

The significance of the multi-isotopic record preserved in K-feldspars is assessed on samples from the Aar metagranite, Central Alps, Switzerland having very tight independent geological constraints. Stepwise leaching reveals that two diachronically grown K-feldspar generations coexist: Kfs-1 (≥ 35 Ma old, Ca-poor, Rb-Cl-rich, with low 87Sr/86Sr and high 206Pb/204Pb) and Kfs-2 (≤ 10 Ma old, antithetic isotopic signatures deriving from external fluids). Microtextures imaged by cathodoluminescence, backscattered electrons, and electron probe microanalysis are patchy and chemically heterogeneous, with pronounced enrichments in Ba in the retrogressed regions. This confirms the simultaneous presence of fluid-dominated retrogression and recrystallization and isotopic inheritance. The staircase-shaped 40Ar/39Ar age spectrum correlates with the Ca/K and Cl/K signatures. This reflects a mixture of heterochemical K-feldspar generations, and not an intracrystalline Ar gradient caused by diffusion. The shape of the age spectrum and the in vacuo release kinetics proceed from entirely different physical and geological phenomena. What K-feldspars can be effectively used for is to constrain the timing of the fluids that interacted with them by multi-isotopic analyses, rather than to model a “cooling history” from 39Ar release alone. The identification of multiple mineral generations by imaging combined with multi-isotopic analysis enables the accurate dating of the events of a multistage evolution after the initial crystallization of the rock in which the minerals occur.

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

Independent variations of CH₄ emissions and isotopic composition over the past 160,000 years

During the last glacial cycle, greenhouse gas concentrations fluctuated on decadal and longer timescales. Concentrations of methane, as measured in polar ice cores, show a close connection with Northern Hemisphere temperature variability, but the contribution of the various methane sources and sinks to changes in concentration is still a matter of debate. Here we assess changes in methane cycling over the past 160,000 years by measurements of the carbon isotopic composition delta C-13 of methane in Antarctic ice cores from Dronning Maud Land and Vostok. We find that variations in the delta C-13 of methane are not generally correlated with changes in atmospheric methane concentration, but instead more closely correlated to atmospheric CO2 concentrations. We interpret this to reflect a climatic and CO2-related control on the isotopic signature of methane source material, such as ecosystem shifts in the seasonally inundated tropical wetlands that produce methane. In contrast, relatively stable delta C-13 values occurred during intervals of large changes in the atmospheric loading of methane. We suggest that most methane sources-most notably tropical wetlands-must have responded simultaneously to climate changes across these periods.

A new methodological protocol for the use of dendrogeomorphological data in flood risk analysis

Dendrogeomorphology uses information sources recorded in the roots, trunks and branches of trees and bushes located in the fluvial system to complement (or sometimes even replace) systematic and palaeohydrological records of past floods. The application of dendrogeomorphic data sources and methods to palaeoflood analysis over nearly 40 years has allowed improvements to be made in frequency and magnitude estimations of past floods. Nevertheless, research carried out so far has shown that the dendrogeomorphic indicators traditionally used (mainly scar evidence), and their use to infer frequency and magnitude, have been restricted to a small, limited set of applications. New possibilities with enormous potential remain unexplored. New insights in future research of palaeoflood frequency and magnitude using dendrogeomorphic data sources should: (1) test the application of isotopic indicators (16O/18O ratio) to discover the meteorological origin of past floods; (2) use different dendrogeomorphic indicators to estimate peak flows with 2D (and 3D) hydraulic models and study how they relate to other palaeostage indicators; (3) investigate improved calibration of 2D hydraulic model parameters (roughness); and (4) apply statistics-based cost–benefit analysis to select optimal mitigation measures. This paper presents an overview of these innovative methodologies, with a focus on their capabilities and limitations in the reconstruction of recent floods and palaeofloods.

The abundance and isotopic composition of Cd in iron meteorites

Cadmium is a highly volatile element and its abundance in meteorites may help better understand volatility-controlled processes in the solar nebula and on meteorite parent bodies. The large thermal neutron capture cross section of 113Cd suggests that Cd isotopes might be well suited to quantify neutron fluences in extraterrestrial materials. The aims of this study were (1) to evaluate the range and magnitude of Cd concentrations in magmatic iron meteorites, and (2) to assess the potential of Cd isotopes as a neutron dosimeter for iron meteorites. Our new Cd concentration data determined by isotope dilution demonstrate that Cd concentrations in iron meteorites are significantly lower than in some previous studies. In contrast to large systematic variations in the concentration of moderately volatile elements like Ga and Ge, there is neither systematic variation in Cd concentration amongst troilites, nor amongst metal phases of different iron meteorite groups. Instead, Cd is strongly depleted in all iron meteorite groups, implying that the parent bodies accreted well above the condensation temperature of Cd (i.e., ≈650 K) and thus incorporated only minimal amounts of highly volatile elements. No Cd isotope anomalies were found, whereas Pt and W isotope anomalies for the same iron meteorite samples indicate a significant fluence of epithermal and higher energetic neutrons. This observation demonstrates that owing to the high Fe concentrations in iron meteorites, neutron capture mainly occurs at epithermal and higher energies. The combined Cd-Pt-W isotope results from this study thus demonstrate that the relative magnitude of neutron capture-induced isotope anomalies is strongly affected by the chemical composition of the irradiated material. The resulting low fluence of thermal neutrons in iron meteorites and their very low Cd concentrations make Cd isotopes unsuitable as a neutron dosimeter for iron meteorites.

Two-phase change in CO2, Antarctic temperature and global climate during Termination II

The end of the Last Glacial Maximum (Termination I), roughly 20 thousand years ago (ka), was marked by cooling in the Northern Hemisphere, a weakening of the Asian monsoon, a rise in atmospheric CO2 concentrations and warming over Antarctica. The sequence of events associated with the previous glacial–interglacial transition (Termination II), roughly 136 ka, is less well constrained. Here we present high-resolution records of atmospheric CO2 concentrations and isotopic composition of N2—an atmospheric temperature proxy—from air bubbles in the EPICA Dome C ice core that span Termination II. We find that atmospheric CO2 concentrations and Antarctic temperature started increasing in phase around 136 ka, but in a second phase of Termination II, from 130.5 to 129 ka, the rise in atmospheric CO2 concentrations lagged that of Antarctic temperature unequivocally. We suggest that during this second phase, the intensification of the low-latitude hydrological cycle resulted in the development of a CO2 sink, which counteracted the CO2 outgassing from the Southern Hemisphere oceans over this period.

Comment on “The phase relation between atmospheric carbon dioxide and global temperature” Humlum et al. [Glob. Planet. Change 100: 51–69.]: Isotopes ignored

A recent study relying purely on statistical analysis of relatively short time series suggested substantial re-thinking of the traditional view about causality explaining the detected rising trend of atmospheric CO2 (atmCO2) concentrations. If these results are well-justified then they should surely compel a fundamental scientific shift in paradigms regarding both atmospheric greenhouse warming mechanism and global carbon cycle. However, the presented work suffers from serious logical deficiencies such as, 1) what could be the sink for fossil fuel CO2 emissions, if neither the atmosphere nor the ocean – as suggested by the authors – plays a role? 2) What is the alternative explanation for ocean acidification if the ocean is a net source of CO2 to the atmosphere? Probably the most provocative point of the commented study is that anthropogenic emissions have little influence on atmCO2 concentrations. The authors have obviously ignored the reconstructed and directly measured carbon isotopic trends of atmCO2 (both δ13C, and radiocarbon dilution) and the declining O2/N2 ratio, although these parameters provide solid evidence that fossil fuel combustion is the major source of atmCO2 increase throughout the Industrial Era.

Can the carbon isotopic composition of methane be reconstructed from multi-site firn air measurements?

Methane is a strong greenhouse gas and large uncertainties exist concerning the future evolution of its atmospheric abundance. Analyzing methane atmospheric mixing and stable isotope ratios in air trapped in polar ice sheets helps in reconstructing the evolution of its sources and sinks in the past. This is important to improve predictions of atmospheric CH4 mixing ratios in the future under the influence of a changing climate. The aim of this study is to assess whether past atmospheric δ13C(CH4) variations can be reliably reconstructed from firn air measurements. Isotope reconstructions obtained with a state of the art firn model from different individual sites show unexpectedly large discrepancies and are mutually inconsistent. We show that small changes in the diffusivity profiles at individual sites lead to strong differences in the firn fractionation, which can explain a large part of these discrepancies. Using slightly modified diffusivities for some sites, and neglecting samples for which the firn fractionation signals are strongest, a combined multi-site inversion can be performed, which returns an isotope reconstruction that is consistent with firn data. However, the isotope trends are lower than what has been concluded from Southern Hemisphere (SH) archived air samples and high-accumulation ice core data. We conclude that with the current datasets and understanding of firn air transport, a high precision reconstruction of δ13C of CH4 from firn air samples is not possible, because reconstructed atmospheric trends over the last 50 yr of 0.3–1.5 ‰ are of the same magnitude as inherent uncertainties in the method, which are the firn fractionation correction (up to ~2 ‰ at individual sites), the Kr isobaric interference (up to ~0.8 ‰, system dependent), inter-laboratory calibration offsets (~0.2 ‰) and uncertainties in past CH4 levels (~0.5 ‰).

Origin and evolution of reactive and noble gases dissolved in matrix pore water

Reactive and noble gases dissolved in matrix pore water of low permeable crystalline bedrock were successfully extracted and characterized for the fist time based on drillcore samples from the Olkiluoto investigation site (SW Finland). Interaction between matrix pore water and fracture groundwater occurs predominately by diffusion. Changes in the chemical and isotopic composition of gases dissolved in fracture groundwater are transmitted and preserved in the pore water. Absolute concentrations, their ratios and the stable carbon isotope signature of hydrocarbon gases dissolved in pore water give valuable indications about the evolution of these gases in the nearby-flowing fracture groundwaters. Inert noble gases dissolved in matrix pore water and their isotopes combined with their in-situ production and accumulation rates deliver information about the residence time of pore water.

Stalagmite water content as a proxy for drip water supply in tropical and subtropical areas

In this pilot study water was extracted from samples of two Holocene stalagmites from Socotra Island, Yemen, and one Eemian stalagmite from southern continental Yemen. The amount of water extracted per unit mass of stalagmite rock, termed "water yield" hereafter, serves as a measure of its total water content. Based on direct correlation plots of water yields and δ18Ocalcite and on regime shift analyses, we demonstrate that for the studied stalagmites the water yield records vary systematically with the corresponding oxygen isotopic compositions of the calcite (δ18Ocalcite). Within each stalagmite lower δ18Ocalcite values are accompanied by lower water yields and vice versa. The δ18Ocalcite records of the studied stalagmites have previously been interpreted to predominantly reflect the amount of rainfall in the area; thus, water yields can be linked to drip water supply. Higher, and therefore more continuous drip water supply caused by higher rainfall rates, supports homogeneous deposition of calcite with low porosity and therefore a small fraction of water-filled inclusions, resulting in low water yields of the respective samples. A reduction of drip water supply fosters irregular growth of calcite with higher porosity, leading to an increase of the fraction of water-filled inclusions and thus higher water yields. The results are consistent with the literature on stalagmite growth and supported by optical inspection of thin sections of our samples. We propose that for a stalagmite from a dry tropical or subtropical area, its water yield record represents a novel paleo-climate proxy recording changes in drip water supply, which can in turn be interpreted in terms of associated rainfall rates.

Stable-Isotopic Composition of Precipitation Over the Northern Slope of the Central Himalaya

Stable-water-isotope data (deltaD and delta(18)O) from three groups of samples (fresh-snow and snow-pit samples collected on Qomolangma (Mount Everest) and Xixabangma during field seasons 1997,1998 and 2001, and precipitation samples collected at Tingri station during summer 2000) are presented and used to survey the isotopic composition of precipitation over the northern slope of the central Himalaya. Multi-year snow-pit samples on Qomolangma have a local meteoric water-line (slope = 8) close to the global value. Deuterium excess (d = deltaD - 8delta(18)O) values at Tingri are much lower than those in fresh snow from Qomolangma, probably due to differences in moisture source and air-mass trajectories as well as local weather conditions. There is no obvious seasonal trend for d values in the Qomolangma region. A negative relationship exists between delta(18)O and d values in both fresh snow on Qomolangma and precipitation at Tingri. Fresh-snow samples collected from different altitudes on Xixabangma allow us to investigate the altitude effect on delta(18)O values in snow. Of four storm events, only one has an obvious altitude effect on delta(18)O variation and a very low gradient of -0.1% per 100 in elevation.

Oxygen Isotopic and Soluble Ionic Composition of a Shallow Firn Core, Inilchek Glacier, Central Tien Shan

Oxygen isotopic and soluble ionic measurements made on snow-pit (2 in depth) and firn-core (12.4 m depth samples recovered from the accumulation zone 5100 m) of Inilchek glacier 43degrees N, 79degrees E) provide information on recent (1992-98) climatic and environmental conditions in the central Tien Shan region of central Asia. The combined 14.4 m snow-pit/firn-core profile lies within the firn zone, arid contains only one observed melt feature (10 m temperature = - 12 degreesC), Although some post-depositional attenuation of the sub-seasonal delta(18)O record is possible, annual cycles are apparent throughout the isotope profile. We therefore use the preserved delta(18)O record to establish a depth/age scale for the core. Mean delta(18)O values for the entire core and for summer periods are consistent with delta(18)O/temperature observations, and suggest the delta(18)O record provides a means to reconstruct past changes in summer surface temperature at the site. Major-ion (Na(+), K(+), Mg(2+), Ca(2+), NH(4)(+), Cl(-), NO(3)(-), SO(4)(2-)) data from the core demonstrate the dominant influence of dust deposition on the soluble chemistry at the site, arid indicate significant interannual variability in atmospheric-dust loading during the 1900s. Anthropogenic impacts oil NH(4)(+) concentrations are observed at the site, and suggest a summer increase in atmospheric NH(4)(+) that may be related to regional agricultural (nitrogen-rich fertilizer use activities.

Sulfur Isotopic Measurements from a West Antarctic Ice Core: Implications for Sulfate Source and Transport

Measurements of delta(34)S covering the years 1935-76 and including the 1963 Agung (Indonesia) eruption were made on a West Antarctic firn core, RIDSA (78.73 degrees S, 116.33 degrees W; 1740m a.s.l.), and results are used to unravel potential source functions in the sulfur cycle over West Antarctica. The delta(34)S values Of SO42- range from 3.1 parts per thousand to 9.9 parts per thousand. These values are lower than those reported for central Antarctica, from near South Pole station, of 9.3-18.1 parts per thousand (Patris and others, 2000). While the Agung period is isotopically distinct at South Pole, it is not in the RIDSA dataset, suggesting differences in the source associations for the sulfur cycle between these two regions. Given the relatively large input of marine aerosols at RIDSA (determined from Na+ data and the seasonal SO42- cycle), there is likely a large marine biogenic SO42- influence. The delta(34)S values indicate, however, that this marine biogenic SO42-, with a well-established delta(34)S of 18 parts per thousand, is mixing with SO42- that has extremely negative delta(34)S values to produce the measured isotope values in the RIDSA core. We suggest that the transport and deposition of stratospheric SO42- in West Antarctica, combined with local volcanic input, accounts for the observed variance in delta(34)S values.