Orbitally Tuned Research Articles

Combining orbital tuning and direct dating approaches to age-depth model development for Chew Bahir, Ethiopia

The directly dated RRMarch2021 age model (Roberts et al., 2021) for the ∼293 m long composite core from Chew Bahir, southern Ethiopia, has provided a valuable chronology for long-term climate changes in northeastern Africa. However, the age model has limitations on shorter time scales (less than 1–2 precession cycles), especially in the time range <20 kyr BP (kiloyears before present or thousand years before 1950) and between ∼155 and 428 kyr BP. To address those constraints we developed a partially orbitally tuned age model. A comparison with the ODP Site 967 record of the wetness index from the eastern Mediterranean, 3300 km away but connected to the Ethiopian plateau via the River Nile, suggests that the partially orbitally tuned age model offers some advantages compared to the exclusively directly dated age model, with the limitation of the reduced significance of (cross) spectral analysis results of tuned age models in cause-effect studies. The availability of this more detailed age model is a prerequisite for further detailed spatiotemporal correlations of climate variability and its potential impact on the exchange of different populations of Homo sapiens in the region.

Orbital tuning for the middle Eocene to early Oligocene Monte Cagnero Section (Central Italy): Paleoenvironmental and paleoclimatic implications

During the middle Eocene to early Oligocene Earth transitioned from a greenhouse to an icehouse climate state. The interval comprises the Middle Eocene Climatic Optimum (MECO; ~40 Ma) and a subsequent long-term cooling trend that culminated in the Eocene-Oligocene transition (EOT; ~34 Ma) with the Oi-1 glaciation. Here, we present a refined calcareous nannofossil biostratigraphy and an orbitally tuned age model for the Monte Cagnero (MCA) section spanning the middle Eocene to the early Oligocene (~41 to ~33 Ma). Spectral analysis of magnetic susceptibility (MS) data displays strong cyclicities in the orbital frequency band allowing us to tune the identified 405 kyr eccentricity minima in the MS record to their equivalents in the astronomical solution. Our orbitally tuned age model allows us to estimate the position and duration of polarity chrons (C18 to C13) and compare them with other standard and orbitally tuned ages. We were also able to constrain the timing and duration of the MECO event, which coincides with a minimum in the 2.4 Myr and 405 kyr eccentricity cycles. Our study corroborates the previous estimated age for the base of the Rupelian stage (33.9 Ma) and estimates the base of the Priabonian stage in the MCA section to be at 37.4 Ma. Finally, calcareous nannofossils with known paleoenvironmental preferences suggest a gradual shift from oligotrophic to meso-eutrophic conditions with an abrupt change at ~36.8 Ma. Besides, nannofossil assemblages suggest that enhanced nutrient availability preceded water cooling at the late Eocene. Altogether, this evidence points to a poorly developed water column stratification prior to the cooling trend.

Age of Matuyama-Brunhes boundary constrained by U-Pb zircon dating of a widespread tephra

ABSTRACTThe youngest geomagnetic polarity reversal, the Matuyama-Brunhes boundary (MBB), provides an important datum plane for sediments, ice cores, and lavas. Its frequently cited age of 780 ka is based on orbital tuning of marine sedimentary records, and is supported by 40 Ar/ 39 Ar dating of Hawaiian lavas using recent age calibrations. Challenging this age, how-ever, are reports of younger astrochronological ages based on oxygen isotope stratigraphy of high-sedimentation-rate marine records, and cosmogenic nuclides in marine sediments and an Antarctic ice core. Here, we present a U-Pb zircon age of 772.7 ± 7.2 ka from a marine-deposited tephra just below the MBB in a forearc basin in Japan. U-Pb dating has a distinct advantage over 40 Ar/ 39 Ar dating in that it is relatively free from assumptions regarding stan-dardization and decay constants. This U-Pb zircon age, coupled with a newly obtained oxygen isotope chronology, yields an MBB age of 770.2 ± 7.3 ka. Our MBB age is consistent with those based on the latest orbitally tuned marine sediment records and on an Antarctic ice core. We provide the first direct comparison between orbital tuning, U-Pb dating, and magnetostratig-raphy for the MBB, fulfilling a key requirement in calibrating the geological time scale.INTRODUCTION

Correlating carbon and oxygen isotope events in early to middle Miocene shallow marine carbonates in the Mediterranean region using orbitally tuned chemostratigraphy and lithostratigraphy.

