(Table S1) Age determination of sediment cores of the Gulf of Alaska

Oxygen isotope data from planktonic and benthic foraminifera, on a high-resolution age model (44 14C dates spanning 17,400 years), document deglacial environmental change on the southeast Alaska margin (59°33.32'N, 144°9.21'W, 682 m water depth). Surface freshening (i.e., d18O reduction of 0.8 per mil) began at 16,650 ± 170 cal years B.P. during an interval of ice proximal sedimentation, likely due to freshwater input from melting glaciers. A sharp transition to laminated hemipelagic sediments constrains retreat of regional outlet glaciers onto land circa 14,790 ± 380 cal years B.P. Abrupt warming and/or freshening of the surface ocean (i.e., additional d18O reduction of 0.9 per mil) coincides with the Bølling Interstade of northern Europe and Greenland. Cooling and/or higher salinities returned during the Allerød interval, coincident with the Antarctic Cold Reversal, and continue until 11,740 ± 200 cal years B.P., when onset of warming coincides with the end of the Younger Dryas. An abrupt 1 per mil reduction in benthic d18O at 14,250 ± 290 cal years B.P. likely reflects a decrease in bottom water salinity driven by deep mixing of glacial meltwater, a regional megaflood event, or brine formation associated with sea ice. Two laminated opal-rich intervals record discrete episodes of high productivity during the last deglaciation. These events, precisely dated here at 14,790 ± 380 to 12,990 ± 190 cal years B.P. and 11,160 ± 130 to 10,750 ± 220 cal years B.P., likely correlate to similar features observed elsewhere on the margins of the North Pacific and are coeval with episodes of rapid sea level rise. Remobilization of iron from newly inundated continental shelves may have helped to fuel these episodes of elevated primary productivity and sedimentary anoxia.

Proglacial river dataset from the Akuliarusiarsuup Kuua River northern tributary, Southwest Greenland, 2008 - 2010, version 1.0

Pressing scientific questions concerning the Greenland ice sheet's climatic sensitivity, hydrology, and contributions to current and future sea level rise require hydrological datasets to resolve. While direct observations of ice sheet meltwater losses can be obtained in terrestrial rivers draining the ice sheet and from lake levels, few such datasets exist. We present a new dataset of meltwater river discharge for the vicinity of Kangerlussuaq, Southwest Greenland. The dataset contains measurements of river stage and discharge for three sites along the Akuliarusiarsuup Kuua (Watson) River's northern tributary, with 30 minute temporal resolution between June 2008 and August 2010. Additional data of water temperature, air pressure, and lake water depth and temperature are also provided. Discharge data were measured at sites with near-ideal properties for such data collection. Regardless, high water bedload and turbulent flow introduce considerable uncertainty. These were constrained and quantified using statistical techniques, thereby providing a high quality dataset from this important site. The greatest data uncertainties are associated with streambed elevation change and measurements. Large portions of stream channels deepened according to statistical tests, but poor precision of streambed depth measurements also added uncertainty. Quality checked data are freely available for scientific use as supplementary online material.

Sedimentology and stable carbon and oxygen isotope record of the Central Arctic

Stable oxygen and carbon isotope and sedimentological-paleontological investigations supported by accelerator mass spectrometry 14C datings were carried out on cores from north of 85°N in the eastern central Arctic Ocean. Significant changes in accumulation rates, provenance of ice-rafted debris (IRD), and planktic productivity over the past 80,000 years are documented. During peak glacials, i.e., oxygen isotope stages 4 and 2, the Arctic Ocean was covered by sea ice with decreased seasonal variation, limiting planktic productivity and bulk sedimentation rates. In early stage 3 and during Termination I, major deglaciations of the circum-Arctic regions caused lowered salinities and poor oxygenation of central Arctic surface waters. A meltwater spike and an associated IRD peak dated to ~14-12 14C ka can be traced over the southern Eurasian Basin of the Arctic Ocean. This event was associated with the early and rapid deglaciation of the marine-based Barents Sea Ice Sheet. A separate Termination Ib meltwater event is most conspicuous in the central Arctic and is associated with characteristic dolomitic carbonate IRD. This lithology suggests an origin of glacial ice from northern Canada and northern Greenland where lower Paleozoic platform carbonates crop extensively out.

