A self-consistent cloud model for brown dwarfs and young giant exoplanets: comparison with photometric and spectroscopic observations

A self-consistent cloud model for brown dwarfs and young giant exoplanets: comparison with photometric and spectroscopic observations

Authors:

Charnay et al

Abstract:

We developed a simple, physical and self-consistent cloud model for brown dwarfs and young giant exoplanets. We compared different parametrisations for the cloud particle size, by either fixing particle radii, or fixing the mixing efficiency (parameter fsed) or estimating particle radii from simple microphysics. The cloud scheme with simple microphysics appears as the best parametrisation by successfully reproducing the observed photometry and spectra of brown dwarfs and young giant exoplanets. In particular, it reproduces the L-T transition, due to the condensation of silicate and iron clouds below the visible/near-IR photosphere. It also reproduces the reddening observed for low-gravity objects, due to an increase of cloud optical depth for low gravity. In addition, we found that the cloud greenhouse effect shifts chemical equilibriums, increasing the abundances of species stable at high temperature. This effect should significantly contribute to the strong variation of methane abundance at the L-T transition and to the methane depletion observed on young exoplanets. Finally, we predict the existence of a continuum of brown dwarfs and exoplanets for absolute J magnitude=15-18 and J-K color=0-3, due to the evolution of the L-T transition with gravity. This self-consistent model therefore provides a general framework to understand the effects of clouds and appears well-suited for atmospheric retrievals.

Whole planet coupling between climate, mantle, and core: Implications for the evolution of rocky planets

Whole planet coupling between climate, mantle, and core: Implications for the evolution of rocky planets

Authors:

Foley et al

Abstract:

Earth's climate, mantle, and core interact over geologic timescales. Climate influences whether plate tectonics can take place on a planet, with cool climates being favorable for plate tectonics because they enhance stresses in the lithosphere, suppress plate boundary annealing, and promote hydration and weakening of the lithosphere. Plate tectonics plays a vital role in the long-term carbon cycle, which helps to maintain a temperate climate. Plate tectonics provides long-term cooling of the core, which is vital for generating a magnetic field, and the magnetic field is capable of shielding atmospheric volatiles from the solar wind. Coupling between climate, mantle, and core can potentially explain the divergent evolution of Earth and Venus. As Venus lies too close to the sun for liquid water to exist, there is no long-term carbon cycle and thus an extremely hot climate. Therefore plate tectonics cannot operate and a long-lived core dynamo cannot be sustained due to insufficient core cooling. On planets within the habitable zone where liquid water is possible, a wide range of evolutionary scenarios can take place depending on initial atmospheric composition, bulk volatile content, or the timing of when plate tectonics initiates, among other factors. Many of these evolutionary trajectories would render the planet uninhabitable. However, there is still significant uncertainty over the nature of the coupling between climate, mantle, and core. Future work is needed to constrain potential evolutionary scenarios and the likelihood of an Earth-like evolution.

NIR-driven Moist Upper Atmospheres of Synchronously Rotating Temperate Terrestrial Exoplanets

NIR-driven Moist Upper Atmospheres of Synchronously Rotating Temperate Terrestrial Exoplanets 

Authors:

Fuji et al

Abstract:
H2O is a key molecule in characterizing atmospheres of temperate terrestrial planets, and observations of transmission spectra are expected to play a primary role in detecting its signatures in the near future. The detectability of H2O absorption features in transmission spectra depends on the abundance of water vapor in the upper part of the atmosphere. We study the three-dimensional distribution of atmospheric H2O for synchronously rotating Earth-sized aquaplanets using the general circulation model (GCM) ROCKE-3D, and examine the effects of total incident flux and stellar spectral type. We observe a more gentle increase of the water vapor mixing ratio in response to increased incident flux than one-dimensional models suggest, in qualitative agreement with the climate-stabilizing effect of clouds around the substellar point previously observed in GCMs applied to synchronously rotating planets. However, the water vapor mixing ratio in the upper atmosphere starts to increase while the surface temperature is still moderate. This is explained by the circulation in the upper atmosphere being driven by the radiative heating due to absorption by water vapor and cloud particles, causing efficient vertical transport of water vapor. Consistently, the water vapor mixing ratio is found to be well-correlated with the near-infrared portion of the incident flux. We also simulate transmission spectra based on the GCM outputs, and show that for the more highly irradiated planets, the H2O signatures may be strengthened by a factor of a few, loosening the observational demands for a H2O detection.

