The Wunugetushan porphyry Cu–Mo deposit is located in northeastern China. The deposit lies within the Mongolia–Erguna metallogenic belt, which is associated with the evolution of the Mongol–Okhotsk Ocean. The multiple episodes of magmatism in the ore district, occurred from 206 to 173 Ma, can be divided into pre-mineralization stage (biotite granite), mineralization stage (monzogranitic porphyry and rhyolitic porphyry), and post-mineralization stage (andesitic porphyry). The biotite granite has (⁸⁷Sr/⁸⁶Sr)_i values of 0.704105–0.704706, ε_Nd(t) values of −0.67 to −0.07, and ε_Hf(t) values of −0.4 to 2.8, yielding Hf two-stage model ages (T_DM2) 1250–1067 Ma, and Nd model ages of 1.04–0.96 Ga, indicating that the pre-mineralization magmas were generated by the remelting of Neoproterozoic juvenile crustal material. The monzogranitic porphyry has (⁸⁷Sr/⁸⁶Sr)_i values of 0.704707–0.706134, ε_Nd(t) values of 0.29–1.33, and ε_Hf(t) values of 1.0–2.9, yielding T_DM2 model ages of 1173–1047 Ma. The rhyolitic porphyry has (⁸⁷Sr/⁸⁶Sr)_i ratio of 0.702129, ε_Nd(t) value of −0.21, and ε_Hf(t) values of −0.5 to 7.1, T_DM2 model ages from 1269 to 782 Ma. These results show that the magmas of mineralization stage were generated by the partial melting of juvenile crust mixed with mantle-derived components. The andesitic porphyry has (⁸⁷Sr/⁸⁶Sr)_i ratio of 0.705284, ε_Nd(t) value of 0.82, and ε_Hf(t) values from 4.1 to 7.4, indicating that the post-mineralization magma source contained more mantle-derived material. The Mesozoic Cu–Mo deposits which genetically related to Mongol–Okhotsk Ocean were temporally distributed in Middle to Late Triassic (240–230 Ma), Early Jurassic (200–180 Ma), and Later Jurassic (160–150 Ma) period. The Middle Triassic to Early Jurassic Cu–Mo mineralization was dominated by Mongol–Okhotsk oceanic plate southeast-directed subducted beneath the Erguna massif. The Later Jurassic Cu–Mo mineralization was controlled by the continent–continent collision between Siberia plate and Erguna massif.

The results of SHRIMP U-Pb ages and in situ Hf isotope of zircons from three granites in the southern Manzhouli region of northeast China, provide new data to understand the subduction process of Mongol-Okhotsk Plate beneath the Erguna massif. SHRIMP U-Pb zircon geochronology results yield an age of 265.5 Ma (middle Permian) for fine-grained monzogranite. Rocks from the Early–Middle Triassic are mainly granodiorite (247.4 ± 4.6 and 249.3 ± 4.9 Ma), the granites are with SiO₂ = 60.0–77.4 wt.%, Al₂O₃ = 12.3–16.8wt.% and Na₂O/K₂O = 0.7–1.9. Chemically, they are metaluminous to peraluminous and belong to the high-K calc-alkaline series. Enrichments in the large ion lithophile elements (e.g., Rb, Ba, and K) and depletions in the high field strength elements (e.g., Nb, Ta, and Ti) are typical for these rock types. The monzogranite (~265 Ma) and granodiorite (~247 Ma) contain zircons with εHf(t) values of 6.3–8.5 and 5.1–7.9, yielding T_DM2 model ages of 888–752 and 958–774 Ma, respectively. These geochemical and zircon Hf isotopic data indicate that primary magmas for Middle Permian–Early Triassic granites crystallized from primary magmas generated by Neoproterozoic crustal materials, formed in an active continental margin setting. The andesite of the Gegenaobao formation is similar with the Izu–Bonin–Mariana arc, relating to subduction initiation. Based on the characteristics of exposed rocks and zircon U-Pb ages of andesite and granitoid rocks in the study area, we conclude the onset subduction of Mongol-Okhotsk Plate beneath the Erguna massif may occur at early-middle Permian.