Supplementary material: XRF core scanning data, CaCO3 values, foraminiferal carbon and oxygen isotope data, Mg/Ca data, clay mineralogy and helium isotope results

The Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) is associated with abrupt climate change, carbon cycle perturbation, ocean acidification, as well as biogeographic shifts in marine and terrestrial biota that were largely reversed as the climatic transient waned. We report a clear exception to the behavior of the PETM as a reversing climatic transient in the eastern North Atlantic (Deep-Sea Drilling Project Site 401, Bay of Biscay) where the PETM initiates a greatly prolonged environmental change compared to other places on Earth where records exist. The observed environmental perturbation extended well past the d13C recovery phase and up to 650 kyr after the PETM onset according to our extraterrestrial 3He-based age-model. We observe a strong decoupling of planktic foraminiferal d18O and Mg/Ca values during the PETM d13C recovery phase, which in combination with results from helium isotopes and clay mineralogy, suggests that the PETM triggered a hydrologic change in western Europe that increased freshwater flux and the delivery of weathering products to the eastern North Atlantic. This state change persisted long after the carbon-cycle perturbation had stopped. We hypothesize that either long-lived continental drainage patterns were altered by enhanced hydrological cycling induced by the PETM, or alternatively that the climate system in the hinterland area of Site 401 was forced into a new climate state that was not easily reversed in the aftermath of the PETM.

Magnetic susceptibility of deep-sea cores from the Scotia Sea

We studied two deep-sea cores from the Scotia Sea to reconstruct past atmospheric circulation in the southern hemisphere and to resolve a long-standing debate on the interpretation of magnetic susceptibility (MS) records in Southern Ocean (SO) sediment. High-sedimentation sites MD07-3134 (0.2 - 1.2 m/kyr) and MD07-3133 (0.3 - 2 m/kyr) cover the last 92.5 kyr and 36 kyr, respectively. Both exhibit a one-to-one coupling of the MS and Ca2+ signal to the non-sea salt (nss) Ca2+ signal of the EDML ice core, clearly identifying atmospheric circulation as means of distribution. Comparison of additional proxies also excludes major influence by volcanic sources, sea-ice, icebergs, or oceanic current transport. The close resemblance of the dust proxies over the last glacial cycle, in turn, allows for the establishment of an age model of unprecedented resolution and precision for SO deep-sea sediment because atmospheric transport involves no major leads or lags. This is of particular importance because MS is routinely measured on deep-sea cores in the SO but the sediments usually lack biogenic carbonate and therefore had only limited stratigraphic control so far. Southern South America (SSA) is the likely source of eolian material because Site MD07-3133, located closer to the continent, has slightly higher MS values than Site MD07-3134, and also the MS record of Patagonian Site SALSA shows comparable variability. Patagonia was the dust source for both the Scotia Sea and East Antarctica. Dust fluxes were several times higher during glacial times, when atmospheric circulation was either stronger or shifted in latitude, sea level was lowered, shelf surfaces were exposed, and environmental conditions in SSA were dominated by glaciers and extended outwash plains. Hence, MS records of SO deep-sea sediment are reliable tracers of atmospheric circulation, allowing for chronologically-constrained reconstructions of the circum Antarctic paleoclimate history.

Radiocarbon ages and stable isotope measurements on sediment cores of the subarctic Pacific and its marginal seas

Under modern conditions only North Pacific Intermediate Water is formed in the northwest Pacific Ocean. This situation might have changed in the past. Recent studies with general circulation models indicate a switch to deep-water formation in the northwest Pacific during Heinrich Stadial 1 (17.5-15.0 ka) of the last glacial termination. Reconstructions of past ventilation changes based on paleoceanographic proxy records are still insufficient to test whether a deglacial mode of deep-water formation in the North Pacific Ocean existed. Here we present deglacial ventilation records based on radiocarbon-derived ventilation ages in combination with epibenthic stable carbon isotopes from the northwest Pacific including the Okhotsk Sea and Bering Sea, the two potential source regions for past North Pacific ventilation changes. Evidence for most rigorous ventilation of the intermediate-depth North Pacific occurred during Heinrich Stadial 1 and the Younger Dryas, simultaneous to significant reductions in Atlantic Meridional Overturning Circulation. Concurrent changes in d13C and ventilation ages point to the Okhotsk Sea as driver of millennial-scale changes in North Pacific Intermediate Water ventilation during the last deglaciation. Our records additionally indicate that changes in the d13C intermediate-water (700-1750 m water depth) signature and radiocarbon-derived ventilation ages are in antiphase to those of the deep North Pacific Ocean (>2100 m water depth) during the last glacial termination. Thus, intermediate- and deep-water masses of the northwest Pacific have a differing ventilation history during the last deglaciation.

The early Eocene global warming events: Implications from ODP Leg 208 Walvis Ridge depth transect

