Article Content

1. Introduction

Situated between the Tarim Block, Siberian Block (SB), and North China Craton (NCC), the Central Asian Orogenic Belt (CAOB) (Figure 1a) is the world’s largest and youngest segment of Phanerozoic crust, formed through complex accretionary processes [1,2]. Its evolution is dominated by the bidirectional subduction of the Paleo-Asian Ocean (PAO): northward beneath the SB during the Ordovician-Silurian, and southward beneath the NCC from the Devonian to Triassic [3,4]. The CAOB represents the most significant region documenting Phanerozoic crustal growth, with subduction and accretion processes initiating as early as the Ordovician along the southern Siberian margin. While these early stages were dominated by progressive accretionary tectonics, our study focuses specifically on the terminal closure of the PAO during the late Permian to Triassic along the northern NCC margin. The CAOB represents the most significant region documenting Phanerozoic crustal growth [1,2], with subduction and accretion processes initiating as early as the Ordovician along the southern Siberian margin. While these early stages were dominated by progressive accretionary tectonics, our study focuses specifically on the terminal closure of the PAO during the late Permian to Triassic along the northern NCC margin. Since the Neoproterozoic, the CAOB has undergone complex subduction–accretion processes, closely linked to the double-sided subduction of the PAO plate during the Paleozoic–Mesozoic. Its unique geographical position and geological history make it a focus of geological research [5,6]. The northeastern segment of the CAOB, located within Northeastern China, is traditionally referred to as the Xingmeng Orogenic Belt (XOR). This belt comprises a series of crustal fragments (microcontinental blocks, island arcs, accretionary wedges, and ophiolites that underwent progressive amalgamation [7,8,9]. Recent studies [4] suggest many of these fragments may have shared a common tectonic history since ~500 Ma, though their final amalgamation records the PAO’s subduction and subsequent Siberian-NCC collision [3,10,11]. During the Mesozoic, the XOR was further modified by the superposition of the circum-Pacific and Mongolia–Okhotsk tectonic domains [10], which resulted in an extremely complex tectonic evolution [3,12,13,14,15,16]. Therefore, a better understanding of the geological evolution of the XOR is essential to reconstruct the tectonic history of the CAOB. Current data indicate several issues and controversies regarding the Late Paleozoic to Early Mesozoic evolution of the XOR, particularly concerning the timing and location of the closure of the PAO and the collision processes along the northern margin of the NCC. Two competing viewpoints exist regarding the timing of the closure of the PAO: one hypothesis suggests that the collision between the NCC and northern blocks occurred between the Middle to Late Devonian and Early Carboniferous [10,17], while the other suggests that the collision occurred during the Late Permian–Early Middle Triassic [3,18]. One view places the collision between the NCC and northern blocks along the Solonker–Hegen–Heihe zone [7,19,20], whereas another locates it along the Solonker–Xar Moron–Changchun–Yanji zone [9,13,15,18,21]. Another major debate centers on whether the eastern margin of the NCC was in a collisional or post-collisional stage during the Triassic. Different interpretations include the following: (1) an extensional environment due to post-collisional gravitational collapse triggered by isostatic rebound following slab breakoff [22]; (2) extension-induced decompression melting within the mantle wedge above the subducted slab [19]; (3) possible continuation of the NCC–northern blocks collision into the Middle Triassic [23]. These interpretations reflect two fundamental debates about Triassic tectonics: (1) The tectonic state—whether the region was in post-collisional extension or still experiencing terminal collision; (2) the driving mechanisms of extension, where some models emphasize (a) slab breakoff-induced collapse [22]; (b) mantle wedge decompression [19]; (c) detachment of PAO fragments [14,16,24,25].
Figure 1. Palinspastic reconstruction of the current ocean–continent framework of China and adjacent regions ((a) modified after Shi et al. [26]; sketch tectonic map of the southern Central Asian orogenic belt and the northern margin of the North China block (b), modified after Xiao et al. [3]); geological sketch map (c) and regional location map ((d) modified after Wu et al. [27]) of the study area.
The northern Liaoning region is located south of the Solonker–Xar Moron–Changchun–Yanji suture zone and spans two major tectonic units: the NCC and XOR (Figure 1b). During the Paleozoic, this region recorded the responses to the retreating subduction zone of the PAO, characterized by back-arc basin extension and fragmentation of the continental margin, culminating in the final closure of the PAO. This complex evolution involved intense magmatic activity and diverse geological processes [26,28,29,30,31,32,33,34,35,36,37]. To address stratigraphic uncertainties in northern Liaoning, our geological fieldwork employed detailed structural–lithological mapping, with particular attention to tectonic mélange zones as key markers of paleo-subduction and collision boundaries [38]. This approach follows established methodologies for resolving complex orogenic belts, where mélange zones—despite their complexity—provide critical evidence for reconstructing tectonic interfaces between distinct crustal blocks [38]. This approach enabled us to redefine and formally separate the previously ambiguous “Fangjiatun Formation” and “Tongjiatun Formation”. These two formations were distinguished by contrasting matrix compositions and rock fragment assemblages. Both formations contain a series of Permian–Triassic intermediate metavolcanic rocks. The Fangjiatun–Tongjiatun contact zone, a suture-parallel melange, preserves relics of the PAO subduction. Volcanic clasts in Fangjiatun directly correlate with our dated magmatic suite, while Tongjiatun’s Proterozoic clasts indicate crustal recycling during the NCC–Songliao collision [18,38]. This study investigates the geochronology and geochemistry of these magmatic rocks and explores the nature of the magma source and tectonic evolution in the area covered by the Songliao Basin along the eastern margin of the NCC, providing essential data for understanding the tectonic history of the southern margin of the XOR.

2. Geological Background and Sample Description

The study area is tectonically complex and bordered by the Solonker–Xar Moron– Changchun–Yanji suture zone, with the eastern segment of the CAOB to the north and the NCC to the south (Figure 1c,d). The research area lies within the PAO domain, encompassing distinct near-EW-trending active continental margin belts of both Early and Late Paleozoic age, primarily developed along the northern margin of the NCC since the Ordovician. Since the Mesozoic, the study area has been within the circum-Pacific tectonic domain, with extensive superposition of rift basins and the development of NE- and NNE-oriented structures. The study area preserves a pre-Mesozoic basement composed primarily of highly metamorphosed and deformed intrusive-volcanic complexes, now partially overlain by Mesozoic sedimentary cover. After the Mesozoic, the area became covered by the sedimentary strata of the Songliao Basin and intermediate basic to acidic volcanic rocks. The exposed strata mainly include Early Cretaceous Yixian Formation volcanic rocks, Early Cretaceous Shahuizhen Formation sedimentary strata, and mid-Cretaceous Quantou Formation sedimentary strata. The study area predominantly comprises Silurian, Permian, Triassic, and Jurassic granitic intrusions (Figure 1d).
For this study, Mesozoic–Triassic magmatic rocks were chosen for geochronological and geochemical analyses. Representative sample collection locations are shown in Figure 1c. The samples include the Gongzhuling rhyolite and dacite, Wangjiadian dacite, Wafangxi rhyolite, Haoguantun rhyolite, and Sheshangou pluton.
The Gongzhuling rhyolite (PM107-9, 123°16′32″ E, 42°30′55″ N) and Gongzhuling dacite (PM107-8, 123°16′32″ E, 42°30′54″ N) are located in the eastern part of the study area, covering a small area, in contact with Early to Middle Permian volcanic rock structures. The Gongzhuling rhyolite has a foliated granoblastic texture with foliate and massive structures. The rock consists of quartz (±60%), muscovite (±35%), and hornblende (±5%), with quartz being the predominant felsic mineral (Figure 2a and Figure 3a). The Gongzhuling dacite has a similar foliated granoblastic texture with a mylonitic structure. It is composed of felsic minerals (±55%), biotite (±20%), muscovite (±10%), and hornblende (±15%). The felsic minerals are mainly quartz, which is elongated and oriented, with some parts showing banded quartz (Figure 2b and Figure 3b).
Figure 2. Macro characteristics of Permian–Triassic magmatic rocks in the study area.
Figure 3. Micro characteristics under CPL of Permian–Triassic magmatic rocks in the study area Qtz—quartz; Pl—plagioclase; Fsp—feldspar; Mus—muscovit.
The Wangjiadian dacite (LWY, 123°15′30″ E, 42°33′12″ N) is exposed in the northeastern part of the study area, covering approximately 0.51 km2. The Gongzhuling dacite is tectonically juxtaposed with Early–Middle Permian volcanic rocks within the NNE-trending shear zone, as evidenced by its foliated granoblastic texture and mylonitic structure. This contact represents a tectonic boundary rather than a magmatic inclusion, confirming its emplacement during post-collisional deformation. The rock is composed of felsic minerals (±57%), muscovite (±25%), biotite (±15%), and hornblende (3%). The felsic minerals include feldspar and quartz, with some feldspar exhibiting a ball structure, indicating a mylonitic structure. Biotite and muscovite are arranged in a continuously oriented manner, forming a foliated structure (Figure 2c and Figure 3c).
The Wafangxi (PM201-5, 123°33′16″ E, 42°39′12″ N) and Haoguantun rhyolites (PM202-6, 123°34′21″ E, 42°37′57″ N) are exposed in the northwestern part of the study area and are in contact with Early to Middle Permian volcanic rocks. The Wafangxi rhyolite covers an area of approximately 1.1 km2. The rock appears grayish-white, with a mylonitic structure, and the matrix has a foliated granoblastic texture containing phenocrysts of plagioclase and quartz. Plagioclase (±7%) develops polysynthetic twinning, with weak clay and carbonate alteration. Quartz (±5%) forms several lens-shaped and short-banded aggregates that are oriented. The matrix is mainly composed of small felsic minerals and muscovite, with felsic minerals (±48%) being granular and embedded with muscovite (±40%), showing an oriented distribution (Figure 2d and Figure 3d).
The Haoguantun rhyolite (PM202-6, 123°34′21″ E, 42°37′57″ N) is exposed in the northwestern part of the study area, covering an area of approximately 1.2 km2. The rock appears grayish-white, with a weak mylonitic structure, predominantly composed of felsic components. The rock consists of muscovite (±5%), feldspar (±15%), alkaline feldspar (±45%), and quartz (±35%). The feldspar and quartz appear elongated or lens-shaped, with muscovite appearing as fibrous sheets, constituting a weak mylonitic structure. Plagioclase exhibits polysynthetic twinning. Alkaline feldspar is composed of microcline and exhibits grid-like twinning (Figure 2e and Figure 3e).
The Sheshangou pluton (SSGC, 123°29′47″ E, 42°30′34″ N) is exposed in the central part of the study area, covering 12.95 km2. It intrudes Middle Triassic granite and Middle Permian quartz diorite bodies. The lithologies include monzogranite and muscovite-bearing monzogranite. The rock appears gray with a faint pink hue, exhibiting a medium-to-coarse-grained granite structure and massive texture. The mineral composition includes muscovite (±5%), plagioclase (±30%), alkaline feldspar (±40%), and quartz (±30%). Plagioclase exhibits polysynthetic twinning, while alkaline feldspar appears as striped feldspar, with some developing Carlsbad twinning (Figure 2f and Figure 3f).
Major elements were analyzed by XRF (1–5% precision) and trace elements by ICP-MS (>10% precision), using contamination-free equipment. Zircon U-Pb dating employed a 193 nm laser ablation system (32 μm spot) coupled with ICP-MS, with data processed using ICPMSDATACAL and Isoplot, following Andersen’s Pb-correction method [39] and Yuan et al.’s protocol [40]. Zircon Lu-Hf isotopes were analyzed via laser ablation MC-ICP-MS (60 μm spot), yielding 176Hf/177Hf values of 0.282316 ± 30 (91500) and 0.282507 ± 50 (MT), with methods detailed in Geng et al. [41].

3. Test Results

3.1. Chronological Characteristics

Zircon age determinations were conducted on the six collected samples. Colorless, transparent zircons without inclusions or fissures were selected for U-Pb isotopic analysis. These zircons exhibit relatively high Th/U values (0.37–1.53) (Table 1) and display euhedral–subhedral shapes, with sizes ranging from 50 to 200 μm and aspect ratios of 1:3. They exhibit visible oscillatory zoning (Figure 4), characteristic of magmatic zircon [42], representing the crystallization age of the volcanic rock.
Figure 4. Cathodoluminescence (CL) images of representative zircons of the Permian–Triassic magmatic rocks in the study areas.
Table 1. Zircon LA-ICP-MS U-Pb dating results for the Permian–Triassic magmatic rocks of the study area.
For the Gongzhuling rhyolite (PM107-9), the Th/U values range from 0.53 to 1.08. They exhibit typical magmatic zircon characteristics, with some zircons showing narrow white bright rims, indicating weak recrystallization. The zircons are elongated, ranging in length from 80 to 120 μm, with aspect ratios of 1:2.5 to 1:4 (Figure 4a). A total of 17 valid test points were analyzed, and 206Pb/238U surface ages ranging from 259.5 to 262.2 Ma (Table 1) were obtained. The weighted average 206Pb/238U age is 260.5 ± 2.2 Ma (MSWD = 0.031, n = 17) (Figure 5a).
Figure 5. Zircon U-Pb concordia diagram from the Permian–Triassic magmatic rocks of the study area.
Zircons from the Gongzhuling dacite (PM107-8) exhibit distinct oscillatory zoning (Figure 4b), with Th and U contents ranging from 61.9 × 10−6 to 404.6 × 10−6 and 86.2 × 10−6 to 665.8 × 10−6, respectively, with Th/U values ranging from 0.37 to 0.85 (Table 1). The 16 measurement points exhibit 206Pb/238U ages falling on or near the concordia curve, ranging from 254.5 to 269.1 Ma (Table 1), with a weighted average age of 260.3 ± 2.4 Ma (MSWD = 1.70, n = 16) (Figure 5b).
For the Wangjiadian dacite (LWY), the zircons are semi-automorphic short columnar, with aspect ratios of 1:1 to 1:3 (Figure 4c). Th/U values range from 0.39 to 1.50. The 206Pb/238U surface ages of the measurement points range from 236.2 to 252.8 Ma and fall on or near the concordia curve (Table 1). The weighted average age is 244.1 ± 2.5 Ma (MSWD = 0.97, n =15) (Figure 5c).
Zircons from the Wafangxi rhyolite (PM201-5) range in length from 80 to 160 μm, with aspect ratios of 1:1 to 1:3 (Figure 4d) and Th/U values between 0.42 and 1.16. The surface ages of the measurement points range from 239.3 to 250.7 Ma (Table 1), falling on or near the concordia curve, yielding a weighted average age of 243.9 ± 3.0 Ma (MSWD = 0.68, n = 11) (Figure 5d).
Zircons from the Haoguantun rhyolite (PM202-6) are euhedral or subhedral, with well-defined internal structures and growth zoning, and some display core–rim structures. The aspect ratios range from 1:1 to 1:2.5 (Figure 4e), and Th/U values range from 0.74 to 1.53. The 11 valid measurement points yielded surface ages ranging from 239.3 to 246.1 Ma (Table 1), falling on or near the U-Pb age concordia, with a weighted average age of 240.9 ± 2.2 Ma (MSWD = 0.15, n = 11) (Figure 5e).
Zircons from the Sheshangou pluton (SSGC) are semi-automorphic, long columnar, with long-axis lengths of 50–100 μm and aspect ratios of 1:1 to 1:3 (Figure 4f). Th/U values range from 0.70 to 1.11. A total of 12 valid measurement points were analyzed, and surface ages ranging from 228.3 to 231.8 Ma (Table 1) were obtained. The ages fall on or near the concordia curve, with a weighted average age of 230.1 ± 1.7 Ma (MSWD = 0.28, n = 12) (Table 1).
While some samples exhibit intense mylonitization, zircon U-Pb ages reflect magmatic crystallization rather than recrystallization events. This is evidenced by (1) concordant ages with magmatic Th/U (>0.4 in 148/159 spots); (2) juvenile εHf(t) values uncorrelated with deformation intensity; (3) exclusion of inherited cores (e.g., 487 Ma in PM202-6) via CL-guided spot placement. Low Th/U values (0.37–0.44) in PM107-8 reflect late-stage magmatic crystallization, not metamorphic overprinting—consistent with their high-U (>200 ppm) and trace element signatures (Ce/Ce* > 10).
While some samples show intense mylonitization, zircon U-Pb ages reflect magmatic crystallization rather than recrystallization events, supported by the following: (1) Concordant ages with magmatic Th/U (>0.4 in 148/159 spots), typical of magmatic zircon, contrasting with metamorphic zircon (Th/U < 0.1). (2) Juvenile εHf(t) values independent of deformation: Zircon Hf isotopes exhibit juvenile characteristics (e.g., εHf(t) = +2.3 to +5.1 in PM107-8) with no correlation to deformation intensity, confirming magmatic origin. (3) CL-guided exclusion of inherited cores: Cathodoluminescence imaging targeted magmatic oscillatory zoning, excluding inherited detrital cores (e.g., 487 Ma in PM202-6). (4) Low Th/U values indicating late magmatic differentiation: PM107-8 zircons with Th/U = 0.37–0.44 show high U (>200 ppm) and Ce/Ce* > 10, consistent with late-stage magmatic crystallization rather than metamorphic overprinting. (5) Greenschist facies constraints: Metamorphic conditions (≤350 °C, <0.3 GPa) were insufficient to reset zircon isotopes, as evidenced by preserved magmatic zoning and lack of metamorphic inclusions. (6) Spatial–temporal separation from younger intrusions: Mylonitized samples lie >10 km from the 230.1 Ma Sheshangou pluton, outside its <2 km contact metamorphic halo, ruling out thermal disturbance [3,18,24,42].

3.2. Geochemical Characteristics

The results of the major and trace element analyses are presented in Table 2 (all data used for plotting were calculated after loss on ignition).
Table 2. Major (wt%) and trace (×10−6) elements for the Permian–Triassic of magmatic rocks of the study area.
The Gongzhuling rhyolite (PM107-9) exhibits a relatively high SiO2 content (75.00%–76.20%), with smaller amounts of CaO (3.11%–3.62%) and Na2O+K2O (3.33%–3.88%) (Table 2). The TAS diagram shows that the samples primarily fall within the rhyolite range (Figure 6a). The Rittmann Index (σ) ranges from 0.34 to 0.47 (Table 2). The aluminum saturation index (A/CNK) is 1.17–1.20 (Table 2), classifying the rock as peraluminous. The ratio of light rare earth elements (LREEs) to heavy rare earth elements (HREEs) for the Gongzhuling rhyolite samples is 4.53–5.04 (Table 2), with a notably high total rare earth element content (183.12 × 10−6–248.76 × 10−6). The samples exhibit a strong negative Eu anomaly (δEu = 0.26–0.42), and the rare earth element distribution pattern is overall right-leaning (Figure 7a). The element spider diagram shows that the samples are relatively enriched in large-ion lithophile elements (LILEs), such as Rb, Th, and K, and depleted in high-field strength elements (HFSEs), such as Nb, P, and Ti (Figure 8a).
Figure 6. SiO2 vs. total alkali (Na2O+K2O) ((a) modified after Middlemost [43]), SiO2 vs. K2O ((b) modified after Peccerillo and Taylor [44]), and A/CNK vs. A/NK ((c) modified after Maniar and Piccoli [45]) diagrams for the Permian–Triassic magmatic rocks of the study area.
Figure 7. Chondrite-normalized REE patterns (normalization values after Taylor and McLennan [46] for the Permian–Triassic magmatic rocks of the study area.
Figure 8. Chondrite-normalized trace element spider diagrams (normalization values after Sun and McDonough [47]) for the Permian–Triassic magmatic rocks of the study area.
The TAS diagram for the Gongzhuling dacite (PM107-8) indicates that samples primarily fall within the dacite range (Figure 6a). Although both rhyolite and dacite occur within a NE-trending shear zone, their distinct geochemical signatures (e.g., SiO2 content, A/CNK ratios, and trace element patterns) are interpreted as primary magmatic features rather than alteration products. These samples have a relatively lower SiO2 content (65.33%–66.19%), MgO content (1.44%–1.75%), Mg# values of 29.6–34.7, and K2O/Na2O ratios of 0.64–0.91 (Table 2). The SiO2-K2O diagram shows that the samples fall within the calc-alkaline series (Figure 6b). The Al2O3 content of the samples ranges from 13.90% to 14.39%, indicating that they are quasi-aluminous rocks (Figure 6c). The total rare earth element content is relatively high (105.25 × 10−6–141.00 × 10−6), with LREE/HREE ratios of 4.53–5.04. The Gongzhuling dacite is enriched in LREEs and relatively depleted in HREEs (Table 2) with a right-leaning rare earth element distribution (Figure 7b). The element spider diagram shows that the samples are relatively enriched in LILEs (Rb, Ba, Th, and K) and depleted in HFSE (Nb, P, and Ti) (Figure 8b).
The Wangjiadian dacite (LWY) samples exhibit a relatively low SiO2 content (63.01%–64.24%) and MgO content (2.15%–2.88%), with Mg# values of 27.0–39.7. These samples have a higher Al2O3 content (15.02%–15.70%) (Table 2). The TAS diagram shows that the samples mainly fall within the andesitic range (Figure 6a). The Rittmann Index (σ) varies from 1.38 to 1.71 (Table 2), indicating that they belong to the calc-alkaline series (Figure 6b). The aluminum saturation index (A/CNK) values range from 1.33 to 1.55 (Table 2), classifying these rocks as peraluminous. These samples have a relatively high total rare earth element content (114.72 × 10−6–141.82 × 10−6), and the rare earth element distribution pattern is overall right-leaning (Figure 7c). The Wangjiadian dacite is enriched in LREEs and relatively depleted in HREEs ((La/Yb)N ranges from 3.12 to 5.15) (Table 2), with a slight negative Eu anomaly (δEu = 0.88–0.96). The element spider diagram shows that the samples are relatively enriched in LILEs (Rb, Th, Ba, and K) but depleted in HFSEs (Nb, P, and Ti) and Sr (Figure 8c).
The Wafangxi rhyolite (PM201-5) has a relatively low SiO2 content (66.99%–68.30%) and MgO content (0.77%–1.20%), with Mg# values of 33.7–42.8 and K2O/Na2O ratios of 0.88–1.19 (Table 2). The TAS diagram shows that the samples mainly fall within the rhyolite range (Figure 6a). The SiO2-K2O diagram further classifies them within the high potassium and calc-alkaline series (Figure 6b). The Al2O3 content of the samples is in the range of 14.03%–15.47% (Table 2), indicating that they are quasi-aluminous rocks (Figure 6c). The rare-earth element distribution pattern is overall right-leaning (Figure 7d). The total rare earth element content is high (ΣREE = 164.14 × 10−6–238.74 × 10−6), with LREE/HREE ratios ranging from 7.07 to 8.19. The Wafangxi rhyolite is enriched in LREEs and relatively depleted in HREEs ((La/Yb)N ranges from 5.68 to 7.58) (Table 2). The element spider diagram shows that the samples are relatively enriched in LILEs (Rb, Ba, Th, and K) and depleted in HFSEs (Nb, P, and Ti) (Figure 8d).
The TAS diagram for the Haoguantun rhyolite (PM202-6) shows that the samples fall within the rhyolite range (Figure 6a). The samples exhibit high SiO2 content (77.77%–78.33%), Na2O+K2O content (8.00%–8.11%), K2O/Na2O ratios of 1.02–1.34 (Table 2), and low FeOT content (0.85%–1.08%). The Rittmann Index (σ) is 1.81–1.89 (Table 2), and the SiO2-K2O diagram shows that the samples primarily fall within the high potassium and calc-alkaline series (Figure 6b). The samples are enriched in LREEs and relatively depleted in HREEs ((La/Yb)N values in the range of 2.84–4.56) (Table 2), with negative Eu anomalies (δEu = 0.31–0.47). The rare earth element distribution pattern is right-leaning (Figure 7e). The element spider diagram shows that the samples are relatively enriched in LILEs (Rb, Th, and K) and depleted in HFSEs (Nb, P, and Ti) (Figure 8e).
The Sheshangou pluton exhibits a high SiO2 content (75.62%–77.49%) and low MgO content (0.03%–0.41%), with Mg# values ranging from 35.0 to 46.9 (Table 2). The TAS diagram classifies the samples within the granite range (Figure 6a). The Rittmann Index (σ) ranges from 1.71 to 2.02 (Table 2). The SiO2-K2O diagram indicates that the samples predominantly fall within the high potassium and calc-alkaline series (Figure 6b). The Sheshangou pluton samples exhibit strong negative Eu anomalies (δEu = 0.04–0.06) (Figure 7f) and a relatively high total rare earth element content (ΣREE = 48.54 × 10−6–66.67 × 10−6), with LREE/HREE ratios of 2.80–3.08. The Sheshan Gou pluton is enriched in LREEs and relatively depleted in HREEs. The trace element spider diagram shows that the samples are relatively enriched in LILEs (Rb, Th, and K) and depleted in HFSEs (Nb, P, and Ti) and Sr (Figure 8f).

3.3. Hf Isotope Characteristics

This study conducted in situ Lu-Hf isotope analyses on the magmatic rocks with determined zircon ages. All analyses yielded fLu/Hf values between −0.98 and −0.90 (Table 3), which are significantly lower than those of mafic (−0.34, Amelin et al. [48]) and silicon-aluminous crust (−0.72, Vervoort et al. [49]), suggesting that the two-stage model age better reflects the time when the source material was extracted from the depleted mantle or the average residence age of the source material in the crust [50].
WhatsApp