Gem-Quality Zircon Megacrysts from Placer Deposits in the Central Highlands , Vietnam — Potential Source and Links to Cenozoic Alkali Basalts

Gem-quality zircon megacrysts occur in placer deposits in the Central Highlands, Vietnam, and have euhedral to anhedral crystal shapes with dimensions of ~3 cm in length. These zircons have primary inclusions of calcite, olivine, and corundum. Secondary quartz, baddeleyite, hematite, and CO2 fluid inclusions were found in close proximity to cracks and tubular channels. LA-ICP-MS U-Pb ages of analyzed zircon samples yielded two age populations of ca. 1.0 Ma and ca. 6.5 Ma, that were consistent with the ages of alkali basalt eruptions in the Central Highlands at Buon Ma Thuot (5.80–1.67 Ma), Pleiku (4.30–0.80 Ma), and Xuan Loc (0.83–0.44 Ma). The zircon geochemical signatures and primary inclusions suggested a genesis from carbonatite-dominant melts as a result of partial melting of a metasomatized lithospheric mantle source, but not from the host alkali basalt. Chondrite-normalized rare earth element patterns showed a pronounced positive Ce, but negligible Eu anomalies. Detailed hyperspectral Dy3+ photoluminescence images of zircon megacrysts revealed resorption and re-growth processes.


