Geochronological and Geochemical Constraints on the Petrogenesis of Permian Dolerite Dyke Swarms in the Beishan Orogenic Belt, NW China

Extensive Early Permian mafic-ultramafic intrusions, doleritic dykes, and basalts crop out within the Beishan area, southern Central Asian Orogenic Belt (CAOB). We present new geochronological and geochemical data for Gubaoquan dolerite dyke swarms in the Beishan orogenic belt. Zircon U-Pb Dating of the Gubaoquan dykes indicates that they were emplaced during the Early Permian (280.7 ± 4 Ma), that was coeval with Yinaoxia and Podong mafic dykes in Beishan area. The dykes are characterized by low Mg# (47–84) in the clinopyroxene crystals, and the content of whole-rock Fe2O3 (t), MgO, and alkali (Na2O + K2O) range from 12.5–17.4, 4.06–5.51, and 2.8–4.4 wt.%, respectively. The samples from the Gubaoquan dykes have high and variable Ba/La (5.93–14.2) and Ba/Nb (15.0–37.3) ratios but low Th/Yb (0.17–0.24) ratios. The rocks show slightly enrichments in LREE, HFSE, Th, and Hf, and depletion in Nb and Ta. The εNd (t = 280 Ma) values and initial 87Sr/86Sr ratios of the Gubaoquan dykes show variations ranging from 6.4 to 6.8 and 0.706240 to 0.707546, respectively. These data suggest that the parental magmas for the Gubaoquan dykes were probably derived from partially decompressed melting of upwelling depleted asthenosphere mantle that was metasomatized previously by subducted fluids.


INTRODUCTION
The Central Asian Orogenic Belt (CAOB) is the important Phanerozoic accretionary orogenic belt in the world ( Figure 1A; Sengör et al., 1993;Jahn et al., 2000). It was formed by the amalgamation of accretionary complexes, magmatic arcs, arc-related basins, ophiolites, seamounts, and microcontinents during the long-lived subduction of the Paleo-Asian Ocean Windley et al., 2007;Xiao et al., 2008Xiao et al., , 2011Xiao et al., , 2019. The Beishan Orogenic Belt is located in the southern margin of the CAOB, and experienced a complex process of tectono-magmatic evolution during the Paleozoic. Therefore, it has become a key part to study the geological evolutionary history of the CAOB. The Early Permian was a crucial period in the evolutionary history of the CAOB, marked by the closure of the Paleo-Asian Ocean, the transition from subductionaccretion to post-collisional tectonism and the occurrence of extensive mantle-derived magmatism (Wartes et al., 2002;Charvet et al., 2007). The Early Permian (290-270 Ma) magmatism and associated mineralization have been considered as the post-subduction byproducts after the closure of the Paleo-Asian Ocean Yuan et al., 2010;Song et al., 2011Song et al., , 2013Zhang et al., 2011). The petrogenesis of the Early Permian mafic-ultramafic complexes, doleritic dykes, and basalts in the Beishan orogenic belt, and is still in controversial. Zhou et al. (2004) believe that the Permian magmatic rocks in the Tarim area and the adjacent Beishan area together constitute a Large Igneous Province (LIP), which is a product of mantle plume activities (Zhou et al., 2004;Qin et al., 2011;Su et al., 2011Su et al., , 2012. Zhang et al. (2011) suggest that the Permian mantle-derived coeval in the Beishan area are both products of the post-subduction extension Li et al., 2013Li et al., , 2020Zheng et al., 2014). Ao et al. (2010) believe that Permian mafic-ultramafic complexes, such as Pobei and Hongshishan, were formed by subduction process of the Paleo-Asian Ocean, representing the Late Paleozoic island arc settings in the Beishan area (Ao et al., 2010).
Previous studies have mostly focused on the investigation of the Early Permian mafic-ultramafic complex and basalt in the Beishan area, but the petrogenesis of the doleritic dykes is poorly constrained and lack of comparison between different regions, due to the limited geochemical investigations of mafic dykes in this area. As a specific product of tectonomagmatism (Hoek and Seitz, 1995), dolerite dyke swarms can be formed in the marginal areas of the mantle plume-related LIP, the back-arc extension environment and the post-subduction extension settings. The doleritic dykes are important signs of tectono-magmatic evolutionary events and can maximally retain information of mantle source (Halls, 1982;Liu et al., 2013). Thus, they are widely used in reconstructing the aggregation, extension, and breakup of ancient plates (Zhou et al., 2000). In addition, mafic dykes can serve as magma conduits, allowing the transportation of mafic magma from the mantle to the crust and potentially preserving primitive magma compositions that provide crucial information on the nature of the underlying mantle (Adams et al., 2005). Accordingly, we present new geochronological and geochemical data for Gubaoquan mafic dyke swarms in the Beishan area ( Figure 1B), and investigate their petrogenesis and the dynamic background This study aims to provide a basis for further understanding the geotectonic background in Beishan area during Early Permian.

