Dating brittle tectonic movements with cleft monazite: Fluid-rock interaction and formation of REE minerals

[1] Two millimeter-sized hydrothermal monazites from an open fissure (cleft) that developed late during a dextral transpressional deformation event in the Aar Massif, Switzerland, have been investigated using electron microprobe and ion probe. The monazites are characterized by high Th/U ratios typical of other hydrothermal monazites. Deformation events in the area have been subdivided into three phases: (D1) main thrusting including formation of a new schistosity, (D2) dextral transpression, and (D3) local crenulation including development of a new schistosity. The two younger deformational structures are related to a subvertically oriented intermediate stress axis, which is characteristic for strike slip deformation. The inferred stress environment is consistent with observed kinematics and the opening of such clefts. Therefore, the investigated monazite-bearing cleft formed at the end of D2 and/or D3, and during dextral movements along NNW dipping planes. Interaction of cleft-filling hydrothermal fluid with wall rock results in rare earth element (REE) mineral formation and alteration of the wall rock. The main newly formed REE minerals are Y-Si, Y-Nb-Ti minerals, and monazite. Despite these mineralogical changes, the bulk chemistry of the system remains constant and thus these mineralogical changes require redistribution of elements via a fluid over short distances (centimeter). Low-grade alteration enables local redistribution of REE, related to the stability of the accessory phases. This allows high precision isotope dating of cleft monazite. 232Th/208Pb ages are not affected by excess Pb and yield growth domain ages between 8.03 ± 0.22 and 6.25 ± 0.60 Ma. Monazite crystallization in brittle structures is coeval or younger than 8 Ma zircon fission track data and hence occurred below 280°C.

Dating prograde fluid pulses during subduction by in situ U–Pb and oxygen isotope analysis

Keywords High-pressure fluids · Whiteschists · U–Pb dating · Oxygen isotopes · Ion microprobe · Metasomatism Introduction The subduction of crustal material to mantle depths and its chemical modification during burial and exhumation contribute to element recycling in the mantle and the formation of new crust through arc magmatism. Crustal rocks that Abstract The Dora-Maira whiteschists derive from metasomatically altered granites that experienced ultrahighpressure metamorphism at ~750 °C and 40 kbar during the Alpine orogeny. In order to investigate the P–T–time– fluid evolution of the whiteschists, we obtained U–Pb ages from zircon and monazite and combined those with trace element composition and oxygen isotopes of the accessory minerals and coexisting garnet. Zircon cores are the only remnants of the granitic protolith and still preserve a Permian age, magmatic trace element compositions and δ18O of ~10 ‰. Thermodynamic modelling of Si-rich and Si-poor whiteschist compositions shows that there are two main fluid pulses during prograde subduction between 20 and 40 kbar. In Si-poor samples, the breakdown of chlorite to garnet + fluid occurs at ~22 kbar. A first zircon rim directly overgrowing the cores has inclusions of prograde phlogopite and HREE-enriched patterns indicating zircon growth at the onset of garnet formation. A second main fluid pulse is documented close to peak metamorphic conditions in both Si-rich and Si-poor whiteschist when talc + kyanite react to garnet + coesite + fluid. A second metamorphic overgrowth on zircon with HREE depletion was observed in the Si-poor whiteschists, whereas a single metamorphic overgrowth capturing phengite and talc inclusions was observed in the Si-rich whiteschists. Garnet rims, zircon rims and monazite are in chemical and isotopic equilibrium for oxygen, demonstrating that they all formed at peak metamorphism at 35 Ma as constrained by the age of monazite (34.7 ± 0.4 Ma) and zircon rims (35.1 ± 0.8 Ma). The prograde zircon rim in Si-poor whiteschists has an age that is within error indistinguishable from the age of peak metamorphic conditions, consistent with a minimum rate of subduction of 2 cm/year for the Dora-Maira unit. Oxygen isotope values for zircon rims, monazite and garnet are equal within error at 6.4 ± 0.4 ‰, which is in line with closed-system equilibrium fractionation during prograde to peak temperatures. The resulting equilibrium Δ18Ozircon-monazite at 700 ± 20 °C is 0.1 ± 0.7 ‰. The in situ oxygen isotope data argue against an externally derived input of fluids into the whiteschists. Instead, fluidassisted zircon and monazite recrystallisation can be linked to internal dehydration reactions during prograde subduction. We propose that the major metasomatic event affecting the granite protolith was related to hydrothermal seafloor alteration post-dating Jurassic rifting, well before the onset of Alpine subduction.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

The dating of dipterocarp tree rings: establishing a record of carbon cycling and climatic change in the tropics

