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

(Table 1) Age control and temporal resolution of paleoceanographic records of the North Atlantic

Short-term changes in sea surface conditions controlling the thermohaline circulation in the northern North Atlantic are expected to be especially efficient in perturbing global climate stability. Here we assess past variability of sea surface temperature (SST) in the northeast Atlantic and Norwegian Sea during Marine Isotope Stage (MIS) 2 and, in particular, during the Last Glacial Maximum (LGM). Five high-resolution SST records were established on a meridional transect (53°N-72°N) to trace centennial-scale oscillations in SST and sea-ice cover. We used three independent computational techniques (SIMMAX modern analogue technique, Artificial Neural Networks (ANN), and Revised Analog Method (RAM)) to reconstruct SST from planktonic foraminifer census counts. SIMMAX and ANN reproduced short-term SST oscillations of similar magnitude and absolute levels, while RAM, owing to a restrictive analog selection, appears less suitable for reconstructing "cold end" SST. The SIMMAX and ANN SST reconstructions support the existence of a weak paleo-Norwegian Current during Dansgaard-Oeschger (DO) interstadials number 4, 3, 2, and 1. During the LGM, two warm incursions of 7°C water to occurred in the northern North Atlantic but ended north of the Iceland Faroe Ridge. A rough numerical estimate shows that the near-surface poleward heat transfer from 53° across the Iceland-Faroe Ridge up to to 72° N dropped to less than 60% of the modern value during DO interstadials and to almost zero during DO stadials. Summer sea ice was generally confined to the area north of 70°N and only rarely expanded southward along the margins of continental ice sheets. Internal LGM variability of North Atlantic (>40°N) SST in the GLAMAP 2000 compilation (Sarnthein et al., 2003, doi:10.1029/2002PA000771; Pflaumann et al., 2003, doi:10.1029/2002PA000774) indicates maximum instability in the glacial subpolar gyre and at the Iberian Margin, while in the Nordic Seas, SST was continuously low.

Marine04 marine radiocarbon age calibration, 0-26 cal kyr BP

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0–26 cal kyr BP (Before Present, 0 cal BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0–10.5 cal kyr BP. Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific 14C reservoir age information to provide a single global marine mixed-layer calibration from 10.5–26.0 cal kyr BP. A substantial enhancement relative to IntCal98 is the introduction of a random walk model, which takes into account the uncertainty in both the calendar age and the 14C age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

Dust flux records from the Subarctic North Pacific

We present a new record of eolian dust flux to the western Subarctic North Pacific (SNP) covering the past 27000 years based on a core from the Detroit Seamount. Comparing the SNP dust record to the NGRIP ice core record shows significant differences in the amplitude of dust changes to the two regions during the last deglaciation, while the timing of abrupt changes is synchronous. If dust deposition in the SNP faithfully records its mobilization in East Asian source regions, then the difference in the relative amplitude must reflect climate-related changes in atmospheric dust transport to Greenland. Based on the synchronicity in the timing of dust changes in the SNP and Greenland, we tie abrupt deglacial transitions in the 230Th-normalized 4He flux record to corresponding transitions in the well-dated NGRIP dust flux record to provide a new chronostratigraphic technique for marine sediments from the SNP. Results from this technique are complemented by radiocarbon dating, which allows us to independently constrain radiocarbon paleoreservoir ages. We find paleoreservoir ages of 745 ± 140 yr at 11653 yr BP, 680 ± 228 yr at 14630 yr BP and 790 ± 498 yr at 23290 yr BP. Our reconstructed paleoreservoir ages are consistent with modern surface water reservoir ages in the western SNP. Good temporal synchronicity between eolian dust records from the Subantarctic Atlantic and equatorial Pacific and the ice core record from Antarctica supports the reliability of the proposed dust tuning method to be used more widely in other global ocean regions.

Bulk geochemical analysis and planktonic foraminifera of sediment core GeoB7112-5

A multivariable approach utilising bulk sediment, planktonic Foraminifera and siliceous phytoplankton has been used to reconstruct rapid variations in palaeoproductivity in the Peru-Chile Current System off northern Chile for the past 19000 cal. yr. During the early deglaciation (19000-16000 cal. yr BP), our data point to strongest upwelling intensity and highest productivity of the past 19 000 cal. yr. The late deglaciation (16000-13000 cal. yr BP) is characterised by a major change in the oceanographic setting, warmer water masses and weaker upwelling at the study site. Lowest productivity and weakest upwelling intensity are observed from the early to the middle Holocene (13000-4000 cal. yr BP), and the beginning of the late Holocene (<4000 cal. yr BP) is marked by increasing productivity, mainly driven by silicate-producing organisms. Changes in the productivity and upwelling intensity in our record may have resulted from a large-scale compression and/or displacement of the South Pacific subtropical gyre during more productive periods, in line with a northward extension of the Antarctic Circumpolar Current and increased advection of Antarctic water masses with the Peru-Chile Current. The corresponding increase in hemispheric thermal gradient and wind stress induced stronger upwelling. During the periods of lower productivity, this scenario probably reversed.

Analyses of lake La Thuile sediment core

Lake La Thuile, in the Northern French Prealps (874 m a.s.l.), provides an 18 m long sedimentary sequence spanning the entire Lateglacial/Holocene period. The high resolution multi-proxy (sedimentological, palynological, geochemical) analysis of the uppermost 6.2 meters reveals the Holocene dynamics of erosion in the catchment in response to landscape modifications. The mountain belt is at relevant altitude to study past human activities and the watershed is sufficiently disconnected from large valleys to capture a local sedimentary signal. From 12,000 to 10,000 cal. BP (10 to 8 ka cal. BC), the onset of hardwood species triggered a drop in erosion following the Lateglacial/Holocene transition. From 10,000 to 4500 cal. BP (8 to 2.5 ka cal. BC), the forest became denser and favored slope stabilization while erosion processes were very weak. A first erosive phase was initiated at ca . 4500 cal. BP without evidence of human presence in the catchment. Then, the forest declined at approximately 3000 cal. BP, suggesting the first human influence on the landscape. Two other erosive phases are related to anthropic activities: approximately 2500 cal. BP (550 cal. BC) during the Roman period and after 1600 cal. BP (350 cal. AD) with a substantial accentuation in the Middle Ages. In contrast, the lower erosion produced during the Little Ice Age, when climate deteriorations are generally considered to result in an increased erosion signal in this region, suggests that anthropic activities dominated the erosive processes and completely masked the natural effects of climate on erosion in the late Holocene.