Constraining the evolution of shear zones in the Himalayan mid crust in Central–Western Nepal: implications for the tectonic evolution of the Himalayan metamorphic core

Abstract Structural analysis, petrochronology and metamorphic petrology enable identification and bracketing of the timing of a newly mapped high-temperature ductile shear zone (Jagat Shear Zone (JSZ)) in the Himalayan metamorphic core in Central-Western Nepal. In situ U-Th-Pb monazite petrochronology constrains the timing of top-to-the-S/SW shearing between 28–27 Ma and 17 Ma. Burial and prograde metamorphisms in footwall rocks were linked to thrust-sense movement along the JSZ, while the hanging wall rocks were retrogressed and exhumed. The identification and age of the JSZ (as part of a regional system of shear zones: the High Himalayan Discontinuity (HHD)) coupled with the localization and timing of activity of the Main Central Thrust (MCT) (i) fills a gap in tracing the HHD along orogenic strike, (ii) supports the identification of the position and timing of the long-debated MCT and (iii) helps to place the boundaries of the Himalayan metamorphic core and its internal architecture. Thus, our study is a significant step towards a precise identification of the burial, assembly and exhumation mechanisms of the Himalayan metamorphic core.


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 firstorder 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 (1939), its position, age and role in the tectonic evolution of the belt are still debated Martin, 2017aMartin, , 2017bMukhopadhyay et al. 2017Montemagni et al. 2018Montemagni et al. , 2020Chakraborty et al. 2019;Iaccarino et al. 2020). According to Godin et al. (2006) and Catlos et al. (2018Catlos et al. ( , 2020, the age of the MCTz spans between 25 and 18-16 Ma down to nearly 3 Ma (Catlos et al. 2001). 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. 2016a(Parsons et al. , 2016b 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 Montemagni et al. 2019Montemagni et al. , 2020. Montomoli et al. (2015) pointed to the occurrence of a regionalscale structure (High Himalayan Discontinuity (HHD)) within the GHS. Carosi et al. (2018) 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. 2006). 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, 1998;Searle & Godin, 2003;Searle, 2010) (Fig. 2a), has been interpreted as a ramp of the HHD (Walters & Kohn, 2017) 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. 2000;Williams & Jercinovic, 2002, 2012Carosi et al. 2022;Dumond et al. 2022) and in situ radiometric dating of mineral (re)crystallization within different fabrics helps to unravel the deformation history (Di Vincenzo et al. 2004;Villa, 2016). 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. 2000;Rubatto et al. 2001Rubatto et al. , 2013Gibson et al. 2004;Kohn et al. 2005;Pyle et al. 2005;Pyle, 2006;Corrie & Kohn, 2011;Kohn, 2016). 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, 2017;Schulz, 2021;Dumond et al. 2022;Imayama et al. 2022). 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. 1999Williams et al. , 2007Gibson et al. 2004;Pyle, 2006;Rubatto et al. 2013;Kohn et al. 2017). In this paper, we focus on laser ablation split stream (LASS) dating (Cottle et al. 2012;Kylander-Clark et al. 2013) 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).
The HHD was first recognized in Western Nepal (Carosi et al. 2010;Montomoli et al. 2013Montomoli et al. , 2015; see Carosi et al. 2018, 2019 for a review), but gaps exist remain in our understanding of its alongstrike continuation to the east Waters, 2019). 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.a. The Himalayas
The Himalayas (Fig. 1) developed after the initial collision between the Indian and Eurasian plates at~59-55 Ma (Le Fort 1975;Green et al. 2008;Hu et al. 2016). 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. 2012;Searle, 2019). 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, 1999;Hodges, 2000). A km-thick ductile to ductile-brittle shear zone, the MCTz (Heim & Gansser, 1939;Gansser, 1964;Searle et al. 2008;Martin, 2017aMartin, , 2017b, separates the mostly greenschistwith minor lower amphibolitefacies 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, 1975;Hodges, 2000;Carosi et al. 2018 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, 2010(Searle, , 2013Visonà et al. 2012, Weinberg, 2016, such as the Manaslu granite (Guillot et al. 1995;Cottle et al. 2019), are abundant in the upper part of the GHS. Tectonically above the GHS, unmetamorphosed sedimentary rocks (Frank & Fuchs, 1970;Garzanti, 1999) with minor amphibolite facies rocks    Montomoli et al. 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. 1984;Burchfiel et al. 1992;Carosi et al. 1998;, Kellett et al. 2010Searle, 2010;Iaccarino et al. 2017b;Nania et al. 2022). The occurrence of the STDS at the top of the GHS and of the MCTz at its bottom, as two coeval northdipping 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. 2010;Martin et al. 2010Martin et al. , 2015Corrie & Kohn, 2011;Imayama et al. 2012Imayama et al. , 2022Montomoli et al. 2013Montomoli et al. , 2015Ambrose et al. 2015;Cottle et al. 2015;He et al. 2015;Khanal et al. 2015;Larson et al. 2015;Wang et al. 2015Wang et al. , 2016Zeiger et al. 2015;Agustsson et al. 2016;Iaccarino et al. 2017a;Walters & Kohn, 2017;Goscombe et al. 2018;Chakraborty et al. 2019;Waters, 2019; 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.

1264
R Carosi et al. Shrestha et al. 2020;Benetti et al. 2021) 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 . 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. 2010;Montomoli et al. 2013Montomoli et al. , 2015.

