Constructing a stacked benthic d18O record

Composite stacks were constructed by superimposing 6 to 13 benthic foraminiferal d18O records covering the period 0-850 ka. An initial timescale for each core was established using radioisotopic age control points and assuming constant sedimentation rates between these points. The average of these records is our 13-core "untuned" stack. Next, we matched the 41 kyr component of each record individually to variations in Earth's obliquity. Four of the 13 records produced timescales that were inconsistent with one or more of the known radioisotopic ages. The nine remaining cores were averaged to create a "minimally tuned" stack. Six of the minimally tuned cores were assembled into a "tropical" stack. For each stack we estimated the uncertainty envelope from the standard deviation of the constituents. Spectral analysis of the three stacks indicates that benthic d18O is dominated by a 100 kyr oscillation that has a narrow spectral peak. The contribution of precession to the total variance is small when compared to prior results from planktic stacks.

Response of nannoplankton to major changes in sea-surface temperature at the South Chatham Rise, southeastern New Zealand

A high-resolution history of paleoceanographic changes in the subpolar waters of the southern margin of the Subtropical Convergence Zone during the last 130 kyr, is present in foraminiferal assemblages of DSDP Site 594. The foraminifera indicate that sea-surface temperatures during the Last Interglacial Climax were warmer than today, and that between substage 5d through to the end of isotope stage 2, temperatures were mostly cooler than Holocene temperatures. The paleotemperatures suggest that (1) the Subtropical Convergence was located over the site during substage 5e, later moving further north, then moving southwards to near the site during the Holocene, and (2) the Polar Front was positioned over the Site during glacial stages 6, 4, 2 and possibly parts of stage 3. Several major events are indicated by the nannofloral assemblages during these large changes in sea-surface temperature and associated reorganization of ocean circulation. First, the time-progressive trends between E. huxleyi and medium to large Gephyrocupsa are unique to this site, with E. huxleyi dominating over medium Gephyrocupsa during stages 5c-a, middle part of stage 4 and after the middle point of stage 3. This unusual trend may (at least partly) be caused by the shift of the Polar Front across the site. Second, upwelling flora (E. huxleyi and small placoliths) increase in abundance during stages 1, 3 and 5, suggesting that upwelling or disturbance of water stratification took place during the interglacials. Thirdly, there are no significant differences between the distribution patterns of the various morphotypes of medium to large Gephyrocupsu, and the combined value of all medium Gephyrocupsu increases in abundance during glacials (stages 2 and 4 and the end of stage 6), similar to the abundance trends in benthic foraminifera. Finally, subordinate nannofossil taxa also show distinctive climatic trends during the last glacial cycle: (1) Syrucosphaera spp. are present in increased abundance during warmer extremes in climate (substages 5e, 5a, and stage 1); (2) Coccolithus pelagicus and Culcidiscus leptoporus dominate the subordinate nannofossil taxa, and their relative proportions seem to provide a useful paleoceanographic index, with C. pelagicus dominating when the Polar Front Zone is over the site (stages 6, 4 and 2), whilst C. leptoporus is relatively more abundant when the STC is positioned over the site (stages 1 and 5e). Increased abundance of C. pelagicus also can indicate intensified coastal upwelling.

CaCO3 content, terrigenous components, foraminifera abundance and age of sediment core MD95-2006, Barra Fan, UK

Lithology, lithic petrology, planktonic foraminiferal abundances, and clastic grain sizes have been determined in a 30 m-long core recovered from the Barra Fan off northwest Scotland. The record extends back to around 45 kyr B.P., with sedimentation rates ranging between 50 and 200 cm/kyr. The abundance of ice-rafted debris indicates 16 glacimarine events, including temporal equivalents to Heinrich events 1-4. Enhanced concentrations of basaltic material derived from the British Tertiary Province suggest that the glacimarine sediments record variations in a glacial source on the Hebrides shelf margin. Glacimarine zones are separated by silty intervals with high planktonic foraminifera concentrations that reflect an interstadial circulation regime in the Rockall Trough. The results suggest that the last British Ice Sheet fluctuated with a periodicity of 2000-3000 years, in common with the Dansgaard-Oeschger climate cycle.

Description and chemical composition of manganese nodules during R/V "Tyro" cruise APNAP 1 to the Madeira Abyssal Plain in 1986

