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.

Studying the atmospheric Planetary Boundary Layer (PBL) is crucial to understand the climate of a planet. The meteorological measurements by the instruments onboard InSight make a uniquely rich dataset to study the active turbulent dynamics of the daytime PBL on Mars. Here we use the high-sensitivity continuous pressure, wind, temperature measurements in the first 400 sols of InSight operations to analyze convective gusts, cells, and vortices in Mars' daytime PBL. We compare the InSight measurements to turbulence-resolving Large-Eddy Simulations (LES). The daytime PBL turbulence at the InSight landing site is found to be very active, with clearly identified signatures of convective cells and a vast population of almost 10,000 recorded vortex encounters, adequately represented by a power-law with a 3.5 exponent in agreement with LES results. While the daily variability of vortex encounters at the InSight landing site can be explained by the statistical nature of turbulence, the seasonal variability is strongly correlated with ambient wind speed, which is supported by LES. However, wind gustiness is more correlated to surface temperature than ambient wind speed, confirming the radiative control of the daytime martian PBL; and fewer convective vortices are forming in LES when the background wind is doubled. Thus, the long-term seasonal variability of vortex encounters at the InSight landing site is mainly controlled by the advection of convective vortices by the ambient wind speed. Typical tracks followed by vortices forming in the LES show a similar distribution in direction and length as orbital imagery of the InSight region.

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.

We extracted a Martian water-ice cloud climatology from OMEGA data covering 7 Martian years (MY 26-32). We derived two products, the Reversed Ice Cloud Index (ICIR) and the Percentage of Cloudy Pixels (PCP), indicating the mean cloud thickness and nebulosity over a regular grid (1{\deg} longitude x 1{\deg} latitude x 1{\deg} Ls x 1 h Local Time). The ICIR has been shown to be a proxy of the water-ice column derived from the Mars Climate Database. The PCP confirms the location of the main cloud structures mapped with the ICIR, and gives a more accurate image of the cloud cover. We observed a denser cloud coverage over Hellas Planitia, Lunae Planum and over large volcanoes in the aphelion belt. For the first time, thanks to the fact that Mars Express is not in Sun-synchronous orbit, we can explore the cloud diurnal cycle at a given season by combining 7 years of observations. However, because of the eccentric orbit, the temporal coverage remains limited. Other limitations of the dataset are its small size, the difficult distinction between ice clouds and frosts, and the impact of surface albedo on data uncertainty. We could nevertheless study the diurnal cloud life cycle by averaging the data over larger regions: from specific topographic features (covering a few degrees in longitude and latitude) up to large climatic bands (all longitudes). We found that in the tropics around northern summer solstice, the diurnal thermal tide modulates the abundance of clouds, which is reduced around noon. At northern midlatitudes, clouds corresponding to the edge of the north polar hood are observed mainly in the morning and around noon during northern winter (Ls=260-30{\deg}). Over Chryse Planitia, low lying morning fogs dissipate earlier and earlier in the afternoon during northern winter. Over Argyre, clouds are present over all daytime during two periods, around Ls = 30 and 160{\deg}.

The impact of gravity waves (GW) on diurnal tides and the global circulation in the middle/upper atmosphere of Mars is investigated using a General Circulation Model (GCM). We have implemented a stochastic parameterization of non-orographic GW into the Laboratoire de M\'et\'eorologie Dynamique (LMD) Mars GCM (LMD-MGCM) following an innovative approach. The source is assumed to be located above typical convective cells ($\sim$ 250 Pa) and the effect of GW on the circulation and predicted thermal structure above 1 Pa ($\sim$ 50 km) is analyzed. We focus on the comparison between model simulations and observations by the Mars Climate Sounder (MCS) on board Mars Reconnaissance Orbiter during Martian Year 29. MCS data provide the only systematic measurements of the Martian mesosphere up to 80 km to date. The primary effect of GW is to damp the thermal tides by reducing the diurnal oscillation of the meridional and zonal winds. The GW drag reaches magnitudes of the order of 1 m/s/sol above 10$^{-2}$ Pa in the northern hemisphere winter solstice and produces major changes in the zonal wind field (from tens to hundreds of m/s), while the impact on the temperature field is relatively moderate (10-20K). It suggests that GW induced alteration of the meridional flow is the main responsible for the simulated temperature variation. The results also show that with the GW scheme included, the maximum day-night temperature difference due to the diurnal tide is around 10K, and the peak of the tide is shifted toward lower altitudes, in better agreement with MCS observations.

