Showing 1–40 of 40 results for author: Pierrehumbert, R

Carbon cycle instability for high-$\mathrm{CO_2}$ exoplanets: Implications for habitability

Authors: R. J. Graham, R. T. Pierrehumbert

Abstract: Implicit in the definition of the classical circumstellar habitable zone (HZ) is the hypothesis that the carbonate-silicate cycle can maintain clement climates on exoplanets with land and surface water across a range of instellations by adjusting atmospheric $\mathrm{CO_2}$ partial pressure ($p\mathrm{CO_2}$). This hypothesis is made by analogy to the Earth system, but it is an open question wheth… ▽ More Implicit in the definition of the classical circumstellar habitable zone (HZ) is the hypothesis that the carbonate-silicate cycle can maintain clement climates on exoplanets with land and surface water across a range of instellations by adjusting atmospheric $\mathrm{CO_2}$ partial pressure ($p\mathrm{CO_2}$). This hypothesis is made by analogy to the Earth system, but it is an open question whether silicate weathering can stabilize climate on planets in the outer reaches of the HZ, where instellations are lower than those received by even the Archean Earth and $\mathrm{CO_2}$ is thought likely to dominate atmospheres. Since weathering products are carried from land to ocean by the action of water, silicate weathering is intimately coupled to the hydrologic cycle, which intensifies with hotter temperatures under Earth-like conditions. Here, we use global climate model (GCM) simulations to demonstrate that the hydrologic cycle responds counterintuitively to changes in climate on planets with $\mathrm{CO_2}$-$\mathrm{H_2O}$ atmospheres at low instellations and high $p\mathrm{CO_2}$, with global evaporation and precipitation decreasing as $p\mathrm{CO_2}$ and temperatures increase at a given instellation. Within the MAC weathering formulation, weathering then decreases with increasing $p\mathrm{CO_2}$ for a range of instellations and $p\mathrm{CO_2}$ typical of the outer reaches of the HZ, resulting in an unstable carbon cycle that may lead to either runaway $\mathrm{CO_2}$ accumulation or depletion of $\mathrm{CO_2}$ to colder (possibly Snowball) conditions. While the behavior of the system has not been completely mapped out, the results suggest that silicate weathering could fail to maintain habitable conditions in the outer reaches of the nominal HZ. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ

JWST Reveals CH$_4$, CO$_2$, and H$_2$O in a Metal-rich Miscible Atmosphere on a Two-Earth-Radius Exoplanet

Authors: Björn Benneke, Pierre-Alexis Roy, Louis-Philippe Coulombe, Michael Radica, Caroline Piaulet, Eva-Maria Ahrer, Raymond Pierrehumbert, Joshua Krissansen-Totton, Hilke E. Schlichting, Renyu Hu, Jeehyun Yang, Duncan Christie, Daniel Thorngren, Edward D. Young, Stefan Pelletier, Heather A. Knutson, Yamila Miguel, Thomas M. Evans-Soma, Caroline Dorn, Anna Gagnebin, Jonathan J. Fortney, Thaddeus Komacek, Ryan MacDonald, Eshan Raul, Ryan Cloutier , et al. (6 additional authors not shown)

Abstract: Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2… ▽ More Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2$\,μ\mathrm{m}$ JWST NIRISS/SOSS+NIRSpec/G395H transmission spectrum of the 2.2$\,R_\oplus$ TOI-270d ($4.78\,M_\oplus$, $T_\mathrm{eq}$=350-380 K), delivering unprecedented sensitivity for atmospheric characterization in the sub-Neptune regime. We detect five vibrational bands of CH$_4$ at 1.15, 1.4, 1.7, 2.3, and 3.3$\,μ$m (9.4$σ$), the signature of CO$_2$ at 4.3$\,μ$m (4.8$σ$), water vapor (2.5$σ$), and potential signatures of SO$_2$ at 4.0$\,μ\mathrm{m}$ and CS$_2$ at 4.6$\,μ\mathrm{m}$. Intriguingly, we find an overall highly metal-rich atmosphere, with a mean molecular weight of $5.47_{-1.14}^{+1.25}$. We infer an atmospheric metal mass fraction of $58_{-12}^{+8}\%$ and a C/O of $0.47_{-0.19}^{+0.16}$, indicating that approximately half the mass of the outer envelope is in high-molecular-weight volatiles (H$_2$O, CH$_4$, CO, CO$_2$) rather than H$_2$/He. We introduce a sub-Neptune classification scheme and identify TOI-270d as a "miscible-envelope sub-Neptune" in which H$_2$/He is well-mixed with the high-molecular-weight volatiles in a miscible supercritical metal-rich envelope. For a fully miscible envelope, we conclude that TOI-270d's interior is $90_{-4}^{+3}\,$wt$\,\%$ rock/iron, indicating that it formed as a rocky planet that accreted a few wt % of H$_2$/He, with the overall envelope metal content explained by magma-ocean/envelope reactions without the need for significant ice accretion. TOI-270d may well be an archetype of the overall population of sub-Neptunes. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 25 pages, 12 figures

Modeling the day-night temperature variations of ultra-hot Jupiters: confronting non-grey general circulation models and observations

Authors: Xianyu Tan, Thaddeus D. Komacek, Natasha E. Batalha, Drake Deming, Roxana Lupu, Vivien Parmentier, Raymond T. Pierrehumbert

Abstract: Ultra-hot Jupiters (UHJs) are natural laboratories to study extreme physics in planetary atmospheres and their rich observational data sets are yet to be confronted with models with varying complexities at a population level. In this work, we update the general circulation model of Tan & Komacek (2019) to include a non-grey radiative transfer scheme and apply it to simulate the realistic thermal s… ▽ More Ultra-hot Jupiters (UHJs) are natural laboratories to study extreme physics in planetary atmospheres and their rich observational data sets are yet to be confronted with models with varying complexities at a population level. In this work, we update the general circulation model of Tan & Komacek (2019) to include a non-grey radiative transfer scheme and apply it to simulate the realistic thermal structures, phase-dependent spectra, and wavelength-dependent phase curves of UHJs. We performed grids of models over a large range of equilibrium temperatures and rotation periods for varying assumptions, showing that the fractional day-night brightness temperature differences remain almost constant or slightly increase with increasing equilibrium temperature from the visible to mid-infrared wavelengths. This differs from previous work primarily due to the increasing planetary rotation rate with increasing equilibrium temperature for fixed host star type. Radiative effects of varying atmospheric compositions become more significant in dayside brightness temperature in longer wavelengths. Data-model comparisons of dayside brightness temperatures and phase curve amplitudes as a function of equilibrium temperature are in broad agreement. Observations show a large scatter compared to models even with a range of different assumptions, indicating significantly varying intrinsic properties in the hot Jupiter population. Our cloud-free models generally struggle to match all observations for individual targets with a single set of parameter choices, indicating the need for extra processes for understanding the heat transport of UHJs. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS; data underlying this article is available in Zenodo https://doi.org/10.5281/zenodo.10121933

Global Chemical Transport on Hot Jupiters: Insights from 2D VULCAN photochemical model

Authors: Shang-Min Tsai, Vivien Parmentier, João M. Mendonça, Xianyu Tan, Russell Deitrick, Mark Hammond, Arjun B. Savel, Xi Zhang, Raymond T. Pierrehumbert, Edward W. Schwieterman

Abstract: The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photoc… ▽ More The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below $\sim$0.1 mbar while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C$_2$H$_4$ absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres. △ Less

Submitted 25 March, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 19 pages, 20 figures, published in ApJ

Journal ref: Tsai et al. (2024), ApJ, 963, 41

Modelling Non-Condensing Compositional Convection for Applications to Super-Earth and Sub-Neptune Atmospheres

Authors: Namrah Habib, Raymond T. Pierrehumbert

