Peer Reviewed Gems & Gemology, Spring 2015, Vol. 51, No. 1

Blue Spinel from the Luc Yen District of Vietnam


Rough blue spinel crystals from Luc Yen
Figure 1. The Luc Yen district of Vietnam has become a major source of top-quality blue spinel, including these two rough crystals (45 and 70 ct) and the 5 ct faceted stone. Photo by J.B Senoble; © Senoble & Bryl.

ABSTRACT

The Luc Yen district of northern Vietnam is a very productive gem province and the leading source of vivid blue spinel. This study characterizes the origin and gemological properties of these spinels, especially the cause of their unusually bright color, which is directly related to their value. Chemical and spectroscopic analyses indicated that the blue color is due to cobalt (Co2+), with some iron contribution. Petrographic examination identified the context of the gem’s formation, which appears to be linked to intense metamorphism during successive orogenies. The carbonate platforms in the ancient Paleo-Tethys Ocean were sandwiched and highly deformed during this orogeny, leading to marble and spinel formation. The authors propose that the cobalt (and to a lesser extent the iron) necessary for the blue color were transported by fluids during metamorphism of the sedimentary sequence.

INTRODUCTION

Blue spinels are mined in Sri Lanka, Tanzania, Myanmar, Pakistan, and Vietnam (Shigley and Stockton, 1984; Delaunay, 2008; Pardieu and Hughes, 2008). Vietnam’s two major spinel deposits, Luc Yen and Quy Chau, were discovered at the end of the 1980s. The Luc Yen deposits have mostly been mined since the 1990s (Pardieu and Hughes, 2008; Senoble, 2010). This area, known for its gem-quality ruby, red spinel, and sapphire (Webster, 1994; Hauzenberger et al., 2003; Long et al., 2004; Senoble, 2010; Huong et al., 2012), has also been a notable producer of vivid blue spinels since the 2000s (figure 1). Bright, saturated blue gems are very popular, as evidenced by the classic appeal of sapphire and the more recent trend of Paraíba tourmaline and bright blue apatite in the same color range. Therefore, the bright blue color of some spinels has increased the popularity of this gem overall (Delaunay, 2008; Senoble, 2010). In this article, we investigate the gemological characteristics of Vietnamese blue spinels and the geology of the deposits to gain a better understanding of this gem source.

LOCATION AND ACCESS

The Luc Yen district is located in the Yen Bai province, in the north of Vietnam. Luc Yen’s capital of Yen The (22°6′38.84″ N, 104°45′57.80″ E) is a five- or six-hour drive from Hanoi on a 160 km expanse of good road. All of the district’s blue spinel mining sites lie within 20 km of Yen The. Several hours of walking or biking are needed to access these mines. The blue spinel deposits are Bai Gou, May Trung, Bai Son, Bah Linh Mot, Khe Khi, Kuoi Ngan, Khao Ka, Lung Thin, Lung Day, Khin Khang, and Chuong Tran (figure 2). Of these, only May Trung, Bai Son, and Bah Linh Mot are primary deposits; the others are secondary placer deposits. Bai Son, reported by Senoble (2010), was no longer being mined as a primary deposit during our visit in February and March 2012. May Trung is divided into two sites located about 150 meters from each other: a marble cliff that is mined for red and lavender spinels, and a second site that is mined for blue spinel from a vein in marble.

Map of Luc Yen district
Figure 2. Most Vietnamese blue spinel deposits are confined to a 30 km² area in the Luc Yen district of northern Vietnam. The circles represent placer deposits, while diamonds indicate primary deposits in marble.

GEOLOGY

The rich tectonic history of Southeast Asia is inherited from several deformation episodes related to the closure of the Paleo-Tethys Ocean and, later, to the Himalayan orogeny. The geology of northern Vietnam is dominated by metamorphic rocks inherited from these two major orogenic events. The first one, the Indosinian orogeny, led to the collision of the main shields (Yangtze and Indochina) during the Permo-Triassic at about 240–245 Ma (Kušnír, 2000; Lepvrier et al., 2008; Huong et al., 2012). In the later orogeny, the Himalayan collision during the Tertiary period, the terrains were strongly reworked. These terrains are primarily composed of metamorphic rocks, mainly medium-grade mica schist, marble, and granulitic gneisses (Kušnír, 2000; Leloup et al., 2001; Hauzenberger et al., 2003).

Northern Vietnam has been studied extensively to understand how a continental collision (in this case, between India and Eurasia) induced crustal wedges to extrude laterally into the surrounding plates. (Tapponnier et al., 1982, 1990; Leloup et al., 1995, 2001; Jolivet et al., 2001; Anckiewicz et al., 2007). During the Oligo-Miocene (from 35 to 17 Ma), the Indo-Eurasian collision induced strong rock deformation over all of Southeast Asia. The Indochinese block was extruded toward the southeast, and this induced the Red River Shear Zone, extending from the Tibetan plateaus to the China Sea for more than 1,000 km (Jolivet et al., 2001; Leloup et al., 2001; Hauzenberger et al., 2003; Anckiewicz et al., 2007). The Yen Bai province is formed by two different geological units separated by a fault that is part of the Red River shear zone. To the northeast lies the Lo Gam zone, and to the southwest the Day Nui Con Voi range (figure 3).

Geological map of Luc Yen district
Figure 3. This geological map of Luc Yen shows two different geological formations: the Day Nui Con Voi Range in the southwest and the Lo Gam zone in the northeast. The blue spinel deposits are located in the Lo Gam zone. Adapted from Garnier (2003) and Long et al. (2004).

All of Luc Yen’s gem deposits are located in the Lo Gam zone (again, see figure 3). The structure of this unit results from the deformation of the Himalayan orogenesis superimposed on the preexisting Indosinian structure (Garnier et al., 2002, 2005). The Lo Gam formation consists of a sedimentary series metamorphosed into marble, gneiss, calc-silicates, micaschist, and amphibolite. These metamorphic rocks are sometimes intruded by granitic and pegmatitic dykes (Leloup et al., 2001; Garnier et al., 2005, 2008). The marbles are mainly calcitic and interlayered with Al-, V-, and Cr-rich amphibolites.

Blue spinel is found in a layer of marble more than 500 meters thick. It occurs in discontinuous series of lenses, tens of millimeters thick and meter-sized in length, roughly following the regional foliation. These marble lenses are remarkably rich in forsterite (magnesian olivine). The gem is often associated with calcite, forsterite, pargasite (sodi-calcic amphibole), sulfides, and chlorites (magnesian chlorite and clinochlore). Remarkably, blue spinel in these primary deposits is not associated with ruby or red spinel.

