Search

Terrestrial exoplanets around M- and K-type stars are important targets for atmospheric characterisation. Such planets are likely tidally locked with the order of spin-orbit resonances (SORs) depending on eccentricity. We explore the impact of SORs on 3D atmospheric dynamics and chemistry, employing a 3D coupled Climate-Chemistry Model to simulate Proxima Centauri b in 1:1 and 3:2 SOR. For a 1:1 SOR, Proxima Centauri b is in the Rhines rotator circulation regime with dominant zonal gradients (global mean surface temperature 229 K). An eccentric 3:2 SOR warms Proxima Centauri b to 262 K with gradients in the meridional direction. We show how a complex interplay between stellar radiation, orbit, atmospheric circulation, and (photo)chemistry determines the 3D ozone distribution. Spatial variations in ozone column densities align with the temperature distribution and are driven by stratospheric circulation mechanisms. Proxima Centauri b in a 3:2 SOR demonstrates additional atmospheric variability, including daytime-nighttime cycles in water vapour of ${+}$55% to ${-}$34% and ozone ($\pm5.2$%) column densities and periastron-apoastron water vapour cycles of ${+}$17% to ${-}$10%. Synthetic emission spectra for the spectral range of the Large Interferometer For Exoplanets fluctuate by up to 36 ppm with orbital phase angle for a 1:1 SOR due to 3D spatial and temporal asymmetries. The homogeneous atmosphere for the 3:2 SOR results in relatively constant emission spectra and provides an observational discriminant from the 1:1 SOR. Our work emphasizes the importance of understanding the 3D nature of exoplanet atmospheres and associated spectral variations to determine habitability and interpret atmospheric spectra., Comment: 31 pages, 19 figures, revised and re-submitted to The Planetary Science Journal

Which rocky exoplanets have atmospheres, and which do not? This presumably simply question is the first that must be answered in the pathway toward understanding the prevalence of nearby habitable planets. A mere 6.9 pc from Earth, the LTT 1445A system is the closest transiting M-dwarf system, and its largest known planet, at $\rm 1.31\; R_{\oplus}$ and 424 K, is one of the best rocky targets in which to search for an atmosphere. We use HST/WFC3 transmission spectroscopy with the G280 and G141 grisms to study the spectrum of LTT 1445Ab between $\rm 0.2-1.65\;\mu m$. In doing so, we uncover a flare on the neighboring star LTT 1445C, and report one of the first simultaneous near-UV/optical spectra of an M dwarf flare. The planet spectrum is consistent with a flat line, though the infrared portion displays potential features that could be explained by known opacity sources such as HCN. Some atmospheric retrievals weakly favor ($\sim2\sigma$) an atmosphere over a flat line, but it remains challenging to discern between stellar contamination, an atmosphere, and a featureless spectrum at this time. We do, however, confidently rule out $\leq100\times$ solar metallicity atmospheres. Although stellar contamination retrievals cannot fit the infrared features well, the overall spectrum is consistent with stellar contamination from hot spots, cold spots, or both. Based on the UV/optical data, we place limits on the extent of stellar variability expected in the near-infrared ($30-40$ ppm), which will be critical for upcoming JWST observations of this highly optimal target., Comment: In Review at AJ. Submitted on 08/06/2024

Transmission spectroscopy has been a workhorse technique over the past two decades to constrain the physical and chemical properties of exoplanet atmospheres. One of its classical key assumptions is that the portion of the atmosphere it probes -- the terminator region -- is homogeneous. Several works in the past decade, however, have put this into question for highly irradiated, hot ($T_{eq}\gtrsim 1000$ K) gas giant exoplanets both empirically and via 3-dimensional modelling. While models predict clear differences between the evening (day-to-night) and morning (night-to-day) terminators, direct morning/evening transmission spectra in a wide wavelength range has not been reported for an exoplanet to date. Under the assumption of precise and accurate orbital parameters on WASP-39 b, here we report the detection of inhomogeneous terminators on the exoplanet WASP-39 b, which allows us to retrieve its morning and evening transmission spectra in the near-infrared ($2-5\ \mu$m) using JWST. We observe larger transit depths in the evening which are, on average, $405 \pm 88$ ppm larger than the morning ones, also having qualitatively larger features than the morning spectrum. The spectra are best explained by models in which the evening terminator is hotter than the morning terminator by $177^{+65}_{-57}$ K with both terminators having C/O ratios consistent with solar. General circulation models (GCMs) predict temperature differences broadly consistent with the above value and point towards a cloudy morning terminator and a clearer evening terminator., Comment: Published in Nature at https://www.nature.com/articles/s41586-024-07768-4. All code to produce plots (with data) can be found at https://github.com/nespinoza/wasp39-terminators

Planetary surface habitability has so far been, in the main, considered in its entirety. The increasing popularity of 3D modelling studies of (exo)planetary climate has highlighted the need for a more nuanced understanding of surface habitability. Using satellite-derived data of photosynthetic life to represent the observed surface habitability of modern Earth, we validate a set of climatologically-defined metrics previously used in exoplanetary habitability studies. The comparison finds that the metrics defined by temperature show spatial patterns of habitability distinct to those defined by aridity, with no metric able to completely replicate the observed. We build upon these results to introduce a new metric defined by the observed thermal limits of modern Earth-based life, along with surface water fluxes as an analogue for water and nutrient availability. Furthermore, we pay attention to not only the thermal bounds of macroscopic complex life, but additionally the limits of microbial life which have been vital to the generation of Earth's own biosignatures, thus expanding considerations of climatic habitability out of a historically binary definition. Repeating the validation for our metric finds a significant improvement in the spatial representation of habitability, laying the groundwork for more accurate assessments of potential life-supporting environments upon other planets., Comment: 26 pages, 2 figures