During the Miocene prominent oxygen isotope events (Mi‐events) reflect major changes in glaciation, while carbonate isotope maxima (CM‐events) reflect changes in organic carbon burial, particularly during the Monterey carbon isotope excursion. However, despite their importance to the global climate history they have never been recorded in shallow marine carbonate successions. The Decontra section on the Maiella Platform (central Apennines, Italy), however, allows to resolve them for the first time in such a setting during the early to middle Miocene. The present study improves the stratigraphic resolution of parts of the Decontra section via orbital tuning of high‐resolution gamma ray (GR) and magnetic susceptibility data to the 405 kyr eccentricity metronome. The tuning allows, within the established biostratigraphic, sequence stratigraphic, and isotope stratigraphic frameworks, a precise correlation of the Decontra section with pelagic records of the Mediterranean region, as well as the global paleoclimatic record and the global sea level curve. Spectral series analyses of GR data further indicate that the 405 kyr orbital cycle is particularly well preserved during the Monterey Event. Since GR is a direct proxy for authigenic uranium precipitation during increased burial of organic carbon in the Decontra section, it follows the same long‐term orbital pacing as observed in the carbon isotope records. The 405 kyr GR beat is thus correlated with the carbon isotope maxima observed during the Monterey Event. Finally, the Mi‐events can now be recognized in the δ18O record and coincide with plankton‐rich, siliceous, or phosphatic horizons in the lithology of the section.

Orbitally tuned timescale and astronomical forcing in the middle Eocene to early Oligocene

Abstract. Deciphering the driving mechanisms of Earth system processes, including the climate dynamics expressed as paleoceanographic events, requires a complete, continuous, and high-resolution stratigraphy that is very accurately dated. In this study, a robust astronomically calibrated age model was constructed for the middle Eocene to early Oligocene interval (31–43 Ma) in order to permit more detailed study of the exceptional climatic events that occurred during this time, including the middle Eocene climate optimum and the Eocene–Oligocene transition. A goal of this effort is to accurately date the middle Eocene to early Oligocene composite section cored during the Pacific Equatorial Age Transect (PEAT, IODP Exp. 320/321). The stratigraphic framework for the new timescale is based on the identification of the stable long eccentricity cycle in published and new high-resolution records encompassing bulk and benthic stable isotope, calibrated XRF core scanning, and magnetostratigraphic data from ODP Sites 171B-1052, 189-1172, 199-1218, and 207-1260 as well as IODP Sites 320-U1333, and 320-U1334 spanning magnetic polarity Chrons C12n to C20n. Subsequently orbital tuning of the records to the La2011 orbital solution was conducted. The resulting new timescale revises and refines the existing orbitally tuned age model and the geomagnetic polarity timescale from 31 to 43 Ma. The newly defined absolute age for the Eocene–Oligocene boundary validates the astronomical tuned age of 33.89 Ma identified at the Massignano, Italy, global stratotype section and point. The compilation of geochemical records of climate-controlled variability in sedimentation through the middle-to-late Eocene and early Oligocene demonstrates strong power in the eccentricity band that is readily tuned to the latest astronomical solution. Obliquity driven cyclicity is only apparent during 2.4 myr eccentricity cycle minima around 35.5, 38.3, and 40.1 Ma.

Chronology of Lake El'gygytgyn sediments – a combined magnetostratigraphic, palaeoclimatic and orbital tuning study based on multi-parameter analyses

Abstract. A 318-metre-long sedimentary profile drilled by the International Continental Scientific Drilling Program (ICDP) at Site 5011-1 in Lake El'gygytgyn, Far East Russian Arctic, has been analysed for its sedimentologic response to global climate modes by chronostratigraphic methods. The 12 km wide lake is sited off-centre in an 18 km large crater that was created by the impact of a meteorite 3.58 Ma ago. Since then sediments have been continuously deposited. For establishing their chronology, major reversals of the earth's magnetic field provided initial tie points for the age model, confirming that the impact occurred in the earliest geomagnetic Gauss chron. Various stratigraphic parameters, reflecting redox conditions at the lake floor and climatic conditions in the catchment were tuned synchronously to Northern Hemisphere insolation variations and the marine oxygen isotope stack, respectively. Thus, a robust age model comprising more than 600 tie points could be defined. It could be shown that deposition of sediments in Lake El'gygytgyn occurred in concert with global climatic cycles. The upper ~160 m of sediments represent the past 3.3 Ma, equivalent to sedimentation rates of 4 to 5 cm ka−1, whereas the lower 160 m represent just the first 0.3 Ma after the impact, equivalent to sedimentation rates in the order of 45 cm ka−1. This study also provides orbitally tuned ages for a total of 8 tephras deposited in Lake El'gygytgyn.