MINING

The three primary deposits at May Trung (22°1′48.9″ N, 104°48′42.7″ E), Bai Son (2159′47.3″ N, 104°40′9.9″ E), and Bah Linh Mot (22°1′23.7″ N, 104°48′42.8″E) are located on a mountain range composed of marble, standing about 600 meters high. Each site is mined by a handful of locals, mainly farmers trying to earn extra income. The blue spinel is extracted from the marble using hand tools (figure 4, top) and a jackhammer.

Mining of blue spinel at Luc Yen
Figure 4. Gem-quality blue spinels are extracted from primary deposits by locals using hand tools, as shown in Bai Nua Doi (top). Secondary deposits are exploited in the valley with a sluice box to sort the minerals according to their density (bottom left). In the karstic environment, caves trap the gem-rich gravel (bottom right). These special secondary deposits are also processed with sluice boxes. Photos by Boris Chauviré.

Secondary deposits (figure 4, bottom left and bottom right) yield most of the blue spinel production. Some lie in the valley to the east of the spinel-rich mountain range. These include Kuoi Ngan (22°0′7.8″ N, 104°50′41.1″ E); Khao Ka (21°59′6.5″ N, 104°50′52.5″ E); Lung Thin (22°0′12.8″ N, 104°49′31.5″ E); Lung Day (21°59′51.3″ N, 104°49′23.3″ E); and Khin Khang (22°1′46.7″ N, 104°50′9.3″ E). Khe Khi (22°1′36.8″ N, 104°48′41.2″ E) and Bai Gou (22°4′43.2″ N, 104°47′5.5″ E) are located in the mountain in a small secondary basin. Miners use a water hose and a sluice to sort the gem-bearing gravels (figure 4, bottom left). Some secondary deposits are localized in karst caves inside marble (figure 4, bottom right). Heavy gravels are washed and sorted inside the cave and brought up to the surface, where they are sorted under daylight. In secondary deposits, blue spinel is found together with ruby, red spinel, sapphire, tourmaline, and occasionally gold.

PRODUCTION AND DISTRIBUTION

In secondary deposits, blue spinel is a by-product of ruby and red spinel mining. Even so, some large parcels contain more than a thousand carats of millimeter-sized, very saturated blue spinel (see Pardieu, 2012). Some dark grayish blue stones weighing approximately 5 ct have been faceted, but far fewer than ruby and red spinel.

Blue spinels from primary deposits are found in two different forms. Usually miners encounter them in “pockets” as centimeter-sized crystals, occasionally with a pleasing, well-defined octahedral shape with some portion containing gem material, or as aggregates of small octahedral crystals of varying quality. These are broken down or cobbed by miners to extract a small amount of gem spinel. These gem blue spinel can reach 5 ct, but they are often fractured. Spinels are also found as millimeter-sized octahedra or twins in the marble.

The blue spinel’s hue, tone, and saturation vary from one deposit to the next. In the secondary Bai Gou deposit, the crystals have a very dark blue color and often reach 10 ct. In Chuong Tran and Bai Son, the spinel has a bright blue color (Senoble, 2010; Overton and Shen, 2011) and can reach 5 ct. Crystals from May Trung and Khe Khi have a very saturated cobalt blue color but are quite small (rarely larger than 1 ct).

Millimeter-sized blue spinels from the primary deposits are often kept in marble so that the whole piece can be carved. This is done directly at the mining site, and the carvings are taken down to the valley to be sold as decorative pieces.

MATERIALS AND METHODS

Sample Collection. In early 2012, two field trips were organized to collect samples and map the main blue spinel deposits within the Luc Yen district. The first expedition enabled us to visit most of the corundum and spinel deposits and to understand the geology of this area. We also made contact with local merchants and miners for the second expedition just one month later. On the second field trip, we visited only blue spinel deposits and collected whole rock and gem-bearing samples. Most of the rock samples were collected at the mining site. Unfortunately, we did not observe gem samples at the mining sites. All gem blue spinels for spectroscopic and gemological measurements were procured from several local merchants in Yen The.
 

Materials. From the 55 carats of blue spinel we collected from local merchants, six representative rough crystals were selected and prepared in parallel-window plates for gemological and spectroscopic investigation. The least included and fractured samples were chosen for spectroscopic analysis. A polished window was prepared on each one to facilitate gemological and microscopic examination. One additional sample similar to SATBLU1 in color and saturation (labeled SATBLUchem) was prepared as a polished section for laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) chemical analysis. Four additional samples from a later field trip by author EF were added to complete this study. These were also purchased from local merchants. The spinels were divided into three different parcels according to their color category (detailed in “Gemological Characteristics” below). The sample names consisted of the color category (GREBLU, SKYBLU and SATBLU) followed by a number; see table 1. Moreover, 73 rock samples were collected in the field from 11 different mining sites. From these, 19 thin sections were prepared for petrographic examination.

Table 1

 

Methods. Gemological Properties. Specific gravity was measured hydrostatically with a Mettler Toledo JB703-c/FACT (with a precision of 0.001 ct). Internal features were observed with a standard gemological microscope. Refractive index was measured with a PФ-І refractometer with Rayner SVLS orange light. All spinel samples were observed under a 6 W A-Krüss Optronic 240 UV light, and we also tested their Chelsea filter reaction. Color was documented under normalized daylight (D65) and a normalized incandescent light (A).

Spectroscopic Measurements. UV-Vis-NIR absorption spectra of each sample (window plates) were taken with a Cary 5G Varian spectrophotometer in the 200–1500 nm range with a sampling interval of 1 nm and a spectral bandwidth of 1 nm maximum (sampling and spectral bandwidth were sometimes reduced to 0.25 nm to obtain better resolution). Raman spectra were collected on gem samples using both a Jobin-Yvon Labram with a 514 nm, 50 mW laser excitation, and a Jobin-Yvon Spex Horiba T64000 with a 647 nm, 50 mW laser excitation. The spectral range extended from 40 to 1500 cm–1 with a two-second exposure.

Chemical Composition. LA-ICP-MS chemical analysis was conducted at Blaise Pascal University (Clermont-Ferrand, France) using an Agilent 7500 spectrometer with a Resonetics M-50E laser (193 nm ablation wavelength, 5 Hz frequency with an energy between 10 and 12 J/cm2). For these analyses, four indentations (about 73 µm in diameter) were ablated on each sample, and 27Al was used as the internal standard. Data was processed with the GLITTER 4.4.2 software. To complement these analyses, we used a RIGAKU NEX CG energy-dispersive X-ray fluorescence (EDXRF) spectrometer operating at 25 kV and 0.10 mA. The detection limit for the major elements (Al, Mg) is about 0.1 wt.%, and below 0.01 wt.% for the minor elements. Each sample was measured for 90 seconds.

