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.

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.

(Table 2) Ba/Ca ratios of Cibicides wuellerstorfi and Orbulina universa from Bahaman sediment cores

This geochemical investigation utilizes Ba/Ca in the benthic foraminifer Cibicides wuellerstorfi from cores taken from the Bahama Banks and the Caribbean Sea to reconstruct changes in basal thermocline ventilation (800-1000 m) and middepth thermohaline circulation (1000-2000 m) in the western North Atlantic during the last glacial period, focusing on the deglacial transition. Previous studies show that an increase in ventilation of the North Atlantic subtropical gyre during the Last Glacial Maximum (LGM) caused a 30-60% decrease in labile nutrients within the thermocline layer. Using foraminiferal Ba/Ca as a proxy of refractory nutrients, increased ventilation during the LGM produced a depletion of less than 20% compared to Holocene values. Following glaciation, the production of Glacial North Atlantic Intermediate Water (GNAIW) shut down owing to the presence of meltwater in the surface ocean, which resulted in a decrease in ventilation, as seen by an enrichment of barium in the basal thermocline. GNAIW was subsequently replaced by barium-rich southern component water in the middepth western North Atlantic. Foraminiferal Ba/Ca data suggest a 38% contribution from southern component water to a depth as shallow as 1475 m and a 14% contribution at 1123 m during deglaciation.

(Table 2) Radiocarbon dating on cold-water corals collected in the Mediterranean Sea

This study presents newly obtained coral ages of the cold-water corals Lophelia pertusa and Madrepora oculata collected in the Alboran Sea and the Strait of Sicily (Urania Bank). These data were combined with all available Mediterranean Lophelia and Madrepora ages compiled from literature to conduct a basin-wide assessment of the spatial and temporal occurrence of these prominent framework-forming scleractinian species in the Mediterranean realm and to unravel the palaeo-environmental conditions that controlled their proliferation or decline. For the first time special focus was placed on a closer examination of potential differences occurring between the eastern and western Mediterranean sub-basins. Our results clearly demonstrate that cold-water corals occurred sparsely in the entire Mediterranean during the last glacial before becoming abundant during the Bølling-Allerød warm interval, pointing to a basin-wide, almost concurrent onset in (re-)colonisation after ~13.5 ka. This time coincides with a peak in meltwater discharge originating from the northern Mediterranean borderlands which caused a major reorganisation of the Mediterranean thermohaline circulation. During the Younger Dryas and Holocene, some striking differences in coral proliferation were identified between the sub-basins such as periods of highly prolific coral growth in the eastern Mediterranean Sea during the Younger Dryas and in the western basin during the Early Holocene, whereas a temporary pronounced coral decline during the Younger Dryas was exclusively affecting coral sites in the Alboran Sea. Comparison with environmental and oceanographic data revealed that the proliferation of the Mediterranean corals is linked with enhanced productivity conditions. Moreover, corals thrived in intermediate depths and showed a close relationship with intermediate water mass circulation in the Mediterranean sub-basins. For instance, reduced Levantine Intermediate Water formation hampered coral growth in the eastern Mediterranean Sea during sapropel S1 event as reduced Winter Intermediate Water formation did in the westernmost part of the Mediterranean (Alboran Sea) during the Mid-Holocene. Overall, this study clearly demonstrates the importance to consider region-specific environmental changes as well as species-specific environmental preferences in interpreting coral chronologies. Moreover, it highlights that the occurrence or decline of cold-water corals is not controlled by one key parameter but rather by a complex interplay of various environmental variables.

Age determinaton and age model of sediments from the Bay of Biscay

