Search

Uranus and Neptune have atmospheres dominated by molecular hydrogen and helium. In the upper troposphere, methane is the third main molecule and condenses, yielding a vertical gradient in CH4. This condensable species being heavier than H2 and He, the resulting change in mean molecular weight due to condensation comes as a factor countering dry and moist convection. As observations also show latitudinal variations in methane abundance, one can expect different vertical gradients from one latitude to another. In this paper, we investigate the impact of this methane vertical gradient on the atmospheric regimes, especially on the formation and inhibition of moist convective storms in the troposphere of ice giants. We develop a 3D cloud-resolving model to simulate convective processes. Using our simulations, we conclude that typical velocities of dry convection in the deep atmosphere are rather low (of the order of 1 m/s) but sufficient to sustain upward methane transport, and that moist convection at methane condensation level is strongly inhibited. Previous studies derived an analytical criterion on the methane vapor amount above which moist convection should be inhibited. We first validate this analytical criterion numerically. We then show that the critical methane abundance governs the inhibition and formation of moist convective storms, and we conclude that the intensity and intermittency of these storms should depend on the methane abundance and saturation. In ice giants, dry convection is weak, and moist convection is strongly inhibited. However, when enough methane is transported upwards, through dry convection and turbulent diffusion, sporadic moist convective storms can form. These storms should be more frequent on Neptune than on Uranus, because of Neptune's internal heat flow. Our results can explain the observed sporadicity of clouds in ice giants.

With the new generation of space telescopes such as the James Webb Space Telescope (JWST), it is possible to better characterize the atmospheres of exoplanets. The atmospheres of Hot and Ultra Hot Jupiters are highly heterogeneous and asymmetrical. The difference between the temperatures on the day-side and the night-side is especially extreme in the case of Ultra Hot Jupiters. We introduce a new tool to compute synthetic lightcurves from 3D GCM simulations, developed in the Pytmosph3R framework. We show how rotation induces a variation of the flux during the transit that is a source of information on the chemical and thermal distribution of the atmosphere. We find that the day-night gradient linked to Ultra Hot Jupiters has an effect close to the stellar limb-darkening, but opposite to tidal deformation. We confirm the impact of the atmospheric and chemical distribution on variations of the central transit time, though the variations found are smaller than that of available observational data, which could indicate that the east-west asymmetries are underestimated, due to the chemistry or clouds., Comment: Accepted for publication in A&A (16th Feb.) 17 pages, 18 figures + 7 in appendices

Radiative transfer (RT) models are critical in the interpretation of exoplanetary spectra, in simulating exoplanet climates and when designing the specifications of future flagship observatories. However, most models differ in methodologies and input data, which can lead to significantly different spectra. In this paper, we present the experimental protocol of the MALBEC (Modeling Atmospheric Lines By the Exoplanet Community) project. MALBEC is an exoplanet model intercomparison project (exoMIP) that belongs to the CUISINES (Climates Using Interactive Suites of Intercomparisons Nested for Exoplanet Studies) framework which aims to provide the exoplanet community with a large and diverse set of comparison and validation of models. The proposed protocol tests include a large set of initial participating RT models, a broad range of atmospheres (from Hot Jupiters to temperate terrestrials) and several observation geometries, which would allow us to quantify and compare the differences between different RT models used by the exoplanetary community. Two types of tests are proposed: transit spectroscopy and direct imaging modeling, with results from the proposed tests to be published in dedicated follow-up papers. To encourage the community to join this comparison effort and as an example, we present simulation results for one specific transit case (GJ-1214 b), in which we find notable differences in how the various codes handle the discretization of the atmospheres (e.g., sub-layering), the treatment of molecular opacities (e.g., correlated-k, line-by-line) and the default spectroscopic repositories generally used by each model (e.g., HITRAN, HITEMP, ExoMol).

Hot Jupiters exhibit large day-night temperature contrasts. Their cooler nightsides are thought to host clouds. However, the exact nature of these clouds, their spatial distribution, and their impact on atmospheric dynamics, thermal structure, and spectra is still unclear. We investigate the atmosphere of WASP-43 b, a short period hot Jupiter recently observed with JWST, to understand the radiative and dynamical impact of clouds on the atmospheric circulation and thermal structure. We aim to understand the impact of different kinds of condensates potentially forming in WASP-43 b, with various sizes and atmospheric metallicities. We used a 3D global climate model (GCM) with a new temperature-dependent cloud model that includes radiative feed-backs coupled with hydrodynamical integrations to study the atmospheric properties of WASP-43 b. We produced observables from our simulations and compared them to spectral phase curves from various observations. We show that clouds have a net warming effect, meaning that the greenhouse effect caused by clouds is stronger than the albedo cooling effect. We show that the radiative effect of clouds has various impacts on the dynamical and thermal structure of WASP-43 b. Depending on the type of condensates and their sizes, the radiative-dynamical feedback will modify the horizontal and vertical temperature gradient and reduce the wind speed. For super-solar metallicity atmospheres, fewer clouds form in the atmosphere, leading to a weaker feedback. Comparisons with spectral phase curves observed with HST, Spitzer, and JWST indicate that WASP-43 b s nightside is cloudy and rule out sub-micron Mg2SiO4 cloud particles as the main opacity source. Distinguishing between cloudy solar and cloudy super-solar-metallicity atmospheres is not straightforward, and further observations of both reflected light and thermal emission are needed., Comment: 23 pages, 15 Figures + 7 Figures in appendix. Accepted for publication in Astronomy & Astrophysics

Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $\mu$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$\sigma$ upper limit of 1-6 parts per million, depending on model assumptions)., Comment: 61 pages, 13 figures, 4 tables. This preprint has been submitted to and accepted in principle for publication in Nature Astronomy without significant changes

