Argon behaviour in an inverted Barrovian sequence, Sikkim Himalaya: the consequences of temperature and timescale on ⁴⁰Ar/³⁹Ar mica geochronology

40Ar/39Ar dating of metamorphic rocks sometimes yields complicated datasets which are difficult to interpret in terms of timescales of the metamorphic cycle. Single-grain fusion and step-heating data were obtained for rocks sampled through a major thrust-sense shear zone (the Main Central Thrust) and the associated inverted metamorphic zone in the Sikkim region of the eastern Himalaya. This transect provides a natural laboratory to explore factors influencing apparent 40Ar/39Ar ages in similar lithologies at a variety of metamorphic pressure and temperature (P–T) conditions. The 40Ar/39Ar dataset records progressively younger apparent age populations and a decrease in within-sample dispersion with increasing temperature through the sequence. The white mica populations span ~ 2–9 Ma within each sample in the structurally lower levels (garnet grade) but only ~ 0–3 Ma at structurally higher levels (kyanite-sillimanite grade). Mean white mica single-grain fusion population ages vary from 16.2 ± 3.9 Ma (2σ) to 13.2 ± 1.3 Ma (2σ) from lowest to highest levels. White mica step-heating data from the same samples yields plateau ages from 14.27 ± 0.13 Ma to 12.96 ± 0.05 Ma. Biotite yield older apparent age populations with mean single-grain fusion dates varying from 74.7 ± 11.8 Ma (2σ) at the lowest structural levels to 18.6 ± 4.7 Ma (2σ) at the highest structural levels; the step-heating plateaux are commonly disturbed. Temperatures > 600 °C at pressures of 0.4–0.8 GPa sustained over > 5 Ma, appear to be required for white mica and biotite ages to be consistent with diffusive, open-system cooling. At lower temperatures, and/or over shorter metamorphic timescales, more 40Ar is retained than results from simple diffusion models suggest. Diffusion modelling of Ar in white mica from the highest structural levels suggests that the high-temperature rocks cooled at a rate of ~ 50–80 °C Ma− 1, consistent with rapid thrusting, extrusion and exhumation along the Main Central Thrust during the mid-Miocene

A novel palaeoaltimetry proxy based on spore and pollen wall chemistry

Understanding the uplift history and the evolution of high altitude plateaux is of major interest to a wide range of geoscientists and has implications for many disparate fields. Currently the majority of palaeoaltimetry proxies are based on detecting a physical change in climate in response to uplift, making the relationship between uplift and climate difficult to decipher. Furthermore, current palaeoaltimetry proxies have a low degree of precision with errors typically greater than 1 km. This makes the calculation of uplift histories and the identification of the mechanisms responsible for uplift difficult to determine. Here we report on advances in both instrumentation and our understanding of the biogeochemical structure of sporopollenin that are leading to the establishment of a new proxy to track changes in the flux of UV-B radiation over geological time. The UV-B proxy is based on quantifying changes in the concentration of UV-B absorbing compounds (UACs) found in the spores and pollen grains of land plants, with the relative abundances of UACs increasing on exposure to elevated UV-B radiation. Given the physical relationship between altitude and UV-B radiation we suggest that the analysis of sporopollenin chemistry, specifically changes in the concentration of UACs, may offer the basis for the first climate independent palaeoaltimetry proxy. Owing to the ubiquity of spores and pollen in the fossil record our proposed proxy has the potential to enable the reconstruction of the uplift history of high altitude plateaux at unprecedented levels of fidelity, both spatially and temporally

Developing an inverted Barrovian sequence; insights from monazite petrochronology

