Search

In a previous work (Hern\'andez-Bernal et al. 2021) we performed an observational analysis of the Arsia Mons Elongated Cloud (AMEC), which stands out due to its impressive size and shape, quick dynamics, and the fact that it happens during the martian dusty season. Observations show that its morphology can be split in a head, on the western slope of the volcano of around 120 km in diameter; and a tail, that expands to the west reaching more than 1000 km in length, making the AMEC the longest orographic cloud observed so far in the solar system. In this work we run the LMD (Laboratoire de M\'et\'eorologie Dynamique) Mesoscale Model to gain insight into the physics of the AMEC. We note that it is coincident in terms of local time and seasonality with the fastest winds on the summit of Arsia Mons. A downslope windstorm on the western slope is followed by a hydraulic-like jump triggering a strong vertical updraft that propagates upwards in the atmosphere, causing a drop in temperatures of down to 30K at 40-50 km in altitude, spatially and temporarily coincident with the observed head of the AMEC. However the model does not reproduce the microphysics of this cloud: the optical depth is too low and the expansion of the tail does not happen in the model. The observed diurnal cycle is correctly captured by the model for the head of the cloud. This work raises new questions that will guide future observations of the AMEC.

Triton and Pluto are believed to share a common origin, both forming initially in the Kuiper Belt but Triton being later captured by Neptune. Both objects display similar sizes, densities, and atmospheric and surface ice composition, with the presence of volatile ices N2, CH4 and CO. Yet their appearance, including their surface albedo and ice distribution strongly differ. What can explain these different appearances? A first disparity is that Triton is experiencing significant tidal heating due to its orbit around Neptune, with subsequent resurfacing and a relatively flat surface, while Pluto is not tidally activated and displays a pronounced topography. Here we present long-term volatile transport simulations of Pluto and Triton, using the same initial conditions and volatile inventory, but with the known orbit and rotation of each object. The model reproduces, to first order, the observed volatile ice surface distribution on Pluto and Triton. Our results unambiguously demonstrate that obliquity is the main driver of the differences in surface appearance and in climate properties on Pluto and Triton, and give further support to the hypothesis that both objects had a common origin followed by a different dynamical history., Comment: 16 pages, 3 figures

Understanding the set of conditions that allow rocky planets to have liquid water on their surface -- in the form of lakes, seas or oceans -- is a major scientific step to determine the fraction of planets potentially suitable for the emergence and development of life as we know it on Earth. This effort is also necessary to define and refine the so-called "Habitable Zone" (HZ) in order to guide the search for exoplanets likely to harbor remotely detectable life forms. Until now, most numerical climate studies on this topic have focused on the conditions necessary to maintain oceans, but not to form them in the first place. Here we use the three-dimensional Generic Planetary Climate Model (PCM), historically known as the LMD Generic Global Climate Model (GCM), to simulate water-dominated planetary atmospheres around different types of Main-Sequence stars. The simulations are designed to reproduce the conditions of early ocean formation on rocky planets due to the condensation of the primordial water reservoir at the end of the magma ocean phase. We show that the incoming stellar radiation (ISR) required to form oceans by condensation is always drastically lower than that required to vaporize oceans. We introduce a Water Condensation Limit, which lies at significantly lower ISR than the inner edge of the HZ calculated with three-dimensional numerical climate simulations. This difference is due to a behavior change of water clouds, from low-altitude dayside convective clouds to high-altitude nightside stratospheric clouds. Finally, we calculated transit spectra, emission spectra and thermal phase curves of TRAPPIST-1b, c and d with H2O-rich atmospheres, and compared them to CO2 atmospheres and bare rock simulations. We show using these observables that JWST has the capability to probe steam atmospheres on low-mass planets, and could possibly test the existence of nightside water clouds., Comment: Accepted for publication in Astronomy & Astrophysics

This chapter of the Planetary Exploration Horizon 2061 Report reviews the way the six key questions about planetary systems, from their origins to the way they work and their habitability, identified in chapter 1, can be addressed by means of solar system exploration, and how one can find partial answers to these six questions by flying to the different provinces to the solar system: terrestrial planets, giant planets, small bodies, and up to its interface with the local interstellar medium. It derives from this analysis a synthetic description of the most important space observations to be performed at the different solar system objects by future planetary exploration missions. These observation requirements illustrate the diversity of measurement techniques to be used as well as the diversity of destinations where these observations must be made. They constitute the base for the identification of the future planetary missions we need to fly by 2061, which are described in chapter 4. Q1- How well do we understand the diversity of planetary systems objects? Q2- How well do we understand the diversity of planetary system architectures? Q3- What are the origins and formation scenarios for planetary systems? Q4- How do planetary systems work? Q5- Do planetary systems host potential habitats? Q6- Where and how to search for life?, Comment: 107 pages, 37 figures, Horizon 2061 is a science-driven, foresight exercise, for future scientific investigations

Temperature profiles retrieved using the first set of data of the Emirates Mars InfraRed Spectrometer (EMIRS) obtained during the science phase of the Emirates Mars Mission (EMM) are used for the analysis of migrating thermal tides in the Martian atmosphere. The selected data cover a solar longitude (LS) range of 60{\deg}-90{\deg} of Martian Year (MY) 36. The novel orbit design of the Hope Probe leads to a good geographic and local time coverage that significantly improves the analysis. Wave mode decomposition suggests dominant diurnal tide and important semi-diurnal tide with maximal amplitudes of 6K and 2K, respectively, as well as the existence of ~0.5K ter-diurnal tide. The results agree well with predictions by the Mars Planetary Climate Model (PCM), but the observed diurnal tide has an earlier phase (3h), and the semi-diurnal tide has an unexpectedly large wavelength (~200km)., Comment: 11 pages, 5 figures, accepted by GRL special issue "The First Results from the Emirates Mars Mission (EMM)"

The Deep Atmosphere Venus Investigation of Noble gases, Chemistry, and Imaging (DAVINCI) mission described herein has been selected for flight to Venus as part of the NASA Discovery Program. DAVINCI will be the first mission to Venus to incorporate science-driven flybys and an instrumented descent sphere into a unified architecture. The anticipated scientific outcome will be a new understanding of the atmosphere, surface, and evolutionary path of Venus as a possibly once-habitable planet and analog to hot terrestrial exoplanets. The primary mission design for DAVINCI as selected features a preferred launch in summer/fall 2029, two flybys in 2030, and descent sphere atmospheric entry by the end of 2031. The in situ atmospheric descent phase subsequently delivers definitive chemical and isotopic composition of the Venus atmosphere during a cloud-top to surface transect above Alpha Regio. These in situ investigations of the atmosphere and near infrared descent imaging of the surface will complement remote flyby observations of the dynamic atmosphere, cloud deck, and surface near infrared emissivity. The overall mission yield will be at least 60 Gbits (compressed) new data about the atmosphere and near surface, as well as first unique characterization of the deep atmosphere environment and chemistry, including trace gases, key stable isotopes, oxygen fugacity, constraints on local rock compositions, and topography of a tessera., Comment: 41 pages, 14 figures, accepted for publication in the Planetary Science Journal

Thermal tides in the Martian atmosphere are analyzed using temperature profiles retrieved from nadir observations obtained by the TIRVIM Fourier-spectrometer, part of the Atmospheric Chemistry Suite (ACS) onboard the ExoMars Trace Gas Orbiter (TGO). The data is selected near the northern summer solstice at solar longitude (LS) 75{\deg}-105{\deg} of Martian Year (MY) 35. The observations have a full local time coverage, which enables analyses of daily temperature anomalies. The observed zonal mean temperature is lower by 4-6K at ~100Pa, but higher towards the summer pole, compared to the LMD Mars General Circulation Model (GCM). Wave mode decomposition shows dominant diurnal tide and important semi-diurnal tide and diurnal Kelvin wave, with maximal amplitudes of 5K, 3K, and 2.5K, respectively, from tens to hundreds of Pa. The results generally agree well with the LMD Mars GCM, but with noticeable earlier phases of diurnal (~1h) and semi-diurnal (~3h) tides., Comment: 7 pages, 4 figures, accepted by GRL

The InSight mission, featuring continuous high-frequency high-sensitivity pressure measurements, is in ideal position to study the active atmospheric turbulence of Mars. Data acquired during 1.25 Martian year allows us to study the seasonal evolution of turbulence and its diurnal cycle. We investigate vortices (abrupt pressure drops), local turbulence (frequency range 0.01-2 Hz) and non-local turbulence often caused by convection cells and plumes (frequency range 0.002-0.01 Hz). Contrary to non-local turbulence, local turbulence is strongly sensitive at all local times and seasons to the ambient wind. We report many remarkable events with the arrival of northern autumn at the InSight landing site: a spectacular burst of daytime vortices, the appearance of nighttime vortices, and the development of nighttime local turbulence as intense as its daytime counterpart. Nighttime turbulence at this dusty season appears as a result of the combination of a stronger low-level jet, producing shear-driven turbulence, and a weaker stability., Comment: 15 pages, 5 figures, accepted version of the manuscript published in GRL

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) is a community project that aims to quantify how dfferences in general circulation models (GCMs) could impact the climate prediction for TRAPPIST-1e and, subsequently its atmospheric characterization in transit. Four GCMs have participated in THAI so far: ExoCAM, LMD-Generic, ROCKE-3D and the UM. This paper, focused on the simulated observations, is the third part of a trilogy, following the analysis of two land planet scenarios (part I) and two aquaplanet scenarios (part II). Here, we show a robust agreement between the simulated spectra and the number of transits estimated to detect the land planet atmospheres. For the aquaplanet ones, using atmospheric data from any of the four GCMs would require at least 17 transits. This prediction corresponds to UM simulated data which produces the lowest and thinnest clouds. Between 35-40% more clouds are predicted by ExoCAM or LMD-G due to higher thick terminator clouds. For the first time this work provides "GCM uncertainty error bars" of 35-40% that need to be considered in future analyses of transmission spectra. We also analyzed the inter-transit variability induced by weather patterns and changes of terminator cloudiness between transits. Its magnitude differs significantly between the GCMs but its impact on the transmission spectra is within the measurement uncertainties. THAI has demonstrated the importance of model intercomparison for exoplanets and also paved the way for a larger project to develop an intercomparison meta-framework, namely the Climates Using Interactive Suites of Intercomparisons Nested for Exoplanet Studies (CUISINES)., Comment: Accepted in the Planetary Science Journal as Part III of a series of 3 THAI papers

