Abstract

The Xiejiagou deposit is a representative medium-sized gold deposit in Jiaodong the Peninsula, which contains gold reserves of 37.5 t. The orebodies are hosted in the Linglong biotite granite with a zircon LA-ICP-MS U–Pb age of (, ) and are characterized by disseminated- or stockwork-style ores. Mineralization and alteration are structurally controlled by the NE-striking fault. Three stages of mineralization were identified with the early stage being represented by (K-feldspar) sericite quartz pyrite, the middle stage by quartz gold polymetallic sulfide, and the late stage by quartz carbonate. Ore minerals and gold mainly occurred in the middle stage. Three types of primary fluid inclusions were distinguished in the Xiejiagou deposit, including carbonic-aqueous, pure carbonic, and aqueous inclusions. The primary fluid inclusions of the three stages were mainly homogenized at temperatures of 262–386°C, 192–347°C, and 137–231°C, with salinities of 2.22–8.82, 1.02–11.60, and 1.22–7.72 wt% NaCl equivalent, respectively. These data indicate that the initial ore-forming fluids were a medium temperature, CO2-rich, and low-salinity H2O–CO2–NaCl homogeneous system, and the ore-forming system evolved from a CO2-rich mesothermal fluid into a CO2-poor fluid. Considering the fluid inclusion characteristics, H–O–S–Pb isotopes, and regional geological events, the ore-forming fluid reservoir was likely metamorphic in origin. Trapping pressures of the first two hydrothermal stages estimated from the carbonic aqueous inclusion assemblages were ~224–302 MPa and ~191–258 MPa, respectively. This suggests that the gold mineralization of the Xiejiagou gold deposit occurred at a lithostatic depth of ~7.2–9.7 km. Au(HS)2 was the most probable gold-transporting complex at the Xiejiagou deposit. Precipitation of gold was caused by a CO2 effervescence of initial auriferous fluids.

1. Introduction

The Jiaodong Peninsula (Figure 1) is an important gold metallogenic belt in China hosting about 4000 t of gold reserves and also considered as one of the most famous gold mineralization belts in the world [111]. The majority of gold deposits in the region are temporally and spatially related to Mesozoic NE-trending faults, and the deposits have been classified as Jiaojia-type and Linglong-type mineralization, where the former refers to a combination of disseminated and stockwork/veinlet gold mineralization and the latter to auriferous quartz veins [1, 1215]. The Jiaodong gold belt is dominated by Jiaojia-type gold deposits, which make up more than 80% of its gold resources [16, 17], and extensive research has been under taken on them ([14, 15, 1822] and the references therein). Gold mineralization in the region is generally considered to have occurred during the Early Cretaceous, with a peak age at [23], although there are still controversies regarding the ore-forming age [20].

Gold deposits are generally related to tectonic, geological, and geochemical processes, in which the nature of the hydrothermal fluids is of fundamental importance [2426]. Several issues related to the nature and sources of the ore-forming fluids remain debated, and a number of hypotheses have been proposed. Initially, Zhang [27] speculated that the ore-forming fluids were meteoric in origin. More recently, Fan et al. [16, 17] suggested that the ore-forming fluids originated from magmatic water, and Zeng et al. [28] suggested that the gold mineralization resulted from magmatic-hydrothermal processes. In addition, because some of the orebodies within these deposits are accompanied by intermediate-mafic dikes, therefore, the ore-forming fluids have been inferred to be derived from the mantle [28], although a combination of crustal and mantle components has also been proposed [16, 17]. Alternatively, based upon ore geological characteristics and mineralization style, the Jiaodong gold mineralization has been treated as an “orogenic type” and it has been proposed that the ore-forming fluids were derived from metamorphic water [21]. Although some geologists classified these as orogenic gold deposits, the nearly 1.7-billion-year gap between gold mineralization and the latest regional metamorphism, together with the tectonic settings, has raised doubt on the orogenic gold deposit model [5, 9, 18, 23]. Therefore, it is difficult to apply the conventional classification schemes for common ore deposit types (such as epithermal gold, orogenic gold, porphyry gold,…) to define the formation of the Jiaodong gold mineralization.

The Xiejiagou deposit is a representative medium-sized Jiaojia-type gold deposit, which contains gold reserves of >35 t. Given that this deposit is an important producer of gold, some geochemical studies have been undertaken with the aim of characterizing the nature of the Au mineralization (e.g., [2934]). Nevertheless, the nature and evolution of the fluid associated with gold mineralization have not been well documented through petrographic, isotopic, and fluid inclusion studies, which limit our understanding of gold mineralization of the Xiejiagou deposit. Based on previous geological information on the Xiejiagou deposit, in particular, petrographic, microthermometric, and Raman spectroscopic analyses for inclusions contained within different stage veins were conducted to investigate the ore-forming P-T conditions and examine the mechanisms of gold transport and deposition. In addition, based on the geologic characteristics and isotopic compositions (H, O, S, and Pb) of the ores from the deposit, we try to trace the origin of ore-forming fluids. We hope that this study can further address the implications for the better understanding of the gold mineralization in the Jiaodong Peninsula.

2. Regional Geology