The porewater and sediment composition of two boxcores and of a small gravity core, taken on a manganese-nodule-covered hill and in the Madeira Abyssal Plain proper respectively, are compared. The pore-water study of the two boxcores indicates that oxic conditions prevail in both cores. In addition, it indicates that no detectable fluxes of Mn or Fe occur from the porewater to the ocean bottom water. Variations in the geochemical composition of the sediments can be explained by fluctuations in the amount of carbonate, which acts as a diluting agent. A clear carbonate minimum is observed at 20-22 cm depth in the two cores. This minimum is likely to be associated with the last glacial period (10-20 kyr B.P.). This association is supported by the sediment accumulation rate of 15 mm/kyr as found by extrapolation from the rate for pelagic sediments in the Madeira Abyssal Plain. The bulk composition of the manganese nodules recovered from the submarine hill is chemically almost identical to the average composition of Atlantic nodules. The trace metal and Rare Earth Elements composition indicate a hydrogenous origin for the manganese nodules of this study. On the basis of the chemical composition, and that of nodules relative to that of the adjacent sediments, an average nodule accretian rate of 2.8-3.3 mm/myr has been calculated. Although the analyses of the entire ferromanganese nodules that have been studied seem to indicate a homogenous composition, internal structures of the nodules reveal great inhomogeneity, both visually and chemically. These fluctuations may be related to variations in the fluxes of Mn and Fe, which in turn could be climate-related.

Mg/Ca and stable oxygen measurement of the past 27 Ma of ODP Hole 120-747A

We explore the applicability of paired Mg/Ca and 18O/16O measurements on benthic foraminifera from Southern Ocean site 747 to paleoceanographic reconstructions on pre-Pleistocene timescales. We focus on the late Oligocene through Pleistocene (27-0 Ma) history of paleotemperatures and the evolution of the d18O values of seawater (d18Osw) at a temporal resolution of ~100-200 kyr. Absolute paleotemperature estimates depend on assumptions of how Mg/Ca ratios of seawater have changed over the past 27 Myr, but relative changes that occur on geologically brief timescales are robust. Results indicate that at the Oligocene to Miocene boundary (23.8 Ma), temperatures lag the increase in global ice-volume deduced from benthic foraminiferal d18O values, but the smaller-scale Miocene glaciations are accompanied by ocean cooling of -1°C. During the mid-Miocene phase of Antarctic ice sheet growth (~15-13 Ma), water temperatures cool by ~3°C. Unlike the benthic foraminiferal d18O values, which remain relatively constant thereafter, temperatures vary (by 3°C) and reach maxima at ~12 and ~8.5 Ma. The onset of significant Northern Hemisphere glaciation during the late Pliocene is synchronous with an ~4°C cooling at site 747. A comparison of our d18Osw curve to the Haq et al. (1987, doi:10.1126/science.235.4793.1156 ) sea level curve yields excellent agreement between sequence boundaries and times of increasing seawater 18O/16O ratios. At ~12-11 Ma in particular, when benthic foraminiferal d18O values do not support a further increase in ice volume, the d18Osw curve comes to a maximum that corresponds to a major mid-Miocene sea level regression. The agreement between the character of our Mg/Ca-based d18Osw curve and sequence stratigraphy demonstrates that benthic foramaniferal Mg/Ca ratios can be used to trace the d18Osw on pre-Pleistocene timescales despite a number of uncertainties related to poorly constrained temperature calibrations and paleoseawater Mg/Ca ratios. The Mg/Ca record also highlights that deep ocean temperatures can vary independently and unexpectedly from ice volume changes, which can lead to misinterpretations of the d18O record.

Neodymium isotope composition of planktonic foraminifera from ODP Hole 121-758A

This study presents neodymium isotope and elemental data for cleaned planktonic foraminifera from ODP site 758 in the southernmost reaches of the Bay of Bengal in the north-east Indian Ocean. Cleaning experiments using oxidative-reductive techniques suggest that diagenetic Fe-Mn oxyhydroxide coatings can be effectively removed, and that the measured Nd isotope composition reflects the composition of seawater from which the foraminiferal calcium carbonate was precipitated. Modern core-top Pulleniatina obliquiloculata and Globorotalia menardii give epsilon-Nd values of 310.12 +/- 0.16 and 310.28 +/- 0.16, respectively, indistinguishable from recent direct measurements of surface seawater in this area. A high-resolution Nd isotope record obtained from G. menardii for the past 150 kyr shows systematic variations (Delta epsilon-Nd = 3) on glacial-interglacial timescales. The timing of those variations shows a remarkable correspondence with the global oxygen isotope record, which suggests a process controlling the Nd isotope composition that responds in phase with global climate cycles. Palaeoclimate reconstruction indicates that during the last glacial maximum changes in monsoon circulation resulted in a reduction in rainfall over the Indian subcontinent, and a decrease in the flux of river water delivered to the Bay of Bengal. Thus, changes in the riverine input of Nd, a change in either flux or composition, most likely caused the isotope variations, although changes in dust source or local ocean circulation may have also played a role. These results clearly establish a link between climate change and variations in radiogenic isotopes in the oceans, and illustrate the potential of Nd isotopes in foraminifera for highresolution palaeoceanographic reconstruction.

