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

: The subduction processes and geodynamic scenarios of the late-stage southward subduction of the Mongol–Okhotsk oceanic slab since the Early Jurassic are subjects of great debate. This contribution presents new U–Pb zircon dating, trace element geochemistry, Ti-in zircon geothermometry, and Lu–Hf isotopes of zircon, as well as bulk-rock geochemical data for Early–Middle Jurassic intrusive rocks in the Erguna Block, NE China. Approximately 181–198 Ma monzogranites and ca. 162–174 Ma quartz monzonites were identiﬁed in the block. The Early Jurassic monzogranites are high-K calc-alkaline I-type granites, which display moderately concave-upward rare earth element (REE) patterns with slightly negative Eu anomalies, and low zircon crystallization temperatures. The Middle Jurassic quartz monzonites have low Yb and Y concentrations, high Sr / Y ratios, and strong high ﬁeld strength elements (HFSEs) depletions, that are in excellent agreement with adakitic rocks. They exhibit right-sloping REE patterns with negligible Eu anomalies, and a wide range of zircon crystallization temperatures. The intrusions yield ε Hf (t) values between − 4.1 to + 4.8 and juvenile two-stage model (T DM2 ) ages varying from 918–1488 Ma. The geochemical and isotopic signatures suggest that the monzogranites were likely derived by the partial melting of K-rich meta-basalts within the lower part of a juvenile crust that had medium-thickness ( ≤ 40 km), with the involvement of minor mantle materials. Whereas, the quartz monzonites were possibly produced by partial melting of a thickened continental lower crust ( ≥ 50 km). The Mongol–Okhotsk tectonic regime played a dominant role in accounting for their formation. An Andean-type continental arc setting was developed during the Early–Middle Jurassic, with gradual thickening of the continental crust. The signiﬁcant crustal thickening may reach its ultimate stage at ca. 162–174 Ma, which marks the tectonic transition from compression to extension. The southward subduction beneath the Erguna Block was continuous and stable during the Early Jurassic. Rollback of the subducted slab occurred at ca. 174–177 Ma, followed by moderate magmatic activities represented by adakitic rocks.


Introduction
The Mongol-Okhotsk orogenic belt stretches across the Russian Transbaikal and Northeastern (NE) Mongolia (Figure 1a). The belt is considered to be the youngest section of the Central Asian Orogenic Belt (CAOB) that was formed through a gigantic accretionary orogeny from the Neoproterozoic-Late Mesozoic [1][2][3][4]. The Mongol-Okhotsk orogenic belt is tectonically situated between the Siberian Craton to the northeast and the collages of Mongolia-NE China continental blocks to the southeast since the Late Paleozoic (Figure 1a,b) [1][2][3][4], and its geodynamic evolution played a dominant role in controlling the stratigraphy, magmatism, mineralization, as well as topography of the eastern part of the CAOB during the Mesozoic [5][6][7][8].
The Erguna Block of NE China is situated immediately to the southeast of the Mongol-Okhotsk orogenic belt (Figure 1a,b), and is characterized by extensive intrusive and volcanic activities during the Early and Late Mesozoic, respectively [7][8][9][10], providing an ideal natural laboratory for exploring the Mesozoic tectonic evolution of the belt in Chinese territory. A large number of geochronological and geochemical studies focusing on the Mesozoic igneous rocks have thus been performed in the past decade (references in Tables S2 and S6), and have made great progress in understanding the early-stage (from Late Permian to Triassic [7,11]) southward subduction of the Mongol-Okhotsk oceanic plate (MOOP). In contrast, the subduction processes and geodynamic scenarios of the late-stage southward subduction of the MOOP since the Early Jurassic remain pending problems.
Although Early Jurassic subduction-related igneous rocks are widespread in the Erguna Block and adjacent areas (Figure 1c) [11][12][13], the associated subduction processes remain debatable. Some studies suggest the rollback of a subducted slab with supra subduction extensional settings was involved [12]; while others argue that a stable arc-trench system accompanied by continuous subduction under compressional environments was dominant [11,13]. In contrast to the Early Jurassic igneous rocks, the Middle Jurassic igneous rocks are poorly exposed in the Erguna Block and adjacent regions ( Figure 1c). Most previous studies suggest that the Middle Jurassic marks the final closure of the Mongol-Okhotsk Ocean (MOO), and the poor exposure of igneous rocks reflects a magmatic quiescent period related to the continent-continent collision [7,11,14]. Other studies argue that some Late Jurassic magmatic rocks were generated by subduction-related processes, and thus the Middle Jurassic tectonic setting of the Erguna Block was controlled by the southward subduction of the MOOP [13,[15][16][17]. These controversies result in an incomplete understanding of the late-stage southward subduction of the MOOP, as well as the tectonic settings in relation to the subduction processes.
In this study, we have first identified two Middle Jurassic intrusions in the Erguna Block and found that the rocks have typical adakitic affinity that indicates an environment of thickened continental crust. Early Jurassic high-K calc-alkaline I-type granites were also identified in the Erguna Block. Here, we present new U-Pb dating, trace element geochemistry, Ti-in zircon geothermometry, and Lu-Hf isotopes of zircons, as well as bulk-rock geochemical data for these intrusive rocks. These data enable to constrain the petrogenesis of the intrusions and contribute to infer the associated tectonic settings. For the widespread Early Jurassic igneous rocks in the Erguna Block, a large geochemical, isotopic, and geochronological database was compiled using previously published data (Tables S2, S6 and S8). Comparisons of our data with the large dataset enable a more convincing interpretation of the results. The findings could contribute to constrain the Early-Middle Jurassic tectonic settings of the Erguna Block, and give a better understanding of the subduction processes and geodynamic scenarios of the late-stage southward subduction of the MOOP.  Table S1) and previous zircon U-Pb age data from references in Table S2.  Table S1) and previous zircon U-Pb age data from references in Table S2.

