# The Fe-Zn Isotopic Characteristics and Fractionation Models: Implications for the Genesis of the Zhaxikang Sb-Pb-Zn-Ag Deposit in Southern Tibet

Authors: Da Wang; Youye Zheng; Ryan Mathur; Song Wu
Journal: Geofluids
 (2018)

---

## Abstract

The genesis of the Zhaxikang Sb-Pb-Zn-Ag deposit remains controversial. Three different geological environments have been proposed to model mineralization: a hot spring, a magmatic-hydrothermal fluid, and a sedimentary exhalative (SEDEX) overprinted by a hot spring. Here, we present the electron probe microanalysis (EPMA) and Fe-Zn isotopic data (microsampled) of four samples from the first pulse of mineralization that show annular textures to constrain ore genesis. The Zn/Cd ratios from the EPMA data of sphalerite range from 296 to 399 and overlap the range of exhalative systems. Theδ
56Fe values of Mn-Fe carbonate and δ
66Zn values of sphalerite gradually decrease from early to late stages in three samples. A combination of the EPMA and isotopic data shows the Fe-Zn contents also have different correlations with δ
66Zn values in sphalerite from these samples. Rayleigh distillation models this isotope and concentration data with the cause of fractionation related to vapour-liquid partitioning and mineral precipitation. In order to verify this Rayleigh distillation model, we combine our Fe-Zn isotopic data with those from previous studies to establish 12 Fe-Zn isotopic fractionation models. These fractionation models indicate the δ
56Fei and δ
66Zni values (initial Fe-Zn isotopic compositions) of the ore-forming system are in the range of -
0.5
‰​
​
~−1
‰ and -
0
.
28
‰
​
~
0
‰, respectively. To conclude, the EPMA data, Fe-Zn isotopic characteristics, and fractionation models support the SEDEX model for the first pulse of mineralization.


---

## Body

## 1. Introduction

To date, the Zhaxikang Sb-Pb-Zn-Ag deposit is the only super large deposit that has been identified within the North Himalayan Polymetallic Metallogenic Belt (NHMB). Although basic research including the geology, petrography, geochronology, and geochemistry studies has been conducted (e.g., [1, 2]), the genesis of this deposit is still debated due to the complicated mineralogy and the presence of multiple stages of mineralization. The main viewpoints involve a hot spring [3, 4], two magmatic-hydrothermal fluids [5, 6], and a SEDEX overprinted by hot spring [7] genetic models. However, most of these genetic models are based on the S, C, O, H, and Si isotopic evidence, which cannot absolutely trace the metal source.The traditional light stable C, H, O, S, and N isotopes have been widely used to constrain fluid evolution and metal sources in ore deposit studies (e.g., [8, 9]). However, the evidences for metal source from these elements are always indirect and putative as they are not the metallogenic elements themselves [10]. For instance, these elements usually have different characteristics with the changing of tectonic settings, and sometimes they may even have different sources from the metallogenic elements [11]. However, the nontraditional transition metal stable isotopes (e.g., Fe, Zn, Cu, Cd, Mg, Cr, Sn, and Mo) are more precise tracers for the metal sources and ore-forming processes in metallogenic systems [10, 12, 13]. The development of Multicollector-Inductively Coupled Plasma Mass Spectrometer (MC-ICP-MS) technology has greatly improved the precision of isotopic analyses [14, 15], which results in the wide application of the nontraditional transition metal stable isotopes in economic geology studies (e.g., [16, 17]).The Fe-Zn isotopes are two of the most representative isotopes applied in ore deposit studies. For example, Mason et al. [18] and Wilkinson et al. [19] both identified that the δ
66Zn values of minerals precipitating from the same hydrothermal fluid become heavier over time by studying the Zn isotopic fractionation of the Alexandrinka volcanic hosted massive sulfide (VHMS) type deposit in Russia and Midlands Irish-type deposit in Ireland, respectively. The gradual increasing δ
66Zn values both from early to late stages and from south to north within the Red Dog ore district in Alaska record the temporal-spatial evolution of the ore-forming fluid and constrain the SEDEX genesis with a single Zn source [20]. In addition, Fe isotopic studies of skarn Cu−Au±Fe deposits in South China excluded the possibility that Fe originated from sedimentary strata [21, 22]. These Fe isotopes matched the igneous source rocks and mineralization, and the δ
56Fe values of sulfides gradually increase both from early to late stages and away from the ore-related igneous rocks. Wang et al. [21, 22] also revealed that the Fe isotopes fractionate during fluid exsolution and that the ore-forming fluid is enriched in light isotopes relative to ore-related igneous rocks. To the contrary, Wawryk and Foden [23] investigated the Fe-isotope fractionation in the Renison Sn-W deposit in Australia and discovered that Fe isotopic compositions of pyrite (0.61
‰
​
~
1.14
‰), chalcopyrite (0.18
‰
​
~
0.71
‰), and magnetite (0.50
‰
​
~
0.70
‰) are isotopically heavier than Renison granite (0.27
‰
​
~
0.45
‰) and thus hypothesized that a magmatic-hydrothermal fluid exsolved from an isotopically heavy reduced magma could deposit isotopically heavy ore minerals whereas oxidized magmas crystallise magmatic magnetite could result in an isotopically lighter melt and fluid. These studies demonstrate the potential of Fe-Zn isotopes to trace the metal source and provide insights into ore-forming evolution.With regard to the Zhaxikang deposit, Duan et al. [5] have investigated the Zn isotope of sphalerite, galena, Fe-Mn carbonates, and igneous rocks and speculated that the δ
66Zn values of the hydrothermal fluid are 0.39 ± 0.10‰. This value is consistent with those of basement rocks (average value of 0.36‰  ± 0.03‰) and Fe-Mn carbonates (average value of 0.27‰  ± 0.15‰), which is identified as the evidence for the magmatic origin. Meanwhile, the contribution of regional sedimentary rocks is conjectured by the Zn-Pb-S isotopes: (
1
) the Zn isotopic variation range of sulfides (-
0.25
‰
​
~
0.03
‰) is larger than basement rocks (0.05
‰
​
~
0.21
‰); (
2
) the radiogenic Pb isotopic compositions of sulfides (e.g., 206Pb/204Pb = 18.727~19.896) is higher than regional igneous rocks (206Pb/204Pb = 18.4~19.2); (
3
) the δ
34S values of sulfides (6
‰
​
~
12
‰) are lighter than regional sedimentary wall rocks (10.94
‰
​
~
11.49
‰) but higher than mantle value (0 ± 2‰). However, the Zn isotopic fractionation during the fluid exsolution and leaching process [27, 28] has been ignored in Duan et al. [5]. In addition, Wang et al. [26] also studied the Fe-Zn isotopes of the pyrite, sphalerite, and Mn-Fe carbonate in Zhaxikang deposit, which successfully constrained the two pulses of mineralization by the temporally increasing δ
56Fe and decreasing δ
66Zn values recorded in the deposit that coincided with an increase in alteration. The Fe-Zn isotopic research also demonstrated the magmatic-hydrothermal fluid origin of the second pulse of mineralization by the heavier δ
56Fe values of stage 3 pyrite and excluded the possibility that slate is the metal source by the similar δ
66Zn values of slate and sphalerite. Nevertheless, the attempt to trace the metal source of the first pulse of mineralization failed in Wang et al. [26]. In this study, we present the Fe-Zn isotopic values and variations within four annular polished section samples from Zhaxikang deposit to provide more credible evidence for the primary and earlier stages of ore genesis.


## 2. Geological Setting

### 2.1. Regional Geology

The NHMB is in the eastern section of the North Himalayan Tectonic Belt (NH), in the Himalayan terrane. From north to south, the Himalayan terrane is divided into four tectonic belts: the North Himalayan Tethys Sedimentary Fold Belt (the North Himalayan Tectonic Belt; NH), the High Himalayan Crystalline Rock Belt (HH), the Low Himalayan Fold Belt (LH), and the Sub-Himalayan Tectonic Belt (SH; Figure1(a)) [29–31]. These belts are separated by four nearly EW-trending faults, including the South Tibet Detachment System (STDS), the Main Central Thrust (MCT), the Main Boundary Thrust (MBT), and the Main Frontal Thrust (MFT; Figure 1(a)) [24, 32]. The NH, composed of a set of Palaeozoic marine sedimentary sequences that formed in a passive continental margin environment within northern India [33], is located between the Indus-Yarlung Zangbo Suture Zone (IYZS) and the HH (Figure 1(a)).Figure 1
(a) Tectonic framework of the Himalayan Terrane (modified from Yin [24]). (b) Regional geological map of the North Himalayan Polymetallic Metallogenic Belt (modified from Zheng et al. [25]). BNSZ: the Bangong-Nujiang Suture Zone; IYZS: the Indus-Yarlung Zangbo Suture Zone; STDS: the South Tibet Detachment System; MCT: the Main Central Thrust Fault; MBT: the Main Boundary Thrust Fault; MFT: the Main Frontal Thrust; NH: the North Himalayan Tethys Sedimentary Fold Belt; HH: the High Himalayan Crystalline Rock Belt; LH: the Low Himalayan Fold Belt; SH: the Sub-Himalayan Tectonic Belt.The sedimentary sequence in the NH records Late Precambrian to Devonian prerift, Carboniferous to Early Jurassic syn-rift and Middle Jurassic to Cretaceous passive continental margin sediments (Figure1(b)) [66–68]. These sediments crop out in an EW and NWW trend and are predominated by the Precambrian Laguigangri Group and a series of Upper Triassic, Jurassic, Lower Cretaceous, and Quaternary sediments. The Laguigangri Group crops out in the core of the Yelaxiangbo dome and is composed of schist, gneiss, and migmatite units (Figure 1(b)). A set of Late Triassic-Early Cretaceous flysch formations deposit in neritic-bathyal environments and crop out across the NH. This set of formations are dominated by turbidite deposits and host the majority of the Au-Sb-Pb-Zn-Ag deposits in the NH. The lithology of these Late Triassic-Early Cretaceous flysch formations is weak-metamorphic slate that is intercalated with metamorphosed fine-grained sandstone, argillaceous limestone, micrite, and siliceous rock that is intercalated with volcanic rocks [7]. Some quaternary sediments also occur in the central and southern area of the NH, which are composed of gravel, sand gravel, sandy loam, clay, and ice boulder.The EW-trending and NS-trending faults, both with the several episodes of motion, are present in the NH. The EW-trending faults, controlling the distribution of intermediate-acid magmatic rocks and ore deposits in the NH [7], are older and cover a larger area than the NS-trending faults. The representative EW-trending faults include the Lazi-Qiongduojiang, Rongbu-Gudui, and Luozha faults as well as the STDS and numerous metamorphic core complexes (Figure 1(b)). A series of rifts that formed from 25 Ma to present are associated with these faults [29, 69–71], mainly including the Sangri-Cuona, Yadong-Gulu, Shenzha-Xietongmen, and Dangreyongcuo-Gucuo rift zones from east to west [72]. In addition, the NS-trending faults that are considered as the result of east-west extension of the Qinghai-Tibet Plateau [73, 74] also formed during this period especially from 18 to 4 Ma [75, 76]. These NS-trending faults are also the important ore-controlling structures in the NH [1].Magmatism in the NH primarily includes the Mesozoic and Cenozoic magmatism. The Mesozoic magmatism generated multiple suites of mafic-intermediate igneous rocks between the Late Triassic and the Early Cretaceous, including basaltic volcanic interlayers, dyke swarms, and subvolcanic dykes. According to the previous geochronological data, the SHRIMP U-Pb ages of the basic dyke swarms from different area in the NH are 134.9 ± 1.8 Ma, 135.5 ± 2.1 Ma [77], and 138.0 ± 3.5 Ma [78], respectively. The SHRIMP U-Pb age of the gabbro is 155.8 Ma [79]. Tong et al. [78], Pan et al. [80], and Zhong et al. [81] regarded these basic dyke swarms as the result of late-stage massive expansion of Neo-Tethys Ocean under the structural environment of the Himalaya passive continental margin intensive stretching and breaking-off, lithosphere extension-thinning, and asthenosphere upwelling. On the contrary, Zhu et al. [82] and Qiu et al. [83] suggested that these basic dyke swarms are the result of interaction between mantle plume and lithospheric mantle material and form in the continental-rift environment. The Cenozoic magmatism is characterized by the formation of monzogranite, leucogranite, diorite, porphyritic diorite, and aplite units [84, 85]. These Cenozoic intermediate-acidic intrusive masses are distributed along the EW-trending faults and in the core of Ranba, Kangma, and Yelaxiangbo dome in the form of batholith, laccolith, and dykes (Figure 1(b)). This phenomenon is considered to be the result of crustal thickening [86] related to the collision of the India Plate and the Eurasia Plate during the postcollision stage (25 to 0 Ma) [87, 88].The NHMB contains many Sb, Au, Sb-Au, Pb-Zn, and Sb-Pb-Zn-Ag deposits, and the Zhaxikang Sb-Pb-Zn-Ag, the Mazhala Au-Sb, the Chalapu Au, the Bangbu Au, the Shalagang Sb, the Cheqiongzhuobu Sb deposits are representative (Figure1(b)) [7, 25]. The geneses and metallogenic age of these deposits are controversial due to the limited research, the genetic models mainly include the SEDEX overprinted by hot spring, carlin and carlin-like, hot spring, subvolcanic magmatic-hydrothermal fluid, and orogenic types [25].


### 2.2. Ore Deposit Geology

The Zhaxikang Sb-Pb-Zn-Ag deposit is located ~48 km west from Longzi County Town within the southeastern Yangzuoyong-Nariyong composite syncline in the NH (Figure1(b)). This deposit has a reserve of 1.268 Mt Pb-Zn with a 3.66% average Zn grade and a 2.45% average Pb grade, 0.138 Mt Sb with an average grade of 1.08%, 1800 t Ag with an average of 99.55 g/t, 3.9 t associated Au, 361 t associated Ga, and 20 Mt Mn-Fe carbonate ores with an average grade of 42% for Fe + Mn [89], which makes it the largest deposit within the NHMB. The majority of mineralization in the orefield is hosted by the Lower Jurassic Ridang formation that consists of epimetamorphic marine clastic rocks. This formation, dipping shallowly to the north and striking eastwest, is divided into five lithologic sections (Figure 2(a)). A few Upper Jurassic Weimei formations composed of fine-grained metamorphic quartzose sandstone, silty slate, and calcarenite as well as Quaternary sediments distributed along valleys also crop out in the orefield (Figure 2(a)) [7].Figure 2
(a) Geological map of the Zhaxikang Sb-Pb-Zn-Ag polymetallic deposit (modified from Zheng et al. [7]). (b) Cross-section along Exploration Line 7. J
1
r
1 to J
1
r
5: the first to fifth lithologic section of the Lower Jurassic Ridang formation.
(a)
(b)The Zhaxikang deposit developed extensive geological structures. A near northsouth striking fault system is prevalent in the orefield, which coexists with a group of northeast-striking faults and some folds. Engineering and geological mapping projects have identified 16 faults, the majority of which are steeply dipping normal faults associated with tensional stress and wrench faults associated with torsional stress. Faults F2, F4, F5, F6, F7, F13, F14, and F16 are the main ore-bearing faults, faults F1 and F10 are partly mineralized, fault F3 was associated with late-stage mineralization, faults F8 and F9 are wrench faults without any mineralization, and fault F12 is a nonmineralized regional fault (Figure2(a)). The orebodies I–VI are hosted by nearly NS-striking faults and orebodies VII–IX are hosted by nearly NE-striking faults (Figure 2). Our samples in this study are all from the orebody V, which is the largest and richest one among these orebodies within the orefield and hosts more than 80% of the reserves. This orebody is >1400 m long, 1 to 30 m wide, and controlled by fault F7 (Figure 2).The magmatism in the orefield is associated with diabase, porphyritic rhyolite, basalt, and leucogranite units as well as some granite porphyry dykes that intruded into the porphyritic rhyolite (Figure2(a)). The diabase is identified by drillholes and footrill in the central part of the orefield as dykes that emplaced into the Jurassic Ridang Formation and has been dated by zircon U-Pb methods to ~133 Ma [7]. The rhyolite porphyry with the zircon SHRIMP U-Pb age of ~135 Ma crops out in the western part of the orefield [90] and the leucogranites crop out in the southern part of the orefield over an area of <1 km2. Additionally, the basalt usually occurs near the orebody in the form of consequent layer or shear layer distributed in slate and the contact region of slate and diabase.Various types of alteration associated with mineralization have occurred in the orefield, including (1) the silicification that is associated with Sb mineralization and generally located in fault zones in the form of quartz veins, radiating quartz, and quartz druse; (2) the carbonatization that is associated with Pb-Zn mineralization in the form of Mn-Fe carbonate veins and also formed the postmineralization calcite; (3) the chlorite alteration that is generally confined to massive and stellated aggregates of chlorite within diabase; (4) the weak sericite alteration that is associated with chlorite alteration and barren quartz; and (5) the clay alteration that developed along the edges of fracture-related crushed zones. Furthermore, the ore-forming elements display a vertical sequence that is zoned from a lowermost Zn (Pb + Ag) zone through a central Zn + Pb + Ag-(Sb) zone to an uppermost Pb + Zn + Sb + Ag zone, although no horizontal zoning is present [89].


### 2.3. Ore Paragenetic Sequence

