Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic
GMD | Articles | Volume 13, issue 2
Geosci. Model Dev., 13, 707–716, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
Geosci. Model Dev., 13, 707–716, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Model experiment description paper 21 Feb 2020

Model experiment description paper | 21 Feb 2020

TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI): motivations and protocol version 1.0

TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI): motivations and protocol version 1.0
Thomas J. Fauchez1,2,3, Martin Turbet4,5, Eric T. Wolf6,7, Ian Boutle8, Michael J. Way9,3, Anthony D. Del Genio9, Nathan J. Mayne10, Konstantinos Tsigaridis9,12, Ravi K. Kopparapu2,3, Jun Yang11, Francois Forget4, Avi Mandell2,3, and Shawn D. Domagal Goldman2,3 Thomas J. Fauchez et al.
  • 1Goddard Earth Sciences Technology and Research (GESTAR), Universities Space Research Association (USRA), Columbia, MD, USA
  • 2NASA Goddard Space Flight Center, Greenbelt, MD, USA
  • 3GSFC Sellers Exoplanet Environments Collaboration, Greenbelt, MD, USA
  • 4Laboratoire de Méteorologie Dynamique, IPSL, Sorbonne Universités, UPMC Univ. Paris 06, CNRS, 4 Place Jussieu, 75005 Paris, France
  • 5Observatoire Astronomique de l'Université de Genève, Université de Genève, Chemin des Maillettes 51, 1290 Versoix, Switzerland
  • 6Laboratory for Atmospheric and Space Physics, Department of Atmospheric and Oceanic Sciences, University of Colorado Boulder, Boulder, CO, USA
  • 7Virtual Planetary Laboratory, NASA Astrobiology Institute, Seattle, WA, USA
  • 8Met Office, Exeter, UK
  • 9NASA Goddard Institute for Space Studies, New York, NY 10025, USA
  • 10University of Exeter, Exeter, UK
  • 11Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871, China
  • 12Center for Climate Systems Research, Columbia University, New York, NY, USA

