Phase relationships between orbital forcing and the composition of air trapped in Antarctic ice cores

Orbital tuning is central for ice core chronologies beyond annual layer counting, available back to 60 ka (i.e. thousands of years before 1950) for Greenland ice cores. While several complementary orbital tuning tools have recently been developed using δ¹⁸Oatm, δO₂⁄N₂ and air content with different orbital targets, quantifying their uncertainties remains a challenge. Indeed, the exact processes linking variations of these parameters, measured in the air trapped in ice, to their orbital targets are not yet fully understood. Here, we provide new series of δO₂∕N₂ and δ¹⁸Oatm data encompassing Marine Isotopic Stage (MIS) 5 (between 100 and 160 ka) and the oldest part (340–800 ka) of the East Antarctic EPICA Dome C (EDC) ice core. For the first time, the measurements over MIS 5 allow an inter-comparison of δO₂∕N₂ and δ¹⁸Oatm records from three East Antarctic ice core sites (EDC, Vostok and Dome F). This comparison highlights some site-specific δO₂∕N₂ variations. Such an observation, the evidence of a 100 ka periodicity in the δO₂∕N₂ signal and the difficulty to identify extrema and mid-slopes in δO2∕N2 increase the uncertainty associated with the use of δO₂∕N₂ as an orbital tuning tool, now calculated to be 3–4 ka. When combining records of δ¹⁸Oatm and δO₂∕N₂ from Vostok and EDC, we find a loss of orbital signature for these two parameters during periods of minimum eccentricity (∼ 400 ka, ∼ 720–800 ka). Our data set reveals a time-varying offset between δO₂∕N₂ and δ¹⁸Oatm records over the last 800 ka that we interpret as variations in the lagged response of δ¹⁸Oatm to precession. The largest offsets are identified during Terminations II, MIS 8 and MIS 16, corresponding to periods of destabilization of the Northern polar ice sheets. We therefore suggest that the occurrence of Heinrich–like events influences the response of δ¹⁸Oatm to precession.

Benthic and planktic foraminiferal stable isotope record and coccolith abundance spanning the last 500,000 years of ODP Hole 121-758C

The Indian Summer Monsoon (ISM) is an inter-hemispheric and highly variable ocean-atmosphere-land interaction that directly affects the densely populated Indian subcontinent. Here, we present new records of palaeoceanographic variability that span the last 500,000 years from the eastern equatorial Indian Ocean, a relatively under-sampled area of ISM influence. We have generated carbon and oxygen stable isotope records from three foraminiferal species from Ocean Drilling Program Site 758 (5°N, 90°E) to investigate the oceanographic history of this region. We interpret our resultant Dd18O (surface-thermocline) record of upper water-column stratification in the context of past ISM variability, and compare orbital phase relationships in our Site 758 data to other climate and monsoon proxies in the region. Results suggest that upper water-column stratification at Site 758, which is dominated by variance at precession and half-precession frequencies (23, 19 and 11 ka), is forced by both local (5°N) insolation and ISM winds. In the precession (23 ka) band, stratification minima at Site 758 lag northern hemisphere summer insolation maxima (precession minima) by 9 ka, which is consistent with Arabian Sea ISM phase estimates and suggests a common wind forcing in both regions. This phase implicates a strong sensitivity to both ice volume and southern hemisphere insolation forcing via latent heat export from the southern subtropical Indian Ocean. Additionally, we find evidence of possible overprinting of millennial-scale events during glacial terminations in our stratification record, which suggests an influence of remote abrupt climate events on ISM dynamics.

Pliocene-Pleistocene stack of globally distributed benthic stable oxygen isotope records

