Geochronology and Geochemistry of Late Paleozoic Volcanic Rocks and Their Relationship With Iron and Molybdenum Deposits in Xilekuduk Area, Northern Margin of Junggar

A large number of intermediate basic volcanic rocks and porphyry Cu-Mo deposits as well as volcanic-hosted magnetite deposit have been recently discovered in the Xilekuduk area. However, no reports concerning petrogenesis and age or their relationship with mineralization have been published to date. The purpose of this study was to make up for the absence of previous studies on Devonian volcanic activities in the area and to confirm the relationship between two stages of volcanic activities and mineralization so as to provide important theoretical basis for mineral exploration. Based on research results of zircon U-Pb geochronology and element geochemistry of volcanic rocks in the area, the ages of dacite, andesite, and stomatal andesite are considered as 375.2 ± 2.9 Ma, 386.5 ± 3.0 Ma, and 317.9 ± 2.9 Ma, respectively, corresponding to the Middle Devonian and Late Carboniferous Period. The Devonian volcanic rocks belong to the high-K calc-alkaline series and island arc volcanic rocks, which are enriched in LREE, strongly enriched in large ion lithophile elements Th, Rb, Ba, and K and relatively depleted in high-field strength elements (HFSEs) Nb, Ta, and Ti. The Carboniferous volcanic rocks are enriched in LREE, as well as the large ion lithophile elements Th, Rb, Ba, and K are strongly enriched, while depleted in the HFSEs Nb, Ta, and Ti; moreover, the contents of TiO2 and V are 0.94–0.97% and 178–183×10–6, which are higher than those of island arc basalts. According to mineralogical typomorphic characteristics and geochemical analysis, magnetite mineralization is divided into two phases. The early stratiform magnetite ore MT1 has magmatic characteristics, forming a volcanic rock type magnetite deposit related to Devonian volcanic eruption and sedimentation (375–386 Ma). The magnetite MT2 in the magnetite-quartz vein is considered as hydrothermal genesis, which is a metal mineral in the early metallogenic stage of Carboniferous (317.1 ± 2.9 Ma) volcanic eruption and subvolcanism, and may be related to porphyry molybdenum mineralization. Therefore, the volcanism and Fe-Cu-Mo mineralization in this area is characterized by multistage superimposed mineralization.


INTRODUCTION
The Xilekuduk area in Fuyun County, Xinjiang, is located in the northern margin of the Junggar Basin, in the collision suture zone of the Siberian Plate and the Junggar Plate (Sengör et al., 1993;Windley et al., 2007). A large number of Paleozoic strata were formed during the process of plate amalgamation and accretion in the Paleozoic, with crisscross faults and strong magmatic activities forming a geological unit with complex coupling of tectonics, magmatism, and mineralization, which gave birth to an important metallogenic belt of precious gold and nonferrous metals such as iron, copper, gold, nickel, and rare metals in China (Yang et al., 2009;Yang et al., 2012;Yang et al., 2014;Li et al., 2015;Wei et al., 2019).
The porphyry Cu-Mo deposit discovered in the Xilekuduk area in 2005 is of great significance in understanding mineralization and the metallogenic series in the Junggar area. The zircon U-Pb method was applied for accurately dating the granitoids of porphyry (copper) molybdenum deposit. The diagenetic age was concentrated at 340-315 Ma, which was considered as the product of tectonic magmatic evolution of plate collision-post collision (You et al., 2016;Long et al., 2009); the Re-Os isochron age of molybdenite was 327.1 ± 2.9 Ma, and the metallogenic age was determined as early stage of Middle-Late Carboniferous . It is considered that the source of ore-forming material and granite porphyry bear the characteristics of homologous magma (Long et al., 2014). The age data of a large number of andesitic, dacitic volcanic rocks, and pyroclastic rocks exposed in the ore district are yet to be reported. It is generally speculated that they belong to the volcanic sedimentary assemblage of the Lower Carboniferous Namingshui Formation (Long et al., 2009;Ding et al., 2011;Long et al., 2011;Long et al., 2014;You et al., 2016;Wei et al., 2019). Some new scientific issues, however, have also been raised in this study: 1) The volcanic sedimentary assemblage of intermediate acid volcanic rocks-pyroclastic rocks-carbonate rocks distributed in the ore district obviously differs from that of the marine continental interactive facies clastic sedimentary rocks of the Lower Carboniferous Namingshui Formation (the Jiangbastao Formation). Are these volcanic rocks all formed in Carboniferous as previously speculated? Are the two types of sedimentary assemblages the product of the same geological event? 2) In 2015, layered magnetite bodies that occurred at the interface between volcanic rocks and clastic rocks were successively discovered more than 100°m southwest of the porphyry copper-molybdenum deposit; their coupling relationship with volcanic rocks, however, remains unclear.
The above issues hinder our understanding of volcanism and the genesis of the deposit in this area and limit the breakthrough of mineral exploration. Therefore, on the basis of detailed field investigation, the author has carried out systematic geochronology and rock geochemistry of volcanic rocks, determined the formation age and genesis of volcanic rocks, and through mineralogical studies of two periods of magnetite, combined with regional data, discussed the volcanic magmatic events and mineralization in this area, so as to provide scientific guidance for mineral exploration in this area.