In a first step to obtain a proxy record of past climatic events (including the El Ni (n) over tildeo-Southern Oscillation) in the normally aseasonal tropical environment of Sabah, a radial segment from a recently fallen dipterocarp (Shorea Superba) was radiocarbon dated and subjected to carbon isotope analysis. The high-precision radiocarbon results fell into the ambiguous modern plateau where several calibrated dates can exist for each sample. Dating was achieved by wiggle matching using a Bayesian approach to calibration. Using the defined growth characteristics of Shorea superba, probability density distributions were calculated and improbable dates rejected. It was found that the tree most likely started growing around AD 1660-1685. A total of 173 apparent growth increments were measured and, therefore, it could be determined that the tree formed one ring approximately every two years. Stable carbon isotope values were obtained from resin-extracted wholewood from each ring. Carbon cycling is evident in the `juvenile effect', resulting from the assimilation of respired carbon dioxide and lower light levels below the canopy, and in the `anthropogenic effect' caused by increased industrial activity in the late-nineteenth and twentieth centuries. This study demonstrates that palaeoenvironmental information can be obtained from trees growing in aseasonal environments, where climatic conditions prevent the formation of well-defined annual rings.

Helium transfer from water into quartz crystals: A new approach for porewater dating

Several important fundamental and applied problems require a quantification of slow rates of groundwater flow. To resolve these problems helium appears to be a promising tracer. In this contribution we discuss a new approach, which gives the helium inventory in a rock – pore water system by using the relevant mineral record, i.e., without extraction and investigation of the porewater samples. Some U- and Th-poor minerals such as quartz (quartz separates from Permo-Carboniferous Formation, sandstone–shale interlayering, Molasses Basin, Northern Switzerland, hereafter PCF, are used in this study) contain excessive helium having migrated into their internal helium-accessible volume (HAV) from the surrounding porewater [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. These volumes are estimated by using helium as a nano-size penetrating tool, i.e., by saturation of the minerals with helium under controlled pressure–temperature conditions and subsequent measurements of the helium-saturated concentrations. In the quartz separates HAV/total volume ratios vary from 0.017% to 0.16%; along with the measured initial (unsaturated) He concentration the HAV gives the internal helium pressure, the mean value obtained for 7 samples (25 sample aliquots) is P=0.45F0.15 atm (1 r). The product of helium pressure and solubility (7.35_10_3 cc STP He/cc H2O for the temperature and salinity of PCF aquifers reported in [F.J. Pearson, W. Balderer, H.H. Loosli, B.E. Lehmann, A. Matter, T. Peters, H. Schmassmann, A. Gautschi, Applied Isotope Hydrogeology–A Case Study in Northern Switzerland, Elsevier Amsterdam, 1991, 439 pp.]) is the mineral-derived He concentration in the respective porewater, CPW=0.0035F0.0017 cc He/cc H2O. This value is in full accord with measured He concentrations in PCF aquifers, CPCF, varying from 0.0045 to 0.0016 cc He/cc H2O. This agreement validates the proposed approach and also shows that the mineral–porewater helium–concentration equilibrium has been established. Indeed, estimates of the He-migration rates through our quartz samples show that in ~6000 years the internal pressure should equilibrate with He-concentration in related porewater of PCF, and this time interval is short compared to independent estimates [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. The helium inventory in the rock–porewater assemblage shows that helium abundance in pore waters is indeed important. In shale samples (with ~15% porosity) porewaters contain more helium than the host minerals altogether. Porewater heliumconcentration profiles, available from the mineral record, along with helium production rates are input parameters allowing model(s) of helium migration through a hydrological structure to be developed. Quite high helium concentrations in PCF porewaters imply slow removal mechanisms, which will be discussed elsewhere.

Investigation of young water inflow in karst aquifers using SF6–CFC–3H/He–85Kr–39Ar and stable isotope components

Karst aquifers are known for their wide distribution of water transfer velocities. From this observation, a multiple geochemical tracer approach seems to be particularly well suited to provide a significant assessment of groundwater flows, but the choice of adapted tracers is essential. In this study, several common tracers in karst aquifers such as physicochemical parameters, major ions, stable isotopes, and d13C to more specific tracers such as dating tracers – 14C, 3H, 3H–3He, CFC-12, SF6 and 85Kr, and 39Ar – were used, in a fractured karstic carbonated aquifer located in Burgundy (France). The information carried by each tracer and the best sampling strategy are compared on the basis of geochemical monitoring done during several recharge events and over longer time periods (months to years). This study’s results demonstrate that at the seasonal and recharge event time scale, the variability of concentrations is low for most tracers due to the broad spectrum of groundwater mixings. The tracers used traditionally for the study of karst aquifers, i.e., physicochemical parameters and major ions, efficiently describe hydrological processes such as the direct and differed recharge, but require being monitored at short time steps during recharge events to be maximized. From stable isotopes, tritium, and Cl� contents, the proportion of the fast direct recharge by the largest porosity was estimated using a binary mixing model. The use of tracers such as CFC-12, SF6, and 85Kr in karst aquifers provides additional information, notably an estimation of apparent age, but they require good preliminary knowledge of the karst system to interpret the results suitably. The CFC-12 and SF6 methods efficiently determine the apparent age of baseflow, but it is preferable to sample the groundwater during the recharge event. Furthermore, these methods are based on different assumptions such as regional enrichment in atmospheric SF6, excess air, and flow models among others. 85Kr and 39Ar concentrations can potentially provide a more direct estimation of groundwater residence time. Conversely, the 3H–3He method is inefficient in the karst aquifer for dating due to 3He degassing.

Groundwater dating with Atom Trap Trace Analysis of 39Ar

We report on the realization of Atom Trap Trace Analysis for39Ar and its first application to dating of groundwater samples. The presented system achieves an atmospheric39Ar count rate as high as 3.58 ± 0.10 atoms/h allowing for the determination of the39Ar concentration in less than a day. We demonstrate that the measured count rates are proportional to the39Ar concentration by intercomparison with Low-Level Counting results and by measurements on prepared argon samples with defined concentration. For a geophysical application, we degas three different groundwater samples and gas chromatographically extract the argon. The39Ar ages inferred from the count rates extend over the accessible dating range and are in agreement with the Low-Level Counting results as well as with complementary isotope data.

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.

Radiometric 81Kr dating identifies 120,000-year-old ice at Taylor Glacier, Antarctica

We present successful 81 Kr-Kr radiometric dating of ancient polarice. Krypton was extracted from the air bubbles in four∼350-kg polar ice samples from Taylor Glacier in the McMurdo Dry Valleys, Antarctica, and dated using Atom Trap Trace Analysis (ATTA). The 81 Kr radiometric ages agree with independent age estimates obtained from stratigraphic dating techniques with a mean abso-lute age offset of 6±2.5 ka. Our experimental methods and sampling strategy are validated by (i) 85 Kr and 39 Ar analyses that show the samples to be free of modern air contamination and (ii)air content measurements that show the ice did not experience gas loss. We estimate the error in the 81 Kr ages due to past geomagnetic variability to be below 3 ka. We show that ice from the previous interglacial period (Marine Isotope Stage 5e, 130–115 ka before present) can be found in abundance near the surface of Taylor Glacier. Our study paves the way for reliable radiometric dating of ancient ice in blue ice areas and margin sites where large samples are available, greatly enhancing their scientific value as archives of old ice and meteorites. At present, ATTA 81Kr analysis requires a 40–80-kg ice sample; as sample requirements continue to decrease, 81 Kr dating of ice cores is a future possibility.

A rapid and efficient ion-exchange chromatography for Lu–Hf, Sm–Nd, and Rb–Sr geochronology and the routine isotope analysis of sub-ng amounts of Hf by MC-ICP-MS

The development and improvement of MC-ICP-MS instruments have fueled the growth of Lu–Hf geochronology over the last two decades, but some limitations remain. Here, we present improvements in chemical separation and mass spectrometry that allow accurate and precise measurements of 176Hf/177Hf and 176Lu/177Hf in high-Lu/Hf samples (e.g., garnet and apatite), as well as for samples containing sub-nanogram quantities of Hf. When such samples are spiked, correcting for the isobaric interference of 176Lu on 176Hf is not always possible if the separation of Lu and Hf is insufficient. To improve the purification of Hf, the high field strength elements (HFSE, including Hf) are first separated from the rare earth elements (REE, including Lu) on a first-stage cation column modified after Patchett and Tatsumoto (Contrib. Mineral. Petrol., 1980, 75, 263–267). Hafnium is further purified on an Ln-Spec column adapted from the procedures of Münker et al. (Geochem., Geophys., Geosyst., 2001, DOI: 10.1029/2001gc000183) and Wimpenny et al. (Anal. Chem., 2013, 85, 11258–11264) typically resulting in Lu/Hf < 0.0001, Zr/Hf < 1, and Ti/Hf < 0.1. In addition, Sm–Nd and Rb–Sr separations can easily be added to the described two-stage ion-exchange procedure for Lu–Hf. The isotopic compositions are measured on a Thermo Scientific Neptune Plus MC-ICP-MS equipped with three 1012 Ω resistors. Multiple 176Hf/177Hf measurements of international reference rocks yield a precision of 5–20 ppm for solutions containing 40 ppb of Hf, and 50–180 ppm for 1 ppb solutions (=0.5 ng sample Hf 0.5 in ml). The routine analysis of sub-ng amounts of Hf will facilitate Lu–Hf dating of low-concentration samples.