Next Article in Journal
Co-Remediation of Acid Mine Drainage and Industrial Effluent Using Passive Permeable Reactive Barrier Pre-Treatment and Active Co-Bioremediation
Next Article in Special Issue
The Li-Bearing Pegmatites from the Pampean Pegmatite Province, Argentina: Metallogenesis and Resources
Previous Article in Journal
Dissolution of Magnetite and Hematite in Mixtures of Oxalic and Nitric Acid: Mechanisms and Kinetics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Minerals
Volume 12
Issue 5
10.3390/min12050564
minerals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geological, Mineralogical and Geochemical Study of the Aquamarine-Bearing Yamrang Pegmatite, Eastern Nepal with Implications for Exploration Targeting

by
Sushmita Bhandari
1,2,3,
Kezhang Qin
1,2,*,
Qifeng Zhou
4 and
Noreen J. Evans
5
1
Key Laboratory of Mineral Resources, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
3
Department of Mines and Geology, Ministry of Industry Commerce and Supplies, Government of Nepal, Kathmandu 44600, Nepal
4
Institute of Mineral Resources Research, China Metallurgical Geology Bureau, Beijing 101300, China
5
School of Earth and Planetary Sciences, John de Laeter Centre, Curtin University, Perth, WA 6845, Australia
*
Author to whom correspondence should be addressed.
Minerals 2022, 12(5), 564; https://doi.org/10.3390/min12050564
Submission received: 14 March 2022 / Revised: 21 April 2022 / Accepted: 21 April 2022 / Published: 30 April 2022
(This article belongs to the Special Issue Petrology and Mineralogy of Pegmatite Deposits)
Browse Figures
Versions Notes

Abstract

:
The Yamrang Pegmatite in the Ikhabu Pegmatite Field, Eastern Nepal is Nepal’s primary source of aquamarine. This paper reports detailed mineralogy and whole rock granite and pegmatite geochemistry, and major and trace element data for alkali feldspar and muscovite in order to classify the aquamarine-bearing Yamrang Pegmatite, elucidate beryl-saturation processes and evaluate potential geochemical exploration tools for beryl-pegmatites. Five internal mineralogical/textural zones were identified in the Yamrang Pegmatite; zone 1 (saccharoidal albite); zone 2 (blocky perthitic microcline); zone 3 (muscovite–microcline–quartz); zone 4 (beryl-quartz), and zone 5 (miarolitic cavities). Zones 1–4 represent the magmatic stage, while zone 5 formed during the hydrothermal stage of pegmatite genesis. Spectacular aquamarines are recovered from miarolitic zone 5, while beryl saturation is found in zones 3, 4, and 5. Based on beryllium (Be) content, Be partition among co-existing minerals at the magmatic stage is beryl > muscovite > tourmaline > alkali feldspar > quartz. In contrast, the sequence at the hydrothermal stage is beryl > muscovite > albite > tourmaline > quartz. The Be content in rock-forming minerals decreases from pegmatite margin to core, and tourmaline could have played a significant role in Be enrichment processes in the marginal pegmatite zone. High temperature, a low degree of fractionation, and the dominance of Be-compatible mineral phases such as muscovite, calcium-rich alkali feldspar and tourmaline resulted in beryl undersaturation in marginal zones. However, low temperature, high fractional crystallization, and low abundance of Be-compatible mineral phases resulted in beryl saturation in inner zones. The strongly peraluminous nature, low total REE content (<500 ppm), mineral assemblage of beryl, tourmaline, spessartine, columbite-tantalite, depletion of Ba, Nb, and enrichment of Pb, Rb, Cs in the primitive mantle normalized multi-element plots suggest that the beryl-bearing Yamrang Pegmatite corresponds to the LCT pegmatite family. Alkali feldspar with K/Rb values of 30–150, Rb ~3000 ppm, Cs >100 ppm, and muscovite, with K/Rb ranging 18–50, Rb ~6000 ppm, Cs > 500 ppm, and Ta > 65 ppm in inner zones (3–5), indicate that the Yamrang Pegmatite is an intermediate-fractionated, beryl-type rare-element (REL) pegmatite. It is probable that whole rock Be content of >10 ppm could be considered an exploration guide to beryl mineralization in the region.

1. Introduction

Beryl (Be3Al2Si6O18) is an important beryllium (Be) ore and a valuable gem. It is the most common Be mineral in the Earth’s crust, occurring in granites, pegmatites, and veins at magmatic and hydrothermal stages. Emerald and aquamarine are the most valuable gem varieties of beryl, formed by the addition of chromophores Cr3+/V3+ for green emerald and Fe3+ for blue aquamarine. Rare-element pegmatites (REL) are a major source of beryl and other rare metals. Accounting for 1–2% of all granitic pegmatites, they may contain significant amounts of beryl, lithium aluminosilicates, phosphates, except for apatite, oxides, except for magnetite or ilmenite, and other rare minerals [1]. The study and classification of pegmatites primarily focuses on their rare-element mineralogy with REL pegmatites widely classified into the LCT (enriched in Li, Cs, Ta) or NYF (enriched in Nb, Y, F) petrogenetic families [2,3]. LCT pegmatites are of great economic interest as they host several economic rare-metal minerals such as spodumene, lepidolite, pollucite, petalite, columbite-tantalite, beryl, cassiterite, etc. Besides being an important Be ore and valuable gem, beryl is an early rare-metal mineral, crystallizing from fertile magma, that may evolve to deposit Nb, Ta, Li, Cs minerals. Hence, establishing the geological, geochemical, and mineralogical characteristics of beryl pegmatite has ramifications for the exploration of other highly economic commodities.
We investigated the field geology, mineralogy, and whole rock chemistry of pegmatites and granites from Eastern Nepal in order to identify barren and Be-rich bodies. We then concentrated our detailed studies on beryl-mineralized pegmatites in the Ikhabu Pegmatite Field (IPF), focusing on the large Yamrang Pegmatite. The Yamrang Pegmatite is intruded on the lower part of the Higher Himalayan Sequence in the Taplejung district of Eastern Nepal (Figure 1).
For the first time, this study reports major and trace element data for alkali feldspar and muscovite from the Yamrang Pegmatite, and whole rock chemistry of barren and Be-rich pegmatites from the region. The Be partition sequence in rock-forming minerals is compared with experimental data. This contribution aims to classify the beryl-bearing pegmatites from the area, investigate the beryl-saturation process and evaluate potential geochemical exploration guides for beryl-pegmatites. Results will provide the first insights into the properties of REL pegmatites from the region and aid in rare-metal exploration in the Himalayan belt in general, and in similar tectonic settings elsewhere.

2. Regional Geological Setting

2.1. Geology

