We present the analysis of the hot-Jupiter KELT-7 b using transmission and emission spectroscopy from the Hubble Space Telescope (HST), both taken with the Wide Field Camera 3 (WFC3). Our study uncovers a rich transmission spectrum which is consistent with a cloud-free atmosphere and suggests the presence of H2O and H-. In contrast, the extracted emission spectrum does not contain strong absorption features and, although it is not consistent with a simple blackbody, it can be explained by a varying temperature-pressure profile, collision induced absorption (CIA) and H-. KELT-7 b had also been studied with other space-based instruments and we explore the effects of introducing these additional datasets. Further observations with Hubble, or the next generation of space-based telescopes, are needed to allow for the optical opacity source in transmission to be confirmed and for molecular features to be disentangled in emission.

Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

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.

The history of the Earth has been marked by major ecological transitions, driven by metabolic innovation, that radically reshaped the composition of the oceans and atmosphere. The nature and magnitude of the earliest transitions, hundreds of million years before photosynthesis evolved, remain poorly understood. Using a novel ecosystem-planetary model, we find that pre-photosynthetic methane-cycling microbial ecosystems are much less productive than previously thought. In spite of their low productivity, the evolution of methanogenic metabolisms strongly modifies the atmospheric composition, leading to a warmer but less resilient climate. As the abiotic carbon cycle responds, further metabolic evolution (anaerobic methanotrophy) may feed back to the atmosphere and destabilize the climate, triggering a transient global glaciation. Although early metabolic evolution may cause strong climatic instability, a low CO:CH4 atmospheric ratio emerges as a robust signature of simple methane-cycling ecosystems on a globally reduced planet such as the late Hadean/early Archean Earth.

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

This paper presents the atmospheric characterisation of three large, gaseous planets: WASP-127 b, WASP-79 b and WASP-62 b. We analysed spectroscopic data obtained with the G141 grism (1.088 - 1.68 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) using the Iraclis pipeline and the TauREx3 retrieval code, both of which are publicly available. For WASP-127 b, which is the least dense planet discovered so far and is located in the short-period Neptune desert, our retrieval results found strong water absorption corresponding to an abundance of log(H$_2$O) = -2.71$^{+0.78}_{-1.05}$, and absorption compatible with an iron hydride abundance of log(FeH)=$-5.25^{+0.88}_{-1.10}$, with an extended cloudy atmosphere. We also detected water vapour in the atmospheres of WASP-79 b and WASP-62 b, with best-fit models indicating the presence of iron hydride, too. We used the Atmospheric Detectability Index (ADI) as well as Bayesian log evidence to quantify the strength of the detection and compared our results to the hot Jupiter population study by Tsiaras et al. 2018. While all the planets studied here are suitable targets for characterisation with upcoming facilities such as the James Webb Space Telescope (JWST) and Ariel, WASP-127 b is of particular interest due to its low density, and a thorough atmospheric study would develop our understanding of planet formation and migration.

We analyse the transmission and emission spectra of the ultra-hot Jupiter WASP-76 b, observed with the G141 grism of the Hubble Space Telescope's Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauRex 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-grey clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyse this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-grey cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterisation with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers.