While small, Neptune-like planets are among the most abundant exoplanets, our understanding of their atmospheric structure and dynamics remains sparse. In particular, many unknowns remain on the way moist convection works in these atmospheres where condensable species are heavier than the non-condensable background gas. While it has been predicted that moist convection could shut-down above some threshold abundance of these condensable species, this prediction is based on simple linear analysis and relies on strong assumptions on the saturation of the atmosphere. To investigate this issue, we develop a 3D cloud resolving model for H2 atmospheres with large amounts of condensable species and apply this model to a prototypical temperate Neptune-like planet -- K2-18b. Our model confirms the shut-down of moist convection and the onset of a stably stratified layer in the atmosphere, leading to much hotter deep atmospheres and interiors. Our 3D simulations further provide quantitative estimates of the turbulent mixing in this stable layer, which is a key driver of the cycling of condensables in the atmosphere. This allows us to build a very simple, yet realistic 1D model that captures the most salient features of the structure of Neptune-like atmospheres. Our qualitative findings on the behavior of moist convection in hydrogen atmospheres go beyond temperate planets and should also apply to the regions where iron and silicates condense in the deep interior of H2-dominated planets. Finally, we use our model to investigate the likelihood of a liquid ocean beneath a H2 dominated atmosphere on K2-18b. We find that the planet would need to have a very high albedo (>0.5-0.6) to sustain a liquid ocean. However, due to the spectral type of the star, the amount of aerosol scattering that would be needed to provide such a high albedo is inconsistent with the latest observational data., Comment: Accepted for publication in Astronomy and Astrophysics

Observed exoplanet transit spectra are usually retrieved using 1D models to determine atmospheric composition while planetary atmospheres are 3D. With the JWST and future space telescopes such as Ariel, we will be able to obtain increasingly accurate transit spectra. The 3D effects on the spectra will be visible, and we can expect biases in the 1D extractions. In order to elucidate these biases, we have built theoretical observations of transit spectra, from 3D atmospheric modeling through transit modeling to instrument modeling. 3D effects are observed to be strongly nonlinear from the coldest to the hottest planets. These effects also depend on the planet's metallicity and gravity. Considering equilibrium chemistry, 3D effects are observed through very strong variations in certain features of the molecule or very small variations over the whole spectrum. We conclude that we cannot rely on the uncertainty of retrievals at all pressures, and that we must be cautious about the results of retrievals at the top of the atmosphere. However the results are still fairly close to the truth at mid-altitudes (those probed). We also need to be careful with the chemical models used for planetary atmosphere. If the chemistry of one molecule is not correctly described, this will bias all the others, and the retrieved temperature as well. Finally, although fitting a wider wavelength range and higher resolution has been shown to increase retrieval accuracy, we show that this could depend on the wavelength range chosen, due to the accuracy on modeling the different features. In any case, 1D retrievals are still correct for the detection of molecules, even in the event of an erroneous abundance retrieval.

Water vapour atmospheres with content equivalent to the Earth's oceans, resulting from impacts or a high insolation, were found to yield a surface magma ocean. This was, however, a consequence of assuming a fully convective structure. Here we report, using a consistent climate model, that pure steam atmospheres are commonly shaped by radiative layers, making their thermal structure strongly dependent on the stellar spectrum and internal heat flow. The surface is cooler when an adiabatic profile is not imposed: melting Earth's crust requires an insolation several times higher than today, which will not happen during the main-sequence of the Sun. Venus' surface can solidify before the steam atmosphere escapes, which is opposite to previous works. Around the reddest stars ($T_{eff}<$3000K), surface magma oceans cannot form by stellar forcing alone, whatever the water content. These findings affect observable signatures of steam atmospheres and exoplanet mass-radius relationships, drastically changing current constraints on the water content of Trappist-1 planets. Unlike adiabatic structures, radiative-convective profiles are sensitive to opacities. New measurements of poorly constrained high-pressure opacities, in particular far from the H$_2$O absorption bands, are thus necessary to refine models of steam atmospheres, which are important stages in terrestrial planet evolution., Comment: 23 pages, 3 figures in the main text, 6 figures and 1 table in the Methods. Published in Nature: 09 August 2023 (received 14 July 2022, accepted 24 May 2023)

The recently discovered super-Earth LP 890-9 c is an intriguing target for atmospheric studies as it transits a nearby, low-activity late-type M-dwarf star at the inner edge of the Habitable Zone. Its position at the runaway greenhouse limit makes it a natural laboratory to study the climate evolution of hot rocky planets. We present the first 3D-GCM exo-Venus model for a modern Venus-like atmosphere (92 bar surface pressure, realistic composition, H$_2$SO$_4$ radiatively-active clouds), applied to the tidally-locked LP 890-9c to inform observations by JWST and future instruments. If LP 890-9 c has developed into a modern exo-Venus, then the modelled temperatures suggest that H$_2$SO$_4$ clouds are possible even in the substellar region. Like on modern Venus, clouds on LP 890-9 c would create a flat spectrum. The strongest CO$_2$ bands in transmission predicted by our model for LP 890-9 c are about 10 ppm, challenging detection, given JWST estimated noise floor. Estimated phase curve amplitudes are 0.9 and 2.4 ppm for continuum and CO$_2$ bands, respectively. While pointing out the challenge to characterise modern exo-Venus analogues, these results provide new insights for JWST proposals and highlight the influence of clouds in the spectrum of hot rocky exoplanet spectra., Comment: 6 pages, 4 figures, MNRAS letter

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