We present a 5.3-Myr stack (the ''LR04'' stack) of benthic d18O records from 57 globally distributed sites aligned by an automated graphic correlation algorithm. This is the first benthic delta18O stack composed of more than three records to extend beyond 850 ka, and we use its improved signal quality to identify 24 new marine isotope stages in the early Pliocene. We also present a new LR04 age model for the Pliocene-Pleistocene derived from tuning the delta18O stack to a simple ice model based on 21 June insolation at 65 N. Stacked sedimentation rates provide additional age model constraints to prevent overtuning. Despite a conservative tuning strategy, the LR04 benthic stack exhibits significant coherency with insolation in the obliquity band throughout the entire 5.3 Myr and in the precession band for more than half of the record. The LR04 stack contains significantly more variance in benthic delta18O than previously published stacks of the late Pleistocene as the result of higher resolution records, a better alignment technique, and a greater percentage of records from the Atlantic. Finally, the relative phases of the stack's 41- and 23-kyr components suggest that the precession component of delta18O from 2.7-1.6 Ma is primarily a deep-water temperature signal and that the phase of d18O precession response changed suddenly at 1.6 Ma.

Parameters of the Earths orbit for the last 5 Million years in 1 kyr resolution

New values for the astronomical parameters of the Earth's orbit and rotation (eccentricity, obliquity and precession) are proposed for paleoclimatic research related to the Late Miocene, the Pliocene and the Quaternary. They have been obtained from a numerical solution of the Lagrangian system of the planetary point masses and from an analytical solution of the Poisson equations of the Earth-Moon system. The analytical expansion developed in this paper allows the direct determination of the main frequencies with their phase and amplitude. Numerical and analytical comparisons with the former astronomical solution BER78 are performed so that the accuracy and the interval of time over which the new solution is valid can be estimated. The corresponding insolation values have also been computed and compared to the former ones. This analysis leads to the conclusion that the new values are expected to be reliable over the last 5 Ma in the time domain and at least over the last 10 Ma in the frequency domain.

Stable isotope record during the past 1 Myr of ODP Site 162-980 in the subpolar North Atlantic

Benthic foraminiferal carbon isotope records from a suite of drill sites in the North Atlantic are used to trace variations in the relative strengths of Lower North Atlantic Deep Water (LNADW), Upper North Atlantic Deep Water (UNADW), and Southern Ocean Water (SOW) over the past 1 Myr. During glacial intervals, significant increases in intermediate-to-deep delta13C gradients (commonly reaching >1.2?) are consistent with changes in deep water circulation and associated chemical stratification. Bathymetric delta13C gradients covary with benthic foraminiferal delta18O and covary inversely with Vostok CO2, in agreement with chemical stratification as a driver of atmospheric CO2 changes. Three deep circulation indices based on delta13C show a phasing similar to North Atlantic sea surface temperatures, consistent with a Northern Hemisphere control of NADW/SOW variations. However, lags in the precession band indicate that factors other than deep water circulation control ice volume variations at least in this band.

Age model of sediment core V20-120 (Table 2)

Time series analyses of atmospheric and oceanic variables in a late Pleistocene record from the northwest Pacific show the complex relationship of the response of various segments of the climate system to changes in the earth's orbit. Most variance spectra of time series from this subarctic record contain frequency peaks with periods corresponding to at least one of the major orbital components of eccentricity, obliquity, or precession. Although the radiolarian faunal (water mass) assemblages have prominent spectral peaks with 41,000-year periods which are coherent with obliquity at this frequency, only the Transitional faunal assemblage contains variance focused at a frequency corresponding to the 100,000-year period of eccentricity. Three of these faunal time series also show variance concentrated at a frequency with a 20,000-year period. These three time series are not coherent at a 20,000-year frequency with either of the dominant spectral peaks of precession. They are coherent, however, with variations in the second harmonic of the obliquity cycle. Changes in obliquity apparently affect siliceous faunal abundances in the northwest Pacific region of this high-latitude site more than variations in eccentricity or precession. Maxima in the time series of quartz abundance occur during low values of eccentricity and high glacial ice volume. Because atmospheric winds serve as the major source of supply of quartz to the sediments at this site, these high quartz values reflect increased aridity at the source region. Except for short periods during interglacials, the characteristics of the surface waters in this region of the subarctic Pacific during much of the last 460,000 years were similar to those which exist today in the Sea of Okhotsk. The spectrum of winter sea surface temperature estimates, derived from siliceous microfaunal abundances, contains dominant frequency peaks at periods of 100,000, 41,000, and 23,000 years which are coherent with eccentricity, obliquity, and precession, respectively. Based upon the relationship of the Subarctic Front with the dominance of specific faunal asemblages, the front was positioned south of its present-day location throughout much of the late Pleistocene.