Zircon Morphology, Trace Element Compositions, and U-Pb Geochronology
After detailed petrographic observations, two monzogranite (14ES14 and 14ER495) and two quartz monzonite (18X19 and 18X20) samples from the Erguna Block were chosen for LA-ICP-MS zircon U-Pb dating and trace element analyses. The results are given in Tables S3 and S4, [36], the analyzed zircons display a linear trend ranging from magmatic fields toward hydrothermal fields. These features may imply that the zircons have undergone moderate primary high-T hydrothermal alteration during the deuteric stage of magma evolution, with a relatively deep intrusion depth.  [36]; chondrite values are from Sun and McDonough [37]). The abbreviation "t" in e-h represents 206 Pb/ 238 U ages yielded from zircons.
Only age data with concordance >87% have been considered in the following. For monzogranite sample 14ES14-1, three analyses with low concordance (81.2%, 58.5%, and 81.1%) were excluded, seven analyses of zircons yielded a weighted mean 206 Pb/ 238 U age of 181 ± 2 Ma (MSWD = 1.2, solid circles in Figure 5e), with the remaining 10 zircons yielding ages ranging from 185 ± 4 to 198 ± 4 Ma (dotted circles in Figure 5e). Among the 20 analyses from monzogranite sample 14ER495, one spot with low concordance (79.5%) was excluded, the remaining zircons recorded two resolvable growth events. Nine analyses yielded a weighted mean 206 Pb/ 238 U age of 177 ± 2 Ma (MSWD = 0.78, solid circles in Figure 5f), with a further nine zircons yielding a weighted mean 206 Pb/ 238 U age of 186 ± 2 Ma (MSWD = 0.70, dotted circles in Figure 5f). The concordia age of ca. 186 Ma (sample 14ER495), and the age range of ca. 185-198 Ma (sample 14ER495) possibly imply the pulse of plutonic activities during the Early Jurassic. The concordia ages of ca. 181 Ma (sample 14ES14-1) and 177 Ma (sample 14ER495) are within error of each other and probably represent the timing of final crystallization of the monzogranite intrusions in the study area.
A total of 11 analyses of zircons from quartz monzonite sample 18X19 yielded a weighted mean 206 Pb/ 238 U age of 174 ± 3 Ma (MSWD = 0.59, Figure 5g), with one rejected spot (zircon was penetrated during ablation). For quartz monzonite sample 18X20, nine zircons yielded a weighted mean 206 Pb/ 238 U age of 162 ± 3 Ma (MSWD = 1.01, solid circles in Figure 5h), with a further three zircons yielding 206 Pb/ 238 U age from 172 ± 5 to 179 ± 5 Ma (dotted circles in Figure 5h). The weighted mean 206 Pb/ 238 U age of ca. 174 Ma (sample 18X19) is indicative of a Middle Jurassic plutonic activity in the study area, and the age of 162 ± 3 Ma obtained from sample 18X20 suggests the final intrusion of the quartz monzonite. The dating results indicate that the quartz monzonites were emplaced during the Middle Jurassic, rather than the Late Triassic-Early Jurassic or Late Jurassic as previously considered by geological mapping [26].

Ti-in Zircon Geothermometry and Zircon Saturation Temperature
The zircon crystallization temperatures of the monzogranites and quartz monzonites were calculated employing the equation for the Ti-in zircon geothermometer established by Ferry and Watson [38] from high-temperature experimental analyses. The Ti concentration in zircon is sensitive to temperature change and could be employed as geothermometer following: where αSiO 2 and αTiO 2 are activity coefficients of SiO 2 and TiO 2 , respectively. For the analyzed monzogranites and quartz monzonites, αSiO 2 is assumed to be 1 since SiO 2 content has reached saturation evidenced by the existence of quartz, and αTiO 2 is assumed to be 0.6 based on the absence of Ti-bearing minerals (e.g., ilmenite and rutile) under optical observations. The calculation results are presented in Table S4. Zircons with Ti concentrations >100 ppm were rejected during the calculation.  [39,40]. The zircon saturation temperature (T Zr ) can be measured using the following equation: where D Zr Zircon/Melt is the ratio of Zr content in zircon to that in the melt, and M is defined as (2Ca + K + Na)/(Si × Al) cation and obtained using host-rock-normalized concentration. The activity coefficient is presumed to be 1 and Zr in zircon is considered to be 496000 × 10 −6 for pure zircons. As shown in Table S5, the obtained T Zr values exhibit a tight range. For the monzogranites and quartz monzonites, T Zr values vary between 771-790 • C and 782-806 • C, respectively. The zircon crystallization temperatures calculated from the Ti-in zircon geothermometer exhibit a wide range from 755-900 • C and 652-852 • C for the monzogranites and quartz monzonites, respectively. These temperatures represent the crystallization of single zircon crystals in the host silicic melt, and generally vary around the corresponding T Zr values calculated from whole-rocks ( Figure S1). Besides, quartz monzonite sample 18X19 has zircon crystallization temperatures much lower than its T Zr , suggesting that zircons did not crystalize immediately when Zr was saturated in the host melt, and the melt probably experienced a complex ascending process.

Whole-Rock Major and Trace Elements
Major and trace element compositions of the monzogranites and quartz monzonites are given in Table S5. All samples have low loss on ignition (LOI) values below 2%, suggesting that the weathering or fluid modification after crystallization can be neglected.
The  The zircon crystallization temperatures calculated from the Ti-in zircon geothermometer exhibit a wide range from 755-900 °C and 652-852 °C for the monzogranites and quartz monzonites, respectively. These temperatures represent the crystallization of single zircon crystals in the host silicic melt, and generally vary around the corresponding TZr values calculated from whole-rocks ( Figure S1). Besides, quartz monzonite sample 18X19 has zircon crystallization temperatures much lower than its TZr, suggesting that zircons did not crystalize immediately when Zr was saturated in the host melt, and the melt probably experienced a complex ascending process.