Petrographic Examination. Thin sections of rocks were observed with a standard Wild Makroscope M420 petrographic microscope, and a JEOL JSL-5800 LV scanning electron microscope (SEM) operating at 20 kV and 0.3 nA electron beam, with a 37° take-off angle of the detector. Mineral compositions of the samples and their inclusions were first determined by energy-dispersive spectroscopy (EDS) using an IMIX-PTS detector. This detector uses a high-resolution (115 eV) Ge crystal and an ultrathin polymer window, detecting elements ideally down to boron, if it is a major component of the material. The calibration standards used were either pure elements or simple compounds. The PGT software applies phi-rho-z data correction for the effect of X-ray absorption in the analyzed material, taking into account all the matrix effects. Oxygen was calculated from the spectrum, not based on stoichiometry.

GEMOLOGICAL CHARACTERISTICS

Visual Appearance. We separated the spinel samples into three categories according to their color descriptions:

  • SATBLU samples: medium to medium dark tone, strong to vivid saturation, and blue to violetish blue hue
  • SKYBLU samples: medium light to very light tone, strong to very vivid saturation, and blue hue
  • GREBLU samples: medium light to light tone, grayish to slightly grayish saturation, and blue to bluish violet hue

All of the rough samples were slightly fractured and contained very few inclusions. Color was homogeneous in each stone, and most showed a subtle color change from blue under daylight-equivalent normalized light (D65) to violetish blue under incandescent light (see table 1). The authors avoid the commonly used term “color shift” (Senoble, 2010), which Manson and Stockton (1984) defined in garnets while observing the combination of two color phenomena, nowadays identified separately: classical color change with lighting, and Usambara effect (change of color with thickness). We observed that the color change is more pronounced in stones with a more saturated color. While examining numerous parcels in Yen The, we observed that most of the grayish blue spinel—and some of the very saturated blue material, contrary to our other observations—did not show any color change.

Blue spinel from secondary deposits (except Khe Khi) is rounded and can reach several tens of carats. In Khe Khi, blue spinels are found as millimeter-sized octahedra.

Optical and Physical Properties. The samples’ refractive index ranged from 1.711 to 1.718, and their specific gravity was from 3.578 to 3.673. They were isotropic, with no anomalous double refringence, and inert under both short- and long-wave UV light. Under the Chelsea filter, all the samples appeared pink to red (see table 1 for details). We observed that the darker the spinel, the redder the Chelsea filter reaction.

Microscopic Characteristics. Conchoidal fractures and “fingerprint” healed fractures were often present in our samples (figure 5, left). Some showed elongated tubes, while others contained groups of parallel tubes. We observed birefringence in some of these tubes, which suggested that they consisted of an anisotropic solid phase (figure 5, middle and right). Black, opaque, irregular to hexagonal crystal inclusions less than 1 mm, reminiscent of graphite, were also found in some samples (figure 5, middle and right). GREBLU1 was the only sample that had yellowish fractures covered by red crystals (probably ferric oxide hematite).

Internal features of Vietnamese blue spinel
Figure 5. Blue spinels (here, GREBLU3) often show healed fractures (fingerprints, left), and some samples contain irregular opaque black crystals associated with elongated tubes (center and right). Photos by Boris Chauviré; field of view 1 mm (under daylight equivalent light on the left, plane-polarized light in the center, and cross-polarized light on the right).

PETROGRAPHY AND CHEMISTRY OF HOST ROCKS

Minerals. The marble that hosts blue spinel is mainly composed of calcite (sometimes magnesian) and dolomite. The major additional phases are olivine and pargasite (figure 6). Several accessory phases were identified using the petrographic microscope or EDS with SEM. These included titanite, rutile, zircon, graphite, apatite, several sulfide minerals (again, see figure 6), and phyllosilicates. The sulfides were mainly pyrrhotite (Fe1–xS; monoclinic) with pentlandite exsolution lamellae ((Fe,Ni)9S8; cubic) and violarite (FeNi2S4; cubic). Raman spectroscopy helped to distinguish between different phyllosilicates, mainly clinochlore and phlogopite. Humite was not observed in the marble, although this mineral is associated with red or purple spinel, as well as ruby (Hauzenberger et al., 2003; Garnier et al, 2008).

Thin sections of spinel-bearing rock
Figure 6. These views of thin sections from rocks bearing blue spinel (under plane-polarized light) show that blue spinel is always associated with olivine (forsterite) and pargasite in calcite matrix. Clinochlore surrounds all main minerals (left, field of view 1.5 mm). In the matrix, graphite and pyrrhothite are common accessory minerals (right, field of view 0.5 mm). Cc = calcite, Clh = clinochlore, Gph = graphite, Ol = olivine, Pg = pargasite, Pyr = pyrrhotite, Sp = spinel. Photomicrographs by Boris Chauviré.

Texture. The marble that hosts blue spinel has a granoblastic texture, with millimeter to centimeter grain size. SEM imaging with a backscattered electron detector showed exsolution features between calcite and dolomite, and intergrown apatite and calcite (figure 7a). Pentlandite lamellae in pyrrhotite present two different morphologies. The first consists of parallel flat lamellae less than 500 nm thick, crossing some pyrrhotite crystals from end to end. The second is lens-shaped, more than 1 µm thick and about 3 µm long, often associated with parallel flat lamellae (figure 7b).

Inclusions in blue spinel from Luc Yen
Figure 7. An inclusion of apatite in pargasite exhibits intergrowth with calcite (slightly magnesian; left, magnification 750×). In most cases, pyrrhotite inclusions have exsolutions of pentlandite, a sulfide with higher nickel content (center, magnified 1200×), which also contains cobalt. Using scanning electron microscopy with backscattered electron imaging, sensitive to the atomic number, a petrographic thin section of marble-bearing blue spinel shows that spinel and olivine are surrounded by clinochlore. The marble is composed of calcite and dolomite (right, magnified 65×). Ap = apatite, Cc = calcite, Clh = clinochlore, Dol = dolomite, Ol = olivine, Pen = pentlandite, Pg = pargasite, Pyr = pyrrhotite, Sp = spinel.