The paragenetic sequence in the Zhaxikang deposit is divided into six stages of ore formation based on the detailed hand specimen and microscopic observations. These six stages are assigned to two clear pulses: the first pulse consists of stages 1 to 2 and is characterized by the assemblages of Mn-Fe carbonates and sulfides, and the second pulse includes stages 3 to 6 and is primarily dominated by quartz, calcite, sulfosalt minerals, and sulfides (Figure3).Figure 3
Mineral paragenesis within the Zhaxikang Sb-Pb-Zn-Ag deposit (modified from Wang et al. [26]).Stage 1, dominated by a Mn-Fe carbonate + sphalerite + pyrite + arsenopyrite assemblage, is the initial stage of ore formation in the Zhaxikang deposit. Majority of the fine-grained sphalerite, pyrite, and arsenopyrite are hosted by fine-grained Mn-Fe carbonate in the form of laminae (Figures4(a)–4(c)), and a few sulfides occur within the Mn-Fe carbonates as stellated aggregates (Figure 4(d)). The fine-grained layered and colloform with synsedimentary features (Figure 4(b)). The laminae and Mn-Fe carbonates in some samples have been cut by later stage 4 quartz-boulangerite veins (Figure 4(c)) or have been affected by the stage 2 coarse-grained sphalerite (Figure 4(a)).Figure 4
Hand specimen photographs and photomicrographs of representative samples from the Zhaxikang deposit. (a) Stage 1 lamellar sphalerite-pyrite-arsenopyrite and stage 2 massive sphalerite-pyrite hosted within fine-grained Mn-Fe carbonate. (b) Stage 1 lamellar and stage 2 banded Mn-Fe carbonate-sphalerite-galena ore with visible synsedimentary features including rhythmic sedimentation in the upper part and angular folding in the lower part of the sample. (c) Stage 1 lamellar sphalerite-pyrite-arsenopyrite and stage 2 massive and banded sphalerite-pyrite hosted by fine-grained Mn-Fe carbonate. The mineral assemblage is in turn cross-cut by stage 4 quartz-boulangerite veins. (d) Coarse-grained stage 2 Mn-Fe carbonate-sphalerite formed by the recrystallization of fine-grained stage 1 Mn-Fe carbonate-sphalerite. (e) Stage 2 massive coarse-grained pyrite hosted by slate. (f) Massive and brecciated stage 2 sphalerite hosted in stage 2 Mn-Fe carbonate. (g) Stage 2 massive, globular, and concentric annular sphalerite-pyrite hosted by coarse-grained Mn-Fe carbonate. (h) Stage 2 massive galena and pyrite. (i) Stage 2 massive and veined sphalerite-pyrite hosted by coarse-grained Mn-Fe carbonate. (j) Stage 1 lamellar sphalerite-pyrite-arsenopyrite and stage 2 massive sphalerite hosted by fine-grained Mn-Fe carbonate. The sample also contains a Mn-Fe carbonate druse dominated by idiomorphic columnar quartz and valentinite. (k) Stage 2 coarse-grained sphalerite-pyrite hosted by stage 2 Mn-Fe carbonate with banded textures. (l) Stage 2 sphalerite and Mn-Fe carbonate ore with typical Dal Matianite texture. (m) Disseminated stage 3 sphalerite and pyrite hosted in stage 3 quartz and calcite, and stage 3 quartz and calcite cut the slate. (n) Stage 3 brecciated sphalerite within stage 3 quartz-calcite. (o) Stage 3 sphalerite-galena veins cross-cut by stage 6 quartz-calcite veins. (p) Stage 3 sphalerite cross-cut by stage 4 quartz-boulangerite veins. (q) Stage 4 massive and needle-like boulangerite hosted by stage 4 quartz. (r) Stage 4 boulangerite-quartz. (s) Stage 5 elongate stibnite hosted by stage 5 quartz. (t) Stage 5 stibnite-cinnabar hosted by stage 5 quartz. (u) Siliceous sinter formed during the Supergene stage. (v) Ferrihydrite and sardinianite formed during the Supergene stage. (w) Stage 2 pyrite containing automorphic stage 1 arsenopyrite is replaced by later stage 2 sphalerite to form a skeletal texture. (x) Stage 3 sphalerite occurs in stage 3 quartz. (y) The emulsion-like and disseminated stage 3 chalcopyrite grains are dotted in the stage 3 sphalerite. (z) The stage 3 chalcopyrite grains are dotted among the grains of stage 3 pyrite. (aa) Stage 3 sphalerite replaced by stage 4 boulangerite that is in turn cross-cut by stage 5 stibnite. (ab) Stage 5 stibnite cross-cut by stage 6 quartz. Mcar1 = stage 1 fine-grained Mn-Fe carbonate; Apy1 = stage 1 lamellar arsenopyrite; Py1 = stage 1 lamellar pyrite; Sp1 = stage 1 lamellar sphalerite; Mcar2 = stage 2 coarse-grained Mn-Fe carbonate; Apy2 = stage 2 arsenopyrite; Py2 = stage 2 pyrite; Sp2 = stage 2 sphalerite; Gn2 = stage 2 coarse-grained galena; Py3 = stage 3 pyrite; Sp3 = stage 3 sphalerite; Gn3 = stage 3 galena; Ccp3 = stage 3 chalcopyrite; Qtz3 = stage 3 quartz; Cal3 = stage 3 calcite; Blr4 = stage 4 boulangerite; Qtz4 = stage 4 quartz; Stb5 = stage 5 stibnite; Ci5 = stage 5 cinnabar; Qtz5 = stage 5 quartz; Cal6 = stage 6 calcite; Qtz6 = stage 6 quartz; Lm = Supergene stage ferrihydrite.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)
(q)
(r)
(s)
(t)
(u)
(v)
(w)
(x)
(y)
(z)
(aa)
(ab)Stage 2 hosts majority of the Pb-Zn mineralization in the Zhaxikang deposit and comprises an assemblage of Mn-Fe carbonate + galena + sphalerite + pyrite ± arsenopyrite. The more abundant and coarser-grained sulfides are hosted by coarse-grained Mn-Fe carbonate and slate to form the banded (Figure4(k)), net-veined (Figure 4(i)), massive (Figure 4(h)), concentric annular (Figures 4(g) and 5), and Dal Matianite (Figure 4(l)) ores. The Mn-Fe carbonate during this stage recrystallized to different degree (Figures 4(f), 4(g) and 4(i)–4(k)), some even formed the druse containing the idiomorphic columnar quartz, needle-like boulangerite, or valentinite (Figure 4(j)). We can also observe that the later sphalerite replaces the earlier pyrite containing automorphic stage 1 arsenopyrite to form a skeletal texture during this stage (Figure 4(w)). Zheng et al. [7] considered that the ore textures in stages 1 and 2 are similar to those of the Red Dog SEDEX-type ore district in Alaska.Figure 5
The sampling points and Fe-Zn isotopic-elemental variations in the annular polished section samples: (a) 9-8; (b) 9-3; (c) ZXK-1; (d) ZXK-2. Mcar = Mn-Fe carbonate; Apy = arsenopyrite; Py = pyrite; Sp = sphalerite; Qtz = quartz.
(a)
(b)
(c)
(d)Stage 3, characterized by the formation of a quartz ± calcite + pyrite + sphalerite + galena ± chalcopyrite ± arsenopyrite assemblage without Mn-Fe carbonate, is the earliest stage of the second pulse of mineralization. The massive, veined, net-veined, and brecciated sphalerite, galena, and pyrite occur in the quartz and calcite (Figures4(l)–4(p), 4(x)), and most of the sulfides form by the modification of sulfides from earlier stages. Some chalcopyrite grains are distributed in sphalerite, galena, and pyrite as stellated aggregates (Figures 4(y) and 4(z)). Some of these sulfides have been cross-cut by the later quartz-boulangerite or quartz-calcite veins (Figure 4(o)).Stage 4 is marked by a mineral assemblage composed of quartz + antimony-lead-silver sulfosalt minerals (Figure4(r)) that prevailingly include boulangerite and jamesonite, as well as minor bournonite, tetrahedrite, and andorite. This stage hosts the majority of the Sb and Ag mineralization and yields the ores with relatively high average Ag grades. The minerals formed in earlier stages are replaced and cross-cut by the quartz-boulangerite veins and boulangerite of this stage (Figures 4(c) and 4(aa)). Some samples also contain quartz druse filled with needle-like boulangerite (Figure 4(q)).Stage 5 is distinguished by the formation of a quartz + stibnite + cinnabar assemblage and hosts part of the Sb mineralization within the deposit. Elongate-radial stibnite and massive stibnite-cinnabar are hosted by quartz (Figures4(s) and 4(t)). Some stibnite cross-cut the stage 4 boulangerite (Figure 4(aa)).Stage 6, representing the youngest stage of mineralization in the Zhaxikang deposit, is identified by a quartz ± calcite assemblage without sulfides. The quartz-calcite veins of this stage cross-cut the earlier formed minerals (Figures4(o) and 4(ab)). Zheng et al. [7] regarded the ore textures in the second pulse of mineralization as typical hot spring type metallogenic features.Supergene stage in the Zhaxikang deposit consists of ferrihydrite, smithsonite, sardinianite, valentinite, travertine, malachite, and siliceous sinter (Figures4(u) and 4(v)).


## 2.1. Regional Geology

The NHMB is in the eastern section of the North Himalayan Tectonic Belt (NH), in the Himalayan terrane. From north to south, the Himalayan terrane is divided into four tectonic belts: the North Himalayan Tethys Sedimentary Fold Belt (the North Himalayan Tectonic Belt; NH), the High Himalayan Crystalline Rock Belt (HH), the Low Himalayan Fold Belt (LH), and the Sub-Himalayan Tectonic Belt (SH; Figure1(a)) [29–31]. These belts are separated by four nearly EW-trending faults, including the South Tibet Detachment System (STDS), the Main Central Thrust (MCT), the Main Boundary Thrust (MBT), and the Main Frontal Thrust (MFT; Figure 1(a)) [24, 32]. The NH, composed of a set of Palaeozoic marine sedimentary sequences that formed in a passive continental margin environment within northern India [33], is located between the Indus-Yarlung Zangbo Suture Zone (IYZS) and the HH (Figure 1(a)).Figure 1
(a) Tectonic framework of the Himalayan Terrane (modified from Yin [24]). (b) Regional geological map of the North Himalayan Polymetallic Metallogenic Belt (modified from Zheng et al. [25]). BNSZ: the Bangong-Nujiang Suture Zone; IYZS: the Indus-Yarlung Zangbo Suture Zone; STDS: the South Tibet Detachment System; MCT: the Main Central Thrust Fault; MBT: the Main Boundary Thrust Fault; MFT: the Main Frontal Thrust; NH: the North Himalayan Tethys Sedimentary Fold Belt; HH: the High Himalayan Crystalline Rock Belt; LH: the Low Himalayan Fold Belt; SH: the Sub-Himalayan Tectonic Belt.The sedimentary sequence in the NH records Late Precambrian to Devonian prerift, Carboniferous to Early Jurassic syn-rift and Middle Jurassic to Cretaceous passive continental margin sediments (Figure1(b)) [66–68]. These sediments crop out in an EW and NWW trend and are predominated by the Precambrian Laguigangri Group and a series of Upper Triassic, Jurassic, Lower Cretaceous, and Quaternary sediments. The Laguigangri Group crops out in the core of the Yelaxiangbo dome and is composed of schist, gneiss, and migmatite units (Figure 1(b)). A set of Late Triassic-Early Cretaceous flysch formations deposit in neritic-bathyal environments and crop out across the NH. This set of formations are dominated by turbidite deposits and host the majority of the Au-Sb-Pb-Zn-Ag deposits in the NH. The lithology of these Late Triassic-Early Cretaceous flysch formations is weak-metamorphic slate that is intercalated with metamorphosed fine-grained sandstone, argillaceous limestone, micrite, and siliceous rock that is intercalated with volcanic rocks [7]. Some quaternary sediments also occur in the central and southern area of the NH, which are composed of gravel, sand gravel, sandy loam, clay, and ice boulder.The EW-trending and NS-trending faults, both with the several episodes of motion, are present in the NH. The EW-trending faults, controlling the distribution of intermediate-acid magmatic rocks and ore deposits in the NH [7], are older and cover a larger area than the NS-trending faults. The representative EW-trending faults include the Lazi-Qiongduojiang, Rongbu-Gudui, and Luozha faults as well as the STDS and numerous metamorphic core complexes (Figure 1(b)). A series of rifts that formed from 25 Ma to present are associated with these faults [29, 69–71], mainly including the Sangri-Cuona, Yadong-Gulu, Shenzha-Xietongmen, and Dangreyongcuo-Gucuo rift zones from east to west [72]. In addition, the NS-trending faults that are considered as the result of east-west extension of the Qinghai-Tibet Plateau [73, 74] also formed during this period especially from 18 to 4 Ma [75, 76]. These NS-trending faults are also the important ore-controlling structures in the NH [1].Magmatism in the NH primarily includes the Mesozoic and Cenozoic magmatism. The Mesozoic magmatism generated multiple suites of mafic-intermediate igneous rocks between the Late Triassic and the Early Cretaceous, including basaltic volcanic interlayers, dyke swarms, and subvolcanic dykes. According to the previous geochronological data, the SHRIMP U-Pb ages of the basic dyke swarms from different area in the NH are 134.9 ± 1.8 Ma, 135.5 ± 2.1 Ma [77], and 138.0 ± 3.5 Ma [78], respectively. The SHRIMP U-Pb age of the gabbro is 155.8 Ma [79]. Tong et al. [78], Pan et al. [80], and Zhong et al. [81] regarded these basic dyke swarms as the result of late-stage massive expansion of Neo-Tethys Ocean under the structural environment of the Himalaya passive continental margin intensive stretching and breaking-off, lithosphere extension-thinning, and asthenosphere upwelling. On the contrary, Zhu et al. [82] and Qiu et al. [83] suggested that these basic dyke swarms are the result of interaction between mantle plume and lithospheric mantle material and form in the continental-rift environment. The Cenozoic magmatism is characterized by the formation of monzogranite, leucogranite, diorite, porphyritic diorite, and aplite units [84, 85]. These Cenozoic intermediate-acidic intrusive masses are distributed along the EW-trending faults and in the core of Ranba, Kangma, and Yelaxiangbo dome in the form of batholith, laccolith, and dykes (Figure 1(b)). This phenomenon is considered to be the result of crustal thickening [86] related to the collision of the India Plate and the Eurasia Plate during the postcollision stage (25 to 0 Ma) [87, 88].The NHMB contains many Sb, Au, Sb-Au, Pb-Zn, and Sb-Pb-Zn-Ag deposits, and the Zhaxikang Sb-Pb-Zn-Ag, the Mazhala Au-Sb, the Chalapu Au, the Bangbu Au, the Shalagang Sb, the Cheqiongzhuobu Sb deposits are representative (Figure1(b)) [7, 25]. The geneses and metallogenic age of these deposits are controversial due to the limited research, the genetic models mainly include the SEDEX overprinted by hot spring, carlin and carlin-like, hot spring, subvolcanic magmatic-hydrothermal fluid, and orogenic types [25].


## 2.2. Ore Deposit Geology

The Zhaxikang Sb-Pb-Zn-Ag deposit is located ~48 km west from Longzi County Town within the southeastern Yangzuoyong-Nariyong composite syncline in the NH (Figure1(b)). This deposit has a reserve of 1.268 Mt Pb-Zn with a 3.66% average Zn grade and a 2.45% average Pb grade, 0.138 Mt Sb with an average grade of 1.08%, 1800 t Ag with an average of 99.55 g/t, 3.9 t associated Au, 361 t associated Ga, and 20 Mt Mn-Fe carbonate ores with an average grade of 42% for Fe + Mn [89], which makes it the largest deposit within the NHMB. The majority of mineralization in the orefield is hosted by the Lower Jurassic Ridang formation that consists of epimetamorphic marine clastic rocks. This formation, dipping shallowly to the north and striking eastwest, is divided into five lithologic sections (Figure 2(a)). A few Upper Jurassic Weimei formations composed of fine-grained metamorphic quartzose sandstone, silty slate, and calcarenite as well as Quaternary sediments distributed along valleys also crop out in the orefield (Figure 2(a)) [7].Figure 2
(a) Geological map of the Zhaxikang Sb-Pb-Zn-Ag polymetallic deposit (modified from Zheng et al. [7]). (b) Cross-section along Exploration Line 7. J
1
r
1 to J
1
r
5: the first to fifth lithologic section of the Lower Jurassic Ridang formation.
(a)
(b)The Zhaxikang deposit developed extensive geological structures. A near northsouth striking fault system is prevalent in the orefield, which coexists with a group of northeast-striking faults and some folds. Engineering and geological mapping projects have identified 16 faults, the majority of which are steeply dipping normal faults associated with tensional stress and wrench faults associated with torsional stress. Faults F2, F4, F5, F6, F7, F13, F14, and F16 are the main ore-bearing faults, faults F1 and F10 are partly mineralized, fault F3 was associated with late-stage mineralization, faults F8 and F9 are wrench faults without any mineralization, and fault F12 is a nonmineralized regional fault (Figure2(a)). The orebodies I–VI are hosted by nearly NS-striking faults and orebodies VII–IX are hosted by nearly NE-striking faults (Figure 2). Our samples in this study are all from the orebody V, which is the largest and richest one among these orebodies within the orefield and hosts more than 80% of the reserves. This orebody is >1400 m long, 1 to 30 m wide, and controlled by fault F7 (Figure 2).The magmatism in the orefield is associated with diabase, porphyritic rhyolite, basalt, and leucogranite units as well as some granite porphyry dykes that intruded into the porphyritic rhyolite (Figure2(a)). The diabase is identified by drillholes and footrill in the central part of the orefield as dykes that emplaced into the Jurassic Ridang Formation and has been dated by zircon U-Pb methods to ~133 Ma [7]. The rhyolite porphyry with the zircon SHRIMP U-Pb age of ~135 Ma crops out in the western part of the orefield [90] and the leucogranites crop out in the southern part of the orefield over an area of <1 km2. Additionally, the basalt usually occurs near the orebody in the form of consequent layer or shear layer distributed in slate and the contact region of slate and diabase.Various types of alteration associated with mineralization have occurred in the orefield, including (1) the silicification that is associated with Sb mineralization and generally located in fault zones in the form of quartz veins, radiating quartz, and quartz druse; (2) the carbonatization that is associated with Pb-Zn mineralization in the form of Mn-Fe carbonate veins and also formed the postmineralization calcite; (3) the chlorite alteration that is generally confined to massive and stellated aggregates of chlorite within diabase; (4) the weak sericite alteration that is associated with chlorite alteration and barren quartz; and (5) the clay alteration that developed along the edges of fracture-related crushed zones. Furthermore, the ore-forming elements display a vertical sequence that is zoned from a lowermost Zn (Pb + Ag) zone through a central Zn + Pb + Ag-(Sb) zone to an uppermost Pb + Zn + Sb + Ag zone, although no horizontal zoning is present [89].


## 2.3. Ore Paragenetic Sequence

