skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: VapoRock: Thermodynamics of Vaporized Silicate Melts for Modeling Volcanic Outgassing and Magma Ocean Atmospheres
Abstract Silicate vapors play a key role in planetary evolution, especially dominating early stages of rocky planet formation through outgassed magma ocean atmospheres. Our open-source thermodynamic modeling software “VapoRock” combines the MELTS liquid model with gas-species properties from multiple thermochemistry tables. VapoRock calculates the partial pressures of 34 gaseous species in equilibrium with magmatic liquid in the system Si–Mg–Fe–Al–Ca–Na–K–Ti–Cr–O at desired temperatures and oxygen fugacities (fO2, or partial pressure of O2). Comparison with experiments shows that pressures and melt-oxide activities (which vary over many orders of magnitude) are reproduced to within a factor of ∼3, consistent with measurement uncertainties. We also benchmark the model against a wide selection of igneous rock compositions including bulk silicate Earth, predicting elemental vapor abundances that are comparable to (Na, Ca, and Al) or more realistic than (K, Si, Mg, Fe, and Ti) those of the closed-source MAGMA code (with maximum deviations by factors of 10–300 for K and Si). Vapor abundances depend critically on the activities of liquid components. The MELTS model underpinning VapoRock was calibrated and extensively tested on natural igneous liquids. In contrast, MAGMA’s liquid model assumes ideal mixtures of a limited set of chemically simplified pseudospecies, which only roughly approximates the nonideal compositional interactions typical of many-component natural silicate melts. Finally, we explore how relative abundances of SiO and SiO2provide a spectroscopically measurable proxy for oxygen fugacity in devolatilized exoplanetary atmospheres, potentially constrainingfO2in outgassed exoplanetary mantles.  more » « less
Award ID(s):
1725025
PAR ID:
10408397
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
947
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 64
Size(s):
Article No. 64
Sponsoring Org:
National Science Foundation
More Like this
  1. Extensive experimental studies show that all major rock-forming elements (e.g., Si, Mg, Fe, Ca, Al, Na, K) dissolve in steam to a greater or lesser extent. We use these results to compute chemical equilibrium abundances of rocky-element-bearing gases in steam atmospheres equilibrated with silicate magma oceans. Rocky elements partition into steam atmospheres as volatile hydroxide gases (e.g., Si(OH)4, Mg(OH)2, Fe(OH)2, Ni(OH)2, Al(OH)3, Ca(OH)2, NaOH, KOH) and via reaction with HF and HCl as volatile halide gases (e.g., NaCl, KCl, CaFOH, CaClOH, FAl(OH)2) in much larger amounts than expected from their vapor pressures over volatile-free solid or molten rock at high temperatures expected for steam atmospheres on the early Earth and hot rocky exoplanets. We quantitatively compute the extent of fractional vaporization by defining gas/magma distribution coefficients and show that Earth's subsolar Si/Mg ratio may be due to loss of a primordial steam atmosphere. We conclude that hot rocky exoplanets that are undergoing or have undergone escape of steam-bearing atmospheres may experience fractional vaporization and loss of Si, Mg, Fe, Ni, Al, Ca, Na, and K. This loss can modify their bulk composition, density, heat balance, and interior structure. 
    more » « less
  2. Abstract High pressure and temperature experiments were carried out on the oxide mixtures corresponding to the bridgmanite stoichiometry under the hydrous shallow lower mantle conditions (24–25 GPa and 1673–1873 K with 5–10 wt. % of water in the starting material). Oxide mixtures investigated correspond to MgSiO3, (Mg, Fe)SiO3, (Mg, Al, Si)O3, and (Mg, Fe, Al, Si)O3. Melting was observed in all runs. Partitioning of various elements, including Mg, Fe, Si, and H is investigated. Melting under hydrous lower mantle conditions leads to increased (Mg + Fe)O/SiO2in the melt compared to the residual solids. The residual solids often contain a large amount of stishovite, and the melt contains higher (Mg,Fe)O/SiO2ratio than the initial material. (Mg + Fe)O‐rich hydrous melt could explain the low‐velocity anomalies observed in the shallow lower mantle and a large amount of stishovite in the residual solid may be responsible for the scattering of seismic waves in the mid‐lower mantle and may explain the “stishovite paradox. Since stishovite‐rich materials are formed only when silica‐rich source rock (MORB) is melted (not a typical peridotitic rock [bulk silicate Earth]), seismic scattering in the lower mantle provides a clue on the circulation of subducted MORB materials. To estimate hydrogen content, we use a new method of estimating the water content of unquenchable melts, and also propose a new interpretation of the significance of superhydrous phase B inclusions in bridgmanite. The results provide revised values of water partitioning between solid minerals and hydrous melts that are substantially higher than previous estimates. 
    more » « less
  3. Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO2, and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models. 
    more » « less
  4. ABSTRACT Observations of planetary material polluting the atmospheres of white dwarfs are an important probe of the bulk composition of exoplanetary material. Medium- and high-resolution optical and ultraviolet spectroscopy of seven white dwarfs with known circumstellar dust and gas emission are presented. Detections or meaningful upper limits for photospheric absorption lines are measured for: C, O, Na, S, P, Mg, Al, Si, Ca, Ti, Cr, Fe, and Ni. For 16 white dwarfs with known observable gaseous emission discs (and measured photospheric abundances), there is no evidence that their accretion rates differ, on average, from those without detectable gaseous emission. This suggests that, typically, accretion is not enhanced by gas drag. At the effective temperature range of the white dwarfs in this sample (16 000–25 000 K) the abundance ratios of elements are more consistent than absolute abundances when comparing abundances derived from spectroscopic white dwarf parameters versus photometric white dwarf parameters. Crucially, this highlights that the uncertainties on white dwarf parameters do not prevent white dwarfs from being utilized to study planetary composition. The abundances of oxygen and silicon for the three hydrogen-dominated white dwarfs in the sample with both optical and ultraviolet spectra differ by 0.62 dex depending on if they are derived from the optical or ultraviolet spectra. This optical/ultraviolet discrepancy may be related to differences in the atmospheric depth of line formation; further investigations into the white dwarf atmospheric modelling are needed to understand this discrepancy. 
    more » « less
  5. Abstract We present results from conducting a theoretical chemical analysis of a sample of benchmark companion brown dwarfs whose primary star is of type F, G, or K. We summarize the entire known sample of these types of companion systems, termed “compositional benchmarks,” that are present in the literature or recently published as key systems of study in order to best understand brown dwarf chemistry and condensate formation. Via mass balance and stoichiometric calculations, we predict a median brown dwarf atmospheric oxygen sink of 17.8 2.3 + 1.7 % by utilizing published stellar abundances in the local solar neighborhood. Additionally, we predict a silicate condensation sequence such that atmospheres with bulk Mg/Si ≲0.9 will form enstatite (MgSiO3) and quartz (SiO2) clouds, and atmospheres with bulk Mg/Si ≳0.9 will form enstatite and forsterite (Mg2SiO4) clouds. The implications of these results on C/O ratio trends in substellar-mass objects and the utility of these predictions in future modeling work are discussed. 
    more » « less