Showing 1–46 of 46 results for author: Manners, J

3D simulations of TRAPPIST-1e with varying CO2, CH4 and haze profiles

Authors: Mei Ting Mak, Denis Sergeev, Nathan Mayne, Nahum Banks, Jake Eager-Nash, James Manners, Giada Arney, Eric Hebrard, Krisztian Kohary

Abstract: 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 temperat… ▽ More 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. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 17 pages, 16 figures, 2 tables

Quenching-driven equatorial depletion and limb asymmetries in hot Jupiter atmospheres: WASP-96b example

Authors: Maria Zamyatina, Duncan A. Christie, Eric Hébrard, Nathan J. Mayne, Michael Radica, Jake Taylor, Harry Baskett, Ben Moore, Craig Lils, Denis Sergeev, Eva-Maria Ahrer, James Manners, Krisztian Kohary, Adina D. Feinstein

Abstract: 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 quenc… ▽ More 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 μm$, could help distinguish atmospheres of different metallicities that are at chemical equilibrium from those with the upper layers at thermochemical disequilibrium. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: accepted to MNRAS, 27 pages, 24 figures

Modeling Atmospheric Lines By the Exoplanet Community (MALBEC) version 1.0: A CUISINES radiative transfer intercomparison project

Authors: Geronimo L. Villanueva, Thomas J. Fauchez, Vincent Kofman, Eleonora Alei, Elspeth K. H. Lee, Estelle Janin, Michael D. Himes, Jeremy Leconte, Michaela Leung, Sara Faggi, Mei Ting Mak, Denis E. Sergeev, Thea Kozakis, James Manners, Nathan Mayne, Edward W. Schwieterman, Alex R. Howe, Natasha Batalha

Abstract: 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 Line… ▽ More 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). △ Less

Submitted 6 February, 2024; originally announced February 2024.

3D simulations of the Archean Earth including photochemical haze profiles

Authors: M. T. Mak, N. J. Mayne, D. E. Sergeev, J. Manners, J. K. Eager-Nash, G. Arney, E. Hebrard, K. Kohary

Abstract: We present results from 3D simulations of the Archean Earth including a prescribed (non-interactive) spherical haze generated through a 1D photochemical model. Our simulations suggest that a thin haze layer, formed when CH4/CO2 = 0.1, leads to global warming of ~10.6 K due to the change of water vapour and cloud feedback, compared to the simulation without any haze. However, a thicker haze layer,… ▽ More We present results from 3D simulations of the Archean Earth including a prescribed (non-interactive) spherical haze generated through a 1D photochemical model. Our simulations suggest that a thin haze layer, formed when CH4/CO2 = 0.1, leads to global warming of ~10.6 K due to the change of water vapour and cloud feedback, compared to the simulation without any haze. However, a thicker haze layer, formed when CH4/CO2 > 0.1, leads to global cooling of up to ~65 K as the scattering and absorption of shortwave radiation from the haze reduces the radiation from reaching the planetary surface. A thermal inversion is formed with a lower tropopause as the CH4/CO2 ratio increases. The haze reaches an optical threshold thickness when CH4/CO2 ~ 0.175 beyond which the atmospheric structure and the global surface temperature do not vary much. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 30 pages, 22 figures, Accepted in Journal of Geophysical Research: Atmospheres

Simulations of idealised 3D atmospheric flows on terrestrial planets using LFRic-Atmosphere

Authors: Denis E. Sergeev, Nathan J. Mayne, Thomas Bendall, Ian A. Boutle, Alex Brown, Iva Kavcic, James Kent, Krisztian Kohary, James Manners, Thomas Melvin, Enrico Olivier, Lokesh K. Ragta, Ben J. Shipway, Jon Wakelin, Nigel Wood, Mohamed Zerroukat

Abstract: 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-… ▽ More 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. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 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)

3D climate simulations of the Archean find that methane has a strong cooling effect at high concentrations

Authors: Jake K. Eager-Nash, Nathan J. Mayne, Arwen E. Nicholson, Janke E. Prins, Oakley C. F. Young, Stuart J. Daines, Denis E. Sergeev, F. Hugo Lambert, James Manners, Ian A. Boutle, Eric T. Wolf, Inga E. E. Kamp, Krisztian Kohary, Tim M. Lenton

Abstract: 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$:… ▽ More 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. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 36 pages, 18 figures

A modern-day Mars climate in the Met Office Unified Model: dry simulations

Authors: Danny McCulloch, Denis E. Sergeev, Nathan Mayne, Matthew Bate, James Manners, Ian Boutle, Benjamin Drummond, Kristzian Kohary

Abstract: 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 accel… ▽ More 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étéorologie 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. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 20 pages, 12 figures with additional appendix figures

Report number: 16, 621--657

Journal ref: JOURNAL = Geoscientific Model Development, VOLUME = 16, YEAR = 2023, NUMBER = 2, PAGES = 621--657

Observability of signatures of transport-induced chemistry in clear atmospheres of hot gas giant exoplanets

Authors: Maria Zamyatina, Éric Hébrard, Benjamin Drummond, Nathan J. Mayne, James Manners, Duncan A. Christie, Pascal Tremblin, David K. Sing, Krisztian Kohary

Abstract: 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,… ▽ More 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. △ Less

Submitted 21 November, 2022; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: accepted to MNRAS, 27 pages, 15 figures, with DOI for the data repository

3D modelling of the impact of stellar activity on tidally locked terrestrial exoplanets: atmospheric composition and habitability

Authors: Robert J. Ridgway, Maria Zamyatina, Nathan J. Mayne, James Manners, F. Hugo Lambert, Marrick Braam, Benjamin Drummond, Éric Hébrard, Paul I. Palmer, Krisztian Kohary

Abstract: 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… ▽ More 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. △ Less

Submitted 25 October, 2022; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 27 pages, 21 figures, accepted to MNRAS

The Impact of Phase Equilibrium Cloud Models on GCM Simulations of GJ~1214b