Five sections drilled in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge (SE Atlantic) during Ocean Drilling Program (ODP) Leg 208 resulted in the first complete early Paleogene deep-sea record. Here we present high-resolution stratigraphic records spanning a ~4.3 million yearlong interval of the late Paleocene to early Eocene. This interval includes the Paleocene-Eocene thermal maximum (PETM) as well as the Eocene thermal maximum (ETM) 2 event. A detailed chronology was developed with nondestructive X-ray fluorescence (XRF) core scanning records and shipboard color data. These records were used to refine the shipboard-derived spliced composite depth for each site and with a record from ODP Site 1051 were then used to establish a continuous time series over this interval. Extensive spectral analysis reveals that the early Paleogene sedimentary cyclicity is dominated by precession modulated by the short (100 kyr) and long (405 kyr) eccentricity cycles. Counting of precession-related cycles at multiple sites results in revised estimates for the duration of magnetochrons C24r and C25n. Direct comparison between the amplitude modulation of the precession component derived from XRF data and recent models of Earth's orbital eccentricity suggests that the onset of the PETM and ETM2 are related to a 100-kyr eccentricity maximum. Both events are approximately a quarter of a period offset from a maximum in the 405-kyr eccentricity cycle, with the major difference that the PETM is lagging and ETM2 is leading a 405-kyr eccentricity maximum. Absolute age estimates for the PETM, ETM2, and the magnetochron boundaries that are consistent with recalibrated radiometric ages and recent models of Earth's orbital eccentricity cannot be precisely determined at present because of too large uncertainties in these methods. Nevertheless, we provide two possible tuning options, which demonstrate the potential for the development of a cyclostratigraphic framework based on the stable 405-kyr eccentricity cycle for the entire Paleogene.

Barium intensities and tie points between holes of Leg 208

The Paleocene - Eocene thermal maximum (PETM) is one of the best known examples of a transient climate perturbation, associated with a brief, but intense, interval of global warming and a massive perturbation of the global carbon cycle from injection of isotopically light carbon into the ocean-atmosphere system. One key to quantifying the mass of carbon released, identifying the source(s), and understanding the ultimate fate of this carbon is to develop high-resolution age models. Two independent strategies have been employed, cycle stratigraphy and analysis of extraterrestrial Helium (HeET), both of which were first tested on Ocean Drilling Program (ODP) Site 690. Both methods are in agreement for the onset of the PETM and initial recovery, or the clay layer ("main body"), but seem to differ in the final recovery phase of the event above the clay layer, where the carbonate contents rise and carbon isotope values return toward background values. Here we present a state-of-the-art age model for the PETM derived from a new orbital chronology developed with cycle stratigraphic records from sites drilled during ODP Leg 208 (Walvis Ridge, Southeastern Atlantic) integrated with published records from Site 690 (Weddell Sea, Southern Ocean, ODP Leg 113). During Leg 208, five Paleocene - Eocene (P-E) boundary sections (Sites 1262 to 1267) were recovered in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge yielding the first stratigraphically complete P-E deep-sea sequence with moderate to relatively high sedimentation rates (1 to 3 cm/kyr). A detailed chronology was developed with non-destructive X-ray fluorescence (XRF) core scanning records on the scale of precession cycles, with a total duration of the PETM now estimated to be ~ 170 kyr. The revised cycle stratigraphic record confirms original estimates for the duration of the onset and initial recovery, but suggests a new duration for the final recovery that is intermediate to the previous estimates by cycle stratigraphy and HeET.

Geochemical and palynological measurements from the ODP Leg 189 at Site 1172 on the East Tasman Plateau

A brief (~150 kyr) period of widespread global average surface warming marks the transition between the Paleocene and Eocene epochs, ~56 million years ago. This so-called "Paleocene-Eocene thermal maximum" (PETM) is associated with the massive injection of 13C-depleted carbon, reflected in a negative carbon isotope excursion (CIE). Biotic responses include a global abundance peak (acme) of the subtropical dinoflagellate Apectodinium. Here we identify the PETM in a marine sedimentary sequence deposited on the East Tasman Plateau at Ocean Drilling Program (ODP) Site 1172 and show, based on the organic paleothermometer TEX86, that southwest Pacific sea surface temperatures increased from ~26 °C to ~33°C during the PETM. Such temperatures before, during and after the PETM are >10 °C warmer than predicted by paleoclimate model simulations for this latitude. In part, this discrepancy may be explained by potential seasonal biases in the TEX86 proxy in polar oceans. Additionally, the data suggest that not only Arctic, but also Antarctic temperatures may be underestimated in simulations of ancient greenhouse climates by current generation fully coupled climate models. An early influx of abundant Apectodinium confirms that environmental change preceded the CIE on a global scale. Organic dinoflagellate cyst assemblages suggest a local decrease in the amount of river run off reaching the core site during the PETM, possibly in concert with eustatic rise. Moreover, the assemblages suggest changes in seasonality of the regional hydrological system and storm activity. Finally, significant variation in dinoflagellate cyst assemblages during the PETM indicates that southwest Pacific climates varied significantly over time scales of 103 - 104 years during this event, a finding comparable to similar studies of PETM successions from the New Jersey Shelf.

