The Composite InfraRed Spectrometer (CIRS) on board Cassini revealed an equatorial oscillation of stratospheric temperature, reminiscent of the Earth's Quasi-Biennial Oscillation (QBO), as well as anomalously high temperatures under Saturn's rings. To better understand these predominant features of Saturn's atmospheric circulation in the stratosphere, we have extended towards higher altitudes the DYNAMICO-Saturn global climate model (GCM), already used in a previous publication to study the tropospheric dynamics, jets formation and planetary-scale waves activity. Firstly, we study the higher model top impact on the tropospheric zonal jets and kinetic energy distribution. Raising the model top prevents energy and enstrophy accumulation at tropopause levels. The reference GCM simulation with 1/2$^{\circ}$ latitude/longitude resolution and a raised model top exhibits a QBO-like oscillation produced by resolved planetary-scale waves. However, the period is more irregular and the downward propagation faster than observations. Furthermore, compared to the CIRS observation retrievals, the modeled QBO-like oscillation underestimates by half both the amplitude of temperature anomalies at the equator and the vertical characteristic length of this equatorial oscillation. This QBO-like oscillation is mainly driven by westward-propagating waves; a significant lack of eastward wave-forcing explains a fluctuating eastward phase of the QBO-like oscillation. At 20$^{\circ}$N and 20$^{\circ}$S latitudes, the DYNAMICO-Saturn GCM exhibits several strong seasonal eastward jets, alternatively in the northern and southern hemisphere. These jets are correlated with the rings' shadowing. Using a GCM simulation without rings' shadowing, we show its impact on Saturn's stratospheric dynamics. Both residual-mean circulation and eddy forcing are impacted by rings' shadowing.

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

The Europlanet-2020 programme, which ended on Aug 31st, 2019, included an activity called VESPA (Virtual European Solar and Planetary Access), which focused on adapting Virtual Observatory (VO) techniques to handle Planetary Science data. This paper describes some aspects of VESPA at the end of this 4-years development phase and at the onset of the newly selected Europlanet-2024 programme starting in 2020. The main objectives of VESPA are to facilitate searches both in big archives and in small databases, to enable data analysis by providing simple data access and online visualization functions, and to allow research teams to publish derived data in an interoperable environment as easily as possible. VESPA encompasses a wide scope, including surfaces, atmospheres, magnetospheres and planetary plasmas, small bodies, helio-physics, exoplanets, and spectroscopy in solid phase. This system relies in particular on standards and tools developed for the Astronomy VO (IVOA) and extends them where required to handle specificities of Solar System studies. It also aims at making the VO compatible with tools and protocols developed in different contexts, for instance GIS for planetary surfaces, or time series tools for plasma-related measurements. An essential part of the activity is to publish a significant amount of high-quality data in this system, with a focus on derived products resulting from data analysis or simulations.

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.

The Cassini mission unveiled the intense and diverse activity in Saturn's atmosphere: banded jets, waves, vortices, equatorial oscillations. To set the path towards a better understanding of those phenomena, we performed high-resolution multi-annual numerical simulations of Saturn's atmospheric dynamics. We built a new Global Climate Model [GCM] for Saturn, named the Saturn DYNAMICO GCM, by combining a radiative-seasonal model tailored for Saturn to a hydrodynamical solver based on an icosahedral grid suitable for massively-parallel architectures. The impact of numerical dissipation, and the conservation of angular momentum, are examined in the model before a reference simulation employing the Saturn DYNAMICO GCM with a $1/2^{\circ}$ latitude-longitude resolution is considered for analysis. Mid-latitude banded jets showing similarity with observations are reproduced by our model. Those jets are accelerated and maintained by eddy momentum transfers to the mean flow, with the magnitude of momentum fluxes compliant with the observed values. The eddy activity is not regularly distributed with time, but appears as bursts; both barotropic and baroclinic instabilities could play a role in the eddy activity. The steady-state latitude of occurrence of jets is controlled by poleward migration during the spin-up of our model. At the equator, a weakly-superrotating tropospheric jet and vertically-stacked alternating stratospheric jets are obtained in our GCM simulations. The model produces Yanai (Rossby-gravity), Rossby and Kelvin waves at the equator, as well as extratropical Rossby waves, and large-scale vortices in polar regions. Challenges remain to reproduce Saturn's powerful superrotating jet and hexagon-shaped circumpolar jet in the troposphere, and downward-propagating equatorial oscillation in the stratosphere.

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 VESPA data access system focuses on applying Virtual Observatory (VO) standards and tools to Planetary Science. Building on a previous EC-funded Europlanet program, it has reached maturity during the first year of a new Europlanet 2020 program (started in 2015 for 4 years). The infrastructure has been upgraded to handle many fields of Solar System studies, with a focus both on users and data providers. This paper describes the broad lines of the current VESPA infrastructure as seen by a potential user, and provides examples of real use cases in several thematic areas. These use cases are also intended to identify hints for future developments and adaptations of VO tools to Planetary Science.

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.

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.

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.

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.

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.

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.