skip to main content

Title: A comparison between the stability fields of a Cl-rich scapolite and the end-member marialite
Abstract Scapolites are pervasive rock-forming aluminosilicates that are found in metamorphic, igneous, and hydrothermal environments; nonetheless, the stability field of Cl-rich scapolite is not well constrained. This experimental study investigated two reactions involving Cl-rich scapolite. First, the anhydrous reaction 1 of plagioclase + halite + calcite to form scapolite [modeled as: 3 plagioclase (Ab80An20) + 0.8 NaCl + 0.2 CaCO3 = scapolite (Ma80Me20)] was investigated to determine the effect of the Ca-rich meionite (Me = Ca4Al6Si6O24CO3) component on the Na end-member marialite (Ma = Na4Al3Si9O24Cl). Second, the effect of water on this reaction was investigated using the hydrothermally equivalent reaction 2, H2O + scapolite (Ma80Me20) = 3 plagioclase (Ab80An20) + CaCO3 + liquid, where the liquid is assumed to be a saline-rich hydrous-silicate melt. Experiments were conducted with synthetic phases over the range of 500–1030 °C and 0.4–2.0 GPa. For reaction 1, intermediate composition scapolite shows a wide thermal stability and is stable relative to plagioclase + halite + calcite at temperatures above 750 °C at 0.4 GPa and 760 °C at 2.0 GPa. For reaction 2, intermediate scapolite appears to be quite tolerant of water; it forms at a minimum bulk salinity [XNaCl = molar ratio of NaCl/(NaCl+H2O)] of more » the brine of approximately 0.2 XNaCl at 830 and 680 °C at pressures of 2.0 and 1.5 GPa, respectively. Based on the study done by Almeida and Jenkins (2017), pure marialite is very intolerant of water when compared to intermediate composition scapolite. Compositional changes in the scapolite and plagioclase were characterized by X-ray diffraction and electron microprobe analysis and found to shift from the nominal bulk compositions to the observed compositions of Ma85Me15 for scapolite and to Ab91An09 for plagioclase. These results were used to model the phase equilibria along the marialitemeionite join in temperature-composition space. This study demonstrates that a small change in the scapolite composition from end-member marialite to Ma85Me15 expands the stability field of marialite significantly, presumably due to the high entropy of mixing in scapolite, as well as increases its tolerance to water. This supports the much more common presence of intermediate scapolites in hydrothermal settings than either end-member meionite or marialite as is widely reported in the literature. « less
Award ID(s):
Publication Date:
Journal Name:
American Mineralogist
Page Range or eLocation-ID:
1788 to 1799
Sponsoring Org:
National Science Foundation
More Like this
  1. Iron-rich phyllosilicates on Mars comprise nearly 90% of the H2O- and OH-bearing phases observed directly by rovers and remotely by orbiters (Chemtob et al., 2017, JGR). Theories concerning the possible origin of Fe-rich smectite during Mars’ earliest history (phyllosian) are hard to test because of limited knowledge of the upper-thermal stability of Fe-rich phyllosilicates. In this study we present data on the upper-thermal stability of a pure-iron smectite to put some minimum constraints on its possible high-temperature origin early in Mars history either from a primordial atmosphere or by hydrothermal activity. Smectite coexisting with quartz and magnetite was synthesized from the oxides in the system Na2O-FeO-Fe2O3-Al2O3-SiO2-H2O at 500°C and 2 kbar and fO2 near FMQ. Reversal experiments involved mixtures with equal portions of the smectite-synthesis and breakdown products (quartz, fayalite, albite, magnetite (mt) treated in the presence of about 10 wt% H2O over the range of 1-3 kbar and 530-640°C. The average composition (electron microprobe) of smectite formed both in synthesis and in reversal experiments was Na0.35(Fe2+2.28Fe3+0.31Al0.41)(Al1.07Si2.93)O10(OH)2·nH2O, where ferric iron was calculated by summing the octahedral cations to 3.0. Reversals for the reaction smec+mt1 = fayalite+albite+mt2+quartz+H2O were obtained at 538±8, 590±10, and 610±10°C at 1, 2, and 3 kbar, respectively,more »where mt1 and mt2 have the approximate compositions Fe2.8Si0.2O4 and Fe2.8Al0.1O4, respectively, with all other phases being pure. This smectite has up to 2 interlayer H2O at 25°C (and high humidity), losing 1 H2O at <50°C, and the second at 125 ± 25°C. Thermodynamic modeling of this reaction was used to extrapolate the upper-thermal stability of this Fe-smectite down to 10 bars and approximately 239°C. Applications of these results indicate the maximum temperature for forming Fe-smectite from a dense primordial atmosphere of 100 bars is 390 ± 25°C. Crustal storage of water in Fe-smectite ranges up to a maximum of 10.7 wt% at ~2 km and 40°C, 7.4 wt% at 6 km and 120°C, and 3.8 wt% H2O at 32 km and 625°C for a Noachian geotherm of 20°C/km. Plain language summary: This study presents experimental limits on the temperatures at which the clay mineral smectite might form on Mars, either from a dense primordial atmosphere (390°C at 100 bars) or by high-temperature hydrothermal activity (625°C at 32 km). Because this study deals with iron end-member clay, these are minimum temperatures; any solid solution with magnesium will increase these temperatures.« less
  2. Teagle, Damon A (Ed.)
    The Cedars ultramafic block hosts alkaline springs (pH > 11) in which calcium carbonate forms upon uptake of atmospheric CO2 and at times via mixing with surface water. These processes lead to distinct carbonate morphologies with ‘‘floes” forming at the atmosphere-water interface, ‘‘snow” of fine particles accumulating at the bottom of pools and terraced constructions of travertine. Floe material is mainly composed of aragonite needles despite CaCO3 precipitation occurring in waters with low Mg/Ca (<0.01). Precipitation of aragonite is likely promoted by the high pH (11.5–12.0) of pool waters, in agreement with published experiments illustrating the effect of pH on calcium carbonate polymorph selection. The calcium carbonates exhibit an extreme range and approximately 1:1 covariation in d13C (9 to 28‰ VPDB) and d18O (0 to 20‰ VPDB) that is characteristic of travertine formed in high pH waters. The large isotopic fractionations have previously been attributed to kinetic isotope effects accompanying CO2 hydroxylation but the controls on the d13C-d18O endmembers and slope have not been fully resolved, limiting the use of travertine as a paleoenvironmental archive. The limited areal extent of the springs (0.5 km2) and the limited range of water sources and temperatures, combined with our sampling strategy, allow usmore »to place tight constraints on the processes involved in generating the systematic C and O isotope variations. We develop an isotopic reaction–diffusion model and an isotopic box model for a CO2-fed solution that tracks the isotopic composition of each dissolved inorganic carbon (DIC) species and CaCO3. The box model includes four sources or sinks of DIC (atmospheric CO2, high pH spring water, fresh creek water, and CaCO3 precipitation). Model parameters are informed by new floe D44Ca data (0.75 ± 0.07‰), direct mineral growth rate measurements (4.8 to 8  107 mol/m2/s) and by previously published elemental and isotopic data of local water and DIC sources. Model results suggest two processes control the extremes of the array: (1) the isotopically light end member is controlled by the isotopic composition of atmospheric CO2 and the kinetic isotope fractionation factor (KFF (‰) = (a  1)  1000) accompanying CO2 hydroxylation, estimated here to be 17.1 ± 0.8‰ (vs. CO2(aq)) for carbon and 7.1 ± 1.1‰ (vs. ‘CO2(aq)+H2O’) for oxygen at 17.4 ± 1.0 C. Combining our results with revised CO2 hydroxylation KFF values based on previous work suggests consistent KFF values of 17.0 ± 0.3‰ (vs. CO2(aq)) for carbon and 6.8 ± 0.8‰ for oxygen (vs. ‘CO2(aq)+H2O’) over the 17–28 C temperature range. (2) The isotopically heavy endmember of calcium carbonates at The Cedars reflects the composition of isotopically equilibrated DIC from creek or surface water (mostly HCO- 3, pH = 7.8–8.7) that occasionally mixes with the high-pH spring water. The bulk carbonate d13C and d18O values of modern and ancient travertines therefore reflect the proportion of calcium carbonate formed by processes (1) and (2), with process (2) dominating the carbonate precipitation budget at The Cedars. These results show that recent advances in understanding kinetic isotope effects allow us to model complicated but common natural processes, and suggest ancient travertine may be used to retrieve past meteoric water d18O and atmospheric d13C values. There is evidence that older travertine at The Cedars recorded atmospheric d13C that predates large-scale combustion of fossil fuels.« less
  3. 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
  4. Abstract Volcán Quizapu, Chile, is an under-monitored volcano that was the site of two historical eruptions: an effusive eruption in 1846–1847 and a Plinian eruption in 1932, both of which discharged ∼5 km3 (dense rock equivalent) of lava and/or tephra. The majority of material erupted in both cases is trachydacite, nearly identical for each event. We present H2O-saturated, phase equilibrium experiments on this end-member dacite magma, using a pumice sample from the 1932 eruption as the main starting material. At an oxygen fugacity (fO2) of ∼NNO + 0·2 (where NNO is the nickel–nickel oxide buffer), the phase assemblage of An25–30 plagioclase + amphibole + orthopyroxene, without biotite, is stable at 865 ± 10 °C and 110 ± 20 MPa H2O pressure (PH2O), corresponding to ∼4 km depth. At these conditions, experiments also reproduce the quenched glass composition of the starting pumice. At slightly higher PH2O and below 860 °C, biotite joins the equilibrium assemblage. Because biotite is not part of the observed Quizapu phase assemblage, its presence places an upper limit on PH2O. At the determined storage PH2O of ∼110 MPa, H2O undersaturation of the magma with XH2Ofluid = 0·87 would align Ptotal to mineral-based geobarometry estimates of ∼130 MPa. However, XH2Ofluid < 1 is not required to reproduce the Quizapu dacite phasemore »assemblage and compositions. A second suite of experiments at lower fO2 shows that the stability fields of the hydrous silicates (amphibole and biotite) are significantly restricted at NNO – 2 relative to NNO + 0·2. Additional observations of Quizapu lava and pumice samples support the existing hypothesis that rapid pre-eruptive heating drove the effusive 1846–1847 eruption, with important refinements. We demonstrate that microlites in the end-member dacite lavas are consistent with in situ crystallization (during ascent), rather than transfer from an andesite. In one end-member dacite lava, newly identified reverse zoning in orthopyroxene and incipient destabilization of amphibole are consistent with small degrees of heating. Our work articulates a clear direction for future Quizapu studies, which are warranted given the active nature of the Cerro Azul–Descabezado Grande volcanic axis.« less
  5. ABSTRACT Formation of the feldspathoid sodalite (Na6Al6Si6O24·2NaCl) by reaction of nepheline (NaAlSiO4) with NaCl-bearing brines was investigated at 3 and 6 kbar and at a constant temperature of 750 °C to determine the brine concentration at which sodalite forms with variation in pressure. The reaction boundary was located by reaction-reversal experiments in the system NaAlSiO4–NaCl–H2O at a brine concentration of 0.16 ± 0.08 XNaCl [= molar ratio NaCl/(NaCl + H2O)] at 3 kbar and at a brine concentration of 0.35 ± 0.03 XNaCl at 6 kbar. Characterization of the sodalite using both X-ray diffraction and infrared spectroscopy after treatment in these brines indicated no obvious evidence of water or hydroxyl incorporation into the cage structure of sodalite. The data from this study were combined with earlier results by Wellman (1970) and Sharp et al. (1989) at lower (1–1.5 kbar) and higher (7–8 kbar) pressures, respectively, on sodalite formation from nepheline and NaCl which models as a concave-down curve in XNaCl – P space. In general, sodalite buffers the concentration of neutral aqueous NaCl° in the brine to relatively low values at P < 4 kbar, but NaCl° increases rapidly at higher pressures. Thermochemical modeling of these data was done tomore »determine the activity of the aqueous NaCl° relative to a 1 molal (m) standard state, demonstrating very low activities (<0.2 m, or 1.2 wt.%) of NaCl° at 3 kbar and lower, but rising to relatively high activities (>20 m, or 54 wt.%) of NaCl° at 6 kbar or higher. The results from this study place constraints on the concentration of NaCl° in brines coexisting with nepheline and sodalite and, because of the relative insensitivity of this reaction to temperature, can provide a convenient geobarometer for those localities where the fluid compositions that formed nepheline and sodalite can be determined independently.« less