1. Introduction
Tracing first-order lithospheric structures in collisional orogens is challenging due to their complex architecture and tectono-metamorphic history. One example is the position of a major structure in the most classical collisional orogen: the Main Central Thrust zone (MCTz), a first-order tectonic boundary running along the Himalayan belt for ∼ 2400 km along strike (Fig. 1). It separates the metamorphic core (Greater Himalayan Sequence (GHS)) in the hanging wall from the Lesser Himalayan Sequence (LHS) in the footwall. Since the MCTz was defined by Heim and Gansser (Reference Heim and Gansser1939), its position, age and role in the tectonic evolution of the belt are still debated (Searle et al. Reference Searle, Law, Godin, Larson, Streule, Cottle and Jessup2008; Martin, Reference Martin2017a, Reference Martin2017b; Mukhopadhyay et al. Reference Mukhopadhyay, Chakraborty, Trepmann, Rubatto, Anczkiewicz, Gaidies, Dasgupta and Chowdhury2017, Carosi et al. Reference Carosi, Montomoli and Iaccarino2018; Montemagni et al. Reference Montemagni, Iaccarino, Montomoli, Carosi, Jain and Villa2018, Reference Montemagni, Carosi, Fusi, Iaccarino, Montomoli, Villa and Zanchetta2020; Chakraborty et al. Reference Chakraborty, Mukul, Mathew and Pande2019; Iaccarino et al. Reference Iaccarino, Montomoli, Montemagni, Massonne, Langone, Jain, Visonà and Carosi2020). According to Godin et al. (Reference Godin, Grujic, Law and Searle2006) and Catlos et al. (Reference Catlos, Lovera, Kelly, Ashley, Harrison and Etzel2018, Reference Catlos, Perez, Lovera, Dubey, Schmitt and Etzel2020), the age of the MCTz spans between 25 and 18–16 Ma down to nearly 3 Ma (Catlos et al. Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001). The localization and timing of the MCTz play a major role in the definition of the tectonic models envisaged for the exhumation of the belt. Several models (Parsons et al. Reference Parsons, Ferré, Law, Lloyd, Phillips and Searle2016a, Reference Parsons, Law, Lloyd, Phillips and Searle2016b and references therein) are based on the contemporaneous activity of the MCT and the South Tibetan Detachment System (STDS), which separates the GHS from the upper Tethyan Himalayan Sequence (THS), so that the age and the structural position of the MCT are fundamental issues (Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013; Montemagni et al. Reference Montemagni, Montomoli, Iaccarino, Carosi, Jain, Massonne and Villa2019, Reference Montemagni, Carosi, Fusi, Iaccarino, Montomoli, Villa and Zanchetta2020). Montomoli et al. (Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015) pointed to the occurrence of a regional-scale structure (High Himalayan Discontinuity (HHD)) within the GHS. Carosi et al. (Reference Carosi, Montomoli and Iaccarino2018) demonstrated that the MCT and HHD are often confused in the literature for the time span from 25–26 Ma to 16 Ma (Godin et al. Reference Godin, Grujic, Law and Searle2006). The univocal position and dating of structural and metamorphic discontinuities in the GHS became key points for unravelling the tectonic evolution of the orogenic belt. A strand of the STDS, the Chame detachment in the Marsyangdi valley (Coleman, Reference Coleman1998; Searle & Godin, Reference Searle and Godin2003; Searle, Reference Searle2010) (Fig. 2a), has been interpreted as a ramp of the HHD (Walters & Kohn, Reference Walters and Kohn2017) so that even the kinematics of shear zones of the upper GHS is challenged. Dating deformation events in tectonites is a major challenge (Müller et al. Reference Müller, Aerden and Halliday2000; Williams & Jercinovic, Reference Williams and Jercinovic2002, Reference Williams and Jercinovic2012; Carosi et al. Reference Carosi, Montomoli, Iaccarino, Benetti, Petroccia and Simonetti2022; Dumond et al. Reference Dumond, Mahan, Goncalves, Williams and Jercinovic2022) and in situ radiometric dating of mineral (re)crystallization within different fabrics helps to unravel the deformation history (Di Vincenzo et al. Reference Di Vincenzo, Carosi and Palmeri2004; Villa, Reference Villa2016). Metamorphic minerals often show chemical zoning of major and trace elements, and a large spread of dates is obtained, even in a single mineral (Villa et al. Reference Villa, Hermann, Muntener and Trommsdorff2000; Rubatto et al. Reference Rubatto, Schaltegger, Lombardo, Colombo and Compagnoni2001, Reference Rubatto, Chakraborty and Dasgupta2013; Gibson et al. Reference Gibson, Carr, Brown and Hamilton2004; Kohn et al. Reference Kohn, Wieland, Parkinson and Upreti2005; Pyle et al. Reference Pyle, Spear, Wark, Daniel and Storm2005; Pyle, Reference Pyle2006; Corrie & Kohn, Reference Corrie and Kohn2011; Kohn, Reference Kohn2016). It is not trivial to link the obtained dates to the different generation of structures. Nevertheless, the complexity in age and mineral chemistry of geochronometers is the key to unravel the history of deformation and metamorphism (Engi, Reference Engi2017; Schulz, Reference Schulz2021; Dumond et al. Reference Dumond, Mahan, Goncalves, Williams and Jercinovic2022; Imayama et al. Reference Imayama, Hoshino, Keewok and Kawabata2022). Approaches combining in situ U-Th-Pb data, with textural position and geochemical composition of minerals, are useful to date geologic events (e.g., Williams et al. Reference Williams, Jercinovic and Terry1999, Reference Williams, Jercinovic and Hetherington2007; Gibson et al. Reference Gibson, Carr, Brown and Hamilton2004; Pyle, Reference Pyle2006; Rubatto et al. Reference Rubatto, Chakraborty and Dasgupta2013; Kohn et al. Reference Kohn, Engi and Lanari2017). In this paper, we focus on laser ablation split stream (LASS) dating (Cottle et al. Reference Cottle, Kylander-Clark and Vrijmoed2012; Kylander-Clark et al. Reference Kylander-Clark, Hacker and Cottle2013) of monazite in metasedimentary rocks affected by high-temperature (amphibolite facies) ductile shear in the metamorphic core of the Himalaya in the Manaslu massif (Central-Western Nepal) (Figs. 1, 2).
Figure 1. (a) Geological map of the Himalaya (modified after Searle, et al. Reference Searle, Law, Godin, Larson, Streule, Cottle and Jessup2008; Weinberg, Reference Weinberg2016; Searle, Reference Searle, Treloar and Searle2019) and (b) its geographic position. MFT, Main Frontal Thrust; MBT, Main Boundary Thrust; MCT, Main Central Thrust; STDS, South Tibetan Detachment System; P, Peshawar basin; S, Sutlej basin.
Figure 2. (a) Geological map of the Marsyangdi river valley. (b) Geological map of the Budhi Gandaki river valley (after Parsons et al. 2016 and our observations) showing the trace of the Jagat Shear Zone (JSZ) and sample locations. On the right-hand side: stereoplot (Wulff net, lower hemisphere) of the mylonitic foliation and grain/object lineation (red dots) of the JSZ.
The HHD was first recognized in Western Nepal (Carosi et al. Reference Carosi, Montomoli, Rubatto and Visonà2010; Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015; see Carosi et al. Reference Carosi, Montomoli and Iaccarino2018, Reference Carosi, Montomoli, Iaccarino and Visonà2019 for a review), but gaps exist remain in our understanding of its along-strike continuation to the east (Carosi et al. Reference Carosi, Montomoli and Iaccarino2018; Waters, Reference Waters, Treloar and Searle2019). In this framework, rocks from the Marsyangdi and Budhi Gandaki valleys in Central Nepal (Fig. 2) were investigated to detect the occurrence and the structural position of a shear zone, which can be correlated with the HHD, and to constrain its age by in situ U-Th-Pb dating of monazite.
2. Geological setting
2.a. The Himalayas
The Himalayas (Fig. 1) developed after the initial collision between the Indian and Eurasian plates at ∼59–55 Ma (Le Fort Reference Le Fort1975; Green et al. Reference Green, Searle, Corfield and Corfield2008; Hu et al. Reference Hu, Garzanti, Wang, Huang, An and Webb2016). The orogenic belt is composed of different tectonic units derived from the northern part of the Indian Plate, which were deformed, metamorphosed and accreted at different structural levels south of the Indus-Tsangpo Suture Zone (Hébert et al. Reference Hébert, Bezard, Guilmette, Dostal, Wang and Liu2012; Searle, Reference Searle, Treloar and Searle2019). The zone contains remnants of the Neo-Tethys Ocean that was located between the Indian and Eurasian plates. Five main litho-tectonic units, forming the backbone of the belt, are recognized. To the south, Cenozoic synorogenic sedimentary rocks of the Sub-Himalayan sequence are tectonically overlaid by metasedimentary and meta-igneous rocks of the LHS (Upreti & Le Fort, Reference Upreti and Le Fort1999; Hodges, Reference Hodges2000). A km-thick ductile to ductile-brittle shear zone, the MCTz (Heim & Gansser, Reference Heim and Gansser1939; Gansser, Reference Gansser1964; Searle et al. Reference Searle, Law, Godin, Larson, Streule, Cottle and Jessup2008; Martin, Reference Martin2017a, Reference Martin2017b), separates the mostly greenschist – with minor lower amphibolite-facies rocks of the LHS from the tectonically overlying GHS, which is mainly composed of middle to upper amphibolite facies metasedimentary and meta-igneous rocks (Le Fort, Reference Le Fort1975; Hodges, Reference Hodges2000; Carosi et al. Reference Carosi, Montomoli and Iaccarino2018 and references therein) and represents the metamorphic core of the orogenic belt. Deformation with a top-to-the-S kinematics, linked to the MCTz, affected both the bottom of the GHS and the upper part of the LHS. In the lower part of the GHS garnet- to kyanite-bearing paragneiss and micaschist with subordinate calc-schist, quartzite, impure marble and migmatite occur. The middle portion of the GHS is made of calcsilicate gneiss and marble with minor pelitic and psammitic levels (Fig. 2). The uppermost part of the GHS consists of orthogneiss and aluminosilicate-bearing migmatites. Crustally derived, Oligo-Miocene leucogranite intrusions (Searle, Reference Searle2010, Reference Searle2013; Visonà et al. Reference Visonà, Carosi, Montomoli, Peruzzo and Tiepolo2012, Weinberg, Reference Weinberg2016), such as the Manaslu granite (Guillot et al. Reference Guillot, Le Fort, Pêcher, Barman and Aprahamian1995; Cottle et al. Reference Cottle, Lederer and Larson2019), are abundant in the upper part of the GHS. Tectonically above the GHS, unmetamorphosed sedimentary rocks (Frank & Fuchs, Reference Frank and Fuchs1970; Garzanti, Reference Garzanti1999) with minor amphibolite facies rocks (Antolín et al. Reference Antolín, Appel, Montomoli, Dunkl, Ding, Gloaguen, El Bay, Poblet and Lisle2011; Dunkl et al. Reference Dunkl, Antolín, Wemmer, Rantitsch, Kienast, Montomoli, Ding, Carosi, Appel, El Bay, Xu, von Eynatten, Gloaguen and Ratschbacher2011; Montomoli et al. Reference Montomoli, Iaccarino, Antolin, Appel, Carosi, Dunkl, Lin and Visonà2017) of the THS are present. The tectonic boundary between GHS and THS is represented by a system of low-angle ductile shear zones to brittle faults with a normal top-to-the-N or -NE kinematics, referred to the STDS (Burg et al. Reference Burg, Brunel, Gapais, Chen and Liu1984; Burchfiel et al. Reference Burchfiel, Chen, Hodges, Liu, Royden, Changrong and Xu1992; Carosi et al. Reference Carosi, Lombardo, Molli, Musumeci and Pertusati1998; Searle et al. Reference Searle, Simpson, Law, Parrish and Waters2003, Kellett et al. Reference Kellett, Grujic, Warren, Cottle, Jamieson and Tenzin2010, Reference Kellett, Cottle and Larson2019; Searle, Reference Searle2010; Iaccarino et al. Reference Iaccarino, Montomoli, Carosi, Montemagni, Massonne, Langone, Jain and Visonà2017b; Nania et al. Reference Nania, Montomoli, Iaccarino, Leiss and Carosi2022). The occurrence of the STDS at the top of the GHS and of the MCTz at its bottom, as two coeval north-dipping shear zones with opposite kinematics, played a key role in shaping the Himalayas and exhuming the GHS. These processes are, however, complicated, particularly since ductile contractional shear zones within the GHS were found in several transects across the belt (Carosi et al. Reference Carosi, Montomoli, Rubatto and Visonà2010; Martin et al. Reference Martin, Ganguly and Decelles2010, Reference Martin, Copeland and Benowitz2015; Corrie & Kohn, Reference Corrie and Kohn2011; Imayama et al. Reference Imayama, Takeshita, Yi, Cho, Kitajima, Tsutsumi, Kayama, Nishido, Okumura, Yagi, Itaya and Sano2012, Reference Imayama, Hoshino, Keewok and Kawabata2022; Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015; Ambrose et al. Reference Ambrose, Larson, Guilmette, Cottle, Buckingham and Rai2015; Cottle et al. Reference Cottle, Larson and Kellett2015; He et al. Reference He, Webb, Larson, Martin and Schmitt2015; Khanal et al. Reference Khanal, Robinson, Kohn and Mandal2015; Larson et al. Reference Larson, Ambrose, Webb, Cottle and Shrestha2015; Wang et al. Reference Wang, Rubatto and Zhang2015, Reference Wang, Zhang, Liu, Zhang, Wang, Rai and Scheltens2016; Zeiger et al. Reference Zeiger, Gordon, Long, Kylander-Clark, Agustsson and Penfold2015; Agustsson et al. Reference Agustsson, Gordon, Long, Seward, Zeiger and Penfold2016; Iaccarino et al. Reference Iaccarino, Montomoli, Carosi, Massonne and Visonà2017a; Walters & Kohn, Reference Walters and Kohn2017; Goscombe et al. Reference Goscombe, Gray and Foster2018; Chakraborty et al. Reference Chakraborty, Mukul, Mathew and Pande2019; Waters, Reference Waters, Treloar and Searle2019; Shrestha et al. Reference Shrestha, Larson, Martin, Guilmette, Smit and Cottle2020; Benetti et al. Reference Benetti, Montomoli, Iaccarino, Langone and Carosi2021) and related to a regional scale tectono-metamorphic discontinuity known as HHD. The HHD occurs in the middle part of the GHS close to the boundary between migmatite and kyanite-bearing paragneiss, with hanging wall rocks showing a higher degree of melting compared to the footwall rocks and with a decoupling in the prograde and retrograde metamorphic history (Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015). The activity of the HHD is constrained between 28–27 and 17 Ma and facilitated exhumation of the upper part of the GHS before the initiation of the MCT (Carosi et al. Reference Carosi, Montomoli, Rubatto and Visonà2010; Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015).
2.b. Localization of the MCT
Mapping, localization, timing and tectono-metamorphic evolution of the MCTz are still debated (Martin, Reference Martin2017a). Researchers have proposed several criteria for defining the MCTz including lithological, timing of activity, structural and metamorphic criteria. Examples are (i) a ‘sharp’ litho-tectonic boundary (e.g. Heim & Gansser, Reference Heim and Gansser1939) between the medium- to high-grade gneiss (i.e. GHS) and the medium- to low-grade quarzite (LHS); (ii) a metamorphic criterion to identify the MCTz close to the kyanite-in isograd (Le Fort, Reference Le Fort1975; Pêcher, Reference Pêcher1989); (iii) a protolith boundary, where the MCTz separates rocks with different Nd and Sr isotopic signatures (e.g. Ahmad et al. Reference Ahmad, Harris, Bickle, Chapman, Bunbory and Prince2000) and with different maximum depositional age based on zircon U-Pb populations (Parrish & Hodges, Reference Parrish and Hodges1996) and (iv) a structural criteria (Searle et al. Reference Searle, Law, Godin, Larson, Streule, Cottle and Jessup2008), where the MCTz is mapped based on strain criteria, that is, the recognition of a strain gradient on distinct deformation localization highlighted by field and microstructural observations (Law et al. Reference Law, Stahr, Francsis, Ashley, Grasemann and Ahmad2013; Parsons et al. Reference Parsons, Ferré, Law, Lloyd, Phillips and Searle2016a, Reference Parsons, Law, Lloyd, Phillips and Searle2016b, Reference Parsons, Law, Searle, Phillips and Lloyd2016c; Larson et al. Reference Larson, Godin and Price2010, Reference Larson, Cottle, Lederer and Rai2017; Hunter et al. Reference Hunter, Weinberg, Wilson, Luzin and Misra2019). Frequently, a contractional top-to-the-S high-strain zone corresponds to the base of an inverted metamorphic sequence (Searle et al. Reference Searle, Law, Godin, Larson, Streule, Cottle and Jessup2008; Iaccarino et al. Reference Iaccarino, Montomoli, Montemagni, Massonne, Langone, Jain, Visonà and Carosi2020). Moreover, several authors mapped two different thrusts, MCT-I and MCT-II (Arita, Reference Arita1983) or Munsiari and Vaikrita thrusts (Valdiya et al. Reference Valdiya, Paul, Tarachandra-Bhakuni and Upadhyay1999) bounding strongly sheared rocks of the MCTz. Probably, a combination of several of them is required to localize the MCTz. The problem of localization of the MCT is particularly evident in our study area as exemplified by the variable locations of the MCT in different studies (Colchen et al. Reference Colchen, Le Fort and Pêcher1986; Larson et al. Reference Larson, Godin and Price2010, Reference Larson, Cottle and Godin2011; Martin et al. Reference Martin, Copeland and Benowitz2015; Parsons et al. Reference Parsons, Law, Searle, Phillips and Lloyd2016c).
2.c. Study area
The Marsyangdi river and Budhi Gandaki river valleys are located in the Annapurna region of Central-Western Nepal, on the western and the eastern side, of the Manaslu massif and expose the upper part of the LHS, a complete section of the GHS and the lower part of the THS (Colchen et al. Reference Colchen, Le Fort and Pêcher1986; Pêcher, Reference Pêcher1989, Reference Pêcher1991; Coleman, Reference Coleman1996, Reference Coleman1998; Coleman & Hodges, Reference Coleman and Hodges1998; Parsons et al. Reference Parsons, Law, Searle, Phillips and Lloyd2016c). In the Marsyangdi river valley, the LHS, separated from the GHS by the MCTz (Fig. 2), is made of chlorite- to garnet-bearing quartzite. The MCTz is represented by mylonitic schist, quarzite and paragneiss (Pêcher, Reference Pêcher1989, Reference Pêcher1991; Catlos et al. Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001, Reference Catlos, Lovera, Kelly, Ashley, Harrison and Etzel2018; Larson et al., Reference Larson, Godin and Price2010; Catlos, Reference Catlos, Cemen and Catlos2021; Martin et al. Reference Martin, Copeland and Benowitz2015; Gibson et al. Reference Gibson, Godin, Kellett, Cottle and Archibald2016) defining a 2 to 4 km-thick high-strain zone with top-to-the-S kinematics. According to Pêcher (Reference Pêcher1989) and Catlos (Reference Catlos, Cemen and Catlos2021), the top of the MCTz is localized just above the kyanite-in isograd. Recent garnet-based pressure-temperature (P-T) paths were presented by Catlos et al. (Reference Catlos, Lovera, Kelly, Ashley, Harrison and Etzel2018), following Catlos et al. (Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001), for rocks across the MCTz suggesting an increase of both pressure and temperature traced structurally upwards. Coleman (Reference Coleman1998) and Catlos et al. (Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001) presented monazite U-Pb and Th-Pb ages, respectively, suggesting that the shearing along the MCTz was active between 25–22 Ma and c. 3 Ma. However, the interpretation of the U-Th-Pb dates, obtained by these authors, is dubious due to lack of detailed characterization of the monazite chemical composition (Kohn, Reference Kohn2016). Recent monazite dating by Gibson et al. (Reference Gibson, Godin, Kellett, Cottle and Archibald2016) (their sample MSY-03) yielded ages between 32 and 26 Ma for prograde conditions and between 18–17 Ma and ∼ 13 Ma for the exhumation of MCTz rocks.
The GHS in the Marsyangdi river valley (Fig. 2a) is classically subdivided into several formations (Le Fort, Reference Le Fort1975; Colchen et al. Reference Colchen, Le Fort and Pêcher1986; Pêcher & Le Fort, Reference Pêcher, Le Fort, Le Fort, Cholchen and Montenat1986), now more properly referred as units (Searle & Godin, Reference Searle and Godin2003). The lowest Unit I consists of garnet- to kyanite-bearing paragneiss and micaschist (Unit Ia;. 2a) and sillimanite-bearing paragneiss and migmatite (Unit Ib; Fig. 2). Close to the village of Jagat, within Unit I, a ductile shear zone, recognized for the first time in this contribution and described in detail below, splits (Fig. 2a) Unit I (as defined by Corrie & Kohn, Reference Corrie and Kohn2011) into Unit Ia and Unit Ib. Layered calcsilicate, marble and amphibole-biotite-bearing gneiss (Fig. 2a) constitute Unit II. Unit III is represented by orthogneiss, migmatite and biotite-gneiss (Fig. 2a). According to Coleman (Reference Coleman1996), a thin interlayered portion of Unit II splits Unit III into structurally upper and lower parts. Unit IV at the top of the GHS is defined by diopside-amphibole calcsilicate and marble (Fig. 2a). In the studied area, the THS is tectonically separated from the GHS by the STDS which is represented by mylonitic calcsilicate and metapelite sandwiched between a lower shear zone, referred as the Chame detachment as it crops out close to the village of Chame (Coleman, Reference Coleman1996; Parsons et al. Reference Parsons, Ferré, Law, Lloyd, Phillips and Searle2016a, Reference Parsons, Law, Lloyd, Phillips and Searle2016b, Reference Parsons, Law, Searle, Phillips and Lloyd2016c; Fig. 2a), and the upper Phu detachment (Fig. 2a) (Searle & Godin, Reference Searle and Godin2003; Searle, Reference Searle2010). Low-grade marble and metapelite (Schneider & Marsh, Reference Schneider and Marsh1993) represent the main lithologies of the THS (Fig. 2a). In the upper part of the Marsyangdi valley, leucogranite dykes and sills, likely linked to the activity of the Manaslu pluton (22–19 Ma, Guillot et al. Reference Guillot, Hodges, Le Fort and Pêcher1994; Cottle et al. Reference Cottle, Lederer and Larson2019), are common. A debate exists about the nature, kinematics and tectonic meaning of the Chame detachment, which has been previously interpreted as a ductile branch of the STDS (Coleman, Reference Coleman1996, Coleman & Hodges, Reference Coleman and Hodges1998; Searle, Reference Searle2010, Parsons et al. Reference Parsons, Ferré, Law, Lloyd, Phillips and Searle2016a, Reference Parsons, Law, Lloyd, Phillips and Searle2016b, Reference Parsons, Law, Searle, Phillips and Lloyd2016c). In contrast, Walters and Kohn (Reference Walters and Kohn2017) proposed a reverse-sense (top-to-the-S) kinematic framework for this structure, mainly based on petrological and titanite petrochronological data and interpreted it as an intra-GHS thrust, active in the time span of 25–17 Ma. Thus, it is crucial for bracketing the geological evolution in the area to unambiguously determine the kinematic nature of the Chame detachment.
A similar setting is recognizable in the Budhi Gandaki river valley (Larson et al. Reference Larson, Godin and Price2010, Reference Larson, Cottle and Godin2011) on the eastern side of the Manaslu massif (Fig. 2b). The main differences to the Marsyangdi valley are the virtual absence of calcsilicate and marble of Unit II (Colchen et al. Reference Colchen, Le Fort and Pêcher1986) and the occurrence of carbonate-bearing lithologies within the MCTz (Colchen et al. Reference Colchen, Le Fort and Pêcher1986; Nania, Reference Nania2021). A ductile shear zone within Unit I, 1 km south of the village of Philim (Fig. 2b), crops out (Pêcher & Fort, Reference Pêcher and Fort2012).
In both valleys, the metamorphic grade increases from biotite-grade in the LHS to sillimanite+K-feldspar grade in the migmatitic part of the GHS (Fig. 2). GHS metamorphism in both areas occurred during the Eocene to Miocene (Pêcher, Reference Pêcher1989; Coleman, Reference Coleman1998; Catlos et al. Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001; Larson et al. Reference Larson, Cottle and Godin2011; Gibson et al. Reference Gibson, Godin, Kellett, Cottle and Archibald2016, Walters & Kohn, Reference Walters and Kohn2017).
3. The Jagat Shear Zone in the GHS
Detailed work in the Marsyangdi and Budhi Gandaki river valleys (Figs. 1, 2) resulted in the recognition and mapping of a contractional, top-to-the-S shear zone close to the sillimanite-in isograd (Figs. 2, 3) near Jagat and Philim (Figs. 2, 3). This shear zone is here referred to the Jagat Shear Zone (JSZ).
Figure 3. 3D geological model of the Marsyangdi river valley showing the location of the Jagat Shear Zone. Same legend as in Fig. 2.
In the Marsyangdi river valley, the JSZ is a few hundred metres thick and affected garnet-aluminosilicate-bearing gneisses of Unit I. On the field, a sharp transition is discernible from tectonites showing tight to isoclinal folds and related axial plane foliation (in both the hanging wall and the footwall rcoks), towards a high strain zone, showing evidences of strain localization with the developing of a mylonitic foliation and evidences of grain-size reduction. Shear planes dip to the N/NE and the associated stretching lineation plunges to the NE. Shear sense indicators (Fig. 4a–d), such as C-S fabrics, shear bands, asymmetric tails around porphyroclasts, mica and mineral fishes, asymmetric boudins and rotated garnet grains, indicate a top-to-the-S/SW sense of shear (Fig. 4a–d) pointing to a thrust-sense shear zone (Figs. 2, 4a–d). Biotite, sillimanite and minor white mica constitute the main anastomosing schistosity, wrapping around grains of feldspar, garnet and kyanite. Sillimanite is observed around garnet and kyanite oriented along the shear bands. Rare staurolite locally rims large kyanite porphyroblasts (Fig. 4e) and is microstructurally associated with sillimanite. Quartz in sheared rocks of the JSZ shows evidence of high-temperature grain boundary migration (GBMII) and recrystallization leading to chessboard extinction patterns (Fig. 4f; Kruhl, Reference Kruhl1996; Law, Reference Law2014). Plagioclase experienced internal ductile deformation discernable by deformation twins (Fig. 4f) and asymmetric myrmekite. Deformation microstructures in the sheared rocks suggest a high-temperature deformation regime (T > 630–650°C; Kruhl, Reference Kruhl1996; Passchier & Trouw, Reference Passchier and Trouw2005) in agreement with the synkinematic mineral assemblage.
Figure 4. Meso- and microstructures of Jagat Shear Zone (JSZ) rocks and footwall rocks. JSZ: (a) Outcrop-scale sheared rocks showing quartz-feldspar sigmoids (arrow) pointing to a top-to-the-S sense of shear. (b) S-C fabric and asymmetric sigmoidal quartz-feldspar lithons (arrow) in sheared garnet-aluminosilicate-bearing paragneiss, pointing to a top-to-the-S sense of shear. (c) Photomicrograph showing a S-C fabric in sheared paragneiss (crossed nicols: XPL). (d) Sillimanite-bearing shear bands in sheared paragneiss (parallel nicols, //P). (e) Large kyanite (Ky) surrounded by staurolite (St), and minor sillimanite (Sil) (//P). (f) Chessboard extintion in quartz (XPL), pointing to deformation temperatures higher than 630–650°C. Footwall of the JSZ: (g) garnet (Grt) and large Ky (i.e. KyI) porphyroblast in paragneiss (//N). (h) Small St associated with small Ky (i.e. KyII) within the paragneiss (//N).
The JSZ separates sillimanite-garnet-bearing gneiss and migmatite (with relict kyanite) in the hanging wall (Unit Ib) from the kyanite-garnet-rutile-bearing gneiss (Unit Ia) (Fig. 4g, h) in the footwall. In the footwall rocks, two texturally different generations of kyanite I (see also Vannay & Hodges, Reference Vannay and Hodges1996) were observed: kyanite forms large (up to centimetre-sized) porphyroblasts (Fig. 4g) and later kyanite II occurs as small crystals along feldspar-quartz (or garnet) grain boundaries, locally associated with small euhedral staurolite (Fig. 4h). Moreover, larger anhedral staurolite (relict of prograde metamorphism) was observed, as already reported by Coleman (Reference Coleman1998).
The JSZ crops out also in the Budhi Gandaki valley just south of Philim, where a zone of intense deformation occurs affecting garnet-aluminosilicate-bearing gneisses of Unit I. In this valley, footwall rocks are represented mainly by garnet-kyanite-bearing paragneiss, with sporadic relicts of staurolite and minor quartzite, whereas hanging wall rocks are made mostly by sillimanite-bearing migmatitic gneiss with minor calcsilicate levels. Also in this location, shear planes dip to the N/NE and the associated stretching lineation plunges to the NE. The highly deformed rocks show abundant sigmoidal quartz (Fig. 5a) and C-S fabric at the mesoscale. The spaced schistosity is defined by biotite and white mica wrapping around garnet and feldspar grains. Garnet shows an internal foliation at high angle with respect to the external schistosity. Thus, we regard garnet as an intertectonic (Passchier & Trouw, Reference Passchier and Trouw2005) mineral. At the microscale, C-S fabric, foliation fishes and rotated garnet grains point to a top-to-the-S sense of shear (Fig. 5b). Quartz shows evidence of dynamic recrystallization through grain boundary migration (GBMII) (Stipp et al. Reference Stipp, Stunitz, Heilbronner and Schmid2002) with the presence of window structures. Plagioclase is ductilely deformed as indicated by the development of deformation twins. Previously, along the Budhi Gandaki valley, Parsons et al. (Reference Parsons, Law, Searle, Phillips and Lloyd2016c) highlighted on their map S foliations and C planes in correspondence of the shear zones. Pêcher and Fort (Reference Pêcher and Fort2012) hypothesized a thrust-sense shear zone by the occurrence of high-strain levels containing sheath folds.
Figure 5. Meso- and microstructures of Jagat Shear Zone rocks in the Budhi Gandaki valley. (a) Sheared gneiss with quartz sigmoids and rootless folds South of the Philim village. (b) C-S fabric and foliation fish (arrow) pointing to a top-to-the-S sense of shear (XPL).
In order to unravel the JSZ tectono-metamorphic history, selected samples were studied to derive P-T-D-time paths. In addition, a study of the petrofabrics of tectonites along the STDs profile was conducted in order to assess the Chame detachment kinematics and the relationship with the JSZ. The position of selected samples, their main assemblages and microstructural features, including the type of analyses applied, are summarized in Table 1.
Table 1. Summary of the main features (e.g. structural position, rock type) of selected samples (see Fig. 2 for sample locations)
Minerals are listed according to decreasing modal amount.
AvePT, THERMOCALC average PT; AveP, THERMOCALC average P; FW, footwall; HW, hanging wall; JSZ, Jagat Shear Zone; STDS, South Tibetan Deatachment System; Zr-in-Rt, zirconium in rutile thermometer.
4. Rocks of the JSZ: mineral chemistry and geothermobarometry
Thin sections from selected field-oriented samples across a structural profile, cut perpendicular to the main foliation and parallel to the main lineation (approximating the XZ plane of the finite strain ellipsoid), were investigated using a CAMECA SX100 EMP equipped with 5 wavelength-dispersive spectrometers (WDS), hosted at the (now closed) Institut fur Mineralogie und Kristallchemie (Universität Stuttgart). The energy-dispersive spectrometer of this EMP was used for qualitative identification of minerals. Synthetic and natural minerals, glasses and pure oxides were used as standards for the characterization of mineral compositions. Counting times were 20 s at the peak and on the background each. The applied accelerating voltage of 15 kV and a beam current of 30 nA and 10 nA were used to analyse garnet and other minerals (micas, feldspar, staurolite, ilmenite), respectively. Details on the analytical procedure and related errors are reported by Massonne (Reference Massonne2012). In order to assess possible garnet zoning, multiple spots and mineral traverses were conducted together with qualitative garnet WDS X-ray elementary (for Mn, Mg, Fe, Ca and Y) maps at 60 nA with the aforementioned EMP. For analysing the Zr content in rutile, an accelerating voltage of 15 kV, a beam current of 100 nA, and 100 s counting time on peak and background each were selected. Structural formulae of minerals, and related end-member activities (used for geothemobarometry), were calculated with the A-X software (downloaded from T. Holland’s personal web-page).
For geothermobarometrical estimates, a multi-equilibrium (MET) method, based on the Average P-T and Average P (at a fixed temperature) subprograms of THERMOCALC (Powell & Holland, Reference Powell and Holland1994, Reference Powell and Holland2008; version 3.33), was applied. Analysed garnets show slightly variable chemical zoning patterns, with flat, or nearly flat, profiles. Minor zoning in Ca was observed, and a sharp increase of Mn at the outermost rims is a common feature (see representative X-ray maps in supplementary material S1). We regard such profiles in garnet as related to diffusion modification/homogenization due to operation of both retrograde exchange and net-transfer reactions (Kohn & Spear, Reference Kohn and Spear2000). Following Yakymchuk and Godin (Reference Yakymchuk and Godin2012), analyses of garnet with the lowest Mn (and highest Mg) contents, corresponding to mantle (or near rim) position, were combined with rim compositions of matrix minerals (not in direct contact with the garnet rim) such as micas, K-feldspar, plagioclase and staurolite (where present, see Table 1 ‘inferred peak mineral assemblage’). Representative mineral data are given in Table S1. For Zr-in-rutile geothermometry, the pressure sensitive calibration by Tomkins et al. (Reference Tomkins, Powell and Ellis2007) was applied considering the pressure range of 0.7–0.9 GPa and 0.9–1.1 GPa for mylonites of the JSZ and rocks of its footwall, respectively. Only pristine rutile grains (i.e. lacking ilmenite overgrowth) in the matrix were analysed. An uncertainity of ±30°C is assigned to the temperature estimates (Tomkins et al. Reference Tomkins, Powell and Ellis2007). P-T paths were derived combining microstructural and petrographic observations and P-T data from geothemobarometry with predicted mineral stability (and sequence of mineral growth) from petrogenetic grids. Metamorphic conditions for the JSZ, its hanging wall and footwall, are summarized in Table 1 and Fig. 6.
Figure 6. (a) P-T results of selected samples taken from different structural positions along the studied transect. Grey boxes refer to P-T estimates given in Catlos et al. (Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001) (see Table 1 and main text for details, Fig. 2 for sample location). KFMASH equilibria plotted with the GIBBS software (Spear & Menard Reference Spear and Menard1989) using the SPaC2007 database. NKFMASH curves are after Spear et al. (Reference Spear, Kohn and Cheney1999). (b) Reconstructed P-T paths of different GHS portions, based on geothermobarometry and petrographic observations.
Geothermobarometric estimates suggest that JSZ mylonites equilibrated at P-T conditions of c. 0.7–0.8 GPa and c. 700°C, close to the kyanite-sillimanite boundary. Footwall rocks equilibrated at higher P of c. 1.0–1.1 GPa and lower T (<700°C) in the kyanite stability field. Hanging wall rocks returned P-T values of c. 1.1 GPa and c. 735°C, falling in the melt- and kyanite-bearing field (Fig. 6), although kyanite could be lacking in thin sections.
5. In situ U-Th-Pb monazite petrochronology
In situ U-Th-Pb monazite dating was applied to obtain the timing of metamorphism and deformation along the JSZ. Microtextural position, internal features, back-scattered electron (BSE) images and chemical zoning of monazite were evaluated with the aforementioned EMP. Monazite from two selected samples (Table 1) was analysed by LASS inductively coupled plasma mass spectrometry. The analytical procedure is reported in Kylander et al. (Reference Kylander-Clark, Hacker and Cottle2013) with modifications as reported in McKinney et al. (Reference McKinney, Cottle and Lederer2015). U-Th/Pb and trace element data are listed in Tables S2 and S3. Representative BSE images, WDS chemical maps, spot locations and corresponding dates, are given in supplementary material S2, S3, S4 and S5. Sample MA16-90 (Fig. 2) is a mylonitic JSZ paragneiss with white and dark micas, garnet, kyanite and sillimanite. Sample MA16-33 is a K-feldspar-sillimanite-biotite-garnet bearing migmatitic gneiss in the hanging wall of the JSZ in the Budhi Gandaki river valley (Fig. 2), showing late white mica (as flakes and as sericite replacing sillimanite) and chlorite after biotite crystallization, likely linked to a late stage of fluid ingress and alteration.
In both samples, monazite sizes range from 30 to 200 µm in diameter. In sample MA16-90, three analysed monazite grains occur in quartz-feldspar domains: one grain is included in kyanite (Mnz 9); other grains are in biotite oriented along the mylonitic foliation. In sample MA16-33, five monazite grains are in quartz-feldspar domains and another five are hosted in mica domains. Xenotime is lacking in sample MA16-90, whereas it is observed as a rare matrix mineral in sample MA16-33. In the latter sample, in one case xenotime is in contact with high-Y monazite (Mnz 8). We would like to stress that, despite not observed, xenotime may have been present in the early-prograde stages and later on totally consumed (e.g., Pyle & Spear, Reference Pyle and Spear1999, Reference Pyle and Spear2003; Spear & Pyle, Reference Spear and Pyle2010).
All monazite grains are mainly concentrically zoned (Catlos, Reference Catlos2013), with minor patchy and oscillatory zoning in sample MA16-33 (Supplementary material S4, S5), with respect to Y, Th and heavy rare-earth elements (HREE). Cores are enriched in Th and variably depleted in Y and HREE, whereas the rims and mantles are enriched in Y and HREE and depleted in Th (Fig. 7). Interestingly, in sample MA16-33, high-Y and HREE rim domains correspond (in some circumstances) to portions where oscillatory zoning in Th is present (Fig. 7, Supplementary material S4, S5), a feature intepreted as linked to melt-crystallization (Rubatto et al. Reference Rubatto, Chakraborty and Dasgupta2013) during cooling. Moreover, in this sample low- to moderate-Y and HREE monazite overgrowth, as discontinuous outermost rims or forming patchy domains, are sometimes observed.
Figure 7. Monazite petrochronological results from sample MA16-90 and MA16-33. Plot of Y (ppm) vs. 208Pb/232Th dates (a) sample MA16-90 (b) sample MA16-33. Plot of Gd/Yb ratio vs. 208Pb/232Th ages (c) sample MA16-90 (d) sample MA16-33. Representative compositional maps (Y and Th) of monazites from both samples are included in the upper part of the figure. Warm colours point to higher concentration.
There is a good correlation between 208Pb/232Th dates, Y content and Gd/Yb ratio. The dates of low-Y (and HREE) and high Gd/Yb cores range from 31 to 27 Ma, whereas high-Y and low Gd/Yb rims yielded dates between 27 and 13 Ma in sample MA16-90 (Figs. 7, 8, Supplementary material S2, S3). In sample MA16-33, we observed in part similar trends of Y, HREE and Gd/Yb in monazite. Monazite cores yielded dates from 39 to 29 Ma; rims and mantles gave dates from 28 Ma to 13–10 Ma. The older group of dates (39–29 Ma) are linked to domains with low-Y and HREE and higher Gd/Yb values (Fig. 7). The younger group of dates shows more complicated trends. An increase in Y and HREE, together with a decreasing in Gd/Yb, is observed until ∼17–15 Ma, whereas youngest dates (13–10 MA) are associated together with decrease in Y (and HREE) and variable to increasing Gd/Yb values.
Figure 8. Histogram of 208Pb/232Th dates vs. number of analyses for samples MA16-90, MA16-33 (this study),and MSY-03 (Gibson et al. Reference Gibson, Godin, Kellett, Cottle and Archibald2016) with arrows marking prograde and retrograde evolutions according to chemical zoning of monazite and age of the HHD (Jagat Shear Zone) and Main Central Thrust (MCT).
Due to the good correlation among dates and monazite chemical composition, we interpret the obtained dates as follows according to Kohn et al. (Reference Kohn, Wieland, Parkinson and Upreti2005), Braden et al. (Reference Braden, Godin and Cottle2017) and Soucy La Roche et al. (Reference Soucy La Roche, Godin, Cottle and Kellett2018). Interpretation of the Y and HREE zoning in monazite rely on the observation assumptions that garnet is the main silicate controlling the budget of these elements during metamorphisms (Pyle & Spear, Reference Pyle and Spear1999, Reference Pyle and Spear2003; Gibson et al. Reference Gibson, Carr, Brown and Hamilton2004; Kohn et al. Reference Kohn, Wieland, Parkinson and Upreti2005; Spear & Pyle, Reference Spear and Pyle2010 and references therein).
The low-Y, high-Th and high Gd/Yb monazite cores are linked to the prograde burial path and formed contemporaneously with garnet, starting at 39 Ma (a minimum age) and ceasing at 28–27 Ma. It is worth noting that the span of time ∼ 29–27 Ma is retrieved also from monazite (Mnz 9) hosted in kyanite in sample MA16-90, showing intermediate Y contents (Fig. 7). The chronological interpretation of such a grain is tricky. Indeed, Gibson et al. (Reference Gibson, Carr, Brown and Hamilton2004) reported the possible formation of coeval kyanite and monazite, with intermediate Y content (see their reaction 4), around peak pressure along a prograde clockwise P-T path.
On the contrary, Soucy La Roche et al. (Reference Soucy La Roche, Dyer, Zagorevski, Cottle and Gaidies2022) demonstrated that monazite inclusions in kyanite may form on the retrograde path during the partial replacement of kyanite by white mica based on the observation that monazite intersects fractures filled with white mica, which could have served as fluid pathways. We documented cracks in kyanite around hosted monazite, but there is no white mica or other evidence for late replacement of kyanite or fluid infiltration; therefore, we interpret the 29–27 Ma Mnz 9 as a primary inclusion that constrains the timing of peak P conditions.
High-Y and low Gd/Yb monazite rims crystallized along the decompression and retrograde path, where the breakdown of garnet to mica and sillimanite and melt crystallization in the migmatite (MA16-33) occurred in the time span of 28–27 Ma down to 13 Ma. Matrix monazite (Mnz 8) in MA16-33 is interpreted to have co-crystallized at 19–14 Ma with xenotime that may have formed due to garnet breakdown (Spear et al. Reference Spear, Kohn and Cheney1999; Gibson et al. Reference Gibson, Carr, Brown and Hamilton2004; Spear & Pyle, Reference Spear and Pyle2010). The youngest group of ages (down to 10 Ma), in light of their different more variable chemistry, and lack of clear chemical trends are interpreted as related to late, locally, fluid-mediated monazite alteration/recrystallization (Kohn et al. Reference Kohn, Wieland, Parkinson and Upreti2005).
6. Chame detachment: petrofabric data
In order to determine the kinematics of the Chame detachment, which is crucial for any tectonic interpretation of the southernmost JSZ, tectonite samples were characterized and compared with rocks from a strand of the STDS, the Phu detachment, which occurs in the northern Marsyangdi river valley (see Table 1, Figs. 2a, 9). Microstructures, integrated with field observations, were investigated with standard polarization microscopy. Calcite and quartz crystallographic preferred orientations (CPO) were studied at the Geowissenschaftliches Zentrum Göttingen through an X-ray Texture Goniometer (X’Pert Pro MRD_DY2139 by PANalytical). X-ray texture analyses were performed on areas of ∼ 3 cm2 in rock slices (∼1cm thick) cut parallel to the XZ and YZ planes of the finite strain ellipsoid. For each sample, crystallographic orientation data from both planes were combined to recalculate the Orientation Distribution Function through MTEX Toolbox (available for download at: https://mtex-toolbox.github.io/) of the Matlab software. The adopted lattice parameters were: a = b = 4.913 Å, c = 5.405 Å for quartz and a = b = 4.988 Å, c = 17.062 Å for calcite. As both mineral phases are trigonal, crystal symmetry ‘312’ (limited to Laue groups) was adopted. The main crystallographic directions of quartz and calcite were displayed on equal area lower hemisphere pole figures with a projection plane normal to the foliation and parallel to the lineation, using the kernel de la Vallée Poussin function (Schaeben, Reference Schaeben1997) with a halfwidth of 10° (Hielscher & Schaeben, Reference Hielscher and Schaeben2008).
Figure 9. Mesoscopic kinematic indicators of the Chame detachment: a) delta-type porphyroclast in orthogneiss in the uppermost part of the GHS; top-down-to-the-NW sense of shear; b) mica fish in sheared leucogranite; top-down-to-the-NW sense of shear. c) Photomicrograph of sample MC17-02 (deformed gneiss); foliation is highlighted by biotite and shape preferred orientation (SPO) of both feldspar and quartz (XPL); d) photomicrograph of sample MC17-17 (XPL), a pure marble with a grain SPO; e) equal area lower hemisphere pole figures of the main crystallographic elements for quartz and calcite textures in samples MC17-02 and MC17-17 (upper and lower row, respectively).
At the outcrop-scale, mica fishes and rotated mantled porphyroclasts, in sheared leucogranite (Fig. 9a, b), indicate a top-down-to-the-NW sense of shear that is cryptic at the microscale. We selected two samples (MC17-02 and MC17-17) for CPO analyses. Sample MC17-02 is a sheared gneiss (Fig. 2) and sample MC17-17 a marble (Fig. 2) collected at the Chame detachment and in the footwall of the Phu detachment (near the village of Brathang), respectively. Sample MC17-02 (Fig. 9c) shows a heterogeneous mineral distribution, where shape preferred orientation (SPO) of biotite and plastically deformed quartz and feldspar define a spaced foliation.
Quartz CPO patterns (Fig. 9e) show a mid-low CPO intensity, displayed by a classic ‘texture index’ of 1.6 (Bunge, Reference Bunge1982) and a M-Index of 0.01 (designed for large ODF datasets, Skemer et al. Reference Skemer, Katayama, Jiang and Karato2005). Quartz [c]-axis maxima (on [0001] pole figure), between the Z- and X-directions of the finite strain ellipsoid, are consistent with type I cross girdle with a fabric opening angle fabric of ca. 72°, corresponding to a deformation temperature of ∼580° ± 50°C (Law, Reference Law2014). The asymmetric pattern suggests a non-coaxial flow (Lister & Hobbs, Reference Lister and Hobbs1980; Wallis, Reference Wallis1992; Law & Johnson, Reference Law, Johnson, Law, Butler, Holdsworth, Krabbendam and Strachan2010). Dextral asymmetry of [c]-axis maxima, compared to the geographical reference system, is coherent with a top-to-the-N non-coaxial deformation.
Marble MC17-17 (Fig. 9d) shows a rough continuous foliation defined by a SPO of fine-grained biotite and slightly elongated calcite. Calcite is well-interconnected, representing the weak matrix (Handy, Reference Handy1994). Lobate grain boundaries and unimodal distribution indicate a grain boundary mobility mechanism (GBM regime) (Fig. 9d). Type II e-twins, typical for temperatures below 300 °C (Burkhard, Reference Burkhard1993; Ferrill et al. Reference Ferrill, Morris, Evans, Burkhard, Groshong and Onasch2004), occur. A texture index of 1.2 (Bunge, Reference Bunge1982), consistent with a M-Index of 0.03, supports a mid-CPO intensity (Skemer et al. Reference Skemer, Katayama, Jiang and Karato2005) for calcite. Well-defined patterns and asymmetric sinistral [c]-axis maxima in the Z-direction, with further maxima on Y-direction (central girdle), are recognized (Fig. 9e). Sinistral asymmetric poles to the {e}-planes (on {
${01\overline{18}}$
} pole figure) support {e}-twinning development in a top-to-the-N non-coaxial flow (Fig. 9e). The {r}-planes maxima (on {
${10\overline{14}}$
}) and the {f}-planes (on {
${01\overline{12}}$
}) define small circles, coherent with the coexistence of intracrystalline deformation. Intracrystalline slip and twinning support a deformation temperature range of 300–800 °C (De Bresser & Spiers, Reference De Bresser and Spiers1997). According to the c-axis-based method of Wenk et al. (Reference Wenk, Takeshita, Bechler, Erskine and Matthies1987), calcite asymmetric texture (of ca. 7°) results from a general flow with low simple shear contribution (∼26%).
7. Discussion
7.a. JSZ: P-T-D-t path and implications
The in situ U-Th-Pb dating of zoned monazite in rocks from the Marsyangdi and Budhi Gandaki valleys yielded dates from ∼ 39 to 13 Ma (Fig. 8). Similar dates have been obtained by Coleman (Reference Coleman1998), Catlos et al. (Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001, Reference Catlos, Lovera, Kelly, Ashley, Harrison and Etzel2018), Gibson et al. (Reference Gibson, Godin, Kellett, Cottle and Archibald2016) and Larson et al. (Reference Larson, Cottle and Godin2011) highlighting younger dates from north to south. The zoned monazite studied in this work allowed us to better unravel the metamorphic and deformation history. Chemical zoning and LASS results indicate that the older monazite dates are related to prograde metamorphism between 39 and 28–27 Ma, whereas the younger dates, obtained from high-Y and HREE rims, constrain the start of the retrograde path at 28–27 Ma. This path recorded by rocks from the hanging wall of the JSZ is a result of the reverse kinematics of the shear zone forcing the hanging wall rocks to move upward and southward. The youngest monazite ages suggest that in a time interval lasting 11–10 My (from 28 Ma to 17 Ma) the footwall of the shear zone moved downward while the hanging wall was simultaneously moved upward and exhumed. The hanging wall underwent retrograde metamorphism and exhumation over a duration of 11–10 My before the movement of the footwall started.
The occurrence of the thrust-sense, top-to-the-S shear zone, located at the boundary between sillimanite-bearing gneiss and underlying kyanite-bearing rocks, as well as the timing of this zone is in agreement with the structural position and age of the HHD determined by Montomoli et al. (Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015) and Iaccarino et al. (Reference Iaccarino, Montomoli, Carosi, Massonne and Visonà2017a) in Western Nepal, Iaccarino et al. (Reference Iaccarino, Montomoli, Carosi, Massonne, Langone and Visonà2015) and Wang et al. (Reference Wang, Rubatto and Zhang2015, Reference Wang, Zhang, Liu, Zhang, Wang, Rai and Scheltens2016) in Central Nepal, Chakraborty et al. (Reference Chakraborty, Mukul, Mathew and Pande2019) and Benetti et al. (Reference Benetti, Montomoli, Iaccarino, Langone and Carosi2021) in India and Waters (Reference Waters, Treloar and Searle2019) in Nepal. Metamorphic conditions of selected samples are consistent (Fig. 6) with both the inferred ‘peak’ mineral assemblage (and related mineral stability constrained from petrogenetic grids) and previous estimates by Coleman (Reference Coleman1998) and Catlos et al. (Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001). Only the hanging wall sample (MA16-86), containing sillimanite, plots in the kyanite stability field. This could be due to the lack of equilibration, the determination of P-T conditions from mineral that formed/equilibrated prior to the decompression and/or the difficulty to consider P-T errors (here: ±0.1 GPa, ±30°C, see parag. 4) adequately. We note that our derived P-T conditions are consistent with P-T data of calcsilicates collected just above MA16-86 (see Walters & Kohn, Reference Walters and Kohn2017). Therefore, we favour the first explanation ‘P-T conditions prior to the decompression’. Our kinematics and P-T-D-t path of the JSZ are in agreement with data on the HHD in nearby areas such as the Kali Gandaki valley (Iaccarino et al. Reference Iaccarino, Montomoli, Carosi, Massonne, Langone and Visonà2015). As a consequence, the occurrence of the JSZ in the Manaslu massif fills a gap in tracing the HHD to the east, confirming its regional extent.
7.b. How is strain accomodated in the HHD?
The occurrence of a regional scale HHD in the GHS showing a thickness ranging from a few dozens of metres up to nearly 4 km poses a problem if it is a single shear zone or a wide zone of deformation concentrated at several shear zones acting at different times. Montomoli et al. (Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015) firstly documented the occurrence of the HHD between 28 and 17 Ma by the presence of sheared rocks, kinematic indicators, geothermobarometry on hanging wall and footwall rocks and U-Th-Pb in situ dating of monazite. In Western Nepal, the shear zone is characterized by thick mylonites (up to 4 km; Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013; Iaccarino et al. Reference Iaccarino, Montomoli, Carosi, Massonne and Visonà2017a). U-Th-Pb dating on monazite yielded a relatively wide time span of shearing from 28 to 17 Ma.
In Central Nepal, Shrestha et al. (Reference Shrestha, Larson, Martin, Guilmette, Smit and Cottle2020) using detailed petrochronology on rocks from both the Banuwa and Sinuwa thrusts were able to document the occurrence of two main in-sequence thrust-sense shear zones active within the time span of the activity of the HHD: starting at ∼ 24 Ma for the Sinuwa thrust and ∼ 21 Ma for the lower Banuwa thrust. However, the kinematics of the Sinuwa thrust is only inferred, whereas the kinematics of the Banuwa thrust is still under discussion: Martin et al. (Reference Martin, Ganguly and Decelles2010, Reference Martin, Copeland and Benowitz2015) proposed normal-sense kinematics, whereas Corrie and Kohn (Reference Corrie and Kohn2011) suggested a thrust-sense kinematics active from 27 to 22 Ma. Shrestha et al. (Reference Shrestha, Larson, Martin, Guilmette, Smit and Cottle2020) proposed an initial thrust-sense, top-to-the-S Banuwa thrust, which was later re-activated as a normal-sense, top-down-to-the-N, shear zone (at ∼10 Ma).
In NE Nepal, Ambrose et al. (Reference Ambrose, Larson, Guilmette, Cottle, Buckingham and Rai2015) proposed the occurrence of five ‘cryptic’ shear zones in the GHS and out of sequence thrusts, indirectly detected on the basis of integrated pseudosection modelling and monazite petrochronology of paragneiss, active from 31 to 12 Ma, mostly falling within the timing of the HHD. However, both the occurrence and the kinematics of the ‘cryptic’ shear zones were only inferred. In addition, they often lack direct evidence of strain localization (high strain zones) such as mylonites or sheared rocks and related meso- and microscale kinematic indicators.
Eastern Nepal is characterized by a major number of (inferred) metamorphic discontinuities with respect to Western and Central Nepal (Cottle et al. Reference Cottle, Larson and Kellett2015; Larson et al. Reference Larson, Ambrose, Webb, Cottle and Shrestha2015) where one or two tectono-metamorphic discontinuities have been documented so far (Carosi et al. Reference Carosi, Montomoli, Iaccarino, Massonne, Rubatto, Langone, Gemignani and Visonà2016, Reference Carosi, Montomoli and Iaccarino2018) on the basis of structural and field observations in addition to petrological and geochronological data.
The GHS was affected by polyphase tectonics with several generations of folds, foliations and intra-GHS discontinuitites so that much caution is requested in attributing a tectonic significance and, moreover, a kinematics of samples showing different P-T data without the help of field evidences and structural analysis.
However, the occurrence of several metamorphic discontinuities in the GHS, if further confirmed on structural basis, does not prevent the occurrence of a major HHD between sillimanite-bearing migmatites, with a higher degree of melting, and kyanite-bearing gneiss below (Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015).
Further investigations are required to better understand if the HHD is a unique shear zone or if it is made of several (in-sequence?) shear zones and what type of shear zone the HHD is (Type I, II, III and IV sensu Fossen & Cavalcante, Reference Fossen and Cavalcante2017). This latter point is particular interesting (to be the target of future research), since it would allow to better understand how strain is distributed within large, regional scale, shear zones affecting medium- to high-grade metamorphic rocks. As suggested by Montomoli et al. (Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013), it is necessary to perform detailed meso- and microstructural analyses on mylonites or sheared rocks combined with detailed thermobarometry of the hanging wall and footwall rocks as well as petrochronology in order to be able to discriminate real ages of shearing in the wide time span of 28–17 Ma. However, to exclude possible ambiguities in the interpretations, we would like to stress that shear zones and their sense of shear should be recognized at the meso- and micro-scale by structural analysis and not only inferred by differences in P-T conditions or ages.
7.c. The Chame detachment
Walters and Kohn (Reference Walters and Kohn2017) identified a strand of the HHD in the upper part of the Marsyangdi valley interpreting the Chame detachment (Coleman, Reference Coleman1998) as a thrust-sense shear zone and a lateral-ramp of the HHD (Fig. 1). However, new data, at a different scale (Fig. 9), allowed us to evaluate kinematic indicators. Both mesoscopic kinematic indicators (e.g. mica fish in the leucogranite and delta-type porphyroclast in orthogneiss in the uppermost part of the GHS, Fig. 9a, b) and asymmetric quartz textures in gneiss (Fig. 9e) suggest that the Chame detachment deformed the orthogneiss and the leucogranite through a non-coaxial flow with a top-down-to-the-N normal-sense of shear. The same kinematics is constrained for the marbles affected by the Phu detachment (Fig. 9e). Therefore, the Chame detachment cannot be directly identified as a strand of the HHD, as far as its entire activity is concerned.
Indeed, considering that the CPOs obtained for both the Chame and the Phu detachments depend on the strain memory of quartz and calcite (i.e. it is not possible to state without doubt that they reflect the entire deformation history of the tectonic discontinuity rather than a part of the plastic history), the Chame detachment can be still interpreted as an early thrust-sense shear zone later reactivated as a normal-sense shear zone, with a full re-orientation of the quartz crystal lattice and a transposition of the kinematic indicators. It is worthy of note that a comparable tectonic history, with a switch in the kinematics, was also proposed for the Annapurna detachment in the nearby Kali Gandaki valley (Vannay & Hodges, Reference Vannay and Hodges1996). In the broad time span of the HHD activity (28–27 to 17 Ma), the deformation was likely localized at different levels complicated by the presence of thrust and lateral ramp or migrated from north to south. Potentially, it is a first thrust-sense shear zone localized (upper HHD) close (or coincident) to the ‘current’ Chame detachment and then migrated some kilometres downward to the level of the JSZ. This is in agreement with the occurrence of higher and older shear zones in the upper part of the GHS such as the Kalopani shear zone in the Kali Gandaki valley active between 41 and 28 Ma (Carosi et al. Reference Carosi, Montomoli, Iaccarino, Massonne, Rubatto, Langone, Gemignani and Visonà2016).
The later reactivation as a normal sense shear zone in the upper GHS (Chame detachment) could be a tectonic event of the orogenic prism in response to the thickening of the GHS caused by the thrust-sense motion of crustal slices of the GHS along the HHD. The occurrence of a normal sense shear zone at the top of the GHS, assisted by erosion, could explain the exhumation of the GHS between the Chame detachment and the HHD.
7.d. Shifting of deformation from the HHD to the MCT
Movements along the JSZ, correlated with the HHD, from 28–27 to 17 Ma caused a discontinuity in the metamorphic paths of the hanging wall and footwall rocks (Fig. 10). Until 28–27 Ma, these rocks experienced prograde burial, but after the activation of the HHD a remarkable difference in the timing of exhumation was introduced (Figs. 8, 10) with the hanging wall exhumed 11–10 My earlier than the footwall. In this sense, P-T paths of rocks across the GHS alone are similar and not sufficient to recognize or to locate the HHD or other tectono-metamorphic discontinuities without studying kinematics and petrochronology.
Figure 10. Sketch of the tectonic evolution on the Greater Himalayan Sequence in the Marsyangdy and Budhi Gandaki river valleys at 28–27 Ma, the time of activation of the Jagat Shear Zone (HHD) (upper) and at 17 Ma, the time of activation of the Main Central Thrust zone (MCTz) (lower). Violet and green dots are reported at depths appropriate for available P-T data.
Sample MSY-03, studied by Gibson et al. (Reference Gibson, Godin, Kellett, Cottle and Archibald2016), is located in the footwall of the JSZ and in the hanging wall of the MCT (Fig. 2a) and constrains, with in situ U-Th-Pb analyses on monazite, prograde metamorphism between 31 and 17 Ma and retrograde metamorphism from 17 to 13 Ma, that is, nearly 10 My later than the cooling and exhumation of the hanging wall of the JSZ. According to Catlos et al. (Reference Catlos, Harrison, Kohn, Grove, Ryerson, Manning and Upreti2001), the mapped metamorphic isograds are not displaced and the monazite ages decrease with no break from the JSZ to the MCT zone (from 22 to 15 Ma). In addition, no discontinuities have been detected in the field from the JSZ to the MCT so, even if the occurrence of minor ‘cryptic’ discontinuities cannot be ruled out, we regard the sample MSY-03 as representative of the behaviour of the footwall of the JSZ as a whole without the occurrence of major discontinuities. The retrograde ages between 17 and 13–10 Ma recorded in the monazite of rocks from the footwall (sample from Gibson et al. Reference Gibson, Godin, Kellett, Cottle and Archibald2016) and hanging wall of the HHD could be caused by the recrystallization/alteration (?) linked to the MCT (Fig. 10). This is in agreement with the ages of the HHD and MCT in Montomoli et al. (Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013), who demonstrated that deformation shifted in space and time from the HHD to the lower MCT at ∼ 17 Ma, and Carosi et al. (Reference Carosi, Montomoli, Rubatto and Visonà2010), who found that HHD shearing ceased at 17 Ma in Western Nepal. Since the position of the MCT is still under debate, spanning for dozens of kilometres, from the chlorite isograd in the LHS up to sillimanite isograd in the GHS along the Himalayan belt (Searle et al. Reference Searle, Law, Godin, Larson, Streule, Cottle and Jessup2008; Martin, Reference Martin2017b; Carosi et al. Reference Carosi, Montomoli and Iaccarino2018) the precise identification and timing of HHD, MCT and shear zones in the upper LHS is necessary to unambiguously locate the MCT in different sections of the belt. Fixing the tectonic boundaries of the GHS and their timing is a fundamental step towards a meaningful identification of the assembly and the exhumation mechanisms of the GHS (Kohn, Reference Kohn2008; Carosi et al. Reference Carosi, Montomoli, Rubatto and Visonà2010, Reference Carosi, Montomoli, Rubatto and Visonà2013; Corrie & Kohn, Reference Corrie and Kohn2011; Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013, Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015; Montemagni et al. Reference Montemagni, Carosi, Fusi, Iaccarino, Montomoli, Villa and Zanchetta2020; Benetti et al. Reference Benetti, Montomoli, Iaccarino, Langone and Carosi2021).
To unambiguously distinguish between the HHD and MCT, it is necessary to integrate structural, metamorphic and geochronological investigations in sections perpendicular to the strike of the belt. The identification and timing of the HHD and MCT constrained in such sections help to better identify the position and age of the MCT, a still debated topic in the Himalayan geology since the last century. In this way, the lower limit of the GHS can be traced.
7.e. Correlation of samples in regional scale shear zones
The Himalayan belt offers the opportunity to directly observe large-scale shear zones both along strike and parallel to the tectonic transport for hundreds of kilometres. Braden et al. (Reference Braden, Godin, Kellet and Yakymchuk2020), who studied two sections of the MCTz from an external klippe to the inner part of the orogenic belt in Western Nepal, found that these sections of the MCTz record different ages of shearing: they are older in the external klippen and younger in the inner portion of the GHS. This led to the interpretation of the occurrence of an out of sequence MCTz in the inner GHS.
Considering the four types of shear zones described by Fossen and Cavalcante (Reference Fossen and Cavalcante2017), except type III that maintains a constant thickness, the active shearing migrates through time towards the inner (core) or towards the external part (wall rocks) of the shear zone depending on the shear zone type. To avoid possible misinterpretations in correlating samples in large-scale regional shear zones, we need to know the geometry, kinematics and evolution type of the investigated shear zone. Therefore, great care should be taken when correlating few samples over long distances because, for example, a younging age from one place to another could be only apparent without the knowledge of the shear zone type.
The framework is complicated by the fact that the same large-scale shear zone (e.g. HHD or equivalents) in a single section shows rejuvenation ages from the upper part to the lower part as documented by detailed petrochronological studies on the Sinuwa and Banuwa thrusts in Central Nepal by Shrestha et al. (Reference Shrestha, Larson, Martin, Guilmette, Smit and Cottle2020). The same happens at a larger scale when considering the position and age of the Kalopani shear zone (Carosi et al. Reference Carosi, Montomoli, Iaccarino, Massonne, Rubatto, Langone, Gemignani and Visonà2016), the HHD (Montomoli et al. Reference Montomoli, Iaccarino, Carosi, Langone and Visonà2013; Carosi et al. Reference Carosi, Montomoli and Iaccarino2018) within the GHS, and the MCT (Montomoli et al. Reference Montomoli, Carosi, Iaccarino, Mukherjee, Carosi, van der Beek, Mukherjee and Robinson2015; Benetti et al. Reference Benetti, Montomoli, Iaccarino, Langone and Carosi2021). So, the problem of a regional and long-lasting shear zone is a 3D problem plus time, for which many data are necessary. In addition to this, Roberts et al. (Reference Roberts, Weinberg, Hunter and Ganade2020) found a syn-orogenic rotation of the MCT in Sikkim, not detected in other sections of the belt. In this case, we should compare geochronological data considering the changing attitude of the stretching lineation.
By considering shear zones occurring in the same vertical section perpendicular to the belt, we can reasonably compare the age of shearing and the kinematics. Much caution is necessary when comparing samples collected along the tectonic transport over long distances, especially in the frontal part of the belt where the GHS is progressively thinned to the point where the THS meets the LHS (Webb et al. Reference Webb, Yin, Harrison, Célérier and Burgess2007). When the GHS is tectonically thinned by components of both pure shear and simple shear during ongoing deformation (Long & Kohn, Reference Long and Kohn2020), shear zones within it could be ‘telescoped’ and the correlation of samples over dozens or hundreds of kilometres can be very difficult so that great caution is needed in their interpretation.
The problem could be overcome by applying a multidisciplinary approach taking into account structural data at the meso- and microscale, P-T estimations and petrochronology, to compare samples collected along different sections. A similar problem is the comparison of shear zones along the strike of arcuate belts. Such belts could be caused by different amounts of displacement along thrusts or affected by inherited lateral discontinuities, striking perpendicular with respect to the main structural trend of the belt and the shear zones (Gibson et al. Reference Gibson, Godin, Kellett, Cottle and Archibald2016; Soucy La Roche & Godin, Reference Soucy La Roche and Godin2019). In this way, two Himalayan sections, along the same parallel, could record different ages of southward propagation of the thrust-sense shear zone.
8. Conclusion
Geological investigations in the Marsyangdi and Budhi Gandaki river valleys allowed us to detect and to map, on a structural basis, a previously unrecognized ductile shear zone (JSZ) in the GHS, separating sillimanite-garnet bearing gneiss from lower kyanite-garnet bearing gneiss. It is a thrust-sense shear zone with a top-to-the-S sense of shear. The JSZ sheared rocks equilibrated at P-T conditions of ∼ 0.7–0.8 GPa and nearly 700°C, whereas footwall rocks equilibrated at higher P of ∼ 1.0–1.1 GPa and a T < 700°C. In situ U-Th-Pb dating on texturally controlled and chemically zoned monazite allowed us to constrain the age of ductile shearing between 28 and 17 Ma. The activity of the JSZ caused the earlier exhumation of its hanging wall joined to the activity of the Chame detachment (a lower strand of the STDS) and/or erosion and tectonic thinning due to the pure shear component of non-coaxial deformation in the GHS. These features confirm that the JSZ is the along-strike continuation of the HHD in Central Nepal and demonstrate the regional extent of the HHD in the Himalayan orogenic belt. To better understand its architecture and evolution, further detailed studies are necessary by integrating structural investigations on sheared rocks, P-T data and detailed petrochronology in sections perpendicular to the grain of the belt.
Detailed structural investigation by meso- and microscale kinematic indicators joined to quartz and calcite CPOs in the Chame detachment and a further strand of the STDS confirmed its extensional, top-to-the-NE kinematics. It firstly acted as a thrust (Walters & Kohn, Reference Walters and Kohn2017) in the upper portion of the Greater Himalan Sequence and then was re-activated as a detachment in response to the motion of the HHD causing an overthickening in the GHS and subsequent exhumation of the hanging wall rocks of the HHD.
The HHD and MCTz are two distict shear zones in the middle and lower GHS which were active at different times within the MCTz marking the lower boundary of the Himalayan metamorphic core. The upper limit of the metamorphic core in the study area is the Chame detachment acting as a normal-sense shear zone (i.e. the local strand of the STDS). Exhumation of the GHS was accomplished by in-sequence shearing with first activation of the HHD (28–17 Ma) and later activation of the MCTz (17−13–10 Ma) at a lower structural position. Precisely determining the boundaries of the metamorphic core, investigating structures within the core and determining the timing of shearing on the boundaries will allow us to better unravel the tectonic evolution and exhumation history of the Himalayan core.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0016756823000365
Acknowledgements
Research funded by PRIN 2015 (University of Torino: R. Carosi and C. Montomoli), funds Ricerca Locale University of Torino (ex-60%, S. Iaccarino and C. Montomoli) and University of California funds (J. Cottle). We are grateful to A. Pêcher and S. Guillot for the fruitful discussions on the geology of the Manaslu massif, to B. Leiss for the support with the XTG, to T. Theye for the support with EMP and to G. Tartaglia for helping with part of the geochronological analysis.
We thank R. D. Law, M.P. Searle and R. Soucy La Roche for their comments and revisions that greatly helped to improve the quality of this paper. We also thank the Editor O. Lacombe for handling this paper.
Author contributions
RC and CM: field work, sampling, formal analysis, methodology conceptualization, writing, project administration and funding acquisition; SI: field work, sampling, formal analysis, methodology, conceptualization, writing and funding acquisition; JC: formal analysis, methodology, writing and funding acquisition; HJM, LN and MS: formal analysis, methodology and writing.
Competing interests
None.
Open access