Abstract: Compositional convection is atmospheric mixing driven by density variations caused by compositional gradients. Previous studies have suggested that compositional gradients of atmospheric trace species within planetary atmospheres can impact convection and the final atmospheric temperature profile. In this work, we employ 3D convection resolving simulations using Cloud Model 1 (CM1) to gain a funda… ▽ More Compositional convection is atmospheric mixing driven by density variations caused by compositional gradients. Previous studies have suggested that compositional gradients of atmospheric trace species within planetary atmospheres can impact convection and the final atmospheric temperature profile. In this work, we employ 3D convection resolving simulations using Cloud Model 1 (CM1) to gain a fundamental understanding of how compositional variation influences convection and the final atmospheric state of exoplanet atmospheres. We perform 3D initial value problem simulations of non-condensing compositional convection for Earth-Air, $\rm H_2$, and $\rm CO_2$ atmospheres. Conventionally, atmospheric convection is assumed to mix the atmosphere to a final, marginally stable state defined by a unique temperature profile. However, when there is compositional variation within an atmosphere, a continuous family of stable end states is possible, differing in the final state composition profile. Our CM1 simulations are used to determine which of the family of possible compositional end states is selected. Leveraging the results from our CM1 simulations, we develop a dry convective adjustment scheme for use in General Circulation Models (GCMs). This scheme relies on an energy analysis to determine the final adjusted atmospheric state. Our convection scheme produces results that agree with our CM1 simulations and can easily be implemented in GCMs to improve modelling of compositional convection in exoplanet atmospheres. △ Less

Submitted 12 October, 2023; originally announced October 2023.

The Runaway Greenhouse Effect on Hycean Worlds

Authors: Hamish Innes, Shang-Min Tsai, Raymond T. Pierrehumbert

Abstract: Hycean worlds are a proposed subset of sub-Neptune exoplanets with substantial water inventories, liquid surface oceans and extended hydrogen-dominated atmospheres that could be favourable for habitability. In this work, we aim to quantitatively define the inner edge of the Hycean habitable zone using a 1D radiative-convective model. As a limiting case, we model a dry hydrogen-helium envelope abov… ▽ More Hycean worlds are a proposed subset of sub-Neptune exoplanets with substantial water inventories, liquid surface oceans and extended hydrogen-dominated atmospheres that could be favourable for habitability. In this work, we aim to quantitatively define the inner edge of the Hycean habitable zone using a 1D radiative-convective model. As a limiting case, we model a dry hydrogen-helium envelope above a surface ocean. We find that 10 to 20 bars of atmosphere produces enough greenhouse effect to drive a liquid surface ocean supercritical when forced with current Earth-like instellation. Introducing water vapour into the atmosphere, we show the runaway greenhouse instellation limit is greatly reduced due to the presence of superadiabatic layers where convection is inhibited. This moves the inner edge of the habitable zone from $\approx$ 1 AU for a G-star to 1.6 AU (3.85 AU) for a Hycean world with a H$_2$-He inventory of 1 bar (10 bar). For an M-star, the inner edge is equivalently moved from 0.17 AU to 0.28 AU (0.54 AU). Our results suggest that most of the current Hycean world observational targets are not likely to sustain a liquid water ocean. We present an analytical framework for interpreting our results, finding that the maximum possible OLR scales approximately inversely with the dry mass inventory of the atmosphere. We discuss the possible limitations of our 1D modelling and recommend the use of 3D convection-resolving models to explore the robustness of superadiabatic layers. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ, comments welcome

The runaway greenhouse on subNeptune waterworlds

Authors: Raymond T. Pierrehumbert

Abstract: The implications of the water vapor runaway greenhouse phenomenon for water-rich subNeptunes are developed. In particular, the nature of the post-runaway equilibration process for planets that have an extremely high water inventory is addressed. Crossing the threshold from sub-runaway to super-runaway conditions leads to a transition from equilibrated states with cold deep liquid oceans and deep i… ▽ More The implications of the water vapor runaway greenhouse phenomenon for water-rich subNeptunes are developed. In particular, the nature of the post-runaway equilibration process for planets that have an extremely high water inventory is addressed. Crossing the threshold from sub-runaway to super-runaway conditions leads to a transition from equilibrated states with cold deep liquid oceans and deep interior ice-X phases to states with hot supercritical fluid interiors. There is a corresponding marked inflation of radius for a given mass, similar to the runaway greenhouse radius inflation effect noted earlier for terrestrial planets, but in the present case the inflation involves the entire interior of the planet. The calculation employs the AQUA equation of state database to simplify the internal structure calculation. Some speculations concerning the effect of $\mathrm{H_2}$ admixture, silicate cores and hot vs. cold start evolution trajectories are offered. Observational implications are discussed, though the search for the mass-radius signature of the phenomena considered is limited by degeneracies and by lack of data. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: Submitted to Astrophysical Journal

CO$_2$ ocean bistability on terrestrial exoplanets

Authors: R. J. Graham, Tim Lichtenberg, Ray Pierrehumbert

Abstract: Cycling of carbon dioxide between the atmosphere and interior of rocky planets can stabilize global climate and enable planetary surface temperatures above freezing over geologic time. However, variations in global carbon budget and unstable feedback cycles between planetary sub-systems may destabilize the climate of rocky exoplanets toward regimes unknown in the Solar System. Here, we perform cle… ▽ More Cycling of carbon dioxide between the atmosphere and interior of rocky planets can stabilize global climate and enable planetary surface temperatures above freezing over geologic time. However, variations in global carbon budget and unstable feedback cycles between planetary sub-systems may destabilize the climate of rocky exoplanets toward regimes unknown in the Solar System. Here, we perform clear-sky atmospheric radiative transfer and surface weathering simulations to probe the stability of climate equilibria for rocky, ocean-bearing exoplanets at instellations relevant for planetary systems in the outer regions of the circumstellar habitable zone. Our simulations suggest that planets orbiting G- and F-type stars (but not M-type stars) may display bistability between an Earth-like climate state with efficient carbon sequestration and an alternative stable climate equilibrium where CO$_2$ condenses at the surface and forms a blanket of either clathrate hydrate or liquid CO$_2$. At increasing instellation and with ineffective weathering, the latter state oscillates between cool, surface CO$_2$-condensing and hot, non-condensing climates. CO$_2$ bistable climates may emerge early in planetary history and remain stable for billions of years. The carbon dioxide-condensing climates follow an opposite trend in $p$CO$_2$ versus instellation compared to the weathering-stabilized planet population, suggesting the possibility of observational discrimination between these distinct climate categories. △ Less

Submitted 30 September, 2022; originally announced October 2022.

Comments: Accepted for publication in JGR-Planets. 19 pages (including references and Tables), 6 figures, corresponding author email: robert.graham@physics.ox.ac.uk

Journal ref: Journal of Geophysical Research: Planets, 127, e2022JE007456

The Climate and Compositional Variation of the Highly Eccentric Planet HD 80606 b -- the rise and fall of carbon monoxide and elemental sulfur

Authors: Shang-Min Tsai, Maria Steinrueck, Vivien Parmentier, Nikole Lewis, Raymond Pierrehumbert

Abstract: The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chem… ▽ More The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chemical mechanisms of HD 80606 b. We first run three-dimensional general circulation models (GCMs) to establish the atmospheric thermal and dynamical structures for different atmospheric metallicities and internal heat. Based on the GCM output, we then adopted a 1D time-dependent photochemical model to investigate the compositional variation along the eccentric orbit. The transition of the circulation patterns of HD 80606 b matched the dynamics regimes in previous works. Our photochemical models show that efficient vertical mixing leads to deep quench levels of the major carbon and nitrogen species and the quenching behavior does not change throughout the eccentric orbit. Instead, photolysis is the main driver of the time-dependent chemistry. While CH$_4$ dominates over CO through most of the orbits, a transient state of [CO]/[CH$_4$}] $>$ 1 after periastron is confirmed for all metallicity and internal heat cases. The upcoming JWST Cycle 1 GO program will be able to track this real-time CH$_4$--CO conversion and infer the chemical timescale. Furthermore, sulfur species initiated by sudden heating and photochemical forcing exhibit both short-term and long-term cycles, opening an interesting avenue for detecting sulfur on exoplanets. △ Less

Submitted 16 March, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 18 pages, 18 figures, published in MNRAS

The impact of ultraviolet heating and cooling on the dynamics and observability of lava planet atmospheres

Authors: T. Giang Nguyen, Nicolas B. Cowan, Raymond T. Pierrehumbert, Roxana E. Lupu, John E. Moores

Abstract: Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we coup… ▽ More Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we couple the two processes by introducing ultraviolet and infrared radiation to a turbulent boundary layer model. We also test the effect of different vertical temperature profiles on atmospheric dynamics. Results from the model show that UV radiation affects the atmosphere much more than IR. UV heating and cooling work together to produce a horizontally isothermal atmosphere away from the sub-stellar point regardless of the vertical temperature profile. We also find that stronger temperature inversions induce stronger winds and hence cool the atmosphere. Our simulated transmission spectra of the bound atmosphere show a strong SiO feature in the UV that would be challenging to observe in the planet's transit spectrum due to the precision required. Our simulated emission spectra are more promising, with significant SiO spectral features at 4.5 and 9 $μ$m that can be observed with the James Webb Space Telescope. Different vertical temperature profiles produce discernible dayside emission spectra, but not in the way one would expect. △ Less