Whole-Rock Major and Trace Elements
Major and trace element compositions of the monzogranites and quartz monzonites are given in Table S5. All samples have low loss on ignition (LOI) values below 2%, suggesting that the weathering or fluid modification after crystallization can be neglected.
The Early Jurassic monzogranites are plotted on granodiorite and granite fields on a total alkali versus silica (TAS) variation diagram (Figure 6a [42]); (c) A/CNK versus A/NK (after Maniar and Piccoli [43]). The solid triangles indicate the data in this study and the + and × symbols indicate published data (same as below) from references in Table S6.  (Figure 7a). In primitive mantle-normalized trace element spidergrams, the samples exhibit enrichment of LREE and large ion lithophile elements (LILEs; e.g., Rb, Ba, and K), and depletion of HFSEs (e.g., Nb, Ta, and Ti) and P (Figure 7b). They are non-adakitic rocks as their geochemical signatures are distinct from typical adakites derived from both subducted oceanic plate and thickened or delaminated lower crust [44][45][46]. This is further  (Figure 7a). In primitive mantle-normalized trace element spidergrams, the samples exhibit enrichment of LREE and large ion lithophile elements (LILEs; e.g., Rb, Ba, and K), and depletion of HFSEs (e.g., Nb, Ta, and Ti) and P (Figure 7b). They are non-adakitic rocks as their geochemical signatures are distinct from typical adakites derived from both subducted oceanic plate and thickened or delaminated lower crust [44][45][46]. This is further supported by plots of (La/Yb) N versus Yb N , and Sr/Y versus Y (Figure 8a,b).    [44]; Data of magmatic rocks in the Andean orogenic arcs in d are from Asadi et al. [49]; chondrite values are from Sun and McDonough [37]). G518 is an Eastern Pontides gabbro employed as the source rock for the batchmelting modeling [50].

Zircon Lu-Hf Isotopes
A total of 21 representative zircons from four selected samples were chosen for in situ Lu-Hf isotopic analyses. The obtained data are presented in Table S7.
The 176 Hf/ 177 Hf ratios of the nine analytical spots from the Early Jurassic monzogranites vary from 0.28255 to 0.28280. Their εHf(t) values and two-stage model ages (TDM2) are between −4.1 to +4.8 and 918 to 1488 Ma, respectively (Figure 9). Twelve spot analyses from Middle Jurassic quartz monzonites yield 176 Hf/ 177 Hf ratios of 0.28257 to 0.28271 and εHf(t) values of −3.6 to +1.6. The corresponding TDM2 ages range from 1111 to 1444 Ma (Figure 9).  [44]; Data of magmatic rocks in the Andean orogenic arcs in d are from Asadi et al. [49]; chondrite values are from Sun and McDonough [37]). G518 is an Eastern Pontides gabbro employed as the source rock for the batch-melting modeling [50].

Zircon Lu-Hf Isotopes
A total of 21 representative zircons from four selected samples were chosen for in situ Lu-Hf isotopic analyses. The obtained data are presented in Table S7.
The 176 Hf/ 177 Hf ratios of the nine analytical spots from the Early Jurassic monzogranites vary from 0.28255 to 0.28280. Their ε Hf (t) values and two-stage model ages (T DM2 ) are between −4.1 to +4.8 and 918 to 1488 Ma, respectively (Figure 9). Twelve spot analyses from Middle Jurassic quartz monzonites yield 176 Hf/ 177 Hf ratios of 0.28257 to 0.28271 and ε Hf (t) values of −3.6 to +1.6. The corresponding T DM2 ages range from 1111 to 1444 Ma (Figure 9).  Table S8. The data for Middle Jurassic igneous rocks are the first presented for the Erguna Block.

Early Jurassic Monzogranites
Generally, granitoids are petrogenetically subdivided into I-, S-, A-and M-types [51][52][53][54]. In the study area, M-type affinity is firstly precluded due to the rare exposure of coeval and associated mafic rocks and ophiolites near the Early Jurassic monzogranites. The typical A-type granite is characterized by high (Na2O + K2O) contents and FeO T /MgO values [51]. In contrast, the studied monzogranites are subalkaline series on a TAS diagram (Figure 6a) with low (Na2O + K2O) (6.78-8.97 wt%) and FeO T /MgO (1.85-3.23) values. Besides, A-type granites commonly exhibit high Zr, Nb, Ta, Ce, Yb, and Y contents, as well as Ga/Al ratios [51,55,56]. The monzogranites from the study area have low concentrations of these trace elements and fall into I-and S-type granites fields when plotted on both (Zr + Nb + Ce + Y) versus ((K2O + Na2O)/CaO)) and (Zr + Nb + Ce + Y) versus (FeO T /MgO) diagrams (Figure 10a,b). Furthermore, the TZr values of the Early Jurassic monzogranites (771-790 °C) are lower than A-type granites (~839 °C), and similar to unfractionated I-type granites (~781 °C) [39,57,58]. The monzogranite samples are mineralogically characterized by the occurrence of amphibole and biotite (typical of I-type granites) and the absence of Al-rich and alkaline minerals, such as muscovite, cordierite, tourmaline, and corundum (typical of S-type granites) [59,60].  Table S8. The data for Middle Jurassic igneous rocks are the first presented for the Erguna Block.

