Additional CaCO3 measurements from IODP Site 306-U1313

We present further %CaCO3 data from Site U1313 across the Pliocene-Pleistocene intensification of Northern Hemisphere glaciation. This data was measured on the U1313 secondary splice. We also present tie points between the primary and secondary splice for this interval based on graphical tuning of L* (sediment lightness).

CaCO3, foraminifer fragmentation, benthic carbon isotopes, coarse lithic counts and radiogenic isotopes (Sr, Nd, Pb) of the Pliocene and earliest Pleistocene terrigenous component (3.3-2.4 Ma) of IODP Site 306-U1313

We present Plio-Pleistocene records of sediment color, %CaCO3, foraminifer fragmentation, benthic carbon isotopes (d13C) and radiogenic isotopes (Sr, Nd, Pb) of the terrigenous component from IODP Site U1313, a reoccupation of benchmark subtropical North Atlantic Ocean DSDP Site 607. We show that (inter)glacial cycles in sediment color and %CaCO3 pre-date major northern hemisphere glaciation and are unambiguously and consistently correlated to benthic oxygen isotopes back to 3.3 million years ago (Ma) and intermittently so probably back to the Miocene/Pliocene boundary. We show these lithological cycles to be driven by enhanced glacial fluxes of terrigenous material (eolian dust), not carbonate dissolution (the classic interpretation). Our radiogenic isotope data indicate a North American source for this dust (~3.3-2.4 Ma) in keeping with the interpreted source of terrestrial plant wax-derived biomarkers deposited at Site U1313. Yet our data indicate a mid latitude provenance regardless of (inter)glacial state, a finding that is inconsistent with the biomarker-inferred importance of glaciogenic mechanisms of dust production and transport. Moreover, we find that the relation between the biomarker and lithogenic components of dust accumulation is distinctly non-linear. Both records show a jump in glacial rates of accumulation from Marine Isotope Stage, MIS, G6 (2.72 Ma) onwards but the amplitude of this signal is about 3-8 times greater for biomarkers than for dust and particularly extreme during MIS 100 (2.52 Ma). We conclude that North America shifted abruptly to a distinctly more arid glacial regime from MIS G6, but major shifts in glacial North American vegetation biomes and regional wind fields (exacerbated by the growth of a large Laurentide Ice Sheet during MIS 100) likely explain amplification of this signal in the biomarker records. Our findings are consistent with wetter-than-modern reconstructions of North American continental climate under the warm high CO2 conditions of the Early Pliocene but contrast with most model predictions for the response of the hydrological cycle to anthropogenic warming over the coming 50 years (poleward expansion of the subtropical dry zones).

Sea surface temperature estimation for the last 70 ka from IODP Site 306-U1313

During the six Heinrich Events of the last 70 ka episodic calving from the circum-Atlantic ice sheets released large numbers of icebergs into the North Atlantic. These icebergs and associated melt-water flux are hypothesized to have led to a shutdown of Atlantic Meridional Overturning Circulation (AMOC) and severe cooling in large parts of the Northern Hemisphere. However, due to the limited availability of high-resolution records the magnitude sea surface temperature (SST) changes related to the impact of Heinrich Events on the mid-latitude North Atlantic is poorly constrained. Here we present a record of UK37'-based SSTs derived from sediments of Integrated Ocean Drilling Project (IODP) Site U1313, located at the southern end of the ice-rafted debris (IRD)-belt in the mid-latitude North Atlantic (41°N). We demonstrate that all six Heinrich Events are associated with a rapid warming of surface waters by 2 to 4°C in a few thousand years. The presence of IRD leaves no doubt about the simultaneous timing and correlation between rapid surface water warming and Heinrich Events. We argue that this warming in the mid-latitude North Atlantic is related to a northward expansion of the subtropical gyre during Heinrich Events. As a wide-range of studies demonstrated that in the central IRD-belt Heinrich Events are associated with low SSTs, these results thus identify an anti-phased (seesaw) pattern in SSTs during Heinrich Events between the mid-latitude (warm) and northern North Atlantic (cold). This highlights the complex response of surface water characteristics in the North Atlantic to Heinrich Events that is poorly reproduced by fresh water hosing experiments and challenges the widely accepted view that within the IRD-belt of the North Atlantic Heinrich Events coincide with periods of low SSTs.

