Mineral Chemistry, S-Pb-O Isotopes, and S / Se Ratios of the Niubiziliang Ni-(Cu) Sulﬁde Deposit in North Qaidam Orogenic Belt, NW China: Constraints on the Parental Magma Composition, Evolution, and Sulfur Saturation Mechanism

: The Niubiziliang Ni-(Cu) deposit is the ﬁrst magmatic Ni-Cu sulﬁde deposit in the North Qaidam Orogenic Belt (NQOB), NW China, and plays a signiﬁcant role in geological evolution, Ni-Cu mineralization, and exploration in the NQOB. Here, we report on the mineral chemistry, S-Pb-O isotopes, and S / Se ratios of the maﬁc-ultramaﬁc complex, which provide insights on the parental magma, evolution, and sulfur saturation mechanism. The Niubiziliang maﬁc-ultramaﬁc intrusion contains four ore blocks and about ten Ni-(Cu) ore / mineralization bodies. Olivines in Niubiziliang belong to the species of chrysolite with Fo values of 88~89, and the pyroxenes are mainly orthopyroxene (En = 79~82) and clinopyroxene (En = 44~40). The olivines and some pyroxenes likely crystallized in a magma chamber at a depth of 35.45~36.55 km at a high temperature (1289~1369 ◦ C) and pressure (9.38~9.67 kbar), whereas the Niubiziliang complex formed at a moderate depth (8.13~8.70 km) with a temperature and pressure of 1159~1253 ◦ C and 2.15~2.30 kbar, respectively. The parental magma was considered to be high-Mg picritic basalt with MgO and NiO contents of 14.95~16.58% and 0.053~0.068%, respectively, which indicated high-degree partial melting of the depleted mantle. The mantle-derived primary magma underwent signiﬁcant fractional crystallization and crustal assimilation and contamination, which was strongly supported by S-Pb-O isotope data and S / Se ratios, resulting in sulfur saturation and sulﬁde immiscibility in the magma. Crustal assimilation and contamination contributed more to sulfur saturation than fractional crystallization. paper attempts to determine the parental magma composition, crystallization conditions, magmatic evolution, sulfur saturation mechanism, and their constraints on mineralization through a systematic study of the mineralogy, O-S-Pb isotope geochemistry, and S / Se ratios in this complex, in order to provide new insights into the metallogenesis of the Niubiziliang Ni-(Cu) sulﬁde deposit.

. Structural location map (a) and geological sketch map (b) of the northern margin of Qaidam Basin (modified after Qian et al. (2015) [15]).

Mining Area Geology
The strata in the Niubiziliang Ni-(Cu) mining area are dominated by the Paleoproterozoic Jinshuikou Group and Quaternary sedimentary successions ( Figure 2). The Jinshuiokou Group is a set of gneiss with intense migmatitization, mainly including biotite plagiogneiss, amphibolite gneiss, and tremolite marble [11]. The Quaternary sediments are mainly composed of sand, gravel, and sand oil.  [15]).
The strata in the study area mainly include Paleoproterozoic Jinshuikou Group, Jurassic Dameigou Formation, Neogene Youshashan Formation, and Quaternary sediment strata. The Jinshuikou Group is the main stratum in the study area, which can be divided into gneiss formation and marble formation. The Jinshuiou Group in the Niubiziiliang Ni-(Cu) mining area is mainly gneiss formation with NW-SE strike, and the lithology is dominated by Biotite plagiogneiss with some sericite-quartz schist, muscovite-bearing quartzite, and tremolite-bearing marble. The Altun south margin fault, north Qaidam margin fault, and Hongsanhan-Niubiziliang fault formed the regional structural framework of the northwestern margin of Qaidam Massif. At the same time, many ductile shear zones, folds, and joints were formed, which had an important influence on the tectonism and magmatism in this area. The faults can be divided into four groups according to the strike, including a nearly SE direction, NE-SW direction, NW-SE direction, and NS direction. The NE-SW faults are related to the formation of the AOB, while the NW-SE faults are likely related to the formation and evolution of the NQOB. The magmatic rocks in this area are widely distributed, large in scale, and mainly intrusive rocks, with poor volcanic rocks. The intrusive rocks were developed from ultrabasic to acidic with many types of rocks and were mainly formed in the Luliang, Jinning, Hercynian, and Indosinian periods, among which Hercynian and Indosinian acid intrusive rocks represent the most developed. Basic-ultrabasic rocks, including Niubiziliang, Yanchangbeishan, Qaidam gate, western Qaidam gate, and western Qingxinjie Mt. rocks, have mainly been produced as complex or vein-like forms, and are distributed sporadically (Figure 1b). The Cu-Ni mineralization in the study area is almost entirely related to the mafic-ultramafic complex. The most Hercynian diorites, granites, and other intrusive rocks intruded at the strong extensional background after the closure of the northern margin of the Qaidam Massif Ocean (NMQMO), which is consistent with that of the mafic-ultramafic rocks in Niubiziliang [18,19].

Mining Area Geology
The strata in the Niubiziliang Ni-(Cu) mining area are dominated by the Paleoproterozoic Jinshuikou Group and Quaternary sedimentary successions ( Figure 2). The Jinshuiokou Group is a set of gneiss with intense migmatitization, mainly including biotite plagiogneiss, amphibolite gneiss, and tremolite marble [11]. The Quaternary sediments are mainly composed of sand, gravel, and sand oil.
The faults in the Niubiziliang mining area were all formed during the post-mineralization stage and some of them had a destructive effect on the deposit. According to the orientation, they can be divided into two groups, as follows: The NE-SW and NW-SE trending fault systems. Among these, the F 1 fault has the largest scale and almost runs through the whole mining area. This fault strikes the NE and dips to the NW, whereas the middle of F 1 dips to the SE. This gives the interpretation that the F 1 is a sinistral translational fault with normal fault characteristics [20]. It divides the original ore-hosted ultramafic intrusion into the No. II and III ore blocks ( Figure 2). The eastern part of F 1 is covered by the Quaternary sediments, forming a large gully. The normal fault effect of F 1 relatively lowers the No. II ore block in the hanging wall and raises the No. III ore block in the footwall. Due to this, under the same denudation conditions after mineralization, the ore-hosted ultramafic intrusion in the No. III ore block was exposed on the surface and the bottom ore body in the No. II ore block was preserved.
The Niubiziliang mafic-ultramafic complex, which formed in the Devonian period (367 Ma [13], 388~402 Ma [15], and 373~389 Ma [16]), is the metallogenic intrusion of the Niubiziliang Ni-(Cu) sulfide deposit. The rocks from the peridotite facies to the gabbro facies are well-developed in this mafic-ultramafic complex, and the harzburgite is the most important ore-hosting rock. Devonian magmatic rocks, mainly including diorite and granite, are the products of magmatic activities after mineralization and have a great influence on the preservation of the deposit [19,21]. Dark inclusions are common in the diorite, the condensation edge of which develops crystalline rings. The monzonitic granite is widely exposed in the central and northeastern part of the mining area and intruded into the mafic-ultramafic complex, having a very destructive effect on the No. I ore block.
Minerals 2020, 10, x FOR PEER REVIEW 5 of 27 complex, and the harzburgite is the most important ore-hosting rock. Devonian magmatic rocks, mainly including diorite and granite, are the products of magmatic activities after mineralization and have a great influence on the preservation of the deposit [19,21]. Dark inclusions are common in the diorite, the condensation edge of which develops crystalline rings. The monzonitic granite is widely exposed in the central and northeastern part of the mining area and intruded into the mafic-ultramafic complex, having a very destructive effect on the No. I ore block.

