JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 PAGES 1673–1719 2000 Contrasting P–T Paths in the Eastern Himalaya, Nepal: Inverted Isograds in a Paired Metamorphic Mountain Belt BEN GOSCOMBE∗ AND MARTIN HAND GROUP FOR INTEGRATED TECTONIC STUDIES AT DEPARTMENT OF GEOLOGY AND GEOPHYSICS, ADELAIDE UNIVERSITY, ADELAIDE, S.A. 5005, AUSTRALIA RECEIVED NOVEMBER 12, 1998; REVISED TYPESCRIPT ACCEPTED APRIL 17, 2000 Petrology and phase equilibria of rocks from two profiles in Eastern Nepal from the Lesser Himalayan Sequences, across the Main Central Thrust Zone and into the Greater Himalayan Sequences reveal a Paired Metamorphic Mountain Belt (PMMB) composed of two thrust-bound metamorphic terranes of contrasting metamorphic style. At the higher structural level, the Greater Himalayan Sequences experienced high-T/moderate-P metamorphism, with an anticlockwise P–T path. Low-P inclusion assemblages of quartz + hercynitic spinel + sillimanite have been overgrown by peak metamorphic garnet + cordierite + sillimanite assemblages that equilibrated at 837 ± 59°C and 6·7 ± 1·0 kbar. Matrix minerals are overprinted by numerous metamorphic reaction textures that document isobaric cooling and re-equilibrated samples preserve evidence of cooling to 600 ± 45°C at 5·7 ±1·1 kbar. Below the Main Central Thrust, the Lesser Himalayan Sequences are a continuous (though inverted) Barrovian sequence of high-P/moderate-T metamorphic rocks. Metamorphic zones upwards from the lowest structural levels in the south are: Zone A: albite + chlorite + muscovite ± biotite; Zone B: albite + chlorite + muscovite + biotite + garnet; Zone C: albite + muscovite + biotite + garnet ± chlorite; Zone D: oligoclase + muscovite + biotite + garnet ± kyanite; Zone E: oligoclase + muscovite + biotite + garnet + staurolite + kyanite; Zone F: bytownite + biotite + garnet + K-feldspar + kyanite ± muscovite; Zone G: bytownite + biotite + garnet + K-feldspar + sillimanite + melt ± kyanite. The Lesser Himalayan Sequences show evidence for a clockwise P–T path. Peak-P conditions from mineral cores average 10·0 ± 1·2 kbar and 557 ± 39°C, and peak-metamorphic conditions from rims average 8·8 ± 1·1 kbar and 609 ± 42°C in Zones ∗Corresponding author. Present address: 18 Cambridge Road, Aldgate, Adelaide, S.A. 5154, Australia. E-mail: [email protected] Extended data set can be found at: http://www.petrology.oupjournals.org D–F. Matrix assemblages are overprinted by decompression reaction textures, and in Zones F and G progress into the sillimanite field. The two terranes were brought into juxtaposition during formation of sillimanite–biotite ± gedrite foliation seams (S3) formed at conditions of 674 ± 33°C and 5·7 ± 1·1 kbar. The contrasting average geothermal gradients and P–T paths of these two metamorphic terranes suggest they make up a PMMB. The upperplate position of the Greater Himalayan Sequences produced an anticlockwise P–T path, with the high average geothermal gradient being possibly due to high radiogenic element content in this terrane. In contrast, the lower-plate Lesser Himalayan Sequences were deeply buried, metamorphosed in a clockwise P–T path and display inverted isograds as a result of progressive ductile overthrusting of the hot Greater Himalayan Sequences during prograde metamorphism. thermobarometry; P–T paths; Himalaya; metamorphism; inverted isograds; paired metamorphic belts KEY WORDS: INTRODUCTION Previous work in the Himalayas interpreted all metamorphic mineral assemblages to have formed in the Himalayan Metamorphic Cycle during collisional orogenesis in the Tertiary (Pecher, 1989; Vannay & Hodges, 1996). Although one of the youngest metamorphic belts, use of the Himalaya as a type example of tectonic and metamorphic processes in collisional orogens is limited because the tectonometamorphic evolution is complex and not entirely understood. This study of the Eastern Oxford University Press 2000 JOURNAL OF PETROLOGY VOLUME 41 Himalaya recognizes two thrust-bound units of contrasting metamorphic style (termed metamorphic ‘terranes’): (1) the Lesser Himalayan Sequences (LHS), which underwent Barrovian high-P/moderate-T metamorphism and display a clockwise P–T path; (2) at higher structural levels, the high-T/moderate-P Greater Himalayan Sequences (GHS), which display an anticlockwise P–T path. Matrix mineral assemblages in both terranes were formed during SSW- to S-directed ductile overthrusting during collisional orogenesis. Thus, the two terranes constitute a Paired Metamorphic Mountain Belt (PMMB), in a continent–continent collisional orogen (Armstrong et al., 1992), distinct from subduction-type paired metamorphic belts in the sense of Miyashiro (1973). Anticlockwise P–T paths and high average geothermal gradients have not been reported from elsewhere in the Himalayan orogen and may be restricted to a portion of the Eastern Himalaya only. East of the Mount Everest–Dudh Kosi area in Eastern Nepal (Fig. 1), very few metamorphic studies have been previously undertaken. All are from the Makalu–Arun River area (Brunel & Kienast, 1986: Lombardo et al., 1993; Meier & Hiltner, 1993; Pognante & Benna, 1993), with none available for the Kangchenjunga region. These previous metamorphic studies document metamorphic assemblages in the different structural units and attribute these to crystallization events on the one clockwise P–T path. In an attempt to characterize the metamorphism of the wider Eastern Nepal region, both with petrological analysis and P–T calculations, two profiles across the LHS and GHS were mapped and sampled by the author in 1992 and 1997 (Fig. 1). These are the Makalu profile along the Arun River and upper Barun Khola (Fig. 2) and the Kangchenjunga profile along the Tamur River and Kabeli Khola (Fig. 3). This study reports new petrological data and a comparable set of P–T estimates using one thermobarometric method [THERMOCALC; Powell & Holland (1988)], these providing a framework for the tectonometamorphic evolution of the Eastern Nepal region. The distribution of metamorphic assemblages does not differ significantly from that given by previous workers. However, in contrast, the present study reports significant differences in interpretative P–T paths and peak metamorphic conditions for both the LHS (P >3 kbar higher) and GHS (T >50–100°C higher). Variation in metamorphic conditions through the Eastern Nepal profiles is shown to be discontinuous, and is interpreted as reflecting two metamorphic terranes with contrasting P–T evolution, although formed within the same metamorphic cycle. REGIONAL GEOLOGY The entire Himalayan front is a metamorphic belt that was metamorphosed and exhumed in the Tertiary NUMBER 12 DECEMBER 2000 Himalayan Orogeny as a result of collision of the Indian subcontinent with Asia. The Himalayan metamorphic belt has been subjected to only one metamorphic cycle. Along the entire length of the Himalayan front, metamorphic field gradients are inverted, with metamorphic grade increasing towards higher structural levels (e.g. LeFort, 1975; Hodges et al., 1988; Pecher, 1989; Inger & Harris, 1992; Vannay & Hodges, 1996). Isograds are parallel to the main metamorphic foliation and vary smoothly across the Main Central Thrust (MCT) Zone (Table 1), without a break in grade (Hubbard, 1989; Metcalfe, 1993; Hubbard, 1996). Thus metamorphic equilibration is interpreted to have been synchronous with ductile, southward thrusting along the MCT (Pecher, 1989). Numerous models for the origin of the inverted isograds in the Himalayas have been discussed in detail by Mohan et al. (1989), Searle & Rex (1989) and Hubbard (1996). Some workers (Hodges et al., 1988; Pecher, 1989; Inger & Harris, 1992; Vannay & Hodges, 1996) proposed that the Himalayan Metamorphic Cycle is polyphase, with two metamorphic events (Table 2). Those workers saw evidence for an Eohimalayan metamorphism spanning 37–24 Ma age (Inger & Harris, 1992; Hodges et al., 1994; Parrish & Hodges, 1996; Coleman & Hodges, 1998). This high-P/moderate-T Barrovian event is thought to be associated with burial and prograde metamorphism (Inger & Harris, 1992; Vannay & Hodges, 1996). Neohimalayan metamorphism at 22–13 Ma (Hubbard & Harrison, 1989; Inger & Harris, 1992; Metcalfe, 1993) is interpreted as a moderate-P/moderate-T peak metamorphic event (Hodges et al., 1988) associated with upthrusting of the GHS (Table 1) along the MCT (Inger & Harris, 1992; Hodges et al., 1996; Vannay & Hodges, 1996). The Himalayan metamorphic belt is composed of laterally extensive sequences bounded by crustal-scale north-dipping thrusts. The basal unit, the LHS (Table 1), is interpreted as Middle Proterozoic in age (Parrish & Hodges, 1996) or alternatively as Upper Proterozoic to Palaeozoic Tethyan meta-sediments (Vannay & Hodges, 1996). The Main Boundary Thrust forms the basal contact, with the LHS thrust onto the Late Miocene molasse of the Siwalik Group (LeFort, 1975). The LHS is overthrust by the GHS at the MCT. The GHS is continuously exposed along the Himalayan front and interpreted as a metamorphosed sequence of Late Proterozoic supracrustal rocks and orthogneisses (Bhanot et al., 1977: Miller & Frank, 1992; Parrish & Hodges, 1996), along with Late Cambrian granitoid rocks (LeFort et al., 1986). Although the GHS is thought to have originally been overthrust by the Tethyan Sequences of Upper Proterozoic to Eocene age, the two terranes are separated by a series of extensional detachment faults, the Southern Tibet Detachment System (STDS) (Table 1; Fig. 1) (Burg 1674 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Fig. 1. Simplified geology of the Eastern Himalaya in Nepal, with location of Makalu (Fig. 2) and Kangchenjunga profiles (Fig. 3) indicated. (a) Tibetan series sediments. Proterozoic to early Palaeozoic Greater Himalayan Sequences (GHS) composed of (c) Barun Gneiss and Jannu–Kangchenjunga Gneiss, (d) Namche Migmatite Orthogneiss, (e) Black Gneisses and (f ) Miocene leucogranite. Uppermost Lesser Himalayan Sequences (LHS) include (b) Himal Group and (g) Kathmandu Group with (h) Num Orthogneiss Formation and (k) Cambrian granites. Palaeozoic Tethyan rocks of the LHS in the Midlands Group are (i) Kushma Formation, ( j) Ulleri Formation and (l) Seti Formation. Below the Main Boundary Thrust (MBT) is (m) Molasse of the Siwalik Formation. STDS, Southern Tibet Detachment System; ITSZ, Intra-Tibetan Suture Zone. Other major thrusts are labelled MCTI and MCTII after Maruo & Kizaki (1981) (MCT, Main Central Thrust). Geology sourced from Bordet (1961), Shrestha et al. (1984, 1985), Lombardo et al. (1993), Morrison & Oliver (1993) and authors. Regional setting inset after Harris et al. (1993). E, Mount Everest; M, Makalu; K, Kangchenjunga. et al., 1984; Burchfiel et al., 1992; Edwards & Harrison, 1997). The sequence of thrusting and reactivation of crustalscale shear zones during collision is complex. All Tethyan oceanic crust had been subducted and the Indian subcontinent had collided with Asia by Middle Eocene times (Sengor, 1990). Prograde burial of the GHS and LHS occurred during the Middle Eocene to Oligocene by overthrusting of Tethyan Sequences at the Eohimalayan Thrust, the precursor to the STDS (Vannay & Hodges, 1996). Peak metamorphic conditions in the Early Miocene resulted in emplacement of leucogranites of 26–17 Ma age into the highest structural levels of the GHS (Scharer, 1984; Deniel et al., 1987; Copeland et al., 1990; Hodges et al., 1996, 1998; Edwards & Harrison, 1997; Coleman, 1998). The MCT was initiated at >22 Ma (Harrison et al., 1995; Hodges et al., 1996, 1998; Edwards & Harrison, 1997; Coleman, 1998) and remained active throughout metamorphism, resulting in the inverted metamorphic sequence (Hubbard & Harrison, 1989; Inger & Harris, 1992). Metamorphism immediately below the MCT may be as young as 6–13 Ma (Macfarlane et al., 1992; Harrison et al., 1997). Isostatic uplift and extensional collapse of the High Himalayas was accommodated by episodic movement on extensional detachments at structural levels above the MCT, at the STDS, from >22 Ma (Burchfiel et al., 1992; Hodges et al., 1996, 1998; Vannay & Hodges, 1996; Edwards & Harrison, 1997; Harrison et al., 1999). In the region investigated, the High Himalayas of Eastern Nepal, the GHS is composed of a variety of distinct gneissic rock units. In the Makalu profile (Fig. 2); the base of the GHS is marked by the principal thrust of the MCT zone, locally called the Barun Thrust. In the Kangchenjunga profile the MCT is represented by the highly sheared informal ‘Biotite Gneiss’ unit. The overthrust Barun Gneiss (Bordet, 1961; Lombardo et al., 1993) is composed of high-grade quartzo-feldspathic gneisses, mafic gneiss and calc-silicates. At higher structural levels is the Namche Migmatite Orthogneiss and above that the Black Gneisses (Lombardo et al., 1993). In the Kangchenjunga profile (Fig. 3); the entire GHS is called the Jannu–Kangchenjunga Gneiss (Mohan et al., 1989). The Jannu–Kangchenjunga Gneiss is 1675 JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 Fig. 2. Simplified geological map of Makalu profile and location of samples (prefixed by M in tables). Geology based on Bordet (1961), Lombardo et al. (1993), Pognante & Benna (1993) and authors. 1676 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Fig. 3. Simplified geological map of Kangchenjunga profile and location of samples (prefixed by K in tables). Geology based on Shrestha et al. (1984) and authors. Table 1: Abbreviations PMMB Paired Metamorphic Mountain Belt; paired metamorphic belt in a continent–continent collisional orogen environment, distinct from subduction-type paired metamorphic belts of Miyashiro (1973) STDS South Tibet Detachment System MCT Main Central Thrust; ductile thrust zone within the upper LHS; principal movement plane at base of GHS (the Barun Trust and Biotite Gneiss unit in the eastern Nepal region investigated) MBT Main Boundary Thrust; basal thrust to Himalayan sequences GHS Greater Himalayan Sequences, also called Higher Himalayan Crystallines; Proterozoic sequences from the north Indian passive margin, including Barun Gneiss, Black Gneisses, Namche Migmatite Orthogneiss and Jannu–Kangchenjunga Gneiss units in the eastern Nepal region LHS Lesser Himalayan Sequences; Palaeozoic sequences mostly of Tethyan origins, including Himal, Kathmandu and Midlands Groups in the eastern Nepal region petrologically similar to the Barun Gneiss of the Makalu profile and both are correlated with each other. At high structural levels, both the Black Gneisses and Jannu– Kangchenjunga Gneiss contain leucogranite intrusions of 22–24 Ma age (Scharer, 1984). Below the Barun Thrust are the Barrovian metamorphic rocks of the LHS composed of numerous groups and formations (Bordet, 1961; Lombardo et al., 1993) (Figs 1 and 2). The uppermost are the undifferentiated Himal Group (Shrestha et al., 1985) and below that the Num Orthogneiss and Kathmandu Group (Shrestha et al., 1985; Lombardo et al., 1993) (Figs 1–3). Lowest structural levels of the LHS are exposed in antiformal tectonic windows (Figs 1–3). These are thrust-bound units of Tethyan meta-sediments of the Midland Group, including the Kushma, Ulleri and Seti Formations (Shrestha et al., 1985). STRUCTURAL GEOLOGY OF THE MAKALU AND KANGCHENJUNGA PROFILES The lowest structural level of the LHS is a series of thrust slices of Seti, Ulleri and Kushma Formations of the Midland Group (Shrestha et al., 1985), exposed in an antiformal dome in the Kangchenjunga profile (Fig. 3). The Ulleri Formation is particularly highly sheared and 1677 JOURNAL OF PETROLOGY VOLUME 41 almost entirely mylonitic. In the Makalu profile, the principal thrust zone of the MCT and upper bounding margin of the LHS is a series of thrust sheets within the Himal Group (Bordet, 1961; Lombardo et al., 1993) (Fig. 2). In the Kangchenjunga profile the principal thrust zone of the MCT is represented by the 900 m thick, informal Biotite Gneiss unit (Fig. 3). Biotite Gneiss is a highly sheared unit of biotite-rich schists and mylonitic NUMBER 12 DECEMBER 2000 gneisses with numerous pegmatitic and granitic veins and tight to isoclinal folds. The movement history within the Biotite Gneiss unit is unknown. Throughout the entire LHS a strong proto-mylonitic and schistose metamorphic foliation (S1) is developed parallel to layering and thrusts, dipping shallowly to the NE to NW (Fig. 4). L1 stretching lineation is defined by quartz- and feldspar-aggregate ribbons and aligned 1678 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA kyanite, biotite and hornblende grains. L1 is shallowly (20–30°) north-plunging throughout most of the LHS and in the Kangchenjunga profile is NE-plunging in the Himal Group (Fig. 4). S–C relationships, shearbands and -type porphyroclasts throughout indicate north over south transport, consistent with shearing in the Himalayan Metamorphic Cycle (Table 2). Aligned micas, amphiboles and kyanite within the matrix assemblage suggest that ductile shearing immediately preceded or accompanied the peak of metamorphism. S1 foliation is crenulated, with new biotite growth, and is isoclinally to closely folded around sub-horizontal east–west-trending axes, normal to L1. F2 folds and crenulations verge to the south and formed in response to the same tectonic transport as the S1–L1 fabric. Leucocratic partial melt segregations occur only at the highest structural levels in the Himal Group, as thin stromatic layers and elongate moats parallel to L1 and enclosing garnet or tourmaline porphyroblasts. Above the MCT, in contrast to the LHS, the GHS is a less schistose sequence of high-grade migmatitic gneissic rock-types. In the Makalu profile the basal GHS unit is the Barun Gneiss, composed of metapelitic, quartzofeldspathic and mafic gneisses with subordinate calcsilicates and carbonates. The overlying Namche Migmatite Orthogneiss (Lombardo et al., 1993) is structurally and petrologically similar to the Barun Gneiss and also contains subordinate carbonate, calc-silicate and mafic layers; differing only in a greater proportion of stromatic partial melt segregations (Pognante & Benna, 1993). The overlying Black Gneisses are quartzo-feldspathic gneisses that are more schistose than the remainder of the GHS (Lombardo et al., 1993), but are otherwise structurally similar. The nature of the contact between the Namche Migmatite Orthogneiss and the Black Gneisses is unknown. In the Kangchenjunga profile, the Jannu– Kangchenjunga Gneiss (Mohan et al., 1989) unit is petrologically and lithologically similar to and occupies the same structural level as the Barun Gneiss (Fig. 1), with which it is tentatively correlated. In all units of the GHS, the majority of partial melt segregations (typically garnet bearing) pre-date the metamorphic foliation (S1) and less common biotite pegmatitic veins post-date S1. Black Gneisses and the highest structural levels of the Jannu–Kangchenjunga Gneiss are intruded by numerous leucogranite plutons and sills. In the Makalu profile, gneisses of the Barun Gneiss have a coarse polygonal granoblastic texture with a strong S1–L1 fabric defined by aligned micas, sillimanite, hornblende and quartz- and feldspar-aggregate ribbons. This high-grade shear fabric constitutes the coarse matrix assemblage that was annealed at the peak of metamorphism (Table 2). S1 dips moderately NW and two stretching lineation directions are apparent, plunging shallowly to the NNE and west (Fig. 4). The NNEplunging lineation is defined by aligned peak-metamorphic sillimanite prisms and mineral aggregate ribbons. The -type porphyroclasts and asymmetrical isoclinal folding document NNE over SSW transport accompanying the dominant S1–L1 shear fabric. Consequently, this fabric is interpreted to have formed coeval with the S1–L1 schistose foliation in the underlying LHS. East-plunging lineations are unique to the Barun Gneisses, not being recognized in the LHS (Figs 2 and 4). This lineation orientation is possibly due to refolding by tight to isoclinal F2 folds with shallow NNE to NE plunges. In the Jannu–Kangchenjunga Gneiss S1 is nearly flat lying (Fig. 4), of variable intensity and stretching lineations are only rarely apparent. S1 is overprinted at low angles, by thin (millimetrescale), fine-grained sillimanite + biotite ± gedrite shearbands (S3) with NE- to ENE-plunging sillimanite and quartz-aggregate lineations. S3 is common in the Barun Gneiss and Jannu–Kangchenjunga Gneiss and less common within the uppermost Himal Group in the Kangchenjunga profile. The S3 foliation is interpreted to be the latest episode of ductile shearing associated with the MCT Zone. Sense of S3 shear is unknown in the area investigated. This episode may correspond to the Neohimalayan event (Vannay & Hodges, 1996) (Table 2). Tight to close mesoscopic F4 folds with shallow NNE to NE plunge and crenulation cleavages with aligned biotite are recognized in both the GHS and LHS. These may be coeval with the large-scale folding giving rise to the antiformal tectonic windows within the LHS in the Eastern Himalaya (Figs 1 and 3). DESCRIPTION OF ROCKS ALONG THE MAKALU PROFILE Greater Himalaya Sequences Black Gneisses The Black Gneisses are composed of biotite-rich schists and schistose quartzo-feldspathic gneisses and are intruded by the tourmaline + biotite ± garnet-bearing Makalu leucogranite. Schists have well-aligned matrix assemblages of biotite + muscovite + quartz + albite + K-feldspar + ilmenite ± sillimanite ± tourmaline (Table 3). Quartzites and foliated quartzo-feldspathic gneisses have similar assemblages but with lower modal proportions of micas and are often rich in tourmaline. Typically biotite and muscovite are in equilibrium, although coarse muscovites both overgrow biotite and are boudinaged with interstitial biotite growth. Muscovite and feldspars both contain fibrolite inclusions and muscovite is overprinted by thin fibrolite seams. Chlorite is typically retrogressive after biotite and muscovite. 1679 JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 Fig. 4. Lower hemisphere equal-area projections of structural data from both the Makalu and Kangchenjunga profiles. Ε, mineral and mineralaggregate lineations (L1); Φ, tight to isoclinal fold axes (F2); Χ, poles to bedding and the sub-parallel metamorphic foliation (S1). Namche Migmatite Orthogneiss The Namche Migmatite Orthogneiss is a mixed unit composed predominantly of augen orthogneiss and highly migmatized quartzo-feldspathic gneisses with subordinate calc-silicate, carbonate and mafic gneiss. Both quartzo-feldspathic and metapelitic gneisses have biotite + K-feldspar + plagioclase + quartz + ilmenite + sillimanite matrix assemblages (Table 3). Some metapelites also contain garnet, perthite, antiperthite and coarse-grained muscovites that overgrow biotite and K-feldspar. Like Barun Gneiss samples, the coarse foliated matrix assemblage is overprinted by thin, fine-grained S3 seams of fibrolite + biotite + muscovite + opaque minerals. Matrix ilmenite is 1680 (foli-QFG) (foli-QFG) (QFG) (quartzite) M73 M74 M70 M68∗ gn X-S S S S S S SK bi 1681 X (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) M11 M13 M16 M17 M18 M21 M22 M34 M38 M39∗ X X XS X X X X X X Barun Gneiss, metapelites (pelite) M78∗ (QFG) (pelite) M75 MX R SL X SR X S S SR S mu IgnXS C IgnXS C IgnXS IgnXS XS Ign S R IgnXS IgnXS XS IgnXS S IgnXS F S-L R F XS XS Namche Migmatite Orthogneiss (schist) (schist) M71∗ (schist) M69∗ M72∗ (schist) M60 Black Gneisses Sample st XS XS X X X XS X XS X XS X X XS X? X X? S S XS kf XS XS X X XS XS X IgnXS X XS X X XS X X X X X-S S S pl XSC S S XS Ign S S XS IgnXSC XS XSF Iq S S S SCF Imu,q S SK sill ky Table 3: Petrology of aluminous rocks from the Makalu profile Ign sp Ign X IgnXS XS SE Ichl XR R Rmu Imu SR SR ilm XS IgnXS Ign X Ign X XS IgnXS Ign X XS Ign S IgnXS XS X Ign S IgnXS Isill, gn IgnXS IgnXS X XS XS X X X X Imu S S SR S q E E Ibi XF R X R E hm XS E E mg R S R SR R SR R R SR R chl XC X X X cd L L L L L gd gn→green bi gn→bi, hm→mg, bi→chl cd→sill-gd sill→fine sill, ilm→hm±mg gn→bi, sill→fine sill cd→gd, cd→sill, gn→sill gd→fine sill sill→fine sill, kf→mu, bi→mu bi→ilm, ilm→hm mu→bi bi→mu, sill→chl, ilm→hm bi→mu, mu→fibr, bi→chl±ilm mu→ilm bi→chl bi→chl Reaction textures act gn→bi-pl gn→q-cd, ap, zr, mon, gn→bi-sill-cd, graphite zr, rut zr, rut ap, tm sph act, ap, tm tm sph tm tm Accessory GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) M47ii M49 M51b∗ M54 M81 M82 M83 M85 1682 M92i X X X X X X X XC Ign S IgnXS XS S S Ign X XS XS IgnXS Ign S C IgnXS bi R R S mu (QFG) (QFG) (QFG) (QFG) (QFG) (QFG) M14 M15 M52 M58∗ M92ii M95 X X IgnXS IgnXS X S IgnXS R S XCR X XS X X X X X XS XS X XS IXS XG X XG XC XS kf S sill X IgnXS XS X X XS IX XS XS X XS IgnXS XG X XCS XS IgnXS Imu Ccd XSC XS S XS XS XS XS S X G S C IgnXS XC XS pl ky Ign sp X IgnXS ilm XS IgnXS Ign X X Ign X IgnXS Ign X IgnXS IgnXS XS XS IgnXS XG IgnXS IXS XS X X Ibi S IgnXS Ign X-S XS XS SC IgnXGS X-S Ign X IgnXS q X C E X-S hm C X-S mg R R R R R R R chl X IgnXC X XC X X cd L L S S L gd ilm→hm sill→ilm, bi→chl bi→chl cd→chl gn→sill-chl, ilm→bi→pl-q←gn gn→chl-mg, gn→bi-pl±sill ilm→bi→myr→gn, Reaction textures zr, tm sph, rut zr tm zr, rut zr, rut, py rut ap, zr, sph, zr, ap sill→gd sill→fine sill, sill, q→g→q→gn cd→sill, sill→fine gn→cd gn→sill-bi-q, pl→mu cd→phengite-sill, cd→sill, gn→cd ilm→mg-hm gn→q→gn gn→chl gn→green bi, bi→chl, gn→chl ap, phengite gn→phengite, zr, rut ap zr zr rut ap, zr, cc, py Accessory NUMBER 12 X X X st VOLUME 41 Barun Gneiss, quartzo-feldspathic gneisses (pelite) M47i X (pelite) (pelite) M44 gn Sample Table 3: continued JOURNAL OF PETROLOGY DECEMBER 2000 (foli-QFG) X M95a∗ 1683 X X X (quartzite) M99 (schist) M111 X S Ign S XS X S Ist S S S S IgnXS S bi S S S XR SL S S S XSC XS R mu X X X X PHmu st Epl S S kf X X X S X XS XS X pl S-L sill X X X present X XS XS ky sp S ilm XS IgnXS S XR IstXS S C hm C Ign,chl X S S Ibi Ign,st S Ign,st S Ist S S XS X Ist,gn X Ist,gn XS XS Ign, pl Ign XS q mg S S R S R S R SR R chl cd gd bi→mu ky→mu, bi→chl, Reaction textures tm ap piemontite ap ap ap ap,tm tm mu across chl gn→chl-bi-q bi→q-chl-ilm, ky→q-mu bi→chl, pl→mu, mu across bi st→chl±mu gn→chl, bi→chl gn→chl-hm mu→chl, st→mu→chl, ky→seri←bi tm, sph, rut ky→chl, ky→mu ap, zr Accessory Samples grouped according to the mapped units (Fig. 2) and arranged from north to south down the page. I, inclusion phase with host indicated; X, peak metamorphic phase; S, aligned or contained within the dominant fabric in equilibrium with peak metamorphic phases; X-S, coarse aligned matrix; C, corona phase; K, crenulation phase; G, sub-grains; F, fine foliated seam; L, late phase; R, retrogressive phase; E, exsolution phase; PH, pseudomorphed phase; ∗, analysed sample. Postscript indicates phase that is retrogressed or contains the inclusion. Reaction textures presented as reactant phases → product phases. (schist) M109∗ Himal Group, Zone C (schist) M97 Himal Group, Zone D (schist) M110 (schist) M106∗ (schist) (schist) M104 M108∗ (schist) M98∗ Himal Group, Zone E (schist) X gn M2∗ Himal Group, Zone F Sample GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA JOURNAL OF PETROLOGY VOLUME 41 overgrown by haematite coronas. Mafic gneisses have quartz–plagioclase + K-feldspar + hornblende + garnet + biotite + sphene + ilmenite assemblages (Table 4). Barun Gneiss The Barun Gneiss is dominated by layered quartzofeldspathic paragneiss and metapelitic gneisses with numerous thin bands of mafic gneiss, carbonate and calcsilicate. Quartzo-feldspathic paragneisses and metapelites have identical assemblages and differ only by higher proportions of feldspars and less garnet and sillimanite in the quartzo-feldspathic gneisses. The coarse-grained, foliated granoblastic matrix consists of garnet + plagioclase + K-feldspar + quartz + sillimanite + biotite + ilmenite ± rutile ± cordierite (Table 3) and most samples contain partial melt segregations. The matrix assemblage comprises the dominant coarse-grained foliation (S1), which is not significantly reworked apart from minor fine-grained S3 foliation seams (Fig. 5d). Plagioclase is oligoclase to andesine and rarely antiperthitic, and Kfeldspar is microcline and commonly perthitic. Kyanite and staurolite are absent and muscovite is rare and invariably a late-stage phase, though aligned with S1 (Table 3). Magnetite is rare and haematite commonly occurs within S1 or as exsolution lamellae within primary ilmenite grains, but is absent from inclusion assemblages. In contrast to the LHS, garnet is commonly nearly devoid of inclusions, which occur only in cores and are absent from rims. Inclusions are typically biotite + quartz + ilmenite ± rutile and rarely also sillimanite and hercynitic spinel (Fig. 5a) (Table 3). Spinel occurs only as an inclusion phase and is absent from the matrix. Fine sillimanite inclusions preserve an earlier, overgrown crenulation cleavage. Matrix sillimanite is medium grained, in textural equilibrium with and aligned with the S1 matrix assemblage. Non-corrosive, fine-grained sillimanite overgrowths occur on the margins of primary sillimanite (Fig. 5d) and garnet. Garnet is consumed by quartz + biotite + sillimanite ± plagioclase symplectites (Fig. 5b) and overgrown by concentric quartz and garnet coronas. Cordierite mostly occurs as stretched augen in the matrix assemblage and rarely as cordierite ± quartz or cordierite + sillimanite + biotite symplectites that consume garnet porphyroblast margins. Fine-grained sillimanite + gedrite ± biotite and sillimanite + phengite aggregates both replace cordierite and comprise the S3 foliation seams that are sub-parallel to, but overprint the coarse-grained S1 fabric (Fig. 5d). Mafic gneisses have a well-aligned granoblastic matrix composed of quartz + labradorite + green hornblende + phlogopite + ilmenite; more rarely clinopyroxene or garnet or both also occur (Table 4). Retrogressive blue–green hornblende, chlorite or biotite occur on clinopyroxene margins. Garnet is rarely enclosed by coronas NUMBER 12 DECEMBER 2000 of clinopyroxene or quartz + hornblende symplectites. Calc-silicate gneisses have a granoblastic matrix of quartz + bytownite + meionite + clinopyroxene + orange phlogopite + sphene + ilmenite ± green hornblende ± calcite (Table 4). Hornblende-free samples contain forsterite. Matrix microcline is common and garnet is rare. Epidote and clinozoisite are interstitial and appear to be post-peak metamorphic phases. Forsterite also occurs as coronas on meionite. Carbonates have very similar matrix assemblages but are devoid of primary hornblende, meionite is less common and wollastonite is present in meionite-free samples. Wollastonite occurs as a matrix phase and as coronas on clinopyroxene. Lesser Himalayan Sequences Himal Group The Himal Group consists of thrust slices of aluminous schists, schistose gneisses and meta-quartzites with subordinate calc-silicate and mafic gneisses. Aluminous schists and meta-quartzites have coarse foliated matrix assemblages of quartz + albite + biotite + muscovite + ilmenite + kyanite ± garnet ± staurolite (Figs 5c and 6a and b) (Table 3). Chlorite and muscovite are typically in equilibrium with the S1 matrix assemblage but also occur as late-stage grains that overgrow S1 (Fig. 6a). Garnet and staurolite are often poikiloblastic (Fig. 6b) and typically have idiomorphic, inclusion-free overgrowths. Inclusions in both garnet and staurolite consist of quartz, biotite and ilmenite. Matrix sillimanite, andalusite, cordierite and hercynitic spinel are entirely absent from the Himal Group. Melt segregations are also absent except thin leucocratic segregations in gneiss sample (M95a) at the highest structural level (Zone F, see below). Kyanite and staurolite margins are retrogressively replaced by fine muscovite, often with an outer corona of chlorite. Metamorphic zones cannot be as accurately mapped as in the Kangchenjunga profile. However, there is a general progression of increasing metamorphic grade up sequence from the Arun River, north to the Barun Thrust at the base of the Barun Gneiss (Fig. 2). The following metamorphic zones are recognized: Zone C: garnet+biotite+muscovite+albite+chlorite; Zone D: garnet+biotite+muscovite+albite+kyanite; Zone E: garnet+biotite+muscovite+oligoclase+ staurolite+kyanite; Zone F: garnet+biotite+oligoclase+K-feldspar+ kyanite±muscovite (late sillimanite). Calc-silicate gneisses have a granoblastic matrix of plagioclase + meionite + grossular + clinopyroxene 1684 q pl (mafic) X (cpx layer) (pl layer) (mafic) M84∗ M84∗ M87∗ XG X X X X X X X X Ign X X X X Ipl,bi XG 1685 (calc-silicate) (calc-silicate) (calc-silicate) (calc-silicate) (bi-rich layer) X (calc-silicate) (calc-silicate) (calc-silicate) M50 M80 M28 M24 M24 M57∗ M57∗ M79 X X X X X X Ign X X X XG X X X X X Ign X X X X (calc-silicate) (calc-silicate) M48 X X X X M42∗ (calc-silicate) (calc-silicate) (calc-silicate) M27 M37a (calc-silicate) M26 M35 (calc-silicate) M25 Barun Gneiss, calc-silicate gneisses (mafic) (hn layer) M84∗ (mafic) M52 M56 (mafic) M51 Barun Gneiss, mafic gneisses M77 Namche Migmatite Orthogneiss Sample X X X X X X X X X XG X X kf X X X X X X X X X X X me X Ipl X X X X X gn X X X X X X X X X X X X Ipl,gn X X C IX cpx XR R R R R X XR X XR R XR C X XR Ign X X-S X X X hn X X X X ol X-L X XC X-L X X-L X ep X X X X X X X X R cc R C R act/tr X X X X Ign X X X X X X X X X X sph Table 4: Petrology of mafic, calc-silicate and carbonate rocks from the Makalu profile X X X X X XC XC X X X Ipl X X XR X Ign X X-S X-S Ipl X X-S X R R R R bi/phl chl X X X X X X X X X X X X X X X ilm X R X hm X woll piedmontite zr rut, ap tm ap ap rut ap ap Accessory fo→ep cpx→hn cpx→hn cpx→hn, ilm→mg cpx→hn phl→act cpx→phl, cpx→act, cpx→hn, cpx→phl ilm→mg, cpx→hn cpx→hn cpx→hn cpx→hn cpx→hn cpx→chl bi→chl, gn→chl, cpx→bi cpx→hn gn→q-hn, gn→cpx bi→chl, ilm→hn Reaction textures GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA q 1686 (mafic) M6 (mafic) X X X X X X X X gn X R XR hn XC X X X Ign X X R X X X X X X X X cpx X X XC ol X X C X X X X X X-L ep X X X X X X X X X cc C R R R act/tr X X Icpx X X X X X X X X X X X sph X X X X X X X R R R bi/phl chl X X X X X ilm X X hm X XC X woll ap, zr rut ap ap ma ma Accessory ep→cpx cpx→hn, gn→ep cpx→act cpx→hn, cpx→chl fo→hn±bi cpx→woll, woll→tr ilm→mg woll→q-tr, woll→act, pl→ma cpx→woll, cpx→act, bi→ma me→fo Reaction textures Samples grouped according to the mapped units (Fig. 2) and arranged from north to south down the page. I, inclusion phase with host indicated; X, peak metamorphic phase; S, aligned or contained within the dominant fabric in equilibrium with peak metamorphic phases; X-S, coarse aligned matrix; C, corona phase; L, late phase; R, retrogressive phase; G, sub-grains or deformed; ∗, analysed sample. Postscript indicates phase that is retrogressed or contains the inclusion. Reaction textures presented as reactant phases → product phases. M103 X X X X X X X me NUMBER 12 Himal Group, Zones C–E (calc-silicate) X M6 X (calc-silicate) Ign, cpx Ign X X X (calc-silicate) X X X M4∗ X X X X X X X kf M1∗ X X X X X X X X pl VOLUME 41 Himal Group, Zone F (carbonate) (carbonate) M57 M53 (carbonate) M43 (carbonate) (carbonate) M41 M79 (carbonate) M37 (carbonate) (carbonate) M33 M79 (carbonate) M24 Barun Gneiss, carbonates Sample Table 4: continued JOURNAL OF PETROLOGY DECEMBER 2000 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Fig. 5. (a) Hercynite, ilmenite and quartz inclusions in garnet. (b) Plagioclase + quartz symplectite corroding garnet in contact with biotite corona around matrix ilmenite. (c) Matrix assemblage of kyanite + biotite + garnet + plagioclase + K-feldspar + quartz all in textural equilibrium, from Zone F. (d) Granoblastic quartz matrix with sillimanite and garnet porphyroblasts, overprinted by fine sillimanite and biotite in S3 foliation seams (indicated). (Note fine sillimanite overgrowths on primary sillimanite.) All microphotographs in plane-polarized light and from Barun Gneisses [except (c)]. Field of view 1 mm across in (a) and (b) and 6·5 mm in (c) and (d). + green hornblende + sphene + haematite ± epidote ± quartz ± K-feldspar. Mafic gneisses have a wellaligned granoblastic matrix of quartz + andesine + hornblende + epidote at lowest structural levels (in Zones C–E) and bytownite + hornblende + biotite + garnet + meionite in a mafic calc-silicate from the highest structural level (Zone F) (Table 4; Fig. 2). Retrogressive blue–green hornblende coronas enclose clinopyroxene and hornblende. Epidote is enclosed by the clinopyroxene coronas, preserving a prograde metamorphic reaction. Retrogressive epidote occurs as coronas on garnet porphyroblasts. DESCRIPTION OF ROCKS ALONG THE KANGCHENJUNGA PROFILE Greater Himalaya Sequences Jannu–Kangchenjunga Gneiss The Jannu–Kangchenjunga Gneiss is composed of migmatized quartzo-feldspathic paragneisses, metapelitic garnet + cordierite gneisses and augen orthogneiss, with rare carbonate, meta-quartzite, calc-silicate and mafic gneiss. Metapelitic and quartzo-feldspathic gneisses have coarse, granoblastic textures and similar matrix assemblages of quartz + perthitic K-feldspar + plagioclase + orange biotite + sillimanite ± garnet ± cordierite (Table 5). Garnet porphyroblasts typically have few or no inclusions, these being ilmenite, biotite and quartz, and inclusion-free garnet overgrowths are recognized. Peak metamorphic sillimanite occurs both as porphyroblasts and as late-stage grains that overprint the main foliation. Sillimanite is a common coronal phase, occurring in aggregates of sillimanite + gedrite ± plagioclase, sillimanite ± biotite or sillimanite + phengite replacing cordierite and sillimanite ± gedrite ± biotite symplectites replacing biotite. As in the Barun Gneiss, porphyroblastic garnet is partially consumed by coronas and symplectites containing plagioclase, quartz, biotite and sillimanite (such as Fig. 5b) (Table 5). Fine-grained S3 seams of sillimanite ± gedrite are sub-parallel to and overprint the coarse-grained S1 foliation. Kyanite and staurolite are entirely absent from the Jannu– Kangchenjunga Gneiss. Muscovite is a rare matrix phase and typically replaces biotite. Mafic gneisses have a coarse, granoblastic matrix assemblage of green hornblende + quartz + plagioclase 1687 JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 Fig. 6. (a) Matrix staurolite + kyanite + quartz + biotite assemblage with late chlorite, from schist in Zone E. (b) Syn-tectonic garnet core with inclusion-free rims in biotite + muscovite + staurolite schist from Zone E. (c) Syn-tectonic garnet core with post-tectonic idiomorphic garnet overgrowth (indicated) in sericite + biotite schist from Zone B in the Seti Formation. (d) Syn-tectonic garnet overgrown by the biotite + sericite S1 metamorphic foliation (indicated), which in turn is overgrown by post-tectonic chlorite (indicated), also from Zone B. All microphotographs in plane-polarized light and field of view 6·5 mm across. + orange biotite. Calc-silicates have assemblages of plagioclase + meionite + clinopyroxene + forsterite + calcite + sphene and late-stage epidote (Table 6). The Jannu–Kangchenjunga Gneiss unit contains an intensely deformed biotite + sillimanite-bearing granitic augen orthogneiss body (Fig. 3). At high structural levels there are numerous, deformed leucogranite sills containing tourmaline + biotite ± muscovite ± garnet ± sillimanite. Lesser Himalayan Sequences Biotite Gneiss Fine-grained, biotite-rich schists dominate the wide ductile shear-zone between the Jannu–Kangchenjunga Gneiss unit and the highest-grade portion of the Himal Group (Zone G) to the south (Fig. 3). These schists have brown biotite + plagioclase + K-feldspar + quartz ± garnet assemblages with fine-grained sillimanite aligned within the foliation and as inclusions in feldspars and garnet (Table 5). Granitic pegmatites and veins are particularly common. Himal Group The Himal Group consists almost entirely of metapelite schists and schistosic gneisses, with uncommon quartzite, amphibolite and calc-silicate units of only 0·5–2 m thickness and few mappable (<100 m thick) calc-silicate units. This monotonous sequence of metapelitic schists permits the mapping of diagnostic metamorphic assemblages, and these display an increase in metamorphic grade to higher structural levels (inverted isograds). The following matrix assemblage zones are recognized, with increasing grade from south to north (Fig. 3; Table 5): Zone D: garnet+biotite+muscovite+oligoclase± kyanite; Zone E: garnet+biotite+muscovite+oligoclase+ staurolite±kyanite; Zone F: garnet+biotite+bytownite+K-feldspar+ kyanite; Zone G: garnet+biotite+bytownite+K-feldspar+ sillimanite±kyanite. The matrix assemblage consists of coarse granoblastic phases in textural equilibrium and is dominated by wellaligned (S1) red–brown biotite sheets (Figs 5c and 6a 1688 gn (pelite) (schist) (pelite) (pelite) (pelite) (ortho) K13 K67a K67b K67c K71a∗ K77 1689 (QFG) (pelite) (pelite) K65b K65c∗ K65d X X X X X X PHchl X (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) K6 K8a K8b K9 K62a K62b∗ (pelite) (pelite) K18a K18b K60a (pelite) (pelite) (pelite) K60b K60c∗ K60d Himal Group, Zone E (pelite) (pelite) K1 Himal Group, Zone D (pelite) K5 X X X X X X X X X X Kushma Formation, Zone D (pelite) K17 Ulleri Formation, Zone C (pelite) K12 Seti Formation, Zones A and B Sample X X-S C X-S X X-S X-S L Ikf S K X-S S S X S X-S X-S Ipl X S XR S SK S S S S mu present st X-S X-S present Ign X-S Ign X-S present X X-S X-S X X-S X-S SK X-S S S X Ign S X-S X-S X S S XS S X-L S bi X G kf X X X X X X X X X X X X X X X XG X pl sill Iky X X X X Ign X Ign X Ign X Ign X Ign X Ign X X X XG Ign X Ign X Ign X X X Ign X q X X Ign X present X present Ign X X present Ign X X X X ky Table 5: Petrology of aluminous rocks from the Kangchenjunga profile X Ign X-S X Ign X Ign X X Ign X X X Ign L Ign X Ign X Ign X-S Ign X Ign X-S X X Ign X SK Ign X-S X L Ign S ilm Eilm X hm X X Xtimg X mg Rbi SK Rgn Rgn,bi Rgn S-L Rgn X-S L R S-L Rgn S S-L Rgn S L Rst L Rgn chl cd gd zr, ap rut zr tm tm tm tm tm zr tm tm tm Accessory ky→mu mu laths across bi foliation ilm across chl across bi-mu foliation bi laths across mu foliation ilm and chl across foliation chl across Reaction textures GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA (pelite) (pelite) K61b∗ K61c (QFG) (pelite) (pelite) (pelite) K15cii K19a K21a K21b K58 (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) K25b K26a K26b K26c K27 K28b K29 1690 (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) (pelite) K35d K52 K53a K53b K53c K53d X X X X X X X X X X X X Ipl X X X X X X-S X-S XF X-S X-S S X X-S X-S S X-S X-S Ign X-S Ign X-S X-S X-S X X-S X-S X-S X-S X-S Rfeld Rbi Rbi,kf SR X-S L X-S st X X X X X X X X X X X X X X X X X X X X X Ipl kf X X X X X X X X X X X X X X X X X X X X X X X X pl Ign S-L C Ign S-L Cbi F feld S Ign, Ifeld X-S C F X-S C X-S X-S C X-S C S X X-S S-L C S-L F S-L F S-L F S-L F sill Ign X q X Ign X X X-S Ign X X X X X X X X X X X X X X X Ign X Ign X X X X Ign X Ign X X X X present Ign X ky X X X Ign X X X X Ign X Ign X X Ign X S Ign X Ign X-L X X Ign X-S Ign X ilm S-L hm X X X Xtimg mg Rbi Rgn Rbi Rbi R Rgn Rgn X-S Rbi X-S L chl cd S-L L gd py tm zr ap ap, tm zr, rut zr zr zr zr zr, tm ap, sph Accessory foliation sill inclusion foliation gn→sill, sill inc bi→sill foliation sill inclusion foliation sill inclusion sill→fine sill sill→fine sill sill→fine sill sill→fine sill bi→sill foliation, gn→chl sill across gn→cc-chl bi-chl mu laths across Reaction textures NUMBER 12 K35 X X mu Ign X-S Ign X-S bi VOLUME 41 Biotite Gneiss, MCT (pelite) K25a Himal Group, Zone G (pelite) (pelite) K15ci Himal Group, Zone F X (pelite) K60e X gn Sample Table 5: continued JOURNAL OF PETROLOGY DECEMBER 2000 (QFG) (QFG) (pelite) (pelite) (pelite) (orthogneiss) (QFG) (pelite) (pelite) (pelite) (QFG) K40b K40c K41a K41b K43 K45a K45b K46b K46d 1691 K47 K48 X-S X-S (schlieren) (schlieren) (leucogranite) K41dii K41di K41e X X X X Rfeld X Rbi, phl X Rbi, feld X Rfeld Rbi X-S S-L Cbi Ign X-S Rbi Ign X-S F X Rbi X C Ign XC X-L Rbi, feld st X X X X X X X X X X X X X X X X X X XRbi X kf X X X X X X X X X X XC XC X X X X X C X pl X S-L S-L F S-L F X S-L C F X S-L C F X S-L F S-L C F Ifeld S F Ifeld C S-L C Ifeld S-L F S-L S-L sill ky X X X X X Ign X X Ign X X X Ign X Ign X Ign X X Ign X X Ign X X IgnXC X q Rbi X X X Ign X Ign Ign X Ign X X X Ign X X ilm F hm Xtimg mg X PH X cd Rbi Rbi Rbi Rbi Rbi Rbi Rbi Rbi Rbi PH? Rbi X? Rbi,gn FRgn chl F F F? R R CF C gd ep tm, ap, phl ap, zr ap, zr ap, zr zr zr, ap zr tm tm, py rut rut zr tm zr Accessory bi→mu bi→sill bi→sill gn→sill-pl-bi bi→mu, gn→bi-pl, gn→gn over cd→sill-gd sill-phengite pseudomorphs bi→sill-gd-bi cd→gd-sill, cd→fine sill±bi cd→sill-gd±pl, gn→pl→q→sill-gd Reaction textures Samples grouped according to the mapped units (Fig. 3) and arranged in increasing metamorphic grade down page, that is, from south to north. I, inclusion phase with host indicated; X, peak metamorphic phase; S, aligned or contained within the dominant fabric in equilibrium with peak metamorphic phases; X-S, coarse aligned matrix; C, corona phase; K, crenulation phase; G, sub-grains; F, fine foliated seam; L, late phase; R, retrogressive phase; PH, pseudomorphed phase; ∗, analysed sample; present, recognized in hand specimen. Postscript indicates phase that is retrogressed or contains the inclusion. Reaction textures presented as reactant phases → product phases. (schlieren) K41c X-L C Rbi X-S Rbi X-L Cbi mu X-S Rbi, cd Ign X-S S-L C Rbi, pl Ign X-S Leucogranite in Kangchenjunga Gneiss X X X X X XC X X X X X-S (pelite) (pelite) (QFG) K37b K39 K40a X X X X (pelite) K37a bi X gn Kangchenjunga Gneiss K36 (pelite) Sample GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA q pl kf me (mafic) K71c X X X X (mafic) K65a X X X X X X (mafic) X X 1692 (mafic) X X (calc-silicate) K24 X X X X (mafic) (calc-silicate) K32a K32b (carbonate) K46c X X X X X X X X X X XC X X X X X C XC X R X-S X-S X-S X X-S X-S XS X hn X X X ol L X-L R R X-L X X-L X ep X X R R X cc act X X X XC C X X X X X sph X X X X X X X-S X-S bi/phl Rbi Rbi Rbi Rbi chl X X X X X X X X X X X X X X X op tm tm tm rut ap Accessory pl→me cpx→hn op→sph, cpx→hn hn→sph cpx→hn hn→gn Reaction textures Samples grouped according to the mapped units (Fig. 3) and arranged in increasing metamorphic grade down page; that is from south to north. I, inclusion phase; X, peak metamorphic phase; S, aligned or contained within the dominant fabric in equilibrium with peak metamorphic phases; X-S, coarse aligned matrix; C, corona phase; R, retrogressive phase; #, out of situ. Postscript indicates phase that is retrogressed or contains the inclusion. Reaction textures presented as reactant phases → product phases. (mafic) (mafic) K38 K46a X X X X X cpx NUMBER 12 Kangchenjunga Gneiss (mafic) K28a X C gn VOLUME 41 Himal Group, associated with Zone G metapelites (carbonate) K15a# Himal Group, associated with Zone F metapelites K61a Himal Group, associated with Zone E metapelites K63 Kushma Formation, associated with Zone D metapelites (mafic) (mafic) K11 K64 Ulleri Formation, associated with Zone C metapelites (calc-silicate) K71b Seti Formation, associated with Zone A-B metapelites Sample Table 6: Petrology of mafic; calc-silicate and carbonate rocks from the Kangchenjunga profile JOURNAL OF PETROLOGY DECEMBER 2000 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA and b). Garnet porphyroblasts typically have few or no inclusions, those observed being typically quartz, biotite and opaque minerals and, at highest grades (Zone G), aligned fine-grained sillimanite. Kyanite occurs as wellaligned porphyroblasts in textural equilibrium with the matrix assemblage. In Zones F and G there is latestage sillimanite aligned with S1 and thin, fine-grained sillimanite ± gedrite ± biotite ± haematite S3 seams. At highest grades (Zone G), sillimanite occurs as porphyroblasts in textural equilibrium with a matrix assemblage that is almost devoid of kyanite (Table 5). Partial melt segregations are recognized in metapelites in Zone G. In the high-grade zones, kyanite occurs both as laths in equilibrium with the matrix assemblage and as inclusions in garnet porphyroblasts. Primary sillimanite and biotite grains are overgrown by coronas of finegrained sillimanite (Table 5). Mafic gneisses also illustrate an increase in metamorphic grade from south to north, with matrix assemblages of quartz + plagioclase + green hornblende + epidote + sphene associated with Zones C–E metapelites and quartz + plagioclase + green hornblende + biotite ± clinopyroxene associated with Zone G metapelites (Table 6). At highest grades epidote and sphene are late-stage phases, with sphene forming coronas on opaque minerals and hornblende. Calc-silicates in Zones F and G have granoblastic matrix assemblages of quartz + plagioclase + clinopyroxene + calcite + sphene + orange biotite ± grossular with late-stage blue–green hornblende and epidote (Table 6). sphene ± biotite assemblages, rarely with small latestage epidote grains (Table 6). The basal Ulleri Formation is a mixed zone containing metapelitic schists with biotite + muscovite + quartz + K-feldspar ± plagioclase ± garnet assemblages and late-stage chlorite laths overprinting the main foliation (Table 5). Seti Formation The Seti Formation is a monotonous sequence of finegrained, pale grey–green quartz–sericite schists with rare mafic schists, calc-silicate and quartzite units. Biotite occurs in the quartz–sericite schists (Zone A); both within the fine-grained foliation and as aligned (L1) blocky laths that overgrow the sericite foliation. Garnet is common in the highest structural levels (Zone B) and occurs as large snowball poikiloblasts with post-tectonic idiomorphic overgrowths [Fig. 6c; see also Meier & Hiltner (1993)]. Chlorite is mostly a late-stage phase, such as coarse laths overprinting the sericite foliation (Fig. 6d) and retrogressive after garnet. Feldspars are nearly absent except for rare albite. Mafic schists have quartz + plagioclase + green hornblende ± epidote assemblages with garnet coronas on hornblende (Table 6). Deformed biotite + muscovite + tourmaline leucogranite bodies are emplaced in the core of the antiformal dome (Fig. 3). Kushma Formation The Kushma Formation is dominated by meta-quartzites and metapelitic schists with subordinate mafic schists. Metapelitic schists have a coarse granoblastic matrix with well-aligned micas and assemblages of garnet + quartz + albite + biotite + muscovite + tourmaline ± kyanite (Table 5), constituting Zone D assemblages. Garnet porphyroblasts are typically poikiloblastic and snowball garnets, commonly with inclusion-free rims. Late-stage crenulation cleavages with axial planar muscovite + biotite ± chlorite are common. All other chlorite is retrogressive after garnet and biotite. Mafic schists have quartz + plagioclase + blue–green hornblende + ilmenite + rutile assemblages. Ulleri Formation The Ulleri Formation is a 3–4 km thick unit of intensely sheared granitic orthogneiss (Fig. 3), with perthite augen in a granoblastic matrix of quartz + plagioclase + Kfeldspar + muscovite with minor biotite and sieve garnet. Mafic schists of 1 m width are at very low angles to the S1 foliation and interpreted as transposed mafic dykes. These have quartz + plagioclase + hornblende ± MINERAL CHEMISTRY Analyses from the Makalu profile were performed on the Cameca SX50 electron microprobe at the University of Tasmania, and from the Kangchenjunga profile on the Camebax microbeam at the University of Cape Town. For both, operating voltage was 15 kV and current was 20 nA for all phases except micas (10 nA) and feldspar (15 nA), and beam radius was 2·2 m for most phases and 4·4 m for micas and feldspars. The following natural silicates were used as standards and checked periodically throughout sessions. For garnet, staurolite and aluminosilicates: Kakanui pyrope, Kakanui hornblende, rutile and Chromite 52NL-11; for micas, amphiboles and pyroxenes: Kakanui hornblende, rutile and rhodonite; for feldspars: Lake County plagioclase, Nuni albite, Orthoclase-1 and Kakanui pyrope. The range in mineral chemistry, within each rock-type group in the Makalu profile, is summarized in Table 7. Representative mineral analyses are presented in Appendix C and the complete dataset of mineral analyses used is in the Journal of Petrology website at http://www.petrology.oupjournals.org. Mineral end-member activities were calculated largely after 1693 JOURNAL OF PETROLOGY VOLUME 41 Powell & Holland (1985); the formulations used are summarized in Appendix D. THERMOCALC RESULTS Average P–T conditions of equilibration of core and rim assemblages were determined by the method of Powell & Holland (1985, 1988). Calculations were performed using the 1992 thermodynamic dataset and computer program THERMOCALC v2·0b (Powell & Holland, 1988); the results are presented in Appendix A. All results satisfy the 2 test and errors from THERMOCALC, incorporating typical uncertainty for each mineral endmember activity and errors in the thermodynamic dataset, average ±44°C and ±1·0 kbar (Appendix A). Where there are insufficient independent reactions to calculate average P–T loci, the intersection of solutions for average P and average T calculations by THERMOCALC have been used to define equilibration conditions (Appendix A; Table 8). The best constrained P–T loci of both the maximum-T and maximum-P conditions preserved in each sample are summarized in Table 8 and represented in Fig. 7. Almost all samples (except M95a) were calculated by THERMOCALC. As a result of using this one method and this one internally consistent thermodynamic dataset, the resultant P–T calculations are considered directly comparable with each other. There were insufficient mineral end-members in the garnet + biotite quartzo-feldspathic gneiss assemblages, from the Black Gneisses, Biotite Gneiss and Jannu– Kangchenjunga Gneiss samples, to calculate a set of independent reactions and obtain average P–T results by THERMOCALC. Estimates of the equilibration conditions of these rock units, based on the stability field of matrix assemblages, are discussed below (Figs 7–9) and presented in Table 8. Average P–T results by THERMOCALC were, however, obtained from all other rock units. These results cluster into three distinct groups, corresponding to peak metamorphic conditions in the Barun Gneiss, Barun Gneiss samples that have re-equilibrated during cooling, and peak metamorphic conditions in the Himal Group samples (Fig. 7; Table 8). Peak metamorphic results from Barun Gneiss metapelites average 837 ± 59°C and 6·7 ± 1·0 kbar (Table 8), indicating an average geothermal gradient of 36 ± 4°C/ km. These peak results are identical to the phase stability field of Barun Gneiss assemblages (Fig. 8), and thus considered an accurate estimate of the peak of metamorphism (Fig. 7). THERMOCALC results from Barun Gneiss rims are of lower T (T = 24–56°C) and slightly lower P (P = 0·3–1·4 kbar) than core results (Fig. 7; Appendix A), although little significance can be attributed to core vs rim results because T and P are of similar magnitude NUMBER 12 DECEMBER 2000 to the errors (typically ±44°C and ±1·0 kbar) on these results (Appendix A). Barun Gneiss garnets preserve zoning of contrasting style to garnets from the LHS; with flat Ca patterns and increasing Fe and Mn and decreasing Mg from core to rims (Fig. 10). These patterns are typical of high-grade garnets with retrograde zoning caused by diffusion (Tracy, 1982; Tuccillo et al., 1990). Consequently, rim calculations may reflect only cation decoupling; where Ca remained relatively immobile and Fe–Mg continued to exchange during cooling (e.g. Tuccillo et al., 1990). As a result, the assemblages evolved along isopleths of the reaction 3 anorthite = grossular + 2 aluminosilicate+quartz, giving rise to the low, positive P/T, core to rim P–T array (Fig. 7). It is important to realize that this P–T array does not necessarily coincide with the P–T path experienced by the rocks (e.g. Frost & Chacko, 1989). Calc-silicate and mafic gneiss samples from the Barun Gneiss give results clustered around an average of 674 ± 33°C and 5·7 ± 1·1 kbar. These results are similar in pressure, but significantly lower in temperature (T = 75–125°C) than the phase stability fields of the matrix assemblages (Fig. 9). Consequently, these samples are interpreted to have been re-equilibrated at some stage subsequent to the peak of metamorphism. The calculated P–T loci represent either conditions at cessation of cation exchange (Harley, 1992) or re-equilibration in a later metamorphic event (see below). Re-equilibration of matrix mineral compositions in a later metamorphic event, without pervasive recrystallization of the matrix phases, has been documented in the Ungava Orogen (St-Onge & Ijewliw, 1996) and Zambezi Belt (Goscombe et al., 1998). Unlike the calcareous rock-types, metapelites and quartzo-feldspathic gneisses in the Barun Gneiss did not experience later re-equilibration and peak metamorphic mineral compositions were apparently preserved. Himal Group samples equilibrated at P–T conditions entirely distinct from all Barun Gneiss samples (Fig. 7). Calculations from syn-kinematic garnet cores vs postkinematic rims commonly document heating (T typically 4–36°C and one sample of 106°C) accompanying decompression (P up to 2·0 kbar) from near peak-P to peak-T conditions (Fig. 7). Himal Group samples from Zones D–F preserve peak metamorphic conditions averaging 609 ± 42°C and 8·8 ± 1·1 kbar, with the highest pressures recorded in the Kangchenjunga profile (Table 8). These results indicate an average geothermal gradient of 20 ± 2°C/km. The maximum-P conditions indicated by garnet core analyses, where significantly different from those for rims, average 557 ± 39°C and 10·0 ± 1·2 kbar in Zones D–F (Table 8). Heating accompanying decompression immediately before the peak T of metamorphism, as suggested by these P–T calculations, is typical of collisional orogens controlled 1694 0·03–0·06 0·00–0·01 X(Mn,M2) X(Fe 0·44–0·72∩ Fe2O3 wt % 0·22–2·21 1695 0·10 0·68–0·93 0·69–0·92∩ X(me) X(cc) Scapolite Carbonate 0·24–0·64 0·10 0·06–0·07 0·59–0·72 0·43–0·51 0·18–0·65∩ 0·30–0·45∩ 1·19–2·09∩ 0·26–0·44ℜ 0·98–0·99 0·80–0·84 0·23–0·54 0·32–0·44 0·43–0·53 0·14–0·98 0·26–0·43 0·40–0·53∩† 0·41–0·49 0·61–0·65ℜ 0·26–0·55ℜ 0·06–0·07ℜ 0·80–0·83∩ 0·21–0·39 0·36–0·43 0·34–0·51 0·23–1·51ℜ 0·27–0·44 0·41–0·49† 0·61 0·43–0·52 0·41–0·56 0·06–0·29 0·23–0·38 0·63–0·67 0·42–0·58 0·11–0·12 0·05–0·22 0·25–0·45 0·41–0·56∩ 0·06–0·29∩ 0·27–0·38ℜ 0·64–0·69ℜ 0·41–0·58∩ 0·01–0·22ℜ 0·21 0·11–0·23ℜ 0·85–0·88ℜ 0·00–0·02∩ 0·02–0·11∩ 0·11–0·15ℜ 0·03–0·13∩ 0·72–0·83ℜ rim 0·78 0·49–0·50 0·61–0·74 0·49 2·53–4·27 0·23–0·42 0·07–0·09 0·79ℜ 0·59–0·60ℜ 0·60–0·69∩ 0·49 1·12–4·32 0·27–0·43ℜ 0·07–0·10 0·48–0·96 0·33–0·35ℜ∗ 0·45–0·96 0·03–0·05 0·02–0·10 0·10–0·12 0·44–0·47∩ 0·36–0·44ℜ rim 0·26–0·34 0·03–0·05 0·02–0·10 0·10–0·11 0·43–0·48 0·37–0·43 core calc-silicate 0·51–0·61 0·08–0·10 0·33–0·46 0·66–0·69 0·54–0·67 0·08–0·13 0·06–0·13 core 0·56–0·60† 0·08–0·10 0·38–0·54ℜ 0·65–0·68∩ 0·54–0·64 0·15 0·06–0·13 rim biotite–muscovite schist Black Gneisses ∗In mafic schists. †Retrograde or late-stage phase; all other phases are primary matrix phases. a Ten per cent of FeOtotal calculated as Fe2O3; b100% of FeOtotal calculated as Fe2O3; c15% of FeOtotal calculated as Fe2O3; d60% of FeOtotal calculated as Fe2O3. Arrows indicate typical direction of change of component from core to rim. For micas; (A1,T1) signifies aluminium cations in T1 tetrahedal sites and (Na,A) signifies sodium cations in A sites (Powell & Holland, 1985). X(ep) Epidote Na+K 0·29–0·46 Al2O3 wt % Hornblendec XFe 0·24–0·44 1·86–2·29 XFe Cpx 0·41–0·53 0·61–0·64 0·25–0·54 XFe 0·05–0·09∩† 0·21–0·55∩† 0·67–0·69 0·46–0·58∩ 0·11–0·17∩ Chloritec 0·05–0·09 X(Na,A) 0·67–0·70 X(Al,T1) 0·27–0·58 0·48–0·58 XFe Biotitec XFe 0·10–0·19 X(ab) K-feldsparb Muscovited 0·12–0·37 Fe2O3 wt % 0·38–1·40 Fe2O3 wt % 0·84–0·88 0·00–0·02 0·00–0·13 0·04–0·16 0·03–0·16 0·69–0·84 core 0·15 0·56–0·86ℜ 0·00–0·03 0·02–0·04ℜ 0·12–0·29∩ 0·18–0·23∩ 0·53–0·65ℜ rim ZnO wt % 0·49–0·76 0·00–0·02 0·02–0·04 0·13–0·31 0·19–0·24 0·50–0·63 core XFe 0·18–0·38 0·27–0·43ℜ 0·00–0·02 0·06–0·12ℜ 0·09–0·18∩ 0·04–0·06∩ 0·70–0·80ℜ 0·26–0·42 ,M1) 0·14–0·27 X(Mg,M2) 3+ 0·03–0·09 XFe 0·60–0·78 X(Fe,M2) X(Ca,M2) rim core core rim aluminous schist mafic metapelite calc-silicate Himal Group Barun Gneiss Plagioclaseb X(an) Sillimaniteb Kyaniteb Staurolitea Cordierite Garnet Phase Table 7: Summary of mineral compositional ranges for rock-type groups in the Makalu profile GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 Table 8: Best estimates of metamorphic conditions of equilibration of samples from both profiles Sample Rock type Peak-T conditions P (kbar) Peak-P conditions T (°C) P (kbar) Method∗ T (°C) Lesser Himalayan Sequences Seti Formation, Zone B K71a gn–bi–mu–chl schist 7·4±1·9 561±11 8·3 assumed 490±32 8·5±1·9 561±12 (1) P–T, T & P 465±30 (1) T Himal Group and Ulleri Formation, Zone C M109 gn–ky schist K65c gn–bi–mu–chl schist 10·3±1·1 10·0 assumed 587±14 (1) P–T Himal Group and Kushma Formation, Zones D and E M106 gn–ky–st schist 8·4±1·0 650±20 9·7±1·1 646±21 (1) P–T M108 gn–ky–st schist 8·0±1·0 662±21 9·7±0·8 556±20 (1) P–T K61b gn–ky schist 9·7±1·1 596±13 K60c gn–ky schist 8·8±1·5 537±73 9·1±1·5 521±69 K62b gn–ky schist 9·6±1·1 596±14 8·9±1·1 608±28 9·5±1·1 574±37 average (1) P–T (1) P–T (1) P–T Himal Group, Zone F M1 calc-silicate 9·8±0·9 570±37 M2 gn–ky schist 8·0±1·2 612±95 9·9±1·9 600±95 (1) P–T (1) P & T, (2) P M4 calc-silicate 8·7±1·0 542±56 10·7±1·1 506±31 (1) P–T & T, (2) P M95a gn–ky gneiss 8·6±1·0 713±50 (2) P & T Biotite Gneiss, MCT in Kangchenjunga profile All samples sill–gn–kf–bi±gd >5·0–7·5 >670–750 (3) >4·0–7·5 >670–710 (3) Greater Himalayan Sequences Black Gneisses All samples sill–mu–kf–bi(±cd) Barun Gneiss and Namche Orthogneiss M39 gn–cd–sill metapelite 6·8±1·0 832±76 (1) P–T M51B gn–cd–sill metapelite 6·5±0·8 823±60 (1) P–T M58 gn–cd–sill QFG 7·0±0·6 877±55 (1) P–T M78 gn–sill metapelite 6·5±1·6 818±47 (1) P–T 6·7±1·0 837±59 Average All samples metapelite & mafic >5·5–7·0 >780–820 (3) Barun Gneiss samples that have been re-equilibrated M42 calc-silicate 5·7 assumed 695±18 (1) T M57 calc-silicate 5·7 assumed 700±30 (1) T M84 gn–hn–cpx–bi mafic 5·4±1·2 664±37 (1) P–T M87 gn–hn–cpx–bi mafic 6·0±1·1 635±47 (1) P–T 5·7±1·1 674±33 Average Jannu–Kangchenjunga Gneiss All samples QFG & calc-silicate >4·0–7·0 >710 (3) S3 assemblages, Barun Gneiss, Jannu–Kangchenjunga Gneiss and Himal Group Zone F All samples sill–gd–bi±gn >4·5–7·0 >610–750 (3) ∗Methods: (1) THERMOCALC results (Appendix A); average P–T calculations used preferentially; with average T or average P calculations used in few samples as indicated. (2) Conventional geothermobarometry (Appendix B). (3) Estimate based on phase stability field of matrix assemblage (see text). Error on average for rock group; is the average of result errors. Makalu profile samples indicated by prefix M and Kangchenjunga profile by prefix K. 1696 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Fig. 7. Summary of results of P–T calculations (Table 8) from both the Makalu and Kangchenjunga profiles. P–T estimates are largely calculated using THERMOCALC v2.0b (Powell & Holland, 1988) (Appendix A), with some by conventional geothermobarometry (Appendix B). Errors not included on figure for clarity; THERMOCALC errors average ±35°C and ±1·05 kbar in the Barun Gneisses and ±39°C and ±1·20 kbar in the Himal Group (Appendix A) and geothermobarometry errors are assumed to be of the order of ±50°C and ±1 kbar (Kohn & Spear, 1990). Closed symbols are cores, open are rims, small symbols are samples constrained by an assumed P and stars are average P–T loci (Table 8). Labelled metamorphic zones are as defined in the text and Fig. 8. Mineral equilibria from literature are: (1) Fe-staurolite stability field in FASH (Bickle & Archibald, 1984); (2 and 4) theoretical relations in KFMASH (Powell & Holland, 1990); (3) experimental reactions in tholeiite compositions with quartz–fayalite–magnetite (QFM) buffer (Spear, 1981); (5) experimental garnet-in reaction for mid-ocean ridge basalt of XFe = 40% composition (Green & Ringwood, 1967); (6) Grant (1985); (7) Bohlen & Dollase (1983); (9) theoretical relations in FMASH (Hensen & Green, 1973; Baker et al., 1987). by uplift and erosion (England & Thompson, 1984; Spear et al., 1984). As is typical at mid-amphibolite facies conditions, all Himal Group and Midland Group garnet porphyroblasts, except those closest to the MCT (i.e. sample M95a, which underwent retrograde diffusion), still preserve growth zoning from syn-kinematic cores to idiomorphic post-kinematic rims (Fig. 10). Garnet compositional profiles are typical of growth zoning with increasing Fe and Mg and decreasing Ca and Mn towards rims (Tracy, 1982; Loomis, 1983; Tuccillo et al., 1990). Thus core to rim P–T trajectories are considered at least broadly representative of changing P–T conditions during garnet growth, although T and P are of magnitude not much greater than the errors on these results (Appendix A). Core and rim calculations are interpreted to represent a portion of the prograde P–T path as argued by Spear et al. (1984) and St-Onge (1987). In this interpretation, garnet cores forming by syn-tectonic crystallization at the same time as the matrix assemblage, possibly near peak-P conditions, and idiomorphic and inclusion-free rims and overgrowths represent post-kinematic garnet growth at peak-T conditions. The continued growth of garnet on a heating and decompressional P–T path, as suggested by the textural relationships and thermobarometry, at first glance 1697 JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 Fig. 8. Petrogenesis of metapelitic and aluminous schist samples from both Makalu and Kangchenjunga profiles. Phase stability fields indicated are of matrix assemblages in the metamorphic zones discussed in the text. Arrows indicate transition to new mineral parageneses as indicated by mineral reaction textures (Tables 3 and 5). Stars are average P–T loci from the Barun Gneiss (Fig. 7; Table 8). Dashed lines indicate mineral equilibria in calcareous rocks for reference (Fig. 9). Dots are invariant points. Mineral equilibria from the literature are as in Fig. 7; (8) Chatterjee & Froese (1975); (10) Windom & Boettcher (1976) and Huckenholz et al. (1981); (11) Spear & Cheney (1989); (12) chlorite–muscovite stability in KMASH (Seifert, 1970; Bird & Fawcett, 1973); (13) K-feldspar solvus in NaKASH (Morse, 1970); (14) theoretical relations in NaKASH (Thompson, 1974). appears contrary to the general notion that increasing P favours garnet growth. This, however, can be rationalized by comparison with the calculations of Vance & Mahar (1998), made for similar bulk compositions to the samples of this study. The low T and high P of formation of LHS samples, and the typical heating and decompression paths documented (P/T = –0·057 to –0·023), would result in >5–10% increase in modal garnet, and thus garnet overgrowths. Anorthite component in matrix plagioclase increases towards rims (Table 7), indicating decompression coupled with garnet growth (e.g. St-Onge, 1987). Phase stability fields for garnet + kyanite + staurolite schist and epidote-bearing, garnet-free amphibolite assemblages typical of Zones D and E define conditions encompassing the THERMOCALC results from these zones (Figs 7–9). Within the errors of the THERMOCALC calculations, all results are consistent with the phase stability field of the matrix assemblage of the respective sample. Consequently, these average P–T results are considered plausible estimates of the equilibration conditions of the matrix assemblages in LHS samples. Further, peak-T results display a spread broadly consistent with an inverted metamorphic field gradient (Fig. 2; Tables 3 and 8). However, because of structural slicing (Fig. 2) and difficulty in assigning calcareous samples to the metapelite metamorphic zones (Table 8), in conjunction with the errors on P–T calculations, the 1698 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Fig. 9. Petrogenesis of mafic, calc-silicate and carbonate samples from both Makalu and Kangchenjunga profiles. Phase stability fields indicated are of calcareous matrix assemblages associated with the metapelitic metamorphic zones discussed in the text. Arrows indicate transition to new mineral parageneses as indicated by mineral reaction textures (Tables 4 and 6). Stars are average P–T loci from the Barun Gneiss (Fig. 7; Table 8). Dashed lines indicate mineral equilibria in aluminous rocks for reference (Fig. 8). Mineral equilibria from the literature are as in Figs 7 and 8 and: (15) Binns (1968) and Moody et al. (1983); (16) experimental relations with QFM buffer (Liou, 1973); (17) greenschist–amphibolite facies transition (Liou et al., 1974; Ghent et al., 1979); (18) Metz (1976) and Kase & Metz (1980); (19) Ellis (1978); (20) plagioclase-out in olivine tholeiite (Green & Ringwood, 1967); (21) Valley & Essene (1980); (22) relations in CaMS (Yoder, 1976). THERMOCALC results on their own are not sufficiently accurate to define an inverted metamorphic field gradient. GEOTHERMOBAROMETRY Conventional geothermometers and geobarometers have been used to constrain, in part, samples with insufficient mineral end-members to calculate average P–T loci by THERMOCALC (Appendix B). Calc-silicate sample (M4) from the Himal Group is constrained by the Fe-end-member garnet–hornblende–plagioclase geobarometer (Kohn & Spear, 1990) in conjunction with THERMOCALC average T results. Similarly, garnet + kyanite schist sample (M2) is constrained by the Mg-end-member garnet–biotite–muscovite–plagioclase geobarometer (Hoisch, 1991) and THERMOCALC average T results (Appendix A). A garnet + kyanite quartzofeldspathic gneiss sample (M95a) in Zone F is constrained by the garnet–plagioclase–kyanite geobarometer (Perchuk et al., 1985) and garnet–biotite geothermometer of Dachs (1990), with activities calculated by method of Hoinkes (1986) (Appendix B). Typical errors for these methods of ±50°C and ±1 kbar are assumed (Dachs, 1990; Kohn & Spear, 1990). These geothermobarometeric results are very similar to THERMOCALC results from other Himal Group samples and are 1699 JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 Fig. 10. Compositional profiles across garnets in metapelitic schists and gneisses. (a) GHS samples; M39 and M58 from Barun Gneiss and M78 from Namche Migmatite Orthogneiss. (b) LHS samples; M109, M108 and M95A from the Himal Group, Zones C, E and F, respectively, and K71a from the Seti Formation (Zone B). Horizontal scale is arbitrary. consistent with the phase stability field of the matrix assemblage (Fig. 8; Tables 3 and 4), and thus considered reliable. P–T EVOLUTION OF THE GREATER HIMALAYAN SEQUENCES Barun Gneiss and Namche Migmatite Orthogneiss Peak metamorphic conditions in the Barun Gneiss are tightly constrained by both the stability field of matrix assemblages and THERMOCALC results. Garnet + cordierite + sillimanite assemblages in metapelites with moderate XMg, such as the Barun Gneiss samples (Table 7), indicate an upper limiting P of 7·5 kbar (Powell & Holland, 1990). Garnet + clinopyroxene amphibolite assemblages indicate temperatures >800°C (Spear, 1981) and pressures >5·5 kbar (Green & Ringwood, 1967) (Figs 8 and 9). Partial melting is pervasive and muscovite is absent from metapelites and quartzo-feldspathic gneisses, suggesting temperatures >685°C (Chatterjee & Froese, 1975; Harte & Hudson, 1979). These constraints are consistent with peak metamorphic THERMOCALC results averaging 837 ± 59°C and 6·7 ± 1·0 kbar (Table 8). Sequential mineral parageneses in the Barun and Jannu–Kangchenjunga Gneisses document an anticlockwise P–T loop restricted to within the sillimanite field (Fig. 11a and b). Quartz + sillimanite + biotite + hercynitic spinel ± ilmenite inclusion assemblages (Table 3) indicate pre-peak metamorphic conditions of <4·0 kbar and >750°C (Fig. 8; Bohlen & Dollase, 1983; Grant, 1985; Vielzeuf & Montel, 1994). Prograde spinel was not stabilized at higher P by high f O2, because prograde and matrix Fe–Ti oxides are typically ilmenite ± rutile pairs (Table 3), suggesting reducing conditions (Powell & Sandiford, 1988). However, stabilization of spinel inclusions at P >4·0 kbar by ZnO cannot be discounted (Nichols et al., 1992). Because no other phase in the assemblage partitions ZnO, the topology of spinel phase relations is not significantly altered by increasing ZnO contents (Fig. 11a and b) (Nichols et al., 1992; Hand et al., 1994). Consequently, for any ZnO content, spinelbearing assemblages remain at lower pressures than the 1700 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA peak metamorphic garnet + sillimanite + cordierite fields (Hand et al., 1994; Scrimgeour & Hand, 1997) (Fig. 11a and b). Thus the overgrowing of spinel + sillimanite + biotite inclusion assemblages by garnet and ultimately peak-metamorphic cordierite + biotite + sillimanite + garnet assemblages is consistent with an up-P prograde path (Fig. 11b). Garnet overgrowths and coronas on garnet are common (Tables 3 and 5). At relatively high temperatures and low pressures, the increase in the proportion of garnet is consistent with increasing P or isobaric cooling (Vance & Mahar, 1998). Immediately post peak-metamorphic reaction textures indicate partial consumption of garnet by matrix cordierite, and subsequent partial replacement of cordierite by sillimanite + biotite (Tables 3 and 5). Such a sequence formed in response to either hydration or, as indicated by P–T pseudosections derived for similar composition metapelites (Scrimgeour & Hand, 1997), a low P/T decompressive cooling path or decompression followed by isobaric cooling (Fig. 11a). Both of these retrograde P–T vectors could potentially be preceded by either a clockwise or anticlockwise P–T loop and so do not constrain the progade P–T vector (Fig. 11a). Sillimanite inclusions in garnet are common, indicating that alumino-silicate was not totally consumed during prograde metamorphism, and because kyanite is entirely absent from inclusion and matrix assemblages, the prograde path was logically restricted to the sillimanite field (Fig. 11a). With no petrological or thermobarometry evidence for prograde metamorphism at high P, and early spinel parageneses and subsequent increases in the proportion of garnet both suggesting an up-P prograde path, an anticlockwise P–T loop is considered the most plausible for the GHS, particularly in conjunction with the isobaric retrograde path documented below (Fig. 11a and b). Peak metamorphic conditions were terminated by nearisobaric cooling. Matrix cordierite is replaced by sillimanite + gedrite-bearing assemblages and sillimanite + biotite assemblages (see above), indicating isobaric or up-P cooling across the following reactions with low, positive P/T (Hensen & Green, 1973; Spear & Rumble, 1986; Baker et al., 1987; Xu et al., 1994) (Figs 8 and 11a and b): cordierite + garnet → gedrite + sillimanite + quartz K-feldspar + cordierite + garnet → sillimanite + biotite. Numerous post-peak metamorphic reaction textures contain sillimanite (Tables 3 and 5). Kyanite is entirely absent, suggesting little or no increase in P during cooling. THERMOCALC results from matrix assemblages document a linear array from peak metamorphic conditions through re-equilibrated samples ranging to 600 ± 45°C Fig. 11. Phase relationships pertinent to Barun Gneiss and Jannu– Kangchenjunga Gneiss metapelites. (a) Qualitative P–T pseudosection in KFMASH with excess K-feldspar + quartz + melt, modified after Scrimgeour & Hand (1997). Solid arrows constrained by sequential mineral parageneses discussed in text; cordierite coronas on garnet replaced by sillimanite + biotite aggregates. Shaded P–T paths indicate alternative hypothetical prograde paths discussed in text: (i) constrained by spinel + sillimanite + biotite inclusions; (ii) highest pressure P–T path possible, constrained by common sillimanite and absence of kyanite inclusions. (b) Qualitative P–T pseudosection in KFMASHT + Zn modified after Hand et al. (1994) is considered to account for the possibility of Zn-bearing hercynite inclusions. Gedrite-forming reaction in FMASH is after Spear & Rumble (1986) and Xu et al. (1994). Solid arrow is P–T path constrained by the sequence of mineral parageneses: from sillimanite + spinel + biotite inclusion assemblages, overgrown by garnet and ultimately garnet + sillimanite + cordierite + biotite matrix assemblages, which are in turn overgrown by S3 gedrite + sillimanite assemblages. and 5·7 ± 1·1 kbar, defining a near-isobaric cooling path with P/T slope of 0·0035 kbar/°C (Fig. 7). S3 foliation seams with sillimanite ± gedrite ± biotite ± haematite assemblages are common in the Barun and Jannu–Kangchenjunga Gneiss units of the GHS, and are also present in the highest-grade Zones (F and G) of the underlying Himal Group. Garnet growth associated with S3 is indicated by overgrowths on garnet porphyroblasts, enveloping S3 with fine, aligned sillimanite and biotite 1701 JOURNAL OF PETROLOGY VOLUME 41 inclusions preserved. Kyanite is absent from S3 assemblages, indicating that the Barun Gneiss was not in the kyanite field during formation of the S3 foliation (Fig. 8). Peak, post-peak and S3 parageneses all formed in the sillimanite field, indicating that the retrograde trajectory of the Barun Gneiss did not pass through the kyanite field, further confirming a near-isobaric cooling path (Fig. 8). S3 assemblages indicate conditions of 610–750°C and 4·5–7 kbar (Hensen & Green, 1973; Spear & Rumble, 1986; Baker et al., 1987) (Fig. 8). These P–T conditions of S3 formation are compatible with THERMOCALC results from re-equilibrated mafic and calc-silicate rocks in the Barun Gneiss, which average 674 ± 33°C and 5·7 ± 1·1 kbar (Fig. 7; Table 8). S3 foliations are interpreted to have formed during the latest phase of ductile movement in the MCT, during the Neohimalayan metamorphic event. Re-equilibration of the mineral chemistry of matrix assemblages in calcsilicate and mafic rocks in the Barun Gneiss is interpreted to be due to this later metamorphic event. S3 assemblages are identical in both the GHS and upper LHS (Fig. 8), suggesting both terranes were juxtaposed and at the same crustal level during these latest ductile movements of the MCT. In the GHS, S3 conditions must have been preceded by a near-isobaric cooling trajectory (see above), and in the LHS by a near-isothermal decompression path from peak metamorphic conditions in the kyanite field (Fig. 8). Subsequent to S3, both metamorphic terranes broadly evolved together along the same decompression and cooling path. Numerous retrogressive chlorite-, epidote- and muscovite-forming reactions document continued cooling subsequent to formation of the S3 fabric (Tables 3–6; Figs 8 and 9). Retrogressive kyanite or andalusite has not been recognized and the exact barometric response during retrogression subsequent to S3 is unknown. Black Gneisses Coexisting K-feldspar and muscovite, sillimanite-bearing assemblages (Table 3) and the presence of partial melt segregations in the Black Gneisses are consistent with conditions centred on 670–710°C and 4–7·5 kbar (Chatterjee & Froese, 1975; Fig. 8). Cordierite is documented in the Black Gneisses (Pognante & Benna, 1993), limiting pressures to <6 kbar (Powell & Holland, 1990; Fig. 8). These are similar pressures, but lower temperatures, than experienced in the Barun Gneiss (Fig. 8). As in the underlying Barun Gneiss and Namche Migmatite Orthogneiss, kyanite is entirely absent from all Black Gneisses parageneses. Further, reaction textures containing secondary sillimanite (Table 3) and late-stage andalusite documented by Pognante & Benna (1993) indicate a retrogressive P–T path with low, negative NUMBER 12 DECEMBER 2000 P/T. Pervasive sillimanite needle inclusions within muscovite and feldspars suggest that conditions immediately before the peak of metamorphism were in the sillimanite field. Although the prograde P–T path in the Black Gneisses is largely unconstrained, these rocks are thought to have experienced an anticlockwise P–T path similar to the underlying Barun Gneiss, but nested within it at lower temperatures (Fig. 12). This is because, collectively, the entire section from Barun Gneiss to Black Gneisses is a contiguous metamorphic terrane, with temperatures decreasing upwards (i.e. not inverted isograds), that remained at moderate pressures throughout much of its metamorphic history. Partial melts generated in the hightemperature basal units of this terrane (i.e. in the Barun Gneiss) were emplaced as leucogranites within the Black Gneisses at higher structural levels (Patiño Douce & Harris, 1998; Harrison et al., 1999). Jannu–Kangchenjunga Gneiss Peak metamorphic conditions in the Jannu– Kangchenjunga Gneiss are tightly constrained by the stability field of matrix assemblages, these being similar to the Barun Gneiss. Garnet + cordierite + sillimanite metapelites indicate an upper limiting P of 7·5 kbar (Powell & Holland, 1990) (Figs 8 and 9). Sphene in mafic and calc-silicate gneisses indicate pressures >5 kbar (Spear, 1981). Partial melting is pervasive, muscovite is absent and K-feldspar is perthitic in metapelites and quartzo-feldspathic gneisses, suggesting temperatures >700°C (Chatterjee & Froese, 1975; Harte & Hudson, 1979). This is supported by meionite + plagioclase + clinopyroxene + forsterite assemblages in calc-silicate, indicating temperatures >750°C (Metz, 1976; Ellis, 1978; Kase & Metz, 1980). Like the Barun Gneiss, Jannu– Kangchenjunga Gneiss samples preserve sillimanite + gedrite + biotite reaction textures replacing cordierite, and all retrograde parageneses are restricted to the sillimanite field (Table 5). Consequently, the Jannu– Kangchenjunga Gneiss is interpreted to have experienced an isobaric cooling path similar to the Barun Gneiss. Spinel inclusions have not been recognized and the prograde P–T path cannot be constrained. P–T EVOLUTION OF THE LESSER HIMALAYAN SEQUENCES Himal and Midland Groups The LHS is a continuous Barrovian sequence that comprises the Himal Group and underlying Midland Group. The petrology and structural evolution of the LHS is very similar in both Makalu and Kangchenjunga profiles. 1702 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Matrix assemblages span a wide range of high-P amphibolite facies conditions typical of the Barrovian series and metamorphic isograds can be mapped (Figs 2 and 3). Metamorphic grade in both profiles increases towards higher structural levels and thus the sequence preserves an inverted metamorphic field gradient as is common in the LHS throughout the length of the Himalayas (LeFort, 1975; Hodges et al., 1988; Pecher, 1989; Inger & Harris, 1992; Vannay & Hodges, 1996). In the Kangchenjunga profile, the following metapelitic matrix assemblage zones are recognized, from south to higher structural levels in the north (Table 5): Zone A: albite+chlorite+muscovite±biotite (Seti Formation); Zone B: albite+chlorite+muscovite+biotite+garnet (Seti Formation); Zone C: albite+muscovite+biotite+garnet±chlorite (Ulleri Formation); Zone D: oligoclase+muscovite+biotite+garnet+ kyanite (Kushma Formation); Zone E: oligoclase+muscovite+biotite+garnet+ staurolite+kyanite; Zone F: bytownite+biotite+garnet+K-feldspar+ kyanite±muscovite; Zone G: bytownite+biotite+garnet+K-feldspar+ sillimanite+melt±kyanite. Only Zones C–F are recognized in the Makalu profile (Fig. 2). Estimates of peak-metamorphic conditions by THERMOCALC from both profiles also confirm a rise in T up-section (Fig. 7; Table 8). The average peakmetamorphic (peak-T) conditions, represented by rim results from garnet + kyanite ± staurolite assemblages in Zones D–F, are 609 ± 42°C and 8·8 ± 1·1 kbar (Table 8). Peak-metamorphic conditions by geothermobarometry from a garnet + kyanite + K-feldspar metapelitic gneiss in Zone F are 713 ± 50°C and 8·6 ± 1·0 kbar (Table 8), and are entirely consistent with the phase stability field of this sample. Zone G is represented only in the Kangchenjunga profile, with sillimanite ± kyanite + garnet + K-feldspar + melt + biotite and muscovite-absent assemblages, suggesting temperatures >700°C (Thompson, 1974; Chatterjee & Froese, 1975) and pressures (8 kbar at the kyanite– sillimanite transition (Fig. 8). The absence of perthitic K-feldspar sets a maximum limiting T of 770°C and absence of cordierite suggests P >6·5 kbar (Fig. 8). These conditions are supported by calc-silicate assemblages (Fig. 9). At lower structural levels, in the Midlands Group, an inverted metamorphic sequence from Zone A to Zone C is preserved, which is contiguous with and progresses up into the inverted metamorphic sequence in the overlying Himal Group (Figs 3 and 8). Matrix assemblages document an increase in metamorphic grade from >475°C above the biotite isograd in Zone A, across the garnet isograd into Zone B and to >550°C (Fig. 8) in Zone D. The only barometric constraints available for Zones A–C are the presence of kyanite and absence of sillimanite (Table 5), indicating pressures >5 kbar (Fig. 8). In contrast to the GHS, the Himal Group contains matrix kyanite and, in the highest-grade zones, kyanite inclusions, thus attesting to prograde metamorphism entirely within the kyanite field. The peak of metamorphism is interpreted to be post-tectonic (post-S1) and coeval throughout all zones in the LHS. In Zones A and B, snowball garnet growth is syn-tectonic (Powell & Vernon, 1979) and garnet is also deformed within the S1 foliation (Fig. 6d). Biotite porphyroblasts (Zone A) and idiomorphic garnet overgrowths in all zones (Fig. 6c) formed at the peak of metamorphism subsequent to S1. A portion of the prograde P–T path is preserved by garnet growth zoning, with core and rim P–T calculations indicating decompression accompanying heating (Fig. 7). Peak pressures experienced were at least 10 kbar as indicated by THERMOCALC results from mineral cores, implying considerable burial of the LHS and crustal thickening during the Himalayan Metamorphic Cycle. Peak-metamorphic conditions indicated by rim THERMOCALC results were at slightly lower pressures, with most samples between 7·4 and 8·8 kbar. Sequential mineral parageneses and reaction textures constrain portions of what are interpreted to be clockwise P–T paths in LHS zones (Fig 8 and 9). Decompression accompanying peak metamorphism is constrained to be wholly within the kyanite field in Zones A–E and to have passed into the sillimanite field only in the higher-grade Zones F and G (Fig. 8). In Zone G, matrix assemblages equilibrated near the kyanite–sillimanite transition (Table 5). In these samples, the near absence of matrix kyanite and presence of kyanite inclusions in garnet imply prograde metamorphism from the kyanite field into the sillimanite field at the peak of metamorphism (Fig. 8). Furthermore, rare aligned sillimanite inclusions in garnet rims are due to continued garnet growth, within the sillimanite field, at the peak of metamorphism. Garnet is also consumed by metamorphic reactions and even entirely enclosed by coronas of matrix plagioclase. Thus most garnet growth occurred early, during prograde metamorphism. This is also evident in the lower-grade zones, where garnet porphyroblasts are enclosed by the S1 foliation, and idiomorphic (peak-metamorphic) garnet overgrowths are relatively thin and rare (Tables 3 and 5). Late-stage S3 foliation seams with sillimanite ± gedrite ± biotite assemblages are developed in metapelites in Zones F and G (Table 5) and formed at conditions of 610–750°C and 4·5–7 kbar, identical to S3 in the GHS. The GHS and all LHS zones are interpreted to have evolved along a broadly similar retrogressive P–T path 1703 JOURNAL OF PETROLOGY VOLUME 41 subsequent to S3. Numerous retrograde biotite-, muscovite- and chlorite-forming reaction textures are recorded in all LHS zones and the GHS, documenting cooling to at least 400°C (Figs 8 and 9; Tables 3–6). No aluminosilicate phases were formed in post-S3 retrograde parageneses, thus the barometric response during cooling is unconstrained. NUMBER 12 DECEMBER 2000 km. The differing style of metamorphism of the upperplate GHS, with respect to the lower-plate LHS, is further supported by the contrasting interpretative P–T paths in these two metamorphic terranes (Fig. 8). Four tectonic models are discussed below, and their compatibility with observations of the structural and metamorphic evolution of Eastern Nepal is evaluated. Model A: inverted metamorphic sequence Biotite Gneiss (Main Central Thrust) Biotite Gneiss is a highly sheared unit of biotite-rich schistose gneisses representing the principal movement zone of the MCT, at the boundary between the Himal Group and the Jannu–Kangchenjunga Gneiss. Assemblages contain K-feldspar and are muscovite free, and partial melt segregations and granitic veins are common, indicating temperatures >670°C (Chatterjee & Froese, 1967; Thompson, 1974). K-feldspar is not perthitic, limiting temperatures to <750°C (Morse, 1970). Sillimanite is common, gedrite rare and kyanite and cordierite are entirely absent, restricting pressures to 5–7·5 kbar (Hensen & Green, 1973; Baker et al., 1987; Powell & Holland, 1990). These conditions are similar to those in the uppermost Himal Group (Zone G), except for the complete absence of kyanite, suggesting slightly lower P in the Biotite Gneiss unit. Sillimanite is found in inclusion, peak and retrograde parageneses, suggesting these rocks remained in the sillimanite field throughout shearing in the MCT. These assemblages and conditions of formation are very similar to S3 parageneses in the underlying Himal Group and throughout the GHS. Consequently, the S3 foliation is correlated with the movement phase in the MCT represented by Biotite Gneiss assemblages. S3 metamorphic conditions are correlated with the moderate-T/ moderate-P Neohimalayan metamorphic event (Hodges et al., 1988; Pecher, 1989; Vannay & Hodges, 1996). Although the MCT may have experienced multiple movement phases, complete recrystallization and a single pervasive foliation suggest that the Biotite Gneiss assemblages represent the main, and latest, phase of movement. The sense of movement during this period is unknown from the region investigated. DISCUSSION There is no petrological evidence for the GHS having ever been buried deeper than the 6·5 kbar experienced at the peak of metamorphism. Furthermore, peak metamorphism was on a high average geothermal gradient of 36 ± 4°C/km. This high-T/moderate-P metamorphism is incompatible with the high-P/moderate-T Barrovian metamorphism in the underlying LHS, with average geothermal gradients centred on >20 ± 2°C/ An inverted metamorphic sequence is documented from throughout the entire LHS, up to the base of the GHS. This feature is well recognized along the entire length of the Central Himalayan metamorphic front (e.g. LeFort, 1975; Hodges et al., 1988; Pecher, 1989; Inger & Harris, 1992; Vannay & Hodges, 1996). Conventional models of Himalayan metamorphism also include the GHS as the uppermost portion of this inverted metamorphic sequence. However, in the Eastern Nepal region investigated, metamorphism of the GHS was due to an entirely distinct average geothermal gradient and the GHS experienced a contrasting P–T path with respect to the underlying LHS. Thus the GHS cannot be considered continuous with, nor part of, the inverted metamorphic sequence displayed in the LHS. In the GHS, peak-T conditions also record the maximum depth of burial attained (Fig. 8). Peak metamorphic conditions in the GHS are interpreted to have been attained at the same time as peak-P conditions, before the peak of metamorphism, in the LHS (Fig. 12; Table 2). This is because matrix assemblages formed at the same time in both metamorphic terranes, in association with south-directed ductile overthrusting. These matrix assemblages crystallized at peak-P conditions (before peak-T) in the LHS and at peak-T conditions in the GHS. There are numerous models to explain the inverted metamorphic sequence in the Himalayas. These have been discussed in detail by LeFort (1975), Mohan et al. (1989), Searle & Rex (1989) and Hubbard (1996) and fall into two classes: (1) overthrusting of a hot slab, or (2) syn- to post-metamorphic overturning of isograds by ductile overthrusting (Mohan et al., 1989). These are not mutually exclusive and both processes are possibly involved in Himalayan metamorphism. Because metamorphic isograds are parallel to both the main metamorphic foliation and the MCT and associated thrusts, most models propose peak metamorphism synchronous with overthrusting in the MCT Zone. In the LeFort (1975) model, lower-crustal, high-grade metamorphic zones are thrust over the upper-crustal lower-grade zones at a rate faster than the relaxation of isotherms and the resultant inverted sequence preserved by rapid uplift. The anticlockwise P–T path in the GHS can potentially be rationalized as being the overthrust upper-plate in a modified version of this model (see Model D below). 1704 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA However, the distinctly higher average geothermal gradient of this metamorphic terrane needs to be explained. Model B: partially reworked pre-Himalayan parageneses Contrasting metamorphic styles of the GHS and LHS can potentially be resolved by a hypothetical scenario in which the GHS preserve pre-Himalayan metamorphic assemblages. For example, the high-T/moderate-P metamorphism of the GHS may have occurred during a Proterozoic metamorphic event in the Indian plate before collision and was juxtaposed against the LHS during the Himalayan Metamorphic Cycle. Such a scenario requires minimal reworking of the older assemblages during the Himalayan Metamorphic Cycle. Collisional belts incorporating basement terranes preserving older assemblages, which experienced only partial re-equilibration with minimal recrystallization of new assemblages, have been recognized (St-Onge & Ijewliw, 1996; Goscombe et al., 1998) and support the plausibility of such a model. Barnicoat & Treloar (1989) reported that relic preHimalayan tectonometamorphic and magmatic events are extensive in the GHS. However, there is no geochronological evidence in the Nepal Himalayas for Proterozoic metamorphic mineral assemblages. Nor is there conclusive structural evidence for reworked older fabrics in the GHS of the area investigated. Matrix mineral assemblages formed in association with the main S1–L1 fabric in both metamorphic terranes. This metamorphic fabric was coeval and tectonically equivalent in both terranes and associated with ductile overthrusting during the Himalayan Metamorphic Cycle. Similarly, S3 is precluded from being interpreted as the only expression of the Himalayan Metamorphic Cycle reworking older parageneses in the GHS, because identical and similarly late-formed S3 parageneses are also developed in the LHS. Further, in these hypothetical scenarios, the GHS would not necessarily be a hot slab during the Himalayan Metamorphic Cycle, rendering development of an inverted metamorphic sequence in the LHS unlikely. Model C: polyphase Himalayan Metamorphic Cycle A similar hypothetical scenario to Model B can be proposed in which the Himalayan Metamorphic Cycle involved two or more metamorphic events, which are variably preserved in different thrust-bound units in the orogen. For example, high-T/moderate-P parageneses in the GHS may have formed in an earlier event in the Himalayan Metamorphic Cycle and were insignificantly reworked by a later event in which the LHS was metamorphosed. Such a scenario is not considered a plausible explanation for the contrasting metamorphic styles recognized, for the same reasons as discussed for Model B. That is, the main metamorphic fabrics in both metamorphic terranes are considered essentially coeval and associated with the same south-directed ductile overthrusting. Two metamorphic events were recognized in the Himalayan Metamorphic Cycle by previous workers from throughout the Himalayas (Hodges et al., 1988; Pecher, 1989; Inger & Harris, 1992; Vannay & Hodges, 1996). The Eohimalayan event involved prograde Barrovian metamorphism in response to overthrusting and burial during collision in the Middle Eocene to Oligocene. The Early–Middle Miocene Neohimalayan event is characterized by moderate-T/moderate-P conditions in the sillimanite field. Together these metamorphic events comprise only portions (crystallization events) of the one metamorphic cycle (Passchier & Trouw, 1996), with a clockwise P–T path. High-P/moderate-T parageneses do not occur in the GHS. Thus metamorphism of the GHS cannot be simply resolved within the above scheme without modification (see Model D), because the prograde path was presumably of high-T/low-P type. Despite the inability of this model to account for the style of metamorphism in the GHS, both metamorphic events are recognized in the higher-grade zones of the LHS. Matrix mineral assemblages formed at peak-P conditions are of high-P/moderate-T type and consistent with the Eohimalayan event. Sillimanite-bearing S3 assemblages formed during decompression through the peak of metamorphism and constitute the Neohimalayan event. Similarly, two crystallization events are recognized in the GHS; peak metamorphic matrix assemblages associated with S1–L1 fabrics and a later event represented by S3 assemblages. Conditions of formation of S3 assemblages are the same in both the GHS and LHS, and correlated with the Neohimalayan event. Thus a twophase Himalayan Metamorphic Cycle is recognized in both metamorphic terranes, although currently accepted metamorphic conditions of both events are consistent with metamorphism of the LHS only. However, the prograde (Eohimalayan) phase is of contrasting metamorphic style in the two metamorphic terranes and metamorphism of the GHS in eastern Nepal is of a type unrecognized elsewhere in the GHS of the Himalaya. New geochronology suggests Barrovian metamorphism in the LHS may be as young as Late Miocene–Pliocene (Harrison et al., 1997), lending support to a polymetamorphic model. However, such a scenario requires an explanation for the absence of pervasive reworking of the overlying GHS in this younger event. It is plausible that the GHS slab may have been minimally structurally reworked in a Late Miocene–Pliocene event, with reworking being restricted to development of the S3 foliation and re-equilibration of calcareous rock-types without significant recrystallization. In such a scenario, the S3 foliation developed in the uppermost LHS must have 1705 JOURNAL OF PETROLOGY VOLUME 41 formed immediately after and be essentially coeval with the proposed Late Miocene–Pliocene main metamorphic foliation in the LHS. Model D: Paired Metamorphic Mountain Belt (preferred model) Most aspects of the juxtaposed, but contrasting, eastern Nepal metamorphic terranes can be resolved as a Paired Metamorphic Mountain Belt (PMMB) model analogous to the Acadian metamorphism in the Appalachians (Armstrong et al., 1992). The bases of this model are as follows: (1) different P–T paths are experienced at different crustal levels in a collisional orogen, particularly with respect to being upper-plate vs lower-plate components in a crustal-scale overthrusting system (Chamberlain & Karabinos, 1987; Shi & Wang, 1987; Armstrong et al., 1992). (2) Metamorphism is dynamic and intimately associated with progressive crustal-scale overthrusting, resulting in an inverted metamorphic sequence (Hubbard, 1996). (3) Different average geothermal gradients are generated by different tectonic processes within the respective metamorphic terranes. These are suggested to be (a) burial of an anomalously radiogenic GHS generating high-T/moderate-P metamorphism as modelled by Sandiford & Hand (1998) and previously suggested for the Himalayas by Pinet & Jaupart (1987), and (b) burial and crustal overthickening (England & Thompson, 1984) in the LHS. A possible model for the tectonometamorphic evolution of the eastern Nepal region investigated, combining these elements, is outlined below (Fig. 12). The GHS was buried to moderate crustal levels by low-angle overthrusting of the Tethyan Sequences from the north, along a proposed Eohimalayan Thrust (Vannay & Hodges, 1996). Prograde metamorphism of the GHS accompanied crustal thickening during collisional orogenesis, with peak-T conditions being reached at the maximum P experienced (Fig. 12a). The high average geothermal gradient of GHS metamorphism cannot be a conductive response to crustal thickening only (Sandiford & Powell, 1991). A possible explanation for the high-T/moderate-P metamorphism is an internally derived heat source resulting from a particularly radiogenic GHS (Sandiford & Hand, 1998). A high content of radiogenic elements in the GHS has previously been suggested as responsible for the high average geothermal gradient in this terrane (Pinet & Jaupart, 1987). The numerous leucogranite sills and plutons in the uppermost GHS are a consequence of the high average geothermal gradient and not the cause, as they are derived by partial melting from within the GHS (Patiño Douce & Harris, 1998). Furthermore, most leucogranite bodies are in the NUMBER 12 DECEMBER 2000 highest structural levels and lowest-grade portion of the GHS. Nelson et al. (1996) have reported deep seismic evidence for a zone of partial melting and granite melt pooling below the STDS, in the mid-crust, 100–200 km north of the currently exposed GHS. This indicates that high-grade metamorphism and partial melting in the GHS is a continuing process in which these hot rocks are being continually fed south like a conveyor belt as a result of overthrusting at the MCT, resulting in burial metamorphism and isograd inversion of the lower-plate LHS rocks. The GHS progressively overrode the lower-plate LHS, giving rise to high-P/moderate-T Barrovian metamorphism, with the burial phase of the LHS occurring at the same time as the peak of metamorphism in the GHS (Fig. 12b). Prograde metamorphism in both metamorphic terranes constitutes the Eohimalayan event and coincides with ductile, north over south thrusting and development of the pervasive main metamorphic fabric (S1–L1) in both terranes. This prograde (burial) phase culminated in peak-T conditions in the upperplate GHS (Fig. 12a) and peak-P conditions in the lowerplate LHS (Fig. 12b), with coarse-grained S1–L1 matrix assemblages crystallized at these respective conditions (Table 2). Progressive southward transport of the upperplate GHS accompanied cooling of this terrane while it remained at essentially the same relative crustal level, and concomitant progressive southward transport in the LHS resulted in inversion of the metamorphic isograds, this being further enhanced by the overriding hot GHS slab. Prograde metamorphism of the LHS was terminated by rapid decompression through the peak of metamorphism and in the highest-grade zones, progressed into the sillimanite field (Fig. 12c). It would be predicted by this PMMB model that peak metamorphic conditions in the GHS should slightly predate peak metamorphic (peak-T) conditions in the underlying LHS (Fig. 12b and c). This indeed seems to be the case, with peak metamorphism in the GHS being generally accepted as occurring at >22 Ma (Guillot et al., 1994; Harrison et al., 1995; Coleman & Hodges, 1998) and associated with emplacement of leucogranites in the High Himal with ages in the range 26–17 Ma (Scharer, 1984; Deniel et al., 1987; Hubbard & Harrison, 1989; Copeland et al., 1990; Inger & Harris, 1992; Metcalfe, 1993; Hodges et al., 1996; Edwards & Harrison, 1997; Coleman, 1998). Peak metamorphic conditions in the LHS are associated with overthrusting at the MCT and variously interpreted to be at 22–16 Ma (Hodges et al., 1996; Coleman, 1998; Coleman & Hodges, 1998), before cooling ages of 15–13 Ma (Macfarlane et al., 1992; Vannay & Hodges 1996) and possibly as late as 13–6 Ma (Macfarlane et al., 1992; Harrison et al., 1997). Postpeak, S3 foliations with identical assemblages in both 1706 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Fig. 12. Interpretative schematic representation of the preferred paired metamorphic belt model for the tectonometamorphic evolution of the Eastern Nepal region investigated. Age range of Eohimalayan and Neohimalayan metamorphic events and leucogranites as reported in the literature (see text). Dashed line is >600°C geotherm at the end of each indicated period. Dots follow the evolution of currently exposed and sampled metamorphic zones. Double arrows indicate pervasive ductile shearing and S1 foliation formation, and single arrow indicates movement on thrusts or extensional detachment faults. Highlighted portions of the interpretative P–T paths are for the period indicated; Β, average P–T loci at specific metamorphic periods in the different terranes (Fig. 7). 1707 JOURNAL OF PETROLOGY VOLUME 41 metamorphic terranes formed in association with shearing in the MCT (Figs 8 and 12c), implying that the terranes were juxtaposed and occupying approximately the same crustal level (>21 km depth; Fig. 8) at this time. These moderate-T/moderate-P parageneses constitute the Neohimalayan event (Pecher, 1989; Vannay & Hodges, 1996), of >22–13 Ma age (Hubbard & Harrison, 1989; Inger & Harris, 1992; Metcalfe, 1993; Table 2). Neohimalayan metamorphic conditions were immediately preceded by decompression from peak-metamorphic conditions in the LHS and by essentially isobaric cooling in the GHS (Fig. 8). Rapid decompression and extrusion of the entire metamorphic core of the Himalayas, contemporaneous with contraction (Fig. 12c and d), has been proposed by Burg et al. (1984), England & Molnar (1993), Hodges et al. (1993, 1996) and Vannay & Hodges (1996). Extrusion of deep-seated rocks is accommodated by extensional detachments, during episodic periods of orogenic collapse in response to lithospheric overthickening (Hodges et al., 1996). The contrasting P–T trajectories, immediately before S3, in the upper- and lower-plate terranes in the area investigated could plausibly be explained by reactivation of the MCT as an extensional detachment (Fig. 12c). However, extensional movements during the Miocene were accommodated by the STDS above the GHS (Burchfiel et al., 1992; Hodges et al., 1996; Vannay & Hodges, 1996; Edwards & Harrison, 1997), and extensional reactivation of the MCT, at the base of the GHS, has not been documented. Although a PMMB model may explain much of the distribution of palaeo-metamorphic conditions in the investigated portion of eastern Nepal, it remains a working hypothesis. In particular, the sequence of thrusting and extensional movements on major structures in the region needs to be established and the proposed mechanism for the high average geothermal gradient in the GHS has not been tested. Furthermore, comparably high average geothermal gradients and anticlockwise P–T paths, such as proposed for the portion of the GHS investigated, have not been reported from elsewhere in the GHS. Consequently, the model discussed here may not be generally applicable throughout the entire Himalayan metamorphic front. Studies of deeply eroded, older metamorphic belts elsewhere have shown them to be heterogeneous, being composed of metamorphic terranes of contrasting metamorphic style [i.e. Appalachians (Armstrong et al., 1992), Zambezi Belt (Goscombe et al., 1998), Damara Orogen (Miller, 1983), Ubendian Belt (Ring, 1993), Grenville Province (Rivers et al., 1989) and Ungava Orogen (St-Onge & Ijewliw, 1996)], suggesting that it may be unlikely to expect comparable metamorphic profiles throughout the entire 1600 km length of the Central Himalaya metamorphic front. NUMBER 12 DECEMBER 2000 ACKNOWLEDGEMENTS Thanks are due to Mel Lambourne for the inspiration, and to Ned Stolz, Wolf Zwolsman, Lia, Prem, Kancha Lama, Ram Lama, Rhim Bardo Lama, Pemba Dike Sherpani and Shiring Nema Sherpani for their great company and help amongst those bonnie hills. Constructive comments by Drs K. Hodges, T. Hoisch, C. M. Wilson, S. Harley and S. Sorensen, and discussions with Drs R. Oliver and M. Sandiford, greatly improved this paper, and their efforts were very much appreciated. Dr Wieslaw Jablonski at the University of Tasmania and Dick Rickards at the University of Cape Town are thanked for their help with electron microprobe analyses. REFERENCES Armstrong, T. R., Tracy, R. J. & Hames, W. E. (1992). Contrasting styles of Taconian, Eastern Acadian and Western Acadian metamorphism, central and western New England. Journal of Metamorphic Geology 10, 415–426. Baker, J., Powell, R., Sandiford, M. & Muhling, J. (1987). Corona textures between kyanite, garnet and gedrite in gneisses from Errabiddy, Western Australia. Journal of Metamorphic Geology 5, 357–370. Barnicoat, A. C. & Treloar, P. J. (1989). Himalayan metamorphism—an introduction. Journal of Metamorphic Geology 7, 3–8. Berman, R. G. (1990). Mixing properties of Ca–Mg–Fe–Mn garnets. American Mineralogist 75, 328–344. Bhanot, V. B., Sigh, V. P., Kansal, A. K. & Thakur, V. C. (1977). Early Proterozoic Rb–Sr whole rock age for Central Crystalline Gneiss of higher Himalaya, Kumaun. Geological Society of India Journal 18, 90–91. Bickle, M. J. & Archibald, N. J. (1984). Chloritoid and staurolite stability: implications for metamorphism in the Archaean Yilgarn Block, Western Australia. Journal of Metamorphic Geology 2, 179–203. Binns, R. A. (1968). Hydrothermal investigations of the amphibolite– granulite facies boundary. Geological Society of Australia Special Publication 2, 341–344. Bird, G. W. & Fawcett, J. J. (1973). Stability relations of Mg-chlorite– muscovite and quartz between 5 and 10 kbar water pressure. Journal of Petrology 14, 415–428. Bohlen, S. R. & Dollase, W. A. (1983). Calibration of hercynite–quartz stability. Geological Society of American Bulletin 15, 529. Bordet, P. (1961). Recherches géologiques dans l’Himalaya du Nepal, région du Makalu. Paris: CNRS. Brunel, M. & Kienast, J.-R. (1986). Étude petro-structurale des chevauchements ductiles himalayens sur la transversale de l’Everest– Makalu (Nepal oriental). Canadian Journal of Earth Science 23, 1117– 1137. Burchfiel, B. C., Zhiliang, C., Hodges, K. V., Yuping, L., Royden, L. H., Changrong, D. & Jiene, X. (1992). The South Tibetan Detachment System, Himalayan Orogen: extension contemporaneous with and parallel to shortening in a collisional mountain belt. Geological Society of America, Special Paper 269, 41 pp. Burg, J. P., Guiraud, M., Chen, G. M. & Li, G. C. (1984). Himalayan metamorphism and deformations in the North Himalayan Belt (southern Tibet, China). Earth and Planetary Science Letters 69, 391–400. Chamberlain, C. P. & Karabinos, P. (1987). Influence of deformation on pressure–temperature paths of metamorphism. Geology 15, 42–44. 1708 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Chatterjee, N. D. & Froese, E. (1975). A thermodynamic study of the pseudobinary join muscovite–paragonite in the system KAlSi3O8– NaAlSi3O8–Al2O3–SiO2–H2O. American Mineralogist 60, 985–993. Coleman, M. E. (1998). U–Pb constraints on Oligocene–Miocene deformation and anatexis within the central Himalaya, Marsyandi Valley, Nepal. American Journal of Science 298, 553–571. Coleman, M. E. & Hodges, K. V. (1998). Contrasting Oligocene and Miocene thermal histories from the hanging wall and footwall of the South Tibetan detachment in the central Himalaya from 40Ar/ 39 Ar thermochronology, Marsyandi Valley, central Nepal. Tectonics 17(5), 726–740. Copeland, P., Harrison, T. M. & LeFort, P. (1990). Age and cooling of the Manaslu granite: implications for Himalayan tectonics. Journal of Volcanology and Geothermal Research 44, 33–50. Dachs, E. (1990). Geothermobarometry in metasediments of the southern Grossvenediger area (Tauern Window, Austria). Journal of Metamorphic Geology 8, 217–230. Deniel, C., Vidal, P., Fernandez, A., Le Fort, P. & Peucat, J.-J. (1987). Isotopic study of the Manaslu granite (Himalaya, Nepal): inferences on the age and source of Himalayan leucogranites. Contributions to Mineralogy and Petrology 96, 78–92. Edwards, M. A. & Harrison, T. M. (1997). When did the roof collapse? Late Miocene north–south extension in the high Himalaya revealed by Th–Pb monazite dating of the Khula Kangri granite. Geology 25, 543–546. Ellis, D. E. (1978). Stability and phase equilibria of chloride and carbonate bearing scapolites at 750°C and 4000 bar. Geology 42, 1271–1281. England, P. & Molnar, P. (1993). Cause and effect among thrust and normal faulting, anatectic melting and exhumation in the Himalaya. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 401–411. England, P. C. & Thompson, A. B. (1984). Pressure–temperature–time paths of regional metamorphism 1. Heat transfer during the evolution of regions of thickened continental crust. Journal of Petrology 25, 894–928. Frost, B. R. & Chacko, T. (1989). The granulite uncertainty principle: limitations on thermobarometry in granulites. Journal of Geology 97, 435–450. Ghent, E. D., Robbins, D. B. & Stout, M. Z. (1979). Geothermometry, geobarometry, and fluid compositions of metamorphosed calc-silicates and pelites, Mica Creek, British Columbia. American Mineralogist 64, 874–885. Goscombe, B. D., Armstrong, R. & Barton, J. M. (1998). Tectonometamorphic evolution of the Chewore Inliers: partial re-equilibration of high-grade basement during the Pan-African Orogeny. Journal of Petrology 39, 1347–1384. Grant, J. A. (1985). Phase equilibria in partial melting of pelitic rocks. In: Ashworth, J. R. (ed.) Migmatites. Glasgow: Blackie, pp. 86–144. Green, D. H. & Ringwood, A. E. (1967). An experimental investigation of the gabbro to eclogite transformation and its petrological applications. Geochimica et Cosmochimica Acta 31, 767–833. Guillot, S., Hodges, K. V., LeFort, P. & Pecher, A. (1994). New constraints on the age of the Manaslu leucogranite: evidence for episodic tectonic denudation in the central Himalayas. Geology 22, 559–562. Hand, M., Scrimgeour, I., Powell, R., Stuwe, K. & Wilson, C. J. L. (1994). Metapelitic granulites from Jetty Peninsula, east Antarctica: formation during a single event or by polymetamorphism? Journal of Metamorphic Geology 12, 557–573. Harley, S. L. (1992). Proterozoic granulite terranes. In: Condie, K. C. (ed.) Proterozoic Crustal Evolution. Developments in Precambrian Geology 10. Amsterdam: Elsevier, pp. 301–359. Harris, N. B. W., Inger, S. & Massey, J. A. (1993). The role of fluids in the formation of High Himalayan leucogranites. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 323–340. Harrison, T. M., McKeegan, K. D. & LeFort, P. (1995). Detection of inherited monazite in the Manaslu leucogranite by 208Pb/232Th ion microprobe dating: crystallization age and tectonic implications. Earth and Planetary Science Letters 133, 271–282. Harrison, T. M., Ryerson, F. J., LeFort, P., Yin, A., Lovera, O. M. & Catlos, E. J. (1997). A Late Miocene–Pliocene origin for the Central Himalayan inverted metamorphism. Earth and Planetary Science Letters 146, E1–E7. Harrison, T. M., Grove, M., McKeegan, K. D., Coath, C. D., Lovera, O. M. & LeFort, P. (1999). Origin and episodic emplacement of the Manaslu Intrusive Complex, Central Himalaya. Journal of Petrology 40, 3–19. Harte, B. & Hudson, N. F. C. (1979). Pelite facies series and the temperatures and pressures of Dalradian metamorphism in E Scotland. In: The Caledonides of the British Isles—Reviewed. London: Geological Society. Hensen, B. J. & Green, D. H. (1973). Experimental study of the stability of cordierite and garnet in pelitic compositions at high pressures and temperatures. 3. Synthesis of experimental data and geological applications. Contributions to Mineralogy and Petrology 38, 151–166. Hodges, K. V., Hubbard, M. S. & Silverberg, D. S. (1988). Metamorphic constraints on the thermal evolution of the central Himalayan Orogen. Philosophical Transactions of the Royal Society of London, Series A 326, 257–280. Hodges, K. V., Burchfiel, B. C., Royden, L. H., Chen, Z. & Liu, Y. (1993). The metamorphic signature of contemporaneous extension and shortening in the central Himalayan orogen: data from the Nyalam transect, southern Tibet. Journal of Metamorphic Geology 11, 721–737. Hodges, K. V., Hames, W. E., Olszewski, W., Burchfiel, B. C., Royden, L. H. & Chen, Z. (1994). Thermobarometric and 40Ar/ 39 Ar geochronologic constraints on Eohimalayan metamorphism in the Dinggye area, southern Tibet. Contributions to Mineralogy and Petrology 117, 151–163. Hodges, K. V., Parrish, R. R. & Searle, M. P. (1996). Tectonic evolution of the central Annapurna Range, Nepalese Himalayas. Tectonics 15, 1264–1291. Hodges, K., Bowring, S., Davidek, K., Hawkins, D. & Krol, M. (1998). Evidence for rapid displacement on Himalayan normal faults and the importance of tectonic denudation in the evolution of mountain ranges. Geology 26(6), 483–486. Hoinkes, G. (1986). Effect of grossular-content in garnet on the partitioning of Fe and Mg between garnet and biotite. Contributions to Mineralogy and Petrology 92, 393–399. Hoisch, T. D. (1990). Empirical calibration of six geobarometers for the mineral assemblage quartz + muscovite + biotite + plagioclase + garnet. Contributions to Mineralogy and Petrology 104, 225–234. Hoisch, T. D. (1991). Equilibria within the mineral assemblage quartz + muscovite + biotite + garnet + plagioclase, and implications for the mixing properties of octahedrally-coordinated cations in muscovite and biotite. Contributions to Mineralogy and Petrology 108, 43–54. Hubbard, M. S. (1989). Thermobarometric constraints on the thermal history of the Main Central Thrust Zone and Tibetan Slab, eastern Nepal Himalaya. Journal of Metamorphic Geology 7, 19–30. Hubbard, M. S. (1996). Ductile shear as a cause of inverted metamorphism: example from the Nepal Himalaya. Journal of Geology 104, 493–499. 1709 JOURNAL OF PETROLOGY VOLUME 41 Hubbard, M. S. & Harrison, T. M. (1989). 40Ar/39Ar age constraints on deformation and metamorphism in the Main Central Thrust Zone and Tibetan Slab, eastern Nepal Himalaya. Tectonics 8, 865–880. Huckenholz, A. G., Lindhuber, W. & Fehr, K. T. (1981). Stability relations of grossular + quartz + wollastonite + anorthite I. The effect of andradite and albite. Neues Jahrbuch für Mineralogie, Abhandlungen 142, 223–247. Inger, S. & Harris, N. B. W. (1992). Tectonothermal evolution of the High Himalayan Crystalline Sequence, Langtang Valley, northern Nepal. Journal of Metamorphic Geology 10, 439–452. Kase, H. R. & Metz, P. (1980). Experimental investigation of the metamorphism of siliceous dolomites. Contributions to Mineralogy and Petrology 73, 151–159. Kohn, M. J. & Spear, F. S. (1990). Two new geobarometers for garnet amphibolites, with applications to southeastern Vermont. American Mineralogist 75, 89–96. LeFort, P. (1975). Himalayas: the collided range. Present knowledge of the continental arc. American Journal of Science 275, 1–44. LeFort, P., Pecher, A. & Upreti, B. N. (1986). A section through the Tibetan Slab in central Nepal (Kali Gandaki Valley): mineral chemistry and thermobarometry of the Main Central Thrust Zone. Sciences de la Terre 47, 211–228. Liou, J. G. (1973). Synthesis and stability relations of epidote, Ca2Al2FeSi3O12(OH). Journal of Petrology 14, 381–413. Liou, J. G., Kuniyoshi, S. & Ito, K. (1974). Experimental studies of the phase relations between greenschist and amphibolite in a basaltic system. American Journal of Science 274, 613–632. Lombardo, B., Pertusati, P & Borghi, S. (1993). Geology and tectonomagmatic evolution of the eastern Himalaya along the Chomolunga–Makalu transect. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 323–340. Loomis, T. P. (1983). Compositional zoning of crystals: a record of growth and reaction history. In: Saxena, S. K. (ed.) Kinetics and Equilibrium in Mineral Reactions. New York: Springer, pp. 1–60. MacFarlane, A. M., Hodges, K. V. & Lux, D. (1992). A structural analysis of the Main Central thrust zone, Langtang National Park, central Nepal Himalaya. Geological Society of America Bulletin 104, 1389–1402. Maruo, Y. & Kizaki, K. (1981). Structure and metamorphism in Eastern Himalaya. In: Saklani, P. S. (ed.) Metamorphic Tectonites of the Himalaya. Delhi: Today and Tomorrows Printers and Publishers, pp. 175–230. Meier, K. & Hiltner, E. (1993). Deformation and metamorphism within the Main Central Thrust zone, Arun Tectonic Window, eastern Nepal. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 511–523. Metcalfe, R. P. (1993). Pressure, temperature and time constraints and metamorphism across the Main Central Thrust zone and High Himalayan Slab in the Garhwal Himalaya. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 485–509. Metz, P. (1976). Experimental investigation of the metamorphism of siliceous dolomites III. Equilibrium data for the reaction: 1 tremolite + 11 dolomite = 8 forsterite + 13 calcite + 9 CO2 + 1 H2O for the total pressures of 3,000 and 5,000 bars. Contributions to Mineralogy and Petrology 58, 137–148. Miller, C. & Frank, W. (1992). Geochemistry and isotope geology of Proterozoic and early Palaeozoic granitoids in the NW Himalayas. 7th Himalay–Karakorum–Tibet Workshop, Oxford University. Miller, R. Mc. G. (1983). The Pan-African Damara Orogen of Namibia. In: Miller, R. Mc. G. (ed.) The Damara Orogen. Geological Society of South Africa, Special Publication 11, 431–515. NUMBER 12 DECEMBER 2000 Miyashiro, A. (1973). Metamorphism and Metamorphic Belts. London: George Allen and Unwin. Mohan, A., Windley, B. F. & Searle, M. P. (1989). Geothermobarometry and the development of inverted metamorphism in the Darjeeling– Sikkim region of the eastern Himalaya. Journal of Metamorphic Geology 7, 95–110. Moody, J. B., Meyer, D. & Jenkins, J. E. (1983). Experimental characterisation of the greenschist/amphibolite boundary in mafic systems. American Journal of Science 283, 48–92. Morrison, C. W. K. & Oliver, G. J. H. (1993). A study of illite crystallinity and fluid inclusions in the Kathmandu Klippe and the Main Central Thrust zone, Nepal. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 525–540. Morse, S. A. (1970). Alkali feldspars with water at 5 kbar pressure. Journal of Petrology 11, 221–251. Nelson, K. D., Zhao, W., Brown, L. D., Kuo, J., Che, J., Liu, X., Klemperer, S. L., Makovsy, Y., Meissner, R., Mechie, J., Kind, R., Wenzel, F., Ni, J., Nabelek, J., Leshou, C., Tan, H., Wei, W., Jones, A. G., Booker, J., Unsworth, M., Kidd, W. S. F., Hauck, M., Alsdorf, D., Ross, A., Cogan, M., Wu, C., Sandvol, E. & Edwards, M. (1996). Partially molten middle crust beneath southern Tibet: synthesis of project INDEPTH results. Science 274, 1684–1688. Nichols, G. T., Berry, R. F. & Green, D. H. (1992). Internally consistent gahnitic spinel–cordierite–garnet equilibria in the FMASHZn system: geothermobarometry and applications. Contributions to Mineralogy and Petrology 111, 362–377. Parrish, R. R. & Hodges, K. V. (1996). Isotopic constraints on the age and provenance of the Lesser and Greater Himalayan sequences, Nepalese Himalaya. Geological Society of America Bulletin 108, 904–911. Passchier, C. W. & Trouw, R. A. J. (1996). Microtectonics. Berlin: Springer, pp. 1–3. Patiño Douce, A. E. & Harris, N. (1998). Experimental constraints on Himalayan anatexis. Journal of Petrology 39, 689–710. Pecher, A. (1989). The metamorphism in the Central Himalaya. Journal of Metamorphic Geology 7, 31–41. Perchuk, L. L., Aranovich, L. Ya., Podlesskii, K. K., Lavrent’eva, I. V., Gerasimov, V. Yu., Fed’kin, V. V., Kitsul, V. I., Karsakov, L. P. & Berdnikov, N. V. (1985). Precambrian granulites of the Aldan Shield, eastern Siberia, USSR. Journal of Metamorphic Geology 3, 265–310. Pinet, C. & Jaupart, C. (1987). A thermal model for the distribution in space and time of the Himalayan granites. Earth and Planetary Science Letters 84, 87–99. Pognante, U. & Benna, P. (1993). Metamorphic zonation, migmatization and leucogranites along the Everest transect of Eastern Nepal and Tibet: record of an exhumation history. In: Treloar, P. J. & Searle, M. P. (eds) Himalayan Tectonics. Geological Society, London, Special Publications 74, 323–340. Powell, C. McA. & Vernon, R. H. (1979). Growth and rotation history of garnet porphyroblasts with inclusion spirals in a Karakoram schist. Tectonophysics 54, 25–43. Powell, R. & Holland, T. J. B. (1985). An internally consistent thermodynamic dataset with uncertainties and correlations: 1. Methods and a worked example. Journal of Metamorphic Geology 3, 327–342. Powell, R. & Holland, T. J. B. (1988). An internally consistent dataset with uncertainties and correlations: 3. Applications to geobarometry, worked examples and a computer program. Journal of Metamorphic Geology 6, 173–204. Powell, R. & Holland, T. (1990). Calculated mineral equilibria in the pelite system, KFMASH (K2O–FeO–MgO–Al2O3–SiO2–H2O). American Mineralogist 75, 367–380. 1710 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Powell, R. & Sandiford, M. (1988). Sapphirine and spinel phase relationships in the system FeO–MgO–Al2O3–SiO2–TiO2–O2 in the presence of quartz and hypersthene. Contributions to Mineralogy and Petrology 98, 64–71. Ring, U. (1993). Aspects of the kinematic history and mechanisms of superposition of Proterozoic mobile belts of eastern Central Africa (northern Malawi and southern Tanzania). Precambrian Research 62(3), 207–226. Rivers, T., Martignole, J., Gower, C. F. & Davidson, A. (1989). New tectonic divisions of the Grenville Province, southeast Canadian Shield. Tectonics 8, 63–84. Sandiford, M. & Hand, M. (1998). Australian Proterozoic high-temperature, low-pressure meatmorphism in the conductive limit. In: Treloar, P. J. & O’Brien, P. J. (eds) What Drives Metamorphism and Metamorphic Reactions? Geological Society, London, Special Publications 138, 109–120. Sandiford, M. & Powell, R. (1991). Some remarks on high-temperature– low-pressure metamorphism in convergent orogens. Journal of Metamorphic Geology 9, 333–340. Scharer, U. (1984). The effect of initial 230Th disequilibrium on young U–Pb ages: the Makalu case, Himalaya. Earth and Planetary Science Letters 67, 191–204. Scrimgeour, I. & Hand, M. (1997). A metamorphic perspective on the Pan African overprint in the Amery area of MacRobertson Land, East Antarctica. Antarctic Science 9, 313–335. Searle, M. P. & Rex, A. J. (1989). Thermal model for the Zanskar Himalaya. Journal of Metamorphic Geology 7, 127–134. Seifert, F. (1970). Low-temperature compatibility relations of cordierite in haplopelites of the system K2O–MgO–Al2O3–SiO2–H2O. Journal of Petrology 11, 73–99. Sengor, A. M. C. (1990). Plate tectonics and orogenic research after 25 years: a Tethyan perspective. Earth-Science Reviews 27, 1–201. Shi, Y. & Wang, C.-Y. (1987). Two-dimensional modeling of the P–T–t paths of regional metamorphism in simple overthrust terrains. Geology 15, 1048–1051. Shrestha, S. B., Shrestha, J. N. & Sharma, S. R. (1984). Geological Map of Eastern Nepal, 1:250 000. Lainchour, Kathmandu: Ministry of Industry, Department of Mines and Geology. Shrestha, S. B., Shrestha, J. N. & Sharma, S. R. (1985). Geological Map of Central Nepal, 1:250 000. Lainchour, Kathmandu: Ministry of Industry, Department of Mines and Geology. Spear, F. S. (1981). Amphibole–plagioclase equilibria: an empirical model for the relation albite + tremolite = edenite + 4 quartz. Contributions to Mineralogy and Petrology 77, 355–364. Spear, F. S. & Cheney, J. T. (1989). A petrogenetic grid for pelitic schists in the system SiO2–Al2O3–FeO–MgO–K2O–H2O. Contributions to Mineralogy and Petrology 101, 149–164. Spear, F. S. & Rumble, D., III (1986). Pressure, temperature and structural evolution of the Orfordville Belt, west central New Hampshire. Journal of Petrology 27, 1071–1093. Spear, F. S., Selverstone, J., Hickmott, D., Crowley, P. & Hodges, K. V. (1984) P–T paths from garnet zoning: a new technique for deciphering tectonic processes in crystalline terranes. Geology 12, 87–90. St-Onge, M. R. (1987). Zoned poikiloblastic garnets: P–T paths and syn-metamorphic uplift through 30 km of structural depth, Wopmay Orogen, Canada. Journal of Petrology 28, 1–21. St-Onge, M. R. & Ijewliw, O. J. (1996). Mineral corona formation during high-P retrogression of granulitic rocks, Ungava Orogen, Canada. Journal of Petrology 37, 553–582. Thompson, A. B. (1974). Calculation of muscovite–paragonite–alkali feldspar phase relations. Contributions to Mineralogy and Petrology 44, 173–194. Tracy, R. J. (1982). Compositional zoning and inclusions in metamorphic minerals. In: Ferry, J. M. (ed.) Characterization of Metamorphism through Mineral Equilibria. Mineralogical Society of America, Reviews in Mineralogy 10, 355–397. Tuccillo, M. E., Essene, E. J. & van der Pluijm, B. A. (1990). Growth and retrograde zoning in garnets from high-grade metapelites: implications for pressure–temperature paths. Geology 18, 839–842. Valley, J. W. & Essene, E. J. (1980). Akermanite in the Cascade Slide Xenolith and its significance for regional metamorphism in the Adirondacks. Contributions to Mineralogy and Petrology 74, 143–152. Vance, D. & Mahar, E. (1998). Pressure–temperature paths from P–T pseudosections and zoned garnets: potential, limitations and examples from the Zanskar Himalaya, NW India. Contributions to Mineralogy and Petrology 132, 225–245. Vannay, J.-C. & Hodges, K. V. (1996). Tectonometamorphic evolution of the Himalayan metamorphic core between the Annapurna and Dhaulagiri, central Nepal. Journal of Metamorphic Geology 14, 635–656. Vielzeuf, D. & Montel, J.-M. (1994). Partial melting of metagreywackes. Part 1. Fluid-absent experiments and phase relationships. Contributions to Mineralogy and Petrology 117, 375–393. Windom, K. E. & Boettcher, A. L. (1976). The effect of reduced activity of anorthite on the reaction grossular + quartz = anorthite + wollastonite; a model for plagioclase in the Earth’s lower crust and upper mantle. American Mineralogist 61, 889–896. Xu, G., Will, T. M. & Powell, R. (1994). A calculated petrogenetic grid for the system K2O–FeO–MgO–Al2O3–SiO2–H2O, with particular reference to contact-metamorphosed pelites. Journal of Metamorphic Geology 12, 99–119. Yoder, H. S., Jr (1976). Relationship of melilite-bearing rocks to kimberlite; a preliminary report on the system akermanite–CO2. First International Conference on Kimberlites. Physics and Chemistry of the Earth 9, 883–894. 1711 T (°C) P (kbar) 766±57 M78-rim 789±56 M51B-rim 5·7±0·7 6·5±0·8 6·5±1·1 6·8±1·0 6·5±0·9 7·0±0·6 5·1±2·0 6·5±1·7 f 1·08 1·14 1·66 1·56 1·37 0·88 1·58 1·32 1712 635±47 600±45 M87-core M87-rim 561±11 587±14 584±13 596±14 592±13 596±13 521±69 537±73 K71a–overgrowth K65c-core K65c-rim K62b-rim K61b-core K61b-rim K60c-core K60c-rim 8·8±1·5 9·1±1·5 9·7±1·1 9·3±1·1 9·6±1·1 10·1±1·1 10·3±1·1 7·4±1·9 8·0±1·0 9·7±0·8 10·7±1·1 0·804 0·795 0·068 0·082 0·104 0·103 0·177 −0·105 0·80 0·36 0·84 0·73 0·47 0·86 0·39 0·71 0·73 1·18 0·89 0·90 0·86 0·83 0·93 0·89 0·97 0·55 0·85 0·087 0·192 0·040 0·018 0·556 0·654 0·652 0·561 0·542 0·270 0·290 alm–py–gr–phl–ann–east–naph–mu–an–ab–q–ky–H2O alm–py–gr–phl–ann–east–naph–mu–an–ab–q–ky–H2O alm–py–gr–phl–ann–east–naph–mu–an–ab–clin–daph–ames–q–ky–H2O alm–py–gr–phl–ann–east–naph–mu–an–ab–clin–daph–ames–q–ky–H2O alm–py–gr–phl–ann–east–naph–mu–an–ab–clin–daph–ames–q–H2O alm–py–gr–phl–ann–east–naph–mu–an–ab–clin–daph–ames–q–H2O alm–py–gr–phl–ann–east–naph–mu–an–ab–clin–daph–q–H2O alm–py–phl–ann–east–naph–mu–cel–pa–clin–daph–ames–q–H2O alm–gr–phl–ann–mu–cel–mst–fst–an–ky–ames–q–H2O alm–gr–east–mu–cel–fst–an–ky–ames–q–H2O alm–gr–phl–ann–mu–cel–mst–fst–an–ky–ames–q–H2O alm–gr–phl–ann–mu–mst–fst–an–ky–ames–q–H2O alm–an–ab–di–hed–cats–hb–parg–ed–me–q–H2O–CO2 alm–py–gr–an–ab–parg–ed–di–hed–cats–q–H2O alm–py–gr–an–ab–parg–ed–di–hed–cats–q–H2O alm–py–gr–an–ab–parg–ed–anth–di–hed–cats–q–H2O alm–py–gr–an–ab–parg–ed–anth–di–hed–cats–q–H2O alm–py–gr–phl–ann–an–ab–hb–parg–ed–di–cats–q–H2O alm–py–gr–phl–ann–an–ab–hb–parg–di–cats–q–H2O alm–py–crd–fcrd–phl–ann–naph–an–ksp–ab–sill–q–H2O alm–py–crd–fcrd–phl–ann–east–naph–ksp–ab–sill–q–H2O alm–py–gr–crd–fcrd–phl–ann–an–ksp–sill–q–H2O alm–py–crd–fcrd–phl–ann–ksp–sill–q–H2O alm–py–gr–crd–fcrd–sill–ksp–an–ab–phl–ann–naph–east–q–H2O alm–py–crd–fcrd–sill–ksp–an–ab–phl–ann–naph–q–H2O alm–py–gr–sill–phl–ann–mu–cel–an–ksp–q–H2O alm–py–gr–sill–phl–ann–mu–cel–an–ksp–q–H2O Mineral end-members used NUMBER 12 Himal Group, Kangenjunga profile 662±21 M108-rim 8·4±1·0 650±20 556±20 646±21 M106-core M106-rim 506±31 M4-core∗ M108-core 9·7±1·1 570±37 9·8±0·9 553±36 M1-rim 9·4±0·9 5·7±1·1 6·0±1·1 5·9±1·3 5·4±1·2 0·853 0·869 0·869 0·833 0·882 0·848 0·920 0·950 cor VOLUME 41 M1-core Himal Group, Makalu profile 645±43 664±37 M84-rim M84-core Barun Gneiss samples that have been re-equilibrated, Makalu profile 806±84 823±60 832±76 M39-core M51B-core 821±72 M58-rim M39-rim 877±55 M58-core Barun Gneiss, Makalu profile 818±47 M78-core Namche Migmatite Orthogneiss, Makalu profile Sample Table A: Pressure–temperature conditions, derived using THERMOCALC v2.0b (Powell & Holland, 1988) APPENDIX A JOURNAL OF PETROLOGY DECEMBER 2000 400°C 5·8±1·02 Average T results (°C) Sample f 7·0 kbar 8·3±1·87 8·0±1·84 6·8±1·17 500°C 1713 683±31 568±34 M57-core∗ M57-rim∗ 643±95 637±94 471±30 492±30 M2-core M2-rim M109-core M109-rim 561±12 K71a-over 1·1 1·1 0·4 0·5 0·6 0·5 1·7 1·6 1·5 1·1 562±12 543±11 492±32 470±31 620±95 625±95 537±56 614±32 726±30 720±18 0·6 0·7 0·5 0·4 0·3 0·2 1·7 1·4 1·4 1·1 f 0·5 0·8 0·4 f 560±12 540±12 489±32 467±32 604±95 608±95 543±60 645±33 760±29 752±17 9·0 kbar 8·6±2·07 8·3±2·03 7·9±1·32 600°C 0·4 0·7 0·7 0·2 0·2 0·3 1·8 1·4 1·3 1·0 f 0·3 0·6 0·1 f alm–phl–ann–east–naph–mu–cel–pa–daph–ames–q–H2O alm–phl–ann–east–naph–mu–cel–pa–daph–ames–q–H2O alm–py–phl–ann–east–ky–ames–q–H2O alm–py–phl–ann–east–ky–ames–q–H2O alm–py–phl–ann–naph–mu–pa–ky–q–H2O alm–py–phl–ann–naph–mu–pa–ky–q–H2O alm–py–an–hb–di–hed–cats–me–q–H2O hb–ftr–anth–di–hed–cats–ann–naph–an–ab–ksp–cc–q–H2O–CO2 hb–ftr–anth–di–hed–ann–naph–an–ab–ksp–cc–q–H2O–CO2 hb–ftr–anth–di–hed–phl–ann–naph–an–ab–ksp–cc–q–H2O–CO2 Mineral end–members used alm–phl–ann–east–naph–mu–cel–pa–daph–ames–q–H2O alm–phl–ann–east–naph–mu–cel–pa–daph–ames–q–H2O alm–gr–ann–an–mu–ky–q–H2O Mineral end-members used f, 2 test; cor, correlation coefficient. Mineral end-member abbreviations after Powell & Holland (1988). ∗Calculated for XCO2 = 0·5 and XH2O = 0·5. 542±12 K71a-core Himal Group, Kangchenjunga profile 522±54 M4-rim∗ Himal Group, Makalu profile 680±18 M42-core∗ Barun Gneiss samples that have been re-equilibrated, Makalu profile 5·0 kbar 7·9±1·67 0·9 7·7±1·76 K71a–over 1·1 0·7 K71a-core Himal Group, Kangchenjunga profile M2-rim f Average P results (kbar) Himal Group, Makalu profile Sample GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Rock-type 1714 calc-silicate gn–ky schist gn-ky gneiss gn–ky gneiss gn–ky gneiss gn–ky gneiss M4-rim M2-core M95a-core M95a-rim M95a-core M95a-rim Himal Group, Makalu profile Sample 641±50 701±50 650±50 710±50 2·7±1 2·3±1 4·2±1 4·6±1 8·8±1 8·3±1 500°C 6·1±1 6·5±1 9·9±1 9·3±1 600°C gn–pl–ky (Perchuk et al., 1995) gn–bi (Dach, 1990; activity, Hoinkes, 1986) gn–bi–mu–pl–q (Hoisch, 1991; Mg) gn–hn–pl (Kohn & Spear, 1990; Fe#B) Calibration VOLUME 41 632±50 692±50 7·7±1 7·3±1 400°C 8·0 kbar 4·0 kbar 6·0 kbar P result (kbar) T result (°C) Table B: Geothermobarometry results from the Makalu profile APPENDIX B JOURNAL OF PETROLOGY NUMBER 12 DECEMBER 2000 1715 0·72 X(Fe) 0·00 0·00 8·02 OH Sum 0·00 0·00 0·00 Na K 0·87 8·02 0·00 0·10 0·37 0·17 0·81 0·10 Mn 2·40 0·00 0·09 2·06 Fe2+ Mg 0·01 Fe3+ 2·00 0·00 0·00 2·98 99·46 0·00 0·00 0·03 1·12 3·06 2·45 35·33 0·03 20·83 Ca 0·00 1·98 Cr Al 2·97 98·68 Total 0·00 0·00 H 2O Si 0·02 K 2O Ti 0·02 1·07 CaO Na2O 1·42 6·83 MnO MgO 0·25 30·79 21·04 Al2O3 Fe2O3 0·05 Cr2O3 FeO 0·01 0·00 0·00 36·60 37·19 SiO2 rim core TiO2 M58 M58 Sample: gn gn Mineral: 0·81 8·01 0·00 0·00 0·00 0·73 0·42 0·07 1·84 0·04 1·95 0·00 0·00 2·96 99·16 0·00 0·00 0·02 8·51 3·55 1·07 27·48 0·70 20·71 0·05 0·05 37·02 core M87 gn 0·84 8·04 0·00 0·00 0·00 0·68 0·36 0·12 1·91 0·05 1·95 0·00 0·00 2·97 99·39 0·00 0·01 0·00 7·86 2·99 1·75 28·41 0·82 20·58 0·01 0·00 36·95 rim M87 gn 0·88 8·04 0·00 0·00 0·01 0·19 0·34 0·01 2·52 0·01 1·99 0·00 0·00 2·97 99·23 0·00 0·01 0·03 2·12 2·76 0·17 36·85 0·23 20·63 0·00 0·04 36·39 core M106 gn 0·89 8·03 0·00 0·00 0·00 0·19 0·32 0·05 2·50 0·00 2·01 0·00 0·00 2·96 99·08 0·00 0·00 0·03 2·17 2·60 0·72 36·52 0·00 20·79 0·03 0·05 36·17 rim M106 gn 0·78 8·05 0·00 0·00 0·00 1·29 0·33 0·29 1·18 0·09 1·91 0·00 0·00 2·96 101·21 0·00 0·00 0·00 15·47 2·84 4·34 18·09 1·49 20·83 0·05 0·07 38·02 core M1 gn Table C: Representative matrix mineral analyses from the Makalu profile APPENDIX C 0·79 8·05 0·00 0·00 0·00 1·31 0·32 0·30 1·17 0·08 1·92 0·00 0·00 2·95 101·10 0·00 0·00 0·01 15·71 2·72 4·50 17·93 1·43 20·87 0·00 0·01 37·93 rim M1 gn 0·84 33·43 4·00 0·01 0·01 0·00 0·62 0·02 3·37 0·04 17·60 0·00 0·11 7·65 98·04 2·11 0·02 0·02 0·00 1·47 0·10 14·17 0·20 52·52 0·00 0·50 26·92 core M106 st st 0·85 33·38 4·00 0·01 0·00 0·00 0·60 0·02 3·36 0·04 17·48 0·01 0·12 7·74 98·30 2·12 0·02 0·01 0·00 1·43 0·10 14·18 0·20 52·34 0·03 0·56 27·32 rim M106 cd 0·40 11·03 0·00 0·00 0·03 0·00 1·19 0·02 0·79 0·03 3·96 0·00 0·00 5·01 96·36 0·00 0·01 0·15 0·01 7·54 0·17 8·91 0·40 31·78 0·00 0·00 47·40 core M58 cd 0·37 11·02 0·00 0·00 0·03 0·00 1·24 0·02 0·72 0·04 3·97 0·00 0·00 5·00 96·56 0·00 0·01 0·14 0·01 7·94 0·22 8·15 0·55 32·00 0·00 0·00 47·54 rim M58 sill — 3·00 0·00 0·00 0·00 0·00 0·00 0·00 0·00 0·01 2·00 0·00 0·00 0·99 98·34 0·00 0·01 0·00 0·02 0·01 0·00 0·00 0·38 61·74 0·01 0·01 36·16 core M58 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA 3·46 3·67 19·12 Al2O3 8·94 0·00 0·02 8·31 0·00 0·19 9·86 3·99 100·33 MnO MgO CaO Na2O K 2O H 2O Total 3·10 1716 0·19 1·05 0·00 Fe2+ Mn 0·53 X(Fe) 2·00 2·00 9·71 OH Sum 0·92 0·03 0·95 Na K 1·01 0·50 9·68 0·02 0·00 0·93 0·00 Mg Ca 2·70 0·50 9·69 2·00 0·80 0·07 0·00 1·08 0·00 1·08 0·19 1·68 0·09 0·49 9·73 2·00 0·84 0·07 0·00 1·10 0·00 1·07 0·19 1·66 0·09 2·71 99·05 3·96 8·73 0·44 0·03 9·72 0·00 16·92 3·32 18·62 1·58 35·74 rim M106 bi 0·67 9·66 2·00 0·91 0·03 0·00 0·62 0·01 1·26 0·22 1·82 0·14 2·65 100·51 3·95 9·40 0·23 0·00 5·49 0·14 19·86 3·90 20·30 2·41 34·83 core M71 bi 0·64 9·68 2·00 0·91 0·02 0·00 0·70 0·01 1·26 0·22 1·76 0·13 2·67 100·54 3·95 9·43 0·16 0·00 6·14 0·22 19·77 3·88 19·60 2·25 35·15 rim M71 bi 0·29 8·97 2·00 0·69 0·28 0·00 0·05 0·00 0·02 0·03 2·81 0·02 3·07 100·12 4·57 8·25 2·21 0·00 0·48 0·00 0·39 0·65 36·36 0·43 46·77 core M106 mu 0·30 8·98 2·00 0·69 0·26 0·00 0·07 0·00 0·03 0·05 2·81 0·02 3·05 99·25 4·52 8·19 2·04 0·00 0·68 0·00 0·60 1·00 35·88 0·31 46·03 rim M106 mu 0·36 8·96 2·00 0·86 0·08 0·00 0·07 0·00 0·04 0·06 2·73 0·03 3·09 100·63 4·55 10·21 0·63 0·03 0·71 0·00 0·74 1·23 35·09 0·52 46·91 core M71 mu 0·45 8·96 2·00 0·86 0·09 0·00 0·06 0·00 0·05 0·07 2·74 0·02 3·07 100·87 4·55 10·25 0·72 0·00 0·62 0·00 0·88 1·46 35·31 0·44 46·63 rim M71 mu 0·47 17·90 8·00 0·00 0·00 0·01 2·22 0·00 1·98 0·35 2·79 0·01 2·54 100·43 11·54 0 0 0·05 14·34 0·04 22·76 4·46 22·77 0·07 24·39 core M106 chl 0·45 17·26 2·00 0·00 0·18 2·70 2·14 0·02 1·74 0·31 1·02 0·18 6·97 99·16 1·99 0·01 0·60 16·72 9·52 0·17 13·80 2·71 5·75 1·60 46·20 core M87 hn 0·51 17·32 2·00 0·14 0·18 2·13 1·99 0·02 2·09 0·37 1·49 0·03 6·88 99·38 1·98 0·75 0·62 13·15 8·85 0·14 16·53 3·24 8·33 0·28 45·52 core M42 hn 0·48 17·79 2·00 0·43 0·31 1·97 2·05 0·09 1·93 0·34 2·50 0·10 6·07 99·22 1·96 2·19 1·04 12·03 8·97 0·68 15·05 2·95 13·86 0·83 39·67 core M1 hn NUMBER 12 0·00 1·01 0·18 1·68 1·70 Al Fe3+ 2·66 0·20 2·65 0·21 Si 99·11 3·97 8·25 0·51 0·05 9·56 0·08 17·02 3·34 18·90 1·63 35·80 core M106 bi VOLUME 41 Ti 98·70 3·94 9·53 0·14 0·00 15·82 3·27 16·69 Fe2O3 FeO 18·74 35·03 35·19 SiO2 rim core TiO2 M58 M58 Sample: bi bi Mineral: Table C: continued JOURNAL OF PETROLOGY DECEMBER 2000 1717 0·00 0·00 0·00 0·20 0·76 0·01 0·00 0·00 0·00 4·98 Fe2+ Mn Mg Ca Na K Sr Ba OH Sum — 0·00 Fe3+ X( Fe) 1·20 Al 99·37 Total 0·00 0·00 H2O Ti 0·00 BaO 2·81 0·00 SrO Si 0·10 0·00 MgO K2O 11·15 0·00 MnO 4·19 0·00 FeO 8·83 0·05 Fe2O3 CaO 0·00 22·86 Al2O3 Na2O 0·00 0·01 — 4·98 0·00 0·00 0·01 0·01 0·42 0·55 0·00 0·00 0·00 0·00 1·51 0·00 2·48 98·67 0·00 0·00 0·26 0·20 4·71 0·04 0·00 0·08 27·90 54·31 63·31 TiO2 core core SiO2 M87 M58 Sample: pl pl Mineral: — 5·65 0·00 0·00 0·01 0·00 0·70 0·97 0·00 0·00 0·00 0·00 1·89 0·00 2·08 99·49 0·00 0·00 0·23 0·02 0·76 19·36 0·01 0·04 0·00 0·08 34·37 0·02 44·57 core M42 pl — 4·97 0·00 0·00 0·01 0·01 0·80 0·15 0·00 0·00 0·00 0·00 1·14 0·00 2·86 99·15 0·00 0·04 0·25 0·09 9·27 3·18 0·00 0·02 0·00 0·01 21·81 0·01 64·47 core M106 pl — 5·00 0·00 0·00 0·01 0·00 0·19 0·81 0·00 0·00 0·00 0·01 1·76 0·00 2·22 99·12 0·00 0·00 0·24 0·05 2·16 16·23 0·00 0·02 0·00 0·17 32·34 0·02 47·88 core M1 pl — 4·97 0·00 0·00 0·00 0·83 0·14 0·00 0·00 0·00 0·00 0·00 1·01 0·00 2·99 98·25 0·00 0·00 0·00 13·96 1·57 0·09 0·00 0·07 0·00 0·03 18·36 0·00 64·18 core M58 kf — 5·01 0·00 0·01 0·01 0·90 0·10 0·00 0·00 0·00 0·00 0·00 1·01 0·00 2·98 98·61 0·00 0·32 0·31 15·01 1·05 0·05 0·00 0·00 0·00 0·03 18·32 0·01 63·48 core M87 kf — 5·00 0·00 0·00 0·00 0·92 0·08 0·00 0·00 0·00 0·00 0·00 1·03 0·00 2·97 99·39 0·00 0·19 0·18 15·42 0·91 0·01 0·00 0·06 0·00 0·04 18·71 0·00 63·86 core M71 kf — 5·12 0·00 0·00 0·01 0·01 0·20 1·09 0·00 0·00 0·00 0·00 1·66 0·00 2·15 94·28 0·00 0·00 0·21 0·16 2·11 20·33 0·00 0·02 0·04 0·00 28·31 0·00 43·09 core M42 me — 8·01 1·00 0·00 0·00 0·00 0·00 2·03 0·01 0·01 0·00 0·50 2·47 0·01 2·98 100·43 1·92 0·00 0·00 0·00 0·00 24·27 0·09 0·18 0·00 8·53 26·93 0·24 38·24 core M1 ep 0·42 3·99 0·00 0·00 0·00 0·00 0·02 0·89 0·56 0·01 0·41 0·08 0·10 0·01 1·91 99·06 0·00 0·00 0·00 0·01 0·31 21·44 9·72 0·31 12·63 2·71 2·29 0·31 49·30 core M87 cpx 0·44 3·99 0·00 0·00 0·00 0·00 0·01 1·01 0·51 0·02 0·40 0·04 0·04 0·00 1·96 99·39 0·00 0·00 0·00 0·00 0·10 24·36 8·89 0·59 12·38 1·52 0·85 0·07 50·54 core M42 cpx 0·23 4·01 0·00 0·00 0·00 0·00 0·04 0·94 0·64 0·02 0·19 0·12 0·18 0·01 1·87 101·04 0·00 0·00 0·00 0·00 0·56 23·62 11·62 0·62 6·00 4·20 4·00 0·22 50·19 core M1 cpx GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA JOURNAL OF PETROLOGY VOLUME 41 NUMBER 12 DECEMBER 2000 APPENDIX D Table D: End-member activities were calculated as follows, largely after Powell & Holland (1985); abbreviations after Powell & Holland (1988) Assuming ideal solid solution Amphibole ahb = 37·9XVA(XCaM2)2(XMgM3)3(XMgM1)(XAlM1) (XSiT1)3(XAlT1) aparg = 64XNaA(XCaM2)2(XMgM3)3(XMgM1)(XAlM1)(XSiT1)2(XAlT1)2 aed = 9·48XNaA(XCaM2)2(XMgM3)3(XMgM1)2(XSiT1)3(XAlT1) aftr = XVA(XCaM2)2(XFeM3)3(XFeM1)2(XSiT1)4 aanth = XVA(XMgM2)2(XMgM3)3(XMgM3)2(XSiT1)4 Cordierite acrd = (XMgM2)2 afcrd = (XFeM2)2 Feldspar aab = (XNaA) aksp = (XKA) White mica amu = 4XKA(XVM2)(XAlM1)2(XAlT1)(XSiT1)(XSiT2)2 acel = 4XKA(XVM2)(XAlM1)(XMgM1)(XSiT1)2(XSiT2)2 apa = 4XNaA(XVM2)(XAlM1)2(XAlT1)(XSiT1) (XSiT2)2 Biotite aphl = 4XKA(XMgM1)(XMgM2)2(XSiT1)(XAlT1)(XSiT2)2 aann = 4XKA(XFeM1)(XFeM2)2(XSiT1)(XAlT1)(XSiT2)2 aeast = XKA(XAlM2)(XMgM2)2(XAlT1)2(XSiT2)2 anaph = 4XNaA(XMgM1)(XMgM2)2(XSiT1)(XAlT1)(XSiT2)2 Chlorite aames = (XMgM2)4(XAlM1)2(XAlT1)2 Pyroxene adi = (XCaM2)(XMgM1)(XSi)2 ahed = (XCaM2)(XFeM1)(XSi)2 acats = (XCaM2)(XAlM1)(XSi)(XAl) Staurolite amst = (XMgM2)2(XAlM1)9(XSiT1)4 afst = (XFeM2)2(XAlM1)9(XSiT1)4 Scapolite ame = (XAlT3)3(XCaA)4 Oxides aq = asill = aky = 1 Assuming non-ideal solid solution Plagioclase = an(XCaA) an after Hoisch (1990) apy = (pyXMgM2)3(XAlM1)2(XSiT1)3 py after Berman (1990) aalm = (almXFeM2)3(XAlM1)2(XSiT1)3 alm after Berman (1990) agr = (grXCaM2)3(XAlM1)2(XSiT1)3 gr after Berman (1990) aan Garnet 1718 GOSCOMBE AND HAND P–T PATHS IN EASTERN HIMALAYA Mineral Abbreviation Mineral Abbreviation Actinolite act Magnetite mg Akermanite ak Margarite ma Almandine alm Meionitic scapolite me Andalusite and Microcline mic Antiperthite aper Monazite mon Apatite ap Monticellite mo Biotite bi Muscovite mu Calcite cc Myrmekite myr Chlorite chl Olivine ol Chloritoid ctd Opaque op Clinopyroxene cpx Perthite per Clinozoisite clz Phlogopite phl Cordierite cd Plagioclase pl Diopside di Pyrite py Dolomite dol Quartz q Epidote ep Rutile rut Ferrosilite fs Scapolite scap Fibrolite fibr Sericite seri Forsterite fo Sillimanite sill Garnet gn Sphene sph Gedrite gd Staurolite st Haematite hm Stilpnomelane stil Hercynitic spinel sp Ti-magnetite timg Hornblende hn Tourmaline tm Ilmenite ilm Tremolite tr K-feldspar kf Wollastonite woll Kyanite ky Zircon zr 1719
© Copyright 2026 Paperzz