Geochemical Study of Cretaceous Magmatic Rocks and Related Ores of the Hucunnan Cu–Mo Deposit: Implications for Petrogenesis and Poly-Metal Mineralization in the Tongling Ore-Cluster Region

: The Hucunnan porphyry- and skarn-type Cu–Mo deposit is located in the south of the central Shizishan ore field of the Tongling ore-cluster region. The intrusive Hucunnan granodiorite, outcropping in this deposit, has adakitic geochemical features, and its magma is proposed to have originated from partial melting of the oceanic crust mixed with mantle-derived materials. The porphyry-type orebody is hosted in the granodiorite, whereas the skarn-type orebody occurs in the contact zones of intrusions and country rocks. The δ 34 S values of pyrite from the skarn orebodies ranged from +3.9 to +4.7‰ (avg. +4.3‰, n = 6), while those of the porphyry orebodies ranged from +5.1 to +6.2‰ (avg. +5.6‰, n = 4). 208 Pb/ 204 Pb, 207 Pb/ 204 Pb, and 206 Pb/ 204 Pb ratios of the pyrites from the skarn orebodies were 38.04–38.45 (avg. 38.26), 15.55–15.66 (avg. 15.59), and 18.16–18.54 (avg. 18.44), respectively ( n = 6). The pyrites in the porphyry orebodies had 208 Pb/ 204 Pb, 207 Pb/ 204 Pb, and 206 Pb/ 204 Pb ratios of 38.24–38.36, 15.51–15.662, and 18.10–18.41, respectively (avg. 38.32, 15.58, 18.22; n = 4), respectively. The metallogenic model ages from Re–Os isotopic dating were 138.7 ± 1.9 and 140.0 ± 2.8 Ma, respectively. Geochemical data indicate that the ore-forming fluids in the skarn stage are characterized by high temperature, low acidity, and high oxygen fugacity, and the ore-forming materials were mainly from magma and partly from stratum, proving that the skarn orebody has more stratum materials than the porphyry orebody.


Introduction
Porphyry and skarn deposits are widely distributed in the world and represent the most important source of Cu, W, Fe, Mo, Pb-Zn, etc. Various petrogenetic and metallogenic models that involved porphyry and skarn deposits have been proposed over the past few decades [1][2][3][4][5][6][7]. Porphyry-skarn metallogenesis is a dynamic process, and the development of deposits depends on parameters such as the host rock types, geochemistry of the causative intrusions, and local tectonic history. These deposits can form in either convergent arc systems (island arcs or continental arcs) or collisional settings and they are closely related to intermediate to acidic intrusions [8,9]. With the development of adakite rocks, some scholars believe that the high oxygen fugacity of adakite magma is closely related to porphyry copper deposits [10].
The Middle-Lower Yangtze River Belt (MLYRB) is a well-known metallogenic belt in eastern China that contains more than 200 polymetallic Cu-Fe-Au-Mo deposits. It hosts seven large-scale ore-cluster regions from northeast to southwest: Ningzhen (Cu, Fe, Pb, Zn), Ningwu (Fe), Tongling (Cu, Au, Pb, Zn), Luzong (Fe, Cu), Anqing-Guichi (Cu, Fe), Jiurui (Cu, Au), and E'dong (Fe, Cu) [11][12][13][14][15][16][17][18]. Among them, the Tongling ore-cluster region has high economic value with a high degree of geological exploration and intense studies ( [19][20][21][22][23][24][25][26][27][28][29]; Figure 1a), and it is the birthplace of the theory of "stratabound skarn deposit" [11]. Skarn ore deposits are one of the most important ore types in the Tongling ore-cluster region, and the geochronological, geochemical, and fluid evolution features of this type have been reported in numerous studies [30][31][32]. However, in recent years, porphyry ore bodies were found in the depth, and there have been relatively few studies on the mineralization of porphyry-skarn type deposits in this area [33][34][35]. As typical porphyry-skarn Cu-Mo polymetallic deposits, the Hucunanan deposit in the Tongling ore-cluster region provides a good example for revealing the magmatic evolution and enrichment and migration of metallogenic materials during the skarn formation process. Previous studies held that the ages of quartz diorites and granodiorite of Hucunan deposit are 125 Ma and 137.5 Ma, respectively [36]. The H-O isotopic signatures and the study of fluid inclusions show that the ore-forming fluids were dominated by magmatic water in the early stages, and gradually mixed with circulating meteoric water in the late stages [37]. However, there is still a lack of research on the magmatic evolution and source of ore-forming materials.
This study systematically investigated the age, physical conditions (temperature, pressure, and oxygen fugacity), source of mineralization, and petrogenesis of the intrusive rocks in the Hucunnan deposit. Here, the geochemical characteristics of intrusive rock were summarized and the evolution and source of magma were discussed. We also discussed the S-Pb isotopic compositions of metallic sulfides and the elemental compositions of pyrite and garnet. They can be used to trace the evolution of the environment and the source of ore-forming materials, and play important roles in judging the physical and chemical conditions of ore-forming hydrothermal solution and determining the genetic types of ore deposits [38][39][40][41]. The new data provide evidence related to the evolutionary history of oreforming hydrothermal system, all of which are necessary for understanding the ore-forming process.

