Geochronology and Petrogenesis of Early Pleistocene Dikes in the Changbai Mountain Volcanic Field (NE China) Based on Geochemistry and Sr-Nd-Pb-Hf Isotopic Compositions

Changbai Mountains intraplate volcanism (NE China) developed above the 500 km deep stagnant portion of the Pacific slab in the last 20 Ma. The more recent activity includes a shield-forming stage (2.8–0.3 Ma), the Tianchi cone construction stage (1.5–0.01 Ma), and a caldera-forming stage (0.2 Ma-present). Detailed studies on the petrogenesis of the volcanic products between the first two stages and the possible role of geodynamics and local tectonics in controlling the volcanism, however, are lacking. Here, we present structural and whole-rock geochemical and zircon Hf isotopic data on Pleistocene dikes of the Changbai Mountains at the transition from the shield-forming to the Tianchi stage with the aim to constrain their age and the source(s) of their parental magma. The dikes represent the shallower feeding system of monogenetic cones and have a NW-SE strike, which is also the preferred strike of the major fault affecting the area and along which the Changbai Mountains monogenetic scoria cones align. The dikes have a potassic affinity and a trachybasaltic composition. Their zircon U–Pb age is 1.19–1.20 Ma (Calabrian). The trachybasalts are enriched in Rb, Ba, Th, U, Nb, Ta, K, Pb, and LREE and slightly depleted in Sr, Zr, Hf, Ti, and HREE with a weak negative Eu/Eu* (δEu = 0.96–0.97). Trace elements and isotopic compositions are compatible with an OIB-type source with an EMI signature. The calculated (87Sr/86Sr)i (=0.705165–0.705324), (143Nd/144Nd)i (= 0.512552–0.512607, εNd(t) = −0.58 to −1.65), and Hf model ages (TDM2) of 1768–1562 Ma suggest that the trachybasaltic dikes were contaminated by a Mesoproterozoic, relatively basic lower crust. The source of the Calabrian trachybasalts consists of asthenospheric melts modified by a subcontinental lithospheric mantle. These melts upwell from depth and stop at the crust-mantle interface where underplating processes favor the assimilation of ancient lower crust material. During the ascent to the surface along deep-seated crustal discontinuities, these magmas weakly differentiate.


