We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

TRAPPIST-1 is a fantastic nearby (~39.14 light years) planetary system made of at least seven transiting terrestrial-size, terrestrial-mass planets all receiving a moderate amount of irradiation. To date, this is the most observationally favourable system of potentially habitable planets. Since the announcement of the discovery of TRAPPIST-1 planets in 2016, a growing number of techniques and approaches have been used and proposed to reveal its true nature. Here we have compiled a state-of-the-art overview of all the observational and theoretical constraints that have been obtained so far using these techniques and approaches. The goal is to get a better understanding of whether or not TRAPPIST-1 planets can have atmospheres, and if so, what they are made of. For this, we surveyed the literature on TRAPPIST-1 about topics as broad as irradiation environment, orbital architecture, transit observations, density measurements, stellar contamination, and numerical climate and escape models. Each of these topics adds a brick to our understanding of the likely atmospheres of the seven planets. We show that (i) HST transit observations, (ii) density measurements, (iii) atmospheric escape modelling, and (iv) gas accretion modelling altogether offer solid evidence against the presence of H2-dominated atmospheres around TRAPPIST-1 planets. This means they likely have either (i) a high molecular weight atmosphere or (ii) no atmosphere at all. There are several key challenges ahead to characterize the bulk compositions of the atmospheres (if present) of TRAPPIST-1 planets. The main one so far is characterizing and correcting for the effects of stellar contamination. Fortunately, a new wave of observations with the James Webb Space Telescope and near-infrared high-resolution ground-based spectrographs on very large telescopes will bring significant advances in the coming decade.

We present self-consistent three-dimensional climate simulations of possible habitable states for the newly discovered Habitable Zone Earth-sized planet, TOI-700 d. We explore a variety of atmospheric compositions, pressures, and rotation states for both ocean-covered and completely desiccated planets in order to assess the planet's potential for habitability. For all 20 of our simulated cases, we use our climate model outputs to synthesize transmission spectra, combined-light spectra, and integrated broadband phase curves. These climatologically-informed observables will help the community assess the technological capabilities necessary for future characterization of this planet - as well as similar transiting planets discovered in the future - and will provide a guide for distinguishing possible climate states if one day we do obtain sensitive spectral observations of a habitable planet around a M-star. We find that TOI-700 d is a strong candidate for a habitable world and can potentially maintain temperate surface conditions under a wide variety of atmospheric compositions. Unfortunately, the spectral feature depths from the resulting transmission spectra and the peak flux and variations from our synthesized phase curves for TOI-700 d do not exceed 10 ppm. This will likely prohibit the James Webb Space Telescope (JWST) from characterizing its atmosphere; however, this motivates the community to invest in future instrumentation that perhaps can one day reveal the true nature of TOI-700 d, and to continue to search for similar planets around less distant stars.

Stellar evolution models predict that the solar luminosity was lower in the past, typically 20-25 % lower during the Archean (3.8-2.5 Ga). Despite the fainter Sun, there is strong evidence for the presence of liquid water on Earth's surface at that time. This "faint young Sun problem" is a fundamental question in paleoclimatology, with important implications for the habitability of the early Earth, early Mars and exoplanets. Many solutions have been proposed based on the effects of greenhouse gases, atmospheric pressure, clouds, land distribution and Earth's rotation rate. Here we review the faint young Sun problem for Earth, highlighting the latest geological and geochemical constraints on the early Earth's atmosphere, and recent results from 3D global climate models and carbon cycle models. Based on these works, we argue that the faint young Sun problem for Earth has essentially been solved. Unfrozen Archean oceans were likely maintained by higher concentrations of CO2, consistent with the latest geological proxies, potentially helped by additional warming processes. This reinforces the expected key role of the carbon cycle for maintaining the habitability of terrestrial planets. Additional constraints on the Archean atmosphere and 3D fully coupled atmosphere-ocean models are required to validate this conclusion.

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well.

