Eclogites and other high-pressure rocks in the Himalaya: a review

Abstract Himalayan high-pressure metamorphic rocks are restricted to three environments: the suture zone; close to the suture zone; and (mostly) far (>100 km) from the suture zone. In the NW Himalaya and South Tibet, Cretaceous-age blueschists (glaucophane-, lawsonite- or carpholite-bearing schists) formed in the accretionary wedge of the subducting Neo-Tethys. Microdiamond and associated phases from suture-zone ophiolites (Luobusa and Nidar) are, however, unrelated to Himalayan subduction–collision processes. Deeply subducted and rapidly exhumed Indian Plate basement and cover rocks directly adjacent to the suture zone enclose eclogites of Eocene age, some coesite-bearing (Kaghan/Neelum and Tso Morari), formed from Permian Panjal Trap, continental-type, basaltic magmatic rocks. Eclogites with a granulite-facies overprint, yielding Oligocene–Miocene ages, occur in the anatectic cordierite ± sillimanite-grade Indian Plate mostly significantly south of the suture zone (Kharta/Ama Drime/Arun, north Sikkim and NW Bhutan) but also directly at the suture zone at Namche Barwa. The sequence carpholite-, coesite-, kyanite- and cordierite-bearing rocks of these different units demonstrates the transition from oceanic subduction to continental collision via continental subduction. The granulitized eclogites in anatectic gneisses preserve evidence of former thick crust as in other wide hot orogens, such as the European Variscides.

