Next Article in Journal
Properties of Sulfur Particles Formed in Biodesulfurization of Biogas
Next Article in Special Issue
Zircon U-Pb Dating and Petrogenesis of Multiple Episodes of Anatexis in the North Dabie Complex Zone, Central China
Previous Article in Journal
A Comparative Study of Porphyry-Type Copper Deposit Mineralogies by Portable X-ray Fluorescence and Optical Petrography
Previous Article in Special Issue
Genesis of Early–Middle Jurassic Intrusive Rocks in the Erguna Block (NE China) in Response to the Late-Stage Southward Subduction of the Mongol–Okhotsk Oceanic Plate: Constraints from Geochemistry and Zircon U–Pb Geochronology and Lu–Hf Isotopes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Petrogenesis and Tectonic Implications of the Early Cretaceous Granitic Pluton in the Sulu Orogenic Belt: The Caochang Granitic Pluton as an Example

1
College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, China
2
Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China
3
Key Laboratory of Deep-Earth Dynamics of Ministry of Natural Resources, Beijing 100037, China
4
Qingdao Institute of Marine Geology, China Geological Survey, Qingdao 266071, China
5
Department of Geology, Northwest University, Xi’an 710069, China
*
Author to whom correspondence should be addressed.
Minerals 2020, 10(5), 432; https://doi.org/10.3390/min10050432
Submission received: 22 April 2020 / Revised: 7 May 2020 / Accepted: 8 May 2020 / Published: 11 May 2020

Abstract

