• We present results of simulations of the climate of the newly discovered planet Proxima Centauri B, performed using the Met Office Unified Model (UM). We examine the responses of both an `Earth-like' atmosphere and simplified nitrogen and trace carbon dioxide atmosphere to the radiation likely received by Proxima Centauri B. Additionally, we explore the effects of orbital eccentricity on the planetary conditions using a range of eccentricities guided by the observational constraints. Overall, our results are in agreement with previous studies in suggesting Proxima Centauri B may well have surface temperatures conducive to the presence of liquid water. Moreover, we have expanded the parameter regime over which the planet may support liquid water to higher values of eccentricity (>= 0.1) and lower incident fluxes (881.7 Wm-2) than previous work. This increased parameter space arises because of the low sensitivity of the planet to changes in stellar flux, a consequence of the stellar spectrum and orbital configuration. However, we also find interesting differences from previous simulations, such as cooler mean surface temperatures for the tidally-locked case. Finally, we have produced high resolution planetary emission and reflectance spectra, and highlight signatures of gases vital to the evolution of complex life on Earth (oxygen, ozone and carbon dioxide).
  • Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. Current observational instruments (on board the Hubble Space Telescope, Spitzer, and on ground-based facilities) have allowed the scientific community to obtain important information on the physical and chemical properties of these planets. However, for a more in-depth characterisation of these worlds, more powerful telescopes are needed. Thanks to the high sensitivity of their instruments, the next generation of space observatories (e.g. JWST, ARIEL) will provide observations of unprecedented quality, allowing us to extract far more information than what was previously possible. Such high quality observations will provide constraints on theoretical models of exoplanet atmospheres and lead to a greater understanding of their physics and chemistry. Important modelling efforts have been carried out during the past few years, showing that numerous parameters and processes (such as the elemental abundances, temperature, mixing, etc.) are likely to affect the atmospheric composition of exoplanets and subsequently the observable spectra. In this manuscript, we review the different parameters that can influence the molecular composition of exoplanet atmospheres. We show that the high-precision of ARIEL observations will improve our view and characterisation of exoplanet atmospheres. We also consider future developments that are necessary to improve atmospheric models, driven by the need to interpret the available observations.
  • Jacob L. Bean, Kevin B. Stevenson, Natalie M. Batalha, Zachory Berta-Thompson, Laura Kreidberg, Nicolas Crouzet, Björn Benneke, Michael R. Line, David K. Sing, Hannah R. Wakeford, Heather A. Knutson, Eliza M.-R. Kempton, Jean-Michel Désert, Ian Crossfield, Natasha E. Batalha, Julien de Wit, Vivien Parmentier, Joseph Harrington, Julianne I. Moses, Mercedes Lopez-Morales, Munazza K. Alam, Jasmina Blecic, Giovanni Bruno, Aarynn L. Carter, John W. Chapman, Leen Decin, Diana Dragomir, Thomas M. Evans, Jonathan J. Fortney, Jonathan D. Fraine, Peter Gao, Antonio García Muñoz, Neale P. Gibson, Jayesh M. Goyal, Kevin Heng, Renyu Hu, Sarah Kendrew, Brian M. Kilpatrick, Jessica Krick, Pierre-Olivier Lagage, Monika Lendl, Tom Louden, Nikku Madhusudhan, Avi M. Mandell, Megan Mansfield, Erin M. May, Giuseppe Morello, Caroline V. Morley, Nikolay Nikolov, Seth Redfield, Jessica E. Roberts, Everett Schlawin, Jessica J. Spake, Kamen O. Todorov, Angelos Tsiaras, Olivia Venot, William C. Waalkes, Peter J. Wheatley, Robert T. Zellem, Daniel Angerhausen, David Barrado, Ludmila Carone, Sarah L. Casewell, Patricio E. Cubillos, Mario Damiano, Miguel de Val-Borro, Benjamin Drummond, Billy Edwards, Michael Endl, Nestor Espinoza, Kevin France, John E. Gizis, Thomas P. Greene, Thomas K. Henning, Yucian Hong, James G. Ingalls, Nicolas Iro, Patrick G.J. Irwin, Tiffany Kataria, Fred Lahuis, Jérémy Leconte, Jorge Lillo-Box, Stefan Lines, Luigi Mancini, Franck Marchis, Nathan Mayne, Enric Palle, Gaël Roudier, Evgenya L. Shkolnik, John Southworth, Johanna Teske, Giovanna Tinetti, Pascal Tremblin, Gregory S. Tucker, Roy vanBoekel, Ingo P. Waldmann, Ian C. Weaver, Tiziano Zingales
    March 13, 2018 astro-ph.EP
    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.
  • To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics and thermal profile. We couple the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self--consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption and, crucially, scattering. Fluxes directly obtained from the UM are used to produce synthetic SEDs and phase curves. Our simulations show extensive cloud formation in both planets, however, cooler temperatures in the HD 189733 b result in higher cloud particle number densities. Sub-micron particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, result in enhanced cloud coverage for latitudes > 30 degrees and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet(s), as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, our simulated optical phase curve for HD 209458 b shows no `offset', in contrast to observations. Efficient scattering by cloud results an atmospheric cooling of up to 250K, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements and suggests both a note of caution of interpretations neglecting such cloud feedback and scattering, and merits further study.
  • We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD~209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model-observation discrepancy in the 4.5\,{\textmu m} {\it Spitzer}/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominant.
  • In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.
  • WASP-39b is a hot Saturn-mass exoplanet with a predicted clear atmosphere based on observations in the optical and infrared. Here we complete the transmission spectrum of the atmosphere with observations in the near-infrared (NIR) over three water absorption features with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) G102 (0.8-1.1 microns) and G141 (1.1-1.7 microns) spectroscopic grisms. We measure the predicted high amplitude H2O feature centered at 1.4 microns, and the smaller amplitude features at 0.95 and 1.2 microns, with a maximum water absorption amplitude of 2.4 planetary scale heights. We incorporate these new NIR measurements into previously published observational measurements to complete the transmission spectrum from 0.3-5 microns. From these observed water features, combined with features in the optical and IR, we retrieve a well constrained temperature Teq = 1030(+30,-20) K, and atmospheric metallicity 151 (+48,-46)x solar which is relatively high with respect to the currently established mass-metallicity trends. This new measurement in the Saturn-mass range hints at further diversity in the planet formation process relative to our solar system giants.
  • We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally significant hot exoplanets (Equilibrium Temperatures of 547-2710 K). This model grid has been developed using a 1D radiative-convective-chemical equilibrium model termed ATMO, with up-to-date high temperature opacities. We present an interpretation of observations of ten exoplanets, including best fit parameters and $\chi^{2}$ maps. In agreement with previous works, we find a continuum from clear to hazy/cloudy atmospheres for this sample of hot Jupiters. The data for all the 10 planets are consistent with sub-solar to solar C/O ratio, 0.005 to 10 times solar metallicity and water rather than a methane dominated infrared spectra. We then explore the range of simulated atmospheric spectra for different exoplanets, based on characteristics such as temperature, metallicity, C/O-ratio, haziness and cloudiness. We find a transition value for the metallicity between 10 and 50 times solar, which leads to substantial changes in the transmission spectra. We also find a transition value of C/O ratio, from water to carbon species dominated infrared spectra, as found by previous works, revealing a temperature dependence of this transition point ranging from $\sim$0.56 to $\sim$1-1.3 for equilibrium temperatures from $\sim$900 to $\sim$2600 K. We highlight the potential of the spectral features of HCN and C$_2$H$_2$ to constrain the metallicities and C/O ratios of planets, using JWST observations. Finally, our entire grid ($\sim$460,000 simulations) is publicly available and can be used directly with the JWST simulator PandExo for planning observations.
  • To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office global circulation model (GCM), the Unified Model (UM), to these exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the coupling of the UM dynamical core to an accurate radiation scheme based on the two-stream approximation and correlated-k method with state-of-the-art opacities from ExoMol. Our first application of this model is devoted to the extensively studied hot Jupiter HD 209458b. We have derived synthetic emission spectra and phase curves, and compare them to both previous models also based on state-of-the-art radiative transfer, and to observations. We find a reasonable a agreement between observations and both our day side emission and hot spot offset, however, our night side emission is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009) with the SPARC/MITgcm, however, we note several quantitative differences: Our simulations show significant variation in the position of the hottest part of the atmosphere with pressure, as expected from simple timescale arguments, and in contrast to the "vertical coherency" found by Showman et al. (2009). We also see significant quantitative differences in calculated synthetic observations. Our comparisons strengthen the need for detailed intercomparisons of dynamical cores, radiation schemes and post-processing tools to understand these differences. This effort is necessary in order to make robust conclusions about these atmospheres based on GCM results.
  • In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure--temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.
  • We present a primary transit observation for the ultra hot (Teq~2400K) gas giant expolanet WASP-121b, made using the Hubble Space Telescope Wide Field Camera 3 in spectroscopic mode across the 1.12-1.64 micron wavelength range. The 1.4 micron water absorption band is detected at high confidence (5.4 sigma) in the planetary atmosphere. We also reanalyze ground-based photometric lightcurves taken in the B, r', and z' filters. Significantly deeper transits are measured in these optical bandpasses relative to the near-infrared wavelengths. We conclude that scattering by high-altitude haze alone is unlikely to account for this difference, and instead interpret it as evidence for titanium oxide and vanadium oxide absorption. Enhanced opacity is also inferred across the 1.12-1.3 micron wavelength range, possibly due to iron hydride absorption. If confirmed, WASP-121b will be the first exoplanet with titanium oxide, vanadium oxide, and iron hydride detected in transmission. The latter are important species in M/L dwarfs, and their presence is likely to have a significant effect on the overall physics and chemistry of the atmosphere, including the production of a strong thermal inversion.