(Table 1) Age determination of sediment core BC5, Golf of Lion

Hitherto unknown abundance peaks of left coiling (l.c.) Neogloboquadrina pachyderma from a Gulf of Lions piston core indicate that abrupt cold spells associated with Atlantic Heinrich events affected the Mediterranean. N. pachyderma (l.c.) is typical of (sub) polar waters in the open ocean. The southern edge of its glacial North Atlantic bioprovince reached south Portugal. Only trace abundances of N. pachyderma (l.c.) are known from Quaternary Mediterranean sediments, suggesting that no significant "invasions" occured via the Strait of Gibraltar. The Gulf of Lions abundance peaks therefore seem to reflect area-specific thriving of a normally rare but indigenous taxon in the western Mediterranean through local favorable habitat development. The general planktonic foraminiferal record suggests that the basic hydrographic regime in the Gulf of Lions, with wintertime deep convective overturn, was relatively stable over the past 60 kyr. Under these conditions, high abundances of N. pachyderma (l.c.) would essentially imply temperature reductions of the order of 5°-8° relative to the present.

(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 2: Age information for the analyzed cores

Five plankton groups, including diatoms, radiolarians, coccolithophores, foraminifers, and dinoflagellate cysts, were synoptically analyzed in six sediment cores and two sediment traps from the Norwegian-Greenland Sea and the North Atlantic in order to provide more detailed insights into the paleoclimatic and paleoceanographic evolution and the development of plankton assemblages of the northern North Atlantic during the last 15,000 years. Based on Q-mode factor analyses, cold, warm, transitional, and relict assemblages were calculated for each of the plankton groups. Data from the different plankton groups complement one another, although they are not always consistent. However, the multiple plankton-group data set is able to bridge intervals in which single groups lack preservation or the ability to react to changes. Synoptically interpreted, the results provide a detailed picture of the response of plankton assemblages to environmental changes during the time period investigated, which includes the B0lling/Aller0d interstadial, the Younger Dryas cold spell, Termination IB, and, in all likelihood, also the "8,200 Event", and the Hypsithermal (approximately 8-4 14C ky BP).

(Table 1) Age comparison between the GPTS and ODP Site 181-1123 age model for each chron base within the Middle Miocene section

Today the deep western boundary current (DWBC) east of New Zealand is the most important route for deep water entering the Pacific Ocean. Large-scale changes in deep water circulation patterns are thought to have been associated with the development of the East Antarctic Ice Sheet (EAIS) close to the main source of bottom water for the DWBC. Here we reconstruct the changing speed of the southwest Pacific DWBC during the middle Miocene from ~15.5-12.5 Ma, a period of significant global ice accumulation associated with EAIS growth. Sortable silt mean grain sizes from Ocean Drilling Program Site 1123 reveal variability in the speed of the Pacific inflow on the timescale of the 41 kyr orbital obliquity cycle. Similar orbital period flow changes have recently been demonstrated for the Pleistocene epoch. Collectively, these observations suggest that a strong coupling between changes in the speed of the deep Pacific inflow and high-latitude climate forcing may have been a persistent feature of the global thermohaline circulation system for at least the past 15 Myr. Furthermore, long-term changes in flow speed suggest an intensification of the DWBC under an inferred increase in Southern Component Water production. This occurred at the same time as decreasing Tethyan outflow and major EAIS growth between ~15.5 and 13.5 Ma. These results provide evidence that a major component of the deep thermohaline circulation was associated with the middle Miocene growth of the EAIS and support the view that this time interval represents an important step in the development of the Neogene icehouse climate.

(Table S1) Tie-points used in age model construction for ODP Leg 172 sites

ix Ocean Drilling Program (ODP) sites, in the Northwest Atlantic have been used to investigate kinematic and chemical changes in the "Western Boundary Undercurrent" (WBUC) during the development of full glacial conditions across the Marine Isotope Stage 5a/4 boundary (~70,000 years ago). Sortable silt mean grain size(sort s) measurements are employed to examine changes in near bottom flow speeds, together with carbon isotopes measured in benthic foraminifera and % planktic foraminiferal fragmentation as proxies for changes in water-mass chemistry. A depth transect of cores, spanning 1.8-4.6 km depth, allows changes in both the strength and depth of the WBUC to be constrained across millennial scale events. Sort s measurements reveal that the flow speed structure of the WBUC during warm intervals ("interstadials") was comparable to modern (Holocene) conditions. However, significant differences are observed during cold intervals, with higher relative flow speeds inferred for the shallow component of the WBUC (~2 km depth) during all cold "stadial" intervals (including Heinrich Stadial 6), and a substantial weakening of the deep component (~3-4 km) during full glacial conditions. Our results therefore reveal that the onset of full glacial conditions was associated with a regime shift to a shallower mode of circulation (involving Glacial North Atlantic Intermediate Water) that was quantitatively distinct from preceding cold stadial events. Furthermore, our chemical proxy data show that the physical response of the WBUC during the last glacial inception was probably coupled to basin-wide changes in the water-mass composition of the deep Northwest Atlantic.

(Table 1) Age markers and errors for the last 1215 ka in IODP Exp302

Despite its importance in the global climate system, age-calibrated marine geologic records reflecting the evolution of glacial cycles through the Pleistocene are largely absent from the central Arctic Ocean. This is especially true for sediments older than 200 ka. Three sites cored during the Integrated Ocean Drilling Program's Expedition 302, the Arctic Coring Expedition (ACEX), provide a 27 m continuous sedimentary section from the Lomonosov Ridge in the central Arctic Ocean. Two key biostratigraphic datums and constraints from the magnetic inclination data are used to anchor the chronology of these sediments back to the base of the Cobb Mountain subchron (1215 ka). Beyond 1215 ka, two best fitting geomagnetic models are used to investigate the nature of cyclostratigraphic change. Within this chronology we show that bulk and mineral magnetic properties of the sediments vary on predicted Milankovitch frequencies. These cyclic variations record ''glacial'' and ''interglacial'' modes of sediment deposition on the Lomonosov Ridge as evident in studies of ice-rafted debris and stable isotopic and faunal assemblages for the last two glacial cycles and were used to tune the age model. Potential errors, which largely arise from uncertainties in the nature of downhole paleomagnetic variability, and the choice of a tuning target are handled by defining an error envelope that is based on the best fitting cyclostratigraphic and geomagnetic solutions.