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.

Mass-radius relationships for water-rich rocky planets are usually calculated assuming most water is present in condensed (either liquid or solid) form. Planet density estimates are then compared to these mass-radius relationships, even when these planets are more irradiated than the runaway greenhouse irradiation limit (around 1.1~times the insolation at Earth for planets orbiting a Sun-like star), for which water has been shown to be unstable in condensed form and would instead form a thick H2O-dominated atmosphere. Here we use the LMD Generic numerical climate model to derive new mass-radius relationships appropriate for water-rich rocky planets that are more irradiated than the runaway greenhouse irradiation limit, meaning planets endowed with a steam, water-dominated atmosphere. For a given water-to-rock mass ratio, these new mass-radius relationships lead to planet bulk densities much lower than calculated when water is assumed to be in condensed form. In other words, using traditional mass-radius relationships for planets that are more irradiated than the runaway greenhouse irradiation limit tends to dramatically overestimate -- possibly by several orders of magnitude -- their bulk water content. In particular, this result applies to TRAPPIST-1 b, c, and d, which can accommodate a water mass fraction of at most 2, 0.3 and 0.08 %, respectively, assuming planetary core with a terrestrial composition. In addition, we show that significant changes of mass-radius relationships (between planets less and more irradiated than the runaway greenhouse limit) can be used to remove bulk composition degeneracies in multiplanetary systems such as TRAPPIST-1. Finally, we provide an empirical formula for the H2O steam atmosphere thickness which can be used to construct mass-radius relationships for any water-rich, rocky planet more irradiated than the runaway greenhouse irradiation threshold.

With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the project. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 $\mu$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $\mu$m to constrain the brightness temperatures of their daysides. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We estimate for TRAPPIST-1 transit depth measurements mean noise floors of $\sim$35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 $\mu$m, and can only set 3$\sigma$ upper limits on their dayside brightness temperatures (611K for b 586K for c).

Terrestrial planets in temperate orbits around very low mass stars are likely to have evolved in a very different way than solar system planets, and in particular Earth. However, because these are the first planets that are and will be accessible for in-depth atmosphere, clouds and surface characterizations with existing and forthcoming telescopes, we need to develop the best possible observational strategies to maximize the scientific return from these characterizations. Here I discuss and expand on the recent works of Bean et al. (2017) and Turbet et al. (2019) to show that terrestrial planets orbiting in temperate orbits around very low mass stars are potentially an excellent sample of planets to test how universal the processes thought to control the habitability of solar system planets and in particular Earth are. Precise measurements of density or atmospheric CO2 concentration for planets located both inside and outside the Habitable Zone could be used to statistically test habitability concepts such as the silicate-weathering feedback, CO2 condensation, or runaway greenhouse, which have been identified as key processes controlling the present and past habitability of Venus, Mars and Earth.

Reducing atmospheres have recently emerged as a promising scenario to warm the surface of early Mars enough to drive the formation of valley networks and other ancient aqueous features that have been detected so far on the surface of Mars. Here we present a series of experiments and calculations to better constrain CO2+CH4 and CO2+H2 collision-induced absorptions (CIAs) as well as their effect on the prediction of early Mars surface temperature. First, we carried out a new set of experimental measurements (using the AILES line of the SOLEIL synchrotron) of both CO2+CH4 and CO2+H2 CIAs. These measurements confirm the previous results of Turbet et al. 2019, Icarus vol. 321, while significantly reducing the experimental uncertainties. Secondly, we fitted a semi-empirical model to these CIAs measurements, allowing us to compute the CO2+CH4 and CO2+H2 CIAs across a broad spectral domain (0-1500cm-1) and for a wide range of temperatures (100-600K). Last, we performed 1-D numerical radiative-convective climate calculations (using the LMD Generic Model) to compute the surface temperature expected on the surface of early Mars for several CO2, CH4 and H2 atmospheric contents, taking into account the radiative effect of these revised CIAs. These calculations demonstrate that thick CO2+H2-dominated atmospheres remain a viable solution for warming the surface of Mars above the melting point of water, but not CO2+CH4-dominated atmospheres. Our calculated CO2+CH4 and CO2+H2 CIA spectra and predicted early Mars surface temperatures are provided to the community for future uses.