HIGHLIGHTS INTRODUCTION
The eastern China block includes the eastern margin of the Central Asian orogenic belt (the Xing'an-Mongolia orogenic belt), the North China Craton (NCC), the Sulu-Dabie high pressure orogeny belt, the Yangtze Craton, and the southeast coastal orogenic belt of China, and it is characterized by Mesozoic volcanic and intrusive rocks ( Figure 1A). The later Cenozoic volcanic rocks mainly consist of basalts emitted in an intraplate setting (e.g., Wu et al., 2011;Zhang and Luo, 2011;Deng et al., 2017) associated to metallogenic and hydrothermal deposits (Sun et al., 2012(Sun et al., , 2013Deng et al., 2017). Changbai Mountains, which consist of basaltic shields, monogenetic vents, and central-type volcanoes, are located above the 500 km deep stagnant slab of the westward subducting Pacific plate ( Figure 1A). Three main geodynamic models have been suggested to explain the Changbai Mountain volcanism: an upwelling mantle plume (Mao et al., 1999), an intraplate continental rifting (Gilder et al., 1991), and the Pacific plate subduction (e.g., Mao et al., 2013;Ouyang et al., 2014;Zhang et al., 2014;Xu W. L. et al., 2020;Xu Z.-t. et al., 2020). Previous studies on the Changbai Mountains area mainly focus on the Holocene silicic magmatism (e.g., Wei et al., 2007;Liu et al., 2009;Xu J. D. et al., 2013a). The widely distributed early Pleistocene rocks are less known. The origin of the Pleistocene magmatism is still debated and different hypotheses have been proposed: 1) low degrees of partial melting of the asthenospheric mantle (Xu, et al., 2020a) possibly contaminated by the addition of crustal material (Sun et al., 2015); 2) partial melting of the subducted slab (Ma et al., 2015); 3) involvement of DM and EM1 end-members (Basu et al., 1991) or OIB and EM1 end-members (Xu et al., 2012); 4) an asthenospheric source metasomatized by recycled subducted silicate sediments or a subcontinental lithospheric mantle (EM1-like) and by carbonated eclogites with a MORB signature (Choi et al., 2020); and 5) an EM1-like source originated from the contamination of an OIB asthenospheric mantle by Paleoproterozoic lithosphere (Lee et al., 2021). Therefore, questions regarding the nature of the mantle and possible crustal sources of the Changbaishan Mountain volcanism and of the dynamics of the Pacific plate subduction remain open. In this study, we carried out systematic LA-MC-ICP-MS zircon U-Pb geochronology as well as whole-rock geochemistry, Sr-Nd-Pb isotopic analyses, and zircon Hf isotopic analyses to investigate the petrogenesis of the dikes related to the early Pleistocene monogenetic activity in area. The results allow us to improve our understanding of the petrogenesis and tectonic setting of the Changbai Mountains magmatism. New data and constraints on the relative role of the asthenosphere, ancient crust, and stagnant slab in the evolution of intraplate volcanism are provided.

GEOLOGICAL SETTING
The NE China magmatism mainly developed during Mesozoic and Cenozoic (Ge et al., 2003;Wu et al., 2005;Sun et al., 2007) eras. Because of continental rifting processes (Gilder et al., 1996;Zhou and Li, 2000;Zhu et al., 2011), the Changbai Mountain volcanic field (CMVF) is bounded by two NE-NW striking regional fault systems, the Dunhua-Mishan faults to the northeast and the Xar Morron-Changchun faults to the southeast although faults with a NW-SE strike are also locally present ( Figure 1B). The origin of the latter faults is unclear, and they may represent the response of the NE China terrains to different rates of rifting associated to the Pacific slab roll-back. The basement of CVMF consists of Precambrian metamorphic rocks including Archean tonalite, trondhjemite, granodiorite (2.4-2.5 Ga) and Proterozoic quartzite, schists, and granulite assemblages (1.7-1.9 Ga) (Jin and Zhang, 1994;Wang et al., 2003). The CMVF volcanic rocks were emitted by three stratovolcanoes located at the China/North Korea boundary: the Tianchi caldera (last eruption in 1903) and the now inactive Wangtian'e in China and Namphothe edifices in North Korea. Several monogenetic scoria cones surround these three large stratovolcanoes. They roughly align along a NW-SE direction (inset in Figure 1). This direction is the same of the main fault affecting the Tianchi caldera (Baishan-Jince fault in Figure 1). After a basaltic plateau forming phase developed between 20 and 3 Ma, the volcanism in the Tianchi area can be roughly divided in three stages: 1) shield-forming stage (2.8-0.3 Ma), 2) the Tianchi cone construction stage (1.5-0.01 Ma), and 3) caldera-forming, explosive eruption stage (0.2 Ma-present) (Fan et al., 1999;Liu et al., 2015;Takeshi et al., 2020). Monogenetic volcanism mainly developed during the shield-forming stage, although some vents formed also during the cone construction stage. The early two stages were characterized by the emission of basaltic lava flows and scorias, while trachytes and rhyolites were emitted as lava flows and pyroclastic fall and flows, and these intermediate-acid magmatism mainly occurred in the last stage of basaltic magmatism.
The study area is located in the northern sector of the Tianchi area ( Figure 1B). Here, the basement rocks include gneiss of the Paleoarchean Sidaolazi formation, amphibolite of the Yangjiadian formation, and quartz schist of the Mesoproterozoic Huashan formation. The volcanic rocks consist of Pleistocene and Pliocene basaltic lava flows and minor scorias (Figure 2). The products studied here are two dikes representing the shallowest feeding system of two neighbors, partly dismantled cones constituted by oxidized scorias and ash ( Figure 2). Dike 1 strikes NNW-SSE; it is15-20 m long and 0.4-1.8 m thick. Dike 2 strikes NW-SE with dip 80°-85°; it is about 30 m long and 0.7-5.2 m thick with some lateral digitations having length of less than 2 m ( Figures  3A,B). These dikes have 5-15 mm wide vesicles testifying the degassing of magma near the surface ( Figures 3C,D).

SAMPLING AND ANALYTICAL TECHNIQUES
The sampling locations are shown in Figure 2. Unaltered, fresh samples were selected for a combined study involving whole-rock major and trace elements analyses, cathodoluminescence (CL) imaging, zircon U-Pb ages, and whole-rock Sr-Nd-Pb isotopes and zircon Hf isotopic analyses.

Sample Descriptions
The dikes are dark gray in color, massive to poorly vesiculated, and show a trachytic texture with plagioclases phenocrysts representing 80-85 vol% of the total minerals ( Figures 3E,F). Clinopyroxene and biotite account for 5-10 vol% of the phenocrysts.
The groundmass includes fine-grained plagioclase, pyroxene, amphibole, zircon ( Figures 3G,H), apatite, and magnetite microlites. No significant differences in mineral assemblage, texture, and crystal content have been recognized between the two dikes. The Japan trench is reported as red solid line with triangles and the depth of the Pacific slab is reported as dashed lines; the number indicates the depth in km (modified from Zhang et al., 2018;Chen et al., 2021;Fan et al., 2021;. (B) Geological map of the Changbai Mountains volcanic field showing distribution of polygenetic (Changbaishan (Tianchi), Wangtian'e, and Namphothe and monogenetic (gray dots) volcanoes. Geological map legend: 1, pre-shield basalts; 2, shield-forming basalts; 3, post-shield basalts; 4, intermediate to silicic volcanic rocks (mainly trachyte with minor trachyandesite); 5, intermediate to silicic volcanic rocks (mainly trachyte with minor trachyandesite); 6, early rhyolite; 7, rhyolite; 8, pyroclastic deposits (primarily trachytic and rhyolitic in composition) erupted by caldera-forming eruptions of Tianchi. The inset on the top right reports the spatial density of monogenetic vents.

LA-MC-ICP-MS Zircon U-Pb Dating
Zircon was mounted in epoxy discs, polished to expose the grains, cleaned ultrasonically in ultrapure water, and then cleaned again prior to the analysis using AR grade methanol. Pre-ablation was conducted for each spot analysis using 5 laser shots (∼0.3 μm in depth) to remove potential surface contamination. The analysis was performed using 30 μm diameter spot at 5 Hz and a fluence of 2 J/cm 2 . Zircon U-Pb dating was conducted using LA-MC-ICP-MS (laser ablation multiple collector inductively coupled plasma mass spectrometry) in Nanjing Hong Chuang Geological Exploration Technology Service Company Limited (NHGETSCL). The Resolution SE model laser ablation system (Applied Spectra, United States) was equipped with ATLEX 300 excimer laser and a two volume S155 ablation cell. The laser ablation system was coupled to a Neptune Plus MC-ICP-MS (Wrexham, Wales, United Kingdom). Detailed tuning parameters can be seen in Thompson et al. (2018). Iolite software package was used for data reduction (Paton et al., 2010). Zircon 91500 and GJ-1 were used as primary and secondary reference materials, respectively. 91500 was analyzed twice and GJ-1 analyzed once every 10-12 analyses of the sample.
The age of GJ-1 (604 ± 6 Ma, 1 σ) measured during the experiment was consistent with the recommended values within the uncertainty range, and the preferred value of 91500 is in Wiedenbeck et al. (1995); the 207 Pb/ 206 Pb and 206 Pb/ 238 U age of 91500 is 1065.4 and 1062.4 Ma. Based on the decay coefficient, Th/U ratio, and the error of the measurements in young zircons, we corrected the age loss caused by 203 Th depletion (Zou et al., 2010). The LA-MC-ICP-MS U-Pb isotopic data are listed in Table 1.

Zircon Lu-Hf Isotopic Analyses
Experiments of in situ Hf isotope ratio analysis were conducted using the Teledyne Cetac Technologies Analyte Excite laser ablation system (Bozeman, Montana, United States) and the Nu Instruments Nu Plasma II MC-ICP-MS (Wrexham, Wales, United Kingdom) that were hosted at the State Key Laboratory for Mineral Deposits Research, Nanjing University. The 193 nm ArF excimer laser, homogenized by a set of beam delivery systems, was focused on zircon surface with fluence of 6.0 J/cm 2 . Ablation protocol employed a spot diameter of 50 um at 8 Hz repetition rate for 40 s. Helium was applied as carrier gas to efficiently transport aerosol to MC-ICP-MS; the results of the zircon Hf isotopic analyses are given in Table 2. The precision of sample test is 0.000005-0.000012 (0.015-0.04‰, 2 s), and the accuracy is better than 0.000 012 (∼0.04‰).

Major and Trace Element Determinations
The geochemical sample analyses were performed at the State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences. The major element oxides were analyzed by a Panalytical PW4400 X-ray fluorescence spectrometer (AxiosmAX, Holland). The detailed sample-digesting procedure was as follows: 1) sample powder (200 mesh) were placed in an oven at 105°C for drying of 12 h. 2) ∼1.0g dried sample was accurately weighted and placed in the ceramic crucible and then heated in a muffle furnace at 1000°C for 2 h. After cooling to 400°C, this sample was placed in the drying vessel and weighted again in order to calculate the loss on ignition (LOI). 3) 0.6 g sample powder was mixed with 6.0 g cosolvent and 0.3 g oxidant in a Pt crucible, which was placed in the furnace at 1150°C for 14 min. Then, the sample was Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 729905 6 quenched with air for 1 min to produce flat discs on the 5fire brick for the XRF analyses. Trace elements and rare earth elements (REEs) were analyzed by Perkin-Elmer ELAN 300D ICP-MS (United States). The detailed sample-digestion procedure refers to Sun et al. (2013). The analytical precisions for major elements and trace elements were better than 3 and 5%, respectively, and the major and trace element compositions of the trachybasalt dikes (CBS-5a, CBS-5b, CBS-5c, CBS-5a-1, CBS-5b-1, CBS-6a, and CBS-6b) are given in Table 3.