2.b. Localization of the MCT
Mapping, localization, timing and tectono-metamorphic evolution of the MCTz are still debated (Martin, 2017a). 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, 1939) 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, 1975;Pêcher, 1989); (iii) a protolith boundary, where the MCTz separates rocks with different Nd and Sr isotopic signatures (e.g. Ahmad et al. 2000) and with different maximum depositional age based on zircon U-Pb populations (Parrish & Hodges, 1996) and (iv) a structural criteria , 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. 2013;Parsons et al. 2016aParsons et al. , 2016bParsons et al. , 2016cLarson et al. 2010Larson et al. , 2017Hunter et al. 2019). Frequently, a contractional top-to-the-S high-strain zone corresponds to the base of an inverted metamorphic sequence Iaccarino et al. 2020). Moreover, several authors mapped two different thrusts, MCT-I and MCT-II (Arita, 1983) or Munsiari and Vaikrita thrusts (Valdiya et al. 1999) 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. 1986;Larson et al. 2010Larson et al. , 2011Martin et al. 2015;Parsons et al. 2016c).

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. 1986;Pêcher, 1989Pêcher, , 1991Coleman, 1996Coleman, , 1998Coleman & Hodges, 1998;Parsons et al. 2016c). 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, 1989(Pêcher, , 1991Catlos et al. 2001Catlos et al. , 2018Larson et al., 2010;Catlos, 2021;Martin et al. 2015;Gibson et al. 2016) defining a 2 to 4 km-thick high-strain zone with top-to-the-S kinematics. According to Pêcher (1989) and Catlos (2021), the top of the MCTz is localized just above the kyanite-in isograd. Recent garnetbased pressure-temperature (P-T) paths were presented by Catlos et al. (2018), following Catlos et al. (2001), for rocks across the MCTz suggesting an increase of both pressure and temperature traced structurally upwards. Coleman (1998) and Catlos et al. (2001) 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, 2016). Recent monazite dating by Gibson et al. (2016) (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, 1975;Colchen et al. 1986;, now more properly referred as units (Searle & Godin, 2003). 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 Coleman (1996), 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 diopsideamphibole 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, 1996;Parsons et al. 2016aParsons et al. , 2016bParsons et al. , 2016cFig. 2a), and the upper Phu detachment (Fig. 2a) (Searle & Godin, 2003;Searle, 2010). Low-grade marble and metapelite (Schneider & Marsh, 1993) 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. 1994;Cottle et al. 2019), 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, 1996, Coleman & Hodges, 1998Searle, 2010, Parsons et al. 2016a, 2016b, 2016c. In contrast, Walters and Kohn (2017) 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. 2010(Larson et al. , 2011 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. 1986) and the occurrence of carbonate-bearing lithologies within the MCTz (Colchen et al. 1986;Nania, 2021). A ductile shear zone within Unit I, 1 km south of the village of Philim (Fig. 2b), crops out (Pêcher & Fort, 2012).

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).
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 , 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 hightemperature grain boundary migration (GBM II ) and recrystallization leading to chessboard extinction patterns ( Fig. 4f; Kruhl, 1996;Law, 2014). 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, 1996;Passchier & Trouw, 2005) in agreement with the synkinematic mineral assemblage.
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, 2005) 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 (GBM II ) (Stipp et al. 2002) 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. (2016c) highlighted on their map S foliations and C planes in correspondence of the shear zones. Pêcher and Fort (2012) hypothesized a thrust-sense shear zone by the occurrence of high-strain levels containing sheath folds.
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.

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 (2012). 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, 1994, 2008version 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, 2000). Following Yakymchuk and Godin (2012), 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. (2007) 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. 2007). 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.
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.

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) 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-biotitegarnet 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, 1999Spear & Pyle, 2010).
All monazite grains are mainly concentrically zoned (Catlos, 2013), 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  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.
There is a good correlation between 208 Pb/ 232 Th 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 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. Due to the good correlation among dates and monazite chemical composition, we interpret the obtained dates as follows according to Kohn et al. (2005), Braden et al. (2017) and Soucy La Roche et al. (2018). 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, 1999Gibson et al. 2004;Kohn et al. 2005;Spear & Pyle, 2010 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. (2004) 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. (2022) 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 Gibson et al. 2004;Spear & Pyle, 2010). 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. 2005).

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 cm 2 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, 1997) with a halfwidth of 10° ( Hielscher & Schaeben, 2008).  Table 1 and main text for details, Fig. 2 for sample location). KFMASH equilibria plotted with the GIBBS software (Spear & Menard 1989) using the SPaC2007 database. NKFMASH curves are after Spear et al. (1999). (b) Reconstructed P-T paths of different GHS portions, based on geothermobarometry and petrographic observations.

