Variety of Iron Oxide Inclusions in Sapphire from Southern Vietnam: Indication of Environmental Change during Crystallization

: Sapphires from alluvial deposits associated with Cenozoic basalts in Southern Vietnam were collected for investigation of mineral inclusions. In this report, primary iron oxide inclusions were focused on, with detailed mineral chemistry using a Raman spectroscope and electron probe micro-analyzer. Consequently, a variety of iron oxide inclusions were recognized as wüstite, hercynite, and ilmenite. Ilmenite falling within an ilmenite–hematite series ranged in composition between Il 24-30 He 36-38 Mt 35-40 and Il 49-54 He 34-40 Mt 7-10 , classified as titanomagnetite and titanohematite, respectively. Wüstite with non-stoichiometry, (Fe 2+0.3-0.9 )(Ti 3+<0.179 Al 3+ ≤ 0.6 Cr 3+<0.1 Fe 3+ ≤ 0.46 ) (cid:3) ≤ 0.23 O, was associated with hercynite inclusions, clearly indicating cogenetic sapphire formation. Wüstite and sapphire appear to have been formed from the breakdown reaction of hercynite (hercynite = sapphire+wüstite) within a reduction magma chamber. Titanohematite and titanomagnetite series might have crystallized during iron–titanium reequilibration via subsolidus exsolution under a slightly oxidized cooling process.


Introduction
Iron oxide minerals have been considered as a significant geothermometer for their host rocks [1][2][3][4][5][6]. Previous investigations of gem sapphire from Southern Vietnam, such as Dak Nong, Di Linh, and Binh Thuan deposits, have reported several iron oxide inclusions [7][8][9]. Most of these iron oxides were identified as ilmenite, magnetite-hercynite, and chromite-hercynite, using a scanning electron microscope-energy dispersive spectrometer (SEM-EDS) and an X-ray diffractometer (XRD) [8]. Moreover, Izokh et al. [7] reported an iron oxide inclusion, namely, Al-Ti-hematite, that was chemically analyzed by an electron probe micro-analyzer (EPMA); subsequently, they proposed that crystallization of host sapphire should relate to iron-rich syenitic melt and metasomatism between crustal rocks and contaminated basaltic melt in the Dak Nong deposit. Recently, Vu et al. [9] reported spinel and other unidentified oxide inclusions in sapphires from many deposits in Southern Vietnam which more details and further investigation are reported herein this manuscript. Besides, the crystallization of sapphire and related zircon from Southern Vietnam may occur in the lithospheric mantle which is related to carbonatite-dominant melts as a result of partial melting of a metasomatized lithospheric mantle source, at over 900 • C [10].
Although iron oxide inclusions were previously reported, their mineral chemistry has never been fully analyzed. This study was therefore designed to analyze most types of iron oxide inclusion in sapphire from Dak Nong, Di Linh, Binh Thuan, Krong Nang, and Pleiku.

