Geochronology and Geochemistry of the Karadaban Bimodal Volcanic Rocks in the Altyn Area, Xinjiang: Implications for the Tectonic Evolution of the Altyn Ocean

Detailed geochronological, geochemical, and Sr-Nd-Hf isotopic data are presented for early Paleozoic volcanic rocks in the Karadaban area from the northern Altyn region, NW China, with the aim to constrain their petrogenesis and tectonic implications. The Karadaban volcanic rocks show a bimodal distribution in composition, with rhyolite and basalt. The LA-ICP-MS zircon U-Pb age indicates that the volcanic rocks were erupted at 512Ma. The ma ﬁ c rocks are calc-alkaline, enriched in light rare earth elements (LREE) and large-ion lithophile elements (LILE; Ba and U) and depleted in high- ﬁ eld strength elements (HFSE; Nb and Ta). These features together with their depleted isotopic signature (initial 87 Sr/ 86 Sr = 0 70413 – 0 70817 , ε Nd t = 2 7 to 3.7) suggest that they were likely derived from a depleted mantle source but mixed with crustal components while upwelling. The felsic rocks show an A-type a ﬃ nity, with high alkalis and Rb/Sr and Ga/Al ratios; enriched in LILE (e.g., Rb, K, Th, U, and REE) and depleted in Ba, Sr, Nb, P, and Ti; and with fractionated REE patterns with strong negative Eu anomalies. The combination of the decoupling of ε Nd t values ( − 2.5 to − 6.3) and ε Hf t values (+5.5 to +14.7) in the setting of subduction indicates that the felsic rocks were generated by partial melting of the juvenile crustal as a result of magma upwelling. The geochemical and Sr-Nd-Hf isotopic characteristics, coupled with regional geology, indicate that the formation of the Karadaban bimodal volcanic rocks involves an extensional regime associated with a subduction-related environment. The rifting of the back arc in response to the retreat of the subducting northern Altyn oceanic lithosphere may account for the Karadaban bimodal volcanic rocks.