The Jiaodong Peninsula, located in east Shandong Province, occupies the southeastern margin of the North China Block (NCB) is and bounded by the Tanlu fault in the west (Figure 1). This region consists of two tectonic units, the Jiaobei terrane in the west and Sulu terrane in the east, bordered by the Wulian-Qingdao-Yantai fault (WQYF) (Figure 1). The Jiaobei terrane comprises the Jiaobei uplift in the north and the Jiaolai basin in the south (Figure 1). The Jiaobei uplift is the primary host of gold on the Jiaodong Peninsula with gold reserves of more than 3600 t ([19], 90% of the proven gold resources). The main lithological units in the Jiaobei uplift include the metamorphosed Precambrian basement, as well as widespread Mesozoic intrusive rocks [35], Figure 1. The Precambrian basement includes the Neoarchean Jiaodong, Paleoproterozoic Jingshan, and Fenzishan groups [36], Figure 1. The Jiaodong group is characterized by widely distributed late Archean tonalitic, trondhjemitic, and granodioritic gneisses, as well as supracrustal amphibolites and gneisses that yield U–Pb zircon ages of mainly between 2.9 and 2.5 Ga [37]. The Paleoproterozoic Jingshan and Fenzishan groups contain widespread schist, marble, and amphibolite, plus minor amounts of mafic granulite. These groups yield metamorphic zircon ages of about 1.8 Ga and underwent amphibolite- to granulite-facies metamorphism and deformation [38]. Mesozoic granitoids that intruded into the Jiaobei uplift have been traditionally divided into the Late Jurassic Linglong granite and the Early Cretaceous Guojialing granite, with both of them being significant hosts for gold mineralization ([3941], Figure 1). Mesozoic volcanic rocks were widespread in the Jiaolai basin (Figure 1). The volcanic sequences mainly include trachybasalt, trachyte, and rhyolite, forming at about 110–95 Ma [42]. The Sulu terrane is formed by the subduction of the Yangtze Block beneath the NCB in the Triassic [43] and is characterized by the presence of high and ultra-high-pressure (HP to UHP) metamorphic rocks as well as Late Triassic, Late Jurassic, and Early Cretaceous granitoids (Figure 1). Subsequently, these rocks were unconformably overlain by Cenozoic basalt and Quaternary sedimentary rocks (Figure 1).

Most of these gold deposits in the Jiaobei uplift are controlled by the NE- to NNE-trending fault systems ([39], Figure 1): from west to east, the Sanshandao fault (SSDF), Jiaojia fault (JJF), Zhaoping fault (ZPF), Muping-Jimo fault (MJF), and Qixia fault (QXF). These faults are subsidiary faults to the continental-scale Tanlu fault zone [44]. Among them, the ZPF is one of the most important tectonic structures controlling second-order ore-hosting faults in the Jiaobei uplift [1, 44, 45].

3. Ore Deposit Geology

The Xiejiagou gold deposit is located about 5 km southwest of Zhaoyuan city in the northwestern part of the Jiaobei uplift and near the middle part of the ZPF (Figure 1). Structurally, the mineralization and hydrothermal alteration are controlled by a series of NE-trending secondary faults (ductile shear structures and brittle faults) located in the hanging wall of the ZPF fault zone (Figure 2). The principal lithologic unit in the Xiejiagou deposit is a pluton of biotite granite, the so-called Linglong granite. Numerous intermediate-mafic dikes are widespread in the Xiejiagou gold deposit. Among them, the lamprophyre dikes (pre-ore, Figure 2(b)) are cross-cut by orebodies or faults and display intense hydrothermal alteration with a whole rock K–Ar age of 123.6 Ma [29], whereas the orebodies are cut by the NW-trending diorite porphyrite, gabbro diorite, and dolerite dikes (post-ore, Figures 2(a) and 2(b)) with whole rock K–Ar ages of 98.6 to 115.2 Ma [29], 97.2 to 98.8 Ma [29], and 100.5 to 107.6 Ma [29], respectively. These geochronological data indicate that the age of Xiejiagou gold mineralization occurred at ~123.6–115.2 Ma, which is identical to , the generally considered main mineralizing period of the major gold deposits within the Jiaodong Peninsula [46]. The formation of the Jiaodong gold mineralization including Xiejiagou deposit was closely related with subduction of the Paleo-Pacific Plate beneath the Eurasian continent accompanied by lithospheric delamination, asthenospheric upwelling, and intense craton destruction [49].

The Xiejiagou deposit currently has proven reserves of about 37.5 t of Au [33]. Most of the orebodies are lenticular or irregular lodes and are controlled by the NE-trending faults, which cut the granitoids (Figure 2). The ore-controlling faults are 200–1000 m long and 2–15 m in width with a strike of 20 to 35° and a dip of 60° to 85° (Figure 2). Four major orebodies have been identified in the ore district, and no. 3 is the largest and most representative (Figure 2), containing roughly 70% of the total proven gold reserves of the deposit [32]. The no. 3 orebody is 600 m in length with a gentle dip of 76° to 85° and contains an ore grade of 3.0 to 6.8 g/t Au [31]. We have focused on the no. 3 orebody, as it is the most productive for Au. The mineralization style of the Xiejiagou deposit appears associated with pyrite sericite silica-altered granites or fine pyrite veinlets, belonging to the disseminated- and stockwork-style gold mineralization. Extensive hydrothermal alteration also occurs in ore-controlling faults and is characterized by K-feldspar, quartz, sericite, pyrite, chlorite, and carbonate (Figures 3(a)3(f)). Major ore minerals in the Xiejiagou deposit include pyrite, galena, sphalerite, and chalcopyrite, along with minor pyrrhotite (Figures 3(g)3(l)). The gangue minerals include 80% quartz with amounts of K-feldspar, sericite, chlorite, calcite, and clay minerals (Figure 3). Au occurs mostly as native gold, followed by electrum. Native gold grains occur mainly as inclusions in pyrite crystals (Figure 3(k)) and less commonly in tiny fissures of pyrite. Electrum occurs mainly as inclusions in sphalerite or galena crystals (Figure 3(l)).