Submitted 9 May, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: Accepted to MNRAS

A Mini-Chemical Scheme with Net Reactions for 3D GCMs I.: Thermochemical Kinetics

Authors: Shang-Min Tsai, Elspeth K. H. Lee, Raymond Pierrehumbert

Abstract: Growing evidence has indicated that the global composition distribution plays an indisputable role in interpreting observational data. 3D general circulation models (GCMs) with a reliable treatment of chemistry and clouds are particularly crucial in preparing for the upcoming observations. In the effort of achieving 3D chemistry-climate modeling, the challenge mainly lies in the expensive computin… ▽ More Growing evidence has indicated that the global composition distribution plays an indisputable role in interpreting observational data. 3D general circulation models (GCMs) with a reliable treatment of chemistry and clouds are particularly crucial in preparing for the upcoming observations. In the effort of achieving 3D chemistry-climate modeling, the challenge mainly lies in the expensive computing power required for treating a large number of chemical species and reactions. Motivated by the need for a robust and computationally efficient chemical scheme, we devise a mini-chemical network with a minimal number of species and reactions for H$_2$-dominated atmospheres. We apply a novel technique to simplify the chemical network from a full kinetics model -- VULCAN by replacing a large number of intermediate reactions with net reactions. The number of chemical species is cut down from 67 to 12, with the major species of thermal and observational importance retained, including H$_2$O, CH$_4$, CO, CO$_2$, C$_2$H$_2$, NH$_3$, and HCN. The size of the total reactions is greatly reduced from $\sim$ 800 to 20. The mini-chemical scheme is validated by verifying the temporal evolution and benchmarking the predicted compositions in four exoplanet atmospheres (GJ 1214b, GJ 436b, HD 189733b, HD 209458b) against the full kinetics of VULCAN. It reproduces the chemical timescales and composition distributions of the full kinetics well within an order of magnitude for the major species in the pressure range of 1 bar -- 0.1 mbar across various metallicities and carbon-to-oxygen (C/O) ratios. The small scale of the mini-chemical scheme permits simple use and fast computation, which is optimal for implementation in a 3D GCM or a retrieval framework. We focus on the thermochemical kinetics of net reactions in this paper and address photochemistry in a follow-up paper. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 9 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 664, A82 (2022)

Cloud-convection feedback in brown dwarfs atmosphere

Authors: Maxence Lefèvre, Xianyu Tan, Elspeth K. H. Lee, R. T. Pierrehumbert

Abstract: Numerous observational evidence has suggested the presence of active meteorology in the atmospheres of brown dwarfs. A near-infrared brightness variability has been observed. Clouds have a major role in shaping the thermal structure and spectral properties of these atmospheres. The mechanism of such variability is still unclear and both 1D and global circulation model cannot fully study this topic… ▽ More Numerous observational evidence has suggested the presence of active meteorology in the atmospheres of brown dwarfs. A near-infrared brightness variability has been observed. Clouds have a major role in shaping the thermal structure and spectral properties of these atmospheres. The mechanism of such variability is still unclear and both 1D and global circulation model cannot fully study this topics due to resolution. In this study, a convective resolving model is coupled to grey-band radiative transfer in order to study the coupling between the convective atmosphere and the variability of clouds over a large temperature range with a domain of several hundreds of kilometers. Six types of clouds are considered, with microphysics including settling. The clouds are radiatively active using Rosseland mean coefficient. Radiative cloud feedback can drive spontaneous atmospheric variability in both temperature and cloud structure, as modeled for the first time in three dimensions. Silicate clouds have the most effect of the thermal structure with the generation of a secondary convective layer in some cases, depending on the assumed particle size. Iron and Aluminum clouds also have a substantial impact on the atmosphere. Thermal spectra were computed, and we find the strongest effect of clouds is the smoothing of spectral features at optical wavelengths. Compared to observed L and T dwarfs on color-magnitude diagram, the simulated atmospheres are redder for most of the cases. The simulations with the presence of cloud holes are closer to the observations. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 33 pages, 23 figures and 1 table. Accepted in The Astrophysical Journal

K2 and Spitzer phase curves of the rocky ultra-short-period planet K2-141 b hint at a tenuous rock vapor atmosphere

Authors: S. Zieba, M. Zilinskas, L. Kreidberg, T. G. Nguyen, Y. Miguel, N. B. Cowan, R. Pierrehumbert, L. Carone, L. Dang, M. Hammond, T. Louden, R. Lupu, L. Malavolta, K. B. Stevenson

Abstract: K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of t… ▽ More K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of the planet. We measure an infrared eclipse depth of 143 +/- 39 ppm and a peak to trough amplitude variation of 121 +/- 43 ppm. The best fit model to the Spitzer data shows no significant thermal hotspot offset, in contrast to the previously observed offset for the well-studied USP planet 55 Cnc e. We also jointly analyze the new Spitzer observations with the photometry collected by Kepler during two separate K2 campaigns. We model the planetary emission with a range of toy models that include a reflective and a thermal contribution. With a two-temperature model, we measure a dayside temperature of 2049 +/- 361 K and a night-side temperature that is consistent with zero (Tp,n < 1712 K at 2 sigma). Models with a steep dayside temperature gradient provide a better fit to the data than a uniform dayside temperature (DeltaBIC = 22.2). We also find evidence for a non-zero geometric albedo Ag = 0.28 +/- 0.07. We also compare the data to a physically motivated, pseudo-2D rock vapor model and a 1D turbulent boundary layer model. Both models fit the data well. Notably, we find that the optical eclipse depth can be explained by thermal emission from a hot inversion layer, rather than reflected light. A thermal inversion may also be responsible for the deep optical eclipse observed for another USP, Kepler-10 b. Finally, we significantly improve the ephemerides for K2-141 b and c, which will facilitate further follow-up observations of this interesting system with state-of-the-art observatories like JWST. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 31 pages, 27 figures, 12 tables, accepted for publication in A&A

Journal ref: A&A 664, A79 (2022)

Atmospheric dynamics of temperate sub-Neptunes. Part I: dry dynamics

Authors: Hamish Innes, Raymond T. Pierrehumbert

Abstract: Sub-Neptunes (planets with radii between 2 and 4 R$_{\oplus}$) are abundant around M-dwarf stars, yet the atmospheric dynamics of these planets is relatively unexplored. In this paper, we aim to provide a basic underpinning of the dry dynamics of general low mean molecular weight, temperate sub-Neptune atmospheres. We use the ExoFMS GCM with an idealised grey gas radiation scheme to simulate plane… ▽ More Sub-Neptunes (planets with radii between 2 and 4 R$_{\oplus}$) are abundant around M-dwarf stars, yet the atmospheric dynamics of these planets is relatively unexplored. In this paper, we aim to provide a basic underpinning of the dry dynamics of general low mean molecular weight, temperate sub-Neptune atmospheres. We use the ExoFMS GCM with an idealised grey gas radiation scheme to simulate planetary atmospheres with different levels of instellation and rotation rates, using the atmosphere of K2-18b as our control. We find that the atmospheres of tidally-locked, temperate sub-Neptunes have weak horizontal temperature gradients owing to their slow rotation rates and hydrogen-dominated composition. The zonal wind structure is dominated by high-latitude cyclostrophic jets driven by the conservation of angular momentum. At low pressures we observe superrotating equatorial jets, which we propose are driven by a Rossby-Kelvin instability similar to the type seen in simulations of idealised atmospheres with axisymmetric forcing. By viewing the flow in tidally-locked coordinates, we find the predominant overturning circulation to be between the dayside and nightside, and we derive scaling relations linking the tidally-locked streamfunction and vertical velocities to instellation. Comparing our results to the only other GCM study of K2-18b, we find significant qualitative differences in dynamics, highlighting the need for further collaboration and investigation into effects of different dynamical cores and physical parameterizations. This paper provides a baseline for studying the dry dynamics of temperate sub-Neptunes, which will be built on in part II with the introduction of moist effects. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJ

Convection modeling of pure-steam atmospheres

Authors: Xianyu Tan, Maxence Lefevre, Raymond Pierrehumbert

