Hurricanes are one of the most extreme storm systems that occur on Earth, characterized by strong rainfall and fast winds. The terrestrial exoplanets that will be characterized with future infrared space telescopes orbit M dwarf stars. As a result, the best observable terrestrial exoplanets have vastly different climates than Earth, with a large dayside-to-nightside irradiation contrast and relatively slow rotation. Hurricanes may affect future observations of terrestrial exoplanets because they enhance the vertical transport of water vapor and could influence ocean heat transport. In this work, we explore how the environment of terrestrial exoplanets orbiting M dwarf stars affects the favorability of hurricane genesis (formation). To do so, we apply metrics developed to understand hurricane genesis on Earth to three-dimensional climate models of ocean-covered exoplanets orbiting M dwarf stars. We find that hurricane genesis is most favorable on intermediate-rotating tidally locked terrestrial exoplanets with rotation periods of $\sim 8-10~\mathrm{days}$. As a result, hurricane genesis is most favorable for terrestrial exoplanets in the habitable zones of late-type M dwarf stars. The peak in the favorability of hurricane genesis at intermediate rotation occurs because sufficient spin is required for hurricane genesis, but the vertical wind shear on fast-rotating terrestrial exoplanets disrupts hurricane genesis. We find that hurricane genesis is less favorable on slowly rotating terrestrial exoplanets, which agrees with previous work. Future work using simulations that resolve hurricane genesis and evolution can test our expectations for how the environment affects the favorability of hurricane genesis on tidally locked terrestrial exoplanets.

The potential habitability of tidally locked planets orbiting M-dwarf stars has been widely investigated in recent work, typically with a non-dynamic ocean and without continents. On Earth, ocean dynamics are a primary means of heat and nutrient distribution. Continents are a critical source of nutrients, strongly influence ocean dynamics, and participate in climate regulation. In this work, we investigate how the size of a substellar land mass affects the oceans ability to transport heat and upwell nutrients on the tidally locked planet Proxima Centauri b using the ROCKE-3D coupled ocean-atmosphere General Circulation Model (GCM). We find that dayside ice-free ocean and nutrient delivery to the mixed layer via upwelling are maintained across all continent sizes. We also find that Proxima Centauri bs climate is more sensitive to differences among atmospheric GCMs than to the inclusion of ocean dynamics in ROCKE-3D. Finally, we find that Proxima Centauri b transitions from a lobster state where ocean heat transport distributes heat away from the substellar point to an eyeball state where heat transport is restricted and surface temperature decreases symmetrically from the substellar point when the continent size exceeds about 20 percent of the surface area. Our work suggests that both a dynamic ocean and continents are unlikely to decrease the habitability prospects of nearby tidally locked targets like Proxima Centauri b that could be investigated with future observations by the James Webb Space Telescope (JWST).

Liquid water oceans are at the center of our search for life on exoplanets because water is a strict requirement for life as we know it. However, oceans are dynamic habitats---and some oceans may be better hosts for life than others. In Earth's ocean, circulation transports essential nutrients such as phosphate and is a first-order control on the distribution and productivity of life. Of particular importance is upward flow from the dark depths of the ocean in response to wind-driven divergence in surface layers. This `upwelling' returns essential nutrients that tend to accumulate at depth via sinking of organic particulates back to the sunlit regions where photosynthetic life thrives. Ocean dynamics are likely to impose constraints on the activity and atmospheric expression of photosynthetic life in exo-oceans as well, but we lack an understanding of how ocean dynamics may differ on other planets. We address this issue by exploring the sensitivity of ocean dynamics to a suite of planetary parameters using ROCKE-3D, a fully coupled ocean-atmosphere GCM. Our results suggest that planets that rotate slower and have higher surface pressure than Earth may be the most attractive targets for remote life detection because upwelling is enhanced under these conditions, resulting in greater nutrient supply to the surface biosphere. Seasonal deepening of the mixed layer on high obliquity planets may also enhance nutrient replenishment from depth into the surface mixed layer. Efficient nutrient recycling favors greater biological activity, more biosignature production, and thus more detectable life. More generally, our results demonstrate the importance of considering oceanographic phenomena for exoplanet life detection and motivate future interdisciplinary contributions to the emerging field of exo-oceanography.