The search for water-rich Earth-sized exoplanets around low-mass stars is rapidly gaining attention because they represent the best opportunity to characterize habitable planets in the near future. Understanding the atmospheres of these planets and determining the optimal strategy for characterizing them through transmission spectroscopy with our upcoming instrumentation is essential in order to constrain their environments. For this study, we present simulated transmission spectra of tidally locked Earth-sized ocean-covered planets around late-M to mid-K stellar spectral types, utilizing GCM modeling results previously published by Kopparapu et al. (2017) as inputs for our radiative transfer calculations performed using NASA's Planetary Spectrum Generator (psg.gsfc.nasa.gov; Villanueva et al. (2018)). We identify trends in the depth of H$_2$O spectral features as a function of planet surface temperature and rotation rate. These trends allow us to calculate the exposure times necessary to detect water vapor in the atmospheres of aquaplanets through transmission spectroscopy with the upcoming James Webb Space Telescope (JWST) as well as several future flagship space telescope concepts under consideration (LUVOIR and OST) for a target list constructed from the TESS Input Catalog (TIC). Our calculations reveal that transmission spectra for water-rich Earth-sized planets around low-mass stars will be dominated by clouds, with spectral features < 20 ppm, and only a small subset of TIC stars would allow for the characterization of an ocean planet in the Habitable Zone. We thus present a careful prioritization of targets that are most amenable to follow-up characterizations with next-generation instrumentation, in order to assist the community in efficiently utilizing precious telescope time.

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.

Planets residing in circumstellar habitable zones (CHZs) offer our best opportunities to test hypotheses of life's potential pervasiveness and complexity. Constraining the precise boundaries of habitability and its observational discriminants is critical to maximizing our chances at remote life detection with future instruments. Conventionally, calculations of the inner edge of the habitable zone (IHZ) have been performed using both 1D radiative-convective and 3D general circulation models. However, these models lack interactive three-dimensional chemistry and do not resolve the mesosphere and lower thermosphere (MLT) region of the upper atmosphere. Here we employ a 3D high-top chemistry-climate model (CCM) to simulate the atmospheres of synchronously-rotating planets orbiting at the inner edge of habitable zones of K- and M-dwarf stars (between $T_{\rm eff} =$ 2600 K and 4000 K). While our IHZ climate predictions are in good agreement with GCM studies, we find noteworthy departures in simulated ozone and HO$_{\rm x}$ photochemistry. For instance, climates around inactive stars do not typically enter the classical moist greenhouse regime even with high ($< 10^{-3}$ mol mol$^{-1}$) stratospheric water vapor mixing ratios, which suggests that planets around inactive M-stars may only experience minor water-loss over geologically significant timescales. In addition, we find much thinner ozone layers on potentially habitable moist greenhouse atmospheres, as ozone experiences rapid destruction via reaction with hydrogen oxide radicals. Using our CCM results as inputs, our simulated transmission spectra show that both water vapor and ozone features on moist greenhouse atmospheres could be detectable by instruments NIRSpec and MIRI LRS onboard the James Webb Space Telescope.

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.

Small exoplanets of nearby M dwarf stars present the possibility to find and characterize habitable worlds within the next decade. TRAPPIST-1, an ultracool M dwarf star, was recently found to have seven Earth-sized planets of predominantly rocky composition. The planets e, f, and g can have a liquid water ocean on their surface given appropriate atmospheres of N2 and CO2. Particularly, climate models have shown that the planets e and f can sustain a global liquid water ocean, for >=0.2 bar CO2 plus 1 bar N2, or >=2 bars CO2, respectively. These atmospheres are irradiated by ultraviolet emission from the star's moderately active chromosphere, and the consequence of this irradiation is unknown. Here we show that chemical reactions driven by the irradiation can produce and maintain more than 0.2 bar O2 and 0.05 bar CO if the CO2 is >=0.1 bar. The abundance of O2 and CO can rise to more than 1 bar under certain boundary conditions. Because of this O2-CO runaway, habitable environments on the TRAPPIST-1 planets entail an O2- and CO-rich atmosphere with coexisting O3. The only process that would prevent the runaway is direct recombination of O2 and CO in the ocean, a reaction that is facilitated biologically. Our results indicate that O2, O3, and CO should be considered together with CO2 as the primary molecules in the search for atmospheric signatures from temperate and rocky planets of TRAPPIST-1 and other M dwarf stars.

