Next Article in Journal
Simulation Technology Development for Dynamic Analysis of Mechanical System in Deep-Seabed Integrated Mining System Using Multibody Dynamics
Next Article in Special Issue
Genesis of the Wuzhutang Granite and Associated W–Sn–Be Mineralization in the Xuebaoding Mining Area, Sichuan Province, China
Previous Article in Journal
A Thermodynamic and Experimental Assessment of the Recovery of Copper, Iron, Zinc, and Lead from Copper Slag
Previous Article in Special Issue
A Combined EMPA and LA-ICP-MS Study of Muscovite from Pegmatites in the Chinese Altai, NW China: Implications for Tracing Rare-Element Mineralization Type and Ore-Forming Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Chemical Characteristics and Metallogenic Mechanism of Beryl from Cuonadong Sn-W-Be Rare Polymetallic Deposit in Southern Tibet, China

Chengdu Center, China Geological Survey, Chengdu 610081, China
*
Author to whom correspondence should be addressed.
Minerals 2022, 12(5), 497; https://doi.org/10.3390/min12050497
Submission received: 21 February 2022 / Revised: 11 April 2022 / Accepted: 15 April 2022 / Published: 19 April 2022
(This article belongs to the Special Issue Rare Metal Ore Formations and Rare Metal Metallogeny)

Abstract

:
The Cuonadong deposit is the first large scale Sn-W-Be rare polymetallic deposit located in southern Tibet, China, where beryl is the main beryllium-bearing mineral. In this paper, the beryl crystals in the pegmatitic and hydrothermal vein orebody from the Xianglin area of the Cuonadong deposit are the research objects, marked as Beryl-I and Beryl-II, and they are investigated by EPMA, LA-ICP-MS and in situ micro-X-ray diffraction (XRD). Data by EPMA and LA-ICP-MS reveal that beryls from this area are alkaline beryls, among which Beryl-I is composed of Li-Cs beryl, and Beryl-II is composed of Na beryl and Na-Li beryl, indicating that beryls have undergone noticeable alkali metasomatism during formation. The Cs/Na ratio in Beryl-I ranges from 0.10 to 0.44, and the Mg/Fe ratio is almost 0, showing that Beryl-I is formed under high-differentiation evolution conditions and is rarely affected by hydrothermal transformation, whereas the Mg/Fe ratio in Beryl-II ranges from 2.73 to 17.31, and the Cs/Na ratio is nearly 0, indicating that Beryl-II has been obviously affected by late hydrothermal metasomatism. In situ XRD analysis shows that both Beryl-I and Beryl-II are t-beryl, and the c/a ratio of Beryl-I (1.0010–1.0012) is slightly higher than that of Beryl-II (1.0005–1.0008), which may also reflect the transition from magmatism to hydrothermal metasomatism in the late stage of pegmatitic magmatism. Based on comprehensive analysis, we believe that the precipitation of Beryl-I is mainly caused by the emplacement of highly fractionated magma containing Be to the top of the rock mass or surrounding rock, the melt-fluid undercooling, and the crystallization of volatile-bearing minerals (such as tourmaline and fluorite). Moreover, the Be-bearing ore-forming fluid has further migrated upward along the near north–south faults formed in the middle Miocene (16–15 Ma), during which Beryl-II precipitates owing to the hydrothermal water mixing, the ore-forming fluid cooling, and large amounts of crystallization of volatile-bearing minerals (mainly fluorite). Therefore, it can be concluded that beryl mineralization largely reflects the process of magmatic–hydrothermal mineralization. Because of a large number of mineralized areas with the similar metallogenic backgrounds to the Cuonadong deposit in the Himalayan region, it has great potential to be a new globally significant rare metal metallogenic belt.

1. Introduction

Beryllium, the lightest alkaline earth metal element, is highly incompatible; thus, it is commonly present in granitic residual melts and enriched in evolved granites and granitic pegmatites [1,2]. The formation of beryllium minerals runs through the whole magmatic–hydrothermal evolution process, including beryl, chrysoberyl, phenacite and bertrandite [3,4,5,6]. Beryl is the most common beryllium mineral, which is mostly present in peraluminous metallogenic systems and can be formed at multiple stages of the magmatic–hydrothermal transition. According to the evolution stage of the magma–hydrothermal system, it is further divided into three types: magmatic, granitic pegmatite and hydrothermal type [1,7]. Beryl is the most important Be-hosting mineral in granitic pegmatite-type beryllium deposits [6]. During the magmatic–hydrothermal evolution of granitic pegmatite, the crystallization of beryl is predominantly controlled by temperature and pressure [6]. In addition, the stability of beryl in the melt is also affected by the activity of Si, Al and other components in the melt, such as Li, B, P, F, Rb, Cs, etc. [4,5]. Small changes in activity of Be2+ can alter the crystallization behavior of beryl with changes in the composition and temperature of the pegmatitic melt [2,5].
In addition to being an important carrier of Be, beryl can also carry a certain amount of Li, Cs, Na, Fe, Mg, and its chemical composition and structure can be used to indicate the degree of differentiation of pegmatites and trace magmatic–hydrothermal evolution [8,9]. Studies have shown that the higher the Cs content in beryl is, the lower the Na/Cs value, Na/Li value, Mg/Fe value and Mg content are, indicating the higher degree of differentiation and evolution of granitic pegmatite [8,9,10,11].
Beryl, ideally Be3Al2Si6O18, has a ring structure consisting of layers of rings with six Si-O tetrahedra(T1), which are linked vertically and laterally by Be-O tetrahedra (T2) and Al-O octahedra(O), and the centers of the Si rings contain channels [12,13]. Various cations may be incorporated in the beryl structure, resulting in complex substitutions [6]. Therefore, the crystal structure and chemical composition of beryl can reflect the physical and chemical conditions of its formation [1,2,9,10,11,12,13].
The Cuonadong W-Sn-Be rare polymetallic deposit, located in southern Tibet, China, was first discovered in the Himalaya region in recent years, with superlarge-scale prospecting potential [14,15]. The mineralization types are complex, including the pegmatite Be orebody in Miocene pegmatitic leucogranite, the skarn W-Sn-Be orebody in the strongly deformed mantle of the dome, and the hydrothermal vein-type Sn-W-Be orebody faulted by north–south-trending rifts (NSTRs). Preliminary exploration of the W-Sn-Be orebody in the skarn around the dome mantle of the resources are WO3 > 50,000 tons (average grade 0.21%), Sn > 80,000 tons (average grade 0.36%), and BeO > 170,000 tons (average grade 0.09%) [14,15], and they have attracted the attention of many scholars, mainly focusing on geological characteristics [16,17], petrology [18,19,20,21], chronology [22,23,24,25,26,27,28,29], geochemistry [30,31] and basic mineralogy [32,33,34]. Beryl is the most important beryllium-bearing mineral in the Cuonadong deposit [17,28]. Although some scholars have conducted some preliminary research on the composition of beryl in this deposit, they have conducted a rough discussion on the element substitution mechanism of beryl and its effect on pegmatites [33,34]. However, the research on the formation mechanism of beryl is incomplete.
In this study, beryl crystals from pegmatitic and hydrothermal mineralization in the Cuonadong deposit as the research objects are investigated by energy spectrum scanning, electron probe, LA-ICP -MS, in situ microzone XRD, with a view to ascertain the physical and chemical characteristics of beryl from the two types of mineralization, and to discuss the Cuonadong deposit in the magmatic–hydrothermal evolution process of the rare metal Be metallogenic mechanism, which can provide evidence of the Himalayan metallogenic belt having further prospects in rare metals.

2. Geologic Setting

2.1. Himalayan Metallogenic Belt

The Tethyan Himalaya Sequence (THS) is separated by the Indus Yarlung Tsangpo suture zone to the north and the South Tibet Detachment System (STDS) to the south [35,36,37], sandwiched between the Great Himalayan Crystalline Complex (GHC) and the Lhasa terrane [38] (Figure 1). As the largest highly differentiated granite belt developed in the Himalayas, the study of Cenozoic tectono-magmatic evolution and metallogenic relationship between highly differentiated leucogranites and their rare metals in this region has been a hot topic at home and abroad in recent years.
Previous studies revealed that the Himalayan Metallogenic Belt is mainly composed of Pb-Zn-Ag-Au-Sb deposits [39,40,41,42,43]. Meanwhile, researchers have been paying attention to the mineralization of rare metals such as beryllium, tungsten, and tin in the Himalayas. Many rare metal minerals have been found in leucogranites and pegmatites, including lithium-rich minerals such as lithium tourmaline and spodumene [44], and other rare metal minerals such as scheelite, cassiterite, beryl, niobium rutile, columbite–tantalite, iron–lepidolite, spodumene, rubidium-containing mica, etc. [44,45,46,47,48]. Beryl, as a beryllium-rich mineral, was found in the Pusila, Wengbo, Cuona, Yadong, Kangma, Gaowu, Dingjie, Xiaru, Peguco and other leucogranite bodies [47]. The latest exploration results show that Be-Rb-Nb-Ta-W-Sn and other rare metal mineralizations are developed in and around the Himalayan domes (Figure 1). Preliminary studies show that this region has great metallogenic potential and is expected to become an important rare metal metallogenic belt in China [47,48,49,50].

2.2. Geology of Cuonadong Mining Area

The spatial distribution of the Cuonadong Sn-W-Be polymetallic deposit is strictly controlled by the structure of the Cuonadong gneiss dome, which is typically representative of rare metal mineralization in the Himalayan gneiss dome belt [14,15,16,17].

2.2.1. Cuonadong Gneiss Dome

