Thermal history of Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex, NE India: New insights into pre-Cretaceous subduction channel tectonics within the Neo-Tethys

Thermal history of Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex, NE India: New insights into pre-Cretaceous subduction channel tectonics within the Neo-Tethys

Accepted Manuscript Thermal history of Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex, NE India: New insights into pre-...

1MB Sizes 0 Downloads 17 Views

Accepted Manuscript Thermal history of Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex, NE India: New insights into pre-Cretaceous subduction channel tectonics within the NeoTethys

M. Rajkakati, S.K. Bhowmik, A. Ao, T.R. Ireland, J. Avila, G.L. Clarke, A. Bhandari, J.C. Aitchison PII: DOI: Article Number: Reference:

S0024-4937(19)30323-8 https://doi.org/10.1016/j.lithos.2019.105166 105166 LITHOS 105166

To appear in:

LITHOS

Received date: Revised date: Accepted date:

31 January 2019 4 August 2019 5 August 2019

Please cite this article as: M. Rajkakati, S.K. Bhowmik, A. Ao, et al., Thermal history of Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex, NE India: New insights into pre-Cretaceous subduction channel tectonics within the NeoTethys, LITHOS, https://doi.org/10.1016/j.lithos.2019.105166

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT Thermal history of Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex, NE India: New insights into pre-Cretaceous subduction channel tectonics within the NeoTethys

SC RI PT

M. Rajkakati1, S.K.Bhowmik1,*, A. Ao2, T. R. Ireland3, J. Avila3, G.L. Clarke4, A. Bhandari5, J.C. Aitchison6 1

Department of Geology & Geophysics, Indian Institute of Technology, Kharagpur 721 302,

NU

India

Wadia Institute of Himalayan Geology, Dehradun, Uttarakhand 248001, India

3

Research School of Earth Sciences, Australian National University, Canberra ACT 0200,

MA

2

Australia

School of Geosciences, University of Sydney, Sydney NSW 2006, Australia

5

Geological Survey of India, Faridabad, Haryana 121001, India

6

School of Earth and Environmental Sciences, The University of Queensland, St Lucia,

PT

ED

4

AC

ABSTRACT

CE

Queensland 4072, Australia

Eclogite from two locations in a fossil accretionary complex now exposed in Nagaland, NE India, at the northern end of the Indo-Myanmar ranges, provides the oldest evidence for Neo-Tethyan subduction along the Yarlung-Tsangpo suture zone. Metamorphic texture, mineral composition, peak metamorphic P-T estimates, reconstructed metamorphic P-T paths, and U-Pb zircon ages from three eclogite samples collected near Mokie and Thewati villages collectively reveal snapshots of a c. 205–172 Ma subduction burialexhumation cycle in the eastern Neo-Tethys. The Thewati eclogite records a clockwise

ACCEPTED MANUSCRIPT metamorphic P-T path involving: (1) prograde metamorphism traversing the epidote blueschist facies and culminated in peak eclogite facies conditions at 25-28 kbar and ~650 o

C; (2) an early retrograde stage involving decompression with cooling in the eclogite facies

at ~18.3 kbar and 630 oC; followed by (3) cooling through the epidote blueschist, transitional lawsonite blueschist and greenschist facies to ~6 kbar and 300 oC. Circa 189-185 Ma peak

SC RI PT

metamorphism is inferred for the Mokie and Thewati eclogites, and was accompanied by an apparently low thermal gradient (~7-8 oC/km). Such conditions and the clockwise P-T path are associated with a tectonic scenario involving cold mature subduction within the NeoTethys in Nagaland. This

is thus the first comprehensive dataset of an Early Jurassic

NU

subduction channel (with respect to present geographic co-ordinates) for the eastern margin of the Neo-Tethyan suture. The Nagaland samples offer a critical record of a pre-Cretaceous

MA

subduction system within the Neo-Tethys, which incorporated an Early-Middle Jurassic Andean-type convergent margin extending along the southern Eurasian margin from Pakistan

ED

to the Indo-Myanmar Ranges.

PT

Key words: Eclogite; U-Pb zircon geochronology; Metamorphism; CW P-T Path;

AC

1. Introduction

CE

Subduction; Neo-Tethys; Yarlung-Tsangpo Suture Zone

Much of our understanding of Mesozoic plate tectonics, including the sequential amalgamation of Gondwanan continental fragments that assembled the Eurasian continent, derives from metamorphic and isotopic studies of three Tethyan suture zones that dissect Tibet. These are: (1) the Longmu Co-Shuanghu (LCS); (2) Bangong-Nujiang (BNS); and the (3) Yarlung-Tsangpo (YTS) sutures (Fig.1). They include blueschist (Anczkiewicz et al., 2000; Groppo et al., 2016), eclogite (Zhang, Y-X et al., 2018) and amphibolite facies rocks

ACCEPTED MANUSCRIPT (Guilmette et al., 2012; Wang et al 2008). These rocks provide information critical to understanding tectonic processes linked to the closure of the Palaeo- (Devonian-Triassic), Meso- (Early Permian to Late Cretaceous) and Neo-Tethyan (Late Triassic-Paleogene) oceans (see Zhu et al., 2013 for a review). The Yarlung-Tsangpo suture zone (YTSZ) is the youngest and southernmost of the

SC RI PT

three sutures that dissect Tibet. It is potentially the most cryptic in terms of understanding the timing, duration, and periodicity of the Neo-Tethyan subduction (see Hébert et al., 2012 for a review). Tectonic models incorporating the northward subduction of the Neo-Tethyan oceanic lithosphere beneath southern Tibet (Lhasa continental block) are well established,

NU

covering much of the Cretaceous and extending to the Early Cenozoic collision of India and Asia (e.g. Yin and Harrison, 2000). However, the pre-Cretaceous geodynamic history of the

MA

Neo-Tethyan lithosphere remains uncertain and subject to on-going discussion. Models include the opening of the Neo-Tethyan ocean as a Late Triassic back-arc basin due to the

ED

separation of the Lhasa block from northern Australia (Zhu et al., 2013), and the initiation of its subduction at c.180 Ma (Hébert et al., 2012). The existence of Late Triassic to Early

PT

Jurassic Neo-Tethyan subduction beneath the Lhasa Terrane has also been inferred on the

CE

basis of the age and geochemistry of arc components of the Gangdese magmatic belt (Yu et al., 2018). However, in the absence of direct evidence of related HP/LT metamorphism, the

AC

subduction channel dynamics of such early events within the Neo-Tethyan sutures are not well understood.

The Nagaland Ophiolite Complex (NOC), at the northern extent of the Indo-Myanmar Ranges (IMR) (Fig. 2a), is considered a southeastern extension of the Yarlung-Tsangpo suture. It exposes a rare Neo-Tethyan metamorphic sequence including pumpellyite-diopside, greenschist, blueschist and eclogite facies rocks (Chatterjee and Ghose, 2010; Ao and Bhowmik, 2014; Bhowmik and Ao, 2016) (Fig. 2b). Preliminary chronological constraints

ACCEPTED MANUSCRIPT indicate Middle to Upper Jurassic protolith ages for elements of the metamorphosed oceanic crust (Sarkar et al., 1996; Baxter et al., 2011; Aitchison et al., 2019). The timing and thermal history of eclogite facies metamorphism and its relationship with other HP metamorphites in the terrane are not known (see also Ao and Bhowmik, 2014; Bhowmik and Ao, 2016). Numerous studies indicate that eclogite facies rocks preserve critical evidence for the timing

SC RI PT

and thermal history of oceanic subduction (Agard et al., 2009 and references cited therein). In particular, low to medium temperature oceanic eclogites commonly preserve excellent textural and mineral compositional (e.g. compositional zoning in garnet, omphacite, epidote, amphibole and phengite) features amenable to the near complete reconstructions of

NU

metamorphic P-T paths, including well-constrained prograde histories (e.g. Franciscan highgrade eclogite block, California, Page et al., 2007; lawsonite eclogite pods, Sivrihisar Massif,

MA

Turkey, Davis and Whitney, 2008; lawsonite eclogite, Alpine Corsica, Brovarone et al., 2011; Zermatt-Saas, Western Alps, Reinecke, 1998, Angiboust et al., 2009; Monviso, Alps,

ED

Groppo & Castelli, 2010; Angiboust et al., 2012; Schistes-Lustres, Alps, Agard et al., 2001; eclogites of Pouebo terrane, Pam Peninsula, NE New Caledonia, Clarke et al., 1997;

PT

Dominican Republic, Krebs et al., 2011; Rio San Juan mélange Complex, Krebs et al., 2011).

CE

Moreover, recent analytical and theoretical advances that combine in situ U-Pb dating and trace element compositional analysis of zircon offer unprecedented insights into the

AC

timescales involved in forming high-pressure/ultra-high-pressure (HP/UHP) metamorphic terranes (Rubatto and Hermann, 2007; Rubatto, 2017). The Neo-Tethyan tectonic models for Southeast Asia, in general, are strongly influenced by data from the Tibetan belts with a lack of detail from the IMR and adjoining ophiolite belts in eastern and north-eastern Myanmar. Seen in this context, the Nagaland Ophiolite Complex presents a unique opportunity to: (a) test existing predictions for the existence of a pre-Cretaceous subduction system within the Neo-Tethys (e.g. Zedong terrane

ACCEPTED MANUSCRIPT - McDermid et al. 2002; Aitchison et al. 2007) and (b) determine the thermal and dynamic history of such a subduction zone using well-constrained metamorphic P-T paths from the eclogites. In this investigation, we establish the thermal history of an Early Jurassic NeoTethyan subduction system, considered to be the oldest within the Yarlung-Tsangpo suture zone, using textural and mineral compositional data, metamorphic reaction history, peak

SC RI PT

metamorphic P-T estimates, reconstructed metamorphic P-T path and SHRIMP U-Pb zircon age data from two new eclogites locations in the NOC. Given the importance of these findings in the tectonic frame work of Southeast Asia, we also review the published literature from the ophiolite belts of IMR and eastern and northeastern Myanmar. These new results of

NU

Neo-Tethyan tectonics from Nagaland and adjoining Myanmar are correlated with those from

MA

the YTSZ in the Tibetan sector.

2. Geological setting of the Nagaland eclogites

ED

Eclogite samples investigated in this work come from the Nagaland Ophiolite Complex at the northern end of the Indo Myanmar Ranges. The IMR are located at the junction between two

PT

orogenic systems: the Alpine-Himalayan continental collision belt in the north, and the active

CE

subduction zone of the Sunda–Andaman arc to the south (Fig. 2a). These orogenic systems resulted from the interaction of three tectonic plates: the Indian Plate in the southwest, the

AC

West Myanmar Block (alternatively Myanmar Microplate) in the east, and the Asian Plate to the north and northeast (Fig. 2a) (Searle et al., 2007). The IMR, along the western margin of the West Myanmar Block, is considered an accretionary wedge that developed in response to eastward subduction of the Neo-Tethys (attached to the Indian plate) beneath the West Myanmar Block (Searle et al. 2007; Acharyya, 2015; Aitchison et al., 2019). Ophiolitic rocks, with or without HP metamorphic rocks, in the IMR and the West Myanmar Block occur in three sub-parallel, N-trending ophiolite belts, called the Western

ACCEPTED MANUSCRIPT (WOB), Central (COB) and Eastern Ophiolite Belts (EOB) (terminology after Htay et al. 2017) (Fig. 2a). The WOB is the longest of the three belts and separates the eastern margin of the IMR from the West Myanmar Block. It can be traced for ~3000 km from Mayodia (Arunachal Pradesh) in the north near the eastern Himalayan syntaxis through the Naga and Manipur Hills, Kalaymyo, Chin Hills, Mt. Victoria-Kawlun belt in the middle section to the

SC RI PT

Andaman-Nicobar Islands in the south (Fig. 2a). Ophiolitic sequences in this belt generally occur as small, rootless and tectonised bodies with shallow-dipping foliations in low-grade metamorphic rocks, which are overlain by Eocene–Oligocene sedimentary rocks of flysch affinity (Acharyya, 2007). A Late Jurassic age for ultramafic rocks of the WOB was posited

NU

from a K–Ar age for hornblende-bearing pegmatite (158 ± 20 Ma) that intruded serpentinites in the Mt Victoria area (Mitchell, 1986). However, recent U-Pb dating of zircon indicates

MA

Early Cretaceous protoliths for ophiolitic components in this belt (e.g. c.116-119 Ma Manipur ophiolite, Singh et al, 2017, Aitchison et al., 2019; 127 Ma Kalaymyo ophiolite, Liu et al.,

ED

2016a). The formation of amphibolite-facies rocks in Kalaymyo has been interpreted as a metamorphic sole to the ophiolite and is dated at 115-119 Ma (Zhang et al., 2017). In

PT

summary, the ophiolites of the WOB are thought to mark a suture zone, where the Neo-

CE

Tethyan Ocean once existed (Searle et al., 2007; Liu et al., 2016a). The Central Ophiolite Belt extends from the Putao area of Kachin State in the north

AC

through the Jade Mines Belt (JMB) to the Sagaing-Minwun ranges in the south (Fig. 2a). It is composed of a dismembered, incompletely preserved ophiolitic sequence of serpentinite and peridotite (Htay et al., 2017). The Jade Mines Belt (JMB) includes jadeitite, eclogite, amphibolite, blueschist, chromitite, and serpentinized peridotite (Goffé et al., 2002) that reflect the effects of high-P/low-T metamorphism. The P–T conditions have been variously estimated at 10–15 kbar and 300–500 °C (Mével & Kiénast, 1986), >14 kbar and 400-450 °C (Goffé et al., 2002), >10 kbar and 250–370 °C (Shi et al., 2003) and ~15 kbar and ~380 °C

ACCEPTED MANUSCRIPT (Oberhänsli et al., 2007). U–Pb Jurassic ages of zircon taken from jadeitite (163 – 160 Ma) have been interpreted as protolith ages for the ophiolite in the JMB (Shi et al., 2008; Yui et al., 2013). The timing of the HP metamorphic event is not well defined, with age estimates ranging from the Late Jurassic (U–Pb zircon age between 158 and 147 Ma, Shi et al., 2008; Qiu et al., 2009) to the Late Cretaceous (U–Pb zircon age of 77 ± 3 Ma, Yui et al., 2013 and 40

Ar/39Ar age of 80 Ma for phengite in eclogite, blueschist, jadeitite, and amphibolite,

SC RI PT

an

Goffé et al., 2002). Zhang, J. et al. (2018) recently reinterpreted all published zircon age data from the JMB and concluded that the main HP metamorphism and accretionary orogenesis in this segment of the COB was of Early to Middle Jurassic age.

NU

The Eastern Ophiolite Belt (EOB) is represented by the Tagaung-Myitkyina belt (Fig. 2a) comprising massive peridotite of harzburgite to lherzolite composition, dacite and

MA

andesite, hornblende gabbro, diorite, granodiorite and veins of gabbro, rodingite and plagiogranite in peridotite (Liu et al., 2016b; Zhang, J. et al., 2018). The volcanic and

ED

gabbroic to dioritic rocks have supra-subduction-zone-type chemical signatures (Yang et al., 2012; Xu et al., 2017). Recent U-Pb dating of zircon from a variety of lithologies in the

PT

Tagaung-Myitkyina belt has yielded a Middle Jurassic protolith age (ca. 177-166 Ma, Yang

CE

et al., 2012, Liu et al., 2016a; Gardiner et al., 2018). The EOB has been interpreted as: (a) the southern continuation of the Meso-Tethyan suture in central Tibet (Liu et al., 2016a); (b) a

AC

klippen of the Medial Myanmar Suture Zone that developed along the eastern margin of the Mogok Metamorphic and Slate Belt (Mitchell et al., 2015); and (c) a Jurassic continental margin arc (Zhang, J. et al., 2018). The 200 km long and 5-15 km wide, NNE-trending Nagaland Ophiolite Complex can be distinguished from other ophiolite bodies in the WOB by: (1) the presence of a variety of HP/LT metamorphic rocks in ophiolitic mélange; and (2) its association with a disrupted oceanic plate sequence of chert, pillowed basalt, plagiogranite, mafic, and ultramafic

ACCEPTED MANUSCRIPT cumulate and serpentinite (Fig. 2b) (Acharyya, 1986). Some of its serpentinite has mantle peridotite lineage (Brunnschweiler, 1966). Ophiolitic mélange structurally overlies the Disang Formation (Fig. 2b), which has Indian continent-derived clastic sedimentary rock of Eocene age (Imchen et al., 2014). The ophiolite mélange is in turn structurally overlain by an Early Paleozoic metasedimentary sequence (Aitchison et al., 2019) forming the Naga

SC RI PT

Metamorphics (Brunnschweiler, 1966, Fig. 2b). The NOC includes accretionary metasedimentary rocks, dominantly of pelagic chert-limestone association, with attenuated metabasalt, metagabbro, serpentinite, and pillow basalt with igneous texture and mineralogy (Acharyya, 1986; Ao and Bhowmik, 2014; Bhowmik and Ao, 2016). Metamorphic

NU

assemblages in the NOC reflect greenschist, pumpellyite-diopside, and blueschist to eclogite facies conditions (Fig. 2b) (Ao and Bhowmik, 2014; Bhowmik and Ao, 2016; Chatterjee and

MA

Ghose, 2010). Ao and Bhowmik (2014) previously established evidence for a disrupted field array in the Mokie-Satuza sequence in the western part of the mapped area, with

ED

metamorphic grade increasing eastward from greenschist through pumpellyite-diopside to lawsonite blueschist facies (Fig. 2b).

PT

Preliminary whole rock K-Ar geochronology of metabasalt (148 ± 4 Ma, Sarkar et al.,

CE

1996) and radiolarian microfossils from elements of the NOC (latest Bathonian to early Callovian, Aitchison et al., 2019) indicate Middle-Upper Jurassic protoliths that are

AC

interpreted as having been scraped off the subducting oceanic lithosphere. Chatterjee and Ghose (2010) estimated peak metamorphic conditions for eclogite in the NOC at 17–20 kbar and 580–610°C. Bhowmik and Ao (2016) recently reported a hornblende eclogite facies metamorphic sole (peak metamorphism at 13.8 kbar, 625 oC) from a nearby area. In this work, we describe two new eclogite locations: (1) near Mokie village (sample N22F, Mokie eclogite from now on), and (2) ~2 km east of Old Thewati village (samples N23B, N40B, hereafter Thewati eclogite) (Fig. 2). In both locations, eclogite occurs as

ACCEPTED MANUSCRIPT small, meter-scale lensoidal bodies within a lawsonite blueschist facies sequence of basalt, chert and limestone (Fig. 3a-b). Thin discontinuous bands rich in tremolite and/or white mica separate the eclogite bodies from the host metabasaltic rocks (Fig. 3a-b).

2.1 Sample description

SC RI PT

Sample N23B is a coarse-grained eclogite, collected from the core region of the Thewati eclogite body (Fig. 3b). It consists of porphyroblastic garnet (5-7 mm grain diameter) (~18 modal%) and omphacite (2-5 mm grain diameter) (~19%) enveloped by a matrix of phengite (~31%), barroisite (Fig. 3c) and glaucophane (collectively ~22%), and chlorite (~8%), with

NU

minor epidote (~3%), rutile, sphene and accessory zircon (Fig. 3a, 3b, 3d). Sample N40B is a composite, omphacite-rich (~37-57%) eclogite (Fig. 3e) collected from the outer margin of

MA

the Thewati eclogite body (Fig. 3b). Pale green layers dominated by foliated omphacite (Fig. 3f) alternate with darker barroisite-rich layers. Porphyroblastic garnet (1-5 mm in diameter),

ED

enveloped by finely foliated omphacite (Fig. 3e, f) occurs throughout the sample (Fig. 3e). Thin discontinuous, layer-parallel segregations of epidote are also present in the rock. Sample

PT

N22F is a fine-grained, schistose rock collected from the outer margin of the Mokie eclogite

CE

body. It consists of garnet and Na-Ca pyroxene of omphacite to aegirine-augite composition, enveloped by Na-Ca amphibole of barroisitic to magnesio-kataphoritic composition, sphene,

AC

minor glaucophane, phengite, epidote and rutile and accessory zircon (Fig. 3g).The foliation is cut by phengite-rich veins. A total of five thin sections (two from N23B, here named N23B-1 and N23B-3 and three from N40B, named N40B-2, N40B-4 and N40B-6) was utilised to recover the metamorphic history of the Thewati eclogite. Sample N23B-1 was used for the geochronology. The Mokie eclogite sample N22F was used for both geochronological and metamorphic studies.

ACCEPTED MANUSCRIPT 3.Methods 3.1. Electron probe micro-analytical technique Major-element mineral chemistry was determined using a CAMECA SX-100 electron microprobe fitted with four spectrometers at the Department of Geology and Geophysics, Indian Institute of Technology (IIT), Kharagpur. Analyses were performed with an

SC RI PT

accelerating voltage of 15 kV, a beam current of 15 nA and a beam size of ~1µm (~5 μm for phengite). Natural mineral standards were used for calibration, including orthoclase (for K and Si), jadeite (Na), wollastonite (Ca), periclase (Mg), hematite (Fe), corundum (Al), rhodonite (Mn) and chromite (Cr). Analytical variation within standards was kept at or below

NU

∼1.5wt% for primary standards. Representative mineral chemical analyses (presented in

3.2. Metamorphic Reconstruction

MA

Tables 1-3) were recalculated using the AX program (THERMOCALC software, v. 3.26).

ED