The TRAPPIST-1 system, consisting of an ultra-cool host star having seven known Earth-size planets will be a prime target for atmospheric characterization with JWST. However, the detectability of atmospheric molecular species may be severely impacted by the presence of clouds and/or hazes. In this work, we perform 3-D General Circulation Model (GCM) simulations with the LMD Generic model supplemented by 1-D photochemistry simulations at the terminator with the Atmos model to simulate several possible atmospheres for TRAPPIST-1e, 1f and 1g: 1) modern Earth, 2) Archean Earth, and 3) CO2-rich atmospheres. JWST synthetic transit spectra were computed using the GSFC Planetary Spectrum Generator (PSG). We find that TRAPPIST-1e, 1f and 1g atmospheres, with clouds and/or hazes, could be detected using JWST's NIRSpec prism from the CO2 absorption line at 4.3 um in less than 15 transits at 3 sigma or less than 35 transits at 5 sigma. However, our analysis suggests that other gases would require hundreds (or thousands) of transits to be detectable. We also find that H2O, mostly confined in the lower atmosphere, is very challenging to detect for these planets or similar systems if the planets' atmospheres are not in a moist greenhouse state. This result demonstrates that the use of GCMs, self-consistently taking into account the effect of clouds and sub-saturation, is crucial to evaluate the detectability of atmospheric molecules of interest as well as for interpreting future detections in a more global (and thus robust and relevant) approach.

Planets that revolve around a binary pair of stars are known as circumbinary planets. The orbital motion of the stars around their center of mass causes a periodic variation in the total instellation incident upon a circumbinary planet. This study uses both an analytic and numerical energy balance model to calculate the extent to which this effect can drive changes in surface temperature on circumbinary terrestrial planets. We show that the amplitude of the temperature variation is largely constrained by the effective heat capacity, which corresponds to the ocean-to-land ratio on the planet. Planets with large ocean fractions should experience only modest warming and cooling of only a few degrees, which suggests that habitability cannot be precluded for such circumbinary planets. Planets with large land fractions that experience extreme periodic forcing could be prone to changes in temperature of tens of degrees or more, which could drive more extreme climate changes that inhibit continuously habitable conditions.

Habitability is a measure of an environment's potential to support life, and a habitable exoplanet supports liquid water on its surface. However, a planet's success in maintaining liquid water on its surface is the end result of a complex set of interactions between planetary, stellar, planetary system and even Galactic characteristics and processes, operating over the planet's lifetime. In this chapter, we describe how we can now determine which exoplanets are most likely to be terrestrial, and the research needed to help define the habitable zone under different assumptions and planetary conditions. We then move beyond the habitable zone concept to explore a new framework that looks at far more characteristics and processes, and provide a comprehensive survey of their impacts on a planet's ability to acquire and maintain habitability over time. We are now entering an exciting era of terrestiral exoplanet atmospheric characterization, where initial observations to characterize planetary composition and constrain atmospheres is already underway, with more powerful observing capabilities planned for the near and far future. Understanding the processes that affect the habitability of a planet will guide us in discovering habitable, and potentially inhabited, planets.

