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 away1. Indeed, the transiting configuration of these planets combined 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 facilities1,2,3. 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 inward4,5. The seven planets have equilibrium temperatures low enough to make possible liquid water on their surfaces6,7,8.

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

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.Comment: 42 pages, 8 figures, 2 table

We present the consistent evolution of short-period exoplanets coupling the tidal and gravothermal evolution of the planet. Contrarily to previous similar studies, our calculations are based on the complete tidal evolution equations of the Hut (1981) model, valid at any order in eccentricity, obliquity and spin. We demonstrate both analytically and numerically that except if the system was formed with a nearly circular orbit (e < ∼ 0.2), consistently solving the complete tidal equations is mandatory to derive correct tidal evolution histories. We show that calculations based on tidal models truncated at 2nd order in eccentricity, as done in all previous studies, lead to quantitatively and sometimes even qualitatively erroneous tidal evolutions. As a consequence, tidal energy dissipation rates are severely underestimated in all these calculations and the characteristic timescales for the various orbital parameters evolutions can be wrong by up to three orders of magnitude. These discrepancies can by no means be justified by invoking the uncertainty in the tidal quality factors. Based on these complete, consistent calculations, we revisit the viability of the tidal heating hypothesis to explain the anomalously large radius of transiting giant planets. We show that even though tidal dissipation does provide a substantial contribution to the planet's heat budget and can explain some of the moderately bloated hot-Jupiters, this mechanism can not explain alone the properties of the most inflated objects, including HD 209 458 b. Indeed, solving the complete tidal equations shows that enhanced tidal dissipation and thus orbit circularization occur too early during the planet's evolution to provide enough extra energy at the present epoch. In that case either a third, so far undetected, low-mass companion must be present to keep exciting the eccentricity of the giant planet, or other mechanisms -stellar irradiation induced surface winds dissipating in the planet's tidal bulges and thus reaching the convective layers, inefficient flux transport by convection in the planet's interior -must be invoked, together with tidal dissipation, to provide all the pieces of the abnormally large exoplanet puzzle.

Because the solar luminosity increases over geological timescales, Earth climate is expected to warm, increasing water evaporation which, in turn, enhances the atmospheric greenhouse effect. Above a certain critical insolation, this destabilizing greenhouse feedback can "runaway" until all the oceans are evaporated 1,2,3,4 . Through increases in stratospheric humidity, warming may also cause oceans to escape to space before the runaway greenhouse occurs 5,6 . The critical insolation thresholds for these processes, however, remain uncertain because they have so far been evaluated with unidimensional models that cannot account for the dynamical and cloud feedback effects that are key stabilizing features of Earth's climate. Here we use a 3D global climate model to show that the threshold for the runaway greenhouse is about 375 W/m 2 , significantly higher than previously thought 6,7 . Our model is specifically developed to quantify the climate response of Earth-like planets to increased insolation in hot and extremely moist atmospheres. In contrast with previous studies, we find that clouds have a destabilizing feedback on the long term warming. However, subsident, unsaturated regions created by the Hadley circulation have a stabilizing effect that is strong enough to defer the runaway greenhouse limit to higher insolation than inferred from 1D models. Furthermore, because of wavelength-dependent radiative effects, the stratosphere remains cold and dry enough to hamper atmospheric water escape, even at large fluxes. This has strong implications for Venus early water history and extends the size of the habitable zone around other stars. Planetary atmospheres naturally settle into a thermal equilibrium state where their outgoing thermal emission balances the heating due to sunlight absorption. The resulting climate is stabilized by the fact that a temperature increase results in an enhanced thermal-emission cooling. When a condensable greenhouse gas is present at the surface, such as water on Earth, this stabilizing feedback is somewhat hampered by the destabilizing greenhouse feedback: evaporation, and thus water vapor greenhouse effect, increases with temperature, reducing the cooling. Fortunately, under present Earth conditions, this greenhouse feedback is both strong enough to maintain clement surface temperatures and weak enough for the climate to remain stable.When solar heating becomes stronger, however, water vapor can become abundant enough to make the atmosphere optically thick at all thermal wavelengths 7,8 . Thermal flux then originates from the upper troposphere only and reaches a maximum, ∼ 282 W/m 2 , independently from the surface temperature 9 . If the planet absorbs more than this critical flux, thermal equilibrium can be restored only by vaporizing all the water available and reaching high surface temperatures at which the surface starts to radiate at visible wavelengths 4,7 . This is the runaway greenhouse state.Because it has mostly been studied through unidimensional atmosphere models, the afor...