GCM Simulations of Unstable Climates in the Habitable Zone

GCM Simulations of Unstable Climates in the Habitable Zone

Authors:

Paradise et al

Abstract:

It has recently been proposed that Earth-like planets in the outer regions of the habitable zone experience unstable climates, repeatedly cycling between glaciated and deglaciated climatic states. While this result has been confirmed and also extended to explain early Mars climate records, all existing work relies on highly idealized low-dimensional climate models. Here, we confirm that the phenomenology of climate cycles remains in 3D Earth climate models with considerably more degrees of freedom. To circumvent the computational barrier of integrating climate on Gyr timescales, we implement a hybrid 0D-3D integrator that uses a general circulation model (GCM) as a short relaxation step along a long evolutionary climate sequence. We find that GCM climate cycles are qualitatively consistent with reported low-dimensional results. This establishes on a firmer ground the notion that outer habitable zone planets may be preferentially found in transiently glaciated states.

The Thermal Phase Curve Offset on Tidally and Nontidally Locked Exoplanets

The Thermal Phase Curve Offset on Tidally and Nontidally Locked Exoplanets: A Shallow Water Model

Authors:

Penn et al

Abstract:
Using a shallow water model with time-dependent forcing, we show that the peak of an exoplanet thermal phase curve is, in general, offset from the secondary eclipse when the planet is rotating. That is, the planetary hot spot is offset from the point of maximal heating (the substellar point) and may lead or lag the forcing; the extent and sign of the offset are functions of both the rotation rate and orbital period of the planet. We also find that the system reaches a steady state in the reference frame of the moving forcing. The model is an extension of the well-studied Matsuno–Gill model into a full spherical geometry and with a planetary-scale translating forcing representing the insolation received on an exoplanet from a host star. The speed of the gravity waves in the model is shown to be a key metric in evaluating the phase curve offset. If the velocity of the substellar point (relative to the planet's surface) exceeds that of the gravity waves, then the hot spot will lag the substellar point, as might be expected by consideration of forced gravity wave dynamics. However, when the substellar point is moving slower than the internal wave speed of the system, the hottest point may lead the passage of the forcing. We provide an interpretation of this result by consideration of the Rossby and Kelvin wave dynamics, as well as, in the very slowly rotating case, a one-dimensional model that yields an analytic solution. Finally, we consider the inverse problem of constraining planetary rotation rate from an observed phase curve.

Radiation Hydrodynamics Simulations of Photoevaporation of Protoplanetary Disks: Metallicity Dependence

Radiation Hydrodynamics Simulations of Photoevaporation of Protoplanetary Disks: Metallicity Dependence

Authors:
Nakatani et al

Abstract:

Protoplanetary disks are thought to have lifetimes of 3−6 million years in the solar neighborhood, but recent observations suggest that the disk lifetimes are shorter in a low metallicity environment. We perform a suite of radiation hydrodynamics simulations of photoevaporation of protoplanetary disks to study the disk structure and its long-term evolution of ∼10000 years, and the metallicity dependence of mass-loss rate. Our simulations follow hydrodynamics, extreme and far ultraviolet radiative transfer, and non-equilibrium chemistry in a self-consistent manner. Dust grain temperatures are also calculated consistently by solving the radiative transfer of the stellar irradiation and grain (re-)emission. We vary the disk gas metallicity over a wide range of 10−4Z⊙≤Z≤10 Z⊙. For our fiducial model with a 0.5 M⊙ central star with solar metallicity, the time-averaged photoevaporation rate is M˙ph=1.38×10−8M⊙yr−1. The photoevaporation rate is lower with higher metallicity in the range of 10−0.5Z⊙≲Z≲10 Z⊙, because dust shielding effectively prevents far-ultra violet (FUV) photons from penetrating into and heating the dense regions of the disk. The photoevaporation rate sharply declines at even lower metallicities in 10−1Z⊙≲Z≲10−0.5Z⊙, because FUV photoelectric heating is not efficient any more to raise the gas temperature and to drive outflows. At 10−4Z⊙≤Z≲10−1Z⊙, HI photoionization heating acts as a dominant gas heating process and drives photoevaporative flows with roughly a constant rate. The typical disk lifetime is shorter at Z=0.3 Z⊙ than at Z=Z⊙, being consistent with recent observations of the extreme outer galaxy.

The Planetary Accretion Shock: I. Framework for Radiation-hydrodynamical Simulations and First Results

The Planetary Accretion Shock: I. Framework for Radiation-hydrodynamical Simulations and First Results

Authors:

Marleau et al

Abstract:
The key aspect determining the post-formation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the planetary surface. We use 1D radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain post-shock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux which escapes the system (roughly the Hill sphere), taking into account the energy recycling which occurs ahead of the shock in a radiative precursor. We focus here on a constant equation of state to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will require a self-consistent treatment including hydrogen dissocation and ionization, the results show the correct qualitative behavior and can be understood semi-analytically. The shock is found to be isothermal and supercritical for a range of conditions relevant to core accretion (CA), with Mach numbers greater than ca. 3. Across the shock, the entropy decreases significantly, by a few entropy units (k_B/baryon). While nearly 100 percent of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40 percent, implying that a meaningful fraction of the total accretion energy is brought into the planet. For realistic parameter combinations in the CA scenario, a non-zero fraction of the luminosity always escapes the system. This luminosity could explain, at least in part, recent observations in the LkCa 15 and HD 100546 systems.

Collisional Fragmentation is not a Barrier to Close-in Planet Formation

Collisional Fragmentation is not a Barrier to Close-in Planet Formation 

Authors:

Wallace et al

Abstract:
Collisional fragmentation is shown to not be a barrier to rocky planet formation at small distances from the host star. Simple analytic arguments demonstrate that rocky planet formation via collisions of homogeneous gravity-dominated bodies is possible down to distances of order the Roche radius (rRoche). Extensive N-body simulations that include plausible models for fragmentation and merging of gravity-dominated bodies confirm this conclusion and demonstrate that rocky planet formation is possible down to ∼1.1 rRoche. At smaller distances, tidal effects cause collisions to be too fragmenting to allow mass build-up to a final, dynamically stable planetary system. We argue that even differentiated bodies can accumulate to form planets at distances that are not much larger than rRoche.

N-body simulations of planet formation via pebble accretion I: First Results

N-body simulations of planet formation via pebble accretion I: First Results

Authors:

Matsumura et al

Abstract:
Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA, along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on a disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle because, although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a range of densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets.

Simulations Fail to Explain SuperEarth Resonant Chains

Breaking the Chains: Hot Super-Earth systems from migration and disruption of compact resonant chains
Authors:

Izidoro et al

Abstract:
“Hot super-Earths” (or “Mini-Neptunes”) between 1 and 4 times Earth’s size with period shorter than 100 days orbit 30-50% of Sun-like type stars. Their orbital configuration – measured as the period ratio distribution of adjacent planets in multi-planet systems – is a strong constraint for formation models. Here we use N-body simulations with synthetic forces from an underlying evolving gaseous disk to model the formation and long-term dynamical evolution of super-Earth systems. While the gas disk is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disk. These resonant chains are far more compact than the observed super-Earth systems. Once the gas dissipates resonant chains may become dynamically unstable. They undergo a phase of giant impacts that spreads the systems out. Disk turbulence has no measurable effect on the outcome. Our simulations match observations if a small fraction of resonant chains remain stable, while most super-Earths undergo a late dynamical instability. Our statistical analysis restricts the contribution of stable systems to less than 25%. Our results also suggest that the large fraction of observed single planet systems does not necessarily imply any dichotomy in the architecture of planetary systems. Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in ∼5% of the cases. This leads to a mystery: in our simulations only 50-60% of resonant chains became unstable whereas at least 75% (and probably 90-95%) must be unstable to match observations.

Zoom-Simulations of Protoplanetary Disks starting from GMC scales

Zoom-Simulations of Protoplanetary Disks starting from GMC scales

Authors:

Kuffmeier et al

Abstract:

We investigate the formation of protoplanetary disks around nine solar mass stars formed in the context of a 40 pc Giant Molecular Cloud model, using \ramses \ adaptive-mesh resolution simulations extending over a scale range of about 4 million, from an outer scale of 40 pc to a smallest cell size of 2 AU. Our most important qualitative result is that the accretion process is heterogeneous in multiple ways; in time, in space, and among protostars of otherwise similar mass. Accretion is heterogeneous in time, in the sense that accretion rates vary in time, with generally decreasing profiles, whose slopes can vary over a wide range, and where accretion can increase again if a protostar enters a region with increased density and low speed. Accretion is heterogeneous in space, because of the distribution of mass in space, with mass approaching the accreting star and disk in filaments and sheets. Finally, accretion is heterogeneous among stars, since the detailed conditions and dynamics in the neighborhood of each star can vary widely. We investigate in particular the sensitivity of disk formation to physical conditions, such as mass-to-flux ratio, and to sink particle recipe parameters. We find that disks frequently form, even when choosing the least favorable sink particle parameters, and that turbulence carried along from larger scales is a main factor in allowing disks to form even when the magnetic field is comparatively strong.

Effects of variable eccentricity on the climate of an Earth-like world

Effects of variable eccentricity on the climate of an Earth-like world

Authors:

Way et al

Abstract:

The Kepler era of exoplanetary discovery has presented the Astronomical community with a cornucopia of planetary systems very different from the one which we inhabit. It has long been known that Jupiter plays a major role in the orbital parameters of Mars and it's climate, but there is also a long-standing belief that Jupiter would play a similar role for Earth if not for its large moon. Using a three dimensional general circulation model (3-D GCM) with a fully-coupled ocean we simulate what would happen to the climate of an Earth-like world if Mars did not exist, but a Jupiter-like planet was much closer to Earth's orbit. We investigate two scenarios that involve evolution of the Earth-like planet's orbital eccentricity from 0--0.283 over 6500 years, and from 0--0.066 on a time scale of 4500 years. In both cases we discover that they would maintain relatively temperate climates over the time-scales simulated and that their regional habitability is larger than present day Earth. More Earth-like planets in multi-planet systems will be discovered as we continue to survey the skies and the results herein show that the proximity of large gas giant planets may play an important role in the habitability of these worlds. These are the first such 3-D GCM simulations using a fully-coupled ocean with a planetary orbit that evolves over time due to the presence of a giant planet.

Linking long-term planetary N-body simulations with periodic orbits: application to white dwarf pollution

Linking long-term planetary N-body simulations with periodic orbits: application to white dwarf pollution

Authors:

Antoniadou et al

Abstract:

Mounting discoveries of debris discs orbiting newly-formed stars and white dwarfs (WDs) showcase the importance of modeling the long-term evolution of small bodies in exosystems. WD debris discs are in particular thought to form from very long-term (0.1-5.0 Gyr) instability between planets and asteroids. However, the time-consuming nature of N-body integrators which accurately simulate motion over Gyrs necessitates a judicious choice of initial conditions. The analytical tools known as \textit{periodic orbits} can circumvent the guesswork. Here, we begin a comprehensive analysis directly linking periodic orbits with N-body integration outcomes with an extensive exploration of the planar circular restricted three-body problem (CRTBP) with an outer planet and inner asteroid near or inside of the 2:1 mean motion resonance. We run nearly 1000 focused simulations for the entire age of the Universe (14 Gyr) with initial conditions mapped to the phase space locations surrounding the unstable and stable periodic orbits for that commensurability. In none of our simulations did the planar CRTBP architecture yield a long-timescale (≳0.25% of the age of the Universe) asteroid-star collision. The pericentre distance of asteroids which survived beyond this timescale (≈35 Myr) varied by at most about 60%. These results help affirm that collisions occur too quickly to explain WD pollution in the planar CRTBP 2:1 regime, and highlight the need for further periodic orbit studies with the eccentric and inclined TBP architectures and other significant orbital period commensurabilities.

Convection in Condensible-rich Atmospheres

Convection in Condensible-rich Atmospheres

Authors:

Ding et al

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

There is Considerable Uncertainty About the Inner Edge of the Habitable Zone

DIFFERENCES IN WATER VAPOR RADIATIVE TRANSFER AMONG 1D MODELS CAN SIGNIFICANTLY AFFECT THE INNER EDGE OF THE HABITABLE ZONE

Authors:

Yang et al

Abstract:

An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G and M stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 μm) and in the region between 0.2 and 1.5 μm. Differences in outgoing longwave radiation increase with surface temperature and reach 10–20 W m−2; differences in shortwave reach up to 60 W m−2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ≈10% of modern Earth's solar constant (i.e., ≈34 W m−2 in global mean) among band models and ≈3% between the two line-by-line models. These comparisons show that future work is needed that focuses on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

