Thermobarometric constraints on the thermal history of the Main Central Thrust Zone and Tibetan Slab, eastern Nepal Himalaya

ABSTRACT The Main Central Thrust (MCT) south of Mt Everest in eastern Nepal is a 3 to 5km thick shear zone separating chlorite-bearing schist in the lower plate from sillimanite-bearing migmatitic gneiss in the overlying Tibetan Slab. The metamorphic grade increases through the MCT zone toward structurally higher levels. Previous workers have suggested that either post- or synmetamorphic thrust movement has caused this inversion of metamorphic isograds. In an effort to quantify the increase in grade and to constrain proposed structural relations between metamorphism and slip on the fault, four well-calibrated thermobarometers were applied to pelitic samples collected along two cross-strike transects through the MCT zone and Tibetan Slab. Results show an increase in apparent temperature up-section in the MCT zone from 778 K to 990 K and a decrease in temperature to ∼850 K in the lower Tibetan Slab, which is consistent with synmetamorphic thrust movement. A trend in calculated pressures across this section is less well-defined but, on average, decreases up-section with a gradient of ∼28MPa/km, resembling a lithostatic gradient. Pressure-temperature paths for zoned garnets from samples within the MCT zone, modelled using the Gibbs' Method, show a significant decrease in temperature and a slight decrease in pressure from core to rim, which might be expected for upper plate rocks during synmetamorphic thrust movement. Samples from the uppermost Tibetan Slab yield higher temperatures and pressures than those from the lower Tibetan Slab, which may be evidence for later‘resetting’ of thermobarometers by intrusion of the large amounts of leucogranite at that structural level.     

Nd isotopic data reveal the material and tectonic nature of the Main Central Thrust zone in Nepal Himalaya

The Main Central Thrust (MCT) is a tectono-metamorphic boundary between the Higher Himalayan crystallines (HHC) and Lesser Himalayan metasediments (LHS), reactivated in the Tertiary, but which had already formed as a collisional boundary in the Early Paleozoic. To investigate the nature of the MCT, we analyzed whole-rock Nd isotopic ratios of rocks from the MCT and surrounding zones in the Taplejung–Ilam area of far-eastern Nepal, Annapurna–Galyang area of central Nepal, and Maikot–Barekot area of western Nepal. We define the MCT zone as a ductile–brittle shear zone between the upper MCT (UMCT) and lower MCT (LMCT). The protoliths of the MCT zone may provide critical constraints on the tectonic evolution of the Himalaya. The LHS is lithostratigraphically divided into the upper and lower units. In the Taplejung–Ilam area, different lithologic units and their ε_Nd (0) values are as follows; HHC (− 10.0 to − 18.1), MCT zone (− 18.5 to − 26.2), upper LHS unit (− 17.2), and lower LHS unit (− 22.0 to − 26.9). There is a distinct gap in the ε_Nd (0) values across the UMCT except for the southern frontal edge of the Ilam nappe. In the Annapurna–Galyang and Maikot–Barekot areas, different lithologic units and their ε_Nd (0) values are as follows; HHC (− 13.9 to − 17.7), MCT zone (− 23.8 to − 26.2 except for an outlier of − 12.4), upper LHS unit (− 15.6 to − 26.8), and lower LHS unit (− 24.9 to − 26.8). These isotopic data clearly distinguish the lower LHS unit from the HHC. Combining these data with the previously published data, the lowest ε_Nd (0) value in the HHC is − 19.9. We regard rocks with ε_Nd (0) values below − 20.0 as the LHS. In contrast, rocks with those above − 19.9 are not always the HHC, and some parts of them may belong to the LHS due to the overlapping Nd isotopic ratio between the HHC and LHS. Most rocks of the MCT zone have Nd isotopic ratios similar to those of the LHS, but very different from those of the HHC. The spatial patterns in the distribution of ε_Nd (0) value around the UMCT suggest no substantial structural mixing of the HHC and LHS during the UMCT activities in the Tertiary. A discontinuity in the spatial distribution of ε_Nd (0) values is laterally continuous along the UMCT throughout the Himalayas. These facts support the theory that the UMCT was originally a material boundary between the HHC and LHS, suggesting the MCT zone was mainly developed with undertaking a role of sliding planes during overthrusting of the HHC in the Tertiary.     

Structural kinematics,metamorphic P–T profiles and zircon geochronology across the Greater Himalayan Crystalline Complex in south‐central Tibet: implication for a revised channel flow

A specific question about the Himalayas is whether the orogeny grew by distributed extrusion or discrete thrusting. To place firm constraints on tectonic models for the orogeny, kinematic, thermobarometric and geochronological investigations have been undertaken across the Greater Himalayan Crystalline Complex (GHC) in the Nyalam region, south‐central Tibet. The GHC in this section is divided into the lower, upper and uppermost GHC, with kinematically top‐to‐the‐south, alternating with top‐to‐the‐north shear senses. A new thrust named the Nyalam thrust is recognized between the lower and upper GHC, with a 3 kbar pressure reversion, top‐to‐the‐south thrust sense, and was active after the exhumation of the GHC. Peak temperature reached ～749 °C in the cordierite zone, and decreased southwards to 633–667 °C in the kyanite and sillimanite‐muscovite zones, and northwards to greenschist facies at the top of the South Tibetan Detachment System (STDS). Pressure at peak temperature reached a maximum value in the kyanite zone of 9.0–12.6 kbar and decreased northwards to ～4.1 kbar in the cordierite zone. Zircon U‐Pb ages of a sillimanite migmatite and an undeformed leucogranite dyke cutting the mylonitized rocks in the STDS reveal a long‐lived partial melting of the GHC, which initiated at 39.7–34 Ma and ceased at 14.1 Ma. Synthesizing the obtained and collected results, a revised channel flow model is proposed by considering the effect of heat advection and convection by melt and magma migration. Our new model suggests that distributed processes like channel flow dominated during the growth of the Himalayan orogen, while discrete thrusting occurred in a later period as a secondary process.     