Authors: D. A. Christie, N. J. Mayne, R. M. Gillard, J. Manners, E. Hébrard, S. Lines, K. Kohary

Abstract: We investigate the impact of clouds on the atmosphere of GJ~1214b using the radiatively-coupled, phase-equilibrium cloud model {\sc EddySed} coupled to the {\sc Unified Model} general circulation model. We find that, consistent with previous investigations, high metallicity ($100\times$ solar) and clouds with large vertical extents (a sedimentation factor of $f_\mathrm{sed} = 0.1$) are required to… ▽ More We investigate the impact of clouds on the atmosphere of GJ~1214b using the radiatively-coupled, phase-equilibrium cloud model {\sc EddySed} coupled to the {\sc Unified Model} general circulation model. We find that, consistent with previous investigations, high metallicity ($100\times$ solar) and clouds with large vertical extents (a sedimentation factor of $f_\mathrm{sed} = 0.1$) are required to best match the observations, although metallicities even higher than those investigated here may be required to improve agreement further. We additionally find that in our case which best matches the observations ($f_\mathrm{sed}=0.1$), the velocity structures change relative to the clear sky case with the formation of a superrotating jet being suppressed, although further investigation is required to understand the cause of the suppression. The increase in cloud extent with $f_\mathrm{sed}$ results in a cooler planet due to a higher albedo, causing the atmosphere to contract. This also results in a reduced day-night contrast seen in the phase curves, although the introduction of cloud still results in a reduction of the phase offset. We additionally investigate the impact the the {\sc Unified Model}'s pseudo-spherical irradiation scheme on the calculation of heating rates, finding that the introduction of nightside shortwave heating results in slower mid-latitude jets compared to the plane parallel irradiation scheme used in previous works. We also consider the impact of a gamma distribution, as opposed to a log-normal distribution, for the distribution of cloud particle radii and find the impact to be relatively minor. △ Less

Submitted 25 September, 2022; originally announced September 2022.

Comments: Accepted to MNRAS

Bistability of the atmospheric circulation on TRAPPIST-1e

Authors: Denis E. Sergeev, Neil T. Lewis, F. Hugo Lambert, Nathan J. Mayne, Ian A. Boutle, James Manners, Krisztian Kohary

Abstract: 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-… ▽ More 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. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 31 pages, 14 figures, accepted to the Planetary Science Journal

Longitudinally asymmetric stratospheric oscillation on a tidally locked exoplanet

Authors: Maureen Cohen, Massimo A. Bollasina, Paul I. Palmer, Denis E. Sergeev, Ian A. Boutle, Nathan J. Mayne, James Manners

Abstract: 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 k… ▽ More 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. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 18 pages, 7 figures

Journal ref: ApJ 930 152 (2022)

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI). Part III: Simulated Observables -- The return of the spectrum

Authors: Thomas J. Fauchez, Geronimo L. Villanueva, Denis E. Sergeev, Martin Turbet, Ian A. Boutle, Kostas Tsigaridis, Michael J. Way, Eric T. Wolf, Shawn D. Domagal-Goldman, Francois Forget, Jacob Haqq-Misra, Ravi K. Kopparapu, James Manners, Nathan J. Mayne

Abstract: 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 observ… ▽ More 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). △ Less

Submitted 15 September, 2022; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: Accepted in the Planetary Science Journal as Part III of a series of 3 THAI papers

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI). Part II: Moist Cases -- The Two Waterworlds

Authors: Denis E. Sergeev, Thomas J. Fauchez, Martin Turbet, Ian A. Boutle, Kostas Tsigaridis, Michael J. Way, Eric T. Wolf, Shawn D. Domagal-Goldman, Francois Forget, Jacob Haqq-Misra, Ravi K. Kopparapu, F. Hugo Lambert, James Manners, Nathan J. Mayne

Abstract: 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… ▽ More 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. △ Less

Submitted 15 September, 2022; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: 35 pages, 22 figures; Published in The Planetary Science Journal

Journal ref: Planet. Sci. J. (2022) 3 212

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI). Part I: Dry Cases -- The fellowship of the GCMs

Authors: Martin Turbet, Thomas J. Fauchez, Denis E. Sergeev, Ian A. Boutle, Kostas Tsigaridis, Michael J. Way, Eric T. Wolf, Shawn D. Domagal-Goldman, François Forget, Jacob Haqq-Misra, Ravi K. Kopparapu, F. Hugo Lambert, James Manners, Nathan J. Mayne, Linda Sohl

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 September, 2022; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: Published in the Planetary Science Journal as Part I of a series of 3 THAI papers

Journal ref: Martin Turbet et al 2022 Planet. Sci. J. 3 211

The Impact of Mixing Treatments on Cloud Modelling in 3D Simulations of Hot Jupiters

Authors: D. A. Christie, N. J. Mayne, S. Lines, V. Parmentier, J. Manners, I. Boutle, B. Drummond, T. Mikal-Evans, D. K. Sing, K. Kohary

Abstract: We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through th… ▽ More We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through the sedimentation factor $f_\mathrm{sed}$ plays a larger role in shaping cloud thickness and its radiative feedback on the local gas temperatures -- e.g. hot spot shift and day-to-night side temperature gradient -- than the switch in mixing treatment. We demonstrate using our new mixing treatments that simulations with cloud scales which are a fraction of the pressure scale height improve agreement with the observed transmission spectra, the emission spectra, and the Spitzer 4.5 $\mathrm{μm}$ phase curve, although our models are still unable to reproduce the optical and UV transmission spectra. We also find that the inclusion of cloud increases the transit asymmetry in the optical between the east and west limbs, although the difference remains small ($\lesssim 1\%$). △ Less

Submitted 12 July, 2021; originally announced July 2021.

Comments: 17 pages. 14 figures. 1 table

Mineral dust increases the habitability of terrestrial planets but confounds biomarker detection