High-pressure metamorphic rocks such as blueschists (glaucophane ± lawsonite-bearing rocks), eclogites and high-pressure granulites are the products of major geodynamic processes such as subduction, collision and crustal thickening. More important, however, is the fact that metabasic blueschists, eclogites and high-pressure granulites often undergo much less deformation and recrystallization during exhumation to the present-day surface than their host rocks and thus provide crucial preserved evidence for these processes in the form of their metamorphic pressure-temperature histories. Extensive, more complete, information over the timing and extremes of metamorphic conditions for a region reduce significantly the uncertainties in thermomechanical models designed to reconstruct geodynamic processes, and thus allow a much more reliable prediction of how ongoing collision belts could respond and develop. Blueschists form in environments where temperature (T ) conditions are relatively low but pressures (P) are relatively high (i.e. high dP/dT )a situation typical for the upper levels of subduction zones, especially in accretionary prisms. The finding of blueschists in orogenic belts, especially when they are associated with ophiolites, is thus a key indicator of the location of former suture zones (i.e. plate boundaries). Eclogites also form at high pressures but require minimum pressures higher than those reached in normal continental crust such that plagioclase is no longer stable. In contrast to blueschists, eclogite can form at a wide range of temperatures and thus, if not completely retrogressed, can yield evidence of high pressures reached through the processes in crustal thickening as well as subduction. A third group of high-pressure rocks, highpressure granulite, is another indicator of extreme metamorphic conditions, conditions compatible with thickening or even short-lived subduction of the continental crust. Intensive investigation of orogenic belts in the last 30 years, aided by improvements in microanalysis methods such as micro-Raman spectroscopy, has allowed us to identify minerals such as coesite and microdiamond in rocks of subduction-collision orogens: minerals formed at mantle depths of 100 km or more (e.g. Schertl & O'Brien 2013). All of these rock types, blueschist, eclogite and high-pressure granulite, as well as rocks containing coesite and diamond occur within the Himalaya. How have they helped us to understand this orogeny?
The discovery of coesite in Kaghan Valley (Pakistan) eclogite (O'Brien et al. 1999) required existing tectonometamorphic models for the Indian Plate, essentially based on crustal thickening, to be significantly modified to include an early (Eocene) deep-crustal subduction and exhumation within the same, already short, timescale. Likewise, the discovery of granulitized eclogites in low-pressurehigh-temperature gneisses, far from the initial subduction zone and with a Miocene metamorphic age , has also caused a major rethink of collision models. Even more spectacular were finds of microdiamond and other ultrahighpressure minerals in remnants of Tethyan oceanic lithosphere trapped between the Indian and Asian plates (Robinson et al. 2004). Does this discovery require a further modification of existing tectonic models?
In this review of the high-pressure rocks of the Himalaya, the large-scale features and subdivision of the orogen will be summarized and the distribution of particular high-pressure rocks with respect to this subdivision outlined. Subsequently, the key geological, petrographical, petrological and geochronological features of the individual high-pressure locations will be presented. The final section will attempt to integrate knowledge gained from the high-pressure rocks with that from the vastly more abundant host series to allow a coherent tectonometamorphic model for the evolution of this spectacular orogenic belt to emerge.

Geological overview and simplified subdivision of the Himalaya
The closure of the Neo-Tethys Ocean and collision between the Indian and Asian plates and intervening arc terranes, starting in Eocene times and continuing to the present day, has produced the world's highest and most spectacular mountain chain, the Himalaya, as well as the largest area of high elevation, the Tibetan Plateau (Argand 1924; see the review by Searle 2015). On the Indian Plate, south of the suture zone, 10 peaks over 8000 m are located in the Himalaya, whereas four further 8000 m peaks are found at the margin of Asia north of the suture in the Karakoram. Fossiliferous marine sediments and segments of former ocean crust are found thousands of metres above sea level and far from any ocean. How can we explain this? In trying to reconstruct the Himalayan orogenic process as a whole from the metamorphic rocks presently exposed, it is necessary to separate the individual metamorphic stages and their respective environments.
The Himalayan-Tibetan orogen can be most easily subdivided into a number of roughly parallel belts, running for over 2500 km along strike of the collision zone, separated by major faults and thrusts ( Fig. 1) (Gansser 1964;Hodges 2000;Yin & Harrison 2000). These are, from north to south, the Trans-Himalayan Batholith, the Indus-Tsangpo Suture Zone, the Tethyan Zone, Higher Himalaya, Lesser Himalaya and Sub-Himalaya. Northward subduction of Neo-Tethys resulted in the development of an Andean-type, calk-alkaline magmatic arc of Upper Cretaceous-Eocene age, the Trans-Himalayan Batholith, running from Afghanistan to Burma (Honegger et al. 1982;Schärer et al. 1984;Wen et al. 2008). In the Lhasa terrane of Tibet this magmatism defines the Gangdese Belt, whereas further west, in India and Pakistan, identical mainly dioritic, tonalitic and granitic bodies have intruded the Kohistan-Ladakh island arc to produce the Ladakh (Honegger et al. 1982;Weinberg & Dunlap 2000) and Kohistan (Petterson & Windley 1985) batholiths. A last peak of magmatic activity at around 50 Ma and only very minor subsequent intrusion to around 46 Ma (Upadhyay et al. 2008;Wen et al. 2008) reflect the transition from oceanic to continental subduction as the Neo-Tethys Ocean closed completely.
The suture zone itself, roughly marking the course of two of the major rivers draining the southern Tibetan Plateau, is appropriately called the Indus-(Yarlung)-Tsangpo Suture Zone (ITSZ). It is well exposed in Ladakh (Fig. 2) where a series of deformed thrust sheets (Garzanti et al. 1987;Searle et al. 1990;Henderson et al. 2010) brings together shelf and slope rocks of the Indian Plate (Lamayuru Group), rocks of the oceanic realm (ophiolitic mélange, volcanic and volcaniclastic rocks of the Dras intra-oceanic arc) and rocks of the Asian margin (slope, shelf and shallow-marine sediments of the Tar and Jurutse groups). In the same area the youngest marine sediments (nummulitic limestone, 50.5 Ma: Green et al. 2008) and the exposed Ladakh granodiorite are directly overlain by terrigenous, fluvial conglomerates (containing cobbles from the batholith) and sandstones of the Indus Molasse (Frank et al. 1977;Brookfield & Andrews-Speed 1984). In many places along the suture zone, large outcrops of ultramafic-mafic ophiolite sequences, interpreted as back-arc basin oceanic crust of Jurassic-Cretaceous age, are a more obvious reminder that the suture zone marks the site of a former ocean (Bédard et al. 2009).
South of the ITSZ are rocks of the Indian Plate, imbricated, stacked and metamorphosed to different degrees, and separated by a series of roughly north-dipping thrusts and faults. In simple terms, the margin of the Indian Plate comprised Archean and Early Proterozoic gneisses, Lower and Mid-Proterozoic metasediments, and Neoproterozoic-Lower Paleozoic series covered by further Paleozoic-Cenozoic sedimentary series, with local basaltic (especially noteworthy being the Permo-Triassic Panjal traps) and granitic (Lower Proterozoic c. 1850 Ma and Cambro-Ordovician c. 480-500 Ma) magmatic events (e.g. Gansser 1964;DeCelles et al. 2001;Kohn et al. 2010;Gehrels et al. 2011). The continental crust entering the orogen has first undergone high-grade metamorphism, forms the present-day high topography and is known as the Higher (or Greater) Himalaya Crystalline Nappes (HHC). Crust incorporated into the thrust stack later, now tectonically below the HHC, is associated with lower topography, shows a dominantly lower grade of metamorphism (even in tectonic windows 100 km or more from the orogenic front) and is generally known as the Lesser Himalaya (LH). Emplacement of the hot HHC on the incoming LH along the ductile shear zone known as the Main Central Thrust (MCT) is associated with an inversion of metamorphic isograds: that is, prograde metamorphism structurally upwards in the LH, retrograde metamorphism downwards in the HHC (LeFort 1975;Stephenson et al. 2000).
For much of the arc of the Himalaya, Cambrian-Eocene marine platform sequences abut the ITSZ. These units, called Tethys Himalaya (TH), have become detached from the underlying Proterozoic basement as the Indian Plate was subducted, and thus, although strongly folded and imbricated, underwent mostly only very-low-to low-grade (up to greenschist-facies) metamorphism and preserve stratigraphic information as well as fossils (Gaetani & Garzanti 1991). The boundary with the underlying HHC later acted as a normal fault, the South Tibetan Detachment (STD), allowing exhumation of the HHC in an overall still compressional orogen (Burg & Chen 1984). Investigation of agedistribution patterns in detrital and magmatic zircon combined with Nd isotopic signatures, mainly in the central part of the Himalaya, has shown that the HHC and LH differ not just in grade of metamorphism but also in their stratigraphy (Parrish & Hodges 1996;Richards et al. 2005;Gehrels et al. 2011;Mottram et al. 2014;Spencer et al. 2018). The LH is dominated by a Paleoproterozoic series with ε Nd (0) values of −20 to −25, a crustal residence time of T DM > 2.0 Ga, and containing 1800 Ma orthogneisses and detrital zircons showing clusters at 1600-2600 Ma. In contrast, the HHC is characterized by Neoproterozoic-Cambrian sediments exhibiting ε Nd (0) values of −15 to −20, a crustal residence time of T DM < 2.0 Ga and additional detrital zircon groups in the range 500-1100 Ma. In fine detail, the parallel belts of the HHC, LH and TH mapped over the length of the Himalaya is an oversimplification, and in reality the correlation of stratigraphy, metamorphic grade and the nature of bounding thrusts, especially the MCT, is not so straightforward. Thus, in Pakistan and Tso Morari, where highgrade rocks generally mapped as the HHC directly abut the suture zone, the stratigraphy (mostly Phanerozoic) differs from that of the high-grade belt in the Central Himalaya (mostly Proterozoic) which is separated from the suture by the low-grade TH (see e.g. DiPietro & Pogue 2004).
Tectonically below the LH, separated by the Main Boundary Thrust, erosion products of the Himalaya are represented by Cenozoic marine and continental sedimentary rocks, deposited in the foreland basin to the south of the orogen, forming the Sub-Himalaya (Burbank et al. 1996). These rocks have been deformed and thrust southwards along the Main Frontal Thrust onto the Indian Craton over the youngest foreland basin sediments. Detailed study of the temporal evolution in detrital components in the Sub-Himalaya (e.g. Najman 2006) has allowed an insight into the erosion of the ITSZ and the pattern of exhumation of the HHC and LH.
High-pressure metamorphic rocks occur in a number of locations within the Himalayan belt but are restricted to three specific structural-tectonic environments, namely: the ITSZ; deeply subducted and quickly exhumed Indian Plate (HHC) close to the suture zone; and high-grade Indian Plate (HHC) in nappes south of the Tethyan Series far (>100 km) from the suture zone. In the NW Himalaya of India and Pakistan, Cretaceous blueschists and associated rocks (glaucophane-and/or lawsonite-and/ or carpholite-bearing schists) are reported from the accretionary wedge marking the subduction of Tethys beneath the Kohistan-Ladakh oceanic island arc. Further east, in the Sangsang region of south central Tibet, comparable blueschists formed in the subduction complex south of the continental arc along the Asian Plate margin. Ophiolites within the suture zone locally contain ultrahigh-pressure indicators (e.g. at Luobusa in Tibet and Nidar in Ladakh). The metamorphosed continental margin of the Indian Plate is evidenced close to this suture zone in the Upper Kaghan valleys of Pakistan and the Tso Morari area of Ladakh in India in the form of eclogites, some of which are coesite-bearing. In both cases, these eclogites formed from dykes, sills and flows of bimodal magmatic rocks located within granitic, pelitic and carbonate series of Paleozoic and Mesozoic age, thus indicating subduction of both the basement and cover units. Other highpressure rocks close to the suture zone indicating Indian Plate subduction are found at Stak (Nanga Parbat syntaxis, Pakistan), in the Lopu Range (south central Tibet) and in the Namche Barwa syntaxis (SE Tibet). A separate group of eclogites, strongly overprinted by granulite-facies assemblages, occur within the ortho-and paragneiss series of the Higher Himalayan Crystalline Nappes in the Kharta/Ama Drime area of Tibet and across the nearby border in NE Nepal, as well as in northern Sikkim and in NW Bhutan (i.e. significantly south of the suture zone).

Himalayan high-pressure rocks
High-pressure rocks in the ITSZ The ITSZ comprises a complex mixture of rocks from the margins of both the Indian and Asian plates (i.e. shelf, trench, slope and accretionary prism environments), as well as segments of the Neo-Tethys oceanic crust and remnants of intra-oceanic arc volcanic rocks. The majority of high-pressure rocks of the ITSZ occur as blueschists within the ophiolitic mélange of the former accretionary prism, such as in the well-known Shangla (Pakistan) and Sapi-Shergol (Ladakh) locations (Honegger et al. 1989). Another unit in the ITSZ yielding evidence for high-and ultrahigh-pressure processes is the ophiolite itself. East of the Himalaya, in NE India (Nagaland) and Burma, the margin of the Indian Plate is also marked with a lawsonite-blueschist-and eclogite-bearing ophiolite belt similar to the ITSZ (Ghose & Singh 1980;Chatterjee & Ghose 2010;Ao & Bhowmik 2014;Bhowmik & Ao 2015). This eastern margin of the Indian Plate is not the subject of this review but these studies are listed for completeness.
The extension of the suture zone further east in Pakistan is very poorly documented but a unit with calcareous and graphitic schists, phyllites, and marbles, the Banna-Parla Sapat, has been mapped between the mafic-ultramafic rocks of the Kohistan Arc and the high-grade units of the Kaghan Valley (Khan et al. 1995). Erosion products suspected as belonging to this unit were found at Dhumduma, north of Naran, where the Bharjali Da Nar joins the Kunhar River. In quartz nodules of graphitic metapelites, relicts of Fe-Mg carpholite, mostly altered to phengite and chlorite, have been identified indicating high-pressure-low-temperature conditions of 1.2-1.6 GPa and 380-410°C (Fig. 3) (Oberhänsli 2013). Recrystallization of this assemblage was at temperatures of up to 500°C as deduced from Raman graphite thermometry (Oberhänsli 2013).
Sangsang (Tibet). Apart from the blueschists in Pakistan and India, there are also reports of rocks of similar character and age in the suture zone, north of Everest, at Sangsang in southern Tibet (Xiao & Gao 1984). Here, thrust-bounded units of ophiolite and ophiolitic mélange with blueschistfacies (≤540°C) metabasic and metasedimentary rocks as part of the ITSZ have been described (Borneman et al. 2014). A geochronological investigation suggested a protolith age of c. 111 Ma (U-Pb zircon) and 65 Ma for high-pressure metamorphism ( 40 Ar/ 39 Ar amphibole) in a metavolcanic rock (Borneman et al. 2014). The ophiolite itself, as in many other cases in the ITSZ, is Early Cretaceous in age (130-120 Ma: Xia et al. 2008).
Ultrahigh-pressure indicators in serpentinites (India and Tibet). Unusual ultrahigh-pressure minerals have been identified within serpentinized peridotites of the ITSZ ophiolite ( Fig. 1) at Nidar, Ladakh (Das et al. 2015(Das et al. , 2017 and Luobusa, Tibet (Robinson et al. 2004;Yang et al. 2007). In both cases microdiamond occurs, as well as other phases (wadsleyite, clinoenstatite, moissanite and coesite) requiring pressures corresponding to mantle depths of 140 km or more. These serpentinized peridotites represent altered segments of deep mantle that has upwelled during seafloor spreading and ocean-crust formation. The ultrahigh-pressure minerals in these rocks, although occurring spatially related to accretionary prism blueschists and eclogites of subducted crustal origin, do not yield information about the Himalayan subduction-collision process. Instead, they are unusual relict phases preserving evidence of unrelated deep-mantle processes. Modern microanalytical tools have made these phases accessible to study and recent summaries have emphasized how common these relicts are in other ophiolites worldwide (Zhang et al. 2016).
High-pressure rocks: Indian Plate close to the suture Tso Morari (India). The first reliable record of eclogites in the Himalaya, and their interpretation as metamorphosed Panjal Trap lavas and dykes (corona dolerite microstructure), was from the Tso Morari  Liou (1971) and Schmidt (1995). Coe = Qtz: Bohlen & Boettcher (1982); Ab = Jd + Qtz: Holland (1980). (Berthelsen 1953) but an even earlier report from Haydon (1904) had already noted garnet in the southern Tso Morari basic dykes. This finding lay largely forgotten until the 1990s when renewed investigation, spurred on by the reports of eclogite in the Kaghan Valley of Pakistan (Pognante & Spencer 1991), also began to unearth a high-pressure history for Indian Plate rocks in Ladakh Guillot et al. 1997). Subsequently, several groups have published detailed petrographical, petrological and geochronological reports on the Tso Morari eclogites, although a surprisingly large number of these (e.g. de Sigoyer et al. 1997;Sachan et al. 1999 Berthelsen (1953, p. 367, fig. 13).

(Figs 1 & 4) area of Ladakh
Geologically, the crystalline rocks of the Tso Morari gneiss dome form a c. 100 km-long, 50 kmwide, NW-SE-orientated body (Fig. 4) emerging from below the ITSZ, along the Zildat detachment fault (Fig. 5b), at the NE boundary and dipping below the Tethyan Himalaya of Zanskar to the SW (Steck et al. 1998;Epard & Steck 2008). The dominant rocks of the dome are variably deformed augen-gneisses derived from Cambro-Ordovician (479 ± 2 Ma) porphyritic S-type granitethe Puga Gneiss of the Tso Morari nappe (Fig. 6a, b) (Girard & Bussy 1999). Locally, such as at Polokongka La, the granite is hardly deformed but mylonitized variants are identical in age (Steck et al. 1998;Girard & Bussy 1999). The granites were intruded into clastic and carbonate sediments of late Proterozoic-Cambrian age and the whole sequence was intruded by basic dykes, geochemically characterized as subalkaline, within-plate basaltic rocks attributed to the Permian Panjal traps which today form folded and stretched layers and boudins (Rao & Rai 2006;Jonnalagadda et al. 2017a). Eclogite-facies relicts occur within this unit. Structurally overlying nappe units, the Tetraogal (Permian and younger shelf deposits with bimodal volcanic rocks) and Mata-Nyimaling-Tsarap (Late Proterozoic-Mesozoic clastic and carbonate metasediments intruded by the Ordovician I-type Rupshu granite) nappes, and intervening serpentinites, metagabbros, chromites and metabasalts of the Karzok ophiolites, show no eclogite-facies relicts and are mostly of amphibolite to greenschist grade (Steck et al. 1998;Epard & Steck 2008).
The best-preserved eclogites (Fig. 6c) contain Grt + Omp + Phe + Coe/Qtz + Rt ± Ky ± Mgs but show multiple reaction stages as for the metasediments Mukherjee & Sachan 2001;Mukherjee et al. 2003;Wilke et al. 2015). A key feature is inclusion-rich (quartz, rutile, calcic and sodic-calcic amphibole, paragonite,  (Wilke et al. 2015). Major and trace element (e.g. REE) patterns of these garnets allowed the combined thermodynamic and trace element modelling of fractionated garnet growth first at low-temperature eclogite facies and then ultrahighpressure conditions, thus substantiating the absence of coesite in the garnet core and its restriction to Mg-rich garnet rims (Konrad-Schmolke et al. 2008).
The matrix to the garnet (Fig. 7c) shows the patchy development of multiple steps along the exhumation path. Important features are: the presence of glaucophane as an early retrograde phase that is itself then mostly consumed by later sodiccalcic and/or calcic amphiboles ; phengite as a peak phase; paragonite in garnet core or as a secondary reaction product : Wilke et al. 2015; and matrix magnesite as cores to later dolomite and calcite (Mukherjee et al. 2003;Wilke et al. 2015). Additionally, a renewed heating phase to produce ilmenite and pargasite further complicates the petrographical story St-Onge et al. 2013).
The first geothermobarometric results for the eclogites presented >2.0 ± 3 GPa and 550 ± 50°C for the peak (  Initial geochronological studies of the Tso Morari high-pressure-ultrahigh-pressure rocks presented ages of 55 ± 12 Ma (Lu-Hf Grt-Omp-WR metabasite), 55 ± 7 Ma (Sm-Nd Grt-Gla-WR metapelite) and 55 ± 17 Ma (U-Pb Aln) for the eclogite stage; 48 ± 2 Ma ( 40 Ar/ 39 Ar Phe) and 45 ± 4 Ma (Rb-Sr Phe-Ap-Rt) in metapelites together with 47 ± 11 Ma (Sm-Nd Grt-Amp-WR) in a metabasite interpreted as the amphibolite-facies overprint; and mica 40 Ar/ 39 Ar ages of around 30 Ma in deformed metapelites interpreted as final cooling (de Sigoyer et al. 2000). Despite the large analytical errors, a spectacular initial exhumation rate of >5 mm a −1 was proposed (de Sigoyer et al. 2000(de Sigoyer et al. , 2004. However, diffusion modelling to determine the timescale allowing preservation of the sharp compositional step between the inner and outer garnet domains indicated even higher rates (23-45 mm a −1 ) which, when fitted to the existing ages, would require <3 myr (from 45 to 48) (Fig. 9) for the whole eclogite to greenschist history (O'Brien & Sachan 2000;Massonne & O'Brien 2003). Some workers (e.g. Leech et al. 2007) continued to propose 55-45 Ma to encompass the eclogite to amphibolite evolution despite the limitations from preserved garnet zoning but subsequent zircon geochronology (Donaldson et al. 2013;St-Onge et al. 2013) has added strong support to the <50 Ma age for the eclogite stage. A further constraint on the thermal evolution is given by zircon and fission-track ages, indicating exhumation to about 10-15 km depth by 45-40 Ma (Schlup et al. 2003). Taking into account the c. 50 Ma age for the end of marine sedimentation in the ITSZ, this is certainly a phenomenal example of exhumation occurring as fast as subduction (e.g. Rubatto & Hermann 2001).
Kaghan and Neelum valleys (Pakistan). In Pakistan, Indian Plate rocks commonly attributed to the HHC (Fig. 10) in most places abut directly the maficultramafic root of the tilted Kohistan Arc (along what is known as the Main Mantle Thrust or MMT) with locally only small slivers of ITSZ rocks in-between. Unlike the arid, high plateau of the Tso Morari with its eclogite-bearing ortho-augengneiss, the dissected (by glacial valleys) former plateau (e.g. van der Beek et al. 2009) of the upper Kaghan and Neelum valleys is far more varied both in terms of vegetation and geology (Figs 10 & 11a). Here a Proterozoic basement comprising pelitic and psammitic metasediments and granitic orthogneisses is overlain by a Paleozoic, mostly clastic, metasedimentary sequence with minor marbles which also contains granitic orthogneisses. A second Late Paleozoic-Mesozoic carbonate-rich (calc-schist, dolomitic marble, diopside marble, mica schist and quartzite) cover sequence contains abundant lenses and bodies of metabasite and felsic intrusive rocks  Greco et al. 1989). Structurally, the area comprises a stack of isoclinally folded thrust sheets, folded again by east-westtrending open folds to give the present structure (Greco & Spencer 1993;Kaneko et al. 2003;Treloar et al. 2003). Eclogites and amphibolitized eclogites occur as dykes, sills and boudins (Fig. 11b-e) in basement orthogneisses, as well as in cover unit paragneisses and marbles (Pognante & Spencer 1991;Fontan et al. 2000;O'Brien et al. 2001;Kaneko et al. 2003: Wilke et al. 2010a. As in Tso Morari, geochemically, the metabasites are subalkaline-alkaline tholeiitic basaltic rocks of within-plate type Spencer et al. 1995;Wilke et al. 2010a;Rehman et al. 2014Rehman et al. , 2016 and belong to the Permian Panjal Trap bimodal basaltic-rhyolitic sequence, as also indicated by zircon ages (Rehman et al. 2016).
The metamorphic conditions (Fig. 14) for the formation of the Kaghan eclogite, derived from equilibria for the Grt + Omp + Phe + Coe/Qtz + Rt peak assemblage, scatter around 2.5-3.2 GPa and 700 ± 70°C O'Brien et al. 2001;Parrish et al. 2006;Wilke et al. 2010a). Retrograde resetting combined with patchy growth zoning makes the choice of equilibrium mineral composition difficult. The post-peak growth of glaucophane requires cooling to 580-630°C (at 1.0-1.7 GPa) during the initial exhumation (Wilke et al. 2010a). The subsequent amphibolite-facies overprint, the metamorphism dominant in the surrounding Grt + Ky-bearing gneisses (Treloar 1995) and reflected as Cpx + Pl ± Amph symplectites consuming omphacite, multiple varieties of green to brown calcic amphibole in different microstructural sites and ilmenite ± titanite replacing rutile indicates heating to 650-720°C (at 1.0-1.  porphyroblasts over the amphibolite-facies assemblages and the related folding at greenschist facies (400-500°C) marks the last major overprint.

Stak (Pakistan). Northeast of the Kaghan Valley,
where the Indus cuts through the Nanga Parbat-Haramosh Massif, eclogite and high-pressure granulite have been described within rapidly exhuming Indian Plate rocks directly at the contact with the Kohistan Arc at the Stak location (LeFort et al. 1997). Mafic rocks occur as small lenses and boudins within the felsic gneisses forming part of the Paleozoic and Mesozoic cover to basement gneisses, just as in the upper HHC further west. Although eclogite garnet and omphacite compositions and zoning are comparable to that of certain Kaghan eclogites (Lanari et al. 2013;Kouketsu et al. 2016), the rocks lack low-temperature retrograde glaucophane and/or barroisite, and instead contain abundant high-temperature secondary brown hornblende. This difference is reflected in deduced P-T conditions (Fig. 14) of 2.5 GPa and 750°C for the peak eclogite stage followed by exhumation through the highpressure granulite and amphibolite field (Cpx + Pl symplectites and amphiboles) at 0.9-1.6 GPa at 650-700°C (Lanari et al. 2013). Geochemically, the metabasites appear in their major and trace element, as well as isotope patterns very similar to Panjal Trap and Kaghan eclogitized Panjal Trap basalts (Kouketsu et al. 2017). The eclogites contain small (<0.05 mm), irregularly shaped zircons, devoid of oscillatory zoning, some of which yielded a sensitive high-resolution ion microprobe (SHRIMP) U-Pb lower intercept age of 32.0 ± 2.6 Ma, whereas other analyses, typically also showing significantly higher Th/U, yielded 206 Pb/ 238 U ages of 93-158 Ma (Kouketsu et al. 2016). Due to the small size of zircons, some ages represent a mix of younger and older domains in grains but clearly a Mesozoic protolith and Cenozoic metamorphic overgrowth are present. The relationship between the c. 32 Ma age and the eclogite-facies stage is not clear. It is possible that symplectitic replacement of eclogite-facies omphacite led to zircon overgrowth but, regardless of the zircon growth process, it is clear that this segment of the HHC stayed longer at high-grade conditions than the nearby Kaghan crustal segment. Further high-pressure mafic granulites (1.0 GPa and 700°C), probably belonging to the 1.6 Ga dyke swarm in the basement (Treloar et al. 2000), were described from within the kyanite-bearing gneisses of the HHC further north in the Stak Valley . A correlation between these granulites and the granulitized eclogites in the Eastern Himalaya (Namche Barwa, Ama Drime, Bhutan) remains speculative.
Lopu Range (Tibet). In the suture zone of southern Tibet, in the Lopu Range roughly north of Kathmandu (Fig. 1), Indian Plate rocks that experienced Eocene age high-pressure metamorphism are exposed. As elsewhere in the Central Himalaya, ophiolite and molasse of the ITSZ occur between the Gangdese Arc and low-grade Phanerozoic TH rocks but a top-to-the-north normal-sense shear zone exposes high-pressure metasediments in its footwall within this sequence (Laskowski et al. 2016). Petrological and geochronological investigations reveal conditions of >1.4 GPa and ≤600°C for the metamorphism of various Phe + Chl ± Sta ± Grt-bearing metasediments, an age of 40.4 ± 1.4 (Lu-Hf Grt-WR) for garnet, and 40 Ar/ 39 Ar phengite ages of between 34 and 39 Ma (Laskowski et al. 2016). Detrital zircon age spectra from these rocks indicate a correlation with Tethyan sequences (Laskowski et al. 2017) and thus, as in Tso Morari and the Kaghan/Neelum valleys, continental India and its sedimentary cover were both subducted and metamorphosed at the start of the orogeny with at least a small fragment quickly returned to shallow levels, thus preserving the petrological evidence.
Namche Barwa (Tibet). At the easternmost end of the main Himalayan arc, centred around the peaks of Namche Barwa and Gyala Peri (7782 and 7151 m, respectively) in Tibet (Fig. 16), an anticlinal dome of HHC protrudes through the Asian margin, comprising Gangdese Belt granites intruded into Phanerozoic sediments and their basement, very much like the situation of Nanga Parbat protruding through the Kohistan Arc far to the west in Pakistan (Burg et al. 1998). A narrow strip of ophiolitic mafic and mafic rocks, thought to be equivalent to those of the ITSZ, separates these continental units and is associated with metasediments (quartzites, metapelites, calc-silicate rocks), attributed to the TH, that are also strongly deformed (Zhang et al. 1992;Geng et al. 2006). The crystalline rocks of the Namche Barwa core are predominantly anatectic gneisses with local garnet-rich felsic and mafic granulites, and a discrete unit containing marble and calcsilicate marble: both Proterozoic and Paleozoic components have been identified (Zhang et al. 1992).
High-pressure rocks: Indian Plate far from the suture In the same year that coesite eclogite was discovered in the Kaghan Valley, Lombardo et al. (1998) reported for the first time eclogites from the Eastern Himalaya with petrographical and petrological details published soon after . These newly discovered eclogites, strongly granulitized but recognizable as former higherpressure rocks, were found east of the Everest-Makalu region, in the fault-bounded Ama Drime Range, close to the Phung Chu River (Arun River when it crosses into Nepal) and near the village of Kharta. Subsequent studies in the same region, just south across the border into Nepal (Cottle et al. 2009;Corrie et al. 2016), in north Sikkim (Rolfo et al. 2008) and in NW Bhutan (Warren et al. , 2012Regis et al. 2014), also resulted in new finds of relict eclogites (Fig. 18).
In Tibet, the metabasic rocks occur as boudinaged dykes and lenses, decimetre to metre sized,  Cottle et al. 2009). These rocks are thus of continental, rift-related origin but, in contrast to the Panjal Trap basalts of the NW Himalayan eclogites, of Proterozoic rather than Permo-Triassic age.
Detailed petrographical descriptions and photomicrographs in the literature concerning the Eastern Himalaya granulitized eclogites note especially a fine-grained symplectitic intergrowth of Jd-poor clinopyroxene and plagioclase (oligoclase), interpreted as a replacement of former omphacite Liu et al. 2007;Cottle et al. 2009;Warren et al. 2011). Matrix omphacite has not been found but exists as rare inclusions in garnet and zircon (Y.H. . The symplectites locally contain isolated amphibole or orthopyroxene but commonly these domains are surrounded by massive amphibole. Garnet is dominantly Alm-Gr-Prp in the ratio 3:2:1, unzoned, in some samples rich in inclusions (Qtz + Rt), and commonly with reactions at rims that produced orthopyroxene + plagioclase (andesine). In more heavily overprinted samples, garnet sits in moats of secondary plagioclase and thus eclogite peak compositions have probably been lost. Former phengite has been replaced by orientated biotite in plagioclase Groppo et al. 2007;Y.H. Wang et al. 2017) and ilmenite (later also titanite) has replaced most of the rutile in the matrix.
Geochronological studies of the granulitized eclogites, their host gneisses and structures linked to their exhumation have been undertaken by several groups. Zircon U-Pb ages of 13.9 ± 1.2, 14.9 ± 0.7  (Corrie et al. 2016) appears to substantiate this young (Miocene-age) eclogitefacies metamorphism. However, garnet Lu-Hf ages of 37.5 ± 0.8, 36 ± 1.9 and 33.9 ± 0.8 Ma from Ama Drime metabasite samples, yielding additionally U-Pb zircon ages of 13-15 Ma, strongly suggest an Eocene age for eclogite-facies metamorphism and a Miocene age for the granulite-facies overprint (Kellett et al. 2014). Ages of monazites in gneisses hosting the granulitized eclogites record well the high-temperature stage (11-15 Ma: Cottle et al.   2009; Warren et al. 2011Warren et al. , 2012, which was followed shortly after by a period of rapid exhumation and cooling of the granulite-facies HHC from the Kharta/Ama Drime across to NW Bhutan related to normal faulting, especially related to the STD (Kali et al. 2010;Grujic et al. 2011;Kellett et al. 2013).

Discussion
In older Phanerozoic eclogite-bearing, subductioncollision orogens, such as the Caledonides and Variscides in Europe, the mountains have been deeply eroded, the crust and lithosphere has returned to a stable thickness, and dynamic tectonic and metamorphic processes are no longer active. In contrast, the Himalayan region, and other equivalent Cenozoic orogens such as the European Alps, are still undergoing active shortening as reflected in the frequency and magnitude of earthquakes, are typified by high mountains and relief contrast, and have overthickened crust reaching 70 km or more. Looking at the present-day structure (Fig. 19, stage 4) of the Himalayan-Tibetan region as visualized by large-scale geophysical programmes (e.g. INDEPTH: Zhao & Nelson 1993;Nelson et al. 1996;and HIMPROBE: Jain et al. 2012;Caldwell et al. 2013), thickened continental crust is certainly an important feature: a continental crust that is at eclogite-and high-pressure granulite-facies conditions today. However, this large-scale structure is a product of collision processes that have been taking place for 50 myr and are still ongoing. Can we use this present-day information to explain the highpressure rocks at the surface today?
The ages and associated rocks of these different high-pressure units within the Himalaya nicely demonstrate the general tectonic history of this subduction-collision orogen: that is, the transition from oceanic subduction to continental collision via a short-lived deep-continental subduction stage (Fig. 19). The presence of lawsonite-bearing blueschists (Groppo et al. 2016) and carpholite-bearing metasediments (Oberhänsli 2013) indicate lowtemperature (i.e. 'cold') subduction but also an exhumation regime that also remained cool or was sufficiently fast to enable transport to shallow crustal levels. This type of tectonometamorphic regime is typical of Pacific-type convergent margins (see the review by Tsujimori & Ernst 2014), although extensive carpholite-bearing series exist in the European Alps and along the Mesozoic collision belt in Turkey (Oberhänsli 2013). This environment represents the stage of subduction of Neo-Tethys under the Asia margin, and the processes taking place in ocean crust-derived volcanosedimentary units within the accretionary prism. The Cretaceous age of blueschist minerals ties in well to the age of ophiolite formation and obduction (Bédard et al. 2009). The blueschists are in a piggyback position, tectonically above major intracrustal thrusts, and with hot anatectic crustal rock sequences to the north and south. In older collision orogens, where only eroded deeper crustal levels are exposed, it is unlikely that this upper tectonic level would have survived, and perhaps only serpentinized peridotites hint at the former presence of a suture zone.
As the last of the oceanic crust was subducted, the initial continental margin was also deeply subducted (Fig. 19, stage 2). The record of this stage is marked by the eclogites in the Kaghan/Neelum valleys in Pakistan and the Tso Morari area of the NW Himalaya O'Brien et al. 2001). The distribution of eclogites in the Proterozoic basement, and Paleozoic metasediments and granitic orthogneisses indicate that the whole of the Indian Plate crust was subducted, thus contrasting with the tectonic separation of basement (involved in intracrustal stacking and emplacement below the Asia margin) and Tethyan cover (deformed but remaining south of the suture) in the central part of the orogeny. The Lopu area is a minor exception to this pattern that requires further investigation. Coesite in the Kaghan and Tso Morari eclogites and gneisses strongly suggests crustal subduction depths of at least 90 km but geochronological data support only a very short time for both the subduction and the majority of the exhumation: that is, exhumation occurred at subduction rates of several centimetres per year (Massonne & O'Brien 2003;Wilke et al. 2010b). The presence of glaucophane post-dating the metamorphic pressure peak in both areas supports exhumation at lower temperatures, as expected for rocks moving upwards between the mantle wedge and the down-going slab. Continental crust is significantly less dense than the mantle, even if partially transformed to high-pressure minerals, thus buoyancy forces would lead to intracrustal shear-zone development with the consequence that detached crustal slivers would force their way back up the subduction channel, at least to typical (i.e. c. 30 km) crustal depths (Chemenda et al. 2000;O'Brien 2001;Treloar et al. 2003). The presence of these ultrahigh-pressure eclogite units now directly adjacent to the ITSZ strongly supports exhumation in the subduction channel (as well as a steep angle of subduction). In the Kaghan area, the HHC is not a single coherent unit but a series of at least three nappes with coesite eclogite structurally uppermost and basic magmatic rocks only partially altered to corona dolerite in the lowest nappe. Again, support for subduction-channel stacking of the same crustal unit that had reached different depths. The location of ultrahigh-pressure eclogites and the rate of subduction and exhumation may seem spectacular but is perfectly compatible with the results of thermomechanical modelling of subducted crust assuming suitable input parameters (e.g. Warren et al. 2008;Beaumont et al. 2009).
Although the majority of petrological studies of the Kaghan and Tso Morari coesite eclogites concluded a hairpin-like P-T path indicating the initial exhumation as near isothermal or accompanied by cooling (de Sigoyer et al. 2000;O'Brien et al. 2001;Mukherjee et al. 2003;Wilke et al. 2010aWilke et al. , b, 2015St-Onge et al. 2013;Palin et al. 2014Palin et al. , 2017, alternative interpretations do exist. The work of Chatterjee & Jagoutz (2015), in contrast to all previous studies, proposed exhumation of the Tso Morari eclogite by diapiric ascent through the overlying mantle wedge, and also interprets the hightemperature overprint at amphibolite facies proposed by most authors, to have been at significantly higher pressures of the eclogite facies. However, first-order observations such as core-rim compositions of inclusion and matrix omphacite, the absence of lawsonite, the presence of coesite exclusively in the Mg-rich zones of overgrowth garnet, exhumation initially producing glaucophane, and the modelling of trace-element and inclusion patterns in garnet (e.g. de Sigoyer et al. 1997;Konrad-Schmolke et al. 2008;Wilke et al. 2015) do not support this alternative. Once all these problems are resolved the P-T path for diapiric exhumation through the mantle disappears and a deep, fast subduction and fast exhumation history, including a marked heating at around 1.0 GPa, emerges: a path directly comparable, also in its timing, to that of the coesite eclogite of the Kaghan Valley, Pakistan.
It is fortunate that the section of Indian Plate containing the Kaghan and Tso Morari eclogites reached shallow levels such that they were located in the hanging wall to the subsequent low-angle emplacement of Indian Plate crust below Asia. What would have happened to rocks that stayed deeper? The Stak eclogite resembles closely the Kaghan eclogites but has suffered a greater degree of overprint, indicating that it was not so well protected from thermal effects during and after exhumation (Lanari et al. 2013). Metamorphic zircons in this rock are also significantly younger than those of the Kaghan eclogite (32 Ma compared to 47 Ma: Kouketsu et al. 2016). Is this a reflection of more zircon being produced during the symplectite breakdown of omphacite in the Stak sample as opposed to ultrahigh-pressure metamorphic zircon being sampled in fresh Kaghan eclogite, or is there a real difference in the age of eclogite formation? The same question has to be asked of the eclogites and high-pressure rocks in the Eastern Himalaya (Ama Drime/Sikkim/Bhutan/Namche Barwa). A critical question to ask is whether these rocks were metamorphosed early in the orogeny, stayed at depth for 20-30 myr, heated during crustal relaxation and partial exhumation to granulite facies, then were extruded in migmatite-rich units to shallow structural levels? Such a model is consistent with the geochronological results of Kellett et al. (2014), whereby Lu-Hf ages of garnet (Eocene) are clearly separated from those of zircon (Miocene) in the same rock. The linking of zircon growth to exact stages in the metamorphic history is a difficult but critical aspect of interpreting such high-grade terranes and is, for good reason, still a major research area (e.g. Harley et al. 2007;Rubatto 2017). Interestingly, in metasediments of the Eastern Himalaya, monazite ages document several other stages in Oligocene-Miocene times making them unusually more informative, geochronologically, than the metabasites despite their weaker petrographical memory and further supporting an extended high-grade metamorphic evolution (Chakraborty et al. 2016;J.M. Wang et al. 2017).
The detailed models designed to reconstruct the crustal architecture of the Himalaya and its metamorphic evolution (e.g. Beaumont et al. 2001Beaumont et al. , 2004Beaumont et al. , 2006Jamieson et al. 2004) based on the large-scale features deduced from the geophysical studies show several features that could help to answer these questions. Once the short-lived subduction and exhumation of the Indian continental margin was over, continued convergence of India and Asia led to the subduction of India Plate material below Asia, the development of a foreland-propagating series of shallow-angle thrusts leading to a widening of the orogen, and thickening of the lithosphere (to produce kyanite-grade rocks) both to the north and the south of the suture zone (Fig. 19, stage 3). In a crust now dominated by stacked slivers of upper-crustal material with respective upper-crustal radioactive heat production, heating produced sillimanite-grade gneisses overprinting earlier kyanite-bearing assemblages. The high-grade rocks would be partially molten, weak and thus liable to flow if squeezed enough. This type of process was invoked to explain the extrusion of anatectic, sillimanite-grade HHC below a detaching roof (Searle & Rex 1989;Beaumont et al. 2004Beaumont et al. , 2006Searle et al. 2006) (Fig. 19, stage 4) but others (e.g. Kohn 2008Kohn , 2014 still question whether a steady-state critical taper rather than channel flow better explains the thermal pattern. The important point is that the stacked and thickened crust of the past, as that of today, was at eclogite-facies pressures at its base. However, as temperature increased, melting occurred and rock masses started to flow to be emplaced at shallower levels, initiated by the 'indentor' effect of incoming Indian lithosphere (Warren et al. 2008); it is highly likely that the dominant metasedimentary gneisses would be unable to retain mineralogical evidence of the high-pressure stage apart from minor kyanite relicts. In contrast, metabasic eclogites in migmatite terranes commonly allow insight into their past history even when multiply overprinted (e.g. Scott et al. 2013).
The discovery of granulitized eclogites within migmatitic gneiss complexes is well known from other orogens. For example, in the Moldanubian Zone of the European Variscides, numerous metabasites with eclogite-facies omphacite and garnet relics partially overprinted by granulite-facies (augitic clinopyroxene, orthopyroxene, spinel or even Fe-rich olivine) assemblages are enclosed by anatectic cordierite + sillimanite-bearing metapelites (O'Brien & Vrána 1995;O'Brien 1997;Scott et al. 2013). In contrast to the Himalayan granulitized eclogites, Variscan examples are closely associated with serpentinized spinel peridotite, thus yielding strong evidence for a tectonic interleaving of metabasiteperidotite and sedimentary units (i.e. of crust and mantle). Regardless of whether the Kharta/Ama Drime, Sikkim or NW Bhutan granulitized eclogite locations are considered, there are no mantle slices and no convincing evidence that metabasic rocks have been tectonically emplaced into the HHC para-and orthogneisses. A logical interpretation of this situation is that the eclogitization and granulitization of the metabasites was in situ: that is, that these particular nappe units of the Higher Himalaya also experienced eclogite-facies conditions as a whole. In contrast to the Tso Morari and Kaghan areas of the HHC directly adjacent to the suture zone, these rocks were not so deeply subducted and lack a rapid exhumation to cold, shallow levels; instead, staying longer at depth, undergoing significant heating and being exhumed as anatectic bodies to become emplaced over other Indian Plate units that were never at eclogite-facies conditions. Results from modelling of the Himalayan-type crustal collision (e.g. Beaumont et al. 2004Beaumont et al. , 2006Jamieson et al. 2007) predicts a long duration for high-grade metamorphism before rocks become extruded. These models incorporated erosion at a mountain front to aid extrusion but the wide, hot crustal belt in the Himalaya and southern Tibet also fits the characteristics of modelled Grenville Orogen crust (Jamieson et al. 2007) which shows breakouts of hot, deep levels into shallow, cooler levels independent of surface erosion. It is clear from the increasing data from the field (e.g. Goscombe et al. 2006;Carosi et al. 2010;Chakungal et al. 2010;Grujic et al. 2011;Warren et al. 2011Warren et al. , 2012Montomoli et al. 2015) that the HHC is not a single coherent slab but different units, with different histories, separated by ductile shear zones. This allows a more reasonable explanation for the granulitized eclogites yielding Miocene ages; namely, that these rocks were a deeper, eclogite-facies part of the crustal stack that became extruded in the Miocene. It has been speculated that the reason for the extrusion of an anatectic, deeper part of the thickened Himalayan nappe stack is the incorporation to the orogen of a cold, rigid 'plunger' of Indian Plate material at this time ('hard' India-Asia collision): a period also associated with marked slowing of the India-Asia convergence (Warren et al. 2008;Molnar & Stock 2009;Grujic et al. 2011;van Hinsbergen et al. 2012). The cluster of ages, dominantly from zircon and monazite, is most likely to have been due to the widespread granulite-facies overprint which led to the breakdown of former higherpressure phases: that is, the end of the stay in the lower part of the thickened crust is documented by datable minerals in the metabasite, and not the start and duration of the thickened-crust episode. In detail, however, monazite ages in the metapelites indicate a long period of residence at high temperatures (or, less likely, multiple metamorphic episodes), whereas the metabasites, although rich in information petrologically, unusually have a poorer geochronological memory (Cottle et al. 2009;Grujic et al. 2011;Warren et al. 2011Warren et al. , 2012Chakraborty et al. 2016;J.M. Wang et al. 2017). The extrusion of granulitized eclogites and their emplacement on lower-grade units all within the HHC may be the best example of all for channel flow.

Summary and conclusions
Despite the predominance of petrological studies in the Himalaya on metasedimentary rocks, it is perhaps once again the minor metabasic rocks that have a reliable memory of the full extent of the metamorphic development and, by extrapolation, thus also the tectonic history of Earth's most spectacular mountain belt. The minor lawsonite-, glaucophaneand even carpholite-bearing rocks testify to the cold subduction environment predating collision. Despite predominant amphibolite-and even greenschist-facies assemblages in Indian Plate basement-cover sequences in the Kaghan/Neelum valleys and Tso Morari, it is the minor preserved eclogite cores to mafic boudins, especially those containing coesite, that elucidate the full extent of the deep, rapid subduction and equally rapid exhumation of the leading edge of India as Neo-Tethys was finally consumed. In addition, despite sillimanite ± cordierite ± spinel formation, along with copious melt, in HHC metasediments and orthogneisses in the Eastern Himalaya, once again it is the resilience of the metabasites to complete resetting that allows the picture of an early eclogite-facies stage to emerge. This last eclogite stage is perfectly consistent with the geophysical evidence of an overthickened crust still present today. With collision still ongoing after 50 myr, it would not be surprising if further pulses of extrusion from deep levels occurtaking a closer look at older, deeply eroded orogens such as the Variscan, Caledonian, Grenville or Trans-Hudson might yield some indicators as to how this will look (Johansson & Möller 1986;Massonne & O'Brien 2003;Jamieson et al. 2010;St-Onge et al. 2006).