Continuous high-resolution mass accumulation rates (MAR) and X-ray fluorescence (XRF) measurements from marine sediment records in the Bay of Biscay (NE Atlantic) have allowed the determination of the timing and the amplitude of the 'Fleuve Manche' (Channel River) discharges during glacial stages MIS 10, MIS 8, MIS 6 and MIS 4-2. These results have yielded detailed insight into the Middle and Late Pleistocene glaciations in Europe and the drainage network of the western and central European rivers over the last 350 kyr. This study provides clear evidence that the 'Fleuve Manche' connected the southern North Sea basin with the Bay of Biscay during each glacial period and reveals that 'Fleuve Manche' activity during the glaciations MIS 10 and MIS 8 was significantly less than during MIS 6 and MIS 2. We correlate the significant 'Fleuve Manche' activity, detected during MIS 6 and MIS 2, with the extensive Saalian (Drenthe Substage) and the Weichselian glaciations, respectively, confirming that the major Elsterian glaciation precedes the glacial MIS 10. In detail, massive 'Fleuve Manche' discharges occurred at ca 155 ka (mid-MIS 6) and during Termination I, while no significant discharges are found during Termination II. It is assumed that a substantial retreat of the European ice sheet at ca 155 kyr, followed by the formation of ice-free conditions between the British Isles and Scandinavia until Termination II, allowed meltwater to flow northwards through the North Sea basin during the second part of the MIS 6. We assume that this glacial pattern corresponds to the Warthe Substage glacial maximum, therefore indicating that the data presented here equates to the Drenthe and the Warthe glacial advances at ca 175-160 ka and ca 150-140 ka, respectively. Finally, the correlation of our records with ODP site 980 reveals that massive 'Fleuve Manche' discharges, related to partial or complete melting of the European ice masses, were synchronous with strong decreases in both the rate of deep-water formation and the strength of the Atlantic thermohaline circulation. 'Fleuve Manche' discharges over the last 350 kyr probably participated, with other meltwater sources, in the collapse of the thermohaline circulation by freshening the northern Atlantic surface water.

(Table 2) Radiocarbon content and 14C AMS ages of planktonic foraminifera from sediment cores Louis_1900 and Louis_2023

North American freshwater runoff records have been used to support the case that climate flickers were caused by shutdowns of the ocean thermohaline circulation (THC) resulting from reversals of meltwater discharges. Inconsistencies in the documentation of these meltwater switches, however, continue to fuel the debate on the cause/s of the oscillatory nature of the deglacial climate. New oxygen and carbon isotope records from the northern Gulf of Mexico depict in exceptional detail the succession of meltwater floods and pauses through the southern routing during the interval 16 to 8.9 ka (14C years BP; ka, kiloannum). The records underscore the bimodal role played by the Gulf of Mexico as a destination of meltwater discharges from the receding Laurentide Ice Sheet. The evidence indicates that the Gulf of Mexico acted as the principal source of superfloods at 13.4, 12.6, and 11.9 ka that reached the North Atlantic and contributed significantly to density stratification, disruption of ocean ventilation, and cold reversals. Gulf of Mexico lapsed into a "relief valve" position in post-Younger Dryas time, when meltwater discharges were rerouted south at 9.9, 9.7, 9.4, and 9.1 ka, thus temporarily interrupting North Atlantic-bound freshwater discharges from Lake Agassiz. The history of meltwater events in the Gulf of Mexico contradicts the model that meltwater flow via the eastern outlets into the North Atlantic disrupted the ocean THC, causing cooling, while diversions to the Gulf of Mexico via the Mississippi River enhanced THC and warming.

(Table 1) Oxygen isotopic composition of interstitial waters from ODP Hole 161-976

Interstitial waters recovered from Ocean Drilling Program, Leg 161, site 976 in the western Mediterranean Sea are used in conjunction with a numerical model to constrain the delta18O of seawater in the basin since the Last Glacial Maximum, including Sapropel Event 1. To resolve the oxygen isotopic composition of the deep Mediterranean, we use a model that couples fluid diffusion with advective transport, thus producing a profile of seawater delta18O variability that is unaffected by glacial-interglacial variations in marine temperature. Comparing our reconstructed seawater delta18O to recent determinations of 1.0 per mil for the mean ocean change in glacial-interglacial delta18O due to the expansion of global ice volume, we calculate an additional 0.2 per mil increase in Mediterranean delta18O caused by local evaporative enrichment. This estimate of delta18O change, due to salinity variability, is smaller than previous studies have proposed and demonstrates that Mediterranean records of foraminiferal calcite delta18O from the last glacial period include a strong temperature component. Paleotemperatures determined in combination with a stacked record of foraminiferal calcite depict almost 9°C of regional cooling for the Last Glacial Maximum. Model results suggest a decrease of ~1.1 per mil in seawater delta18O relative to the modern value caused by increased freshwater input and reduced salinity accompanying the formation of the most recent sapropel. The results additionally indicate the existence of isotopically light water circulating down to bottom water depths, at least in the western Mediterranean, supporting the existence of an 'anti-estuarine' thermohaline circulation pattern during Sapropel Event 1.

