skip to main content


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
NSF-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

    For rocky exoplanets, knowledge of their geologic characteristics such as composition and mineralogy, surface recycling mechanisms, and volcanic behavior are key to determining their suitability to host life. Thus, determining exoplanet habitability requires an understanding of surface chemistry, and understanding the composition of exoplanet surfaces necessitates applying methods from the field of igneous petrology. Piston‐cylinder partial melting experiments were conducted on two hypothetical rocky exoplanet bulk silicate compositions. HEX1, a composition with molar Mg/Si = 1.42 (higher than bulk silicate Earth's Mg/Si = 1.23) yields a solidus similar to that of Earth's undepleted mantle. However, HEX2, a composition with molar Ca/Al = 1.07 (higher than Earth Ca/Al = 0.72) has a solidus with a slope of ∼10°C/kbar (vs. ∼15°C/kbar for Earth) and as result, has much lower melting temperatures than Earth. The majority of predicted adiabats point toward the likely formation of a silicate magma ocean for exoplanets with a mantle composition similar to HEX2. For adiabats that do intersect HEX2's solidus, decompression melting initiates at pressures more than 4x greater than in the modern Earth's undepleted mantle. The experimental partial melt compositions for these exoplanet mantle analogs are broadly similar to primitive terrestrial magmas but with higher CaO, and for the HEX2 composition, higher SiO2for a given degree of melting. This first of its kind exoplanetary experimental data can be used to calibrate future exoplanet petrologic models and predict volatile solubilities, volcanic degassing, and crust compositions for exoplanets with bulk compositions and ƒO2similar to those explored herein.

     
    more » « less
  3. 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
  4. 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
  5. Abstract

    Beckettite (Ca2V6Al6O20; IMA 2015‐001) is a newly discovered refractory mineral, occurring as micrometer‐sized grains intergrown with hibonite and perovskite, and surrounded by secondary grossular, anorthite, coulsonite, hercynite, and corundum. It occurs within highly altered areas in a V‐rich, Type A Ca‐Al‐rich inclusion (CAI), A‐WP1, from the Allende CV3 carbonaceous chondrite. The type beckettite has an empirical formula of (Ca1.99Na0.01)(V3+3.47Al1.40Ti4+0.57Mg0.25Sc0.08Fe2+0.04)(Al5.72Si0.28)O20, with a triclinic structure in space groupPand cell parametersa= 10.367 Å,b= 10.756 Å,c= 8.895 Å, α = 106.0°, β = 96.0°, γ = 124.7°,V= 739.7 Å3, andZ= 2, which leads to a calculated density of 3.67 g cm−3. Beckettite’s general formula is Ca2(V,Al,Ti,Mg)6Al6O20and the endmember formula is Ca2V6Al6O20. Beckettite is slightly16O‐depleted (Δ17O = −16 ± 2‰) compared to the coexisting hibonite and spinel −24 ± 2‰. Beckettite is a primary high‐temperature mineral resulting from igneous crystallization of an16O‐rich V‐rich CAI melt together with V‐bearing hibonite, perovskite, burnettite, spinel, and paqueite. Subsequently, beckettite experienced an incomplete isotope exchange with an16O‐poor aqueous fluid (Δ17O = −3 ± 2‰) on the Allende parent asteroid.

     
    more » « less