We use a hierarchy of numerical models (a 3-D Global Climate Model, a 1-D radiative-convective model and a 2-D Mantle Dynamics model) to explore the environmental effects of very large impacts on the atmosphere, surface and interior of early Mars. Using a combination of 1-D and 3-D climate simulations, we show that the environmental effects of the largest impact events recorded on Mars are characterized by: (i) a short impact-induced warm period; (ii) a low amount of hydrological cycling of water; (iii) deluge-style precipitation; and (iv) precipitation patterns that are uncorrelated with the observed regions of valley networks. We show that the impact-induced stable runaway greenhouse state predicted by Segura et al. 2012 is physically inconsistent. We confirm the results of Segura et al. 2008 and Urata & Toon 2013 that water ice clouds can significantly extend the duration of the post-impact warm period, and even for cloud coverage lower than predicted in Ramirez & Kasting 2017. However, the range of cloud microphysical properties for which this scenario works is very narrow. Using 2-D Mantle Dynamics simulations we find that large impacts can raise the near-surface internal heat flux up to several hundreds of mW/m$^2$ (i.e. up to $\sim$ 10 times the ambient flux) for several millions years at the edges of the impact crater. However, such internal heat flux is insufficient to keep the martian surface above the melting point of water. Our numerical results support the prediction of Palumbo & Head 2018 that very large impact-induced rainfall could have caused degradation of craters and formed smooth plains, potentially erasing much of the previously visible morphological surface history. Such hot rainfalls may have also led to the formation of aqueous alteration products on Noachian-aged terrains.

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 Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission's primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS's observations were focused on the southern ecliptic hemisphere, resulting in the discovery of three mini-Neptunes orbiting the star TOI-125, a V=11.0 K0 dwarf. We present intensive HARPS radial velocity observations, yielding precise mass measurements for TOI-125b, TOI-125c and TOI-125d. TOI-125b has an orbital period of 4.65 days, a radius of $2.726 \pm 0.075 ~\mathrm{R_{\rm E}}$, a mass of $ 9.50 \pm 0.88 ~\mathrm{M_{\rm E}}$ and is near the 2:1 mean motion resonance with TOI-125c at 9.15 days. TOI-125c has a similar radius of $2.759 \pm 0.10 ~\mathrm{R_{\rm E}}$ and a mass of $ 6.63 \pm 0.99 ~\mathrm{M_{\rm E}}$, being the puffiest of the three planets. TOI-125d, has an orbital period of 19.98 days and a radius of $2.93 \pm 0.17~\mathrm{R_{\rm E}}$ and mass $13.6 \pm 1.2 ~\mathrm{M_{\rm E}}$. For TOI-125b and TOI-125d we find unusual high eccentricities of $0.19\pm 0.04$ and $0.17^{+0.08}_{-0.06}$, respectively. Our analysis also provides upper mass limits for the two low-SNR planet candidates in the system; for TOI-125.04 ($R_P=1.36 ~\mathrm{R_{\rm E}}$, $P=$0.53 days) we find a $2\sigma$ upper mass limit of $1.6~\mathrm{M_{\rm E}}$, whereas TOI-125.05 ( $R_P=4.2^{+2.4}_{-1.4} ~\mathrm{R_{\rm E}}$, $P=$ 13.28 days) is unlikely a viable planet candidate with upper mass limit $2.7~\mathrm{M_{\rm E}}$. We discuss the internal structure of the three confirmed planets, as well as dynamical stability and system architecture for this intriguing exoplanet system.

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.