Mafic-Ultramafic Intrusions
The Niubiziliang Ni-(Cu) deposit can be divided into four ore blocks according to the location of ultramafic intrusion (Figures 2 and 3), among which the Ni-(Cu) ore bodies are mainly hosted in the No. II and III ultramafic intrusions [20]. Gabbros in this deposit are widely distributed, forming the main part of the complex, and the ultramafic intrusions output in stocks into the gabbro or Jinshuikou Group gneiss. The exposed areas of ultramafic intrusions in the four ore blocks are different and the overall areas are small, mostly ranging between 0.08 and 0.12 km 2 . There are a few banded cumulate gabbros in the No. I ore block and the No. IV ore block is composed of widely distributed gabbros and three small separated outcrops of ultramafic rocks. The ore-hosting ultramafic intrusions in No. II and III ore blocks originally belonged to the same intrusion, which was divided into two isolated intrusions by the F1 fault. The ultramafic rocks in No. I and IV ore blocks are mainly pyroxenite, websterite, and olivine pyroxenite, whereas the ultramafic rocks in the No. II and III ore blocks are dominated by harzburgite and lherzolite ( Figure 3).

Mafic-Ultramafic Intrusions
The Niubiziliang Ni-(Cu) deposit can be divided into four ore blocks according to the location of ultramafic intrusion (Figures 2 and 3), among which the Ni-(Cu) ore bodies are mainly hosted in the No. II and III ultramafic intrusions [20]. Gabbros in this deposit are widely distributed, forming the main part of the complex, and the ultramafic intrusions output in stocks into the gabbro or Jinshuikou Group gneiss. The exposed areas of ultramafic intrusions in the four ore blocks are different and the overall areas are small, mostly ranging between 0.08 and 0.12 km 2 . There are a few banded cumulate gabbros in the No. I ore block and the No. IV ore block is composed of widely distributed gabbros and three small separated outcrops of ultramafic rocks. The ore-hosting ultramafic intrusions in No. II and III ore blocks originally belonged to the same intrusion, which was divided into two isolated intrusions by the F 1 fault. The ultramafic rocks in No. I and IV ore blocks are mainly pyroxenite, websterite, and olivine pyroxenite, whereas the ultramafic rocks in the No. II and III ore blocks are dominated by harzburgite and lherzolite ( Figure 3).  The Niubiziliang ore-forming pluton is a complex which was formed by the multiple emplacement of mafic-ultramafic magmas in the Devonian period [16], and was then intruded by post-metallogenic diorite and granite as dyke or stock. It is inferred that there are two possible shapes of the deep ultramafic intrusion based on a large number of drill holes and magnetic data, as follows: (a) Intruding the gabbros as stock, such as the No. I and IV ore blocks. The ultramafic intrusion in the No. IV ore block is a semi-steep lopolith in space, and the ultramafic rocks on the surface are the outcrops of the lopolith edge. The second is (b) intruding the gabbros or gneiss as a steeply dipping and gentle-lateral-trending apophysis, such as the No. II and III ore blocks. They are controlled by the NW trending fault with striking NW and dipping SW ( Figure  4). The shape of the ultramafic intrusion has a strong effect on the location of the ore body. The steeply dipping apophysis usually forms the hanging orebody, while the semi-steep lopolith tends to form the bottom ore body. The Niubiziliang ore-forming pluton is a complex which was formed by the multiple emplacement of mafic-ultramafic magmas in the Devonian period [16], and was then intruded by post-metallogenic diorite and granite as dyke or stock. It is inferred that there are two possible shapes of the deep ultramafic intrusion based on a large number of drill holes and magnetic data, as follows: (a) Intruding the gabbros as stock, such as the No. I and IV ore blocks. The ultramafic intrusion in the No. IV ore block is a semi-steep lopolith in space, and the ultramafic rocks on the surface are the outcrops of the lopolith edge. The second is (b) intruding the gabbros or gneiss as a steeply dipping and gentle-lateral-trending apophysis, such as the No. II and III ore blocks. They are controlled by the NW trending fault with striking NW and dipping SW ( Figure 4). The shape of the ultramafic intrusion has a strong effect on the location of the ore body. The steeply dipping apophysis usually forms the hanging orebody, while the semi-steep lopolith tends to form the bottom ore body.

Characteristics of Ni-(Cu) Ore Body
The mineralization in the four ore blocks is quite different. In general, the gabbro is rarely mineralized, with starspot pyrite occurring locally, whereas the ultramafic rocks mainly include starspot, disseminated, sideronitic, and crumby/massive sulfides. Most ultramafic intrusions contain starspot and disseminated sulfides, but their grades of Cu, Ni, and Co are generally low, with no economic value. The amount of sulfide is closely related to the content of olivine, for example, the harzburgite and lherzolite in the No. II and III ore blocks contain more sulfides than pyroxenite and websterite in the No. I and IV ore blocks. The ore body is usually formed by ultramafic rocks with sideronitic or crumby/massive sulfides, with a generally high grade of Cu, Ni, and Co. About 10 ore bodies (including concealed ore bodies) are defined in the No. II and III ore blocks at present, and their characteristics are listed in Table 1. However, there are no ore bodies in the No. I and IV ore blocks, in which ultramafic rocks only contain starspot or disseminated sulfides.

Characteristics of Ni-(Cu) Ore Body
The mineralization in the four ore blocks is quite different. In general, the gabbro is rarely mineralized, with starspot pyrite occurring locally, whereas the ultramafic rocks mainly include starspot, disseminated, sideronitic, and crumby/massive sulfides. Most ultramafic intrusions contain starspot and disseminated sulfides, but their grades of Cu, Ni, and Co are generally low, with no economic value. The amount of sulfide is closely related to the content of olivine, for example, the harzburgite and lherzolite in the No. II and III ore blocks contain more sulfides than pyroxenite and websterite in the No. I and IV ore blocks. The ore body is usually formed by ultramafic rocks with sideronitic or crumby/massive sulfides, with a generally high grade of Cu, Ni, and Co. About 10 ore bodies (including concealed ore bodies) are defined in the No. II and III ore blocks at present, and their characteristics are listed in Table 1. However, there are no ore bodies in the No. I and IV ore blocks, in which ultramafic rocks only contain starspot or disseminated sulfides.
Five ore/mineralization Ni-(Cu) bodies can be found in the No. II ore block, including one Ni-Cu industrial ore body (G-M1), two low-grade ore bodies (D-M1 and D-M2), and two mineralization bodies (M3 and M4). The M1~M3 ore bodies occur in the cracked harzburgite and lherzolite in the F 3 fracture zone dipping to S at 72 • , in which the hanging wall is ore-bearing harzburgite and the footwall is Jinshuikou Group gneiss ( Figure 4). The other five ore/mineralization bodies have the same ore rock, sulfide assemblage, and texture. They are a group of ore bodies of the same type, and their ore rocks are all harzburgite and lherzolite. The G-M1 ore body contains massive ores, whereas the others mainly contain disseminated sulfides. Prryhotite, pentlandite, and chalcopyrite are the main sulfides of the ore bodies, and the more sulfides, the higher the ore grade.
Seven Ni-(Cu) ore/mineralization bodies have been found in the No. III ore block, all of which are concealed ore bodies, including two industrial ore bodies (M7 and G-M10), four low-grade ore bodies (M6, M8, M9, and D-M10), and one mineralization body (M5). The rock core of ZKIII-0801 and ZKIII-1601 in the G-M1 ore body contains many massive ores and has a high content of sulfides, with a maximum value of about 80%.

