Texture and Geochemistry of Scheelites in the Tongshankou Deposit in Daye, Hubei, China: Implication for REE Substitution Mechanism and Multistage W Mineralization Processes

: The Tongshankou skarn deposit in the Edong ore district is a typical metasomatic deposit associated with adjacent granodiorite porphyry and carbonate rocks. Using comprehensive micro-scopic observations, mineralogical and geochemical analysis, scheelite grains in the skarns can be classiﬁed into three major types, showing multi-stage mineralization characteristics. In the redox ﬂuid environment, scheelites that occur with garnets usually have afﬁnity to garnets, while in later skarn phases others exist with oxides and sulﬁdes. They can be subdivided by trace elements, such as Nb and Eu, to discuss the detailed ore-forming process. Scheelites have three typical substitution mechanisms including: 2Ca 2 + (cid:10) REE 3 + +Na + (1); Ca 2 + + W 6 + (cid:10) REE 3 + +Nb 5+ (2); and 3Ca 2 + (cid:10) 2REE 3+ + (cid:3) Ca ( (cid:3) Ca = Ca site vacancy) (3). Plagioclase or various hydrothermal stages can cause Eu anomalies to ﬂuctuate from positive to negative, and these processes can cause particular zonation in W and Mo contents in scheelites. This study highlights the use of texture and geochemistry of scheelites in skarn deposits, depicting the W mineralization processes. Contributions: Conceptualization, R.-Z.Z. and M.-F.W.; methodology, R.-Z.Z. and J.-P.W.; software, R.-Z.Z. and X.-Y.S.; validation, M.-F.W., H.L.; formal analysis, R.-Z.Z. and X.-Y.S.; investigation, R.-Z.Z., M.-F.W. and X.-Y.S.; resources, M.-F.W. and H.L.; data curation, M.-F.W., H.L. and Z.U.; writing—original draft preparation, R.-Z.Z. and X.-Y.S.; writing—review and editing, R.-Z.Z., M.-F.W., H.L. and Z.U.; visualization, M.-F.W., J.-P.W., and H.L.; supervision, M.-F.W., H.L. and Z.U.; funding acquisition, M.-F.W. All authors read and agreed the published version of the manuscript.


Introduction
Tungsten has been listed as one of the most important industrial metals in alloy manufacturing, with a wide range of applications, and has also been listed as one of the critical metals by nearly all developed countries due to its strategic importance [1]. The Edong ore district, located in the Middle-Lower Yangtze River belt (MLYRB), is a multi-metal metallogenic belt containing mostly porphyry and skarn type deposits. The ore district includes the Tongshankou porphyry-skarn Cu-Mo-(W) deposit, and has a comparable geological environment to the adjacent Ruanjiawan and Fujiashan-Longjiaoshan W deposits [2][3][4]. Scheelite is the main tungsten mineral in both the porphyry and skarn systems in the Tongshankou area. There has been, however, no comprehensive study on the genesis of scheelites and the W transportation mechanism in the Tongshankou deposits. A cold cathodoluminescence study revealed pores and inner texture in scheelite, indicating multi-stage ore formation [5]. Tungsten is mostly concentrated in mafic enclaves compared to diorite intrusions and limestone stratum, which may have originated from a different magma source due to mantle upwelling [6]. The focus of this research is to describe the pared to diorite intrusions and limestone stratum, which may have originated from a different magma source due to mantle upwelling [6]. The focus of this research is to describe the field investigation, geological relationship, textural mineralogy and geochemistry of scheelites in a skarn system in the Tongshankou deposit, for better understanding of their REE substitution mechanisms and formation processes of scheelite. In addition, this research provides insight into the ore prospecting based on the ore-forming processes and geochemical characteristics. The purpose of this research is to investigate the tungsten metallogenic process in order to illustrate the genesis of scheelite mineralization in Tongshankou deposit.

