(Table S1) Age determination of sediment core RAPiD-21-03K

A decadal resolution time series of sea surface temperature (SST) spanning the last two millennia is reconstructed by combining a proxy record from a new sediment sequence with previously published data from core MD99-2275, north of Iceland. The alkenone based SST reconstruction is validated with historic observational data and compared to a new similar temporal resolution reconstruction obtained from sediment core RAPiD21-3K, in the subpolar North Atlantic. The two SST paleorecords show consistent multidecadal scale coolings throughout the interval and similar expressions during the contrasted climatic periods 'colloquially known' as the Medieval Climatic Anomaly (MCA) and Little Ice Age (LIA). In order to further understand the temporal and spatial SST variations and investigate the influence of natural forcings on the observed SST changes during the last millennium, we compare our time series to simulations using the Institut Pierre-Simon Laplace IPSLCM4-v2 climate model. This comparison highlights the potential importance of volcanism as a natural forcing driving coherent abrupt cooling events captured in the subpolar North Atlantic records.

(Table T1) Age determination of seven ODP Leg 201 sites

As part of a wider paleoclimate and paleoceanographic study of Holocene-upper Pleistocene laminated sediments from the eastern equatorial Pacific and Peru continental margin, we completed 32 accelerator mass spectrometry (AMS) 14C dates from cores recovered during Ocean Drilling Program (ODP) Leg 201. Sample preparation and measurement were carried out at the ANTARES AMS facility, Australian Nuclear Science and Technology Organisation (ANSTO), in Sydney, Australia (Lawson et al., 2000, doi:10.1016/S0168-583X(00)00276-7; Fink et al., 2004, doi:10.1016/j.nimb.2004.04.025). Although the sediments are predominantly diatomaceous oozes (D'Hondt, Jørgensen, Miller, et al., 2003, doi:10.2973/odp.proc.ir.201.2003), they contain sufficient inorganic (e.g., foraminifer tests and nannofossil plates) and organic (Meister et al., 2005, doi:10.2973/odp.proc.sr.201.105.2005) carbon to allow 14C dating. These dates permitted us to reconstruct a history of sediment accumulation over the past 20 k.y., particularly on the Peru continental margin. In this report we present 14C AMS dates and other pertinent data from cores from Sites 1227, 1228, and 1229 collected during Leg 201 at the Peru continental margin.

Age determination and stable isotope records of sediment core M78/1_235-1

Paleoenvironmental studies and climate models demonstrate that fluvial runoff and moisture availability in the Caribbean hinterland react very sensitively to climatic variations. Late Pleistocene and Holocene climate records document pronounced dry and wet periods over tropical South America mainly caused by shifts of the Intertropical Convergence Zone (ITCZ). However, forcing mechanisms for changes in the ITCZ position remain controversial. Here we present high-resolution foraminiferal Ba/Ca and d18Oseawater records from a core located within the Orinoco River outflow documenting abrupt hydrological changes in the Orinoco catchment area during the deglacial and Holocene. Our data, obtained from the surface-dwelling foraminifera Globigerinoides ruber (pink), show an abrupt increase in Ba/Ca ratios in the early Holocene, starting ~600 yr after the end of the Younger Dryas (YD) cold interval at ca. 10.8 ka and suggesting a massive reorganization of moisture sources in northern South America. In contrast, the salinity dependent d18Oseawater from the same samples shows a gradual decrease starting at the end of the YD. The offset of our Ba/Ca peak excludes meltwater release in conjunction with the northern Andean glacier retreat well before the end of the YD as a forcing mechanism. We suggest that the Ba/Ca record documents an abrupt increase in Ba-rich waters of a northern Andean source caused by the insolation-driven shift of the ITCZ and/or enhanced monsoon activity.

Age determination of sediment core HZM, lake Holzmar (Tables 1, 2)

AMS radiocarbon ages have been determined on terrestrial macrofossils selected from the annually laminated sediments of lake Holzmaar (Germany). The radiocarbon chronology of this lake covers the last 12.6 ka. Comparison of the radiocarbon dated varve chronology with tree ring data shows that an additional 878 years have to be added to the varve chronology. The corrected 14C varve chronology of Holzmaar reaches back to ca. 13.8 ka cal. BP and compares favourably with the results from Soppensee (Switzerland) (Hajdas et al., 1993, doi:10.1007/BF00209748). The corrected ages for the onset and the end of the Younger Dryas biozone are 11,940 cal. BP and 11,490 cal. BP, respectively. The ash layer of the Laacher See volcanic eruption is dated at 12,201 ± 224 cal. BP and the Ulmener Tephra layer is dated at 10,904 cal. BP.

(Table 2) Age determination of Southern Ocean sediments

The silicic acid leakage hypothesis (SALH) predicts that during glacial periods excess silicic acid was transported from the Southern Ocean to lower latitudes, which favored diatom production over coccolithophorid production and caused a drawdown of atmospheric CO2. Downcore records of 230Th-normalized opal (biogenic silica) fluxes from 31 cores in the Pacific sector of the Southern Ocean were used to compare diatom productivity during the last glacial period to that of the Holocene and to examine the evidence for increased glacial Si export to the tropics. Average glacial opal fluxes south of the modern Antarctic Polar Front (APF) were less than during the Holocene, while average glacial opal fluxes north of the APF were greater than during the Holocene. However, the magnitude of the increase north of the APF was not enough to offset decreased fluxes to the south, resulting in a decrease in opal burial in the Pacific sector of the Southern Ocean during the last glacial period, equivalent to approximately 15 Gt opal/ka1. This is consistent with the work of Chase et al. (2003, doi:10.1016/S0967-0645(02)00595-7), and satisfies the primary requirement of the SALH, assuming that the upwelled supply of Si was approximately equivalent during the Holocene and the glacial period. However, previous results from the equatorial oceans are inconsistent with the other predictions of the SALH, namely that either the Corg:CaCO3 ratio or the rate of opal burial should have increased during glacial periods. We compare the magnitudes of changes in the Southern Ocean and the tropics and suggest that Si escaping the glacial Southern Ocean must have had an alternate destination, possibly the continental margins. There is currently insufficient data to test this hypothesis, but the existence of this sink and its potential impact on glacial pCO2 remain interesting topics for future study.

(Table 1) Age determination of DSDP Hole 94-607 and sediment cores V30-101 and V23-81

Oxygen isotope measurements in Greenland ice demonstrate that a series of rapid warm-cold oscillations -called Dansgaard-Oeschger events- punctuated the last glaciation (Dansgard et al., 1993, doi:10.1038/364218a0). Here we present records of sea surface temperature from North Atlantic sediments spanning the past 90 kyr which contain a series of rapid temperature oscillations closely matching those in the ice-core record, confirming predictions that the ocean must bear the imprint of the Dansgaard-Oeschger events (Broecker et al., 1988, doi:10.1016/0033-5894(88)90082-8; 1990, doi:10.1029/PA005i004p00469). Moreover, we show that between 20 and 80 kyr ago, the shifts in ocean-atmosphere temperature are bundled into cooling cycles, lasting on average 10 to 15 kyr, with asymmetrical saw-tooth shapes. Each cycle culminated in an enormous discharge of icebergs into the North Atlantic (a 'Hein-rich event' (Bond et al., 1992, doi:10.1038/360245a0; Broecker et al., 1992, doi:10.1007/BF00193540), followed by an abrupt shift to a warmer climate. These cycles document a previously unrecognized link between ice sheet behaviour and ocean-atmosphere temperature changes. An important question that remains to be resolved is whether the cycles are driven by external factors, such as orbital forcing, or by inter-nal ice-sheet dynamics.

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.

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

Vertical profiles of environmental parameters measured on discrete water samples collected with Niskin bottles at station TARA_209 during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) was a global survey of ocean ecosystems aboard the Sailing Vessel Tara. It carried out extensive measurements of evironmental conditions and collected plankton (viruses, bacteria, protists and metazoans) for later analysis using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set includes properties of seawater, particulate matter and dissolved matter that were measured from discrete water samples collected with Niskin bottles during the 2009-2013 Tara Oceans expedition. Properties include pigment concentrations from HPLC analysis (10 depths per vertical profile, 25 pigments per depth), the carbonate system (Surface and 400m; pH (total scale), CO2, pCO2, fCO2, HCO3, CO3, Total alkalinity, Total carbon, OmegaAragonite, OmegaCalcite, and dosage Flags), nutrients (10 depths per vertical profile; NO2, PO4, N02/NO3, SI, quality Flags), DOC, CDOM, and dissolved oxygen isotopes. The Service National d'Analyse des Paramètres Océaniques du CO2, at the Université Pierre et Marie Curie, determined CT and AT potentiometrically. More than 200 vertical profiles of these properties were made across the world ocean. DOC, CDOM and dissolved oxygen isotopes are available only for the Arctic Ocean and Arctic Seas (2013).

Vertical profiles of environmental parameters measured on discrete water samples collected with Niskin bottles at station TARA_131 during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) was a global survey of ocean ecosystems aboard the Sailing Vessel Tara. It carried out extensive measurements of evironmental conditions and collected plankton (viruses, bacteria, protists and metazoans) for later analysis using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set includes properties of seawater, particulate matter and dissolved matter that were measured from discrete water samples collected with Niskin bottles during the 2009-2013 Tara Oceans expedition. Properties include pigment concentrations from HPLC analysis (10 depths per vertical profile, 25 pigments per depth), the carbonate system (Surface and 400m; pH (total scale), CO2, pCO2, fCO2, HCO3, CO3, Total alkalinity, Total carbon, OmegaAragonite, OmegaCalcite, and dosage Flags), nutrients (10 depths per vertical profile; NO2, PO4, N02/NO3, SI, quality Flags), DOC, CDOM, and dissolved oxygen isotopes. The Service National d'Analyse des Paramètres Océaniques du CO2, at the Université Pierre et Marie Curie, determined CT and AT potentiometrically. More than 200 vertical profiles of these properties were made across the world ocean. DOC, CDOM and dissolved oxygen isotopes are available only for the Arctic Ocean and Arctic Seas (2013).