The Gaia hypothesis postulates that life regulates its environment to be favorable for its own survival. Most planets experience numerous perturbations throughout their lifetimes such as asteroid impacts, volcanism, and the evolution of their host star's luminosity. For the Gaia hypothesis to be viable, life must be able to keep the conditions of its host planet habitable, even in the face of these challenges. ExoGaia, a model created to investigate the Gaia hypothesis, has been previously used to demonstrate that a randomly mutating biosphere is in some cases capable of maintaining planetary habitability. However, those model scenarios assumed that all non-biological planetary parameters were static, neglecting the inevitable perturbations that real planets would experience. To see how life responds to climate perturbations to its host planet, we created three climate perturbations in ExoGaia: one rapid cooling of a planet and two heating events, one rapid and one gradual. The planets on which Gaian feedbacks emerge without climate perturbations are the same planets on which life is most likely to survive each of our perturbation scenarios. Biospheres experiencing gradual changes to the environment are able to survive changes of larger magnitude than those experiencing rapid perturbations, and the magnitude of change matters more than the sign. These findings suggest that if the Gaia hypothesis is correct, then typical perturbations that a planet would experience may be unlikely to disrupt Gaian systems.

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.

We are on the verge of characterizing the atmospheres of terrestrial exoplanets in the habitable zones of M dwarf stars. Due to their large planet-to-star radius ratios and higher frequency of transits, terrestrial exoplanets orbiting M dwarf stars are favorable for transmission spectroscopy. In this work, we quantify the effect that water clouds have on the amplitude of water vapor transmission spectral features of terrestrial exoplanets orbiting M dwarf stars. To do so, we make synthetic transmission spectra from general circulation model (GCM) experiments of tidally locked planets. We improve upon previous work by considering how varying a broad range of planetary parameters affects transmission spectra. We find that clouds lead to a 10-100 times increase in the number of transits required to detect water features with the James Webb Space Telescope (JWST) with varying rotation period, incident stellar flux, surface pressure, planetary radius, and surface gravity. We also find that there is a strong increase in the dayside cloud coverage in our GCM simulations with rotation periods $\gtrsim 12 \ \mathrm{days}$ for planets with Earth's radius. This increase in cloud coverage leads to even stronger muting of spectral features for slowly rotating exoplanets orbiting M dwarf stars. We predict that it will be extremely challenging to detect water transmission features in the atmospheres of terrestrial exoplanets in the habitable zone of M dwarf stars with JWST. However, species that are well-mixed above the cloud deck (e.g., CO$_2$ and CH$_4$) may still be detectable on these planets with JWST.

Terrestrial planets orbiting within the habitable zones of M-stars are likely to become tidally locked in a 1:1 spin:orbit configuration and are prime targets for future characterization efforts. An issue of importance for the potential habitability of terrestrial planets is whether they could experience snowball events (periods of global glaciation). Previous work using an intermediate complexity atmospheric Global Climate Model (GCM) with no ocean heat transport suggested that tidally locked planets would smoothly transition to a snowball, in contrast with Earth, which has bifurcations and hysteresis in climate state associated with global glaciation. In this paper, we use a coupled ocean-atmosphere GCM (ROCKE-3D) to model tidally locked planets with no continents. We chose this configuration in order to consider a case that we expect to have high ocean heat transport. We show that including ocean heat transport does not reintroduce the snowball bifurcation. An implication of this result is that a tidally locked planet in the habitable zone is unlikely to be found in a snowball state for a geologically significant period of time.