Regional and Deposit Geology
The MLYRB is one of the most important polymetallic mineralization belts of the Yangtze plate. The Lower Yangtze depression is dominated by sedimentary cap rocks in the north of Yangtze plate, and Dabie orogenic belt in the south [7,8] (Figure 1a). The region is bounded by the Xiangfan-Guangji fault to the northwest, the Tancheng-Lujiang regional strike-slip fault to the northeast and the Yangxin-Changzhou fault to the south [9] (Figure 1b). From Echeng to Yinzu intrusions, a significant number of north-dipping thrust faults are linked together, forming complicated tectonics, thrust folds that controlled the majority of mineralized structures in this area, including Tongshankou and Daye [10]. The Yanshanian and Himalayan orogeny in Mesozoic time caused the main magmatic intrusions in the Edong district, and played the major role of orogeny, which is directly connected with Yinzu and Lingxiang intrusions. Echeng, Tieshan, Jinshandian, Yinzu, Lingxiang and Yangxin are the six main magmatic intrusions exposed in the close vicinity of the study area ( Figure 2). The Daye syncline and Yinzu anticline, which exist between the Yinzu and Lingxiang intrusions, are key structures that influenced the metallogeny of the Tongshankou region. At 143-145 Ma, the early formed porphyry system was an oxidation type Cu (Mo) deposit [11] that played the role of intermediate-acid intrusion rocks in skarn mineralization. As mentioned, complicated tectonics within this area included seven composite folds: Echeng anticline, Huajiahu syncline, Tieshan anticline, Huangjinshan syncline, Baoan-Wangren anticline, Daye syncline and Yinzu anticline, most of which were formed during early Yanshanian. There are over 200 faults within the Edong district with different trending that can be classified relatively into 5 types: NWW, NE, NNE, NW and E-W [9][10][11][12].  [9][10][11][12][13]. Abbreviations: TLF-Tancheng-Lujiang Fault, XGF-Xiangfan-Guangji Fault, YCF-Changzhou-Yangxin Fault; (b) detail location of the Edong ore district and simple illustration of surrounding geology environment (modified after [9,12,13]). Changzhou-Yangxin Fault; (b) detail location of the Edong ore district and simple illustration of surrounding geology environment (modified after [9,12,13]).

Figure 2.
The Edong ore district geologic map, illustrating the main magmatic intrusions in this area (modified after [13]).
The nearby Early Triassic carbonate rocks strata offered a suitable geochemical environment for Ca-skarn formation. K-feldspar, biotite, garnet, apatite, wollastonite, diopside, chlorite, hornblende, quartz, scheelite, pyrite, molybdenite, chalcopyrite, calcite and galena are the most common skarn minerals. Scheelite is a common mineral associated with garnet, oxides, sulfides and quartz in skarn deposits [14]. There is a zonation of intensive skarn alteration, weak skarn alteration and mild chlorite alteration within the deposit district from intrusions to the carbonate strata [15].
According to previous Sr-Nd and U-Pb zircon studies, the Tongshankou deposits were formed around 145 Ma and are linked to skarn mineralization [12]. The Tongshankou skarn deposit is a typical metasomatic deposit with layers of granodiorite porphyry and carbonate rocks. The Edong ore district, which is part of the largest MLYRB deposits, has multi-metal resources such as Cu, Fe, Ni, Au, W, Mo, Mn, Pb, Zn and Co. The Tongshankou deposit has an association of economically Cu-Mo-W (Au) ore bodies that have been prospected, mostly in the porphyry portions that are localized towards endoskarn rocks [7,8] (Figure 3). The nearby Early Triassic carbonate rocks strata offered a suitable geochemical environment for Ca-skarn formation. K-feldspar, biotite, garnet, apatite, wollastonite, diopside, chlorite, hornblende, quartz, scheelite, pyrite, molybdenite, chalcopyrite, calcite and galena are the most common skarn minerals. Scheelite is a common mineral associated with garnet, oxides, sulfides and quartz in skarn deposits [14]. There is a zonation of intensive skarn alteration, weak skarn alteration and mild chlorite alteration within the deposit district from intrusions to the carbonate strata [15].
According to previous Sr-Nd and U-Pb zircon studies, the Tongshankou deposits were formed around 145 Ma and are linked to skarn mineralization [12]. The Tongshankou skarn deposit is a typical metasomatic deposit with layers of granodiorite porphyry and carbonate rocks. The Edong ore district, which is part of the largest MLYRB deposits, has multi-metal resources such as Cu, Fe, Ni, Au, W, Mo, Mn, Pb, Zn and Co. The Tongshankou deposit has an association of economically Cu-Mo-W (Au) ore bodies that have been prospected, mostly in the porphyry portions that are localized towards endoskarn rocks [7,8] (Figure 3).