Rock magnetic, element, and color data of three sediment cores along the Senegalese continental margin

We present a suite of new high-resolution records (0-135 ka) representing pulses of aeolian, fluvial, and biogenic sedimentation along the Senegalese continental margin. A multiproxy approach based on rock magnetic, element, and color data was applied on three cores enclosing the present-day northern limit of the ITCZ. A strong episodic aeolian contribution driven by stronger winds and dry conditions and characterized by high hematite and goethite input was revealed north of 13°N. These millennial-scale dust fluxes are synchronous with North Atlantic Heinrich stadials. Fluvial clay input driven by the West African monsoon predominates at 12°N and varies at Dansgaard-Oeschger time scales while marine productivity is strongly enhanced during the African humid periods and marine isotope stage 5. From latitudinal signal variations, we deduce that the last glacial ITCZ summer position was located between core positions at 12°26' and 13°40'N. Furthermore, this work also shows that submillennial periods of aridity over northwest Africa occurred more frequently and farther south than previously thought.

Climatic records of the late Holocene from sediment core GeoB9504-3 off Senegal

High-resolution climatic records of the late Holocene along the north-west African continental margin are scarce. Here we combine sediment grain size, elemental distribution and mineral assemblage data to trace dust and riverine sources at a shallow-marine sediment depocentre in the vicinity of the Senegal River mouth. The aim is to understand how these terrigenous components reflect climate variability during the late Holocene. Major element contents were measured and mineral identification was performed on three sub-fractions of our sediment core: (i) fluvial material <2 µm, (ii) aeolian material of 18-63 µm and (iii) a sub-fraction of dual-origin material of 2-18 µm. Results show that more than 80% of the total Al and Fe terrigenous bulk content is present in the fluviogenic fraction. In contrast, Ti, K and Si cannot be considered as proxies for one specific source off Senegal. The Al/Ca ratio, recording the continental river runoff, reveals two dry periods from 3010 to 2750 cal a BP and from 1900 to 1000 cal a BP, and two main humid periods from 2750 to 1900 cal a BP and from 1000 to 700 cal a BP. The match between (i) intervals of low river runoff inferred by low Al/Ca values, (ii) reduced river discharge inferred by integrated palynological data from offshore Senegal and (iii) periods of enhanced dune reactivation in Mali confirms this interpretation.

Geochemistry measured on sediment core GeoB12309-5

Laminated sediment records from the oxygen minimum zone in the Arabian Sea offer unique ultrahigh-resolution archives for deciphering climate variability in the Arabian Sea region. Although numerous analytical techniques are available it has become increasingly popular during the past decade to analyze relative variations of sediment cores' chemical signature by non-destructive X-ray fluorescence (XRF) core scanning. We carefully selected an approximately 5 m long sediment core from the northern Arabian Sea (GeoB12309-5: 24°52.3' N; 62°59.9' E, 956 m water depth) for a detailed, comparative study of high-resolution techniques, namely non-destructive XRF core scanning (0.8 mm resolution) and ICP-MS/OES analysis on carefully selected, discrete samples (1 mm resolution). The aim of our study was to more precisely define suitable chemical elements that can be accurately analyzed and to determine which elemental ratios can be interpretated down to sub-millimeter-scale resolutions. Applying the Student's t-test our results show significantly correlating (1% significance level) elemental patterns for all S, Ca, Fe, Zr, Rb, and Sr, as well as the K/Ca, Fe/Ti and Ti/Al ratios that are all related to distinct lithological changes. After careful consideration of all errors for the ICP analysis we further provide respective factors of XRF Core Scanner software error's underestimation by applying Chi-square-tests, which is especially relevant for elements with high count rates. As demonstrated by these new, ultra-high resolution data core scanning has major advantages (high-speed, low costs, few sample preparation steps) and represents an increasingly required alternative over the time consuming, expensive, elaborative, and destructive wet chemical analyses (e.g., by ICP-MS/OES after acid digestions), and meanwhile also provides high-quality data in unprecedented resolution.

Geochemistry of sediment core GeoB7920-2

