This study presents the development of a TKE-l parameterization of the diffusion coefficients for the representation of turbulent diffusion in neutral and stable conditions in large-scale atmospheric models. The parameterization has been carefully designed to be completely tunable in the sense that all adjustable parameters have been clearly identified and their number minimized as much as possible to help the calibration and to thoroughly assess the parametric sensitivity. We choose a mixing length formulation that depends on both static stability and wind shear to cover the different regimes of stable boundary layers. We follow a heuristic approach for expressing the stability functions and turbulent Prandlt number in order to guarantee the versatility of the scheme and its applicability for planetary atmospheres composed of an ideal and perfect gas such as that of Earth and Mars. Particular attention has also been paid to the numerical stability at typical time steps used in General Circulation Models. Test, parametric sensitivity assessment and preliminary tuning are performed on single-column idealized simulations of the weakly stable boundary layer. The robustness and versatility of the scheme are also assessed through its implementation in the LMDZ General Circulation Model and the Mars Planetary Climate Model and by running simulations of the Antarctic and Martian nocturnal boundary layers.

The heat flow and physical properties package measured soil thermal conductivity at the landing site in the 0.03 to 0.37 m depth range. Six measurements spanning solar longitudes from 8.0$^\circ$ to 210.0$^\circ$ were made and atmospheric pressure at the site was simultaneously measured using InSight’s Pressure Sensor. We find that soil thermal conductivity strongly correlates with atmospheric pressure. This trend is compatible with predictions of the pressure dependence of thermal conductivity for unconsolidated soils under martian atmospheric conditions, indicating that heat transport through the pore filling gas is a major contributor to the total heat transport. This implies that any cementation or induration of the soil sampled by the experiments must be minimal and that the soil surrounding the mole at depths below the duricrust is unconsolidated. Thermal conductivity data presented here are the first direct evidence that the atmosphere interacts with the top most meter of material on Mars.

Observations of the South Polar Residual Cap suggest a possible erosion of the cap, leading to an increase of the global mass of the atmosphere. We test this assumption by making the first comparison between Viking 1 and InSight surface pressure data that have been recorded with ~40 years of difference. Such a comparison also allows us to determine changes in the dynamics of the seasonal ice caps between these two periods. To do so, we first had to recalibrate the InSight pressure data because of their unexpected sensitivity to the sensor temperature. Then, we had to design a procedure to compare distant pressure measurements. We propose two surface pressure interpolation methods at the local and global scale to do the comparison. The comparison of Viking and InSight seasonal surface pressure variations does not show major changes in the CO2 cycle. Such conclusions are also supported by an analysis of the Mars Science Laboratory (MSL) pressure data. Further comparisons with images of the south seasonal cap taken by the Viking 2 orbiter and MARCI camera do not display significant changes in the dynamic of this cap within ~40 years. Only a possible larger extension of the North Cap after the global storm of MY 34 is observed, but the physical mechanisms behind this anomaly are not well determined. Finally, the first comparison of MSL and InSight pressure data suggests a pressure deficit at Gale crater during southern summer, possibly resulting from a large presence of dust suspended within the crater.

The SEIS seismometer deployed at the surface of Mars in the framework of the NASA-InSight mission has been continuously recording the ground motion at Elysium Planitia for more than one martian year. In this work, we investigate the seasonal variation of the near surface properties using both background vibrations and a particular class of high-frequency seismic events. We present measurements of relative velocity changes over one martian year and show that they can be modeled by a thermoelastic response of the Martian regolith. Several families of high-frequency seismic multiplets have been observed at various periods of the martian year. These events exhibit repeatable waveforms with an emergent character and a coda that is likely composed of scattered waves. Taking advantage of these properties, we use coda waves interferometry to measure relative travel-time changes as a function of the date of occurrence of the quakes. While in some families a stretching of the coda waveform is clearly observed, in other families we observe either no variation or a clear contraction of the waveform. Measurements of velocity changes from the analysis of background vibrations above 5Hz are consistent with the results from coda wave interferometry. We identify a frequency band structure in the power spectral density, that can be tracked over hundreds of days. This band structure is the equivalent in the frequency domain of an autocorrelogram and can be efficiently used to measure relative travel-time changes as a function of frequency. The observed velocity changes can be adequately modeled by the thermoelastic response of the regolith to the time-dependent incident solar flux at the seasonal scale. In particular, the model captures the time delay between the surface temperature variations and the velocity changes in the sub-surface. Our observations could serve as a basis for a joint inversion of the seismic and thermal properties in the first meters below InSIght.

Before dawn on the dustiest regions of Mars, surfaces measured at or below ∼ 148 K are common. Thermodynamics principles indicate that these terrains must be associated with the presence of CO2 frost, yet visible wavelength imagery does not display any ice signature. We interpret this systematic absence as an indication of CO2 crystal growth within the surficial regolith, not on top of it, forming hard-to-distinguish intimate mixtures of frost and dust, i.e., dirty frost. This particular ice/regolith relationship unique to the low thermal inertia regions is enabled by the large difference in size between individual dust grains and the peak thermal emission wavelength of any material nearing 148 K (1-2 μm vs. 18 μm), allowing radiative loss (and therefore ice formation) to occur deep within the pores of the ground, below several layers of grains. After sunrise, sublimation-driven winds promoted by direct insolation and conduction create an upward drag within the surficial regolith that can be comparable in strength to gravity and friction forces combined. This drag displaces individual grains, possibly preventing their agglomeration, induration, and compaction, and can potentially initiate or sustain downslope mass movement such as slope streaks. If confirmed, this hypothesis introduces a new form of CO2-driven geomorphological activity occurring near the equator on Mars and explains how large units of mobile dust are currently maintained at the surface in an otherwise soil-encrusting world.