Abstract: Condensable species are crucial in shaping planetary climate. A wide range of planetary climate systems involve understanding non-dilute condensable substances and their influence on climate dynamics. There has been progress on large-scale dynamical effects and on 1D convection parameterization, but resolved 3D moist convection remains unexplored in non-dilute conditions, though it can have a prof… ▽ More Condensable species are crucial in shaping planetary climate. A wide range of planetary climate systems involve understanding non-dilute condensable substances and their influence on climate dynamics. There has been progress on large-scale dynamical effects and on 1D convection parameterization, but resolved 3D moist convection remains unexplored in non-dilute conditions, though it can have a profound impact on temperature/humidity profiles and cloud structure. We tackle this problem for pure-steam atmospheres using three-dimensional, high-resolution numerical simulations of convection in post-runaway atmospheres where the water reservoir at the surface has been exhausted. We show that the atmosphere is comprised of two characteristic regions, an upper condensing region dominated by gravity waves and a lower noncondensing region characterized by convective overturning cells. Velocities in the condensing region are much smaller than those in the lower noncondensing region, and the horizontal temperature variation is small overall. Condensation in the thermal photosphere is largely driven by radiative cooling and tends to be statistically homogeneous. Some condensation also happens deeper, near the boundary of the condensing region, due to triggering by gravity waves and convective penetrations and exhibit random patchiness. This qualitative structure is insensitive to varying model parameters, but quantitative details may differ. Our results confirm theoretical expectations that atmospheres close to the pure-steam limit do not have organized deep convective plumes in the condensing region. The generalized convective parameterization scheme discussed in Ding & Pierrehumbert (2016) is appropriate to handle the basic structure of atmospheres near the pure-steam limit but is difficult to capture gravity waves and their mixing that appear in 3D convection-resolving models. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 10 pages, 4 figures, accepted to ApJL

Inferring Shallow Surfaces on sub-Neptune Exoplanets with JWST

Authors: Shang-Min Tsai, Hamish Innes, Tim Lichtenberg, Jake Taylor, Matej Malik, Katy Chubb, Raymond Pierrehumbert

