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

Dinocyst-derived sea-surface characteristics, palynomorph abundances, ice-rafted debris, and planktonic foraminifera abundances and d18O of the 35-41 ka BP interval of core MD99-2285

The last glacial millennial climatic events (i.e. Dansgaard-Oeschger and Heinrich events) constitute outstanding case studies of coupled atmosphere-ocean-cryosphere interactions. Here, we investigate the evolution of sea-surface and subsurface conditions, in terms of temperature, salinity and sea ice cover, at very high-resolution (mean resolution between 55 and 155 years depending on proxies) during the 35-41 ka cal BP interval covering three Dansgaard-Oeschger cycles and including Heinrich event 4, in a new unpublished marine record, i.e. the MD99-2285 core (62.69°N; -3.57s°E). We use a large panel of complementary tools, which notably includes dinocyst-derived sea-ice cover duration quantifications. The high temporal resolution and multiproxy approach of this work allows us to identify the sequence of processes and to assess ocean-cryosphere interactions occurring during these periodic ice-sheet collapse events. Our results evidence a paradoxical hydrological scheme where (i) Greenland interstadials are marked by a homogeneous and cold upper water column, with intensive winter sea ice formation and summer sea ice melting, and (ii) Greenland and Heinrich stadials are characterized by a very warm and low saline surface layer with iceberg calving and reduced sea ice formation, separated by a strong halocline from a less warm and saltier subsurface layer. Our work also suggests that this stadial surface/subsurface warming started before massive iceberg release, in relation with warm Atlantic water advection. These findings thus support the theory that upper ocean warming might have triggered European ice-sheet destabilization. Besides, previous paleoceanographic studies conducted along the Atlantic inflow pathways close to the edge of European ice-sheets suggest that such a feature might have occurred in this whole area. Nonetheless, additional high resolution paleoreconstructions are required to confirm such a regional scheme.

Results of modelled Arctic river discharge trends since 7 ka BP

The Arctic hydrological cycle throughout the Holocene is analyzed based on the results of transient simulations with the coupled atmosphere-ocean circulation model ECHO-G. The results suggest a ~ 2 % increase of mid-Holocene to preindustrial Arctic river discharges for the Eurasian continent. However, rivers of the North America Arctic realm show a moderate runoff decline of approximately 4 to 5 % for the same period. The total river discharge into the Arctic Ocean has remained at an approximately constant preindustrial level since the mid Holocene. The positive discharge trend within Eurasia is caused by a more rapid decrease in local net evaporation compared to a smaller decline in advected moisture and hence precipitation. This effect is neither recognized within the North American Arctic domain nor in the far eastern part of the Eurasian Arctic realm. A detailed comparison of these model findings with a variety of proxy studies is conducted. The collected proxy records show trends of continental surface temperatures and precipitation rates that are consistent with the simulations. A continuation of the transient Holocene runs for the 19th and 20th century with increased greenhouse gases indicates an increase of the total river influx into the Arctic Ocean of up to 7.6 %. The Eurasian river discharges increase by 7.5 %, the North American discharges by up to 8.4 %. The most rapid increases have been detected since the beginning of the 20th century. These results are corroborated by the observed rising of Arctic river discharges during the last century which is attributed to anthropogenic warming. The acceleration of the Arctic hydrological cycle in the 20th century is without precedence in the Holocene.

Planktonic foraminifera oxygen isotope and trace metal data ODP Hole 172-1056A, 46-54 ka BP