The paragenetic sequence in the Zhaxikang deposit is divided into six stages of ore formation based on the detailed hand specimen and microscopic observations. These six stages are assigned to two clear pulses: the first pulse consists of stages 1 to 2 and is characterized by the assemblages of Mn-Fe carbonates and sulfides, and the second pulse includes stages 3 to 6 and is primarily dominated by quartz, calcite, sulfosalt minerals, and sulfides (Figure3).Figure 3
Mineral paragenesis within the Zhaxikang Sb-Pb-Zn-Ag deposit (modified from Wang et al. [26]).Stage 1, dominated by a Mn-Fe carbonate + sphalerite + pyrite + arsenopyrite assemblage, is the initial stage of ore formation in the Zhaxikang deposit. Majority of the fine-grained sphalerite, pyrite, and arsenopyrite are hosted by fine-grained Mn-Fe carbonate in the form of laminae (Figures4(a)–4(c)), and a few sulfides occur within the Mn-Fe carbonates as stellated aggregates (Figure 4(d)). The fine-grained layered and colloform with synsedimentary features (Figure 4(b)). The laminae and Mn-Fe carbonates in some samples have been cut by later stage 4 quartz-boulangerite veins (Figure 4(c)) or have been affected by the stage 2 coarse-grained sphalerite (Figure 4(a)).Figure 4
Hand specimen photographs and photomicrographs of representative samples from the Zhaxikang deposit. (a) Stage 1 lamellar sphalerite-pyrite-arsenopyrite and stage 2 massive sphalerite-pyrite hosted within fine-grained Mn-Fe carbonate. (b) Stage 1 lamellar and stage 2 banded Mn-Fe carbonate-sphalerite-galena ore with visible synsedimentary features including rhythmic sedimentation in the upper part and angular folding in the lower part of the sample. (c) Stage 1 lamellar sphalerite-pyrite-arsenopyrite and stage 2 massive and banded sphalerite-pyrite hosted by fine-grained Mn-Fe carbonate. The mineral assemblage is in turn cross-cut by stage 4 quartz-boulangerite veins. (d) Coarse-grained stage 2 Mn-Fe carbonate-sphalerite formed by the recrystallization of fine-grained stage 1 Mn-Fe carbonate-sphalerite. (e) Stage 2 massive coarse-grained pyrite hosted by slate. (f) Massive and brecciated stage 2 sphalerite hosted in stage 2 Mn-Fe carbonate. (g) Stage 2 massive, globular, and concentric annular sphalerite-pyrite hosted by coarse-grained Mn-Fe carbonate. (h) Stage 2 massive galena and pyrite. (i) Stage 2 massive and veined sphalerite-pyrite hosted by coarse-grained Mn-Fe carbonate. (j) Stage 1 lamellar sphalerite-pyrite-arsenopyrite and stage 2 massive sphalerite hosted by fine-grained Mn-Fe carbonate. The sample also contains a Mn-Fe carbonate druse dominated by idiomorphic columnar quartz and valentinite. (k) Stage 2 coarse-grained sphalerite-pyrite hosted by stage 2 Mn-Fe carbonate with banded textures. (l) Stage 2 sphalerite and Mn-Fe carbonate ore with typical Dal Matianite texture. (m) Disseminated stage 3 sphalerite and pyrite hosted in stage 3 quartz and calcite, and stage 3 quartz and calcite cut the slate. (n) Stage 3 brecciated sphalerite within stage 3 quartz-calcite. (o) Stage 3 sphalerite-galena veins cross-cut by stage 6 quartz-calcite veins. (p) Stage 3 sphalerite cross-cut by stage 4 quartz-boulangerite veins. (q) Stage 4 massive and needle-like boulangerite hosted by stage 4 quartz. (r) Stage 4 boulangerite-quartz. (s) Stage 5 elongate stibnite hosted by stage 5 quartz. (t) Stage 5 stibnite-cinnabar hosted by stage 5 quartz. (u) Siliceous sinter formed during the Supergene stage. (v) Ferrihydrite and sardinianite formed during the Supergene stage. (w) Stage 2 pyrite containing automorphic stage 1 arsenopyrite is replaced by later stage 2 sphalerite to form a skeletal texture. (x) Stage 3 sphalerite occurs in stage 3 quartz. (y) The emulsion-like and disseminated stage 3 chalcopyrite grains are dotted in the stage 3 sphalerite. (z) The stage 3 chalcopyrite grains are dotted among the grains of stage 3 pyrite. (aa) Stage 3 sphalerite replaced by stage 4 boulangerite that is in turn cross-cut by stage 5 stibnite. (ab) Stage 5 stibnite cross-cut by stage 6 quartz. Mcar1 = stage 1 fine-grained Mn-Fe carbonate; Apy1 = stage 1 lamellar arsenopyrite; Py1 = stage 1 lamellar pyrite; Sp1 = stage 1 lamellar sphalerite; Mcar2 = stage 2 coarse-grained Mn-Fe carbonate; Apy2 = stage 2 arsenopyrite; Py2 = stage 2 pyrite; Sp2 = stage 2 sphalerite; Gn2 = stage 2 coarse-grained galena; Py3 = stage 3 pyrite; Sp3 = stage 3 sphalerite; Gn3 = stage 3 galena; Ccp3 = stage 3 chalcopyrite; Qtz3 = stage 3 quartz; Cal3 = stage 3 calcite; Blr4 = stage 4 boulangerite; Qtz4 = stage 4 quartz; Stb5 = stage 5 stibnite; Ci5 = stage 5 cinnabar; Qtz5 = stage 5 quartz; Cal6 = stage 6 calcite; Qtz6 = stage 6 quartz; Lm = Supergene stage ferrihydrite.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)
(q)
(r)
(s)
(t)
(u)
(v)
(w)
(x)
(y)
(z)
(aa)
(ab)Stage 2 hosts majority of the Pb-Zn mineralization in the Zhaxikang deposit and comprises an assemblage of Mn-Fe carbonate + galena + sphalerite + pyrite ± arsenopyrite. The more abundant and coarser-grained sulfides are hosted by coarse-grained Mn-Fe carbonate and slate to form the banded (Figure4(k)), net-veined (Figure 4(i)), massive (Figure 4(h)), concentric annular (Figures 4(g) and 5), and Dal Matianite (Figure 4(l)) ores. The Mn-Fe carbonate during this stage recrystallized to different degree (Figures 4(f), 4(g) and 4(i)–4(k)), some even formed the druse containing the idiomorphic columnar quartz, needle-like boulangerite, or valentinite (Figure 4(j)). We can also observe that the later sphalerite replaces the earlier pyrite containing automorphic stage 1 arsenopyrite to form a skeletal texture during this stage (Figure 4(w)). Zheng et al. [7] considered that the ore textures in stages 1 and 2 are similar to those of the Red Dog SEDEX-type ore district in Alaska.Figure 5
The sampling points and Fe-Zn isotopic-elemental variations in the annular polished section samples: (a) 9-8; (b) 9-3; (c) ZXK-1; (d) ZXK-2. Mcar = Mn-Fe carbonate; Apy = arsenopyrite; Py = pyrite; Sp = sphalerite; Qtz = quartz.
(a)
(b)
(c)
(d)Stage 3, characterized by the formation of a quartz ± calcite + pyrite + sphalerite + galena ± chalcopyrite ± arsenopyrite assemblage without Mn-Fe carbonate, is the earliest stage of the second pulse of mineralization. The massive, veined, net-veined, and brecciated sphalerite, galena, and pyrite occur in the quartz and calcite (Figures4(l)–4(p), 4(x)), and most of the sulfides form by the modification of sulfides from earlier stages. Some chalcopyrite grains are distributed in sphalerite, galena, and pyrite as stellated aggregates (Figures 4(y) and 4(z)). Some of these sulfides have been cross-cut by the later quartz-boulangerite or quartz-calcite veins (Figure 4(o)).Stage 4 is marked by a mineral assemblage composed of quartz + antimony-lead-silver sulfosalt minerals (Figure4(r)) that prevailingly include boulangerite and jamesonite, as well as minor bournonite, tetrahedrite, and andorite. This stage hosts the majority of the Sb and Ag mineralization and yields the ores with relatively high average Ag grades. The minerals formed in earlier stages are replaced and cross-cut by the quartz-boulangerite veins and boulangerite of this stage (Figures 4(c) and 4(aa)). Some samples also contain quartz druse filled with needle-like boulangerite (Figure 4(q)).Stage 5 is distinguished by the formation of a quartz + stibnite + cinnabar assemblage and hosts part of the Sb mineralization within the deposit. Elongate-radial stibnite and massive stibnite-cinnabar are hosted by quartz (Figures4(s) and 4(t)). Some stibnite cross-cut the stage 4 boulangerite (Figure 4(aa)).Stage 6, representing the youngest stage of mineralization in the Zhaxikang deposit, is identified by a quartz ± calcite assemblage without sulfides. The quartz-calcite veins of this stage cross-cut the earlier formed minerals (Figures4(o) and 4(ab)). Zheng et al. [7] regarded the ore textures in the second pulse of mineralization as typical hot spring type metallogenic features.Supergene stage in the Zhaxikang deposit consists of ferrihydrite, smithsonite, sardinianite, valentinite, travertine, malachite, and siliceous sinter (Figures4(u) and 4(v)).


## 3. Sampling and Analytical Methods

### 3.1. Sampling

The sampling points for EPMA and Fe-Zn isotopic analyses (the powders are sampled by microdrill) are all in the annular polished section samples 9-3, 9-8, ZXK-1, and ZXK-2. The specific numbers, locations, and photomicrographs of these sampling points are given in Figures5 and 6, respectively. These four samples are all from the first pulse of mineralization, only the sample ZXK-2 has been cut by a stage 3 quartz vein (Figure 5(d)).Figure 6
The photomicrographs and electron probe micrographs of sampling area in the annular polished section samples for EPMA and Fe-Zn isotopic analyses. (a) 9-3-2; (b) 9-8-8; (c) 9-8-2; (d) 9-8-3 and 9-8-5; (e) ZXK-1-4; (f) ZXK-1-3; (g) ZXK-1-2; (h) ZXK-2-1; (i) ZXK-2-5; (j) ZXK-2-2; (k) ZXK-2-3; (l) ZXK-2-4; (m) ZXK-1-3; (n) ZXK-2-5; (o) ZXK-2-2; (p) ZXK-2-3. Abbreviations are as in Figure5.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)


### 3.2. EPMA

Chemical compositions of sulfide, Mn-Fe carbonate, and quartz were determined on a JEOL (Japan Electron Optics Laboratory) JXA-8100 electron microprobe (EMP) at the Second Institute of Oceanography, State Oceanic Administration of China. The accelerating voltage is 15 kV for Mn-Fe carbonate and quartz and 20 kV for sulfide, the beam current is 10 nA, the beam diameter is 1μm, the secondary electronic resolution is 6 nm with the operating distance of 11 mm, and the repeat accuracy of the sample stage is within 1 nm. The standards are natural minerals and synthetic oxides as those of Sun et al. [2]. The correction program supplied by the manufacturer is used for matrix corrections [91, 92].


### 3.3. Fe-Zn Isotopic Analyses

Approximately 10–50 milligrams of sample powders was placed in 15 ml Teflon jars and the solids were dissolved in 4 ml of heated ultrapure aqua regia. The solutions were dried and then Fe and Zn were purified using the BioRad MP-1 anion exchange resin using the protocol from Maréchal et al. [14]. Yields from the columns were tested volumetrically on the ICP-OES at Pennsylvania State University and were all greater than 95%. Isotope values are reported in the traditional per mil values (‰).The Fe isotopes were measured on the Neptune MC-ICP-MS at Pennsylvania State University. The instrument setup, sample introduction, and running conditions are discussed in greater detail in Yesavage et al. [93]. Samples were diluted to a 3 ppm Fe solution which produced approximately a 10 V signal on the shoulder to the argon interference peak (56Fe and A
r
40
16O). Sample intensities matched the intensity of the bracketing standard within 10%. Mass bias was corrected for by standard-sample-standard bracketing. In-house and international standards were measured throughout the sessions and yielded overlapping values of SRM-3126a δ
56Fe = 0.33 ± 0.08‰, n = 8 (accepted values δ
56Fe = 0.34 ± 0.1‰ 2σ) [93], and HPS-WU δ
56Fe = 0.62 ± 0.11‰, n = 8 (accepted values δ
56Fe = 0.60 ± 0.07‰ 2σ) [34]. The samples are reported relative to the international standard IRMM-014 (δ
56Fe (‰) = [
(
F
e
56
/
F
e
54
)
s
a
m
p
l
e
/
(
F
e
56
/
F
e
54
)
I
R
M
M
–
014
-
1
]
×
1000). Reported values are an average of two different measurements and the errors fall within the range 0.1‰ 2σ of the standards.The Zn isotopes were measured on Neptune MC-ICP-MS at Rutgers University. Correction of mass bias for Zn using Cu (NIST 976) was employed for these samples as suggested in [94–97] and the corrected values were then bracketed by the standards. The samples are reported relative to the newly developed Zn isotope standard (AA-ETH; δ
66Zn (‰) = [(
Z
n
66
/
Z
n
64
)
s
a
m
p
l
e
/
(
Z
n
66
/
Z
n
64
)
A
A
–
E
T
H
-
1] × 1000) and all the quoted data from previous literatures in this paper are converted relative to the AA-ETH standard (δ
66
Z
n
A
A
–
E
T
H
=
δ
66
Z
n
J
M
C
3
–
0749
L
-
0.28
‰) [98]. We also compared the new standard relative to IRMM 3702 and obtained a δ
66Zn = 0.03‰, which is within the error reported in Archer et al. [98]. Solutions were kept at 100 ppb Cu and 150 ppb Zn which generated 63Cu = 7 V and 66Zn = 4 V. One block of 30 ratios is reported and the average error for the standard compared to itself throughout the session is 0.05‰ 2σ.


## 3.1. Sampling

The sampling points for EPMA and Fe-Zn isotopic analyses (the powders are sampled by microdrill) are all in the annular polished section samples 9-3, 9-8, ZXK-1, and ZXK-2. The specific numbers, locations, and photomicrographs of these sampling points are given in Figures5 and 6, respectively. These four samples are all from the first pulse of mineralization, only the sample ZXK-2 has been cut by a stage 3 quartz vein (Figure 5(d)).Figure 6
The photomicrographs and electron probe micrographs of sampling area in the annular polished section samples for EPMA and Fe-Zn isotopic analyses. (a) 9-3-2; (b) 9-8-8; (c) 9-8-2; (d) 9-8-3 and 9-8-5; (e) ZXK-1-4; (f) ZXK-1-3; (g) ZXK-1-2; (h) ZXK-2-1; (i) ZXK-2-5; (j) ZXK-2-2; (k) ZXK-2-3; (l) ZXK-2-4; (m) ZXK-1-3; (n) ZXK-2-5; (o) ZXK-2-2; (p) ZXK-2-3. Abbreviations are as in Figure5.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)


## 3.2. EPMA

Chemical compositions of sulfide, Mn-Fe carbonate, and quartz were determined on a JEOL (Japan Electron Optics Laboratory) JXA-8100 electron microprobe (EMP) at the Second Institute of Oceanography, State Oceanic Administration of China. The accelerating voltage is 15 kV for Mn-Fe carbonate and quartz and 20 kV for sulfide, the beam current is 10 nA, the beam diameter is 1μm, the secondary electronic resolution is 6 nm with the operating distance of 11 mm, and the repeat accuracy of the sample stage is within 1 nm. The standards are natural minerals and synthetic oxides as those of Sun et al. [2]. The correction program supplied by the manufacturer is used for matrix corrections [91, 92].


## 3.3. Fe-Zn Isotopic Analyses

Approximately 10–50 milligrams of sample powders was placed in 15 ml Teflon jars and the solids were dissolved in 4 ml of heated ultrapure aqua regia. The solutions were dried and then Fe and Zn were purified using the BioRad MP-1 anion exchange resin using the protocol from Maréchal et al. [14]. Yields from the columns were tested volumetrically on the ICP-OES at Pennsylvania State University and were all greater than 95%. Isotope values are reported in the traditional per mil values (‰).The Fe isotopes were measured on the Neptune MC-ICP-MS at Pennsylvania State University. The instrument setup, sample introduction, and running conditions are discussed in greater detail in Yesavage et al. [93]. Samples were diluted to a 3 ppm Fe solution which produced approximately a 10 V signal on the shoulder to the argon interference peak (56Fe and A
r
40
16O). Sample intensities matched the intensity of the bracketing standard within 10%. Mass bias was corrected for by standard-sample-standard bracketing. In-house and international standards were measured throughout the sessions and yielded overlapping values of SRM-3126a δ
56Fe = 0.33 ± 0.08‰, n = 8 (accepted values δ
56Fe = 0.34 ± 0.1‰ 2σ) [93], and HPS-WU δ
56Fe = 0.62 ± 0.11‰, n = 8 (accepted values δ
56Fe = 0.60 ± 0.07‰ 2σ) [34]. The samples are reported relative to the international standard IRMM-014 (δ
56Fe (‰) = [
(
F
e
56
/
F
e
54
)
s
a
m
p
l
e
/
(
F
e
56
/
F
e
54
)
I
R
M
M
–
014
-
1
]
×
1000). Reported values are an average of two different measurements and the errors fall within the range 0.1‰ 2σ of the standards.The Zn isotopes were measured on Neptune MC-ICP-MS at Rutgers University. Correction of mass bias for Zn using Cu (NIST 976) was employed for these samples as suggested in [94–97] and the corrected values were then bracketed by the standards. The samples are reported relative to the newly developed Zn isotope standard (AA-ETH; δ
66Zn (‰) = [(
Z
n
66
/
Z
n
64
)
s
a
m
p
l
e
/
(
Z
n
66
/
Z
n
64
)
A
A
–
E
T
H
-
1] × 1000) and all the quoted data from previous literatures in this paper are converted relative to the AA-ETH standard (δ
66
Z
n
A
A
–
E
T
H
=
δ
66
Z
n
J
M
C
3
–
0749
L
-
0.28
‰) [98]. We also compared the new standard relative to IRMM 3702 and obtained a δ
66Zn = 0.03‰, which is within the error reported in Archer et al. [98]. Solutions were kept at 100 ppb Cu and 150 ppb Zn which generated 63Cu = 7 V and 66Zn = 4 V. One block of 30 ratios is reported and the average error for the standard compared to itself throughout the session is 0.05‰ 2σ.


## 4. Results

### 4.1. EPMA

All the EPMA data are given in Tables1 and 2. The Mn-Fe carbonate contains 23.562~31.806 wt% Fe and 27.514~32.232 wt% Mn, with a negative correlation between Fe and Mn contents (Figure 7(a)), which indicates that the Mn-Fe carbonates form by the isomorphic substitution of Fe2+ and Mn2+ ions and have a molecular formula of (Mn0.5Fe0.5)CO3. The Fe contents are around 36 wt% for the arsenopyrite samples and 46 wt% for the pyrite samples. The sphalerite samples have 56.941~60.552 wt% Zn and 5.375~9.424 wt% Fe with a negative correlation between these two elements (Figure 7(b)).Table 1
EPMA data for Mn-Fe carbonate and quartz from the Zhaxikang deposit.
Sample number
SiO2 (%)
Na2O (%)
K2O (%)
TiO2 (%)
Al2O3 (%)
MgO (%)
CaO (%)
FeO (%)
MnO (%)
SrO (%)
BaO (%)
CO2 (%)
Total (%)
Mineral
9-3-4
0.015
0.005
0.007
0.011
0.003
1.231
2.061
28.598
29.840
0.007
0.057
38.164
99.999
Mcar
9-3-5
0.024
0.109
0.027
–
0.010
1.080
2.151
28.372
30.582
–
0.002
37.643
100.000
Mcar
9-8-8
0.002
0.031
0.007
–
–
1.276
1.978
27.341
31.338
0.010
0.022
37.993
99.998
Mcar
9-8-9
–
0.025
–
–
0.008
1.189
1.144
31.806
27.514
–
–
38.316
100.002
Mcar
ZXK-1-6
0.001
0.028
0.001
0.022
–
1.207
5.178
24.070
30.717
–
–
38.775
99.999
Mcar
ZXK-1-7
0.028
0.011
0.019
0.036
0.014
1.995
2.940
29.639
27.943
–
0.017
37.358
100.000
Mcar
ZXK-2-7
0.005
0.039
–
0.055
–
1.236
5.542
27.235
27.604
–
–
38.285
100.001
Mcar
ZXK-2-8
–
0.051
–
0.012
0.010
1.140
4.298
24.287
32.694
–
0.008
37.500
100.000
Mcar
ZXK-2-9
0.009
0.049
0.008
–
0.037
1.552
1.946
26.739
30.968
0.007
0.013
38.673
100.001
Mcar
ZXK-2-10
0.015
–
–
0.023
0.012
1.119
4.849
23.562
32.232
–
0.003
38.185
100.000
Mcar
ZXK-2-11
97.000
0.026
0.025
–
0.219
0.028
0.018
0.020
—
–
–
2.664
100.000
Qtz
Notes. (1) Mcar = Mn-Fe carbonate; Qtz = quartz; (2) for the sampling points, please see Figure 5.Table 2
EPMA data for sulfides from the Zhaxikang deposit.
Sample number
As (%)
S (%)
Se (%)
Cu (%)
Co (%)
Zn (%)
Sb (%)
Ag (%)
Au (%)
Cd (%)
Bi (%)
Pb (%)
NiO (%)
FeO (%)
Total (%)
Mineral
9-3-1
–
33.165
–
0.133
0.009
59.095
–
–
–
0.171
–
–
–
6.985
99.558
Sp
9-3-2
–
33.423
–
0.015
0.012
58.224
–
0.002
–
0.151
–
–
–
8.060
99.887
Sp
9-3-3
–
33.262
–
–
0.015
57.114
–
–
–
0.159
–
–
–
9.405
99.955
Sp
9-3-6
0.784
52.884
–
0.005
0.044
–
–
–
0.016
–
–
0.134
–
46.585
100.452
Py
9-3-7
2.244
51.934
–
0.002
0.049
–
–
–
0.004
–
–
0.132
–
46.044
100.409
Py
9-8-1
0.002
33.740
–
0.009
0.010
57.705
–
0.004
–
0.160
–
0.018
0.008
8.827
100.483
Sp
9-8-2
–
33.649
–
0.011
0.005
59.584
–
0.005
–
0.150
–
–
–
7.112
100.516
Sp
9-8-3
–
33.500
–
0.006
0.003
56.941
–
–
–
0.160
–
–
0.002
9.424
100.036
Sp
9-8-5
0.009
34.064
–
0.089
0.004
57.475
–
–
–
0.151
–
–
–
8.901
100.693
Sp
9-8-6
0.006
53.867
–
0.033
0.050
0.058
–
0.005
0.012
–
–
0.300
0.002
46.539
100.872
Py
9-8-7
0.015
53.359
–
0.471
0.043
0.044
–
–
–
–
–
0.174
–
46.518
100.624
Py
ZXK-1-1
–
33.502
–
0.012
0.009
57.497
–
–
–
0.172
–
–
–
9.142
100.334
Sp
ZXK-1-2
–
33.428
–
0.036
–
60.190
0.002
0.002
–
0.201
–
–
–
5.810
99.669
Sp
ZXK-1-3(1)
39.757
24.131
–
0.006
0.043
0.009
–
–
0.027
–
–
0.063
0.002
36.340
100.378
Apy
ZXK-1-3(2)
2.435
51.666
–
0.018
0.067
0.185
–
0.007
0.013
–
–
0.113
0.004
45.877
100.385
Py
ZXK-1-3(3)
–
33.337
0.002
0.014
0.012
58.623
–
0.015
–
0.147
–
–
0.003
7.747
99.900
Sp
ZXK-1-4
1.942
51.982
0.001
0.008
0.046
0.015
–
–
0.003
–
–
0.196
–
46.004
100.197
Py
ZXK-1-5
2.056
52.248
–
0.008
0.054
–
–
–
0.013
–
–
0.137
–
46.163
100.679
Py
ZXK-2-1
–
33.684
–
0.061
0.010
59.726
–
–
–
0.170
–
–
0.005
6.904
100.560
Sp
ZXK-2-2
39.643
23.789
–
0.002
0.050
0.008
0.196
0.006
–
–
–
0.081
0.098
36.245
100.118
Apy
ZXK-2-3
–
33.868
–
0.025
0.012
58.034
–
0.001
–
0.196
–
–
0.002
8.714
100.852
Sp
ZXK-2-4
39.444
24.229
–
–
0.041
–
0.173
–
0.026
–
–
0.061
–
36.424
100.398
Apy
ZXK-2-5
2.417
51.755
–
–
0.045
0.027
–
–
–
–
–
0.131
–
46.227
100.602
Py
ZXK-2-6
0.031
53.630
–
0.107
0.047
–
–
–
0.009
–
–
0.168
–
46.654
100.646
Py
Notes. (1) Py = pyrite, Sp = sphalerite, and Apy = arsenopyrite; (2) for the sampling points, please see Figure 5.Figure 7
(a) The negative correlation between Mn and Fe contents within Mn-Fe carbonate. (b) The negative correlation between Zn and Fe contents within sphalerite. (c) The negative correlation betweenδ
66Zn values and Fe contents within sphalerite from 9-3, 9-8, and ZXK-1. (d) The positive correlation between δ
66Zn values and Zn contents within sphalerite from 9-3, 9-8, and ZXK-1.
(a)
(b)
(c)
(d)