Geological Setting
Southern Vietnam geologically consists of a large-scale structure of the Da Lat active continental margin (Da Lat zone), Indosinian polyepisodic orogenic belt (Srepok orogenic belt), and Kontum massif [11], which is part of the Indosinian craton [12] (Figure 1). This region is composed of Archean-Proterozoic basement rocks, Early to Middle Paleozoic cover rocks, Jurassic sediments, late Mesozoic rocks, and Cenozoic basaltic rocks [12] (Figure 2). The basement rocks in this area are characterized by metamorphic complexes of granulites which are usually covered by volcanogenic sedimentary rocks, metamorphosed sedimentary rocks within greenschist facies of the Early to Middle Paleozoic, as well as sandstone, siltstone, and shale of Jurassic sedimentary formation. These basements and upper rocks are intruded by a number of late Mesozoic igneous rocks, including Triassic granite, granodiorite, and granosyenite (results of the Indosinian-Yangtze collision during the Permo-Triassic about 245-240 Ma) and Cretaceous diorite and granodiorite (result of the Paleo-Pacific plate subduction) [13][14][15]. Late Cenozoic basalts associated with sapphire deposits in Southern Vietnam have been mapped, overlining the older rock formations reported above.  [11] and Hoa et al. [12].
These Late Cenozoic basalts ( Figure 2) are related to regional tectonic terranes, particularly after the end of East Sea opening in the Middle Miocene [16][17][18]. Paleo-Pacific oceanic crustal material, previously subducted into the lower mantle, was subsequently entrained into the Hainan plume which was the main cause of basaltic magmatism in Southern Vietnam [19]. According to Hoang and Flower [20], these volcanic activities were extended more than 100 km in diameter with a thickness up to several hundred meters and covered a total area of approximately 23,000 km 2 . The centers of volcanism appear to have developed during two main eruptive episodes. The early phases generated mainly quartz and olivine tholeiites with rare alkali basalt whereas the later phases produced  [11] and Hoa et al. [12].
These Late Cenozoic basalts ( Figure 2) are related to regional tectonic terranes, particularly after the end of East Sea opening in the Middle Miocene [16][17][18]. Paleo-Pacific oceanic crustal material, previously subducted into the lower mantle, was subsequently entrained into the Hainan plume which was the main cause of basaltic magmatism in Southern Vietnam [19]. According to Hoang and Flower [20], these volcanic activities were extended more than 100 km in diameter with a thickness up to several hundred meters and covered a total area of approximately 23,000 km 2 . The centers of volcanism appear to have developed during two main eruptive episodes. The early phases generated mainly quartz and olivine tholeiites with rare alkali basalt whereas the later phases produced olivine tholeiite, alkali basalt, basanite, and rare nephelinite. Tholeiite eruptions occurred significantly in the centers associated with extensional rifting. On the other hand, alkali basalt, olivine tholeiite, and basanite appear to have erupted along the conjugate strike-slip faults [18]. It should be noted that sapphire and zircon occurrences mainly discovered in Quaternary and Upper-Pleistocene alluvial deposits derived from the alkali basalts [8][9][10]21]. They would originate in the deep-seated formations before being transported as megacryst onto the Earth's surface by alkali basaltic magmas.  [18]. It should be noted that sapphire and zircon occurrences mainly discovered in Quaternary and Upper-Pleistocene alluvial deposits derived from the alkali basalts [8][9][10]21]. They would originate in the deep-seated formations before being transported as megacryst onto the Earth's surface by alkali basaltic magmas.