As JWST begins to return observations, it is more important than ever that exoplanet climate models can consistently and correctly predict the observability of exoplanets, retrieval of their data, and interpretation of planetary environments from that data. Model intercomparisons play a crucial role in this context, especially now when few data are available to validate model predictions. The CUISINES Working Group of NASA's Nexus for Exoplanet System Science (NExSS) supports a systematic approach to evaluating the performance of exoplanet models, and provides here a framework for conducting community-organized exoplanet Model Intercomparison Projects (exoMIPs). The CUISINES framework adapts Earth climate community practices specifically for the needs of exoplanet researchers, encompassing a range of model types, planetary targets, and parameter space studies. It is intended to help researchers to work collectively, equitably, and openly toward common goals. The CUISINES framework rests on five principles: 1) Define in advance what research question(s) the exoMIP is intended to address. 2) Create an experimental design that maximizes community participation, and advertise it widely. 3) Plan a project timeline that allows all exoMIP members to participate fully. 4) Generate data products from model output for direct comparison to observations. 5) Create a data management plan that is workable in the present and scalable for the future. Within the first years of its existence, CUISINES is already providing logistical support to 10 exoMIPs, and will continue to host annual workshops for further community feedback and presentation of new exoMIP ideas., Comment: 14 pages, 2 figures

Life on Earth has experienced numerous upheavals over its approximately 4 billion year history. In previous work we have discussed how interruptions to stability lead, on average, to increases in habitability over time, a tendency we called Entropic Gaia. Here we continue this exploration, working with the Tangled Nature Model of co-evolution, to understand how the evolutionary history of life is shaped by periods of acute environmental stress. We find that while these periods of stress pose a risk of complete extinction, they also create opportunities for evolutionary exploration which would otherwise be impossible, leading to more populous and stable states among the survivors than in alternative histories without a stress period. We also study how the duration, repetition and number of refugia into which life escapes during the perturbation affects the final outcome. The model results are discussed in relation to both Earth history and the search for alien life., Comment: 12 pages, 7 figures

In order to assess observational evidence for potential atmospheric biosignatures on exoplanets, it will be essential to test whether spectral fingerprints from multiple gases can be explained by abiotic or biotic-only processes. Here, we develop and apply a coupled 1D atmosphere-ocean-ecosystem model to understand how primitive biospheres, which exploit abiotic sources of H2, CO and O2, could influence the atmospheric composition of rocky terrestrial exoplanets. We apply this to the Earth at 3.8 Ga and to TRAPPIST-1e. We focus on metabolisms that evolved before the evolution of oxygenic photosynthesis, which consume H2 and CO and produce potentially detectable levels of CH4. O2-consuming metabolisms are also considered for TRAPPIST-1e, as abiotic O2 production is predicted on M-dwarf orbiting planets. We show that these biospheres can lead to high levels of surface O2 (approximately 1-5 %) as a result of \ch{CO} consumption, which could allow high O2 scenarios, by removing the main loss mechanisms of atomic oxygen. Increasing stratospheric temperatures, which increases atmospheric OH can reduce the likelihood of such a state forming. O2-consuming metabolisms could also lower O2 levels to around 10 ppm and support a productive biosphere at low reductant inputs. Using predicted transmission spectral features from CH4, CO, O2/O3 and CO2 across the hypothesis space for tectonic reductant input, we show that biotically-produced CH4 may only be detectable at high reductant inputs. CO is also likely to be a dominant feature in transmission spectra for planets orbiting M-dwarfs, which could reduce the confidence in any potential biosignature observations linked to these biospheres., Comment: 29 pages, 19 figures

Using a 3D General Circulation Model, the Unified Model, we present results from simulations of a tidally-locked TRAPPIST-1e with varying carbon dioxide CO2 and methane CH4 gas concentrations, and their corresponding prescribed spherical haze profiles. Our results show that the presence of CO2 leads to a warmer atmosphere globally due to its greenhouse effect, with the increase of surface temperature on the dayside surface reaching up to ~14.1 K, and on the nightside up to ~21.2 K. Increasing presence of CH4 first elevates the surface temperature on the dayside, followed by a decrease due to the balance of tropospheric warming and stratospheric cooling. A thin layer of haze, formed when the partial pressures of CH4 to CO2 (pCH4/pCO2) = 0.1, leads to a dayside warming of ~4.9K due to a change in the water vapour H2O distribution. The presence of a haze layer that formed beyond the ratio of 0.1 leads to dayside cooling. The haze reaches an optical threshold thickness when pCH4/pCO2 ~0.4 beyond which the dayside mean surface temperature does not vary much. The planet is more favourable to maintaining liquid water on the surface (mean surface temperature above 273.15 K) when pCO2 is high, pCH4 is low and the haze layer is thin. The effect of CO2, CH4 and haze on the dayside is similar to that for a rapidly-rotating planet. On the contrary, their effect on the nightside depends on the wind structure and the wind speed in the simulation., Comment: 17 pages, 16 figures, 2 tables

Convective processes are crucial in shaping exoplanetary atmospheres but are computationally expensive to simulate directly. A novel technique of simulating moist convection on tidally locked exoplanets is to use a global 3D model with a stretched mesh. This allows us to locally refine the model resolution to 4.7 km and resolve fine-scale convective processes without relying on parameterizations. We explore the impact of mesh stretching on the climate of a slowly rotating TRAPPIST-1e-like planet, assuming it is 1:1 tidally locked. In the stretched-mesh simulation with explicit convection, the climate is 5 K colder and 25% drier than that in the simulations with parameterized convection (with both stretched and quasi-uniform meshes)}. This is due to the increased cloud reflectivity - because of an increase of low-level cloudiness - and exacerbated by the diminished greenhouse effect due to less water vapor. At the same time, our stretched-mesh simulations reproduce the key characteristics of the global climate of tidally locked rocky exoplanets, without any noticeable numerical artifacts. Our methodology opens an exciting and computationally feasible avenue for improving our understanding of 3D mixing in exoplanetary atmospheres. Our study also demonstrates the feasibility of a global stretched mesh configuration for LFRic-Atmosphere, the next-generation Met Office climate and weather model., Comment: 15 pages, 7 figures; accepted for publication in The Astrophysical Journal