Ice rafted debris, planktic foraminifera and reconstructed sea surface temperatures from site GIK23414

Planktic foraminiferal census data, faunal sea surface temperatures (SSTs) and oxygen isotopic and lithic records from a site in the northeast Atlantic were analyzed to study the interglacial dynamics of Marine Isotope Stage (MIS) 11, a period thought to closely resemble the Holocene on the basis of orbital forcing. Interglacial conditions during MIS 11 persisted for approximately 26 ka. After the main deglacial meltwater processes ceased, a 10- to 12-ka-long transitional period marked by significant water mass circulation changes occurred before surface waters finally reached their thermal maximum. This SST peak occurred between 400 and 397 ka, inferred from the abundance of the most thermophilic foraminiferal species and was coincident with lowest sea level according to benthic isotope values. The ensuing stepwise SST decrease characterizes the overall climate deterioration preceding the increase in global ice volume by ab. 3 ka. This cooling trend was followed by a more pronounced cold event that began at 388 ka, and that terminated in the recurrence of icebergs at the site around 382 ka. Because the water mass configuration of early MIS 11 evolved quite differently from that of the early Holocene, there is little evidence that MIS 11 can serve as an appropriate analogue for a future Holocene climate, despite the similarity in some orbital parameters.

(Appendix 1) Coral analysis and isostatic rebound effects from different Holes of IODP Expedition 310

Past sea-level records provide invaluable information about the response of ice sheets to climate forcing. Some such records suggest that the last deglaciation was punctuated by a dramatic period of sea-level rise, of about 20 metres, in less than 500 years. Controversy about the amplitude and timing of this meltwater pulse (MWP-1A) has, however, led to uncertainty about the source of the melt water and its temporal and causal relationships with the abrupt climate changes of the deglaciation. Here we show that MWP-1A started no earlier than 14,650 years ago and ended before 14,310 years ago, making it coeval with the Bølling warming. Our results, based on corals drilled offshore from Tahiti during Integrated Ocean Drilling Project Expedition 310, reveal that the increase in sea level at Tahiti was between 12 and 22 metres, with a most probable value between 14 and 18 metres, establishing a significant meltwater contribution from the Southern Hemisphere. This implies that the rate of eustatic sea-level rise exceeded 40 millimetres per year during MWP-1A.

Sedimentology and benthic foraminiferal assemblages of the late Younger Dryas and early Holocene in the Skagerrak

A high-resolution study of palaeoenvironmental changes through the late Younger Dryas and early Holocene in the Skagerrak, the eastern North Atlantic, is based on multi-proxy analyses of core MD99-2286 combined with palaeo-water depth modelling for the area. The late Younger Dryas was characterized by a cold ice-distal benthic foraminiferal fauna. After the transition to the Preboreal (c. 11 650 cal. a BP) this fauna was replaced by a Cassidulina neoteretis dominated fauna, indicating the influence of chilled Atlantic Water at the sea floor. Persisting relatively cold bottom-water conditions until c. 10 300 cal. a BP are presumably a result of an outflow of glacial meltwater from the Baltic area across south-central Sweden, which develops a strong stratification of the water column at MD99-2286. A short-term peak in the C/N ratio at c. 10 200 cal. a BP is suggested to indicate input of terrestrial material, which may represent the drainage of an ice-dammed lake in southern Norway, the Glomma event. After the last drainage route across south-central Sweden closed, c. 10 300 cal. a BP, the meltwater influence diminished, and the Skagerrak resembled a fjord with stable inflow of waters from the North Atlantic through the Norwegian Channel and a gradual increase in boreal species. Full interglacial conditions were established at the sea floor from c. 9250 cal. a BP. Subsequent warm stable conditions were interrupted by a short-term cooling around 8300-8200 cal. a BP, representing the 8.2 ka event.

Clay mineralogy of sediment cores from the Arctic Ocean