Materials and Methods
Sapphire collections of about 5 mm crystals were sampled from Dak Nong, Di Linh, Binh Thuan, Krong Nang, and Pleiku in Southern Vietnam (see Figure 2 and Table 1 and Table S1). The sapphire crystals were detected for metallic opaque inclusions for further study. After being mounted in epoxy, they were ground by a diamond wheel until inclusions were exposed; subsequently, they were polished using 6 µm, 3 µm, and 1 µm diamond pastes. These inclusions were initially characterized using a laser Raman spectroscope, inVia model, Renishaw equipped with a Leica optical microscope at the Gem and Jewelry Institute of Thailand (Public Organization) (GIT). The laser beam was generally set at about 5µm with standard conditions including 532 nm radiation of the NIR diode laser emitting at 785 nm and power of 15.7 mW (about 5 mW on the sample), at a resolution of approximately 0.5 cm −1 . It should be noted that Raman patterns of iron minerals may be transformed rapidly by laser induction, leading to the modification of Raman shift, and the same effect can also be induced by the natural processes such as oxidation, recrystallization, order-disorder transitions (cation redistribution), phase transition, or decomposition [22][23][24]. Therefore, low laser power (0.5 mW) and normal operation power (5

Materials and Methods
Sapphire collections of about 5 mm crystals were sampled from Dak Nong, Di Linh, Binh Thuan, Krong Nang, and Pleiku in Southern Vietnam (see Figure 2 and Table 1 and Table S1). The sapphire crystals were detected for metallic opaque inclusions for further study. After being mounted in epoxy, they were ground by a diamond wheel until inclusions were exposed; subsequently, they were polished using 6 µm, 3 µm, and 1 µm diamond pastes. These inclusions were initially characterized using a laser Raman spectroscope, inVia model, Renishaw equipped with a Leica optical microscope at the Gem and Jewelry Institute of Thailand (Public Organization) (GIT). The laser beam was generally set at about 5µm with standard conditions including 532 nm radiation of the NIR diode laser emitting at 785 nm and power of 15.7 mW (about 5 mW on the sample), at a resolution of approximately 0.5 cm −1 . It should be noted that Raman patterns of iron minerals may be transformed rapidly by laser induction, leading to the modification of Raman shift, and the same effect can also be induced by the natural processes such as oxidation, recrystallization, order-disorder transitions (cation redistribution), phase transition, or decomposition [22][23][24]. Therefore, low laser power (0.5 mW) and normal operation power (5 mW) on the sample were applied in this study to observe the Raman pattern and its alteration in wüstite inclusions because its structure is sensitive to laser power. However, ilmenite and hercynite spinel were analyzed using normal operation laser power of 5 mW, as suggested by Wang et al. [25]. Each spectrum was recorded within the spectral range of 200 cm −1 to 1500 cm −1 for 20 s of exposure time, 6 accumulation, and 50× magnification, at a laboratory temperature of about 22 • C.
After identification using the Raman technique, these inclusion samples were carbon coated prior to major and minor analyses using an electron probe microanalyzer (EPMA, JEOL model JXA-8100) at the Department of Geology, Faculty of Science, Chulalongkorn University. The analytical condition was set at 15 kV acceleration voltages with about 24 nA probe current. Appropriate standards and analytical crystals were selected for analyses with 30 s for peak counts and background counts of each element before automatic ZAF correction, requested for three main effects (i.e., atomic number, absorption, and fluorescence excitation) influencing spectroscopic analysis of characteristic X-rays, was applied to report in oxide contents. Finally, atomic proportions of these analyses were recalculated on the basis of proper oxygen in their formular before Fe 2+ and Fe 3+ ratios were estimated accordingly using the equation of Droop [26].

Results
Although all of the studied iron oxide minerals looked similarly black opaque crystals, their morphological affinities were classified into two groups, octahedral and rhombohedral shapes, indicating significantly primary inclusions. Octahedral iron oxides formed as single and aggregate crystals ranging in size from 20-500 µm ( Figure 3). In contrast, rhombohedral iron oxides occurred typically as single crystals, about 80-600 µm long, and 40-300 µm wide ( Figure 4). Raman spectroscopic identification of these iron oxide inclusions is reported below. mW) on the sample were applied in this study to observe the Raman pattern and its alteration in wüstite inclusions because its structure is sensitive to laser power. However, ilmenite and hercynite spinel were analyzed using normal operation laser power of 5 mW, as suggested by Wang et al. [25]. Each spectrum was recorded within the spectral range of 200 cm −1 to 1500 cm −1 for 20 s of exposure time, 6 accumulation, and 50× magnification, at a laboratory temperature of about 22 °C.
After identification using the Raman technique, these inclusion samples were carbon coated prior to major and minor analyses using an electron probe microanalyzer (EPMA, JEOL model JXA-8100) at the Department of Geology, Faculty of Science, Chulalongkorn University. The analytical condition was set at 15 kV acceleration voltages with about 24 nA probe current. Appropriate standards and analytical crystals were selected for analyses with 30 s for peak counts and background counts of each element before automatic ZAF correction, requested for three main effects (i.e., atomic number, absorption, and fluorescence excitation) influencing spectroscopic analysis of characteristic X-rays, was applied to report in oxide contents. Finally, atomic proportions of these analyses were recalculated on the basis of proper oxygen in their formular before Fe 2+ and Fe 3+ ratios were estimated accordingly using the equation of Droop [26].

