Received 7 December 2015; accepted 19 January 2016; published 22 January 2016
The Mushgai khudag alkaline volcanic-plutonic complex with REE-Sr-Fe-F-P mineralization is located in South-Gobi, Mongolia. Russian and Mongolian Joint Geological Expedition discovered the Mushgai khudag complex in early 1970’s. Since that, number of new studies about alkaline rocks, carbonatites, and REE mineralization the deposit has been completed  - .
Geological and geochemical surveys demonstrated a wide occurrence of ore-bearing rocks containing rare earth elements with economic grades, especially light lanthanides and also strontium and barium  . The alkaline complex is composed of various types of LREE enriched rocks significantly apatite rocks and monazite-rich apatite rock. LREE-rich apatites contain up to 14 % LREE2O3, while monazite up to 68% LREE2O3 (LREE2O3 represent sum of La2O3, Ce2O3, Nd2O3, and Pr2O3).
In this study, we discuss the chemistry of LREE-rich apatite crystals from the Mushgai khudag complex. Apatite from the Monazite zone, which has unusual crystal texture with high LREE composition, is main contribution of the study because such type of apatite with individual monazite is rarely found in naturally. Three types of apatite crystals have been examined. Crystal-1 is euhedral and sometimes associated with celestine grains, and crystal-2 is subhedral and brecciated, whereas crystal-3 has an unusual idiomorphic crystal that demonstrates compositional zoning texture and associated with sulfur-rich monazite grains.
Monazite associated with apatite was studied in natural occurrences  . Moreover, experimental studies conducted that monazite forms from the host apatite via the action of fluids   . In the apatite structure, (Y + REE) are principally charge balanced through the coupled substitutions: 1) Si4+ + (Y + REE)3+ = P5+ + Ca2+, 2) Na+ + (Y + REE)3+ = 2Ca2+  . Apatite from the monazite zone exhibits the compositional zoning texture. The individual monazite band formed in apatite. Interestingly, we note that monazite found in this study is enriched in S (0.93 - 14.72), Sr and Ca (≤7.21 and ≤8.25 wt.% oxides, respectively). Incorporation of these elements in the monazite structure involves the complex substitution: Ce3+ + P5+ = Ca2+ (Sr2+) + S6+ to form of sulfur-rich monazite  .
2. Geological Background
The Mushgai khudag alkaline complex is associated with late Jurassic alkaline magmatic activity close to the Main Mongolian lineament along where Mesozoic rift structures have been developed. The Mushgai khudag alkaline complex is hosted in the Silurian-Devoniun Mandal-oboo Formation’s turbidite sediments, lower Devoniun Orgol Formation’s carbonate sediments, and mid Devoniun Bot-uuli Formation’s volcanic sediments  (Figure 1).
Figure 1. Geological map of the Mushgai khudag deposit (modified from exploration report, “Mongol gazar” LLC, 2010).
Late Jurassic alkaline complex consists of trachytic volcanics and alkaline syenites, melanephelinite, and melaleucitite. K-Ar ages of alkaline rocks vary from 161 to 121 Ma  and Rb-Sr isochron ages is 139 Ma, respectively  .
REE mineralization zones occur peripheries of the late Jurassic alkaline intrusives. The peripheries of intrusives and their roofs are characterized by distribution of quartz-syenite, micro-syenites and porphyritic syenite dykes. Large rectangular-shaped porphyritic phenocrysts in porphyritic syenite and quartzsyenite demonstrate a rapakivi texture where phenocrysts of K-feldspar (Ab5-27 An0-3 Or70-94) are rimmed by plagioclase (Ab68-97 An2-20 Or0-11). Coarse-grained K-feldspar phenocrysts compose main groundmass of quartzsyenite and porphyritic syenite. Round-shaped porphyritic phenocrysts of micro-syenite show an antirapakivi texture with phenocrysts of plagioclase (Ab65-73 An18-25 Or5-10) rimmed by K-feldspar (Ab3-14 Or86-94). Fine-grained mostly K-feldspar (Ab7-8 Or79-93) phenocrysts composed groundmass of micro-syenite. The quartz-syenite consist of (<5 vol.%) quartz grains that shows poiklitic texture  .
The carbonatite zone in central part of alkaline complex has approximately 1 km length and 500 meter width. The largest mineralized zone is composed of low-grade REE carbonatite rocks (total LREE 1.64 - 2.11 wt.%). The Magnetite zone is composed of massive magnetite rocks, and low-grade LREE-bearing apatite veins occur in massive magnetite rocks’s cracks (LREE up to 0.11 wt.%). The Apatite zone is composed of various type apatite rocks with high-grade LREE (5 - 10 wt.%). Monazite zone composed of high-grade apatite rocks that contain LREE up to 10 wt.%. The apatite rocks from monazite zone demonstrate zoning apatite with individual sulfur-rich monazite grain (Figure 1).
3. Mineral Chemistry
3.1. Analytical Method
Mineral compositions were determined using an electron microprobe (JEOL JXA-8600 SX) with energy dispersive spectroscopy (INCA-6025, Oxford Instruments) at Korea University, Seoul, South Korea. The microprobe was operated with an accelerating voltage = 15.0 kV, beam current = 3.0 nA, live time = 100 seconds, and probe diameter = 3 μm. Used standards materials: e.g. quartz (SiO2) for Si, corundum (Al2O3) for Al, magnetite (Fe3O4) for Fe, rutile (TiO2) for Ti, manganosite (MnO) for Mn, wollastonite (CaSiO3) for Ca, jadeite (NaAlSi2O6) for Na, for sanidine (KAlSi3O8) for K, Ni-oxide (NiO) for Ni, fluorite (CaF2) for F, apatite (Ca5(PO4)8(OH)) for P, REEs (LaB6, NdB6, CaB6) for LREE with SPI and JEOL Company.
Three types of apatite crystals have been analyzed. In BSE image, euhedral elongated large apatite crystal-1 contains numerous of celestine inclusions (Figure 3(a)). Apatite crystal-1 from Apatite zone contains LREE (up to 3 wt.% Ce2O3). Such type of apatite crystals are mostly euhedral and contained LREE (0.72 - 5.52 wt.%), (Table 1).
Subhedral, brecciated apatite crystals are included in apatite crystal-2 type. In BSE image, apatite crystal-2 shows round shape and discerns to dark and light areas (Figure 3(b)). Dark areas of the apatite crystal-2 contain LREE2O3 (11.44 wt.%), whereas light areas are contained LREE2O3 (13.27 wt.%), (Table 1).
Apatite crystal-3 demonstrates euhedral compositional zoning textures. In BSE image, the most compositional zoning apatite ranges from LREE depleted in core1 (moderate), core2 (dark) areas (up to 10 wt.% Ce2O3 + La2O3 + Nd2O3), and LREE strong enriched in rim1 (light) area (up to 14.64 wt.% Ce2O3 + La2O3 + Nd2O3) (Figure 3(c), Figure 3(d)). Representative apatite compositions from the REE ores and host rocks are given in Table 1. The compositional variations of apatite are plotted in the solid solution among F-OH-Cl diagram  (Figure 2). The REE-bearing apatites from the Mushgai khudag deposit are mostly fluorapatite to hydroxyl-bearing fluorapatite with variable REE contents.
The fine-grained subhedral monazite grains, associated with apatite, mostly occur as a rim of apatite crystals at Monazite zone (Figure 3(c), Figure 3(d)). Representative monazite compositions from the REE ores are given
Table 1.Repressentative chemical composition of LREE-rich apatite from the monazite zone,apatite zone,carbonatite zone,and sysenite.
Figure 2. Apatite is showing solid solution among F-OH-Cl end-members, (after Pan and Fleet, 1995). The apatite re- presents fluorapatite and hydroxyl-bearing fluorapatite.
in Table 2. The monazite associated with fluorapatite or hydroxyl-bearing fluorapatite. The apatites show idiomorphic crystal grains, and compositional zoning texture. They have commonly LREE-poor core2 (dark area), LREE-rich core1 (moderate), rim1 (light areas). The apatites in the syenites, apatite zone and carbonatite zone have a high content of phosphate and calcium with depletion in Si, Na, and LREE, whereas apatite from monazite zone has a low content of P and Ca with increasing of Si, Na, and LREE (Figure 4).
The substitution of Si and/or Na from the apatite structure without coupled removal of (Y + REE)3+ the resulting charge imbalance, due to the coupled substitution reactions in apatite from Monazite zone namely:
As pointed by Harlov (2013) and Pan (1997), the substitutions encouraged the removal of Ca2+ and the subsequent growth of monazite and/or xenotime inclusions within the apatite grain body via dissolution-precipita- tion from the P5+ available in the immediate vicinity of the growing crystal and (Y + REE)3+ transported in from over a much broader area    .
The intensive coupled substitutions with (Y + REE)3+ contributed the nucleation of the individual sulfur-rich monazite between the rim1 and the rim 2 area in the apatite from Monazite zone (Figure 3(d)).
Principally, in the monazite structure, rare earth elements are charge balanced by substitution scheme: 2Ce3+ = Ca2+(Sr2+) + Th4 that known as Th-bearing barbantite and Th, Ca, Sr-bearing strontiobrabantite    .
This substitution scheme is not applicable for monazite from the Mushgai khudag, Mongolia, because relative depletion of Th content. Th was not detected due to below detection limit. Nevertheless, samples from the monazite-bearing apatite at the Monazite zone demonstrate Th content around (300 - 400) ppm in XRF data  . Therefore another mechanism of substitution as inferred by Chakhmouradain (1999) may account for the monazite from the Mushgai khudag complex, that Ce3+ + P5+ = Ca2+(Sr2+) + S6+   .
The monazite in the Mushgai khudag complex contains high S in nature (1.5 - 14.72 wt.% SO3), and enriched in Sr and Ca (up to 7.21 and ≤8.25 wt.% oxides, respectively). Our data show S and P, S, and LREE moderate negative trend (Figure 5(a), Figure 5(c)), whereas S and Ca2+, Sr2+ divalent cations represent weak positive pattern (Figure 5(b)), implying that Ce3+ + P5+ = Ca2+(Sr2+) + S6+ scheme of substitution should be a major mechanism of incorporation of Ca, Sr and S in the composition of monazite from the Mushgai khudag complex, South Mongolia.
Figure 3. (a) Euhedral apatite crystal-1 from apatite zone and celestine inclusions associated with apatite crystal; (b) Subhedral apatite crystal-2 from apatite zone and celestine. Apatite crystal-2 can be discern to dark and light areas; (c) Euhedral apatite crystal-3 associated with monazite grains; (d) Apatite crystal-3 can be discern to core1 (c1) moderate, core2 (c2) dark, and rim1 (r1), rim2 (r2) light areas. Monazite grain occurs between rim1 (r1) and rim2 (r2); (e) Subhedral very fine-grained monazites associated with Celestine; (f) Elongated euhedral very fine grained two monazite grains associated with celestine.
Table 2. Representative chemical composition of monaztite from the Mushgai khudag complex.
Note: Cations in crystal-chemical formulae normalized to 4 oxygen atoms.
The mobility of the rare earth elements in hydrothermal fluid is strongly depend on the relative concentrations of possible ligands, such as F−, Cl−, OH−, , and, and the pH of the solution. Also, hydrothermal fluid’s temperature and pressure are important  .
It has been reported that (= 200 ppb) concentrations, in acid waters at Yellowstone hydrothermal system (USA)  . REE content tends to be increased in Low-pH waters than neutral or alkaline solutions. Low-pH, sulfate-rich fluids with the highest REE concentrations in hydrothermal fluids have been found at the Valles-caldera (= 339 ppb),  .
Moreover, Saline magmatic solutions with high-temperature can effective carry REE. Also, the solutions can transport a large amount of REE (1290 ppm) in the ancient hydrothermal system  . Wood’s (1990) experimental data and theoretical predictions suggest that the sulfate ions form relatively strong complexes with REE. At the temperatures ≥200˚C and in the absence of other ligands, sulfate complexes will be the dominant form of REE transport, even at free sulfate environment  .
Nevertheless, Brookins (1989) noted that the REE with sulfate complexes is considerably important only in the absence of other potentially strong ligands. If strong complex-forming ligands such as fluoride, carbonate, or possibly the phosphate were present, the suggested sulfate complexes with REE would probably not be important because of the strong affinity of the hard REE ions for hard ligands such as F−, OH−, and, especially at higher temperatures  . It has been suggested that the REE cations are classified as hard ions, as ions with a small size and high charge. Hard cations bind preferentially with hard ligands  .
Figure 4. Plots of cation values as atoms per formula unit (a.p.f.u) reflecting coupled substitution reaction (a) Si4+ + LREE3+ = P5+ + Ca2+, (b) Na+ + LREE3+ = 2Ca2+, for the apatite from monazite zone, that discern to core2 (dark), core1 (moderate), rim1 (light) areas. Apatite from carbonatite zone, apatite zone, and in syenite (after Harlov et al., 2003).
Figure 5. Variation in minor components “in Ca-Sr-rich monazite from the Mushgai khudag deposit” (a) P5+ versus S6+, (b) (Ca2+ + Sr2+) versus S6+, and (c) LREE3+ versus S6+ (after Chakhmourdain and Mitchell, 1999).
In this study, we suggest sulfate ligand, which encouraged the rare earth elements mobility, based on the chemical characteristics of monazite-bearing apatite in the Mushgai khudag complex. It is no doubt that the hard ligand phosphate is abundant with a high alkaline condition of the first stage of magmatism in the Mushgai khudag complex. However, the hard ligands fluoride and carbonate were dominant with the carbonatite magmatism. In case, the sulfate ligand complexed with REE in the Mushgai khudag complex, it could be the other ligands abundance would have been minor. In addition, the pH condition should be low. Thus, the sulfate ligand complexing with REE may have been related to the last stage of hydrothermal activity (metasomatism) in the Mushgai khudag complex. Our result show that the monazite in the Mushgai khudag complex. Our result show that the monazite contains CaSO4 (anhydrite), SrSO4 (celestine) component. In other words, the monazite was formed by sulfate-rich fluid induced metasomatism. The S6+ versus LREE3+ diagram show the S6+ moderate positive correlation with LREE3+. The S6+ content is increased with LREE3+ increase in the apatite from the Mushgai khudag complex. Especially the sulfur-rich monazite forming (core2), (core1), (rim1) areas are strongly enriched in S content suggesting that the sulfur-rich monazite in the Mushgai khudag complex strongly depends on the sulfate ligand (Figure 6).
In the Mushgai khudag deposit, LREE-bearing minerals are dominantly fluorapatite, hydroxyl-bearing fluorapatite, and monazite. Fluorapatite and hydroxyl-bearing fluorapatite can be found in all mineralization zones, whereas monazite distributed limited area. S and Sr-rich individual monazite grains have been found in the rim of fluorapatite and hydroxyl-bearing fluorapatite in the samples from the Monazite zone. It can be explained that the Mushgai khudag rocks were affected by strong REE mineralization due to various stages of magmatic and hydrothermal segregations. REEs have a strong mobility in the last stage of metasomatism due to sulfate-rich fluid complexed with REE. Apatite from monazite zone shows diverse zoning texture, and individual high S content monazite grains are formed along the rim of apatite crystals, and is well described in the present study.
Figure 6. Plots of cation values as atoms per formula unit (a.p.f.u), S6+ versus LREE3+ for the apatite from monazite zone, discern to core2 (dark), core1 (moderate), rim1 (light) areas. Apatite from carbonatite zone, apatite zone, and in syenite.
Such zonation of apatite with monazite is a relatively rare case in terms of REE mineralization with phosphate.
This work was supported by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government Ministry of Knowledge Economy (20152510101890) and the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP), (2014R1AEAEA01007294). An earlier manuscript of this paper was improved by R. Giere and a number of referees who suggested valuable recommendations. First stage of the manuscript improvement was made by Professor B. Munkhtsengel and Mr, D. Batbold. We wish to thank Professor O. Gerel for last critical reading of the English content and greatly improved this paper. Finally, thanks to the fellows of the “Mongol Gazar” company for their help to provide research opportunity, especially Mrs, M. Delgertsetseg for her help in field work.
 Baatar, M., Ochir, G., Kynicky, J., Iizumi, S. and Comi-Ciaramonti, P. (2013) Some Notes on the Lugiin Gol, Mushgai-Khudag- and Bayan-Khoshuu Alkaline Complexes, Southern Mongolia. International Journal of Geosciences, 4, 1200-1214.
 Harlov, D.E. and Forster, H.J. (2003) Fluid-Induced Nucleation of (Y+REE)-Phosphate Mineral within Apatite: Nature Experiment. Part II. Fluorapatite. American Mineralogist, 88, 1209-1229.
 Harlov, D.E., Forster, H.J. and Schmidt, C. (2003) High P-T Experimental Metasomatism of a Fluorapatite with Significant Britholite and Fluorsellestadite Components: Implications for LREE Mobility during Granulite-Facies Metamorphism. Mineralogical Magazine, 67, 61-72.
 Harlov, D.E., Forster, H.J. and Nijland, T.G. (2002) Fluid-Induced Nucleation of (Y + REE) Phosphate Minerals in Apatite: Nature and Experiment. Part I. Chlorapatite. American Mineralogist, 87, 245-261. http://dx.doi.org/10.2138/am-2002-2-306
 Pan, Y. and Fleet, M.E. (2002) Compositions of the Apatite-Group Minerals: Subtitution Mechanisms and Controlling Factors. Reviews in Mineralogy and Geochemistry, 48, 13-49.
 Chakhmouradian, A.R. and Mitchell, R.H. (1999) Niobianilmenite, Hydroxylapatite and Sulfur-Rich Monazite: Alternative Hosts for Incompatible Elements in Calcite Kimberlite from Internatsional’Naya, Yakutia. The Canadian Mineralogist, 37, 1177-1189.
 Watt, G.R. (1995) High-Thorium Monazite-(Ce) Formed during Disequilbrium Melting of Metapelties under Granulite-Facies Conditions. Mineralogical Magazine, 59, 735-743.
 Chakhmouradian, A.R. and Mitchel, R.H. (1998) Lueshite, Pyrochlore and Monazite-(Ce) from Apatite-Dolomite Carbonatite, Lesnaya Varaka Complex, Kola Peninsula, Russia. The Canadian Mineralogist, 62, 769-782.
 Jaroslav, P., Martin, O., Peter, B. and Bartosz, B. (2010) Methamorphic-Hydrothermal REE Minerals in the Bacuch Magnetite Deposit, Western Carpathias, Slovakia: (Sr, S)-Rich Monazite-(Ce) and Nd-Dominant Hingganite. The Canadian Mineralogist, 48, 81-94. http://dx.doi.org/10.3749/canmin.48.1.81
 Giere, R. (1993) Transport and Deposition of REE in H2S-Rich Fluids: Evidence from Accessory Mineral Assemblages. Chemical Geology, 110, 251-268. http://dx.doi.org/10.1016/0009-2541(93)90257-J
 Palmer, M.A. (1998) Rare Earth Element Speciation in Geothermal Fluids from Yellowstone National Park, Wyoming, USA. Geochimica et Cosmochimica Acta, 62, 657-663. http://dx.doi.org/10.1016/S0016-7037(97)00367-0
 Michard, A. (1989) Rare Earth Element Systematics in Hydrothermal Fluids. Geochimica et Cosmochimica Acta, 53, 745-750. http://dx.doi.org/10.1016/0016-7037(89)90017-3
 Banks, D.Y. (1994) REE Composition of an Aqueous Magmatic Fluid: A Fluid Inclusion Study from the Captain Pluton, New Mexico. Chemical Geology, 113, 22. http://dx.doi.org/10.1016/0009-2541(94)90070-1
 Wood, A. (1990) The Aquesous Geochemistry of the Rare-Earth Elements and Yttrium. Theoretical Predictions of Speciation in Hydrothermal Solutions to 350?C at Saturation Water Vapor Pressure. Chemical Geology, 88, 99-125.