Regional Geological Background
The tectonic location of the study area is located in the Gabosar Island Arc Belt of the East Junggar Orogenic Belt (You et al., 2016) ( Figure 1A). The regional exposed strata are predominantly Paleozoic, Devonian, and Carboniferous, followed by Cenozoic tertiary and Quaternary, with local occurrence of Permian and Jurassic. Among these, the Devonian is the most widely distributed, which is a set of littoral and neritic clastic rock, volcanic lava, pyroclastic rock with carbonate rock, and terrigenous clastic rock. It is the main ore-bearing layer of iron and copper polymetallic in the area. Regional fold and fault structures are developed, mainly including the Yitieke aketas composite anticline, Sarbulak karatongke composite syncline, and Jiaposar composite anticline as well as the northwest trending Erqisi deep fault, Ulunguhe deep fault, and Kalashinger fault. Among these, the Kalashinger fault has been formed most recently as the main rock-controlling and orecontrolling fault structure in the region. The magmatic activity in the area is intense and various intrusive bodies are formed. The rock types are mainly granite. The diagenetic age is concentrated in the Early-Middle Devonian of Late Paleozoic, followed by the Early Carboniferous. The volcanic rocks mainly occur in the Devonian Tuorangekuduke, Beitashan, Yundukala, and Kashion formations, and Carboniferous Jiangbastao and Batamayineishan formations (Zhang et al., 2006;Li et al., 2015) ( Figure 1B).