Early Jurassic Monzogranites
Generally, granitoids are petrogenetically subdivided into I-, S-, A-and M-types [51][52][53][54]. In the study area, M-type affinity is firstly precluded due to the rare exposure of coeval and associated mafic rocks and ophiolites near the Early Jurassic monzogranites. The typical A-type granite is characterized by high (Na 2 O + K 2 O) contents and FeO T /MgO values [51]. In contrast, the studied monzogranites are subalkaline series on a TAS diagram (Figure 6a) with low (Na 2 O + K 2 O) (6.78-8.97 wt%) and FeO T /MgO (1.85-3.23) values. Besides, A-type granites commonly exhibit high Zr, Nb, Ta, Ce, Yb, and Y contents, as well as Ga/Al ratios [51,55,56]. The monzogranites from the study area have low concentrations of these trace elements and fall into I-and S-type granites fields when plotted on both (Zr + Nb + Ce + Y) versus ((K 2 O + Na 2 O)/CaO)) and (Zr + Nb + Ce + Y) versus (FeO T /MgO) diagrams (Figure 10a,b). Furthermore, the T Zr values of the Early Jurassic monzogranites (771-790 • C) are lower than A-type granites (~839 • C), and similar to unfractionated I-type granites (~781 • C) [39,57,58]. The monzogranite samples are mineralogically characterized by the occurrence of amphibole and biotite (typical of I-type granites) and the absence of Al-rich and alkaline minerals, such as muscovite, cordierite, tourmaline, and corundum (typical of S-type granites) [59,60].
Early Jurassic monzogranites are widespread in the Erguna Block (Figure 1c) [7,11,13]. As shown on the petrogenetic discrimination diagrams, these monzogranites and other coeval felsic rocks are dominated by I-type geochemical affinity ( Figure 10). These I-type granites are possibly generated by diverse petrogenetic mechanisms, including (1) highly fractional crystallization of mantle-originated mafic magma or mixing of mantle-originated magma and crustal materials [61,62]; or (2) partial melting of metaigneous rocks in the crust, followed by fractionation [63,64]. The formation of felsic rocks via the fractional crystallization of mantle-derived mafic magma usually demands the parental magma being voluminous [64]. The subordinate Early Jurassic mafic members in the Erguna Block preclude the possibility that the widespread monzogranites were originated by a simple fractional history. Besides, the La/Sm and Zr/Nb ratios of the studied samples increase with elevating La and Zr concentrations, respectively (Table S5), suggesting crustal partial melting played a significant role during the magma evolution. Therefore, we suggest that the partial melting of crustal metaigneous rocks was the dominant process responsible for the formation of the monzogranites.  [51]; c is after Chappell [52]; and f is after Chappell and White [55]).
Early Jurassic monzogranites are widespread in the Erguna Block ( Figure 1c) [7,11,13]. As shown on the petrogenetic discrimination diagrams, these monzogranites and other coeval felsic rocks are dominated by I-type geochemical affinity ( Figure 10). These I-type granites are possibly generated by diverse petrogenetic mechanisms, including (1) highly fractional crystallization of mantle-originated mafic magma or mixing of mantle-originated magma and crustal materials [61,62]; or (2) partial melting of metaigneous rocks in the crust, followed by fractionation [63,64]. The formation of felsic rocks via the fractional crystallization of mantle-derived mafic magma usually demands the parental magma being voluminous [64]. The subordinate Early Jurassic mafic members in the Erguna Block preclude the possibility that the widespread monzogranites were originated by a simple fractional history. Besides, the La/Sm and Zr/Nb ratios of the studied samples increase with elevating La and Zr concentrations, respectively (Table S5), suggesting crustal partial melting played a significant role during the magma evolution. Therefore, we suggest that the partial melting of crustal metaigneous rocks was the dominant process responsible for the formation of the monzogranites.
As discussed above, the monzogranites from the study area are high-K calc-alkaline I-type granites. Melting experiments have demonstrated that high-K silicic melts can be generated by partial melting of tonalitic or granodioritic rocks [65]. The resulting melts geochemically would show A-type affinity, which is significantly inconsistent with our samples. Besides, if the rocks were derived by dehydration melting of tholeiitic amphibolites, they should exhibit low K2O concentrations and  [51]; c is after Chappell [52]; and f is after Chappell and White [55]).
As discussed above, the monzogranites from the study area are high-K calc-alkaline I-type granites. Melting experiments have demonstrated that high-K silicic melts can be generated by partial melting of tonalitic or granodioritic rocks [65]. The resulting melts geochemically would show A-type affinity, which is significantly inconsistent with our samples. Besides, if the rocks were derived by dehydration melting of tholeiitic amphibolites, they should exhibit low K 2 O concentrations and Na 2 O/K 2 O values [66], which are distinct from the studied rocks. Experimental works have exhibited that the dehydration melting of slightly hydrous medium-to high-K basaltic rocks could generate high-K silicic melts [64]. The derived silicic melts are characterized by metaluminous to slightly peraluminous signatures, which are in good agreement with the monzogranites from the study area. Moreover, the monzogranites are depleted in HFSEs and HREEs, and exhibit high Rb/Sr values (0.19-0.27) and low zircon crystallization temperatures, implying a mica-bearing hydrous source [67,68]. They display moderate negative Eu anomalies with flat HREE patterns, and negative anomalies in Sr, Yb, Y, and Ti. These geochemical characteristics indicate the involvement of plagioclase, hornblende, and garnet as major residual phases [62,69]. The leaving behind of a hornblende dominated residue by dehydration melting suggests that the magma source was possibly under 8-13 kbar and 700-800 • C conditions (≤40 km depth) [70,71]. This is inconsistent with our calculated zircon saturation temperatures (771-790 • C) and Ti-in zircon geothermometry (755-900 • C with most < 800 • C). The zircons from the monzogranites experienced moderate primary high-T hydrothermal alteration during the deuteric stage of magma evolution (Figure 5c,d), with relatively good crystallinity of the rocks (Figure 2a,b), implying a deep intrusion depth of~5 km. In a Sm/Yb versus La/Sm diagram (Figure 8d), the monzogranites and coeval felsic magmatic rocks in the Erguna Block plot in the field with a thickness thinner than~45-km-thick Andean orogenic arcs. Furthermore, batch-melting modeling displays that the parental magma of the monzogranites was formed in an environment where the P-T conditions are similar to garnet-amphibole-facies metamorphism, corresponding to a depth of ≤40 km [70] (Figure 8a,c).
The monzogranites have relatively high Mg# values (49)(50)(51)(52)(53)(54)(55)(56) and MgO concentrations (0.76-2.13 wt%), which are indicative of interactions between the primary granitic melts and mantle materials. Minor mafic microgranular enclaves (MMEs) are locally found within the monzogranites, but mingling textures and complex compositional zoning in mafic minerals are seldom found, implying that the interactions were not extensive. In addition, the rocks also record a moderately scattered distribution of zircon ε Hf (t) values between −4.1 to +4.8 with corresponding T DM2 ages of 918 to 1488 Ma (Figure 9). The above field, geochemical, and isotopic signatures together indicate that the primary magma of the monzogranites was likely produced by partial melting of an accreted Meso-Neoproterozoic lower crust, with the injection of minor mantle components during the formation of the silicic magma. In conclusion, the Early Jurassic High-K calc-alkaline I-type monzogranites in the Erguna Block were likely derived by partial melting of K-rich meta-basalts within the lower part of a juvenile crust with medium-thickness.
The quartz monzonites have a narrow range of ε Hf (t) values (−3.6 to +1.6), low Mg# values (36)(37)(38)(39)(40)(41)(42)(43)(44) and MgO concentrations (0.76-1.24 wt%), with lacking MMEs and mingling textures, suggesting that magma mixing is negligible. The classic delaminated lower crust-originated adakites commonly have high Mg# values due to the involvement of mantle peridotite during magma ascending [79,80]. The low Mg# values of our samples preclude the petrogenetic processes of partial melting of the delaminated lower crust. Minerals that account for AFC-generated adakitic rocks usually have high partition coefficients (K D ) for Y and HREE, such as amphibole [75,76]. Amphibole fractionation is commonly coupled with the residue of plagioclase [81]. Besides, K D for middle-REEs is higher than HREE for amphibole [82,83]. Therefore, fractional crystallization of amphibole and plagioclase under a natural system would result in decreased Dy/Yb ratios, and negative Eu anomalies with concave-upward REE trends for the leaving melts [75,76]. The studied quartz monzonites have constant Dy/Yb ratios, negligible Eu anomalies with right-sloping REE patterns, indicating that AFC of primary basaltic magmas is insignificant. Besides, xenocrystic zircons, which are indicative of crustal assimilation during magma ascent and emplacement, were not found in the quartz monzonite samples, suggesting that the primary magma did not experience noteworthy crustal assimilation. The representative oceanic slab-originated adakitic rocks are characterized by low Rb/Sr values (0.01-0.04) and MORB-like Lu-Hf isotopic signatures [44,84]. The studied quartz monzonites yield higher Rb/Sr ratios of 0.05-0.09, and distinctively lower ε Hf (t) values between −1.3 to +1.6 in comparison with adakitic rocks originated from the oceanic slab. Moreover, slab-derived adakitic rocks are medium-K calc-alkaline series with high Na 2 O/K 2 O values (>2) [13,84]. In contrast, the samples are high-K calc-alkaline rocks and show lower Na 2 O/K 2 O ratios between 1.55-2.02. Thus, the partial melting of a subducted oceanic slab is unlikely.
Experimental works have indicated that basaltic melts in thickened regions of the continental lower crust commonly exhibit low Mg# values (<43) and MgO contents [45,63]. Similarly, the studied quartz monzonites show low Mg# values (36)(37)(38)(39)(40)(41)(42)(43)(44) and MgO concentrations (0.76-1.24 wt%). The batch melting model illustrates that the quartz monzonites likely originated from~10%-20% partial melting of an assumed 10% garnet-bearing amphibolite source (Figure 8a). Besides, they have high (Dy/Yb) N ratios of 1.22-1.30 (>1) and Sr concentrations of 847-1000 ppm (>400 ppm), but low (La/Sm) N ratios of 4.08-4.80 (<6) and Yb contents of 0.86-1.03 ppm (<2 ppm), with negligible Eu anomalies, which strongly indicate that amphibolite and garnet are dominant residual phase [85]. The significant depletion of Nb and Ta implies the existence of rutile with garnet in the residual phase [77,86]. This is further supported by the Y versus Sr/Y plots (Figure 8c), which show that the parental provinces of the primary magma under P-T conditions similar to amphibole-eclogite-facies metamorphism. Therefore, the P-T conditions for the quartz monzonites are probably >15 kbar and 750-950 • C, corresponding to a depth of ≥50 km [70,86]. This is in good agreement with our obtained zircon saturation temperatures (782-806 • C). The Ti-in zircon geothermometry shows a relatively wide range of temperatures (652-852 • C), which may imply complicated ascending and crystallization processes of the parental magma. Besides, the quartz monzonites exhibit a cumulate texture (Figure 2d). The thickened crust makes the magma to experience a long-distance ascending process before its final emplacement, with early crystallizing of amphiboles and plagioclases, followed by the protracted crystallizing of quartzs as the inter-cumulus phase. The zircons probably underwent moderate primary high-T hydrothermal alteration during the deuteric stage of magma evolution (Figure 5c,d). Combined with the relatively good crystallinity of the quartz monzonites (Figure 2c,d), a deep intrusion depth of~5 km may be constrained. In a Sm/Yb versus La/Sm diagram (Figure 8d), the quartz monzonites plot between the~45-and~60-km-thick Andean orogenic arcs. Besides, the samples have relatively clustered ε Hf (t) values and T DM2 ages range from −3.6 to +1.6 and 1111 to 1444 Ma, respectively (Figure 9), suggesting that they were possibly derived from the partial melting of an accreted Meso-Neoproterozoic lower crust. Thus, we suggest that the Middle Jurassic quartz monzonites with adakitic affinity in the Erguna Block were possibly generated by the partial melting of a thickened juvenile continental lower crust (≥50 km).

