(Table 1) Stable isotope ratios of Globigerinoides ruber and pore water, and Cl measurement at ODP Hole 133-817C

The oxygen isotopic composition of pore waters squeezed from sediments in Hole 817C co-varies with the oxygen isotopic composition of Globigerinoides ruber below 8 mbsf. The magnitude of the variation in the pore water d18O is approximately 30% of the variation in the foraminifers. Overall, the d18O of the pore waters increases down the core, a trend that is also present in the Cl- concentrations. The variations in the d18O of pore waters may be the result of either of two phenomena. First, these may reflect original variations in the waters, the magnitude of which has been subsequently reduced by process of diffusion. Second, these may reflect recrystallization of the precursor sediment and isotopic exchange between the fluids and the recrystallized sediment. At the moment data are not available to ascertain which process is responsible although the correlation between the Cl- and the d18O data suggests that these values reflect the original composition modified by diffusion.

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

A 5 million year reconstruction of sea level, temperature, and sea water d18O

Marine sediment records from the Oligocene and Miocene reveal clear 400,000-year (400-kyr) climate cycles related to variations in orbital eccentricity. These cycles are also observed in the Plio-Pleistocene records of the global carbon cycle. However they are absent in the Late Pleistocene ice-age record over the past 1.5 million years. Here, we present a simulation of global ice volume over the past 5 million years with a coupled system of four 3-D ice-sheet models. Our simulation shows that the 400-kyr long eccentricity cycles of Antarctica vary coherently with d13C records during the Pleistocene suggesting that they drive the long-term carbon cycle changes throughout the past 35 million years. The 400-kyr response of Antarctica is eventually suppressed by the dominant 100-kyr glacial cycles of the large ice sheets in the Northern Hemisphere (NH).

(Table 1) Age control points and of sediment core MD01-2444

The marine isotopic stage 3 (MIS3) at Ocean Drilling Program (ODP) Site 1060 (Gulf Stream) shows both sharp onset and end of interstadials, the existence of very short lived warm events during stadials, and points to differences in detail between the sea surface temperature (SST) record from the western North Atlantic and the atmospheric temperature record inferred from d18O in Greenland ice. Investigating MIS3 and obtaining comparable data from other locations appears crucial. The eastern Atlantic provides well-documented records of climate changes. We have selected a core from off Portugal and use it to examine Dansgaard/Oeschger events (D/O) at centennial-scale resolution (139 years on average between two data points). We have obtained a faunal data set for core MD01-2444, 37°N, 10°W, 2600 m water depth and use a group of species (Globigerina bulloides + Globigerinita glutinata) as a proxy of upwelling intensity driven by trade winds intensity changes. We tentatively relate the variation of this group to a North Atlantic Oscillation-like phenomenon (NAO) off Portugal. We observe that it resembles the rainfall index in the Caribbean as recorded at ODP Site 1002 (Cariaco Basin) which traces the Intertropical Convergence Zone (ITCZ) location through changes of terrigenous inputs. The driest intervals (ITZC to the south) at Site 1002 correspond to intervals of increased upwelling in MD01-2444 as well as the driest periods identified during stadials on similar cores in the area. Because the ITZC to the south is consistent with an El Niño-Southern Oscillation (ENSO+) situation, our study suggests a positive correlation between ENSO-like conditions and NAO-like conditions at a millennial timescale. During interstadial intervals when increased wetness over Cariaco is recorded (ITCZ to the north) and the upwelling in MD01-2444 is decreased, we see from both SSTs and faunal tropical indicators that MD01-2444 site is warm. In addition, interstadials are equally warm through each so-called Bond cycle. This contrasts with the Greenland Ice Core Project (GRIP) record where interstadial peaks are successively cooler through each Bond cycle. This record confirms a link between tropical climate linked to ITCZ position and the climate of southern Europe at millennial timescales, in spite of showing a very good correlation with polar latitudes (GRIP) through d18O on Globigerina bulloides. In addition, because the warmest SSTs and the d18O on G. bulloides are so remarkably different, our work points to changes in seasonality as a strong control over the climatic pattern of the North Atlantic area and the marked influence of winter conditions.