Geological Characteristics of the Deposit
According to previous studies, the volcanic and pyroclastic rock strata exposed in the ore district are assigned to the Lower Carboniferous Nanmingshui Formation (Long et al., 2009;Ding et al., 2011;Long et al., 2011;Long et al., 2014;You et al., 2016;Wei et al., 2019). Two groups of volcanic rock-clastic rock formations were decomposed in this field investigation: one group was the combination of intermediatebasic volcanic lava-sedimentary rock, while the other was the combination of intermediate-basic volcanic lava-volcanic clastic rock-clastic sedimentary rock. The two were in a fault contact relationship ( Figure 1C, Figure 2). The intermediate basic volcanic lava is mainly distributed in the northeast of the ore district, with an area of about 1.1 km 2 and a lithology thickness of about 520 m. Intermediate basic volcanic lava-pyroclastic rocks-clastic sedimentary rocks are mainly distributed in the northwest of the ore district, with an area of 8.9 km 2 .
It can be divided into five volcanic rhythms with a thickness of 3018 m, showing an integrated contact relationship of volcanic eruption: 1) Overflow facies andesite-erupted andesite lithic tuff, about 57 m thick. 2) Eruptive facies stomatal and esite-falling accumulation facies volcanic sedimentary tuff, with a thickness of about 110 m. 3) Overflow facies andesite-erupted andesite lithic tuff-sedimentary facies tuffaceous sandstone, about 1540 m thick. 4) Overflow facies dacite-sedimentary facies tuffaceous sandstone, about 810 m thick. 5) Overflow facies andesitesedimentary facies carbonate rock, sandstone, with a thickness of about 510 m, magnetite ore bodies are distributed in the skarnized zone of contact between pyroclastic rocks and carbonate rock. On the basis of type of rock association, the volcanic rock and clastic rock association exposed in the ore district is not completely Carboniferous.
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 657083 3 diorite occur successively from the north to the south of the ore district. Fault structures are extremely developed, comprising three main groups of NNW-NS, NW-EW, and NEE trending faults. Among them, NNW-NS fault is the most developed in the central part of the ore district, which is the main rock-controlling and ore-controlling fault in the ore district. It belongs to the secondary fault of the NW-trending regional Erqis south fault (Ablimit, 2017) ( Figure 1C).
The Cu-Mo deposit is a concealed porphyry deposit located in the northern part of the area. It mainly occurs in the internal and external contact zone of granite porphyry, quartz diorite, and tuffaceous siltstone. The ore body is stratoid and lenticular, strikes 350°-360°, inclines to the east with dips of 25°-40°, and occurs in parallel with local pinch of outbranches. It is 100-300 m long, 90-380 m wide, and 0.84-30.37 m thick. The average grade of Mo ore body is 0.061-0.097%. The metal minerals are mainly molybdenite and pyrite, with a small amount of chalcopyrite and magnetite. The gangue minerals mainly include quartz, sericite, chlorite, and epidote. The iron ore body is located in the skarn zone at the lithologic interface between pyroclastic sedimentary rock and carbonate rock in the south of the ore district. The surface mineralization chiefly comprises a beaded or lenticular shape. The exposed ore body is 1 m thick with a Cu grade of 1.57% and TFe content of 22.6%; the magnetite ore body with a thickness of 3 m is delineated by deep drilling (Figure 3), and the average TFe content is 18.54%. The main mineralization types are skarn type copper (gold) ore, skarn type copperbearing magnetite, and skarn type copper-bearing pyrrhotite magnetite.

1) Characteristics of Volcanic Rock Samples
According to the field observation and the measured section, the volcanic lava in the southwest of the ore district includes dacite and andesite, and the volcanic lava in the northeast is mainly stomatal andesite. The specific petrographic characteristics are as follows: Northeast stomatal andesite: the weathered surface is gray, while the fresh surface is cyan gray, in porphyritic structure; the matrix is felsic and in fine interwoven structure and massive structure; the rock is composed of phenocrysts and the matrix ( Figure 4A). The phenocryst minerals are plagioclase (PL) and hornblende (HBL), accounting for 15-25% of the total rock mass; the matrix is microcrystalline plagioclase, hornblende, and biotite; plagioclase is arranged in a directional manner; microcrystalline quartz is in felsic interlaced structure, accounting for about 80-85% of the total rock, and the accessory mineral is magnetite ( Figures 4B,C).
Southwest dacite: the weathered surface is gray, while the fresh surface is grayish white, in porphyritic structure; the matrix is in felsic interlaced structure and massive structure; the rock is composed of phenocrysts and matrix ( Figure 4D). The phenocryst minerals are plagioclase (PL), biotite, and quartz (QTZ), accounting for 8-10% of the total rock mass; the matrix is microcrystalline plagioclase, quartz, and biotite. Plagioclase is directionally arranged with quartz in felsic and fine texture, accounting for about 90-92% of the total rock mass ( Figures 4E,F). The accessory minerals are magnetite, while the altered minerals are sericite and chlorite.
Southwest andesite: the weathered surface is grayish green while the fresh surface is greenish green ( Figure 4G). The rocks are composed of phenocrysts and matrix. The phenocryst minerals are plagioclase (PL) and hornblende (HBL), accounting for about 10-20% of the total rock mass; the matrix is microcrystalline plagioclase, hornblende, and quartz. Plagioclase is oriented and interwoven, accounting for 85-90% of the total rock mass ( Figures 4H,I). The accessory mineral is magnetite.

