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.

Age determination and sedimentology of sediment core Lz1101 from the Melles Lake

A lacustrine sediment core from Store Koldewey, northeast Greenland, was biogeochemically, biologically and sedimentologically investigated in order to reconstruct long- and short-term climatic and environmental variability. The chronology of the uppermost 189 cm of the record is based on ten 14C AMS age determinations of aquatic mosses. The record covers almost the entire Holocene and revealed changes on multidecadal to centennial scales. Dating of the oldest mosses shows that lacustrine biogenic productivity already began at around 11 cal. kyr BP. This age pre-dates the onset of biogenic productivity in other lakes on Store Koldewey by about 2 kyr. In spite of the early onset of biogenic production organic matter accumulation remained low and minerogenic sedimentation dominated. At about 9.5 cal. kyr BP moss, sulphur, organic carbon and biogenic silica content started to increase, indicating that the environment stabilized and the biogenic production in the lake adjusted to more preferable conditions. Subsequently, the biogenic productivity experienced repeated changes and varied both on long- and short-term scales. The long-term trend shows a maximum during the early Holocene thus responding to increased temperatures during the Holocene Thermal Maximum. Superimposed on the long-term trend, biogenic productivity also experienced repeated short-term fluctuations that match partly the NGRIP temperatures. The most pronounced decrease of biogenic productivity occurred at around 8.2 cal. kyr BP. Perennial lake ice coverage resulting from low temperatures is supposed to have caused decreased lacustrine biogenic productivity. From the middle Holocene to the present repeated decreases of productivity occurred that could be related to periods with severe sea-ice conditions of the East Greenland Current. Besides the dependence on air temperature it therefore demonstrates the sensitivity of lacustrine biogenic productivity in coastal high arctic areas to short-term cold spells that are mediated by the currents emanating from the Arctic Ocean. However, the data also emphasize the difficulties associated with the interpretation of lacustrine records.

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.