(Table 1) Age determination of sediment core PC17

We have investigated glacial-interglacial differences in sea surface temperature (SST) near Hawaii using two relatively high deposition rate, shallow-water piston cores collected near Oahu, Hawaii. Modern hydrographic data show that local surface water temperatures are broadly consistent with the regional pattern of SSTs in the southern subtropical North Pacific. Past SSTs were estimated on the basis of three independently measured parameters: (1) UK'37 values of alkenones, (2) d18O of Globigerinoides ruber, and (3) assemblages of planktonic foraminifera using the modern analog technique (MAT). The two cores yield similar SST records, and if differences in the ecology of foraminifera and coccolithophores are considered, the three different approaches to estimating SSTs yield consistent results. UK'37-based temperatures, which may represent winter values at this location, were ~2.5°C colder during the Last Glacial Maximum than today, which is consistent with the February MAT estimates. The d18O-based temperature estimates, likely biased toward summer temperatures, indicate that the glacial SSTs were at least 1°C cooler than today, which is comparable to the results of MAT August estimates.

Age determination of sediment core GIK16523-1

High-resolution studies of a planktonic foraminifer core record from the South China Sea (SCS) (31KL: 18°45.4'N, 115°52.4'E, water depth 3360 m) reveal changes driven by ice-volume forcings in the climate of the East Asian monsoon in the western Pacific marginal sea during the late Quaternary. The analyses of planktonic foraminifer faunal abundance data from the core indicate significant variations in the relative abundances of the dominant taxa over the past 100,000 years since the isotope stage 5. The transfer function estimates of faunal sea surface temperatures (SST) correlate well with a long-term (104-105 years) trend of global glaciation. About 65,000 years ago, there was an observable change in the mode of SST variability as many low-latitude records have shown. These findings suggest that the SCS surface circulation and the East Asian winter monsoon systems are mainly driven by variations in global glaciation levels. The association of surface ocean cooling in the SCS with global climatic events suggests that fluctuations in the strength of the East Asian winter monsoon may be linked to shifts in the latitudinal position of the westerly winds and the Siberian high-pressure system.

Sedimentology, age models and stable isotope ratios of sediment cores from the equatorial Pacific

Based on detailed reconstructions of global distribution patterns, both paleoproductivity and the benthic d13C record of CO2, which is dissolved in the deep ocean, strongly differed between the Last Glacial Maximum and the Holocene. With the onset of Termination I about 15,000 years ago, the new (export) production of low- and mid-latitude upwelling cells started to decline by more than 2-4 Gt carbon/year. This reduction is regarded as a main factor leading to both the simultaneous rise in atmospheric CO2 as recorded in ice cores and, with a slight delay of more than 1000 years, to a large-scale gradual CO2 depletion of the deep ocean by about 650 Gt C. This estimate is based on an average increase in benthic d13C by 0.4-0.5 per mil. The decrease in new production also matches a clear 13C depletion of organic matter, possibly recording an end of extreme nutrient utilization in upwelling cells. As shown by Sarnthein et al., [1987], the productivity reversal appears to be triggered by a rapid reduction in the strength of meridional trades, which in turn was linked via a shrinking extent of sea ice to a massive increase in high-latitude insolation, i.e., to orbital forcing as primary cause.

(Table 2) Age determination of sediment cores from the Gulf of Cadiz

We present micropalaeontological and grain-size records for a set of sediment cores from the Gulf of Cadiz (southwest Spain) that reflect changes in the position and strength of the Mediterranean Outflow (MO) current. The cores sample a sediment drift (the Gil Eanes Drift) that is positioned lower on the slope in the Gulf of Cadiz than the position of the main current today. The data indicate that the drift is of glacial age and that the glacial MO current was positioned lower on the slope than today but also that it was active over a considerably reduced area of the slope. We argue that this observation is consistent with physical constraints on the Gibraltar Exchange and on the likely settling and spreading behavior of the MO plume along the Iberian Margin under glacial environmental and sea level conditions. The deeper settling of the MO is likely to have influenced the formation of glacial North Atlantic Intermediate Water and also may have exerted indirect influence on the formation of glacial North Atlantic Deep Water.