The Himalayan orogenic belt resulted from the Cenozoic Indo-Asian collision along the Yarlung–Zangbo Suture Zone and extends 2400 km from Namche Barwa in the east to Nanga Parbat in the west. It comprises the Sub-Himalayan sequence, Lesser Himalayan sequence (LHS), Higher Himalayan sequence (HHS), and Tethyan Himalayan sequence (THS) from south to north (Figure 1).
The Eastern Nepal Higher Himalaya (Tamor–Kanchenjunga region) is bounded by the Main Central Thrust (MCT) in the south and South Tibetan Detachment System (STDS) in the north, and is divided into three lithological units; (a) Junbesi Paragneiss, (b) Kanchenjunga Migmatite, and (c) Kanchenjunga Paragneiss (Sil-Bt quartzo-feldspathic paragneiss in Figure 1) formation from the bottom to top [4,5,6]. They are correlated with Formation I, II, and III of Le Fort (1975) [7], respectively.
The Precambrian to Cambrian Junbesi Paragneiss (Formation I) is a 4–7 km thick unit with St-Ky gneiss at the lower part and Sil-grade migmatitic gneiss, concordant leucosomes, and discrete leucogranite dikes and sills at the top [4,5,6] (Figure 1). The gradational boundary between Junbesi Paragneiss and Kanchenjunga Migmatite (Formation II) is marked by an increased proportion of Sil-migmatite and roughly coincides with the position of the High Himal Thrust (HHT 1) [4,8]. The Kanchenjunga Migmatite is a 5–18 km thick unit comprising a sillimanite ± cordierite-bearing migmatitic gneiss, and quartzite with subordinate calc-silicates metabasites [6]. Cordierite-bearing migmatites appear in the mid-portion of the Kanchenjunga Migmatite and are restricted to the north of HHT 2 [8,9] (Figure 1).
Minerals 12 00564 g001 550
Figure 1. Regional Geological Map of Tamor–Kanchenjunga region, Eastern Nepal, showing the Ikhabu Pegmatite Field (IPF) and sample locations (compiled after [4,5,6,8]. Lower data clusters show the Chyanthapu area, and upper data clusters show the Kanchenjunga area. Abbreviations: LHS, Lesser Himalayan Sequence; HHS, Higher Himalayan Sequence; LMCT, Lower Main Central Thrust; MCT, Main Central Thrust; HHT, High Himal Thrust (HHT-1, HHT of [4]; HHT-2, HHT of [5]); STDS, South Tibet Detachment System; Bt, biotite; Sil, sillimanite. Mineral abbreviations are after [10].
Figure 1. Regional Geological Map of Tamor–Kanchenjunga region, Eastern Nepal, showing the Ikhabu Pegmatite Field (IPF) and sample locations (compiled after [4,5,6,8]. Lower data clusters show the Chyanthapu area, and upper data clusters show the Kanchenjunga area. Abbreviations: LHS, Lesser Himalayan Sequence; HHS, Higher Himalayan Sequence; LMCT, Lower Main Central Thrust; MCT, Main Central Thrust; HHT, High Himal Thrust (HHT-1, HHT of [4]; HHT-2, HHT of [5]); STDS, South Tibet Detachment System; Bt, biotite; Sil, sillimanite. Mineral abbreviations are after [10].
Minerals 12 00564 g001
The Kangchenjunga Migmatite has a gradational boundary with the overlying Kanchenjunga Paragneiss (Formation III) (Schelling, 1992). The Kanchenjunga Paragneiss formation is a 3–6 km thick unit consisting of Bt-Sil-gneiss/schist, quartzite, and widespread intrusions of Cambrian–Ordovician granites (510–480 Ma, now orthogneiss) [7,11,12,13] and Miocene (23–12 Ma) leucogranite sills and dikes [14,15,16] (Figure 1 and Figure 2a–d).
The studied beryl-mineralized pegmatites belong to the Ikhabu Pegmatite Field of Eastern Nepal. They are intruded into the middle part of the Junbesi Paragneiss in the vicinity of the sillimanite-in isograd (Figure 1). The beryl pegmatites are characterized by a mineral assemblage of quartz, K-feldspar, albite, muscovite, tourmaline, garnet, beryl, and ± rutile or columbite-tantalite (Figure 2h). The host rock (kyanite to sillimanite grade, locally migmatitic gneiss; Figure 2f,g) strikes NW–SE and dips NE at 25°–35°. The studied barren pegmatites are from the lower part of the Junbesi Paragneiss, closer to the Main Central Thrust at the Ankhop–Chyanthapu area. These pegmatites vary from 10 cm to 7 m in thickness. These are common pegmatites with a mineral assemblage of quartz, K-feldspar, albite, biotite, muscovite, tourmaline, and garnet (Figure 2e).

2.2. Granitic Plutons

Two types of granitic plutons are identified in the Eastern Nepal Higher Himalayan Sequence: (1) pre-Himalayan granites (orthogneiss) (510–480 Ma) of Kanchenjunga Paragneiss (formation III) formed by Cambrian–Ordovician magmatism [7,11,12,13] and (2) Himalayan leucogranites (23–12 Ma) formed by Miocene magmatism [14,15,16]. Leucogranites are essentially seen above HHT 1 as thin dykes and sills. The proportion of leucogranites increases from 10 to 20% at 4000 m to >60% above 6000 m altitude and they comprise the Kanchenjunga and Jannu mountain peaks [4,17]. The leucogranites at the Mt. Kanchenjunga-Mt Jannu area are known as the Kanchenjunga–Jannu leucogranites [17]. The plutons studied in this work are from the Kanchenjunga Paragneiss formation, upper part of HHS around the Lonak–Kanchenjunga–Jannu area, about 35–40 km northeast of the IPF, and from the Junbesi Paragneiss formation, lower part of HHS around Chyanthapu, about 35–40 km southeast from IPF (Figure 1 and Figure 2).
The coarse-grained orthogneiss consists of quartz, augen-shaped K-feldspar, plagioclase, garnet, tourmaline, biotite, muscovite, and occasionally cordierite and sillimanite (Figure 2a,c). The foliation is defined by biotite, tourmaline, and augen feldspar. The leucogranite sills and dikes were intruded into paragneiss and orthogneiss (Figure 2b,c). Himalayan leucogranites are divided into four types; biotite-, two-mica (Ms > Bt)-, tourmaline-, and garnet-leucogranites, based on dominant mineralogy [18,19,20,21,22]. Fine to medium-grained tourmaline-leucogranite consisting of quartz, feldspar, muscovite and tourmaline is dominant in the Kanchenjunga–Jannu area (Figure 2d); however, garnet-leucogranite fragments were also found around the base camp at Mt. Kanchenjunga and Mt. Jannu. All orthogneiss samples were collected from the Kanchenjunga area, while leucogranite samples were from the Kanchenjunga–Jannu and Chyanthapu areas.

3. Geology of Yamrang Pegmatite

Aquamarines were first discovered in Nepal in 1934 [23] and they are now famous worldwide. The Ikhabu pegmatites had been previously mentioned as a source of aquamarine in Nepal [24,25,26]; however, no further information is available in the literature. The present work represents the first detailed geological information on the beryl (aquamarine)-bearing pegmatites from Nepal Himalaya.
The beryl-bearing Ikhabu pegmatites are intruded into Ky-Sil grade biotite-schist to migmatitic gneiss in the middle part of the Junbesi paragneiss formation (Figure 2 and Figure 3). They are tabular to lenticular dikes with an undifferentiated, homogeneous to well-zoned character, ranging in width from 1 to 27 m and with a visible length from a few meters to 130 m. They show angular and sharp contact with country rocks (Figure 2f,g) and exhibit little to no deformation, implying that they are younger than Himalayan metamorphism. The long axes of larger dikes frequently follow the strike of the foliation, with pinched and inflated edges. The major minerals are quartz, alkali feldspar, muscovite >> biotite, with tourmaline and garnet as common accessories. Rare-metal minerals include beryl, columbite–tantalite, rutile, zircon, monazite, and xenotime (Figure 3 and Figure 4). Big granitic bodies were not exposed in the area; however, migmatites and thin to thick leucogranite sills are frequent in the upper part of the Junbesi Paragneiss [4,5,8].
The large Yamrang Pegmatite is located at the NE aspect of Yamrang hill at Sangsabu village in the Taplejung district, Eastern Nepal (Figure 3). The hill is steep and covered with dense vegetation. The Yamrang Pegmatite-host rock contact was not exposed in the field, however, a small section of sharp contact was observed in the roof of the Yamrang Pegmatite manual mining site. The pegmatite intruded kyanite grade, medium-grained, highly deformed, strongly foliated, biotite-rich schistose gneiss and was exposed at 1840 m elevation. Its exposed dimension is 27 m wide and 130 m long (Figure 3). The Yamrang Pegmatite shows axial-symmetric zonation and can be differentiated into five mineralogical–textural zones (Figure 3c and Figure 4, Table 1). All zones have gradational contact except zone 5, which is overprinted on zones 2 and 3. Zones 1–4 represent the magmatic stage of the Yamrang Pegmatite, while zone 5 represents the hydrothermal stage.
Zone 1 is the marginal saccharoidal albite zone, with a 1.5–2 m thick saccharoidal albite and quartz mineral assemblage, containing subordinate muscovite and accessory fine tourmaline (Figure 3c and Figure 4c). Zone 2 is the blocky perthitic microcline zone, 4.5–7 m thick with a mineral assemblage of perthitic microcline + quartz ± platy albite with accessory garnet and tourmaline (Figure 4d). Zone 3 is the muscovite–microcline-quartz zone, 5–6 m thick and containing coarse perthitic-microcline, quartz, albite, muscovite, subordinate tourmaline and garnet, and rare mineral beryl, columbite-tantalite, and rutile (Figure 3b). Zone 4 is the core zone, measuring 1–1.5 m wide and consisting of a coarse beryl–quartz assemblage (Figure 4f). Zone 5 comprises a group of miarolitic cavities overprinted on zones 2 and 3, mainly along the eastern boundaries of zones 2 and 3, and 3 and 4 (Figure 3c and Figure 4d). Aquamarine, translucent white beryls associated with platy albite, euhedral smoky quartz, book muscovite, and fine-grained columbite-tantalite are found in the miarolitic cavities (Figure 4h,i).

4. Analytical Methods

Petrographic observation and visual determination of modal mineral abundance in the thin section were performed under a polarizing microscope. Based on petrographic observation, representative samples were carefully selected for further analysis. Thin sections were used for electron probe microanalysis (EPMA), major element determination and in-situ laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) trace element analysis of minerals.

4.1. Whole-Rock Major and Trace Element Analysis

We divided our whole rock samples into five clusters: beryl-bearing pegmatite (Ikhabu pegmatites), paragneiss (host rock of Ikhabu pegmatites), barren pegmatite, orthogneiss, and leucogranite. Fresh whole rock samples (11 beryl pegmatites, 10 barren pegmatites, 3 orthogneiss and 9 leucogranites) were selected for whole rock major and trace element analyses (Supplementary Table S1). Studying pegmatite geochemistry from representative bulk samples is challenging given the heterogeneous texture and large mineral crystals. However, when samples are chosen carefully and combined with mineral chemistry, we can still enhance our understanding of pegmatite chemistry. In the case of the zoned Yamrang Pegmatite, large samples (2–5 kg) from internal zones (zone 1-YRP-1, YRP-11, YRP-13; zone 2-YRP-2; zone 3-YRP-6, YRP-7, YRP-9; zone 5-IW-8) were selected, each representative of the texture and mineral assemblage/mineral distribution of that zone. Similarly, from the barren pegmatites and beryl-bearing Wairung (IW-1) and Sangsabu (IW-5, IW-6) pegmatites that lack internal zonation and contain giant crystals, large samples (2–5 kg) representative of texture and pegmatite mineral assemblage were selected.
Analyses were conducted at Beijing Createch Testing Technology Co. Ltd. The whole rock was crushed to 200 mesh using a tungsten carbide ball mill. Major element compositions were determined on fused glass disks using an X-ray fluorescence spectrometer (XRF-1800, SHIMADZU, Kyoto, Japan) with an analytical precision within ±2% for the oxides present in an abundance > 0.5 wt.% and within ±5% for oxides > 0.1 wt.%. Loss on ignition (LOI) was measured for each sample after heating to 1000 °C for three hours in a muffle furnace.
Trace element analysis of an aliquot of whole rocks (40 mg) was performed using inductively coupled plasma mass spectrometry (Agilent 7500, Santa Clara, CA, USA) after acid digestion of samples in a HF + HNO3 flux for 48 h at about 190 °C. The solution was evaporated to incipient dryness, dissolved in concentrated HNO3, and evaporated at 150 °C to dispel the fluorides. The samples were diluted to about 80 g for analysis after being re-dissolved in 30% HNO3 overnight. An internal standard solution containing Rh monitored signal drift during ion counting. Analytical results for USGS standards indicated that the uncertainties for most elements were within 5%. The analytical procedures for whole-rock major and trace elements were similar to those described by [27].

4.2. Electron Microprobe Analysis

Major element compositions and BSE images of alkali feldspar and muscovite in carbon-coated thin sections were acquired using an electron probe microanalyzer, JEOL JXA-8100 at the State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China. An accelerating voltage of 15 kV, beam current of 10 nA, counting time of 10–30 s on each peak, and a beam diameter of 5 μm were used to carry out routine quantitative analyses. Natural samples and synthetic oxides were used for standard calibration: albite for Na, fluorite for F, BaSO4 for Ba, apatite for P, diopside for Ca, Mg, and Si, SrSO4 for Sr, bustamite for Mn, K-feldspar for K, Al2O3 for Al, Fe2O3 for Fe, rutile for Ti, Rb glass for Rb and pollucite for Cs. Raw data reduction was undertaken using a routine ZAF procedure. All elements were analyzed at greater than 98% precision. The structural formula for feldspar was calculated based on eight oxygen atoms, and end members were calculated as Or, Ab, and An. The structural formula for the muscovite was calculated based on 24 atoms of O, OH, and F. Li2O calculation was done after [28], and H2O calculation was done after [29].

4.3. Laser Ablation Inductively Coupled Plasma Mass Spectrometry

The trace element content of the alkali feldspar and muscovite was determined after removing the carbon coating from the thin sections used for EPMA. The experiment was carried out on an Agilent 7500a Q-ICP-MS instrument (Agilent Technologies, Santa Clara, CA, USA) coupled to a 193 nm ArF excimer laser system (Geolas HD, Lambda Physik, Göttingen, Germany) or an Analyte G2 193 nm ArF excimer laser ablation system at the State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences. The approach was similar to [30,31], with isotopes measured in peak-hopping mode with a laser beam diameter of ca. 44 μm and a 5 Hz repetition rate. The laser energy density was ~6 J cm−2. Helium was employed as the ablation gas to improve the transporting efficiency of ablated aerosols. NIST SRM 610 reference glass was used as calibration material, and ARM-1 [32] and BCR-2G were analyzed for data quality control. Silicon (29Si) was used as an internal standard, using the Si value determined from EPMA analysis. The raw data were reduced using GLITTER [33]. For most trace elements (>0.10 μg/g), the accuracy was better than ±10%, with analytical precision (1 RSD) of ±10%.

5. Results

5.1. Whole-Rock Chemistry

The major and trace element results for whole rock analysis are presented in Table 2, Supplementary Table S1. Most of the studied rocks plot in the granite field. Some samples from the alkali feldspar-rich zones of Yamrang Pegmatite plot in quartz-monzonite and syenite areas (Figure 5a). Since alkali feldspar has a high K and Na content, the bulk composition falls in the syenite, rather than the granite, field. All rocks analyzed have a peralkaline and peraluminous nature and plot in the S-type field with ASI index > 1.1 (1.2–3.3) (Figure 5b). The total REE content of studied pegmatite and granite samples ranges from 1 to 161 ppm, while paragneiss samples contain 121–178 ppm total REE (Supplementary Table S1).
The beryl-bearing Ikhabu pegmatites have generally high total alkali contents (14 wt.%) and a high Fe (1.3 wt.%) and Be content (4–531 ppm). Other notable elements include Ca (2.2 wt.%), P (0.5 wt.%), Nb up to 60 ppm, and Ta up to 17 ppm (Table 2). The barren pegmatites have a relatively lower Be content (<4.5 ppm), and a high Nb content (73 ppm) with 13 ppm Ta, 2.4 wt.% Ca, 1.5 wt.% Fe and 0.5 wt.% P. The abundance of Be, Li, Cs, Ta, and Nb in the whole rocks is shown in Figure 6. Enrichments of Be, Cs, and Nb is noted in beryl-bearing Ikhabu pegmatites. Lithium is enriched in orthogneiss (average-34.3, range 12–72 ppm), leucogranite (average-23.3, range 0.7–66 ppm), and barren pegmatite (average-22, range 6–69 ppm) relative to the Ikhabu pegmatites (average-13.3, range-0.7–33 ppm) (Supplementary Table S1). All studied rocks show a higher Nb than Ta content (Figure 6). Among the granitic plutons, two samples of tourmaline-leucogranite (KGP-18 and KGP-20) are enriched in Be (up to 11 ppm), Nb (21 ppm) and Ta (8 ppm) (Figure 6). These leucogranites are also enriched in Fe (~1.2 wt.%) and Ca (~1.2 wt.%), with P (0.1 wt.%), Mg (0.3 wt.%), and Ti (0.1 wt.%). The Be content in other plutons is <3 ppm, although some are richer in Nb (~33 ppm) and Ta (12 ppm). The paragneiss host rocks of Ikhabu pegmatites are rich in Fe (2.3 wt.%), with Ca (1.6 wt.%), Mg (3.2 wt.%), Ti (0.5 wt.%), and P (~0.1 wt.%). The host rock of the Yamrang Pegmatite (YRP-14) contain ~2.4 ppm Be and ~12.5 ppm Nb with low Ta (1 ppm) (Table 2).
Chondrite-normalized REE distribution patterns for the pegmatite samples show LREE enrichments and both negative and positive Eu anomalies (Figure 7a,c). Positive Eu anomalies in some pegmatite samples were probably due to wall rock interaction or weak plagioclase fractionation; however, this does not affect the overall REE trend [34]. Beryl pegmatites have a gentler REE chondrite pattern with a flatter HREE pattern than barren pegmatites (Figure 7a,c). Some beryl pegmatite samples show a positive Tm anomaly.
Minerals 12 00564 g005 550
Figure 5. Plots of (a) total alkali (Na2O + K2O) vs. SiO2 after [35]. Due to bulk samples being taken from the alkali feldspar-rich zones of the Yamrang Pegmatite, some data plot in the syenite field. (b) A/NK vs. A/CNK after [36] of pegmatites and granites from the Higher Himalaya of Eastern Nepal.
Figure 5. Plots of (a) total alkali (Na2O + K2O) vs. SiO2 after [35]. Due to bulk samples being taken from the alkali feldspar-rich zones of the Yamrang Pegmatite, some data plot in the syenite field. (b) A/NK vs. A/CNK after [36] of pegmatites and granites from the Higher Himalaya of Eastern Nepal.
Minerals 12 00564 g005
Minerals 12 00564 g006 550
Figure 6. Concentration profiles of (a) Li-, Cs-, Nb-, Ta-; (b) Be- in pegmatites, gneiss, and granites from Eastern Nepal, Higher Himalaya.
Figure 6. Concentration profiles of (a) Li-, Cs-, Nb-, Ta-; (b) Be- in pegmatites, gneiss, and granites from Eastern Nepal, Higher Himalaya.
Minerals 12 00564 g006
Minerals 12 00564 g007 550
Figure 7. Chondrite-normalized REE patterns and primitive mantle-normalized multi-element patterns for (a,b) Be-rich pegmatite; (c,d) Barren pegmatites, and (e,f) Orthogneiss and leucogranites. Chondrite values are from [37], and primitive mantle values are from [38].
Figure 7. Chondrite-normalized REE patterns and primitive mantle-normalized multi-element patterns for (a,b) Be-rich pegmatite; (c,d) Barren pegmatites, and (e,f) Orthogneiss and leucogranites. Chondrite values are from [37], and primitive mantle values are from [38].
Minerals 12 00564 g007
Chondrite-normalized REE distribution patterns for the leucogranites and orthogneiss samples are similar, showing enrichment of light REE and negative Eu anomalies (Figure 7e). However, the orthogneiss and Chyanthapu leucogranites show a flatter HREE pattern than the Kanchenjunga leucogranites (Figure 7e).
The primitive mantle normalized multi-element diagrams for beryl-pegmatites, barren pegmatites, leucogranites, and orthogneiss show a depletion of Ba, Nb, Sr, and moderate to strong enrichment of Cs, Rb, U, K, Pb (Figure 7b,d,f). Beryl pegmatites have a steeper normalized multi-element pattern than barren pegmatites, leucogranites, and orthogneiss, with stronger element enrichments and depletions (Figure 7b).

5.2. Feldspar Chemistry

Feldspars were sampled and analyzed from all five zones of the Yamrang Pegmatite. Samples YRP-1A, YRP-11A, YRP-13A were from zone 1; YRP-2B, YRP-12A from zone 2, YRP-4, YRP-9A from zone 3, and YRP-7E was from zone 5. YRP-8B is the microcline sample from the core margin associated with quartz and beryl. We have assigned it to core zone 4. The feldspar in samples from zones 1 and 5 is albite, while zones 2, 3, and 4 contain perthitic microcline (Figure 8). Feldspar samples show a homogeneous texture in BSE images.
Endmembers were determined in terms of Ab, An, and Or using feldspar major element data from various zones. The representative major element contents are shown in Table 3, and a full set of data is provided in Supplementary Table S2. The alkali feldspars from zones 1 and 5 were albite (Na-alkali feldspar), and those from zones 2, 3, and 4 were perthitic microcline (K-alkali feldspar). Alkali feldspar from zone 1 has a higher An component (An9–12), while the innermost part (core) has the lowest An (An0.1–0.7) (Figure 9). Albite from zone 1 consists of Ab87–100, An0–12, Or0–2. The albite sample from zone 5 is Ab96–98, An2–4, Or0.5–1, whereas the microcline composition from zone 2 is Ab10–98, An0–2, Or1–91. Zone 3 feldspar is Ab10–99, An0–2, Or1–90 and zone 4 is Ab12–99, An0–1, Or1–88 (Figure 9).
The average trace element data for the alkali feldspar (excluding the albite phase from perthitic microcline) is presented in Table 4, and the full set of data is given in Supplementary Table S3. Alkali feldspar from the Yamrang Pegmatite consists of average contents of K-42,479 ppm (range 114–165,414 ppm), Rb-543 ppm (range 1–2925 ppm), Cs-46 ppm (range 0–323 ppm), P-293 ppm (range 22–906 ppm), Ga-11.4 ppm (range 1–30 ppm), Nb < 1 ppm, Ba-179 ppm (range 1–462 ppm), and Ta below the limit of detection. The K/Rb values range from 46 to 1265, with K/Cs, Rb/Sr and Nb/Ta values of 375–13,259, 0–2259, and 3–20, respectively (Figure 10, Supplementary Table S3). The Cs vs. K/Rb plot for alkali feldspar shows an overall negative correlation between K/Rb and Cs, and data reach the Be-pegmatite field (Figure 10a). The median K/Cs value for alkali feldspar decreases from zones 2 to 5, with zone 1 and 3 being similar. The median K/Rb value also decreases from zones 1 to 5, with zone 5 being slightly higher than zone 4. The median Nb/Ta and Rb/Sr increase from zone 1 to 4 but Rb/Sr decreases in zone 5 (Figure 10).
At the magmatic stage, the Be content in alkali feldspar roughly decreases from outer to inner zones (zone 1–4), with an average of 5.6 ppm (range 1–10 ppm) in zone 1, 5.3 ppm (range 4–7 ppm) in zone 2, 5.5 ppm (range 3.5–8.5 ppm) in zone 3, and approximately 4 ppm in zone 4 (Figure 10b, Supplementary Table S3). The miarolitic albite from zone 5 (representing the hydrothermal stage) has the highest Be content of 11–16 ppm, higher than that in the zones representing the magmatic stage (zones 1–4) (Figure 10b). The Be content is higher in alkali feldspar from zone 1 (up to 10 ppm), and this zone also has the highest An component (An9–12) among all magmatic zones (zone 1–4). Be in alkali feldspar shows a strong positive correlation with Cl, Cr, Ga (r > 0.7) and a strong negative correlation with V, Ni (r < −0.7).

5.3. Muscovite Chemistry

Muscovite was sampled from different zones of the Yamrang Pegmatite (Figure 8); however, muscovite was not obtained from the core zone 4. Sample YRP-10A represents the innermost part of zone 3, closest to the core zone. Sample YRP-13B was from zone 1, YRP-2B was from zone 2, YRP-6A and YRP-10A were from zone 3, and YRP-5B and YRP-7F were from zone 5. The zone 1 sample shows secondary muscovite overgrowing albite (Figure 8a), while other samples show primary muscovite. Muscovite has a homogeneous texture in BSE images.
Representative muscovite major and trace element data from different zones of the Yamrang Pegmatite are shown in Table 5 and Table 6. The full set of data are given in Supplementary Tables S4 and S5. The Mg-Li vs. Fe+Mg+Ti-AlVI diagram for mica [41] in Figure 11a shows the composition of muscovite from the Yamrang Pegmatite. In Figure 11b, the Fe/(Fe+Mg) vs. Al plot shows that all analyzed muscovite was primary.
The K/Rb value ranges from 18 to 140, K/Cs from 199 to 7900, Rb/Sr from 170 to 19,660, and Nb/Ta from 2 to 17 (Table 6). In Figure 12a, the Cs vs. K/Rb plot of muscovite shows an overall negative correlation between K/Rb and Cs, and data points are distributed between Mus-class and REL-class pegmatites. In Figure 12b, the Be vs. K/Rb plot of muscovite shows few data points falling in the REL-class pegmatite field. Muscovites from the miarolitic zone are the most enriched in Be and Cs. The median Be content of muscovite fluctuates less in all zones, ranging from 23 to 30 ppm (Figure 12c). The median muscovite K/Rb and K/Cs value decreases from zones 2 to 5 (Figure 12d,e). The Nb/Ta ratio is highest at zone 3, corresponding to columbite-tantalite mineralization (Figure 12f).

6. Discussion

6.1. Classification of the Yamrang Pegmatite

Geoscientists have made many attempts to find a suitable classification scheme for pegmatites [2,3,24,39,44,45,46,47]. Ginsburg et al. (1979) [45] divided pegmatites into four categories based on the depth of emplacement, relationship to metamorphism, and granitic plutons: abyssal, muscovite, rare-element (REL) and miarolitic classes. Later this was revised into five classes: abyssal, muscovite, muscovite-rare-element, rare-element, and miarolitic, along with the introduction of two distinct petrogenetic families (LCT and NYF) by Černý and Ercit (2005) [2,3], which are widely used today (Figure 13).
Economically important, rare-element class pegmatites containing an economic concentration of lithophile rare elements are divided into two subclasses: REL-REE, corresponding to the NYF family, and REL-Li, corresponding to the LCT family. REL-REE subclass pegmatites are derived from post to anorogenic metaluminous to peraluminous granites with allanite, monazite, euxenite, and gadolinite as distinctive minerals. In contrast, REE-Li pegmatites are derived chiefly from syn to late orogenic peraluminous granites and intrude upper-greenschist to amphibolite facies host rocks (Figure 13) with beryl, columbite, tantalite, spodumene, petalite as distinctive minerals [2,3]. The Yamrang Pegmatite intruded into amphibolite facies Ky-Sil grade schistose gneiss, and contains beryl, with columbite-tantalite in inner zones (3–5). These features are consistent with the REL class, beryl type, beryl-columbite subtype pegmatite (Figure 13). The alkali feldspar, having K/Rb 30–150, Rb ~3000 ppm and Cs > 100 ppm, and muscovite having K/Rb 18–50, Rb ~6000 ppm, Cs > 500 ppm and Ta > 65 ppm indicates that the Yamrang Pegmatite is a low to intermediately fractionated REL pegmatite (Figure 10 and Figure 12) [48]. Barton and Young (2002) [49] divided Be-mineralization based on alumina and silica saturation into four types: strongly to weakly peraluminous, metaluminous to weakly peraluminous, peralkaline to metaluminous quartz-saturated, and silica-undersaturated peralkaline types. According to their classification, the pegmatitic Be-mineralization of the Ikhabu pegmatite Field with a high ASI index > 1.2 can be classified as a strong to weak peraluminous type.
Černý (1991) and Černý and Ercit (2005) [2,3] split REL class pegmatites into two distinct petrogenetic families: LCT and NYF types, based on trace element signatures. The LCT family is related to S-type sources in an orogenic setting with a peraluminous nature, abundant rare alkalis, Be, Bo, P and Sn, with a low total REE content (<500 ppm), and a low Nb/Ta ratio. In contrast, the NYF family is related to A-type sources commonly in an anorogenic setting, with a subaluminous to metaluminous nature, enriched HREE and F, and high Nb/Ta ratios [1,2,3]. The LCT and NYF families were proposed using trace-element signatures expressed by exotic minerals such as spodumene, pollucite, tantalite, pyrochlore, etc., which are only exhibited by the most fractionated pegmatites [2,3,50,51]. Hence, the LCT pegmatites are expected to be enriched in distinctive elements, Li, Cs, and Ta; however, this is not always the case with less fractionated pegmatites. Besides the enrichment of distinctive minerals, several other characteristics are used to classify the pegmatite into the LCT family, such as peraluminous nature, high Be, Rb, Sn, Ga, B, P, F, low total REE content (<500 ppm), low U, Y, Zr, Ti, a mineral assemblage of abundant beryl, tourmaline, muscovite, spessartine, minor presence of topaz, apatite and also association with an orogenic setting [3,50,52,53].
The whole rock geochemistry of the studied pegmatites and granites shows a strong peraluminous nature (ASI index of 1.2–3.3) and plots in the S-type field (Figure 5). The total REE content of the studied pegmatite and granite samples is very low, ranging from 1 to 161 ppm (Supplementary Table S1). The REE chondrite-normalized patterns show a negative slope with a negative Eu anomaly, although some pegmatite samples show a positive Eu anomaly. The primitive mantle normalized multi-element diagrams for the Yamrang Pegmatite show strong enrichment of Cs, Rb, Pb, U, K, and depletion of Ba, Nb, Sr, indicating a continental crustal origin. The Yamrang Pegmatite displays characteristics of low to intermediate fractionation, and does not exhibit an enrichment of distinctive minerals (Li, Cs, Ta) typical of LCT pegmatite; however, its strongly peraluminous nature, low total REE content, mineral assemblage of beryl, tourmaline, spessartine and columbite–tantalite, depletion of Ba and Nb, and enrichment of Pb, Rb and Cs in primitive mantle normalized multi-element plots suggest that the beryl-bearing Yamrang Pegmatite corresponds to the LCT pegmatite family.
This study also shows that the studied orthogneiss and leucogranite are S-type plutons. The S-type source of the orthogneiss and leucogranite agrees with previous research on these Himalaya plutons [54,55,56]. The S-type leucogranites with Be enrichment (up to 8–11 ppm) support the possibility of rare-metal mineralization prospects in the terrain, as discovered in other parts of the Himalayas [20,57,58,59,60].

6.2. Beryl Saturation Process in Yamrang Pegmatite

6.2.1. Magmatic Stage

The rare-element pegmatite whole-rock Be contents are reported in the range of ~32–575 ppm [61,62]. The whole-rock sample from the zone 1 of the Yamrang Pegmatite is enriched in Be (up to 531 ppm; sample YRP-11), but is not beryl saturated with a mineral assemblage of albite (60%), quartz (35%), with subordinate muscovite and tourmaline (5%). The Be content in quartz reported from worldwide REL pegmatites is the lowest (average, 0.1–0.44 ppm) [63,64] of all rock-forming minerals discussed in this paper. Hence, we assigned quartz to the bottom of the Be partition sequence. Albite from zone 1 has a relatively higher calcium component (An9–12) and a higher Be content (~10.4 ppm) than alkali feldspar from other zones (2–4). In addition, the Be content of muscovite (~31 ppm) and tourmaline (~15 ppm) from this zone is relatively higher than from later zones (zones 3, 4) (Be contents of beryl and tourmaline discussed in this paper are taken from our unpublished data). Accordingly, the Be partition sequence as determined by the Be content in coexisting minerals from zone 1 is muscovite > tourmaline > albite > quartz (Figure 14). The Be content in muscovite from zone 2 is (~28 ppm), with tourmaline and perthitic microcline having ~11 ppm and ~7 ppm, respectively. This yields a Be partition sequence of muscovite > tourmaline > microcline > quartz for this zone (Figure 14). Experimental models of beryl saturation show that at high temperatures (~700 °C), 99% of the Be in the melt can be accommodated by rock-forming minerals [65]. While beryl may not crystallize from a melt with 140 ppm Be at 700 °C, at a lower temperature (~400 °C), even 40 ppm Be in the melt can trigger beryl crystallization [65]. Beryl undersaturation in zones 1 and 2 probably occurred due to high-temperature conditions at the early pegmatite intrusion and consolidation stage. Additionally, Be is most compatible in plagioclase with a composition of An31, and Be compatibility decreases with digression from An30, on either side [65]. In agreement with this study, alkali feldspar enriched in the calcic component (An9–12) from marginal zone 1 had a higher Be content (up to 10 ppm) than the inner magmatic zones (zones 2–4). The ASI index of marginal zone 1 is lower (1.7) than that of inner zones (2.2), which may have increased the solubility of Be in the melt, favoring retention in the melt over partitioning into crystals under the higher temperature conditions [66].
Zone 3 is the first zone of magmatic beryl saturation, where beryl coexists with quartz and K-feldspar (microcline). Beryl appears during the late magmatic stages in zone 3 and exhibits a single episode of beryl saturation. This zone has low whole-rock Be content (5–31 ppm), with 4–9 ppm Be in microcline and 4–10 ppm Be in tourmaline, but the zone 3 muscovite has a similar Be content (26–33 ppm) as muscovite in other zones. The Be partition sequence for zone 3 is beryl > muscovite > tourmaline > microcline > quartz (Figure 14). In zone 4, beryl coexists with K-feldspar and quartz at the core margin and forms a beryl and quartz assemblage at the center. The microcline from the zone 4 has the lowest Be content (~4 ppm) among all zones, while muscovite is enriched in Be (21–47 ppm) relative to other zones. It should be noted that microcline and muscovite samples were not precisely taken from the core but were as close to the core as possible. The Be partition sequence for zone 4 is beryl > muscovite > microcline > quartz (Figure 14).
Evensen and London (2002) calculated Be partition coefficients for a variety of rock-forming minerals, demonstrating that Be is compatible in the order of cordierite >> calcic oligoclase (An31) > white mica > dark mica > albite > alkali feldspar ≅ quartz [65]. Similarly, Hezel et al. (2011) studied the Be partition sequence among cogenetic minerals in natural samples based on their Be concentration showing the order as staurolite > plagioclase > muscovite > tourmaline > biotite > K-feldspar > garnet [67]. The Be partition sequences among coexisting minerals at the Yamrang Pegmatite magmatic stage is in agreement with these experimental and natural sample data.
Low temperatures were likely prevalent during the late pegmatite crystallization stage, significantly increasing beryl saturation and many other factors likely contributed to beryl saturation in zones 3–5. As crystallization progressed towards the core of the pegmatite, calcic alkali feldspars were deposited in the outer zones, and incompatible elements and fluxes were concentrated in the residual melt. Hence, Be was saturated in the inner zones. In addition, interior zones with a low abundance of minerals compatible with Be, such as calcic albite and tourmaline, may have meant the melt became increasingly enriched in Be, facilitating beryl saturation late in pegmatite crystallization.

6.2.2. Hydrothermal Stage

Miarolitic cavities are empirical evidence of aqueous fluid and record fluid saturation conditions during pegmatite crystallization [50,51,68,69]. Texturally late, miarolitic cavities (designated as zone 5) represent the hydrothermal stage of Yamrang Pegmatite (Figure 14). Zone 5 is not geometrically conformable with previous magmatic zones; rather, it comprises a group of pockets overprinted on magmatic zones along the core margin and along the boundary of zone 2 and 3 in the east portion. Not all pockets have the same mineral assemblages or evolution stages. Beryl-bearing pockets show the highest stage of evolution with a typical mineral assemblage of albite-quartz-muscovite-tourmaline-beryl, which was also obtained experimentally from vapor precipitation [70]. The Be partition sequence among the coexisting minerals from zone 5 is beryl > muscovite > albite > tourmaline > quartz (Figure 14). In contrast to the magmatic zones, we observed Be to be more compatible in albite than in tourmaline in the miarolitic cavities. The pocket-forming P-T conditions in the pegmatite system are constrained as 2–3 kbar and 550 °C to as low as 300 °C [45,50,51], suggesting low-temperature conditions during crystallization of zone 5. Low temperatures can induce beryl saturation at a low concentration of Be in the melt [65]; however, the whole-rock Be content is high (114 ppm) in the miarolitic zone. This suggests Be was enriched in the fluid during the late stages of magmatic crystallization, although Be is less compatible in the vapor phase than in silicate melt with a melt/vapor partition coefficient of 0.2 at 650 °C [70]. A high concentration of Be in the fluid and low-temperature conditions probably facilitated beryl saturation in miarolitic cavities.

6.3. Exploration Implications

Whole rock samples of barren pegmatite yielded < 4.5 ppm of Be, while beryl-pegmatites predominantly had >10 ppm Be (Figure 6), which agrees with the range of Be concentration in fertile granites as listed in [48] (1–604 ppm). The chondrite-normalized REE patterns of beryl-bearing pegmatites have a gentler negative slope than barren pegmatites. The HREE pattern of beryl-pegmatites is flatter than that of barren pegmatites, likely due to the fractionation of garnets. The primitive mantle normalized multi-element diagrams for beryl-pegmatites and barren pegmatites look similar, however, the pattern for beryl pegmatites is steeper with stronger enrichment and depletion of elements, and barren pegmatites show more extreme Ti depletion than beryl pegmatites (Figure 7).
The median K/Rb, K/Cs, and Be content of alkali feldspar and muscovite roughly decreases from outer to inner zones in the Yamrang Pegmatite, reflecting increasing fractionation [2,48,50,71]. However, it should be noted that several fractionation indices break their trend in zone 2 (Figure 10 and Figure 12) and further research will investigate this feature.
The Be content of alkali feldspar ranges from 1 to 16 ppm in the Yamrang Pegmatite. A similar range of alkali feldspar Be content is reported from the Koktokay No. 3 (2–18 ppm) and Luumäki gem beryl pegmatites (7–24 ppm) [72,73]. The median Be content of muscovite shows little variation within zones, ranging from 23 to 30 ppm, similar to the Be concentration range in muscovite from REL-pegmatites such as the Koktokay No. 3 (8–71 ppm), Kenticha (11–38 ppm), and Renli pegmatites (10–23 ppm) [72,74,75]. While the ~10 ppm of Be in alkali feldspar and ~30 ppm of Be muscovite could be a useful indicator of proximal beryl mineralization, the lack of data on Be content in barren pegmatite feldspar and mica make this supposition inconclusive. Nonetheless, whole-rock Be contents of >10 ppm may serve as a potential exploration vector to beryl mineralization in the region. In addition, knowledge of the beryl concentration in zones 3, 4, and 5 may help refine selective mining of the Yamrang Pegmatite.
Beryllium shows a strong positive correlation with Cl, Cr, Ga (r > 0.7) and a strong negative correlation with V, Ni (r < −0.7) in alkali feldspar. These elements in alkali feldspar can also be potential indicators of Be enrichment. We also observed that tourmaline has a significant Be concentration (up to 15 ppm), similar to alkali feldspar (1–16 ppm) in the Yamrang Pegmatite. Tourmaline comprises a significant volume in the marginal zones in pegmatites and likely plays an important role in the Be enrichment process.
Future work should focus on comparing the mineral chemistry of the Yamrang Pegmatite with barren pegmatites and Be-rich leucogranites. More data on mineral geochemistry from different zones will help elucidate Be-enrichment processes during the magmatic and hydrothermal stages.

7. Conclusions

The conclusions of this study can be summarized as follows:
  • The Yamrang Pegmatite can be classified as an intermediate-fractionated REL class, beryl type, beryl-columbite subtype pegmatite of the LCT family.
  • High temperature, low fractionation, dominance of Be-compatible mineral phases such as muscovite, calcium-rich alkali feldspar, and tourmaline resulted in beryl undersaturated marginal zones. In contrast, low temperature, high fractional crystallization, and low abundance of Be-compatible mineral phases resulted in beryl saturated inner zones.
  • The Be partition sequence among coexisting minerals at the magmatic stage is beryl > muscovite > tourmaline > alkali feldspar > quartz, whereas at the hydrothermal stage, the sequence is beryl > muscovite > albite > tourmaline > quartz.
  • A whole rock Be content of >10 ppm may indicate beryl mineralization in the region.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/min12050564/s1, Table S1: The whole-rock major and trace element composition of pegmatites, gneisses and granites, Eastern Nepal; Table S2: The EPMA major element composition of alkali feldspar from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal; Table S3: The LA-ICP-MS trace element composition of alkali feldspar from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal; Table S4: The EPMA major element composition of muscovite from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal; Table S5: The LA-ICP-MS trace element composition of muscovite from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.

Author Contributions

Conceptualization, S.B. and Q.Z.; methodology, S.B. and Q.Z.; software, formal analysis, data curation, S.B.; design and resources, K.Q.; writing—original draft preparation, S.B.; writing—review and editing, S.B., K.Q., Q.Z. and N.J.E.; supervision, K.Q. funding acquisition, K.Q. All authors have read and agreed to the published version of the manuscript.

Funding

The first author was awarded with the CAS-TWAS President’s Fellowship. This research was funded by the Second Tibetan Plateau Scientific Expedition and Research (2019QZKK0806) and the Key Deployment Project of CAS (ZDRW-ZS-2020-4-1).

Data Availability Statement

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

Acknowledgments

The first author thanks profusely to family and friends; Gezahegn Abdissa Bayessa, Kabita Karki, Sabina Koirala, and Suman Panday for their enormous support, love, and motivation to overcome hardship and carry on with work during the COVID-19 Pandemic. We acknowledge Baburam Gyawali, Diwakar Lamichhane, He Changtong, Mahendra Maharjan, and Sun Zhenghao for fieldwork assistance. We are grateful to the editors and anonymous reviewers for their constructive comments. Finally, we recognize the Department of Mines and Geology of Nepal for providing continuous administrative support and study leave to the first author to conduct this research work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. London, D. A Petrologic Assessment of Internal Zonation in Granitic Pegmatites. Lithos 2014, 184–187, 74–104. [Google Scholar] [CrossRef]
  2. Černý, P. Rare-Element Granitic Pegmatites. Part I: Anatomy and Internal Evolution of Pegmatite Deposits. Geosci. Can. 1991, 18, 49–67. [Google Scholar]
  3. Černý, P.; Ercit, T.S. The Classification of Granitic Pegmatites Revisited. Can. Mineral. 2005, 43, 2005–2026. [Google Scholar] [CrossRef] [Green Version]
  4. Goscombe, B.; Gray, D.; Hand, M. Crustal Architecture of the Himalayan Metamorphic Front in Eastern Nepal. Gondwana Res. 2006, 10, 232–255. [Google Scholar] [CrossRef]
  5. Imayama, T.; Takeshita, T.; Arita, K. Metamorphic P-T Profile and P-T Path Discontinuity across the Far-Eastern Nepal Himalaya: Investigation of Channel Flow Models. J. Metamorph. Geol. 2010, 28, 527–549. [Google Scholar] [CrossRef]
  6. Schelling, D. The Tectonostratigraphy and Structure of the Eastern Nepal Himalaya. Tectonics 1992, 11, 925–943. [Google Scholar] [CrossRef]
  7. Le Fort, P. Himalayas: The Collided Range. Present Knowledge of the Continental Arc. Am. J. Sci. 1975, 275-A, 1–44. [Google Scholar]
  8. Larson, K.P.; Camacho, A.; Cottle, J.M.; Coutand, I.; Buckingham, H.M.; Ambrose, T.K.; Rai, S.M. Cooling, Exhumation, and Kinematics of the Kanchenjunga Himal, Far East Nepal. Tectonics 2017, 36, 1037–1052. [Google Scholar] [CrossRef]
  9. Imayama, T.; Takeshita, T.; Yi, K.; Cho, D.-L.; Kitajima, K.; Tsutsumi, Y.; Kayama, M.; Nishido, H.; Okumura, T.; Yagi, K.; et al. Two-Stage Partial Melting and Contrasting Cooling History within the Higher Himalayan Crystalline Sequence in the Far-Eastern Nepal Himalaya. Lithos 2012, 134–135, 1–22. [Google Scholar] [CrossRef]
  10. Whitney, D.L.; Evans, B.W. Abbreviations for Names of Rock-Forming Minerals. Am. Mineral. 2010, 95, 185–187. [Google Scholar] [CrossRef]
  11. DeCelles, P.G.; Gehrels, G.E.; Quade, J.; LaReau, B.; Spurlin, M. Tectonic Implications of U-Pb Zircon Ages of the Himalayan Orogenic Belt in Nepal. Science 2000, 288, 497–499. [Google Scholar] [CrossRef] [PubMed]
  12. Hodges, K.V. Tectonics of the Himalaya and Southern Tibet from Two Perspectives. GSA Bull. 2000, 112, 324–350. [Google Scholar] [CrossRef]
  13. Parrish, R.R.; Hodges, V. Isotopic Constraints on the Age and Provenance of the Lesser and Greater Himalayan Sequences, Nepalese Himalaya. Geol. Soc. Am. Bull. 1996, 108, 904–911. [Google Scholar] [CrossRef]
  14. Harrison, M.T.; Grove, M.; Mckeegan, K.D.; Coath, C.D.; Lovera, O.M.; Fort, P.L. Origin and Episodic Emplacement of the Manaslu Intrusive Complex, Central Himalaya. J. Petrol. 1999, 40, 3–19. [Google Scholar] [CrossRef]
  15. Searle, M.P.; Simpson, R.L.; Law, R.D.; Parrish, R.R.; Waters, D.J. The Structural Geometry, Metamorphic and Magmatic Evolution of the Everest Massif, High Himalaya of Nepal-South Tibet. J. Geol. Soc. London. 2003, 160, 345–366. [Google Scholar] [CrossRef] [Green Version]
  16. Simpson, R.L.; Parrish, R.R.; Searle, M.P.; Waters, D.J. Two Episodes of Monazite Crystallization during Metamorphism and Crustal Melting in the Everest Region of the Nepalese Himalaya. Geology 2000, 28, 403. [Google Scholar] [CrossRef]
  17. Searle, M.P.; Szulc, A.G. Channel Flow and Ductile Extrusion of the High Himalayan Slab-the Kangchenjunga-Darjeeling Profile, Sikkim Himalaya. J. Asian Earth Sci. 2005, 25, 173–185. [Google Scholar] [CrossRef]
  18. Liu, X.C.; Wu, F.Y.; Yu, L.J.; Liu, Z.C.; Ji, W.Q.; Wang, J.G. Emplacement Age of Leucogranite in the Kampa Dome, Southern Tibet. Tectonophysics 2016, 667, 163–175. [Google Scholar] [CrossRef]
  19. Liu, Z.C.; Wu, F.Y.; Ji, W.Q.; Wang, J.G.; Liu, C.Z. Petrogenesis of the Ramba Leucogranite in the Tethyan Himalaya and Constraints on the Channel Flow Model. Lithos 2014, 208, 118–136. [Google Scholar] [CrossRef]
  20. Wu, F.; Liu, X.-C.; Liu, Z.-C.; Wang, R.-C.; Xie, L.; Wang, J.; Ji, W.; 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]
  21. Visona, D.; Carosi, R.; Montomoli, C.; Tiepolo, M.; Peruzzo, L. Lithos Miocene Andalusite Leucogranite in Central-East Himalaya (Everest–Masang Kang Area): Low-Pressure Melting during Heating. Lithos 2012, 144–145, 194–208. [Google Scholar] [CrossRef]
  22. Le Fort, P. Metamorphism and Magmatism during the Himalayan Collision. Geol. Soc. Lond. Spec. Publ. 1986, 19, 159–172. [Google Scholar] [CrossRef]
  23. Smith, C.P.; Gübelin, E.J.; Bassett, A.M.; Manandhar, M.N. Rubies and Fancy-Color Sapphires from Nepal. Gems Gemol. 1997, 33, 24–41. [Google Scholar] [CrossRef]
  24. Dill, H.G. Pegmatites and Aplites: Their Genetic and Applied Ore Geology. Ore Geol. Rev. 2015, 69, 417–561. [Google Scholar] [CrossRef]
  25. ESCAP. Atlas of Mineral Resources of the ESCAP Region. Geology and Mineral Resources of Nepal. Explanatory Brochure. V. 9; United Nations: New York, NY, USA, 1993; Volume 9.
  26. DMG. Mineral Resources of Nepal, 2nd ed.; Department of Mines and Geology: Kathmandu, Nepal, 2017.
  27. Li, X.; Liu, Y.; Tu, X.; Hu, G.; Zeng, W. Precise Determination of Chemical Compositions in Silicate Rocks Using ICP-AES and ICP-MS: A Comparative Study of Sample Digestion Techniques of Alkali Fusion and Acid Dissolution. Geochimica 2002, 31, 289–294. [Google Scholar] [CrossRef]
  28. Monier, G.; Robert, J.L. Muscovite Solid Solutions in the System K2O-MgO-FeO-Al2O3-SiO2-H2O: An Experimental Study at 2 Kbar PH2O and Comparison with Natural Li-Free White Micas. Mineral. Mag. 1986, 50, 257–266. [Google Scholar] [CrossRef]
  29. Tindle, A.G.; Webb, P.C. Estimation of Lithium Contents in Trioctahedral Micas Using Microprobe Data: Application to Micas from Granitic Rocks. Eur. J. Mineral. 1990, 2, 595–610. [Google Scholar] [CrossRef]
  30. Xie, L.; Zhang, Y.; Zhang, H.; Sun, J.; Wu, F. In Situ Simultaneous Determination of Trace Elements, U-Pb and Lu-Hf Isotopes in Zircon and Baddeleyite. Chin. Sci. Bull. 2008, 53, 1565–1573. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, S.; Karius, V.; Schmidt, B.C.; Simon, K.; Wörner, G. Comparison of Ultrafine Powder Pellet and Flux-Free Fusion Glass for Bulk Analysis of Granitoids by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry. Geostand. Geoanalytical Res. 2018, 42, 575–591. [Google Scholar] [CrossRef]
  32. Wu, S.; Wörner, G.; Jochum, K.P.; Stoll, B.; Simon, K.; Kronz, A. The Preparation and Preliminary Characterisation of Three Synthetic Andesite Reference Glass Materials (ARM-1, ARM-2, ARM-3) for In Situ Microanalysis. Geostand. Geoanalytical Res. 2019, 43, 567–584. [Google Scholar] [CrossRef] [Green Version]
  33. Griffin, W.L. GLITTER: Data Reduction Software for Laser Ablation ICP-MS. Laser Ablation ICP-MS Earth Sci. Curr. Pract. Outst. Issues 2008, 308–311. [Google Scholar]
  34. Rushmer, T.; Knesel, K. Defining Geochemical Signatures and Timescales of Melting Processes in the Crust: An Experimental Tale of Melt Segregation, Migration and Emplacement. In Timescales of Magmatic Processes; John Wiley & Sons, Ltd.: Chichester, UK, 2010; pp. 181–211. [Google Scholar] [CrossRef]
  35. Middlemost, E.A.K. Naming Materials in the Magma/Igneous Rock System. Earth-Sci. Rev. 1994, 37, 215–224. [Google Scholar] [CrossRef]
  36. Maniar, P.D.; Piccoli, P.M. Tectonic Discrimination of Granitoids. GSA Bull. 1989, 101, 635–643. [Google Scholar] [CrossRef]
  37. Boynton, W.V. Cosmochemistry of the Rare Earth Elements: Meteorite Studies. Dev. Geochem. 1984, 2, 63–114. [Google Scholar] [CrossRef]
  38. Sun, S.-S.; McDonough, W.F. Chemical and Isotopic Systematics of Oceanic Basalts: Implications for Mantle Composition and Processes. Geol. Soc. Lond. Spec. Publ. 1989, 42, 313–345. [Google Scholar] [CrossRef]
  39. Trueman, D.L.; Černý, P. Exploration for Rare-Element Granitic Pegmatites. In Proceedings of the Short Course in Granitic Pegmatites in Science and Industry; Mineralogical Association of Canada: Quebec, QC, Canada, 1982; pp. 463–493. [Google Scholar]
  40. Morteani, G.; Preinfalk, C.; Horn, A.H. Classification and Mineralization Potential of the Pegmatites of the Eastern Brazilian Pegmatite Province. Miner. Depos. 2000, 35, 638–655. [Google Scholar] [CrossRef]
  41. Tischendorf, G.; Gottesmann, B.; Förster, H.-J.; Trumbull, R.B. On Li-Bearing Micas: Estimating Li from Electron Microprobe Analyses and an Improved Diagram for Graphical Representation. Mineral. Mag. 1997, 61, 809–834. [Google Scholar] [CrossRef]
  42. Miller, C.F.; Stoddard, E.F.; Bradfish, L.J.; Dollase, W.A. Composition of Plutonic Muscovite: Genetic Implications. Can. Mineral. 1981, 19, 25–34. [Google Scholar]
  43. Cerny, P.; Burt, D.M. Paragenesis, Crystallochemical Characteristics, and Geochemical Evolution of the Micas in Granite Pegmatites. Rev. Mineral. Geochem. 1984, 13, 257–297. [Google Scholar]
  44. Cameron, E.N.; Jahns, R.H.; Mcnair, A.H.; Page, L.R. Internal Structure of Granitic Pegmatites, No. 2. Econ. Geol. Urbana. Ill. 1949, 2, 115. [Google Scholar]
  45. Ginsburg, A.I.; Timofeyev, I.N.; Feldman, L.G. Principles of Geology of the Granitic Pegmatites. Nedra. Moscow 1979, 296. [Google Scholar]
  46. Zagorsky, V.Y.; Makagon, V.M.; Shmakin, B.M. The Systematics of Granitic Pegmatites. Can. Mineral. 1999, 37, 800–802. [Google Scholar]
  47. Černý, P. Anatomy and Classification of Granitic Pegmatites. In Granitic Pegmatites in Science and Industry; Černý, P., Ed.; Mineralogical Association of Canada: Québec, QC, Canada, 1982; Volume 8, pp. 1–39. [Google Scholar]
  48. Selway, J.B.; Breaks, F.W.; Tindle, A.G. A Review of Rare-Element (Li-Cs-Ta) Pegmatite Exploration Techniques for the Superior Province, Canada, and Large Worldwide Tantalum Deposits. Explor. Min. Geol. 2005, 14, 1–30. [Google Scholar] [CrossRef]
  49. 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]
  50. London, D. Pegmatites; Martin, R.F., Ed.; Mineralogical Association of Canada: Quebec, QC, Canada, 2008; Volume 10. [Google Scholar]
  51. London, D. Ore-Forming Processes within Granitic Pegmatites. Ore Geol. Rev. 2018, 101, 349–383. [Google Scholar] [CrossRef]
  52. Černý, P.; London, D.; Novák, M. Granitic Pegmatites as Reflections of Their Sources. Elements 2012, 8, 289–294. [Google Scholar] [CrossRef]
  53. Duc-Tin, Q.; Keppler, H. Monazite and Xenotime Solubility in Granitic Melts and the Origin of the Lanthanide Tetrad Effect. Contrib. Mineral. Petrol. 2015, 169, 1–26. [Google Scholar] [CrossRef]
  54. Gehrels, G.E.; DeCelles, P.G.; Martin, A.; Ojha, T.P.; Pinhassi, G.; Upreti, B.N. Initiation of the Himalayan Orogen as an Early Paleozoic Thin-Skinned Thrust Belt. GSA Today 2003, 13, 4. [Google Scholar] [CrossRef]
  55. Cawood, P.A.; Johnson, M.R.W.; Nemchin, A.A. Early Palaeozoic Orogenesis along the Indian Margin of Gondwana: Tectonic Response to Gondwana Assembly. Earth Planet. Sci. Lett. 2007, 255, 70–84. [Google Scholar] [CrossRef]
  56. Hopkinson, T.N.; Harris, N.B.W.; Warren, C.J.; Spencer, C.J.; Roberts, N.M.W.; Horstwood, M.S.A.; Parrish, R.R.; EIMF. The Identification and Significance of Pure Sediment-Derived Granites. Earth Planet. Sci. Lett. 2017, 467, 57–63. [Google Scholar] [CrossRef] [Green Version]
  57. Wang, J.-G.; Liu, X.-C.; Liu, Z.-C.; Ji, W.-Q.; Wu, F.-Y.; Ding, L. Highly Fractionated Late Eocene (~35 Ma) Leucogranite in the Xiaru Dome, Tethyan Himalaya, South Tibet. Lithos 2015, 240–243, 337–354. [Google Scholar] [CrossRef]
  58. 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]
  59. Qin, K.; Zhao, J.; He, C.; Shi, R. Discovery of the Qongjiagang Giant Lithium Pegmatite Deposit in Himalaya, Tibet, China. Acta Petrol. Sin. 2021, 37, 3277–3286. [Google Scholar] [CrossRef]
  60. Zhou, Q.; Qin, K.; He, C.; Wu, H.; Liu, Y.; Niu, X.; Mo, L.; Liu, X.; Zhao, J. Li-Be-Nb-Ta Mineralogy of the Kuqu Leucogranite and Pegmatite in the Eastern Himalaya, Tibet, and Its Implication. Acta Petrol. Sin. 2021, 37, 3305–3324. [Google Scholar] [CrossRef]
  61. London, D.; Evensen, J.M. Beryllium in Silicic Magmas and the Origin of Beryl-Bearing Pegmatites. Rev. Mineral. Geochem. 2002, 50, 445–486. [Google Scholar] [CrossRef]
  62. Zhao, J.; He, C.; Qin, K.; Shi, R.; Liu, X.; Hu, F.; Yu, K.; Sun, Z. Geochronology, Source Features and the Characteristics of Fractional Crystallization in Pegmatite at the Qongjiagang Giant Pegmatite-Type Lithium Deposit, Himalaya, Tibet. Acta Petrol. Sin. 2021, 37, 3325–3347. [Google Scholar] [CrossRef]
  63. Müller, A.; Keyser, W.; Simmons, W.B.; Webber, K.; Wise, M.; Beurlen, H.; Garate-Olave, I.; Roda-Robles, E.; Galliski, M.Á. Quartz Chemistry of Granitic Pegmatites: Implications for Classification, Genesis and Exploration. Chem. Geol. 2021, 584, 120507. [Google Scholar] [CrossRef]
  64. Breiter, K.; Ďurišová, J.; Dosbaba, M. Chemical Signature of Quartz from S- and A-Type Rare-Metal Granites – A Summary. Ore Geol. Rev. 2020, 125, 103674. [Google Scholar] [CrossRef]
  65. Evensen, J.M.; London, D. Experimental Silicate Mineral/Melt Partition Coefficients for Beryllium and the Crustal Be Cycle from Migmatite to Pegmatite. Geochim. Cosmochim. Acta 2002, 66, 2239–2265. [Google Scholar] [CrossRef]
  66. 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]
  67. Hezel, D.C.; Angelika, K.; Marschall, H.R.; Ludwig, T.; Meyer, H.P. Major-Element and Li, Be Compositional Evolution of Tourmaline in an S-Type Granite-Pegmatite System and Its Country Rocks: An Example from Ikaria, Aegean Sea, Greece. Can. Mineral. 2011, 49, 321–340. [Google Scholar] [CrossRef]
  68. Simmons, W.B.; Pezzotta, F.; Shigley, J.E.; Beurlen, H. Granitic Pegmatites as Sources of Colored Gemstones. Elements 2012, 8, 281–287. [Google Scholar] [CrossRef]
  69. Jahns, R.H.; Burnham, C.W. Experimental Studies of Pegmatite Genesis; l, A Model for the Derivation and Crystallization of Granitic Pegmatites. Econ. Geol. 1969, 64, 843–864. [Google Scholar] [CrossRef]
  70. London, D.; Hervig, R.L.; Morgan, G.B. Melt-Vapor Solubilities and Elemental Partitioning in Peraluminous Granite-Pegmatite Systems: Experimental Results with Macusani Glass at 200 MPa. Contrib. Mineral. Petrol. 1988, 99, 360–373. [Google Scholar] [CrossRef]
  71. Cerny, P.; Meintzer, R.E.; Anderson, A.J. Extreme Fractionation in Rare-Element Granitic Pegmatites: Selected Examples of Data and Mechanisms. Can. Mineralogist. 1985, 23, 381–421. [Google Scholar]
  72. Zhou, Q.; Qin, K.; Tang, D.; Ding, J.; Guo, Z. Mineralogy and Significance of Micas and Feldspars from the Koktokay No. 3 Pegmatitic Rare-Element Deposit, Altai. Acta Petrol. Sin. 2013, 29, 3004–3022. [Google Scholar]
  73. Michallik, R.M.; Wagner, T.; Fusswinkel, T.; Heinonen, J.S.; Heikkilä, P. Chemical Evolution and Origin of the Luumäki Gem Beryl Pegmatite: Constraints from Mineral Trace Element Chemistry and Fractionation Modeling. Lithos 2017, 274–275, 147–168. [Google Scholar] [CrossRef]
  74. Küster, D.; Romer, R.L.; Tolessa, D.; Zerihun, D.; Bheemalingeswara, K.; Melcher, F.; Oberthür, T. The Kenticha Rare-Element Pegmatite, Ethiopia: Internal Differentiation, U–Pb Age and Ta Mineralization. Miner. Depos. 2009, 44, 723–750. [Google Scholar] [CrossRef]
  75. Li, P.; Li, J.; Chen, Z.; Liu, X.; Huang, Z.; Zhou, F. Compositional Evolution of the Muscovite of Renli Pegmatite-Type Rare-Metal Deposit, Northeast Hunan, China: Implications for Its Petrogenesis and Mineralization Potential. Ore Geol. Rev. 2021, 138, 104380. [Google Scholar] [CrossRef]