R Carosi et al.
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  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, 1982) and a M-Index of 0.01 (designed for large ODF datasets, Skemer et al. 2005). 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, 2014). The asymmetric pattern suggests a non-coaxial flow (Lister & Hobbs, 1980;Wallis, 1992;Law & Johnson, 2010). Dextral asymmetry of [c]-axis maxima, compared to the geographical reference system, is coherent with a top-to-the-N non-coaxial deformation.
According to the c-axis-based method of Wenk et al. (1987), calcite asymmetric texture (of ca. 7°) results from a general flow with low simple shear contribution (~26%).

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 (1998), Catlos et al. (2001, Gibson et al. (2016) and Larson et al. (2011) 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. (2013Montomoli et al. ( , 2015 and Iaccarino   Iaccarino et al. (2015) and Wang et al. (2015Wang et al. ( , 2016 in Central Nepal, Chakraborty et al. (2019) and Benetti et al. (2021) in India and Waters (2019) 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 (1998) and Catlos et al. (2001). 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, 2017). 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 . 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. (2013Montomoli et al. ( , 2015 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. 2013;Iaccarino et al. 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. (2020) using detailed petrochronology on rocks from both the Banuwa and Sinuwa thrusts were able to document the occurrence of two main insequence 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. (2010Martin et al. ( , 2015 proposed normal-sense kinematics, whereas Corrie and Kohn (2011) suggested a thrust-sense kinematics active from 27 to 22 Ma. Shrestha et al. (2020) 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. (2015) 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 Larson et al. 2015) where one or two tectono-metamorphic discontinuities have been documented so far (Carosi et al. 2016 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 sillimanitebearing migmatites, with a higher degree of melting, and kyanitebearing gneiss below .
Further investigations are required to better understand if the HHD is a unique shear zone or if it is made of several (insequence?) shear zones and what type of shear zone the HHD is (Type I, II, III and IV sensu Fossen & Cavalcante, 2017). 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 highgrade metamorphic rocks. As suggested by Montomoli et al. (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 (2017) identified a strand of the HHD in the upper part of the Marsyangdi valley interpreting the Chame detachment (Coleman, 1998) 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, 1996). 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 1274 R Carosi et al.
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. 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.
Sample MSY-03, studied by Gibson et al. (2016), 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. (2001), 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. 2016) 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. (2013), who demonstrated that deformation shifted in space and time from the HHD to the lower MCT at~17 Ma, and Carosi et al. (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 Martin, 2017b;Carosi et al. 2018) 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, 2008;Carosi et al. 2010Carosi et al. , 2013Corrie & Kohn, 2011;Montomoli et al. 2013Montomoli et al. , 2015Montemagni et al. 2020;Benetti et al. 2021).
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. (2020), 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 (2017), 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 largescale 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. (2020). The same happens at a larger scale when considering the position and age of the Kalopani shear zone (Carosi et al. 2016), the HHD Carosi et al. 2018) within the GHS, and the MCT Benetti et al. 2021). 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. (2020) 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. 2007). When the GHS is tectonically thinned by components of both pure shear and simple shear during ongoing deformation (Long & Kohn, 2020), 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. 2016;Soucy La Roche & Godin, 2019). In this way, two Himalayan sections, along the same parallel, could record different ages of southward propagation of the thrust-sense shear zone.

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 kyanitegarnet bearing gneiss. It is a thrust-sense shear zone with a top-tothe-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, 2017) 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