Paragenesis. Blue spinel is observed only in olivine-rich lenses, associated with dolomite and calcite (figure 8). No blue spinel is observed in the marble when olivine is absent. The spinel-rich lenses are elongated nearly parallel to the regional foliation: roughly 45° toward the northeast. Spinel and pargasite show inclusions of apatite and sulfides similar in shape and composition for both host minerals. This suggests that apatite and sulfides preexisted spinel and pargasite. In some titanite crystals, SEM imaging revealed inclusions of zircon and pargasite. Therefore, titanite probably represents a later stage of mineralization. Clinochlore crystals surround all the other minerals (figure 7c), meaning it probably crystallized later during a hydration phase, and possibly during exhumation.

Mable-hosted lenses in which blue spinel deposits appear]
Figure 8. In Vietnam, primary blue spinel deposits appear as approximately lens-shaped bodies rich in olivine. These lenses are hosted in marble, and pargasite is found throughout the surrounding marble. Photo and drawing by Boris Chauviré.

CHEMICAL COMPOSITION

Spinel. The composition of the three types of spinel crystals was measured in thin sections using EDS, and all rough samples were analyzed by EDXRF. These analyses identified them as spinel sensu stricto (MgAl2O4). Table 2 presents LA-ICP-MS chemical analyses on representative samples of the three color categories. The main impurities detected were Li, Fe, and Zn. Significant traces of Be, Ti, V, Cr, Mn, Ga, Ni, and Co were also detected. All analyzed spinels had concentrations of Ga, Zn, and Li, consistent with those observed only in natural blue spinel (Muhlmeister et al., 1993; Krzemnicki, 2008; Saeseaw et al., 2009). All the samples presented nearly uniform concentration in iron and in copper: around 10,000 ppma (equal to 1 atomic percent) and 4 ppma, respectively. The other elements showed strong variation among samples. Sample SATBLUchem (saturated blue) was enriched in Ti, V, Cr, Mn, Co, and Ni compared to the other samples. Samples SKYBLU1 and GREBLU2 were enriched in Be and Zn compared to SATBLUchem. GREBLU2 is also enriched in Li and Ga compared to the two others.

TABLE 2. Trace-element composition of three spinel samples, measured by LA-ICP-MS.
Element
(ppma)
Detection
Limit (ppma)
SKYBLU1 GREBLU2 SATBLUchem
Li 4 2778 6030 2120
Be 10 552 946 32
Ti 2 3 3 202
V 1 11 6 362
Cr 3 16 8 1111
Mn 5 238 106 287
Fe 20 11,009 9362 12,794
Co 0.2 84 22 1236
Ni 1 85 29 2514
Cu 1 4 4 4
Zn 4 7242 4887 1047
Ga 0.2 299 1088 234
 

Host Rocks. The chemical composition and characteristics of associated minerals were also examined with EDS analysis. Olivine is 99% pure forsterite (Mg2SiO4). Pargasite is rich in titanium, sodium, and chlorine. Apatites are fluorapatites with up to 20% chlorine in substitution of fluorine. In two thin sections, we analyzed one REE-rich unknown mineral and several molybdenum- and tungsten-rich unknown minerals. Cobalt was found in sulfides, as high as 1.5 wt.% in pentlandite and 3.5 wt.% in violarite.

SPECTROSCOPIC PROPERTIES OF BLUE SPINEL

UV-Vis Absorption Spectra. All UV-visible spectra showed a broad, intense absorption band between 500 and 670 nm composed of several narrower bands at about 545, 550, 560, 580, 590, and 625 nm (figure 9). Two transmission windows were seen in the visible part of the spectra, in the violet to blue region (400–500 nm) and in the red region (670–700 nm). We also observed several weak peaks between 300 and 500 nm at about 371, 386, 418, 427, 455, 460, and 480 nm. The bands at 427 and 460 nm are not visible on the spectra that show the most intense main band between 500 and 670 nm (samples SKYBLU2 and SATBLU1). Additionally, we noted a large, weak band centered at about 440 nm only on the SKYBLU samples. For samples GREBLU1 and SATBLU1, we also note an increasing absorption from 450 nm toward the UV.

UV-visible spectra of blue spinel from Vietnam
Figure 9. The UV-Visible spectra of typical blue spinels from Vietnam show a major composite absorption band between 500 and 650 nm, a large transmission window in the blue to violet region, and a smaller one in the red. The bands at 371, 386, 418, 455, 460, 480, 560, and 590 nm are due to Fe2+.The bands at 545, 550, 580, and 625 nm are due to Co2+. The band at 427 nm is not allocated.

Raman and Luminescence. The Raman spectra were typical of spinel, with weak peaks at 405, 665, and 766 cm–1 (figure 10a; Fraas et al., 1973). The 405 cm–1 peak was 9 cm–1 wide, evidence that the analyzed spinels were natural and unheated (Krzemnicki, 2008; Saeseaw et al., 2009). However, this Raman signal of spinel was overwhelmed by luminescence with the two available excitation wavelengths (514 or 647 nm). The luminescence band was centered at 107 cm–1 for the 647 nm excitation wavelength, corresponding to a 650 nm emission (figure 10b). In this case, the sample showed a strong red luminescence (figure 10b, inset) consistent with a broad band emission centered at 650 nm. In addition, many weak peaks between 673 and 710 nm, grouped in apparent triplets, were visible: 685, 687, and 689 nm; 696, 697, and 700 nm; and 704, 707, and 709 nm (figure 10b).

Raman spectra of blue spinel from Vietnam
Figure 10. The samples in this study displayed the typical Raman signal for spinel (top), with a strong continuum due to cobalt luminescence. But when the concentration in cobalt was too high, luminescence overwhelmed the signal for spinel. In some samples, spectra acquired using a 647 nm excitation showed additional luminescence peaks of chromium (the so-called organ pipe spectrum, bottom). Under green laser excitation (514 nm), the sample reacted with a red luminescence (bottom photo).

DISCUSSION

Primary Geological Origin. Red and blue spinels are always found in marble (figure 11). Garnier et al. (2008) proposed that this marble originated from an old carbonate platform (considered Precambrian to Permo-Triassic), which later metamorphosed. Graphite crystals in these Vietnamese marbles likely derive from metamorphism of organic matter (Giuliani et al., 2003; Garnier et al., 2008). As already mentioned, blue spinels are always associated with olivine (nearly pure forsterite). This paragenesis is typical of the granulitic metamorphic facies (high temperature above 550°C for a CO2-rich system; Bucher and Frey, 1994; Janardhan et al., 2001; Proyer et al., 2008). Pargasite is ubiquitous in marble, also representing a high-temperature phase. Pargasite, olivine, and spinel are nearly contemporaneous, and they may have crystallized from the destabilization of diopside with increasing pressure and temperature in a prograde reaction (Proyer et al., 2008; Ferry et al., 2011).