Neptune's tropospheric winds are among the most intense in the Solar System, but the dynamical mechanisms that produce them remain uncertain. Measuring wind speeds at different pressure levels may help understand the atmospheric dynamics of the planet. The goal of this work is to directly measure winds in Neptune's stratosphere with ALMA Doppler spectroscopy. We derived the Doppler lineshift maps of Neptune at the CO(3-2) and HCN(4-3) lines at 345.8 GHz ($\lambda$~0.87 mm) and 354.5 GHz (0.85 mm), respectively. For that, we used spectra obtained with ALMA in 2016 and recorded with a spatial resolution of ~0.37" on Neptune's 2.24" disk. After subtracting the planet solid rotation, we inferred the contribution of zonal winds to the measured Doppler lineshifts at the CO and HCN lines. We developed an MCMC-based retrieval methodology to constrain the latitudinal distribution of wind speeds. We find that CO(3-2) and HCN(4-3) lines probe the stratosphere of Neptune at pressures of $2^{+12}_{-1.8}$ mbar and $0.4^{+0.5}_{-0.3}$ mbar, respectively. The zonal winds at these altitudes are less intense than the tropospheric winds based on cloud tracking from Voyager observations. We find equatorial retrograde (westward) winds of $-180^{+70}_{-60}$ m/s for CO, and $-190^{+90}_{-70}$ m/s for HCN. Wind intensity decreases towards mid-latitudes, and wind speeds at 40$^\circ$S are $-90^{+50}_{-60}$ m/s for CO, and $-40^{+60}_{-80}$ m/s for HCN. Wind speeds become 0 m/s at about 50$^\circ$S, and we find that the circulation reverses to a prograde jet southwards of 60$^\circ$S. Overall, our direct stratospheric wind measurements match previous estimates from stellar occultation profiles and expectations based on thermal wind equilibrium. These are the first direct Doppler wind measurements performed on the Icy Giants, opening a new method to study and monitor their stratospheric dynamics., Comment: Accepted for publication in A&A Letters on 10/05/2023. 8 pages, 1 Table, 5 Figures

Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information content of the data resulted in high sensitivity to the varying assumptions made in the treatment of instrument systematics and the atmospheric retrieval analysis. Here we present a dayside thermal emission spectrum of the ultra-hot Jupiter WASP-18b obtained with the NIRISS instrument on JWST. The data span 0.85 to 2.85 $\mu$m in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at $>$6$\sigma$ confidence) and evidence for optical opacity, possibly due to H$^-$, TiO, and VO (combined significance of 3.8$\sigma$). Models that fit the data require a thermal inversion, molecular dissociation as predicted by chemical equilibrium, a solar heavy element abundance (''metallicity'', M/H = 1.03$_{-0.51}^{+1.11}$ $\times$ solar), and a carbon-to-oxygen (C/O) ratio less than unity. The data also yield a dayside brightness temperature map, which shows a peak in temperature near the sub-stellar point that decreases steeply and symmetrically with longitude toward the terminators., Comment: JWST ERS bright star observations. Uploaded to inform JWST Cycle 2 proposals. Manuscript under review. 50 pages, 14 figures, 2 tables

