DeepMIP: model intercomparison of early Eocene climatic optimum (EECO) large-scale climate features and comparison with proxy data
Authors: Lunt Daniel J. , Bragg Fran , Chan Wing-Le , Hutchinson David K. , Ladant Jean-Baptiste , Morozova Polina , Niezgodzki Igor , Steinig Sebastian , Zhang Zhongshi , Zhu Jiang , Abe-Ouchi Ayako , Anagnostou Eleni , de Boer Agatha M. , Coxall Helen K. , Donnadieu Yannick , Foster Gavin , Inglis Gordon N. , Knorr Gregor , Langebroek Petra M. , Lear Caroline H. , Lohmann Gerrit , Poulsen Christopher J. , Sepulchre Pierre , Tierney Jessica E. , Valdes Paul J. , Volodin Evgeny M. , Dunkley Jones Tom , Hollis Christopher J. , Huber Matthew , Otto-Bliesner Bette L. 
Published: January 2021

Abstract. We present results from an ensemble of eight climate models, each
of which has carried out simulations
of the early Eocene climate optimum (EECO, ∼ 50 million years
ago). These simulations have been carried out in the framework of the Deep-Time Model Intercomparison Project
(DeepMIP; http://www.deepmip.org, last access: 10 January 2021); thus, all models have been configured with the same
paleogeographic and vegetation boundary conditions. The results indicate that
these non-CO2 boundary conditions contribute between 3 and
5 ∘C to Eocene warmth. Compared with
results from previous studies, the DeepMIP simulations generally show a reduced spread of the global mean surface temperature response across the ensemble for a given atmospheric CO2 concentration as well as an increased climate sensitivity on average. An energy balance analysis of the model ensemble indicates that global mean warming in the Eocene compared with the preindustrial period mostly arises from decreases in emissivity due to the elevated CO2 concentration (and associated water vapour and long-wave cloud feedbacks), whereas the reduction in the Eocene in terms of the meridional temperature gradient is primarily due to emissivity and albedo changes owing to the non-CO2 boundary conditions (i.e. the removal of the Antarctic ice sheet and changes in vegetation). Three of the models (the Community Earth System Model, CESM; the Geophysical Fluid Dynamics Laboratory, GFDL, model; and the Norwegian Earth System Model, NorESM) show results that are consistent with the proxies in terms of the global mean temperature, meridional SST gradient, and CO2, without prescribing changes to model parameters. In addition, many of the models agree well with the first-order spatial patterns in the SST proxies. However, at a more regional scale, the models lack skill. In particular, the modelled anomalies are substantially lower than those indicated by the proxies in the southwest Pacific; here, modelled continental surface air temperature anomalies are more consistent with surface air temperature proxies, implying a possible inconsistency between marine and terrestrial temperatures in either the proxies or models in this region. Our aim is that the documentation of the large-scale features and model–data comparison presented herein will pave the way to further studies that explore aspects of the model simulations in more detail, for example the ocean
circulation, hydrological cycle, and modes of variability, and encourage sensitivity studies to aspects such as paleogeography, orbital configuration, and aerosols.