Geological Background
The Beishan orogenic belt is located in the southernmost margin of CAOB, adjacent to the Central Tianshan orogenic belt in the north, and the Tarim-North China Craton in the south ( Figure 1B). The boundary between CAOB and the North China Craton is the South Tianshan-Liuyuan-Solonker suture zone ( Figure 1A), which is considered to represent the final suture zone on the southern margin of CAOB. High-pressure to ultrahigh-pressure metamorphic rock belts related to the subduction of the Paleo-Asian Ocean are distributed across the suture zone, including Akeyazi blueschist and eclogite in the southwestern Tianshan (Ar-Ar age: 315-319 Ma) (Gao and Klemd, 2003;Klemd, 2003;Klemd et al., 2015;Li et al., 2015), Gubaoquan eclogite in the south of Liuyuan (age of peak metamorphism: ∼465 Ma; Liu et al., 2012;Saktura et al., 2017), and the Wenduermiao metamorphic rocks of blueschist facies in Solonker (Ar-Ar age: ∼450 Ma) (Gao et al., 1990).
Available zircon U-Pb chronological data has shown that there were at least two phases of mafic-ultramafic magmatism in the Beishan area: one from the Late Ordovician to Devonian (360-450 Ma) and the other from the Late Carboniferous to Early Permian (260-305 Ma). The widely distributed ophiolite belts formed during the first stage in the Beishan area includes high-Mg diorite, Nb-rich basalt, normal-mid-oceanic ridge basalt (N-MORB), enriched-MORB (E-MORB), and island arc tholeiite, indicating that the Beishan area may have experienced oceanic crust subduction from the Late Ordovician to Early Devonian (Mao et al., 2012;Zheng et al., 2018). In the later phase, mafic-ultramafic magmatism was frequent, and several maficultramafic complexes and dolerite dyke swarms were developed. These dolerite dykes are extensively distributed across the entire Beishan area and were once associated with the Tarim mantle plume (Xiao et al., 2006).

Permian Mafic Dyke Swarm
The dolerite dyke swarms in the Beishan area are concentrated in the Pobei area in the west, Gubaoquan area in the middle, and Yinwaxia area in the east. Moreover, there are scattered outcrops in other regions ( Figure 1B). The dykes are generally parallel, striking toward NW (310-355 • ) and NE (10-80 • ), and are nearly vertical (dips of 80 • -90 • ; Figure 2E). The dyke swarms predominantly intrude the Silurian-Devonian (380-442 Ma) granodiorite and quartz diorite, as well as the Precambrian to Permian strata (Figures 2A,C,E; Mao, 2008;Mao et al., 2010) in the Liuyuan area. The contact boundaries between the dykes and country rocks and strata are conspicuous (Figures 2A,C,E). Baked edges can be seen between individual dykes and the granodiorite.
Gubaoquan-Hongliuyuan Fault. These dykes generally strike between NNE-NE and NW-SE, and are various in size, with widths ranging from 0.5 to 6 m. Dykes extend from hundreds of meters to thousands of meters, and are nearly vertical ( Figure 2E). The Gubaoquan doleritic dykes mainly intrude the Ordovician-Devonian granodiorite, and some of them intrude the metamorphic strata of the Early Paleozoic Huaniushan Group ( Figure 1B). The dykes are interspersed in the rock mass and cut through the strata, mostly in apophyses or stocks. The doleritic dykes exposed in the strata generally occur parallel to each other, with greenish-black to dark gray color, and have undergone weak chlorite and sericite alteration ( Figure 2F). The dolerite dykes in the granodiorite rock mass are well preserved, mostly in dark green color, and the boundary between the dyke and country rocks is clear and flat ( Figure 2E). These mafic dykes contain clinopyroxene (∼ 15 vol.%) and plagioclase (∼ 10 vol.%) phenocrysts ( Figures 2B,D). Matrix minerals are composed of clinopyroxene and plagioclase with minor hornblende and Fe-Ti oxides ( Figure 2D). The Gubaoquan doleritic dykes are porphyritic and fresh, with the ophitic textures ( Figures 2B,D,F) and comprise the low proportion of Fe-Ti oxides (2 vol.% ∼3 vol.%).