Petrological constraints on the tectonic setting of the Kathmandu Nappe in the Langtang–Gosainkund–Helambu regions,Central Nepal Himalaya

The Gosainkund–Helambu region in central Nepal occupies a key area for the development of Himalayan kinematic models, connecting the well‐investigated Langtang area to the north with the Kathmandu Nappe (KN), whose interpretation is still debated, to the south. In order to understand the structural and metamorphic architecture of the Greater Himalayan Sequence (GHS) in this region, a detailed petrological study was performed, focusing on selected metapelite samples from both the Gosainkund–Helambu and Langtang transects. The structurally lowest sample investigated belongs to the Lesser Himalayan Sequence; its metamorphic evolution is characterized by a narrow hairpin P–T path with peak P–T conditions of 595 ± 25 °C, 7.5 ± 1 kbar. All of the other samples here investigated belong to the GHS. Along the Langtang section, two tectono‐metamorphic units have been distinguished within the GHS: the Lower Greater Himalayan Sequence (L‐GHS), characterized by peak P–T conditions at 728 ± 11 °C, 10 ± 0.5 kbar (corresponding to a T/depth ratio of 22 ± 1 °C km^?1), and the structurally higher Upper Greater Himalayan Sequence, with peak metamorphic conditions at 780 ± 20 °C, 7.8 ± 0.8 kbar (corresponding to a T/depth ratio of 31 ± 4 °C km^?1). This confirms the existence of a main tectono‐metamorphic discontinuity within the GHS, as previously suggested by other authors. The results of petrological modelling of the metapelites from the Gosainkund–Helambu section show that this region is entirely comprised within a sub‐horizontal and thin L‐GHS unit: the estimated peak metamorphic conditions of 734 ± 19 °C, 10 ± 0.8 kbar correspond to a uniform T/depth ratio of 23 ± 3 °C km^?1. The metamorphic discontinuity identified along the Langtang transect and dividing the GHS in two tectono‐metamorphic units is located at a structural level too high to be intersected along the Gosainkund–Helambu section. Our results have significant implications for the interpretation of the KN and provide a contribution to the more general discussion of the Himalayan kinematic models. We demonstrate that the structurally lower unit of the KN (known as Sheopuri Gneiss) can be correlated with the L‐GHS unit; this result strongly supports those models that correlate the KN to the Tethyan Sedimentary Sequence and that suggest the merging of the South Tibetan Detachment System and the Main Central Thrust on the northern side of the KN. Moreover we speculate that, in this sector of the Himalayan chain, the most appropriate kinematic model able to explain the observed tectono‐metamorphic architecture of the GHS is the duplexing model, or hybrid models which combine the duplexing model with another end‐member model.     

Petrology and tectonic significance of metabasite slivers in the Lesser and Higher Himalayan domains of Sikkim,India

The metamorphic history of the Himalayas has been constrained mostly through studies of the ubiquitous metapelitic rocks. Non‐eclogitic metabasite rock lenses that occur intercalated with the metapelites have received little attention and it is not clear whether they share a common metamorphic history. This study reports the results of a petrological study of the metabasite lenses (dm³–m³) from the Lesser Himalayan (LH) and the Higher Himalayan (HH) domains in Sikkim. These have similar bulk chemical compositions and chemical affinities (sub‐alkaline tholeiitic basalts), with plagioclase and amphibole as the dominant mineralogical constituents. Garnet and clinopyroxene occur in some samples depending on small variations in bulk chemistry; and orthopyroxene is developed as a retrograde phase in some rocks. Minor phases are ilmenite, chlorite, titanite and rutile. The rocks were metamorphosed at similar conditions (～9–12 kbar, 800 °C). Minor differences in bulk chemical composition lead to different phase assemblages and mineral chemistry in adjacent metabasite lenses, a feature that is used to demonstrate that metamorphic conditions (peak P–T as well as retrograde P–T path) can be reliably retrieved through a combination of pseudosection analysis and kinetically constrained individual thermobarometry. The peak P–T conditions of the metabasites from this region are independent of the present geographic or tectonic (i.e. within the LH or HH) location of the samples and they differ from the conditions at which the regional metapelites (i.e. metapelites not immediately adjacent to the metabasite lenses) were metamorphosed. Metapelites that are immediately adjacent to the metabasite lenses differ in their appearance, phase assemblage and recorded P–T history from those of the regional metapelites, either because they were emplaced as slivers along with the metabasites, or because they were modified when they came in contact with the metabasites. The retrograde P–T paths of the LH and HH metabasites are different: the HH samples underwent steep decompression whereas the LH followed a more gentle exhumation path. The P–T conditions of peak metamorphism (9–12 kbar, 800 °C) are commensurate with a thermal perturbation at the base of a crust of average thickness and may be the signature of a widespread (samples found across different regions in the Himalaya) and long‐lasting (e.g. homogeneous garnet compositions) crustal underplating event that occurred during the early stages (?subduction) of the Himalayan orogeny, or earlier if the metamorphism was pre‐Himalayan.     