Stable oxygen isotope record of benthic foraminifera from the tropical Atlantic

Benthic foraminiferal oxygen isotope ratios from two sediment cores recovered at 426 and 1299 m water depth in the eastern and western tropical Atlantic show that a slowdown of the thermohaline circulation (THC) during Heinrich event H1 and the Younger Dryas was accompanied by rapid and intense warming of intermediate depth waters. Millennial-scale covariations of low paleosalinities in the subpolar North Atlantic with decreased benthic oxygen isotope ratios in the eastern tropical Atlantic throughout the past 10,000 years suggest that THC weakening might be related to middepth warming during the Holocene period as well. Climate model experiments simulating a strong reduction of the THC in the Atlantic Ocean under present-day and glacial conditions reveal that the increase of temperature in the middepth tropical and South Atlantic is a common feature for both climatic states, caused by a reduced ventilation of cold intermediate and deep waters in conjunction with downward mixing of heat from the thermocline. From the similarity of the paleoclimatic records with the model simulations, we infer that the characteristic pattern of temperature change in the Atlantic Ocean related to weakened thermohaline circulation can serve as an indicator of present-day and future THC slowdown.

Sea surface temperature reconstruction for sediment core M35003-4 in the western tropical North Atlantic

Planktonic foraminiferal census counts are used to construct high-resolution sea surface temperature (SST) and subsurface (thermocline) temperature records at a core site in the Tobago Basin, Lesser Antilles. The record is used to document climatic variability at this tropical site in comparison to middle- and high-latitude sites and to test current concepts of cross-equatorial heat transports as a major player in interhemispheric climate variability. Temperatures are estimated using transfer function and modern analog techniques. Glacial - maximum cooling of 2.5°-3°C is indicated; maximum cooling by 4°C is inferred for isotope stage 3. The SST record displays millennial-scale variability with temperature jumps of up to 3°C and closely tracks the structure of ice-core Dansgaard/Oeschger cycles. SST variations in part of the record run opposite to the SST evolution at high northern latitude sites, pointing to thermohaline circulation and marine heat transport as an important factor driving SST in the tropical and high-latitude Atlantic, both on orbital and suborbital timescales.

Stable isotopes, radionuclides, and calculated sea surface temperature of two sediment cores in the Sicily Channel

The Eastern Mediterranean Transient (EMT) occurred in the Aegean Sea from 1988 to 1995 and is the most significant intermediate-to-deep Mediterranean overturning perturbation reported by instrumental records. The EMT was likely caused by accumulation of high salinity waters in the Levantine and enhanced heat loss in the Aegean Sea, coupled with surface water freshening in the Sicily Channel. It is still unknown whether similar transients occurred in the past and, if so, what their forcing processes were. In this study, sediments from the Sicily Channel document surface water freshening (SCFR) at 1910±12, 1812±18, 1725±25 and 1580±30 CE. A regional ocean hindcast links SCFR to enhanced deep-water production and in turn to strengthened Mediterranean thermohaline circulation. Independent evidence collected in the Aegean Sea supports this reconstruction, showing that enhanced bottom water ventilation in the Eastern Mediterranean was associated with each SCFR event. Comparison between the records and multi-decadal atmospheric circulation patterns and climatic external forcings indicates that Mediterranean circulation destabilisation occurs during positive North Atlantic Oscillation (NAO) and negative Atlantic Multidecadal Oscillation (AMO) phases, reduced solar activity and strong tropical volcanic eruptions. They may have recurrently produced favourable deep-water formation conditions, both increasing salinity and reducing temperature on multi-decadal time scales.

