More Like this

1. Experimental Study of Halogen Partitioning Between Amphibole and Melt in Subduction Magmas

   Dawson, H.G.; Krawczynski, M.J. (December 2022, AGU Fall meeting)

   Subduction of oceanic plates releases large amounts of chlorine and fluorine into the mantle. These halogens are transported into the crust through hydrous melting, where they may be incorporated into minerals such as biotite, apatite, and amphibole. Halogen concentrations are measured in volcanic or plutonic material, while the concentration of Cl and F released from the subducting slab must be calculated based upon the amount of crystallized material and the partition coefficients of each mineral. As amphibole is the most common halogen bearing igneous mineral, it is commonly studied as a bearer of Cl and F. However, the partition coefficient of F between amphibole and a hydrous melt has not been agreed upon by previous studies. Here we show that F is moderately to highly compatible in amphibole, in agreement with other experiments performed at crustal conditions. As amphibole may be able to incorporate a large amount of F, cryptic amphibole crystallization may raise the Cl/F ratio of residual magma, which will then be transported to the surface bearing this geochemical signature, even with little crystallized amphibole present in erupted material. This provides further evidence for the occurrence of cryptic amphibole crystallization, previously predicted based on REE studies and phase equilibria. A better understanding of the halogen reservoirs present in the crust will allow for more accurate estimates of the amount of Cl and F released by subducting slabs. 

   more » « less
2. The Geochemistry of Magnetite and Apatite from the El Laco Iron Oxide-Apatite Deposit, Chile: Implications for Ore Genesis

   La Cruz, N. (November 2020, Economic geology)

   null (Ed.)

   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, and 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. 

   more » « less
3. An ab-initio study on the thermodynamics of disulfide, sulfide, and bisulfide incorporation into apatite and the development of a more comprehensive temperature, pressure, pH, and composition-dependent model for ionic substitution in minerals

   https://doi.org/10.2138/am-2022-8250  

   Kim, YoungJae; Konecke, Brian; Simon, Adam; Fiege, Adrian; Becker, Udo (November 2022, American Mineralogist)

   Abstract 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 S 6+/Σ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, the 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 is 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 is also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS(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 to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures. 

   more » « less
4. An ab-initio study on the thermodynamics of disulfide, sulfide, and bisulfide incorporation into apatite and the development of a more comprehensive temperature, pressure, pH, and composition-dependent model for ionic substitution in minerals

   Young Jae Kim, Brian Konecke (September 2022, The American mineralogist)

   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, the 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. 

   more » « less
5. Constraining Halogen Budgets in Continental Arcs at the Whole-Rock Scale

   Bucheli, César; Segee-Wright, George; Grabiec, Justine G; Cooperdock, Emily HG; Barnes, Jaime D; Lackey, Jade Star; Stockli, Lisa D; Stockli, Daniel F; Ardill, Katie (June 2024, 2024 V.M. Goldschmidt Conference Abstracts)

   Halogens (F, Cl, Br, I) are primary components of volcanic gas emissions and play an essential role in continental arc magmatic environments due to their solubility in fluids that generate metallic ore deposits. Despite their ubiquity, the behavior and budget of halogens in continental arc environments are poorly constrained. We investigated the plutonic and volcanic halogen budgets in intermediate-to-felsic igneous rocks (56–77 wt% SiO2) from the Sierra Nevada (California) - a Mesozoic continental arc where plutonic and volcanic outcrops can be correlated via their geographic, compositional, and geochronologic framework. We measured the halogen concentrations of bulk rock powders and their leachates via ion chromatography (F, Cl) and ICP-MS (Br, I). Halogen concentrations in our rock powders range between 107–727 μg/g F, 13–316 μg/g Cl, 2–323 ng/g Br, and 1–69 ng/g I. In contrast, leachates yielded 3–4 orders of magnitude less Cl and F, one order of magnitude less I, and similar amounts of Br compared to their corresponding bulk rocks. Preliminary data show no significant differences between volcanic and plutonic samples, suggesting that halogen concentrations in these rocks are insensitive to shallow fractionation. Although F and I exhibit no correlation with major element compositions, Cl and Br display negative trends with increasing SiO2 and K2O, and positive trends with increasing Fe2O3T, MnO, MgO, CaO, and TiO2, suggesting mafic minerals as important hosts of structurally bound halogens. Overall, Sierran plutonic rocks display low halogen contents (max. F, Cl = 727, 315 μg/g), consistent with biotite- and apatite-bearing granitoids reported in [1]. This work suggests that halogens do not preferentially enrich in shallow plutonic or volcanic portions of a continental arc system and that mafic mineral phases likely serve as primary reservoirs of these elements in intermediate-to-felsic igneous rocks. These hypotheses will be further investigated in future work through in-situ analysis of halogen concentrations in crystals. [1] Teiber, Marks, Wenzel, Siebel, Altherr & Markl (2014), Chemical Geology, vol. 374–375, pp. 92–109, doi: 10.1016/j.chemgeo.2014.03.006 

   more » « less