Ore-forming mechanisms and spatio-temporal framework for intrusion-related deposits in NE China

Northeastern China is composed of the eastern part of the Central Asian Orogenic Belt and the northeastern margin of the North China Craton. It underwent two major sets of orogenic events, including the pre-Mesozoic amalgamation of several micro-continents and the Mesozoic subduction of the Paleo-Pacific plate. It hosts numerous ore deposits of dominantly porphyry and skarn types, most of which have Mesozoic ages. In this study, both mineralization types were studied, aiming to improve the understanding of ore genesis, hydrothermal evolution, mineralization mechanism, regional metallogeny as well as the geodynamic setting. Systematic zircon U-Pb and/or molybdenite Re-Os dating on five porphyry deposits (i.e., Aolunhua, Haisugou, Shabutai, Banlashan, and Yangchang) in the northern Xilamulun district indicates that the timing of the magmatism and the Mo mineralization is broadly coeval, mainly at 130-140 Ma. Major and trace element geochemistry reveals the intrusions hosting Mo-only deposits (e.g., Haisugou) have stronger crystal fractionation than intrusions hosting porphyry Cu and Cu-Mo deposits (e.g., Aolunhua), indicating that such a process may have played a role in selective enrichment of Mo. A comparison of zircon Ce/Nd ratios as a proxy for the oxidation state of magmas between mineralized and barren intrusions shows that the mineralized intrusions are associated with more oxidized magmas than the cospatial barren granites, and therefore it is proposed that higher oxygen fugacity may also be important to produce economic Mo mineralization. Whole rock Sr-Nd-Pb and zircon Hf isotopes show that these mineralized granites in Xilamulun are associated with magmas generated from three different source regions (i.e., remelting of old crust material, mixing of old crust material with depleted mantle component, and juvenile mantle-derived magmas). The variation in the origin of the magmas from which the porphyry Mo systems were generated suggests that the composition of magma sources is unlikely to have played a major role in the formation of Mo deposits. The compilation of existing geochronological data on Mo deposits in NE China, including the newly obtained data from this study, shows that Mesozoic Mo deposits (~250 to 90 Ma) widely occur in this region and are linked to three tectonic-magmatic events: (1) Triassic Mo deposits (250–220 Ma) are mainly distributed along the east-west Xilamulun fault and are related to post-collisional crustal extension following the final closure of the Paleo-Asian ocean; (2) Jurassic to Early Cretaceous Mo mineralization (200–130 Ma) displays a clear younging trend from southeast to northwest and is interpreted to be related to the northwestward flat-slab subduction of the Paleo-Pacific plate beneath the Eurasian continent that started from Early Jurassic (ca. 200 Ma); (3) Cretaceous Mo mineralization (130–90 Ma) shows a distinctly reversed migration trend from northwest to southeast, and can be explained by the coastward migration of slab rollback related lower crust delamination, asthenospheric upwelling and lithospheric thinning. For skarns, the Baiyinnuo'er Zn-Pb deposit was selected as a representative example for detailed study in this study. It is one of the largest Zn-Pb deposits in China, with 32.74 Mt resources averaging 5.44% Zn, 2.02% Pb and 31.36 g/t Ag. Several phases of igneous rocks, including Permian, Triassic and Early Cretaceous intrusions, are exposed in the mining areas, and among them the Early Cretaceous granites, which intruded into limestone of the early Permian Huanggangliang Formation, are interpreted to be the source of ore, since their Pb isotope compositions (²⁰⁶Pb/²⁰⁴Pb = 18.25–18.35, ²⁰⁷Pb/²⁰⁴Pb = 15.50–15.56 and ²⁰⁸Pb/²⁰⁴Pb = 38.14–38.32) are highly consistent with the sulfides including sphalerite, galena and chalcopyrite (²⁰⁶Pb/²⁰⁴Pb = 18.23–18.37, ²⁰⁷Pb/²⁰⁴Pb = 15.47–15.62 and ²⁰⁸Pb/²⁰⁴Pb = 37.93–38.44). Sulfur isotope values of the sulfides fall in a narrow δ³⁴S interval of -6.1 to -4.6‰ (mean = -5.4‰, n = 15), suggesting the ore-forming fluid is of magmatic origin. The deposit formed in three stages: the pre-ore stage (prograde skarn minerals with minor magnetite), the syn-ore stage (sulfides and retrograde skarn minerals including calcite and minor quartz), and the post-ore stage (late veins composed of calcite, quartz, fluorite and chlorite; cutting the above mineral assemblages).The pre-ore stage fluids trapped in pyroxene have higher temperatures (471 ± 31 °C), higher salinity (43.0 ± 3.1 wt. % NaCl eq.), and higher concentrations of Zn (~1.1 wt. %), Pb (~1.7 wt. %), and other elements (e.g., Na, K, Li, As, Rb, Sr, Cs, Ba, Cl and Br) than syn-ore mineralizing fluids (<400 °C, <12 wt. % NaCl eq., ~0.05 wt. % Zn and ~0.03 wt. % Pb). The post-ore fluids are much cooler (<270 °C; averaging ~210°C), with much lower salinity (<5.1 wt. % NaCl eq.), Zn (~38 ppm) and Pb (~19 ppm). Geochemically, the fluids of all paragenetic stages in Baiyinnuo'er are characterized by magmatic signatures based on the element ratios, which are distinctively different from basin brines. The inclusion fluids in pre-ore stage show little variation in composition between ~520 °C and ~420 °C, indicative of a closed cooling system. In contrast, the major components of the syn- and post-ore stage fluids including Cl, Na and K decrease with the temperature dropping from ~350 to <200 °C, indicating a dilution by mixing with groundwater. The metal contents in pre-ore fluid are significantly higher than in syn-ore fluid, but no mineralization occurred. This confirms that the early fluid was, although enriched in metal elements, not responsible for ore precipitation, likely due to its high temperature high salinity nature. The metal deposition was mostly due to mixing with groundwater, which caused temperature decrease and dilution that significantly reduced the metal solubility, thereby promoting metal deposition. The deposition was probably accompanied and facilitated by carbonate dissolution that buffered the acidity generated during the breakdown of Zn (Pb)-Cl complexes and the formation of sulfides. Boiling occurred in both pre-ore and early part of the syn-ore stages, but no evidence indicates that it was related to metal deposition. The current Baiyinnuo'er massive skarns contain both prograde and retrograde minerals (including ore minerals). Paragenetically, they were not formed at the same time, but could be attributed to two (or more) successive pluses of hydrothermal fluids released episodically from residual melts of a progressively downward crystallizing magma. The prograde alteration increased the permeability and porosity, and created sufficient spaces, which was essential for later metal deposition