Salinity increase in the subtropical gyre system may have pre-conditioned the North Atlantic Ocean for a rapid return to stronger overturning circulation and high-latitude warming following meltwater events during the Last Glacial period. Here we investigate the Gulf Stream - subtropical gyre system properties over Dansgaard-Oeschger (DO) cycles 14 to 12, including Heinrich ice-rafting event 5. During the Holocene and Last Glacial Maximum a positive gradient in surface dwelling planktonic foraminifera d18O (Globigerinoides ruber) can be observed between the Gulf Stream and subtropical gyre, due to decreasing temperature, increasing salinity, and a change from summer to year-round occurrence of G. ruber. We assess whether this gradient was a common feature during stadial-interstadial climate oscillations of Marine Isotope Stage 3, by comparing existing G. ruber d18O from ODP Site 1060 (subtropical gyre location) and new data from ODP Site 1056 (Gulf Stream location) between 54 and 46 ka. Our results suggest that this gradient was largely absent during the period studied. During the major warm DO interstadials 14 and 12 we infer a more zonal and wider Gulf Stream, influencing both ODP Sites 1056 and 1060. A Gulf Stream presence during these major interstadials is also suggested by the large vertical d18O gradient between shallow dwelling planktonic foraminifera species, especially G. ruber, and the deep dwelling species Globorotalia inflata at site 1056, which we associate with strong summer stratification and Gulf Stream presence. A major reduction in this vertical d18O gradient from 51 ka until the end of Heinrich event 5 at 48.5 ka suggests site 1056 was situated within the subtropical gyre in this mainly cold period, from which we infer a migration of the Gulf Stream to a position nearer to the continental shelf, indicative of a narrower Gulf Stream with possibly reduced transport.

Reference datasets and pollen record from sediment core YAK, Yakumo, SW Hokkaido

This study reconstructs middle and late Holocene vegetation and climate dynamics in the Oshima Peninsula, SW Hokkaido, using the published method of biome reconstruction and modern analogue technique applied to the Yakumo pollen record (42°17'03''N, 140°15'34''E) spanning the last 5500 years. Two previously published matrices assigning Japanese plant/pollen taxa to the major vegetation types (biomes) are tested using a newly compiled dataset of 78 surface pollen spectra from Hokkaido. With both matrices showing strengths and weaknesses in reconstructing cool mixed and temperate deciduous forests of Hokkaido, the results suggest the necessity to consider the whole list of identified terrestrial pollen taxa for generating robust vegetation reconstructions for northern Japan. Applied to the fossil pollen data, both biome-reconstruction approaches demonstrate consistently that oak-dominated cool mixed forest spread in the study region between 5.5 and 3.6 cal ka BP and was subsequently replaced by beech-dominated temperate deciduous forest. The pollen-based climate reconstruction suggests this change in the vegetation composition was caused by a shift from cooler and drier than present climate to warmer and wetter, similar to modern conditions about 3.6 cal ka BP. Comparing the pollen-based reconstruction results with the published marine records from the NW Pacific, the reconstructed vegetation and climate dynamics can be satisfactorily explained by the greater role played by the warm Tsushima Current in the Sea of Japan and in the Tsugaru Strait during the middle and late Holocene. An increase in sea surface temperatures west and south of the study site would favour air temperature rise and moisture uptake and cause an increase in precipitation and snow accumulation in the western part of Hokkaido during the late Holocene.

Lake level reconstructions at Laguna Potrok Aike (southern Patagonia, Argentina) for the last 49 ka

Four seismic surveys and a stratigraphic record from southernmost Patagonia (Argentina) based on 51 AMS-14C dates obtained in the framework of ICDP expedition 5022 "Potrok Aike Maar Lake Sediment Archive Drilling Project" (PASADO) provide a database to compare the 106 m composite profile from the lake centre with piston cores from the littoral and outcrops in the catchment area. Based on event correlation using distinct volcanic ash layers with unique geochemical composition and optically stimulated luminescence (OSL) dates on feldspars, sediment records are firmly linked. This approach allows to match the sediment record with water levels during the past ca. 49 ka providing evidence for lake level variations. Reconstructed lake levels were 20 m higher than today during the last Glacial until the early Holocene. With the migration of the Southern Hemispheric Westerlies over this site the lake level dropped ca. 55 m for a period of two millennia. Thereupon the water balance was more positive again causing a stepwise rise of the lake level until the maximum was reached during the Little Ice Age with a subsequent lowering since the 20th century. We suggest that the mid- to late-Holocene lake level variation is caused by intensity changes of the Southern Hemispheric Westerlies.