### 4.2. Fe-Zn Isotopes

All the Fe-Zn isotopic data are given in Table3. The annular polished section samples have δ
56
F
e
I
R
M
M
-
014 of -
1.95
‰
​
​
~
0.43
‰, with an average of −0.50‰  ± 1.09‰ (2SD, n = 19), and δ
66
Z
n
A
A
-
E
T
H of -
0.38
‰
​
~
0.07
‰ with an average of −0.25‰  ± 0.19‰ (2SD, n = 12). The Mn-Fe carbonate and pyrite show the δ
56Fe values range from −1.95‰ to −0.59‰ (average value of −0.97‰  ± 0.86‰; 2SD, n = 8) and from −0.26‰ to 0.23‰ (average value of −0.07‰  ± 0.35‰; 2SD, n = 9), respectively. The two arsenopyrite samples exhibit the δ
56Fe values of −1.01‰ (ZXK-1-3; Figure 5(c)) and −0.18‰ (ZXK-2-2; Figure 5(d)).Table 3
Fe-Zn isotopic data for annular polished section samples from the Zhaxikang deposit.
Sample number
Mineral
δ
56
F
e
I
R
M
M
–
014
2σ
δ
66
Z
n
A
A
–
E
T
H
2σ
9-3-1
Sp
−0.12
0.05
9-3-2
Sp
−0.23
0.05
9-3-3
Sp
−0.31
0.05
9-3-4
Mcar
−0.59
0.08
9-3-5
Mcar
−1.95
0.08
9-3-7
Py
−0.26
0.08
9-8-1
Sp
−0.17
0.05
9-8-2
Sp
−0.09
0.05
9-8-3
Sp
−0.35
0.05
9-8-4
Py
0.23
0.08
9-8-5
Sp
−0.23
0.05
9-8-6
Py
−0.26
0.08
9-8-7
Py
0.09
0.08
9-8-8
Mcar
−0.78
0.08
9-8-9
Mcar
−0.92
0.08
ZXK-1-1
Sp
−0.38
0.05
ZXK-1-2
Sp
−0.15
0.05
ZXK-1-3(1)
Apy
−1.01
0.08
ZXK-1-3(2)
Py
−0.22
0.08
ZXK-1-3(3)
Sp
−0.18
0.05
ZXK-1-4
Py
0.12
0.08
ZXK-1-5
Py
−0.05
0.08
ZXK-1-6
Mcar
−0.66
0.09
ZXK-1-7
Mcar
−0.78
0.08
ZXK-2-1
Sp
−0.32
0.05
ZXK-2-2
Apy
−0.18
0.08
ZXK-2-3
Sp
−0.28
0.05
ZXK-2-5
Py
−0.07
0.08
ZXK-2-6
Py
−0.17
0.08
ZXK-2-7
Mcar
−1.06
0.08
ZXK-2-8
Mcar
−1.04
0.08
Notes. (1) Abbreviations and sampling points are as in Tables 1 and 2; (2) 2σ is two times the standard deviation.In sample 9-3, the 9-3-1, 9-3-2 (Figure6(a)), and 9-3-3 sphalerite have the δ
66Zn values of −0.12‰, −0.23‰, and −0.31‰ with a gradually decreasing trend, the 9-3-4 and 9-3-5 Mn-Fe carbonate show the same trend with the δ
56Fe values of −0.59‰ and −1.95‰, and the δ
56Fe value of 9-3-7 pyrite is −0.26‰ (Figure 5(b)). Similarly, in sample 9-8, from 9-8-2 (−0.09‰; Figure 6(c)) through 9-8-5 (−0.23‰) to 9-8-3 (−0.35‰; Figure 6(d)), the δ
66Zn values of sphalerite also present a gradually decreasing trend; meanwhile the δ
56Fe values also decrease from −1.06‰ (9-8-8; Figure 6(b)) to −1.38‰ (9-8-9) for Mn-Fe carbonate and from 0.23‰ (9-8-4) to 0.09‰ (9-8-7) for pyrite, although the 9-8-1 sphalerite and 9-8-6 pyrite have the δ
66Zn value of −0.17‰ and δ
56Fe values of −0.26‰ (Figure 5(a)). In sample ZXK-1, the δ
66Zn values of ZXK-1-3 (−0.18‰; Figures 6(f) and 6(m)) and ZXK-1-2 (−0.15‰; Figure 6(g)) sphalerite are almost the same and heavier than ZXK-1-1 (−0.38‰) sphalerite; the δ
56Fe values for pyrite increase from ZXK-1-3 (−0.22‰) to ZXK-1-4 (0.12‰; Figure 6(e)) and ZXK-1-5 (−0.05‰); and, as for Mn-Fe carbonate, the ZXK-1-6 (−0.66‰) has heavier δ
56Fe values than ZXK-1-7 (−0.78‰; Figure 5(c)). By contrast, in sample ZXK-2, the ZXK-2-1 (−0.32‰; Figure 6(h)) and ZXK-2-3 (−0.28‰; Figures 6(k) and 6(p)) sphalerite have almost the same δ
66Zn values; the ZXK-2-7 (−1.06‰) and ZXK-2-8 (−1.04‰) Mn-Fe carbonate also yield approximate δ
56Fe values; also the ZXK-2-5 (Figures 6(i) and 6(n)) and ZXK-2-6 pyrite show the similar δ
56Fe values of −0.07‰ and −0.17‰ (Figure 5(d)).


## 4.1. EPMA

All the EPMA data are given in Tables1 and 2. The Mn-Fe carbonate contains 23.562~31.806 wt% Fe and 27.514~32.232 wt% Mn, with a negative correlation between Fe and Mn contents (Figure 7(a)), which indicates that the Mn-Fe carbonates form by the isomorphic substitution of Fe2+ and Mn2+ ions and have a molecular formula of (Mn0.5Fe0.5)CO3. The Fe contents are around 36 wt% for the arsenopyrite samples and 46 wt% for the pyrite samples. The sphalerite samples have 56.941~60.552 wt% Zn and 5.375~9.424 wt% Fe with a negative correlation between these two elements (Figure 7(b)).Table 1
EPMA data for Mn-Fe carbonate and quartz from the Zhaxikang deposit.
Sample number
SiO2 (%)
Na2O (%)
K2O (%)
TiO2 (%)
Al2O3 (%)
MgO (%)
CaO (%)
FeO (%)
MnO (%)
SrO (%)
BaO (%)
CO2 (%)
Total (%)
Mineral
9-3-4
0.015
0.005
0.007
0.011
0.003
1.231
2.061
28.598
29.840
0.007
0.057
38.164
99.999
Mcar
9-3-5
0.024
0.109
0.027
–
0.010
1.080
2.151
28.372
30.582
–
0.002
37.643
100.000
Mcar
9-8-8
0.002
0.031
0.007
–
–
1.276
1.978
27.341
31.338
0.010
0.022
37.993
99.998
Mcar
9-8-9
–
0.025
–
–
0.008
1.189
1.144
31.806
27.514
–
–
38.316
100.002
Mcar
ZXK-1-6
0.001
0.028
0.001
0.022
–
1.207
5.178
24.070
30.717
–
–
38.775
99.999
Mcar
ZXK-1-7
0.028
0.011
0.019
0.036
0.014
1.995
2.940
29.639
27.943
–
0.017
37.358
100.000
Mcar
ZXK-2-7
0.005
0.039
–
0.055
–
1.236
5.542
27.235
27.604
–
–
38.285
100.001
Mcar
ZXK-2-8
–
0.051
–
0.012
0.010
1.140
4.298
24.287
32.694
–
0.008
37.500
100.000
Mcar
ZXK-2-9
0.009
0.049
0.008
–
0.037
1.552
1.946
26.739
30.968
0.007
0.013
38.673
100.001
Mcar
ZXK-2-10
0.015
–
–
0.023
0.012
1.119
4.849
23.562
32.232
–
0.003
38.185
100.000
Mcar
ZXK-2-11
97.000
0.026
0.025
–
0.219
0.028
0.018
0.020
—
–
–
2.664
100.000
Qtz
Notes. (1) Mcar = Mn-Fe carbonate; Qtz = quartz; (2) for the sampling points, please see Figure 5.Table 2
EPMA data for sulfides from the Zhaxikang deposit.
Sample number
As (%)
S (%)
Se (%)
Cu (%)
Co (%)
Zn (%)
Sb (%)
Ag (%)
Au (%)
Cd (%)
Bi (%)
Pb (%)
NiO (%)
FeO (%)
Total (%)
Mineral
9-3-1
–
33.165
–
0.133
0.009
59.095
–
–
–
0.171
–
–
–
6.985
99.558
Sp
9-3-2
–
33.423
–
0.015
0.012
58.224
–
0.002
–
0.151
–
–
–
8.060
99.887
Sp
9-3-3
–
33.262
–
–
0.015
57.114
–
–
–
0.159
–
–
–
9.405
99.955
Sp
9-3-6
0.784
52.884
–
0.005
0.044
–
–
–
0.016
–
–
0.134
–
46.585
100.452
Py
9-3-7
2.244
51.934
–
0.002
0.049
–
–
–
0.004
–
–
0.132
–
46.044
100.409
Py
9-8-1
0.002
33.740
–
0.009
0.010
57.705
–
0.004
–
0.160
–
0.018
0.008
8.827
100.483
Sp
9-8-2
–
33.649
–
0.011
0.005
59.584
–
0.005
–
0.150
–
–
–
7.112
100.516
Sp
9-8-3
–
33.500
–
0.006
0.003
56.941
–
–
–
0.160
–
–
0.002
9.424
100.036
Sp
9-8-5
0.009
34.064
–
0.089
0.004
57.475
–
–
–
0.151
–
–
–
8.901
100.693
Sp
9-8-6
0.006
53.867
–
0.033
0.050
0.058
–
0.005
0.012
–
–
0.300
0.002
46.539
100.872
Py
9-8-7
0.015
53.359
–
0.471
0.043
0.044
–
–
–
–
–
0.174
–
46.518
100.624
Py
ZXK-1-1
–
33.502
–
0.012
0.009
57.497
–
–
–
0.172
–
–
–
9.142
100.334
Sp
ZXK-1-2
–
33.428
–
0.036
–
60.190
0.002
0.002
–
0.201
–
–
–
5.810
99.669
Sp
ZXK-1-3(1)
39.757
24.131
–
0.006
0.043
0.009
–
–
0.027
–
–
0.063
0.002
36.340
100.378
Apy
ZXK-1-3(2)
2.435
51.666
–
0.018
0.067
0.185
–
0.007
0.013
–
–
0.113
0.004
45.877
100.385
Py
ZXK-1-3(3)
–
33.337
0.002
0.014
0.012
58.623
–
0.015
–
0.147
–
–
0.003
7.747
99.900
Sp
ZXK-1-4
1.942
51.982
0.001
0.008
0.046
0.015
–
–
0.003
–
–
0.196
–
46.004
100.197
Py
ZXK-1-5
2.056
52.248
–
0.008
0.054
–
–
–
0.013
–
–
0.137
–
46.163
100.679
Py
ZXK-2-1
–
33.684
–
0.061
0.010
59.726
–
–
–
0.170
–
–
0.005
6.904
100.560
Sp
ZXK-2-2
39.643
23.789
–
0.002
0.050
0.008
0.196
0.006
–
–
–
0.081
0.098
36.245
100.118
Apy
ZXK-2-3
–
33.868
–
0.025
0.012
58.034
–
0.001
–
0.196
–
–
0.002
8.714
100.852
Sp
ZXK-2-4
39.444
24.229
–
–
0.041
–
0.173
–
0.026
–
–
0.061
–
36.424
100.398
Apy
ZXK-2-5
2.417
51.755
–
–
0.045
0.027
–
–
–
–
–
0.131
–
46.227
100.602
Py
ZXK-2-6
0.031
53.630
–
0.107
0.047
–
–
–
0.009
–
–
0.168
–
46.654
100.646
Py
Notes. (1) Py = pyrite, Sp = sphalerite, and Apy = arsenopyrite; (2) for the sampling points, please see Figure 5.Figure 7
(a) The negative correlation between Mn and Fe contents within Mn-Fe carbonate. (b) The negative correlation between Zn and Fe contents within sphalerite. (c) The negative correlation betweenδ
66Zn values and Fe contents within sphalerite from 9-3, 9-8, and ZXK-1. (d) The positive correlation between δ
66Zn values and Zn contents within sphalerite from 9-3, 9-8, and ZXK-1.
(a)
(b)
(c)
(d)


## 4.2. Fe-Zn Isotopes

All the Fe-Zn isotopic data are given in Table3. The annular polished section samples have δ
56
F
e
I
R
M
M
-
014 of -
1.95
‰
​
​
~
0.43
‰, with an average of −0.50‰  ± 1.09‰ (2SD, n = 19), and δ
66
Z
n
A
A
-
E
T
H of -
0.38
‰
​
~
0.07
‰ with an average of −0.25‰  ± 0.19‰ (2SD, n = 12). The Mn-Fe carbonate and pyrite show the δ
56Fe values range from −1.95‰ to −0.59‰ (average value of −0.97‰  ± 0.86‰; 2SD, n = 8) and from −0.26‰ to 0.23‰ (average value of −0.07‰  ± 0.35‰; 2SD, n = 9), respectively. The two arsenopyrite samples exhibit the δ
56Fe values of −1.01‰ (ZXK-1-3; Figure 5(c)) and −0.18‰ (ZXK-2-2; Figure 5(d)).Table 3
Fe-Zn isotopic data for annular polished section samples from the Zhaxikang deposit.
Sample number
Mineral
δ
56
F
e
I
R
M
M
–
014
2σ
δ
66
Z
n
A
A
–
E
T
H
2σ
9-3-1
Sp
−0.12
0.05
9-3-2
Sp
−0.23
0.05
9-3-3
Sp
−0.31
0.05
9-3-4
Mcar
−0.59
0.08
9-3-5
Mcar
−1.95
0.08
9-3-7
Py
−0.26
0.08
9-8-1
Sp
−0.17
0.05
9-8-2
Sp
−0.09
0.05
9-8-3
Sp
−0.35
0.05
9-8-4
Py
0.23
0.08
9-8-5
Sp
−0.23
0.05
9-8-6
Py
−0.26
0.08
9-8-7
Py
0.09
0.08
9-8-8
Mcar
−0.78
0.08
9-8-9
Mcar
−0.92
0.08
ZXK-1-1
Sp
−0.38
0.05
ZXK-1-2
Sp
−0.15
0.05
ZXK-1-3(1)
Apy
−1.01
0.08
ZXK-1-3(2)
Py
−0.22
0.08
ZXK-1-3(3)
Sp
−0.18
0.05
ZXK-1-4
Py
0.12
0.08
ZXK-1-5
Py
−0.05
0.08
ZXK-1-6
Mcar
−0.66
0.09
ZXK-1-7
Mcar
−0.78
0.08
ZXK-2-1
Sp
−0.32
0.05
ZXK-2-2
Apy
−0.18
0.08
ZXK-2-3
Sp
−0.28
0.05
ZXK-2-5
Py
−0.07
0.08
ZXK-2-6
Py
−0.17
0.08
ZXK-2-7
Mcar
−1.06
0.08
ZXK-2-8
Mcar
−1.04
0.08
Notes. (1) Abbreviations and sampling points are as in Tables 1 and 2; (2) 2σ is two times the standard deviation.In sample 9-3, the 9-3-1, 9-3-2 (Figure6(a)), and 9-3-3 sphalerite have the δ
66Zn values of −0.12‰, −0.23‰, and −0.31‰ with a gradually decreasing trend, the 9-3-4 and 9-3-5 Mn-Fe carbonate show the same trend with the δ
56Fe values of −0.59‰ and −1.95‰, and the δ
56Fe value of 9-3-7 pyrite is −0.26‰ (Figure 5(b)). Similarly, in sample 9-8, from 9-8-2 (−0.09‰; Figure 6(c)) through 9-8-5 (−0.23‰) to 9-8-3 (−0.35‰; Figure 6(d)), the δ
66Zn values of sphalerite also present a gradually decreasing trend; meanwhile the δ
56Fe values also decrease from −1.06‰ (9-8-8; Figure 6(b)) to −1.38‰ (9-8-9) for Mn-Fe carbonate and from 0.23‰ (9-8-4) to 0.09‰ (9-8-7) for pyrite, although the 9-8-1 sphalerite and 9-8-6 pyrite have the δ
66Zn value of −0.17‰ and δ
56Fe values of −0.26‰ (Figure 5(a)). In sample ZXK-1, the δ
66Zn values of ZXK-1-3 (−0.18‰; Figures 6(f) and 6(m)) and ZXK-1-2 (−0.15‰; Figure 6(g)) sphalerite are almost the same and heavier than ZXK-1-1 (−0.38‰) sphalerite; the δ
56Fe values for pyrite increase from ZXK-1-3 (−0.22‰) to ZXK-1-4 (0.12‰; Figure 6(e)) and ZXK-1-5 (−0.05‰); and, as for Mn-Fe carbonate, the ZXK-1-6 (−0.66‰) has heavier δ
56Fe values than ZXK-1-7 (−0.78‰; Figure 5(c)). By contrast, in sample ZXK-2, the ZXK-2-1 (−0.32‰; Figure 6(h)) and ZXK-2-3 (−0.28‰; Figures 6(k) and 6(p)) sphalerite have almost the same δ
66Zn values; the ZXK-2-7 (−1.06‰) and ZXK-2-8 (−1.04‰) Mn-Fe carbonate also yield approximate δ
56Fe values; also the ZXK-2-5 (Figures 6(i) and 6(n)) and ZXK-2-6 pyrite show the similar δ
56Fe values of −0.07‰ and −0.17‰ (Figure 5(d)).


## 5. Discussion

### 5.1. The Fe-Zn Isotopic and Elemental Variations

