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

Benthic foraminiferal record and age dating of sediment cores from the Gullmar Fjord at the west coast of Sweden

A high-resolution study of benthic foraminiferal assemblages was performed on a ca. eight metre long sediment core from Gullmar Fjord on the west coast of Sweden. The results of 210Pb- and AMS 14C-datings show that the record includes the two warmest climatic episodes of the last 1500 years: the Medieval Warm Period (MWP) and the recent warming of the 20th century. Both periods are known to be anomalously warm and associated with positive NAO winter indices. Benthic foraminiferal successions of both periods are compared in order to find faunal similarities and common denominators corresponding to past climate changes. During the MWP, Adercotryma glomerata, Cassidulina laevigata and Nonionella iridea dominated the assemblages. Judging from dominance of species sensitive to hypoxia and the highest faunal diversity for the last ca. 2400 years, the foraminiferal record of the MWP suggests an absence of severe low oxygen events. At the same time, faunas and d13C values both point to high primary productivity and/or increased input of terrestrial organic carbon into the fjord system during the Medieval Warm Period. Comparison of the MWP and recent warming revealed different trends in the faunal record. The thin-shelled foraminifer N. iridea was characteristic of the MWP, but became absent during the second half of the 20th century. The recent Skagerrak-Kattegat fauna was rare or absent during the MWP but established in Gullmar Fjord at the end of the Little Ice Age or in the early 1900s. Also, there are striking differences in the faunal diversity and absolute abundances of foraminifera between both periods. Changes in primary productivity, higher precipitation resulting in intensified land runoff, different oxygen regimes or even changes in the fjord's trophic status are discussed as possible causes of these faunal differences.

Estimates of annual mean surface temperatures in the late Quaternary of the tropical Atlantic

At least two modes of glacial-interglacial climate change have existed within the tropical Atlantic Ocean during the last 20,000 years. The first mode (defined by cold glacial and warm interglacial conditions) occurred symmetrically north and south of the equator and dominated the eastern boundary currents and tropical upwelling areas. This pattern suggests that mode 1 is driven by a glacial modification of surface winds in both hemispheres. The second mode of oceanic climate change, defined by temperature extremes centered on the deglaciation, was hemispherically asymmetrical, with the northern tropical Atlantic relatively cold and the southern tropical Atlantic relatively warm during deglaciation. A likely cause for this pattern of variation is a reduction of the presently northward cross-equatorial heat flux during deglaciation. No single mechanism accounts for all the data. Potential contributors to oceanic climate changes are linkage to high-latitude climates, modification of monsoonal winds by ice sheet and/or insolation changes, atmospheric CO2 and greenhouse effects, indirect effects of glacial meltwater, and variations in thermohaline overturn of the oceans.

(Table 1) Age determination of sediment cores HU89038_8 and OCE326-GGC5

Models indicate that a complete shutdown of deep and intermediate water production is a possible consequence of extreme climate conditions in the northern North Atlantic, and the high ratio of 231Pa to 230Th on Bermuda Rise is evidence that this might have happened ?17 ka during Heinrich event 1 (H1). However, new radiocarbon data from bivalves that lived at ?4.6 km on the Bermuda Rise during H1 lead to a different conclusion. The bivalve data do indeed indicate ventilation of the deep western North Atlantic was suppressed during H1 but not as much as it was during the last glacial maximum. We propose that high diatom flux to the Bermuda Rise during H1 is at least in part responsible for increased 231Pa/230Th at that time. Although we cannot say for sure why opal production was so high in a gyre center location at that time, increased leakage of silica rich waters from the Southern Ocean to the North Atlantic is one possibility.

Multivariate statistical analysis on the kaolinite/chlorite ratios from 20 South Atlantic sediment cores

Multivariate statistical analysis on the kaolinite/chlorite ratios from 20 South Atlantic sediment cores allowed for the extraction of two processes controlling the fluctuations of the kaolinite/chlorite ratio during the last 130,000 yrs, (1) the relative strength of North Atlantic Deep Water (NADW) inflow into the South Atlantic Ocean and (2) the influx of aeolian sediments from the south African continent. The NADW fluctuation can be traced in the entire deep South Atlantic while the dust signal is restricted to the vicinity of South Africa. Our data indicate that NADW formation underwent significant changes in response to glacial/interglacial climate changes with enhanced export to the Southern Hemisphere during interglacials. The most pronounced phases with Enhanced South African Dust Export (ESADE) occurred during cold Marine Isotope Stage (MIS) 5d and across the Late Glacial/Holocene transition from 16 ka to 4 ka (MIS 2 to 1). This particular pattern is attributed to the interaction of Antarctic Sea Ice extent, the position of the westerlies and the South African monsoon system.

Magnetic susceptibility dataset from the Danian Contessa Highway section, Gubbio, Italy