Massic magnetic susceptibility and CaCO3 content at the P. ohmi and T. hugii zones from the Rio Argos X.Ag-1 section, southeastern Spain

An orbital floating time scale of the Hauterivian-Barremian transition (Early Cretaceous) is proposed using high-resolution magnetic susceptibility measurements. Orbital tuning was performed on the Río Argos section (southeast Spain), the candidate for a Global boundary Stratotype Section and Point (GSSP) for the Hauterivian-Barremian transition. Spectral analyses of MS variations, coupled with the frequency ratio method, allow the recognition of precession, obliquity and eccentricity frequency bands. Orbitally-tuned magnetic susceptibility provides minimum durations for ammonite biozones. The durations of well-constrained ammonite zones are assessed at 0.78 myr for Pseudothurmannia ohmi (Late Hauterivian) and 0.57 myr for Taveraidiscus hugii (Early Barremian). These results are consistent with previous estimates from the other reference section (Angles, southeast France) and tend to show that the Río Argos section displays a complete succession for this time interval. They differ significantly from those proposed in the Geologic Time Scale 2008 and may help to improve the next compilation. The Faraoni Oceanic Anoxic Event, a key Early Cretaceous oceanographic perturbation occurring at the P. ohmi/P. catulloi subzone boundary has a duration estimated at 0.10-0.15 myr, which is similar to previous assessments.

(Table 1) Age control points of ODP Site 160-968

The astronomical timescale of the Eastern Mediterranean Plio-Pleistocene builds on tuning of sapropel layers to Northern Hemisphere summer insolation maxima. A 3000-year precession lag has become instrumental in the tuning procedure as radiocarbon dating revealed that the midpoint of the youngest sapropel, S1, in the early Holocene occurred approximately 3000 years after the insolation maximum. The origin of the time lag remains elusive, however, because sapropels are generally linked to maximum African monsoon intensities and transient climate modeling results indicate an in-phase behavior of the African monsoon relative to precession forcing. Here we present new high-resolution records of bulk sediment geochemistry and benthic foraminiferal oxygen isotopes from ODP Site 968 in the Eastern Mediterranean. We show that the 3000-year precession time lag of the sapropel midpoints is consistent with (1) the global marine isotope chronology, (2) maximum (monsoonal) precipitation conditions in the Mediterranean region and China derived from radiometrically dated speleothem records, and (3) maximum atmospheric methane concentrations in Antarctica ice cores. We show that the time lag relates to the occurrence of precession-paced North Atlantic cold events, which systematically delayed the onset of strong boreal summer monsoon intensity. Our findings may also explain a non-stationary behavior of the African monsoon over the past 3 million years due to more frequent and intensive cold events in the Late Pleistocene.

Construction of a precise early Eocene orbital cyclostratigraphy for DSDP Site 80-550 from the Goban Spur, North Atlantic

We constructed a precise early Eocene orbital cyclostratigraphy for DSDP Site 550 (Leg 80, Goban Spur, North Atlantic) utilizing precession related cycles as represented in a high resolution X-Ray Fluorescence based Barium core log. Based on counting of those cycles, we constrain the exact timing of two volcanic ash layers in Site 550 which correlate to ashes +19 and -17 of the Fur Formation in Denmark. The ashes, relative to the onset of the Paleocene/Eocene Thermal Maximum (PETM), are offset by 862 kyr and 672 kyr, respectively. When combined with published absolute ages for ash -17, the absolute age for the onset of the PETM is consistent with astronomically calibrated ages. Using the current absolute age of 28.02 Ma for the Fish Canyon Tuff (FCT) standard for calibrating the absolute age of ash -17 is consistent with tuning option 2 in the astronomically calibrated Paleocene time scale of Westerhold et al. (2008) [Westerhold, T., Röhl, U., Raffi, I., Fornaciari, E., Monechi, S., Reale, V., Bowles, J., and Evans, H.F., 2008, Astronomical calibration of the Paleocene time: Palaeogeography, Palaeoclimatology, Palaeoecology, v. 257, p. 377-403]. Using the recently recalibrated absolute age of 28.201 Ma for the FCT standard is consistent with tuning option 3 in the astronomically calibrated Paleocene time scale. The new results do not support the existence of any additional 405-kyr cycle in the early Paleocene astronomically tuned time scale.