Sr-Nd-Pb Isotopic Analyses
Whole-rock Pb isotopic measurements were conducted by a MAT-261 thermal ionization mass spectrometer. All chemical preparations were performed on class 100 work benches within a class 1000 over-pressured clean laboratory. Sample digestion is as follows: 1) sample powder (200 mesh) was placed in an oven at 105°C for drying of 12 h; 2) 50-200 mg sample powder was accurately weighed and placed in a Teflon bomb; 3) 1-3 ml HNO 3 and 1-3 ml HF were added into the Teflon bomb; 4) Teflon bomb was put in a stainless steel pressure jacket and heated to 190°C in an oven for >24 h; 5) after cooling, the Teflon bomb was opened and placed on a hotplate at 140°C and evaporated to incipient dryness, and then 1 ml HNO3 was added and evaporated to dryness again; and 6) the sample was dissolved in 1.0 ml of 1.0 M HBr. After centrifugation, the supernatant solution was loaded into an ion-exchange column packed with AG resin. After complete draining of the sample solution, columns were rinsed with 1.0 M HBr to remove undesirable matrix elements. Finally, the Pb fraction was eluted using 6.0 M HCl and gently evaporated to dryness prior to massspectrometric measurement. The sample analyses were performed at the Analytical Laboratory of the Beijing Research Institute of Uranium Geology (BRIUG), the precision of 20x Pb/ 204 Pb is (0.002-0.025%, 1σ), and the accuracy is better than 0.03%. The results of the wholerock Pb isotopic analyses are given in Table 4.
The results of Sr-Nd isotope analyses are reported in Table 5. All chemical preparations were performed on class 100 work benches within a class 1000 over-pressured clean laboratory; the experimental process refers to Pan et al. (2020).
The REE solution was used to separate the Nd fraction by the Nd-column method. Sr and Nd isotopes analyses were performed on a Nu Plasma II MC-ICP-MS (Wrexham, Wales, UK) at the Analytical Laboratory of the BRIUG, Beijing, China. The precision of the test sample 0.000005-0.00002 (0.01-0.05‰, 1σ), better than 0.000025 (∼0.05‰). Raw data of isotopic ratios were corrected for mass fractionation by normalizing to 86 Sr/ 88 Sr 0.1194 for Sr, 146 Nd/ 144 Nd 0.7219 for Nd, and 205 Tl/ 203 Tl 2.3885 for Pb with exponential law. International isotopic standards (NIST SRM 987 for Sr, JNdi-1 for Nd, and NIST SRM 981 for Pb) were periodically analyzed to correct instrumental drift. Geochemical reference materials of USGS BCR-2, BHVO-2, AVG-2, and RGM-2 were treated as quality control.