Blue spinel in marble host
Figure 11. Blue spinel from Luc Yen in its marble host. Photo by Vincent Pardieu/GIA.

We detected some fluorine and chlorine in apatite and pargasite, and some sodium, lithium, and beryllium in blue spinel. These elements are indications that evaporitic rocks played a role during metamorphism (Proyer et al., 2008). Giuliani et al. (1993) and Garnier et al. (2005, 2008) also proposed this hypothesis from the study of fluid inclusions in gem ruby from the Luc Yen area.

Different areas yielding ruby and red or blue spinel show distinct characteristics. Red spinel and rubies have a very similar paragenesis. Forsterite is only associated with blue spinel, and clinohumite is only associated with red spinel. Clinohumite can also grow from diopside in a prograde reaction with dolomite and water (Proyer et al., 2008). Ruby-bearing rocks are very different from those containing blue spinel, as they underwent different metamorphic histories. Because of the intense tectonic activity in Luc Yen, it is possible that two rocks with very different geological histories have been brought in contact.

Garnier et al. (2008) did not observe evidence of a fluid circulating through the marble. They proposed that aluminum and chromium originally sedimented within the carbonate platform. These elements were mobilized due to the presence of halogen elements (fluorine and chlorine) from evaporitic rocks. We propose that the mobilization of Ni and Co happened through the same process. Another hypothesis is that Ni and Co were mobilized from amphibolitic rocks interlayered in the marble (observed by Garnier et al., 2006) via halogen-rich fluids. Fluids can be formed by the metamorphism of clay minerals, evaporate, and organic matter (Giuliani et al., 2003; Garnier et al., 2008).

Proposed Geological History. The ancient Paleo-Tethys Ocean (possibly Proterozoic to Permo-Triassic) separated the China (Yangtze) and Indochina cratons (now Vietnam, Laos, Cambodia, Thailand, and Myanmar). In this ocean, a carbonaceous platform developed by sedimentation. Tectonic movements caused the closing of this ocean, and evaporitic minerals were deposited. The two main blocks (Yangtze and Indochina) subsequently collided, and all the sedimentary and magmatic rocks of the oceanic crust underwent intense deformation and metamorphism. During the collision, the metamorphism of the mix of former carbonate platform minerals and some detritic material (such as clays) deposited with it may have led to the formation of diopside through the following reactions:

tremolite + calcite → diopside + dolomite
dolomite + quartz → diopside + carbon dioxide

By increasing metamorphism, diopside destabilized into olivine, spinel, and clinohumite. The reduction of evaporitic minerals such as sulfates formed chlorine- and fluorine-rich fluids. These fluids were involved in the mobilization of aluminum and other elements such as chromium (Giuliani et al., 2003; Garnier et al., 2008).

Some processes remain poorly understood. Why do some areas show clinohumite with red spinel while others show olivine and blue spinel? What is the main difference responsible for mobilizing more chromium (red spinel) or more cobalt (blue spinel) in the marble?

Origin of Blue Color and Color Change in Vietnamese Blue Spinels. The main absorption band between 500 and 670 nm, the dominant origin of color in these blue spinels, is composed of a series of bands at approximately 545, 550, 560, 580, 590, and 625 nm (again, see figure 9). Bands at 545, 550, 580, and 625 nm are due to cobalt (Co2+) substituting for Mg2+ in tetrahedral sites of the spinel structure (Wherry, 1929; Pappalardo et al., 1961; Shigley and Stockton, 1984; Kuleshov et al., 1993; Muhlmeister et al., 1993; Delaunay et al., 2008; Duan et al., 2012; Bosi et al., 2012; D’Ippolito et al., 2015). The remaining absorption bands (at 371, 386, 418, 455, 460, 480, 560, and 590 nm) are allocated to iron (Fe2+) in tetrahedral sites of the spinel structure (Gaffney, 1973; Dickson and Smith, 1976; Muhlmeister et al., 1993; Delaunay et al., 2008; D’Ippolito et al., 2015). A weak band observed at 427 nm is not attributed but may be linked with other measurable elements such as Ni. Consequently, the spectra show transmission windows between 300 and 500 nm and between 700 and 900 nm that explain the blue color. As expected, the spectra show that iron (Fe2+) and cobalt (Co2+) are the main chromophore elements. The other trace elements detected either do not give rise to absorption in the visible range or are much less efficient absorbers than cobalt. Chromium, which is the main chromophore for red and pink spinel, makes a significant contribution to color if the concentration is above 1000 ppma. (Muhlmeister et al., 1993; T. Häger, pers. comm, 2014). Cr concentration in SATBLU samples is borderline, but the Co concentration is higher, too. The contribution is considered negligible.

We observed that the SATBLU samples, which had the most saturated color, also had the most important cobalt optical absorption. In addition, the main band had an absorption coefficient greater than 10 cm–1, and the iron optical absorptions observed were weak. For the parcel classified as SKYBLU, iron and cobalt optical absorption seemed to have a similar importance in the optical spectra, reaching a maximum of 2 cm–1 as compared to SATBLU samples. GREBLU samples had the most significant iron band, but the main absorption band only reached 0.5 cm–1. We also observed a correlation between the cobalt absorption bands and the color saturation.

SATBLU2 had a Fe/Co value of approximately 10 (table 2), and the SATBLU samples had the most saturated color (table 1). GREBLUE2, with a Fe/Co ratio of about 425 (table 2), had a visible gray hue component (table 1). For intermediate Fe/Co ratios of about 130 (measured on SKYBLU1; see table 2), the spinel had a sky-blue color (table 1). Moreover, sample GREBLU1, which had the grayest color, showed more significant bands due to Fe2+. We propose that the GREBLU samples are colored mainly by iron and the SATBLU samples by cobalt. The SKYBLU samples’ colors arise from both iron and cobalt absorption. Hue differences are more significantly controlled by iron (with different species), while saturation is largely dictated by cobalt (D’Ippolito et al., 2015).