2) Characteristics of Magnetite Samples
There are two main occurrence types of magnetite in the study area: 1) MT1 type: stratiform magnetite ore occurring along the strata ( Figure 5A). The surrounding rocks are mainly composed of andesitic volcanic tuff and tuffaceous siltstone ( Figures 3A,B, 5A). Magnetite is mainly of euhedral and subhedral texture ( Figures 5B,C). Intergrowth granular pyrite can be seen in Some MT2 magnetite veins cut through the contact zone between MT1 magnetite and volcanic rock stratum, and are cut by late pyrite vein ( Figure 5H). MT1 magnetite is often wrapped by MT2 magnetite, which indicates that the formation of MT2 magnetite is later than that of MT1.

1) Zircon U-Pb Dating
Zircon U-Pb dating was carried out in the Chinese Academy of Sciences' Guiyang Institute of Geochemistry. The test instrument was the GeoLas2005 laser ablation system produced by Agilent, company of the United States, and Elan 6100DRC ICP-MS. The laser spot diameter was 30 μm, the laser pulse was 10 Hz, the energy was 32-36 mJ, and the depth of laser ablation sample was 20-40 μm. The international standard zircon 91,500 was used as the external reference material for zircon age determination. For detailed analysis steps and data processing methods, please refer to the relevant literature (Yuan et al., 2003). The isotopic ratio and element content of zircon were calculated by Gitter (ver 4.0, Macquarie University) program and Andersen Tom (Andersen, 2002). LA-ICP-MS Common Lead Correction (ver 3.15) software was used for correction of common lead, and Isopolt (ver 3.0) software was used for age calculation and concordance diagram drawing.

2) Analysis of Elements
The major rare earth and trace elements were determined in the State Key Laboratory of Institute of Geology and Geophysics of the Chinese Academy of Sciences. The major elements were Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 657083 6 determined with application of the melting XRF method (national standard GB/t14506.28-1993) on X-ray fluorescence spectrometer 3080 E. The contents of rare earth and trace elements were determined by double-focus high-resolution ICP-MS produced by Finnigan MAT.
The zircon U-Pb isotopic analysis data of volcanic rocks ( Table 1) show that they are more consistent within the error range in the values of 207 Pb/ 235 U and 206 Pb/ 238 U, and the data fall on and near the concordant curve, the 206 Pb/ 238 U age of andesite in the northeast area was between (311 ± 4.6) to (319 ± 4.6) Ma, and its weighted average was (317.9 ± 2.9) Ma, with the weighted mean variance (MSWD) being 0.13, which can categorize the corresponding diagenetic age to the Lower Carboniferous ( Figure 7A). The 206 Pb/ 238 U ages of the dacite in the southwestern area ranged from 373 ± 5.1 Ma to 393 ± 10.7 Ma, and its weighted average was 375.2 ± 2.9 Ma (MSWD 0.20) ( Figure 7B). The 206 Pb/ 238 U ages of the southern porous andesite were between 369 ± 9.2 Ma to 382 ± 9.6 Ma, and its weighted average was 386.5 ± 3.0 Ma (MSWD 0.30), with the diagenetic age corresponding to the Middle and Late Devonian ( Figure 7C).

1) Major Element Geochemistry
The results of major elements of volcanic rocks are shown in Table 2. The SiO 2 content of Carboniferous stomatal andesite was between 55.01 and 55.80%, and the SiO 2 content was relatively low. The Al 2 O 3 content was between 17.08 and 17.48%, the CaO content was between 5.86 and 6.11%, and the content of MgO was   16.62-17.89%; the CaO content 2.66-3.56%; the MgO content 1.18-1.51%; the Fe 2 O 3 content 0.84-1.04%; and the TiO 2 content 0.34-0.35%. The total alkali content is between 5.78 and 6.68%, and it is sodium-rich and potassium-poor (Na 2 O/K 2 O 3.77-5.29), Mg number is 42.50-53. 45 and 32.31-38.74, respectively. The consolidation index (SI), the felsic index (FL), and the ironmagnesium index (MF) increase compared with the nonamygdaloid ones, indicating a high degree of magmatic crystallization differentiation (Table 3).