Sample Description
Representative rock and ore samples were collected from two drill holes in the Tongshankou deposit to determine mineral composition and scheelite trace element geochemistry. Exoskarn samples were collected from molybdenite and scheelite-garnet zone (ZK2002-1) at a depth of 760 m, and from the ore-bearing quartz-sulfides zone (ZK2002-2) at a depth of 1250 m, while the endoskarn samples of diopside-garnet (ZKB20S) were taken at a depth of 650-700 m.
Scheelites in the Tongshankou deposit are widely distributed in both endoskarn and exoskarn, and are associated with garnet. Three types of scheelite (SchA, B and C) are identified under the microscope on the basis of its paragenesis. Epidote, chlorite, actinolite and serpentine are the most common retrograde minerals here, which formed after skarn stage. SchA is subhedral to anhedral scheelite within the ZK2002 exoskarn sample. Some grains of SchA are associated with sulfides, which coexist in the highly fractured and metasomatized garnet (Figure 4a,b); SchB is granular scheelite within both ZK2002 and ZKB20S samples (Figure 4c,d); while SchC mostly occurring in ZKB20S sample is subhedral to anhedral, associated with calcite and clustered stockwork of quartz grains ( Figure  4e-h). Therefore, the SchA and SchB scheelites are distributed in exoskarn deposits, with skarn stage minerals and SchC scheelites identified in endoskarn deposit and associated with calcite filled by the later stage minerals.

Sample Description
Representative rock and ore samples were collected from two drill holes in the Tongshankou deposit to determine mineral composition and scheelite trace element geochemistry. Exoskarn samples were collected from molybdenite and scheelite-garnet zone (ZK2002-1) at a depth of 760 m, and from the ore-bearing quartz-sulfides zone (ZK2002-2) at a depth of 1250 m, while the endoskarn samples of diopside-garnet (ZKB20S) were taken at a depth of 650-700 m.
Scheelites in the Tongshankou deposit are widely distributed in both endoskarn and exoskarn, and are associated with garnet. Three types of scheelite (SchA, B and C) are identified under the microscope on the basis of its paragenesis. Epidote, chlorite, actinolite and serpentine are the most common retrograde minerals here, which formed after skarn stage. SchA is subhedral to anhedral scheelite within the ZK2002 exoskarn sample. Some grains of SchA are associated with sulfides, which coexist in the highly fractured and metasomatized garnet (Figure 4a,b); SchB is granular scheelite within both ZK2002 and ZKB20S samples (Figure 4c,d); while SchC mostly occurring in ZKB20S sample is subhedral to anhedral, associated with calcite and clustered stockwork of quartz grains (Figure 4e-h). Therefore, the SchA and SchB scheelites are distributed in exoskarn deposits, with skarn stage minerals and SchC scheelites identified in endoskarn deposit and associated with calcite filled by the later stage minerals.

Microprobe Analysis
Major element compositions and elemental mapping were analyzed by electron microprobe analysis (EMPA) at the State Key Laboratory of Geological Processes and Mineral Resources (GPMR), China University of Geosciences (Wuhan), with a JEOL JXA-8100 EMPA equipped with four wavelength-dispersive spectrometers (WDS). An accelerating voltage of 15 kV, beam current of 20 nA and <5 µm-diameter focused beam was used in the analyses.
The raw data were corrected online by using the modified ZAF (atomic number, absorption, fluorescence) correction procedure. Element peaks and backgrounds were measured for all the elements with counting times of 10 s and 5 s, respectively, except for

Optical Microscopy Cathodoluminescence (OM-CL)
Polished thin sections of representative samples were examined by OM-CL, using a Leica DM2700P microscope coupled with a CITLMK 5-2 system, at GPMR, China University of Geosciences (Wuhan). The system was operated at 13 kV accelerating voltage and a current density of about 350 цA.

Scanning Electron Microscopy Cathodoluminescence (SEM-CL)
The BSE and CL imaging were performed on polished and carbon-coated samples to characterize the morphology and internal texture of scheelite. The imaging was conducted with a Quanta 450 FEG SEM, equipped with a SDD Inca X-Max 50 and a MonoCL 4+ detector (20 kV accelerating voltage and 20 nA beam current) at GPMR, China University of Geosciences (Wuhan).