Introduction
Zircon (ZrSiO 4 , tetragonal, I4 1 /amd) is an accessory mineral in most types of igneous and metamorphic rocks [1].Zircon megacrysts are often found in placer deposits derived from intraplate basaltic fields as xenocrysts or xenolith debris in alkali basaltic rocks [2].Alkali basalt fields in South-East Asia commonly contain mantle xenoliths that include garnet lherzolites, spinel lherzolites, and harzburgite, as well as mantle-and/or crust-derived megacrysts of pyroxene, olivine, plagioclase, garnet, zircon, and corundum [3].Therefore, zircon placers that are related to basaltic magmatism are often associated with megacrysts of other important gem materials (e.g., sapphire, garnet, and spinel), as reported from various localities, including Australia, Vietnam, Cambodia, Thailand, China, Myanmar, Sri Lanka, and Tanzania [4][5][6][7][8][9][10][11].Zircon material from the Central Highlands, Vietnam, shows brownish-red colors and blue, upon heat-treatment under reducing conditions [12].This color effect is similar to zircon found in placers in the Ratanakiri district, Cambodia [13].The Vietnamese zircon deposits are known since the late 1980s, and sporadic mining activities are conducted by local people at numerous small sites.These deposits are exploited by digging to a depth of approximately 1-2 m, then picking the gems by hand after washing the alluvial material.The exploration of placer deposits of gem-quality zircons and other gem megacrysts is of high economic importance in many South-East Asian countries that supply the global gem market [14].The source characteristics and forming processes of placer zircon play an important role for further exploration of this kind of gem deposits in geologically related regions.
The genesis of zircon megacrysts associated with alkali basalt has been addressed previously, but it is still vigorously debated [2,[15][16][17][18].Some authors agree that large zircon and corundum crystals represent xenocrysts that are transported by alkaline-basalt magmas and do not crystallize primarily from them (e.g., Hinton and Upton [2]), while others consider crystallization directly from basaltic magma [19,20].Various hypotheses have been reported for the origin of zircon megacrysts from various sites in different geological contexts: (1) crystallization from melts derived from a metasomatized upper mantle, as supported by O-isotope studies [15,17,18,21]; (2) formation within a late-stage fractional crystallization process of oceanic island basalt magma [16,19]; and (3) crystallization from a primitive alkaline mafic magma, which later evolved to a less alkaline host magma [20].According to Cong et al. [10], zircon megacrysts from Cenozoic basalts in northeastern Cambodia crystallized in the mantle during metasomatic events caused by phosphate-rich fluids and/or silicate melts enriched in zirconium.Piilonen et al. [22] recently reported findings of megacrystic zircons as single crystals enclosed in the alkali basalt from Ratanakiri, Cambodia, and suggested a single stage growth in a carbonate-influenced environment.
In the Central Highlands of Vietnam, gem quality zircon is commonly accompanied by corundum (sapphire) in alluvial deposits that are considered eroded from Cenozoic alkali basalt fields in close proximity (e.g., Garnier et al. [17]).In this study, zircon megacrysts from various alluvial deposits in the Central Highlands were systematically studied to investigate their origin and potential links to Cenozoic basalt eruptions.In addition to trace element data and U-Pb geochronological data, photoluminescence (PL) imaging methods were applied to reveal internal growth and potential secondary alteration textures.In this paper, we used mineral abbreviations as given by Whitney and Evans [23].

Geological Setting
South-East Asia was formed by an amalgamation of several crustal blocks, including South China, Indochina, Siumasu, Inthanon, the West Burma block, and the Trans Vietnam Orogenic Belt (TVOB) (Figure 1A) [24].The TVOB was first proposed by Osanai et al. [24] as a zone of Permo-Triassic metamorphic rocks in Vietnam, which were formed by continent-continent collision between South China and Indochina blocks.This orogenic belt is characterized by numerous shear zones with strong deformation, such as the Red River shear zone, Song Ma suture zone, Tam Ky-Phuoc Son shear zone, and the Dak To Kan shear zone (Kontum Massif).The northern extension of these shear zones reaches the Yuan Nan Province through Ailaoshan in China.The Central Highlands lies entirely within the Indochina block and are located in the southern part of the TVOB, partly extending into eastern Vietnam and Laos (e.g., Hutchison [25]) (Figure 1A).
Miocene-Pliocene alkaline and subalkaline basalts are exposed over a vast region of Thailand, Laos, and Vietnam [26].In southern and central Vietnam, this province stretches across an area of approximately 23,000 km 2 with a thickness up to several hundred meters [3] (Figure 1A).The basalt plateau is accompanied by pull-apart structures composed of short extensional rifts bounded by strike-slip faults [27].At least two eruptive episodes ("early" and "late") have been reported by Hoang et al. [28].Tholeiitic, and rarely, alkali basaltic flows represent the "early" episode at Dalat (17.60-7.90Ma), whereas in the "late" episodes, olivine tholeiite, alkali basalt, basanite, and (rarely) nephelinite erupted at Phuoc Long (<8.00-3.40Ma), Buon Ma Thuot (5.80-1.67Ma), Pleiku (4.30-0.80Ma), Xuan Loc (0.83-0.44 Ma), and the Re Island centers (0.80-0 Ma) [27] (Figure 1B).Hoang et al. [28] proposed that tholeiitic basalts are the most common basalt type in the region and build up much more volumetric mass in comparison to the alkali basalts and the rarely occurring nephelinites which erupted from small volcanoes.

Sample Description and Methods
In this study, representative samples of zircon megacrysts collected from several alluvial deposits in the Central Highlands, Vietnam, were investigated in more detail.A representative zircon sample (Rata) from the Ratanakiri district, Cambodia was included for comparison.Large, gemquality zircon specimens have dimensions up to several centimeters (Figure 2).Their color varies from colorless, orange, brownish-orange to dark brown and dark red (Figure 2).Some crystals have anhedral to subhedral shapes with rounded termination, while most grains are euhedral crystals with a typical combination of bipyramid and tetragonal prism (Figure 2B).Occasionally, internal colored zones that comprise an oscillatory change of orange and brown bands are visible to the naked eye (Figure 2B, sample C16).
Selected samples were cut along their long axis for petrographic analyses.Inclusions in zircon samples were observed and identified using an optical microscope attached to a confocal laser Raman

Sample Description and Methods
In this study, representative samples of zircon megacrysts collected from several alluvial deposits in the Central Highlands, Vietnam, were investigated in more detail.A representative zircon sample (Rata) from the Ratanakiri district, Cambodia was included for comparison.Large, gem-quality zircon specimens have dimensions up to several centimeters (Figure 2).Their color varies from colorless, orange, brownish-orange to dark brown and dark red (Figure 2).Some crystals have anhedral to subhedral shapes with rounded termination, while most grains are euhedral crystals with a typical combination of bipyramid and tetragonal prism (Figure 2B).Occasionally, internal colored zones that comprise an oscillatory change of orange and brown bands are visible to the naked eye (Figure 2B, sample C16).
Selected samples were cut along their long axis for petrographic analyses.Inclusions in zircon samples were observed and identified using an optical microscope attached to a confocal laser Raman spectrometer system (JASCO NSR-3100, JASCO (USA)) at Kyushu University.A 100× objective (with a numerical aperture NA = 0.90) and a green, continuous 532 nm frequency-doubled YAG:Nd laser (with an energy output of approx.10 mW at the sample surface) was used to perform the Raman spot analyses.Note that potential bias of Raman spectra from the photoluminescence (PL) of Er 3+ in zircon has been reported when using this laser wavelength [29].A very minor PL contribution of Er 3+ was identified in most recorded Raman spectra, but did not hamper the identification of typical Raman bands as obtained from various inclusion phases.Inclusions that were accessible via the polished sample surface were further identified qualitatively using a JEOL JSM-5310S-JED2140 scanning electron microprobe (SEM, JEOL, Ltd., Tokyo, Japan) equipped with an energy-dispersive X-Ray (EDX) detector (Li-doped Si semiconductor) at Kyushu University.
The internal texture of the prepared zircon crystals was investigated using a back-scattered electron (BSE) detector system and a GATAN MiniCL panchromatic cathodoluminescence (CL) detector attached to the above-mentioned SEM system at Kyushu University.In addition, we applied laser-induced photoluminescence spectroscopy using a HORIBA LabRam HR800 Evolution spectrometer equipped with a Peltier-cooled Si detector and a grating of 600 lines per millimeter (Horiba, Ltd., Kyoto, Japan).Confocal PL spot measurements were performed using an Olympus BX80 microscope (manufactured by Olympus Corporation, Tokyo, Japan) with a solid-state continuous laser, operating at 473 nm, to excite the most prominent REE 3+ photoluminescence emissions in zircon (e.g., Dy 3+ ; Lenz et al. [29]).Hyperspectral images of polished, large zircon single-crystals were produced by measuring a multitude of spots in a point-by-point raster with a step-width of 10-20 µm using a software-controlled Maerzhauser mechanical x-y table.A software-based, automated data treatment procedure for each single spectrum was applied, and the plotted spectral parameters of interest were color-coded.With respect to the large size of investigated zircon crystals, and the comparatively large mapping step-widths, we used a 50× objective (NA = 0.5) and a confocal hole of 200 µm to adjust the lateral spatial resolution to be approximately 8-10 µm.However, we note that a diffraction-limited maximum spatial resolution of ~1 µm (planar) and 2-3 µm in depth may be achieved using the combination of a small confocal hole and a 100× objective.
Zircon U-Pb dating of samples C05, C07, and C16 was carried out by applying laser-ablation inductively coupled plasma mass-spectrometry (LA-ICP-MS) using an Agilent 7500cx quadrupole ICP-MS with a New Wave Research UP-213 YAG:Nd laser at the Kyushu University, that produced a laser-ablation spot of 100 µm.Details of the analytical procedure are presented by Adachi et al. [30].The zircon standards Temora (417 Ma; [31]) and FC-1 (1099 Ma; [32]) were used for calibration and accuracy checks, respectively.The NIST SRM-611 glass standard was applied to determine the Th/U ratios.Raw data were treated using GLITTER software (Version 4.4.2,Glitter sold through Access Macquarie Ltd., Sydney, Australia) [33] and plotted in concordia diagrams using the Isoplot/Ex 3.7 software (Berkeley Geochronology Center, Berkeley, CA, USA) [34].Trace element concentrations were obtained from similar zones in close proximity to spots used for U-Pb dating.Details of the analytical procedure are given in Nakano et al. [35] using the theoretical Si concentration in zircon as an internal standard.
Dating of zircon samples DL1, DL2, and DL3 was carried out using an Agilent 7700 quadrupole ICP-MS instrument attached to a Photon Machines Excimer 193 nm laser system at GEMOC, Macquarie University, using a beam diameter of ca.50 µm with a 5 Hz repetition rate, and energy of around 0.06 J/cm 2 to 8 J/cm 2 .Ablation was carried out in He gas to improve the sample transport efficiency, and to provide stable signals with reproducible Pb/U fractionation.Sample analyses were run along with analyses of the GEMOC GJ-1 zircon standard [36].This standard is slightly discordant, and has a TIMS 207 Pb/ 206 Pb age of 608.5 Ma [37].The other well-characterized zircon standard 91500 and Mud Tank were analysed within the run as an independent control on reproducibility and instrument stability.Individual time-resolved data analysis provides isotopically homogeneous segments of the signal to be selected for integration.We corrected the integrated ratios for ablation related fractionation and instrumental mass bias by the calibration of each selected time segment against the identical time segments for the standard zircon analyses.Furthermore, we employed the common-Pb correction procedure described by Andersen [38].The analyses presented here were corrected assuming recent lead-loss with a common-lead composition corresponding to the present-day average orogenic lead, as given by the second-stage growth curve of Stacey and Kramers [39] for 238 U/ 204 Pb = 9.74.No correction was applied to the analyses that were concordant within a 2σ analytical error in 206 Pb/ 238 U and 207 Pb/ 235 U, or which had less than 0.2% common lead.Trace element concentrations of samples DL1, DL2, and DL3 were obtained using the same LA-ICP-MS system.Calibration of the relative element sensitivities were done using the NIST-610 standard glass as external calibration.Zircon BCR-2g and GJ-1 reference material were analyzed along the measurement runs as an independent control on reproducibility and instrument stability.Theoretical Zr content of the zircon was used for internal calibration of unknown zircon samples.The precision and accuracy of the NIST-610 analyses were 1-2% for REE, Y, Nb, Hf, Ta, Th, and U at the ppm concentration level, and at 5-10% for Ti (further details see in Belousova et al. [16]).
corrected assuming recent lead-loss with a common-lead composition corresponding to the presentday average orogenic lead, as given by the second-stage growth curve of Stacey and Kramers [39] for 238 U/ 204 Pb = 9.74.No correction was applied to the analyses that were concordant within a 2σ analytical error in 206 Pb/ 238 U and 207 Pb/ 235 U, or which had less than 0.2% common lead.Trace element concentrations of samples DL1, DL2, and DL3 were obtained using the same LA-ICP-MS system.Calibration of the relative element sensitivities were done using the NIST-610 standard glass as external calibration.Zircon BCR-2g and GJ-1 reference material were analyzed along the measurement runs as an independent control on reproducibility and instrument stability.Theoretical Zr content of the zircon was used for internal calibration of unknown zircon samples.The precision and accuracy of the NIST-610 analyses were 1-2% for REE, Y, Nb, Hf, Ta, Th, and U at the ppm concentration level, and at 5-10% for Ti (further details see in Belousova et al. [16]).

Inclusion Features
Zircon crystals investigated in this study contained various fluid and mineral inclusions throughout the grains, especially along the cracks (Figure 3A-G).Tubular channels, bifurcating cavities, and vesicles in association with fluid inclusions were observed along the 2-dimensional planes (Figure 3A-G).Typical fluid inclusions along the fractures or fissures consisted of two phases, liquid (H 2 O) with CO 2 or O 2 bubbles.Mineral inclusions were identified as calcite, hematite, corundum, olivine, baddeleyite, quartz, and feldspar.We identified two different paragenetic types of mineral inclusions: (1) the inclusions that were distributed throughout the grains without fluid inclusions in close proximity; and (2) others that were distributed along the cracks or were associated with tubular channels, bifurcating cavities, and vesicles with fluid inclusions nearby.Calcite inclusions were found in both areas: ovoid and droplet shaped ones sat inside the tubular channels or were associated with fluid inclusions (Figure 3B); whilst square or rectangular shaped calcite inclusions were accumulated along the chains or in groups that did not occur next to the fluid inclusions or within the tubular channels (Figure 3I).We rarely found corundum (Figure 3H) and olivine inclusions.They were of an euhedral and subhedral shape, and appeared in areas where no cracks, fluid inclusions, or tubular channels were present.Quartz, baddeleyite, hematite, and feldspar were found with irregular or ovoid shapes with rounded termination.They were located along the fissures, within tubular channels, cavities or associated with fluid inclusions (Figure 3C-G).

Inclusion Features
Zircon crystals investigated in this study contained various fluid and mineral inclusions throughout the grains, especially along the cracks (Figure 3A-G).Tubular channels, bifurcating cavities, and vesicles in association with fluid inclusions were observed along the 2-dimensional planes (Figure 3A-G).Typical fluid inclusions along the fractures or fissures consisted of two phases, liquid (H2O) with CO2 or O2 bubbles.Mineral inclusions were identified as calcite, hematite, corundum, olivine, baddeleyite, quartz, and feldspar.We identified two different paragenetic types of mineral inclusions: (1) the inclusions that were distributed throughout the grains without fluid inclusions in close proximity; and (2) others that were distributed along the cracks or were associated with tubular channels, bifurcating cavities, and vesicles with fluid inclusions nearby.Calcite inclusions were found in both areas: ovoid and droplet shaped ones sat inside the tubular channels or were associated with fluid inclusions (Figure 3B); whilst square or rectangular shaped calcite inclusions were accumulated along the chains or in groups that did not occur next to the fluid inclusions or within the tubular channels (Figure 3I).We rarely found corundum (Figure 3H) and olivine inclusions.They were of an euhedral and subhedral shape, and appeared in areas where no cracks, fluid inclusions, or tubular channels were present.Quartz, baddeleyite, hematite, and feldspar were found with irregular or ovoid shapes with rounded termination.They were located along the fissures, within tubular channels, cavities or associated with fluid inclusions (Figure 3C-G).

Internal Texture of Zircon Megacrysts
Backscattered-electron images obtained from the analysis of polished sections revealed no internal textural features although high contrast levels were applied (not shown).The latter indicated a homogeneous distribution and/or low concentration levels of trace elements.In contrast, cathodoluminescence images shown in Figure 4, revealed complex internal textures comprised of oscillatory zoning and sector zoning, although overall CL intensities were comparably low, except for the colorless zircon (Zrn-C05, see Figure 4).Anhedral grains were found to be split or broken into sections of former even larger grains (e.g., F03, F04), whereas the euhedral grains were characterized by multiple CL growth bands that retraced multiple distinct growth stages within a single crystal (e.g., C16 in Figure 4).To obtain minor differences in the CL zonation, high intensity contrast-levels were applied.This caused white stripes to be present in the CL images of some samples that did not correlate with the crystal's internal zonation (Sample C01, C07, and F01).They were most prominent along cracks and break-outs of fine material that produced mechanically-induced defects on the sample surface during polishing.Those structural defect centers mostly caused broad and intense CL signals in the UV spectral region, which were detected efficiently using panchromatic CL detectors.

Internal Texture of Zircon Megacrysts
Backscattered-electron images obtained from the analysis of polished sections revealed no internal textural features although high contrast levels were applied (not shown).The latter indicated a homogeneous distribution and/or low concentration levels of trace elements.In contrast, cathodoluminescence images shown in Figure 4, revealed complex internal textures comprised of oscillatory zoning and sector zoning, although overall CL intensities were comparably low, except for the colorless zircon (Zrn-C05, see Figure 4).Anhedral grains were found to be split or broken into sections of former even larger grains (e.g., F03, F04), whereas the euhedral grains were characterized by multiple CL growth bands that retraced multiple distinct growth stages within a single crystal (e.g., C16 in Figure 4).To obtain minor differences in the CL zonation, high intensity contrast-levels were applied.This caused white stripes to be present in the CL images of some samples that did not correlate with the crystal's internal zonation (Sample C01, C07, and F01).They were most prominent along cracks and break-outs of fine material that produced mechanically-induced defects on the sample surface during polishing.Those structural defect centers mostly caused broad and intense CL signals in the UV spectral region, which were detected efficiently using panchromatic CL detectors.We applied laser-induced PL hyperspectral mapping to the large euhedral crystal grains with the aim of visualizing the internal luminescence distribution patterns as specifically caused by the emissions of REEs (Figure 5).The advantage of applying the latter technique is that the PL of specific REE species may be excited effectively using the appropriate laser sources in the visible spectral range without being obscured by intensive broad-band luminescence features as observed by CL spectroscopy and imaging [29].Figure 5 shows hyperspectral images of the zircon samples from the Central Highlands (C16, DL2) and from Ratanakiri, Cambodia (Rata).The integrated intensity of the  We applied laser-induced PL hyperspectral mapping to the large euhedral crystal grains with the aim of visualizing the internal luminescence distribution patterns as specifically caused by the emissions of REEs (Figure 5).The advantage of applying the latter technique is that the PL of specific REE species may be excited effectively using the appropriate laser sources in the visible spectral range without being obscured by intensive broad-band luminescence features as observed by CL spectroscopy and imaging [29].Figure 5 shows hyperspectral images of the zircon samples from the Central Highlands (C16, DL2) and from Ratanakiri, Cambodia (Rata).The integrated intensity of the most prominent PL emission of Dy 3+ ( 4 F 9/2 → 6 H 13/2 ) was used as a plotted spectral parameter (see Figure 6A; and Figure 1 in Lenz et al. [29]).Note that the Dy sub-level bands of Vietnamese and Cambodian (Rata) zircon samples were characterized by exceptionally narrow band-widths (sublevel at 581 nm had a band-width of 11.0 cm −1 ; see Figure 6A), which indicated a high degree of crystallinity and the absence of structural radiation-damage that might be caused by decay of radioactive U and Th.This was in accordance to findings by Lenz and Nasdala [40] and Zeug et al. [13] who reported very low PL and Raman band widths for zircon from Ratanakiri, Cambodia (Raman ν 3 [SiO 4 ] ~1.8 cm −1 ).This has been interpreted to be the result of very low α-doses due to low U and Th concentrations in combination with a very young age (ca. 1 Ma).
In this study, we found a systematic linear correlation of the integrated PL intensity of Dy with its concentration as obtained from multiple LA-ICP-MS measurement spots (Figure 6B and Table S1).Furthermore, the latter laser-ablation reference spots of Dy concentrations and their correlation with PL intensity were used as an external calibration for the PL hyperspectral images to infer Dy concentrations from the PL signal in regions of unknown trace-element chemistry (see color-coded scales in Figure 5).Especially in the euhedral crystals (e.g., C16 and DL3), we identified several growth stages that were characterized by abrupt changes in concentrations of trace Dy (REE).Typically, the crystal's cores were found to be enriched in Dy, followed by a decreasing oscillatory REE substitutional budget (see inset of sample C16 in Figure 5).PL patterns further revealed episodical resorption of pre-existent crystals and fast re-growth.Growth zones that were characterized by low REE concentrations cross-cut former zonation patterns and occurred along with large holes and cavities (see arrows in Figure 5, samples C16 and DL3).
Minerals 2019, 9, 89 8 of 18 most prominent PL emission of Dy 3+ ( 4 F9/2 → 6 H13/2) was used as a plotted spectral parameter (see Figure 6A; and Figure 1 in Lenz et al. [29]).Note that the Dy sub-level bands of Vietnamese and Cambodian (Rata) zircon samples were characterized by exceptionally narrow band-widths (sublevel at 581 nm had a band-width of 11.0 cm −1 ; see Figure 6A), which indicated a high degree of crystallinity and the absence of structural radiation-damage that might be caused by decay of radioactive U and Th.This was in accordance to findings by Lenz and Nasdala [40] and Zeug et al. [13] who reported very low PL and Raman band widths for zircon from Ratanakiri, Cambodia (Raman ν3 [SiO4] ~1.8 cm −1 ).This has been interpreted to be the result of very low α-doses due to low U and Th concentrations in combination with a very young age (ca. 1 Ma).
In this study, we found a systematic linear correlation of the integrated PL intensity of Dy with its concentration as obtained from multiple LA-ICP-MS measurement spots (Figure 6B and Table S1).Furthermore, the latter laser-ablation reference spots of Dy concentrations and their correlation with PL intensity were used as an external calibration for the PL hyperspectral images to infer Dy concentrations from the PL signal in regions of unknown trace-element chemistry (see color-coded scales in Figure 5).Especially in the euhedral crystals (e.g., C16 and DL3), we identified several growth stages that were characterized by abrupt changes in concentrations of trace Dy (REE).Typically, the crystal's cores were found to be enriched in Dy, followed by a decreasing oscillatory REE substitutional budget (see inset of sample C16 in Figure 5).PL patterns further revealed episodical resorption of pre-existent crystals and fast re-growth.Growth zones that were characterized by low REE concentrations cross-cut former zonation patterns and occurred along with large holes and cavities (see arrows in Figure 5, samples C16 and DL3). ) and from Ratanakiri, Cambodia (Rata).The integrated intensity of the Dy 3+ ( 4 F9/2 → 6 H13/2) emission (compare Figure 6) is the plotted spectral parameter (grey-scale).Traceelement concentrations of Dy from multiple LA-ICP-MS spots on the grains, were used for external calibration to correlate the Dy concentrations with its PL response.Concentrations given color-coded were back-calculated using the latter correlation (see Figure 6) to extrapolate the Dy distribution in regions of unknown trace-element chemistry.Besides oscillatory and sector zoning, several growth zones may be distinguished that were interpreted to result from resorption and re-growth events, as indicated by the cross-cut zonation patterns, abruptly changing (decreasing) Dy concentrations, and the appearance of large holes and cavities (see arrows). ) and from Ratanakiri, Cambodia (Rata).The integrated intensity of the Dy 3+ ( 4 F 9/2 → 6 H 13/2 ) emission (compare Figure 6) is the plotted spectral parameter (grey-scale).Trace-element concentrations of Dy from multiple LA-ICP-MS spots on the grains, were used for external calibration to correlate the Dy concentrations with its PL response.Concentrations given color-coded were back-calculated using the latter correlation (see Figure 6) to extrapolate the Dy distribution in regions of unknown trace-element chemistry.Besides oscillatory and sector zoning, several growth zones may be distinguished that were interpreted to result from resorption and re-growth events, as indicated by the cross-cut zonation patterns, abruptly changing (decreasing) Dy concentrations, and the appearance of large holes and cavities (see arrows).S1).Linear regression from the data pairs PL vs. Dy conc.(R² = 0.95) was used to calibrate hyperspectral maps presented in Figure 5.

Geochemical Characteristics
Results from trace-element analyses of Vietnamese zircon samples, in addition to a Cambodian sample from Ratanakiri (Rata) are summarized in Table S1.Representative samples that comprised a colorless (C05), a dark-red (C07), an orange (C16), and three heavily zoned samples (DL1, DL2, DL3), were selected for detailed analyses.All zircon megacrysts had low total Y + REE concentrations typically within a 50-550 ppm range.Somewhat higher concentrations were detected at measurement spots placed in the core region of the heavily zoned samples DL2 and DL3 that were up to 2300 ppm Y + REE (see Table S1).The chondrite-normalized REE concentrations revealed a steep slope from La to Lu (Figure 7).Note that the REE geochemical characteristics of samples from the Central Highlands closely coincided with those of the samples from Ratanakiri (see samples Rata and data from Cong et al. [10]).All samples were characterized by a pronounced positive Ce-anomaly and a lack of Euanomaly, which is typical of magmatic zircon crystallized from more evolved magma sources with co-crystallization of feldspars that substitute Eu 2+ on Ca-sites.The average Eu/Eu* ratio, that reflects the deviation of the measured Eu concentration from a theoretical Eu* concentration-extrapolated from chondrite-normalized Sm and Gd abundance (Eu* = EuN / [SmN + GdN] 0.5 )-was found to be close to 1 (Figure 8A).The obtained REE compositional characteristics excluded crustal derived rocks or granitoids as a potential source of the zircon samples studied, and were more typical of the zircon of kimberlitic and/or carbonatitic origin (Figure 8A).Similar to REEs, the overall U and Th concentrations were found to be comparably low, with typical variability at 10-60 ppm and 3-30 ppm, respectively.Core regions of samples DL2 and 3 were found to be more enriched in U and Th with concentrations up to 500 ppm.Despite the pronounced variation in U and Th concentrations among the different growth zones within the samples, their U/Th ratios showed a clear correlation (Figure 8B).Likewise, the Y and U data pairs in a Y vs. U discrimination diagram scattered appreciably amongst the different growth zones within and among the samples, but coincided well with the compositional fields that were reported for basic, kimberlitic, and carbonatitic rocks (Figure 8C; Belousova et al. [16]).A source rock discrimination based on the Hf composition as proposed by Shnukov et al. [41], is presented in Figure 8D.The samples' Hf composition was found to be very low and shared compositional similarities with zircon from carbonatitic rocks, which had the lowest reported Hf concentrations (<0.7 wt %).S1).Linear regression from the data pairs PL vs. Dy conc.(R 2 = 0.95) was used to calibrate hyperspectral maps presented in Figure 5.

Geochemical Characteristics
Results from trace-element analyses of Vietnamese zircon samples, in addition to a Cambodian sample from Ratanakiri (Rata) are summarized in Table S1.Representative samples that comprised a colorless (C05), a dark-red (C07), an orange (C16), and three heavily zoned samples (DL1, DL2, DL3), were selected for detailed analyses.All zircon megacrysts had low total Y + REE concentrations typically within a 50-550 ppm range.Somewhat higher concentrations were detected at measurement spots placed in the core region of the heavily zoned samples DL2 and DL3 that were up to 2300 ppm Y + REE (see Table S1).The chondrite-normalized REE concentrations revealed a steep slope from La to Lu (Figure 7).Note that the REE geochemical characteristics of samples from the Central Highlands closely coincided with those of the samples from Ratanakiri (see samples Rata and data from Cong et al. [10]).All samples were characterized by a pronounced positive Ce-anomaly and a lack of Eu-anomaly, which is typical of magmatic zircon crystallized from more evolved magma sources with co-crystallization of feldspars that substitute Eu 2+ on Ca-sites.The average Eu/Eu* ratio, that reflects the deviation of the measured Eu concentration from a theoretical Eu* concentration-extrapolated from chondrite-normalized Sm and Gd abundance (Eu* = Eu N / [Sm N + Gd N ] 0.5 )-was found to be close to 1 (Figure 8A).The obtained REE compositional characteristics excluded crustal derived rocks or granitoids as a potential source of the zircon samples studied, and were more typical of the zircon of kimberlitic and/or carbonatitic origin (Figure 8A).Similar to REEs, the overall U and Th concentrations were found to be comparably low, with typical variability at 10-60 ppm and 3-30 ppm, respectively.Core regions of samples DL2 and 3 were found to be more enriched in U and Th with concentrations up to 500 ppm.Despite the pronounced variation in U and Th concentrations among the different growth zones within the samples, their U/Th ratios showed a clear correlation (Figure 8B).Likewise, the Y and U data pairs in a Y vs. U discrimination diagram scattered appreciably amongst the different growth zones within and among the samples, but coincided well with the compositional fields that were reported for basic, kimberlitic, and carbonatitic rocks (Figure 8C; Belousova et al. [16]).A source rock discrimination based on the Hf composition as proposed by Shnukov et al. [41], is presented in Figure 8D.The samples' Hf composition was found to be very low and shared compositional similarities with zircon from carbonatitic rocks, which had the lowest reported Hf concentrations (<0.7 wt %).[42].All samples showed a steep slope rise from La to Lu, towards heavy REEs.A positive pronounced Ce-anomaly, but no Eu-anomaly was detected.

Geochronology
Results of the measured U-Th-Pb isotopic ratios and calculated ages of samples DL1, DL2, C05, C07, C16 are summarized in Tables S2 and S3 (see supplementary materials) and Figure 9. Vietnamese zircon samples from the Central Highlands were found to have Cenozoic, fairly concordant (~50% of discordance in average), and U−Pb ages (Figure 9A).Detailed comparison of the more robust 206 Pb/ 238 U ages indicated two distinct events of zircon megacrysts formation.Samples DL1, DL2, and C16 had a U-Pb age of around ca. 6.5 Ma, whereas samples C05 and C07 were even younger with an age of ca. 1 Ma.Note that the LA-ICP-MS measurement spots of individual single crystals were from different growth zones among the entire grains, which were clearly discernable in CL imaging, as well as PL hyperspectral mapping, and they had significant differences in traceelement concentrations (results above).We, however, found no systematic difference in U-Pb ages between the growth zones of individual samples within the standard errors.

Discussion
Zircon megacrysts from the Central Highlands that were investigated in this study had various colors, ranging from colorless, orange, brownish-orange, and dark brown to dark red.Some were characterized by colored zonation visible to the naked eye.Most of the selected grains showed a euhedral crystal shape with typical combinations of bipyramid and tetragonal prisms.Other grains had a subhedral to anhedral shape with rounded termination (Figure 2B).The latter might result from magmatic resorption on crystal faces and/or erosion during weathering and transportation into the placer.The internal texture of single crystals revealed by CL imaging and PL hyperspectral mapping was characterized by wide to narrow oscillatory and/or sector zoning (Figures 4 and 5).Typical crystal shape and textural features were conclusive hints that the megacrysts had a magmatic origin [1].In heavily zoned samples (DL1−3, C16), dark-brownish colored zones correlated with the elevated concentrations of REE, U, and Th.On the other hand, a colorless zircon sample (C05, Figure 2B) was found to be more enriched in REEs, Th, and U, than the non-transparent dark-red sample C07 (Figures 2B, 7 and Table S2).Hence, any rigorous interpretation based on color information should be avoided.For example, we found no corroboration for the interpretation that dark-brownish coloration of zircon was related to the presence of structural radiation damage.Multiple zircon samples obtained from placer deposits in Vietnam, including samples from Ratanakiri, Cambodia, had various colors ranging from transparent to non-transparent dark-brown, but all were characterized by very low Raman and PL spectral band-widths that indicated no radiation damage to be present [13].

Geochronology
Results of the measured U-Th-Pb isotopic ratios and calculated ages of samples DL1, DL2, C05, C07, C16 are summarized in Tables S2 and S3 (see Supplementary Materials) and Figure 9. Vietnamese zircon samples from the Central Highlands were found to have Cenozoic, fairly concordant (~50% of discordance in average), and U-Pb ages (Figure 9A).Detailed comparison of the more robust 206 Pb/ 238 U ages indicated two distinct events of zircon megacrysts formation.Samples DL1, DL2, and C16 had a U-Pb age of around ca. 6.5 Ma, whereas samples C05 and C07 were even younger with an age of ca. 1 Ma.Note that the LA-ICP-MS measurement spots of individual single crystals were from different growth zones among the entire grains, which were clearly discernable in CL imaging, as well as PL hyperspectral mapping, and they had significant differences in trace-element concentrations (results above).We, however, found no systematic difference in U-Pb ages between the growth zones of individual samples within the standard errors.

Discussion
Zircon megacrysts from the Central Highlands that were investigated in this study had various colors, ranging from colorless, orange, brownish-orange, and dark brown to dark red.Some were characterized by colored zonation visible to the naked eye.Most of the selected grains showed a euhedral crystal shape with typical combinations of bipyramid and tetragonal prisms.Other grains had a subhedral to anhedral shape with rounded termination (Figure 2B).The latter might result from magmatic resorption on crystal faces and/or erosion during weathering and transportation into the placer.The internal texture of single crystals revealed by CL imaging and PL hyperspectral mapping was characterized by wide to narrow oscillatory and/or sector zoning (Figures 4 and 5).Typical crystal shape and textural features were conclusive hints that the megacrysts had a magmatic origin [1].In heavily zoned samples (DL1−3, C16), dark-brownish colored zones correlated with the elevated concentrations of REE, U, and Th.On the other hand, a colorless zircon sample (C05, Figure 2B) was found to be more enriched in REEs, Th, and U, than the non-transparent dark-red sample C07 (Figure 2B, Figure 7 and Table S2).Hence, any rigorous interpretation based on color information should be avoided.For example, we found no corroboration for the interpretation that dark-brownish coloration of zircon was related to the presence of structural radiation damage.Multiple zircon samples obtained from placer deposits in Vietnam, including samples from Ratanakiri, Cambodia, had various colors ranging from transparent to non-transparent dark-brown, but all were characterized by very low Raman and PL spectral band-widths that indicated no radiation damage to be present [13].S2 and  S3).We used PL hyperspectral imaging of the Dy 3+ emission as an effective tool to visualize the distribution of REEs across the zircon samples (Figure 5).The latter technique was more reliable with respect to specifically exciting and detecting the emissions of REEs, because the panchromatic CL imaging is often strongly affected by other luminescence emissions that are caused by defects that are not coupled to the substitution of REEs (see mechanically induced bright CL striations caused by polishing; Figure 4).We further found that PL integrated intensities of Dy 3+ used for PL hyperspectral imaging, correlated well with the Dy concentrations determined using spot LA-ICP-MS analyses (Figure 6).We, therefore, used PL hyperspectral images to quantitatively visualize Dy concentrations in regions of unknown chemical compositions as extrapolated from Dy 3+ integrated emission intensities.Note, however, that the correlation of Dy concentrations with their PL response may be strongly hampered in other zircon samples that show a much higher accumulation of structural radiation damage, e.g., Lenz and Nasdala [40] reported that the presence of radiation damage resulted in substantial quenching of the PL intensities.
In this study, we identified multiple growth stages that included the resorption and re-growth of large volumes of the zircon megacrysts (Figure 5, samples C16 and DL3).Outer growth zones irregularly cut the former zones repeatedly and were characterized by lower REE concentrations, which indicated that the growth condition of megacrysts may have changed spontaneously to a more REE undersaturated growth environment that resulted in re-growth of the zones with depleted REE concentrations.Although absolute REE concentrations changed considerably across the different growth zones, their general chondrite-normalized concentration patterns were qualitatively similar.S2 and  S3).We used PL hyperspectral imaging of the Dy 3+ emission as an effective tool to visualize the distribution of REEs across the zircon samples (Figure 5).The latter technique was more reliable with respect to specifically exciting and detecting the emissions of REEs, because the panchromatic CL imaging is often strongly affected by other luminescence emissions that are caused by defects that are not coupled to the substitution of REEs (see mechanically induced bright CL striations caused by polishing; Figure 4).We further found that PL integrated intensities of Dy 3+ used for PL hyperspectral imaging, correlated well with the Dy concentrations determined using spot LA-ICP-MS analyses (Figure 6).We, therefore, used PL hyperspectral images to quantitatively visualize Dy concentrations in regions of unknown chemical compositions as extrapolated from Dy 3+ integrated emission intensities.Note, however, that the correlation of Dy concentrations with their PL response may be strongly hampered in other zircon samples that show a much higher accumulation of structural radiation damage, e.g., Lenz and Nasdala [40] reported that the presence of radiation damage resulted in substantial quenching of the PL intensities.
In this study, we identified multiple growth stages that included the resorption and re-growth of large volumes of the zircon megacrysts (Figure 5, samples C16 and DL3).Outer growth zones irregularly cut the former zones repeatedly and were characterized by lower REE concentrations, which indicated that the growth condition of megacrysts may have changed spontaneously to a more REE undersaturated growth environment that resulted in re-growth of the zones with depleted REE concentrations.Although absolute REE concentrations changed considerably across the different growth zones, their general chondrite-normalized concentration patterns were qualitatively similar.
Chondrite-normalized REE patterns in the core and overgrowth regions had a steep slope from light REE (LREE) to heavy REE (HREE), with pronounced positive Ce-, but with a lack of an Eu-anomaly (Figure 7).The positive Ce anomaly is typically caused by the significant difference between Ce 4+ and LREE 3+ to fit into the zircon's structure, as Ce 4+ has the same charge and similar ionic radius compared to Zr 4+ .The Eu anomaly of zircon in REE patterns is generally explained by Eu 2+ fractionation in the plagioclase and/or alkali feldspar that crystallizes before or during zircon formation from the magma [1,43,44].Therefore, a pronounced negative Eu anomaly is commonly found in zircons from crustal-derived felsic rocks [16,44], whereas no Eu anomaly is found in zircon from feldspar free rocks, such as little to non-fractionated or mantle-derived rocks [45].Note also that zircons of mantle origin usually have low REE and Y concentrations (zircon in kimberlite from southern Africa with TREE = 5-39 ppm and Y = 11-74 ppm [46]; zircon from Jwaneng kimberlite TREE up to 12 ppm and Y up to 23 ppm [47]), whereas zircons in crustal rocks are more enriched in REEs ranging from 250 ppm to 5000 ppm, with a 1500-2000 ppm average [1,16,48].Trace-element concentrations of zircon crystals investigated in this study (TREE = 25-309 ppm, Y = 26-392 ppm) shared chemical characteristics comparable to mantle-derived zircon, and had no Eu anomaly (Figure 7).Moreover, the provenance discrimination plots proposed by Belousova et al. [16] and Shnukov et al. [41] revealed that the trace-element composition of Vietnamese zircon megacrysts were consistent with zircon typically found in kimberlite, syenite, carbonatites (Figure 8A-C), or alkaline rocks and alkaline metasomatites of alkaline complexes (Figure 8D).As the observed trace-element chemical characteristics ambiguously meet with those of several discrimination fields, it is hard to define an exact source composition.In fact, zircon samples from carbonatites often fall into various discrimination fields [49,50].For example, Saava et al. [49] reported that only 45 out of 100 analyzed zircon grains extracted from a carbonatite body meet the trace-element signatures that indicate a carbonatitic source, while all others suggest an alkaline or ultramafic source.This compositional variation has been interpreted to result from the interaction between carbonatitic melt and co-magmatic silicate sources that broaden the chemical variation of carbonatitic zircon.In this study, total REE contents of zircon typically varied in the range of 25 to 309 ppm.These concentrations were slightly higher than those found in the kimberlitic zircon that had typically less than 50 ppm [46,47], but had a much lower concentration than the one from zircon of corundum bearing syenites (up to 3500 ppm [2]).Low concentrations of U and Th (Figure 8B), and low Nb and Ta concentrations are also typical of the zircon that originates from very little, to non-fractionated, Si-poor melts [16].Trace-element chemical signatures of the zircon megacrysts of south-central Vietnam were found to be very comparable to those of zircon from Ratanakiri, Northeast Cambodia (Figure 7 in Cong et al. [10]; Figure 7 in this study) and other places such as New Zealand [11], East Australia [18], eastern and northeast China [5,8].Results of a detailed study by Cong et al. [10] and Piilonen et al. [22] using δ 18 O and 176 Hf/ 177 Hf isotopes clearly demonstrated that zircon from Ratanakiri, Cambodia, had a mantle origin and may be derived from metasomatized, partially melted lithospheric mantle material (such as peridotite, harzburgite), with strong carbonatitic geochemical fingerprint.In fact, experimental studies by Foley et al. [51], have demonstrated that mantle metasomatism, as induced by the presence of CO 2 and H 2 O, lowers the solidus temperature at upper mantle conditions and promotes partial melting of the peridotites that results in melts with carbonatitic compositions at low degree, and carbonated silicate melts with higher degrees of melting.Early carbonate-dominant melts were characterized by a low Hf compatibility, that is consistent with low Hf concentrations obtained from zircon in this study (Figure 8D) [5,8,10,11].With increasing degrees of partial melting, Hf is more compatible in the carbonated silicate melt [51].
Owing to the close geographical, geochemical, and geochronological relation of Cambodian zircon from Ratanakiri with material from the Central Highlands in Vietnam, studied here, we considered a very similar genetic origin for the latter.This was further supported by the presence of primary inclusions.Euhedral inclusions of calcite and corundum, and subhedral olivine are present in the zircon host, that is free of fluid inclusions, tubular channels, which shows no indication of visible cracks or alteration (Figure 3).These findings were in accordance to those of Le Bas [52], who interpreted the presence of abundant primary carbonate and olivine inclusions in the zircon megacrysts to be indicative of the growth in a silica-undersaturated, carbonatite-like melt.The experimental work of Baldwin [53] demonstrates that corundum and olivine likely crystallize from carbonatitic melts during early phases.However, no carbonatite bodies were reported from the studied area in the Central Highlands, Vietnam.In fact, the only known occurrence of carbonatites in Vietnam has been found in South Nam Xe, Northwestern Vietnam.Those rift-related carbonatites have reported ages of 28-44 Ma (biotite K-Ar ages) and 30-32 Ma U-Th-Pb isochron ages [54,55], and hence, may not be considered as a potential source.Instead, we found U-Pb zircon ages that were consistent with periods of basaltic volcanic activities in southern Vietnam.Two distinct zircon populations with U-Pb ages of ca.6.5 Ma and ca.1.0 Ma were identified.Although zircon megacrysts may be potentially derived from carbonatitic, metasomatized mantle material and do not primarily relate to alkali basalts, the zircon ages obtained were consistent with bi-episodal eruptive events at Buon Ma Thuot, Dak Lak (5.8-1.67Ma), Pleiku (4.3-0.8Ma), and the younger basalt eruptions exposed at Xuan Loc (0.83-0.44 Ma) (Figure 1B) [4].Zircon megacrysts from North-East Cambodia not only occur geographically close to the megacrysts found in the present study, but also share a similar U-Pb age connected to the younger period of volcanic activity in the region (0.98 ± 0.04 Ma, Cong et al. [10]; 0.88 ± 0.22 Ma to 1.56 ± 0.21 Ma, Piilonen et al. [22]).Obtained zircon U-Pb ages were just slightly older than the eruption events that were dated based on the K-Ar ages for the whole rock [4].This might be due to the difference in closure temperatures for U-Pb ages of zircon and K-Ar ages from the basalt, that indicated rapid extraction of zircon from its formation source, and the comparably fast uplift by the alkali-basalt eruption events [17,22].Note, however, that the rim zones of resorption and re-growth (visualized by PL images in Figure 5) gave the very same ages as the zircon core regions within statistical errors.We excluded a potential alteration process to explain this texture.Late-stage, fluid-driven dissolution-reprecipitation typically results in geochronological resetting [56], and geochemical signatures typically indicate a progressive REE fractionation, e.g., a more pronounced negative Eu-anomaly.As both criteria are not verified for the overgrowth zones found in the zircon megacrysts, in this study, we assumed that multiple resorption and re-growth processes took place in the deep sub-crustal "carbonatitic" source magma chambers at temperatures higher than U-Pb closure (>900 • C; Cherniak et al. [57]).Multiple, late-Cenozoic, relatively fast uplifts of basaltic melts induced by tectonic weakening along shear-zones in the tectonically active region, may have either caught up with xenolithic zircon megacrysts from and/or mixed with pre-existent carbonatite-dominated melts generated by CO 2 and/or H 2 O metasomatization of the upper mantle material.Inclusions of hematite, baddeleyite, quartz, and feldspar, as well as some ovoid or droplet carbonate inclusions along fractures or along channels, that are associated with fluid inclusions of liquid H 2 O and vapor bubbles of CO 2 , are considered to be of secondary origin (Figure 3).These phases are well known to occur upon zircon alteration.Particularly, the formation of baddeleyite due to the metasomatic alteration of zircon induced by Ca-bearing fluids has been proven by the experimental works of Lewerentz [58].Therefore, the occurrences of these phases along fluid inclusions in close proximity to fractures, within tubular channels and/or vesicles are interpreted to result from late stage circulation of oxidative fluids in the ascending alkali-basaltic magma.Finally, weathering processes eroded zircon megacrysts from the alkali-basalt hosts in the Central Highlands and resulted in the deposition and enrichment of the nearby alluvial zircon placers found today.