To identify promising exoplanets for atmospheric characterization and to make the best use of observational data, a thorough understanding of their atmospheres is needed. 3D general circulation models (GCMs) are one of the most comprehensive tools available for this task and will be used to interpret observations of temperate rocky exoplanets. Due to parameterization choices made in GCMs, they can produce different results, even for the same planet. Employing four widely-used exoplanetary GCMs -- ExoCAM, LMD-G, ROCKE-3D and the UM -- we continue the TRAPPIST-1 Habitable Atmosphere Intercomparison by modeling aquaplanet climates of TRAPPIST-1e with a moist atmosphere dominated by either nitrogen or carbon dioxide. Although the GCMs disagree on the details of the simulated regimes, they all predict a temperate climate with neither of the two cases pushed out of the habitable state. Nevertheless, the inter-model spread in the global mean surface temperature is non-negligible: 14 K and 24 K in the nitrogen and carbon dioxide dominated case, respectively. We find substantial inter-model differences in moist variables, with the smallest amount of clouds in LMD-Generic and the largest in ROCKE-3D. ExoCAM predicts the warmest climate for both cases and thus has the highest water vapor content and the largest amount and variability of cloud condensate. The UM tends to produce colder conditions, especially in the nitrogen-dominated case due to a strong negative cloud radiative effect on the day side of TRAPPIST-1e. Our study highlights various biases of GCMs and emphasizes the importance of not relying solely on one model to understand exoplanet climates., Comment: 35 pages, 22 figures; Published in The Planetary Science Journal