Recently, gas disks have been discovered around main sequence stars well beyond the usual protoplanetary disk lifetimes (i.e., > 10 Myrs), when planets have already formed. These gas disks, mainly composed of CO, carbon, and oxygen seem to be ubiquitous in systems with planetesimal belts (similar to our Kuiper belt), and can last for hundreds of millions of years. Planets orbiting in these gas disks will accrete a large quantity of gas that will transform their primordial atmospheres into new secondary atmospheres with compositions similar to that of the parent gas disk. Here, we quantify how large a secondary atmosphere can be created for a variety of observed gas disks and for a wide range of planet types. We find that gas accretion in this late phase is very significant and an Earth's atmospheric mass of gas is readily accreted on terrestrial planets in very tenuous gas disks. In slightly more massive disks, we show that massive CO atmospheres can be accreted, forming planets with up to sub-Neptune-like pressures. Our new results demonstrate that new secondary atmospheres with high metallicities and high C/O ratios will be created in these late gas disks, resetting their primordial compositions inherited from the protoplanetary disk phase, and providing a new birth to planets that lost their atmosphere to photoevaporation or giant impacts. We therefore propose a new paradigm for the formation of atmospheres on low-mass planets, which can be tested with future observations (JWST, ELT, ARIEL). We also show that this late accretion would show a very clear signature in Sub-Neptunes or cold exo-Jupiters. Finally, we find that accretion creates cavities in late gas disks, which could be used as a new planet detection method, for low mass planets a few au to a few tens of au from their host stars.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect a planetary mass companion, PDS70b, within the disk cavity. Moreover, observations in H_alpha with MagAO and MUSE revealed emission associated to PDS70b and to another new companion candidate, PDS70c, at a larger separation from the star. Aims. Our aim is to confirm the discovery of the second planet PDS70c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g. ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS70c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 M Jup and with a T_eff around 900 K. Also, it has a low gravity with log g between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12") from the star and gravitationally bound to the star. Its spectrum is however very blue, so that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We, however, cannot completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of few tens of M_Earth. Moreover, we propose a possible structure for the planetary system based on our data that, however, cannot be stable on a long timescale.

In this white paper, we present a cross-section of important scientific questions that remain partially or completely unanswered, ranging from Titan exosphere to the deep interior, and we detail which instrumentation and mission scenarios should be used to answer them. Our intention is to formulate the science goals for the next generation of planetary missions to Titan in order to prepare the future exploration of the moon. The ESA L-class mission concept that we propose is composed of a Titan orbiter and at least an in situ element (lake lander and/or drone(s)).

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

Three dimensional models that account for chemistry are useful tools to predict the chemical composition of (exo)planet and brown dwarf atmospheres and interpret observations of future telescopes. Recent Juno observations of the NH3 tropospheric distribution in Jupiter also indicate that 3D chemical modelling may be necessary to constrain the deep composition of the giant planets of the Solar System. However, due to the high computational cost of chemistry calculations, 3D chemical modelling has so far been limited. A solution to include chemical kinetics in a 3D model is to use a reduced chemical scheme. In this view, we have developed a reduced scheme containing 30 species and 181 reversible reactions, from the full chemical scheme of Venot et al. 2012. This scheme reproduces accurately the vertical profiles of the observable species (H2O, CH4, CO, CO2, NH3, and HCN) and has a large range of validity. It can be used to study all kind of warm atmospheres, except hot C-rich ones, which contains a high amount of C2H2. It can be used in 1D models, for fast computations, but also in 3D models for hot giant (exo)planet and brown dwarf atmospheres.

Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus around the planet. Consequently, studies have focused on the effect of atmospheric horizontal heterogeneities along the limb. Here we demonstrate that the region probed in transmission -- the limb -- extends significantly toward the day and night sides of the planet. We show that the strong day-night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb to result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures. We first produce synthetic JWST observations from a simulation of GJ 1214b and demonstrate the necessity of a real 3D approach to model data for such precise observatories. Second, we investigate how day-night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesize a large set of forward spectra for prototypical HD 209458 b and GJ 1214 b type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrieval analyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside. These biases are difficult to detect as the 1D profiles used in the retrieval are found to provide an excellent match to the observed spectra. This needs to be kept in mind when interpreting real data.