(Table 1) Nd isotope ratios of ODP Leg 208 samples

We present new Nd isotope records from Walvis Ridge Ocean Drilling Program (ODP) sites 1262-1264 (southeastern Atlantic) spanning the past 24 Ma to investigate the Neogene evolution of Atlantic thermohaline circulation. The new data indicate that deepwater epsilon-Nd(t) values from ODP Site 1262 decrease from -11.0 at 10.6 Ma to -12.5 by 7.3 Ma. This decrease parallels the Nd isotope trends contained in Fe-Mn crust records from the northwestern Atlantic; however, the shift at ODP Site 1262 (4755 m water depth) occurred up to ~6 Ma earlier than the corresponding changes in crust records from the Atlantic and earlier than any Atlantic site shallower than 2700 m paleowater depth. Recent interpretations of the rapid decrease in Fe-Mn crust Nd epsilon-Nd(t)ss values invoke changes in weathering inputs to the Labrador Sea region rather than a fundamental change in deepwater convection in the Labrador Sea. However, the new evidence for significant depth stratification of the Nd isotope signal in the southeastern Atlantic between 10.6 and 7.3 Ma suggests that the onset of deepwater convection in the Labrador Sea may have played a role in the deepwater decrease in Nd isotopic composition. Climatic conditions during the middle to late Miocene likely favored an increase in the importance of glacially induced mechanical weathering, while at the same time promoting deep convection in the Labrador Sea.

Age determination of sediment core W8709A-13

Northeast Pacific benthic foraminiferal d18O and d13 reveal repeated millennial-scale events of strong deep-sea ventilation (associated with nutrient depletion and/or high gas exchange) during stadial (cool, high ice volume) episodes from 10 to 60 ka, opposite the pattern in the deep North Atlantic. Two climate mechanisms may explain this pattern. North Pacific surface waters, chilled by atmospheric transmission from a cold North Atlantic and made saltier by reduced freshwater vapor transports, could have ventilated the deep Pacific from above. Alternatively, faster turnover of Pacific bottom and mid-depth waters, driven by Southern Ocean winds, may have compensated for suppressed North Atlantic Deep Water production during stadial intervals. During the Younger Dryas event (~11.6-13.0 cal ka), ventilation of the deep NE Pacific (~2700 m) lagged that in the Santa Barbara Basin (~450 m) by >500 years, suggesting that the NE Pacific was first ventilated at intermediate depth from above and then at greater depth from below. This apparent lag may reflect the adjustment time of global thermohaline circulation.

(Table 1) Age comparison between the GPTS and ODP Site 181-1123 age model for each chron base within the Middle Miocene section

Today the deep western boundary current (DWBC) east of New Zealand is the most important route for deep water entering the Pacific Ocean. Large-scale changes in deep water circulation patterns are thought to have been associated with the development of the East Antarctic Ice Sheet (EAIS) close to the main source of bottom water for the DWBC. Here we reconstruct the changing speed of the southwest Pacific DWBC during the middle Miocene from ~15.5-12.5 Ma, a period of significant global ice accumulation associated with EAIS growth. Sortable silt mean grain sizes from Ocean Drilling Program Site 1123 reveal variability in the speed of the Pacific inflow on the timescale of the 41 kyr orbital obliquity cycle. Similar orbital period flow changes have recently been demonstrated for the Pleistocene epoch. Collectively, these observations suggest that a strong coupling between changes in the speed of the deep Pacific inflow and high-latitude climate forcing may have been a persistent feature of the global thermohaline circulation system for at least the past 15 Myr. Furthermore, long-term changes in flow speed suggest an intensification of the DWBC under an inferred increase in Southern Component Water production. This occurred at the same time as decreasing Tethyan outflow and major EAIS growth between ~15.5 and 13.5 Ma. These results provide evidence that a major component of the deep thermohaline circulation was associated with the middle Miocene growth of the EAIS and support the view that this time interval represents an important step in the development of the Neogene icehouse climate.

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