From the Archean toward the Proterozoic, the Earth's atmosphere underwent a major shift from anoxic to oxic conditions, around 2.4 to 2.1 Gyr, known as the Great Oxidation Event (GOE). This rapid transition may be related to an atmospheric instability caused by the formation of the ozone layer. Previous works were all based on 1D photochemical models. Here, we revisit the GOE with a 3D photochemical-climate model to investigate the possible impact of the atmospheric circulation and the coupling between the climate and the dynamics of the oxidation. We show that the diurnal, seasonal and transport variations do not bring significant changes compared to 1D models. Nevertheless, we highlight a temperature dependence for atmospheric photochemical losses. A cooling during the late Archean could then have favored the triggering of the oxygenation. In addition, we show that the Huronian glaciations, which took place during the GOE, could have introduced a fluctuation in the evolution of the oxygen level. Finally, we show that the oxygen overshoot which is expected to have occurred just after the GOE, was likely accompanied by a methane overshoot. Such high methane concentrations could have had climatic consequences and could have played a role in the dynamics of the Huronian glaciations.

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $\mu$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $\mu$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$\sigma$) and G395H (4.5$\sigma$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations., Comment: 39 pages, 14 figures, accepted to be published in Nature

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $\mu$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $\mu$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere., Comment: 35 pages, 13 figures, 3 tables, Nature, accepted

Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the medium-resolution (R$\sim$600) transmission spectrum of an exoplanet atmosphere between 3-5 $\mu$m covering multiple absorption features for the Saturn-mass exoplanet WASP-39b, obtained with JWST NIRSpec G395H. Our observations achieve 1.46x photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO$_2$ (28.5$\sigma$) and H$_2$O (21.5$\sigma$), and identify SO$_2$ as the source of absorption at 4.1 $\mu$m (4.8$\sigma$). Best-fit atmospheric models range between 3 and 10x solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO$_2$, underscore the importance of characterising the chemistry in exoplanet atmospheres, and showcase NIRSpec G395H as an excellent mode for time series observations over this critical wavelength range., Comment: 44 pages, 11 figures, 3 tables. Resubmitted after revision to Nature

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

Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 {\mu}m in wavelength and show a prominent CO2 absorption feature at 4.3 {\mu}m (26{\sigma} significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 {\mu}m that is not reproduced by these models., Comment: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

New-generation spectrographs dedicated to the study of exoplanetary atmospheres require a high accuracy in the atmospheric models to better interpret the input spectra. Thanks to space missions, the observed spectra will cover a large wavelength range from visible to mid-infrared with an higher precision compared to the old-generation instrumentation, revealing complex features coming from different regions of the atmosphere. For hot and ultra hot Jupiters (HJs and UHJs), the main source of complexity in the spectra comes from thermal and chemical differences between the day and the night sides. In this context, one-dimensional plane parallel retrieval models of atmospheres may not be suitable to extract the complexity of such spectra. In addition, Bayesian frameworks are computationally intensive and prevent us from using complete three-dimensional self-consistent models to retrieve exoplanetary atmospheres. We propose the TauREx 2D retrieval code, which uses two-dimensional atmospheric models as a good compromise between computational cost and model accuracy to better infer exoplanetary atmospheric characteristics for the hottest planets. TauREx 2D uses a 2D parametrization across the limb which computes the transmission spectrum from an exoplanetary atmosphere assuming azimuthal symmetry. It also includes a thermal dissociation model of various species. We demonstrate that, given an input observation, TauREx 2D mitigates the biases between the retrieved atmospheric parameters and the real atmospheric parameters. We also show that having a prior knowledge on the link between local temperature and composition is instrumental in inferring the temperature structure of the atmosphere. Finally, we apply such a model on a synthetic spectrum computed from a GCM simulation of WASP-121b and show how parameter biases can be removed when using two-dimensional forward models across the limb., Comment: 16 pages, 16 figures. Accepted for publication in Astronomy & Astrophysics

Considering the relatively high precision that will be reached by future observatories, it has recently become clear that one dimensional (1D) atmospheric models, in which the atmospheric temperature and composition of a planet are considered to vary only in the vertical, will be unable to represent exoplanetary transmission spectra with a sufficient accuracy. This is particularly true for warm to (ultra-) hot exoplanets because the atmosphere is unable to redistribute all the energy deposited on the dayside, creating a strong thermal and often compositional dichotomy on the planet. This situation is exacerbated by transmission spectroscopy, which probes the terminator region. This is the most heterogeneous region of the atmosphere. However, if being able to compute transmission spectra from 3D atmospheric structures (from a global climate model, e.g.) is necessary to predict realistic observables, it is too computationally expensive to be used in a data inversion framework. For this reason, there is a need for a medium-complexity 2D approach that captures the most salient features of the 3D model in a sufficiently fast implementation. With this in mind, we present a new open-source documented version of Pytmosph3R that handles the computation of transmission spectra for atmospheres with up to three spatial dimensions and can account for time variability. Taking the example of an ultra hot Jupiter, we illustrate how the changing orientation of the planet during the transit can allow us to probe the horizontal variations in the atmosphere. We further implement our algorithm in TauREx to allow the community to perform 2D retrievals. We describe our extensive cross-validation benchmarks and discuss the accuracy and numerical performance of each model., Comment: Accepted for publication in A&A. Pytmosph3R is available at http://perso.astrophy.u-bordeaux.fr/~jleconte/pytmosph3r-doc/index.html; Table 2 is corrected in this version

Hot Jupiters are very good targets for transmission spectroscopy analysis. Their atmospheres have a large scale height implying a high signal to noise ratio. As these planets orbit close to their stars, they often present strong thermal and chemical hetereogeneities between the day and the night side of their atmosphere. For the hottest ones, the thermal dissociation of several species occurs in their atmospheres which leads to a stronger chemical dichotomy between the two hemispheres. It has already been shown that the current retrieval algorithms, which are using 1D forward models, find biased molecular abundances in ultra hot Jupiters. Here, we quantify the effective temperature domain over which these biases are present. We use a set of 12 simulations of typical Hot Jupiters from Teq = 1000 K to Teq = 2100 K performed with the Substellar and Planetary Atmospheric Radiation and Circulation global climate model and generate transmission spectra that fully account for the 3D structure of the atmosphere with Pytmosph3R. These spectra are then analyzed using the 1D TauREx retrieval code. We find that for JWST-like data, accounting for non-isothermal vertical temperature profiles is required over the whole temperature range. We further find that 1D retrieval codes start to estimate wrong parameter values for planets with equilibrium temperatures greater than 1400 K if there are absorbers in the visible able to create a hot stratosphere. The high temperatures at low pressures indeed entail a thermal dissociation of species which creates a strong chemical day-night dichotomy. As a by-product, we demonstrate that when using synthetic observations to assess the detectability of a given feature or process using a Bayesian framework, it is valid to use non-randomized input data, as long as the anticipated observational uncertainties are correctly taken into account., Comment: 18 pages, 13 figures, 3 tables, 1 appendix. Received 3 August 2021; revised 17 August 2021; accepted 14 October 2021 in A&A

Earth has had oceans for nearly four billion years and Mars had lakes and rivers 3.5-3.8 billion years ago. However, it is still unknown whether water has ever condensed on the surface of Venus because the planet - now completely dry - has undergone global resurfacing events that obscure most of its history. The conditions required for water to have initially condensed on the surface of Solar System terrestrial planets are highly uncertain, as they have so far only been studied with one-dimensional numerical climate models that cannot account for the effects of atmospheric circulation and clouds, which are key climate stabilizers. Here we show using three-dimensional global climate model simulations of early Venus and Earth that water clouds - which preferentially form on the nightside, owing to the strong subsolar water vapour absorption - have a strong net warming effect that inhibits surface water condensation even at modest insolations (down to 325 W/m2, that is, 0.95 times the Earth solar constant). This shows that water never condensed and that, consequently, oceans never formed on the surface of Venus. Furthermore, this shows that the formation of Earth's oceans required much lower insolation than today, which was made possible by the faint young Sun. This also implies the existence of another stability state for present-day Earth: the 'Steam Earth', with all the water from the oceans evaporated into the atmosphere., Comment: Published in Nature (13 October 2021)

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution., Comment: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

We present a study on the spatially scanned spectroscopic observations of the transit of GJ 1132 b, a warm ($\sim$500 K) Super-Earth (1.13 R$_\oplus$) that was obtained with the G141 grism (1.125 - 1.650 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. We used the publicly available Iraclis pipeline to extract the planetary transmission spectra from the five visits and produce a precise transmission spectrum. We analysed the spectrum using the TauREx3 atmospheric retrieval code with which we show that the measurements do not contain molecular signatures in the investigated wavelength range and are best-fit with a flat-line model. Our results suggest that the planet does not have a clear primordial, hydrogen-dominated atmosphere. Instead, GJ 1132 b could have a cloudy hydrogen-dominated envelope, a very enriched secondary atmosphere, be airless, or have a tenuous atmosphere that has not been detected. Due to the narrow wavelength coverage of WFC3, these scenarios cannot be distinguished yet but the James Webb Space Telescope may be capable of detecting atmospheric features, although several observations may be required to provide useful constraints., Comment: 17 pages, 11 figures. Accepted for publication in AJ

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)

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.

Uranus and Neptune are the archetypes of "ice giants", a class of planets that may be among the most common in the Galaxy. They hold the keys to understand the atmospheric dynamics and structure of planets with hydrogen atmospheres inside and outside the solar system; however, they are also the last unexplored planets of the Solar System. Their atmospheres are active and storms are believed to be fueled by methane condensation which is both extremely abundant and occurs at low optical depth. This means that mapping temperature and methane abundance as a function of position and depth will inform us on how convection organizes in an atmosphere with no surface and condensates that are heavier than the surrounding air, a general feature of giant planets. Owing to the spatial and temporal variability of these atmospheres, an orbiter is required. A probe would provide a reference atmospheric profile to lift ambiguities inherent to remote observations. It would also measure the abundances of noble gases which can be used to reconstruct the history of planet formation in the Solar System. Finally, mapping the planets' gravity and magnetic fields will be essential to constrain their global composition, atmospheric dynamics, structure and evolution. An exploration of Uranus or Neptune will be essential to understand these planets and will also be key to constrain and analyze data obtained at Jupiter, Saturn, and for numerous exoplanets with hydrogen atmospheres., Comment: arXiv admin note: substantial text overlap with arXiv:1908.02092

With the major increase in the volume of the spectroscopic line lists needed to perform accurate radiative transfer calculations, disseminating accurate radiative data has become almost as much a challenge as computing it. Considering that many planetary science applications are only looking for heating rates or mid-to-low resolution spectra, any approach enabling such computations in an accurate and flexible way at a fraction of the computing and storage costs is highly valuable. For many of these reasons, the correlated-k approach has become very popular. Its major weakness has been the lack of ways to adapt the spectral grid/resolution of precomputed k-coefficients, making it difficult to distribute a generic database suited for many different applications. Currently, most users still need to have access to a line-by-line transfer code with the relevant line lists or high-resolution cross sections to compute k-coefficient tables at the desired resolution. In this work, we demonstrate that precomputed k-coefficients can be binned to a lower spectral resolution without any additional assumptions, and show how this can be done in practice. We then show that this binning procedure does not introduce any significant loss in accuracy. Along the way, we quantify how such an approach compares very favorably with the sampled cross section approach. This opens up a new avenue to deliver accurate radiative transfer data by providing mid-resolution k-coefficient tables to users who can later tailor those tables to their needs on the fly. To help with this final step, we briefly present Exo_k, an open-access, open-source Python library designed to handle, tailor, and use many different formats of k-coefficient and cross-section tables in an easy and computationally efficient way., Comment: In press at Astronomy and Astrophysics. The Exo_k library can be found at http://perso.astrophy.u-bordeaux.fr/~jleconte/exo_k-doc/index.html along with tutorials, example scripts and notebooks, and an extensive online documentation

Hubble (HST) spectroscopic transit observations of the temperate sub-Neptune K2-18b were interpreted as the presence of water vapour with potential water clouds. 1D modelling studies also predict the formation of water clouds at some conditions. However, such models cannot predict the cloud cover, driven by atmospheric dynamics and thermal contrasts, and thus their real impact on spectra. The main goal of this study is to understand the formation, distribution and observational consequences of water clouds on K2-18b and other temperate sub-Neptunes. We simulated the atmospheric dynamics, water cloud formation and spectra of K2-18b for H2-dominated atmosphere using a 3D GCM. We analysed the impact of atmospheric composition (with metallicity from 1*solar to 1000*solar), concentration of cloud condensation nuclei and planetary rotation rate. Assuming that K2-18b has a synchronous rotation, we show that the atmospheric circulation in the upper atmosphere essentially corresponds to a symmetric day-to-night circulation. This regime preferentially leads to cloud formation at the substellar point or at the terminator. Clouds form for metallicity >100*solar with relatively large particles. For 100-300*solar metallicity, the cloud fraction at the terminators is small with a limited impact on transit spectra. For 1000*solar metallicity, very thick clouds form at the terminator. The cloud distribution appears very sensitive to the concentration of CCN and to the planetary rotation rate. Fitting HST transit data with our simulated spectra suggests a metallicity of ~100-300*solar. In addition, we found that the cloud fraction at the terminator can be highly variable, leading to a potential variability in transit spectra. This effect could be common on cloudy exoplanets and could be detectable with multiple transit observations., Comment: Accepted for publication in A&A. 16 pages, 14 figures

We present an atmospheric characterization study of two medium sized planets bracketing the radius of Neptune: HD 106315 c (R$_{\rm{P}}$=4.98 $\pm$ 0.23 R$_{\oplus}$) and HD 3167 c (R$_{\rm{P}}$=2.740$_{-0.100}^{+0.106}$ R$_{\oplus}$). We analyse spatially scanned spectroscopic observations obtained with the G141 grism (1.125 - 1.650 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. We use the publicly available Iraclis pipeline and TauREx3 atmospheric retrieval code and we detect water vapor in the atmosphere of both planets with an abundance of $\log_{10}[\mathrm{H_2O}]=-2.1^{+0.7}_{-1.3}$ ($\sim$5.68$\sigma$) and $\log_{10}[\mathrm{H_2O}]=-4.1^{+0.9}_{-0.9}$ ($\sim$3.17$\sigma$) for HD 106315 c and HD 3167 c, respectively. The transmission spectrum of HD 106315 c shows also a possible evidence of ammonia absorption ($\log_{10}[\mathrm {NH_3}]=-4.3^{+0.7}_{-2.0}$, $\sim$1.97$\sigma$ -even if it is not significant-), whilst carbon dioxide absorption features may be present in the atmosphere of HD 3167 c in the $\sim$1.1-1.6~$\mu$m wavelength range ($\log_{10}[\mathrm{CO_{2}}]= -2.4^{+0.7}_{-1.0}$, $\sim$3.28$\sigma$). However the CO$_2$ detection appears significant, it must be considered carefully and put into perspective. Indeed, CO$_2$ presence is not explained by 1D equilibrium chemistry models, and it could be due to possible systematics. The additional contribution of clouds, CO and CH$_4$ are discussed. HD 106315 c and HD 3167 c will be interesting targets for upcoming telescopes such as the James Webb Space Telescope (JWST) and the Atmospheric Remote-Sensing Infrared Exoplanet Large-Survey (Ariel)., Comment: Accepted for publication in AJ

We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets' densities may be described with a single rocky mass-radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition, but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm/sec, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small; the orbits are extremely coplanar; and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations, and speculate on possible implications of the planet densities for the formation, migration and evolution of the planet system., Comment: Final version to be published in the Planetary Sciences Journal. 56 pages, 30 figures. Data from the paper and a complete table of forecast JWST times may be found at https://github.com/ericagol/TRAPPIST1_Spitzer/

TRAPPIST-1 is a fantastic nearby (~39.14 light years) planetary system made of at least seven transiting terrestrial-size, terrestrial-mass planets all receiving a moderate amount of irradiation. To date, this is the most observationally favourable system of potentially habitable planets. Since the announcement of the discovery of TRAPPIST-1 planets in 2016, a growing number of techniques and approaches have been used and proposed to reveal its true nature. Here we have compiled a state-of-the-art overview of all the observational and theoretical constraints that have been obtained so far using these techniques and approaches. The goal is to get a better understanding of whether or not TRAPPIST-1 planets can have atmospheres, and if so, what they are made of. For this, we surveyed the literature on TRAPPIST-1 about topics as broad as irradiation environment, orbital architecture, transit observations, density measurements, stellar contamination, and numerical climate and escape models. Each of these topics adds a brick to our understanding of the likely atmospheres of the seven planets. We show that (i) HST transit observations, (ii) density measurements, (iii) atmospheric escape modelling, and (iv) gas accretion modelling altogether offer solid evidence against the presence of H2-dominated atmospheres around TRAPPIST-1 planets. This means they likely have either (i) a high molecular weight atmosphere or (ii) no atmosphere at all. There are several key challenges ahead to characterize the bulk compositions of the atmospheres (if present) of TRAPPIST-1 planets. The main one so far is characterizing and correcting for the effects of stellar contamination. Fortunately, a new wave of observations with the James Webb Space Telescope and near-infrared high-resolution ground-based spectrographs on very large telescopes will bring significant advances in the coming decade., Comment: 55 pages, 11 figures, 2 tables. Review article accepted for publication in Space Science Reviews. Part of ISSI special collection on 'Understanding the Diversity of Planetary Atmospheres'. Abstract abridged to meet ArXiv size limit

We present the analysis of the hot-Jupiter KELT-7b using transmission and emission spectroscopy from the Hubble Space Telescope (HST), both taken with the Wide Field Camera 3 (WFC3). Our study uncovers a rich transmission spectrum which is consistent with a cloud-free atmosphere and suggests the presence of H2O and H-. In contrast, the extracted emission spectrum does not contain strong absorption features and, although it is not consistent with a simple blackbody, it can be explained by a varying temperature-pressure profile, collision induced absorption (CIA) and H-. KELT-7 b had also been studied with other space-based instruments and we explore the effects of introducing these additional datasets. Further observations with Hubble, or the next generation of space-based telescopes, are needed to allow for the optical opacity source in transmission to be confirmed and for molecular features to be disentangled in emission., Comment: 19 pages, 11 figures, accepted in AJ on June 23, 2020

This paper presents the atmospheric characterisation of three large, gaseous planets: WASP-127b, WASP-79b and WASP-62b. We analysed spectroscopic data obtained with the G141 grism (1.088 - 1.68 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) using the Iraclis pipeline and the TauREx3 retrieval code, both of which are publicly available. For WASP-127 b, which is the least dense planet discovered so far and is located in the short-period Neptune desert, our retrieval results found strong water absorption corresponding to an abundance of log(H$_2$O) = -2.71$^{+0.78}_{-1.05}$, and absorption compatible with an iron hydride abundance of log(FeH)=$-5.25^{+0.88}_{-1.10}$, with an extended cloudy atmosphere. We also detected water vapour in the atmospheres of WASP-79 b and WASP-62 b, with best-fit models indicating the presence of iron hydride, too. We used the Atmospheric Detectability Index (ADI) as well as Bayesian log evidence to quantify the strength of the detection and compared our results to the hot Jupiter population study by Tsiaras et al. 2018. While all the planets studied here are suitable targets for characterisation with upcoming facilities such as the James Webb Space Telescope (JWST) and Ariel, WASP-127 b is of particular interest due to its low density, and a thorough atmospheric study would develop our understanding of planet formation and migration., Comment: Accepted for publication in AJ

We analyse the transmission and emission spectra of the ultra-hot Jupiter WASP-76b, observed with the G141 grism of the Hubble Space Telescope's Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauRex 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-grey clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyse this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-grey cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterisation with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers., Comment: Accepted for publication in AJ

The distant ice giants of the Solar System, Uranus and Neptune, have only been visited by one space mission, Voyager 2. The current knowledge on their composition remains very limited despite some recent advances. A better characterization of their composition is however essential to constrain their formation and evolution, as a significant fraction of their mass is made of heavy elements, contrary to the gas giants Jupiter and Saturn. An in situ probe like Galileo would provide us with invaluable direct ground-truth composition measurements. However, some of the condensibles will remain out of the grasp of a shallow probe. While additional constraints could be obtained from a complementary orbiter, thermochemistry and diffusion modeling can further help us to increase the science return of an in situ probe., Comment: Accepted for publication in Space Science Reviews

Most planets currently amenable to transit spectroscopy are close enough to their host star to exhibit a relatively strong day to night temperature gradient. For hot planets, this leads to cause a chemical composition dichotomy between the two hemispheres. In the extreme case of ultra hot jupiters, some species, such as molecular hydrogen and water, are strongly dissociated on the day-side while others, such as carbon monoxide, are not. However, most current retrieval algorithm rely on 1D forward models that are unable to model this effect. We thus investigate how the 3D structure of the atmosphere biases the abundances retrieved using commonly used algorithms. We study the case of Wasp-121b as a prototypical ultra hot Jupiter. We use the simulations of this planet performed with the SPARC/MIT global climate model (GCM) and generate transmission spectra that fully account for the 3D structure of the atmosphere with Pytmopsh3R. These spectra are then analyzed using the \taurex retrieval code. We find that such ultra hot jupiter's transmission spectra exhibit muted H$_2$O features that originate in the night-side where the temperature, hence the scale-height, is smaller than on the day-side. However, the spectral features of molecules present on the day-side are boosted by both its high temperature and low mean molecular weight. As a result, the retrieved parameters are strongly biased compared to the ground truth. In particular the [CO]/[H$_2$O] is overestimated by one to three orders of magnitude. This must be kept in mind when using such retrieval analysis to infer the C/O ratio of a planet's atmosphere. We also discuss whether indicators can allow us to infer the 3D structure of an observed atmosphere. Finally we show that HST/WFC3 transmission data of Wasp-121b are compatible with the day-night thermal and compositional dichotomy predicted by models., Comment: 13 pages, 10 figures, accepted on March 12 2020 in Astronomy & Astrophysics, section 10. Planets and planetary systems

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

Mass-radius relationships for water-rich rocky planets are usually calculated assuming most water is present in condensed (either liquid or solid) form. Planet density estimates are then compared to these mass-radius relationships, even when these planets are more irradiated than the runaway greenhouse irradiation limit (around 1.1~times the insolation at Earth for planets orbiting a Sun-like star), for which water has been shown to be unstable in condensed form and would instead form a thick H2O-dominated atmosphere. Here we use the LMD Generic numerical climate model to derive new mass-radius relationships appropriate for water-rich rocky planets that are more irradiated than the runaway greenhouse irradiation limit, meaning planets endowed with a steam, water-dominated atmosphere. For a given water-to-rock mass ratio, these new mass-radius relationships lead to planet bulk densities much lower than calculated when water is assumed to be in condensed form. In other words, using traditional mass-radius relationships for planets that are more irradiated than the runaway greenhouse irradiation limit tends to dramatically overestimate -- possibly by several orders of magnitude -- their bulk water content. In particular, this result applies to TRAPPIST-1 b, c, and d, which can accommodate a water mass fraction of at most 2, 0.3 and 0.08 %, respectively, assuming planetary core with a terrestrial composition. In addition, we show that significant changes of mass-radius relationships (between planets less and more irradiated than the runaway greenhouse limit) can be used to remove bulk composition degeneracies in multiplanetary systems such as TRAPPIST-1. Finally, we provide an empirical formula for the H2O steam atmosphere thickness which can be used to construct mass-radius relationships for any water-rich, rocky planet more irradiated than the runaway greenhouse irradiation threshold., Comment: Final version published in A&A. (small typos in Tables 1 and 2 have been corrected)

Eccentricity tides generate a torque that can drive an ocean planet towards asynchronous rotation states of equilibrium when enhanced by resonances associated with the oceanic tidal modes. We investigate the impact of eccentricity tides on the rotation of rocky planets hosting a thin uniform ocean and orbiting cool dwarf stars such as TRAPPIST-1, with orbital periods ~1-10 days. Combining the linear theory of oceanic tides in the shallow water approximation with the Andrade model for the solid part of the planet, we develop a global model including the coupling effects of ocean loading, self-attraction, and deformation of the solid regions. We derive from this model analytic solutions for the tidal Love numbers and torque exerted on the planet. These solutions are used with realistic values of parameters provided by advanced models of the internal structure and tidal oscillations of solid bodies to explore the parameter space both analytically and numerically. Our model allows us to fully characterise the frequency-resonant tidal response of the planet, and particularly the features of resonances associated with the oceanic tidal modes (eigenfrequencies, resulting maxima of the tidal torque and Love numbers) as functions of the planet parameters (mass, radius, Andrade parameters, ocean depth and Rayleigh drag frequency). Resonances associated with the oceanic tide decrease the critical eccentricity beyond which asynchronous rotation states distinct from the usual spin-orbit resonances can exist. We provide an estimation and scaling laws for this critical eccentricity, which is found to be lowered by roughly one order of magnitude, switching from ~0.3 to ~0.06 in typical cases and to ~0.01 in extremal ones., Comment: Accepted for publication in Astronomy & Astrophysics, 28 pages, 9 figures

Thermal atmospheric tides have a strong impact on the rotation of terrestrial planets. They can lock these planets into an asynchronous rotation state of equilibrium. We aim at characterizing the dependence of the tidal torque resulting from the semidiurnal thermal tide on the tidal frequency, the planet orbital radius, and the atmospheric surface pressure. The tidal torque is computed from full 3D simulations of the atmospheric climate and mean flows using a generic version of the LMDZ general circulation model (GCM) in the case of a nitrogen-dominated atmosphere. Numerical results are discussed with the help of an updated linear analytical framework. Power scaling laws governing the evolution of the torque with the planet orbital radius and surface pressure are derived. The tidal torque exhibits i) a thermal peak in the vicinity of synchronization, ii) a resonant peak associated with the excitation of the Lamb mode in the high frequency range, and iii) well defined frequency slopes outside these resonances. These features are well explained by our linear theory. Whatever the star-planet distance and surface pressure, the torque frequency spectrum -- when rescaled with the relevant power laws -- always presents the same behaviour. This allows us to provide a single and easily usable empirical formula describing the atmospheric tidal torque over the whole parameter space. With such a formula, the effect of the atmospheric tidal torque can be implemented in evolutionary models of the rotational dynamics of a planet in a computationally efficient, and yet relatively accurate way., Comment: Accepted for publication in Astronomy & Astrophysics, 23 pages, 9 figures

Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus around the planet. Consequently, studies have focused on the effect of atmospheric horizontal heterogeneities along the limb. Here we demonstrate that the region probed in transmission -- the limb -- extends significantly toward the day and night sides of the planet. We show that the strong day-night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb to result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures. We first produce synthetic JWST observations from a simulation of GJ 1214b and demonstrate the necessity of a real 3D approach to model data for such precise observatories. Second, we investigate how day-night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesize a large set of forward spectra for prototypical HD 209458 b and GJ 1214 b type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrieval analyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside. These biases are difficult to detect as the 1D profiles used in the retrieval are found to provide an excellent match to the observed spectra. This needs to be kept in mind when interpreting real data., Comment: Accepted for publication at Astronomy and Astrophysics. Abstract abridged to meet ArXiv size limit