Tectonic Regime
NE China is tectonically situated in the eastern CAOB, which archives the Paleozoic-Mesozoic complex tectonic evolution of the Palaeo-Asian, Palaeo-Pacific, and Mongol-Okhotsk tectonic regimes [1,23,86]. The overprinting of multiple tectonic events results in debatable models for the tectonic history. The Early-Middle Jurassic tectonic settings for the intrusive rocks in the Erguna Block remain ambiguous, especially the Middle Jurassic, due to the poor exposure of magmatic rocks. The final closure of the Palaeo-Asian Ocean was probably finished before the Early Triassic along the Solonker-Xar Moron-Changchun-Yanji suture belt (Figure 1b) [1,9]. The resulting post-collisional Late Triassic bimodal magmatism and intrusions of A-type granitoids mainly occurred in the Lesser Xing'an and Zhangguangcai ranges, indicating an extensional setting [28,87]. The initiation of the westward subduction of the Palaeo-Pacific oceanic slab beneath the continental region of NE China was likely by the Late Triassic or Early Jurassic [88,89]. Subduction-related magmatism such as the Early Jurassic calc-alkaline volcanism in eastern Jilin-Heilongjiang provinces dominantly happened in eastern NE China [8]. The Erguna Block has geographically situated more than 900 km and 1000 km away from the Solonker-Xar Moron-Changchun-Yanji suture belt and the Palaeo-Pacific subduction zone, respectively. These long distances are generally considered to exceed the greatest extent of influence caused by far-field tectonics. On the contrary, the Erguna Block is located immediately to the southeast of the Mongol-Okhotsk suture belt, with a distance less than 200 km. Moreover, the Paleozoic-Mesozoic magmatic activities in the Erguna Block show a NE-SW-trending, which well parallels to the Mongol-Okhotsk suture belt [10,13,90]. Consequently, the spatial distribution of magmatic rocks and tectonic locations of the Erguna Block indicate that the magmatism within the block was associated with neither the Paleo-Asian nor the Paleo-Pacific tectonic regimes. In addition, the impact of the Mongol-Okhotsk tectonic regime is suggested to extend at least as far as the eastern margin of the Xing'an Block, which lies in the southeast of the Erguna Block (Figure 1b) [12,91]. Thus, we consider that the Mongol-Okhotsk tectonic regime played a dominant role in accounting for the generation of the Early-Middle Jurassic intrusive rocks in the Erguna Block.