Examining the Accuracy of Astrophysical Disk Simulations With a Generalized Hydrodynamical Test Problem

Examining the Accuracy of Astrophysical Disk Simulations With a Generalized Hydrodynamical Test Problem

Authors:

Raskin et al

Abstract:

We discuss a generalization of the classic Keplerian disk test problem allowing for both pressure and rotational support, as a method of testing astrophysical codes incorporating both gravitation and hydrodynamics. We argue for the inclusion of pressure in rotating disk simulations on the grounds that realistic, astrophysical disks exhibit non-negligible pressure support. We then apply this test problem to examine the performance of various smoothed particle hydrodynamics (SPH) methods incorporating a number of improvements proposed over the years to help SPH better address problems noted in modeling the classical gravitation only Keplerian disk. We also apply this test to a newly developed extension of SPH based on reproducing kernels called CRKSPH. Counterintuitively, we find that pressure support worsens the performance of traditional SPH on this problem, causing unphysical collapse away from the steady-state disk solution even more rapidly than the purely gravitational problem, whereas CRKSPH greatly reduces this error.

LONG-TERM AND LARGE-SCALE HYDRODYNAMICAL SIMULATIONS OF MIGRATING PLANETS

LONG-TERM AND LARGE-SCALE HYDRODYNAMICAL SIMULATIONS OF MIGRATING PLANETS

