Centennial-scale evolution of Dansgaard-Oeschger events in sediment core MD95-2006

There is much uncertainty surrounding the mechanisms that forced the abrupt climate fluctuations found in many palaeoclimate records during Marine Isotope Stage (MIS)-3. One of the processes thought to be involved in these events is the Atlantic Meridional Overturning Circulation (MOC), which exhibited large changes in its dominant mode throughout the last glacial period. Giant piston core MD95-2006 from the northeast Atlantic Ocean records a suite of palaeoceanographic proxies related to the activity of both surface and deep water masses through a period of MIS-3 when abrupt climate fluctuations were extremely pronounced. A two-stage progression of surface water warming during interstadial warm events is proposed, with initial warming related to the northward advection of a thin warm surface layer within the North Atlantic Current, which only extended into deeper surface layers as the interstadial progressed. Benthic foraminifera isotope data also show millennial-scale oscillations but of a different structure to the abrupt surface water changes. These changes are argued to partly be related to the influence of low-salinity deepwater brines. The influence of deepwater brines over the site of MD95-2006 reached a maximum at times of rapid warming of surface waters. This observation supports the suggestion that brine formation may have helped to destabilize the accumulation of warm, saline surface waters at low latitudes, helping to force the MOC into a warm mode of operation. The contribution of deepwater brines relative to other mechanisms proposed to alter the state of the MOC needs to be examined further in future studies.

Sedimentary and carbonate preservation of ODP Sites 1262, 1263 and 1266

Rapid carbon input into the ocean-atmosphere system caused a dramatic shoaling of the lysocline during the Paleocene-Eocene thermal maximum (PETM), a transient (~170 kyr) global warming event that occurred roughly 55 Ma. Carbon cycle models invoking an accelerated carbonate-silicate feedback mechanism to neutralize ocean acidification predict that the lysocline would subsequently deepen to depths below its original position as the marine carbonate system recovered from such a perturbation. To test this hypothesis, records of carbonate sedimentation and preservation for PETM sections in the Weddell Sea (ODP Site 690) and along the Walvis Ridge depth transect (ODP Sites 1262, 1263, and 1266) were assembled within the context of a unified chronostratigraphy. The meridional gradient of undersaturation delimited by these records shows that dissolution was more severe in the subtropical South Atlantic than in the Weddell Sea during the PETM, a spatiotemporal pattern inconsistent with the view that Atlantic overturning circulation underwent a transient reversal. Deepening of the lysocline following its initial ascent is signaled by increases in %CaCO3 and coarse-fraction content at all sites. Carbonate preservation during the recovery period is appreciably better than that seen prior to carbon input with carbonate sedimentation becoming remarkably uniform over a broad spectrum of geographic and bathymetric settings. These congruent patterns of carbonate sedimentation confirm that the lysocline was suppressed below the depth it occupied prior to carbon input, and are consistent with the view that an accelerated carbonate-silicate geochemical cycle played an important role in arresting PETM conditions.

Atlantic 231Pa/230Th and biogenic opal compilation of the Holocene and the LGM

The strength and geometry of the Atlantic meridional overturning circulation is tightly coupled to climate on glacial-interglacial and millennial timescales, but has proved difficult to reconstruct, particularly for the Last Glacial Maximum. Today, the return flow from the northern North Atlantic to lower latitudes associated with the Atlantic meridional overturning circulation reaches down to approximately 4,000 m. In contrast, during the Last Glacial Maximum this return flow is thought to have occurred primarily at shallower depths. Measurements of sedimentary 231Pa/230Th have been used to reconstruct the strength of circulation in the North Atlantic Ocean, but the effects of biogenic silica on 231Pa/230Th-based estimates remain controversial. Here we use measurements of 231Pa/230Th ratios and biogenic silica in Holocene-aged Atlantic sediments and simulations with a two-dimensional scavenging model to demonstrate that the geometry and strength of the Atlantic meridional overturning circulation are the primary controls of 231Pa/230Th ratios in modern Atlantic sediments. For the glacial maximum, a simulation of Atlantic overturning with a shallow, but vigorous circulation and bulk water transport at around 2,000 m depth best matched observed glacial Atlantic 231Pa/230Th values. We estimate that the transport of intermediate water during the Last Glacial Maximum was at least as strong as deep water transport today.