A candidate explanation for Early Mars rivers is atmospheric warming due to surface release of H$_2$ or CH$_4$ gas. However, it remains unknown how much gas could be released in a single event. We model the CH$_4$ release by one mechanism for rapid release of CH$_4$ from clathrate. By modeling how CH$_4$-clathrate release is affected by changes in Mars' obliquity and atmospheric composition, we find that a large fraction of total outgassing from CH$_4$ clathrate occurs following Mars' first prolonged atmospheric collapse. This atmosphere-collapse-initiated CH$_4$-release mechanism has three stages. (1) Rapid collapse of Early Mars' carbon dioxide atmosphere initiates a slower shift of water ice from high ground to the poles. (2) Upon subsequent CO$_2$-atmosphere re-inflation and CO$_2$-greenhouse warming, low-latitude clathrate decomposes and releases methane gas. (3) Methane can then perturb atmospheric chemistry and surface temperature, until photochemical processes destroy the methane. Within our model, we find that under some circumstances a Titan-like haze layer would be expected to form, consistent with transient deposition of abundant complex abiotic organic matter on the Early Mars surface. We also find that this CH$_4$-release mechanism can warm Early Mars, but special circumstances are required in order to uncork 10$^{17}$ kg of CH$_4$, the minimum needed for strong warming. Specifically, strong warming only occurs when the fraction of the hydrate stability zone that is initially occupied by clathrate exceeds 10%, and when Mars' first prolonged atmospheric collapse occurs for atmospheric pressure > 1 bar.

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 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.

The long-standing debate on the existence of ancient oceans on Mars has been recently revived by evidence for tsunami resurfacing events that date from the Late Hesperian geological era. It has been argued that these tsunami events originated from the impact of large meteorites on a deglaciated or nearly deglaciated ocean present in the northern hemisphere of Mars. Here we show that the presence of such a late ocean faces a paradox. If cold, the ocean should have been entirely frozen shortly after its formation, thus preventing the formation of tsunami events. If warm, the ice-free ocean should have produced fluvial erosion of Hesperian Mars terrains much more extensively than previously reported. To solve this apparent paradox, we suggest a list of possible tests and scenarios that could help to reconcile constraints from climate models with tsunami hypothesis. These scenarios could be tested in future dedicated studies.

Transmission spectra of H$_2$O+CO$_2$ mixtures have been recorded, at 296, 325 and 366 K, for various pressures and mixture compositions using two experimental setups. Their analysis enables to retrieve values of the 'continuum' absorption by the CO$_2$-broadened H$_2$O line wings between 100 and 1500 cm$^{-1}$. The results are in good agreement with those, around 1300 cm$^{-1}$, of the single previous experimental study available. Comparisons are also made with direct predictions based on line-shape correction factors $\chi$ calculated, almost thirty years ago, using a quasistatic approach and an input H$_2$O-CO$_2$ intermolecular potential. They show that this model quite nicely predicts, with slightly overestimated values, the continuum over a spectral range where it varies by more than three orders of magnitude. An empirical correction is proposed, based on the experimental data, which should be useful for radiative transfer and climate studies in CO$_2$ rich planetary atmospheres.

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.

We present an experimental study of the absorption, between 40 and 640 cm$^{-1}$, by CO$_2$, CH$_4$ and H$_2$ gases as well as by H$_2$+CO$_2$ and CH$_4$+CO$_2$ mixtures at room temperature. A Fourier transform spectrometer associated to a multi-pass cell, whose optics were adjusted to obtain a 152 m pathlength, were used to record transmission spectra at total pressures up to about 0.98 bar. These measurements provide information concerning the collision-induced absorption (CIA) bands as well as about the wing of the CO$_2$ 15 $\mu$m band. Our results for the CIAs of pure gases are, within uncertainties, in agreement with previous determinations, validating our experimental and data analysis procedures. We then consider the CIAs by H$_2$+CO$_2$ and CH$_4$+CO$_2$ and the low frequency wing of the pure CO$_2$ 15 $\mu$m band, for which there are, to our knowledge, no previous measurements. We confirm experimentally the theoretical prediction of Wordsworth et al. 2017 that the H$_2$+CO$_2$ and CH$_4$+CO$_2$ CIAs are significantly stronger in the 50-550 cm$^{-1}$ region than those of H$_2$+N$_2$ and CH$_4$+N$_2$, respectively. However, we find that the shape and the strength of these recorded CIAs differ from the aforementioned predictions. For the pure CO$_2$ line-wings, we show that both the $\chi$-factor deduced from measurements near 4 $\mu$m and a line-mixing model very well describe the observed strongly sub-Lorentzian behavior in the 500-600 cm$^{-1}$ region. These experimental results open renewed perspectives for studies of the past climate of Mars and extrasolar analogues.

