(Appendix A) Relative abundance of Emiliania huxleyi in sediment cores of the North Atlantic

The coccolithophore species Emiliania huxleyi is characterized by a wide range of sizes, which can be easily distinguished in the light microscope. In this study we have quantified the abundance of large (coccoliths > 4 µm in maximum length) E. huxleyi specimens during the last 25 kyr in sedimentary records from eleven cores and drill sites in the NE Atlantic and W Mediterranean Sea, to prove its usefulness in the reconstruction of water mass dynamics and biostratigraphic potential. During the Last Glacial Maximum this large form, a cold-water indicator, was common in the NE Atlantic and Mediterranean, and its regional variation in abundance indicates a displacement of the climatic zones southwards in agreement with the development of ice sheets and sea ice in the Northern Hemisphere during this period. On the other hand, the gradient between northern and southern surface water masses in the Subtropical Gyre appears to have been more pronounced than at present, while the Portugal and Canary Currents were more intense. In the western Mediterranean basin temperatures were cooler than in the adjacent Atlantic, provoking a quasi-endemism of these specimens until the end of Heinrich Event 1. This may have been due to a restriction in the communication between the Atlantic and Mediterranean through the Strait of Gibraltar, the arrival of cold surface water and the amplification of cooling after the development of ice sheets in the Northern Hemisphere. During the deglaciation, large E. huxleyi specimens decreased in abundance at medium and low latitudes, but were still numerous close to the Subarctic region during the Holocene. In transitional waters this decrease to present day abundances occurred after Termination Ib. The abrupt change in abundance of this large E. huxleyi form is proposed as a new biostratigraphic event to characterize the Holocene in mid- to low-latitude water masses in the North Atlantic, although this horizon seems to be diachronous by 5 kyr from tropical to subarctic regions, in agreement with the gradual onset of warm conditions.

Micropaleontological and geochemical proxy studies of sediment core GeoB7926-2

[1] The low-latitude upwelling regime off the Mauritanian coast in the subtropical NE Atlantic accounts for a significant part of global export production. Although productivity variations in coastal upwelling areas are usually attributed to changes in wind stress and upwelling intensity, productivity dynamics off Mauritania are less straightforward because of the complex atmospheric and hydrographic setting. Here we integrate micropaleontological (diatoms) and geochemical (bulk biogenic sediment components, X-ray fluorescence, and alkenones) proxies to examine on submillennial-to-millennial changes in diatom production that occurred off Mauritania, NW Africa, for the last 25 ka. During the Last Glacial Maximum (LGM, 19.0-23.0 ka B.P.), moderate silicate content of upwelled waters coupled with weakened NE trade winds determined moderate diatom productivity. No significant cooling is observed during the LGM, suggesting that our alkenone-based SST reconstruction represents a local, upwelling-related signal rather than a global insolation related one. Extraordinary increases in diatom and opal concentrations during Heinrich event 1 (H1, 15.5-18.0 ka B.P.) and the Younger Dryas (YD, 13.5-11.5 ka B.P.) are attributed to enhanced upwelling of silica-rich waters and an enlarged upwelling filament, due to more intense NE trade winds. The synchronous increase of CaCO3 and K intensity and the decreased opal and diatoms values mark the occurrence of the Bølling/Allerød (BA, 13.5-15.5 ka B.P.) due to weakened eolian input and more humid conditions on land. Although the high export of diatoms is inextricably linked to upwelling intensity off Mauritania, variability in the nutrient content of the thermocline also plays a decisive role.

Paleotemperature reconstructions on sediment core GIK161603-3

