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.

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.

Discoveries of terrestrial, Earth-sized exoplanets that lie within the habitable zone (HZ) of their host stars continue to occur at increasing rates. Transit spectroscopy can potentially enable the detection of molecular signatures from such worlds, providing an indication of the presence of an atmosphere and its chemical composition, including gases potentially indicative of a biosphere. Such planets around nearby M-dwarf stars - such as TRAPPIST-1 - provide relatively good signal, high signal/noise, and frequent transits for follow-up spectroscopy. However, even with these advantages, transit spectroscopy of terrestrial planets in the HZ of nearby M-stars will still be a challenge. Herein, we examine the potential for future space observatories to conduct such observations, using a Global Climate Model (GCM), a photochemical model, and a radiative transfer suite to simulate modern-Earth-like atmospheric boundary conditions on TRAPPIST-1e. The detectability of biosignatures on such an atmosphere via transmission spectroscopy is modeled for JWST, LUVOIR, HabEx, and Origins. We show that for any of these observatories, only CO2 would be detectable at the 3 sigma level in transmission spectroscopy, when clouds are included in our simulations. This is because the impacts of clouds on scale height strongly limits the detectability of molecules in the atmosphere. Synergies between space- and ground-based spectroscopy may be essential in order to overcome these difficulties.

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.

The collision-induced fundamental vibration-rotation band at 6.4 um is the most significant absorption feature from O2 in the infrared (Timofeyev and Tonkov, 1978; Rinslandet al., 1982, 1989), yet it has not been previously incorporated into exoplanet spectral analyses for several reasons. Either CIAs were not included or incomplete/obsolete CIA databases were used. Also, the current version of HITRAN does not include CIAs at 6.4 um with other collision partners (O2-X). We include O2-X CIA features in our transmission spectroscopy simulations by parameterizing the 6.4 um O2-N2 CIA based on Rinsland et al.(1989) and the O2-CO2 CIA based on Baranov et al. (2004). Here we report that the O2-X CIA may be the most detectable O2 feature for transit observations. For a potentialTRAPPIST-1e analogue system within 5 pc of the Sun, it could be the only O2 detectable signature with JWST (using MIRI LRS) for a modern Earth-like cloudy atmosphere with biological quantities of O2. Also, we show that the 6.4 um O2-X CIA would be prominent for O2-rich desiccated atmospheres (Luger and Barnes, 2015) and could be detectable with JWST in just a few transits. For systems beyond 5 pc, this feature could therefore be a powerful discriminator of uninhabited planets with non-biological "false positive" O2 in their atmospheres - as they would only be detectable at those higher O2 pressures.

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.

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 similar to Earth - but slightly more irradiated - are expected to enter into a runaway greenhouse state, where all surface water rapidly evaporates, forming an optically thick H2O-dominated atmosphere. For Earth, this extreme climate transition is thought to occur for a ~6% increase only of the solar luminosity, though the exact limit at which the transition would occur is still a highly debated topic. In general, the runaway greenhouse is believed to be a fundamental process in the evolution of Earth-size, temperate planets. Using 1-D radiative-convective climate calculations accounting for thick, hot water vapour-dominated atmospheres, we evaluate the transit atmospheric thickness of a post-runaway greenhouse atmosphere, and find that it could possibly reach over a thousand kilometers (i.e., a few tens of % of Earth radius). This abrupt radius inflation - resulting from the runaway-greenhouse-induced transition - could be detected statistically by ongoing and upcoming space missions such as TESS, CHEOPS and PLATO (combined with precise radial velocity mass measurements with ground-based spectrographs such as ESPRESSO, CARMENES or SPIRou), or even in particular cases in multiplanetary systems such as TRAPPIST-1 (when masses and radii will be known with good enough precision). This result provides the community with an observational test of (1) the concept of runaway greenhouse, that defines the inner edge of the traditional Habitable Zone, and the exact limit of the runaway greenhouse transition. In particular, this could provide an empirical measurement of the irradiation at which Earth analogs transition from a temperate to a runaway greenhouse climate state. This astronomical measurement would make it possible to statistically estimate how close Earth is from the runaway greenhouse. (2) the presence (and statistical abundance) of water in temperate, Earth-size exoplanets.

We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system.

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets.

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.

1D climate models are less sophisticated than 3D global circulation models (GCMs), however their computational time is much less expensive, allowing a large number of runs in a short period of time to explore a wide parameter space. Exploring parameter space is particularly important for predicting the observable properties of exoplanets, for which few parameters are known with certainty. Therefore, 1D climate models are still very useful tools for planetary studies. In most of these 1D models, clouds are not physically represented in the atmosphere, despite having a well-known, significant impact on a planetary radiative budget. This impact is simulated by artificially raising surface albedo, in order to reproduce the observed-averaged surface temperature (i.e. 288 K for modern Earth) and a radiative balance at the top of the atmosphere. This non-physical representation of clouds, causes atmospheric longwave and shortwaves fluxes to not match observational data. Additionally, this technique represents a parameter that is highly-tuned to modern Earth's climate, and may not be appropriate for planets that deviate from modern Earth's climate conditions. In this paper, we present an update to the climate model within the Atmos 1D atmospheric modeling package with a physical representation of clouds. We show that this physical representation of clouds in the atmosphere allows both longwave and shortwave fluxes to match observational data. This improvement will allow us to study the energy fluxes for a variety of cloudy rocky planets, and increase our confidence in future simulations of temperature profile and net energy balance.

Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmosphere and surface can be derived by analyzing light of the parent star reflected by the planet. We investigate the influence of the surface albedo $A_{\rm s}$, the optical thickness $b_{\rm cloud}$ and altitude of water clouds, and the mixing ratio $\eta$ of biosignature O$_2$ on the strength of the O$_2$ A-band (around 760 nm) in flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ($\eta$ < 0.4) will yield moderately deep bands in flux and moderate to small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on their properties such as optical thickness, top altitude, particle phase, coverage fraction, horizontal distribution. Depending on the surface albedo, and cloud properties, different O$_2$ mixing ratios $\eta$ can give similar absorption band depths in flux and band strengths in polarization, in particular if the clouds have moderate to high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, in particular not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O$_2$ or any other absorbing gas.

Asteroid modeling efforts in the last decade resulted in a comprehensive dataset of almost 400 convex shape models and their rotation states. This amount already provided a deep insight into physical properties of main-belt asteroids or large collisional families. We aim to increase the number of asteroid shape models and rotation states. Such results are an important input for various further studies such as analysis of asteroid physical properties in different populations, including smaller collisional families, thermophysical modeling, and scaling shape models by disk-resolved images, or stellar occultation data. This provides, in combination with known masses, bulk density estimates, but constrains also theoretical collisional and evolutional models of the Solar System. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as an input for the convex inversion method, and derive 3D shape models of asteroids, together with their rotation periods and orientations of rotation axes. The key ingredient is the support of more that one hundred observers who submit their optical data to publicly available databases. We present updated shape models for 36 asteroids, for which mass estimates are currently available in the literature or their masses will be most likely determined from their gravitational influence on smaller bodies, which orbital deflection will be observed by the ESA Gaia astrometric mission. This was achieved by using additional optical data from recent apparitions for the shape optimization. Moreover, we also present new shape model determinations for 250 asteroids, including 13 Hungarias and 3 near-Earth asteroids.