Petrogenesis of Early Triassic Felsic Volcanic Rocks in the East Kunlun Orogen, Northern Tibet: Implications for the Paleo-Tethyan Tectonic and Crustal Evolution

: The felsic volcanic rocks in orogenic belts are vital probes to understand the tectonic evolution and continental crust growth. Here, we present a comprehensive study on the zircon U–Pb geochronology, whole-rock geochemistry, and zircon Lu-Hf isotopes of Early felsic volcanic rocks from the Hongshuichuan Formation, East Kunlun Orogen, Northern Tibet, aiming to explore their petrogenesis and implications for the Paleo-Tethyan orogeny and crustal evolution. The studied felsic volcanics comprise rhyolite porphyry and rhyolite, exhibiting coeval zircon U–Pb ages of ca. 247–251 Ma. Rhyolite porphyries show metaluminous to peraluminous nature (A/CNK = 0.88–1.24) with high SiO 2 contents (72.1–78.9 wt%) and moderate Mg # values (22–40), and they display enrichment of LREE with (La/Yb) N ratios of 6.02–17.9 and depletion of high ﬁeld strength elements. In comparison, the rhyolites are strongly peraluminous (A/CNK = 1.09–1.74) with high SiO 2 contents (71.7–74.3 wt%) and high Mg # values (43–52) and are also enriched in LREE ((La/Yb) N of 6.65–18.4) and depleted in HFSE (e.g., Nb, Ta, Ti). Combining with their different zircon Lu-Hf isotopes, i.e., enriched isotopes for the rhyolite porphyries ( ε Hf(t) = − 7.3 to − 3.8) and depleted Hf isotopes for the rhyolites ( ε Hf = − 0.6 to +3.0), we interpret that the studied rhyolite porphyries and rhyolites were derived by partial melting of Mesoproterozoic metagreywacke sources followed by plagioclase-dominated fractional crystallization, but the latter shows the signiﬁcant contribution of crust–mantle magma mixing. The mixed mantle-derived magma comes from an enriched lithospheric mantle source that had been metasomatized by subduction-related ﬂuids. Combining with other geological evidence, we propose that the studied Early Triassic felsic volcanic rocks were formed in a subduction arc setting, and the reworking of ancient continental crust with crust–mantle magma mixing is the major mechanism of crustal evolution in the East Kunlun Paleo-Tethyan orogenic belt.


