Chronology and petrology in MIS4-2 of DSDP Site 94-609

Ice-rafting evidence for a '1500-year cycle' sparked considerable debate on millennial-scale climate change and the role of solar variability. Here, we reinterpret the last 70,000 years of the subpolar North Atlantic record, focusing on classic DSDP Site 609, in the context of newly available raw data, the latest radiocarbon calibration (Marine09) and ice core chronology (GICC05), and a wider range of statistical methodologies. A ~1500-year oscillation is primarily limited to the short glacial Stage 4, the age of which is derived solely from an ice flow model (ss09sea), subject to uncertainty, and offset most from the original chronology. Results from the most well-dated, younger interval suggest that the original 1500 ± 500 year cycle may actually be an admixture of the ~1000 and ~2000 cycles that are observed within the Holocene at multiple locations. In Holocene sections these variations are coherent with 14C and 10Be estimates of solar variability. Our new results suggest that the '1500-year cycle' may be a transient phenomenon whose origin could be due, for example, to ice sheet boundary conditions for the interval in which it is observed. We therefore question whether it is necessary to invoke such exotic explanations as heterodyne frequencies or combination tones to explain a phenomenon of such fleeting occurrence that is potentially an artifact of arithmetic averaging.

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.

The Antarctic ice core chronology (AICC2012)

Abstract of Bazin et al. (2013): An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120-800 ka. In this framework, new measurements of d18Oatm over Marine Isotope Stage (MIS) 11-12 on EDC and a complete d18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from d18Oatm, dO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one. Abstract of Veres et al. (2013): The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on d15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores.

Age model of firn core DML05C98_32 (B32)

In the framework of the European Project for Ice Coring in Antarctica (EPICA), a comprehensive glaciological pre-site survey has been carried out on Amundsenisen, Dronning Maud Land, East Antarctica, in the past decade. Within this survey, four intermediate-depth ice cores and 13 snow pits were analyzed for their ionic composition and interpreted with respect to the spatial and temporal variability of volcanic sulphate deposition. The comparison of the non-sea-salt (nss)-sulphate peaks that are related to the well-known eruptions of Pinatubo and Cerro Hudson in AD 1991 revealed sulphate depositions of comparable size (15.8 ± 3.4 kg/km**2) in 11 snow pits. There is a tendency to higher annual concentrations for smaller snow-accumulation rates. The combination of seasonal sodium and annually resolved nss-sulphate records allowed the establishment of a time-scale derived by annual-layer counting over the last 2000 years and thus a detailed chronology of annual volcanic sulphate deposition. Using a robust outlier detection algorithm, 49 volcanic eruptions were identified between AD 165 and 1997. The dating uncertainty is ±3 years between AD 1997 and 1601, around ±5 years between AD 1601 and 1257, and increasing to ±24 years at AD 165, improving the accuracy of the volcanic chronology during the penultimate millennium considerably.

Chronology and petrology of sediment cor DAPC2

Evidence from paleoclimatic archives suggests that Earth's climate experienced rapid temperature changes associated with pronounced interhemispheric asymmetry during the last glacial period. Explanations for these climate excursions have converged on nonlinear interactions between ice sheets and the ocean's thermohaline circulation, but the driving mechanism remains to be identified. Here we use multidecadal marine records of faunal, oxygen isotope, and sediment proxies from the northeast Atlantic proximal to the western margins of the last glacial British Ice Sheet (BIS) to document the coupling between ice sheet dynamics, ocean circulation, and insolation changes. The core data reveal successions of short-lived (80-100 years), high-amplitude ice-rafted debris (IRD) events that were initiated up to 2000 years before the deposition of detrital carbonate during Heinrich events (HE) 1 and 2. Progressive disintegration of the BIS 19-16 kyr before present (B.P.) occurred in response to abrupt ocean-climate warmings that impinged on the northeast Atlantic during the early deglaciation. Peak IRD deposition recurs at 180-220 year intervals plausibly involving repeated breakup of glacial tidewater margins and fringing marine ice shelves. The early deglaciation culminated in a major meltwater pulse at ~16.3 kyr B.P. followed by another discharge associated with HE1 some 300 years after. We conclude that temperature changes related to external forcing and marine heat transport caused a rapid response of the BIS and possibly other margins of the Eurasian Ice Sheet. Massive but short-lived meltwater surges influenced the Atlantic meridional overturning circulation thereby contributing to North Atlantic climate variability and bipolar climatic asymmetry.

Age determination of late Quaternary sediment cores from the South Tasman Rise

We use sediment cores from the South Tasman Rise (STR) to reconstruct deep- water circulation in the southwest Pacific sector of the Southern Ocean. Sediment cores MD972106 (45° 09' S, 146° 17' E, 3310 m water depth) and GC34 (45° 06' S, 147° 45' E, 4002 m water depth) preserve records covering the last 160 kyr, with chronology controlled by calibrated accelerator mass spectrometry radiocarbon dates and benthic foraminiferal d18O tied to SPECMAP. The STR benthic foraminiferal d13C records provide new d13C values for Southern Ocean deep water spanning the last 160 kyr at sites unlikely to be affected by variations in productivity. The records establish that glacial benthic foraminifera (Cibicidoides spp.) d13C values are lower relative to interglacial values and are comparable to previous glacial benthic d13C records in the Indian and Pacific sectors of the Southern Ocean. Comparisons of the benthic foraminiferal d13C time series at the STR are made with the equatorial Pacific (V19-30 and Site 846) and the equatorial Atlantic (GeoB1115). The similarity of benthic d13C records at the STR to the equatorial Pacific suggest the Southern Ocean deep-water mass closely tracked those of the deep Pacific, and the presence of a d13C gradient between the STR and the equatorial Atlantic suggests there was continual production of northern source deep water over the past 160 kyr.