Transport-induced quenching in hot Jupiter atmospheres is a process that determines the boundary between the part of the atmosphere at chemical equilibrium and the part of the atmosphere at thermochemical (but not photothermochemical) disequilibrium. The location of this boundary, the quench level, depends on the interplay between the dynamical and chemical timescales in the atmosphere, with quenching occurring when these timescales are equal. We explore the sensitivity of the quench level position to an increase in the planet's atmospheric metallicity using aerosol-free 3D GCM simulations of a hot Jupiter WASP-96b. We find that the temperature increase at pressures of $\sim$$10^{4}-10^{7}$ Pa that occurs when metallicity is increased could shift the position of the quench level to pressures dominated by the jet, and cause an equatorial depletion of $CH_4$, $NH_3$ and $HCN$. We discuss how such a depletion affects the planet's transmission spectrum, and how the analysis of the evening-morning limb asymmetries, especially within $\sim3-5 {\mu}m$, could help distinguish atmospheres of different metallicities that are at chemical equilibrium from those with the upper layers at thermochemical disequilibrium., Comment: accepted to MNRAS, 27 pages, 24 figures

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 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

We demonstrate that LFRic-Atmosphere, a model built using the Met Office's GungHo dynamical core, is able to reproduce idealised large-scale atmospheric circulation patterns specified by several widely-used benchmark recipes. This is motivated by the rapid rate of exoplanet discovery and the ever-growing need for numerical modelling and characterisation of their atmospheres. Here we present LFRic-Atmosphere's results for the idealised tests imitating circulation regimes commonly used in the exoplanet modelling community. The benchmarks include three analytic forcing cases: the standard Held-Suarez test, the Menou-Rauscher Earth-like test, and the Merlis-Schneider Tidally Locked Earth test. Qualitatively, LFRic-Atmosphere agrees well with other numerical models and shows excellent conservation properties in terms of total mass, angular momentum and kinetic energy. We then use LFRic-Atmosphere with a more realistic representation of physical processes (radiation, subgrid-scale mixing, convection, clouds) by configuring it for the four TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) scenarios. This is the first application of LFRic-Atmosphere to a possible climate of a confirmed terrestrial exoplanet. LFRic-Atmosphere reproduces the THAI scenarios within the spread of the existing models across a range of key climatic variables. Our work shows that LFRic-Atmosphere performs well in the seven benchmark tests for terrestrial atmospheres, justifying its use in future exoplanet climate studies., Comment: 34 pages, 9(12) figures; Submitted to Geoscientific Model Development; Comments are welcome (see Discussion tab on the journal's website: https://egusphere.copernicus.org/preprints/2023/egusphere-2023-647)

Determining the habitability and interpreting future atmospheric observations of exoplanets requires understanding the atmospheric dynamics and chemistry from a 3-D perspective. Previous studies have shown significant spatial variability in the ozone layer of synchronously rotating M-dwarf planets, assuming an Earth-like initial atmospheric composition. We use a 3-D Coupled Climate-Chemistry model to understand this distribution of ozone and identify the mechanism responsible for it. We document a previously unreported connection between the ozone production regions on the photochemically active dayside hemisphere and the nightside devoid of stellar radiation and thus photochemistry. We find that stratospheric dayside-to-nightside overturning circulation can advect ozone-rich air to the nightside. On the nightside, ozone-rich air subsides at the locations of two quasi-stationary Rossby gyres, resulting in an exchange between the stratosphere and troposphere and the accumulation of ozone at the gyre locations. We identify the hemispheric contrast in radiative heating and cooling as the main driver of this ozone circulation. Dynamically-driven chemistry also impacts other tracer species in the atmosphere (gaseous and non-gaseous phase) as long as chemical lifetimes exceed dynamical lifetimes. These findings illustrate the 3-D nature of planetary atmospheres, predicting spatial and temporal variability that will impact spectroscopic observations of exoplanet atmospheres., Comment: 13 pages, 12 figures, submitted to MNRAS

Methane is thought to have been an important greenhouse gas during the Archean, although its potential warming has been found to be limited at high concentrations due to its high shortwave absorption. We use the Met Office Unified Model, a general circulation model, to further explore the climatic effect of different Archean methane concentrations. Surface warming peaks at a pressure ratio CH$_4$:CO$_2$ of approximately 0.1, reaching a maximum of up to 7 K before significant cooling above this ratio. Equator-to-pole temperature differences also tend to increase up to pCH$_4$ $\leq$300 Pa, which is driven by a difference in radiative forcing at the equator and poles by methane and a reduction in the latitudinal extend of the Hadley circulation. 3D models are important to fully capture the cooling effect of methane, due to these impacts of the circulation., Comment: 36 pages, 18 figures

We present results from the Met Office Unified Model (UM), a world-leading climate and weather model, adapted to simulate a dry Martian climate. We detail the adaptation of the basic parameterisations and analyse results from two simulations, one with radiatively active mineral dust and one with radiatively inactive dust. These simulations demonstrate how the radiative effects of dust act to accelerate the winds and create a mid-altitude isothermal layer during the dusty season. We validate our model through comparison with an established Mars model, the Laboratoire de M\'et\'eorologie Dynamique planetary climate model (PCM), finding good agreement in the seasonal wind and temperature profiles but with discrepancies in the predicted dust mass mixing ratio and conditions at the poles. This study validates the use of the UM for a Martian atmosphere, highlighting how the adaptation of an Earth general circulation model (GCM) can be beneficial for existing Mars GCMs and provides insight into the next steps in our development of a new Mars climate model., Comment: 20 pages, 12 figures with additional appendix figures

We present the first exoplanet phase curve measurement made with the JWST NIRSpec instrument, highlighting the exceptional stability of this newly-commissioned observatory for exoplanet climate studies. The target, WASP-121b, is an ultrahot Jupiter with an orbital period of 30.6 hr. We analyze two broadband light curves generated for the NRS1 and NRS2 detectors, covering wavelength ranges of 2.70-3.72 micron and 3.82-5.15 micron, respectively. Both light curves exhibit minimal systematics, with approximately linear drifts in the baseline flux level of 30 ppm/hr (NRS1) and 10 ppm/hr (NRS2). Assuming a simple brightness map for the planet described by a low-order spherical harmonic dipole, our light curve fits suggest that the phase curve peaks coincide with orbital phases $3.36 \pm 0.11$ deg (NRS1) and $2.66 \pm 0.12$ deg (NRS2) prior to mid-eclipse. This is consistent with the strongest dayside emission emanating from eastward of the substellar point. We measure planet-to-star emission ratios of $3,924 \pm 7$ ppm (NRS1) and $4,924 \pm 9$ ppm (NRS2) for the dayside hemisphere, and $136 \pm 8$ ppm (NRS1) and $630 \pm 10$ ppm (NRS2) for the nightside hemisphere. The latter nightside emission ratios translate to planetary brightness temperatures of $926 \pm 12$ K (NRS1) and $1,122 \pm 10$ K (NRS2), which are low enough for a wide range of refractory condensates to form, including enstatite and forsterite. A nightside cloud deck may be blocking emission from deeper, hotter layers of the atmosphere, potentially helping to explain why cloud-free 3D general circulation model simulations systematically over-predict the nightside emission for WASP-121b., Comment: Accepted for publication in Astrophysical Journal Letters on December 29, 2022

Cloud cover at the planetary limb of water-rich Earth-like planets is likely to weaken chemical signatures in transmission spectra, impeding attempts to characterize these atmospheres. However, based on observations of Earth and solar system worlds, exoplanets with atmospheres should have both short-term weather and long-term climate variability, implying that cloud cover may be less during some observing periods. We identify and describe a mechanism driving periodic clear sky events at the terminators in simulations of tidally locked Earth-like planets. A feedback between dayside cloud radiative effects, incoming stellar radiation and heating, and the dynamical state of the atmosphere, especially the zonal wavenumber-1 Rossby wave identified in past work on tidally locked planets, leads to oscillations in Rossby wave phase speeds and in the position of Rossby gyres and results in advection of clouds to or away from the planet's eastern terminator. We study this oscillation in simulations of Proxima Centauri b, TRAPPIST 1-e, and rapidly rotating versions of these worlds located at the extreme inner edge of their stars' habitable zones. We simulate time series of the transit depths of the 1.4 {\mu}m water feature and 2.7 {\mu}m carbon dioxide feature. The impact of atmospheric variability on the transmission spectra is sensitive to the structure of the dayside cloud cover and the location of the Rossby gyres, but none of our simulations have variability significant enough to be detectable with current methods., Comment: 21 pages, 11 figures

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

Transport-induced quenching, i.e., the homogenisation of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical modelling of this process exists, the three-dimensional nature of transport-induced chemistry is underexplored. Here we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b, HD 189733b, HD 209458b, and WASP-17b including coupled hydrodynamics, radiative transfer and chemistry. Our simulations were performed with two chemical schemes: a chemical kinetics scheme, which is capable of capturing transport-induced quenching, and a simpler, more widely used chemical equilibrium scheme. We find that transport-induced quenching is predicted to occur in atmospheres of all planets in our sample; however, the extent to which it affects their synthetic spectra and phase curves varies from planet to planet. This implies that there is a "sweet spot" for the observability of signatures of transport-induced quenching, which is controlled by the interplay between the dynamics and chemistry., Comment: accepted to MNRAS, 27 pages, 15 figures, with DOI for the data repository

With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors., Comment: Accepted to PSJ

Stellar flares present challenges to the potential habitability of terrestrial planets orbiting M dwarf stars through inducing changes in the atmospheric composition and irradiating the planet's surface in large amounts of ultraviolet light. To examine their impact, we have coupled a general circulation model with a photochemical kinetics scheme to examine the response and changes of an Earth-like atmosphere to stellar flares and coronal mass ejections. We find that stellar flares increase the amount of ozone in the atmosphere by a factor of 20 compared to a quiescent star. We find that coronal mass ejections abiotically generate significant levels of potential bio-signatures such as N$_2$O. The changes in atmospheric composition cause a moderate decrease in the amount of ultraviolet light that reaches the planets surface, suggesting that while flares are potentially harmful to life, the changes in the atmosphere due to a stellar flare act to reduce the impact of the next stellar flare., Comment: 27 pages, 21 figures, accepted to MNRAS

Determining the habitability and interpreting atmospheric spectra of exoplanets requires understanding their atmospheric physics and chemistry. We use a 3-D Coupled Climate-Chemistry Model, the Met Office Unified Model with the UK Chemistry and Aerosols framework, to study the emergence of lightning and its chemical impact on tidally-locked Earth-like exoplanets. We simulate the atmosphere of Proxima Centauri b orbiting in the Habitable Zone of its M-dwarf star, but the results apply to similar M-dwarf orbiting planets. Our chemical network includes the Chapman ozone reactions and hydrogen oxide (HO$_{\mathrm{x}}$=H+OH+HO$_2$) and nitrogen oxide (NO$_{\mathrm{x}}$=NO+NO$_2$) catalytic cycles. We find that photochemistry driven by stellar radiation (177-850 nm) supports a global ozone layer between 20-50 km. We parameterise lightning flashes as a function of cloud-top height and the resulting production of nitric oxide (NO) from the thermal decomposition of N$_2$ and O$_2$. Rapid dayside convection over and around the substellar point results in lightning flash rates of up to 0.16 flashes km$^{-2}$yr$^{-1}$, enriching the dayside atmosphere below altitudes of 20 km in NO$_{\mathrm{x}}$. Changes in dayside ozone are determined mainly by UV irradiance and the HO$_{\mathrm{x}}$ catalytic cycle. ~45% of the planetary dayside surface remains at habitable temperatures (T$_{\mathrm{surf}}$>273.15 K) and the ozone layer reduces surface UV radiation levels to 15%. Dayside-nightside thermal gradients result in strong winds that subsequently advect NO$_{\mathrm{x}}$ towards the nightside, where the absence of photochemistry allows NO$_{\mathrm{x}}$ chemistry to involve reservoir species. Our study also emphasizes the need for accurate UV stellar spectra to understand the atmospheric chemistry of exoplanets., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

Using a 3D general circulation model, we demonstrate that a confirmed rocky exoplanet and a primary observational target, TRAPPIST-1e presents an interesting case of climate bistability. We find that the atmospheric circulation on TRAPPIST-1e can exist in two distinct regimes for a 1~bar nitrogen-dominated atmosphere. One is characterized by a single strong equatorial prograde jet and a large day-night temperature difference; the other is characterized by a pair of mid-latitude prograde jets and a relatively small day-night contrast. The circulation regime appears to be highly sensitive to the model setup, including initial and surface boundary conditions, as well as physical parameterizations of convection and cloud radiative effects. We focus on the emergence of the atmospheric circulation during the early stages of simulations and show that the regime bistability is associated with a delicate balance between the zonally asymmetric heating, mean overturning circulation, and mid-latitude baroclinic instability. The relative strength of these processes places the GCM simulations on different branches of the evolution of atmospheric dynamics. The resulting steady states of the two regimes have consistent differences in the amount of water content and clouds, affecting the water absorption bands as well as the continuum level in the transmission spectrum, although they are too small to be detected with current technology. Nevertheless, this regime bistability affects the surface temperature, especially on the night side of the planet, and presents an interesting case for understanding atmospheric dynamics and highlights uncertainty in 3D GCM results, motivating more multi-model studies., Comment: 31 pages, 14 figures, accepted to the Planetary Science Journal

We present new analysis of infrared transmission spectroscopy of the cloud-free hot-Saturn WASP-96b performed with the Hubble and Spitzer Space Telescopes (HST and Spitzer). The WASP-96b spectrum exhibits the absorption feature from water in excellent agreement with synthetic spectra computed assuming a cloud-free atmosphere. The HST-Spitzer spectrum is coupled with Very Large Telescope (VLT) optical transmission spectroscopy which reveals the full pressure-broadened profile of the sodium absorption feature and enables the derivation of absolute abundances. We confirm and correct for a spectral offset of $\Delta R_{{\rm p}}/R_{\ast}=(-4.29^{+0.31}_{-0.37})\,\times10^{-3}$ of the VLT data relative to the HST-Spitzer spectrum. This offset can be explained by the assumed radius for the common-mode correction of the VLT spectra, which is a well-known feature of ground-based transmission spectroscopy. We find evidence for a lack of chromospheric and photometric activity of the host star which, therefore, make a negligible contribution to the offset. We measure abundances for Na and O that are consistent with solar to supersolar, with abundances relative to solar values of $21^{+27}_{-14}$ and $7^{+11}_{-4}$, respectively. We complement the transmission spectrum with new thermal emission constraints from Spitzer observations at 3.6 and $4.5\mu$m, which are best explained by the spectrum of an atmosphere with a temperature decreasing with altitude. A fit to the spectrum assuming an isothermal blackbody atmosphere constrains the dayside temperature to be $T_{\rm{p}}$=$1545$$\pm$$90$K., Comment: Accepted for publication in MNRAS

The temperature profile of a planetary atmosphere is a key diagnostic of radiative and dynamical processes governing the absorption, redistribution, and emission of energy. Observations have revealed dayside stratospheres that either cool or warm with altitude for a small number of gas giant exoplanets, while other dayside stratospheres are consistent with constant temperatures. Here we report spectroscopic phase curve measurements for the gas giant WASP-121b, which constrain stratospheric temperatures throughout the diurnal cycle. Variations measured for a water vapour spectral feature reveal a temperature profile that transitions from warming with altitude on the dayside hemisphere to cooling with altitude on the nightside hemisphere. The data are well explained by models assuming chemical equilibrium, with water molecules thermally dissociating at low pressures on the dayside and recombining on the nightside. Nightside temperatures are low enough for perovskite (CaTiO3) to condense, which could deplete titanium from the gas phase and explain recent non-detections at the day-night terminator. Nightside temperatures are also consistent with the condensation of refractory species such as magnesium, iron, and vanadium. Detections of these metals at the day-night terminator suggest, however, that if they do form nightside clouds, cold trapping does not efficiently remove them from the upper atmosphere. Horizontal winds and vertical mixing could keep these refractory condensates aloft in the upper atmosphere of the nightside hemisphere until they are recirculated to the hotter dayside hemisphere and vaporised., Comment: Published in Nature Astronomy (publisher version is Open Access)

Thermal emission has now been observed from many dozens of exoplanet atmospheres, opening the gateway to population-level characterization. Here, we provide theoretical explanations for observed trends in $\textit{Spitzer}$ IRAC channel 1 (3.6 $\mu m$) and channel 2 (4.5 $\mu m$) photometric eclipse depths (EDs) across a population of 34 hot Jupiters. We apply planet-specific, self-consistent atmospheric models, spanning a range of recirculation factors, metallicities, and C/O ratios, to probe the information content of $\textit{Spitzer}$ secondary eclipse observations across the hot-Jupiter population. We show that most hot Jupiters are inconsistent with blackbodies from $\textit{Spitzer}$ observations alone. We demonstrate that the majority of hot Jupiters are consistent with low energy redistribution between the dayside and nightside (hotter dayside than expected with efficient recirculation). We also see that high equilibrium temperature planets (T$_{eq}$ $\ge$ 1800 K) favor inefficient recirculation in comparison to the low temperature planets. Our planet-specific models do not reveal any definitive population trends in metallicity and C/O ratio with current data precision, but more than 59 % of our sample size is consistent with the C/O ratio $\leq$ 1 and 35 % are consistent with whole range (0.35 $\leq$ C/O $\leq$ 1.5). We also find that for most of the planets in our sample, 3.6 and 4.5 $\mu m$ model EDs lie within $\pm$1 $\sigma$ of the observed EDs. Intriguingly, few hot Jupiters exhibit greater thermal emission than predicted by the hottest atmospheric models (lowest recirculation) in our grid. Future spectroscopic observations of thermal emission from hot Jupiters with the James Webb Space Telescope will be necessary to robustly identify population trends in chemical compositions with its increased spectral resolution, range and data precision., Comment: 24 pages and 10 Figures. Published in the Astrophysical Journal (ApJ). Supplementary figures available in the online version of the journal

Using a three-dimensional general circulation model, we show that the atmospheric dynamics on a tidally locked Earth-like exoplanet, simulated with the planetary and orbital parameters of Proxima Centauri b, support a longitudinally asymmetric stratospheric wind oscillation (LASO), analogous to Earth's quasi-biennial oscillation (QBO). In our simulations, the LASO has a vertical extent of 35--55 km, a period of 5--6.5 months, and a peak-to-peak wind speed amplitude of -70 to +130 m/s with a maximum at an altitude of 41 km. Unlike the QBO, the LASO displays longitudinal asymmetries related to the asymmetric thermal forcing of the planet and to interactions with the resulting stationary Rossby waves. The equatorial gravity wave sources driving the LASO are localised in the deep convection region at the substellar point and in a jet exit region near the western terminator, unlike the QBO, for which these sources are distributed uniformly around the planet. Longitudinally, the western terminator experiences the highest wind speeds and undergoes reversals earlier than other longitudes. The antistellar point only experiences a weak oscillation with a very brief, low-speed westward phase. The QBO on Earth is associated with fluctuations in the abundances of water vapour and trace gases such as ozone which are also likely to occur on exoplanets if these gases are present. Strong fluctuations in temperature and the abundances of atmospheric species at the terminators will need to be considered when interpreting atmospheric observations of tidally locked exoplanets., Comment: 18 pages, 7 figures

Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for general use in planetary atmosphere RT problems. We aim to automate the post-processing pipeline, starting from direct global circulation model (GCM) output to synthetic spectra. We develop albedo, emission and transmission spectra modes for 3D and 1D input structures. We include capability to use correlated-k and high-resolution opacity tables, the latter of which can be Doppler shifted inside the model. We post-process results from several GCM groups including ExoRad, SPARC/MITgcm THOR, UK Met Office UM, Exo-FMS and the Rauscher model. Users can therefore take advantage of desktop and HPC GPU computing solutions. gCMCRT is well suited for post-processing large GCM model grids produced by members of the community and for high-resolution 3D investigations., Comment: 20 pages, 9 figures, submitted to ApJ 29/10/21, revised post-reviewer comments 18/3/22

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

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

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

We present ground-based optical transmission spectroscopy of the low-density hot Jupiter WASP-88b covering the wavelength range 4413-8333 {\AA} with the FORS2 spectrograph on the Very Large Telescope. The FORS2 white light curves exhibit a significant time-correlated noise which we model using a Gaussian Process and remove as a wavelength-independent component from the spectroscopic light curves. We analyse complementary photometric observations from the Transiting Exoplanet Survey Satellite and refine the system properties and ephemeris. We find a featureless transmission spectrum with increased absorption towards shorter wavelengths. We perform an atmospheric retrieval analysis with the AURA code, finding tentative evidence for haze in the upper atmospheric layers and a lower likelihood for a dense cloud deck. Whilst our retrieval analysis results point toward clouds and hazes, further evidence is needed to definitively reject a clear-sky scenario., Comment: 19 pages, 12 figures, 3 tables, Accepted for publication in MNRAS

We present a ground-based optical transmission spectrum for the warm Saturn-mass exoplanet WASP-110b from two transit observations made with the FOcal Reducer and Spectrograph (FORS2) on the Very Large Telescope (VLT). The spectrum covers the wavelength range from 4000 to 8333\AA, which is binned in 46 transit depths measured to an averaged precision of 220 parts per million (ppm) over an averaged 80\AA~bin for a Vmag=12.8 star. The measured transit depths are unaffected by a dilution from a close A-type field dwarf, which was fully resolved. The overall main characteristic of the transmission spectrum is an increasing radius with wavelength and a lack of the theoretically predicted pressure-broadened sodium and potassium absorption features for a cloud-free atmosphere. We analyze archival high-resolution optical spectroscopy and find evidence for low to moderate activity of the host star, which we take into account in the atmospheric retrieval analysis. Using the AURA retrieval code, we find that the observed transmission spectrum can be best explained by a combination of unocculted stellar faculae and a cloud deck. Transmission spectra of cloud-free and hazy atmospheres are rejected at a high confidence. With a possible cloud deck at its terminator, WASP-110b joins the increasing population of irradiated hot-Jupiter exoplanets with cloudy atmospheres observed in transmission., Comment: Accepted for publication in The Astronomical Journal

Thermal dissociation and recombination of molecular hydrogen, H_2, in the atmospheres of ultra-hot Jupiters (UHJs) has been shown to play an important role in global heat redistribution. This, in turn, significantly impacts their planetary emission, yet only limited investigations on the atmospheric effects have so far been conducted. Here we investigate the heat redistribution caused by this dissociation/recombination reaction, alongside feedback mechanisms between the atmospheric chemistry and radiative transfer, for a planetary and stellar configuration typical of UHJs. To do this, we have developed a time-dependent pseudo-2D model, including a treatment of time-independent equilibrium chemical effects. As a result of the reaction heat redistribution, we find temperature changes of up to $\sim$400 K in the atmosphere. When TiO and VO are additionally considered as opacity sources, these changes in temperature increase to over $\sim$800 K in some areas. This heat redistribution is found to significantly shift the region of peak atmospheric temperature, or hotspot, towards the evening terminator in both cases. The impact of varying the longitudinal wind speed on the reaction heat distribution is also investigated. When excluding TiO/VO, increased wind speeds are shown to increase the impact of the reaction heat redistribution up to a threshold wind speed. When including TiO/VO there is no apparent wind speed threshold, due to thermal stabilisation by these species. We also construct pseudo-2D phase curves from our model, and highlight both significant spectral flux damping and increased phase offset caused by the reaction heat redistribution., Comment: Accepted in MNRAS

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)