4. Hydrothermal Quartz Vein Sequences and Mineral Assemblages

Three hydrothermal stages and a supergene stage have been identified in the Xiejiagou deposit on the basis of mineralogical assemblage, textures, and cross-cutting relationships observed in hand specimens and thin sections (Figures 3 and 4). The hydrothermal ore-forming processes occurred during an early (K-feldspar) sericite quartz pyrite stage, a middle quartz gold polymetallic sulfide stage, and a late quartz carbonate stage (Figure 5). The gold mineralization occurred in the middle stage.

The early hydrothermal stage is characterized by (K-feldspar) sericite pyrite silica alteration granites (Figure 4(a)) or in ore-barren quartz veinlet places (Figure 4(b)). Minerals comprise milky quartz, sericite, K-feldspar, and sparsely distributed euhedral cube pyrite (Figures 3(a), 3(b), 3(g), 4(a), and 4(b)). The main mineral assemblage in this stage is milky subhedral/anhedral quartz (Qz1), sericite, K-feldspar, and sparsely distributed euhedral cube pyrite, with variables of albite and muscovite (Figures 3(a), 3(b), 3(g), 4(a), 4(b), and 5). Not any gold was found in our observation. The early-stage ores were locally cut or cemented by middle- or late-stage veinlets (Figures 4(c) and 4(f)).

The middle stage is the main gold-producing stage and is characterized by the quartz gold polymetallic sulfide veinlets (Figures 4(c) and 4(d)). During this stage, large amounts of sulfide minerals were precipitated including pyrite, galena, sphalerite, chalcopyrite, and minor pyrrhotite (Figures 3(h)3(j)). Quartz in the middle stage (Qz2) is commonly dark-gray (smoky) in color coexisting with sulfide minerals and minor calcite (Figures 4(c)4(e)). Gold mainly occurs as native gold and that coexists with other sulfide minerals (Figures 4(e) and 4(f)).

The late stage is characterized by quartz carbonate veinlets. This stage mainly contains calcite, quartz, and siderite, with trace amounts of pyrite (Figures 4(e) and 4(f)). No gold has been identified under reflected light in this stage. The paragenetic sequence of the Xiejiagou gold deposit is summarized in Figure 5.

5. Sampling and Analytical Methodology

5.1. Zircon U–Pb Dating

Zircons were separated by an unaltered Linglong biotite granite sample (C156-13) from the south of the Xiejiagou gold deposit (Figure 2(a)). The main minerals of the medium-grained biotite granite sample (C156-13) are plagioclase (30%), K-feldspar (30%), quartz (35%), and biotite (5%) (Figure 6(a)), with accessory titanite, zircon and apatite. Zircon sample preparation, cathodoluminescence (CL), and back-scattered electron imaging were completed at the Nanjing Hongchuang Exploration Technology Service Company Limited. U–Pb dating was performed on a laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) for constraining the age of the wall-rocks of the deposit.

LA-ICP-MS U–Pb isotope analyses were conducted using an Agilent 7500 mass spectrometer connected to a 193 nm ArF excimer laser ablation system at the Wuhan Sample Solution Analytical Technology Company Limited. The synthetic silicate glass standard reference material NIST SMR610 was used to calibrate the instrument, and the international reference standard zircon 91500 was used as the external standard for age calibration. During the analyses, high-purity He gas was used to transfer the ablated materials. The spot size was 32 mm. For detailed procedures and methods, as well as the instrument settings, see Yuan et al. [47]. The ICPMSDataCal [48, 49] and Isoplot 3.0 software [50] were used for data reduction. Common Pb corrections were undertaken using the approach of Andersen [51]. Fifteen zircon analyses were performed via LA-ICP-MS. Uncertainties on individual LA-ICP-MS analyses are quoted at the 1σ level, with results given in Table 1.

5.2. Fluid Inclusion Studies