Different proxies for sea surface temperature (SST) often exhibit divergent trends for deglacial warming in tropical regions, hampering our understanding of the phase relationship between tropical SSTs and continental ice volume at glacial terminations. To reconcile divergent SST trends, we report reconstructions of two commonly used paleothermometers (the foraminifera G. ruber Mg/Ca and the alkenone unsaturation index) from a marine sediment core collected in the southwestern tropical Indian Ocean encompassing the last 37,000 years. Our results show that SSTs derived from the alkenone unsaturation index (UK'37) are consistently warmer than those derived from Mg/Ca by ~2-3°C except for the Heinrich Event 1. In addition, the initial timing for the deglacial warming of alkenone SST started at ~15.6 ka, which lags behind that of Mg/Ca temperatures by 2.5 kyr. We argue that the discrepancy between the two SST proxies reflects seasonal differences between summer and winter rather than post-depositional processes or sedimentary biases. The UK'37 SST record clearly mimics the deglacial SST trend recorded in the North Atlantic region for the earlier part of the termination, indicating the early deglacial warming trend attributed to local summer temperatures was likely mediated by changes in the Atlantic Meridional Overturning Circulation at the onset of the deglaciation, In contrast, the glacial to interglacial SST pattern recorded by G. ruber Mg/Ca probably reflects cold season SSTs. This indicates that the cold season SSTs was likely mediated by climate changes in the southern hemisphere, as it closely tracks the Antarctic timing of deglaciation. Therefore our study reveals that the tropical southwestern Indian Ocean seasonal SST was closely linked to climate changes occurring in both hemispheres. The austral summer and winter recorded by each proxy is further supported with seasonal SST trends modeled by AOGCMs for our core site. Our interpretation that the alkenone and Mg/Ca SSTs are seasonally biased may also explain similar proxy mismatches observed in other tropical regions at the onset of the last termination.

Nd isotopes of sediments, foraminifera and fish debris of the western North Atlantic

Understanding changes in ocean circulation during the last deglaciation is crucial to unraveling the dynamics of glacial-interglacial and millennial climate shifts. We used neodymium isotope measurements on postdepositional iron-manganese oxide coatings precipitated on planktonic foraminifera to reconstruct changes in the bottom water source of the deep western North Atlantic at the Bermuda Rise. Comparison of our deep water source record with overturning strength proxies shows that both the deep water mass source and the overturning rate shifted rapidly and synchronously during the last deglacial transition. In contrast, any freshwater perturbation caused by Heinrich event 1 could have only affected shallow overturning. These findings show how changes in upper-ocean overturning associated with millennial-scale events differ from those associated with whole-ocean deglacial climate events.

Age determination and geochemical composition of sediment core DA04-31P

There is limited knowledge pertaining to the history of the Greenland Ice Sheet (GIS) during the last glacial-interglacial transition as it retreated from the continental margins to an inland position. Here we use multiproxy data, including ice-rafted debris (IRD); planktonic isotopes; alkenone temperatures; and tephra geochemistry from the northern Labrador Sea, off southwest Greenland, to investigate the deglacial response of the GIS and evaluate its implications for the North Atlantic deglacial development. The results imply that the southern GIS retreated in three successive stages: (1) early deglaciation of the East Greenland margins, by tephra-rich IRD that embrace Heinrich Event 1; (2) progressive retreat during Allerød culminating in major meltwater releases (d18O depletion of 1.2 per mil) at the Allerød-Younger Dryas transition (12.8-13.0 kyr B.P.); and (3) a final stage of glacial recession during the early Holocene (~9-11 kyr B.P.). Rather than indicating local temperatures of ambient surface water, the alkenones likely were transported to the core site by the Irminger Current. We attribute the timing of GIS retreat to the incursion of warm intermediate waters along the base of grounded glaciers and below floating ice shelves on the continental margin. The results lend support to the view that GIS meltwater presented a forcing factor for the Younger Dryas cooling.

Glacial-interglacial sediment record from the North Atlantic

Studies from the subtropical western and eastern Atlantic Ocean, using the 231Pa/230Th ratio as a kinematic proxy for deep water circulation, provided compelling evidence for a strong link between climate and the rate of meridional overturning circulation (MOC) over the last deglaciation. In this study, we present a compilation of existing and new sedimentary 231Pa/230Th records from North Atlantic cores between 1710 and 4550 m water depth. Comparing sedimentary 231Pa/230Th from different depths provides new insights into the evolution of the geometry and rate of deep water formation in the North Atlantic during the last 20,000 years. The 231Pa/230Th ratio measured in upper Holocene sediments indicates slow water renewal above ?2500 m and rapid flushing below, consistent with our understanding of modern circulation. In contrast, during the Last Glacial Maximum (LGM), Glacial North Atlantic Intermediate Water (GNAIW) drove a rapid overturning circulation to a depth of at least ?3000 m depth. Below ~4000 m, water renewal was much slower than today. At the onset of Heinrich event 1, transport by the overturning circulation declined at all depths. GNAIW shoaled above 3000 m and significantly weakened but did not totally shut down. During the Bølling-Allerød (BA) that followed, water renewal rates further decreased above 2000 m but increased below. Our results suggest for the first time that ocean circulation during that period was quite distinct from the modern circulation mode, with a comparatively higher renewal rate above 3000 m and a lower renewal rate below in a pattern similar to the LGM but less accentuated. MOC during the Younger Dryas appears very similar to BA down to 2000 m and slightly slower below.

XRF-scanned elemental concentrations from sediment cores off Peru

We present a high-resolution marine record of sediment input from the Guayas River, Ecuador, that reflects changes in precipitation along western equatorial South America during the last 18ka. We use log (Ti/Ca) derived from X-ray Fluorescence (XRF) to document terrigenous input from riverine runoff that integrates rainfall from the Guayas River catchment. We find that rainfall-induced riverine runoff has increased during the Holocene and decreased during the last deglaciation. Superimposed on those long-term trends, we find that rainfall was probably slightly increased during the Younger Dryas, while the Heinrich event 1 was marked by an extreme load of terrigenous input, probably reflecting one of the wettest period over the time interval studied. When we compare our results to other Deglacial to Holocene rainfall records located across the tropical South American continent, different modes of variability become apparent. The records of rainfall variability imply that changes in the hydrological cycle at orbital and sub-orbital timescales were different from western to eastern South America. Orbital forcing caused an antiphase behavior in rainfall trends between eastern and western equatorial South America. In contrast, millennial-scale rainfall changes, remotely connected to the North Atlantic climate variability, led to homogenously wetter conditions over eastern and western equatorial South America during North Atlantic cold spells. These results may provide helpful diagnostics for testing the regional rainfall sensitivity in climate models and help to refine rainfall projections in South America for the next century.

Late glacial and deglacial age model and biomarker analyses of sediment core MSM05/5_712-2 from the western continental slope of Svalbard

The transition from last glacial to deglacial and subsequently to modern interglacial climate conditions was accompanied by abrupt shifts in the palaeoceanographic setting in the subpolar North Atlantic. Knowledge about the role that sea ice coverage played during these rapid climate reversals is limited since most marine sediment cores from the higher latitudes provide only a coarse temporal resolution and often poorly preserved microfossils. Here we present a highly resolved reconstruction of sea ice conditions that characterised the eastern Fram Strait - a key area for water mass exchange between the Arctic Ocean and the North Atlantic - for the past 30 ka BP. This reconstruction is based on the distribution of the sea ice biomarker IP25 and phytoplankton derived biomarkers in a sediment core from the continental slope of western Svalbard. During the late glacial (30 ka to 19 ka BP), recurrent advances and retreats of sea ice characterised the study area and point to a hitherto less considered oceanic (and/or atmospheric) variability. A long-lasting perennial sea ice coverage in eastern Fram Strait persisted only at the very end of the Last Glacial Maximum (i.e. from 19.2 to 17.6 ka BP) and was abruptly reduced at the onset of Heinrich Event 1 - coincident with or possibly even inducing the collapse of the Atlantic Meridional Overturning Circulation (AMOC). Further maximum sea ice conditions prevailed during the Younger Dryas cooling event and support the assumption of an AMOC reduction due to increased formation and export of Arctic sea ice through Fram Strait. A significant retreat of sea ice and sea surface warming are observed for the Early Holocene.

Stepwise consolidation of glacigenic sediments of the Prydz Bay

A series of samples from the five sites drilled across the continental shelf and upper slope in Prydz Bay during ODP Leg 119 were consolidation tested in an oedometer. Preconsolidation stresses increase downcore at Sites 739 and 742 in a stepwise manner, and the steps are interpreted to represent periods of increased action of grounded glaciers covering the entire shelf. By the use of theoretical ice sheet surface profiles giving the range of possible ice thicknesses, sediment loading and subsequent erosion seem to be the most important factor for increasing the overconsolidation ratios, and a total glacial erosion exceeding 1 km is possible. Four separate steps in consolidation, here termed "load events" have been identified. The lowermost load event, 1, is correlated to the onset of glaciations reaching the shelf edge and an early period of extensive glaciations, starting in early Oligocene or possibly earlier. Glacial activity related to the buildup of ice in West Antarctica in the late Miocene is tentatively correlated to load event 2. Event 3 is the trace of relatively extensive glacial erosion probably in the Pliocene, whereas the upper step in preconsolidation stress, load event 4, results from the last glaciation reaching the shelf edge, possibly during the late Weichselian. Correlations to other data related to Antarctic glacial history are, however, hampered by the poor age control of the cored diamictites. Consolidation tests may provide a tool for finding the position for hiatuses and unconformities formed subglacially and obscured by subglacial reworking.

Sea-surface temperature reconstruction for sediment core GeoB1023-5

Sediment core GeoB 1023-5 from the eastern South Atlantic was investigated at high temporal resolution for variations of sea-surface temperature (SST) during the past 22 kyr, using the alkenone (UK'37) method. SSTs increased by 3.5°C from about 18°C during the Last Ice Age (21±2 cal kyr BP) to about 21.5°C at 14.5 cal kyr BP. This warming trend associated with the deglaciation phase was followed by a cooling event with lowest SSTs near 20°C, persisting for about 1000 years between 13 and 12 cal kyr BP. The SSTs then continued to increase to about 22.5°C at the Holocene climatic optimum at 7 cal kyr BP, and decreased again during the Late Holocene to a core-top value of 19.8°C that is comparable to modern annual mean SST values. When compared with alkenone SST records from the eastern North Atlantic, our SST record indicates continuous warming throughout the deglaciation phase in the Benguela Current, while its northern counterpart, the Canary Current, experienced prominent cooling during 'Heinrich Event 1' (H1). On the other hand, for the time period corresponding to the 'Younger Dryas' (YD) cooling event, the Benguela SST record exhibits a cold-temperature interval that corresponds to that observed in the eastern North Atlantic SST records. This observation suggests that interhemispheric climate response in Atlantic eastern boundary current systems was different with respect to the two abrupt climate events associated with Termination I. For the H1, the eastern South Atlantic SST record strongly supports the hypothesis that an 'anti-phase' thermal behavior in South Atlantic surface waters was forced by the slowdown of the North Atlantic Deep Water formation during cold spells in the North Atlantic. In contrast, the abrupt cooling in the eastern South Atlantic coincident with the YD period was probably induced by more vigorous global atmospheric circulation, enhancing the upwelling intensity in both eastern boundary current systems. This atmospheric control may have overridden any effect caused by changes in thermohaline circulation on the South Atlantic SSTs during the YD, which leads to the assumption that the thermohaline circulation was already much closer to its interglacial mode during the YD than during the H1.

Uranium, Thorium and Protactinium isotopes of ODP Site 172-1063

Ocean circulation may have undergone reductions and reinvigorations in the past closely tied to regional climate changes. Measurements of 231Pa/230Th ratios in a sediment core from the Bermuda Rise have been interpreted as evidence that the Atlantic Meridional Overturning Circulation (AMOC) was weakened or completely eliminated during a period of catastrophic iceberg discharges (Heinrich-Event 1, H1). Here we present new data from the Bermuda Rise that show further 231Pa/230Th peaks during Heinrich-2 (H2) and Heinrich-3 (H3). Additionally, a tight correlation between diatom abundances (biogenic silica) and 231Pa/230Th is discovered in this core. Our results redirect the interpretation of 231Pa/230Th from the Bermuda Rise as a proxy for ocean circulation towards a proxy that reacts highly sensitive to changes of particle composition and water mass properties.