P–T evolution of the Precambrian Metamorphic Complex,NW Iran: a study of metapelitic rocks

The Mahneshan Metamorphic Complex (MMC) is one of the Precambrian terrains exposed in the northwest of Iran. The MMC underwent two main phases of deformation (D₁ and D₂) and at least two metamorphic events (M₁ and M₂). Critical metamorphic mineral assemblages in the metapelitic rocks testify to regional metamorphism under amphibolite‐facies conditions. The dominant metamorphic mineral assemblage in metapelitic rocks (M₁) is muscovite, biotite I, Garnet I, staurolite, Andalusite I and sillimanite. Peak metamorphism took place at 600–620°C and ∼7 kbar, corresponding to a depth of ca. 24 km. This was followed by decompression during exhumation of the crustal rocks up to the surface. The decrease of temperature and pressure during exhumation produced retrograde metamorphic assemblages (M₂). Secondary phases such as garnet II biotite II, Andalusite II constrain the temperature and pressure of M₂ retrograde metamorphism to 520–560°C and 2.5–3.5 kbar, respectively. The geothermal gradient obtained for the peak of metamorphism is 33°C km⁻¹, which indicates that peak metamorphism was of Barrovian type and occurred under medium‐pressure conditions. The MMC followed a ‘clockwise’ P–T path during metamorphism, consistent with thermal relaxation following tectonic thickening. The bulk chemistry of the MMC metapelites shows that their protoliths were deposited at an active continental margin. Together with the presence of palaeo‐suture zones and ophiolitic rocks around the high‐grade metamorphic rocks of the MMC, these features suggest that the Iranian Precambrian basement formed by an island‐arc type cratonization. Copyright © 2010 John Wiley & Sons, Ltd.     

P–T evolution of Precambrian eclogite in the Sveconorwegian orogen,SW Sweden

Conditions of the prograde, peak‐pressure and part of the decompressional P–T path of two Precambrian eclogites in the eastern Sveconorwegian orogen have been determined using the pseudosection approach. Cores of garnet from a Fe–Ti‐rich eclogite record a first prograde and syn‐deformational stage along a Barrovian gradient from ~670 °C and 7 kbar to 710 °C and 8.5 kbar. Garnet rims grew during further burial to 16.5–19 kbar at ~850–900 °C, along a steep dP/dT gradient. The pseudosection model of a kyanite‐bearing eclogite sample of more magnesian bulk composition confirms the peak conditions. Matrix reequilibration associated with subsequent near‐isothermal decompression and partial exhumation produced plagioclase‐bearing symplectites replacing kyanite and clinopyroxene at an estimated 850–870 °C and 10–11 kbar. The validity of the pseudosections is discussed in detail. It is shown that in pseudosection modelling the fractionation of FeO in accessory sulphides may cause a significant shift of field boundaries (here displaced by up to 1.5 kbar and 70 °C) and must not be neglected. Fast burial, exhumation and subsequent cooling are supported by the steepness of both the prograde and the decompressional P–T paths as well as the preservation of garnet growth zoning and the symplectitic reaction textures. These features are compatible with deep tectonic burial of the eclogite‐bearing continental crust as part of the underthrusting plate (Eastern Segment, continent Baltica) in a collisional setting that led to an effectively doubled crustal thickness and subsequent exhumation of the eclogites through tectonic extrusion. Our results are in accordance with regional structural and petrological relationships, which demonstrate foreland‐vergent partial exhumation of the eclogite‐bearing nappe along a basal thrust zone and support a major collisional stage at c. 1 Ga. We argue that the similarities between Sveconorwegian and Himalayan eclogite occurrences emphasize the modern style of Grenvillian‐aged tectonics.     

Inverted metamorphism and the Main Central Thrust: field relations and thermobarometric constraints from the Kishtwar Window, NW Indian Himalaya

Following the early Eocene collision of the Indian and Asian plates, intracontinental subduction occurred along the Main Central Thrust (MCT) zone in the High Himalaya. In the Kishtwar–Zanskar Himalaya, the MCT is a 2 km thick shear zone of high strain, distributed ductile deformation which emplaces the amphibolite facies High Himalayan Crystalline (HHC) unit south‐westwards over the lower greenschist facies Lesser Himalaya. An inverted metamorphic field gradient, mapped from the first appearance of garnet, staurolite and kyanite index minerals, is coincident with the high strain zone. Petrography and garnet zoning profiles indicate that rocks in the lower MCT zone preserve a prograde assemblage, whereas rocks in the HHC unit show retrograde equilibration. Thermobarometric results derived using THERMOCALC indicate a P–T increase of c. 180 °C and c. 400 MPa across the base of the MCT zone, which is a consequence of the syn‐ to postmetamorphic juxtaposition of M1 kyanite grade rocks of the HHC unit on a cooling path over biotite grade footwall rocks, which subsequently attain their peak (M2) during thrusting. Inclusion thermobarometry from the lower MCT zone reveals that M2 was accompanied by loading, and peak conditions of 537±38 °C and 860±120 MPa were attained. M1 kyanite assemblages in the HHC unit, which have not been overprinted by M2 fibrolitic sillimanite, were not significantly affected by M2, and conditions of equilibration are estimated as 742±53 °C and 960±180 MPa. There is no evidence for dissipative or downward conductive heating in the MCT zone. Instead, the primary control on the distribution of peak assemblages, represented by the index minerals, is postmetamorphic ductile thrusting in a downward propagating shear zone. Polymetamorphism and diachroneity of equilibration are also important controls on the thermal profile through the MCT zone and HHC unit.     