Samples for the fluid inclusion study were collected from the main orebody (lode no. 3) of the Xiejiagou gold deposit to determine the nature of the fluids associated with mineralization. (Figure 4). Doubly polished thin sections (about 0.20 mm thick) were made from 45 quartz samples associated with different stages. Fluid inclusion petrography included careful observation of the shapes, characteristics of spatial distribution, genetic and composition types, and vapor/liquid ratios. Twenty-nine typical samples with abundant and representative fluid inclusions were selected for microthermometric measurements and laser Raman spectroscopy analyses. Various inclusions in quartz were selected for microthermometry based on the distribution, size, types, and textural relationship. The microthermometric study was carried out at the Institute of Mineral Resources, Chinese Academy of Geological Sciences, and Institute of Geology, Chinese Academy of Geological Science, using a Linkam THMS-600 heating-freezing stage with a temperature range of −195°C to +600°C. The equipment was calibrated with synthetic samples of fluid inclusions: pure water inclusions (0°C), pure CO2 inclusions (−56.6°C), and potassium bichromate (398°C). The estimated precision of the measurements was ±0.2°C for temperatures lower than 31°C, ±1°C for the interval of 31–300°C, and ±2°C for temperatures higher than 300°C. The salinities of NaCl–H2O inclusions were calculated using the final melting temperatures of ice [52]. The salinities of CO2-bearing fluid inclusions were calculated using the melting temperatures of clathrate [53]. Densities and pressure were calculated using FLINCOR software according to the microthermometry data [54, 55]. Representative fluid inclusion volatiles were analyzed using a Renishaw System 2000 Raman Microspectrometer at the Institute of Mineral Resources, Chinese Academy of Geological Sciences, with a blazing power of 20 MW and blazing wave length of 514.5 nm according to the method of Burke [56].

5.3. Oxygen Isotope and Hydrogen Isotope Analyses

Six quartz samples that were intergrown with ores during the different stages of mineralization were collected for detailed hydrogen and oxygen isotope analyses. All samples were handpicked and checked under a binocular microscope to achieve >95% purity. The isotopic compositions of oxygen and hydrogen were analyzed using a MAT-253 stable isotope ratio mass spectrometer in the analytical laboratory of the Beijing Research Institute of Uranium Geology. Oxygen was extracted for analysis using the BrF5 technique [57], and hydrogen from H2O was released from fluid inclusions [58]. The isotopic ratios are reported in standard δ notation (‰) relative to standard mean ocean water (SMOW) for oxygen and hydrogen. The analytical precision was better than ±0.2‰ for δ18O and ±2‰ for δD. The oxygen isotope ratios of water in equilibrium with quartz were calculated by the equation [59]. is the average homogenization temperature of fluid inclusions from different stages.

5.4. Sulfur Isotope Analysis

Sulfur isotope analyses were carried out from 50 to 100 mg of sulfide minerals including pyrite, sphalerite, and galena. The sulfide grains were mixed with cuprous oxide and crushed into 200 mesh powder. SO2 was produced through the reaction of sulfide and cuprous oxide at 980°C under a vacuum pressure of . The SO2 was then measured by MAT-252 mass spectrometer for sulfur isotope. The sulfur isotope analysis was carried out at the analytical laboratory of the Beijing Research Institute of Uranium Geology. All the analytical uncertainties were better than ±0.2‰.

5.5. Lead Isotope Analysis

For lead isotope ratios, approximately 10 to 50 mg of five sulfide samples was first leached in acetone to remove surface contamination and then washed by distilled water and dried at 60°C in the oven. Washed sulfides were dissolved in dilute mix solution of nitric acid and hydrofluoric acid. Following ion exchange chemistry, the lead in the solution was loaded onto rhenium filaments using a phosphoric acid silica gel emitter. The lead isotopic compositions were measured on MAT-261 thermal ionization mass spectrometer with the standard sample NBS 981. The lead isotope analysis was carried out at the analytical laboratory of the Beijing Research Institute of Uranium Geology. The analytical precision of Pb isotope is better than ±0.09‰.

6. Results

6.1. Zircon U–Pb Age of Biotite Granite

Zircons from the biotite granite (C156-13) are subhedral to euhedral, with a crystal length of 100 to 300 μm and aspect ratios of 2 : 1 to 4 : 1. Most of the grains display well-developed oscillatory zoning (Figure 6(b)), indicating a typical magmatic origin [60, 61]. The zircons contain variable moderate concentrations of U (1445 – 2493 ppm) and Th (125 – 971 ppm). The Th/U ratios range from 0.09 to 0.85. They have concordant 206Pb/238U ages of to , yielding a concordia age of (, ; Figure 6(c)) and a weighted mean 206Pb/238U age of (, ; Figure 6(d)). The weighted mean 206Pb/238U age with a low mean square of weighted deviation (MSWD) of 1.2 is interpreted as the crystallization age of the biotite granite.

6.2. Type and Occurrence of Fluid Inclusions

Fluid inclusions (FIs) occurring as isolated inclusions, randomly distributed groups, or in clusters were considered to be primary or pseudosecondary in origin, whereas those aligned along microfractures in transgranular trails were interpreted as secondary [62, 63]. Abundant primary, secondary, and pseudosecondary fluid inclusions are observed (Figure 7(a)).

Three compositional types of fluid inclusions (Figure 7) were identified based on their phases at room temperature, their observed phase transitions during heating and cooling runs, and laser Raman spectroscopy results. The FI types, in decreasing order of abundance, are aqueous carbonic (Figures 7(b)7(d); CO2–H2O–NaCl; type Ι), pure carbonic (Figure 7(e); CO2; type II), and aqueous (Figures 7(f) and 7(g); H2O–NaCl; type IΙΙ).