Conclusions
Crystal morphology and internal texture of gem-quality, placer zircon megacrysts of the Central Highlands, Vietnam, are indicative of magmatic origin.The geochemical signatures and primary inclusions of calcite, olivine, and corundum indicate that zircon megacrysts might have crystallized from a carbonatite-dominant melt caused by the low-grade partial melting of metasomatized lithospheric mantle.Zircon resorption and re-growth textural features, observed with PL imaging, indicate an extended residence time in the sub-crustal magma chambers at variable REE-saturation levels, but at temperatures higher than closure of the U-Pb system (>900 • C).Zircon megacrysts were subsequently incorporated into ascending alkali basalts as xenocrysts.Two distinct populations with U-Pb ages at ca. 6.5 Ma and ca.1.0 Ma that correlate with the eruption of alkali basalt fields in the Central Highlands, were identified.Recent reports of the direct findings of megacrystic zircon in the geographically and genetically related basalts of Ratanakiri, Cambodia [22] strongly supports this close relationship of placer zircon with Cenozoic alkali basalts as a potential host.The basalt magma and/or late stage circulation of carbonate-rich fluids must have resulted in the entrapment of further secondary inclusions like CO 2 -H 2 O fluids, baddeleyite, quartz, hematite, and feldspar along the cracks and tubular channels.

