Introduction
In continental high-pressure metamorphic terranes, it is commonly observed that quartzofeldspathic rocks such as granites, gneisses and metasedimentary rocks, which host eclogite-facies mafic rocks, record lower-pressure amphibolite-facies conditions (e.g. Chopin, Reference Chopin1984; Smith, Reference Smith1984; Young and Kylander-Clark, Reference Young and Kylander-Clark2015). The discovery of minor ultrahigh-pressure (UHP) mafic rocks within lower pressure felsic terranes led to an early debate as to whether the whole tectonic unit experienced the same conditions or whether the mafic eclogites were introduced tectonically (Lappin, Reference Lappin1977; Smith, Reference Smith1980; Cuthbert and Carswell, Reference Cuthbert and Carswell1982). The observation that, in some cases, the host quartzofeldspathic rocks also preserve some transformation at eclogite-facies conditions led to the understanding that the eclogites were metamorphosed in situ, and that water availability, residence time at depth, bulk composition and deformation were all factors that influence the equilibration of felsic rocks under eclogite-facies pressures (e.g. Austrheim, Reference Austrheim1986; Koons et al., Reference Koons, Rubie and Frueh-Green1987; Rubie, Reference Rubie1990; Cuthbert et al., Reference Cuthbert, Carswell, Krogh-Ravna and Wain2000; Wain et al., Reference Wain, Waters, Jephcoat and Olijynk2000). Debate continues as to whether the pre-existing mineral assemblages in quartzofeldspathic rocks persist metastably throughout their high-pressure metamorphic evolution, or are transformed and retrogressed pervasively leaving little record of high-pressure mineral assemblages (e.g. Wain et al., Reference Wain, Waters and Austrheim2001; Young and Kylander-Clark, Reference Young and Kylander-Clark2015).
The occurrence and timing of mineral transformations during continental subduction is relevant to geophysical models that either neglect phase transitions, or assume complete transformation over large volumes. Transformation of felsic rocks during subduction influences the density and strength of crustal material and is therefore a key component related to the exhumation of UHP material from depth (Warren, Reference Warren2013). Incomplete transformation at eclogite-facies conditions may be a common feature in ‘dry’ felsic rocks, allowing strong metastable crust to persist during subduction (e.g. Peterman et al., Reference Peterman, Hacker and Baxter2009; Krabbendam et al., Reference Krabbendam, Wain and Andersen2000; Wain et al., Reference Wain, Waters and Austrheim2001).
Zones of low strain at regional and outcrop scales provide a unique window into the early tectonic history of metamorphic rocks in these high-pressure terranes and can preserve an extensive and protracted history (e.g. Krabbendam and Wain, Reference Krabbendam and Wain1997; Wain et al., Reference Wain, Waters and Austrheim2001). Eclogite bodies, granulite-facies intermediate–mafic metaplutonic rocks (Krabbendam and Wain, Reference Krabbendam and Wain1997) and metagranitoids (Biino and Compangnoni, Reference Biino and Compangnoni1992; Gebauer et al., Reference Gebauer, Schertl, Brix and Schreyer1997; Koons et al., Reference Koons, Rubie and Frueh-Green1987) occur as low-strain enclaves preserving eclogite-facies mineral assemblages within pervasively deformed gneisses in terranes such as the Bergen Arcs (Austrheim and Griffin, Reference Austrheim and Griffin1985) and Western Gneiss Region of Norway (Mørk, Reference Mørk1985; Austrheim and Engvik, Reference Austrheim and Engvik1997); the Sulu orogen (Hirajima et al., Reference Hirajima, Wallis, Zhai and Ye1993; Wallis et al., Reference Wallis, Ishiwatari, Hirajima, Ye, Guo, Nakamura, Kato, Zhai, Enami, Cong and Banno2007; Zhang and Liou, Reference Zhang and Liou1997), the Sesia Zone (Oberhänsli et al., Reference Oberhänsli, Hunziker, Martinotti and Stern1985) and Dora Maira massif (Koons et al., Reference Koons, Rubie and Frueh-Green1987; Biino and Compagnoni, Reference Biino and Compangnoni1992) in the European Alps; and the Tso Morari Complex (Sachan et al., Reference Sachan, Mukherjee, Ogasawara, Maruyama, Ishida, Muko and Yoshioka2004) and Kaghan valley (Pognante and Spencer, Reference Pognante and Spencer1991) in the Himalaya.
This investigation describes the microstructure, mineral composition, timing and conditions of the metamorphic evolution of a suite of low-strain metagranites from the Tso Morari Complex, an ultrahigh-pressure metamorphic dome in Ladakh, NW Himalaya (Fig. 1). Whereas eclogite-facies metamorphism has been recorded in mafic bodies, very little evidence has survived in the volumetrically dominant amphibolite-facies felsic rocks. The peak pressure estimates of 26–42 kbar in this complex (de Sigoyer et al., Reference de Sigoyer, Guillot, Lardeaux and Mascle1997; Guillot et al., Reference Guillot, Sigoyer and Lardeaux1997; St-Onge et al., Reference St-Onge, Rayner, Palin, Searle and Waters2013; Chatterjee and Jagoutz, Reference Chatterjee and Jagoutz2015; Wilke et al., Reference Wilke, OBrien, Schmidt and Ziemann2015) lie within the coesite stability field, and were first confirmed by the discovery of coesite preserved in an eclogite boudin (Sachan et al., Reference Sachan, Mukherjee, Ogasawara, Maruyama, Ishida, Muko and Yoshioka2004). However, evidence of ultrahigh-pressure metamorphism in the felsic rocks of the Tso Morari Complex has been documented only as pseudomorphs after coesite in a low-strain metagranite (Bidgood et al., Reference Bidgood, Parsons, Lloyd, Waters and Goddard2021). Accompanying high-pressure metamorphic reactions have not been previously recorded in the host Puga gneiss, owing to the pervasive amphibolite-facies overprint.
The Tso Morari metagranites form metre- to kilometre-scale zones of low strain within the high-strain amphibolite-facies Puga Gneiss, which is itself of broadly granitic composition. They preserve evidence of a history of emplacement, alteration and metamorphism that spans c. 480 Myr. Observations of mineral textures and compositions in the granite are interpreted as a record of magmatic and magmatic–hydrothermal processes, overprinted by high-pressure metamorphism, and ultimately followed by Barrovian metamorphism with attendant regional deformation. Many of these early features are preserved in discrete textural domains in the rock, where local equilibrium was important, with dry conditions limiting reaction progress and aiding preservation potential.
Geological setting
The Tso Morari dome (Fig. 1) is composed of a set of stacked nappes, forming a doubly-plunging northwest–southeast elongate dome, with the stratigraphically lowest units cropping out in the centre of the dome (Steck et al., Reference Steck, Epard, Vannay, Hunziker, Girard, Morard and Robyr1998). The lowermost nappe is known as the Tso Morari Complex and contains Ordovician granite intruded into Cambrian sediments atop the Indian plate. The dominant lithology within the complex is the Puga Gneiss. Low-strain granitic domains are present throughout the complex, ranging in size from approximately one metre in diameter to map-scale domains that are over a square kilometre in area, such as in the Polokongka La region, where large boulders litter the glacial landscape, and Nuruchen, southwest of Tso Kar. The Polokongka La granite is a coarse-grained peraluminous S-type granite with an Al2O3/(CaO+Na2O+K2O) ratio in the range 0.88–1.98 (Girard and Bussy, Reference Girard and Bussy1999). Early investigators mapped the Polokongka La granite as a discrete body, presumed intrusive to the Puga Gneiss (Thakur and Virdi, Reference Thakur and Virdi1979; Fuchs and Linner., Reference Fuchs and Linner1996; Guillot et al., Reference Guillot, Sigoyer and Lardeaux1997). However, Girard and Bussy (Reference Girard and Bussy1999) and Epard and Steck (Reference Epard and Steck2008) emphasised the lack of observed intrusive relationships, the heterogeneous deformation, and the continuous transition between low-strain granite and gneiss. Near identical zircon morphology and U–Pb ages (Girard and Bussy, Reference Girard and Bussy1999) suggest that these lithologies actually belong to the same intrusive package and represent relatively low- and high-strain versions.
Both the Tso Morari metagranite and the deformed Puga Gneiss host mafic eclogite boudins, which indicates that the entirety of the Tso Morari Complex experienced eclogite-facies metamorphism (de Sigoyer et al., Reference de Sigoyer, Guillot, Lardeaux and Mascle1997; Guillot et al., Reference Guillot, Sigoyer and Lardeaux1997; Sachan et al., Reference Sachan, Mukherjee, Ogasawara, Maruyama, Ishida, Muko and Yoshioka2004). However, there is currently little reported evidence, other than prograde garnet zoning (Bidgood, Reference Bidgood2020) and quartz textures after coesite (Bidgood et al., Reference Bidgood, Parsons, Lloyd, Waters and Goddard2021), that the Puga Gneiss and granitoids underwent mineralogical transformation to eclogite-facies assemblages. The majority of the Tso Morari complex contains regional amphibolite facies Barrovian assemblages (Epard and Steck, Reference Epard and Steck2008; St-Onge et al., Reference St-Onge, Rayner, Palin, Searle and Waters2013; Palin et al., Reference Palin, St-Onge, Waters, Searle and Dyck2014).
Analytical methods
Scanning electron microscopy
Back-scattered electron imaging of textures and semi-quantitative maps of corona and pseudomorphic textures were collected from carbon coated thin sections using an FEI Quanta 650 FEG E-Scanning electron microscope at the University of Oxford, equipped with an energy dispersive spectrometer (EDS). Maps were acquired with an accelerating voltage of 15 kV using Oxford Instruments AZtec (version 3.3) acquisition software.
Electron probe microanalysis
Thin sections were made from 11 samples from the Tso Morari metagranites (Table 1). Samples TMG01, 02, 06 and 11 were selected for detailed petrological analysis. Major-element compositions of minerals that exhibit solid solutions were measured using a Cameca SX-5 field emission electron probe microanalyser (EPMA) at the Department of Earth Sciences, University of Oxford. All EPMA were conducted with 15 kV acceleration potential and 20 nA beam current. For all major elements, the on-peak count time was 30 s with 30 s of background collection. Titanium was measured with a 60 s on-peak count time and 60 s of background collection. A range of natural and synthetic oxide standards were used including: albite (Na, Al, Si); MgO (Mg); wollastonite (Ca); andradite (Fe); Mn metal (Mn); and TiO2 (Ti), and analyses were verified against secondary mineral standards. Raw analyses were automatically corrected for atomic number, absorption and fluorescence using the φρz matrix correction routine (Merlet, Reference Merlet1994). Mineral spot analyses and line profiles were taken across micas and garnet to determine the extent of intracrystalline compositional variation. Representative compositions are in Table 2, and a full list is in the Supplementary materials.
*Mineral abbreviations mostly according to Warr (Reference Warr2021): Ab = albite; Ap = apatite; Bt = biotite; Czo = Clinozoisite; Ep = epidote; Grt = garnet; Ilm = ilmenite; Kfs = K-feldspar; Ky = kyanite; Mnz = monazite; Ms/Ph = muscovite/phengite; Qz = quartz; Rt = rutile; and Zrn = zircon.
*Abbreviations – Ab = albite; Bt = biotite; Crd = cordierite; Grt = garnet; Kfs = K-feldspar; Ms = Muscovite; and Pl = plagioclase. '–' = not detected.
Laser ablation mass spectrometry
Monazite grains were imaged using back-scattered electron imaging to identify any zoning present using the FEI Quanta 650 FEG E-Scanning electron microscope (E-SEM) at the University of Oxford. The U–Pb age of monazite was measured in situ using a Nu Instruments Attom single-collector inductively coupled plasma mass spectrometer (SC-ICP-MS) at the Geochronology and Tracers Facility, British Geological Survey, using the method described in Mottram et al. (Reference Mottram, Warren, Regis, Roberts, Harris, Argles and Parrish2014). Analytical conditions and results are provided in Supplementary tables. Laser parameters were a static spot of 12 μm diameter, ablated for 30 s, at a repetition of 10 Hz, with a fluence of 2.5 J/cm2. The primary reference material used was Stern monazite (Palin et al., Reference Palin, Searle, Waters, Parrish, Roberts, Horstwood, Yeh, Chung and Anh2013; Horstwood et al., Reference Horstwood, Košler, Gehrels, Jackson, McLean, Paton, Pearson, Sircombe, Sylvester, Vermeesch and Bowring2016), and secondary reference materials Manangotry monazite (Paquette et al., Reference Paquette, Nédélec, Moine and Rakotondrazafy1994; Palin et al., Reference Palin, Searle, Waters, Parrish, Roberts, Horstwood, Yeh, Chung and Anh2013; Horstwood et al., Reference Horstwood, Košler, Gehrels, Jackson, McLean, Paton, Pearson, Sircombe, Sylvester, Vermeesch and Bowring2016) and 44069 (Aleinikoff et al., Reference Aleinikoff, Schenck, Plank, Srogi, Fanning, Kamo and Bosbyshell2006; Buick et al., Reference Buick, Lana and Gregory2011). Data reduction used an in-house spreadsheet for normalisation and IsoplotR (Vermeesch, Reference Vermeesch2018) for plotting and age calculation. Uncertainty propagation follows Horstwood et al. (Reference Horstwood, Košler, Gehrels, Jackson, McLean, Paton, Pearson, Sircombe, Sylvester, Vermeesch and Bowring2016), with the final age quoted as ± α/β, where α does not propagate systematic uncertainties, and β does. The secondary reference materials gave weighted mean 206Pb/238U ages of 553.0 ± 5.6 Ma (Manangotry) and 430.8 ± 4.3 Ma (44069), which are within 1.5% of the accepted ages.
Phase equilibrium modelling
Constraints on the pressure–temperature evolution of sample TMG01 were determined by scanning a representative area within a garnet + mica domain using the SEM-EDS system described above. The integrated bulk composition in weight % oxides, normalised to 100% anhydrous (Table 3), is SiO2 45.26; TiO2 0.20; Al2O3 30.93; FeO 12.75; MnO 0.48; MgO 3.53; CaO 0.12; Na2O 0.33; and K2O 6.40. H2O was set in excess, to fully saturate all assemblages in the investigated P–T space. The phase equilibrium model was calculated in the ten-component system MnO–Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2 (MnNCKFMASHT) using Theriak/Domino v.2015 with thermodynamic dataset ds55 (Holland and Powell, Reference Holland and Powell1998; updated August 2004), adapted by Doug Tinkham. The following activity composition models were used: Holland and Powell (Reference Holland and Powell2003) for ternary plagioclase and K-feldspar; White et al. (Reference White, Powell and Holland2007) for garnet with added spessartine and biotite with added mnbi phase component from the software list; Mahar et al. (Reference Mahar, Baker, Powell, Holland and Howell1997) and Holland and Powell (Reference Holland and Powell1998) for chlorite with added mnchl phase component, as in Tinkham et al. (Reference Tinkham, Zuluaga and Stowell2001); Coggon and Holland (Reference Coggon and Holland2002) for white mica with their pyrophyllite and margarite members excluded due to the consistent overestimation of the stability of these end-members; White et al. (Reference White, Powell, Holland and Worley2000, Reference White, Pomroy and Powell2005) for ilmenite (Mn system); and Holland and Powell (Reference Holland and Powell1998) for epidote. The following pure phases that are not part of a converted solution model were added: H2O, quartz, kyanite, albite, rutile and sphene (titanite).
Description and composition of textural domains
The Tso Morari metagranites are coarse-grained porphyritic granitic rocks that consist of quartz, alkali feldspar, white mica, biotite and albite ± garnet with accessory apatite, zircon, monazite and rutile with ilmenite rims. During deformation of the complex, the major part of the granitic protolith developed a pervasive fabric, forming the Puga Gneiss, which is dominated by assemblages of mica, albite and quartz, with less abundant K-feldspar. Quartz- and mica-rich folia wrap K-feldspar augen. The state of strain is variable, but is commonly highest towards the margins of eclogite dykes and boudins, where the gneiss may also contain abundant garnet.
Primary igneous minerals and textures are best preserved in the low-strain metagranites together with evidence for mineral reactions at different stages of its evolution, allowing a more comprehensive understanding of the geological history of the complex to be decoded. The Tso Morari metagranites contain coarse-grained phenocrysts of inclusion-rich K-feldspar reaching up to 15 mm in length, with a matrix grain size of 5 mm. The coarse-grained minerals, as seen in thin section (Fig. 2a) or in hand specimen (Fig. 2g), are best characterised as mineral domains that retain relict igneous textures. The metagranites also contain numerous enclaves, with a rounded shape and diffuse margins, reaching up to 50 cm in diameter. Although compositionally distinct, they contain mineral domains closely comparable to those in the porphyritic metagranite. The most common type consists of domain 1 (quartz, Fig. 2b) and domain 5 (garnet + mica, Fig. 2f).
There are five textural domains common to all of the granitic samples investigated (Fig. 2b–f). These are described below, in non-genetic terms.
Quartz-rich domains
The quartz-rich domain consists of both individual and clusters of bipyramidal quartz grains that are ~3–4 mm in diameter. Most quartz-rich domains display undulose extinction and/or fine deformation lamellae (Fig. 2b). Quartz clusters are made up of idiomorphic crystals that are typically joined along prism and rhombohedral facets (Fig. 3a). Both the clustered quartz grains and individual idiomorphic grains are surrounded by other mineral/textural domains of similar size, retaining primary igneous relationships with these domains (Fig. 2a,b).
Single crystals of quartz exhibit low-strain internal features such as subgrain formation and incipient recrystallisation, which is commonly confined to grain boundaries (Fig. 3b). Chessboard subgrains (~0.5 mm in diameter) are present within some of these grains (Fig. 3c). Locally, smaller subgrains have developed on a scale of <0.1 mm within the chessboard subgrains (Fig. 3c). Quartz grain boundaries are commonly gently lobate at a scale of ~0.1 mm as a result of minor grain-boundary migration (Fig. 3d,e). In contrast to the single quartz crystals preserved in all of the samples investigated from the Tso Morari Complex quartz domains in sample TMG06 from the Polokongka La region preserve palisade microstructures (Fig. 3f).
K-feldspar + white mica domains
The K-feldspar + white mica domains form the largest domains (up to 15 cm in length in outcrop) and are commonly observed to enclose the other four domain types (Fig. 2c, 4a–d). These domains commonly contain a single large K-feldspar megacryst (Or93), with microcline twinning and fine string and bleb perthitic microtextures, that has undergone near-pervasive sericitisation (Fig. 2a and 4a). Fine-grained white mica forms micaceous patches within the K feldspar (Fig. 4c) which plot along the celadonite exchange vector, with increasing Si and Mg per formula unit (pfu) (Fig. 5a), and the main cluster ranges from 6.5 to 6.9 Si pfu The K-feldspar + white mica domains also contain coarse-grained white micas replacing K-feldspar, ~0.5–1 mm in length (Fig. 2c, Fig. 4b–d) with a low phengite and high paragonite content (Fig. 5a–c). A finer-grained recrystallised rim of the coarser generation (Fig. 4b–d, 5d,f) has Si pfu > 6.7 (Fig. 5a–c).
Certain areas within K-feldspar phenocrysts represent inclusions of other minerals. Commonly these are rectangular areas in which albite is intergrown with a fine-grained white mica (Fig. 4a). This mica is texturally distinct from the even finer grained sericite-like mica that forms irregular patches after K-feldspar (Fig. 4c). These areas strongly resemble the white mica + albite domains described in the following section, although they are smaller in size.
White mica + albite domains
The third type of domain is the white mica + albite domain which consists of a fine-grained aggregate of white mica and albite that, under plane polarised light, appears darker than the K-feldspar + white mica domain and commonly exhibits zoning between dark core regions and pale rims (Fig. 2d, 4e–f). The composition of the white mica plots along the celadonite exchange vector with > 6.5 Si pfu (Fig. 5a). In some samples (e.g. TMG11) radial clusters of clinozoisite needles occur in turbid regions within this domain (Fig. 4f). In less-altered samples (e.g. TMG01), islands of plagioclase are also present with an anorthite content of 12–25 mol.%, preserved among micaceous alteration.
Micaceous domains
The micaceous domains vary in character across the samples. In the least altered samples (e.g. TMG01), this domain forms rectangular regions, commonly occupied by a large single mica crystal of muscovite or biotite (Fig. 2e, 6b, c). In other samples, the micaceous domains manifest as polycrystalline clusters of biotite and white mica. Coarse-grained muscovite in the micaceous domains has the lowest phengite and highest paragonite content of the white micas (Fig. 5a–c), comparable to coarse-grained white micas in the K-feldspar + white mica domains (Fig. 5d–f). The rims of these coarse-grained micas have recrystallised to form fine-grained phengitic white mica.
The coarse grained biotites are compositionally homogeneous with a dark brown-red colour and a Ti content of 0.25–0.38 pfu and X Mg = Mg/(Mg + Fe) of 0.21–0.39 (Fig. 6b). Geometrical arrays of rutile needles form a sagenitic texture (c.f. Shau et al., Reference Shau, Yang and Peacor1991) parallel to the cleavage planes (Fig. 6c). Trails of densely packed zircon inclusions and apatite crosscut the cleavage of the grains. The coarse-grained biotites have reaction coronas ~30 μm wide (Fig. 6a–c) consisting of phengite (6.6–6.8 Si pfu) intergrown with fingers of ilmenite, surrounded by a continuous corona of granular garnet. The garnet coronas preserve a range of compositions from Alm87Py6Grs2Sps5 to Alm55Py3Grs40Sps2 with a pronounced gradient in grossular content from 2 mol.% in contact with the biotite to 40 mol.% in contact with albite. In some instances, the white mica + ilmenite and garnet coronas are surrounded by a second rim of white mica with a higher phengite content of 6.9–7.1 Si pfu.
In some samples, coarse-grained biotites are replaced by aggregates of abundant red-brown fine-grained biotite (Fig. 6e–f) which form distinct, non-rectangular regions surrounded by thin rims of white mica and ilmenite (e.g. samples TMG06, 08 and 09). The boundaries of the micaceous aggregates are commonly defined by small (~0.5 mm) idioblastic garnets with distinct core and rim zones (Fig. 6f). In sample TMG06, large garnet crystals completely replace coarse red biotite domains and have a rim of brown biotite (Fig. 6d, e).
Garnet + mica domains
The garnet + mica domains have well-defined boundaries and are essentially composed of fine-grained white mica with small irregular clusters of garnet and kyanite (Fig. 7a–e). Coarse aggregates of white mica and biotite have occasionally grown perpendicular to the domain boundaries (Fig. 2f). Biotite is generally colourless, with lower Ti and higher Mg (X Mg = 0.41–0.55) than observed in Ti-rich biotites of the micaceous domain. Coarse-grained white micas have a phengite content of c. 6.2 Si pfu which correlates closely with the lower phengite content of the coarse-grained white mica in both the micaceous domains and the K feldspar + white mica domains.
Small garnets (0.1 mm) contain randomly oriented kyanite needles and have notably grossular-poor compositions of Alm72–79Py20–26Grs0.06–0.2Sps1.3–2.5. Garnets may be partially- or completely-replaced by fine-grained brown biotite (Fig. 7f) with a low Ti content (0.09 pfu) and a relatively high X Mg of 0.46–0.61. The replacement of garnet is most commonly observed in granite enclaves (e.g. TMG03) in which garnet + mica domains are intergrown with radial quartz rosettes. Kyanite appears to be less abundant in the fine micaceous matrix than it is in garnet, and it is absent from the coarse micas.
The composition of the garnet + mica domain from EDS is given in Table 3. Notably, the domain-scale composition represents that of cordierite (basic formula (Fe,Mg)2Al4Si5O18) with additional K and H2O. Accordingly, we interpret the garnet + mica domains as hydrated and chemically altered pseudomorphs after cordierite.
U–Pb Monazite Geochronology
Monazite is commonly associated with the garnet + mica domains, though is also found as inclusions in igneous micas, particularly biotite (Fig. 8a–d). Monazite occurs as rounded grains, ~50 μm in diameter. Back-scattered electron imaging reveals no zoning. The weighted mean 206Pb/238U age of the concordant analyses (40% of total analyses) is taken as the time of monazite formation and yields an age of 474.0 ± 11.6 Ma, MSWD = 1.5 (Fig. 8e). The high mean square weighted deviation (MSWD) indicates that even this age estimate, derived from concordant analyses, involves a variable amount of discordance, probably due to both incorporation of small amounts of common lead, and to minor disturbance of the isotopic system during younger events. The latter is compatible with the fact the Tso Morari Complex underwent high pressure metamorphism during the Eocene (e.g. de Sigoyer et al. Reference de Sigoyer, Guillot, Lardeaux and Mascle1997, St-Onge et al., Reference St-Onge, Rayner, Palin, Searle and Waters2013).
Pressure-temperature constraints
The phase equilibrium model (Fig. 9) was constructed using the bulk composition of a representative area in a garnet + mica domain, approximating a homogeneous and chemically isolated domain with abundant hydrous phases, which might be expected to record prograde metamorphism. Several features support restricted chemical exchange between this domain and others: the preservation of cordierite stoichiometry in the less mobile elements Si, Al, Fe and Mg, the very low concentration of Ca and Mn in the garnets, and the Mg-rich nature of biotite. In contrast, the microstructures in other domains, such as the narrow coronas on igneous micas, display textural and compositional disequilibrium, and are not suitable for P–T determination.
The bulk composition (Table 3) falls almost entirely in the simpler KFMASH system, with very minor amounts of MnO, Na2O, CaO and TiO2. The uncertainties on these small amounts are relatively large, however they have almost no effect on the overall phase relationships that are critical for determining pressure–temperature conditions. MnO is mainly accommodated in garnet, and TiO2 in biotite. Outside the stability field of biotite, the model predicts minor ilmenite at lower P, rutile at higher P. CaO is a very minor component in garnet, and plagioclase is predicted only at low P in the sillimanite field. Na2O is accommodated in muscovite, except at the highest pressures, where a tiny amount of Na-rich pyroxene or Na-mica is predicted (Fig. 9). Isopleths were determined for critical mineral properties, i.e. the modal abundance and pyrope content of garnet, the Si content of phengitic muscovite, and the stability limit and modal abundance of kyanite, using mineral composition ranges determined in the same domain that was used to determine the bulk composition.
Garnet is present throughout the P–T region of interest on Fig. 9, but has an interval of rapid prograde growth, shown as a pink band, accounting for at least half of the final garnet volume. Almost no further growth is predicted in the broad HP–HT fields above 600°C and 15 kbar. The observed pyrope mole fraction in garnet lies between 0.18 and 0.26. This variation corresponds to the predicted garnet growth spurt, up to the maximum volume proportion seen in the domain analysed. Essentially all kyanite growth is predicted to occur, at the expense of chloritoid, over a narrow (c. 6°C) temperature interval. This band cuts obliquely across the interval of rapid garnet growth, such that above ~26 kbar, kyanite has completed its growth at the beginning of the rapid increase in garnet, thereby allowing garnet to trap abundant kyanite inclusions. In contrast, at P <~22 kbar, kyanite does not appear until garnet is almost fully grown. Thus, the entrapment of kyanite inclusions in garnet at high pressure is consistent with the P–T path deduced from metabasic eclogites (grey arrows a and b). The interval of kyanite and garnet growth corresponds to a predicted silica content in phengite of ~6.6 per formula unit of 22 (O), matching the upper part of the range in fine-grained white mica in the domain. Therefore, we conclude that the garnet-bearing assemblage in this domain type formed during prograde eclogite-facies metamorphism at c. 550–650°C, >25 kbar.
Interpretation of magmatic-metamorphic evolution in domains
The five textural domains found in the Tso Morari granitoids can be interpreted as originating from the distribution of primary igneous minerals, i.e: (1) magmatic quartz; (2) magmatic K-feldspar, with its inclusions; (3) magmatic plagioclase; (4) magmatic muscovite and biotite; and (5) cordierite. When taken together with the record of subsequent chemical and microstructural modification, the five microdomains collectively preserve a sequence of events that includes the following four stages (Fig. 10): (1) magmatic crystallisation; (2) late-magmatic hydrothermal alteration; (3) (ultra)high-pressure metamorphism; and (4) Barrovian metamorphism with attendant regional deformation. In this section we discuss the evidence that supports a protracted magmatic–metamorphic evolution of the Tso Morari granites. The inferred mineral assemblage associated with each event is recorded in Table 4.
Magmatic crystallisation
Primary (igneous) mineral microstructures of quartz, K-feldspar, biotite and muscovite are partially preserved in most samples. The idiomorphic and bipyramidal microstructures in quartz domains indicate that many of the primary igneous quartz grains have been generally unaffected by subsequent events. Euhedral quartz grains cluster into aggregates with interstitial feldspar, suggesting that quartz was one of the first phases on the liquid line of descent. That the quartz–quartz grain boundaries form along well-developed crystal facets (Fig. 3a) indicates quartz underwent synneusis; the mutual attachment of crystals in a melt (Vance, Reference Vance1969). These sintering microstructures imply that: (1) quartz grains were crystallised early during the magma solidification process when their growth would not have been impinged by other solid phases; and (2) quartz crystals were able to move freely (i.e. in suspension) prior to attachment in low energy grain-boundary configurations (Swanson and Fenn, Reference Swanson and Fenn1986; Dyck and Holness, Reference Dyck and Holness2022).
K-feldspar megacrysts are also interpreted here as being magmatic in origin. The megacrysts contain inclusions of both single crystals and aggregates of phases which correspond to the different primary igneous domain types described in this paper. These mineral inclusions are distinctly smaller than the phases and aggregates outside the megacrysts, and they may reflect the incipient stages of crystallisation of primary igneous minerals from the magma over which the megacrysts initially grew. The now preserved K-feldspar is perthitic, occasionally showing microcline twinning (Fig. 4c), indicating that the magma was hydrous (Parsons et al., Reference Parsons and Boyd1971).
Further evidence for an earlier magmatic history is found in the coarse-grained muscovite which occurs as clusters in the matrix and is associated with coarse-grained Ti-rich biotite and polycrystalline biotite aggregates. The relatively high Ti content of biotite is indicative of a high-temperature magmatic origin (Forbes and Flower, Reference Forbes and Flower1974; Robert, Reference Robert1976; Tronnes et al., Reference Tronnes, Edgar and Arima1985; Abrecht et al., Reference Abrecht and Hewitt1988). Though the polycrystalline aggregates might not necessarily be primary igneous in nature, their high Ti content is consistent with them forming by the recrystallisation of primary coarse-grained biotite while the granite was still at elevated temperatures.
Late-magmatic hydrothermal alteration
Plagioclase and cordierite are both interpreted as being part of the original igneous assemblage, however, their compositions no longer reflect the primary mineral composition expected in igneous rocks and thus we interpret them as having undergone late- to post-magmatic hydration/metasomatic alteration. Several features are analogous to textures of late-magmatic alteration that are commonly observed in unmetamorphosed peraluminous granites. An example is the partial replacement of K-feldspar with laths or skeletal networks of muscovite, as in the Land's End granite, Cornwall, UK (Booth, Reference Booth1968) and the anatectic biotite-granite from the Trois Seigneurs Massif, Pyrenees (Fig. 11a–b) (Wickham, Reference Wickham1987). Equivalent microstructures from the Tso Morari metagranite are shown in Fig. 11c–d, indicating that these experienced a similar alteration prior to metamorphism. Moreover, the composition of these muscovites is comparable with that of magmatic white micas in the same samples (Fig. 5a–c).
The bulk composition of the garnet + mica domain resembles a cordierite composition with additional K and H2O. It is common in peraluminous granites for cordierite to destabilise upon cooling, causing it to break down during hydrothermal alteration to muscovite, phlogopite/biotite and/or chlorite (Booth, Reference Booth1968; Farina et al., Reference Farina, Stevens and Villaros2012). Unmetamorphosed examples of this can be found in the Peninsula granite of the Cape Granite Suite, South Africa (Farina et al., Reference Farina, Stevens and Villaros2012, their figure 4a), and biotite granite from the Trois Seigneurs Massif (Fig. 12a–b) (Wickham, Reference Wickham1987). In the Land's End granite of the Cornubian batholith in Cornwall, England both pinitised cordierite and secondary muscovite are associated with late- to post-magmatic alteration (Booth, Reference Booth1968; Exley et al., Reference Exley, Stone and Lees1982; Hawkes and Dangerfield, Reference Hawkes and Dangerfield1978). In the Peninsula granite, cordierite is a common phase in the enclaves within the granite, where it has been replaced by a mineral assemblage involving “aggregates of green-brown Ti-poor biotite and white mica as well as by very fine-grained mats of phyllosilicates varying under the polar microscope from grey to the typical third order blue of muscovite” (Farina et al., Reference Farina, Stevens and Villaros2012), which we interpret as being replaced partly by coarser-grained micas at higher temperature, and subsequently by pinite.
In the Tso Morari metagranites, coarse-grained muscovite and Ti-poor, Mg-rich biotite commonly form rims around the edge of the cordierite pseudomorphs, with fine-grained aggregates of phengite at the core (Fig. 7a–e, 12c–d). This texture is also commonly observed in coarse-grained enclaves, indicating that these were composed originally of cordierite and quartz. The similarity in composition between the coarse-grained white mica in the cordierite pseudomorphs, compared with the igneous white mica (Fig. 5a–c), suggests that these also were formed during the late-stage hydrothermal alteration of the granite. The biotite within the cordierite pseudomorph has a high-Mg, low-Ti composition, consistent with formation by replacement of a relatively magnesian cordierite. Subsequent metamorphism of the fine-grained micaceous pinite in the core of this domain resulted in the observed assemblage of garnet, kyanite and phengite.
Close association of monazite with pseudomorphs after cordierite may suggest that monazite formed either during crystallisation of the granite or else is associated with the breakdown of cordierite. However, cordierite does not typically contain rare earth elements (REE), so monazite is not likely to form by liberation of REE from cordierite, but is more probably associated with crystallisation of the granite. The U–Pb age of monazite of sample TMG01 is 474.0 ± 11.6 Ma (2s), and is within uncertainty of a single grain concordant U–Pb zircon age of 479 ± 2 Ma from the Polokongka La granite, which is interpreted as the age of crystallisation from a granitic melt (Girard and Bussy, Reference Girard and Bussy1999).
(Ultra)high-pressure metamorphism
The presence of palisade quartz microstructures in sample TMG06 is particularly noteworthy as this microstructure has been demonstrated recently to have formed during the coesite–quartz transition (Bidgood et al., Reference Bidgood, Parsons, Lloyd, Waters and Goddard2021) and is a cryptic indicator of ultrahigh-pressure metamorphism in felsic rocks. The fact that palisade quartz is found in only a limited number of Tso Morari granite samples, even though there are abundant quartz domains in all samples, probably demonstrates the rarity of forming and preserving this microstructure rather than a lower peak-pressure history for the Tso Morari granites.
Coronas of garnet + phengite around igneous biotite preserve evidence of biotite breakdown at high metamorphic pressures. In the Polokongka La granite this corona texture was interpreted by Girard and Bussy (Reference Girard and Bussy1999) to form from the breakdown of igneous biotite at eclogite-facies conditions, and has been documented elsewhere in the Monte Mucrone metaquartz diorites of the Sesia Zone (Oberhänsli et al., Reference Oberhänsli, Hunziker, Martinotti and Stern1985; Koons et al., Reference Koons, Rubie and Frueh-Green1987) and the Brossasco metagranite (Biino and Compagnoni, Reference Biino and Compangnoni1992) of the European Alps, as well as the Yankou metagranite of the Sulu UHP province, China (Hirajima et al., Reference Hirajima, Wallis, Zhai and Ye1993). In the Tso Morari metagranites, garnet + phengite coronas around biotite are common, although volumetrically minor, indicating that all of these samples underwent a similar eclogite-facies metamorphic history.
The fine-grained parts of the garnet + mica domains, i.e. pseudomorphs after cordierite, contain the assemblage phengite + garnet + kyanite, expected to belong to the eclogite facies. This is confirmed by the results of phase equilibrium modelling discussed above in ‘Pressure-temperature constraints’, which indicate prograde metamorphism at >25 kbar pressure. Fine-grained white mica associated with coarse micas of igneous origin have compositions that plot along the celadonite exchange vector, suggesting recrystallisation at high pressure. The mineral assemblage and textures in this domain are comparable to the assemblage preserved in metamorphosed pseudomorphs after cordierite in the French Massif Central (Lotout et al., Reference Lotout, Pitra and Pujol2017, their figure 3f).
The presence of clinozoisite needles in white mica + albite domains, i.e. former plagioclase, in TMG01 and TMG11, suggests a high-pressure metamorphic transformation. Although these domains now contain white mica rather than jadeite, the pattern and scale of the clinozoisite needles strongly resembles the distribution of (clino)zoisite within the jadeite-rich pseudomorphs after plagioclase in Monte Mucrone metaquartz-diorites (Koons et al., Reference Koons, Rubie and Frueh-Green1987).
In the Puga Gneiss, large idioblastic garnets are occasionally present in which a prograde to high-pressure metamorphic zoning pattern is preserved (Bidgood, Reference Bidgood2020). The progression from thin coronitic rims of garnet (sample TMG01), to granular or skeletal replacement of red-brown biotite (sample TMG04) to euhedral garnets (sample TMG06) in increasingly metamorphosed samples, could represent a more complete transformation of biotite at eclogite-facies conditions, facilitated by the availability of fluid at high pressure.
Barrovian metamorphism with attendant deformation
Pressure and temperature paths from mafic eclogites in the Tso Morari Complex predict steep, approximately isothermal decompression from the eclogite-facies peak to amphibolite-facies conditions (Chatterjee and Jagoutz, Reference Chatterjee and Jagoutz2015; de Sigoyer et al., Reference de Sigoyer, Guillot, Lardeaux and Mascle1997; Guillot et al., Reference Guillot, Sigoyer and Lardeaux1997; St-Onge et al., Reference St-Onge, Rayner, Palin, Searle and Waters2013; Wilke et al., Reference Wilke, OBrien, Schmidt and Ziemann2015). The Tso Morari metagranites, in contrast to the Puga Gneiss, show limited evidence for mineral growth and re-equilibration under amphibolite-facies conditions. The extent of recrystallisation (e.g. of micas) appears to correlate with increasing strain.
Some insight can be obtained from the phase equilibrium modelling of the former cordierite domains (Fig. 9). During decompression to c. 15 kbar, very little change is predicted in mineral proportions or composition, however some re-equilibration may have occurred in the sample to account for Si contents in phengite decreasing to c. 6.4 pfu, and the corresponding increase in Al is perhaps coupled with an observed (and predicted) decrease in the amount of kyanite in the mica matrix, compared with that trapped in garnet. Under amphibolite-facies conditions (<15 kbar), biotite becomes stable, however almost no devolatilisation is predicted on decompression at c. 650°C. The implied reaction is essentially Grt + Ms = Bt + Ky + Qz, which (at constant K2O) conserves H2O, however there is no evidence that this reaction occurred. Instead, the observed replacement of garnet by biotite seen in many samples is not accompanied by new growth of an aluminium silicate (Ky or Sil). This simple replacement would require addition of external fluid, and a further implication would be that K2O was also introduced in this fluid. In the absence of deformation, relatively little fluid is involved, as replacements are not always complete and the microstructures are not significantly modified.
The breakdown of garnet and its replacement by biotite could also suggest that the Puga Gneiss, strongly deformed under amphibolite-facies conditions and potentially fluid-saturated, given the extent of amphibolite-facies retrogression shown by the basic bodies it contains, has lost much of any garnet it developed at higher pressures.
The observation that plagioclase (white mica + albite) domains contain only hints that they were transformed to jadeite-bearing aggregates suggests either that much plagioclase was not converted to jadeite, or that jadeite and quartz were transformed to albite (and white mica) on decompression. It is also possible that plagioclase was extensively sericitised at low temperature during post-magmatic alteration. The latest features developed in the Tso Morari metagranites may relate to the low-strain deformation recorded in many samples. Deformation temperatures can be estimated by comparing observed quartz microstructures to qualitative temperature ranges for quartz deformation regimes (e.g. bulging, subgrain rotation and grain-boundary migration; Stipp et al., Reference Stipp, Stunitz, Heilbronner and Schmid2002). High-temperature quartz microstructures in the form of chessboard subgrains are observed in granitoid samples from the Polokongka La (Fig. 3b,c,e) region of the Tso Morari Complex. These microstructures are indicative of low-strain deformation at temperatures >550–600°C (Kruhl, Reference Kruhl1996; Wallis et al., Reference Wallis, Parsons and Hansen2019). Amoeboid shapes along grain boundaries are evidence for low-strain features of active, though limited, deformation in the grain-boundary migration regime at ~550°C (Stipp et al., Reference Stipp, Stunitz, Heilbronner and Schmid2002). Metagranite samples from the Polokongka La region also display tabular subgrains internal to the grain, suggesting recovery associated with dislocation glide and climb at somewhat lower temperatures (Passchier and Trouw, Reference Passchier and Trouw2005). Peak deformation temperatures of >600°C are consistent with electron back-scattered diffraction analysis of quartz in the Puga Gneiss across the Tso Morari Complex by Dutta and Mukherjee (Reference Dutta and Mukherjee2021).
Discussion
The incomplete transformation of metagranitoids at (ultra)high pressures
The scarcity of mineralogical evidence preserved in felsic rocks that have experienced (ultra)high pressures implies that either: (1) (ultra)high-pressure assemblages are not readily produced (i.e. low-pressure minerals persist in a metastable manner); or (2) (ultra)high-pressure assemblages form but are readily transformed to lower pressure assemblages upon decompression. The evidence we present here for the preservation of both igneous and late-magmatic hydrothermal alteration indicates the importance of hypothesis 1 in the absence of penetrative deformation, regardless of uncertainty over the extent of high-P breakdown of plagioclase to jadeite.
The Tso Morari metagranites underwent partial transformation in the eclogite facies. Examples of igneous textural domains, preserved despite later replacement of the igneous minerals, are found in low-strain igneous rocks such as the previously mentioned Monte Mucrone metaquartz diorite (Dal Piaz et al., Reference Dal Piaz, Hunziker and Martinotti1972; Oberhänsli et al., Reference Oberhänsli, Hunziker, Martinotti and Stern1985; Koons et al., Reference Koons, Rubie and Frueh-Green1987, Zucali et al., Reference Zucali, Corti, Roda, Ortolano, Visalli and Zanoni2021), the Brossasco metagranites (Schertl et al., Reference Schertl, Schreyer and Chopin1991; Biino and Compangnoni., Reference Biino and Compangnoni1992) and the Yankou metagranite (Wallis et al., Reference Wallis, Ishiwatari, Hirajima, Ye, Guo, Nakamura, Kato, Zhai, Enami, Cong and Banno2007), as well as in the Kokchetav Massif in Kazakhstan (Zhang et al., Reference Zhang, Liou, Ernst, Coleman, Sobolev and Shatsky1997). Additionally, in the Flaktraket mangerite–granulite complex in the Western Gneiss region of Norway, >95% of the complex is thought to have preserved metastable igneous textures and coronas around the magmatic orthopyroxene during eclogite-facies conditions (Wain et al., Reference Wain, Waters and Austrheim2001).
There are some marked differences between the Tso Morari metagranites and high-pressure metagranites from other terranes. For example, omphacite is expected to be the equilibrium pyroxene for felsic bulk compositions, and in the case of local equilibrium, plagioclase domains may be transformed to jadeite-rich aggregates (e.g. Koons et al., Reference Koons, Rubie and Frueh-Green1987). However, sodic pyroxene was not observed in the feldspar domains of the Tso Morari granites. Nevertheless, textures preserved in the plagioclase domains, including the presence of arrays of clinozoisite needles similar to those in transformed plagioclase, and the lack of primary twin lamellae, suggest that most of the primary plagioclase did not survive high-pressure metamorphism. Therefore, on balance, we consider that clinopyroxene was present at high pressures in at least some of these rocks, however it was subsequently overprinted during exhumation.
Static transformation in granite can be limited to within igneous mineral-scale domains, where a local effective bulk composition is dominant and there is very little communication between domains. For example, Koons et al. (Reference Koons, Rubie and Frueh-Green1987) observed jadeite cores and omphacite rims within pyroxene, zoisite and quartz-bearing domains replacing igneous plagioclase in the Monte Mucrone metaquartz diorites. Initially, the formation of jadeitic clinopyroxene after plagioclase was controlled by the effective bulk composition of the igneous Na-rich plagioclase. Subsequent mineral changes during growth of clinopyroxene were interpreted to be a result of partial equilibration of local mineral domains to the bulk composition as the equilibration volume enlarged and recrystallisation was promoted by progressive deformation, rather than directly due to changes in pressure and temperature (Koons et al., Reference Koons, Rubie and Frueh-Green1987).
Garnet + phengite + ilmenite coronas around biotite are a common feature of partially transformed granites in high-pressure terranes. Koons et al. (Reference Koons, Rubie and Frueh-Green1987) provided a detailed interpretation of this texture in the Monte Mucrone metaquartz diorites, where the asymmetrical compositional zoning of grossular in garnet coronas surrounding biotite is attributed to the influence of the breakdown of plagioclase into jadeite + quartz, which acts as a source for CaO and Al2O3. The rate-determining step in this reaction is thought to be the diffusion of Al2O3 and CaO from the plagioclase domain to the biotite rim and also accounts for the asymmetric zoning pattern seen in the garnet coronas. Koons et al. (Reference Koons, Rubie and Frueh-Green1987) suggest that the non-completion of the biotite breakdown reaction could be due to the decrease in diffusive flux owing to the growth of garnet and muscovite coronas, lengthening diffusive pathways and therefore reducing the chemical potential gradients. Alternatively, the production of Al2O3 and CaO may have ceased as the plagioclase breakdown reaction ran to completion. The presence of asymmetric compositional zoning in garnet may therefore be an indication of plagioclase breakdown in the eclogite facies.
Partial transformation within the Tso Morari metagranites was limited by the availability of H2O at high pressure, leading to corona textures where reactions did not run to completion. There was little communication between igneous domains, indicated by the essentially isochemical transformation of hydrous domains such as metamorphosed cordierite pseudomorphs. The absence of ‘typical’ eclogite-facies mineral assemblages such as jadeite + quartz within plagioclase domains could be attributed to lack of preservation owing to the re-equilibration during Barrovian metamorphism, or to extensive pre-metamorphic sericitisation.
Conclusions
Low-strain metagranite domains are found within the Tso Morari Complex, and are observed in this investigation to preserve textural and compositional evidence of early and late magmatic–hydrothermal features, followed by later eclogite-facies metamorphism. The Tso Morari metagranites therefore provide a unique window into these early histories. The microstructure of the igneous domains, the assemblages and the mineral compositions in the Tso Morari metagranites reveal a complex evolution in which a peraluminous S-type granite was emplaced during the Ordovician at 474.0 ± 11.6 Ma and hydrothermally altered during the late stages of emplacement.
Eclogite-facies metamorphism led to the partial transformation of igneous biotite to form garnet + phengite coronas and the complete transformation of hydrous pseudomorphs after cordierite to garnet, phengite and kyanite during high-pressure metamorphism, indicating prograde metamorphism at >25 kbar. The relatively low abundance of diagnostic high-pressure textures in the Tso Morari metagranites and preservation of igneous minerals such as mica indicates that metamorphic transformations were limited by fluid availability, and controlled by mineral-scale compositional domains in the rock. In cases where transformation occurred in the eclogite facies, the absence of characteristic minerals such as jadeite may be attributed to overprinting during Barrovian metamorphism where H2O was more readily available, leading to partial re-equilibration of original feldspar and cordierite domains.
The metagranites stand in marked contrast to their host rock, the high-strain Puga Gneiss, which contains amphibolite-facies mineral assemblages, fabric development and recrystallisation having erased almost all of its earlier history. The Tso Morari metagranites draw attention to the role of low strain and lack of fluid availability in preserving the original features and revealing the P–T history of continental basement rocks.
Acknowledgements
This work was funded by the Natural Environmental Research Council, grant number NE/L002612/1 awarded to AKB. Fieldwork to Ladakh was undertaken in 2016 and 2017 as part of the PhD of AKB and was partially funded by the Geological Society of London Mike Coward fund, the Mineralogical Society of Great Britain and Ireland, Edinburgh Geological Society, University College Oxford and the Royal Geographical Society.
The illustrated biotite granite sample from the Trois Seigneurs massif, eastern Pyrenees, was collected in 1984 by DJW on a field excursion led by R.S.J. Sparks and S.M. Wickham, supported by the Department of Earth Sciences, University of Cambridge. The sample L08-04 (TMG11) was collected in 2008 by M.P. Searle. Thanks to Jon Wade and Phil Gopon from Oxford University for help and training on the SEM and EMPA.
We would like to thank Michael Anenburg and two anonymous reviewers for their constructive comments on this manuscript.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1180/mgm.2022.121