LA-ICP-MS Trace Element Analysis
Trace elemental compositions of the Tongshankou scheelite (SchA to SchC) were determined at the National Research Center for Geo-analysis (NRCG), the Chinese Academy of Geological Sciences, Beijing, China by using in situ laser ablation inductively coupled plasma-mass spectrometry (LA-ICP-MS), followed by New Wave NWR-193 laser ablation system coupled to a Thermo Element II mass spectrometer. Prior to entering the ICP, helium (carrying gas) and argon (makeup gas) were combined through a T-connector. Helium (carrier gas) and argon (makeup gas) were mixed via a T-connector prior to the entering of the ICP. The analyzing parameters of LA-ICP-MS demonstrated an accelerating frequency of 10 Hz, beam energy of~2 mJ and a beam density of 23-25 J/cm 2 with 35 µm laser spot size. Each analysis comprises a 15 s gas blank measurement (laser-off), followed by a 40 s sample data acquisition (laser-on), and a 15 s cleaning between successive measurements. The ablated materials were flushed in a continuous helium (carrying gas) and argon (makeup gas) combination flow through a T-connector by using an external standard of NIST 612 to correct for linear drift.

Cold Cathodoluminescence Analysis of Scheelite
The inner textural features such as oscillatory zoning and core, mantle and rim texture are difficult to study under a microscope. Therefore, OM-CL analysis has been widely used in research in recent years to examine the inner texture of minerals. However, fine interior texture, such as various phases of mineral formation is difficult to identify. However, scheelite can be analyzed because of its fluorescence characteristics by using cold cathodoluminescence method [5,16]. By using this technique, the detailed texture, filling, and cutting relationships of several mineral phases are noticed, which can be explained furtherly in discussion part.

EMPA Mapping of Garnet
Garnets were chosen because of their larger grains and potential association with endoskarn scheelite to explore the major element variation of early-stage skarn minerals. Comparing the proportion of grossular and andradite in garnet, we can apparently demonstrate the redox environment of primary skarn minerals formation [17][18][19]. The scheelite grains coexisted with garnet and seemed to have a similar growing method from core to rim, so the brighter part of scheelite ( Figure 5) can be seen as the faster growing part of scheelite, which was influenced by a high disturbed fluid environment. Therefore, REE analyses are performed to better understand the mineralization of scheelite in Tongshankou deposit. interior texture, such as various phases of mineral formation is difficult to identify. However, scheelite can be analyzed because of its fluorescence characteristics by using cold cathodoluminescence method [5,16]. By using this technique, the detailed texture, filling, and cutting relationships of several mineral phases are noticed, which can be explained furtherly in discussion part.

EMPA Mapping of Garnet
Garnets were chosen because of their larger grains and potential association with endoskarn scheelite to explore the major element variation of early-stage skarn minerals. Comparing the proportion of grossular and andradite in garnet, we can apparently demonstrate the redox environment of primary skarn minerals formation [17][18][19]. The scheelite grains coexisted with garnet and seemed to have a similar growing method from core to rim, so the brighter part of scheelite ( Figure 5) can be seen as the faster growing part of scheelite, which was influenced by a high disturbed fluid environment. Therefore, REE analyses are performed to better understand the mineralization of scheelite in Tongshankou deposit. According to the mapping of garnet, it is so clear that garnet has an evolution tendency from grossular to andradite. When the results of garnet EMPA mapping, major elements, garnets crystallography, and dynamic recrystallization are combined, then distal garnet rims have higher Al and Fe contents ( Figure 6), which can be seen as regeneration from grossular to andradite. Therefore, these crystallographic processes reveal the hydrothermal evolution that has effects on scheelites. According to the mapping of garnet, it is so clear that garnet has an evolution tendency from grossular to andradite. When the results of garnet EMPA mapping, major elements, garnets crystallography, and dynamic recrystallization are combined, then distal garnet rims have higher Al and Fe contents ( Figure 6), which can be seen as regeneration from grossular to andradite. Therefore, these crystallographic processes reveal the hydrothermal evolution that has effects on scheelites. SchB has an evolution tendency from inner grossular zones (which has relative lower Fe content) to andradite zones within distal rims (which has obvious higher Fe content). This can result from a highly disturbed growing environment, and andradite crystalized faster than grossular in this environment [20]. SchB has an evolution tendency from inner grossular zones (which has relative lower Fe content) to andradite zones within distal rims (which has obvious higher Fe content). This can result from a highly disturbed growing environment, and andradite crystalized faster than grossular in this environment [20].