Figure 1 .
Figure 1.(A) Location of various alluvial zircon deposits in Cambodia and Vietnam within South-East Asia [12,24].(B) Distribution of Neogene and Quaternary basalts with their K-Ar and Ar-Ar ages (after [4,6]).Reproduced with permission from all authors.Green circles indicate the sampling localities.

Figure 1 .
Figure 1.(A) Location of various alluvial zircon deposits in Cambodia and Vietnam within South-East Asia [12,24].(B) Distribution of Neogene and Quaternary basalts with their K-Ar and Ar-Ar ages (after [4,6]).Reproduced with permission from all authors.Green circles indicate the sampling localities.

Figure 2 .
Figure 2. (A) Typical large, gem-quality, cut zircon samples from the Central Highlands are available on the gem market.(B) Selected zircon megacrysts used in this study from the same region.Samples C05, C07, and C16 are from the Buon Ma Thuot and Xuan Loc alkali-basaltic field.Samples DL1, DL2, and DL3 (not shown) are from deposits in the Dak Lak alkali-basaltic field (Figure 1B).

Figure 2 .
Figure 2. (A) Typical large, gem-quality, cut zircon samples from the Central Highlands are available on the gem market.(B) Selected zircon megacrysts used in this study from the same region.Samples C05, C07, and C16 are from the Buon Ma Thuot and Xuan Loc alkali-basaltic field.Samples DL1, DL2, and DL3 (not shown) are from deposits in the Dak Lak alkali-basaltic field (Figure 1B).

