Accepted Manuscript Title: Neoproterozoic crustal evolution in Sri Lanka: Insights from petrologic, geochemical and zircon U-Pb and Lu-Hf isotopic data and implications for Gondwana assembly Author: M. Santosh T. Tsunogae Sanjeewa P.K. Malaviarachchi Zeming Zhang Huixia Ding Li Tang P.L. Dharmapriya PII: DOI: Reference:
S0301-9268(14)00333-7 http://dx.doi.org/doi:10.1016/j.precamres.2014.09.017 PRECAM 4092
To appear in:
Precambrian Research
Received date: Revised date: Accepted date:
2-8-2014 13-9-2014 15-9-2014
Please cite this article as: Santosh, M., Tsunogae, T., Malaviarachchi, S.P.K., Zhang, Z., Ding, H., Tang, L., Dharmapriya, P.L.,Neoproterozoic crustal evolution in Sri Lanka: Insights from petrologic, geochemical and zircon U-Pb and Lu-Hf isotopic data and implications for Gondwana assembly, Precambrian Research (2014), http://dx.doi.org/10.1016/j.precamres.2014.09.017 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.
Page |1
Neoproterozoic crustal evolution in Sri Lanka: Insights from petrologic, geochemical and zircon U‐
ip t
Pb and Lu‐Hf isotopic data and implications for
cr
Gondwana assembly
us
an
M. Santosh1,2*, T. Tsunogae3,4, Sanjeewa P.K. Malaviarachchi5,6,
M
Zeming Zhang7, Huixia Ding7, Li Tang1, P.L. Dharmapriya6
1
School of Earth Science and Resources, China University of Geosciences Beijing, No. 29
d
Xueyuan Road, Haidian District, Beijing, 100083, China 2
780-8520, Japan 3
te
Division of Interdisciplinary Science, Faculty of Science, Kochi University, Akebono-cho, Kochi Faculty of Life and Environmental Sciences, University of Tsukuba, Ibaraki 305-8572, Japan
4
Ac ce p
Department of Geology, University of Johannesburg, Auckland Park 2006, South Africa
5
Department of Geology, Faculty of Science, University of Peradeniya, Peradeniya 20400, Sri Lanka 6 Postgraduate Institute of Science, University of Peradeniya, Peradeniya 20400, Sri Lanka 7
State Key Laboratory of Continental Tectonic and Dynamics, Institute of Geology, Academy of
Geological Science, No. 26 Baiwanzhuang Road, Beijing 100037, China
Corresponding author e-mail:
[email protected]
Page 1 of 111
Page |2
ip t
Research Highlights
Neoproterozoic arc magmatism prior to Gondwana assembly in Sri
cr
Lanka
us
Crustal growth through both juvenile and reworked components Highland Complex is an accretionary suture
an
Double-sided subduction in Neoproterozoic and collision in
te
Abstract
d
M
Cambrian
Ac ce p
Sri Lanka, the ‘pendant’ of Gondwana, is a collage of distinct crustal blocks that preserve important records of major Neoproterozoic tectonothermal events.
Here we present the
petrology, geochemistry, zircon U-Pb geochronology and Lu-Hf isotopes on a suite of metaigneous rocks including granodiorite, diorite and garnet amphibolite from the Kadugannawa Complex (KC), granodiorite from the Wanni Complex (WC) and mafic granulites, gabbros and garnet-bearing charnockite from the Highland Complex (HC) along a NW – SE transect. The regional metamorphic peak P-T conditions were estimated from garnet-clinopyroxeneplagioclase-quartz assemblage in the metagabbro as 830-860°C and 9.4-9.8 kbar. Slightly lower temperature ranges of 700-780°C were obtained from garnet amphibolite, metagranodiorite and metadiorite, corresponding retrograde conditions. Trace element and rare earth element patterns as well as Rb-Y-Nb and Rb-Yb-Ta discrimination plots show volcanic arc affinity for the granodiorite, diorite and garnet charnockite suggesting that the protoliths of the rocks were formed from felsic to intermediate arc magmas. The mafic granulites and magnesian metagabbro also suggest volcanic arc affinity and indicate subduction-related mafic magmatism and magma underplating. The garnet-bearing metagabbro shows N-MORB signature, whereas the garnet
Page 2 of 111
Page |3 amphibolite displays oceanic island alkali basalt affinity. These rocks therefore represent accretion of the remnants of oceanic lithosphere during the subduction-collisional event. Zircons in a metadiorite and the surrounding metagranodiorite from the KC yield ages of 980 ± 16 Ma to 916 ± 57 Ma marking early Neoproterozoic magmatism followed by metamorphism at 532 ± Zircons in the garnet amphibolite from this complex show extensive metamorphic
recrystallization yielding a weighted mean
206
ip t
18 Ma.
Pb/238U age of 520.7 ± 6.6 Ma. From the WC,
zircons in a metagranodiorite define three groups of weighted mean
206
Pb/238U ages at 805 ± 12
cr
Ma (emplacement of the magmatic protolith), 734.0 ± 4.6 Ma (Cryogenian thermal event) and 546.0 ± 5.7 Ma (latest Neoproterozoic-Cambrian metamorphism). Zircons from the HC record
us
multiple late Neoproterozoic – Cambrian thermal events with weighted mean
206
Pb/238U ages of
576.8 ± 9.3 Ma and 523.0 ± 7.1 Ma (metagabbro) and 579 ± 10 Ma, 540.4 ± 6.0 Ma and 511.1 ± 5.9 Ma (garnet charnockite), 553.0 ± 3.2 Ma (mafic granulite), and 539.1 ± 4.4 Ma (mafic granulite Lu-Hf data reveal dominantly positive εHf(t) values for zircons in the metadiorite and
an
sill).
metagranodiorite from the KC (-1.1 to 7.2) and Hf crustal model ages (TDMC) in the range of 1206 to 1733 Ma suggesting a mixed source from both juvenile and Paleo-Mesoproterozoic
M
components. However, zircons in the garnet amphibolite from this complex show dominantly negative εHf(t) values (mean -16.5) with TDMC in the range of 2356 to 2828 Ma suggesting reworked Neoarchean-Paleoproterozoic crustal source. Zircons in metagranodiorite of the WC
d
also possess negative εHf(t) values (mean -6.2) with TDMC in the range of 1799 to 2498 Ma
te
suggesting reworked Paleoproterozoic crust as the magma source. The HC rocks also preserve distinct imprints of reworking of older crust. The zircon εHf(t) values in mafic granulite show a tight
Ac ce p
cluster from -2.2 to 0.1 with TDMC in the range of 1501-1651 Ma suggesting a mixed source from both juvenile Neoproterozoic and reworked Mesoproterozoic components. Zircons in the metagabbro from this complex show negative εHf(t) values (mean -6.3) and TDMC of 1847-1978 Ma. Zircons in the garnet charnockite also display highly negative εHf(t) values (mean -17.7) and older TDMC) (mean 2614 Ma) suggesting reworked Paleoproterozoic crust. Zircons in the mafic granulite sample show negative εHf(t) values (mean -14.1) and TDMC between 2263 to 2790 Ma indicating that the source material for the magma evolved from the Neoarchean-Paleoproterozoic crust. In summary, the 916 to 980 Ma ages from the KC represent arc magmatism during early Neoproterozoic, followed by the 805 Ma granodioritic magma emplacement in the WC. Repeated thermal events during mid and late Neoproterozoic are also recorded from the Wanni and Highland Complexes, culminating in Cambrian high-grade metamorphism that reached ultrahightemperature conditions. We propose a model of double-sided subduction during the Neoproterozoic, where the Wanni Complex to the west and the Vijayan Complex to the east represent continental arcs, culminating in collision along the HC during late NeoproterozoicCambrian.
Page 3 of 111
Page |4
Key words: Petrology and geochemistry; Zircon U-Pb geochronology and Lu-Hf isotopes; Arc magmatism; Neoproterozoic; Sri Lanka.
ip t
1. Introduction
cr
Crustal evolution models on both Archean and Phanerozoic Earth identify the
us
role of arc magmas derived from both juvenile and reworked older components (e.g., Condie and Kröner, 2013; Condie and Aster, 2013; Santosh et al., 2013a; In convergent margins, continental growth occurs
an
Xiao and Santosh, 2014).
through vertical addition of magmas derived by the melting of downgoing oceanic
M
slabs, as well as lateral accretion of oceanic and trench material, together with other exotic oceanic and crustal fragments.
For instance, the Central Asian
te
d
Orogenic Belt (CAOB) representing the largest Phanerozoic accretionary orogen on the globe is a typical example of continental growth through lateral accretion
Ac ce p
of young arc complexes and old micro-continents together with vertical growth through underplating of mantle-derived magmas (e.g., Xiao and Santosh, 2014, and references therein). Many of the major Precambrian collisional orogenic belts were also constructed through prolonged subduction-accretion processes, similar to the modern analogue of the ongoing subduction and accretion in the Western Pacific (Santosh et al., 2009). Accretionary orogens also characterize the margins of the major supercontinents built through the amalgamation of continental fragments (Cawood and Buchan, 2007).
The geology of Sri Lanka provides important insights into continental growth in
Page 4 of 111
Page |5 the Neoproterozoic Earth. Tectonics of Sri Lanka has also been in focus in relation to the history of assembly of supercontinents, particularly because of its central position within the India – Madagascar – Africa – East Antarctica collage
ip t
of the late Neoproterozoic Gondwana supercontinent (Yoshida et al., 1992). The Sri Lankan Precambrian basement has been subdivided into four major terrains
cr
(e.g. Cooray, 1994): the Wanni Complex (WC) to the west together with the
us
Kadugannawa Complex (KC), the Highland Complex (HC) at the middle, and the Vijayan Complex (VC) to the east (Fig. 1), largely based on Nd- model age
an
mapping and regional structural interpretations (Milisenda et al., 1988, 1994; Kröner et al., 1987, 1991; Kehelpannala, 1991, 2003, 2004). The WC is
M
considered to represent a higher crustal level than that of the HC although there
d
is no clear structural break between the rocks of the two complexes, and the
te
contact between these two has been obliterated by later events (Voll and Kleinschrodt, 1991). The boundary between the VC and the HC is well defined as
Ac ce p
a thrust/shear contact (Hatherton et al., 1975; Vitanage, 1985; Voll and Kleinschrodt, 1991; Kriegsman, 1994). Three tectonic klippens of similar petrological and structural features to the HC are exposed in the south eastern part of Sri Lanka as inliers, namely Kataragama Klippe, Buttala Klippe and Kuda Oya Klippe. Some authors (e.g. Kehelpannala, 1991, 1997) have suggested that the KC forms part of the WC based on geochronological, geochemical and structural evidence. The HC domain is interpreted to be a part of a supracrustal basin developed in a suture zone with the Lützow-Holm Complex in East Antarctica during the final phase of Gondwana assembly (Shiraishi et al. 1994).
Page 5 of 111
Page |6 The WC, KC and VC domains are considered as Grenville-aged terranes of arcrelated settings at the outer margin of the Rodinia supercontinent (Kehelpannala, 2003, 2004; Kröner et al., 2003; Willbold et al., 2004; Kröner et al., 2013).
ip t
Although Sri Lanka is a relatively small continental fragment, mimicking a pendant on the heart of Gondwana, the several distinct tectonic units in this
cr
island and their history of amalgamation offer important clues for evaluating
us
Precambrian crustal evolution history.
an
In this study, we present results from our field investigations and petrologic, geochemical, zircon U-Pb and Lu-Hf isotopic studies on a suite of
Our data reveal Neoproterozoic crustal growth through protracted
d
Complexes.
M
metamorphosed magmatic rocks from the Wanni, Kadugannawa and Highland
te
arc magmatism derived from both juvenile and reworked sources, and lateral accretionary growth by the incorporation of oceanic and continental material.
Ac ce p
The subduction-accretion process culminated in continental collision and highgrade metamorphism during late Neoproterozoic – Cambrian.
2. Geological background
2.1 Overview of geology and petrology
The HC is dominantly composed of granulitic meta-quartzites, marbles, calc-silicates and metapelitic gneisses, in association with charnockites (Cooray,
Page 6 of 111
Page |7 1962, 1984; 1994; Perera, 1984; Mathavan et al. 1999; Mathavan and Fernando, 2001; Dharmapriya et al. 2014). The metasedimentary rocks such as marble and quartzite can be traced for more than 40 km in the central and northeastern part
ip t
of the HC, whereas marble and quartzite are scarce in the southwestern part of this complex (Mathavan et al. 1999). Bands of wollastonite-scapolite, diopside-
cr
and scapolite-bearing calc-granulites and cordierite-bearing gneisses occur in the
us
south and southwestern parts of the HC (Cooray, 1962, 1984; Hapuarachchi, 1968; Perera, 1984; Prame, 1991a; Mathavan et al., 1999; Mathavan and
an
Fernando, 2001).
Metamorphism under granulite facies conditions has been well established
M
in the HC from meta-basaltic and gabbroic to intermediate rocks (Sandiford et al.,
d
1988; Schumacher et al., 1990; Schumacher and Faulhabar, 1994), charnockites
te
(Prame, 1991b) and metamorphosed pelitic rocks (Perera, 1984, 1987; Prame, 1991a, Hiroi et al., 1994, Raase and Schenk, 1994). A P-T zonation across the
Ac ce p
HC has been established where pressures and temperatures decrease from 9-10 kbar and 830 °C in the east and southeast to 5-6 kbar and 700 °C in the northwest (Faulhaber and Raith, 1991; Schumacher and Faulhaber, 1994). Sajeev and Osanai (2005) identified a thermal gradient of 591 - 996 °C across south western to the central region of the HC. A clockwise P-T path is inferred for the pelitic rocks based on the sequence of kyanite and staurolite (occurring as inclusions in garnet) followed by sillimanite and andalusite (Perera, 1987; Hiroi et al., 1994, Raase and Schenk, 1994; Mathavan et al., 1999, Kriegsman and Schumacher, 1999). In contrast, reaction textures involving pyroxenes,
Page 7 of 111
Page |8 plagioclase, garnet and quartz in some meta-igneous rocks (Perera, 1987; Schumacher et al., 1990; Prame, 1991b) and high temperatures (>900 °C) from pyroxene exsolution textures (Schenk et al., 1988) have been used to suggest
ip t
isobaric cooling prior to the uplift. Although isobaric cooling is apparently not recorded in pelitic rocks of the HC apart from the report of Karunaratne et al.
cr
(2002), Dharmapriya et al. (2014) showed a rare evidence from a khondalite as
us
evidence for isobaric cooling from T of ~900 down to 770 °C at P of 7.5 kbar before uplift. Voll et al. (1994) derived peak temperatures of metamorphism at
an
850-900 °C using two-feldspar thermometry. Ultrahigh temperature (UHT) granulites have been reported in few localities in the central HC and rarely from
M
the southwestern part of the HC. These rocks witnessed extreme crustal
d
metamorphism at temperature of 925 °C to 1150 °C and pressures of 9-12.5
te
kbar (e.g. Kriegsman and Schumacher, 1999; Osanai et al., 2006; Sajeev and
Ac ce p
Osanai, 2004a, 2004b; Sajeev et al., 2007).
The WC occurs to the west-northwest of the HC and is dominantly
composed of upper amphibolite to granulite facies orthogneisses and minor meta-sediments.
The meta-igneous rocks show varied protolith-chemistry
ranging from granitic, granodioritic, monzonitic, tonalitic, charnockitic and enderbitic composition (e.g. Pohl and Emmermann, 1991). Garnet sillimanite gneisses, cordierite gneisses, quartzites and calc-silicate rocks occur as minor bands, particularly close to the inferred boundary with the HC. Metacarbonates are absent in the WC although apatite bearing unmetamorphosed carbonatite
Page 8 of 111
Page |9 body is exposed in the northern WC. The western part of the WC is mainly composed of less-deformed granites, including the unmetamorphosed posttectonic K-feldspar rich granite at Tonigala (e.g. Hölzl et al., 1991, Cooray, 1994).
ip t
Migmatization is also widespread within the unit. The peak metamorphic conditions of the WC are estimated to be 700-830 °C and 5-7 kbar (Schenk et al.,
cr
1991; Raase and Schenk, 1994). Classic outcrops of incipient charnockites
us
represent granulite formation within the host amphibolite facies gneisses in the central regions of the WC. In these zones, well-foliated light gray hornblende-
an
gneisses show veins and patches of dark, greenish-brown and coarse-grained incipient charnockite formed along mesoscopic shear zones and foliation planes,
M
representing pervasive fluid influx in the Sri Lankan lower crust (e.g. Hanson et
d
al., 1987; Burton and O’Nions, 1990; Baur et al., 1991; Milisenda et al., 1991;
te
Santosh et al., 1991; Kehelpannala, 1999; Perchuk et al., 2000).
Ac ce p
The KC, often referred as ‘Arenas’ (Vitanage, 1972; Almond, 1991), a
synonym used in Sri Lanka for doubly plunging synforms, is exposed in the northwestern part of Kandy. The complex includes upper amphibolite to granulite facies basement rocks (Kröner et al., 1991; Cooray, 1994; Malaviarachchi and Takasu, 2005, 2011a). Perera (1983) argued a sedimentary origin of the wellbanded hornblende biotite gneisses in the KC, whereas Cooray (1984) regarded their possible protoliths as basaltic lavas or volcanic ash. Relics of primary magmatic layering are preserved in these synforms mainly composed of hornblende gneisses, amphibolites and minor pyroxenites (Kröner et al., 1991).
Page 9 of 111
P a g e | 10 Hornblende- and biotite-bearing orthogneissic rocks with gabbroic, dioritic and trondhjemitic rocks of calc-alkaline affinity (Pohl and Emmermann, 1991) are interlayered within granodioritic to granitic gneisses, charnockites, enderbites and
ip t
minor metasediments (Kröner et al., 2003). Minor layered mafic and ultramafic rocks have also been reported (Stosch, 1991). The KC is considered by some
cr
authors as a migmatized and metamorphosed equivalent of a layered mafic-
us
ultramafic complex, named the ‘Kandy Layered Intrusion’ (KLI: Vol and Kleinschrodt, 1991; Stosch, 1991; Kleinschrodt et al., 1991).
an
The VC exposed in eastern Sri Lanka consists predominantly of metaigneous rocks including upper amphibolite facies calc-alkaline suite of granitoid
M
gneisses, augen-gneisses and minor amphibolites possibly derived from mafic
d
dykes. This granitoid suite was interpreted as a product of subduction-related
te
magmatism (Milisenda, 1991; Pohl and Emmermann, 1991). However, in contrary to previous interpretations, Kröner et al. (2013) reported that the Vijayan
Ac ce p
gneisses show evidence for granulite-facies metamorphism. Minor sedimentary enclaves such as metaquartzite, calc-silicate rocks and marble also occur within the VC (Dahanayake, 1982; Dahanayake and Jayasena, 1983; Kröner et al., 1987,
1991).
The
granitoid
gneisses
show
TTG
(tonalite-trondhjemite-
granodiorite)-dominated signature (Milisenda, 1991; Pohl and Emmermann, 1991) whereas the hornblende gneisses display compositions between diorite and leucogranite with the dominant composition being granodiorite and granite (Kröner et al., 1991; Milisenda, 1991; Kleinschrodt, 1994; Milisenda et al., 1994). Minor metagabbroic rocks and mafic dykes occur interlayered and infolded with
Page 10 of 111
P a g e | 11 the granitoid gneisses. Many exposures of the Vijayan gneisses show evidence of migmatization. Prame (1997) noted that many of the granitic rocks in the VC have A-type characteristics and inferred their formation through melting of felsic
ip t
crust. Igneous charnockite bodies occurring at the SE margin of the VC are interpreted as tectonically infolded or inter-sliced fragments of rocks of the HC
cr
similar to those of the HC-klippens present in the VC, during the thrusting of the
us
HC over the VC (Kröner et. al. 1991).
an
2.2 Summary of previous geochronological studies and tectonic models
d
M
2.2.1 Geochronology
te
The Nd model ages of 3.0-2.2 Ga reported from the HC (Milisenda et al., 1988, 1994) suggest that the rocks in this complex were derived from late
Ac ce p
Archean sources. Hölzl et al. (1991, 1994) and Köhler et al. (1991) presented Rb-Sr whole rock, Rb-Sr biotite and Sm-Nd garnet indicating ages of ~2 Ga. Cordani and Cooray (1990) reported Rb-Sr isochron age of 1100 Ma and De Maesschalck et al. (1990) interpreted a Rb-Sr whole rock isochron age of 1930 Ma as the timing of high-grade metamorphism. Rb-Sr and Sm-Nd isochron ages of 2.6-2.3 Ga and 500-450 Ma were also reported from the HC (Kagami et al. 1990). Garnet-whole rock Sm-Nd ages for widely separated HC paragneisses and metabasites show metamorphic ages of ca. 600 Ma (Hölzl et al., 1991; Milisenda, 1991). Ion microprobe U-Pb study on zircons by Kröner et al. (1987)
Page 11 of 111
P a g e | 12 documented Archean to early Proterozoic ages (3.2-2.4 Ga) from detrital zircons in the HC metasediments. They also reported an intrusion age of 1.1 Ga and new zircon growth at about 550 Ma from meta-igneous rocks. Köhler et al. (1991)
ip t
reported U-Pb zircon and monazite data indicating depositional ages of ~2 Ga and ~550 Ma ages for the high-grade metamorphism. Baur et al. (1991) reported
cr
U-Pb zircon crystallization ages of 1940 Ma and a Pb loss event at 560-550 Ma
us
indicating high-grade metamorphism. Kröner and Williams (1993) reported SHRIMP zircon ages for granitic gneiss from the HC as 1.9 Ga for the
an
crystallization of the igneous protolith and 531 Ma as the age of metamorphism. The lower intercept ages in concordia plots fall between 550-610 Ma for
M
discordant single zircon and zircon fractions in the granulites from the HC, whereas multigrain fractions of metamorphic zircon, monazite and rutile gave Pb/206Pb ages of 541-608 Ma (Baur et al., 1991, Hölzl et
d
207
te
concordant U-Pb and
al., 1991, 1994; Kröner et al., 1994). Burton and O’Nions (1990) analysed U and
Ac ce p
Pb isotopes in ilmenite giving an isochron age of 1.1 Ga which they interpreted as the metamorphic age. Hölzl et al. (1994) used U-Pb method of zircon and monazite from the same rock and reported ages of 555-592 Ma. In a recent study, Malaviarachchi and Takasu (2011b) obtained ca. 728-460 Ma ages for metamorphism using electron microprobe monazite CHIME method. Sajeev et al. (2003) obtained the Sm-Nd isochron whole-rock age of ca. 1.5 Ga from the UHT granulites of the HC and inferred that the granulites could be relics of a pre-Pan African metamorphism. In another study, Sajeev et al. (2007) reported U-Pb zircon metamorphic age of ca. 580 Ma from relatively
Page 12 of 111
P a g e | 13 higher pressure mafic granulites which is in agreement with zircon and monazite ages of ca. 570 Ma, determined from sapphirine granulites of the central HC
HC with the assembly of the Gondwana supercontinent.
ip t
(Sajeev et al. 2010), correlating the ultrahigh-temperature metamorphism in the
From the WC, Nd-modal ages of ca. 700 Ma -1.1 Ga have been reported,
cr
with multiple magmatic events as documented by zircon ages of 1100-530 Ma
us
(Burton and O’Nions, 1990; Kröner et al., 1991, 2003). Massive charnockites from the northern part of the WC have also older protolith emplacement ages as
an
from zircon evaporation (Kröner et al. 1994). Deposition of detrital zircons in metapelites and the emplacement of magmatic zircons of charnockitic and
M
enderbitic rocks yield ion microprobe U-Pb concordia upper intercept age of ~1.1
d
Ga considered as the timing of emplacement of the precursor magmas. The
te
lower intercept ages range between 540-590 Ma correlate with the high-grade metamorphism (e.g. Hölzl et al., 1991). Typical arrested (in-situ) charnockites and
Ac ce p
its host granitic/tonalitic gneiss from the Udadigana quarry in Kurunegala show ages on a common regression line with concordia intercepts at 771 Ma and 563 Ma, considered as the intrusion age and the metamorphic age respectively (Baur et al. 1991). Burton and O’Nions (1990) dated zircons from charnockitized patches and inferred that the in-situ charnockitization in the WC has occurred at 535 Ma. Thus, the arrested or in-situ charnockitization of the WC is considered to post-date the peak granulite facies metamorphism, and principally controlled by fluid influx along shear zones (Kehelpannala, 1999). Nd-model ages of the rocks from the KC range between ca. 900 Ma and
Page 13 of 111
P a g e | 14 1.6 Ga (Milisenda et al., 1988, 1994). Pb-Pb zircon evaporation ages from the KC show ca. 770 Ma -1.1 Ga indicating multiple calc-alkaline magmatic events covering a period of approximately 300 million years (Kröner et al., 2003; Willbold
ip t
et al., 2004). Most of the dioritic to granodioritic rocks of the KC display ca. 900 Ma ages. Xenomorphic zircons of gabbroic gneisses show relatively younger
cr
ages of ~700 Ma whereas idiomorphic zircons from the same sample gave ages
us
of ca. 900 Ma. Perera and Kagami (2001) argued that the regional high grade
isotope systematics on charnockitic rocks.
an
metamorphism of the KC must be older than 1100 Ma based on a study of Nd-Sr
Whole-rock Nd isotope analyses from widely separated localities show
M
mean crustal residence ages of 1.0–1.8 Ga for the VC rocks. Sr isotope data are
d
in line with the Nd isotope results suggesting that the majority of the VC rocks
te
were juvenile and derived from a typical LREE-depleted, mantle-derived precursor and not from anatexis of older crustal lithologies (Milisenda et al.,
Ac ce p
1988, 1994; Kröner et al., 1991). Rb–Sr whole-rock isochron ages of ca. 800 Ma reported by De Maesschalck et al. (1990) were ascribed to metamorphic resetting whereas the magmatic zircons reflect emplacement of the protoliths of the Vijayan gneisses at 1000–1100 Ma (Kröner et al., 2003, 2013). Zircon U-Pb ages of ca. 591- 456 Ma (Hölzl et al., 1991; Kröner et al., 1987, 1991, 2013) are now regarded as the age of metamorphism of the VC.
2.2.2. Structural and tectonic models
Page 14 of 111
P a g e | 15 The HC granulites have been subjected to a major folding and thrusting event (e.g. Berger and Jayasinghe, 1976; Kriegsman, 1994; Kleinschrodt, 1994). The large folds with amplitudes of about 7-10 km and an exposed length parallel
ip t
to their axes of up to 50 km in the central Highlands formed after the peak of metamorphism (800-850 °C, 8-9 kbar) at slightly lower temperature (700-750 °C)
cr
(Kleinschrodt, 1994). These large upright folds were generated during the last
us
deformation event following high-grade metamorphism (e.g. Berger and Jayasinghe, 1976; Kriegsman, 1994). During this event, horizontal shorting and
an
vertical thickening of crust took place with unroofing of the crustal rocks into upper level (Kleinschrodt, 1994; Kriegsman, 1994). Hiroi et al. (2014) reported
M
evidence for fast exhumation of lower-crustal rocks to andalusite-stable upper-
d
crustal conditions by channel flow in a continental collision orogeny, based on
of the HC.
te
felsic/granitic inclusions within high grade ordinary granulites and UHT granulites
Ac ce p
Transposition of structures in the WC is less severe, compared to those in
the HC. In ductile shear zones, gneissic foliations are severely flattened resulting in the obliteration of the previous intrusive and deformation histories. However, in some regions, multiple phases of gabbroic to granitoid intrusions and their deformation are preserved in multiple foliations (Kriegsman, 1991; Kröner et al., 1994). Macroscopic folds of the open to tight type with gently plunging hinges are well developed in the WC similar to those in the HC and their axial traces continue for several kilometers (Tani and Yoshida, 1996). The contact between the WC and the HC is still controversial due to absence of any structural
Page 15 of 111
P a g e | 16 discontinuity in between. The combined HC – WC unit is considered to represent a tilted section of the former lower–middle crust, with the WC representing the
ip t
higher level.
The entire rock suite of the KC is folded into large-scale doubly plunging
cr
five NW-trending synforms (Arenas). These macroscopic plunging folds with
us
upright axial surfaces were generated after the formation of the major compositional layering. The interior domains of the KC preserve prograde
an
amphibolite facies conditions compared to the surrounding HC (Perera, 1983; Sandiford et al., 1988; Voll and Kleinschrodt, 1991). Hence, the KC is considered
M
to be at a higher crustal level than the HC (Sandiford et al., 1988; Kriegsman,
d
1994) whereas Kehelpannala (1997) interpreted that the KC represents a deeper
te
crustal level compared to the WC. Although the boundary between the HC and the KC is well defined as a structural break, that between the WC and the KC is
Ac ce p
not clear e.g. Kröner et al., 2003).
Kehelpannala (2003) noted that the homogeneous to compositionally well-
layered gneisses experienced strong ductile deformation during more than one phase of folding. The boundary between the VC and the HC is considered to be a sub-horizontal ductile shear zone/thrust (Hatherton et al., 1975; Vitanage, 1985; Voll and Kleinschrodt, 1991; Kriegsman, 1995), formed when the HC collided with the VC and produced late upright folds in the former and ductile deformation in the latter (Kehelpannala, 1997, 2003, 2004). Shear sense indicators and sheath folds found at the suture between the VC and the HC show
Page 16 of 111
P a g e | 17 that the HC, namely the hanging wall, apparently moved from N to S over the VC (Kehelpannala,
2003)
accompanying
widespread
mylonitization
and
ip t
migmatization along the shear zones (Kleinschrodt, 1994).
cr
2.3 Present study area and description of sampling locations
us
The present study focuses on petrology, geochemistry and geochronology of meta-igneous rocks collected from a traverse running from the WC through the
an
KC into the HC in a direction of increasing elevation (Fig. 1). Fresh samples were chosen from quarries and road-cut outcrops along this traverse. The NW-SE
M
sampling in our study is oblique to the general trend of foliation so that a
d
meaningful stratigraphic section is traversed.
te
In this paper, we present analysis of three samples from the KC, four samples from the HC, and one sample from the WC (Fig. 2). A summary of the
Ac ce p
location and salient features of the samples analyzed for geochemistry and geochronology are given in Table 1. Samples of the KC represent the Dumbara Synform (KDG-1), Huluganga Synform (KDG-2) and Metibokka Synform (KDG4). The HC samples (HC-1, HC-7, HC-8, HC-10 and HC-11) were selected from the areas which represent strong ductile deformation and relatively higher grade metamorphism in the HC. The sample of the WC (WC-2) comes from a classical outcrop of in-situ/incipient charnockitization at Udadigana in Kurunegala. Representative field photographs are shown in Figures 3 and 4.
Page 17 of 111
P a g e | 18 2.3.1 Kadugannawa Complex
ip t
Metagranodiorite (KDG-4) and metadiorite (KDG-1)
The metagranodiorite sample KDG-4 was collected from a quarry at Maussawa
cr
along Harankahawa Road in Galagedara. The quarry exposes evidence for mafic
us
magma intrusion into a felsic magma chamber and subsequent mixing, mingling and fractionation. The presence of mafic-felsic layers in the outcrop may indicate
an
lower crustal melting in a granodiorite-granite magma chamber where the restite portion is represented by the hornblendite-rich mafic layers. The felsic magma
M
chamber was intruded by mafic/intermediate (dioritic) magma as dykes/sills
d
resulting in magma mixing and mingling as evidenced from the mafic magmatic
te
enclaves (MMEs) of various sizes within the granodiorite. The rocks display intense deformation as seen by refolded isoclinal folds and eye-like sheath or
Ac ce p
doubly plunging folds while preserving the magmatic history.
The metadiorite KDG-1 was collected from a road side exposure at Kundasale (about 5 km from Thennekumbura) along the Kandy-Mahiyangana road. The outcrop lies at the lower south-western part of the doubly plunging Dumbara synform of the KC. The rock shows strongly flattened and stretched compositional layers, and bedding parallel with the surrounding TTG gneisses that show compositional layering. Cross cutting dykes of aplitic and granitic composition are also common in the highly deformed parts of the metadiorite.
Page 18 of 111
P a g e | 19
Garnet amphibolite (KDG-2)
ip t
This sample was collected from a newly excavated road cut at the Mahaberiyatenna Airport site. The rock occurs as a several tens of meter thick-
cr
layer within migmatitic hornblende-biotite gneiss in the south western margin of
us
the doubly plunging Huluganga synform of the KC. Coarse-grained calcic amphibole and porphyroblastic red-eye garnets are visible in the outcrop, where garnets
show
breakdown
textures
forming
an
some
plagioclase-bearing
symplectites. Hornblende-rich cumulate layers also occur within the magmatic
d
M
layers associated with the amphibolite.
te
2.3.2. Wanni Complex
Ac ce p
Metagranodiorite (WC-2)
The sample was collected from an abandoned quarry at Udadigana in Kurunegala, tectonically located on the western limb of the Kurunegala Synform of the WC. The quarry exposes excellent evidence for in-situ/incipient charnockitization with greenish patches and veins of orthopyroxene-bearing anhydrous zones overprinting the host gneiss. The sample collected for this study is from the tonalitic to granodioritic host rock of the incipient charnockite. A gradational transition of the host gneiss to charnockite is evidenced at the quarry
Page 19 of 111
P a g e | 20 where some parts are totally converted into a coarse-grained charnockite, suggesting three dimensional networks of fluid- filled veins or tubes.
ip t
2.3.3. Highland Complex
us
cr
Mafic granulite (HC-1, HC-11)
Mafic granulite bands, boudins and enclaves of variable thickness (0.1 to 3m)
an
occur within marbles at the Digana quarry. Our sample HC-1 is from one of these bands. The marble is composed of dolomitic and calc-silicate layers parallel to
M
the regional compositional layering of the HC. Mafic granulite enclaves are
d
located parallel to bedding, and are fragmented with angular edges. They show
te
strong deformation including evidence for flattening and stretching of minerals like pyroxene and hornblende. In some parts of the quarry, the mafic granulite
Ac ce p
shows boudinaged features.
The sample HC 11 is from a mafic sill along a roadside exposure ~2 km from Kivullinda Junction towards Randenigala, along the Kandy-Randenigala road. This location is close to the hinge on the west limb of the Digana Antiform of the HC. Several metadioritic sills ranging in thickness from few tens of cm to several meters occur within alternating layers (5 to 40 cm thickness) of quartzite and psammitic gneiss, parallel to bedding with sharp contacts. Some of the mafic sills contain garnet.
Page 20 of 111
P a g e | 21
Metagabbro (HC-7 and HC-8)
ip t
The samples of metagabbro were collected from road side outcrops along the Kandy-Hantana Road. These rocks occur as discontinuous lenses (few tens of
cr
metres thick) within quartzite layers of approx. 2 km thickness and biotite gneiss
us
(partly charnockitic at places), close to the hinge on the western limb of the Kandy Antiform of the HC. Reddish porphyroblastic garnets are extensively
an
distributed in the rock, partly reflecting a cumulate-textured mafic granulite. Minor quartz-bearing domains (sample HC-8) occur within the metagabbroic host rock
te
Garnet charnockite (HC-10)
d
M
(HC-7).
Ac ce p
The garnet- and orthopyroxene- bearing charnockite sample was collected from a road outcrop near Hakurutale, between 12 and 13 km along the KandyRandenigala Road. The 0.5 – 1 km thick layer is interbanded between a granitic gneiss and a biotite gneiss layer. In some domains, coarse deep red almandinerich garnets occur and some of the garnet grains show plagioclase-bearing symplectites.
Symplectitic
garnets
and
coarse
grained
dark
brownish
hypersthene grains occur along the foliation plane. The rock is generally medium to coarse grained showing brownish to gray-green greasy charnockitic luster.
Page 21 of 111
P a g e | 22 3. Analytical techniques
ip t
3.1 Petrography and mineral chemistry
Polished thin sections were prepared for petrographic study at the State Key
cr
Laboratory of Geological Processes of Mineral Resources of China University of
us
Geosciences Beijing, and at the University of Tsukuba, Japan. Mineral chemical analyses were carried out using an electron microprobe analyzer (JEOL
an
JXA8530F) at the Chemical Analysis Division of the Research Facility Center for Science and Technology, the University of Tsukuba. The analyses were
M
performed under conditions of 15 kV accelerating voltage and 10 nA sample
te
supplied by JEOL.
d
current, and the data were regressed using an oxide-ZAF correction program
Ac ce p
3.2 Whole-rock geochemistry
Representative samples were chosen for whole-rock geochemical analyses from the WC, HC, and KC which are devoid of surface alteration or weathering. The size of the samples were initially reduced in a jaw crusher, and then manually fine-powdered in an agate mortar. Major oxides were analyzed by Thermo JarrellAsh ENVIRO II ICP, and minor and trace elements by Perkin Elmer Sciex ELAN 9000 ICP/MS at Actlabs, Canada, using lithium metaborate/tetraborate fusion technique.
Page 22 of 111
P a g e | 23
3.3 Zircon U-Pb and Lu-Hf isotopes
ip t
Zircon grains were separated by gravimetric and magnetic separation from crushed rock samples, and then purified by hand picking under a binocular
cr
microscope at the Yu’neng Geological and Mineral Separation Survey Centre,
us
Langfang City, Hebei Province, China. The grains were mounted in epoxy resin discs and polished to reveal mid-sections, followed by gold sputter coating.
an
Zircons were imaged under both transmitted and reflected light, and were imaged using cathodoluminescence (CL) to identify internal structures and choose
M
potential target sites for U-Pb analyses. The CL imaging was carried out at the Beijing Geoanalysis Centre. Individual grains were mounted along with the Pb/238U age of 417 Ma (Black et al., 2003), onto
d
206
te
standard TEMORA 1, with
double-sided adhesive tape and enclosed in epoxy resin discs. The discs were
Ac ce p
polished to a mid depth and gold coated for cathodoluminescence (CL) imaging and U-Pb isotope analysis. Zircon morphology, inner structure and texture were examined by using a JSM-6510 Scanning Electron Microscope (SEM) equipped with a backscatter probe and Chroma CL probe. The zircon grains were also examined under transmitted and reflected light images using a petrological microscope.
Zircon U-Pb dating and trace element analysis were simultaneously conducted by LA-ICP-MS at the Stage Key Laboratory of Geological Processes
Page 23 of 111
P a g e | 24 and Mineral Resources, China University of Geosciences. Detailed operating conditions for the laser ablation system and the ICP-MS instrument are as described by Liu et al. (2008, 2010). Laser sampling was conducted using an
ip t
Excimer 193 nm GeoLas 2005 System with a spot size of 32 μm. An Agilent 7500a ICP-MS instrument was used to acquire ion-signal intensities. Nitrogen
cr
was introduced into the central gas flow (ArtHe) of the Ar plasma in the LA-ICP-
us
MS analysis, which increases the sensitivity for most elements by a factor of ~2 compared with the results without adding nitrogen (Hu et al., 2008). To keep time-
an
dependent elemental fractionation at a low level, a laser frequency of 4 Hz and a laser energy of 60 mJ were applied. The uncertainty of analysis is within ~1% for
M
the zircon standard.
d
Zircon 91500 was used as external standard for U-Pb dating, and was
te
analyzed twice every five analyses. Time-dependent drifts of U-Th-Pb isotopic ratios were corrected using a linear interpolation (with time) for every five
Ac ce p
analyses according to the variations of 91500 (i.e. two zircon 91500 + five samples + two zircon 91500). Preferred U-Th-Pb isotopic ratios used for 91500 are from Wiedenbeck et al. (1995). Uncertainty in the preferred values for the external standard 91500 was propagated to the ultimate results for the samples. Concordia diagrams and weighted mean calculations were made using Isoplot/Ex_ver3 (Ludwig, 2003). Trace element compositions of zircons were calibrated against multiple-reference materials (BCR-2G and BIR-1G) combined with internal standardization (Liu et al., 2010). The preferred values of element concentrations for the USGS reference glasses used are from the GeoReM
Page 24 of 111
P a g e | 25 database (http://georem.mpch-mainz.gwdg.de/). In situ Hf isotope compositions of zircon were determined using a Neptune Plus multicollector (MC)-ICP-MS system (Thermo Fisher Scientific, Germany)
ip t
and a Geolas 2005 Excimer ArF laser ablation system (Lambda Physik, Götingen, Germany) at the State Key Laboratory of Geological Processes and
cr
Mineral Resources, China University of Geosciences, Wuhan. The energy
us
density of laser ablation used was 5.3 J cm-2. Helium was used as the carrier gas in the ablation cell and was merged with argon (makeup gas) after the ablation
an
cell. All data for zircon in this study were acquired in the single-spot ablation mode at a spot size of 32 mm. Each measurement consisted of 20 s of
M
acquisition of the background signal followed by 50 s of ablation signal
d
acquisition. The operating conditions for the laser ablation system and the MC-
te
ICP-MS instrument, and the analytical method, are the same as those described
Ac ce p
in detail by Hu et al. (2012).
4. Results
4. 1 Petrology and mineral chemistry
4.1.1 Petrography
A summary of the petrographic features of each of the rock types analyzed in this study is given below. Representative photomicrographs are shown in Fig. 5.
Page 25 of 111
P a g e | 26
Metagranodiorite and metadiorite
ip t
The metagranodiorite and metadiorite in our study are generally coarse-grained with quartz, plagioclase, biotite and calcic amphibole (hornblende) as the
cr
dominant minerals. Sample WC-2, a typical metagranodiorite, is composed of
us
plagioclase (35-45%), K-feldspar (20-30%), quartz (20-25%), biotite (5-10%), and hornblende (2-5%), with accessory Fe-Ti oxides (1-2%), zircon, and apatite (Fig.
an
5a). Plagioclase and K-feldspar are coarse-grained (grain size range of 0.3-1.2 mm and 0.4-1.3 mm, respectively) showing granoblastic texture. Fine, rounded
M
grains of quartz (0.2-0.7 mm) are present along their grain boundaries, although
d
elongated coarse-grained (~4.2 mm) quartz also occurs aligned along the rock
te
foliation. Myrmekite textures can be commonly observed in some of the samples.
Ac ce p
The metagranodiorite (sample KDG-4) comprises plagioclase (30-40%), Kfeldspar (15-25%), hornblende (15-25%), biotite (10-15%), quartz (10-15%), and magnetite (1-2%). Weak foliation of the rock is defined by aligned flakes of hornblende (0.2-1.9 mm) and biotite (0.2-1.8 mm). Plagioclase and K-feldspar are coarse-grained (0.2-1.2 mm), sub-equigranular, and granoblastic (Fig. 5b). Quartz is relatively fine grained (0.2-0.8 mm) compared to that of the sample WC2.
The sample KDG-1 is an orthopyroxene-bearing and K-feldspar-free metadiorite
Page 26 of 111
P a g e | 27 showing textures similar to those of the sample KDG-4. The mineral assemblage is plagioclase (50-60%), hornblende (15-25%), biotite (10-15%), quartz (5-10%), and ilmenite (2-3%), with accessory orthopyroxene, pyrite, apatite, and zircon
ip t
(Fig. 5c). The rock does not have a distinct foliation, and is characterized by granoblastic and semi-equigranular texture of plagioclase (0.6-2.2 mm) and
cr
quartz (0.2-0.8 mm). Both hornblende (0.5-1.5 mm) and biotite (0.2-0.8 mm) are
us
subhedral and randomly oriented along the grain boundaries of plagioclase and quartz. Orthopyroxene is also fine-grained (0.2-0.4 mm) and intergrown with
an
biotite or hornblende.
d
M
Garnet amphibolite
te
Sample KDG-2 is composed of very coarse-grained (~9 mm) poikiloblastic garnet (5-15%) in the matrix of aggregates of hornblende (20-30%) and plagioclase (30-
Ac ce p
40%) (Fig. 5d). Foliation is not obvious in thin section. Accessory minerals are ilmenite-magnetite intergrowth (2-5%), apatite (2-5%) and zircon, and quartz is present only as rare inclusions in garnet. The garnet contains numerous inclusions of ilmenite, apatite, hornblende, and quartz. Hornblende (0.3-2.4 mm) and plagioclase (0.5-1.7 mm) are coarse grained, but is fine-grained compared to the garnet. The garnet is surrounded by medium-grained intergrowth of hornblende + plagioclase, suggesting the progress of the following CFMASH continuous reaction (1); Grt + Qtz + H2O => Hbl + Pl
(1)
Page 27 of 111
P a g e | 28
ip t
Mafic granulite
The sample HC-1 is a medium-grained garnet-free mafic granulite that is
cr
composed of plagioclase (30-40%), clinopyroxene (15-25%), hornblende (15-
us
25%), orthopyroxene (10-20%), ilmenite (2-5%) and magnetite (<1%) (Fig. 5e). Calcite occurs as an accessory mineral, in accordance with the field occurrence
an
of the rock as enclave within marbles. The rock shows a weak foliation defined by aligned aggregates of ferromagnesian minerals. All the dominant minerals are
d
M
medium to fine grained (0.2-0.8 mm).
te
Sample HC-11 is also a garnet-free mafic granulite without any marked foliation and composed of plagioclase (30-40%), clinopyroxene (20-30%), orthopyroxene
Ac ce p
(10-20%), hornblende (10-20%), and accessory ilmenite and zircon (Fig. 5f) in an equigranular matrix (0.2-0.8 mm). Fine-grained aggregates of orthopyroxene and plagioclase occur as grayish spots of about 4 mm in diameter, typically associated with quartz (Fig. 5g), suggesting the progress of the following decomposition reaction of garnet: Grt + Qtz => Opx + Pl
(2)
This reaction involved the near total consumption of garnet possibly associated with post-peak decompression (e.g., Harley, 1989).
Page 28 of 111
P a g e | 29 Metagabbro
The sample HC-7 is a coarse-grained garnet-bearing metagabbro with variable
ip t
quartz-bearing and absent domains. The rock is composed of clinopyroxene (4050%), garnet (30-40%), orthopyroxene (5-10%), and plagioclase (5-10%), with
cr
accessory Fe-Ti oxide (2-5%) and hornblende (1-2%) (Fig. 5h). No obvious
us
foliation can be seen. Garnet and clinopyroxene are coarse grained (~1.7 mm) and semi-equigranular. Garnet is subhedral, and almost free from inclusions
an
except for minor apatite, hornblende, ilmenite, and quartz. Clinopyroxene is also subhedral, and contains thin exsolution lamellae of orthopyroxene. Plagioclase
M
(0.2-0.7 mm) is anhedral and fills the matrix of garnet and pyroxenes. Fine-
d
grained and rounded Fe-Ti oxide (intergrowth of ilmenite and magnetite) grains
te
are scattered throughout the thin section. These textures are similar to metagabbros reported from the Palghat-Cauvery Suture Zone in Southern India
Ac ce p
(e.g., Sajeev et al., 2009; Santosh et al., 2010; Saitoh et al., 2011; Koizumi et al., 2014).
Sample HC-8 is a quartz-rich domain of sample HC-7. It comprises plagioclase (30-40%), clinopyroxene (30-40%), garnet (10-20%), quartz (5-10%), and orthopyroxene (2-5%) with accessory ilmenite, magnetite, rutile, and zircon. Foliation is defined by elongated quartz grains (~4 mm) and garnet aggregates. Garnet is coarse-grained (0.3-3.3 mm), subhedral to anhedral, and often surrounded by fine-grained (~0.2 mm) symplectitic intergrowth of vermicular
Page 29 of 111
P a g e | 30 orthopyroxene and plagioclase or thin corona of orthopyroxene + plagioclase (Fig. 5i). Such intergrowths can be readily observed in quartz-rich portion of the rock, confirming the progress of reaction (2) probably related to rapid
ip t
decompression. Clinopyroxene (0.2-1.1 mm) is subhedral, and often occur as aggregates elongated along the foliation particularly at plagioclase-rich portion of
cr
the rock. Ilmenite occurs as fine-grained (<0.3 mm) rounded to sub-rounded
us
grains, or irregular mineral intergrowing with clinopyroxene.
an
Garnet charnockite
M
Sample HC-10 is characterized by coarse-grained (~2.5 mm) porphyroblastic
d
garnet in hand specimen (Fig. 3), although the modal abundance varies from 0 to
te
5 %. The examined sample corresponds to a garnet-absent portion of the rock, and comprises plagioclase (30-40%), quartz (10-20%), clinopyroxene (10-15%),
Ac ce p
orthopyroxene (5-10%), hornblende (5-10%), with accessory ilmenite, pyrite, calcite, zircon, and retrograde biotite (Fig. 5j). K-feldspar is absent in the rock, which is consistent with its dioritic composition as discussed later. The rock is also characterized by ribbon quartz elongated parallel to the rock foliation defined by thin alternation of pyroxene-rich and plagioclase-rich layers. Plagioclase is coarse-grained (0.7-2.5 mm), granoblastic, and semi-equigranular in plagioclaserich layers, while it is fine-grained (0.2-0.7 mm), subhedral, and locally intergrowing with fine-grained orthopyroxene (~0.2 mm) in pyroxene-rich layers as a product of reaction (2) (Fig. 5j). The lack of garnet in the sample probably
Page 30 of 111
P a g e | 31 implies complete consumption of the mineral during decompression. Mediumgrained (0.2-0.8 mm) and subhedral ortho- and clinopyroxenes, and fine-grained (0.1-0.5 mm) hornblende are dominant minerals in the pyroxene-rich layers,
ip t
whereas their grain size decreases in plagioclase-rich layers. Individual pyroxene and plagioclase grains show no obvious alignment, whereas hornblende is
cr
weakly oriented along the rock foliation. Minor biotite mantles pyroxene and
an
us
hornblende grains as a product of post-peak retrograde hydration.
M
4.1. 2 Mineral chemistry
d
Representative compositions of minerals from electron microprobe analyses are
Ac ce p
below.
te
given in Supplementary Tables 2 to 4, plotted in Figure 6, and briefly discussed
Garnet
Garnet in the metagabbros (HC-7, HC-8) and garnet amphibolite (KDG-2) is essentially almandine-rich with minor pyrope, grossular, and spessartine contents whose abundance varies slightly depending on the assemblage and core-rim positions (Supplementary Table 2, Figs. 6a,b). Garnet in sample HC-7 shows the highest XMg (= Mg/(Mg+Fe)) content of 0.24-0.28. Its almandine content increases slightly from core (Alm57-58 Prp22-23 Grs18-19 Sps1-2) toward rim (Alm59-60
Page 31 of 111
P a g e | 32 Prp19-20 Grs18-19 Sps1-2), which is a typical zoning pattern for garnet grains in mafic granulites. Garnet in sample HC-8 shows nearly consistent core-rim variation, but the core is slightly grossular rich as Alm59-60 Prp14-15 Grs24-25 Sps1-2 (core) and
ip t
Alm64-65 Prp13-14 Grs20-21 Sps1-2 (rim). In contrast, garnet in sample KDG-2 shows slightly Mn-rich rim (Alm58-59 Prp15-16 Grs19-20 Sps5-6) than core (Alm60-61 Prp15-16
us
cr
Grs20-21 Sps3-4), although their XMg values are nearly consistent (0.20-0.23).
an
Pyroxenes
Clinopyroxene in the mafic suites of present study is Mg-rich (XMg = 0.53-0.82),
M
and all are classified as augite (Fig. 6c), although their XMg varies depending on
d
samples. Clinopyroxene in the metagabbros (samples HC-7 and HC-8) shows
te
lower XMg of 0.61-0.64 and 0.53-0.59, than that in the mafic granulite (samples HC-1 and HC-11, XMg = 0.69-0.82). Single grains do not display any
Ac ce p
compositional zoning.
Orthopyroxene also shows a wide compositional range of XMg = 0.39-0.72. Orthopyroxene in mafic granulite (sample HC-1) shows the highest XMg of 0.710.72. Orthopyroxene in sample HC-11 is slightly less magnesian as XMg = 0.600.61, whereas that associated with symplectite in the same sample shows lower XMg of 0.58-0.59. A similar XMg range was also obtained for orthopyroxene in metadiorite (sample KDG1; XMg = 0.57-0.58). Orthopyroxene in metagabbro (sample HC-8) displays the lowest XMg of 0.39-0.42. Al content in the
Page 32 of 111
P a g e | 33 orthopyroxene is generally low (Al <0.1 pfu).
ip t
Plagioclase
Plagioclase in the examined rocks shows significant compositional variations
cr
depending on lithologies (Supplementary Table 4, Fig. 6d). Plagioclase in
us
symplectite-bearing metagabbro and mafic granulite (e.g., samples HC-8 and HC-11, respectively) shows higher anorthite contents of An61-77. Particularly, the
an
coarse-grained matrix plagioclase in sample HC-11 (An68-77) displays the highest anorthite contents of An68-77 as compared to that associated with symplectite
M
(An64-65). Plagioclase in metadiorite (sample KDG-1) shows the highest albite
Ac ce p
Amphibole
te
the samples as Or80-86.
d
content of An24-27. K-feldspar is compositionally nearly homogeneous throughout
Amphibole in metagranodiorite (samples WC-2 and KDG-4), metadiorite (sample KDG-1), mafic granulite (samples HC-1 and HC-11), and garnet amphibolite (sample KDG-2) shows notable compositional variation depending on its occurrence (Supplementary Table 2, Fig. 6e). The mineral is characterized by high Ca and (Na+K) contents (1.79-1.92 pfu and 0.54-0.94 pfu, respectively). Calcic amphibole in metagranodiorite and metadiorite is Si-rich, and classified as edenite – pargasite (sample KDG-1), ferro-edenite (sample KDG-4), and ferro-
Page 33 of 111
P a g e | 34 edenite – ferro-pargasite (sample WC-2). Amphibole in garnet amphibolite is slightly Al-rich and classified as ferro-pargasite. Calcic amphiboles in garnet-free mafic granulite (samples HC-1 and HC-11) are Mg- and Al-rich as pargasite. Ti
ip t
content of amphiboles in metagranodiorite, metadiorite, and amphibolite are low (0.18-0.25 pfu), but that in mafic granulites is slightly enriched (0.27-0.41 pfu).
cr
Calcic amphibole in metagabbro (sample HC-7), which occurs as inclusions in
us
poikiloblastic garnet, is significantly Ti-rich (0.52-0.54 pfu) and classified as
an
kaersutite.
M
Biotite
d
The biotite in metagranodiorite (samples WC-2 and KDG-4) and metadiorite
te
(sample KDG-1) exhibits very high TiO2 content of 3.9-5.6 wt.% (Fig. 6f). XMg of biotite in samples WC-2 and KDG-4 (metagranodiorite) is 0.44-0.51, which is
Ac ce p
lower than that in the sample KDG-1 (metadiorite) (XMg = 0.56-0.61).
4.2 Geothermobarometry
The garnet-clinopyroxene-plagioclase-quartz geothermobarometers were applied to porphyroblastic garnet and clinopyroxene with plagioclase and quartz assemblage in samples HC-7 and HC-8 (metagabbro). The estimated temperature range for garnet-clinopyroxene pairs are 830-860°C for sample HC7 and 840-860°C for sample HC-8 (at 10 kbar) based on the method of Ellis and
Page 34 of 111
P a g e | 35 Green (1979), which is a widely applied thermometer for mafic granulites and relies on experimental calibration of Fe–Mg fractionation between garnet and clinopyroxene at 750°C to 1350°C and 24 to 30 kbar. Calculated temperature
ip t
ranges using the method of Ganguly et al. (1996), which is based on revised solution model of garnet, are slightly lower but nearly consistent, 820-850°C
cr
(sample HC-7) and 790-800°C (sample HC-8). The temperatures were calculated
us
at 10 kbar, a reference pressure inferred from the peak pressure condition of the samples estimated by garnet-clinopyroxene-plagioclase-quartz geobarometer of
M
nearly consistent among all the samples.
an
Moecher et al. (1988) as discussed below. The estimated temperature range is
d
Metamorphic pressure for the metagabbro was calculated using garnet-
te
clinopyroxene-plagioclase-quartz assemblage in the samples HC-7 and HC-8 based on experimental calibration of Perkins and Newton (1981). The estimated
Ac ce p
results are 7.1-7.6 kbar for sample HC-7 and 6.6-6.8 kbar for sample HC-8 at 850ºC. These results are, however, significantly low if compared to the petrographic observation of quartz in the rock because the stability of quartz in mafic granulite possibly indicates high-pressure condition of P > 9 kbar at 850°C (e.g., Green and Ringwood, 1967). We therefore adopted the method of Moecher et al. (1988), which is based on improved garnet-clinopyroxene-plagioclasequartz geobarometer using new thermodynamic and experimental data. Application of the method yielded a pressure range of 9.5-9.8 kbar (sample HC7) and 9.4-9.6 kbar (sample HC-8) at 850°C, suggesting relatively high-pressure
Page 35 of 111
P a g e | 36 conditions.
The Fe-Mg cation exchange between garnet and hornblende has been calibrated
ip t
as a geothermometer by Graham and Powell (1984). We applied this method to garnet amphibolite (sample KDG-2) and obtained a temperature range of 700-
cr
720°C. The estimated temperatures are significantly lower than the results
us
obtained from garnet-clinopyroxene assemblage, and probably record retrograde
an
conditions.
As hornblende coexists with plagioclase in metagranodiorite (samples WC-2 and
M
KDG-4) and metadiorite (sample KDG-1), the mineral pair has been used for
d
geothermometry. Based on hornblende solid-solution models and well-
te
constrained natural and experimental studies, two geothermometers were calibrated by Holland and Blundy (1994) based on the edenite-tremolite reaction,
Ac ce p
which is applicable to quartz-bearing rocks, and edenite-richterite reaction, which is applicable to both quartz-bearing and quartz-free rocks. As quartz is present in the samples, the two thermometers were used for the calculation of temperature of the felsic rocks. The calculated results are 730−770°C for sample KDG-1, 700−740°C for sample KDG-4, and 720−780°C for sample WC-2, probably indicate retrograde re-equilibration.
4.3 Geochemistry
Page 36 of 111
P a g e | 37 Whole rock geochemical data, including major, minor, trace and rare earth elements, on nine representative samples (two metagranodiorites, one metadiorite, one garnet amphibolite, two mafic granulites, two metagabbros and
ip t
one garnet charnockite) are given in Supplementary Table 5. The salient
cr
geochemical features of the various rock types are summarized below.
us
4.3.1 Major and trace elements
an
Major element contents show that the metagranodiorite has high SiO2 content in a restricted range (71.0-71.3 wt.%) (Fig. 7), whereas the metadiorite displays
M
intermediate SiO2 of 55.8 wt.%. The lowest SiO2 content is displayed by the
d
garnet amphibolite (41.5%). The two mafic granulites have SiO2 content in the
te
range of 49.86-49.94 wt.%, whereas the metagabbro is slightly SiO2 poor as 44.7-48.8 % (Fig. 7). Interestingly, the metagranodiorites in the two different
Ac ce p
crustal units (KDG-4 and WC-2) show almost identical composition. P2O5 contents of the garnet amphibolite shows relatively high values, which is consistent with the abundant occurrence of apatite in the sample (~5%). Garnet amphibolite (KDG-2), metagabbro (HC-7) and metagranodiorites clearly plot away from the general trends displayed by other meta-mafic rocks except for their Na2O-SiO2 relationship. These meta-mafic rocks have systematic positive correlation of Al2O3, K2O and P2O5 and a negative correlation of MgO, Fe2O3, MnO, CaO, TiO2 with increasing silica content. This geochemical variability is clearly reflected by the SiO2 versus Fe2O3T/MgO diagram (Fig. 7j) suggesting
Page 37 of 111
P a g e | 38 three distinct geochemical groups of rocks: magnesian (Fe2O3T/MgO <3), ferroan (Fe2O3T/MgO >3) and felsic (SiO2 >70 %). The three compositional types can also be distinguished in terms of their SiO2 (magnesian: 48.8-54.5 wt.%; ferroan:
ip t
41.5-44.7 wt.%; felsic: >70 wt. %) and TiO2 (magnesian: 0.55-1.2 wt.%; ferroan:
cr
3.0-3.4 wt.%; felsic: 0.31-0.41 wt.%) contents.
us
The major element behaviour shown by the three groups of samples is also reflected in their trace elements (Fig. 8). In the primitive-mantle normalized (Sun
an
and McDonough, 1989) trace element plot, magnesian and felsic samples are characterized by enrichment of large ion lithophile elements (LILE) like Rb, Ba
M
and K, negative Nb, Ta, P and Ti anomalies, and relatively constant high-field
d
strength elements (HFSE: Zr, Hf) (Figs. 8a,c). Nevertheless, the ferroan
te
metagabbro (sample HC-7) is depleted in LILE and does not show Nb and Ta negative anomalies although Zr, Hf, Y and Yb are depleted (Fig. 8b). The garnet
Ac ce p
amphibolite (sample KDG-2) shows a characteristic depletion of LILE whereas the other trace elements are relatively enriched and yield a smooth pattern (Fig. 8b). Another striking feature is marked Sr enrichment in the ferroan samples and no systematic behaviour in the magnesian samples. The felsic samples however, are enriched in LILE with obvious negative anomalies of Ta and Nb, P, Ti and Sr. The felsic samples show approximately one magnitude higher contents of LILE as compared to the other two groups. Rb-Y-Nb and Rb-Yb-Ta diagrams of Pearce et al. (1984) for granitic rocks (Fig. 9) suggest that metagranodiorite and metadiorite are characterized by low Y, Yb, Nb, and Ta, and moderate Rb
Page 38 of 111
P a g e | 39 contents.
In the chondrite normalized (McDonough and Sun, 1995) REE plot (Figs. 8d-f),
ip t
the magnesian samples are characterized by light REE (LREE, e.g. La) enriched patterns with smooth constant heavy REE (HREE, e.g. Yb), except for the HREE
cr
depleted KDG-1 sample. There is no obvious Eu anomaly in the magnesian
us
samples although the garnet charnockite shows a slight positive Eu anomaly (Fig. 8d). The ferroan sample of the KC (KDG-2) has LREE enriched and HREE
an
depleted pattern whereas that of the HC (HC-7) has a LREE depleted pattern with smooth HREE. The garnet amphibolite (sample KDG-2) suggests
M
continuous fractionation from LREE to HREE ((La/Sm)cn = 2.0 and (Gd/Yb)cn =
d
2.4) (Fig. 8e). The felsic samples are highly enriched in LREE and show clear
te
negative Eu anomaly. Further, these samples have relatively flat HREE pattern. Although the samples of the felsic group are from two complexes (WC-1 and
Ac ce p
KDG-4), whole rock chemistry, trace element and REE patterns are remarkably similar.
4.3.2. Petrogenetic implications
The metagranodiorite and metadiorite fall in the field of volcanic-arc granites (VAG) in the Rb-Y-Nb and Rb-Yb-Ta diagrams of Pearce et al. (1984) (Fig. 9). The arc magmatic signature is also supported by the primitive-mantle normalized trace element plot (Figs. 8a,c) characterized by the enrichment of LILE, negative
Page 39 of 111
P a g e | 40 Nb, Ta, P and Ti anomalies, and relatively constant HFSE of the metagranodiorite and metadiorite. LREE enrichment relative to middle REE and relatively constant HREE of the metagranodiorite and metadiorite in the chondrite normalized REE
ip t
plot (Figs. 8d,f) also support volcanic arc affinity of the felsic to intermediate
cr
rocks.
us
In the TiO2-MnO-P2O5 triangular diagram (Fig. 10a) for basaltic rocks, the mafic granulite, magnesian metagabbro, and garnet charnockite (samples HC-1, HC-8,
an
HC-10, and HC-11) show island-arc tholeiite affinity, whereas the ferroan metagabbro (sample HC-7) and garnet amphibolite (sample KDG-2) are
M
classified as MORB and oceanic island alkali basalt, respectively. The arc
d
signature of the magnesian rocks is also supported by enrichment of LILE
te
elements in Fig. 8a, and LREE enrichment relative to middle REE and relatively constant HREE in Fig. 8d. On the other hand, continuous fractionation from
Ac ce p
LREE to HREE, high P2O5 content in TiO2-MnO-P2O5 plot, and high Nb content in Nb-Zr-Y plot support oceanic-island affinity of the garnet amphibolite. The MORB signature of ferroan metagabbro is supported by depleted in LILE and positive Nb and Ta anomalies in Fig. 8b, and slightly LREE-depleted and nearly flat HREE pattern in Fig. 8e.
4.4. Zircon U-Pb geochronology and Lu-Hf isotopes
4.4.1. U-Pb geochronology
Page 40 of 111
P a g e | 41 The zircon U-Pb analytical results are given in Supplementary Table 6, and zircon REE data are given in Supplementary Table 7. Representative CL images of the zircons from different rock types in this study are shown in Figs. 11-13, and
ip t
the age data are plotted in Figs. 14-17. Chondrite-normalized REE patterns are
cr
shown in Figs. 18 and 19.
us
Mafic granulite
Zircon grains in the mafic granulite sample HC-1 are translucent and colorless to
an
pale pink. Most of the grains are subhedral, displaying spherical to weakly elongate features with lengths ranging from 30 to 100 μm and an aspect ratio of
M
1.5:1 to 1:1. In CL images, the zircon grains are dark gray and structureless. A
d
total of 21 spots from 21 zircons were analyzed from this rock. Their Th, U and
te
Th/U values are in the range of 474-933 ppm, 1559-2280 ppm and 0.24-0.37, respectively (Supplementary Table 7). They exhibit fractionated REE patterns
Ac ce p
with LREE depletion and HREE enrichment, and weakly negative Eu anomalies (Fig. 18a). All of the analyzed spots form a coherent group on the concordia, yielding a weighted mean
206
Pb/238U age of 553.0 ± 3.2 Ma (MSWD = 2.4) (Fig
14a,b), which we interpret as the timing of metamorphism.
Metagabbro The zircon grains in metagabbro sample HC-7 are translucent, colorless or light brown. The euhedral to anhedral grains show ellipsoidal or prismatic morphology with a length of 50-200 μm and a length to width ratio of 2:1 to 1:1. In
Page 41 of 111
P a g e | 42 CL images, the grains show slight patchy zoning patterns or structureless domains with homogeneous gray color. A total of 15 spots from 15 zircons were analyzed from this rock. Their Th, U and Th/U contents show ranges of 5-95
ip t
ppm, 32-298 ppm, and 0.14-0.33, respectively (Supplementary Table 7). They exhibit slight fractionated REE patterns with LREE depletion and HREE
cr
enrichment, and variable negative Eu anomalies (Fig. 18b). All the spots are
us
distributed near the concordia line and can be divided into two groups (Fig 14c,d). Seven spot ages range from 555 ± 16 Ma to 589 ± 12 Ma and define a
an
weighted mean 206Pb/238U age of 576.8 ± 9.3 Ma (MSWD = 0.92) with Th/U ratios (0.15-0.33). Another eight spots form a coherent tight group on the concordia and 206
Pb/238U age of 523.0 ± 7.1 Ma (MSWD = 0.54) with
M
yield a weighted mean
Ac ce p
Garnet charnockite
te
distinct thermal events.
d
Th/U ratios (0.14-0.32). These two groups of ages are regarded to represent two
The zircon grains in garnet charnockite sample HC-10 are translucent,
colorless or light brown. The subhedral to anhedral grains show ellipsoidal or irregular morphology with a length of 50-300 μm and a length to width ratio of 3:1 to 1:1. In CL images, some grains possess dark inherited core and bright rim texture. A total of 29 spots from 29 zircons were analyzed from this rock. Their Th, U and Th/U contents show a range of 5-289 ppm, 36-758 ppm, and 0.021.42, respectively (Supplementary Table 7). The inherited cores of zircon exhibit distinct fractionated REE patterns and with negative Eu anomalies (Fig. 18b),
Page 42 of 111
P a g e | 43 whereas the bright rims of zircon shows slightly fractionated and even flat REE patterns with variable negative Eu anomalies (Fig. 18c). The majority of U-Pb data are distributed along the concordia or along the Pb loss line and yield a
ip t
lower intercept age of 533 ± 33 Ma (MSWD = 4.6) which is taken to indicate the timing of metamorphism of this rock. Twenty spots show ages ranging from 509.0
206
Pb/238U ages of 511.1 ±
us
be divided into three groups yielding weighted mean
cr
± 5.8 Ma to 624.0 ± 11.4 Ma and form a cluster near the lower intercept and can
5.9 Ma (MSWD = 0.13), 540.4 ± 6.0 Ma (MSWD = 1.5) and 579 ± 10 Ma (MSWD
an
= 0.13) (Fig 15a,b) with Th/U ratios 0.21-0.52, 0.24-1.42, and 0.02-0.24, respectively. These ages might represent three distinct thermal events of zircon
d te
Mafic granulite sill
M
growth.
The zircons from mafic-granulite sill sample HC-11 are colorless or light
Ac ce p
brown, transparent to translucent. The subhedral to anhedral grains show ellipsoidal, irregular or elongate morphology with a length of 40-300 μm and a length to width ratio of 3:1 to 1:1. In CL images, most zircons are bright with structureless or slight patchy zoning, and a few zircons show dark core and bright rim textures, suggesting their complex origin. A total of 22 spots from 22 zircons were analyzed from this rock. Their Th, U and Th/U contents show a range of 11-205 ppm, 47-720 ppm, and 0.09-0.91, respectively (Supplementary Table 7). Most spots exhibit slightly fractionated or flat REE patterns with variable negative Eu anomalies (Fig. 18d). Sixteen analyses of CL-bright, structureless
Page 43 of 111
P a g e | 44 domains exhibit coherent age and yield a weighted mean
206
Pb/238U age of 539.1
± 4.4 Ma (MSWD = 1.6) (Fig 15 c, d). Their Th/U ratios range between 0.14 and 0.91 with 10 of the 16 spot showing Th/U>0.3. The remaining six analyses of CL-
ip t
dark core or bright patchy domains show older ages between 563 ± 4 Ma and 665 ± 10 Ma. These spots have variable Th/U ratios from 0.09 to 0.35. The
cr
oldest core age of ca. 665 Ma might represent the timing of magmatic
us
emplacement of these rocks, and the weighted mean age of ca. 539 Ma from the
an
bright structureless zircons marks the metamorphic age.
M
Metadiorite
d
The zircons from metadiorite sample KDG-1 are colorless or light brown,
te
transparent to translucent. The euhedral to subhedral grains show near-spherical or prismatic morphology with a length of 40-200 μm and a length to width ratio in
Ac ce p
the range of 3:1 to 1:1. In CL images, most zircons show oscillatory-zoned core and bright thin rim texture. A total of 34 spots from 33 zircons were analyzed from this rock. Their Th, U and Th/U contents show a range of 119-806 ppm, 226-1417 ppm, and 0.13-1.06, respectively (Supplementary Table 7). They exhibit distinct fractionated patterns with variable negative Eu anomalies (Fig. 19a). All of the data when regressed define an upper intercept age of 916 ± 57 Ma and a lower intercept age of 462 ± 180 Ma (MSWD = 2.7) (Fig. 16 a,b). Among these, thirty spots exhibit a wide range of 206Pb/238U ages between 659 ± 5 Ma and 1086 ± 11 Ma suggesting variable Pb loss after the emplacement of the protolith in early
Page 44 of 111
P a g e | 45 Neoproterozoic. The remaining four spots form a coherent group near the lower intercept and yield a weighted mean
206
Pb/238U age of 532 ± 18 Ma (MSWD =
ip t
6.2), representing the metamorphic age.
cr
Garnet amphibolite
us
Zircons in the garnet amphibolite sample KDG-2 are colorless or light brown, transparent to translucent. The subhedral to anhedral grains show near-spherical
an
or irregular morphology with a length of 50-250 μm and a length width ratio of 3:1 to 1:1. Most grains are structureless and light gray in CL images. A total of 21
M
spots from 21 zircons were analyzed from this rock. Their Th, U and Th/U
d
contents show a range of 1-96 ppm, 37-195 ppm, and 0.03-0.60, respectively
te
(Supplementary Table 7). Their REE patterns can be divided into two types, the first type is characterized by LREE depletion and HREE enrichment, whereas the
Ac ce p
second type by flat HREE (Fig. 19b). Two spot ages show 923 Ma (Th/U=0.03) and 848 Ma (Th/U=0.07). Nineteen spots with
206
Pb/238U age ranging from 501 ±
9 Ma to 549 ± 14 Ma form a coherent tight group and yield a weighted mean 206
Pb/238U age of 520.7 ± 6.6 Ma (MSWD = 1.6) with Th/U ratios 0.32-0.60
(Supplementary Table 7, Fig. 16 c,d ). The Cambrian ages might mark the timing of metamorphism of this rock.
Meta-granodiorite
Page 45 of 111
P a g e | 46 Zircons in meta-granodiorite sample KDG-4 are colorless or light brown, and transparent to translucent. The euhedral to subhedral grains show near-spherical or prismatic morphology with a length of 30-200 μm and a length to width ration
ip t
of 3:1 to 1:1. In CL images, most zircons show clear oscillatory zoning and some of the grains display heterogeneous fractured domains. A total of 33 spots from
of
182-2242
ppm,
385-2485
ppm,
and
0.35-1.51,
respectively
us
range
cr
31 zircons were analyzed from this rock. Their Th, U and Th/U contents show a
(Supplementary Table 7), suggesting magmatic crystallization. Moreover, these
anomalies (Fig. 19c). Spot 4.4 yields
206
an
spots exhibit distinct fractionated REE patterns with strong negative Eu Pb/238U spot age of 538 ± 6 Ma, probably
M
corresponding to metamorphism. The remaining 32 spots are distributed along 207
Pb/206Pb age of 980 ± 16 Ma
d
the concordia and yield a weighted mean
this rock.
te
(MSWD = 1.13) (Fig. 16 e,f), which is taken to indicate the emplacement age of
Ac ce p
Zircons in the meta-granodiorite sample WC-2 are colorless or light brown
and translucent. The euhedral to subhedral grains show near-spherical or prismatic morphology with a length of 50-250 μm and a length to width ratio of 3:1 to 1:1. In CL images, most zircons show core-rim texture with weakly light rim and clear oscillatory zoning or show heterogeneous fractured core. A total of 37 spots from 37 zircons were analyzed from this rock. Their Th, U and Th/U contents show a range of 56-776 ppm, 55-1550 ppm, and 0.08-1.05, respectively (Supplementary Table 7). They have variable LREE contents but similar and fractionated HREE patterns with negative Eu anomalies (Fig. 19d). Six spots
Page 46 of 111
P a g e | 47 form a tight group (Group 1) with and a weighted mean
206
206
Pb/238U ages between 789 Ma and 822 Ma
Pb/238U age of 805 ± 12 Ma (MSWD = 4.0) with high
Th/U ratios 0.34-0.96. Six spots define the second group (Group 2) with a narrow
206
206
Pb/238U ages between 723 Ma and 745 Ma and a weighted mean
ip t
range of
Pb/238U age of 734.0 ± 4.6 Ma (MSWD = 1.1) (Th/U ratios 0.19-1.05) (Fig. 17). 206
Pb/238U ages
cr
Eighteen spots constitute another tight group (Group 3) with
us
between 524 Ma and 565 Ma and a weighted mean 206Pb/238U age of 546.0 ± 5.7 Ma (MSWD = 7.8) and relatively low Th/U ratios 0.08-0.48 (15 of the 18 spots
an
show Th/U<0.27, 12 of the 18 spots show Th/U<0.18). The result suggests magmatic emplacement of the protolith at ca. 805 Ma and another thermal event
Ac ce p
te
4.4.2 Zircon Lu-Hf isotopes
d
M
at ca. 734 Ma, followed by metamorphism ca. 546 Ma.
In situ Hf isotope compositions of the dated zircons are given in Table 8, and
plotted in Fig 20. In the eight samples analyzed, the fLu/Hf values of zircons vary from -0.92 to -1.0, which are obviously lower than the fLu/Hf values of mafic crust (0.34, Amelin et al., 2000) and sialic crust (-0.72, Vervoort and Patchett, 1996). Therefore, the two-stage model age is considered to be a better representation than the single stage model age to evaluate the time of source material extraction from the depleted mantle or the residence time of the source material in the crust (Blichert-Toft and Albarede, 1997).
Page 47 of 111
P a g e | 48 Mafic granulite
Seventeen zircons were analyzed from mafic granulite sample HC-1 for Lu176
Hf/177Hf ratios in the range of 0.282365 to
ip t
Hf isotopes. The result show initial
0.282432. The εHf(t) values show a tight cluster from -2.2 to 0.1 with an average
cr
value of -1.1, falling close the line of the Chondrite Uniform Reservoir (CHUR)
us
(Fig. 20). The Hf depleted model ages (TDM) are between 1135 Ma and 1227 Ma (mean 1184 Ma), and Hf crustal model ages (TDMC) show a range of 1501-1651
an
Ma (mean 1582 Ma). The results suggest a mixed source from both juvenile
M
Neoproterozoic and reworked Mesoproterozoic components.
te
d
Metagabbro
Twelve zircons were analyzed from metagabbro sample HC-7 for Lu-Hf
Ac ce p
isotopes. The result shows initial
176
Hf/177Hf ratios ranging from 0.282221 to
0.282280. The εHf(t) values are in the range of -7.5 to -5.5 with an average value of -6.3. The Hf depleted model ages (TDM) are between 1338 Ma and 1418 Ma (mean 1372 Ma), and Hf crustal model ages (TDMC) in the range of 1847-1978 Ma (mean 1902 Ma). The predominantly negative εHf(t) values suggest that the source material involved reworked Paleoproterozoic crust.
Garnet charnockite
Page 48 of 111
P a g e | 49 Twenty-five zircons were analyzed from the garnet charnockite sample HC10 for Lu-Hf isotopes. The result show initial
176
Hf/177Hf ratios ranging from
0.281500 to 0.282144. The εHf(t) values are in the range of -33.3 to -10.5 with an
ip t
average value of -17.7. The Hf depleted model ages (TDM) are between 1522 Ma and 2401 Ma (mean 1809 Ma), and Hf crustal model ages (TDMC) varies 2156-
cr
3580 Ma (mean 2614 Ma averagely). The highly negative εHf(t) values suggest
us
that the source material for the charnockite was reworked Paleoproterozoic crust. Twenty zircons were analyzed from mafic granulite sample HC-11 for Lu-Hf 176
Hf/177Hf ratios in the range of 0.281857 to
an
isotopes. The result show initial
0.282095. The εHf(t) values are in the range of -20.5 to -12.1 with an average
M
value of -14.1. The Hf depleted model ages (TDM) are between 1590 Ma and
d
1922 Ma (mean 1666 Ma), and Hf crustal model ages (TDMC) varies from 2263 to
te
2790 Ma (mean 2386 Ma). The dominantly negative εHf(t) values suggest that the source material for the magma evolved from the Neoarchean-Paleoproterozoic
Ac ce p
crust.
Metadiorite
Twenty-six zircons were analyzed from metadiorite sample KDG-1 for Lu-Hf
isotopes. The result show initial
176
Hf/177Hf ratios ranging from 0.282270 to
0.282505. The εHf(t) values are in the range of -1.1 to 7.2 with an average value of 1.5. The Hf depleted model ages (TDM) are between 1031 Ma and 1360 Ma (mean 1259 Ma), and Hf crustal model ages (TDMC) varies from 1206 to 1733 Ma
Page 49 of 111
P a g e | 50 (mean 1569). The results suggest a mixed source from both juvenile and Proterozoic components.
ip t
Garnet amphibolite
176
Hf/177Hf ratios in the range of
us
for Lu-Hf isotopes. The result show initial
cr
Twenty zircons were analyzed from the garnet amphibolite sample KDG-2
0.281845 to 0.282058. The εHf(t) values are in the range of -21.3 to -13.8 with an
an
average value of -16.5. The Hf depleted model ages (TDM) are between 1641 Ma and 1931 Ma (mean 1744 Ma), and Hf crustal model ages (TDMC) varies from
M
2356 to 2828 Ma (mean 2526 Ma). The predominantly negative εHf(t) values
Ac ce p
Meta-granodiorite
te
Paleoproterozoic crust.
d
suggest that the source material was evolved from the Neoarchean-
Twenty-six zircons were analyzed from meta-granodiorite sample KDG-4 for
Lu-Hf isotopes. The result show initial
176
Hf/177Hf ratios ranging from 0.282230 to
0.282341. The εHf(t) values are in the range of 2.5 to 6.4 with an average value of 3.9. The Hf depleted model ages (TDM) are between 1266 Ma and 1417 Ma (mean 1365 Ma), and Hf crustal model ages (TDMC) varies from 1431 to 1679 Ma (mean 1592 Ma). The results suggest a mixed source from both juvenile and Paleo-Mesoproterozoic components. The dominantly positive εHf(t) values
Page 50 of 111
P a g e | 51 suggest that the source material was evolved from the Mesoproterozoic juvenile material. Thirty-two zircons were analyzed from metagranodiorite sample WC-2 for 176
Hf/177Hf ratios range from 0.281968 to
ip t
Lu-Hf isotopes. The result show initial
0.282283. The εHf(t) values are in the range of -15.0 to -3.9 with an average
cr
value of -6.2. The Hf depleted model ages (TDM) are between 1344 Ma and 1789
us
Ma (mean 1440 Ma), and Hf crustal model ages (TDMC) varies from 1799 to 2498 Ma (mean 1944 Ma). The dominantly negative εHf(t) values suggest that the
an
source material was evolved from reworked Paleoproterozoic crust.
d
M
5. Discussion
te
The samples analyzed in this study are petrographically classified as metagranodiorite (samples KDG-4 and WC-2; Pl + Kfs + Qtz + Bt + Hbl),
Ac ce p
metadiorite (sample KDG-1; Pl + Qtz + Bt + Hbl), mafic granulite (samples HC-1 and HC-11; Cpx + Opx + Pl + Hbl), metagabbro (samples HC-7 and HC-8; Grt + Cpx + Opx + Pl), garnet amphibolite (sample KDG-2; Grt + Hbl + Pl) and garnet charnockite (sample HC-10; Grt + Pl + Qtz + Opx + Cpx). Pressure and temperature conditions of the study area were estimated based on conventional geothermobarometry on mafic granulite, amphibolite, metagranodiorite, and metadiorite. The peak P-T conditions were estimated from garnet-clinopyroxeneplagioclase-quartz assemblages in samples HC-7 and HC-8 as 830-860°C and 9.4-9.8 kbar. Slightly lower temperature ranges of 700-720°C obtained for garnet-
Page 51 of 111
P a g e | 52 hornblende pairs in garnet amphibolite (sample KDG-2) and 720−780°C for hornblende-plagioclase pairs in metagranodiorite and metadiorite probably correspond to retrograde conditions. The peak metamorphic conditions of ~860
ip t
°C and ~9.8 kbar estimated in this study are consistent with previous P-T estimates from the Highland Complex (~850-900 °C and ~8-10 kbar; Schenk et
cr
al., 1988; Schumacher et al., 1990; Kriegsman, 1991; Hiroi et al., 1994, Raase
us
and Schenk, 1994; Schumacher and Faulhaber, 1994; Kriegsman and Schumacher, 1999; Dharmapriya et al., 2014). The UHT metamorphic conditions
an
of >900 °C or even >1100 °C estimated based on the stability of diagnostic mineral assemblages (e.g., sapphirine + quartz) and phase equilibria (e.g.,
M
Kriegsman and Schumacher, 1999; Osanai et al., 2006; Sajeev and Osanai,
te
Highland Complex.
d
2004a, 2004b; Sajeev et al., 2007) is mainly restricted to some localities in the
Ac ce p
The petrographic and geochemical features discussed in previous sections suggest distinct petrogenetic history for the protoliths of the meta-igneous suite from the crustal blocks of Sri Lanka. The REE and trace element patterns (Fig. 8) as well as Rb-Y-Nb and Rb-Yb-Ta diagrams (Fig. 9) reveal volcanic arc affinity for the metagranodiorite, metadiorite and garnet charnockite, suggesting that the protoliths of the rocks were derived from felsic to intermediate arc magmas. The two mafic granulites (HC-1 and HC-11) and magnesian metagabbro (HC-8) also suggest volcanic arc affinity. The garnet-free mafic granulite (sample HC-1) occurring within marble and magnesian metagabbro (sample HC-8) enclosed in
Page 52 of 111
P a g e | 53 quartzite shows evidence for subduction-related mafic magmatism in the crustal blocks of WC, KC and the HC.
ip t
In contrast, the protolith of sample HC-7 (ferroan metagabbro from the Highland Complex) is inferred as N-MORB, suggesting the accretion of remnant accreted
cr
oceanic lithosphere accreted along with the arc components and trench
us
sediments during the subduction-collision event. The garnet amphibolite (KDG-2)
an
displays a signature of oceanic island alkali basalt.
From the Kadugannawa Complex, a metadiorite (KDG-1), garnet amphibolite
M
(KDG-2), and a meta-granodiorite (KDG-4) were analyzed in this study. Zircons in
d
the metadiorite define an upper intercept age of 916 ± 57 Ma and a lower
te
intercept age of 462 ± 180 Ma (MSWD = 2.7). Metamorphic zircons in this rock form a coherent group and yield a weighted mean
206
Pb/238U age of 532 ± 18 Ma
Ac ce p
(MSWD = 6.2). Two zircon grains from the garnet amphibolite show spot ages of 923 Ma and 848 Ma (Th/U=0.07). The majority of zircons in this rock are of metamorphic origin and define a weighted mean
206
Pb/238U age of 520.7 ± 6.6
Ma (MSWD = 1.6). Zircons in the meta-granodiorite sample yield a weighted mean
207
Pb/206Pb age of 980 ± 16 Ma (MSWD = 1.13) marking the timing of
emplacement of the protolith. From the Wanni Complex, zircons in a metagranodiorite (WC-2) define three groups with a weighted mean
206
Pb/238U
age of 805 ± 12 Ma (MSWD = 4.0) (emplacement of the magmatic protolith) 734.0 ± 4.6 Ma (MSWD = 1.1) (Cryogenian thermal event) and 546.0 ± 5.7 Ma
Page 53 of 111
P a g e | 54 (MSWD = 7.8) (latest Neoproterozoic-Cambrian metamorphism). Three samples were analyzed from the Highland Complex. All of the analyzed spots in a mafic granulite (HC-1) form a coherent group on the 206
Pb/238U age of 553.0 ± 3.2 Ma (MSWD =
ip t
concordia, yielding a weighted mean
2.4) marking the timing of metamorphism. Zircons in a metagabbro sample (HC206
Pb/238U ages of 576.8 ± 9.3 Ma
cr
7) form two groups with weighted mean
us
(MSWD = 0.92) and 523.0 ± 7.1 Ma (MSWD = 0.54) suggesting two distinct thermal events. Zircons from a garnet charnockite (HC-10) show Pb loss with a
an
lower intercept age of 533 ± 33 Ma (MSWD = 4.6) marking the timing of metamorphism. Another group of zircons from this rock show spot ages in the
M
range of 509.0 ± 5.8 Ma to 624.0 ± 11.4 Ma with weighted mean
206
Pb/238U ages
d
of 511.1 ± 5.9 Ma (MSWD = 0.13), 540.4 ± 6.0 Ma (MSWD = 1.5) and 579 ± 10
te
Ma (MSWD = 0.13), suggesting different stages of zircon growth during the late Neoproterozoic. Zircons in a mafic granulite sill (HC-11) yield a weighted mean Pb/238U age of 539.1 ± 4.4 Ma (MSWD = 1.6). Another group with CL-dark
Ac ce p
206
core or bright patchy domains shows older ages between 563 ± 4 Ma and 665 ± 10 Ma, and a weighted mean metamorphic age of 539 Ma from the bright structureless zircons.
Zircons in the metadiorite (KDG-1) and metagranodiorite (KDG-4) from the
Kadugannawa complex show mostly positive εHf(t) values in the range of -1.1 to 7.2 and Hf crustal model ages (TDMC) varies from 1206 to 1733 Ma suggesting a mixed source from both juvenile and Paleo-Mesoproterozoic components. However, zircons in the garnet amphibolite (KDG-2) show dominantly negative
Page 54 of 111
P a g e | 55 εHf(t) values (mean -16.5) with TDMC in the range of 2356 to 2828 Ma suggesting reworked Neoarchean-Paleoproterozoic crustal source. Zircons in the Wanni Complex metagranodiorite (WC-2) also possess negative εHf(t) values (mean -
ip t
6.2) with TDMC in the range of from 1799 to 2498 Ma suggesting reworked Paleoproterozoic crust as the magma source.
cr
The Highland Complex rocks preserve distinct imprints of reworking of
us
older crust. The zircon εHf(t) values in mafic granulite show a tight cluster from 2.2 to 0.1 with TDMC in the range of 1501-1651 Ma suggesting a mixed source
an
from both juvenile Neoproterozoic and reworked Mesoproterozoic components. Zircons in the metagabbro (HC-7) from this complex show negative εHf(t) values
crust as the source.
M
(mean -6.3) and TDMC of 1847-1978 Ma suggesting reworked Paleoproterozoic Zircons in the garnet charnockite (HC-10) also display
te
d
highly negative εHf(t) values (mean -17.7) and older TDMC) (mean 2614 Ma) suggesting reworked Paleoproterozoic crust. Zircons in the mafic granulite sill
Ac ce p
(HC-11) also show negative εHf(t) values (mean -14.1) and TDMC between 2263 to 2790 Ma suggesting that the source material for the magma evolved from the Neoarchean-Paleoproterozoic crust. The εNd(t) and εHf(t) isotopic data reported in previous studies and those
from the present study from the Sri Lankan crustal blocks provide evidence for magma generation from mixed sources involving both juvenile and reworked components. The Nd isotope data reported elsewhere for the KC and the WC show highly heterogeneous εNd(t) values, ranging from +2 to -10 (Willbold et al., 2004), and +4 to -8 (Millisenda et al., 1994). These data suggest variable crustal
Page 55 of 111
P a g e | 56 contamination of mantle-derived melts in a continental magmatic arc setting. However, in clear contrast, the HC rocks show dominantly negative εNd(t) values from -7 to -25 suggesting a different tectonic genetic setting from that of the KC
ip t
and the WC, with substantial contribution from older reworked crustal components.
cr
In summary the 916 to 980 Ma ages from the Kadugannawa Complex
us
represent arc magmatism during early Neoproterozoic. Zircon in the diorites from this complex show positive εHf(t) values suggesting juvenile addition. This was
an
followed by the 805 Ma tonalitic magma emplacement in the Wanni Complex. Multiple thermal events during Neoproterozoic are recorded from the Wanni and
M
Highland Complexes, with subduction-related signature, culminating in Cambrian
te
d
high-grade metamorphism (~550 Ma) that overprinted all the three complexes.
We propose a tectonic model where we envisage that the Wanni and Vijayan
Ac ce p
complexes as continental arcs underlain by Neoarchean – Paleoproterozoic basement (Fig. 21). The similarity in the arc-related geochemical features of the magmatic rocks emplaced broadly and coevally in both Wanni and Vijayan complexes prompt us to suggest a double-sided subduction of an oceanic lithosphere during early to mid Neoproterozoic. The broadly negative εNd(t) and εHf(t) values of these magmatic suites, together with the Neoarchean – Paleoproterozoic two stage Hf crustal model ages are also in accordance with melting of underlying older basement in both of these blocks, as well as the presence of older crust within the WC. The mafic/intermediate suites in the KC
Page 56 of 111
P a g e | 57 complex
(metadiorites)
show
input
of
juvenile
components
within
a
suprasubduction setting. In our model, the WC is envisaged as an accretionary suture with both oceanic and continental components, together with fragments of
ip t
older continental crust, analogous to the tectonic milieu in modern Phanerozoic accretionary belts such as the Central Asian Orogenic Belt (e.g., Xiao and The N-MORB signature of the garnet-bearing metagabbro
cr
Santosh, 2014).
us
(sample HC-7) and the oceanic alkali basalt chemistry of the garnet amphibolite (sample KDG-2) support our model of the accretion of remnants of the oceanic
an
lithosphere along with trench sediments onto the arc (continent). Occurrence of granitic-ring structures, serpentinites, anorthosites, and Cu-magnetite deposits in
M
association with chert bands and apatite in the VC at the boundary with the HC
d
(Pathirana, 1980; Cooray, 1984) indicates active continental margin volcanism
te
and supra-subduction zone magmatism at the junction of the HC and the VC. The submarine canyon with wall heights >1 km in the north east VC and hot
Ac ce p
springs are contiguous with the HC-VC contact zone. These features lend further support our subduction model.
Importantly, our study identifies juvenile magmatic additions during mid
Neoproterozoic (KDG-1, associated with subduction) and latest Neoproterozoic – Cambrian (HC-1; magmatism in post-collisional asthenospheric upwelling). The post-collisional slab-break off and asthenospheric upwelling is also envisaged as the principal contributor of heat that triggered widespread ultrahigh-temperature metamorphism within the HC, similar to the scenario envisaged for the formation
Page 57 of 111
P a g e | 58 of the UHT rocks elsewhere (e.g., Santosh et al., 2013).
As discussed in earlier sections, previous studies had identified the HC, WC and
ip t
VC as distinct tectonic units based largely on Nd model age mapping, combined with regional structural interpretations. The Nd model ages are as follows: HC=2-
cr
3.4 Ga; WC=1-2 Ga; VC=1-1.8Ga; and KC=1.4-1.8Ga. The zircon ages from
us
these three different tectonic units (either magma emplacement ages or depositional ages of sediments) are HC =~ 2.3Ga; WC=750-1.1Ga; VC=650Clearly the HC preserves records of the oldest
an
1.9Ga; and KC=0.8-1.0Ga.
sources of sedimentary units in Sri Lanka, and/or the oldest recycled continental
M
crust in this island. In the present study also, the oldest TDMC (2790 Ma) is Such older
d
obtained from zircons in a garnet charnockite within the HC.
te
components are common in many younger accretionary belts where ancient micro-continents or arcs are accreted and admixed during the final collision In the present case, the subduction-accretion continued until late
Ac ce p
stage.
Neoproterozoic and the final architecture of the HC was built during the Cambrian collisional event.
The ‘flower-style’ extrusion of the HC is clearly
marked by the thrust zone mapped in previous studies at the border of the HC with the VC; the structural break is also corroborated by geochronological studies that show a distinct boundary between the HC and the VC. However, the boundary is unclear on the western margin of the HC at its contact with the WC, where the large arc magma chamber of the KC is found emplaced at the inferred boundary. Away from the KC to the north and south-west also the HC-WC
Page 58 of 111
P a g e | 59 boundary is not clear. However, the doubly plunging synforms defined by the KC are considered to mark the structural boundary in the central region along the HC-WC inferred boundary. Such large-scale sheath folds are considered as the
ip t
hallmarks of major accretionary orogens elsewhere (e.g., Chetty et al., 2012). According our model, KC is not a discrete crustal block, but part of a disrupted
cr
huge arc magma chamber that was exhumed and transposed along the margin
us
of the WC.
an
Some of the existing models suggest that the WC, HC and VC were discrete terranes (Vitanage, 1972, 1985). The WC is proposed to have collided with HC
M
first, and subsequently, both WC and HC together as a unified block collided
d
with, and was thrust over, the VC (Vitanage, 1972; Voll and Kleinschrodt, 1991;
te
Kröner and Jaeckel, 1994; Kriegsman, 1995) during the amalgamation of Gondwana (Kriegsman, 1995; Kehelpannala 1997, 2003, 2004; Kröner et al.
Ac ce p
2003). However, there is no robust geochronologic and tectonic evidence to support this speculation, particularly the two discrete docking events, which should produce two separate metamorphic orogens with different ages along both margins of the HC. In fact, our data, coupled with available information clearly show that the WC and the VC were coeval Neoproterozoic arcs which were welded along the HC during late Neoproterozoic-Cambrian. The Cambrian collision event was accompanied by metamorphic overprint in all the rocks, formation of UHT metamorphic orogen along the HC, and post-collisional mantlederived magmatism within the HC (e.g. sample HC-1 in our study), all suggesting
Page 59 of 111
P a g e | 60 the HC as a collisional suture.
Whether the older continental basement
underlying WC and VC were contiguous prior to Neoproterozoic, or whether these bocks were derived independently from other continental blocks, remains
ip t
to be investigated in future studies.
cr
In previous studies, the magmatic episodes ranging from the Grenvillian
us
(ca. 1.1 Ga) to ca. 882 Ma in the KC (e.g. Kröner et al., 2003; Willbold et al., 2004) were considered to mark the beginning of the breakup of the Rodinia
an
supercontinent. However, the Neoproterozoic magmatic pulses in the Sri Lankan crustal complexes are clearly related to convergent margin setting as shown in
M
our study, as against the extensional setting required for Rodinia breakup. Except
d
for the Late Neoproterozoic-Cambrian mafic dykes (sample HC-1), and some
te
latest Neoproterozoic syn- and post-collisional granitoids in the HC, there are no other rock records of widespread extension-related mafic or bimodal magmatism
Ac ce p
in Sri Lanka during mid Neoproterozoic, and therefore a Rodinia break-up model to explain the magmatic suites is untenable. Accretion of oceanic components and arc magmatism in a convergent margin setting are also clearly indicated from our results. The global assembly of the Rodinia supercontinent with South China at its center is considered to have occurred at ca. 900 Ma (Li et al., 2008). Initial rifting of Rodinia started at ca. 860 with widespread breakup during 830800 Ma. From the Eastern Ghats Belt in India, Korhonen et al. (2011) reported SHRIMP ages of 980 Ma in monazites associated with garnet and cordierite from ultrahigh-temperature Mg–Al granulites, marking the timing of collisional
Page 60 of 111
P a g e | 61 amalgamation of this terrane within the Rodinia assembly. They also reported a range of
207
Pb/235U monazite ages (from ca. 1014 Ma to 959 Ma for one sample
and ca. 1043 Ma to 922 Ma for the second sample) suggesting protracted
ip t
monazite growth during the high-temperature retrograde evolution, which would broadly correspond with the global assembly of the Neoproterozoic Rodinia However, the data from crustal blocks in Sri Lanka indicate
cr
supercontinent.
us
protracted arc magmatism in a convergent margin setting, commencing in early Neoproterozoic, continuing through mid Neoproterozoic and culminating in
suture (HC).
an
regional metamorphism in all the blocks, with UHT imprint along the collisional The arc-related signature of these magmatic suites typical of
M
convergent margin magmas prompt us to suggest that the Neoproterozoic
d
magmatic events culminating in Cambrian collisional metamorphism might be
te
related to the prolonged Pacific-type orogeny associated with the assembly of Gondwana supercontinent, as also proposed for the Southern Granulite Terrane
Ac ce p
in India (Santosh et al., 2009).
Thus, our model deviates from the existing
concepts that link Sri Lanka with Rodinia rifting, and correlates the Neoproterozoic arc-related magmatic suites with convergent margin processes associated with the assembly of Gondwana.
8. Conclusions
Petrological and geochemical analyses and zircon U-Pb and Lu-Hf data from a suite of metamorphosed magmatic rocks from the major tectonic units in Sri
Page 61 of 111
P a g e | 62 Lanka, coupled with the available information from the crustal blocks, lead us to the following conclusions.
ip t
1. Regional metamorphism in the crustal bocks reached up to 860°C and 9.8 kbar, together with the extruded metamorphic orogen in the Highland
cr
Complex showing ultrahigh-temperature signature, which we correlate to
us
heat input from asthenospheric sources in a post-collisional slab-break off
an
setting.
2. Trace element and rare earth element patterns as well as Rb-Y-Nb and
M
Rb-Yb-Ta discrimination plots show volcanic arc affinity and protolith
te
d
crystallization from felsic to intermediate arc magmas.
3. The mafic magmatic units including those incorporated within the
Ac ce p
metasedimentary belt of Highland Complex show N-MORB signature, or oceanic island alkali basalt affinity, suggesting accretion of the remnants of oceanic lithosphere during the subduction-collisional event.
4. Zircon U-Pb ages ranging from 916 to 980 Ma from the Kadugannawa Complex represent arc magmatism during early Neoproterozoic, followed by the 805 Ma granodioritic magma emplacement in the Wanni Complex. Repeated thermal events during mid and late Neoproterozoic are also recorded from the Wanni and Highland Complexes, culminating in
Page 62 of 111
P a g e | 63 Cambrian high-grade metamorphism that reached ultrahigh-temperature conditions.
ip t
5. Zircon Lu-Hf data reveal mixed sources with the involvement of both juvenile and reworked Neoarchean – Paleoproterozoic basement. There
cr
is also clear indication of juvenile addition during both Neoproterozoic and
us
Cambrian tectonothermal events.
6. According to our data, the KC is part of a disrupted huge arc magma
an
chamber that was exhumed and transposed along the margin of the WC
d
simply a part of the WC.
M
which was previously regarded to be either an exotic crustal block or
te
7. A model of double-sided subduction during the Neoproterozoic, where the Wanni Complex to the west and the Vijayan Complex to the east represent
Ac ce p
continental arcs during the Neoproterozoic, culminating in collision along the Highland Complex during late Neoproterozoic-Cambrian is proposed. The occurrence of supra-subduction zone magmatic complexes and volcanic suites formed in a convergent margin setting along the contact of the VC with the HC lend support to our model.
8. The data presented in this study are consistent with a prolonged convergent margin setting from early to late Neoproterozoic, followed by final collisional amalgamation in the Cambrian, with the Highland Complex
Page 63 of 111
P a g e | 64 defining an accretionary belt as well as the collisional suture.
9. Both vertical (arc magmatic) and lateral (accretionary) growth of
ip t
continental crust during the Neoproterozoic is well recorded in the crustal
cr
blocks of Sri Lanka.
us
Acknowledgments
an
We thank Prof. Guochun Zhao, Editor and two referees for constructive and helpful comments. We extend our thanks to Qiongyan Yang and Xueming Teng at
M
the China University of Geosciences Beijing for help with zircon analyses. This study contributes to the Talent Award to M. Santosh from the Chinese
te
d
Government under the 1000 Talent Plan. Partial funding for this project was produced by a Grant-in-Aid for Scientific Research (B) from Japan Society for the
Ac ce p
Promotion of Science (JSPS) to Tsunogae (No. 26302009).
References
Amelin, Y., Lee, D.C., Halliday, A.N., 2000. Early-Middle Archaean crustal evolution deduced from Lu-Hf and U-Pb isotopic studies of single zircon grains. Geochim. Cosmochim. Acta 64, 4205-4225. Almond, D.C., 1991. Arena Gneiss and Kandy Gneiss-a proposed subdivision of the Highland Series around Kandy, and its significance. J. Geol. Soc. Sri Lanka 3, 41-50. Blichert-Toft, J., Albarede, F., 1997. The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle–crust system. Earth and Planetary Science Letters 148, 243-258.
Page 64 of 111
P a g e | 65
Baur, N, Kröner, A, Todt, W, Liew, T.C., Hofman, A.W., 1991. U-Pb isotopic systematics of zircons from prograde and retrograde transition zones in high grade orthogneisses, Sri Lanka. J. Geol. 99, 527-545.
ip t
Berger, A.R., Jayasinghe, N.R., 1976. Precambrian structure and chronology in the Highland Series of Sri Lanka. Precambrian Res. 3, 559-576.
cr
Burton, K.W., O’Nions, R.K., 1990. The timescale and mechanism of granulite formation at Kurunegala, Sri Lanka. Contrib. Mineral. Petrol. 106, 66-89.
us
Cawood, P.A., Buchan, C., 2007. Linking accretionary orogenesis with supercontinent assembly. Earth-Science Reviews 82, 217-256.
an
Chetty, T.R.K., Yellappa, T., Mohanty, D.P., Nagesh, P., Sivappa, V.V., Santosh, M., Tsunogae, T., 2012. Mega Sheath Fold of the Mahadevi Hills, Cauvery Suture Zone, Southern India: Implication for Accretionary Tectonics. J. Geol. Soc. India 80, 747-758.
M
Condie, K.C., Aster, R.C., 2013. Refinement of the supercontinent cycle with Hf, Refinement of the supercontinent cycle with Hf, Nd and Sr isotopes: Geoscience Frontiers 4, 667-680.
te
d
Condie, K.C., Kröner, A., 2013. The building blocks of continental crust: Evidence for a major change in the tectonic setting of continental growth at the end of the Archean: Gondwana Res. 23, 394-402.
Ac ce p
Cordani, U., Cooray, P.G., 1990. Rb-Sr ages of granulites and gneisses from the Precambrian of Sri Lanka. J. Geol. Soc. Sri Lanka 2, 35-43. Cooray, P.G., 1962. Charnockites and their associated gneisses in the Precambrian of Ceylon. J. Geol. Soc. London 118, 239– 273. Cooray, P.G., 1984. An introduction to the geology of Sri Lanka. National Museum Sri Lanka Publication, 340. Cooray, P.G., 1994. The Precambrian of Sri Lanka: a historic review. Precambrian Res. 66, 3-18. De Maesschalck, A.A, Oen, I.S., Hebeda, E.H., Verschure, R., Arps, E.S., 1990. Rb-Sr whole rock ages of Kataragama and Pottuvil Charnockites and east Vijayan gneiss: indication of a 2 Ga metamorphism in the highlands of south east Sri Lanka. J. Geol. 98, 772-779. Dahanayake, K., Jayasena, H.A.K., 1983. General geology and petrology of some Precambrian crystalline rocks from the Vijayan Complex of Sri Lanka.
Page 65 of 111
P a g e | 66 Precambrian Res. 19, 301–315. Dahanayake, K., 1982. Structural and petrological studies on the Precambrian Vijayan Complex of Sri Lanka. Revista Brasiliano Geociencias 12, 89–93.
cr
ip t
Dharmapriya, P.L., Malaviarachchi, S.P.K., Galli, A., Su, B., Subasinghe, N.D., Dissanayake, C.B., Nimalsiri, T.B., Zhu, B., 2014. P-T evolution of a spinel + quartz bearing khondalite from the Highland Complex, Sri Lanka: Implications for non-UHT metamorphism, J. Asian Earth Sci., Doi: http://dx.doi.org/10.1016/j.jseaes.2014.05.003.
us
Ellis, D.J., Green, D.D., 1979. An experimental study of the effect of Ca upon garnet-clinopyroxene Fe-Mg exchange equilibria. Contrib. Mineral. Petrol 71, 13-22.
an
Faulhaber, S., Raith, M., 1991. Geothermometry and geobarometry of high grade rocks: a case study on garnet pyroxene granulites in southern Sri Lanka. Min. Mag. 55, 33-56.
M
Ganguly, J., Cheng, W., Tirone, M., 1996. Thermodynamics of aluminosilicate garnet solid solution: new experimental data, an optimized model, and thermometric applications. Contrib. Mineral. Petrol. 126, 137-151.
te
d
Graham, C.M., Powell, R., 1984. A garnet-hornblende geothermometer: calibration, testing, and application to Pelona Schist, Southern California. J. Metamorph. Geol. 2, 13-31.
Ac ce p
Green, T.H., Ringwood, A.E., 1967. An experimental investigation of the gabbro to eclogite transformation and its petrological applications. Geochim. Cosmochim. Acta 31, 767-833. Hansen, E.C., Janardan, A.S., Newton, R.C., Prame, W.K.B.N., Ravindrakumar, G.R., 1987. Arrested charnockite formation in southern India and Sri Lanka. Contrib. Mineral. Petrol. 196, 66-89. Hapuarachchi, D.J.A.C., 1968. Cordierite and wollastonite-bearing rocks of south-west of Ceylon. Geol. Mag. 105, 317-324. Harley, S.L., 1989. The origin of granulites: a metamorphic perspective. Geol. Mag. 126, 215-247. Hatherton, T., Pattiarachchi, D.B., Ranasinghe, V.V.C., 1975. Gravity map of Sri Lanka, 1:1,000,000. Geol. Surv. Dept., Sri Lanka, Prof. Paper 3, 39. Hiroi, Y., Ogo, Y., Namba, L., 1994. Evidence for prograde metamorphic evolution of Sri Lankan politic granulites and implications for the development of
Page 66 of 111
P a g e | 67 continental crust. Precambrian Res. 66, 245-263.
ip t
Hiroi, Y., Yanagi, A., Kato, M., Kobayashi, T., Prame, B., Hokada, T., SatishKumar, M., Ishikawa, M., Adachi, T., Osanai, Y., Motoyoshi, Y., Shiraishi, K., 2014. Super cooled melt inclusions in lower-crustal granulites as a consequence of rapid exhumation by channel flow. Gondwana Res. 25, 226234.
cr
Holland, T., Blundy, J., 1994. Non-ideal interactions in calcic amphiboles and their bearing on amphibole-plagioclase thermometry. Contrib. Mineral. Petrol. 116, 433-447.
us
Hölzl, S., Köhler, H., Kröner, A., Jaeckel, P., Liew, T.C., 1991. Geochronology of the Sri Lankan basement. Geological Survey Department, Sri Lanka, Professional Paper 5, 237-257.
an
Hölzl, S., Hofmann, A.W., Todt, W., Köhler, H., 1994. U–Pb Geochronology of the Sri Lankan basement. Precambrian Res. 66, 123–149.
M
Hu, Z.C., Gao, S., Liu,Y.S., Hu, S.H., Chen, H.H., Yuan, H.L., 2008. Signal enhancement in laser ablation ICP-MS by addition of nitrogen in the central channel gas. Journal of Analytical Atomic Spectrometry 23, 1093–1101.
Ac ce p
te
d
Hu, Z.C., Liu,Y.S., Chen, L., Zhou, L., Li, M., Zong, K.Q., Zhu, L., Gao, S., 2011. Contrasting matrix induced elemental fractionation in NIST SRM and rock glasses during laser ablation ICPMS analysis at high spatial resolution. Journal of Analytical Atomic Spectrometry 26, 425–430. Kagami, H., Owada, M., Osanai, Y., Hiroi, Y., 1990. Preliminary geochronological study of Sri Lankan rocks. Interim Report of Japan-Sri Lanka Joint Research, p. 55-70. Karunaratne, D.N., Perera, L.R.K. and Fernando, G.W.A.R., 2002. The missing link between Prograde decompression (?) and Retrograde Isobaric Cooling in the Highland Complex Granulites of Sri Lanka. Proc. 18th Ann. Technical Sessions of the Geol. Soc. of Sri Lanka, p. 3 (Abstract). Kehelpannala, K.V.W., 1991. Structural evolution of high-grade terrains in Sri Lanka with special reference to the areas around Dodanhaslanda and Kandy. Geological Survey Department, Sri Lanka, Professional Paper 5, 69-88. Kehelpannala, K.V.W., 1997. Deformation of a high-grade Gondwana fragment, Sri Lanka. Gondwana Res. 1, 47–68. Kehelpannala, K.V.W., 1999. Shear-zone controlled arrested charnockitization,
Page 67 of 111
P a g e | 68 retrogression and metasomatism of high-grade rocks. Gondwana Res. 2, 573-577. Kehelpannala, K.V.W., 2003. Structural evolution of the middle to lower crust in Sri Lanka—a review. J. Geol. Soc. Sri Lanka 11, 45–86.
ip t
Kehelpannala, K.V.W., 2004. Arc accretion around Sri Lanka during the assembly of Gondwana. Gondwana Res. 7, 1323–1328.
cr
Kleinschrodt, R., Voll, G., Kehelpannala, W., 1991. A layered basic intrusion, deformed and metamorphosed in granulite facies of the Sri Lanka basement. Geol. Rundschau 80, 779-800.
us
Kleinschrodt, R., 1994. Large-scale thrusting in the lower crustal basement of Sri Lanka. Precambrian Res. 66, 39-57.
an
Köhler, H., Dratch, V., Holzl, S., Fehr, T., 1991. Geochronology of Sri Lanka basement. Part II: Rb-Sr and Sm-Nd dating on mineral and whole rock samples. Terra Abstracts 3, 504.
te
d
M
Koizumi, T., Tsunogae, T., Santosh, M., Tsutsumi, Y., Chetty, T.R.K., Saitoh, Y., 2014. Petrology and zircon U-Pb geochronology of metagabbros from a mafic-ultramafic suite at Aniyapuram: Neoarchean to Early Paleoproterozoic convergent margin magmatism and Middle Neoproterozoic high-grade metamorphism in southern India. J. Asian Earth Sci., doi: 10.1016/j.jseaes.2014.04.013.
Ac ce p
Korhonen, F.J., Saw, A.K., Clark, C., Brown, M., Bhattacharya, S., 2011. New constraints on UHT metamorphism in the Eastern Ghats Province through the application of phase equilibria modelling and in situ geochronology. Gondwana Res. 20, 764-781. Kriegsman, L.M., 1991. Structural geology of the Sri Lankan basement – a preliminary review. Geological Survey Department, Sri Lanka, Professional Paper 5, 52-68. Kriegsman, L., 1994. Evidence for a fold nappe in the high-grade basement of central Sri Lanka: terrane assembly in the Pan-African lower crust? Precambrian Res. 66, 59-76. Kriegsman, L., 1995. The Pan-African events in East Antarctica: a review from Sri Lanka and the Mozambique Belt. Precambrian Res. 75, 263–277. Kriegsman, L.M., Schumacher, J.C., 1999. Petrology of sapphirine-bearing and associated granulites from central Sri Lanka. J. Petrol. 40, 1211–1239.
Page 68 of 111
P a g e | 69 Kroner, A and Jaeckel, P., 1994. Zircon ages from rocks of the Wanni Complex, Sri Lanka. J.Geol. Soc. Sri Lanka 5, 41-57.
ip t
Kröner, A., Williams, I.S., Compston, W., Baur, N., Vithanage, P.W., Perera., L.R.K., 1987. Zircon ion microprobe dating of high grade rocks in Sri Lanka. J. Geol. 95, 775-791.
cr
Kröner, A., Williams, I.S., Compston, W., Baur, N., Vithanage, P.W., Perera., L.R.K., 1987. Zircon ion microprobe dating of high grade rocks in Sri Lanka. J. Geol. 95, 775-791.
an
us
Kröner, A., Cooray, P.G., Vitanage, P.W., 1991. Lithotectonic subdivision of the Precambrian basement in Sri Lanka. Geological Survey Department, Sri Lanka, Professional Paper 5, 5-21. Kröner, A., Williams, I.S., 1993. Age of metamorphism in the high-grade rocks of Sri Lanka. J. Geol. 101, 513–521.
te
d
M
Kröner, A., Jaeckel, P., Williams, I.S., 1994. Pb-loss patterns in zircons from a high-grade metamorphic terrain as revealed by different dating methods: UPb and Pb-Pb ages for igneous and metamorphic zircons from northern Sri Lanka. Precambrian Res. 66, 151-181.
Ac ce p
Kröner, A., Kehelpannala, K.V.W., Hegner, A., 2003. Ca. 750–1100 Ma magmatic events and Grenville-age deformation in Sri Lanka: relevance for Rodinia supercontinent formation and dispersal, and Gondwana amalgamation. J. Asian Earth Sci. 22, 279-300. Kröner, A., Rojas-Agramonte, Y., Kehelpannala, K. V. W., Zack, T., Hegner, E., Geng, H. Y., Wong, J., Barth, M., 2013. Age, Nd–Hf isotopes, and geochemistry of the Vijayan Complex of eastern and southern Sri Lanka: A Grenville-age magmatic arc of unknown derivation, Precambrian Res. 234, 288-321. Li, Z.X., Bogdanova, S.V., Collins, A.S., Davidson, A., De Waele, B., Ernst, R.E., Fitzsimons, I.C.W., Fuck, R.A., Gladkochub, D.P., Jacobs, J., Karlstrom, K.E., Lu, S., Natapov, L.M., Pease, V., Pisarevsky, S.A., Thrane, K., Vernikovsky, V., 2008. Assembly, configuration, and break-up history of Rodinia: a synthesis. Precambrian Res. 160, 179–210. Liu, Y.S., Hu, Z.C., Gao, S., Gunther, D., Xu, J., Gao, C.G., Chen, H.H., 2008. In situ analysis of major and trace elements of anhydrous minerals by LA-ICPMS without applying an internal standard. Chemical Geology 257, 34–43.
Page 69 of 111
P a g e | 70
Liu, Y.S., Gao, S., Hu, Z., Gao, C., Zong, K., Wang, D., 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen: U-Pb dating, Hf isotopes and trace elements in zircons of mantle xenoliths. Journal of Petrology 51, 537–571.
ip t
Ludwig, K.R., 2003. ISOPLOT 3.00: A Geochronological Toolkit for Microsoft Excel. Berkeley Geochronology Center, Special Publications 1a.
cr
Malaviarachchi, S., Takasu, A. 2005. Petrographic study on the high-grade metamorphic rocks from the Highland and Kadugannawa Complexes, central Sri Lanka. Geoscience Reports of Shimane University, Japan, 24, 31-46.
an
us
Malaviarachchi, S.P.K., Takasu, A., 2011a. Petrology of metamorphic rocks from the Highland and Kadugannawa Complexes, Sri Lanka. J. Geol. Soc. Sri Lanka 14 103-122. Malaviarachchi, S.P.K., Takasu, A., 2011b. Electron Microprobe dating of Monazites from Sri Lanka. J. Geol. Soc. Sri Lanka 14, 81-90.
d
M
Mathavan, V, W. K. B. N., Prame, W.K.B.N., Cooray, P.G., 1999. Geology of the high grade Proterozoic terrains of Sri Lanka and the assembly of Gondwana: an update on recent developments, Gondwana Res. 2, 237-250.
te
Mathavan, V., Fernando G.W.A.R., 2001. Reactions and textures in grossularwollastonite-scapolite calc-silicate granulites from Maligawila, Sri Lanka. Lithos, 59, 217-232.
Ac ce p
McDonough, W.F., Sun, S.-S., 1995. The composition of the Earth. Chemical Geology 120, 223-253. Milisenda, C.C., 1991. Compositional characteristics of the Vijayan Complex. Geological Survey Department, Sri Lanka, Professional Paper 5, 135-140. Milisenda, C.C., Pohl, J.R., Hofmann, A.W., 1991. Charnockite formation at Kurunegala, Sri Lanka. Geological Survey Department, Sri Lanka, Professional Paper 5, 141-149. Milisenda, C.C, Leiw, T.C., Hofmann, A.W., Kröner, A., 1988. Isotopic mapping of age provinces in Precambrian high grade terrains: Sri Lanka. J. Geol. 96, 608-615. Milisenda, C.C., Leiw, T.C., Hofmann, A.W., Köhler, H., 1994. Nd isotopic mapping of the Sri Lanka basement: update, and additional constraints from Sr isotopes. Precambrian Res. 66, 95-110.
Page 70 of 111
P a g e | 71 Moecher, D.P., Essene, E.J., Anovitz, L.M., 1988. Calculation and application of clinopyroxene-garnet-plagioclase-quartz geobarometers. Contrib. Mineral. Petrol. 100, 92–106.
ip t
Osanai, Y., Sajeev, K., Owada, M., Kehelpannala, K.V.W., Prame, W.K.B. Nakano, N., Jayatileke, S., 2006. Metamorphic evolution of ultrahightemperature and high-pressure granulites from Highland Complex, Sri Lanka. J. Asian Earth Sci. 28, 20-37.
cr
Pathirana, H.D.N.C., 1980. Geology of Sri Lanka in relation to plate tectonics. J. National Science Council Sri Lanka. 8, 75-85.
us
Pearce, J.A., Harris, N.B.W., Tindle, A.G., 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. J. Petrol. 25, 956– 983.
an
Perchuk, L.L., Safonov, O.G., Gerya, T.V., Fu, B., Harlov, D.E., 2000. Mobility of components in metasomatic transformation and partial melting of gneisses: an example from Sri Lanka. Contrib. Mineral. Petrol. 140, 212-232.
M
Perkins, D., Newton, R.C., 1981. Charnockite geobarometers based on coexistence garnet-pyroxene-plagioclase-quartz. Nature 292, 144–146.
te
d
Perera, L.R.K., 1983. The origin of pink granites of Sri Lanka – another view. Precambrian Res. 20, 17-37.
Ac ce p
Perera, L.R.K., 1984. Co-existing cordierite-almandine-a key to the metamorphic history of Sri Lanka. Precambrian Res. 25, 349-364. Perera, L.R.K., 1987, Petrogenesis of granulite-facies metamorphic rocks in Sri Lanka. University of Peradeniya (Mphil. thesis) (unpubl.). Perera, L.R.K., Kagami, H., 2011. Centimetre- and metre-scale Nd and Sr isotopic homogenization in Kadugannawa Complex, Sri Lanka. J. Geol. Soc. Sri Lanka 14, 129-141. Pohl, J.G., Emmermann, R., 1991. Chemical composition of the Sri Lanka Precambrian Basement. Geological Survey Department, Sri Lanka, Professional Paper 5, 94-123. Prame, W.K.B.N, 1991a. Metamorphism and nature of the granulite facies crust in south western Sri Lanka: Characterization by pelitic/ psammopelitic rocks and associated granulites. Geological Survey Department, Sri Lanka, Professional Paper 5, 188-200. Prame, W.K.B.N., 1991b. Petrology of the Kataragama Complex, Sri Lanka:
Page 71 of 111
P a g e | 72 evidence for high pressure granulite facies metamorphism and subsequent isobaric cooling. Geological Survey Department, Sri Lanka, Professional Paper 5, 200-224.
ip t
Prame, W.K.B.N., 1997. Geochemistry and genesis of granitoid rocks from southern Sri Lanka. Min. Petrol. 60, 245–265.
cr
Raase, P., Schenk, V., 1994. Petrology of granulite facies metapelites of the Highland Complex, Sri Lanka: implication for the metamorphic zonation and the P-T path. Precambrian Res. 66, 265-294.
an
us
Saitoh, Y., Tsunogae, T., Santosh, M., Chetty, T.R.K., Horie, K., 2011. Neoarchean high-pressure metamorphism from the northern margin of the Palghat-Cauvery Suture Zone, southern India: petrology and zircon SHRIMP geochronology. J. Asian Earth Sci. 42, 268-285.
M
Sajeev, K., Osanai, Y., 2004a. Ultrahigh-temperature metamorphism (1150° C,12 kbar) and multistage evolution of Mg-, Al-rich granulites from the Central Highland Complex Sri Lanka. J. Petrol. 45, 1821–1844.
te
d
Sajeev, K., Osanai, Y., 2004b. 'Osumilite' and 'spinel+quartz' from Highland Complex, Sri Lanka: a case of cooling and decompression after ultrahightemperature metamorphism. J. Min. Petrol. Sci. 99, 320–327.
Ac ce p
Sajeev, K., Osanai, Y., 2005. Thermal gradients in the Sri Lankan granulite terrane: a garnet–biotite thermometric approach, J. Metamorph. Geol. 23, 383-397. Sajeev, K., Osanai, Y., Suzuki, S., Kagami, H., 2003. Geochronological evidence for multistage metamorphic events in ultra-high temperature granulites from central Highland Complex, Sri Lanka. Polar Geosci. 16, 137-148. Sajeev, K., Osanai, Y., Connolly, J.A.D., Suzuki, S., Ishioka, J., Kagami, H., Rino, S., 2007. Extreme Crustal Metamorphism during a Neoproterozoic Event in Sri Lanka: A study of dry mafic granulites. J. Geol. 115, 563–582. Sajeev, K., Windley, B.F., Connolly, J.A.D., Kon, Y., 2009. Retrogressed eclogite (20 kbar, 1020ºC) from the Neoproterozoic Palghat-Cauvery suture zone, southern India. Precambrian Res. 171, 23-36. Sajeev, K., Williams, I.S., Osanai, Y., 2010. Sensitive high-resolution ion microprobe U–Pb dating of prograde and retrograde ultrahigh-temperature metamorphism as exemplified by Sri Lankan granulites. Geology 38, 971– 974.
Page 72 of 111
P a g e | 73
Sandiford, N., Powell, R., Martin, S.F., Perera, L.R.K., 1988. Thermal and baric evolution of garnet granulites from Sri Lanka. J. Metamorph. Geol. 6, 351364.
cr
ip t
Santosh, M., Yoshida, M., Nanda-Kumar, V., 1991. Fluid characteristics across a gneiss-charnockite reaction front in Sri Lanka: Implication for granulite formation in Gondwana deep crust. J. Mineral. Petrol. Econnomic Geol. 86, 27-44.
us
Santosh, M., Tsunogae, T., Shimizu, H., Dubessy, J., 2010. Fluid characteristics of retrogressed eclogites and mafic granulites from the Cambrian Gondwana suture zone in southern India. Contrib. Mineral. Petrol. 159, 349-369.
an
Santosh, M., Maruyama, S., Sato, K., 2009. Anatomy of a Cambrian suture in Gondwana: Pacific-type orogeny in southern India? Gondwana Res. 16, 321341.
d
M
Santosh, M., Yang, Q.Y., Shaji, E., Tsunogae, T., Ram Mohan, M., Satyanarayanan, M., 2013a. An exotic Mesoarchean microcontinent: The Coorg Block, southern India. Gondwana Res., DOI: 10.1016/j.gr.2013.10.005
Ac ce p
te
Santosh, M., Liu, S.J., Tsunogae, T., Li, J.H., 2013b. Paleoproterozoic ultrahightemperature granulites in the North China Craton: Implications for tectonic models on extreme crustal metamorphism. Precambrian Res. 222-223, 77106. Schenk, V., Raase, P., Schumacher, R., 1988. Very high temperatures and isobaric cooling before tectonic uplift in the Highland Series of Sri Lanka. Terra Cognita 8, 265. Schenk, V., Raase, P., Schumacher, R., 1991. Metamorphic zonation and P–T history of the Highland Complex in Sri Lanka. Geological Survey Department, Sri Lanka, Professional Paper 5, pp. 150–163. Schumacher, R., Shenk, V., Raase, P., Vithanage, P.W., 1990. Granulite facies metamorphism of metabasic and intermediate rocks in the Highland Series of Sri Lanka. In: Ashworth, J.R., Brown, M. (Eds.), High temperature metamorphism and Crustal Anatexis. The Mineralogical Society series 2, 235271. Schumacher ,R., Faulhaber, S., 1994. Summary and discussion of P-T estimates from garnet-pyroxene-plagioclase-quartz bearing granulite facies rocks from Sri Lanka. Precambrian Res. 66, 295-308.
Page 73 of 111
P a g e | 74
Shiraishi, K., Ellis, D.J. ,Hiroi, Y., Fanning, C.M., Motoyoshi, Y., Nakai, Y., 1994. Cambrian orogenic belt in east Antarctica and Sri Lanka: implications for Gondwana assembly. J. Geol. 102, 47-65.
ip t
Stosch, H.G., 1991. Geochemistry of the mafic intrusion in the synforms of the Kandy area. Geological Survey Department, Sri Lanka, Professional Paper 5, 125-134.
us
cr
Sun, S.-S., McDonough, W.F., 1989. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In Saunders, A.D., Norry M.J. (Eds.), Magmatism in the Ocean Basins. Geol. Soc. London, Spec. Publ. 42, 313-345
an
Tani, Y., Yoshida, M., 1996. The structural evolution of the Arena Gneisses and its bearing on Proterozoic tectonics of Sri Lanka. J. Southeast Asian Earth Sci. 14, 309-329.
M
Vervoort, J.D., Patchett, P.J., Gehrels, G.E., Nutman, A. P., 1996. Constraints on early Earth differentiation from hafnium and neodymium isotopes. Nature 379, 624-627.
te
d
Voll, G., Evangelakakisb, C., Kroll, H., 1994. Revised two-feldspar geothermometry applied to Sri Lankan feldspars. Precambrian Res. 66, 351367.
Ac ce p
Voll, G., Kleinschrodt, R., 1991. Sri Lanka: Structural, magmatic and metamorphic development of a Gondwana fragment. Geological Survey Department, Sri Lanka, Professional Paper 5, 22-52. Vitanage, P.W., 1972. Post-Precambrian uplift and regional neotectonic movements in Ceylon. 24th Int. Geol. Congr., Montreal, Sect. 3, 642-654. Vitanage, P.W., 1985. Tectonics and mineralization in Sri Lanka. Bull. Geol. Soc. Finland 57, 157-168. Wiedenbeck, M., Alle, P., Corfu, F., Griffin, W.L., Meier, M., Oberli, F., Quadt, A.V., Roddick, J.C., Spiegel, W., 1995. Three natural zircon standards for U-Th-Pb, Lu-Hf, trace element and REE analyses. Geostandards and Geoanalytical Research 19, 1-23. Willbold, M., Hegner, E., Kleinschrodt, R., Stosch, H.-G., Kehelpannala, K.V.W., Dulski, P., 2004. Geochemical evidence for a Neoproterozoic magmatic continental margin in Sri Lanka-relevence for the Rodinia-Gondwana
Page 74 of 111
P a g e | 75 supercontinental cycle. Precambrian Res. 130, 185-198. Xiao, W.J., Santosh, M., 2014. The western Central Asian Orogenic Belt: A window to accretionary orogenesis and continental growth. Gondwana Res. 25, 1429-1444.
an
us
cr
ip t
Yoshida, M., Funaki, M., Vitanage, P.W., 1992. Proterozoic to Mesozoic East Gondwana: The juxtaposition of India, Sri Lanka, and Antarctica. Tectonics 11, 381-391.
M
Figure captions
Fig. 1 Generalized geological and tectonic framework of Sri Lanka showing the
d
major crustal blocks and their boundaries (after Cooray, 1994). The study area is
te
shown in the Fig. 2 is marked by box. Sample localities are shown by open
Ac ce p
circles. (Sanjeewa)
Fig. 2 Detailed geological map and sample locations of the study areas (Figure modified after Kröner et al. 2003).
Fig. 3 Representative field photographs of locations KDG-4, KDG-1, KDG-2 from the Kadugannawa Complex and WC-2 from the Wanni Complex.
(a)
Metagranodiorite;
(d)
(b)
metadiorite;
(c)
garnet-bearing
amphibolite;
metagranodiorite with patches of incipient charnockite. See text for discussion.
Page 75 of 111
P a g e | 76 Fig. 4 Representative field photographs of locations HC-1, HC-11, HC-7 and HC10 from the Highland Complex. (a) Blocks and boudins of mafic granulite within
garnet-bearing charnockite. See text for discussion.
ip t
marble; (b) mafic sills in metasediments; (c) garnet-bearing metagabbro; (d)
cr
Fig. 5. Photomicrographs showing textures of representative samples. (a) A
us
typical granoblastic texture of hornblende + biotite + quartz + plagioclase + Kfeldspar in metagranodiorite (sample WC-2). (b) Metagranodiorite from the
an
Kadugannawa Complex (sample KDG-4). (c) Orthopyroxene-bearing and Kfeldspar-absent metadiorite (sample KDG-1). (d) Porphyroblastic garnet in the
M
matrix of fine-grained hornblende and plagioclase in garnet amphibolite (sample
d
KDG-2). (e) Medium-grained foliated pyroxene + hornblende + plagioclase
te
assemblage in mafic granulite (sample HC-1). (f) Equigranular texture of
Ac ce p
pyroxene + hornblende + plagioclase in mafic granulite (sample HC-11).
Fig. 5 (continued). (g) Orthopyroxene + plagioclase symplectite present as grayish spots in sample HC-11 as an evidence of decomposition of garnet during exhumation. (h) Coarse-grained massive metagabbro with subidioblastic garnet, clinopyroxene, and minor orthopyroxene (sample HC-7). Plagioclase and ilmenite fill the matrix of the minerals. (i) Medium-grained garnet and quartz which are separated by orthopyroxene + plagioclase corona in metagabbro (sample HC-8). (j) Granoblastic plagioclase + clinopyroxene + orthopyroxene + hornblende with symplectitic aggregates of orthopyroxene + plagioclase garnet charnockite
Page 76 of 111
P a g e | 77 (sample HC-10).
Fig. 6. Compositional diagrams showing chemistry of representative minerals. (a) (b)
Compositional
diagrams
showing
Ca/(Fe+Mg+Ca+Mn)
and
ip t
and
Mn/(Fe+Mg+Ca+Mn) versus XMg of garnet. (c) Wollastonite-enstatite-ferrosilite
cr
diagram showing compositions of clinopyroxene and clinopyroxene. (d) Anorthite-
us
albite-orthoclase diagram showing compositions of plagioclase. (e) Si (pfu) versus XMg diagram showing compositions of calcic amphibole. (f) XMg versus
an
TiO2 (wt.%) diagram showing biotite chemistry. GR: metagranodiorite, GD:
M
metadiorite, GA: garnet amphibolite, MG: mafic granulite, MGB: metagabbro.
d
Fig. 7. Major element variation diagrams showing concentrations of major
te
elements in the examined samples.
Ac ce p
Fig. 8. (a) to (c) Primitive mantle-normalized multi-element variation diagrams for magnesian (a), ferroan (b), and felsic (c) groups. Normalizing values are from Sun and McDonough (1989). (d) to (f) Chondrite-normalized REE spider plots for magnesian (d), ferroan (e), and felsic (f) groups. Normalizing values are from McDonough and Sun (1995).
Fig. 9. Discrimination diagrams for granites based on Rb-Y-Nb and Rb-Yb-Ta variations after Pearce et al. (1984). (a) Nb-Y diagram. (b) Ta-Yb diagram. (c) RbY+Nb diagram. (d) Rb-Yb+Ta diagram. VAG: volcanic-arc granites, syn-COLG:
Page 77 of 111
P a g e | 78 syn-collisional granites, WPG: within plate granites, ORG: ocean-ridge granites
Fig. 10. Triangular diagrams showing compositions of mafic granulites and
ip t
amphibolite. (a) TiO2-MnOx10-P2O5x10 diagram. OIT: oceanic-island tholeiite or seamount tholeiite, OIA: oceanic-island alkali basalt or seamount alkali basalt,
cr
CAB: island-arc calc-alkaline basalt, IAT: island-arc tholeiite, Bon: boninite. (b)
us
2Nb-Zr/4-Y diagram. AI: within-plate alkali basalt, AII: within-plate alkali basalt and within-plate alkali tholeiite, B: E-type MORB, C: within-plate tholeiite and
an
volcanic-arc basalt, D: N-type MORB and volcanic-arc basalt. (c) Y/15-La/10Nb/8 diagram. 1A: calc-alkali basalt, 1C: volcanic arc tholeiite, 1B: an area of
M
overlap between 1A and 1C. 2A: continental basalt, 2B: back-arc basin basalt,
d
3A: alkali basalt from intercontinental rift, 3B: enriched E-MORB, 3C: weakly
te
enriched E-MORB, 3D: N-MORB.
Ac ce p
Fig. 11 Cathodoluminescence (CL) images of representative zircons in samples HC -1, HC- 7, and HC- 10. Analytical points (white circles for U-Pb and red circles for Lu-Hf), ages in Ma and εHf (t) values are also shown.
Fig. 12 Cathodoluminescence (CL) images of representative zircons in samples HC-11, KDG-1, and KDG-2. Analytical points (white circles for U-Pb and red circles for Lu-Hf), ages in Ma and εHf (t) values are also shown.
Fig. 13 Cathodoluminescence (CL) images of representative zircons in samples
Page 78 of 111
P a g e | 79 KDG-4 and WC-2. Analytical points (white circles for U-Pb and red circles for LuHf), ages in Ma and εHf (t) values are also shown.
ip t
Fig. 14 Zircon U-Pb concordia plots, age data histograms and probability curves
cr
for samples HC-1, HC-7.
us
Fig. 15 Zircon U-Pb concordia plots, age data histograms and probability curves
an
for sample HC-10, HC-1.
Fig. 16 Zircon U-Pb concordia plots, age data histograms and probability curves
d
M
for samples KDG-1, KDG-2, KDG-4.
te
Fig. 17 Zircon U-Pb concordia plots, age data histograms, probability curves and bar charts for sample WC-2.
Ac ce p
Fig. 18 REE patterns of zircons in samples HC-1, HC- 7, HC-10, HC-11. Fig. 19 REE patterns of zircons in samples KDG-1, KDG-2, KDG-4, WC-2. Fig. 20 εHf (t) versus
207
Pb/206Pb age diagram of zircons from the magmatic
suites of Sri Lanka.
Fig. 21 Plate tectonic cartoons illustrating the model of Neoproterozoic doublesided subduction (a) culminating in Cambrian collision (b) proposed in this study. See text for details.
Page 79 of 111
P a g e | 80 Table captions
Table 1 Location, rock type and mineralogy of samples from the crustal blocks in
ip t
Sri Lanka analysed in this study.
Table 8 LA-MC-ICPMS Lu-Hf isotope data on zircons from the magmatic suites
Ac ce p
te
d
M
an
us
cr
of Sri Lanka.
Page 80 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Graphical Abstract (for review)
Page 81 of 111
Table-1
HC-11 WC-2
cr
metagabbro
Cpx Grt Opx Pl Ilm Mag Hbl Zrn
o
o
metagabbro
Pl Cpx Grt Qtz Opx Ilm Mag Rt Zrn
o
o
garnet charnockite
Pl Qtz Cpx Opx Hbl Ilm Py Cal Zrn Bt
o
o
mafic granulite
Pl Cpx Opx Hbl Ilm Zrn
o
o
metagranodiorite
Pl Kfs Qtz Bt Hbl Ilm Zrn Ap
N 7 15' 59.9"; E 80 37' 58.4"
Hantana
N 7 15' 00.1"; E 80 37' 59.0"
Hantana
N 7 12' 22.3"; E 80 47' 54.9"
Hakurutale 2km SE from Kivullinda JCT Udadigana
Pl Cpx Hbl Opx Ilm Mag Zrn
o
us
HC-10
mafic granulite
o
N 7 14' 21.0"; E 80 43' 48.1" N 7 30' 09.3"; E 80 20' 43.3"
an
HC-8
Pl Kfs Hbl Bt Qtz Mag Zrn
o
M
HC-7
metagranodiorite
o
N 7 18' 07.4"; E 80 44' 12.1"
Diagana
Hbl Pl Grt Ilm Mag Ap Zrn Qtz
ed
HC-1
garnet amphibolite
pt
KDG-4
o o Mahaberiyatenna Airport site N 7 16' 36.8"; E 80 45' 34.3" N 7o 22' 47.2"; E 80o 31' 51.8" Maussawa
Ac ce
KDG-2
ip t
Table 1 Location, rock type and mineralogy of samples from the crustal blocks in Sri Lanka analysed in this study. U-Pb geochronology and Lu-Hf isotope analysis Sample No. Locality Co-ordinates Rock type Assemblage o o N 7 16' 44.1"; E 80 40' 49.1" KDG-1 Kundasale metadiorite Pl Hbl Bt Qtz Ilm Opx Py Ap Zrn
Page 82 of 111
Yb/177Hf
176
Lu/177Hf
176
Hf/177Hf
WC-2-01
608.0
0.034894
0.001226
0.282199
WC-2-02
608.0
0.021266
0.000788
0.282249
WC-2-03
608.0
0.020237
0.000744
0.282257
WC-2-04
608.0
0.018627
0.000675
0.282168
WC-2-05
608.0
0.027675
0.000978
WC-2-06
608.0
0.016021
WC-2-07
608.0
WC-2-08 WC-2-09
1s
176
Hf/177Hfi
εHf(0)
εHf(t)
TDM (Ma)
TDMC (Ma)
fLu/Hf
0.282185
-20.3
-7.4
1495
2018
-0.96
0.000011
0.282240
-18.5
-5.4
1409
1895
-0.98
0.000011
0.282249
-18.2
-5.1
1395
1875
-0.98
0.000010
0.282161
-21.3
-8.2
1516
2071
-0.98
0.282292
0.000009
0.282281
-17.0
-4.0
1355
1803
-0.97
0.000569
0.282262
0.000012
0.282256
-18.0
-4.9
1382
1860
-0.98
0.015693
0.000541
0.281985
0.000012
0.281979
-27.8
-14.7
1762
2475
-0.98
608.0
0.027203
0.000977
0.282197
0.000012
0.282186
-20.3
-7.3
1487
2014
-0.97
608.0
0.028383
0.001019
0.282255
0.000010
0.282243
-18.3
-5.3
1409
1888
-0.97
WC-2-10
608.0
0.028674
0.000987
0.282172
0.000011
0.282161
-21.2
-8.2
1523
2070
-0.97
WC-2-12
608.0
0.018053
0.000683
0.282233
0.000010
0.282225
-19.1
-6.0
1427
1928
-0.98
WC-2-13
608.0
0.035602
0.001250
0.282257
0.000011
0.282243
-18.2
-5.3
1414
1888
-0.96
WC-2-14
608.0
0.014895
0.000568
0.282273
0.000008
0.282267
-17.6
-4.5
1367
1834
-0.98
WC-2-15
608.0
0.045200
0.001550
0.282143
0.000014
0.282125
-22.2
-9.5
1587
2149
-0.95
WC-2-16
608.0
0.019193
0.000718
0.282256
0.000012
0.282248
-18.2
-5.1
1396
1876
-0.98
WC-2-17
608.0
0.029196
0.001037
0.282291
0.000011
0.282279
-17.0
-4.0
1359
1807
-0.97
WC-2-18
608.0
0.020072
0.000761
0.282205
0.000010
0.282196
-20.1
-7.0
1469
1993
-0.98
608.0
0.033156
0.001157
0.282253
0.000011
0.282240
-18.4
-5.4
1417
1895
-0.97
608.0
0.019062
0.000709
0.282288
0.000010
0.282280
-17.1
-4.0
1351
1805
-0.98
608.0
0.015912
0.000597
0.282251
0.000009
0.282244
-18.4
-5.3
1398
1885
-0.98
608.0
0.030750
0.001084
0.282255
0.000010
0.282243
-18.3
-5.3
1410
1887
-0.97
608.0
0.013215
0.000543
0.282282
0.000009
0.282276
-17.3
-4.2
1354
1814
-0.98
WC-2-24
608.0
0.018073
0.000667
0.282273
0.000011
0.282265
-17.6
-4.5
1371
1838
-0.98
WC-2-25
608.0
0.013856
0.000551
0.282289
0.000009
0.282283
-17.1
-3.9
1344
1799
-0.98
WC-2-26
608.0
0.021998
0.000825
0.282281
0.000012
0.282272
-17.4
-4.3
1365
1824
-0.98
WC-2-20 WC-2-21 WC-2-22 WC-2-23
d
ep te
Ac c
WC-2-19
an
0.000009
M
Age (Ma)
us
176
No.
cr
ip t
Table-8
Page 83 of 111
ip t 0.001009
0.282174
WC-2-28
608.0
0.025933
0.000933
0.281979
WC-2-29
608.0
0.017038
0.000656
0.282261
WC-2-31
608.0
0.016344
0.000586
0.282266
WC-2-32
608.0
0.024456
0.000898
0.282258
WC-2-33
608.0
0.018576
0.000708
0.282240
HC-11-01
539.1
0.000998
0.000031
HC-11-02
539.1
0.000446
0.000013
HC-11-04
539.1
0.000477
0.000014
HC-11-05
539.1
0.000516
HC-11-06
539.1
HC-11-07
0.000011
cr
0.027214
0.282163
us
608.0
-21.1
-8.2
1521
2067
-0.97
0.281968
-28.0
-15.0
1789
2498
-0.97
0.000010
0.282254
-18.1
-4.9
1387
1864
-0.98
0.000010
0.282259
-17.9
-4.8
1378
1852
-0.98
0.000015
0.282247
-18.2
-5.2
1400
1878
-0.97
0.000012
0.282232
-18.8
-5.7
1417
1912
-0.98
0.282058
0.000012
0.282057
-25.3
-13.4
1641
2346
-1.00
0.282028
0.000010
0.282028
-26.3
-14.5
1681
2412
-1.00
0.282059
0.000010
0.282059
-25.2
-13.4
1639
2343
-1.00
0.000015
0.282043
0.000010
0.282043
-25.8
-13.9
1660
2377
-1.00
0.001572
0.000054
0.282058
0.000013
0.282057
-25.3
-13.4
1642
2346
-1.00
539.1
0.000503
0.000015
0.282061
0.000010
0.282060
-25.2
-13.3
1636
2339
-1.00
HC-11-08
539.1
0.000441
0.000013
0.282057
0.000010
0.282057
-25.3
-13.4
1641
2346
-1.00
HC-11-09
539.1
0.000124
0.000003
0.282004
0.000010
0.282004
-27.2
-15.3
1712
2464
-1.00
HC-11-10
539.1
0.000987
0.000031
0.282073
0.000012
0.282072
-24.7
-12.9
1621
2313
-1.00
HC-11-11
539.1
0.000487
0.000017
0.282081
0.000010
0.282081
-24.4
-12.6
1608
2293
-1.00
HC-11-12
539.1
0.000689
0.000020
0.282071
0.000010
0.282071
-24.8
-12.9
1622
2315
-1.00
HC-11-14
539.1
0.000309
0.000009
0.282057
0.000009
0.282057
-25.3
-13.4
1641
2346
-1.00
539.1
0.000939
0.000029
0.282067
0.000010
0.282067
-24.9
-13.1
1628
2324
-1.00
539.1
0.000741
0.000021
0.282079
0.000010
0.282079
-24.5
-12.7
1611
2298
-1.00
539.1
0.000060
0.000002
0.281999
0.000010
0.281999
-27.3
-15.5
1719
2475
-1.00
539.1
0.000902
0.000028
0.282095
0.000011
0.282095
-23.9
-12.1
1590
2263
-1.00
539.1
0.000130
0.000004
0.282024
0.000011
0.282024
-26.5
-14.6
1686
2421
-1.00
HC-11-20
539.1
0.001070
0.000035
0.281947
0.000011
0.281947
-29.2
-17.3
1791
2591
-1.00
HC-11-21
539.1
0.000737
0.000022
0.282072
0.000010
0.282072
-24.8
-12.9
1621
2314
-1.00
HC-11-22
539.1
0.009487
0.000325
0.281860
0.000010
0.281857
-32.2
-20.5
1922
2790
-0.99
HC-11-16 HC-11-17 HC-11-18 HC-11-19
d
ep te
Ac c
HC-11-15
an
0.000012
M
WC-2-27
Page 84 of 111
ip t 0.000496
0.282002
HC-10-03
533.0
0.000208
0.000006
0.282106
HC-10-04
533.0
0.004298
0.000152
0.282137
HC-10-05
533.0
0.004251
0.000164
0.282015
HC-10-06
533.0
0.000561
0.000015
0.282145
HC-10-07
533.0
0.004830
0.000186
HC-10-08
533.0
0.006191
HC-10-09
533.0
0.004883
HC-10-11
533.0
0.000161
0.000005
HC-10-12
533.0
0.012997
HC-10-13
533.0
HC-10-14
0.000019
cr
0.013494
0.281997
us
533.0
-27.2
-15.7
1737
2483
-0.99
0.282106
-23.6
-11.8
1575
2242
-1.00
0.000013
0.282136
-22.4
-10.8
1537
2175
-1.00
0.000010
0.282013
-26.8
-15.1
1705
2447
-1.00
0.000011
0.282144
-22.2
-10.5
1522
2156
-1.00
0.281993
0.000011
0.281991
-27.6
-15.9
1736
2497
-0.99
0.000241
0.281993
0.000011
0.281990
-27.6
-15.9
1738
2498
-0.99
0.000178
0.281933
0.000019
0.281931
-29.7
-18.0
1816
2629
-0.99
0.282017
0.000011
0.282017
-26.7
-15.0
1694
2438
-1.00
0.000450
0.281765
0.000022
0.281760
-35.6
-24.1
2059
3007
-0.99
0.002461
0.000099
0.282142
0.000018
0.282141
-22.3
-10.6
1529
2163
-1.00
533.0
0.008187
0.000293
0.282114
0.000018
0.282111
-23.3
-11.6
1575
2229
-0.99
HC-10-15
533.0
0.000393
0.000011
0.282120
0.000010
0.282119
-23.1
-11.4
1556
2212
-1.00
HC-10-16
533.0
0.006434
0.000236
0.281929
0.000015
0.281927
-29.8
-18.2
1824
2639
-0.99
HC-10-17
533.0
0.004681
0.000178
0.281927
0.000011
0.281925
-29.9
-18.3
1825
2644
-0.99
HC-10-18
533.0
0.034608
0.001208
0.281676
0.000013
0.281664
-38.8
-27.5
2222
3218
-0.96
HC-10-19
533.0
0.005998
0.000237
0.281502
0.000016
0.281500
-44.9
-33.3
2401
3580
-0.99
HC-10-20
533.0
0.007711
0.000309
0.281694
0.000013
0.281691
-38.1
-26.5
2147
3160
-0.99
533.0
0.030751
0.001051
0.281673
0.000011
0.281663
-38.9
-27.5
2217
3221
-0.97
533.0
0.008653
0.000329
0.281826
0.000012
0.281823
-33.5
-21.9
1969
2870
-0.99
533.0
0.018805
0.000679
0.281850
0.000011
0.281844
-32.6
-21.1
1953
2823
-0.98
533.0
0.002831
0.000107
0.281977
0.000010
0.281976
-28.1
-16.5
1754
2531
-1.00
HC-10-22 HC-10-23 HC-10-24 HC-10-25
d
ep te
Ac c
HC-10-21
an
0.000011
M
HC-10-01
533.0
0.004084
0.000161
0.281972
0.000010
0.281970
-28.3
-16.6
1763
2543
-1.00
HC-10-26
533.0
0.003430
0.000131
0.282035
0.000010
0.282034
-26.0
-14.4
1675
2401
-1.00
HC-10-27
533.0
0.003593
0.000141
0.282029
0.000011
0.282028
-26.3
-14.6
1684
2415
-1.00
HC-7-01
543.0
0.001018
0.000031
0.282260
0.000011
0.282260
-18.1
-6.2
1366
1892
-1.00
Page 85 of 111
ip t 0.000022
0.282279
HC-7-05
543.0
0.001113
0.000033
0.282251
HC-7-06
543.0
0.000665
0.000018
0.282254
HC-7-07
543.0
0.000549
0.000016
0.282262
HC-7-08
543.0
0.001852
0.000056
0.282280
HC-7-09
543.0
0.000577
0.000017
0.282255
HC-7-11
543.0
0.000873
0.000024
HC-7-13
543.0
0.000688
0.000020
HC-7-14
543.0
0.000992
0.000030
HC-7-15
543.0
0.001226
HC-7-16
543.0
HC-1-01
0.000010
cr
0.000762
0.282279
us
543.0
-17.4
-5.5
1340
1850
-1.00
0.282251
-18.4
-6.5
1378
1912
-1.00
0.000010
0.282254
-18.3
-6.4
1374
1905
-1.00
0.000011
0.282262
-18.0
-6.1
1363
1888
-1.00
0.000013
0.282279
-17.4
-5.5
1340
1849
-1.00
0.000011
0.282254
-18.3
-6.4
1373
1904
-1.00
0.282280
0.000010
0.282280
-17.4
-5.5
1338
1847
-1.00
0.282255
0.000010
0.282255
-18.3
-6.4
1372
1903
-1.00
0.282233
0.000012
0.282232
-19.1
-7.1
1403
1953
-1.00
0.000040
0.282236
0.000012
0.282236
-19.0
-7.0
1399
1946
-1.00
0.000874
0.000026
0.282221
0.000012
0.282221
-19.5
-7.5
1418
1978
-1.00
553.0
0.005983
0.000268
0.282368
0.000014
0.282365
-14.3
-2.2
1227
1651
-0.99
HC-1-02
553.0
0.006039
0.000256
0.282374
0.000010
0.282371
-14.1
-2.0
1218
1637
-0.99
HC-1-03
553.0
0.006708
0.000282
0.282373
0.000016
0.282370
-14.1
-2.0
1220
1639
-0.99
HC-1-04
553.0
0.007236
0.000298
0.282375
0.000014
0.282372
-14.1
-2.0
1218
1636
-0.99
HC-1-05
553.0
0.005477
0.000231
0.282380
0.000015
0.282378
-13.8
-1.8
1208
1621
-0.99
HC-1-06
553.0
0.005387
0.000228
0.282395
0.000014
0.282392
-13.3
-1.3
1188
1589
-0.99
HC-1-07
553.0
0.003809
0.000138
0.282417
0.000013
0.282416
-12.5
-0.4
1155
1537
-1.00
553.0
0.005521
0.000238
0.282390
0.000015
0.282388
-13.5
-1.4
1195
1600
-0.99
553.0
0.005921
0.000250
0.282393
0.000011
0.282390
-13.4
-1.3
1192
1594
-0.99
553.0
0.004044
0.000142
0.282416
0.000015
0.282414
-12.6
-0.5
1157
1540
-1.00
553.0
0.006902
0.000294
0.282427
0.000011
0.282424
-12.2
-0.1
1146
1519
-0.99
HC-1-09 HC-1-10 HC-1-11 HC-1-12
d
ep te
Ac c
HC-1-08
an
0.000012
M
HC-7-04
553.0
0.005637
0.000229
0.282404
0.000011
0.282401
-13.0
-0.9
1176
1569
-0.99
HC-1-13
553.0
0.007241
0.000308
0.282435
0.000012
0.282432
-11.9
0.1
1135
1501
-0.99
HC-1-14
553.0
0.006267
0.000270
0.282387
0.000012
0.282384
-13.6
-1.5
1200
1607
-0.99
HC-1-15
553.0
0.006233
0.000268
0.282407
0.000012
0.282404
-12.9
-0.8
1172
1562
-0.99
Page 86 of 111
ip t 0.000243
0.282429
HC-1-18
553.0
0.005871
0.000241
0.282399
KDG-4-01
980.0
0.045121
0.001538
0.282369
KDG-4-02
980.0
0.030845
0.001091
0.282290
KDG-4-03
980.0
0.012870
0.000482
0.282338
KDG-4-04
980.0
0.030585
0.001150
KDG-4-05
980.0
0.042430
KDG-4-06
980.0
0.030300
KDG-4-07
980.0
0.024680
0.000915
KDG-4-08
980.0
0.035258
KDG-4-09
980.0
KDG-4-10
0.000015
cr
0.005770
0.282426
us
553.0
-12.1
-0.1
1142
1513
-0.99
0.282396
-13.2
-1.1
1183
1580
-0.99
0.000014
0.282341
-14.3
6.4
1266
1431
-0.95
0.000013
0.282270
-17.1
3.9
1363
1590
-0.97
0.000012
0.282329
-15.4
6.0
1275
1458
-0.99
0.282272
0.000011
0.282251
-17.7
3.2
1389
1631
-0.97
0.001499
0.282335
0.000014
0.282307
-15.5
5.2
1314
1506
-0.95
0.001114
0.282278
0.000011
0.282257
-17.5
3.5
1380
1618
-0.97
0.282246
0.000011
0.282230
-18.6
2.5
1417
1679
-0.97
0.001267
0.282295
0.000012
0.282271
-16.9
4.0
1362
1586
-0.96
0.027611
0.001038
0.282276
0.000013
0.282257
-17.5
3.5
1379
1618
-0.97
980.0
0.033971
0.001216
0.282278
0.000018
0.282255
-17.5
3.4
1384
1622
-0.96
KDG-4-11
980.0
0.051090
0.001764
0.282350
0.000017
0.282317
-14.9
5.6
1301
1483
-0.95
KDG-4-12
980.0
0.025094
0.000901
0.282274
0.000012
0.282257
-17.6
3.5
1378
1618
-0.97
KDG-4-13
980.0
0.035424
0.001379
0.282270
0.000012
0.282245
-17.7
3.0
1401
1645
-0.96
KDG-4-14
980.0
0.010098
0.000429
0.282289
0.000011
0.282281
-17.1
4.3
1340
1564
-0.99
KDG-4-15
980.0
0.029547
0.001075
0.282261
0.000012
0.282242
-18.1
2.9
1402
1653
-0.97
KDG-4-16
980.0
0.030448
0.001091
0.282273
0.000012
0.282253
-17.7
3.3
1386
1628
-0.97
980.0
0.021464
0.000856
0.282307
0.000010
0.282291
-16.5
4.7
1330
1542
-0.97
980.0
0.044247
0.001600
0.282279
0.000014
0.282250
-17.4
3.2
1396
1635
-0.95
980.0
0.043887
0.001672
0.282303
0.000012
0.282272
-16.6
4.0
1365
1584
-0.95
980.0
0.033846
0.001309
0.282305
0.000013
0.282281
-16.5
4.3
1349
1565
-0.96
KDG-4-18 KDG-4-19 KDG-4-20 KDG-4-21
d
ep te
Ac c
KDG-4-17
an
0.000012
M
HC-1-16
980.0
0.028177
0.001111
0.282271
0.000011
0.282250
-17.7
3.2
1390
1633
-0.97
KDG-4-22
980.0
0.033561
0.001188
0.282263
0.000012
0.282241
-18.0
2.9
1404
1653
-0.96
KDG-4-23
980.0
0.024356
0.000897
0.282253
0.000013
0.282236
-18.4
2.7
1407
1665
-0.97
KDG-4-24
980.0
0.075508
0.002571
0.282345
0.000013
0.282297
-15.1
4.9
1338
1529
-0.92
Page 87 of 111
ip t 0.001550
0.282287
KDG-4-26
980.0
0.021552
0.000854
0.282265
KDG-2-01
520.7
0.002531
0.000086
0.281977
KDG-2-02
520.7
0.000244
0.000006
0.282023
KDG-2-03
520.7
0.001201
0.000042
0.282020
KDG-2-04
520.7
0.000201
0.000005
KDG-2-05
520.7
0.003280
KDG-2-06
520.7
0.001328
KDG-2-08
520.7
0.001784
0.000070
KDG-2-09
520.7
0.004295
KDG-2-10
520.7
KDG-2-12
0.000017
cr
0.041104
0.282259
us
980.0
-17.1
3.5
1383
1615
-0.95
0.282249
-17.9
3.2
1389
1636
-0.97
0.000016
0.281976
-28.1
-16.7
1753
2538
-1.00
0.000012
0.282022
-26.5
-15.1
1687
2435
-1.00
0.000014
0.282019
-26.6
-15.2
1693
2442
-1.00
0.281996
0.000011
0.281996
-27.4
-16.0
1723
2493
-1.00
0.000136
0.281896
0.000011
0.281895
-31.0
-19.6
1864
2718
-1.00
0.000048
0.281933
0.000011
0.281933
-29.7
-18.3
1810
2634
-1.00
0.281892
0.000013
0.281892
-31.1
-19.7
1866
2725
-1.00
0.000162
0.281996
0.000015
0.281995
-27.4
-16.1
1730
2497
-1.00
0.000281
0.000007
0.281994
0.000012
0.281994
-27.5
-16.1
1726
2498
-1.00
520.7
0.001944
0.000065
0.281997
0.000013
0.281996
-27.4
-16.0
1725
2493
-1.00
KDG-2-13
520.7
0.000934
0.000030
0.282022
0.000011
0.282022
-26.5
-15.1
1689
2436
-1.00
KDG-2-14
520.7
0.001584
0.000048
0.282058
0.000016
0.282058
-25.2
-13.8
1641
2356
-1.00
KDG-2-15
520.7
0.001525
0.000047
0.282006
0.000009
0.282005
-27.1
-15.7
1712
2473
-1.00
KDG-2-16
520.7
0.002683
0.000108
0.281846
0.000015
0.281845
-32.7
-21.3
1931
2828
-1.00
KDG-2-17
520.7
0.000766
0.000021
0.282039
0.000013
0.282039
-25.9
-14.5
1666
2399
-1.00
KDG-2-18
520.7
0.000173
0.000004
0.282012
0.000011
0.282012
-26.9
-15.4
1701
2458
-1.00
520.7
0.003531
0.000133
0.281983
0.000013
0.281981
-27.9
-16.5
1747
2526
-1.00
520.7
0.000428
0.000011
0.281989
0.000013
0.281989
-27.7
-16.3
1733
2510
-1.00
755.0
0.003763
0.000131
0.282389
0.000012
0.282387
-13.5
3.0
1193
1471
-1.00
755.0
0.002483
0.000095
0.282474
0.000010
0.282473
-10.5
6.1
1075
1279
-1.00
755.0
0.003083
0.000120
0.282408
0.000011
0.282406
-12.9
3.7
1167
1429
-1.00
KDG-1-05
755.0
0.015456
0.000596
0.282344
0.000011
0.282336
-15.1
1.2
1270
1586
-0.98
KDG-1-06
755.0
0.028044
0.001131
0.282370
0.000014
0.282353
-14.2
1.9
1252
1547
-0.97
KDG-1-07
755.0
0.003021
0.000115
0.282507
0.000011
0.282505
-9.4
7.2
1031
1206
-1.00
KDG-2-20 KDG-1-01 KDG-1-03 KDG-1-04
d
ep te
Ac c
KDG-2-19
an
0.000011
M
KDG-4-25
Page 88 of 111
ip t 0.000677
0.282322
KDG-1-09
755.0
0.024695
0.000976
0.282313
KDG-1-10
755.0
0.013257
0.000540
0.282341
KDG-1-11
755.0
0.025974
0.001031
0.282348
KDG-1-12
755.0
0.020893
0.000864
0.282342
KDG-1-13
755.0
0.010112
0.000385
0.282353
KDG-1-14
755.0
0.014504
0.000592
KDG-1-15
755.0
0.014082
0.000570
KDG-1-16
755.0
0.021951
0.000844
KDG-1-17
755.0
0.017256
KDG-1-18
755.0
KDG-1-19
0.000014
cr
0.017363
0.282313
us
755.0
-15.9
0.4
1303
1638
-0.98
0.282300
-16.2
-0.1
1325
1667
-0.97
0.000011
0.282334
-15.2
1.2
1272
1591
-0.98
0.000013
0.282333
-15.0
1.1
1279
1592
-0.97
0.000012
0.282330
-15.2
1.0
1281
1599
-0.97
0.000013
0.282348
-14.8
1.6
1251
1560
-0.99
0.282337
0.000013
0.282329
-15.4
1.0
1279
1602
-0.98
0.282278
0.000013
0.282270
-17.5
-1.1
1360
1733
-0.98
0.282340
0.000012
0.282328
-15.3
0.9
1284
1605
-0.97
0.000737
0.282341
0.000012
0.282330
-15.3
1.0
1279
1599
-0.98
0.018840
0.000774
0.282336
0.000012
0.282325
-15.4
0.9
1287
1610
-0.98
755.0
0.007692
0.000315
0.282320
0.000014
0.282315
-16.0
0.5
1294
1633
-0.99
KDG-1-20
755.0
0.022287
0.000900
0.282337
0.000015
0.282324
-15.4
0.8
1290
1613
-0.97
KDG-1-21
755.0
0.015511
0.000630
0.282291
0.000012
0.282282
-17.0
-0.7
1344
1705
-0.98
KDG-1-22
755.0
0.025541
0.001046
0.282350
0.000014
0.282335
-14.9
1.2
1276
1587
-0.97
KDG-1-23
755.0
0.019186
0.000755
0.282370
0.000014
0.282360
-14.2
2.1
1238
1533
-0.98
KDG-1-24
755.0
0.017021
0.000682
0.282355
0.000012
0.282345
-14.7
1.6
1257
1565
-0.98
KDG-1-25
755.0
0.013953
0.000561
0.282291
0.000013
0.282283
-17.0
-0.6
1342
1704
-0.98
755.0
0.005558
0.000218
0.282343
0.000012
0.282340
-15.2
1.4
1259
1578
-0.99
755.0
0.021774
0.000862
0.282363
0.000012
0.282351
-14.5
1.7
1253
1553
-0.97
KDG-1-27
d
ep te
Ac c
KDG-1-26
an
0.000014
M
KDG-1-08
Page 89 of 111
Ac ce p
te
d
M
an
us
cr
ip t
Figure-1
Page 90 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-2
Page 91 of 111
Ac ce p
te
d
M
an
us
cr
ip t
Figure-3
Page 92 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-4
Page 93 of 111
Ac ce p
te
d
M
an
us
cr
ip t
Figure-5 part I
Page 94 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-5 part 2
Page 95 of 111
Figure-6
(b)
0.26
rim
KDG-1(GD) KDG-2 (GA) KDG-4 (GR) HC-1 (MG) HC-7 (MGB) HC-8 (MGB) HC-11 (MG) WC-2 (GR)
core
0.24
0.22
0.06
0.05
0.04
ip t
(a)
mantle
core
0.03 rim rim
0.02 rim
0.18 0.15
0.2
core
0.25
core
0.01 0.15
0.3
(d)
M
50
augite
Wo
endiopside
d
pigeonite
0.25
an
diopside hedenbergite 50
0.2
core
0.3
Mg/(Fe+Mg)
Mg/(Fe+Mg)
(c)
rim
rim
us
core
cr
0.2
mantle
50
50
Or
orthopyroxene 50
(e)
Ac ce pt e
En
Fs
(f)
0.65
50
Ab
An
6
CaB > 1.5, (Na+K)A > 0.5
5.5
0.6
0.55
0.5
0.45
0.4
5
pargasite
edenite
4.5
kaersutite
ferro-edenite
6.8
6.6
4
ferro-pargasite 6.4
Si (pfu)
6.2
6
3.5
0.4
0.45
0.5
0.55
0.6
0.65
Mg/(Fe+Mg)
Page 96 of 111
Figure-7
Al2O3
16
2
14
1
40
25
45
50
55
60
65
70
(c)
75
40
15
45
50
55
60
65
70
(d)
0.3
Fe2O3
20
0
MnO
0.2
10
10
45
50
55
60
65
70
75
0
40
14
(e)
MgO
8
45
50
8
40
45
50
55
60
65
75
2
40
45
5
(g)
4
70
55
60
(h)
Na2O
4
Ac ce pt e
3
3
50
70
75
65
70
75
70
75
70
75
K2O
d
0
4
M
6
2
65
CaO
10
4
60
(f)
12
6
55
an
40
us
0.1 5 0
75
cr
12
TiO2
(b)
3
ip t
(a)
18
2
2
1
1
40
2
45
(i)
50
55
Metagranodiorite
1.5
Metagabbro
HC-7 (ferroan)
KDG-4 (felsic)
HC-8 (magnesian)
70
75
P2O5
0
40
45
(j)
4
Mafic granulite
50
55
Grt amphibolite
Fe2O3/MgO
Ferroan
Felsic
Grt charnockite
2
HC-10 (magnesian)
Magnesian
1
40
45
50
65
HC-11 (magnesian)
KDG-2 (ferroan)
0
60
3
HC-1 (magnesian)
KDG-1 (magnesian)
0.5
65
WC-2 (felsic)
Metadiorite
1
60
55
60
SiO2 (wt.%)
65
70
75
40
45
50
55
60
65
SiO2 (wt.%)
Page 97 of 111
Figure-8
10 2
10 3 (d) Magnesian
10 2
1
10
10
Sr
K
Rb Ba Th Ta Nb Ce
P
Zr
Hf Sm Ti
Y
Yb
2
1
Ce
Pr
Nd Sm Eu Gd Tb
(e) Ferroan
10 2
Metagabbro
1
Grt amphibolite
Rb Ba Th Ta Nb Ce
P
Zr
Hf Sm Ti
Y
1
M
KDG-2 K
Er Tm Yb
Lu
10
HC-7
Sr
Dy Ho
an
10
10
La
10 3 (b) Ferroan
0.1
HC-8
us
0.1
Grt charnockite
HC-1 HC-10 HC-11 Metadiorite KDG-1 Metagabbro
ip t
10
Mafic granulite
cr
(a) Magnesian
Yb
3
La
Ce
Pr
Nd Sm Eu Gd Tb
Dy Ho
Er Tm Yb
Lu
Pr
Nd Sm Eu Gd Tb
Dy Ho
Er Tm Yb
Lu
10 3
(c) Felsic
(f) Felsic
Metagranodiorite
10
1
Sr
K
Ac ce pt e
10
d
WC-2 KDG-4
2
Rb Ba Th Ta Nb Ce
P
Zr
Hf Sm Ti
Y
Yb
10 2
10
1
La
Ce
Page 98 of 111
Figure-9
100
1000
Metagranodiorite WC-2 KDG-4
WPG
VAG + syn-COLG
syn-COLG 1
10
ORG
1
10
Y (ppm)
100
1000
0.1
0.1
1
1000
1000
syn-COLG WPG
10
10
100
an 1
1
10
100
Yb + Ta (ppm)
d
Y + Nb (ppm)
1000
ORG
Ac ce pt e
1
VAG
M
ORG
100
us
VAG
10
10
WPG
100
100
Yb (ppm)
cr
syn-COLG
1
ORG
VAG
ip t
1
WPG
Metadiorite KDG-1
10
100
Page 99 of 111
TiO2
ip t
Figure-10
(a)
cr
2Nb
MORB
M an
AI OIT
3D 1C
1B
2B
AII
3C
3B
B
IAT OIA
1A
C
Bon CAB
2A
3A
Y
La/10
Nb/8
ce pt
P2O5x10 Zr/4
ed
D
Ac
MnOx10
Y/15
(c)
us
(b)
Mafic granulite HC-1 (magnesian) HC-11 (magnesian) Metagabbro HC-7 (ferroan) HC-8 (magnesian) Grt amphibolite KDG-2 (ferroan) Grt charnockite HC-10 (magnesian)
Page 100 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-11
Page 101 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-12
Page 102 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-13
Page 103 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-14
Page 104 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-15
Page 105 of 111
Ac ce p
te
d
M
an
us
cr
ip t
Figure-16
Page 106 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-17
Page 107 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-18
Page 108 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-19
Page 109 of 111
Ac
ce
pt
ed
M
an
us
cr
i
Figure-20
Page 110 of 111
Ac ce p
te
d
M
an
us
cr
ip t
Figure-21 revised
Page 111 of 111