Petrology and geochronology of "muscovite age standard" B4M

Muscovite B4M, distributed in 1961 as an age standard, was ground under ethanol. Five grain size fractions were obtained and characterized by X-ray diffraction. They display a mixing trend between a phengitic (enriched in the fraction <0.2 µm) and a muscovitic component (predominant in the fraction >20 µm). High-pressure phengite is preserved as a relict in retrograde muscovite. Electron microprobe analyses of the distributed mineral separate reveal at least four white mica populations based on Si, Al, Mg, Na, Fe and F. Rb/K ratios vary by one order of magnitude. Rb–Sr analyses link the mineralogical heterogeneity to variable Rb/Sr and 87Sr/86Sr ratios. The grain size fractions define no internal isochron. Relict fine-grained phengite gives older ages than coarse-grained retrograde greenschist facies muscovite. The inverse grain size–age relationship also characterizes 39Ar/40Ar analyses. Cl/K anticorrelates with step ages: Cl-rich coarse muscovite is younger than Cl-poor fine relict phengite. Sr and Ar preserve a similar isotopic inheritance despite peak metamorphism reaching 635±20 °C. A suitable mineral standard requires that its petrological equilibrium first be demonstrated. Relicts and retrograde reaction textures are a guarantee of isotopic disequilibrium and heterogeneous ages within single crystal at the micrometre scale.

Crystal chemistry and hydrogen bonding of rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2 with variable OH, Cl, F

Three samples of the skarn mineral rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2, space group C2/c, a ≈7.6, b ≈ 18.5, c ≈ 15.5 Å, β ≈ 104°, with variable OH, Cl, F content were investigated by electron microprobe, single-crystal X-ray structure refinements, and Raman spectroscopy. “Rust1LCl” is a low chlorine rustumite Ca10(Si2O7)2(SiO4)(OH1.88F0.12)(Cl1.28,OH0.72) from skarns associated with the Rize batholith near Ikizedere, Turkey. “Rust2F” is a F-bearing rustumite Ca10(Si2O7)2(SiO4)(OH1.13F0.87) (Cl1 96OH0.04) from xenoliths in ignimbrites of the Upper Chegem Caldera, Northern Caucasus, Russia. “Rust3LClF” represents a low-Cl, F-bearing rustumite Ca10(Si2O7)2(SiO4)0.87(H4O4)0.13(OH1.01F0.99) (Cl1.00 OH1.00) from altered merwinite skarns of the Birkhin massif, Baikal Lake area, Eastern Siberia, Russia. Rustumite from Birkhin massif is characterized by a significant hydrogarnet-like or fluorine substitution at the apices of the orthosilicate group, leading to specific atomic displacements. The crystal structures including hydrogen positions have been refined from single-crystal X-ray data to R1 = 0.0205 (Rust1_LCl), R1 = 0.0295 (Rust2_F), and R1 = 0.0243 (Rust3_LCl_F), respectively. Depletion in Cl and replacement by OH is associated with smaller unit-cell dimensions. The substitution of OH by F leads to shorter hydrogen bonds O-H⋯F instead of O-H⋯OH. Raman spectra for all samples have been measured and confirm slight strengthening of the hydrogen bonds with uptake of F.This study discusses the complex crystal chemistry of the skarn mineral rustumite and may provide a wider understanding of the chemical reactions related to contact metamorphism of limestones.

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.

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.

REE and Hf distribution among mineral phases in the CV-CK clan : a way to explain present-day Hf isotopic variations in chondrites

Chondrites are among the most primitive objects in the Solar System and constitute the main building blocks of telluric planets. Among the radiochronometers currently used for dating geological events, Sm–Nd and Lu–Hf are both composed of refractory, lithophile element. They are thought to behave similarly as the parent elements (Sm and Lu) are generally less incompatible than the daughter elements (Nd and Hf) during geological processes. As such, their respective average isotopic compositions for the solar system should be well defined by the average of chondrites, called Chondritic Uniform Reservoir (CHUR). However, while the Sm–Nd isotopic system shows an actual spread of less than 4% in the average chondritic record, the Lu–Hf system shows a larger variation range of 28% [Bouvier A., Vervoort J. D. and Patchett P. J. (2008) The Lu–Hf and Sm–Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth Planet. Sci. Lett.273, 48–57]. To better understand the contrast between Sm–Nd and Lu–Hf systems, the REE and Hf distribution among mineral phases during metamorphism of Karoonda (CK) and Vigarano-type (CV) carbonaceous chondrites has been examined. Mineral modes were determined from elemental mapping on a set of five CK chondrites (from types 3–6) and one CV3 chondrite. Trace-element patterns are obtained for the first time in all the chondrite-forming minerals of a given class (CK chondrites) as well as one CV3 sample. This study reveals that REE are distributed among both phosphates and silicates. Only 30–50% of Sm and Nd are stored in phosphates (at least in chondrites types 3–5); as such, they are not mobilized during early stages of metamorphism. The remaining fraction of Sm and Nd is distributed among the same mineral phases; these elements are therefore not decoupled during metamorphism. Of the whole-rock total of Lu, the fraction held in phosphate decreases significantly as the degree of metamorphism increases (30% for types 3 and 4, less than 5% in type 6). In contrast to Lu, Hf is mainly hosted by silicates with little contribution from phosphates throughout the CK metamorphic sequence. A significant part of Sm and Nd are stored in phosphates in types 3–5, and these elements behave similarly during CK chondrite metamorphism. That explains the robustness of the Sm/Nd ratios in chondrites through metamorphism, and the slight discrepancies observed in the present-day isotopic Nd values in chondrites. On the contrary, Lu and Hf are borne by several different minerals and consequently they are redistributed during metamorphism–induced recrystallization. The Lu/Hf ratios are therefore significantly disturbed during chondrites metamorphism, leading to the high discrepancies observed in present-day Hf isotopic values in chondrites.

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.