Early Eocene benthic foraminifera of DSDP Hole 48-401

Within the last decade, several early Eocene hyperthermals have been detected globally. These transient warming events have mainly been characterized geochemically - using stable isotopes, carbonate content measurements or XRF core scanning - yet detailed micropaleontological records are sparse, limiting our understanding of the driving forces behind hyperthermals and of the contemporaneous paleoceanography. Here, detailed geochemical and quantitative benthic foraminiferal records are presented from lower Eocene pelagic sediments of Deep Sea Drilling Project Site 401 (Bay of Biscay, northeast Atlantic). In calcareous nannofossil zone NP11, several clay-enriched levels correspond to negative d13C and d18O bulk-rock excursions with amplitudes of up to ~0.75 per mil, suggesting that significant injections of 12C-enriched greenhouse gasses and small temperature rises took place. Coeval with several of these hyperthermal events, the benthic foraminiferal record reveals increased relative abundances of oligotrophic taxa (e.g. Nuttallides umbonifera) and a reduction in the abundance of buliminid species followed by an increase of opportunistic taxa (e.g. Globocassidulina subglobosa and Gyroidinoides spp.). These short-lived faunal perturbations are thought to be caused by reduced seasonality of productivity resulting in a decreased Corg flux to the seafloor. Moreover, the sedimentological record suggests that an enhanced influx of terrigenous material occurred during these events. Additionally, the most intense d13C decline (here called level d) gives rise to a small, yet pronounced long-term shift in the benthic foraminiferal composition at this site, possibly due to the reappraisal of upwelling and the intensification of bottom water currents. These observations imply that environmental changes during (smaller) hyperthermal events are also reflected in the composition of deep-sea benthic communities on both short (<100 kyr) and longer time scales. We conclude that the faunal patterns of the hyperthermals observed at Site 401 strongly resemble those observed in other deep-sea early Paleogene hyperthermal deposits, suggesting that similar processes have driven them.

Abrupt shifts of the Sahara-Sahel boundary during Heinrich stadials demonstrated on a sediment core transect

Relict dune fields that are found as far south as 14° N in the modern-day African Sahel are testament to equatorward expansions of the Sahara desert during the Late Pleistocene. However, the discontinuous nature of dune records means that abrupt millennial-timescale climate events are not always resolved. High-resolution marine core studies have identified Heinrich stadials as the dustiest periods of the last glacial in West Africa although the spatial evolution of dust export on millennial timescales has so far not been investigated. We use the major-element composition of four high-resolution marine sediment cores to reconstruct the spatial extent of Saharan-dust versus river-sediment input to the continental margin from West Africa over the last 60 ka. This allows us to map the position of the sediment composition corresponding to the Sahara-Sahel boundary. Our records indicate that the Sahara-Sahel boundary reached its most southerly position (13° N) during Heinrich stadials and hence suggest that these were the periods when the sand dunes formed at 14° N on the continent. Heinrich stadials are associated with cold North Atlantic sea surface temperatures which appear to have triggered abrupt increases of aridity and wind strength in the Sahel. Our study illustrates the influence of the Atlantic meridional overturning circulation on the position of the Sahara-Sahel boundary and on global atmospheric dust loading.

(Tabls S1) Bulk geochemistry for sediment samples collected from the Cenomanian-Turonian boundary section at Pont d'Issole

