Search

Discoveries of giant planet candidates orbiting white dwarf stars and the demonstrated capabilities of the James Webb Space Telescope bring the possibility of detecting rocky planets in the habitable zones of white dwarfs into pertinent focus. We present simulations of an aqua planet with an Earth-like atmospheric composition and incident stellar insolation orbiting in the habitable zone of two different types of stars - a 5000 K white dwarf and main-sequence K-dwarf star Kepler-62 with a similar effective temperature - and identify the mechanisms responsible for the two differing planetary climates. The synchronously-rotating white dwarf planet's global mean surface temperature is 25 K higher than that of the synchronously-rotating planet orbiting Kepler-62, due to its much faster (10-hr) rotation and orbital period. This ultra-fast rotation generates strong zonal winds and meridional flux of zonal momentum, stretching out and homogenizing the scale of atmospheric circulation, and preventing an equivalent build-up of thick, liquid water clouds on the dayside of the planet compared to the synchronous planet orbiting Kepler-62, while also transporting heat equatorward from higher latitudes. White dwarfs may therefore present amenable environments for life on planets formed within or migrated to their habitable zones, generating warmer surface environments than those of planets with main-sequence hosts to compensate for an ever shrinking incident stellar flux., Comment: 18 pages, 8 figures, 2 tables, accepted for publication in The Astrophysical Journal

This paper highlights methods from geostatistics that are relevant to the interpretation, intercomparison, and synthesis of atmospheric model data, with a specific application to exoplanet atmospheric modeling. Climate models are increasingly used to study theoretical and observational properties of exoplanets, which include a hierarchy of models ranging from fast and idealized models to those that are slower but more comprehensive. Exploring large parameter spaces with computationally-expensive models can be accomplished with sparse sampling techniques, but analyzing such sparse samples can pose challenges for conventional interpolation functions. Ordinary kriging is a statistical method for describing the spatial distribution of a data set in terms of the variogram function, which can be used to interpolate sparse samples across any number of dimensions. Variograms themselves may also be useful diagnostic tools for describing the spatial distribution of model data in exoplanet atmospheric model intercomparison projects. Universal kriging is another method that can synthesize data calculated by models of different complexity, which can be used to combine sparse samples of data from slow models with larger samples of data from fast models. Ordinary and universal kriging can also provide a way to synthesize model predictions with sparse samples of exoplanet observations and may have other applications in exoplanet science., Comment: Accepted for publication in the Planetary Science Journal

Potentially habitable exoplanets are targets of great interest for the James Webb Space Telescope and upcoming mission concepts such as the Habitable Worlds Observatory. Clouds strongly affect climate and habitability, but predicting their properties is difficult. In Global Climate Models (GCMs), especially those aiming at simulating Earth, cloud microphysics is often crudely approximated by assuming that all cloud particles have a single, constant size or a prescribed size distribution and that all clouds in a grid cell are identical. For exoplanets that range over a large phase space of planetary properties, this method could result in large errors. In this work, our goal is to determine how cloud microphysics on terrestrial exoplanets, whose condensable is mainly water vapor, depend on aerosol properties and planetary parameters such as surface pressure, surface gravity, and incident stellar radiation. We use the Community Aerosol and Radiation Model for Atmospheres as a 1D microphysical model to simulate the formation and evolution of clouds including the processes of nucleation, condensation, evaporation, coagulation, and vertical transfer. In these 1D idealized experiments, we find that the parameters that determine the macrophysical thermal structure of the atmospheres, including surface pressure and stellar flux, impact cloud radiative effect (CRE) most significantly. Parameters such as gravity and number density of aerosols working as cloud condensation nuclei affect the microphysical processes of cloud formation, including activation and vertical transfer. They also have a significant, though weaker effect on CRE. This work motivates the development of more accurate GCM cloud schemes and should aid in the interpretation of future observations., Comment: 17 pages, 10 figures