Authors:

Benítez-Llambay et al

Abstract:

We present a new method that allows for long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists of a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiency problems related to the migration of multi-planet systems in gaseous disks and to improve the current results of the migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.

Particle dynamics in discs with turbulence generated by the vertical shear instability

Particle dynamics in discs with turbulence generated by the vertical shear instability

Authors:

Stoll et al

Abstract:

Among the candidates for generating turbulence in accretion discs in situations with low intrinsic ionization the vertical shear instability (VSI) has become an interesting candidate, as it relies purely on a vertical gradient in the angular velocity. Existing simulations have shown that α-values a few times 10−4 can be generated. The particle growth in the early planet formation phase is determined by the dynamics of dust particles. Here, we address in particular the efficiency of VSI-turbulence in concentrating particles in order to generate overdensities and low collision velocities. We perform 3D numerical hydrodynamical simulations of accretion discs around young stars that include radiative transport and irradiation from the central star. The motion of particles within a size range of a fraction of mm up to several m is followed using standard drag formula. We confirm that under realistic conditions the VSI is able to generate turbulence in full 3D protoplanetary discs. The irradiated disc shows turbulence within 10 to 60au. The mean radial motion of the gas is such that it is directed inward near the midplane and outward in the surface layers. We find that large particles drift inward with the expected speed, while small particles can experience phases of outward drift. Additionally, the particles show bunching behaviour with overdensities reaching 5 times the average value, which is strongest for dimensionless stopping times around unity. Particles in a VSI-turbulent discs are concentrated in large scale turbulent eddies and show low relative speeds that allow for growing collisions. The reached overdensities will also allow for the onset streaming instabilities further enhancing particle growth. The outward drift for small particles at higher disk elevations allows for the transport of processed high temperature material in the Solar System to larger distances.