Results
Although all of the studied iron oxide minerals looked similarly black opaque crystals, their morphological affinities were classified into two groups, octahedral and rhombohedral shapes, indicating significantly primary inclusions. Octahedral iron oxides formed as single and aggregate crystals ranging in size from 20-500 µm ( Figure 3). In contrast, rhombohedral iron oxides occurred typically as single crystals, about 80-600 µm long, and 40-300 µm wide ( Figure 4). Raman spectroscopic identification of these iron oxide inclusions is reported below.   Three distinctive types were recognized by Raman spectra, including wüstite, hercynite, and ilmenite groups. In comparison with morphological features, wüstite and hercynite inclusions typically showed octahedral shape in forms of single and aggregate crystals. Moreover, wüstite-hercynite composite inclusions ( Figure 5) could be found. On the other hand, ilmenite inclusions were mostly characterized by rhombohedral single crystals.
At the low laser power (0.5 mW), all studied wüstite inclusions clearly showed only a sharp peak at 670 cm −1 , which closely resembled magnetite spectra ( Figure 6a). On the other hand, at high laser power (5 mW), some wüstite inclusions (e.g., DL50, DL56, PT17) showed characteristic patterns of magnetite (weaker peak in the range of 650-670 cm −1 ), hematite with higher intensity peaks at 218, 285, and 388, and wüstite assigned by the 595 cm −1 peak, as suggested by Hanesch [23] (Figure 6b). These spectra were matched well with wüstite spectra observed at low and high laser powers reported by Faria et al. [22]. In addition, some wüstite inclusions showed a broader band in the 650-670 cm −1 region (strongest peak at approximately 667 cm −1 ) of magnetite with higher intensity (Figure 6c). The Raman spectrum of hercynite is representatively displayed in Figure 7. The spectrum yielded a strong band at 753 cm −1 and a weak peak at 701 cm −1 which indicate the Three distinctive types were recognized by Raman spectra, including wüstite, hercynite, and ilmenite groups. In comparison with morphological features, wüstite and hercynite inclusions typically showed octahedral shape in forms of single and aggregate crystals. Moreover, wüstite-hercynite composite inclusions ( Figure 5) could be found. On the other hand, ilmenite inclusions were mostly characterized by rhombohedral single crystals. Three distinctive types were recognized by Raman spectra, including wüstite, hercynite, and ilmenite groups. In comparison with morphological features, wüstite and hercynite inclusions typically showed octahedral shape in forms of single and aggregate crystals. Moreover, wüstite-hercynite composite inclusions ( Figure 5) could be found. On the other hand, ilmenite inclusions were mostly characterized by rhombohedral single crystals.
At the low laser power (0.5 mW), all studied wüstite inclusions clearly showed only a sharp peak at 670 cm −1 , which closely resembled magnetite spectra (Figure 6a). On the other hand, at high laser power (5 mW), some wüstite inclusions (e.g., DL50, DL56, PT17) showed characteristic patterns of magnetite (weaker peak in the range of 650-670 cm −1 ), hematite with higher intensity peaks at 218, 285, and 388, and wüstite assigned by the 595 cm −1 peak, as suggested by Hanesch [23] (Figure 6b). These spectra were matched well with wüstite spectra observed at low and high laser powers reported by Faria et al. [22]. In addition, some wüstite inclusions showed a broader band in the 650-670 cm −1 region (strongest peak at approximately 667 cm −1 ) of magnetite with higher intensity (Figure 6c). The Raman spectrum of hercynite is representatively displayed in Figure 7. The spectrum yielded a strong band at 753 cm −1 and a weak peak at 701 cm −1 which indicate the At the low laser power (0.5 mW), all studied wüstite inclusions clearly showed only a sharp peak at 670 cm −1 , which closely resembled magnetite spectra (Figure 6a). On the other hand, at high laser power (5 mW), some wüstite inclusions (e.g., DL50, DL56, PT17) showed characteristic patterns of magnetite (weaker peak in the range of 650-670 cm −1 ), hematite with higher intensity peaks at 218, 285, and 388, and wüstite assigned by the 595 cm −1 peak, as suggested by Hanesch [23] (Figure 6b). These spectra were matched well with wüstite spectra observed at low and high laser powers reported by Faria et al. [22]. In addition, some wüstite inclusions showed a broader band in the 650-670 cm −1 region (strongest peak at approximately 667 cm −1 ) of magnetite with higher intensity (Figure 6c). peak, 683 cm −1 , was due to ilmenite, as suggested by Wang et al. [25], whereas 299 and 498 cm −1 peaks were caused by hematite [22]. Additionally, some ilmenite inclusions revealed broadened peaks at 220, 293, and 613 cm −1 which matched well with the hematite's pattern; moreover, an additional characteristic peak of ilmenite was present around 683 cm -1 as well as weak peaks at 399, 600, and 1298 cm −1 , likely indicating oxidized titanomagnetite (Figure 8b) [24,25]. These different Raman patterns indicated that ilmenite inclusions in this study may have a variety of chemical compositions.   The Raman spectrum of hercynite is representatively displayed in Figure 7. The spectrum yielded a strong band at 753 cm −1 and a weak peak at 701 cm −1 which indicate the vibration of AlO 4 tetrahedra, characteristic of a spinel structure, as suggested by Cynn et al. [27]. This spectrum was similar to the hercynite spinel reported by Ospitali et al. [28].
vibration of AlO4 tetrahedra, characteristic of a spinel structure, as suggested by Cynn et al. [27]. This spectrum was similar to the hercynite spinel reported by Ospitali et al. [28].
An ilmenite inclusion showed Raman characteristics of mixture minerals of which Raman peaks were recognized at 299, 498, and 683 cm −1 (Figure 8a). The highest intensity peak, 683 cm −1 , was due to ilmenite, as suggested by Wang et al. [25], whereas 299 and 498 cm −1 peaks were caused by hematite [22]. Additionally, some ilmenite inclusions revealed broadened peaks at 220, 293, and 613 cm −1 which matched well with the hematite's pattern; moreover, an additional characteristic peak of ilmenite was present around 683 cm -1 as well as weak peaks at 399, 600, and 1298 cm −1 , likely indicating oxidized titanomagnetite (Figure 8b) [24,25]. These different Raman patterns indicated that ilmenite inclusions in this study may have a variety of chemical compositions.   An ilmenite inclusion showed Raman characteristics of mixture minerals of which Raman peaks were recognized at 299, 498, and 683 cm −1 (Figure 8a). The highest intensity peak, 683 cm −1 , was due to ilmenite, as suggested by Wang et al. [25], whereas 299 and 498 cm −1 peaks were caused by hematite [22]. Additionally, some ilmenite inclusions revealed broadened peaks at 220, 293, and 613 cm −1 which matched well with the hematite's pattern; moreover, an additional characteristic peak of ilmenite was present around 683 cm −1 as well as weak peaks at 399, 600, and 1298 cm −1 , likely indicating oxidized titanomagnetite (Figure 8b) [24,25]. These different Raman patterns indicated that ilmenite inclusions in this study may have a variety of chemical compositions. Mineral chemistry of these iron oxide inclusions, based on the EPMA, clearly supported Raman spectroscopic identification. Representative EPMA analyses of each iron oxide type are shown in Tables 1-3.  Mineral chemistry of these iron oxide inclusions, based on the EPMA, clearly supported Raman spectroscopic identification. Representative EPMA analyses of each iron oxide type are shown in Tables 1-3.   (Figure 9). It should be noted again that some hercynite was found in close association with wüstites in sapphire from Binh Thuan (sample PT55, as shown in Figure 5, with EPMA analyses presented in Tables 1 and 2). Note: * previously reported by Vu et al. [9]. Fe 2+ and Fe 3+ were recalculated from total FeO after the method of Droop [26]. ΣR 2+ = Fe 2+ + Mn + Mg + Zn + Ca + Ni. ΣR 3+ = Si + Ti + Al + Cr + Fe 3+ .