Geological Setting
The Tongling ore-cluster region is along the lower Yangtze depression in the northern Yangtze block and belongs to the secondary uplift depression, which is also known as the Tongling uplift [11]. Apart from the absence of Lower and Middle Devonian strata, the main stratigraphic system of this region is from Silurian to Triassic (Figure 1b) [42]. The Cu-Au orebodies are mainly located in neriticfacies and littoral-facies carbonate sequences, e.g., the Carboniferous Huanglong and Chuanshan Formations and the Lower Permian Nanlinghu Formation. There are more than 70 magmatic plutons in the Tongling ore-cluster region. They are distributed from east to west for nearly 25 km in the Tongling-Daijiahui magmatic belt (Figure 1b). These magmatic rocks are divided into four rock assemblages [43,44]: (1) Pyroxene diorite-pyroxene monzodiorite combination, (2) quartz dioritequartz monzodiorite combination, (3) granodiorite combination, and (4) K-feldspar granite. The ages of granodiorites are concentrated in the range of 144-137 Ma, whereas pyroxene diorite-pyroxene monzodiorite and quartz diorite-quartz monzodiorite range from 143-137 Ma and 140-136 Ma, respectively [19]. The quartz monzodiorite and granodiorite are closely related to copper mineralization, and belong to high potassium calc-alkaline and olivine anorthite [25]. There are three groups of faults in the basement of the Tongling ore-cluster region with trending directions of EW, SN, and N-NE. In addition, there are three groups of shallow fault systems, i.e., NE, NS, and NW. These structures controlled the distribution of rocks and deposits in this area.
There are more than 200 ore deposits in the Tongling ore-cluster region, which are divided into five ore fields: Tongguanshan (Cu, Au, Fe), Shizishan (Cu, Au, S, Pb, Zn), Xinqiao-Shujiaidan (Cu, S, Au), Fenghuangshan (Cu, Au, Fe), and Shatanjiao (Pb, Zn, Au, Cu) ( Figure 1b) [11]. The metallogenic age gradually changes from east to west, the metallogenic age difference of the adjacent ore fields is about 2 Ma, and the metallogenic ages of the whole ore-cluster region are between 141 and 137 Ma [19]. According to different mineralization, the deposits in the area can be divided into porphyry, sedimentation-hydrothermal superposition reworked type, skarn, hydrothermal, and weathering types (Table 1).