Temperature Structure and Atmospheric Circulation of Dry, Tidally Locked Rocky Exoplanets

Temperature Structure and Atmospheric Circulation of Dry, Tidally Locked Rocky Exoplanets

Authors:

Koll et al

Abstract:

Next-generation space telescopes will observe the atmospheres of rocky planets orbiting nearby M-dwarfs. Understanding these observations will require well-developed theory in addition to numerical simulations. Here we present theoretical models for the temperature structure and atmospheric circulation of dry, tidally locked rocky exoplanets with grey radiative transfer and test them using a general circulation model (GCM). First, we develop a radiative-convective model that captures surface temperatures of slowly rotating and cool atmospheres. Second, we show that the atmospheric circulation acts as a global heat engine, which places strong constraints on large-scale wind speeds. Third, we develop a radiative-convective-subsiding model which extends our radiative-convective model to hot and thin atmospheres. We find that rocky planets develop large day-night temperature gradients at a ratio of wave-to-radiative timescales up to two orders of magnitude smaller than the value suggested by work on hot Jupiters. The small ratio is due to the heat engine inefficiency and asymmetry between updrafts and subsidence in convecting atmospheres. Fourth, we show using GCM simulations that rotation only has a strong effect on temperature structure if the atmosphere is hot or thin. Our models let us map out atmospheric scenarios for planets such as GJ 1132b and show how thermal phase curves could constrain them. Measuring phase curves of short-period planets will require similar amounts of time on the James Webb Space Telescope as detecting molecules via transit spectroscopy, so future observations should pursue both techniques.