According to the Snowball Earth hypothesis, Earth has experienced periods of low-latitude glaciation in its deep past. Prior studies have used general circulation models (GCMs) to examine the effects such an extreme climate state might have on the structure and dynamics of Earth's troposphere, but the behavior of the stratosphere has not been studied in detail. Understanding the snowball stratosphere is important for developing an accurate account of the Earth's radiative and chemical properties during these episodes. Here we conduct the first analysis of the stratospheric circulation of the Snowball Earth using ECHAM6 general circulation model simulations. In order to understand the factors contributing to the stratospheric circulation, we extend the Statistical Transformed Eulerian Mean framework. We find that the stratosphere during a snowball with prescribed modern ozone levels exhibits a weaker meridional overturning circulation, reduced wave activity, stronger zonal jets, and is extremely cold relative to modern conditions. Notably, the snowball stratosphere displays no sudden stratospheric warmings. Without ozone, the stratosphere displays slightly weaker circulation, a complete lack of polar vortex, and even colder temperatures. We also explicitly quantify for the first time the cross-tropopause mass exchange rate and stratospheric mixing efficiency during the snowball and show that our values do not change the constraints on CO$_2$ inferred from geochemical proxies during the Marinoan glaciation ($\sim$635 Ma), unless the O$_2$ concentration during the snowball was orders of magnitude less than the CO$_2$ concentration.

The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari (2013), and determines how unstable the atmosphere is to baroclinic instability (the baroclinic "criticality") through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes.

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.

These notes contain everything necessary to run a flipped course on "The Atmosphere" at an introductory undergraduate level. There are notes for the students to read before each course meeting and problems for them to work on in small groups during course meetings. Topics include (1) atmospheric composition, structure, and thermodynamics; (2) solar and terrestrial radiation in the atmospheric energy balance; (3) atmospheric dynamics and circulation. I include 10 problem sets, six practice midterms, and three practice finals. Problems are drawn from the atmospheres of modern and past Earth, solar system planets, and extrasolar planets. I can provide solutions to the in-class problems and problem sets to teachers upon request. You are free to use these notes in your classes, and to expand them as you please. If you catch any typos or errors, please send them to me. Enjoy!

A rigorous definition of the habitable zone and its dependence on planetary properties is part of the search for habitable exoplanets. In this work, we use the general circulation model ExoCAM to determine how the inner edge of the habitable zone of tidally locked planets orbiting M dwarf stars depends on planetary radius, surface gravity, and surface pressure. We find that the inner edge of the habitable zone for more massive planets occurs at higher stellar irradiation, as found in previous one-dimensional simulations. We also determine the relative effects of varying planetary radius and surface gravity. Increasing the planetary radius leads to a lower planetary albedo and warmer climate, pushing the inner edge of the habitable zone to lower stellar irradiation. This results from a change in circulation regime that leads to the disruption of the thick, reflective cloud deck around the substellar point. Increasing gravity increases the outgoing longwave radiation, which moves the inner edge of the habitable zone to higher stellar irradiation. This is because the column mass of water vapor decreases with increasing gravity, leading to a reduction in the greenhouse effect. The effect of gravity on the outgoing longwave radiation is stronger than the effect of radius on the planetary albedo, so that increasing gravity and radius together causes the inner edge of the habitable zone to move to higher stellar irradiation. Our results show that the inner edge of the habitable zone for more massive terrestrial planets occurs at a larger stellar irradiation.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote sensing of planetary spectroscopic and temporal signals, along with careful fitting of this data to advanced models of planets and their atmospheres. Terrestrial planets are small and significantly more challenging to observe compared to their larger gaseous brethren; however observatories coming on-line in the coming decade will begin to allow their characterization. Still, it is not enough to invest only in observational endeavors. Comprehensive modeling of planetary atmospheres is required in order to fully understand what it is that our grand telescopes see in the night-sky. In our quest to characterize habitable, and possibly inhabited worlds, 3D general circulation models (GCMs) should be used to evaluate potential climate states and their associated temporal and spatial dependent observable signals. 3D models allow for coupled, self-consistent, multi-dimensional simulations, which can realistically simulate the climates of terrestrial extrasolar planets. A complete theoretical understanding of terrestrial exoplanetary atmospheres, gained through comprehensive 3D modeling, is critical for interpreting spectra of exoplanets taken from current and planned instruments, and is critical for designing future missions that aim to measure spectra of potentially habitable exoplanets as one of their key science goals. We recommend continued institutional support for 3D GCM modeling teams that focus on planetary and exoplanetary applications.