2) Geochemistry of Rare Earth and Trace Elements
The results of rare earth elements (REEs) and trace elements are shown in Table 4. The weight abundance (w ( REE)) of the Carboniferous stomatal andesite samples and their values of europium anomaly (δ Eu) and the values of (La) N / (Yb) N are (95.54-116.07)×10 -6 , 1.14-1.16, and 5.93-8.35, respectively. The diagram of chondrite-normalized REE distribution curve ( Figure 9A) shows that the rocks are slightly enriched in LREE and have no obvious Eu anomaly, implying that there was no obvious fractional crystallization of  plagioclase during the magmatic evolution. The primitive mantle standard diagram of incompatible elements ( Figure 9C) shows that large ion lithophile elements (LILEs) are relatively enriched and the high-field strength elements (HFSEs) are relatively depleted. The w ( REE) of Devonian dacite samples varies slightly 108.16-130.96 × 10 -6 , d Eu is 0.97-1.16, and w (La) N /(Yb) N is 5.67-7.65. The w ( REE) of Devonian andesite is 108.26-130.02 × 10 -6 , d Eu is 0.84-1.09, and w (La) N /(Yb) N is 5.12-7.73. The normalized REE distribution curve of chondrite ( Figure 9B) shows right-leaning LREE enrichment type, and there is no Eu anomaly in the Devonian dacites and andesites. On the primitive mantle standard diagram of incompatible elements ( Figure 9D), the LILEs are relatively enriched and the LILEs are relatively depleted.

Diagenetic Age and Tectonic Environment
As predecessors have speculated, andesitic, dacitic volcanic rocks and pyroclastic rocks exposed in the ore district belong to the volcanic sedimentary assemblage of the Lower Carboniferous Namingshui Formation (Long et al., 2009;Ding et al., 2011;Long et al., 2011;Long et al., 2014;You et al., 2016;Wei et al., 2019). However, in this study, LA-ICP-MS zircon U-Pb data show that the diagenetic age of dacite is 375.2 ± 2.9 Ma and that of andesite is 386.5 ± 3.0 Ma, indicating that the volcanic rock was formed in the Middle Devonian and the basaltic andesite is 317.9 ± 2.9 Ma, corresponding to the early Late Carboniferous. According to the 1:250000 and 1:50000 regional geological survey data and previous studies, there is a set of Middle Devonian Beitashan Formation mainly composed of shallow and coastal volcanic rocks, pyroclastic rocks, and carbonate rocks, and a set of Late Carboniferous Jiangbastao Formation mainly composed of marine and continental facies clastic rocks Zhang et al., 2008). Based on this, the strata in this area can be divided into two sets: the dacite-pyroclastic rocks in the southeast belong to the Beitashan Formation of Middle Devonian, and the clastic sedimentary rocks with stomatal andesite volcanic rocks in the northwest belong to the Lower Carboniferous Jiangbastao Formation.
The Devonian volcanic rock assemblage consists of andesite and dacite. The content of TiO 2 varies from 0.34 to 0.35% and from 0.42 to 0.47%, which belongs to high-K calc-alkaline series. The REE standard distribution model belongs to the LREE enrichment type and the large ion lithophile elements Rb, Ba, and K are strongly enriched, and the high-field strength elements (HFSEs) Nb, Ta, and Ti are relatively depleted. The Zr/Nb ratio is 15.38-23.90, which is close to the Zr/Nb ratio range of MORB (10-60, Davisin, 1996), and the Nb/Ta ratio is 9.41-12.14, which are lower than those of the original mantle (−17), indicating that Nb and Ta are separated by subducting fluid metasomatism of the original mantle (Othman et al., 1989). The geochemical characteristics of the volcanic rocks are similar to those of the Beitashan Formation in the region and are typical island arc volcanic rocks (Deng et al., 2009). The Carboniferous volcanic rocks are mainly stomatal andesite the REE standard distribution pattern belongs to LREE enrichment type. There is no obvious negative Ce anomaly, showing obvious continental volcanic characteristics (Zhang et al., 2006). The large ion lithophile elements Th, Rb, Ba, and K are strongly enriched, and the HFSEs Nb, Ta, and Ti are relatively depleted, which are similar to island arc basalt. However, the variation of TiO 2 content is 0.94-0.97%, and that of V content is 178-183×10 -6 , which is higher than that of island arc basalt. The metamorphic basicity of  the Lower Carboniferous Jiangbastao Formation in the eastern part of the study area shows the characteristics of back arc basalts (Woodhead et al., 1993;Zhou et al., 2005). The cause of this phenomenon may be linked to the transitional environment and complex formation process; consequently, it is inferred that the diagenetic environment belongs to the early stage of plate collision. On the Nb-Zr-Y tectonic environment discrimination diagram, Devonian and Carboniferous volcanic rock samples fall mainly in the volcanic arc basalt and the intraplate tholeiitic basalt area ( Figure 10A), and all samples fall in island arc volcanic rocks on the Hf-Th-Ta diagram, showing the products of island arc environment ( Figure 10B). This indicates that the diagenesis has certain transitional characteristics, and this may be related to the complex plate subduction and collision process.