The Cuonadong gneiss dome is shaped as a rhombus, which is 20 km in length, 20 km in width and approximately 400 km2 in area. The dome can be divided into three parts: core, mantle and rim [15,19,22] (Figure 2). First, the core is located under the lower detachment fault, consisting of early Paleozoic orthogneiss and Miocene leucogranite and pegmatite. Furthermore, the mantle and core are separated by the lower ductile detachment fault. The mantle shows typical characteristics of the ductile shear zone, formed from a set of strongly deformed quartz–mica schists with mylonitized carbonate rocks metamorphosed at Paleogene. Strong contact metasomatism occurred between the marble and leucogranite or granitic pegmatite, thus forming a discontinuous skarn zone around the dome where W-Sn-Be polymetallic ore bodies are found. In addition, the rim and mantle are separated by the upper brittle detachment fault, and the rim is dominated by a set of low-grade metamorphic rocks composed of phyllite, sandy slate and carbonaceous slate.
The Cuonadong dome has experienced three strong magmatism stages. The first stage occurred under the background of the dome activity of the south Tibet detachment system (STDS) in southern Tibet in the Oligocene (34~26 Ma) and was mainly distributed in the mantle and rim of the dome [50]. The lithology is mainly composed of strong deformed two-mica granites and granitic pegmatite. The second stage occurred in the lower detachment fault in the Early Miocene (24~18 Ma) [8,10]. The lithology consists of two-mica granites with weak deformation. Particularly in the vicinity of the lower detachment fault, the second-phase two-mica granite is directed with schist and marble and shows strong mineral orientations and lenticular shapes [18,20,30]. ③ The third stage occurred in the Middle and Late Miocene (16~14 Ma) with undeformed massive muscovite leucogranite and highly differentiated garnet–tourmaline-bearing muscovite leucogranite and beryl-bearing pegmatite [27,31]. The main body of the dikes is mainly distributed in a north–south direction, intrudes into various geological bodies and structural layers, and cuts the detachment fault in the form of veinlets, which is probably related to the north–south-trending rifts (NSTRs) (Figure 2).

2.2.2. Mineralization Characteristics

Four regions with important metallogenic potential are distributed around the dome, including Xianglin, Yunazha, Rina and Qugaer (Figure 2). Among them, the Xianglin area in the northwestern part of Cuonadong has the most abundant mineralization types and the largest inferred resources; thus, it is the key research area in this study and the main area of recent mineral exploration (Figure 3).
The main mineralization types at Xianglin include pegmatite Be, skarn W-Sn-Be and hydrothermal vein Sn-W-Be ores. Preliminary exploration of Xianglin shows that the resources of Sn, WO3 and BeO are over 80,000 t at 0.36%, 50,000 t at 0.21% and 34,000 t at 0.08% [27,28].
The pegmatite Be orebody is spatially closed to muscovite granite (16 ± 1 Ma), and a large amount of beryl-bearing pegmatite is developed on the top of the muscovite granite. The gradual and transitional contact between the two implies an evolutionary relationship. Beryl is the main beryllium-bearing mineral in the disseminated form, with crystals up to 10 cm and colors ranging from blue–green to pale green, which co-exists with quartz, plagioclase, K-feldspar, muscovite, tourmaline and fluorite (Figure 4a). The content of beryl in pegmatite dikes is not evenly distributed, and the scale of individual veins is tens of centimeters to a few meters, which has not yet been evaluated [14,17].
Skarn ore bodies are mainly W-Sn mineralization, followed by Be. Cassiterite is the dominant tin-bearing mineral, and scheelite is the main tungsten-bearing mineral (Figure 4b). Beryllium mainly occurs in silicate minerals, such as vesuvianite, albite and garnet, through element substitutions and a small amount of Be-containing minerals such as phenacite, bertrandite [32]. However, as far as the current mineral processing technology is concerned, beryllium in these skarn minerals is difficult to fully utilize [51].
Hydrothermal vein ore bodies are produced in high-angle and nearly north–south faults and detached faults, mainly exhibiting greisenization, albitization, silicification, fluoridization and pyritization. The mineralization includes cassiterite–quartz veins (Figure 4c) and cassiterite sulfide veins (Figure 4d), which is superimposed on skarns, or occurs on the top or in the edge of granitic pegmatites, which developed later than weakly oriented mica granites, indicating that hydrothermal veins may be formed in garnet-bearing mica granites. The latest exploration results show that 11 new hydrothermal vein Sn-W-Be ore bodies have been discovered in the Xianglin area, among which one main ore body (No. Z6 orebody) has the predicted potential resources of Sn 73,800 tons at 1.75%, WO3 19,900 tons at 0.36% and BeO 3100 tons at 0.14% [51]. Beryllium occurs mainly in beryl, and tin and tungsten occur in cassiterite and scheelite, respectively. The hydrothermal vein Sn-W-Be ore bodies are the most significant mineralization in the Cuonadong mining area [17,51].

3. Analytical Methods

In this paper, two beryl crystals from pegmatitic and four from hydrothermal mineralization of the Xianglin area are taken as the key research objects (Figure 5), and the location of the samples is shown in Figure 3.
The samples under investigation were characterized through a combination of optical examination, combined with scanning electron microscopy, electron microprobe analyses and laser ablation inductively coupled plasma mass spectro metry.
Major- and minor-element compositions of beryls were conducted on a JEOL JXA-8230 electron microprobe at the Key Laboratory of Metallogeny and Mineral Resource Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences. The analytical conditions were: an accelerating voltage of 15 Kv, a beam current of 20 nA, and a beam size of 5μm. The following standards were used: albite (Na), sanidine (Si,Al,K), pyrope (Mg), titanite (Ca,Ti), almandine (Fe), rhodonite (Mn), Rb-leucite (Rb) and pollucite (Cs). The analytical data were corrected using the ZAF methods.
Trace-element compositions of beryls were performed using an X-Series ICP-MS (Thermo Fisher Scientific, Braunschweig, Germany) coupled with a J-200 343 nm Yb: Fiber femto-second laser ablation system (Applied Spectra, West Sacramento, CA, USA) housed at the National Research Center for Geoanalysis, Chinese Academy of Geological Sciences. Helium gas carrying the ablated sample aerosol from the chamber was mixed with argon make-up gas and nitrogen as an additional diatomic gas to enhance sensitivity. A baffled-type smoothing device was used in front of the ICP-MS to reduce fluctuation effects induced by laser-ablation pulses and to improve the quality of the data [52]. The beryl samples were ablated for 50 s at a repetition rate of 8 Hz at 6 J/cm2, and ablation pits were ~50 μm in diameter. The calibration was performed externally using two NIST SRM 610 and one NIST SRM 612 for every 10 samples with 29Si, measured by electron microprobe, as the internal standards to correct for instrument drift. Data reduction was carried out with the commercial software ICPMS Data Cal 10.8 [53,54]. The detection limits of LA-ICP-MS ranged from 0.05 to 0.1 ppm for REE. Repeated analyses of the standards SRM 610 and SRM 612 indicate that both precision and accuracy were better than 10% for most analyzed elements.
The mapping of beryl was determined using a NWR 193 nm ArF Excimer laser-ablation system coupled to an iCAP RQ (ICPMS) at the Guangzhou Tuoyan Analytical Technology Co., Ltd., Guangzhou, China. The ICPMS was tuned using NIST 610 standard glass to yield low oxide production rates. Then, 0.7 L/min He carrier gas was fed into the cup, and the aerosol was subsequently mixed with 0.89 l/min Ar make-up gas. The laser fluence was 5 J/cm2, with a repetition rate of 20 or 30 Hz and a 5, 8 or 10 μm spot size. The laser scan speed was 15, 20 or 35 um/s. The raw isotope data were reduced using the “TRACE ELEMENTS” data reduction scheme (DRS). The DRS runs within the freeware IOLITE package of Paton [55]. In IOLITE, user-defined time intervals are established for the baseline correction procedure to calculate session-wide baseline-corrected values for each isotope. Blocks of two standards (one NIST 610 and one GSE-2G) and one MASS-1 sulfide standard analyses were followed by 5 to 8 unknown samples. 29Si in beryl measured by electron microprobe was used as internal standard to correct instrument drift.
In situ micro-X-ray diffraction (XRD) experiment was conducted with a Rigaku D/max Rapis IIR micro-XRD system at the Guangzhou Institute of Geochemistry, Chinese Academy of Sciences. All measurements were carried out at the 40 kV and 250 mA (Cu Kα radiation) with a collection time of 22 min, 2θ range is 20~110 °C, and continuous scan. The scan speed was 4°/min. The X-ray beam was about 40 μm in diameter, focused on the selected spots on the thin sections. The data point interval was 0.045°. The software of MDI Jade 6.0 (International Centre for Diffraction Data, Denver, CO, America) was used to analyze the obtained XRD patterns for the derivation of unit-cell parameters by whole pattern fitting.

4. Results

4.1. Occurrence

Beryl from pegmatite (Beryl-I) is mostly symbiotic with rock-forming minerals such as quartz, potash feldspar, albite, muscovite and fluorite, with well-developed hexagonal column shape and particle sizes of 0.1–5 cm (Figure 5a,b). BSE images of these beryl particles show a homogeneous structure and locally developed altered edges, which is metasomatized by quartz, muscovite, K-feldspa, and albite along the edges or fissures (Figure 5b,c).
Beryl from the hydrothermal veins (Beryl-II) is mainly associated with cassiterite, scheelite, quartz, muscovite, phlogopite, fluorite, albite, sphalerite, and arsenopyrite, with the sizes ranging from 0.1 to 20 mm (Figure 5d,e). It is typically characterized by the development of complex zonal structures, which are locally replaced or dissolved by later K-feldspar, albite, and scorodite (Figure 5e,f).