Constraining the conditions of Barrovian metamorphism in Sikkim,India: P–T–t paths of garnet crystallization in the Lesser Himalayan Belt

Garnet crystallization in metapelites from the Barrovian garnet and staurolite zones of the Lesser Himalayan Belt in Sikkim is modelled utilizing Gibbs free energy minimization, multi‐component diffusion theory and a simple nucleation and growth algorithm. The predicted mineral assemblages and garnet‐growth zoning match observations remarkably well for relatively tight, clockwise metamorphic P–T paths that are characterized by prograde gradients of ～30 °C kbar^?1 for garnet‐zone rocks and ～20 °C kbar^?1 for rocks from the staurolite zone. Estimates for peak metamorphic temperature increase up‐structure toward the Main Central Thrust. According to our calculations, garnet stopped growing at peak pressures, and protracted heating after peak pressure was absent or insignificant. Almost identical P–T paths for the samples studied and the metamorphic continuity of the Lesser Himalayan Belt support thermo‐mechanical models that favour tectonic inversion of a coherent package of Barrovian metamorphic rocks. Time‐scales associated with the metamorphism were too short for chemical diffusion to substantially modify garnet‐growth zoning in rocks from the garnet and staurolite zones. In general, the pressure of initial garnet growth decreases, and the temperature required for initial garnet growth was reached earlier, for rocks buried closer toward the MCT. Deviations from this overall trend can be explained by variations in bulk‐rock chemistry.     

Age and P–T conditions of the Gridino‐type eclogite in the Belomorian Province,Russia

Although eclogites in the Belomorian Province have been regarded as Archean in age and among the oldest in the world, there are also multiple studies that have proposed a Paleoproterozoic age. Here, we present new data for the Gridino‐type eclogites, which occur as boudins and metamorphosed dykes within tonalite–trondhjemite–granodiorite gneisses. Zircon from these eclogites has core and rim structures. The cores display high Th/U ratios (0.18–0.45), negative Eu anomalies and strong enrichment in HREE, and have Neoarchean U–Pb ages of c. 2.70 Ga; they are interpreted to be magmatic in origin. Zircon cores have δ¹⁸O of 5.64–6.07‰ suggesting the possibility of crystallization from evolved mantle‐derived magmas. In contrast, the rims, which include the eclogite facies minerals omphacite and garnet, are characterized by low Th/U ratios (<0.035) and flat HREE patterns, and yield U–Pb ages of c. 1.90 Ga; they are interpreted to be metamorphic in origin. Zircon rims have elevated δ¹⁸O of 6.23–6.80‰, which was acquired during eclogite facies metamorphism. Based on petrography and phase equilibria modelling, we recognize a prograde epidote amphibolite facies mineral assemblage, the peak eclogite facies mineral assemblage and a retrograde high‐P amphibolite facies mineral assemblage. The peak metamorphic conditions of 695–755°C at >18 kbar for the Gridino‐type eclogites suggest an apparent thermal gradient of <39–42°C/kbar for the Lapland–Kola collisional orogeny.     

The significance of Mn‐rich ilmenite and the determination of P–T paths from zoned garnet in metasedimentary rocks from the western Cape Breton Highlands,Nova Scotia

The Jumping Brook Metamorphic Suite in the western Cape Breton Highlands of Nova Scotia is part of an inverted Barrovian sequence that formed during a Late Silurian–Early Devonian promontory–promontory collision in the Canadian Appalachians. In this study, systematic discrepancies between geochemical observations and thermodynamic model predictions led to the discovery of a systematic relationship linking the style of garnet core isopleth intersection (GCII) to the pyrophanite (MnTiO₃) component of co‐existing ilmenite. Samples that yielded tight GCIIs at or near the garnet‐in curve were found to contain ilmenite with negligible pyrophanite components, whereas samples yielding GCIIs far removed (up to 105°C) from the garnet‐in curve were found to contain ilmenite with significant pyrophanite and/or ecandrewsite (ZnTiO₃) components. Based on petrographic and geochemical observations, Mn(±Zn)‐rich ilmenite are interpreted to have sequestered Mn throughout prograde metamorphism due to sluggish intracrystalline diffusion. The amount of reactive Mn input into the thermodynamic models from whole‐rock analyses were, in some cases, overestimated, resulting in garnet‐in curve topologies that extend to erroneously low P–T conditions. Modifications to the whole‐rock chemistry that account for Mn sequestration into ilmenite, however, yielded robust model results. Our results show that, in addition to uncertainties in thermodynamic data sets and phenomenon related to reaction kinetics, Mn‐rich ilmenite may superimpose additional complexities related to the interpretation of predicted equilibria involving garnet. Numerical simulations of garnet crystallization were used to infer P–T paths of metamorphism for one sample from the garnet zone (Mn corrected) and two samples from the staurolite zone (Mn uncorrected) of the inverted sequence. Model results are remarkably similar among the three samples and indicate that garnet crystallization occurred along relatively steep (31–37°C/km) clockwise P–T paths. The peak conditions of garnet crystallization and metamorphism (560–590°C, 7.4–8.0 kbar) are interpreted to have been attained approximately simultaneously, such that the paths are characterized by tight prograde‐to‐retrograde transitions. The hairpin nature of the P–T paths is interpreted to represent the onset of thrust‐related exhumation and isograd inversion along ductile shear zones, consistent with available field and geochronological constraints.     

