Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

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

Isotope and Mg/Ca measurements from Cores ME0005A-43JC and 202-1242 of tropical Pacific

Multiproxy geologic records of d18O and Mg/Ca in fossil foraminifera from sediments under the Eastern Pacific Warm Pool (EPWP) region west of Central America document variations in upper ocean temperature, pycnocline strength, and salinity (i.e., net precipitation) over the past 30 kyr. Although evident in the paleotemperature record, there is no glacial-interglacial difference in paleosalinity, suggesting that tropical hydrologic changes do not respond passively to high-latitude ice sheets and oceans. Millennial variations in paleosalinity with amplitudes as high as 4 practical salinity units occur with a dominant period of 3-5 ky during the glacial/deglacial interval and 1.0-1.5 ky during the Holocene. The amplitude of the EPWP paleosalinity changes greatly exceeds that of published Caribbean and western tropical Pacific paleosalinity records. EPWP paleosalinity changes correspond to millennial-scale climate changes in the surface and deep Atlantic and the high northern latitudes, with generally higher (lower) paleosalinity during cold (warm) events. In addition to Intertropical Convergence Zone (ITCZ) dynamics, which play an important role in tropical hydrologic variability, changes in Atlantic-Pacific moisture transport, which is closely linked to ITCZ dynamics, may also contribute to hydrologic variations in the EPWP. Calculations of interbasin salinity average and interbasin salinity contrast between the EPWP and the Caribbean help differentiate long-term changes in mean ITCZ position and Atlantic-Pacific moisture transport, respectively.

(Table 1) Age determination on monospecies planktonic foraminifera of sediment core AN25-934

Reconstruction of regional climate and the Okhotsk Sea (OS) environment for the Last Glacial Maximum (LGM), deglaciation and Holocene were performed on the basis of high-resolution records of ice rafted debris (IRD), CaCO3, opal, total organic carbon (TOC), biogenic Ba (Ba_bio) and redox sensitive element (Mn, Mo) content, and diatom and pollen results of four cores that form a north-southern transect. Age models of the studied cores were earlier established by AMS 14C data, oxygen - isotope chronostratigraphy and tephrochronology. According to received results, since 25 ka the regional climate and OS environmental conditions have changed synchronously with LGM condition, cold Heinrich event 1, Bølling -Allerød (BA) warming, Younger Dryas (YD) cooling and Pre-Boreal (PB) warming recorded in the Greenland ice core, North Atlantic sediment, and China cave stalagmites. Calculation of IRD MAR in sediment of north-south transect cores indicate an increase of sea ice formation several times in the glacial OS as compared to the Late Holocene. Accompanying ice formation, increased brine rejection and the larger potential density of surface water at the north shelf due to a drop of glacial East Asia summer monsoon precipitation and Amur River run off, led to strong enhancement of the role of the OS in glacial North Pacific Intermediate Water (NPIW) formation. The remarkable increase in OS productivity during BA and PB warming was probably related with significant reorganisation of the North Pacific deep water ventilation and nutrient input into the NPIW and OS Intermediate Water (OSIW). Seven Holocene OS millennial cold events based on the elevated values of the detrended IRD stack record over the IRD broad trend in the sediments of the studied cores have occurred synchronously with cold events recorded in the North Atlantic, Greenland ice cores and China cave stalagmites after 9 ka. Diatom production in the OS were mostly controlled by sea ice cover changes and surface water stratification induced by sea-ice melting; therefore significant opal accumulation in sediments of this basin begin from 4-6 ka ago simultaneously with a remarkable decrease of sea ice cover.

Relative contents of foraminifera species in bottom sediments of northern and southern seas

Correlation of paleoceanographic events in several key regions of the World Ocean: North Atlantic, Antarctic, West Arctic Seas, North Pacific and tropical Indo-Pacific has been carried out for the last 135 ka based on micropaleontological, stable isotope, geochronological (AMS-14C) and other data. It has been shown that the global thermohaline circulation controls remote climatic teleconnections on millennial-scale and partly on centennial-scale, while short-term climate changes are mainly transferred by the atmosphere. The basic information is given about the recent thermohaline circulation and stages of its development during Neogene.