Major and Trace Elements in Scheelites
The 62 EMPA results are summarized in Table 1. SchA has the highest WO 3 contents (average of 78.5 wt.%) and SchB has around 75.8 wt.% in average, while the SchC has the lowest WO 3 contents (74.6 wt.% in average), respectively. Molybdenum is more susceptible to redox conditions than W element due to its two distinct electric valences. SchB has the highest average MoO 3 concentration (4.85 wt.%), as well as the highest single grain content (approximately 25 wt.%), while MoO 3 contents of SchC are often greater than those of SchA ( Figure 7 and Table 1). Molybdenum is more susceptible to redox conditions than W element due to its two distinct electric valences. SchB has the highest average MoO3 concentration (4.85 wt.%), as well as the highest single grain content (approximately 25 wt.%), while MoO3 contents of SchC are often greater than those of SchA ( Figure 7 and Table 1).
SchA, mainly within exoskarn, has the highest WO3 contents and lowest MoO3 contents, indicating a reducing crystallization state (Figure 7a). The MoO3 concentration and occurrence of 15 scheelite grains were categorized based on their substantial differences ( Figure 7b). Furthermore, under a microscope and in OM-CL images, these distinct grains exhibit obvious mineralogical features that may be identified. The MoO3 concentrations in SchC (inside the garnet skarn) are higher compared with SchA and SchB formed within the sulfides stage and were mostly observed in the molybdenum ore sample, so it has the highest MoO3 content. Strontium, Pb, Nb, Na and Ba are the trace elements observed in scheelite with considerable abundances. There are many spots where Na concentrations are below the detection limit, while within the same skarn Nb and Sr have relatively similar concentrations. In Tongshankou deposit, Sr contents of scheelites range variously and have an average of 340.4 ppm. Niobium abundances are about 13 ppm (average) within exoskarn, whereas they are about 28 ppm (average) in endoskarn ( Table 2). ΣREE substitution mechanism analyses can be applied for the correlation of Nb and Sr elements abundance.  SchA, mainly within exoskarn, has the highest WO 3 contents and lowest MoO 3 contents, indicating a reducing crystallization state (Figure 7a). The MoO 3 concentration and occurrence of 15 scheelite grains were categorized based on their substantial differences (Figure 7b). Furthermore, under a microscope and in OM-CL images, these distinct grains exhibit obvious mineralogical features that may be identified. The MoO 3 concentrations in SchC (inside the garnet skarn) are higher compared with SchA and SchB formed within the sulfides stage and were mostly observed in the molybdenum ore sample, so it has the highest MoO 3 content.
Strontium, Pb, Nb, Na and Ba are the trace elements observed in scheelite with considerable abundances. There are many spots where Na concentrations are below the detection limit, while within the same skarn Nb and Sr have relatively similar concentrations. In Tongshankou deposit, Sr contents of scheelites range variously and have an average of 340.4 ppm. Niobium abundances are about 13 ppm (average) within exoskarn, whereas they are about 28 ppm (average) in endoskarn (Table 2). ΣREE substitution mechanism analyses can be applied for the correlation of Nb and Sr elements abundance.
Several points have some elements below the detection limit, which made Lu seem to have abnormal positive anomalies. Consequently, the "bdl" points have been wiped out to weaken the test errors.
Compared with REE characteristics of scheelites from other skarn deposits [36][37][38][39][40][41], analyzed results of samples in the Tongshankou deposit are mostly Eu positive. SchA and B have higher ΣREE and stronger Eu positive anomalies and SchC usually has similar LREE-enriched and HREE-depleted REE patterns.
Scheelites in the Tongshankou skarn deposit have the evolution process from "flat" and relatively high ΣREE contents to enrichment in LREE and depletion in HREE patterns, which is consistent with bulk rock background value.
whereas La, Ce and Eu anomalies represent magmatic or hydrothermal fluid redox conditions [28][29][30][31][32][33]. The Tongshankou deposit also has a porphyry system that is related to oxidized magma [2,3,6,12,34]. The diorite and skarn are both enriched in LREE and depleted in HREE with a relatively strong slope of LREE type and a weak negative Eu anomaly [34,35]. To better examine the scheelites REE geochemistry, the Eu anomaly was focused to understand the evolutionary process ( Figure 8).  (2); 3Ca 2+ ⇌ 2REE 3 + + □Ca (□Ca = Ca site vacancy) (3) [36,42,43] involved in schee tallization. Nb vs REE-Eu + Y and Eu vs Eu/Eu* ratios have been plotted to dem the substitution mechanism of scheelite in the Tongshankou deposit (Figure 9c) Nb vs ΣREE-Eu + Y (Figure 9a), ΣREE-Eu + Y are used to eliminate Eu anomaly i EuN vs Eu*N/Eu* were used to preclude Ca content influence in the mechanism (3 9b,c) [44], because Eu 2+ ions tend to enter the scheelite lattice to fill in the Ca vaca In the substitution mechanism (2), SchC exhibits a clear linear trend, indicating that it is the main mechanism in endoskarn mineralization (Figure 9a). The majority of SchA and a portion of SchB exhibit a correlation between Eu N and Eu/Eu*, implying that Eu anomalies in SchA and SchB should be triggered by 3Ca 2+ 2REE 3+ + Ca ( Ca = Ca site vacancy) (Figure 9b). The fact that Eu anomaly has increased together with ΣREE, indicates that Eu has a more distinct entry mechanism than other REEs (Eu 2+ ions enter the Ca vacancy). During crystallization, some SchC grains, particularly the "double linear grain" may have numerous substitution mechanisms (Figure 9c). SchA has a clear linear pattern, and the data was plotted to find if it was corresponding to mechanism (3), and it has a strong correlation with a coefficient of 0.96.

