Stable isotope analysis and temperature reconstruction data from DSDP Hole 94-609B and ODP Hole 162-984B

The Greenland ice sheet is accepted as a key factor controlling the Quaternary glacial scenario. However, the origin and mechanisms of major Arctic glaciation starting at 3.15 Ma and culminating at 2.74 Ma are still controversial. For this phase of intense cooling Ravelo et al. proposed a complex gradual forcing mechanism. In contrast, our new submillennial-scale paleoceanographic records from the Pliocene North Atlantic suggest a far more precise timing and forcing for the initiation of northern hemisphere glaciation (NHG), since it was linked to a 2-3 °C surface water warming during warm stages from 2.95 to 2.82 Ma. These records support previous models, claiming that the final closure of the Panama Isthmus (3.0- ~2.5 Ma induced an increased poleward salt and heat transport. Associated strengthening of North Atlantic Thermohaline Circulation and in turn, an intensified moisture supply to northern high latitudes resulted in the build-up of NHG, finally culminating in the great, irreversible climate crash at marine isotope stage G6 (2.74 Ma). In summary, there was a two-step threshold mechanism that marked the onset of NHG with glacial-to-interglacial cycles quasi-persistent until today.

Neogene calcareous nannofossil abundance and datums of ODP Site 186-1150 and 186-1151 sediments

Calcareous nannofossil assemblages were studied from Ocean Drilling Program Holes 1150A, 1150B, 1151A, 1151C, and 1151D in order to estimate the age of sediments drilled in the Japan Trench of the western Pacific Ocean. The abundance and species diversity of nannofossil flora are generally low but are sufficient to show that the sedimentary sequences range from Quaternary to Miocene in age (nannofossil Zones CN15-CN3). The abundance of Coccolithus pelagicus, a cold-water indicator, was studied from sediments younger than 3.83 Ma from both Holes 1150A and 1151A in order to elucidate past climate conditions. Between 3.83 and 2.82 Ma, the abundance of C. pelagicus was generally low, but abundance increased significantly after 2.82 Ma. In agreement with previous studies, this increase appears to be related to a change in the current system around the western Pacific Ocean and eastern Atlantic Ocean that occurred in response to the final elevation of the Isthmus of Panama.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

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.

Stable isotope analysis of benthic foraminifera and sedimentology on ODP Site 162-984

Sediments recovered at ODP Site 984 on the Reykjanes Ridge provided multicentennial-scale records (SST, planktic and benthic delta18O, IRD and magnetic susceptibility) of Late Pliocene climate change over the onset of Northern Hemisphere glaciation (NHG), 2.95-2.82 Ma. Short-term climate variations prior and after the onset of continent-wide glaciation were compared to test the hypothesis, whether Dansgaard-Oeschger (DO) cycles may have been triggered by continental ice breakouts. By means of spectral analyses for two selected interglacial stages prior to and after NHG (G15 and G1), we found that climate variability resembled that of the Holocene and the mid-Pliocene warm period. In contrast, DO-like periodicities near 1470, 2900, and 4400 yr indeed only occurred in glacial stages after the onset of NHG (G14, G6, and 104), but hardly in stage G20 prior to the onset. These results suggest a causal link between DO cycles and the Late Pliocene onset of major NHG and ice breakouts in the North Atlantic. This data set provides all primary data and spectral analysis related to this scientific work.