HIP 65426 b is a recently discovered exoplanet imaged during the course of the SPHERE-SHINE survey. Here we present new $L'$ and $M'$ observations of the planet from the NACO instrument at the VLT from the NACO-ISPY survey, as well as a new $Y-H$ spectrum and $K$-band photometry from SPHERE-SHINE. Using these data, we confirm the nature of the companion as a warm, dusty planet with a mid-L spectral type. From comparison of its SED with the BT-Settl atmospheric models, we derive a best-fit effective temperature of $T_{\text{eff}}=1618\pm7$ K, surface gravity $\log g=3.78^{+0.04}_{-0.03}$ and radius $R=1.17\pm0.04$ $R_{\text{J}}$ (statistical uncertainties only). Using the DUSTY and COND isochrones we estimate a mass of $8\pm1$ $M_{\text{J}}$. Combining the astrometric measurements from our new datasets and from the literature, we show the first indications of orbital motion of the companion (2.6$\sigma$ significance) and derive preliminary orbital constraints. We find a highly inclined orbit ($i=107^{+13}_{-10}$ deg) with an orbital period of $800^{+1200}_{-400}$ yr. We also report SPHERE sparse aperture masking observations that investigate the possibility that HIP 65426 b was scattered onto its current orbit by an additional companion at a smaller orbital separation. From this data we rule out the presence of brown dwarf companions with masses greater than 16 $M_{\text{J}}$ at separations larger than 3 AU, significantly narrowing the parameter space for such a companion.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

We report on results of an observing campaign to support the Juno mission. At the beginning of 2016, using TEXES (Texas Echelon cross-dispersed Echelle Spectrograph), mounted on the NASA Infrared Telescope Facility (IRTF), we obtained data cubes of Jupiter in the 1930--1943 cm$^{-1}$ spectral ranges (around 5 $\mu$m), which probe the atmosphere in the 1--4 bar region, with a spectral resolution of $\approx$ 0.15 cm$^{-1}$ and an angular resolution of $\approx$ 1.4". This dataset is analysed by a code that combines a line-by-line radiative transfer model with a non-linear optimal estimation inversion method. The inversion retrieves the vertical abundance profiles of NH$_3$ - which is the main contributor at these wavelengths - with a maximum sensitivity at $\approx$ 1--3 bar, as well as the cloud transmittance. This retrieval is performed on over more than one thousand pixels of our data cubes, producing maps of the disk, where all the major belts are visible. We present our retrieved NH$_3$ abundance maps which can be compared with the distribution observed by Juno's MWR (Bolton et al., 2017; Li et al., 2017) in the 2 bar region and discuss their significance for the understanding of Jupiter's atmospheric dynamics. We are able to show important latitudinal variations - such as in the North Equatorial Belt (NEB), where the NH$_3$ abundance is observed to drop down to 60 ppmv at 2 bar - as well as longitudinal variability. In the zones, we find the NH$_3$ abundance to increase with depth, from 100 $\pm$ 15 ppmv at 1 bar to 500 $\pm$ 30 ppmv at 3 bar. We also display the cloud transmittance--NH$_3$ abundance relationship, and find different behaviour for the NEB, the other belts and the zones. Using a simple cloud model, we are able to fit this relationship, at least in the NEB, including either NH$_3$-ice or NH$_4$SH particles with sizes between 10 and 100 $\mu$m.

Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball states within the habitable zone. The large amplitude changes in obliquity and eccentricity cause the ice edge, the lowest latitude extent of the ice caps, to become unstable and grow to the equator. We apply an analytical theory of the ice edge latitude to show that obliquity is the primary driver of the instability. The thermal inertia of the ice sheets and the spectral energy distribution of the G dwarf star increase the sensitivity of the model to triggering runaway glaciation. Finally, we apply a machine learning algorithm to demonstrate how this technique can be used to extend the power of climate models. This work illustrates the importance of orbital evolution for habitability in dynamically rich planetary systems. We emphasize that as potentially habitable planets are discovered around G dwarfs, we need to consider orbital dynamics.

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.

The obliquity of the Earth, which controls our seasons, varies by only ~2.5 degrees over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth's ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semi-analytic model that allows us to map broad regions of parameter space. We find that in general, the largest drivers of obliquity variations are secular spin-orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.