Zircon U-Pb Dating
Zircons from the trachybasalt dikes (sample CBS-5a and CBS-6a) are light gray with no inclusions. They are euhedral-subhedral columnar crystals with size of 50-200 μm (aspect ratio of 1:1 to 3:1). In cathodoluminescence (CL) images, the zircons exhibit low, homogeneous luminescence and occasionally display oscillatory zoning (Figure 4). Th/U is 0.8-2.5, a range of values consistent with a magmatic origin (Hoskin and Schaltegger, 2003). The 25 analyses of zircon grains from CBS-5a and CBS-6a give weighted mean 206 Pb/ 238 U ages of 1.19 ± 0.0057 Ma (MSWD 2.00, n 015) and 1.20 ± 0.0048 Ma (MSWD 0.89, n 10), respectively (Figures  5B,D). Because the magmatism in the Tianchi area occurred in the late Cenozoic, the content of common Th, U, and Pb in the Th-U-Pb system of zircons may be considered high, and the content of 235 U decays to 207 Pb is low. Therefore, 206 Pb/ 208 U is selected as the weighted average age. Tera-Wasserburg concordia diagrams were used to calculate the concordant age (1.19 ± 0.0030 Ma and 1.20 ± 0.0038, respectively, Figures 5A,C). This age is taken to represent the emplacement age of the studied dikes.