Geology of Deposit
The Hucunnan Cu-Mo deposit, a medium-sized porphyry-and skarn-type deposit, is located in the south of Shizishan orefield in the Tongling ore-cluster region (Figure 1b). The exposed surface strata, exposed in the mining area, are mainly Middle Triassic Dongmaanshan Formation (T2d) limestone and dolomite (Figure 2a), and the deep part is Early Triassic (T1) and Permian (P) sedimentary rocks.
The fold structures in the mining area are Qingshanjiao anticline and Zhucun syncline. The main faults can be classified into three groups, i.e., SN, NW-W, and N-NE directions, which are important rock-and ore-controlling structures in this area.
There are two types of orebodies in the Hucunnan deposit, i.e., skarn and porphyry type. The porphyry-type orebody is located in the granodiorite, and the skarn-type orebody is located in the contacts of intrusions and country rocks (Qixia and Gufeng Formation) (Figure 2b). The Cu mineralization is dominant in the adjacent contact zones of the southwest and northeast sides of the intrusion. Generally, the orebodies strike N 30° and dip 10°-25° to the southeast. The Cu orebody is 6. Metal minerals are mainly chalcopyrite, pyrite, and molybdenite, followed by magnetite, sphalerite, and galena. Gangue minerals are mainly quartz, garnet, pyroxene, calcite, plagioclase, and potassium feldspar, followed by biotite, sericite, epidote, and chlorite. The ore structure is diverse, including automorphic-semi-automorphic, heteromorphic, metasomatic, and solid solution separation structures. The structures are mainly disseminated, veined, and reticulated (Figure 3d-f).

Sample Collection
All samples were collected from drills (Nos. ZK907 and ZK927) in the Hucunnan Cu-Mo deposits from −1400 m to −1600 m.
The five samples of granodiorite were from No. ZK927. The rocks were relatively fresh, graywhite, with idiomorphic-semi-idiomorphic granular and block structure. In this work, whole-rock major and trace elements were tested, and the rock-forming minerals (biotite) were analyzed by electron probe micro-analysis (EPMA, JEOL, Beijing, China).
Skarn of different depths and two metallogenic types of pyrite-bearing samples were selected to grind probe plates. Seventeen garnets of skarn, ten pyrites of skarn type (PySK 1-10; Figure 4a-c), and seven pyrites of porphyry type (PyGD 1-7; Figure 4d-f) were selected for EPMA. Ten samples were selected for S-Pb isotope analysis. Six samples were skarn type (TW 01-06; Figure 4a) and the others were porphyry type (TW 07-10; Figure 4d).
The garnets were usually reddish brown, with mainly rhombic dodecahedron or tetragonal trihedron self-shaped/semi-self-shaped granular structures. The particle size was generally 3-6 mm, with a small part reaching more than 1 cm, and the band structure was obvious. The pyrites were usually bright yellow, automorphic/semi-automorphic, disseminated, and reticulated.

Whole-Rock Major and Trace Elements
Whole-rock major and trace elements of fresh samples were analyzed by ALS Minerals-ALS Chemex (Guangzhou, China).
Major elements were determined by ME-XRF26 spectrometry. A prepared sample was fused with lithium metaborate-lithium tetraborate flux, which also included an oxidizing agent (lithium nitrate), and then poured into a platinum mold. The resultant glass disk was analyzed by X-ray fluorescence (XRF) spectrometry. The XRF analysis was determined in conjunction with a loss-onignition at 1000 °C. The resulting data from both determinations were combined to produce a total.
Trace elements and rare earth elements (REE) were determined by inductively coupled plasmamass spectrometry (ICP-MS) and inductively coupled plasma-atomic emission spectroscopy (ICP-AES). First, we prepared two pulp subsamples. One subsample was digested with perchloric, nitric, and hydrofluoric acids. The residue was leached with diluted hydrochloric acid and diluted to volume. The solution was then analyzed by ICP-MS for ultra-trace-level elements. The same digestion solution was also analyzed by ICP-AES for trace-level elements. Results were corrected for spectral interelement interference. The other subsample was added to lithium metaborate/lithium tetraborate flux, mixed well, and fused in a furnace at 1025 °C. The resulting melt was then cooled and dissolved in an acid mixture containing nitric, hydrochloric, and hydrofluoric acids. This solution was then analyzed by ICP-MS. According to the actual situation of the sample and the digestion effect, the comprehensive value was the final test result.

