Geochronology and Genesis of the Xitian WSn Polymetallic Deposit in Eastern Hunan Province , South China : Evidence from Zircon UPb and Muscovite ArAr Dating , Petrochemistry , and Wolframite Sr-Nd-Pb Isotopes

The recently explored Xitian tungsten-tin (W-Sn) polymetallic ore field, located in Hunan province, South China, is one of the largest ore fields in the Nanling Range (NLR). Two major metallogenic types appeared in this ore field, skarnand quartz vein-type. They are distributed within Longshang, Heshuxia, Shaiheling, Hejiangkou, Goudalan, and so on. Hydrothermal zircons from two altered granites yielded U-Pb ages of 152.8 ± 1.1 Ma, and 226.0 ± 2.8 Ma, respectively. Two muscovite samples from ore-bearing quartz vein yielded 40Ar/39Ar plateau ages of 156.6 ± 0.7 Ma, 149.5 ± 0.8 Ma, respectively. Combined with the geological evidence, two metallogenic events are proposed in the Xitian ore field, with skarn-type W-Sn mineralization in Late Triassic (Indosinian) and quartz vein/greisen type W-Sn mineralization in Late Jurassic (Yanshanian). The relatively low Ce/Ce* ratios and high Y/Ho ratios in zircons from two altered granites indicate that the hydrothermal fluids of two metallogenic events are characterized by low oxygen fugacities and enrichment in F. The similar chondrite-normalized patterns between the skarn and Xitian Indosinian granites and Sr-Nd-Pb isotopic compositions of wolframite suggest that the metal sources for both types W-Sn mineralization are derived from a crustal source.