Sample 9-3 with typical concentric annular texture has the gradually decreasingδ
66Zn values of sphalerite and δ
56Fe values of Mn-Fe carbonate from core to edge (early to late; Figure 5(b)). Similarly, sample 9-8 shows the same isotopic variation trend of pyrite, sphalerite, and Mn-Fe carbonate except the sulfides in the core (Figure 5(a)). As the core in 9-8 consists of crushed, earlier formed slate breccias and sulfides, the annular sulfides and Mn-Fe carbonate formed around the core, and thus this core is not included in the decreasing trend. In sample ZXK-1, as the laminae (ZXK-1-3; Figure 5(c)) is regarded as the earliest, the Fe-Zn isotopic values of sphalerite and Mn-Fe carbonate still present gradually decreasing trend from ZXK-1-3 (−0.18‰) to ZXK-1-1 (−0.38‰) and ZXK-1-6 (−0.66‰) to ZXK-1-7 (−0.78‰), respectively (Figure 5(c)), whereas the δ
56Fe values of pyrite show the different variation, which might be related to the influence of Fe isotopic fractionation within pyrite-arsenopyrite-sphalerite mineral pair (Figures 6(f) and 6(m)). However, there is no similar decreasing Fe-Zn isotopic variation trend in sample ZXK-2, and the modification by the stage 3 quartz vein must be the main cause (Figures 5(d), 6(h)–6(l) and 6(n)–6(p)).Previous studies [12, 62, 99] have proposed a Rayleigh distillation model to explain an increasing trend in δ
66Zn values within precipitates over time for the hydrothermal fluid. This Rayleigh distillation model is as follows: the ore-forming materials derived from a single source would be subjected to kinetic Rayleigh fractionation that would lead to the early formed mineral precipitants being preferentially enriched in light isotopes, as well as residual fluids and later precipitants with heavier isotopic values, causing an increasing trend in isotopic values within precipitants over time. Several previous studies have used this distillation model to explain the Zn isotopic variation within different types of deposits (e.g., VHMS: [18]; Irish-type: [19, 100]; SEDEX: [20, 53]). Likewise, this Rayleigh distillation model is also applicable to the Fe isotopic variation in skarn deposits [21, 22]. However, the δ
66Zn and δ
56Fe values gradually decrease from early to late stages within Zhaxikang deposit, which cannot be explained by this distillation model. Another Rayleigh distillation mechanism models this decreasing trend: the metallogenic elements are transported by the ore-forming system consisting of vapour and liquid phases, and there is partitioning between vapour-liquid phases and the ratios change with the temperature decreasing. Then the minerals precipitate from the liquid phase of the ore-forming system. During this period, the vapour-liquid partitioning and mineral precipitation cause the Rayleigh fractionation, and this Rayleigh fractionation leads to the mineral precipitation being preferentially enriched in heavy isotopes relative to the ore-forming system. Thus, the isotopic values of subsequent minerals are lighter and lighter [101, 102].This Rayleigh distillation model is supported by the following evidence: (1) The vapour-liquid partitioning and related isotopic fractionation for transition metal elements (e.g., Cu and Mo) have been confirmed by previous research in the Dahutang W-Cu-Mo ore field [103]; (2) minerals typically precipitate from the liquid phase; however, according to the previous literature [58], in unique cases the vapour phase containing metal can even directly condensate to form solid phases from high-temperature ore-forming system (e.g., VMS and volcano related system), which can demonstrate the existence of the vapour phase for metal elements. As for the Zhaxikang deposit, the theoretical calculated ore-forming temperature from Fe isotopic data is 500~800°C [26], and thus there should be a transitory high-temperature period, making the vapour-liquid partitioning possible; (3) the fluid inclusion data demonstrate that Mn-Fe carbonates and sulfides exist in three types of inclusions: ① gas-liquid two-phase water inclusions (W type, more than 90%), ② pure liquid inclusions (L type), and ③ pure CO2 inclusions (PG type) [104] from Zhaxikang deposit.On the other hand, for sphalerite, there are positive correlations between Zn contents andδ
66Zn values, with negative correlations between Fe contents and δ
66Zn values in samples 9-8 (Figure 5(a)), 9-3 (Figure 5(b)), and ZXK-1 (Figure 5(c)), respectively. Moreover, plotting all the data from these three samples on a diagram together, the correlations are also good with R
2 > 0.6 (Figures 7(c) and 7(d)). Zinc and iron usually show similar geochemical behaviour as both of them are highly mobile in chloride-bearing hydrothermal fluids [23, 100]. Thus, the Zn2+ and Fe2+ ions are preferentially enriched in the liquid phase relative to the vapour phase before precipitation [23, 100], which cause the decreasing Zn contents of sphalerite over time. As the total content of Zn2+ and Fe2+ is constant in sphalerite, the Fe contents of sphalerite gradually increase with the decreasing Zn contents. These correlations further support the hypothesis that the ore-forming system is the mixture of vapour and liquid phases. In addition, the sample ZXK-2 cut by a later stage 3 quartz vein does not present the same correlations and variations, which is a new evidence for two pulses of mineralization proposed by Zheng et al. [7] and Wang et al. [26].All of the evidence above reveals that the ore-forming elements are transported by the ore-forming system that consists of vapour and liquid phases. The vapour-liquid partitioning and mineral precipitation are the main cause of Fe-Zn isotopic and elemental variations. Afterwards, the overprint by the second pulse of mineralization has also partly modified the Fe-Zn isotopic and elemental compositions of some earlier samples (Figure5(d)) [26].


### 5.2. The Fe-Zn Isotopic Fractionation Models for the Ore-Forming System

In order to verify the Rayleigh distillation model in Section5.1 and obtain more information of the ore-forming system, we use the following equations to establish Fe-Zn isotopic fractionation models for the ore-forming system and mineral precipitation:(1)
δ
M
i
n
e
r
a
l
s
‰
=
α
α
i
×
δ
i
+
1000
×
F
α
m
-
1
-
1000
,
δ
O
r
e
-
f
o
r
m
i
n
g
S
y
s
t
e
m
‰
=
1
α
i
×
δ
i
+
1000
×
F
α
m
-
1
-
1000
.
δ
M
i
n
e
r
a
l
s, δ
O
r
e
-
f
o
r
m
i
n
g
S
y
s
t
e
m, and δ
i are the δ
56Fe-δ
66Zn values of momentary mineral precipitation and momentary and initial ore-forming system, respectively; α, α
i, and α
m are the isotopic fractionation factors between ore-forming system and mineral precipitation that refer to the momentary condensation temperature T, the initial temperature of ore-forming system T
i, and (
T
+
T
i
)
/
2, respectively; and F is the fraction of remaining ore-forming system that consists of vapour and liquid phases [102, 105]. In addition, the equations for fractionation factors are approximately using ln α
F
e
=
0.4432
×
10
6
/
T
2 for Fe [106] and ln α
Z
n
=
0.2853
×
10
6
/
T
2
+
0.0535 for Zn (T is absolute temperature in K) [107].Wang et al. [26] have calculated the ore-forming temperature (500~800°C) of the first pulse of mineralization in Zhaxikang deposit using the Fe isotopic fractionation between pyrite and Mn-Fe carbonate. Although this temperature range is a little wide, the highest temperature of such ore-forming system can reach around 500°C according to the previous ore deposit studies [6, 22, 23]; hence we regard 500°C as the initial temperature of the ore-forming system. The homogenization temperature (240°C) of the fluid inclusions [104] from the first pulse of mineralization in Zhaxikang deposit is considered as the momentary condensation temperature. Furthermore, for the purpose of making the fractionation models more comprehensive and exact, we also quote the Fe-Zn isotopic data of pyrite (δ
56Fe: stage 1: −0.33‰ to −0.09‰; stage 2: −0.30‰ to 0.19‰; stage 3: 0.16‰ to 0.43‰), sphalerite (δ
66Zn: −0.31‰ to 0.07‰), and Mn-Fe carbonate (δ
56Fe: −0.80‰ to −0.55‰; δ
66Zn: −0.11‰ to 0.04‰) from Wang et al. [26], as well as the δ
66Zn values of sphalerite (−0.25‰ to 0.03‰) and Mn-Fe carbonate (−0.01‰) from Duan et al. [5]. Finally, we set up 12 Fe-Zn isotopic fractionation models for pyrite, sphalerite, and Mn-Fe carbonate (Figure 8) with different δ
56Fei values of 0‰ (mean value of magma) [51], −0.5‰, and −1‰, as well as δ
66
Z
n
i values of −0.28‰ (the lightest value of seafloor hydrothermal fluids) [62], 0‰ (mean value of bulk earth) [54], and 0.23‰ (the mean value of deep sea water) [63, 108, 109].Figure 8
The Fe-Zn isotopic fractionation models for the ore-forming system and minerals (pyrite, sphalerite, and Mn-Fe carbonate) with differentδ
56Fei and δ
66Zni values.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)These fractionation models show that theF ranges for ore-forming system highly depend on the δ
i values (Figure 8). The details are as follows: the pyrite covers the F ranges of 30.1
%
​
~
62.5
% (δ
56Fei = 0‰; Figure 8(a)), 48.8
%
​
~
99.0
% (δ
56Fei = −0.5‰; Figure 8(b)), and more than 77.7% (δ
56Fei = −1‰; Figure 8(c)). The F ranges for sphalerite are 38.7
%
​
~
70.9
% (δ
66Zni = −0.28‰; Figure 8(d)), 26.7
%
​
~
48.7
% (δ
66Zni = 0‰; Figure 8(e)), and 19.8
%
​
~
36.3
% (δ
66Zni = 0.23‰; Figure 8(f)), respectively. In comparison, the Mn-Fe carbonates has the F ranges of 7.1
%
​
~
24.9
% (δ
56Fei = 0‰; Figure 8(g)), 10.9
%
​
~
39.7
% (δ
56Fei = −0.5‰; Figure 8(h)), and 17.3
%
​
~
63.3
% (δ
56Fei = −1‰; Figure 8(i)) for Fe isotope, as well as 55.6
%
​
~
68.1
% (δ
66Zni = −0.28‰; Figure 8(j)), 38.2
%
​
~
46.8
% (δ
66Zni = 0‰; Figure 8(k)), and 28.5
%
​
~
34.8
% (δ
66Zni = 0.23‰; Figure 8(l)) for Zn isotope. All of these results suggest that the Fe-Zn isotopic data of Zhaxikang deposit fit these Rayleigh fractionation models well.However, we need to consider the following facts in the Zhaxikang deposit: (1) the second pulse of mineralization has brought some Fe to form the stage 3 pyrite with heavierδ
56Fe values (0.23
‰
​
~
0.43
‰) [26]; thus the stage 3 pyrite does not fit the fractionation models; (2) as most of the sphalerite and pyrite are paragenetic during the first pulse of mineralization (Figures 4(g), 4(i), 4(k), 5(a)–5(c) and 6(c)), especially in the earliest lamina (Figures 4(a)–4(c), 5(c), 6(f) and 6(m)), the sphalerite and pyrite should overlap on F ranges; (3) in theory, the Mn-Fe carbonates would have the same F range for Fe-Zn isotopes. Nonetheless, in view of the tight Zn isotopic variation range, the F range for Fe isotope should cover that for Zn isotope; (4) during the earlier period, as the ore-forming system consisting of vapour and liquid phases is dominant, F values in fractionation models should be large. Taking all these facts and 12 fractionation models into consideration, the δ
56Fei value is supposed to be in the range of -
0.5
‰
​
~−1‰, and the δ
66
Z
n
i value should be between −0.28‰ and 0‰.


### 5.3. Implications for the Genesis of Zhaxikang Deposit

#### 5.3.1. Excluding the Possibility of Hot Spring Genetic Model

The hot spring model predicts that metals (e.g., Zn, Pb, Sb, Ag, and Fe) are leached from sedimentary wall rocks, which is supported by the following evidence: (1) theδ
34S values of the sulfides (4.5
‰
​
~
12
‰) are similar to those of sedimentary wall rocks (4.93
‰
​
~
11.49
‰); (2) the δ
30Si values of quartz (-
0.90
‰
​
~−0.40‰) are the same as those of siliceous rocks with hot spring genesis; (3) the δ
D
V
–
S
M
O
W (-
135
‰
​
~−127‰) and δ
18
O
H
2
O values (-
13.7
‰
​
~
12.4
‰) of fluid inclusions trapped in quartz are similar to those of the south Tibetan hot spring; (4) the Pb isotopes (P
b
206
/
P
b
204: 18.474~19.637; P
b
207
/
P
b
204: 15.649~15.774; P
b
208
/
P
b
204: 39.660~40.010) show the characteristics of radiogenic Pb; (5) the He-Ar isotopes demonstrate the contribution of crustal fluid and meteoric water [3, 4].However, this genetic model is inconsistent with textural and Fe-Zn isotopic evidence presented here. Firstly, the primary sedimentary wall rocks in the orefield are slate. As suggested by the continuous batch experimental research of Fernandez and Borrok [27], the ore-forming fluid would preferentially leach out the heavy Zn isotopes. Likewise, Chen et al. [110] have analyzed the Zn isotopic compositions of samples from 8 hot springs, and the results show most of the hot springs have relatively constant and heavier δ
66Zn values (approximately 0.42‰) than host rocks (−0.42‰ to 0.14‰). Therefore, if the metallogenic elements are leached from the slate by hot spring, there would be some sphalerite with heavier δ
66Zn values than these slates in Zhaxikang deposit. Nevertheless, the δ
66Zn values of the slate from Zhaxikang orefield range from −0.23‰ to 0.10‰ that are similar to sphalerite, especially the unmodified slate sample with the δ
66Zn value of 0.10‰ that is even a little heavier than that of sphalerite (−0.38‰ to 0.07‰; Figure 9(b)) [26]. As the slate has heavier Zn isotopic compositions than those of host rocks from Chen et al. [110], the δ
66Zn values of the hot spring in Zhaxikang orefield would even be heavier than 0.42‰, which are much heavier than the δ
66Zni values (-
0.28
‰
​
~
0
‰) of the fractionation models in Section 5.2. Secondly, Sharam et al. [46] have measured the δ
56Fe values (−0.59‰ to −0.12‰) of hot springs in Juan de Fuca Ridge, which is much heavier than the δ
56Fei value (-
0.5
‰
​
~−1‰) gained from the fractionation models. Moreover, the marine fluids usually have lighter Fe isotopic compositions (Figure 9(a)); thus the hot springs in Tibet would have heavier δ
56Fe values than those of Sharam et al. [46].Figure 9
(a) Fe isotopic compositions of pyrite, Mn-Fe carbonate, and arsenopyrite in the Zhaxikang deposit (part of the data is quoted from Wang et al. [26]). Other Fe isotopic data for the Bulk Silicate Earth [10], shale [34, 35], igneous rocks [21–23, 28, 36–45], seafloor hydrothermal fluid system [34, 46–49], midoceanic ridges pyrite [48, 50], the Bayan Obo Fe-REE magmatic-type deposit in China [51], the skarn-type deposits in Tongling ore district in China [21, 22], the Schwarzwald hydrothermal vein deposit in Germany [52], the Renison Sn-W deposit in Australia [23], and the Dongshengmiao SEDEX-type deposit in China [53] are also plotted for comparison. (b) Zn isotopic compositions of sphalerite, Mn-Fe carbonate, and slate in the Zhaxikang deposit (part of the data is quoted from Wang et al. [26]), compared with the Bulk Silicate Earth [54], sedimentary rocks [55–57], igneous rocks [28, 54, 57–60], deep-sea carbonates [61], seafloor hydrothermal fluid system [62], deep sea water [51, 54, 63], and other deposits with different geneses: the Gorno and Raibl magmatic-type deposit in Italy [14], the skarn-type deposits in the Tongling ore district in China [12], the Tianqiao and Bangbangqiao carbonated-hosted Pb-Zn sulfide deposits in China [64], the Irish-type deposit in Ireland [19], the Cévennes MVT deposit in France [65], the Alexandrinka VHMS-type deposit in Russia [18], the Red Dog SEDEX-type ore district in Alaska [20], and the Dongshengmiao SEDEX-type deposit in China [53].
(a)
(b)The possibility of hot spring genesis for the first pulse of mineralization event can be excluded by Fe-Zn isotopic data. Evidence for the second pulse of mineralization is based on the fact that the Fe-Zn isotope values do not follow similar concentric patterns with the ore textures as seen in sample ZXK-2. Meanwhile, the evidence from Si-H-O isotopes demonstrates that the second pulse of mineralization may be related to hot spring. Additionally, the Fe-Zn isotopic data demonstrate that the sedimentary wall rocks have not provided significant amounts of metals, although the S-Pb isotopic data show that these wall rocks constitute some contribution for S-Pb [4, 5], whereas the similar Zn isotopic compositions of slate and sphalerite suggest that they share the same Zn origin.


#### 5.3.2. Inconsistency with the Magmatic-Hydrothermal Fluid Genetic Models