Direct correlation of Chinese speleothem δ18O and South China Sea planktonic δ18O: Transferring a speleothem chronology to the benthic marine chronology

Orbital tuning of benthic δ18O is a common approach for assigning ages to ocean sediment records. Similar environmental forcing of the northern South China Sea and the southeast Asian cave regions allows for transfer of the speleothem δ18O radiometric chronology to the planktonic and benthic δ18O records from Ocean Drilling Program Site 1146, yielding a new chronology with 41 radiometrically calibrated datums, spanning the past 350 kyr. This approach also provides for an independent assessment of the accuracy of the orbitally tuned benthic δ18O chronology for the last 350 kyr. The largest differences relative to the latest chronology occur in marine isotope stages (MIS) 5.4, 5.5, 6, 7, and 9.3. Prominent suborbital‐scale structure believed to be global in nature is identified within MIS 5.4 and MIS 7.2. On the basis of the radiometrically calibrated chronology, the time constant of the ice sheet is found to be 5.4 kyr at the precession band (light δ18O lags precession minima by −55.4°) and 10.4 kyr at the obliquity band (light δ18O lags obliquity maxima by 57.4°). These values are significantly shorter than the single 17 kyr time constant originally estimated by Imbrie et al. (1984), based primarily on the timing of terminations I and II and the 15 kyr time constant used by Lisiecki and Raymo (2005).

To tune or not to tune: Detecting orbital variability in Oligo-Miocene climate records

We address the problem of detecting quasi-periodic variability at orbital frequencies within pre-Pleistocene climate records using depth-derived and orbitally tuned chronologies. Many studies describing orbital variability in pre-Pleistocene sediment hosted isotope records employ climatic records that are set on orbitally tuned chronologies, without accounting for the bias in spectral power estimates introduced by orbital tuning. In this study we develop a method to quantify the effects of tuning upon spectral estimates and, in particular, to more properly determine the statistical significance of spectral peaks associated with orbital frequencies. We apply this method to two marine sediment δ18O records spanning the Oligo-Miocene, from ODP cores 1090 and 1218. We find that using linear age–depth relationships reveals statistically significant spectral peaks matching eccentricity in core 1090, and obliquity and precession in core 1218, where the last appears most significant. Tuning the chronologies to the orbital solutions of Laskar et al. (2004) increases the statistical significance of the precession peak, whereas the obliquity and eccentricity peaks become indistinguishable from those expected from tuning noise. This result can be understood in that tuning records with high signal to noise ratios tends to lead to more significant spectral peaks, whereas a linear age–depth relationship is better suited for detecting peaks when signal to noise ratios are low. We also demonstrate this concept using synthetic records.

Towards orbital dating of the EPICA Dome C ice core using δO&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt;/N&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt;

Abstract. Based on a composite of several measurement series performed on ice samples stored at −25 °C or −50 °C, we present and discuss the first δO2/N2 record of trapped air from the EPICA Dome C (EDC) ice core covering the period between 300 and 800 ka (thousands of years before present). The samples stored at −25 °C show clear gas loss affecting the precision and mean level of the δO2/N2 record. Two different gas loss corrections are proposed to account for this effect, without altering the spectral properties of the original datasets. Although processes at play remain to be fully understood, previous studies have proposed a link between surface insolation, ice grain properties at close-off, and δO2/N2 in air bubbles, from which orbitally tuned chronologies of the Vostok and Dome Fuji ice core records have been derived over the last four climatic cycles. Here, we show that limitations caused by data quality and resolution, data filtering, and uncertainties in the orbital tuning target limit the precision of this tuning method for EDC. Moreover, our extended record includes two periods of low eccentricity. During these intervals (around 400 ka and 750 ka), the matching between δO2/N2 and the different insolation curves is ambiguous because some local insolation maxima cannot be identified in the δO2/N2 record (and vice versa). Recognizing these limitations, we restrict the use of our δO2/N2 record to show that the EDC3 age scale is generally correct within its published uncertainty (6 kyr) over the 300–800 ka period.

Orbital tuning, eccentricity, and the frequency modulation of climatic precession