3.2.1. Methods of determination of metamorphic pressures (P) and temperatures (T) We adopt a conventional thermobarometric approach to recover the P–T conditions of

PT

eclogite facies metamorphism at the two locations. This utilizes simultaneous solutions of

CE

garnet-clinopyroxene Fe2+-Mg exchange thermometric (reaction R1, 2 pyrope + 3 ferrosilite = 2 almandine + 3 enstatite, after Ravna, 2000) and garnet-clinopyroxene-phengite

AC

barometric (reaction R2, 6 diopside + 3 muscovite = 2 grossular + 1 pyrope + 3 celadonite, after Ravna and Terry, 2004) equations. This methodology has been widely applied to eclogite localities from around the world. P–T estimates made using these calibrations (from now on R00_RT04 and Method 1) are also compared with the results of computations of the same equilibria using the internally consistent thermodynamic dataset of Holland and Powell (1998) and THERMOCALC software (version 3.31, thermodynamic database tc-ds55, 22 November 2003) (Method 2). Mineral end-member activities and their uncertainties were

ACCEPTED MANUSCRIPT calculated using activity-composition (a-X) datafile coding for pseudosection calculations, following the approach outlined by Powell and Holland (2008). Solid solution models used for this purpose are clinopyroxene (Green et al., 2007), garnet (White et al., 2007) and phengite (Coggon and Holland, 2002). For the spessartine-rich garnet composition (Sps17) in N22F, the garnet activity model of Ganguly et al. (1996) was adopted for the Method 2. The

SC RI PT

model incorporates Mg-Mn interaction relevant to the peak mineral assemblage in the sample and for the computations of garnet end-member activities. Finally, we have also estimated peak conditions using the THERMOCALC average P-T method (Method 3) (see Bhowmik and Ao, 2016 for additional information on the methodology). The average P-T computations

NU

were carried out using the same garnet-clinopyroxene-phengite assemblage and mineral chemistry. The same THERMOCALC software, thermodynamic database and mineral end

MA

member activity models as adopted in the Method 2 has been applied here. For all the three methodologies, we have estimated ferric content in clinopyroxene using charge-balance

ED

criterion and in garnet following the methodology of Ravna (2000). The P-T results are

PT

presented in Table 4.

CE

3.2.2 P-T Pseudosection modelling

Estimates of peak metamorphic conditions and the likely reaction history of the Nagaland

AC

eclogite were also evaluated using calculated metamorphic phase diagrams. Sample N40B preserves excellent reaction textures (Fig. 4) and mineral compositional features (Fig. 5-6, Tables 1-3) consistent with evidence for prograde, peak and retrograde conditions being preserved in selected domains of the rock (Section 4). For this reason, we calculated P-T pseudosections using effective bulk rock compositions for two domains in the N40B: domain N40B-4 represents the garnet-pyroxene-rich peak assemblage; and domain N40B-6 is rich in retrograde minerals such as barroisite and epidote and represents the retrograde stage.

ACCEPTED MANUSCRIPT Effective bulk rock compositions were calculated using XMapTools software, following the methodology of Lanari (2014). For each domain, a mineral mask file was prepared representing the modes of constituent minerals (Fig. 7a-b, Table 5). Mineral modes were then converted to their mass fractions and then combined with the calculated average oxide weight percentage of all major elements of minerals present in the chosen domain to obtain effective

SC RI PT

rock compositions (Table 5).

The chemical system Na2O - CaO - K2O - FeO - MgO - Al2O3 - SiO2 - H2O- TiO2 - O2 (NCKFMASHTO) was used for the computation of the pseudosections, together with the same version of the THEMOCALC software and thermodynamic dataset used in the

NU

thermobarometry. The component MnO was omitted because all minerals other than garnet in the mafic eclogites have very low manganese content, and there are currently no manganese

MA

solid-solution models for several of the rock-forming minerals. The inclusion of manganese in the model system has the main effect of subtly expanding garnet-bearing equilibria to

ED

lower-grade conditions. However, as the spessartine content of garnet grain cores can be comparatively high (Sps11) the proportion of MnO was combined with the FeO content in the

PT

rock. The following mixing models were applied: clinopyroxene (Green et al., 2007),

CE

amphibole (Diener et al., 2007), garnet (White et al., 2007), epidote (Holland and Powell, 1998), plagioclase (Holland and Powell, 2003) and phengite (Coggon and Holland, 2002).

AC

Quartz, sphene, rutile and lawsonite are treated as pure phases. Fe3+ contents in the effective rock compositions have estimated on the basis of average Fe3+ contents in the constituent minerals and their proportions. Previous studies have demonstrated a strong influence of the redox state of the rock in the P-T stability of epidote-, pyroxene- and amphibole-bearing mineral assemblages in eclogite facies metamorphic rocks (e.g. Rebay et al., 2010). Although, we are not in a position to provide an accurate estimate of the O2 content in the modelled rock compositions, the approach adopted here provides a reasonable approximation

ACCEPTED MANUSCRIPT of phase stabilities in HP metabasites as in this study. The metamorphic phase relationships for water-saturated conditions are presented in Figs. 8-9.

3.3. U-Pb zircon geochronology Zircon grains from the two Nagaland eclogite samples (N22F and N23B-1) were analysed for

SC RI PT

U-Th-Pb content using the Sensitive High Resolution Ion Microprobe (SHRIMP-II) instrument at the Research School of Earth Sciences (RSES), The Australian National University (ANU). Zircon grains in the samples were initially identified in polished thin sections using a JEOL JSM 6490 scanning electron microscope (SEM) with EDAX

NU

attachment at the EPMA National facility in IIT Kharagpur. Although small (~10-20 μm), zircon is common throughout the two samples in a variety of textural settings. Back-

MA

scattered-electron (BSE) imaging revealed diverse zircon morphology, internal structure, and mineral inclusions. Because of their small grain size, parts of the thin sections with the most

ED

zircon were cut for in-situ analysis using an ultrasonic disk cutter. These thin section pieces were then mounted on double-sided tape together with Temora 2 reference zircon prior to

PT

casting in 25mm diameter epoxy mounts. The Temora 2 grains were grounded to expose their

CE

mid-sections and polished with 1µm diamond paste prior to mounting together with the thin section pieces. Mounts were then mapped with an automated Leica D6000 microscope in

AC

reflected and transmitted light in order to locate the previously identified zircon grains based on textural positioning in the original thin section. Subsequently, Cathodoluminescence (CL) imaging of zircon grains in the mounts was completed using a JEOL-JSM-6610A SEM at RSES-ANU to characterise their internal structure. For U-Pb isotopic analysis on SHRIMP-II, the grain mounts were coated with ca. of 10 nm of Au. After coating, the epoxy mounts were loaded into a steel holder for insertion into SHRIMP. A reflected light map of each mount was uploaded onto the SHRIMP-II

ACCEPTED MANUSCRIPT computer. Every individual spot location was checked for correct positioning. SHRIMP measurements were performed with a O2- primary beam of ca. 1.5nA, with an accelerating voltage of 10 kV and focused to sputter an area of ~ 10 μm x 15 μm (Köhler illumination). Secondary ions were extracted using an accelerating voltage of 10 kV and measured by single collector analysis on an ETP electron multiplier. Masses were selected by magnetic

SC RI PT

peak switching. SHRIMP-II is operated at a mass resolving power of

~ 6500 (at 10%

peak height). The acquisition time for each spot was ca. 16 min, which consists of six cycles through the following peaks: 248

ThO+, and

254

196

Zr2O+,

204

UO+. The peak position of

Pb+, background, 204

206

Pb+,

207

Pb+,

Pb is fixed relative to

196

208

Pb+,

238

U+,

Zr2O, and the peak

NU

positions of 207Pb and 208Pb are fixed relative to 206Pb.

Uranium concentrations are determined by normalizing to SL13 zircon (238 ppm U;

MA

Roddick and van Breemen, 1994) which were measured at the beginning of every session. The analysis of every five unknown spots was both preceded and followed by the analysis of

ED

Temora2 standard and used as a U/Pb normalizing reference (416.8 ± 1.3 Ma; Black et al., 2003). All analyses were corrected for common Pb using the measured

204

Pb content. Data

PT

were reduced using POXI, an in-house ANU program that follows the data reduction

CE

methods of SQUID 1. The program Isoplot/Ex version 3.0 (Ludwig, 2003) was used to generate probability density and concordia plots and explore possible mixed age populations

AC

through the "mixture modelling" approach of Sambridge and Compston (1994). Single analyses in the data table (Table 6) and concordia plots are presented with 1σ errors.

4. Petrography and mineral chemistry 4.1 Petrography 4.1.1. Thewati eclogite

ACCEPTED MANUSCRIPT In the Thewati eclogite sampleN23B, garnet (Fig. 4a-b), omphacite (Fig. 4c) and epidote (Fig. 4a) occur as large to medium-sized porphyroblasts, enveloped by foliated phengite (Fig. 3d), barroisite (Fig. 4d), epidote and sphene. Rutile, quartz and omphacite occur as inclusions in garnet (Fig. 4a) and omphacite (Fig. 4c). Evidence for an early prograde stage is indicated by foliated, fine-grained glaucophane (glaucophane1), epidote (epidote1), sphene (sphene1),

SC RI PT

phengite and quartz inclusions concentrated in garnet grain cores (Fig. 4a-b), and also rarely in omphacite (Fig. 4c). When present, inclusions in garnet grain rims are comparatively coarse-grained and dominated by omphacite and rutile, and only very rare epidote (Fig. 4a). The peak metamorphic assemblage for this sample is inferred to have involved garnet,

NU

omphacite, epidote, rutile, phengite and quartz.

The composite Thewati eclogite (sample N40B) contains garnet similar to that in

MA

sample N23B (Fig. 4f), and a second variety of garnet that is characteristically smaller (1-2 mm grain diameter) and has an idioblastic habit (Fig. 4e). These garnet types are named

ED

garnet-A and garnet-B, respectively. The cores of garnet-B grains have inclusions of epidote (epidote1) and omphacite (Fig. 4e), but in smaller proportions than occur in garnet-A (Fig.

PT

4f). The outer rims of garnet-B grains are inclusion free. Matrix pyroxene, amphibole and

CE

epidote envelop both garnet-A and garnet-B grains. The post-peak history of the Thewati eclogite involved a series of retrograde stages,

AC

that successively formed: (a) structurally-late epidote (epidote2) with inclusions of garnet and omphacite as aggregates and discontinuous layer-parallel segregations (Fig. 4g, h); (b) the syn- (Fig. 4i), to post-tectonic (with respect to the pervasive matrix foliation) growth of barroisite that partially pseudomorphs garnet, epidote2, and omphacite (Fig. 4d, j–k); (c) a late generation of idioblastic glaucophane (glaucophane2) and sphene (sphene2) that form rims on barroisite (Fig. 4d, l) and rutile, respectively; (d) post-tectonic rims of chlorite and phengite on garnet (Fig. 4m) and (e) cross-cutting veins of lawsonite, albite, winchite-

ACCEPTED MANUSCRIPT ferrowinchite, aegirine-augite, and sphene2 (Fig. 4n-q). The minerals constituting individual veins varies along their length and they can be classified into three types: (i) albite + winchite + aegirine-augite + sphene veins cutting garnet-A grain cores of (Fig. 4n-o); (ii) lawsonite + aegirine-augite veins cutting grains of garnet-B (Fig. 4p) and (iii) lawsonite + albite + ferro-

SC RI PT

winchite veins cutting barroisite porphyroblasts (Fig. 4q).

4.1.2. Mokie eclogite

The Mokie eclogite records prograde, peak and retrograde stages (Fig.2g, 3r) similar to the Thewati eclogite. The rock lacks the texturally late epidote, and late lawsonite, albite, sphene,

prograde history of the Mokie

NU

(ferro-)winchite, aegirine-augite-bearing veins observed in the Thewati eclogite. The sample is preserved by inclusions of glaucophane1 in

MA

omphacite (Fig.4s) and sphene1 in rutile (Fig. 4t). The peak mineral assemblage is inferred to have involved garnet, omphacite/aegirine-augite, rutile, phengite and epidote (Fig.4r). In

ED

contrast to the porphyroblastic garnet in the Thewati eclogite, garnet in the Mokie eclogite occurs as small idioblastic grains, and occasionally as clusters within pyroxene and the rock

PT

matrix (Fig. 4r-s). An intergrowth of barroisite and sphene partially rims pyroxene and

CE

encloses garnet, rutile and epidote (Fig. 4s-u). Barroisite in such intergrowths is in turn

AC

partially pseudomorphed by a thin rim of glaucophane2 (Fig. 4v) and chlorite (Fig. 4u).

4.2. Mineral chemistry Garnet in the Thewati eclogite is dominantly a solid solution among spessartine (Sps 11 to Sps00-01), pyrope (Prp06 to Prp17), almandine (Alm54 to Alm64) and grossular (Grs19 to Grs28) with minor andradite (Adr02 to Adr05) contents (Table 1). Mn-, Mg-, Fe- and Ca-X-ray element images of three representative large garnet grains in samples N23B and N40B (cf. garnet-A in N40B-6) show a uniform feature of Mn-rich core and Mn-poor rim and a

ACCEPTED MANUSCRIPT continuous rim-ward enrichment in Mg contents (Fig. 5a-b, e-f, i-j). The Fe and Ca images show a feature of complementary increase and fall in their concentrations up to the inner rim followed by fall and rise towards the outer rim respectively (Fig. 5c-d, g-h, k-l). The compositional zonation of these elements, inferred from the X-ray element images is also confirmed in detailed compositional profiles (Fig. 5q-ae). The rim-ward depletions in

SC RI PT

spessartine (11-06→00-01 mol%) and enrichments in pyrope (06-08→15-17) and XMg (1113→20-23) contents in garnet (Table 1) are attributed to growth during prograde metamorphism. Characteristic features of elemental zoning in these garnet grains include discontinuities in almandine [54-58 (core)→60-64 (inner rim)→55-61 (outer rim)] and

NU

combined grossular + andradite [26-30 (core)→22-24 (inner rim)→24-27 (outer rim)] profiles (Fig. 5s-t, x-y, ac-ad). On the basis of these compositional features, we consider

MA

garnets in N23B to be comparable with garnet-A in N40B. Type-B (smaller) garnet grains in N40B also show features of growth zoning in the

ED

X-ray element images (Fig. 5m-p) and compositional profiles (Fig. 5af-aj). However, their

PT

grain cores are deficient in spessartine (Sps02 vs. Sps11-06) but enriched in pyrope (Prp13 vs. Prp06-08) and XMg (17 vs. 11-13) contents compared to that in the garnet-A cores (Table 1).

CE

These distinctions are attributed to the growth of garnet-B in sample N40B from a

AC

fractionated bulk rock composition (e.g. in Mn) following the early growth of (the cores of) large garnet grains, and at elevated temperature conditions. However, the spessartine-poor and pyrope-rich compositions in the garnet outer rims are the same in all the textural varieties of garnet. Garnet-B, is also distinct from the other garnet grains by its developments of a narrow (~50μm wide) rims marked by sharp fall in pyrope and XMg and complementary rise in almandine contents (Fig. 5ah-aj) at its outer edge in contact with the matrix omphacite (Fig.4e). These are attributed to localized diffusive re-equilibration.

ACCEPTED MANUSCRIPT Garnet in the Mokie eclogite is characteristically Mn-rich (Sps14-26) and growthzoned [Sps26 Prp07 Alm45 Grs19 (C)→Sps14-17 Prp10-12 Alm52-57 Grs13-23 (R)] as in the Thewati eclogite (Table 1). The XPrp and XMg contents in its outer rim are, however, lower than that in the Thewati eclogite. Pyroxenes in the Nagaland eclogites show a variation in composition from aegirine-

SC RI PT

augite to omphacite (Fig. 6a-b, Table 2), which can be related to their textural setting. In the Thewati eclogites, pyroxene inclusions (Type 1 pyroxene) in the cores of large garnet grains straddle the boundary of the aegirine-augite and omphacite fields (Fig. 6a). Pyroxene in the garnets rims is omphacite (Fig. 6a). Pyroxene inclusions in the smaller garnet-B grains in

NU

N40B are omphacite (Fig. 6a). Pyroxene grain in the matrix (Type 2A pyroxene) or its relict within the texturally late barroisite grains (Type 2B pyroxene) is omphacite (Fig. 6b). The

MA

pyroxene porphyroblasts in the Mokie eclogite are omphacite to aegirine-augite (Fig. 6b). The vein-type pyroxenes are characteristically more oxidised and plot in the field of aegirine-

ED

augite in the ternary Q (Quad)-Jd (Jadeite)-Aeg (Aegirine) diagram (Fig. 6c). Na-Ca amphibole in the Thewati eclogite is barroisite. The Mokie eclogite has

PT

amphibole ranging in composition from barroisite to magnesio-kataphorite (Fig. 6d-e, Table

CE

3). The Na-Ca amphiboles in the veins are ferro-winchite (Fig. 6d). The texturally early and later sodic amphiboles in the Nagaland eclogites are glaucophane (Fig. 6f). In the binary NaB

AC

vs. (Al+Fe3+)VI plot, both the sodic and sodic-calcic amphibole groups show compositional spreads along the glaucophane substitution vector (Fig. 6g). The sodic-calcic amphibole group additionally shows minor tschermak substitution (Fig. 6g). Ferro-winchite in the veins is comparatively enriched in NaB content and lies very close to the boundary between sodiccalcic and sodic amphibole groups (Fig. 6g).

ACCEPTED MANUSCRIPT Phengite in the Thewati eclogite sample N23B is more siliceous [Si (c.p.f.u.) = 3.533.54] than that in sample N40B (Si = 3.46). Phengite of even higher Si-content is found in the Mokie eclogite sample N22F (Si = 3.63). In the Thewati eclogite, epidote1 inclusions within garnet are more pistacite-rich

the Mokie eclogite is pistacite-poor (XPs = 0.17).

SC RI PT

(XPs=0.25) than the foliation-defining epidote2 (XPs= 0.16-0.20). Texturally early epidote in

Chlorite coronae around garnet in the Thewati eclogite is clinochlore (XMg=0.68). Sphene1 and sphene2 in the Thewati and Mokie eclogites contain small amount of alumina (Al2O3 = 1.04-1.43 wt %). Albite, lawsonite and rutile are near pure phases.

NU

Summarising, the textural and mineral compositional features of the Nagaland eclogites reveal a multistage mineral assemblage evolution comprising: (1) early epidote

MA

blueschist facies conditions that accompanied prograde garnet growth (Fig. 5) together with Na-Ca pyroxene (Fig. 6a), glaucophane1, epidote1, sphene1, phengite and quartz; (2) peak

ED

conditions represented by garnet outer rims of spessartine-poor and pyrope-rich compositions, jadeite-rich omphacite in the matrix and inclusions in the garnet outer rim,

PT

rutile and phengite with or without epidote and quartz; and (3) a protracted post-peak

CE

evolution marked by sequential appearances of epidote2 (Fig.4g-h), barroisite + sphene2 ± epidote2 (Fig. 4i-k, s-u) and glaucophane2 + chlorite (Fig.4k-l, m, t-v). The Thewati eclogite

AC

additionally records the development of lawsonite with or without albite, ferro-winchite and aegirine-augite in cross-cutting veins (Fig. 4n-q).

4.3. Peak P-T conditions of metamorphism Mineral compositions from the pyrope-rich, growth-zoned garnet rims were paired with jadeite-rich omphacite from either the cores of larger matrix grains or inclusions in garnet rims together with the high-Si matrix phengite to estimate peak metamorphic

ACCEPTED MANUSCRIPT conditions. According to equations R1 and R2, these chemical parameters most likely reflect the maximum pressure conditions experienced by eclogite facies basites. For the Thewati eclogites, the applications of Method 1 yielded peak P-T estimates of 28.4 ± 1.6 kbar (weighted mean value with errors at 2 sigma) and 630 ± 50°C for the sample N23B, which are similar to those obtained by the Method 2 (28.1 ± 2.8kbar, 645 ± 80 °C) (Table 4). For

SC RI PT

the sample N40B, the estimated pressures (P = 26.1 ± 2.0 kbar, Method 1 and 25.5 ± 3.3 kbar, Method 2) are lower by ~2.3 to 2.6 kbar at similar TMax (T= 635±60°C, Method 1 and T= 650±100 °C, Method 2). We attribute the higher P-estimates in N23B to a more siliceous composition of phengite in the sample [Si = 3.54 (c.p.f.u.) vs. 3.46]. Compared to the

NU

Methods 1 and 2, the Method 3 yielded temperatures, which exceed by ~25 oC in sample N23B and ~ 100 oC for N40B. In particular, the peak estimate of ~755 oC in N40B is much

MA

higher than the more robust garnet-clinopyroxene Fe2+-Mg exchange thermometry results. For this reason, we discount the results of the Method 3 from further discussions on the P-T

ED

constraints of peak metamorphism in the Nagaland eclogites. Although the results from the first two methods are very close, we favour the computations by the Method 2. This is

PT

because these P-T values can be directly compared with the results of the phase equilibria

CE

modelling as the two methodologies adopt the same thermodynamic dataset and mineral solid solution models. Accordingly, we consider P ~ 25.5 to 28.1 kbar, T ~650 oC as the peak P-T

AC

estimates of the Thewati eclogites. Using the same reasoning, the peak P-T estimate of the Mokie eclogite has been constrained at P ~ 23.8 ± 5.1 kbar, T~ 555 ± 90 oC (Table 4).