In the Himalayan region of Sikkim, the well-developed inverted metamorphic sequence of the Main Central Thrust (MCT) zone is folded, thus exposing several transects through the structure that reached similar metamorphic grades at different times. In-situ LA-ICP-MS U–Th–Pb monazite ages, linked to pressure–temperature conditions via trace-element reaction fingerprints, allow key aspects of the evolution of the thrust zone to be understood for the first time. The ages show that peak metamorphic conditions were reached earliest in the structurally highest part of the inverted metamorphic sequence, in the Greater Himalayan Sequence (GHS) in the hanging wall of the MCT. Monazite in this unit grew over a prolonged period between ~37 and 16 Ma in the southerly leading-edge of the thrust zone and between ~37 and 14.5 Ma in the northern rear-edge of the thrust zone, at peak metamorphic conditions of ~790 ◦C and 10 kbar. Monazite ages in Lesser Himalayan Sequence (LHS) footwall rocks show that identical metamorphic conditions were reached ~4–6 Ma apart along the ~60 km separating samples along the MCT transport direction. Upper LHS footwall rocks reached peak metamorphic conditions of ~655 ◦C and 9 kbar between ~21 and 16 Ma in the more southerly-exposed transect and ~14.5–12 Ma in the northern transect. Similarly, lower LHS footwall rocks reached peak metamorphic conditions of ~580 ◦C and 8.5 kbar at ~16 Ma in the south, and 9–10 Ma in the north. In the southern transect, the timing of partial melting in the GHS hanging wall (~23–19.5 Ma) overlaps with the timing of prograde metamorphism (~21 Ma) in the LHS footwall, confirming that the hanging wall may have provided the heat necessary for the metamorphism of the footwall. Overall, the data provide robust evidence for progressively downwards-penetrating deformation and accretion of original LHS footwall material to the GHS hanging wall over a period of ~5 Ma. These processes appear to have occurred several times during the prolonged ductile evolution of the thrust. The preserved inverted metamorphic sequence therefore documents the formation of sequential ‘paleothrusts’ through time, cutting down from the original locus of MCT movement at the LHS–GHS protolith boundary and forming at successively lower pressure and temperature conditions. The petrochronologic methods applied here constrain a complex temporal and thermal deformation history, and demonstrate that inverted metamorphic sequences can preserve a rich record of the duration of progressive ductile thrusting

Li and ?7Li in Himalayan rivers: Proxies for silicate weathering?

This paper presents the first systematic survey of lithium and its isotopes in the dissolved load and suspended and bed sediments of tributaries of the Ganges, both before and after the monsoon. Locations were chosen in order to cover catchments draining both silicates and carbonates, at high (2000–4000 m) and low (550–1300 m) altitudes. Modelling of the dissolved composition shows that the Li / Ca ratio of the silicate endmember in Himalayan rivers is at least an order-of-magnitude higher than that of the carbonate endmember. Most of the dissolved Li (> 90%) is derived from silicates even in carbonate-dominated catchments. While the Sr-isotope composition of the dissolved load reflects that of the bedrock, the main control on its Li-isotope composition is fractionation during weathering. Fractionation between the dissolved and suspended load in silicate-dominated catchments is greatest at high altitude and lower at low altitude where weathering is more intense. Tributaries draining silicates have lower dissolved ?7Li values (by 2.3‰ to 4.2‰) following the monsoon when weathering is more intense because of higher runoff and elevated temperatures. Our data suggest that riverine Li fluxes largely reflect silicate weathering rates, while riverine ?7Li varies with weathering intensity. As rivers presently contribute ?50% of the Li input to the oceans, seawater Li concentrations and ?7Li show potential as proxies for global silicate weathering processes

Relative contributions of silicate and carbonate rocks to riverine Sr fluxes in the headwaters of the Ganges

Exhumation of the Himalayan-Tibetan orogen is implicated in the marked rise in seawater 87Sr/86Sr ratios since 40 Ma. However both silicate and carbonate rocks in the Himalaya have elevated 87Sr/86Sr ratios and there is disagreement as to how much of the 87Sr flux is derived from silicate weathering. Most previous studies have used element ratios from bedrock to constrain the proportions of silicate- and carbonate-derived Sr in river waters. Here we use arrays of water compositions sampled from the head waters of the Ganges in the Indian and Nepalese Himalaya to constrain the end-member element ratios. The compositions of tributaries draining catchments restricted to a limited range of geological units can be described by two-component mixing of silicate and carbonate-derived components and lie on a plane in multicomponent composition space. Key elemental ratios of the carbonate and silicate components are determined by the intersection of the tributary mixing plane with the planes Na = 0 for carbonate and constant Ca/Na for silicate. The fractions of Sr derived from silicate and carbonate sources are then calculated by mass-balance in Sr-Ca-Mg-Na composition space. Comparison of end-member compositions with bedrock implies that secondary calcite deposition may be important in some catchments and that dissolution of low-Mg trace calcite in silicate rocks may explain discrepancies in Sr-Ca-Na-Mg covariation. Alternatively, composition-dependent precipitation or incongruent dissolution reactions may rotate mixing trends on cation-ratio diagrams. However the calculations are not sensitive to transformations of the compositions by incongruent dissolution or precipitation processes provided that the transformed silicate and carbonate component vectors are constrained. Silicates are calculated to provide 50% of the dissolved Sr flux from the head waters of the Ganges assuming that discrepancies between Ca-Mg-Na covariation and the silicate rock compositions arise from addition of trace calcite. If the Ca-Mg-Na mixing plane is rotated by composition-dependent secondary calcite deposition, this estimate would be increased. Moreover, when 87Sr/86Sr ratios of the Sr inputs are considered, silicate Sr is responsible for 70 ± 16% (1σ) of the 87Sr flux forcing changes in seawater Sr-isotopic composition. Since earlier studies predict that silicate weathering generates as little as 20% of the total Sr flux in Himalayan river systems, this study demonstrates that the significance of silicate weathering can be greatly underestimated if the processes that decouple the water cation ratios from those of the source rocks are not properly evaluated