Methane is thought to have been an important greenhouse gas during the Archean, although its potential warming has been found to be limited at high concentrations due to its high shortwave absorption. We use the Met Office Unified Model, a general circulation model, to further explore the climatic effect of different Archean methane concentrations. Surface warming peaks at a pressure ratio CH$_4$:CO$_2$ of approximately 0.1, reaching a maximum of up to 7 K before significant cooling above this ratio. Equator-to-pole temperature differences also tend to increase up to pCH$_4$ $\leq$300 Pa, which is driven by a difference in radiative forcing at the equator and poles by methane and a reduction in the latitudinal extend of the Hadley circulation. 3D models are important to fully capture the cooling effect of methane, due to these impacts of the circulation., Comment: 36 pages, 18 figures

With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors., Comment: Accepted to PSJ

Planets in synchronous rotation around low-mass stars are the most salient targets for current ground- and space-based missions to observe and characterize. Such model calculations can help to prioritize targets for observation with current and future missions; however, intrinsic differences in the complexity and physical parameterizations of various models can lead to different predictions of a planet's climate state. Understanding such model differences is necessary if such models are to guide target selection and aid in the analysis of observations. This paper presents a protocol to intercompare models of a hypothetical planet with a 15 day synchronous rotation period around a 3000 K blackbody star across a parameter space of surface pressure and incident instellation. We conduct a sparse sample of 16 cases from a previously published exploration of this parameter space with the ExoPlaSim model. By selecting particular cases across this broad parameter space, the SAMOSA intercomparison will identify areas where simpler models are sufficient as well as areas where more complex GCMs are required. Our preliminary comparison using ExoCAM shows general consistency between the climate state predicted by ExoCAM and ExoPlaSim except in regions of the parameter space most likely to be in a steam atmosphere or incipient runaway greenhouse state. We use this preliminary analysis to define several options for participation in the intercomparison by models of all levels of complexity. The participation of other GCMs is crucial to understand how the atmospheric states across this parameter space differ with model capabilities., Comment: Accepted for publication in the Planetary Science Journal

A key item of interest for planetary scientists and astronomers is the habitable zone, or the distance from a host star where a terrestrial planet can maintain necessary temperatures in order to retain liquid water on its surface. However, when observing a system's habitable zone, it is possible that one may instead observe a Venus-like planet. We define "Venus-like" as greenhouse-gas-dominated atmosphere occurring when incoming solar radiation exceeds infrared radiation emitted from the planet at the top of the atmosphere, resulting in a runaway greenhouse. Our definition of Venus-like includes both incipient and post-runaway greenhouse states. Both the possibility of observing a Venus-like world and the possibility that Venus could represent an end-state of evolution for habitable worlds, requires an improved understanding of the Venus-like planet; specifically, the distances where these planets can exist. Understanding this helps us define a "Venus zone", or the region in which Venus-like planets could exist, and assess the overlap with the aforementioned "Habitable Zone". In this study, we use a 1D radiative-convective climate model to determine the outer edge of the Venus zone for F0V, G2V, K5V, and M3V and M5V stellar spectral types. Our results show that the outer edge of the Venus zone resides at 3.01, 1.36, 0.68, 0.23, and 0.1 AU, respectively. These correspond to incident stellar fluxes of 0.8, 0.55, 0.38, 0.32, and 0.3 S, respectively, where stellar flux is relative to Earth (1.0). These results indicate that there may be considerable overlap between the habitable zone and the Venus zone., Comment: Accepted for publication in the Planetary Science Journal

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) is a community project that aims to quantify how dfferences in general circulation models (GCMs) could impact the climate prediction for TRAPPIST-1e and, subsequently its atmospheric characterization in transit. Four GCMs have participated in THAI so far: ExoCAM, LMD-Generic, ROCKE-3D and the UM. This paper, focused on the simulated observations, is the third part of a trilogy, following the analysis of two land planet scenarios (part I) and two aquaplanet scenarios (part II). Here, we show a robust agreement between the simulated spectra and the number of transits estimated to detect the land planet atmospheres. For the aquaplanet ones, using atmospheric data from any of the four GCMs would require at least 17 transits. This prediction corresponds to UM simulated data which produces the lowest and thinnest clouds. Between 35-40% more clouds are predicted by ExoCAM or LMD-G due to higher thick terminator clouds. For the first time this work provides "GCM uncertainty error bars" of 35-40% that need to be considered in future analyses of transmission spectra. We also analyzed the inter-transit variability induced by weather patterns and changes of terminator cloudiness between transits. Its magnitude differs significantly between the GCMs but its impact on the transmission spectra is within the measurement uncertainties. THAI has demonstrated the importance of model intercomparison for exoplanets and also paved the way for a larger project to develop an intercomparison meta-framework, namely the Climates Using Interactive Suites of Intercomparisons Nested for Exoplanet Studies (CUISINES)., Comment: Accepted in the Planetary Science Journal as Part III of a series of 3 THAI papers