Sr-Nd Isotope Composition
Five samples of the trachybasalt dikes were analyzed for their whole-rock Sr-Nd isotopic composition. They have evolved radiogenic isotope compositions with ( 87 Sr/ 86 Sr) i 0.705165-0.705324 and ( 143 Nd/ 144 Nd) i 0.512552-0.512607 corresponding to εNd (t) values of −0.58 to −1.65. According to the data reported in Figure 8C, the Sr-Nd isotope composition of the trachybasaltic dikes approaches that of the coeval basaltic rocks in the CMVF region.    Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 729905 11

Timing and Tectonic Control of Volcanism at CMVF
Application of U-Pb disequilibrium dating of zircons with contrasting stabilities and compositions can provide a unique perspective to study the evolution of magma within the framework of absolute time (Wu et al., 2004). All dike outcrops in the field consist of hypabyssal basalts ( Figures  3A,B). All dated zircons from the dikes are euhedral-subhedral columnar crystals and display oscillatory zoning; their Th/U is 0.8-2.5, a value consistent with a magmatic origin. Although the zircon crystals of the two dikes could be interpreted as antecrysts, however, their petrographic features do not support such hypothesis. The growth of zircons in poorly evolved magmas may be due to different processes including (Borisova et al., 2020) the partial melting of hydrated peridotite, the felsic crustal involvement/recycling or zircon sealing, total zircon dissolution, and reprecipitation and/or a contribution of metasomatic Si-F-rich fluid. Our data do not allow us to better constrain the origin of the CMVF zircons studied here. As preliminary hypothesis, crustal recycling (see Discussion) or metasomatic processes in a hypabyssal environment could be responsible for the formation of the zircon crystals of CMVF dikes. Field observations suggest that dikes were formed by basaltic magma hypabyssal process ( Figures 3A,B), consistent with the euhedral-subhedral columnar morphology and oscillatory zoning of all the dated zircon crystals. Single grain LA-MC-ICP-MS zircon U-Pb concordance ages of 1.19-1.20 Ma (n 25) and 206 Pb/ 238 U ages (1.18-1.20 Ma) are here obtained from the trachybasaltic dikes and confirm that the age is 1.18-1.20 Ma. We conclude that studied dikes emplaced in the Calabrian stage of Pleistocene. As previously reported, the last three main stages of the Tianchi   (Boynton, 1984;Sun and McDonough, 1989). Other data sources are as for Figure 7C; 143 Nd/ 144 Nd versus 87 Sr/ 86 Sr plot of the CMVF trachybasalts. The field of the Changbaishan basalts is from Zhang et al. (2018). Reference ranges for DM, enriched mantle (EM1 and EM2), MORB, and OIB are from Zindler and Hart (1986 Zou et al., 2021). Therefore, our new geochronological data show that the trachybasalt dikes emplaced between the shield and cone forming stages ( Figures  9A-C). The trachybasalt dikes, which are coeval with the less evolved basaltic magmatism, evidence an early phase of evolution processes within the CMVF magmatic system, at least in the Tianchi norther sector. Our field data better constrain the role of regional tectonics in controlling the CMVF volcanism during Calabrian. In fact, the NW-SE strike of dikes mirrors that of the major structural discontinuity affecting the CMVF area, i.e., the Baishan-Jince fault, and the preferred alignment of the monogenetic volcanoes ( Figure 1). As a result, the CVMF Calabrian magmatism developed along a deep discontinuity probably affecting the whole crust along which basaltic magma upwelled and partly differentiated. The SiO 2 content and the values of A/CNK, MgO, Fe 2 O 3 , TiO 2 , P 2 O 5 , and V of the dikes confirm that they represent poorly evolved products of mantle melts.