Sea surface temperature estimation for MIS 16 of IODP Site 306-U1313

Hudson Strait (HS) Heinrich Events, ice-rafting events in the North Atlantic originating from the Laurentide ice sheet (LIS), are among the most dramatic examples of millennial-scale climate variability and have a large influence on global climate. However, it is debated as to whether the occurrence of HS Heinrich Events in the (eastern) North Atlantic in the geological record depends on greater ice discharge, or simply from the longer survival of icebergs in cold waters. Using sediments from Integrated Ocean Drilling Program (IODP) Site U1313 in the North Atlantic spanning the period between 960 and 320 ka, we show that sea surface temperatures (SSTs) did not control the first occurrence of HS Heinrich(-like) Events in the sedimentary record. Using mineralogy and organic geochemistry to determine the characteristics of ice-rafting debris (IRD), we detect the first HS Heinrich(-like) Event in our record around 643 ka (Marine Isotope Stage (MIS) 16), which is similar as previously reported for Site U1308. However, the accompanying high-resolution alkenone-based SST record demonstrates that the first HS Heinrich(-like) Event did not coincide with low SSTs. Thus, the HS Heinrich(-like) Events do indicate enhanced ice discharge from the LIS at the end of the Mid-Pleistocene Transition, not simply the survivability of icebergs due to cold conditions in the North Atlantic.

Biomarker and XRD results for the last 3.4 Ma from IODP Site 306-U1313

Various types of abrupt/millennial-scale climate variability such as Dansgaard/Oeschger and Heinrich Events characterized the last glacial period. Over the last decade, a number of studies demonstrated that such millennial-scale climate variability was not limited to the last glacial but inherent to Quaternary climate. Here we review the occurrence and origin of millennial ice-rafting events in the North Atlantic during the late Pliocene and Pleistocene (last 3.4 Ma) with a special focus on North Atlantic Hudson Strait (HS) Heinrich(-like) Events. Besides a clear biomarker signature, we show that Heinrich Layers 5, 4, 2, and 1 in marine sediment cores from across the North Atlantic all bear the organic geochemical fingerprint of the Hudson area. Using this framework and combining previously published results, detailed investigations into the organic and inorganic chemistry of ice-rafted debris (IRD) found across the North Atlantic demonstrate that prior to MIS 16 (~ 650 ka) IRD in the North Atlantic did not originate from the Hudson area of northern Canada. The signature of this early IRD is distinctly different compared to that of HS Heinrich Layers. Rather ice-rafting events during the late Pliocene and early Pleistocene predominantly emanated from the calving of the Greenland and Fennoscandian ice sheets and possibly minor contributions from local ice streams from the North American and British ice sheets. Compared to North Atlantic HS Heinrich Events, these early Pleistocene IRD-events had a limited impact on surface water characteristics in the North Atlantic. North Atlantic HS Heinrich(-like) Events first occurred during MIS 16. At the same time, the dominant frequency in silicate-rich IRD accumulation shifted from the obliquity (41-ka) to a 100-ka frequency across the North Atlantic. Iceberg survivability or a change in iceberg trajectory likely did not control this change in IRD-regime. These results lend further support for the existing hypothesis that an increase in size (thickness) of the Laurentide ice sheet controls the occurrence of North Atlantic HS Heinrich Events, favoring an internal dynamic mechanism for their occurrence.

Stable carbon and oxygen isotope record od IODP Hole 306-U1313