To identify promising exoplanets for atmospheric characterization and to make the best use of observational data, a thorough understanding of their atmospheres is needed. 3D general circulation models (GCMs) are one of the most comprehensive tools available for this task and will be used to interpret observations of temperate rocky exoplanets. Due to parameterization choices made in GCMs, they can produce different results, even for the same planet. Employing four widely-used exoplanetary GCMs -- ExoCAM, LMD-G, ROCKE-3D and the UM -- we continue the TRAPPIST-1 Habitable Atmosphere Intercomparison by modeling aquaplanet climates of TRAPPIST-1e with a moist atmosphere dominated by either nitrogen or carbon dioxide. Although the GCMs disagree on the details of the simulated regimes, they all predict a temperate climate with neither of the two cases pushed out of the habitable state. Nevertheless, the inter-model spread in the global mean surface temperature is non-negligible: 14 K and 24 K in the nitrogen and carbon dioxide dominated case, respectively. We find substantial inter-model differences in moist variables, with the smallest amount of clouds in LMD-Generic and the largest in ROCKE-3D. ExoCAM predicts the warmest climate for both cases and thus has the highest water vapor content and the largest amount and variability of cloud condensate. The UM tends to produce colder conditions, especially in the nitrogen-dominated case due to a strong negative cloud radiative effect on the day side of TRAPPIST-1e. Our study highlights various biases of GCMs and emphasizes the importance of not relying solely on one model to understand exoplanet climates., Comment: 35 pages, 22 figures; Published in The Planetary Science Journal

With the commissioning of powerful, new-generation telescopes such as the James Webb Space Telescope (JWST) and the ground-based Extremely Large Telescopes, the first characterization of a high molecular weight atmosphere around a temperate rocky exoplanet is imminent. Atmospheric simulations and synthetic observables of target exoplanets are essential to prepare and interpret these observations. Here we report the results of the first part of the TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) project, which compares 3D numerical simulations performed with four state-of-the-art global climate models (ExoCAM, LMD-Generic, ROCKE-3D, Unified Model) for the potentially habitable target TRAPPIST-1e. In this first part, we present the results of dry atmospheric simulations. These simulations serve as a benchmark to test how radiative transfer, subgrid-scale mixing (dry turbulence and convection), and large-scale dynamics impact the climate of TRAPPIST-1e and consequently the transit spectroscopy signature as seen by JWST. To first order, the four models give results in good agreement. The intermodel spread in the global mean surface temperature amounts to 7K (6K) for the N2-dominated (CO2-dominated) atmosphere. The radiative fluxes are also remarkably similar (intermodel variations less than 5%), from the surface (1 bar) up to atmospheric pressures around 5 mbar. Moderate differences between the models appear in the atmospheric circulation pattern (winds) and the (stratospheric) thermal structure. These differences arise between the models from (1) large-scale dynamics, because TRAPPIST-1e lies at the tipping point between two different circulation regimes (fast and Rhines rotators) in which the models can be alternatively trapped, and (2) parameterizations used in the upper atmosphere such as numerical damping., Comment: Published in the Planetary Science Journal as Part I of a series of 3 THAI papers

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison (THAI) workshop (14-16 September 2020). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements and limitations, as well as direction for future model development was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community., Comment: Accepted in the AAS Planetary Science Journal (PSJ)