K2-18 b is a transiting mini-Neptune that orbits a nearby (38 pc) cool M3 dwarf and is located inside its region of temperate irradiation. We report on the search for hydrogen escape from the atmosphere K2-18 b using Lyman-$\alpha$ transit spectroscopy with the Space Telescope Imaging Spectrograph (STIS) instrument installed on the Hubble Space Telescope (HST). We analyzed the time-series of the fluxes of the stellar Lyman-$\alpha$ emission of K2-18 in both its blue- and redshifted wings. We found that the average blueshifted emission of K2-18 decreases by $67\% \pm 18\%$ during the transit of the planet compared to the pre-transit emission, tentatively indicating the presence of H atoms escaping vigorously and being blown away by radiation pressure. This interpretation is not definitive because it relies on one partial transit. Based on the reconstructed Lyman-$\alpha$ emission of K2-18, we estimate an EUV irradiation between $10^1-10^2$ erg s$^{-1}$ cm$^{-2}$ and a total escape rate in the order of $10^8$ g s$^{-1}$. The inferred escape rate suggests that the planet will lose only a small fraction (< 1%) of its mass and retain its volatile-rich atmosphere during its lifetime. More observations are needed to rule out stellar variability effects, confirm the in-transit absorption and better assess the atmospheric escape and high-energy environment of K2-18 b.

Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an inter-comparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (20-30 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance.

By measuring the regular oscillations of the density of CO$_{2}$ in the upper atmosphere (between 120 and 190~km), the mass spectrometer MAVEN/NGIMS (Atmosphere and Volatile EvolutioN/Neutral Gas Ion Mass Spectrometer) reveals the local impact of gravity waves (GWs). This yields precious information on the activity of GWs and the atmospheric conditions in which they propagate and break. The intensity of GWs measured by MAVEN in the upper atmosphere has been shown to be dictated by saturation processes in isothermal conditions. As a result, GWs activity is correlated to the evolution of the inverse of the background temperature. Previous data gathered at lower altitudes ($\sim$95 to $\sim$150~km) during aerobraking by the accelerometers on board MGS (Mars Global Surveyor), ODY (Mars Odyssey) and MRO (Mars Reconnaissance Orbiter) are analyzed in the light of those recent findings with MAVEN. The anti-correlation between GW-induced density perturbations and background temperature is plausibly found in the ODY data acquired in the polar regions, but not in he MGS and MRO data. MRO data in polar regions exhibit a correlation between the density perturbations and the Brunt-V\"{a}is\"{a}l\"{a} frequency, obtained from Global Climate Modeling. At lower altitude levels (between 100 and 120~km), although wave saturation might still be dominant, isothermal conditions are no longer verified. In this case, theory predicts that the intensity of GWs is no more correlated to background temperature, but to static stability. At other latitudes in the three aerobraking datasets, the GW-induced relative density perturbations are correlated with neither inverse temperature nor static stability; in this particular case, this means that the observed activity of GWs is not only controlled by saturation, but also by the effects of GWs sources and wind filtering through critical levels.

We have reconstructed longitude-latitude maps of column dust optical depth (CDOD) for Martian year (MY) 34 (May 5, 2017 --- March 23, 2019) using observations by the Mars Climate Sounder (MCS) aboard NASA's Mars Reconnaissance Orbiter spacecraft. Our methodology works by gridding standard and newly available estimates of CDOD from MCS limb observations, using the "iterative weighted binning" methodology. In this work, we reconstruct four gridded CDOD maps per sol, at different Mars Universal Times. Together with the seasonal and day-to-day variability, the use of several maps per sol allows to explore also the daily variability of CDOD in the MCS dataset, which is shown to be particularly strong during the MY 34 equinoctial Global Dust Event (GDE). Regular maps of CDOD are then produced by daily averaging and spatially interpolating the irregularly gridded maps using a standard "kriging" interpolator, and can be used as "dust scenario" for numerical model simulations. In order to understand whether the daily variability of CDOD has a physical explanation, we have carried out numerical simulations with the "Laboratoire de M\'et\'eorologie Dynamique" Mars Global Climate Model. Using a "free dust" run initiated at $L_s \sim 210^\circ$ with the corresponding kriged map, but subsequently free of further CDOD forcing, we show that the model is able to account for some of the observed daily variability in CDOD. The model serves also to confirm that the use of the MY 34 daily-averaged dust scenario in a GCM produces results consistent with those obtained for the MY 25 GDE.

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.

Pluto's atmospheric haze settles out rapidly compared with geological timescales. It needs to be accounted for as a surface material, distinct from Pluto's icy bedrock and from the volatile ices that migrate via sublimation and condensation on seasonal timescales. This paper explores how a steady supply of atmospheric haze might affect three distinct provinces on Pluto. We pose the question of why they each look so different from one another if the same haze material is settling out onto all of them. Cthulhu is a more ancient region with comparatively little present-day geological activity, where the haze appears to simply accumulate over time. Sputnik Planitia is a very active region where glacial convection, as well as sublimation and condensation rapidly refresh the surface, hiding recently deposited haze from view. Lowell Regio is a region of intermediate age featuring very distinct coloration from the rest of Pluto. Using a simple model haze particle as a colorant, we are not able to match the colors in both Lowell Regio and Cthulhu. To account for their distinct colors, we propose that after arrival at Pluto's surface, haze particles may be less inert than might be supposed from the low surface temperatures. They must either interact with local materials and environments to produce distinct products in different regions, or else the supply of haze must be non-uniform in time and/or location, such that different products are delivered to different places.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

New Horizons observations suggest that CH4 on Pluto has a complex history, involving reservoirs of different composition, thickness and stability controlled by volatile processes occurring on different timescales. In order to interpret these observations, we use a Pluto volatile transport model able to simulate the cycles of N2 and CH4 ices over millions of years. By assuming fixed solid mixing ratios, we explore how changes in surface albedos, emissivities and thermal inertias impact volatile transport. This work is therefore a direct and natural continuation of the work by Bertrand et al. (2018), which only explored the N2 cycles. Results show that bright CH4 deposits can create cold traps for N2 ice outside Sputnik Planitia, leading to a strong coupling between the N2 and CH4 cycles. Depending on the assumed albedo for CH4 ice, the model predicts CH4 ice accumulation (1) at the same equatorial latitudes where the Bladed Terrain Deposits are observed, supporting the idea that these CH4-rich deposits are massive and perennial, or (2) at mid-latitudes (25{\deg}N-70{\deg}N), forming a thick mantle which is consistent with New Horizons observations. In our simulations, both CH4 ice reservoirs are not in an equilibrium state and either one can dominate the other over long timescales, depending on the assumptions made for the CH4 albedo. This suggests that long-term volatile transport exists between the observed reservoirs. The model also reproduces the formation of N2 deposits at mid-latitudes and in the equatorial depressions surrounding the Bladed Terrain, as observed by New Horizons. At the poles, only seasonal CH4 and N2 deposits are obtained in Pluto's current orbital configuration. Finally, we show that Pluto's atmosphere always contained, over the last astronomical cycles, enough gaseous CH4 to absorb most of the incoming Lyman-flux.

We investigate the global seasonal variations of near-surface relative humidity and relevant attributes, like temperature and water vapor volume mixing ratio on Mars using calculations from modelled and measurement data. We focus on 2 am local time snapshots to eliminate daily effects related to differences in insolation, and to be able to compare calculations based on modelling data from the LMDZ GCM with the observations of MGS TES. We study the seasonal effects by examining four specific dates in the Martian year, the northern spring equinox, summer solstice, autumn equinox and winter solstice. We identify three specific zones, where the near-surface relative humidity levels are systematically higher than in their vicinity regardless of season. We find that these areas coincide with low thermal inertia features, which control surface temperatures on the planet, and are most likely covered with unconsolidated fine dust with grain sizes less than $\sim$ 40$\mu$m. By comparing the data of relative humidity, temperature and water vapor volume mixing ratio at two different heights (near-surface, $\sim$ 23 m above the surface), we demonstrate that the thermal inertia could play an important role in determining near-surface humidity levels. We also notice that during the night the water vapor levels drop at $\sim$ 4 m above the surface. This, together with the temperature and thermal inertia values, shows that water vapor likely condenses in the near-surface atmosphere and on the ground during the night at the three aforementioned regions. This condensation may be in the form of brines, wettening of the fine grains or deliquescence. This study specifies areas of interest on the surface of present day Mars for the proposed condensation, which may be examined by in-situ measurements in the future.

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.

An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional (1D) radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G- and M-stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 um) and in the region between 0.2 and 1.5 um. Differences in outgoing longwave radiation increase with surface temperature and reach 10-20 Wm^-2; differences in shortwave reach up to 60 Wm^-2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ~10% of modern Earth's solar constant (i.e., ~34 Wm^-2 in global mean) among band models and ~3% between the two line-by-line models. These comparisons show that future work is needed focusing on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