Fluid Inclusion and Oxygen Isotope Constraints on the Origin and Hydrothermal Evolution of the Haisugou Porphyry Mo Deposit in the Northern Xilamulun District, NE China

The Haisugou porphyry Mo deposit is located in the northern Xilamulun district, northeastern China. Based on alteration and mineralization styles and crosscutting relationships, the hydrothermal evolution in Haisugou can be divided into three stages: an early potassic alteration stage with no significant metal deposition, a synmineralization sericite-chlorite alteration stage with extensive Mo precipitation, and a postmineralization stage characterized by barren quartz and minor calcite and fluorite. The coexistence of high-salinity brine inclusions with low-salinity inclusions both in potassic alteration stage (~440°C) and locally in the early time of mineralization stage (380–320°C) indicates the occurrence of fluid boiling. The positive correlations between the homogenization temperatures and the salinities of the fluids and the low oxygen isotopic compositions (δ18Ofluid < 3‰) of the syn- to postmineralization quartz together suggest the mixing of magmatic fluids with meteoric water, which dominated the whole mineralization process. The early boiling fluids were not responsible for ore precipitation, whereas the mixing with meteoric water, which resulted in temperature decrease and dilution that significantly reduced the metal solubility, should have played the major role in Mo mineralization. Combined fluid inclusion microthermometry and chlorite geothermometer results reveal that ore deposition mainly occurred between 350 and 290°C in Haisugou

Mesozoic Mo mineralization in northeastern China did not require regional-scale pre-enrichment