Hot Jupiters are submitted to an intense stellar heating. The resulting thermal tides can torque their atmospheres into asynchronous rotation, while these planets are usually assumed to be locked into spin-orbit synchronization with their host star. Particularly, the thermal atmospheric torque can be greatly enhanced by the dynamical component of the tidal response, that is the component associated with the propagation of internal waves. Owing to the involved complex dynamics, semi-analytical approaches are crucial to understand the physical mechanisms that are responsible for the frequency-resonant behavior of thermal tides, and quantify the atmospheric tidal torque. In this work, we revisit the early works by Arras \& Socrates (2010) and present an improved modeling of thermal tides taking into account rotation and radiative cooling. Using this new modeling, we compute analytically the atmospheric tidal response of hot Jupiters and show that resonances associated with low-frequency internal gravity waves are able to drive asynchronous zonal flows in the range 1-30 days., Comment: 5 pages, 3 figures, proceedings of the SF2A annual meeting

Oceanic tides are a major source of tidal dissipation. They are a key actor for the orbital and rotational evolution of planetary systems, and contribute to the heating of icy satellites hosting a subsurface ocean. Oceanic tides are characterized by a highly frequency-resonant behavior, which is mainly due to the propagation of surface gravity waves in the case of thin oceans, and internal waves when they are deeper. In this work, we derive self-consistent ab initio expressions of the oceanic tidal torque as a function of the key physical parameters of the system (the ocean depth, the Brunt-V\"ais\"al\"a stratification frequency, the rotation rate, the tidal frequency, the Rayleigh friction). These solutions include the coupled mechanisms of internal and surface gravito-inertial waves, which allows us to study the case of planets hosting deep oceans and offer interesting prospects for the coupling between subsurface oceans and ice shells in the case of icy satellites., Comment: 3 pages, 2 figures, proceedings of the SF2A annual meeting

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., Comment: 68 pages, 23 figures, published in Icarus https://doi.org/10.1016/j.icarus.2019.07.011

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

A large fraction of known rocky exoplanets are expected to have been spun-down to a state of synchronous rotation, including temperate ones. Studies about the atmospheric and surface processes occurring on such planets thus assume that the day/night sides are fixed with respect to the surface over geological timescales. Here we show that this should not be the case for many synchronous exoplanets. This is due to True Polar Wander (TPW), a well known process occurring on Earth and in the Solar System that can reorient a planet without changing the orientation of its rotational angular momentum with respect to an inertial reference frame. As on Earth, convection in the mantle of rocky exoplanets should continuously distort their inertia tensor, causing reorientation. Moreover, we show that this reorientation is made very efficient by the slower rotation rate of synchronous planets. This is due to the weakness of their combined rotational/tidal bulge---the main stabilizing factor limiting TPW. Stabilization by an elastic lithosphere is also shown to be inefficient. We thus expect the axes of smallest and largest moment of inertia to change continuously over time but to remain closely aligned with the star-planet and orbital axes, respectively., Comment: Published in Nature Geoscience. Accepted version of the manuscript before editing by the publisher

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission., Comment: PASP in press

Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting to study the formation and evolution of terrestrial planets that formed in the same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all known to better than 5% precision, their densities have significant uncertainties (between 28% and 95%) because of poor constraints on the planet's masses. Aims.The goal of this paper is to improve our knowledge of the TRAPPIST-1 planetary masses and densities using transit-timing variations (TTV). The complexity of the TTV inversion problem is known to be particularly acute in multi-planetary systems (convergence issues, degeneracies and size of the parameter space), especially for resonant chain systems such as TRAPPIST-1. Methods. To overcome these challenges, we have used a novel method that employs a genetic algorithm coupled to a full N-body integrator that we applied to a set of 284 individual transit timings. This approach enables us to efficiently explore the parameter space and to derive reliable masses and densities from TTVs for all seven planets. Results. Our new masses result in a five- to eight-fold improvement on the planetary density uncertainties, with precisions ranging from 5% to 12%. These updated values provide new insights into the bulk structure of the TRAPPIST-1 planets. We find that TRAPPIST-1\,c and e likely have largely rocky interiors, while planets b, d, f, g, and h require envelopes of volatiles in the form of thick atmospheres, oceans, or ice, in most cases with water mass fractions less than 5%., Comment: in press in Astronomy & Astrophysics, 24 pages, 13 figures