Relationship Between Volcanism and Iron and Molybdenum Mineralization
Magnetite is widely distributed in various types of deposits, such as the skarn type deposit, the Kiruna type deposit, the iocg type deposit, the magmatic type deposit, the BIF type deposit, and the porphyry type deposit (Dupuis and Beaudoin, 2011;Huberty et al., 2012;Hu et al., 2014). The content and variation range of major elements, trace elements, and rare earth elements of magnetite with different genesis differ: this can indicate different formation temperature, oxygen fugacity, sulfur fugacity, ore-forming material source, and tectonic environment, and so on (Pearce and Gale., 1977;Nadoll et al., 2012). Therefore, the mineral chemical characteristics of magnetite have certain typomorphic significance for distinguishing different types of ore deposits. For example, magnetite in skarn-type iron ore contains lower V and Ti, higher Ca, Al, and Mn contents than that in magmatic iron ore, and their REE characteristics are also obviously different (Frietsch and Perdahl, 1995;Dupuis and Beaudoin, 2011;Hu et al., 2014;Nadoll et al., 2014). The iron-bearing minerals in the Xilekuduk area are mainly magnetite. Through microscopic mineral research, it is found that MT1 magnetite is generally wrapped and metasomatized by MT2, and the contact boundary between them is clear, indicating that there are at least two stages of mineralization in this area. In recent years, some scholars have also identified two phases of magnetite in adjacent areas. For example, the Qiaoxiahala and Laoshankou iron-copper-gold  Boynton (1984) and primitive mantle data from Sun and MCDonough (1989)) Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 657083 deposits in the northern margin of Junggar, Xinjiang, were divided into two phases of magnetite, and it was noted that "dissolution reprecipitation" may have occurred in the mineralization process (Shang et al., 2017). The EPMA data of two typical magnetite samples (to be published in another article) show that MT1 magnetite has higher Ti and V and lower Ca, Al, Cr than MT2 magnetite. In the genetic diagram of TiO 2 -Al 2 O 3 -MgO ( Figure 11A), most of the early crystallized magnetites fall in the acid-alkaline magmatic magnetite region. The magnetite crystallized in the late stage is dispersed, and some of them fall in the area of sedimentary metamorphism-contact metasomatism.
In the genetic diagram of TiO 2 -Al 2 O 3 -MgO+ MnO genetic diagram ( Figure 11B), magnetite crystallized in the early stage mainly falls in the volcanic rock type and the magmatic melt separated titanomagnetite type area. The magnetite crystallized in the late stage mainly occurs in contact metasomatic and skarn type areas. In the genetic diagram of TiO 2 -Al 2 O 3 -MgO ( Figure 11C), most of the early crystallized magnetites fall in the granite area (acidic magmatic pegmatite). The magnetite crystallized in the late stage mainly occurs in hydrothermal and calc skarn type areas. Therefore, the early crystallized magnetite MT1 is mostly of magmatic origin, while the late crystallized magnetite MT2 shows the characteristics of hydrothermal origin. This further confirms the existence of at least two stages of mineralization of magnetite. According to the banded occurrence characteristics of MT1 magnetite and the occurrence along with volcanic rocks in the Middle Devonian, it is inferred that the formation of MT1 magnetite is related to the   (Xu and Shao, 1979): I: accessory mineral type; II: magmatic fused titanomagnetite type; III: volcanic rock type; IV: contact metasomatic type; V: skarn type; VI: sedimentary metamorphic type. (C) (Lin, 1982):I: granite area (acid magmatic pegmatite); II: basalt area (tholeiite, etc.); III: gabbro area (gabbro peridotite, monzonite, and anorthosite-accessory minerals and iron ore); IV: peridotite area (peridotite, peridotite, pyroxene, and other accessory minerals and iron ore); V1: amphibolite area (including clinopyroxenite); V2: diorite area; VI: kimberlite area; VII: hydrothermal type and calc skarn type (mainly plutonic hydrothermal type above dotted line, hydrothermal type and calc skarn type below); VIII: hydrothermal type and magnesium skarn type (plutonic hydrothermal type, part of hydrothermal metasomatism; magnesium skarn type); IX: sedimentary metamorphism and hydrothermal superimposition type; X: carbonate rock area (×1 is related to ultrabasic rocks, X2 is related to wall rock replacement); XI: transitional area.
eruption and deposition of Devonian volcanoes (375-386 Ma). However, MT2 magnetite occurs in the form of magnetite quartz vein. Microscopically, it can be observed that MT2 magnetite is metasomatized by molybdenite and pyrite, which has hydrothermal characteristics. This may be similar to the "dissolution-reprecipitation" of magnetite in the surrounding the Qiaoxiahala and Laoshankou iron-copper-gold deposits. The relationship between magnetite and molybdenite in this stage is relatively close. Combined with the porphyry molybdenum metallogenic event in the ore district, it is speculated that MT2 magnetite may be related to porphyry molybdenum mineralization. According to regional data, the study area experienced plate subduction, collision, and post-collision regional tectonic events in the late Paleozoic (Wang et al., 2006;Cai et al., 2011;Xiao and Santosh, 2014), accompanied by strong volcanism, such as volcanic eruption, eruption, and subvolcanism (Liu, 2001). According to the two stages of Late Paleozoic volcanic rock formation in the Xilekuduk area, the coupling relationship between two stages of volcanism and iron and molybdenum mineralization is sorted out: In the early Devonian, the Junggar plate began to subduction northward and was in the tectonic environment of the plate margin. Volcanic rocks transited from tholeiitic (low-potassium) series to calc-alkaline series and then to high-potassium series of dacite, indicating that the island arc experienced the evolution process from early stage to mature stage. There are a large number of andesite dacite assemblages in the Xilekuduk area. The ages of these volcanic rock assemblages are between 375 and 385 Ma, showing typical island arc magmatic characteristics. During the intermittent period of volcanic eruption, the magmatic fluid rich in iron and copper first crystallized and precipitated between the layers of volcanic rocks, or interposed in the fissures of volcanic rocks, forming magnetite ore bodies mainly in strip and vein (Wilkinson, 2013), and a large number of banded and vein magnetite (MT1) were developed. With the crystallization differentiation of magma and the precipitation of copper element, the high background value of copper in volcanic rocks of the Beitashan Formation is formed. The Cu content of basalt is 95 × 10 −6 -25 × 10 -6 , which is about 4 times of the Clark value of crust and that of andesite is 20 × 10 -6 , which is roughly equivalent to the Clark value of crust, which provides material source for the formation of later copper deposits (Yang et al., 2010). In the Qiaoxiahara iron-copper-gold deposit on the northern side of the study area, iron ore body mainly occurs in the basalt-andesite volcanic rock of the Beitashan Formation of the Middle Devonian, which is derived from the early intermediate-basic volcanism (Ying, 2007). Therefore, the basalt andesite formation in the Xilekuduk area is an important iron-bearing area.
In the Early Carboniferous, the Junggar plate continued to subduction to form the Devonian-Early Carboniferous back arc ocean basin (Li, 2009). Volcanism in this stage was mainly subvolcanic facies, which was composed of intermediate-intermediate acid hypabyssal and ultrahypabyssal intrusive facies igneous rocks with different depths. The diagenetic age was concentrated in 340-318 Ma. Basaltic andesite porphyrite xenoliths are distributed in the granodiorite in the Xilekuduk area, suggesting that the source of crust mantle mixing of acid rock mass: the upper mantle material upwelling of asthenosphere induces the melting of lower crust material, and the mantle-derived magma rises and mixes with the crust source material, resulting in largescale magma mixing (Wang et al., 2015), forming ore-bearing magma. Subsequently, the ore-bearing magma intruded and emplaced along the extensional structure, metasomatized the Devonian overflow andesite-sedimentary carbonate formation, forming a skarn belt. At the same time, a large amount of magmatic segregation formed a small amount of magnetite quartz veins (MT2), superimposed or wrapped early stratiform magnetite ore (MT1). With the gradual evolution of fluid, a large amount of Cu-Mo precipitates to form the Xilekuduk porphyry Cu-Mo deposit.
Volcanism and Fe-Cu-Mo mineralization in this area, therefore, bear the characteristics of multistage superimposed mineralization. The andesite-dacite formation in the Middle Devonian is closely related to the formation of early-stage magnetite, and the Early Carboniferous homologous heterogeneous volcanic rock formation, such as basaltic andesite-diorite-granite porphyry, is closely related to the porphyry copper mineralization and late-stage magnetite.