Northeastern China is one of the richest Mo-mineralized regions in China, with 95 Mesozoic Mo-bearing deposits and a total metal resource of >12.2 Mt Mo. The reasons behind the large-scale Mesozoic Mo mineralization in NE China remain unclear, and whether or not there was any regional-scale pre-enrichment of the source region (e.g., a Mo-rich lower crust) is still a matter of debate. In this study, whole-rock Nd and zircon Hf isotope compositions of the intrusions related to the Mo mineralization have been compiled. The results show that the isotopic compositions are highly heterogeneous among these deposits with different ages and Mo tonnages, indicating that Mo-related magmas could be derived from either ancient lower crust or relatively juvenile lower crust or via mixing of mantle-derived magmas with varying proportions of crustal melts. This suggests that different magmas, independently from their sources, can produce Mo mineralization in NE China, and, therefore, that there was probably not an unusually Mo rich basement underlying NE China. In addition, Monte Carlo simulations have been carried out to explore the magmatic processes potentially associated with the formation of magmatic-hydrothermal (porphyry- or skarn-type) Mo deposits. The results reveal that a variably large magma volume (e.g., >150 km3) was required to provide enough Mo metal to form the large (>0.1 Mt Mo) deposits. In NE China, the Jurassic Mo deposits are much more abundant and larger than the Triassic and Cretaceous deposits, which could be attributed to the specific Jurassic tectonic regime. In the Jurassic, the subduction of the Paleo-Pacific oceanic plate led to regional compression, which favors accumulation of larger amounts of magmas at depth, ultimately resulting in larger Mo deposit formation. In this study, we highlight the importance of magma volume, rather than magma source, in the formation of the Mesozoic Mo deposits in NE China

Skarn deposits of China

Skarn deposits are one of the most common deposit types in China. The 386 skarns summarized in this review contain ~8.9 million tonnes (Mt) Sn (87% of China’s Sn resources), 6.6 Mt W (71%), 42 Mt Cu (32%), 81 Mt Zn-Pb (25%), 5.4 Mt Mo (17%), 1,871 tonnes (t) Au (11%), 42,212 t Ag (10%), and ~8,500 Mt Fe ore (~9%; major source of high-grade Fe ore). Some of the largest Sn, W, Mo, and Zn-Pb skarns are world-class. The abundance of skarns in China is related to a unique tectonic evolution that resulted in extensive hydrous magmas and widespread belts of carbonate country rocks. The landmass of China is composed of multiple blocks, some with Archean basements, and oceanic terranes that have amalgamated and rifted apart several times. Subduction and collisional events generated abundant hydrous fertile magmas. The events include subduction along the Rodinian margins, closures of the Proto-Tethys, Paleo-Asian, Paleo-Tethys, and Neo-Tethys Oceans, and subduction of the Paleo-Pacific plate. Extensive carbonate platforms developed on the passive margins of the cratonic blocks during multiple periods from Neoarchean to Holocene also facilitated skarn formation. There are 231 Ca skarns replacing limestone, 15 Ca skarns replacing igneous rocks, siliciclastic sedimentary rocks, or metamorphic silicate rocks, 113 Ca-Mg skarns replacing dolomitic limestone or interlayered dolomite and limestone, and 28 Mg skarns replacing dolomite in China. The Ca and Ca-Mg skarns host all types of metals, as do Mg skarns, except for major Cu and W mineralization. Boron mineralization only occurs in Mg skarns. The skarns typically include a high-temperature prograde stage, iron oxide-rich higher-temperature retrograde stage, sulfide-rich lower-temperature retrograde stage, and a latest barren carbonate stage. The zoning of garnet/pyroxene ratios depends on the redox state of both the causative magma and the wall rocks. In an oxidized magma-reduced wall-rock skarn system, such as is typical of Cu skarns in China, the garnet/pyroxene ratio decreases, and garnet color becomes lighter away from the intrusion. In a reduced intrusion-reduced wall-rock skarn system, such as a cassiterite- and sulfide-rich Sn skarn, the skarn is dominated by pyroxene with minor to no garnet. Manganese-rich skarn minerals may be abundant in distal skarns. Metal associations and endowment are largely controlled by the magma redox state and degree of fractionation and, in general, can be grouped into four categories. Within each category there is spatial zonation. The first category of deposits is associated with reduced and highly fractionated magma. They comprise (1) greisen with Sn ± W in intrusions, grading outward to (2) Sn ± Cu ± Fe at the contact zone, and farther out to (3) Sn (distal) and Zn-Pb (more distal) in veins, mantos, and chimneys. The second category is associated with oxidized and poorly to moderately fractionated magma. Ores include minor porphyry-style Mo and/or porphyry-style Cu mineralization ± Cu skarns replacing xenoliths or roof pendants inside intrusions, zoned outward to major zones of Cu and/or Fe ± Au ± Mo mineralization at the contact with and in adjacent country rocks, and farther out to local Cu (distal) + Zn-Pb (more distal) in veins, mantos, and chimneys. Oxidized and highly fractionated magma is associated with porphyry Mo or greisen W inside an intrusion, outward to Mo and/or W ± Fe ± Cu skarns at the contact zone, and farther to Mo or W ± Cu in distal veins, mantos, and chimneys. The final category is associated with reduced and poorly to moderately fractionated magma. No major skarns of this type have been recognized in China, but outside China there are many examples of such intrusions related to Au-only skarns at the contact zone. Reduced Zn-Au skarns in China are inferred to be distal parts of such systems. Tungsten and Sn do not occur together as commonly as was previously thought. The distal part of a skarn ore system may transition to carbonate replacement deposits. Distal stratabound mantos and crosscutting veins/chimneys may contain not only Zn-Pb but also major Sn, W, Cu, Mo, and Au mineralization. The Zn-Pb mineralization may be part of either an oxidized system (e.g., Cu, Mo, Fe) or a reduced system (e.g., Sn). In China, distal Zn-Pb is more commonly related to reduced magmas. Gold and W may also be related to both oxidized and reduced magmas, although in China they are more typically related to oxidized magma. There are numerous examples of distal mantos/chimneys that continuously transition to proximal skarns at intrusion-wall-rock contact zones, and this relationship strongly supports the magmatic affiliation of such deposits and suggests that distal skarns/carbonate replacement deposits systems should be explored to find more proximal mineralization. Carbonate xenoliths or roof pendants may host the majority of mineralization in some deposits. In contact zones, skarns are better developed where the intrusion shape is complicated. The above two skarn positions imply that there may be multiple skarn bodies below drill interceptions of intrusive rocks. Many of the largest skarns for all commodities in China are related to small or subsurface intrusions (except for Sn skarns), have multiple mineralization centers, are young (<~160 Ma), and have the full system from causative intrusion(s) to distal skarns or carbonate replacement extensions discovered. Chinese skarn deposits fall in several age groups: ~830, ~480 to 420, ~383 to 371, ~324 to 314, ~263 to 210, ~200 to 83, ~80 to 72, and ~65 to 15 Ma. They are typically associated with convergent plate boundaries, mostly in subduction settings but also in collisional settings. Seven major skarn metallogenic belts are recognized based on skarn geographic location and geodynamic background. In subduction settings, skarns may form in a belt up to 4,000 km long and 1,000 km inland, with skarns continuously forming for up to 120 m.y., e.g., the eastern China belt. In most other belts, skarns form in 5- to 20-m.y. episodes similar to the situation in South America. In collisional settings, skarns may form up to 50 m.y. after an ocean closure, and the distance to the collisional/accretionary boundary may extend to ~150 km inland. The size of collision-related skarns may be as large as the largest skarns related to oceanic crust subduction. Older suture zones may be favorable sites for younger mineralization, for example, the Triassic Paleo-Tethys suture between the North and South China blocks for the younger and largest skarn cluster of the Middle-Lower Yangtze belt in the eastern China belt, and the Triassic sutures in southwestern China for Cretaceous to Tertiary mineralization