The habitable zone (HZ) is the region around a star(s) where standing bodies of water could exist on the surface of a rocky planet. The classical HZ definition makes a number of assumptions common to the Earth, including assuming that the most important greenhouse gases for habitable planets are CO2 and H2O, habitable planets orbit main-sequence stars, and that the carbonate-silicate cycle is a universal process on potentially habitable planets. Here, we discuss these and other predictions for the habitable zone and the observations that are needed to test them. We also, for the first time, argue why A-stars may be interesting HZ prospects. Instead of relying on unverified extrapolations from our Earth, we argue that future habitability studies require first principles approaches where temporal, spatial, physical, chemical, and biological systems are dynamically coupled. We also suggest that next-generation missions are only the beginning of a much more data-filled era in the not-too-distant future, when possibly hundreds to thousands of HZ planets will yield the statistical data we need to go beyond just finding habitable zone planets to actually determining which ones are most likely to exhibit life.

This is a white paper in response to the National Academy of Sciences "Exoplanet Science Strategy" call. We summarize recent advances in theoretical habitability studies and argue that such studies will remain important for guiding and interpreting observations. Interactions between 1-D and 3-D climate modelers will be necessary to resolve recent discrepancies in model results and improve habitability studies. Observational capabilities will also need improvement. Although basic observations can be performed with present capabilities, technological advances will be necessary to improve climate models to the level needed for planetary habitability studies.

The recent detections of temperate terrestrial planets orbiting nearby stars and the promise of characterizing their atmospheres motivates a need to understand how the diversity of possible planetary parameters affects the climate of terrestrial planets. In this work, we investigate the atmospheric circulation and climate of terrestrial exoplanets orbiting both Sun-like and M-dwarf stars over a wide swath of possible planetary parameters, including the planetary rotation period, surface pressure, incident stellar flux, surface gravity, planetary radius, and cloud particle size. We do so using a general circulation model (GCM) that includes non-grey radiative transfer and the effects of clouds. The results from this suite of simulations generally show qualitatively similar dependencies of circulation and climate on planetary parameters as idealized GCMs, with quantitative differences due to the inclusion of additional model physics. Notably, we find that the effective cloud particle size is a key unknown parameter that can greatly affect the climate of terrestrial exoplanets. We confirm a transition between low and high dayside cloud coverage of synchronously rotating terrestrial planets with increasing rotation period. We determine that this cloud transition is due to eddy-driven convergence near the substellar point and should not be parameterization-dependent. Finally, we compute full-phase light curves from our simulations of planets orbiting M-dwarf stars, finding that changing incident stellar flux and rotation period affect observable properties of terrestrial exoplanets. Our GCM results can guide expectations for planetary climate over the broad range of possible terrestrial exoplanets that will be observed with future space telescopes.

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.

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.

The ice-albedo feedback on rapidly-rotating terrestrial planets in the habitable zone can lead to abrupt transitions (bifurcations) between a warm and a snowball (ice-covered) state, bistability between these states, and hysteresis in planetary climate. This is important for planetary habitability because snowball events may trigger rises in the complexity of life, but could also endanger complex life that already exists. Recent work has shown that planets tidally locked in synchronous rotation states will transition smoothly into the snowball state rather than experiencing bifurcations. Here we investigate the structure of snowball bifurcations on planets that are tidally influenced, but not synchronously rotating, so that they experience long solar days. We use PlaSIM, an intermediate-complexity global climate model, with a thermodynamic mixed layer ocean and the Sun's spectrum. We find that the amount of hysteresis (range in stellar flux for which there is bistability in climate) is significantly reduced for solar days with lengths of tens of Earth days, and disappears for solar days of hundreds of Earth days. These results suggest that tidally influenced planets orbiting M and K-stars that are not synchronously rotating could have much less hysteresis associated with the snowball bifurcations than they would if they were rapidly rotating. This implies that the amount of time it takes them to escape a snowball state via CO$_2$ outgassing would be greatly reduced, as would the period of cycling between the warm and snowball state if they have a low CO$_2$ outgassing rate.

