Magnetic susceptibility dataset from the Danian Contessa Highway section, Gubbio, Italy

We studied a high-resolution multiproxy data set, including magnetic susceptibility (MS), CaCO3 content, and stable isotopes (d18O and d13C), from the stratigraphic interval covering the uppermost Maastrichtian and the lower Danian, represented by the pelagic limestones of the Scaglia Rossa Formation continuously exposed in the classic sections of the Bottaccione Gorge and the Contessa Highway near Gubbio, Italy. Variations in all the proxy series are periodic and reflect astronomically forced climate changes (i.e., Milankovitch cycles). In particular, the MS proxy reflects variations in the terrigenous dust input in this pelagic, deep-marine environment. We speculate that the dust is mainly eolian in origin and that the availability and transport of dust are influenced by variations in the vegetation cover on the Maastrichtian-Paleocene African or Asian zone, which were respectively located at tropical to subtropical latitudes to the south or far to the east of the western Tethyan Umbria-Marche Basin, and were characterized by monsoonal circulation. The dynamics of monsoonal circulation are known to be strongly dependent on precession-driven and obliquity-driven changes in insolation. We propose that a threshold mechanism in the vegetation coverage may explain eccentricity-related periodicities in the terrigenous eolian dust input. Other mechanisms, both oceanic and terrestrial, that depend on the precession amplitude modulated by eccentricity, can be evoked together with the variation of dust influx in the western Tethys to explain the detected eccentricity periodicity in the d13C record. Our interpretations of the d18O and MS records suggest a warming event ~400 k.y. prior to the Cretaceous-Paleogene (K-Pg) boundary, and a period of climatic and environmental instability in the earliest Danian. Based on these multiproxy phase relationships, we propose an astronomical tuning for these sections; this leads us to an estimate of the timing and duration of several late Maastrichtian and Danian biostratigraphic and magnetostratigraphic events.

Late Cretaceous to Eocene calcareous nannofossil distribution near Gubbio, Italy

Using a modified sample preparation technique, we have been able to establish a detailed lower Campanian to upper Eocene nannofossil stratigraphy in the Bottaccione and Contessa Highway sections near Gubbio. Appearance and extinction levels of virtually all the commonly used calcareous nannofossil zonal markers have been recognized and can now be closely correlated with the planktonic foraminifera zonation and the magnetic reversal stratigraphy previously established in these sections. Comparisons with the nannofossil calibrations of the oceanic magnetic anomaly sequence in Deep Sea Drilling Project (DSDP) sites suggest that magmetic Subchrons C17N and C25N are missing in the Bottaccione section. The observed variability of the relative stratigraphic position of most plankton events is confirmed to less than one magnetic subchron. Absolute abundance, paleobiogeographic restriction, and differential preservation render some of the traditionally used biostratigraphic events less reliable than others.

Natural remanent magnetization (NRM) from ODP Site 1263 covering magnetochron C20r and C21n and high-resolution bulk carbon isotope (d13C) records from Sites 702 and 1263

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.

Carbon isotope record of 3 profiles from Campanian-Maastrichtian

Carbon-isotope stratigraphy has proven to be a powerful tool in the global correlation of Cretaceous successions. Here we present new, high-resolution carbon-isotope records for the Global Boundary Stratotype Section and Point (GSSP) of the Maastrichtian stage at Tercis les Bains (France), the Bottaccione and Contessa sections at Gubbio (Italy), and the coastal sections at Norfolk (UK) to provide a global d13C correlation between shelf-sea and oceanic sites. The new d13C records are correlated with d13C-stratigraphies of the boreal chalk sea (Trunch borehole, Norfolk, UK, Lägerdorf-Kronsmoor-Hemmoor section, northern Germany, Stevns-1 core, Denmark), the tropical Pacific (ODP-Hole 1210B, Shatsky Rise) and the South Atlantic and Southern Ocean (DSDP Hole 525A, ODP Hole 690C) by using an assembled Gubbio d13C record as a reference curve. The global correlation allows the identification of significant high-frequency d13C variations that occur superimposed on prominent Campanian-Maastrichtian events, namely the Late Campanian Event (LCE), the Campanian-Maastrichtian Boundary Event (CMBE), the mid-Maastrichtian Event (MME), and the Cretaceous-Paleogene transition (KPgE). The carbon-isotope events are correlated with the geomagnetic polarity scale recalculated using the astronomical 40Ar/39Ar calibration of the Fish Canyon sanidine. This technique allows the evaluation of the relative timing of base occurrences of stratigraphic index fossils such as ammonites, planktonic foraminifera and calcareous nannofossils. Furthermore, the Campanian-Maastrichtian boundary, as defined in the stratotype at Tercis, can be precisely positioned relative to carbon-isotope stratigraphy and the geomagnetic polarity timescale. The average value for the age of the Campanian-Maastrichtian boundary is 72.1 ± 0.1 Ma, estimated by three independent approaches that utilize the Fish Canyon sanidine calibration and Option 2 of the Maastrichtian astronomical timescale. The CMBE covers a time span of 2.5 Myr and reflects changes in the global carbon cycle probably related to tectonic processes than to glacioeustasy. The duration of the high-frequency d13C variations instead coincides with the frequency band of long eccentricity, indicative of orbital forcing of changes in climate and the global carbon cycle.

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.

Magnetic suceptibility, carbon and oxygen stable isotope ratios and reflectance from the Bottaccione Gorge and Furlo section, Italy

The oceans at the time of the Cenomanian-Turonian transition were abruptly perturbed by a period of bottom-water anoxia. This led to the brief but widespread deposition of black organic-rich shales, such as the Livello Bonarelli in the Umbria-Marche Basin (Italy). Despite intensive studies, the origin and exact timing of this event are still debated. In this study, we assess leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world, by providing a 6-Myr-long astronomically-tuned timescale across the Cenomanian-Turonian boundary. We procure insights in the relationship between orbital forcing and the Late Cretaceous carbon cycle by deciphering the imprint of astronomical cycles on lithologic, geophysical, and stable isotope records, obtained from the Bottaccione, Contessa and Furlo sections in the Umbria-Marche Basin. The deposition of black shales and cherts, as well as the onset of oceanic anoxia, is related to maxima in the 405-kyr cycle of eccentricity-modulated precession. Correlation to radioisotopic ages from the Western Interior (USA) provides unprecedented age control for the studied Italian successions. The most likely tuned age for the Livello Bonarelli base is 94.17 ± 0.15 Ma (tuning #1); however, a 405-kyr older age cannot be excluded (tuning #2) due to uncertainties in stratigraphic correlation, radioisotopic dating, and orbital configuration. Our cyclostratigraphic framework suggests that the exact timing of major carbon cycle perturbations during the Cretaceous may be linked to increased variability in seasonality (i.e. a 405-kyr eccentricity maximum) after the prolonged avoidance of seasonal extremes (i.e. a 2.4-Myr eccentricity minimum). Volcanism is probably the ultimate driver of oceanic anoxia, but orbital periodicities determine the exact timing of carbon cycle perturbations in the Late Cretaceous. This unites two leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world.