HR2562 is an F5V star located at ~33 pc from the Sun hosting a substellar companion that was discovered using the GPI instrument. The main objective of the present paper is to provide an extensive characterisation of the substellar companion, by deriving its fundamental properties. We observed HR 2562 with the near-infrared branch (IFS and IRDIS) of SPHERE at the VLT. During our observations IFS was operating in the YJ band, while IRDIS was observing with the $H$ broad-band filter. The data were reduced with the dedicated SPHERE GTO pipeline, which is custom-designed for this instrument. On the reduced images, we then applied the post-processing procedures that are specifically prepared to subtract the speckle noise. The companion is clearly detected in both IRDIS and IFS datasets. We obtained photometry in three different spectral bands. The comparison with template spectra allowed us to derive a spectral type of T2-T3 for the companion. Using both evolutionary and atmospheric models we inferred the main physical parameters of the companion obtaining a mass of 32+/-14$~MJup, T_{eff}=1100+/-200 K and log(g)=4.75+/-0.41.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

We present simultaneous HST WFC3 + Spitzer IRAC variability monitoring for the highly-variable young ($\sim$20 Myr) planetary-mass object PSO J318.5-22. Our simultaneous HST + Spitzer observations covered $\sim$2 rotation periods with Spitzer and most of a rotation period with HST. We derive a period of 8.6$\pm$0.1 hours from the Spitzer lightcurve. Combining this period with the measured $v sin i$ for this object, we find an inclination of 56.2$\pm 8.1^{\circ}$. We measure peak-to-trough variability amplitudes of 3.4$\pm$0.1$\%$ for Spitzer Channel 2 and 4.4 - 5.8$\%$ (typical 68$\%$ confidence errors of $\sim$0.3$\%$) in the near-IR bands (1.07-1.67 $\mu$m) covered by the WFC3 G141 prism -- the mid-IR variability amplitude for PSO J318.5-22 one of the highest variability amplitudes measured in the mid-IR for any brown dwarf or planetary mass object. Additionally, we detect phase offsets ranging from 200--210$^{\circ}$ (typical error of $\sim$4$^{\circ}$) between synthesized near-IR lightcurves and the Spitzer mid-IR lightcurve, likely indicating depth-dependent longitudinal atmospheric structure in this atmosphere. The detection of similar variability amplitudes in wide spectral bands relative to absorption features suggests that the driver of the variability may be inhomogeneous clouds (perhaps a patchy haze layer over thick clouds), as opposed to hot spots or compositional inhomogeneities at the top-of-atmosphere level.

We developed a simple, physical and self-consistent cloud model for brown dwarfs and young giant exoplanets. We compared different parametrisations for the cloud particle size, by either fixing particle radii, or fixing the mixing efficiency (parameter fsed) or estimating particle radii from simple microphysics. The cloud scheme with simple microphysics appears as the best parametrisation by successfully reproducing the observed photometry and spectra of brown dwarfs and young giant exoplanets. In particular, it reproduces the L-T transition, due to the condensation of silicate and iron clouds below the visible/near-IR photosphere. It also reproduces the reddening observed for low-gravity objects, due to an increase of cloud optical depth for low gravity. In addition, we found that the cloud greenhouse effect shifts chemical equilibriums, increasing the abundances of species stable at high temperature. This effect should significantly contribute to the strong variation of methane abundance at the L-T transition and to the methane depletion observed on young exoplanets. Finally, we predict the existence of a continuum of brown dwarfs and exoplanets for absolute J magnitude=15-18 and J-K color=0-3, due to the evolution of the L-T transition with gravity. This self-consistent model therefore provides a general framework to understand the effects of clouds and appears well-suited for atmospheric retrievals.

