Genesis of the Zhijiadi Ag-Pb-Zn Deposit , Central North China Craton : Constraints from Fluid Inclusions and Stable Isotope Data

The Zhijiadi Ag-Pb-Zn deposit is located in the central North China Craton. Fluid inclusions (FIs) studies indicate three types of FIs, including aqueous, aqueous-carbonic, and daughter mineral-bearing multiphase inclusions. The daughter minerals in FIs are mainly composed of marcasite, chalcopyrite, calcite, and dolomite. Microthermometric data show that the homogenization temperature and salinity of FIs decrease gradually from early to late stages. Homogenization temperatures from early to main to late stages span from 244 to 334C, from 164 to 298C, and from 111 to 174C, respectively, while their salinities are 4.0–9.9 wt.% NaCl equiv., 0.5–12.7 wt.% NaCl equiv., and 0.2–8.8 wt.% NaCl equiv., respectively. Trapping pressures drop from 203–299MPa (the early stage) to 32–158MPa (the main stage). The dropping of pressure and temperature and mixing and/or dilution of oreforming fluids result in the formation of ore deposit. Combined with C-O-S-Pb isotopic compositions, the initial ore-forming fluids andmaterials were likely derived from amagmatic system. As a whole, we proposed that this deposit belongs to medium-low temperature hydrothermal deposit related to volcanic and subvolcanic magmatism strictly controlled by the fault zones.