The TRAPPIST-1, Proxima Centauri, and LHS 1140 systems are the most exciting prospects for future follow-up observations of potentially inhabited planets. All orbit nearby M-stars and are likely tidally locked in 1:1 spin-orbit states, which motivates the consideration of the effects that tidal locking might have on planetary habitability. On Earth, periods of global glaciation (snowballs) may have been essential for habitability and remote signs of life (biosignatures) because they are correlated with increases in the complexity of life and in the atmospheric oxygen concentration. In this paper we investigate the snowball bifurcation (sudden onset of global glaciation) on tidally locked planets using both an energy balance model and an intermediate-complexity global climate model. We show that tidally locked planets are unlikely to exhibit a snowball bifurcation as a direct result of the spatial pattern of insolation they receive. Instead they will smoothly transition from partial to complete ice coverage and back. A major implication of this work is that tidally locked planets with an active carbon cycle should not be found in a snowball state. Moreover, this work implies that tidally locked planets near the outer edge of the habitable zone with low CO2 outgassing fluxes will equilibrate with a small unglaciated substellar region rather than cycling between warm and snowball states. More work is needed to determine how the lack of a snowball bifurcation might affect the development of life on a tidally locked planet.

The search for habitable exoplanets and life beyond the Solar System is one of the most compelling scientific opportunities of our time. Nevertheless, the high cost of building facilities that can address this topic and the keen public interest in the results of such research requires the rigorous development of experiments that can deliver a definitive advance in our understanding. Most work to date in this area has focused on a "systems science" approach of obtaining and interpreting comprehensive data for individual planets to make statements about their habitability and the possibility that they harbor life. This strategy is challenging because of the diversity of exoplanets, both observed and expected, and the limited information that can be obtained with astronomical instruments. Here we propose a complementary approach that is based on performing surveys of key planetary characteristics and using statistical marginalization to answer broader questions than can be addressed with a small sample of objects. The fundamental principle of this comparative planetology approach is maximizing what can be learned from each type of measurement by applying it widely rather than requiring that multiple kinds of observations be brought to bear on a single object. As a proof of concept, we outline a survey of terrestrial exoplanet atmospheric water and carbon dioxide abundances that would test the habitable zone hypothesis and lead to a deeper understanding of the frequency of habitable planets. We also discuss ideas for additional surveys that could be developed to test other foundational hypotheses is this area.

Terrestrial exoplanets in the canonical habitable zone may have a variety of initial water fractions due to random volatile delivery by planetesimals. If the total planetary water complement is high, the entire surface may be covered in water, forming a "waterworld." On a planet with active tectonics, competing mechanisms act to regulate the abundance of water on the surface by determining the partitioning of water between interior and surface. Here we explore how the incorporation of different mechanisms for the degassing and regassing of water changes the volatile evolution of a planet. For all of the models considered, volatile cycling reaches an approximate steady-state after $\sim 2 \ \mathrm{Gyr}$. Using these steady-states, we find that if volatile cycling is either solely dependent on temperature or seafloor pressure, exoplanets require a high abundance ($\gtrsim 0.3\%$ of total mass) of water to have fully inundated surfaces. However, if degassing is more dependent on seafloor pressure and regassing mainly dependent on mantle temperature, the degassing rate is relatively large at late times and a steady-state between degassing and regassing is reached with a substantial surface water fraction. If this hybrid model is physical, super-Earths with a total water fraction similar to that of the Earth can become waterworlds. As a result, further understanding of the processes that drive volatile cycling on terrestrial planets is needed to determine the water fraction at which they are likely to become waterworlds.

The habitable zone concept is important because it focuses the scientific search for extraterrestrial life and aids the planning of future telescopes. Recent work has shown that planets near the outer edge of the habitable zone might not actually be able to stay warm and habitable if CO$_2$ outgassing rates are not large enough to maintain high CO$_2$ partial pressures against removal by silicate weathering. In this paper I use simple equations for the climate and CO$_2$ budget of a planet in the habitable zone that can capture the qualitative behavior of the system. With these equations I derive an analytical formula for an effective outer edge of the habitable zone, including limitations imposed by the CO$_2$ outgassing rate. I then show that climate cycles between a Snowball state and a warm climate are only possible beyond this limit if the weathering rate in the Snowball climate is smaller than the CO$_2$ outgassing rate (otherwise stable Snowball states result). I derive an analytical solution for the climate cycles including a formula for their period in this limit. This work allows us to explore the qualitative effects of weathering processes on the effective outer edge of the habitable zone, which is important because weathering parameterizations are uncertain.

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.