Figure 4 .
Figure 4. Cathodoluminescence images of the selected zircon samples from the Central Highlands, resemble zircon internal textural features that are typical for zircon of magmatic growth, e.g., oscillatory, sector, and growth zoning.Note, however, that whitish striations (indicated by black arrows) in CL contrasts obtained in the images of samples C01, C07, and F01 are due to mechanicallyinduced structural defect centers emanating from cracks and break-outs during polishing.

Figure 4 .
Figure 4. Cathodoluminescence images of the selected zircon samples from the Central Highlands, resemble zircon internal textural features that are typical for zircon of magmatic growth, e.g., oscillatory, sector, and growth zoning.Note, however, that whitish striations (indicated by black arrows) in CL contrasts obtained in the images of samples C01, C07, and F01 are due to mechanically-induced structural defect centers emanating from cracks and break-outs during polishing.

Figure 5 .
Figure5.Laser-induced photoluminescence hyperspectral images of zircon samples from the Central Highlands (C16 and DL3) and from Ratanakiri, Cambodia (Rata).The integrated intensity of the Dy 3+ ( 4 F9/2 → 6 H13/2) emission (compare Figure6) is the plotted spectral parameter (grey-scale).Traceelement concentrations of Dy from multiple LA-ICP-MS spots on the grains, were used for external calibration to correlate the Dy concentrations with its PL response.Concentrations given color-coded were back-calculated using the latter correlation (see Figure6) to extrapolate the Dy distribution in regions of unknown trace-element chemistry.Besides oscillatory and sector zoning, several growth zones may be distinguished that were interpreted to result from resorption and re-growth events, as indicated by the cross-cut zonation patterns, abruptly changing (decreasing) Dy concentrations, and the appearance of large holes and cavities (see arrows).

