Garnet–monazite rare earth element relationships in sub-solidus metapelites: a case study from Bhutan

A key aim of modern metamorphic geochronology is to constrain precise and accurate rates and timescales of tectonic processes. One promising approach in amphibolite and granulite-facies rocks links the geochronological information recorded in zoned accessory phases such as monazite to the pressure–temperature information recorded in zoned major rock-forming minerals such as garnet. Both phases incorporate rare earth elements (REE) as they crystallize and their equilibrium partitioning behaviour potentially provides a useful way of linking time to temperature. We report REE data from sub-solidus amphibolite-facies metapelites from Bhutan, where overlapping ages, inclusion relationships and Gd/Lu ratios suggest that garnet and monazite co-crystallized. The garnet–monazite REE relationships in these samples show a steeper pattern across the heavy (H)REE than previously reported. The difference between our dataset and the previously reported data may be due to a temperature-dependence on the partition coefficients, disequilibrium in either dataset, differences in monazite chemistry or the presence or absence of a third phase that competed for the available REE during growth. We urge caution against using empirically-derived partition coefficients from natural samples as evidence for, or against, equilibrium of REE-bearing phases until monazite–garnet partitioning behaviour is better constrained

Kyanite petrogenesis in migmatites: resolving melting and metamorphic signatures

Aluminosilicates (kyanite, sillimanite and andalusite) are useful pressure–temperature (P–T) indicators that can form in a range of rock types through different mineral reactions, including those that involve partial melting. However, the presence of xenocrystic or inherited grains may lead to spurious P–T interpretations. The morphologies, microtextural positions, cathodoluminescence responses and trace element compositions of migmatite-hosted kyanite from Eastern Bhutan were investigated to determine whether sub-solidus kyanite could be distinguished from kyanite that crystallised directly from partial melt, or from kyanite that grew peritectically during muscovite dehydration reactions. Morphology and cathodoluminescence response were found to be the most reliable petrogenetic indicators. Trace element abundances generally support petrographic evidence, but protolith bulk composition exerts a strong control over absolute element abundance in kyanite. Sample-normalised concentrations show distinctive differences between petrogenetic types, particularly for Mg, Ti, V, Cr, Mn, Fe and Ge. LA-ICP-MS element maps, particularly combined to show Cr/V, provide additional information about changing geochemical environments during kyanite growth. Most kyanite in the studied migmatitic leucosomes is of sub-solidus origin, with less widespread evidence for peritectic crystallisation. Where present, grain rims commonly crystallised directly from the melt; however, entire grains crystallised exclusively from melt are rare. The presence of kyanite in leucosomes does not, therefore, necessarily constrain the P–T conditions of melting, and the mechanism of growth should be determined before using kyanite as a P–T indicator. This finding has significant implications for the interpretation of kyanite-bearing migmatites as representing early stages of melting during Himalayan evolution

Using U-Th-Pb petrochronology to determine rates of ductile thrusting: time windows into the Main Central Thrust, Sikkim Himalaya

Quantitative constraints on the rates of tectonic processes underpin our understanding of the mechanisms that form mountains. In the Sikkim Himalaya, late structural doming has revealed time-transgressive evidence of metamorphism and thrusting that permit calculation of the minimum rate of movement on a major ductile fault zone, the Main Central Thrust (MCT), by a novel methodology. U-Th-Pb monazite ages, compositions, and metamorphic pressure-temperature determinations from rocks directly beneath the MCT reveal that samples from ~50 km along the transport direction of the thrust experienced similar prograde, peak, and retrograde metamorphic conditions at different times. In the southern, frontal edge of the thrust zone, the rocks were buried to conditions of ~550°C and 0.8 GPa between ~21 and 18 Ma along the prograde path. Peak metamorphic conditions of ~650°C and 0.8–1.0 GPa were subsequently reached as this footwall material was underplated to the hanging wall at ~17–14 Ma. This same process occurred at analogous metamorphic conditions between ~18–16 Ma and 14.5–13 Ma in the midsection of the thrust zone and between ~13 Ma and 12 Ma in the northern, rear edge of the thrust zone. Northward younging muscovite 40Ar/39Ar ages are consistently ~4 Ma younger than the youngest monazite ages for equivalent samples. By combining the geochronological data with the >50 km minimum distance separating samples along the transport axis, a minimum average thrusting rate of 10 ± 3 mm yr−1 can be calculated. This provides a minimum constraint on the amount of Miocene India-Asia convergence that was accommodated along the MCT

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

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

Allanite U–Pb dating places new constraints on the high‐pressure to high‐temperature evolution of the deep Himalayan crust

During continental collision, crustal rocks are buried, deformed, transformed and exhumed. The rates, timescales and tectonic implications of these processes are constrained through the sequence and conditions of metamorphic reactions in major and accessory phases. Petrographic, isotopic and elemental data from metabasite samples in NW Bhutan, eastern Himalaya, suggest initial equilibration under high‐pressure (plagioclase‐absent and rutile‐present) conditions, followed by decompression to lower pressure conditions at high‐temperatures that stabilized plagioclase, orthopyroxene and ilmenite. Field observations and chemical indicators suggest equilibration under the lower pressure conditions is likely linked to the infiltration of melt from the host metasedimentary rocks. The metabasites preserve two metamorphic growth stages of chemically‐and petrographically distinct allanite that temporally overlap two stages of zircon growth. Allanite cores and zircon mantles grew at c. 19 ± 2 and 17–15.5 Ma respectively, linked texturally and chemically to the high‐pressure evolution. Symplectitic rims on embayed allanite cores, wholly symplectized Aln–Ilm and Aln–Cpx grains, and high U zircon rims grew at c. 15.5–14.5 Ma, linked chemically to the presence of melt and lower pressure, high‐temperature conditions. A single garnet Lu–Hf date is interpreted as geologically meaningless, with the bulk rock composition modified by melt infiltration after garnet formation. The open system evolution of these rocks precludes precise determination of the reactive bulk composition during metamorphic evolution and thus absolute conditions, especially during the early high‐pressure evolution. Despite these limitations, we show that combined geochemical and petrographic datasets are still able to provide insights into the rates and timescales of deep orogenic processes. The data suggest a younger and shallower evolution for the NW Bhutan metabasites compared to similar rocks in the central and eastern Himalayas

Figure S8: Garnet–monazite rare earth element relationships in sub-solidus metapelites: a case study from Bhutan

Supplementary field photographs, thin section information, monazite chemical maps and laser-hole photographs for LG-09-5