Type Ι inclusions are as follows: CO2-bearing aqueous solution inclusions comprise 60% of the total inclusions and are mainly three-phase CO2-bearing fluid inclusions () at room temperature (Figure 7(b)). Occasionally, two phases can be observed (; ) (Figure 7(c)). Besides, two subtypes of aqueous carbonic FIs were identified (type Ιa and type Ιb). The aqueous carbonic type Ιa and type Ιb inclusions contain two or three phases, but the volumes of the carbonic phases are 10–50% and 50–80%, respectively (Figures 7(b)7(d)). These inclusions have ellipsoidal or irregular shapes with a long axis ranging from 5 to 20 μm (mainly 10–15 μm). Type Ιa inclusions commonly coexist with type ΙΙ inclusions in the early stage. Types Ιa and Ιb inclusions commonly coexist with type ΙΙ and type ΙΙΙ inclusions in the middle stage.

Type ΙΙ inclusions are as follows: these types of inclusions consist of pure CO2 (Figure 7(e)). The pure CO2 inclusions are less than 5% of the total population and mainly occur in the early and middle stages. They are composed of single or two phases (LCO2, VCO2, or ). These types of inclusions are oval to negative crystal morphologies with a long axis of about 5–15 μm, with most measuring between 5 and 10 μm. Type ΙΙ FIs are distributed in clusters with type Ι in quartz grains from the early stage and type Ι and type ΙΙΙ FIs in quartz grains from the middle stage.

Type ΙΙΙ inclusions are as follows: NaCl-H2O-type inclusions consist of H2O liquid and H2O vapor at room temperature and comprise 35% of total inclusions. Three subtypes of aqueous FIs were identified (type ΙΙΙa, type ΙΙΙb, and type ΙΙΙc). The aqueous type ΙΙΙa and type ΙΙΙb inclusions contain two phases, but the volumes of the vapor phase are 10–50% and 50–80%, respectively (Figures 7(f) and 7(g)). Besides, a small number of one-phase aqueous inclusions (type ΙΙΙc) contain liquid water at room temperature; these inclusions are confined to the late stage and occur commonly in clusters with the two-phase type ΙΙΙa inclusions (Figure 8(g)). The type ΙΙΙ inclusions are oval, nearly round, or irregular. The long axis of this type of inclusion ranges from 5 to 20 μm, with most measuring between 8 and 15 μm. Type ΙΙΙ FIs commonly coexist with type Ι and type ΙΙ inclusions in the middle stage.

6.3. Fluid Inclusion Microthermometry

Microthermometric data obtained from primary, pseudosecondary, and secondary type Ι, ΙΙ, and ΙΙΙ FIs of the three mineralization stages are summarized in Table 2 and presented in Figure 9. Microthermometric measurements were not possible for some FIs due to their small sizes (i.e., for type Ι and type ΙΙ, for type Ι, and for type ΙΙΙ).

6.4. Laser Raman Microprobe Analysis

To constrain the fluid inclusion compositions, representative samples of the three different hydrothermal stages were examined using laser Raman microspectroscopy (Figure 10). The type Ι and ΙΙ FIs from the early and middle hydrothermal stages contain a vapor phase of CO2 with minor amounts of CH4 (Figures 10(a)10(d)). The main volatile component of the type ΙΙΙ FIs for the last two stages is H2O with minor amounts of CO2 and H2S (Figures 10(e) and 10(f)).

6.5. Oxygen and Hydrogen Isotopes

The oxygen and hydrogen isotopic results of Xiejiagou gold deposit are given in Table 3, including the data from previous studies [30]. The measured δ18O values for quartz crystals from the gold mineralization at Xiejiagou have intervals between 3.1‰ and 13.2‰. The δ18Owater was calculated according to the oxygen isotopic compositions for quartz and the corresponding homogenization temperatures obtained from our fluid inclusion study (Table 3); average temperatures for early stage, middle stage, and late stage are estimated to be 327, 268, and 183°C, respectively. The calculated δ18Owater values of the fluid range from −9.7‰ to 5.7‰, and the δD values range from −101.8‰ to −83.1‰.

6.6. Sulfur Isotopes

The sulfur isotopic results of Xiejiagou gold deposit are given in Table 4, including the data from previous studies [34]. The sulfur isotope compositions of hydrothermal sulfide have a narrow range of δ34S between 4.7‰ and 7.8‰; the δ34S values of pyrite samples are between 5.9‰ and 7.8‰; the δ34S value of one sphalerite sample is 6.4‰; and the δ34S value of one galena sample is 4.7‰.

6.7. Lead Isotope

Lead isotope data were obtained for sulfides in ore samples from Xiejiagou deposit (Table 5). The 206Pb/204Pb (17.251 to 17.315), 207Pb/204Pb (15.486 to 15.519), and 208Pb/204Pb (37.904 to 38.029) ratios are relatively constant.

7. Discussion

7.1. Evolution of Ore-Forming Fluids

In this study, fluid inclusion microthermometry and laser Raman spectroscopy indicate that dominant type Ιa and rare type ΙΙ FIs are developed in the early-stage quartz crystals. The homogenization temperatures of the early-stage FIs are from 262 to 386°C with salinities of 2.22–8.82 wt% NaCl equivalent. These factors indicate that the initial ore-forming fluids were characterized by medium-high temperatures, low salinities, and enrichment in CO2 with minor CH4 (Table 1, Figures 9 and 11). Therefore, the initial ore fluids belonged to a homogeneous fluid system. Quartz crystals formed during the middle stage of mineralization contain type Ι, ΙΙ, and ΙΙΙ FIs that yield moderate homogenization temperatures (192°C–347°C) and variable salinities (1.02–11.60 wt% NaCl equivalent) (Table 1, Figures 9 and 11). These features indicate that the ore-forming fluids evolved into a system at moderate temperatures and variable salinities. The ore-forming fluid then became an H2O–NaCl system at the end of mineralization as only type ΙΙΙ FIs existed in the late stage. In addition, FIs trapped during the late stage of mineralization have low temperatures (137°C–231°C) and salinities (1.22–7.72 wt% NaCl equivalent) similar to that of meteoric water [7], which indicated that the ore-forming fluids experienced cooling and diluting (Figure 11).