A terrestrial planet in an orbit far outside of the standard habitable zone could maintain surface liquid water as a result of H2-H2 collision-induced absorption by a thick H2 atmosphere. Without a stabilizing climate feedback, however, habitability would be accidental and likely brief. In this letter I propose stabilizing climate feedbacks for such a planet that require only that biological functions have an optimal temperature and operate less efficiently at other temperatures. For example, on a planet with a net source of H2 from its interior, H2-consuming life (such as methanogens) could establish a stable climate. If a positive perturbation is added to the equilibrium temperature, H2 consumption by life will increase (cooling the planet) until the equilibrium climate is reestablished. The potential existence of such feedbacks makes H2-warmed planets more attractive astrobiological targets.

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.

Surface liquid water is essential for standard planetary habitability. Calculations of atmospheric circulation on tidally locked planets around M stars suggest that this peculiar orbital configuration lends itself to the trapping of large amounts of water in kilometers-thick ice on the night side, potentially removing all liquid water from the day side where photosynthesis is possible. We study this problem using a global climate model including coupled atmosphere, ocean, land, and sea-ice components as well as a continental ice sheet model driven by the climate model output. For a waterworld we find that surface winds transport sea ice toward the day side and the ocean carries heat toward the night side. As a result, night-side sea ice remains O(10 m) thick and night-side water trapping is insignificant. If a planet has large continents on its night side, they can grow ice sheets O(1000 m) thick if the geothermal heat flux is similar to Earth's or smaller. Planets with a water complement similar to Earth's would therefore experience a large decrease in sea level when plate tectonics drives their continents onto the night side, but would not experience complete day-side dessication. Only planets with a geothermal heat flux lower than Earth's, much of their surface covered by continents, and a surface water reservoir O(10 %) of Earth's would be susceptible to complete water trapping.

Planetary rotation rate is a key parameter in determining atmospheric circulation and hence the spatial pattern of clouds. Since clouds can exert a dominant control on planetary radiation balance, rotation rate could be critical for determining mean planetary climate. Here we investigate this idea using a three-dimensional general circulation model with a sophisticated cloud scheme. We find that slowly rotating planets (like Venus) can maintain an Earth-like climate at nearly twice the stellar flux as rapidly rotating planets (like Earth). This suggests that many exoplanets previously believed to be too hot may actually be habitable, depending on their rotation rate. The explanation for this behavior is that slowly rotating planets have a weak Coriolis force and long daytime illumination, which promotes strong convergence and convection in the substellar region. This produces a large area of optically thick clouds, which greatly increases the planetary albedo. In contrast, on rapidly rotating planets a much narrower belt of clouds form in the deep tropics, leading to a relatively low albedo. A particularly striking example of the importance of rotation rate suggested by our simulations is that a planet with modern Earth's atmosphere, in Venus' orbit, and with modern Venus' (slow) rotation rate would be habitable. This would imply that if Venus went through a runaway greenhouse, it had a higher rotation rate at that time.

In the spirit of minimal modeling of complex systems, we develop an idealized two-column model to investigate the climate of tidally locked terrestrial planets with Earth-like atmospheres in the habitable zone of M-dwarf stars. The model is able to approximate the fundamental features of the climate obtained from three-dimensional (3D) atmospheric general circulation model (GCM) simulations. One important reason for the two-column model's success is that it reproduces the high cloud albedo of the GCM simulations, which reduces the planet's temperature and delays the onset of a runaway greenhouse state. The two-column model also clearly illustrates a secondary mechanism for determining the climate: the nightside acts as a ``radiator fin'' through which infrared energy can be lost to space easily. This radiator fin is maintained by a temperature inversion and dry air on the nightside, and plays a similar role to the subtropics on modern Earth. Since 1D radiative-convective models cannot capture the effects of the cloud albedo and radiator fin, they are systematically biased towards a narrower habitable zone. We also show that cloud parameters are most important for determining the day--night thermal emission contrast in the two-column model, which decreases and eventually reverses as the stellar flux increases. This reversal is important because it could be detected by future extrasolar planet characterization missions, which would suggest that the planet has Earth-like water clouds and is potentially habitable.

