Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an inter-comparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (20-30 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance.

The James Webb Space Telescope (JWST) will open up the possibility of comprehensively measuring the emission spectra of rocky exoplanets orbiting M dwarfs to detect and characterize their atmospheres. In preparation for this opportunity, we present model atmospheres for three M-dwarf planets particularly amenable to secondary eclipse spectroscopy -- TRAPPIST-1b, GJ 1132b, and LHS 3844b. Using three limiting cases of candidate atmospheric compositions (pure H2O, pure CO2 and solar abundances) we calculate temperature-pressure profiles and emission and reflection spectra in radiative-convective equilibrium, including the effects of a solid surface. We find that the atmospheric radiative transfer is significantly influenced by the cool M-star irradiation; H2O and CO2 absorption bands in the near-infrared are strong enough to absorb a sizeable fraction of the incoming stellar light at low pressures, which leads to temperature inversions in the upper atmosphere. The non-gray band structure of gaseous opacities in the infrared is hereby an important factor. Opacity windows are muted at higher atmospheric temperatures, so we expect temperature inversions to be common only for sufficiently cool planets. We also find that pure CO2 atmospheres exhibit lower overall temperatures and stronger reflection spectra compared to models of the other two compositions. We estimate that for GJ 1132b and LHS 3844b we should be able to distinguish between different atmospheric compositions with JWST. The emission lines from the predicted temperature inversions are currently hard to measure, but high resolution spectroscopy with future Extremely Large Telescopes may be able to detect them.

The upcoming launch of the James Webb Space Telescope (JWST) means that we will soon have the capability to characterize the atmospheres of rocky exoplanets. However, it is still unknown whether such planets orbiting close to M dwarf stars can retain their atmospheres, or whether high-energy irradiation from the star will strip the gaseous envelopes from these objects. We present a new method to detect an atmosphere on a synchronously rotating rocky exoplanet around a K/M dwarf, by using thermal emission during secondary eclipse to infer a high dayside albedo that could only be explained by bright clouds. Based on calculations for plausible surface conditions, we conclude that a high albedo could be unambiguously interpreted as a signal of an atmosphere for planets with substellar temperatures of $T_{sub}=$ 410-1250 K. This range corresponds to equilibrium temperatures of $T_{eq}=$ 300-880 K. We compare the inferred albedos of eight possible planet surface compositions to cloud albedo calculations. We determine that a layer of clouds with optical depths greater than $\tau=0.5$-$7$, would have high enough albedos to be distinguishable from a bare rock surface. This method of detecting an atmosphere on a rocky planet is complementary to existing methods for detecting atmospheres, because it provides a way to detect atmospheres with pressures below 1 bar (e.g. Mars), which are too tenuous to transport significant heat but thick enough to host high-albedo clouds.

We study the red-optical photometry of the ultra-hot Jupiter WASP-121 b as observed by the Transiting Exoplanet Survey Satellite (TESS) and model its atmosphere through a radiative transfer simulation. Given its short orbital period of 1.275 days, inflated state and bright host star, WASP-121 b is exceptionally favorable for detailed atmospheric characterization. Towards this purpose, we use allesfitter to characterize its full red-optical phase curve, including the planetary phase modulation and the secondary eclipse. We measure the day and nightside brightness temperatures in the TESS passband as 2940+38-41 K and 2190+294-146 K, respectively, and find no phase shift between the brightest and substellar points. This is consistent with an inefficient heat recirculation on the planet. We then perform an atmospheric retrieval analysis to infer the dayside atmospheric properties of WASP-121 b such as its bulk composition, albedo and heat recirculation. We confirm the temperature inversion in the atmosphere and suggest H-, TiO and VO as potential causes of the inversion, absorbing heat at optical wavelengths at low pressures. Future HST and JWST observations of WASP-121 b will benefit from its first full phase curve measured by TESS.

Most known terrestrial planets orbit small stars with radii less than 60% that of the Sun. Theoretical models predict that these planets are more vulnerable to atmospheric loss than their counterparts orbiting Sun-like stars. To determine whether a thick atmosphere has survived on a small planet, one approach is to search for signatures of atmospheric heat redistribution in its thermal phase curve. Previous phase curve observations of the super-Earth 55 Cancri e (1.9 Earth radii) showed that its peak brightness is offset from the substellar point $-$ possibly indicative of atmospheric circulation. Here we report a phase curve measurement for the smaller, cooler planet LHS 3844b, a 1.3 Earth radius world in an 11-hour orbit around a small, nearby star. The observed phase variation is symmetric and has a large amplitude, implying a dayside brightness temperature of $1040\pm40$ kelvin and a nightside temperature consistent with zero kelvin (at one standard deviation). Thick atmospheres with surface pressures above 10 bar are ruled out by the data (at three standard deviations), and less-massive atmospheres are unstable to erosion by stellar wind. The data are well fitted by a bare rock model with a low Bond albedo (lower than 0.2 at two standard deviations). These results support theoretical predictions that hot terrestrial planets orbiting small stars may not retain substantial atmospheres.

Young terrestrial planets can capture or outgas hydrogen-rich atmospheres with tens to hundreds of bars of H2, which persist for 100 Myrs or longer. Although the earliest habitable conditions on Earth and terrestrial exoplanets could thus arise while the atmosphere is still dominated by H2, the climatic effects of H2 remain poorly understood. Previous work showed that H2 induces strong greenhouse warming at the outer edge of the habitable zone. Here we use a 1D radiative-convective model to show that H2 also leads to strong warming near the inner edge of the habitable zone. Unlike H2's greenhouse warming at the outer edge, however, its effect near the inner edge is driven by thermodynamics: H2's large thermal scale height allows the atmosphere to store more water vapor than either a pure-H2O atmosphere or an atmosphere with a heavy background gas, such as N2 or CO2, thereby amplifying the greenhouse effect of H2O. Using idealized grey calculations, we then present a general argument for how different background gases affect the inner edge of the habitable zone. H2 stands out for its ability to induce novel "Souffle" climates, which further support its warming effect. Our results show that if the earliest conditions on a planet near the inner edge of the habitable zone were H2-rich, they were likely also hot: 1 bar of H2 is sufficient to raise surface temperatures above 340 K, and 50 bar of H2 are sufficient to raise surface temperatures above 450 K.

Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of secondary eclipses could quickly identify "candidate" atmospheres, by searching for rocky planets with atmospheres thick enough that atmospheric heat transport noticeably reduces their dayside thermal emission compared to that of a bare rock. For a planet amenable to atmospheric follow-up, we find that JWST should be able to confidently detect the heat redistribution signal of an O(1) bar atmosphere with one to two eclipses. One to two eclipses is generally much less than the effort needed to infer an atmosphere via transmission or emission spectroscopy. Candidate atmospheres can be further validated via follow-up spectroscopy or phase curves. In addition, because this technique is fast it could enable a first atmospheric survey of rocky exoplanets with JWST. We estimate that the TESS mission will find ~100 planets that are too hot to be habitable but that can be quickly probed via eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have atmospheres is important not only for understanding the evolution of uninhabitable worlds: if atmospheres are common on hot planets, then cooler, potentially habitable planets around M dwarfs are also likely to have atmospheres.

Atmospheric heat redistribution shapes the remote appearance of rocky exoplanets but there is currently no easy way to predict a planet's heat redistribution from its physical properties. In this paper I derive an analytical scaling theory for the heat redistribution on tidally locked rocky exoplanets. The main parameters of the theory are a planet's equilibrium temperature, its surface pressure, and its broadband longwave optical thickness. I validate the theory against general circulation model simulations of TRAPPIST-1b, GJ1132b, and LHS 3844b. I find that heat redistribution becomes efficient, and a planet's observable thermal phase curve and secondary eclipse start to deviate significantly from that of a bare rock, once surface pressure exceeds O(1) bar. These results thus bridge the gap between theory and imminent observations with the James Webb Space Telescope. They can also be used to parameterize the effect of 3D atmospheric dynamics in 1D models, thereby improving the self-consistency of such models.

Recent studies have shown that ocean dynamics can have a significant warming effect on the permanent night sides of 1 to 1 tidally locked terrestrial exoplanets with Earth-like atmospheres and oceans in the middle of the habitable zone. However, the impact of ocean dynamics on the habitable zone's boundaries (inner edge and outer edge) is still unknown and represents a major gap in our understanding of this type of planets. Here we use a coupled atmosphere-ocean global climate model to show that planetary heat transport from the day to night side is dominated by the ocean at lower stellar fluxes and by the atmosphere near the inner edge of the habitable zone. This decrease in oceanic heat transport (OHT) at high stellar fluxes is mainly due to weakening of surface wind stress and a decrease in surface shortwave energy deposition. We further show that ocean dynamics have almost no effect on the observational thermal phase curves of planets near the inner edge of the habitable zone. For planets in the habitable zone's middle range, ocean dynamics moves the hottest spot on the surface eastward from the substellar point. These results suggest that future studies of the inner edge may devote computational resources to atmosphere-only processes such as clouds and radiation. For studies of the middle range and outer edge of the habitable zone, however, fully coupled ocean-atmosphere modeling will be necessary. Note that due to computational resource limitations, only one rotation period (60 Earth days) has been systematically examined in this study; future work varying rotation period as well as other parameters such as atmospheric mass and composition is required.

A key legacy of the recently launched TESS mission will be to provide the astronomical community with many of the best transiting exoplanet targets for atmospheric characterization. However, time is of the essence to take full advantage of this opportunity. JWST, although delayed, will still complete its nominal five year mission on a timeline that motivates rapid identification, confirmation, and mass measurement of the top atmospheric characterization targets from TESS. Beyond JWST, future dedicated missions for atmospheric studies such as ARIEL require the discovery and confirmation of several hundred additional sub-Jovian size planets (R_p < 10 R_Earth) orbiting bright stars, beyond those known today, to ensure a successful statistical census of exoplanet atmospheres. Ground-based ELTs will also contribute to surveying the atmospheres of the transiting planets discovered by TESS. Here we present a set of two straightforward analytic metrics, quantifying the expected signal-to-noise in transmission and thermal emission spectroscopy for a given planet, that will allow the top atmospheric characterization targets to be readily identified among the TESS planet candidates. Targets that meet our proposed threshold values for these metrics would be encouraged for rapid follow-up and confirmation via radial velocity mass measurements. Based on the catalog of simulated TESS detections by Sullivan et al. (2015), we determine appropriate cutoff values of the metrics, such that the TESS mission will ultimately yield a sample of $\sim300$ high-quality atmospheric characterization targets across a range of planet size bins, extending down to Earth-size, potentially habitable worlds.

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 (1D) 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 um) and in the region between 0.2 and 1.5 um. Differences in outgoing longwave radiation increase with surface temperature and reach 10-20 Wm^-2; differences in shortwave reach up to 60 Wm^-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 Wm^-2 in global mean) among band models and ~3% between the two line-by-line models. These comparisons show that future work is needed focusing on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

