Magnetic susceptibility, XRF and color reflectance data of DSDP Hole 72-516F, cores 113 and 114

Deep marine successions of early Campanian age from DSDP site 516F drilled at low paleolatitudes in the South Atlantic reveal distinct sub-Milankovitch variability in addition to precession and eccentricity related variations. Elemental abundance ratios point to a similar 5 climatic origin for these variations and exclude a quadripartite structure - as observed in the Mediterranean Neogene - of the precession related cycles as an explanation for the inferred semi-precession cyclicity in MS. However, the semi-precession cycle itself is likely an artifact, reflecting the first harmonic of the precession signal. The sub-Milankovitch variability is best approximated by a ~ 7 kyr cycle as shown by 10 spectral analysis and bandpass filtering. The presence of sub-Milankovitch cycles with a period similar to that of Heinrich events of the last glacial cycle is consistent with linking the latter to low-latitude climate change caused by a non-linear response to precession induced variations in insolation between the tropics.

Age determination of sediment cores from the northern Mid-Atlantic Ridge

We construct age models for a suite of cores from the northeast Atlantic Ocean by means of accelerator mass spectrometer dating of a key core, BOFS 5K, and correlation with the rest of the suite. The effects of bioturbation and foraminiferal species abundance gradients upon the age record are modeled using a simple equation. The degree of bioturbation is estimated by comparing modeled profiles with dispersal of the Vedde Ash layer in core 5K, and we find a mixing depth of roughly 8 cm for sand-sized material. Using this value, we estimate that age offsets between unbioturbated sediment and some foraminifera species after mixing may be up to 2500 years, with lesser effect on fine carbonate (< 10 µm) ages. The bioturbation model illustrates problems associated with the dating of 'instantaneous' events such as ash layers and the 'Heinrich' peaks of ice-rafted detritus. Correlations between core 5K and the other cores from the BOFS suite are made on the basis of similarities in the downcore profiles of oxygen and carbon isotopes, magnetic susceptibility, water and carbonate content, and via marker horizons in X radiographs and ash beds.

Sediment, isotopic and geochemical record during Marine Isotope Stages 29 to 34 of IODP Site 339-U1387

Centennial-to-millennial scale records from IODP Site U1387, drilled during IODP Expedition 339 into the Faro Drift at 558 m water depth, now allow evaluating the climatic history of the upper core of the Mediterranean Outflow (MOW) and of the surface waters in the northern Gulf of Cadiz during the early Pleistocene. This study focuses on the period from Marine Isotope Stage (MIS) 29 to 34, i.e. the interval surrounding extreme interglacial MIS 31. Conditions in the upper MOW reflect obliquity, precession and millennial-scale variations. The benthic d18O signal follows obliquity with the exception of an additional, smaller d18O peak that marks the MIS 32/31 transition. Insolation maxima (precession minima) led to poor ventilation and a sluggish upper MOW core, whereas insolation minima were associated with enhanced ventilation and often also increased bottom current velocity. Millennial-scale periods of colder sea-surface temperatures (SST) were associated with short-term maxima in flow velocity and better ventilation, reminiscent of conditions known from MIS 3. A prominent contourite layer, coinciding with insolation cycle 100, was formed during MIS 31 and represents one of the few contourites developing within an interglacial period. MIS 31 surface water conditions were characterized by an extended period (1065-1091 ka) of warm SST, but SST were not much warmer than during MIS 33. Interglacial to glacial transitions experienced 2 to 3 stadial/interstadial cycles, just like their mid-to-late Pleistocene counterparts. Glacial MIS 30 and 32 recorded periods of extremely cold (< 12°C) SST that in their climatic impact were comparable to the Heinrich events of the mid and late Pleistocene. Glacial MIS 34, on the other hand, was a relative warm glacial period off southern Portugal. Overall, surface water and MOW conditions at Site U1387 show strong congruence with Mediterranean climate, whereas millennial-scale variations are closely linked to North Atlantic circulation changes.

Age deterimation of sediment cores from the Atlantic Ocean

High resolution benthic oxygen isotope records combined with radiocarbon datings, from cores retrieved in the North, Equatorial, and South Atlantic are used to establish a reliable cronostratigraphy for the last 60 ky. This common temporal framework enables us to study the timing of the sub-Milankovitch climate variability in the entire surface Atlantic during this period, as reflected in planktonic oxygen isotope records. Variations in sea surface temperatures in the Equatorial and South Atlantic reveal two warm periods during the mid-stage 3 which are correlated to the warming observed in the North Atlantic after Heinrich events (HL) 5 and 4. However, the records show that the warming started about 1500 y earlier in the South Atlantic. A zonally averaged ocean circulation model simulates a similar north-south thermal antiphasing between the latitudes of our coring sites, when pertubated by a freshwater flux anomaly. We infer that the observed phase relationship between the northern and the southern Atlantic is related to periods of reduced NADW production in the North Atlantic, such as during HL5 and HL4.

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.

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.