We studied a high-resolution multiproxy data set, including magnetic susceptibility (MS), CaCO3 content, and stable isotopes (d18O and d13C), from the stratigraphic interval covering the uppermost Maastrichtian and the lower Danian, represented by the pelagic limestones of the Scaglia Rossa Formation continuously exposed in the classic sections of the Bottaccione Gorge and the Contessa Highway near Gubbio, Italy. Variations in all the proxy series are periodic and reflect astronomically forced climate changes (i.e., Milankovitch cycles). In particular, the MS proxy reflects variations in the terrigenous dust input in this pelagic, deep-marine environment. We speculate that the dust is mainly eolian in origin and that the availability and transport of dust are influenced by variations in the vegetation cover on the Maastrichtian-Paleocene African or Asian zone, which were respectively located at tropical to subtropical latitudes to the south or far to the east of the western Tethyan Umbria-Marche Basin, and were characterized by monsoonal circulation. The dynamics of monsoonal circulation are known to be strongly dependent on precession-driven and obliquity-driven changes in insolation. We propose that a threshold mechanism in the vegetation coverage may explain eccentricity-related periodicities in the terrigenous eolian dust input. Other mechanisms, both oceanic and terrestrial, that depend on the precession amplitude modulated by eccentricity, can be evoked together with the variation of dust influx in the western Tethys to explain the detected eccentricity periodicity in the d13C record. Our interpretations of the d18O and MS records suggest a warming event ~400 k.y. prior to the Cretaceous-Paleogene (K-Pg) boundary, and a period of climatic and environmental instability in the earliest Danian. Based on these multiproxy phase relationships, we propose an astronomical tuning for these sections; this leads us to an estimate of the timing and duration of several late Maastrichtian and Danian biostratigraphic and magnetostratigraphic events.

Sedimentology, geochemistry and ostracodes from Lake Ohrid

Ancient Lake Ohrid, located in the southern Balkan Peninsula in Macedonia and Albania is characterized by a high degree of endemism and it is considered to be the oldest lake in Europe. But its exact age (between one and ten million years) and also its origin are so far not known. To unravel these uncertainties an ICDP (International Continental Scientific Drilling Program) drilling project (Scientific Collaboration On Past Speciation Conditions in Ohrid (SCOPSCO)), started in April 2013. In addition to the investigations about the age and origin, other paleolimnological studies, e.g., the reconstruction of past climate and of past lake level changes, should be performed with the drilled cores. Used proxies in such paleolimnological studies are, e.g., ostracodes because they respond sensitively to environmental changes but an accurate knowledge of their preferences and tolerances to specific environmental conditions is necessary for this purpose. So far, this knowledge about the, mostly endemic, Ohrid ostracodes was limited. Thus, within the framework of this thesis, ostracodes and a multiplicity of environmental data were collected in Lake Ohrid and its adjacent waters during four field campaigns. In a total of 47 ostracode species could be detected in the entire study area and 32 of them were found alive in Lake Ohrid. Multivariate statistic identified that water depth, salinity, conductivity, pH, and dissolved oxygen were the main determining factors for ostracode distribution in the entire study area. In Lake Ohrid, the distribution was mainly controlled by water depth, water temperature, and pH. Some ostracodes were identified as strong indicator species for important environmental variables, e.g., water temperature and water depth. A distinctive feature of Lake Ohrid was the finding of the ostracode genus Amnicythere whose species normally inhabit oligo-(meso-)haline waters and this could point to a marine origin of the lake. So far, the specialized endemic ostracodes show the highest abundances and the greatest spatial distribution in Lake Ohrid but during the sampling eight widespread species were found for the first time in the lake. They inhabited mainly the northern part of the lake, where two cities are located and industry and agriculture play a major role, and they were limited to water depths above 50 m and this could be an evidence for an increasing anthropogenic pressure because widespread ostracode species often replace endemic species. To unravel the human impact on Lake Ohrid during the last decades short sediment cores were taken and the multi-proxy study indicated that the lake productivity between the early 1920s and the late 1980s was relatively low. Diatom assemblages indicate a rising productivity in the southern part of Lake Ohrid since the mid 1970s and geochemical proxies and ostracodes point to an increasing productivity since the late 1980s in the southern and in the northern part. A slight increase in the productivity continued until 2009. Noticeable is the fact that since the early 1990s, the increasing productivity and the increasing concentrations of heavy metals correspond to a decreasing number of ostracodes in the northern part of Lake Ohrid. Perhaps, this indicates that living conditions in this lake part became less favorable for the mostly endemic ostracode species. Furthermore, the sediment samples from the cores show relatively high concentrations of arsenic, iron, and nickel. Fluctuations in ostracode assemblages from three longer sediment cores, the longest spans approximately 136 ka, taken in Lake Ohrid, correspond to fluctuations in the productivity, in the carbonate content, of the lake level, and of climate changes. Between the marine isotope stage (MIS) 6 and MIS 2 the number of ostracode valves is very low or the valves were completely absent. This corresponds to a low lake productivity, a low carbonate content, and a low lake level. At the onset of the Holocene, the number of valves increased markedly and this correlates with an increased productivity and carbonate content and a warmer climate. But during the Little Ice Age (LIA), the number of valves dropped again and species which prefer warmer waters disappeared completely. This drop corresponds also to a low productivity. After the LIA, the number of species increased again but since 1895 AD a strong and abrupt decrease is visible. A reason for this could be an increase in the heavy metal concentrations.