ANALYTICAL METHODS
Zircon crystals were separated from a fresh doleritic dyke sample (∼25 kg, GBQ-04, 41 • 00 43.09 N, 95 • 00 21.52 E) using the conventional heavy liquid and magnetic separation techniques in a commercial laboratory at Langfang, Hebei Province, China. Representative zircon crystals those were euhedral, colorless, transparent and underformed were mounted in an epoxy resin disk and then polished to analysis. Cathodoluminescence (CL) and back-scatted electron (BSE) images were obtained at the Beijing SHRIMP Center, Institute of Geology, Chinese Academy of Geological Sciences. CL and BSE images were used together to help select the zircon grains that have oscillatory or simple zoning patterns indicative of magmatic crystallization for U-Pb dating. Zircon U-Pb dating for mafic dykes from the Gubaoquan areas was carried out using LA-ICP-MS at the Key Laboratory of Western China's Mineral Resources and Geological Engineering, Ministry of Education. Detailed descriptions of the instrumentation and analytical procedures can be found in Luan et al. (2019). A Plešovice zircon ( 206 Pb/ 238 U age = 337 Ma, Sláma et al., 2008) and Qinghu zircon ( 206 Pb/ 238 U age = 159 Ma, Li et al., 2009) standard were used as external standard for correcting mass discrimination and isotope fractionation, and Zircon 91500 ( 207 Pb/ 206 Pb age = 1065 Ma, Wiedenbeck, 1995) was analyzed alternately with unknown zircons. Concordia diagrams and weighted mean calculations were carried out using the Isoplot 3.0 program (Ludwig, 2003). Common Pb correction was made using the measured 204 Pb and the values of Stacey and Kramers (1975).
Major element concentrations for clinopyroxene and plagioclase were determined by a JEOL JXA-8100 electron microprobe at Chang'an University, China. The analytical conditions were: 15-kV accelerating voltage, 20-nA beam current, 1-µm beam size, and a 20-s peak counting time. International mineral standards were used for calibration. Analytical reproducibility was within ± 2% relative.
Whole-rock major and trace element compositions were determined at Chang'an University, China. The major element concentrations were measured by X-ray fluorescence (XRF) spectrometer (Shimadzu 1800) on fused glass disks, and Chinese national rock standards GBW07112 (gabbro) was used for calibration. Analytical precision for major elements was less than 3%. Whole-rock trace elements analyses were performed using a American Thermoelectric X-7 ICP-MS instrument, at Chang'an University. Sample powder (40mg) were acid-digested with mixed HNO 3 and HF in steel-bomb coated Teflon beakers for 48 h in order to ensure complete dissolution of refractory minerals. Take the USGS basalt standard BCR-2 as external standard sample to calibrate the trace elemental concentrations of the measured samples, and analytical precision was generally better than 5%.