There are two genetic models for the magmatic-hydrothermal fluid genesis. In the first model, Duan et al. [5] considered that the genesis of Zhaxikang deposit relates to the mid-low temperature magma-related hydrothermal activity and that the metallogenic elements are mainly sourced from the mixing of basement and the sedimentary wall rocks. The evidence is mainly from the Zn-S-Pb isotopes that we mentioned in Section 1.Duan et al. [5] have analyzed the δ
66Zn values of the sulfides (-
0.25
‰
​
~
0.03
‰) and basement rocks (0.05
‰
​
~
0.21
‰). The dominating sedimentary wall rocks in the orefield are slate (δ
66Zn values: -
0.23
‰
​
~
0.10
‰) [26], which is formed by the epimetamorphism of shale and sandstone. Meanwhile, combining the data from Wang et al. [26] with this study, the sphalerite should have a range from −0.38‰ to 0.07‰ in Zhaxikang deposit. Both the basement and sedimentary wall rocks have heavier Zn isotopic compositions than sphalerite in Zhaxikang deposit. However, just like we discussed in Section 5.3.1, if the Zn is sourced from mixing of the basement and sedimentary wall rocks, there should be some sphalerite with heavier δ
66Zn values than these rocks. This inference is also evidenced by the research of Zhou et al. [64]: the Paleozoic carbonate host rocks and Precambrian basements are considered to be the origin of metals, and these rocks have lighter δ
66Zn values (−0.52‰ to 0.16‰) than the sphalerite from the Tianqiao (−0.54‰ to 0.30‰) and Bangbangqiao (−0.21‰ to 0.43‰) deposits in the Sichuan-Yunnan-Guizhou Pb-Zn metallogenic province (Figure 9(b)). Additionally, in respect of the Fe isotope, the Schwarzwald hydrothermal vein deposit in Germany can be used as an analogy [52]. Iron in this deposit originates from the basement consisting of granites and gneisses, as well as sedimentary rocks including shale and sandstone. The basement of Zhaxikang is composed of dolerite, quartz diorite, rhyolite porphyry, pyroclastics, and porphyritic monzogranite [5], and these crust-derived igneous rocks have the similar Fe isotopic composition with granites according to the data from previous studies [21–23, 28, 36–45]. Likewise, the δ
56Fe values of shale and sandstone in the Schwarzwald deposit are −0.21‰, 0.03‰, and 0.22‰, all of which fall into the Fe isotopic variation range of shale (-
0.39
‰
​
~
0.71
‰) from Beard et al. [34] and Rouxel et al. [35]. Therefore, we regard the fact that the metal-sourced rocks in these two deposits have the similar Fe isotopic compositions. And, yet, Markl et al. [52] suggested that fluid-rock interaction make the ore-forming fluid have the δ
56Fe value of -
0.5
‰
​
~
0
‰, which is heavier than the δ
56Fei value (-
0.5
‰
​
~−1‰) of the fractionation models in Section 5.2. Furthermore, it is generally known that the equilibrium isotope fractionation is a function of temperature, with larger fractionation generated at lower temperatures [111]. Although the temperature of ore-forming fluid (100~200°C) [52] in the Schwarzwald deposit is lower than that in Zhaxikang deposit (~250°C) [5], the lightest δ
56Fe value (−1.36‰) of siderite from this deposit is much heavier than that of Mn-Fe carbonate (−1.95‰) from Zhaxikang deposit (Figure 9(a)). All of these inconsistencies from Fe-Zn isotopic data indicate that this genetic model is not appropriate for Zhaxikang deposit.In the second genetic model, Xie et al. [6] suggested the mineralization in Zhaxikang deposit was genetically related to the Miocene dome-related magmatism. This magmatism generated the pegmatite and two-mica granite in the core of the regional domes (Figure 1(b); e.g., Cuonadong, Yalaxiangbo, Ranba, and Kangma), and the high-temperature ore-forming fluids were derived from magmatic melts exsolution. This hypothesis is based on the research of field geology, petrography, melt and fluid inclusions, and C−H−O isotopes: (1) the δ
13
C
V
–
P
D
B (-
6.1
‰
​
~−6.9‰) and δ
18
O
V
–
S
M
O
W (+
9.9
‰
​
~+11.8‰) values of rhodochrosite show the mantle origin; (2) the δ
13
D
V
–
S
M
O
W ( -
144.8
‰
​
~−110‰) and δ
18
O
V
–
S
M
O
W (-
9.85
‰
​
~−8.89‰) values, the low salinity (0.2~7.9 wt.% NaCleqv), high temperature (298~457°C), and rich CO2 content with minor CH4, N2, C2H6, C3H8, and C6H6 of the melt and fluid inclusions trapped in quartz and beryl within pegmatite indicate the magmatic origin; (3) the 40Ar-39Ar plateau age of pegmatite (18.93 ± 0.27 Ma) is similar to stage 5 quartz-pyrite-stibnite (17.9 ± 0.5 Ma).The pegmatite and two-mica granite in the domes have the typical characteristics of S-type granite: (1) containing abundant Al-rich minerals, A/CNK: 1.07~1.24; (2) w(SiO2): 73.26
%
​
~
74.33
%, w(K2O)/w(Na2O): 1.1~1.2, w(FeO)/w(Fe + Mn): 0.64~0.76; (3) the content of corundum molecules > 1% in the CIPW standard minerals [112, 113]. Meanwhile, based on the facts that all the Fe-bearing minerals in Zhaxikang deposit are ferrous, and the magmatic-hydrothermal fluid is CO2-rich with minor CH4, N2, C2H6, C3H8, and C6H6 [6], and the parent magma is most likely S-type reduced magma. In consideration of the situation in the Renison Sn-W deposit that we mentioned in Section 1, the ore-forming fluid is considered to exsolve from S-type reduced magma and the sulfides have heavier Fe isotopic compositions than ore-related igneous rocks [23]. The case in Zhaxikang deposit is contrary yet the pyrite from the first pulse of mineralization has the δ
56Fe value of -
0.33
‰
​
~
0.23
‰ that is much lighter than those of granitoids (−0.08% to 0.59%) [22, 28]. Moreover, Heimann et al. [43] proposed that we would not expect Fe isotopic compositions of high F/Cl magmatic/hydrothermal systems to significantly deviate from the average of igneous rocks, hence as Xie et al. [6] considered that the Zhaxikang parent magma has high F contents, the magmatic-hydrothermal fluid in this genetic model would have the similar Fe isotopic compositions with granitoids, which is not in line with our fractionation models. Even if the parent magma in Zhaxikang orefield is oxidized-type and is similar to the magma related to I-type granitoids with δ
56Fe values of -
0.04
%
​
~
0.59
% in Tongling ore district (Figure 9(a)) [21, 22], it is still hard to produce so light δ
56Fei value (-
0.5
‰
​
~−1‰). Similarity, Zn isotopic data do not support this genetic model, either. Chen et al. [54] studied the Zn isotopic fractionation during igneous process and suggested that the maximum Zn isotopic variation induced by high-temperature igneous processes is no larger than δ
66Zn~0.10%. Besides this, Telus et al. [28] measured the δ
66Zn values of pegmatite (0.25
‰
​
~
0.59
‰) and some other granitoids (-
0.16
‰
​
~
0.21
‰) and then found there is even no variation in δ
66Zn values during fluid exsolution in some cases. As the granitoids in the regional domes are principally dominated by pegmatite and the magmatic fluids which have high temperature of 298~457°C [6], it is also hard to generate δ
66Zni value between −0.28‰ and 0‰ as yielded from the Fe-Zn isotopic fractionation models. Consequently, the Fe-Zn isotopic data are also not in favor of the second genetic model.On the other hand, neither of these two genetic models can explain the different Fe-Zn isotopic and elemental variations in sample ZXK-2, which is considered to result from the overprint by the second pulse of mineralization. Overall, all of the Fe-Zn isotopic and elemental evidences are inconsistent with both of these magmatic-hydrothermal fluid genetic models. Nevertheless, the evidence for these two genetic models may prove that the second pulse of mineralization is related to magmatic-hydrothermal fluids.


#### 5.3.3. Constraints on SEDEX Modified by Hydrothermal Fluid Genetic Model

Zheng et al. [7, 25] considered that the first pulse of mineralization (Pb-Zn) has the SEDEX genesis, and the second pulse of mineralization (Sb-Ag) is related to hot spring that overprints the earlier mineralization. The previous evidences are mainly as follows: (1) the evidence from ore textures and obviously late Sb mineralization compared to Pb-Zn mineralization in Section 2.3; (2) the stage 1 and 2 ores have high Mn, Fe, Ba, and B contents, Ga ≫ In (Ga/In: 1.49~4.47), Pb + Zn ≫ Cu, and host rocks are exhalative rocks and Mn-Fe carbonates; (3) the δ
34S values of sphalerite and galena (7.7
‰
​
~
12
‰) are different from those of stibnite (4.5
‰
200
B
~
7.1
‰); (4) the data from valentinite (δ
30Si values: -
0.9
‰
​
~−0.4‰) and fluid inclusions trapped in quartz (δ
D
V
–
S
M
O
W: -
162
‰
​
~−142‰, δ
18
O
H
2
O: -
12.9
‰
​
~−1.9‰, homogenization and freezing temperatures: 164~313°C and −3.2​
~−0.7°C, salinity: 0.7
%
​
~
5.3
%, and density: 0.74~0.93 g/cm3) both demonstrate the Sb mineralization is related to the south Tibetan hot spring; (5) the Rb-Sr isotopic isochrone age of stage 2 sphalerite is 147.2 ± 3.2 Ma that is similar to the Jurassic sedimentary wall rocks, whereas a later stage quartz-pyrite-stibnite vein has the A
r
40-A
r
39 plateau age of 17.9 ± 0.5 Ma.Wen et al. [17] investigated several Pb-Zn deposits with different geneses in China, and the results show that the Zn/Cd ratios of sphalerite vary with different geneses: (1) high-temperature systems including the porphyry, magmatic hydrothermal, skarn, and volcanic hosted massive sulfide (VMS)-type deposits: 155~223; (2) low-temperature systems that include the Mississippi Valley-type (MVT) deposits: 17~201; (3) SEDEX-type deposits of exhalative systems: 316~368; (4) seafloor hydrothermal sulfides of exhalative systems: 211~510. According to the EPMA data, the Zn/Cd ratios of sphalerite range from 296 to 399 in Zhaxikang deposit, which overlap the range of exhalative systems and much higher than those of high-temperature and low-temperature systems.The Fe-Zn isotopic data also conform to the SEDEX modified by hydrothermal fluid genetic model. Firstly, as discussed in Sections5.3.1 and 5.3.2, neither the hot spring nor the magmatic-hydrothermal fluids can have the δ
56Fei (-
0.5
‰
​
~
-
1
‰) and δ
66Zni (-
0.28
‰
​
~
0
‰) values to satisfy the Fe-Zn isotopic fractionation models. However, the Fe-Zn isotopic compositions of seafloor hydrothermal fluid system covers the range of -
1.79
‰
​
~
0.04
‰ for Fe isotope and -
0.28
‰
​
~
0.96
‰ for Zn isotope according to previous studies (Figure 9) [38, 46–49, 62], which can generate the ore-forming system with δ
56Fei and δ
66Zni values to meet the fractionation models. Secondly, although there are large overlaps in Zn isotopic compositions among deposits with different geneses, the Zn isotopic compositions of Zhaxikang deposit is most similar to the Alexandrinka VHMS-type and Red Dog SEDEX-type deposits with marine origin, as well as obviously distinguishing from the narrow range of magmatic-hydrothermal deposits (Figure 9(b)). Thirdly, the lightest δ
56Fe value in Zhaxikang deposit is −1.95‰, and only the minerals with marine origin (mid-oceanic ridges pyrite) or SEDEX genesis (sphalerite and pyrrhotite in Dongshengmiao SEDEX-type deposit) can have so light δ
56Fe values (Figure 9(a)). Fourthly, under the conditions of high P
C
O
2 and low PH (<8), the Zn precipitated as sulfides is isotopically nearly unfractionated with respect to the parent hydrothermal fluid, whereas, under the conditions of high P
C
O
2 and high-PH (>9), negative δ
66Zn values down to 0.6‰ can be expected in sulfides precipitated from hydrothermal fluid [107]. In Zhaxikang deposit, P
C
O
2 would be high as there are plenty of Mn-Fe carbonates; meanwhile, owing to the facts that the modern seawater has the PH around 8 and there is more CO2 in Jurassic atmosphere than present, the Jurassic seawater would have a lower PH than modern seawater. This can well explain that the δ
66Zn values of sphalerite (-
0.38
‰
​
~
0.07
‰) slightly fractionate with the ore-forming system (δ
66Zni value: -
0.28
‰
​
~
0
‰). Besides these, the overprinting of earlier ores by second pulse of mineralization is not only proved by the different Fe-Zn isotopic and elemental variations in sample ZXK-2 from the other 3 samples in this research (Figure 5) but also evidenced by the temporally increasing δ
56Fe and decreasing δ
66Zn values recorded in this deposit that coincided with an increase in alteration [26]. Nevertheless, further research is required to confirm whether this hydrothermal fluid is hot spring or magmatic-hydrothermal fluid.To sum up, among the various genetic models, the Fe-Zn isotopic and EPMA evidence indicate the SEDEX modified by hydrothermal fluid genetic model is the most plausible. Our research also demonstrates that the Fe-Zn isotopes have the potential to trace the metal source and provide insights into ore-forming processes.


## 5.1. The Fe-Zn Isotopic and Elemental Variations

Sample 9-3 with typical concentric annular texture has the gradually decreasingδ
66Zn values of sphalerite and δ
56Fe values of Mn-Fe carbonate from core to edge (early to late; Figure 5(b)). Similarly, sample 9-8 shows the same isotopic variation trend of pyrite, sphalerite, and Mn-Fe carbonate except the sulfides in the core (Figure 5(a)). As the core in 9-8 consists of crushed, earlier formed slate breccias and sulfides, the annular sulfides and Mn-Fe carbonate formed around the core, and thus this core is not included in the decreasing trend. In sample ZXK-1, as the laminae (ZXK-1-3; Figure 5(c)) is regarded as the earliest, the Fe-Zn isotopic values of sphalerite and Mn-Fe carbonate still present gradually decreasing trend from ZXK-1-3 (−0.18‰) to ZXK-1-1 (−0.38‰) and ZXK-1-6 (−0.66‰) to ZXK-1-7 (−0.78‰), respectively (Figure 5(c)), whereas the δ
56Fe values of pyrite show the different variation, which might be related to the influence of Fe isotopic fractionation within pyrite-arsenopyrite-sphalerite mineral pair (Figures 6(f) and 6(m)). However, there is no similar decreasing Fe-Zn isotopic variation trend in sample ZXK-2, and the modification by the stage 3 quartz vein must be the main cause (Figures 5(d), 6(h)–6(l) and 6(n)–6(p)).Previous studies [12, 62, 99] have proposed a Rayleigh distillation model to explain an increasing trend in δ
66Zn values within precipitates over time for the hydrothermal fluid. This Rayleigh distillation model is as follows: the ore-forming materials derived from a single source would be subjected to kinetic Rayleigh fractionation that would lead to the early formed mineral precipitants being preferentially enriched in light isotopes, as well as residual fluids and later precipitants with heavier isotopic values, causing an increasing trend in isotopic values within precipitants over time. Several previous studies have used this distillation model to explain the Zn isotopic variation within different types of deposits (e.g., VHMS: [18]; Irish-type: [19, 100]; SEDEX: [20, 53]). Likewise, this Rayleigh distillation model is also applicable to the Fe isotopic variation in skarn deposits [21, 22]. However, the δ
66Zn and δ
56Fe values gradually decrease from early to late stages within Zhaxikang deposit, which cannot be explained by this distillation model. Another Rayleigh distillation mechanism models this decreasing trend: the metallogenic elements are transported by the ore-forming system consisting of vapour and liquid phases, and there is partitioning between vapour-liquid phases and the ratios change with the temperature decreasing. Then the minerals precipitate from the liquid phase of the ore-forming system. During this period, the vapour-liquid partitioning and mineral precipitation cause the Rayleigh fractionation, and this Rayleigh fractionation leads to the mineral precipitation being preferentially enriched in heavy isotopes relative to the ore-forming system. Thus, the isotopic values of subsequent minerals are lighter and lighter [101, 102].This Rayleigh distillation model is supported by the following evidence: (1) The vapour-liquid partitioning and related isotopic fractionation for transition metal elements (e.g., Cu and Mo) have been confirmed by previous research in the Dahutang W-Cu-Mo ore field [103]; (2) minerals typically precipitate from the liquid phase; however, according to the previous literature [58], in unique cases the vapour phase containing metal can even directly condensate to form solid phases from high-temperature ore-forming system (e.g., VMS and volcano related system), which can demonstrate the existence of the vapour phase for metal elements. As for the Zhaxikang deposit, the theoretical calculated ore-forming temperature from Fe isotopic data is 500~800°C [26], and thus there should be a transitory high-temperature period, making the vapour-liquid partitioning possible; (3) the fluid inclusion data demonstrate that Mn-Fe carbonates and sulfides exist in three types of inclusions: ① gas-liquid two-phase water inclusions (W type, more than 90%), ② pure liquid inclusions (L type), and ③ pure CO2 inclusions (PG type) [104] from Zhaxikang deposit.On the other hand, for sphalerite, there are positive correlations between Zn contents andδ
66Zn values, with negative correlations between Fe contents and δ
66Zn values in samples 9-8 (Figure 5(a)), 9-3 (Figure 5(b)), and ZXK-1 (Figure 5(c)), respectively. Moreover, plotting all the data from these three samples on a diagram together, the correlations are also good with R
2 > 0.6 (Figures 7(c) and 7(d)). Zinc and iron usually show similar geochemical behaviour as both of them are highly mobile in chloride-bearing hydrothermal fluids [23, 100]. Thus, the Zn2+ and Fe2+ ions are preferentially enriched in the liquid phase relative to the vapour phase before precipitation [23, 100], which cause the decreasing Zn contents of sphalerite over time. As the total content of Zn2+ and Fe2+ is constant in sphalerite, the Fe contents of sphalerite gradually increase with the decreasing Zn contents. These correlations further support the hypothesis that the ore-forming system is the mixture of vapour and liquid phases. In addition, the sample ZXK-2 cut by a later stage 3 quartz vein does not present the same correlations and variations, which is a new evidence for two pulses of mineralization proposed by Zheng et al. [7] and Wang et al. [26].All of the evidence above reveals that the ore-forming elements are transported by the ore-forming system that consists of vapour and liquid phases. The vapour-liquid partitioning and mineral precipitation are the main cause of Fe-Zn isotopic and elemental variations. Afterwards, the overprint by the second pulse of mineralization has also partly modified the Fe-Zn isotopic and elemental compositions of some earlier samples (Figure5(d)) [26].


## 5.2. The Fe-Zn Isotopic Fractionation Models for the Ore-Forming System

In order to verify the Rayleigh distillation model in Section5.1 and obtain more information of the ore-forming system, we use the following equations to establish Fe-Zn isotopic fractionation models for the ore-forming system and mineral precipitation:(1)
δ
M
i
n
e
r
a
l
s
‰
=
α
α
i
×
δ
i
+
1000
×
F
α
m
-
1
-
1000
,
δ
O
r
e
-
f
o
r
m
i
n
g
S
y
s
t
e
m
‰
=
1
α
i
×
δ
i
+
1000
×
F
α
m
-
1
-
1000
.
δ
M
i
n
e
r
a
l
s, δ
O
r
e
-
f
o
r
m
i
n
g
S
y
s
t
e
m, and δ
i are the δ
56Fe-δ
66Zn values of momentary mineral precipitation and momentary and initial ore-forming system, respectively; α, α
i, and α
m are the isotopic fractionation factors between ore-forming system and mineral precipitation that refer to the momentary condensation temperature T, the initial temperature of ore-forming system T
i, and (
T
+
T
i
)
/
2, respectively; and F is the fraction of remaining ore-forming system that consists of vapour and liquid phases [102, 105]. In addition, the equations for fractionation factors are approximately using ln α
F
e
=
0.4432
×
10
6
/
T
2 for Fe [106] and ln α
Z
n
=
0.2853
×
10
6
/
T
2
+
0.0535 for Zn (T is absolute temperature in K) [107].Wang et al. [26] have calculated the ore-forming temperature (500~800°C) of the first pulse of mineralization in Zhaxikang deposit using the Fe isotopic fractionation between pyrite and Mn-Fe carbonate. Although this temperature range is a little wide, the highest temperature of such ore-forming system can reach around 500°C according to the previous ore deposit studies [6, 22, 23]; hence we regard 500°C as the initial temperature of the ore-forming system. The homogenization temperature (240°C) of the fluid inclusions [104] from the first pulse of mineralization in Zhaxikang deposit is considered as the momentary condensation temperature. Furthermore, for the purpose of making the fractionation models more comprehensive and exact, we also quote the Fe-Zn isotopic data of pyrite (δ
56Fe: stage 1: −0.33‰ to −0.09‰; stage 2: −0.30‰ to 0.19‰; stage 3: 0.16‰ to 0.43‰), sphalerite (δ
66Zn: −0.31‰ to 0.07‰), and Mn-Fe carbonate (δ
56Fe: −0.80‰ to −0.55‰; δ
66Zn: −0.11‰ to 0.04‰) from Wang et al. [26], as well as the δ
66Zn values of sphalerite (−0.25‰ to 0.03‰) and Mn-Fe carbonate (−0.01‰) from Duan et al. [5]. Finally, we set up 12 Fe-Zn isotopic fractionation models for pyrite, sphalerite, and Mn-Fe carbonate (Figure 8) with different δ
56Fei values of 0‰ (mean value of magma) [51], −0.5‰, and −1‰, as well as δ
66
Z
n
i values of −0.28‰ (the lightest value of seafloor hydrothermal fluids) [62], 0‰ (mean value of bulk earth) [54], and 0.23‰ (the mean value of deep sea water) [63, 108, 109].Figure 8
The Fe-Zn isotopic fractionation models for the ore-forming system and minerals (pyrite, sphalerite, and Mn-Fe carbonate) with differentδ
56Fei and δ
66Zni values.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)These fractionation models show that theF ranges for ore-forming system highly depend on the δ
i values (Figure 8). The details are as follows: the pyrite covers the F ranges of 30.1
%
​
~
62.5
% (δ
56Fei = 0‰; Figure 8(a)), 48.8
%
​
~
99.0
% (δ
56Fei = −0.5‰; Figure 8(b)), and more than 77.7% (δ
56Fei = −1‰; Figure 8(c)). The F ranges for sphalerite are 38.7
%
​
~
70.9
% (δ
66Zni = −0.28‰; Figure 8(d)), 26.7
%
​
~
48.7
% (δ
66Zni = 0‰; Figure 8(e)), and 19.8
%
​
~
36.3
% (δ
66Zni = 0.23‰; Figure 8(f)), respectively. In comparison, the Mn-Fe carbonates has the F ranges of 7.1
%
​
~
24.9
% (δ
56Fei = 0‰; Figure 8(g)), 10.9
%
​
~
39.7
% (δ
56Fei = −0.5‰; Figure 8(h)), and 17.3
%
​
~
63.3
% (δ
56Fei = −1‰; Figure 8(i)) for Fe isotope, as well as 55.6
%
​
~
68.1
% (δ
66Zni = −0.28‰; Figure 8(j)), 38.2
%
​
~
46.8
% (δ
66Zni = 0‰; Figure 8(k)), and 28.5
%
​
~
34.8
% (δ
66Zni = 0.23‰; Figure 8(l)) for Zn isotope. All of these results suggest that the Fe-Zn isotopic data of Zhaxikang deposit fit these Rayleigh fractionation models well.However, we need to consider the following facts in the Zhaxikang deposit: (1) the second pulse of mineralization has brought some Fe to form the stage 3 pyrite with heavierδ
56Fe values (0.23
‰
​
~
0.43
‰) [26]; thus the stage 3 pyrite does not fit the fractionation models; (2) as most of the sphalerite and pyrite are paragenetic during the first pulse of mineralization (Figures 4(g), 4(i), 4(k), 5(a)–5(c) and 6(c)), especially in the earliest lamina (Figures 4(a)–4(c), 5(c), 6(f) and 6(m)), the sphalerite and pyrite should overlap on F ranges; (3) in theory, the Mn-Fe carbonates would have the same F range for Fe-Zn isotopes. Nonetheless, in view of the tight Zn isotopic variation range, the F range for Fe isotope should cover that for Zn isotope; (4) during the earlier period, as the ore-forming system consisting of vapour and liquid phases is dominant, F values in fractionation models should be large. Taking all these facts and 12 fractionation models into consideration, the δ
56Fei value is supposed to be in the range of -
0.5
‰
​
~−1‰, and the δ
66
Z
n
i value should be between −0.28‰ and 0‰.