(Table 2) Radiocarbon content and 14C AMS ages of planktonic foraminifera from sediment cores Louis_1900 and Louis_2023

North American freshwater runoff records have been used to support the case that climate flickers were caused by shutdowns of the ocean thermohaline circulation (THC) resulting from reversals of meltwater discharges. Inconsistencies in the documentation of these meltwater switches, however, continue to fuel the debate on the cause/s of the oscillatory nature of the deglacial climate. New oxygen and carbon isotope records from the northern Gulf of Mexico depict in exceptional detail the succession of meltwater floods and pauses through the southern routing during the interval 16 to 8.9 ka (14C years BP; ka, kiloannum). The records underscore the bimodal role played by the Gulf of Mexico as a destination of meltwater discharges from the receding Laurentide Ice Sheet. The evidence indicates that the Gulf of Mexico acted as the principal source of superfloods at 13.4, 12.6, and 11.9 ka that reached the North Atlantic and contributed significantly to density stratification, disruption of ocean ventilation, and cold reversals. Gulf of Mexico lapsed into a "relief valve" position in post-Younger Dryas time, when meltwater discharges were rerouted south at 9.9, 9.7, 9.4, and 9.1 ka, thus temporarily interrupting North Atlantic-bound freshwater discharges from Lake Agassiz. The history of meltwater events in the Gulf of Mexico contradicts the model that meltwater flow via the eastern outlets into the North Atlantic disrupted the ocean THC, causing cooling, while diversions to the Gulf of Mexico via the Mississippi River enhanced THC and warming.

The impact of fluctuating light on the dinoflagellate Prorocentrum micans depends on NO3- and CO2 availability

Increasing atmospheric pCO2 and its dissolution into oceans leads to ocean acidification and warming, which reduces the thickness of upper mixing layer (UML) and upward nutrient supply from deeper layers. These events may alter the nutritional conditions and the light regime to which primary producers are exposed in the UML. In order to better understand the physiology behind the responses to the concomitant climate changes factors, we examined the impact of light fluctuation on the dinoflagellate Prorocentrum micans grown at low (1 µmol/L) or high (800 µmol/L) [NO3(-)] and at high (1000 µatm) or low (390 µatm, ambient) pCO2. The light regimes to which the algal cells were subjected were (1) constant light at a photon flux density (PFD) of either 100 (C100) or 500 (C500) µmol/m**2/s or (2) fluctuating light between 100 or 500 µmol photons/m**2/s with a frequency of either 15 (F15) or 60 (F60) min. Under continuous light, the initial portion of the light phase required the concomitant presence of high CO2 and NO3(-) concentrations for maximum growth. After exposure to light for 3h, high CO2 exerted a negative effect on growth and effective quantum yield of photosystem II (F'(v)/F'(m)). Fluctuating light ameliorated growth in the first period of illumination. In the second 3h of treatment, higher frequency (F15) of fluctuations afforded high growth rates, whereas the F60 treatment had detrimental consequences, especially when NO3(-) concentration was lower. F'(v)/F'(m) respondent differently from growth to fluctuating light: the fluorescence yield was always lower than at continuous light at 100 µmol/m**2/s, and always higher at 500 µmol/m**2/s. Our data show that the impact of atmospheric pCO2 increase on primary production of dinoflagellate depends on the availability of nitrate and the irradiance (intensity and the frequency of irradiance fluctuations) to which the cells are exposed. The impact of global change on oceanic primary producers would therefore be different in waters with different chemical and physical (mixing) properties.

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.

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.

Pollen analysis of sediment core GeoB1023-5

A high resolution marine pollen record from site GeoB1023, west of the northern Namib desert provides data on vegetation and climate change for the last 21 ka at an average resolution of 185 y. Pollen and spores are mainly delivered to the site by the Cunene river and by surface and mid-tropospheric wind systems. The main pollen source areas are located between 13°S and 21°S, which includes the northern Namib desert and semi-desert, the Angola-northern Namibian highland, and the north-western Kalahari. The pollen spectra reflect environmental changes in the region. The last glacial maximum (LGM) was characterised by colder and more arid conditions than at present, when a vegetation with temperate elements such as Asteroideae, Ericaceae, and Restionaceae grew north of 21°S. At 17.5 ka cal. B.P., an amelioration both in temperature and humidity terminated the LGM but, in the northern Kalahari, mean annual rainfall in the interval 17.5-14.4 ka cal. B.P. was probably 100-150 mm lower than at present (400-500 mm/y). The Late-glacial to early Holocene transition includes two arid periods, i.e. 14.4-12.5 and 10.9-9.3 ka cal. B.P. The last part of the former period may be correlated with the Younger Dryas. The warmest and most humid period in the Holocene occurred between 6.3 and 4.8 ka cal. B.P. During the last 2000 years, human impact, as reflected by indications of deforestation, enhanced burning and overgrazing, progressively intensified.