Zircon U-Pb Age
Zircon U-Pb data for the Gubaoquan mafic dykes are listed in Table 1. Sample GBQ-04 was collected from a mafic dyke in the Gubaoquan area that was intruded into Silurian-Devonian granodiorite (424-380 Ma; Mao et al., 2010) and Precambrian orthogneisses (Figures 2A,C; Yang et al., 2006;Mao, 2008;Saktura et al., 2017). The zircons separated from this sample are colorless and transparent in subhedral granular or columnar shapes, with a length of 50-100 µm and a length to width ratio of 1:1-2:1. In cathodoluminescence (CL) images, the zircon grains have flat crystal planes, and the majority of these zircons have well-developed oscillatory zoning, indicating a magmatic origin (Figure 3a). The majority of the older zircon grains have recrystallized blurred zoning and are stubby and subhedral, whereas the three youngest zircons are elongated grain and have more developed oscillatory zoning (Figure 3a). The fourteen measured zircon 206 Pb/ 238 U age values vary in a wide range (Figure 3b). Four zircons yielded Precambrian ages (2044 to 897 Ma), and seven zircons yielded Silurian ages (441 to 425 Ma). In addition, the three youngest ages yield a weighted average age of 280.7 ± 4 Ma (Figure 3b). Accordingly, we interpret the age of 280.7 ± 4 Ma as the emplacement age of the Bubaoquan doleritic dykes, and all older grains as inherited zircons captured from the wall rocks.

Mineral Chemistry
The compositions of clinopyroxene and plagioclase in the Gubaoquan doleritic dykes are given as Supplementary , with a small amount of diopside and endiopside (Figure 4). Clinopyroxene from Gubaoquan mafic dykes have the variable Mg # values (47-84), and mainly concentrated in 66 to 84.

Major and Trace Elements
Whole-rock major and trace element compositions of the mafic dykes from the Gubaoquan are given as Supplementary Table 2. The content of whole-rock SiO 2 and Fe 2 O 3 T from Gubaoquan dykes ranges from 47.7 to 50.1 wt.% and from 12.5 to 17.4 wt.%, respectively, with relatively lower concentrations of MgO (4.06-5.51 wt.%) and total alkali (Na 2 O + K 2 O = 2.83-4.37 wt.%). The Mg# values (magnesium index of rock) of the Gubaoquan dolerite are between 41.5 and 57. These dykes are relatively rich  in iron and poor in magnesium thereby belonging to iron-rich mafic rock group, and are also sub-alkaline and basaltic in a TAS diagram (Figure 5; Le Bas et al., 1986), that is similar to the major element characteristics of dolerites in the Yinwaxia and Podong areas (Zhang et al., 2015;Xue et al., 2016), whereas the type of dolerite formed in the same period in the Tarim Bachu area is dominately alkaline basalt (Yu, 2009). The chondrite-normalized rare-earth element (REE) patterns and primitive mantle-normalized immobile incompatible element patterns are illustrated in Figures 6A,B. The Gubaoquan dykes show slightly depleted to enriched LREE patterns ([La/Yb] N = 0.93-1.32) and weakly negative Eu anomalies (δEu = 0.87-0.99; Figure 6A), which may be caused by the  Sun, 1995). As show in Figure 6B, the samples are characterized by pronounced Nb-Ta negative anomalies and slightly light REE enrivhments relative to heavy REE, that is similar to the trace element distribution characteristics of the Mesozoic island arc basalts. Compared with the regional contemporaneous dykes and basalt, the Liuyuan basalt is relatively more enriched in LREE and HFSE, while the Gubaoquan dykes have slightly higher trace element content than that of doleritic dykes in the Yinwaxia and Podong areas. However, they have similar distribution patterns of rare earth elements and trace elements which are significantly different from the geochemical characteristics of Bachu mafic dykes of the same period in the Tarim area.
primary magmatic zircons are difficult to crystalize in dolerite.
Considering the fact that mafic dykes are able to capture zircons from their wall rocks during ascent, the majority of these zircons can be interpreted as xenocrysts captured from wall rocks. Previous data show that Precambrian metamorphic basement were extensively developed in Beishan area (China University of Geosciences Press, 1997;Niu et al., 2020). The doleritic dykes from Gubaoquan area mainly intruded the Silurian-Devonian granodiorite (424-380 Ma; Mao et al., 2010) and Precambrian orthogneisses (Figures 2A,C; Yang et al., 2006;Mao, 2008;Saktura et al., 2017), respectively. It is inferred that the formation age of the Gubaoquan doleritic dykes should be later than the Devonian. This means that the three youngest zircons (279.1, 280.5, and 282.8 Ma) are most likely to crystallize from the mafic magma, and their weighted average 206 Pb /238 U age of 280.7 ± 4 Ma can be taken as the best estimation of the emplacement age of the Gubaoquan dykes. Whereas, the older age groups (2044-897 Ma and 441-425 Ma) represent the ages of zircons captured from the basement strata (Precambrian orthogneisses) and country rocks (granodiorite) during the ascent of magma. In addition, According to the statistics of the age of doleritic dykes (Table 3) in the Beishan area, the emplacement age range of the dolerite dyke swarms is 271.2 ± 2.9-313.6 ± 3.3 Ma, and concentrated in 285-270 Ma, similar to that of the mafic dykes in Podong, Huaniushan, and Yinwaxia area. Moreover, the age of doleritic dykes obtained by Zhang et al. (2015) in the Liuyuan Gubaoquan area is 282 Ma, indicating that there was extensive mafic magmatism from the Late Carboniferous to Early Permian in the Beishan area.