EPMA
The element composition of pyrite was determined using a JEOL JXA-8230M electron microprobe at the School of Resources and Environmental Engineering, Hefei University of Technology, Hefei, China. The analytical conditions were as follows: Accelerating voltage of 15 kV, probe current of 20 nA, and peak diameter of 5 μm.

Sulfur Isotope
The S isotope was tested at the Chinese Academy of Geological Sciences, Beijing, China. Samples contained no more than 100 ug of sulfur and three to five times the amount of V2O5, wrapped in 9-mm × 5-mm tin cups. They were automatically sampled in the combustion reactor with an oxygen supply and sufficient combustion under 1020 °C. All gas was combusted under helium load flow into and through the WO3 filler and redox reaction of Cu wire tube, and all gas was oxidized adequately. At the same time, a small amount of SO3 was generated through the reduction of SO2 by the Cu wire. Gas was separated from SO2 and other impure gases by a column (PN 260 070 80 Sulphur Separation Column for IRMS/HT) and then tested by mass spectrometry. Three international standard substances, IAEA-SO-5, IAEA-SO-6, and NBS 127, were used. The accuracy of standard sample analysis was better than 0.2‰.

Pb Isotope
The Pb isotope was tested at the Chinese Academy of Geological Sciences, Beijing, China. For determination of Pb isotope composition, all the methods were static. 202 Hg + was measured to monitor the interference of 204 Hg + to 204 Pb + in argon. 205 Tl + and 203 Tl + were measured to use 203 Tl/ 205 Tl as an external standard calibration instrument for mass fractionation effect. Before testing the samples, the parameters of Neptune Plus were optimized using NIST 981 200 ug/L standard solution, including plasma parameters (rectangular tube position and carrier gas velocity) and ion lens parameters, in order to achieve maximum sensitivity. The sample after chemical separation underwent mass spectrometry with 2% HNO3 solution, and the signal intensity of 208 Pb + was about 8 V (the concentration of Pb in solution was about 200 ug/L). The sample was injected by free atomizer. The lead isotope ratio was normalized by 203 Tl/ 205 Tl = 0.418922.

Re-Os Isotope
The Re-Os isotope analysis of molybdenite was completed at the Chinese Academy of Geological Sciences, Beijing, China. The isotope dilution Carius tube reverse water sample technique was used to seal and melt the sample for 24 h at 220 °C [46], direct distillation was used to separate and enrich Os [47], acetone extraction was used to separate and purify Re [48], and four-stage rod mass spectrometry (thermoelectric company x-series) was used to determine isotope ratios.

Whole-Rock Geochemistry
The major and trace element results are listed in Supplementary Table S1
The S isotope composition of the different orebody types at Hucunnan differed slightly, but the two generally had little change and were close, which indicates that different types of orebodies may have the same S source but were slightly different in the evolution process.

Pb Isotope
The Pb isotope composition of different types of orebodies is shown in Supplementary

