Abstract. We present results from an ensemble of seven 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 DeepMIP (www.deepmip.org), and as such all models have been configured with identical paleogeographic and vegetation boundary conditions. The results indicate that these non-CO2 boundary conditions contribute between 3 and 5 °C to Eocene warmth. Compared to results from previous studies, the DeepMIP simulations show reduced spread of global mean surface temperature response across the ensemble, for a given atmospheric CO2 concentration. In a marked departure from the results from previous simulations, at least two of the DeepMIP models (CESM and GFDL) are consistent with proxy indicators of global mean temperature, and atmospheric CO2, and meridional SST gradients. The best agreement with global SST proxies from these models occurs at CO2 concentrations of around 2400 ppmv. At a more regional scale the models lack skill in reproducing the proxy SSTs, in particular in the southwest Pacific, around New Zealand and south Australia, where the modelled anomalies are substantially less than indicated by the proxies. However, in these regions modelled continental surface air temperature anomalies are consistent with surface air temperature proxies, implying an 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 and aerosols.