Atmospheric convection plays a key role in the climate of tidally locked planets:

insights from high-resolution simulations

Department of Meteorology | University of Reading | 27 April 2020

Denis Sergeev

University of Exeter, UK

This presentation is available at slides.com/denissergeev/2020-04-27-uor-met

Exeter Exoplanet Theory Group

Group leader: Nathan Mayne

Staff: Eric Hebrard, F. Hugo Lambert

Postdocs: Duncan Christie, Denis Sergeev,   Maria Zamyatina

PhD students: Michelle Bieger, Jake Eager, Robert Ridgway

IT support: Krisztian Kohary

UK Met Office collaborators: Ian Boutle, James Manners, Benjamin Drummond, Stefan Lines

Talk outline

  1. Challenges and importance of simulating convection on exoplanets
  2. Global simulations of Trappist-1e and Proxima Centauri b
  3. Regional high-resolution simulations and potential implications for the global climate

Introduction

Exoplanets

  • >4000 exoplanets discovered
  • Most detected are gas giants (like Jupiter)
  • More and more rocky planets are detected

... are awesome!

Credit: ESO

Terrestrial tidally locked planets

  • Habitable ~ have surface liquid water 
  • M-dwarf are relatively cool
  • Potentially habitable planets around M-dwarfs have to orbit closer
  • Strong tidal forces
  • Planet can become tidally locked
  • Orbital period = rotation period
  • Permanent day-side and night-side

 

Credit: Engine House Animation Studio

Credit: ESA/Hubble

Climate of Earth-like tidally locked exoplanets

Climate of Earth-like tidally locked exoplanets

Climate of Earth-like tidally locked exoplanets

Challenges in convection modelling

  • Wide scale separation between convection and planetary-scale circulation 
  • Resolving convection is computationally expensive
  • Convection is subgrid-scale and parameterised
  • All existing parameterisations are tested against observations on Earth and still are a big uncertainty for future climate projections
  • No in-situ measurements on exoplanets

Credit: NASA

Modelling framework

Met Office Unified Model

  • Full deep-atmosphere non-hydrostatic Navier-Stokes equations
  • Full suite of parameterisations
  • Seamless modelling: the same model is used for global climate & local high-res modelling
  • Adapted for exoplanet atmospheres

Mayne+ (2014)

Drummond+ (2018)

Lines+ (2019)

Boutle+ (2017)

Adaptation of the dynamical core

Climate modelling of a terrestrial exoplanet

Chemistry on Hot Jupiters 

Clouds on Hot Jupiters 

  • Orbital parameters of Trappist-1e and Proxima Centauri b
  • Planets are tidally locked to M-dwarf stars
  • \(N_2\)-dominated atmosphere with traces of \(CO_2\) and \(H_2O\)
  • Mean surface pressure of 1 bar
  • Aquaplanet regime (a slab ocean)
  • Two modes
    • Global coarse-resolution (\(2^\circ\times 2.5^\circ\)) model with parameterised convection
    • Regional high-resolution (4 km) model with explicit convection

Credit: NASA/JPL-Caltech

Credit: ESA

Trappist-1e

Proxima b

Representation of convection in our global setup

  • Standard operational scheme in the UM
  • Most sophisticated class of convection schemes
  • Involves a lot of "tuning" parameters
  • Convective adjustment to a reference state
  • Only 2 free parameters
  • Similar to Betts-Miller convection scheme

NoCnvPm

Adjust

MassFlux (Control)

  • Convection parameterisation switched off

Mean global climate and its sensitivity to the convection parameterisation

Global simulation of Trappist-1e

Shown are:

  • Surface temperature
  • Cloud condensate
  • Upper troposphere wind vectors

Surface temperature (\(T_s\))

  • \(T_s\) day-night contrast larger on Proxima b than on Trappist-1e

MassFlux

Surface temperature (\(T_s\))

  • \(T_s\) day-night contrast larger on Proxima b than on Trappist-1e

MassFlux

  • \(T_s\) increases globally, especially in cold traps on the night side
  • \(\Delta T_{s,dn}\) reduces

Adjust

Adjust minus MassFlux

Adjust minus MassFlux

Surface temperature (\(T_s\))

  • \(T_s\) day-night contrast larger on Proxima b than on Trappist-1e

MassFlux

  • \(T_s\) increases globally, especially in cold traps on the night side
  • \(\Delta T_{s,dn}\) reduces

NoCnvPm

Adjust

  • Small global-mean change
  • Different response on Trappist1-e vs Proxima b

NoCnvPm minus MassFlux

NoCnvPm minus MassFlux

Adjust minus MassFlux

Adjust minus MassFlux

Surface temperature (\(T_s\)) and free troposphere winds

  • \(T_s\) day-night contrast larger on Proxima b than on Trappist-1e
  • \(T_s\) increases globally, especially in cold traps on the night side
  • \(\Delta T_{s,dn}\) reduces

NoCnvPm

Adjust

MassFlux

  • Small global-mean change
  • Different response on Trappist1-e vs Proxima b

NoCnvPm minus MassFlux

NoCnvPm minus MassFlux

Adjust minus MassFlux

Adjust minus MassFlux

Circulation regime (eddy components)

  • Planetary-scale equatorial Rossby wave-like pattern

MassFlux

Circulation regime (eddy components)

  • Planetary-scale equatorial Rossby wave-like pattern

MassFlux

Circulation regime (eddy components)

  • Planetary-scale equatorial Rossby wave-like pattern

MassFlux

  • Almost a mirror-opposite of the MassFlux
  • Dominated by extratropical baroclinic gyres