Hot Jupiters provide valuable natural laboratories for studying potential contributions of high-energy radiation to prebiotic synthesis in the atmospheres of exoplanets. In this fourth paper of the MOVES (Multiwavelength Observations of an eVaporating Exoplanet and its Star) programme, we study the effect of different types of high-energy radiation on the production of organic and prebiotic molecules in the atmosphere of the hot Jupiter HD 189733b. Our model combines X-ray and UV observations from the MOVES programme and 3D climate simulations from the 3D Met Office Unified Model to simulate the atmospheric composition and kinetic chemistry with the STAND2019 network. Also, the effects of galactic cosmic rays and stellar energetic particles are included. We find that the differences in the radiation field between the irradiated dayside and the shadowed nightside lead to stronger changes in the chemical abundances than the variability of the host star's XUV emission. We identify ammonium (NH4+) and oxonium (H3O+) as fingerprint ions for the ionization of the atmosphere by both galactic cosmic rays and stellar particles. All considered types of high-energy radiation have an enhancing effect on the abundance of key organic molecules such as hydrogen cyanide (HCN), formaldehyde (CH2O), and ethylene (C2H4). The latter two are intermediates in the production pathway of the amino acid glycine (C2H5NO2) and abundant enough to be potentially detectable by JWST., Comment: 15 pages, 11 figures, 4 tables