The X-ray fluorescence (XRF) core scanner provides bulk-sediment chemistry data measured nondestructively at the split core sediment surface. Although this method is widely accepted, there is little known about the effects of physical properties such as density and water content on XRF core scanner data. Comparison of XRF scanner measurements from the sediment surface and dry powder samples of sediment core GeoB7920 indicates strongly reduced element intensities for the lighter elements Al and Si. We relate the lower element intensities of the measurements taken at the sediment surface to the amount of water in the sample volume analyzed by the XRF core scanner. The heavier elements K, Ca, Ti, and Fe remain relatively unaffected by the variation of any physical property within sediment core GeoB7920. Additionally, we successfully use the elemental intensity of Cl as a proxy for the seawater content in the sample volume analyzed by the XRF core scanner. This enables the establishment of a correction function for the elements Al and Si that corrects for the radiation absorption of the water content in sediment core GeoB7920 off Cape Blanc, NW Africa.

Bulk stable carbon and oxygen isotopes from ODP Site 1258, and Fe intensities from ODP Site 1262, and minimal tuning age model options for ODP Site 1258 and 1262

Timing is crucial to understanding the causes and consequences of events in Earth history. The calibration of geological time relies heavily on the accuracy of radioisotopic and astronomical dating. Uncertainties in the computations of Earth's orbital parameters and in radioisotopic dating have hampered the construction of a reliable astronomically calibrated time scale beyond 40 Ma. Attempts to construct a robust astronomically tuned time scale for the early Paleogene by integrating radioisotopic and astronomical dating are only partially consistent. Here, using the new La2010 and La2011 orbital solutions, we present the first accurate astronomically calibrated time scale for the early Paleogene (47-65 Ma) uniquely based on astronomical tuning and thus independent of the radioisotopic determination of the Fish Canyon standard. Comparison with geological data confirms the stability of the new La2011 solution back to ~54 Ma. Subsequent anchoring of floating chronologies to the La2011 solution using the very long eccentricity nodes provides an absolute age of 55.530 {plus minus} 0.05 Ma for the onset of the Paleocene/Eocene Thermal Maximum (PETM), 54.850 {plus minus} 0.05 Ma for the early Eocene ash -17, and 65.250 {plus minus} 0.06 Ma for the K/Pg boundary. The new astrochronology presented here indicates that the intercalibration and synchronization of U/Pb and 40Ar/39Ar radiometric geochronology is much more challenging than previously thought.

Sedimentary climate indicators of sediment core MD08-3167 and GeoB1711-4, southwestern Africa

The past climate evolution of southwestern Africa is poorly understood and interpretations of past hydrological changes are sometimes contradictory. Here we present a record of leaf-wax dD and View the MathML source taken from a marine sediment core at 23°S off the coast of Namibia to reconstruct the hydrology and C3 versus C4 vegetation of southwestern Africa over the last 140 000 years (140 ka). We find lower leaf-wax dD and higher View the MathML source (more C4 grasses), which we interpret to indicate wetter Southern Hemisphere (SH) summer conditions and increased seasonality, during SH insolation maxima relative to minima and during the last glacial period relative to the Holocene and the last interglacial period. Nonetheless, the dominance of C4 grasses throughout the record indicates that the wet season remained brief and that this region has remained semi-arid. Our data suggest that past precipitation increases were derived from the tropics rather than from the winter westerlies. Comparison with a record from the Congo Basin indicates that hydroclimate in southwestern Africa has evolved in antiphase with that of central Africa over the last 140 ka.

Stable oxygen isotope, UK'37, productivity, CaCO3 and XRF records from sediment core MD06-3075

Proxy records of hydrologic variability in the West Pacific Warm Pool (WPWP) have revealed wide-scale changes in past convective activity in response to orbital and sub-orbital climate forcings. However, attributing proxy responses to regional changes in WPWP hydrology versus local variations in precipitation requires independent records linking the terrestrial and marine realms. We present high-resolution stable isotope, UK'37 sea-surface temperature, X-ray fluorescence (XRF) core scanning and coccolithophore-derived paleoproductivity records covering the past 120 ka from International Marine Global Change (IMAGES) Program Core MD06-3075 (6°29' N, 125°50' E, water depth 1878 m), situated in the Davao Gulf on the southern side of Mindanao. XRF-derived log(Fe/Ca) records provide a robust proxy for runoff-driven sedimentary discharge from Mindanao, whilst past changes in local productivity are associated with variable freshwater runoff and stratification of the surface layer. Significant precessional-scale variability in sedimentary discharge occurred during Marine Isotope Stage (MIS) 5, with peaks in discharge contemporaneous with Northern Hemisphere summer insolation minima. We attribute these changes to the latitudinal migration of the Intertropical Convergence Zone (ITCZ) over the WPWP together with variability in the strength of the Walker circulation acting on precessional timescales. Between 60 and 15 ka sedimentary discharge at Mindanao was muted, displaying little orbital- or millennial-scale variability, likely in response to weakened precessional insolation forcing and lower sea level driving increased subsidence of air masses over the exposed Sunda Shelf. These results highlight the high degree of local variability in the precipitation response to past climate changes in the WPWP.