The accuracy of geologic chronologies can, in principle, be improved through orbital tuning, the systematic adjustment of a chronology to bring the associated record into greater alignment with an orbitally derived signal. It would be useful to have a general test for the success of orbital tuning, and one proposal has been that eccentricity ought to covary with the amplitude envelope associated with precession variability recorded in tuned geologic records. A common procedure is to filter a tuned geologic record so as to pass precession period variability and compare the amplitude modulation of the resulting signal against eccentricity. There is a reasonable expectation for such a relationship to be found in paleoclimate records because the amplitude of precession forcing depends upon eccentricity. However, there also exists a relationship between eccentricity and the frequency of precession such that orbital tuning generates eccentricity-like amplitude modulation in filtered signals, regardless of the accuracy of the chronology or the actual presence of precession. This relationship results from the celestial mechanics governing eccentricity and precession and from the interaction between frequency modulation and amplitude modulation caused by filtering. When the eccentricity of Earth's orbit is small, the frequency of climatic precession undergoes large variations and less precession energy is passed through a narrow-band filter. Furthermore, eccentricity-like amplitude modulation is routinely obtained from pure noise records that are orbitally tuned to precession and then filtered. We conclude that the presence of eccentricity-like amplitude modulation in precession-filtered records does not support the accuracy of orbitally tuned time scales.

An ice core perspective on the age of the Matuyama–Brunhes boundary

Two intervals of enhanced 10Be flux thought to be associated with periods of low dipole intensity and identified as the Matuyama–Brunhes transition and a precursor event have been observed in the bottom section of the EPICA Dome C ice core. The peaks span 764–776 ka and 788–798 ka on the new EDC3 chronology with a stated absolute age uncertainty of 6 ka (2 σ). This chronology uses orbital tuning of atmospheric oxygen-18 ( δ 18O atm) to correct for anomalies in ice flow in the bottom 500 m of the core. An additional 28 δ 18O atm data points have been measured to improve resolution and verify the accuracy of the tuning and the stated timescale uncertainty. Both the dating of the increased 10Be, and that relative to climatic records, are compared to paleointensity records found in orbitally tuned marine sediments. The mid-point of the 10Be peak associated with the M-B is approximately 10 ka younger than the age determined radioisotopically from lavas with transitional orientations, taking into account recent revisions to the 40Ar/ 39Ar dating standard and improved precision. Climatic constraints on the EDC3 agescale make an error of this magnitude in the ice chronology implausible. This age difference, however, is consistent with recent modeling suggesting that directional changes are spatially asynchronous, and may precede the dipole intensity minimum in some locations. Although formally less precise than the published age from astrochronologically dated marine sediments, ice core ages are potentially more accurate because they are not subject to lock-in depth uncertainties.

A radiometric calibration of the SPECMAP timescale

The astronomical theory of climate change asserts that Earth's climate is affected by changes in its orbit, which vary the seasonal and latitudinal distribution of solar radiation. This theory is the basis of the orbitally tuned SPECMAP timescale. A key constraint for this important chronology was the mid-point of the Penultimate Deglaciation, initially dated to 127,000 years ago. Recent work suggests this event may be considerably older, casting doubt on the astronomical theory, the SPECMAP timescale, and the accuracy of orbitally tuned chronologies. Difficulties with U/Th coral dating of sea-level events have impeded progress on this problem, because most corals are not closed systems. Here, we use a new approach to U/Th dating that corrects for open-system behavior and produces a sea-level curve of sufficient resolution to confidently correlate with SPECMAP over the last 240,000 years, permitting a reassessment of both this critical chronology and a central tenet of climate change theory. High-precision ages for 24 oxygen isotope events provide a 240,000-year chronology for marine δ 18O records that is independent of orbital tuning assumptions. Although there appear to be significant differences between the radiometric and orbitally tuned timescales near the lastglacial maximum and at the Marine Isotope Stage 7/6 boundary, a comparison of radiometric and SPECMAP ages for identical isotope events suggest that the SPECMAP timescale is quite accurate and that its errors were, in general, overestimated. Despite suborbital complexity, orbital cyclicity is clearly evident in our record. High-amplitude sea-level oscillations at periods greater than ∼20,000 years are very close in phase to summer insolation in the Northern Hemisphere. Although sea-level changes cannot be uniquely tied to a specific season or latitude of insolation forcing, the simplest explanation is that long-period, high-amplitude sea-level change is linked to Northern Hemisphere insolation forcing. These results validate the principles of orbital tuning and suggest such timescales are generally robust.

A stratigraphic network across the Subtropical Front in the central South Atlantic: Multi-parameter correlation of magnetic susceptibility, density, X-ray fluorescence and δ18O records