Classification of Iron Oxide Inclusions
Based on morphological and Raman spectroscopic characteristics combined with EPMA analyses, iron oxide minerals in sapphire from Southern Vietnam can be grouped into wüstite, hercynite, and ilmenite. These EPMA analyses are newly reported for the mineral chemistry of these oxide inclusions in sapphires from Southern Vietnam. This is the first report of wüstite found as an inclusion in gem sapphire. Its composition, Fe 1−x O with 0.04 < x < 0.12, should be stable at < 570 • C and can be transformed between α-Fe and Fe 3 O 4 [31]. As reported above, low power laser (0.5 mW) excitation can decrease alterations on the wüstite surface and then generate a Raman spectrum very similar to the pattern of magnetite (Fe 3 O 4 ) (Figure 6a). This indicates that the structure of wüstite is made up of non-equivalent sites, called non-stoichiometric wüstite (FeO), because the structure of non-equivalent sites of non-stoichiometric wüstite is similar to that of magnetite (Fe 3 O 4 ) [32]. However, this spectrum changes to a strong decrease in half-intensity peak positions of the magnetite structure (Fe 3 O 4 ) and clearly displays new peaks of hematite structure (Fe 2 O 3 ) plus a peak at 595 cm −1 attributed to wüstite, as noticed by Hanesch [23] (Figure 6b) at high laser power (5 mW). This effect supports that wüstite's surface with a magnetite structure might be transformed into a hematite structure by the activation of high laser power. The irreversible transformation of magnetite to hematite caused by a lowering of the temperature in the natural process is commonly known as matitization [33]: Besides, a strong increase in the peak around 667 cm -1 (magnetite) at high laser power was visible in the other wüstite inclusions (Figure 6c), which may be attributed to the Raman active vibrations of spinel groups [27], indicating more spinel components in the structure. They are comparable to the EPMA analyses that yielded non-stoichiometric wüstite (FeO) and wüstite-hercynite ( Table 1). The hercynite component (e.g., 23 percent in sample PT49) obtained in non-stoichiometric wüstites (Table 1) should affect the intensity of such a magnetite peak (Figure 6c).
Raman spectroscopic features of hercynite spinel inclusions clearly belong to the spinel group. Additionally, EPMA analyses yielded the composition of hercynite (Table 2 and Figure 9). This is the first report of wüstite found as an inclusion in gem sapphire. Its composition, Fe1−xO with 0.04 < x < 0.12, should be stable at < 570 °C and can be transformed between α-Fe and Fe3O4 [31]. As reported above, low power laser (0.5 mW) excitation can decrease alterations on the wüstite surface and then generate a Raman spectrum very similar to the pattern of magnetite (Fe3O4) (Figure 6a). This indicates that the structure of wüstite is made up of non-equivalent sites, called non-stoichiometric wüstite (FeO), because the structure of non-equivalent sites of non-stoichiometric wüstite is similar to that of magnetite (Fe3O4) [32]. However, this spectrum changes to a strong decrease in halfintensity peak positions of the magnetite structure (Fe3O4) and clearly displays new peaks of hematite structure (Fe2O3) plus a peak at 595 cm −1 attributed to wüstite, as noticed by Hanesch [23] (Figure 6b) at high laser power (5 mW). This effect supports that wüstite's surface with a magnetite structure might be transformed into a hematite structure by the activation of high laser power. The irreversible transformation of magnetite to hematite caused by a lowering of the temperature in the natural process is commonly known as matitization [33]: Besides, a strong increase in the peak around 667 cm -1 (magnetite) at high laser power was visible in the other wüstite inclusions (Figure 6c), which may be attributed to the Raman active vibrations of spinel groups [27], indicating more spinel components in the structure. They are comparable to the EPMA analyses that yielded non-stoichiometric wüstite (FeO) and wüstite-hercynite ( Table 1). The hercynite component (e.g., 23 percent in sample PT49) obtained in non-stoichiometric wüstites (Table 1) should affect the intensity of such a magnetite peak (Figure 6c).
Raman spectroscopic features of hercynite spinel inclusions clearly belong to the spinel group. Additionally, EPMA analyses yielded the composition of hercynite (Table 2 and Figure 9).  Raman spectroscopic features of ilmenite inclusions only exhibited the mixture patterns of ilmenite and hematite, called ilmenite solid solutions. In addition, the broadened peaks at 399, 600, and 1298 cm −1 of hematite ( Figure 8b) were caused by the oxidation of titanomagnetite. This oxidation seems to be influenced by ilmenite exsolution in the Ti-rich magnetite (titanomagnetite) which was reported previously in igneous intrusive from Bijigou, Pazhihua, and Xinjie in China [35]. Meanwhile, the EPMA analyses of ilmenite solid solutions clearly defined titanohematite and titanomagnetite (Table 3 and Figure 10). peaks at 399, 600, and 1298 cm −1 of hematite (Figure 8b) were caused by the oxidation of titanomagnetite. This oxidation seems to be influenced by ilmenite exsolution in the Tirich magnetite (titanomagnetite) which was reported previously in igneous intrusive from Bijigou, Pazhihua, and Xinjie in China [35]. Meanwhile, the EPMA analyses of ilmenite solid solutions clearly defined titanohematite and titanomagnetite (Table 3 and Figure 10). From the results, Raman spectroscopic features of ilmenite inclusions clearly exhibited ilmenite solid solutions (e.g., ilmenohematite and titanomagnetite), which are characteristic of rhombohedral solid solutions [25]; on the other hand, hercynite spinel and wüstite have cubic and octahedral shapes [30]. The EPMA shows attributes of wüstite and hercynite in octahedral iron oxide inclusions, whereas rhombohedral iron oxides favor ilmenite-hematite (titanohematite and titanomagnetite). The morphological characteristics of the iron oxide minerals in sapphire from Southern Vietnam match the results of Raman and EPMA analyses.