Authors: Ian A. Boutle, Manoj Joshi, F. Hugo Lambert, Nathan J. Mayne, Duncan Lyster, James Manners, Robert Ridgway, Krisztian Kohary

Abstract: 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… ▽ More 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. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 15 pages, 4 figures, 3 tables

Journal ref: Nature Communications, 2020

Implications of different stellar spectra for the climate of tidally-locked Earth-like exoplanets

Authors: Jake K. Eager, David J. Reichelt, Nathan J. Mayne, F. Hugo Lambert, Denis E. Sergeev, Robert J. Ridgway, James Manners, Ian A. Boutle, Timothy M. Lenton, Krisztian Kohary

Abstract: 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. Spe… ▽ More 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. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 12 pages, 8 figures

Journal ref: A&A 639, A99 (2020)

Atmospheric convection plays a key role in the climate of tidally-locked terrestrial exoplanets: insights from high-resolution simulations

Authors: Denis E. Sergeev, F. Hugo Lambert, Nathan J. Mayne, Ian A. Boutle, James Manners, Krisztian Kohary

Abstract: 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 tem… ▽ More 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. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 26 pages, 15 figures, 5 tables. Accepted to ApJ

Implications of three-dimensional chemical transport in hot Jupiter atmospheres: results from a consistently coupled chemistry-radiation-hydrodynamics model

Authors: Benjamin Drummond, Eric Hebrard, Nathan J. Mayne, Olivia Venot, Robert J. Ridgway, Quentin Changeat, Shang-min Tsai, James Manners, Pascal Tremblin, Nathan Luke Abraham, David Sing, Krisztian Kohary

Abstract: 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 t… ▽ More 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. △ Less

Submitted 14 April, 2020; v1 submitted 30 January, 2020; originally announced January 2020.

Comments: Accepted in A&A

Ozone chemistry on tidally locked M dwarf planets

Authors: Jack S. Yates, Paul I. Palmer, James Manners, Ian Boutle, Krisztian Kohary, Nathan Mayne, Luke Abraham

Abstract: 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 m… ▽ More 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. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: Accepted for publication in MNRAS. 25 pages, 17 figures

Overcast on Osiris: 3D radiative-hydrodynamical simulations of a cloudy hot Jupiter using the parameterised, phase-equilibrium cloud formation code EddySed

Authors: S. Lines, N. J. Mayne, J. Manners, I. A. Boutle, B. Drummond, T. Mikal-Evans, K. Kohary, D. K. Sing

Abstract: We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effect… ▽ More We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effects, suggesting this important mechanism can not be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed $f_{\textrm{sed}}$ in EddySed), but also the temperature-pressure profile of the deeper atmosphere. We briefly discuss the large difference between the resolved cloud structures of this work, adopting a phase-equilibrium and parameterised cloud model, and our previous work incorporating a cloud microphysical model, although a fairer comparison where, for example, the same list of constituent condensates is included in both treatments, is reserved for a future work. Our results underline the importance of further study into the potential condensate size distributions and vertical structures, as both strongly influence the radiative impact of clouds on the atmosphere. Finally, we present synthetic observations from our simulations reporting an improved match, over our previous cloud-free simulations, to the observed transmission, HST WFC3 emission and 4.5 $μ$m Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations. △ Less

Submitted 27 June, 2019; originally announced June 2019.

Comments: 26 pages, accepted for publication in MNRAS

The Limits of the Primitive Equations of Dynamics for Warm, Slowly Rotating Small Neptunes and Super Earths

Authors: N. J. Mayne, B. Drummond, F. Debras, E. Jaupart, J. Manners, I. A. Boutle, I. Baraffe, K. Kohary