Introduction
As the most archetypal products of a long-lived orogeny, felsic volcanics record the continental crust growth and tectonic evolution of Paleo-orogenic belts [1][2][3].Several petrogenetic models have been proposed for the origin of felsic volcanic rocks in orogenic belts: (1) advanced degrees of fractional crystallization of mantle-derived mafic magmas [4,5], (2) partial melting of continental crust triggered by the underplating of mantle-derived mafic magmas [6][7][8], and (3) hybridization of crustal felsic melts with mantle-derived mafic melts [9][10][11][12].Thus, identifying the petrogenesis of felsic volcanics is the key to reveal the mechanism of crust-mantle interaction as well as continental evolution and tectonic history of Paleo-orogenic belts.
The East Kunlun Orogen Belt (EKOB), one typical tectonic-magmatic belt in the northern Tibet Plateau, has experienced Paleo-Tethyan subduction collision orogeny as recorded by large-scale Permian-Triassic felsic magmatism.Although much research has been carried out on the Paleo-Tethyan felsic magmatism in the EKOB, its petrogenesis and geodynamic background are still controversial.Some studies propose that most Early to Middle Triassic felsic volcanics are formed in a subduction setting and the Kunlun Paleo-Tethyan ocean was finally closed in Late Triassic, while other studies emphasize that the Kunlun Paleo-Tethyan ocean was closed in Late Permian, and the Triassic felsic volcanics all are generated in a post-collisional setting.Meanwhile, some scholars even proposed that the Paleo-Tethyan oceanic subduction lasted until the Late Triassic or Jurassic [2,[13][14][15].Due to the rare occurrence of Early Triassic felsic volcanic rocks in the EKOB, previous studies mainly focus on the large-scale developed Late Triassic volcanic rocks, including abundant andesitic to rhyolitic lavas and volcaniclastic rocks [16,17].However, the petrogenesis and geodynamic background of the Early Triassic felsic volcanic magmatism in the EKOB are still unclear, which restricts our understanding of the East Kunlun Paleo-Tethyan orogeny.
Herein, we carried out systematic studies on zircon U-Pb ages, bulk-rock major and trace element compositions, and zircon Lu-Hf isotopes of Early Triassic felsic volcanic rocks in the EKOB to discuss their petrogenesis.These data, as well as the published dataset of the Late Triassic felsic rocks, help us reveal the tectonic setting of the Early Triassic data for volcanic rocks and reconstruct the tectonic and crustal evolutionary history of the East Kunlun Paleo-Tethyan orogen.

Regional Geologic Background
The EKOB is located in the northern Tibet Plateau, NW China, adjoins the Qaidam block in the north, the Bayan Har Orogen in the south, and the West Qinling Orogen in the east (Figure 1a).The orogen extends for nearly 1500 km from east to west along the northern margin of the Tibetan Plateau, with a width between 50 km and 200 km.Based on the central and south Kunlun faults, the EKOB can be divided into the Northeast Kunlun terrane and the Southeast Kunlun terrane.Due to Paleo-Tethyan oceanic subduction and subsequent collision, large-scale Late Permian to Triassic felsic magmatism occurred in the EKOB (Figure 1b) [2,14].Permian-Triassic felsic intrusive rocks are dominated by granites and granodiorites with abundant mafic magmatic enclaves [13,[18][19][20][21], while the felsic volcanic rocks are mainly Late Triassic in age (e.g., Elashan Formation) [17], with smaller volumes of Late Permian to Early Triassic age [22].The crustal basement of the EKOB is represented by the Paleo-and Mesoproterozoic Jinshuikou Group comprising the lower Baishahe and upper Xiaomiao formations.The lower Baishahe Formation consists of marbles, gneisses, migmatites, and amphibolites.The upper Xiaomiao Formation comprises marbles, gneisses, greenschists, and quartzites [24][25][26][27].The basement is covered by the Permian-Triassic volcanic-sedimentary strata, such as the Late Permian Gequ Formation, the Early Triassic Hongshuichuan Formation, and the Middle Triassic Naocangjiangou Formation.Among them, the Early Triassic Hongshuichuan Formation is mainly composed of quartz sandstone, bioclastic limestone, lavas, and pyroclastic rocks.This set of strata is mainly composed of limestone and sandstone; lavas mainly occur in it in the form of volcanic strata in different stratigraphic units of the Hongshui Chuan Formation, which is the product of a marine eruptive environment.The volcanic lavas in this study are sampled from the Hongshuichuan Formation, which mainly comprises rhyolite and rhyolite porphyry and is interbedded with contemporaneous early Mesozoic sedimentary strata (Figure 1c).

Sampling and Petrography
The Early Triassic felsic volcanic rocks in the Hongshuichuan Formation overlie limestone.Rhyolite is located above rhyolite porphyry and is separated by quartz sandstone and bioclastic limestone.The rhyolite porphyries are greyish-green and have spherulitic or porphyritic texture (Figure 2c).The phenocrysts are K-feldspar (20-25 vol%), quartz (5- The crustal basement of the EKOB is represented by the Paleo-and Mesoproterozoic Jinshuikou Group comprising the lower Baishahe and upper Xiaomiao formations.The lower Baishahe Formation consists of marbles, gneisses, migmatites, and amphibolites.The upper Xiaomiao Formation comprises marbles, gneisses, greenschists, and quartzites [24][25][26][27].The basement is covered by the Permian-Triassic volcanic-sedimentary strata, such as the Late Permian Gequ Formation, the Early Triassic Hongshuichuan Formation, and the Middle Triassic Naocangjiangou Formation.Among them, the Early Triassic Hongshuichuan Formation is mainly composed of quartz sandstone, bioclastic limestone, lavas, and pyroclastic rocks.This set of strata is mainly composed of limestone and sandstone; lavas mainly occur in it in the form of volcanic strata in different stratigraphic units of the Hongshui Chuan Formation, which is the product of a marine eruptive environment.The volcanic lavas in this study are sampled from the Hongshuichuan Formation, which mainly comprises rhyolite and rhyolite porphyry and is interbedded with contemporaneous early Mesozoic sedimentary strata (Figure 1c).

LA-ICP-MS Zircon U-Pb Dating
The studied samples were collected from fresh outcrops and zircons were separated by heavy liquid and magnetic methods.Zircon grains were photographed with an optical microscope and the internal structures were analyzed by cathodoluminescence (CL).Zircon U-Pb isotope and trace element concentration were determined using a 193 nm Geo Las Pro and an Agilent 7500a ICP-MS with a laser spot beam of 32 µm at the State Key Laboratory of Geological Processes and Mineral Resources (GPMR), China University of Geosciences, Wuhan.For the external standards and calibration of mass discrimination and isotope fractionation, zircon 91,500 and zircon GJ-1 were used, respectively.A background acquisition time of 18-20 s and a data acquisition time of 50 s was adopted for each measurement.Helium was applied as a carrier gas, and Ar was used as the makeup gas and mixed with the carrier gas via a T-connector before entering the ICP source.Nitrogen was added to the central gas flow (Ar + He) of the Ar plasma to lower the detection limits and improve precision.Concordia diagrams and weighted mean ages were made by Isoplot/Ex_ver3 [28].Data were processed using ICPMSDataCal [29].The detailed operating conditions for the laser system and ICP-MS instrument are described by Liu et al. [29].

LA-ICP-MS Zircon U-Pb Dating
The studied samples were collected from fresh outcrops and zircons were separated by heavy liquid and magnetic methods.Zircon grains were photographed with an optical microscope and the internal structures were analyzed by cathodoluminescence (CL).Zircon U-Pb isotope and trace element concentration were determined using a 193 nm Geo Las Pro and an Agilent 7500a ICP-MS with a laser spot beam of 32 µm at the State Key Laboratory of Geological Processes and Mineral Resources (GPMR), China University of Geosciences, Wuhan.For the external standards and calibration of mass discrimination and isotope fractionation, zircon 91,500 and zircon GJ-1 were used, respectively.A background acquisition time of 18-20 s and a data acquisition time of 50 s was adopted for each measurement.Helium was applied as a carrier gas, and Ar was used as the makeup gas and mixed with the carrier gas via a T-connector before entering the ICP source.Nitrogen was added to the central gas flow (Ar + He) of the Ar plasma to lower the detection limits and improve precision.Concordia diagrams and weighted mean ages were made by Isoplot/Ex_ver3 [28].Data were processed using ICPMSDataCal [29].The detailed operating conditions for the laser system and ICP-MS instrument are described by Liu et al. [29].

Whole-Rock Major and Trace Element Analyses
Whole-rock samples were crushed in a corundum jaw crusher (to 60 mesh).About 60 g of this material was powdered in an agate ring mill to <200 mesh for whole-rock geochemical analysis.The major element analysis was conducted by standard X-ray fluorescence (XRF) methods, using a Shimadzu Sequential 1800 spectrometer at the GPMR.Precision is <3% and accuracy is <4% for the major elements.The detailed techniques for major element analysis were described by Ma et al. [30].The trace element compositions were determined using an Agilent 7700a ICP-MS with analytical precision ranging from 5 to 10%.The samples were digested by HF +HNO 3 in Teflon bombs.Analyses of USGS standards (AGV-2, BHVO-2, BCR-2 and RGM-2) indicate accuracy better than 5%-10% for most trac elements.The analysis was conducted using ICP-MS (Agilent 7700) and ICP-AES (Agilent VISIA) methods after the rock powder (200 mesh) had been dissolved by HF-HNO 3 -HClO 4 -HCl acid digestion.The relative standard deviation of the elemental concentrations is close to or below 10%.The sample digestion procedures and ICP-MS methods are as described by Liu et al. [31].

Zircon Lu-Hf Isotope Analyses
In-situ Lu-Hf isotope analysis was conducted at the State Key Laboratory of GPMR, China, using a Neptune Plus MC-ICP-MS equipped with a Geolas 193 nm ArF Laser Ablation System (LAS).Measurement was completed on the ablation spots close to those for the U-Pb dating with a spot beam of 44 µm.Acquisition of background and ablation signals were completed for 20 s and 50 s, respectively.Detailed operating conditions for the laser ablation system and the MC-ICP-MS instrument and analytical method are the same as described by Hu et al. [32].Off-line selection and integration of analytical signals and mass bias calibrations were performed using ICPMSDataCal [31].

Zircon Trace Elements
The data for zircon trace elements are shown in Table S2.As shown in a chondritenormalized REE plot, all analyses are LREE depleted and HREE enriched and show

Zircon Trace Elements
The data for zircon trace elements are shown in Table S2.As shown in a chondritenormalized REE plot, all analyses are LREE depleted and HREE enriched and show positive Ce and negative Eu anomalies, chemically identical to typical magmatic zircons (Figure 5a) [34].Zircons from rhyolite contain higher Hf (~9066 ppm), Th/U (~0.

Petrogenetic Type of the Felsic Volcanic Rocks
Granitoids and their eruptive equivalents are usually divided into I-, S-, M-, and Atype granitic rocks based on their petrogenetic mechanisms.In general, I-type granitic rocks are usually derived by the crust-mantle interaction or partial melting of metamorphic igneous rocks [14,44], S-type granitic rocks are formed by melting of aluminum-rich metasedimentary rocks associated with collision or post-collision setting [45], M-type granitic rocks are mainly derived by fractional crystallization of mantle-derived magma, while A-type granitic rocks mainly originated from various magma sources in extensional environment [16].Revealing the petrogenetic types of granitoids and their volcanic equivalents is the key to understanding their petrogenesis and geodynamic implications.Based on petrography and geochemical studies, we propose that the rhyolite and rhyolite porphyry are I-type granitic rocks.The critical supporting evidence is as follows: Firstly, the studied rhyolite and rhyolite porphyry all display low Ce and HFSEs (Zr + Nb + Y + Ce) contents, low Ga/Al and FeO T /MgO ratios, contrary to the Late Triassic felsic volcanics, and all plot in the I-type granitic field according to various chemical classification diagrams (Figure 9a,b) [46].This identification is also corroborated by the anticorrelation between SiO2 and Al2O3 contents (Figure 9c,d).Additionally, the rhyolite and rhyolite porphyry show higher CaO contents (average of 1.11 wt% and 1.28 wt%) relative to the A-type granitic rocks (0.75 wt%) but closely match the I-type granitic rocks (1.71 wt%, [46][47][48]).

Petrogenetic Type of the Felsic Volcanic Rocks
Granitoids and their eruptive equivalents are usually divided into I-, S-, M-, and A-type granitic rocks based on their petrogenetic mechanisms.In general, I-type granitic rocks are usually derived by the crust-mantle interaction or partial melting of metamorphic igneous rocks [14,44], S-type granitic rocks are formed by melting of aluminum-rich metasedimentary rocks associated with collision or post-collision setting [45], M-type granitic rocks are mainly derived by fractional crystallization of mantle-derived magma, while A-type granitic rocks mainly originated from various magma sources in extensional environment [16].Revealing the petrogenetic types of granitoids and their volcanic equivalents is the key to understanding their petrogenesis and geodynamic implications.Based on petrography and geochemical studies, we propose that the rhyolite and rhyolite porphyry are I-type granitic rocks.The critical supporting evidence is as follows: Firstly, the studied rhyolite and rhyolite porphyry all display low Ce and HFSEs (Zr + Nb + Y + Ce) contents, low Ga/Al and FeO T /MgO ratios, contrary to the Late Triassic felsic volcanics, and all plot in the I-type granitic field according to various chemical classification diagrams (Figure 9a,b) [46].This identification is also corroborated by the anti-correlation between SiO 2 and Al 2 O 3 contents (Figure 9c,d).Additionally, the rhyolite and rhyolite porphyry show higher CaO contents (average of 1.11 wt% and 1.28 wt%) relative to the A-type granitic rocks (0.75 wt%) but closely match the I-type granitic rocks (1.71 wt%, [46][47][48]).
The studied rhyolite porphyry has high contents of SiO2, low contents of MgO, TiO2, and FeO T , and shows enrichment of large ion lithophile elements (Rb, Th, U, and K) and depletion of high field strength elements (Nb, Ta, Zr, and Hf) with remarkable depletion of Sr and Eu (Figure 10), which are typical features of crustal origin.In addition, the K2O value of rhyolite porphyry varies widely, which we believe is due to the hydration of glass [49].The rhyolite may have undergone weathering after its formation.This is likely to lead to alter in the original glass matrix, and most samples absorbed K and lost Na.The ratios of incompatible trace elements (e.g., Th/Nb, Th/La, La/Nb, and Nb/Ta) usually remain constant during magma differentiation; thus, they are reflective of the nature of magma sources [50].The average ratios of incompatible trace elements for the studied rhyolite porphyry (e.g., Th/Nb = 1.47;Th/La = 0.47; La/Nb = 3.21; Nb/Ta = 12.4) are higher than those of the primitive mantle (Th/Nb = 0.18; Th/La = 0.13; La/Nb = 0.94; Nb/Ta = 16.0) but closely identical to the average ratios of the continental crust (Th/Nb = 0.44; Th/La = 0.20; La/Nb = 2.20; Nb/Ta = 11.0-12.0)[51,52].In addition, the studied rhyolite porphyry exhibits elemental geochemistry identical to those derived from pure crustal melting in a normal arc setting (Figure 11a,b).The rhyolite porphyry has high values of Al2O3/(FeOtot + MgO) and CaO/(MgO + FeOtot), which is consistent with the melt derived from the metagreywacke source as expected by the experimental petrology (Figure 11c,d).Isotopically, their negative zircon εHf(t) values (−7.3 to −3.8) are significantly lower than those of the mantle-derived rocks or the juvenile crustal rocks (Figure 8).These lines of evidence, in combination with their Mesoproterozoic two-stage Hf model ages (1.4-1.6 Ga), suggest that the rhyolite porphyry was derived by partial melting of a Mesoproterozoic metagreywacke source in the lower crust.
The studied rhyolite porphyry has high contents of SiO 2 , low contents of MgO, TiO 2 , and FeO T , and shows enrichment of large ion lithophile elements (Rb, Th, U, and K) and depletion of high field strength elements (Nb, Ta, Zr, and Hf) with remarkable depletion of Sr and Eu (Figure 10), which are typical features of crustal origin.In addition, the K 2 O value of rhyolite porphyry varies widely, which we believe is due to the hydration of glass [49].The rhyolite may have undergone weathering after its formation.This is likely to lead to alter in the original glass matrix, and most samples absorbed K and lost Na.The ratios of incompatible trace elements (e.g., Th/Nb, Th/La, La/Nb, and Nb/Ta) usually remain constant during magma differentiation; thus, they are reflective of the nature of magma sources [50].The average ratios of incompatible trace elements for the studied rhyolite porphyry (e.g., Th/Nb = 1.47;Th/La = 0.47; La/Nb = 3.21; Nb/Ta = 12.4) are higher than those of the primitive mantle (Th/Nb = 0.18; Th/La = 0.13; La/Nb = 0.94; Nb/Ta = 16.0) but closely identical to the average ratios of the continental crust (Th/Nb = 0.44; Th/La = 0.20; La/Nb = 2.20; Nb/Ta = 11.0-12.0)[51,52].In addition, the studied rhyolite porphyry exhibits elemental geochemistry identical to those derived from pure crustal melting in a normal arc setting (Figure 11a,b).The rhyolite porphyry has high values of Al 2 O 3 /(FeO tot + MgO) and CaO/(MgO + FeO tot ), which is consistent with the melt derived from the metagreywacke source as expected by the experimental petrology (Figure 11c,d).Isotopically, their negative zircon εHf(t) values (−7.3 to −3.8) are significantly lower than those of the mantle-derived rocks or the juvenile crustal rocks (Figure 8).These lines of evidence, in combination with their Mesoproterozoic two-stage Hf model ages (1.4-1.6 Ga), suggest that the rhyolite porphyry was derived by partial melting of a Mesoproterozoic metagreywacke source in the lower crust.Similarly, the relatively high ratios of Al2O3/(FeOtot + MgO) and CaO/(MgO + FeOtot) suggest that the studied rhyolite is mainly derived from a metagreywacke source (Figure 11c,d).However, the rhyolite has high contents of SiO2 (71.7-74.3wt%) and high values of Mg # (43)(44)(45)(46)(47)(48)(49)(50)(51)(52)(53).This, as well as their juvenile or depleted zircon Lu-Hf isotopes, suggests that mantle materials, rather than pure crustal melts, contributed to their origins (Figure 11a).The εHf(t) values of zircons from the rhyolite (−0.6 to 3.0) are much higher than the Triassic granitoids derived from ancient continental crust in the EKOB (average εHf(t) values of −5.5; [62,63]) but close to the rocks derived from the enriched mantle or crust-mantle   Similarly, the relatively high ratios of Al2O3/(FeOtot + MgO) and CaO/(MgO + FeOtot) suggest that the studied rhyolite is mainly derived from a metagreywacke source (Figure 11c,d).However, the rhyolite has high contents of SiO2 (71.7-74.3wt%) and high values of Mg # (43)(44)(45)(46)(47)(48)(49)(50)(51)(52)(53).This, as well as their juvenile or depleted zircon Lu-Hf isotopes, suggests that mantle materials, rather than pure crustal melts, contributed to their origins (Figure 11a).The εHf(t) values of zircons from the rhyolite (−0.6 to 3.0) are much higher than the Triassic granitoids derived from ancient continental crust in the EKOB (average εHf(t) values of −5.5; [62,63]) but close to the rocks derived from the enriched mantle or crust-mantle Similarly, the relatively high ratios of Al 2 O 3 /(FeO tot + MgO) and CaO/(MgO + FeO tot ) suggest that the studied rhyolite is mainly derived from a metagreywacke source (Figure 11c,d).However, the rhyolite has high contents of SiO 2 (71.7-74.3wt%) and high values of Mg # (43)(44)(45)(46)(47)(48)(49)(50)(51)(52)(53).This, as well as their juvenile or depleted zircon Lu-Hf isotopes, suggests that mantle materials, rather than pure crustal melts, contributed to their origins (Figure 11a).The εHf(t) values of zircons from the rhyolite (−0.6 to 3.0) are much higher than the Triassic granitoids derived from ancient continental crust in the EKOB (average εHf(t) values of −5.5; [62,63]) but close to the rocks derived from the enriched mantle or crust-mantle magma mixing, further indicating a significant contribution of mantle melts probably derived from the enriched lithospheric mantle (Figure 8).Petrographic studies show that most plagioclase phenocrysts in rhyolite exhibit zoning texture, which is contrasted to the rhyolite porphyry.The widespread presence of zoning textures of plagioclase phenocryst in volcanic rocks is commonly due to the magma mixing between mafic and felsic magmas in open magma systems [63,64].
The rhyolites have similar Hf isotopic compositions with Early Triassic mafic rocks derived from enriched mantle in East Kunlun.(~240Ma gabbro:εHf(t) = −1.8-−0.6 [44]), indicating a significant addition of mantle material probably derived from enriched lithospheric mantle.Rhyolites are LILE enriched and HFSE depleted and show obvious Nb-Ta-Ti troughs in primitive mantle-normalized trace element variation diagrams, exhibiting geochemical features of arc-related magmas, which implies an enrichment mechanism that was related to subducted components.The relative enrichments of Rb, Th, Ta, and Ba (Figure 7), and the obvious depletion in Yb and Y, indicate that mantle below a continental arc, enriched in fluid, could account for these geochemical features.Therefore, mantle-derived magma comes from an enriched lithospheric mantle source that had been metasomatized by subduction-related fluids.
Thus, we propose that the studied rhyolite was derived from the melting of Mesoproterozoic metagraywacke in the lower continental crust with significant involvement of enriched lithospheric mantle -derived magma.

Fractional Crystallization
The rhyolite porphyry and rhyolite display typical geochemical characteristics of high silica rhyolite with extremely high SiO 2 contents (>70 wt%; [65]), thus suggesting that the magma has undergone significant evolution.Fractional crystallization during magma evolution is attested by the marked depletions concerning Ba, Nb, Sr, P, Ti, and Eu on the trace and REE distribution diagrams (Figure 7).Negative Nb-Ti anomalies mirror fractionation of titanium-bearing phases, such as ilmenite or titanite, while negative P anomalies are consistent with apatite fractionation, as also evident by the negative relationship with P 2 O 5 and SiO 2 .Plagioclase fractionation is demonstrated by the negative correlation between SiO 2 and Al 2 O 3 (Figure 10c).Plagioclase is an Al-rich mineral; when it is removed from melt, the magma is depleted in Al 2 O 3 .These scenarios explain the pronounced negative Eu anomalies (Figure 9b) because Ca 2+ is substituted by Eu 2+ in plagioclase [66].A decrease in Al 2 O 3 as a function of plagioclase separation has been proven experimentally [67].The linear relationship in (Figure 10c,d) further demonstrates the separation of plagioclase during magma evolution.
Zircon geochemistry may allow for a clear picture of magma evolutionary processes [68].Zircon from the studied rhyolite and rhyolite porphyry exhibit a significant fractionation trend in U/Yb vs. Eu/Eu* diagram (Figure 5d).The Eu/Eu* ratios negatively correlate with Hf (Figure 5c), which can be explained by the crystallization of plagioclase.The diverse zircon compositions along the observed evolutionary path thus reflect amphibole, apatite, and/or titanite fractionation (Figure 5b).

Geodynamic Implications
Recently, Paleozoic ophiolites with ages of ca.355-308 Ma have been identified in the East Kunlun Paleo-Tethyan oceanic suture zone, indicating the opening of the Tethyan ocean not later than Carboniferous [69-72].The regional angular unconformity between the Upper Permian Gequ Formation and the Upper Carboniferous Haoteluowa Formation in the EKOB indicates the occurrence of regional tectonic disturbance during the Late Permian [3].This regional angular unconformity, as well as the gradually developed Late Permian magmatic rocks with arc-like geochemical affinities, indicates that the Kunlun Paleo-Tethyan oceanic subduction in the EKOB initiated at Late Permian [25,73,74].However, the subduction duration and the closure time of the Paleo-Tethyan Ocean remain poorly constrained, which restricts our understanding of the Early Triassic magmatism and tectonics in the EKOB.
The low Sr/Y (1.59-15.5)and La/Yb (8.38-25.7)ratios as well as high Th/Yb ratios (3.76-14.8)and elevated Yb and Y content suggests that the studied Early Triassic felsic volcanic rocks were derived from a normal arc-related magmatic source (Figure 11b, [75]).The studied felsic volcanics are enriched in LILE and depleted in HFSE with negative Nb, Ta, P, and Ti anomalies (Figure 7b poorly constrained, which restricts our understanding of the Early Triassic magmatism and tectonics in the EKOB. The low Sr/Y (1.59-15.5)and La/Yb (8.38-25.7)ratios as well as high Th/Yb ratios (3.76-14.8)and elevated Yb and Y content suggests that the studied Early Triassic felsic volcanic rocks were derived from a normal arc-related magmatic source (Figure 11b, [75]).The studied felsic volcanics are enriched in LILE and depleted in HFSE with negative Nb, Ta, P, and Ti anomalies (Figure 7b  The above conclusion is also supported by stratigraphic evidence.Previous studies on the sedimentary strata in the EKOB show that the Late Permian Gequ Formation comprises conglomerate, sandstone, and calcilutite; the Early Triassic Hongshuichuan Formation is mainly composed of bioclastic limestone, lavas, and pyroclastic rocks; and the Middle Triassic Naocangjiangou Formation is mainly composed of sandstone and limestone [3,78,79].The studies on sedimentary stratigraphy suggest that these Late Permian to Middle Triassic strata containing marine sedimentary rocks all crop out along the southern margin of the proposed East Kunlun continental magmatic arc and are originally deposited in a forearc basin [78,79], recording the northward subduction of the East Kunlun Paleo-Tethyan Ocean.The Late Triassic sedimentary strata are composed of continental volcanic rocks and clastic rocks and are unconformably overlaid on the Late Permian to Middle Triassic strata [3,11,80], representing the closure of the East Kunlun Paleo-Tethyan Ocean and the following continental collision.The studied rhyolite porphyry and rhyolite are geochemically distinct from Late Triassic volcanic rocks in the EKOB, because Late Triassic volcanics have chemical affinity similar to typical within-plate rather than The above conclusion is also supported by stratigraphic evidence.Previous studies on the sedimentary strata in the EKOB show that the Late Permian Gequ Formation comprises conglomerate, sandstone, and calcilutite; the Early Triassic Hongshuichuan Formation is mainly composed of bioclastic limestone, lavas, and pyroclastic rocks; and the Middle Triassic Naocangjiangou Formation is mainly composed of sandstone and limestone [3,78,79].The studies on sedimentary stratigraphy suggest that these Late Permian to Middle Triassic strata containing marine sedimentary rocks all crop out along the southern margin of the proposed East Kunlun continental magmatic arc and are originally deposited in a forearc basin [78,79], recording the northward subduction of the East Kunlun Paleo-Tethyan Ocean.The Late Triassic sedimentary strata are composed of continental volcanic rocks and clastic rocks and are unconformably overlaid on the Late Permian to Middle Triassic strata [3,11,80], representing the closure of the East Kunlun Paleo-Tethyan Ocean and the following continental collision.The studied rhyolite porphyry and rhyolite are geochemically distinct from Late Triassic volcanic rocks in the EKOB, because Late Triassic volcanics have chemical affinity similar to typical within-plate rather than arc-related rocks (Figure 12, [11, 16,17]).Zhu et al. [16] and Shao et al. [17] propose that the Late Triassic volcanic rocks were more likely to generate in a within-plate extension setting.Thus, we propose that the East Kunlun Paleo-Tethyan Ocean subducted northward beneath the East Kunlun-Qaidam terrane during Late Permian-Middle Triassic.
The geochemical and zircon Lu-Hf isotopic studies show that the studied rhyolite porphyry and rhyolite originated from the melting of pure ancient continental crust, but the origin of the latter one has a certain contribution of crust-mantle magma mixing.Therefore, the reworking of the ancient continental lithosphere and magma mixing are important mechanisms for continental crust growth.In summary, we propose a tectonic and crustal evolution model for the EKOB during Early Triassic (Figure 13).Owing to the sinking slab of the Kunlun Paleo-Tethyan oceanic crust and lithospheric mantle during the Early Triassic, the melting of the metasomatized mantle wedge generated the mafic magma.The following underplating of the mafic magma heated the lower continental crust and triggered its melting to form the felsic magma.The felsic magma underwent fractionation and then directly erupted upward to form low-Mg # felsic volcanic rocks (i.e., rhyolite porphyry) with enriched zircon Hf isotopic compositions, while some felsic magmas mixed with mantle-derived mafic magmas to form high-Mg # felsic rocks (i.e., rhyolite) with plagioclase zoning and depleted zircon Hf isotopes.
plate rather than arc-related rocks (Figure 12, [11, 16,17]).Zhu et al. [16] and Shao et al. [17] propose that the Late Triassic volcanic rocks were more likely to generate in a within-plate extension setting.Thus, we propose that the East Kunlun Paleo-Tethyan Ocean subducted northward beneath the East Kunlun-Qaidam terrane during Late Permian-Middle Triassic.
The geochemical and zircon Lu-Hf isotopic studies show that the studied rhyolite porphyry and rhyolite originated from the melting of pure ancient continental crust, but the origin of the latter one has a certain contribution of crust-mantle magma mixing.Therefore, the reworking of the ancient continental lithosphere and magma mixing are important mechanisms for continental crust growth.In summary, we propose a tectonic and crustal evolution model for the EKOB during Early Triassic (Figure 13).Owing to the sinking slab of the Kunlun Paleo-Tethyan oceanic crust and lithospheric mantle during the Early Triassic, the melting of the metasomatized mantle wedge generated the mafic magma.The following underplating of the mafic magma heated the lower continental crust and triggered its melting to form the felsic magma.The felsic magma underwent fractionation and then directly erupted upward to form low-Mg # felsic volcanic rocks (i.e., rhyolite porphyry) with enriched zircon Hf isotopic compositions, while some felsic magmas mixed with mantle-derived mafic magmas to form high-Mg # felsic rocks (i.e., rhyolite) with plagioclase zoning and depleted zircon Hf isotopes.Elemental and isotopic geochemistry studies indicate that the rhyolite porphyry and rhyolite have typical lower crust-derived characteristics, and they are generated by partial melting of Mesoproterozoic metagreywacke sources with a certain amount of mantle-derived mafic magmas involved in the origin of the high Mg # rhyolite.

3.
The studied ca.250 Ma felsic volcanic rocks were formed during the northward subduction of the East Kunlun Paleo-Tethyan oceanic slab.This study suggests that the reworking of ancient continental crust and crust-mantle magma mixing are important mechanisms for the evolution of the orogenic continental crust.

Figure 1 .
Figure 1.(a) Tectonic outline of the Tibetan Plateau (modified after [23]); (b) simplified geological map of the eastern section of East Kunlun showing the distribution of Permian-Triassic granites (modified after [14]); and (c) simplified distribution map of the Hongshuichuan volcanics, showing the studied sample's location.

Figure 1 .
Figure 1.(a) Tectonic outline of the Tibetan Plateau (modified after [23]); (b) simplified geological map of the eastern section of East Kunlun showing the distribution of Permian-Triassic granites (modified after [14]); and (c) simplified distribution map of the Hongshuichuan volcanics, showing the studied sample's location.

Figure 5 .
Figure 5. (a) Chondrite-normalized REE patterns of zircons from the Hongshuichuan volcanics and (b-d) variations in selected trace elements for zircons from the Hongshuichuan volcanics.

Figure 5 .
Figure 5. (a) Chondrite-normalized REE patterns of zircons from the Hongshuichuan volcanics and (b-d) variations in selected trace elements for zircons from the Hongshuichuan volcanics.

Figure 10 .
Figure 10.Selected major and trace element diagrams of the Hongshuichuan volcanics showing fractionation trends of major and accessory mineral phases.(a) Sr versus Rb/Sr; (b) Ba versus Sr; (c) SiO2 versus Al2O3; and (d) Eu/Eu* versus Sr.

Figure 10 .
Figure 10.Selected major and trace element diagrams of the Hongshuichuan volcanics showing fractionation trends of major and accessory mineral phases.(a) Sr versus Rb/Sr; (b) Ba versus Sr; (c) SiO 2 versus Al 2 O 3 ; and (d) Eu/Eu* versus Sr.

Figure 10 .
Figure 10.Selected major and trace element diagrams of the Hongshuichuan volcanics showing fractionation trends of major and accessory mineral phases.(a) Sr versus Rb/Sr; (b) Ba versus Sr; (c) SiO2 versus Al2O3; and (d) Eu/Eu* versus Sr.
), showing affinity to arc-related magmatism.The felsic volcanics plot in the field of volcanic arc granitoids (VAG) in the Ta vs. Yb, Rb vs. (Y + Ta), Rb vs. (Y + Nb), and Nb vs. Y diagrams (Figure 12) further indicate their derivation in a subduction arc setting.The conclusion of the Early Triassic subduction arc setting is consistent with other geological evidence, including the Early Triassic arc-related adakitic quartz diorites and granodiorites (U-Pb age, ca.243 Ma; [14]) and the subduction-related I-type granites (ca.244-240 Ma; [44,76]).
), showing affinity to arc-related magmatism.The felsic volcanics plot in the field of volcanic arc granitoids (VAG) in the Ta vs. Yb, Rb vs. (Y + Ta), Rb vs. (Y + Nb), and Nb vs. Y diagrams (Figure 12) further indicate their derivation in a subduction arc setting.The conclusion of the Early Triassic subduction arc setting is consistent with other geological evidence, including the Early Triassic arc-related adakitic quartz diorites and granodiorites (U-Pb age, ca.243 Ma; [14]) and the subduction-related I-type granites (ca.244-240 Ma; [44,76]).

Figure 13 .
Figure 13.(a)Schematic cartoons showing the tectonic magmatic evolution of the Early Triassic igneous rocks in the EKOB; (b) Schematic diagram of the genesis of rhyolite and rhyolite porphyry; (c) Zircon εHf(t) and Mg # value distribution of the studied felsic volcanic rocks.

1 .
The studied rhyolite and rhyolite porphyry in the Hongshuichuan Formation from the East Kunlun orogenic belt, North Tibet, are synchronously formed in Early Triassic (ca.251-247 Ma). 2. Elemental and isotopic geochemistry studies indicate that the rhyolite porphyry and rhyolite have typical lower crust-derived characteristics, and they are generated by partial melting of Mesoproterozoic metagreywacke sources with a certain amount of mantle-derived mafic magmas involved in the origin of the high Mg # rhyolite.3. The studied ca.250 Ma felsic volcanic rocks were formed during the northward subduction of the East Kunlun Paleo-Tethyan oceanic slab.This study suggests that

Figure 13 .
Figure 13.(a) Schematic cartoons showing the tectonic magmatic evolution of the Early Triassic igneous rocks in the EKOB; (b) Schematic diagram of the genesis of rhyolite and rhyolite porphyry; (c) Zircon εHf(t) and Mg # value distribution of the studied felsic volcanic rocks.
and rhyolite porphyry in the Hongshuichuan Formation from the East Kunlun orogenic belt, North Tibet, are synchronously formed in Early Triassic (ca.251-247 Ma). 2.