## 5.3. Implications for the Genesis of Zhaxikang Deposit

### 5.3.1. Excluding the Possibility of Hot Spring Genetic Model

The hot spring model predicts that metals (e.g., Zn, Pb, Sb, Ag, and Fe) are leached from sedimentary wall rocks, which is supported by the following evidence: (1) theδ
34S values of the sulfides (4.5
‰
​
~
12
‰) are similar to those of sedimentary wall rocks (4.93
‰
​
~
11.49
‰); (2) the δ
30Si values of quartz (-
0.90
‰
​
~−0.40‰) are the same as those of siliceous rocks with hot spring genesis; (3) the δ
D
V
–
S
M
O
W (-
135
‰
​
~−127‰) and δ
18
O
H
2
O values (-
13.7
‰
​
~
12.4
‰) of fluid inclusions trapped in quartz are similar to those of the south Tibetan hot spring; (4) the Pb isotopes (P
b
206
/
P
b
204: 18.474~19.637; P
b
207
/
P
b
204: 15.649~15.774; P
b
208
/
P
b
204: 39.660~40.010) show the characteristics of radiogenic Pb; (5) the He-Ar isotopes demonstrate the contribution of crustal fluid and meteoric water [3, 4].However, this genetic model is inconsistent with textural and Fe-Zn isotopic evidence presented here. Firstly, the primary sedimentary wall rocks in the orefield are slate. As suggested by the continuous batch experimental research of Fernandez and Borrok [27], the ore-forming fluid would preferentially leach out the heavy Zn isotopes. Likewise, Chen et al. [110] have analyzed the Zn isotopic compositions of samples from 8 hot springs, and the results show most of the hot springs have relatively constant and heavier δ
66Zn values (approximately 0.42‰) than host rocks (−0.42‰ to 0.14‰). Therefore, if the metallogenic elements are leached from the slate by hot spring, there would be some sphalerite with heavier δ
66Zn values than these slates in Zhaxikang deposit. Nevertheless, the δ
66Zn values of the slate from Zhaxikang orefield range from −0.23‰ to 0.10‰ that are similar to sphalerite, especially the unmodified slate sample with the δ
66Zn value of 0.10‰ that is even a little heavier than that of sphalerite (−0.38‰ to 0.07‰; Figure 9(b)) [26]. As the slate has heavier Zn isotopic compositions than those of host rocks from Chen et al. [110], the δ
66Zn values of the hot spring in Zhaxikang orefield would even be heavier than 0.42‰, which are much heavier than the δ
66Zni values (-
0.28
‰
​
~
0
‰) of the fractionation models in Section 5.2. Secondly, Sharam et al. [46] have measured the δ
56Fe values (−0.59‰ to −0.12‰) of hot springs in Juan de Fuca Ridge, which is much heavier than the δ
56Fei value (-
0.5
‰
​
~−1‰) gained from the fractionation models. Moreover, the marine fluids usually have lighter Fe isotopic compositions (Figure 9(a)); thus the hot springs in Tibet would have heavier δ
56Fe values than those of Sharam et al. [46].Figure 9
(a) Fe isotopic compositions of pyrite, Mn-Fe carbonate, and arsenopyrite in the Zhaxikang deposit (part of the data is quoted from Wang et al. [26]). Other Fe isotopic data for the Bulk Silicate Earth [10], shale [34, 35], igneous rocks [21–23, 28, 36–45], seafloor hydrothermal fluid system [34, 46–49], midoceanic ridges pyrite [48, 50], the Bayan Obo Fe-REE magmatic-type deposit in China [51], the skarn-type deposits in Tongling ore district in China [21, 22], the Schwarzwald hydrothermal vein deposit in Germany [52], the Renison Sn-W deposit in Australia [23], and the Dongshengmiao SEDEX-type deposit in China [53] are also plotted for comparison. (b) Zn isotopic compositions of sphalerite, Mn-Fe carbonate, and slate in the Zhaxikang deposit (part of the data is quoted from Wang et al. [26]), compared with the Bulk Silicate Earth [54], sedimentary rocks [55–57], igneous rocks [28, 54, 57–60], deep-sea carbonates [61], seafloor hydrothermal fluid system [62], deep sea water [51, 54, 63], and other deposits with different geneses: the Gorno and Raibl magmatic-type deposit in Italy [14], the skarn-type deposits in the Tongling ore district in China [12], the Tianqiao and Bangbangqiao carbonated-hosted Pb-Zn sulfide deposits in China [64], the Irish-type deposit in Ireland [19], the Cévennes MVT deposit in France [65], the Alexandrinka VHMS-type deposit in Russia [18], the Red Dog SEDEX-type ore district in Alaska [20], and the Dongshengmiao SEDEX-type deposit in China [53].
(a)
(b)The possibility of hot spring genesis for the first pulse of mineralization event can be excluded by Fe-Zn isotopic data. Evidence for the second pulse of mineralization is based on the fact that the Fe-Zn isotope values do not follow similar concentric patterns with the ore textures as seen in sample ZXK-2. Meanwhile, the evidence from Si-H-O isotopes demonstrates that the second pulse of mineralization may be related to hot spring. Additionally, the Fe-Zn isotopic data demonstrate that the sedimentary wall rocks have not provided significant amounts of metals, although the S-Pb isotopic data show that these wall rocks constitute some contribution for S-Pb [4, 5], whereas the similar Zn isotopic compositions of slate and sphalerite suggest that they share the same Zn origin.


### 5.3.2. Inconsistency with the Magmatic-Hydrothermal Fluid Genetic Models

There are two genetic models for the magmatic-hydrothermal fluid genesis. In the first model, Duan et al. [5] considered that the genesis of Zhaxikang deposit relates to the mid-low temperature magma-related hydrothermal activity and that the metallogenic elements are mainly sourced from the mixing of basement and the sedimentary wall rocks. The evidence is mainly from the Zn-S-Pb isotopes that we mentioned in Section 1.Duan et al. [5] have analyzed the δ
66Zn values of the sulfides (-
0.25
‰
​
~
0.03
‰) and basement rocks (0.05
‰
​
~
0.21
‰). The dominating sedimentary wall rocks in the orefield are slate (δ
66Zn values: -
0.23
‰
​
~
0.10
‰) [26], which is formed by the epimetamorphism of shale and sandstone. Meanwhile, combining the data from Wang et al. [26] with this study, the sphalerite should have a range from −0.38‰ to 0.07‰ in Zhaxikang deposit. Both the basement and sedimentary wall rocks have heavier Zn isotopic compositions than sphalerite in Zhaxikang deposit. However, just like we discussed in Section 5.3.1, if the Zn is sourced from mixing of the basement and sedimentary wall rocks, there should be some sphalerite with heavier δ
66Zn values than these rocks. This inference is also evidenced by the research of Zhou et al. [64]: the Paleozoic carbonate host rocks and Precambrian basements are considered to be the origin of metals, and these rocks have lighter δ
66Zn values (−0.52‰ to 0.16‰) than the sphalerite from the Tianqiao (−0.54‰ to 0.30‰) and Bangbangqiao (−0.21‰ to 0.43‰) deposits in the Sichuan-Yunnan-Guizhou Pb-Zn metallogenic province (Figure 9(b)). Additionally, in respect of the Fe isotope, the Schwarzwald hydrothermal vein deposit in Germany can be used as an analogy [52]. Iron in this deposit originates from the basement consisting of granites and gneisses, as well as sedimentary rocks including shale and sandstone. The basement of Zhaxikang is composed of dolerite, quartz diorite, rhyolite porphyry, pyroclastics, and porphyritic monzogranite [5], and these crust-derived igneous rocks have the similar Fe isotopic composition with granites according to the data from previous studies [21–23, 28, 36–45]. Likewise, the δ
56Fe values of shale and sandstone in the Schwarzwald deposit are −0.21‰, 0.03‰, and 0.22‰, all of which fall into the Fe isotopic variation range of shale (-
0.39
‰
​
~
0.71
‰) from Beard et al. [34] and Rouxel et al. [35]. Therefore, we regard the fact that the metal-sourced rocks in these two deposits have the similar Fe isotopic compositions. And, yet, Markl et al. [52] suggested that fluid-rock interaction make the ore-forming fluid have the δ
56Fe value of -
0.5
‰
​
~
0
‰, which is heavier than the δ
56Fei value (-
0.5
‰
​
~−1‰) of the fractionation models in Section 5.2. Furthermore, it is generally known that the equilibrium isotope fractionation is a function of temperature, with larger fractionation generated at lower temperatures [111]. Although the temperature of ore-forming fluid (100~200°C) [52] in the Schwarzwald deposit is lower than that in Zhaxikang deposit (~250°C) [5], the lightest δ
56Fe value (−1.36‰) of siderite from this deposit is much heavier than that of Mn-Fe carbonate (−1.95‰) from Zhaxikang deposit (Figure 9(a)). All of these inconsistencies from Fe-Zn isotopic data indicate that this genetic model is not appropriate for Zhaxikang deposit.In the second genetic model, Xie et al. [6] suggested the mineralization in Zhaxikang deposit was genetically related to the Miocene dome-related magmatism. This magmatism generated the pegmatite and two-mica granite in the core of the regional domes (Figure 1(b); e.g., Cuonadong, Yalaxiangbo, Ranba, and Kangma), and the high-temperature ore-forming fluids were derived from magmatic melts exsolution. This hypothesis is based on the research of field geology, petrography, melt and fluid inclusions, and C−H−O isotopes: (1) the δ
13
C
V
–
P
D
B (-
6.1
‰
​
~−6.9‰) and δ
18
O
V
–
S
M
O
W (+
9.9
‰
​
~+11.8‰) values of rhodochrosite show the mantle origin; (2) the δ
13
D
V
–
S
M
O
W ( -
144.8
‰
​
~−110‰) and δ
18
O
V
–
S
M
O
W (-
9.85
‰
​
~−8.89‰) values, the low salinity (0.2~7.9 wt.% NaCleqv), high temperature (298~457°C), and rich CO2 content with minor CH4, N2, C2H6, C3H8, and C6H6 of the melt and fluid inclusions trapped in quartz and beryl within pegmatite indicate the magmatic origin; (3) the 40Ar-39Ar plateau age of pegmatite (18.93 ± 0.27 Ma) is similar to stage 5 quartz-pyrite-stibnite (17.9 ± 0.5 Ma).The pegmatite and two-mica granite in the domes have the typical characteristics of S-type granite: (1) containing abundant Al-rich minerals, A/CNK: 1.07~1.24; (2) w(SiO2): 73.26
%
​
~
74.33
%, w(K2O)/w(Na2O): 1.1~1.2, w(FeO)/w(Fe + Mn): 0.64~0.76; (3) the content of corundum molecules > 1% in the CIPW standard minerals [112, 113]. Meanwhile, based on the facts that all the Fe-bearing minerals in Zhaxikang deposit are ferrous, and the magmatic-hydrothermal fluid is CO2-rich with minor CH4, N2, C2H6, C3H8, and C6H6 [6], and the parent magma is most likely S-type reduced magma. In consideration of the situation in the Renison Sn-W deposit that we mentioned in Section 1, the ore-forming fluid is considered to exsolve from S-type reduced magma and the sulfides have heavier Fe isotopic compositions than ore-related igneous rocks [23]. The case in Zhaxikang deposit is contrary yet the pyrite from the first pulse of mineralization has the δ
56Fe value of -
0.33
‰
​
~
0.23
‰ that is much lighter than those of granitoids (−0.08% to 0.59%) [22, 28]. Moreover, Heimann et al. [43] proposed that we would not expect Fe isotopic compositions of high F/Cl magmatic/hydrothermal systems to significantly deviate from the average of igneous rocks, hence as Xie et al. [6] considered that the Zhaxikang parent magma has high F contents, the magmatic-hydrothermal fluid in this genetic model would have the similar Fe isotopic compositions with granitoids, which is not in line with our fractionation models. Even if the parent magma in Zhaxikang orefield is oxidized-type and is similar to the magma related to I-type granitoids with δ
56Fe values of -
0.04
%
​
~
0.59
% in Tongling ore district (Figure 9(a)) [21, 22], it is still hard to produce so light δ
56Fei value (-
0.5
‰
​
~−1‰). Similarity, Zn isotopic data do not support this genetic model, either. Chen et al. [54] studied the Zn isotopic fractionation during igneous process and suggested that the maximum Zn isotopic variation induced by high-temperature igneous processes is no larger than δ
66Zn~0.10%. Besides this, Telus et al. [28] measured the δ
66Zn values of pegmatite (0.25
‰
​
~
0.59
‰) and some other granitoids (-
0.16
‰
​
~
0.21
‰) and then found there is even no variation in δ
66Zn values during fluid exsolution in some cases. As the granitoids in the regional domes are principally dominated by pegmatite and the magmatic fluids which have high temperature of 298~457°C [6], it is also hard to generate δ
66Zni value between −0.28‰ and 0‰ as yielded from the Fe-Zn isotopic fractionation models. Consequently, the Fe-Zn isotopic data are also not in favor of the second genetic model.On the other hand, neither of these two genetic models can explain the different Fe-Zn isotopic and elemental variations in sample ZXK-2, which is considered to result from the overprint by the second pulse of mineralization. Overall, all of the Fe-Zn isotopic and elemental evidences are inconsistent with both of these magmatic-hydrothermal fluid genetic models. Nevertheless, the evidence for these two genetic models may prove that the second pulse of mineralization is related to magmatic-hydrothermal fluids.


### 5.3.3. Constraints on SEDEX Modified by Hydrothermal Fluid Genetic Model

Zheng et al. [7, 25] considered that the first pulse of mineralization (Pb-Zn) has the SEDEX genesis, and the second pulse of mineralization (Sb-Ag) is related to hot spring that overprints the earlier mineralization. The previous evidences are mainly as follows: (1) the evidence from ore textures and obviously late Sb mineralization compared to Pb-Zn mineralization in Section 2.3; (2) the stage 1 and 2 ores have high Mn, Fe, Ba, and B contents, Ga ≫ In (Ga/In: 1.49~4.47), Pb + Zn ≫ Cu, and host rocks are exhalative rocks and Mn-Fe carbonates; (3) the δ
34S values of sphalerite and galena (7.7
‰
​
~
12
‰) are different from those of stibnite (4.5
‰
200
B
~
7.1
‰); (4) the data from valentinite (δ
30Si values: -
0.9
‰
​
~−0.4‰) and fluid inclusions trapped in quartz (δ
D
V
–
S
M
O
W: -
162
‰
​
~−142‰, δ
18
O
H
2
O: -
12.9
‰
​
~−1.9‰, homogenization and freezing temperatures: 164~313°C and −3.2​
~−0.7°C, salinity: 0.7
%
​
~
5.3
%, and density: 0.74~0.93 g/cm3) both demonstrate the Sb mineralization is related to the south Tibetan hot spring; (5) the Rb-Sr isotopic isochrone age of stage 2 sphalerite is 147.2 ± 3.2 Ma that is similar to the Jurassic sedimentary wall rocks, whereas a later stage quartz-pyrite-stibnite vein has the A
r
40-A
r
39 plateau age of 17.9 ± 0.5 Ma.Wen et al. [17] investigated several Pb-Zn deposits with different geneses in China, and the results show that the Zn/Cd ratios of sphalerite vary with different geneses: (1) high-temperature systems including the porphyry, magmatic hydrothermal, skarn, and volcanic hosted massive sulfide (VMS)-type deposits: 155~223; (2) low-temperature systems that include the Mississippi Valley-type (MVT) deposits: 17~201; (3) SEDEX-type deposits of exhalative systems: 316~368; (4) seafloor hydrothermal sulfides of exhalative systems: 211~510. According to the EPMA data, the Zn/Cd ratios of sphalerite range from 296 to 399 in Zhaxikang deposit, which overlap the range of exhalative systems and much higher than those of high-temperature and low-temperature systems.The Fe-Zn isotopic data also conform to the SEDEX modified by hydrothermal fluid genetic model. Firstly, as discussed in Sections5.3.1 and 5.3.2, neither the hot spring nor the magmatic-hydrothermal fluids can have the δ
56Fei (-
0.5
‰
​
~
-
1
‰) and δ
66Zni (-
0.28
‰
​
~
0
‰) values to satisfy the Fe-Zn isotopic fractionation models. However, the Fe-Zn isotopic compositions of seafloor hydrothermal fluid system covers the range of -
1.79
‰
​
~
0.04
‰ for Fe isotope and -
0.28
‰
​
~
0.96
‰ for Zn isotope according to previous studies (Figure 9) [38, 46–49, 62], which can generate the ore-forming system with δ
56Fei and δ
66Zni values to meet the fractionation models. Secondly, although there are large overlaps in Zn isotopic compositions among deposits with different geneses, the Zn isotopic compositions of Zhaxikang deposit is most similar to the Alexandrinka VHMS-type and Red Dog SEDEX-type deposits with marine origin, as well as obviously distinguishing from the narrow range of magmatic-hydrothermal deposits (Figure 9(b)). Thirdly, the lightest δ
56Fe value in Zhaxikang deposit is −1.95‰, and only the minerals with marine origin (mid-oceanic ridges pyrite) or SEDEX genesis (sphalerite and pyrrhotite in Dongshengmiao SEDEX-type deposit) can have so light δ
56Fe values (Figure 9(a)). Fourthly, under the conditions of high P
C
O
2 and low PH (<8), the Zn precipitated as sulfides is isotopically nearly unfractionated with respect to the parent hydrothermal fluid, whereas, under the conditions of high P
C
O
2 and high-PH (>9), negative δ
66Zn values down to 0.6‰ can be expected in sulfides precipitated from hydrothermal fluid [107]. In Zhaxikang deposit, P
C
O
2 would be high as there are plenty of Mn-Fe carbonates; meanwhile, owing to the facts that the modern seawater has the PH around 8 and there is more CO2 in Jurassic atmosphere than present, the Jurassic seawater would have a lower PH than modern seawater. This can well explain that the δ
66Zn values of sphalerite (-
0.38
‰
​
~
0.07
‰) slightly fractionate with the ore-forming system (δ
66Zni value: -
0.28
‰
​
~
0
‰). Besides these, the overprinting of earlier ores by second pulse of mineralization is not only proved by the different Fe-Zn isotopic and elemental variations in sample ZXK-2 from the other 3 samples in this research (Figure 5) but also evidenced by the temporally increasing δ
56Fe and decreasing δ
66Zn values recorded in this deposit that coincided with an increase in alteration [26]. Nevertheless, further research is required to confirm whether this hydrothermal fluid is hot spring or magmatic-hydrothermal fluid.To sum up, among the various genetic models, the Fe-Zn isotopic and EPMA evidence indicate the SEDEX modified by hydrothermal fluid genetic model is the most plausible. Our research also demonstrates that the Fe-Zn isotopes have the potential to trace the metal source and provide insights into ore-forming processes.