Constraints from monazite and xenotime growth modelling in the MnCKFMASH‐PYCe system on the P–T path of a metapelite from Shillong‐Meghalaya Plateau: implications for the Indian shield assembly

The Palaeo‐Mesoproterozoic metapelite granulites from northern Garo Hills, western Shillong‐Meghalaya Gneissic Complex (SMGC), northeast India, consist of resorbed garnet, cordierite and K‐feldspar porphyroblasts in a matrix comprising shape‐preferred aggregates of biotite±sillimanite+quartz that define the penetrative gneissic fabric. An earlier assemblage including biotite and sillimanite occurs as inclusions within the garnet and cordierite porphyroblasts. Staurolite within cordierite in samples without matrix sillimanite is interpreted to have formed by a reaction between the sillimanite inclusion and the host cordierite during retrogression. Accessory monazite occurs as inclusions within garnet as well as in the matrix, whereas accessory xenotime occurs only in the matrix. The monazite inclusions in garnet contain higher Ca, and lower Y and Th/U than the matrix monazite outside resorbed garnet rims. On the other hand, matrix monazite away from garnet contains low Ca and Y, and shows very high Th/U ratios. The low Th/U ratios (<10) of the Y‐poor garnet‐hosted monazite indicate subsolidus formation during an early stage of prograde metamorphism. A calculated P–T pseudosection in the MnCKFMASH‐PYCe system indicates that the garnet‐hosted monazite formed at <3 kbar/600 °C (Stage A). These P–T estimates extend backward the previously inferred prograde P–T path from peak anatectic conditions of 7–8 kbar/850 °C based on major mineral equilibria. Furthermore, the calculated P–T pseudosections indicate that cordierite–staurolite equilibrated at ~5.5 kbar/630 °C during retrograde metamorphism. Thus, the P–T path was counterclockwise. The Y‐rich matrix monazite outside garnet rims formed between ~3.2 kbar/650 °C and ~5 kbar/775 °C (Stage B) during prograde metamorphism. If the effect of bulk composition change due to open system behaviour during anatexis is considered, the P–T conditions may be lower for Stage A (<2 kbar/525 °C) and Stage B (~3 kbar/600 °C to ~3.5 kbar/660 °C). Prograde garnet growth occurred over the entire temperature range (550–850 °C), and Stage‐B monazite was perhaps initially entrapped in garnet. During post‐peak cooling, the Stage‐B monazite grains were released in the matrix by garnet dissolution. Furthermore, new matrix monazite (low Y and very high Th/U ≤80, ~8 kbar/850–800 °C, Stage C), some monazite outside garnet rims (high Y and intermediate Th/U ≤30, ~8 kbar/800–785 °C, Stage D), and matrix xenotime (<785 °C) formed through post‐peak crystallization of melt. Regardless of textural setting, all monazite populations show identical chemical ages (1630–1578 Ma, ±43 Ma). The lithological association (metapelite and mafic granulites), and metamorphic age and P–T path of the northern Garo Hills metapelites and those from the southern domain of the Central Indian Tectonic Zone (CITZ) are similar. The SMGC was initially aligned with the southern parts of CITZ and Chotanagpur Gneissic Complex of central/eastern India in an ENE direction, but was displaced ~350 km northward by sinistral movement along the north‐trending Eastern Indian Tectonic Zone in Neoproterozoic. The southern CITZ metapelites supposedly originated in a back‐arc associated with subducting oceanic lithosphere below the Southern Indian Block at c. 1.6 Ga during the initial stage of Indian shield assembly. It is inferred that the SMGC metapelites may also have originated contemporaneously with the southern CITZ metapelites in a similar back‐arc setting.     

An anticlockwise P–T–t path at high‐pressure,high‐temperature conditions for a migmatitic gneiss from the island of Fjørtoft,Western Gneiss Region,Norway, indicates two burial events during the Caledonian orogeny