The Mid-Pleistocene transition (MPT) was the time when quasi-periodic (? 100 kyr), high-amplitude glacial variability developed in the absence of any significant change in the character of orbital forcing, leading to the establishment of the characteristic pattern of late Pleistocene climate variability. It has long been known that the interval around 900 ka stands out as a critical point of the MPT, when major glaciations started occurring most notably in the northern hemisphere. Here we examine the record of climatic conditions during this significant interval, using high-resolution stable isotope records from benthic and planktonic foraminifera from a sediment core in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313). We have considered the time interval from late in Marine Isotope Stage (MIS) 23 to MIS 20 (910 to 790 ka). Our data indicate that interglacial MIS 21 was a climatically unstable period and was broken into four interstadial periods, which have been identified and correlated across the North Atlantic region. These extra peaks tend to contradict previous studies that interpreted the MIS 21 variability as consisting essentially of a linear response to cyclical changes in orbital parameters. Cooling events in the surface record during MIS 21 were associated with low benthic carbon isotope excursions, suggesting a coupling between surface temperature changes and the strength of the Atlantic meridional overturning circulation. Time series analysis performed on the whole interval indicates that benthic and planktonic oxygen isotopes have significant concentrations of spectral power centered on periods of 10.7 kyr and 6 kyr, which is in agreement with the second and forth harmonic of precession. The excellent correspondence between the foraminifera d18O records and insolation variations at the Equator in March and September suggests that a mechanism related to low-latitude precession variations, advected to the high latitudes by tropical convective processes, might have generated such a response. This scenario accounts for the presence of oscillations at frequencies equal to precession harmonics at Site U1313, as well as the occurrence of higher amplitude oscillations between the MIS22/21 transition and most of MIS 21, times of enhanced insolation variability.

Stable carbon and oxygen isotope record of IODP Site 306-U1313, MIS 23-19

Since the seminal work by Hays et al. (1976), a plethora of studies has demonstrated a correlation between orbital variations and climatic change. However, information on how changes in orbital boundary conditions affected the frequency and amplitude of millennial-scale climate variability is still fragmentary. The Marine Isotope Stage (MIS) 19, an interglacial centred at around 785 ka, provides an opportunity to pursue this question and test the hypothesis that the long-term processes set up the boundary conditions within which the short-term processes operate. Similarly to the current interglacial, MIS 19 is characterised by a minimum of the 400-kyr eccentricity cycle, subdued amplitude of precessional changes, and small amplitude variations in insolation. Here we examine the record of climatic conditions during MIS 19 using high-resolution stable isotope records from benthic and planktonic foraminifera from a sedimentary sequence in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313) in order to assess the stability and duration of this interglacial, and evaluate the climate system's response in the millennial band to known orbitally induced insolation changes. Benthic and planktonic foraminiferal d18O values indicate relatively stable conditions during the peak warmth of MIS 19, but sea-surface and deep-water reconstructions start diverging during the transition towards the glacial MIS 18, when large, cold excursions disrupt the surface waters whereas low amplitude millennial scale fluctuations persist in the deep waters as recorded by the oxygen isotope signal. The glacial inception occurred at ~779 ka, in agreement with an increased abundance of tetra-unsaturated alkenones, reflecting the influence of icebergs and associated meltwater pulses and high-latitude waters at the study site. After having combined the new results with previous data from the same site, and using a variety of time series analysis techniques, we evaluate the evolution of millennial climate variability in response to changing orbital boundary conditions during the Early-Middle Pleistocene. Suborbital variability in both surface- and deep-water records is mainly concentrated at a period of ~11 kyr and, additionally, at ~5.8 and ~3.9 kyr in the deep ocean; these periods are equal to harmonics of precession band oscillations. The fact that the response at the 11 kyr period increased over the same interval during which the amplitude of the response to the precessional cycle increased supports the notion that most of the variance in the 11 kyr band in the sedimentary record is nonlinearly transferred from precession band oscillations. Considering that these periodicities are important features in the equatorial and intertropical insolation, these observations are in line with the view that the low-latitude regions play an important role in the response of the climate system to the astronomical forcing. We conclude that the effect of the orbitally induced insolation is of fundamental importance in regulating the timing and amplitude of millennial scale climate variability.