Pluto's landscape is shaped by the cycles of the volatile ices covering its surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large diversity of terrains, with bright and dark plains, pits, topographic depressions and evidences of recent and past glacial flows. Outside SP, New Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at mid-northern latitudes. These observations suggest a complex history involving volatile and glacial processes on different timescales. We present numerical simulations of volatile transport on Pluto performed with a model able to simulate the N2 cycle over millions of years (Myrs), taking into account the changes of obliquity and orbital parameters as experienced by Pluto. Results show that over one obliquity cycle, the latitudes of SP between 25{\deg}S-30{\deg}N are dominated by N2 condensation, while the latitudes between 30-50{\deg}N are dominated by N2 sublimation. We find that a net amount of 1 km of ice has sublimed at the northern edge of SP during the last 2 Myrs. By comparing these results with the observed geology of SP, we can relate the formation of the pits and the brightness of the ice to the ice flux occurring at the annual timescale, while the glacial flows at its eastern edge and the erosion of the water ice mountains all around the ice sheet are related to the astronomical timescale. We also perform simulations with a glacial flow scheme which shows that SP is currently at its minimum extent. We also explore the stability of N2 ice outside SP. Results show that it is not stable at the poles but rather in the equatorial regions, in particular in depressions, where thick deposits may persist over tens of Myrs, before being trapped in SP. Finally, another key result is that the minimum and maximum surface pressures obtained over the simulated Myrs remain in the range of mm-Pa and Pa, respectively.

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.

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.