(Table S1) Stable carbon and oxygen isotope ratios of benthic foraminifera of the western Atlantic Ocean

Oxygen and carbon isotopic data were produced on the benthic foraminiferal taxa Cibicidoides and Planulina from 25 new piston cores, gravity cores, and multicores from the Brazil margin. The cores span water depths from about 400 to 3000 m and intersect the major water masses in this region. These new data fill a critical gap in the South Atlantic Ocean and provide the motivation for updating the classic glacial western Atlantic d13C transect of Duplessy et al. (1988). The distribution of 13C of SumCO2 requires the presence of three distinct water masses in the glacial Atlantic Ocean: a shallow (~1000 m), southern source water mass with an end-member d13C value of about 0.3-0.5 per mil VPDB, a middepth (~1500 m), northern source water mass with an end-member value of about 1.5 per mil, and a deep (>2000 m), southern source water with an end-member value of less than -0.2 per mil, and perhaps as low as the -0.9 per mil values observed in the South Atlantic sector of the Southern Ocean (Ninnemann and Charles, 2002, doi:10.1016/S0012-821X(02)00708-2). The origins of the water masses are supported by the meridional gradients in benthic foraminiferal d18O. A revised glacial section of deep water d13C documents the positions and gradients among these end-member intermediate and deep water masses. The large property gradients in the presence of strong vertical mixing can only be maintained by a vigorous overturning circulation.

Nd, Pb and Sr isotopes in sediment cores from the Gulf of Cadiz and the Portuguese margin

Mediterranean Outflow Water (MOW) is characterised by higher temperatures and salinities than other ambient water masses. MOW spreads at water depths between 500 and 1500 m in the eastern North Atlantic and has been a source of salinity for the Atlantic Meridional Overturning Circulation in the North Atlantic. We used high-resolution Nd and Pb isotope records of past ambient seawater obtained from authigenic ferromanganese coatings of sediments in three gravity cores at 577, 1745 and 1974 m water depth in the Gulf of Cadiz and along the Portuguese margin complemented by a selection of surface sediments to reconstruct the extent and pathways of MOWover the past 23 000 years. The surface and downcore Nd isotope data from all water depths exhibit only a very small variability close to the present day composition of MOW but do not reflect the present day Nd isotopic stratification of the water column as determined from a nearby open ocean hydrographic station. In contrast, the Pb isotope records show significant and systematic variations, which provide evidence for a significantly different pattern of the MOW pathways between 20 000 and 12 000 years ago compared with the subsequent period of time.

Methane and stable isotope record of ice core EDML from Dronning Maud Land, Antarctica

Precise knowledge of the phase relationship between climate changes in the two hemispheres is a key for understanding the Earth's climate dynamics. For the last glacial period, ice core studies have revealed strong coupling of the largest millennial-scale warm events in Antarctica with the longest Dansgaard-Oeschger events in Greenland through the Atlantic meridional overturning circulation. It has been unclear, however, whether the shorter Dansgaard-Oeschger events have counterparts in the shorter and less prominent Antarctic temperature variations, and whether these events are linked by the same mechanism. Here we present a glacial climate record derived from an ice core from Dronning Maud Land, Antarctica, which represents South Atlantic climate at a resolution comparable with the Greenland ice core records. After methane synchronization with an ice core from North Greenland, the oxygen isotope record from the Dronning Maud Land ice core shows a one-to-one coupling between all Antarctic warm events and Greenland Dansgaard-Oeschger events by the bipolar seesaw. The amplitude of the Antarctic warm events is found to be linearly dependent on the duration of the concurrent stadial in the North, suggesting that they all result from a similar reduction in the meridional overturning circulation.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