We present a publicly available library of model atmospheres with radiative-convective equilibrium Pressure-Temperature ($P$-$T$) profiles fully consistent with equilibrium chemical abundances, and the corresponding emission and transmission spectrum with R$\sim$5000 at 0.2 $\mu$m decreasing to R$\sim$35 at 30 $\mu$m, for 89 hot Jupiter exoplanets, for four re-circulation factors, six metallicities and six C/O ratios. We find the choice of condensation process (local/rainout) alters the $P$-$T$ profile and thereby the spectrum substantially, potentially detectable by JWST. We find H$^-$ opacity can contribute to form a strong temperature inversion in ultra-hot Jupiters for C/O ratios $\geq$ 1 and can make transmission spectra features flat in the optical, alongside altering the entire emission spectra. We highlight how adopting different model choices such as thermal ionisation, opacities, line-wing profiles and the methodology of varying the C/O ratio, effects the $P$-$T$ structure and the spectrum. We show the role of Fe opacity to form primary/secondary inversion in the atmosphere. We use WASP-17b and WASP-121b as test cases to demonstrate the effect of grid parameters across their full range, while highlighting some important findings, concerning the overall atmospheric structure, chemical transition regimes and their observables. Finally, we apply this library to the current transmission and emission spectra observations of WASP-121b, which shows H$_2$O and tentative evidence for VO at the limb, and H$_2$O emission feature indicative of inversion on the dayside, with very low energy redistribution, thereby demonstrating the applicability of library for planning and interpreting observations of transmission and emission spectrum., Comment: 26 pages, 19 figures in the main paper. 13 pages, 6 figures, 3 tables in the supplementary material attached with the main paper here. Accepted for Publication in MNRAS. Full grid of model P-T profiles, chemical abundances, transmission and emission spectra, contribution functions are available here, https://drive.google.com/drive/folders/1zCCe6HICuK2nLgnYJFal7W4lyunjU4JE