W-Mo Substitution Mechanism
The substitution of W-Mo contents in scheelite has been studied in detail, which is one of the most important influencing factors for the formation of tungsten ore [45,46] (Figure 7a). Scheelites have been shown to exhibit various mineralization processes even within the same sample [5,47]. Although W and Mo show a clear negative correlation (W 6+ and Mo 6+ replacement) (Figure 7a), redox influences in the concentration of W and Mo are not clearly demonstrated. In general, oxidized scheelite with more Mo content coexists with skarn minerals that formed earlier than minerals in oxidation stage, and reducing scheelite coexists with sulfides, quartz and calcite. Scheelite within endoskarn contains higher Mo contents and lower W contents, which may indicate an early oxidized fluid. When the distal skarn is concerned, reducing scheelite with a higher W content and a lower Mo content is crystallized with sulfides and ore-bearing quartz veins (Figure 7b).
The substation mechanisms (3) of Ca 2+ + W 6+ REE 3+ + Nb 5+ (2) and 3Ca 2 + 2REE 3+ + Ca ( Ca = Ca site vacancy) were compared (Figure 9c), and element plots were used to determine the major mechanism during ore formation. The results show that mechanism (2) controls the progression of SchC, which is mostly based on linear and fine correlation, while mechanism (3) controls the formation of SchA and SchB. Scheelites may have undergone numerous substitutions during the ore formation as a result of the changing redox conditions of environment fluid. Actually, substitution mechanisms (2) and (3) have shown to be active (or present) in a single grain of SchC scheelite type. Further investigation of its texture under SEM-CL images and OM-CL images revealed that it originally had distinct later filling and cutting texture with different optical properties (another generation of scheelite) and its mottled oscillatory zoning should be evidence of this process.

