(Table S1) Calcium isotope ratios of Pacific Ocean sediments

Multiple lines of evidence have shown that the isotopic composition and concentration of calcium in seawater have changed over the past 28 million years. A high-resolution, continuous seawater calcium isotope ratio curve from marine (pelagic) barite reveals distinct features in the evolution of the seawater calcium isotopic ratio suggesting changes in seawater calcium concentrations. The most pronounced increase in the d44/40Ca value of seawater (of 0.3 per mil) occurred over roughly 4 million years following a period of low values around 13 million years ago. The major change in marine calcium corresponds to a climatic transition and global change in the carbon cycle and suggests a reorganization of the global biogeochemical system.

(Table S1) Neodymium isotopic results from ODP Site 172-1063

The Atlantic Meridional Overturning Circulation (AMOC) plays an important role in the Northern Hemisphere climate system. Significant interest went into the question of how excessive freshwater input through melting of continental ice can affect its overturning vigor and, hence, heat supply, to higher northern latitudes. Such forcing can be tested by investigating its behavior during extreme iceberg discharge events into the open North Atlantic during the last glacial period, the so-called Heinrich events (HE). Here we present neodymium (Nd) isotope compositions of past seawater, a sensitive chemical water mass tag, extracted from sediments of Ocean Drilling Program Site 1063 in the western North Atlantic (Bermuda Rise), covering the period surrounding HE 2, the Last Glacial Maximum, and the early deglaciation. These data are compared with a record of the kinematic circulation tracer (231Pa/230Th)xs extracted from the same sediment core. Both tracers indicate significant circulation changes preceding intense ice rafting during HE 2 by almost 2 kyr. Moreover, the Nd isotope record suggests the presence of deeply ventilating North Atlantic Deep Water early during Marine Isotope Stage 2 until it was replaced by Southern Source Water at ~27 ka. The early switch to high (Pa/Th)xs and radiogenic epsilon-Nd in relation to intensified ice rafting during HE 2 suggests that ice rafting into the open North Atlantic during major HE 2 was preceded by an early change of the AMOC. This opens the possibility that variations in AMOC contributed to or even triggered the ice sheet instability rather than merely responding to it.

(Table S1) Uranium and Thorium isotopes and U/Th ages of IODP Exp310 cores from offshore Tahiti

The timing of sea-level change provides important constraints on the mechanisms driving Earth's climate between glacial and interglacial states. Fossil corals constrain the timing of past sea level by their suitability for dating and their growth position close to sea level. The coral-derived age for the last deglaciation is consistent with climate change forced by Northern Hemisphere summer insolation (NHI), but the timing of the penultimate deglaciation is more controversial. We found, by means of uranium/thorium dating of fossil corals, that sea level during the penultimate deglaciation had risen to ~85 meters below the present sea level by 137,000 years ago, and that it fluctuated on a millennial time scale during deglaciation. This indicates that the penultimate deglaciation occurred earlier with respect to NHI than the last deglacial, beginning when NHI was at a minimum.

(Table S1) Age determination of Atlantic Ocean sediment cores

Well-dated benthic foraminifer oxygen isotopic records (d18O) from different water depths and locations within the Atlantic Ocean exhibit distinct patterns and significant differences in timing over the last deglaciation. This has two implications: on the one hand, it confirms that benthic d18O cannot be used as a global correlation tool with millennial-scale precision, but on the other hand, the combination of benthic isotopic records with independent dating provides a wealth of information on past circulation changes. Comparing new South Atlantic benthic isotopic data with published benthic isotopic records, we show that (1) circulation changes first affected benthic d18O in the 1000-2200 m range, with marked decreases in benthic d18O taking place at ~17.5 cal. kyr B.P. (ka) due to the southward propagation of brine waters generated in the Nordic Seas during Heinrich Stadial 1 (HS1) cold period; (2) the arrival of d18O-depleted deglacial meltwater took place later at deeper North Atlantic sites; (3) hydrographic changes recorded in North Atlantic cores below 3000 m during HS1 do not correspond to simple alternations between northern- and southern-sourced water but likely reflect instead the incursion of brine-generated deep water of northern as well as southern origin; and (4) South Atlantic waters at ~44°S and ~3800 m depth remained isolated from better-ventilated northern-sourced water masses until after the resumption of North Atlantic Deep Water (NADW) formation at the onset of the Bølling-Allerod, which led to the propagation of NADW into the South Atlantic.