Identification of habitable planets beyond our solar system is a key goal of current and future space missions. Yet habitability depends not only on the stellar irradiance, but equally on constituent parts of the planetary atmosphere. Here we show, for the first time, that radiatively active mineral dust will have a significant impact on the habitability of Earth-like exoplanets. On tidally-locked planets, dust cools the day-side and warms the night-side, significantly widening the habitable zone. Independent of orbital configuration, we suggest that airborne dust can postpone planetary water loss at the inner edge of the habitable zone, through a feedback involving decreasing ocean coverage and increased dust loading. The inclusion of dust significantly obscures key biomarker gases (e.g. ozone, methane) in simulated transmission spectra, implying an important influence on the interpretation of observations. We demonstrate that future observational and theoretical studies of terrestrial exoplanets must consider the effect of dust., Comment: 15 pages, 4 figures, 3 tables

The majority of potentially habitable exoplanets detected orbit stars cooler than the Sun, and therefore are irradiated by a stellar spectrum peaking at longer wavelengths than that incident on Earth. Here, we present results from a set of simulations of tidally-locked terrestrial planets orbiting three different host stars to isolate the effect of the stellar spectra on the simulated climate. Specifically, we perform simulations based on TRAPPIST-1e, adopting an Earth-like atmosphere and using the UK Met Office Unified Model in an idealised 'aqua-planet' configuration. Whilst holding the planetary parameters constant, including the total stellar flux (900 W/m$^2$) and orbital period (6.10 Earth days), we compare results between simulations where the stellar spectrum is that of a quiescent TRAPPIST-1, Proxima Centauri and the Sun. The simulations with cooler host stars had an increased proportion of incident stellar radiation absorbed directly by the troposphere compared to the surface. This, in turn, led to an increase in the stability against convection, a reduction in overall cloud coverage on the dayside (reducing scattering), leading to warmer surface temperatures. The increased direct heating of the troposphere also led to more efficient heat transport from the dayside to the nightside and, therefore, a reduced day-night temperature contrast. We inferred that planets with an Earth-like atmosphere orbiting cooler stars had lower dayside cloud coverage, potentially allowing habitable conditions at increased orbital radii, compared to similar planets orbiting hotter stars for a given planetary rotation rate., Comment: 12 pages, 8 figures

