Bias associated with detrital zircon geochronology and thermochronometry

Detrital studies that utilize zircon U–Pb geochronology and fission-track (FT) thermochronometry are subject to a range of potential sources of bias that should be properly evaluated and minimized. Some of them are common to any single-grain mineral analysis (e.g., variable bedrock mineral fertility, hydraulic sorting during transport, selective grain loss during sample processing), whereas others are intrinsic to zircon, and are related to radiation damage and age discordance. In this article, we quantify the impact of intrinsic bias on detrital studies thanks to the analysis of modern detritus shed from the European Alps, and illustrate the general implications on geological interpretations. We show that detrital zircon U–Pb age distributions based on statistically robust datasets are highly reproducible and representative of the parent bedrock ages in the catchment. Arbitrary or selective removal of discordant grain ages can be minimized by using the Kolmogorov–Smirnov test to identify an appropriate cutoff level. Loss of metamict (α-damaged) zircon has a minor impact on data representativeness, and is mainly controlled by regional metamorphism rather than by mechanical abrasion during river transport. Zircon FT grain-age distributions were found to have poor reproducibility, although age spectra are consistent with bedrock data. However, unlike the U–Pb datasets, U-rich zircon grains (> 1000 ppm) are systematically missed, and undatable grains may exceed 50%. We identify two major sources of distribution bias specific to zircon FT datasets: (i) sediment sources dominated by U-rich zircon grains are markedly underrepresented in the detrital record, because such grains often have uncountable high densities of fission tracks (“U concentration bias”); (ii) sediment sources that shed zircon grains with high levels of α-damage are underrepresented, because these grains are lost during chemical etching for FT revelation (“etching bias”). In the case of multimethod dating on the same grains (e.g., FT and U–Pb double dating), bias affecting detrital zircon FT dating propagates to the entire dataset. These effects may not impact on exhumation-rate studies that utilize the youngest grain ages (i.e., lag-time approach). However, they represent a limiting factor for conventional provenance studies, and generally preclude application of zircon FT dating to sediment budget calculations.

Dating polygenetic metamorphic assemblages along a transect across the Western Alps

Multichronometric analyses were performed on samples from a transect in the French-Italian Western Alps crossing nappes derived from the Briançonnais terrane and the Piemonte-Liguria Ocean, in an endeavour to constrain the high-pressure (HP) metamorphism and the retrogression history. 12 samples of white mica were analysed by 39Ar-40Ar stepwise heating, complemented by 2 samples from the Monte Rosa 100 km to the NE and also attributed to the Briançonnais terrane. One Sm-Nd and three Lu-Hf garnet ages from eclogites were also obtained. White mica ages decrease from ca. 300 Ma in the westernmost samples (Zone Houillère), reaching ca. 300 °C during Alpine metamorphism, to < 48 Ma in the internal units to the East, which reached ca. 500 °C during Alpine orogeny. The conventional “thermochronological” interpretation postulates Cretaceous Eo-Alpine HP metamorphism and younger “cooling ages” in the higher-temperature samples. However, Eocene Lu-Hf and Sm-Nd ages from the same samples cannot be interpreted as post-metamorphic cooling ages, which makes a Cretaceous eclogitization untenable. The age date from this transect require instead to replace conventional “thermochronology” by an approach combining age dating with detailed geochemical, petrological and microstructural investigations. Petrology reveals important mineralogical differences along the transect. Samples from the Zone Houillère mostly contain detrital mica. White mica with Si > 6.45 atoms per formula unit becomes more abundant eastward. Across the whole traverse, HP phengitic mica forms the D1 foliation. Syn-D2 mica is Si-poorer and associated with nappe stacking, exhumation, and hydrous retrogression under greenschist facies conditions. D1 phengite is very often corroded, overgrown or intergrown by syn-D2 muscovite. Most importantly, syn-D2 recrystallization is not limited to S2 schistosity domains; microchemical fingerprinting shows that it also can form pseudomorphs after crystals that could be mistaken to have formed during D1 based on microstructural arguments alone. Thereby the Cl concentration in white mica is a useful discriminator, since D2 retrogression was associated with a less saline fluid than eclogitization. Once the petrological stage is set, geochronology is straightforward. All samples contain mixtures of detrital, syn-D1 and syn-D2 mica, and retrogression phases (D3) in greatly varying proportions according to local pressure-temperature-fluid activity-deformation conditions. The correlation of age vs. Cl/K clearly identifies 47 ± 1 Ma as the age of formation of syn-D1 mica along the entire transect, including the Monte Rosa nappe samples. The inferred age of the greenschist-facies low-Si syn-D2 mica generation ranges within 39-43 Ma, with local variations. Coexistence of D1 and D2 ages, and the constancy of non-reset D1 ages along the entire transect, are strong evidence that the D1 white mica ages are very close to formation ages. Volume diffusion of Ar in white mica (activation energy E = 250 kJ/mol; pressure-adjusted diffusion coefficient D’0 < 0.03 cm2 s-1) has a subordinate effect on mineral ages compared to both prograde and retrograde recrystallization in most samples. Eocene Lu-Hf and Sm-Nd garnet ages are prograde and predate the HP peak.