Abstract: We present significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practicall… ▽ More We present significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practically identical numerical setups, within the same model. The differences arise due to the breakdown of the so-called `shallow-fluid' and traditional approximations, which worsens when rotation rates are slowed, and day-night temperature contrasts are increased. The changes in the zonal advection between simulations solving the full and simplified equations, give rise to significant differences in the atmospheric redistribution of heat, altering the position of the hottest part of the atmosphere and temperature contrast between the day and night sides. The implications for the atmospheric chemistry and, therefore, observations need to be studied with a model including a more detailed treatment of the radiative transfer and chemistry. Small Neptunes and super Earths are extremely abundant and important, potentially bridging the structural properties (mass, radius, composition) of terrestrial and gas giant planets. Our results indicate care is required when interpreting the output of models solving the primitive equations of motion for such planets. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 26 pages, 11 figures, Accepted for Publication in ApJ

The 3D thermal, dynamical and chemical structure of the atmosphere of HD 189733b: implications of wind-driven chemistry for the emission phase curve

Authors: Benjamin Drummond, Nathan J. Mayne, James Manners, Isabelle Baraffe, Jayesh Goyal, Pascal Tremblin, David K. Sing, Krisztian Kohary

Abstract: In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind-… ▽ More In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind--driven advection on the chemical composition is qualitatively similar to our previous results for the warmer atmosphere of HD~209458b, found using the same model. However, we find more significant alterations to both the thermal and dynamical structure for the cooler atmosphere of HD~189733b, with changes in both the temperature and wind velocities reaching $\sim10\%$. We also present the contribution function, diagnosed from our simulations, and show that wind--driven chemistry has a significant impact on its three--dimensional structure, particularly for regions where methane is an important absorber. Finally, we present emission phase curves from our simulations and show the significant effect of wind--driven chemistry on the thermal emission, particularly within the 3.6 \textmu m Spitzer/IRAC channel. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: Accepted in ApJ

Exonephology: Transmission spectra from a 3D simulated cloudy atmosphere of HD209458b

Authors: S. Lines, J. Manners, N. J. Mayne, J. Goyal, A. L. Carter, I. A. Boutle, E. K. H. Lee, Ch. Helling, B. Drummond, D. M. Acreman, D. K. Sing

Abstract: We present high resolution transmission spectra, calculated directly from a 3D radiative-hydrodynamics simulation that includes kinetic cloud formation, for HD209458b. We find that the high opacity of our vertically extensive cloud deck, composed of a large number density of sub-micron particles, flattens the transmission spectrum and obscures spectral features identified in observed data. We use… ▽ More We present high resolution transmission spectra, calculated directly from a 3D radiative-hydrodynamics simulation that includes kinetic cloud formation, for HD209458b. We find that the high opacity of our vertically extensive cloud deck, composed of a large number density of sub-micron particles, flattens the transmission spectrum and obscures spectral features identified in observed data. We use the PandExo simulator to explore features of our HD209458b spectrum which may be detectable with the James Webb Space Telescope (JWST). We determine that an 8 - 12 micron absorption feature attributed to the mixed-composition, predominantly silicate cloud particles is a viable marker for the presence of cloud. Further calculations explore, and trends are identified with, variations in cloud opacity, composition heterogeneity and artificially scaled gravitational settling on the transmission spectrum. Principally, by varying the upper extent of our cloud decks, rainout is identified to be a key process for the dynamical atmospheres of hot-Jupiters and shown to dramatically alter the resulting spectrum. Our synthetic transmission spectra, obtained from the most complete, forward atmosphere simulations to--date, allow us to explore the model's ability to conform with observations. Such comparisons can provide insight into the physical processes either missing, or requiring improvement. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: Accepted to MNRAS, 14 Pages and 8 Figures

Simulating the cloudy atmospheres of HD 209458 b and HD 189733 b with the 3D Met Office Unified Model

Authors: S. Lines, N. J. Mayne, Ian A. Boutle, James Manners, Graham K. H. Lee, Ch. Helling, Benjamin Drummond, David S. Amundsen, Jayesh Goyal, David M. Acreman, Pascal Tremblin, Max Kerslake

Abstract: 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… ▽ More 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. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Comments: 28 pages, 23 figures and accepted for publication in A&A

Journal ref: A&A 615, A97 (2018)

Observable signatures of wind--driven chemistry with a fully consistent three dimensional radiative hydrodynamics model of HD 209458b

Authors: Benjamin Drummond, N. J. Mayne, James Manners, Aarynn L. Carter, Ian A. Boutle, Isabelle Baraffe, Eric Hebrard, Pascal Tremblin, David K. Sing, David S. Amundsen, Dave Acreman

Abstract: 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 r… ▽ More 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. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: Accepted for publication in ApJL

The influence of a sub-stellar continent on the climate of a tidally-locked exoplanet

Authors: Neil T. Lewis, F. Hugo Lambert, Ian A. Boutle, Nathan J. Mayne, James Manners, David M. Acreman

Abstract: Previous studies have demonstrated that continental carbon-silicate weathering is important to the continued habitability of a terrestrial planet. Despite this, few studies have considered the influence of land on the climate of a tidally-locked planet. In this work we use the Met Office Unified Model, coupled to a land surface model, to investigate the climate effects of a continent located at th… ▽ More Previous studies have demonstrated that continental carbon-silicate weathering is important to the continued habitability of a terrestrial planet. Despite this, few studies have considered the influence of land on the climate of a tidally-locked planet. In this work we use the Met Office Unified Model, coupled to a land surface model, to investigate the climate effects of a continent located at the sub-stellar point. We choose to use the orbital and planetary parameters of Proxima Centauri B as a template, to allow comparison with the work of others. A region of the surface where $T_{\text{s}} > 273.15\,\text{K}$ is always retained, and previous conclusions on the habitability of Proxima Centauri B remain intact. We find that sub-stellar land causes global cooling, and increases day-night temperature contrasts by limiting heat redistribution. Furthermore, we find that sub-stellar land is able to introduce a regime change in the atmospheric circulation. Specifically, when a continent offset to the east of the sub-stellar point is introduced, we observe the formation of two mid-latitude counterrotating jets, and a substantially weakened equatorial superrotating jet. △ Less

Submitted 1 February, 2018; originally announced February 2018.

Comments: 19 pages, 8 figures, Accepted for publication in The Astrophysical Journal

The effect of metallicity on the atmospheres of exoplanets with fully coupled 3D hydrodynamics, equilibrium chemistry, and radiative transfer

Authors: Benjamin Drummond, N. J. Mayne, Isabelle Baraffe, Pascal Tremblin, James Manners, David S. Amundsen, Jayesh Goyal, Dave Acreman

Abstract: 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 hea… ▽ More 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. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Comments: Accepted for publication in Astronomy and Astrophysics

A library of ATMO forward model transmission spectra for hot Jupiter exoplanets

Authors: Jayesh M. Goyal, Nathan Mayne, David K. Sing, Benjamin Drummond, Pascal Tremblin, David S. Amundsen, Thomas Evans, Aarynn L. Carter, Jessica Spake, Isabelle Baraffe, Nikolay Nikolov, James Manners, Gilles Chabrier, Eric Hebrard

Abstract: 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 o… ▽ More 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 $χ^{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. △ Less

Submitted 29 July, 2020; v1 submitted 27 October, 2017; originally announced October 2017.

Comments: No changes in replacement, main paper now made visible on arXiv as earlier the published Erratum associated with the paper was instead visible. 34 pages, 28 figures. Accepted for Publication in MNRAS. Full grid of model transmission spectra and chemical abundances are available here, https://drive.google.com/drive/folders/1tvlgSWyEAA0cX_yxnA5URD4Zo3Zj4GcH

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 474, Issue 4, March 2018, Pages 5158-5185

Advection of potential temperature in the atmosphere of irradiated exoplanets: a robust mechanism to explain radius inflation

Authors: P. Tremblin, G. Chabrier, N. J. Mayne, D. S. Amundsen, I. Baraffe, F. Debras, B. Drummond, J. Manners, S. Fromang

Abstract: The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosp… ▽ More The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosphere or oceans. At depth, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD209458b, but also the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non uniform atmospheric heating thus provides a robust mechanism explaining the inflated radii of irradiated hot Jupiters. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: accepted in ApJ

Results from a set of three-dimensional numerical experiments of a hot Jupiter atmosphere

Authors: N. J. Mayne, F. Debras, I. Baraffe, John Thuburn, David S. Amundsen, David M. Acreman, Chris Smith, Matthew K. Browning, James Manners, Nigel Wood

Abstract: We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom bou… ▽ More We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow. △ Less

Submitted 3 April, 2017; originally announced April 2017.

Comments: 26 pages, 22 Figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 604, A79 (2017)

Exploring the climate of Proxima B with the Met Office Unified Model

Authors: Ian A. Boutle, Nathan J. Mayne, Benjamin Drummond, James Manners, Jayesh Goyal, F. Hugo Lambert, David M. Acreman, Paul D. Earnshaw

Abstract: 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 plane… ▽ More 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). △ Less

Submitted 24 August, 2018; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: Astronomy and Astrophysics, in press

Journal ref: A&A 601, A120 (2017)

Treatment of overlapping gaseous absorption with the correlated-k method in hot Jupiter and brown dwarf atmosphere models

Authors: David S. Amundsen, Pascal Tremblin, James Manners, Isabelle Baraffe, Nathan J. Mayne

Abstract: The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models,… ▽ More The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: Accepted for publication in Astronomy & Astrophysics. 10 pages, 13 figures, 4 tables

Journal ref: A&A 598, A97 (2017)

The UK Met Office GCM with a sophisticated radiation scheme applied to the hot Jupiter HD 209458b

Authors: David S. Amundsen, Nathan J. Mayne, Isabelle Baraffe, James Manners, Pascal Tremblin, Benjamin Drummond, Chris Smith, David M. Acreman, Derek Homeier

Abstract: 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 c… ▽ More 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. △ Less

Submitted 9 September, 2016; v1 submitted 30 August, 2016; originally announced August 2016.

Comments: Accepted for publication in Astronomy & Astrophysics. 11 pages, 11 figures, 1 table

Journal ref: A&A 595, A36 (2016)

The mineral clouds on HD 209458b and HD189733b

Authors: Ch. Helling, G. Lee, I. Dobbs-Dixon, N. Mayne, D. S. Amundsen, J. Khaimova, A. A. Unger, J. Manners, D. Acreman, C. Smith

Abstract: 3D atmosphere model results are used to comparatively study the kinetic, non-equilibrium cloud formation in the atmospheres of two example planets guided by the giant gas planets HD209458b and HD189733b. Rather independently of hydrodynamic model differences, our cloud modelling suggests that both planets are covered in mineral clouds throughout the entire modelling domain. Both planets harbour ch… ▽ More 3D atmosphere model results are used to comparatively study the kinetic, non-equilibrium cloud formation in the atmospheres of two example planets guided by the giant gas planets HD209458b and HD189733b. Rather independently of hydrodynamic model differences, our cloud modelling suggests that both planets are covered in mineral clouds throughout the entire modelling domain. Both planets harbour chemically complex clouds that are made of mineral particles that have a height-dependent material composition and size. The remaining gas-phase element abundances strongly effects the molecular abundances of the atmosphere in the cloud forming regions. Hydrocarbon and cyanopolyyne molecules can be rather abundant in the inner, dense part of the atmospheres of HD189733b and HD209458b. No one value for metallicity and the C/O ratio can be used to describe an extrasolar planet. Our results concerning the presence and location of water in relation to the clouds explain some of the observed differences between the two planets. In HD189733b, strong water features have been reported while such features appear less strong for HD209458b. By considering the location of the clouds in the two atmospheres, we see that obscuring clouds exist high in the atmosphere of HD209458b, but much deeper in HD189733b. We further conclude that the (self-imposed) degeneracy of cloud parameters in retrieval methods can only be lifted if the cloud formation processes are accurately modelled in contrast to prescribing them by independent parameters. △ Less

Submitted 23 March, 2016; v1 submitted 13 March, 2016; originally announced March 2016.

Comments: 32 pages, accepted for publication in MNRAS (references corrected)

Accuracy tests of radiation schemes used in hot Jupiter global circulation models

Authors: David Skålid Amundsen, Isabelle Baraffe, Pascal Tremblin, James Manners, Wolfgang Hayek, N. J. Mayne, David M. Acreman

Abstract: The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have be… ▽ More The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution. △ Less

Submitted 18 February, 2014; v1 submitted 4 February, 2014; originally announced February 2014.

Comments: Accepted for publication in Astronomy & Astrophysics. 16 pages, 17 figures, 11 tables

Complete Multiwavelength Characterization of Faint Chandra X-ray Sources Seen in the Spitzer Wide-Area IR Extragalactic (SWIRE) Survey

Authors: A. Franceschini, J. Manners, M. Polletta, C. J. Lonsdale, E. Gonzalez-Solares, J. Surace, D. Shupe, F. Fan, C. K. Xu, D. Farrah, S. Berta, G. Rodighiero, I. Perez-Fournon, E. Hatziminaoglou, H. E. Smith, B. Siana, M. Rowan-Robinson, K. Nandra, T. Babbedge, M. Vaccari, S. Oliver, B. Wilkes, F. Owe, D. Padgett, D. Frayer , et al. (2 additional authors not shown)

Abstract: We exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at l… ▽ More We exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at least 2 IRAC bands and 988 sources detected at 24micron with MIPS brighter than 0.1 mJy. Of these, 102 IRAC and 59 MIPS sources have Chandra counterparts, out of a total of 122 X-ray sources present in the area with S(0.5-8 kev)>10^(-15) erg/cm^2/s. We have constructed SEDs for each source using data from the 4 IRAC wavebands, Chandra fluxes, and optical follow-up data in the wavebands U, g', r', i', Z, and H. We fit a number of spectral templates to the SEDs at optical and infrared wavelengths to determine photometric redshifts and spectral categories, and also make use of diagnostics based on the X-ray luminosities, hardness ratios, X-ray to infrared spectral slopes and optical morphologies. Although we have spectroscopic redshifts for only a minority of the Chandra sources, the available SEDs constrain the redshifts for most of the sample sources, which turn out to be typically at 0.5<z<2. We find that 39% of the Chandra sources are dominated by type-1 AGN emission, 23% display optical/IR spectra typical of type-2 AGNs, while the remaining 38% fraction show starburst-like or even normal galaxy spectra. Since we prove that all these galaxies are dominated by AGN emission in X-rays this brings the fraction of type-1 AGNs to be 80% of the type-2: even assuming that all the Chandra sources undetected by Spitzer are type-2 AGNs, the type-1 fraction would exceed 1/3 of the total population (abridged). △ Less

Submitted 17 December, 2004; originally announced December 2004.

Comments: Accepted for publication in AJ, March 2005 issue

Mid-infrared sources in the ELAIS Deep X-ray Survey

Authors: J. C. Manners, S. Serjeant, S. Bottinelli, M. Vaccari, A. Franceschini, I. Perez-Fournon, E. Gonzalez-Solares, C. J. Willott, O. Johnson, O. Almaini, M. Rowan-Robinson, S. Oliver

Abstract: We present a cross-correlation of the European Large Area ISO survey (ELAIS) with the ELAIS Deep X-ray Survey of the N1 and N2 fields. There are 7 Chandra point sources with matches in the ELAIS Final Analysis 15um catalogue, out of a total of 28 extragalactic ISO sources present in the Chandra fields. Five of these are consistent with AGN giving an AGN fraction of ~19 per cent in the 15um flux… ▽ More We present a cross-correlation of the European Large Area ISO survey (ELAIS) with the ELAIS Deep X-ray Survey of the N1 and N2 fields. There are 7 Chandra point sources with matches in the ELAIS Final Analysis 15um catalogue, out of a total of 28 extragalactic ISO sources present in the Chandra fields. Five of these are consistent with AGN giving an AGN fraction of ~19 per cent in the 15um flux range 0.8-6 mJy. We have co-added the hard X-ray fluxes of the individually-undetected ISO sources and find a low significance detection consistent with star formation in the remaining population. We combine our point source cross-correlation fraction with the XMM-Newton observations of the Lockman Hole and Chandra observations of the Hubble Deep Field North to constrain source count models of the mid-infrared galaxy population. The low dust-enshrouded AGN fraction in ELAIS implied by the number of cross-identifications between the ELAIS mid-infrared sample and the Chandra point sources is encouraging for the use of mid-infrared surveys to constrain the cosmic star formation history, provided there are not further large undetected populations of Compton-thick AGN. △ Less

Submitted 29 September, 2004; originally announced September 2004.

Comments: 9 pages, 6 figures, MNRAS accepted

Dust and gas obscuration in ELAIS Deep X-ray Survey reddened quasars

Authors: Chris J. Willott, Chris Simpson, Omar Almaini, Olivia Johnson, Andrew Lawrence, James S. Dunlop, Nathan D. Roche, Robert G. Mann, James C. Manners, Eduardo Gonzalez-Solares, Ismael Perez-Fournon, Rob J. Ivison, Stephen Serjeant, Seb J. Oliver, Richard G. McMahon, Michael Rowan-Robinson

Abstract: Hard X-ray surveys have uncovered a large population of heavily obscured AGN. They also reveal a population of quasars with moderate obscuration at both visible and X-ray wavelengths. We use Chandra selected samples of quasars from the ELAIS Deep X-ray Survey (EDXS) and CDF-N to investigate the obscuration towards the nuclei of moderately obscured AGN. We find an inverse correlation between the… ▽ More Hard X-ray surveys have uncovered a large population of heavily obscured AGN. They also reveal a population of quasars with moderate obscuration at both visible and X-ray wavelengths. We use Chandra selected samples of quasars from the ELAIS Deep X-ray Survey (EDXS) and CDF-N to investigate the obscuration towards the nuclei of moderately obscured AGN. We find an inverse correlation between the optical to X-ray flux ratio and the X-ray hardness ratio which can be interpreted as due to obscuration at visible and X-ray wavelengths. We present detailed optical and near-infrared data for a sample of optically-faint (R>23) quasars from the EDXS. These are used to constrain the amount of rest-frame UV/optical reddening towards these quasars. It is found that optically-faint quasars are mostly faint due to obscuration, not because they are intrinsically weak. After correcting for reddening, the optical magnitudes of most of these quasars are similar to the brighter quasars at these X-ray fluxes. Combining with gas column densities inferred from the X-ray observations we consider the gas-to-dust ratios of the obscuring matter. We find that the quasars generally have higher gas-to-dust absorption than that seen in the Milky Way - similar to what has been found for nearby Seyfert galaxies. We consider the possible existence of a large population of X-ray sources which have optical properties of Type 1 (unobscured) quasars, but X-ray properties of Type 2 (obscured) quasars. We show that the observed distribution of optical-to-X-ray flux ratios of quasars at z>1 is skewed to low values compared to the intrinsic distribution... △ Less

Submitted 18 March, 2004; originally announced March 2004.

Comments: 12 pages, 10 figures, ApJ, in press

Journal ref: Astrophys.J. 610 (2004) 140-150

Obscured AGN from the ELAIS Deep X-ray Survey

Authors: C. J. Willott, C. Simpson, O. Almaini, J. C. Manners, O. Johnson, A. Lawrence, J. S. Dunlop, R. J. Ivison, S. Rawlings, E. Gonzalez-Solares, I. Perez-Fournon, S. Serjeant, S. J. Oliver, N. D. Roche, R. G. Mann, M. Rowan-Robinson

Abstract: The sources discovered in deep hard X-ray surveys with 2-8 keV fluxes of 10^-14 erg cm^-2 s^-1 make up the bulk of the X-ray background at these energies. We present here detailed multi-wavelength observations of three such sources from the ELAIS Deep X-ray Survey. The observations include sensitive near-infrared spectroscopy with the Subaru Telescope and X-ray spectral information from the Chan… ▽ More The sources discovered in deep hard X-ray surveys with 2-8 keV fluxes of 10^-14 erg cm^-2 s^-1 make up the bulk of the X-ray background at these energies. We present here detailed multi-wavelength observations of three such sources from the ELAIS Deep X-ray Survey. The observations include sensitive near-infrared spectroscopy with the Subaru Telescope and X-ray spectral information from the Chandra X-ray Observatory. The sources observed all have optical-to-near-IR colours redder than an unobscured quasar and comprise a reddened quasar, a radio galaxy and an optically-obscured AGN. The reddened quasar is at a redshift z=2.61 and shows a very large X-ray absorbing column of N_H approx 3.10^23 cm^-2. This contrasts with the relatively small amount of dust reddening, implying a gas-to-dust ratio along the line-of-sight a hundred times greater than that of the Milky Way. The radio galaxy at z=1.57 shows only narrow emission lines, but has a surprisingly soft X-ray spectrum. The softness of this spectrum either indicates an unusually low gas-to-dust ratio for the absorbing medium or X-ray emission related to the young radio source. The host galaxy is extremely red (R-K=6.4) and its optical/near-IR spectrum is best fit by a strongly reddened (A_V~2) starburst. The third X-ray source discussed is also extremely red (R-K=6.1) and lies in a close grouping of three other R-K>6 galaxies. No emission or absorption lines were detected from this object, but its redshift (and that of one of the nearby galaxies) are constrained by SED-fitting to be just greater than z=1. The extremely red colours of these two galaxies can be accounted for by old stellar populations. These observations illustrate the diverse properties of hard X-ray selected AGN. △ Less

Submitted 10 October, 2002; originally announced October 2002.

Comments: 14 pages, 11 figures, MNRAS in press

Journal ref: Mon.Not.Roy.Astron.Soc. 339 (2003) 397

The ELAIS Deep X-ray Survey I: Chandra Source Catalogue and First Results

Authors: J. C. Manners, O. Johnson, O. Almaini, C. J. Willott, E. Gonzalez-Solares, A. Lawrence, R. G. Mann, I. Perez-Fournon, J. S. Dunlop, R. G. McMahon, S. J. Oliver, M. Rowan-Robinson, S. Serjeant

Abstract: We present an analysis of two deep (75 ks) Chandra observations of the European Large Area ISO Survey (ELAIS) fields N1 and N2 as the first results from the ELAIS deep X-ray survey. This survey is being conducted in well studied regions with extensive multi-wavelength coverage. Here we present the Chandra source catalogues along with an analysis of source counts, hardness ratios and optical clas… ▽ More We present an analysis of two deep (75 ks) Chandra observations of the European Large Area ISO Survey (ELAIS) fields N1 and N2 as the first results from the ELAIS deep X-ray survey. This survey is being conducted in well studied regions with extensive multi-wavelength coverage. Here we present the Chandra source catalogues along with an analysis of source counts, hardness ratios and optical classifications. A total of 233 X-ray point sources are detected in addition to 2 soft extended sources, which are found to be associated with galaxy clusters. An over-density of sources is found in N1 with 30% more sources than N2, which we attribute to large-scale structure. A similar variance is seen between other deep Chandra surveys. The source count statistics reveal an increasing fraction of hard sources at fainter fluxes. The number of galaxy-like counterparts also increases dramatically towards fainter fluxes, consistent with the emergence of a large population of obscured sources. △ Less

Submitted 1 April, 2003; v1 submitted 29 July, 2002; originally announced July 2002.

Comments: 17 pages, 15 figures, MNRAS accepted; fig 2 corrected, differential counts added, text revised and updated. Further information including electronic source catalogues can be found at http://www.roe.ac.uk/~jcm/edxs/index.html

Journal ref: Mon.Not.Roy.Astron.Soc. 343 (2003) 293

The X-ray variability of high redshift QSOs

Authors: J. Manners, O. Almaini, A. Lawrence

Abstract: We present an analysis of X-ray variability in a sample of 156 radio quiet quasars taken from the ROSAT archive, covering a redshift range 0.1 < z < 4.1. A maximum likelihood method is used to constrain the amplitude of variability in these low signal to noise light curves. Through combining these data in ensembles we are able to identify trends in variability amplitude with luminosity and with… ▽ More We present an analysis of X-ray variability in a sample of 156 radio quiet quasars taken from the ROSAT archive, covering a redshift range 0.1 < z < 4.1. A maximum likelihood method is used to constrain the amplitude of variability in these low signal to noise light curves. Through combining these data in ensembles we are able to identify trends in variability amplitude with luminosity and with redshift. The decline in variability amplitude with luminosity identified in local AGN (z < 0.1) is confirmed out to z=2. There is tentative evidence for an increase in QSO X-ray variability amplitude towards high redshifts (z > 2) in the sense that QSOs of the same X-ray luminosity are more variable at z > 2. We discuss possible explanations for this effect. The simplest explanation may be that high redshift QSOs are accreting at a larger fraction of the Eddington limit than local AGN. △ Less

Submitted 22 October, 2001; originally announced October 2001.

Comments: 9 pages, 6 figures, MNRAS accepted

Journal ref: Mon.Not.Roy.Astron.Soc. 330 (2002) 390

The coincidence and angular clustering of Chandra and SCUBA sources

Authors: O. Almaini, S. E. Scott, J. S. Dunlop, J. C. Manners, C. J. Willott, A. Lawrence, R. J. Ivison, O. Johnson, A. W. Blain, J. A. Peacock, S. J. Oliver, M. J. Fox, R. G. Mann, I. Perez-Fournon, E. Gonzalez-Solares, M. Rowan-Robinson, S. Serjeant, F. Cabrera-Guerra, D. H. Hughes

Abstract: We explore the relationship between the hard X-ray and sub-mm populations using deep Chandra observations of a large, contiguous SCUBA survey. In agreement with other recent findings, we confirm that the direct overlap is small. Of the 17 sub-mm sources detected in this field at 850 microns, only one is coincident with a Chandra source. The resulting limits imply that the majority of SCUBA sourc… ▽ More We explore the relationship between the hard X-ray and sub-mm populations using deep Chandra observations of a large, contiguous SCUBA survey. In agreement with other recent findings, we confirm that the direct overlap is small. Of the 17 sub-mm sources detected in this field at 850 microns, only one is coincident with a Chandra source. The resulting limits imply that the majority of SCUBA sources are not powered by AGN, unless the central engine is obscured by Compton-thick material with a low (<1 per cent) scattered component. Furthermore, since Chandra detects only ~5 per cent of SCUBA sources, the typical obscuration would need to be almost isotropic. The X-ray upper limits are so strong that in most cases we can also rule out a starburst SED at low redshift, suggesting that the majority of SCUBA sources lie at z>1 even if they are purely starburst galaxies. Despite the low detection rate, we find evidence for strong angular clustering between the X-ray and sub-mm populations. The implication is that AGN and SCUBA sources trace the same large-scale structure but do not generally coincide. If bright sub-mm sources represent massive elliptical galaxies in formation, we suggest that (for a given galaxy) the major episode of star-formation must be distinct from the period of observable quasar activity. △ Less

Submitted 28 August, 2002; v1 submitted 24 August, 2001; originally announced August 2001.

Comments: 12 pages, 8 figures, MNRAS in press. Minor revision, including a comparison with the SED of Arp-220 and an X-ray stacking analysis for the SCUBA sources

Journal ref: Mon.Not.Roy.Astron.Soc. 338 (2003) 303-311

The ELAIS Deep X-ray Survey

Authors: C. J. Willott, O. Almaini, J. Manners, O. Johnson, A. Lawrence, J. S. Dunlop, R. G Mann, I. Perez-Fournon, E. Gonzalez-Solares, F. Cabrera-Guerra, S. Serjeant, S. J. Oliver, M. Rowan-Robinson

Abstract: We present initial follow-up results of the ELAIS Deep X-ray Survey which is being undertaken with the Chandra and XMM-Newton Observatories. 235 X-ray sources are detected in our two 75 ks ACIS-I observations in the well-studied ELAIS N1 and N2 areas. 90% of the X-ray sources are identified optically to R=26 with a median magnitude of R=24. We show that objects which are unresolved optically (i.… ▽ More We present initial follow-up results of the ELAIS Deep X-ray Survey which is being undertaken with the Chandra and XMM-Newton Observatories. 235 X-ray sources are detected in our two 75 ks ACIS-I observations in the well-studied ELAIS N1 and N2 areas. 90% of the X-ray sources are identified optically to R=26 with a median magnitude of R=24. We show that objects which are unresolved optically (i.e. quasars) follow a correlation between their optical and X-ray fluxes, whereas galaxies do not. We also find that the quasars with fainter optical counterparts have harder X-ray spectra, consistent with absorption at both wavebands. Initial spectroscopic follow-up has revealed a large fraction of high-luminosity Type 2 quasars. The prospects for studying the evolution of the host galaxies of X-ray selected Type 2 AGN are considered. △ Less

Submitted 31 May, 2001; originally announced May 2001.

Comments: 9 pages, 5 figures, To appear in Proceedings of XXI Moriond Conference: "Galaxy Clusters and the High Redshift Universe Observed in X-rays", edited by D. Neumann, F.Durret, & J. Tran Thanh Van

Starburst activity in a ROSAT Narrow Emission-Line Galaxy

Authors: K. F. Gunn, I. M. McHardy, O. Almaini, T. Shanks, T. J. Sumner, T. W. B. Muxlow, A. Efstathiou, L. R. Jones, S. M. Croom, J. C. Manners, A. M. Newsam, K. O. Mason, S. B. G. Serjeant, M. Rowan-Robinson, .

Abstract: We present multiwaveband photometric and optical spectropolarimetric observations of the R=15.9 narrow emission line galaxy R117_A which lies on the edge of the error circle of the ROSAT X-ray source R117 (from McHardy et al 1998). The overall spectral energy distribution of the galaxy is well modelled by a combination of a normal spiral galaxy and a moderate-strength burst of star formation. Th… ▽ More We present multiwaveband photometric and optical spectropolarimetric observations of the R=15.9 narrow emission line galaxy R117_A which lies on the edge of the error circle of the ROSAT X-ray source R117 (from McHardy et al 1998). The overall spectral energy distribution of the galaxy is well modelled by a combination of a normal spiral galaxy and a moderate-strength burst of star formation. The far infra-red and radio emission is extended along the major axis of the galaxy, indicating an extended starburst. On positional grounds, the galaxy is a good candidate for the identification of R117 and the observed X-ray flux is very close to what would be expected from a starburst of the observed far infra-red and radio fluxes. Although an obscured high redshift QSO cannot be entirely ruled out as contributing some fraction of the X-ray flux, we find no candidates to K=20.8 within the X-ray errorbox and so conclude that R117_A is responsible for a large fraction, if not all, of the X-ray emission from R117. Searches for indicators of an obscured AGN in R117_A have so far proven negative; deep spectropolarimetric observations show no signs of broad lines to a limit of one per cent and, for the observed far infra-red and radio emission, we would expect a ten times greater X-ray flux if the overall emission were powered by an AGN. We therefore conclude that the X-ray emission from R117 is dominated by starburst emission from the galaxy R117_A. △ Less

Submitted 24 October, 2000; originally announced October 2000.

Comments: 8 pages, 9 figures, MNRAS accepted

Journal ref: Mon.Not.Roy.Astron.Soc. 324 (2001) 305