Adjust

Circulation regime (eddy components)

  • Planetary-scale equatorial Rossby wave-like pattern
  • Almost a mirror-opposite of the MassFlux
  • Dominated by extratropical baroclinic gyres

NoCnvPm

Adjust

MassFlux

  • Broadly similar response as in Adjust

Circulation regime (eddy components)

  • Planetary-scale equatorial Rossby wave-like pattern
  • Almost a mirror-opposite of the MassFlux
  • Dominated by extratropical baroclinic gyres

NoCnvPm

Adjust

MassFlux

  • Broadly similar response as in Adjust
  • Different shape of latent heating on the day side
  • Change in baroclinicity

Clouds and precipitation

MassFlux

Clouds and precipitation

MassFlux

Clouds and precipitation

MassFlux

Clouds and precipitation

  • Large reduction of clouds on the day side
  • More clouds in polar regions of Trappist-1e

Adjust

MassFlux

Clouds and precipitation

  • Large reduction of clouds on the day side
  • More clouds in polar regions of Trappist-1e

NoCnvPm

Adjust

  • Day side: closer to MassFlux
  • Night side: closer to Adjust

MassFlux

Summary of global simulations: the effects of convection parameterisation

  • The global climate of TL exoplanets is altered (and not just the day side!)
  • The effect is planet-dependent
  • Caution should be taken when interpreting results of GCMs with one parameterisation or another.

Extra

Circulation regime

What's happening in the Trappist-1e case?

  • Different shape of latent heating on the day side
  • Change in baroclinicity
\lambda_{Rh} = \frac{\text{Rhines length}}{\text{planet's radius}}
\lambda_{Ro} = \frac{ \text{equatorial radius of deformation} }{ \text{planet's radius} }

Differences in circulation regimes

  • Leconte+ (2013)
  • Haqq-Misra+ (2017)

Energy redistribution

  • Using the horizontal divergence of the MSE flux, \(\int_0^{z_{top}} \rho \nabla\cdot\left(\vec u (c_p T + g z + L_v q) \right)dz\)
  • The largest difference is in the moist component (\(Lq\))
  • Adjust results in a more intense energy transport

Vertical profiles of temperature and humidity

Vertical cross-sections of heat fluxes at \(90^\circ E\)

High-resolution simulations

"global" (MassFlux)

"HighRes"

Nested grid setup

Spatial variability: top-of-atmosphere outgoing LW radiation

Spatial variability: top-of-atmosphere outgoing LW radiation

Spatial variability: top-of-atmosphere outgoing LW radiation

Spatial variability: top-of-atmosphere outgoing LW radiation

Histograms of instantaneous TOA OLR

Spatial variability: top-of-atmosphere outgoing LW radiation

Histograms of instantaneous TOA OLR

Vertical cloud structure (averaged profiles)

Vertical cloud structure (averaged profiles)

Vertical cloud structure (averaged profiles)

  • Biases in liquid water clouds
  • Thicker layer of mixed-phase clouds
  • Smaller bias in ice clouds

Latent heating (averaged profiles)

  • Latent heating is stronger in the global model, compared to HighRes
  • Impact on the upper-level flow divergence?

Extra

Global impact of resolved convection

  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?
  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?
  1. Run an ensemble of global experiments with perturbed convection parameters
  1. Run an ensemble of global experiments with perturbed convection parameters
  2. Integrate a metric of convection strength over the substellar region
  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?
  1. Run an ensemble of global experiments with perturbed convection parameters
  2. Integrate a metric of convection strength over the substellar region
  3. Correlate with a metric of day-night heat redistribution
  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?
  1. Run an ensemble of global experiments with perturbed convection parameters
  2. Integrate a metric of convection strength over the substellar region
  3. Correlate with a metric of day-night heat redistribution
  4. Extrapolate for the convection metric from HighRes simulations
  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?
  1. Run an ensemble of global experiments with perturbed convection parameters
  2. Integrate a metric of convection strength over the substellar region
  3. Correlate with a metric of day-night heat redistribution
  4. Extrapolate for the convection metric from HighRes simulations
  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?
  1. Run an ensemble of global experiments with perturbed convection parameters
  2. Integrate a metric of convection strength over the substellar region
  3. Correlate with a metric of day-night heat redistribution
  4. Extrapolate for the convection metric from HighRes simulations
  • Caveat: one-way nesting of the model
  • How to assess the impact of HighRes simulations on the global climate?

Summary

  • Representation of convection is important for the global climate of exoplanets
  • Using a simpler adjustment scheme leads to
    • on the day side: reduction of clouds, decrease in albedo (by ~60%)
    • on the night side: warming of 20-30 K
    • reduction in the day-night temperature contrast
    • change in the global circulation
  • The sensitivity to convection scheme is planet-dependent
  • The choice of convection scheme does not push the simulated climate out of habitable state
  • A potential global convection-permitting simulation may enhance the day-night temperature contrast

Future work

  • Extending this analysis to broader range of planetary and stellar configurations
  • Exploring the sensitivity to the lower boundary condition, e.g. a presence of a continent on the day side
  • Coupling the high-resolution simulation with a chemical cycle and space weather (e.g. stellar flares)
  • Global resolved convection on tidally locked exoplanets (LFRic project at the Met Office)

Thank you!

My twitter: meteodenny

Accepted preprint: arxiv.org/abs/2004.03007

Seminar talk at the Department of Meteorology, University of Reading

By Denis Sergeev

Seminar talk at the Department of Meteorology, University of Reading

  • 340