Precambrian to Paleozoic zircon record in the Siviez-Mischabel basement (western Swiss Alps)

U–Pb zircon analyses from three meta-igneous and two metasedimentary rocks from the Siviez-Mischabel nappe in the western Swiss Alps are presented, and are used to derive an evolutionary history spanning from Paleoarchean crustal growth to Permian magmatism. The oldest components are preserved in zircons from metasedimentary albitic schists. The oldest zircon core in these schists is 3.4 Ga old. Detrital zircons reveal episodes of crustal growth in the Neoarchean (2.7–2.5 Ga), Paleoproterozoic (2.2–1.9 Ma) and Neoproterozoic (800–550 Ma, Pan-African event). The maximum age of deposition for the metasedimentary rocks is given by the youngest detrital zircons within both metasedimentary samples dated at ~490 Ma (Cambrian-Ordovician boundary). This is in the age range of two granitoid samples dated at 505 ± 4 and 482 ± 7 Ma, and indicates sedimentation and magmatism in an extensional setting preceding an Ordovician orogeny. The third felsic meta-igneous rock gives a Permian age of intrusion, and is part of a long-lasting Variscan to post-Variscan magmatic activity. The zircons record only minor disturbance of the U–Pb system during the Alpine orogeny.

Age of cleft monazites in the eastern Tauern Window: constraints on crystallization conditions of hydrothermal monazite

Monazite-bearing Alpine clefts located in the Sonnblick region of the eastern Tauern Window, Austria, are oriented perpendicular to the foliation and lineation. Ion probe (SIMS) Th–Pb and U–Pb dating of four cleft monazites yields crystallization ages of different growth domains and aggregate regions ranging from 18.99 ± 0.51 to 15.00 ± 0.51 Ma. The crystallization ages obtained are overlapping or slightly younger than zircon fission track ages but older than zircon (U–Th)/He cooling ages from the same area. This constrains cleft monazite crystallization in this area to *300–200 �C. LA-ICP-MS data of dated hydrothermal monazites indicate that in graphite-bearing, reduced host lithologies, cleft monazite is poor in As and has higher La/Yb values and U concentrations, whereas in oxidised host rocks opposite trends are observed. Monazites show negative Eu anomalies and variable La/Yb values ranging from 520 to 6050. The positive correlation between Ca and Sr concentration indicates dissolution of plagioclase or carbonates as the source of these elements. The data show that early exhumation and cleft formation in the Tauern is related to metamorphic dome formation caused by the collision of the Adriatic with the European plate and that monazite crystallization in the clefts occurred later. Our data also demonstrate that hydrothermal monazite ages offer great potential in helping to constrain the chronology of exhumation in collisional orogens.

Structural, metamorphic and geochronological relations between the Zanskar Shear Zone and the Miyar Shear Zone (NW Indian Himalaya): Evidence for two distinct tectonic structures and implications for the evolution of the High Himalayan Crystalline of Zanskar

The geologic structures and metamorphic zonation of the northwestern Indian Himalaya contrast significantly with those in the central and eastern parts of the range, where the high-grade metamorphic rocks of the High Himalayan Crystalline (HHC) thrust southward over the weakly metamorphosed sediments of the Lesser Himalaya along the Main Central Thrust (MCT). Indeed, the hanging wall of the MCT in the NW Himalaya mainly consists of the greenschist facies metasediments of the Chamba zone, whereas HHC high-grade rocks are exposed more internally in the range as a large-scale dome called the Gianbul dome. This Gianbul dome is bounded by two oppositely directed shear zones, the NE-dipping Zanskar Shear Zone (ZSZ) on the northern flank and the SW-dipping Miyar Shear Zone (MSZ) on the southern limb. Current models for the emplacement of the HHC in NW India as a dome structure differ mainly in terms of the roles played by both the ZSZ and the MSZ during the tectonothermal evolution of the HHC. In both the channel flow model and wedge extrusion model, the ZSZ acts as a backstop normal fault along which the high-grade metamorphic rocks of the HHC of Zanskar are exhumed. In contrast, the recently proposed tectonic wedging model argues that the ZSZ and the MSZ correspond to one single detachment system that operates as a subhorizontal backthrust off of the MCT. Thus, the kinematic evolution of the two shear zones, the ZSZ and the MSZ, and their structural, metamorphic and chronological relations appear to be diagnostic features for discriminating the different models. In this paper, structural, metamorphic and geochronological data demonstrate that the MSZ and the ZSZ experienced two distinct kinematic evolutions. As such, the data presented in this paper rule out the hypothesis that the MSZ and the ZSZ constitute one single detachment system, as postulated by the tectonic wedging model. Structural, metamorphic and geochronological data are used to present an alternative tectonic model for the large-scale doming in the NW Indian Himalaya involving early NE-directed tectonics, weakness in the upper crust, reduced erosion at the orogenic front and rapid exhumation along both the ZSZ and the MSZ.