We use a one-dimensional (1-D) cloud-free climate model to estimate habitable zone (HZ) boundaries for terrestrial planets of masses 0.1 M$_{E}$ and 5 M$_{E}$ around circumbinary stars of various spectral type combinations. Specifically, we consider binary systems with host spectral types F-F, F-G, F-K, F-M, G-G, G-K, G-M, K-K, K-M and M-M. Scaling the background N2 atmospheric pressure with the radius of the planet, we find that the inner edge of the HZ moves inwards towards the star for 5ME compared to 0.1ME planets for all spectral types. This is because the water-vapor column depth is smaller for larger planets and higher temperatures are needed before water vapor completely dominates the outgoing longwave radiation. The outer edge of the HZ changes little due to competing effects of the albedo and greenhouse effect. While these results are broadly consistent with the trend of single star HZ results for different mass planets, there are significant differences between single star and binary star systems for the inner edge of the HZ. Interesting combinations of stellar pairs from our 1-D model results can be used to explore for in-depth climate studies with 3-D climate models. We identify a common HZ stellar flux domain for all circumbinary spectral types

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.

This White Paper presents a mission concept called MIRECLE - the Mid-InfraRed Exoplanet CLimate Explorer. With a moderately sized aperture of 2 meters, broad wavelength coverage (4 - 25 um), and next generation instruments, MIRECLE will be capable of efficiently characterizing a statistically significant sample of terrestrial planets, many of which will be in their host stars's habitable zones. Spectroscopic characterization of terrestrial atmospheres will provide constraints for the distribution of planets with tenuous vs. substantial atmospheres, on the inner and outer edges of the habitable zone, and climate models to assess the potential for habitability. For the few brightest targets, the detection of specific combinations of molecules would provide evidence of biosignatures. For all other targets, this comprehensive survey would filter out the airless, desiccated, or lifeless worlds, thus providing a subset of potentially habitable worlds ready for in-depth atmospheric characterization using a larger aperture telescope.

Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location of water vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.

The lunar history of water deposition, loss, and transport post-accretion has become an important consideration in relation to the possibility of a human outpost on the Moon. Very recent work has shown that a secondary primordial atmosphere of up to 10 mb could have been emplaced ~3.5 billion years ago due to volcanic outgassing from the maria. Using a zero dimensional chemistry model we demonstrate the temperature dependence of the resulting major atmospheric components (CO or CO$_2$). We use a three dimensional general circulation model to test the viability of such an atmosphere and derive its climatological characteristics. Based on these results we then conjecture on its capability to transport volatiles outgassed from the maria to the permanently shadowed regions at the poles. Our preliminary results demonstrate that atmospheres as low as 1 mb are viable and that permanent cold trapping of volatiles is only possible at the poles.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Rotation and orbital eccentricity both strongly influence planetary climate. Eccentricities can often be measured for exoplanets, but rotation rates are currently difficult or impossible to constrain. Here we examine how the combined effects of rotation and eccentricity on observed emission from ocean-rich terrestrial planets can be used to infer their rotation rates in circumstances where their eccentricities are known. We employ an Earth climate model with no land and a slab ocean, and consider two eccentricities ($e=0.3$ and 0.6) and two rotation rates: a fast Earth-like period of 24 hours, and a slower pseudo-synchronous period that generalizes spin synchronization for eccentric orbits. We adopt bandpasses of the Mid-Infrared Instrument on the James Webb Space Telescope as a template for future photometry. At $e=0.3$ the rotation rates can be distinguished if the planet transits near periastron, because slow rotation produces a strong day-night contrast and thus an emission minimum during periastron. However, light curves behave similarly if the planet is eclipsed near periastron, as well as for either viewing geometry at $e=0.6$. Rotation rates can nevertheless be distinguished using ratios of emission in different bands, one in the water vapor window with another in a region of strong water absorption. These ratios vary over an orbit by $\lesssim\!0.1$ dex for Earth-like rotation, but by 0.3--0.5 dex for pseudo-synchronous rotation because of large day-night contrast in upper-tropospheric water. For planets with condensible atmospheric constituents in eccentric orbits, rotation regimes might thus be distinguished with infrared observations for a range of viewing geometries.

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.

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.

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.