Petrogenesis of the Hucunnan Granodiorite
The granodiorite of the Hucunnan Cu-Mo deposit is consistent with the magmatic rock related to copper ore in the Shizishan ore field, both of which belong to the high-K calcium alkaline series.
The chemical compositions of biotite can also reflect oxidation conditions during magma crystallization. The concentrations of Fe 3+ , Fe 2+ , and Mg 2+ in biotite, which is paragenetic with K-feldspar and magnetite, can be used to calculate magmatic oxygen fugacity [54][55][56][57][58]. From the Fe 3+ -Fe 2+ -Mg diagram of biotite in the intrusion (Figure 8a), it can be seen that the sample points of biotites fall between the two buffer lines of Ni-NiO and Fe2O3-Fe3O4, indicating that biotite was crystallized under the condition of high oxygen fugacity [57]. The contents of MgO and FeO in biotite can determine the source of magmatic material. In the MgO-FeO/(FeO + MgO) diagram for biotites (Figure 8b), data are mainly plotted in the crust-mantle mixed-source area, indicating that Hucunnan granodiorite has the characteristics of crust-mantle mixed source.
The Hucunnan granodiorite was enriched in K, Al, and Na; enriched in large ion lithophile and light rare earth elements; depleted in high-field-strength elements Nb, Ta, and Ti; and obviously enriched in Sr and Ba, but Rb was low. The Rb/Y-Nb/Y diagram (Figure 9a) shows that the distribution of Hucunnan granodiorite was basically along the line of the evolution of crustal contamination or subduction zones. These characteristics reflect that the original magma came from the lithospheric mantle, and the lithospheric mantle source area may have been caused by plate subduction and may also reflect crust mixing during the magma rising process.
In the La/Sm-La diagram (Figure 9b), the La/Sm ratio and La content of the sample show an obvious linear increasing trend on the partial melting evolution line, indicating that partial melting occurred in the process of magma evolution. In the Harker diagram ( Figure 10), Al2O3, CaO, P2O5, TiO2, and MgO decreased with the increase of SiO2, indicating that with the evolution of magma, hornblende, pyroxene, sphene, biotite, and other minerals in the magma gradually underwent precipitation crystallization separation, and crystallization differentiation occurred during the evolution of magma. However, Na2O, did not change with SiO2.   MgO and SiO2 of granodiorite in the area were negatively correlated, reflecting the typical characteristics of subducted oceanic crust melting adakite (Figure 10f) [60]. In the Sr/Y-Y diagram (Figure 11a), all samples show a negative correlation trend and fall within the area of adakite with an affinity of melting of oceanic subduction, which is obviously different from that of melting of thickened lower continental crust (Figure 11b).
To sum up, we believe that the Hucunnan granodiorite is adakite that originated from the partial melting of oceanic crust, mixed with mantle-derived materials. After melting, the magma underwent crystallization differentiation of hornblende, pyroxene, sphene, biotite, and other minerals, and the crystallization environment had high oxygen fugacity.

Evolutionary Characteristics of Garnet
Garnet is the main skarn mineral in the Hucunnan deposit. In this study, the main elements of garnet at different locations were measured by electron probe. According to the distance from intrusion, the elements from far to near were zk927-b24, zk907-b59, and zk927-b26, respectively.
As one of the common minerals in skarn-type polymetallic deposits, garnet can be used to indicate the diagenetic and metallogenic process of skarn and serve as a direct indicator of mineralized species and genera [62,63]. In the Hucunnan Cu-Mo deposit, the end-member compositions were andradite (And = 61.91-97.88%), grossularite (Gro = 1.19-34.16%), spessartite (Spe = 0.43-2.11%), pyrope (Pyr = 0.00-1.53%), and almandine (Alm = 0.00-1.74%), which belong to the series of andradite and agrossularite solid solutions. The molecular contents of the andradite are similar to the end components of garnet in skarn-type Cu-Mo deposits in the world (Figure 12a).
As skarn gradually moves away from the contact zone and the influence of alteration, its mineral composition also shows a certain change trend. The contents of andradite and grossularite in the Hucunnan deposit were linearly correlated (Figure 12b), and the composition of garnet gradually changed from andradite to grossularite with the direction close to the intrusion. This indicates that the metasomatism of fluid may have decreased gradually with the change of spatial position. There was a negative correlation between FeO and Al2O3 (Figure 12c), and a slight positive correlation between SiO2 and CaO (Figure 12d). With the increase of Fe content, Al content in garnet decreases obviously, which reflects the variation characteristics of composition in the Ca-Fe-Ca-Al garnet solid-melt system and was formed under nonequilibrium condition [64]. Skarn mineral composition is closely related to intrusive rock composition, surrounding rock composition, depth, oxygen fugacity, temperature, pressure, etc. [66]. Based on comprehensive research on the main mineral composition and symbiotic combination of skarn in skarn deposits, the changes in physical and chemical conditions in the process of diagenesis and mineralization of the deposit can be revealed. Changes in the mineral compositions of garnet have important indicative significance for the formation environment of skarn [63]. The Fe 2+ /Fe 3+ ratios of garnet mineral in the Hucunnan deposit were 0.00-0.04, indicating that it was formed in an oxidized environment. The end members of the minerals are composed of Gro1.19-34.16And61.91-97.88Spe + Pyr + Alm0.62-3.85, which is dominated by andradite, indicating that it was formed under more oxidized conditions [67]. Previous experiments on the physical and chemical conditions for the formation of andradite showed that its main formation temperature was 450-600 °C, and it was in an oxidization-weak alkaline environment with pH 4.0-11.0 [68], indicating that the ore-forming fluids in the skarn stage of the Hucunnan deposit can be characterized by high temperature, low acidity, and high oxygen fugacity.

