Wind is the process that connects Mars' climate system. Measurements of Mars atmospheric winds from orbit would dramatically advance our understanding of Mars and help prepare for human exploration of the Red Planet. Multiple instrument candidates are in development and will be ready for flight in the next decade. We urge the Decadal Survey to make these measurements a priority for 2023-2032.

The authors and co-signers of the Terrestrial Planets Comparative Climatology (TPCC) mission concept white paper advocate that planetary science in the next decade would greatly benefit from comparatively studying the fundamental behavior of the atmospheres of Venus and Mars, contemporaneously and with the same instrumentation, to capture atmospheric response to the same solar forcing, and with a minimum of instrument-related variability. Thus, this white paper was created for the 2023-2032 Planetary Science Decadal Survey process. It describes the science rationale for such a mission, and a mission concept that could achieve such a mission.

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.

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.

Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.

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.

The Akatsuki instrument LIR measured an unprecedented wave feature at the top of Venusian cloud layer. Stationary bow-shape waves of thousands of kilometers large lasting several Earth days have been observed over the main equatorial mountains. Here we use for the first time a mesoscale model of the Venus's atmosphere with high-resolution topography and fully coupled interactive radiative transfer computations. Mountain waves resolved by the model form large-scale bow shape waves with an amplitude of about 1.5 K and a size up to several decades of latitude similar to the ones measured by the Akatsuki spacecraft. The maximum amplitude of the waves appears in the afternoon due to an increase of the near-surface stability. Propagating vertically the waves encounter two regions of low static stability, the mixed layer between approximately 18 and 30 km and the convective layer between 50 and 55 km. Some part of the wave energy can pass through these regions via wave tunneling. These two layers act as wave filter, especially the deep atmosphere layer. The encounter with these layers generates trapped lee waves propagating horizontally. No stationary waves is resolved at cloud top over the polar regions because of strong circumpolar transient waves, and a thicker deep atmosphere mixed layer that filters most the mountain waves

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.

Jupiter's upper troposphere and stratosphere is host to a rich dynamical and chemical activity, modulating the thermal structure and distribution of trace species and aerosols, which in turn impact the radiative budget and dynamics. We detail a computationally efficient 1-D radiative-convective model of Jupiter's atmosphere taking into account radiative forcings from the main hydrocarbons (methane, ethane, acetylene), ammonia, collision-induced absorption, several tropospheric and stratospheric cloud and haze layers and an internal heat flux. Sensitivity studies to several parameters are presented and we discuss the expected seasonal, vertical and meridional temperature variations in Jupiter's atmosphere. We compare the resulting thermal structure with that derived from thermal infrared observations by Fletcher et al., 2016. The modeled temperature is in excellent agreement with observations at and near the 10-mbar pressure level, where auroral aerosols warm the atmosphere by up to 20K at latitudes 45-60{\deg}. Those aerosols act to shorten the radiative timescales, estimated to 100 (Earth) days at the 10-mbar level. At pressures lower than 3 mbar, our modeled temperatures are underestimated by ~5K compared to observations, suggesting that other processes, such as dynamical heating by gravity wave breaking, eddies, or thermospheric circulation, play an important role. We then exploit knowledge of heating and cooling rates to diagnose the residual-mean circulation in Jupiter's stratosphere, under the assumption that the eddy heat flux convergence term is negligible. In the lower stratosphere (1-30 mbar), the residual-mean circulation is characterized by upwelling at high latitudes with vertical velocities of 0.3 mm.s-1 and subsidence at 10{\deg}S. In the middle and upper stratosphere, the circulation features rising motion at the equator and subsidence at both poles.

The atmospheric pressure fluctuations on Mars induce an elastic response in the ground that creates a ground tilt, detectable as a seismic signal on the InSight seismometer SEIS. The seismic pressure noise is modeled using Large Eddy Simulations of the wind and surface pressure at the InSight landing site and a Green's function ground deformation approach that is subsequently validated via a detailed comparison with two other methods based on Sorrells' theory (Sorrels 1971; Sorrels et al. 1971). The horizontal acceleration as a result of the ground tilt due to the LES turbulence-induced pressure fluctuations are found to be typically ~2 - 40 nm/s^2 in amplitude, whereas the direct horizontal acceleration is two orders of magnitude smaller and is thus negligible in comparison. The vertical accelerations are found to be ~0.1 - 6 nm/s^2 in amplitude. We show that under calm conditions, a single-pressure measurement is representative of the large-scale pressure field (to a distance of several kilometers), particularly in the prevailing wind direction. However, during windy conditions, small-scale turbulence results in a reduced correlation between the pressure signals, and the single-pressure measurement becomes less representative of the pressure field. Nonetheless, the correlation between the seismic signal and the pressure signal is found to be higher for the windiest period because the seismic pressure noise reflects the atmospheric structure close to the seismometer. In the same way that we reduce the atmospheric seismic signal by making use of a pressure sensor that is part of the InSight APSS (Auxiliary Payload Sensor Suite), we also the use the synthetic noise data obtained from the LES pressure field to demonstrate a decorrelation strategy.

We report on spectroscopic observations of Saturn's stratosphere in July 2011 with the Texas Echelon Cross Echelle Spectrograph (TEXES) mounted on the NASA InfraRed Telescope Facility (IRTF). The observations, targeting several lines of the CH$_4$ $\nu_4$ band and the H$_2$ S(1) quadrupolar line, were designed to determine how Saturn's stratospheric thermal structure was disturbed by the 2010 Great White Spot. A study of Cassini Composite Infrared Spectrometer (CIRS) spectra had already shown the presence of a large stratospheric disturbance centered at a pressure of 2~hPa, nicknamed the beacon B0, and a tail of warm air at lower pressures (Fletcher et al. 2012. Icarus 221, 560--586). Our observations confirm that the beacon B0 vertical structure determined by CIRS, with a maximum temperature of $180\pm1$K at 2~hPa, is overlain by a temperature decrease up to the 0.2-hPa pressure level. Our retrieved maximum temperature of $180\pm1$K is colder than that derived by CIRS ($200\pm1$K), a difference that may be quantitatively explained by terrestrial atmospheric smearing. We propose a scenario for the formation of the beacon based on the saturation of gravity waves emitted by the GWS. Our observations also reveal that the tail is a planet-encircling disturbance in Saturn's upper stratosphere, oscillating between 0.2 and 0.02~hPa, showing a distinct wavenumber-2 pattern. We propose that this pattern in the upper stratosphere is either the signature of thermal tides generated by the presence of the warm beacon in the mid-stratosphere, or the signature of Rossby wave activity.

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.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

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.

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.

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.