Description: 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 (MIS), MIS KM5c and K1.

Description: Model output of the intermediate complexity climate model CLIMBER-2 over the past 5 million years. The simulations were forced with insolation data (O), insolation and land ice data (OI), insolation and carbon dioxide data (OC) and with insolation, land ice and carbon dioxide data (OIC). Sheet 1 contains the main results: northern hemispheric (30-90 deg N), southern hemispheric (30-90 deg S) and global temperatures. Sheet 2 contains the land ice and carbon dioxide forcing in terms of globally averaged radiative forcing. Details are given in the publication. More information or data can be obtained by contacting L.B. Stap (lennert.stap@awi.de).

Description: It is still an open question how equilibrium warming in response to increasing radiative forcing – the specific equilibrium climate sensitivity S – depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6–5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Description: Marine sediment records from the Oligocene and Miocene reveal clear 400,000-year (400-kyr) climate cycles related to variations in orbital eccentricity. These cycles are also observed in the Plio-Pleistocene records of the global carbon cycle. However they are absent in the Late Pleistocene ice-age record over the past 1.5 million years. Here, we present a simulation of global ice volume over the past 5 million years with a coupled system of four 3-D ice-sheet models. Our simulation shows that the 400-kyr long eccentricity cycles of Antarctica vary coherently with d13C records during the Pleistocene suggesting that they drive the long-term carbon cycle changes throughout the past 35 million years. The 400-kyr response of Antarctica is eventually suppressed by the dominant 100-kyr glacial cycles of the large ice sheets in the Northern Hemisphere (NH).

Description: Model output of a coupled ice sheet-climate model, inversely forced by benthic d18O over the past 38 million years. Sheet 1 contains the main results from the reference simulation: benthic d18O, CO2, ice-volume-equivalent sea level and global temperature. Sheet 2 contains global, Northern Hemisphere (40-80 deg N), and Antarctic (60-90 deg S) temperatures, from the reference run and the run with ice uncoupled, only albedo coupled, and only surface height coupled. Sheet 3 contains global temperature, from the reference run, and the runs with fixed PD ice, fixed LGM ice, and no ice. Details are given in the publication. More information or data can be obtained by contacting L.B. Stap (lennert.stap@awi.de).

Description: During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.