The New Horizons spacecraft, which flew by Pluto on July 14, 2015, revealed the presence of haze in Pluto's atmosphere that were formed by CH4/N2 photochemistry at high altitudes in Pluto's atmosphere, as on Titan and Triton. In order to help the analysis of the observations and further investigate the formation of organic haze and its evolution at global scales, we have implemented a simple parametrization of the formation of organic haze in our Pluto General Circulation Model. The production of haze in our model is based on the different steps of aerosol formation as understood on Titan and Triton: photolysis of CH4 in the upper atmosphere by Lyman-alpha UV radiation, production of various gaseous species, and conversion into solid particles through accumulation and aggregation processes. We compared two reference simulations ran with a particle radius of 50 nm: with, and without South Pole N2 condensation. We discuss the impact of the particle radius and the lifetime of the precursors on the haze distribution. We simulate CH4 photolysis and the haze formation up to 600 km above the surface. Results show that CH4 photolysis in Pluto's atmosphere in 2015 occured mostly in the sunlit summer hemisphere with a peak at an altitude of 250 km, though the interplanetary source of Lyman-alpha flux can induce some photolysis even in the Winter hemisphere. We obtained an extensive haze up to altitudes comparable with the observations, and with non-negligible densities up to 500 km altitude. In both reference simulations, the haze density is not strongly impacted by the meridional circulation. With no South Pole N2 condensation, the maximum nadir opacity and haze extent is obtained at the North Pole. With South Pole N2 condensation, the descending parcel of air above the South Pole leads to a latitudinally more homogeneous haze density with a slight density peak at the South Pole.

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.

We use a 3D general circulation model to compare the primitive Martian hydrological cycle in "warm and wet" and "cold and icy" scenarios. In the warm and wet scenario, an anomalously high solar flux or intense greenhouse warming artificially added to the climate model are required to maintain warm conditions and an ice-free northern ocean. Precipitation shows strong surface variations, with high rates around Hellas basin and west of Tharsis but low rates around Margaritifer Sinus (where the observed valley network drainage density is nonetheless high). In the cold and icy scenario, snow migration is a function of both obliquity and surface pressure, and limited episodic melting is possible through combinations of seasonal, volcanic and impact forcing. At surface pressures above those required to avoid atmospheric collapse (~0.5 bar) and moderate to high obliquity, snow is transported to the equatorial highland regions where the concentration of valley networks is highest. Snow accumulation in the Aeolis quadrangle is high, indicating an ice-free northern ocean is not required to supply water to Gale crater. At lower surface pressures and obliquities, both H2O and CO2 are trapped as ice at the poles and the equatorial regions become extremely dry. The valley network distribution is positively correlated with snow accumulation produced by the cold and icy simulation at 41.8 degrees obliquity but uncorrelated with precipitation produced by the warm and wet simulation. Because our simulations make specific predictions for precipitation patterns under different climate scenarios, they motivate future targeted geological studies.