Low-mass stars show evidence of vigorous magnetic activity in the form of large flares and coronal mass ejections. Such space weather events may have important ramifications for the habitability and observational fingerprints of exoplanetary atmospheres. Here, using a suite of three-dimensional coupled chemistry-climate model (CCM) simulations, we explore effects of time-dependent stellar activity on rocky planet atmospheres orbiting G-, K-, and M-dwarf stars. We employ observed data from the MUSCLES campaign and Transiting Exoplanet Satellite Survey and test a range of rotation period, magnetic field strength, and flare frequency assumptions. We find that recurring flares drive K- and M-dwarf planet atmospheres into chemical equilibria that substantially deviate from their pre-flare regimes, whereas G-dwarf planet atmospheres quickly return to their baseline states. Interestingly, simulated O$_2$-poor and O$_2$-rich atmospheres experiencing flares produce similar mesospheric nitric oxide abundances, suggesting that stellar flares can highlight otherwise undetectable chemical species. Applying a radiative transfer model to our CCM results, we find that flare-driven transmission features of bio-indicating species, such as nitrogen dioxide, nitrous oxide, and nitric acid, show particular promise for detection by future instruments., Comment: 27 pages, 14 figures total, 6 main text figures, 8 extended data figures, updated to the published version in the March 2021 Issue of Nature Astronomy; Volume 5, Pages 298-310; the published version is available at https://www.nature.com/articles/s41550-020-01264-1

Much attention has been given to the climate dynamics and habitable boundaries of synchronously rotating planets around low mass stars. However, other rotational states are possible, particularly when higher eccentricity orbits can be maintained in a system, including spin-orbit resonant configurations. Additionally, the oscillating strain as a planet moves from periastron to apoastron results in friction and tidal heating, which can be an important energy source. Here, we simulate the climate of ocean-covered planets near the inner edge of the habitable zone around M to solar stars with ROCKE-3D, and leverage the planetary evolution software package, VPLanet, to calculate tidal heating rates for Earth-sized planets orbiting 2600 K and 3000 K stars. This study is the first to use a 3-D General Circulation Model that implements tidal heating to investigate habitability for multiple resonant states. We find that in the absence of tidal heating, the resonant state has little impact on the inner edge, because for a given stellar flux, higher-order states tend to be warmer than synchronous rotators, but for a given temperature, have drier upper atmospheres. However, when strong tidal heating is present, the rotational component implies a strong dependence of habitable conditions on the system evolution and rotational state. Since tidal and stellar heating both decrease with orbital distance, this results in a compact orbital width separating temperate and uninhabitable climates. We summarize these results and also compare ROCKE-3D to previously published simulations of the inner edge that used a modified version of the NCAR CAM4 model.

The discovery of terrestrial exoplanets is uncovering increasingly diverse architectures. Of particular interest are those systems that contain exoplanets at a variety of star-planet separations, allowing direct comparison of exoplanet evolution (comparative planetology). The Kepler-1649 system contains two terrestrial planets similar both in size and insolation flux to Venus and Earth, although their eccentricities remain largely unconstrained. Here we present results of dynamical studies of the system and the potential effects on climate. The eccentricities of the Kepler-1649 system are poorly constrained, and we show that there are dynamically viable regions for further terrestrial planets in between the two known planets for a limited range of eccentricities. We investigate the effect of eccentricity of the outer planet on the dynamics of both planets and show that this results in high-frequency (1000-3000 year) eccentricity oscillations in long-term stable configurations. We calculate the resulting effect of these eccentricity variations on insolation flux and present the results of 3D climate simulations for the habitable zone planet. Our simulations demonstrate that, despite large eccentricity variations, the planet can maintain stable climates with relatively small temperature variations on the substellar hemisphere for a variety of initial climate configurations. Such systems thus provide key opportunities to explore alternative Venus/Earth climate evolution scenarios., Comment: 12 pages, 6 figures, accepted for publication in the Astronomical Journal