On the basis of geological evidence, it is often stated that the early martian climate was warm enough for liquid water to flow on the surface thanks to the greenhouse effect of a thick atmosphere. We present 3D global climate simulations of the early martian climate performed assuming a faint young sun and a CO 2 atmosphere with surface pressure between 0.1 and 7 bars. The model includes a detailed representation of radiative transfer using revised CO 2 gas collision induced absorption properties, and a parameterisation of CO2 ice cloud microphysical and radiative properties. A wide range of possible climates is explored using various values of obliquities, orbital parameters, cloud microphysic parameters, atmospheric dust loading, and surface properties.Unlike on present-day Mars, for pressures higher than a fraction of a bar surface temperatures vary with altitude because of adiabatic cooling / warming of the atmosphere. In most simulations, CO 2 ice clouds cover a major part of the planet. Previous studies suggested that they could have warmed the planet thanks to their scattering greenhouse effect. However, even assuming parameters that maximize this effect, it does not exceed +15 K. Combined with the revised CO 2 spectroscopy and the impact of surface CO 2 ice on the planetary albedo, we find that a CO 2 atmosphere could not have raised the annual mean temperature above 0 • C anywhere on the planet. The collapse of the atmosphere into permanent CO 2 ice caps is predicted for pressures higher than 3 bar, or conversely at pressure lower than one bar if the obliquity is low enough. Summertime diurnal mean surface temperatures above 0 • C (a condition which could have allowed rivers and lakes to form) are predicted for obliquity larger than 40 • at high latitudes but not in locations where most valley networks or layered sedimentary units are observed. In the absence of other warming mechanisms, our climate model results are thus consistent with a cold early Mars scenario in which non climatic mechanisms must occur to explain the evidence for liquid water. In a companion paper by Wordsworth et al., we simulate the hydrological cycle on such a planet and discuss how this could have happened in more detail.

The inner edge of the classical habitable zone is often defined by the critical flux needed to trigger the runaway greenhouse instability. This 1D notion of a critical flux, however, may not be all that relevant for inhomogeneously irradiated planets, or when the water content is limited (land planets). Based on results from our 3D global climate model, we present general features of the climate and large-scale circulation on close-in terrestrial planets. We find that the circulation pattern can shift from super-rotation to stellar/anti stellar circulation when the equatorial Rossby deformation radius significantly exceeds the planetary radius, changing the redistribution properties of the atmosphere. Using analytical and numerical arguments, we also demonstrate the presence of systematic biases among mean surface temperatures and among temperature profiles predicted from either 1D or 3D simulations. After including a complete modeling of the water cycle, we further demonstrate that two stable climate regimes can exist for land planets closer than the inner edge of the classical habitable zone. One is the classical runaway state where all the water is vaporized, and the other is a collapsed state where water is captured in permanent cold traps. We identify this "moist" bistability as the result of a competition between the greenhouse effect of water vapor and its condensation on the night side or near the poles, highlighting the dynamical nature of the runaway greenhouse effect. We also present synthetic spectra showing the observable signature of these two states. Taking the example of two prototype planets in this regime, namely Gl 581 c and HD 85512 b, we argue that depending on the rate of water delivery and atmospheric escape during the life of these planets, they could accumulate a significant amount of water ice at their surface. If such a thick ice cap is present, various physical mechanisms observed on Earth (e.g., gravity driven ice flows, geothermal flux) should come into play to produce long-lived liquid water at the edge and/or bottom of the ice cap. Consequently, the habitability of planets at smaller orbital distance than the inner edge of the classical habitable zone cannot be ruled out. Transiting planets in this regime represent promising targets for upcoming exoplanet characterization observatories, such as EChO and JWST.

While conventional interior models for Jupiter and Saturn are based on the simplistic assumption of a solid core surrounded by a homogeneous gaseous envelope, we have derived new models with an inhomogeneous distribution of heavy elements within these planets. Such a compositional gradient hampers large-scale convection that turns into double-diffusive convection, yielding an inner thermal profile that departs from the traditionally assumed adiabatic interior and affecting these planets heat content and cooling history.To address this problem, we have developed an analytical approach to describe layered double-diffusive convection and apply this formalism to solar system gaseous giant planet interiors. These models satisfy all observational constraints and yield values for the metal enrichment of our gaseous giants that are up to 30% to 60% higher than previously thought. The models also constrain the size of the convective layers within the planets. Because the heavy elements tend to be redistributed within the gaseous envelope, the models predict smaller than usual central cores inside Saturn and Jupiter, with possibly no core for the latter. These models open a new window and raise new challenges to our understanding of the internal structure of giant (solar and extrasolar) planets, in particular on how to determine their heavy material content, a key diagnostic for planet formation theories.