The Blåhø Nappe on the island of Fjørtoft, which represents an isolated portion of the Seve Nappe Complex in the Western Gneiss Region, Norway, has been suggested to have experienced two deep burial cycles during the Caledonian orogeny. However, evidence on this multiple burial process by the derivation of a pressure–temperature–time (P–T–t) path has never been given in the literature. In this study, the ‘diamondiferous’ kyanite–garnet gneiss from the Blåhø Nappe on Fjørtoft was revisited to determine if such a process was correct. Two types of garnet, porphyroblastic garnet‐1 and fine‐grained garnet‐2, were recognized in the gneiss. The core of garnet‐1 is poor in Ca and documents P–T conditions of 1.2–1.3 GPa at c. 880°C based on pseudosection modelling. The inner rims of garnet‐1 and the core of garnet‐2 are both richer in Ca, recording P–T conditions of 1.35–1.45 GPa and 770–820°C. Application of conventional geothermobarometry on the outer rim of garnet‐1 and the rim of garnet‐2 yielded retrograde P–T conditions of 0.75–0.90 GPa and 610–685°C. These estimates define an anticlockwise P–T path at pressures below 1.5 GPa. Accessory monazite was dated with the electron microscope. Relicts of detrital monazite in the gneiss point to Sveconorwegian and possibly also Cryogenian provenance for the detritus of the sedimentary protolith. Metamorphic monazite in the gneiss records a wide age range from 460 to 380 Ma, with a peak c. 435 Ma and a shoulder at 395 Ma. These data suggest that the original (Ediacaran?) Baltica margin sediment (gneiss protolith) was transported to the base of an overlying plate during the early Caledonian (pre‐Scandian) orogeny. A long residence time of the metasedimentary rock at this base caused its heating to 880°C and homogenization of the early garnet chemistry. The late Caledonian (Scandian) collision between Baltica and Laurentia led to further burial, during which the studied gneiss was close to the former surface of the downgoing continental plate and, thus, cooled. The reconstructed P–T–t path confirms the multiple burial history of the Blåhø Nappe but contradicts previous ideas of deep burial of the Fjørtoft gneiss to more than 100 km.     

Petrology,geochemistry and P–T–t path of lawsonite‐bearing retrograded eclogites in the Changning–Menglian orogenic belt,southeast Tibetan Plateau

The Changning–Menglian orogenic belt (CMOB) in the southeastern Tibetan Plateau, is considered as the main suture zone marking the closure of the Palaeo‐Tethys Ocean between the Indochina and Sibumasu blocks. Here, we investigate the recently discovered retrograded eclogites from this suture zone in terms of their petrological, geochemical and geochronological features, with the aim of constraining the metamorphic evolution and protolith signature. Two types of metabasites are identified: retrograded eclogites and mafic schists. The igneous precursors of the retrograded eclogites exhibit rare earth element distribution patterns and trace element abundance similar to those of ocean island basalts, and are inferred to have been derived from a basaltic seamount in an intra‐oceanic tectonic setting. In contrast, the mafic schists show geochemical affinity to arc‐related volcanics with the enrichment of Rb, Th and U, and depletion of Nb, Ta, Zr, Hf and Ti, and their protoliths possibly formed at an active continental margin tectonic setting. Retrograded eclogites are characterized by peak metamorphic mineral assemblages of garnet, omphacite, white mica, lawsonite and rutile, and underwent five‐stage metamorphic evolution, including pre‐peak prograde stage (M₁) at 18–19 kbar and 400–420°C, peak lawsonite‐eclogite facies (M₂) at 24–26 kbar and 520–530°C, post‐peak epidote–eclogite facies decompression stage (M₃) at 13–18 kbar and 530–560°C, subsequent amphibolite facies retrogressive stage (M₄) at 8–10 kbar and 530–600°C, and late greenschist facies cooling stage (M₅) at 5–8 kbar and 480–490°C. Laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) U–Pb spot analyses of zircon show two distinct age groups. The magmatic zircon from both the retrograded eclogite and mafic schist yielded protolith ages of 451 ± 3 Ma, which is consistent with the ages of Early Palaeozoic ophiolitic complexes and ocean island sequences in the CMOB reported in previous studies. In contrast, metamorphic zircon from the retrograded eclogite samples yielded consistent Triassic metamorphic ages of 246 ± 2 and 245 ± 2 Ma, which can be interpreted as the timing of closure of the Palaeo‐Tethys Ocean. The compatible peak metamorphic mineral assemblages, P–T–t paths and metamorphic ages, as well as the similar protolith signatures for the eclogites in the CMOB and Longmu Co–Shuanghu suture (LCSS) suggest that the two belts formed part of a cold oceanic subduction system in the Triassic. The main suture zone of the Palaeo‐Tethyan domain extends at least 1,500 km in length from the CMOB to the LCSS in the Tibetan Plateau. The identification of lawsonite‐bearing retrograded eclogites in the CMOB provides important insights into the tectonic framework and complex geological evolution of the Palaeo‐Tethys.     

Subduction initiation in the Neo‐Tethys: constraints from counterclockwise P–T paths in amphibolite rocks of the Nagaland Ophiolite Complex,India

