(Table 2) Age determination of sediment core NIOP-C2_929

In this study we present a sea surface temperature (SST) record from the western Arabian Sea for the last 20,000 years. We produced centennial-scale d18O and Mg/Ca SST time series of core NIOP929 with focus on the glacial-interglacial transition. The western Arabian Sea is influenced by the seasonal NE and SW monsoon wind systems. Lowest SSTs occur during the SW monsoon season because of upwelling of cold water, and highest SSTs can be found in the low-productivity intermonsoon season. The Mg/Ca-based temperature record reflects the integrated SST of the SW and NE monsoon seasons. The results show a glacial-interglacial SST difference of ~2°C, which is corroborated by findings from other Arabian Sea cores. At 19 ka B.P. a yet undescribed warm event of several hundred years duration is found, which is also reflected in the d18O record. A second centennial-scale high SST/low d18O event is observed at 17 ka B.P. This event forms the onset of the stepwise yet persistent trend toward Holocene temperatures. Highest Mg/Ca-derived SSTs in the NIOP929 record occurred between 13 and 10 ka B.P. Interglacial SST is ~24°C, indicating influence of upwelling. The onset of Arabian Sea warming occurs when the North Atlantic is experiencing minimum temperatures. The rapid temperature variations at 19, 17, and 13 ka B.P. are difficult to explain with monsoon changes alone and are most likely also linked to regional hydrographic changes, such as trade wind induced variations in warm water advection.

Age determination and stable carbon isotope ratios of Neogloboquadrina pachyderma of IODP Hole 302-M0004C

The Integrated Ocean Drilling Program (IODP) Arctic Coring Expedition (ACEX) Hole 4C from the Lomonosov Ridge in the central Arctic Ocean recovered a continuous 18 m record of Quaternary foraminifera yielding evidence for seasonally ice-free interglacials during the Matuyama, progressive development of large glacials during the mid-Pleistocene transition (MPT) ~1.2-0.9 Ma, and the onset of high-amplitude 100-ka orbital cycles ~500 ka. Foraminiferal preservation in sediments from the Arctic is influenced by primary (sea ice, organic input, and other environmental conditions) and secondary factors (syndepositional, long-term pore water dissolution). Taking these into account, the ACEX 4C record shows distinct maxima in agglutinated foraminiferal abundance corresponding to several interglacials and deglacials between marine isotope stages (MIS) 13-37, and although less precise dating is available for older sediments, these trends appear to continue through the Matuyama. The MPT is characterized by nearly barren intervals during major glacials (MIS 12, 16, and 22-24) and faunal turnover (MIS 12-24). Abundant calcareous planktonic (mainly Neogloboquadrina pachyderma sin.) and benthic foraminifers occur mainly in interglacial intervals during the Brunhes and very rarely in the Matuyama. A distinct faunal transition from calcareous to agglutinated foraminifers 200-300 ka in ACEX 4C is comparable to that found in Arctic sediments from the Lomonosov, Alpha, and Northwind ridges and the Morris Jesup Rise. Down-core disappearance of calcareous taxa is probably related to either reduced sea ice cover prior to the last few 100-ka cycles, pore water dissolution, or both.

Age determinations and stable isotope records of four sediment profiles off Nova Scotia

Continental margin sediments off Nova Scotia accumulate at high rates (up to 360 cm/kyr) and contain a history of millennial-scale environmental changes which are dominated by the proximity of the Laurentide ice sheet during the latest Quaternary. Using stable isotope ratios of oxygen, accelerator mass spectrometer radiocarbon dating, micropaleontology, and sedimentology, we document these changes in six piston cores ranging in water depth from ab. 450 to ab. 4300 m. We find that maximum d18O in N. pachyderma occurred about 15 ka and preceded the maximum abundance of this species in these cores by ab. 1000 years. Between 13 and 14 ka we find a second peak in abundance of N. pachyderma, minimum d18O, and two pulses of ice rafting. The sediment lithology supports terrestrial studies which indicate that there was a general withdrawal of ice beyond the upper Paleozoic and Mesozoic red beds by 14 ka in southeastern Canada, so the ice rafting events between 13 and 14 ka probably reflect ice stream activity in the St. Lawrence valley. The Younger Dryas event is recognized as a peak in abundance of N. pachyderma and ice rafting (dated as ab. 11.3 ka), but meltwater discharge to the Gulf of St. Lawrence was either too small or occurred over too long a time to leave a distinct d18O minimum off Nova Scotia. At 7.1 ka, in the middle of Holocene warming, we find a third peak in abundance of N. pachyderma and another d18O minimum but no ice rafting. We interpret these data as evidence of a late-occurring meltwater event which, if correct, could have originated in the Great Lakes, in the Labrador-Ungava region, or in both. The final millennial-scale phenomenon off Nova Scotia is the onset of "Neoglaciation," marked by increased ice rafting and increased % N. pachyderma beginning about 5 kyr ago.