Minerals 12 00564 g002 550
Figure 2. Field photographs of pegmatites and granites. (a) Orthogneiss exposure at Lonak (Kanchenjunga area); (b) Field view of Kanchenjunga–Jannu leucogranite at Mt. Jannu base camp; (c) Tour-leucogranite crosscutting the orthogneiss; (d) Tour-leucogranite; (e) Barren pegmatite intruded on banded gneiss at Chyanthapu area; (f,g) Be-rich pegmatite showing angular and sharp contact with host-rock gneiss at Ikhabu area; (h) Beryl seen in the outcrop of Be-rich pegmatite at Ikhabu area. Mineral abbreviations: Brl, beryl; Qz, quartz; Tur, tourmaline. Mineral abbreviations are after [10].
Figure 2. Field photographs of pegmatites and granites. (a) Orthogneiss exposure at Lonak (Kanchenjunga area); (b) Field view of Kanchenjunga–Jannu leucogranite at Mt. Jannu base camp; (c) Tour-leucogranite crosscutting the orthogneiss; (d) Tour-leucogranite; (e) Barren pegmatite intruded on banded gneiss at Chyanthapu area; (f,g) Be-rich pegmatite showing angular and sharp contact with host-rock gneiss at Ikhabu area; (h) Beryl seen in the outcrop of Be-rich pegmatite at Ikhabu area. Mineral abbreviations: Brl, beryl; Qz, quartz; Tur, tourmaline. Mineral abbreviations are after [10].
Minerals 12 00564 g002
Minerals 12 00564 g003 550
Figure 3. Field photographs and sketch map of the Yamrang Pegmatite. (a) Location of Yamrang Pegmatite at Yamrang hill; (b) Exposure of Yamrang Pegmatite, zone 3; (c) Schematic sketch of internal zonation of the Yamrang Pegmatite. Mineral abbreviations: Brl, beryl; Mc, microcline; Ms, muscovite; Qz, quartz. Mineral abbreviations are after [10].
Figure 3. Field photographs and sketch map of the Yamrang Pegmatite. (a) Location of Yamrang Pegmatite at Yamrang hill; (b) Exposure of Yamrang Pegmatite, zone 3; (c) Schematic sketch of internal zonation of the Yamrang Pegmatite. Mineral abbreviations: Brl, beryl; Mc, microcline; Ms, muscovite; Qz, quartz. Mineral abbreviations are after [10].
Minerals 12 00564 g003
Minerals 12 00564 g004 550
Figure 4. Field photographs and photomicrographs of the Yamrang Pegmatite. (a) Yamrang Pegmatite and host-rock gneiss contact zone; (b) Biotite-rich schistose gneiss hosting Yamrang Pegmatite under plane-polarized light; (c) Sample from saccharoidal albite zone showing albite-quartz assemblage; (d) Blocky perthitic microcline zone with miarolitic cavities; (e) Beryl and aquamarine in association with K-feldspar; (f) Quartz-beryl association from core zone; (g) Pegmatite sample showing graphic texture; (h) Aquamarine, quartz, and tourmaline gems from Yamrang Pegmatite; (i) Smoky quartz, muscovite books and albite assemblage from miarolitic cavities. Mineral abbreviations: Ab, albite; Aqm, aquamarine; Brl, beryl; Bt, biotite; Grt, garnet; Kfs, K-feldspar; Ms, muscovite; Pl, plagioclase; Qz, quartz; Tur, tourmaline. Mineral abbreviations are after [10].
Figure 4. Field photographs and photomicrographs of the Yamrang Pegmatite. (a) Yamrang Pegmatite and host-rock gneiss contact zone; (b) Biotite-rich schistose gneiss hosting Yamrang Pegmatite under plane-polarized light; (c) Sample from saccharoidal albite zone showing albite-quartz assemblage; (d) Blocky perthitic microcline zone with miarolitic cavities; (e) Beryl and aquamarine in association with K-feldspar; (f) Quartz-beryl association from core zone; (g) Pegmatite sample showing graphic texture; (h) Aquamarine, quartz, and tourmaline gems from Yamrang Pegmatite; (i) Smoky quartz, muscovite books and albite assemblage from miarolitic cavities. Mineral abbreviations: Ab, albite; Aqm, aquamarine; Brl, beryl; Bt, biotite; Grt, garnet; Kfs, K-feldspar; Ms, muscovite; Pl, plagioclase; Qz, quartz; Tur, tourmaline. Mineral abbreviations are after [10].
Minerals 12 00564 g004
Minerals 12 00564 g008 550
Figure 8. Photomicrographs and BSE images of alkali feldspar, muscovite, and quartz from the Yamrang Pegmatite. (a) Assemblage of albite, primary muscovite, and secondary mica overgrowing albite from zone 1 under cross-polarized light; (b) BSE image of homogeneous muscovite and quartz from zone 2; (c) Perthitic microcline and tourmaline from zone 2 under plane-polarized light; (d) Assemblage of perthitic microcline, anhedral quartz, and tourmaline from zone 3 under cross-polarized light; (e) BSE image of perthitic microcline from zone 4; (f) BSE image of homogeneous quartz from the core; (g) BSE image of muscovite sampled near core zone 4; (h) BSE image of miarolitic homogeneous muscovite and quartz; (i) Platy albite and muscovite books assemblage from miaroles under cross-polarized light. Mineral abbreviations: Ab, albite; Mc, microcline; Ms, muscovite; Qz, quartz; Tur, tourmaline. Mineral abbreviations are after [10].
Figure 8. Photomicrographs and BSE images of alkali feldspar, muscovite, and quartz from the Yamrang Pegmatite. (a) Assemblage of albite, primary muscovite, and secondary mica overgrowing albite from zone 1 under cross-polarized light; (b) BSE image of homogeneous muscovite and quartz from zone 2; (c) Perthitic microcline and tourmaline from zone 2 under plane-polarized light; (d) Assemblage of perthitic microcline, anhedral quartz, and tourmaline from zone 3 under cross-polarized light; (e) BSE image of perthitic microcline from zone 4; (f) BSE image of homogeneous quartz from the core; (g) BSE image of muscovite sampled near core zone 4; (h) BSE image of miarolitic homogeneous muscovite and quartz; (i) Platy albite and muscovite books assemblage from miaroles under cross-polarized light. Mineral abbreviations: Ab, albite; Mc, microcline; Ms, muscovite; Qz, quartz; Tur, tourmaline. Mineral abbreviations are after [10].
Minerals 12 00564 g008
Minerals 12 00564 g009 550
Figure 9. The ternary diagram of Ab-An-Or composition of alkali feldspar from five zones of the Yamrang Pegmatite.
Figure 9. The ternary diagram of Ab-An-Or composition of alkali feldspar from five zones of the Yamrang Pegmatite.
Minerals 12 00564 g009
Minerals 12 00564 g010 550
Figure 10. (a) Plot of K/Rb vs. Cs (ppm); Box and whisker plots of (b) Be; (c) K/Cs; (d) K/Rb; (e) Nb/Ta; (f) Rb/Sr of alkali feldspar from five zones of the Yamrang Pegmatite. Connecting lines join the median value. Ratios are calculated as ppm/ppm. Discrimination lines of barren, Be-, and Li, Be- pegmatites are from [39] as presented in [40].
Figure 10. (a) Plot of K/Rb vs. Cs (ppm); Box and whisker plots of (b) Be; (c) K/Cs; (d) K/Rb; (e) Nb/Ta; (f) Rb/Sr of alkali feldspar from five zones of the Yamrang Pegmatite. Connecting lines join the median value. Ratios are calculated as ppm/ppm. Discrimination lines of barren, Be-, and Li, Be- pegmatites are from [39] as presented in [40].
Minerals 12 00564 g010
Minerals 12 00564 g011 550
Figure 11. (a) Fe+Mg+T-AlVI vs. Mg-Li diagram of mica showing the position of muscovite of Yamrang Pegmatite [41]; (b) Total Al (apfu) vs. Fe/(Fe+Mg) diagram of muscovite from Yamrang Pegmatite. The discrimination line of primary and secondary muscovite is from [42].
Figure 11. (a) Fe+Mg+T-AlVI vs. Mg-Li diagram of mica showing the position of muscovite of Yamrang Pegmatite [41]; (b) Total Al (apfu) vs. Fe/(Fe+Mg) diagram of muscovite from Yamrang Pegmatite. The discrimination line of primary and secondary muscovite is from [42].
Minerals 12 00564 g011
Minerals 12 00564 g012 550
Figure 12. Plots of (a) K/Rb vs. Cs (ppm); (b) K/Rb vs. Be (ppm); Box and whisker plots of (c) Be; (d) K/Rb; (e) K/Cs; (f) Nb/Ta of muscovite from five zones of Yamrang Pegmatite. Connecting lines join the median value. Ratios are calculated from ppm/ppm value. Discrimination lines of Mus-class pegmatites and REL-class pegmatites are from [43], as presented in [40].
Figure 12. Plots of (a) K/Rb vs. Cs (ppm); (b) K/Rb vs. Be (ppm); Box and whisker plots of (c) Be; (d) K/Rb; (e) K/Cs; (f) Nb/Ta of muscovite from five zones of Yamrang Pegmatite. Connecting lines join the median value. Ratios are calculated from ppm/ppm value. Discrimination lines of Mus-class pegmatites and REL-class pegmatites are from [43], as presented in [40].
Minerals 12 00564 g012
Minerals 12 00564 g013 550
Figure 13. Schematic P-T fields of regional host rocks of five classes of granitic pegmatites: abyssal (AB), muscovite (MS), muscovite - rare-element (MSREL), rare-element (REL), and miarolitic (MI) (redrawn from [3]). Yellow star denotes the Yamrang Pegmatite in the REL pegmatite field.
Figure 13. Schematic P-T fields of regional host rocks of five classes of granitic pegmatites: abyssal (AB), muscovite (MS), muscovite - rare-element (MSREL), rare-element (REL), and miarolitic (MI) (redrawn from [3]). Yellow star denotes the Yamrang Pegmatite in the REL pegmatite field.
Minerals 12 00564 g013
Minerals 12 00564 g014 550
Figure 14. Box and whisker plot of Be-concentration in minerals and whole rock samples from five zones of Yamrang Pegmatite. Zone 1–4 represents the magmatic stage, and zone 5 represents the hydrothermal stage. Beryl saturation occurred at the late magmatic to the hydrothermal stage. Be content of beryl and tourmaline are taken from unpublished data of this author.
Figure 14. Box and whisker plot of Be-concentration in minerals and whole rock samples from five zones of Yamrang Pegmatite. Zone 1–4 represents the magmatic stage, and zone 5 represents the hydrothermal stage. Beryl saturation occurred at the late magmatic to the hydrothermal stage. Be content of beryl and tourmaline are taken from unpublished data of this author.
Minerals 12 00564 g014
Table
Table 1. Description of internal zonation of Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Table 1. Description of internal zonation of Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Zone NumberZone Name (Sequential)Zone Name (Mineralogical)ThicknessMain Mineral Assemblage (Relative Abundances wt.%)Accessory and (Rare) MineralsDescriptive FeaturesTexture
Zone 1Marginal ZoneSaccharoidal albite2 mAb (56), Mc (8), Qz (30), Ms (4)Drv-Srl, (Fl, Zrn, Mnz)saccharoidal albite, inwardly oriented tourmaline layersFine to medium-grained
Zone 2Outer Intermediate ZoneBlocky perthitic microcline4–7 mAb (10), Mc (70), Qz (15), Ms (3)Sps, Srl-Drv, (Zrn, Mnz)Massive microcline, graphic pegmatite, platy albite, quartz lenses, muscovite clustersmonomineralic, massive
Zone 3Inner Intermediate ZoneMuscovite-microcline-quartz5–6 mAb (15), Mc (40), Qz (25), Ms (15)Brl, Sps, Rt, Srl, (Col-Tan, Zrn, Mnz, Xtm)coarse graphic pegmatite, platy albite, book muscovite veins, quartz-beryl pods, quartz-muscovite nestscoarse-grained to massive
Zone 4Core ZoneBeryl-quartz1–1.5 mMc (5), Qz (60), Ms (1), Brl (33)Srl, (Zrn, Mnz)quartz-beryl 90%, microcline patches at core marginmedium to coarse-grained
Zone 5Miarolitic ZoneAlbite-quartzfew to tens of cmAb (40), Mc (15), Qz (25), Ms (15)Brl, Aqm, Col-Tan, Rt, Srl, Hem, (Zrn, Mnz, Xtm)clay filled pockets, platy albite, dipyramidal smoky quartz, aquamarine, muscovite books, overprinted on zone 2 and 3,fine to coarse-grained
Note: Ab, albite; Aqm, aquamarine; Mc, microcline; Qz, quartz; Ms, muscovite; Brl; beryl; Col-Tan, columbite-tantalite; Drv-Srl, dravite-schorl; Fl, fluorite; Zrn, zircon; Mnz, monazite; Sps, spessartine; Rt, rutile; Xtm, xenotime; Hem, hematite. Mineral abbreviations are after [10].
Table
Table 2. Whole-rock major and trace element composition of selected pegmatite, gneiss and granite samples.
Table 2. Whole-rock major and trace element composition of selected pegmatite, gneiss and granite samples.
LithologyBeryl PegmatiteParagneissBarren PegmatiteLeucograniteOrthogneiss
SampleZone 1Zone 2Zone 3Zone 5WPSPYRP-14IW-7WPT-88WPT-112KGP-20WPT-205KGP-6
wt.%
SiO273.1579.1670.2275.1073.2472.3671.8381.7378.2873.3973.1175.7676.49
Al2O316.2811.7616.1415.0614.7516.6312.117.9213.0515.1514.5214.1013.62
TiO20.020.050.010.010.030.050.500.220.060.020.140.010.11
TFe2O30.390.580.190.270.150.162.212.320.670.471.130.450.47
MgO0.070.050.000.010.060.102.743.150.220.130.230.110.22
MnO0.030.020.090.260.000.000.030.000.020.050.030.030.01
CaO0.510.700.130.470.362.211.580.460.550.571.100.300.56
Na2O7.433.982.738.312.636.521.712.460.923.393.424.226.32
K2O1.593.1810.110.287.460.542.691.124.356.244.744.470.86
P2O50.160.120.120.130.140.090.050.060.370.170.120.240.15
LOI0.360.390.250.080.510.474.420.511.470.290.840.270.27
TOTAL100.0100.0100.0100.099.399.199.9100.0100.099.999.4100.099.1
A/CNK1.711.501.251.661.411.792.031.962.241.51.571.571.76
μg/g
Li16.9632.3918.867.06bdl2.2014.2911.3312.4368.7667.5127.0512.07
Be5318.5617.371154.0111.262.420.452.851.7110.821.410.57
B0.000.000.000.000.000.000.000.0020.4676.070.001870.00
Sc1.640.660.310.310.090.6911.854.5718.311.851.472.954.05
Ti0.000.000.000.0011356.320.000.000.000.007210.00437
V0.250.270.110.264.257.0929524.530.801.679.830.445.39
Cr3.311.690.932.190.000.0058.4325.382.834.230.005.280.00
Mn0.000.000.000.0021.1929.870.000.000.000.002000.0056.14
Co0.330.700.080.090.440.182.193.710.500.480.920.180.30
Ni1.560.750.630.670.750.827.729.850.470.730.780.470.70
Cu2.327.781.132.161.121.084.771.362.272.370.432.800.15
Zn12.2521.782.753.540.075.3712314.2310.404.6344.007.867.65
Ga12.467.529.3214.153.9811.9517.2610.1020.4011.3812.9010.5113.16
Rb160161108317.20289.109.5813341.0118128544014131.46
Sr2.0221.141.950.6136.4988.4376.7422.5835.1615465.2119.9452.06
Zr1703.514.8610.7120.8511.8056.4657.1216.5916.9569.8916.3754.16
Nb7.908.858.2059.613.381.5312.487.7620.390.8920.511.2022.88
Mo0.050.010.010.020.050.0211.730.070.020.060.160.030.19
Sn0.000.000.000.006.061.640.000.000.000.0021.230.007.28
Cd0.240.120.090.340.000.000.100.010.020.090.000.090.00
In0.010.010.010.010.000.000.070.020.150.040.000.030.00
Cs53.6323.1093.326.2257.450.606.721.846.735.6633.950.181.00
Ba2.9332.975.513.13106928.9643893.56165.47711.3119347.0624.26
Hf8.630.150.280.711.010.501.611.740.570.443.400.443.00
Ta2.832.111.1917.122.060.660.860.691.670.137.180.1011.52
W0.861.330.691.348.743.381.460.3722.660.236.170.113.26
Tl1.100.835.250.201.520.061.400.150.961.742.630.870.19
Pb29.0066.3277.2523.0568.4993.7827.653.0717.6439.6576.5225.875.18
Bi2.5211.994.180.760.000.000.440.010.100.060.000.060.00
Th3.111.250.802.221.010.6212.4815.080.470.6414.880.158.20
U25.393.0914.765.879.3928.667.053.450.930.5850.100.8611.85
μg/g
Y37.275.060.454.144.597.4917.2810.7114.382.4114.881.9812.03
La6.701.460.331.852.244.1338.2325.622.423.3417.830.933.26
Ce12.252.380.513.604.476.4465.3246.628.005.5035.601.766.94
Pr1.400.250.050.370.470.658.956.220.960.434.180.160.85
Nd4.500.810.161.041.612.1529.8420.044.191.3715.220.523.09
Sm1.610.270.040.440.400.544.943.602.190.273.790.161.11
Eu0.040.260.020.010.291.130.890.680.190.780.550.040.14
Gd1.790.300.030.360.510.714.192.972.080.244.370.151.70
Tb0.630.090.010.120.110.200.610.410.570.040.710.040.43
Dy5.700.770.060.750.741.233.051.983.340.293.530.302.44
Ho1.190.150.010.090.170.280.600.350.530.060.630.050.47
Er4.240.490.030.200.590.941.640.961.370.181.640.171.36
Tm1.110.120.030.060.100.150.250.120.250.040.220.040.22
Yb9.180.720.040.260.811.061.550.631.570.271.440.311.63
Lu1.610.120.010.040.100.130.220.090.200.040.160.040.21
∑REE89.213.21.7613.317.227.217812142.315.31056.6535.9
Note: bdl indicates signal below the detection limit; Zone 1, 2, 3, 4 are from the Yamrang Pegmatite; WP, Wairung Pegmatite; SP, Sang-sabu Pegmatite; Major mineral associations are: Zone 1, albite, quartz, tourmaline; Zone 2, blocky perthitic microcline; Zone 3, musco-vite, microcline, quartz; Zone 5, miarolitic albite, quartz, beryl; WP, graphic quartz & K-feldspar, muscovite; SP, graphic quartz & feld-spar, tourmaline.
Table
Table 3. The representative EPMA major element composition of alkali feldspar from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Table 3. The representative EPMA major element composition of alkali feldspar from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Zone 1 (Marginal Zone)Zone 2 (Outer IZ)Zone 3 (Inner IZ)Zone 4 (Core)Zone 5 (Miarolitic Zone)
AlbiteMicroclineMicroclineMicroclineAlbite
SampleYRP-1A-F1YRP-1A-F2YRP-1A-F3YRP-12A-F1YRP-12A-F5YRP-6A-F1YRP-4-F2YRP-9A-F1YRP-8B-F2YRP-8B-F5YRP-7B-F2YRP-7B-F3
SiO2 (wt%)65.1665.0464.8663.8264.4464.6065.3664.4764.3064.3267.8668.22
TiO20.000.000.000.000.000.030.000.040.000.000.000.00
Al2O320.9321.0220.9117.7218.2317.9918.1917.7118.0117.9119.5618.86
FeO0.010.000.060.000.000.020.010.000.010.020.010.00
MnO0.000.000.000.000.020.030.000.000.020.030.000.00
MgO0.000.010.000.000.000.000.000.000.000.000.030.02
CaO2.562.502.490.010.000.040.040.000.050.010.780.36
Na2O10.6510.5910.661.471.001.061.161.321.581.3011.6712.03
K2O0.200.130.3013.6014.5714.5414.1214.0413.8713.940.170.09
Total99.5299.2899.2896.6198.2698.3098.8997.5897.8497.54100.0799.59
Oxygen atoms = 8
Si (apfu)2.8702.8732.8633.0423.0273.0333.0513.0453.0233.0402.9602.984
Ti0.0000.0000.0000.0000.0000.0010.0000.0010.0000.0000.0000.000
Al1.0871.0941.0880.9951.0090.9951.0010.9860.9980.9981.0050.972
Fe0.0000.0000.0020.0000.0000.0010.0010.0000.0000.0010.0000.000
Mn0.0000.0000.0000.0000.0010.0010.0000.0000.0010.0010.0000.000
Mg0.0000.0010.0000.0000.0000.0000.0000.0000.0000.0000.0020.001
Ca0.1210.1180.1180.0010.0000.0020.0020.0000.0020.0010.0360.017
Na0.9100.9070.9120.1360.0910.0970.1050.1210.1440.1190.9871.020
K0.0110.0070.0170.8270.8720.8700.8410.8460.8310.8400.0090.005
Total5.0005.0005.0005.0005.0005.0005.0005.0005.0005.0005.0005.000
End members
Ab87.2987.8587.1214.089.489.9711.0612.5114.7512.3895.5997.91
An11.6111.4411.260.060.010.190.200.000.250.053.511.63
Or1.100.701.6285.8590.5189.8488.7487.4985.0087.570.900.46
Note: IZ, Intermediate zone; structural formulae were calculated based on 8 No. of oxygen.
Table
Table 4. The average LA-ICP-MS trace element composition of alkali feldspar from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Table 4. The average LA-ICP-MS trace element composition of alkali feldspar from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Zone 1 (Marginal Zone)Zone 2 (Outer IZ)Zone 3 (Inner IZ)Zone 4 (Core)Zone 5 (Miarolitic Zone)
PhaseAlbiteMicroclineMicroclineMicroclineAlbite
n = 16n = 5n = 5n = 3n = 5
ElementAvgMinMaxσAvgMinMaxσAvgMinMaxσAvgMinMaxσAvgMinMaxσ
Li (ppm)2.311.274.270.816.303.6611.12.8824.23.4748.719.848.619.969.025.6-bdlbdl-
Be5.611.2710.43.425.294.367.151.195.513.488.522.324.123.944.402.2813.311.315.71.71
B22.511.443.88.726.721.137.57.538.127.349.015.4-bdlbdl--bdlbdl-
P84.541.521552.337630740540.247930859910869567170926384179490640.5
K/1032.330.653.911.181481431544.1913290.816530.713913414365.56.191.1915.45.94
Ti98.25.515640.96.85.109.71.9920.23.647.724.05.915.915.9149.84.463.815.100.91
Ga8.241.2118.75.8112.611.614.00.9413.010.314.91.718.614.021.16.9328.827.630.21.13
Rb5.870.9411.33.1193788497733.81272816166932026972417292598072.82.0817778.9
Sr1196.3323783.722.717.034.36.717.01.069.229.25.011.077.3564.60.940.521.740.49
Nb0.550.400.790.140.250.220.280.04-bdlbdl--bdlbdl--bdlbdl-
Sn2.950.766.921.794.303.464.890.754.564.144.990.377.765.579.062.182.162.162.16-
Sb1.060.352.700.73-bdlbdl--bdlbdl-0.92bdlbdl-1.631.471.780.22
Cs0.480.450.540.0419.112.740.612.195.413.532313110682.912210422.75.9541.017.6
Ba2170.924621517143.49519.041.11.7214559.629.929.830.011221.1bdlbdl-
La5.042.808.001.900.140.120.170.020.680.112.110.96-bdlbdl-0.710.321.660.55
Ce7.933.7613.03.240.130.110.150.030.880.063.811.640.070.070.073.230.850.252.210.78
Ta0.080.050.190.05-bdlbdl--bdlbdl--bdlbdl--bdlbdl-
Pb23.51.0215344.718616723628.415811819435.813311316764.158.738.998.324.1
K/Rb54034512652271581531623.9810587.712112.652.145.959.042.633886.71017454
K/Cs40671439723825549426368312,0523315409138767312824136111011721230945537556699.5
Nb/Ta7.162.7010.02.23----------------
Rb/Sr0.060.010.230.0543.527.652.09.3942511.811964521024398225962765.04.0016971.9
Note: IZ, Intermediate zone; σ, one standard deviation; bdl, below detection limit for the utilized method; -, value could not be determined. The data of albite phase of perthitic microcline from zone 2–4 are excluded in this table but are provided in Supplementary Table S3.
Table
Table 5. The representative EPMA major element composition of muscovite from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Table 5. The representative EPMA major element composition of muscovite from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
SampleZone 1 (Marginal Zone)Zone 2 (Outer IZ)Zone 3 (Inner IZ)Zone 4 (Core)Zone 5 (Miarolitic Zone)
YRP-13B-M1YRP-13B-M2YRP-2B-M1YRP-2B-M2YRP-6A-M1YRP-6A-M2YRP-10A-M1YRP-10A-M2YRP-5B-M1YRP-7F-M2
SiO247.0847.9546.9847.1548.2547.8746.1947.0047.2947.77
TiO20.130.220.200.240.110.120.160.130.130.51
Al2O332.1932.1033.1232.8834.8834.9032.9433.3533.5532.69
FeO2.412.292.362.302.011.941.911.931.902.15
MnO0.120.070.050.090.050.090.050.050.040.04
MgO0.871.171.121.100.470.420.940.830.850.94
CaO0.010.010.000.000.010.000.040.000.040.04
Na2O0.330.420.740.640.870.820.840.710.570.56
K2O9.539.399.089.139.099.139.269.409.309.34
F0.250.510.180.210.230.300.410.230.300.19
Cl0.000.010.010.000.010.000.000.000.000.00
Cr2O30.020.000.000.010.010.010.000.000.000.02
NiO0.020.000.010.010.000.050.000.000.010.03
Li2O*0.000.080.000.000.000.020.050.000.020.00
H2O*4.274.214.364.344.454.404.194.334.314.38
Subtotal97.2498.4398.2298.10100.45100.0896.9897.9698.3198.66
O=F,Cl0.110.220.080.090.100.120.170.100.130.08
Total97.1398.2198.1598.02100.3599.9596.8197.8698.1898.58
Structural formulae based on O, OH, F = 24
Si6.436.466.346.376.346.326.326.356.366.41
Al iv1.571.541.661.631.661.681.681.651.641.59
Y
Al vi3.613.563.613.603.743.753.633.673.683.58
Ti0.010.020.020.020.010.010.020.010.010.05
Cr0.000.000.000.000.000.000.000.000.000.00
Fe0.280.260.270.260.220.210.220.220.210.24
Mn0.010.010.010.010.010.010.010.010.000.00
Mg0.180.230.230.220.090.080.190.170.170.19
Ni0.000.000.000.000.000.000.000.000.000.00
Li*0.000.040.000.000.000.010.030.000.010.00
X
Ca0.000.000.000.000.000.000.010.000.010.01
Na0.090.110.190.170.220.210.220.190.150.15
K1.661.611.561.571.521.541.621.621.601.60
OH*3.893.783.923.913.903.883.823.903.873.92
F0.110.220.080.090.100.120.180.100.130.08
TOTAL17.8417.8517.8817.8617.8217.8317.9417.8817.8417.82
Y total4.094.134.134.124.074.084.094.074.094.07
X total1.751.731.761.741.751.751.841.811.751.75
Al total5.185.105.275.235.405.435.315.315.325.17
Fe/Fe+Mg0.610.520.540.540.700.720.530.570.560.56
Fe+Mg+Aliv+Ti4.074.074.124.114.074.064.064.074.084.06
Mg-Li0.180.190.230.220.090.070.160.170.160.19
Fe+Mg+Ti-Alvi−3.14−3.04−3.09−3.10−3.42−3.44−3.21−3.27−3.28−3.10
Note: IZ, Intermediate zone; *, calculated from stoichiometry; Li2O* calculation after [28]; H2O* calculation after [29].
Table
Table 6. The average LA-ICP-MS trace element composition of muscovite from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
Table 6. The average LA-ICP-MS trace element composition of muscovite from the Yamrang Pegmatite, Ikhabu Pegmatite Field (IPF), Eastern Nepal.
ElementZone 1 (Marginal Zone)Zone 2 (Outer IZ)Zone 3 (Inner IZ)Zone 4 (Core)Zone 5 (Miarolitic Zone)
n = 16n = 5n = 5n = 3n = 5
AvgMinMaxσAvgMinMaxσAvgMinMaxσAvgMinMaxσAvgMinMaxσ
Li (ppm)52643259384.036032939220.128225632724.353551955816.556748069662.9
Be30.329.930.80.4923.218.328.02.9429.425.633.42.8429.821.046.910.628.824.332.02.71
B15410919944.716711419529.1200.5193.5205.44.5716213120028.514011916714.5
Na42373562492868341453473448832554595367554663.6543641146869981425932385081599
Mg45213939484750548834553521924822332159239375.7507945435738486518445125851498
P77.567.388.510.665.255.172.89.1263.252.681.811.082.657.510122.487.971.611314.6
K/1031101101100.291091081111.221031001052.04104991135.511271221365.02
Sc14.213.914.90.5518.413.520.32.315.615.016.20.4815.011.318.02.8810.52.6420.87.95
Ti1009749124825011971128128657.070667873216.1101690212401462574117644911485
Mn46041048944.027524829722.942140643910.542836052470.5674515986171
Fe15,33714,08515,997108615,39414,50016,07254415,22914,52115,75240214,69512,83316,168157521,38416,70329,5754730
Zn73.154.489.317.555.648.465.46.310495.01207.862.618.975.324.498.859.720341.3
Ga12811215926.590.587.492.51.81461411513.11171061258.2915513418118.2
Rb247515404057137792077413541962125205622024915251480156140.84398217268812288
Sr0.730.420.970.281.50.92.00.40.470.180.770.30-bdlbdl-0.440.350.560.07
Zr1.600.543.281.472.51.43.10.61.291.001.850.292.751.844.200.932.221.742.550.25
Nb29618636796.587.668.714224.42452362526.6415813618218.0450213701234
In1.150.451.950.750.50.30.80.20.730.500.960.15-bdlbdl-2.050.613.581.41
Sn24673.547320587.463.315431.41401331444.4868.458.784.512.8491113880385
Cs18847.639017937.113.812139.963.661.966.21.4832.630.734.41.4631054.0621264.8
Ba29.321.134.17.1626.318.849.610.5-bdlbdl-10.62.6817.77.853.321.146.772.14
Ta10343.617666.912.44.945.014.615.113.915.90.8014.59.822.45.1390.817.917472.1
W38.328.444.48.6846.642.649.82.742.138.144.62.2641.532.261.811.649.045.453.42.7
Pb9.76.911.92.613.712.215.20.97.526.977.950.3612.19.715.52.3416.312.820.52.4
K/Rb52.827.271.222.912280.714019.648.447.849.90.7168.063.272.93.5538.918.059.921.2
K/Cs11502822305104252709067900280116181573167536.5318429613423199120719923751044
Nb/Ta3.41.94.31.311.03.115.94.416.215.317.00.6911.68.1513.82.427.34.012.93.4
Rb/Sr3369225841839976243899371826434275111,8374781----10,070430919,6605700
Note: IZ, Intermediate zone; σ, one standard deviation; bdl, below detection limit; -, value could not be determined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bhandari, S.; Qin, K.; Zhou, Q.; Evans, N.J. Geological, Mineralogical and Geochemical Study of the Aquamarine-Bearing Yamrang Pegmatite, Eastern Nepal with Implications for Exploration Targeting. Minerals 2022, 12, 564. https://doi.org/10.3390/min12050564

AMA Style

Bhandari S, Qin K, Zhou Q, Evans NJ. Geological, Mineralogical and Geochemical Study of the Aquamarine-Bearing Yamrang Pegmatite, Eastern Nepal with Implications for Exploration Targeting. Minerals. 2022; 12(5):564. https://doi.org/10.3390/min12050564

Chicago/Turabian Style

Bhandari, Sushmita, Kezhang Qin, Qifeng Zhou, and Noreen J. Evans. 2022. "Geological, Mineralogical and Geochemical Study of the Aquamarine-Bearing Yamrang Pegmatite, Eastern Nepal with Implications for Exploration Targeting" Minerals 12, no. 5: 564. https://doi.org/10.3390/min12050564

APA Style

Bhandari, S., Qin, K., Zhou, Q., & Evans, N. J. (2022). Geological, Mineralogical and Geochemical Study of the Aquamarine-Bearing Yamrang Pegmatite, Eastern Nepal with Implications for Exploration Targeting. Minerals, 12(5), 564. https://doi.org/10.3390/min12050564

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