-
Distinguishing oceans of water from magma on mini-Neptune K2-18b
Authors:
Oliver Shorttle,
Sean Jordan,
Harrison Nicholls,
Tim Lichtenberg,
Dan J. Bower
Abstract:
Mildly irradiated mini-Neptunes have densities potentially consistent with them hosting substantial liquid water oceans (`Hycean' planets). The presence of CO2 and simultaneous absence of ammonia (NH3) in their atmospheres has been proposed as a fingerprint of such worlds. JWST observations of K2-18b, the archetypal Hycean, have found the presence of CO2 and the depletion of NH3 to <100 ppm; hence…
▽ More
Mildly irradiated mini-Neptunes have densities potentially consistent with them hosting substantial liquid water oceans (`Hycean' planets). The presence of CO2 and simultaneous absence of ammonia (NH3) in their atmospheres has been proposed as a fingerprint of such worlds. JWST observations of K2-18b, the archetypal Hycean, have found the presence of CO2 and the depletion of NH3 to <100 ppm; hence, it has been inferred that this planet may host liquid water oceans. In contrast, climate modelling suggests that many of these mini-Neptunes, including K2-18b, may likely be too hot to host liquid water. We propose a solution to this discrepancy between observation and climate modelling by investigating the effect of a magma ocean on the atmospheric chemistry of mini-Neptunes. We demonstrate that atmospheric NH3 depletion is a natural consequence of the high solubility of nitrogen species in magma at reducing conditions; precisely the conditions prevailing where a thick hydrogen envelope is in communication with a molten planetary surface. The magma ocean model reproduces the present JWST spectrum of K2-18b to < 3 sigma, suggesting this is as credible an explanation for current observations as the planet hosting a liquid water ocean. Spectral areas that could be used to rule out the magma ocean model include the >4um region, where CO2 and CO features dominate: Magma ocean models suggest a systematically lower CO2/CO ratio than estimated from free chemistry retrieval, indicating that deeper observations of this spectral region may be able to distinguish between oceans of liquid water and magma on mini-Neptunes.
△ Less
Submitted 21 February, 2024; v1 submitted 11 January, 2024;
originally announced January 2024.
-
Interior dynamics of super-Earth 55 Cancri e
Authors:
Tobias G. Meier,
Dan J. Bower,
Tim Lichtenberg,
Mark Hammond,
Paul J. Tackley
Abstract:
The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot…
▽ More
The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot enough to harbour a magma ocean. We use results from general circulation models (GCMs) of atmospheres to constrain the surface temperature contrasts. There is still a large uncertainty on the vigour and style of mantle convection in super-Earths, especially those that experience stellar irradiation large enough to harbour a magma ocean. In this work, we aim to constrain the mantle dynamics of the tidally locked lava world 55 Cancri e. Using the surface temperature contrasts as boundary condition, we model the mantle flow of 55 Cancri e using 2D mantle convection simulations and investigate how the convection regimes are affected by the different climate models. We find that large super-plumes form on the dayside if that hemisphere is covered by a magma ocean and the nightside remains solid or only partially molten. Cold material descends into the deep interior on the nightside, but no strong downwellings form. In some cases, the super-plume also moves several tens of degrees towards the terminator. A convective regime where the upwelling is preferentially on the dayside might lead to preferential outgassing on that hemisphere which could lead to the build-up of atmospheric species that could be chemically distinct from the nightside.
△ Less
Submitted 1 August, 2023;
originally announced August 2023.
-
Investigating the visible phase-curve variability of 55 Cnc e
Authors:
E. A. Meier Valdés,
B. M. Morris,
B. -O. Demory,
A. Brandeker,
D. Kitzmann,
W. Benz,
A. Deline,
H. -G. Florén,
S. G. Sousa,
V. Bourrier,
V. Singh,
K. Heng,
A. Strugarek,
D. J. Bower,
N. Jäggi,
L. Carone,
M. Lendl,
K. Jones,
A. V. Oza,
O. D. S. Demangeon,
Y. Alibert,
R. Alonso,
G. Anglada,
J. Asquier,
T. Bárczy
, et al. (65 additional authors not shown)
Abstract:
55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-cur…
▽ More
55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-curve modulation in 55 Cnc e. To this end, we used the CHaracterising ExOPlanet Satellite (CHEOPS), whose exquisite photometric precision provides an opportunity to characterise minute changes in the phase curve from one orbit to the next. CHEOPS observed 29 individual visits of 55 Cnc e between March 2020 and February 2022. Based on these observations, we investigated the different processes that could be at the origin of the observed modulation. In particular, we built a toy model to assess whether a circumstellar torus of dust driven by radiation pressure and gravity might match the observed flux variability timescale. We find that the phase-curve amplitude and peak offset of 55 Cnc e do vary between visits. The sublimation timescales of selected dust species reveal that silicates expected in an Earth-like mantle would not survive long enough to explain the observed phase-curve modulation. We find that silicon carbide, quartz, and graphite are plausible candidates for the circumstellar torus composition because their sublimation timescales are long. The extensive CHEOPS observations confirm that the phase-curve amplitude and offset vary in time.We find that dust could provide the grey opacity source required to match the observations. However, the data at hand do not provide evidence that circumstellar material with a variable grain mass per unit area causes the observed variability. Future observations with the James Webb Space Telescope promise exciting insights into this iconic super-Earth.
△ Less
Submitted 27 July, 2023; v1 submitted 12 July, 2023;
originally announced July 2023.
-
Diverse Carbonates in Exoplanet Oceans Promote the Carbon Cycle
Authors:
Kaustubh Hakim,
Meng Tian,
Dan J. Bower,
Kevin Heng
Abstract:
Carbonate precipitation in oceans is essential for the carbonate-silicate cycle (inorganic carbon cycle) to maintain temperate climates. By considering the thermodynamics of carbonate chemistry, we demonstrate that the ocean pH decreases by approximately 0.5 for a factor of 10 increase in the atmospheric carbon dioxide content. The upper and lower limits of ocean pH are within 1-4 of each other, w…
▽ More
Carbonate precipitation in oceans is essential for the carbonate-silicate cycle (inorganic carbon cycle) to maintain temperate climates. By considering the thermodynamics of carbonate chemistry, we demonstrate that the ocean pH decreases by approximately 0.5 for a factor of 10 increase in the atmospheric carbon dioxide content. The upper and lower limits of ocean pH are within 1-4 of each other, where the upper limit is buffered by carbonate precipitation and defines the ocean pH when the carbon cycle operates. If the carbonate compensation depth (CCD) resides above the ocean floor, then carbonate precipitation and the carbon cycle cease to operate. The CCD is deep (>40 km) for high ocean temperature and high atmospheric carbon dioxide content. Key divalent carbonates of magnesium, calcium and iron produce an increasingly wider parameter space of deep CCDs, suggesting that chemical diversity promotes the carbon cycle. The search for life from exoplanets will benefit by including chemically more diverse targets than Earth twins.
△ Less
Submitted 9 January, 2023; v1 submitted 6 January, 2023;
originally announced January 2023.
-
Solubility of water in peridotite liquids and the prevalence of steam atmospheres on rocky planets
Authors:
Paolo A. Sossi,
Peter M. E. Tollan,
James Badro,
Dan J. Bower
Abstract:
Atmospheres are products of time-integrated mass exchange between the surface of a planet and its interior. On Earth and other planetary bodies, magma oceans likely marked significant atmosphere-forming events, during which both steam- and carbon-rich atmospheres may have been generated. However, the nature of Earth's early atmosphere, and those around other rocky planets, remains unclear for lack…
▽ More
Atmospheres are products of time-integrated mass exchange between the surface of a planet and its interior. On Earth and other planetary bodies, magma oceans likely marked significant atmosphere-forming events, during which both steam- and carbon-rich atmospheres may have been generated. However, the nature of Earth's early atmosphere, and those around other rocky planets, remains unclear for lack of constraints on the solubility of water in liquids of appropriate composition. Here we determine water solubility in 14 peridotite liquids, representative of Earth's mantle, synthesised in a laser-heated aerodynamic levitation furnace. We explore oxygen fugacities (fO$_2$) between -1.9 and +6.0 log units relative to the iron-wüstite buffer at constant temperature (2173$\pm$50 K) and total pressure (1 bar). The resulting fH$_2$O ranged from 0 to 0.027 bar and fH$_2$ from 0 to 0.064 bar. Total H$_2$O contents were determined by transmission FTIR spectroscopy from the absorption band at 3550cm$^{-1}$ and applying the Beer-Lambert law. The mole fraction of water in the liquid is $\propto$ (fH$_2$O)$^{0.5}$, attesting to its dissolution as OH. The data are fit by a solubility coefficient of 524$\pm$16 ppmw/bar$^{0.5}$, for a molar absorption coefficient, $ε_{3550}$, of 6.3$\pm$0.3 m$^2$/mol in basaltic glasses or 647$\pm$25 ppmw/bar$^{0.5}$, with $ε_{3550}$ = 5.1$\pm$0.3m$^2$/mol for peridotitic glasses. These solubility constants are 10-25 % lower than those for basaltic liquids at 1623 K and 1 bar. Higher temperature lowers water solubility, offsetting the greater depolymerisation of peridotite melts that would otherwise increase H$_2$O solubility relative to basaltic liquids. Because the solubility of water remains high relative to that of CO$_2$, steam atmospheres are rare, although they may form under oxidising conditions on telluric bodies, provided high H/C ratios prevail.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
Interior heating of rocky exoplanets from stellar flares with application to TRAPPIST-1
Authors:
Alexander Grayver,
Dan J. Bower,
Joachim Saur,
Caroline Dorn,
Brett M. Morris
Abstract:
Many stars of different spectral types with planets in the habitable zone are known to emit flares. Until now, studies that address the long-term impact of stellar flares and associated Coronal Mass Ejections (CMEs) assumed that the planet's interior remains unaffected by interplanetary CMEs, only considering the effect of plasma/UV interactions on the atmosphere of planets. Here, we show that the…
▽ More
Many stars of different spectral types with planets in the habitable zone are known to emit flares. Until now, studies that address the long-term impact of stellar flares and associated Coronal Mass Ejections (CMEs) assumed that the planet's interior remains unaffected by interplanetary CMEs, only considering the effect of plasma/UV interactions on the atmosphere of planets. Here, we show that the magnetic flux carried by flare-associated CMEs results in planetary interior heating by Ohmic dissipation and leads to a variety of interior--exterior interactions. We construct a physical model to study this effect and apply it to the TRAPPIST-1 star whose flaring activity has been constrained by Kepler observations. Our model is posed in a stochastic manner to account for uncertainty and variability in input parameters. Particularly for the innermost planets, our results suggest that the heat dissipated in the silicate mantle is both of sufficient magnitude and longevity to drive geological processes and hence facilitate volcanism and outgassing of the TRAPPIST-1 planets. Furthermore, our model predicts that Joule heating can further be enhanced for planets with an intrinsic magnetic field compared to those without. The associated volcanism and outgassing may continuously replenish the atmosphere and thereby mitigate the erosion of the atmosphere caused by the direct impact of flares and CMEs. To maintain consistency of atmospheric and geophysical models, the impact of stellar flares and CMEs on atmospheres of close-in exoplanetary systems needs to be studied in conjunction with the effect on planetary interiors.
△ Less
Submitted 8 December, 2022; v1 submitted 11 November, 2022;
originally announced November 2022.
-
VapoRock: Thermodynamics of vaporized silicate melts for modeling volcanic outgassing and magma ocean atmospheres
Authors:
Aaron S. Wolf,
Noah Jäggi,
Paolo A. Sossi,
Dan J. Bower
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 (Ghiorso et al., 1995) with gas-species properties from multiple thermochemistry tables (e.g., Chase et al., 1998). VapoRock calculates the partial p…
▽ More
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 (Ghiorso et al., 1995) with gas-species properties from multiple thermochemistry tables (e.g., Chase et al., 1998). 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 (Na, Ca, & Al) or more realistic (K, Si, Mg, Fe, & Ti) than those of the closed-source MAGMA code (with maximum deviations by factors of 10-300 for K & 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 pseudo-species, which only roughly approximates the non-ideal compositional interactions typical of many-component natural silicate melts. Finally, we explore how relative abundances of SiO and SiO2 provide a spectroscopically measurable proxy for oxygen fugacity in devolatilized exoplanetary atmospheres, potentially constraining fO2 in outgassed exoplanetary mantles.
△ Less
Submitted 20 December, 2022; v1 submitted 19 August, 2022;
originally announced August 2022.
-
Evolution of Mercury's Earliest Atmosphere
Authors:
Noah Jäggi,
Diana Gamborino,
Dan J. Bower,
Paolo A. Sossi,
Aaron S. Wolf,
Apurva V. Oza,
Audrey Vorburger,
André Galli,
Peter Wurz
Abstract:
MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, a…
▽ More
MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, and atmospheric composition. Consideration of two initial Mercury sizes and four magma ocean compositions determine the atmospheric speciation at a given surface temperature. A coupled interior-atmosphere model determines the cooling rate and therefore the lifetime of surficial melt. Combining the melt lifetime and escape flux calculations provide estimates for the total mass loss from early Mercury. Loss rates by Jeans escape are negligible. Plasma heating and photoionization are limited by homopause diffusion rates of $\sim10^{6}$ kg/s. Loss by photoevaporation depends on the timing of Mercury formation and assumed heating efficiency and ranges from $\sim10^{6.6}$ to $\sim10^{9.6}$ kg/s. The material for photoevaporation is sourced from below the homopause and is therefore energy-limited rather than diffusion-limited. The timescale for efficient interior-atmosphere chemical exchange is less than ten thousand years. Therefore, escape processes only account for an equivalent loss of less than 2.3 km of crust ($0.3\%$ of Mercury's mass). Accordingly, $\leq0.02\%$ of the total mass of H$_2$O and Na is lost. Therefore, cumulative loss cannot significantly modify Mercury's bulk mantle composition during the magma ocean stage. Mercury's high core:mantle ratio and volatile-rich surface may instead reflect chemical variations in its building blocks resulting from its solar-proximal accretion environment.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
Retention of water in terrestrial magma oceans and carbon-rich early atmospheres
Authors:
Dan J. Bower,
Kaustubh Hakim,
Paolo A. Sossi,
Patrick Sanan
Abstract:
Massive steam and CO$_2$ atmospheres have been proposed for magma ocean outgassing of Earth and terrestrial planets. Yet formation of such atmospheres depends on volatile exchange with the molten interior, governed by volatile solubilities and redox reactions. We determine the evolution of magma ocean--atmosphere systems for a range of oxygen fugacities, C/H ratios and hydrogen budgets that includ…
▽ More
Massive steam and CO$_2$ atmospheres have been proposed for magma ocean outgassing of Earth and terrestrial planets. Yet formation of such atmospheres depends on volatile exchange with the molten interior, governed by volatile solubilities and redox reactions. We determine the evolution of magma ocean--atmosphere systems for a range of oxygen fugacities, C/H ratios and hydrogen budgets that include redox reactions for hydrogen (H$_2$--H$_2$O), carbon (CO--CO$_2$), methane (CH$_4$), and solubility laws for H$_2$O and CO$_2$. We find that small initial budgets of hydrogen, high C/H ratios, and oxidizing conditions, suppress outgassing of hydrogen until the late stage of magma ocean crystallization. Hence early atmospheres in equilibrium with magma oceans are dominantly carbon-rich, and specifically CO-rich except at the most oxidizing conditions. The high solubility of H$_2$O limits its outgassing to melt fractions below $\sim$30\%, the fraction at which the mantle transitions from vigorous to sluggish convection with melt percolation. Sluggish melt percolation could enable a surface lid to form, trapping water in the interior and thereby maintaining a carbon-rich atmosphere (equilibrium crystallization). Alternatively, efficient crystal settling could maintain a molten surface, promoting a transition to a water-rich atmosphere (fractional crystallization). However, additional processes, including melt trapping and H dissolution in crystallizing minerals, further conspire to limit the extent of H outgassing, even for fractional crystallization. Hence, much of the water delivered to planets during their accretion can be safely harbored in their interiors during the magma ocean stage, particularly at oxidizing conditions.
△ Less
Submitted 14 April, 2022; v1 submitted 15 October, 2021;
originally announced October 2021.
-
Hemispheric Tectonics on super-Earth LHS 3844b
Authors:
Tobias G. Meier,
Dan J. Bower,
Tim Lichtenberg,
Paul J. Tackley,
Brice-Olivier Demory
Abstract:
The tectonic regime of rocky planets fundamentally influences their long-term evolution and cycling of volatiles between interior and atmosphere. Earth is the only known planet with active plate tectonics, but observations of exoplanets may deliver insights into the diversity of tectonic regimes beyond the solar system. Observations of the thermal phase curve of super-Earth LHS 3844b reveal a soli…
▽ More
The tectonic regime of rocky planets fundamentally influences their long-term evolution and cycling of volatiles between interior and atmosphere. Earth is the only known planet with active plate tectonics, but observations of exoplanets may deliver insights into the diversity of tectonic regimes beyond the solar system. Observations of the thermal phase curve of super-Earth LHS 3844b reveal a solid surface and lack of a substantial atmosphere, with a temperature contrast between the substellar and antistellar point of around 1000 K. Here, we use these constraints on the planet's surface to constrain the interior dynamics and tectonic regimes of LHS 3844b using numerical models of interior flow. We investigate the style of interior convection by assessing how upwellings and downwellings are organized and how tectonic regimes manifest. We discover three viable convective regimes with a mobile surface: (1) spatially uniform distribution of upwellings and downwellings, (2) prominent downwelling on the dayside and upwellings on the nightside, and (3) prominent downwelling on the nightside and upwellings on the dayside. Hemispheric tectonics is observed for regimes (2) and (3) as a direct consequence of the day-to-night temperature contrast. Such a tectonic mode is absent in the present-day solar system and has never been inferred from astrophysical observations of exoplanets. Our models offer distinct predictions for volcanism and outgassing linked to the tectonic regime, which may explain secondary features in phase curves and allow future observations to constrain the diversity of super-Earth interiors.
△ Less
Submitted 3 March, 2021;
originally announced March 2021.
-
Vertically resolved magma ocean-protoatmosphere evolution: H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, and N$_2$ as primary absorbers
Authors:
Tim Lichtenberg,
Dan J. Bower,
Mark Hammond,
Ryan Boukrouche,
Patrick Sanan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositi…
▽ More
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically-resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, or N$_2$. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N$_2$, and O$_2$ with minimal effect, H$_2$O, CO$_2$, and CH$_4$ with intermediate influence, and H$_2$ with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multi-wavelength astronomical observations.
△ Less
Submitted 26 January, 2021;
originally announced January 2021.
-
Lithologic Controls on Silicate Weathering Regimes of Temperate Planets
Authors:
Kaustubh Hakim,
Dan J. Bower,
Meng Tian,
Russell Deitrick,
Pierre Auclair-Desrotour,
Daniel Kitzmann,
Caroline Dorn,
Klaus Mezger,
Kevin Heng
Abstract:
Weathering of silicate rocks at a planetary surface can draw down CO$_2$ from the atmosphere for eventual burial and long-term storage in the planetary interior. This process is thought to provide an essential negative feedback to the carbonate-silicate cycle (carbon cycle) to maintain clement climates on Earth and potentially similar temperate exoplanets. We implement thermodynamics to determine…
▽ More
Weathering of silicate rocks at a planetary surface can draw down CO$_2$ from the atmosphere for eventual burial and long-term storage in the planetary interior. This process is thought to provide an essential negative feedback to the carbonate-silicate cycle (carbon cycle) to maintain clement climates on Earth and potentially similar temperate exoplanets. We implement thermodynamics to determine weathering rates as a function of surface lithology (rock type). These rates provide upper limits that allow estimating the maximum rate of weathering in regulating climate. This modeling shows that the weathering of mineral assemblages in a given rock, rather than individual minerals, is crucial to determine weathering rates at planetary surfaces. By implementing a fluid-transport controlled approach, we further mimic chemical kinetics and thermodynamics to determine weathering rates for three types of rocks inspired by the lithologies of Earth's continental and oceanic crust, and its upper mantle. We find that thermodynamic weathering rates of a continental crust-like lithology are about one to two orders of magnitude lower than those of a lithology characteristic of the oceanic crust. We show that when the CO$_2$ partial pressure decreases or surface temperature increases, thermodynamics rather than kinetics exerts a strong control on weathering. The kinetically- and thermodynamically-limited regimes of weathering depend on lithology, whereas, the supply-limited weathering is independent of lithology. Our results imply that the temperature-sensitivity of thermodynamically-limited silicate weathering may instigate a positive feedback to the carbon cycle, in which the weathering rate decreases as the surface temperature increases.
△ Less
Submitted 1 February, 2021; v1 submitted 26 August, 2020;
originally announced August 2020.
-
The influence of bulk composition on long-term interior-atmosphere evolution of terrestrial exoplanets
Authors:
Rob J. Spaargaren,
Maxim D. Ballmer,
Dan J. Bower,
Caroline Dorn,
Paul J. Tackley
Abstract:
Aims: The secondary atmospheres of terrestrial planets form and evolve as a consequence of interaction with the interior over geological time. We aim to quantify the influence of planetary bulk composition on the interior--atmosphere evolution for Earth-sized terrestrial planets to aid in the interpretation of future observations of terrestrial exoplanet atmospheres. Methods: We used a geochemical…
▽ More
Aims: The secondary atmospheres of terrestrial planets form and evolve as a consequence of interaction with the interior over geological time. We aim to quantify the influence of planetary bulk composition on the interior--atmosphere evolution for Earth-sized terrestrial planets to aid in the interpretation of future observations of terrestrial exoplanet atmospheres. Methods: We used a geochemical model to determine the major-element composition of planetary interiors (MgO, FeO, and SiO2) following the crystallization of a magma ocean after planet formation, predicting a compositional profile of the interior as an initial condition for our long-term thermal evolution model. Our 1D evolution model predicts the pressure-temperature structure of the interior, which we used to evaluate near-surface melt production and subsequent volatile outgassing. Volatiles are exchanged between the interior and atmosphere according to mass conservation. Results: Based on stellar compositions reported in the Hypatia catalog, we predict that about half of rocky exoplanets have a mantle that convects as a single layer (whole-mantle convection), and the other half exhibit double-layered convection due to the presence of a mid-mantle compositional boundary. Double-layered convection is more likely for planets with high bulk planetary Fe-content and low Mg/Si-ratio. We find that planets with low Mg/Si-ratio tend to cool slowly because their mantle viscosity is high. Accordingly, low-Mg/Si planets also tend to lose volatiles swiftly through extensive melting. Moreover, the dynamic regime of the lithosphere (plate tectonics vs.\ stagnant lid) has a first-order influence on the thermal evolution and volatile cycling. These results suggest that the composition of terrestrial exoplanetary atmospheres can provide information on the dynamic regime of the lithosphere and the thermo-chemical evolution of the interior.
△ Less
Submitted 5 October, 2020; v1 submitted 17 July, 2020;
originally announced July 2020.
-
Linking the evolution of terrestrial interiors and an early outgassed atmosphere to astrophysical observations
Authors:
Dan J. Bower,
Daniel Kitzmann,
Aaron S. Wolf,
Patrick Sanan,
Caroline Dorn,
Apurva V. Oza
Abstract:
A terrestrial planet is molten during formation and may remain so if subject to intense insolation or tidal forces. Observations continue to favour the detection and characterisation of hot planets, potentially with large outgassed atmospheres. We aim to determine the radius of hot Earth-like planets with large outgassed atmospheres and explore differences between molten and solid silicate planets…
▽ More
A terrestrial planet is molten during formation and may remain so if subject to intense insolation or tidal forces. Observations continue to favour the detection and characterisation of hot planets, potentially with large outgassed atmospheres. We aim to determine the radius of hot Earth-like planets with large outgassed atmospheres and explore differences between molten and solid silicate planets and their influence on the mass-radius relationship and transmission and emission spectra. An interior-atmosphere model, combined with static structure calculations, tracks the evolving radius of a rocky mantle that is outgassing CO$_2$ and H$_2$O. Synthetic emission and transmission spectra are generated for CO$_2$ and H$_2$O dominated atmospheres. Atmospheres dominated by CO$_2$ suppress the outgassing of H$_2$O to a greater extent than previously realised, as previous studies have applied an erroneous relationship between volatile mass and partial pressure. We therefore predict more H$_2$O can be retained by the interior during the later stages of magma ocean crystallisation. Furthermore, formation of a lid at the surface can tie outgassing of H$_2$O to the efficiency of heat transport through the lid, rather than the atmosphere's radiative timescale. Contraction of the mantle as it solidifies gives $\sim5\%$ radius decrease, which can partly be offset by addition of a relatively light species to the atmosphere. We conclude that a molten silicate mantle can increase the radius of a terrestrial planet by around $5\%$ compared to its solid counterpart, or equivalently account for a $13\%$ decrease in bulk density. An outgassing atmosphere can perturb the total radius according to its speciation. Atmospheres of terrestrial planets around M-stars that are dominated by CO$_2$ or H$_2$O can be distinguished by observing facilities with extended wavelength coverage (e.g., JWST).
△ Less
Submitted 7 November, 2019; v1 submitted 17 April, 2019;
originally announced April 2019.
-
Direct imaging of molten protoplanets in nearby young stellar associations
Authors:
Irene Bonati,
Tim Lichtenberg,
Dan J Bower,
Miles L Timpe,
Sascha P Quanz
Abstract:
During their formation and early evolution, rocky planets undergo multiple global melting events due to accretionary collisions with other protoplanets. The detection and characterization of their post-collision afterglows (magma oceans) can yield important clues about the origin and evolution of the solar and extrasolar planet population. Here, we quantitatively assess the observational prospects…
▽ More
During their formation and early evolution, rocky planets undergo multiple global melting events due to accretionary collisions with other protoplanets. The detection and characterization of their post-collision afterglows (magma oceans) can yield important clues about the origin and evolution of the solar and extrasolar planet population. Here, we quantitatively assess the observational prospects to detect the radiative signature of forming planets covered by such collision-induced magma oceans in nearby young stellar associations with future direct imaging facilities. We have compared performance estimates for near- and mid-infrared instruments to be installed at ESO's Extremely Large Telescope (ELT), and a potential space-based mission called Large Interferometer for Exoplanets (LIFE). We modelled the frequency and timing of energetic collisions using \textit{N}-body models of planet formation for different stellar types, and determine the cooling of the resulting magma oceans with an insulating atmosphere. We find that the probability of detecting at least one magma ocean planet depends on the observing duration and the distribution of atmospheric properties among rocky protoplanets. However, the prospects for detection significantly increase for young and close stellar targets, which show the highest frequencies of giant impacts. For intensive reconnaissance with a K band (2.2 $μm$) ELT filter or a 5.6 $μm$ LIFE filter, the $β$ Pictoris, Columba, TW Hydrae, and Tucana-Horologium associations represent promising candidates for detecting a molten protoplanet. Our results motivate the exploration of magma ocean planets using the ELT and underline the importance of space-based direct imaging facilities to investigate and characterize planet formation and evolution in the solar vicinity.
△ Less
Submitted 18 November, 2018;
originally announced November 2018.
-
Numerical solution of a non-linear conservation law applicable to the interior dynamics of partially molten planets
Authors:
Dan J. Bower,
Patrick Sanan,
Aaron S. Wolf
Abstract:
The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. In this formulation the effective or eddy diffusivity depends on the entropy gradient, $\partial S/\partial r$, as well as entropy. First we present a simplified model with semi-an…
▽ More
The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. In this formulation the effective or eddy diffusivity depends on the entropy gradient, $\partial S/\partial r$, as well as entropy. First we present a simplified model with semi-analytical solutions, highlighting the large dynamic range of $\partial S/\partial r$, around 12 orders of magnitude, for physically-relevant parameters. It also elucidates the thermal structure of a magma ocean during the earliest stage of crystal formation. This motivates the development of a simple, stable numerical scheme able to capture the large dynamic range of $\partial S/\partial r$ and provide a flexible and robust method for time-integrating the energy equation.
We then consider a full model including energy fluxes associated with convection, mixing, gravitational separation, and conduction that all depend on the thermophysical properties of the melt and solid phases. This model is discretised and evolved by applying the finite volume method (FVM), allowing for extended precision calculations and using $\partial S/\partial r$ as the solution variable. The FVM is well-suited to this problem since it is naturally energy conserving, flexible, and intuitive to incorporate arbitrary non-linear fluxes that rely on lookup data. Special attention is given to the numerically challenging scenario in which crystals first form in the centre of a magma ocean.
Our computational framework is immediately applicable to modelling high melt fraction phenomena in Earth and planetary science research. Furthermore, it provides a template for solving similar non-linear diffusion equations arising in other disciplines, particularly for non-linear functional forms of the diffusion coefficient.
△ Less
Submitted 20 November, 2017;
originally announced November 2017.
-
Assessing the Interior Structure of Terrestrial Exoplanets with Implications for Habitability
Authors:
Caroline Dorn,
Dan J. Bower,
Antoine B. Rozel
Abstract:
Astrophysical observations reveal a large diversity of radii and masses of exoplanets. It is important to characterize the interiors of exoplanets to understand planetary diversity and further determine how unique, or not, Earth is. Assessing interior structure is challenging because there are few data and large uncertainties. Thus, for a given exoplanet a range of interior structure models can sa…
▽ More
Astrophysical observations reveal a large diversity of radii and masses of exoplanets. It is important to characterize the interiors of exoplanets to understand planetary diversity and further determine how unique, or not, Earth is. Assessing interior structure is challenging because there are few data and large uncertainties. Thus, for a given exoplanet a range of interior structure models can satisfy available data. Typically, interior models aim to constrain the radial structure and composition of the core and mantle, and additionally ice, ocean, and gas layer if appropriate. Constraining the parameters of these layers may also inform us about interior dynamics. However, it remains challenging to constrain interior dynamics using interior structure models because structure models are relatively insensitive to the thermal state of a planet. Nevertheless, elucidating interior dynamics remains a key goal in exoplanetology due to its role in determining surface conditions and hence habitability. Thus far, Earth-like habitability can be excluded for super-Earths that are in close proximity to their stars and therefore have high surface temperatures that promote local magma oceans.
△ Less
Submitted 16 October, 2017;
originally announced October 2017.