Transit spectroscopy of terrestrial planets around nearby M dwarfs is a primary goal of space missions in coming decades. 3-D climate modeling has shown that slow-synchronous rotating terrestrial planets may develop thick clouds at the substellar point, increasing the albedo. For M dwarfs with Teff > 3000 K, such planets at the inner habitable zone (IHZ) have been shown to retain moist greenhouse conditions, with enhanced stratospheric water vapor (fH2O > 1E-3) and low Earth-like surface temperatures. However, M dwarfs also possess strong UV activity, which may effectively photolyze stratospheric H2O. Prior modeling efforts have not included the impact of high stellar UV activity on the H2O. Here, we employ a 1-D photochemical model with varied stellar UV, to assess whether H2O destruction driven by high stellar UV would affect its detectability in transmission spectroscopy. Temperature and water vapor profiles are taken from published 3-D climate model simulations for an IHZ Earth-sized planet around a 3300 K M dwarf with an N2-H2O atmosphere; they serve as self-consistent input profiles for the 1-D model. We explore additional chemical complexity within the 1-D model by introducing other species into the atmosphere. We find that as long as the atmosphere is well-mixed up to 1 mbar, UV activity appears to not impact detectability of H2O in the transmission spectrum. The strongest H2O features occur in the JWST MIRI instrument wavelength range and are comparable to the estimated systematic noise floor of ~50 ppm.

A planet's atmospheric constituents (e.g., O$_2$, O$_3$, H$_2$O$_v$, CO$_2$, CH$_4$, N$_2$O) can provide clues to its surface habitability, and may offer biosignature targets for remote life detection efforts. The plethora of rocky exoplanets found by recent transit surveys (e.g., the Kepler mission) indicates that potentially habitable systems orbiting K- and M-dwarf stars may have very different orbital and atmospheric characteristics than Earth. To assess the physical distribution and observational prospects of various biosignatures and habitability indicators, it is important to understand how they may change under different astrophysical and geophysical configurations, and to simulate these changes with models that include feedbacks between different subsystems of a planet's climate. Here we use a three-dimensional (3D) Chemistry-Climate model (CCM) to study the effects of changes in stellar spectral energy distribution (SED), stellar activity, and planetary rotation on Earth-analogs and tidally-locked planets. Our simulations show that, apart from shifts in stellar SEDs and UV radiation, changes in illumination geometry and rotation-induced circulation can influence the global distribution of atmospheric biosignatures. We find that the stratospheric day-to-night side mixing ratio differences on tidally-locked planets remain low ($<20\%$) across the majority of the canonical biosignatures. Interestingly however, secondary photosynthetic biosignatures (e.g., C$_2$H$_6$S) show much greater (${\sim}67\%$) day-to-night side differences, and point to regimes in which tidal-locking could have observationally distinguishable effects on phase curve, transit, and secondary eclipse measurements. Overall, this work highlights the potential and promise for 3D CCMs to study the atmospheric properties and habitability of terrestrial 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.

The TRAPPIST-1 system provides an extraordinary opportunity to study multiple terrestrial extrasolar planets and their atmospheres. Here we use the National Center for Atmospheric Research Community Atmosphere Model version 4 to study the possible climate and habitability of the planets in the TRAPPIST-1 system. We assume ocean-covered worlds, with atmospheres comprised of N2, CO2, and H2O, and with orbital and geophysical properties defined from observation. Model results indicate that the inner three planets (b, c, and d) presently reside interior to the inner edge of the traditional liquid water habitable zone. Thus if water ever existed on the inner planets, they would have undergone a runaway greenhouse and lost their water to space, leaving them dry today. Conversely the outer 3 planets (f, g, and h) fall beyond the maximum CO2 greenhouse outer edge of the habitable zone. Model results indicate that the outer planets cannot be warmed despite as much as 30 bar CO2 atmospheres, instead entering a snowball state. The middle planet (e) represents the best chance for a presently habitable ocean-covered world in the TRAPPIST-1 system. Planet e can maintain at least some habitable surface area with 0 - 2 bar CO2, depending on the background N2 content. Near present day Earth surface temperatures can be maintained for an ocean-covered planet e with either 1 bar N2 and 0.4 bar CO2, or a 1.3 bar pure CO2 atmosphere.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

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