(Table S1) Age determination of sediment core OCE326-GGC5

The Atlantic meridional overturning circulation is widely believed to affect climate. Changes in ocean circulation have been inferred from records of the deep water chemical composition derived from sedimentary nutrient proxies (Boyle and Keigwin, 1987, doi:10.1038/330035a0), but their impact on climate is difficult to assess because such reconstructions provide insufficient constraints on the rate of overturning (LeGrand and Wunsch, 1995, doi:10.1029/95PA01455). Here we report measurements of 231Pa/230Th, a kinematic proxy for the meridional overturning circulation, in a sediment core from the subtropical North Atlantic Ocean. We find that the meridional overturning was nearly, or completely, eliminated during the coldest deglacial interval in the North Atlantic region, beginning with the catastrophic iceberg discharge Heinrich event H1, 17,500 yr ago, and declined sharply but briefly into the Younger Dryas cold event, about 12,700 yr ago. Following these cold events, the 231Pa/230Th record indicates that rapid accelerations of the meridional overturning circulation were concurrent with the two strongest regional warming events during deglaciation. These results confirm the significance of variations in the rate of the Atlantic meridional overturning circulation for abrupt climate changes.

(Table S1) Chronostratigraphic constrains for sediment core AND1-1B

Thirty years after oxygen isotope records from microfossils deposited in ocean sediments confirmed the hypothesis that variations in the Earth's orbital geometry control the ice ages (Hays et al., 1976, doi:10.1126/science.194.4270.1121), fundamental questions remain over the response of the Antarctic ice sheets to orbital cycles (Raymo and Huybers, 2008, doi:10.1038/nature06589). Furthermore, an understanding of the behaviour of the marine-based West Antarctic ice sheet (WAIS) during the 'warmer-than-present' early-Pliocene epoch (~5-3 Myr ago) is needed to better constrain the possible range of ice-sheet behaviour in the context of future global warming (Solomon et al., 2007). Here we present a marine glacial record from the upper 600 m of the AND-1B sediment core recovered from beneath the northwest part of the Ross ice shelf by the ANDRILL programme and demonstrate well-dated, ~40-kyr cyclic variations in ice-sheet extent linked to cycles in insolation influenced by changes in the Earth's axial tilt (obliquity) during the Pliocene. Our data provide direct evidence for orbitally induced oscillations in the WAIS, which periodically collapsed, resulting in a switch from grounded ice, or ice shelves, to open waters in the Ross embayment when planetary temperatures were up to ~3° C warmer than today ( Kim and Crowley, 2000, doi:10.1029/1999PA000459) and atmospheric CO2 concentration was as high as ~400 p.p.m.v. (van der Burgh et al., 1993, doi:10.1126/science.260.5115.1788, Raymo et al., 1996, doi:10.1016/0377-8398(95)00048-8). The evidence is consistent with a new ice-sheet/ice-shelf model (Pollard and DeConto, 2009, doi:10.1038/nature07809) that simulates fluctuations in Antarctic ice volume of up to +7 m in equivalent sea level associated with the loss of the WAIS and up to +3 m in equivalent sea level from the East Antarctic ice sheet, in response to ocean-induced melting paced by obliquity. During interglacial times, diatomaceous sediments indicate high surface-water productivity, minimal summer sea ice and air temperatures above freezing, suggesting an additional influence of surface melt (Huybers, 2006, doi:10.1126/science.1125249) under conditions of elevated CO2.

(Table S1) Opal content of sediment core OXMZ01MV-GC31