Abstract: Planets smaller than Neptune and larger than Earth make up the majority of the discovered exoplanets. Those with H$_2$-rich atmospheres are prime targets for atmospheric characterization. The transition between the two main classes, super-Earths and sub-Neptunes, is not clearly understood as the rocky surface is likely not accessible to observations. Tracking several trace gases (specifically the… ▽ More Planets smaller than Neptune and larger than Earth make up the majority of the discovered exoplanets. Those with H$_2$-rich atmospheres are prime targets for atmospheric characterization. The transition between the two main classes, super-Earths and sub-Neptunes, is not clearly understood as the rocky surface is likely not accessible to observations. Tracking several trace gases (specifically the loss of ammonia (NH$_3$) and hydrogen cyanide (HCN)) has been proposed as a proxy for the presence of a shallow surface. In this work, we revisit the proposed mechanism of nitrogen conversion in detail and find its timescale on the order of a million years. NH$_3$ exhibits dual paths converting to N$_2$ or HCN, depending on the UV radiation of the star and the stage of the system. In addition, methanol (CH$_3$OH) is identified as a robust and complementary proxy for a shallow surface. We follow the fiducial example of K2-18b with a 2D photochemical model (VULCAN) on an equatorial plane. We find a fairly uniform composition distribution below 0.1 mbar controlled by the dayside, as a result of slow chemical evolution. NH$_3$ and CH$_3$OH are concluded to be the most unambiguous proxies to infer surfaces on sub-Neptunes in the era of the James Webb Space Telescope (JWST). △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJL

A multispecies pseudoadiabat for simulating condensable-rich exoplanet atmospheres

Authors: R. J. Graham, Tim Lichtenberg, Ryan Boukrouche, Ray Pierrehumbert

Abstract: Central stages in the evolution of rocky, potentially habitable planets may play out under atmospheric conditions with a large inventory of non-dilute condensable components. Variations in condensate retention and accompanying changes in local lapse rate may substantially affect planetary climate and surface conditions, but there is currently no general theory to effectively describe such atmosphe… ▽ More Central stages in the evolution of rocky, potentially habitable planets may play out under atmospheric conditions with a large inventory of non-dilute condensable components. Variations in condensate retention and accompanying changes in local lapse rate may substantially affect planetary climate and surface conditions, but there is currently no general theory to effectively describe such atmospheres. In this article, expanding on the work by Li et al. (2018), we generalize the single-component moist pseudoadiabat derivation in Pierrehumbert (2010) to allow for multiple condensing components of arbitrary diluteness and retained condensate fraction. The introduction of a freely tunable retained condensate fraction allows for a flexible, self-consistent treatment of atmospheres with non-dilute condensable components. To test the pseudoadiabat's capabilities for simulating a diverse range of climates, we apply the formula to planetary atmospheres with compositions, surface pressures, and temperatures representing important stages with condensable-rich atmospheres in the evolution of terrestrial planets: a magma ocean planet in a runaway greenhouse state; a post-impact, late veneer-analogue planet with a complex atmospheric composition; and an Archean Earth-like planet near the outer edge of the classical circumstellar habitable zone. We find that variations in the retention of multiple non-dilute condensable species can significantly affect the lapse rate and in turn outgoing radiation and the spectral signatures of planetary atmospheres. The presented formulation allows for a more comprehensive treatment of the climate evolution of rocky exoplanets and early Earth analogues. △ Less

Submitted 29 August, 2021; originally announced August 2021.

Comments: Accepted for publication in Planetary Science Journal on 8/24/2021. For a video summary see https://bit.ly/3ztVhfF

Journal ref: The Planetary Science Journal (2021), 2, 207

Beyond runaway: initiation of the post-runaway greenhouse state on rocky exoplanets

Authors: Ryan Boukrouche, Tim Lichtenberg, Raymond T. Pierrehumbert

Abstract: The runaway greenhouse represents the ultimate climate catastrophe for rocky, Earth-like worlds: when the incoming stellar flux cannot be balanced by radiation to space, the oceans evaporate and exacerbate heating, turning the planet into a hot wasteland with a steam atmosphere overlying a possibly molten magma surface. The equilibrium state beyond the runaway greenhouse instellation limit depends… ▽ More The runaway greenhouse represents the ultimate climate catastrophe for rocky, Earth-like worlds: when the incoming stellar flux cannot be balanced by radiation to space, the oceans evaporate and exacerbate heating, turning the planet into a hot wasteland with a steam atmosphere overlying a possibly molten magma surface. The equilibrium state beyond the runaway greenhouse instellation limit depends on the radiative properties of the atmosphere and its temperature structure. Here, we use 1-D radiative-convective models of steam atmospheres to explore the transition from the tropospheric radiation limit to the post-runaway climate state. To facilitate eventual simulations with 3-D global circulation models, a computationally efficient band-grey model is developed, which is capable of reproducing the key features of the more comprehensive calculations. We analyze two factors which determine the equilibrated surface temperature of post-runaway planets. The infrared cooling of the planet is strongly enhanced by the penetration of the dry adiabat into the optically thin upper regions of the atmosphere. In addition, thermal emission of both shortwave and near-IR fluxes from the hot lower atmospheric layers, which can radiate through window regions of the spectrum, is quantified. Astronomical surveys of rocky exoplanets in the runaway greenhouse state may discriminate these features using multi-wavelength observations. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ; 19 pages, 11 figures; Summary video available at https://exoplanet-talks.org/talk/379

Journal ref: The Astrophysical Journal (2021), 919, 130

Simulating gas giant exoplanet atmospheres with Exo-FMS: Comparing semi-grey, picket fence and correlated-k radiative-transfer schemes

Authors: Elspeth K. H. Lee, Vivien Parmentier, Mark Hammond, Simon L. Grimm, Daniel Kitzmann, Xianyu Tan, Shang-Min Tsai, Raymond T. Pierrehumbert

Abstract: Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM us… ▽ More Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM using these RT schemes to hot Jupiter simulation results from the literature. We perform a HD 209458b-like simulation with the three schemes and compare their results. These simulations are then post-processed to compare their observable differences. The semi-grey scheme results show qualitative agreement with previous studies in line with variations seen between GCM models. The real gas model reproduces well the temperature and dynamical structures from other studies. After post-processing our non-grey picket fence scheme compares very favourably with the real gas model, producing similar transmission spectra, emission spectra and phase curve behaviours. Exo-FMS is able to reliably reproduce the essential features of contemporary GCM models in the hot gas giant regime. Our results suggest the picket fence approach offers a simple way to improve upon RT realism beyond semi-grey schemes. △ Less

Submitted 28 June, 2021; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: MNRAS accepted 22 June 2021 - V2, typos fixed

3D convection-resolving model of temperate, tidally-locked exoplanets

Authors: Maxence Lefèvre, Martin Turbet, Raymond Pierrehumbert

Abstract: A large fraction of known terrestrial-size exoplanets located in the Habitable Zone of M-dwarfs are expected to be tidally-locked. Numerous efforts have been conducted to study the climate of such planets, using in particular 3-D Global Climate Models (GCM). One of the biggest challenges in simulating such an extreme environment is to properly represent the effects of sub-grid convection. Most GCM… ▽ More A large fraction of known terrestrial-size exoplanets located in the Habitable Zone of M-dwarfs are expected to be tidally-locked. Numerous efforts have been conducted to study the climate of such planets, using in particular 3-D Global Climate Models (GCM). One of the biggest challenges in simulating such an extreme environment is to properly represent the effects of sub-grid convection. Most GCMs use either a simplistic convective-adjustment parametrization or sophisticated (e.g., mass flux scheme) Earth-tuned parametrizations. One way to improve the representation of convection is to study convection using Convection Resolving numerical Models (CRMs), with an fine spatial resolution . In this study, we developed a CRM coupling the non-hydrostatic dynamical core WRF with the radiative transfer and cloud/precipitation models of the LMD-Generic climate model to study convection and clouds on tidally-locked planets, with a focus on Proxima b. Simulations were performed for a set of 3 surface temperatures (corresponding to three different incident fluxes) and 2 rotation rates, assuming an Earth-like atmosphere. The main result of our study is that while we recover the prediction of GCMs that (low-altitude) cloud albedo increases with increasing stellar flux, the cloud feedback is much weaker due to transient aggregation of convection leading to low partial cloud cover. △ Less

Submitted 12 April, 2021; originally announced April 2021.

Comments: Accepted for publications in The Astrophysical Journal. 42 pages, 14 figures, 1 table

Evidence for disequilibrium chemistry from vertical mixing in hot Jupiter atmospheres. A comprehensive survey of transiting close-in gas giant exoplanets with warm-Spitzer/IRAC

Authors: C. Baxter, J-M. Désert, S-M. Tsai, K. O. Todorov, J. L. Bean, D. Deming, V. Parmentier, J. J. Fortney, M. Line, D. Thorngren, R. T. Pierrehumbert, A. Burrows, A. P. Showman

Abstract: [Abridged] Aims. We present a large atmospheric study of 49 gas giant exoplanets using infrared transmission photometry with Spitzer/IRAC at 3.6 and 4.5um. Methods. We uniformly analyze 70 photometric light curves of 33 transiting planets using our custom pipeline, which implements pixel level decorrelation. We use this survey to understand how infrared photometry traces changes in atmospheric che… ▽ More [Abridged] Aims. We present a large atmospheric study of 49 gas giant exoplanets using infrared transmission photometry with Spitzer/IRAC at 3.6 and 4.5um. Methods. We uniformly analyze 70 photometric light curves of 33 transiting planets using our custom pipeline, which implements pixel level decorrelation. We use this survey to understand how infrared photometry traces changes in atmospheric chemical properties as a function of planetary temperature. We compare our measurements to a grid of 1D radiative-convective equilibrium forward atmospheric models which include disequilibrium chemistry. We explore various strengths of vertical mixing (Kzz = 0 - 10^12 cm2/s) as well as two chemical compositions (1x and 30x solar). Results. We find that, on average, Spitzer probes a difference of 0.5 atmospheric scale heights between 3.6 and 4.5um, which is measured at 7.5sigma level of significance. We find that the coolest planets show a lack of methane compared to expectations, which has also been reported by previous studies of individual objects. We show that the sample of coolest planets rule out 1x solar composition with >3sigma confidence while supporting low vertical mixing (Kzz = 10^8 cm2/s). On the other hand, we find that the hot planets are best explained by models with 1x solar metallicity and high vertical mixing (Kzz = 10^12 cm2/s). We interpret this as the lofting of CH4 to the upper atmospheric layers. Changing the interior temperature changes the expectation for equilibrium chemistry in deep layers, hence the expectation of disequilibrium chemistry higher up. We also find a significant scatter in the transmission signatures of the mid-temperate and ultra-hot planets, likely due to increased atmospheric diversity, without the need to invoke higher metallicities. Additionally, we compare Spitzer transmission with emission for the same planets and find no evidence for correlation. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 43 pages, 17 Figures. Accepted on 9 Feb 2021 in Astronomy & Astrophysics

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.

Comments: Accepted for publication in JGR Planets; 23 pages, 10 figures; summaries available at https://bit.ly/MO21_blog (blog) and http://bit.ly/MO21_video (video)

Journal ref: Journal of Geophysical Research: Planets (2021), Exoplanets: The Nexus of Astronomy and Geoscience

The Phase-curve Signature of Condensible Water-rich Atmospheres on Slowly Rotating Tidally Locked Exoplanets

Authors: Feng Ding, Raymond T. Pierrehumbert

Abstract: We use an idealized three-dimensional general circulation model to study condensible-rich atmospheres with an ineffective cold trap on slowly rotating tidally locked terrestrial planets. In particular, we show the climate dynamics in a thin and temperate atmosphere with condensible water vapor. The similarities between our thin and temperate atmosphere and the warm and thick atmosphere approaching… ▽ More We use an idealized three-dimensional general circulation model to study condensible-rich atmospheres with an ineffective cold trap on slowly rotating tidally locked terrestrial planets. In particular, we show the climate dynamics in a thin and temperate atmosphere with condensible water vapor. The similarities between our thin and temperate atmosphere and the warm and thick atmosphere approaching the water vapor runaway greenhouse in previous works are discussed, including the reversal of the thermal emission between the day and night hemispheres. Different from the transit spectroscopy of water vapor that depends on the absolute amount of atmospheric water vapor, the contrast between the dayside and nightside thermal emission provides information regarding the relative ratio of water vapor to the background atmosphere as well as the atmospheric pressure near the substellar tropopause and the emission level on the nightside on potentially habitable worlds. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: 15 pages, 4 figures

Journal ref: ApJL, 901:L33, 2020 October

The Equatorial Jet Speed on Tidally Locked Planets: I -- Terrestrial Planets

Authors: Mark Hammond, Shang-Min Tsai, Raymond T. Pierrehumbert

Abstract: The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and… ▽ More The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and the vertical motion of the atmosphere. The three-dimensional structure of the zonal acceleration is vital to this mechanism. We demonstrate this mechanism in a hierarchy of models. We calculate the three-dimensional stationary waves induced by the forcing on these planets, and show the vertical structure of the zonal acceleration produced by these waves, which we use to suggest a mechanism for how the jet forms. GCM simulations are used to confirm the equilibrium state predicted by this mechanism, where the acceleration from these waves is balanced by an interaction between the zonal-mean vertical velocity and the jet. We derive a simple model of this using the "Weak Temperature Gradient" approximation, which gives an estimate of the jet speed on a terrestrial tidally locked planet. We conclude that the proposed mechanism is a good description of the formation of an equatorial jet on a terrestrial tidally locked planet, and should be useful for interpreting observations and simulations of these planets. The mechanism requires assumptions such as a large equatorial Rossby radius and weak acceleration due to transient waves, and a different mechanism may produce the equatorial jets on gaseous tidally locked planets. △ Less

Submitted 1 September, 2020; originally announced September 2020.

Comments: Accepted by ApJ

Thermodynamic and Energetic Limits on Continental Silicate Weathering Strongly Impact the Climate and Habitability of Wet, Rocky Worlds

Authors: R. J. Graham, R. T. Pierrehumbert

Abstract: The "liquid water habitable zone" (HZ) concept is predicated on the ability of the silicate weathering feedback to stabilize climate across a wide range of instellations. However, representations of silicate weathering used in current estimates of the effective outer edge of the HZ do not account for the thermodynamic limit on concentration of weathering products in runoff set by clay precipitatio… ▽ More The "liquid water habitable zone" (HZ) concept is predicated on the ability of the silicate weathering feedback to stabilize climate across a wide range of instellations. However, representations of silicate weathering used in current estimates of the effective outer edge of the HZ do not account for the thermodynamic limit on concentration of weathering products in runoff set by clay precipitation, nor for the energetic limit on precipitation set by planetary instellation. We find that when the thermodynamic limit is included in an idealized coupled climate/weathering model, steady-state planetary climate loses sensitivity to silicate dissolution kinetics, becoming sensitive to temperature primarily through the effect of temperature on runoff and to pCO$_2$ through an effect on solute concentration mediated by pH. This increases sensitivity to land fraction, CO$_2$ outgassing, and geological factors such as soil age and lithology, all of which are found to have a profound effect on the position of the effective outer edge of the HZ. The interplay between runoff sensitivity and the energetic limit on precipitation leads to novel warm states in the outer reaches of the HZ, owing to the decoupling of temperature and precipitation. We discuss strategies for detecting the signature of silicate weathering feedback through exoplanet observations in light of insights derived from the revised picture of weathering. △ Less

Submitted 20 May, 2020; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: 27 pages excluding references, 30 pages total, 8 figures. Submitted to ApJ, accepted 15 May 2020. (reuploaded to arxiv 5/20/2020 to include corresponding author's email address and news of acceptance)

Simplified 3D GCM modelling of the irradiated brown dwarf WD0137-349B

Authors: Elspeth K. H. Lee, Sarah L. Casewell, Katy L. Chubb, Mark Hammond, Xianyu Tan, Shang-Min Tsai, Raymond T. Pierrehumbert

Abstract: Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical… ▽ More Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical properties of the Brown dwarf WD0137-349B. Methods: We use the three dimensional GCM model Exo-FMS, with a dual-band grey radiative-transfer scheme to model the atmosphere of WD0137-349B. The results of the GCM model are post-processed using the three dimensional Monte Carlo radiative-transfer model \textsc{cmcrt}. Results: Our results suggest inefficient day-night energy transport and a large day-night temperature contrast for WD0137-349B. Multiple flow patterns are present, shifting energy asymmetrically eastward or westward depending on their zonal direction and latitude. Regions of overturning are produced on the western terminator. We are able to reproduce the start of the system near-IR emission excess at $\gtrsim$ 1.95 $μ$m as observed by the GNIRS instrument. Our model over predicts the IR phase curve fluxes by factors of $\approx$1-3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using \textsc{vulcan} suggests a highly ionised region at high altitudes can form on the dayside of the Brown dwarf. Conclusions: We present a first attempt at simulating the atmosphere of a short period White dwarf - Brown dwarf binary in a 3D setting. Further studies into the radiative and photochemical heating from the UV irradiation is required to more accurately capture the energy balance inside the Brown dwarf atmosphere. Cloud formation may also play an important role in shaping the emission spectra of the Brown dwarf. △ Less

Submitted 25 June, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: Version 2 - Accepted June 19 2020 (MNRAS)

The Need for Laboratory Measurements and Ab Initio Studies to Aid Understanding of Exoplanetary Atmospheres

Authors: Jonathan J. Fortney, Tyler D. Robinson, Shawn Domagal-Goldman, Anthony D. Del Genio, Iouli E. Gordon, Ehsan Gharib-Nezhad, Nikole Lewis, Clara Sousa-Silva, Vladimir Airapetian, Brian Drouin, Robert J. Hargreaves, Xinchuan Huang, Tijs Karman, Ramses M. Ramirez, Gregory B. Rieker, Jonathan Tennyson, Robin Wordsworth, Sergei N Yurchenko, Alexandria V Johnson, Timothy J. Lee, Chuanfei Dong, Stephen Kane, Mercedes Lopez-Morales, Thomas Fauchez, Timothy Lee , et al. (63 additional authors not shown)

Abstract: We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as c… ▽ More We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Submitted as an Astro2020 Science White Paper

Wave-Mean Flow Interactions in the Atmospheric Circulation of Tidally Locked Planets

Authors: Mark Hammond, Raymond T. Pierrehumbert

Abstract: We use a linear shallow-water model to investigate the global circulation of the atmospheres of tidally locked planets. Simulations, observations, and simple models show that if these planets are sufficiently rapidly rotating, their atmospheres have an eastward equatorial jet and a hot-spot east of the substellar point. We linearize the shallow-water model about this eastward flow and its associat… ▽ More We use a linear shallow-water model to investigate the global circulation of the atmospheres of tidally locked planets. Simulations, observations, and simple models show that if these planets are sufficiently rapidly rotating, their atmospheres have an eastward equatorial jet and a hot-spot east of the substellar point. We linearize the shallow-water model about this eastward flow and its associated height perturbation. The forced solutions of this system show that the shear flow explains the form of the global circulation, particularly the hot-spot shift and the positions of the cold standing waves on the night-side. We suggest that the eastward hot-spot shift seen in observations and 3D simulations of these atmospheres is caused by the zonal flow Doppler-shifting the stationary wave response eastwards, summed with the height perturbation from the flow itself. This differs from other studies which explained the hot-spot shift as pure advection of heat from air flowing eastward from the substellar point, or as equatorial waves travelling eastwards. We compare our solutions to simulations in our climate model Exo-FMS, and show that the height fields and wind patterns match. We discuss how planetary properties affect the global circulation, and how they change observables such as the hot-spot shift or day-night contrast. We conclude that the wave-mean flow interaction between the stationary planetary waves and the equatorial jet is a vital part of the equilibrium circulation on tidally locked planets. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: Accepted by ApJ

Global or Local Pure-Condensible Atmospheres: Importance of Horizontal Latent Heat Transport

Authors: Feng Ding, Raymond T. Pierrehumbert

Abstract: The distribution of a pure condensible planetary atmosphere in equilibrium with a surface reservoir is revisited employing the energy budget of the climate system, emphasizing the atmospheric horizontal latent heat transport. This configuration is applicable to icy Solar System bodies such as Triton as well as a range of possible exoplanet atmospheres, including water or $\mathrm{CO_2}$ iceballs o… ▽ More The distribution of a pure condensible planetary atmosphere in equilibrium with a surface reservoir is revisited employing the energy budget of the climate system, emphasizing the atmospheric horizontal latent heat transport. This configuration is applicable to icy Solar System bodies such as Triton as well as a range of possible exoplanet atmospheres, including water or $\mathrm{CO_2}$ iceballs or ocean worlds, and lava planets with mineral vapor atmospheres. Climate regimes for slowly rotating planets with the hot-spot near the substellar point, and for rapidly rotating planets with a warm equatorial belt, are both treated. A non-dimensional parameter controlling the fractional variation of the surface pressure is derived; it measures whether the pure condensible atmosphere is global or localized. The global pure condensible atmosphere with the non-dimensional parameter much less than order of unity is maintained by the strong horizontal latent heat transport associated with an "evaporation/sublimation-driven flow" from warm to cold places that compensates for the incoming differential radiative forcing. We show that the variation of surface temperature can be estimated in terms of this non-dimensional parameter if it is not too large. In the case of a pure water-vapor atmosphere with an ice or liquid surface, we show that the atmosphere is thick enough to maintain nearly isothermal surface conditions even when the substellar surface temperature is around the freezing point. Finally, it is proposed that the evaporation/sublimation-driven flow regime for global atmospheres could be detected via its effect on the inhomogeneous distribution of minor non-condensible components in the atmosphere. △ Less

Submitted 28 September, 2018; originally announced September 2018.

Comments: accepted for publication in ApJ

Exploring the Atmosphere of Neoproterozoic Earth: The Effect of O$_{2}$ on Haze Formation and Composition

Authors: Sarah M. Hörst, Chao He, Melissa S. Ugelow, A. Mark Jellinek, Raymond T. Pierrehumbert, Margaret A. Tolbert

Abstract: Previous studies of haze formation in the atmosphere of the Early Earth have focused on N$_{2}$/CO$_{2}$/CH$_{4}$ atmospheres. Here, we experimentally investigate the effect of O$_{2}$ on the formation and composition of aerosols to improve our understanding of haze formation on the Neoproterozoic Earth. We obtained in situ size, particle density, and composition measurements of aerosol particles… ▽ More Previous studies of haze formation in the atmosphere of the Early Earth have focused on N$_{2}$/CO$_{2}$/CH$_{4}$ atmospheres. Here, we experimentally investigate the effect of O$_{2}$ on the formation and composition of aerosols to improve our understanding of haze formation on the Neoproterozoic Earth. We obtained in situ size, particle density, and composition measurements of aerosol particles produced from N$_{2}$/CO$_{2}$/CH$_{4}$/O$_{2}$ gas mixtures subjected to FUV radiation (115-400 nm) for a range of initial CO$_{2}$/CH$_{4}$/O$_{2}$ mixing ratios (O$_{2}$ ranging from 2 ppm to 0.2\%). At the lowest O$_{2}$ concentration (2 ppm), the addition increased particle production for all but one gas mixture. At higher oxygen concentrations (20 ppm and greater) particles are still produced, but the addition of O$_{2}$ decreases the production rate. Both the particle size and number density decrease with increasing O$_{2}$, indicating that O$_{2}$ affects particle nucleation and growth. The particle density increases with increasing O$_{2}$. The addition of CO$_{2}$ and O$_{2}$ not only increases the amount of oxygen in the aerosol, but it also increases the degree of nitrogen incorporation. In particular, the addition of O$_{2}$ results in the formation of nitrate bearing molecules. The fact that the presence of oxygen bearing molecules increases the efficiency of nitrogen fixation has implications for the role of haze as a source of molecules required for the origin and evolution of life. The composition changes also likely affect the absorption and scattering behavior of these particles but optical properties measurements are required to fully understand the implications for the effect on the planetary radiative energy balance and climate. △ Less

Submitted 8 April, 2018; originally announced April 2018.

Comments: 15 pages, 3 tables, 8 figures, accepted in Astrophysical Journal

Journal ref: Astrophysical Journal, 858:199, 2018

Exoplanet Diversity in the Era of Space-based Direct Imaging Missions

Authors: Ravi Kopparapu, Eric Hebrard, Rus Belikov, Natalie M. Batalha, Gijs D. Mulders, Chris Stark, Dillon Teal, Shawn Domagal-Goldman, Dawn Gelino, Avi Mandell, Aki Roberge, Stephen Rinehart, Stephen R. Kane, Yasuhiro Hasegawa, Wade Henning, Brian Hicks, Vardan Adibekyan, Edward W. Schwieterman, Erika Kohler, Johanna Teske, Natalie Hinkel, Conor Nixon, Kevin France, William Danchi, Jacob Haqq-Misra , et al. (33 additional authors not shown)

Abstract: This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This wil… ▽ More This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux. △ Less

Submitted 10 March, 2018; originally announced March 2018.

Comments: A white paper submitted to the National Academy of Sciences Exoplanet Science Strategy

Linking the Climate and Thermal Phase Curve of 55 Cancri e

Authors: Mark Hammond, Raymond Pierrehumbert

Abstract: The thermal phase curve of 55 Cancri e is the first measurement of the temperature distribution of a tidally locked Super-Earth, but raises a number of puzzling questions about the planet's climate. The phase curve has a high amplitude and peak offset, suggesting that it has a significant eastward hot-spot shift as well as a large day-night temperature contrast. We use a general circulation model… ▽ More The thermal phase curve of 55 Cancri e is the first measurement of the temperature distribution of a tidally locked Super-Earth, but raises a number of puzzling questions about the planet's climate. The phase curve has a high amplitude and peak offset, suggesting that it has a significant eastward hot-spot shift as well as a large day-night temperature contrast. We use a general circulation model to model potential climates, and investigate the relation between bulk atmospheric composition and the magnitude of these seemingly contradictory features. We confirm theoretical models of tidally locked circulation are consistent with our numerical model of 55 Cnc e, and rule out certain atmospheric compositions based on their thermodynamic properties. Our best-fitting atmosphere has a significant hot-spot shift and day-night contrast, although these are not as large as the observed phase curve. We discuss possible physical processes which could explain the observations, and show that night-side cloud formation from species such as SiO from a day-side magma ocean could potentially increase the phase curve amplitude and explain the observations. We conclude that the observations could be explained by an optically thick atmosphere with a low mean molecular weight, a surface pressure of several bar and a strong eastward circulation, with night-side cloud formation a possible explanation for the difference between our model and the observations. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

Dynamics of atmospheres with a non-dilute condensible component

Authors: Raymond T. Pierrehumbert, Feng Ding

Abstract: The diversity of characteristics for the host of recently discovered exoplanets opens up a great deal of fertile new territory for geophysical fluid dynamics, particularly when the fluid flow is coupled to novel thermodynamics, radiative transfer or chemistry. In this paper, we survey one of these new areas-the climate dynamics of atmospheres with a non-dilute condensible component, defined as the… ▽ More The diversity of characteristics for the host of recently discovered exoplanets opens up a great deal of fertile new territory for geophysical fluid dynamics, particularly when the fluid flow is coupled to novel thermodynamics, radiative transfer or chemistry. In this paper, we survey one of these new areas-the climate dynamics of atmospheres with a non-dilute condensible component, defined as the situation in which a condensible component of the atmosphere makes up a substantial fraction of the atmospheric mass within some layer. Non-dilute dynamics can occur for a wide range of condensibles, generically applying near both the inner and the outer edges of the conventional habitable zone and in connection with runaway greenhouse phenomena. It also applies in a wide variety of other planetary circumstances. We first present a number of analytical results developing some key features of non-dilute atmospheres, and then show how some of these features are manifest in simulations with a general circulation model adapted to handle non-dilute atmospheres. We find that non-dilute atmospheres have weak horizontal temperature gradients even for rapidly rotating planets, and that their circulations are largely barotropic. The relative humidity of the condensible component tends towards 100% as the atmosphere becomes more non-dilute, which has important implications for runaway greenhouse thresholds. Non-dilute atmospheres exhibit a number of interesting organized convection features, for which there is not yet any adequate theoretical understanding. △ Less

Submitted 10 August, 2016; originally announced August 2016.

Comments: Accepted for publication in Proc. R. Soc. A

Convection in Condensible-rich Atmospheres

Authors: Feng Ding, Raymond T. Pierrehumbert

Abstract: Condensible substances are nearly ubiquitous in planetary atmospheres. For the most familiar case-water vapor in Earth's present climate-the condensible gas is dilute, in the sense that its concentration is everywhere small relative to the noncondensible background gases. A wide variety of important planetary climate problems involve nondilute condensible substances. These include planets near or… ▽ More Condensible substances are nearly ubiquitous in planetary atmospheres. For the most familiar case-water vapor in Earth's present climate-the condensible gas is dilute, in the sense that its concentration is everywhere small relative to the noncondensible background gases. A wide variety of important planetary climate problems involve nondilute condensible substances. These include planets near or undergoing a water vapor runaway and planets near the outer edge of the conventional habitable zone, for which CO2 is the condensible. Standard representations of convection in climate models rely on several approximations appropriate only to the dilute limit, while nondilute convection differs in fundamental ways from dilute convection. In this paper, a simple parameterization of convection valid in the nondilute as well as dilute limits is derived and used to discuss the basic character of nondilute convection. The energy conservation properties of the scheme are discussed in detail and are verified in radiative-convective simulations. As a further illustration of the behavior of the scheme, results for a runaway greenhouse atmosphere for both steady instellation and seasonally varying instellation corresponding to a highly eccentric orbit are presented. The latter case illustrates that the high thermal inertia associated with latent heat in nondilute atmospheres can damp out the effects of even extreme seasonal forcing. △ Less

Submitted 10 August, 2016; originally announced August 2016.

Comments: Accepted for publication in ApJ with minor revisions

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Authors: Jonathan J. Fortney, Tyler D. Robinson, Shawn Domagal-Goldman, David Skålid Amundsen, Matteo Brogi, Mark Claire, David Crisp, Eric Hebrard, Hiroshi Imanaka, Remco de Kok, Mark S. Marley, Dillon Teal, Travis Barman, Peter Bernath, Adam Burrows, David Charbonneau, Richard S. Freedman, Dawn Gelino, Christiane Helling, Kevin Heng, Adam G. Jensen, Stephen Kane, Eliza M. -R. Kempton, Ravi Kumar Kopparapu, Nikole K. Lewis , et al. (9 additional authors not shown)

Abstract: Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements… ▽ More Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr) △ Less

Submitted 23 February, 2016; v1 submitted 19 February, 2016; originally announced February 2016.

Comments: This is a white paper originally initiated within NASA's Nexus for Exoplanet System Science (NExSS), with additional community input. Community comment is welcome. Please submit comments and suggestions by going to the Google Doc version of the paper: https://docs.google.com/document/d/1109uy_5n9VWuOOn-huTVv-gJgBFRO9QH0Ojkao1a6Lk/edit?usp=sharing

Comparison of "warm and wet" and "cold and icy" scenarios for early Mars in a 3D climate model

Authors: Robin D. Wordsworth, Laura Kerber, Raymond T. Pierrehumbert, Francois Forget, James W. Head

Abstract: We use a 3D general circulation model to compare the primitive Martian hydrological cycle in "warm and wet" and "cold and icy" scenarios. In the warm and wet scenario, an anomalously high solar flux or intense greenhouse warming artificially added to the climate model are required to maintain warm conditions and an ice-free northern ocean. Precipitation shows strong surface variations, with high r… ▽ More We use a 3D general circulation model to compare the primitive Martian hydrological cycle in "warm and wet" and "cold and icy" scenarios. In the warm and wet scenario, an anomalously high solar flux or intense greenhouse warming artificially added to the climate model are required to maintain warm conditions and an ice-free northern ocean. Precipitation shows strong surface variations, with high rates around Hellas basin and west of Tharsis but low rates around Margaritifer Sinus (where the observed valley network drainage density is nonetheless high). In the cold and icy scenario, snow migration is a function of both obliquity and surface pressure, and limited episodic melting is possible through combinations of seasonal, volcanic and impact forcing. At surface pressures above those required to avoid atmospheric collapse (~0.5 bar) and moderate to high obliquity, snow is transported to the equatorial highland regions where the concentration of valley networks is highest. Snow accumulation in the Aeolis quadrangle is high, indicating an ice-free northern ocean is not required to supply water to Gale crater. At lower surface pressures and obliquities, both H2O and CO2 are trapped as ice at the poles and the equatorial regions become extremely dry. The valley network distribution is positively correlated with snow accumulation produced by the cold and icy simulation at 41.8 degrees obliquity but uncorrelated with precipitation produced by the warm and wet simulation. Because our simulations make specific predictions for precipitation patterns under different climate scenarios, they motivate future targeted geological studies. △ Less

Submitted 15 June, 2015; originally announced June 2015.

Comments: 20 pages, 14 figures, 2 tables. In press at JGR (Planets)

Abiotic oxygen-dominated atmospheres on terrestrial habitable zone planets

Authors: Robin Wordsworth, Raymond Pierrehumbert

Abstract: Detection of life on other planets requires identification of biosignatures, i.e., observable planetary properties that robustly indicate the presence of a biosphere. One of the most widely accepted biosignatures for an Earth-like planet is an atmosphere where oxygen is a major constituent. Here we show that lifeless habitable zone terrestrial planets around any star type may develop oxygen-domina… ▽ More Detection of life on other planets requires identification of biosignatures, i.e., observable planetary properties that robustly indicate the presence of a biosphere. One of the most widely accepted biosignatures for an Earth-like planet is an atmosphere where oxygen is a major constituent. Here we show that lifeless habitable zone terrestrial planets around any star type may develop oxygen-dominated atmospheres as a result of water photolysis, because the cold trap mechanism that protects H2O on Earth is ineffective when the atmospheric inventory of non-condensing gases (e.g., N2, Ar) is low. Hence the spectral features of O2 and O3 alone cannot be regarded as robust signs of extraterrestrial life. △ Less

Submitted 20 March, 2014; v1 submitted 11 March, 2014; originally announced March 2014.

Comments: 5 pages, 3 figures, accepted for publication in ApJ Letters. Minor changes from previous version

Water loss from terrestrial planets with CO2-rich atmospheres

Authors: Robin Wordsworth, Raymond Pierrehumbert

Abstract: Water photolysis and hydrogen loss from the upper atmospheres of terrestrial planets is of fundamental importance to climate evolution but remains poorly understood in general. Here we present a range of calculations we performed to study the dependence of water loss rates from terrestrial planets on a range of atmospheric and external parameters. We show that CO2 can only cause significant water… ▽ More Water photolysis and hydrogen loss from the upper atmospheres of terrestrial planets is of fundamental importance to climate evolution but remains poorly understood in general. Here we present a range of calculations we performed to study the dependence of water loss rates from terrestrial planets on a range of atmospheric and external parameters. We show that CO2 can only cause significant water loss by increasing surface temperatures over a narrow range of conditions, with cooling of the middle and upper atmosphere acting as a bottleneck on escape in other circumstances. Around G-stars, efficient loss only occurs on planets with intermediate CO2 atmospheric partial pressures (0.1 to 1 bar) that receive a net flux close to the critical runaway greenhouse limit. Because G-star total luminosity increases with time but XUV/UV luminosity decreases, this places strong limits on water loss for planets like Earth. In contrast, for a CO2-rich early Venus, diffusion limits on water loss are only important if clouds caused strong cooling, implying that scenarios where the planet never had surface liquid water are indeed plausible. Around M-stars, water loss is primarily a function of orbital distance, with planets that absorb less flux than ~270 W/m^2 (global mean) unlikely to lose more than one Earth ocean of H2O over their lifetimes unless they lose all their atmospheric N2/CO2 early on. Because of the variability of H2O delivery during accretion, our results suggest that many 'Earth-like' exoplanets in the habitable zone may have ocean-covered surfaces, stable CO2/H2O-rich atmospheres, and high mean surface temperatures. △ Less

Submitted 12 October, 2013; v1 submitted 13 June, 2013; originally announced June 2013.

Comments: Accepted for publication in ApJ with minor revisions

The Effect of Host Star Spectral Energy Distribution and Ice-Albedo Feedback on the Climate of Extrasolar Planets

Authors: Aomawa L. Shields, Victoria S. Meadows, Cecilia M. Bitz, Raymond T. Pierrehumbert, Manoj M. Joshi, Tyler D. Robinson

Abstract: Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. Here we explore this effect using a one dimensional (1-D), line-by-line, radiative-transfer model to calculate broadband planetary albedos as input to a seasonally varying, 1-D energy-balance climate model. A three-dimensional general circul… ▽ More Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. Here we explore this effect using a one dimensional (1-D), line-by-line, radiative-transfer model to calculate broadband planetary albedos as input to a seasonally varying, 1-D energy-balance climate model. A three-dimensional general circulation model is also used to explore the atmosphere's response to changes in incoming stellar radiation, or instellation, and surface albedo. Using this hierarchy of models we simulate planets covered by ocean, land, and water ice of varying grain size, with incident radiation from stars of different spectral types. Terrestrial planets orbiting stars with higher near-UV radiation exhibit a stronger ice-albedo feedback. We find that ice-covered conditions occur on an F-dwarf planet with only a 2% reduction in instellation relative to the present instellation on Earth, assuming fixed CO2 (present atmospheric level on Earth). A similar planet orbiting the Sun at an equivalent flux distance requires an 8% reduction in instellation, while a planet orbiting an M-dwarf star requires an additional 19% reduction in instellation to become ice-covered, equivalent to 73% of the modern solar constant. The surface ice-albedo feedback effect becomes less important at the outer edge of the habitable zone, where atmospheric CO2 can be expected to be high in order to maintain clement conditions for surface liquid water. We show that 3-10 bars of CO2 will entirely mask the climatic effect of ice and snow, leaving the outer limits of the habitable zone unaffected by the spectral dependence of water ice and snow albedo. However, less CO2 is needed to maintain open water for a planet orbiting an M-dwarf star, than would be the case for hotter main-sequence stars. △ Less

Submitted 29 May, 2013; originally announced May 2013.

Comments: 63 pages, 2 tables, 20 figures; accepted for publication in Astrobiology

Hydrogen Greenhouse Planets Beyond the Habitable Zone

Authors: Raymond Pierrehumbert, Eric Gaidos

Abstract: We show that collision-induced absorption allows molecular hydrogen to act as an incondensible greenhouse gas, and that bars or tens of bars of primordial H2-He mixtures can maintain surface temperatures above the freezing point of water well beyond the "classical" habitable zone defined for CO2 greenhouse atmospheres. Using a 1-D radiative-convective model we find that 40 bars of pure H2 on a 3 E… ▽ More We show that collision-induced absorption allows molecular hydrogen to act as an incondensible greenhouse gas, and that bars or tens of bars of primordial H2-He mixtures can maintain surface temperatures above the freezing point of water well beyond the "classical" habitable zone defined for CO2 greenhouse atmospheres. Using a 1-D radiative-convective model we find that 40 bars of pure H2 on a 3 Earth-mass planet can maintain a surface temperature of 280K out to 1.5AU from an early-type M dwarf star and 10 AU from a G-type star. Neglecting the effects of clouds and of gaseous absorbers besides H2, the flux at the surface would be sufficient for photosynthesis by cyanobacteria (in the G star case) or anoxygenic phototrophs (in the M star case). We argue that primordial atmospheres of one to several hundred bars of H2-He are possible, and use a model of hydrogen escape to show that such atmospheres are likely to persist further than 1.5 AU from M stars, and 2 AU from G stars, assuming these planets have protecting magnetic fields. We predict that the microlensing planet OGLE-05-390L could have retained a H2-He atmosphere and be habitable at ~2.6 AU from its host M star. △ Less

Submitted 29 April, 2011; originally announced May 2011.

Comments: Accepted to Astrophysical Journal Letters