Introduction
Bimodal volcanism typically characterizes an extensional environment, which can occur within various tectonic settings, including continental rifts [1,2], within-plate extensional settings [3][4][5], intraoceanic islands [6], ocean island arcs [7,8], incipient back-arc depressions [9], mature island/active continental margins [7][8][9], and back arcs [9][10][11][12]. In each of these modern tectonic environments, the volcanic activity may give rise to specific features, such as lithological assemblage, geochemical signature, and type of associated ore deposits [12,13]. Moreover, of all these various associations of mafic and felsic rocks, it is perhaps the A-type rhyolites with basalts that are the most enigmatic but may also prove to be the most instructive. Therefore, the recognition of such distinctive bimodal igneous suites in ancient orogenic belts can not only provide diagnostic geodynamic tracers for constraining the tectonic evolution of these belts but also document important information on crustal growth by magmatic underplating in continental interiors.
The Altyn Tagh Belt (ATB) marks the northern margin of the Qinghai-Tibet Plateau [14], lying between the Tarim block to the north, the Qaidam block, the Qilian orogen, and the Kunlun orogenic belt to the south [15,16]. As one of the most important tectonic zones in northwestern China, the Altyn Tagh carved and dominated several tectonic units to match the eastern extrusion of the Tibetan Plateau [17][18][19][20][21]. In the last decade, the North Altyn Tagh region was recognized as an ophiolite [14,[22][23][24], with HP/UHP metamorphic rocks [25][26][27] and igneous rocks [28][29][30][31][32]. Previous studies show that three stages can be divided about the collisional orogeny in the Altyn area: the oceanic-continent collision stage [31,32], the oceanic-continent subduction stage [33][34][35], and the postcollisional extension stage [36,37]. At present, the mechanisms of the tectonic setting of the last two stages of magmatism are almost clear due to the zircon U-Pb ages from igneous rocks and metamorphic ages from HP/UHP metamorphic rocks [14,31,34,35,37,38]. However, the oceanic-continental collision stage of magmatism in the north Altyn region is still unclear, which is subject to few systematic studies. In addition, the oceaniccontinental subduction setting is always accompanied by the mixture of crustal and mantle magmatism, but most of the studies are focused on granite only. Recently, there was a set of bimodal volcanic rocks discovered in the NE of the North Altyn Tagh, which is helpful for studying the magmatic events and determining the timing of oceaniccontinental subduction. In this paper, we present U-Pb zircon dating and geochemical and Sr-Nd-Hf isotopic composition for Karadaban bimodal volcanic rocks from the north Altyn region. Our aims are (1) to constrain petrogenesis and magma sources of bimodal volcanic rocks in the subduction setting and (2) to discuss the early Paleozoic evolution of the Altyn Tagh region.

Geological Background
The Altyn Mountain is located at the juncture of the Tarim Basin, the Qaidam Basin, and the Qilian orogen in Western China [19,[39][40][41][42][43], which plays a significant role in accommodating convergence between India and Eurasia [20,44,45]. The Altyn Tagh fault system is the first-order tectonic unit that marks the northern margin of the Qinghai-Tibet Plateau and influences the geological evolution of surrounding areas. Pin and Paquette proposed that the sedimentation in the southern Tarim Basin was associated with this large-scale (>1600 km long) intracontinental fault system [8] (Figure 1(a)).
The tectonic units of the ATB have been divided into four different schemes (Figure 1(b)) [19,[47][48][49]: the North Altyn Tagh Archean complex (NAAC), the North Altyn Tagh subduction-collision complex (NASC), the Milanhe-Jinyanshan block (MJB), and the South Altyn Tagh subduction-collision complex (SATSC) (Figures 1(b) and 1(c)) [19,25,50]. In the NAAC region, the HP argillaceous rocks, glaucophane schist, and eclogite occurred belonging the HP-LT metamorphic belt. Previous studies showed that the phengite in the eclogite and the paragonite in the glaucophane schist have a yield of 39 Ar- 40 Ar isochron ages of 513 ± 5 Ma and 497 ± 10 Ma, respectively [51]. These HP-LT metamorphic rocks are an important indicator of the palaeosubduction zone [52]. Moreover, there are many mafic-ultramafic plutons that are developed in the NASC, in which the ophiolite belt is dated at 524-437 Ma, with the features of the ocean island basalt, indicating that the oceanic crust may have developed probably in the Early Cambrian [14,53,54]. Researches about the age of felsic rocks in NASC show that they mainly ranged from 500 Ma to 420 Ma, which is associated with the setting of subduction to collision to postcollision [55][56][57]. The NASC is about 20 km wide, which extends from Hongliugou eastward to Lapeiquan for 150 km, and is controlled by the NEE-Altyn faults. In the Altyn region, the Paleozoic formation of the Zhuo'abulake and Simierbulake volcanic sedimentary rock series is located at the center of the Karadaban fault, which is fault-contacted to the Cenozoic and Archean strata. In addition, there is a suite of Karadaban bimodal volcanic rocks that developed in the north of NASC, which occurred in the Zhuo'abulake Formation ( Figure 2). Magmatism peaked at the early Paleozoic and stretched along the NE faults, which produced a series of granodioritic and granitic rocks. Chen et al. reported that the age of Kaladawan felsic volcanic rocks (zircon SHRIMP U-Pb method) was 485-477 Ma, and the characteristic of geochemistry indicated a setting of active continental margin [34]. Cui et al. (2010) verified the geochemistry of basic volcanic rocks in the Karadaban region. These rocks showed the features of tholeiitic basalt and conjectured an ocean island arc setting, which developed at the hanging wall of subduction in the early Paleozoic [58]. However, the timing of the subduction was still disputable.

Petrology and Sampling
In this study, 12 samples within the Karadaban bimodal volcanic rocks were collected for geochemical analysis, 8 of them were used for Sr-Nd isotopic analysis, and analyses of Hf isotopes were made on zircon grains from samples KB-26 and KB-81 that were previously dated by laser ablation multiple collector inductively coupled plasma mass spectrometry (LA-ICP-MS). Lithology of these samples are basalt and rhyolite ( Figure 3).
The rhyolite samples are gray to light yellow in color, massive, and porphyritic and rhyolitic in texture, which are interbedded with the basalt (Figure 3(a)). The felsic samples contain 10-15% plagioclase and quartz phenocrysts in a groundmass consisting of glassy cryptocrystalline material and feldspar and quartz (crystal sizes < 0 1 mm). The quartz is subhedral-anhedral, some of which has been rounded by dissolution (Figures 3(b) and 3(c)). The plagioclase is subhedral and lath-shaped (up to 0.5 mm) and hypidiomorphic in texture. The groundmass is typically a fine-grained assemblage of plagioclase, quartz, and muscovite. Muscovite is present as microscopic flakes (generally <0.1 mm), preferentially oriented aggregates, and disseminated single crystals with a preferred orientation (Figure 3(d)). The dark greyishgreen basalts occur as massive layers with a porphyritic texture. The main phenocrysts are plagioclase and pyroxene (Figures 3(e) and 3(f)). The plagioclase is euhedral, tabular, and up to 0.3-0.4 mm, and the pyroxene occurs as thin 2 Geofluids tabular crystals up to 0.5-0.6 mm in size. The groundmass includes fine-grained plagioclase, pyroxene, and magnetite, which is characterized by the fasciculate variolitic structure of fibroid pyroxene and plagioclase (Figure 3(g)).

Analytical Methods
The Karadaban igneous samples were collected from adits and outcrops ( Figure 2

Zircon U-Pb
Dating. Cathodoluminescence (CL) images were obtained by using a scanning electron microscope (SEM, Leo 1450VP, Germany) at the Chinese Academy of Geological Sciences before in situ U-Pb and Hf isotopes analyses. Zircon LA-ICP-MS U-Pb dating was performed using a laser ablation inductively coupled plasma mass spectrometer (LA-ICP-MS) at the MLR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences (CAGS), Beijing. The laser spot diameter and frequency were 40 μm and 10 Hz, respectively. The Harvard zircon 91500 was used as an external standard for zircon U-Th-Pb analyses and NIST610 as an external standard to calculate the contents of U, Th, and Pb. The 207 Pb/ 206 Pb and 206 Pb/ 238 U ratios were calculated using the GLITTER program [60], and common Pb was corrected using the method of Andersen [61]. Age calculations and concordia plots were done using Isoplot (ver. 3.0) [62]. The details of the analytical techniques were documented by Hou et al. [63].

Lu-Hf Isotope Analyses.
In situ zircon Hf isotopic analyses were carried out at the Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, using a Neptune MC-ICP-MS with an ArF excimer laser ablation system. During analyses, the spot sizes of 32 and 63 μm and a laser repetition rate of 10 Hz with 100 mJ were used. During analyses, the 176 Hf/ 177 Hf and 176 Lu/ 177 Hf ratios of standard zircon (91500) were 0 282294 ± 15 (2σn, n = 20) and 0.00031, similar to the commonly accepted 176 Hf/ 177 Hf ratio of 0 282284 ± 3 (1σ) measured using the solution method [64,65].

Sr-Nd Isotopic Analyses.
Samples for isotopic analysis were dissolved in Teflon bombs after being spiked with 84 Sr, 87 Rb, 150 Nd, and 147 Sm tracers prior to HF+HNO 3 (with a ratio of 2 : 1) dissolution. Rubidium, Sr, Sm, and Nd were separated using conventional ion exchange procedures and measured using multicollector mass spectrometer (IsoProbe-T) at the Analytical Laboratory of BRIUG [66]. Procedural blanks were <100 pg for Sm and Nd and <500 pg for Rb and Sr. 143 Nd/ 144 Nd were corrected for mass fractionation by normalization to 146 Nd/ 144 Nd = 0 7219, and 87 Sr/ 86 Sr ratios were normalized to 86 Sr/ 88 Sr = 0 1194. Typical within-run precision (2σ) for Sr and Nd was estimated to be ±0.000015. The measured values for the La Jolla and BCR-1 Nd standards and the NBS-607 Sr standard were 143 Nd/ 144 Nd = 0 511853 ± 7 (2σ) and 0 512604 ± 7 (2σ), respectively, and the measured value during the period of data acquisition is 87 Sr/ 86 Sr = 1 20042 ± 2 (2σ).

Major and Trace
Elements. Twelve whole-rock major and trace element data are presented in Table 1. The loss on ignition (LOI) for the rocks was in the generally range from 0.74 to 3.68 wt.% (one sample at 6.87 wt.%), which probably reflects the presence of alteration-related hydration. Considering the potential modification from the concentrations of mobile components [67,68], the effects of alteration and the mobility of the major and trace elements were evaluated prior to the geochemical and    4 Geofluids petrological investigation. Only immobile elements were used to classify rock samples and discuss their tectonic settings and petrogenesis.
The SiO 2 content of the volcanic rocks lie within a range from 44.74 to 76.14 wt.% and with a compositional gap at 46-68 wt.% SiO 2 , which shows the bimodal feature (Table 1).    6 Geofluids Note: major elements were obtained by X-ray fluorescence spectrometer (XRF) at the Analytical Laboratory of BRIUG, Beijing, using fused lithium-tetraborate glass pellets. The analytical precision as determined on the Chinese national standard GSR-3 was generally around 1-5%. Trace elements were analyzed using a PerkinElmer ELAN-DRC-e ICP-MS at the Analytical Laboratory of BRIUG, Beijing. The powdered samples (50 mg) were dissolved in high-pressure Teflon bombs using a HF+HNO 3 mixture for 48 h, at~195°C. Rh was used as an internal standard to monitor signal drift during counting. The international standards GBPG-1 and OU-6 and the Chinese national standard GSR-3 were used for analytical quality control, with the precision generally better than 5%.

Geofluids
Combining with the Zr/TiO 2 -Nb/Y plot, the samples were identified as basalt and rhyolite (Figure 4(a)). Mafic rocks, on the TiO 2 -Zr/P 2 O 5 diagram (Figure 4 Table 2. Zircons selected from the medium-to coarsegrained (sample KB-81) and medium-to fine-grained (sample KB-26) rhyolites are mostly display euhedral crystals with clear oscillatory zoning, euhedral, and up to 50 μm in length, with length to width ratios of nearly 2 : 1. Their CL images commonly show oscillatory zoning ( Figure 6), which indicated a typical feature of magmatic zircons [73].  Figure 5: Chondrite-normalized REE element pattern and primitive-mantle normalized trace element spider grams for the bimodal volcanic rocks. The normalized values of chondrite are from [70], and primitive mantle values are from [71]. Data of the A-type granite are from [72], Lapeiquan granite are from [24], Hongliugou granite are from [57], Kaladawan rhyolite are from [34], and Highland 4337 granodiorite are from [32].  10 Geofluids  According to the relationship between rhyolite and basalt in the outcrop, the two rocks are the contemporary volcanic rocks, which belong to the early Paleozoic.

In Situ Zircon Hf
Isotopes. Zircon Lu-Hf isotopic data are listed in Table 3. Their initial εHf t values were calculated using their U-Pb zircon ages.   Table 4. The whole-rock initial 87 Sr/ 86 Sr ratios and ε Nd t values have been calculated at t = 512 Ma on the basis of the zircon U-Pb ages of this study ( Figure 6). The initial 87     1.79 14 Geofluids the northern Altyn region, and the age of volcanic rocks was 512 4 ± 2 8 Ma and 512 2 ± 3 Ma, which belongs to early Cambrian. It is an important stage for researchers to study the evolution of the northern Altyn region. Previous studies show that the magmatic activity in the Altyn area generated a series of magma evolution from basic to felsic rocks, which mainly occurred from 520~480 Ma [29,30,56,76,77]. The early Paleozoic tectonic belt of Hongliugou-Lapeiquan was recognized, and the results showed that the basic volcanic rocks are mainly with the features of MORB and OIB [78]. In this belt, the LA-ICP-MS zircon U-Pb age of the eutectic gabbro is 513 ± 3 Ma, which can be classified as a subduction zone ophiolite, and mainly developed in the oceanic subduction belt [79]. Liu et al. (1998) reported that the whole-rock isochron age of the pillow basalt in Hongliuquan area is 524 ± 44 Ma [80], which probably represented the forming age of ocean island basalt. Gao et al. (2011) obtained a LA-ICP-MS zircon U-Pb age of 513 ± 2 Ma from the plagiogranite, which is intrusive in the ophiolite mélange.  obtained the SHRIMP zircon U-Pb age of 482~477 Ma from the felsic volcanic rocks in the Kaladawan area, which indicated that the characteristic of magmatism continued from 520 Ma to 477 Ma. In addition, there are high pressure-ultrahigh pressure (HP/LP) metamorphic zones in the tectonic belt, and the 39 Ar- 40 Ar isochron age of 513 ± 5 Ma from muscovite in the blueschist showed that the HP/LP was effected by the oceanic subduction. In this paper, we verified that the age of felsic volcanic rocks in Kaladaban is 512 Ma, corresponding with the setting of the north Altyn region. Combining with previous studies, we inferred that there was an ocean basin that existed at the north Altyn region in the late Neoproterozoic and the subduction activities probably continued until the late Ordovician. Since the oceanic crust subduction and collision, a series of subduction-related volcanic rocks were developed in the north Altyn region [34]. The magmatic activity in the northern Altyn region is multistage and probably successive from Neo-Proterozoic to Late Paleozoic.

Petrogenesis of Bimodal Volcanic Rocks.
Bimodal volcanic rocks occur within a variety of tectonic settings, such as continental rifts, oceanic islands, zones of continental break-up, back-arc basins, orogenic belts, intraoceanic island arcs, and mature island arcs/active margins [9,81,82]. Recent studies have proposed two models of the genesis of bimodal volcanic rocks, depending on the cognate of the rhyolites and basalts (e.g., [83]). Since the mafic and felsic magmas originated from different sources, the trace element contents and the Sr and Nd isotopic compositions should differ between the basalts and the rhyolites [84,85], thereby providing a reliable indicator for the genesis of bimodal volcanic rocks.
Samples from the northern Altyn region have experienced less alteration. The loss of ignition (LOI) values for the mafic and felsic volcanic rocks in the study area are 1.92-6.87 and 0.74-3.68 wt.%, respectively. This alteration is a result of low-temperature metamorphism related to intense deformation and fluid infiltration, whose effects are principally evident in the various types of abundance of the mobile elements such as K 2 O, Rb, Ba, and Cs. However, the Fe, Al, Ca, Mg, REE, and HFSE elements are considered to be relatively immobile during the low-temperature alteration [65,[85][86][87]. Accordingly, the following discussion of the magma geochemical composition will focus on the immobile element and their ratios, as well as isotopic compositions.
6.2.1. Petrogenesis of the Mafic Rocks. The low silica content (SiO 2 = 44 74 -46 65 wt %) and relatively high concentrations of T Fe 2 O 3 (11.49-17.12 wt.%), Cr (175-215 ppm), and Ni (94-118 ppm) content of the mafic rocks suggest that a mantle component played a prominent role in their genesis [88]. The slight lack of negative Eu anomalies on the chondrite-normalized REE diagrams indicates an absence of plagioclase fractionation during the process of magmatic evolution (Figures 5(a) and 5(b)). Meanwhile, the basalt displays a calc-alkaline nature and is enriched in LREEs and is also depleted in HFSEs (e.g., Nb and Ta). These features are similar to the subduction-related magmas [89]. However, the basaltic samples show MgO, Cr, Ni, V, and Co are negatively correlated with SiO 2 , indicating that the basaltic magma is not primitive and possibly experienced some fractional crystallization of pyroxene. This is consistent with the presence of pyroxene as the dominant phenocryst in the basalt. The Karadaban basalt is enriched in Th (4.08-11.90 ppm) and positively correlated the plots of Nb/La and Nb/Th (Figure 7(a)). So we conjectured the involvement of Th-enriched component in its petrogenesis, which is due to the crustal contamination [37]. The continental crust is typically depleted in Nb and Ta [90], and the upper continental crust is enriched in La and Th, while the lower continental crust is not always enriched in Th [91]. Combining with the (Th/Ta) PM vs. (La/Ta) PM plot, the basalt samples define an UCC trend (Figure 7(b); [92]), implying that the contaminator is the highly evolved upper continental crust. It is noted that mafic magmas derived from the asthenosphere typically have La/Nb ratios of <1.5 and La/Ta ratios of <22, whereas those rocks from lithosphere have La/Nb ratios > 1 5 and La/Ta ratios > 22 [93,94]. High La/Nb ratios (2.15-3.36), La/Ta ratios (36.57-53.49), and positive εNd t values (+2.7 to +3.7) indicate that the basalt was derived from a depleted mantle source but mixed with crustal components while upwelling. Mantlederived magmas generally assimilate crustal materials during their ascent or their storage in a crustal magma chamber [95]. The Nb/U values (2. 16-9.13) are lower than those of the source mantle (34), MORB (47 ± 11), or OIB (52 ± 15), further indicating that the basalt was affected by the crustal contamination [96]. Moreover, the mafic rocks from the Karadaban area are Nb-Ta depleted but enriched in Zr, Hf, and Ti ( Figure 5(b)), yielding patterns similar to those expected for crustal material which suggests that these magmas assimilated crustal material during ascent [97].
Excluding the garnet residue in the source, the Nb/Yb and Ta/Yb ratios are effective in constraining the relative depletion of the mantle source when the degree of partial melting is even lower [70,98]. Meanwhile, the Th/Yb ratio is a sensitive indicator of the inclusion of subduction zone materials in magma [87,91]. The Th/Yb vs. Nb/Yb diagram ( Figure 8) indicates that the basalt samples in the present study developed in an active continental margin, mainly close to the boundary between the E-MORB and OIB arrays. In summary, the geochemical and Sr-Nd isotope data demonstrate that the basalt of the Karadaban bimodal volcanic rocks was derived from the depleted mantle source and experienced fractional crystallization, which was dominated by the fractionation of pyroxene and crustal contamination possibly by upper crustal materials. These compositional features suggest that the samples in the present study originated from a source that lacked garnet but contained plagioclase in the residue, which indicates partial melting under low-pressure conditions [100,101]. The absence of Ba, Nb, Ta, Sr, P, and Ti, as shown in the primitive-mantle-normalized trace-element spider diagram ( Figure 5(d)), indicates similar trace element composition of A-type granites [101], and the felsic rocks generally display similar patterns to the basalts. However, the Nb/Ta and Th/Ta ratios of the felsic rocks are low, suggesting that the rhyolite is contaminated by upwelling basaltic magma. Furthermore, the Karadaban felsic samples fall into the A-type granite area on a Ga/Al-Zr plot (Figure 9(a)), and the Yb/Ta and Y/Nb ratios of the felsic samples are higher than those of OIB and close to IAB (Figure 9(b)). This indicated that the felsic rocks are A-type granites and exhibit similarities to the crust between OIB and IAB. Based on the subdivision of A-type granites [102], the Karadaban felsic samples are type A2 ( Figure 10); however, the other granites (e.g., Highland 4337 and Hongliugou) show I-type features. It coincides where all data fall into the A-type granite and (b) the Yb/Ta vs. Y/Nb diagram (after [102]) exhibiting similarities to average crustal ratios between OIB (ocean island basalt) and IAB (island-arc basalt). Data of the Hongliugou granite are from [57], Highland 4337 granodiorite are from [32], and south Kaladawan granite are from [46].  Figure 10: Plots of the volcanic rocks from the North Altyn Tagh region in Nb-Y-3Ga and Nb-Y-Ce (after [102]). A1: anorogenic granite; A2: postorogenic granite. Data of the Hongliugou granite are from [57], Highland 4337 granodiorite are from [32], and south Kaladawan granite are from [46]. 17 Geofluids with the amount of island arc granitoid developed in the north of Karadaban volcanic rocks [32,57]. Hence, we suggest the Karadaban felsic rocks probably generated in an extensional environment.
Two types of process have been proposed for the origin of the silicic end-member in the bimodal magmatic suite, including (1) extensive fractional crystallization from a common mantle-derived magma parental to the mafic endmember, coupled with crustal contamination [103], and (2) crustal anatexis caused by mantle-derived mafic magma with distinct isotopic compositions [104][105][106]. In the case of the felsic rocks, the distinct difference documented by the whole-rock ε Nd t and zircon εHf t values between the basalt and felsic rocks indicated that the basalt could not have been produced by fractional crystallization of the coeval basalt. However, the negative ε Nd t (−6.3 to −2.5) and positive εHf t (4.5 to 14.7) of Karadaban felsic rocks showed Nd-Hf isotopic decoupling. The decoupling was probably effected by fluid/melt metasomatism in the subduction zone; the solubility of Nd was higher than that of Hf and led to the higher Nd/Hf ratios in the fluid/melt derived from the subducted plate [107,108]. Coincidentally, from Cambrian to Ordovician, the whole north Altyn region was in a subduction setting, and a great deal of crustal growth occurred through arc-related magmatism (Table 5). For the Nd-Hf isotopic decoupling, the Hf isotope composition of zircon is more realistic than Nd. The positive zircon εHf t (4.5 to 14.7) and T DM2 ages of 551 to 1186 Ma suggest a juvenile source rather than an old basement rocks. In addition, the wide range of εHf t indicated the mixture of the end-member during the magma evolution [109]. When the subduction metasomatized mantle-derived felsic crustal materials were rejuvenated, the newly formed magmas would inherit the decoupled Nd-Hf characteristics such as Karadaban rhyolite in accordance with A 2 -type granites representing magmas derived from the continental crust which has been through a cycle of island-arc magmatism (Figures 11 and 12; [102]). As presented above, the continued  18 Geofluids extension of the crust and repeated underplating of subsequent mafic magma partial melting of the juvenile crust to generate felsic melts will ascend and erupt on the surface forming the Karadaban rhyolite. Based on the discussion above, we conclude that the felsic volcanic rocks produced mainly by partial melting of the juvenile crustal and the origin melt source itself likely resulted from mantle-derived magmatic upwelling in an extensional setting, probably near the subduction-related active continental margin.
6.3. Tectonic Implications. The north Altyn region is a part of the Proto-Tethys tectonic domain. The Cambrian-Ordovician basic rocks and calc-alkaline granite in this region are interpreted to be island arc volcanic rocks [34,35]. SHRIMP zircon ages of 506 ± 2 Ma for a granodiorite and 481 ± 5 Ma for arc granite indicate that the Cambrian to early Ordovician age is related to arc magmatism [30,31]. In addition, bimodal volcanic rocks indicate an extensional setting in which they form, and several types of the extensional setting can account for the generation of bimodal volcanic suite, such as within-plate (continental rifting), passive continental marginal rifting, and incipient back-arc or continental active margins [9]. Hence, we propose that the Karadaban bimodal volcanism could be regarded as a significant indication of arc-related extensional tectonics.
On the one hand, according to the feature of trace element, the Karadaban volcanic suite displays some geochemical features of subduction-related melts, such as enrichment of LILEs but depletion in HFSEs and Sr-Nd-Hf isotope compositions of magmas. This suggests that the mixing of juvenile crust and mantle-derived magma, as well as the Karadaban bimodal volcanic suite, was formed in a subduction-related environment. The conclusion is further supported by a series of tectonomagmatic discrimination diagrams. Based on the Hf-Th-Nb and Zr-Ti-Y ternary diagrams, the Karadaban rocks are classified as volcanic arc calc-alkaline basalts (Figure 13), which support the hypothesis that magmatic traces existed in the study area [99]. The Karadaban basalts exhibit both arc-like and within-plate basalt affinities, implying that they are probably formed within a back-arc tectonic setting. Generally, the Kaladaban bimodal volcanic suite in the north Altyn region is likely an analogue of the Dzungaria Ocean in the southern Central Asian Orogenic Belt during the late Paleozoic [70]. The Karadaban bimodal volcanic rocks resemble those of the back-arc side of the Dzungaria arc, which are believed to result from the initiation opening of the Dzungaria Ocean [70,112]. On the other hand, regional geological data suggest that the final closure of the Altyn Ocean did not occur until the Early Ordovician and even until the later Ordovician in the north Altyn region [24,35]. This demonstrates that the Kaladaban bimodal volcanism was produced in a back-arc environment (Lapeiquan back-arc basin), rather than in the postorogenic or intracontinental rifting setting.
The Altyn region is of great significance to realize the evolution of the Tibetan Plateau and the tectonic pattern of Western China. It developed as a unified back-arc basin system between the northern Altyn Complex and MJB before the formation of the Altyn Tagh, which is the northern boundary of the Proto-Tethys in China [113]. Three stages can be identified by previous studies: north Altyn Oceanic rifting (~750 Ma) [47,114], north Altyn oceanic subduction (520~460 Ma) [34,76,77,115], and oceanic consumption followed by continental collision (440~420 Ma) [29,109,116,117]. This paper mainly focuses on the study of the bimodal volcanic rocks developed during the subduction of the northern Altyn Ocean. Combined with previous studies, the subduction is presumed to have occurred around 520~460 Ma (Table 5) [19,77,117] from the plentiful research of regional magmatic sequences  to the mafic to intermediate to felsic rocks [31,46,109]. Based on our new data and previous studies, we propose a model for the tectonic evolutions of the north Altyn region during the Early Cambrian. This model emphasizes a southward subduction system of the oceanic lithosphere (Altyn Ocean) between the NASC and MJB. Beneath the MJB, its consolidated NASC margin during the early Paleozoic period and an arc and back-arc system developed, resembling the present western Pacific continental margin, appearing during the early Cambrian [24,35,46]. The Karadaban bimodal volcanic rocks formed in the incipient stage of the Cambrian back-arc basin/extension. The system might have lasted into the Late Ordovician, during which the back-arc basin gradually became mature and finally closed during the Early Silurian. In conclusion, the north Altyn oceanic crust subducted towards the MJB and transmitted the latter into the active continental margin. The magma, which derived from the upwelling and mixing caused by subduction, eventually emplaced into the crust as Kaladaban bimodal volcanic rocks.

Conclusions
The following conclusions can be drawn from this study.
(1) New LA-ICP-MS zircon U-Pb dating indicates that the newly discovered bimodal volcanic rocks in the northern Altyn Tagh area formed at ca. 512 Ma (2) Whole-rock geochemical and Sr-Nd-Hf isotopic compositions suggest that the mafic rocks were derived mainly from a depleted mantle source but mixed with crustal components. The felsic rocks were probably generated by partial melting of the juvenile crustal as a result of basic magma upwelling (3) The generation of the Karadaban bimodal volcanic rocks was not only related to an active continental margin setting but also associated with the subduction-related back-arc extension environment

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.