(Table 1) Age control and temporal resolution of paleoceanographic records of the North Atlantic

Short-term changes in sea surface conditions controlling the thermohaline circulation in the northern North Atlantic are expected to be especially efficient in perturbing global climate stability. Here we assess past variability of sea surface temperature (SST) in the northeast Atlantic and Norwegian Sea during Marine Isotope Stage (MIS) 2 and, in particular, during the Last Glacial Maximum (LGM). Five high-resolution SST records were established on a meridional transect (53°N-72°N) to trace centennial-scale oscillations in SST and sea-ice cover. We used three independent computational techniques (SIMMAX modern analogue technique, Artificial Neural Networks (ANN), and Revised Analog Method (RAM)) to reconstruct SST from planktonic foraminifer census counts. SIMMAX and ANN reproduced short-term SST oscillations of similar magnitude and absolute levels, while RAM, owing to a restrictive analog selection, appears less suitable for reconstructing "cold end" SST. The SIMMAX and ANN SST reconstructions support the existence of a weak paleo-Norwegian Current during Dansgaard-Oeschger (DO) interstadials number 4, 3, 2, and 1. During the LGM, two warm incursions of 7°C water to occurred in the northern North Atlantic but ended north of the Iceland Faroe Ridge. A rough numerical estimate shows that the near-surface poleward heat transfer from 53° across the Iceland-Faroe Ridge up to to 72° N dropped to less than 60% of the modern value during DO interstadials and to almost zero during DO stadials. Summer sea ice was generally confined to the area north of 70°N and only rarely expanded southward along the margins of continental ice sheets. Internal LGM variability of North Atlantic (>40°N) SST in the GLAMAP 2000 compilation (Sarnthein et al., 2003, doi:10.1029/2002PA000771; Pflaumann et al., 2003, doi:10.1029/2002PA000774) indicates maximum instability in the glacial subpolar gyre and at the Iberian Margin, while in the Nordic Seas, SST was continuously low.

(Table 2) Age determination of sediment cores from the Portuguese Margin

We demonstrate that changes in the behavior of the Mediterranean Outflow Water (MOW) prior to and through the last deglaciation played an important role in promoting Meridional Overturning Circulation (MOC). Estimation of past MOW salt and heat fluxes indicates that they gradually increased through the last deglaciation. Between 17.5 and 14.6 thousand years ago (ka B.P., where B.P. references year 1950), net evaporation from the Mediterranean exported sufficient fresh water from the North Atlantic catchment to cause an average salinity increase of 0.5 psu throughout the upper 2000 m of the entire North Atlantic to the north of 25°N. Combined with rapid intensification and shoaling of the MOW plume, which we identify around 15-14.5 ka B.P., this deglacial MOW-related salt accumulation preconditioned the North Atlantic for abrupt resumption of the MOC at 14.6 ka B.P.

Age determination of sediment cores from the Laptev Sea

To establish a chronology of the Holocene transgression in Arctic Siberia, a total of 14 sediment cores from the Laptev Sea continental slope and shelf were studied covering the water depth range between 983 and 21 m. The age models of the cores were derived from 119 radiocarbon datings, which were all analyzed on marine biogenic calcite (mainly bivalve shells). The oldest shell sample was found at the slope and dates back to about 15.3 cal. ka, indicating that the time interval investigated starts prior to the onset of the meltwater pulse 1A (~14.2 cal. ka) when global sea-level rose dramatically. The inundation history was reconstructed mainly on the basis of major changes in average sedimentation rates (ASR), but also other sedimentological parameters were incorporated. A diachronous reduction in ASR from the outer to the inner shelf region is recognized, which was related to the southward migration of the coastline as the primary sediment source. We estimate that the flooding of the 50-, 43-, and 31-m isobaths was completed by approximately 11.1, 9.8, and 8.9 cal. ka, and that Holocene sea-level highstand was approached near 5 cal. ka. Between these time intervals, sea level in the Laptev Sea rose by 5.4, 13.3, and 7.9 mm/year, respectively.