The substellar companion HD206893b has recently been discovered by direct imaging of its disc-bearing host star with the SPHERE instrument. We investigate the atypical properties of the companion, which has the reddest near-infrared colours among all known substellar objects, either orbiting a star or isolated, and we provide a comprehensive characterisation of the host star-disc-companion system. We conducted a follow-up of the companion with adaptive optics imaging and spectro-imaging with SPHERE, and a multiinstrument follow-up of its host star. We obtain a R=30 spectrum from 0.95 to 1.64 micron of the companion and additional photometry at 2.11 and 2.25 micron. We carried out extensive atmosphere model fitting for the companions and the host star in order to derive their age, mass, and metallicity. We found no additional companion in the system in spite of exquisite observing conditions resulting in sensitivity to 6MJup (2MJup) at 0.5" for an age of 300 Myr (50 Myr). We detect orbital motion over more than one year and characterise the possible Keplerian orbits. We constrain the age of the system to a minimum of 50 Myr and a maximum of 700 Myr, and determine that the host-star metallicity is nearly solar. The comparison of the companion spectrum and photometry to model atmospheres indicates that the companion is an extremely dusty late L dwarf, with an intermediate gravity (log g 4.5-5.0) which is compatible with the independent age estimate of the system. Though our best fit corresponds to a brown dwarf of 15-30 MJup aged 100-300 Myr, our analysis is also compatible with a range of masses and ages going from a 50 Myr 12MJup planetary-mass object to a 50 MJup Hyades-age brown dwarf...

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.

Oxygen isotopes in marine cherts have been used to infer hot oceans during the Archean with temperatures between 60{\deg}C (333 K) and 80{\deg}C (353 K). Such climates are challenging for the early Earth warmed by the faint young Sun. The interpretation of the data has therefore been controversial. 1D climate modeling inferred that such hot climates would require very high levels of CO2 (2-6 bars). Previous carbon cycle modeling concluded that such stable hot climates were impossible and that the carbon cycle should lead to cold climates during the Hadean and the Archean. Here, we revisit the climate and carbon cycle of the early Earth at 3.8 Ga using a 3D climate-carbon model. We find that CO2 partial pressures of around 1 bar could have produced hot climates given a low land fraction and cloud feedback effects. However, such high CO2 partial pressures should not have been stable because of the weathering of terrestrial and oceanic basalts, producing an efficient stabilizing feedback. Moreover, the weathering of impact ejecta during the Late Heavy Bombardment (LHB) would have strongly reduced the CO2 partial pressure leading to cold climates and potentially snowball Earth events after large impacts. Our results therefore favor cold or temperate climates with global mean temperatures between around 8{\deg}C (281 K) and 30{\deg}C (303 K) and with 0.1-0.36 bar of CO2 for the late Hadean and early Archean. Finally, our model suggests that the carbon cycle was efficient for preserving clement conditions on the early Earth without necessarily requiring any other greenhouse gas or warming process.

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.

As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is $0.00 \leqslant e \leqslant 0.32$, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO$_2$ in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60$^\circ$ and 90$^\circ$). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO$_2$ is increased to 5-bar levels. In a low-CO$_2$ case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO$_2$ levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet's substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.

The next generation of ground and space-based telescopes will image habitable planets around nearby stars. A growing literature describes how to characterize such planets with spectroscopy, but less consideration has been given to the usefulness of planet colors. Here, we investigate whether potentially Earth-like exoplanets could be identified using UV-visible-to-NIR wavelength broadband photometry (350-1000 nm). Specifically, we calculate optimal photometric bins for identifying an exo-Earth and distinguishing it from uninhabitable planets including both Solar System objects and model exoplanets. The color of some hypothetical exoplanets - particularly icy terrestrial worlds with thick atmospheres - is similar to Earth's because of Rayleigh scattering in the blue region of the spectrum. Nevertheless, subtle features in Earth's reflectance spectrum appear to be unique. In particular, Earth's reflectance spectrum has a 'U-shape' unlike all our hypothetical, uninhabitable planets. This shape is partly biogenic because O2-rich, oxidizing air is transparent to sunlight, allowing prominent Rayleigh scattering, while ozone absorbs visible light, creating the bottom of the 'U'. Whether such uniqueness has practical utility depends on observational noise. If observations are photon limited or dominated by astrophysical sources (zodiacal light or imperfect starlight suppression), then the use of broadband visible wavelength photometry to identify Earth twins has little practical advantage over obtaining detailed spectra. However, if observations are dominated by dark current then optimized photometry could greatly assist preliminary characterization. We also calculate the optimal photometric bins for identifying extrasolar Archean Earths, and find that the Archean Earth is more difficult to unambiguously identify than a modern Earth twin.