With the commissioning of powerful, new-generation telescopes such as the James Webb Space Telescope (JWST) and the ground-based Extremely Large Telescopes, the first characterization of a high molecular weight atmosphere around a temperate rocky exoplanet is imminent. Atmospheric simulations and synthetic observables of target exoplanets are essential to prepare and interpret these observations. Here we report the results of the first part of the TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) project, which compares 3D numerical simulations performed with four state-of-the-art global climate models (ExoCAM, LMD-Generic, ROCKE-3D, Unified Model) for the potentially habitable target TRAPPIST-1e. In this first part, we present the results of dry atmospheric simulations. These simulations serve as a benchmark to test how radiative transfer, subgrid-scale mixing (dry turbulence and convection), and large-scale dynamics impact the climate of TRAPPIST-1e and consequently the transit spectroscopy signature as seen by JWST. To first order, the four models give results in good agreement. The intermodel spread in the global mean surface temperature amounts to 7K (6K) for the N2-dominated (CO2-dominated) atmosphere. The radiative fluxes are also remarkably similar (intermodel variations less than 5%), from the surface (1 bar) up to atmospheric pressures around 5 mbar. Moderate differences between the models appear in the atmospheric circulation pattern (winds) and the (stratospheric) thermal structure. These differences arise between the models from (1) large-scale dynamics, because TRAPPIST-1e lies at the tipping point between two different circulation regimes (fast and Rhines rotators) in which the models can be alternatively trapped, and (2) parameterizations used in the upper atmosphere such as numerical damping., Comment: Published in the Planetary Science Journal as Part I of a series of 3 THAI papers

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison (THAI) workshop (14-16 September 2020). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements and limitations, as well as direction for future model development was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community., Comment: Accepted in the AAS Planetary Science Journal (PSJ)

