-
Asymmetries in the simulated ozone distribution on TRAPPIST-1e due to orography
Authors:
Anand Bhongade,
Daniel R Marsh,
Felix Sainsbury-Martinez,
Gregory J Cooke
Abstract:
TRAPPIST-1e is a tidally locked rocky exoplanet orbiting the habitable zone of an M dwarf star. Upcoming observations are expected to reveal new rocky exoplanets and their atmospheres around M dwarf stars. To interpret these future observations we need to model the atmospheres of such exoplanets. We configured CESM2-WACCM6, a chemistry climate model, for the orbit and stellar irradiance of TRAPPIS…
▽ More
TRAPPIST-1e is a tidally locked rocky exoplanet orbiting the habitable zone of an M dwarf star. Upcoming observations are expected to reveal new rocky exoplanets and their atmospheres around M dwarf stars. To interpret these future observations we need to model the atmospheres of such exoplanets. We configured CESM2-WACCM6, a chemistry climate model, for the orbit and stellar irradiance of TRAPPIST-1e assuming an initial Earth-like atmospheric composition. Our aim is to characterize the possible ozone (O$_3$) distribution and explore how this is influenced by the atmospheric circulation shaped by orography, using the Helmholtz wind decomposition and meridional mass streamfunction. The model included Earth-like orography and the substellar point was located over the Pacific Ocean. For such a scenario, our analysis reveals a North-South asymmetry in the simulated O$_3$ distribution. The O$_3$ concentration is highest below 10 hPa level (below $\sim$30 km) near the South pole. This asymmetry results from the presence of land masses on the night side that cause drag in near-surface flows and lead to an asymmetric meridional overturning circulation. Catalytic species were roughly symmetrically distributed and were not found to be primary driver for the O$_3$ asymmetry. The total ozone column (TOC) density was higher for TRAPPIST-1e compared to Earth, with 8000 Dobson Units (DU) near the South pole and 2000 DU near the North pole. The results emphasize the sensitivity of O$_3$ to model parameters, illustrating how incorporating Earth-like orography can affect atmospheric dynamics and O$_3$ distribution. This link between surface features and atmospheric dynamics underlines the importance of how changing model parameters used to study exoplanet atmospheres can influence the interpretation of observations.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
How Land-Mass Distribution Influences the Atmospheric Dynamics of Tidally Locked Terrestrial Exoplanets
Authors:
F. Sainsbury-Martinez,
C. Walsh,
G. J. Cooke,
D. R. Marsh
Abstract:
Interpretation of the ongoing efforts to simulate the atmospheres of potentially-habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally-locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terres…
▽ More
Interpretation of the ongoing efforts to simulate the atmospheres of potentially-habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally-locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terrestrial planets: the inclusion of an Earth-like land-ocean distribution with orography. To do this we use the Earth-system model WACCM6/CESM2 to run a pair of TRAPPIST-1e models with N$_2$-O$_2$ atmospheres and with the sub-stellar point fixed over either land or ocean. The presence of orography shapes atmospheric transport, and in the case of Earth-like orography, breaks the symmetry between the northern and southern hemispheres which was previously found in slab ocean models. For example, peak zonal jet speeds in the southern hemisphere are $50\rightarrow100\%$ faster than similar jets in the northern hemisphere. This also affects the meridional circulation, transporting equatorial material towards the south-pole. As a result we also find significant changes in the atmospheric chemistry, including the accumulation of potentially lethal quantities of ozone at both the south pole and the surface. Future studies which investigate the effects of land-mass distribution on the dynamics of exoplanetary atmospheres should pay close attention to both the day-side land-fraction as well as the orography of the land. Simply modelling a flat land-mass will not give a complete picture of its dynamical impact.
△ Less
Submitted 2 July, 2024; v1 submitted 1 July, 2024;
originally announced July 2024.
-
Lethal surface ozone concentrations are possible on habitable zone exoplanets
Authors:
G. J. Cooke,
D. R. Marsh,
C. Walsh,
F. Sainsbury-Martinez
Abstract:
Ozone ($\textrm{O}_3$) is important for the survival of life on Earth because it shields the surface from ionising ultraviolet (UV) radiation. However, the existence of $\textrm{O}_3$ in Earth's atmosphere is not always beneficial. Resulting from anthropogenic activity, $\textrm{O}_3$ exists as a biologically harmful pollutant at the surface when it forms in the presence of sunlight and other poll…
▽ More
Ozone ($\textrm{O}_3$) is important for the survival of life on Earth because it shields the surface from ionising ultraviolet (UV) radiation. However, the existence of $\textrm{O}_3$ in Earth's atmosphere is not always beneficial. Resulting from anthropogenic activity, $\textrm{O}_3$ exists as a biologically harmful pollutant at the surface when it forms in the presence of sunlight and other pollutants. As a strong oxidiser, $\textrm{O}_3$ can be lethal to several different organisms; thus, when assessing the potential habitability of an exoplanet, a key part is determining whether toxic gases could be present at its surface. Using the Whole Atmosphere Community Climate Model version 6 (WACCM6; a three-dimensional chemistry-climate model), twelve atmospheric simulations of the terrestrial exoplanet TRAPPIST-1 e are performed with a variety of $\textrm{O}_2$ concentrations and assuming two different stellar spectra proposed in the literature. Four atmospheric simulations of the exoplanet Proxima Centauri b are also included. Some scenarios for both exoplanets exhibit time-averaged surface $\textrm{O}_3$ mixing ratios exceeding harmful levels of 40 ppbv, with 2200 ppbv the maximum concentration found in the cases simulated. These concentrations are toxic and can be fatal to most life on Earth. In other scenarios $\textrm{O}_3$ remains under harmful limits over a significant fraction of the surface, despite there being present regions which may prove inhospitable. In the case that $\textrm{O}_3$ is detected in a terrestrial exoplanet's atmosphere, determining the surface concentration is an important step when evaluating a planet's habitability.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
The Impact of Cometary 'impacts' on the Chemistry, Climate, and Spectra of Hot Jupiter Atmospheres
Authors:
Felix Sainsbury-Martinez,
Catherine Walsh
Abstract:
Impacts from icy and rocky bodies have helped shape the composition of solar system objects, for example the Earth-Moon system, or the recent impact of comet Shoemaker-Levy 9 with Jupiter. It is likely that such impacts also shape the composition of exoplanetary systems. Here we investigate how cometary impacts might affect the atmospheric composition/chemistry of hot Jupiters, which are prime tar…
▽ More
Impacts from icy and rocky bodies have helped shape the composition of solar system objects, for example the Earth-Moon system, or the recent impact of comet Shoemaker-Levy 9 with Jupiter. It is likely that such impacts also shape the composition of exoplanetary systems. Here we investigate how cometary impacts might affect the atmospheric composition/chemistry of hot Jupiters, which are prime targets for characterisation. We introduce a parametrised cometary impact model that includes thermal ablation and pressure driven breakup, which we couple with the 1D `radiative-convective' atmospheric model ATMO, including disequilibrium chemistry. We use this model to investigate a wide range of impactor masses and compositions, including those based on observations of Solar System comets, and interstellar ices (with JWST). We find that even a small impactor (R = 2.5 km) can lead to significant short-term changes in the atmospheric chemistry, including a factor $>10$ enhancement in H$_2$O, CO, CO$_2$ abundances, and atmospheric opacity more generally, and the near complete removal of observable hydrocarbons, such as CH$_4$, from the upper atmosphere. These effects scale with the change in atmospheric C/O ratio and metallicity. Potentially observable changes are possible for a body that has undergone significant/continuous bombardment, such that the global atmospheric chemistry has been impacted. Our works reveals that cometary impacts can significantly alter or pollute the atmospheric composition/chemistry of hot Jupiters. These changes have the potential to mute/break the proposed link between atmospheric C/O ratio and planet formation location relative to key snowlines in the natal protoplanetary disc.
△ Less
Submitted 8 February, 2024;
originally announced February 2024.
-
Characterising The Atmospheric Dynamics Of HD209458b-like Hot Jupiters Using AI Driven Image Recognition/Categorisation
Authors:
F. Sainsbury-Martinez,
P. Tremblin,
M. Mancip,
S. Donfack,
E. Honore,
M. Bourenane
Abstract:
In-order to understand the results of recent observations of exoplanets, models have become increasingly complex. Unfortunately this increases both the computational cost and output size of said models. We intend to explore if AI-image-recognition can alleviate this burden. We used DYNAMICO to run a series of HD209458-like models with different orbital-radii. Training data for a number of features…
▽ More
In-order to understand the results of recent observations of exoplanets, models have become increasingly complex. Unfortunately this increases both the computational cost and output size of said models. We intend to explore if AI-image-recognition can alleviate this burden. We used DYNAMICO to run a series of HD209458-like models with different orbital-radii. Training data for a number of features of interest was selected from the initial outputs of these models. This was used to train a pair of multi-categorisation convolutional-neural-networks (CNN), which we applied to our outer-atmosphere-equilibrated models. The features detected by our CNNs revealed that our models fall into two regimes: models with a shorter orbital-radii exhibit significant global mixing which shapes the entire atmospheres dynamics. Whereas, models with longer orbital-radii exhibit negligible mixing except at mid-pressures. Here, the initial non-detection of any trained features revealed a surprise: a night-side hot-spot. Analysis suggests that this occurs when rotational influence is sufficiently weak that divergent flows from the day-side to the night-side dominate over rotational-driven transport, such as the equatorial jet. We suggest that image-classification may play an important role in future, computational, atmospheric studies. However special care must be paid to the data feed into the model, from the colourmap, to training the CNN on features with enough breadth and complexity that the CNN can learn to detect them. However, by using preliminary-studies and prior-models, this should be more than achievable for future exascale calculations, allowing for a significant reduction in future workloads and computational resources.
△ Less
Submitted 19 September, 2023;
originally announced September 2023.
-
Evidence of Radius Inflation in Radiative GCM Models of WASP-76b due to the Advection of Potential Temperature
Authors:
Felix Sainsbury-Martinez,
Pascal Tremblin,
Aaron David Schneider,
Ludmila Carone,
Isabelle Baraffe,
Gilles Chabrier,
Christiane Helling,
Leen Decin,
Uffe Gråe Jørgensen
Abstract:
Understanding the discrepancy between the radii of observed hot Jupiters and standard 'radiative-convective' models remains a hotly debated topic in the exoplanet community. One mechanism which has been proposed to bridge this gap, and which has recently come under scrutiny, is the vertical advection of potential temperature from the irradiated outer atmosphere deep into the interior, heating the…
▽ More
Understanding the discrepancy between the radii of observed hot Jupiters and standard 'radiative-convective' models remains a hotly debated topic in the exoplanet community. One mechanism which has been proposed to bridge this gap, and which has recently come under scrutiny, is the vertical advection of potential temperature from the irradiated outer atmosphere deep into the interior, heating the deep, unirradiated, atmosphere, warming the internal adiabat, and resulting in radius inflation. Specifically, a recent study which explored the atmosphere of WASP-76b using a 3D, non-grey, GCM suggested that their models lacked radius inflation, and hence any vertical enthalpy advection. Here we perform additional analysis of these, and related, models, focusing on an explicit analysis of vertical enthalpy transport and the resulting heating of the deep atmosphere compared with 1D models. Our results indicate that, after any evolution linked with initialisation, all the WASP-76b models considered here exhibit significant vertical enthalpy transport, heating the deep atmosphere significantly when compared with standard 1D models. Furthermore, comparison of a long time-scale (and hence near steady-state) model with a Jupiter-like internal-structure model suggests not only strong radius-inflation, but also that the model radius, $1.98 \mathrm{R_{J}}$, may be comparable with observations ($1.83\pm0.06 \mathrm{R_{J}}$). We thus conclude that the vertical advection of potential temperature alone is enough to explain the radius inflation of WASP-76b, and potentially other irradiated gas giants, albeit with the proviso that the exact strength of the vertical advection remains sensitive to model parameters, such as the inclusion of deep atmospheric drag.
△ Less
Submitted 21 June, 2023;
originally announced June 2023.
-
Exploring deep and hot adiabats as a potential solution to the radius inflation problem in brown dwarfs: Long-timescale models of the deep atmospheres of KELT-1b, Kepler-13Ab, and SDSS1411B
Authors:
F. Sainsbury-Martinez,
S. L. Casewell,
J. D. Lothringer,
M. W. Phillips,
P. Tremblin
Abstract:
The anomalously large radii of highly-irradiated gaseous exoplanets has long been a mystery. One mechanism suggested as a solution for hot Jupiters is the heating of the deep atmosphere via the vertical advection of potential temperature resulting in an increased internal entropy. Here we intend to explore if this mechanism can also explain the observed brown dwarf radii trend: a general increase…
▽ More
The anomalously large radii of highly-irradiated gaseous exoplanets has long been a mystery. One mechanism suggested as a solution for hot Jupiters is the heating of the deep atmosphere via the vertical advection of potential temperature resulting in an increased internal entropy. Here we intend to explore if this mechanism can also explain the observed brown dwarf radii trend: a general increase in radius with irradiation, with an exception for highly-irradiated brown dwarfs orbiting white dwarfs. We use a 3D GCM, DYNAMICO, to run a series of long-timescale models of the atmospheres of Kepler-13Ab, KELT-1b, and SDSS1411B. These models allow us to explore not only if a stable advective adiabat can develop, but also the associated dynamics. We find that our models fall into two distinct regimes: Kepler-13Ab and KELT-1b both show signs of significant deep heating and hence maintain adiabats that are hotter than 1D models predict. On the other hand, SDSS1411B exhibits a much weaker downward heating profile which not only struggles to heat the interior under ideal conditions, but is highly sensitive to the presence of deep radiative dynamics. We find that the vertical advection of potential temperature by large-scale atmospheric circulations represents a robust mechanism to explain the trend of increasing inflation with irradiation, including the exception for highly irradiated brown dwarfs orbiting white dwarfs. This can be understood as occurring due to the role that increasing rotational influence plays on mid- to-high latitude advective dynamics. Furthermore, when paired with a suitable parametrisation of the outer atmosphere irradiation profile, this mechanism alone could potentially provide a complete explanation for the observed levels of inflation in our brown dwarfs.
△ Less
Submitted 27 September, 2021; v1 submitted 17 September, 2021;
originally announced September 2021.
-
Idealised simulations of the deep atmosphere of hot jupiters: Deep, hot, adiabats as a robust solution to the radius inflation problem
Authors:
F. Sainsbury-Martinez,
P. Wang,
S. Fromang,
P. Tremblin,
T. Dubos,
Y. Meurdesoif,
A. Spiga,
J. Leconte,
I. Baraffe,
G. Chabrier,
N. Mayne,
B. Drummond,
F. Debras
Abstract:
Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scena…
▽ More
Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
Radial velocity eclipse mapping of exoplanets
Authors:
Nikolay Nikolov,
Felix Sainsbury-Martinez
Abstract:
Planetary rotation rates and obliquities provide information regarding the history of planet formation, but have not yet been measured for evolved extrasolar planets. Here we investigate the theoretical and observational perspective of the Rossiter-McLauglin effect during secondary eclipse (RMse) ingress and egress for transiting exoplanets. Near secondary eclipse, when the planet passes behind th…
▽ More
Planetary rotation rates and obliquities provide information regarding the history of planet formation, but have not yet been measured for evolved extrasolar planets. Here we investigate the theoretical and observational perspective of the Rossiter-McLauglin effect during secondary eclipse (RMse) ingress and egress for transiting exoplanets. Near secondary eclipse, when the planet passes behind the parent star, the star sequentially obscures light from the approaching and receding parts of the rotating planetary surface. The temporal block of light emerging from the approaching (blue-shifted) or receding (red-shifted) parts of the planet causes a temporal distortion in the planet's spectral line profiles resulting in an anomaly in the planet's radial velocity curve. We demonstrate that the shape and the ratio of the ingress-to-egress radial velocity amplitudes depends on the planetary rotational rate, axial tilt and impact factor (i.e. sky-projected planet spin-orbital alignment). In addition, line asymmetries originating from different layers in the atmosphere of the planet could provide information regarding zonal atmospheric winds and constraints on the hot spot shape for giant irradiated exoplanets. The effect is expected to be most-pronounced at near-infrared wavelengths, where the planet-to-star contrasts are large. We create synthetic near-infrared, high-dispersion spectroscopic data and demonstrate how the sky-projected spin axis orientation and equatorial velocity of the planet can be estimated. We conclude that the RMse effect could be a powerful method to measure exoplanet spins.
△ Less
Submitted 16 June, 2015;
originally announced June 2015.