The Xitian W-Sn polymetallic ore field, located in the middle of the Nanling region, is one of the largest newly discovered ore fields in recent years (Figure 2).A large number of studies have been carried out in the Xitian area by geochemical and isotopic methods, and these studies have shown that the formation of this deposit is genetically related to the Xitian pluton [19,28,38,39].Many dating technologies have been applied to study the emplacement age of this pluton and metallogenic age of this deposit, including LA-ICPMS, and ion probe by either SHRIMP or CAMECA zircon U-Pb, mica and cassiterite 40 Ar- 39 Ar, and molybdenite Re-Os isotopic techniques [19, [40][41][42].These precise data provide detailed chronological constraints on the emplacement age of the Xitian pluton and for the time interval between W-Sn mineralization of the Xitian deposit.The majority of these chronological data show that the Xitian pluton could be subdivided into the Late Triassic (Indosinian) granites (230-220 Ma) and the Late Jurassic (Yanshanian) granites (160-140 Ma), and the time interval for W-Sn  [2]), showing the distribution of granitic plutons, basalts, and related W-Sn deposits and their geochronological data compiled from [17][18][19]21,29,[31][32][33][34][35][36][37].TB: Tarim block; CAB: Cathaysian Block; NCB: North China Block; SCB: South China Block; YZB: Yangtze Block.
The Xitian W-Sn polymetallic ore field, located in the middle of the Nanling region, is one of the largest newly discovered ore fields in recent years (Figure 2).A large number of studies have been carried out in the Xitian area by geochemical and isotopic methods, and these studies have shown that the formation of this deposit is genetically related to the Xitian pluton [19,28,38,39].Many dating technologies have been applied to study the emplacement age of this pluton and metallogenic age of this deposit, including LA-ICP-MS, and ion probe by either SHRIMP or CAMECA zircon U-Pb, mica and cassiterite 40 Ar- 39 Ar, and molybdenite Re-Os isotopic techniques [19, [40][41][42].These precise data provide detailed chronological constraints on the emplacement age of the Xitian pluton and for the time interval between W-Sn mineralization of the Xitian deposit.The majority of these chronological data show that the Xitian pluton could be subdivided into the Late Triassic (Indosinian) granites (230-220 Ma) and the Late Jurassic (Yanshanian) granites (160-140 Ma), and the time interval for W-Sn mineralization is 160-150 Ma.From the above-mentioned evidence, it seems that the formation of this ore field is attributed to Yanshanian magmatic activity.Deng et al. [43] obtained a molybdenite Re-Os age of 225.5 ± 3. 6 Ma from altered granites in the Indosinian granitic batholith, indicating a possibility for the Indosinian mineralization event in the Xitian ore field; however, this age may not represent the age of the large-scale skarn-type W-Sn mineralization in the Xitian ore field.It was proposed that the skarn occurred in the contact zone between the Yanshanian granites and the Devonian dolomitic limestone [19].However, recent studies show that the granites belong to Indosinian granites rather than Yanshanian granites [34].Therefore, is the skarn type ore body related to the Indosinian granites, rather than to the Yanshanian granites?Ore-forming age of some deposits in the Xitian ore field is still in doubt; for example, the Hejiangkou deposit.In this study, we display the results of zircon U-Pb dating, zircon compositions, muscovite 40 Ar-39 Ar dating, skarn geochemistry, and wolframite Sr-Nd-Pb isotopic compositions, with the aims of constraining the time interval between mineralization and the emplacement of associated granitic rocks, outlining the genetic relationship between two episodes of granitic magmatism and two types of W-Sn mineralization, and probing into the genesis of the two types of W-Sn mineralization in the Xitian ore field.
Minerals 2018, 8, x FOR PEER REVIEW 3 of 23 mineralization is 160-150 Ma.From the above-mentioned evidence, it seems that the formation of this ore field is attributed to Yanshanian magmatic activity.Deng et al. [43] obtained a molybdenite Re-Os age of 225.5 ± 3. 6 Ma from altered granites in the Indosinian granitic batholith, indicating a possibility for the Indosinian mineralization event in the Xitian ore field; however, this age may not represent the age of the large-scale skarn-type W-Sn mineralization in the Xitian ore field.It was proposed that the skarn occurred in the contact zone between the Yanshanian granites and the Devonian dolomitic limestone [19].However, recent studies show that the granites belong to Indosinian granites rather than Yanshanian granites [34].Therefore, is the skarn type ore body related to the Indosinian granites, rather than to the Yanshanian granites?Ore-forming age of some deposits in the Xitian ore field is still in doubt; for example, the Hejiangkou deposit.In this study, we display the results of zircon U-Pb dating, zircon compositions, muscovite 40 Ar-39 Ar dating, skarn geochemistry, and wolframite Sr-Nd-Pb isotopic compositions, with the aims of constraining the time interval between mineralization and the emplacement of associated granitic rocks, outlining the genetic relationship between two episodes of granitic magmatism and two types of W-Sn mineralization, and probing into the genesis of the two types of W-Sn mineralization in the Xitian ore field.

Geological Setting
The South China Block (SCB) was formed by the amalgamation of the Yangtze Block (YZB) to the northwest and the Cathaysian Block (CAB) to the southeast at ca. 820 Ma (Figure 1a; [45]).The Nanling region, located in the central section of SCB, is comprised of Guangxi, Guangdong, Hunan, and Jiangxi province occupying an area of 170,000 km 2 [46].This region has undergone several significant tectonic-magmatic events, the most famous of which are the Indosinian and Yanshanian tectonic events during the Mesozoic [47][48][49][50][51][52][53].Due to the superior metallogenic geotectonic setting, it is characterized by widespread igneous rocks and numerous large-scale W-Sn polymetallic deposits [2,3,16,20,54,55].The basement of the Nanling region consists of weakly metamorphosed Precambrian, late Paleozoic sedimentary strata which are mainly Devonian and Carboniferous carbonate rocks, and lesser amounts of Upper Triassic to Tertiary sandstone and siltstone [2,46].The regional fault is the NE-trending Chenzhou-Linwu fault which controls the spatial distribution of the granitic intrusions and numerous W-Sn polymetallic deposits associated with the granitic magma activities (Figure 1b; [56,57]).Numerous granitic intrusions were emplaced in this region, most of which are Indosinian and Yanshanian pluton, and these granites are mostly peraluminous, calc-alkaline and remelted granites [12,23,27,28,58,59].
The Xitian ore field, located in Chaling, Hunan province, is characterized by intensive and widely distributed granitoids associated with numerous non-ferrous and rare-metal minerals of Mesozoic age [19,39,44,60].

Sedimentary Rocks
The strata outcropping in the Xitian area are Ordovician, Devonian, Carboniferous, and Cretaceous sedimentary rocks, among which the middle to upper Devonian and Carboniferous rocks are dominant (Figure 2).Lying unconformably on Ordovician metasedimentary rocks, the Devonian strata can be subdivided into the Middle Devonian Tiaomajian and Qiziqiao Formations, and the upper Devonian Shetianqiao and Xikuangshan Formations [61].The Tiaomajian formation, 35-42 m in thickness, consists of conglomerate-bearing quartzite.The Qiziqiao formation, over 200 m thick, comprises impure carbonate rocks and arenaceous shale, and is the typical ore reservoirs of the Xitian ore field.The Shetianqiao formation is up to 500 m thick, composed of quartz sandstone and argillaceous siltstone.The Xikuangshan formation, 110-130 m in thickness, are mainly quartz sandstone, arenaceous shale and nodular limestones.The Carboniferous Yanguan formation which is about 275 m in thickness, mainly consists of sandshale and siltstone.

Structure
The Xitian ore field is located to the east of the NE-trending Chenzhou-Linwu deep fault (Figure 1b), which is considered to be the boundary between the Yangtze Block and Cathaysian Block [56,57,62].The tectonic framework of this region is controlled mainly by two trends of faults which are approximately NE-, nearly SN-and NW-trending.The NE-trending faults are the larger in scale, and some of these faults are truncated by the NW-and/or SN trending-faults (Figure 2).The NE-trending faults are the main ore-controlling faults in Xitian ore field, with 2-13 km in outcropped length, 60-70 • in angle of trend, and 60-85 • in angle of dip [60].The nearly SN-trending faults are also important ore-bearing structures, including a series of NNW-, SN-and NNE-trending small faults [60].The NW-trending faults are about 1.5-8.0km in outcropped length, with a dip of NNE and large inclined angle [60].

Geology of the Ore Deposits
Previous studies have revealed that four types of W-Sn polymetallic ore bodies were exploited in the Xitian ore field consisting of skarn-, quartz vein-, greisen-and structurally altered rock-types [19,60].The skarn-type ore bodies, occurred mainly in Longshang, Hejiangkou, and Shaiheling, are characterized by W-Sn mineralization (Figure 2).The quartz vein-and greisen-type ore bodies are also characterized by W-Sn mineralization, distributed in Longshang, Hejiangkou, Heshuxia, and Goudalan (Figure 2).The structurally altered rock-type ore bodies are mainly found in Shaiheling featured by Lead (Pb)-Zinc (Zn) mineralization (Figure 2).The morphology of ore body, specimen and micrographs were presented in Figures 3 and 4, respectively.

Geology of the Ore Deposits
Previous studies have revealed that four types of W-Sn polymetallic ore bodies were exploited in the Xitian ore field consisting of skarn-, quartz vein-, greisen-and structurally altered rock-types [19,60].The skarn-type ore bodies, occurred mainly in Longshang, Hejiangkou, and Shaiheling, are characterized by W-Sn mineralization (Figure 2).The quartz vein-and greisen-type ore bodies are also characterized by W-Sn mineralization, distributed in Longshang, Hejiangkou, Heshuxia, and Goudalan (Figure 2).The structurally altered rock-type ore bodies are mainly found in Shaiheling featured by Lead (Pb)-Zinc (Zn) mineralization (Figure 2).The morphology of ore body, specimen and micrographs were presented in Figures 3 and 4, respectively.

Longshang
The Longshang deposit, located in the western part of the Xitian ore field, is the largest W-Sn deposit in this area (Figure 2).Two types of ore body are exploited in this deposit, which are skarnand quartz vein type.The skarn-type W-Sn ore bodies are stratiform, stratoid and lentoid, mainly occurring at the endo-or exo-contact zone between the Xitian Triassic granites and the Devonian dolomitic limestone (Figure 3a).In addition, the N-S-trending and E-dipping skarn-type ore bodies are 2700 m long and 4.5-33.1 m thick with ore grade of 0.14-0.77%Sn and 0.038-0.83%WO3 [35].The quartz vein-type ore bodies, occurring in the Xitian Triassic granites, are E-W-trending and S-dipping with ore veins of 3-20 cm in thickness (Figure 3b).A complex skarn-vein type W-Sn ore bodies, newly found in the Longshang deposit, have relatively high grade with 2.3% WO3 and 1.7% Sn in some parts of the ore bodies (Figure 3c).
The major ore minerals of skarn-type ore bodies comprise scheelite, pyrite, and some other minerals, with gangue mineral mainly consisting of garnet, epidote, idocrase, quartz, and other skarn minerals (Figure 4a,b,f,g).The major ore minerals in quartz vein-type ore bodies are wolframite, cassiterite, scheelite, pyrite, and some other minerals, with gangue mineral mainly consisting of quartz, fluorite, calcite, and sericite.

Longshang
The Longshang deposit, located in the western part of the Xitian ore field, is the largest W-Sn deposit in this area (Figure 2).Two types of ore body are exploited in this deposit, which are skarn-and quartz vein type.The skarn-type W-Sn ore bodies are stratiform, stratoid and lentoid, mainly occurring at the endo-or exo-contact zone between the Xitian Triassic granites and the Devonian dolomitic limestone (Figure 3a).In addition, the N-S-trending and E-dipping skarn-type ore bodies are 2700 m long and 4.5-33.1 m thick with ore grade of 0.14-0.77%Sn and 0.038-0.83%WO 3 [35].The quartz vein-type ore bodies, occurring in the Xitian Triassic granites, are E-W-trending and S-dipping with ore veins of 3-20 cm in thickness (Figure 3b).A complex skarn-vein type W-Sn ore bodies, newly found in the Longshang deposit, have relatively high grade with 2.3% WO 3 and 1.7% Sn in some parts of the ore bodies (Figure 3c).
The major ore minerals of skarn-type ore bodies comprise scheelite, pyrite, and some other minerals, with gangue mineral mainly consisting of garnet, epidote, idocrase, quartz, and other skarn minerals (Figure 4a,b,f,g).The major ore minerals in quartz vein-type ore bodies are wolframite, cassiterite, scheelite, pyrite, and some other minerals, with gangue mineral mainly consisting of quartz, fluorite, calcite, and sericite.

Hejiangkou
The Hejiangkou deposit is also located in the western part of the Xitian ore field (Figure 2).Similar to the Longshang deposit, Hejiangkou is also characterized by skarn-type Sn-Cu mineralization and quartz vein type W-Sn mineralization.The stratiform and lentoid skarn-type ore bodies are the mainly mining target in this deposit, which occur in the endo-or exo-contact zone between the Xitian Triassic granites and the Devonian limestone, with 500-1200 m long and 30-50 m thick.The quartz vein-type W-Sn ore veins, hosted in the Xitian Triassic granites and Devonian sandstone, are NEE and/or NNW-trending and N-dipping with ore veins of 50-200 m long and 0.1-0.8m thick (Figure 3d).
Ore minerals of the skarn-type ore bodies are mainly cassiterite, chalcopyrite, scheelite, sphalerite, and pyrite, with gangue minerals consisting of quartz, calcite, garnet, epidote, idocrase and other skarn minerals.Ore minerals of the quartz vein-type ore bodies are mainly composed of wolframite, cassiterite, scheelite, and pyrite, and gangue minerals are quartz, fluorite, and muscovite.

Heshuxia
The Heshuxia deposit is in the eastern part of the Xitian ore field (Figure 2).Unlike the Longshang and Hejiangkou deposits in the western part of the Xitian ore field, the Heshuxia deposit is mainly characterized by quartz vein type of W mineralization.The NW-trending and SW-dipping ore veins are mainly hosted in the Xitian Triassic granites, with length of 100-1000 m, thickness of 30-100 cm, and ore grade of 0.172-0.700%WO 3 (Figure 3e,f).Ore minerals are mainly wolframite, with a small quantity of scheelite, molybdenite, pyrite, and chalcopyrite.Gangue minerals are mainly composed of quartz, fluorite, and muscovite (Figure 4c,d,e,h,i).

Goudalan
The Goudalan deposit, located in the southeast part of the Xitian ore field, is also characterized by quartz vein type of W-Sn mineralization.The NEE trending ore veins are also hosted in the Xitian Triassic granites with length of 100-1000 m, thickness of 0.1-1 m and ore grade of 0.086-0.762%WO 3 (Figure 3g).Ore minerals mainly consist of wolframite, with a small quantity of scheelite, molybdenite, pyrite, and chalcopyrite.Gangue minerals are mainly composed of quartz, fluorite, and muscovite.

Shaiheling
The Shaheling deposit, located in the northeast part of the Xitian ore field, are characterized by the structurally altered rock-type of Pb-Zn mineralization and skarn-type of W-Sn mineralization.The structurally altered rock-type of Pb-Zn ore bodies are mainly hosted in the fracture zone of the skarn and/or carbonate formations, with NW trending and NE dipping (Figure 3h).The skarn-type W-Sn ore bodies are stratiform, stratoid and lentoid, mainly occurring at the endo-or exo-contact zone between the Xitian Triassic granites and the Devonian dolomitic limestone as in the Longshang and Hejiangkou deposits, with average length of 1.4 km and thickness of 3.3 m (Figure 3i).The ore minerals of the structurally altered rock-type ore mainly comprise of sphalerite, galenite, scheelite, pyrite and chalcopyrite, with gangue minerals consisting of quartz, feldspar, chlorite, and so on.Ore minerals of the skarn-type W-Sn ore are mainly composed of scheelite, cassiterite, pyrite and chalcopyrite, with gangue minerals consisting of quartz, garnet, epidote, and so on.

Sampling and Analytical Methods
The analyzed samples were collected from underground mines (Figure 2).Zircon grains used for LA-ICP-MS U-Pb dating were separated from a sericitic coarse-grained biotite granite (sample No. 19-4s1, Figure 5a,b) and a sericitic fine-grained two mica granite (sample No. 24-15s1, Figure 5c,d which were collected from the Longshang deposit.The muscovite used for Ar-Ar dating were extracted from quartz vein-type ore in Hejiangkou (Figure 3d, sample 1401-3-1) and Heshuxia (Figure 3e, sample 1401-5-3).The skarn used for geochemical analyses were all endo-skarn, some of which contained sulfides, collected from the Longshang, Shaiheling and Huamu deposits.The wolframite, separated from the quartz vein type ore and used for Sr-Nd-Pb analysis, were collected from Longshang, Goudalan and Heshuxia.The sampling locations were marked in Figure 2.

In Situ LA-ICPMS Zircon U-Pb Dating and Trace Element Compositions
Zircon grains from these samples were separated by conventional magnetic and heavy liquid techniques before they were hand-picked under a binocular microscope.They were then mounted into epoxy resin blocks and polished to obtain flat surfaces.Cathodoluminescence (CL) imaging technique was used to visualize the internal structures of individual zircon grains, with a scanning electron microscope (TESCAN MIRA 3 LMH FE-SEM, TESCAN, Brno, Czech Republic) at the Sample Solution Analytical Technology Co., Ltd., Wuhan, China.Zircon U-Pb dating was undertaken with an Agilent 7700 inductively coupled plasma-mass spectrometer (ICP-MS, Agilent, Santa Clara, CA, USA), combined with a Coherent 193 laser ablation (LA) system at Sample Solution Analytical Technology Co., Ltd., Wuhan, China.Two zircon standards, 91500 (1062 ± 4 Ma; [64] and GJ-1 (610.0 ± 1.7 Ma; [65], were used as external standards for dating.Standard silicate glass (NIST SRM610) was used for external standardization for trace element analysis, and 29 Si was used for internal standardization (32.8% SiO2 in zircon).The standard protocol correction method was used in analyzing the 91500 and GJ-1 standard zircons twice and once, respectively, after every five analyses.The raw ICP-MS data were processed using ICPMSDataCal software [66], and common Pb was corrected following [67].Concordia diagrams and weighted mean calculations were processed using Isoplot (version 3.0; [68]).

In Situ LA-ICP-MS Zircon U-Pb Dating and Trace Element Compositions
Zircon grains from these samples were separated by conventional magnetic and heavy liquid techniques before they were hand-picked under a binocular microscope.They were then mounted into epoxy resin blocks and polished to obtain flat surfaces.Cathodoluminescence (CL) imaging technique was used to visualize the internal structures of individual zircon grains, with a scanning electron microscope (TESCAN MIRA 3 LMH FE-SEM, TESCAN, Brno, Czech Republic) at the Sample Solution Analytical Technology Co., Ltd., Wuhan, China.Zircon U-Pb dating was undertaken with an Agilent 7700 inductively coupled plasma-mass spectrometer (ICP-MS, Agilent, Santa Clara, CA, USA), combined with a Coherent 193 laser ablation (LA) system at Sample Solution Analytical Technology Co., Ltd., Wuhan, China.Two zircon standards, 91500 (1062 ± 4 Ma; [64] and GJ-1 (610.0 ± 1.7 Ma; [65], were used as external standards for dating.Standard silicate glass (NIST SRM610) was used for external standardization for trace element analysis, and 29 Si was used for internal standardization (32.8% SiO 2 in zircon).The standard protocol correction method was used in analyzing the 91500 and GJ-1 standard zircons twice and once, respectively, after every five analyses.The raw ICP-MS data were processed using ICPMSDataCal software [66], and common Pb was corrected following [67].Concordia diagrams and weighted mean calculations were processed using Isoplot (version 3.0; [68]).

Muscovite 40 Ar-39 Ar Dating
The Muscovite grains were carefully handpicked using a binocular microscope from the crushed sample to ensure purity up to 99.9%, then these grains were washed repeatedly in an ultrasonic bath using deionized water and acetone.Aliquots of approximately 10 mg were wrapped in Al foil and stacked in quartz vials.After samples had been stacked, the sealed quartz vials were put in a quartz canister, which was wrapped with cadmium foil (0.5 mm in thickness) to act as a slow neutron shield thereby preventing interface reactions during irradiation.The irradiation procedure was put the samples in channel B4 of Beijing 49-2 reactor for 50 h at the Chinese Academy of Nuclear Energy Sciences.During irradiation, the vials were rotated at a speed of two cycles per minute to ensure uniformity of the irradiation.The biotite standard ZBH-2506 (132.5 Ma; [69]) was used to monitor the neutron flux. 40Ar/ 39 Ar stepwise heating analyses were performed at the Key Laboratory of Tectonics and Petroleum Resources, China University of Geosciences, Wuhan, China.Analyses were carried out using an Argus VI mass spectrometer combined with Coherent 50 W CO 2 laser system.The time of heating was 60 s for every single stage with a laser beam diameter of 2.5 mm, and the time of gas purification was 400 s with two Zr-Al scavenger.The detailed analytical procedures were given by [70].K 2 SO 4 and CaF 2 crystals were analyzed to calculate Ca, K correction factors: ( 39 Ar/ 37 Ar) Ca = 8.984 × 10 −4 , ( 36 Ar/ 37 Ar) Ca = 2.673 × 10 −4 , ( 40 Ar/ 39 Ar) K = 5.97 × 10 −3 .The data-processing software and diagrams of plateau age we used was the ArArCALC 2.52 software by [71].

Skarn Major and Trace Elements Analysis
The skarn samples were crushed in a milling machine to 200 mesh before elemental analyses were conducted.The major and trace element compositions of skarn were analyzed at ALS Chemex, Guangzhou, China.The major element contents were measured using a Panalytical Axios Max X-ray fluorescence (XRF, Panalytical, Almelo, The Netherlands) instrument, with analytical accuracy of about 1-5%.Trace element compositions were measured using ICP-MS (Perkin Elmer Elan 9000, Perkin, Waltham, MA, USA), with analytical accuracy of better than 5%.

Wolframite Sr-Nd-Pb Isotopic Composition Analysis
Sr-Nd-Pb isotopic analyses were carried out at the Key Laboratory of Crust-Mantle Materials and Environments, Chinese Academy of Sciences, University of Science and Technology of China, Hefei, China, using a Finnigan MAT-262 multicollector thermal ionization mass spectrometer (MC-TIMS).Rb-Sr was separated and purified using conventional cation exchange (AG50W-X12, 200-400 resin), whereas Sm and Nd were separated and purified using Teflon and a Power resin, respectively.The correction for mass fractionation of the Sr-Nd isotopic ratio was undertaken by normalizing to 86 Sr/ 88 Sr = 0.1194 and 146 Nd/ 144 Nd = 0.7219.The Sr standard (NBS987, 87 Sr/ 86 Sr = 0.710249 ± 0.000012 (2σ)) and the Nd standard (La Jolla, 143 Nd/ 144 Nd = 0.511869 ± 0.000006 (2σ)) were used as the standard solution in this study.The analytical accuracy of the Sr and Nd isotope data are superior to 0.003%.The 208 Pb/ 206 Pb, 207 Pb/ 206 Pb, and 204 Pb/ 206 Pb ratios of the Pb standard (NBS981) are 2.1681 ± 0.0008 (2σ), 0.91464 ± 0.00033 (2σ), and 0.059042 ± 0.000037(2σ), respectively.The analytical accuracy of the Pb isotope data is better than 0.01%.Specific procedures of the Sr-Nd-Pb isotopic analytical techniques are given by [72].

Zircon U-Pb Dating
Most of zircon grains of sample 19-4s1 are euhedral or subhedral in shape and black in CL imaging.The length of zircons ranges from 80 to 180 µm with length-to-width ratios from 1:1 to 3:1.CL images show that these zircons are with weak internal oscillatory zoning, and/or irregular, patchy to granular internal structures, and growth zonings can be found in some zircons (zircon No. 8, 9, 12, and 15, Figure 6a), indicating that they might not be the typical magmatic zircons [73].Th and U contents of these zircon grains vary from 509 ppm to 1495 ppm (mean = 931 ppm) and 663 ppm to 3129 ppm (mean = 1705 ppm), respectively.The Th/U ratios are variable ranging from 0.18 to 1.73 (mean = 0.66).Several isotopic data of the analyzed zircon grains have relatively big errors, which are eliminated in the process of dating calculation.The 206 Pb/ 238 U ages of thirteen zircons ranges from 216. 4 6b).Most of the zircon grains from sample 24-15s1 are xenotopic, with small amounts of idiomorph, and the aspect ratios are ranging from 1:1 to 4:1 with lengths of 50-200 µm.CL imaging indicates that most of the zircon grains are black in color with weak internal oscillatory zoning, and/or granular internal texture (Figure 6c).These grains show abnormally high contents of U (3277-59,113 ppm; mean = 30,823 ppm), and Th (3347-14,922 ppm; mean = 6389 ppm, Table S1 of Supplementary Materials), which is much higher than the granites without alteration in Xitian pluton with U (354-7047 ppm) and Th (192-1257 ppm) contents [28].Their relatively low Th/U ratios (0.209-0.67, mean = 0.25) indicate a hydrothermal origin [73].The 206 Pb/ 238 U ages of 13 zircon grains from this sample range from 150.6 Ma to 156.2 Ma (Table S1 of Supplementary Materials) and are plotted on or close to the concordia curve, with a weighted mean 206 Pb/ 238 U age of 152.8 ± 1.1 Ma (MSWD = 0.31; Figure 6d).Most of the zircon grains from sample 24-15s1 are xenotopic, with small amounts of idiomorph, and the aspect ratios are ranging from 1:1 to 4:1 with lengths of 50-200 µm.CL imaging indicates that most of the zircon grains are black in color with weak internal oscillatory zoning, and/or granular internal texture (Figure 6c).These grains show abnormally high contents of U (3277-59,113 ppm; mean = 30,823 ppm), and Th (3347-14,922 ppm; mean = 6389 ppm, Table S1 of Supplementary Materials), which is much higher than the granites without alteration in Xitian pluton with U (354-7047 ppm) and Th (192-1257 ppm) contents [28].Their relatively low Th/U ratios (0.209-0.67, mean = 0.25) indicate a hydrothermal origin [73].The 206 Pb/ 238 U ages of 13 zircon grains from this sample range from 150.6 Ma to 156.2 Ma (Table S1 of Materials) and are plotted on or close to the concordia curve, with a weighted mean 206 Pb/ 238 U age of 152.8 ± 1.1 Ma (MSWD = 0.31; Figure 6d).

Trace Element Compositions of Zircons
The zircon grains used for trace element analysis are the same as those which were dated in this study.Their trace element compositions and related parameters are given in Table S2 of Supplementary Materials.
Zircon grains of sample 19-4s1 are characterized by high contents of the heavy rare earth elements (HREEs) and relatively low contents of light rare earth elements (LREEs), with LREE/HREE ratios ranging from 0.04 to 0.24 (mean = 0.1).They also have relatively variable and high contents of the rare earth elements (REEs) with ΣREE ranging from 623 ppm to 2058 ppm (mean = 1071 ppm).Chondrite normalized REE patterns of these zircon grains are characterized by steep slopes, elevated heavy rare earth elements (HREEs), positive Ce anomalies (most of the Ce/Ce* ratios range from 1.61 to 74.14, with average = 10.84), and negative Eu anomalies (Eu/Eu* = 0.17-0.41,mean = 0.28, Figure 7a).

Trace Element Compositions of Zircons
The zircon grains used for trace element analysis are the same as those which were dated in this study.Their trace element compositions and related parameters are given in Table S2 of Supplementary Materials.
Zircon grains of sample 19-4s1 are characterized by high contents of the heavy rare earth elements (HREEs) and relatively low contents of light rare earth elements (LREEs), with LREE/HREE ratios ranging from 0.04 to 0.24 (mean = 0.1).They also have relatively variable and high contents of the rare earth elements (REEs) with ΣREE ranging from 623 ppm to 2058 ppm (mean = 1071 ppm).Chondrite normalized REE patterns of these zircon grains are characterized by steep slopes, elevated heavy rare earth elements (HREEs), positive Ce anomalies (most of the Ce/Ce* ratios range from 1.61 to 74.14, with average = 10.84), and negative Eu anomalies (Eu/Eu* = 0.17-0.41,mean = 0.28, Figure 7a).
Zircon grains from sample 24-15s1 are also characterized by high contents of the heavy rare earth elements (HREEs) and relatively low contents of light rare earth elements (LREEs), with LREE/HREE ratios ranging from 0.01 to 0.13 (mean = 0.02).The REEs contents are tremendously high with ΣREE ranging from 7073 ppm to 39,062 ppm (mean = 15,241 ppm) and Y ranging from 3995 ppm to 43,198 ppm (mean = 18,180 ppm).Chondrite normalized REE patterns that are also characterized by steep slopes, elevated heavy rare earth elements (HREEs), relatively positive Ce anomalies (most of the Ce/Ce* ratios range from = 1.33-8.88,mean = 3.68), and significant negative Eu anomalies (Eu/Eu* = 0.01-0.04,mean = 0.01, Figure 7b).

Muscovite 40 Ar-39 Ar Dating
The Ar-Ar isotopic data of two muscovite samples are given in Table S3 of Supplementary Materials.Thirteen laser-heating stages were carried out for sample 1401-3-1, and stages from 6 to 12 had been obtained the flat age spectrum.This sample yields a plateau age of 156.6 ± 0.7 Ma (MSWD = 0.51, Figure 9a) with 60% release of 39 Ar, consistent with the inverse isochron age of 156.0 ± 0.7 Ma (MSWD = 0.50, Figure 9b).
Seventeen laser-heating stages were carried out for sample 1401-5-3, and stages from 4 to 16 had been obtained the flat age spectrum.This sample yields a well-defined plateau age of 150.0 ± 0.6 Ma (MSWD = 0.39 Ma, Figure 9c) with 97% release of 39 Ar, consistent with the inverse isochron age of 149.5 ± 0.8 Ma (MSWD = 0.17, Figure 9d).

Muscovite 40 Ar-39 Ar Dating
The Ar-Ar isotopic data of two muscovite samples are given in Table S3 of Supplementary Materials.Thirteen laser-heating stages were carried out for sample 1401-3-1, and stages from 6 to 12 had been obtained the flat age spectrum.This sample yields a plateau age of 156.6 ± 0.7 Ma (MSWD = 0.51, Figure 9a) with 60% release of 39 Ar, consistent with the inverse isochron age of 156.0 ± 0.7 Ma (MSWD = 0.50, Figure 9b).
Seventeen laser-heating stages were carried out for sample 1401-5-3, and stages from 4 to 16 had been obtained the flat age spectrum.This sample yields a well-defined plateau age of 150.0 ± 0.6 Ma (MSWD = 0.39 Ma, Figure 9c) with 97% release of 39 Ar, consistent with the inverse isochron age of 149.5 ± 0.8 Ma (MSWD = 0.17, Figure 9d).
The trace and rare earth elements (REEs) contents of skarns are given in Table S4 of Supplementary Materials.The total REE (ΣREE) contents range from 40.4 to 184.5 ppm (mean = 103.2ppm), with high (La/Yb)N ratios (5.2-13.2),and significant negative Eu anomalies (Eu/Eu* = 0.36-0.82).The REE distribution patterns display right-dipping V-type curves, which are exactly similar to those of the Indosinian granites in Xitian (Figure 10a).These samples also show similar primitivemantle-normalized trace element patterns to those of Indosinian granites in Xitian, characterized by enrichment in Rb, K, U, Zr, Hf and REE, and depletion in Ti, P, Sr, Ba, and Nb, which are also similar to the Indosinian granites in the Xitian ore field (Figure 10b).
The trace and rare earth elements (REEs) contents of skarns are given in Table S4 of Supplementary Materials.The total REE (ΣREE) contents range from 40.4 to 184.5 ppm (mean = 103.2ppm), with high (La/Yb) N ratios (5.2-13.2),and significant negative Eu anomalies (Eu/Eu* = 0.36-0.82).The REE distribution patterns display right-dipping V-type curves, which are exactly similar to those of the Indosinian granites in Xitian (Figure 10a).These samples also show similar primitive-mantle-normalized trace element patterns to those of Indosinian granites in Xitian, characterized by enrichment in Rb, K, U, Zr, Hf and REE, and depletion in Ti, P, Sr, Ba, and Nb, which are also similar to the Indosinian granites in the Xitian ore field (Figure 10b).ppm), with high (La/Yb)N ratios (5.2-13.2),and significant negative Eu anomalies (Eu/Eu* = 0.36-0.82).The REE distribution patterns display right-dipping V-type curves, which are exactly similar to those of the Indosinian granites in Xitian (Figure 10a).These samples also show similar primitivemantle-normalized trace element patterns to those of Indosinian granites in Xitian, characterized by enrichment in Rb, K, U, Zr, Hf and REE, and depletion in Ti, P, Sr, Ba, and Nb, which are also similar to the Indosinian granites in the Xitian ore field (Figure 10b).Normalizing factors are from [76].Data for the Xitian Indosinian and Yanshanian granites are from [28,38].[76].Data for the Xitian Indosinian and Yanshanian granites are from [28,38].

Wolframite Sr-Nd-Pb Isotopic Composition Analysis
The wolframite Sr-Nd isotopic compositions of the four samples analyzed in this study are given in Table S5 of Supplementary Materials.The initial 87 Sr/ 86 Sr and 143 Nd/ 144 Nd ratios were calculated using the muscovite 40 Ar-39 Ar dating result of 150.0 ± 0.6 Ma.These samples have high initial ( 87 Sr/ 86 Sr) i ratios (0.71282-0.72003),low initial ( 143 Nd/ 144 Nd) i ratios (0.511644-0.512155),negative ε Nd (t) values ranging from −15.6 to −5.6, and old model ages (T DM2 ) of 2168-1396 Ma.
The Pb isotopic compositions of the samples are shown in Table S5 of Supplementary Materials.All samples are enriched in radiogenic Pb, with 206 Pb/ 204 Pb ratios of 18.489-18.569, 207Pb/ 204 Pb ratios of 15.724-15.877,and 208 Pb/ 204 Pb ratios of 39.055-39.335,respectively.The Pb isotopic values were calculated with the muscovite 40 Ar-39 Ar dating result of 150.0 ± 0.6 Ma using single-stage Pb isotopic evolution model [77], with values of ∆β ranging from 26.09 to 36.06, ∆γ values from 48.95 to 58.62, and µ values from 9.81 to 9.99, respectively.

Timing of Mineralization and Granitic Magmatism
In order to constrain the time interval of mineralization in the Xitian W-Sn deposit, several studies have been carried out using various dating technologies [19, [40][41][42]78].It was first reported the muscovite 40 Ar-39 Ar isotopic ages of ore-bearing quartz vein in skarn and greisen in the Longshang deposit, with ages of 155.6 ± 1.3 Ma and 157.2 ± 1.4 Ma, respectively [41].Then, some authors obtained the metallogenic age of other deposits in the Xitian ore field, such as Heshuxia (molybdenite Re-Os age of 150.0 ± 2.7 Ma; [78]), Shantian (molybdenite Re-Os age of 158.9 ± 2.2 Ma; [40]), and Hejiangkou (molybdenite Re-Os age of 225.5 ± 3.6 Ma; [43]).An overwhelming majority of dating minerals are collected from ore-bearing quartz vein and/or greisen with ages ranging from 159 Ma to 149 Ma which has a congruent relationship with the early Yanshanian granitic magmatism [19, [40][41][42]78].It seems that the quartz vein and/or greisen type ore bodies in the Xitian W-Sn deposit is close to the early Yanshanian magmatic hydrothermal activities. Itwas reported a molybdenite Re-Os age of 225.5 ± 3.6 Ma of altered granite in the Hejiangkou deposit, regarded as the proof for the Indosinian mineralization in Xitian ore field [43].Here, we reported a muscovite 40 Ar-39 Ar isotopic age of the quartz vein type ore body (156.0 ± 0.7 Ma) in the Hejinagkou deposit of the Xitian ore field, indicating that the quartz vein type W-Sn mineralization in this deposit is closely linked to Yanshanian granitic magmatic activities.However, there is still no direct chronologic evidence to seek an answer for the skarn-type mineralization and restricting the genesis of the skarn-type mineralization in this area.We reported the muscovite 40 Ar-39 Ar isotopic ages of ore-bearing quartz vein for the Hejiangkou deposit and Heshuxia deposit, and two hydrothermal zircon U-Pb ages for the altered granites are 156.6 ± 0.7 Ma, 149.5 ± 0.8 Ma, 152.8 ± 1.1 Ma, and 226.0 ± 2.8 Ma, respectively.Together with the ages published and obtained in this study, we draw a conclusion that there two phases of metallogenic events have been recorded with the time interval of the quartz vein and greisen type W-Sn mineralization ranging from 159 to 149 Ma and skarn-type W-Sn mineralization ca.225 Ma in the Xitian ore field, respectively.
The Yanshanian period is a significant time interval of granitic magmatism and W-Sn metallogenesis in the Nanling region, which was regarded as the period of mineralization explosion in South China [79].However, compared with the Yanshanian W-Sn mineralization event in the Nanling region, the Indosinian W-Sn mineralization event is inconspicuous (Table S6 of Supplementary Materials).Although only a few deposits were determined to be the products of Indosinian magmatic hydrothermal events, such as Shuiyuanshan (220.6 ± 1.1 Ma; [31]), Yeziwo (227.2 ± 1.5; [31]), Hehuaping (224.0 ± 1.9 Ma; [32]), and Limu (214.1 ± 1.9 Ma; [33]), it indicates that the Indosinian is also an important epoch for the W-Sn mineralization in South China.

Physico-Chemical Conditions of the Ore-Forming Fluids
Recently, several studies proved that chemical compositions of zircons could be used as a valid tracer to reflect the physico-chemical conditions of the magmatic melt and/or hydrothermal fluid [80][81][82][83][84]. Specifically, the positive Ce anomalies in zircon are the result of the oxidation of Ce 3+ to Ce 4+ .Ce 4+ is compatible in zircon and can substitute for Zr 4+ , Hf 4+ and other tetravalent elements in zircon lattice.Thus, Ce anomalies can provide information for the oxidation state of magma and/or related fluid, where higher Ce/Ce* ratios are in accordance with the high oxygen fugacity (f O 2 ; [83][84][85].However, zircon Eu anomalies are not the efficient tracer, because Eu anomalies in zircon are controlled not only by the redox state of the fluids but also by the crystallization history of plagioclase [86,87].Due to the relatively low ratios of Ce/Ce* for the zircons from sample 19-4s1 and 24-15s1, almost all the zircon grains plot in the field of low f O 2 indicating that these formed in a lower f O 2 environment (Figure 11a).Using the model proposed by [83], we also estimate the redox conditions for the samples 19-4s1 and 24-15s1 (Table S2 of

Physico-Chemical Conditions of the Ore-Forming Fluids
Recently, several studies proved that chemical compositions of zircons could be used as a valid tracer to reflect the physico-chemical conditions of the magmatic melt and/or hydrothermal fluid [80][81][82][83][84]. Specifically, the positive Ce anomalies in zircon are the result of the oxidation of Ce 3+ to Ce 4+ .Ce 4+ is compatible in zircon and can substitute for Zr 4+ , Hf 4+ and other tetravalent elements in zircon lattice.Thus, Ce anomalies can provide information for the oxidation state of magma and/or related fluid, where higher Ce/Ce* ratios are in accordance with the high oxygen fugacity (fO2; [83][84][85].However, zircon Eu anomalies are not the efficient tracer, because Eu anomalies in zircon are controlled not only by the redox state of the fluids but also by the crystallization history of plagioclase [86,87].Due to the relatively low ratios of Ce/Ce* for the zircons from sample 19-4s1 and 24-15s1, almost all the zircon grains plot in the field of low fO2 indicating that these formed in a lower fO2 environment (Figure 11a).Using the model proposed by [83], we also estimate the redox conditions for the samples 19-4s1 and 24-15s1 (Table S2 of Supplementary Materials), yielding the relatively low oxygen fugacities with log fO2 values ranging from −19 to −15 (mean = −17) and from −19 to −13 (mean = −15), respectively.This data further confirmed that both the Indosinian and Yanshanian hydrothermal fluids are reducing fluids which are favorable for the W-Sn mineralization [88][89][90][91].
As an efficient tracer, Y/Ho ratios can provide evidence about the chemical characteristics of source fluids [81].Fractionation between the Y and Ho occurs in highly evolved granitic melts or hydrothermally altered granites with participation of F-rich fluids which contain high concentrations of Y, Li, B, and/or P [81,82,92,93].Zircon grains of the sample 19-4s1 and 24-15s1 are ranging from 28 to 34 (mean = 31) and from 29 to 31 (mean = 31), respectively, which are higher than the chondritic value of 28 (Figure 11a; [94]).The high ratios of zircons from these two samples suggest that they were crystallized in F-rich fluids, which are consistent with the existence of abundant fluorite ore bodies (Figure 11b).

Source of Ore-Forming Metals
On account of the extremely similar geochemical behavior, the REEs are always involved in the geological process in group, and the hydrothermal metamorphism will not change the composition mode and distribution mode of REEs in minerals or rocks, making them efficient tracers for determining the source rocks and element migration mechanism in ore-forming processes [95][96][97][98].According to the chondrite-normalized REE patterns of Xitian skarns, they all show a good consistency, exhibiting right-dipping V-type curves with obvious negative Eu anomalies.Furthermore, the REE patterns of Xitian skarns are extremely similar to those of Xitian Indosinian granites, which are distinctly different to Yanshanian granites (Figure 10b; [28,38]).It indicates that the origin of skarn has a genetic relationship with Indosinian granites rather than Yanshanian granites.The Indosinian granitoids of the Xitian pluton are high-Si, high-K, weakly to strongly peraluminous, and highly fractionated S-type granites with high initial 87 Sr/ 86 Sr isotope ratios (0.71397-0.71910),low εNd(t) values ranging from −10.1 to −9.4, and old Nd model ages (1858-1764 Ma), indicating that the Xitian Indosinian granites were mainly originated from partial melting of Paleoproterozoic metamorphic basement with small amounts of mantle-derived magma involved [38].Since the Indosinian granites are closely related to the skarn type W-Sn mineralization, it can be inferred that the source of ore-forming metals from skarn type ore bodies should be mainly originated from a crustal source.
Radiogenic Sr, Nd, and Pb isotopes are powerful tools not only to determine magma sources, but also to determine ore-forming metals in minerals [99][100][101][102][103].As the most important tungstenbearing mineral in the tungsten deposit, geochemical compositions and mineralogical features of wolframite can provide abundant metallogenetic information, i.e., fluid, environment, and metal source [1,5,44,104].The samples of wolframite analyzed in this study have high initial 87 Sr/ 86 Sr ratios (0.71282-0.72003),indicating that the wolframites in quartz vein type ore bodies are characterized by a crustal source.Furthermore, these samples also have negative εNd(t) values ranging from −15.9 to −5.6 and cover a wide range of model ages (TDM2 = 2168 − 1399 Ma).In addition, almost all the samples are plotted near the evolution line of the Upper Continental Crust (UCC), and two samples are plotted in the field of the Meso-Paleoproterozoic low mature basement (Figure 12a), indicating that the ore-forming metals are mainly originated from the crustal source, which are proved by the diagram of age versus εNd(t) (Figure 12b).Lead isotopic compositions of wolframites from the Xitian deposit are characterized by high radiogenic Pb isotope values with 206 Pb/ 204 Pb ratios of 18.489-18.569,[82].Symbols are the same as those in Figure 7.
As an efficient tracer, Y/Ho ratios can provide evidence about the chemical characteristics of source fluids [81].Fractionation between the Y and Ho occurs in highly evolved granitic melts or hydrothermally altered granites with participation of F-rich fluids which contain high concentrations of Y, Li, B, and/or P [81,82,92,93].Zircon grains of the sample 19-4s1 and 24-15s1 are ranging from 28 to 34 (mean = 31) and from 29 to 31 (mean = 31), respectively, which are higher than the chondritic value of 28 (Figure 11a; [94]).The high ratios of zircons from these two samples suggest that they were crystallized in F-rich fluids, which are consistent with the existence of abundant fluorite ore bodies (Figure 11b).

Source of Ore-Forming Metals
On account of the extremely similar geochemical behavior, the REEs are always involved in the geological process in group, and the hydrothermal metamorphism will not change the composition mode and distribution mode of REEs in minerals or rocks, making them efficient tracers for determining the source rocks and element migration mechanism in ore-forming processes [95][96][97][98].According to the chondrite-normalized REE patterns of Xitian skarns, they all show a good consistency, exhibiting right-dipping V-type curves with obvious negative Eu anomalies.Furthermore, the REE patterns of Xitian skarns are extremely similar to those of Xitian Indosinian granites, which are distinctly different to Yanshanian granites (Figure 10b; [28,38]).It indicates that the origin of skarn has a genetic relationship with Indosinian granites rather than Yanshanian granites.The Indosinian granitoids of the Xitian pluton are high-Si, high-K, weakly to strongly peraluminous, and highly fractionated S-type granites with high initial 87 Sr/ 86 Sr isotope ratios (0.71397-0.71910),low ε Nd (t) values ranging from −10.1 to −9.4, and old Nd model ages (1858-1764 Ma), indicating that the Xitian Indosinian granites were mainly originated from partial melting of Paleoproterozoic metamorphic basement with small amounts of mantle-derived magma involved [38].Since the Indosinian granites are closely related to the skarn type W-Sn mineralization, it can be inferred that the source of ore-forming metals from skarn type ore bodies should be mainly originated from a crustal source.
Radiogenic Sr, Nd, and Pb isotopes are powerful tools not only to determine magma sources, but also to determine ore-forming metals in minerals [99][100][101][102][103].As the most important tungsten-bearing mineral in the tungsten deposit, geochemical compositions and mineralogical features of wolframite can provide abundant metallogenetic information, i.e., fluid, environment, and metal source [1,5,44,104].The samples of wolframite analyzed in this study have high initial 87 Sr/ 86 Sr ratios (0.71282-0.72003),indicating that the wolframites in quartz vein type ore bodies are characterized by a crustal source.Furthermore, these samples also have negative ε Nd (t) values ranging from −15.9 to −5.6 and cover a wide range of model ages (T DM2 = 2168 − 1399 Ma).In addition, almost all the samples are plotted near the evolution line of the Upper Continental Crust (UCC), and two samples are plotted in the field of the Meso-Paleoproterozoic low mature basement (Figure 12a), indicating that the ore-forming metals are mainly originated from the crustal source, which are proved by the diagram of age versus ε Nd (t) (Figure 12b).Lead isotopic compositions of wolframites from the Xitian deposit are characterized by high radiogenic Pb isotope values with 206 Pb/ 204 Pb ratios of 18.489-18.569, 207Pb/ 204 Pb ratios of 15.724-15.877,and 208 Pb/ 204 Pb ratios of 39.055-39.335,and high values of µ ranging from 9.70 to 9.99, indicating that Pb was derived from the upper crust.Additionally, the Pb isotopic compositions of the wolframites are similar to the sulfides of the quartz vein type ore bodies, and almost all the samples are plotted primarily toward the upper crust evolution field and/or line in the ∆β versus ∆γ and 206 Pb/ 204 Pb versus 207 Pb/ 204 Pb diagrams (Figure 12c,d).Briefly, based on the diagenetic and metallogenic geochronology, the relationship between the metallogenic and magmatic activities, and the evidence of mineral isotopes, we deduce that the ore-forming metals of quartz vein/greisen type ore bodies were also derived from a crustal source.[59,78,105,106].Data of sulfides are from [107].MORB: Mid Ocean Ridge Basalt; DMM: depleted end-member; CHUR: chondritic uniform reservoir.

Conclusions
1. Hydrothermal zircon U-Pb and muscovite 40 Ar/ 39 Ar dating suggests that there are two epochs of W-Sn mineralization in the Xitian ore field, with skarn-type W-Sn mineralization at ca. 226 Ma and quartz vein/greisen type W-Sn mineralization at ca. 156 Ma.
2. The ore-forming fluids for the two metallogenic events are both characterized by enrichment in F and low oxygen fugacities.
3. The ore-forming metals for the skarn-type and quartz vein/greisen type W-Sn mineralization

Figure 2 .
Figure 2. Schematic geological map of the Xitian W-Sn ore field showing the location of samples (modified from [44]).

Figure 2 .
Figure 2. Schematic geological map of the Xitian W-Sn ore field showing the location of samples (modified from [44]).

Figure 3 .
Figure 3. Photographs of the related ore bodies showing the morphology of skarn-and quartz veintype ore body.(a) The stratiform and stratoid ore bodies occurring in the contact zone between the Indosinian granites and the Devonian dolomitic limestone in Longshang ore district; (b) Quartz vein type of ore body in Longshang ore district; (c) Skarn-type ore bodies are superimposed by calcite-and fluorite-veins in Longshang ore district; (d) Quartz vein type of ore body in Hejiangkou ore district; (e) Quartz vein type of ore body in Heshuxia ore district; (f) Quartz vein type of ore body in Heshuxia ore district; (g) Quartz vein type of ore body in Goudalan ore district; (h) The structurally altered rock-type of Pb-Zn ore body in Shaiheling ore district; (i) The stratiform and stratoid ore bodies occurring in the contact zone between the Indosinian granites and the Devonian dolomitic limestone in Shaiheling ore district; Cal: calcite; Fl: fluorite; Gn: galenite; Mo: molybdenite; Ms: muscovite; Wol: wolframite.

Figure 3 .
Figure 3. Photographs of the related ore bodies showing the morphology of skarn-and quartz vein-type ore body.(a) The stratiform and stratoid ore bodies occurring in the contact zone between the Indosinian granites and the Devonian dolomitic limestone in Longshang ore district; (b) Quartz vein type of ore body in Longshang ore district; (c) Skarn-type ore bodies are superimposed by calcite-and fluorite-veins in Longshang ore district; (d) Quartz vein type of ore body in Hejiangkou ore district; (e) Quartz vein type of ore body in Heshuxia ore district; (f) Quartz vein type of ore body in Heshuxia ore district; (g) Quartz vein type of ore body in Goudalan ore district; (h) The structurally altered rock-type of Pb-Zn ore body in Shaiheling ore district; (i) The stratiform and stratoid ore bodies occurring in the contact zone between the Indosinian granites and the Devonian dolomitic limestone in Shaiheling ore district; Cal: calcite; Fl: fluorite; Gn: galenite; Mo: molybdenite; Ms: muscovite; Wol: wolframite.
Minerals 2018, 8, x FOR PEER REVIEW 8 of 23 contained sulfides, collected from the Longshang, Shaiheling and Huamu deposits.The wolframite, separated from the quartz vein type ore and used for Sr-Nd-Pb analysis, were collected from Longshang, Goudalan and Heshuxia.The sampling locations were marked in Figure 2.

Figure 5 .
Figure 5. Hand specimen and micrographs of the altered granite.(a) Hand specimen of the altered Indosinian coarse-grained biotite granite; (b) The K-feldspar is altered and replaced by sericite; (c) Hand specimen of the altered Yanshanian fine-grained two-mica granite; (d) Almost the feldspars are altered and replaced by sericite.Kfs: K-feldspar; Pl: plagioclase; Ser: sericite.

Figure 5 .
Figure 5. Hand specimen and micrographs of the altered granite.(a) Hand specimen of the altered Indosinian coarse-grained biotite granite; (b) The K-feldspar is altered and replaced by sericite; (c) Hand specimen of the altered Yanshanian fine-grained two-mica granite; (d) Almost the feldspars are altered and replaced by sericite.Kfs: K-feldspar; Pl: plagioclase; Ser: sericite.

Figure 6 .
Figure 6.Cathodoluminescence (CL) images of representative zircon grains and concordia diagrams of zircon U-Pb geochronological data for the samples taken from the Xitian ore field.(a) Cathodoluminescence (CL) images of the zircons from altered Indosinian granites; (b) Concordia diagram of zircon U-Pb data for the altered Indosinian granites; (c) Cathodoluminescence (CL) images of the zircons from altered Yanshanian granites; (b) Concordia diagram of zircon U-Pb data for the altered Yanshanian granites.

Figure 6 .
Figure 6.Cathodoluminescence (CL) images of representative zircon grains and concordia diagrams of zircon U-Pb geochronological data for the samples taken from the Xitian ore field.(a) Cathodoluminescence (CL) images of the zircons from altered Indosinian granites; (b) Concordia diagram of zircon U-Pb data for the altered Indosinian granites; (c) Cathodoluminescence (CL) images of the zircons from altered Yanshanian granites; (d) Concordia diagram of zircon U-Pb data for the altered Yanshanian granites.

Figure 7 .
Figure 7. Chondrite-normalized REE chemistry of zircon grains for the samples taken from the Xitian ore field, with normalizing factors from [74].(a) Chondrite-normalized REE chemistry of zircon grains for altered Indosinian granites; (b) Chondrite-normalized REE chemistry of zircon grains for altered Yanshanian granites.

Figure 7 . 23 Figure 8 .
Figure 7. Chondrite-normalized REE chemistry of zircon grains for the samples taken from the Xitian ore field, with normalizing factors from [74].(a) Chondrite-normalized REE chemistry of zircon grains for altered Indosinian granites; (b) Chondrite-normalized REE chemistry of zircon grains for altered Yanshanian granites.
Supplementary Materials), yielding the relatively low oxygen fugacities with log f O 2 values ranging from −19 to −15 (mean = −17) and from −19 to −13 (mean = −15), respectively.This data further confirmed that both the Indosinian and Yanshanian hydrothermal fluids are reducing fluids which are favorable for the W-Sn mineralization [88-91].Minerals 2018, 8, x FOR PEER REVIEW 15 of 23

Figure 11 .
Figure 11.(a) Ce/Ce* versus Eu/Eu*; and (b) Y/Ho versus Y plots of zircon grains from the Xitian ore field.(b) is modified from[82].Symbols are the same as those in Figure7.

Figure 11 .
Figure 11.(a) Ce/Ce* versus Eu/Eu*; and (b) Y/Ho versus Y plots of zircon grains from the Xitian ore field.(b) is modified from[82].Symbols are the same as those in Figure7.
5a to 233.5Ma which plot on or near the concordia curve (TableS1of Supplementary Materials), yielding a weighted mean 206 Pb/ 238 U age of 226.0 ± 2.8 Ma (MSWD = 2.1, Figure6b).Several isotopic data of the analyzed zircon grains have relatively big errors, which are eliminated in the process of dating calculation.The 206 Pb/ 238 U ages of thirteen zircons ranges from 216.4 Ma to 233.5 Ma which plot on or near the concordia curve (TableS1of Supplementary Materials), yielding a weighted mean 206 Pb/ 238 U age of 226.0 ± 2.8 Ma (MSWD = 2.1, Figure (mean = 0.66).