CO2+H2 greenhouse warming has recently emerged as a promising scenario to sufficiently warm the early martian surface to allow the formation of valley networks and lakes. Here we present numerical 3-D global climate simulations of the early martian climate that we have performed assuming dense CO2+H2 atmospheres. Our climate model, derived from earlier works by Forget et al. (2013) and Wordsworth et al. (2013), is coupled to an asynchronous model of the long-term evolution of martian glaciers and lakes. Simulations were carried out at 40{\deg} obliquity to investigate how (i) water content and (ii) H2 content (added to 1 or 2 bars of CO2) can shape the climate and hydrologic cycle of early Mars. We show that the adiabatic cooling mechanism (Wordsworth et al. 2013) that leads to the accumulation of ice deposits in the southern highlands in cold climate (the so called 'icy highland scenario') also works in warm climates, with impact crater lakes acting as the main water reservoirs. This produces rainfall mainly localized in the southern highlands of Mars. If one adjust (i) the amount of CO2 and H2, (ii) the size and location of the water reservoirs, and (iii) the ancient topography (i.e. by removing Tharsis), the spatial patterns of surface runoff (from rainfall or snowmelt) in the simulations can match -- with a few exceptions -- the observed distribution of valley networks and impact crater lakes. Although our results are obtained for CO2-dominated atmospheres enriched with H2, they should also apply to assess the impact of any combination of powerful long-lived greenhouse gases on early Mars., Comment: Submitted for publication to Icarus. Comments and suggestions are welcome!