In another way, the average δ18O value of ore-forming fluids slightly decreases from 4.8‰ to 4.2‰ and to −9.7‰ in stages 1–3, as well as the δD values for each stage (median, −86.3‰ → −96.2‰ → −101.8‰). Thus, the O–H isotopic compositions obviously have a shift towards meteoric water from the early to late hydrothermal stage (Figure 12), indicating that significant amount of meteoric water infiltrated into the fluid system during the late mineralization stage.

7.2. Source of Metallic and Hydrothermal Components of the Ore-Forming Fluids

Sulfur isotopic composition of sulfide minerals can be used to trace the source of sulfur in ore fluids [64, 65]. The δ34S values of the sulfides in Xiejiagou deposit display a narrow range (4.7‰–7.8‰, with an average of 6.6‰, Figure 13), suggesting an extremely homogeneous source; this is similar to the majority of gold deposits in Jiaodong [2, 3, 66]. The average δ34S values of those sulfide minerals from Xiejiagou deposit show a trend (Figure 13, Table 4), suggesting equilibrium fractionation of sulfur isotopes among the sulfide minerals. Based on the average homogenization temperatures of the corresponding fluid inclusions from the first two stages (327°C and 268°C), the δ34SH2S values of the mineralizing fluid are 4.5 to 6.9‰ using the equations of Ohmoto and Rye [67]. These δ34SH2S values are obviously higher than that of chondrite () indicating that the sulfur is unlikely to have originated from the mantle. In addition, the average δ34SH2S value is within the typical range of the orogenic gold deposits worldwide (Figure 13, 0 to 10‰) and also within the ranges of whole rocks and pyrite from the Archean Jiaodong group, the Mesozoic Linglong and Guojialing granites, and the Early Cretaceous mafic-intermediate dikes. Even though the average δ34SH2S value of the Xiejiagou deposit is close to that of the Late Jurassic Linglong granite (host rock, pre-ore), it is still difficult to establish whether the ore-forming fluids had a magmatic source or other sources.

Lead isotope compositions (particularly of sulfide minerals) are an excellent proxy for constraining the sources of lead in ore-forming fluids [68, 69]. Sulfide minerals in the Xiejiagou deposit show Pb isotope compositions of , , and . In a similar manner to S isotope data, all of the Pb isotope data overlap with the compositions of whole rocks from Mesozoic granites and dikes, as well as Precambrian metamorphic rocks (Figure 14), and are consistent with the lead isotope compositions of most gold deposits in Jiaodong [19]. These values overlap the field of Late Jurassic Linglong granite, which is the main host rock of the Xiejiagou deposit. Thus, the lead isotopic compositions of sulfides can be interpreted to suggest that the host rocks maybe one of the lead sources for sulfides but can not be used to trace the source of gold. In another way, on plumbotectonic diagrams (Figure 14, Table 5), the Pb isotope compositions of sulfide minerals fall within a range that spans within the orogene and lower crust lead isotope evolution curves [70], indicating that the lead in ore-forming fluids may be originated from the lower crust reservoir in an orogenic regime. The Jiaodong gold mineralization had a genetic relationship with continental extension that was induced by the upwelling of hot felsic magma during the Mesozoic [8, 71, 72]. Therefore, the lead in the ore-forming fluids, and included lead in the Mesozoic granites and dikes, may be originated from the Mesozoic orogenic lead reservoir.