Variations in the reflective properties of the bulk material that comprises the surface of land-dominated planets will affect the planetary energy balance by interacting differently with incident radiation from the host star. Furthermore, low-mass cool stars, such as nearby M8V dwarf TRAPPIST-1, emit a significant fraction of their flux in longer wavelengths relative to the Sun in regions where terrestrial materials may exhibit additional variability in albedo. Using the Community Earth System Model (CESM) we investigate the effect of the composition of the land surface and its albedo on planetary climate in the context of spatially homogeneous, entirely land-covered planets with dry atmospheres at the orbital separation of TRAPPIST-1d, TRAPPIST-1e, and TRAPPIST-1f. We use empirically derived spectra of four terrestrial compositional endmembers (granite, calcite, aridisol, and dune sand) and a composite spectrum of TRAPPIST-1 for these simulations and compare these model output to an aquaplanet and several Sol-spectrum control cases. We report a difference of approximately 50 K in global mean surface temperature, variations in atmospheric rotational features, and a reduction in cross-equatorial heat transport between scenarios in which materials with higher albedo in the infrared (calcite and dune sand) were used and those with more absorptive crustal material, such as granite or dry soils. An aquaplanet TRAPPIST-1d scenario results in an unstable runaway greenhouse regime. Therefore, we demonstrate that determining the composition and albedo of continental landmasses is crucial for making accurate determinations of the climate of terrestrial exoplanets., Comment: 18 pages, 11 figures. Accepted for publication in the Astrophysical Journal (ApJ)

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., Comment: 55 pages, 11 figures, 2 tables. Review article accepted for publication in Space Science Reviews. Part of ISSI special collection on 'Understanding the Diversity of Planetary Atmospheres'. Abstract abridged to meet ArXiv size limit

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., Comment: 22 pages, 7 figures, 1 table. Accepted for publication in Space Science Reviews. Part of ISSI special collection on Diversity of Atmospheres

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., Comment: 35 pages, 21 figures, 6 tables

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., Comment: Accepted and published in Geoscientific Model Development (GMD)

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., Comment: Accepted to the Astronomical Journal. 21 pages, 11 figures

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., Comment: Accepted for publication in AJ (30 pages, 13 figures, 4 tables, 2 appendices)

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., Comment: 34 pages, 17 figures

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., Comment: 7 pages, 3 figures, accepted at ApJ Letters

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., Comment: Accepted to the Astrophysical Journal. 25 pages, 10 figures

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., Comment: ApJ accepted

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., Comment: 36 pages, 19 figures

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., Comment: Accepted for publication in JGR-Planets

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., Comment: Book chapter for "Planetary Astrobiology" (Space science series)

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., Comment: Accepted at ApJ, 8 pages, 4 figures

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., Comment: 23 pages, 11 figures, and 2 tables. Published in the Astrophysical Journal on 14 November 2019

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., Comment: 26 pages, 15 figures, 3 tables. Accepted for publication in the Astronomical Journal

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., Comment: White paper submitted in response to the solicitation of feedback for the "2020 Decadal Survey" by the National Academy of Sciences

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., Comment: NAS ASTRO 2020, White Paper! 2 figures, and 10 pages

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., Comment: Astrophysical Journal in review; 11 pages; 4 figures; 2 tables

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., Comment: 12 pages, 4 figures, 1 table; published in the Astrophysical Journal Letters on 15 November 2018; corrected two small typos since published version on 26 December 2018

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., Comment: 25 pages, 1 table, and 10 figures

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., Comment: A white paper submitted to the National Academy of Sciences Exoplanet Science Strategy

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., Comment: This is a white paper submitted to the National Academies 2018 "Exoplanet Science Strategy" call(6 pages total, including cover page). Corrected references section in this version. http://sites.nationalacademies.org/SSB/CurrentProjects/SSB_180659

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., Comment: 33 pages, 11 figures. The update corrects the definition in Eq. (8) and of beta. The results are unaffected as these were transcription errors only and did not influence the calculations. Errata are published at doi:10.3847/1538-4357/ab9a4b and doi:10.3847/1538-4357/abb19a

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., Comment: Accepted to Astrophysical Journal

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., Comment: Published in International Journal of Astrobiology

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., Comment: Submitted to ApJL, 5 figures, New version includes errata appended to the end

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., Comment: accepted at Astrophysical Journal, 3 tables, 6 figures, 43 pages

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., Comment: 21 pages, 10 figures, published in ApJ (836, 49)

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., Comment: 111 pages, 15 figures, 4 tables, accepted for publication in Astrobiology

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., Comment: Submitted to Monthly Notices of the Royal Astronomical Society

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., Comment: Accepted to Astrophysical Journal. 31 pages, 11 figures