(Table S1) Brominated flame retardant concentrations in an ice core from Holtedahlfonna, Svalbard

Brominated flame retardants (BFRs) have been found in Arctic wildlife, lake sediment, and air. To identify the atmospheric BFR deposition history on Svalbard, Norway, we analyzed 19 BFRs, including hexabromocyclododecane (HBCD), 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE), decabromodiphenyl ethane (DBDPE), pentabromoethylbenzene (PBEB),and 15 polybrominated diphenyl ether congeners (PBDE) in the upper 34 m of an ice core (representing 1953-2005) from Holtedahlfonna, the western-most ice sheet on Svalbard. All of the non-PBDE compounds were detected in nearly continuous profiles in the core. Seven PBDEs were not observed above background (28,47,66,100,99,154,153), while 4 were found in 1 or 2 of 6 segments (17,85,138,183). BDEs-49,71,190,209 had nearly continuous profiles but only BDE-209 in large amounts. The greatest inputs were HBCD and BDE-209, 910, and 320 pg/cm**2/yr from 1995-2005. DBDPE, BTBPE, and PBEB show nearly continuous input growth in recent core segments, but all were <6 pg/cm**2/yr. Long-range atmospheric processes may have moved these particle-bound BFRs to the site, probably during the Arctic haze season. Average air mass trajectories over 10 years show >75% of atmospheric flow to Holtedahlfonna coming from Eurasia during haze periods (March and April).

(Table S1) Stable carbon and oxygen isotope ratios of benthic foraminifera of the western Atlantic Ocean

Oxygen and carbon isotopic data were produced on the benthic foraminiferal taxa Cibicidoides and Planulina from 25 new piston cores, gravity cores, and multicores from the Brazil margin. The cores span water depths from about 400 to 3000 m and intersect the major water masses in this region. These new data fill a critical gap in the South Atlantic Ocean and provide the motivation for updating the classic glacial western Atlantic d13C transect of Duplessy et al. (1988). The distribution of 13C of SumCO2 requires the presence of three distinct water masses in the glacial Atlantic Ocean: a shallow (~1000 m), southern source water mass with an end-member d13C value of about 0.3-0.5 per mil VPDB, a middepth (~1500 m), northern source water mass with an end-member value of about 1.5 per mil, and a deep (>2000 m), southern source water with an end-member value of less than -0.2 per mil, and perhaps as low as the -0.9 per mil values observed in the South Atlantic sector of the Southern Ocean (Ninnemann and Charles, 2002, doi:10.1016/S0012-821X(02)00708-2). The origins of the water masses are supported by the meridional gradients in benthic foraminiferal d18O. A revised glacial section of deep water d13C documents the positions and gradients among these end-member intermediate and deep water masses. The large property gradients in the presence of strong vertical mixing can only be maintained by a vigorous overturning circulation.

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