The oxygen and hydrogen isotopic compositions of hydrothermal quartz crystals are useful tracers for determining the source of ore-forming fluids. Unlike hydrothermal minerals that form during the late stage of mineralization, those that form during the early and middle stages preserve O and H isotope ratios that should not have been significantly influenced by meteoric water and thus have the potential to accurately reflect the nature of the original ore-forming fluids [26, 73]. The O–H compositions of ore-forming fluids are shown in Figure 12. The calculated δ18OH2O (1.9‰ to 5.7‰) and δD (–99.4‰ to –83.1‰) values for fluid of the first two stages are slightly lower than the magmatic water composition ([74], −50‰ to −85‰ for δD; 5.5‰ to 9.0‰ for δ18OH2O). As such, it is difficult to establish whether the ore-forming fluids had a magmatic source. Mesozoic granites are widespread on the Jiaodong Peninsula (Figure 1), and many researchers have therefore suggested that magmatic water is likely to be the dominant source of the fluids that formed the Jiaodong gold deposits [16, 17, 28]. However, FIs of the Xiejiagou deposit are characteristic of low-salinity mesothermal CO2-rich fluids, which are widely accepted to be of metamorphic origin and differ markedly from typical high-temperature, high-salinity exsolved magmatic fluids [73, 75]. Furthermore, because the mineralization of the Xiejiagou deposit occurred at about 123.6–115.2 Ma, it is reasonable to narrow the field of magmatic water compositions to include only those of Early Cretaceous Guojialing granite (132–126 Ma; [39]) in the Jiaodong Peninsula. None of the analyzed compositions fall within this narrower field ([76]; magmatic water of Guojialing granite: −50‰ to −85‰ for δD; 5.5‰ to 9.0‰ for δ18OH2O), indicating that magmatic water may not to be the source of the fluids that formed the Xiejiagou deposit. On the other hand, it is worth noting that the ore-forming fluids of reduced intrusion-related gold deposit are also characterized by CO2-rich fluids with variable temperature and salinity [77, 78]. This raises the possibility that the Xiejiagou deposit is a reduced intrusion-related gold deposit. Unfortunately, as with previous researchers, we were unable to find any syn-ore-reduced igneous rocks in the Jiaodong Peninsula and we did not observe any of the exsolved textures (e.g., unidirectional solidification textures and quartz eyes) that are typical of primary magmatic water [77]. There is also no evidence of sheet-formed gold-bearing veins ([78]; representative mineralization style of reduced intrusion-related gold deposit) in the Xiejiagou deposit or, for that matter, in any other gold deposits within the Jiaodong Peninsula. In combination, these factors indicate that Xiejiagou is not a reduced intrusion-related deposit.

Representatively, because the nature of ore-forming fluids of the gold deposits in the Jiaodong Peninsula was mesothermal, has low salinity, and enriched in CO2, it has been suggested that the Jiaodong gold deposits including Xiejiagou might be orogenic gold deposit [21] and that the ore-forming fluids were of metamorphic origin. In addition, the secondary FIs within quartz crystals of the early and middle stages (Table 1; Figure 9) probably lead to inaccurate δD values for the hydrothermal fluids in each stage [22, 7981], making the calculated O–H values slightly less than those of the typical range of the majority of orogenic gold deposit (Figure 12; [82]). Therefore, the true δD values for the hydrothermal fluids are higher than the δD values of FIs and may come to lie in the metamorphic water field.

In summary, we speculate that the most probable explanation is that the fluids that formed the Xiejiagou deposit were dominated by a metamorphic source. It is commonly believed that orogenic gold deposits, which are predominantly hosted by greenschist facies or low-amphibolite facies metamorphic terranes, have gold-forming ages that are either synchronous with, or slightly younger than, terrane metamorphism and deformational events [82]. However, the Jiaodong gold deposits including Xiejiagou are of the granite-hosted “gold-only” type [1], with only a limited number of metamorphic rocks within the vicinity of the ore deposit (Figure 2). In addition, the ages of mineralization at the Jiaodong Peninsula postdate peak at high amphibolite to granulite facies metamorphism in the area by at least 1.7 Ga, indicating that the source of ore-forming fluids can not come from the Neoarchean and Paleoproterozoic basement rocks [22]. Therefore, the original ore fluids may in fact have formed via metamorphic dehydration, decarbonization, and desulfidation of the subducting Paleo-Pacific slab within an accretionary orogeny [14, 15, 21, 22].

7.3. Trapping Pressure of FIs and Metallogenic Depth

The assemblage of the fluid inclusions from the different hydrothermal stages can estimate the pressure conditions of the fluid inclusions trapped during ore forming [25, 26]. Given that the type Ι FIs are common in the early and middle hydrothermal stages, these FIs exhibit values of approximately -56.6°C; as such, the pressure can be estimated based on the CO2–H2O–NaCl system with a range of isochores calculated using FLINCOR software [54] and the formula of Burke [56].

Based on the homogenization temperatures and salinity of the type Ι FIs in the first two stages, a representative isochore calculated for the H2O–CO2–NaCl system containing 6 wt% NaCl equivalent and 15 mol% CO2 ([55], Figure 15) was constructed to estimate the trapping pressures of type Ι FIs from the early and middle stages. Given that the homogenization temperatures of FIs from the first two stages have wide ranges and can not provide the accurate temperature and pressure of fluid entrapment, upper and lower quartiles are more suitable than the average values to reflect the data distribution [80, 83]. Therefore, intersections of isochores of FIs with the upper and lower quartiles of homogenization temperatures and bulk densities are used to constrain the trapping pressure range of the fluids.

The trapping pressure for FIs from quartz veinlets of the early hydrothermal stage were calculated to be 224–302 MPa using homogenization temperatures ranging from 305 to 359°C and CO2 densities ranging from 0.574 to 0.751 g/cm3 (Figure 15(a)). The pressure of the gold polymetallic quartz veinlets from the middle hydrothermal stage was estimated ranging from 191 to 258 MPa by homogenization temperatures of 246 to 289°C and CO2 densities ranging from 0.589 to 0.767 g/cm3 (Figure 15(b)). Applying an ancient ground pressure gradient of 0.0265 GPa/km (mean density of upper crust), the calculated trapping pressure for the gold polymetallic quartz veinlets corresponds to a depth of 7.2 to 9.7 km within a lithostatic fluid system. The thermochronologic study indicates that the erosion depth of the Jiaobei uplift has been calculated as about 3.0 to 7.0 km between the Early Cretaceous and the present [84]. Consequently, it is inferred that the Xiejiagou gold deposit has undergone very less total denudation since its formation. Accordingly, it can be predicted that large gold reserves could still exist at deep levels within the Xiejiagou deposit.