## 5.3.1. Excluding the Possibility of Hot Spring Genetic Model

The hot spring model predicts that metals (e.g., Zn, Pb, Sb, Ag, and Fe) are leached from sedimentary wall rocks, which is supported by the following evidence: (1) theδ
34S values of the sulfides (4.5
‰
​
~
12
‰) are similar to those of sedimentary wall rocks (4.93
‰
​
~
11.49
‰); (2) the δ
30Si values of quartz (-
0.90
‰
​
~−0.40‰) are the same as those of siliceous rocks with hot spring genesis; (3) the δ
D
V
–
S
M
O
W (-
135
‰
​
~−127‰) and δ
18
O
H
2
O values (-
13.7
‰
​
~
12.4
‰) of fluid inclusions trapped in quartz are similar to those of the south Tibetan hot spring; (4) the Pb isotopes (P
b
206
/
P
b
204: 18.474~19.637; P
b
207
/
P
b
204: 15.649~15.774; P
b
208
/
P
b
204: 39.660~40.010) show the characteristics of radiogenic Pb; (5) the He-Ar isotopes demonstrate the contribution of crustal fluid and meteoric water [3, 4].However, this genetic model is inconsistent with textural and Fe-Zn isotopic evidence presented here. Firstly, the primary sedimentary wall rocks in the orefield are slate. As suggested by the continuous batch experimental research of Fernandez and Borrok [27], the ore-forming fluid would preferentially leach out the heavy Zn isotopes. Likewise, Chen et al. [110] have analyzed the Zn isotopic compositions of samples from 8 hot springs, and the results show most of the hot springs have relatively constant and heavier δ
66Zn values (approximately 0.42‰) than host rocks (−0.42‰ to 0.14‰). Therefore, if the metallogenic elements are leached from the slate by hot spring, there would be some sphalerite with heavier δ
66Zn values than these slates in Zhaxikang deposit. Nevertheless, the δ
66Zn values of the slate from Zhaxikang orefield range from −0.23‰ to 0.10‰ that are similar to sphalerite, especially the unmodified slate sample with the δ
66Zn value of 0.10‰ that is even a little heavier than that of sphalerite (−0.38‰ to 0.07‰; Figure 9(b)) [26]. As the slate has heavier Zn isotopic compositions than those of host rocks from Chen et al. [110], the δ
66Zn values of the hot spring in Zhaxikang orefield would even be heavier than 0.42‰, which are much heavier than the δ
66Zni values (-
0.28
‰
​
~
0
‰) of the fractionation models in Section 5.2. Secondly, Sharam et al. [46] have measured the δ
56Fe values (−0.59‰ to −0.12‰) of hot springs in Juan de Fuca Ridge, which is much heavier than the δ
56Fei value (-
0.5
‰
​
~−1‰) gained from the fractionation models. Moreover, the marine fluids usually have lighter Fe isotopic compositions (Figure 9(a)); thus the hot springs in Tibet would have heavier δ
56Fe values than those of Sharam et al. [46].Figure 9
(a) Fe isotopic compositions of pyrite, Mn-Fe carbonate, and arsenopyrite in the Zhaxikang deposit (part of the data is quoted from Wang et al. [26]). Other Fe isotopic data for the Bulk Silicate Earth [10], shale [34, 35], igneous rocks [21–23, 28, 36–45], seafloor hydrothermal fluid system [34, 46–49], midoceanic ridges pyrite [48, 50], the Bayan Obo Fe-REE magmatic-type deposit in China [51], the skarn-type deposits in Tongling ore district in China [21, 22], the Schwarzwald hydrothermal vein deposit in Germany [52], the Renison Sn-W deposit in Australia [23], and the Dongshengmiao SEDEX-type deposit in China [53] are also plotted for comparison. (b) Zn isotopic compositions of sphalerite, Mn-Fe carbonate, and slate in the Zhaxikang deposit (part of the data is quoted from Wang et al. [26]), compared with the Bulk Silicate Earth [54], sedimentary rocks [55–57], igneous rocks [28, 54, 57–60], deep-sea carbonates [61], seafloor hydrothermal fluid system [62], deep sea water [51, 54, 63], and other deposits with different geneses: the Gorno and Raibl magmatic-type deposit in Italy [14], the skarn-type deposits in the Tongling ore district in China [12], the Tianqiao and Bangbangqiao carbonated-hosted Pb-Zn sulfide deposits in China [64], the Irish-type deposit in Ireland [19], the Cévennes MVT deposit in France [65], the Alexandrinka VHMS-type deposit in Russia [18], the Red Dog SEDEX-type ore district in Alaska [20], and the Dongshengmiao SEDEX-type deposit in China [53].
(a)
(b)The possibility of hot spring genesis for the first pulse of mineralization event can be excluded by Fe-Zn isotopic data. Evidence for the second pulse of mineralization is based on the fact that the Fe-Zn isotope values do not follow similar concentric patterns with the ore textures as seen in sample ZXK-2. Meanwhile, the evidence from Si-H-O isotopes demonstrates that the second pulse of mineralization may be related to hot spring. Additionally, the Fe-Zn isotopic data demonstrate that the sedimentary wall rocks have not provided significant amounts of metals, although the S-Pb isotopic data show that these wall rocks constitute some contribution for S-Pb [4, 5], whereas the similar Zn isotopic compositions of slate and sphalerite suggest that they share the same Zn origin.


## 5.3.2. Inconsistency with the Magmatic-Hydrothermal Fluid Genetic Models

There are two genetic models for the magmatic-hydrothermal fluid genesis. In the first model, Duan et al. [5] considered that the genesis of Zhaxikang deposit relates to the mid-low temperature magma-related hydrothermal activity and that the metallogenic elements are mainly sourced from the mixing of basement and the sedimentary wall rocks. The evidence is mainly from the Zn-S-Pb isotopes that we mentioned in Section 1.Duan et al. [5] have analyzed the δ
66Zn values of the sulfides (-
0.25
‰
​
~
0.03
‰) and basement rocks (0.05
‰
​
~
0.21
‰). The dominating sedimentary wall rocks in the orefield are slate (δ
66Zn values: -
0.23
‰
​
~
0.10
‰) [26], which is formed by the epimetamorphism of shale and sandstone. Meanwhile, combining the data from Wang et al. [26] with this study, the sphalerite should have a range from −0.38‰ to 0.07‰ in Zhaxikang deposit. Both the basement and sedimentary wall rocks have heavier Zn isotopic compositions than sphalerite in Zhaxikang deposit. However, just like we discussed in Section 5.3.1, if the Zn is sourced from mixing of the basement and sedimentary wall rocks, there should be some sphalerite with heavier δ
66Zn values than these rocks. This inference is also evidenced by the research of Zhou et al. [64]: the Paleozoic carbonate host rocks and Precambrian basements are considered to be the origin of metals, and these rocks have lighter δ
66Zn values (−0.52‰ to 0.16‰) than the sphalerite from the Tianqiao (−0.54‰ to 0.30‰) and Bangbangqiao (−0.21‰ to 0.43‰) deposits in the Sichuan-Yunnan-Guizhou Pb-Zn metallogenic province (Figure 9(b)). Additionally, in respect of the Fe isotope, the Schwarzwald hydrothermal vein deposit in Germany can be used as an analogy [52]. Iron in this deposit originates from the basement consisting of granites and gneisses, as well as sedimentary rocks including shale and sandstone. The basement of Zhaxikang is composed of dolerite, quartz diorite, rhyolite porphyry, pyroclastics, and porphyritic monzogranite [5], and these crust-derived igneous rocks have the similar Fe isotopic composition with granites according to the data from previous studies [21–23, 28, 36–45]. Likewise, the δ
56Fe values of shale and sandstone in the Schwarzwald deposit are −0.21‰, 0.03‰, and 0.22‰, all of which fall into the Fe isotopic variation range of shale (-
0.39
‰
​
~
0.71
‰) from Beard et al. [34] and Rouxel et al. [35]. Therefore, we regard the fact that the metal-sourced rocks in these two deposits have the similar Fe isotopic compositions. And, yet, Markl et al. [52] suggested that fluid-rock interaction make the ore-forming fluid have the δ
56Fe value of -
0.5
‰
​
~
0
‰, which is heavier than the δ
56Fei value (-
0.5
‰
​
~−1‰) of the fractionation models in Section 5.2. Furthermore, it is generally known that the equilibrium isotope fractionation is a function of temperature, with larger fractionation generated at lower temperatures [111]. Although the temperature of ore-forming fluid (100~200°C) [52] in the Schwarzwald deposit is lower than that in Zhaxikang deposit (~250°C) [5], the lightest δ
56Fe value (−1.36‰) of siderite from this deposit is much heavier than that of Mn-Fe carbonate (−1.95‰) from Zhaxikang deposit (Figure 9(a)). All of these inconsistencies from Fe-Zn isotopic data indicate that this genetic model is not appropriate for Zhaxikang deposit.In the second genetic model, Xie et al. [6] suggested the mineralization in Zhaxikang deposit was genetically related to the Miocene dome-related magmatism. This magmatism generated the pegmatite and two-mica granite in the core of the regional domes (Figure 1(b); e.g., Cuonadong, Yalaxiangbo, Ranba, and Kangma), and the high-temperature ore-forming fluids were derived from magmatic melts exsolution. This hypothesis is based on the research of field geology, petrography, melt and fluid inclusions, and C−H−O isotopes: (1) the δ
13
C
V
–
P
D
B (-
6.1
‰
​
~−6.9‰) and δ
18
O
V
–
S
M
O
W (+
9.9
‰
​
~+11.8‰) values of rhodochrosite show the mantle origin; (2) the δ
13
D
V
–
S
M
O
W ( -
144.8
‰
​
~−110‰) and δ
18
O
V
–
S
M
O
W (-
9.85
‰
​
~−8.89‰) values, the low salinity (0.2~7.9 wt.% NaCleqv), high temperature (298~457°C), and rich CO2 content with minor CH4, N2, C2H6, C3H8, and C6H6 of the melt and fluid inclusions trapped in quartz and beryl within pegmatite indicate the magmatic origin; (3) the 40Ar-39Ar plateau age of pegmatite (18.93 ± 0.27 Ma) is similar to stage 5 quartz-pyrite-stibnite (17.9 ± 0.5 Ma).The pegmatite and two-mica granite in the domes have the typical characteristics of S-type granite: (1) containing abundant Al-rich minerals, A/CNK: 1.07~1.24; (2) w(SiO2): 73.26
%
​
~
74.33
%, w(K2O)/w(Na2O): 1.1~1.2, w(FeO)/w(Fe + Mn): 0.64~0.76; (3) the content of corundum molecules > 1% in the CIPW standard minerals [112, 113]. Meanwhile, based on the facts that all the Fe-bearing minerals in Zhaxikang deposit are ferrous, and the magmatic-hydrothermal fluid is CO2-rich with minor CH4, N2, C2H6, C3H8, and C6H6 [6], and the parent magma is most likely S-type reduced magma. In consideration of the situation in the Renison Sn-W deposit that we mentioned in Section 1, the ore-forming fluid is considered to exsolve from S-type reduced magma and the sulfides have heavier Fe isotopic compositions than ore-related igneous rocks [23]. The case in Zhaxikang deposit is contrary yet the pyrite from the first pulse of mineralization has the δ
56Fe value of -
0.33
‰
​
~
0.23
‰ that is much lighter than those of granitoids (−0.08% to 0.59%) [22, 28]. Moreover, Heimann et al. [43] proposed that we would not expect Fe isotopic compositions of high F/Cl magmatic/hydrothermal systems to significantly deviate from the average of igneous rocks, hence as Xie et al. [6] considered that the Zhaxikang parent magma has high F contents, the magmatic-hydrothermal fluid in this genetic model would have the similar Fe isotopic compositions with granitoids, which is not in line with our fractionation models. Even if the parent magma in Zhaxikang orefield is oxidized-type and is similar to the magma related to I-type granitoids with δ
56Fe values of -
0.04
%
​
~
0.59
% in Tongling ore district (Figure 9(a)) [21, 22], it is still hard to produce so light δ
56Fei value (-
0.5
‰
​
~−1‰). Similarity, Zn isotopic data do not support this genetic model, either. Chen et al. [54] studied the Zn isotopic fractionation during igneous process and suggested that the maximum Zn isotopic variation induced by high-temperature igneous processes is no larger than δ
66Zn~0.10%. Besides this, Telus et al. [28] measured the δ
66Zn values of pegmatite (0.25
‰
​
~
0.59
‰) and some other granitoids (-
0.16
‰
​
~
0.21
‰) and then found there is even no variation in δ
66Zn values during fluid exsolution in some cases. As the granitoids in the regional domes are principally dominated by pegmatite and the magmatic fluids which have high temperature of 298~457°C [6], it is also hard to generate δ
66Zni value between −0.28‰ and 0‰ as yielded from the Fe-Zn isotopic fractionation models. Consequently, the Fe-Zn isotopic data are also not in favor of the second genetic model.On the other hand, neither of these two genetic models can explain the different Fe-Zn isotopic and elemental variations in sample ZXK-2, which is considered to result from the overprint by the second pulse of mineralization. Overall, all of the Fe-Zn isotopic and elemental evidences are inconsistent with both of these magmatic-hydrothermal fluid genetic models. Nevertheless, the evidence for these two genetic models may prove that the second pulse of mineralization is related to magmatic-hydrothermal fluids.


## 5.3.3. Constraints on SEDEX Modified by Hydrothermal Fluid Genetic Model

Zheng et al. [7, 25] considered that the first pulse of mineralization (Pb-Zn) has the SEDEX genesis, and the second pulse of mineralization (Sb-Ag) is related to hot spring that overprints the earlier mineralization. The previous evidences are mainly as follows: (1) the evidence from ore textures and obviously late Sb mineralization compared to Pb-Zn mineralization in Section 2.3; (2) the stage 1 and 2 ores have high Mn, Fe, Ba, and B contents, Ga ≫ In (Ga/In: 1.49~4.47), Pb + Zn ≫ Cu, and host rocks are exhalative rocks and Mn-Fe carbonates; (3) the δ
34S values of sphalerite and galena (7.7
‰
​
~
12
‰) are different from those of stibnite (4.5
‰
200
B
~
7.1
‰); (4) the data from valentinite (δ
30Si values: -
0.9
‰
​
~−0.4‰) and fluid inclusions trapped in quartz (δ
D
V
–
S
M
O
W: -
162
‰
​
~−142‰, δ
18
O
H
2
O: -
12.9
‰
​
~−1.9‰, homogenization and freezing temperatures: 164~313°C and −3.2​
~−0.7°C, salinity: 0.7
%
​
~
5.3
%, and density: 0.74~0.93 g/cm3) both demonstrate the Sb mineralization is related to the south Tibetan hot spring; (5) the Rb-Sr isotopic isochrone age of stage 2 sphalerite is 147.2 ± 3.2 Ma that is similar to the Jurassic sedimentary wall rocks, whereas a later stage quartz-pyrite-stibnite vein has the A
r
40-A
r
39 plateau age of 17.9 ± 0.5 Ma.Wen et al. [17] investigated several Pb-Zn deposits with different geneses in China, and the results show that the Zn/Cd ratios of sphalerite vary with different geneses: (1) high-temperature systems including the porphyry, magmatic hydrothermal, skarn, and volcanic hosted massive sulfide (VMS)-type deposits: 155~223; (2) low-temperature systems that include the Mississippi Valley-type (MVT) deposits: 17~201; (3) SEDEX-type deposits of exhalative systems: 316~368; (4) seafloor hydrothermal sulfides of exhalative systems: 211~510. According to the EPMA data, the Zn/Cd ratios of sphalerite range from 296 to 399 in Zhaxikang deposit, which overlap the range of exhalative systems and much higher than those of high-temperature and low-temperature systems.The Fe-Zn isotopic data also conform to the SEDEX modified by hydrothermal fluid genetic model. Firstly, as discussed in Sections5.3.1 and 5.3.2, neither the hot spring nor the magmatic-hydrothermal fluids can have the δ
56Fei (-
0.5
‰
​
~
-
1
‰) and δ
66Zni (-
0.28
‰
​
~
0
‰) values to satisfy the Fe-Zn isotopic fractionation models. However, the Fe-Zn isotopic compositions of seafloor hydrothermal fluid system covers the range of -
1.79
‰
​
~
0.04
‰ for Fe isotope and -
0.28
‰
​
~
0.96
‰ for Zn isotope according to previous studies (Figure 9) [38, 46–49, 62], which can generate the ore-forming system with δ
56Fei and δ
66Zni values to meet the fractionation models. Secondly, although there are large overlaps in Zn isotopic compositions among deposits with different geneses, the Zn isotopic compositions of Zhaxikang deposit is most similar to the Alexandrinka VHMS-type and Red Dog SEDEX-type deposits with marine origin, as well as obviously distinguishing from the narrow range of magmatic-hydrothermal deposits (Figure 9(b)). Thirdly, the lightest δ
56Fe value in Zhaxikang deposit is −1.95‰, and only the minerals with marine origin (mid-oceanic ridges pyrite) or SEDEX genesis (sphalerite and pyrrhotite in Dongshengmiao SEDEX-type deposit) can have so light δ
56Fe values (Figure 9(a)). Fourthly, under the conditions of high P
C
O
2 and low PH (<8), the Zn precipitated as sulfides is isotopically nearly unfractionated with respect to the parent hydrothermal fluid, whereas, under the conditions of high P
C
O
2 and high-PH (>9), negative δ
66Zn values down to 0.6‰ can be expected in sulfides precipitated from hydrothermal fluid [107]. In Zhaxikang deposit, P
C
O
2 would be high as there are plenty of Mn-Fe carbonates; meanwhile, owing to the facts that the modern seawater has the PH around 8 and there is more CO2 in Jurassic atmosphere than present, the Jurassic seawater would have a lower PH than modern seawater. This can well explain that the δ
66Zn values of sphalerite (-
0.38
‰
​
~
0.07
‰) slightly fractionate with the ore-forming system (δ
66Zni value: -
0.28
‰
​
~
0
‰). Besides these, the overprinting of earlier ores by second pulse of mineralization is not only proved by the different Fe-Zn isotopic and elemental variations in sample ZXK-2 from the other 3 samples in this research (Figure 5) but also evidenced by the temporally increasing δ
56Fe and decreasing δ
66Zn values recorded in this deposit that coincided with an increase in alteration [26]. Nevertheless, further research is required to confirm whether this hydrothermal fluid is hot spring or magmatic-hydrothermal fluid.To sum up, among the various genetic models, the Fe-Zn isotopic and EPMA evidence indicate the SEDEX modified by hydrothermal fluid genetic model is the most plausible. Our research also demonstrates that the Fe-Zn isotopes have the potential to trace the metal source and provide insights into ore-forming processes.


## 6. Conclusions

The EPMA and Fe-Zn isotopic data allow us to make the following conclusions:(1)
The ore-forming elements are transported by the ore-forming system that is the mixture of vapour and liquid phases; the vapour-liquid partitioning and mineral precipitation are the main cause of Fe-Zn isotopic and elemental variations.
(2)
The Fe-Zn isotopic fractionation models demonstrate that theδ
56Fei and δ
66Zni values of the ore-forming system are in the range of -
0.5
‰
​
~−1
‰ and −0.28‰~0‰, respectively.
(3)
Based on the evidence from the EPMA data, Fe-Zn isotopic characteristics, and fractionation models, the SEDEX modified by hydrothermal fluid genetic model is most plausible for the Zhaxikang deposit.
(4)
There are two pulses of mineralization in the Zhaxikang deposit; the overprint by the second pulse of mineralization has also partly modified the Fe-Zn isotopic and elemental compositions of some earlier samples.



---

*Source: 2197891-2018-03-20.xml*