:
The Sulu orogenic belt is the source of information on important magmatic events associated with the collision of the Yangtze craton and North China craton (NCC) and the destruction of the NCC during the Mesozoic in eastern China. In this study, we have, for the first time, identified a monzonitic granitic pluton. We hereby present petrological, geochemical, and zircon U-Pb-Hf-O isotopic data, shedding new light on the petrogenesis and tectonic implications for the granitic pluton in the Sulu belt. LA-ICP-MS and SHRIMP II analyses of zircon grains suggest that the monzonitic granitic pluton was crystallized in the Early Cretaceous (ca. 120 Ma). Geochemically, the granitic pluton shows sub-alkaline, high-K calc-alkaline, and metaluminous signatures, and is genetically of I-type granite, excluding the possibility of S-type granite, as evidenced by mantle-like zircon oxygen isotopic features. In addition, the pluton is enriched in light REE and large-ion lithophile elements (LILE) (e.g., La, Cs, Ba, K, and Pb), but depleted in high-field-strength elements (HFSE) (e.g., Nb, Ta, P, and Ti), suggesting an arc-related affinity. Zircon Hf isotopes (εHf(t) = −27.51~−32.35; TDM2 = 2979~3175 Ma) and mantle-like δ18O values (5.12–6.24‰) together indicate that the identified granitic pluton is derived from the partial melting (reworking) of the ancient mafic lower crustal material, with no supra-crustal material participation. Moreover, high Magnesium number (Mg# = 42–49) values and mafic micro-granular enclaves suggest that mantle-derived magma participated in the evolution of the granitic pluton in this study. Integrating the findings of this study and previous work, we propose that the Caochang granitic pluton is derived from the partial melting of the deep Yangtze basaltic lower crust during the Early Cretaceous, and that the large-scale thinning of the lithospheric mantle was the main factor that led to Early Cretaceous magmatic flare-up in the Sulu orogenic belt.

1. Introduction

The Sulu orogenic belt is an ideal place for studying continental deep subduction and the crust–mantle interaction [1,2,3,4,5,6,7]. Multiple studies have reported that the Sulu orogenic belt experienced two stages of evolution, namely the deep-subduction stage and the subsequent exhumation stage [4,8,9]. In addition, multiple geodynamic scenarios have been constructed in the Sulu belt [1,2,10,11,12,13]. Early Cretaceous granitoid rocks (130–110 Ma), which are widely distributed throughout the Sulu orogenic belt, have been regarded as post-collisional products [4,5,11,14,15]. Competing models were proposed to explain the petrogenesis of the Early Cretaceous granitoid rocks in the Sulu belt [16,17,18,19,20]. Niu [19] and Gao [20] proposed that the granites were most probably derived from anataxis of the continental lower crust (ca. 20–25% partial melting of the mafic granulite), and further pointed out that the stagnant paleo-Pacific slab beneath the subcontinental lithospheric mantle was the main inducing factor that brought about convective asthenosphere, leading to lower crustal partial melting. However, Dai [21] and Wang [12] argued that rollback of the western paleo-Pacific slab caused the strongest mantle perturbation, resulting in large-scale partial melting and strong Early Cretaceous (130–110 Ma) magmatic events in the Dabie–Sulu orogenic belt and other regions of eastern China. In addition to above models, Ling [22] argued that Early Cretaceous magmatism was related to ridge subduction, and stressed that the ridge subduction was key driving force leading to the destruction of the North China craton (NCC) in ca. 125 Ma. Additionally, the magma sources of Early Cretaceous granites are in dispute. Yan [23] and Zhao [5] concluded that the granites were produced by partial melting of upper continental crust, while Gao [20] argued that granites were products of lower continental crust with injection of mantle magma, and reported that upper crustal material played a minor or negligible role. Together with zircon oxygen isotopes (δ18O), Meng [11] considered the felsic rocks to be the products of melting of deep crustal material which did not experience high-temperature water–rock interactions. With regard to the magma source, it was proved that granites from the Sulu belt were products of partial melting of the Yangtze crust [4,18,24,25]. According to a recent new study, the under-thrust lower crust of the Yangtze was the source for the granites in both the Jiaobei terrane (akin to the NCC affinity) and the Sulu belt [26]. This study challenged traditional notions, and implied that crustal material of the NCC was not the sole source of Late Jurassic to Early Cretaceous granites in the Jiaobei terrane. In order to evaluate related debates and issues, we identified a new granitic outcrop (artificial exposure—a quarry) in order carry out new petrological, geochronological, and geochemical studies for the first time in the central segment of the Sulu belt. Based on these datasets, we discuss its petrogenesis and the geodynamic scenario, and evaluate the tectonic setting of granites in the study area.

2. Regional Geological Setting and Petrological Features

Owing to the left-lateral movement of the Tanlu fault in the Early Jurassic, the Dabie–Sulu orogenic belt was separated into two parts, the Sulu belt to the east and Dabie belt to the west (Figure 1a). According to regional faults and mélange, the Sulu belt is bound by the Jiashan–Xiangshui fault (JXF) to the south and by the Baichihe (Wulian)–Qingdao-Yantai fault (WQYF) to the north [27,28]. Furthermore, Liu [3] proposed that the whole Sulu orogenic belt could be divided into the ultra-high pressure (UHP) and high-pressure (HP) sub-zones that are characterized in the Shuyang–Jinping suture zone. The UHP and HP rocks experienced widely retrograde metamorphism (amphibolite facies), and subsequently were intruded upon by Late Jurassic–Early Cretaceous granites [3] (Figure 1b). It is accepted that the framework of the Dabie–Sulu orogenic belt was well constructed. According to metamorphism and magmatism in the Dabie–Sulu belt, four important stages were recognized: (1) Continental deep subduction in the Triassic between the Yangtze and north China craton (ultra-high pressure magmatism); (2) Syn-exhumation anataxis in the Late Triassic (syn-exhumation magmatism); (3) Anataxis of the thickened crust in the Late Jurassic (post-collisional magmatism); and (4) Anataxis of orogenic crust in the Early Cretaceous (post-collisional magmatism) [1,2,4,8,10,14]. The four stages correspond with three different magmatic activities in the Sulu belt [1,2,4,5,14,24,25,29], with the stages containing: (1) 225–205 Ma; (2) 160–150 Ma; and (3) 130–110 Ma. The Late Triassic magmatic rocks, composed of granites and mafic complexes, only outcrop in the Shidao area. Geochemical and isotopic evidence suggests them to be related to syn-exhumation tectonism [14]. Compared to the Late Triassic magmatic rocks, the Late Jurassic magmatic rocks have more outcrops, which are mainly distributed in Rushan and Wendeng areas [24,25]. Multiple studies demonstrated that the Late Jurassic magmatic rocks were generated from partial melting of thickened continental crust [4]. The Early Cretaceous magmatic activities were significantly sharp (magmatic flare-up), and cropped out throughout the Sulu orogenic belt [4,18]. Early Cretaceous magmatic rocks are generally considered to be products of the NCC destruction. The Late Jurassic and Early Cretaceous magmatism had a post-collisional background and was formed in the within-plate setting.
The study region, located in the middle segment of the Sulu belt, is adjacent to the village of Caochang. Thus, we denominated the newly identified pluton as Caochang granite (Figure 1c,d and Figure 2a,b). Although the simplified map shows that the study region is characterized by granites and UHP rocks, the study region is covered by quaternary covers, as suggested by Google Earth and a field survey (Figure 1c and Figure 2a). Therefore, new identified exposures could provide new clues to evaluate the debates, as described in the introduction. Using 1:50,000 and 1:250,000 mapping and a geological survey, the study region was found to be characterized by granitoid rocks and granitic gneisses (Figure 1d).
The newly identified pluton showed good artificial exposures (Figure 2a,b). We collected representative samples from the limited rock exposures for research. The field observations revealed that the granitic pluton showed typical magmatic structures, and few mafic micro-granular enclaves were found in the granitic pluton (Figure 2c). The granitic samples showed medium-grained and equigranular textures, mainly comprising euhedral K-feldspar (30~35%), plagioclase (~35%), biotite (~15%), and quartz (10~20%), as well as a small amounts of accessory minerals (~5%) (Figure 2d–f). Equivalent contents (volume) of plagioclase and K-feldspar suggested that the granitic pluton was composed of monzonitic granite. K-feldspar crystals were marked by earthy appearances and Carlsbad twinning, with few crystals showing envelope and sector-shaped twinning (Figure 2d,e). Plagioclase grains had typical zonal textures and polysynthetic twinning, and few crystals were subjected to sericite (Figure 2d,e). Compared to the host granite, the mafic micro-granular enclaves (MMEs) were crystallite textures marked by sub-euhedral to anhedral biotite, amphibole, and plagioclase crystals (Figure 2f–i). The MMEs had a clear boundary with the host pluton, and few sanidine crystals were identified near the contact boundary (Figure 2f). The sanidine crystals, which were marked by clear appearances, suggested a quenching process. Additionally, acicular apatite crystals also suggested a quick cooling process (Figure 2h). As such, feldspar and quartz megacrysts indicated a magma mixing process between the MME and host granite (Figure 2g,i).
Field observations and petrological-mineral analyses suggested that the granitic pluton was composed of monzonitic granite, corresponding to CIPW calculation results.

3. Methods

3.1. Zircon U-Pb Geochronology and Oxygen Isotope Analyses

Zircon separation and cathode luminescence (CL) images were undertaken at the Institute of Geology, Chinese Academy of Geological Sciences (Beijing, China). Zircon LA-ICP-MS U-Pb dating was completed at the Wuhan Sample Solution Analytical Technology Co., Ltd (Wuhan, China) using a GeolasPro laser ablation system comprising a COMPexPro 102 ArF excimer laser and a MicroLas optical system. In addition, ICP-MS (Agilent 7700e) was adopted to obtain ion signal intensities. During the experiment, helium (He) was applied as a carrier gas. Standard zircon 91500 and glass NIST610 were used to correct the U-Th-Pb isotope system and trace elements, respectively. Due to weak signals of 204Pb and interference of 204Hg, the traditional method could not be used to determine accurate contents of U and Pb. Therefore, we used Com Pb Corr#_315G software embedded in EXCEL to carry out common Pb corrections. The detailed analytical procedures and principles were described by Liu [30] and Zong [31]. Zircon SHRIMP II analyses were performed at the Institute of Geology, Chinese Academy of Geological Sciences. A primary O2− ion bean of 3~6 nA was utilized to bombard the surface of zircons with a 20-µm spot size in diameter. The standard zircons 91500 (U = 91 ppm), SL13 (U = 238 ppm), and M257 (U = 840 ppm) were used to correct the U-Th-Pb isotope system. TEMORA, dated at ca. 417 Ma, was analyzed for calibration of 206Pb/238U ratios per three analyses. The related techniques and processes were described by Williams [32]. Zircon oxygen isotopes were also completed at the Institute of Geology, Chinese Academy of Geological Sciences (Beijing, China), using the SHRIMP II MC technique (Beijing SHRIMP Center, Beijing, China). Detailed oxygen isotopic procedures and analytical principles were described by Black [33].

3.2. Zircon Lu–Hf Isotopes

In situ Zircon Lu–Hf isotopes were obtained at the Wuhan Sample Solution Analytical Technology Co., Ltd. using a Neptune Plus MC-ICP-MS (Thermo Fisher Scientific, Germany) in combination with a Geolas HD excimer ArF laser ablation system. Helium (He) was used as the carrier gas within the ablation cell and was merged with argon (makeup gas) after the ablation cell. During the experiment, standard zircons GJ-1 (internal standard) and 91500 (external standard) were used to verify and correct measured data. The assayed Hf values were consistent with the recommended values within the uncertainty. Detailed analytical conditions and procedures, as well as processes, are described by Hu [34].

3.3. Major and Trace Whole-Rock Elements

Major and trace whole-rock element compositions were assayed at the Testing center of Shandong Bureau of China Metallurgical Geology Bureau (Jinan, China). Major elements were determined using the ARL 9900XP XRF technique, and trace elements were assayed by the ICP-MS method. The analytical uncertainties were controlled within 0.1–1% (RSD), and loss on ignition (LOI) was obtained by means of a 500-µg rock powder that was heated up to ~1100 °C for one hour. The precision of the ICP-MS was ±10%, and some trace elements were found to represent ±5%.

4. Results

4.1. Zircon U-Pb Dating Results

Three samples from the different exposures of the granitic pluton were selected for zircon U-Pb dating. Zircon CL images revealed that zircons from the investigated area showed prismatic crystal morphologies, with ratios of length and width ranging from 1:1 to 4:1 (Figure 3). High ratios of Th/U (>0.4) and clear zircon zonal textures suggested an igneous origin [35]. The zircon LA-ICP-MS and SHRIMP II U-Pb dating results are listed in Table 1.
In total, 19 and 26 spots were assayed using LA-ICP-MS for samples My121 and My122, respectively. Excluding a few dating spots, the My121 sample had ages of 206Pb/238U ranging from 113.9 ± 2.1 to 126.2 ± 3.2 Ma, with a weighted average age of 119.6 ± 1.8 Ma (MSWD = 2.7). The sample My122 yielded ages of 206Pb/238U varying from 114.3 ± 1.5 to 126.9 ± 2.4 Ma, with a weighted average age of 118.7 ± 1.3 Ma (MSWD = 2.2). In order to further address absolute emplacement age of the pluton, the third sample, My1211, which contained 14 effective spots, was analyzed using the SHRIMP II method. The analytical results suggested the absolute crystallization age of the pluton was ca. 121.9 ± 1.4 Ma (MSWD = 1.4). Figure 4 demonstrates that all the dating spots fell along the concordant curves, suggesting no radiogenic Pb losses.
Despite uncertainty, two analytical methods were consistent, and both indicated that the Caochang monzonitic granitic pluton was crystallized at the Aptian stage of the Early Cretaceous, corresponding to the volcanic activities of the Qingshan Group in the study area [11].

4.2. Zircon Hf–O Isotopic Results.

A total of 40 in situ zircon Lu–Hf analyses were conducted on samples of My121 and My122 (one spot was excluded). The analytical data of zircon Lu–Hf are shown in Table 2. The analytical results revealed that values of 176Lu/177Hf ranged from 0.000621 to 0.00314 (mean value = 0.001215), with the majority of Lu–Hf spots having a low value (low accumulation) (<0.002), corresponding to a high quality of data (Figure 5a). The εHf(t) values of assayed zircons ranged from −27.51 to −32.35, with a mean value of −30.21 (Figure 5b,c), and two-stage model ages (TDM2) varied from 2979 to 3175 Ma, with a mean age of 3053 Ma (Figure 5d).
Thirty representative zircon grains were selected for in situ oxygen analyses. The oxygen isotopic results are listed in Table 3. Analytical results revealed that zircon grains had a narrow range of δ18O values (5.12–6.24‰) that were consistent with the normal mantle range (Figure 5e), reflecting partial melting of normal mantle material or rocks. The decoupling of Hf–O isotopes suggests that the granitic pluton originated from an ancient mafic lower crust-dominated source (Figure 5f).

4.3. Whole-Rock Geochemical Results

A total of nine representative samples were selected for whole-rock geochemical analyses. The analytical results are listed in Table 4. In the following diagrams, the major elements were recalculated to 100% on a volatile-free basis.
In the TAS diagram, the samples are plotted in the granodiorite to quartz monzonite fields (Figure 6a). According to the CIPW results and mineral features, we identified the identified pluton as monzonitic granite. The studied samples showed calc-alkaline features, as suggested by the A.R. versus SiO2 diagram (Figure 6b). A further classification diagram suggested that the monzonitic granite belonged to high-K clac-alkaline series, but one sample fell in the medium-K calc-alkaline field that might be related to magma mixing (with greater MME composition participation) (Figure 6c). In the molar Al2O3/CaO + Na2O+ K2O (A/CNK) versus molar Al2O3/Na2O + K2O (A/NK) diagram (Figure 6d), the studied samples demonstrated metaluminous affinity corresponding to I-type granite. Detailed petrogenetic classifications are discussed below.
The total rare earth element (∑REE) contents ranged from 143 to 188 ppm, and light REE (LREE)/heavy REE (HREE) and (La/Yb)N ratios varied from 11.28 to 14.59 and from 14.34 to 21.33, respectively, both suggesting an enrichment of LREE compared to HREE (Figure 7a). In addition, the samples showed weak Eu anomalies, with Eu/Eu* ratios varying from 0.76 to 0.82 (mean ratio = 0.78). In the primitive-normalized spider diagram (Figure 7b), the majority of samples were enriched in large-ion lithophile elements (e.g., Cs, K, Ba, and LREE), but depleted in high-field-strength elements (HSFE), suggesting an arc-type or crustal magma nature [37]. Additionally, the samples were enriched in Pb, Zr, and Hf, also corresponding to crustal affinity. Combined with above geochemical features, we tentatively argue that the Caochang granitic pluton is most probably a result of the partial melting of crustal material; detailed petrogenetic process of the granitic pluton and deep insights are addressed below.

5. Discussion

5.1. Classification for the Granitic Pluton

According to petrogenetic setting and magma source, granitoids can be divided into the M-, I-, A-, and S-types [47]. A-type granites are usually produced in high-temperature anorogenic settings (e.g., the extensional setting), and are marked by high alkaline contents and anhydrous features [48,49]. It is well recognized that eastern China experienced strongly lithospheric thinning (foundering) during the Early Cretaceous (ca. 125 Ma) [50], leading to massive mantle magma underplating and extensional tectonic formation (e.g., the metamorphic core complex). Thus, the underplating of mantle magma with the heat flux and extensional tectonics provided an advantage for the formation of A-type granites in the study region during the Early Cretaceous. Considering this, we will evaluate the possibility that is akin to A-type granite in the study region.
Discrimination diagrams revealed that the studied samples belonged to ordinary granite, excluding the A-type possibility (Figure 8a–c). Additionally, the combination of whole-rock Zr and zircon Ti thermometers revealed a low crystallization temperature [51] (Figure 8d,e) that is inconsistent with typical high-temperature A-type granite [49]. The normalized-REE patterns were in agreement with I-type granitic rocks that were reported in the Sulu orogenic belt [20,52] (Figure 7a). Typical A-type granitic rocks, which were reported in the Sulu belt, showed significant a Eu-negative anomaly corresponding to highly-fractionated granite (Figure 7a) [20,23]. These genetic features suggest that the new identified granitic pluton belongs to I- or S-type granite rather than A-type and high-fractionated granite.
Additionally, in the SiO2 versus P2O5 diagram (Figure 8f), P2O5 shows a decreasing trend with increasing SiO2, corresponding to I-type granite [54]. In addition, Rb versus Y and Rb versus Th discriminate diagrams also suggest an I-type granitic trend for the studied samples in this study (Figure 8g,h). In addition, low zircons δ18O values of 5.12–6.24‰ exclude the possibility of S-type granite (Figure 5e). It is well known that granites sourced from partial melting of sedimentary rocks or reworking of S-type granites will have high zircon δ18O values (>8‰), indicating a supra-crustal interaction [56,57,58]. In addition, microstructural analyses revealed that the granitic samples correspond with I-type granite, comprising plagioclase, K-feldspar, quartz, and biotite (Figure 2c–e). Compared to I-type granite, S-type granite is rich in high-aluminum (Al) minerals (e.g., muscovite, cordierite, and garnet). In conclusion, the granitic pluton identified in the Caochang region is composed of typical I-type granite that originated from the partial melting of meta-igneous rocks.

5.2. Magma Source and Petrogenesis

Multiple studies have demonstrated that the granites are products of partial melting of the Yangtze crustal material, as suggested by low values of oxygen isotopes and Neoproterozoic (600–800 Ma) inherited zircons in the Sulu orogenic belt [4,5,11]. However, magmatic rocks collected from the North China craton (NCC) have high values of oxygen isotopes and Archean inherited zircons. The Wulian (Baichihe)–Qingdao-Yantai fault (WQYF) is regarded as the boundary of the Yangtze and the NCC. It is clear that granites collected from the west of the WQYF show a genetic association with partial melting of the NCC crustal material, whereas those from the east of the WQYF have a close relation with the Yangtze craton [4,27,28]. However, the nature of granites from the fault zone and neighboring regions (tectonic amalgamation region) is poorly understood. Recently, Zhao [59] provided a summary and proposed that the granites across the fault zone and neighboring regions all belong to the Yangtze craton. This statement implies that granites distributed in the Sulu belt and boundary region have a close affinity with the Yangtze craton. According to Zhao’s [59] proposal, our samples, which are located in the east of the fault zone (Figure 1b), have a close affinity with the Yangtze craton. Additionally, the absence of the Archean inherited zircons also rules out the possibility of the NCC (Figure 4).
The zircon Lu–Hf isotope is a useful tool in discriminating magma sources and petrogenetic processes [60]. The analytical results revealed that the granitic samples have homogeneous Hf isotopic compositions, with narrow variations of εHf(t) values ranging from −27.51 to −32.35. In addition, most of the zircons have low values of 176Lu/177Hf (Figure 5a) (<0.002), suggesting a good quality of data [61]. It is well established that positive εHf(t) values suggest a crustal growth or recycling of juvenile crust, and negative εHf(t) values represent ancient crustal recycling or reworking. The newly identified granite showed strongly negative εHf(t) values (−27.51 to −32.35, mean value: −30.21) (Figure 5b,c), indicating that the magma source of the studied pluton was derived from partial melting of ancient crust. The two-stage Hf model ages ranged from 2979 to 3175 Ma (mean age: 3053 Ma) (Figure 5d), also indicating partial melting of the Mesoarchean crustal material. In conclusion, we argue that the granitic pluton is derived from partial melting of the Yangtze ancient crustal material. The strongly negative εHf(t) values suggest multiple reworking of the magma source.
As described above, the granitic pluton originated from partial melting of the Yangtze crustal material. The effective geochemical diagrams indicate that the newly identified pluton was derived from partial melting of mafic material (Figure 9a,b). Compared to the middle-upper crust, the lower crust is mafic in composition [62,63]. Considering the previous argument, we argue that the granitic pluton originated due to the partial melting of the Yangtze mafic lower crust. Zhao [4] concluded that most of the Early Cretaceous granites occurring in the Sulu belt had normal oxygen (δ18O) values close to those of normal mantle zircons (5.3 ± 0.3‰), suggesting a deep magma source with no significant water–rock interaction. If the rocks were derived from the middle-upper crust, they would experience significant influence of atmospheric precipitation, demonstrating low oxygen isotopes. In this study, our samples had δ18O values of 5.12–6.25‰ that were consistent with juvenile basaltic crust (5.5–6.5‰) being formed by partial melting of mantle-derived magma [64] (δ = 18Ozrn = 5.3 ± 0.3‰) (Figure 5e), suggesting that the granite was probably derived from partial melting of basaltic lower crust and experienced negligible supra-crustal contamination in the parental magma. This was also indicated by coeval felsic volcanic rocks (δ18O = 4.82–5.86‰) on Lingshan Island [11].
Compared to zircon oxygen isotopic information, zircon Hf isotope is also a necessary requirement for constraining the formation and evolution of continental crust [60,65]. The combination of Hf–O isotopes is a very powerful tool for deciphering the nature of magma and the magmatic evolution process of host granites and mafic micro-granular enclaves [57,66]. Unlike the δ18O results, the εHf(t) values of the studied pluton displayed ancient crustal characteristics (Figure 5b, εHf(t) = −27.51~−32.35). Furthermore, the zircon Hf two-stage model ages (TDM2) showed Mesoarchean rocks from the Yangtze craton [67,68] (Figure 5d). Therefore, the studied zircons have mantle-derived oxygen isotopic information, whereas Hf isotopic compositions suggest ancient crustal features (Figure 5f). Under most circumstances, significant negative εHf(t) values are associated with crustal material input, especially the supra-crustal sediments. In contrast, the zircon δ18O isotopic values are very sensitive to supra-crustal material input. Crustal contamination will lead to sharp deviations from the mantle-like δ18O values (Figure 5f). As a consequence, the correlation (Figure 5f) is not suitable via the assimilation fractional crystallization (AFC) process.
A tenable explanation for the correlation of δ18O and εHf(t) values is that the Caochang granitic pluton potentially originated from partial melting of an ancient mafic lower crust. It has been demonstrated that a mafic lower crust derived from a depleted mantle can show enriched Hf isotopic compositions through a considerable period of time [60,69,70]. If supra-crustal materials were not injected and did not experience low- and high-temperature water–rock interactions, the mafic lower crust would have negligible variations of oxygen isotope values that retain the mantle-like δ18O characteristics. Based on the above discussion, we propose that the Caochang granitic pluton derived from the partial melting of the mafic lower crust. This model is in agreement with the geochemical discrimination diagrams that suggest partial melting of the mafic lower crust (Figure 9a,b).
Furthermore, experimental petrology demonstrated that partial melting of mafic or basaltic lower crust could not produce high Mg-number (Mg/Mg + TFe) magmas [71]. Therefore, high Mg-number magmas suggest injection of mantle-derived magma, excluding partial melting of the pure crust (Figure 9c). Moreover, the granite pluton shows an arc-type geochemical feature characterized by enrichment of light REE and Pb as well as Nb–Ta–Ti depletion (Figure 7a,b), demonstrating that the granite pluton originated from the recycling of arc-related rocks and preserved arc-related features. Although supra-crustal materials also display arc-type geochemical features, this possibility is precluded owing to mantle-like oxygen characteristics (Figure 5e,f). This is also supported by the Nb versus Nb/Th diagram (Figure 9d).

5.3. Mineral Fractional Crystallization

In Figure 10, TiO2 and MgO demonstrate non-linear trends with the variations of SiO2 contents, suggesting that mafic minerals (e.g., amphibole, pyroxene, and biotite) have a minor role in the magma evolution (Figure 10a,c). Although minor amounts of biotite were found in the thin sections and field exposures (Figure 2e), biotite is usually considered as volatile, representing a later magma evolution process. Na2O, TFe2O3, Al2O3, and CaO showed good linear relations with SiO2 (R2 > 0.5) (Figure 10b,d–f). It is notable that Na2O and Al2O3 showed decreasing trends with increasing SiO2, suggesting plagioclase and alkali feldspar crystallization fractionation (Figure 10b,e). However, steady linear relations, as suggested by TFe2O3 and CaO, suggest the absence of basic feldspar and mafic minerals in the evolution process (Figure 10d,f). In Figure 11, trace element discrimination diagrams correspond well with the Harker diagrams (Figure 10). Figure 11a,b also suggests that the combination of K-feldspar and plagioclase played a key role in the granitic pluton during the magma evolution process. In addition, Figure 11c,d indicates monazite crystallization and fractionation. However, Zr/Hf ratios had consistent ratios, ranging from 39 to 42 (Table 4), suggesting that zircon fractional crystallization was negligible [76,77].

5.4. Evaluation on the Adakite-Like Rocks

Previous studies revealed that Late Jurassic granites show high Sr/Y and (La/Yb)N ratios, suggesting that they originated from partial melting of thickened crust. As mentioned above, Early Cretaceous granites show complicated geochemical features. Some granitic plutons, distributed in northeastern segment of the Jiaodong peninsula, have high Sr/Y and (La/Yb)N ratios that are consistent with adakite-like rocks [79,80,81], whereas the other granites are characterized by low Sr/Y and (La/Yb)N ratios in the Sulu belt. Compared to the Late Jurassic adakite-like rocks, the petrogenesis and setting of the Early Cretaceous adakite-like rocks are poorly understood. In this study, our samples show arc-related geochemical characteristics (Figure 12a). However, the YbN versus (La/Yb)N diagram indicates an adakite-like feature corresponding to thickened crust (Figure 12b). It is well known that eastern China experienced significant lithospheric thinning and foundering [11,50,82], so it is untenable that the Early Cretaceous adakite-like rocks were formed in a thickened setting.
It has been attested that adakite-like rocks, as defined by Sr/Y and (La/Yb)N discriminant diagrams, can be formed in different settings: (1) Partial melting of thickened crust or foundering crust [83,84,85]; (2) Partial melting of the oceanic slab [86]; (3) Crystallization and differentiation [87,88]; (4) A high Sr/Y magma source [43]; (5) Magma mixing of different end-members [89,90]; and (6) Partial melting of granulite rocks [91]. It is generally accepted that (1) and (2) occur for true adakite rocks, while the other types (3 to 6) occur for pseudo-adakites [85,92]. He [93] proposed new discriminant diagrams that can better discriminate true and pseudo adakites. Figure 12c–f reveals that the samples plot into the normal granitic fields rather than normal genetic adakites. Combined with regional background and geochemical characteristics of early Cretaceous granites, we argue that early Cretaceous granites originated from partial melting of normal crustal material rather than thickened crust in the Sulu belt. However, our samples are characterized by high (La/Yb)N ratios that might be related to magma mixing, as suggested by magma mixing lines [78] (Figure 10g,h).

5.5. Tectonic Implications

It is well known that underplating of mantle-derived, basaltic magma at the crust–mantle boundary might be an important mechanism for crustal growth and recycling [71,95]. The significantly negative εHf(t) values indicate multiple crustal recycling for the studied pluton. It was also proved that mantle-derived material and basaltic magma supplied sufficient heat to result in partial melting of crust material at the crust–mantle boundary. During the interaction process, the mantle material probably participates in the formation of granitic rocks. Multiple studies have demonstrated that mafic microgranular enclaves (MME) and high-Mg-number granites reflect input and participation of mantle material [71,92].
MMEs and the high-Mg-number values (>40) reveal that the mantle material probably participated in the evolution of the granitic pluton in study area (Figure 2a–c and Figure 9c). Multiple studies have shown that formational processes of MMEs can be classified as follows: (1) fragments or xenoliths of mafic country rocks [96]; (2) refractory solid residues [97,98]; (3) early products from co-genetic magma or mafic mineral cumulates [99]; and (4) mixed mafic–felsic magmas in the chamber [92,100]. Compared to other models, mafic–felsic magmas mixing is the dominant model in producing high-Mg-number granites. Based on microtectonics and geochemical features, we argue that magma mixing is the best candidate for forming the granitic pluton with the MMEs in the Sulu belt.
Multiple studies have proven that the eastern China passed through significant lithosphere thinning during the Mesozoic. Rifted basins and multiple extensional tectonics (e.g., metamorphic core complex) are the shallow responses of the lithosphere thinning in eastern China [11,101,102]. In addition, Early Cretaceous large-scale mineralization and magmatic flare-up also suggest tectonic mechanism transfer during the Jurassic to Early Cretaceous [11,103]. The Late Jurassic granitoid rocks mainly originated from partial melting of pure crust material, with no significant mantle material input in the Sulu belt. However, the Early Cretaceous granitoid rocks have complex petrogenetic mechanisms, including pure crust-derived granitoids, crust–mantle mixing-derived granitoids, and large amounts of intermediate-mafic dykes. Particular rock assemblages and extensional tectonics suggest that eastern China experienced large-scale extension and thinning during the Early Cretaceous. Zhu [104] proposed that the Pacific plate subduction beneath the eastern China was the main factor causing lithosphere destruction and thinning of eastern China during the Early Cretaceous. Furthermore, Zhu [104] reported that the stagnant slab lying in the mantle transition zone resulted in an unsteady mantle flow (convective asthenosphere), leading to overlying lithosphere mantle thinning and destruction. Combined above field observations and geochemical and Hf–O isotopes, we propose that underplating of mantle-derived magma was the main factor leading to partial melting of the Yangtze mafic lower crust in the Sulu orogenic belt. During the evolution process, the mantle-derived magma also participated in the formation of the Caochang granitic pluton. Mantle-like δ18O values reveal that the Caochang granitic pluton originated from a deep crustal setting, with no water–rock interaction and strong tectonic hot or cold events.

6. Concluding Remarks

(1) The new identified granitic pluton was crystallized during the Early Cretaceous (ca. 120 Ma), as suggested by zircon U–Pb dating.
(2) In-situ zircon Hf-O isotope data from the Caochang granitic pluton demonstrate a decoupling correlation of significant negative εHf(t) and mantle-like δ18O values, suggesting a basaltic lower crust origin with no supra-crustal material input and water–rock interaction.
(3) High Mg-number values and mafic micro-granular enclaves suggest mantle material participation during the granitic magma evolution.
(4) In this study, the newly identified granitic pluton was the product of the destruction of the NCC. The NCC destruction resulted in decompression and mantle-derived magma upwelling that led to partial melting of Yangtze lower crust material, forming the Caochang granitic pluton.

Author Contributions

Y.M. carried out a detailed geological survey and wrote this manuscript; Z.W. provided ideas and experimental funds; B.G. made a contribution to tectonic models and geochemical discrimination diagrams; J.L. made a contribution to experimental measurements and data processing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Young innovative projects of Shandong Province (2019KJH004), the National Natural Science Foundation of China (41902230, 41876059), and the Shandong Provincial Natural Science Foundation (ZR2019QD002).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, J.H.; Chung, S.L.; Wilde, S.A.; Wu, F.Y.; Chu, M.F.; Lo, Q.H.; Fan, H.R. Petrogenesis of post-orogenic syenites in the Sulu Orogenic Belt, East China: Geochronological, geochemical and Nd-Sr isotopic evidence. Chem. Geol. 2005, 214, 99–125. [Google Scholar] [CrossRef]
  2. Yang, J.H.; Wu, F.Y.; Chung, S.L.; Wilde, S.A.; Chu, M.F.; Lo, C.H.; Song, B. Petrogenesis of Early Cretaceous intrusions in the Sulu ultrahigh-pressure orogenic belt, East China and their relationship to lithospheric thinning. Chem. Geol. 2005, 222, 200–231. [Google Scholar] [CrossRef]
  3. Liu, F.L.; Liu, L.S.; Liu, P.H.; Wang, F.; Cai, J.; Liu, J.H.; Wang, W.; Ji, L. A relic slice of Archean-early Paleoproterozoic basement of Jiaobei Terrane identified within the Sulu UHP belt: Evidence from protolith and metamorphic ages from meta-mafic rocks, TTG-granitic gneisses, and metasedimentary rocks in the Haiyangsuo region. Precambrian Res. 2017, 303, 117–152. [Google Scholar] [CrossRef]
  4. Zhao, Z.F.; Liu, Z.B.; Chen, Q. Melting of subducted continental crust: Geochemical evidence from Mesozoic granitoids in the Dabie-Sulu orogenic belt, east-central China. J. Asian Earth Sci. 2017, 145, 260–277. [Google Scholar] [CrossRef]
  5. Zhao, Z.F.; Zheng, Y.F. Remelting of subducted continental lithosphere: Petrogenesis of Mesozoicmagmatic rocks in the Dabie-Sulu orogenic belt. Sci. China Ser. D 2009, 52, 1295–1318. [Google Scholar] [CrossRef]
  6. Xu, Y.; Yang, M.L.; Polat, A.; Yang, K.G.; Meng, Y.K.; Li, R.H.; Wang, Z.B. Zircon U-Pb ages and Hf isotopic compositions of two types of supracrustal rocks in the Northeastern Sulu UHP terrane: Constraints on the surface boundary between South China and North China. Int. Geol. Rev. 2018, 61, 539–561. [Google Scholar] [CrossRef]
  7. Wang, S.J.; Schertl, H.P.; Pang, Y.M. Geochemistry, geochronology and Sr-Nd-Hf isotopes of two types of early cretaceous granite porphyry dykes in the Sulu Orogenic Belt, Eastern China. Can. J. Earth Sci. 2019, 57, 249–266. [Google Scholar] [CrossRef]
  8. Zheng, Y.F. A perspective view on ultrahigh-pressure metamorphism and continental collision in the Dabie-Sulu orogenic belt. Chin. Sci. Bull. 2008, 53, 3081–3104. [Google Scholar] [CrossRef] [Green Version]
  9. Liou, J.G.; Ernst, W.G.; Zhang, R.Y.; Tsujimori, T.; Jahn, B.M. Ultrahigh-pressure minerals and metamorphic terranes-the view from China. J. Asian Earth Sci. 2009, 35, 199–231. [Google Scholar] [CrossRef]
  10. Zheng, Y.F.; Fu, B.; Gong, B.; Li, L. Stable isotope geochemistry of ultrahigh pressure metamorphic rocks from the Dabie–Sulu orogeny in China: Implications for geodynamics and fluid regime. Earth Sci. Rev. 2003, 62, 105–161. [Google Scholar] [CrossRef]
  11. Meng, Y.K.; Santosh, M.; Li, R.H.; Xu, Y.; Hou, F.H. Petrogenesis and tectonic implications of Early Cretaceous volcanic rocks from Lingshan Island in the Sulu Orogenic Belt. Lithos 2018, 312–313, 244–257. [Google Scholar] [CrossRef]
  12. Wang, H.; Wang, S.H.; Xu, Z.W.; Fu, B.; Zhao, Z.X.; Li, Z.G.; Dong, Y.H.; Tang, L.M.; Li, J. Geochemical and Sr-Nd-Pb-Hf-O isotopic compositions of the Tiezhai complex: Implications for lithosphere destruction of the North China Craton. Gondwana Res. 2018, 61, 203–221. [Google Scholar] [CrossRef]
  13. Wang, J.; Chang, S.C.; Chen, Y.; Yan, S.Y. Early Cretaceous transpressional and transtensional tectonics straddling the Sulu orogenic belt, East China. Geosci. Front. 2019, 10, 2287–2300. [Google Scholar] [CrossRef]
  14. Xu, H.J.; Zhang, J.F.; Wang, Y.F.; Liu, W.L. Late Triassic alkaline complex in the Sulu UHP terrane: Implication for post-collisional magmatism and subsequent fractional crystallization. Gondwana Res. 2016, 35, 390–410. [Google Scholar] [CrossRef] [Green Version]
  15. Zhao, Z.F.; Zheng, Y.F.; Zhang, J.; Dai, L.Q.; Liu, X.M. Syn-exhumation magmatism during continental collision: Evidence from alkaline intrusives of Triassic age in the Sulu orogen. Chem. Geol. 2012, 328, 70–88. [Google Scholar] [CrossRef]
  16. Jahn, B.M.; Wu, F.Y.; Lo, C.H.; Tsai, C.H. Crust-mantle interaction induced by deep subduction of the continental crust: Geochemical and Sr-Nd isotopic evidence from post-collisional mafic-ultramafic intrusions of the northern Dabie complex, central China. Chem. Geol. 1999, 157, 119–146. [Google Scholar] [CrossRef]
  17. Zhou, X.M.; Li, W.X. Origin of Late Mesozoic igneous rocks of southeastern China: Implications for lithosphere subduction and underplating of mafic magma. Tectonophysics 2000, 326, 269–287. [Google Scholar] [CrossRef]
  18. Zhao, Z.F.; Zheng, Y.F.; Chen, Y.X.; Sun, G.C. Partial melting of subducted continental crust: Geochemical evidence from synexhumation granite in the Sulu orogen. Geol. Soc. Am. Bull. 2017, 129, 1692–1707. [Google Scholar] [CrossRef]
  19. Niu, Y.L.; Liu, Y.; Xue, Q.Q.; Shao, F.L.; Chen, S.; Duan, M.; Guo, P.Y.; Gong, H.M.; Hu, Y.; Hu, Z.X.; et al. The exotic origin of the continental shelf of East and South China Seas-a testable hypothesis with consequences on the tectonic evolution of the western Pacific and eastern China since the Mesozoic. Sci. Bull. 2015, 18, 1598–1616. [Google Scholar] [CrossRef] [Green Version]
  20. Gao, Y.J.; Niu, Y.L.; Duan, M.; Xue, Q.Q.; Sun, P.; Chen, S.; Xiao, Y.Y.; Guo, P.Y.; Wang, X.H.; Chen, Y.H. The petrogenesis and tectonic significance of the Early Cretaceous intraplate granites in eastern China: The Laoshan granite as an example. Lithos 2019, 328–329, 200–211. [Google Scholar] [CrossRef] [Green Version]
  21. Dai, L.Q.; Zheng, Y.F.; Zhao, Z.F. Termination time of peak decratonization in North China: Geochemical evidence from mafic igneous rocks. Lithos 2016, 240–243, 327–336. [Google Scholar] [CrossRef]
  22. Ling, M.X.; Li, Y.; Ding, X.; Teng, F.Z.; Yang, X.Y.; Fan, W.M.; Xu, Y.G.; Sun, W.D. Destruction of the North China Craton induced by ridge subductions. J. Geol. 2013, 121, 197–213. [Google Scholar] [CrossRef]
  23. Yan, Q.S.; Shi, X.F. Geochemistry and petrogenesis of the Cretaceous A-type granites in the Laoshan granitic complex, eastern China. Island Arc 2014, 23, 221–235. [Google Scholar] [CrossRef]
  24. Zhang, J.; Zhao, Z.F.; Zheng, Y.F.; Dai, M.N. Postcollisional magmatism: Geochemical constraints on the petrogenesis of Mesozoic granitoids in the Sulu orogen, China. Lithos 2010, 119, 512–536. [Google Scholar] [CrossRef]
  25. Zhang, J.; Zhao, Z.F.; Zheng, Y.F.; Liu, X.M.; Xie, L.W. Zircon Hf–O isotope and whole-rock geochemical constraints on origin of postcollisional mafic to felsic dykes in the Sulu orogen. Lithos 2012, 136, 225–245. [Google Scholar] [CrossRef]
  26. Meng, Y.K.; Li, R.H.; Xu, Y.; Hou, F.H. U-Pb-Hf isotopes and tectonic significance of early Cretaceous detrital zircons on Lingshan Island, Qingdao of Shandong Province. Earth Sci. 2018, 43, 3302–3323. [Google Scholar]
  27. Zhou, J.B.; Wilde, S.A.; Zhao, G.C.; Zhang, X.Z.; Zheng, C.Q.; Jin, W.; Cheng, H. SHRIMP U–Pb zircon dating of theWulian complex: Defining the boundary between the north and south China cratons in the Sulu orogenic belt, China. Precambrian Res. 2008, 162, 559–576. [Google Scholar] [CrossRef]
  28. Zhou, J.B.; Wilde, S.; Zhao, G.C.; Zheng, C.Q.; Jin, W.; Zhang, X.Z.; Cheng, H. SHRIMP U-Pb zircon dating of the Neoproterozoic Penglai Group and Archean gneisses from the Jiaobei terrane, north China, and their tectonic implications. Precambrian Res. 2008, 160, 323–340. [Google Scholar] [CrossRef]
  29. Guo, J.H.; Chen, F.K.; Zhang, X.M.; Siebel, W.; Zhai, M.G. Evolution of syn- to post-collisional magmatism from north Sulu UHP belt, eastern China: Zircon U-Pb geochronology. Acta Petrol. Sin. 2005, 21, 1281–1301. [Google Scholar]
  30. Liu, Y.S.; Hu, Z.C.; Zong, K.Q.; Gao, C.G.; Gao, S.; Xu, J.; Chen, H.H. Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS. Chin. Sci. Bull. 2010, 55, 1535–1546. [Google Scholar] [CrossRef]
  31. Zong, K.Q.; Liu, Y.S.; Gao, C.G.; Hu, Z.H.; Gao, S.; Gong, H.J. In situ U-Pb dating and trace element analysis of zircons in thin sections of eclogite: Refining constraints on the UHP metamorphism of the Sulu terrane, China. Chem. Geol. 2010, 269, 237–251. [Google Scholar] [CrossRef]
  32. Williams, I.S. U-Th-Pb geochronology by ion microprobe. In Applications of Microanalytical Techniques to Understanding Mineralizing Processes; McKibben, M.A., Shanks, W.C., Ridley, W.I., Eds.; Review of Economic Geology; Society of Economic Geologists: Littleton, CO, USA, 1998; Volume 7, pp. 1–35. [Google Scholar]
  33. Black, L.P.; Kamo, S.L.; Allen, C.M.; Davis, D.W.; Aleinikoff, J.N.; Valley, J.W.; Mundil, R.; Campbell, I.H.; Korsch, R.J.; Williams, I.S.; et al. Improved 206Pb/238U microprobe geochronology by the monitoring of a trace-element related matrix effect; SHRIMP, ID-TIMS, ELA-ICP-MS and oxygen isotope documentation for a series of zircon standards. Chem. Geol. 2004, 205, 115–140. [Google Scholar] [CrossRef]
  34. Hu, Z.C.; Liu, Y.S.; Gao, S.; Liu, W.G.; Yang, L.; Zhang, W.; Tong, X.R.; Lin, L.; Zong, K.Q.; Li, M.; et al. Improved in situ Hf isotope ratio analysis of zircon using newly designed X skimmer cone and Jet sample cone in combination with the addition of nitrogen by laser ablation multiple collector ICP-MS. J. Anal. At. Spectrom. 2012, 27, 1391–1399. [Google Scholar] [CrossRef]
  35. Hoskin, P.W.O.; Schaltegger, U. The composition of zircon and igneous and metamorphic petrogenesis. Rev. Mineral. Geochem. 2003, 53, 27–62. [Google Scholar] [CrossRef]
  36. Andersen, T. Correction of common lead in U-Pb analyses that do not report 204Pb. Chem. Geol. 2002, 192, 59–79. [Google Scholar] [CrossRef]
  37. Xiong, F.H.; Meng, Y.K.; Yang, J.S.; Liu, Z.; Xu, X.Z.; Eslami, A.; Zhang, R. Geochronology and petrogenesis of the mafic dykes from the Purang ophiolite: Implications for evolution of the western Yarlung-Tsangpo suture zone, Tibet. Geosci. Front. 2020, 11, 277–292. [Google Scholar] [CrossRef]
  38. Middlemost, E.A.K. Naming materials in the magma/igneous rock system. Earth-Sci. Rev. 1994, 37, 215–224. [Google Scholar] [CrossRef]
  39. Irvine, T.N.; Baragar, W.R.A. A guide to the chemical classification of the common volcanic rocks. Can. J. Earth Sci. 1971, 8, 523–548. [Google Scholar] [CrossRef]
  40. Wright, J.B. A simple alkalinity ratio and its application to questions of non-orogenic granite genesis. Geol. Mag. 1969, 106, 370–384. [Google Scholar] [CrossRef]
  41. Peccerillo, A.; Taylor, S.R. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, northern Turkey. Contrib. Miner. Petrol. 1976, 58, 63–81. [Google Scholar] [CrossRef]
  42. Maniar, P.D.; Piccoli, P.M. Tectonic discrimination of granitoids. Geol. Soc. Am. Bull. 1989, 101, 635–643. [Google Scholar] [CrossRef]
  43. Moyen, J.F. High Sr/Y and La/Yb ratios: The meaning of the adakitic signature. Lithos 2009, 112, 556–574. [Google Scholar] [CrossRef]
  44. Sun, T.; Zhou, X.M.; Chen, P.R.; Li, H.M.; Zhou, H.Y.; Wang, Z.C.; Shen, W.Z. Mesozoic strong peraluminous granite and its tectonic significance in the eastern Nanling Mountains, China. Sci. China Ser. D 2003, 33, 1209–1218. [Google Scholar]
  45. Boynton, W.V. Geochemistry of the rare earth elements: Meteorite studies. In Rare Earth Element Geochemistry; Henderson, P., Ed.; Elsevier: Amsterdam, The Netherlands, 1984; pp. 63–114. [Google Scholar]
  46. Sun, S.S.; McDonough, W.F. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In Magmatism in the Ocean Basins; Saunders, A.D., Norry, M.J., Eds.; Geological Society of London Special Publications: London, UK, 1989; Volume 42, pp. 313–345. [Google Scholar]
  47. Chappell, B.W.; White, A.J.R. Two contrasting granite types. Pac. Geol. 1974, 8, 173–174. [Google Scholar]
  48. Loiselle, M.C.; Wones, D.R. Characteristics and origin of anorogenic granites. Geol. Soc. Am. Abstr. 1979, 11, 468. [Google Scholar]
  49. Hao, L.L.; Wang, Q.; Wyman, D.A.; Ma, L.; Wang, J.; Xia, X.P.; Qu, Q. First identification of postcollisional A-type magmatism in the Himalayan-Tibetan orogen. Geology 2019, 47, 187–190. [Google Scholar] [CrossRef]
  50. Zhu, R.X.; Yang, J.H.; Wu, F.Y. Timing of destruction of the North China Craton. Lithos 2012, 149, 51–60. [Google Scholar] [CrossRef]
  51. Watson, E.B.; Wark, D.A.; Thomas, J.B. Crystallization thermometers for zircon and rutile. Contrib. Mineral. Petrol. 2006, 151, 413–433. [Google Scholar] [CrossRef]
  52. Zhao, G.T.; Wang, D.Z.; Cao, Q.C. The geochemistry and genesis of the Laoshan Granitoids, Shandong Province. Geol. J. China Univ. 1997, 3, 1–15. [Google Scholar]
  53. Whalen, J.B.; Currie, K.L.; Chappell, B.W. A-type granites: Geochemical characteristics, discrimination and petrogenesis. Contrib. Miner. Petrol. 1987, 95, 407–419. [Google Scholar] [CrossRef]
  54. Chappell, B.W.; White, A.J.R. I- and S-type granites in the Lachlan Fold Belt. Transactions of the Royal Society of Edinburgh: Earth Sci. 1992, 83, 1–26. [Google Scholar] [CrossRef]
  55. Chappell, B.W. Aluminium saturation in I- and S-type granites and the characterization of fractionated haplogranites. Lithos 1999, 46, 535–551. [Google Scholar] [CrossRef]
  56. Zhang, S.B.; Zheng, Y.F. Time and space of Neoproterozoic low δ18O magmatic rocks in South China. Chin. Sci. Bull. 2013, 58, 2344–2350. [Google Scholar]
  57. Gao, P.; Zheng, Y.F.; Zhao, Z.F. Distinction between S-type and peraluminous I-type granites: Zircon versus whole-rock geochemistry. Lithos 2016, 258, 77–91. [Google Scholar] [CrossRef]
  58. Lu, Y.H.; Gao, P.; Zhao, Z.F.; Zheng, Y.F. Whole-rock geochemical and zircon Hf-O isotopic constraints on the origin of granitoids and their mafic enclaves from the Triassic Mishuling pluton in West Qinling, central China. J. Asian Earth Sci. 2019, 189, 104136. [Google Scholar] [CrossRef]
  59. Zhao, R.; Wang, Q.F.; Liu, X.F.; Wang, W.; Pan, R.G. Architecture of the Sulu crustal suture between the North China Craton and Yangtze Craton: Constraints from Mesozoic granitoids. Lithos 2016, 266, 348–361. [Google Scholar] [CrossRef]
  60. Wu, F.Y.; Li, X.H.; Zheng, Y.F.; Gao, S. Lu-Hf isotopic systematics and their applications in petrology. Acta Petrol. Sin. 2007, 23, 185–220. [Google Scholar]
  61. Griffin, W.L.; Wang, X.; Jackson, S.E.; Pearson, N.J.; O’Reilly, S.Y.; Xu, X.; Zhou, X. Zircon chemistry and magmamixing, SE China: In-situ analysis of Hf isotopes, Tonglu and Pingtan igneous complexes. Lithos 2002, 61, 237–269. [Google Scholar] [CrossRef]
  62. Dawson, J.B.; Carswell, D.A.; Hall, J.; Wedepohl, K.H. The nature of the Lower Continental Crust; Geological Society Special Publication: London, UK, 1986; Volume 24, pp. 363–374. [Google Scholar]
  63. Rudnick, R.L.; Gao, S. Composition of the continental crust. Treatise Geochem. 2003, 3, 1–64. [Google Scholar]
  64. Valley, J.W.; Lackey, J.S.; Cavosie, A.J.; Clechenko, C.C.; Spicuzza, M.J.; Basei, M.A.S.; Bindeman, I.N.; Ferreira, V.P.; Sial, A.N.; King, E.M.; et al. 4.4 billion years of crustal maturation: Oxygen isotope ratios of magmatic zircon. Contrib. Mineral. Petrol. 2005, 150, 561–580. [Google Scholar] [CrossRef]
  65. Ma, Q.; Xu, Y.G.; Huang, X.L.; Zheng, J.P.; Ping, X.Q.; Xia, X.P. Eoarchean to Paleoproterozoic crustal evolution in the North China Craton: Evidence from U-Pb and Hf-O isotopes of zircons from deep-crustal xenoliths. Geochimica et Cosmochimica Acta. 2019, 278, 94–109. [Google Scholar] [CrossRef]
  66. Kemp, A.I.S.; Hawkesworth, C.J.; Foster, G.L.; Paterson, B.A.; Woodhead, J.D.; Hergt, J.M.; Gray, C.M.; Whitehouse, M.J. Magmatic and crustal differentiation history of granitic rocks from Hf-O isotopes in zircon. Science 2007, 315, 980–983. [Google Scholar] [CrossRef]
  67. Gao, S.; Yang, J.; Zhou, L.; Li, M.; Hu, Z.; Guo, J.; Yuan, H.; Gong, H.; Xiao, G.; Wei, J. Age and growth of the Archean Kongling terrain, South China, with Emphasis on 3.3 Ga granitoid gneisses. Am. J. Sci. 2011, 311, 153–182. [Google Scholar] [CrossRef]
  68. Guo, J.L.; Wu, Y.B.; Gao, S.; Jin, Z.M.; Zong, K.Q.; Hu, Z.C.; Chen, K.; Chen, H.H.; Liu, Y.S. Episodic Paleoarchean-Paleoproterozoic (3.3–2.0 Ga) granitoid magmatism in Yangtze craton, South China: Implications for Late Archean Tectonics. Precambrian Res. 2015, 270, 246–266. [Google Scholar] [CrossRef]
  69. Taylor, S.R.; Campbell, I.H.; McCulloch, M.T.; McLennan, S.M. A lower crustal origin for massif-type anorthosites. Nature 1984, 311, 372. [Google Scholar] [CrossRef]
  70. He, H.L.; Song, X.Y.; Zhai, M.G.; Yu, S.Y.; Du, Z.S. Lower crustal contribution to the magma formation of the Damiao massif-type anorthosite, North China Craton: Evidence from zircon Hf-O isotopes. Precambrian Res. 2019, 332, 105396. [Google Scholar] [CrossRef]
  71. Rapp, R.P.; Watson, E.B. Dehydration melting of metabasalt at 8–32 kbar: Implications for continental growth and crust-mantle recycling. J. Petrol. 1995, 36, 891–931. [Google Scholar] [CrossRef]
  72. Patiño Douce, A.E. What do experiments tell us about the relative contributions of crust and mantle to the origin of granitic magmas? In Understanding Granites: IntergratingNew and Classical Techniques; Castro, A., Fernandez, C., Vigneresse, J.L., Eds.; Special Publications, Geological Society London: London, UK, 1999; Volume 168, pp. 55–75. [Google Scholar]
  73. Altherr, R.; Holl, A.; Hegner, E. High-potassium, calc-alkaline I-type plutonism in the European Variscides: Northern Vosges (France) and northern Schwarzwald (Germany). Lithos 2000, 50, 51–73. [Google Scholar] [CrossRef]
  74. Hofmann, A.W. Chemical differentiation of the Earth: The relationship between mantle, continental crust, and oceanic crust. Earth Planet. Sci. Lett. 1988, 90, 297–314. [Google Scholar] [CrossRef] [Green Version]
  75. Schmidberger, S.S.; Hegner, E. Geochemistry and isotope systematics of calc-alkaline volcanic rocks from the Saar-Nahe basin (SW Germany): Implications for Late-Variscan orogenic development. Contrib. Mineral. Petrol. 1999, 135, 373–385. [Google Scholar] [CrossRef]
  76. Dostal, J.; Chatterjee, A.K. Constrasting behaviour of Nb/Ta and Zr/Hf ratios in a peraluminous granitic pluton (Nova Scotia, Canada). Chem. Geol. 2000, 163, 207–218. [Google Scholar] [CrossRef]
  77. Linnen, R.L.; Keppler, H. Columbite solubility in granitic melts: Consequence for the enrichment and fractionation of Nb and Ta in the Earth’s crust. Contrib. Mineral. Petrol. 1997, 128, 213–227. [Google Scholar] [CrossRef]
  78. Langmuir, C.H.; Vocke, R.D.; Hanson, G.N.; Hart, S.R. A general mixing equation with applications to Icelandic basalts. Earth Planet. Sci. Lett. 1978, 37, 380–392. [Google Scholar] [CrossRef]
  79. Yang, J.H.; Zhu, M.F.; Liu, W.; Zhai, M.G. Geochemistry and petrogenesis of Guojialing granodiorites from the northwestern Jiaodong Peninsula, eastern China. Acta Petrol. Sin. 2003, 19, 692–700. [Google Scholar]
  80. Zhang, H.F.; Zhai, M.G.; Tong, Y.; Peng, P.; Xu, B.L.; Guo, J.H. Petrogenesis of the Sanfoshan high-Ba-Sr granite, Jiaodong Peninsula, eastern China. Geol. Rev. 2006, 52, 43–53. [Google Scholar]
  81. Du, L.; Zhang, Y.Y.; Huang, Z.Y.; Li, X.-P.; Yuan, C.; Wu, B.; Long, X.P. Devonian to Carboniferous tectonic evolution of the Kangguer Ocean in the Eastern Tianshan, NW China: Insights from three episodes of granitoids. Lithos 2019, 350, 105243. [Google Scholar] [CrossRef]
  82. Liu, J.G.; Cai, R.H.; Pearson, D.G.; Scott, J.M. Thinning and destruction of the lithospheric mantle root beneath the North China Craton: A review. Earth Sci. Rev. 2019, 196, 102873. [Google Scholar] [CrossRef]
  83. Petford, N.; Atherton, M. Na–rich partial melts from newly underplated basaltic crust: The Cordillera Blanca Batholith. Peru. J. Petrol. 1996, 37, 1491–1521. [Google Scholar] [CrossRef] [Green Version]
  84. Xu, J.F.; Shinjo, R.; Defant, M.J.; Wang, Q.; Rapp, R.P. Origin of Mesozoic adakitic intrusive rocks in the Ningzhen area of east China: Partialmelting of delaminated lower continental crust? Geology 2002, 30, 1111–1114. [Google Scholar] [CrossRef]
  85. Xu, H.J.; Ma, C.Q.; Ye, K. Early cretaceous granitoids and their implications for the collapse of the Dabie orogen, eastern China: SHRIMP zircon U-Pb dating and geochemistry. Chem. Geol. 2007, 240, 238–259. [Google Scholar] [CrossRef]
  86. Defant, M.J.; Drummond, M.S. Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature 1990, 347, 662–665. [Google Scholar] [CrossRef]
  87. Richards, J.P.; Kerrich, R. Adakie-like rocks: Their diverse origins and questionable role in metallogenesis. Econ. Geol. 2007, 102, 537–576. [Google Scholar] [CrossRef]
  88. Xu, W.C.; Zhang, H.F.; Luo, B.J.; Guo, L.; Yang, H. Adakite-like geochemical signature produced by amphibole–dominated fractionation of arc magmas: An example from the Late Cretaceous magmatism in Gangdese belt, southern Tibet. Lithos 2015, 232, 197–210. [Google Scholar] [CrossRef]
  89. Chen, B.; Jahn, B.M.; Arakawa, Y.; Zhai, M.G. Petrogenesis of the Mesozoic intrusive complexes from the southern Taihang Orogen, North China Craton: Elemental and Sr-Nd-Pb isotopic constraints. Contrib. Mineral. Petrol. 2004, 148, 489–501. [Google Scholar] [CrossRef]
  90. Guo, F.; Nakamuru, E.; Fan, W.M.; Kobayoshi, K.; Li, C.W. Generation of Palaeocene adakitic andesites by magma mixing; Yanji Area, NE China. J. Petrol. 2007, 48, 661–692. [Google Scholar] [CrossRef]
  91. Jiang, N.; Liu, Y.S.; Zhou, W.G.; Yang, J.H.; Zhang, S.Q. Derivation of Mesozoic adakitic magmas from ancient lower crust in the North China craton. Geochim. Cosmochim. Acta 2007, 71, 2591–2608. [Google Scholar] [CrossRef]
  92. Meng, Y.K.; Xiong, F.H.; Xu, Z.Q.; Ma, X.X. Petrogenesis of Late Cretaceous mafic enclaves and their host granites in the Nyemo region of southern Tibet: Implications for the tectonic-magmatic evolution of the Central Gangdese Belt. J. Asian Earth Sci. 2019, 176, 27–41. [Google Scholar] [CrossRef]
  93. He, Y.S.; Li, S.G.; Hoefs, J.; Huang, F.; Liu, S.A.; Hou, Z.H. Post-collisional granitoids from the Dabie orogen: New evidence for partial melting of thickened continental crust. Geochim. Cosmochim. Acta 2011, 75, 3815–3838. [Google Scholar] [CrossRef]
  94. Defant, M.; Xu, J.; Kepezhinskas, P.; Wang, Q.; Zhang, Q.; Xiao, L.; Defant, M. Adakites: Some variations on a theme. Acta Petrol. Sin. 2002, 18, 129–142. [Google Scholar]
  95. Rudnick, R. Growing from below. Nature 1990, 347, 711–712. [Google Scholar] [CrossRef]
  96. Maas, R.; Nicholls, I.A.; Legg, C. Igneous and metamorphic enclaves in the S-type Deddick granodiorite, Lachlan Fold Belt, SE Australia: Petrographic, geochemical and Nd-Sr isotopic evidence for crustal melting and magma mixing. J. Petrol. 1997, 38, 815–841. [Google Scholar] [CrossRef]
  97. White, R.V.; Tarney, J.; Kerr, A.C.; Saunders, A.D.; Kempton, P.D.; Pringle, M.S.; Klaver, G.T. Modification of an oceanic plateau, Aruba. Dutch Caribbean: Implications for the generation of continental crust. Lithos 1999, 46, 43–68. [Google Scholar] [CrossRef]
  98. Kazemi, K.; Kananian, A.; Xiao, Y.L.; Sarjoughian, F. Petrogenesis of Middle-Eocene granitoids and their mafic microgranular enclaves in central Urmia-Dokhtar Magmatic Arc (Iran): Evidence for interaction between felsic and mafic magmas. Geosci. Front. 2018, 10, 705–723. [Google Scholar] [CrossRef]
  99. Shellnutt, J.G.; Jahn, B.M.; Dostal, J. Elemental and Sr-Nd isotope geochemistry of microgranular enclaves from peralkaline A-type granitic plutons of the Emeishan large igneous province, SW China. Lithos 2010, 119, 34–46. [Google Scholar] [CrossRef]
  100. Barbarin, B. Mafic magmatic enclaves and mafic rocks associated with some granitoids of the central Sierra Nevada batholith, California: Nature, origin, and relations with the hosts. Lithos 2005, 88, 155–177. [Google Scholar] [CrossRef]
  101. Charles, N.; Gumiaux, C.; Augier, R.; Chen, Y.; Zhu, R.; Lin, W. Metamorphic core complexes vs. synkinematic plutons in continental extension setting: Insights from key structures (Shandong Province, eastern China). J. Asian Earth Sci. 2011, 40, 261–278. [Google Scholar] [CrossRef] [Green Version]
  102. Charles, N.; Augier, R.; Gumiaux, C.; Monié, P.; Chen, Y.; Faure, M.; Zhu, R. Timing, duration and role of magmatism in wide rift systems: Insights from the Jiaodong Peninsula (China, East Asia). Gondwana Res. 2013, 24, 412–428. [Google Scholar] [CrossRef] [Green Version]
  103. Wang, J.; Chang, S.C.; Wang, K.L.; Lu, H.B.; Zhang, H.C. Geochronology and geochemistry of early cretaceous igneous units from the central Sulu orogenic belt: Evidence for crustal delamination during a shift in the regional tectonic regime. J. Asian Earth Sci. 2015, 112, 49–59. [Google Scholar] [CrossRef]
  104. Zhu, R.X.; Xu, Y.G. The subduction of the west Pacific plate and the destruction of the north China craton. Sci. China Earth Sci. 2019, 62, 1340–1350. [Google Scholar] [CrossRef]
Figure 1. (a) Simplified tectonic framework of the Dabie–Sulu orogenic belt (from Yang [1,2] and Xu [14]). (b) Simplified geological map of the Sulu orogenic belt. (c) Modification of the Google Earth map. (d) Geological map of the study and adjacent areas (modified from 1:250,000 geological map of the PRC). WQYF: Wulian–Qingdao-Yantai fault; JXF: Jiashan–Xiangshui fault; MF: Milashan Fault; CCSD: China continental scientific drilling.
Figure 1. (a) Simplified tectonic framework of the Dabie–Sulu orogenic belt (from Yang [1,2] and Xu [14]). (b) Simplified geological map of the Sulu orogenic belt. (c) Modification of the Google Earth map. (d) Geological map of the study and adjacent areas (modified from 1:250,000 geological map of the PRC). WQYF: Wulian–Qingdao-Yantai fault; JXF: Jiashan–Xiangshui fault; MF: Milashan Fault; CCSD: China continental scientific drilling.
Minerals 10 00432 g001
Figure 2. (a,b) Field artificial exposures of the identified granitic pluton. (c) Granitic rocks with the mafic micro-granular enclave. (d,e) Micro-structural textures of the identified granitic pluton. (f) The boundary of mafic micro-granular enclaves (MMEs) and host granite pluton. (g) Plagioclase metacryst of MME. (h) Acicular apatite crystal. (i) Quartz metacryst of MMEs.
Figure 2. (a,b) Field artificial exposures of the identified granitic pluton. (c) Granitic rocks with the mafic micro-granular enclave. (d,e) Micro-structural textures of the identified granitic pluton. (f) The boundary of mafic micro-granular enclaves (MMEs) and host granite pluton. (g) Plagioclase metacryst of MME. (h) Acicular apatite crystal. (i) Quartz metacryst of MMEs.
Minerals 10 00432 g002
Figure 3. Representative zircon CL images of the granitic pluton (red circle for the U-Pb dating field, yellow circle for Lu–Hf field, white dotted circle for the oxygen isotopic field).
Figure 3. Representative zircon CL images of the granitic pluton (red circle for the U-Pb dating field, yellow circle for Lu–Hf field, white dotted circle for the oxygen isotopic field).
Minerals 10 00432 g003
Figure 4. Zircon U-Pb concordia diagrams for the Caochang granitic pluton in the Sulu belt ((ad) for the LA-ICP-MS technique, (e,f) for the SHRIMP II technique).
Figure 4. Zircon U-Pb concordia diagrams for the Caochang granitic pluton in the Sulu belt ((ad) for the LA-ICP-MS technique, (e,f) for the SHRIMP II technique).
Minerals 10 00432 g004
Figure 5. Zircon Hf-O isotopic results. (a) 176Lu/177Hf versus 176Hf/177Hf diagram. (b) Zircon U-Pb age versus εHf(t) diagram. (c) Histogram of εHf(t) values. (d) Histogram of the two-stage model age (TDM2). (e) Histogram of zircon oxygen isotopes. (f) Zircon δ18O values versus εHf(t) values. CC: continental crust; CHUR: chondritic uniform reservoir; red dot for zircon Hf isotopic value, one red dot representing one datum; red diamond for zircon oxygen value, one diamond for one oxygen datum
Figure 5. Zircon Hf-O isotopic results. (a) 176Lu/177Hf versus 176Hf/177Hf diagram. (b) Zircon U-Pb age versus εHf(t) diagram. (c) Histogram of εHf(t) values. (d) Histogram of the two-stage model age (TDM2). (e) Histogram of zircon oxygen isotopes. (f) Zircon δ18O values versus εHf(t) values. CC: continental crust; CHUR: chondritic uniform reservoir; red dot for zircon Hf isotopic value, one red dot representing one datum; red diamond for zircon oxygen value, one diamond for one oxygen datum
Minerals 10 00432 g005
Figure 6. (a) TAS diagram (from Middlemost [38]) (red dot for the studied samples). The dotted line shows the Irvine line (from Irvine and Baragar [39]), while the area below the dotted line is sub-alkaline, and that above the dotted line is alkaline. (b) Alkalinity ratio (A.R.) versus SiO2 plot (after Wright [40]). (c) SiO2 versus K2O plot (from Peccerillo [41]). (d) A/CNK versus A/NK plot (from Maniar and Piccoli [42]).
Figure 6. (a) TAS diagram (from Middlemost [38]) (red dot for the studied samples). The dotted line shows the Irvine line (from Irvine and Baragar [39]), while the area below the dotted line is sub-alkaline, and that above the dotted line is alkaline. (b) Alkalinity ratio (A.R.) versus SiO2 plot (after Wright [40]). (c) SiO2 versus K2O plot (from Peccerillo [41]). (d) A/CNK versus A/NK plot (from Maniar and Piccoli [42]).
Minerals 10 00432 g006
Figure 7. (a) Chondrite-normalized REE patterns (garnet and hornblende equilibrium melt modeling curves are from Moyen [43]; highly-fractionated magma curve are from Sun [44]) (red dot for the studied samples). (b) Primitive mantle-normalized multiple trace element diagram (chondrite normalization values are from Boynton [45]; primitive mantle normalization values are from Sun [46]).
Figure 7. (a) Chondrite-normalized REE patterns (garnet and hornblende equilibrium melt modeling curves are from Moyen [43]; highly-fractionated magma curve are from Sun [44]) (red dot for the studied samples). (b) Primitive mantle-normalized multiple trace element diagram (chondrite normalization values are from Boynton [45]; primitive mantle normalization values are from Sun [46]).
Minerals 10 00432 g007
Figure 8. (ac) Discrimination diagrams of A-type granite (from Whalen [53]) (red dot for the studied samples). (d,e) Whole-rock Zr and zircon Ti thermometers. (f) SiO2 versus P2O5 diagram (from Chappell [54]). (g,h) Rb versus Y and Rb versus Th diagrams (from Chappell [55]).
Figure 8. (ac) Discrimination diagrams of A-type granite (from Whalen [53]) (red dot for the studied samples). (d,e) Whole-rock Zr and zircon Ti thermometers. (f) SiO2 versus P2O5 diagram (from Chappell [54]). (g,h) Rb versus Y and Rb versus Th diagrams (from Chappell [55]).
Minerals 10 00432 g008
Figure 9. Compositional discrimination diagrams of granites (red dot for the studied samples). (a) Na2O + K2O + FeOT + MgO + TiO2 versus molar (Na2O + K2O)/(FeOT + MgO + TiO2) (from Patiño Douce [72]). (b) Molar Al2O3/(MgO + Fe2O3T) versus molar CaO/(MgO + Fe2O3T) (from Altherr [73]). (c) SiO2 versus Mg# diagram (Mg# = molecular Mg2+/(Mg2+ + Fe2+)) (from Rapp and Watson [71]). (d) Nb versus Nb/Th (primitive mantle data are after Hofmann [74], MORB + OIB and arc-type rocks fields are after Schmidberger and Hegner [75])
Figure 9. Compositional discrimination diagrams of granites (red dot for the studied samples). (a) Na2O + K2O + FeOT + MgO + TiO2 versus molar (Na2O + K2O)/(FeOT + MgO + TiO2) (from Patiño Douce [72]). (b) Molar Al2O3/(MgO + Fe2O3T) versus molar CaO/(MgO + Fe2O3T) (from Altherr [73]). (c) SiO2 versus Mg# diagram (Mg# = molecular Mg2+/(Mg2+ + Fe2+)) (from Rapp and Watson [71]). (d) Nb versus Nb/Th (primitive mantle data are after Hofmann [74], MORB + OIB and arc-type rocks fields are after Schmidberger and Hegner [75])
Minerals 10 00432 g009
Figure 10. (af) Major elements versus SiO2 Harker diagrams (red dot for the studied samples). (g,h) Diagrams of SiO2/Al2O3 versus CaO/Al2O3 and Al2O3/MgO versus CaO/MgO (the magma mixing model is from Langmuir [78]).
Figure 10. (af) Major elements versus SiO2 Harker diagrams (red dot for the studied samples). (g,h) Diagrams of SiO2/Al2O3 versus CaO/Al2O3 and Al2O3/MgO versus CaO/MgO (the magma mixing model is from Langmuir [78]).
Minerals 10 00432 g010
Figure 11. Mineral fractional evolution curves for the studied samples; (a,b) Sr versus Rb/Sr and Sr versus Ba diagrams (modeling curves are from Yang [2]). (c) La versus (La/Yb)N; (d) Th versus LREE. Allan: allanite; Ap: apatite; Bt: biotite; Cpx: clinopyroxene; Hb: hornblende; Kfs or Kf: K-feldspar; Mon: monazite; Opx: orthopyroxene; Plag or Pl: plagioclase; Sph: sphene; Zr: zircon
Figure 11. Mineral fractional evolution curves for the studied samples; (a,b) Sr versus Rb/Sr and Sr versus Ba diagrams (modeling curves are from Yang [2]). (c) La versus (La/Yb)N; (d) Th versus LREE. Allan: allanite; Ap: apatite; Bt: biotite; Cpx: clinopyroxene; Hb: hornblende; Kfs or Kf: K-feldspar; Mon: monazite; Opx: orthopyroxene; Plag or Pl: plagioclase; Sph: sphene; Zr: zircon
Minerals 10 00432 g011
Figure 12. Discrimination diagrams for adakitic rocks (red dot for studied samples); (a,b) Y versus Sr/Y and YbN versus (La/Yb)N diagrams for adakite rocks (from Petford [83] and Defant [94]). (c) SiO2 versus Sr. (d) CaO versus Sr. (e) SiO2 versus (Dy/Yb)N, (f) SiO2 versus Sr/CaO (cf are from He [93]).
Figure 12. Discrimination diagrams for adakitic rocks (red dot for studied samples); (a,b) Y versus Sr/Y and YbN versus (La/Yb)N diagrams for adakite rocks (from Petford [83] and Defant [94]). (c) SiO2 versus Sr. (d) CaO versus Sr. (e) SiO2 versus (Dy/Yb)N, (f) SiO2 versus Sr/CaO (cf are from He [93]).
Minerals 10 00432 g012
Table 1. Zircon U-Pb dating results of the granitic pluton in the study area.
Table 1. Zircon U-Pb dating results of the granitic pluton in the study area.
Spot NumberThUTh/UTi contents°C207Pb/235U206Pb/238U206Pb/238UMethods
RatioRatioAge (Ma)
My1-2-1-011691421.194.686790.128100.012570.018600.00045118.82.8I
My1-2-1-021681181.427.707220.135580.013690.018550.00034118.52.1I
My1-2-1-042291731.326.767100.117450.009550.018740.00032119.72.1I
My1-2-1-054052002.024.126690.128900.011530.019750.00030126.11.9I
My1-2-1-062523750.677.447190.448940.015700.054280.00058340.83.5I
My1-2-1-081891681.123.556570.119940.009220.018320.00030117.01.9I
My1-2-1-091691151.472.916420.144380.013950.019780.00051126.23.2I
My1-2-1-101761371.287.317170.144200.014350.019430.00041124.02.6I
My1-2-1-111751331.323.366530.116030.012810.018440.00033117.82.1I
My1-2-1-125122472.075.876980.131400.007870.018410.00027117.61.7I
My1-2-1-132881951.473.796620.108970.013150.017830.00034113.92.1I
My1-2-1-141591371.174.136690.142100.010810.019140.00034122.22.1I
My1-2-1-151981271.5638.988920.126860.011490.018630.00034119.02.2I
My1-2-1-162131181.805.056850.118500.011700.017890.00036114.32.3I
My1-2-1-172161501.443.426540.124530.009840.019190.00033122.52.1I
My1-2-1-183001871.603.746610.125570.009210.018610.00028118.91.8I
My1-2-1-192121751.225.056850.123920.010020.018960.00029121.11.8I
My1-2-2-012101421.485.786970.127840.009260.019070.00035121.82.2I
My1-2-2-021521281.186.217030.135030.010140.018250.00034116.62.2I
My1-2-2-031361171.165.106860.128190.013010.018600.00033118.82.1I
My1-2-2-041841431.295.266890.135670.011090.018340.00035117.22.2I
My1-2-2-053561801.972.536310.120470.009560.018780.00031120.01.9I
My1-2-2-061841501.236.407050.125270.010130.018700.00032119.42.0I
My1-2-2-072871661.739.757430.138830.012320.018750.00033119.72.1I
My1-2-2-0876661.162.826390.125120.011450.018560.00044118.62.8I
My1-2-2-093852171.783.756610.132040.008920.018420.00026117.71.6I
My1-2-2-101161101.055.466920.118260.010860.018870.00039120.52.5I
My1-2-2-112671841.455.156870.122330.009260.018540.00029118.41.9I
My1-2-2-128733092.824.236710.130810.007750.018610.00023118.91.5I
My1-2-2-136253032.065.726960.121590.006680.018560.00022118.51.4I
My1-2-2-151001020.995.506930.140100.011420.019880.00039126.92.4I
My1-2-2-161481141.309.527410.145990.014920.019660.00041125.52.6I
My1-2-2-185352252.384.796810.133950.008590.019440.00032124.12.0I
My1-2-2-195672632.156.367050.136670.010360.018670.00033119.22.1I
My1-2-2-204632312.003.926650.117820.007200.018130.00026115.81.6I
My1-2-2-211331201.123.546570.114090.010770.017890.00037114.32.3I
My1-2-2-221781401.276.337050.124650.013780.018410.00035117.62.2I
My1-2-2-232313560.652.766380.119780.006290.017890.00023114.31.5I
My1-2-2-262581811.435.046850.121050.008480.018630.00032119.02.0I
MY1211-1.11381331.04//0.163800.000560.020210.00018129.02.3II
MY1211-1.22211981.11//0.122000.00190.018960.00019121.12.3II
MY1211-1.32001821.10//0.152000.00160.02020.00019128.92.5II
MY1211-1.41331510.88//0.126000.00180.018510.00028118.23.3II
MY1211-1.53372311.46//0.120000.00150.019440.00018124.12.3II
MY1211-1.61521530.99//0.098000.0030.018750.00036119.84.3II
MY1211-1.74272751.55//0.103000.00180.018630.00021119.02.5II
MY1211-1.81691501.13//0.129000.00290.019050.00025121.63.0II
MY1211-1.92472101.18//0.131000.00110.019310.00019123.32.3II
MY1211-1.101721621.06//0.105000.00310.019280.00029123.13.5II
MY1211-1.112582111.23//0.079000.00180.018840.00017120.32.1II
MY1211-1.124022381.69//0.134000.00150.02010.00018128.32.3II
MY1211-1.132191891.16//0.152000.00130.020190.00019128.92.4II
MY1211-1.141391550.90//0.119000.00180.019180.0002122.52.4II
MY1211-1.152281701.34//0.087000.00320.019150.00029122.33.5II
MY1211-1.161361520.90//0.144000.00270.019660.00026125.53.3II
MY1211-1.174202711.55//0.110000.00210.018230.00023116.42.6II
Notes: common Pb corrections are from Andersen [36]. I: method for LA-ICP-MS; II: method for the SHRIMP II technique. ’’/’’ for null for SHRIMP II method.
Table 2. Zircon Lu–Hf isotopic results of the granitic pluton in the study area.
Table 2. Zircon Lu–Hf isotopic results of the granitic pluton in the study area.
Spot NumberAge176Hf/177Hf176Lu/177Hf176Yb/177HfIHfεHf(t)fLu/HfTDM1TDM2
My121-011200.281840 0.000009 0.001198 0.000007 0.027346 0.000178 0.28184 −30.43 0.32 −0.96 1992 12 3084 19
My121-021200.281804 0.000009 0.001066 0.000004 0.024296 0.000079 0.28180 −32.13 0.31 −0.97 2036 12 3163 19
My121-041200.281848 0.000008 0.001230 0.000005 0.028240 0.000122 0.28185 −30.57 0.30 −0.96 1983 12 3066 18
My121-051200.281869 0.000009 0.000914 0.000014 0.021804 0.000406 0.28187 −29.80 0.32 −0.97 1938 13 3019 20
My121-063410.281801 0.000008 0.001190 0.000039 0.031458 0.001135 0.28179 −27.51 0.30 −0.96 2047 12 3042 18
My121-081200.281859 0.000009 0.000965 0.000006 0.021987 0.000117 0.28186 −29.74 0.33 −0.97 1954 13 3042 20
My121-091200.281875 0.000009 0.000785 0.000008 0.018712 0.000174 0.28187 −29.59 0.33 −0.98 1923 13 3006 20
My121-101200.281845 0.000009 0.000973 0.000013 0.022548 0.000282 0.28184 −30.66 0.32 −0.97 1974 12 3072 20
My121-111200.281847 0.000009 0.001342 0.000012 0.032351 0.000303 0.28184 −30.60 0.33 −0.96 1990 13 3068 20
My121-121200.281860 0.000010 0.001369 0.000007 0.032400 0.000172 0.28186 −30.16 0.37 −0.96 1974 15 3041 23
My121-131200.281800 0.000010 0.001151 0.000005 0.026790 0.000123 0.28180 -31.82 0.35 −0.97 2045 14 3171 22
My121-141200.281831 0.000009 0.001071 0.000001 0.024793 0.000044 0.28183 −31.15 0.32 −0.97 1998 12 3103 19
My121-151200.281850 0.000008 0.001030 0.000006 0.024307 0.000174 0.28185 −30.49 0.29 −0.97 1970 11 3062 18
My121-161200.281800 0.000011 0.002338 0.000043 0.056820 0.001104 0.28180 −32.35 0.38 −0.93 2111 15 3175 23
My121-171200.281847 0.000010 0.001233 0.000020 0.028555 0.000421 0.28184 −30.61 0.37 −0.96 1984 15 3069 23
My121-181200.281858 0.000009 0.000813 0.000002 0.019332 0.000029 0.28186 −29.78 0.31 −0.98 1948 12 3044 19
My121-191200.281862 0.000009 0.001128 0.000006 0.027226 0.000113 0.28186 −30.06 0.31 −0.97 1958 12 3035 19
My122-011190.281832 0.000009 0.001638 0.000021 0.039715 0.000564 0.28183 −31.18 0.33 −0.95 2027 13 3103 20
My122-021190.281860 0.000010 0.000915 0.000004 0.021709 0.000130 0.28186 −30.16 0.35 −0.97 1951 14 3040 21
My122-031190.281851 0.000009 0.001102 0.000009 0.025715 0.000200 0.28185 −30.49 0.32 −0.97 1973 13 3061 20
My122-041190.281847 0.000008 0.001058 0.000006 0.025040 0.000122 0.28184 −30.19 0.29 −0.97 1975 12 3069 18
My122-051190.281854 0.000009 0.001105 0.000018 0.026502 0.000496 0.28185 −30.38 0.32 −0.97 1968 13 3054 20
My122-061190.281884 0.000008 0.001285 0.000006 0.030369 0.000140 0.28188 −29.34 0.29 −0.96 1937 11 2989 18
My122-071190.281848 0.000010 0.003147 0.000051 0.078701 0.001352 0.28184 −30.74 0.37 −0.91 2088 16 3074 23
My122-081190.281880 0.000011 0.001113 0.000067 0.026345 0.001534 0.28188 −29.46 0.38 −0.97 1933 15 2997 23
My122-091190.281871 0.000009 0.001041 0.000015 0.024160 0.000300 0.28187 −29.34 0.31 −0.97 1941 12 3016 19
My122-101190.281865 0.000008 0.000753 0.000009 0.017240 0.000179 0.28186 −29.97 0.30 −0.98 1936 12 3029 18
My122-111190.281862 0.000009 0.001206 0.000021 0.028586 0.000564 0.28186 −30.11 0.32 −0.96 1963 13 3037 20
My122-121190.281871 0.000009 0.001453 0.000071 0.037634 0.001977 0.28187 −29.82 0.34 −0.96 1963 14 3019 21
My122-131190.281861 0.000008 0.001040 0.000005 0.024390 0.000134 0.28186 −30.13 0.29 −0.97 1956 12 3039 18
My122-151190.281874 0.000009 0.001114 0.000004 0.026054 0.000099 0.28187 −29.24 0.30 −0.97 1941 12 3010 19
My122-161190.281887 0.000009 0.000621 0.000024 0.014994 0.000592 0.28189 −29.16 0.30 −0.98 1898 12 2979 18
My122-181190.281875 0.000008 0.001066 0.000013 0.024323 0.000330 0.28187 −29.62 0.29 −0.97 1937 11 3007 18
My122-191190.281847 0.000012 0.001792 0.000014 0.043603 0.000350 0.28184 −30.67 0.43 −0.95 2014 17 3071 27
My122-201190.281858 0.000009 0.001045 0.000006 0.024099 0.000134 0.28186 −30.23 0.31 −0.97 1960 12 3044 19
My122-211190.281876 0.000009 0.000720 0.000007 0.016596 0.000204 0.28187 −29.14 0.32 −0.98 1919 12 3004 20
My122-221190.281866 0.000008 0.001126 0.000014 0.026669 0.000330 0.28186 −29.96 0.27 −0.97 1953 11 3028 17
My122-231190.281809 0.000010 0.002324 0.000056 0.057767 0.001431 0.28180 −32.06 0.37 −0.93 2098 15 3157 23
My122-261190.281885 0.000009 0.000911 0.000012 0.020823 0.000306 0.28188 −29.25 0.31 −0.97 1915 12 2984 19
Table 3. Zircon oxygen isotopic results of the granitic pluton in the study area.
Table 3. Zircon oxygen isotopic results of the granitic pluton in the study area.
Sample NumberZircon δ18O (‰)
MY1211-15.530.20
MY1211-106.000.21
MY1211-115.500.16
MY1211-126.240.17
MY1211-135.850.15
MY1211-146.010.24
MY1211-155.160.26
MY1211-165.680.20
MY1211-175.990.18
MY1211-185.930.16
MY1211-195.680.14
MY1211-25.680.18
MY1211-205.660.14
MY1211-215.790.24
MY1211-225.440.18
MY1211-235.500.29
MY1211-245.270.21
MY1211-255.920.16
MY1211-265.940.20
MY1211-275.760.14
MY1211-285.970.21
MY1211-295.120.18
MY1211-35.790.20
MY1211-305.830.23
MY1211-45.720.12
MY1211-55.780.22
MY1211-65.700.16
MY1211-76.040.20
MY1211-85.430.18
MY1211-95.620.23
Table 4. Whole-rock geochemical results of the granitic pluton in the study area.
Table 4. Whole-rock geochemical results of the granitic pluton in the study area.
Sample NumberMy1-2-1My1-2-2My1-2-3My1-2-4My1-2-4-1My1-2-5My1-2-6My1-2-7My1-2-8
SiO264.3765.4365.1364.9564.2260.6365.4665.4558.80
TiO20.570.530.520.550.540.610.550.550.55
Al2O315.5215.7015.5816.0215.7516.0315.7615.6616.51
TFe2O34.483.974.124.214.165.074.204.254.59
Na2O3.633.673.743.793.783.983.873.606.80
K2O3.563.253.653.313.643.293.373.581.49
MgO2.111.761.762.001.812.602.022.021.66
CaO3.393.312.813.732.983.733.273.524.47
MnO0.110.120.150.070.140.180.120.070.17
P2O50.210.190.190.200.210.230.210.210.20
LOI1.622.162.750.572.63.471.510.894.51
Total99.57100.09100.4099.4099.8399.82100.3499.8099.75
Al838848487584467858168575888080844238381891774
Li13.0015.8019.5010.4015.7025.7016.9010.205.07
Be1.602.092.751.542.432.661.821.442.91
Sc16.8014.1015.6014.2017.9020.6016.6017.3013.50
Ti354731953207340833783792340634303504
V74.7067.1067.8069.2067.7087.3072.4070.6061.50
Cr24.1020.8020.4023.3020.8060.8023.0021.7019.60
Mn9219721194598108214239676131291
Co11.7010.4010.5011.3010.8014.4011.0011.2011.20
Ni13.5011.7012.1013.1012.7027.2013.3012.7013.00
Cu8.206.966.879.236.959.626.649.414.66
Zn56.2052.9079.4051.8067.1072.6064.1050.7067.60
Ga17.9017.2018.0018.4018.0019.2018.2018.1019.90
Rb77.9064.2097.5062.5081.1078.2076.6067.0047.20
Sr547469423580495518581565317
Y16.5014.2015.0015.1014.6018.0015.3015.7013.20
Zr205201197198195197198208202
Nb10.109.849.749.419.2411.708.859.6110.40
Mo0.520.701.860.350.441.051.110.410.26
Cd<0.0300.030.14<0.030<0.030<0.0300.040.03<0.030
Sn1.181.081.221.131.341.481.271.136.03
Cs1.371.532.290.981.241.451.191.050.70
Ba17411518139515401736194415831614640
Hf5.085.104.944.764.665.004.775.134.91
Ta0.740.720.720.670.650.750.640.710.77
W0.220.300.770.170.770.400.270.150.87
Tl0.440.380.610.450.620.480.470.360.28
Pb26.9027.0050.0022.6026.0032.5026.3023.1013.10
Bi0.100.120.130.040.200.130.070.040.51
Th8.5810.0011.7010.108.088.629.189.658.67
U1.831.862.032.141.923.071.592.013.38
La38.4036.4048.7044.3042.7040.0048.3049.7042.40
Ce69.2056.9079.4073.9069.6069.4078.1079.9072.30
Pr8.247.338.608.368.248.488.729.107.90
Nd29.7025.5029.0028.9028.3030.4029.4031.0026.80
Sm4.744.114.354.474.475.044.474.694.12
Eu1.211.091.091.141.191.231.171.180.99
Gd4.743.994.434.544.364.884.554.614.01
Tb0.600.520.540.550.540.620.550.580.50
Dy3.022.582.722.722.633.142.712.852.39
Ho0.650.540.580.570.540.670.570.580.49
Er1.711.551.641.591.541.901.661.681.41
Tm0.280.251.060.250.240.300.250.260.22
Yb1.751.571.601.571.521.881.591.641.34
Lu0.290.260.270.260.250.310.260.270.23
10000Ga/Al2.132.032.132.142.102.182.162.162.17
FeOT/MgO1.912.032.111.892.071.761.871.892.49
Zr+Nb+Ce+Y300.8281.9301.1296.4288.4296.1300.3313.2297.9
A/CNK0.971.011.020.961.010.950.990.970.79
A/NK1.581.641.541.631.551.581.571.601.29
Nb/Th1.180.980.830.931.141.360.961.001.20
(La/Yb)N14.7915.6320.5219.0218.9414.3420.4820.4321.33
Mg#494746494751494942
(K2O+Na2O)/CaO2.122.092.631.902.491.952.212.041.85
K2O+Na2O7.196.927.397.107.427.277.247.188.29
Sr/CaO16114215115516613917816171
(Dy/Yb)N1.121.071.101.121.121.081.111.131.16
Zr/Hf40.3539.4139.8841.6041.8539.4041.5140.5541.14
Rb/Sr0.140.140.230.110.160.150.130.120.15
Sr/Y33.1533.0328.2038.4133.9028.7837.9735.9924.02
SiO2/Al2O34.154.174.184.054.083.784.154.183.56
Al2O3/MgO7.368.928.858.018.706.177.807.759.95
CaO/Al2O30.220.210.180.230.190.230.210.220.27
A.R.1.381.331.401.331.391.321.351.371.11
Al2O3/(MgO+FeOT)1.401.641.601.531.591.221.501.481.64
CaO/(MgO+FeOT)0.560.630.530.650.550.520.570.610.81
Molar (N+K)/(F+M+T)0.830.940.970.880.960.730.890.871.19
N+K+F+M+T14.1913.0513.7013.6013.9015.6113.7513.7215.36
YbN8.377.517.667.517.279.007.617.856.41
LREE/HREE11.6211.6613.3313.3713.3011.2814.0214.0814.59
∑REE164.53142.59183.98173.12166.12168.25182.30188.04165.10
LREE156.23135.32175.57165.61158.86159.43174.71180.18158.52
Eu/Eu*0.780.820.760.770.820.760.790.780.74
CIPW calculated results (%)
Ab31.0031.1631.6632.4132.1933.9232.7930.6757.96
Or21.2319.2721.5819.7721.6519.5819.9421.308.87
An15.6815.2312.7017.1013.5016.2715.7016.1010.20
Plagioclase An number323227332831313314
Abbreviation notes. A.R. = (Al2O3 + CaO + Na2O + K2O)/(Al2O3 + CaO − Na2O − K2O); (N + K)/(F + M + T) = (Na2O + K2O)/(FeOT + MgO + TiO2); N + K + F + M + T = Na2O + K2O + FeOT + MgO + TiO2. ∑REE: total rare earth element; LREE: light REE; HREE: heavy REE; LOI: loss on ignition. Eu/Eu* representing Eu element anomaly; An: anorthite.

Share and Cite

MDPI and ACS Style

Meng, Y.; Wang, Z.; Gan, B.; Liu, J. Petrogenesis and Tectonic Implications of the Early Cretaceous Granitic Pluton in the Sulu Orogenic Belt: The Caochang Granitic Pluton as an Example. Minerals 2020, 10, 432. https://doi.org/10.3390/min10050432

AMA Style

Meng Y, Wang Z, Gan B, Liu J. Petrogenesis and Tectonic Implications of the Early Cretaceous Granitic Pluton in the Sulu Orogenic Belt: The Caochang Granitic Pluton as an Example. Minerals. 2020; 10(5):432. https://doi.org/10.3390/min10050432

Chicago/Turabian Style

Meng, Yuanku, Zhongbo Wang, Baoping Gan, and Jinqing Liu. 2020. "Petrogenesis and Tectonic Implications of the Early Cretaceous Granitic Pluton in the Sulu Orogenic Belt: The Caochang Granitic Pluton as an Example" Minerals 10, no. 5: 432. https://doi.org/10.3390/min10050432

APA Style

Meng, Y., Wang, Z., Gan, B., & Liu, J. (2020). Petrogenesis and Tectonic Implications of the Early Cretaceous Granitic Pluton in the Sulu Orogenic Belt: The Caochang Granitic Pluton as an Example. Minerals, 10(5), 432. https://doi.org/10.3390/min10050432

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop