(Table 1) Color, Fe-C-S chemistry and carbonate content from ODP Sites 172-1062 and 172-1063

Reflectance spectra collected during ODP Leg 172 were used in concert with solid phase iron chemistry, carbonate content, and organic carbon content measurements to evaluate the agents responsible for setting the color in sediments. Factor analysis has proved a valuable and rapid technique to detect the local and regional primary factors that influence sediment color. On the western North Atlantic drifts, sediment color is the result of primary mineralogy as well as diagenetic changes. Sediment lightness is controlled by the carbonate content while the hue is primarily due to the presence of hematite and Fe2+/Fe3+ changes in clay minerals. Hematite, most likely derived from the Permo-Carboniferous red beds of the Canadian Maritimes, is differentially preserved at various sites due to differences in reductive diagenesis and dilution by other sedimentary components. Various intensities for diagenesis result from changes in organic carbon content, sedimentation rates, and H2S production via anaerobic methane oxidation. Iron monosulfides occur extensively at all high sedimentation sites especially in glacial periods suggesting increased high terrigenous flux and/or increased reactive iron flux in glacials.

Uranium, Thorium and Protactinium isotopes of ODP Site 172-1063

Ocean circulation may have undergone reductions and reinvigorations in the past closely tied to regional climate changes. Measurements of 231Pa/230Th ratios in a sediment core from the Bermuda Rise have been interpreted as evidence that the Atlantic Meridional Overturning Circulation (AMOC) was weakened or completely eliminated during a period of catastrophic iceberg discharges (Heinrich-Event 1, H1). Here we present new data from the Bermuda Rise that show further 231Pa/230Th peaks during Heinrich-2 (H2) and Heinrich-3 (H3). Additionally, a tight correlation between diatom abundances (biogenic silica) and 231Pa/230Th is discovered in this core. Our results redirect the interpretation of 231Pa/230Th from the Bermuda Rise as a proxy for ocean circulation towards a proxy that reacts highly sensitive to changes of particle composition and water mass properties.

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

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.

(Table S1) Tie-points used in age model construction for ODP Leg 172 sites

ix Ocean Drilling Program (ODP) sites, in the Northwest Atlantic have been used to investigate kinematic and chemical changes in the "Western Boundary Undercurrent" (WBUC) during the development of full glacial conditions across the Marine Isotope Stage 5a/4 boundary (~70,000 years ago). Sortable silt mean grain size(sort s) measurements are employed to examine changes in near bottom flow speeds, together with carbon isotopes measured in benthic foraminifera and % planktic foraminiferal fragmentation as proxies for changes in water-mass chemistry. A depth transect of cores, spanning 1.8-4.6 km depth, allows changes in both the strength and depth of the WBUC to be constrained across millennial scale events. Sort s measurements reveal that the flow speed structure of the WBUC during warm intervals ("interstadials") was comparable to modern (Holocene) conditions. However, significant differences are observed during cold intervals, with higher relative flow speeds inferred for the shallow component of the WBUC (~2 km depth) during all cold "stadial" intervals (including Heinrich Stadial 6), and a substantial weakening of the deep component (~3-4 km) during full glacial conditions. Our results therefore reveal that the onset of full glacial conditions was associated with a regime shift to a shallower mode of circulation (involving Glacial North Atlantic Intermediate Water) that was quantitatively distinct from preceding cold stadial events. Furthermore, our chemical proxy data show that the physical response of the WBUC during the last glacial inception was probably coupled to basin-wide changes in the water-mass composition of the deep Northwest Atlantic.

Revision of the early-middle Pleistocene calcareous nannofossil biochronology

The extant nannofossil biostratigraphic and biochronologic framework for the early-middle Pleistocene time interval has been tested through the micropaleontological analysis of globally distributed high-quality low- to mid-latitude deep-sea successions. The quantitative temporal distribution patterns of relative abundances of selected taxa were reconstructed in critical intervals, and the following biohorizons were defined: first occurrence of medium-sized Gephyrocapsa spp. (bmG); last occurrence of Calcidiscus macintyrei (tCm); first occurrence of large Gephyrocapsa spp. (blG); last occurrence of large Gephyrocapsa spp. (tlG); first occurrence of Reticulofenestra asanoi (bRa); re-entrance of medium-sized Gephyrocapsa spp. (reemG) and last occurrence of Reticulofenestra asanoi (tRa). The detailed patterns of abundance change at these biohorizons were used to generate a detailed biostratigraphy, and the biostratigraphic data were transformed into a precise biochronology by means of correlation to isotope stratigraphies and astronomical timescales. The degree of isochrony or diachrony of the biohorizons was evaluated. Biohorizons tlG and tRa are isochronous occurring close to marine isotope stages (MIS)55 and MIS 22, respectively, and bmG and blG are slightly diachronous on the order of 30-40 kyr, whereas biohorizons tCm, reemG and bRa are confirmed as diachronous on the order of 100, 80 and 60 kyr, respectively. Some of the events are clearly controlled by environmental conditions, e.g. the last occurrence of R. asanoi, related to significant environmental changes associated with the first large-amplitude glaciation of the late Quaternary, MIS 22.

