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