This study focuses on sedimentological investigations of sediment cores recovered during the international Arctic'91, expeditions with the German research ice breaker RV "Polarstern" to the European sector of the Arctic Ocean. Here, we deduce the last glacial/interglacial changes in transport mechanism and sedimentation from the clay mineral group smectite. We choose the smectites as an example of how sediment mineralogy can be linked with particular source regions (the Kara and Laptev seas), distinct transport mechanism (sea ice and surface currents) and sedimentation processes. Smectite contents in Arctic sediments discussed for two time slices, including the Last Glacial Maximum (LGM), and the last deglaciation (Termination I), reveal the highest variability subsequent to the retreat of the Eurasian ice sheets. Our results show that smectite anomalies in the Eurasian Basin are associated with distinct meltwater pulses and occurred around 13.5-13.0 14C ka B.P. Compelling evidence is provided that these anomalies are deduced from sea-ice entrained sediments from the eastern Kara Sea that entered the Arctic Ocean after ice-sheet break-up and eventually flooding of the Kara Sea. We propose that smectite anomalies in sediments of the eastern Arctic Ocean can be utilized to identify deglacial events and to help decipher configurations of the Eurasian ice sheets. The identification of smectite maxima along the modern sea-ice edge in the Eurasian Basin further indicates biologically enhanced sedimentation from melting sea ice allowing the reconstruction of seasonally open water in the region. Hence, considering the poor preservation conditions of primary paleoceanographic proxies in the Arctic Ocean, the clay mineral contents, particularly the smectite group, may be one alternative tool for paleoclimatic reconstruction in the Eurasian Basin.

Sea surface temperature calculation based on alkenones in sediments off northwest Africa

