In this study, zircon grains are applied for U–Pb dating, Hf isotopes and trace elements to reveal the origin of magmatism and tectonic evolution of Late Paleozoic rocks of the Indian plate, Northern Pakistan. Most of the zircons are characterized by oscillatory zoning, depletion of light rare earth elements (LREE) and enrichment of heavy rare earth elements (HREE) with Ce and Eu anomalies. The yielded ages for these rocks are 256 ± 1.9 Ma and are plotted in the zones defined for the continental setting with few deviated toward the mid-oceanic ridge and the oceanic arc setting. Deviated zircons are recognized as inherited zircons by displaying a high concentration of normalized primitive La and Pr values, while others are plotted in the continental zones. Rare earth elements (REE) and trace elements including Th, Hf, U, Nb, Sc and Ti discriminate Swat orthogneisses into the within plate setting and the inherited zircons are plotted in the orogenic or the arc-related setting. The LREE discriminated these zircons into a magmatic zone with inherited zircons deviated toward the hydrothermal zone. The temperature calculated for these rocks based on the Ti content in zircon ranges from 679 to 942°C. The εHf(t) ranging from −11.1 to +1.4 reveals that the origin is the continental crust with the minute input of the juvenile mantle.
Traditionally, the importance of zircon (ZrSiO4) is known for its geochronological dating, but zircon trace element geochemistry is an emerging field to interpret tectonic settings [1,2,3]. Titanium present in the zircon can help calculate the temperature of the rocks, and the Hf isotopic study can play a significant role to unravel crustal recycling [4,5,6,7]. LREE depletion and HREE enrichment with Eu and Ce anomalies are used to know the magmatic and metamorphic origins of the rocks from the zircon study [8,9]. Zircon occurs in several igneous, sedimentary and metamorphic rocks as an accessory mineral . It is due to its high durability and stability to harsh physical and chemical weathering as well as due to its dominant role in the growth with granitoid petrology that it is regarded as the most common mineral used in igneous, metamorphic and sedimentary petrology [10,11]. Trace element geochemistry in zircon is an emerging field to interpret the origin and the tectonic setting of igneous rocks [2,12,13], and this technique is applied for the first time on the Swat orthogneisses, Northern Pakistan. The REEs like Ce and HREE present in zircons display unique enrichment normalized trends, which are used to find the protolith of the rock [3,14]. As hafnium is most compatible with zircon, the ages determined by the U–Pb dating are also used for the Lu–Hf isotope system, which is the best choice for determining the evolution of igneous and metamorphic rocks in collisional zones . Lu is least compatible with the zircon, while hafnium and zirconium are geochemically similar, and 1–4% of hafnium is present in the zircon, which is geotectonically very important. The zircon is a very resistant mineral to harsh chemical and physical weathering, and during physical weathering, the zircon moves to the sand portion along with hafnium, while the Lu is discarded into the clay portion, which is the basis for finding the geochemical settings [5,16]. Although the importance of Hf as a geochemical tracer was recognized in the early 1980s, Hf separation was very difficult, and later the MC-ICP-MS technique was adopted to determine the accurate ratios of Hf isotopes in the zircon, which is a revolution to study the mantle crust processes in detail as a geochemical tracer in subduction settings [17,18]. The Hf isotopic record in the zircons is the best decipher for source rocks as well as plays important roles in the magma mixing, assimilation and petrogenesis [9,19,20,21,22]. It is noted that the geochemical behavior of Zr and Hf during magma crystallization is the same, and all the zircons have 0.4–4 wt% Hf, which is a part of the zircon lattice . The Lu is not compatible with zircon, and hence, the initial Lu/Hf ratio is much less (<0.001) compared to the continental crust with the high Lu/Hf ratio. The crust/mantle differentiation in the context of 176Hf/177Hf is different, and this specific characteristic of Hf is used as the trace for differentiation [5,23]. The technique to determine the Hf isotopic ratio from the same spot that is used for finding the age is a new development in the petrochronology [24,25,26]. The previous ages determined for Swat orthogneiss were 268 ± 7 Ma, 265 ± 2, 278 ± 4 by U–Pb and Rb/Sr techniques [27,28]. Crustal and mantle processes are particularly studied in detail by using the evolution of hafnium in zircon . The morphology and the chemistry of trace elements of the zircon combined with the Hf isotope composition help interpret petrogenetic information . In this article, Swat orthogneisses (Figure 1) of zircons are analyzed for Hf isotopes, U–Pb dating and trace elements to find the origin and the evolution of magma. The zircon morphology, trace elements in zircon and Hf isotopes from the zircon were determined by using cathodoluminescence (CL) imaging, laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) and laser ablation multicollector inductively coupled plasma mass spectrometry (LA-MC-ICPMS).
2 Geological setting
Northern Pakistan is tectonically divided into three terranes from the south to the north, i.e., Indian plate rocks in the south, Kohistan–Ladakh Island arc in the middle and Karakoram plate rocks in the north [30,31]. The main mantle thrust (MMT) separates the Indian plate rocks (Figure 1) from the Kohistan island arc in the north [32,33,34]. The Kohistan Island in the north is bounded by the main Karakoram thrust or the northern suture zone. The part of the Indian plate in Pakistan is further divided into three tectonic bodies by the main central thrust (MCT), main boundary thrust and salt range thrust. The zone between the MCT and the MMT (Figure 1) is a part of greater Himalaya [35,36,37,38]. The Indian plate rocks in the south are metamorphosed as a result of initial subduction phases and later on collision with the Kohistan-Ladakh island arc in Cretaceous to Cenozoic Himalayan orogeny . The Greater Himalayan zone consists of pelitic-psammitic schists, gneisses, amphibolites, meta-carbonates and intrusion of plutonic bodies ranging in age from Archean to recent (ca. 1,847–47 Ma) [28,39,40]. The Himalayan orogeny has great effects in the context of sedimentation, metamorphism and plutonism on the Indian plate rocks, which are formed as early as the Archean era. Granitic rocks are the most abundant intrusive bodies ranging in size from plutons to smaller size sills and dykes in northern Pakistan . The granitic rocks from northern Pakistan were grouped based on their radiometric ages, but the data available for these extensive rock bodies are limited [35,42]. Swat orthogneisses were metamorphosed to greenschist and amphibolite facies during Eocene–Oligocene [43,44]. These orthogneisses are divided into seven groups starting from Early Proterozoic (ca. >2,175 Ma) to Late Quaternary (ca. 10 Ma). This study focuses on Late Paleozoic (ca. 350–268 Ma) orthogneisses by interpreting the trace elements from zircon grains. Swat orthogneisses were explored by using whole-rock major and trace element geochemistry, as well as geochronological studies were reported from this terrane, but zircon trace element geochemistry is limited [3,5,14]. The felsic plutonic activity is the main cause of crustal differentiation, but to correlate felsic magma with the crust development, knowledge of magmatic evolution is required, which is difficult from the whole rock composition. Hf isotopes extracted from the zircons are the best tracers for finding the magma evolution in continental granitoids, continental arcs and subduction zones [19,45]. Hf isotopes in the zircons are selected as it is highly compatible with the zircon. The precise ages determined by the U–Pb system in zircon can accurately provide the isotopic evolution of Hf. Both Hf and Lu are incompatible elements, but Hf segregated more strongly into the melts than Lu during magma differentiation stages. The concentration of 176Hf/177Hf is increased in the mantle compared to the crust. Granitic magmas with high 176Hf/177Hf are the indication of magma intruded directly from mantle or by melting of the lower crust. The lower concentration of 176Hf/177Hf is determined for magma from crustal reworking [19,45].
3 Analytical techniques
3.1 Zircon U–Pb dating and Hf isotope
Zircons grains were extracted from coarsely crushed orthogneisses (80 mesh). These zircons were separated by using the combined technique of magnetic separation and heavy liquid before handpicking under a binocular microscope. These separated zircons were mounted in epoxy resins and were polished to half size to best expose grain interiors. The CL images and transmitted and reflected images of the polished zircons were obtained (Figure 2a–f). Detailed interpretation of all these images was carried out to select those areas that have no cracks and inclusions for laser ablation. The zircon U–Pb dating was performed at the state key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Beijing (GPMR-CUGB). Zircon standard 91500 was used as an external standard as described in Ref. . The standard was run two times after every eight samples. The first 20 s were used for blank gas followed by 50 s for sample. The laser ablation was performed by using a GerlasPro laser ablation system comprising a COMPexPro 102 ArF excimer laser with a wavelength of 193 nm and a maximum energy of 200 mJ and a MicroLas optical system. The laser ablation was done with a ablation spot size of 35µm, and ion signal intensities were taken by the ICPMS instrument (Agilent 7700×). Helium gas was applied as a carrier gas, and Argon was used as a make-up gas that mixes with the carrier gas through a T-connector prior entering the ICP. The standards of the U–Pb dating analysis are Plešovice, 91500, GJ-1 and SRM 610, which were used as external standards . The concentrations of U–Pb and trace elements obtained from the LA-ICPMS analysis were interpreted for results on the ICPMS-Data-Cal 10.8 software by using the method defined in Ref. . The common lead correction was carried out by following the method designed in Ref. . U–Pb Concordia and the weighted average mean were obtained through Iso-plot version 2003 by using the method described in Ref. . The REE and trace elements obtained during this research work are listed in Table 1.
|Trace elements (ppm)|
The in situ Hf isotopic analysis of zircon was carried out by using LA-MC-ICPMS at GPMR-CUGB. The Lu–Hf zircon isotopic analysis was performed by using the 35 µm spot size. The geo-standards follow those given in Ref. , including 91500 (176Hf/177Hf value 0.282306 ± 30) and Plešovice (176Hf/177Hf value 0.282477 ± 30). The isotopic concentrations of εHf(i) and εHf(t) were determined by using the chondrite uniform reservoir standards defined in Ref. [19,51,52]. The concentrations of these are as follows: (176Lu/177Hf)CHUR = 0.0336, (176Hf/177Hf)CHUR = 0.282785 and (176Lu/177Hf)mean crust = 0.015. Initially, the zircons were run for U–Pb dating, and the best concordance results were then interpreted for the Hf isotopic study. The calculation of the initial and final values of the Hf isotopic composition followed the method defined in Ref.  and is presented in Table 2.
|Zircon spots||Age (Ma)||176Lu/177Hf||176Hf/177Hf||176Hf/177Hf i||εHf(0)||εHf(t)||TDM||TDM(Hf)C||fLu/Hf|
4.1 Rock petrography
The Swat orthogneisses are essentially composed of potash feldspar, quartz and muscovite with accessory minerals including epidote, apatite, plagioclase, clinopyroxenes and rutile (Figure 3a–d). The modal percentage of potash feldspar including both orthoclase and microcline ranges from 45% to 65% in different views. The microcline displays the coarsed grain texture with anhedral to subhedral shape. The grain size of the microcline ranges from 1 to 2.5 mm with hatched twinning. The microcline shares its boundaries with the muscovite and recrystallized quartz. The quartz is a second dominant phase, and it ranges in modal composition from 15% to 25%. Quartz is mostly anhedral, and the smaller grains display strains that are easily visible from the undulose extinction. The orthoclase displays simple twins and mostly represented by a cloudy surface in the plane light due to sericitization. Muscovite is less than 10% in these rocks and is characterized by the tabular appearance. The dominancy of muscovite over biotite in these rocks represents its peraluminous nature. Biotite is less abundant in most views, and its percentage is <5% but mostly associated with muscovite grains. The laths of tabular muscovite exhibit truncated contact with microcline. Megacrysts of potash feldspar (microcline) are the common mineral in these rocks. Inclusions of biotite grains are present in the potash feldspars. The graphic texture is observed where quartz shares a boundary with a microcline. Quartz grains with the anhedral to subhedral shape are clear, while the alkali feldspars have minute input of sericitization and saussuritization in the form of sericite, muscovite and epidote. Biotite and muscovite display subhedral to euhedral shape, and topotaxial growth is also observed in small grains where muscovite is crystallized at the expense of biotite. Clusters of epidotes occur in these rocks with dominant zoning, which are bluish and orange that represent metamorphic stages in these rocks. Biotite is fresh looking with perfect one set of cleavage with some spots altered to greenish color chlorite in the center. Recrystallized quartz grains mostly of the anhedral shape are present along the boundary of microcline, and these are developed during high-grade metamorphism. The rock is characterized by hypidiomorphic texture with most of the minerals represented by anhedral with few euhedral shapes. The biotites are less than 0.5 mm, and muscovite is coarse grained with an average size of 0.5–1 mm. The quartz grains are less than 0.5 mm to more than 2 mm in different views. The recrystallized quartz grains are represented by smaller anhedral shapes, and the grain boundaries are sutured. Traces of rutile as an euhedral mineral are also observed in these rocks. The regional stress direction is represented by the orientation of muscovite, biotite and quartz grains in these rocks.
4.2 Zircon morphology and age
The zircons extracted were mostly 100–150 µm in length and exhibit cracks as examined in the CL images. Those zones that display igneous zoning in CL images and display no inclusions in transmitted light were selected, and the care was also taken to select those areas, which were not diluted by Pb loss due to internal stresses and fractures. These zircons have preserved the oscillatory zoning, which is the characteristic of igneous zircons, and the broad zoning without visible cores reveals that these zircons were developed during one stage of the magmatic growth . Prismatic grains in these zircons display weak oscillatory zoning, while nonprismatic subhedral elliptical grains are characterized by good oscillatory zones. There are lighter zones that display their metamorphic history, while the magmatic sector and oscillatory zones are preserved in almost all zircons (Figure 2). Most of them are subhedral to euhedral, and due to the long tectonic history, some are split into pieces as a result of tectonic forces. The uranium content in these zircons ranges from 74 to 3,192 ppm, while the Th/U content of the analyzed zircons ranges from 0.087 to 1.15, which display its magmatic age . The concordant age (Figure 4a and Supplementary Table 3) determined for this sample is 256 ± 1.9 Ma (mean square weighted deviation (MSWD) = 2.7, n = 15). The age is calculated based on 15 zircons in which the concordance of most of the grains is more than 95%. The range of time calculated for the granitic rocks from this study is 249–256 Ma based on the 206Pb/238U ratio. The Lu(n) values in these zircons range from 2,940 to 13,010, similar to the range defined in Ref.  for the igneous-type granitoids. The Y concentration in these zircons is 937–4,220 ppm, which is similar to granitoids (500–4,534 ppm) reported in Ref. .
4.3 Zircon trace elements and Hf isotopes
Zircons were analyzed for REE and trace elements including Sc, Ti, Y, Zr, Nb, Hf and Ta to find out the tectonic setting of Swat orthogneisses (Table 1). The concentrations of trace elements are as follows: Sc, 274 to 547 ppm; Ti, 4.6–57.4 ppm; Y, 937–4,220 ppm; Nb, 1.6–34.4 ppm; Hf, 10,869–17,994 ppm; and Ta, 1.8–14.5 ppm. All the zircons were characterized by negative Eu and positive Ce anomaly (Figure 4b) as plotted on the normalized primitive values in Ref. . The concentration of REE in these zircons was dominated by the enrichment of HREE and depletion of LREE. The Ti content was variable and converted to a temperature range of 679–942°C. The highest temperature of >800°C was displayed by the metamorphic region that has a high concentration of Ti, and the zircons having a temperature range of <750°C are all the magmatic zircons with a low titanium content . The composition of Hf in the context of the 176Hf/177Hf ratio ranges from 0.282315 to 0.282684, and εHf(t) for this composition is −11.15 to +1.47.
5.1 Tectonic setting of Swat orthogneisses
The origin of granitic rocks based on multiple trace elements like Hf, U, Th, Y and HREE easily distinguished the island arc setting, within plate orogeny and convergent plate orogeny [3,14]. Trace elements in zircons, i.e., U, Y and Hf, are used to distinguish the continental and oceanic setting . The REE normalized pattern of zircon in this study is shown in Figure 4b. Zircons from Swat orthogneisses are characterized by continental granitoids as plotted in Ref.  the normalized pattern of REE. The LREE depletion with prominent positive Ce anomaly and the enrichment of HREE discriminated these zircons into typical magmatic origin . The pattern is compared with the standard given in Ref.  for zircons, and it is found that the zircons are plotted on the continental zone with less or no contamination of the oceanic island setting. Few of the zircons are noted with high La (>5 ppm) and Pr (>10 ppm) from the LREE, and these are included in the inherited zircons as the enrichment of La and Pr is a dominant feature in the inherited zircons . Granitic rocks are further classified into I-type, S-type and A-type on the basis of REE patterns, and in this study, zircons show the I-type concentration of REE . The Eu-negative anomaly in the REE pattern normalized to  primitive values is the characteristic feature of plagioclase crystallization in these rocks. The Th/U and Nb/Hf  ratios are critical to distinguish orogenic within the plate setting from the arc-related orogenic setting (Figure 5a). The zircons from the Swat orthogneisses have a high concentration of the Nb/Hf ratio, which is a clue for these zircons to have originated in the anorogenic stable within the plate setting as a result of plutonic magma in the continental crust. The inherited zircons with a low concentration of Nb/Hf are plotted in the zone for arc or orogenic setting. The high concentration of Th than Nb and similarly the high concentration of Hf than Nb are characteristic features of zircon from the continental arc-related orogenic setting. The zircons in this study are plotted in the anorogenic within the plate zone  with the low concentration of Hf and Nb compared to the arc-related setting, while the inherited zircons are deviated toward the zone or the arc setting (Figure 5b). The Hf content in the zircons in this study is higher (>10,000 ppm), which reveals that the origin of magma for these rocks is from more evolved felsic magma . The heavy REE element Yb in the zircon has a characteristic to be differentiated into the oceanic zone and the continental zone correlated with the uranium content . The zircons in this study have a high concentration of uranium and are plotted in the zone defined for the continental survey. The rocks are plotted in the zone of anorogenic within the plate setting in the figures of Refs. [61,62] on the basis of Th/Nb versus Hf/Th and Th/U versus Nb/Hf ratios (Figure 5a and b).
It is noted that multiple trace element ratios are more efficient to discriminate the tectonic setting compared to the single element. The ratio of trace elements like Nb/Yb versus U/Yb (Figure 5c) is helpful to define the tectonic setting for Swat orthogneisses to have originated in the continental arc-type setting, and the similar setting is determined for these rocks as plotted on the binary (Figure 5d) Sc versus U . U/Yb is used for distinguishing the zircons from the continental and oceanic settings as major changes are observed in the ratio for these tectonic settings [64,65]. The average U/Yb values for N-MORB, primitive arc basalts, average lower and upper crust are ∼0.01, 0.3, 0.13 and 1.35, respectively. Based on these values, the tectonic settings are easily distinguished. The Swat orthogneisses have the U/Yb ratio ranging from 0.18 to 3.47 and the high values indicating the enriched crust, which have influence on the mature arc. The higher ratio (>0.1) discarded the phenomenon of these rocks to have originated from the mid-ocean ridge basalt (MORB) setting . Zircon is characterized by partitioning Yb compared to U, and thus, the concentration of U in the magma increased in the later stages of magma is increased . This increasing trend is used to differentiate the primary magma from the MORB and the later secondary stages of magmas.
Niobium  is the tracer for interpreting the subduction zones by its normalized negative anomaly, but in the zircons trace element geochemistry, the enrichment of Nb is plotted for felsic igneous within plate rocks and alkali-rich magmas. The zircons are plotted in the field of the continental arc setting based on Nb/Yb and U/Yb discrimination (Figure 5c), and the crystallization of these granitic rocks in the continental setting is confirmed. Trace elements from the zircons are used to correlate with rock types, and the concentration is documented for classifying the igneous rocks . The dominant trace elements used for discriminating igneous rocks into different fields include Y, U, Yb, Sm, Nb and Ta, and the number of rocks recognized based on these elements from zircons ranges from ultramafic to felsic compositions, i.e., kimberlites to leucogranites. The Yb/Sm versus Y composition of zircons from Swat orthogneisses is used to know the rock type on the plot , and the zircons are plotted in the field of granodiorites (Figure 6). The Ta and Nb as well as Ce/Ce* versus Y plot discriminate the zircons into the same granodiorites. Similarly, Nb/Ta versus Y plot also discriminates the zircons from Swat orthogneisses into the granodiorite field. It is confirmed from the trace element concentration that the original rock in which these zircons grow is granodiorite as shown in all the discriminatory figures. The granitoid field in the previous study  is further divided into three zones, namely, (1) aplites and leucogranites, (2) granites and (3) granodiorites and tonalites. The rocks are plotted in the granodiorites, which is subfield 3 from the discrimination diagram . The rock classification figures of the trace elements are shown in Figure 6. The ratio of Nb/Ta analyzed in this study ranges from 0.9 to 2.73 with an average value of 2 and the chondritic Nb/Ta ratio  is 17. The high ratio is found in most fractionated rocks like granites, while the low ratio is reported in kimberlites. The Nb/Ta and Y are positively correlated as the ionic size of Nb is larger than Ta, and thus, Ta is most compatible with zircon during magma differentiation. The normalized pattern of trace elements is used to distinguish between the magmatic and the hydrothermal zircons . The LREE (La and Sm) normalized values are plotted versus La , and most of the zircons in this study are plotted in the magmatic zone with few deviated toward the hydrothermal zone (Figure 7a). The deviated zircons also have high contents of La and Pr on the normalized spider pattern (Figure 4b) and are included in the inherited zircons. The deviation of magmatic zircons toward the hydrothermal zone is also documented in Ref.  as the inclusion of apatite, monazite and titanite, which can greatly enhance the composition of LREE and thus deviates from the magmatic zone to the hydrothermal zone. The hafnium concentration of this study zircons plotted versus U/Yb has the characteristic of the enriched mantle zone or crustal input in the continental zone (Figure 7b). The Hf concentration is >10,000 ppm, and the ratio of U/Yb is more than 0.1. The Yb plotted versus Ti (Figure 7c) shows the input of the mid-oceanic ridge setting with most of the zircons plotted in the continental zone compilation with the arrow marking toward the amphibole absent zone. Most of the zircons have Ti content less than 10 ppm plotted in the continental compilation with few having a Ti content greater than 10 ppm in the mid-oceanic ridge compilation .
The boundary between the continental and oceanic crust setting is best explained by U/Yb and Th/Yb ratios from zircons [3,14]. The average U/Yb ratios for N-MORB, primitive basalts, lower and upper crust are 0.01, 0.3, 0.13 and 1.35, respectively [64,65]. This ratio calculated for the Swat orthogneisses ranges from 0.18 to 3.47 with an average value of 1.72, which is for the upper crustal setting. Few zircons are recognized as hydrothermal based on a flattened normalized light REE pattern . The Sc/Yb ratio from the zircon is used as a proxy for tectonic settings including continental arc, mid-ocean ridge and ocean island arc. The high silica crustal island arcs have a high concentration of Sc/Yb (>0.1), and the zircons extracted from the Swat orthogneisses have the same concentration, i.e., 2.6–7.4 with an average value of 5.2. The mid-ocean ridge and oceanic arcs have less than 0.1 concentration reported for Sc/Yb [3,14]. LREE enrichment is commonly related to hydrothermal alteration, and Swat orthogneisses also have this pattern for few zircons as shown in Figure 4b. The roles of inclusions like apatite, titanite and allanite can also enhance the LREE concentration [69,70]. Th enrichment pattern of LREE is also associated with sub-solidus condition during metamorphism, and Swat orthogneisses displayed the normalized enrichment pattern of few zircons due to the regional metamorphism . The Swat orthogneisses are metamorphosed from granitic protolith, which are plotted in the anorogenic plot in Ref. , but few zircons deviated toward the arc-related orogenic setting. The deviated zircons that display arc-related settings are characterized by the high concentration of LREE.
5.2 Temperature of Swat orthogneisses
Temperature calculation of the Ti content present in the zircon is a technique to evaluate the minimum crystallization temperature of the rock [6,7,72,73]. For high-grade metamorphic rocks, the techniques defined in Refs. [74,75] are used to determine the temperature of Ti in the zircon. The titanium content in the zircons from the Swat orthogneisses ranges from 4.6 to 57.4 ppm, and there is a temperature variation from 679 to 942°C. The temperatures calculated for the lowest and highest values in Ref.  are 678.8 and 942.4°C. The variation in temperature calculated for these zircons is due to the magmatic origin and metamorphism, and the inherited zircons also have a high concentration of Ti, which is converted to high temperature. Igneous zircons have Ti concentration in the range of 4–10 ppm, and displayed temperature ranges from 700 to 800°C. These zircons display a temperature range above granitic rocks solidus and thus are characterized by the actual temperature of the zircon crystallization . These zircons are also plotted in the continental zones on different variation figures as discussed in the discrimination portion of trace element. Inherited zircons have Ti content ranging from 20 to 30 ppm, and they have a temperature range of >800°C. The inherited zircons with the variable Hf content and the high Ti content display very high temperature, and these values are not related to the bulk rock crystallization conditions, but either they later on mixed with the magma or the variation is also due to inclusion of Ti bearing minerals in the zircon like ilmenite and rutile. The Hf versus the increasing temperature trend is shown in Figure 7d.
5.3 Hafnium isotopic evolution
The Lu–Hf system in zircon is the best choice for interpreting the geological setting as a geochemical tracer. Hf isotope from the zircon is used to know the magma mixing and the composition of magma in the crystallization of granitoids. The Lu–Hf system is also characterized by tracking back the chemical composition of crust and mantle as Lu is fractionated from Hf during crust mantle differentiation [61,62]. The half-life of this system is approximately 35.7 billion years, and the initial 176Hf/177Hf ratio in the zircon is preserved due to the high incompatibility of Hf than Lu during mantle melting [18,77]. The εHf value for the rocks derived from mantle obtained in Ref.  is +14, while the εHf value of the current oceanic basaltic magma is +23. These values show high variation, and the magma sources for these mantle rocks are characterized into different sources . The negative values of εHf (−10 and −12) are recognized for the lower crust by the previous pioneer works [78,79]. Another source for the enriched εHf is the subduction of the oceanic crust beneath the continental or oceanic plate. The variation in Hf isotopic composition is highly variable, and it is noted from the results that there is an input of the juvenile crust with time in the Swat orthogneisses. Two sources are concluded from the isotopic composition, i.e., initially, the crust derived magma is later on mixed with the juvenile crustal input material with the positive Hf values from negative ones. The zircons from this study have dominantly negative values with some partially positive input, i.e., εHf(t) −11.5 to +1.47. The high variation in the values of εHf(t) is due to inherent heterogeneity in the source magma from which these rocks are crystallized. The age versus εHf(t) plot in Figure 8 shows that most of the zircons are clustered at εHf(t) = −4 to +1 and shows the igneous evolution of magma, which is dominant from the crustal input toward the juvenile mantle. The age versus 176Hf/177Hf ratio displays that these zircons are formed below the depleted mantle zone within the zone of the lower crust in the ranges of >0.282 and <0.283. The inherited zircons are deviated in the εHf(t) plot and deviated from the main magmatic cluster showing highly variable values. The εHf(t) values of the inherited zircons show that these are from the upper crust of the juvenile mantle. The Hf evolution reveals that the Swat orthogneisses rocks have negative values, which is confirmed for its origin from the crustal section, but partial addition of positive values is either from the magma mixing or exposure of these rocks to continental collision. The Hf composition from the zircon spots revealed that the recycled crustal material played a role to generate plutonic magma in this terrane, as well as with time, the Hf composition changes to positive values, which is a clue for mixing of juvenile mantle material in the already generated crustal section [23,77,80].
On the basis of trace elements in zircon, Hf isotopic analysis and geochronology, the following conclusions are presented:
The age (256 ± 1.9 Ma) determined for these rocks in the current research work is younger from the previous calculated ages (i.e., 268 Ma).
The Swat orthogneisses have akin their geochemical behavior similar to upper crustal level rocks based on concentration determined by U/Yb and Sc/Yb ratios.
These rocks are formed in the anorogenic/within plate setting based on trace element concentration reported in this study on zircons.
The average temperature calculated for these rocks ranges from 700 to 800°C based on the Ti content present in the zircons.
εHf(t) has mainly negative values, which confirmed the fact that these rocks originated in the continental setting, and slightly positive εHf(t) addition is either due to magma mixing or from the juvenile mantle zone.
The authors are grateful to Di-Cheng Zhu, Xin Tong and Dong Liu for their help in the revised manuscript. This study was funded by grant from the National Key Research and Development Program of China (2016YFC0600304). The authors are grateful for careful reviews by anonymous reviewers, which greatly helped us improve the paper.
 Watson EB, Chemiak D, Hanchar JM, Harrison T, Wark D. The incorporation of Pb into zircon. Chem Geol. 1997;141:19–31. Search in Google Scholar
 Belousova E, Griffin WL, O’Reilly SY, Fisher N. Igneous zircon: trace element composition as an indicator of source rock type. Contribut Mineral Petrol. 2002;143:602–22. Search in Google Scholar
 Grimes C, Wooden J, Cheadle M, John B. “Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contribut Mineral Petrol. 2015;170:46. Search in Google Scholar
 Amelin Y, Lee D-C, Halliday AN, Pidgeon RT. Nature of the Earth’s earliest crust from hafnium isotopes in single detrital zircons. Nature. 1999;399:252. Search in Google Scholar
 Hoskin PW, Schaltegger U. The composition of zircon and igneous and metamorphic petrogenesis. Rev Mineral Geochem. 2003;53:27–62. Search in Google Scholar
 Watson E, Wark D, Thomas J. Crystallization thermometers for zircon and rutile. Contribut Mineral Petrol. 2006;151:413. Search in Google Scholar
 Ferry J, Watson E. New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers. Contribut Mineral Petrol. 2007;154:429–37. Search in Google Scholar
 Trail D, Watson EB, Tailby ND. Ce and Eu anomalies in zircon as proxies for the oxidation state of magmas. Geochim Cosmochim Acta. 2012;97:70–87. Search in Google Scholar
 Kemp A, Hawkesworth C, Foster G, Paterson B, Woodhead J, Hergt J, et al. Magmatic and crustal differentiation history of granitic rocks from Hf–O isotopes in zircon. Science. 2007;315:980–3. Search in Google Scholar
 Corfu F, Hanchar JM, Hoskin PW, Kinny P. Atlas of zircon textures. Rev Mineral Geochem. 2003;53:469–500. Search in Google Scholar
 Belousova E, Griffin W, O’Reilly SY. Zircon crystal morphology, trace element signatures and Hf isotope composition as a tool for petrogenetic modelling: examples from Eastern Australian granitoids. J Petrology. 2005;47:329–53. Search in Google Scholar
 Chapman JB, Gehrels GE, Ducea MN, Giesler N, Pullen A. A new method for estimating parent rock trace element concentrations from zircon. Chem Geol. 2016;439:59–70. Search in Google Scholar
 Lei W, Shi G, Santosh M, Ng Y, Liu Y, Wang J, et al. Trace element features of hydrothermal and inherited igneous zircon grains in mantle wedge environment: a case study from the Myanmar jadeitite. Lithos. 2016;266:16–27. Search in Google Scholar
 Grimes CB, John BE, Kelemen P, Mazdab F, Wooden J, Cheadle MJ, et al. Trace element chemistry of zircons from oceanic crust: a method for distinguishing detrital zircon provenance. Geology. 2007;35:643–6. Search in Google Scholar
 Tang M, Wang X-L, Shu X-J, Wang D, Yang T, Gopon P. Hafnium isotopic heterogeneity in zircons from granitic rocks: Geochemical evaluation and modeling of “zircon effect” in crustal anatexis. Earth Planet Sci Lett. 2014;389:188–99. Search in Google Scholar
 Patchett PJ. Importance of the Lu–Hf isotopic system in studies of planetary chronology and chemical evolution. Geochim Cosmochim Acta. 1983;47:81–91. Search in Google Scholar
 Yuan H-L, Gao S, Dai M-N, Zong C-L, Günther D, Fontaine GH, et al. Simultaneous determinations of U–Pb age, Hf isotopes and trace element compositions of zircon by excimer laser-ablation quadrupole and multiple-collector ICP-MS. Chem Geol. 2008;247:100–18. Search in Google Scholar
 Gerdes A, Zeh A. Zircon formation versus zircon alteration – new insights from combined U–Pb and Lu–Hf in situ LA-ICP-MS analyses, and consequences for the interpretation of Archean zircon from the Central Zone of the Limpopo Belt. Chem Geol. 2009;261:230–43. Search in Google Scholar
 Griffin W, Wang X, Jackson S, Pearson N, O’Reilly SY, Xu X, et al. Zircon chemistry and magma mixing, SE China: in situ analysis of Hf isotopes, Tonglu and Pingtan igneous complexes. Lithos. 2002;61:237–69. Search in Google Scholar
 Flowerdew M, Millar IL, Vaughan A, Horstwood M, Fanning C. The source of granitic gneisses and migmatites in the Antarctic Peninsula: a combined U–Pb SHRIMP and laser ablation Hf isotope study of complex zircons. Contribut Mineral Petrol. 2006;151:751–68. Search in Google Scholar
 Zeh A, Gerdes A, Barton Jr J, Klemd R. U–Th–Pb and Lu–Hf systematics of zircon from TTG’s, leucosomes, meta-anorthosites and quartzites of the Limpopo Belt (South Africa): constraints for the formation, recycling and metamorphism of Palaeoarchaean crust. Precambrian Res. 2010;179:50–68. Search in Google Scholar
 Belousova E, Kostitsyn Y, Griffin WL, Begg GC, O’Reilly SY, Pearson NJ. The growth of the continental crust: constraints from zircon Hf-isotope data. Lithos. 2010;119:457–66. Search in Google Scholar
 Kinny PD, Maas R. Lu–Hf and Sm–Nd isotope systems in zircon. Rev Mineral Geochem. 2003;53:327–41. Search in Google Scholar
 Thirlwall MF, Walder AJ. In situ hafnium isotope ratio analysis of zircon by inductively coupled plasma multiple collector mass spectrometry. Chem Geol. 1995;122:241–7. Search in Google Scholar
 Harrison T, Blichert-Toft J, Müller W, Albarede F, Holden P, Mojzsis SJ. Heterogeneous Hadean hafnium: evidence of continental crust at 4.4 to 4.5 Ga. Science. 2005;310:1947–50. Search in Google Scholar
 Kemp A, Hawkesworth C, Paterson B, Kinny P. Episodic growth of the Gondwana supercontinent from hafnium and oxygen isotopes in zircon. Nature. 2006;439:580. Search in Google Scholar
 Anczkiewicz R. Structural and geochronological study of the India-Kohistan arc collision, lower Swat region of Pakistan, NW Himalaya. Unpublished PhD thesis, # 12771, ETH Zürich, Switzerland; 1998. Search in Google Scholar
 DiPietro JA, Isachsen CE. U–Pb zircon ages from the Indian plate in northwest Pakistan and their significance to Himalayan and pre-Himalayan geologic history. Tectonics. 2001;20:510–25. Search in Google Scholar
 Taylor SR, McLennan SM. The continental crust: its composition and evolution. Oxford, UK: Blackwell; 1985. p. 1–312. Search in Google Scholar
 Tahirkheli RK. Geology of Kohistan and adjoining Eurasian and Indo-Pakistan continents, Pakistan. Geol Bull Univ Peshawar. 1979;11:1–30. Search in Google Scholar
 Kazmi AH, Jan MQ. Geology and tectonics of Pakistan. Karachi, Pakistan: Graphic Publishers; 1997. p. 1–554. Search in Google Scholar
 Searle M, Khan MA, Fraser J, Gough S, Jan MQ. The tectonic evolution of the Kohistan Karakoram collision belt along the Karakoram Highway transect, north Pakistan. Tectonics. 1999;18:929–49. Search in Google Scholar
 Searle M, Treloar P. Was Late Cretaceous–Paleocene obduction of ophiolite complexes the primary cause of crustal thickening and regional metamorphism in the Pakistan Himalaya? vol. 338, London: Geological Society; Special Publications, 2010. p. 345–59. Search in Google Scholar
 Burg J-P. The Asia–Kohistan–India collision: review and discussion. In Arc-continent collision. Germany: Springer-Verlag Berlin Heidelberg; 2011. p. 279–309. Search in Google Scholar
 DiPietro J, Lawrence R. Himalayan structure and metamorphism south of the Main Mantle thrust, Lower Swat, Pakistan. J Metamorphic Geol. 1991;9:481–95. Search in Google Scholar
 Greco A, Spencer DA. A section through the Indian plate, Kaghan valley, NW Himalaya, Pakistan. vol. 74. London: Geological Society; Special Publications, 1993. p. 221–36. Search in Google Scholar
 DiPietro JA, Pogue KR. Tectonostratigraphic subdivisions of the Himalaya: a view from the west. Tectonics. 2004;23(5):1–20. Search in Google Scholar
 Hussain A, DiPietro JA, Pogue AR, Ahmad I. Geological map of the 43B degree sheet, NWFP. Pakistan: Geological Survey of Pakistan; 2004. Search in Google Scholar
 Argles T, Foster G, Whittington A, Harris N, George M. Isotope studies reveal a complete Himalayan section in the Nanga Parbat syntaxis. Geology. 2003;31:1109–12. Search in Google Scholar
 Kaneko Y, Katayama I, Yamamoto H, Misawa K, Ishikawa M, Rehman H, et al. Timing of Himalayan ultrahigh pressure metamorphism: Sinking rate and subduction angle of the Indian continental crust beneath Asia. J Metamorphic Geol. 2003;21:589–99. Search in Google Scholar
 Ahmad I, Jan M, DiPietro J. Age and Tectonic implications of granitoid rocks from the Indian plate of Northern Pakistan. J Virtual Explorer. 2003;11:13–34. Search in Google Scholar
 Shams FA. Granites of Himalayas, Karakorum, and Hindu Kush. Institute of Geology, Punjab, University; 1983. Search in Google Scholar
 Maluski H, Matte P. Ages of alpine tectonometamorphic events in the northwestern Himalaya (northern Pakistan) by 39Ar/40Ar method. Tectonics. 1984;3:1–18. Search in Google Scholar
 DiPietro JA. Metamorphic pressure temperature conditions of Indian Plate rocks south of the main mantle thrust, lower swat, Pakistan. Tectonics. 1991;10:742–57. Search in Google Scholar
 Miller JS, Matzel JEP, Miller CF, Burgess SD, Miller RB. Zircon growth and recycling during the assembly of large, composite arc plutons. J Volcanol Geotherm Res. 2007;167(1–4):282–99. Search in Google Scholar
 Wiedenbeck M, Alle P, Corfu F, Griffin W, Meier M, Oberli F, et al. Three natural zircon standards for U–Th–Pb, Lu–Hf, trace element and REE analyses. Geostand Newslett. 1995;19:1–23. Search in Google Scholar
 Liu Y, Gao S, Hu Z, Gao C, Zong K, Wang D. 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 from mantle xenoliths. J Petrol. 2010;51:537–71. Search in Google Scholar
 Andersen T. Correction of common lead in U–Pb analyses that do not report 204Pb. Chem Geol. 2002;192:59–79. Search in Google Scholar
 Ludwig KR. User’s manual for isoplot 3.00, a geochronlogical toolkit for microsoft excel. Berkeley Geochronol Cent Spec Publ. 2003;4:25–32. Search in Google Scholar
 Woodhead JD, Hergt JM. A preliminary appraisal of seven natural zircon reference materials for in situ Hf isotope determination. GeostGeoanalytical Res. 2005;29:183–95. Search in Google Scholar
 Söderlund U, Patchett PJ, Vervoort JD, Isachsen CE. The 176Lu decay constant determined by Lu–Hf and U–Pb isotope systematics of Precambrian mafic intrusions. Earth Planet Sci Lett. 2004;219:311–24. Search in Google Scholar
 Bouvier A, Vervoort JD, Patchett PJ. The Lu–Hf and Sm–Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth Planet Sci Lett. 2008;273(1–2):48–57. Search in Google Scholar
 Williams I, Claesson S. Isotopic evidence for the Precambrian provenance and Caledonian metamorphism of high grade paragneisses from the Seve Nappes, Scandinavian Caledonides. Contribut Mineral Petrol. 1987;97:205–17. Search in Google Scholar
 Hoskin PW, Ireland TR. Rare earth element chemistry of zircon and its use as a provenance indicator. Geology. 2000;28:627–30. Search in Google Scholar
 Poller U, Huth J, Hoppe P, Williams IS. REE, U, Th, and Hf distribution in zircon from western Carpathian Variscan granitoids: a combined cathodoluminescence and ion microprobe study. Am J Sci. 2001;301:858–67. Search in Google Scholar
 Sun SS, McDonough WF. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. London: Geological Society; Special Publications, 1989. vol. 42, p. 313–45. Search in Google Scholar
 Hoskin PW. Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the Jack Hills, Australia. Geochim Cosmochim Acta. 2005;69:637–48. Search in Google Scholar
 Chen RX, Zheng YF, Xie L. Metamorphic growth and recrystallization of zircon: distinction by simultaneous in situ analyses of trace elements, U–Th–Pb and Lu–Hf isotopes in zircons from eclogite-facies rocks in the Sulu orogen. Lithos. 2010;114:132–54. Search in Google Scholar
 Wang Q, Zhu DC, Zhao Z, Guan Q, Zhang XQ, Sui QL, et al. Magmatic zircons from I-, S- and A-type granitoids in Tibet: trace elements characteristics and their application to detrital zircon provenance study. J Asian Earth Sci. 2012;53:59–66. Search in Google Scholar
 McDonough WF, Sun S-S. The composition of the Earth. Chem Geol. 1995;120:223–53. Search in Google Scholar
 Hawkesworth C, Kemp A. Using hafnium and oxygen isotopes in zircons to unravel the record of crustal evolution. Chem Geol. 2006;226:144–62. Search in Google Scholar
 Yang J, Cawood PA, Du Y, Huang H, Huang H, Tao P. Large Igneous Province and magmatic arc sourced Permian–Triassic volcanogenic sediments in China. Sediment Geol. 2012;261:120–31. Search in Google Scholar
 Claiborne LL, Miller CF, Wooden JL. Trace element composition of igneous zircon: a thermal and compositional record of the accumulation and evolution of a large silicic batholith, Spirit Mountain, Nevada. Contribut Mineral Petrol. 2010;160:511–31. Search in Google Scholar
 Kelemen PB, Hanghøj K, Greene A. One view of the geochemistry of subduction-related magmatic arcs, with an emphasis on primitive andesite and lower crust. Treatise Geochem. 2003;3:659. Search in Google Scholar
 Rudnick RL, Gao S. Composition of the continental crust. Treatise Geochem. 2003;3:659. Search in Google Scholar
 Bea F. Residence of REE, Y, Th and U in granites and crustal protoliths, implications for the chemistry of crustal melts. J Petrol. 1996;37:521–52. Search in Google Scholar
 Pearce JA, Harris NB, Tindle AG. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. J Petrol. 1984;25:956–83. Search in Google Scholar
 Anders E, Grevesse N. Abundances of the elements: meteoritic and solar. Geochim Cosmochim Acta. 1989;53:197–214. Search in Google Scholar
 Cavosie AJ, Valley JW, Kita NT, Spicuzza MJ, Ushikubo T, Wilde SA. The origin of high δ18O zircons: marbles, megacrysts, and metamorphism. Contribut Mineral Petrol. 2011;162:961–74. Search in Google Scholar
 Bindeman IN, Serebryakov NS, Schmitt AK, Vazquez JA, Guan Y, Azimov PY, et al. Field andmicroanalytical isotopic investigation of ultradepleted in 18O Paleoproterozoic “Slushball Earth” rocks from Karelia, Russia. Geosphere. 2014;10:308–39. Search in Google Scholar
 Xia QX, Zheng YF, Hu Z. Trace elements in zircon and coexisting minerals from low-T/UHP metagranite in the Dabie orogen: implications for action of supercritical fluid during continental subduction-zone metamorphism. Lithos. 2010;14(3):385–412. Search in Google Scholar
 Fu B, Mernagh TP, Kita NT, Kemp AI, Valley JW. Distinguishing magmatic zircon from hydrothermal zircon: a case study from the Gidginbung high-sulphidation Au–Ag–(Cu) deposit, SE Australia. Chem Geol. 2009;259:131–42. Search in Google Scholar
 Moecher D, McDowell S, Samson SD, Miller C. Ti-in-zircon thermometry and crystallization modeling support hot Grenville granite hypothesis. Geology. 2014;42:267–70. Search in Google Scholar
 Ewing TA, Hermann J, Rubatto D. The robustness of the Zr-in-rutile and Ti-in-zircon thermometers during high-temperature metamorphism (Ivrea-Verbano Zone, northern Italy). Contribut Mineral Petrol. 2013;165:757–79. Search in Google Scholar
 Liu Y-C, Deng L-P, Gu X-F, Groppo C, Rolfo F. Application of Ti-in-zircon and Zr-in-rutile thermometers to constrain high-temperature metamorphism in eclogites from the Dabie orogen, central China. Gondwana Res. 2015;27:410–23. Search in Google Scholar
 Clemens J. S-type granitic magmas – petrogenetic issues, models and evidence. Earth-Sci Rev. 2003;61:1–18. Search in Google Scholar
 Amelin Y, Lee D-C, Halliday A. Early-middle Archaean crustal evolution deduced from Lu–Hf and U–Pb isotopic studies of single zircon grains. Geochim Cosmochim Acta. 2000;64:4205–25. Search in Google Scholar
 Corfu F, Stott GM. Age and petrogenesis of two late Archean magmatic suites, northwestern Superior Province, Canada: zircon U–Pb and Lu–Hf isotopic relations. J Petrol. 34, 1993. 817–38. Search in Google Scholar
 Vervoort J, Patchett P, Gehrels GE, Nutman A. Constraints on early Earth differentiation from hafnium and neodymium isotopes. Nature. 1996;379:624. Search in Google Scholar
 Vervoort JD, Blichert-Toft J. Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochim Cosmochim Acta. 1999;63:533–56. Search in Google Scholar
© 2020 Lawangin Sheikh et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.