Observations show that the lower thermosphere of Mars ($\sim$100--140 km) is up to 40 K colder than the current general circulation models (GCMs) can reproduce. Possible candidates for physical processes missing in the models are larger abundances of atomic oxygen facilitating stronger CO$_2$ radiative cooling, and thermal effects of gravity waves. Using two state-of-the-art Martian GCMs, the Laboratoire de M\'et\'eorologie Dynamique and Max Planck Institute models that self-consistently cover the atmosphere from the surface to the thermosphere, these physical mechanisms are investigated. Simulations demonstrate that the CO$_2$ radiative cooling with a sufficiently large atomic oxygen abundance, and the gravity wave-induced cooling can alone result in up to 40 K colder temperature in the lower thermosphere. Accounting for both mechanisms produce stronger cooling at high latitudes. However, radiative cooling effects peak above the mesopause, while gravity wave cooling rates continuously increase with height. Although both mechanisms act simultaneously, these peculiarities could help to further quantify their relative contributions from future observations.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We have produced a multiannual climatology of airborne dust from Martian year 24 to 31 using multiple datasets of retrieved or estimated column optical depths. The datasets are based on observations of the Martian atmosphere from April 1999 to July 2013 made by different orbiting instruments: the Thermal Emission Spectrometer (TES) aboard Mars Global Surveyor, the Thermal Emission Imaging System (THEMIS) aboard Mars Odyssey, and the Mars Climate Sounder (MCS) aboard Mars Reconnaissance Orbiter (MRO). The procedure we have adopted consists of gridding the available retrievals of column dust optical depth (CDOD) from TES and THEMIS nadir observations, as well as the estimates of this quantity from MCS limb observations. Our gridding method calculates averages and uncertainties on a regularly spaced, but possibly incomplete, spatio-temporal grid, using an iterative procedure weighted in space, time, and retrieval uncertainty. In order to evaluate strengths and weaknesses of the resulting gridded maps, we validate them with independent observations of CDOD. We have statistically analyzed the irregularly gridded maps to provide an overview of the dust climatology on Mars over eight years, specifically in relation to its interseasonal and interannual variability. Finally, we have produced multiannual, regular daily maps of CDOD by spatially interpolating the irregularly gridded maps using a kriging method. These synoptic maps are used as dust scenarios in the Mars Climate Database version 5, and are useful in many modelling applications in addition to forming a basis for instrument intercomparisons. The derived dust maps for the eight available Martian years are publicly available and distributed with open access.

Radiative effects of water ice clouds have noteworthy consequences on the Martian atmosphere, its thermal structure and circulation. Accordingly, the inclusion of such effects in the LMD Mars Global Climate Model (GCM) greatly modifies the simulated Martian water cycle. The intent of this paper is to address the impact of radiatively active clouds on atmospheric water vapor and ice in the GCM and improve its representation. We propose a new enhanced modeling of the water cycle, consisting of detailed cloud microphysics with dynamic condensation nuclei and a better implementation of perennial surface water ice. This physical modeling is based on tunable parameters. This new version of the GCM is compared to the Thermal Emission Spectrometer observations of the water cycle. Satisfying results are reached for both vapor and cloud opacities. However, simulations yield a lack of water vapor in the tropics after Ls=180{\deg} which is persistent in simulations compared to observations, as a consequence of aphelion cloud radiative effects strengthening the Hadley cell. Every year, our GCM simulations indicate that permanent surface water ice on the north polar cap increases at latitudes higher than 80{\deg}N and decreases at lower latitudes. Supersaturation above the hygropause is predicted in line with SPICAM observations. The model also shows for the first time that the scavenging of dust by water ice clouds alone fails to fully account for observed dust detached layers.

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.

High-resolution spectra of Pluto in the 1.66 um region, recorded with the VLT/CRIRES instrument in 2008 (2 spectra) and 2012 (5 spectra), are analyzed to constrain the spatial and vertical distribution of methane in Pluto's atmosphere and to search for mid-term (4 year) variability. A sensitivity study to model assumptions (temperature structure, surface pressure, Pluto's radius) is performed. Results indicate that (i) no variation of the CH4 atmospheric content (column-density or mixing ratio) with Pluto rotational phase is present in excess of 20 % (ii) CH4 column densities show at most marginal variations between 2008 and 2012, with a best guess estimate of a ~20 % decrease over this time frame. As stellar occultations indicate that Pluto's surface pressure has continued to increase over this period, this implies a concomitant decrease of the methane mixing ratio (iii) the data do not show evidence for an altitude-varying methane distribution; in particular, they imply a roughly uniform mixing ratio in at least the first 22-27 km of the atmosphere, and high concentrations of low-temperature methane near the surface can be ruled out. Our results are also best consistent with a relatively large (> 1180 km) Pluto radius. Comparison with predictions from a recently developed global climate model GCM indicates that these features are best explained if the source of methane occurs in regional-scale CH4 ice deposits, including both low latitudes and high Northern latitudes, evidence for which is present from the rotational and secular evolution of the near-IR features due to CH4 ice. Our "best guess" predictions for the New Horizons encounter in 2015 are: a 1184 km radius, a 17 ubar surface pressure, and a 0.44 % CH4 mixing ratio with negligible longitudinal variations.