Crystallization Environment
Iron oxide inclusions in sapphire from Southern Vietnam are mostly characterized by euhedral crystals (rhombohedral and octahedral) which appear to be primary inclusions [36]. Therefore, their chemical and physical conditions can be used to reconstruct the crystallization environment of the host sapphire.
Non-stoichiometric wüstite appears to have formed under a strongly reducing condition [29,37]. Furthermore, the co-existing wüstite and hercynite inclusions in sapphire from Southern Vietnam ( Figure 5) may be involved by the hercynite's breakdown reaction process [38,39]: FeAl2O4 (hercynite) = Al2O3 (sapphire/ruby) + FeO (wüstite) Therefore, cogenetic wüstite and hercynite inclusions should be formed in a strongly reducing environment related to the hercynite breakdown reaction process. On the other hand, titanomagnetite and titanohematite inclusions should be generated by sub-solidus re-equilibration [1,40] which suggests an oxidizing environment. The compositions of titanomagnetite inclusions vary towards hematite components without ulvӧspinel components (Table 3) which should be re-equilibrated under the slow cooling oxidized sub-solidus environment [1,41]. This supports the complete oxidation reaction of spinelss by the From the results, Raman spectroscopic features of ilmenite inclusions clearly exhibited ilmenite solid solutions (e.g., ilmenohematite and titanomagnetite), which are characteristic of rhombohedral solid solutions [25]; on the other hand, hercynite spinel and wüstite have cubic and octahedral shapes [30]. The EPMA shows attributes of wüstite and hercynite in octahedral iron oxide inclusions, whereas rhombohedral iron oxides favor ilmenitehematite (titanohematite and titanomagnetite). The morphological characteristics of the iron oxide minerals in sapphire from Southern Vietnam match the results of Raman and EPMA analyses.