Correspondence: Thomas J. Fauchez (

Back to toptop

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

1 Introduction
Back to toptop

M dwarfs are the most common type of stars in our galaxy, and rocky exoplanets orbiting M-dwarf stars will likely be the first to be characterized with upcoming astronomical facilities such as the James Webb Space Telescope (JWST). Ultra-cool dwarfs (T<2700 K) are a substellar class of late M dwarfs and represent nearly 15 % of astronomical objects in the stellar neighborhood of the Sun (Cantrell et al.2013). Their smaller size compared to other stellar types allows for easier detection of rocky exoplanets in close orbits, and this potential was recently realized by the discovery of the TRAPPIST-1 system (Gillon et al.2016, 2017). Located about 12 pc (parsec) away, TRAPPIST-1 has seven known planets, and is one of the most promising rocky-planet systems for follow-up observations due to the depths of the transit signals (Gillon et al.2017; Luger et al.2017). Transit timing variation (TTV) measurements of the TRAPPIST-1 planets suggest a terrestrial composition likely enriched in volatiles, and possibly water (Grimm et al.2018). Also, it has been found that three planets (TRAPPIST-1e, f and g) are in the habitable zone (HZ; Kopparapu et al.2013) where surface temperatures would allow surface water to exist (Gillon et al.2017; Wolf2017, 2018; Turbet et al.2018).

TRAPPIST-1 is an active M-dwarf star (O'Malley-James and Kaltenegger2017; Wheatley et al.2017; Vida and Roettenbacher2018) which offers an environment very hostile to the survival of planetary atmospheres. However, Bolmont et al. (2017) and Bourrier et al. (2017) argued that depending on their initial water contents, the TRAPPIST-1 planets could have retained some water presently. Assuming that this water has remained in sufficient quantity, TRAPPIST-1e may be able to maintain habitable conditions (locally or globally around the planet) through a very large set of atmospheric configurations (Wolf2017; Turbet et al.2018; Grootel et al.2018, and references therein). The first attempt to characterize those planets through transmission spectroscopy has been conducted by de Wit et al. (2016, 2018) using the Hubble Space Telescope (HST) for the six innermost planets. Their analysis suggests that the TRAPPIST-1 planets do not have a cloud- or haze-free H2-dominated atmosphere, and that a large set of high-mean-molecular-weight atmospheres are possible, such as thick N2-, O2-, H2O-, CO2-, or CH4-dominated atmospheres. Using laboratory measurements and models Moran et al. (2018) have also shown that H2 dominated atmospheres with cloud/haze can also be ruled out. Note that the uncertainties of these HST observations were very large, on the order of hundreds of parts per million (ppm), and further investigations with future facilities such as the JWST (Barstow and Irwin2016; Morley et al.2017) will be needed to determine the nature of atmospheres heavier than hydrogen.

Upstream of future JWST characterization of TRAPPIST-1e, it is important to derive constraints on its possible atmosphere to serve as a guideline for the observations. For this purpose, 3-D global climate models (GCMs) are the most advanced tools (Wolf et al.2019). However, GCMs are very complex models and their outputs can vary from one model to another for a variety of reasons. GCM intercomparisons have been widely used by the Earth science community. For instance the Coupled Model Intercomparison Project (CMIP), initiated in 1995 and currently at version 6 (Eyring et al.2016), focuses on the differences in GCM responses to forcings for anthropogenic climate change. While exoplanets receive considerable attention from climate modelers, and atmospheric data from Earth-like worlds may be imminent, to our knowledge only one intercomparison of planetary GCMs has been published (Yang et al.2019). They found significant differences in global surface temperature between the models for planets around M-dwarf stars due to differences in atmospheric dynamics, clouds and radiative transfer. However, Yang et al. (2019) concerns planets near the inner edge of the HZ and focuses on highly idealized planetary configurations. Note that another model intercomparison have been run for the exoplanet community: the Palaeoclimate and Terrestrial Exoplanet Radiative Transfer Model Intercomparison Project (PALAEOTRIP). The protocol of this experiment is described in Goldblatt et al. (2017) and aims to compare a large variety of radiation codes used for paleoclimate or exoplanet sciences and to identify the limit conditions for which each model can produce accurate results. The timeline and information about PALAEOTRIP can be found at (last access: 8 February 2020). The motivation behind the TRAPPIST Habitable Atmosphere Intercomparison (THAI) is to highlight differences among GCM simulations of a confirmed exoplanet, TRAPPIST-1e, that is potentially characterizable in the near term (with JWST or ground-based facilities) and to evaluate how these differences may impact our interpretations of retrievals of its atmospheric properties from delivered observables. Our objective is also to provide a clear protocol intended for other GCMs to join the intercomparison, which is therefore not only limited to the GCMs presented in this paper. Results of the intercomparison will be presented in a following paper. In this paper, the motivations, including a presentation of TRAPPIST-1e and of the GCMs, are presented in Sect. 2. In Sect. 3 we present the THAI protocol describing all the parameters to be set up in the GCM. In Sect. 4, we list the model parameters to be provided in order for a given model to be comparable to other GCM simulations. A summary is given in Sect. 5.

2 TRAPPIST-1e climate simulation and motivations
Back to toptop

2.1 Motivations for a planetary GCM intercomparison

Global climate models (GCMs) are 3-dimensional numerical models designed to represent physical processes at play in planetary atmospheres and surfaces. They are the most sophisticated way to model the atmospheres and oceans of real planets. GCMs can be seen as a complex network of 1-D time-marching climate models connected together through a dynamical core (see description below). Each 1-D column contains physical parameterizations for radiative transfer, convection, boundary layer processes, cloud macroscale and microscale physics, aerosols, precipitation, surface snow and sea ice accumulation, and other processes at varying levels of complexity.

The motivation behind this experimental protocol is to evaluate how some of the differences between the models can impact the assessment of the planet's habitability and its observables through transmission spectroscopy and thermal phase curves with upcoming observatories such as JWST. The intercomparison protocol was designed to evaluate three possible sources of differences between the models listed.

  1. The dynamical core. The dynamical core is a numerical solver of the hydrodynamic equations on the (rotating) planetary sphere. It calculates the winds that transport atmospheric gases, clouds, aerosols, sensible and latent heat, and momentum from one atmospheric column to another.

  2. The radiative transfer. Each model has its own radiative transfer working assumptions and may use different spectroscopic databases and even different versions of the same spectroscopic database (e.g., HITRAN), collision-induced absorption (CIA), line-by-line versus correlated-k distribution (Lacis and Oinas1991), line cutoff, spectral resolution, etc.

  3. The moist physics. The treatment of water in all of its thermodynamic phases is critical for the simulation of habitable planets. In particular, cloud and convection process are a significant source of differences between climate models, and these differences are often exacerbated when modeling exoplanets around M-dwarf stars (Yang et al.2014, 2019).

Note that a particular emphasis will be given to the differences of cloud properties between the models, because they may have a large impact on the strength of the spectral signatures simulated by current radiative transfer tools (Fauchez et al.2019). Yet a sufficient understanding of 3-D cloud fields is needed to provide realistic observational constraints to observers. It is therefore crucial to address these potential differences between the GCMs.

Four GCMs (in their planetary version) are initially included with THAI:

  1. the Laboratoire de Météorologie Dynamique – Generic model (LMDG; Wordsworth et al.2011, a review paper on the model is currently under preparation),

  2. the Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D, Planet 1.0 version derived from the NASA GISS Model E; Way et al.2017),

  3. the Exoplanet Community Atmospheric Model (ExoCAM, available on GitHub:, last access: 8 February 2020, available from the National Center for Atmospheric Research (NCAR):, last access: 8 February 2020, derived from the CAM4 NCAR model; Neale2010),

  4. the Met Office Unified Model (UM; Mayne et al.2014; Boutle et al.2017).

By publishing our protocols in advance of the intercomparison work, we hope that other teams will also use this protocol to compare their own GCM with the four GCMs of this study.

2.2 The TRAPPIST-1e benchmark

TRAPPIST-1e is up to now one of the best habitable planet candidates for atmospheric characterization through transmission spectroscopy with the JWST. Therefore, it is also an obvious candidate for an experimental protocol for GCM intercomparison. In Table 1 we summarize the TRAPPIST-1e parameters used in the THAI project based on Grimm et al. (2018).

Table 1TRAPPIST-1 stellar spectrum and TRAPPIST-1e planetary parameters from Grimm et al. (2018).

Download Print Version | Download XLSX

3 The THAI protocol
Back to toptop

3.1 Atmospheric configurations

For THAI, we have chosen a set of four planetary configurations with increasing complexity. We have chosen to start with benchmark cases of dryland planets with N2- and CO2-dominated atmospheres, which will allow us to assess atmospheric dynamical plus boundary layer processes and CO2 radiative transfer differences. Next we conduct aqua-planet simulations of N2- and CO2-dominated atmospheres, providing characteristic cold and warm, respectively, habitable states for TRAPPIST-1e. By gradually increasing the complexity of our simulations, we hope to be able to parse out meaningful differences between atmospheric dynamical plus boundary layer, radiative transfer, and moist physical processes. The motivation for each of these cases is described below.

  • Benchmark case 1 (Ben1). In this case, composed of 1 bar of N2 and 400 ppm of CO2, the purpose is to test the differences of the planetary boundary layer (PBL) schemes, the dynamical core and the associated heat redistribution between the different models. Note that N2N2 CIA should be included.

  • Benchmark case 2 (Ben2). In this case, composed of 1 bar of CO2, we test the PBL schemes and dynamical core differences as well as the CO2 radiative transfer.

  • Habitable case 1 (Hab1). In this case, composed of a modern Earth-like atmosphere of 1 bar of N2 and 400 ppm of CO2, the dynamical core, the clouds and atmospheric processes are tested together. It is also the most widespread benchmark for habitable planets in the literature (Barstow and Irwin2016; Morley et al.2017; Lincowski et al.2018).

  • Habitable case 2 (Hab2). In this case, composed of 1 bar of CO2, the dynamical core, CO2 radiative transfer assumption and the clouds and atmospheric processes are tested. This case is likely representative of early Earth (during the Hadean epoch), early Venus, and early Mars at a time when Martian valley networks and lakes were formed (Haberle et al.2017; Kite2019).

We have therefore two lands planets (Ben1 and Ben2) and two aqua planets (Hab1 and Hab2). Note that Ben1 and Hab1 share the same atmospheric composition of 1 bar of N2 and 400 ppm of CO2, and Ben2 and Hab2 share the same atmospheric composition of 1 bar of CO2. The protocol is therefore symmetrical with respect to the atmospheric composition between the land planets and aqua planets.

Figure 1Surface contours for surface temperature, thermal emitted radiation (top of atmosphere, TOA) and reflected stellar radiation (TOA) for Hab1 simulated by the four GCMs: the UK Met Office United Model (UM), the Laboratoire Météorologie Dynamique Generic model (LMDG), the Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) and the National Center for Atmospheric Research Community Atmosphere Model version 4 modified for exoplanets (ExoCAM).


In each case, it is crucial to start each simulation with the same initial conditions. The simplest approach is then to start with an isothermal atmosphere. For THAI, we fixed the initial surface and atmosphere temperature at 300 K. The atmospheric configurations for the two benchmark (dry land) cases and two habitable cases are listed in Table 2, upper block. Note that for Ben2 initial results indicate that some models feature cold trap temperatures on the nightside slightly below the CO2 condensation point at 1 bar (194 K). However, because not all the models include CO2 condensation, it should be disabled in the models that allow it. Ben2 is thus to be viewed as an idealization for the sake of study. Initial results indicate that Hab1 is representative of a cool, largely ice covered world but with liquid water in the substellar region. Hab2 is significantly warmer than Hab1, owing to a strong CO2 greenhouse effect and the water vapor greenhouse feedback, and it is representative of a temperate habitable world. The amount and variability of clouds and the strength of the atmospheric processes should be enhanced, providing a more challenging comparison than in Hab1.

In Fig. 1 we show results from preliminary simulations for case Hab1 conducted with four different GCMs; UM, LMDG, ROCKE-3D and ExoCAM. We show surface contours for surface temperature, thermal emitted radiation (top of atmosphere (TOA)) and reflected stellar radiation (TOA). We can see significant differences in the maximum, minimum and mean values of these parameters between the models. For such a complex atmosphere, it is difficult to disentangle the effects leading to these differences. However, it seems clear that the patterns of thermal emitted and reflected stellar TOA fluxes are strongly influenced by the cloud patterns produced by each respective model. Here we have shown preliminary outputs to demonstrate the feasibility of the described experiments. In-depth analysis of these simulations will be discussed in a following article in preparation.

