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.
Hot Jupiters (HJs) 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 heterogeneities between the day- and nightside of their atmosphere. For the hottest of these planets, 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 ultrahot Jupiters. Here, we quantify the effective temperature domain over which these biases are present. We used a set of 12 simulations of typical HJs 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 were then analyzed using the 1D TauREx retrieval code. We find that for James Webb Space Telescope like data, accounting for nonisothermal vertical temperature profiles is required over the whole temperature range. We further find that 1D retrieval codes start to estimate incorrect parameter values for planets with equilibrium temperatures greater than 1400 K if there are absorbers in the visible (such as TiO and VO, e.g.) that are able to create a hot stratosphere. The high temperatures at low pressures indeed entail a thermal dissociation of species that creates a strong chemical day-night dichotomy. As a byproduct, we demonstrate that when synthetic observations are used to assess the detectability of a given feature or process using a Bayesian framework (e.g., by comparing the Bayesian evidence of retrievals with different model assumptions), it is valid to use nonrandomized input data as long as the anticipated observational uncertainties are correctly taken into account.
Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal ...tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than ∼0.5 to 0.7 solar mass. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.
Most planets currently amenable to transit spectroscopy are close enough to their host stars to exhibit a relatively strong day to night temperature gradient. For hot planets this leads to 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 algorithms rely on 1D forward models that are unable to reproduce 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 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 TauREx retrieval code. We find that the ultra-hot Jupiter transmission spectra exhibit muted H
2
O features that originate on 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 the retrieval analysis to infer the C/O 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 Wide Field Camera 3 from
Hubble
Space Telescope transmission data of Wasp-121b are compatible with the day–night thermal and compositional dichotomy predicted by models.
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.
In an atmosphere, a cloud condensation region is characterized by a strong vertical gradient in the abundance of the related condensing species. On Earth, the ensuing gradient of mean molecular ...weight has relatively few dynamical consequences because N 2 is heavier than water vapor, so that only the release of latent heat significantly impacts convection. On the contrary, in a hydrogen dominated atmosphere (e.g., giant planets), all condensing species are significantly heavier than the background gas. This can stabilize the atmosphere against convection near a cloud deck if the enrichment in the given species exceeds a critical threshold. This raises two questions. What is transporting energy in such a stabilized layer, and how affected can the thermal profile of giant planets be? To answer these questions, we first carry out a linear analysis of the convective and double-diffusive instabilities in a condensable medium showing that an efficient condensation can suppress double-diffusive convection. This suggests that a stable radiative layer can form near a cloud condensation level, leading to an increase in the temperature of the deep adiabat. Then, we investigate the impact of the condensation of the most abundant species (water) with a steady-state atmosphere model. Compared to standard models, the temperature increase can reach several hundred degrees at the quenching depth of key chemical tracers. Overall, this effect could have many implications for our understanding of the dynamical and chemical state of the atmosphere of giant planets, for their future observations (with Juno for example), and for their internal evolution.
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 watts per square metre, 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.
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 H
2
O-dominated atmosphere. Here we use a 1-D radiative-convective inverse version of the LMD generic numerical climate model to derive new theoretical 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. As a result of the runaway greenhouse radius inflation effect introduced in previous work, these new mass-radius relationships significantly differ from those traditionally used in the literature. 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. Broadly speaking, our results demonstrate that non-H
2
/He-dominated atmospheres can have a first-order effect on the mass-radius relationships, even for rocky planets receiving moderate irradiation. Finally, we provide an empirical formula for the H
2
O steam atmosphere thickness as a function of planet core gravity and radius, water content, and irradiation. This formula can easily be used to construct mass-radius relationships for any water-rich, rocky planet (i.e., with any kind of interior composition ranging from pure iron to pure silicate) more irradiated than the runaway greenhouse irradiation threshold.
The increase in solar luminosity over geological timescales should warm the Earth's climate, increasing water evaporation, which will in turn enhance the atmospheric greenhouse effect. Above a ...certain critical insolation, this destabilizing greenhouse feedback can 'run away' until the oceans have completely evaporated. Through increases in stratospheric humidity, warming may also cause evaporative loss of the oceans to space before the runaway greenhouse state occurs. The critical insolation thresholds for these processes, however, remain uncertain because they have so far been evaluated using one-dimensional models that cannot account for the dynamical and cloud feedback effects that are key stabilizing features of the Earth's climate. Here we use a three-dimensional global climate model to show that the insolation threshold for the runaway greenhouse state to occur is about 375 W m(-2), which is significantly higher than previously thought. Our model is specifically developed to quantify the climate response of Earth-like planets to increased insolation in hot and extremely moist atmospheres. In contrast with previous studies, we find that clouds have a destabilizing feedback effect on the long-term warming. However, subsident, unsaturated regions created by the Hadley circulation have a stabilizing effect that is strong enough to shift the runaway greenhouse limit to higher values of insolation than are inferred from one-dimensional models. Furthermore, because of wavelength-dependent radiative effects, the stratosphere remains sufficiently cold and dry to hamper the escape of atmospheric water, even at large fluxes. This has strong implications for the possibility of liquid water existing on Venus early in its history, and extends the size of the habitable zone around other stars.