Downcore rock-magnetic properties of ICDP Core 5045-1, Lake Ohrid

The bulk magnetic mineral record from Lake Ohrid, spanning the past 637 kyr, reflects large-scale shifts in hydrological conditions, and, superimposed, a strong signal of environmental conditions on glacial-interglacial and millennial timescales. A shift in the formation of early diagenetic ferrimagnetic iron sulfides to siderites is observed around 320 ka. This change is probably associated with variable availability of sulfide in the pore water. We propose that sulfate concentrations were significantly higher before ~320 ka, due to either a higher sulfate flux or lower dilution of lake sulfate due to a smaller water volume. Diagenetic iron minerals appear more abundant during glacials, which are generally characterized by higher Fe/Ca ratios in the sediments. While in the lower part of the core the ferrimagnetic sulfide signal overprints the primary detrital magnetic signal, the upper part of the core is dominated by variable proportions of high- to low-coercivity iron oxides. Glacial sediments are characterized by high concentration of high-coercivity magnetic minerals (hematite, goethite), which relate to enhanced erosion of soils that had formed during preceding interglacials. Superimposed on the glacial-interglacial behavior are millennial-scale oscillations in the magnetic mineral composition that parallel variations in summer insolation. Like the processes on glacial-interglacial timescales, low summer insolation and a retreat in vegetation resulted in enhanced erosion of soil material. Our study highlights that rock-magnetic studies, in concert with geochemical and sedimentological investigations, provide a multi-level contribution to environmental reconstructions, since the magnetic properties can mirror both environmental conditions on land and intra-lake processes.

Temperature variabilty in the southeastern Black-Sea during 64-20ka BP

The Eurasian inland propagation of temperature anomalies during glacial millennial-scale climate variability is poorly understood but this knowledge is crucial to understanding hemisphere-wide atmospheric teleconnection patterns and climate mechanisms. Based on biomarkers and geochemical paleothermometers, a pronounced continental temperature variability between 64,000 and 20,000 years ago, coinciding with the Greenland Dansgaard-Oeschger cycles, was determined in a well-dated sediment record from the formerly enclosed Black Sea. Cooling during Heinrich events was not stronger than during other stadials in the Black Sea. This is corroborated by modeling results showing that regular Dansgaard-Oeschger cycles penetrated deeper into the Eurasian continent than Heinrich events. The pattern of coastal ice-rafted detritus suggests a strong dependence on the climate background state, with significantly milder winters during periods of reduced Eurasian ice sheets and an intensified meridional atmospheric circulation.

Pollen profile and age determination for sediment core M72/5_619-1 (22-GC3)

Sediments from the Black Sea, a region historically dominated by forests and steppe landscapes, are a valuable source of detailed information on the changes in regional terrestrial and aquatic environments at decadal to millennial scales. Here we present multi-proxy environmental records (pollen, dinoflagellate cysts, Ca, Ti and oxygen isotope data) from the uppermost 305 cm of the core 22-GC3 (42°13.53' N, 36°29.55' E) collected from a water depth of 838 m in the southern part of the Black Sea in 2007. The records span the last ~ 18 kyr (all ages are given in cal kyr BP). The pollen data reveal the dominance of the Artemisia-steppe in the region, suggesting rather dry/cold environments ~ 18-14.5 kyr BP. Warming/humidity increase during melt-water pulses (~ 16.1-14.5 kyr BP), indicated by d18O records from the 22-GC3 core sediment and from the Sofular Cave stalagmite, is expressed in more negative d13C values from the Sofular Cave, usually interpreted as the spreading of C3 plants. The records representing the interstadial complex (~ 14.5-12.9 kyr BP) show an increase in temperature and moisture, indicated by forest development, increased primary productivity and reduced surface run-off, whereas the switch from primary terrigenous to primary authigenic Ca origin occurs ~ 500 yr later. The Younger Dryas cooling is clearly demonstrated by more negative d13C values from the Sofular Cave and a reduction of pines. The early Holocene (11.7-8.5 kyr BP) interval reveals relatively dry conditions compared to the mostly moist and warm middle Holocene (8.5-5 kyr BP), which is characterized by the establishment of the species-rich warm mixed and temperate deciduous forests in the low elevation belt, temperate deciduous beech-hornbeam forests in the middle and cool conifer forest in upper mountain belt. The border between the early and middle Holocene in the vegetation records coincides with the opening of the Mediterranean corridor at ~ 8.3 kyr BP, as indicated by a marked change in the dinocyst assemblages and in the sediment lithology. Changes in the pollen assemblages indicate a reduction in forest cover after ~ 5 kyr BP, which was likely caused by increased anthropogenic pressure on the regional vegetation.

