Astronomical calibration of the geological timescale: closing the middle Eocene gap

To explore cause and consequences of past climate change, very accurate age models such as those provided by the astronomical timescale (ATS) are needed. Beyond 40 million years the accuracy of the ATS critically depends on the correctness of orbital models and radioisotopic dating techniques. Discrepancies in the age dating of sedimentary successions and the lack of suitable records spanning the middle Eocene have prevented development of a continuous astronomically calibrated geological timescale for the entire Cenozoic Era. We now solve this problem by constructing an independent astrochronological stratigraphy based on Earth's stable 405 kyr eccentricity cycle between 41 and 48 million years ago (Ma) with new data from deep-sea sedimentary sequences in the South Atlantic Ocean. This new link completes the Paleogene astronomical timescale and confirms the intercalibration of radioisotopic and astronomical dating methods back through the Paleocene–Eocene Thermal Maximum (PETM, 55.930 Ma) and the Cretaceous–Paleogene boundary (66.022 Ma). Coupling of the Paleogene 405 kyr cyclostratigraphic frameworks across the middle Eocene further paves the way for extending the ATS into the Mesozoic

Orbital forcing of the Paleocene and Eocene carbon cycle

This is the final version of the article. Available from the publisher via the DOI in this record.Multimillion-year proxy records across the Paleocene and Eocene show prominent variations on orbital time scales. The cycles, which have been identified at various sites across the globe, preferentially concentrate spectral power at eccentricity and precessional frequencies. It is evident that these cycles are an expression of changes in global climate and carbon cycling paced by astronomical forcing. However, little is currently known about the link between orbital forcing and the carbon cycle-climate system and the amplitude of associated atmospheric CO2 variations. Here we use simple and complex carbon cycle models to explore the basic effect of different orbital forcing schemes and noise on the carbon cycle. Our primary modeling target is the high-resolution, ∼7.7 Myr long, benthic isotope record at Ocean Drilling Program Site 1262 in the South Atlantic. For direct insolation forcing (as opposed to artificial eccentricity-tilt-precession), one major challenge is understanding how the system transfers spectral power from high to low frequencies. We discuss feasible solutions, including insolation transformations analogous to electronic AC-DC conversion (DC’ing). Regarding mechanisms, we focus on tropical insolation and a long-term carbon imbalance in terrestrial organic burial/oxidation but do not rule out other scenarios. Our analysis shows that high-latitude mechanisms are unlikely drivers of orbitally paced changes in the late Paleocene-early Eocene (LPEE) Earth system. Furthermore, we provide constraints on the origin and isotopic composition of a possible LPEE cyclic carbon imbalance/source responding to astronomical forcing. Our simulations also reveal a mechanism for the large 13C-eccentricity lag at the 400 kyr period observed in Paleocene, Oligocene, and Miocene sections. We present the first estimates of orbital-scale variations in atmospheric CO2 during the late Paleocene and early EoceneThis research was supported by U.S. NSF grants OCE12-20615 and OCE16-58023 to R.E.Z. and J.C.Z. and the Deutsche Forschungsgemeinschaft (DFG) to T.W

A high‐resolution benthic stable‐isotope record for the South Atlantic: implications for orbital scale changes in Late Paleocene–Early Eocene climate and circulation