Sulfides and Ni-(Cu) Mineralization
Pyrrhotite, pentlandite, chalcopyrite, and pyrite are the main sulfides of the Niubiziliang Ni-(Cu) deposit, with minor amounts of violarite, bornite, magnetite, and ilmenite ( Figure 5), and there are a few supergene minerals (e.g., malachite, annabergite, and limonite) in the surface. They are often disseminated, while sideronitic and massive sulfides occur locally. Pentlandite usually occurs in the inner or outer margin of pyrrhotite in the form of xenomorphic aggregates, showing a typical exsolution texture and sulfide assemblage in the magmatic Ni-Cu sulfide deposit.
Pyrite is the most widely distributed sulfide in Niubiziliang, with a variable content of 5~20%. It generally has a xenomorphic-hypautomorphic granularity, is vein-like, and is crumby in shape, with a size of 0.05~2.50 mm (Figure 5a). The vein-like pyrites formed in the late stage often cut the early rocks and sulfides.
The pyrrhotite is distributed between the gangue minerals in the form of a xenomorphic or irregular granularity, with a size of 0.02~1.00 mm. It generally coexists with pentlandite and chalcopyrite, and the contact boundary between them is smooth and straight, showing a common symbiotic texture (Figure 5b-d). means that its Ni grade is higher than that of other ore/mineralization bodies. A crumby and massive ore usually contains >60% sulfides and is mainly hosted in the GM-1 Ni-Cu ore body, which provides the majority of Ni resources of Niubiziliang, with the highest grade of Ni of 1.57%. The GM-1 Ni-Cu ore body contains many pyrites, in addition to the Po + Pn + Ccp assemblage, and chalcopyite is often associated with pyrite (Figure 5a), which is quite different from other ore/mineralization bodies.

Samples and Petrography
The Niubiziliang mafic-ultramafic complex can be divided into mafic to ultramafic rocks, and is comprised of peridotites, pyroxenites, and gabbros. The gabbro is widely distributed, constituting the main part of the complex, and is then intruded by the peridotite, pyroxenite, and postmineralization diorite and granite. The probable emplaced sequence of the mafic-ultramafic rocks is gabbro→pyroxenite→peridotite, according to a large number of rock core data. Harzburgite is the most important ore-hosted rock, while the mineralization of gabbro is weak, with a few pyrites occurring locally. The mafic-ultramafic rocks in Niubiziliang mainly contain harzburgite, olivine pyroxenite, lherzolite, olivine websterite, hornblende lherzolite, pyroxenite, websterite, olivine gabbro, melagabbro, bojite, and gabbro ( Figure 6).
The gabbro is gray in color, has a gabbro texture, and is massive, and mainly contains clinopyroxene (40%), plagioclase (55%), and hornblende (5%) (Figure 6a). The plagioclase is euhedral clintheriform, measuring 1.0~1.5 mm, with a weak alteration. The pyroxene is mainly filled between Pentlandite is the most important Ni-bearing mineral in Niubiziliang and provides the majority of Ni resources. It usually occurs in the inner or outer margin of pyrrhotite in the form of xenomorphic-hypidiomorphic aggregates with a minor hypidiomorphic granularity in gangue minerals, with a size of about 0.1~0.5 mm (Figure 5b-d).
Chalcopyrite in Niubiziliang mainly occurs in the following two forms: (a) Coexisting with pyrite ( Figure 5a), and (b) coexisting with pentlandite and pyrrhotite (Figure 5b-d). The former chalcopyrite is mainly hosted in the GM-1 Ni-Cu ore body in disseminated, crumby, or massive forms, whereas the latter is mainly distributed in the remaining Ni-(Cu) ore/mineralization bodies in dissemination, dropwise, or veinlet forms.
The original sulfide ore is the most important ore in the Niubiziliang Ni-(Cu) deposit, with a few oxidized ores. It is generally composed of pyrite, pyrrhotite, chalcopyrite, and pentlandite, and can be divided into disseminated (<20% sulfides), sideronitic (20~40% sulfides), vein-like (10~30% sulfides), crumby (40~60% sulfides), and massive (>60% sulfides) types. The ore texture is mainly a xenomorphic granular texture, xenomorphic intersertal texture, and sideronitic texture, and the disseminated structure, vein structure, and massive structure are the main ore structures in the Niubiziliang deposit. In fact, the formation of different types of ore is closely related to the magma exhibiting variable sulfide melting.
A disseminated ore is usually composed of <20% sulfides, being the main type of ore in Niubiziliang. The sulfides in the disseminated ore are mainly pyrrhotite, pentlandite, and chalcopyrite, with a few pyrites and magnetites (Figure 5c), in which pyrrhotite, pentlandite, and chalcopyrite often coexist.
A disseminated ore becomes a sideronitic ore when the sulfide content reaches 20~40% and sulfides are often filled between the silicate minerals in the form of a xenomorphic assemblage. The sideronitic ore mainly occurs in the M10 Ni-(Cu) ore body, which means that its Ni grade is higher than that of other ore/mineralization bodies.
A crumby and massive ore usually contains >60% sulfides and is mainly hosted in the GM-1 Ni-Cu ore body, which provides the majority of Ni resources of Niubiziliang, with the highest grade of Ni of 1.57%. The GM-1 Ni-Cu ore body contains many pyrites, in addition to the Po + Pn + Ccp assemblage, and chalcopyite is often associated with pyrite ( Figure 5a), which is quite different from other ore/mineralization bodies.

Samples and Petrography
The Niubiziliang mafic-ultramafic complex can be divided into mafic to ultramafic rocks, and is comprised of peridotites, pyroxenites, and gabbros. The gabbro is widely distributed, constituting the main part of the complex, and is then intruded by the peridotite, pyroxenite, and post-mineralization diorite and granite. The probable emplaced sequence of the mafic-ultramafic rocks is gabbro→pyroxenite→peridotite, according to a large number of rock core data. Harzburgite is the most important ore-hosted rock, while the mineralization of gabbro is weak, with a few pyrites occurring locally. The mafic-ultramafic rocks in Niubiziliang mainly contain harzburgite, olivine pyroxenite, lherzolite, olivine websterite, hornblende lherzolite, pyroxenite, websterite, olivine gabbro, melagabbro, bojite, and gabbro ( Figure 6). the plagioclases in the form of a xenomorphic-hypidiomorphic granularity, with a size of about 1.0 mm. In addition, there are a few banded gabbros in the No. I ore block, and the pyroxenes and plagiaclases occur in band and alternation, which indicates the in situ fractional crystallization and cumulation of mafic minerals in the diagenic process.
The olivine pyroxenite is black in color, has a granular texture, and is massive, and is mainly composed of pyroxene (50%), olivine (40%), hornblende (5%), and plagioclase (5%) (Figure 6d). The pyroxenes are mainly composed of orthopyroxene (70%) and clinopyroxene (30%), and they have a xenomorphic-hypidiomorphic granularity, with a size of 0.5~1.5 mm. The olivine is rotund in shape and 0.3~1.0 mm in size, with cleavage and high protuberance, and has been serpentinized and talcolized. The hornblende is brownish in color, 0.5~2.0 mm in size, and filled between the pyroxene and olivine grains in the form of xenomorphic granularity.
The harzburgite is black in color, granular, and massive, and is mainly composed of olivine (60%), orthopyroxene (30%), clinopyroxene (5%), and opaque mineral (5%), with strong alteration (Figure 6e,f). The olivine is rotund in shape and 1~2 mm in size, with cleavage and high protuberance, and has been serpentinized and talcolized. The orthopyroxene has a xenomorphic-hypidiomorphic granularity, with a size of 1.5~3.0 mm, and most of it is altered into tremolite and talc. The opaque minerals are mainly pyrrhotite, chalcopyrite, and pentlandite.  The gabbro is gray in color, has a gabbro texture, and is massive, and mainly contains clinopyroxene (40%), plagioclase (55%), and hornblende (5%) (Figure 6a). The plagioclase is euhedral clintheriform, measuring 1.0~1.5 mm, with a weak alteration. The pyroxene is mainly filled between the plagioclases in the form of a xenomorphic-hypidiomorphic granularity, with a size of about 1.0 mm. In addition, there are a few banded gabbros in the No. I ore block, and the pyroxenes and plagiaclases occur in band and alternation, which indicates the in situ fractional crystallization and cumulation of mafic minerals in the diagenic process.
The olivine pyroxenite is black in color, has a granular texture, and is massive, and is mainly composed of pyroxene (50%), olivine (40%), hornblende (5%), and plagioclase (5%) (Figure 6d). The pyroxenes are mainly composed of orthopyroxene (70%) and clinopyroxene (30%), and they have a xenomorphic-hypidiomorphic granularity, with a size of 0.5~1.5 mm. The olivine is rotund in shape and 0.3~1.0 mm in size, with cleavage and high protuberance, and has been serpentinized and talcolized. The hornblende is brownish in color, 0.5~2.0 mm in size, and filled between the pyroxene and olivine grains in the form of xenomorphic granularity.
The harzburgite is black in color, granular, and massive, and is mainly composed of olivine (60%), orthopyroxene (30%), clinopyroxene (5%), and opaque mineral (5%), with strong alteration (Figure 6e,f). The olivine is rotund in shape and 1~2 mm in size, with cleavage and high protuberance, and has been serpentinized and talcolized. The orthopyroxene has a xenomorphic-hypidiomorphic granularity, with a size of 1.5~3.0 mm, and most of it is altered into tremolite and talc. The opaque minerals are mainly pyrrhotite, chalcopyrite, and pentlandite.