4.2. Chemical Composition

As Li and Be are light elements (i.e., the X-ray yield through electron excitation is poor), and it is not possible to obtain precise Li2O and BeO analyses on an electron microprobe. Furthermore, Be is a major element in beryl; thus, the results measured by LA-ICP-MS are inaccurate [33]. Assuming that the crystallochemical formula of beryl is X3Y2[Z6O18], the anion method is adopted for calculation, and the theoretical number of oxygen atoms is 18 [56]. The theoretical atomic number of Be, herein, is calculated by taking the sum of the cations at the X site as 3. Meanwhile, assuming that Be at the X site is only replaced by Li, the theoretical atomic number of Be can be calculated by Be = 3-Li apfu [56]. The quantification of H2O was estimated as H2O wt.% = (0.845958 × Na2O in wt.%) + 0.08373 [57].
In this paper, two methods are used to obtain the Li content in beryl, one is through chemical calculation of EPMA (Li = Ca + Na + K + Cs, the number of atoms obtained by stoichiometric calculation) [58], the other is direct determination by LA-ICP-MS. By comparing the data results of these two methods, we found that Li content in beryl calculated by EPMA (0.04–0.87 wt.%) is significantly higher than that directly determined by LA-ICP-MS (1.23–2.48 wt.%). The reason is that there is complex metal replacement in beryl [59]. Li+ replacement of Be2+ results in a charge deficit, which requires monovalent alkali metal to enter the channel to balance the charge. This is the same as the replacement process of Al2+ by Fe2+ and Mg2+ [34].
Beryls from the Xianglin area contain different amounts of FeO (0.17–1.31 wt.%) and MgO (0~2.06 wt.%), leading to the high content of Li obtained by chemical calculation of EPMA. It implies that the content of Li in beryls is not completely consistent with the sum of other alkali metals (Figure 6). Thus, the content of Li calculated by using the EPMA component may lead to deviation. Moreover, the deviation of the Li content calculation value will further affect the content of BeO. Therefore, in this paper, the content of Li directly measured by LA-ICP-MS is used to calculate the content of BeO. By calculation, the BeO content of beryl in pegmatite is 11.69~12.05 wt.%, and that in the hydrothermal vein is 12.27~13.75 wt.%. The H2O content of the two beryls is 1.41~1.60 and 1.26~2.25 wt.%.
The major, minor and trace elements of beryls from the Xiangling area are listed in Table 1 and Table 2. The beryl is essentially composed of SiO2, Al2O3, BeO, and some typically minor elements, such as Li2O, Na2O, Cs2O, MgO and FeO (total Fe in the paper). Based on the chemical composition obtained for EPMA and LA-ICP-MS (Table 1 and Table 2), the formulas for the two types of beryl from the Xianglin area, using average values of the most important components, are:
(1) Beryl-I (beryl in pegmatite):
Ch(Na0.15,Cs0.04)0.19 T2(Be2.70, Li0.30)3.00 O(Al1.94,Fe0.03)1.97 T1Si6.03O18
(2) Beryl-II (beryl in the hydrothermal veins):
ChNa0.25 T2(Be2.89, Li0.11)3.00 O(Al1.76, Mg0.22, Fe0.05)2.03 T1Si6.00O18
where T2 and T1 are the tetrahedral sites for Be and Si, respectively, O is the octahedral site, and Ch is the “channel”.
A thorough analysis is presented in Section 5.1 and Section 5.2.

4.3. LA-ICP-MS Mapping

LA-ICP-MS mapping is used to analyze the complex zonal structure of beryl in the hydrothermal vein of Xianglin area. Judging from the mapping, the complex bands of beryl-II are mainly caused by the difference in the contents of Be, Al, Na, Cs, and Cr (Figure 7).

4.4. Unit Cell Parameters

The in situ micro-X-ray diffraction (XRD) analysis results (Table 3) show that the absolute values and average of the unit cell parameter a,b and the c/a ratio have a slight difference between the two types of beryls. Beryl-I has the slightly higher c value (9.2209–9.2243, avg.c = 9.2228) and c/a ratio (1.0010–1.0012, avg. 1.0011) than beryl-II (c = 9.2180–9.2206, avg.c = 9.2197; c/a =1.0004–1.0010, avg. c/a =1.0008).

5. Discussion

5.1. Crystal Chemical Features

According to the crystallochemical formula of beryl, various cations are incorporated in the beryl structure, resulting in complex substitutions. As divalent cations, such as Fe2+ and Mg2+, substitute trivalent Al3+ in the Y site (octahedral), leading to a charge deficit, the long-radius monovalent cations and divalent cation Ca2+ enter the channel to balance the charge [59,60,61,62,63,64]. Since the Ca2+ and K+ contents of the two types of beryls from the Xianglin area are negligible (Table 1), the occupying position of Ca2+ and K+ is not discussed in this paper. The atoms per formula unit of Be2+ in the X site of Beryl-I and Beryl-II is 2.70 and 2.89 apfu, respectively, showing obvious unsaturation. Li+ is the main substitute cation in the X site, and the atoms per formula unit are 0.30 and 0.11 apfu, respectively, indicating that the replacement degree of the X site ion of Beryl-I is higher than that of Beryl-II. The atoms per formula unit of Al3+ in the Y site of Beryl-I is 1.94 apfu, and the isomorphic substitution is 0.03 apfu, manifesting that the substitution degree of the Y site of Beryl-I is low, whereas the atoms per formula unit of Al3+ is 1.76 apfu, and the substitution of the Y site is 0.27 apfu, indicating that the substitution degree of the Y site of Beryl-II is relatively higher than that of Beryl-I.
Ideally, the sum of divalent cations in the Y site (R2+) and monovalent alkali cations in channel (A+) of beryl is close to 1:1 [65]. It can be illustrated with Figure 8 that all the points of Beryl-I lie at a high position above the 1:1 slope line, showing an obvious substitution in the channel. The points of Beryl-II are mostly located near the slope of R2+ = A+, and some points lie at a high position above the slope, indicating a certain amount of A+ entering the X and Y sites.
Linear fitting is performed on the main cation contents of the two types of beryl at the X site, Y site and channel (Figure 9). In Beryl-I, the negative correlation between Li and Na, Li and Cs, Na and Cs, and Fe and Al indicates that the mechanism of element substitution in Beryl-I includes (Na,Cs)Li□-1Be-1 channel-tetrahedral substitution, (Na,Cs)Fe2+-1Al-1 channel-octahedral substitution and NaCs-1 the mutual substitution of alkali metal ions in the “channel”. By comparison, in Beryl-II, Na and Li are weakly negatively correlated, Al is negatively correlated with Mg + Fe2+, Mg, Fe2+, and Mg is also negatively correlated with Fe2+, which can be inferred that the element substitution mechanism of Beryl-II is composed of NaLi□-1Be-1 channel-tetrahedral substitution and Na(Fe2+,Mg)□-1Al-1 channel-octahedral substitution.

5.2. Classification

As the ions are replaced, the unit cell parameters of the minerals also change. Based on the difference in the ratio of c/a, beryl is divided into three types: ① n-beryl (normal beryl, c/a = 0.997~0.998), ideal pure beryl without replacement; ② o-beryl (octahedral beryl, c/a = 0.991~0.996), of which Al3+ in the center is partially replaced; ③ t-beryl (tetrahedral beryl, c/a = 0.999~1.003), of which a large number of metal ions entering the structure channel to balance the electricity price [4]. The c/a ratios of Beryl-I and Beryl-II from the Xianglin area range from 1.0004 to 1.0012 (Table 3), which shows that the beryls from the Xianglin area are t-beryl. Significantly, the c/a ratios of Beryl-I and -II show a slight decreasing trend, which is probably caused by the accommodation of alkali elements and substitution of Be2+ and Li+ in the beryl structure along the c-axis [66]. It may reflect the hydrothermal process in the late stage of pegmatite evolution [9,10].
According to the contents of alkali metal elements, beryl can be classified into five types [3]: (1)alkali-free beryl (total R2O < 0.5 wt.%, where R = Na, Li, Cs, K); (2) alkali beryl (total R2O > 0.5 wt%); (3) Na beryl (total R2O > 0.5 wt%, Na2O = 0.5–2 wt%,Li2O = 0.1–0.5 wt%); (4) Na–Li beryl (Li2O = 0.5–1.5 wt% and Na2O = 1–2.5 wt%); and (5) Li–Cs beryl (Li2O = 0.1–1 wt%, Cs2O < 3 wt%, Na2O = 0.3–1 wt%). Therefore, the beryls at the Xianglin area are alkali beryls, of which beryl-I is Li–Cs beryl, with R2O = 1.97–3.35 wt%, Na2O = 0.68–0.90 wt%, Li2O = 0.74–0.87 wt%, Cs2O = 0.31–1.77 wt% (Figure 10a); while beryl-II consists of Na and Na-Li beryl, with R2O = 1.40–2.20 wt%, Na2O = 1.25–1.66 wt%, Li2O = 0.04–0.56 wt%, Cs2O = 0.00–0.39 wt% (Figure 10b).
Furthermore, the contents of alkali metal elements can indicate the magmatic–hydrothermal evolution trend of granitic pegmatite [67]; that is, the increasing content of Cs and Cs/Na ratio of beryl has been identified as the process of evolution and differentiation of pegmatite from low to high, with the possibility of decreasing content of Mg and Mg/Fe ratio. The high Mg/Fe ratio is generally considered to be affected by hydrothermal alkali-metasomatism [8,34].
Judging from the plot of Mg/Fe vs. Cs/Na (Figure 11), there are obvious evolutionary differences between the two beryls. The content of Cs2O of beryl-I in pegmatite ranges from 0.31 to 1.77 wt% (avg. = 1.01 wt%), and the ratio of Cs/Na(apfu) = 0.10–0.44 (avg. = 0.28), with almost no MgO detected. It is indicated that beryl-I was formed under highly differentiated evolution conditions [4]. By contrast, the content of MgO in beryl-II is higher, ranging from 0.46 to 2.25 wt.%, avg. = 1.48 wt.%, and the ratio of Mg/Fe(apfu) ranging from 2.73 to 17.31, avg. = 6.33, with the low content of Cs2O inferring that Beryl-II probably experienced hydrothermal alkali-metasomatism [10].