Oceanic Anoxic Event 2 (OAE2), spanning the Cenomanian-Turonian boundary (CTB), represents one of the largest perturbations in the global carbon cycle in the last 100 Myr. The d13Ccarb, d13Corg, and d18O chemostratigraphy of a black shale-bearing CTB succession in the Vocontian Basin of France is described and correlated at high resolution to the European CTB reference section at Eastbourne, England, and to successions in Germany, the equatorial and midlatitude proto-North Atlantic, and the U.S. Western Interior Seaway (WIS). Delta13C (offset between d13Ccarb and d13Corg) is shown to be a good pCO2 proxy that is consistent with pCO2 records obtained using biomarker d13C data from Atlantic black shales and leaf stomata data from WIS sections. Boreal chalk d18O records show sea surface temperature (SST) changes that closely follow the Delta13C pCO2 proxy and confirm TEX86 results from deep ocean sites. Rising pCO2 and SST during the Late Cenomanian is attributed to volcanic degassing; pCO2 and SST maxima occurred at the onset of black shale deposition, followed by falling pCO2 and cooling due to carbon sequestration by marine organic productivity and preservation, and increased silicate weathering. A marked pCO2 minimum (~25% fall) occurred with a SST minimum (Plenus Cold Event) showing >4°C of cooling in ~40 kyr. Renewed increases in pCO2, SST, and d13C during latest Cenomanian black shale deposition suggest that a continuing volcanogenic CO2 flux overrode further drawdown effects. Maximum pCO2 and SST followed the end of OAE2, associated with a falling nutrient supply during the Early Turonian eustatic highstand.

Age models for ODP Leg 172 sites

Ten ODP sites drilled in a depth transect (2164-4775 m water depth) during Leg 172 recovered high-deposition rate (>20 cm/kyr) sedimentary sections from sediment drifts in the western North Atlantic. For each site an age model covering the past 0.8-0.9 Ma has been developed. The time scales have a resolution of 10-20 kyr and are derived by tuning variations of estimated carbonate content to the orbital parameters precession and obliquity. Based on the similarity in the signature of proxy records and the spectral character of the time series, the sites are divided into two groups: precession cycles are better developed in carbonate records from a group of shallow sites (2164-2975 m water depth, Sites 1055-1058) while the deeper sites (2995-4775 m water depth, Sites 1060-1063) are characterized by higher spectral density in the obliquity band. The resulting time scales show excellent coherence with other dated carbonate and isotope records from low latitudes. Besides the typical Milankovitch cyclicity significant variance of the resulting carbonate time series is concentrated at millennial-scale changes with periods of about 12, 6, 4, 2.5, and 1.5 kyr. Comparisons of carbonate records from the Blake Bahama Outer Ridge and the Bermuda Rise reveal a remarkable similarity in the time and frequency domain indicating a basin-wide uniform sedimentation pattern during the last 0.9 Ma.

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.

Neodymium isotopes from fossil fish teeth and detrital sediment, and oxygen isotopes from benthic forams from ODP Hole 119-738 B

Ephemeral polar glaciations during the middle-to-late Eocene (48-34 Ma) have been proposed based on far-field ice volume proxy records and near-field glacigenic sediments, although the scale, timing, and duration of these events are poorly constrained. Here we confirm the existence of a transient cool event within a new high-resolution benthic foraminiferal d18O record at Ocean Drilling Program (ODP) Site 738 (Kerguelen Plateau; Southern Ocean). This event, named the Priabonian oxygen isotope maximum (PrOM) Event, lasted ~140 kyr and is tentatively placed within magnetochron C17n.1n (~37.3 Ma) based on the correlation to ODP Site 689 (Maud Rise, Southern Ocean). A contemporaneous change in the provenance of sediments delivered to the Kerguelen Plateau occurs at the study site, determined from the <63 µm fraction of decarbonated and reductively leached sediment samples. Changes in the mixture of bottom waters, based on fossil fish tooth epsilon-Nd, were less pronounced and slower relative to the benthic d18O and terrigenous epsilon-Nd changes. Terrigenous sediment epsilon-Nd values rapidly shifted to less radiogenic signatures at the onset of the PrOM Event, indicating an abrupt change in provenance favoring ancient sources such as the Paleoproterozoic East Antarctic craton. Bottom water epsilon-Nd reached a minimum value during the PrOM Event, although the shift begins much earlier than the terrigenous epsilon-Nd excursion. The origin of the abrupt change in terrigenous sediment provenance is compatible with a change in Antarctic terrigenous sediment flux and/or source as opposed to a reorganization of ocean currents. A change in terrigenous flux and/or source of Antarctic sediments during the oxygen isotope maximum suggests a combination of cooling and ice growth in East Antarctica during the early late Eocene.

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.