The North American monsoon (NAM), an onshore wind shift occurring between July and September, has evolved in character during the Holocene largely due to changes in Northern Hemisphere insolation. Published paleoproxy and modeling studies suggest that prior to ~8000 cal years BP, the NAM affected a broader region than today, extending westward into the Mojave Desert of California. Holocene proxy SST records from the Gulf of California (GoC) and the adjacent Pacific provide constraints for this changing NAM climatology. Prior to ~8000 cal years BP, lower GoC SSTs would not have fueled northward surges of tropical moisture up the GoC, which presently contribute most of the monsoon precipitation to the western NAM region. During the early Holocene, the North Pacific High was further north and SSTs in the California Current off Baja California were warmer, allowing monsoonal moisture flow from the subtropical Pacific to take a more direct, northwesterly trajectory into an expanded area of the southwestern U.S. west of 114°W. A new upwelling record off southwest Baja California reveals that enhanced upwelling in the California Current beginning at ~7500 cal year BP may have triggered a change in NAM climatology, focusing the geographic expression of NAM in the southwest USA into its modern core region east of ~114°W, in Arizona and New Mexico. Holocene proxy precipitation records from the southwestern U.S. and northwestern Mexico, including lakes, vegetation/pollen, and caves are reviewed and found to be largely supportive of this hypothesis of changing Holocene NAM climatology.

(Table S1) Age determination of sediment cores SK99_17 and MD76-131

Productivity in the Arabian Sea is one of the highest in the world. It is controlled by seasonally reversing monsoonal wind-driven upwelling of nutrient-rich deeper waters which fuel phytoplankton growth. The detailed history of upwelling-induced productivity in the eastern Arabian Sea is unknown. Here we present paleoproductivity records from a composite sediment core at the millennial scale during the last 80 kyr B.P. These records are based on relative abundance counts of planktonic foraminifera and organic carbon contents, which are shown to mainly vary in concert. The eastern Arabian Sea upwelling-induced productivity was higher in the glacial period than in the Holocene, but it fell repeatedly on millennial timescales. These productivity declines occurred during cold events in the North Atlantic region, with the most pronounced changes prevailing during the Heinrich events. Hence, seasonal monsoon winds that drive upwelling-induced productivity in the east were weak when the North Atlantic was cold. These weak winds resulted in stratification of the water column, comparable to today's Arabian Sea stratification in the intermonsoonal period. Combining the new eastern with published western Arabian Sea results shows that the entire biological factory was severely diminished during the North Atlantic Heinrich events, and the seasonal productivity change in the Arabian Sea monsoon system was reduced with year-round low productivity.

(Table S1) Diatom abundance according to ice cycles in ANDRILL AND1-1B drill core

The influence of Antarctica and the Southern Ocean on Late Pliocene global climate reconstructions has remained ambiguous due to a lack of well-dated Antarctic-proximal, paleoenvironmental records. Here we present ice sheet, sea-surface temperature, and sea ice reconstructions from the ANDRILL AND-1B sediment core recovered from beneath the Ross Ice Shelf. We provide evidence for a major expansion of an ice sheet in the Ross Sea that began at ~3.3 Ma, followed by a coastal sea surface temperature cooling of ~2.5°C, a stepwise expansion of sea ice, and polynya-style deep mixing in the Ross Sea between 3.3 and 2.5 Ma. The intensification of Antarctic cooling resulted in strengthened westerly winds and invigorated ocean circulation. The associated northward migration of Southern Ocean fronts has been linked with reduced Atlantic Meridional Overturning Circulation by restricting surface water connectivity between the ocean basins, with implications for heat transport to the high latitudes of the North Atlantic. While our results do not exclude low-latitude mechanisms as drivers for Pliocene cooling, they indicate an additional role played by southern high-latitude cooling during development of the bipolar world.

(Table S1) Radiocarbon dates of organic matter from Lena Delta ice wedges

Relative to the past 2,000 years, the Arctic region has warmed significantly over the past few decades. However, the evolution of Arctic temperatures during the rest of the Holocene is less clear. Proxy reconstructions, suggest a long-term cooling trend throughout the mid- to late Holocene, whereas climate model simulations show only minor changes or even warming. Here we present a record of the oxygen isotope composition of permafrost ice wedges from the Lena River Delta in the Siberian Arctic. The isotope values, which reflect winter season temperatures, became progressively more enriched over the past 7,000 years, reaching unprecedented levels in the past five decades. This warming trend during the mid- to late Holocene is in opposition to the cooling seen in other proxy records. However, most of these existing proxy records are biased towards summer temperatures. We argue that the opposing trends are related to the seasonally different orbital forcing over this interval. Furthermore, our reconstructed trend as well as the recent maximum are consistent with the greenhouse gas forcing and climate model simulations, thus reconciling differing estimates of Arctic and northern high-latitude temperature evolution during the Holocene.