Super-Earths with solid surfaces such CoRot-7b and Kepler-10b are expected to be extremely hot. It has been suggested that they must be atmosphere-free and that a lava ocean is present on their hot dayside. Here, we use several dedicated thermal models to study how observations with NIRSPEC on the JWST could further confirm and constrain, or reject the atmosphere-free lava ocean planet model for very hot super-Earths. Using CoRoT-7b as a working case, we explore the consequences on the phase-curve of a non tidal-locked rotation, with/without an atmosphere, and for different values of the albedo. We simulate future observations of the reflected light and thermal emission from CoRoT-7b with NIRSPEC-JWST and look for detectable signatures, such as time lag, as well as the possibility to retrieve the latitudinal surface temperature distribution. We demonstrate that we should be able to constrain several parameters after observations of two orbits (42 h) with a broad range of wavelengths: i)The Bond albedo is retrieved to within +/- 0.03 in most cases. ii) The lag effect allows us to retrieve the rotation period within 3 hours for a planet, with rotation half the orbital period. iii) Any spin period shorter than a limit in the range 30 - 800 h, depending on the thickness of the thermal layer in the soil, would be detected. iv) The presence of a thick gray atmosphere with a pressure of one bar, and a specific opacity higher than 1.0E-5 m-2 kg-1 is detectable. v) With spectra up to 4.5 mu, the latitudinal temperature profile can be retrieved to within 30 K with a risk of a totally wrong solution in 5 % of the cases. We conclude that it should thus be possible to distinguish the reference situation from other cases. Surface temperature map and the albedo brings important constraints on the nature or the physical state of the soil of hot super-Earths.

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.

What kind of environment may exist on terrestrial planets around other stars? In spite of the lack of direct observations, it may not be premature to speculate on exoplanetary climates, for instance to optimize future telescopic observations, or to assess the probability of habitable worlds. To first order, climate primarily depends on 1) The atmospheric composition and the volatile inventory; 2) The incident stellar flux; 3) The tidal evolution of the planetary spin, which can notably lock a planet with a permanent night side. The atmospheric composition and mass depends on complex processes which are difficult to model: origins of volatile, atmospheric escape, geochemistry, photochemistry. We discuss physical constraints which can help us to speculate on the possible type of atmosphere, depending on the planet size, its final distance for its star and the star type. Assuming that the atmosphere is known, the possible climates can be explored using Global Climate Models analogous to the ones developed to simulate the Earth as well as the other telluric atmospheres in the solar system. Our experience with Mars, Titan and Venus suggests that realistic climate simulators can be developed by combining components like a "dynamical core", a radiative transfer solver, a parametrisation of subgrid-scale turbulence and convection, a thermal ground model, and a volatile phase change code. On this basis, we can aspire to build reliable climate predictors for exoplanets. However, whatever the accuracy of the models, predicting the actual climate regime on a specific planet will remain challenging because climate systems are affected by strong positive destabilizing feedbacks (such as runaway glaciations and runaway greenhouse effect). They can drive planets with very similar forcing and volatile inventory to completely different states.

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 Martian Planetary Boundary Layer [PBL] is a crucial component of the Martian climate system. Global Climate Models [GCMs] and Mesoscale Models [MMs] lack the resolution to predict PBL mixing which is therefore parameterized. Here we propose to adapt the "thermal plume" model, recently developed for Earth climate modeling, to Martian GCMs, MMs, and single-column models. The aim of this physically-based parameterization is to represent the effect of organized turbulent structures (updrafts and downdrafts) on the daytime PBL transport, as it is resolved in Large-Eddy Simulations [LESs]. We find that the terrestrial thermal plume model needs to be modified to satisfyingly account for deep turbulent plumes found in the Martian convective PBL. Our Martian thermal plume model qualitatively and quantitatively reproduces the thermal structure of the daytime PBL on Mars: superadiabatic near-surface layer, mixing layer, and overshoot region at PBL top. This model is coupled to surface layer parameterizations taking into account stability and turbulent gustiness to calculate surface-atmosphere fluxes. Those new parameterizations for the surface and mixed layers are validated against near-surface lander measurements. Using a thermal plume model moreover enables a first order estimation of key turbulent quantities (e.g. PBL height, convective plume velocity) in Martian GCMs and MMs without having to run costly LESs.