Multistage Ore-Forming Process-"Double Linear Grain"
"Double linear grain" refers to the progression of two of three substitution mechanisms in one grain, which is most common in endoskarn. The grain is assumed to have undergone multi-stage mineralization (SchI-early stage and SchII-late stage) by analyzing its internal texture in cold cathodoluminescence ( Figure 10). It is clear that two substitution mechanisms played an important role in this "double linear grain" (Figure 9c). These kind of grains can be evidence of at least two substitutions occurring during the formation of a single scheelite. In fact, the other two of the three selected SchC grains have a lower slope linear trend such as SchA (Figure 9c), indicating disseminated, cutting texture and oscillatory zoning. These textures may represent different stages of scheelite These kind of grains can be evidence of at least two substitutions occurring during the formation of a single scheelite. In fact, the other two of the three selected SchC grains have a lower slope linear trend such as SchA (Figure 9c), indicating disseminated, cutting texture and oscillatory zoning. These textures may represent different stages of scheelite mineralization and their own crystallization process.
Eu N vs. Eu/Eu* plots are used in combination with REE partition mapping to determine if the entry of Eu 2+ ions was the primary factor for Eu/Eu* formation (Figures 8a-d and 9a-c). Higher linear relationship indicates that Eu ions have more than "Eu 3+ " form when entering scheelite lattice compared to other REE (as shown in mechanism (3)) (Figure 9b). At a late stage, scheelites have similar substitution and forming mechanism, which can be seen as a transition from endoskarn to exoskarn.
The slope of linear trend in plots would be lower (at the lowest it would be a horizonal line) if the Eu N was mostly made up of the Eu 3+ ions. However, it was significantly greater in Figure 9b, indicating that Eu 2+ ions penetrate the lattice mostly in exoskarn SchB and SchA, rather than endoskarn SchC (relatively oxidating).
In fact, this could have originated from the geochemical characteristic of the Tongshankou porphyry source region: High Sr/Yb ratio at the average of 73.5, enrichment in LREE, depletion in HREE, high silica-alkali ratios [34,48,49], which illustrates the typical adakite signal in Tongshankou intrusions. The Sr/Yb ratio can be considered as an overview of the P-T conditions source region [50], and the high Sr/Yb and weak negative Eu anomaly revealed very little fractional crystallization of plagioclase in the source region of Tongshankou magma [36,51,52]. Therefore, the relatively weak Eu anomaly in SchC (lower values compared with SchA and SchB) could be linked to the geochemical characteristics of intrusions. Eu ions were mostly induced to enter plagioclase during early skarn mineralization. When early formed plagioclase dissolved and was released into the fluid, then Eu re-entered the hydrothermal fluid with other metallogenic elements and was transported to the exoskarn. In the exoskarn, relative reducing conditions caused Eu 2+ to be more abundant than Eu 3+ , and mechanism (3) is dominated by crystallization. Therefore, the positive Eu anomalies are widely distributed in REE patterns of Sch A and B (Figure 8a-d).
However, the uncoordinated spots continued to appear in the REE distribution pattern, and these spots were more prominent in those grains with distinct core, middle and rim partition and had strong positive Eu anomalies (Figure 8c,d).

Internal Texture of "Uncoordinate" Grains
These uncoordinated spots points are mostly observed in OM-CL images as zones with lighter optical properties as compared to others, but there is no obvious cutting and filling texture under a microscope ( Figure 11). It revealed that a random distribution of crystal-growing environment (hydrothermal fluid) eventually led to the formation of composite scheelite grains. Therefore, these uncoordinated spots can be seen as a special footprint indicating multi-stage hydrothermal mineralization of skarn deposit. In relation to previous discussion of substitution mechanisms, the Eu anomaly should be an indication of late-stage reducing hydrothermal fluid activity [44,46].
The tungsten concentration of reduced scheelite bearing skarns is greater than that of oxidized scheelite bearing skarns. [16,46,48]. The reducing hydrothermal fluid formed different types of scheelite grains in the early skarn stage mineralization within the endoskarn. The reducing W skarn deposit has significantly more tungsten concentration because of the fluid transportation process. Tungsten should be carried better in this process than under oxidized conditions because tungsten can be carried in ionic complex form during transportation, and competes with Mo, which has the same 6+ electronic valence. This may elucidate why scheelite is commonly associated with sulfides and reduced minerals, particularly molybdenite. The mafic enclaves in the Tongshankou district have the highest concentration of W [8]. Such mafic enclaves hosted in the ore-forming porphyries were most likely generated by a magma mixing process [6]. Therefore, W has a different magma source from the oxidized porphyry magma.