Figure 5 .
Figure5.Laser-induced photoluminescence hyperspectral images of zircon samples from the Central Highlands (C16 and DL3) and from Ratanakiri, Cambodia (Rata).The integrated intensity of the Dy 3+ ( 4 F 9/2 → 6 H 13/2 ) emission (compare Figure6) is the plotted spectral parameter (grey-scale).Trace-element concentrations of Dy from multiple LA-ICP-MS spots on the grains, were used for external calibration to correlate the Dy concentrations with its PL response.Concentrations given color-coded were back-calculated using the latter correlation (see Figure6) to extrapolate the Dy distribution in regions of unknown trace-element chemistry.Besides oscillatory and sector zoning, several growth zones may be distinguished that were interpreted to result from resorption and re-growth events, as indicated by the cross-cut zonation patterns, abruptly changing (decreasing) Dy concentrations, and the appearance of large holes and cavities (see arrows).

Figure 6 .
Figure 6.(A) Laser-induced PL spectra (laser-excitation wavelength λexc = 473 nm) obtained from three different measurement spots on the zircon sample DL3 (spots 4, 5, and 7; Figure 5).(B) The integrated intensity of the most intensive emission of Dy 3+ in the spectral range 17,000-17,450 cm −1 (see dotted lines in sub-image (A)) from multiple spots on samples DL3, C16, and Rata plotted against the Dy trace-element concentration as obtained from LA-ICP-MS (TableS1).Linear regression from the data pairs PL vs. Dy conc.(R² = 0.95) was used to calibrate hyperspectral maps presented in Figure5.