Fluid compositions reveal fluid nature, metal deposition mechanisms, and mineralization potential: an example at the Haobugao Zn-Pb skarn, China

Fluid inclusion compositions obtained from laser ablation–inductively coupled plasma–mass spectrometry at the Haobugao Zn-Pb skarn in northeastern China provide constraints on fluid origin, evolution, and metal deposition mechanisms and an example of evaluating mineralization potential. Metal concentrations in the prograde fluids were high (up to 1.4 wt% Zn and 1.8 wt% Pb) but remained in solution, likely due to the high temperatures (440–575 °C) and salinities (35.4–45.3 wt % NaCl equivalent). Absolute concentrations of elements (e.g., Rb and Na) and mass ratios (e.g., Zn/Na and K/Na) reveal that the early, prograde fluids were magmatic, consistent with the oxygen isotope composition of fluids (δ18OH2O = 5.5‰–8.5‰). Later mixing with a meteoric fluid caused dilution and Zn-Pb deposition, as revealed by lowered element concentrations and Pb/(Na + K) and Zn/(Na + K) ratios in the sulfide-stage fluid inclusions. Elevated Ca/K ratios in sphalerite-hosted inclusions indicate fluid-carbonate reactions that buffered fluid pH, also facilitating Zn-Pb precipitation. Although cassiterite and molybdenite occur locally at Haobugao, mass balance calculation shows low metal endowment (maximum 2900 t Sn and 2200 t Mo) of the system. Furthermore, the generally unchanged Sn/(Na + K) and Mo/(Na + K) ratios from pre- to late-mineralization fluids suggest that the fluids were never saturated in Sn and Mo. Therefore, finding much Sn or Mo at Haobugao is unlikely. This demonstrates a potential tool for evaluating the metal endowment of a mineral prospect, which may guide exploration