(Table 2) Age determination of Florida Straits sediment cores

The waters passing through the Florida Straits today reflect both the western portion of the wind-driven subtropical gyre and the northward flow of the upper waters which cross the equator, compensating North Atlantic Deep Water export as part of the large-scale Atlantic meridional overturning circulation. It has been postulated from various lines of evidence that the overturning circulation was weaker during the Younger Dryas cold event of the last deglaciation. We show here that the contrast in the oxygen isotopic composition of benthic foraminiferal tests across the Florida Current is reduced during the Younger Dryas. This most likely reflects a decrease in the density gradient across the channel and a decrease in the vertical shear of the Florida Current. This reduced shear is consistent with the postulated reduction in the Atlantic meridional overturning circulation. We find that the onset of this change in density structure and flow at the start of the Younger Dryas is very abrupt, occurring in less than 70 years.

(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 2) Age determination of Aegean Sea sediments

The modern Aegean Sea is an important source of deep water for the eastern Mediterranean. Its contribution to deep water ventilation is known to fluctuate in response to climatic variation on a decadal timescale. This study uses marine micropalaeontological and stable isotope data to investigate longer-term variability during the late glacial and Holocene, in particular that associated with the deposition of the early Holocene dysoxic/anoxic sapropel S1. Concentrating on the onset of sapropel-forming conditions, we identify the start of 'seasonal' stratification and highlight a lag in d18O response of the planktonic foraminifer N. pachyderma to termination T1b as identified in the d18O record of G. ruber. By use of a simple model we determine that this offset cannot be a function of bioturbation effects. The lag is of the order of 1 kyr and suggests that isolation of intermediate/deep water preceded the start of sapropel formation by up to 1.5 kyr. Using this discovery, we propose an explanation for the major unresolved problem in sapropel studies, namely, the source of nutrient supply required for export productivity to reach levels needed for sustained sapropel deposition. We suggest that nutrients had been accumulating in a stagnant basin for 1-1.5 kyr and that these accumulated resources were utilized during the deposition of S1. In addition, we provide a first quantitative estimate of the diffusive (1/e) mixing timescale for the eastern Mediterranean in its "stratified" sapropel mode, which is of the order of 450 years.

(Table 1) Age determination of sediment core AAS-62/2

Variations in the Indian summer monsoon (ISM) intensity during the last 16.7?ka have been studied using organic carbon (Corg), d15N of sedimentary organic matter, CaCO3, sediment texture, d18OC, and Mg/Ca-derived sea surface temperature, d18O of sea water and sea surface salinity, in a 14C-dated sediment core from the eastern Arabian Sea. The d18O in water and planktonic foraminifera shells off the central west coast of India may be controlled by the ISM intensity as this area receives high precipitation and land runoff. Also, the Corg and CaCO3 contents of sediments and d15N of sedimentary organic matter may be linked to ISM-induced productivity and denitrification. The results of the present study reveal that between 16 and 15.2 ka BP, the ISM was weak with minor fluctuations and started intensifying around 15.2 ka BP, at the onset of the Bølling-Ållerød (B-A) event. The B-A event is characterized by higher water column denitrification rates comparable to the present day. The ISM signatures observed in the d18OC record of B-A event compare well with those from Timta cave of the western Himalayas and also the Asian summer monsoon signatures from the Hulu caves in China and warming signatures in Greenland Ice Sheet Project 2 (GISP2) suggesting atmospheric teleconnections through Intertropical Convergence Zone. The boundary between the Younger Dryas and the Holocene is discernible with small episodes of abrupt events of increased ISM intensity. This decrease in d18OC values at ~11.8 ka BP is contemporary with June solar insolation maximum at 30° north and the increase in methane in the GISP2 ice core supporting episodes of warmer climate and increase in ISM intensity. The ISM seems to have been most stable between 7 and 5.6 ka BP. The core exhibits periodicity of 500 years that is comparable to the Atlantic water formation and the Chinese monsoon.