Figure 6 .
Figure 6.(A) Laser-induced PL spectra (laser-excitation wavelength λ exc = 473 nm) obtained from three different measurement spots on the zircon sample DL3 (spots 4, 5, and 7; Figure 5).(B) The integrated intensity of the most intensive emission of Dy 3+ in the spectral range 17,000-17,450 cm −1 (see dotted lines in sub-image (A)) from multiple spots on samples DL3, C16, and Rata plotted against the Dy trace-element concentration as obtained from LA-ICP-MS (TableS1).Linear regression from the data pairs PL vs. Dy conc.(R 2 = 0.95) was used to calibrate hyperspectral maps presented in Figure5.

Figure 7 .
Figure 7. Chondrite-normalized plot of REE concentrations of zircon samples from the Central Highlands compared to samples from Ratanakiri, Cambodia.REE composition of C1 chondrite is described by Sun and McDonough [42].All samples showed a steep slope rise from La to Lu, towards heavy REEs.A positive pronounced Ce-anomaly, but no Eu-anomaly was detected.

Figure 7 .
Figure 7. Chondrite-normalized plot of REE concentrations of zircon samples from the Central Highlands compared to samples from Ratanakiri, Cambodia.REE composition of C1 chondrite is described by Sun and McDonough [42].All samples showed a steep slope rise from La to Lu, towards heavy REEs.A positive pronounced Ce-anomaly, but no Eu-anomaly was detected.