Large terrestrial planets are expected to have muted topography and deep oceans, implying that most super-Earths should be entirely covered in water, so-called waterworlds. This is important because waterworlds lack a silicate weathering thermostat so their climate is predicted to be less stable than that of planets with exposed continents. In other words, the continuously habitable zone for waterworlds is much narrower than for Earth-like planets. A planet's water is partitioned, however, between a surface reservoir, the ocean, and an interior reservoir, the mantle. Plate tectonics transports water between these reservoirs on geological timescales. Degassing of melt at mid-ocean ridges and serpentinization of oceanic crust depend negatively and positively on seafloor pressure, respectively, providing a stabilizing feedback on long-term ocean volume. Motivated by Earth's approximately steady-state deep water cycle, we develop a two-box model of the hydrosphere and derive steady-state solutions to the water partitioning on terrestrial planets. Critically, hydrostatic seafloor pressure is proportional to surface gravity, so super-Earths with a deep water cycle will tend to store more water in the mantle. We conclude that a tectonically active terrestrial planet of any mass can maintain exposed continents if its water mass fraction is less than ~0.2%, dramatically increasing the odds that super-Earths are habitable. The greatest source of uncertainty in our study is Earth's current mantle water inventory: the greater its value, the more robust planets are to inundation. Lastly, we discuss how future missions can test our hypothesis by mapping the oceans and continents of massive terrestrial planets.

Planets in M dwarf stars' habitable zones are likely to be tidally locked with orbital periods of order tens of days. This means that the effects of rotation on atmospheric dynamics will be relatively weak, which requires small horizontal temperature gradients above the boundary layer of terrestrial atmospheres. An analytically solvable and dynamically consistent model for planetary climate with only three free parameters can be constructed by making the weak temperature gradient (WTG) approximation, which assumes temperatures are horizontally uniform aloft. The extreme numerical efficiency of a WTG model compared to a 3D general circulation model (GCM) makes it an optimal tool for Monte Carlo fits to observables over parameter space. Additionally, such low-order models are critical for developing physical intuition and coupling atmospheric dynamics to models of other components of planetary climate. The objective of this paper is to determine whether a WTG model provides an adequate approximation of the effect of atmospheric dynamics on quantities likely to be observed over the next decade. To do this we first tune a WTG model to GCM output for an Earth-like tidally locked planet with a dry, 1 bar atmosphere, then generate and compare the expected phase curves of both models. We find that differences between the two models would be extremely difficult to detect from phase curves using JWST. This result demonstrates the usefulness of the WTG approximation when used in conjunction with GCMs as part of a modeling hierarchy to understand the climate of remote planets.

The habitable zone (HZ) is the circumstellar region where a planet can sustain surface liquid water. Searching for terrestrial planets in the HZ of nearby stars is the stated goal of ongoing and planned extrasolar planet surveys. Previous estimates of the inner edge of the HZ were based on one-dimensional radiative-convective models. The most serious limitation of these models is the inability to predict cloud behavior. Here we use global climate models with sophisticated cloud schemes to show that due to a stabilizing cloud feedback, tidally locked planets can be habitable at twice the stellar flux found by previous studies. This dramatically expands the HZ and roughly doubles the frequency of habitable planets orbiting red dwarf stars. At high stellar flux, strong convection produces thick water clouds near the substellar location that greatly increase the planetary albedo and reduce surface temperatures. Higher insolation produces stronger substellar convection and therefore higher albedo, making this phenomenon a stabilizing climate feedback. Substellar clouds also effectively block outgoing radiation from the surface, reducing or even completely reversing the thermal emission contrast between dayside and nightside. The presence of substellar water clouds and the resulting clement surface conditions will therefore be detectable with the James Webb Space Telescope.