5.3. Metallogenic Mechanism

The solution and saturation of beryllium in granitic magma are mainly controlled by five factors [4,5,6,7]: temperature, pressure, melt composition, fractionation, and the volatile components such as B, P, and F.
The geological characteristics, metallogenic chronology, stable isotopes and fluid inclusions of the Cuonadong Be-Sn-W polymetallic deposit show that the ore-forming materials and fluids of the deposit are derived from Cuonadong leucogranite [14,15,16,17,18,19,20,21,22,23,24]. Cuonadong leucogranite, a typical representative of the Himalayan leucogranite belt, is a complex pluton formed by the multi-stage evolution of homologous magma, including the early strongly deformed two-mica granite, the middle weakly deformed two-mica leucogranite and the late undeformed muscovite leucogranite, of which the chemical composition gradually evolved to be of higher silicon, more alkalis and more volatile [19,20,21,22,23,24]. Beryllium, an incompatible element, would preferentially enter the melt phase during the differentiation and evolution of magma. With the continuous process of this interaction, volatile compounds and ore-forming fluids can often be released from the late magmatism [68]. Highly differentiated leucogranites generally have a high alumina saturation index (ASI), of which beryllium content is positively proportional to ASI value [5,69]. In addition, the beryllium content has a weak to moderate negative correlation with the potassium content, and the reduction of potassium is due to the homogeneous reaction of volatile anions such as F, P, and B with the Na-containing magmatic melt [5]. That is, the more potassium phases (e.g., K-feldspar, mica) that gradually dominate the crystallization mode, the more sodium remains in the melt [70]. Flux components (F, B, P) in pegmatite melts can form some stable fluoride, borate and phosphate migration with rare metals [4,71]; thereby, the transitional fluids with the coexistence of magmatic fluid and residual magma, rich in Be, Sn, W, B, F, H2O, migrate to the top of the granitic or surrounding rocks [72] (Figure 12). Due to the intrusion of pegmatite magma into the surrounding rock, there is a significant temperature difference between the two, resulting in undercooling of the melt-liquid phase and a decrease in the solubility of rare metal minerals [69]. Meanwhile, some minerals rich in flux components (such as fluorite, tourmaline, and apatite) crystallize [73], leading to a great decrease in contents of flux components such as F, B, and P in the ore-forming system and Beryl-I crystallization. It is also the reason why Beryl-I is mainly produced in the contact between albitite pegmatite and the surrounding rock.
With the continuous differentiation and crystallization of the leucogranitic magma, a large number of ore-forming materials such as Be, W and Sn are brought out. When the ore-bearing fluid migrates to the marble layer with pressure releasing, extensive water–rock reactions occur between the ore-bearing fluid and marble layer [72], and then the skarn be-W-Sn ore body is formed (Figure 12). Beryllium mainly occurs in skarn minerals such as vesuvianite, scapolite and garnet modified by late hydrothermal solution, in the form of isomorphisms [32].
Influenced by the east–west extensional structure, some nearly north–south normal faults and secondary faults were developed in the middle Miocene (16–15 Ma) [27]. Ore-forming fluids migrated along these faults, and with the decrease in temperature and the mixing of some hydrothermal water, the hydrothermal vein-type Be-Sn-W ore (chemical) bodies were formed [72] (Figure 12). It can be seen from the above discussion (Section 5.2) that beryl in hydrothermal vein-type ore bodies is related to hydrothermal activity. The microthermometry results of fluid inclusions show that the temperature of ore-forming fluid of hydrothermal vein orebodies ranges from 197 to 343 °C [28]. Under such mineralization conditions, beryllium mostly exists in the form of fluoride, chloride, and complexes containing fluorine or chlorine, while beryllium tends to migrate as beryllium carbonate complexes if the activity of fluorine in the solution is low and the activity of carbonic acid is high [74,75]. There are many wide fluorite veins in the Cuonadong mining area [17,28], indicating that the ore-forming fluid is extremely rich in fluorine. Therefore, it can be inferred that beryllium is most likely transported as beryllium fluoride (BeF) or fluoride complex ions (e.g., [BeF4]2−) [74], in this area. The result of beryl-II’s mapping shows that the development of these complex zonings is mainly caused by the differences among the contents of elements such as Na, Cs, Be, Al, Cr, and As (Figure 7). It probably reflects the unbalanced crystallization environment and the heterogeneity of ore-forming fluids, indicating that extensive fluid immiscibility and metasomatism has occurred in the metallogenic system [72], and it has undergone the superimposed transformation of later ore-forming fluids, such as strong alkaline metasomatism in the late period of beryl mineralization [10,76]. Field geological survey studies have shown that a variety of alterations, such as greisenization and albitization, often develop in the hydrothermal vein orebodies [77]. During the process, liquid H2O is bound to be generated, which cannot only reduce the temperature, but can also react with beryllium complexes in the ore-forming fluid and combine with Si and Al, thereby forming hydrothermal beryl. The probable reaction formula is as follows [77]:
mK Al Si 3 O 8   +   n H 2 O   +   p Na OH   +   q Na 2 Be F 4 a Na Al 3 O 8   +   bK Al 3 Si 3 O 10 OH 2   + ( K - feldspar )                                                   ( albite )           ( muscovite ) c Be 3 Al 2 Si 6 O 18   +   d Si O 2   +   e H 2 O   +   f F 2 ( beryl )       ( quartz )
Beryl occurrences in the pegmatites and hydrothermal vein ore bodies in the Cuonadong mining area indicate that beryllium mineralization is characterized by magmatic–hydrothermal transition. Some granitic rare metal deposits in South China have similar characteristics; that is, the pegmatic crusts usually develop on the top of the granitic rock mass, while hydrothermal mineralization occurs at the surroundings of the rock mass, and the hydrothermal mineralization pattern is closely related to the surrounding rocks [10,78,79]. A large number of leucogranites similar to the Cuonadong leucogranite are developed in the Himalayas, and the Tethys Himalayan sedimentary sequence is mainly a set of clastic rocks and carbonate rocks originating from passive Indian continental margin [32]. In conclusion, there is further great potential and prospects for the pegmatite, skarn and hydrothermal-vein rare metal deposits related to leucogranite in the Himalayas.

6. Conclusions

The results of this study allow for the following conclusions for beryls from the Cuonadong deposit.
(1)
In terms of composition, beryls in the Cuonadong deposit are alkaline beryls, among which beryl-I is Li-Cs beryl, and beryl-II consists of Na and Na-Li beryl. Structurally, they are t-beryls. This indicates that beryl has experienced magmatism to hydrothermal alkali-metasomatism in the late stage of pegmatitic magmatism during formation.
(2)
The mechanism of element substitution in Beryl-I includes (Na,Cs)Li□-1Be-1 channel-tetrahedral substitution, (Na,Cs)Fe2+-1Al-1 channel-octahedral substitution and NaCs-1 the mutual substitution of alkali metal ions in the “channel”, whereas that in Beryl-II consists of NaLi□-1Be-1 channel-tetrahedral substitution and Na(Fe2+,Mg)□-1Al-1 channel-octahedral substitution.
(3)
The precipitation of Beryl-I is mostly caused by the emplacement of highly fractionated magma containing Be to the top of the rock mass or surrounding rock, the melt-fluid undercooling, and the crystallization of volatile minerals (such as tourmaline and fluorite).
(4)
Beryl-II precipitates owe to the ore-forming fluid mixing with the hydrothermal water and cooling and large amounts of crystallization of volatile minerals (mainly fluorite).

Author Contributions

Y.W.: Investigation, Conceptualization, Supervision, Formal Analysis, Writing—Original Draft, Writing—Review and Editing. G.L.: Investigation, Formal Analysis, Funding Acquisition. W.L.: Investigation, Resources. Z.Z.: Formal Analysis, Data Curation. All authors have read and agreed to the published version of the manuscript.

Funding

This work was sponsored by the National Key R&D Program of China (2019YFC0605201) and the National Natural Science Foundation of China (91955208 and 42002097).

Data Availability Statement

All data derived from this research are presented in the enclosed figures and tables.

Acknowledgments