Figure 9 .
Figure 9. Summary of U-Pb geochronology data obtained by LA−ICP−MS (see details in TableS2 and S3).(A) Representative concordia plots of samples C05, C07, and C16.Ellipses give a 2σ error, and red ellipses correspond to measurement spots placed in the core region of single crystals.(B)206 Pb/ 238 U ages of the samples analyzed with the mean weighted age given for each sample.Error bars represent 2σ.Note that the ages obtained do not differ systematically amongst the measurement spots placed in the different growth zones within the individual samples.Two distinct sample populations were clearly differentiated (ca.6.5 Ma and ca. 1 Ma).
Figure 9. Summary of U-Pb geochronology data obtained by LA−ICP−MS (see details in TableS2 and S3).(A) Representative concordia plots of samples C05, C07, and C16.Ellipses give a 2σ error, and red ellipses correspond to measurement spots placed in the core region of single crystals.(B)206 Pb/ 238 U ages of the samples analyzed with the mean weighted age given for each sample.Error bars represent 2σ.Note that the ages obtained do not differ systematically amongst the measurement spots placed in the different growth zones within the individual samples.Two distinct sample populations were clearly differentiated (ca.6.5 Ma and ca. 1 Ma).

Figure 9 .
Figure 9. Summary of U-Pb geochronology data obtained by LA-ICP-MS (see details in TablesS2 and S3).(A) Representative concordia plots of samples C05, C07, and C16.Ellipses give a 2σ error, and red ellipses correspond to measurement spots placed in the core region of single crystals.(B)206 Pb/ 238 U ages of the samples analyzed with the mean weighted age given for each sample.Error bars represent 2σ.Note that the ages obtained do not differ systematically amongst the measurement spots placed in the different growth zones within the individual samples.Two distinct sample populations were clearly differentiated (ca.6.5 Ma and ca. 1 Ma).
Figure 9. Summary of U-Pb geochronology data obtained by LA-ICP-MS (see details in TablesS2 and S3).(A) Representative concordia plots of samples C05, C07, and C16.Ellipses give a 2σ error, and red ellipses correspond to measurement spots placed in the core region of single crystals.(B)206 Pb/ 238 U ages of the samples analyzed with the mean weighted age given for each sample.Error bars represent 2σ.Note that the ages obtained do not differ systematically amongst the measurement spots placed in the different growth zones within the individual samples.Two distinct sample populations were clearly differentiated (ca.6.5 Ma and ca. 1 Ma).