Erosion history of the Tibetan Plateau since the last Interglacial: constraints from the first studies of cosmogenic 10Be from Tibetan bedrock

The cosmogenic 10Be exposure histories of in situ bedrock surfaces from the Tibetan Plateau indicate low erosion rates of <30 mm/ka in southern and central Tibet during the last interglacial–glacial cycle that contrast strongly with unusually rapid erosion rates (60–2000 mm/ka) for Kunlun in northern Tibet during the Holocene, comparable with published values from the Himalaya. By comparing apatite fission-track ages with cosmogenic data, erosion rates in southern Tibet appear to be decelerating since the Miocene, whereas in the Kunlun, erosion rates have accelerated over the same timescale. Such secular changes suggest that the southern and central regions of the plateau had formed their present flat relief by the Pleistocene. Unusually high erosion rates along the northern margin of the plateau may reflect intense tectonic activity during the Holocene. These findings indicate that over much of the high plateau erosion rates are exceptionally low, and therefore the sources of detritus carried by the great Asian rivers that rise in Tibet lie overwhelmingly in bedrocks at lower altitudes. This study illustrates the potential of cosmogenic studies for unraveling the most recent phase of the erosion/exhumation history of orogenic belts that cannot be resolved by either Ar-isotope or fission-track thermochronometers. Autho

The geology and tectonics of central Bhutan

Lithotectonic mapping, metamorphic observations and U–Pb zircon ages underpin a substantial revision of central Bhutan geology, notably a more extensive and continuous outcrop of the Tethyan Sedimentary Series (TSS) than previously mapped. Metamorphic grade in the TSS increases downward towards a basal north-vergent tectonic contact with the underlying Greater Himalayan Series (GHS), interpreted as a southward continuation of the South Tibetan Detachment (STD). Miocene (c. 17–20 Ma) leucogranite sheets are associated with the STD in this region but appear to diminish southwards. Two leucogranite dykes that cross-cut TSS structures yield ages of 17.8 ± 0.2 and 17.9 ± 0.5 Ma. A 500 ± 4 Ma (U–Pb zircon) metamorphosed ash bed in the Pele La Group within the psammite-dominated lower TSS yields the first direct isotopic age for the TSS in the eastern Himalaya, confirming existing age constraints from detrital zircon and fossil studies. A continuation of the Paro metasedimentary unit underlying the GHS was mapped near Wangdue Phodrang. Our observations, notably the exposure of a wholly ductile STD so far south and the significance of large nappe-like structures in the TSS, prompt a major revision to the geological map of the Bhutan Himalaya and require a reassessment of tectonic interpretations of the Bhutan Himalaya

The identification and significance of pure sediment-derived granites

© 2017 The Author(s)The characterization of the geochemical reservoirs of the Earth's continental crust, including the determination of representative upper and lower crustal compositions, underpins our understanding of crustal evolution. The classic I- and S-type granite classification has often been invoked to distinguish between melts derived from igneous protoliths and those derived from the melting of a sedimentary source. Recent geochemical studies suggest that most granites, even those cited as typical examples of ‘S-type’, show evidence for a mixture of mantle and upper crustal sources, thereby implying that granite formation is evidence for overall crustal growth. We have examined the source of leucogranite bodies in one of the world's youngest collisional orogens using novel zircon techniques that can resolve the presence of even minor mantle contributions. 232 zircons from 12 granites from the Bhutan Himalaya were analysed by in-situ techniques for O, Hf and U–Pb isotopic signatures. In combination with data from the granite host rocks, our data show that the Himalayan leucogranites were derived solely from metamorphosed crustal sediments, and do not record any mantle contribution. This finding is consistent with the time-lag between crustal thickening and widespread crustal melting, and the heat-producing capacities of the pelitic source rocks. We conclude that Himalayan leucogranites provide a more suitable type locality for ‘S-type’ granites than the Lachlan area in South-East Australia where the term was first defined. The Himalayan leucogranites therefore provide evidence that syn-orogenic melting during collisional events does not necessarily result in crustal growth. Importantly, crustal growth models should not always assume that crustal growth is achieved during collisional orogenesis