(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.

(Table S1) Age determination of sediment cores from the Labrador Sea

High resolution flow speed reconstructions of two core sites located on Gardar Drift in the northeast Atlantic Basin and Orphan Knoll in the northwest Atlantic Basin reveal a long-term decrease in flow speed of Northeast Atlantic Deep Water (NEADW) after 6,500 years. Benthic foraminiferal oxygen isotopes of sites currently bathed in NEADW show a 0.2per mil depletion after 6,500 years, shortly after the start of the development of a carbon isotope gradient between NEADW and Norwegian Sea Deep Water. We consider these changes in near-bottom flow vigor and benthic foraminiferal isotope records to mark a significant reorganization of the Holocene deep ocean circulation, and attribute the changes to a weakening of NEADW flow during the mid to late Holocene that allowed the shoaling of Lower Deep Water and deeper eastward advection of Labrador Sea Water into the northeast Atlantic Basin.

(Table S1) TEX86 and BIT record from sediment core PBBC-1

A molecular organic geochemical proxy (TEX86) for sea surface temperature (SST) is compared with a foraminifera-based SST proxy (Mg/Ca) in a decadal-resolution marine sedimentary record spanning the last 1000 years from the Gulf of Mexico. We assess the relative strengths of the organic and inorganic paleoceanographic techniques for reconstructing high-resolution SST variability during recent climate events, including the Little Ice Age (LIA) and the Medieval Warm Period (MWP). SST estimates based on the molecular organic proxy TEX86 show a similar magnitude and pattern of SST variability to foraminiferal Mg/Ca-SST estimates but with some important differences. For instance, both proxies show a cooling (1°C-2°C) of Gulf of Mexico SSTs during the LIA. During the MWP, however, Mg/Ca-SSTs are similar to near-modern SSTs, while TEX86 indicates SSTs that were cooler than modern. Using the respective SST calibrations for each proxy results in TEX86-SST estimates that are 2°C-4°C warmer than Mg/Ca-SST throughout the 1000 year record. We interpret the TEX86-SST as a summer-weighted SST signal from the upper mixed layer, whereas the Mg/Ca-SST better reflects the mean annual SST. Downcore differences in the SST estimates between the two proxies (DeltaT = TEX86 - Mg/Ca) are interpreted in the context of varying seasonality and/or changing water column temperature gradients.

Data S1. Results of the analyses considering developmental time, supplementary figures and tables & Video S1

Our ability to project the impact of global change on marine ecosystem is limited by our poor understanding on how to predict species sensitivity. For example, the impact of ocean acidification is highly species-specific, even in closely related taxa. The aim of this study was to test the hypothesis that the tolerance range of a given species to decreased pH corresponds to their natural range of exposure. Larvae of the green sea urchin Strongylocentrotus droebachiensis were cultured from fertilization to metamorphic competence (29 days) under a wide range of pH (from pHT = 8.0/pCO2 ~ 480 ?atm to pHT = 6.5/pCO2 ~ 20 000 ?atm) covering present (from pHT 8.7 to 7.6), projected near-future variability (from pHT 8.3 to 7.2) and beyond. Decreasing pH impacted all tested parameters (mortality, symmetry, growth, morphometry and respiration). Development of normal, although showing morphological plasticity, swimming larvae was possible as low as pHT >= 7.0. Within that range, decreasing pH increased mortality and asymmetry and decreased body length (BL) growth rate. Larvae raised at lowered pH and with similar BL had shorter arms and a wider body. Relative to a given BL, respiration rates and stomach volume both increased with decreasing pH suggesting changes in energy budget. At the lowest pHs (pHT <= 6.5), all the tested parameters were strongly negatively affected and no larva survived past 13 days post fertilization. In conclusion, sea urchin larvae appeared to be highly plastic when exposed to decreased pH until a physiological tipping point at pHT = 7.0. However, this plasticity was associated with direct (increased mortality) and indirect (decreased growth) consequences for fitness.