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 (
- Award ID(s):
- 1412175
- NSF-PAR ID:
- 10036394
- Date Published:
- Journal Name:
- Astrophysical journal
- Volume:
- 824
- ISSN:
- 1538-4357
- Page Range / eLocation ID:
- id 103
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract f O2, 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 constrainingf O2in outgassed exoplanetary mantles. -
more » « less
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.
-
Abundances of minor and trace elements in olivine are increasingly used as petrogenetic indicators for mantle source lithologies, mantle metasomatism history, mantle potential temperatures, and magmatic differentiation. As it is common for olivine to be complexly zoned on a fine-scale, high precision analytical methods for EPMA (electron microprobe microanalyzer, or Electron Microprobe) trace element analysis under high spatial resolution have been developed. However, previous studies have focused more on analytical precision with fewer efforts in examining the accuracy of the data. In this study, we used the Cameca SXFive field emission (FE) EPMA to fully evaluate the effects of beam settings, background offsets and background regression models, and primary calibration standards on the data accuracy of 10 key petrogenetic elements (Na, Al, P, Ca, Ti, Cr, Mn, Co, Ni, and Zn) using MongOlSh11–2 olivine as a reference material. Our results indicate that high voltage, high beam current and long counting time not only improve data precision, but also improve data accuracy, especially on elements with low P/B (peak/background) ratios such as Zn and Cr. Importantly, careful background offsets and background regression models need to be obtained via high resolution WDS relative scans or step scans on each target element. Special care needs to be paid to Co element analysis to avoid or correct for peak interference of Fe Kβ. Among 10 minor and trace elements, exponential background regression models need to be applied to Al, Mn, and Ti elements, whereas other elements require linear background regression. Furthermore, to avoid Al and Zn surface contamination due to alumina polishing or brass presence, ultrasonic cleaning between each intermediate polishing steps and plasma cleaning immediately prior to EPMA experiments is highly recommended. Micro-inclusions such as chromite and spinel in olivine or adjacent Ca-rich phases need to be avoided to minimize primary or secondary fluorescence-related contamination on Al, Cr, or Ca. As a volatile element, Na element needs to be analyzed first with appropriate counting time to minimize the Na loss under high beam conditions. It needs mentioning that major elements (Mg, Fe, and Si) are best analyzed using MongOlSh11–2 or San Carlos olivine as primary standards for calibrations, which can yield more accurate data for both major elements and trace elements because of the improved matrix- corrections. Using our recommended analytical protocols, we have successfully discriminated “depleted” mantle olivine cores from an EMORB in northern East Pacific Rise (EPR) via Ca, Ti, Ni, Co, and Mn abundances. Our olivine data from Siqueiros Transform and the nearby 8◦20′ N seamounts also help reveal a metasomatized peridotite mantle beneath the northern EPR. Overall, the protocols proposed in this study can serve as a guide for accurate EPMA olivine trace element analyses, which potentially contributes to the efforts of fostering a comparable olivine database worldwide.more » « less
-
We provide estimates of atmospheric pressure and surface composition on short-period, rocky exoplanets with dayside magma pools and silicate-vapor atmospheres. Atmospheric pressure tends toward vapor-pressure equilibrium with surface magma, and magma-surface composition is set by the competing effects of fractional vaporization and surface-interior exchange. We use basic models to show how surface-interior exchange is controlled by the planet’s temperature, mass, and initial composition. We assume that mantle rock undergoes bulk melting to form the magma pool, and that winds flow radially away from the substellar point. With these assumptions, we find that: (1) atmosphere-interior exchange is fast when the planet’s bulk-silicate FeO concentration is low, and slow when the planet’s bulk-silicate FeO concentration is high; (2) magma pools are compositionally well mixed for substellar temperatures ≲2400 K, but compositionally variegated and rapidly variable for substellar temperatures ≳2400 K; (3) currents within the magma pool tend to cool the top of the solid mantle (“tectonic refrigeration”) (4) contrary to earlier work, many magma planets have time-variable surface compositions.more » « less
-
more » « less
Abstract The transition between blueschist and eclogite plays an important role in subduction zones via dehydration and densification processes in descending oceanic slabs. There are a number of previous petrological studies describing potential mineral reactions taking place at the transition. An experimental determination of such reactions could help constrain the pressure–temperature conditions of the transition as well as the processes of dehydration. However, previous experimental contributions have focused on the stability of spontaneously formed hydrous minerals in basaltic compositions rather than on reactions among already formed blueschist facies minerals. Therefore, this study conducted three groups of experiments to explore the metamorphic reactions among blueschist facies minerals at conditions corresponding to warm subduction, where faster reaction rates are possible on the time scale of laboratory experiments. The first group of experiments was to establish experimental reversals of the reaction glaucophane+paragonite to jadeite+pyrope+quartz+H2O over the range of 2.2–3.5 GPa and 650–820°C. This reaction has long been treated as key to the blueschist–eclogite transition. However, only the growth of glaucophane+paragonite was observed at the intersectional stability field of both paragonite and jadeite+quartz, confirming thermodynamic calculations that the reaction is not stable in the system Na2O–MgO–Al2O3–SiO2–H2O. The second set of experiments involved unreversed experiments using glaucophane+zoisite ±quartz in low‐Fe and Ca‐rich systems and were run at 1.8–2.4 GPa and 600–780°C. These produced omphacite+paragonite/kyanite+H2O accompanied by compositional shifts in the sodium amphibole, glaucophane, towards sodium–calcium amphiboles such as winchite (☐(CaNa)(Mg4Al)Si8O22(OH)2) and barroisite (☐(CaNa)(Mg3Al2)(AlSi7)O22(OH)2). This suggests that a two‐step dehydration occurs, first involving the breakdown of glaucophane+zoisite towards a paragonite‐bearing assemblage, then the breakdown of paragonite to release H2O. It also indicates that sodium–calcium amphibole can coexist with eclogite phases, thereby extending the thermal stability of amphibole to greater subduction zone depths. The third set of experiments was an experimental investigation at 2.0–2.4 GPa and 630–850°C involving a high‐Fe (Fe#=Fetotal/(Fetotal+Mg)≈0.36) natural glaucophane, synthetic paragonite and their eclogite‐forming reaction products. The results indicated that garnet and omphacite grew over most of these pressure–temperature conditions, which demonstrates the importance of Fe‐rich glaucophane in forming the key eclogite assemblage of garnet+omphacite, even under warm subduction zone conditions. Based on the experiments of this study, reaction between glaucophane+zoisite is instrumental in controlling dehydration processes at the blueschist–eclogite transition during warm subduction.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/0004-637X/824/2/103
