- 1School of Resources and Environment Engineering, Anhui University, Hefei, China
- 2Exploration Research Institute, Anhui Provincial Bureau of Coal Geology, Hefei, China
- 3Faculty of Science, University of Brunei Darussalam, Bandar Seri Begawan, Brunei
- 4Third Exploration Team of Anhui Provincial Bureau of Coal Geology, Suzhou, China
In the eastern North China Craton (NCC), Mesozoic tectonics was dominated by the Paleo-Pacific subduction rollback and the Tanlu crustal-scale fault movement. The regional transtension had generated extensive adakitic magmatism, some Cu-Au ore-forming but others not. To establish the geodynamic setting and any metallogenic link for the adakites from the southeastern NCC margin, we analyzed the ore-barren adakitic rocks from underground mines in the Huaibei-Linhuan coalfield (where surface igneous outcrops are scarce), and compared their ages and geochemistry with other mineralized and ore-barren adakites across Eastern China. Zircon U-Pb dating reveals two magmatic episodes in the Huaibei-Linhuan coalfield: 1) early-Early Cretaceous (ca. 130–129 Ma) (quartz-)diorite and granodiorite, and 2) late-Early Cretaceous (ca. 115.8 and 105.8 Ma) microgabbro and dolerite. Whole-rock geochemistry indicates that the (quartz-)diorite and granodiorite are high-Mg adakitic, featured by low K2O/Na2O (avg. 0.33), high Sr/La (avg. 44.3), and lack of correlation between SiO2 (fractionation index) and Sr/Y (avg. 56.55) and MREE/HREE (avg. 1.09), resembling typical adakites derived from oceanic-slab partial melting. Geochronological correlation with the regional tectonic events suggests that the slab-melting may have been caused by the Paleo-Pacific subduction rollback. Further extension and crustal thinning in the late-Early Cretaceous along the southern Tanlu fault may have formed the gabbro-dolerite in the coalfield. Geochemical comparison suggests that parental magma of the Huaibei-Linhuan adakites may have had similar water content [similar zircon 10,000*(Eu/Eu*)/Y and Eu/Eu* ratios] to typical porphyry Cu-Au ore-forming magmas, yet the former may have been considerably more reduced (lower zircon Ce/Nd and whole-rock V/Sc ratios). We considered that the assimilation of Carboniferous-Permian coal seams in the area may have further lowered the magma fO2 and thus its potential to form Cu-Au mineralization.
Introduction
The Mesozoic eastern North China craton (NCC) is tectonically complex due to the co-influence of a number of tectonic events, including the North China-South China collision, the Paleo-Pacific subduction and rollback, the gradual NCC decratonization, together with the ∼500-km sinistral movement of the regional Tanlu fault zone and the accompanied pull-apart basin formation (e.g., Li et al., 2012, 2017; Hong et al., 2020; Zhang et al., 2020; Sun et al., 2021; Ye et al., 2021). The complex tectonic interactions resulted in the formation of important Cu-Au mineralization across the region, notably the world-class Jiaodong-Liaodong orogenic gold province (e.g., Goldfarb and Santosh, 2014; Zhang Z. et al., 2019; Liu et al., 2019, 2020; Li J. et al., 2020; Deng et al., 2020), and the Cu-Au-Fe polymetallic Middle-Lower Yangtze River Metallogenic Belt (MLYRB) (e.g., Pirajno and Zhou, 2015; Zhou et al., 2015; Zhang Y. et al., 2019; Lü et al., 2021). Many of these deposits are closely related in space-time to adakitic magmatism, yet many adakitic units in the region are also ore-barren (e.g., Liu S.-A. et al., 2010; Li C. et al., 2020; Jiang et al., 2020). Genetic link between adakitic magmatism and Cu-Au mineralization is rather well-established (e.g., Reich et al., 2003; Richards and Kerrich, 2007; Chiaradia et al., 2012), yet understanding why some adakites are ore-barren would also benefit both ore deposit modelling and regional Cu-Au exploration (Table 1).
TABLE 1. Summary of ages, intrusions and lithology of major barren and fertile adakite units in Eastern China.
In the southeastern NCC margin (Anhui province) where the North China Alluvial Plain is currently situated, intrusive rocks are rarely exposed and thus less studied. In this study, we collected samples from various underground coal mines across the Huaibei-Linhuan coalfield (Figure 1). Through zircon U-Pb dating and whole-rock geochemical analyses, and by comparing geochemically with the coeval Cu-Au mineralized/barren adakites across Eastern China (esp. in the MLYRB and along/near the southern Tanlu fault), we discuss why some adakites are ore-forming and some are not.
FIGURE 1. Tectonic maps of China (A) and North China (B), and (C) geologic map of the Huaibei-Linhuan coalfield, showing the distribution of basement rocks, sampling sites and major structures (modified after Zhao et al., 2016; Zhu et al., 2005, 2009; An et al., 2018).
Geological Background
Archean to Mesozoic geological record is well-preserved in the NCC, except locally in its eastern part due to later decratonization. The craton was largely formed by extensive Archean tonalite-trondhjemite-granodiorite (TTG) emplacement, and remained relatively stable throughout the Paleozoic-Early Triassic. Afterward, the convergence between the Yangtze and North China plates (collided in the Early-Middle Jurassic) had caused widespread Mesozoic magmatism, forming the Cretaceous intermediate-mafic and bimodal (sub-)volcanic rocks in the southern NCC margin (Li et al., 2017). Some intrusive rocks were suggested to have formed by mantle underplating and had a crustal-mantle mixed source. Geochemistry of these rocks implies that the lithospheric mantle beneath the NCC is highly heterogenous (Yang et al., 2008; Zhao et al., 2009; Xia et al., 2013).
Although igneous rock exposure is generally poor in the southeastern NCC margin, many intrusions are encountered in the underground coal mines there, notably those in the Huaibei-Linhuan coalfield. The Huaibei coalfield is located close to the Xuzhou–Suzhou fault system and between the EW-trending Fengpei and Bangpu uplift zones (Figure 1B). Faults of various sizes and types were documented in the coalfield, but consist mainly of high-angle normal faults and rare reverse faults. The regional Huangji and Subei faults were interpreted to have formed in the same stage as the regional Carboniferous-Permian folding (which formed the Tongting syncline and Nanping anticline) (Figure 1C) (An et al., 2018). The Linhuan coalfield is bounded by the N-S-trending Taihe-Wuhe and Subei faults, and by the E-W-trending Nanping and Fengwo faults, giving a rectangular shape to the coalfield. The Taihe-Wuhe fault is a steeply-dipping normal fault, which connects to the Tanlu fault zone in the east. Secondary faults and interlayer structures are well-developed in the coalfield: the former is mainly NE-trending and less-commonly NW-trending, whilst the latter is developed on both hanging-wall and footwall of the coal seams, which locally damaged the coal seams. Major fault activities likely commenced in the Jurassic (early-stage Yanshanian orogeny), during which the regional stress transition (to NW-SE-trending) may have formed the Xuzhou-Suzhou imbricated structure. Tectonic regime in Eastern China changed again to active extension during the middle- and late-stage Yanshanian orogeny, forming many normal faults and fault-bounded basins. Normal faulting likely reached its peak during the Cenozoic Himalayan orogeny, during which many ancient structures were also reactivated/modified (Yao and Liu, 2012). Intermediate volcanism and plutonism were active during the Yanshanian movement, and were mainly controlled by the NNE-trending and EW-trending faults (Xu et al., 2004; Yang et al., 2010).
Previous geological survey documented that stratigraphy in the Linhuan coalfield comprises Ordovician and Carboniferous-Permian sequences (#3 Exploration Team of Anhui Bureau of Coal Geology, 2007). The coal-bearing sequences include the Upper Carboniferous Benxi and Taiyuan formations, Lower Permian Shanxi Formation, Middle Permian Shihezi Formation, and the Middle-Upper Permian Shiqianfeng Formation. The Carboniferous sequences comprise mainly greyish mudstone in the lower part, and crystalline limestone and sandy shale-mudstone in the upper part. The Permian sequences comprise mainly sandstone and mudstone (An et al., 2018, and ref. therein). Plutonic emplacement is locally controlled by ancient faults and interlayer fractures, which likely served as magma conduits. Many Jurassic-Cretaceous intermediate-mafic dykes were identified in the Huaibei coalfield, and are controlled mainly by NNE-trending faults and minor by EW-trending and NW-trending ones. Middle Jurassic-Early Cretaceous intermediate-felsic plutons are widespread in the southern NCC margin, and eclogite xenoliths and adakitic rocks were reported in the Xuhuai area as a rectangle tectonic region in the Southern Tan-Lu Fault (STLF) showing in Figure 1B (Xu et al., 2006; Xu et al., 2009). Plutonic rocks intruding the coal-bearing sequence in the Huaibei coalfield include mainly diorite porphyry stock, dykes, and sills. All coal mines in the Linhuan coalfield were affected (metamorphosed) by various degrees of plutonism, which intensified toward the Subei and Daliujia faults, e.g., at the Yangliu and Haizi mines (#3 Exploration Team of Anhui Bureau of Coal Geology, 2007; Yangliu Coal Industry Co. Ltd, 2017; Zhao et al., 2019).
Intrusive rock samples in this study were collected from coal mines in the Huaibei-Linhuan coalfield. Plutons are mainly dioritic in the coal mines along the NE-trending Qiluzhou fault (Taiping, Haizi, Yangliu), whilst those in the coal mines along the NS-trending Daliujia fault (Yuandian-1, Yuandian-2, Qingdong) are mainly gabbroic and minor peridotitic.
Sampling and Analysis Methods
In this study, we analyzed igneous rock samples from six mines in the Huaibei-Linhuan coalfield, with the sampling locations shown in Figure 1C. Sample lithologies include (quartz-)diorite-granodiorite (n = 6) from Yangliu (YL-03, YL-07), Haizi (HZ-08, HZ-13), and Taiping (TP-01, TP-04), and microgabbro-dolerite (n = 6) from Qingdong (TT-07, TT-12), Yuandian-1 (Y1-29, Y1-32), Yuandian-2 (Y2-08, Y2-23). The rock samples were prepared into petrographic thin sections and rock powder (200 mesh) for microscopic (including SEM) observations and whole-rock geochemical analyses, respectively. Sample preparation and petrographic observation were conducted at the Geological Laboratory of the School of Resources and Environment Engineering, Anhui University. Optical microscopic analysis was performed with an Olympus BX53PLM to identify textural features and rock-forming minerals. Whole-rock major and trace elemental compositions were measured on all samples by X-ray fluorescence (XRF) and inductively coupled plasma mass spectrometry (ICP-MS), respectively.
The XRF analysis follows the Chinese national standard procedure GB/T 14506.28-2010. The powdered samples were dried at 105°C for 2–4 h and then cooled down in a desiccator. For each sample, 0.700 g was weighted and mixed with 5.200 g Li2B4O7, 0.400 g LiF, and 0.300 g NH4NO3. The mixture was placed in platinum crucibles (95% Pt + 5% Au). Afterward, 1 ml LiBr (1.5%) was added to the crucibles, dried on a hot-plate, and then heat (with lid closed) at 1,200°C for 10 min to form a glass disc.
Solution ICP-MS sample preparation follows the Chinese National Standard Procedure GB/T 30714-2014. Weighted dry sample powder (0.1000 g) was placed in polytetrafluoroethylene (PTFE) bottles, mixed successively with 10 ml HNO3, 10 ml HF, and 2 ml HCl, and then heated on hot-plate for 2–3 h to completely dissolve the samples. After cooling, the residue was washed with ultra-pure water, dried again, and then redissolved in 2.0 ml aqua reqia and 1.0 ml LH3BO3 on hot-plate for around 16 min. The solution was then transferred to 50 ml bottles, diluted with water and then measured (with two blanks) on an Agilent 7700e ICP-MS. The analysis follows the Chinese national standard procedure DZ/T0223-2001.
LA-ICP-MS zircon U-Pb dating was performed at the Wuhan Sample Solution Analytical Technology Co. Ltd. (China) with an Agilent 7900 ICPMS, coupled to a GeoLas HD 193 nm excimer laser ablation system. Zircon 91500 and glass NIST610 were used as external standards for U-Pb dating and trace element calibration, respectively. Each analysis incorporated a background acquisition of approximately 20–30 s followed by 50 s of data acquisition from the sample. An Excel-based software ICPMSDataCal was used to perform off-line selection and integration of background and analyzed signals, time-drift correction and quantitative calibration for trace element analysis and U-Pb dating (Liu et al., 2008; Liu et al., 2010 Y.). Concordia diagrams and weighted mean calculations were made using Isoplot/Ex_ver3 (Ludwig, 2003). Analytical conditions include: 80 mJ laser energy, 5 Hz frequency, and 32 µm laser beam size. The standards were analyzed twice for every six samples. The data processing and concordia plot production were completed with ICPMSDataCal 10.9 and Isoplot 3.0, respectively.
Results
Petrographic Features
Taiping mine: Diorite samples were collected from a dyke (4–6 m thick) at 300 m depth. The rocks are dark-grey, medium-fine-grained, and massive structure. The rocks contain mainly plagioclase (75%) and hornblende (20%) and rare pyroxene (<5%). Some samples are porphyritic and have zoned plagioclase phenocrysts in a dark grey to greyish-green microcrystalline groundmass. The plagioclase grains are largely euhedral-subhedral, whilst hornblende grains are subhedral (Figures 2A,D).
FIGURE 2. Hand-specimen and thin-section microscopic photos of subvolcanic rocks from the Huaibei-Linhuan coalfield: (A) Taiping (TP-01) diorite; (B) Haizi (HZ-08) diorite; (C) Yangliu (YL-03) diorite; (D) Taiping granodiorite with feldspar phenocrysts; (E) Haizi porphyritic diorite with altered euhedral feldspar phenocrysts; (F) Yangliu granodiorite porphyry with quartz and plagioclase phenocrysts; (G) Qingdong (TT-07) microgabbro; (H) Yuandian-1 (Y1-29) dolerite; (I) Yuandian-2 (Y2-23) dolerite; (J) Pyroxene and plagioclase in the Qingdong microgabbro (XPL); (K) Plagioclase and pyroxene in the ophitic-textured Yuandian-1 dolerite (XPL); (L) Plagioclase and fine-grained Fe-Ti oxides in the Yuandian-2 dolerite (XPL).
Yangliu mine: The diorite-granodiorite samples were collected from a sill at 463–504 m depth (drill-hole 2017-02). The rocks are weakly porphyritic, with subrounded quartz and alkali feldspar phenocrysts and (minor) hornblende phenocrysts set in a dark-grey hornblende groundmass. Feldspar phenocrysts (2 mm) are generally larger than subhedral hornblende phenocrysts (0.5 mm) (Figures 2C,F).
Haizi mine: The diorite porphyry samples were collected from a dyke (8–10 m thick) at 800 m depth. The rocks are dark grey and medium-fine-grained massive, and contain mainly phenocrysts of plagioclase (60%), hornblende (20%) and pyroxene (<5%), which set in a dark-grey hornblende and (minor) pyrite groundmass (Figures 2B,E). Partially assimilated porous and rounded-/irregular-shaped natural coke fragments are included in the diorite (Figure 3A).
FIGURE 3. SEM image of plutonic rock samples, showing partially-assimilated porous natural coke inclusions with rounded corroded margin: (A) HZ-13; (B) TT-07; (C) Y1-32.
Qingdong mine: The microgabbro samples were collected from a sill at 478–483 m depth (drill-hole 2016-03). The rocks are greyish-green, fine-grained massive, and contain mainly plagioclase feldspars (55%), pyroxene (35%) and olivine (10%) (Figures 2G,J). Partially assimilated natural coke fragments, featured by being porous and rounded/irregular, are observed under SEM in the microgabbro (Figure 3B).
Yuandian-1 mine: The dolerite samples were collected from a sill at 295–320 m depth (drill-hole 2016-06). The rocks are greyish-green, fine-medium-grained massive (weakly porphyritic), and contain euhedral-subhedral plagioclase (65%) and pyroxene (25%) phenocrysts set in a groundmass (10%) of similar mineral content (Figures 2H,K). Some dolerite samples contain partially assimilated natural coke fragments (Figure 3C).
Yuandian-2 mine: The dolerite samples were collected from a sill at 278–305 m depth (drill-hole 2016-20). The rocks are greyish-green, fine-medium-grained massive, and weakly porphyritic. The rocks have mainly plagioclase (60%), pyroxene (20%), and minor (10% each) hornblende and Fe-Ti oxides. The pyroxenes are commonly altered to fine-grained amphibole and chlorite (Figures 2I,L).
Zircon U-Pb Geochronology and Trace Element Geochemistry
In this study, one Taiping diorite (TP-01) and Haizi diorite porphyry (HZ-08), and one Qingdong microgabbro (TT-07) and Yuandian-1 dolerite (Y1-32) were zircon U-Pb dated. Zircon grains from all rock samples are largely similar in size (50–150 μm, mostly ∼100 μm). Zircons from the diorite samples (TP-01 and HZ-08) are featured by rounded, partially resorbed inherited cores (20–120 μm, mostly 25–90 μm) surrounded by oscillatory zoned rims, giving an overall short prismatic or subrounded shape. Most zircon cores have markedly different CL reflectance (darker or brighter) from their growth rims. The cores are mostly unzoned (homogenous) but sector zoning is also found in certain grains. For the microgabbro-dolerite samples (TT-07 and Y1-32), most zircon grains do not have inherited core, and most of them are elongated prismatic (some occur as broken fragments), and are oscillatory- and/or sector-zoned (Figure 4).
FIGURE 4. Representative zircon CL images for the Huaibei-Linhuan intrusive rock samples: (A) Taiping diorite (TP-01); (B) Haizi diorite porphyry (HZ-08); (C) Qingdong microgabbro (TT-07); (D) Yuandian-1 dolerite (Y1-32).
Zircon samples from the Taiping diorite (TP-01) have high Th/U ratios for both the inherited cores (avg. 0.62) and rims (avg. 0.39), resembling typical igneous zircons (Hoskin and Black, 2000). Most age data (n = 19) from the cores cluster closely around 2,500 Ma, and yielded an upper intercept 206Pb/238U age of 2,463 ± 23 Ma (2σ, MSWD = 1.2). Meanwhile, age data (n = 21) from the rims yielded a weighted average 206Pb/238U age of 129.2 ± 1.4 Ma (2σ, MSWD = 2.0), which we interpreted to be the diorite emplacement age (Figures 4A, 5B; Table 2).
FIGURE 5. Zircon LA-MC-ICP-MS U-Pb Concordia diagrams of the Huaibei-Linhuan intrusive rock samples: (A) inherited zircon cores from Haizi diorite porphyry (HZ-08); (B) zircon growth rims from Haizi diorite porphyry (HZ-08); (C) inherited zircon cores from Taiping diorite (TP-01); (D) zircon growth rims from Taiping diorite (TP-01); (E) Qingdong microgabbro (TT-07); (F) Yuandian-1 dolerite (Y1-32).
TABLE 2. Zircon U-Pb dating results of the plutonic rock samples from the Huaibei-Linhuan coalfield.
Similarly, zircons from the Haizi diorite porphyry (HZ-08) have high Th/U ratios for both the inherited cores (avg. 0.47) and rims (avg. 0.84), resembling typical igneous zircons (Hoskin and Black, 2000). Age data (n = 19) from the cores also cluster closely around 2,500 Ma, and yielded an upper intercept 206Pb/238U age of 2,475 ± 33 Ma (2σ, MSWD = 2.3). Analyses of the zircon rims (n = 15) yielded a weighted average 206Pb/238U age of 130.0 ± 2.4 Ma (2σ, MSWD = 2.4), which we interpreted to be the diorite crystallization age (Figures 4B, 5A,C; Table 2).
For the Qingdong microgabbro (TT-07), the zircon ages range narrowly from 110.48 to 120.38 Ma and there are no old inherited zircons. The zircons are likely igneous due to their high Th/U ratios (avg. 0.50) and the ubiquitous igneous oscillatory zoning. Twenty-five analyses yielded a weighted average 206Pb/238U age of 115.8 ± 1.1 Ma (2σ, MSWD = 2.0), which likely reflects the timing of the microgabbro formation (Figures 4C, 5B,E; Table 2).
For the Yuandian-1 dolerite (Y1-32), most zircon ages range narrowly from 104.38 to 124.48 Ma. Again, the zircons are likely igneous due to their high Th/U ratios (avg. 0.93) and their igneous oscillatory zoning. Fifteen analyses yielded a weighted average 206Pb/238U age of 105.8 ± 1.8 Ma (2σ, MSWD = 2.8), which is largely the same (within error) as the upper intercept 206Pb/238U age (106.6 ± 2.4 Ma) and likely represents the dolerite formation age (Figures 4D, 5F; Table 2).
Most zircons, regardless of their host-rock types and ages, have similar left-inclining REE patterns and positive Ce anomalies (Ce4+/Ce3+ = mostly 0.21–90.78), resembling typical igneous zircons (Belousova et al., 2002). LREE contents are mostly low (often below the detection limit), and the few zircons with abnormally high La content are likely caused by apatite inclusions (Lu et al., 2016). Different rock samples, however, have different zircon Eu anomalies: Zircon rims of samples HZ-08 and TP-01 have no discernible anomalies [δEu: (Eu*/Eu) = 0.42–0.86], whereas their Neoarchean-Paleoproterozoic cores have distinct negative anomalies (δEu = 0.03–0.48). Varying negative Eu anomalies are present for both samples TT-07 (δEu = 0.03–0.84) and Y1-32 (δEu = 0.12–0.55) (Figure 6; Table 3).
FIGURE 6. Chondrite-normalized zircon REE patterns for the Huaibei-Linhuan intrusive rock samples: Chondrite normalizing values are from Sun and McDonough (1989).
Whole-Rock Elemental Compositions
For the 12 samples analyzed, their whole-rock major element contents were back-calculated to 100 wt% to remove the volatile influence. A few samples yielded very high LOI (e.g., TT07 = 12.36 wt%, Y2-23 = 10.19 wt%), which we considered to be caused by the natural coke inclusions (instead of alteration/weathering), as demonstrated by the fresh igneous texture and partially-assimilated natural coke fragments observed under the SEM (Figure 3). The microgabbro-dolerite samples have SiO2 = 42.58–50.88 wt%, Al2O3 = 12.22–14.55 wt%, MgO = 5.44–8.63 wt%, FeOT = 6.99–14.68 wt%, and (Na2O + K2O) = 2.62–4.08 wt%, whilst the diorite-granodiorite (Haizi, Taiping, and Yangliu) samples have SiO2 = 57.45–61.27 wt%, Al2O3 = 14.53–16.12 wt%, MgO = 2.87–4.63 wt%, FeOT = 4.0–6.71 wt%, and (Na2O + K2O) = 6.45–9.15 wt% (Table 4).
TABLE 4. | Major and trace element compositions of plutonic rock samples from the Huaibei-Linhuan coalfield.
In the Total Alkali-Silica (TAS) diagram, the gabbro-dolerite samples fall inside the (olivine-)gabbro fields (some marginally alkali), whilst the diorite-granodiorite samples fall inside the (quartz-)monzonite and syenite fields (Figure 7A). Similar results are shown in the immobile Zr/Ti vs. Nb/Y diagram, suggesting that alteration/weathering has no major influence on the major element geochemistry (Figure 7B). Both the gabbro-dolerite (ASI = 0.94–1.38) and diorite-granodiorite (ASI = 1.09–1.42) are peraluminous (Figure 7C), with the former being calc-alkaline and latter high-K calc-alkaline (Figure 7D). The diorite-granodiorite samples have high Sr (>677.7 ppm) and low Y (<15.48 ppm) (Sr/Y = 49.97–106.95), showing adakitic affinity. In contrast, the gabbro-dolerite samples have lower Sr (<384.2 ppm) and higher Y (>24.13 ppm)(Sr/Y = 15.26–83.19, avg. 47), resembling normal arc magmatic rocks (Figure 8A). The gabbro-dolerite and diorite-granodiorite samples have similar total (Σ) REE contents [73.14–193.10 ppm (avg. 112.88) and 50.51–113.45 ppm (avg. 84.99), respectively] and chondrite-normalized REE patterns, which are featured by enrichments in LREE/MREE [(La/Yb)N = 4.27–17.85 and 8.63–10.62, respectively] and MREE/HREE [(Dy/Yb)N = 1.76–2.16 and 1.81–1.98, respectively] (Tables 4; Figures 7–8). In primitive mantle-normalized diagrams, both the gabbro-dolerite and diorite-granodiorite show large ion lithophile element (LILE) over high field strength element (HFSE) enrichments, as well as negative Ta-Nb and Ti anomalies and distinct positive Pb and Sr anomalies.
FIGURE 7. Petrological classification plots for the igneous rock samples from the Huaibei-Linhuan coalfield: (A) Total Alkali Silica (TAS) (base map after Middlemost, 1994). The curve denotes the Ir (Irvine) boundary, separating the alkaline (above) and sub-alkaline (below). Number: 1) olivine gabbro; 2a) alkaline gabbro; 2b) sub-alkaline gabbro; 3) gabbroic diorite; 4) diorite; 5) granodiorite; 6) granite; 7) quartzolite; 8) monzogabbro; 9) monzodiorite; 10) monzonite; 11) quartz monzonite; 12) syenite; 13) foid gabbro; 14) foid monzogabbro; 15) foid monzosyenite; 16) foid syenite; 17) foidolite; 18) tawite/urtite/italite; (B) Zr/Ti vs. Nb/Y (after Pearce, 1996); (C) A/NK vs. A/CNK (base map after Maniar and Piccoli, 1989); (D) K2O vs. SiO2 (after Peccerillo and Taylor, 1976). For Figures 5–6, 9–10, data from the adakitic rocks in the Bengbu-Wuhe gold belt (Li C. et al., 2020), Banjing pluton (Zhou et al., 2019), Eastern Dabie Orogen (Zhang et al., 2013), Lower Yangtze River Belt (Liu S.-A. et al., 2010) and southern Tanlu fault (Liu S.-A. et al., 2010) are also shown for comparison.
FIGURE 8. Adakite discrimination plots for the igneous rock samples from the Huaibei-Linhuan coalfield: (A) Sr/Y vs. Y (after Defant and Drummond, 1993), andesites, dacites and rhyolites (ADR); (B) K2O/Na2O vs. Al2O3 (after Liu et al., 2010a, and ref. therein); (C) MgO vs. SiO2 (base map after Liu et al., 2010a, and ref. therein); (D) Sr/Y vs. SiO2 (after Richards et al., 2012).
Discussion
Multiphase Magmatism in the Huaibei-Linhuan Coalfield
Intrusive rocks dated in this study yielded at least two magmatic phases: The Early Cretaceous diorite-granodiorite ages (Haizi: 130.0 ± 2.0 Ma; Taiping: 129.2 ± 1.4 Ma), whilst the microgabbro-dolerite yielded Late Cretaceous ages (Qingdong: 115.8 ± 1.1 Ma; Yuandian-1: 105.8 ± 1.8 Ma). The Early Cretaceous Huaibei-Linhuan diorite-granodiorite are coeval with the high-Mg adakite-like diorite-granodiorite in the Liujiawa pluton (∼128 Ma), which is located further southwest in the eastern Dabie Orogen close to the Tanlu fault zone (Figure 1) (Jiang et al., 2020). Also along the Tanlu fault, the ore-related Nvshan syenogranite and the Huaiguang granodiorite in the Bengbu-Wuhe goldfield (Anhui Province) were both dated to be ∼130 Ma (Li C. et al., 2020). The Huaibei-Linhuan (quartz-)diorite-granodiorite are also of similar age (slightly younger) than the Tongguanshan Cu-Au ore-forming adakites (∼136.7 Ma; Tongling orefield), and the ore-barren high-Mg adakitic plutons at Fangjiangzhuang (∼129.1 Ma) and Qiaotouji (∼131.7 Ma) in the MLYRB (Liu S.-A. et al., 2010; Jiang et al., 2020). This suggests that the emplacement of the Huaibei-Linhuan (quartz-)diorite and granodiorite was part of the extensive Early Cretaceous (high-Mg-)adakitic magmatism, which extended from the along the eastern MLYRB northward along the Tanlu fault. From Eastern China, the Tanlu fault zone likely continues NE-ward across the Bohai sea basin into NE China, although its northeastern-most end is yet to be well identified (Ye et al., 2021). Early Cretaceous adakitic granitoids were documented at Xiuyan (ca. 129–126 Ma) in the Liaodong Peninsula (NE China) (Dong et al., 2020), and Krinichnoye (131.4 Ma) in the southern Sikhote-Alin (SE Russia Far East). As to be discussed in the later section, we attributed this major regional adakitic magmatism as a product of Paleo-Pacific subduction and rollback and/or tearing, facilitated by the reactivation of the Tanlu fault zone movement.
The younger late-Early Cretaceous magmatic phase in Huaibei-Linhuan (Qingdong: 115.8 Ma; Yuandian-1: 105.8 Ma) was also reported elsewhere along the Tanlu fault. For instance, normal-arc-type granite- and diorite porphyries and monzogranite from the Bengbu-Wuhe goldfield were dated at ca. 112–114 Ma, and the associated Bengbu mafic dykes yielded similar Ar-Ar age (plateau: 111.8 ± 0.6 Ma; isochron: 109.6 ± 1.5 Ma) (Liu et al., 2012; Li C. et al., 2020).
Petrogenesis of the Huaibei-Linhuan Magmatic Rocks
In this study, we have identified an earlier adakitic (quartz-)diorite-granodiorite phase and a younger gabbro-dolerite phase for the Early Cretaceous magmatism in Huaibei-Linhuan. For the adakitic (quartz-)diorite-granodiorite, they have relatively high Mg (MgO = 2.87–4.63; Mg# = 0.65–0.76) and low K2O/Na2O (0.35–0.54), resembling typical adakites derived from oceanic slab melting. Such features are distinct from the adakites derived from melting of the thickened lower crust, notably those from the Bengbu-Wuhe goldfield situated in the interior of the North China Craton (Figure 8) (Li C. et al., 2020). It is not to say that there is no crustal influence in the Huaibei-Linhuan adakitic rocks, as evidence from the presence of Precambrian inherited cores (2,317–2,632 Ma) in the Cretaceous zircons. The weighted average ages of the inherited cores for the (quartz-)diorite-granodiorite samples from Haizi (2,475 ± 33 Ma) and Taiping (2,463 ± 23 Ma) are similar to the relict magmatic zircons in (ca. 2,455 and 2,716 Ma) from the Cenozoic (∼1 Ma) Nushan basalts in the southeastern North China Craton (Ping et al., 2019), suggesting clear NCC crustal input for the Huaibei-Linhuan adakite formation. The lack of the clear correlation between the adakitic signature (Sr/Y) and fractional crystallization (SiO2) suggest that the adakites were not formed by high-pressure fractionation (garnet and/or amphibole) (Figure 7D). This conclusion is also supported by the lack of upward-concave MREE-HREE pattern in the chondrite-normalized REE plots (Figure 9), or the lack of correlation between MREE/HREE (Dy/Yb) and fractionation (SiO2), which would be expect for amphibole-garnet fractionation-derived adakites or normal calc-alkaline rocks (such as the Huaibei-Linhuan gabbro-dolerite) (Castillo et al., 1999; Kamvong et al., 2014) (Figure 10A). Previous studies suggested that if the adakites were formed by fractional crystallization, substantial amphibole (85%) and garnet (15%) fractionation are needed (Mori et al., 2007). Our (quartz-)diorite-granodiorite samples do not fall on these fractionation trends (Figures 10A–B), and no garnet (and very few amphibole) was found in the rocks, suggesting that the adakites were not formed from fractional crystallization. In the Rb/Y vs. Nb/Y and Ba/Th vs. Th discrimination plots, the Huaibei-Linhuan adakite samples fall between the slab-derived and melt-derived enrichment trends, suggesting that the magmas were sourced from both slab-derived fluids and melts (Figures 10C–D). The Sr/La ratio has been used to distinguish slab-derived adakites from lower-continental crustal (LCC)-derived ones, because altered oceanic crust (with MORB-type LREE depletion and Sr enrichment by seawater) has much higher Sr/La than the LCC (Liu S.-A. et al., 2010). Our adakitic (quartz-)diorite-granodiorite samples have elevated Sr/La ratios (Figure 11A), indicating a slab-derived origin. The slab-derived origin is also supported by the zircon trace element compositions. The zircons show in general increasing Yb/Gd with Hf content, which demonstrate a magma cooling and fractionation trend (Figure 12A) (Barth and Wooden, 2010). The Cretaceous zircons (or zircon rims) have elevated U/Yb ratio (a proxy of LILE-enrichment), suggesting a more enriched magma source resembling that of typical continental arc (Grimes et al., 2015). This is markedly different from the Precambrian cores, which have lower U/Yb ratio that mimics zircons from typical ocean islands (e.g., Iceland). The Cretaceous and Precambrian zircons have similar Hf and Nb/Yb, showing their host magmas have similar degree of fractionation and alkalinity, respectively (Figures 12B–C) (Grimes et al., 2015).
FIGURE 9. Chondrite-normalized whole-rock REE patterns for the Huaibei-Linhuan intrusive rock samples: (A) gabbro-dolerite; (B) (quartz-)diorite and granodiorite. Chondrite normalizing values are from Sun and McDonough (1989).
FIGURE 10. (A–B) Trace element plots for identifying amphibole–garnet fractionation: (A) Dy/Yb vs. SiO2. Model fractionation paths are from Davidson et al. (2007), Macpherson (2008); (B) Dy vs. Nb; (C) Nb/Ta vs. Nb. Fractional crystallization (FC) and assimilation-fractional crystallization (AFC) trends are from Gómez-Tuena et al. (2007); (C–D) Trace element plots for differentiating fluid-related vs. melt-related enrichments (after Zhao et al., 2020): (C) Rb/Y vs. Nb/Y; and (D) Ba/Th vs. Th. Abbreviations: amp = amphibole; cpx = clinopyroxene; grt = garnet, pl = plagioclase.
FIGURE 11. Petrogenetic discrimination plots for the Huaibei-Linhuan intrusive rock samples: (A) Ce/Pb vs. Sr/La (after Liu S.-A. et al., 2010); (B) V/Sc vs. MgO (after Lee et al., 2005; Kamvong et al., 2014). V/Sc = proxy for fO2 of parental magma. Data from the ore-related porphyries from Quimsacocha (Ecuador) (Beate et al., 2001) and Dexing (Eastern China) (Wang et al., 2006) and ore-barren calc-alkaline rocks from the Sunda arc (Reubi and Nicholls, 2004) and Aleutian arc (Jicha and Singer, 2006) are also shown for comparison.
FIGURE 12. Plots of (A) Yb/Gd vs. Hf (base map from Barth and Wooden, 2010), (B) U/Yb vs. Hf (base map from Grimes et al., 2015), and (C) U/Yb vs. Nb/Yb (base map from Grimes et al., 2015) for zircons from the Huaibei-Linhuan intrusive rock samples.
After the magma formation, the lack of distinct Ce or Eu anomalies (some even with weakly positive Eu anomaly), indicates that plagioclase fractionation was insignificant. In the chondrite-normalized REE and primitive mantle-normalized multi-element plots (Sun and McDonough, 1989) (Figures 9, 13), the adakitic (quartz-)diorite-granodiorite and the non-adakitic gabbro-dolerite have generally similar arc-type patterns, e.g., LREE enrichments (between EMORB and OIB), HREE depletions (relative to MORB), negative Nb-Ta-Ti anomalies, and positive Pb-Sr anomalies, suggesting that both rock suites were formed in a subduction-related setting. The only difference between them (apart from the adakitic characters such as the Sr-Y contents), i.e., higher LILE and LREE contents for the (quartz-)diorite-granodiorite, is likely attributed to its more fractionated nature compared with the gabbro-dolerite.
FIGURE 13. Primitive mantle-normalized spider diagrams for the Huaibei-Linhuan intrusive rock samples: (A) gabbro-dolerite; (B) (quartz-)diorite and granodiorite. Primitive mantle normalizing values are from Sun and McDonough (1989).
Tectonic Evolution of Cretaceous Eastern China
Cretaceous adakitic magmatism in Eastern China has been variably attributed to the progressive decratonization of North China, and/or the Paleo-Pacific subduction and rollback. As discussed in the previous section, our Huaibei-Linhuan (quartz-)diorite-granodiorite samples show geochemical features typical of oceanic slab melting-derived adakites, but also contain Precambrian inherited zircon core (Haizi: 2,475 ± 33 Ma; Taiping: 2,463 ± 23 Ma) that reflect NCC crustal inheritance (Cenozoic Nushan basalt: ca. 2,455 and 2,716 Ma) (Ping et al., 2019). An oceanic slab-melting model was also proposed for the adakites (some high-Mg) in the nearby Banjing pluton (Xuhuai area), which is slightly younger (zircon U-Pb: 126.4 ± 2.1 Ma) and has similar structural setting (proximal to the Tanlu fault) to the Huaibei-Linhuan adakites (Zhou et al., 2019). The Huaibei-Linhuan and Banjing are geochemically similar in terms of their low Th content and low Rb/Y and Nb/Y ratios (suggestive of minor crustal melt-derived enrichment), which are in marked contrast with the thickened lower crust-derived Bengbu-Wuhe gold ore-forming granitoids (Liu S.-A. et al., 2010) (Figures 10C–D).
Although not without controversies, many studies suggest that the Tanlu fault zone was first formed as a Triassic-Jurassic syn-collisional sinistral fault between the North China and South China (phengite Ar-Ar: ca. 221–210 and 198–181 Ma) (Zhu et al., 2009), and had displaced the Dabie–Sulu UHP metamorphic belt across a ∼500 km distance (Zhu et al., 2005). The fault (esp. the southern part where the Huaibei-Linhuan coalfield is located) has likely reactivated several times during the Early Cretaceous (e.g., hornblende Ar-Ar = 143.3 Ma; phengite Ar-Ar = 138.8 Ma) (Zhu et al., 2005). The sinistral movement in the southern Tanlu fault (138.8 Ma; Zhu et al., 2005) was accompanied (or soon followed) by emplacement of granodiorite (∼131 Ma), diorite (∼130 Ma), and lamprophyre dykes (∼128 Ma) of the Guandian intrusive complex (Zhao et al., 2020), whose ages are largely coeval (within error) with the ca. 130–129 Ma ages obtained for the Huaibei-Linhuan (quartz-)diorite-granodiorite. Based on whole-rock geochemical and Sr-Nd isotope geochemical data, Zhao et al. (2020) suggested that these Early Cretaceous intrusive rocks at Guandian were formed in an extensional setting, possibly related to the Paleo-Pacific subduction rollback. Such tectonic setting would fit with our suggestion that the Huaibei-Linhuan high-Mg adakitic (quartz-)diorite-granodiorite was formed by partial melting of the subducted slab. A slab-tearing may also be possible (Zhao et al., 2016) to generate the adakites, although more work will be needed to constrain when and how it actually took place.
A detailed review of the Early Cretaceous tectonics across the North Pacific by Liu et al. (2021) has revealed that the Eastern Eurasia, as facilitated by the Paleo-Pacific subduction rollback (via eastward migration of the shallow mantle convection system), has undergone three major extensional events at around 160–145 Ma, 135–120 Ma (peak), and post-120 Ma. These extensional events were accompanied by synextensional magmatism at around 160 Ma, 130–120 Ma (magmatic flare-up), and 100–80 Ma. It is thus clear that the Huaibei-Linhuan high-Mg adakitic magmatism, like many Cu-Au ore-fertile or ore-barren adakitic rocks in the southern Tanlu fault and the LMYRB, were formed during this 130–120 Ma magmatic flare-up during the peak extension episode. Meanwhile, the late-Early Cretaceous gabbro-dolerite magmatism in Huaibei-Linhuan (Qingdong: 115.8 Ma; Yuandian-1: 105.8 Ma) was formed in the post-peak extension. The highly attenuated North China continental crust around the southern Tanlu fault may have enabled higher degree of partial melting with lower degree of crustal assimilation and fractionation, thereby forming more mafic rocks in Huaibei-Linhuan. The presence of highly-attenuated crust is also supported by the complete absence of inherited zircon core or xenocrystic zircon in the Huaibei-Linhuan gabbro-dolerite.
Metallogenic Implications in Fertile Versus Barren Adakites
From a metallogenic perspective, it is still unclear on why some Early Cretaceous adakites in Eastern China are Cu-Au ore-fertile (e.g., those in the MLYRB), whereas some others are ore-barren (e.g., those along the southern Tanlu fault). Based on whole-rock elemental and Sr-Nd-Pb isotope data, Liu S.-A. et al. (2010) argued that the ore-fertile MLYRB adakites were formed from partial melting of the seawater-metasomatized Paleo-Pacific slab, whereas the ore-barren adakites along the southern Tanlu fault were derived from the delaminated North China eclogitic lower continental crust. This explanation, however, is unlikely applicable to the ore-barren Huaibei-Linhuan adakitic (quartz-)diorite-granodiorite, as their formation was likely slab-melting-related similar to the MLYRB adakites. This suggestion is supported by the elevated Sr/La values, which are comparable to the ore-fertile MLYRB adakites but considerably higher than the ore-barren adakites along the southern Tanlu fault (Figure 12A).
Porphyry-type Cu-Au fertility is controlled by a number of factors, including the contents of metals, water and sulfur of the magma, the magma oxygen fugacity (fO2), as well as thickness of the crust into which the magma was emplaced (e.g., Audétat et al., 2008; Richards et al., 2012; Park et al., 2019, 2021). Comparing with our Huaibei-Linhuan samples, both the mineralized adakites in the MLYRB and the barren adakites in the southern Tanlu fault are more fractionated (Liu S.-A. et al., 2010) (Figure 7D). The MLYRB adakites have distinctly higher Sr/Y (67.96–121.89) than their barren counterparts in Huaibei-Linhuan (15.26–106.95) and the southern Tanlu fault (42.49–81.19). The “more-adakitic” (higher Sr/Y) character in the MLYRB rocks may have caused by the higher magma water content (i.e., “wetter magma”), which can suppress plagioclase fractionation and promote early amphibole ± garnet fractionation (Richards, 2011, and ref. therein).
Previous studies suggested that elevated zircon Eu/Eu* (>0.3) can reflect elevated magma water content, as widely reported from ore-forming porphyries at El Salvador and Chuquicamata-El Abra (Chile), Yanacocha (Peru), Yerington (Nevada, US), Batu Hijau (Indonesia), Tampakan (Philippines), and Dexing, Qulong, and Jiama (China) (Ballard et al., 2002; Dilles et al., 2015; Lu et al., 2016). Lu et al. (2016) proposed that wet fertile magmas have commonly high 10,000*(Eu/Eu*)/Y and Eu/Eu* ratios, with both ratios being positively correlated with (Ce/Nd)/Y. Zircon Ce/Nd ratio can be used as a proxy of magma fO2 in replace of Ce4+/Ce3+, which is hard to be accurately determined due to the very low (often below the detection limit) zircon La and Pr contents, and the zircon La content is also susceptible to contamination by minute apatite inclusions (Lu et al., 2016; Loucks et al., 2020). Yttrium is added to the denominator of 10,000*(Eu/Eu*)/Y and (Ce/Nd)/Y because its content is commonly low in fertile magmas, caused probably by early hornblende fractionation (Lu et al., 2016, and ref. therein).
From Figures 14A–B, it is seen that the Cretaceous Huaibei-Linhuan adakites have comparable 10,000*(Eu/Eu*)/Y and Eu/Eu* ratios to those of the ore-forming granodiorite porphyries in the MLYRB (Jiurui orefield; Xu et al., 2021), as well as to many other fertile adakitic suites (Lu et al., 2016). Similar to the MLYRB ore-forming porphyries, the zircon 10,000*(Eu/Eu*)/Y and Eu/Eu* ratios of the Huaibei-Linhuan adakites are correlated positively with (Ce/Nd)/Y. This suggests that the Huaibei-Linhuan adakitic magma was not H2O-poor, as also supported by the considerable amount of hornblende present in our samples (Figure 14C).
FIGURE 14. (A) 10000*(Eu/Eu*)/Y vs. (Ce/Nd)/Y and (B) (Eu/Eu*)/Y vs. (Ce/Nd)/Y plots for zircons from the Huaibei-Linhuan intrusive rock samples (base map from Lu et al., 2016); (C) Thin-section microscopic photos showing the abundance of igneous hornblende in Taiping diorite (TP-01); (D) Zircon Ce/Ce* vs. Ce/Nd plot for the Huaibei-Linhuan intrusive rock samples, showing that Ce/Nd is a good proxy for Ce anomalies.
Therefore, we considered that the cause of Cu-Au infertility for the Huaibei-Linhuan adakites may lie in the lower fO2 of their parental magma. In Figure 14, the Huaibei-Linhuan adakites have generally lower (Ce/Nd)/Y (avg. 0.035) than their MLYRB counterpart (avg. 0.051), suggesting lower fO2 in the former. For the case of Huaibei-Linhuan, the Ce/Nd ratio can serve as good proxy for the Ce anomaly, as indicated by certain positive correlation (R2 = 0.29) between the two (Figure 14D). The more reducing nature of the Huaibei-Linhuan adakitic magma is further supported by the whole-rock V/Sc content, which serves as a proxy for the original magma fO2 (Lee et al., 2005). In the V/Sc vs. Mg discrimination plot (Figure 11B), the Huaibei-Linhuan adakites have distinctly lower V/Sc values than the Cu-Au ore-related porphyries from Quimsacocha (Ecuador) (Beate et al., 2001) and Dexing (Eastern China) (Wang et al., 2006). The V/Sc values remain constant with decreasing MgO, indicating that the V/Sc is not dependent on fractionation. Elevated magma fO2 has been suggested to be critical in forming porphyry-style Cu-Au mineralization, as it would inhibit early sulfide crystallization and allow chalcophile metals to concentrate in late-stage hydrothermal fluids (e.g., Mungall, 2002; Sillitoe, 2010; Zhang et al., 2017). Mungall et al. (2015) further proposed that sulfide liquids can ascend to the shallow crust through attaching with vapor bubbles, rather than only residing in the deep crust as traditionally believed. The reasons why the parental magma of the Huaibei-Linhuan adakites was particularly low-fO2 (even lower than that of the ore-barren adakites along the southern Tanlu fault) is unclear, but may partly be attributed to the NCC lithospheric heterogeneity, especially along/near the southern Tanlu fault (e.g., Lü, 2019; Wang et al., 2020). In Huaibei-Linhuan, assimilation of the Carboniferous-Permian coal seams (as shown by the partially-assimilated natural coke fragments in our samples; Figure 3) would likely further reduce the fO2 of the intruding dioritic-granodioritic magma, and hence further decrease its potential to form porphyry-style/-related Cu-Au mineralization.
Conclusions
1) Petrography, zircon U-Pb dating and whole-rock geochemistry suggests that the Huaibei-Linhuan intrusive rocks fall into two types, early-Early Cretaceous (ca. 130–129 Ma) high-Mg adakitic (quartz-)diorite and granodiorite, and late-Early Cretaceous (ca. 115.8 and 105.8 Ma) calc-alkaline microgabbro and dolerite.
2) The Huaibei-Linhuan high-Mg adakitic (quartz-)diorite and granodiorite are featured by low K2O/Na2O, high Sr/La, and lack of correlation between SiO2 (fractionation index) and Sr/Y and MREE/HREE. This suggests that the adakites were not derived from high-pressure garnet-amphibole fractionation or partial melting of the thickened lower crust, but from partial melting of the subducted slab. Age correlation with the regional tectonic events suggests that the slab-melting may have caused by rollback of the subducting Paleo-Pacific slab instead. Further extension and crustal delamination along the southern Tanlu fault may have led to higher-degree partial melting (and lower-degree assimilation-fractionation) that formed the gabbro-dolerite in the Huaibei-Linhuan coalfield.
3) Parental magma of the Huaibei-Linhuan high-Mg adakites was considerably less oxidized than typical porphyry Cu-Au ore-forming magmas, although their water content may have been similarly high. The low magma fO2 nature is supported by both the low zircon Ce/Nd and whole-rock V/Sc ratios. The assimilation of the Carboniferous-Permian coal-bearing sequences in the area may have further decreased the magma fO2, and thus its potential to form Cu-Au 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.
Author Contributions
SZ: investigation, funding and writing. YA: supervision and writing. CL: data interpretation and writing. HW: analysis and academic support. YL: analysis and academic support. All authors have contributed to the article and approved the submitted version.
Funding
The project is jointly funded by the Natural Science Foundation of China (Grant No. 41802244) to the corresponding author.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We sincerely thank Yuangao Sun, Shuai Zhang, and Qiang Fu for their field and logistic support. We would also like to thank Zhongqing Lin and Zheng Liu for assisting with the zircon U-Pb dating and whole-rock geochemical analysis. We are grateful to the Chief Editor Prof. David Lentz, the Handling Editor Prof. Hassan Mirnejad, and two knowledgeable reviewers for their insightful comments.
References
Aeolus Lee, C.-T., Leeman, W. P., Canil, D., and Li, Z.-X. A. (2005). Similar V/Sc systematics in MORB and arc basalts: Implications for the oxygen fugacities of their mantle source regions. J. Petrol. 46, 2313–2336. doi:10.1093/petrology/egi056
An, Y., Liu, L., Wang, M., Zheng, S., Guo, Y., Zhang, S., et al. (2018). Source and enrichment of toxic elements in coal seams around mafic intrusions: Constraints from pyrites in the Yuandian Coal Mine in Anhui, Eastern China. Minerals 8 (4), 164. doi:10.3390/min8040164
Audétat, A., Pettke, T., Heinrich, C. A., and Bodnar, R. J. (2008). Special Paper: The Composition of Magmatic-Hydrothermal Fluids in Barren and Mineralized Intrusions. Econ. Geology. 103, 877–908. doi:10.2113/gsecongeo.103.5.877
Ballard, J. R., Palin, M. J., and Campbell, I. H. (2002). Relative oxidation states of magmas inferred from Ce(IV)/Ce(III) in zircon: Application to porphyry copper deposits of northern Chile. Contrib. Mineral. Petrol. 144, 347–364. doi:10.1007/s00410–002–0402–5
Barth, A. P., and Wooden, J. L. (2010). Coupled elemental and isotopic analyses of polygenetic zircons from granitic rocks by ion microprobe, with implications for melt evolution and the sources of granitic magmas. Chem. Geology. 277, 149–159. doi:10.1016/j.chemgeo.2010.07.017
Beate, B., Monzier, M., Spikings, R., Cotten, J., Silva, J., Bourdon, E., et al. (2001). Mio-Pliocene adakite generation related to flat subduction in southern Ecuador: the Quimsacocha volcanic center. Earth Planet. Sci. Lett. 192, 561–570. doi:10.1016/S0012-821X(01)00466–6
Belousova, E., Griffin, W., O'Reilly, S. Y., and Fisher, N. (2002). Igneous zircon: trace element composition as an indicator of source rock type. Contrib. Mineral. Petrol. 143, 602–622. doi:10.1007/s00410–002–0364–7
Castillo, P. R., Janney, P. E., and Solidum, R. U. (1999). Petrology and geochemistry of Camiguin Island, southern Philippines: Insights to the source of adakites and other lavas in a complex arc setting. Contrib. Mineralogy Pet. 134, 33–51. doi:10.1007/s004100050467
Chiaradia, M., Ulianov, A., Kouzmanov, K., and Beate, B. (2012). Why large porphyry Cu deposits like high Sr/Y magmas. Sci. Repep 2, 685. doi:10.1038/srep00685
Davidson, J., Turner, S., Handley, H., Macpherson, C., and Dosseto, A. (2007). Amphibole “sponge” in arc crust. Geol 35, 787–790. doi:10.1130/G23637A.110.1130/g23637a.1
Defant, M. J., and Drummond, M. S. (1993). Mount St. Helens: Potential example of the partial melting of the subducted lithosphere in a volcanic arc. Geol 21, 547–550. doi:10.1130/0091-7613(1993)021<0547:mshpeo>2.3.co;2
Deng, J., Yang, L.-Q., Groves, D. I., Zhang, L., Qiu, K.-F., and Wang, Q.-F. (2020). An integrated mineral system model for the gold deposits of the giant Jiaodong province, eastern China. Earth-Science Rev. 208, 103274. doi:10.1016/j.earscirev.2020.103274
Dilles, J. H., Kent, A. J. R., Wooden, J. L., Tosdal, R. M., Koleszar, A., Lee, R. G., et al. (2015). Zircon Compositional Evidence for Sulfur-degassing from Ore-forming Arc Magmas. Econ. Geology. 110, 241–251. doi:10.2113/econgeo.110.1.241
Dong, Y., Liu, J., Zhang, Y., Dou, S., Li, Y., Liu, S., et al. (2020). Early Cretaceous Xiuyan adakitic granitoids in the Liaodong Peninsula, eastern China: petrogenesis and implications for lithospheric thinning of the North China Craton. Can. J. Earth Sci. 58, 50–66. doi:10.1139/cjes–2020–0033
Goldfarb, R. J., and Santosh, M. (2014). The dilemma of the Jiaodong gold deposits: Are they unique. Geosci. Front. 5, 139–153. doi:10.1016/j.gsf.2013.11.00110.1016/j.gsf.2013.11.001
Gómez-Tuena, A., Langmuir, C. H., Goldstein, S. L., Straub, S. M., and Ortega-Gutiérrez, F. (2007). Geochemical evidence for slab melting in the trans-Mexican volcanic belt. J. Petrol. 48, 537–562. doi:10.1093/petrology/egl071
Grimes, C. B., Wooden, J. L., Cheadle, M. J., and John, B. E. (2015). "Fingerprinting" tectono-magmatic provenance using trace elements in igneous zircon. Contrib. Mineral. Petrol. 170, 46. doi:10.1007/s00410–015–1199–3
Hong, L., Xu, Y., Zhang, L., Liu, Z., Xia, X., and Kuang, Y. (2020). Oxidized Late Mesozoic subcontinental lithospheric mantle beneath the eastern North China Craton: A clue to understanding cratonic destruction. Gondwana Res. 81, 230–239. doi:10.1016/j.gr.2019.11.012
Hoskin, P. W. O., and Black, L. P. (2000). Metamorphic zircon formation by solid-state recrystallization of protolith igneous zircon. J. Metamorphic Geology. 18 (4), 423–439. doi:10.1046/j.1525–1314.2000.00266.x
Jiang, X.-Y., Deng, J.-H., Luo, J.-C., Zhang, L.-P., Luo, Z.-B., Yan, H.-B., et al. (2020). Petrogenesis of Early Cretaceous adakites in Tongguanshan Cu-Au polymetallic deposit, Tongling region, Eastern China. Ore Geology. Rev. 126, 103717. doi:10.1016/j.oregeorev.2020.103717
Jicha, B. R., and Singer, B. S. (2006). Volcanic history and magmatic evolution of Seguam Island, Aleutian Island arc, Alaska. Geol. Soc. America Bull. 118, 805–822. doi:10.1130/B2586110.1130/b25861.1
Kamvong, T., Khin Zaw, K., Meffre, S., Maas, R., Stein, H., and Lai, C.-K. (2014). Adakites in the Truong Son and Loei fold belts, Thailand and Laos: Genesis and implications for geodynamics and metallogeny. Gondwana Res. 26, 165–184. doi:10.1016/j.gr.2013.06.011
Li, C., Yan, J., Wang, A.-G., Liu, J.-M., and Li, Z.-S. (2020a). Petrogenesis of Cretaceous granitoids in the Bengbu-Wuhe area, southeastern North China Craton: Implications for gold mineralization. Ore Geology. Rev. 126, 103740. doi:10.1016/j.oregeorev.2020.103740
Li, J., Wang, K.-y., Cai, W.-y., Sun, F.-y., Liu, H.-l., Fu, L.-j., et al. (2020b). Triassic gold-silver metallogenesis in Qingchengzi orefield, North China Craton: Perspective from fluid inclusions, REE and H-O-S-Pb isotope systematics. Ore Geology. Rev. 121, 103567. doi:10.1016/j.oregeorev.2020.103567
Li, S., Jahn, B.-m., Zhao, S., Dai, L., Li, X., Suo, Y., et al. (2017). Triassic southeastward subduction of North China Block to South China Block: Insights from new geological, geophysical and geochemical data. Earth-Science Rev. 166, 270–285. doi:10.1016/j.earscirev.2017.01.009
Li, S., Zhao, G., Dai, L., Liu, X., Zhou, L., Santosh, M., et al. (2012). Mesozoic basins in eastern China and their bearing on the deconstruction of the North China Craton. J. Asian Earth Sci. 47, 64–79. doi:10.1016/j.jseaes.2011.06.008
Liu, J., Liu, F.-X., Li, S.-H., and Lai, C.-K. (2019). Formation of the Baiyun gold deposit, Liaodong gold province, NE China: Constraints from zircon U-Pb age, fluid inclusion, and C-H-O-Pb-He isotopes. Ore Geol. Rev., 104. doi:10.1016/j.oregeorev.2018.12.006
Liu, J., Liu, F. X., Li, S. H., and Lai, C. K. (2020). Genesis of the Xiaotongjiapuzi gold deposit of the Liaodong gold province, Northeast China: Fluid inclusion thermometry and S-Pb-H-O-He isotope constraints. Geol. J. 55, 1023–1040. doi:10.1002/gj.3454
Liu, J., Ni, J., Chen, X., Craddock, J. P., Zheng, Y., Ji, L., et al. (2021). Early Cretaceous tectonics across the North Pacific: New insights from multiphase tectonic extension in Eastern Eurasia. Earth-Science Rev. 217, 103552. doi:10.1016/j.earscirev.2021.103552
Liu, S.-A., Li, S., Guo, S., Hou, Z., and He, Y. (2012). The Cretaceous adakitic-basaltic-granitic magma sequence on south-eastern margin of the North China Craton: Implications for lithospheric thinning mechanism. Lithos 134-135 (135), 163–178. doi:10.1016/j.lithos.2011.12.015
Liu, S.-A., Li, S., He, Y., and Huang, F. (2010a). Geochemical contrasts between early Cretaceous ore-bearing and ore-barren high-Mg adakites in central-eastern China: Implications for petrogenesis and Cu-Au mineralization. Geochimica et Cosmochimica Acta 74, 7160–7178. doi:10.1016/j.gca.2010.09.003
Liu, Y., Gao, S., Hu, Z., Gao, C., Zong, K., and Wang, D. (2010b). Continental and Oceanic Crust Recycling-induced Melt-Peridotite Interactions in the Trans-North China Orogen: U-Pb Dating, Hf Isotopes and Trace Elements in Zircons from Mantle Xenoliths. J. Pet. 51, 537–571. doi:10.1093/petrology/egp082
Liu, Y., Hu, Z., Gao, S., Günther, D., Xu, J., Gao, C., et al. (2008). In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard. Chem. Geology. 257, 34–43. doi:10.1016/j.chemgeo.2008.08.004
Loucks, R. R., Fiorentini, M. L., and Henríquez, G. J. (2020). New Magmatic Oxybarometer Using Trace Elements in Zircon. J. Petrol., 61. doi:10.1093/petrology/egaa034
Lü, Q., Meng, G., Zhang, K., Liu, Z., Yan, J., Shi, D., et al. (2021). The lithospheric architecture of the Lower Yangtze Metallogenic Belt, East China: Insights into an extensive Fe-Cu mineral system. Ore Geology. Rev. 132, 103989. doi:10.1016/j.oregeorev.2021.103989
Lu, Y.-J., Loucks, R. R., Fiorentini, M., McCuaig, T. C., Evans, N. J., Yang, Z.-M., et al. (2016). Zircon Compositions as a Pathfinder for Porphyry Cu ± Mo ± Au Deposits. Tectonics Metallog. Tethyan Orog. Belt 19, 0. doi:10.5382/SP.19.13
Lü, Z. (2019). Uppermost mantle velocity and anisotropy structure beneath the North China Craton and its adjacent regions. Tectonophysics 754, 45–55. doi:10.1016/j.tecto.2019.01.014
Ludwig, K. R. (2003). ISOPLOT 3.00: A Geochronological Toolkit for Microsoft Excel. Berkeley: Berkeley Geochronology CenterCalifornia.
Macpherson, C. G. (2008). Lithosphere erosion and crustal growth in subduction zones: Insights from initiation of the nascent East Philippine Arc. Geol 36, 311–314. doi:10.1130/G24412A.1
Maniar, P. D., and Piccoli, P. M. (1989). Tectonic discrimination of granitoids. Geol. Soc. Am. Bull. 101 (5), 635–643. doi:10.1130/0016-7606(1989)101<0635:tdog>2.3.co;2
Mao, J., Wang, Y., Lehmann, B., Yu, J., Du, A., Mei, Y., et al. (2006). Molybdenite Re-Os and albite 40Ar/39Ar dating of Cu-Au-Mo and magnetite porphyry systems in the Yangtze River valley and metallogenic implications. Ore Geology. Rev. 29, 307–324. doi:10.1016/j.oregeorev.2005.11.001
Middlemost, E. A. K. (1994). Naming materials in the magma/igneous rock system. Earth-Science Rev. 37, 215–224. doi:10.1016/0012–8252(94)90029–9
Mori, L., Gómez-Tuena, A., Cai, Y., and Goldstein, S. L. (2007). Effects of prolonged flat subduction on the Miocene magmatic record of the central Trans-Mexican Volcanic Belt. Chem. Geology. 244, 452–473. doi:10.1016/j.chemgeo.2007.07.002
Mungall, J. E., Brenan, J. M., Godel, B., Barnes, S. J., and Gaillard, F. (2015). Transport of metals and sulphur in magmas by flotation of sulphide melt on vapour bubbles. Nat. Geosci 8, 216–219. doi:10.1038/ngeo2373
Mungall, J. E. (2002). Roasting the mantle: Slab melting and the genesis of major Au and Au-rich Cu deposits. Geol 30, 915–918. doi:10.1130/0091-7613(2002)030<0915:rtmsma>2.0.co;2
Pearce, J. (1996). A user's Guide to basalt discrimination diagrams. Geol. Assoc. Can. Short Course Notes 12, 79–113. doi:10.1111/j.1438–8677.1980.tb03374.x10.1093/oq/12.3.83
Peccerillo, A., and Taylor, S. R. (1976). Geochemistry of Eocene Calc-Alkaline Volcanic Rocks from the Kastamonu Area, Northern Turkey. Contr. Mineral. Petrol. 58, 63–81. doi:10.1007/BF00384745
Ping, X., Zheng, J., Xiong, Q., Griffin, W. L., Yu, C., and Su, Y. (2019). Downward rejuvenation of the continental lower crust beneath the southeastern North China Craton. Tectonophysics 750, 213–228. doi:10.1016/j.tecto.2018.11.012
Pirajno, F., and Zhou, T. (2015). Intracontinental Porphyry and Porphyry-Skarn Mineral Systems in Eastern China: Scrutiny of a Special Case "Made-in-China". Econ. Geology. 110, 603–629. doi:10.2113/econgeo.110.3.603
Reich, M., Parada, M. A., Palacios, C., Dietrich, A., Schultz, F., and Lehmann, B. (2003). Adakite-like signature of Late Miocene intrusions at the Los Pelambres giant porphyry copper deposit in the Andes of central Chile: metallogenic implications. Mineralium Deposita 38, 876–885. doi:10.1007/s00126–003–0369–9
Reubi, O., and Nicholls, I. A. (2004). Magmatic evolution at Batur volcanic field, Bali, Indonesia: petrological evidence for polybaric fractional crystallization and implications for caldera-forming eruptions. J. Volcanology Geothermal Res. 138, 345–369. doi:10.1016/j.jvolgeores.2004.07.009
Richards, J. P. (2011). High Sr/Y Arc Magmas And Porphyry Cu Mo Au Deposits: Just Add Water. Econ. Geology. 106, 1075–1081. doi:10.2113/econgeo.106.7.1075
Richards, J. P., and Kerrich, R. (2007). Special Paper: Adakite-Like Rocks: Their Diverse Origins and Questionable Role in Metallogenesis. Econ. Geology. 102 (4), 537–576. doi:10.2113/gsecongeo.102.4.537
Richards, J. P., Spell, T., Rameh, E., Razique, A., and Fletcher, T. (2012). High Sr/Y Magmas Reflect Arc Maturity, High Magmatic Water Content, and Porphyry Cu ± Mo ± Au Potential: Examples from the Tethyan Arcs of Central and Eastern Iran and Western Pakistan. Econ. Geol. 107, 295–332. doi:10.2113/econgeo.107.2.295
Sillitoe, R. H. (2010). Porphyry copper systems. Econ. Geology. 105, 3–41. doi:10.2113/gsecongeo.105.1.3
Sun, P., Guo, P., and Niu, Y. (2021). Eastern China continental lithosphere thinning is a consequence of paleo-Pacific plate subduction: A review and new perspectives. Earth-Science Rev. 218, 103680. doi:10.1016/j.earscirev.2021.103680
Sun, S.-s., and McDonough, W. F. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geol. Soc. Lond. Spec. Publications 42, 313–345. doi:10.1144/GSL.SP.1989.042.01.19
Sun, W., Xie, Z., Chen, J., and Zhang, X. (2003). Os-Os dating of copper and molybdenum deposits along the Middle and Lower reaches of the Yangtze River, China. Econ. Geology. 98, 175–180. 0361–0128/01/3323/175–610.2113/98.1.175.
Wang, K., Chen, L., Xiong, X., Yan, Z., and Xie, R. (2020). The role of lithospheric heterogeneities in continental rifting: Implications for rift diversity in the North China Craton. J. Geodynamics 139, 101765. doi:10.1016/j.jog.2020.101765
Wang, Q., Xu, J.-F., Jian, P., Bao, Z.-W., Zhao, Z.-H., Li, C.-F., et al. (2006). Petrogenesis of Adakitic Porphyries in an Extensional Tectonic Setting, Dexing, South China: Implications for the Genesis of Porphyry Copper Mineralization. J. Petrol. 47, 119–144. doi:10.1093/petrology/egi070
Xia, Q.-K., Hao, Y.-T., Liu, S.-C., Gu, X.-Y., and Feng, M. (2013). Water contents of the Cenozoic lithospheric mantle beneath the western part of the North China Craton: Peridotite xenolith constraints. Gondwana Res. 23 (1), 108–118. doi:10.1016/j.gr.2012.01.010
Xie, G.-Q., Mao, J.-W., Li, R.-L., Qü, W.-J., Pirajno, F., and Du, A.-D. (2007). Re-Os molybdenite and Ar-Ar phlogopite dating of Cu-Fe-Au-Mo (W) deposits in southeastern Hubei, China. Mineralogy Pet. 90, 249–270. doi:10.1007/s00710–006–0176–y
Xu, W.-L., Gao, S., Yang, D.-B., Pei, F.-P., and Wang, Q.-H. (2009). Geochemistry of eclogite xenoliths in Mesozoic adakitic rocks from Xuzhou-Suzhou area in central China and their tectonic implications. Lithos 107 (3), 269–280. doi:10.1016/j.lithos.2008.11.004
Xu, W., Gao, S., Wang, Q., Wang, D., and Liu, Y. (2006). Mesozoic crustal thickening of the eastern North China craton: Evidence from eclogite xenoliths and petrologic implications. Geol 34 (9), 721–724. doi:10.1130/G22551.1
Xu, W. L., Wang, Q. H., Liu, X. C., Wang, D. Y., and Guo, J. H. (2004). Chronology and Sources of Mesozoic Intrusive Complexes in the Xuzhou-Huainan Region, Central China: Constraints from SHRIMP Zircon U-Pb Dating. Acta Geologica Sinica (English Edition) 78 (1), 96–106. doi:10.1111/j.1755–6724.2004.tb00679.x
Xu, Y.-M., Jiang, S.-Y., and Zhu, J.-X. (2021). Factors controlling the formation of large porphyry Cu deposits: A case study from the Jiurui ore district of Middle-Lower Yangtze River Metallogenic Belt using in situ zircon and apatite chemistry from syn-mineralization intrusions. Ore Geology. Rev. 133, 104082. doi:10.1016/j.oregeorev.2021.104082
Yang, D.-B., Xu, W.-L., Wang, Q.-H., and Pei, F.-P. (2010). Chronology and geochemistry of Mesozoic granitoids in the Bengbu area, central China: Constraints on the tectonic evolution of the eastern North China Craton. Lithos 114 (1-2), 200–216. doi:10.1016/j.lithos.2009.08.009
Yang, J.-H., Wu, F.-Y., Wilde, S. A., Belousova, E., and Griffin, W. L. (2008). Mesozoic decratonization of the North China block. Geol 36 (6), 467–470. doi:10.1130/G24518A.1
Yangliu Coal Industry Co. Ltd (2017). Resources verification report on the reserves of Yangliu Coal Mine in Suixi County, Anhui Province. Huaibei, Anhui, 62–65.
Yao, Y., and Liu, D. (2012). Effects of igneous intrusions on coal petrology, pore-fracture and coalbed methane characteristics in Hongyang, Handan and Huaibei coalfields, North China. Int. J. Coal Geology. 96-97, 72–81. doi:10.1016/j.coal.2012.03.007
Ye, G., Liu, C., Luo, X., Jin, S., Wei, W., Dong, H., et al. (2021). Dynamical significance of the Tanlu Fault Zone in the destruction of the North China Craton: The evidence provided by the three-dimensional Magnetotelluric array study. Tectonophysics 813, 228910. doi:10.1016/j.tecto.2021.228910
Zhang, C.-c., Sun, W.-d., Wang, J.-t., Zhang, L.-p., Sun, S.-j., and Wu, K. (2017). Oxygen fugacity and porphyry mineralization: A zircon perspective of Dexing porphyry Cu deposit, China. Geochimica et Cosmochimica Acta 206, 343–363. doi:10.1016/j.gca.2017.03.013
Zhang, C., Ma, C., Holtz, F., Koepke, J., Wolff, P. E., and Berndt, J. (2013). Mineralogical and geochemical constraints on contribution of magma mixing and fractional crystallization to high-Mg adakite-like diorites in eastern Dabie orogen, East China. Lithos 172-173, 118–138. doi:10.1016/j.lithos.2013.04.011
Zhang, J., Qu, J., Zhang, B., Zhao, H., Niu, P., Zhao, S., et al. (2020). Mesozoic intraplate deformation of the central North China Craton: Mechanism and tectonic setting. J. Asian Earth Sci. 192, 104269. doi:10.1016/j.jseaes.2020.104269
Zhang, Y., Cheng, J., Tian, J., Pan, J., Sun, S., Zhang, L., et al. (2019a). Texture and trace element geochemistry of quartz in skarn system: Perspective from Jiguanzui Cu-Au skarn deposit, Eastern China. Ore Geol. Rev., 109. doi:10.1016/j.oregeorev.2019.05.007
Zhang, Z., Wang, Y., Li, D., and Lai, C. (2019b). Lithospheric architecture and metallogenesis in Liaodong peninsula, North China Craton: Insights from zircon Hf-Nd isotope mapping. Minerals 9, 179. doi:10.3390/min9030179
Zhao, L., Allen, R. M., Zheng, T., and Hung, S.-H. (2009). Reactivation of an Archean craton: Constraints from P- and S-wave tomography in North China. Geophys. Res. Lett. 36 (17), 367–389. doi:10.1029/2009GL039781
Zhao, M., An, Y., Wang, M., Ding, M., and Lai, C. (2019). New Genesis of Natural Coke around Magmatic Intrusion at the Shitai Coalmine of Huaibei City, North China. Acta Geologica Sinica ‐ English Edition 93 (4), 1158–1159. doi:10.1111/1755-6724.13827
Zhao, T., Zhu, G., Lin, S., and Wang, H. (2016). Indentation-induced tearing of a subducting continent: Evidence from the Tan-Lu Fault Zone, East China. Earth-Science Rev. 152, 14–36. doi:10.1016/j.earscirev.2015.11.003
Zhao, Z., Liang, S., Santosh, M., and Wei, J. (2020). Lithospheric extension associated with slab rollback: Insights from early Cretaceous magmatism in the southern segment of Tan-Lu fault zone, central-eastern China. Lithos 362-363, 105487. doi:10.1016/j.lithos.2020.105487
Zhou, H., Shang, D. F., Chan, S. W., and Chen, J. (2019). Zircon U-Pb dating, petrogenesis and geological significance of the Banjing pluton in Xuhuai area. J. Mineral. Petrol. 39, 9–16. (in Chinese with English abstract).
Zhou, T., Wang, S., Fan, Y., Yuan, F., Zhang, D., and White, N. C. (2015). A review of the intracontinental porphyry deposits in the Middle-Lower Yangtze River Valley metallogenic belt, Eastern China. Ore Geology. Rev. 65, 433–456. doi:10.1016/j.oregeorev.2014.10.002
Zhu, G., Liu, G. S., Niu, M. L., Xie, C. L., Wang, Y. S., and Xiang, B. (2009). Syn-collisional transform faulting of the Tan-Lu fault zone, East China. Int. J. Earth Sci. (Geol Rundsch) 98, 135–155. doi:10.1007/s00531–007–0225–8
Keywords: high-Mg adakites, North China Craton, Paleo-Pacific subduction, Cretaceous, zircon U-Pb dating
Citation: Zheng S, An Y, Lai C, Wang H and Li Y (2021) Genesis of High-Mg Adakites in the Southeastern Margin of North China Craton: Geochemical and U-Pb Geochronological Perspectives. Front. Earth Sci. 9:731233. doi: 10.3389/feart.2021.731233
Received: 26 June 2021; Accepted: 21 September 2021;
Published: 11 October 2021.
Edited by:
Hassan Mirnejad, Miami University, United StatesReviewed by:
Wayne Powell, Brooklyn College (CUNY), United StatesMohsen Rezaei, Shahid Chamran University of Ahvaz, Iran
Copyright © 2021 Zheng, An, Lai, Wang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanfei An, YW55QGFodS5lZHUuY24=