Mineralization and Material Source
The metallogenic model ages of the Hucunnan Cu-Mo deposit are 138.70 ± 1.87 Ma and 140.03 ± 2.8 Ma, which are consistent with the main metallogenic ages of the Tongling ore-cluster region [19]. The diagenetic age of Hucunnan granodiorite measured by predecessors is 137.5 ± 1.2 Ma [36], indicating that the age of diagenesis and mineralization is basically the same. Therefore, the formation of the Hucunnan Cu-Mo deposit may be closely related to the Hucunnan granodiorite.
Variations of Fe and S contents (Figure 13a) are predominantly related to trace metal substitution in pyrite [69]. The well-defined negative correlation between Cu and Fe in the Hucunnan deposit (Figure 13b) indicates Cu 2+ substitution for Fe 2+ in pyrite [70]. Previous studies suggested that Cu can be structurally bound in pyrite, reaching weight percent levels [70,71], and can also be concentrated in pyrite due to micro-or nanoparticle inclusions of chalcopyrite [72]. We propose that the negative correlation of Cu and Fe is due to Cu 2+ solid solution in pyrite because of the low Cu contents in the studied pyrite samples and pyrite coexisting with chalcopyrite. The elemental compositions of pyrite in skarn type and porphyry type have similar evolutionary trends, indicating that they may have had the same source of ore-forming materials. The sulfur isotopic composition of metal sulfide is widely used to trace the source of sulfur in ore deposits. δ 34 S will vary greatly from one source to another. From mantle sulfur, δ 34 S will be 0 ± 3‰, from sea sulfur, it will be +20‰, and from sedimentary sulfur, it will be negative [73].
The S isotopes show that the δ 34 S values of sulfide were generally less than 10‰ in the Tongling ore-cluster region, mainly concentrated in the 2‰-8‰ range, while the carboniferous strata of pyrite in the δ 34 S values were 0.6‰-29.5‰ [74,75]. The δ 34 S values of the Hucunnan deposit were consistent with other deposits in the area, indicating that most of the sulfur in this deposit did not originate from the strata (Figure 14), but from magma of acidic intrusive rocks in the area. Meanwhile, the δ 34 S value of skarn ore was lower than that of porphyry ore, indicating that the ore-forming hydrothermal underwent further evolution during the intrusion process and more stratigraphic materials were involved in the skarn mineralization stage. Figure 14. Distribution of sulfur isotope (data after [76][77][78][79][80][81][82][83][84][85][86][87][88][89][90][91][92]).
By projecting the isotopic composition of lead in various types of ore onto the diagram of 206 Pb/ 204 Pb-207 Pb/ 204 Pb, and 206 Pb/ 204 Pb-208 Pb/ 204 Pb (Figure 15a), it can be seen that the projection points of the isotopic composition of lead in two types of ore are relatively concentrated, suggesting that they may have had the same origin or evolutionary process [93]. The Pb isotope is located between the upper crust and the mantle. It is inferred that Pb in the Hucunnan deposit is mainly mantlederived, which may be contaminated by crustal material (Figure 15b).  [94]; Geochron, zero isochron line; NHRL, northern hemisphere reference line [95]. Lines of mantle, lower crust, and upper crust of China are after [96].
The 206 Pb/ 204 Pb values of ore lead in the Hucunnan deposit were 18.10-18.54, with an average value of 18.35, which is comparable to the 206 Pb/ 204 Pb value of rock lead in the Yanshan period (18.30) and the value of post-magmatic hydrothermal deposits in the Middle and Lower Yangtze River [97,98], showing that the ore-forming materials were originated from magma. In the diagram of Pb isotope of ore, the ore sample point of Hucunnan deposit also mainly falls in the magmatic area, and a small amount of skarn ore falls in the boundary area between sedimentary and magmatic rock, which shows that ore ore-forming materials have the characteristics of magma source, and some stratigraphic materials were added in the skarn stage.
In conclusion, we believe that the formation of the Hucunnan Cu-Mo deposit was closely related to the intrusion of granodiorite. The ore-forming materials were mainly from magma and partly from strata. The skarn orebody has more stratigraphic materials than the porphyry orebody.