We present four new secondary eclipse observations for the ultrahot Jupiter WASP-121b acquired using the Hubble Space Telescope Wide Field Camera 3. The eclipse depth is measured to a median precision of 60ppm across 28 spectroscopic channels spanning the 1.12-1.64 micron wavelength range. This is a considerable improvement to the 90ppm precision we achieved previously for a single eclipse observation using the same observing setup. Combining these data with those reported at other wavelengths, a blackbody spectrum for WASP-121b is ruled out at >6-sigma confidence and we confirm the interpretation of previous retrieval analyses that found the data is best explained by a dayside thermal inversion. The updated spectrum clearly resolves the water emission band at 1.3-1.6 micron, with higher signal-to-noise than before. It also fails to reproduce a bump in the spectrum at 1.25 micron derived from the first eclipse observation, which had tentatively been attributed to VO emission. We conclude the latter was either a statistical fluctuation or a systematic artefact specific to the first eclipse dataset., Comment: Accepted in MNRAS

Using a 3D general circulation model (GCM), we investigate the sensitivity of the climate of tidally-locked Earth-like exoplanets, Trappist-1e and Proxima Centauri b, to the choice of a convection parameterization. Compared to a mass-flux convection parameterization, a simplified convection adjustment parameterization leads to a $>$60% decrease of the cloud albedo, increasing the mean day-side temperature by $\approx$10 K. The representation of convection also affects the atmospheric conditions of the night side, via a change in planetary-scale wave patterns. As a result, using the convection adjustment scheme makes the night-side cold traps warmer by 17-36 K for the planets in our simulations. The day-night thermal contrast is sensitive to the representation of convection in 3D GCM simulations, so caution should be taken when interpreting emission phase curves. The choice of convection treatment, however, does not alter the simulated climate enough to result in a departure from habitable conditions, at least for the atmospheric composition and planetary parameters used in our study. The near-surface conditions both in the Trappist-1e and Proxima b cases remain temperate, allowing for an active water cycle. We further advance our analysis using high-resolution model experiments, in which atmospheric convection is simulated explicitly. Our results suggest that in a hypothetical global convection-permitting simulation the surface temperature contrast would be higher than in the coarse-resolution simulations with parameterized convection. In other words, models with parameterized convection may overestimate the inter-hemispheric heat redistribution efficiency., Comment: 26 pages, 15 figures, 5 tables. Accepted to ApJ

Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well., Comment: Accepted and published in Geoscientific Model Development (GMD)