Magma Source and Crustal Contamination
The results of the whole-rock Sr-Nd isotopic simulation reveal that the magma of the Gubaoquan doleritic dykes experienced contamination with upper-crust materials during the emplacement. In the simulation, the two-component of the magma isotopic composition are taken from theprimitive mantle that is the intersection point between the sample and the mantle trend line [εNd (t) = 9.5, ( 87 Sr/ 86 Sr) i = 0.703], and the upper crust [Rudnick and Gao, 2014; εNd (t) = 6, ( 87 Sr/ 86 Sr) i = 0.72]. The magma evolved from the depleted mantle has Sr and Nd concentrations of 90 ppm and 7.3 ppm, respectively. The Sr and Nd concentrations of upper-crust contaminant end members are assumed to be 320 ppm and 27 ppm, respectively. These values lie within the mean value ranges of oceanic basalts and upper crust, respectively (Sun and McDonough, 1989;Rudnick and Gao, 2014). The simulation calculation results show that the Gubaoquan mafic dykes have a crustal contamination degree of 5∼10% (Figure 7), which is close to the Liuyuan basalts and higher than the doleritic dykes in Yinwaxia and Pobei areas. In addition, the presence of xenocrystic zircons from the Precambrian strata and granodiorite provided the evidence for crustal contamination during the magma emplacement.
Mafic magma generally originates from the asthenosphere or lithospheric mantle (Fiona et al., 2017). The Fe/Mn ratios of mafic rocks can be used to identify their parental magma composition (Liu et al., 2008). In general, MORB have the lower Fe/Mn ratios (55-58; Zhang et al., 2011) than the plume-derived basalts such as Hawaiian OIBs (about 65-71; Humayun et al., 2004) and continental flood basalts (∼ 63.7; Wang et al., 2011). Olivine and clinopyroxene fractionation will decrease and increase Fe/Mn ratios, respectively (Liu et al., 2008). The relatively uniform Fe/Mn ratios (54.1-59.7) of the Gubaoquan doleritic dykes is similar to the MORB, and suggest that the crystallization of olivine and pyroxene has little overall impact on the Fe/Mn ratio of the magma. This reveals that the Fe/Mn ratios of the Gubaoquan dykes reflect the composition of their mantle source and imply a MORB-type parental magma. In addition, the εNd (t) value of the Gubaoquan doleritic dykes varies between 6.4 and 6.8, showing the characteristics of a depleted mantle. All of the mafic dykes plot above the MORB-OIB array on a Th/Yb versus Nb/Yb diagram, that are similar to Liuyuan basalts, Podong and Yinaoxia doleritic dykes (Figure 8), further suggesting the involvement of a subduction-related component (Pearce, 2008). However, the Bachu doleritic dykes and Keping basalts plot close to the OIB in the MORB-OIB array, with obvious different sources. Subduction components mainly include siliceous melts and slab-derived fluids, both can impose distinctive geochemical signatures on mantle sources (Pearce, 2005). Highly mobile elements (e.g., Rb, Ba, Cs, and U) are tend to partition into aqueous fluids, whereas Th and LREE are strongly partitioned into sediment-derived siliceous melts (Plank and Wade, 2005). The high and variable Ba/La (5.93-14.2) and Ba/Nb (15.0-37.3) ratios but low Th/Yb (0.17-0.24) ratios of the Gubaoquan doleritic dykes strongly suggest they were derived from a mantle source that was metasomatized by slab-derived fluids (Woodhead   (2008); The CAB, OIB, MORB fields are from Li et al. (2015); Date sources: Yinaoxia dykes (Zhang et al., 2015), Podong dykes (Xue et al., 2016), Tarim Bachu dolerite dykes (Yu, 2009), and Tarim Keping basalt (Jiang et al., 2004a), and Beishan metamorphic rock (Xie et al., 2012(Xie et al., ). et al., 2001Plank and Wade, 2005). Therefore, based on the above results, we propose that the parental magma for the Gubaoquan mafic dykes was derived from the depleted Asthenospheric mantle that metasomatized previously by the subducted fluids, and then contaminated by crustal materials.

Geodynamic Background of Magmatism
The parent magma of these dykes underwent assimilation and contamination with upper-crust materials during the emplacement process, which may have caused of the negative Nb-Ta anomalies in the Gubaoquan dolerite, because the continental crust is severely depleted in Nb and Ta (Rudnick and Gao, 2003). However, there are other possibilities to cause negative Nb-Ta anomalies. Island arc basaltic magma and the mantle-derived magma that was previously metasomatized by a subductionrelated component both have the characteristics of negative Nb-Ta anomalies (Rudnick and Gao, 2003;Pearce et al., 2013). In order to further understand the genesis of negative Nb-Ta anomalies in the Gubaoquan doleritic dykes, we use the whole-rock Nb/Yb and Th/Yb ratios to evaluate the explanations (Pearce, 2008). As shown in Figure 8, the trace element ratios represent the parental magma composition, if we assume that the negative Nb-Ta anomalies in dykes are entirely caused by the contamination with crustal, up to ∼15% crustal contamination is required (Figure 8), indicating that the assimilation and contamination with more crustal materials is required to reach the current level of Nb-Ta depletion. As mentioned above, the Sr-Nd isotope compositions of the mafic dykes show significantly lower degrees of crustal contamination (< 10%, Figure 7). The contradiction rules out the possibility of crustal contamination as a major cause for the negative Nb-Ta anomalies in the dykes. And the Nb/Yb and Th/Yb ratios of Gubaoquan dykes are within the field of Cenozoic continental arc basalts. The characteristics of trace element and Sr-Nd isotopes indicate that the negative Nb-Ta anomalies probably reflect the primary signature of a mantle source. The parent magma of the Gubaoquan doleritic dykes had a certain degree of Nb-Ta depletion. This show that the parent magma of the Gubaoquan dykes probably was the basaltic magma in island arc setting or the mantle-derived magma whose source had been metasomatized by subduction fluids before its formation.
The Beishan orogenic belt was formed by the suturing of blocks located in different places in the Paleo-Asian Ocean during the Paleozoic (Zheng et al., 2014(Zheng et al., , 2016Chen et al., 2017;Wang et al., 2018). Available geological data have shown that the ophiolite belt in the Beishan area was formed between 330 Ma and 516 Ma (Zhang and Guo, 2008;Li et al., 2012;Zheng et al., 2013;Wang et al., 2014), which represent the upper time limit of the oceanic crust subduction between different blocks. The metamorphic age of the high-temperature and high-pressure metamorphic rock belt distributed across the South Tianshan-Liuyuan-Solonker suture zone on the southernmost margin of CAOB was 465-315 Ma (Gao and Klemd, 2003;Klemd, 2003;Liu et al., 2012;Li et al., 2015;Saktura et al., 2017), whick represents the end of oceanic crust subduction, i.e., the subduction of Paleo-Asian oceanic crust ended before the Late Carboniferous (∼315 Ma), and subsequently, the stage of arc-continent collision and post-subduction extension was started.
The Gubaoquan dolerite was formed in the Early Permian (∼280 Ma). After the subduction of the Paleo-Asian Ocean, the Beishan area has entered the stage of intracontinental evolution. Moreover, considering that there was no arc-related magmatism in the Huaniushan arc or Shibanshan arc since Early Permian. In other words, the significant negative Nb-Ta anomalies in the dykes probably were not formed by the basaltic magma in the island arc setting. Instead, they are the combined results of metasomatism of the mantle-derived magma by fluids from ancient subducting slab and crustal contamination. The trace rare earth elements and Sr-Nd isotopic characteristics of the Gubaoquan doleritic dykes indicate that its primary magma originated from the depleted asthenosphere mantle, and the magma experienced 5%∼10% of assimilation and contamination with upper-crust materials during the emplacement. The petrological and geochemical characteristics of the Gubaoquan doleritic dykes also indicate that they formed under conditions of low water fugacity that are different from the conditions of arc-related basaltic rocks (Feig et al., 2006;Metcalf and Shervais, 2008). Therefore, during the post-subduction extension stage (∼280 Ma), lithospheric delamination system was most likely responsible for the genesis of the Permian Gubaoquan mafic dykes in the Beishan area (Figure 9).
Consequently, we suggest that the Gubaoquan doleritic dykes were formed during the tectonic period after the plate subduction. The lithospheric delamination in an extensional environment led to the partial melting of upwelling asthenospheic mantle and the overlying metasomatized mantle, which consequently generated the mafic magma of dolerite. In addition, on a larger scale in Beishan, the formation age of the Gubaoquan doleritic dykes was similar to that of the widely distributed doleritic dykes in the Beishan area, and the Gubaoquan mafic dykes have similar geochemical characteristics to the dykes in the Yinwaxia in the east and Podong in the west. This indicates that the Late Carboniferous to Early Permian doleritic dykes in the Beishan orogenic belt may be controlled by a unified continental geodynamic background, i.e., the dykes are related to the post-subduction extensional setting.

CONCLUSION
Zircon U-Pb dating indicates that the Gubaoquan doleritic dykes in the Beishan orogenic belt were emplaced at ca. 280 Ma, that was the same time as the formation of the dolerite dyke FIGURE 9 | Schematic diagrams for the Permian tectonic evolution of the southern Beishan area. Lithospheric delamination leading to extensional tectonism and partial melting of upwelling asthenospheric mantle and overlying lithospheric mantle. swarms extensively distribution in Beishan and the Tarim LIP. The geochemical characteristics of these dykes are obviously different from the mafic dykes related to mantle plume activity characterized by OIB-type mantle in Bachu area, Tarim. On the contrary, they are similar to the dykes in Yinaoxia and Podong area. Trace elements and Sr-Nd isotopes indicate that the parental magma for these mafic dykes was derived from the partial melting of the depleted asthenosphere mantle that metasomatized previously by the fluids from ancient subducting slab. This suggests that the dolerite dyke swarms in Beishan area were formed after the closure of the Paleo-Asian Ocean. It was a product of basaltic magmatism formed by the partial melting of the decompressed asthenosphere caused by the lithospheric delamination in a post-subduction extensional setting.

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 authors.

AUTHOR CONTRIBUTIONS
GX: field geological investigation, ideas, formal analysis, and writing -original draft. JD: field geological investigation, validation, and writing -review and editing. WG: data curation and visualization. RW: project administration and validation. ZS: field geological investigation and visualization. BM: data collection and software. JS: experimental test assistance and software. All authors contributed to the article and approved the submitted version.