A review of air–ice chemical and physical interactions (AICI): liquids, quasi-liquids, and solids in snow

Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termed quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale. Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions. Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

Ediacaran 2,500-km-long synchronous deep continental subduction in the West Gondwana Orogen

The deeply eroded West Gondwana Orogen is a major continental collision zone that exposes numerous occurrences of deeply subducted rocks, such as eclogites. The position of these eclogites marks the suture zone between colliding cratons, and the age of metamorphism constrains the transition from subduction-dominated tectonics to continental collision and mountain building. Here we investigate the metamorphic conditions and age of high-pressure and ultrahigh-pressure eclogites from Mali, Togo and NE-Brazil and demonstrate that continental subduction occurred within 20 million years over at least a 2,500-km-long section of the orogen during the Ediacaran. We consider this to be the earliest evidence of large-scale deep-continental subduction and consequent appearance of Himalayan-scale mountains in the geological record. The rise and subsequent erosion of such mountains in the Late Ediacaran is perfectly timed to deliver sediments and nutrients that are thought to have been necessary for the subsequent evolution of sustainable life on Earth.

Multiple metamorphic stages within an eclogite-facies terrane (Sesia Zone, Western Alps) revealed by Th-U-Pb petrochronology

Convergent plate margins typically experience a transition from subduction to collision dynamics as massive continental blocks enter the subduction channel. Studies of high-pressure rocks indicate that tectonic fragments are rapidly exhumed from eclogite facies to midcrustal levels, but the details of such dynamics are controversial.To understand the dynamics of a subduction channel we report the results of a petrochronological study from the central Sesia Zone, a key element of the internalWestern Alps.This comprises two polymetamorphic basement complexes (Eclogitic Micaschist Complex and Gneiss Minuti Complex) and a thin, dismembered cover sequence (Scalaro Unit) associated with pre-Alpine metagabbros and metasediments (Bonze Unit). Structurally controlled samples from three of these units (Eclogitic Micaschist Complex and Scalaro-Bonze Units) yield unequivocal petrological and geochronological evidence of two distinct high-pressure stages. Ages (U-Th-Pb) of growth zones in accessory allanite and zircon, combined with inclusion and textural relationships, can be tied to the multi-stage evolution of single samples.Two independent tectono-metamorphic ‘slices’ showing a coherent metamorphic evolution during a given time interval have been recognized: the Fondo slice (which includes Scalaro and Bonze rocks) and the Druer slice (belonging to the Eclogitic Micaschist Complex).The new data indicate separate stages of deformation at eclogite-facies conditions for each recognized independent kilometer-sized tectono-metamorphic slice, between ~85 and 60 Ma, with evidence of intermittent decompression (∆P~0.5 GPa) within only the Fondo slice. The evolution path of the Druer slice indicates a different P-T-time evolution with prolonged eclogite-facies metamorphism between ~85 and 75Ma. Our approach, combining structural, petrological and geochronological techniques, yields field-based constraints on the duration and rates of dynamics within a subduction channel.

The mobility of Nb in rutile-saturated NaCl- and NaF-bearing aqueous fluids from 1-6.5 GPa and 300-800 °C

Rutile (TiO2) is an important host phase for high field strength elements (HFSE) such as Nb in metamorphic and subduction zone environments. The observed depletion of Nb in arc rocks is often explained by the hypothesis that rutile sequesters HFSE in the subducted slab and overlying sediment, and is chemically inert with respect to aqueous fluids evolved during prograde metamorphism in the forearc to subarc environment. However, field observations of exhumed terranes, and experimental studies, indicate that HFSE may be soluble in complex aqueous fluids at high pressure (i.e., >0.5 GPa) and moderate to high temperature (i.e., >300 degrees C). In this study, we investigated experimentally the mobility of Nb in NaCl- and NaF-bearing aqueous fluids in equilibrium with Nb-bearing rutile at pressure-temperature conditions applicable to fluid evolution in arc environments. Niobium concentrations in aqueous fluid at rutile saturation were measured directly by using a hydrothermal diamond-anvil cell (HDAC) and synchrotron X-ray fluorescence (SXRF) at 2.1 to 6.5 GPa and 300-500 degrees C, and indirectly by performing mass loss experiments in a piston-cylinder (PC) apparatus at similar to 1 GPa and 700-800 degrees C. The concentration of Nb in a 10 wt% NaCl aqueous fluid increases from 6 to 11 mu g/g as temperature increases from 300 to 500 degrees C, over a pressure range from 2.1 to 2.8 GPa, consistent with a positive temperature dependence. The concentration of Nb in a 20 wt% NaCl aqueous fluid varies from 55 to 150 mu g/g at 300 to 500 degrees C, over a pressure range from 1.8 to 6.4 GPa; however, there is no discernible temperature or pressure dependence. The Nb concentration in a 4 wt% NaF-bearing aqueous fluid increases from 180 to 910 mu g/g as temperature increases from 300 to 500 degrees C over the pressure range 2.1 to 6.5 GPa. The data for the F-bearing fluid indicate that the Nb content of the fluid exhibits a dependence on temperature between 300 and 500 degrees C at >= 2 GPa, but there is no observed dependence on pressure. Together, the data demonstrate that the hydrothermal mobility of Nb is strongly controlled by the composition of the fluid, consistent with published data for Ti. At all experimental conditions, however, the concentration of Nb in the fluid is always lower than coexisting rutile, consistent with a role for rutile in moderating the Nb budget of arc rocks.

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.

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.