Controls of sediment dynamics on the Galician slope reconstructed from three sediment cores (GeoB11035-1, GeoB130206-1, GeoB13071-1)

Controls of sediment dynamics at the Galician continental slope (NW Iberia) during the past 30 ka were reconstructed from three new gravity cores (GeoB11035-1, 130206-1, 13071-1) based on sedimentological (e.g. sortable silt, IRD), micropalaeontological (e.g. coccoliths), geochemical (AMS 14C, XRF) and geophysical (e.g. magnetic susceptibility) diagnostics. The data are consistent with existing regional knowledge that, during marine isotope stages 3-1, variations in detrital input, marine productivity and sea level were the essential drivers of sediment availability on the slope, whereas deep-water current velocities controlled sediment deposition: (1) the period prior to 30 cal ka BP is characterized by minor but systematic variations in various proxies which can be associated with D-O cycles; (2) between 30 and 18 cal ka BP, high detrital input and steady slope-parallel currents led to constant sedimentation; (3) from the LGM until 10 cal ka BP, the shelf-transgressive sea-level rise increased the detrital particle flux; sedimentation was influenced by significantly enhanced deep-water circulation during the Bølling/Allerød, and subsequent slowing during the Younger Dryas; (4) an abrupt and lasting change to hemipelagic sedimentation at ca. 10 cal ka BP was probably due to Holocene warming and decelerated transgression; (5) after 5 cal ka BP, additional input of detrital material to the slope is plausibly linked to the evolution of fine-grained depocentres on the Galician shelf, this being the first report of this close shelf-slope sedimentary linkage off NW Iberia. Furthermore, there is novel evidence of the nowadays strong outer shelf Iberian Poleward Current becoming established at about 15.5 cal ka BP. The data also demonstrate that small-scale morphologic features and local pathways of sediment export from the neighbouring shelf play an important role for sediment distribution on the NW Iberian slope, including a hitherto unknown sediment conduit off the Ría de Arousa. By implication, the impact of local morphology on along- and down-slope sediment dynamics is more complex than commonly considered, and deserves future attention.

Sediment composition in NE Greenland Fjords, Shelf, and Slope

In order to document changes in Holocene glacier extent and activity in NE Greenland (~73° N) we study marine sediment records that extend from the fjords (PS2631 and PS2640), across the shelf (PS2623 and PS2641), to the Greenland Sea (JM07-174GC). The primary bedrock geology of the source areas is the Caledonian sediment outcrop, including Devonian red beds, plus early Neoproterozoic gneisses and early Tertiary volcanics. We examine the variations in colour (CIE*), grain size, and bulk mineralogy (from X-ray diffraction of the <2 mm sediment fraction). Fjord core PS2640 in Sofia Sund, with a marked red hue, is distinct in grain size, colour and mineralogy from the other fjord and shelf cores. Five distinct grain-size modes are distinguished of which only one is associated with a coarse ice-rafting signal - this mode is rare in the mid- and late Holocene. A sediment unmixing program (SedUnMixMC) is used to characterize down-core changes in sediment composition based on the upper late Holocene sediments from cores PS2640 (Sofia Sund), PS2631 (Kaiser Franz Joseph Fjord) and PS2623 (south of Shannon Is), and surface samples from the Kara Sea (as an indicator of transport from the Russian Arctic shelves). Major changes in mineral composition are noted in all cores with possible coeval shifts centred c. 2.5, 4.5 and 7.5 cal. ka BP (±0.5 ka) but are rarely linked with changes in the grain-size spectra. Coarse IRD (>2 mm) and IRD-grain-size spectra are rare in the last 9-10 cal. ka BP and, in contrast with areas farther south (~68° N), there is no distinct IRD signal at the onset of neoglaciation. Our paper demonstrates the importance of the quantitative analysis of sediment properties in clarifying source to sink changes in glacial marine environments.