Petrogenesis
As known, intraplate magmatism can be due to different processes including 1) asthenosphere upwelling carried by mantle plume associated or not to a hot spot, 2) lithosphere thinning, continental stretching, and consequent upwelling of the lithospheric mantle, 3) fluid release from a subducting plate and mantle metasomatism, and 4) crust-mantle interaction (Defant and Drummond., 1990;Puffer, 2001;Sisson et al., 2005). The Sr, Nd, Pb, and Hf isotopic rations ( Figure 7C; Figures Figure 8B) are consistent with an OIBtype mantle source ( Figure 10B). In particular, the values of the Sr and Nd isotopic ratios ( Figure 7C) overlap those of the CVMF basaltic rocks, for which an OIB source has been recognized (Zhang et al., 2018;Choi et al., 2020). However, the trace element pattern indicates that the source is more complex. The (Na  Wu et al., 2002;Sun et al., 2014;Wang et al., 2017); the modeling results of assimilation of lower and upper crust material by DM and basaltic melts are reported as black lines; the fractions (in %) of assimilated crustal material are reported as tick marks; (D) εHf(t) vs. εNd(t) diagram (diagram modified after Dobosi et al., 2003 andZhu, 2007). Other data sources are as for Figure 7.
Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 729905 trace element abundances such as Th, U, Nb, and Ta suggest the involvement of an enriched mantle type I (EMI) (Figures 10C,D) and some interactions with crustal material ( Figure 10A). The CVMF Pleistocene basaltic rocks consist of early sodium-rich rocks to later potassic rocks, and the data obtained in this study indeed indicate an alkaline potassic affinity of the rocks emplaced during Calabrian. According to Wang et al. (2017) and Lei et al. (2020), the EMI signature of the NE China potassic volcanism could be due to fluids released by the Pacific sediments of the subducted Pacific slab. More in general, potassic rocks as those studied here are generally considered to be derived from the partial melting of a mantle modified by slab-derived fluids or melts from the subcontinental lithospheric mantle (SCLM) (Foley and Peccerillo, 1992). The last stage of the Pacific subduction plate between NE China and Japan may favor the dehydration of the oceanic crust allowing the partial melting of the mantle. In this framework, the interaction between hydrous fluid and the mantle can produce a phlogopite-bearing source whose partial melting may produce potassic magmas (Wyllie and Sekine, 1982;Condamine and Médard, 2014). Generally, melts in equilibrium with phlogopite are expected to have higher Rb/Sr (>0.1) and lower Ba/Rb (<20) values than those formed from amphibolebearing mantle sources (Furman and Graham, 1999), and the early Pleistocene basaltic rocks of the study area have low Rb/Sr (0.07-0.09) and  values. This suggests that potassic amphibole is the main K-rich phase. Because the mantle material is depleted of Zr and Hf, it is impossible for the mantle derived magma to crystallize zircon without the addition of crustal material. As a result, the occurrence of zircon in the mantle derived magma indicates that some crust-mantle interactions occurred during the magma emplacement (Zheng et al., 2006 Figure 6A). This indicates that there is some crustal material involved in the magma generation processes, i.e., crustal contamination. On the other hand, the consistent negative Nd isotopic compositions (εNd(t) −0.58 to −1.67) of the trachybasalt dikes preclude their derivation from the melting of an old lower continental crust (T DM1 1021-1139 Ma). The variation ranges of the ( 87 Sr/ 87 Sr) I and εNd(t) values in the mafic rocks in the region are 0.70485-0.70619 and +0.02 to −1.65, respectively. These values indicate an initial magmatic source consisting of OIB-EMI type melts ( Figure 8C). Furthermore, the Pb-Pb isotopic data ( Figure 8B) of the trachybasalt dikes and CVMF basaltic rocks support a mixed source consisting of mantle and ancient lower crust components. In the εHf(t)-εNd(t) diagram ( Figure 8D), all the samples fall into the overlapping area of the global lower crust and OIB. The results of the modeling of assimilation of lower crust material in Figure 8C suggest that the studied trachybasalts assimilated about 10-20% of crustal material.