Composition and evolution of fluids forming the Baiyinnuo'er Zn-Pb Skarn Deposit, Northeastern China: insights from laser ablation ICP-MS study of fluid inclusions

The Baiyinnuo'er skarn deposit is one of the largest Zn-Pb deposits in northeastern China, with 32.74 million metric tons (Mt) resources averaging 5.44% Zn, 2.02% Pb, and 31.36 g/t Ag. The deposit formed in three stages: the preore stage (prograde skarn minerals with minor magnetite), the synore stage (sulfides and retrograde skarn minerals including calcite and minor quartz), and the postore stage (late veins composed of calcite, quartz, fluorite, and chlorite; cutting the above mineral assemblages). In this study we analyzed the composition of single fluid inclusions using laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) to (1) determine the composition of the fluids and the evolution through the stages, (2) infer the fluid and metal sources, and (3) explore the metal deposition mechanisms. The preore fluids trapped in pyroxene have higher homogenization temperatures (432°–504°C), higher salinity (36.5–46.1 wt % NaCl equiv), and higher concentrations of Zn (~0.9 wt %), Pb (~1.4 wt %), and other elements (e.g., Na, K, Li, As, Rb, Sr, Cs, Ba, Cl, and Br) than synore mineralizing fluids (<370°C, <10 wt % NaCl equiv, ~450 ppm Zn, and ~290 ppm Pb). The postore fluids show lower temperatures (<250°C) and a rather dilute composition (<4 wt % NaCl equiv, ~33 ppm Zn, and ~24 ppm Pb). Geochemically, the fluids of all paragenetic stages in Baiyinnuo'er have magmatic signatures based on the element mass ratios, including elevated K/Na, Zn/Na, and Rb/Na ratios, lower Ca/K ratios, and combined Cl/Br-Na/K ratios, which are distinctively different from basinal brines. Inclusion fluids in preore stage show little variation in composition between ~510° and ~430°C, indicative of a closed cooling system. In contrast, the major components of the syn- and postore fluids, including Cl, Na, and K, decrease and correlate with a drop of homogenization temperatures from ~370° to ~200°C, indicating a dilution by mixing with groundwater. The Baiyinnuo'er mineralizing fluids (trapped in sphalerite) have higher Ca/K mass ratios (avg ~0.78) than other proximal magmatic hydrothermal systems (0.06–0.52) but lower than that of the distal El Mochito skarn (avg ~6.4), probably reflecting a gradually weakened magmatic signal away from the causative intrusions. The metal contents in preore fluids are significantly higher than those in synore fluids, but no mineralization occurred. This confirms that the early fluids were, although enriched in metals, not responsible for ore precipitation, most likely due to their high temperature and high salinities. One important factor controlling Zn-Pb mineralization was mixing with groundwater, which resulted in temperature decrease and dilution that significantly reduced the metal solubility, thereby promoting metal deposition. Another main driving force was the interaction with carbonate wall rock that buffered the acidity generated during the breakdown of Zn and (Pb)-Cl complexes and the precipitation of sulfides. Phase separation occurred in both the preore and the early part of the synore stages, but no evidence indicated that it caused metal deposition. The prograde minerals and retrograde minerals (including ore minerals) coexisting in the same samples could have been caused by two (or more) successive pulses of hydrothermal fluids released from residual melts of a progressively downward crystallizing magma. Each fluid produced a series of proximal high-temperature prograde to distal low-temperature assemblages, with the lower temperature assemblages of later fluids overprinting the higher temperature assemblages at most locations

Zircon U–Pb geochronology and Sr–Nd–Pb–Hf isotopic constraints on the timing and origin of Mesozoic granitoids hosting the Mo deposits in northern Xilamulun district, NE China