The authors appreciate two anonymous reviewers for the scientific suggestions and Zhen-Yu Chen at the Chinese Academy of Geological Sciences and Li-Min Zhou at the National Research Center for Geoanalysis for their laboratory assistance and thoughtful discussions. This work also benefited from helpful discussions with Yong Huang, Jian-Gang Fu, Anping Xiang, Cheng-Shi Qing, Jian-Yang Wu, Li Ying-Xu, Guo-Tao Ma, Lei Dong, and Ke Gao.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. London, D.; Evensen, J.M. Beryllium in silicic magmas and the origin of beryl-bearing pegmatites. Rev. Mineral. Geochem. 2002, 50, 445–485. [Google Scholar] [CrossRef]
  2. Evensen, J.M.; London, D. Experimental partitioning of Be, Cs, and other trace elements between cordierite and felsic melt, and the chemical signature of S-type granite. Contrib. Mineral. Petrol. 2003, 144, 739–757. [Google Scholar] [CrossRef]
  3. Černý, P. Alkali variations in pegmatitic beryls and their petrogenetic implications. Neues Jahrb. Mineral. Abh. 1975, 123, 198–212. [Google Scholar]
  4. Charoy, B. Beryllium speciation in evolved granitic magmas: Phosphates versus silicates. Eur. J. Mineral. 1999, 11, 135–148. [Google Scholar] [CrossRef]
  5. Evensen, J.M.; London, D.; Wendlandt, R.F. Solubility and stability of beryl in granitic melts. Am. Mineral. 1999, 84, 733–745. [Google Scholar] [CrossRef]
  6. London, D. Reading Pegmatites: Part 1—What Beryl Says. Rocks Min. 2015, 90, 138–153. [Google Scholar] [CrossRef]
  7. Barton, M.D.; Young, S. Non-pegmatitic Deposits of Beryllium: Mineralogy, Geology, Phase Equilibria and Origin. Rev. Mineral. Geochem. 2002, 50, 591–691. [Google Scholar] [CrossRef]
  8. Zhou, Q.F.; Qin, K.Z.; Tang, D.M.; Wang, C.L.; Ma, L.S. Mineralogical characteristics and significance of beryl from the rare-element pegmatites in the Lushi County, East Qinling, China. Acta Petrol. Sin. 2019, 35, 1999–2012. [Google Scholar]
  9. Černý, P.; Anderson, A.J.; Tomascak, P.B.; Chapman, R. Geochemical and morphological features of beryl from the Bikita granitic pegmatite, Zimbabwe. Can. Mineral. 2003, 41, 1003–1011. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, R.C.; Che, X.D.; Zhang, W.L.; Zhang, A.C.; Zhang, H. Geochemical evolution and late re-equilibration of Na-Cs-rich beryl from the Koktokay #3 pegmatite (Altai, NW China). Eur. J. Mineral. 2009, 21, 795–809. [Google Scholar]
  11. Zhou, Q.; Qin, K.; Tang, D.; Wang, C.; Tian, Y.; Sakyi, P.A. Mineralogy of the Koktokay No.3 Pegmatite, Altai, NW China: Implications for Evolution and Melt-Fluid Processes of Rare-metalpegmatites. Eur. J. Mineral. 2015, 27, 433–457. [Google Scholar] [CrossRef]
  12. Bragg, W.L.; West, J. The structure of beryl, Be3Al2Si6O18. Proc. R. Soc. London. Ser. A Contain. Pap. Math. Phys. Character. 1926, 111, 691–714. [Google Scholar]
  13. Gibbs, G.V.; Breck, D.W.; Meagher, E.P. Structural refinement of hydrous and anhydrous synthetic beryl, Al2(Be3Si6)O18 and emerald, Al1.9Cr0.1(Be3Si6)O18. Lithos 1968, 1, 275–285. [Google Scholar] [CrossRef]
  14. Li, G.M.; Zhang, L.K.; Jiao, Y.J.; Xia, X.B.; Dong, S.L.; Fu, J.G.; Liang, W.; Zhang, Z.; Wu, J.Y.; Dong, L.; et al. First discovery and implications of Cuonadong superlarge Be-W-Sn polymetallic deposit in Himalayan metallogenic belt, souther Tibet. Miner. Depos. 2017, 36, 1003–1008. [Google Scholar]
  15. Zhang, Z.; Zhang, L.K.; Li, G.M.; Liang, W.; Xia, X.B.; Fu, J.G.; Dong, S.L.; Ma, G.T. The Cuonadong Gneiss Dome of North Himalaya: A New Member of Gneiss Dome and a New Proposition for the Ore-controlling Role of North Himalaya Gneiss Domes. Acta Petrol. Sin. 2017, 38, 754–766. [Google Scholar]
  16. Liang, W.; Zhang, L.K.; Xia, X.B.; Ma, G.T.; Huang, Y.; Zhang, Z.; Fu, J.G.; Cao, H.W.; Miao, H.Q.; Li, G.M. Geology and preliminary mineral genesis of the Cuonadong W-Sn polymetallic deposit, Southern Tibet, China. Earth Sci. 2018, 43, 2742–2754. [Google Scholar]
  17. Xia, X.B.; Li, G.M.; Zhang, L.K.; Zhang, Z.; Cao, H.W.; Liang, W. Geological characteristics of and prospecting strategy for the Xianglin Be-Sn polymetallic ore deposit in the Cuonadong gneiss dome in southern Tibet. Earth Sci. Front. 2022, 29, 93–106. [Google Scholar]
  18. Gao, L.E.; Gao, J.H.; Zhao, L.H.; Hou, K.J.; Tang, S.H. The Miocene leucogranite in the Nariyongcuo gneiss dome, southern Tibet: Products from melting metapelite and fractional crystallization. Acta Petrol. Sin. 2017, 33, 2395–2411. [Google Scholar]
  19. Zhang, L.K.; Zhang, Z.; Li, G.M.; Dong, S.L.; Xia, X.B.; Liang, W.; Fu, J.G.; Cao, H.W. Rock Assemblage, Structural Characteristics and Genesis Mechanism of the Cuonadong Dome, Tethys Himalaya. Earth Sci. 2018, 43, 2664–2683. [Google Scholar]
  20. Dong, H.W.; Xu, Z.Q.; Meng, Y.K.; Yi, Z.Y. Geochronology of leucogranites in the Cuonadong dome, southern Tibet and limitation of the timing of the Southern Tibet etachment System (STDS). Acta Petrol. Sin. 2017, 33, 3741–3752. [Google Scholar]
  21. Fu, J.G.; Li, G.M.; Wang, G.H.; Zhang, L.K.; Liang, W.; Zhang, Z.; Dong, S.L.; Huang, Y. Timing of E-W extension deformation in North Himalaya: Evidences from ArAr age in the Cuonadong dome, South Tibet. Earth Sci. 2018, 43, 2638–2650. [Google Scholar]
  22. Fu, J.G.; Li, G.M.; Wang, G.H.; Zhang, L.K.; Liang, W.; Zhang, Z.; Zhang, X.Q.; Huang, Y. Synchronous granite intrusion and E–W extension in the Cuonadong dome, southern Tibet, China: Evidence from field observations and thermochronologic results. Int. J. Earth Sci. 2018, 107, 2023–2041. [Google Scholar] [CrossRef]
  23. Fu, J.; Li, G.; Wang, G.; Zhang, L.; Liang, W.; Zhang, X.; Jiao, Y.; Dong, S.; Huang, Y. Structural analysis of sheath folds and geochronology in the Cuonadong Dome, southern Tibet, China: New constraints on the timing of the South Tibetan Detachment System and its relationship to North Himalayan Gneiss Domes. Terra Nova 2020, 32, 300–323. [Google Scholar] [CrossRef]
  24. Xie, L.; Qiu, H.; Bai, X.; Zhang, W.; Wang, Q.; Xia, X. Geochronological and geochemical constraints on the Cuonadong leucogranite, eastern Himaaya. Acta Geochim. 2018, 37, 347–359. [Google Scholar] [CrossRef]
  25. Xia, X.B.; Li, G.M.; Cao, H.W.; Liang, W.; Fu, J.G. Petrogenic Age and Geochemical Characteristics of the Mother Rock of Skarn Type Ore Body in the Cuonadong Be-W-Sn Polymetallic Deposit, Southern Tibet. Earth Sci. 2019, 44, 2207–2223. [Google Scholar]
  26. Xia, X.B.; Xiang, A.P.; Li, G.M.; Zhang, L.K.; Cao, H.W.; Zhang, Z.; Liang, W. Neoproterozoic magmatism of Cuonadong dome and its tectonic implication, Tibet, China. J. Chengdu Univ. Technol. (Sci. Technol. Ed.) 2019, 46, 435–448. [Google Scholar]
  27. Cao, H.W.; Li, G.M.; Zhang, Z.; Zhang, L.K.; Dong, S.L.; Xia, X.B.; Liang, W.; Fu, J.G.; Huang, Y.; Xiang, A.P.; et al. Miocene Sn polymetallic mineralization in the Tethyan Himalaya, southeastern Tibet: A case study of the Cuonadong deposit. Ore Geol. Rev. 2020, 119, 103403. [Google Scholar] [CrossRef]
  28. Cao, H.W.; Li, G.M.; Zhang, R.Q.; Zhang, Y.H.; Zhang, L.K.; Dai, Z.W.; Zhang, Z.; Liang, W.; Dong, S.L.; Xia, X.B. Genesis of the Cuonadong tin polymetallic deposit in the Tethyan Himalaya: Evidence from geology, geochronology, fluid inclusions and multiple isotopes. Gondwana Res. 2021, 92, 72–101. [Google Scholar] [CrossRef]
  29. Liang, W.; Li, G.M.; Zhang, L.K.; Fu, J.G.; Huang, Y.; Zhang, Z. Cuonadong Be-rare polymetallic metal deposit: Constraints from Ar-Ar age of hydrothermal muscovite. Sediment. Geol. Tethyan Geol. 2020, 40, 76–81. [Google Scholar]
  30. Lin, B.; Tang, J.X.; Zheng, W.B.; Leng, Q.F.; Lin, X.; Wang, Y.Y.; Meng, Z.; Tang, P.; Ding, S.; Xu, Y.F.; et al. Geochemical characteristics, age and genesis of Cuonadong leucogranite, Tibet. Acta Petrol. Miner. 2016, 35, 391–406. [Google Scholar]
  31. Huang, C.M. Petrogenesis of the Cuonadong-Lhozag Leucogranites and Implication for Tectonic Evolution and Be-W-Sn Metallogeny in Southern Tibet. Ph.D. Thesis, China University of Geosciences, Beijing, China, 2019. [Google Scholar]
  32. He, C.T.; Qin, K.Z.; Li, J.X.; Zhou, Q.F.; Zhao, J.X.; Li, G.M. Preliminary study on occurrence status of beryllium and genetic mechanism in Cuonadong tungsten-tin-beryllium deposit, eastern Himalaya. Acta Petrol. Sin. 2020, 36, 3593–3606. [Google Scholar]
  33. Shen, J.Q.; Hu, Z.K.; Cui, S.Y.; Zhang, Y.F.; Li, E.Q.; Liang, W.; Xu, B. A Study on Beryl in the Cuonadong Be-W-Sn Polymetallic Deposit, Longzi County, Tibet, China. Crystals 2021, 11, 777. [Google Scholar] [CrossRef]
  34. Tao, X.Y.; Xie, L.; Wang, R.C.; Zhang, R.Q.; Hu, H.; Liu, C. Mineralogical characteristics of beryl: A case study of the beryls from Cuona and Qomolangma district in Himalaya. J. Nanjing Univ. (Nat. Sci.) 2020, 56, 815–829. [Google Scholar]
  35. Burchfiel, B.C.; Chen, Z.; Hodges, K.V.; Liu, Y.; Royden, L.H.; Deng, C.; Xu, J. The south Tibetan detachment system, Himalayan orogen: Extension contemporaneous with and parallel to shortening in a collisional mountain belt. Geol. Soc. Am. Spec. Pap. 1992, 269, 1–41. [Google Scholar]
  36. Pan, G.T.; Wang, L.Q.; Li, R.S.; Yuan, S.H.; Ji, W.H.; Yin, F.G.; Zhang, W.P.; Wang, B.D. Tectonic evolution of the Qinghai-Tibet Plateau. J. Asian Earth Sci. 2012, 53, 3–14. [Google Scholar] [CrossRef]
  37. Burg, J.-P.; Bouilhol, P. Timeline of the South-Tibet-Himalayan belt: The geochronological record of subduction, collision, and underthrusting from zircon and monazite U-Pb ages. Can. J. Earth Sci. 2019, 56, 1318–1332. [Google Scholar] [CrossRef]
  38. Yin, A.; Harrison, T.M. Geologic evolution of the Himalayan-Tibetan orogen. Annu. Rev. Earth Planet. Sci. 2000, 28, 211–280. [Google Scholar] [CrossRef] [Green Version]
  39. Yang, Z.S.; Hou, Z.Q.; Meng, X.J.; Liu, Y.C.; Fei, H.C.; Tian, S.H.; Li, Z.Q.; Gao, W. Post-collisional Sb and Au mineralization related to the South Tibetan detachment system, Himalayan orogen. Ore Geol. Rev. 2009, 36, 194–212. [Google Scholar] [CrossRef]
  40. Zhai, W.; Sun, X.; Yi, J.; Zhang, X.; Mo, R.; Zhou, F.; Wei, H.; Zeng, Q. Geology, geochemistry, and genesis of orogenic gold–antimony mineralization in the Himalayan Orogen, South Tibet, China. Ore Geol. Rev. 2014, 58, 68–90. [Google Scholar] [CrossRef]
  41. Hou, Z.Q.; Zhang, H.R. Geodynamics and metallogeny of the eastern Tethyan metallogenic domain. Ore Geol. Rev. 2015, 70, 346–384. [Google Scholar] [CrossRef]
  42. Sun, X.M.; Wei, H.X.; Zhai, W.; Shi, G.Y.; Liang, Y.H.; Mo, R.W.; Han, M.X.; Yi, J.Z.; Zhang, X.G. Fluid inclusion geochemistry and Ar–Ar geochronology of the Cenozoic Bangbu orogenic gold deposit, southern Tibet, China. Ore Geol. Rev. 2016, 74, 196–210. [Google Scholar] [CrossRef]
  43. Wang, D.; Zheng, Y.; Mathur, R.; Jiang, J.; Zhang, S.; Zhang, J.; Yu, M. Multiple mineralization events in the Zhaxikang Sb–Pb–Zn–Ag deposit and their relationship with geodynamic evolution in the North Himalayan Metallogenic Belt, South Tibet. Ore Geol. Rev. 2019, 105, 201–215. [Google Scholar] [CrossRef]
  44. Liu, C.; Wang, R.C.; Wu, F.Y.; Xie, L.; Liu, X.C.; Li, X.K.; Yang, L.; Li, X.J. Spodumene pegmatites from the Pusila pluton in the higher Himalaya, South Tibet: Lithium mineralization in a highly fractionated leucogranite batholith. Lithos 2020, 358–359, 105421. [Google Scholar] [CrossRef]
  45. Visonà, D.; Lombardo, B. Two mica and tourmaline leucogranites from the Everest⁃Makalu region (NepalTibet). Himalayan leucogranite genesis by isobaric heating? Lithos 2002, 62, 125–150. [Google Scholar] [CrossRef]
  46. Xie, L.; Tao, X.Y.; Wang, R.C. Highly fractionated leucogranites in the eastern Himalayan Cuonadong dome and related magmatic Be-Nb-Ta and hydrothermal Be-W-Sn mineralization. Lithos 2020, 354–355, 105286. [Google Scholar] [CrossRef]
  47. Wang, R.C.; Wu, F.Y.; Xie, L.; Liu, X.C.; Wang, J.M.; Yang, L.; Lai, W.; Liu, C. A preliminary study of rare-metal mineralization in the Himalayan leucogranite belts, South Tibet. China Earth Sci. 2017, 60, 1655–1663. [Google Scholar] [CrossRef]
  48. Wu, F.Y.; Liu, X.C.; Liu, Z.C.; Wang, R.C.; Xie, L.; Wang, J.M.; Ji, W.Q.; Yang, L.; Liu, C.; Khanal, G.P.; et al. Highly fractionated Himalayan leucogranites and associated rare-metal mineralization. Lithos 2020, 352–353, 105319. [Google Scholar] [CrossRef]
  49. Wu, F.Y.; Wang, R.C.; Liu, X.C.; Xie, L. New breakthroughs in the studies of Himalayan rare-metal mineralization. Acta Petrol. Sin. 2021, 37, 3261–3276. [Google Scholar]
  50. Zhang, Z.; Li, G.M.; Zhang, L.K.; Dong, S.L.; Liang, W.; Fu, J.G.; Huang, Y.; Cao, H.W.; Xia, X.B. The early Oligocene beryl-bearing pegmatite in the Cuonadong dome, southern Tibet: Its forming mechanism and geological significances. Sedimentary Geol. Tethyan Geol. 2020, 40, 14–30. [Google Scholar]
  51. Li, G.M.; Zhang, L.K.; Zhang, Z.; Xia, X.B.; Liang, W.; Hou, C.Q. New exploration progresses, resource potentials and prospecting targets of strategic minerals in the southern Qinghai-Tibet Plateau. Sedimentary Geol. Tethyan Geol. 2021, 41, 351–360. [Google Scholar]
  52. Tunheng, A.; Hirata, T. Development of signal smoothing device for precise elemental analysis using laser ablation-ICP-mass spectrometry. J. Anal. At. Spectrom. 2004, 19, 932–934. [Google Scholar] [CrossRef]
  53. Liu, S.B.; Wang, D.H.; Chen, Y.C.; Li, J.K.; Ying, L.J.; Xu, J.X.; Zeng, Z.L. 40Ar/39Ar ages of muscovite from different types tungsten-bearing quartz veins in the ChongYu-You concentrated mineral area in Gannan region and its geological significance. Acta Petrol. Sin. 2008, 82, 932–940. [Google Scholar]
  54. Liu, Y.; Hu, Z.; Gao, S.; Günther, D.; Xu, J.; Gao, C. In situ, analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard. Chem. Geol. 2008, 257, 34–43. [Google Scholar] [CrossRef]
  55. Paton, C.; Hellstrom, J.; Paul, B.; Woodhead, J.; Hergt, J. Iolite: Freeware for the visualisation and processing of mass spectrometric data. J. Anal. At. Spectrom. 2011, 26, 2508–2518. [Google Scholar] [CrossRef]
  56. Sardi, F.G.; Heimann, A. Pegmatitic beryl as indicator of melt evolution: Example from the Velasco district, Pampeana pegmatite province, Argentina, and review of worldwide occurrences. Can. Mineral. 2014, 52, 809–836. [Google Scholar] [CrossRef]
  57. Giuliani, G.; Cheilletz, A.; Zimmermann, J.L.; Cheiletz, A.; Arboleda, C.; Charoy, B.; Goget, P.; Frontan, F.; Giars, D. Fluid composition, δD of channel H2O, and δ18O of lattice oxygen in beryls: Genetic implications for Brazilian, Colombian, and Afghanistani emerald deposits. Int. Geol. Rev. 1997, 39, 400–424. [Google Scholar] [CrossRef]
  58. Zhao, T.X.; Chen, W.D.; Yin, T. EPMA quantitative analysis of beryllium mineral beryl. Acta Miner. Sin. 2019, 39, 169–175. [Google Scholar]
  59. Turner, D.; Groat, L.A.; Hart, C.J.R.; Mortensen, J.K.; Linnen, R.L.; Giuliani, G.; Wengzynowski, W. Mineralogical and geochemical study of the True Blue aquamarine showing, southern Yukon. Can. Mineral. 2007, 45, 203–227. [Google Scholar] [CrossRef]
  60. Hawthorne, F.C.; Černý, P. The alkali-metal positions in Cs-Li beryl. Can. Mineral. 1977, 15, 414–421. [Google Scholar]
  61. Sherriff, B.L.; Grundy, H.D.; Hartman, J.S.; Hawthorne, F.C.; Černý, P. The incorporation of alkalis in beryl: Multi-nuclear MAS NMR and crystal-structure study. Can. Mineral. 1991, 29, 271–285. [Google Scholar]
  62. Artioli, G.; Rinaldi, R.; Stahl, K.; Zanazzi, P.F. Structure refinements of beryl by single-crystal neutron and X-ray diffraction. Am. Mineral. 1993, 78, 762–768. [Google Scholar]
  63. Černý, P. Mineralogy of beryllium in granitic pegmatites. Mineral. Geochem. 2002, 50, 405–444. [Google Scholar] [CrossRef]
  64. Hawthorne, F.C.; Huminicki, D.M.C. The crystal chemistry of beryllium. Rev. Mineral. Geochem. 2002, 50, 333–403. [Google Scholar] [CrossRef]
  65. Lum, J.E.; Viljoen, F.; Cairncross, B.; Frei, D. Mineralogical and geochemical characteristics of BERYL (AQUAMARINE) from the Erongo Volcanic Complex, Namibia. J. Afr. Earth Sci. 2016, 124, 104–125. [Google Scholar] [CrossRef]
  66. Liu, Y.; Deng, J.; Sun, D.; Zhou, Y. Morphology and gensis typomorphism of minerals in W-Sn-Be deposit of Huya, Sichuan. Geosci. J. 2007, 32, 75–81. [Google Scholar]
  67. Uher, P.; Chudík, P.; Bacik, P. Beryl composition and evolution trends: An example from granitic pegmatites of the beryl⁃columbite subtype,Western Carpathians, Slovakia. J. Geosci. 2010, 55, 69–80. [Google Scholar] [CrossRef] [Green Version]
  68. Qin, K.Z.; Zhou, Q.F.; Tang, D.M.; Wang, R.C. Types, interal structure patterns, mineralization and prospects of rare-element pegmatites in East Qinling Mountain in comparison with feature of Chinese Altay. Miner. Depos. 2019, 38, 970–982. [Google Scholar]
  69. Tang, Y.; Qin, S.X.; Zhao, J.Y.; Lü, Z.H.; Liu, X.Q.; Wang, H.; Chen, J.Z.; Zhang, H. Solubility ofrare metals as a constraint on mineralization ofgranitic pegmatite. Earth Sci. Front. 2022, 29, 81–92. [Google Scholar]
  70. Zhao, K.H. Mineralogy and Metallogenic Characteristics of Beryl in Pingjiang, Hunan Province, China. Master’s Thesis, China University of Geosciences, Beijing, China, 2021. [Google Scholar]
  71. Xu, X.W.; Hong, T.; Li, H.; Ke, Q.; Chen, J.Z.; Liu, S.K.; Zhai, M.G. Concept of high-tempreture granite-pegmatite Li-Be metallogenic system with a primary study in the middle Altyn-Tagh. Acta Petrol. Sin. 2020, 36, 3572–3592. [Google Scholar]
  72. Dai, Z.W. Study on Mineralization of the Cuonadong Be-Sn-W Polymetallic Deposit, Tibet, China. Ph.D. Thesis, Chengdu University of Technology, Chengdu, China, 2020. [Google Scholar]
  73. Dai, Z.W.; Li, G.M.; Ding, J.; Zhang, L.K.; Cao, H.W.; Zhang, Z.; Liang, W. Chemical and Boron Isotopic Composition, and Significance of Tourmaline from the Cuonadong Tourmaline Granite, Tibet. Earth Sci. 2019, 44, 1849–1859. [Google Scholar]
  74. Wood, S.A. Theoretical prediction of speciation and solubility of beryllium in hydrothermal solution to 300 °C at saturated vapor-pressure-Application to Bertrandite Phenakite Deposits. Ore Geol. Rev. 1992, 7, 249–278. [Google Scholar] [CrossRef]
  75. Long, Z.Y.; Yu, X.Y.; Zheng, Y.Y. Ore formation of the Dayakou emerald deposit (Southwest China) constrained by chemical and boron isotopic composition of tourmaline. Ore Geol. Rev. 2021, 135, 104208. [Google Scholar] [CrossRef]
  76. Rao, C.; Wang, R.C.; Hu, H. Paragenetic assemblage of beryllium silicates and phosphates from the Nanping granite pegmatite dyke; Fujian province southeastern China. Can. Mineral. 2011, 49, 1175–1187. [Google Scholar] [CrossRef]
  77. Ge, Z.B. Emerald orebody in metaggranite-greisen beryl deposit. Geol. Explor. 1983, 6, 24–31. [Google Scholar]
  78. Rao, C.; Wang, R.C.; Hu, H. Electron-Microprobe Compositions and Genesis of Beryls from the Nanping No. 31 Granitic Pegmatite (Fujian Province, Southeastern China). Geol. J. China Univ. 2009, 15, 496–505. [Google Scholar]
  79. Wu, F.Y.; Liu, Z.C.; Liu, X.C. Himalayan leucogranite: Petrogenesis and implications to orogenesis and plateau uplift. Acta Petrol. Sin. 2015, 31, 1–36. [Google Scholar]