Pluto's atmospheric profiles (temperature and pressure) have been studied for decades from stellar occultation lightcurves. In this paper, we look at recent Pluto Global Climate Model (GCM) results (3D temperature, pressure, and density fields) from Bertrand et al. (2020) and use the results to generate model observer's plane intensity fields (OPIF) and lightcurves by using a Fourier optics scheme to model light passing through Pluto's atmosphere (Young, 2012). This approach can accommodate arbitrary atmospheric structures and 3D distributions of haze. We compared the GCM model lightcurves with the lightcurves observed during the 15-AUG-2018 Pluto stellar occultation. We find that the climate scenario which best reproduces the observed data includes an N2 ice mid-latitude band in the southern hemisphere. We have also studied different haze and P/T ratio profiles: the haze effectively reduces the central flash strength, and a lower P/T ratio both reduces the central flash strength and incurs anomalies in the shoulders of the central flash.

Exoplanetary science continues to excite and surprise with its rich diversity. We discuss here some key aspects potentially influencing the range of exoplanetary terrestrial-type atmospheres which could exist in nature. We are motivated by newly emerging observations, refined approaches to address data degeneracies, improved theories for key processes affecting atmospheric evolution and a new generation of atmospheric models which couple physical processes from the deep interior through to the exosphere and consider the planetary-star system as a whole. Using the Solar System as our guide we first summarize the main processes which sculpt atmospheric evolution then discuss their potential interactions in the context of exoplanetary environments. We summarize key uncertainties and consider a diverse range of atmospheric compositions discussing their potential occurrence in an exoplanetary context.

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., Comment: 22 pages, 7 figures, 1 table. Accepted for publication in Space Science Reviews. Part of ISSI special collection on Diversity of Atmospheres

Studying the atmospheric Planetary Boundary Layer (PBL) is crucial to understand the climate of a planet. The meteorological measurements by the instruments onboard InSight at a latitude of 4.5$^{\circ}$N 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 (from northern late winter to midsummer) to analyze wind gusts, convective cells and vortices in Mars' daytime PBL. We compare InSight measurements to turbulence-resolving Large-Eddy Simulations (LES). The daytime PBL turbulence at the InSight landing site is very active, with clearly identified signatures of convective cells and a vast population of 6000 recorded vortex encounters, adequately represented by a power-law with a 3.4 exponent. While the daily variability of vortex encounters at InSight can be explained by the statistical nature of turbulence, the seasonal variability is positively correlated with ambient wind speed, which is supported by LES. However, wind gustiness is positively correlated to surface temperature rather than ambient wind speed and sensible heat flux, 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 ambient wind speed. Typical tracks followed by vortices forming in the LES show a similar distribution in direction and length as orbital imagery., Comment: 44 pages, 19 figures, revised manuscript after peer-reviewed comments for consideration in Journal of Geophysical Research Planets (InSight special issue)

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., Comment: Accepted and published in Geoscientific Model Development (GMD)

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., Comment: 5 pages, 4 figures, accepted for publication in A&A Letters, V2 with editorial corrections

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., Comment: 34 pages, 17 figures

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., Comment: 30 pages, 20 figures, submitted to JGR Planets

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}., Comment: PDF, 55 pages, 22 figures, 2 tables. Accepted article (journal : Icarus)

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., Comment: Nature Scientific Reports, in press

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., Comment: Accepted for publication in Icarus. Abstract significantly abridged to meet ArXiv size limit

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., Comment: 25 pages, 1 table, and 10 figures

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., Comment: Accepted for publication in Icarus. 22 pages, 7 figures, 2 tables

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., Comment: 23 pages, 11 figures, in revision for publication in Astronomy & Astrophysics, abstract abridged for arXiv submission

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., Comment: 21 pages, 7 figures, 2 tables. Accepted for publication in Earth and Planetary Science Letters

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., Comment: 45 pages, 7 figures, accepted for publication in EPSL

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.