High‐P metamorphic rocks that are formed at the onset of oceanic subduction usually record a single cycle of subduction and exhumation along counterclockwise (CCW) P–T paths. Conceptual and thermo‐mechanical models, however, predict multiple burial–exhumation cycles, but direct observations of these from natural rocks are rare. In this study, we provide a new insight into this complexity of subduction channel dynamics from a fragment of Middle‐Late Jurassic Neo‐Tethys in the Nagaland Ophiolite Complex, northeastern India. Based on integrated textural, mineral compositional, metamorphic reaction history and geothermobarometric studies of a medium‐grade amphibolite tectonic unit within a serpentinite mélange, we establish two overprinting metamorphic cycles (M₁–M₂). These cycles with CCW P–T trajectories are part of a single tectonothermal event. We relate the M₁ metamorphic sequence to prograde burial and heating through greenschist and epidote blueschist facies to peak metamorphism, transitional between amphibolite and hornblende‐eclogite facies at 13.8 ± 2.6 kbar, 625 ± 45 °C (error 2σ values) and subsequent cooling and partial exhumation to greenschist facies. The M₂ metamorphic cycle reflects epidote blueschist facies prograde re‐burial of the partially exhumed M₁ cycle rocks to peak metamorphism at 14.4 ± 2 kbar, 540 ± 35 °C and their final exhumation to greenschist facies along a relatively cooler exhumation path. We interpret the M₁ metamorphism as the first evidence for initiation of subduction of the Neo‐Tethys from the eastern segment of the Indus‐Tsangpo suture zone. Reburial and final exhumation during M₂ are explained in terms of material transport in a large‐scale convective circulation system in the subduction channel as the latter evolves from a warm nascent to a cold and more mature stage of subduction. This Neo‐Tethys example suggests that multiple burial and exhumation cycles involving the first subducted oceanic crust may be more common than presently known.     

Coupled role of deformation and metamorphism in the construction of inverted metamorphic sequences: an example from far‐northwest Nepal

In the Greater Himalayan sequence of far northwestern Nepal, detailed mapping, thermobarometry, and microstructure analysis are used to test competing models of the construction of Himalayan inverted metamorphism. The inverted Greater Himalayan sequence, which is characterized by an increase in peak metamorphic temperatures up structural section from 580 to 720 °C, is divided into two tectonometamorphic domains. The lower domain contains garnet‐ to kyanite‐zone rocks whose peak metamorphic assemblages suggest a metamorphic field pressure gradient that increases up structural section from 8 to 11 kbar, and which developed during top‐to‐the‐south directed shearing. The upper portion of the Greater Himalayan sequence is composed of kyanite‐ and sillimanite‐zone migmatitic gneisses that contain a metamorphic pressure gradient that decreases up structural section from 10 to 5 kbar. The lower and upper portions of the Greater Himalayan sequence are separated by a metamorphic discontinuity that spatially coincides with the base of the lowest migmatite unit. Temperatures inferred from quartz recrystallization mechanisms and the opening angles of quartz c‐axis fabrics increase up section through the Greater Himalayan sequence from ～530 to >700 °C and yield similar results to peak metamorphic temperatures determined by thermometry. The observations from the Greater Himalayan sequence in far northwestern Nepal are consistent with numerical predictions of channel‐flow tectonic models, whereby the upper hinterland part evolved as a ductile southward tunnelling mid‐crustal channel and the lower foreland part ductily accreted in a critical‐taper system at the leading edge of the extruding channel. The boundary between the upper and lower portions of the Greater Himalayan sequence is shown to represent a foreland–hinterland transition zone that is used to reconcile the different proposed tectonic styles documented in western Nepal.     

The P–T–t evolution of the exhumed Himalayan metamorphic core in the Likhu Khola region,East Central Nepal

The recent identification of multiple strike‐parallel discontinuities within the exhumed Himalayan metamorphic core has helped revise the understanding of convergence accommodation processes within the former mid‐crust exposed in the Himalaya. Whilst the significance of these discontinuities to the overall development of the mountain belt is still being investigated, their identification and characterization has become important for potential correlations across regions, and for constraining the kinematic framework of the mid‐crust. The result of new phase equilibria modelling, trace element analysis and high‐precision Lu–Hf garnet dating of the metapelites from the Likhu Khola region in east central Nepal, combined with the previously published monazite petrochronology data confirms the presence of one of such cryptic thrust‐sense tectonometamorphic discontinuities within the lower portion of the exhumed metamorphic core and provides new constraints on the P–T estimates for that region. The location of the discontinuity is marked by an abrupt change in the nature of P–T–t paths of the rocks across it. The rocks in the footwall are characterized by a prograde burial P–T path with peak metamorphic conditions of ~660°C and ~9.5 kbar likely in the mid‐to‐late Miocene, which are overlain by the hanging wall rocks, that preserve retrograde P–T paths with P–T conditions of >700°C and ~7 kbar in the early Miocene. The occurrence of this thrust‐sense structure that separates rock units with unique metamorphic histories is compatible with orogenic models that identify a spatial and temporal transition from early midcrustal deformation and metamorphism in the deeper hinterland to later deformation and metamorphism towards the shallower foreland of the orogen. Moreover, these observations are comparable with those made across other discontinuities at similar structural levels along the Himalaya, confirming their importance as important orogen‐scale structures.     

Influence of tectonic overpressure on P–T paths of HP–UHP rocks in continental collision zones: thermomechanical modelling