Metallogenic Geological Background
According to previous studies, during the Late Jurassic period (165 ± 5 Ma to 145 Ma), the ancient Pacific plate began to subduct, resulting in the eastern part of China being extruded and uplifted as a whole to form a plateau, suffering from erosion and missing the deposition of the Late Jurassic, and the lithosphere rapidly thickened [99,100]. A fluid or melt formed by the dehydration of the oceanic crust produced metasomatism of the mantle wedge to a low degree of partial melting of the mantle wedge or mixing with mantle magma to form alkaline basaltic magma. Then, it entered the postcollisional orogenic stress transition period (145-130 Ma). The tectonic stress transitioned from compression to tension. The alkaline basaltic magma from the enriched mantle bottom invaded the lower crust, which not only caused the crust to thicken, but also caused partial melting of the crust to form a deep magma chamber, with a part of the deep magma acidic like adakite. The adakite magma mixed with the lower crust source magma, evolved, and transgressed upward, forming intrusive rocks, mainly composed of intermediate acid in the Tongling area. Later, it entered the extensional period (130-80 Ma), during which the lithosphere was desecrated and thinned, the asthenosphere underwent upwelling, and large-scale volcanic eruption occurred in the region [27]. In the tectonic discrimination diagrams (Figure 16), the data of Hucunnan granodiorite are mostly plotted in the field of volcanic arc granite (VAG). Combined with its diagenetic age, it would have been against the dynamic background of the transformation from compression to tension.  [101]. VAG, volcanic arc granite; ORG, oceanic ridge granite; Syn-COLG, syn-collision granite; WPG, intra-plate granite.