Andean-Type Arc Setting for the Early Jurassic Monzogranites
Huge volumes of Early Jurassic igneous sequences crop out in the Erguna Block (Figure 1c) [11,12]. The sequences are dominated by felsic and intermediate rocks, with minor mafic rocks, and belong to high-K calc-alkaline and calc-alkaline series ( Figure 6). The studied monzogranites are geochemically high-K calc-alkaline I-type granites. Previous studies have shown that such granites may be produced in volcanic arc settings akin to the Andes, or post-collisional settings like the Caledonides [92]. The MOO is thought to have closed diachronously in a scissor-like manner from west to east due to the counterrotation of the Siberian Craton relative to the Mongolia-NE China continental blocks [4,93,94]. Some previous researches suggested that the closure of the MOO was by the Middle Jurassic [4,7,11,14], while others argued that the completion was during the Late Jurassic-Early Cretaceous [13,[15][16][17]. Although the timing of final closure is still controversial, the Erguna Block is considered to have been situated within an active margin setting dominated by the southward subduction of the MOOP during the Early Jurassic [7,10,12,95]. Moreover, recent studies have identified several subduction-related Early Jurassic Cu-Mo deposits in the Erguna Block, such as the Wunugetushan deposit (ca. 180 Ma) [94]. Therefore, the high-K calc-alkaline I-type monzogranites in the Erguna Block were not generated in a post-collisional setting such as the Caledonides. On tectonic setting discrimination diagrams of (Yb + Nb) versus Rb, (Yb + Ta) versus Rb, Y versus Nb, and the ternary plot of Ta×3 -Rb/30 -Hf (Figure 11a-d), all samples were plotted on the fields of volcanic granites, indicating a volcanic arc setting. Furthermore, previous works demonstrated that continental arc granites have relatively elevated Th and Ta concentrations in contrast to island arc granites [96][97][98][99][100][101][102][103][104][105][106][107][108][109][110]. The studied monzogranites contain Th and Ta similar to those felsic rocks from representative continental arc settings, such as Northern Andes [96,97], Colombian Caribbean area [98], and Colombian Andes [98], with compelling evidence from the Ta/Yb versus Th/Yb plot (Figure 11h), suggesting a continental arc setting. This is further evidenced by discriminant-function-based multi-dimensional robust diagrams (Figure 11e-g), where all samples were plotted on continental arc fields. The coeval felsic rocks within the Erguna Block also display identical features (Figure 11a-h). Thus, the Erguna Block was located in an Andean-type arc setting. A compressional environment was produced by continuous southward subduction of the MOOP beneath the trench. On a log[CaO/(K 2 O + Na 2 O)] versus SiO 2 diagram, all studied samples and most of the coexistent felsic rocks within the Erguna Block fall into the compressional fields (Figure 11i). The continental crust was relatively shortened to a medium thickness of appr. ≤40 km. The heat bonus provided by the mantle wedge triggered the partial melting of K-rich meta-basalts within the lower part of a juvenile crust, resulting in the widespread high-K calc-alkaline magmatism in the Erguna Block. Notably, an Early Jurassic bimodal volcanism was observed in the east Songliao Basin, implying an extensional environment for the Palaeo-Pacific tectonic regime [8]. This is significantly different from the contemporary tectonic setting of the Erguna Block, further indicating that the Erguna Block was subjected to the Mongol-Okhotsk tectonic regime. part of a juvenile crust, resulting in the widespread high-K calc-alkaline magmatism in the Erguna Block. Notably, an Early Jurassic bimodal volcanism was observed in the east Songliao Basin, implying an extensional environment for the Palaeo-Pacific tectonic regime [8]. This is significantly different from the contemporary tectonic setting of the Erguna Block, further indicating that the Erguna Block was subjected to the Mongol-Okhotsk tectonic regime.  [100]; (d-f) are after Verma et al. [101]; (g) is after Gorton et al. [102]; (h) is after Harris et al. [103]; (i) is after Brown [104]). Felsic rocks of representative island arc are from Izu-Bonin-Mariana arc [105,106], Kermadec arc [107], Lesser Antilles arc [108], South Sandwich arc [109], and Tonga arc [110]. Felsic rocks of the representative continental include Northern Andes [96,97], Colombian Caribbean area [98], and Colombian Andes [99].