The warm sub-Neptune GJ1214b has a featureless transit spectrum which may be due to the presence of high and thick clouds or haze. Here, we simulate the atmosphere of GJ1214b with a 3D General Circulation Model for cloudy hydrogen-dominated atmospheres, including cloud radiative effects. We show that the atmospheric circulation is strong enough to transport micrometric cloud particles to the upper atmosphere and generally leads to a minimum of cloud at the equator. By scattering stellar light, clouds increase the planetary albedo to 0.4-0.6 and cool the atmosphere below 1 mbar. However, the heating by ZnS clouds leads to the formation of a stratospheric thermal inversion above 10 mbar, with temperatures potentially high enough on the dayside to evaporate KCl clouds. We show that flat transit spectra consistent with HST observations are possible if cloud particle radii are around 0.5 micron, and that such clouds should be optically thin at wavelengths > 3 microns. Using simulated cloudy atmospheres that fit the observed spectra we generate transit, emission and reflection spectra and phase curves for GJ1214b. We show that a stratospheric thermal inversion would be readily accessible in near and mid-infrared atmospheric spectral windows. We find that the amplitude of the thermal phase curves is strongly dependent on metallicity, but only slightly impacted by clouds. Our results suggest that primary and secondary eclipses and phase curves observed by the James Webb Space Telescope in the near to mid-infrared should provide strong constraints on the nature of GJ1214b's atmosphere and clouds.

GJ1214b is a warm sub-Neptune transiting in front of a nearby M dwarf star. Recent observations indicate the presence of high and thick clouds or haze whose presence requires strong atmospheric mixing. In order to understand the transport and distribution of such clouds/haze, we study the atmospheric circulation and the vertical mixing of GJ1214b with a 3D General Circulation Model for cloud-free hydrogen-dominated atmospheres (metallicity of 1, 10 and 100 times the solar value) and for a water-dominated atmosphere. We analyze the effect of the atmospheric metallicity on the thermal structure and zonal winds. We also analyze the zonal mean meridional circulation and show that it corresponds to an anti-Hadley circulation in most of the atmosphere with upwelling at mid-latitude and downwelling at the equator in average. This circulation must be present on a large range of synchronously rotating exoplanets with strong impact on cloud formation and distribution. Using simple tracers, we show that vertical winds on GJ1214b can be strong enough to loft micrometric particles and that the anti-Hadley circulation leads to a minimum of tracers at the equator. We find that the strength of the vertical mixing increases with metallicity. We derive 1D equivalent eddy diffusion coefficients and find simple parametrizations from Kzz=7x10^2xP_{bar}^{-0.4} m^2/s for solar metallicity to Kzz=3x10^3xP_{bar}^{-0.4} m^2/s for the 100xsolar metallicity. These values should favor an efficient formation of photochemical haze in the upper atmosphere of GJ1214b.