4.4. Calculations of metamorphic phase diagrams and deduction of metamorphic P-T path Phase diagrams calculated for the Thewati eclogite are presented in Figs. 8 & 9. For the pyroxene-garnet-rich domains of eclogite sample N40B (cf. N40B-4), the phase diagram predicts garnetiferous assemblages at T>475-550 oC, and the loss of chlorite at T>550-575 oC

ACCEPTED MANUSCRIPT and amphibole at T>500-725 oC, depending on metamorphic pressures (Fig. 8a). The positive dP/dT slope of the lawsonite-out curve and the stability of the lawsonite-bearing assemblages at T<525-650 oC is nearly identical to the experimentally determined lawsonite-out curve (e.g. Liou, 1971; Schmidt, 1995). Epidote-bearing assemblages define an inverted V-shaped field at P<25 kbar and at T conditions between 525 and 675 oC. The early prograde loss of

SC RI PT

lawsonite before epidote during heating at P<25 kbar is consistent with phase diagram predictions made for metabasic bulk rock composition elsewhere (Wei and Clarke, 2011). The phase diagram also predicts peak eclogite facies conditions of T>625-675oC, which are consistent with the estimates for the sample N40B (P~ 25.5 kbar, T~650 oC) (Fig. 8a).

NU

We now evaluate the progressive changes in garnet composition (Fig. 8b) and garnet (Fig.8c), pyroxene (Fig. 8d) and amphibole (Fig. 8e) modes along three model P-T paths:

MA

Path1, a prograde path across the epidote-blueschist facies; Path 2, a counter-clockwise (CCW) P-T trajectory involving cooling and decompression; and Path 3, a clockwise (CW)

ED

P-T trajectory involving decompression followed by cooling. Path P1 predicts (a) the sequential intersections of the garnet-in, chlorite-out, amphibole-out and epidote-out curves;

PT

(b) prograde garnet growth with garnet composition becoming progressively magnesian

CE

(decreasing x(g) isopleths, see Table 4 for the terminology) and Ca-poor (decreasing z(g) isopleths) and (c) increasing omphacite and decreasing amphibole contents in the rock. This

AC

is consistent with the observed textural and garnet compositional features. Assuming the rock was not appreciably fractionated prior to the growth of garnet, the P–T condition at which garnet is predicted to first formed can be estimated using garnet core compositions and the isopleth thermobarometry approach (after Vance and Mahar, 1998). For this purpose, we have used the core composition of the largest garnet-A grain in sample N40B, which has the highest spessartine, lowest pyrope and highest grossular + andradite content. The x(g) =0.85 and z(g) =0.30) garnet isopleths intersect at ~18.8kbar, 555oC (Fig. 8f). This P-T estimate

ACCEPTED MANUSCRIPT falls in the epidote blueschist facies, in good agreement with the reconstructed mineral reaction history, and can be considered a reasonable approximation of the conditions accompanying early garnet growth for the Thewati eclogite. An expanded MnO-bearing system would likely displace this result to subtly lower grade (see, for example, Marmo et al., 2002). Both of the retrograde paths P2 and P3 (Fig. 8c) can be expected to induce the

SC RI PT

decomposition of garnet and clinopyroxene, and growth of epidote and amphibole. Path 2 posits the growth of lawsonite before amphibole (Fig. 8f), inconsistent with the deduced mineral reaction history. In contrast, Path 3 posits the sequential growth of epidote, then epidote with amphibole, then chlorite (Fig. 8f), in good agreement with the deduced sequence

NU

of retrograde minerals.

The phase diagram for the amphibole-rich compositions of sample N40B resembles

MA

that for the pyroxene-rich compositions in having similar P-T stability fields for lawsoniteand chlorite-bearing assemblages (Fig. 9a. However, the P-T range of epidote-bearing

ED

mineral assemblages is significantly enlarged (Fig. 9a), due to the higher calcic and alumina content of the modelled composition (Table 5). The phase diagram is contoured for garnet

PT

composition and garnet, pyroxene and amphibole modes (Fig. 9b-e). Whereas the phase

CE

diagram broadly predicts prograde and retrograde growth sequences similar to that of the garnet-pyroxene-rich composition (compare Figs 8f and 9f), the amphibole-rich bulk

AC

composition better matches discontinuities in grossular zoning during the prograde garnet growth. The prediction of a progressive decrease in z(g) content from core to the inner rim followed by an increase towards the garnet outer rim (along P1, Fig. 9b) resembles the measured garnet composition (Fig. 5t, y, ad). The phase diagram is also contoured with isopleths of NaB content in amphibole for values 0.5 and 1.5 [where the numbers refer to cations per formula unit (c.p.f.u.) in the B structural site of amphibole] (Fig. 9e). When considered together with the positions of

ACCEPTED MANUSCRIPT chlorite-in and sphene-in curves, the phase diagram with the modelled amphibole compositions predict the observed growth of glaucophane after barroisite, chlorite after garnet and sphene after rutile (Fig. 4d, l-m) along the retrograde path P2 (Fig. 9f). This gives confidence in using the phase diagram together with detail from chemically-zoned amphibole grains to assess the validity of the retrograde paths P2 and P3 (Fig. 9f).

SC RI PT

We now combine the reconstructed metamorphic evolutionary pathways, key details of the phase equilibria modelling, estimated peak P-T conditions of metamorphism and P-T positions of the experimentally determined lawsonite-out curve and calculated mineral equilibria, corresponding to greenschist-blueschist-facies transition and transformation of

NU

albite to jadeite + quartz to establish the metamorphic P-T path of evolution from the Thewati eclogite sample (Fig. 10). We additionally calculate and plot four activity-corrected end

MA

member mineral reactions (reactions R3 to R6) to supplement phase equilibria predictions of mineral stabilities, particularly involving retrograde minerals that formed in localised

ED

domains. These are:10 (jadeite)Omphacite + (4 grossular + 5 pyrope)Garnet + 11 quartz + 8 H2O = 5 (glaucophane)Barroisite+ 6 (clinozoisite)Epidote (R3).

PT

5 (glaucophane)Barroisite+12 (clinozoisite)Epidote + 14 quartz + 28 H2O = 3 (tremolite)Fe-winchite +

CE

10 albite + 18 lawsonite (R4).

3 daphnite + 5 pyrope = 3 clinochlore + 5 almandine (R5).

AC

Ab = Jd17 + Qtz (R6).

These reactions reveal the thermal and baric conditions for the stability of retrograde barroisite + epidote, chlorite and lawsonite, albite, ferro-winchite and aegirine-augite in late veins. In calculating R3-R5, we have applied the same thermodynamic dataset and activitycomposition models of mineral end members as adopted during the phase equilibria modelling. This makes the P-T positions of these equilibria to be internally consistent and comparable with the results of phase equilibria modelling.

ACCEPTED MANUSCRIPT The disposition of the sphene-out, garnet-out, chlorite-out, amphibole-out and epidote-out curves in sequence and the P-T conditions of formation of the prograde garnet core at 18.8 kbar, 540 oC, when considered together provide a tight constraint for the prograde P-T path (cf. MPR metamorphic stage, where PR represents prograde) that culminated in peak eclogite facies metamorphism at 25.5 kbar, 650 oC (MP metamorphism as

SC RI PT

P represents peak, Fig. 10). The intersections of epidote-out and amphibole-out curves in an opposite sense and reaction R3, leading to the successive appearances of retrograde epidote2 and barroisite2 + epidote2 constrain the first stage of retrogression (cf. MR1 metamorphic stage, where R refers to retrograde) as one of decompression with minor cooling to ~18 kbar,

NU

~620 oC (at amphibole eclogite facies condition) (Fig. 10). The subsequent crossing of the NaB =1.5 isopleth in amphibole and intersections of the chlorite-out and sphene-out curves in

MA

the reverse sense and R5 indicates further decompression and cooling in the epidote blueschist facies at ~10kbar, 510 oC (cf. MR2 metamorphic stage) (Fig. 10). The dispositions

ED

of the reactions, R4 and R6, reflecting the stability of aegirine-augite + lawsonite + albite + ferro-winchite assemblage indicate a terminal phase of retrogression to one of protracted

PT

cooling (∆T ~200 oC) at low-pressures (P~6 kbar) (cf. M3R metamorphic stage, Fig. 10). The

CE

development of a narrow diffusion zoning in garnet-B (Fig.5aj) in contact with omphacite provides additional support for post-peak cooling. Summarising, our study reveals that this

AC

progressive change in P-T vectors from loading accompanying heating (MPR) culminating in peak metamorphism (MP) and three successive stages of decompression and cooling (MR1 to MR3) defines a clockwise P-T path of evolution of the Thewati eclogites (Fig. 10). Although we have not quantified the metamorphic P-T path of evolution of the Mokie eclogites, the broadly analogous metamorphic reaction pathways in the Mokie eclogite with the Thewati eclogite raise the possibility that the former also records a broadly similar subduction burialexhumation history in the Nagaland subduction channel.

ACCEPTED MANUSCRIPT

5. Results of SHRIMP U-Pb zircon geochronology

Analytical data for zircon are presented in Table 6. Zircon textural and cathodoluminescence (CL) features are summarized in Fig. 11 and chemical composition and Th/U vs. age

SC RI PT

distribution shown in Fig. 12. A total of 34 spot analyses from 33 zircon grains were completed on the two eclogite samples (17 grains in N22F and 16 grains in N23B-1). In both samples, zircon grains are common but small (grain diameter range ~ 8-20 μm). They occur in various textural settings: as inclusions in garnet (Fig. 11a), Na-Ca pyroxene (Fig. 11b) and

NU

rutile, and also in the matrix foliation-defining barroisite, phengite and sphene. Zircon in N22F generally occurs as elongate grains with pyramidal terminations and has grain cores

MA

dominated by ovoid inclusions (Fig. 11c). Zircon in N23B-1 mostly occurs as equidimensional grains with well-developed crystal faces that are free of mineral inclusions (Fig.

ED

11d). Where inclusions are relatively coarse they are determined to be garnet (Fig. 11e), omphacite (Fig. 11f), and rutile.

PT

Their small size partially inhibits the high resolution CL images of zircon in the two

CE

Nagaland eclogites. However, they have a variety of CL structures (Fig. 11g-x) that can be used to identify three types of CL domains that reflect stages of grain growth,

AC

recrystallisation and alteration (Table 7). Type 1 zircon is marked by darkly luminescent, faintly banded (Fig. 11g, i, n-p) elongate domains that characteristically have euhedral shapes and pyramidal terminations (Fig. 11j-l, r-t). These features are consistent with zircon growth in the presence of a fluid or melt (e.g. Rubatto, 2017). Type 2 zircon is marked by (a) structure-less grain domains with morphology similar to type-1 zircon (Fig. 11h, j-m, r, u), commonly with patchy, irregular zones of bright luminescence (Fig. 11q, s-t); (b) textural modification of zircon with type-1 grain form involving grain rounding and/or the truncation

ACCEPTED MANUSCRIPT of pyramidal faces (Fig. 11j, l, p, r-s, u-v); and (c) development of bright luminescent narrow rims to type-1 zircon (Fig. 11h, j-m, p-v), commonly with equant and euhedral external shape oblique to the enclosed type-1 zircon core (Fig. 11j-k). These features are attributed to the recrystallisation and/or replacement of type 1 zircon. Type 3 zircon only occurs in N22F and is marked by irregular bright-luminescent domains that partially replace type-1 zircon cores

SC RI PT

(e.g. domain 1 in grain Z-9) (Fig. 11v). The zones of replacement have serrated to curvilinear inner boundaries (e.g. in grains Z-9, 33) (Fig. 11v-w). In some zircon grains, the partial replacement of the older core is extensive, and gives rise to ghost and fuzzy zoning (e.g. grain Z-28) (Fig. 11x). Features of this third stage are consistent with grain modification during the

NU

localised fluid-mediated alteration of type 1 zircon.

Types 1 and 2 zircon grains (both from primary and recrystallised varieties) in N23B-

MA

1 have limited compositional variation with low Th (0.05-1.62 ppm) and U (29-200 ppm) contents (Fig. 12a). Equivalent zircon domains in N22F are characteristically uraniferous (U:

ED

1110-2478 ppm) and more thorian (Th: 8-66 ppm) (Fig. 12a). A Type 2 zircon rim in N23B-1 is depleted in U relative to the recrystallised grain core (e.g. 13 ppm in rim vs. 29 ppm in

PT

core, grain Z-66, N23B-1) (Fig. 12a). Type 3 zircon is chemically distinct by being more

CE

thorian (Th: 94-271 ppm) and has high Th/U ratios (0.09 to 0.20) (Table 6, Fig.12a-b). Inter-sample compositional variability of Th-U in the Nagaland zircon grains is a

AC

feature common in many HP metamorphic rocks. It is generally attributed to zircon growth in the absence of a dominant buffering fluid (Rubatto and Hermann, 2007). Type 1 zircon in the Nagaland samples also has low Th/U ratios (≤0.02-0.03) (Table 7,Fig. 12b), that are common in metamorphic zircon from blueschist and eclogite facies rocks (Rubatto and Gebauer, 2000; Rubatto and Hermann, 2007). Given the small size of Nagaland zircon grains, it is not possible to analyse most type2 zircon rims in the two samples. One analysis was made of a comparatively thick rim in

ACCEPTED MANUSCRIPT N23B-1 (see below). The SHRIMP analysis sessions thus targeted zircon core domains, including the type1 zircon, and additionally the type3 zircon in N22F (See Fig. 11 for the analytical spots in representative zircon grains). Among the 17 analysed spots in N22F, six spots were excluded from the calculation of concordia and mean ages based on the following considerations: (a) these analyses reveal discordance in the concordia plot (Fig. 13a); (b) four

SC RI PT

of the analyses are from the Th-rich altered domain (type 3 zircon); and (c) all have yielded older dates (between 337 and 234 Ma). Following previous studies that have indicated the isotopic disturbance of altered zircon grains in HP metamorphic terranes, we are of the opinion that zircon ages from such texturally late, fluid-mediated alteration in N22F are

NU

likely geologically spurious and lack geological significance (e.g. Rubatto and Hermann 2003).

MA

The remaining 11 analyses show a spread in zircon ages from 216 to 171 Ma (Table 6) and reveal three groups at 206 ± 4 Ma (error 2s), 189 ± 3 Ma, and a minor one at 172 ± 4

ED

Ma (Fig. 13b). Concordia ages of the two dominant age groups are 205± 6 Ma [n=4, Mean Square Weighted Deviates (MSWD) =0.49, probability of concordance =0.48] (Fig. 13c) and

PT

189 ±1Ma (n=5,MSWD=0.05, probability=0.82) (Fig. 13d).

CE

In N23B-1, 18 analyses show a spread in zircon dates from ~217 Ma to ~152 Ma, with two broad age groups at 204 ± 10 Ma and 187 ± 7 Ma (Fig. 13e) similar to that observed

AC

in N22F. Two analyses from core-rim domains of one zircon grain (grain Z-66, Fig. 11h) give dates at 190 ± 18 Ma for the grain core (type 1 zircon) and 152 ± 42 Ma for the grain rim (type 2 zircon). Although this is consistent with zircon recrystallisation being younger than the growth of the core, the large age error due to the presence of high common Pb does not allow this age to be useful for estimating the timing of the recrystallisation event. For concordia age calculations, we have excluded 9 analyses because of negative or indeterminate radiogenic

207

Pb/235U ratios. The remaining 9 analyses have yielded two

ACCEPTED MANUSCRIPT concordia ages at 205± 4 Ma (n=5, MSWD = 0.10, probability = 0.75) (Fig. 13f) and 185 ± 4Ma (n=4, MSWD = 0.08, probability=0.78) (Fig. 13g) that are near identical to that obtained from N22F.

6. Discussion

SC RI PT

6.1. Early Jurassic eclogite facies metamorphism in the Nagaland Ophiolite Complex Nearly identical zircon age groupings of ~205 Ma, ~185-189 Ma and (minor) ~172 Ma have been identified in two separated eclogite locations in Nagaland. We recognise two stages of zircon growth and/or recrystallisation, in small zircon grains that challenge the analytical

NU

method. Despite minor overlap, the data largely define two end member textural types (1 and 2, below). In the context of petrological observations, the concordia ages of the first two age

MA

groups are most probably geologically meaningful. The oldest zircon age group was recovered from type 1 zircon, and reflects the earliest stage of zircon growth (e.g. in N23B).

ED

There are three possible interpretations of the c. 205 Ma date: (a) the growth of zircon in the eclogitic protolith; (b) a partially reset growth age due to the effects of post-growth heating

PT

during later prograde conditions; and (c) zircon growth during prograde metamorphism. In

CE

the absence of zircon trace element data it is difficult to distinguish between these possibilities, but the characteristic low Th/U ratios (<0.01) of these zircons are most

AC

consistent with the c. 205 Ma date to being a minimum age for the eclogitic protolith [option (b) above]. The c.185-189 Ma zircon ages were recovered mostly from stage 2 recrystallised grain cores, which include eclogite facies mineral inclusions. Stage 2 zircon is thus interpreted to have formed during or shortly after peak conditions in the eclogitic facies, at c. 189-185 Ma. The limited dataset of the c. 172 Ma concordant zircon ages is interpreted as recording the timing of retrograde cooling and exhumation. However, more detailed work is needed to firmly establish its significance. An important outcome of these data is the

ACCEPTED MANUSCRIPT identification of a Triassic age (≥ c. 205 Ma) for this segment of the Neo-Tethyan oceanic crust at Nagaland.

6.2. Thermal history and accretionary processes in an Early Mesozoic Neo-Tethyan subduction channel

SC RI PT

Eclogite from two locations in the Nagaland Ophiolite Complex preserve a record of broadly coeval eclogite facies conditions in the Early Jurassic: ~28.1 to 25.5 kbar, ~ 650 oC for the Thewati body and 23.8 kbar, 555 oC for the Mokie body. Whereas the Mokie eclogite records low-temperature (LT) eclogite facies conditions (terminology and P-T field

NU

of global LT eclogites are after Fig. 3a in Agard et al., 2016), the Thewati eclogites underwent burial to deeper levels (~80-90 km depths) that were transitional to ultra-high

MA

pressure metamorphic conditions (Fig. 14a). In a previous study, Chatterjee and Ghose (2010) reported eclogite from a location near the Satuza area, and inferred a CW P-T path

ED

with peak conditions of 17–20 kbar and 580–610 °C (Path 3, Fig.14a). Whereas the prograde and retrograde segments of the Satuza P-T path resemble that of the Thewati

PT

eclogite (Path 1, Fig. 14a), the Satuza locality records shallower levels and higher thermal

CE

gradients (>10 oC/km) (Fig. 14a). Our findings, however, suggest that the Thewati eclogite falls under the select category of oceanic eclogites that are recovered from the maximum

AC

depth in the subduction zone (PMax ~28 kbar) and with a low apparent peak thermal gradient of ≤6–9 °C/km (Fig. 15). In contrast to the Nagaland samples, jadeite-glaucophane-bearing and garnetphengite-bearing schists from the adjoining Central Ophiolite Belt of Myanmar preserve a much larger spread in peak conditions, generally in blueschist to low-T eclogite facies (reviewed in Nyunt et al., 2017; Fig. 14a). A well-constrained metamorphic P-T path for the evolution of the terrane is lacking. Earlier studies of the Himalayan blueschists (see O'Brien,

ACCEPTED MANUSCRIPT 2018 for the latest review on the Himalayan HP/LT metamorphism) inferred low-T conditions and CW P-T paths (Paths 4 & 6, Fig. 14a) (Guiraud, 1982; Jan, 1985; Honegger et al., 1989; Anczkiewicz et al., 2000; Oberhänsli, 2013). On the basis of a detailed metamorphic investigation of Sapi-Shergol garnetiferous lawsonite blueschists, Groppo et al. (2016) recently established a significantly deeper and cold subduction system of possible

SC RI PT

Cretaceous age (cf. Anczkiewicz et al., 2000) within the Neo-Tethys from the NW Himalayas (peak P-T at 19 kbar, 470 oC) and a clockwise hairpin P-T loop along low thermal gradients (< 8-9 °C/km; Path 5, Fig. 14a). While the geometry of the P-T path and lower thermal gradients of the Shangla lawsonite blueschist resemble with that of the

NU

Nagaland blueschist (Path 2, Fig. 14a), our new geochronological constraints indicate the subduction zone metamorphism in the eastern Neo-Tethys initiated much earlier in the Early

MA

Jurassic.

The lower apparent thermal gradient at metamorphic peak in the Nagaland eclogite

ED

samples (7-8 oC/km) (Fig. 14a) is comparable to the gradient recovered from the associated lawsonite blueschist in Mokie (~9 oC/km, Ao and Bhowmik, 2014) and other HP

PT

metamorphic belts in the Zagros suture zone, Southern Iran (Angiboust et al., 2016), Sapi-

CE

Shergol lawsonite blueschist, Western Himalaya in the YTSZ (Groppo et al., 2016), Dominican Republic and Rio San Juan Mélange Complex (Krebs et al., 2011), the Western

AC

Alps (Angiboust et al., 2009; Agard et al., 2009, 2016 and reference cited therein), the Bantimala Complex, South Sulawesi, Indonesia (Setiawan et al., 2016) and also predicted by thermo-mechanical models (Gerya et al., 2002) (Fig. 15). The inferred Nagaland P-T path is similar to that inferred for the Sulawesi eclogite (Path 6, Fig. 15). For the latter, subduction of an old and cold oceanic crust was suggested (Setiawan et al., 2016). For the present study, we relate this characteristic lower thermal gradient during peak eclogite and blueschist facies metamorphism, the (hairpin) clockwise P-T paths of blueschists and

ACCEPTED MANUSCRIPT eclogites (Fig. 14a) and the geological setting with the development of a cold mature stage of an intra-oceanic subduction system within the Neo-Tethys in Nagaland (see also Ao and Bhowmik, 2014; Bhowmik and Ao, 2016). In many subduction-related HP-LT terranes, different localities of adjacent tectonometamorphic units record a systematic increase of peak-T and associated P, producing a

SC RI PT

positive dP/dT slope of the P-T array (e.g. Fig. 13b of Groppo et al., 2016). In the NOC, our studies (e.g. Ao and Bhowmik, 2014 and this work) demonstrate a near-complete preservation of a progressive metamorphic sequence from greenschist through pumpellyitediopside, blueschist to eclogite facies from the Mokie-Satuza area (Figs. 2c, 16). The peak P-

NU

T estimates of these rocks reveal a metamorphic field gradient (Fig. 14b) that is a characteristic feature of many HP/LT metamorphic terranes worldwide (e.g. Oman: Yamato

MA

et al., 2007; Turkey: Plunder et al., 2015; New Caledonia: Vitale Brovarone and Agard, 2013; Corsica: Vitale Brovarone et al., 2014; Schistes Lustres of the Western Alps: Plunder et al.,

ED

2012), also predicted by thermo-mechanical models (Fig. 14b). We attribute this metamorphic zonation in the NOC to repeated accretion of the top layer of the oceanic crust

PT

subducted at different depths and subsequently off-scraped at the base of the accretionary

CE

prism (e.g. Agard et al., 2009, 2018 and references therein). In the NOC, these HP metamorphic rocks are tectonically imbricated with serpentinites and mafic and ultramafic

AC

cumulate rocks, giving rise to an apparent chaotic internal arrangement in the ophiolitic mélange, though we note the coherence of individual members of the HP metamorphic sequence, in particular from the Mokie-Satuza area. The Nagaland P-T path reveals a two stage exhumation of the Nagaland eclogites: (a) initial exhumation (from 25.5 kbar, 650 oC to 10 kbar, 510 oC, and combining MR1 and MR2 metamorphic stages) took place along a steep dP/dT gradient (Fig. 14a). Such an exhumation path resembles that deduced for cold oceanic eclogites from different metamorphic belts

ACCEPTED MANUSCRIPT worldwide (e.g. Fig. 15); (b) later exhumation of the eclogites took place along a gentler dP/dT gradient (cf. MR3 metamorphic stage, Fig. 14a), and is distinctive for the Nagaland rocks. Three tectonic scenarios have been invoked to account for the exhumation of deeply buried eclogite along a steep dP/dT gradient (see Agard et al., 2009 for a general review). (1) In the first scenario that involves high-grade blocks (amphibolites/eclogites along CCW P-T

SC RI PT

path), the exhumation takes place within the first 5–15 Myrs of subduction infancy, in rapidly cooling subduction zones (e.g. Agard et al., 2009). In this scenario, exhumation is triggered by rheological weakening of the subducted oceanic crust due to partial melting in warm subduction zones (Sorensen and Barton, 1987). It has been proposed that this leads to the

NU

boudinage and detachment of the high-grade blocks in a weaker, sedimentary mud-matrix mélange (Cloos and Shreve, 1988). (2) This scenario involves buoyancy-driven exhumation

MA

along the subduction channel as a consequence of serpentinisation of the mantle wedge, controlled by fluid released during dehydration of hydrous phases in the subducting oceanic

ED

crust (Cloos, 1982, 1985; Gerya et al., 2002; Hermann et al., 2000; Guillot et al., 2009). In this setting, oceanic blueschist and eclogite are thought to occur as thin, disrupted bodies in

PT

association with serpentinite and/or a mechanically weak mud matrix (e.g. Gerya et al.,

CE

2002). The scenario predicts the exhumation of eclogite as part of material transport in a large-scale convective circulation system in the subduction channel that is evolving to a cold

AC

and mature stage (Cloos, 1982; Peacock & Wang, 1999; Gerya et al., 2002). The exhumation is inferred to be facilitated by the buyoancy and low viscosity of the serpentinites in the channel (Guillot et al., 2009). (3) A third scenario that involves the exhumation of eclogite as a consequence of the positively buoyant continental crust entering the subduction channel towards the end of the oceanic subduction, allowing the extrusion of large volumes of deeply buried continental crust (e.g. Norway, Andersen et al., 1991; Oman, Chemenda et al., 1996; Dabie Shan, Hacker et al., 1996; Zermatt-Saas eclogite, Western Alps, Angiboust et al.,

ACCEPTED MANUSCRIPT 2009; New Caledonia, Fitzherbert et al., 2004). In this scenario, oceanic eclogite is exhumed as part of a forced return flow of buyoant continental material with or without slab break-off (Von Blanckenburg and Huw Davies, 1995). For the Nagaland eclogites, features such as the (a) characteristic LT eclogite facies metamorphism along a CW P-T path, (b) occurrence of eclogite as rare, small lensoidal bodies within blueschist, (c) presence of thin, fluid-assisted

SC RI PT

metasomatic bands rich in amphibole + white mica as rims around the eclogitic cores and (d) absence of deeply buried continental crustal fragment in the ophiolitic mélange are collectively most consistent with exhumation via the second scenario.

The stage 2 exhumation that involved protracted cooling at shallower levels (Fig. 14a)

NU

can be linked with thrust stacking of partially exhumed eclogites with the cooler, prograde blueschists at shallower crustal levels. Despite differences in peak P-T conditions of

MA

metamorphism in the eclogites and blueschists, the P-T conditions of terminal retrogression in these rocks are broadly similar (P~5-6 kbar, T~300 oC, Fig. 14a). This suggest that the

ED

eclogites and blueschists were finally tectonically juxtaposed along their exhumation paths. On the basis of this deduction, we propose that the final mixing of rocks of contrasting

PT

lithologies, originally placed at different depths in the subduction channel, and producing a

CE

tectonic collage took place at a shallower crustal level during the development of the Nagaland accretionary complex. Future work should be directed to provide structural details

AC

in support of the suggested model above.

6.3. Regional correlations and Early Mesozoic tectonic framework of SE Asia The Nagaland eclogite samples dated in this work are older than the Early Cretaceous NeoTethyan ophiolites of the WOB (e.g. Manipur and Kalaymyo ophiolite belts) and its amphibolite-facies metamorphic sole (119-115 Ma, Liu et al., 2016a; Zhang et al., 2017). The Nagaland Ophiolite also predates the Mid-Jurassic aged Myitkyina ophiolite of the EOB by

ACCEPTED MANUSCRIPT more than ~30 Myr. It is also older than Cretaceous blueschists (e.g. 100-80 Ma Sapi-Shergol and Shangla blueschist, NW Himalayas, Anczkiewicz et al., 2000) and amphibolites (e.g. Saga amphibolite metamorphic sole between 135 and 127 Ma, Guilmette et al., 2012) (locations in Fig. 1 of the YTSZ). Nevertheless, it is contemporary with Late Triassic-Early Jurassic arc magmatism in the Gangdese magmatic belt (e.g. Yu et al., 2018). We propose

SC RI PT

that this cryptic record for the earliest phase of Neo-Tethyan subduction can be traced westwards from Nagaland-Myanmar, along the southern Eurasian margin to Pakistan as part of an Early to Middle-Jurassic Andean-type convergent margin (Searle et al., 1999; Yu et al., 2018).

NU

Several recent studies present interpretations for the Early Mesozoic tectonic framework of northern Myanmar and southeastern Tibet (see Searle et al., 2017; Gardiner et

MA

al., 2018; Zhang, J. et al., 2018 for general reviews). However, there remains considerable uncertainty regarding the: (a) polarity (easterly, Zhang, J. et al., 2018 or westerly, Mitchell,

ED

1993) and type of the subduction (single long-lived, e.g. Zhang, J. et al., 2018; or multiple, Mitchell, 1993; Acharyya, 2007, 2015; Mitchell et al., 2015; Searle et al., 2017; Gardiner et

PT

al., 2018); (b) status of the Western Myanmar block [e.g. a western extension of the Shan

CE

Plateau, an extension of the Lhasa Block, through Tengchong and Baoshan in Yunnan or as a part of the west Sumatra block (e.g. Zhang, J. et al., 2018)] and (c) any continuity of the

AC

three major Mesozoic Tibetan sutures into Myanmar. Liu et al. (2016a) interpret the EOB as a Jurassic suture of the Meso-Tethys (equivalent to the BN suture in central Tibet). This extended suture would pass between the Sibumasu block to the east and the West Myanmar block to the west (Fig. 17a). Gardiner et al. (2018) recently proposed a southerly termination of the Meso-Tethyan BN suture in NE Myanmar, extending the suture between the Tengchong block to the west and the coalesced Sibumasu–Baoshan Blocks to the east (Fig. 17b). In this model, an east-subducting Meso-Tethyan sea floor developed an Andean-type

ACCEPTED MANUSCRIPT active continental margin setting above the Sibumasu block between the Late Triassic and Early Jurassic. The closure of this Meso-Tethyan subduction through the amalgamation of the two continental blocks is inferred to have occurred before the Late Cretaceous. Alternatively, Zhang, J. et al. (2018) propose a long-lived east-dipping subduction within the Neo-Tethyan Ocean that progressively migrated westwards from the Early Jurassic – Middle Eocene and

SC RI PT

the present day. In this model, the subduction of the Neo-Tethyan ocean beneath the Sibumasu continental block during Early to Middle Jurassic time resulted in an accretionary complex of HP metamorphic rocks (now represented by the Jade Mines Belt and Kumon range), and a broadly coeval continental margin Jurassic magmatism in the Myitkyina

NU

ophiolite and Mogok metamorphic belt (Fig. 17c). However, this model predicts peak and retrograde metamorphism in the NOC, which is situated west of the Jade Mines Belt, to be

MA

the product of a second accretionary event in the Late Cretaceous. The reconstruction of an Early Jurassic accretionary complex in the NOC (this study)

ED

raises the possibility that it was continuous with the accretionary complex of the Jade Mines Belt further northeast (with respect to current geographic co-ordinates; Fig. 17c). The two

PT

ophiolites have been disrupted by younger strike-slip faulting along a structure parallel to the

CE

N-trending Sagaing fault. This proposition is supported by the analogous spatial association of continental slivers of Gondwanan affinity to the east of the HP accretionary complexes in

AC

both the locations: continental Naga Metamorphics of Lower Palaeozoic age (Aitchison et al., 2019) in the NOC and Mogok Metamorphic and Slate Belt as part of the Sibumasu Continental block in northern Myanmar (cf. Metcalfe, 1995; Zhang, J. et al., 2018) (Fig. 17c). Future studies should be directed to understand the Neo-Tethyan tectonics that was responsible for the present association of Jurassic and Cretaceous accretionary complexes in the WOB.

ACCEPTED MANUSCRIPT 7. Conclusions 1. Oceanic eclogite from two new locations near Mokie and Thewati in the Nagaland Ophiolite Complex, NE India, occurs as meter-scale lensoidal bodies within a lawsoniteblueschist facies metamorphosed oceanic sequence of basalt-chert-limestone. 2. SHRIMP U-Pb zircon dating of the eclogites yielded age groupings at ~205 Ma,

SC RI PT

~185-189 Ma and ~172 Ma (minor). The c. 205 Ma date is interpreted as a minimum age for the oceanic protolith, the c.185-189 Ma age is thought to reflect zircon growth at eclogite facies conditions and the c.172 Ma age possibly represents the onset of retrograde cooling and exhumation.

NU

3. Maximum P-T conditions recorded by the eclogite samples (~28.1 to 25.5 kbar and ~ 650 oC for the Thewati eclogite and ~23.8 kbar and ~555 oC for the Mokie eclogite)

MA

involved a low apparent thermal gradient (~7-8 oC/km).

4. The Thewati eclogite was buried to depths of ~80-90 km, transitional to UHP

ED

conditions, thus falling in the select category of oceanic eclogites worldwide that are recovered from the maximum depths in the subduction zone.

PT

5. The Thewati eclogite records a near complete clockwise metamorphic P-T path of

CE

evolution that involved an epidote blueschist facies prograde stage, peak metamorphism in the epidote eclogite facies and multistage retrogression through amphibole-eclogite facies at

AC

~18.3 kbar, 630 oC, epidote blueschist facies to the transitional lawsonite blueschist and greenschist facies metamorphic conditions at ~6 kbar, 300 oC. 6. The low thermal gradients (~7-9 oC/km) that accompanied peak conditions, together with the clockwise P-T paths of the blueschists and eclogites are consistent with a Lower Jurassic-aged, cold intra-oceanic subduction system within the Neo-Tethys in Nagaland.

ACCEPTED MANUSCRIPT 7. The Nagaland eclogites were exhumed in two stages: (a) the initial exhumation (from 25.5 kbar, 650 oC to 10 kbar, 510 oC, and along a steep dP/dT gradient) was likely to have taken place in the subduction channel as a consequence of buoyancy-driven material transport in a rheologically weak and fluidised mantle; (b) in the latter phase, the exhumation of the eclogites along a gentler dP/dT gradient was accompanied by extended cooling at

SC RI PT

shallower levels. This exhumation stage is linked with thrust stacking of the partially exhumed eclogites with the cooler, prograde blueschists at shallower crustal levels. 8. The current disposition of the Nagaland eclogites in lower grade blueschist likely reflects stage 2 of their exhumation in the subduction channel, when the Nagaland

NU

accretionary complex was assembled.

10. An Early Jurassic accretionary complex in the NOC was likely continuous with

current geographic co-ordinates).

