skip to main content

Title: Cl-bearing fluorcalciobritholite in high-Ti basalts from Apollo 11 and 17: Implications for volatile histories of late-stage lunar magmas.
Abstract We report the occurrence of a previously unidentified mineral in lunar samples: a Cl-,F-,REE-rich silico-phosphate identified as Cl-bearing fluorcalciobritholite. This mineral is found in late-stage crystallization assemblages of slowly cooled high-Ti basalts 10044, 10047, 75035, and 75055. It occurs as rims on fluorapatite or as a solid-solution between fluorapatite and Cl-fluorcalciobritholite. The Cl-fluorcalciobritholite appears to be nominally anhydrous. The Cl and Fe2+ of the lunar Cl fluorcalciobritholite distinguishes it from its terrestrial analog. The textures and chemistry of the Clfluorcalciobritholite argue for growth during the last stages of igneous crystallization, rather than by later alteration/replacement by Cl-, REE-bearing metasomatic agents in the lunar crust. The igneous growth of this Cl- and F-bearing and OH-poor mineral after apatite in the samples we have studied suggests that the Lunar Apatite Paradox model (Boyce et al. 2014) may be inapplicable for high-Ti lunar magmas. This new volatile-bearing mineral has important potential as a geochemical tool for understanding Cl isotopes and REE chemistry of lunar samples.
; ;
Award ID(s):
Publication Date:
Journal Name:
American Mineralogist
Page Range or eLocation-ID:
255 to 261
Sponsoring Org:
National Science Foundation
More Like this
  1. The textures of outcrop and near-surface exposures of the massive magnetite orebodies (>90 vol % magnetite) at the Plio-Pleistocene El Laco iron oxide-apatite (IOA) deposit in northern Chile are similar to basaltic lava flows and have compositions that overlap high- and low-temperature hydrothermal magnetite. Existing models— liquid immiscibility and complete metasomatic replacement of andesitic lava flows—attempt to explain the genesis of the orebodies by entirely igneous or entirely hydrothermal processes. Importantly, those models were developed by studying only near-surface and outcrop samples. Here, we present the results of a comprehensive study of samples from outcrop and drill core that require a new model for the evolution of the El Laco ore deposit. Backscattered electron (BSE) imaging, electron probe microanalysis (EPMA), and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) were used to investigate the textural and compositional variability of magnetite and apatite from surface and drill core samples in order to obtain a holistic understanding of textures and compositions laterally and vertically through the orebodies. Magnetite was analyzed from 39 surface samples from five orebodies (Cristales Grandes, Rodados Negros, San Vicente Alto, Laco Norte, and Laco Sur) and 47 drill core samples from three orebodies (Laco Norte, Laco Sur, andmore »Extensión Laco Sur). The geochemistry of apatite from eight surface samples from three orebodies (Cristales Grandes, Rodados Negros, and Laco Sur) was investigated. Minor and trace element compositions of magnetite in these samples are similar to magnetite from igneous rocks and magmatic-hydrothermal systems. Magnetite grains from deeper zones of the orebodies contain >1 wt % titanium, as well as ilmenite oxyexsolution lamellae and interstitial ilmenite. The ilmenite oxyexsolution lamellae, interstitial ilmenite, and igneous-like trace element concentrations in titanomagnetite from the deeper parts of the orebodies are consistent with original crystallization of titanomagnetite from silicate melt or high-temperature magmatic-hydrothermal fluid. The systematic decrease of trace element concentrations in magnetite from intermediate to shallow depths is consistent with progressive growth of magnetite from a cooling magmatic-hydrothermal fluid. Apatite grains from surface outcrops are F rich (typically >3 wt %) and have compositions that overlap igneous and magmatic-hydrothermal apatite. Magnetite and fluorapatite grains contain mineral inclusions (e.g., monazite and thorite) that evince syn- or postmineralization metasomatic alteration. Magnetite grains commonly meet at triple junctions, which preserve evidence for reequilibration of the ore minerals with hydrothermal fluid during or after mineralization. The data presented here are consistent with genesis of the El Laco orebodies via shallow emplacement and eruption of magnetite-bearing magmatic-hydrothermal fluid suspensions that were mobilized by decompression- induced collapse of the volcanic edifice. The ore-forming magnetite-fluid suspension would have rheological properties similar to basaltic lava flows, which explains the textures and presence of cavities and gas escape tubes in surface outcrops.« less
  2. The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X‑ray absorption near-edge structure spectroscopy (S-XANES) where S6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms. One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, themore »energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel. The thermodynamics of incorporation of disulfide and bisulfide into apatite are evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S are also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS2(aq)(2–n) and HnS(aq)(2–n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite and decreases for hydroxylapatite as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended also to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.« less
  3. The Gallinas Mountains district located in Lincoln and Torrance Counties, New Mexico, is host to hydrothermal REE-bearing fluorite veins and breccia deposits. The rare earth elements (REE) are found in bastnäsite-(Ce) ([La,Ce]CO3F) which is also the primary ore mineral mined in several important carbonatite deposits (e.g. Mountain Pass in California; Bayan Obo in China). Minor production of REE, fluorite, Cu, Pb, Zn, Ag, and Fe has been recorded in the Gallinas Mountains district between the early 1900s and the 1950s. The REE-bearing fluorite veins and breccias are hosted in Permian sedimentary rocks as well as genetically related trachyte/syenite sills and dikes emplaced between 28-30 Ma. Previous studies have described the REE occurrences in the Gallinas Mountains but the controls of hydrothermal processes on the transport and deposition of REE in the district remain unclear. In this study, we combine microtextural observations with mineral and bulk rock chemistry of hydrothermal REE-bearing fluorite veins and breccias to determine the vein types, alteration styles and establish a detailed mineral paragenesis. The goal of this study is to determine lithogeochemical vectors towards REE enriched zones in the district by linking thin section and deposit scale observations with mineral and bulk rock geochemistry. This districtmore »is an exceptional natural laboratory for studying the role of hydrothermal processes for transport/deposition of REE in an alkaline F-rich magmatic- hydrothermal system because very few deposits worldwide have such well-preserved and exposed geology. Hand samples of hydrothermal veins and breccias containing fluorite ± calcite ± barite ± bastnäsite-(Ce) were collected from outcrops, prospect pits, and mine dumps. Optical microscopy was used to identify minerals and determine the textural features and crosscutting relationships of the different fluorite veins. The veins were classified into: i) hematite-fluorite veins; ii) barite- bearing bastnäsite-fluorite veins; iii) barite-bearing (fluorite)-calcite veins. Nearly all of the barite crystals in the fluorite veins display dissolution textures (skeletonized crystals) with infilling of mostly fluorite and minor calcite, suggesting that barite is part of an earlier paragenetic mineral assemblage. Bastnäsite-(Ce) is commonly found in veins containing barite and occurs either as disseminated crystals in the fluorite veins or together with fluorite infills around large barite crystals. A few of the barren fluorite-calcite veins display an intergrowth with euhedral barite crystals indicating that these could be part of an earlier barite paragenesis. These textural observations suggest a key control of REE mineralization in the Gallinas Mountains district by a coupled dissolution of barite-bearing (fluorite)-calcite veins and precipitation of later bastnäsite- fluorite veins. Geochemical bulk rock data collected from the New Mexico Bureau of Geology and Mineral Resources database were analyzed using the IMDEX ioGASTM program to definegeochemical signatures of rock types, alteration styles, and vein types. Preliminary data analysis indicates a positive correlation between Ba, F, and total rare earth oxides (TREO). These trends corroborate with the observed vein microtextures, suggesting that the interaction of a hydrothermal fluids with the barite-bearing (fluorite)-calcite veins represents a key process for defining geochemical vectors in the district.« less
  4. Societal demand for critical metals used in the high-tech and green industries has led to an increased interest in REEs associated with ore deposits. Hydrothermal mineralization of monazite (CePO4) in various REE deposits can display significant variations in REE mineralogy and rock chemistry. Monazite displays textural and REE compositional variations, such as those observed in the giant Bayan Obo REE deposit in China and in the Pea Ridge iron-oxide-apatite (IOA) deposit in Missouri. The coupling of compositional variations of monazite with thermodynamic modeling of fluid-rock interaction processes may provide a useful vectoring tool in these ore deposits. However, interpreting these geochemical fingerprints requires building an internally consistent thermodynamic dataset for REE minerals and their relevant aqueous complexes. In this study, a series of hydrothermal solubility experiments were carried out using synthetic monazite crystals (i.e., LaPO4, PrPO4, NdPO4, and EuPO4) to assess the consistency of reported mineral calorimetric data and the thermodynamic data of the aqueous REE complexes. The solubility experiments were conducted in aqueous HClO4-H3PO4–bearing solutions at temperatures between 100° and 250°C and at saturated water vapor pressure. Equilibrium constants (Ks0) for the dissolution reaction of monazite end members were retrieved as a function of temperature and extrapolated to standardmore »conditions of 25°C and 1 bar. Results indicate significant differences between the new solubility data and those reported in the literature. We demonstrate the impact of these new thermodynamic data in a series of fluid-rock interaction models using the GEMS code package ( and the MINES thermodynamic database ( The simulated monazite stability can be correlated to field observations and allows for the prediction of the behavior of REE in hydrothermal fluids and their association to alteration observed in ore deposits.« less
  5. Abstract Garnet is a common metamorphic and igneous mineral with extensive solid solution that can be stable to mantle depths ≥400 km. High-T and/or high-P garnet may contain oriented lamellae of other minerals, most commonly simple oxides (e.g., rutile, ilmenite), apatite, and, in ultrahigh-P cases, silicates including pyroxene and amphibole. Lamellae have classically been considered to be precipitation features preserving a record of former garnet chemistry richer in the lamellae nutrients (e.g., Ti4+). Such microtextural origins in precipitation systems (e.g., alloys) have long been studied via the crystallographic orientation relationships (COR) that form between a host and a separating phase, and by the shape-preferred orientation (SPO) of the lamellae. Recently, however, alternative hypotheses to precipitation have been suggested that require emplacement of lamellae in garnet by fluids, or co-growth, overgrowth, or inheritance mechanisms. These hypotheses posit that lamellae cannot be used to study former garnet chemistry. Moreover, they predict that lamellae phases, SPO, and COR should differ widely between localities, as lamellae formation will be controlled by various local rock-specific factors such as fluid presence, fluid chemistry, or mineral growth sequence. On the other hand, if lamellae characteristics are largely consistent between localities, it likely reflects control by precipitation energetics,more »rather than external factors. There have been few comparative COR studies in geologic systems, but the integrative assessment of COR, SPO, and lamellae assemblages should fingerprint lamellae growth process. To test the precipitation and alternative hypotheses, we collected large electron backscatter diffraction (EBSD) data sets for rutile, ilmenite, and apatite lamellae in garnet from the Brimfield Schist, Connecticut (≥1000 °C metamorphism; Central Maine Terrane, U.S.A.). We analyzed these data alongside published EBSD data for rutile, ilmenite, and corundum from metapegmatites metamorphosed in the eclogite facies from the Austrian Alps (Griffiths et al. 2016). The apatite data set is the first of its kind, and reveals that apatite preferentially aligns its close-packed direction parallel to that of garnet (c-axisapatite//<111>garnet). We also recognize a rutile-garnet COR related to those in meteorites with Widmanstätten patterns that are unequivocal products of exsolution. This is the first identification of direct similarities between silicate-oxide and metal-metal COR of which we are aware. Significantly, this rutile-garnet COR is found in diverse geologic settings including Connecticut and Idaho (U.S.A.), Austria, Germany, Greece, and China over a broad range of bulk-rock compositions. Results for all lamellae minerals show that COR are largely consistent between localities and, furthermore, are shared between apatite, ilmenite, and corundum. Moreover, between 70% and 95% of lamellae have COR and there is a dominant COR for each lamellae phase. Calculations show that d-spacing ratios of host-lamellae pairs can successfully predict the most commonly observed specific COR (those COR with two or more axial alignments with the host). These results, especially similarity of COR from markedly different geologic settings and a low diversity of lamellae minerals, are fully consistent with lamellae formation by precipitation (likely via exsolution). In contrast, the alternative hypotheses remain unsupported by COR results as well as by mineralogical and petrological evidence. Lamellae with similar traits as those in this work should thus be considered precipitates formed during unmixing of garnet compositions originally stable at elevated or extreme pressures and temperatures.« less