Using chemical and spectroscopic analysis from eight of our samples, we calculated the molar absorptivity of cobalt in spinel (sensu stricto) for three absorption bands. At wavelengths of 545, 580, and 625 nm, we took the apparent maximum of each band. We determined a molar absorptivity of 530±29, 664 ±18, and 586±11 L·mol–1·cm–1, respectively. For the sake of comparison, we calculated the molar absorptivity of iron in spinel (sensu stricto) for bands 371, 386, 480, and 590 nm. These bands are attributed to ferrous iron (Fe2+) in the tetrahedral site. Our chemical analysis measured only the total iron content. Assuming all iron was in the form Fe2+, we propose that the molar absorptivity of ferrous iron in the tetrahedral site had an order of magnitude of about 30 L·mol–1·cm–1 for each band. With this method, values of molar absorptivity are not very accurate but provide a working assumption for our preliminary study. We recognize that further investigation is needed to fully understand the color in blue spinel. A Gaussian decomposition of spectra can improve the precision of these values. Nevertheless, it is apparent that in spinel, Co2+ is approximately 20 times more efficient at absorbing light, and thus creating color, than Fe2+ (consistent with D’Ippolito et al., 2015).

Spectra have two transmission windows between 350 and 500 nm (in the blue region) and between 670 and 900 nm (in the red region). This explains the pink to red reaction under the Chelsea filter and the color change. Indeed, the Chelsea filter probes a transmission window in the red. The color change is also explained when the spectral composition of the lighting environment is compared with the absorption spectra of spinel, although this change is not observed in every example. Compact fluorescent light emits more in the blue region than in the red, and therefore the spinel appears blue. Under incandescent light, which is richer in red, spinel displays a violetish blue color that is mostly blue with minor red.

Origin of Red Luminescence. Under laser excitation, our samples showed a strong red luminescence. In spectra acquired using a 647 nm excitation, we observed several peaks (in groups of three) between 673 and 710 nm (at about 673, 675, 685, 687, 689, 696, 697, 700, 704, 707, and 709 nm; see figure 10b). These peaks are known to be due to trivalent chromium (Cr3+) substituting for aluminum in the octahedral site (Burns et al., 1965; Wood et al., 1968; Skvortsova et al., 2011). The broad band centered at 650 nm is allocated to divalent cobalt in the tetrahedral site of the spinel structure (Abritta and Blak, 1991; Kuleshov et al., 1993). These luminescence behaviors are consistent with our chemical analysis, as the strongest luminescence was observed in the SATBLU samples, which had higher concentrations of Cr3+ (1111 ppma) and Co2+ (1236 ppma).

CONCLUSION

We confirmed that the saturated “cobalt-blue” color of Vietnamese spinels (figures 12 and 13) is due predominantly to Co2+ substituting for Mg2+ in the tetrahedral site of the spinel structure. For the most saturated blue spinel, cobalt is the main coloring agent, even if iron is more abundant. Indeed, cobalt is a powerful coloring agent, with a molar absorptivity between 500 and 700 L·mol–1·cm–1 depending on wavelength, whereas iron (Fe2+ in the tetrahedral site) has a molar absorptivity of about 30 L·mol–1·cm–1. The higher the iron/cobalt ratio is, the grayer the color. The red transmission window of these gems explains both their pink to red Chelsea filter reaction and their slight change of color from blue to “lavender” with a change of lighting environment. The red luminescence is due to both Cr3+ and Co2+, and it may have a minor influence on the perceived color.

Examining a blue spinel crystal
Figure 12. A gem merchant examines a blue spinel from the Luc Yen district. Photo by Vincent Pardieu/GIA.
2.59 ct cobalt blue spinel
Figure 13. Vietnam’s spinel production yielded this 2.59 ct cobalt blue gem. Photo by Robert Weldon/GIA, courtesy of Palagems.com.

This study offers clues to the definition of “cobalt-blue” spinel. Cobalt is actually the main chromophore, but the presence of iron is also significant. The term “cobalt-blue” can be clarified by further investigations on the significance of each chromophore elements (iron and cobalt). These investigations can propose a limit on the ratio of iron/cobalt above which the term “cobalt-blue” cannot be used.

Spinels from Luc Yen contain few inclusions. Fractures and fingerprints were the most common inclusions found. Sometimes, we observed parallel elongated tubes with black, irregular solid inclusions associated.

From a geological standpoint, gem-quality blue spinels are associated with intense metamorphism. Their marble host results from the metamorphism of an ancient carbonaceous platform. This platform was located in the Paleo-Tethys Ocean, which separated Indochina and China. During the convergence of these “paleo-continents,” the ocean closed off, accompanied by the formation of evaporitic rocks. The ocean crust, associated with the carbonaceous platform and evaporitic platform, was sandwiched between the two continents. The collision led to the metamorphism of the evaporite rocks, in turn producing fluids mobilizing some elements, possibly including cobalt. Spinel grew in the marble during this intense metamorphism. These processes of metamorphism and fluid interaction led to the crystallization of attractive blue spinels in the marble mountains of Luc Yen.

Mr. Chauviré (boris.chauvire@univ-nantes.fr) is a PhD student at the Laboratoire de Planétologie et Géodynamique de Nantes, France, and Dr. Rondeau is an assistant professor at the same laboratory (CNRS Team 6112). Dr. Fritsch (CNRS Team 6502) is a professor of physics at the University of Nantes, Institut des Matériaux Jean Rouxel. Mr. Ressigeac is product manager for Montepuez Ruby Mining, Mozambique. Mr. Devidal is an engineer specialist of ICP-MS-LA at the Laboratoire Magmas et Volcans, Clermont Ferrand, France.

We are grateful to Vincent Pardieu, senior manager of field gemology at GIA’s Bangkok laboratory, for his valuable aid during the preparation and progress of the expeditions. We thank Vincent’s contacts for acquiring blue spinel samples. We are grateful to Mr. Chuãn, our guide, for his knowledge of the field and his logistical support. We also thank Pham Van Long, director of the Center for Gem and Gold Research and Identification in Hanoi, for his logistical support and for exporting the samples collected. GIA’s laboratory in Bangkok and its director, Kenneth Scarratt, provided technical and logistical support. Jean-Pierre Lorand (LPGN-CNRS) generously shared his knowledge about sulfides. We thank Alexandre Droux from the Laboratoire Français de Gemmologie for EDXRF measurements. We also thank Tobias Häger of Johannes Gutenberg University in Mainz, Germany, for his help with interpreting UV-Vis spectra. We also thank reviewers that participate to improve this study.

Abritta T., Blak F.H. (1991) Luminescence study of ZnGa2O4:Co. Journal of Luminescence, Vol. 48–49, pp. 558–560.
 
Anczkiewicz R., Viola, G., Müntener O., Thirlwall M.F., Villa I.M., Qong N.Q. (2007) Structure and shearing conditions in the Day Nui Con Voi massif: Implications for the evolution of the Red River shear zone in northern Vietnam. Tectonics, Vol. 26, pp. 1–21.
 