Nd isotopes of sediments, foraminifera and fish debris of the western North Atlantic

Understanding changes in ocean circulation during the last deglaciation is crucial to unraveling the dynamics of glacial-interglacial and millennial climate shifts. We used neodymium isotope measurements on postdepositional iron-manganese oxide coatings precipitated on planktonic foraminifera to reconstruct changes in the bottom water source of the deep western North Atlantic at the Bermuda Rise. Comparison of our deep water source record with overturning strength proxies shows that both the deep water mass source and the overturning rate shifted rapidly and synchronously during the last deglacial transition. In contrast, any freshwater perturbation caused by Heinrich event 1 could have only affected shallow overturning. These findings show how changes in upper-ocean overturning associated with millennial-scale events differ from those associated with whole-ocean deglacial climate events.

Stable isotope ratios and abundances of selected foraminifera taxa from ODP Leg 172 sites

This data report describes the results of post-Leg 172 sampling of Sites 1054, 1055, and 1063 for two purposes: to investigate the climatic significance of red-colored intervals in the hemipelagic sediments cored during Leg 172 and to better understand the stratigraphy and chronology of Carolina Slope Sites 1054 and 1055. Gravity cores collected from the Carolina Slope on site survey cruise Knorr 140/2 show very high rates of sedimentation during the Holocene and lower rates during the last glacial maximum (LGM). Because of the high rates, many of the sediments in the recovered cores never reached the LGM. In other cores, it is possible that deglacial oscillations have been mistaken for the LGM. Although radiocarbon dating could solve that problem, some of the gravity cores are at or very close to the Ocean Drilling Program (ODP) sites, and it is useful to compare the isotope stratigraphies among them before proceeding with dating. Furthermore, some of the site survey cores have red-colored intervals and others do not, even though there is some indication they are time equivalent. Either the stratigraphy is wrong, diagenesis has affected the color of the sediment, or red sediment is carried to some sites but not to others that differ in depth by only a few hundred meters.

Mineralogy and major element oxide composition of sediments from ODP Hole 172-1063D

Dansgaard-Oeschger (D-O) cycles in sediment at Site 1063 are characterized by distinct fluctuations in physical properties. Stadials are marked by low bulk density and interstadials by high bulk density. Compressional (P-)wave velocity is in phase with bulk density over some but not all depth intervals. Four of the D-O cycles straddling the oxygen isotope Stage 4/5 boundary have been studied in detail to understand the origin of the physical properties changes. Sediment on the Bermuda Rise is comprised of three main components: calcite, aluminosilicate minerals, and biogenic silica. Calcite concentrations vary from 1% to 43% of bulk sediment and are highest during interstadials. Aluminosilicate concentrations vary from 52% to 92% of bulk sediment and are highest during stadials. The major element ratios Al2O3/TiO2 and K2O/Al2O3 show increases across bulk density cycles, suggesting a change in the composition of aluminosilicates. This interpretation is supported by mineralogical analyses, which show a subtle change in clay composition. Biogenic silica concentrations vary from 0% to 23% of bulk sediment and are also highest during stadials. However, the abundance of silica varies significantly from one D-O cycle to another. Silt and fine sand abundance also increase during the first of the four stadials. This coarsening of sediment coincides with the increase in biogenic silica. The low grain density and high porosity associated with biogenic silica result in intervals of low bulk-sediment density. The abundance of biogenic silica closely matches P-wave velocity, suggesting that silica imparts a greater rigidity to the sediment.

(Table T1) Calcium carbonate content of sediments from ODP Leg 172 sites

Visual-domain diffuse reflectance data collected aboard the JOIDES Resolution with the Minolta spectrometer CM-2002 during Ocean Drilling Program Leg 172 have been used to estimate successfully the carbonate content of sediments. Calibration equations were developed for each site and for each lithostratigraphic unit (or subunit at Site 1063) using multiple linear regression on raw as well as pretreated reflectance spectra (i.e., first-order derivation and squaring of raw reflectance spectra) for a total of 4141 direct carbonate measurements. The root-mean-square errors of 4% to 7% are within the range of previous estimates using diffuse reflectance data and are acceptable for the general extensive range of carbonate contents (i.e., 0-70 wt%) that characterize sedimentation at Leg 172 sites.