Electron Microprobe Analysis
A total of eight samples (e.g., harzburgite, websterite, olivine pyroxenite, olivine gabbro, gabbro, and bojite) were selected to conduct the electron microprobe analysis for the main silicate minerals in this study. Eight polished thin sections were analyzed on a JEOL JXA8230 electron probe microanalyzer at the MLR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences. The operating conditions were executed with an accelerating voltage of 20 kV and a beam current of 20 nA, with a 5 µm diameter beam size. The spectral lines, standards, and count times that were used and averages of detection limits are listed in Supplementary Materials Table S1, and the EMPA data of olivines, pyroxenes, plagioclases, and hornblendes presented in this paper are all shown in Tables S2-S5.

Sulfide S-Pb Isotope Analysis
A total of twelve sulfides (i.e., pyrite, pentlandite, and pyrrhotite) were selected from the massive and disseminated ores in the No. II and III ore blocks for S and Pb isotope analyses, and were tested at the Beijing Research Institute of Uranium Geology. A Finnigan MAT 251 mass spectrometer was utilized to analyze the S isotopic compositions and the uncertainty was ±0.2% for δ 34 S V-CDT . The reported data are relative to Vienna Canon Diablo Troilite (V-CDT) sulfide. As for the Pb-isotopic compositions of sulfides, 50 mg of powder was dissolved in a mixed solution of hydrofluoric and perchloric acid, followed by elution through anion exchange resin, in order to separate Pb. The Pb isotopic compositions were measured on an ISOPROBE-T Thermal Ionization Mass Spectrometer with the NBS SRM 981 standard. The precisions for 204 Pb/ 206 Pb and 208 Pb/ 206 Pb ratios were better than 0.005, and the results are listed in Tables 2 and 3.

Whole-Rock O Isotope Analysis
Four mafic-ultramafic rock samples (i.e., harzburgite and gabbro) were obtained from the No. II and III ore blocks to conduct the O isotopic analysis, and were tested at the Beijing Research Institute of Uranium Geology. The whole-rock samples were reacted with pure bromine pentafluoride at 500~680 • C for 14 h, in order to release the O 2 and impurity under a vacuum condition of 10 −3 Pa in the sample preparation unit. After this, the impurity components, such as SiF 4 and BrF 3 , were separated using the freezing method, the pure O 2 was reacted with graphite at a constant temperature under the condition of 700 • C, and a platinum catalyst was employed to produce CO 2 . CO 2 was collected by using the freezing method, and the O isotopic composition of the whole-rock samples was analyzed by MAT253 gas isotopic mass spectrometry. The measurement results are defined as SMOW and marked as δ 18 O V-SMOW , and the analysis accuracy is better than ±0.2% . The reference standard for an oxygen isotope is GBW-04409 and GBW-04410 of quartz, of which δ 18 O is 11.11 ± 0.06% and −1.75 ± 0.08% , respectively, and the results are listed in Table 4. Olivine-pyroxene hornblendite Disseminated sulfides 6.5

Pyroxene
Pyroxene mainly includes orthopyroxene and clinopyroxene (Table S3). The orthopyroxene is mainly distributed in the ultramafic rocks, being an important and symbol mineral of the ore-hosting ultramafic rocks. The En value of the orthopyroxene in Niubiziliang is 79~82, which belongs to clinoenstatite (Figure 7a), and is close to that of Xiarihamu (En = 78~86 [24]).

Pyroxene
Pyroxene mainly includes orthopyroxene and clinopyroxene (Table S3). The orthopyroxene is mainly distributed in the ultramafic rocks, being an important and symbol mineral of the ore-hosting ultramafic rocks. The En value of the orthopyroxene in Niubiziliang is 79~82, which belongs to clinoenstatite (Figure 7a), and is close to that of Xiarihamu (En = 78~86 [24]).

Plagioclase
Plagioclase is widely distributed in all kinds of mafic-ultramafic rocks, and mainly falls into the labradorite, with some falling into bytownite and oligoclase, in the An-Ab-Or diagram (Figure 7c and Table S4). The An values for plagiaclase in ultramafic rocks are usually large, and are dominated by bytownite and labradorite. However, the An values for plagioclase in gabbro vary greatly, and are mainly concentrated at 50~60 (labradorite), with the minority being as low as 10~14 (oligoclase). The An values for plagioclase in ultramafic rocks are generally higher than those of mafic rocks, indicating that the An values of plagioclase are closely related to the degree of basicity of rocks.
The clinopyroxene mainly occurs in gabbro, bojite, and other mafic rocks often coexisting with hornblende, and the websterite and olivine pyroxenite also contain some clinopyroxenes. The clinopyroxenes mainly fall into the augite and diopside in the Wo-En-Fs diagram (Figure 7a) and fall into the tholeiite in the Al 2 O 3 -Na 2 O-TiO 2 diagram (Figure 7b).

Plagioclase
Plagioclase is widely distributed in all kinds of mafic-ultramafic rocks, and mainly falls into the labradorite, with some falling into bytownite and oligoclase, in the An-Ab-Or diagram (Figure 7c and Table S4). The An values for plagiaclase in ultramafic rocks are usually large, and are dominated by bytownite and labradorite. However, the An values for plagioclase in gabbro vary greatly, and are mainly concentrated at 50~60 (labradorite), with the minority being as low as 10~14 (oligoclase). The An values for plagioclase in ultramafic rocks are generally higher than those of mafic rocks, indicating that the An values of plagioclase are closely related to the degree of basicity of rocks.

Hornblende
The hornblende in the Niubiziliang complex mainly belongs to calcareous hornblende (Table S5). In the classification diagram of Leake et al. (1997) [30], the hornblende of olivine pyroxenite was pargasite, whereas the hornblende of basic rocks (e.g., gabbro, bojite, and melagabbro) was pargasite and tremolite ( Figure 8). Hornblendes in olivine pyroxenite all have a mantle origin, while hornblendes from the basic rocks mainly have a crust origin, with some mantle origin, in the diagrams of Al 2 O 3 -TiO 2 and Si-Al atoms (Figure 9). This could indicate that crustal contamination and the ultrabasic magma and mantle-origin hornblende evolved to the basic magma and crust-origin hornblende due to the addition of crustal felsic components.