Bosi F., Hålenius U., D’Ippolito V., Andreozzi G.B. (2012) Blue spinel crystals in the MgAl2O4-CoAl2O4 series: Part II. Cation ordering over short-range and long-range scales. American Mineralogist, Vol. 97, No. 11–12, pp. 1834–1840, http://dx.doi
.org/10.2138/am.2012.4139.
 
Bucher K, Frey M. (1994) Petrogenesis of Metamorphic Rocks, 6th ed. Springer-Verlag, 308 pp.
 
Burns G., Geiss E. A., Jenkins B.A., Nathan M. I. (1965) Cr3+ fluorescence in garnet and other crystals. Physical Review A, Vol. 139, No. 5A, pp. A1687–A1693, http://dx.doi.org/10.1103/PhysRev.139.A1687.
 
Delaunay A. (2008) Les spinelles bleus. http://www.geminterest.
com/.
 
Dickson B.L., Smith G. (1976) Low-temperature optical absorption and Mössbauer spectra of staurolite and spinel. Canadian Mineralogist, Vol. 14, pp. 206–215.
 
D’IppolitoV., Andreozzi G. B., Hålenius U., Skogby H., Hametner K., Günther D. (2015) Color mechanisms in spinels: cobalt and iron interplay for the blue color. Physics and Chemistry of Materials, http://dx.doi.org/10.1007/s00269-015-0734-0.
 
Duan X., Wang X., Yu F., Liu X. (2012) Effects of Co content and annealing temperature on the structure and optical properties of CoxMg1–xAl2O4 nanoparticles. Materials Chemistry and Physics, Vol. 137, No. 2, pp. 652–659, http://dx.doi.org/10.1016/j.matchemphys.2012.10.016.
 
Ferry J.M., Ushikubo T., Valley J.W. (2011) Formation of forsterite by silicification of dolomite during contact metamorphism. Journal of Petrology, Vol. 52, No. 9, pp. 1619–1640, http://dx.doi.org/10.1093/petrology/egr021.
 
Fraas L.M., Moore J.E., Salzberg J.B. (1973) Raman characterization studies of synthetic and natural MgAl2O4 crystals. The Journal of Chemical Physics, Vol. 58. No. 9. pp 3585–3592, http://dx.doi.org/10.1063/1.1679704.
 
Gaffney E.S. (1973) Spectra of tetrahedral Fe2+ in MgA12O4. Physical Review B, Vol. 8, No. 7, pp. 3484–3486.
 
Garnier V., Giuliani G., Maluski H., Ohnenstetter D., Phan Trong T., Hoàng Quang V., Pham Van L., Vu Van T., Schwarz D. (2002) Ar-Ar ages in phlogopites from marble-hosted ruby deposits in northern Vietnam: evidence for Cenozoic ruby formation. Chemical Geology, Vol. 188, pp. 33–49, http://dx.doi.org/10.1016/S0009-2541(02)00063-3.
 
Garnier V., Ohnenstetter D., Giuliani G., Maluski H., Deloule E., Phan Trong, T., Pham Van L., Hoàng Quang V. (2005) Age and significance of ruby-bearing marble from the Red River Shear Zone, northern Vietnam. The Canadian Mineralogist, Vol. 43, pp. 1315–1329.
 
Garnier V., Giuliani G., Ohnenstetter D., Fallick A.E., Dubessy J., Banks D., Vinh H.Q., Lhomme T., Maluski H., Pêcher A., Bakhsh K.A., Long P.V., Trinh P.T., Schwarz D. (2008) Marble-hosted ruby deposits from Central and South-East Asia: Toward a new genetic model. Ore Geology Reviews, Vol. 34, pp. 169–191, http://dx.doi.org/10.1016/j.oregeorev.2008.03.003.
 
Giuliani G., Hoàng Quang V., Phan Trong T., France-Lanord C., Coget P. (1999) Carbon isotopes study on graphite and coexisting calcite-graphite pairs in marbles from the Luc Yen and Yen Bai district, North of Vietnam. Bulletin de Liaison S.F.M.C., Vol. 11, pp. 80–82.
 
Giuliani G., Dubessy J., Banks D., Hoàng Quang V., Lhomme T., Pironon J., Garnier V., Phan Trong T., Pham Van L., Ohnenstetter D., Schwarz D. (2003) CO2-H2S-COS-S8-AlO(OH)-bearing fluid inclusions in ruby from marble-hosted deposits in Luc Yen area, North Vietnam. Chemical Geology, Vol. 194, pp. 167–165, http://dx.doi.org/10.1016/S0009-2541(02)00276-0.
 
Hauzenberger C.A., Häger T., Hofmeister W, Quang V.X., Rohan Fernando G.W.A. (2003) Origin and formation of gem quality corundum from Vietnam. Geo- and Materials-Science on Gem-Minerals of Vietnam, Proceedings of the International Workshop, Hanoi, October 1–8.
 
Huong L.T.-T., Häger T., Hofmeister W., Hauzenberger C., Schwarz D., Long P.V., Wehmeister U., Khoi N.N., Nhung N.T. (2012) Gemstones from Vietnam: An update. G&G, Vol. 48, No. 3, pp. 158–176, http://dx.doi.org/10.5741/GEMS.48.3.158.
 
Janardhan A.S., Sriramguru K., Basava S., Shankara M.A. (2001) Geikielite-Mg.Al.Spinel-titanoclinohumite association from a marble quarry near Rajapalaiyam area, part of the 550 Ma of southern granulite terrain, southern India. Gondwana Research, Vol. 4, No. 3, pp. 359–366, http://dx.doi.org/10.1016/S1342-937X(05)70335-X.
 
Jolivet L., Beyssac O., Goffé B., Avigad D., Lepvrier C., Maluski H., Than T.T. (2001) Oligo-miocène midcrustal subhorizontal shear zone in Indochina. Tectonics, Vol. 20, No. 1, pp. 46–57.
 
Krzemnicki M.S. (2008) Trade Alert: Flux grown synthetic spinels again on the market. SSEF Newsletter, October.
 
Kuleshov N.V., Mikhailov V.P., Sherbitsky V.G., Prokoshin P.V., Yumashev K.V. (1993) Absorption and luminescence of tetrahedral Co2+ ion in MgAl2O4. Journal of Luminescence, Vol. 55, pp. 265–269, http://dx.doi.org/10.1016/0022-2313(93)
90021-E.
 
Kušnír I. (2000) Mineral resources of Vietnam. Acta Montanistica Slovaca, Vol. 5, No. 2, pp. 165–172.
 