Table 2THAI experimental protocol.

Download Print Version | Download XLSX

3.2 Surface

The surfaces considered in THAI (Table 2, lower block) are simple. The land planets (Ben1 and Ben2) are covered by sand with a subsurface depth of at least 3 m with a constant albedo of 0.3. The ocean planets (Hab1 and Hab2) are fully covered by a 100 m deep slab (no horizontal heat transport) ocean. The ocean albedo is fixed at 0.06, and the ice and snow albedos are fixed at a constant value of 0.25. Note that the sea ice or snow albedo parameterization is a common source of discrepancy between the models. Some models, like ROCKE-3D, account for the spectral dependence of the sea ice albedo over multiple bands and variations due to snowfall, aging, depth and melt ponds, while other models, such as LMDG, compute the wavelength-dependent albedo of water ice or snow from a simplified albedo spectral law, calibrated to get an ice or snow bolometric albedo of ∼0.25 around an ultra-cool star like TRAPPIST-1 (Joshi and Haberle2012; von Paris et al.2013; Shields et al.2013). Differences in sea ice albedo have been found to have a large impact on planetary climate and habitability (Turbet et al.2018). However, for the sake of this intercomparison, this discrepancy can be easily avoided by fixing the sea ice and snow albedo at a constant bolometric value of 0.25.

Figure 2Globally averaged surface temperature (a) and cloud water column (b) as a function of the number of orbits for Hab1 and Hab2 simulated with ExoCAM (Wolf and Toon2015). Surface temperature for Hab1 and cloud water column for the two cases vary by a couple of tens of percent on a timescale of 10 orbits due to weather patterns.


3.3 Model spatial resolutions and time steps

The model spatial resolution is an important parameter because every process taking place at a subgrid level would be parameterized and those parameterizations often diverge between the models. Similarly, the model time steps control the numerical stability and accuracy. However, the choice of time steps is fundamental to how each model operates under a given parameterization, and arbitrarily fixing these parameters may prevent some models from correctly and fairly performing the intercomparison. In addition, models should be compared using the specifications that they commonly use for exoplanet studies. Therefore, for the sake of the THAI, we do not impose the model spatial resolution nor time steps. However, note that we recommend (but this is not a requirement) the radiative time step (a parameter much more flexible than the others among the models) to be set to 1800 s. This value should provide a good coupling of the radiation with temporal changes to the atmosphere without slowing down the model too much.

We also ask the contributing scientists to disable the subgrid gravity wave parameterizations in their model. Indeed, not all of the models have implemented a gravity wave parameterization, and some have prescribed or predicted gravity wave formation, tuned for Earth topography and meteorology. Therefore, to avoid differences in atmospheric dynamics, especially above the tropopause, we recommend to not include the subgrid gravity wave parameterizations in this intercomparison. Gravity waves whose wavelengths are greater than the model grid are explicitly resolved in the models and do not need to be modified.

Note that under the requirements of the protocol, the atmospheric simulation of TRAPPIST-1e may actually not represent what each individual model can simulate with all their parameterizations fully activated. This is especially true for the sea ice and snow albedo parameterization. Therefore, complementary to the Hab1 case, we propose the Hab1* which should be simulated with the commonly used model parameterizations fully activated. Therefore, only the requirements on the atmospheric composition (1 bar of N2 and 400 ppm of CO2) and the planet and star properties of Table 1 are constrained for Hab1*.

Table 3Instantaneous fields to be output by the GCM. For each diagnostic, the mean value and the standard deviation are computed from data output at the specified frequency and number of orbits for the case. OLR and ASR correspond to outgoing longwave radiation (at TOA) and absorbed shortwave radiation (at TOA), respectively; SW and LW correspond to shortwave and longwave radiation, respectively; CF corresponds to cloud fraction and MMR corresponds to mass mixing ratio. The empty set symbol () represents parameters that are not relevant for a specific case.

Download Print Version | Download XLSX

4 Outputs
Back to toptop

To compare the difference between models of a particular (instantaneous) output variable, both the average and standard deviation over the specified frequency and number of orbits for the case will be computed. Four categories of outputs frequently used in climate simulations have been selected: radiation, surface, atmospheric profiles and clouds. The radiation outputs are the outgoing longwave radiation (OLR) and absorbed shortwave radiation (ASR) for clear and cloudy skies, also commonly known as emitted thermal and absorbed stellar fluxes, respectively, both at the top of the atmosphere (TOA). The surface outputs are the temperature map, the downward total SW flux and net LW flux, and the open ocean fraction (for Hab1/Hab1* and Hab2 only). The atmosphere outputs are the temperature and the U, V and W wind speed profiles. Finally, the cloud outputs for Hab1/Hab1* and Hab2 are the water vapor and cloud condensed water and ice integrated columns; the cloud profiles of the cloud fraction; and the mass mixing ratio for the liquid, ice, and both combined. Also in these two cases, the spatial and temporal variability is much weaker than in Hab1/Hab1* and Hab2. Therefore, to mitigate the amount of data, we choose to only output data for 10 consecutive orbits (with a 6 h output frequency). Concerning Hab1/Hab1* and Hab2, we can see in Fig. 2 that weather patterns modulate the surface temperature and cloud water column of Hab1 on a period nearly equal to 10 orbits. Also Hab2 (hotter than Hab1) has a more efficient heat transport and is therefore more homogeneous in temperature, but the cloud variability is very important. Therefore, more orbits (100) are needed in order to smooth out this variability. A summary of the output parameters is given in Table 3.