S Isotopes
The δ 34 S value of the fifteen sulfides in Niubiziliang is −8.4~9.47‰, with a wide range of variation (  (Table 3), with a wide range of variation. The discriminant diagrams show that the Pb in Niubiziliang is mainly derived from the orogenic belt,

S Isotopes
The δ 34 S value of the fifteen sulfides in Niubiziliang is −8.4~9.47‰, with a wide range of variation (  (Table 3), with a wide range of variation. The

S Isotopes
The δ 34 S value of the fifteen sulfides in Niubiziliang is −8.4~9.47% , with a wide range of variation (Table 2), among which eight samples are −1.8~1.4% , and thus close to that of the upper mantle (0 ± 2% [32]). However, the δ 34 S values of the remaining seven samples are −8.4~−2.8% and 3.7~9.47% , which are significantly different from that of the upper mantle.

Pb Isotopes
The values of 206 (Table 3), with a wide range of variation. The discriminant diagrams show that the Pb in Niubiziliang is mainly derived from the orogenic belt, with some upper crust Pb (Figure 10). The Pb isotope of pentlandite in the No. II ore block (falling outside of Figure 10) is quite different from that of the sulfides in the No. III ore block and is abnormal Pb, which could be related to the later hydrothermal activities.

O Isotopes
The δ 18 O value of the four samples employed in this study is 6.7~8.1‰, with a wide range of variation (Table 4), which is close to the data obtained by   [13] (δ 18 O = 5.1~7.8‰). Therefore, all of the samples have higher values than the δ 18 O value of the mantle, except for one sample, which is close to the mantle.

S/Se Ratios
The S/Se value of five samples from the Niubiziliang complex is 1864~5890, most of which are higher than that of the mantle (Table 5). One sample is in the mantle region, another is in the S loss region, and those remaining are in the S addition region (Figure 11).

O Isotopes
The δ 18 O value of the four samples employed in this study is 6.7~8.1% , with a wide range of variation (Table 4), which is close to the data obtained by   [13] (δ 18 O = 5.1~7.8% ). Therefore, all of the samples have higher values than the δ 18 O value of the mantle, except for one sample, which is close to the mantle.

S/Se Ratios
The S/Se value of five samples from the Niubiziliang complex is 1864~5890, most of which are higher than that of the mantle (Table 5). One sample is in the mantle region, another is in the S loss region, and those remaining are in the S addition region (Figure 11). sample, which is close to the mantle.

S/Se Ratios
The S/Se value of five samples from the Niubiziliang complex is 1864~5890, most of which are higher than that of the mantle (Table 5). One sample is in the mantle region, another is in the S loss region, and those remaining are in the S addition region (Figure 11).

Parental Magma
Previous studies have shown that the distribution coefficient of Mg-Fe between peridotite and melt was a relatively stable value, that is, Kd Ol-Melt = (TFeO/MgO) Ol /(TFeO/MgO mag ) = 0.3~0.33 [35,36]. The MgO/TFeO value of the melt coexisting with olivine in this study was obtained by using the above formula, resulting in a value of 1.37~1.51 (choosing Fo max = 89 for olivine), which was lower than that of the whole-rock geochemistry (MgO/TFeO = 2.36~2.98); selecting the whole-rock geochemistry data from Yu et al. (2019) [16] in the Niubiziliang complex. Therefore, the analyzed results in this study could not represent the coexisting melt components, and there were excess olivines adding to the pre-emplaced magma, implying that the pre-emplaced magma is a derivative magma that has undergone deep olivine crystallization by the parental magma. Meanwhile, the NiO content of olivine (0.10~0.23%) in this study was significantly lower than that of olivine (≈0.4% [24]) crystallized from the primitive basaltic magma, which equilibrated with pyrolite. This indicates that the magma of the Niubiziliang complex is derivative magma from which the parental magma underwent differentiation and resulted in strong NiO depletion. Additionally, all of the samples are located below the curve representing the parental magma composition in the MgO-Fo-FeO T diagram (Figure 12), which also indicates the addition of extra olivine. 6. Discussion

Parental Magma
Previous studies have shown that the distribution coefficient of Mg-Fe between peridotite and melt was a relatively stable value, that is, KdOl-Melt = (TFeO/MgO)Ol/(TFeO/MgOmag) = 0.3~0.33 [35,36]. The MgO/TFeO value of the melt coexisting with olivine in this study was obtained by using the above formula, resulting in a value of 1.37~1.51 (choosing Fomax = 89 for olivine), which was lower than that of the whole-rock geochemistry (MgO/TFeO = 2.36~2.98); selecting the whole-rock geochemistry data from Yu et al. (2019) [16] in the Niubiziliang complex. Therefore, the analyzed results in this study could not represent the coexisting melt components, and there were excess olivines adding to the pre-emplaced magma, implying that the pre-emplaced magma is a derivative magma that has undergone deep olivine crystallization by the parental magma. Meanwhile, the NiO content of olivine (0.10~0.23%) in this study was significantly lower than that of olivine (≈0.4% [24]) crystallized from the primitive basaltic magma, which equilibrated with pyrolite. This indicates that the magma of the Niubiziliang complex is derivative magma from which the parental magma underwent differentiation and resulted in strong NiO depletion. Additionally, all of the samples are located below the curve representing the parental magma composition in the MgO-Fo-FeO T diagram ( Figure 12), which also indicates the addition of extra olivine.  [35]). Meanwhile, the MgO and FeO T contents of the parental magma could also be estimated to be 14.9% and 13.43%, respectively, by using the extrapolation method of MgO-FeO T (Figure 13a). Therefore, we believe that the MgO content of the parental magma is about 14.95~16.58%, which is higher than the value (MgO = 10.8%) estimated previously by Ling et al. (2014a) [13], and the corresponding FeO T content is about 13.43~14.89%. The difference mainly reflects the different Fo value of olivine in the samples used in  [35]). Meanwhile, the MgO and FeO T contents of the parental magma could also be estimated to be 14.9% and 13.43%, respectively, by using the extrapolation method of MgO-FeO T (Figure 13a). Therefore, we believe that the MgO content of the parental magma is about 14.95~16.58%, which is higher than the value (MgO = 10.8%) estimated previously by Ling et al. (2014a) [13], and the corresponding FeO T content is about 13.43~14.89%. The difference mainly reflects the different Fo value of olivine in the samples used in the calculations. The Fo value of olivine in harzburgite samples chosen by Ling et al. (2014a) [13] was 84, which was five lower than in the samples used in this study. It resulted in a low MgO content, a large error in the composition of the parental magma, and a low degree of partial melting of the mantle. The MgO content of parental magma in Niubiziliang is significantly higher than that of Jinchuan (12.6% [38]) and Xiarihamu (12.48% [24] and 12.74% [39]), and the NiO content in it could be estimated to be 0.053~0.068% by using the discrimination diagram of MgO-NiO (Figure 13b). Meanwhile, the Mg # value of parental magma was calculated to be 71 by the formula of Mg # = 1/(1 + (1 − Fo)/(Kd × Fo)) (Kd = 0.3, Fo = 89), which is close to the Mg # value of primary basaltic magma. Therefore, the parental magma of the Niubiziliang complex is considered to be high-Mg picritic basalt, which also suggests high-degree melting of depleted mantle (ε Hf(t) = 4.2~10.9 [16]).

Olivine
The crystallization temperature of olivine is generally calculated by the geological thermometer formula of T(°C) = 1056.6 + 17.3 × MgO ± 26 [42], where the MgO, here, presents its content in parental magma. Therefore, the crystallization temperature of olivine in Niubiziliang was calculated to be 1289~1369 °C by the geological thermometer formula above. Meanwhile, the CaO content of olivine was low (0~0.35%), indicating a high crystallization pressure [43], that it had crystallized in the magma chamber, and that olivine is frequently surrounded by pyroxene in a round shape to form a poikilitic texture.