Crystallization Environment
Iron oxide inclusions in sapphire from Southern Vietnam are mostly characterized by euhedral crystals (rhombohedral and octahedral) which appear to be primary inclusions [36]. Therefore, their chemical and physical conditions can be used to reconstruct the crystallization environment of the host sapphire.
Non-stoichiometric wüstite appears to have formed under a strongly reducing condition [29,37]. Furthermore, the co-existing wüstite and hercynite inclusions in sapphire from Southern Vietnam ( Figure 5) may be involved by the hercynite's breakdown reaction process [38,39]: Therefore, cogenetic wüstite and hercynite inclusions should be formed in a strongly reducing environment related to the hercynite breakdown reaction process. On the other hand, titanomagnetite and titanohematite inclusions should be generated by sub-solidus re-equilibration [1,40] which suggests an oxidizing environment. The compositions of titanomagnetite inclusions vary towards hematite components without ulvöspinel components (Table 3) which should be re-equilibrated under the slow cooling oxidized subsolidus environment [1,41]. This supports the complete oxidation reaction of spinel ss by the re-equilibration: 4Fe 3 O 4 (in exsolution) + O 2 = 6Fe 2 O 3 (in ilmenite) to ilmenite exsolution in titanomagnetite.
This information indicates an environmental change during the formation of these oxide inclusions as well as their host sapphire. Initial crystallization should take place in a strongly reducing (indicated by wüstite occurrences) magma chamber prior to the slow cooling sub-solidus stage under a low oxidizing condition (based on phase transformation from titanomagnetite to titanohematite). This process may be related to oxygen fugacity and temperature changing towards the slow cooling process.
Inclusions and their host sapphire from Southern Vietnam were suggested to have formed in the lower crust [7,9] and in the lithospheric mantle [10]. In this study, wüstite with a cubic shape should be typified by crystallization at a temperature of about 570 • C [31]. This thermal process should be located in the continental crust (≤800 • C) [42]. Furthermore, wüstite appears to have occurred in the continental crust, as suggested by Seifert et al. [29], who reported wüstite in fluorapatite crystallized from S-type granite melts. Additionally, hercynite has been recognized in magmatic sapphires also formed in crust [43,44]. In addition, titanomagnetite inclusion with Il 24-30 components seems to be formed in plutonic rocks under a crustal environment, as suggested by Buddington et al. [45]. Therefore, coexisting wüstite, hercynite, titanohematite, and titanomagnetite clearly indicate that these oxide inclusions and their host sapphire should have crystallized directly from magma in crust instead of mantle.
Under the silica-saturated condition, iron oxide minerals (i.e., wüstite, hercynite, titanomagnetite, and titanohematite) are not stable; all Fe 2+ atoms preferentially enter silicate structures instead. Therefore, wüstite, hercynite, titanomagnetite, and titanohematite may only be formed in the silica-undersaturated environment. Titanomagnetite indicated a higher temperature (at about 600-700 • C) than the ilmenite-hematite miscibility gap [46]. Titanomagnetite should crystallize directly from magma prior to sub-solidus reequilibration of iron-titanium oxides during a slow cooling process within a low oxidizing environment. This late state re-equilibration led to a decrease in magnetite in spinel ss , with fO 2 increasing slightly.

Conclusions
Iron oxide mineral inclusions provide useful indicators of crystallization conditions of their host sapphire from Southern Vietnam. These iron oxide inclusions include wüstite, hercynite, titanomagnetite, and titanohematite series. Crystal morphology, Raman spectroscopy, and mineral chemical signatures of these iron oxides indicate silicaundersaturated magmatic origin in the lower crust. Wüstite might have crystallized from the hercynite's breakdown reaction (hercynite = sapphire + wüstite), whereas titanomagnetite and titanohematite series should have formed in sub-solidus re-equilibration within the slow cooling process. These results indicate an environmental change during the crystallization process of sapphire, wüstite, and hercynite in the reducing magma chamber prior to a slow cooling sub-solidus re-equilibrium of titanomagnetite to titanohematite with low oxidizing conditions. Supplementary Materials: The following are available online at https://www.mdpi.com/2075-163 X/11/3/241/s1, Table S1: Sapphire samples used in this study. Funding: The first author was supported for her PhD study by "the Scholarship Program for ASIAN Countries Chulalongkorn University". Institutional Review Board Statement: "Not applicable" for studies not involving humans or animals.