Abstract Computer aided multi-parameter signal correlation is used to develop a common high-precision age model for eight gravity cores from the subtropical and subantarctic South Atlantic. Since correlations between all pairs of multi-parameter sequences are used, and correlation errors between core pairs (A, B) and (B, C) are controlled by comparison with (A, C), the resulting age model is called a stratigraphic network. Precise inter-core correlation is achieved using high-resolution records of magnetic susceptibility κ, wet bulk density ρ and X-ray fluorescence scans of elemental composition. Additional δ18O records are available for two cores. The data indicate nearly undisturbed sediment series and the absence of significant hiatuses or turbidites. After establishing a high-precision common depth scale by synchronously correlating four densely measured parameters (Fe, Ca, κ, ρ), the final age model is obtained by simultaneously fitting the aligned δ18O and κ records of the stratigraphic network to orbitally tuned oxygen isotope [J. Imbrie, J. D. Hays, D. G. Martinson, A. McIntyre, A. C. Mix, J. J. Morley, N. G. Pisias, W. L. Prell, N. J. Shackleton, The orbital theory of Pleistocene climate: support from a revised chronology of the marine δ18O record, in: A. Berger, J. Imbrie, J. Hays, G. Kukla, B. Saltzman (Eds.), Milankovitch and Climate: Understanding the Response to Orbital Forcing, Reidel Publishing, Dordrecht, 1984, pp. 269-305; D. Martinson, N. Pisias, J. Hays, J. Imbrie, T. C. Moore Jr., N. Shackleton, Age dating and the orbital theory of the Ice Ages: development of a high-resolution 0 to 300.000-Year chronostratigraphy, Quat. Res. 27 (1987) 1-29.] or susceptibility stacks [T. von Dobeneck, F.Schmieder, Using rock magnetic proxy records for orbital tuning and extended time series analyses into the super-and sub-Milankovitch Bands, in: G. Fischer, G. Wefer (Eds.), Use of proxies in paleoceanography: Examples from the South Atlantic, Springer-Verlag, Berlin (1999), pp. 601-633.]. Besides the detection and elimination of errors in single records, the stratigraphic network approach allows to check the intrinsic consistency of the final result by comparing it to the outcome of more restricted alignment procedures. The final South Atlantic stratigraphic network covers the last 400 kyr south and the last 1200 kyr north of the Subtropical Front (STF) and provides a highly precise age model across the STF representing extremely different sedimentary regimes. This allows to detect temporal shifts of the STF by mapping δMn / Fe. It turns out that the apparent STF movements by about 200 km are not directly related to marine oxygen isotope stages.

Climato-hydrological environment inferred from Lake Baikal sediments based on an automatic orbitally tuned age model

We propose a new automatic orbital tuning algorithm to adjust climatic signals to insolation. This can tune two signals with periodicity only without shape similarity. A Genetic Algorithm (GA) is used as an optimizing method. The new age model for the Brunhes epoch in Lake Baikal core BDP98 defined climatic shifts at about 250 kyr B.P., 350 kyr B.P. and 700 kyr B.P. The sedimentation rate for the interval from 350 kyr B.P. to 700 kyr B.P. was comparatively low and stable. This new model also indicates that there was a stadial during the super interglacial period (MIS 11), and that the terrestrial response to climate change was complex in this period.

A comparison of the Vostok ice deuterium record and series from Southern Ocean core MD 88-770 over the last two glacial-interglacial cycles

Taking advantage of the fact that the Vostok deuterium (δD) record now covers almost two entire climatic cycles, we have applied the orbital tuning approach to derive an age-depth relation for the Vostok ice core, which is consistent with the SPECMAP marine time scale. A second age-depth relation for Vostok was obtained by correlating the ice isotope content with estimates of sea surface temperature from Southern Ocean core MD 88-770. Both methods lead to a close correspondence between Vostok and MD 88-770 time series. However, the coherence between the correlated δD and insolation is much lower than between the orbitally tuned 8D and insolation. This reflects the lower accuracy of the correlation method with respect to direct orbital tuning. We compared the ice and marine records, set in a common temporal framework, in the time and frequency domains. Our results indicate that changes in the Antarctic air temperature quite clearly lead variations in global ice volume in the obliquity and precession frequency bands. Moreover, the average phase we estimated between the filtered δD and insolation signals at precessional frequencies indicates that variations in the southern high latitude surface temperature could be induced by changes in insolation taking place during a large period of the summer in northern low latitudes or winter in southern low latitudes. The relatively large lag found between Vostok δD variations and obliquity-driven changes in insolation suggests that variations in the local radiative balance are not the only mechanism responsible for the variability in surface temperature at those frequencies. Finally, in contrast to the cross-spectral analysis method used in previous studies, the method we use here to estimate the phases can reveal errors in cross-correlations with orbitally tuned chronologies.