Middle to Late Quaternary sedimentation and rock-magnetic of ODP Site 105-646

We examine rock-magnetic, carbonate, and planktonic foraminiferal fluxes to identify climatically controlled changes of terrigenous and pelagic sedimentation at Ocean Drilling Program (ODP) Site 646 (the Labrador Sea). Terrigenous sediments are brought to the site principally by bottom currents. We use a rock-magnetic parameter sensitive to changes in magnetic mineral grain size, the ratio of anhysteretic susceptibility to low-field magnetic susceptibility (XARM/X), to monitor changes in bottom-current intensity over time, with large values of XARM/X (finer-grained magnetic minerals) indicating weaker bottom currents. A second rock-magnetic parameter, magnetic mineral accumulation rate (KaT) was used to indicate variations in terrigenous flux. Planktonic foraminiferal and carbonate accumulation rates (Pfar and CaC03ar) are used as indicators of pelagic flux. Absolute age assignments are based on correlation between the planktonic foraminiferal oxygen-isotope variations for Site 646 and the SPECMAP master oxygen-isotope curve. Cross-correlation analyses of the parameters that we studied with respect to the SPECMAP curve suggest that from oxygen-isotope stages 21 to 11, sedimentation rate, KaT, X, CaCO3ar, and Pfar were at their maximums, whereas XARM/X was at its minimum during peak interglacials (i.e., 0 k.y. lag time with respect to minimum ice volume). However, all parameters we examined lag behind minimum ice volume from stages 11 to 1, indicating a change in timing of both pelagic and terrigenous fluxes at approximately 400 k.y. BP. The negative correlation coefficient between XARM/X and the SPECMAP curve further suggest that finer-grained magnetic minerals are deposited during glacial periods, which probably reflects weaker bottom currents. The shift observed in the lag times of parameters examined with respect to the SPECMAP record is attributed to a change in significance of orbital parameters. Spectral results exhibit strong power in eccentricity (about 100 k.y.) throughout the record. Kap X, CaCO3flr, and Pfar show significant power in obliquity (about 41 k.y.), whereas XARM/X shows significant power at 73 k.y. from stages 21 to 11. The 73-k.y. period in XARM/X is near the difference tone of obliquity and eccentricity: 1/43-1/102 = 1/69. Kar and XARM/X show power only in eccentricity from stages 11 to 1. X and Pfar show significant power in precession (about 18 and 22 k.y.) whereas CaC03ar has power at 34 k.y, which could be a combination of precession and obliquity. The shift in power of orbital parameters may by attributed to the effect of the about 413-k.y. signal of eccentricity.

Benthic oxygen isotope and dust flux record from ODP Site 108-659 in the North Atlantic

High-resolution benthic oxygen isotope and dust flux records from Ocean Drilling Program site 659 have been analyzed to extend the astronomically calibrated isotope timescale for the Atlantic from 2.85 Ma back to 5 Ma. Spectral analysis of the delta18O record indicates that the 41-kyr period of Earth's orbital obliquity dominates the Pliocene record. This is shown to be true regardless of fundamental changes in the Earth's climate during the Pliocene. However, the cycles of Sahelian aridity fluctuations indicate a shift in spectral character near 3 Ma. From the early Pliocene to 3 Ma, the periodicities were dominantly precessional (19 and 23 kyr) and remained strong until 1.5 Ma. Subsequent to 3 Ma, the variance at the obliquity period (41 kyr) increased. The timescale tuned to precession suggests that the Pliocene was longer than previously estimated by more than 0.5 m.y. The tuned ages for the magnetic boundaries Gauss/Gilbert and Top Cochiti are about 6-8% older than the ages of the conventional timescale. A major phase of Pliocene northern hemisphere ice growth occurred between 3.15 Ma and 2.5 Ma. This was marked by a gradual increase in glacial Atlantic delta18O values of 1per mil and an increase in amplitude variations by up to 1.5 per mil, much larger than in the Pacific deepwater record (site 846). The first maxima occured in cold stages G6-96 between 2.7 Ma and 2.45 Ma. Prior to 3 Ma, the isotope record is characterized by predominantly low amplitude fluctuations (< 0.7 per mil). When obliquity forcing was at its minimum between 4.15 and 3.6 Ma and during the Kaena interval, delta18O amplitude fluctuations were minimal. From 4.9 to 4.3 Ma, the delta18O values decreased by about 0.5 per mil, reaching a long-term minimum at 4.15 Ma, suggesting higher deepwater temperatures or a deglaciation. Deepwater cooling and/or an increase in ice volume is indicated by a series of short-term delta18O fluctuations between 3.8 and 3.6 Ma.