(Supporting Information) 14Cage plateau tuning of local Atlantic sediment cores

Modeling studies predict that changes in radiocarbon (14C) reservoir ages of surface waters during the last deglacial episode will reflect changes in both atmospheric 14C concentration and ocean circulation including the Atlantic Meridional Overturning Circulation. Tests of these models require the availability of accurate 14C reservoir ages in well-dated late Quaternary time series. We here test two models using plateau-tuned 14C time series in multiple well-placed sediment core age-depth sequences throughout the lower latitudes of the Atlantic Ocean. 14C age plateau tuning in glacial and deglacial sequences provides accurate calendar year ages that differ by as much as 500-2500 years from those based on assumed global reservoir ages around 400 years. This study demonstrates increases in local Atlantic surface reservoir ages of up to 1000 years during the Last Glacial Maximum, ages that reflect stronger trades off Benguela and summer winds off southern Brazil. By contrast, surface water reservoir ages remained close to zero in the Cariaco Basin in the southern Caribbean due to lagoon-style isolation and persistently strong atmospheric CO2 exchange. Later, during the early deglacial (16 ka) reservoir ages decreased to a minimum of 170-420 14C years throughout the South Atlantic, likely in response to the rapid rise in atmospheric pCO2 and Antarctic temperatures occurring then. Changes in magnitude and geographic distribution of 14C reservoir ages of peak glacial and deglacial surface waters deviate from the results of Franke et al. (2008) but are generally consistent with those of the more advanced ocean circulation model of Butzin et al. (2012).

(Table 1) Age determination of sediment core MD99-2284

Dansgaard-Oeschger (D-O) cycles are the most dramatic, frequent, and wide-reaching abrupt climate changes in the geologic record. On Greenland, D-O cycles are characterized by an abrupt warming of 10 ± 5°C from a cold stadial to a warm interstadial phase, followed by gradual cooling before a rapid return to stadial conditions. The mechanisms responsible for these millennial cycles are not fully understood but are widely thought to involve abrupt changes in Atlantic Meridional Overturning Circulation due to freshwater perturbations. Here we present a new, high-resolution multiproxy marine sediment core monitoring changes in the warm Atlantic inflow to the Nordic seas as well as in local sea ice cover and influx of ice-rafted debris. In contrast to previous studies, the freshwater input is found to be coincident with warm interstadials on Greenland and has a Fennoscandian rather than Laurentide source. Furthermore, the data suggest a different thermohaline structure for the Nordic seas during cold stadials in which relatively warm Atlantic water circulates beneath a fresh surface layer and the presence of sea ice is inferred from benthic oxygen isotopes. This implies a delicate balance between the warm subsurface Atlantic water and fresh surface layer, with the possibility of abrupt changes in sea ice cover, and suggests a novel mechanism for the abrupt D-O events observed in Greenland ice cores.

Radionuclides, Nd isotopes and biogenic Opal in Atlanic Ocean sediments

Reconstructing past modes of ocean circulation is an essential task in paleoclimatology and paleoceanography. To this end, we combine two sedimentary proxies, Nd isotopes (epsilon-Nd) and the 231Pa/230Th ratio, both of which are not directly involved in the global carbon cycle, but allow the reconstruction of water mass provenance and provide information about the past strength of overturning circulation, respectively. In this study, combined 231Pa/230Th and epsilon-Nd down-core profiles from six Atlantic Ocean sediment cores are presented. The data set is complemented by the two available combined data sets from the literature. From this we derive a comprehensive picture of spatial and temporal patterns and the dynamic changes of the Atlantic Meridional Overturning Circulation over the past ~25 ka. Our results provide evidence for a consistent pattern of glacial/stadial advances of Southern Sourced Water along with a northward circulation mode for all cores in the deeper (>3000 m) Atlantic. Results from shallower core sites support an active overturning cell of shoaled Northern Sourced Water during the LGM and the subsequent deglaciation. Furthermore, we report evidence for a short-lived period of intensified AMOC in the early Holocene.