The principle of lithostatic pressure is habitually used in metamorphic geology to calculate burial/exhumation depth from pressure given by geobarometry. However, pressure deviation from lithostatic, i.e. tectonic overpressure/underpressure due to deviatoric stress and deformation, is an intrinsic property of flow and fracture in all materials, including rocks under geological conditions. In order to investigate the influences of tectonic overpressure on metamorphic P–T paths, 2D numerical simulations of continental subduction/collision zones were conducted with variable brittle and ductile rheologies of the crust and mantle. The experiments suggest that several regions of significant tectonic overpressure and underpressure may develop inside the slab, in the subduction channel and within the overriding plate during continental collision. The main overpressure region that may influence the P–T paths of HP–UHP rocks is located in the bottom corner of the wedge‐like confined channel with the characteristic magnitude of pressure deviation on the order of ～0.3 GPa and 10–20% from the lithostatic values. The degree of confinement of the subduction channel is the key factor controlling this magnitude. Our models also suggest that subducted crustal rocks, which may not necessarily be exhumed, can be classified into three different groups: (i) UHP‐rocks subjected to significant (≥0.3 GPa) overpressure at intermediate subduction depth (50–70 km, P = 1.5–2.5 GPa) then underpressured at depth ≥100 km (P ≥ 3 GPa); (ii) HP‐rocks subjected to ≥0.3 GPa overpressure at peak P–T conditions reached at 50–70 km depth in the bottom corner of the wedge‐like confined subduction channel (P = 1.5–2.5 GPa); (iii) lower‐pressure rocks formed at shallower depths (≤40 km depth, P ≤ 1 GPa), which are not subjected to significant overpressure and/or underpressure.     

Metamorphic evolution of medium‐temperature ultra‐high pressure (MT‐UHP) eclogites from the South Dabie orogen,Central China: an insight from phase equilibria modelling

Medium‐temperature ultrahigh pressure (MT‐UHP) eclogites from the south Dabie orogen, as represented by samples from the Jinheqiao, Shuanghe and Bixiling areas, consist of garnet, omphacite, phengite, epidote, hornblendic amphibole, quartz/coesite and rutile with or without kyanite and talc. Garnet is mostly anhedral and unzoned, but a few porphyroblasts are weakly zoned with core–mantle increasing grossular (X_gr) and decreasing pyrope (X_py) contents. Garnet compositions are closely correlated with the bulk compositions. For instance, the X_py and X_gr contents are positively correlated with the bulk MgO and CaO contents. Phengite is occasionally zoned with core–rim deceasing Si content, and phengite grains as inclusions in garnet show higher Si than in the matrix, suggesting differently resetting during post‐peak stages. The maximum Si contents are mostly 3.60–3.63 p.f.u. for the three areas. Pseudosections calculated using THERMOCALC suggest that the MT‐UHP eclogites should have a peak assemblage of garnet + omphacite + lawsonite + phengite + coesite in most rocks of higher MgO content. In this assemblage, the X_py in garnet mostly depends on bulk compositions, whereas the X_gr in garnet and the Si contents in phengite regularly increase, respectively, as temperature and as pressure rise, and thus, can provide robust thermobarometric constraints. Using the X_gr and Si isopleths in pseudosections, the peak P–T conditions were estimated to be 40 kbar/730 °C for the Jinheqiao, 41 kbar/726 °C for the Shuanghe, and 37–52 kbar and 700–830 °C for the Bixiling eclogites. Some eclogites with higher FeO are predicted to have a peak assemblage of garnet + omphacite + coesite ± phengite without lawsonite, where the garnet and phengite compositions highly depend on bulk compositions and generally cannot give available thermobarometric constraints. Decompression of the eclogites with lawsonite in the peak stage is inferred to be accompanied with cooling and involves two stages: an early‐stage decompression is dominated by lawsonite dehydration, resulting in increase in the mode of anhydrous minerals, or further eclogitization, and formation of epidote porphyroblasts and kyanite‐bearing quartz veins in eclogite. As lawsonite dehydration can facilitate evolution of assemblages under fluid‐present conditions, it is difficult to recover real peak P–T conditions for UHP eclogites with lawsonite. This may be a reason why the P–T conditions estimated for eclogites using thermobarometers are mostly lower than those estimated for the coherent ultramafic rocks, and lower than those suggested from the inclusion assemblages in zircon from marble. A late‐stage decompression is dominated by formation of hornblendic amphibole and plagioclase with fluid infiltration. The lawsonite‐absent MT‐UHP eclogites have only experienced a decompression metamorphism corresponding to the later stage and generally lack the epidote overprinting.     

Metamorphic history of glaucophane‐paragonite‐zoisite eclogites from the Shanderman area,northern Iran

The Shanderman eclogites and related metamorphosed oceanic rocks mark the site of closure of the Palaeotethys ocean in northern Iran. The protolith of the eclogites was an oceanic tholeiitic basalt with MORB composition. Eclogite occurs within a serpentinite matrix, accompanied by mafic rocks resembling a dismembered ophiolite. The eclogitic mafic rocks record different stages of metamorphism during subduction and exhumation. Minerals formed during the prograde stages are preserved as inclusions in peak metamorphic garnet and omphacite. The rocks experienced blueschist facies metamorphism on their prograde path and were metamorphosed in eclogite facies at the peak of metamorphism. The peak metamorphic mineral paragenesis of the rocks is omphacite, garnet (pyrope‐rich), glaucophane, paragonite, zoisite and rutile. Based on textural relations, post‐peak stages can be divided into amphibolite and greenschist facies. Pressure and temperature estimates for eclogite facies minerals (peak of metamorphism) indicate 15–20 kbar at ~600 °C. The pre‐peak blueschist facies assemblage yields <11 kbar and 400–460 °C. The average pressure and temperature of the post‐peak amphibolite stage was 5–6 kbar, ~470 °C. The Shanderman eclogites were formed by subduction of Palaeotethys oceanic crust to a depth of no more than 75 km. Subduction was followed by collision between the Central Iran and Turan blocks, and then exhumation of the high pressure rocks in northern Iran.