Assessing the Habitability of Exoplanetary Earth-like Atmospheres Through Simualtions

Assessing the habitability of planets with Earth-like atmospheres with 1D and 3D climate modeling

Authors:

Godolt et al

Abstract:

The habitable zone (HZ) describes the range of orbital distances around a star where the existence of liquid water on the surface of an Earth-like planet is in principle possible. While 3D climate studies can calculate the water vapor, ice albedo, and cloud feedback self-consistently and therefore allow for a deeper understanding and the identification of relevant climate processes, 1D model studies rely on fewer model assumptions and can be more easily applied to the large parameter space possible for exoplanets. We evaluate the applicability of 1D climate models to estimate the potential habitability of Earth-like exoplanets by comparing our 1D model results to those of 3D climate studies in the literature. We applied a cloud-free 1D radiative-convective climate model to calculate the climate of Earth-like planets around different types of main-sequence stars with varying surface albedo and relative humidity profile. These parameters depend on climate feedbacks that are not treated self-consistently in most 1D models. We compared the results to those of 3D model calculations in the literature and investigated to what extent the 1D model can approximate the surface temperatures calculated by the 3D models. The 1D parameter study results in a large range of climates possible for an Earth-sized planet with an Earth-like atmosphere and water reservoir at a certain stellar insolation. At some stellar insolations the full spectrum of climate states could be realized, i.e., uninhabitable conditions as well as habitable surface conditions, depending only on the relative humidity and surface albedo assumed. When treating the surface albedo and the relative humidity profile as parameters in 1D model studies and using the habitability constraints found by recent 3D modeling studies, the same conclusions about the potential habitability of a planet can be drawn as from 3D model calculations.