Age determination and stable isotope record of planktonic foraminifera from sediments of the Baffin Shelf

Stable isotopic values on planktonic foraminifera in a suite of cores from basins across the SE Baffin Shelf are used to extract a record of meltwater events during Termination I deglaciation. Resolution and Hatton basins lie on the SE Baffin Shelf at water depths > 500 m, seaward of major conduits for ice drainage from the eastern sector of the Laurentide Ice Sheet (LIS). Accelerator mass spectrometry 14C dates are used to constrain our chronology of events in ten cores. In Resolution Basin, three cores have 14C AMS dates on foraminifera of > 20 ka at their bases; whereas Hatton Basin cores terminate in sediments < 13 kyr. Sedimentation rates varied between 0.1 to 4.5 m/ka. Stable oxygen and carbon isotopic ratios were obtained on 146 samples of the planktonic foraminifera Neogloboquadrina pachyderma (Ehrenberg) sinistral, from seven of the ten cores. No evidence was found to indicate that test morphology or size affected delta18O. Between 7 and 13.5 ka the surface water on the shelf was on average 1 per mil lower than the open ocean signal. Significant temporal variations were found in both delta18O and delta13C. Evidence for significant low delta18O events occurred between 13 and 8 ka. The delta13C record from the planktonic foraminifera suggests a threefold division of events between 13 and 7 ka, with positive values between 10.8 and 13.0 ka, negative values between 9 and 10.8 ka, and positive values from 7 to 9 ka. The delta18O data suggest the presence of meltwater on the shelf some 3,000 years prior to the first late glacial dates on terrestrial deglaciation (at circa 10.4 ka). "Hudson Strait must be the real key to the importance of the calving process during deglaciation, because it is potentially the largest marine outlet for the Laurentide Ice Sheet and because it leads into the very center of the ice sheet.....the rates of calving through Hudson Strait during the period of initial ?18O rise unfortunately are unknown." W. F. Ruddiman (1987, p. 151)

Age models and stable isotope ratios of SO12 and SO35 sediment cores

Based on organic carbon accumulation rates, nine time slices of oceanic export paleoproductivity (Pnew) are presented which depict the variability of Pnew on a global scale through the last 30,000 years and document that the basic distribution patterns did not change through glacial and interglacial times. However, the glacial ocean shows an increased contrast of high- versus low-productivity zones. d13C values of near-surface-dwelling planktonic foraminifera Globigerinoides ruber suggest that the same contrast applies to the glacial nutrient inventories of the ambient surface waters, with a significant glacial transfer of PO4 from low- to high-productivity zones. In this way, glacial Pnew increased by a global average of about 2-4 Gt C/yr and led, via an enhanced CaCO3 dissolution and alkalinity in the deep ocean, to a significant extraction of CO2 from the surface water and the atrnosphere.

Age determinations and analyses of organic matter of sediments from the Argentine Basin

The first radiocarbon chronology for sediments of the Argentine basin has been determined using accelerator mass spectrometer (AMS) analyses of 54 total organic carbon samples from four box and two piston cores collected from the downstream and upstream sides of two central Argentine Basin mudwaves. Throughout the Holocene, sediment from the geomorphically defined upstream side of each wave accumulated at rates of 30 to 105 cm/1000 years. Sediments from the downstream side of each wave accumulated at rates of 2 to 10 cm/1000 years in the late and early Holocene, while the mid Holocene is characterized by sedimentation rates less than 1.0 cm/1000 years. During the mid-Holocene, increased aridity reduced chemical weathering and the flow of the rivers draining to the continental shelf, causing a concomitant decrease in fine-grained terrigenous input to the basin as evidenced by decreased sedimentation rates, lower N/C ratios, and depleted delta13Corg values. It is estimated that all of the organic carbon deposited in the central basin during the mid-Holocene was of a marine origin. During the late and early Holocene, however, approximately 35% of the organic carbon deposited was of terrestrial origin. Bottom water flow speeds in the late Holocene were estimated using a lee-wave model and found to average 14 cm/s. This estimate is comparable to 10 cm/s mean and 15-20 cm/s maximum flow speeds measured by current meters deployed within the basin. Flow speeds in the Argentine Basin were 10% higher than today from 8000 to 2000 B.P., and are consistent with a general invigoration of thermohaline circulation that began between 9000 and 8000 B.P. It is proposed that the introduction of warm, salty Indian Ocean water into the northern North Atlantic at 9000 B.P. was the mechanism that provided the excess salt needed to stabilize the North Atlantic Deep Water thermohaline circulation system in its present mode.