CONCLUSION
1) The diagenetic age of dacite is 375.2 ± 2.9 Ma and that of andesite is 386.5 ± 3.0 Ma, which indicates that the volcanic rock was formed in the Middle Devonian, and the diagenetic age of stomatal andesite is 317.9 ± 2.9 Ma, it was formed in the Late Carboniferous. 2) The Devonian volcanic rock assemblage is composed of andesite and dacite. The content of TiO 2 varies from 0.34 to 0.35% and 0.42-0.47%, respectively. It belongs to the high-K calc-alkaline series, enriched in LREE, enriched in large ion lithophile elements Th, Rb, Ba, and K, and relatively depleted in HFSEs Nb, Ta, and Ti. The Carboniferous volcanic rocks are mainly stomatal andesite. The REE standard distribution pattern belongs to LREE enrichment type. The large ion lithophile elements Rb, Ba, and K are strongly enriched. The HFSEs Nb, Ta, and Ti are relatively depleted. The contents of TiO 2 and V are 0.94-0.97% and 178-183×10 -6 , which are higher than those of island arc basalts. They are characterized by transitional diagenetic environment and are close to the stage of plate collision. 3) There are two periods of magnetite in the ore district. The stratiform magnetite ore MT1 produced along the early strata is mostly of magmatic origin. It is inferred that its formation is related to the eruption and deposition of Devonian volcano (375-386 Ma), while the magnetite MT2 of late magnetite quartz vein shows the characteristics of hydrothermal mineral, indicating the metal minerals precipitation in the early stage of the mineralization of the Carboniferous (317.1 ± 2.9 ma) porphyry copper molybdenum deposit. Therefore, the volcanism and Fe-Cu-Mo mineralization in this area have the characteristics of multistage superimposed mineralization.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.