Strong Crustal Thickening for the Middle Jurassic Quartz Monzonites
Middle Jurassic igneous rocks are rarely observed in the Erguna Block and adjacent areas. The studied quartz monzonites have geochemical characteristics that are similar to adakitic rocks that originated from a thickened continental crust. The P-T conditions for the quartz monzonites correspond to a depth of ≥50 km, indicating significant crustal shortening and thickening during the Middle Jurassic.
The quartz monzonites were produced in a volcanic arc setting, as evidenced by tectonic setting discrimination plots of (Yb + Nb) versus Rb, (Yb + Ta) versus Rb, Y versus Nb, and ternary plot of Ta×3 -Rb/30 -Hf (Figure 11a-d). Their elevated Th and Ta concentrations with respect to island arc igneous rocks imply a continental arc setting similar to the Early Jurassic monzogranites (Figure 11e-g).
Moreover, all studied samples and most of the coeval felsic rocks within the Xing'an Block plot on the continental arc setting fields on the multi-dimensional robust diagrams (Figure 11e-g). We invoke, therefore, the Middle Jurassic quartz monzonites were generated in an Andean-type arc setting, similar to the Early Jurassic monzogranites. However, the continental crust of the Erguna Block was much uplifted during the Middle Jurassic in comparison with the Early Jurassic, as recorded by the distinctively different P-T conditions of the monzogranites and quartz monzonites.
Notably, on a log[CaO/(K 2 O + Na 2 O)] versus SiO 2 diagram, the quartz monzonite samples plot on the transitional zone between the compressional and extensional fields, while the coexistent felsic rocks from the Xing'an Block fall into the extensional fields, which indicate a tectonic inversion. Moreover, in a Sm/Yb versus La/Sm diagram (Figure 8d), these felsic rocks from the Xing'an Block plot on the field with a thickness thinner than~45-km-thick Andean orogenic arcs, which may imply back-arc extension in the southeast of the Erguna Block. In central Erguna Block, intermediate volcanic rock successions with zircon U-Pb ages of ca. 162 Ma exhibit a geochemical transition from sub-alkaline to alkaline series [8,111]. Besides, in the northern Erguna Block, A-type granites with zircon U-Pb ages of ca. 155 Ma was identified in the Badaguan area [4]. The above observations demonstrate that the Erguna Block has experienced a tectonic transition from compression to extension dominated by the southward subduction of the MOOP beneath the Erguna Block during the Middle Jurassic. Our quartz monzonites yielded zircon U-Pb ages between ca. 162-174 Ma. Thus, we suggest that the significant crustal thickening may reach its ultimate stage at this time. The subsequent extensional environment dominated the formation of alkaline volcanic rocks and A-type granites.