7.4. Mechanisms of Gold Transport and Deposition

Gold is commonly dissolved and transported in the form of gold bisulfide complexes (Au(HS)0 and Au(HS)2) in hydrothermal solutions at temperatures of <400°C [85]. The aqueous speciation of gold(I) sulfide complexes is sensitive to physical and chemical conditions of solution [86]. For example, when the temperature is below 400°C, Au(HS)0 dominates under acidic conditions, whereas Au(HS)2– is more likely to form under weakly acidic to neutral conditions [8688]. The alteration mineral assemblages of the Xiejiagou deposit consist of quartz sericite K-feldspar chlorite association (Figures 3 and 4), indicating that the pH of the ore solution was near-neutral to weakly acidic [89, 90]. Thus, gold was most probably transported as the Au(HS)2 complex at the Xiejiagou deposit. It is also supported by the observation that gold is generally accompanied by pyrite and other sulfides in the Xiejiagou deposit (Figures 3(k) and 3(l)), as is typical in orogenic gold deposits worldwide, in which gold(I) sulfide complexes (especially Au(HS)2) were considered the most likely species for transporting gold [82]. Besides, the abundance of type Ι and ΙΙ FIs in quartz from the early hydrothermal stage demonstrates that the ore-forming solution was enriched in CO2. CO2 can buffer or neutralize the fluid within the pH range in which the gold bisulfide complexes are stable and increase its solubility in the ore-forming fluids [91], which provides favorable conditions for gold bisulfide migration. Therefore, the initial auriferous hydrothermal fluids need changes of physical and chemical conditions to reduce solubility of Au(HS)2 so that gold can be precipitated.

In the pressure-temperature isochore diagram (Figure 15(b)), the P–T window with the yellow-shaded area of the middle stage is mainly below the solvus (6 wt% NaCl equivalent with 15 mol% CO2), indicating that the initial homogeneous fluids entered the two-phase field and experienced phase separation process in this stage [55]. Besides, the existence of the types Ι, ΙΙ and ΙΙΙ FIs within a single thin-section of the middle stage with different compositions and densities is inferred to be caused by the CO2 effervescence from an initial homogeneous fluid system [92, 93]. In addition, the homogenization temperatures of type ΙΙΙ FIs are slightly lower than those of coexisting type Ι FIs within middle-stage samples, also because the homogenization temperatures of aqueous inclusions can thus be significantly lowered by CO2 effervescence resulting from pressure fluctuations [94]. Such CO2 effervescence in auriferous hydrothermal systems could have led to the exsolution of the volatiles such as H2S (Figure 10(e)). The exsolution of carbon dioxide of the original ore fluid led to an increase in pH values and a decrease in temperature. This process, in conjunction with the exsolution of H2S from the ore-forming fluid, which may have made gold bisulfide complexes unstable and reduced their solubility, resulted in the precipitation of abundant gold.

8. Conclusions

(1)Three types of primary fluid inclusion were identified in the Xiejiagou gold deposit: type Ι (aqueous carbonic), type ΙΙ (pure carbonic), and type ΙΙΙ (aqueous). Type Ι and minor type ΙΙ inclusions exist in quartz from the early stage of mineralization, whereas all three types of inclusions are observed in quartz from the middle period of mineralization. In contrast, quartz from the late stage of mineralization only contains type ΙΙΙ inclusions(2)The initial ore-forming fluids of the Xiejiagou deposit belonged to a medium-high temperature (262–386°C), CO2-rich, low salinity (2.22–8.82 wt% NaCl), and homogeneous fluid system. During mineralization, the fluid finally evolved into a medium-low temperature NaCl–H2O system as a result of an influx of meteoric water. Considering the fluid inclusion characteristics, H–O–S–Pb isotopes, and regional geological events, the ore-forming fluid reservoir of the Xiejiagou deposit was likely metamorphic in origin(3)Trapping pressures of the first two hydrothermal stages estimated from the type Ι inclusion assemblages were ~224–302 MPa and ~191–258 MPa, respectively. The gold mineralization of the Xiejiagou gold deposit occurred at a lithostatic depth of ~7.2–9.7 km(4)Au(HS)2 was the most probable gold-transporting complex at the Xiejiagou deposit. CO2 effervescence in auriferous hydrothermal systems resulted in the precipitation of gold

Data Availability

The manuscript is a data self-contained article, whose results were obtained from the laboratory analysis, and the entire data is presented within the article.

Conflicts of Interest

We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the submitted paper.

Acknowledgments

This research work was jointly supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant no. XDA20070304), the China Geological Survey (DD20160024 and grant no. 121201102000150011), the National Natural Science Foundation of China (grant nos. 41320104004 and 41602084), and the China Postdoctoral Science Foundation funded project (grant no. 2016M590119). We especially thank four reviewers for their kind and critically constructive comments and suggestions, which greatly improved the quality of our manuscript. We are grateful to the Institute of Mineral Resources, Chinese Academy of Geological Science, for laser Raman spectroscopic analyses. We are also grateful to the Beijing Research Institute of Uranium Geology for H–O–S–Pb isotopic analyses.