In the past 15 years, astronomers have revealed that a significant fraction of the stars should harbor planets and that it is likely that terrestrial planets are abundant in our galaxy. Among these planets, how many are habitable, i.e. suitable for life and its evolution? These questions have been discussed for years and we are slowly making progress. Liquid water remains the key criterion for habitability. It can exist in the interior of a variety of planetary bodies, but it is usually assumed that liquid water at the surface interacting with rocks and light is necessary for the emergence of a life able to modify its environment and evolve. A first key issue is thus to understand the climatic conditions allowing surface liquid water assuming a suitable atmosphere. This have been studied with global mean 1D models which has defined the "classical habitable zone", the range of orbital distances within which worlds can maintain liquid water on their surfaces (Kasting et al. 1993). A new generation of 3D climate models based on universal equations and tested on bodies in the solar system is now available to explore with accuracy climate regimes that could locally allow liquid water. A second key issue is now to better understand the processes which control the composition and the evolution of the atmospheres of exoplanets, and in particular the geophysical feedbacks that seems to be necessary to maintain a continuously habitable climate. From that point of view, it is not impossible that the Earth's case may be special and uncommon.

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.

Airborne dust is the main climatic agent in the Martian environment. Local dust storms play a key role in the dust cycle; yet their life cycle is poorly known. Here we use mesoscale modeling that includes the transport of radiatively active dust to predict the evolution of a local dust storm monitored by OMEGA on board Mars Express. We show that the evolution of this dust storm is governed by deep convective motions. The supply of convective energy is provided by the absorption of incoming sunlight by dust particles, rather than by latent heating as in moist convection on Earth. We propose to use the terminology "rocket dust storm", or conio-cumulonimbus, to describe those storms in which rapid and efficient vertical transport takes place, injecting dust particles at high altitudes in the Martian troposphere (30 to 50 km). Combined to horizontal transport by large-scale winds, rocket dust storms produce detached layers of dust reminiscent of those observed with Mars Global Surveyor and Mars Reconnaissance Orbiter. Since nighttime sedimentation is less efficient than daytime convective transport, and the detached dust layers can convect during the daytime, these layers can be stable for several days. The peak activity of rocket dust storms is expected in low-latitude regions at clear seasons (late northern winter to late northern summer), which accounts for the high-altitude tropical dust maxima unveiled by Mars Climate Sounder. Dust-driven deep convection have strong implications for the Martian dust cycle, thermal structure, atmospheric dynamics, cloud microphysics, chemistry, and robotic and human exploration.

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.

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 CO2 atmosphere with pressure between 0.1 and 7 bars. The model includes a detailed radiative transfer model using revised CO2 gas collision induced absorption properties, and a parameterisation of the CO2 ice cloud microphysical and radiative properties. A wide range of possible climates is explored by 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 the adiabatic cooling and warming of the atmosphere when it moves vertically. In most simulations, CO2 ice clouds cover a major part of the planet but greenhouse effect does not exceed +15 K. We find that a CO2 atmosphere could not have raised the annual mean temperature above 0{\deg}C anywhere on the planet. The collapse of the atmosphere into permanent CO2 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{\deg}C (a condition which could have allowed rivers to form) are predicted for obliquity larger than 40{\deg} at high latitudes but not in locations where most valley networks 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.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

Although transit spectroscopy is a powerful method for studying the composition, thermal properties and dynamics of exoplanet atmospheres, only a few transiting terrestrial exoplanets will be close enough to allow significant transit spectroscopy. Thermal phase curves (variations of the apparent infrared emission of the planet with its orbital phase) have been observed for hot Jupiters in both transiting and nontransiting configurations, and could be observed for hot terrestrial exoplanets. We study the wavelength and phase changes of the thermal emission of a tidally-locked terrestrial planet as atmospheric pressure increases, and address the observability of these multiband phase-curves and the ability to use them to detect atmospheric constituents. We used a 3D climate model (GCM) to simulate the CO2 atmosphere of a terrestrial planet on an 8-day orbit around a M3 dwarf and its apparent infrared emission as a function of its orbital phase. We estimated the signal to photon-noise ratio in narrow bands between 2.5 and 20 microns for a 10 pc target observed with a 6 and 1.5 m telescope (respectively the sizes of JWST and EChO). We find that tmospheric absorption bands produce signatures in what we call the variation spectrum. Planets with no atmosphere produce large variations and can be easily distinguished from dense absorbing atmospheres. Photon-noise limited spectro-photometry of nearby systems could allow us to detect and characterize the atmosphere of nontransiting terrestrial planets known from radial velocity surveys. Two obvious impediments for these types of observations are the required photometric sensitivity (1E-5) over the duration of at least one orbit and the intrinsic stellar variability. However, overcoming these obstacles would give access to one order of magnitude more targets than transit spectroscopy.

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.