Geodynamic Implications
The late Cenozoic basaltic rocks in the NE China continental margin have geochemical compositions consistent with OIB-type and EMI sources (Zhou and Armstrong, 1982;Xu et al., 2005;Zhi et al., 1990;Liu et al., 2015;Wang et al., 2015;Mao and Wang, 2000;Liu et al., 2008;Chen et al., 2009;Andreeva et al., 2014;Zhang et al., 2009;Ge et al., 2001;Wang et al., 2011). However, the CMVF Calabrian dikes studied here have moderate Mg # (27.85-28.09) and higher Na 2 O and K 2 O contents with respect to the less evolved CMVF basaltic rocks ( Figure 6). Therefore, the dikes do not represent primary magmas and experienced some crystal fractionation. Despite these features, they show geochemical features, e.g., values of Na 2 O/K 2 O and isotopic ratios, and some trace element contents, comparable with FIGURE 9 | (A) A summary of ages of the late Cenozoic volcanic rocks in the Changbai Mountain volcanic area and adjacent areas (e.g., Fan et al., 1999;Nakamura et al., 2000); (B) age frequency diagram of volcanic rocks during the cone construction stage in the Changbaishan Mountains volcanic area (e.g., Liu, 1981Wang et al., 1999Wang et al., , 2001); (C) age frequency diagram of volcanic rocks during the shield-forming stage (other data sources are as for Figure 7). those of the CMVF less evolved rocks ( Figure 6, Figure 7B, and Figures 8B-D). This allows us to discuss the significance of the OIB and EMI signatures of the Pleistocene studied dike in the intraplate geodynamic context of NE China. Seismic tomography studies show that there is a high velocity anomaly in the mantle transition zone (MTZ) below NE China at 410-660 km depth, which is interpreted as the slab of the Pacific plate in the MTZ, and a low velocity body above the flat portion of the slab interpreted as a big mantle wedge (BMW) (Li and Hilst, 2010;Tian et al., 2016;Wang et al., 2020). From the middle Cenozoic, the subducted plate gradually subsided into the MTZ and partly melts. The resulting water-rich melt and the slab-derived fluids may metasomatize the overlying mantle. However, the EMI signature of the studied CMVF dikes cannot be explained by slab-derived sediments and associated metasomatizing fluids (Xu et al., 2012) because of their moderate Sr and low Pb isotopic ratios. In accordance with the geochemical signature of the gas emitted in NE China (Xu S. et al., 2013b) and of the Cenozoic basalts (Zhang et al., 2009), the EM1 source could be the subcontinental lithospheric mantle (SCLM) ( Figure 10D; Zhu et al., 2011). Regarding the source of the late Cenozoic initial OIB magma in the continental margin of NE China, this source could be a response to the back-arc extension and rifting caused by the roll-back of the Pacific plate (Li and Hilst, 2010;Zheng et al., 2018). Some geophysical studies suggest, however, that the Pacific slab is broken and an OIB-type mantle plume upwells from the lower mantle crossing the slab and entering MTZ (Tang et al., 2014). This model is not supported by shear-wave splitting studies (Lynner and Long, 2014) which show trench-parallel sub-slab fast splitting directions and not, as expected for a mantle plume from the lower mantle, a vertical direction. Based on geochemical data, Kuritani et al. (2009) propose that a back-arc type rifting associated to the slab retreat and  (Li, 1992;Deng, 2004;Zhao, 2016), showing that the eastern Jilin Province basaltic rocks (early Pleistocene, including the Changbai Mountain trachybasaltic dikes) mainly plot in the intraplate field. Data of the trachybasalt dikes are from this study, and other data sources are as for Figure 7.
Frontiers in Earth Science | www.frontiersin.org October 2021 | Volume 9 | Article 729905 progressive flattening may favor the decompression melting of the asthenosphere. This model could explain the OIB signature of the CMVF basalts including the trachybasaltic dikes studied here. In this scenario, the CMVF OIB melts upwell from the asthenosphere to the lithosphere, where they mix to SCLM melts. These melts accumulate at the base of the lower crust in the early Pleistocene, where they may assimilate lower crust material of the ancient Asia continent ( Figure 8A and Figure 10A). However, assimilation processes could affect the mantle melts during their upwelling and storage in the crust. Tomographic data in NE China  clearly show a thickening of the Moho due to the accumulation of mantle melts at the base of the crust below CMVF. These data suggest active underplating process where these mantle melts may promote the melting of the lower crust. Therefore, we propose that most of the assimilation of crustal material underwent by the CMVF trachybasalts occurs at the Moho depth and not during their ascent and storage in the crust.

CONCLUSION
Out results indicate that trachybasaltic dikes emplaced in the Calabrian stage of Pleistocene at 1.19-1.20 Ma, i.e., in a time period marked by the transition from a monogenetic-type volcanism to the Tianchi cone construction, central-type activity. The Calabrian volcanism is associated to the emplacement of scoria cones whose shallower feeding systems are controlled by regional tectonics and, in particular, by NW-SE striking structures. Geochemical and Sr-Nd-Pb-Hf isotopic data show that the trachybasaltic dikes have a potassic affinity and their parental magma is constituted by an OIB-type asthenospheric source with an EM I signature related to SCLM melts associated to back-arc and rifting processes linked the slab retreat ( Figure 11). The magma upraising from the mantle stops at the Moho and assimilates Mesoproterozoic lower crust material by underplating ( Figure 11).

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

AUTHOR CONTRIBUTIONS
ML: formal analysis, and writing-original draft. ZX: investigation and data curation. GV: writing-review and editing. XP: project administration and formal analysis. DH: supervision and writingreview and editing. GG: writing-review and editing. DY: investigation and formal analysis. BP: data and formal analysis. JF: drafting the work and revising manuscript critically.