Geochemistry of Fe-Mn oxyhydroxide samples of ODP Hole 172-1060A

The strength of the North Atlantic Meridional Overturning Circulation during climatically highly variable Marine Isotope Stage (MIS) 3 has attracted much attention in recent years. Here we present high-resolution Nd isotope compositions of past seawater derived from authigenic Fe-Mn oxyhydroxides recovered from drift sediments on the Blake Ridge in the deep western North Atlantic (ODP Leg 172, Site 1060, 3481 m water depth). The data cover the period from 45 to 35 ka BP, tracing circulation changes during major Heinrich iceberg discharge event 4 (H4, ~40-39 ka BP). The Nd isotope record suggests that there was no northern-source water (NSW) mass like modern NADW at the deeper part of Blake Ridge at any time between 45 and 35 ka. This is fundamentally different from the hydrographic situation during the Holocene where NADW extends below 4500 m at this location. The epsilon-Nd of past deep water recorded in the Blake Ridge sediments was least radiogenic during Dansgaard/Oeschger (D/O) Interstadial (IS) 8 (epsilon-Nd = -11.3) and most radiogenic immediately preceding IS 9 (epsilon-Nd = -9.8). More radiogenic compositions were also recorded during H4 (-10.2 <= epsilon-Nd <= -9.9). The Nd isotope variability in MIS 3 matches that of a physical bottom current strength reconstruction from the same location. Neither record follows the pattern of Northern Hemisphere D/O climatic cycles. In our record, reduced mixing with northern source waters started in stadial 12 and lasted until after H4 in stadial 9, followed by a rapid increase in NSW contribution thereafter. This major change in the Nd isotope record predates the iceberg discharge event Heinrich 4 by more than 3 ka indicating a shallowing of the water mass boundary between Glacial North Atlantic Intermediate Water and Southern Source Water beneath. This early change in bottom water properties at the deep Blake Ridge suggests that North Atlantic deep water advection may already have decreased several thousand years before the actual iceberg discharge event and associated freshening of the surface waters in the North Atlantic. The change can thus not be attributed to climatic events in the North Atlantic but may be related to changes in flux of deep water from the South.

Nd, Sr and Pb isotopes and clay minerals in sediment cores from the Gulf of Cadiz and the Portuguese margin

The assemblages of marine sediments on the SW Iberian shelf have been controlled by contributions from distinct sources, which have varied in response to environmental changes since the Last Glacial Maximum (LGM). The rapid, decadal scale Mediterranean overturning circulation permits mixing of suspended particles from the entire Mediterranean Sea. They are entrained into the suspended particulate matter (SPM) carried by Mediterranean Outflow Water (MOW), which enters the eastern North Atlantic through the Strait of Gibraltar and spreads at intermediate depths in the Gulf of Cadiz and along the Portuguese continental margin. Other major sediment sources that have contributed to the characteristics and budget of SPM along the flow path of MOW on the SW Iberian shelf are North African dust and river-transported particles from the Iberian Peninsula. To reconstruct climate- and circulation-driven changes in the supply of sediments over the past ~23000 cal yr B.P., radiogenic Nd, Sr and Pb isotope records of the clay-size sediment fraction were obtained from one gravity core in the Gulf of Cadiz (577 m water depth) and from two gravity cores on the Portuguese shelf (1745 m, 1974 m water depth). These records are supplemented by time series analyses of clay mineral abundances from the same set of samples. Contrary to expectations, the transition from the LGM to the Holocene was not accompanied by strong changes in sediment provenance or transport, whereas Heinrich Event 1 (H1) and the African Humid Period (AHP) were marked by significantly different isotopic signatures reflecting changes in source contributions caused by supply of ice rafted material originating from the North American craton during H1 and diminished supply of Saharan dust during the AHP. The data also reveal that the timing of variations in the clay mineral abundances was decoupled from that of the radiogenic isotope signatures.