We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the large-scale atmospheric flow in hot Jupiter atmospheres. We couple a flexible chemical kinetics scheme to the Met Office Unified Model which enables the study of the interaction of chemistry, radiative transfer and fluid dynamics. We use a newly-released "reduced" chemical network comprising 30 chemical species that has been specifically developed for application in 3D atmosphere models. We simulate the atmospheres of the well-studied hot Jupiters HD~209458b and HD~189733b which both have dayside--nightside temperature contrasts of several hundred Kelvin and superrotating equatorial jets. We find qualitatively quite different chemical structures between the two planets, particularly for methane (CH$_4$), when advection of chemical species is included. Our results show that consideration of 3D chemical transport is vital in understanding the chemical composition of hot Jupiter atmospheres. 3D mixing leads to significant changes in the abundances of absorbing gas-phase species compared with what would be expected by assuming local chemical equilibrium, or from models including 1D - and even 2D - chemical mixing. We find that CH$_4$, carbon dioxide (CO$_2$) and ammonia (NH$_3$) are particularly interesting as 3D mixing of these species leads to prominent signatures of out-of-equilibrium chemistry in the transmission and emission spectra, detectable with near-future instruments., Comment: Accepted in A&A

We use the Met Office Unified Model to explore the potential of a tidally locked M dwarf planet, nominally Proxima Centauri b irradiated by a quiescent version of its host star, to sustain an atmospheric ozone layer. We assume a slab ocean surface layer, and an Earth-like atmosphere of nitrogen and oxygen with trace amounts of ozone and water vapour. We describe ozone chemistry using the Chapman mechanism and the hydrogen oxide (HO$_x$, describing the sum of OH and HO$_2$) catalytic cycle. We find that Proxima Centauri radiates with sufficient UV energy to initialize the Chapman mechanism. The result is a thin but stable ozone layer that peaks at 0.75 parts per million at 25 km. The quasi-stationary distribution of atmospheric ozone is determined by photolysis driven by incoming stellar radiation and by atmospheric transport. Ozone mole fractions are smallest in the lowest 15 km of the atmosphere at the sub-stellar point and largest in the nightside gyres. Above 15 km the ozone distribution is dominated by an equatorial jet stream that circumnavigates the planet. The nightside ozone distribution is dominated by two cyclonic Rossby gyres that result in localized ozone hotspots. On the dayside the atmospheric lifetime is determined by the HO$_x$ catalytic cycle and deposition to the surface, with nightside lifetimes due to chemistry much longer than timescales associated with atmospheric transport. Surface UV values peak at the substellar point with values of 0.01 W/m$^2$, shielded by the overlying atmospheric ozone layer but more importantly by water vapour clouds., Comment: Accepted for publication in MNRAS. 25 pages, 17 figures

Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres., Comment: 28 pages, 20 pages of text - 8 of appendices, 9 figures in text - 6 in appendices

Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results. Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states. Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode. Results. Our code is benchmarked against classical wave and instability test cases in multiple geometries. The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes. Conclusions. We have developed and benchmarked a code for the study of linear processes in planetary atmospheres, with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably to further study the initial spin up of superrotation of GCM simulations of hot Jupiter., Comment: 18 pages (6 of appendix), 9 figures

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system., Comment: White paper for ESA Voyage 2050; 24 pages