Vertical profiles of environmental parameters measured on discrete water samples collected with Niskin bottles at station TARA_132 during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) was a global survey of ocean ecosystems aboard the Sailing Vessel Tara. It carried out extensive measurements of evironmental conditions and collected plankton (viruses, bacteria, protists and metazoans) for later analysis using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set includes properties of seawater, particulate matter and dissolved matter that were measured from discrete water samples collected with Niskin bottles during the 2009-2013 Tara Oceans expedition. Properties include pigment concentrations from HPLC analysis (10 depths per vertical profile, 25 pigments per depth), the carbonate system (Surface and 400m; pH (total scale), CO2, pCO2, fCO2, HCO3, CO3, Total alkalinity, Total carbon, OmegaAragonite, OmegaCalcite, and dosage Flags), nutrients (10 depths per vertical profile; NO2, PO4, N02/NO3, SI, quality Flags), DOC, CDOM, and dissolved oxygen isotopes. The Service National d'Analyse des Paramètres Océaniques du CO2, at the Université Pierre et Marie Curie, determined CT and AT potentiometrically. More than 200 vertical profiles of these properties were made across the world ocean. DOC, CDOM and dissolved oxygen isotopes are available only for the Arctic Ocean and Arctic Seas (2013).

Vertical profiles of environmental parameters measured on discrete water samples collected with Niskin bottles at station TARA_139 during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) was a global survey of ocean ecosystems aboard the Sailing Vessel Tara. It carried out extensive measurements of evironmental conditions and collected plankton (viruses, bacteria, protists and metazoans) for later analysis using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set includes properties of seawater, particulate matter and dissolved matter that were measured from discrete water samples collected with Niskin bottles during the 2009-2013 Tara Oceans expedition. Properties include pigment concentrations from HPLC analysis (10 depths per vertical profile, 25 pigments per depth), the carbonate system (Surface and 400m; pH (total scale), CO2, pCO2, fCO2, HCO3, CO3, Total alkalinity, Total carbon, OmegaAragonite, OmegaCalcite, and dosage Flags), nutrients (10 depths per vertical profile; NO2, PO4, N02/NO3, SI, quality Flags), DOC, CDOM, and dissolved oxygen isotopes. The Service National d'Analyse des Paramètres Océaniques du CO2, at the Université Pierre et Marie Curie, determined CT and AT potentiometrically. More than 200 vertical profiles of these properties were made across the world ocean. DOC, CDOM and dissolved oxygen isotopes are available only for the Arctic Ocean and Arctic Seas (2013).

Vertical profiles of environmental parameters measured on discrete water samples collected with Niskin bottles at station TARA_124 during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) was a global survey of ocean ecosystems aboard the Sailing Vessel Tara. It carried out extensive measurements of evironmental conditions and collected plankton (viruses, bacteria, protists and metazoans) for later analysis using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set includes properties of seawater, particulate matter and dissolved matter that were measured from discrete water samples collected with Niskin bottles during the 2009-2013 Tara Oceans expedition. Properties include pigment concentrations from HPLC analysis (10 depths per vertical profile, 25 pigments per depth), the carbonate system (Surface and 400m; pH (total scale), CO2, pCO2, fCO2, HCO3, CO3, Total alkalinity, Total carbon, OmegaAragonite, OmegaCalcite, and dosage Flags), nutrients (10 depths per vertical profile; NO2, PO4, N02/NO3, SI, quality Flags), DOC, CDOM, and dissolved oxygen isotopes. The Service National d'Analyse des Paramètres Océaniques du CO2, at the Université Pierre et Marie Curie, determined CT and AT potentiometrically. More than 200 vertical profiles of these properties were made across the world ocean. DOC, CDOM and dissolved oxygen isotopes are available only for the Arctic Ocean and Arctic Seas (2013).