Introduction
The North China Craton (NCC), a major Precambrian craton in China, has been in focus in recent years with regard to craton destruction and the mechanism models [1][2][3].There is consensus that the NCC was completed by the amalgamation of the western and eastern blocks bounded by the central zone named as the Trans-North China Orogen (TNCO) (Figure 1(a)) [2,[4][5][6].Previous studies on ore deposits in the NCC are mainly focused along the margins of the craton, such as the Jiao-Liao gold province (eastern margin), Xiaoqinling gold province (southern margin), and Yan-Liao gold mineralization belt (northern margin) [7,8].There are only few systematic investigations on the mineral deposits in the central NCC, especially for the Wutai-Hengshan terrain [9,10].Most of the deposits in the Wutai-Hengshan region were related to Mesozoic magmatic-hydrothermal activity (e.g., [9]), and the types of deposits are primarily skarn, hydrothermal, and auriferous quartz vein deposit.The deposits discovered in the Wutai-Hengshan region include mainly Au, Cu, Ag, Pb, Zn, and Mo deposits, including the Diaoquan Ag-Cu deposit [10], the Zhijiadi Ag-Pb-Zn deposit [11], the Yixingzhai Au deposit [9], and the Boqiang Cu-Mo deposit [12].
ore-forming fluids were with medium-low temperature and low salinity on the basis of microthermometric data from only the aqueous (W-type) fluid inclusions [14,19,21,22].Only Wang et al. [23] reported a few microthermometric data of CO 2 -bearing fluid inclusions (C-type) in the quartz veins, but the petrographic features of this type fluid inclusion are unclear.Fluid inclusions are the "fossils" of ore-forming fluid system, and only primary FIs can be used to determine the genesis of the original fluids and the metallogenic mechanism of ore deposits (e.g., [24]).The study of CO 2bearing fluid inclusions can provide an effective way to calculate the trapping pressure and to trace the ore-forming processes (e.g., [25]).There is no systematic petrographic and microthermometric data for the C-type FIs in the Zhijiadi deposit.In this study, we observed abundant CO 2 -bearing FIs and daughter mineral-bearing FIs in fluorite and sphalerite, and hence here we present a detailed fluid inclusion study from different minerals to find out the nature of the initial ore-forming fluids of the Zhijiadi deposit and to reconstruct the temperature-pressure condition for the mineralization.In addition, based on the geologic characteristics and isotopic compositions (C-O-S-Pb), we try to trace the ore material sources and to unravel the ore-forming conditions and ore genesis.

Regional Geology
The Wutai-Hengshan region is located at the central (NCC), to the west of the Taihang Mountains.The Precambrian basement rocks at this region from south to north are termed as the Fuping (FP) Complex, Wutai (WT) Complex, and Hengshan (HS) Complex, respectively (Figure 1(b)) [9].The FP and HS Complexes are composed of 2.7-2.5 Ga amphibolite, granitic gneisses, mafic granulites, and metasedimentary rocks in granulite facies [26].The WT Complex consists of Late Archean to Paleoproterozoic granitoid plutons, metamorphosed volcanic and sedimentary rocks (Figure 1(b)), which could be subdivided into three segments, namely, the lower, middle, and upper Wutai subgroups [27][28][29].The lower part of the WT Complex is mainly composed of amphibolite facies peridotites, oceanic tholeiite, cherts, BIF (banded iron formation), sandstones, shales, calc-silicate rocks, siltstones, and minor limestones [29,30].The middle part of the WT Complex contains felsic volcanic rocks and tholeiitic basalts which are metamorphosed to greenschist facies [30].The upper part of the WT Complex mainly comprises quartz wackes, conglomerates, and mafic to intermediate and felsic volcanic rocks, metamorphosed into lower greenschist facies (Figure 1(b)).The WT Complex is considered to be unconformably overlain and also structurally interleaved by the Paleoproterozoic Huotuo Group metasedimentary rocks [29].During Mesoproterozoic to Paleozoic era, the Wutai-Hengshan region also witnessed a prolonged stage of quiescence with deposition of sedimentary cover rocks on the metamorphic basement [12].In this region, there are two intense stages of tectonic activities.The early stage NNW-and EW-trending faults are associated with Mesoarchean to Paleoproterozoic tectonism.The later stage is a series of NE-and NW-trending faults as a result of the compression of NW-SE trending horizontal stress in the Yanshan tectonic movements, which are associated with the Mesozoic magmatism and mineralization [31].During the Late Jurassic to Early Cretaceous, the tectonic setting in this region changed from early compressional setting gradually into later extensional setting [32].
The magmatic activities are very extensive during Mesozoic.Mesozoic volcanic rocks are mainly distributed in the eastern part of the Wutai-Hengshan region, which are mainly composed of tuff, dacite, and rhyolite, accompanying a series of pyroclastic rocks.There are minor Proterozoic granite and diabase as dikes and Mesozoic intermediatefelsic hypabyssal rocks in the region.These intermediatefelsic hypabyssal rocks are mainly composed of quartz porphyry, granite porphyry, granite, felsite, granodiorite, quartz monzonite, and diorite porphyry.Widespread magmatic and tectonic activities also lead to hydrothermal mineralization and the formation of volcanic-sedimentary basins.Up to now, numerous Au, Ag, Cu, Pb, Zn, and Mo ore deposits are mostly distributed in and around the intermediate-felsic hypabyssal rocks along the NW-trending faults.The timing of magmatism and metallogeny in the Wutai-Hengshan region is concentrated at 130-142 Ma based on the LA-ICP-MS zircon U-Pb ages (e.g., [12]).

Deposit Geology
The Zhijiadi deposit is located at the center of the Wutai-Hengshan region in the central NCC.The strata outcropped in the mining area are composed of basement and cover rocks (Figure 2(a)).The basement comprises the Paleoproterozoic metamorphic rocks of the Wutai Group that is mainly made up of gneiss, amphibolite, leptite, and mica quartz schist.The cover rocks in the area include Neoproterozoic carbonate rocks, Paleozoic carbonate, and clastic sedimentary rocks.Most of the area in the studied region is occupied by Later Jurassic pyroclastic rocks, and they could be subdivided into the Baiqi Formation and the Zhangjiakou Formation.The Baiqi Formation is composed of andesitic, dacitic, and tuffaceous breccia, while the Zhangjiakou Formation consists of volcanic suite including rhyolite, rhyolitic, and tuffaceous breccias.
Only a small number of the Paleoproterozoic granites and Early Cretaceous subvolcanic rocks are distributed in the northwest of the mining area.The Early Cretaceous subvolcanic rocks include quartz porphyry and granite porphyry with a small amount of sulfides disseminated in these rocks observed in the drill cores, such as sphalerite and galena.The quartz porphyry and granite porphyry show typical porphyritic texture.Plagioclase and quartz are the main phenocrysts in quartz porphyry, while K-feldspar and quartz occur as phenocrysts in granite porphyry.Pyrite, zircon, and apatite are the main accessory minerals in both of them.The quartz porphyry and granite porphyry yield LA-ICP-MS zircon U-Pb ages of 135.4 ± 1.3 Ma and 132.1 ± 1.8 Ma (our unpublished data) and its geochemistry indicates that it belongs to the highly potassic and strongly peraluminous series [17].
The occurrence of orebodies in the Zhijiadi deposit is strictly controlled by a series of NW-, EW-, and NE-trending faults (Figure 2(b)).The F 1 fault is an EW-trending tensionshear fault in the northern part of mining area, which inclines 62-80 ∘ to southeast and extends up to 1800 m.The F 2 fault belongs to a NW-trending compression-shear fault, which is approximately 1300 meters long and 2-3 meters wide at the central part of the mining area, with the inclination of the southwest at an angle of 50-70 ∘ .The NNW-trending F 3 fault is 1500 meters long and approximately 2 meters wide in the southwest side of the mining area.It is a tension-shear fault and dips at 50-70 ∘ toward the southwest.Among these faults, the F 2 fault is the most important ore-controlling structure (Figure 3), along which abundant hydrothermal breccias and cataclastic rocks are well developed (Figure 3(a)).In plain view, the hydrothermal breccia has a spindle or elliptical shape with long axes of approximately 600 m and vertical length of less than 200 m [11].The breccia fragments are composed of quartz porphyry, tuff, and dacite, with weak chloritization.
Twelve orebodies have been identified in the mining area.Most of orebodies are restricted to distribution along the F 2 fault zone, and their thickness varies from 1 to 80 m along the northeastern direction with an average thickness at 8 m (Figure 3(b)).The largest orebody (number A) expands more than 1000 m in length and 0.5-10 m width, and the average grades of which are 270 g/t Ag, 1.2% Pb, and 1.7% Zn [15].Recent exploration survey shows that the orebodies do not pinch in the controlled depth, suggesting this deposit has great potential for further exploration at depth [13].The distribution of orebodies in the Zhijiadi deposit (modified from [14]).
The hydrothermal alterations associated with Ag-Pb-Zn mineralization in the Zhijiadi deposit mainly include silicification, carbonation, sericitization, and chloritization.On the basis of mineral assemblages and the cross-cutting relationship, the mineralization process at Zhijiadi deposit can be divided into primary hydrothermal stages and a supergene epoch (Figure 5).The primary hydrothermal stages could be subdivided into three stages, namely, (1) the early stage characterized by the assemblages of quartz + pyrite veins (Figure 6  Figure 3: (a) Simplified level plan 1320 (modified from [13]); (b) geological cross section of number 5 exploration of the Zhijiadi deposit (modified from [17]).

Analytical Methods
stages were examined under transmitted light to identify the phase, shape, size, and distribution.Eighteen of doublepolished thin sections were selected for microthermometric measurements and laser Raman spectroscopic analyses.Fluid inclusion petrographic studies were carried out in National Demonstration Center for Experimental Mineral Exploration Education, China University of Geosciences, Wuhan.Microthermometric measurements were conducted by a Linkam THMSG 600 heating-freezing stage (from −196 to 550 ∘ C) mounted on a Leica DM2700P microscope at the China University of Geosciences, Wuhan.The heatingfreezing stages calibration is carried out at −56.6 ∘ C, 0.0 ∘ C, and 374.1 ∘ C using the standard of synthetic FIs supplied by FLUID INC, USA.During heating/freezing runs, the heating/freezing rate was generally constrained at 0.2-5 ∘ C/min but was reduced to <0.2 ∘ C/min near the phase transformation.The estimated accuracies of the heating and freezing measurements were ±0.5 ∘ C, ±0.2 ∘ C, and ±2 ∘ C for runs in the range of −120 ∘ C to −70 ∘ C, −70 ∘ C to +100 ∘ C, and +100 ∘ C to +600 ∘ C, respectively.Ice-melting temperatures were obtained at a heating rate of less than 0.1 ∘ C/min and homogenization temperatures at a rate of ≤1 ∘ C/min.Salinities were calculated using the equation of Hall et al. [40] for aqueous inclusions and of Roedder [41] for aqueouscarbonic inclusions.Densities (g/cm 3 ) and pressures were estimated using the Flincor software package according to microthermometric data for H 2 O-NaCl-CO 2 system [42].
Laser Raman spectroscopic analyses of FIs were carried out at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan.The analyses were conducted with a DXR laser Raman spectrometer manufactured by the Thermo Fisher Scientific Company, using a 514.5 nm Ar + ion laser as the excitation source operating at 20 mW.The exposure time was 20-40 s and ranged from 50 to 3500 cm −1 .Spectrograph aperture was defined as 50 m pinhole.Instrumental settings were kept constant during all analyses.

Carbon and Oxygen Isotope
Analyses.Rhodochrosite generally coexists with fluorite and sphalerite in the main stage (Figure 4).Eight samples are collected from the orebodies in level 1320.Carbon and oxygen isotope analyses were conducted at Key Laboratory of Tectonics and Petroleum Resources, China University of Geosciences, Wuhan.The carbon and oxygen isotopic compositions were measured as follows: dry samples were powdered and reacted overnight with 100% phosphoric acid at 72 ∘ C to extract CO 2 .The trapped CO 2 was purified and analyzed on a MAT-253 isotope ratio mass spectrometer for the carbon and oxygen isotopic compositions.Isotopic values are reported to Vienna Pee Dee Belemnite standard (V-PDB), but oxygen isotopes have been converted to V-SMOW standard, using the formula  18 O V-SMOW = 1.03086 ×  18 O V-PDB + 30.86 [43].analyses.Minerals include pyrite, sphalerite, galena, and sphalerite covering the early and main stages.
In situ sulfur isotope analyses on sulfides were conducted by the Nu Plasma II MC-ICP-MS equipped with the Resonetics-S155 excimer ArF laser ablation system at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan.The diameter of laser beam was set to 33 m with a laser repetition rate of 10 Hz, and the ablation process was set to last for 40 s.High-purity helium gas, mixing with Ar gas and a small amount of N 2 gas, was used to transport the ablated materials into the mass spectrometer.Standard-sample bracket method was used to determine the  34 S values of samples throughout the whole MC-ICP-MS analytical process.The true sulfur isotope ratio was calculated by the correction for instrumental mass bias by linear interpolation between the biases calculated from two neighboring standard-sample analyses.The primary standards used in this study were NBS-123 for sphalerite and WS-1 (in-house standard) for pyrite, galena, and chalcopyrite.The  34 S V-CDT value for NBS-123 is recommended to be 17.0‰ [44], and the  34 S V-CDT values for WS-1 natural pyrite were determined to be 0.9‰ through the isotope ratio mass spectrometer, as well as 1.1 ± 0.2‰ by SIMS at the Chinese Academy of Geochemistry, Guangzhou [45].The  34 S V-CDT value of NBS-123 calibrated through the WS-1 is 17.0 ± 0.5‰, which indicated that the matrix effect of different sulfides is not obvious.The sulfur isotopic data are reported in delta notation (‰) relative to Vienna Canyon Diablo Troilite (V-CDT).The analytical precision (1) was about ±0.1 per mil and detailed analytical conditions and procedures can be seen in [45,46].

Lead Isotope Analysis.
The lead isotopic compositions of four sphalerite, ten galena samples, and four quartz porphyry samples have been analyzed on a Neptune Plus MC-ICP-MS at State Key Laboratory for Mineral Deposits Research, Nanjing University.Sulfide samples were completely dissolved by HNO 3 and purified by cation-exchange resin.The standard errors of data were calculated by the software of Neptune through analyses of 40 data types (4 blocks, 10 cycles).In 200 ppb Pb solution, the signal intensity of 208 Pb = 6V, whereas in the blank solution 208 Pb = 20 mv.Tl NBS997 ( 203 Tl/ 205 Tl = 0.418922) is added to sample solution for correcting mass bias of Pb isotope ratios, with analytical reproducibility of 0.1% (2) for 206 Pb/ 204 Pb, 207 Pb/ 204 Pb, and 0.2% (2) for 208 Pb/ 204 Pb, and detailed analytical conditions and procedures can be seen in [47].

Electron Probe Microanalyses (EPMA) of Sulfides.
Mineral compositions were determined at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan, with a JEOL JXA-8100 Electron Probe Microanalyzer equipped with four wavelength-dispersive spectrometers (WDS).The samples were firstly coated with a thin conductive carbon film prior to analysis.The precautions suggested by [48] were used to minimize the difference of carbon film thickness between samples and obtain a ca.20 nm approximately uniform coating.During the analysis, an accelerating voltage of 20 kV, a beam current of 20 nA, and a 2 m spot size were used to analyze minerals.Data were corrected online using a modified ZAF (atomic number, absorption, and fluorescence) correction procedure.The peak counting time was 10 s for S, Fe, Cd, Cu, Pb, and Zn and 20 s for Ge and Ga.The background counting time was one-half of the peak counting time on the high-and low-energy background positions.The following standards were used: germanium (Ge), galena (S), pyrite (Fe), gallium arsenide (Ga), cadmium (Cd), copper (Cu), galena (Pb), and zinc (Zn).The minimum hydrothermal fluid temperatures were calculated by using the correlation of the Fe/Zn ratios in the linear equation [49].

Type and Occurrence of Fluid Inclusions.
Petrographic study of fluid inclusion was conducted with optical microscope.All samples contain abundant primary and pseudosecondary inclusions with ovoid, elliptical, irregular, and negative crystal shapes.In general, the primary FIs were identified by their occurrence in crystal growth zones of the host mineral or as the isolated inclusions (Figure 7(a)), while secondary FIs present as penetrating crystal boundaries (Figures 7(a) and 7(b)).These FIs vary from 4 m to over 60 m; most of them are 5 m to 30 m in size.On the basis of phase proportion at room temperature (25 ∘ C) and compositional components, three types of fluid inclusions were identified from quartz, rhodochrosite, fluorite, and sphalerite, as follows: (1) W-type: primary aqueous (W-type) inclusions consist of two phases (liquid phase and vapour phase) (Figures 7(c)-7(e)).They are the most common type, which occupy 60% of the total primary FIs in different stages and minerals (Figures 7(c)-7(e)).They vary from negative crystal to irregular shapes.The diameter of them range from 4 to 20 m, and their vapour volume spans from 5 to 45 percent.

Microthermometry.
All the C-type and W-type consistently homogenize into liquid phase.The microthermometric results, calculated salinities and densities of W-type and C-type of the early, main, and late stages, are shown in Table 1 and Figures 9 and 10.When these S-type FIs were   heated, the vapour phase disappeared with the temperature of 200-260 ∘ C; however these daughter minerals remained to be stable in size with the temperature up to 450 ∘ C. Considering that the S-type FIs were not completely homogenized, the microthermometric data of S-type FIs are not included in Table 1.
In the early stage, the W-and C-type FIs are recognized as primary FIs in the quartz grain.The melting temperatures ( m-CO2 ) of solid CO 2 in the C-type FIs range from −57.4 to −57.1 ∘ C, which are similar to the triple point for pure CO 2 (−56.6 ∘ C).Clathrate melting temperatures ( m-cla ) vary from 5.5 to 6.8 ∘ C, with which the salinities are calculated to be 6.1 to 8.3 wt.% NaCl equivalent by the equation of Hall et al. [40].Homogenization temperatures ( h-CO2 ) of the carbonic phase to the liquid phase are between 25.6 and 29.2 ∘ C. The total homogenization temperatures ( h ) of Ctype FIs to the liquid phase range from 294 to 298 ∘ C. The CO 2 densities are estimated to be from 0.65 to 0.70 g/cm 3 , and their bulk densities range from 0.98 to 1.00 g/cm 3 .Wtype FIs yield the ice-melting temperatures ( m-ice ) from −7.5 to −2.4 ∘ C, and their salinities are estimated to be 4.0 to 9.9 wt.% NaCl equivalent by the equation of Roedder [41].These FIs completely homogenize to the liquid phase at temperatures ( h ) of 244-334 ∘ C. Their densities are also acquired, fluctuating from 0.74 to 0.87 g/cm 3 (Table 1).
In the main stage, three transparent minerals are measured for the microthermometry, including rhodochrosite, fluorite, and sphalerite.In rhodochrosite crystals, ice-melting temperatures ( m-ice ) of W-type FIs range from −8.9 to −1.9 ∘ C, with corresponding salinities from 3.2 to 12.7 wt.% NaCl equivalent.They are homogenized into the liquid phase at the temperatures of 188 to 298 ∘ C, and their densities are 0.86 to 0.96 g/cm 3 .
Fluorite crystals from the main stage contain abundant C-and W-type FIs.The C-type FIs have narrow  m-CO2 range from −57.4 to −56.4 ∘ C. The CO 2 clathrate was melted at the temperatures of 7.1 to 7.9 ∘ C, corresponding to salinities of 4.3 to 5.6 wt.% NaCl equivalent.The homogenization temperatures ( h-CO2 ) range from 29.5 to 31.1 ∘ C. The CO 2 densities and total densities are estimated to be 0.47-0.62g/cm 3 and 0.92-0.99g/cm 3 , respectively.The C-type FIs show a total homogenization (all homogenized to liquid) at temperatures from 197 to 268 ∘ C.These W-type FIs yield final ice-melting temperatures ( m-ice ) of −5.3 to −0.3 ∘ C and the salinities of 0.5 to 8.3 wt.% NaCl equivalent.They homogenized to the liquid phase at the temperatures of 194 to 278 ∘ C.
FIs in sphalerite are mainly composed of C-and Wtype FIs.The solid CO 2 yields the melting temperatures ( m-CO2 ) of −57.3 to −55.8 ∘ C. The  m-cla are between 6.2 and 8.9 ∘ C, which correspond to the salinities from 2.2 to 3.9 wt.% NaCl equivalent.The homogenization temperatures ( h-CO2 ) vary from 20.7 to 27.7 ∘ C, and the total homogenization temperatures ( h ) of these C-type FIs into the liquid phase span from 190 to 224 ∘ C. Their CO 2 densities are estimated to be 0.66-0.77g/cm 3 , and their bulk densities are from 0.90 to 0.96 g/cm 3 . m-ice of the W-type range from −5.9 to −1.3 ∘ C, corresponding to salinities of 2.2-9.1 wt.% NaCl equivalent.The W-type FIs totally homogenize to the liquid phase at 164-245 ∘ C, with the calculated densities of 0.86-0.96g/cm 3 .
As a result of the extremely small size of quartz and carbonate grains, unfortunately, none of FIs has been measured in quartz-calcite-dolomite veins.Lu et al. [50] proposed that the secondary fluid inclusions could provide clues to study the nature of fluids after mineral crystallization.Furthermore, this view was also supported by Peng et al. [39], who reported that the characteristics of secondary FIs at the main stage (quartz and calcite veins) are similar to these of primary FIs at the late stage (calcite veins) in the Hongniu-Hongshan Cu skarn deposit.So we chose secondary W-type FIs in the fluorite and sphalerite in the main stage try to represent the characteristics of the late stage in the Zhijiadi deposit.Many secondary W-type FIs observed in the grains of fluorite and sphalerite have been measured in this study.All these kinds of FIs homogenize into the liquid phase.These W-type FIs in fluorite completely homogenize at the temperatures of 115-174 ∘ C. Their  m-ice vary from −4.6 to −0.2 ∘ C with the calculated salinities of 0.2-7.3wt.% NaCl equivalent.Similar to them, those FIs in the sphalerite have the icemelting temperatures ( m-ice ) of −5.7 to −1.0 ∘ C, and the salinities are estimated to be 1.7-8.8wt.% NaCl equivalent.The homogenization temperatures ( h ) range from 111 to 152 ∘ C. the Zhijiadi deposit are given in Table 2.The  13 C V-PDB values vary from −2.8‰ to −4.4‰, with an average of −3.0‰.The  18 O V-PDB values range from −14.4‰ to −17.5‰, with an average of −16.6‰.The  18 O V-SMOW values are calculated to be 12.8‰ to 16.03‰, and their average value is −13.7‰.The C and O isotopic compositions are relatively homogenous in the rhodochrosite.

Sulfur Isotopic Compositions.
One hundred and eighteen in situ LA-MC-ICP-MS sulfur isotopic compositions were determined from one quartz porphyry sample (3 in pyrite) and thirteen ore samples (44 in pyrite, 24 in sphalerite, 44 in galena, and 3 in chalcopyrite) covering the early and main stages.The results are shown in Table 3 and Figure 11.
As a whole, the  34 S values of sulfides in the Zhijiadi deposit show a decreasing sequence of  34 S Py (2.8‰) >  34 S Sp (1.8‰) >  34 S Ccp (1.7‰) >  34 S Gn (−1.8‰), which is consistent with the  34 S enrichment at the equilibrium conditions of sulfur isotope fractionation [51].Hence, the sulfur isotopic equilibrium temperatures of mineral pair for the sulfides (sphalerite and galena) and quartz porphyry (subvolcanic rocks) from the Zhijiadi deposit.The base map is from [34].These ranges of subvolcanic rocks (quartz porphyry and granite porphyry) are quoted from the data of [35].
Table 2: Carbon and oxygen isotope data of rhodochrosite from the Zhijiadi Ag-Pb-Zn deposit.(sphalerite and galena) can be calculated using the fractionation equation [51].The equilibrium temperatures for mineral pair of sphalerite and galena are estimated to be 237 to 267 ∘ C.  [35].In addition, all lead isotope data of sulfides are plotted into the area of subvolcanic rocks (Figure 12).

EPMA Data of Sphalerite.
A total of seven electron microprobe spot analyses were performed in sphalerite.The EPMA data for the main stage are summarized in Table 5.
Among all the analytical elements, Ga, Cd, and Pb are mostly below the detection limit.

Discussion
6.1.Nature and Evolution of Ore-Forming Fluids.Carbon dioxide (CO 2 ) is a common component of ore-forming fluids in orogenic-type and granitic intrusion-related ore deposits [52][53][54].It is widely accepted that CO 2 plays a critical role during metal transport by buffering the fluid in a pH range where elevated metal concentration can be maintained by complexation with reduced sulfur [55,56].On the other hand, the abundance of C-type FIs and the content of CO 2 in ore-forming fluids give us an effective way to understand the nature of the ore-forming fluids and the genesis of some deposit; for example, Chi and Xue [25] confirmed the sedimentary basin-hosted Jinman Cu deposit in Yunnan was formed in a hydrothermal system rather than basinal mineralization system on the basis of the existence of a large number of CO 2 -rich FIs in this deposit; Lu et al. [50] proposed the content of CO 2 in ore-forming fluids can be a marker for gold exploration; Chen et al. [24] also regard the content of CO 2 as an important criterion for the division of different types of hydrothermal gold deposits.There are a large number of CO 2 -bearing (C-type) FIs having been observed and measured in the Zhijiadi Ag-Pb-Zn deposit, so it is possible to systemically and comprehensively constrain the nature of initial ore-forming fluids.C-type FIs are distributed in the early quartz-pyrite veins.These C-type FIs are composed of two-phase and minor three-phase FIs and generally coexist with the W-type.Laser Raman spectroscopy shows that the gas phase of C-type FIs contains CO 2 .Their  m-CO2 temperatures range from −57.7 to −57.1 ∘ C, which are similar to the triple point for pure CO 2 (−56.6 ∘ C), showing that the carbonic phase in these FIs is nearly pure CO 2 .The peak homogenization temperature of these FIs (C-type and W-type FIs) in the early stage ranges from 280 to 300 ∘ C (Figure 10(a)), and their salinities range from 4.0 to 9.9 wt.% NaCl equivalent (Figure 10(b)).Hence, the initial ore-forming fluids are characterized by mediumlow temperature and low salinity and belong to a NaCl-H 2 O-CO 2 system.
Compared with the early stage, the main stage minerals mainly contain abundant C-type, W-type, and S-type FIs, especially fluorite and sphalerite crystals.These transparent daughter minerals in S-type FIs for fluorite and sphalerite are calcite and dolomite, determined by laser Raman spectroscopy, while a few opaque daughter minerals are detected as chalcopyrite in S-type FIs of the sphalerite.Although only W-type FIs are observed in rhodochrosite, laser Raman spectroscopy results have identified the peak of CO 2 in the vapour bubble of some W-type FIs.The low content of CO 2 in FIs of rhodochrosite crystals could be the product of the reaction: . In fluorite and sphalerite grains, these C-type FIs generally coexist with these transparent daughter-bearing S-type FIs.The CO 2 melting temperatures ( m-CO2 ) range from −57.3 to −55.8 ∘ C, closer to the triple point for pure CO 2 (−56.6 ∘ C).The salinities of C-type FIs in the main stage are estimated to be 0.2-5.6 wt.% NaCl equivalent, which are lower than that of C-type FIs in the early stage (6.1-8.3 wt.% NaCl equivalent).As a whole, the ore-forming fluids at this stage are characterized by relatively lower temperature (peak temperature between 200 and 240 ∘ C) (Figure 10(c)) and salinity (0.5-12.7 wt.% NaCl equivalent) (Figure 10(d)) than the early stage.Additionally, a lot of fluorite crystallizing at this stage illustrates that the ore-forming fluids are enriched in F − and belong to an F-rich NaCl-H 2 O-CO 2 system.
In the late stage, as a result of the extremely small size of quartz and carbonate grains, none of FIs has been observed in quartz-carbonate veins.However, abundant secondary Wtype FIs in fluorite and sphalerite have been measured with the features of low temperature (111-174 ∘ C) (Figure 10(e)) and low salinity (0.2-8.8 wt.% NaCl equivalent) (Figure 10(f)).To a certain extent, petrographic characteristics and compositional types of FIs can be considered to stand for the nature of ore-forming fluids at the late stage [50].Hence, the ore-forming fluids in the late stage are likely to be characterized by low temperature (Figure 10(e)) and low salinity (Figure 10(f)).
In summary, the ore-forming fluids in the Zhijiadi deposit are characterized by medium-low temperature, low salinity, and an F-rich NaCl-H 2 O-CO 2 system.

Mineralization Trapping
Pressure.The estimation of trapping pressure and mineralization depth are of great significance for targeting of deep-seated blind orebodies and distinguishing deposit type [57].The existence of abundant C-type and W-type FIs in the Zhijiadi deposit makes it possible to obtain the trapping pressure and mineralization depth [58,59].Isochores of fluid inclusions are defined on the P-T diagram through the Flincor program with the equation proposed by Bowers and Helgeson [42].
In the early stage, C-type FIs from the quartz-pyrite veins were selected for the pressure estimation.The maximum and minimum densities of these C-type FIs are 0.70 and 0.65 g/cm 3 , respectively, and they have been completely homogenized into the liquid phase at the temperature of 188-298 ∘ C, with the peak temperature of 280-300 ∘ C. Their peak homogenization temperatures are used to represent the lower limit of fluid temperatures, and these isochores of the C-type FIs imply that the pressure of ore-forming fluids in this stage is approximately 203 to 299 MPa at the temperature of 280-300 ∘ C (Figure 13(a)).
Abundant W-type and C-type FIs in the fluorite provide a great opportunity to acquire reliable trapping pressure of ore-forming fluids in the main stage.The intersection of isochores in different types of inclusions (i.e., W-and C-type) is one important method for pressure estimation [41,50,59].The maximum and minimum densities of Ctype FIs in fluorite are 0.62 g/cm 3 and 0.47 g/cm 3 , and those of W-type FIs are 0.91 g/cm 3 and 0.82 g/cm 3 , respectively.The trapping pressure based on intersections among the isochores in different types of FIs is from 32 to 158 MPa at the temperature of 210-279 ∘ C (Figure 13(b)).This intersection temperature is consistent with the independent temperature calculated from EPMA data of sphalerite (236-275 ∘ C) and/or the sulfur isotope equilibrium temperatures of sphalerite and galena pairs (237-267 ∘ C) within the analytic error.Hence, the initial trapping pressure of ore-forming fluids in the main stage can be obtained to be 32 to 158 MPa at temperature of 210-279 ∘ C. As is shown in Figure 13(c), there is an obviously downward trend where the trapping pressures drop from 203-299 MPa in the early stage to 32 to 158 MPa in the main stage.
On the basis of the deposit geology and regional geology, the pressure estimated from this study can stand for the pressure evolution of initial ore-forming fluids from the early to main stage by the following reasons: (1) the occurrence of most of orebodies in the Zhijiadi deposit is controlled by the fault F 2 .The F 2 fault belongs to a compression-shear type fault (Figure 4 process in the F 2 fault zone was carried out under a relatively sealing condition and high pressure; (2) it is important to note that the early stage quartz veins are typically quite narrow, and the central part of the early stage quartz veins is refilled by the main stage polymetallic sulfides (mainly galena and sphalerite) (Figure 6(a)), which supports the point where there is a tectonic shift from the close system in early stage (203-299 MPa) to the relatively open system in the main stage (32 to 158 MPa); (3) geologic observations show that hydrothermal breccia, cataclastic rocks, and fault gouge are well developed along the F 2 fault zone.This may indicate the existence of multiple deformation events under different structural regimes [60]; (4) in general, the orogenic-type deposits are characterized by the high-content of CO 2 (e.g., [24]) and have an obvious decrease from the early high pressure to the late low pressure [61,62].The ore-forming fluids in the Zhijiadi deposit have also characteristics by the existence of abundant CO 2 ; (5) Zhang et al. [17] proposed that the metallogenic timing of the Zhijiadi deposit is considered to be the Late Jurassic to Early Cretaceous.During this period, there existed an obviously tectonic transformation in the Wutai-Hengshan region from the compression to the extension [32], which corresponds to the pressure evolution of ore-forming fluids in the Zhijiadi deposit from the early stage to the main stage.We suggest that such circumstances may be caused by tectonic activity, which means the metallogenic pressure (203-299 MPa in the early stage and 32 to 158 MPa in the main stage) of the Zhijiadi deposit may record the tectonic activities rather than the upper crustal rock conditions.In general, hydrostatic pressure and lithostatic pressure are used to estimate the mineralization depth.Given that 2.7 g/cm 3 is the average density of upper crustal rocks, the pressure of the early and main stages corresponds to the lithostatic depth of 7.5-11.1 km and 1.2-8.1 km, respectively.However, these estimated depth results do not conform to the geological characteristics that abundant volcanic rocks (e.g., tuff) occurred in the mining area which may be produced at relatively shallow environment.The existence  [36,37]. of multistage breccia in the Zhijiadi deposit indicated that fluid pressure locally exceeded lithostatic pressure [63].In consideration of early stage quartz-pyrite vein crossed by the main stage galena + sphalerite vein (Figure 6(a)), it may imply that the metallogenetic environment changed from compression (shear) to extension (open) dynamics.As a result, the estimated depth, 7.5-11.1 km for early stage and 1.2-8.1 km for main stage, by using 2.7 g/cm 3 as the average density of upper crustal rocks, could not represent the true metallogenic depth.

Sources of the
Ore-Forming Fluids and Materials 6.3.1.Source of Ore-Forming Fluids.The C-O isotopic compositions of rhodochrosite from the Zhijiadi deposit suggest that they formed in a magmatic-hydrothermal system.All of their  13 C V-PDB values (−2.8 to −4.4‰) are significantly higher than organic material in sedimentary rocks (−30 to −10‰ [64]) but are similar to or slightly higher than those of mantle (−5 to −7‰ [65]) and igneous rocks and magmatic systems (−3 to −30‰ [65]).Moreover, the plotted data (Figure 14) indicates that CO 2 of ore-forming fluids might be main source from a magmatic-hydrothermal system, sharing similar source with the H 2 O in the Zhijiadi deposit [20].

Source of Ore-Forming Materials.
Previous studies reported that the sulfur isotopic compositions of sulfides in the Zhijiadi deposit were relatively stable around zero ( 34 S ≈ 0‰), demonstrating that sulfur of this deposit was provided by the deep-sourced magmatic system [14,19,20].However, the  34 S values of the sulfides in this study were found to vary in a wide range (−17.0‰ to +5.5‰) (Figure 11) by means of in situ sulfur isotope analyses.Compared with the  34 S values of local Early Cretaceous subvolcanic intrusions (quartz porphyry, 2.3‰ to 3.3‰) and magmatic source ( 34 S = 0 ± 3‰ [64,66]), the S isotopic compositions for a majority of sulfides in the Zhijiadi deposit seem to be similar to them despite the slight difference caused by galena at the later main stage.Particularly, the  34 S values of pyrite (+0.2‰ to +5.5‰) in the early stage, sphalerite (0‰ to +3.3‰), and chalcopyrite (+1.1‰ to +3.0‰) in the main stage yield similar range to magmatic source, whereas the  34 S values of galena in the latest crystallization of the main stage are the most negative with galena sulfur isotope data fluctuating in a large span from −17.0‰ to 2.9‰.
Sakai [67] and Ohmoto [38] proposed that the S isotopic compositions of sulfur minerals that formed in the hydrothermal system are influenced by not only the  34 S value of source materials, but also the physicochemical conditions (temperature, pH value, oxygen fugacity, and ion activity) of the ore-forming fluids.In general, there are two main probably factors resulting in the negative sulfur isotopic value: (1) biologic processes (e.g., reduction of sulfate to sulfide by sulfur-reducing bacteria during the late stage mineralization and the mixture of biogenic and magmatic sulfur lead to the isotopic values as low as −17.0‰); (2) the sharp change of oxygen fugacity from the main stage to the late stage [38].Firstly, there are not any sulfide-bearing sedimentary rocks or biogenetic texture of these sulfides.Secondly, the C-O isotopic compositions of rhodochrosite also suggest that there is no addition of the sedimentary organic matter into the ore-forming fluids.Therefore, the former explanation for the negative sulfur isotopic values of these sulfides can be excluded.Moreover, previous H-O isotopic analyses of the Zhijiadi deposit show that its ore-forming fluids mainly came from magmatic system with a significant input of oxidized meteoric water [20].Hence, the large variation and negative values in the sulfur isotopic compositions of galena can be interpreted as the shift in oxidation state, which was caused by progressive mixture with oxidized meteoric waters.In general, it can be concluded that the sulfur of the Zhijiadi deposit can be attributed to a magmatic sulfur reservoir with variation caused by the change in oxygen fugacity between these ore stages (Figure 15) [38].
The Pb isotope compositions of the sulfides from the Zhijiadi deposit show a narrow range, with all data plotted into the area of the Early Cretaceous subvolcanic rocks.The ratios of 206   12(a)), while the lead isotope data are plotted into the field between the lower crust and orogenic belt evolution lines in the 208 Pb/ 204 Pb versus 206 Pb/ 204 Pb diagram (Figure 12(b)).This suggests that the sulfides and subvolcanic rocks share the same Pb source, namely, the lower crust, further demonstrating that the ore-forming materials were of magmatic origin.systems (e.g., [68]).In combination with geological characteristics and stable isotope studies, fluid inclusion study can shed light on the evolution of ore-forming fluids and provide useful constraints on the metallogenic process.
In natural systems, metals in an ore-forming fluid are generally found in complexes with aqueous Cl − and HS − and are transported in the form of metal-ligand complexes (e.g., [69]).Temperature is an important factor for the stability of metal complexes and the solubility of ore-forming materials.A certain degree of cooling for the ore-forming fluids could not only destroy the stability of the metal complexes but also reduce the solubility of the ore-forming materials, and it finally causes the precipitation of base metals from the hydrothermal system [70][71][72].The fluid evolution diagram for the Zhijiadi deposit, represented by a homogenization temperature versus salinity plot, is shown in Figure 16.There is a clear trend that the homogenization temperature and salinity of the hydrothermal fluids gradually decreased from the early stage, through the main stage, to the late stage.
The trapping pressures that we obtained also share similar variation, dropping from 203-299 MPa to 32-158 MPa during the early and main stages of the hydrothermal process.Therefore, remarkable reduction of pressure could lead to the ascension of the magmatic fluid along fault zones to higher levels, and the magmatic-hydrothermal fluids then were mixed with circulating meteoric water, finally precipitating sulfides and forming ore deposit (e.g., [63]).Hence, the reduction of pressure has a significant effect on the reduction of temperature and mixing and/or dilution for ore-forming fluids, as well as the formation of ore deposit.Furthermore, previous H-O isotopic analyses show that the ore-forming fluids are sourced from a magmatic system with a significant input of meteoric water into it [20].Progressive input of meteoric water caused the gradual dilution of ore-forming fluid and reduced its temperature and salinity (Figure 15).
A great number of researches were carried on the Zhijiadi deposit, whereas its genesis remains controversial.On the basis of geological characteristics, stable isotopic analysis, and the mineralization temperature, there are two main views that were proposed to interpret its genesis, for example, (1) subvolcanic rock-hosted medium-low temperature hydrothermal deposit [14] and (2) epithermal deposit [17,73].However, our new data from the estimation of trapping pressure from fluid inclusions and stable isotope analyses indicated that the Zhijiadi deposit has many differences from typical epithermal deposits.
Up to now, it has been widely accepted that the epithermal deposit generally forms at low temperatures (<300 ∘ C), medium-low salinity fluids, and low content of CO 2 , with relatively low mineralization pressure [24,[74][75][76].However, the microthermometric data in our work suggests that the ore-forming fluids in the Zhijiadi deposit are characterized by medium-low temperature (117-334 ∘ C), low salinity, and abundant CO 2 .The trapping pressures range from 203-209 MPa in the early stage to 32-158 MPa in the main stage, which are obviously different from the characteristic of epithermal deposits (<10 MPa) (e.g., [76]).
In addition, the C-O-S isotope data indicate that the oreforming fluids and materials are sourced from a magmatic system.The Pb isotopes further confirm the sulfides and subvolcanic rocks have the same sources, namely, the lower crust, demonstrating that the ore-forming materials were of magmatic origin.Furthermore, the occurrence of orebodies in the Zhijiadi deposits are strictly controlled by a series of NW-, EW-, and NE-trending faults, especially for the F 2 fault.The hydrothermal alterations associated with Ag-Pb-Zn mineralization in the Zhijiadi deposit are mainly composed of silicification, carbonation, sericitization, and chloritization.The gangue minerals are quartz, rhodochrosite, fluorite, calcite, and dolomite, which are not accordant with the typical characteristics of alternation in epithermal deposits (e.g., adularization or alunitization [24,76,77]), which can further rule out the possibility that the Zhijiadi deposit is a kind of epithermal deposit.Based on all available geological and geochemical evidence, we propose that the Zhijiadi deposit belongs to medium-low temperature hydrothermal deposit related to volcanic and subvolcanic magmatism strictly controlled by fault zones instead of epithermal deposit.

Conclusions
(1) The Zhijiadi deposit, located in the central North China Craton, is hosted in volcanic and subvolcanic rocks and fault zones.According to the mineral assemblages and crosscutting relationship, hydrothermal metallogenic processes can be divided into the early stage, the main stage, and the late stage, represented by quartz + pyrite veins, rhodochrosite + fluorite + galena + sphalerite + chalcopyrite + multiple Agbearing sulfides, and quartz + carbonate veins.
(2) The fluid inclusions of the Zhijiadi deposit are primarily aqueous (W-type), aqueous-carbonate (C-type), and daughter mineral-bearing (S-type) inclusions.Marcasite is recognized as the daughter minerals in the early stage quartzpyrite vein, whereas chalcopyrite, calcite, and dolomite occur as daughter minerals in the main stage fluorite and sphalerite.The ore-forming fluids are characterized by medium-low temperature, low salinity, and abundant CO 2 and belong to an F-rich H 2 O-NaCl-CO 2 system.From the early, through the main, to the late stage, the homogenization temperatures and salinities of the hydrothermal fluids gradually decrease.
(3) The C-O isotopic compositions indicate that the oreforming fluids of the Zhijiadi deposit were sourced from a magmatic system in origin.Analyses of S isotopes suggest that the sulfur of Zhijiadi deposit was attributed to the magmatic sulfur reservoir with variation caused by the change in oxygen fugacity.The Pb isotopes further confirm the sulfides and subvolcanic rocks share the same sources, namely, the lower crust.
(4) The Zhijiadi deposit belongs to medium-low temperature hydrothermal deposit strictly controlled by the fault zones.There is an obviously downward trend that the trapping pressures of ore-forming fluids in the Zhijiadi deposit drop from 203-299 MPa in the early stage to 32 to 158 MPa in the main stage.The reduction of pressure has a significant effect on the reduction of temperature, mixing, and/or dilution for ore-forming fluids, as well as the formation of ore deposit.

Figure 2 :
Figure 2: (a) Simplified geologic map (modified from unpublished internal map of the third geological exploration institute of China Metallurgical Geology Bureau).(b)The distribution of orebodies in the Zhijiadi deposit (modified from[14]).

4. 1 .
Fluid Inclusion Studies.Samples were collected from the orebodies in level 1320 around F 2 fault zone in the Zhijiadi deposit.Thirty-seven double-polished thin sections (about 200 m thick) associated with three hydrothermal

4. 3 .
Sulfur Isotope Analysis.Thirteen standard polished sections (one pyrite-bearing quartz porphyry sample and 12 ore samples) are collected to carry out in situ sulfur isotope 20.0 kV ×40 100 m WD 11 mm

Figure 6 :
Figure 6: Photomicrographs showing three stages in the Zhijiadi deposit.(a) The early stage quartz-pyrite vein.(b) The main stage rhodochrosite polymetallic sulfides vein cross-cutting quartz porphyry and some sphalerite and galena disseminated in quartz porphyry.(c) The main stage fluorite-galena filling in the dacite breccia.(d) Coexistence of sphalerite and rhodochrosite in the hydrothermal breccia.(e) Brecciated sphalerite distributing in the altered quartz porphyry breccia.(f) The late stage quartz-carbonate vein cross-cutting the tuffaceous breccia.Py, pyrite; Gn, galena; Sp, sphalerite; Rds, rhodochrosite; Fl, fluorite; Cal, calcite; Dol, dolomite.

5. 2 .
Laser Raman Spectroscopy.Laser Raman spectroscopy analyses were carried out to constrain the fluid inclusion compositions of representative samples from three different hydrothermal stages, and the results are shown in Figure8.The results show that the vapour phase of aqueous-carbonic inclusions in the early stage quartz grains was confirmed as CO 2 by the peak of 1281 and 1385 cm −1 (Figure8(a)).The gas phase of two-phase fluid inclusions in the main stage rhodochrosite contains minor CO 2 with the peak of 1283 and 1386 cm −1 (Figure8(b)).The main stage fluorite and sphalerite also contain abundant CO 2 in the vapour phase of C-type FIs, as indicated by the peak of 1280 and 1384 cm −1 and 1282 and 1384 cm −1 , respectively (Figures8(c) and 8(d)).Transparent daughter minerals in S-type FIs of the main stage fluorite and sphalerite are calcite and/or dolomite at the peaks of 1081 cm −1 and 1098 cm −1 , respectively (Figures8(e) and 8(f)).Opaque daughter minerals in the early stage quartz are marcasite, which are indicated by the peak of 339 cm −1 (Figure8(g)).The 289 cm −1 peak demonstrates that some opaque daughter minerals in the main stage sphalerite are chalcopyrite (Figure8(h)).
S e c o n d a r y F ls Sp

Figure 7 :Figure 8 :Figure 9 :
Figure 7: Photomicrographs showing the types and distribution of fluid inclusions observed in the Zhijiadi deposit.(a) Primary aqueous (W-type) inclusions coexisting with daughter mineral-bearing (S-type) inclusions in the main stage fluorite and these secondary W-type inclusions in a linear distribution penetrating crystal boundaries of fluorite.(b) Cluster of primary W-type inclusions in the main stage sphalerite and these secondary W-type inclusions penetrating crystal boundaries of sphalerite.(c) W-type inclusions clustered in the early stage quartz.(d) Primary W-type inclusions clustered in the main stage rhodochrosite.(e) Isolation of W-type inclusions from the main stage sphalerite.(f) Coexistence of aqueous-carbonic (C-type) inclusions and W-type inclusions in the early stage quartz crystals; C-type inclusions in the main stage fluorite (g) and sphalerite (h); (i) S-type inclusions in main stage fluorite (daughter minerals including calcite and dolomite), showing an obvious difference from the secondary W-type inclusions; (j) calcite-bearing S-type inclusions with CO 2 vapour bubble in the main stage fluorite; (k) marcasite-bearing S-type inclusions in the early stage quartz; (l) chalcopyrite-bearing S-type inclusions in the main stage sphalerite.V CO2 , vapour CO 2 ; L CO2 , liquid CO 2 ; V H2O , vapour H 2 O; L H2O , liquid H 2 O; Qz, quartz; Rds, rhodochrosite; Fl, fluorite; Sp, sphalerite; Cal, calcite; Dol, dolomite; Mrc, marcasite.

GeofluidsFigure 10 :
Figure 10: Histograms of homogenization temperatures ( h ) and salinities of C-type and W-type fluid inclusions in the early stage (a, b), the main stage (c, d), and secondary W-type fluid inclusions from fluorite and sphalerite which represent the late stage (e, f).Qz, quartz; Rds, rhodochrosite; Sp, sphalerite; Fl, fluorite.

Figure 13 :
Figure 13: Pressure-temperature conditions for the early and main stages in the Zhijiadi deposit.(a) Minimum trapping pressure for the early stage based on isochores of C-type FIs in the quartz and the peak homogenization temperatures.(b) Trapping pressures for the main stage based on intersections of isochores of the minimum and maximum densities from the C-and W-type FIs in fluorite.(c) Model for pressure-temperature evolution of the ore-forming fluids from the early to main stages.Qz, quartz; Sp, sphalerite; Fl, fluorite.

Figure 14 :
Figure 14:  18 O versus  13 C diagram of rhodochrosite from the Zhijiadi deposit.The base map is from[36,37].

6. 4 .Figure 15 :
Figure15: log  O2 versus pH isothermal diagram showing stability relationships in the system Fe-O-S for early and main stage (A) and late stage (B) of the paragenesis (modified from[38]).Given Fe-O-S mineral boundaries at ΣS = 0.1 moles/kg H 2 O,  = 250 ∘ C, and  = 1.0 condition, thinner red lines represent the stability field boundaries for pyrite, pyrrhotite, magnetite, and hematite with the following reasons: (1) the sulfur isotopic equilibrium temperatures are estimated to be 237-267 ∘ C; (2) previous study shows that pH value of the ore-forming fluid is near-neutral pH[19]; (3) the sericitization occurs in the Zhijiadi deposit.Kaol, kaolinite; Ms, muscovite; K-feld, K-feldspar.

Figure 16 :
Figure 16: Homogenization temperature versus salinity plot of fluid inclusions showing fluid evolution in the Zhijiadi deposit.

Table 1 :
Microthermometric data of fluid inclusions from different stages in the Zhijiadi Ag-Pb-Zn deposit. m,CO2 , melting temperature of solid CO 2 ;  m,cla , temperature of CO 2 -clathrate dissociation;  h,CO2 , homogenization temperature of CO 2 ;  m,ice , ice-melting temperature;  h , total homogenization temperature.All of fluid inclusions including W-type and Ctype FIs homogenize into the liquid phase.

Table 3 :
Sulfur isotope data of sulfide minerals from the Zhijiadi Ag-Pb-Zn deposit.

Table 4 :
Lead isotope data of sulfides and rocks from the Zhijiadi Ag-Pb-Zn deposit.
Pb/ 204 Pb, 207 Pb/ 204 Pb, and 208 Pb/ 204 Pb for these sulfides vary from 16.5323 to 16.6984, 15.2634 to 15.2712, and 36.6673 to 36.9998, respectively.In the 207 Pb/ 204 Pb versus 206 Pb/ 204 Pb diagram, all of the samples are plotted near to the lower crust line (Figure