Olivine
The crystallization temperature of olivine is generally calculated by the geological thermometer formula of T( • C) = 1056.6 + 17.3 × MgO ± 26 [42], where the MgO, here, presents its content in parental magma. Therefore, the crystallization temperature of olivine in Niubiziliang was calculated to be 1289~1369 • C by the geological thermometer formula above. Meanwhile, the CaO content of olivine was low (0~0.35%), indicating a high crystallization pressure [43], that it had crystallized in the magma chamber, and that olivine is frequently surrounded by pyroxene in a round shape to form a poikilitic texture.

Pyroxene
Since orthopyroxene and clinopyroxene widely existed in the Niubiziliang mafic-ultramafic complex, two-pyroxene thermobarometry proposed by Putirka (2008) [44] was very appropriate for estimating their crystallization temperature and pressure. Firstly, we used the T-independent equation of P(kbar) = −279.8 + 293X opx Al(VI) + 455X opx Na + 299X opx Cr + 519X opx Fm 2 Si 2 O 6 − 563X opx En + 371X opx Di + 372α opx En + 1.19/K f to obtain the crystallization pressure of ortho-and clinopyroxene for the two groups, as follows: 1.76~2.89 kbar (average 2.30 kbar) and 9.38~9.67 kbar. The K f = X opx Ca / (1 − X cpx Ca ) is as in Mercier et al. (1984) [45]; Fm 2 Si 2 O 6 opx = EnFs opx , X opx En = (X opx Fm 2 Si 2 O 6 )(X opx Mg /[X opx Mg , and here, X opx Fe 2+ = X opx Fe − X opx Fe 3+ , where Fe 3+ was calculated as in Papike et al. (1974) [46], i.e., Fe 3+ = Al IV + Na − Al VI -Cr − 2Ti. Secondly, their corresponding crystallization depth was calculated to be 6.65~10.92 km (average 8.7 km) and 35.45~36.55 km using the equation d(km) = 100 × P(kbar)/(ρg), where ρ is 2.7 g/cm 3 for the average basement of Jinshuikou Group gneiss. In addition, the Al IV /(Al IV + Al VI ) value of clinopyroxene was 0.32~0.78, with an average of 0.56, indicating the relatively low crystallization pressure [47]. This further suggested that clinopyroxene likely crystallized after the magma intruded the existing space and its fractional crystallization was not remarkable. Therefore, we suggest that 6.65~10.92 km could be the invaded depth of the Niubiziliang mafic-ultramafic complex, while 34.45~36.55 km likely presents the depth of the magma chamber. Finally, we calculated the crystallization temperatures of ortho-and clinopyrexene to be 1185~1335 • C, with an average of 1253 • C, by using the thermometer of 10 4 /T( [44], which is slightly lower than that of olivine.

Hornblende
The thermometer of T( • C) = −151.487 × Si* + 2041 proposed by Ridolfi et al. (2009) [37] is one of the most widely used and precise thermometers for estimating the crystallization temperature of hornblende. Here, Si* = Si + Al IV /15 − 2Ti − Al/2 + Ti/1.8 + Fe 3+ /9 + Fe 2+ /3.3 + Mg/26 + Ca B /5 + Na B /1.3 − Na A /15 + K A /2.3, where cations in this thermometer are calculated using 23 standard oxygen atoms. The crystallization temperature of hornblende, except for tremolite, in Niubiziliang was calculated to be 1098~1234 • C, with an average of 1159 • C. Secondly, we used the calcareous amphibole formula of P(kbar) = 19.209 × e 1.438Al(Tot) /100 [48], where Al(Tot) = Al IV + Al VI , to estimate the crystallization pressure of hornblende, except for tremolite, resulting in a value of 1.68~2.77 kbar, with an average of 2.15 kbar. Moreover, we also calculated the crystallization depth by applying the physical equation for the gravitational pressure of d(km) = 100 × P(kbar)/(ρg), where ρ is 2.7 g/cm 3 for the average basement of Jinshuikou Group gneiss. Therefore, it could be estimated that the crystallization depth of hornblende in the Niubiziliang complex is about 6.35~10.47 km, with an average of 8.13 km.
According to the crystallization temperature mentioned above, the crystallization sequence of silicate minerals in Niubiziliang is as follows: Olivine→orthopyroxene→clinopyroxene-(plagioclase)→hornblende. This is consistent with the observation results under the microscope, and reflects the mineral crystallization at different temperatures, pressures, and depths. Olivines and some orthopyroxenes first crystallized in the magma chamber at a depth of 35.45~36.55 km, with a pressure and temperature of 9.38~9.67 kbar and 1289~1369 • C, respectively, and they are the optimum mineral for restoring the parental magma. However, most of the silicate minerals crystallized in the existing space at a depth of 8.13~8.70 km, with a pressure and temperature of 2.15~2.30 kbar and 1159~1253 • C, respectively. Finally, the Niubiziliang complex and different kinds of mafic-ultramafic rocks formed through magmatic crystallization and crustal contamination.

Fractional Crystallization
Fractional crystallization is a geological process that occurs extensively in the process of mafic-ultramafic diagenesis and mineralization, and mainly takes place in the early stage of magmatic crystallization. Niubiziliang is a representative mafic-ultramafic complex with multiple invasions by magma in NQOB. It has complete lithofacies, diverse rock types, clear lithofacies boundaries, and a high degree of differentiation. Cumulates such as peridotite, pyroxenite, and anorthosite are all formed by the strong fractional crystallization of primitive basaltic magma in the magma chamber. Meanwhile, some banded gabbros in the No. I ore block (Figure 3a) indicate the local and relatively weak fractional crystallization of magma after invasion. The olivine-encircled texture and interstitial texture were developed in the Niubiziliang complex, which showed that the residual magma between olivine grains formed clinopyroxene to enclose olivine or filled spaces between olivine grains. This is the crystallized result of olivine, which crystallized in the magma chamber at 35.45~36.55 km, and residual magma in the existing space at 8.13~8.70 km.
At the same time, the geochemistry of mafic-ultramafic rocks also shows the significant fractional crystallization in the diagenetic process. The Mg # value of ultramafic rocks in Niubiziliang is 81~84 [16], which is higher than that of primary basaltic magma (68~73 [49,50]), suggesting that fractional crystallization occurred at an early stage. Meanwhile, the NiO content of olivine (0.10~0.23%) is significantly lower than that of olivine crystallized by primitive basaltic magma, which is balanced with pyrolite (≈0.4%), also supporting magma differentiation and leading to a strong deficit of Ni in the parental magma. SiO 2 correlates negatively with MgO (Figure 14a), and Ni, Co, and Cr show positive correlations with MgO (Figure 14d-f), indicative of the crystallization of olivine, and the strong negative correlation between Al 2 O 3 , CaO, and MgO reflects the fractional crystallization of plagioclase or clinopyroxene (Figure 14b,c). If the magma has a significant fractional crystallization of clinopyroxene, in general, the CaO/Al 2 O 3 value and Mg # value show a positive correlation, while the CaO content and CaO/Al 2 O 3 value in the magma also decrease accordingly [51]. However, these characteristics are not clear in the Niubiziliang mafic-ultramafic rocks (Figure 14g,h), indicating that the fractional crystallization of clinopyroxene did not dominate the magma evolution. Despite this, the fractional crystallization of clinopyroxene in gabbro is quite obvious, similar to what occurs for olivine and orthopyroxene in ultramafic rocks during magma evolution (Figure 14i). Simkin et al. (1970) [43] suggested that the CaO content typically increased at a faster cooling rate and low pressures during olivine crystallization. However, the low CaO content of olivine (<0.35%) in Niubiziliang suggests a slow rate and high-pressure (about 9.38~9.67 kbar) conditions [52], which is consistent with what was concluded due to the wide variation of the Fo value and Mg # of olivine. This further suggests that the magma chamber is deep (about 35.45~36.55 km) and sufficiently differentiates, which is favorable for mineralization. However, these characteristics are not clear in the Niubiziliang mafic-ultramafic rocks (Figure 14g,h), indicating that the fractional crystallization of clinopyroxene did not dominate the magma evolution. Despite this, the fractional crystallization of clinopyroxene in gabbro is quite obvious, similar to what occurs for olivine and orthopyroxene in ultramafic rocks during magma evolution (Figure 14i). Simkin et al. (1970) [43] suggested that the CaO content typically increased at a faster cooling rate and low pressures during olivine crystallization. However, the low CaO content of olivine (<0.35%) in Niubiziliang suggests a slow rate and high-pressure (about 9.38~9.67 kbar) conditions [52], which is consistent with what was concluded due to the wide variation of the Fo value and Mg # of olivine. This further suggests that the magma chamber is deep (about 35.45~36.55 km) and sufficiently differentiates, which is favorable for mineralization.

Crustal Assimilation and Contamination
Crustal assimilation and contamination and fractional crystallization typically occur simultaneously during the evolution and emplacement of magma [53], and the crustal assimilation and contamination, which often occur during the formation of a large-scale Ni-Cu deposit, are considered to be some of the key factors in the formation of a magmatic Ni-Cu sulfur deposit (Voisey's Bay [4], Jinchuan [5], Noril'sk [9], and Xiarihamu [54]). First of all, the presence of many xenoliths in the Niubiziliang complex indicates that crustal contamination inevitably occurs during the emplacement of magma. Then, the crustal contamination is not only manifested in the geology, but also usually gives several geochemical fingerprints, such as an increase of SiO 2 , K 2 O, Rb, Ba, Th, Zr, S, 87 Sr/ 86 Sr, 207 Pb/ 204 Pb, δ 18 O, S/Se, and 187 Os/ 188 Os, and a decrease of P 2 O 5 , TiO 2 , Ti/Yb, Ce/Pb, and 143 Nd/ 144 Nd in magma [55]. The sulfur solubility in magma decreases rapidly with an increase of, for example, SiO 2 , CaO, and K 2 O contents, which plays a significant role in the sulfur saturation.
The mantle and crust have quite different values of δ 34 S, δ 18 O, S/Se, 207 Pb/ 204 Pb, 208 Pb/ 204 Pb, 187 Os/ 188 Os, γ Os , ( 87 Sr/ 86 Sr) i , and ε Nd (t), and therefore, they can be effectively used to judge the crustal assimilation and contamination of magma [56,57]. It is generally believed that the metallogenic materials (e.g., S, Pb, Cu, Ni, Co, O, and Se) of magmatic Ni-Cu sulfide deposits are mostly from the mantle, but they change significantly with great variation if contaminated with crustal compositions. Firstly, the S/Se value of the mantle is about 2850~4350, whereas that of the crust is about 3500~10,000, with a varied δ 34 S value of <-40~>30% [58]. Yamamoto (1976) [59] suggested that the increase of the S/Se value was the result of the assimilation of crustal sulfur, while the decrease was likely related to the desulfurization of hydrothermal fluids. The S/Se value of the Niubiziliang mafic-ultramafic rocks is 1864~5890 (Table 5), among which three samples are notably higher than that of the mantle and fall into the S additional area ( Figure 11). Secondly, the δ 34 S value of the mantle is generally within ±2% , while δ 34 S values with a wide range of variation are generally considered to be associated with metamorphism and deposition [21]. The sulfide δ 34 S value in Niubiziliang is -8.4~9.47% , with a significantly wide range of variation (Table 2), which indicates a mixture of sulfur with different origins. It is reasonable to attribute the δ 34 S values of −1.8~1.4% to a mantle origin, but not the δ 34 S values of −8.4~−2.8% and 3.7~9.47% . They are significantly different to that of the upper mantle and are considered to be the result of the sulfur having assimilated from the Paleoproterozoic Jinshuikou Group in NQOB. Thirdly, the sulfides in Niubiziliang have a great variability in terms of the Pb isotopic compositions (  (Figure 10), which also provides strong support for crustal contamination. Last but not least, the whole-rock δ 18 O value of mantle-derived magma is usually modest, at 6.0% , even with the influence of mantle heterogeneity or magmatic crystallization. It is generally believed that there are two mechanisms for increasing the value of δ 18 O in mantle-derived magmatic hydrothermal alteration and crustal contamination. The δ 18 O value of mafic-ultramafic rocks in Niubiziliang is 6.5~8.1% (except for one sample of 5.1% ), with a wide range of variation (Table 4), all of which are higher than that of the mantle. Therefore, it is more likely that the high δ 18 O values of mafic-ultramafic rocks are the result of contamination with the Jinshuikou Group metamorphic rocks, considering that there is no notable hydrothermal alteration in Niubiziliang. In conclusion, they all give strong evidence of crustal assimilation and contamination resulting in greatly elevated S/Se, δ 34 S, 207 Pb/ 204 Pb, 208 Pb/ 204 Pb, and δ 18 O values, with a wide range of variation, and play an important role in the sulfur saturation of parental magma in Niubiziliang. Moreover, hornblendes in Niubiziliang have mantle and crust origins, which could indicate crustal contamination, and the ultrabasic magma and mantle-origin hornblende evolved to the basic magma and crust-origin hornblende due to the addition of crustal felsic components.

Sulfur Saturation Mechanism
Many investigations have indicated that sulfide liquation led to a sharp loss of Ni content in magma under sulfur saturation, and thus the olivine also showed a loss of Ni [3,9]. In other words, the Ni content of olivine was significantly controlled by sulfide liquation, and the Ni deficit in olivine was determined by the volume ratio of olivine to sulfide, in magma. Therefore, a sharp decrease of the Ni content in olivine was considered to be a symbol of deep sulfide liquation, and the greater the Ni deficit, the more complete the sulfide liquation and the more favorable the conditions for mineralization. The Ni content of olivine in Niubiziliang mafic-ultramafic rocks was 0.08~0.18%, which was lower than the normal content of olivine by 0.25%, indicating that, in the magma, sulfide liquation had occurred before olivine crystallization, leading to a loss of Ni content in olivine. Stanley and Karleen (1978) [63] believed that the Ni of olivine crystallized in S unsaturated magma was positively correlated with MgO; on the contrary, it represented that the S of magma had reached saturation. However, the projection points of olivine exhibit a great dispersion and no linear relationship in the Ni-MgO diagram (Figure 15a). At the same time, all of the samples fall in the region of sulfide liquation in the Fo-NiO diagram ( Figure 15b) and both indicate that sulfide liquation occurred in the magma when the olivine crystallized, leading to the NiO deficit in the olivine.
Minerals 2020, 10, x FOR PEER REVIEW 22 of 27 which was lower than the normal content of olivine by 0.25%, indicating that, in the magma, sulfide liquation had occurred before olivine crystallization, leading to a loss of Ni content in olivine. Stanley and Karleen (1978) [63] believed that the Ni of olivine crystallized in S unsaturated magma was positively correlated with MgO; on the contrary, it represented that the S of magma had reached saturation. However, the projection points of olivine exhibit a great dispersion and no linear relationship in the Ni-MgO diagram (Figure 15a). At the same time, all of the samples fall in the region of sulfide liquation in the Fo-NiO diagram (Figure 15b) and both indicate that sulfide liquation occurred in the magma when the olivine crystallized, leading to the NiO deficit in the olivine. The factors affecting the sulfur solubility in magma mainly include the temperature, pressure, oxygen fugacity, sulfur fugacity, and FeO content [2,64], among which the solubility of sulfur in magma is negatively correlated with the pressure, but positively correlated with the temperature [8,65,66]. However, the temperature has little influence on the sulfur solubility compared with the pressure; that is, even if the sulfur in the mantle source is saturated, it remains unsaturated after invading the crust [3]. Therefore, the sulfur in magma tends to be unsaturated, with a decrease of pressure during the rising process of magma, and sulfur in magma has difficulty reaching saturation and forming an Ni-Cu deposit if there is no interference from external factors. A key factor in the formation of magmatic sulfide deposits is the saturation of sulfur in magma, causing immiscibility between the silicate magma and sulfides, and enrichment in chalcophile elements in the sulfide melt [3]. Sulfur saturation in magma can generally be achieved by any combination of the following processes: (a) Rapid cooling, such as the sulfides at the bottom of the Bushveld complex [67]; (b) magma mixing, such as the Merensky Reef in Bushveld [68]; (c) fractional crystallization, such as the sulfides in chromitites in Bushveld, and the Hongqiling, Huangshan, and Kalatongke deposits [6]; (d) the addition of external sulfur, such as the Voisey's Bay, Duluth, Tianyu, and Baishiquan deposits [7]; and (e) the addition of crustal compositions, such as the Noril'sk, Huangshannan, and Kalatongke deposits [9]. It should be pointed out that the addition of external sulfur and crustal compositions are attributed to the crustal assimilation and contamination, which play a key role in Ni-Cu sulfide deposit formation and often occur during the formation of large-super large Ni-Cu deposits (Voisey's Bay [4], Jinchuan [5], Noril'sk [9], and Xiarihamu [54]). There could be one or several sulfur saturation mechanisms for the Niubiziliang Ni-Cu deposit. First of all, the low CaO content of olivine (<0.01%) and sufficient differentiation of lithofacies suggest that the magma is crystallizing at a slow rate and higher pressures. In other words, the rate of magma temperature decline is low, which suggests that the rapid cooling of the magma has a limited role in S saturation. Secondly, there is no magmatic inclusion of other constituents found in the Niubiziliang mafic-ultramafic complex, which implies a low probability of magma mixing for S saturation. Thirdly, as previously mentioned, the complete lithofacies of complex and discrimination diagrams suggest the intense fractional crystallization of magma. They lead to the crystallization of large amounts of mafites (e.g., ilmenite, magnetite, olivine, and pyroxene), which alters the magma composition, in particular, rapidly bringing down the FeO content, resulting in a decrease in the sulfur solubility. The factors affecting the sulfur solubility in magma mainly include the temperature, pressure, oxygen fugacity, sulfur fugacity, and FeO content [2,64], among which the solubility of sulfur in magma is negatively correlated with the pressure, but positively correlated with the temperature [8,65,66]. However, the temperature has little influence on the sulfur solubility compared with the pressure; that is, even if the sulfur in the mantle source is saturated, it remains unsaturated after invading the crust [3]. Therefore, the sulfur in magma tends to be unsaturated, with a decrease of pressure during the rising process of magma, and sulfur in magma has difficulty reaching saturation and forming an Ni-Cu deposit if there is no interference from external factors. A key factor in the formation of magmatic sulfide deposits is the saturation of sulfur in magma, causing immiscibility between the silicate magma and sulfides, and enrichment in chalcophile elements in the sulfide melt [3]. Sulfur saturation in magma can generally be achieved by any combination of the following processes: (a) Rapid cooling, such as the sulfides at the bottom of the Bushveld complex [67]; (b) magma mixing, such as the Merensky Reef in Bushveld [68]; (c) fractional crystallization, such as the sulfides in chromitites in Bushveld, and the Hongqiling, Huangshan, and Kalatongke deposits [6]; (d) the addition of external sulfur, such as the Voisey's Bay, Duluth, Tianyu, and Baishiquan deposits [7]; and (e) the addition of crustal compositions, such as the Noril'sk, Huangshannan, and Kalatongke deposits [9]. It should be pointed out that the addition of external sulfur and crustal compositions are attributed to the crustal assimilation and contamination, which play a key role in Ni-Cu sulfide deposit formation and often occur during the formation of large-super large Ni-Cu deposits (Voisey's Bay [4], Jinchuan [5], Noril'sk [9], and Xiarihamu [54]).
There could be one or several sulfur saturation mechanisms for the Niubiziliang Ni-Cu deposit. First of all, the low CaO content of olivine (<0.01%) and sufficient differentiation of lithofacies suggest that the magma is crystallizing at a slow rate and higher pressures. In other words, the rate of magma temperature decline is low, which suggests that the rapid cooling of the magma has a limited role in S saturation. Secondly, there is no magmatic inclusion of other constituents found in the Niubiziliang mafic-ultramafic complex, which implies a low probability of magma mixing for S saturation. Thirdly, as previously mentioned, the complete lithofacies of complex and discrimination diagrams suggest the intense fractional crystallization of magma. They lead to the crystallization of large amounts of mafites (e.g., ilmenite, magnetite, olivine, and pyroxene), which alters the magma composition, in particular, rapidly bringing down the FeO content, resulting in a decrease in the sulfur solubility. Therefore, fractional crystallization is considered to be one of the factors in S saturation in the Niubiziliang Ni-Cu deposit. Last but not least, the crustal assimilation and contamination in Niubiziliang determine the S saturation in at least three ways, as follows: (a) Absorbing the sulfur in Jinshuikou Group and increasing the S concentration in magma supported by the S/Se and δ 34 S values; (b) altering the composition of magma and decreasing the solubility of sulfur in magma, such as the increase of SiO 2 , K 2 O, Al 2 O 3 , and CaO content; and (c) enhancing the f O2 of magma, resulting in a decrease of the sulfur solubility supported by the δ 18 O values. Therefore, we conclude that the fractional crystallization of magma and crustal assimilation and contamination are the main mechanisms of sulfur saturation in the parental magma of the Niubiziliang Ni-Cu deposit. However, the effect of crustal assimilation and contamination on sulfur saturation is much greater than that of fractional crystallization.
As the primary magma evolved, the fractional crystallization of olivine and pyroxene, together with crustal assimilation and contamination, caused sulfur saturation and strong sulfide immiscibility. Therefore, the parental magma differentiated into a series of magma with various sulfide melting in the magma chamber at a depth of 35.45~36.55 km, including barren, ore-bearing, ore-rich, and ore magma. Then, these magmas were emplaced in multiple stages under the geo-background of continuous extension during the middle-late Devonian in NQOB, and the Niubiziliang mafic-ultramafic complex with Ni-Cu mineralization finally formed in the existing space at a depth of 8.13~8.70 km.

Conclusions
(1) The parental magma of the Niubiziliang complex is considered to be high-Mg picritic basalt, with MgO and NiO contents of 14.95~16.58% and 0.053~0.068%, respectively, indicating a high degree of partial melting of the depleted mantle.
(2) The Niubiziliang complex formed at a moderate depth (8.13~8.70 km) with the temperature and pressure of 1159~1253 • C and 2.15~2.30 kbar, respectively, and the crystallization sequence of the main silicate minerals is olivine→pyroxene-(plagioclase)→hornblende, which is consistent with the observation made under the microscope. However, the olivines and some pyroxenes likely crystallized in the magma chamber at a depth of 35.45~36.55 km, with a high temperature (1289~1369 • C) and pressure (9.38~9.67 kbar).
(3) The primary magma of Niubiziliang underwent significant fractional crystallization and crustal assimilation and contamination, which resulted in sulfur saturation and sulfide immiscibility in the magma. However, the contribution of crustal assimilation and contamination to sulfur saturation is much greater than that of fractional crystallization.