Human influence on the biosphere has been evident at least since the development of widespread agriculture, and some stratigraphers have suggested that the activities of modern civilization indicate a geological epoch transition. The study of the anthropocene as a geological epoch, and its implication for the future of energy-intensive civilizations, is an emerging transdisciplinary field in which astrobiology can play a leading role. Habitability research of Earth, Mars, and exoplanets examines extreme cases relevant for understanding climate change as a planetary process. Energy-intensive civilizations will also face thermodynamic limits to growth, which provides an important constraint for estimating the longevity of human civilization and guiding the search for extraterrestrial intelligence. We recommend that missions concepts such as LUVOIR, HabEx, and OST be pursued in order to make significant progress toward understanding the future evolution of life on our planet and the possible evolution of technological, energy-intensive life elsewhere in the universe.

We investigate the atmospheric dynamics of terrestrial planets in synchronous rotation within the habitable zone of low-mass stars using the Community Atmosphere Model (CAM). The surface temperature contrast between day and night hemispheres decreases with an increase in incident stellar flux, which is opposite the trend seen on gas giants. We define three dynamical regimes in terms of the equatorial Rossby deformation radius and the Rhines length. The slow rotation regime has a mean zonal circulation that spans from day to night side, with both the Rossby deformation radius and the Rhines length exceeding planetary radius, which occurs for planets around stars with effective temperatures of 3300 K to 4500 K (rotation period > 20 days). Rapid rotators have a mean zonal circulation that partially spans a hemisphere and with banded cloud formation beneath the substellar point, with the Rossby deformation radius is less than planetary radius, which occurs for planets orbiting stars with effective temperatures of less than 3000 K (rotation period < 5 days). In between is the Rhines rotation regime, which retains a thermally-direct circulation from day to night side but also features midlatitude turbulence-driven zonal jets. Rhines rotators occur for planets around stars in the range of 3000 K to 3300 K (rotation period ~ 5 to 20 days), where the Rhines length is greater than planetary radius but the Rossby deformation radius is less than planetary radius. The dynamical state can be observationally inferred from comparing the morphology of the thermal emission phase curves of synchronously rotating planets.

Terrestrial planets in the habitable zones (HZs) of low-mass stars and cool dwarfs have received significant scrutiny recently because their shorter orbital periods increase their chances of detection and characterization compared to planets around G-dwarfs. As these planets are likely tidal-locked, improved 3D numerical simulations of such planetary atmospheres are needed to guide target selection. Here we use a 3-D climate system model, updated with new water-vapor absorption coefficients derived from the HITRAN 2012 database, to study ocean covered planets at the inner edge of the HZ around late-M to mid-K stars ($2600$ K <= Teff <= 4500K). Our results indicate that these updated water-vapor coefficients result in significant warming compared to previous studies, so the inner HZ around M-dwarfs is not as close as suggested by earlier work. Assuming synchronously rotating planets with background 1 bar N2 atmospheres, we find that planets at the inner HZ of stars with Teff > 3000K undergo the classical "moist-greenhouse" (H2O mixing ratio > 10-3 in the stratosphere) at significantly lower surface temperature (~ 280K) in our 3-D model compared with 1-D climate models (~ 340K). This implies that some planets around low mass stars can simultaneously undergo water-loss and remain habitable. However, for star with Teff <= 3000K, planets at the inner HZ may directly transition to a runaway state, while bypassing the moist greenhouse water-loss entirely. We analyze transmission spectra of planets in a moist green-house regime, and find that there are several prominent H2O features within JWST instruments range. Thus, relying only upon standard Earth-analog spectra with 24-hour rotation period around M-dwarfs for habitability studies will miss the strong H2O features that one would expect to see on synchronously rotating planets around M-dwarf stars, with JWST.