Near-surface perennial water ice on Mars has been previously inferred down to latitudes of about 45{\deg} and could result from either water vapor diffusion through the regolith under current conditions or previous ice ages precipitations. In this paper we show that at latitudes as low as 25{\deg} in the southern hemisphere buried water ice in the shallow (< 1 m) subsurface is required to explain the observed surface distribution of seasonal CO2 frost on pole facing slopes. This result shows that possible remnants of the last ice age, as well as water that will be needed for the future exploration of Mars, are accessible significantly closer to the equator than previously thought, where mild conditions for both robotic and human exploration lie.

In this paper, we analyze water ice occurrences at the surface of Mars using near-infrared observations, and we study their distribution with a climate model. Latitudes between 45{\deg}S and 50{\deg}N are considered. Data from the Observatoire pour la Min\'eralogie, l'Eau, les Glaces et l'Actitit\'e and the Compact Reconnaissance Imaging Spectrometer for Mars are used to assess the presence of surface water ice as a function of location and season. A modeling approach combining the 1-D and 3-D versions of the General Circulation Model of the Laboratoire de M\'et\'eorologie Dynamique de Jussieu is developed and successfully compared to observations. Ice deposits 2-200 \mu m thick are observed during the day on pole facing slopes in local fall, winter and early spring. Ice extends down to 13{\deg} latitude in the Southern Hemisphere but is restricted to latitudes higher than 32{\deg} in the north. On a given slope, the pattern of ice observations at the surface is mainly controlled by the global variability of atmospheric water (precipitation and vapor), with local ground properties playing a lower role. Only seasonal surface ice is observed: no exposed patches of perennial ground ice have been detected. Surface seasonal ice is however sensitive to subsurface properties: the results presented in this study are consistent with the recent discovery of low latitude subsurface ice obtained through the analysis of CO2 frost.

The recently discovered exoplanet Gl581d is extremely close to the outer edge of its system's habitable zone, which has led to much speculation on its possible climate. We have performed a range of simulations to assess whether, given simple combinations of chemically stable greenhouse gases, the planet could sustain liquid water on its surface. For best estimates of the surface gravity, surface albedo and cloud coverage, we find that less than 10 bars of CO2 is sufficient to maintain a global mean temperature above the melting point of water. Furthermore, even with the most conservative choices of these parameters, we calculate temperatures above the water melting point for CO2 partial pressures greater than about 40 bar. However, we note that as Gl581d is probably in a tidally resonant orbit, further simulations in 3D are required to test whether such atmospheric conditions are stable against the collapse of CO2 on the surface.

The polar condensation/sublimation of CO2, that involve about one fourth of the atmosphere mass, is the major Martian climatic cycle. Early observations in visible and thermal infrared have shown that the sublimation of the Seasonal South Polar Cap (SSPC) is not symmetric around the geographic South Pole. Here we use observations by OMEGA/Mars Express in the near-infrared to detect unambiguously the presence of CO2 at the surface, and to estimate albedo. Second, we estimate the sublimation of CO2 released in the atmosphere and show that there is a two-step process. From Ls=180{\deg} to 220{\deg}, the sublimation is nearly symmetric with a slight advantage for the cryptic region. After Ls=220{\deg} the anti-cryptic region sublimation is stronger. Those two phases are not balanced such that there is 22% +/- 9 more mass the anti-cryptic region, arguing for more snow precipitation. We compare those results with the MOLA height measurements. Finally we discuss implications for the Martian atmosphere about general circulation and gas tracers, e.g. Ar.

Characterizing the Martian atmosphere is an essential objective to understand its meteorology and its climate. The lower atmosphere (< 40 km) and middle atmosphere (40-80 km) of Mars appear dynamically coupled at much higher levels than in the case of the Earth. The vertical extension of the weather phenomena is considerable with for example Hadley's cells reaching the top of the neutral atmosphere (120 km). The circulation in the middle atmosphere modifies the meteorology of the lower atmosphere, affecting the transport and climatic processes Observations of the CO rotational lines at millimeter (mm) wavelengths (Clancy et al 1990) have strongly contributed in the study of the vertical distribution of this compound and the thermal profile in the atmosphere of Mars over 0-70 km. Singledish observations of the CO Doppler lineshift have allowed direct wind measurements in the martian middle atmosphere near 50 km altitude (Lellouch et al 1991), but at a low spatial resolution (12'') enabling only an essentially hemispheric resolution of the martian disk. The use of mm interferometry has allowed us to better spatially resolve the Martian disk, in order to obtain wind maps of the middle atmosphere (Moreno et al 2001).