Late Pliocene variations of the Mediterranean outflow

Late Pliocene changes in the advection of Mediterranean Outflow Water (MOW) derivates were reconstructed at northeast Atlantic DSDP/ODP sites 548 and 982 and compared to records of WMDW at West Mediterranean Site 978. Neodymium isotope (epsilon-Nd) values more positive than ~10.5/~ 11 reflect diluted MOW derivates that spread almost continuously into the northeast Atlantic from 3.7 to 2.55 Ma, reaching Rockall Plateau Site 982 from 3.63 to 2.75 Ma. From 3.4 to 3.3 Ma average MOW temperature and salinity increased by 2°-4 °C and ~1 psu both at proximal Site 548 and distal Site 982. The rise implies a rise in flow strength, coeval with a long-term rise in both west Mediterranean Sea surface salinity by almost 2 psu and average bottom water salinity (BWS) by ~1 psu, despite inherent uncertainties in BWS estimates. The changes were linked with major Mediterranean aridification and a drop in African monsoon humidity. In contrast to model expectations, the rise in MOW salt discharge after 3.4 Ma did not translate into improved ventilation of North Atlantic Deep Water, since it possibly was too small to significantly influence Atlantic Meridional Overturning Circulation. Right after ~2.95 Ma, with the onset of major Northern Hemisphere Glaciation, long-term average bottom water temperature (BWT) and BWS at Site 548 dropped abruptly by ~5 °C and ~1-2 psu, in contrast to more distal Site 982, where BWT and BWS continued to oscillate at estimates of ~2 °C and 1.5-2.5 psu higher than today until ~2.6 Ma. We relate the small-scale changes both to a reduced MOW flow and to enhanced dilution by warm waters of a strengthened North Atlantic Current temporarily replacing MOW derivates at Rockall Plateau.

(Table S1) Neodymium isotopic results from ODP Site 172-1063

The Atlantic Meridional Overturning Circulation (AMOC) plays an important role in the Northern Hemisphere climate system. Significant interest went into the question of how excessive freshwater input through melting of continental ice can affect its overturning vigor and, hence, heat supply, to higher northern latitudes. Such forcing can be tested by investigating its behavior during extreme iceberg discharge events into the open North Atlantic during the last glacial period, the so-called Heinrich events (HE). Here we present neodymium (Nd) isotope compositions of past seawater, a sensitive chemical water mass tag, extracted from sediments of Ocean Drilling Program Site 1063 in the western North Atlantic (Bermuda Rise), covering the period surrounding HE 2, the Last Glacial Maximum, and the early deglaciation. These data are compared with a record of the kinematic circulation tracer (231Pa/230Th)xs extracted from the same sediment core. Both tracers indicate significant circulation changes preceding intense ice rafting during HE 2 by almost 2 kyr. Moreover, the Nd isotope record suggests the presence of deeply ventilating North Atlantic Deep Water early during Marine Isotope Stage 2 until it was replaced by Southern Source Water at ~27 ka. The early switch to high (Pa/Th)xs and radiogenic epsilon-Nd in relation to intensified ice rafting during HE 2 suggests that ice rafting into the open North Atlantic during major HE 2 was preceded by an early change of the AMOC. This opens the possibility that variations in AMOC contributed to or even triggered the ice sheet instability rather than merely responding to it.

CO2 concentration and stable isotope ratios of three Antarctic ice cores covering the period from 149.4 - 1.5 kyr before 1950

We present new d13C measurements of atmospheric CO2 covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in d13C(atm) of 0.5 permil occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in d13C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO2]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of d13C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.