This is the author accepted manuscript. The final version is available from Elsevier via the DOI in this record.The Late Paleocene and Early Eocene were characterized by warm greenhouse climates, punctuated by a series of rapid warming and ocean acidification events known as “hyperthermals”, thought to have been paced or triggered by orbital cycles. While these hyperthermals, such as the Paleocene Eocene Thermal Maximum (PETM), have been studied in great detail, the background low-amplitude cycles seen in carbon and oxygen-isotope records throughout the Paleocene–Eocene have hitherto not been resolved. Here we present a 7.7 million year (myr) long, high-resolution, orbitally-tuned, benthic foraminiferal stable-isotope record spanning the late Paleocene and early Eocene interval (∼52.5–60.5 Ma) from Ocean Drilling Program (ODP) Site 1262, South Atlantic. This high resolution (∼2–4 kyr) record allows the changing character and phasing of orbitally-modulated cycles to be studied in unprecedented detail as it reflects the long-term trend in carbon cycle and climate over this interval. The main pacemaker in the benthic oxygen-isotope (δ18O) and carbon-isotope (δ13C) records from ODP Site 1262, are the long (405 kyr) and short (100 kyr) eccentricity cycles, and precession (21 kyr). Obliquity (41 kyr) is almost absent throughout the section except for a few brief intervals where it has a relatively weak influence. During the course of the Early Paleogene record, and particularly in the latest Paleocene, eccentricity-paced negative carbon-isotope excursions (δ13C, CIEs) and coeval negative oxygen-isotope (δ18O) excursions correspond to low carbonate (CaCO3) and coarse fraction (%CF) values due to increased carbonate dissolution, suggesting shoaling of the lysocline and accompanied changes in the global exogenic carbon cycle. These negative CIEs and δ18O events coincide with maxima in eccentricity, with changes in δ18O leading changes in δ13C by ∼6 (±5) kyr in the 405-kyr band and by ∼3 (±1) kyr in the higher frequency 100-kyr band on average. However, these phase lags are not constant, with the lag in the 405-kyr band extending from ∼4 (±5) kyr to ∼21 (±2) kyr from the late Paleocene to the early Eocene, suggesting a progressively weaker coupling of climate and the carbon-cycle with time. The higher amplitude 405-kyr cycles in the latest Paleocene are associated with changes in bottom water temperature of 2–4 °C, while the most prominent 100 kyr-paced cycles can be accompanied by changes of up to 1.5 °C. Comparison of the 1262 record with a lower resolution, but orbitally-tuned benthic record for Site 1209 in the Pacific allows for verification of key features of the benthic isotope records which are global in scale including a key warming step at 57.7 Ma.Thanks to Alexis Kersey for picking foraminifera and assisting with sample processing (at UCSC), and to Walker Weir, Alejandro Aguilar, and Phillip Staudigel for lab assistance, and to Dyke Andreasen and Chih-Ting Hsieh for stable-isotope support (UCSC). Thanks to Barbara Donner (MARUM) for coordinating foraminifera picking, and to Monika Segl and her team (MARUM) for stable-isotope analyses. We thank Roy Wilkens (Hawaii) for core images analysis. Sediment samples were supplied by the Ocean Drilling Program. Funding for this project was provided by NSF grant (grant number EAR-0628719) to J.Z. and DFG grants (RO 1113/2 through RO 1113/4) to U.R

Data report: depths of Site U1553 off-splice data adjusted to the Site U1553 splice, IODP Expedition 378

A near-complete sedimentary sequence was spliced together for the upper part of International Ocean Discovery Program (IODP) Holes U1553A, U1553B, and U1553E. Poor core recovery precluded a complete splice for the deeper section cored in Holes U1553C and U1553D. The history of Deep Sea Drilling Project Site 277, which was cored nearby, suggests that the Site U1553 splice will be heavily sampled and that eventually samples will be taken from intervals of core that are not included in the splice (i.e., off-splice). Although the depths of all cores have been shifted to a common scale during the splicing process by aligning significant features shared by cores from the different holes, core disturbance and natural variability often lead to misalignment between features in the splice and the same features in off-splice data. To remedy this problem for future sampling, data from off-splice intervals are squeezed or stretched to match spliced intervals using a set of tie points between the splice and off-splice data. The difference in depths can be significant when considering sedimentation rates and orbital periods of precession, obliquity, and eccentricity and sometimes even change the phase relationship compared to the splice. Results are presented as tables of tie points between each hole and the splice that can be used to interpolate the proper splice depth of off-splice samples

Data report: splice adjustment for Site U1553

Postcruise examination of the data splice for International Ocean Discovery Program Expedition 378 Site U1553, in light of new X-ray fluorescence data, revealed three cores from Hole U1553E that were misaligned. These cores have been shifted to fill in some gaps in the original splice

High resolution cyclostratigraphy of the early Eocene – new insights into the origin of the Cenozoic cooling trend