M-dwarf stars are more abundant than G-dwarf stars, so our position as observers on a planet orbiting a G-dwarf raises questions about the suitability of other stellar types for supporting life. If we consider ourselves as typical, in the anthropic sense that our environment is probably a typical one for conscious observers, then we are led to the conclusion that planets orbiting in the habitable zone of G-dwarf stars should be the best place for conscious life to develop. But such a conclusion neglects the possibility that K-dwarfs or M-dwarfs could provide more numerous sites for life to develop, both now and in the future. In this paper we analyze this problem through Bayesian inference to demonstrate that our occurrence around a G-dwarf might be a slight statistical anomaly, but only the sort of chance event that we expect to occur regularly. Even if M-dwarfs provide more numerous habitable planets today and in the future, we still expect mid G- to early K-dwarfs stars to be the most likely place for observers like ourselves. This suggests that observers with similar cognitive capabilities as us are most likely to be found at the present time and place, rather than in the future or around much smaller stars.

Conventional definitions of habitability require abundant liquid surface water to exist continuously over geologic timescales. Water in each of its thermodynamic phases interacts with solar and thermal radiation and is the cause for strong climatic feedbacks. Thus, assessments of the habitable zone require models to include a complete treatment of the hydrological cycle over geologic time. Here, we use the Community Atmosphere Model from the National Center for Atmospheric Research to study the evolution of climate for an Earth-like planet at constant CO2, under a wide range of stellar fluxes from F-, G-, and K-dwarf main sequence stars. Around each star we find four stable climate states defined by mutually exclusive global mean surface temperatures (Ts); snowball (Ts < 235 K), waterbelt (235 K < Ts < 250 K), temperate (275 K < Ts < 315 K), and moist greenhouse (Ts > 330 K). Each is separated by abrupt climatic transitions. Waterbelt, temperate, and cooler moist greenhouse climates can maintain open-ocean against both sea-ice albedo and hydrogen escape processes respectively, and thus constitute habitable worlds. We consider the warmest possible habitable planet as having Ts ~ 355 K, at which point diffusion limited water-loss could remove an Earth ocean in ~1 Gyr. Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO2, habitability can be maintained for an upper limit of ~2.2, ~2.4 and ~4.7 Gyr around F-, G- and K-dwarf stars respectively due to main sequence brightening.

Hazes are common in known planet atmospheres, and geochemical evidence suggests early Earth occasionally supported an organic haze with significant environmental and spectral consequences. The UV spectrum of the parent star drives organic haze formation through methane photochemistry. We use a 1D photochemical-climate model to examine production of fractal organic haze on Archean Earth-analogs in the habitable zonesof several stellar types: the modern and early Sun, AD Leo (M3.5V), GJ 876 (M4V), $\epsilon$ Eridani (K2V), and $\sigma$ Bo\"otis (F2V). For Archean-like atmospheres, planets orbiting stars with the highest UV fluxes do not form haze due to the formation of photochemical oxygen radicals that destroy haze precursors. Organic hazes impact planetary habitability via UV shielding and surface cooling, but this cooling is minimized around M dwarfs whose energy is emitted at wavelengths where organic hazes are relatively transparent. We generate spectra to test the detectability of haze. For 10 transits of a planet orbiting GJ 876 observed by the James Webb Space Telescope, haze makes gaseous absorption features at wavelengths $<$ 2.5 $\mu$m 2-10$\sigma$ shallower compared to a haze-free planet, and methane and carbon dioxide are detectable at $>$5$\sigma$. A haze absorption feature can be detected at 5$\sigma$ near 6.3 $\mu$m, but higher signal-to-noise is needed to distinguish haze from adjacent absorbers. For direct imaging of a planet at 10 parsecs using a coronagraphic 10-meter class ultraviolet-visible-near infrared telescope, a UV-blue haze absorption feature would be strongly detectable at $>$12$\sigma$ in 200 hours.