Age determination of sediment core NP90-52 (Table 1)

Data on glacial erosion have been compiled and synthesised using a wide range of sediment budget and sediment yield studies from the Svalbard-Barents Sea region. The data include studies ranging in timescale from 1 to 10**6 yr, and in size of drainage basin from 101 to 105 km**2. They show a clear dependence of sediment yield on the mode of glacierization. Polar glaciers erode at rates comparable to those found in Arctic fluvial basins, or about 40 t/km**-2/ yr or 0.02 mm/yr. In contrast, rates of erosion by polythermal glaciers are 800-1000 t/km**2/ yr (or ca 0.3-0.4 mm/yr), while rates from fast-flowing glaciers are slightly more than twice this: 2100 t/km**2/yr (or 1 mm/yr). Similar rates are also found for large glacierized basins like those in the southwestern parts of the Barents Sea. In contrast to the situation in fluvial basins, in which sediment yield typically decreases with increasing basin size, the tendency in glacierized basins is for erosion to be independent of basin size. In studies of sediment yield from glaciers it is sometimes difficult to distinguish between material actually dislodged from the bedrock by glaciers and material dislodged by other processes in interglacial times and simply transported to a depocenter by a glacier. Our data suggest that pulses of sediment resulting from advance of a glacier over previously-dislodged material last on the order of 10**3 yr, and result in inferred erosion rates that are approximately 25% higher than long-term average rates of glacial erosion. The maximum sediment fluxes from the large Storfjorden and Bear Island drainage basins occurred in mid-Pleistocene. The onset of this period of high sediment yield coincided with the onset of the 100 kyr glacial cycle. We presume that this was the beginning of a period of increased glacial activity, but one in which glaciers still advanced and retreated frequently. During the last two to four 100 kyr cycles, however, sediment yields appear to have decreased. This decrease may be the result of the submergence of the Barents Sea. Glacier erosion would be much higher for a subaerial Barents Sea setting than it would be for a present day subsea Barents Sea. A classical question in Quaternary Geology is whether glaciers are more erosive than rivers. We surmise that if factors such as the lithology and the available potential energy (mgh) of the precipitation falling at a given altitude, whether in liquid or solid form, are held constant, then glaciers are vastly more effective agents of erosion than rivers.

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

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.

Diatom oxygen and silicon isotopes data from IODP Hole 323-U1341B in the Bering Sea through the late Pliocene and over the onset of major Northern Hemisphere Glaciation

The Bering Sea is one of the most biologically productive regions in the marine system and plays a key role in regulating the flow of waters to the Arctic Ocean and into the subarctic North Pacific Ocean. Cores from Integrated Ocean Drilling Program (IODP) Expedition 323 to the Bering Sea provide the first opportunity to obtain reconstructions from the region that extend back to the Pliocene. Previous research at Bowers Ridge, south Bering Sea, has revealed stable levels of siliceous productivity over the onset of major Northern Hemisphere Glaciation (NHG) (circa 2.85-2.73 Ma). However, diatom silica isotope records of oxygen (d18Odiatom) and silicon (d30Sidiatom) presented here demonstrate that this interval was associated with a progressive increase in the supply of silicic acid to the region, superimposed on shift to a more dynamic environment characterized by colder temperatures and increased sea ice. This concluded at 2.58 Ma with a sharp increase in diatom productivity, further increases in photic zone nutrient availability and a permanent shift to colder sea surface conditions. These transitions are suggested to reflect a gradually more intense nutrient leakage from the subarctic northwest Pacific Ocean, with increases in productivity further aided by increased sea ice- and wind-driven mixing in the Bering Sea. In suggesting a linkage in biogeochemical cycling between the south Bering Sea and subarctic Northwest Pacific Ocean, mainly via the Kamchatka Strait, this work highlights the need to consider the interconnectivity of these two systems when future reconstructions are carried out in the region.

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.