Leloup P.H., Lacassin R., Tapponnier P., Schärer U., Dalai Z., Xiaohan L., Liangshang Z., Shaocheng J, Phan Trong T. (1995) The Ailao Shan-Red River shear zone (Yunnan, China), Tertiary transform boundary of Indochina. Tectonophysics, Vol. 251, pp. 3–84, http://dx.doi.org/10.1016/0040-1951(95)00070-4.
 
Leloup P.H., Arnaud N., Lacassin R., Kienast J.R., Harrison T.M., Phan Trong T.T., Replumaz A., Taponnier P. (2001) New constraints on the structure, thermochronology, and timing of the Ailo Shan-Red River shear zone, SE Asia. Journal of Geophysical Research, Vol. 106, No. B4, pp. 6683–6732.
 
Lepvrier C., Vuong N.V., Maluski H., Thi P.T., Vu T.V. (2008) Indosinian tectonics in Vietnam. Comptes Rendus Geoscience, Vol. 340, No. 2–3, pp. 94–111, http://dx.doi.org/10.1016/j.crte.2007.10.005.
 
Long P.V., Giuliani G., Garnier V., Ohnenstetter D. (2004) Gemstones in Vietnam. The Australian Mineralogist, Vol. 22, No. 4. pp. 162–168.
 
Manson D.V., Stockton C.M. (1984) Pyrope-spessartine garnets with unusual color behavior. G&G, Vol. 20, No. 4, pp. 200–207, http://dx.doi.org/10.5741/GEMS.20.4.200.
 
Muhlmeister S., Koivula J.I., Kammerling R.C., Smith C.P., Fritsch E., Shigley J.E. (1993) Flux-grown synthetic red and blue spinels from Russia. G&G, Vol. 29, No. 2, pp. 81–98, http://dx.doi.org/10.5741/GEMS.29.2.81.
 
Nassau K. (1980) Gems Made by Man. Chilton Book Co., Radnor, PA, 364 pp.
 
Overton T.W., Shen A.H. (2011) Gem News International: Cobalt blue-colored spinel from Khuoi Ngan, Vietnam. G&G, Vol. 47, No. 4, pp. 328–329.
 
Pappalardo R., Wood D.L., Linares R.C. (1961) Optical absorption study of Co-doped oxide systems. II. The Journal of Chemical Physics, Vol. 35, No. 6, pp. 2041–2059, http://dx.doi.org/10.1063/1.1732208.
 
Pardieu V. (2012) FE34 Vietnam blue spinels. http://www.field
gemology.org
 
Pardieu V., Hughes R.W. (2008) Spinel: The resurrection of a classic. InColor, Summer, pp. 10–12 (also on www.ruby-sapphire.com and www.fieldgemology.org).
 
Proyer A., Mposkos E., Baziotis I., Hoinkes G. (2008) Tracing high-pressure metamorphism in marbles: Phase relations in high-grade aluminous calcite-dolomite marbles from the Greek Rhodope massif in the system CaO-MgO-Al2O3-SiO2-CO2 and indications of prior aragonite. Lithos, Vol. 104, No. 1–4, pp. 119–130, http://dx.doi.org/10.1016/j.lithos.2007.12.002.
 
Rice J.M. (1980) Phase equilibria involving humite minerals in impure dolomitic limestones: part I. Calculated stability of clinohumite. Contributions to Mineralogy and Petrology, Vol. 71, No. 3, pp. 219–235, http://dx.doi.org/10.1007/BF00371664.
 
Saeseaw S., Wang W., Scarratt K., Emmett J.L., Douthit T.R. (2009) Distinguishing heated spinels from unheated natural spinels and from synthetic spinels. http://www.gia.edu/gia-news-research-NR32209A.
 
Satish-Kumar M. (1999) An overview of petrology of calc-silicate granulites from the Trivandrum block, southern India. Journal of Geosciences, Vol. 42. Osaka City University, pp. 127–159.
 
Senoble J.B. (2010) Beauty and rarity – A quest for Vietnamese blue spinels. InColor, Summer, pp. 2–7.
 
Shigley J.E., Stockton C.M. (1984) “Cobalt-blue” gem spinels. G&G, Vol. 20, No. 1, pp. 34–41, http://dx.doi.org/10.5741/GEMS.20.1.34.
 
Skvortsova V., Mironova-Ulmane N., Riekstina D. (2011) Structure and phase changes in natural and synthetic magnesium aluminum spinel. Proceedings of the 8th International Scientific and Practical Conference, Vol. II, pp. 100–106.
 
Tapponnier P., Peltzer G., Le Dain A.Y., Armijo R., Cobbold P. (1982) Propagating extrusion tectonics in Asia: New insights from simple experiments with plasticine. Geology, Vol. 10, No. 12, pp. 611–616, http://dx.doi.org/10.1130/0091-7613(1982)10611:PETIAN2.0.CO;2.
 
Tapponnier P., Lacassin R., Leloup P.H., Shärer U., Dalai Z., Haiwei W., Xiaohan Liu, Shaocheng J., Lianshang Z., Jiayou Z. (1990) The Ailao Shan/Red River metamorphic belt: Tertiary left-lateral shear between Indochina and South China. Nature, Vol. 343, No. 6257, pp. 431–437, http://dx.doi.org/10.1038/343431a0.
 
Ulmer P., Luth R.W. (1991) The graphite-COH fluid equilibrium in P, T, fO2 space. Contributions to Mineralogy and Petrology, Vol. 106, pp. 265–272.
 
Webster R. (1994) Gems: Their Sources, Descriptions and Identification, 5th ed. Butterworth-Heinemann, Oxford, edited by P.G. Read, 1026 pp.
 
Wherry E.T. (1929) Mineral determination by absorption spectra II. The American Mineralogist, Vol. 14, No. 9, pp. 323–328.
 
Wood D.L., Imbusch G.F., Macfarlane R.M., Kisliuk P., Larkin D.M. (1968) Optical spectrum of Cr3+ ions in spinels. The Journal of Chemical Physics, Vol. 48, No. 11, pp. 5255–5263, http://dx.doi.org/10.1063/1.1668202.
 
Yu J., Mao J., Chen F., Wang Y., Che L., Wang T., Liang J. (2014) Metallogeny of the Shilu Fe-Co-Cu deposit, Hainan Island, South China: Constraints from fluid inclusions and stable isotopes. Ore Geology Reviews, Vol. 57, pp. 351–362, http://dx.doi.org/10.1016/j.oregeorev.2013.08.018.