Total organic carbon content and redox-sensitiv trace metal concentrations of selected late Cenoman claystones from DSDP Hole 93-603B

The Cenomanian/Turonian (C/T) intervals at DSDP Sites 105 and 603B from the northern part of the proto-North Atlantic show high amplitude, short-term cyclic variations in total organic carbon (TOC) content. The more pronounced changes in TOC are also reflected by changes in lithology from green claystones (TOC<1%) to black claystones (TOC>1%). Although their depositional history was different, the individual TOC cycles at Sites 105 and 603B can be correlated using stable carbon isotope stratigraphy. Sedimentation rates obtained from the isotope stratigraphy and spectral analyses indicate that these cycles were predominately precession controlled. The coinciding variations in HI, OI, delta13Corg and the abundance of marine relative to terrestrial biomarkers, as well as the low abundance of lignin pyrolysis products generated from the kerogen of the black claystones, indicate that these cyclic variations reflect changes in the contribution of marine organic matter (OM). The cooccurrence of lamination, enrichment of redox-sensitive trace metals and presence of molecular fossils of pigments from green sulfur bacteria indicate that the northern proto-North Atlantic Ocean water column was periodically euxinic from the bottom to at least the base of the photic zone (<150 m) during the deposition of the black claystones. In contrast, the green claystones are bioturbated, are enriched in Mn, do not show enrichments in redox-sensitive trace metals and show biomarker distributions indicative of long oxygen exposure times, indicating more oxic water conditions. At the same time, there is evidence (e.g., abundance of biogenic silica and significant 13C-enrichment for OC of phytoplanktic origin) for enhanced primary productivity during the deposition of the black claystones. We propose that increased primary productivity periodically overwhelmed the oxic OM remineralisation potential of the bottom waters resulting in the deposition of OM-rich black claystones. Because the amount of oxygen used for OM remineralisation exceeded the amount supplied by diffusion and deep-water circulation, the northern proto-North Atlantic became euxinic during these periods. Both Sites 105 and 603B show trends of continually increasing TOC contents and HI values of the black claystones up section, which most likely resulted from both enhanced preservation due to increased anoxia and increased production of marine OM during oceanic anoxic event 2 (OAE2).

Magnetic susceptibility (MS) measurements on the Sancerre (France) drill-core

The Toarcian Oceanic Anoxic Event (T-OAE) of the early Jurassic period involves one of the largest perturbations of the carbon cycle in the past 250 Ma, recorded by a pronounced negative carbon-isotope excursion (CIE). Numerous studies have focused on potential causes of the T-OAE and CIE, but are hampered by an uncertain timescale. Here we present high-resolution (~2 kyr) magnetic susceptibility (MS) measurements from the marine marls of the Sancerre-Couy drill-core, southern Paris Basin, spanning the entire Toarcian Stage. The MS variations document a rich series of sub-Milankovitch to Milankovitch frequencies (precession, obliquity and eccentricity) with the periodic g2-g5 (405 kyr) and quasi-periodic g4-g3 (~2.4 Myr Cenozoic mean periodicity) eccentricity terms being the most prominent. The MS-related g4-g3 variation reflects third-order eustatic sequences, and constrains the sequence stratigraphic framework of the Toarcian Stage. In addition, MS variations reveal a modulation of g2-g5 by g4-g3 eccentricity related cycles, suggesting that sea-level change was the main control on the deposition of the Toarcian Sancerre marls, in tune with the astro-climatic frequencies. The stable 405 kyr cyclicity constrains a minimum duration of the Toarcian Stage to ~8.3 Myr, and the well documented CIE, associated with the T-OAE, to ~300 to 500 kyr. The 405 kyr MS timescale calibrates the periodicity of the prominent high-frequency d13C cycles that occur in the decreasing part of the CIE to 30 to 34 kyr, consistent with the Toarcian obliquity period predicted for an Earth experiencing sustained tidal dissipation.