During accretion, the young rocky planets are so hot that they become endowed with a magma ocean. From that moment, the mantle convective thermal flux control the cooling of the planet and an atmosphere is created by outgassing. This atmosphere will then play a key role during this cooling phase. Studying this cooling phase in details is a necessary step to explain the great diversity of the observed telluric planets and especially to assess the presence of surface liquid water. We used here a radiative-convective 1D atmospheric model (H2O, CO2) to study the impact of the Bond albedo on the evolution of magma ocean planets. We derived from this model the thermal emission spectrum and the spectral reflectance of these planets, from which we calculated their Bond albedos. Compared to Marcq et al. (2017), the model now includes a new module to compute the Rayleigh scattering, and state of the art CO2 and H2O gaseous opacities data in the visible and infrared spectral ranges. We show that the Bond albedo of these planets depends on their surface temperature and results from a competition between Rayleigh scattering from the gases and Mie scattering from the clouds. The colder the surface temperature is, the thicker the clouds are, and the higher the Bond albedo is. We also evidence that the relative abundances of CO2 and H2O in the atmosphere strongly impact the Bond albedo. The Bond albedo is higher for atmospheres dominated by the CO2, better Rayleigh scatterer than H2O. Finally, we provide the community with an empirical formula for the Bond albedo that could be useful for future studies of magma ocean planets.

Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting to study the formation and evolution of terrestrial planets that formed in the same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all known to better than 5% precision, their densities have significant uncertainties (between 28% and 95%) because of poor constraints on the planet's masses. Aims.The goal of this paper is to improve our knowledge of the TRAPPIST-1 planetary masses and densities using transit-timing variations (TTV). The complexity of the TTV inversion problem is known to be particularly acute in multi-planetary systems (convergence issues, degeneracies and size of the parameter space), especially for resonant chain systems such as TRAPPIST-1. Methods. To overcome these challenges, we have used a novel method that employs a genetic algorithm coupled to a full N-body integrator that we applied to a set of 284 individual transit timings. This approach enables us to efficiently explore the parameter space and to derive reliable masses and densities from TTVs for all seven planets. Results. Our new masses result in a five- to eight-fold improvement on the planetary density uncertainties, with precisions ranging from 5% to 12%. These updated values provide new insights into the bulk structure of the TRAPPIST-1 planets. We find that TRAPPIST-1\,c and e likely have largely rocky interiors, while planets b, d, f, g, and h require envelopes of volatiles in the form of thick atmospheres, oceans, or ice, in most cases with water mass fractions less than 5%.

In this work, we measured the absorption by CO2+H2O mixtures from 2400 to 2600 cm-1 which corresponds to the spectral region beyond the nu3 band head of CO2. Transmission spectra of CO2 mixed with water vapor were recorded with a high-resolution Fourier-transform spectrometer for various pressure, temperature and concentration conditions. The continuum absorption by CO2 due to the presence of water vapor was determined by subtracting from measured spectra the contribution of local lines of both species, that of the continuum of pure CO2 as well as of the self- and CO2-continua of water vapor induced by the H2O-H2O and H2O-CO2 interactions. The obtained results are in very good agreement with the unique previous measurement (in a narrower spectral range). They confirm that the H2O-continuum of CO2 is significantly larger than that observed for pure CO2. This continuum thus must be taken into account in radiative transfer calculations for media involving CO2+H2O mixture. An empirical model, using sub-Lorentzian line shapes based on some temperature-dependent correction factors chi is proposed which enables an accurate description of the experimental results.