Internal Texture of "Uncoordinate" Grains
These uncoordinated spots points are mostly observed in OM-CL images as zones with lighter optical properties as compared to others, but there is no obvious cutting and filling texture under a microscope ( Figure 11). It revealed that a random distribution of crystal-growing environment (hydrothermal fluid) eventually led to the formation of composite scheelite grains. Therefore, these uncoordinated spots can be seen as a special footprint indicating multi-stage hydrothermal mineralization of skarn deposit. In relation to previous discussion of substitution mechanisms, the Eu anomaly should be an indication of late-stage reducing hydrothermal fluid activity [44,46]. Figure 11. The characteristics of red curves, which are "uncoordinated" among SchA (a), SchB (b), SchC (c): diagrams on the left are REE distribution patterns of scheelite, the middle photos are cold cathodoluminescence images and the photos on the right are corresponding SEM-CL images of the grains. The normalization values were taken from Taylor and McLennan (1985) [33].
The tungsten concentration of reduced scheelite bearing skarns is greater than that of oxidized scheelite bearing skarns. [16,46,48]. The reducing hydrothermal fluid formed different types of scheelite grains in the early skarn stage mineralization within the endoskarn. The reducing W skarn deposit has significantly more tungsten concentration because of the fluid transportation process. Tungsten should be carried better in this process than under oxidized conditions because tungsten can be carried in ionic complex form during transportation, and competes with Mo, which has the same 6+ electronic valence. This may elucidate why scheelite is commonly associated with sulfides and reduced minerals, particularly molybdenite. The mafic enclaves in the Tongshankou district have the highest concentration of W [8]. Such mafic enclaves hosted in the ore-forming porphyries

Conclusions
This research presents the first detailed investigation of scheelites in skarns at the Tongshankou deposit where three types of scheelites (SchA, B and C) have been identified. The following geochemical and CL features of different scheelites can be concluded:

1.
The geochemistry and OM-CL show that the W mineralization has been highly disrupted during the formation of scheelite. The skarn mineralization system has an apparent evolution in the redox environment from endoskarn (SchC) to the distal exoskarn (SchA). W-Mo concentrations have fluctuation due to the crystallization of multi-stage minerals. SchC is considered to have similar ΣREE features as intrusions with the negative Eu anomaly, whereas SchA and SchB have higher ΣREE characteristics and a comparatively high positive Eu anomaly. 2.
"Double linear grain", postulated for the first time in this work, could be one of the evidences of two substitution mechanisms in a single scheelite grain, which is the typical characteristic of SchC. In addition, the cutting and filling texture in scheelites are also prominent in OM-CL images. Early endoskarn scheelites were dominated by mechanism Ca 2 + + W 6 + REE 3 + +Nb 5+ (2), in scheelite crystallization and there was no appropriate reducing environment to sustain mechanism (3) 3Ca 2 + 2REE 3 + + Ca ( Ca = Ca site vacancy). The reducing fluid metasomatized the carbonate strata, traversed and regenerated some earlier scheelites following mechanism (3), and caused a positive Eu anomaly in the late skarn stage.

3.
The ore formation at the Tongshankou area is a multistage process, including a porphyry system at~143 Ma and a skarn system at~140 Ma. The porphyry system has an obvious priority relationship in time scale compared with the skarn system and it was probably the source of both the metallogenic element and the hydrothermal fluid. According to this study, the oxidation mineralization of Cu and Mo may have different magma sources for W. Scheelites formed with mechanism (2) during the early skarn stage mineralization within endoskarn, when there was no appropriate reducing environment for fluid transportation. Both the mechanism Ca 2 + + W 6 + REE 3 + +Nb 5+ (2) and mechanism 3Ca 2 + 2REE 3 + + Ca ( Ca = Ca site vacancy) (3), indicate fluctuation in ΣREE in Tongshankou scheelite, suggesting multi-stage mineralization: At the first stage, mechanism (2) hosted within endoskarn proceeded through mechanism (3), and the ΣREE increased throughout the single systematical skarn mineralization. Those grains with uncoordinated REE characteristics might be indicative of a reducing metallogenic fluid entering the skarn system during the late stage of evolution.

4.
The association of skarn W mineralization with an oxidized porphyry Cu (Mo) deposit in Tongshankou could be explained by magma sources coming from various locations and being transported with distinct redox fluids. Endoskarn and exoskarn scheelites exhibit distinct growing environments and W-Mo concentrations, indicating that W may have a favorable transportation mechanism from an oxidation to a reducing environment. Mo-rich scheelites are formed from an early metallogenic fluid and when it evolved at a later stage, magma from different source regions upwelled, generating a more reducing environment, though contamination by some earlierformed intrusions. Data Availability Statement: All data generated or analyzed during this study are included in this published article.