Geochemistry of sediments from the upwelling off Morocco

The high-productive upwelling area off Morocco is part of one of the four major trade-wind driven continental margin upwelling zones in the world oceans. While coastal upwelling occurs mostly on the shelf, biogenic particles derived from upwelling are deposited mostly at the upper continental slope. Nutrient-rich coastal water is transported within the Cape Ghir filament region at 30°N up to several hundreds of kilometers offshore. Both upwelling intensity and filament activity are dependent on the strength of the summer Trades. This study is aimed to reconstruct changes in trade wind intensity over the last 250,000 years by the analysis of the productivity signal contained in the sedimentary biogenic particles of the continental slope and beneath the Cape Ghir filament. Detailed geochemical and geophysical analyses (TOC, carbonate, C/N, delta13Corg, delta15N, delta13C of benthic foraminifera, delta18O of benthic and planktic foraminifera, magnetic susceptibility) have been carried out at two sites on the upper continental slope and one site located further offshore influenced by the Cape Ghir filament. A second offshore site south of the filament was analyzed (TOC, magnetic susceptibility) to distinguish the productivity signal related to the filament signal from the general offshore variability. Higher productivity during glacial times was observed at all four sites. However, the variability of productivity during glacial times was remarkably different at the filament-influenced site compared to the upwelling-influenced continental slope sites. In addition to climate-related changes in upwelling intensity, zonal shifts of the upwelling area due to sea-level changes have impacted the sedimentary productivity record, especially at the continental slope sites. By comparison with other proxies related to the strength and direction of the prevailing winds (Si/Al ratio as grain-size indicator, pollen) the productivity record at the filament-influenced site reflects mainly changes in trade-wind intensity. Our reconstruction reveals that especially during glacial times trade-wind intensity was increased and showed a strong variability with frequencies related to precession.

Susceptibility measurements at Gubbio sections

The complex interplay between extraterrestrial events and earth-bound processes that triggered one of the greatest biological crises of the Phanerozoic requires a high resolution timescale. Detailed magnetic susceptibility measurements at the Contessa Highway and Bottaccione sections (Italy) span the Cretaceous-Paleogene boundary and reveal clear orbital signatures in the sedimentary record. Identification of precession and 405 kyr eccentricity cycles allows an estimate of 324+/-40 kyr for the duration of the Maastrichtian part of Chron C29r. We present in the same high resolution time frame sites in Spain and the North and South Atlantic and bio-horizons, biotic changes, stable isotopic excursions and the decrease in Osmium isotopes recorded in these sections. The onset of 187 Os/ 188 Os decrease coincides with the d13 C negative excursion K-PgE1, thus suggesting a first pulse in Deccan volcanism at 66.64 Ma. The K-PgE3 d13 C negative excursion is possibly the expression of a second pulse at 66.26 Ma. Late Maastrichtian d13 C negative excursions are of low intensity and span durations of one to two eccentricity cycles, whereas early Danian excursions are brief (about 30 kyr) and acute. In Biotic response to late Maastrichtian perturbations occurred with a delay of ca. 200 kyr after the beginning of K-PgE1 shortly before K-PgE3. The biotic perturbation could be thus either a delayed response to K-PgE1, or a direct response to K-PgE3, and possibly, a threshold response to the stepwise buildup of CO2 atmospheric injections. No delay is evident in response to early Danian hyperthermal events. These differences suggest that short-lived, volcanically-derived environmental perturbations were buffered within the stable late Maastrichtian oceanic realm whereas they were amplified by the more sensitive and highly disturbed early Danian oceanic ecosystem.