(Table 1) Age determination of Barents Sea sediment cores

he separate roles of oceanic heat advection and orbital forcing on influencing early Holocene temperature variability in the eastern Nordic Seas is investigated. The effect of changing orbital forcing on the ocean temperatures is tested using the 1DICE model, and the 1DICE results are compared with new and previously published temperature reconstructions from a transect of five cores located underneath the pathway of Atlantic water, from the Faroe-Shetland Channel in the south to the Barents Sea in the north. The stronger early Holocene summer insolation at high northern latitudes increased the summer mixed layer temperatures, however, ocean temperatures underneath the summer mixed layer did not increase significantly. The absolute maximum in summer mixed layer temperatures occurred between 9 and 6 ka BP, representing the Holocene Thermal Maximum in the eastern Nordic Seas. In contrast, maximum in northward oceanic heat transport through the Norwegian Atlantic Current occurred approximately 10 ka BP. The maximum in oceanic heat transport at 10 ka BP occurred due to a major reorganization of the Atlantic Ocean circulation, entailing strong and deep rejuvenation of the Atlantic Meridional Overturning Circulation, combined with changes in the North Atlantic gyre dynamic causing enhanced transport of heat and salt into the Nordic Seas.

(Table S1) Age determination of sediment cores SK99_17 and MD76-131

Productivity in the Arabian Sea is one of the highest in the world. It is controlled by seasonally reversing monsoonal wind-driven upwelling of nutrient-rich deeper waters which fuel phytoplankton growth. The detailed history of upwelling-induced productivity in the eastern Arabian Sea is unknown. Here we present paleoproductivity records from a composite sediment core at the millennial scale during the last 80 kyr B.P. These records are based on relative abundance counts of planktonic foraminifera and organic carbon contents, which are shown to mainly vary in concert. The eastern Arabian Sea upwelling-induced productivity was higher in the glacial period than in the Holocene, but it fell repeatedly on millennial timescales. These productivity declines occurred during cold events in the North Atlantic region, with the most pronounced changes prevailing during the Heinrich events. Hence, seasonal monsoon winds that drive upwelling-induced productivity in the east were weak when the North Atlantic was cold. These weak winds resulted in stratification of the water column, comparable to today's Arabian Sea stratification in the intermonsoonal period. Combining the new eastern with published western Arabian Sea results shows that the entire biological factory was severely diminished during the North Atlantic Heinrich events, and the seasonal productivity change in the Arabian Sea monsoon system was reduced with year-round low productivity.

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.

Chemistry and stable isotope record of sediment core GeoB4905-4

