Stable isotope, UK'37, SST, chlorin, alkenone and TOC data of sediment cores GIK17927-2 and GIK17928-3

Upwelling intensity in the South China Sea has changed over glacial-interglacial cycles in response to orbital-scale changes in the East Asian Monsoon. Here, we evaluate new multi-proxy records of two sediment cores from the north-eastern South China Sea to uncover millennial-scale changes in winter monsoondriven upwelling over glacial Terminations I and II. On the basis of U/Th-based speleothem chronology, we compare these changes with sediment records of summer monsoondriven upwelling east of South Vietnam. Ocean upwelling is traced by reduced (UK'37-based) temperature and increased nutrient and productivity estimates of sea surface water (d13C on planktic foraminifera, accumulation rates of alkenones, chlorins, and total organic carbon). Accordingly, strong winter upwelling occurred north-west of Luzon (Philippines) during late Marine Isotope Stage 6.2, Heinrich (HS) and Greenland stadials (GS) HS-11, GS-26, GS-25, HS-1, and the Younger Dryas. During these stadials, summer upwelling decreased off South Vietnam and sea surface salinity reached a maximum suggesting a drop in monsoon rains, concurrent with speleothem records of aridity in China. In harmony with a stadial-to-interstadial see-saw pattern, winter upwelling off Luzon in turn was weak during interstadials, in particular those of glacial Terminations I and II, when summer upwelling culminated east of South Vietnam. Most likely, this upwelling terminated widespread deep-water stratification, coeval with the deglacial rise in atmospheric CO2. Yet, a synchronous maximum in precipitation fostered estuarine overturning circulation in the South China Sea, in particular as long as the Borneo Strait was closed when sea level dropped below -40 m.

Rockmagnetic investigations of IODP Site U1302-U1303

Integrated Ocean Drilling Program (IODP) Sites U1302-U1303, drilled on the SE flank of Orphan Knoll (Labrador Sea), preserve a record of detrital layers and other proxies of hydrographic change that extend the record of ice-sheet/ocean interactions through most of the Brunhes Chron. The age model is built by tandem matching of relative paleointensity (RPI) and oxygen isotope data (d18O) from Neogloboquadrina pachyderma (sin.) to reference records, indicating a mean Brunhes sedimentation rate of 14 cm/kyr. Sedimentation back to marine isotope stage (MIS) 18 is characterized by detrital layers that are detected by higher than background gamma-ray attenuation (GRA) density, peaks in X-ray fluorescence (XRF) indicators for detrital carbonate (Ca/Sr) and detrital silicate (Si/Sr), and an ice-rafted debris (IRD) proxy (wt.% >106 µm). The age model enables correlation of Site U1302/03 to IODP Site U1308 in the heart of the central Atlantic IRD belt where an age model and a similar set of detrital-layer proxies have already been derived. Ages of Heinrich (H) layers H1, H2, H4, H5 and H6 are within ~2 kyr at the two sites (H0, H3 and H5a are not observed at Site U1308), and agree with previous work at Orphan Knoll within ~3 kyr. At Site U1308, Brunhes detrital layers are restricted to peak glacials and glacial terminations back to marine isotope stage (MIS) 16 and have near-synchronous analogs at Site U1302/03. Detrital layers at Site U1302/03 are distributed throughout the record in both glacial and most interglacial stages. We distinguish Heinrich-like layers associated with IRD from detrital layers marked by multiple detrital-layer proxies (including Ca/Sr) but usually not associated with IRD, that may be attributed to lofted sediment derived from drainage and debris-flow events funneled down the nearby Northwest Atlantic Mid-Ocean Channel (NAMOC). The prominent detrital layers at Sites U1302/03 and U1308 can be correlated to millennial scale features in the Chinese speleothem (monsoon) record over the last 400 kyr, implying a link between monsoon precipitation and Laurentide Ice Sheet (LIS) instability. The detrital-layer stratigraphy at Site U1302/03 provides a long record of LIS dynamics against which other terrestrial and marine records can be compared.

Sedimentology of core GeoTu_SL110 from the eastern Mediterranean Sea

Clay mineral assemblages in a sediment core from the distal Nile discharge plume off Israel have been used to reconstruct the late Quaternary Nile sediment discharge into the Eastern Mediterranean Sea (EMS). The record spans the last ca. 140 ka. Smectite abundances indicate the influence of the Blue Nile and Atbara that have their headwaters in the volcanic rocks of the Ethiopian highlands. Kaolinite abundances indicate the influence of wadis, which contribute periodically to the suspension load of the Nile. Due to the geographical position, the climate and the sedimentary framework of the EMS is controlled by two climate systems. The long-term climate regime was governed by the African monsoon that caused major humid periods with enhanced sediment discharge at 132 to <126 ka (AHP5), 116 to 99 ka (AHP4), and 89 to 77 ka (AHP3). They lasted much longer than the formation of the related sapropel layers S5 (>2 ka), S4 (3.5 ka) and S3 (5 ka). During the last glacial period (MIS 4-2) the long-term changes of the monsoonal system were superimposed by millennial-scale changes of an intensified mid-latitude glacial system. This climate regime caused short but pronounced drought periods in the Nile catchment, which are linked to Heinrich Events and alternate with more humid interstadials. The clay mineral record further implies that feedback mechanisms between vegetation cover and sediment discharge of the Nile are detectable but of minor importance for the sedimentary record in the southeastern Mediterranean Sea during the investigated African Humid Periods.

Age determination of sediment core W8709A-13

Northeast Pacific benthic foraminiferal d18O and d13 reveal repeated millennial-scale events of strong deep-sea ventilation (associated with nutrient depletion and/or high gas exchange) during stadial (cool, high ice volume) episodes from 10 to 60 ka, opposite the pattern in the deep North Atlantic. Two climate mechanisms may explain this pattern. North Pacific surface waters, chilled by atmospheric transmission from a cold North Atlantic and made saltier by reduced freshwater vapor transports, could have ventilated the deep Pacific from above. Alternatively, faster turnover of Pacific bottom and mid-depth waters, driven by Southern Ocean winds, may have compensated for suppressed North Atlantic Deep Water production during stadial intervals. During the Younger Dryas event (~11.6-13.0 cal ka), ventilation of the deep NE Pacific (~2700 m) lagged that in the Santa Barbara Basin (~450 m) by >500 years, suggesting that the NE Pacific was first ventilated at intermediate depth from above and then at greater depth from below. This apparent lag may reflect the adjustment time of global thermohaline circulation.

Radiocarbon age, Mg/Ca and d18O measurements on planktonic foraminifera of sediment core GeoB12605-3

The sea surface temperature (SST) of the tropical Indian Ocean is a major component of global climate teleconnections. While the Holocene SST history is documented for regions affected by the Indian and Arabian monsoons, data from the near-equatorial western Indian Ocean are sparse. Reconstructing past zonal and meridional SST gradients requires additional information on past temperatures from the western boundary current region. We present a unique record of Holocene SST and thermocline depth variations in the tropical western Indian Ocean as documented in foraminiferal Mg/Ca ratios and d18O from a sediment core off northern Tanzania. For Mg/Ca and thermocline d18O, most variance is concentrated in the centennial to bicentennial periodicity band. On the millennial time scale, an early to mid-Holocene (~7.8-5.6 ka) warm phase is followed by a temperature drop by up to 2°C, leading to a mid-Holocene cool interval (5.6-4.2 ka). The shift is accompanied by an initial reduction in the difference between surface and thermocline foraminiferal d18O, consistent with the thickening of the mixed layer and suggestions of a strengthened Walker circulation. However, we cannot confirm the expected enhanced zonal SST gradient, as the cooling of similar magnitude had previously been found in SSTs from the upwelling region off Sumatra and in Flores air temperatures. The SST pattern probably reflects the tropical Indian Ocean expression of a large-scale climate anomaly rather than a positive Indian Ocean Dipole-like mean state.

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

Sea surface temperature (TEX86), marine GDGTs and BIT index from ODP Hole 161-977A

The TEX86H temperature proxy is a relatively new proxy based on crenarchaeotal lipids and has rarely been applied together with other temperature proxies. In this study, we applied the TEX86H on a sediment core from the Alboran Sea (western Mediterranean, core ODP-977A) covering the penultimate climate cycle, that is, from 244 to 130 ka, and compared this with previously published sea surface temperatures derived from the Uk'37 of alkenones of haptophyta and Mg/Ca records of planktonic foraminifera. The TEX86H temperature record shows remarkably similar stadial-interstadial patterns and abrupt temperature changes to those observed with the Uk'37 palaeothermometer. Absolute TEX86H temperature estimates are generally higher than those of Uk'37, though this difference (<3°C in 81% of the data points) is mainly within the temperature calibration error for both proxies, suggesting that crenarchaeota and haptophyta experienced similar temperature variations. During occasional events (<5% of the analyzed time span), however, the TEX86H exhibits considerably higher absolute temperature estimates than the Uk'37. Comparison with Mg/Ca records of planktonic foraminifera as well as other Mediterranean TEX86 and Uk'37 records suggests that part of this divergence may be attributed to seasonal differences, that is, with TEX86H reflecting mainly the warm summer season while Uk'37 would show annual mean. Biases in the global calibration of both proxies or specific biases in the Mediterranean are an alternative, though less likely, explanation. Despite differences between absolute TEX86H and Uk'37 temperatures, the correlation between the two proxies (r**2 = 0.59, 95% significance) provides support for the occurrence of abrupt temperature variations in the western Mediterranean during the penultimate interglacial-to-glacial cycle.

Fractional abundances of isoprenoid and branched glycerol dialkyl glycerol tetraethers (GDGT) of cores SO201-2-12KL and SO201-2-114KL

It has been proposed that North Pacific sea surface temperature (SST) evolution was intimately linked to North Atlantic climate oscillations during the last glacial-interglacial transition. However, during the early deglaciation and the Last Glacial Maximum, the SST development in the subarctic northwest Pacific and the Bering Sea is poorly constrained as most existing deglacial SST records are based on alkenone paleothermometry, which is limited prior to 15 ka B.P. in the subarctic North Pacific realm. By applying the TEXL86 temperature proxy we obtain glacial-Holocene-SST records for the marginal northwest Pacific and the Western Bering Sea. Our TEXL86-based records and existing alkenone data suggest that during the past 15.5 ka, SSTs in the northwest Pacific and the Western Bering Sea closely followed millennial-scale climate fluctuations known from Greenland ice cores, indicating rapid atmospheric teleconnections with abrupt climate changes in the North Atlantic. Our SST reconstructions indicate that in the Western Bering Sea SSTs drop significantly during Heinrich Stadial 1 (HS1), similar to the known North Atlantic climate history. In contrast, progressively rising SST in the northwest Pacific is different to the North Atlantic climate development during HS1. Similarities between the northwest Pacific SST and climate records from the Gulf of Alaska point to a stronger influence of Alaskan Stream waters connecting the eastern and western basin of the North Pacific during this time. During the Holocene, dissimilar climate trends point to reduced influence of the Alaskan Stream in the northwest Pacific.

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

Sea-surface temperature and salinity reconstruction for sediment core GeoB5844-2

We present high-resolution paleoceanographic records of surface and deep water conditions within the northern Red Sea covering the last glacial maximum and termination I using alkenone paleothermometry, stable oxygen isotopes, and sediment compositional data. Paleoceanographic records in the restricted desert-surrounded northern Red Sea are strongly affected by the stepwise sea level rise and appear to record and amplify well-known millennial-scale climate events from the North Atlantic realm. During the last glacial maximum (LGM), sea surface temperatures were about 4°C cooler than the late Holocene. Pronounced coolings associated with Heinrich event 1 (~2°C below the LGM level) and the Younger Dryas imply strong atmospheric teleconnections to the North Atlantic. Owing to the restricted exchange with the Indian Ocean, Red Sea salinity is particularly sensitive to changes in global sea level. Paleosalinities exceeded 50 psu during the LGM. A pronounced freshening of the surface waters is associated with the meltwater peaks MWP1a and MWP1b owing to an increased surface-near inflow of "normal" saline water from the Indian Ocean. Vertical delta18O gradients are also increased during these phases, indicating stronger surface water stratification. The combined effect of deglacial changes in sea surface temperature and salinity on water column stratification initiated the formation of two sapropel layers, which were deposited under almost anoxic condition in a stagnant water body.

Age model and palaeoclimate records over the past 22,000 years of sediment core GeoB10053-7

The Australian-Indonesian monsoon is an important component of the climate system in the tropical Indo-Pacific region. However, its past variability, relation with northern and southern high-latitude climate and connection to the other Asian monsoon systems are poorly understood. Here we present high-resolution records of monsoon-controlled austral winter upwelling during the past 22,000 years, based on planktic foraminiferal oxygen isotopes and faunal composition in a sedimentary archive collected offshore southern Java. We show that glacial-interglacial variations in the Australian-Indonesian winter monsoon were in phase with the Indian summer monsoon system, consistent with their modern linkage through cross-equatorial surface winds. Likewise, millennial-scale variability of upwelling shares similar sign and timing with upwelling variability in the Arabian Sea. On the basis of element composition and grain-size distribution as precipitation-sensitive proxies in the same archive, we infer that (austral) summer monsoon rainfall was highest during the Bølling-Allerød period and the past 2,500 years. Our results indicate drier conditions during Heinrich Stadial 1 due to a southward shift of summer rainfall and a relatively weak Hadley cell south of the Equator. We suggest that the Australian-Indonesian summer and winter monsoon variability were closely linked to summer insolation and abrupt climate changes in the northern hemisphere.