Located in the east section of the Central Asian orogen in northeastern China, the Xilamulun district comprises several newly discovered molybdenum deposits, primarily of porphyry type and Mesozoic ages. This district is divided by the Xilamulun fault into the southern and the northern parts. In this paper, we present new zircon U–Pb dating, trace elements and Hf isotope, and/or whole rock Sr–Nd–Pb isotopic results for the host granitoids from three Mo deposits (Yangchang, Haisugou and Shabutai) in northern Xilamulun. Our aim is to constrain the age and petrogenesis of these intrusions and their implications for Mo mineralization. Zircon U–Pb LA-ICP-MS dating shows that the monzogranites from the Shabutai and Yangchang deposits formed at 138.4 ± 1.5 and 137.4 ± 2.1 Ma, respectively, which is identical to the molybdenite Re–Os ages and coeval well with the other Mo deposits in this region, thereby indicating an Early Cretaceous magmatism and Mo mineralization event. Zircon Ce/Nd ratios from the mineralized intrusions are significantly higher than the barren granites, implying that the mineralization-related magmas are characterized by higher oxygen fugacity. These mineralized intrusions share similar zircon in-situ Hf and whole rock Sr–Nd isotopic compositions, with slightly negative to positive εHf(t) ranging from − 0.8 to + 10.0, restricted εNd(t) values from − 3.7 to + 1.6 but a little variable (87Sr/86Sr)i ratios between 0.7021 and 0.7074, indicative of formation from primary magmas generated from a dominantly juvenile lower crust source derived from depleted mantle, despite diverse consequent processes (e.g., magma mixing, fractional crystallization and crustal contamination) during their evolution. The Pb isotopes (whole rock) also show a narrow range of initial compositions, with (206Pb/204Pb)i = 18.03–18.88, (207Pb/204Pb)i = 15.48–15.58 and (208Pb/204Pb)i = 37.72–38.28, in agreement with Sr–Nd–Hf isotopes reflecting the dominance of a mantle component. An integration of geochronological and isotopic data points to three stages of Mo mineralization events (i.e., Triassic, Late Jurassic, and Early Cretaceous) associated with magmas generated from three different source regions in the Xilamulun district, NE China. The variation in the origin of the magmas from which the porphyry Mo systems were generated suggests that the formation of Mo deposit lies not in the composition of magma sources. Other factors, including high magma oxygen fugacity, may have taken a fundamental role in Mo enrichment and subsequent mineralization

Regional metallogeny of Mo-bearing deposits in northeastern China, with new Re-Os Dates of porphyry Mo deposits in the northern Xilamulun district

With many new discoveries of giant Mo deposits in the past decade, northeastern China has become the largest molybdenum mineralization region in China. The Xilamulun district is in the southern part of northeastern China and contains 25 Mo deposits and prospects, dominantly of porphyry type. These deposits are mostly along the E-striking Xilamulun fault. Re-Os dating of hydrothermal molybdenite from four deposits in the northern Xilamulun district revealed mineralization ages of 137 to 129 Ma (129.4 ± 3.4 Ma for Aolunhua, 135.3± 2.6 Ma for Shabutai, 136.4 ± 0.8 Ma for Haisugou, and 136.1 ± 6.6 Ma for Banlashan), in general agreement with the crystallization ages of their host granitoids.\ud \ud The compilation of existing data on Mo-bearing deposits in northeastern China, including the new data of this study, shows that Mesozoic Mo deposits in this region have a wide range of ages, from ~250 to 90 Ma. We propose that they are linked to three tectonic-magmatic events: (1) Triassic Mo deposits (250–220 Ma) are mainly distributed along the E-trending Xilamulun fault and are related to postcollisional crustal extension following the final closure of the Paleo-Asian ocean; (2) Jurassic Mo mineralization (200–140 Ma) displays a clear younging trend from southeast to northwest, coincides with the regional magmatism trend, and is interpreted to be related to the northwestward flat-slab subduction of the Paleo-Pacific plate beneath the Eurasian continent that started in the Early Jurassic (ca. 200 Ma); (3) Cretaceous Mo mineralization (140–90 Ma) shows a reversed migration trend from northwest to southeast, and can be explained by the coastward migration of slab rollback related to lower crust delamination, asthenospheric upwelling, and lithospheric thinning in eastern China. The spatial-temporal distribution of Mesozoic Mo mineralization is important for regional metallogeny and exploration. Recently numerous epithermal Au (Cu, Mo) deposits have been recognized in the southeast margin of northeastern China. According to this study, we predict that hidden porphyry/skarn systems could be located beneath these epithermal deposits