Evidence of rapid climatic oscillations like those observed in the Greenland ice cores and sediments from high latitudes of the northern Atlantic have been recognized in the pulses of terrigenous material to continental margin sediments off Cameroon. Fe/Ca ratios used as a parameter to quantify the relative proportions of terrigenous fluxes versus marine carbonate monitor the variability of the west African monsoon. They reveal the history of abrupt changes in precipitation over western and central Africa during the past 52 kyr. These rapid changes are particularly pronounced during the last glacial period and occur at timescales of a few thousand years. Stable oxygen isotope (delta18O) records of Globigerinoides ruber (pink) show high negative values reflecting periods of high monsoon precipitation. The Fe/Ca pattern is very similar to the Dansgaard-Oeschger cycles from the Greenland ice cores. The good correspondence between the warm interstadials of the Dansgaard-Oeschger cycles from the GISP2 ice core records and the high pulses of Fe/Ca sedimentation in our core suggest a strong teleconnection between the low-latitude African climate and the high-latitude northern hemisphere climate oscillations during the last glacial. This climatic link is probably vested in the west African monsoonal fluctuation that alters tropical sea surface temperatures, thermohaline circulations and in turn net export of heat from the south to the north Atlantic, coupled with the variability of the low-latitude southeast trade winds.

(Table 2) Age determination of Irminger Basin sediment cores

High-resolution records of coarse lithic content and oxygen isotope have been obtained in a piston core from the Irminger Basin. The last glacial period is characterized by numerous periods of increased iceberg discharges originating partly from Iceland and corresponding to millennial-scale instabilities of the coastal ice sheets and ice shelves in the Nordic area. A comparison with midlatitude sediment cores shows that ice-rafted material corresponding to the Heinrich events was deposited synchronously from 40° to 60°N. There are thus two oscillating systems: every 5-10 kyr massive iceberg armadas are released from large continental ice caps, whereas more frequent instabilities of the coastal ice sheets in the high latitude regions occur every 1.2-3.8 kyr. At the time of the Heinrich events the synchroneity of the response from all the northern hemisphere ice sheets attests the existence of strong interactions between the two systems.

Laminated sediment cores in the Bering Sea

During the last glacial termination, the upper North Pacific Ocean underwent dramatic and rapid changes in oxygenation that lead to the transient intensification of oxygen minimum zones (OMZs), recorded by the widespread occurrence of laminated sediments on circum-Pacific continental margins. We present a new laminated sediment record from the mid-depth (1100 m) northern Bering Sea margin that provides insight into these deglacial OMZ maxima with exceptional, decadal-scale detail. Combined ultrahigh-resolution micro-X-ray-fluorescence (micro-XRF) data and sediment facies analysis of laminae reveal an alternation between predominantly terrigenous and diatom-dominated opal sedimentation. The diatomaceous laminae are interpreted to represent spring/summer productivity events related to the retreating sea ice margin.We identified five laminated sections in the deglacial part of our site. Lamina counts were carried out on these sections and correlated with the Bølling-Allerød and Preboreal phases in the North Greenland Ice Core (NGRIP) oxygen isotope record, indicating an annual deposition of individual lamina couplets (varves). The observed rapid decadal intensifications of anoxia, in particular within the Bølling-Allerød, are tightly coupled to short-term warm events through increases in regional export production. This dependence of laminae formation on warmer temperatures is underlined by a correlation with published Bering Sea sea surface temperature records and d18O data of planktic foraminifera from the Gulf of Alaska. The rapidity of the observed changes strongly implies a close atmospheric teleconnection between North Pacific and North Atlantic regions.We suggest that concomitant increases in export production and subsequent remineralization of organic matter in the Bering Sea, in combination with oxygen-poor waters entering the Being Sea, drove down oxygen concentrations to values below 0.1ml/l and caused laminae preservation. Calculated benthic-planktic ventilation ages show no significant variations throughout the last deglaciation, indicating that changes in formation rates or differing sources of North Pacific mid-depth waters are not prime candidates for strengthening the OMZ at our site. The age models established by our correlation procedure allow for the determination of calendar age control points for the Bølling-Allerød and the Preboreal that are independent of the initial radiocarbon-based chronology. Resulting surface reservoir ages range within 730-990 yr during the Bølling-Allerød, 800-1100 yr in the Younger Dryas, and 765-775 yr for the Preboreal.