Recognizing whether a planet can support life is a primary goal of future exoplanet spectral characterization missions, but past research on habitability assessment has largely ignored the vastly different conditions that have existed in our planet's long habitable history. This study presents simulations of a habitable yet dramatically different phase of Earth's history, when the atmosphere contained a Titan-like organic-rich haze. Prior work has claimed a haze-rich Archean Earth (3.8-2.5 billion years ago) would be frozen due to the haze's cooling effects. However, no previous studies have self-consistently taken into account climate, photochemistry, and fractal hazes. Here, we demonstrate using coupled climate-photochemical-microphysical simulations that hazes can cool the planet's surface by about 20 K, but habitable conditions with liquid surface water could be maintained with a relatively thick haze layer (tau ~ 5 at 200 nm) even with the fainter young sun. We find that optically thicker hazes are self-limiting due to their self-shielding properties, preventing catastrophic cooling of the planet. Hazes may even enhance planetary habitability through UV shielding, reducing surface UV flux by about 97% compared to a haze-free planet, and potentially allowing survival of land-based organisms 2.6.2.7 billion years ago. The broad UV absorption signature produced by this haze may be visible across interstellar distances, allowing characterization of similar hazy exoplanets. The haze in Archean Earth's atmosphere was strongly dependent on biologically-produced methane, and we propose hydrocarbon haze may be a novel type of spectral biosignature on planets with substantial levels of CO2. Hazy Archean Earth is the most alien world for which we have geochemical constraints on environmental conditions, providing a useful analog for similar habitable, anoxic exoplanets.

The population of known extrasolar planets includes giant and terrestrial planets that closely orbit their host star. Such planets experience significant tidal distortions that can force the planet into synchronous rotation. The combined effects of tidal deformation and centripetal acceleration induces significant asphericity in the shape of these planets, compared to the mild oblateness of Earth, with maximum gravitational acceleration at the poles. Here we show that this latitudinal variation in gravitational acceleration is relevant for modeling the climate of oblate planets including Jovian planets within the solar system, closely-orbiting hot Jupiters, and planets within the habitable zone of white dwarfs. We compare first- and third-order approximations for gravitational acceleration on an oblate spheroid and calculate the geostrophic wind that would result from this asphericity on a range of solar system planets and exoplanets. Third-order variations in gravitational acceleration are negligible for Earth but become significant for Jupiter, Saturn, and Jovian exoplanets. This latitudinal variation in gravitational acceleration can be measured remotely, and the formalism presented here can be implemented for use in general circulation climate modeling studies of exoplanet atmospheres.

Terrestrial planets at the inner edge of the habitable zone of late-K and M-dwarf stars are expected to be in synchronous rotation, as a consequence of strong tidal interactions with their host stars. Previous global climate model (GCM) studies have shown that, for slowly-rotating planets, strong convection at the substellar point can create optically thick water clouds, increasing the planetary albedo, and thus stabilizing the climate against a thermal runaway. However these studies did not use self-consistent orbital/rotational periods for synchronously rotating planets placed at different distances from the host star. Here we provide new estimates of the inner edge of the habitable zone for synchronously rotating terrestrial planets around late-K and M-dwarf stars using a 3-D Earth-analog GCM with self-consistent relationships between stellar metallicity, stellar effective temperature, and the planetary orbital/rotational period. We find that both atmospheric dynamics and the efficacy of the substellar cloud deck are sensitive to the precise rotation rate of the planet. Around mid-to-late M-dwarf stars with low metallicity, planetary rotation rates at the inner edge of the HZ become faster, and the inner edge of the habitable zone is farther away from the host stars than in previous GCM studies. For an Earth-sized planet, the dynamical regime of the substellar clouds begins to transition as the rotation rate approaches ~10 days. These faster rotation rates produce stronger zonal winds that encircle the planet and smear the substellar clouds around it, lowering the planetary albedo, and causing the onset of the water-vapor greenhouse climatic instability to occur at up to ~25% lower incident stellar fluxes than found in previous GCM studies. For mid-to-late M-dwarf stars with high metallicity and for mid-K to early-M stars, we agree with previous studies.