Titan's equatorial regions are covered by eastward propagating linear dunes. This direction is opposite to mean surface winds simulated by Global Climate Models (GCMs), which are oriented westward at these latitudes, similar to trade winds on Earth. Different hypotheses have been proposed to address this apparent contradiction, involving Saturn's gravitational tides, large scale topography or wind statistics, but none of them can explain a global eastward dune propagation in the equatorial band. Here we analyse the impact of equinoctial tropical methane storms developing in the superrotating atmosphere (i.e. the eastward winds at high altitude) on Titan's dune orientation. Using mesoscale simulations of convective methane clouds with a GCM wind profile featuring superrotation, we show that Titan's storms should produce fast eastward gust fronts above the surface. Such gusts dominate the aeolian transport, allowing dunes to extend eastward. This analysis therefore suggests a coupling between superrotation, tropical methane storms and dune formation on Titan. Furthermore, together with GCM predictions and analogies to some terrestrial dune fields, this work provides a general framework explaining several major features of Titan's dunes: linear shape, eastward propagation and poleward divergence, and implies an equatorial origin of Titan's dune sand.

Several clues indicate that Titan's atmosphere has been depleted in methane during some period of its history, possibly as recently as 0.5-1 billion years ago. It could also happen in the future. Under these conditions, the atmosphere becomes only composed of nitrogen with a range of temperature and pressure allowing liquid or solid nitrogen to condense. Here, we explore these exotic climates throughout Titan's history with a 3D Global Climate Model (GCM) including the nitrogen cycle and the radiative effect of nitrogen clouds. We show that for the last billion years, only small polar nitrogen lakes should have formed. Yet, before 1 Ga, a significant part of the atmosphere could have condensed, forming deep nitrogen polar seas, which could have flowed and flooded the equatorial regions. Alternatively, nitrogen could be frozen on the surface like on Triton, but this would require an initial surface albedo higher than 0.65 at 4 Ga. Such a state could be stable even today if nitrogen ice albedo is higher than this value. According to our model, nitrogen flows and rain may have been efficient to erode the surface. Thus, we can speculate that a paleo-nitrogen cycle may explain the erosion and the age of Titan's surface, and may have produced some of the present valley networks and shorelines. Moreover, by diffusion of liquid nitrogen in the crust, a paleo-nitrogen cycle could be responsible of the flattening of the polar regions and be at the origin of the methane outgassing on Titan.

The analyses of Titan's gravity field obtained by Cassini space mission suggest the presence of an internal ocean beneath its icy surface. The characterization of the geophysical parameters of the icy shell and the ocean is important to constrain the evolution models of Titan. The knowledge of the librations, that are periodic oscillations around a uniform rotational motion, can bring piece of information on the interior parameters. The objective of this paper is to study the librational response in longitude from an analytical approach for Titan composed of a deep atmosphere, an elastic icy shell, an internal ocean, and an elastic rocky core perturbed by the gravitational interactions with Saturn. We start from the librational equations developed for a rigid satellite in synchronous spin-orbit resonance. We introduce explicitly the atmospheric torque acting on the surface computed from the Titan IPSL GCM (Institut Pierre Simon Laplace General Circulation Model) and the periodic deformations of elastic solid layers due to the tides. We investigate the librational response for various interior models in order to compare and to identify the influence of the geophysical parameters and the impact of the elasticity. The main librations arise at two well-separated forcing frequency ranges: low forcing frequencies dominated by the Saturnian annual and semi-annual frequencies, and a high forcing frequency regime dominated by Titan's orbital frequency around Saturn. We find that internal structure models including an internal ocean with elastic solid layers lead to the same order of libration amplitude than the oceanless models, which makes more challenging to differentiate them by the interpretation of librational motion.

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. Through increases in stratospheric humidity, warming may also cause oceans to escape to space before the runaway greenhouse occurs. 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. 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.