Metallogenic Model
According to previous studies, the diagenetic ages of granodiorite in the Tongling ore cluster are mainly from 147-137 Ma, and the metallogenic ages are mainly from 141-137 Ma [19]. The diagenetic age of Hucunnan granodiorite is 137.5 ± 1.2 Ma [36], and the ages of the metallogenic model are 138.70 ± 1.87 Ma and 140.03 ± 2.8 Ma, which are consistent with the overall diagenetic and metallogenic age of the region, indicating that magmatism and mineralization are inseparable. Research on the major elements of biotite and garnet showed that the ore-forming environment of the Hucunnan Cu-Mo deposit was characterized by high oxygen fugacity. This is consistent with the view that large-scale Cu-Mo polymetallic mineralization is closely related to high oxygen fugacity magmas [102][103][104][105][106][107][108][109]. Previous studies on the ore-forming fluid of the deposit showed that the time changes of redox conditions, acid balance, and temperature in the ore-forming fluid led to the temporal separation of copper and molybdenum in Hucunnan skarn deposit [37]. Research on the major elements and S-Pb isotopes of pyrite showed that the ore-forming materials of skarn and porphyry orebody in this deposit have similar sources. They are mainly from magma and are mixed with some stratigraphic materials. Compared with porphyry orebody, skarn orebody has more stratigraphic materials. Combined with the geochemical characteristics of the rocks, the compositional patterns are summarized as follows ( Figure 17).
During 145-135 Ma, the regional tectonic stress field changed to a stress relaxation stage of transition from compression to tension. Before that, the collision and welding of the North China plate and Yangtze plate from north to south and the subsequent subduction northwest to the ancient Pacific plate resulted in lithospheric thickening and melting due to decompression. A fluid or melt formed by the dehydration of the oceanic crust made metasomatism of the mantle wedge to a low degree of partial melting of the mantle wedge or mixed with mantle magma to form alkaline basaltic magma, and the underplating of the lower crust rock formed adakitic magma. The two kinds of magma show mixed coalescence, evolution, and random rise and emplacement along the east-west deep fault. During the upwelling process, crystallization of apatite, sphene, amphibole, pyroxene, and plagioclase and assimilation and contamination of crustal materials took place, forming the Hucunnan granodiorite. Magmatism in the area not only brought ore-forming materials, but also provided a lot of heat energy for groundwater circulation. The ore-forming hydrothermal fluid rose to the shallow intrusion, forming lenticular orebodies and vein and veinlet disseminated ore in the intrusive body and its silicate surrounding rock, developing slightly zonal K-silicified hydrothermal alteration and forming porphyry deposits. Skarnization formed in the contact zone between the rock body and carbonate formation, with strong contact thermal metamorphism. Under the influence of thermal activity and volatilization, the ore-forming materials in the strata were activated and migrated into the ore-forming hydrothermal solution, and were enriched in favorable spaces, such as interlayer fissures and fold ends, to form skarn orebodies, usually occurring along the strata. In addition, a small amount of ore-forming hydrothermal fluids continued to intrude along the fractures, forming a small-scale hydrothermal orebody near the surface, which is the end metallogenic effect of the metallogenic system.

Conclusions
(1) The geochemical affinity of this granodiorite related to Cu-Mo mineralization is adakite, the metallogenic model ages of the Hucunnan Cu-Mo deposit are 138.70 ± 1.87 Ma and 140.03 ± 2.8 Ma, and the ages of diagenesis and mineralization are basically the same. Magma originated in the mantle and subducted oceanic crust and mixed with some crustal material. It formed under the dynamic background of a transformation from compression to extension. (2) The chemical compositions of biotites show that the magma formed in a high oxygen fugacity environment, and garnets show that the ore-forming fluids in the skarn stage are characterized by high temperature, low acidity, and high oxygen fugacity. (3) The formation of the Hucunnan Cu-Mo deposit was closely related to the intrusion of granodiorite. The ore-forming materials were mainly originated from magma and partly from strata. The skarn orebody has more stratigraphic materials than the porphyry orebody. The porphyry/skarn metallogenic system was formed under the joint action of stratum structure magmatic rocks.

Supplementary Materials:
The following are available online at www.mdpi.com/2075-163X/10/2/107/s1, Table  S1: Major (wt.%) and trace element (ppm) results of Hucunnan granodiorites; Table S2: Chemical compositions of biotites in the Hucunnan deposit; Table S3: Chemical composition of garnets in the Hucunnan deposit; Table  S4: Chemical composition of pyrites in the Hucunnan deposit; Table S5: S isotopes of ores in the Hucunnan deposit; Table S6: Pb isotopes of ores in the Hucunnan deposit; Table S7: Re-Os isotopes of molybdenites in the Hucunnan deposit.
Author Contributions: K.S. wrote the paper; X.Y., J.D., and J.C. designed the experiments; Q.W. and Y.C. took part in the field investigation. All authors have read and agreed to the published version of the manuscript.