All the simulations should have reached radiative equilibrium at TOA at ±1 W m−2. If such a limit can not be achieved by the model, the radiative equilibrium can be established if no discernible trends are observable in the last 10-year average global mean temperature.

To facilitate comparison between each GCM, we ask the contributing scientists to provide their outputs in netCDF format. The contributing scientist will be able to upload their data to a public permanent repository at (last access: 8 February 2020) after requesting an IP address authorization from Thomas Fauchez (

The main objective of THAI is to highlight how differences in atmospheric profiles produced by each GCM are going to impact the predictions of atmosphere detectability and observational constraints for habitable planet targets such as TRAPPIST-1e (Morley et al.2017; Fauchez et al.2019). Therefore, in addition to the parameters of Table 3, we will emphasize the differences between the models in terms of the planet's climate and habitability with a particular attention on the cloud coverage. Also, the objective will be to identify and quantify the differences in the simulated JWST observations, through simulated transmission spectra (in NIRSpec (near-infrared spectrograph) prism and MIRI LRS (mid-infrared instrument, low-resolution spectrometer) ranges) and thermal phase curves (in MIRI LRS range) due to the differences in atmospheric profiles (temperature, pressure and gas mixing ratios) output by each GCM. The planetary spectrum generator (PSG; Villanueva et al.2018) will be used to simulate transmission and emission spectra. The comparison of the spectra for Hab1 and Hab2 cases will therefore highlight the sensitivity of model characteristics to predict transmission spectra of habitable planets.

Note that while additional simulations with a simple Newton cooling model, a 1-D column model, or with cloud radiative effects disabled would help to better understand the differences due to the dynamical cores and cloud physics, they would also dramatically increase the computational time, amount of data and effort. THAI aims to be easily reproducible and not time consuming in order to reach many GCM user groups. The five simulations proposed in THAI should be enough to understand the main differences between the GCMs and their impact on the observables. THAI could also be used as a benchmark for future GCM intercomparisons that specifically aim to understand each differences between the models.

5 Summary
Back to toptop

THAI is an intercomparison project of planetary GCMs focused on the exciting new habitable planet candidate, TRAPPIST-1e. Because rocky exoplanets in the habitable zone of nearby M dwarfs have the highest chance to be the first Earth-size exoplanets to be characterized with future observatories, TRAPPIST-1e is currently the best benchmark we could think of to compare the capability of planetary GCMs. In this first paper we have presented the planet and GCM parameters to be used in this experiment which already has four GCMs included (LMDG, ROCKE-3D, ExoCAM and UM), but we hope more GCMs will join the project. The results of the comparison of these four models will be given in a second paper and a THAI workshop is planned for fall 2020.

Code availability
Back to toptop
Code availability. 

ExoCAM (Wolf and Toon2015) is available on GitHub: (last access: 8 February 2020). The Met Office Unified Model is available for use under license; see (Met Office2020, last access: 8 February 2020). ROCKE-3D is public domain software and available for download for free from (last access: 8 February 2020, NASA Goddard Institute for Space Studies2020a). Annual tutorials for new users take place annually, whose recordings are freely available online at (last access: 8 February 2020b, NASA Goddard Institute for Space Studies2020b). LMDG is obtainable upon request from Martin Turbet ( and François Forget (

Author contributions
Back to toptop
Author contributions. 

TJF lead the THAI project and has written the article. ETW ran the simulation for Fig. 1 and plotted the figures. Every author contributed to the development of the THAI protocol, to the discussions and to the editing of the article.

Competing interests
Back to toptop
Competing interests. 

The author declares that there is no conflict of interest.

Back to toptop

Martin Turbet acknowledges the use of the computing resources of OCCIGEN (CINES, French national high-performance computing center).

The THAI GCM intercomparison team is grateful to Anong's Thai Cuisine restaurant in Laramie, WY, USA, for hosting its first meeting on 15 November 2017. We would like to thank the reviewer Daniel D. B. Koll, the anonymous reviewer and the topical editor Julia Hargreaves for comments that greatly improved our article.

Financial support
Back to toptop
Financial support. 

This project has received funding from the GSFC Sellers Exoplanet Environments Collaboration (SEEC), which is funded by the NASA Planetary Science Division's Internal Scientist Funding Model. This project has also received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie actions (grant agreement no. 832738/ESCAPE). Michael J. Way and Anthony D. Del Genio acknowledge funding from the NASA Astrobiology Program through participation in the Nexus for Exoplanet System Science (NExSS).

Review statement
Back to toptop
Review statement. 

This paper was edited by Julia Hargreaves and reviewed by Daniel D. B. Koll and one anonymous referee.

Back to toptop

Barstow, J. K. and Irwin, P. G. J.: Habitable worlds with JWST: transit spectroscopy of the TRAPPIST-1 system?, Mon. Not. R. Astron. Soc., 461, L92–L96,, 2016. a, b

Bolmont, E., Selsis, F., Owen, J. E., Ribas, I., Raymond, S. N., Leconte, J., and Gillon, M.: Water loss from terrestrial planets orbiting ultracool dwarfs: implications for the planets of TRAPPIST-1, Mon. Not. R. Astron. Soc., 464, 3728–3741,, 2017. a

Bourrier, V., de Wit, J., Bolmont, E., Stamenkovic, V., Wheatley, P. J., Burgasser, A. J., Delrez, L., Demory, B.-O., Ehrenreich, D., Gillon, M., Jehin, E., Leconte, J., Lederer, S. M., Lewis, N., Triaud, A. H. M. J., and Grootel, V. V.: Temporal Evolution of the High-energy Irradiation and Water Content of TRAPPIST-1 Exoplanets, Astron. J., 154, 121,, 2017. a

Boutle, I. A., Mayne, N. J., Drummond, B., Manners, J., Goyal, J., Hugo Lambert, F., Acreman, D. M., and Earnshaw, P. D.: Exploring the climate of Proxima B with the Met Office Unified Model, Astron. Astrophys., 601, A120,, 2017. a

Cantrell, J. R., Henry, T. J., and White, R. J.: The Solar Neighborhood XXIX: The Habitable Real Estate of Our Nearest Stellar Neighbors, Astron. J., 146, 99,, 2013. a

de Wit, J., Wakeford, H. R., Gillon, M., Lewis, N. K., Valenti, J. A., Demory, B.-O., Burgasser, A. J., Burdanov, A., Delrez, L., Jehin, E., Lederer, S. M., Queloz, D., Triaud, A. H. M. J., and Van Grootel, V.: A combined transmission spectrum of the Earth-sized exoplanets TRAPPIST-1 b and c, Nature, 537, 69–72,, 2016. a

de Wit, J., Wakeford, H. R., Lewis, N. K., Delrez, L., Gillon, M., Selsis, F., Leconte, J., Demory, B.-O., Bolmont, E., Bourrier, V., Burgasser, A. J., Grimm, S., Jehin, E., Lederer, S. M., Owen, J. E., Stamenković, V., and Triaud, A. H. M. J.: Atmospheric reconnaissance of the habitable-zone Earth-sized planets orbiting TRAPPIST-1, Nature Astronomy, 2, 214–219,, 2018. a

Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958,, 2016. a

Fauchez, T. J., Turbet, M., Villanueva, G. L., Wolf, E. T., Arney, G., Kopparapu, R. K., Mandell, A., de Wit, J., Pidhorodetska, D., and Domagal-Goldman, S. D.: Impact of Clouds and Hazes on the Simulated JWST Transmission Spectra of Habitable Zone Planets in the TRAPPIST-1 System, Astrophys. J., 887, 194,, 2019. a, b

Gillon, M., Jehin, E., Lederer, S. M., Delrez, L., de Wit, J., Burdanov, A., Van Grootel, V., Burgasser, A. J., Triaud, A. H. M. J., Opitom, C., Demory, B.-O., Sahu, D. K., Bardalez Gagliuffi, D., Magain, P., and Queloz, D.: Temperate Earth-sized planets transiting a nearby ultracool dwarf star, Nature, 533, 221–224,, 2016. a

Gillon, M., Triaud, A. H. M. J., Demory, B.-O., Jehin, E., Agol, E., Deck, K. M., Lederer, S. M., de Wit, J., Burdanov, A., Ingalls, J. G., Bolmont, E., Leconte, J., Raymond, S. N., Selsis, F., Turbet, M., Barkaoui, K., Burgasser, A., Burleigh, M. R., Carey, S. J., Chaushev, A., Copperwheat, C. M., Delrez, L., Fernandes, C. S., Holdsworth, D. L., Kotze, E. J., Van Grootel, V., Almleaky, Y., Benkhaldoun, Z., Magain, P., and Queloz, D.: Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1, Nature, 542, 456–460,, 2017. a, b, c

Goldblatt, C., Kavanagh, L., and Dewey, M.: The Palaeoclimate and Terrestrial Exoplanet Radiative Transfer Model Intercomparison Project (PALAEOTRIP): experimental design and protocols, Geosci. Model Dev., 10, 3931–3940,, 2017. a

Grimm, S. L., Demory, B.-O., Gillon, M., Dorn, C., Agol, E., Burdanov, A., Delrez, L., Sestovic, M., Triaud, A. H. M. J., Turbet, M., Bolmont, É., Caldas, A., Wit, J. D., Jehin, E., Leconte, J., Raymond, S. N., Grootel, V. V., Burgasser, A. J., Carey, S., Fabrycky, D., Heng, K., Hernandez, D. M., Ingalls, J. G., Lederer, S., Selsis, F., and Queloz, D.: The nature of the TRAPPIST-1 exoplanets, Astron. Astrophys., 613, A68,, 2018. a, b, c

Grootel, V. V., Fernandes, C. S., Gillon, M., Jehin, E., Manfroid, J., Scuflaire, R., Burgasser, A. J., Barkaoui, K., Benkhaldoun, Z., Burdanov, A., Delrez, L., Demory, B.-O., de Wit, J., Queloz, D., and Triaud, A. H. M. J.: Stellar Parameters for Trappist-1, Astrophys. J., 853, 30, doi:10.3847/1538-4357/aaa023, 2018. a

Haberle, R. M., Zahnle, K., and Barlow, N.: Impact Delivery of Reduced Greenhouse Gases on Early Mars, in: Fourth International Conference on Early Mars: Geologic, Hydrologic, and Climatic Evolution and the Implications for Life, vol. 2014 of LPI Contributions, p. 3022, 2017. a

Joshi, M. M. and Haberle, R. M.: Suppression of the Water Ice and Snow Albedo Feedback on Planets Orbiting Red Dwarf Stars and the Subsequent Widening of the Habitable Zone, Astrobiology, 12, 3–8,, 2012. a

Kite, E. S.: Geologic Constraints on Early Mars Climate, Space Sci. Rev., 215, 10,, 2019. a

Kopparapu, R. K., Ramirez, R., Kasting, J. F., Eymet, V., Robinson, T. D., Mahadevan, S., Terrien, R. C., Domagal-Goldman, S., Meadows, V., and Deshpande, R.: Habitable Zones around Main-sequence Stars: New Estimates, Astrophys. J., 765, 131,, 2013. a

Lacis, A. A. and Oinas, V.: A description of the correlated k distribution method for modeling nongray gaseous absorption, thermal emission, and multiple scattering in vertically inhomogeneous atmospheres, J. Geophys. Res.-Atmos., 96, 9027–9063,, 1991. a

Lincowski, A. P., Meadows, V. S., Crisp, D., Robinson, T. D., Luger, R., Lustig-Yaeger, J., and Arney, G. N.: Evolved Climates and Observational Discriminants for the TRAPPIST-1 Planetary System, Astrophys. J., 867, 76,, 2018. a

Luger, R., Sestovic, M., Kruse, E., Grimm, S. L., Demory, B.-O., Agol, E., Bolmont, E., Fabrycky, D., Fernandes, C. S., Van Grootel, V., Burgasser, A., Gillon, M., Ingalls, J. G., Jehin, E., Raymond, S. N., Selsis, F., Triaud, A. H. M. J., Barclay, T., Barentsen, G., Howell, S. B., Delrez, L., de Wit, J., Foreman-Mackey, D., Holdsworth, D. L., Leconte, J., Lederer, S., Turbet, M., Almleaky, Y., Benkhaldoun, Z., Magain, P., Morris, B. M., Heng, K., and Queloz, D.: A seven-planet resonant chain in TRAPPIST-1, Nature Astronomy, 1, 0129,, 2017. a

Mayne, N. J., Baraffe, I., Acreman, D. M., Smith, C., Wood, N., Amundsen, D. S., Thuburn, J., and Jackson, D. R.: Using the UM dynamical cores to reproduce idealised 3-D flows, Geosci. Model Dev., 7, 3059–3087,, 2014. a

Met Office: Unified Model, available at:, last access: 8 February 2020. a

Moran, S. E., Hörst, S. M., Batalha, N. E., Lewis, N. K., and Wakeford, H. R.: Limits on Clouds and Hazes for the TRAPPIST-1 Planets, Astron. J., 156, 252,, 2018. a

Morley, C. V., Kreidberg, L., Rustamkulov, Z., Robinson, T., and Fortney, J. J.: Observing the Atmospheres of Known Temperate Earth-sized Planets with JWST, Astrophys. J., 850, 121,, 2017. a, b, c

NASA Goddard Institute for Space Studies: GISS Model E, available at:, last access: 8 February 2020a. a

NASA Goddard Institute for Space Studies: ROCKE-3D workshop tutorials, available at:, last access: 8 February 2020b. a

Neale, R. B.: Description of the NCAR Community Atmosphere Model (CAM 4.0), Tech. rep., NCAR/TN-486+STR NCAR TECHNICAL NOTE, 2010. a

O'Malley-James, J. T. and Kaltenegger, L.: UV surface habitability of the TRAPPIST-1 system, Mon. Not. R. Astron. Soc., 469, L26–L30,, 2017. a

Shields, A. L., Meadows, V. S., Bitz, C. M., Pierrehumbert, R. T., Joshi, M. M., and Robinson, T. D.: The Effect of Host Star Spectral Energy Distribution and Ice-Albedo Feedback on the Climate of Extrasolar Planets, Astrobiology, 13, 715–739,, 2013. a

Turbet, M., Bolmont, E., Leconte, J., Forget, F., Selsis, F., Tobie, G., Caldas, A., Naar, J., and Gillon, M.: Modeling climate diversity, tidal dynamics and the fate of volatiles on TRAPPIST-1 planets, Astron. Astrophys., 612, A86,, 2018. a, b, c

Vida, K. and Roettenbacher, R. M.: Finding flares in Kepler data using machine-learning tools, Astron. Astrophys., 616, A163,, 2018. a

Villanueva, G. L., Smith, M. D., Protopapa, S., Faggi, S., and Mandell, A. M.: Planetary Spectrum Generator: An accurate online radiative transfer suite for atmospheres, comets, small bodies and exoplanets, J. Quant. Spectrosc. Ra., 217, 86–104,, 2018. a

von Paris, P., Selsis, F., Kitzmann, D., and Rauer, H.: The Dependence of the Ice-Albedo Feedback on Atmospheric Properties, Astrobiology, 13, 899–909,, 2013. a

Way, M. J., Aleinov, I., Amundsen, D. S., Chandler, M. A., Clune, T. L., Genio, A. D. D., Fujii, Y., Kelley, M., Kiang, N. Y., Sohl, L., and Tsigaridis, K.: Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets, Astrophys. J. Suppl. S., 231, 12 pp.,, 2017. a

Wheatley, P. J., Louden, T., Bourrier, V., Ehrenreich, D., and Gillon, M.: Strong XUV irradiation of the Earth-sized exoplanets orbiting the ultracool dwarf TRAPPIST-1, Mon. Not. R. Astron. Soc., 465, L74–L78,, 2017. a

Wolf, E. T.: Assessing the Habitability of the TRAPPIST-1 System Using a 3D Climate Model, Astrophys. J. Lett., 839, L1,, 2017. a, b

Wolf, E. T.: Erratum: “Assessing the Habitability of the TRAPPIST-1 System Using a 3D Climate Model” (2017, ApJL, 839, L1), Astrophys. J., 855, L14,, 2018. a

Wolf, E. T. and Toon, O. B.: The evolution of habitable climates under the brightening Sun, J. Geophys. Res.-Atmos., 120, 5775–5794,, 2015.  a, b

Wolf, E. T., Kopparapu, R., Airapetian, V., Fauchez, T., Guzewich, S. D., Kane, S. R., Pidhorodetska, D., Way, M. J., Abbot, D. S., Checlair, J. H., Davis, C. E., Del Genio, A., Dong, C., Eggl, S., Fleming, D. P., Fujii, Y., Haghighipour, N., Heavens, N., Henning, W. G., Kiang, N. Y., Lopez-Morales, M., Lustig-Yaeger, J., Meadows, V., Reinhard, C. T., Rugheimer, S., Schwieterman, E. W., Shields, A. L., Sohl, L., Turbet, M., and Wordsworth, R. D.: The Importance of 3D General Circulation Models for Characterizing the Climate and Habitability of Terrestrial Extrasolar Planets, NAS ASTRO 2020, White Paper, arXiv e-prints, arXiv:1903.05012, 2019. a

Wordsworth, R. D., Forget, F., Selsis, F., Millour, E., Charnay, B., and Madeleine, J.-B.: Gliese 581d is the First Discovered Terrestrial-mass Exoplanet in the Habitable Zone, Astrophys. J. Lett., 733, L48,, 2011. a

Yang, J., Boué, G., Fabrycky, D. C., and Abbot, D. S.: Strong Dependence of the Inner Edge of the Habitable Zone on Planetary Rotation Rate, Astrophys. J., 787, L2,, 2014. a

Yang, J., Leconte, J., Wolf, E. T., Merlis, T., Koll, D. D. B., Forget, F., and Abbot, D. S.: Simulations of Water Vapor and Clouds on Rapidly Rotating and Tidally Locked Planets: A 3D Model Intercomparison, Astrophys. J., 875, 46,, 2019. a, b, c

Publications Copernicus
Short summary
Atmospheric characterization of rocky exoplanets orbiting within the habitable zone of nearby M dwarf stars is around the corner with the James Webb Space Telescope (JWST), expected to be launch in 2021. Global climate models (GCMs) are powerful tools to model exoplanet atmospheres and to predict their habitability. However, intrinsic differences between the models can lead to various predictions. This paper presents an experiment protocol to evaluate these differences.
Atmospheric characterization of rocky exoplanets orbiting within the habitable zone of nearby M...