Here we present a high-resolution cyclostratigraphy based on X-ray fluorescence (XRF) core scanning data from a new record retrieved from the tropical western Atlantic (Demerara Rise, ODP Leg 207, Site 1258). The Eocene sediments from ODP Site 1258 cover magnetochrons C20 to C24 and show well developed cycles. This record includes the missing interval for reevaluating the early Eocene part of the Geomagnetic Polarity Time Scale (GPTS), also providing key aspects for reconstructing high-resolution climate variability during the Early Eocene Climatic Optimum (EECO). Detailed spectral analysis demonstrates that early Eocene sedimentary cycles are characterized by precession frequencies modulated by short (100 kyr) and long (405 kyr) eccentricity with a generally minor obliquity component. Counting of both the precession and eccentricity cycles results in revised estimates for the duration of magnetochrons C21r through C24n. Our cyclostratigraphic framework also corroborates that the geochronology of the Eocene Green River Formation (Wyoming, USA) is still questionable mainly due to the uncertain correlation of the "Sixth tuff" to the GPTS. <br><br> Right at the onset of the long-term Cenozoic cooling trend the dominant eccentricity-modulated precession cycles of ODP Site 1258 are interrupted by strong obliquity cycles for a period of ~800 kyr in the middle of magnetochron C22r. These distinct obliquity cycles at this low latitude site point to (1) a high-latitude driving mechanism on global climate variability from 50.1 to 49.4 Ma, and (2) seem to coincide with a significant drop in atmospheric CO<sub>2</sub> concentration below a critical threshold between 2- and 3-times the pre-industrial level (PAL). The here newly identified orbital configuration of low eccentricity in combination with high obliquity amplitudes during this ~800-kyr period and the crossing of a critical <i>p</i>CO<sub>2</sub> threshold may have led to the formation of the first ephemeral ice sheet on Antarctica as early as ~50 Ma ago

Calcareous nannofossils across the Eocene-Oligocene transition: Preservation signals and biostratigraphic remarks from ODP Site 1209 (NW Pacific, Shatsky Rise) and IODP Hole U1411B (NW Atlantic Ocean, Newfoundland Ridge)

This work provides a detailed biostratigraphic correlation through the Eocene-Oligocene Transition (EOT), based on an integrated stratigraphic approach and the study of calcareous nannofossils, between two disparate sites, one in the NW Atlantic (IODP Hole U1411B) and one in the NW Pacific (ODP Site 1209). The precise site-to-site correlation provided by these data allows for a comparison of carbonate preservation across the EOT including identification of the main post-depositional processes that impact the calcareous nannofossil ooze at Site 1209. The main aim of this work is to understand the extent to which the bulk δ18O and δ13C records and their sources (mainly calcareous nannofossils) are altered by diagenesis. Our detailed SEM study highlights some differences before, during and after the EOT, suggesting local diagenetic dynamics. At Site 1209, a distinctive change, both in nannofossil assemblage composition and preservation state, is observed from the pre-EOT phase to the Late Eocene Event (LEE), with a shift in the dominant process from dissolution to recrystallisation. Surprisingly, despite the overall poor preservation, only the interval between 141 and 142.4 (adj. rmcd) was compromised in term of isotopic values and assemblage diversity and abundance. This interval, recorded in the upper Eocene, was characterized by severe dissolution, concomitant with deposition of secondary calcite on solution-resistant forms. Diagenetic processes have strongly biased the δ18O isotopic signal, resulting in a positive oxygen isotope anomaly through the upper Eocene that is difficult to reconcile with other published trends. For the remaining time intervals, diagenesis seems not to have altered the bulk δ18O profile, which closely resembles that of other sites across the world, and is particularly consistent with other data from the Pacific Ocean. In summary, the impact of diagenesis on nannofossil preservation even if clearly visible both in SEM and optical microscope observations does not always cause a pervasive alteration of the primary isotopic signal and can instead provide important clues on local depositional dynamics

Late Miocene to Holocene high-resolution eastern equatorial Pacific carbonate records: stratigraphy linked by dissolution and paleoproductivity