MA

the accretionary complex of the Jade Mines Belt of northern Myanmar (with respect to

ED

11. The Nagaland eclogite samples present the oldest Neo-Tethyan remnant (cf. Late Triassic age) and evidence of the earliest subduction (cf. Early Jurassic) within the Yarlung-

PT

Tsangpo Suture Zone. The cryptic record of this pre-Cretaceous subduction system within the

CE

Neo-Tethys that incorporated an Early-Middle Jurassic Andean-type convergent margin can

Pakistan.

AC

be traced westwards from the Indo-Myanmar Ranges along the southern Eurasian margin to

Acknowledgements The paper constitutes a part of the doctoral thesis work of MR for which she acknowledges funding from the Indian Institute of Technology Kharagpur (IIT KGP) for an Institute doctoral fellowship and to carry out geological fieldwork. S.K.B. acknowledges financial support from the Indian Space Research Organization (Grant No. IIT/KCSTC/Chairman/New

ACCEPTED MANUSCRIPT Approval/15-16/09) and Australia-India Strategic Research Fund (DST/INT/AUS/P-55/2013 and AISRF07021) for fieldwork and analyses. A.A. thanks the Director, Wadia Institute of Himalayan Geology. We would also like to acknowledge the generous support of the Nagaland State Mineral Development Corporation Ltd. and their support staff at the Weziho Field Hostel without which, it wouldn't have been possible for us to undertake geological

SC RI PT

fieldwork in this remote and inaccessible part of the Nagaland state bordering Myanmar. The manuscript has benefitted from earlier reviews by S.A. Wilde and T. Johnson. We thank two anonymous journal reviewers for detailed reviews, and Marco Scambelluri for efficient

NU

editorial handling.

References

MA

Acharyya, S.K., 1986. Tectono-stratigraphic history of Naga Hills Ophiolites. In: Geology of the Nagaland Ophiolite (ed. Ghosh, D.B.), Geological Survey of India Memoirs

ED

119,94–103.

Acharyya, S.K., 2007. Collisional emplacement history of the Naga-Andaman ophiolites and

PT

the position of the eastern Indian suture. Journal of Asian Earth Sciences 29, 229-242.

CE

Acharyya, S.K., 2015. Indo-Burma Range: a belt of accreted microcontinents, ophiolites and Mesozoic–Paleogene flyschoid sediments. International Journal of Earth Sciences 104,

AC

1235-1251.

Agard, P., Jolivet, L., Goffé, B., 2001. Tectonometamorphic evolution of the Schistes Lustrés complex: implications for the exhumation of HP and UHP rocks in the Western Alps. Bulletin de la Société Géologique de France 172 (5), 617–636. Agard, P., Yamato, P., Jolivet, L., Burov, J.E., 2009. Exhumation of oceanic blueschists and eclogites in subduction zones: Timing and mechanisms. Earth Science Reviews 92,5379.

ACCEPTED MANUSCRIPT Agard, P., Yamato, P., Soret, M., Prigent, C., Guillot, S., Plunder, A., Dubacq, B., Chauvet,A., Monié, P., 2016. Plate interface rheological switches during subduction infancy: control on slab penetration and metamorphic sole formation. Earth and Planetary Science Letters 451, 208–220. Agard, P., Plunder, A., Angiboust, S., Bonnet, G., Ruh, J., 2018. The subduction plate

SC RI PT

interface: Rock record and mechanical coupling (from long to short time scales). Lithos 320-321, 537-566. https://doi.org/10.1016/j.lithos.2018.09.029.

Aitchison, J.C., McDermid, I.R.C., Ali, J.R., Davis, A.M., Zyabrev, S.V., 2007. Shoshonites in southern Tibet record Late Jurassic rifting of a Tethyan intra-oceanic island arc.

NU

Journal of Geology, 115: 197-213.

Aitchison, J.C., Ao, A., Bhowmik, S., Clarke, G.L., Ireland, T.R., Kachovich, S., Lokho,

MA

K.,Stojanovic, D., Roeder, T., Truscott, N., Zhen, Y., Zhou, R., 2019. Tectonic evolution of the western margin of the Myanmar microplate based on new fossil and

ED

radiometric age constraints. Tectonics 38, 1718–1741. https://doi.org/10.1029/2018TC005049.

PT

Anczkiewicz, R., Burg, J.P., Villa, I.M., Meier, M., 2000. Late Cretaceous blueschist

CE

metamorphism in the Indus suture zone, Shangla region, Pakistan Himalaya. Tectonophysics 324, 111-134.

AC

Andersen, T. B., Jamtveit, B., Dewey, J. F., Swensson, E., 1991. Subduction and eduction of continental crust: major mechanisms during continent‐continent collision and orogenic extensional collapse, a model based on the south Norwegian Caledonides. Terra Nova 3(3), 303-310. Angiboust, S., Agard, P., Jolivet, L., Beyssac, O., 2009. The Zermatt‐Saas ophiolite: the largest (60km wide) and deepest (c. 70–80 km) continuous slice of oceanic lithosphere detached from a subduction zone? Terra Nova 21(3), 171-180.

ACCEPTED MANUSCRIPT Angiboust, S., Langdon, R., Agard, P., Waters, D., Chopin, C., 2012. Eclogitization of the Monviso ophiolite (W. Alps) and implications on subduction dynamics. Journal of Metamorphic Geology 30, 37-61. Angiboust, S., Agard, P., Glodny, J., Omrani, J., Oncken, O., 2016. Zagros blueschists: Episodic underplating and long-lived cooling of a subduction zone. Earth and Planetary

SC RI PT

Science Letters 443, 48-58. https://doi.org/10.1016/j.epsl.2016.03.017. Anon., 1986. Geology of Nagaland ophiolite. Geological Survey of India Memoir, 119, 113pp.

Ao, A., Bhowmik, S.K., 2014. Cold subduction of the Neotethys: The metamorphic record

NU

from finely banded lawsonite and epidote blueschists and associated metabasalts of the Nagaland Ophiolite Complex, India. Journal of Metamorphic Geology 32(8), 829-860.

MA

Baxter, A.T., Aitchison, J.C., Zyabrev, S.V., Ali, J.R., 2011. Upper Jurassic radiolarians from the Naga ophiolite, Nagaland, northeast India. Gondwana Research 20, 638-644.

ED

Bhowmik, S.K., Ao, A., 2016. Subduction initiation in the Neo‐Tethys: Constraints from counterclockwise P–T paths in amphibolite rocks of the Nagaland Ophiolite Complex,

PT

India. Journal of Metamorphic Geology 34(1), 17-44.

CE

Black, L.P., Kamo, S.L., Allen, C.M., Aleinikoff, J.N., 2003. TEMORA 1: a new zircon standard for Phanerozoic U-Pb geochronology. Chemical Geology 200, 155–17.

AC

Bohlen, S.R., Boettcher, A.L., 1982. The quartz–coesite transformation: a pressure determination and the effects of other components. Journal of Geophysical Research104,208–224. Brovarone, A.V., Agard, P., 2013. True metamorphic isograds or tectonically sliced metamorphic sequence? New high-spatial resolution petrological data for the New Caledonia case study. Contributions to Mineralogy and Petrology 166(2), 451-469.

ACCEPTED MANUSCRIPT Brovarone, A. V., Groppo, C., Hetenyi, G., Compagnoni, R., Malavieille, J., 2011. Coexistence of lawsonite‐bearing eclogite and blueschist: phase equilibria modelling of Alpine Corsica metabasalts and petrological evolution of subducting slabs. Journal of Metamorphic Geology 29(5), 583-600. Brovarone, A. V., Picatto, M., Beyssac, O., Lagabrielle, Y., Castelli, D. 2014. The

SC RI PT

blueschist–eclogite transition in the Alpine chain: P–T paths and the role of slowspreading extensional structures in the evolution of HP–LT mountain belts.Tectonophysics 615, 96-121.

Brunnschweiler, R.O., 1966. On the geology of Indo-Burman Ranges. Geological Society of

NU

Australia 13, 137–195.

Chatterjee, N., Ghose, N.C., 2010. Metamorphic evolution of the Naga Hills eclogite and

MA

blueschist, Northeast India: Implications for early subduction of the Indian plate under the Burma microplate. Journal of Metamorphic Geology 28(2), 209-225.

ED

Chemenda, A. I., Mattauer, M., Bokun, A. N., 1996. Continental subduction and a mechanism for exhumation of high-pressure metamorphic rocks: new modelling and

PT

field data from Oman. Earth and Planetary Science Letters 143(1-4), 173-182.

CE

Clarke, G. L., Aitchison, J. C., Cluzel, D., 1997. Eclogites and blueschists of the Pam Peninsula, NE New Caledonia: a reappraisal. Journal of Petrology 38(7), 843-876.

AC

Cloos, M., 1982. Flow melanges: numerical modelling and geologic constraints on their origin in the Franciscan subduction complex, California. Bulletin of Geological Society of America 93, 330–345. Cloos, M., 1985. Thermal evolution of convergent plate margins: Thermal modeling and reevaluation of isotopic Ar‐ages for blueschists in the Franciscan Complex of California. Tectonics 4(5), 421-433.

ACCEPTED MANUSCRIPT Cloos, M., Shreve, R., 1988. Subduction channel mode of prism accretion, mélange formation, sediment subduction, and subduction erosion at convergent plate margins: 2. Implications and discussions. Pure and Applied Geophysics 128, 501–541. Coggon, R., Holland, T.J.B., 2002. Mixing properties of phengitic micas and revised garnet‐phengite thermobarometers. Journal of Metamorphic Geology 20(7), 683-696.

SC RI PT

Davis, P.B. Whitney, D.L., 2008. Petrogenesis and structural petrology of high-pressure metabasalt pods, Sivrihisar, Turkey. Contributions to Mineralogy and Petrology 156, 217– 241.

Diener, J.F.A., Powell, R., White, R.W. Holland, T.J.B., 2007. A new thermodynamic model

NU

for clino- and orthoamphiboles in the system Na2O-CaO-FeO-MgO-Al2O3-SiO2-H2OO. Journal of Metamorphic Geology 25, 631–656.

MA

Fitzherbert, J. A., Clarke, G. L., Marmo, B.,Powell, R., 2004. The origin and P–T evolution of peridotites and serpentinites of NE New Caledonia: prograde interaction between

ED

continental margin and the mantle wedge. Journal of Metamorphic Geology 22(4), 327344. https://doi.org/10.1111/j.1525-1314.2004.00517.x.

PT

Ganguly, J., Cheng, W., Tirone, M., 1996. Thermodynamics of aluminosilicate garnet solid

CE

solution: new experimental data, an optimized model, and thermometric applications. Contributions to Mineralogy and Petrology 126, 137–151.

AC

Gardiner, N.J., Searle, M.P., Morley, C.K., Robb, L.J., Whitehouse, M.J., Roberts, N.M., Kirkland, C.L., Spencer, C.J., 2018. The crustal architecture of Myanmar imaged through zircon U-Pb, Lu-Hf and O isotopes: Tectonic and metallogenic implications. Gondwana Research 62, 27-60. doi:10.1016/j.gr.2018.02.008. Gerya, T.V., Stoeckhert, B., Perchuk, A.L., 2002. Exhumation of high-pressure metamorphic rocks in a subduction channel: a numerical simulation. Tectonics 21, 61-69, 1056. doi:10.1029/2002TC001406.

ACCEPTED MANUSCRIPT Goffé, B., Rangin, C., Maluski, H., 2000. Jade and associated rocks from Jade Mines area, northern Myanmar, as record of a polyphased high-pressure metamorphism. EOS Transactions AGU 81, Fall Meeting Supplement, F1365. Goffé, B., Rangin, C., Maluski, H., 2002. Jade and associated rocks from the Jade Mines area, northern Myanmar as record of a polyphased high-pressure metamorphism.

SC RI PT

Himalaya Karakoram-Tibet Workshop Meeting (Abstracts). Journal of Asian Earth Sciences 20, 16–17.

Green, E.C.R., Holland, T.J.B., Powell, R., 2007. An order disorder model for omphacitic pyroxenes in the system jadeite- diopside-hedenbergite-acmite, with applications to

NU

eclogite rocks. American Mineralogist 92, 1181–1189.

Groppo, C., Castelli, D., 2010. Prograde P-T evolution of a lawsonite eclogite from the

MA

Monviso metaophiolite (Western Alps): dehydration and redox reactions during subduction of oceanic FeTi-oxide gabbro. Journal of Petrology 51, 2489–2514.

ED

Groppo, C., Rolfo, F., Sachan, H.K., Rai, S.K., 2016. Petrology of blueschist from the Western Himalaya (Ladakh, NW India): Exploring the complex behavior of a

PT

lawsonite-bearing system in a paleo-accretionary setting. Lithos 252, 41-56.

CE

Guillot, S., Hattori, K., Agard, P., Schwartz, S., Vidal, O., 2009. Exhumation processes in oceanic and continental subduction contexts: a review. In: Lallemand, S., Funcinello, F.

AC

(Eds.), Subduction Zone Geodynamics: pp. 175–205. http://dx.doi.org/10.1007/978-3540-87974-9.

Guilmette, C., Hébert, R., Dostal, J., Indares, A., Ullrich, T., Bédard, É., Wang, C., 2012. Discovery of a dismembered metamorphic sole in the Saga ophiolitic mélange, South Tibet: Assessing an Early Cretaceous disruption of the Neo-Tethyan supra-subduction zone and consequences on basin closing. Gondwana Research 22(2), 398-414.

ACCEPTED MANUSCRIPT Guiraud, M., 1982. Géothermométrie du faciés schistes vert à glaucophane. Modélisation et applications (Afghanistan, Pakistan, Corse, Bohème). Ph.D. thesis, University of Montpellier, France. Hacker, B.R., Wang, X., Eide, E.A., Ratsbacher, L., 1996. The Qinling‐Dabie ultrahigh‐pressure collisional orogen. In: Tectonic Evolution of Asia (Yin, A.,

SC RI PT

Harrison, T.M. eds), Cambridge University Press, Cambridge, pp. 345-370. Hébert, R., Bezard, R., Guilmette, C., Dostal, J., Wang, C.S., Liu, Z.F., 2012. The Indus– Yarlung Zangbo ophiolites from Nanga Parbat to Namche Barwa syntaxes, southern Tibet: First synthesis of petrology, geochemistry, and geochronology with incidences

NU

on geodynamic reconstructions of Neo-Tethys. Gondwana Research 22, 377-397. Hermann, J., Müntener, O., Scambelluri, M., 2000. The importance of serpentinite mylonites

MA

for subduction and exhumation of oceanic crust. Tectonophysics 327(3-4), 225-238. Holland, T.J.B., 1980. The reaction albite = jadeite+ quartz determined experimentally in the

ED

range 600–1200 oC. American Mineralogist 65(1-2), 129-134. Holland, T.J.B., Powell, R., 1998. An internally consistent thermodynamic dataset for phases

PT

of petrological interest. Journal of Metamorphic Geology 16, 309–343.

CE

Holland, T.J.B., Powell, R., 2003. Activity–composition relations for phases in petrological calculations: an asymmetric multicomponent formulation. Contributions to Mineralogy

AC

and Petrology 145, 492–501. Honegger, K., Le Fort, P., Mascle, G., Zimmerman, J.L., 1989. The blueschists along the Indus Suture Zone in Ladakh, N.W. Himalaya. Journal of Metamorphic Geology 7, 57– 72. Htay, H., Zaw, K., Oo, T.T., 2017. The mafic–ultramafic (ophiolitic) rocks of Myanmar. In: Barber, A.J., Khin, Z., & Crow, M.J. (eds) Myanmar: Geology, Resources and Tectonics. Geological Society, London, Memoirs 48, 117–141.

ACCEPTED MANUSCRIPT Imchen, W., Thong, G.T., Pongen, T., 2014. Provenance, tectonic setting and age of the sediments of the Upper Disang Formation in the Phek District, Nagaland. Journal of Asian Earth Sciences 88, 11–27. Jan, M.Q., 1985. High-P rocks along the suture zone around Indo-Pakistan plate and phase chemistry of blueschists from eastern Ladakh. Geological Bulletin University of

SC RI PT

Peshawar 18, 1–40.

Krebs, M., Schertl, H.P., Maresch, W.V., Draper, G., 2011. Mass flow in serpentinite-hosted subduction channels: P-T-t path patterns of metamorphic blocks in the Rio San Juan melange (Dominican Republic). Journal of Asian Earth Sciences 42, 569–595.

NU

Lanari, P., Vidal, O., De Andrade, V., Dubacq, B., Lewin, E., Grosch, E.G., Schwartz, S.,2014. XMapTools: A MATLAB©-based program for electron microprobe X-ray

MA

image processing and geothermobarometry. Computers & Geosciences 62, 227-240. Leake, B.E., Wooley, A.R., Arps, C.E., Birch, W.D., Gilbert, M.C., Grice, J.D., Hawthorne,

ED

F.C., Kato, A., Kisch, H.J., Krivovichec, V.G. and Linthout, K., 1997. Nomenclature of amphiboles: report of the subcommittee on amphiboles of the International

PT

Mineralogical Association, commission on new minerals and mineral names. American

CE

Mineralogist 82, 1019–1037.

Liou, J.G., 1971. P-T Stabilities of laumontite, wairakite, lawsonite, and related minerals in

AC

the system CaAl2Si2O8-SiO2-H2O. Journal of Petrology 12(2), 379-411. Liu, C.Z., Chung, S.L., Wu, F.Y., Zhang, C., Xu, Y., Wang, J.G., Chen, Y., Guo, S., 2016a. Tethyan suturing in Southeast Asia: zircon U-Pb and Hf-O isotopic constraints from Myanmar ophiolites. Geology 44, 311–314. Liu, C.Z., Zhang, C., Xu, Y., Wang, J.G., Chen, Y., Guo, S., Wu, F.Y., Sein, K., 2016b. Petrology and geochemistry of mantle peridotites from the Kalaymyo and Myitkyina ophiolites (Myanmar): implications for tectonic settings. Lithos 264, 495–508.

ACCEPTED MANUSCRIPT Ludwig, K.R., 2003. User’s manual for a geochronological toolkit for Microsoft Excel (Isoplot/Ex version 3.0). Berkeley Geochronology Center, Special Publication 4, 1-71. Mahéo, G., Fayoux, X., Guillot, S., Garzanti, E., Capiez, P., Mascle, G., 2006. Relicts of an intra-oceanic arc in the Sapi-Shergol mélange zone (Ladakh, NW Himalaya, India): implications for the closure of the Neo-Tethys Ocean. Journal of Asian Earth Sciences

SC RI PT

26, 695–707.

Maruyama, S., Cho, M., Liou, J. G., 1986. Experimental investigations of blueschistgreenschist transition equilibria: Pressure dependence of Al2O3 contents in sodic amphiboles- A new geobarometer. Blueschists and Eclogites 164, 1-16.

NU

McDermid, I.R.C., Aitchison, J.C., Davis, A.M., Harrison, T.M., Grove, M., 2002. The Zedong terrane: A Late Jurassic intra-oceanic magmatic arc within the Yarlung-Zangbo

MA

suture zone, southeastern Tibet. Chemical Geology 187, 267-277. Metcalfe, I., 1995. Gondwana dispersion and Asian accretion. Journal of Geology 5-6, 223–

ED

266.

Mével, C., Kiénast, J.R., 1986. Jadeite-kosmochlor solid solution and chromian sodic

PT

amphiboles in jadeitites and associated rocks from Tawmaw (Burma). Bulletin de

CE

Mineralogie 109, 617–633.

Mitchell, A.H.G., 1986. Ophiolite and associated rocks in four settings: relationships to

AC

subduction and collision. Tectonophysics 125, 269-285. Mitchell, A.H.G., 1993. Cretaceous-Cenozoic Tectonic events in the western Myanmar (Burma)-Assam region. Journal of Geological Society, London 150, 1089–1102. Mitchell, A.H.G., Htay, M.T., Htun, K.M., 2015. The medial Myanmar suture zone and the western Myanmar Mogok foreland. Journal of the Myanmar Geosciences Society 6,73– 88.

ACCEPTED MANUSCRIPT Nyunt, T.T., Massonne, H.J., Sun, T.T., 2017. Jadeitite and other high-pressure metamorphic rocks from the Jade Mines Belt, Tawmaw area, Kachin State, northern Myanmar. Geological Society, London, Memoirs 48(1), 295-315. Oberhhänsli, R., 2013. High-pressure – low-temperature evolution in the Indus-Tsangpo Suture along the Kohistan Arc (Kaghan Valley, NE Pakistan). Episodes 36, 87–93.

SC RI PT

Oberhänsli, R., Bousquet, R., Moinzadeh, H., Moazzen, M., Arvin, M., 2007. The field of stability of blue jadeite: a new occurrence of jadeitite from Sorkhan, Iran, as a case study. Canadian Mineralogist 45, 1501–1509.

O'Brien, P.J., 2018. Eclogites and other high-pressure rocks in the Himalaya: a review.

NU

Geological Society, London, Special Publications 483,https://doi.org/10.1144/SP483.13 Okamoto, K., Maruyama, S., 1999. The high-pressure synthesis of lawsonite in the

MA

MORB+H2O system. American Mineralogist 84(3), 362-373. Page, F.Z., Armstrong, L.S., Essene, E.J., Mukasa, S.B., 2007. Prograde and retrograde

ED

history of the Junction School eclogite, California, and an evaluation of garnet-

153, 533–555.

PT

phengite-clinopyroxene thermobarometry. Contributions to Mineralogy and Petrology

CE

Peacock, S.M., Wang, K., 1999. Seismic consequences of warm versus cool subduction metamorphism: examples from southwest and northeast Japan. Science 286, 937–939.

AC

Plunder, A., Agard, P., Dubacq, B., Chopin, C., Bellanger, M., 2012. How continuous and precise is the record of P–T paths? Insights from combined thermobarometry and thermodynamic modelling into subduction dynamics (Schistes Lustrés, W. Alps). Journal of Metamorphic Geology 30(3), 323-346. Plunder, A., Agard, P., Chopin, C., Pourteau, A., Okay, A. I., 2015. Accretion, underplating and exhumation along a subduction interface: From subduction initiation to continental subduction (Tavşanlı zone, W. Turkey). Lithos 226, 233-254.

ACCEPTED MANUSCRIPT Powell, R., Holland, T., 2008. On thermobarometry. Journal of Metamorphic Geology 26, 155–179. Qiu, Z.L., Wu, F.Y., Yang, S.F., Zhu, M., Sun, J.F., Yang, P., 2009. Age and genesis of the Myanmar jadeite: constraints from U-Pb ages and Hf isotopes of zircon inclusions. Chinese Science Bulletin 54, 658–668.

SC RI PT

Ravna, K., 2000. The garnet–clinopyroxene Fe2+–Mg geothermometer: an updated calibration. Journal of Metamorphic Geology 18, 211-219.

Ravna, E.J., Terry, M.P., 2004. Geothermobarometry of UHP and HP eclogites and schists– an evaluation of equilibria among garnet–clinopyroxene–kyanite–phengite–

NU

coesite/quartz. Journal of Metamorphic Geology 22, 57-592.

Rebay, G., Powell, R., Diener, J.F.A., 2010. Calculated phase equilibria for a MORB

MA

composition in a P–T range, 450–650°C and 18–28 kbar: the stability of eclogite. Journal of Metamorphic Geology 28, 635–645.

ED

Reinecke, T., 1998. Prograde high- to ultrahigh-pressure metamorphism and exhumation of oceanic sediments at Lago di Cignana, Zermatt-Saas Zone, western Alps. Lithos 42,

PT

147–189.

CE

Roddick, J.C., van Breemen, O., 1994. U-Pb zircon dating: A comparison of ion microprobe and single grain conventional analyses. In Radiogenic Age and Isotopic Studies,

AC

Report8. Geological Survey of Canada Current Res 1994-F, p 1-9. Rubatto, D., 2017. Zircon: the metamorphic mineral. Reviews in Mineralogy and Geochemistry 83, 261-295. Rubatto, D., Gebauer, D., 2000. Use of cathodoluminescence for U-Pb zircon dating by ion microprobe: some examples from the Western Alps. In Cathodoluminescence in geosciences. Springer, Berlin, Heidelberg. pp. 373-400.

ACCEPTED MANUSCRIPT Rubatto, D., Hermann, J., 2003. Zircon formation during fluid circulation in eclogites (Monviso, Western Alps): implications for Zr and Hf budget in subduction zones. Geochimica et Cosmochimica Acta 67, 2173-2187. Rubatto, D., Hermann, J., 2007. Zircon behaviour in deeply subducted rocks. Elements 3, 3135.

SC RI PT

Sambridge, M.S., Compston, W., 1994. Mixture modeling of multi-component data sets with application to ion-probe zircon ages. Earth and Planetary Science Letters 128, 373–390. Sarkar, A., Datta, A.K., Poddar, B.C., Bhattacharyya, B.K., Kollapuri, V.K., Sanwal, R.,1996. Geochronological studies of Mesozoic igneous rocks from eastern India.

NU

Journal of Southeast Asian Earth Sciences 13, 77–81.

Schmidt, M.W., 1995. Lawsonite: upper pressure stability and formation of higher density

MA

hydrous phases. American Mineralogist 80 (11-12), 1286-1292. Searle, M.P., Kahn, M.A., Fraser, J.E., Gough, S.J., 1999. The tectonic evolution of the

ED

Kohistan-Karakoram collision belt along the Karakoram Highway transect, north Pakistan. Tectonics 18, 929–949.

PT

Searle, M.P., Noble, S.R., Cottle, J.M. et al., 2007. Tectonic evolution of the Mogok

CE

Metamorphic belt, Burma (Myanmar) constrained by U-Th-Pb dating of metamorphic and magmatic rocks. Tectonics 26, TC3014. https://doi:10.1029/2006TC002083.

AC

Searle, M.P., Morley, C.K., Waters, D.J., Gardiner, N.J., Htun, U. K., Nu, T.T., Robb, L.J., 2017.Tectonic and metamorphic evolution of the Mogok Metamorphic and Jade Mines belts and ophiolitic terranes of Burma (Myanmar). In: Barber, A.J., Khin, Z., & Crow, M.J. (eds) Myanmar: Geology, Resources and Tectonics. Geological Society, London, Memoirs 48, 261–293, https://doi.org/10.1144/M48.12.

ACCEPTED MANUSCRIPT Setiawan, N.I., Osanai, Y., Nakano, N., Adachi, T., Yonemura, K., Yoshimoto, A., 2016. Prograde and retrograde evolution of eclogites from the Bantimala Complex in South Sulawesi, Indonesia. Journal of Mineralogical and Petrological Sciences 111, 211–225. Shi, G.H., Cui, W.Y., Tropper, P., Wang, C.Q., Shu, G.M., Yu, H.X., 2003. The petrology of a complex sodic and sodic-calcic amphibole association and its implications for the

SC RI PT

metasomatic processes in the jadeitite area in northwestern Myanmar, formerly Burma. Contributions to Mineralogy and Petrology 145, 355–376.

Shi, G.H., Cui, W.Y., Cao, S.M., Jiang, N., Jian, P., Liu, D.Y., Miao, L.C., Chu, B.B., 2008. Ion microprobe zircon U–Pb age and geochemistry of the Myanmar jadeitite. Journal of

NU

the Geological Society, London 165, 221–234.

Singh, A.K., Chung, S.-L., Bikramaditya, R.K., Lee, H.Y., 2017. New U–Pb zircon ages of

MA

plagiogranites from the Nagaland-Manipur Ophiolites, Indo-Myanmar Orogenic Belt, NE India. Journal of the Geological Society, London 74, 170-179.

ED

Sorensen, S. S., Barton, M. D., 1987. Metasomatism and partial melting in a subduction complex: Catalina Schist, southern California. Geology 15(2), 115-118.

PT

Thet Tin Nyunt., 2009. Petrological and geochemical contribution to the origin of jadeitite

CE

and associated rocks of the Tawmaw area, Kachin State, Myanmar. Doctoral dissertation, Institut für Mineralogie und Kristallchemie der Universität Stuttgart.

AC

Vance, D., Mahar, E., 1998. Pressure-temperature paths from P-T pseudosections and zoned garnets: potential, limitations and examples from the Zanskar Himalaya, NW India. Contributions to Mineralogy and Petrology 132(3), 225-245. Von Blanckenburg, F., Davies, J. H., 1995. Slab break off: a model for syncollisional magmatism and tectonics in the Alps. Tectonics 14(1), 120-131.

ACCEPTED MANUSCRIPT Wang, W.L., Aitchison, J.C., Lo, C.H., Zeng, Q.G., 2008. Geochemistry and geochronology of the amphibolite blocks in ophiolitic mélanges along Bangong-Nujiang suture, central Tibet. Journal of Asian Earth Sciences 33, 122-138. Wei, C.J., Clarke, G.L., 2011. Calculated phase equilibria for MORB compositions: a reappraisal of the metamorphic evolution of lawsonite eclogite. Journal of

SC RI PT

Metamorphic Geology 29 (9), 939-952.

White, R.W., Powell, R., Holland, T.J.B., 2007. Progress relating to calculation of partial melting equilibria for metapelites. Journal of Metamorphic Geology 25, 511–527. Whitney, D.L., Evans, B.W., 2010. Abbreviations for names of rock-forming minerals.

NU

American Mineralogist 95, 185-187.

Xu, Y., Liu, C.Z., Guo, S., Wang, J., Sein, K., 2017. Petrogenesis and tectonic implications

MA

of gabbro and plagiogranite intrusions in mantle peridotites of the Myitkyina ophiolite, Myanmar. Lithos 284-285, 180–193.

ED

Yang, J.S., Xu, Z.Q., Duan, X.D., Li, J., Xiong, F.H., Liu, Z., Cai, Z.H., Li, H.Q., 2012. Discovery of a Jurassic SSZ ophiolite in the Myitkyina region of Myanmar. Acta

PT

Petrologica Sinica 28, 1710–1730 (in Chinese with English abstract).

CE

Yamato, P., Agard, P., Burov, E., Le Pourhiet, L., Jolivet, L., Tiberi, C., 2007. Burial and exhumation in a subduction wedge: mutual constraints from thermo-mechanical

AC

modelling and natural P-T-t data (Schistes Lustres, W. Alps). Journal of Geophysical Research 112, B07410. doi:10.1029/2006JB004441. Yin, A., Harrison, T.M., 2000. Geologic evolution of the Himalayan Tibetan orogen. Annual Review of Earth and Planetary Sciences 28, 211–80. Yu, Y-P., Xie, C-M., Fan, J-J., Wang, M., Dong, Y-C, Wang, B., Hao, Y-J., 2018. Zircon U– Pb geochronology and geochemistry of Early Jurassic granodiorites in Sumdo area,

ACCEPTED MANUSCRIPT Tibet: Constraints on petrogenesis and the evolution of the Neo-Tethyan Ocean. Lithos 320-321, 134-143. Yui, T.F., Fukoyama, M., Iizuka, Y. et al., 2013. Is Myanmar jadeitite of Jurassic age? A result from incompletely recrystallized inherited zircon. Lithos 160–161, 268–282. Zhang, J., Xiao, W., Windley, B.F., Cai, F., Sein, K., Naing, S., 2017. Early Cretaceous

SC RI PT

wedge extrusion in the Indo-Burma Range accretionary complex: Implications for the Mesozoic subduction of Neotethys in SE Asia. International Journal of Earth Sciences 106, 1391–1408.

Zhang, J., Xiao, W., Windley, B., Wakabayashi, J., Cai, F., Sei, K., Wu, H., Naing, S., 2018.

NU

Multiple alternating forearc- and backarc-ward migration of magmatism in the IndoMyanmar Orogenic Belt since the Jurassic: Documentation of the orogenic architecture

MA

of eastern Neotethys in SE Asia. Earth Science Reviews 185, 704-731. Zhang, Y.X., Jin, X., Zhang, K.J., Sun, W.D., Liu, J.M., Zhou, X.Y., Yan, L.L., 2018. Newly

ED

discovered Late Triassic Baqing eclogite in central Tibet indicates an anticlockwise West–East Qiangtang collision. Scientific Reports. 8(1), 966.

PT

https://doi:10.1038/s41598-018-19342-w.

CE

Zhu, D-C., Zhao, Z-D., Niu, Y., Dilek, Y., Hou, Z-Q., Mo, X-X., 2013. The origin and pre-

AC

Cenozoic evolution of the Tibetan Plateau. Gondwana Research 23, 1429-1454.

Figure Captions

Fig. 1: A tectonic map of eastern Tibet, SE China, Myanmar and NE India showing the occurrences of different Tethyan sutures (LCS, Longmu Co-Shuanghu suture; BNS, Bangong-Nujiang suture; YTS, Yarlung-Tsangpo suture (after Liu et al., 2016a) and locations of oceanic blueschists and amphibolites (1-4) along the YTS zone. N/M: Nagaland/Myktyina ophiolites. Other abbreviations used: IMR, Indo-Myanmar Ranges;

ACCEPTED MANUSCRIPT WMB, West Myanmar Block; MCT, Main Central thrust; STDS, South Tibetan Detachment System. Fig. 2: (a) Locations of the Indo-Myanmar Ranges (IMR) and western, central and eastern ophiolite belts (WOB, COB and EOB) in the tectonic framework of Indian, Asian and Burmese plates (modified after Searle et al., 2007). The outlined box in the northern margin

SC RI PT

of the IMR marks the location of the present study in the Nagaland Ophiolite Complex (NOC), shown in detail in Fig. 2b. Geochronological data sources for magmatic and metamorphic events are superscripted and subscripted as follows: Superscript:1. Shi et al.(2008); 2. Yui et al.(2013); 3. Goffé et al. (2002); 4. Singh et al. (2017); 5. Liu et al.

NU

(2016a); 6. Yang et al. (2012); 7. Zhang et al. (2017). Subscript: a, Jadeitite; b, Eclogite; c, Plagiogranite; d, Diorite and Gabbro; e, Rodingites; f, Amphibolite; g, Andesitic basalt,

MA

Leucogabbro and Plagiogranite; h, Trondhjemite. (b) Geological map of the NOC (modified after Anon., 1986; Ao and Bhowmik, 2014). Also shown are the distributions of different

ED

metamorphic facies and locations of eclogite samples [(Mokie eclogite, N22F and Thewati eclogites, N23B and N40B, this study and of previous study, abbreviated as CG10 (after

PT

Chatterjee & Ghose, 2010), along with their peak P-T estimates], blueschist occurrences (BS-

CE

A86, after Anon., 1986) and dated basalt sample (S-96, after Sarkar et al., 1996).

Eclogite.

AC

Abbreviations used: GS, Greenschist; PD, Pumpellyite-diopside; BS, Blueschist, EC,

Fig. 3: (a-b) Outcrop-scale maps showing the occurrences of Mokie (Sample N22F) (a) and Thewati (Samples N23B and N40B) (b) eclogites in a blueschist facies ensemble of metabasalt, marble and metachert. Photographs of hand specimens (c, e, g), photomicrograph (f) and false-colour thin section map (d) of Thewati(c-f) and Mokie (g) eclogite samples showing the mineralogy, texture and micro-structure of these rocks. See text for

ACCEPTED MANUSCRIPT more details. Mineral abbreviations in this and other figures and tables are after Whitney & Evans (2010). Fig. 4: Back-scattered electron images (BSEI, a-b, e, h, n-q, s-u) and photo-micrographs (c-d, f-g, i-m, r) of Thewati (a-q) and Mokie eclogites showing textural features of prograde, peak and multi-stage retrograde metamorphism. (a-c): Porphyroblastic Grt (a-b) and omphacitic

SC RI PT

Cpx (c) contain prograde inclusions of Ep1, Gln1 and Ph. (d) Relics of Omp and Rt occur within Brs porphyroblast. Brs is rimmed by a thin corona of Gln2. (e-f) The Thewati eclogite sample N40B shows a second textural variety of Grt (cf. Grt-B) (e-f), which is smaller in size compared to the larger Grt porphyroblasts (cf. Grt-A) (f). Grt-A shows inclusion-rich core

NU

and inclusion-poor rim. (g-h) The occurrence of Ep2 as very coarse aggregates in layer parallel segregations (g) or in the matrix foliation domain (h). Ep2 contains relics of Grt and

MA

Omp (g-h). (i-l) Syn- (i) to post-tectonic (j-l) growth of Brs. Brs includes Rt, Ep2, Grt and Ttn. Thin idioblastic corona of Gln2 occurs around Brs porphyroblast. (m) Composite corona

ED

of Chl + Ph occurs around Grt porphyroblast. (n-q) The developments of late Lws-bearing assemblages with or without Ab, (F)Wnc, Agt and Spn as veins in Grt (n-p) and Brs (q)

PT

porphyroblasts. (r-v) Sample N22F showing peak minerals of Omp porphyroblast and small

CE

Rt and Grt (r-s). Grt crystals, with an idioblastic habit and often in clusters occur in a matrix of retrograde Brs, Gln and Spn. Prograde inclusion of Gln1 occurs within the Omp

AC

porphyroblast. An intergrowth of Brs + Spn occurs around Omp and includes relics of Grt, Rt and Ep (s-u). Brs is later rimmed by Chl and Gln2 (t-v). Fig. 5: Mn-, Mg-, Fe-, and Ca- X-ray element images of representative porphyroblastic Grt crystals in Thewati eclogite samples N23B (a-h) and N40B (Grt-A in Fig. i-l and Grt-B in Fig. m-p) and corresponding compositional profiles along selected traverses (q-aj). The compositional features are attributed to growth zoning in Grt.

ACCEPTED MANUSCRIPT Fig. 6: (a-c) Ternary quadrilateral pyroxene-jadeite-aegirine diagram showing the compositional features of three textural varieties of Na-Ca Px crystals in the investigated samples: (i) Type 1 (T1) pyroxene as inclusions within Grt (a); (ii) Type 2 pyroxenes in the matrix (T2A type) or as relics in barroisite (T2B type) (b); (iii) T3 pyroxenes in veins (c). (df) Compositional varieties of Na-Ca (d-e) and Na amphiboles (f). (g) Binary NaB vs.

SC RI PT

(Al+Fe3+)VI diagram showing coupled substitution vectors in Na-Ca (dark-grey field) and Na(grey field) amphiboles.

Fig. 7: Colour coded mineral maps of two selected compositional layers in the banded Thewati eclogite sample N40B: Pyroxene-garnet-rich layer (a) and Barroisite-epidote-rich

NU

layer (b) using the software XMapTools (after Lanari et al., 2014). These two compositional layers are chosen for the computations of effective bulk rock compositions for phase

MA

equilibria modelling.

Fig. 8:(a) NCKFMASHTO P-T pseudosection for the Px-Grt-rich bulk rock compositional

ED

domain in sample N40B, showing stability of equilibrium mineral assemblage fields and key index minerals. The P-T estimate (with 2σ error) plotted in the diagram refers to the peak Ep-

PT

absent, eclogite facies metamorphism for the Thewati eclogite. (b-e) The mineral assemblage

CE

fields are contoured with isopleths of Grt compositions [x(g) and z(g) parameters, see Table 4 for the nomenclature] and modes of Grt [abbreviated as (Grt(m)], Omp [O(m)] and Amp

AC

[Amp(m)]. The reconstructed metamorphic reaction pathways in the Thewati eclogites are evaluated with compositional (e.g. Grt composition) and modal variations of these minerals along three representative metamorphic P-T paths of evolution (P1, P2 and P3). (f) A summary P-T diagram showing the stability fields of key index minerals, the P-T condition of formation

of

the

prograde

Grt

core

and

the

possible

sequences

of

mineral

appearances/disappearances along the three model P-T paths (see text for details). The numbers in Fig. 9a refer to different mineral assemblage fields, listed as: (1) Cpx Ph Ep Grt

ACCEPTED MANUSCRIPT Lws Amp Qz Rt, (2) Cpx Ph Ep Grt Lws Amp Rt, (3) Cpx Chl Ph Ep Grt Lws Amp Rt, (4) Cpx Ph Ep Grt Amp Rt, (5) Cpx Chl Ph Ep Amp Rt, (6) Cpx Chl Ph Amp Lws Spn Rt, (7) Cpx Ph Ep Grt Lws Qz Rt, (8) Cpx Ph Grt Amp Lws Qz Rt, (9) Cpx Ph Grt Amp Lws Rt, (10) Cpx Chl Ph Amp Lws Spn, (11) Cpx Ph Grt Lws Rt, (12) Cpx Pl Ph Grt Qz Rt, (13) Cpx Chl Ph Ep Lws Amp Rt, (14) Cpx Ph Grt Lws Coe Rt, (15) Cpx Ph Grt Amp Qz Rt, (16) Cpx

SC RI PT

Chl Ph Grt Lws Rt.

Fig. 9: (a) NCKFMASHTO P-T pseudosection for the Amp-Ep-rich bulk rock compositional domain in sample N40B, showing the stability of equilibrium mineral assemblage fields and key index minerals. (b-e) The mineral assemblage fields are contoured with compositional

NU

isopleths of Grt and Amp (NaB values at 0.5 and 1.5, where the values represent cations per formula unit) and modes of Grt, Omp and Amp. The prograde P-T path P1 predicts a Grt

MA

compositional zoning with progressive depletion in Grs + Adr content from core to inner rim of Grt followed by its enrichment toward the Grt outer rim, consistent with the measured Grt

ED

composition. The Amp compositional isopleths reveal the P-T windows at which Na-Ca (cf. Brs) and Na- (cf. Gln) amphiboles are likely to become stable along the retrograde path, P3.

PT

(f) The summary P-T diagram showing the stability fields of key index minerals and the

CE

possible sequences of mineral appearances along the two model retrograde P-T paths (see text for details). The numbers in Fig. 9a refer to different mineral assemblage fields, listed as:

AC

(1) Cpx Chl Ph Grt Lws Ep Spn, (2) Cpx Chl Ph Amp Lws Ep Spn, (3) Cpx Chl Ph Grt Amp Lws Ep Rt, (4) Cpx Chl Ph Grt Lws Ep Rt, (5) Cpx Ph Amp Pl Ep Rt, (6) Cpx Ph Amp Grt Pl Ep Rt, (7) Cpx Ph Grt Lws Qz Rt, (8) Cpx Ph Grt Lws Coe Rt, (9) Cpx Ph Grt Ep Lws Qz Rt, (10) Cpx Ph Grt Ep Lws Rt, (11) Cpx Ph Grt Amp Ep Rt, (12) Cpx Ph Amp Ep Rt, (13) Cpx Ph Amp Ep Qz Rt, (14) Cpx Ph Amp Pl Ep GrtQz Rt, (15) Cpx Chl Ph Amp Ep Lws Rt, (16) Cpx Chl Ph Ep Lws Rt, (17) Cpx Chl Ph Lws Ep Spn, (18) Cpx Ph Grt Ep Coe Rt, (19) Cpx Ph Grt Ep Lws Coe Rt.

ACCEPTED MANUSCRIPT Fig. 10: Reconstruction of the metamorphic P-T path of evolution of the Thewati eclogite using the following approaches: (a) plotting the P-T locus of the first appearance (shown by the symbol +) or disappearance (symbol -) of key minerals calculated through phase equilibria modelling, (b) calculating the P-T positions of activity corrected mineral reactions R3 to R6 (see text for details for the calculation procedures), (c) placing the mineral

SC RI PT

equilibria, R7 (after Maruyama et al., 1986), R8 (after Bohlen and Boettcher, 1982), R9 (after Holland, 1980) and experimentally determined Lws-out curve after Liou (1971) and Schmidt (1995) (Curve R10) in the P-T space and (d) the estimated P-T conditions of formation of prograde Grt core following the isopleth thermobarometry method and peak metamorphism

NU

through conventional thermobarometry. The deduced metamorphic evolutionary pathways when evaluated with (a) to (d) indicate a CW P-T path of evolution involving a prograde

MA

segment (MPR stage) in the epidote blueschist facies, and culminating in peak metamorphism in the eclogite facies (cf. MP stage) and three successive stages of retrograde evolution in

ED

amphibole eclogite facies (cf. MR1 stage), epidote blueschist facies (cf. MR2 stage) and transitional between greenschist and lawsonite blueschist facies (cf. M R3 stage). Mineral

PT

reactions: R3, 10(Jd)Omp+ (4Gr+5Prp)Grt+ 11Qtz + 8H2O → 5(Gl)Brs+6(Cz)Ep; R4, 5(Gl)Brs+

CE

12(Cz)Ep+ 14Qtz + 28H2O → 3(Tr)Fe-Wnc+ 10 Ab + 18 Lws; R5, 3(Dph)Chl + 5(Prp)Grt= 3(Clc)Chl + 5(Alm)Grt; R6, Ab = Jd17 + Qtz; R7, 25Gln + 6Cz + 7 Qtz + 14H2O = 6Tr + 9Chl

AC

+ 50Ab; R8; Qtz = Coe; R9; LAb = Jd + Qtz. Fig. 11: (a-b) BSE images showing the textural locations of zircon in N22F. (c-d) BSE images showing the distribution patterns of mineral inclusions and shape of host zircon crystal in N22F (c) and N23B (d). (e-f) Cathodoluminescence (CL) (e) and BSE (f) images of zircon showing the presence of Grt and Omp inclusions. (g-x) CL structures of representative analysed zircon grains in samples N23B (g-n) and N22F (o-x). 1-2-3 represent different CL domains (see text for details).

ACCEPTED MANUSCRIPT Fig. 12: Binary plots of Th vs. U compositions (in ppm) (a) and Th/U ratio vs. age (in Ma) (b) of analysed zircon grains in the two eclogite samples. Fig. 13: Concordia plots (a, c-d, f-g) and probability distribution diagrams (b, e) of 206Pb/238U ages (in Ma) in zircon grains from samples N22F (a-d) and N23B (e-g). Compositional ellipses in the concordia diagrams are linked with individual zircon grains (abbreviated as Z-

SC RI PT

10 etc.).

Fig. 14: Peak P-T conditions of metamorphism and metamorphic P–T path of the Nagaland eclogites (e.g. the Thewati eclogite path 1, this study) are compared with those from the Mokie blueschist (path 2) and Satuza eclogite (Path 3), all from the NOC, Himalayan

NU

blueschists (paths 4-6) and jadeitite and associated HP metamorphic rocks in the Central Ophiolite Belt (COB) of Myanmar. The P–T results/paths are plotted in the metamorphic

MA

facies diagram (after Okamoto and Maruyama, 1999). The apparent peak thermal gradients of the Nagaland eclogites are compared with model geotherms. The P-T field of global LT

ED

eclogites, and outlined by a red circle is adopted from Fig. 3a of Agard et al. (2016). See text for details. Abbreviations used: AM, Amphibolites facies; Ep-EC, Epidote eclogite facies;

PT

Amp-EC, Amphibole eclogite facies; BS, Blueschist facies; Dry-EC, Dry eclogite facies; EA,

CE

Epidote amphibolites facies; GS, Greenschist facies. References for metamorphic P-T paths and P-T conditions of metamorphism: 1-6: Metamorphic P-T paths of Neo-Tethyan HP/LT

AC

rocks; 1, Thewati EC (N40B), NOC (this study); 2, Mokie BS, NOC (after Ao & Bhowmik, 2014); 3, Satuza EC, NOC (Chatterjee & Ghose, 2010); 4,Sapi-Shergol-Zildat BS (Mahéo et al., 2006); 5, Sapi-Shergol LBS (Groppo et al., 2016); 6, Shangla BS (Guiraud, 1982, Jan, 1985); 7-11: P-T conditions of metamorphism of jadeitite and associated HP metamorphic rocks in the COB; 7, Mével et al., 1986; 8, Goffé et al., 2000; 9, Shi et al., 2003; 10, Oberhänsli et al., 2007; 11a/11b, Gln schist/Grt-Mica schist (after Thet Tin Nyunt, 2009). (b) P-T diagram comparing the peak P-T conditions of pumpellyite-diopside (PD), BS and

ACCEPTED MANUSCRIPT eclogite facies metamorphic rocks from the Mokie-Satuza (red circles) and Thewati (green circles) areas in the NOC. For a comparison, the P-T gradient of the subducted oceanic crust in a thermo-mechanical model-simulated intra-oceanic subduction system (after Gerya et al., 2002) is also plotted. The Nagaland HP/LT metamorphic rocks lie on the P-T gradient predicted by the thermo-mechanical model.

SC RI PT

Fig. 15: P–T path of the Nagaland eclogite (Path 1, NOC, this study) is compared with that from oceanic eclogites in metamorphic belts of the world. The grey P-T band refers to the average low T/P gradients for mature subduction (after Agard et al., 2018). The curve G02 (after Gerya et al., 2002) is the numerical simulation of prograde paths computed for

NU

subduction interface rocks. References for other paths: 2, Zermatt-Saas, Western Alps (Reinecke, 1998; Angiboust et al., 2009); 3, Dominican Republic (Krebs et al., 2011); 4,

MA

Monviso, Alps (Groppo and Castelli, 2010, Angiboust et al., 2012); 5: Rio San Juan mélange Complex (Krebs et al., 2011); 6: Sulawesi, Indonesia (Setiawan et al., 2016).

ED

Fig. 16: A simplified geological map of the NOC showing a progressive metamorphic sequence from GS through PD, BS to EC facies ensemble of basalt-limestone and chert. Also

PT

shown are the locations of Middle-Upper Jurassic-aged basalt and radiolarian chert samples

CE

of previous studies and PD, BS and eclogite facies samples with peak P-T estimates (Ao and Bhowmik, 2014 and this study). Ages of protolith formation and peak eclogite facies

AC

metamorphism (abbreviated as P and M respectively) of the eclogite samples (this study) are also indicated. To the east, the ophiolite terrane is bounded by Early Paleozoic Naga Metamorphics, a Lhasa/Tengchong-equivalent continental slice. Fig. 17: Cartoon diagrams showing the different tectonic models that correlate the Western, Central and Eastern Ophiolite belts in NE India and adjoining Myanmar and SE China with the Tibetan sutures (modified after Zhang, J.et al., 2018). We suggest that the HP metamorphic rocks of the Nagaland Ophiolite Complex and Jade Mines Belt together

ACCEPTED MANUSCRIPT constitute once continuous Early-Middle Jurassic accretionary complex within the NeoTethys.

Table 1: Representative mineral chemical analyses of garnet in eclogites from the NOC. Table 2: Representative mineral chemical analyses of pyroxene in eclogites from the NOC.

sphene, lawsonite, albite and rutile from the NOC.

SC RI PT

Table 3: Representative mineral chemical analysis of phengite, epidote, amphibole, chlorite,

Table 4: Results of thermobarometry of peak eclogite facies metamorphism in the NOC. Table 5: Calculations of effective bulk rock compositions of two representative chemical

NU

domains in sample N40B.

Table 6: SHRIMP U-Pb analyses of representative zircon grains in Mokie and Thewati

MA

eclogites from the NOC.

AC

CE

PT

ED

Table 7: A summary of characteristics of zircon from the Nagaland eclogite samples.

ACCEPTED MANUSCRIPT Table 1. Representative mineral chemical analysis of garnet in eclogite from the Nagaland Ophiolite complex, India N23B-1 133_I1 Grt C 37.14 0.17 21.03 0.01 27.08 2.90 2.21 9.56 0.09 b.d.l. 100.18 12 2.95 1.95 0.01 — 0.04 1.74 0.19 0.26 0.81 0.01 — 7.96 0.06 0.09 0.58 0.25 0.02 0.13

N23B-1 40_23_I1 Grt IR 37.13 0.06 21.13 0.13 30.11 0.49 3.34 7.99 0.04 0.02 100.42 12 2.93 1.94 — 0.01 0.05 1.91 0.03 0.39 0.67 0.01 — 7.93 0.01 0.13 0.64 0.20 0.02 0.17

N23B-1 57_1F2 Grt OR 37.78 0.12 21.13 b.d.l. 28.39 0.19 3.81 8.43 0.02 0.01 99.88 12 2.98 1.96 0.01 — 0.04 1.83 0.01 0.45 0.71 — — 7.98 b.d.l. 0.15 0.61 0.22 0.02 0.20

N23B-1 10_1F1 Grt OR 38.57 0.05 21.52 0.08 28.24 0.37 4.17 8.96 0.03 b.d.l. 101.98 12 2.97 1.95 — — 0.05 1.76 0.02 0.48 0.74 — — 7.96 0.01 0.16 0.59 0.22 0.02 0.21

N23B-3 1_1_I2 Grt C 37.76 0.62 20.71 0.04 25.60 4.40 2.09 9.86 0.06 b.d.l. 101.14 12 2.97 1.92 0.04 — 0.05 1.63 0.29 0.24 0.83 0.01 — 7.98 0.10 0.08 0.54 0.25 0.02 0.13

N23B-3 1_34_I2 Grt IR 37.92 0.02 21.35 b.d.l. 28.75 0.46 4.06 8.13 b.d.l. 0.02 100.71 12 2.97 1.95 — — 0.04 1.82 0.03 0.47 0.68 — — 7.96 0.01 0.16 0.61 0.20 0.02 0.20

SC RI PT

N23B-1 8_1F2 Grt C 37.32 0.11 20.82 0.05 25.64 3.55 1.70 10.67 0.01 b.d.l. 99.84 12 2.97 1.95 0.01 — 0.05 1.65 0.24 0.20 0.91 — — 7.97 0.08 0.07 0.55 0.28 0.02 0.11

NU

N23B-1 40_2_I1 Grt C 37.85 0.16 21.35 b.d.l. 27.87 3.14 2.17 9.29 0.03 b.d.l. 101.85 12 2.96 1.95 0.01 — 0.04 1.77 0.21 0.25 0.77 — — 7.96 0.07 0.08 0.59 0.24 0.02 0.12

MA

Sample Anal no. Mineral Site SiO₂ TiO₂ Al₂O₃ Cr₂O₃ FeO MnO MgO CaO Na₂O K₂O Total Oxygen Si Al Ti Cr Feᵌ⁺ Fe²⁺ Mn Mg Ca Na K Sum XSps XPrp XAlm XGrs XAdr XMg

N23B-3 1_39_I2 Grt OR 38.15 0.06 21.44 b.d.l. 27.73 0.40 3.97 9.57 0.08 b.d.l. 101.41 12 2.95 1.94 — — 0.05 1.73 0.03 0.45 0.79 0.01 — 7.96 0.01 0.15 0.58 0.24 0.03 0.21

N40B-4 5 _I4 GrtA C 37.85 0.25 20.81 b.d.l. 26.02 3.62 2.56 9.27 0.04 b.d.l. 100.42 12 2.99 1.94 0.01 — 0.05 1.67 0.24 0.30 0.78 0.01 — 8.00 0.08 0.10 0.56 0.24 0.02 0.15

N40B-4 30_I4 GrtA C 38.06 0.25 20.81 0.02 25.66 4.06 2.48 9.47 0.05 b.d.l. 100.86 12 3.00 1.93 0.02 — 0.05 1.64 0.27 0.29 0.80 0.01 — 8.00 0.09 0.10 0.55 0.24 0.03 0.15

N40B 35_I4 GrtA OR 38.03 0.11 21.17 b.d.l 27.49 0.36 4.20 9.17 0.04 b.d.l 100.5 12 2.97 1.94 0.01 — 0.06 1.73 0.02 0.49 0.76 0.01 — 7.97 0.01 0.16 0.58 0.23 0.03 0.21

AC

CE

PT

ED

Calculation of Fe3+ in Grt after Ravna & Terry (2004) and in pyroxenes, phengite, amphiboles, epidote and chlorite are bas abbreviations are after Whitney and Evans (2010); other abbreviations used: C/IR/OR/R: Core/Inner Rim/Outer Rim/Rim; I/M: In Porphyroblast/Vein; Fol/Cor: Foliation domain/Corona; Garnet and Pyroxene end members are in mole fraction; X Mg= Mg/( b.d.l: below detection limit.

Table 1. Continued Sample Anal no. Mineral Site SiO₂ TiO₂ Al₂O₃ Cr₂O₃ FeO MnO MgO CaO Na₂O K₂O Total Oxygen

N40B-4 P114_I3 GrtB Edge 38.81 0.04 21.81 b.d.l. 27.11 0.32 3.87 9.16 b.d.l. b.d.l. 101.11 12

N40B-4 23_I3 GrtB OR 38.30 0.04 20.95 0.08 28.52 0.31 4.32 8.41 0.05 b.d.l. 100.99 12

N40B-4 66_I3 GrtB OR 37.89 0.07 20.82 b.d.l. 27.47 0.34 4.20 8.72 b.d.l. b.d.l. 99.53 12

N40B-4 P96_I3 GrtB OR 38.43 0.01 21.23 b.d.l. 27.00 0.43 4.30 8.64 0.03 b.d.l. 100.06 12

N40B-4 15_I3 GrtB OR 38.27 0.04 21.13 b.d.l. 26.42 0.38 4.16 9.02 b.d.l. 0.02 99.44 12

N40B-6 P3_I3 GrtB C 37.08 0.20 20.19 b.d.l. 26.00 4.72 2.10 9.09 0.09 b.d.l. 99.46 12

N40B-6 P56_I3 GrtB IR 37.19 0.17 20.00 b.d.l. 28.28 0.43 3.88 8.34 0.02 b.d.l. 98.31 12

N40B-6 58_I3 GrtB IR 37.73 0.04 20.82 b.d.l. 26.48 0.49 4.25 9.16 b.d.l. b.d.l. 98.96 12

N40B-6 P61_I3 GrtB OR 37.69 0.07 20.38 b.d.l. 26.35 0.43 4.22 9.50 0.01 b.d.l. 98.63 12

N22F 69_F2 Grt C 37.47 0.22 20.16 b.d.l. 21.05 11.27 1.81 7.77 0.10 0.01 99.86 12

N22F 53_F2 Grt OR 37.63 0.06 20.88 0.02 23.29 6.23 2.54 8.11 0.10 b.d.l. 98.85 12

N22F 3_I1 Grt I(R)^Ttn 37.25 0.35 20.77 b.d.l. 26.21 7.78 3.02 5.34 0.06 b.d.l. 100.78 12

ACCEPTED MANUSCRIPT 2.98 1.92 — — 0.07 1.78 0.02 0.50 0.70 0.01 — 7.98 0.01 0.17 0.59 0.20 0.04 0.22

2.99 1.93 — — 0.06 1.75 0.02 0.49 0.74 — — 7.99 0.01 0.16 0.58 0.21 0.03 0.22

3.01 1.97 — — 0.03 1.74 0.03 0.50 0.73 — — 8.00 0.01 0.17 0.58 0.23 0.02 0.22

3.02 1.97 — — 0.03 1.72 0.03 0.49 0.76 — — 8.00 0.01 0.16 0.57 0.24 0.01 0.22

2.97 1.90 0.01 — 0.09 1.65 0.32 0.25 0.78 0.01 — 7.99 0.11 0.08 0.55 0.22 0.04 0.13

2.99 1.89 0.01 — 0.10 1.79 0.03 0.46 0.72 — — 7.99 0.01 0.15 0.60 0.19 0.05 0.21

2.99 1.94 — — 0.06 1.69 0.03 0.50 0.78 — — 7.99 0.01 0.17 0.56 0.23 0.03 0.23

SC RI PT

3.01 2.00 — — — 1.76 0.02 0.45 0.76 — — 8.00 0.01 0.15 0.59 0.25 0.00 0.20

NU

Si Al Ti Cr Feᵌ⁺ Fe²⁺ Mn Mg Ca Na K Sum XSps XPrp XAlm XGrs XAdr XMg

3.00 1.91 — — 0.09 1.66 0.03 0.50 0.81 — — 7.99 0.01 0.17 0.55 0.23 0.04 0.23

3.00 1.91 0.01 — 0.07 1.35 0.77 0.22 0.67 0.02 — 8.00 0.26 0.07 0.19 0.45 0.04 0.14

3.02 1.99 — — 0.01 1.57 0.43 0.31 0.70 0.02 — 8.00 0.14 0.10 0.23 0.52 0.00 0.16

2.96 1.93 0.02 — 0.05 1.68 0.52 0.35 0.45 0.01 — 7.96 0.17 0.12 0.13 0.56 0.03 0.17

Table 2. Representative mineral chemical analysis of pyroxenes in eclogite from the Nagaland Ophiolite Complex, India. N23B-1 9_F1 Omp I^Grt (R) 55.80 0.11 8.01 0.01 8.57 b.d.l. 7.93 13.68 6.45 0.01 100.57 6 2.00 0.34 — — 0.10 0.16 — 0.42 0.53 0.45 — 4.00 0.55 0.35 0.10 0.73

MA

N23B-1 10_1F2 Agt I^Grt (C) 53.62 0.04 4.50 0.03 11.01 0.17 9.10 16.20 4.82 0.01 99.49 6 1.98 0.20 — — 0.20 0.14 0.01 0.50 0.64 0.34 — 4.00 0.65 0.16 0.19 0.78

ED

N23B-1 7_1F2 Omp I^Grt (C) 53.67 0.08 4.26 0.01 10.17 0.26 9.67 16.71 4.17 b.d.l. 99.00 6 1.99 0.19 — — 0.13 0.19 0.01 0.53 0.66 0.30 — 4.00 0.70 0.17 0.13 0.74

PT

N23B-1 15_F1 Omp I^Grt (C) 54.89 0.04 5.32 b.d.l. 10.85 0.13 8.73 16.15 5.05 0.02 101.16 6 1.99 0.23 — — 0.15 0.18 — 0.47 0.63 0.35 — 4.00 0.64 0.21 0.15 0.73

CE

Site SiO₂ TiO₂ Al₂O₃ Cr₂O₃ FeO MnO MgO CaO Na₂O K₂O Total Oxygen Si Al Ti Cr Feᵌ⁺ Fe²⁺ Mn Mg Ca Na K Sum XQuad XJd XAeg XMg

N23B-1 14_F1 Omp I^Grt (C) 52.57 0.01 5.53 0.01 12.95 0.19 9.30 14.47 4.38 0.08 99.48 6 1.95 0.24 — — 0.18 0.22 0.01 0.51 0.57 0.31 — 4.00 0.67 0.16 0.16 0.70

AC

Sample Anal no. Mineral

N23B-1 53_1F2 Omp I^Grt (R) 55.59 0.16 8.58 0.03 8.42 0.03 7.26 12.76 6.83 b.d.l. 99.66 6 2.01 0.37 — — 0.09 0.17 — 0.39 0.49 0.48 — 3.99 0.52 0.39 0.09 0.70

N23B-1 13_4F1 Omp

N23B-1 16_4F1 Omp

P(C)

P(R)

55.81 0.08 6.68 0.04 10.47 0.31 7.39 11.82 7.39 0.01 99.99 6 2.02 0.28 — — 0.20 0.12 0.01 0.40 0.46 0.52 — 4.00 0.48 0.30 0.21 0.77

55.51 0.15 5.78 b.d.l. 11.31 0.23 8.03 12.06 6.94 0.01 100.03 6 2.01 0.25 — — 0.21 0.13 0.01 0.43 0.47 0.49 — 4.00 0.51 0.26 0.23 0.99

N23B-1 59 _IF2 Omp Relics^ Brs(P) 54.45 b.d.l. 7.97 0.04 7.75 b.d.l. 8.45 13.93 6.16 0.02 98.76 6 1.98 0.34 — — 0.13 0.11 — 0.46 0.54 0.44 — 4.00 0.56 0.32 0.12 0.81

N40B-4 4_I4 Omp I^GrtA (C) 54.90 0.06 6.04 0.06 9.61 0.18 8.66 15.09 5.62 b.d.l. 100.21 6 1.99 0.26 — — 0.15 0.14 0.01 0.47 0.59 0.40 — 4.00 0.60 0.25 0.15 0.77

N40B-4 29_I4 Omp I^GrtA (C) 54.69 0.05 6.46 0.01 10.10 0.15 8.45 14.97 5.85 0.01 100.74 6 1.97 0.28 — — 0.19 0.12 0.01 0.45 0.58 0.41 — 4.00 0.58 0.24 0.18 0.80

N40B40_I4 Omp I^GrtA (R) 54.19 0.08 8.37 0.05 9.57 0.13 7.36 12.99 6.77 0.01 99.50 6 1.97 0.36 — — 0.18 0.11 — 0.40 0.51 0.48 — 3.99 0.52 0.31 0.17 0.79

ACCEPTED MANUSCRIPT

Table 2. Continued N40B-4 42_I3 Omp

Fol(C)

Fol(R)

54.94 0.01 8.16 b.d.l. 6.81 b.d.l. 8.64 14.45 5.94 0.03 98.98 6 2.00 0.35 — — 0.08 0.13 — 0.47 0.56 0.42 — 4.00 0.58 0.34 0.08 0.79

55.23 0.10 7.93 0.06 7.17 0.08 8.43 14.37 6.01 b.d.l. 99.37 6 2.00 0.34 — — 0.07 0.14 — 0.46 0.56 0.42 — 3.99 0.58 0.35 0.08 0.76

N40B-4 84_I3 Agt V(+Lws)^ GrtB 53.23 0.06 3.79 b.d.l. 19.92 0.33 3.81 6.77 9.47 b.d.l. 97.38 6 2.01 0.17 — — 0.50 0.13 0.01 0.21 0.27 0.69 — 4.00 0.31 0.17 0.52 0.63

N40B-4 23_I4 Agt V(+Ab)^ GrtA 53.74 0.22 4.27 0.02 21.82 0.33 2.23 4.84 10.97 0.01 98.44 6 2.00 0.19 0.01 — 0.59 0.10 0.01 0.12 0.19 0.79 — 4.00 0.21 0.19 0.60 0.57

N40B-4 74_I3 Agt V^MOmp 53.92 0.03 5.32 0.01 16.25 0.25 4.63 7.91 9.21 0.03 97.56 6 2.01 0.23 — — 0.41 0.10 0.01 0.26 0.32 0.67 — 4.00 0.33 0.24 0.42 0.73

SC RI PT

N40B-4 41_I3 Omp

NU

N40B-4 14_I3 Omp I^GrtB (OR) 54.67 b.d.l. 7.77 0.07 8.70 0.06 7.73 13.77 6.25 b.d.l. 99.01 6 2.00 0.33 — — 0.12 0.15 — 0.42 0.54 0.44 — 4.00 0.56 0.33 0.12 0.74

MA

Site

N40B-6 11_I3 Omp I^GrtB (OC) 54.78 0.10 7.51 b.d.l. 8.57 0.20 8.08 14.38 6.08 0.02 99.72 6 1.99 0.32 — — 0.13 0.13 0.01 0.44 0.56 0.43 — 4.00 0.57 0.30 0.13 0.77

N40B-6 32_I3 Omp I^GrtB (IR) 54.34 0.08 8.40 b.d.l. 8.99 0.19 7.12 13.18 6.71 0.01 99.02 6 1.98 0.36 — — 0.15 0.13 0.01 0.39 0.52 0.47 — 4.00 0.52 0.33 0.14 0.75

N40B-6 39_I3 Omp I^GrtB (IR) 54.82 0.05 9.25 b.d.l. 8.58 0.10 7.27 12.56 6.92 b.d.l. 99.54 6 1.98 0.39 — — 0.13 0.13 — 0.39 0.49 0.48 — 4.00 0.51 0.37 0.12 0.75

N40B-6 48_I3 Omp Fol 54.55 0.01 7.63 0.01 9.31 0.02 7.75 13.37 6.55 0.03 99.23 6 1.98 0.33 — — 0.17 0.12 — 0.42 0.52 0.46 — 4.00 0.53 0.30 0.16 0.78

CE

PT

SiO₂ TiO₂ Al₂O₃ Cr₂O₃ FeO MnO MgO CaO Na₂O K₂O Total Oxygen Si Al Ti Cr Feᵌ⁺ Fe²⁺ Mn Mg Ca Na K Sum XQuad XJd XAeg XMg

N40B-4 65_I3 Omp I^GrtB (OR) 55.23 0.06 9.21 0.02 7.90 0.07 7.46 12.58 6.93 b.d.l. 99.46 6 1.99 0.39 — — 0.10 0.14 — 0.40 0.49 0.49 — 4.00 0.51 0.38 0.10 0.75

ED

Sample Anal no. Mineral

AC

Table 3. Representative mineral chemical analysis of phengite, epidote, amphiboles, chlorite, sphene, lawsonite, albite and Nagaland Ophiolite Complex, India. Sample Anal no. Mineral

N23B-1 103_I1 Ph

N23B-1 56_2_I2 Ph

N40B-4 26 _I4 Ph

N40B-4 67_I3 Ph

N22F 44_I1 Ph

N22F 1_I4 Ep

N23B-1 2_4F1 Ep1

N40B-4 3_I4 Ep1

N40B-4 29_I3 Ep2

N40B-6 49_I3 Ep2

N23B-1 42_2_I2 Brs

N40B-4 93_I3 Brs

Site

I^Grt(R)

M

I^GrtA(C)

M(C)

M

I^Omp

I^Grt(C)

I^Grt(C)

Fol

I^Brs(P)

Fol

P

SiO₂ TiO₂ Al₂O₃ Cr₂O₃ FeO MnO MgO CaO Na₂O K₂O Total Oxygen

53.57 0.10 23.25 0.18 5.06 0.05 4.13 0.04 0.17 10.93 97.48 11

52.12 0.25 23.25 0.02 3.74 b.d.l. 4.22 0.01 0.14 11.12 94.87 11

51.85 0.07 22.11 b.d.l. 5.32 0.15 3.49 0.08 0.12 11.09 94.28 11

52.66 0.43 26.03 b.d.l. 3.87 b.d.l. 3.18 b.d.l. 0.67 10.49 97.33 11

55.59 0.02 21.86 b.d.l. 2.55 0.09 5.87 b.d.l. 0.12 11.60 97.70 11

38.13 0.17 26.62 0.13 8.35 0.21 0.11 23.23 0.01 0.00 96.95 12.50

37.80 0.07 23.86 0.01 11.51 0.24 0.04 22.78 b.d.l. b.d.l. 96.31 12.5

38.27 0.04 24.03 0.02 11.50 0.31 0.03 23.37 0.05 b.d.l. 97.62 12.5

38.06 0.16 24.46 b.d.l. 9.51 0.11 0.04 23.77 0.02 b.d.l. 96.13 12.5

38.26 0.06 25.13 b.d.l. 8.74 0.18 0.01 24.16 b.d.l. b.d.l. 96.54 12.5

51.70 0.17 7.43 b.d.l. 10.72 0.16 14.99 9.42 2.85 0.29 97.70 23

50.97 0.12 6.80 b.d.l. 12.89 0.15 13.36 8.18 3.22 0.20 95.88 23

ACCEPTED MANUSCRIPT 3.54 1.81 0.01 0.01 0.05 0.23 — 0.41 — 0.02 0.92 6.99 0.64 –

3.53 1.86 0.01 — — 0.21 — 0.43 — 0.02 0.96 7.02 0.67 –

3.57 1.79 — — — 0.31 0.01 0.36 0.01 0.02 0.97 7.03 0.46 –

3.46 2.02 0.02 — — 0.21 — 0.31 — 0.09 0.88 6.99 0.59 –

3.63 1.68 — — — 0.14 0.01 0.57 — 0.02 0.97 7.02 0.80 –

2.99 2.46 0.01 0.01 0.52 0.03 0.01 0.01 1.95 0.00 0.00 8.00 – 0.17

3.01 2.24 b.d.l. b.d.l. 0.76 0.01 0.02 b.d.l. 1.94 b.d.l. b.d.l. 7.99 – 0.25

3.01 2.23 — — 0.75 0.01 0.02 — 1.97 0.01 — 8.00 – 0.25

3.05 2.31 0.01 — 0.44 0.20 0.01 0.01 2.04 — — 8.07 – 0.16

3.03 2.34 — — 0.57 0.01 0.01 — 2.05 — — 8.01 – 0.20

7.34 1.24 0.02 — 0.22 1.06 0.02 3.17 1.43 0.78 0.05 15.33 0.75 –

7.42 1.17 0.01 — 0.30 1.27 0.02 2.90 1.28 0.91 0.04 15.31 0.70 –

NU

SC RI PT

Si Al Ti Cr Feᵌ⁺ Fe²⁺ Mn Mg Ca Na K Sum XMg XPs

I^Grt(C)

COR
55.81 b.d.l. 7.21 0.13 15.88 0.21 9.33 1.90 6.30 0.01 96.79 23 7.95 1.21 — 0.01 0.39 1.51 0.03 1.98 0.29 1.74 — 15.11 0.57 –

57.52 0.01 6.22 b.d.l. 14.52 0.01 11.23 0.96 6.81 b.d.l. 97.27 23 8.00 1.02 — — 0.81 0.88 — 2.33 0.14 1.84 — 15.01 0.73 –

ED

N23B-1 50_I2 Gln2

PT

SiO₂ TiO₂ Al₂O₃ Cr₂O₃ FeO MnO MgO CaO Na₂O K₂O Total Oxygen Si Al Ti Cr Feᵌ⁺ Fe²⁺ Mn Mg Ca Na K Sum XMg XPs

N23B-1 12_I1 Gln1

CE

Site

N40B-4 101_I3 Fwn V(+Lws) ^Brs 52.88 0.03 1.41 b.d.l. 25.59 0.57 7.21 2.88 5.19 0.04 95.79 23 7.95 0.25 — — 1.17 2.05 0.07 1.62 0.46 1.51 0.01 15.10 0.44 –

AC

Sample Anal no. Mineral

N40B-2 72_I3 Gln2

N22F 101_I2 Gln1

N22F 63_I2 Gln2

N23B-1 38_I2 Chl

70_I1 Chl

N40B-4 6_I4 Spn1

N23B-1 85_I1 Spn2

COR
I^Agt

Cor
Cor
Cor
I^Grt(C)

Cor
56.63 0.05 7.67 0.06 12.79 0.21 11.17 1.85 6.08 0.02 96.53 23 7.97 1.27 0.01 0.01 0.31 1.20 0.03 2.34 0.28 1.66 — 15.06 0.66 –

57.58 0.02 9.79 0.02 11.09 0.24 11.52 1.09 7.13 0.02 98.49 23 7.87 1.58 — — 0.30 0.97 0.03 2.35 0.16 1.89 — 15.13 0.71 –

57.95 b.d.l. 8.80 b.d.l. 10.71 0.16 12.05 1.28 6.97 0.05 97.97 23 7.95 1.42 — — 0.28 0.95 0.02 2.46 0.19 1.85 0.01 15.13 0.72 –

28.37 0.02 18.48 b.d.l. 18.33 0.61 21.65 0.04 0.05 0.01 87.56 14 2.88 2.21 — — 0.03 1.53 0.05 3.28 0.01 0.01 — 10.00 0.68 –

29.06 0.05 18.97 0.00 14.75 0.33 24.86 0.05 0.00 0.00 88.09 14.00 2.87 2.21 — — 0.04 1.18 0.03 3.66 0.01 0.00 0.00 10.00 0.76 –

30.10 38.43 1.07 0.04 1.12 0.14 b.d.l. 28.35 0.03 0.02 99.30 5 0.99 0.04 0.95 — — 0.03 — — 1.00 — — 3.02 – –

30.70 38.25 1.43 b.d.l. 0.54 0.11 0.01 29.45 0.06 0.01 100.57 5 1.00 0.05 0.94 — — 0.01 — — 1.03 — — 3.03 – –

MA

Table 3. Continued

N22F

ACCEPTED MANUSCRIPT

ts of thermobarometry of peak eclogite facies metamorphism in the NOC. Grt

Cpx

Simultaneous solution of Grt-Cpx thermometry and Grt-Cpx-Phe barometry

Phe

Thermocalc average Cpx-Ph ass

t)

z(grt)

f(grt)

m(grt)

x(cpx)

j(cpx)

f(cpx)

Si(mu)

x(mu)

y(mu)

na(mu)

PRT , TRAV

PTC , TTC

PTC(Av.)(SdP),TTC(Av.)(S

0

0.24, ¹0.24

0.02

0.00

0.30

0.48

0.19

3.54

0.36

0.39

0.02

29.0±2, 650±60

28.7±3.7, 665±110

29.8±3.7, 685±70

29.1±3.7, 665±110

29.2±3.5, 650±65

0.24, ¹0.24

0.02

0.00

0.30

0.48

0.19

3.53

0.33

0.39

0.02

0

0.24, ¹0.24

0.02

0.00

0.27

0.45

0.23

3.54

0.36

0.39

0.02

27.3±2, 600±60

27.1±3.5, 615±90

29.3±3.6, 675±65

0

0.24, ¹0.24

0.02

0.00

0.27

0.45

0.23

3.53

0.33

0.39

0.02

27.6±2, 605±60

27.5±3.5, 620±90

28.7±3.5, 645±60

0.39

0.02

28.5±2, 640±60

28.1±3.4, 655±100

30.5±3.7, 710±70

0.33 0.39 Weighted average

0.02

28.8±2, 640±60

28.5±3.4, 660±100 28.1±2.8, 645±80

29.7±3.5, 675±65

28.4±1.6, 630±50 0.09

26.5±2, 675±60

26.0±3.3, 695±110

28.6±3.5, 770±75

0.09

25.6±2, 605±60

25.0±3.6, 625±100

30.2±4.7, 735±85

0.09

26.6±2, 645±60

26.0±3.1, 660±100

28.9±3.3, 765±75

9

0.25, ¹0.25

0.02

0.01

0.27

0.45

0.23

3.54

9

0.25, ¹0.25

0.02

0.01

0.27

0.45

0.23

3.53

0.36

8

0.24, ¹0.25

0.02

0.01

0.29

0.45

0.22

3.46

0.41

0.49

8

0.23, ¹0.24

0.04

0.01

0.26

0.51

0.25

3.46

0.41

0.49

8

0.26, ¹0.26

0.01

0.01

0.26

0.44

0.26

3.46

0.41

0.49

8

0.25, ¹0.25

0.03

0.01

0.25

0.49

0.21

3.46

0.41

0.49

Weighted average

7

0.26, ¹0.26

0.03

0.01

0.25

0.48

0.25

7

0.26, ¹0.26

0.03

0.01

0.25

0.47

0.30 Weighted average

0.15, ¹0.18

0.03

0.17

0.23

0.50

0.38

3.63

0.20

0.31

3

0.15, ¹0.18

0.03

0.17

0.24

0.46

0.41

3.63

0.20

0.31

4

0.16, ¹0.19

0.02

0.16

0.24

0.46

0.41

3.63

0.20

0.09

0.31

Weighted average

25.6±2, 605±60

24.8±3.4, 620±100

29.7±4.6, 739±90

25.5±3.3, 650±100

29.2±3.8, 755±80

²635±60

²645±95 ²640±90 ²640±130

0.02

²630±85 24.2±2, 550±60

24.1±4.8, 550±80

0.02

23.6±2, 565±60

23.4±4.4, 570±80

0.02

29.5±2.9, 670±50

26.1±2.0, 635±60 ²625±60

NU

3

SC RI PT

0

29.4±2, 650±60

24.0±2, 540±60 23.9±2.3, 550±70

24.0±4.3, 545±75 23.8±5.1, 555±90 2+

2+

#

24.1±4.2, #475±80

#

22.9±4.1, #470±80

#

23.1±3.6, #455±75 23.3±4.5, 465±90

AC

CE

PT

ED

MA

used: Chemical compositions of garnet (grt), clinopyroxene (cpx) and phengite (mu), x(mineral)=Fe /(Fe +Mg); z(grt)=Ca/(Ca+M m(grt)=Mn/(Ca+Mg+Fe2++Mn); j(cpx)=Na(total); f(cpx)=Fe3+/(Fe3++Al(VI)); y(mu)=XAlM2A; na (mu)=XNaA ; TRAV/PRT , temperatur ed using Grt-Cpx thermometry of Ravna (2000) & Grt-Cpx-Ph barometry of Ravna and Terry (2004); TTC/PTC: temperature (⁰C)/ g Grt-Cpx thermometric and Grt-Cpx-Ph barometric reactions calculated by using Thermocalc software (v.3.31); T TC(Av)/PTC(Av): Therm C)/pressure (kbar) calculated using Thermocalc software (v.3.31) (see text for details); #Average P-T computations with Grt activity mo , z(grt)=Ca/(Ca+Mg+Fe2+); 2,temperature calculated at a PRef=26 kbar.

ACCEPTED MANUSCRIPT

Table 5: Calculations of effective bulk rock compositions of two representative chemical domains in sample N40B NOC. XMapTools Data

Mole %

52.44±1.57 0.11±0.07 5.97±1.96 — 15.09±3.01 0.22±0.15 12.31±1.39 6.81±2.09 3.88±0.83 0.18±0.08 —

50.31±0.88 0.24±0.08 12.96±1.66 0.03±0.03 12.77±2.52 0.19±0.17 6.29±0.93 10.76±1.26 4.04±0.58 1.38±0.05 —

52.56 0.19 7.98 0.00 11.27 0.17 9.79 12.04 4.09 0.92 0.99

97.00±4.73

98.97±3.70

100.00

Calculated Bulk

Mole%

48.50±0.45 0.37±0.22 13.57±0.67 0.02±0.04 11.97±0.84 0.36±0.26 6.40±0.39 13.77±0.47 3.34±0.21 0.71±0.07 — 99.01±1.56

50.48 0.29 8.34 0.00 10.52 0.32 9.93 15.36 3.37 0.47 0.92 100.00

Grt

Ep

Ttn

Ph

Omp

Brs

17.04 4.00 0.34 15.40 57.11 6.11 100.00

SiO2 TiO2 Al2O3 Cr2O3 FeO MnO MgO CaO Na2O K2O O2

38.07±0.46 0.07±0.05 20.82±0.34 0.02±0.03 27.49±0.88 0.53±0.17 3.88±0.32 8.80±0.46 0.04±0.03 0.00±0.01 —

38.40±0.30 0.09±0.03 24.77±0.64 0.04±0.04 10.11±1.03 0.06±0.06 0.10±0.11 23.40±0.50 0.03±0.02 0.00±0.00 —

30.96±0.09 36.92±1.64 2.09±0.97 — 0.86±0.47 0.04±0.06 — 29.09±0.11 0.02±0.03 0.02±0.01 —

52.84±1.25 0.34±0.12 25.12±1.80 0.01±0.02 3.96±0.21 0.03±0.03 3.56±0.51 0.01±0.02 0.45±0.22 10.68±0.29 —

54.78±0.93 0.07±0.07 7.40±2.13 0.03±0.04 9.65±3.68 0.12±0.22 7.65±1.25 13.20±1.80 6.56±0.87 0.01±0.01 —

Total

99.73±1.20

97.00±1.35

100.00±1.97

96.99±2.30

99.46±4.96

N40B-6

SC RI PT

Mineral

N40B-6

Sample

Average mineral composition(EPMA)

Mineral

Grt

Ep

Ttn

13.48 20.94 0.62 7.81 37.39 19.76 100.00

SiO2 TiO2 Al2O3 Cr2O3 FeO MnO MgO CaO Na2O K2O O2 Total

37.80±0.45 0.24±0.39 20.29±0.38 0.02±0.03 27.01±1.17 1.97±1.42 3.19±0.66 9.38±0.52 0.05±0.06 0.05±0.19 — 100±2.15

38.60±0.32 0.23±0.64 24.70±0.65 0.04±0.08 10.34±0.92 0.10±0.08 0.08±0.08 24.38±0.28 0.03±0.05 0.01±0.01 — 98.50±1.37

29.85±0.15 39.00±0.33 1.26±0.10 0.04±0.05 0.72±0.29 0.01±0.01 — 29.05±0.16 0.05±0.03 0.02±0.01 — 100.00±0.51

PT CE

Ph

Omp

Brs

51.92±1.93 0.35±0.25 23.72±2.73 0.01±0.02 4.27±0.16 0.06±0.04 3.89±0.80 0.07±0.10 0.43±0.27 11.28±0.50 — 96.00±3.50

54.96±0.34 0.07±0.04 7.93±0.67 0.02±0.03 8.48±0.99 0.08±0.07 8.08±0.45 14.13±0.65 6.26±0.35 0.01±0.01 — 100.00±1.52

51.24±0.38 0.15±0.04 6.67±0.21 0.03±0.05 13.53±0.30 0.11±0.08 13.41±0.18 8.35±0.29 3.27±0.14 0.24±0.03 — 97.00±0.65

NU

Mode

AC

Grt Ep Spn Ph Omp Brs Total

Calculated Bulk

Mode

XMapTools Data Phases

Average mineral composition(EPMA)

MA

Grt Ep Spn Ph Omp Brs Total

N40B-4

Sample

ED

Phases

N40B-4

ACCEPTED MANUSCRIPT Table 6: SHRIMP U-Pb analyses of representative zircon grains in Mokie and Thewati eclogites from the Nagaland Ophiolite Complex ²⁰⁴Pb /²⁰⁶Pb

% comm 206 Pb

4corr ²⁰⁷Pb* /²⁰⁶Pb*

% err

4corr ²⁰⁷Pb* /²³⁵U

% err

4corr ²⁰⁶Pb* /²³⁸U

% err

204corr ²⁰⁶Pb* /²³⁸U Age(Ma)

31.23

4.6E-4

0.84

.0477

5.5

0.177

5.7

.0269

1.4

171

0.15

35.05

7.0E-4

1.28

.0701

5.8

0.519

6.0

.0536

1.5

337

66

0.03

60.67

5.0E-4

0.92

.0476

4.9

0.212

5.1

.0322

1.4

204

1267

30

0.02

40.33

2.1E-4

0.38

.0600

7.5

0.305

7.6

.0369

1.4

234

Zrn-13(O)

1110

36

0.03

26.37

9.0E-4

1.64

.0446

9.3

0.167

9.4

.0272

1.6

173

Zrn-22(N1)

1705

13

0.01

46.00

2.3E-4

0.42

.0491

2.9

0.212

3.3

.0313

1.6

199

Zrn-28(N1)

1380

271

0.20

52.09

3.1E-4

0.57

.0528

5.2

0.318

5.5

.0437

1.6

276

Zrn-33(N1)

1091

94

0.09

41.31

1.4E-3

2.58

.0566

13.3

0.335

13.4

.0429

1.8

271

Zrn-35(N1)

1511

88

0.06

44.94

8.2E-4

1.51

.0496

7.9

0.233

8.1

.0341

1.7

216

Zrn-36(N1)

2387

15

0.01

63.24

1.0E-3

1.90

.0495

7.9

0.207

8.1

.0302

1.6

192

Zrn-40(N1)

1197

10

0.01

33.80

6.8E-4

1.24

.0565

6.4

0.253

6.6

.0325

1.7

206

Zrn-45(N1)

1939

24

0.01

50.06

3.5E-4

0.65

.0480

3.4

0.198

3.8

.0299

1.7

190

Zrn-46(N1)

2478

8

0.003

63.27

1.1E-4

0.19

.0501

1.8

0.205

2.4

.0297

1.6

188

Zrn-47(N1)

1469

40

0.03

38.31

1.9E-3

3.57

.0591

11.6

0.238

11.8

.0293

1.8

186

Zrn-49(N1)

1561

29

0.02

57.64

1.3E-4

0.24

.0555

3.4

0.328

3.8

.0429

1.7

271

Zrn-50(N1)

2051

36

0.02

52.05

1.8E-4

0.34

.0497

2.8

0.202

3.2

.0294

1.6

187

Zrn-55(N1) Sample N23B

1765

163

0.09

61.25

3.5E-3

6.49

.0553

16.6

0.288

16.7

.0378

1.7

239

Zrn-10

106

0.88

0.01

3.03

4.8E-3

8.77

.0339

91.2

0.142

91.4

.0305

5.6

194

Zrn-16

97

0.85

0.01

2.48

1.2E-3

2.11

.0488

21.1

0.196

21.3

.0291

2.9

185

Zrn-21

45

0.16

0.004

1.19

7.7E-3

14.05

.0694

84.9

0.253

85.3

.0265

7.8

168

Zrn-22

200

0.57

0.003

6.22

1.2E-2

22.72









.0280

5.3

178

Zrn-23

92

0.24

0.003

2.98

1.4E-2

25.27









.0281

7.0

178

Zrn-34

136

0.39

0.003

3.78

6.7E-4

1.22

.0479

19.0

0.212

19.3

.0320

3.5

203

Zrn-46

61

0.002

1.97

6.5E-3

11.20

.0756

65.2

0.347

65.6

.0333

7.0

211

²⁰⁴Pb /²⁰⁶Pb

% comm 206

4corr ²⁰⁷Pb* /²⁰⁶Pb*

% err

4corr ²⁰⁷Pb* /²³⁵U

% err

4corr ²⁰⁶Pb* /²³⁸U

% err

204corr ²⁰⁶Pb* /²³⁸U Age(Ma)

U (ppm)

Th (ppm)

Th/U

²⁰⁶Pb* (ppm)

Zrn-8(O)

1340

27

0.02

Zrn-9(O)

751

113

Zrn-10(O)

2172

Zrn-11(O)

Spot No.

NU

MA

ED

PT

CE

AC

Table 6: contd.

0.14

SC RI PT

Sample N22F

U (ppm)

Th (ppm)

Th/U

²⁰⁶Pb* (ppm)

57

0.18

0.003

1.84

6.9E-3

12.72

.0235

79.5

0.106

79.6

.0328

4.1

208

Zrn-53

126

0.23

0.002

3.43

1.0E-3

1.87

.0408

29.2

0.175

29.4

.0311

3.7

198

Zrn-54

55

0.12

0.002

1.52

3.6E-3

6.60









.0298

5.2

189

Zrn-62

59

0.11

0.002

1.70

6.4E-3

11.80









.0297

5.5

189

Zrn-64

42

0.15

0.004

1.32

1.5E-2

27.27









.0266

10.8

169

Zrn-66

29

0.04

0.001

0.88

7.9E-3

14.41









.0299

9.7

190

Zrn-66.2

13

0.06

0.004

0.47

2.2E-2

41.12









.0239

27.7

152

Zrn-69

55

0.24

0.004

1.70

7.8E-3

14.28









.0307

6.8

195

Zrn-71

83

0.27

0.003

2.52

1.4E-3

2.66

.0828

19.4

0.391

19.8

.0343

4.0

217

Zrn-76

93

1.62

0.02

2.48

2.4E-3

4.37

.0411

62.6

0.168

62.7

.0297

4.8

188

Spot No.

Zrn-49

ACCEPTED MANUSCRIPT

AC

CE

PT

ED

MA

NU

SC RI PT

Notes: Pb* indicate the radiogenic Pb corrected using measured ²⁰⁴Pb & ²⁰⁷Pb . The “-” means undetected or imponderable.

ACCEPTED MANUSCRIPT Table 7: A summary of characteristics of zircon from two Nagaland eclogite samples Stage-1

Stage-2

Stage-3

Dark luminescent, faintly banded elongated domain in Zrn cores (cf. Type 1 Zrn); Euhedral external shape, often with pyramidal termination

Only in N22F; Irregular bright luminiscent domain that replaces stage 2 Zrn cores (cf. Type 3 Zrn)

Composition

Low Th (0.12-0.27 ppm) and U (55-126 ppm) concentrations in N23B; More uraniferous (15112387 ppm) and thoriferous (15-88 ppm) in N22F; Low Th/U ratios in both sample groups : <0.01 (N23B); 0.01-0.06 (N22F) ~205 Ma

Commonly structure-less dark luminescent core, locally with patchy, irregular zones of bright luminescence (cf. Type 2 Zrn); Change of regular shape of Type-1 Zrn; Bright luminescent narrow rim around Zrn cores; Contains inclusions of Grt, Omp and Rt Low Th (0.05-1.62 ppm) and U (29-200 ppm) concentrations in N23B; uraniferous (1110-2478 ppm) and thoriferous (8-66 ppm) in N-22F; Low Th/U ratios in both sample groups : ≤0.02 (N-23B); ≤0.03 (N22F) ~189-185 Ma

Older 206Pb/238U dates between 337 and 239 Ma; Discordant in concordia plot Localised fluid-mediated alteration of recrystallised Zrn core

Statistically Th-rich (94271) and U-poor (751-1765 ppm) than other textural domains in N22F; High Th/U ratios (0.09-0.20)

MA

Recrystallization/replacement of growth-zoned Zrn due to in situ dissolution–reprecipitation in sub-solidus metamorphic conditions Metamorphic/ Minimum age of protolith Peak eclogite facies Timing unknown; Older Magmatic formation metamorphism at ~189-185 discordant dates context Ma geologically meaningless Mineral abbreviations after Whitney & Evans (2010); Other abbreviations used: CL, Cathodoluminescence; Soln, Solution.

CE

PT

ED

Primary growth due to precipitation from fluid/melt

AC

Origin

NU

Age

SC RI PT

Zrn evolution stages CL Zoning and texture

ACCEPTED MANUSCRIPT Highlights  Nagaland eclogites record the oldest evidence for Neo-Tethyan subduction in the YTSZ.  Minimum age of Nagaland eclogitic protolith is Late Triassic (≥ c. 205 Ma).  Peak eclogite facies metamorphism occurred between 189 and 185 Ma.

SC RI PT

 Metamorphic peak P-T conditions range from 25.5-28.1 kbar, ~650 oC to 23.8 kbar, ~555oC.

 Eclogites record low thermal gradient (~7-8 oC/km) at metamorphic peak and a CW

AC

CE

PT

ED

MA

NU

P-T path of evolution.