As the world warms due to rising greenhouse gas concentrations, the Earth system moves toward climate states without societal precedent, challenging adaptation. Past Earth system states offer possible model systems for the warming world of the coming decades. These include the climate states of the Early Eocene (ca. 50 Ma), the Mid-Pliocene (3.3–3.0 Ma), the Last Interglacial (129–116 ka), the Mid-Holocene (6 ka), preindustrial (ca. 1850 CE), and the 20th century. Here, we quantitatively assess the similarity of future projected climate states to these six geohistorical benchmarks using simulations from the Hadley Centre Coupled Model Version 3 (HadCM3), the Goddard Institute for Space Studies Model E2-R (GISS), and the Community Climate System Model, Versions 3 and 4 (CCSM) Earth system models. Under the Representative Concentration Pathway 8.5 (RCP8.5) emission scenario, by 2030 CE, future climates most closely resemble Mid-Pliocene climates, and by 2150 CE, they most closely resemble Eocene climates. Under RCP4.5, climate stabilizes at Pliocene-like conditions by 2040 CE. Pliocene-like and Eocene-like climates emerge first in continental interiors and then expand outward. Geologically novel climates are uncommon in RCP4.5 (<1%) but reach 8.7% of the globe under RCP8.5, characterized by high temperatures and precipitation. Hence, RCP4.5 is roughly equivalent to stabilizing at Pliocene-like climates, while unmitigated emission trajectories, such as RCP8.5, are similar to reversing millions of years of long-term cooling on the scale of a few human generations. Both the emergence of geologically novel climates and the rapid reversion to Eocene-like climates may be outside the range of evolutionary adaptive capacity.

Given the inherent uncertainties in predicting how climate and environments will respond to anthropogenic emissions of greenhouse gases, it would be beneficial to society if science could identify geological analogues to the human race’s current grand climate experiment . This has been a focus of the geological and palaeoclimate communities over the last 30 years, with many scientific papers claiming that intervals in Earth history can be used as an analogue for future climate change. Using a coupled ocean–atmosphere modelling approach, we test this assertion for the most probable pre-Quaternary candidates of the last 100 million years: the Mid- and Late Cretaceous, the Palaeocene–Eocene Thermal Maximum (PETM), the Early Eocene, as well as warm intervals within the Miocene and Pliocene epochs. These intervals fail as true direct analogues since they either represent equilibrium climate states to a long-term CO ₂ forcing—whereas anthropogenic emissions of greenhouse gases provide a progressive (transient) forcing on climate—or the sensitivity of the climate system itself to CO ₂ was different. While no close geological analogue exists, past warm intervals in Earth history provide a unique opportunity to investigate processes that operated during warm (high CO ₂ ) climate states. Palaeoclimate and environmental reconstruction/modelling are facilitating the assessment and calculation of the response of global temperatures to increasing CO ₂ concentrations in the longer term (multiple centuries); this is now referred to as the Earth System Sensitivity, which is critical in identifying CO ₂ thresholds in the atmosphere that must not be crossed to avoid dangerous levels of climate change in the long term. Palaeoclimatology also provides a unique and independent way to evaluate the qualities of climate and Earth system models used to predict future climate.

The characteristics of the mid-Pliocene warm period (mPWP: 3.264–3.025 Ma BP) have been examined using geological proxies and climate models. While there is agreement between models and data, details of regional climate differ. Uncertainties in prescribed forcings and in proxy data limit the utility of the interval to understand the dynamics of a warmer than present climate or evaluate models. This uncertainty comes, in part, from the reconstruction of atime slabrather than atime slice, where forcings required by climate models can be more adequately constrained. Here, we describe the rationale and approach for identifying a time slice(s) for Pliocene environmental reconstruction. A time slice centred on 3.205 Ma BP (3.204–3.207 Ma BP) has been identified as a priority for investigation. It is a warm interval characterized by a negative benthic oxygen isotope excursion (0.21–0.23‰) centred on marine isotope stage KM5c (KM5.3). It occurred during a period of orbital forcing that was very similar to present day. Climate model simulations indicate that proxy temperature estimates are unlikely to be significantly affected by orbital forcing for at least a precession cycle centred on the time slice, with the North Atlantic potentially being an important exception.

In modern environmental and climate science it is necessary to assimilate observational datasets collected over decades with outputs from numerical models, to enable a full understanding of natural systems and their sensitivities. During the twentieth and twenty-first centuries, numerical modelling became central to many areas of science from the Bohr model of the atom to the Lorenz model of the atmosphere. In modern science, a great deal of time and effort is devoted to developing, evaluating, comparing and modifying numerical models that help us synthesise our understanding of complex natural systems. Here we provide an assessment of the contribution of past (palaeo) climate modelling to multidisciplinary science and to society by answering the following question: What can palaeoclimate modelling do for you? We provide an assessment of how palaeoclimate modelling can develop in the future to further enhance multidisciplinary research that aims to understand Earth’s evolution, and what this may tell us about the resilience of natural and social systems as we enter the Anthropocene.

Accurate understanding of permafrost dynamics is critical for evaluating and mitigating impacts that may arise as permafrost degrades in the future; however, existing projections have large uncertainties. Studies of how permafrost responded historically during Earth’s past warm periods are helpful in exploring potential future permafrost behavior and to evaluate the uncertainty of future permafrost change projections. Here, we combine a surface frost index model with outputs from the second phase of the Pliocene Model Intercomparison Project to simulate the near‐surface (~3 to 4 m depth) permafrost state in the Northern Hemisphere during the mid-Pliocene warm period (mPWP, ~3.264 to 3.025 Ma). This period shares similarities with the projected future climate. Constrained by proxy-based surface air temperature records, our simulations demonstrate that near‐surface permafrost was highly spatially restricted during the mPWP and was 93 ± 3% smaller than the preindustrial extent. Near‐surface permafrost was present only in the eastern Siberian uplands, Canadian high Arctic Archipelago, and northernmost Greenland. The simulations are similar to near‐surface permafrost changes projected for the end of this century under the SSP5-8.5 scenario and provide a perspective on the potential permafrost behavior that may be expected in a warmer world.

Examining the nature of ice sheet and sea level response to past episodes of enhanced greenhouse gas forcing may help constrain future sea level change. Here, for the first time, we present the transient nature of ice sheets and sea level during the late Pliocene. The transient ice sheet predictions are forced by multiple climate snapshots derived from a climate model set up with late Pliocene boundary conditions, forced with different orbital forcing scenarios appropriate to two Marine Isotope Stages (MISs), MIS KM5c, and K1. Our results indicate that during MIS KM5c both the Antarctic and Greenland ice sheets contributed to sea level rise relative to present and were relatively stable. Insolation forcing between the hemispheres was out of phase during MIS K1 and led to an asynchronous response of ice volume globally. Therefore, when variations of precession were high, inferring the behavior of ice sheets from benthic isotope or sea level records is complex.