Different solutions have been proposed to solve the "faint young Sun problem", defined by the fact that the Earth was not fully frozen during the Archean despite the fainter Sun. Most previous studies were performed with simple 1-D radiative convective models and did not account well for the clouds and ice-albedo feedback or the atmospheric and oceanic transport of energy. We apply a global climate model (GCM) to test the different solutions to the faint young Sun problem. We explore the effect of greenhouse gases (CO2 and CH4), atmospheric pressure, cloud droplet size, land distribution, and Earth's rotation rate. We show that neglecting organic haze, 100 mbar of CO2 with 2 mbar of CH4 at 3.8 Ga and 10 mbar of CO2 with 2 mbar of CH4 at 2.5 Ga allow a temperate climate (mean surface temperature between 10{\deg}C and 20{\deg}C). Such amounts of greenhouse gases remain consistent with the geological data. Removing continents produces a warming lower than +4{\deg}C. The effect of rotation rate is even more limited. Larger droplets (radii of 17 microns versus 12 microns) and a doubling of the atmospheric pressure produce a similar warming of around +7{\deg}C. In our model, ice-free water belts can be maintained up to 25{\deg}N/S with less than 1 mbar of CO2 and no methane. An interesting cloud feedback appears above cold oceans, stopping the glaciation. Such a resistance against full glaciation tends to strongly mitigate the faint young Sun problem.

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 so relevant for inhomogeneously irradiated planets, or when the water content is limited (land planets). Here, based on results from our 3D global climate model, 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. Using analytical and numerical arguments, we also demonstrate the presence of systematic biases between mean surface temperatures or temperature profiles predicted from either 1D or 3D simulations. Including a complete modeling of the water cycle, we further demonstrate that for land planets closer than the inner edge of the classical habitable zone, two stable climate regimes can exist. One is the classical runaway state, 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. We also present synthetic spectra showing the observable signature of these two states. Taking the example of two prototype planets in this regime, namely Gl581c and HD85512b, we argue that they could accumulate a significant amount of water ice at their surface. If such a thick ice cap is present, gravity driven ice flows and 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 observatories like EChO and JWST.

We discuss 3D global simulations of the early Martian climate that we have performed assuming a faint young Sun and denser CO2 atmosphere. We include a self-consistent representation of the water cycle, with atmosphere-surface interactions, atmospheric transport, and the radiative effects of CO2 and H2O gas and clouds taken into account. We find that for atmospheric pressures greater than a fraction of a bar, the adiabatic cooling effect causes temperatures in the southern highland valley network regions to fall significantly below the global average. Long-term climate evolution simulations indicate that in these circumstances, water ice is transported to the highlands from low-lying regions for a wide range of orbital obliquities, regardless of the extent of the Tharsis bulge. In addition, an extended water ice cap forms on the southern pole, approximately corresponding to the location of the Noachian/Hesperian era Dorsa Argentea Formation. Even for a multiple-bar CO2 atmosphere, conditions are too cold to allow long-term surface liquid water. Limited melting occurs on warm summer days in some locations, but only for surface albedo and thermal inertia conditions that may be unrealistic for water ice. Nonetheless, meteorite impacts and volcanism could potentially cause intense episodic melting under such conditions. Because ice migration to higher altitudes is a robust mechanism for recharging highland water sources after such events, we suggest that this globally sub-zero, `icy highlands' scenario for the late Noachian climate may be sufficient to explain most of the fluvial geology without the need to invoke additional long-term warming mechanisms or an early warm, wet Mars.

It has been suggested that the recently discovered exoplanet GJ581d might be able to support liquid water due to its relatively low mass and orbital distance. However, GJ581d receives 35% less stellar energy than Mars and is probably locked in tidal resonance, with extremely low insolation at the poles and possibly a permanent night side. Under such conditions, it is unknown whether any habitable climate on the planet would be able to withstand global glaciation and / or atmospheric collapse. Here we present three-dimensional climate simulations that demonstrate GJ581d will have a stable atmosphere and surface liquid water for a wide range of plausible cases, making it the first confirmed super-Earth (exoplanet of 2-10 Earth masses) in the habitable zone. We find that atmospheres with over 10 bar CO2 and varying amounts of background gas (e.g., N2) yield global mean temperatures above 0 degrees Celsius for both land and ocean-covered surfaces. Based on the emitted IR radiation calculated by the model, we propose observational tests that will allow these cases to be distinguished from other possible scenarios in the future.