Because of their intense incident stellar irradiation and likely tidally locked spin states, hot Jupiters are expected to have wind speeds that approach or exceed the speed of sound. In this work we develop a theory to explain the magnitude of these winds. We model hot Jupiters as planetary heat engines and show that hot Jupiters are always less efficient than an ideal Carnot engine. Next, we demonstrate that our predicted wind speeds match those from three-dimensional numerical simulations over a broad range of parameters. Finally, we use our theory to evaluate how well different drag mechanisms can match the wind speeds observed with Doppler spectroscopy for HD 189733b and HD 209458b. We find that magnetic drag is potentially too weak to match the observations for HD 189733b, but is compatible with the observations for HD 209458b. In contrast, shear instabilities and/or shocks are compatible with both observations. Furthermore, the two mechanisms predict different wind speed trends for hotter and colder planets than currently observed. As a result, we propose that a wider range of Doppler observations could reveal multiple drag mechanisms at play across different hot Jupiters.

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.

Next-generation space telescopes will allow us to characterize terrestrial exoplanets. To do so effectively it will be crucial to make use of all available data. We investigate which atmospheric properties can, and cannot, be inferred from the broadband thermal phase curve of a dry and tidally locked terrestrial planet. First, we use dimensional analysis to show that phase curves are controlled by six nondimensional parameters. Second, we use an idealized general circulation model (GCM) to explore the relative sensitivity of phase curves to these parameters. We find that the feature of phase curves most sensitive to atmospheric parameters is the peak-to-trough amplitude. Moreover, except for hot and rapidly rotating planets, the phase amplitude is primarily sensitive to only two nondimensional parameters: 1) the ratio of dynamical to radiative timescales, and 2) the longwave optical depth at the surface. As an application of this technique, we show how phase curve measurements can be combined with transit or emission spectroscopy to yield a new constraint for the surface pressure and atmospheric mass of terrestrial planets. We estimate that a single broadband phase curve, measured over half an orbit with the James Webb Space Telescope, could meaningfully constrain the atmospheric mass of a nearby super-Earth. Such constraints will be important for studying the atmospheric evolution of terrestrial exoplanets as well as characterizing the surface conditions on potentially habitable planets.