TRAPPIST-1 planets are invaluable for the study of comparative planetary science outside our Solar System and possibly habitability. First, we derive from N-body simulations possible planetary evolution scenarios, and show that each of the planets are likely to be in synchronous rotation. We then use a 3-D Global Climate Model to explore the possible climates of cool planets of the TRAPPIST-1 system. In particular, we look at the conditions required for cool planets to prevent possible volatile species to be lost by permanent condensation, irreversible burying or photochemical destruction. We also explore the resilience of the same volatiles (when in condensed phase) to a runaway greenhouse process. We find that background atmospheres made of N2, CO or O2 are resistant to atmospheric collapse. However, it should be difficult for TRAPPIST-1 planets to accumulate significant greenhouse gases like CO2, CH4, or NH3. CO2 can easily condense on the nightside, forming glaciers that would flow toward the substellar region. A complete CO2 ice cover is possible on TRAPPIST-1g and h only, although CO2 ice deposits could be gravitationally unstable and get buried beneath the water ice shell in geologically short timescales. Given TRAPPIST-1 planets large EUV irradiation (at least 1000x Titan's flux), CH4 and NH3 should be photodissociated rapidly and thus be hard to accumulate in the atmosphere. Photochemical hazes could then sedimentate and form a surface layer of tholins. Regarding habitability, we confirm that few bars of CO2 would suffice to warm the surface of TRAPPIST-1f and g above the melting point of water. We also show that TRAPPIST-1e is a remarkable candidate for surface habitability. If the planet is today synchronous and abundant in water, then it should always sustain surface liquid water at least in the substellar region, whatever the atmosphere considered.

Ozak et al. [2016] claimed that explicitly including the effect of CO$_2$ collisional line mixing (LM) in their radiative transfer calculations yield CO$_2$ atmospheres that are more transparent to infrared radiation than when spectra calculations are made using sub-Lorentzian line shapes. This would in particular imply significantly colder surface temperatures (up to 15K) for early Mars than estimated in previous studies. Here we show that the relative cooling that Ozak et al. [2016] associated to the effect of collisional line mixing is in fact due to a wrong choice of broadening species (air instead of CO$_2$). We then calculated Line-by-Line spectra of pure CO$_2$ atmospheres using a line-mixing model developed for self-broadened CO$_2$. Using the LMD Generic model (in 1D radiative-convective mode), we find that calculations made with the proper collisional line mixing model and with sub-Lorentzian line shapes lead to differences between early Mars surface temperatures smaller than 2 Kelvins only.

It is widely believed that the carbonate-silicate cycle is the main agent to trigger deglaciations by CO$_2$ greenhouse warming on Earth and on Earth-like planets when they get in frozen state. Here we use a 3D Global Climate Model to simulate the ability of frozen planets to escape from glaciation by accumulating enough gaseous CO$_2$. We find that Earth-like planets orbiting a Sun-like star may never be able to escape from glaciation if their orbital distance is greater than $\sim$ 1.27 AU (Flux $<$ 847 W m$^{-2}$ or 62$\%$ of the Solar constant), because CO$_2$ would condense at the poles forming permanent CO$_2$ ice caps. This limits the amount of CO$_2$ in the atmosphere and thus its greenhouse effect. The amount of CO$_2$ that can be trapped in the polar caps depends on the efficiency of CO$_2$ ice to flow laterally as well as its graviational stability relative to subsurface water ice. The flow of CO$_2$ ice from poles to equator is mostly controlled by the bottom temperature, and hence by the internal heat flux. We find that a frozen Earth-like planet located at 1.30 AU of a Sun-like star could store as much as 1.5/4.5/15 bars of dry ice at the poles, for internal heat fluxes of 100/30/10 mW m$^{-2}$. But these amounts are lower limits. For planets with a significant water ice cover, we show that CO$_2$ ice deposits should be gravitationnally unstable. They get buried beneath the water ice cover in geologically short timescales of $\sim$~10$^4$ yrs, mainly controlled by the viscosity of water ice. CO$_2$ would be permanently sequestered underneath the water ice cover, in the form of CO$_2$ liquids, CO$_2$ clathrate hydrates and/or dissolved in subglacial water reservoirs (if any). This would considerably increase the amount of CO$_2$ trapped and further reduce the probability of deglaciation.

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf.

One focus of modern astronomy is to detect temperate terrestrial exoplanets well-suited for atmospheric characterisation. A milestone was recently achieved with the detection of three Earth-sized planets transiting (i.e. passing in front of) a star just 8% the mass of the Sun 12 parsecs away. Indeed, the transiting configuration of these planets with the Jupiter-like size of their host star - named TRAPPIST-1 - makes possible in-depth studies of their atmospheric properties with current and future astronomical facilities. Here we report the results of an intensive photometric monitoring campaign of that star from the ground and with the Spitzer Space Telescope. Our observations reveal that at least seven planets with sizes and masses similar to the Earth revolve around TRAPPIST-1. The six inner planets form a near-resonant chain such that their orbital periods (1.51, 2.42, 4.04, 6.06, 9.21, 12.35 days) are near ratios of small integers. This architecture suggests that the planets formed farther from the star and migrated inward. The seven planets have equilibrium temperatures low enough to make possible liquid water on their surfaces.

Mars was characterized by cataclysmic groundwater-sourced surface flooding that formed large outflow channels and that may have altered the climate for extensive periods during the Hesperian era. In particular, it has been speculated that such events could have induced significant rainfall and caused the formation of late-stage valley networks. We present the results of 3-D Global Climate Model simulations reproducing the short and long term climatic impact of a wide range of outflow channel formation events under cold ancient Mars conditions. We find that the most intense of these events (volumes of water up to 10$^7 $km$^3$ and released at temperatures up to 320 Kelvins) cannot trigger long-term greenhouse global warming, regardless of how favorable are the external conditions (e.g. obliquity and seasons). In any case, outflow channel formation events at any atmospheric pressure are unable to produce rainfall or significant snowmelt at latitudes below 40$^{\circ}$N. On the long term, for an obliquity of $\sim $45$^{\circ}$ and atmospheric pressures $>$ 80 mbar, we find that the lake ice (formed quickly after the outflow event) is transported progressively southward through the mechanisms of sublimation and adiabatic cooling. At the same time, and as long as the initial water reservoir is not entirely sublimated, ice deposits remain in the West Echus Chasma Plateau region where hints of hydrological activity contemporaneous with outflow channel formation events have been observed. However, because the high albedo of ice drives Mars to even colder temperatures, snowmelt produced by seasonal solar forcing is difficult to attain.

Radial velocity monitoring has found the signature of a $M \sin i = 1.3$~M$_\oplus$ planet located within the Habitable Zone (HZ) of Proxima Centauri \citep{Anglada16}. Despite a hotter past and an active host star the planet Proxima~b could have retained enough volatiles to sustain surface habitability \citep{Ribas2016}. Here we use a 3D Global Climate Model (GCM) to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect. We find that a broad range of atmospheric compositions allow surface liquid water. On a tidally-locked planet with sufficient surface water inventory, liquid water is always present, at least in the substellar region. With a non-synchronous rotation, this requires a minimum greenhouse warming ($\sim$10~mbar of CO$_2$ and 1~bar of N$_2$). If the planet is dryer, $\sim$0.5~bar/1.5~bars of CO$_2$ (respectively for asynchronous/synchronous rotation) suffice to prevent the trapping of any arbitrary small water inventory into polar/nightside ice caps. We produce reflection/emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible by direct imaging with forthcoming large telescopes. The angular separation of $7 \lambda/D$ at 1~$\mu$m (with the E-ELT) and a contrast of $\sim$10$^{-7}$ will enable high-resolution spectroscopy and the search for molecular signatures, including H$_2$O, O$_2$, and CO$_2$. The observation of thermal phase curves can be attempted with JWST, thanks to a contrast of $2\times10^{-5}$ at 10~$\mu$m. Proxima~b will also be an exceptional target for future IR interferometers. Within a decade it will be possible to image Proxima~b and possibly determine whether this exoplanet's surface is habitable.

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.