In order to understand the climate on terrestrial planets orbiting nearby Sun-like stars, one would like to know their thermal inertia. We use a global climate model to simulate the thermal phase variations of Earth-analogs and test whether these data could distinguish between planets with different heat storage and heat transport characteristics. In particular, we consider a temperate climate with polar ice caps (like modern Earth), and a snowball state where the oceans are globally covered in ice. We first quantitatively study the periodic radiative forcing from, and climatic response to, rotation, obliquity, and eccentricity. Orbital eccentricity and seasonal changes in albedo cause variations in the global-mean absorbed flux. The responses of the two climates to these global seasons indicate that the temperate planet has 3 times the bulk heat capacity of the snowball planet due to the presence of liquid water oceans. The temperate obliquity seasons are weaker than one would expect based on thermal inertia alone; this is due to cross-equatorial oceanic and atmospheric energy transport. Thermal inertia and cross-equatorial heat transport have qualitatively different effects on obliquity seasons, insofar as heat transport tends to reduce seasonal amplitude without inducing a phase lag. For an Earth-like planet, however, this effect is masked by the mixing of signals from low thermal inertia regions (sea ice and land) with that from high thermal inertia regions (oceans), which also produces a damped response with small phase lag. We then simulate thermal lightcurves as they would appear to a high-contrast imaging mission (TPF-I/Darwin) and consider the inverse problem of estimating thermal inertia based solely on time-resolved photometry. [Abridged]

It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be discovered. Water content helps determine surface land fraction, which influences planetary weathering behavior. This is important because the silicate weathering feedback determines the width of the habitable zone in space and time. Here a low-order model of weathering and climate, useful for gaining qualitative understanding, is developed to examine climate evolution for planets of various land-ocean fractions. It is pointed out that, if seafloor weathering does not depend directly on surface temperature, there can be no weathering-climate feedback on a waterworld. This would dramatically narrow the habitable zone of a waterworld. Results from our model indicate that weathering behavior does not depend strongly on land fraction for partially ocean-covered planets. This is powerful because it suggests that previous habitable zone theory is robust to changes in land fraction, as long as there is some land. Finally, a mechanism is proposed for a waterworld to prevent complete water loss during a moist greenhouse through rapid weathering of exposed continents. This process is named a "waterworld self-arrest," and it implies that waterworlds can go through a moist greenhouse stage and end up as planets like Earth with partial ocean coverage. This work stresses the importance of surface and geologic effects, in addition to the usual incident stellar flux, for habitability.

Identifying liquid water on the surface of planets is a high priority, as this traditionally defines habitability. One proposed signature of oceans is specular reflection ("glint"), which increases the apparent albedo of a planet at crescent phases. We post-process a global climate model of an Earth-like planet to simulate reflected lightcurves. Significantly, we obtain glint-like phase variations even though we do not include specular reflection in our model. This false positive is the product of two generic properties: 1) for modest obliquities, a planet's poles receive less orbit-averaged stellar flux than its equator, so the poles are more likely to be covered in highly reflective snow and ice, and 2) we show that reflected light from a modest-obliquity planet at crescent phases probes higher latitudes than at gibbous phases, therefore a planet's apparent albedo will naturally increase at crescent phase. We suggest that this "latitude-albedo effect" will operate even for large obliquities: in that case the equator receives less orbit-averaged flux than the poles, and the equator is preferentially sampled at crescent phase. Using rotational and orbital color variations to map the surfaces of directly imaged planets and estimate their obliquity will therefore be a necessary pre-condition for properly interpreting their reflected phase variations. The latitude-albedo effect is a particularly convincing glint false positive for zero-obliquity planets, and such worlds are not amenable to latitudinal mapping. This effect severely limits the utility of specular reflection for detecting oceans on exoplanets.

Rogue planets have been ejected from their planetary system. We investigate the possibility that a rogue planet could maintain a liquid ocean under layers of thermally-insulating water ice and frozen atmosphere as a result of geothermal heat flux. We find that a rogue planet of Earth-like composition and age could maintain a subglacial liquid ocean if it were ~3.5 times more massive than Earth. If a rogue planet had about ten times higher water mass fraction or a thick cryo-atmospheric layer, it would need to be only ~0.3 times the mass of Earth to maintain a liquid ocean. Such a planet could be detected from reflected solar radiation and its thermal emission could be characterized in the far-IR if it passed within O(1000) AU of Earth.