Coherent variation in CaCO3 burial is a feature of the Cenozoic eastern equatorial Pacific. Nevertheless, there has been a long-standing ambiguity in whether changes in CaCO3 dissolution or changes in equatorial primary production might cause the variability. Since productivity and dissolution leave distinctive regional signals, a regional synthesis of data using updated age models and high-resolution stratigraphic correlation is an important constraint to distinguish between dissolution and production as factors that cause low CaCO3. Furthermore, the new chronostratigraphy is an important foundation for future paleoceanographic studies. The ability to distinguish between primary production and dissolution is also important to establish a regional carbonate compensation depth (CCD). We report late Miocene to Holocene time series of XRF-derived (X-ray fluorescence) bulk sediment composition and mass accumulation rates (MARs) from eastern equatorial Pacific Integrated Ocean Drilling Program (IODP) sites U1335, U1337, and U1338 and Ocean Drilling Program (ODP) site 849, and we also report bulk-density-derived CaCO3 MARs at ODP sites 848, 850, and 851. We use physical properties, XRF bulk chemical scans, and images along with available chronostratigraphy to intercorrelate records in depth space. We then apply a new equatorial Pacific age model to create correlated age records for the last 8 Myr with resolutions of 1–2 kyr. Large magnitude changes in CaCO3 and bio-SiO2 (biogenic opal) MARs occurred within that time period but clay deposition has remained relatively constant, indicating that changes in Fe deposition from dust is only a secondary feedback to equatorial productivity. Because clay deposition is relatively constant, ratios of CaCO3 % or biogenic SiO2 % to clay emulate changes in biogenic MAR. We define five major Pliocene–Pleistocene low CaCO3 % (PPLC) intervals since 5.3 Ma. Two were caused primarily by high bio-SiO2 burial that diluted CaCO3 (PPLC-2, 1685–2135 ka, and PPLC-5, 4465–4737 ka), while three were caused by enhanced dissolution of CaCO3 (PPLC-1, 51–402 ka, PPLC-3, 2248–2684 ka, and PPLC-4, 2915–4093 ka). Regional patterns of CaCO3 % minima can distinguish between low CaCO3 caused by high diatom bio-SiO2 dilution versus lows caused by high CaCO3 dissolution. CaCO3 dissolution can be confirmed through scanning XRF measurements of Ba. High diatom production causes lowest CaCO3 % within the equatorial high productivity zone, while higher dissolution causes lowest CaCO3 percent at higher latitudes where CaCO3 production is lower. The two diatom production intervals, PPLC-2 and PPLC-5, have different geographic footprints from each other because of regional changes in eastern Pacific nutrient storage after the closure of the Central American Seaway. Because of the regional variability in carbonate production and sedimentation, the carbonate compensation depth (CCD) approach is only useful to examine large changes in CaCO3 dissolution

Solar total and spectral irradiance reconstruction over the last 9000 years

Changes in solar irradiance and in its spectral distribution are among the main natural drivers of the climate on Earth. However, irradiance measurements are only available for less than four decades, while assessment of solar influence on Earth requires much longer records. The aim of this work is to provide the most up-to-date physics-based reconstruction of the solar total and spectral irradiance (TSI/SSI) over the last nine millennia. The concentrations of the cosmogenic isotopes 14C and 10Be in natural archives have been converted to decadally averaged sunspot numbers through a chain of physics-based models. TSI and SSI are reconstructed with an updated SATIRE model. Reconstructions are carried out for each isotope record separately, as well as for their composite. We present the first ever SSI reconstruction over the last 9000 years from the individual 14C and 10Be records as well as from their newest composite. The reconstruction employs physics-based models to describe the involved processes at each step of the procedure. Irradiance reconstructions based on two different cosmogenic isotope records, those of 14C and 10Be, agree well with each other in their long-term trends despite their different geochemical paths in the atmosphere of Earth. Over the last 9000 years, the reconstructed secular variability in TSI is of the order of 0.11%, or 1.5 W/m2. After the Maunder minimum, the reconstruction from the cosmogenic isotopes is consistent with that from the direct sunspot number observation. Furthermore, over the nineteenth century, the agreement of irradiance reconstructions using isotope records with the reconstruction from the sunspot number by Chatzistergos et al. (2017) is better than that with the reconstruction from the WDC-SILSO series (Clette et al. 2014), with a lower chi-square-value

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

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