The molecular stratigraphy of Biogeochemical Oceanic Flux Study core 31K (19°N, 20°10'W) and Ocean Drilling Program Hole 658C (20°45'N, 18°35'W) has been studied for C37 alkenone abundances over the past 80 ka at high resolution (~circa 200-500 years). The derived Uk 37' sea surface temperature record for both cores shows a range of temperatures from about 18°C during the last glacial to 21.5°C during the early Holocene. Both records also reveal changes in sea surface temperature as much as 2°-4°C over a few hundred years, which correlate well with similar abrupt climatic changes observed in cores from elsewhere in the NE Atlantic, associated with 'Heinrich events'. Our results indicate that meltwater produced by these ice-rafting events was transmitted southward by the Canary Current, where it had considerable impact on sea surface temperatures in the subtropical eastern Atlantic.

Age determination and stable isotope composition of sediment cores from the Arctic Ocean

On the basis of 52 sediment cores, analyzed and dated at high resolution, the paleoceanography and climate of the Last Glacial Maximum (LGM) were reconstructed in detail for the Fram Strait and the eastern and central Arctic Ocean. Sediment composition and stable isotope data suggest three distinct paleoenvironments: (1) a productive region in the eastern to central Fram Strait and along the northern Barents Sea continental margin characterized by Atlantic Water advection, frequent open water conditions, and occasional local meltwater supply and iceberg calving from the Barents Sea Ice Sheet; (2) an intermediate region in the southwestern Eurasian Basin (up to 84-85°N) and the western Fram Strait characterized by subsurface Atlantic Water advection and recirculation, a moderately high planktic productivity, and a perennial ice cover that breaks up only occasionally; and (3) a central Arctic region (north of 85°N in the Eurasian Basin) characterized by a low-salinity surface water layer and a thick ice cover that strongly reduces bioproduction and bulk sedimentation rates. Although the total inflow of Atlantic Water into the Arctic Ocean may have been reduced during the LGM, its impact on ice coverage and halocline structure in the Fram Strait and southwestern Eurasian Basin was strong.

Calcareous plankton response to orbital and millennial-scale climate changes from ODP Leg 161-975

The paleoenvironmental conditions through MIS 15-9 at the Mediterranean Ocean Drilling Program (ODP) Site 975 were interpreted by high resolution study of calcareous plankton assemblages compared with available d18O and d13C records and high resolution paleoclimate proxies from the Atlantic Ocean. Sea Surface Temperatures (SSTs) have been estimated from planktonic foraminiferal assemblages using the artificial neural networks method. Calcareous plankton varied dominantly on a glacial-interglacial scale as testified by the SST record, foraminiferal diversity, total coccolith abundance and changes in warm-water calcareous nannofossil taxa. A general increase in foraminiferal diversity and of total coccolith abundance is observed during interglacials. Warmest SSTs are reached during MIS 11, while MIS 12 and MIS 10 represent the coldest intervals of the studied record. During MIS 12, one of the most extreme glacials of the last million years, occurrence of Globorotalia inflata and of neogloboquadrinids indicates a shoaling of the interface between Atlantic inflowing and Mediterranean outflowing waters. Among calcareous nannofossils the distribution of Gephyrocapsa margereli-G. muellerae > 4 µm also supports a reduced Atlantic-Mediterranean exchange during MIS 12. Superimposed on glacial-interglacial variability, six short-terms coolings are recognized during MIS 12 and 10, which appear comparable in their distribution and amplitude to the Heinrich - type events documented in the Atlantic Ocean in the same interval. During these H-type events, N. pachyderma (s) and G. margereli-G. muellerae > 4 µm increase as a response to the enhanced inflow of cold Atlantic water into the Mediterranean via the Strait of Gibraltar. Mediterranean surface water hydrography appears to have been most severely affected at Termination V during the H-type event Ht4, possibly as a response to a large volume of Atlantic meltwater inflow via the Strait of Gibraltar and/or to freshwater/terrigenous input deriving from local mountain glaciers. Three additional SST coolings are recorded through MIS 14-16, but these are not well correlated with Heinrich - type events documented in the Atlantic Ocean in the same interval; during these cooling episodes only the subpolar Turborotalita quinqueloba increases. These results highlight the sensitive response of the Mediterranean basin to millennial-scale climate variations related to Northern Hemisphere ice-sheet instability and support the hypothesis that the tight connection between high latitude climate dynamics and Mediterranean sea surface water features can be traced through the Middle Pleistocene.

(Table 1) Age determination of Arctic Ocean sediment cores

Accelerator mass spectrometer 14C dated stable isotope data from Neogloboquadrina pachyerma in cores raised from the Mendeleyev Ridge and slope provide evidence for significant influx of meltwater to the western Arctic Ocean during the early part of marine oxygen isotope stage 1 (OIS 1) and during several intervals within OIS 3. The strongest OIS 3 meltwater event occurred before ca. 45 ka (conventional radiocarbon age) and was probably related to the deglaciation at the beginning of OIS 3. Major meltwater input to the western Arctic Ocean during the last deglaciation coincides closely with the maximum rate of global sea-level rise as determined from the Barbados sea-level record, demonstrating a strong link between the global record and changes in the central Arctic Ocean. OIS 2, which includes the last glacial maximum, is very condensed or absent in the cores. Abundance and d13C values for N. pachyderma in the middle part of OIS 3 are similar to modern values, indicating high productivity and seasonal ice-free areas along the Arctic margin at that time. These records indicate that the Arctic Ocean was a source of heat and moisture to the northern polar atmosphere during parts of OIS 3.

(Table 1) Age determination of sediment cores south of the Faeroe Islands

A study was made of three cores from the Faeroe-Shetland gateway, based on planktonic foraminifera, oxygen isotopes, accelerator mass spectrometry 14C dates, magnetic susceptibility, and counts of ice rafted debris (IRD). The data, covering the period 30-10 ka, show that during the Last Glacial Maximum the Arctic Front occupied a position close to the Faeroes, allowing a persisting inflow of Atlantic surface water into the Faeroe-Shetland Channel. The oceanographic environment during deposition of two IRD layers is influenced by Atlantic surface water masses during the lower IRD layer, with transport of icebergs from N-NW. Polar surface water conditions prevailed only during deposition of the upper IRD layer. There is no indication of surface meltwater influence in the region during the deglaciation, but there is a persistent influence of Atlantic surface water masses in the region. Thus we conclude that during almost the entire period (30-10 ka) the Faeroe-Shetland Channel was a gateway for transport of Atlantic surface water toward the Norwegian Sea.

(Table 1) Age determination and isotopic signature of last glacial Patagonian caly/silt

Ice cores provide a record of changes in dust flux to Antarctica, which is thought to reflect changes in atmospheric circulation and environmental conditions in dust source areas (Forster et al., 2007; Diekmann et al. 2000, doi:10.1016/S0031-0182(00)00138-3; Winckler et al., 2008, doi:10.1126/science.1150595; Reader et al., 1999, doi:10.1029/1999JD900033; Mahowald et al., 1999, doi:10.1029/1999JD900084; Petit et al., 1999, doi:10.1038/20859; 1990, doi:10.1038/343056a0 Delmonte et al., 2009, doi:10.1029/2008GL033382; Lambert et al., 2008, doi:10.1038/nature06763). Isotopic tracers suggest that South America is the dominant source of the dust (Grousset et al., 1992, doi:10.1016/0012-821X(92)90177-W; Basile et al., 1997, doi:10.1016/S0012-821X(96)00255-5; Gaiero et al., 2007, doi:10.1016/j.chemgeo.2006.11.003), but it is unclear what led to the variable deposition of dust at concentrations 20-50 times higher than present in glacial-aged ice (Petit et al., 1990, doi:10.1038/343056a0; Lambert et al., 2008, doi:10.1038/nature06763). Here we characterize the age and composition of Patagonian glacial outwash sediments, to assess the relationship between the Antarctic dust record from Dome C (refs Lambert et al., 2008, doi:10.1038/nature06763; Wolff et al., 2006, doi:10.1038/nature04614) and Patagonian glacial fluctuations (Sugden et al., 2005; McCulloch et al., 2005, doi:10.1111/j.0435-3676.2005.00260.x; Kaplan et al., 2008, doi:10.1016/j.quascirev.2007.09.013) for the past 80,000 years. We show that dust peaks in Antarctica coincide with periods in Patagonia when rivers of glacial meltwater deposited sediment directly onto easily mobilized outwash plains. No dust peaks were noted when the glaciers instead terminated directly into pro-glacial lakes. We thus propose that the variable sediment supply resulting from Patagonian glacial fluctuations may have acted as an on/off switch for Antarctic dust deposition. At the last glacial termination, Patagonian glaciers quickly retreated into lakes, which may help explain why the deglacial decline in Antarctic dust concentrations preceded the main phase of warming, sea-level rise and reduction in Southern Hemisphere sea-ice extent (Wolff et al., 2006, doi:10.1038/nature04614).

Analysis of minerals from Arctic Ocean sediments

The Arctic Ocean System is a key player regarding the climatic changes of Earth. Its highly sensitive ice Cover, the exchange of surface and deep water masses with the global ocean and the coupling with the atmosphere interact directly with global climatic changes. The output of cold, polar water and sea ice influences the production of deep water in the North Atlantic and controls the global ocean circulation ("the conveyor belt"). The Arctic Ocean is surrounded by the large Northern Hemisphere ice sheets which not only affect the sedimentation in the Arctic Ocean but also are supposed to induce the Course of glacials and interglacials. Terrigenous sediment delivered from the ice sheets by icebergs and meltwater as well as through sea ice are major components of Arctic Ocean sediments. Hence, the terrigenous content of Arctic Ocean sediments is an outstanding archive to investigate changes in the paleoenvironment. Glazigenic sediments of the Canadian Arctic Archipelago and surface samples of the Arctic Ocean and the Siberian shelf regions were investigated by means of x-ray diffraction of the bulk fraction. The source regions of distinct mineral compositions were to be deciphered. Regarding the complex circumpolar geology stable christalline shield rocks, active and ancient fold belts including magmatic and metamorphic rocks, sedimentary rocks and wide periglacial lowlands with permafrost provide a complete range of possible mineral combinations. Non- glaciated shelf regions mix the local input from a possible point source of a particular mineral combination with the whole shelf material and function as a sampler of the entire region draining to the shelf. To take this into account, a literature research was performed. Descriptions of outcropping lithologies and Arctic Ocean sediments were scanned for their mineral association. The analyses of glazigenic and shelf sediments yielded a close relationship between their mineral composition and the adjacent source region. The most striking difference between the circumpolar source regions is the extensive outcrop of carbonate rocks in the vicinity of the Canadian Arctic Archipelago and in N Greenland while siliciclastic sediments dominate the Siberian shelves. In the Siberian shelf region the eastern Kara Sea and the western Laptev Sea form a destinct region defined by high smectite, (clino-) pyroxene and plagioclase input. The source of this signal are the extensive outcrops of the Siberian trap basalt in the Putorana Plateau which is drained by the tributaries of the Yenissei and Khatanga. The eastern Laptev Sea and the East Siberian Sea can also be treated as one source region containing a feldspar, quartz, illite, mica, and chlorite asscciation combined with the trace minerals hornblende and epidote. Franz Josef Land provides a mineral composition rich in quartz and kaolinite. The diverse rock suite of the Svalbard archipelago distributes specific mineral compositions of highly metamorphic christalline rocks, dolomite-rich carbonate rocks and sedimentary rocks with a higher diagenetic potential manifested in stable newly built diagenetic minerals and high organic maturity. To reconstruct the last 30,000 years as an example of the transition between glacial and interglacial conditions a profile of sediment cores, recovered during the RV Polarstern" expedition ARK-VIIIl3 (ARCTIC '91), and additional sediment cores around Svalbard were investigated. Besides the mineralogy of different grain size fractions several additional sedimentological and organo-geochemical Parameterswere used. A detailed stratigraphic framework was achieved. By exploiting this data set changes in the mineral composition of the Eurasian Basin sediments can be related to climatic changes. Certain mineral compositions can even be associated with particular transport processes, e.g. the smectitel pyroxene association with sea ice transport from the eastern Kara Sea and the western Laptev Sea. Hence, it is possible to decipher the complex interplay between the influx of warm Atlantic waters into the Southwest of the Eurasian Basin, the waxing and waning of the Svalbard1Barents- Sea- and Kara-Sea-Ice-Sheets, the flooding of the Siberian shelf regions and the surface and deep water circulation. Until now the Arctic Ocean was assumed to be a rather stable System during the last 30,000 years which only switched from a completely ice covered situation during the glacial to seasonally Open waters during the interglacial. But this work using mineral assemblages of sediment cores in the vicinity of Svalbard revealed fast changes in the inflow of warm Atlantic water with the Westspitsbergen Current (< 1000 years), short periods of advances and retreats of the marine based Eurasian ice sheets (1000-3000 years), and short melting phases (400 years?). Deglaciation of the marine-based Eurasian and the land-based north American and Greenland ice sheets are not simultaneous. This thesis postulates that the Kara Sea Ice Sheet released an early meltwater signal prior to 15,000 14C years leading the Barents Sea Ice Sheet while the western land-based ice sheets are following later than 13,500 14C years. The northern Eurasian Basin records the shift between iceberg and sea-ice material derived from the Canadian Arctic Archipelago and N-Greenland and material transported by sea-ice and surface currents from the Siberian shelf region. The phasing of the deglaciation becomes very obvious using the dolomite and quartd phyllosilicate record. It is also supposed that the flooding of the Laptev Sea during the Holocene is manifested in a stepwise increase of sediment input at the Lomonosov Ridge between the Eurasian and Amerasian Basin. Depending on the strength of meltwater pulses from the adjacent ice sheets the Transpolar Drift can probably be relocated. These movements are traceable by the distribution of indicator minerals. Based on the outcome of this work the feasibility of bulk mineral determination can be qualified as excellent tool for paleoenvironmental reconstructions in the Arctic Ocean. The easy preparation and objective determination of bulk mineralogy provided by the QUAX software bears the potential to use this analyses as basic measuring method preceding more time consuming and highly specialised mineralogical investigations (e.g. clay mineralogy, heavy mineral determination).