Much of our understanding of Earth’s past climate comes from the measurement of oxygen and carbon isotope variations in deep-sea benthic foraminifera. Yet, long intervals in existing records lack the temporal resolution and age control needed to thoroughly categorize climate states of the Cenozoic era and to study their dynamics. Here, we present a new, highly resolved, astronomically dated, continuous composite of benthic foraminifer isotope records developed in our laboratories. Four climate states—Hothouse, Warmhouse, Coolhouse, Icehouse—are identified on the basis of their distinctive response to astronomical forcing depending on greenhouse gas concentrations and polar ice sheet volume. Statistical analysis of the nonlinear behavior encoded in our record reveals the key role that polar ice volume plays in the predictability of Cenozoic climate dynamics.

Global warming during the Palaeocene-Eocene Thermal Maximum1,2 (PETM, ~56 Ma) is commonly interpreted as being primarily driven by the destabilization of carbon from surficial sedimentary reservoirs such as methane hydrates3. However, the source(s) of carbon remain controversial1,3–5. Resolving this is key to understanding the proximal cause, as well as quantifying the roles of triggers versus feedbacks in driving the event. Here we present new boron isotope data – a proxy for seawater pH – that demonstrate the occurrence of persistently suppressed surface ocean pH across the PETM. Our pH data, alongside a paired carbon isotope record, are assimilated in an Earth system model to reconstruct the unfolding carbon cycle dynamics across the event6,7. We find strong evidence for a much larger (>10,000 PgC) and on average isotopically heavier carbon source than considered previously8,9. This leads us to identify volcanism associated with the North Atlantic Igneous Province, rather than carbon from a surficial reservoir, as the main driver of the PETM10,11. We also find that, although amplifying organic carbon feedbacks with climate likely played only a subordinate role in driving the event, enhanced organic matter burial was important in ultimately sequestering the released carbon and accelerating the recovery of the Earth system12.

The Early Eocene Climate Optimum (EECO, which occurred about 51 to 53 million years ago) 1 , was the warmest interval of the past 65 million years, with mean annual surface air temperature over ten degrees Celsius warmer than during the pre-industrial period 2-4 . Subsequent global cooling in the middle and late Eocene epoch, especially at high latitudes, eventually led to continental ice sheet development in Antarctica in the early Oligocene epoch (about 33.6 million years ago). However, existing estimates place atmospheric carbon dioxide (CO 2 ) levels during the Eocene at 500-3,000 parts per million [5][6][7] , and in the absence of tighter constraints carbon-climate interactions over this interval remain uncertain. Here we use recent analytical and methodological developments 8-11 to generate a new high-fidelity record of CO 2 concentrations using the boron isotope (δ 11 Β) composition of well preserved planktonic foraminifera from the Tanzania Drilling Project, revising previous estimates 6 . Although species-level uncertainties make absolute values difficult to constrain, CO 2 concentrations during the EECO were around 1,400 parts per million. The relative decline in DOI: 10.1038/nature17423 Page 2 of 38 CO 2 concentration through the Eocene is more robustly constrained at about fifty per cent, with a further decline into the Oligocene 12 . Provided the latitudinal dependency of sea surface temperature change for a given climate forcing in the Eocene was similar to that of the late Quaternary period 13 , this CO 2 decline was sufficient to drive the well documented high-and low-latitude cooling that occurred through the Eocene 14 . Once the change in global temperature between the pre-industrial period and the Eocene caused by the action of all known slow feedbacks (apart from those associated with the carbon cycle) is removed 2-4 , both the EECO and the late Eocene exhibit an equilibrium climate sensitivity relative to the pre-industrial period of 2.1 to 4.6 degrees Celsius per CO 2 doubling (66 per cent confidence), which is similar to the canonical range (1.5 to 4.5 degrees Celsius 15 ), indicating that a large fraction of the warmth of the early Eocene greenhouse was driven by increased CO 2 concentrations, and that climate sensitivity was relatively constant throughout this period.Over the past 540 million years, Earth's climate has oscillated between a globally warm 'greenhouse state' and an 'icehouse state' with substantial continental glaciation 16 . The most recent of these transitions occurred between the warmest time interval of the last 65 million years-the EECO (about 14 ± 3 °C warmer than preindustrial times 2 )-and the rapid growth of ice on Antarctica in the earliest icehouse state of the Oligocene (~33.6 Myr ago 1 ). It has been suggested that variations in the concentration of the greenhouse gas CO 2 were responsible for both the overall warmth of the Eocene and the subsequent cooling 17 . Recent studies have documented the importance of CO 2 decline for the final step into the icehous...

Abstract. The early Eocene (56 to 48 million years ago) is inferred to have been the most recent time that Earth's atmospheric CO2 concentrations exceeded 1000 ppm. Global mean temperatures were also substantially warmer than those of the present day. As such, the study of early Eocene climate provides insight into how a super-warm Earth system behaves and offers an opportunity to evaluate climate models under conditions of high greenhouse gas forcing. The Deep Time Model Intercomparison Project (DeepMIP) is a systematic model–model and model–data intercomparison of three early Paleogene time slices: latest Paleocene, Paleocene–Eocene thermal maximum (PETM) and early Eocene climatic optimum (EECO). A previous article outlined the model experimental design for climate model simulations. In this article, we outline the methodologies to be used for the compilation and analysis of climate proxy data, primarily proxies for temperature and CO2. This paper establishes the protocols for a concerted and coordinated effort to compile the climate proxy records across a wide geographic range. The resulting climate “atlas” will be used to constrain and evaluate climate models for the three selected time intervals and provide insights into the mechanisms that control these warm climate states. We provide version 0.1 of this database, in anticipation that this will be expanded in subsequent publications.

Abstract. Past warm periods provide an opportunity to evaluate climate models under extreme forcing scenarios, in particular high ( >  800 ppmv) atmospheric CO2 concentrations. Although a post hoc intercomparison of Eocene ( ∼  50  Ma) climate model simulations and geological data has been carried out previously, models of past high-CO2 periods have never been evaluated in a consistent framework. Here, we present an experimental design for climate model simulations of three warm periods within the early Eocene and the latest Paleocene (the EECO, PETM, and pre-PETM). Together with the CMIP6 pre-industrial control and abrupt 4 ×  CO2 simulations, and additional sensitivity studies, these form the first phase of DeepMIP – the Deep-time Model Intercomparison Project, itself a group within the wider Paleoclimate Modelling Intercomparison Project (PMIP). The experimental design specifies and provides guidance on boundary conditions associated with palaeogeography, greenhouse gases, astronomical configuration, solar constant, land surface processes, and aerosols. Initial conditions, simulation length, and output variables are also specified. Finally, we explain how the geological data sets, which will be used to evaluate the simulations, will be developed.

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.

No abstract

Despite recent advances, the link between the evolution of atmospheric CO2 and climate during the Eocene greenhouse remains uncertain. In particular, modelling studies suggest that in order to achieve the global warmth that characterised the early Eocene, warmer climates must be more sensitive to CO2 forcing than colder climates. Here, we test this assertion in the geological record by combining a new high-resolution boron isotope-based CO2 record with novel estimates of Global Mean Temperature. We find that Equilibrium Climate Sensitivity (ECS) was indeed higher during the warmest intervals of the Eocene, agreeing well with recent model simulations, and declined through the Eocene as global climate cooled. These observations indicate that the canonical IPCC range of ECS (1.5 to 4.5 °C per doubling) is unlikely to be appropriate for high-CO2 warm climates of the past, and the state dependency of ECS may play an increasingly important role in determining the state of future climate as the Earth continues to warm.