Geodynamic Scenario
The closure of the MOO was completed by bidirectional-subduction beneath the Siberian Craton in the northwest and the Erguna and Xing'an blocks in the southeast [6,13]. The northward subduction of the MOOP could trace back to the Devonian [112]. The subduction-related magmatic activities in the Transbaikalia such as the Selenge arc archived persistent north-directed subduction from the Middle Carboniferous to the Triassic or Jurassic [92,112]. The southward subduction was initiated since the Carboniferous as evidenced by the Middle Gobi arc [113], followed by widespread Permo-Triassic subduction-related magmatic activities in the Central Mongolia, Erguna, and Xing'an blocks [12,13,31]. During the Early Jurassic, the continued subduction of the MOOP triggered the initiation of crustal thickening, accompanied by the intrusion of high-K calc-alkaline I-type monzogranites in the Erguna Block at ca. 177-198 Ma ( Figure 12). Thus, the petrogenesis of the huge volumes of Early Jurassic magmatic rocks in the Erguna Block and adjacent areas were possibly dominated by an Andean-type continental arc environment (Figure 13a). A significant magmatic quiescent is observed in the Great Xing'an Range, NE China, spanning from ca. 172-177 Ma ( Figure 12). Several tectonic scenarios may account for such a magmatic gap, including (1) final closure of the ocean basin; (2) low-angle plate subduction; (3) mid-ocean ridge subduction; or (4) slab-rollback. Previous studies suggested that the MOO remained several thousand kilometers in width in the northwest of the Erguna Block during the Early Jurassic [16]. Besides, the studied Middle Jurassic quartz monzonites are typical adakitic rocks derived from an active continental margin. The Middle Jurassic granites within the Xing'an Block are high-K calc-alkaline I-type granites, indicating a subduction-related origin [12]. Moreover, a Late Jurassic (ca. 150 Ma) subduction-related Fukeshan Cu-Mo deposit was identified within the Erguna Block recently [13]. Thus, the scenario (1) could be precluded.  Table S1. A significant magmatic quiescent is observed in the Great Xing'an Range, NE China, spanning from ca. 172-177 Ma ( Figure 12). Several tectonic scenarios may account for such a magmatic gap, including (1) final closure of the ocean basin; (2) low-angle plate subduction; (3) mid-ocean ridge subduction; or (4) slab-rollback. Previous studies suggested that the MOO remained several thousand kilometers in width in the northwest of the Erguna Block during the Early Jurassic [16]. Besides, the studied Middle Jurassic quartz monzonites are typical adakitic rocks derived from an active continental margin. The Middle Jurassic granites within the Xing'an Block are high-K calc-alkaline Itype granites, indicating a subduction-related origin [12]. Moreover, a Late Jurassic (ca. 150 Ma) subduction-related Fukeshan Cu-Mo deposit was identified within the Erguna Block recently [13]. Thus, the scenario (1) could be precluded.  Table S1.  Table S1. A significant magmatic quiescent is observed in the Great Xing'an Range, NE China, spanning from ca. 172-177 Ma ( Figure 12). Several tectonic scenarios may account for such a magmatic gap, including (1) final closure of the ocean basin; (2) low-angle plate subduction; (3) mid-ocean ridge subduction; or (4) slab-rollback. Previous studies suggested that the MOO remained several thousand kilometers in width in the northwest of the Erguna Block during the Early Jurassic [16]. Besides, the studied Middle Jurassic quartz monzonites are typical adakitic rocks derived from an active continental margin. The Middle Jurassic granites within the Xing'an Block are high-K calc-alkaline Itype granites, indicating a subduction-related origin [12]. Moreover, a Late Jurassic (ca. 150 Ma) subduction-related Fukeshan Cu-Mo deposit was identified within the Erguna Block recently [13]. Thus, the scenario (1) could be precluded. Subduction with a steep angle usually leads to a thick mantle wedge which motivates intense interactions between mantle materials, slab melts, and crustal rocks, and provides high thermal input [77,114]. In contrast, subduction with a low angle commonly produce an uneffective mantle wedge, and the far extended subducted-slab may hinder the conduction of heat from the mantle to the overriding crust [115,116]. The lithosphere would be consequently cooled, with suppression of magmatic activities. Mid-ocean ridge subduction may greatly shallow the subduction angle and cause a magmatic lull [117]. Subsequently, the buoyant oceanic topographic rises would result in the formation of a slab window, leading to a protracted magmatic activity represented by mafic and adakitic lavas near the trench [118,119]. Thus, both the low-angle plate subduction and mid-ocean ridge subduction are characterized by a shallowed subduction angle. Previous studies have demonstrated that low-angle subduction would greatly hamper the crustal thickening by forming an "eroding" margin at given stable conditions (medium rate orthogonal subduction), such as the Central American trench [18,120]. This is hard to reconcile with the case of the Erguna Block. As discussed above, the continental crust in the Middle Jurassic was significantly thicker than in the Early Jurassic, strongly indicating steep subduction. Besides, the slab window formed by mid-ocean ridge subduction would cause extensive magmatism with mafic rocks immediately after the magmatic quiescent period, contrasting with the observations in the Erguna Block. Therefore, the rollback of the subducted MOOP is the most plausible scenario that accounts for the magmatic gap.
The slab-rollback may start at ca. 177 Ma, the long-term persistent subduction resulted in the continuous growth of the continental crust in the Erguna Block. On the other hand, the gradual removal of thermal input from the mantle wedge cooled the lithosphere. The slab-rollback was terminated at ca. 174 Ma, accompanied by moderate magmatic activities that are characterized by adakitic rocks and high-K calc-alkaline I-type granites [12] (Figure 13b). The subduction-related crustal thickening may reach its ultimate stage at the same time. The consequent collapse triggered an extensional environment, resulting in the formation of alkaline volcanic sequences and A-type granitoids in the Erguna and Xing'an blocks [8,111]. The final closure of the MOO may be accomplished later at ca. 143-145 Ma [13,90], corresponding to a magmatic gap followed by intense post-collisional magmatism observed in the Great Xing'an Range ( Figure 12).

2.
The Early Jurassic I-type monzogranites were likely originated by partial melting of K-rich meta-basalts from the lower part of a juvenile crust with medium-thickness (≤40 km), with the injection of minor mantle materials. The Middle Jurassic quartz monzonites were probably produced by partial melting of a thickened juvenile continental lower crust (≥50 km).

3.
The Mongol-Okhotsk tectonic regime played a dominant role in accounting for the generation of the Early-Middle Jurassic intrusive rocks within the Erguna Block. An Andean-type continental arc setting was developed during the Early-Middle Jurassic, with continuous thickening of the continental crust. The significant crustal thickening may reach its peak during ca. 162-174 Ma, which marks the tectonic transition from compression to extension. 4.
The MOOP was subducted southward beneath the Erguna Block during the Early Jurassic, followed by slab-rollback since ca. 178 Ma. The slab-rollback was terminated at ca. 174 Ma, accompanied by moderate magmatic activities represented by adakitic rocks.  Table S9: Metadata for LA-ICPMS zircon U-Pb and trace element analyses of this study; Table S10: Metadata for whole-rock trace element analyses of this study; Table S11. Metadata for zircon Lu-Hf isotopic analyses of this study; Figure S1: Plot of zircon crystallization temperature and zircon saturation temperature.