Figure 1. Map of the Himalayan orogenic belt showing the distribution of leucogranites, gneiss domes, tectonic zones and the Cuonadong deposit [27].
Figure 1. Map of the Himalayan orogenic belt showing the distribution of leucogranites, gneiss domes, tectonic zones and the Cuonadong deposit [27].
Minerals 12 00497 g001
Figure 2. Geological map of the Cuonadong dome [16].
Figure 2. Geological map of the Cuonadong dome [16].
Minerals 12 00497 g002
Figure 3. Geological map and sample location of Xianglin in the Cuonadong deposit.
Figure 3. Geological map and sample location of Xianglin in the Cuonadong deposit.
Minerals 12 00497 g003
Figure 4. Typical mineralization at Cuonadong deposit. (a): The contact zone between muscovite granite and pegmatite, in which a beryl crystal has a hexagonal section. (b): Skarn-type ores under ultraviolet light showing scheelite (distinct light blue color). (c): Ore-bearing hydrothermal veins interspersed in viritic rock and skarn. (d): In the greisen,the cassiserite-quartz vien was cut by the cas-siserite sulfide vein.
Figure 4. Typical mineralization at Cuonadong deposit. (a): The contact zone between muscovite granite and pegmatite, in which a beryl crystal has a hexagonal section. (b): Skarn-type ores under ultraviolet light showing scheelite (distinct light blue color). (c): Ore-bearing hydrothermal veins interspersed in viritic rock and skarn. (d): In the greisen,the cassiserite-quartz vien was cut by the cas-siserite sulfide vein.
Minerals 12 00497 g004
Figure 5. Occurrence characteristics of beryl in the Xianglin area. (a): Beryl-bearing pegmatite. (b): beryl co-existing with quartz, plagioclase, K-feldspar, muscovite and fluorite (plane-polarized light). (c): beryl co-existing with quartz, plagioclase, K-feldspar, muscovite (BSE image); (d) Beryl- and cassiterite-bearing hydrothermal vein. (e) Beryl with band structure, associated with cassiterite, quartz, muscovite and albite (plane-polarized light). (f): Beryl with complex zonal structures, associated with cassiterite, K-feldspar, quartz, muscovite, albite and scorodite (BSE image).
Figure 5. Occurrence characteristics of beryl in the Xianglin area. (a): Beryl-bearing pegmatite. (b): beryl co-existing with quartz, plagioclase, K-feldspar, muscovite and fluorite (plane-polarized light). (c): beryl co-existing with quartz, plagioclase, K-feldspar, muscovite (BSE image); (d) Beryl- and cassiterite-bearing hydrothermal vein. (e) Beryl with band structure, associated with cassiterite, quartz, muscovite and albite (plane-polarized light). (f): Beryl with complex zonal structures, associated with cassiterite, K-feldspar, quartz, muscovite, albite and scorodite (BSE image).
Minerals 12 00497 g005
Figure 6. Li(apfu) vs. Na + Cs + K(apfu) of the beryls from the Xianglin area.
Figure 6. Li(apfu) vs. Na + Cs + K(apfu) of the beryls from the Xianglin area.
Minerals 12 00497 g006
Figure 7. Backscattered electron image (a) and LA-ICP-MS mapping images (bl) of beryl-II from the hydrothermal vein.
Figure 7. Backscattered electron image (a) and LA-ICP-MS mapping images (bl) of beryl-II from the hydrothermal vein.
Minerals 12 00497 g007
Figure 8. Plot of R2+ versus A+. R2+ represents divalent cations in the Y (octahedral) site, including Fe, Mg; A+ represents monovalent alkali cations in the channel, including Na, Cs.
Figure 8. Plot of R2+ versus A+. R2+ represents divalent cations in the Y (octahedral) site, including Fe, Mg; A+ represents monovalent alkali cations in the channel, including Na, Cs.
Minerals 12 00497 g008
Figure 9. Plots of main cations in the X site, Y site and channel of beryl. (a): Li+ vs. Na+ plot. (b): Cs+ vs. Na+ plot. (c): Fe2+ vs. Mg2+ plot. (d): Fe2+ vs. Al3+ plot. (e): Mg2+ vs. Al3+ plot. (f): (Fe2+ + Mg2+) vs. Al3+ plot. Fe is one of the most important chromogenic ions of beryl. In addition to replacing Al3+ in the Y (octahedral) site, Fe appears in the channel, making it appear light green, blue–green, green–blue and light blue [56]. Emerald is chromatized by Fe, V and Cr [59], and a certain amount of Fe, V and Cr is detected in Beryl-I and Beryl-II. Thereby, there is likely to be a potential for further emerald prospecting at the Xianglin area.
Figure 9. Plots of main cations in the X site, Y site and channel of beryl. (a): Li+ vs. Na+ plot. (b): Cs+ vs. Na+ plot. (c): Fe2+ vs. Mg2+ plot. (d): Fe2+ vs. Al3+ plot. (e): Mg2+ vs. Al3+ plot. (f): (Fe2+ + Mg2+) vs. Al3+ plot. Fe is one of the most important chromogenic ions of beryl. In addition to replacing Al3+ in the Y (octahedral) site, Fe appears in the channel, making it appear light green, blue–green, green–blue and light blue [56]. Emerald is chromatized by Fe, V and Cr [59], and a certain amount of Fe, V and Cr is detected in Beryl-I and Beryl-II. Thereby, there is likely to be a potential for further emerald prospecting at the Xianglin area.
Minerals 12 00497 g009
Figure 10. Frequency histogram of alkali elements of beryl from the Xianglin area. (a):Alkali elements of Beryl-I. (b): Alkali elements of Beryl-II.
Figure 10. Frequency histogram of alkali elements of beryl from the Xianglin area. (a):Alkali elements of Beryl-I. (b): Alkali elements of Beryl-II.
Minerals 12 00497 g010
Figure 11. Mg/Fe vs. Cs/Na diagram of beryl from the Xianglin area.
Figure 11. Mg/Fe vs. Cs/Na diagram of beryl from the Xianglin area.
Minerals 12 00497 g011
Figure 12. Metallogenic model of the Cuonadong deposit, Tibet (after [72]).
Figure 12. Metallogenic model of the Cuonadong deposit, Tibet (after [72]).
Minerals 12 00497 g012
Table 1. Representative electronic microprobe analysis of beryls from the Xianglin area.
Table 1. Representative electronic microprobe analysis of beryls from the Xianglin area.
Elt (wt.%)Beryl-I n = 9Beryl-II n = 23
Max.Min.Avg.CND21-01CND21-02CND21-03Max.Min.Avg.XL20-03XL20-05XL20-08
SiO264.4263.3363.7263.4863.4064.1065.4163.5764.3365.1365.4164.25
Al2O317.6317.1917.3817.2617.6317.2917.4314.2415.9716.6116.7015.66
MgO0.030.000.010.030.020.012.061.071.621.371.321.66
CaO0.020.000.010.010.000.000.060.000.030.020.010.02
FeO0.590.260.410.480.260.461.310.170.610.230.250.67
MnO0.040.000.010.020.000.010.050.000.010.000.030.00
TiO20.060.000.020.000.000.000.040.000.010.000.010.03
Na2O0.900.680.810.840.680.901.661.251.431.331.321.25
K2O0.050.020.030.040.030.030.080.010.040.050.020.02
Cs2O1.770.311.011.291.770.420.390.000.030.050.000.00
Rb2O0.000.000.000.000.000.000.000.000.000.000.000.00
Li2O-LA0.870.740.780.760.870.750.560.040.310.050.050.49
Li2Ocal2.481.231.852.182.481.351.731.291.531.451.351.29
BeOcal12.0511.6911.8911.8811.6911.9813.5512.2712.8613.5113.5512.46
H2Ocal1.601.411.521.551.411.602.251.262.001.971.961.90
Total100.7596.1398.2498.1998.3398.04100.9798.0499.56100.51100.8498.60
apfu based on 18 oxygen atoms
Al3+1.971.921.941.931.971.921.901.601.751.801.801.73
Mg2+0.000.000.000.000.000.000.290.150.230.190.180.23
Fe2+0.050.020.030.040.020.040.100.010.050.020.020.05
Mn2+0.000.000.000.000.000.000.000.000.000.000.000.00
Ti2+0.000.000.000.000.000.000.000.000.000.000.000.00
∑O2.001.951.971.972.001.962.091.962.032.002.002.02
Ca2+0.000.000.000.000.000.000.010.000.000.000.000.00
Na+0.160.120.150.150.120.160.300.230.260.240.230.23
K+0.010.000.000.000.000.000.010.000.000.010.000.00
Cs+0.070.010.040.050.070.020.020.000.000.000.000.00
Rb+0.000.000.000.000.000.000.000.000.000.000.000.00
∑Ch0.210.180.190.210.200.180.310.230.270.250.240.23
Li+0.330.280.300.290.330.290.210.020.120.020.020.19
Be2+2.722.672.702.712.672.712.982.792.882.982.982.81
∑(T2)3.003.003.003.003.003.003.003.003.003.003.003.00
Si (T1)6.056.026.036.036.026.056.085.946.005.995.996.04
Table 2. Representative trace element compositions of beryls from the Xianglin area, as determined through LA-ICP-MS.
Table 2. Representative trace element compositions of beryls from the Xianglin area, as determined through LA-ICP-MS.
Let (ppm)Max.Min.Avg.Beryl-I n = 10Max.Min.Avg.Beryl-II n = 12
CND21-01CND21-02CND21-03XL20-05XL20-06XL20-07
Li4253.133429.603676.103519.224033.133495.121106.652291.52207.70207.701913.45215.32
Be53,207.2348,756.0550,687.4750,161.0251,617.3251,450.7752,146.6456,414.4949,970.4653,987.7150,709.4252,193.78
Rb215.31111.98167.14215.31141.68111.98147.24246.5341.68240.99139.4541.68
Cs16,734.255499.5410,339.285499.5412,632.125587.391182.122313.76472.711953.34695.12472.71
V3.180.000.780.000.943.1836.1565.953.3518.053.3565.95
Cr634.5813.63136.4820.2432.65243.24206.25642.6017.5399.3142.9659.71
Co4.960.000.500.000.000.000.645.060.000.000.001.03
Ni6.030.000.910.006.030.008.4265.860.0015.5565.860.00
Cu18.400.003.990.0015.780.0020.98111.860.000.0021.3936.21
Zn482.89231.49396.58463.16346.16231.49280.301781.6839.6339.63188.7469.41
Ga41.3824.3234.7237.3830.4139.0843.9794.9418.3594.9441.3140.44
As43.170.009.637.496.6732.091511.086432.480.00278.44917.086432.48
Sr1.770.000.201.770.000.002.5016.880.000.000.0016.88
Nb0.330.000.090.140.330.000.130.510.000.000.350.12
Ta1.020.000.140.110.010.000.110.270.000.020.000.10
W0.150.000.030.000.000.000.080.320.000.000.000.18
Bi0.120.000.030.040.000.042.229.170.000.252.409.17
Pb2.560.000.410.010.090.361.335.000.000.040.032.06
Th1.730.000.180.000.000.042.6216.460.000.591.7916.46
Table 3. In situ micro-X-ray diffraction (XRD) analysis results of beryls from the Xianglin area.
Table 3. In situ micro-X-ray diffraction (XRD) analysis results of beryls from the Xianglin area.
Samplea (Å)b (Å)c (Å)V (Å3)c/a
CND21-19.21039.21039.2209677.29201.0012
CND21-29.21419.21419.2243677.92901.0011
CND21-39.21369.21369.2231677.91801.0010
Avg.9.21279.21279.2228677.71301.0011
Max.9.21419.21419.2243677.92901.0012
Min.9.21039.21039.2209677.29201.0010
XL20-19.21409.21409.2180677.74101.0004
XL20-29.21119.21119.2202677.47501.0010
XL20-39.21399.21349.2206677.91301.0007
XL20-49.21169.21309.2200677.71901.0009
Avg.9.21279.21299.2197677.71201.0008
Max.9.21409.21409.2206677.91301.0010
Min.9.21119.21119.2180677.47501.0004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Li, G.; Liang, W.; Zhang, Z. The Chemical Characteristics and Metallogenic Mechanism of Beryl from Cuonadong Sn-W-Be Rare Polymetallic Deposit in Southern Tibet, China. Minerals 2022, 12, 497. https://doi.org/10.3390/min12050497

AMA Style

Wang Y, Li G, Liang W, Zhang Z. The Chemical Characteristics and Metallogenic Mechanism of Beryl from Cuonadong Sn-W-Be Rare Polymetallic Deposit in Southern Tibet, China. Minerals. 2022; 12(5):497. https://doi.org/10.3390/min12050497

Chicago/Turabian Style

Wang, Yiyun, Guangming Li, Wei Liang, and Zhi Zhang. 2022. "The Chemical Characteristics and Metallogenic Mechanism of Beryl from Cuonadong Sn-W-Be Rare Polymetallic Deposit in Southern Tibet, China" Minerals 12, no. 5: 497. https://doi.org/10.3390/min12050497

APA Style

Wang, Y., Li, G., Liang, W., & Zhang, Z. (2022). The Chemical Characteristics and Metallogenic Mechanism of Beryl from Cuonadong Sn-W-Be Rare Polymetallic Deposit in Southern Tibet, China. Minerals, 12(5), 497. https://doi.org/10.3390/min12050497

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop