Perturbations of the global carbon cycle across the Cretaceous–Palaeogene boundary

doi:10.21203/rs.3.rs-3893195/v1

Download PDF

Article

Perturbations of the global carbon cycle across the Cretaceous–Palaeogene boundary

https://doi.org/10.21203/rs.3.rs-3893195/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 10 May, 2024

Read the published version in Communications Earth & Environment →

Version 1

posted

You are reading this latest preprint version

Two hyperthermal events with different carbon cycle perturbations occurred across the Cretaceous-Palaeogene (K-Pg) boundary, i.e., the late Maastrichtian Warming Event (LMWE) and the early Danian Dan-C2 event. However, the roles played by Deccan volcanism and orbital forcing in these two hyperthermals are still debated. Here, we obtain a new terrestrial δ¹³C_carb record in the Nanxiong Basin (southeastern China) and compare it with marine records. The results show that both the LMWE and Dan-C2 event can be well distinguished in the terrestrial record and that the Dan-C2 event is characterized by a typical hyperthermal event; however, the specificity of the context under which this event occurred has resulted in inconsistencies in the marine records. In addition, the δ¹³C excursion during the LMWE was more muted and prolonged than that during the Dan-C2 event, and the short-eccentricity cycle disappeared in the marine record during the LMWE, indicating that Deccan volcanism perturbed the carbon cycle during the LMWE, while the Dan-C2 event was less influenced by volcanic perturbation.

Earth and environmental sciences/Climate sciences/Palaeoclimate

Earth and environmental sciences/Solid Earth sciences/Geology

Terrestrial record

Cretaceous-Palaeogene boundary

Global carbon cycle

Deccan Traps

Orbital forcing

A series of hyperthermal events characterized by transient warming and significant negative carbon isotope excursion (CIE) occurred during the late Cretaceous–early Palaeogene^1,2. Among them, the Late Maastrichtian Warming Event (LMWE) and the early Danian Dan-C2 event occurred across the Cretaceous-Palaeogene (K-Pg) boundary. Both were related to significant climatic and biotic changes; for instance, the LMWE was probably related to mass extinction at the end of the Cretaceous^3–5, while the Dan-C2 event promoted biotic recovery after the extinction^4,6.

The Dan-C2 event has been distinguished in the Atlantic Ocean^1,6–9, the Tethys Ocean^10–12, and the middle and lower latitudes of Eurasia^4,13–16. Compared to other hyperthermal events, the δ¹³C negative excursions during the Dan-C2 event were restricted to planktonic foraminifera and bulk records in parts of the Atlantic and Tethys Oceans, while benthic foraminifera rarely recorded this event^1,9. Furthermore, the warming indicated by oxygen isotopes (δ¹⁸O) during this event was also limited to surface waters in parts of the North Atlantic, with evidence of warming in bottom waters generally lacking ^1,9,17(supporting information Text S1, Fig. S1). Therefore, the global significance of the Dan-C2 event and even whether it should be regarded as a true hyperthermal are still controversial^1,7,17.

Two main sources of the massive ¹²C-rich CO₂ triggering of the carbon cycle during the LMWE and Dan-C2 event were proposed: Deccan volcanism and carbon pools controlled by orbital forcing. The extremely negative values of the CIE for the LMWE and Dan-C2 event are located at the maxima of the 405-kyr long-eccentricity cycle^18,4,12. Moreover, the two CIEs of the Dan-C2 event correspond to two 100-kyr short-eccentricity maxima¹, indicating that orbital forcing contributed to the CIE during both the LMWE and the Dan-C2 event. The overlapping occurrences of Deccan volcanism and two hyperthermals suggest that the large amount of volcanic carbon emitted could have also played a role in driving hyperthermals, especially the LMWE^3,5,6,18. However, this speculation is plausible and contemporaneous and lacks causal connection¹⁹. Recent modelling efforts to simulate CO₂ emission scenarios from Deccan volcanism have yielded conflicting results, with either more CO₂²⁰, half of the CO₂⁶, or one-third of the CO₂ ²¹released before the K-Pg boundary. Therefore, the roles that Deccan volcanism played in perturbing the carbon cycle in the LMWE and the Dan-C2 event are still unclear.

To date, most of the carbon isotope records are from marine sediments, whereas terrestrial records are lacking. Our previous work showed that the LMWE and Dan-C2 event were recorded by red strata in the Nanxiong Basin (Southeastern China, Fig. 1). Hg and its isotopes indicate that they formed during Deccan volcanism; moreover, both hyperthermals were related to the 405-kyr long eccentricity cycle^4,15,16. However, due to the low resolution of the δ¹³C_carb data, detailed comparisons with marine records and in-depth analyses of carbon cycle processes are limited. Here, a total of 274 fresh samples from the upper part of the Zhenshui Formation to the lower Xiahui part of the Shanghu Formation (Fig. 1b) were collected, then we obtain a new δ¹³C_carb record across the K-Pg boundary in the basin and compare it with other published records to reveal 1) the terrestrial records of the LMWE and Dan-C2 event, as well as the global significance of the Dan-C2 event; and 2) the relative contributions of Deccan volcanism and orbital forcing to the carbon cycle during the two hyperthermal events.

The δ¹³C value of the bulk carbonate from the Nanxiong Basin (Fig. 2c) is consistent with the previous δ¹³C values of pedogenic carbonate nodules that were not affected by diagenesis ^22,23(Fig. S2a); moreover, there is no significant correlation between δ¹³C_carb and δ¹⁸O_carb (Fig. S2b), suggesting that the influence of diagenesis on our bulk carbonate δ¹³C can be ruled out. δ¹³C_carb shows a steady decrease from 66.4 to 66.2 Ma, reaches a minimum at approximately 66.2 Ma, with a negative CIE of ~ 1.5‰, and then slowly increases to the K-Pg boundary, corresponding to the LMWE. δ¹³C_carb decreases sharply immediately after the K-Pg boundary and then shows an overall increase, with several negative excursions, such as the double CIEs at ~ 65.8 and ~ 65.7 Ma, which indicate that the Dan-C2 event characterized by the first CIE (~ 3‰) was larger than the second CIE (~ 2‰). In addition, another excursion of δ¹³C_carb (~ 2‰) occurred at ~ 65.3 Ma, which should correspond to the Lower C29n event. Overall, the δ¹³C_carb values in the Nanxiong Basin can be well compared with the δ¹³C_bulk values at ODP site 1262 (Fig. 2e); both show more muted CIEs during the LMWE than during the Dan-C2 event. In addition, the total duration of the onset, peak, and recovery of the LMWE (~ 200–300 kyr) is significantly longer than that of each of the CIEs during the Dan-C2 event (~ 100 kyr) (Figs. 2c and 2e). However, the magnitudes of the CIEs in the Nanxiong Basin (1.5-3‰) are greater than those in ODP Site 1262 (0.5-1‰), which is also consistent with the findings of previous studies showing that the CIE magnitude is greater in the terrestrial record than in the marine record for the same hyperthermal event^24,25.

Evolutionary power spectral analysis revealed that both the δ¹³C_carb of the Nanxiong Basin and the δ¹³C_bulk of ODP site 1262 exhibited significant 405-kyr long eccentricity cycles from 66.6 to 65 Ma (Fig. 3a). The dramatic negative excursion of global δ¹³C immediately after the K-Pg boundary was probably related to the large amount of CO₂ injected into the atmosphere by the Chicxulub impact³⁰, as well as the vital effects caused by mass extinction³¹. This dramatic excursion could lead to the insignificance of the 100-kyr short eccentricity cycle in both terrestrial and marine records under a 500 kyr sliding window during 66.2 to 65.8 Ma (Fig. 3a). To eliminate this bias, two discrete time windows without dramatic excursions were divided and then analysed separately with a sliding window of 150 kyr. The results reveal significant 100-kyr short eccentricity cycles both below and above the K-Pg boundary in the Nanxiong Basin (Fig. 3b); however, this cycle is complicated at ODP site 1262: the short eccentricity cycle was significant from 66.6 to 66.3 Ma, after which the significance level decreased and even disappeared from 66.2–66.0 Ma (Fig. 3c, right), after which the significance returned again after the K-Pg event (Fig. 3c, left).

Was the Dan-C2 event a global hyperthermal event?

Our study provides a new terrestrial δ¹³C record of the LMWE and Dan-C2 event, expanding their global distribution, especially the Dan-C2 event. However, there are still inconsistencies in marine records of the Dan-C2 event (supporting information, Text S1, Fig. S1), which are summarized in the following aspects: 1) the negative excursion of δ¹³C is mainly recorded by bulk samples and planktonic foraminifera, while benthic foraminifera are rarely recorded; 2) the distribution is restricted to parts of the Atlantic Ocean and the Tethys Ocean; and 3) evidence of bottom-seawater warming is generally lacking.

To explore the possible reasons for the inconsistencies in the Dan-C2 event, we compared it to the Palaeocene-Eocene Thermal Maximum (PETM), which is the most studied hyperthermal event. The CIE of the PETM can reach 7‰ in terrestrial records and only 3‰ in marine records^24,25,32. Two mechanisms have been proposed: 1) marine CIE records are commonly truncated by carbonate dissolution, so the most extreme values are not represented, leading to an incomplete CIE^32,33; and 2) the terrestrial CIE is amplified by environmental changes and fractionation effects caused by photosynthesis in plants³⁴. A similar phenomenon is also shown with the Dan-C2 event. The CIEs of δ¹³C_organic and δ¹³C_carb in Boltysh crater (Fig. S1l) and the Nanxiong Basin (Fig. 1c) can reach ~ 3.0‰, whereas the marine CIEs are less than 1.5‰ (Fig. S1). Differences in the carbon isotope records during the PETM are not only manifested in marine and terrestrial records but also among different types of carbonate records (bulk samples, planktonic foraminifera, or benthic foraminifera) and even within the same type of carbonate but in different marine regions³³ (supporting information Fig. S3 and Fig. 4). Ocean acidification during the PETM led to carbonate dissolution and shoaling of the carbonate compensation depth (CCD), causing the sediments to be clay-rich. The thickness of the clay layer increases with increasing palaeodepth^32,33, indicating enhanced carbonate dissolution. The CIEs of δ¹³C_bulk and δ¹³C_benthic have significant negative correlations with palaeodepth (Figs. 4a and 4b), suggesting that the dissolution of carbonate suppressed the amplitude of the CIE. In addition to carbonate dissolution, the sedimentation rate also affects the CIE. The CIEs of δ¹³C_bulk and δ¹³C_benthic have positive correlations with sediment thickness (Figs. 4c and 4d). The greater the sediment thickness is, the greater the sedimentation rate, the more complete the carbonate isotope record, and the greater the CIE. Similarly, for the Dan-C2 event, the CIEs of δ¹³C_bulk show a significant negative correlation with palaeodepth (Fig. 4e) but a positive correlation with sediment thickness (Fig. 4f) and carbonate concentration (Fig. 4g), indicating that the CIEs of the Dan-C2 event are also influenced by the water depth, sedimentation rate, and carbonate dissolution, leading to some differences in the CIE records from different regions.

However, compared with the PETM, the Dan-C2 event was a short-lived and muted warming event that occurred under a specific environmental background, for instance, biotic turnover and drastic ecosystem changes caused by mass extinction. Under normal oceanic conditions, phytoplankton preferentially convert ¹²CO₂ into organic matter through photosynthesis (primary productivity), and then, ¹²C-rich organic carbon is transported to the bottom water through a biological pump (export productivity), thus promoting carbon exchange between the surface and the deep ocean³⁵. After the K-Pg boundary, a collapse occurred in the surface-bottom δ¹³C gradient³⁶, which was initially thought to reflect either a collapse in primary productivity (“Strangelove Ocean”³⁶) or export productivity (“Living Ocean”³⁷) after the mass extinction. However, benthic faunal records show a lack of significant extinction of phytoplankton-dependent benthic foraminifera and an increased food flux to the seafloor in the southeastern Atlantic Ocean and the Pacific Ocean³¹. In addition, biogenic barium records indicate geographic heterogeneity in export productivity, an increase in the central Pacific Ocean, no changes in upwelling or shelf Atlantic sites, and decreases in the northeast and southwest Atlantic Ocean, Southern Ocean, and Indian Ocean³⁸. Thus, the “Resilient Ocean” and “Heterogeneous Ocean” have been proposed to explain these phenomena³⁹.

The spatial heterogeneity of primary/export productivity could cause heterogeneous carbon cycling processes between the surface and deep ocean, leading to inconsistent marine δ¹³C records from the Dan-C2 event. However, this heterogeneity was proposed to be due to the limited number of sites, which is insufficient to reveal a robust pattern, and the mechanism responsible for this heterogeneity is still unclear³⁸. The geographic location⁴⁰, circulation, nutrient runoff from land, and stratification³⁸ are potential drivers of the spatially heterogeneous ocean and inconsistent records of the Dan-C2 event. For example, ODP site 1049C is located at a productive coastal upwelling site⁴¹ in the western North Atlantic, where a great deal of terrestrial materials are transported by surface ocean currents⁴²; in addition, the export productivity is stable⁴³, leading to strong negative excursions from all δ¹³C_bulk, δ¹³C_planktic and δ¹³C_benthic (Fig. S1a), especially negative excursions from δ¹³C_benthic. Although increased export productivity was recorded at ODP sites 1209 and 1210 in the Central Pacific Ocean^24,38, there are no records of Dan-C2 event, probably due to pelagic and oligotrophic environments, as well as the lack of nutrients from terrestrial sources. Evidently, the mass extinction and the dramatic changes it brought are the indispensable causes of the spatial heterogeneity of the Dan-C2 records, but proposing a reasonable model to explain this mechanism requires additional records in the future.

The roles of Deccan volcanism and orbital forcing in the carbon cycle

The large amount of CO₂ released by large igneous provinces (LIPs) can cause perturbations to the global carbon cycle. In addition to LIPs, other carbon pools, such as peat and methane hydrates, which are modulated by eccentricity forcing, could also contribute to the carbon cycle during hyperthermals. For instance, during eccentricity minima, seasonally uniform annual precipitation is more suitable for carbon burial, whereas during eccentricity maxima, short wet seasons and prolonged dry seasons caused by “monsoon-like” precipitation could promote carbon release². In addition, methane hydrates buried in the marine shelf become unstable and decompose in response to orbital-driven warming, leading to large quantities of light carbon being emitted to the atmosphere-ocean system, further perturbing the global carbon cycle⁴⁴. Previous work has shown that the total mercury (Hg) concentration in the Nanxiong Basin has been anomalous from 66.4 to 65.6 Ma (Fig. 2b); combined with Hg isotope data, Ma et al.⁴ attributed these anomalies to volcanism in the central Deccan Traps. Both the LMWE and the Dan-C2 event temporally overlapped with the central Deccan volcanism (Fig. 2). Moreover, the LMWE and Dan-C2 event occurred during the last 405-kyr long eccentricity of the Maastrichtian and the first 405-kyr long eccentricity of the Danian, respectively (Figs. 2c, 2d, and 2e), and their CIE maxima were all within the maxima of the 405-kyr eccentricity cycle¹². Moreover, 100-kyr short eccentricity cycles were significant in both terrestrial and marine records (except for the LMWE of ODP 1262; Fig. 3). These findings imply that both Deccan volcanism and orbital forcing contributed to the LMWE and Dan-C2 event.

However, there are noticeable differences between the LMWE and Dan-C2 event: 1) the magnitude of the CIE during the LMWE (~ 1.5‰ in the Nanxiong Basin,<0.5‰ in ODP 1262) is more muted than in the Dan-C2 event (2–3‰ in the Nanxiong Basin, 0.6‰ in ODP 1262); 2) the LMWE was characterized by both surface and deep sea warming^1,18, while the Dan-C2 event was characterized by surface ocean warming, accompanied by little appreciable deep sea warming; 3) each of the double CIEs of the Dan-C2 event corresponds to a maxima of the 100-kyr eccentricity cycle, while the total duration of the onset, peak and recovery of the LMWE (200 ~ 300 kyr) was significantly longer than each CIE of the Dan-C2 event (Fig. 2); 4) the short eccentricity cycles are significant during the Dan-C2 event, whereas they are insignificant and even disappear during the LMWE in the marine record (Fig. 3); and 5) although the CIE maxima of both the Dan-C2 and the LMWE were all within the maxima of the 405-kyr eccentricity cycle, the onset of the LMWE occurred at the minima of the 405-kyr eccentricity cycle^12,18 (Figs. 2c and 2e). These apparent differences suggest that Deccan volcanism and orbital forcing played different roles in driving the LMWE and Dan-C2 event, as well as in the global carbon cycle.

High-precision chronologies indicate that both the eruption rate and volume were low during the early stage of the central Deccan Traps^26,27 (Fig. 2a). However, CO₂ release has the potential to decouple from rates of surface volcanism because large amounts of CO₂ can be released through passive degassing^6,27, especially from intrusive magmas²⁰. The reconstructed atmospheric CO₂ concentration based on the pedogenic carbonate nodules showed higher pCO₂ values during the LMWE than during the Dan-C2 event^5,45 (Fig. S4), which is consistent with direct measurements of melt-inclusion CO₂ concentrations, suggesting that early Deccan magmas were enriched with more CO₂²⁰. Although the onset of the LMWE occurred at the minima of the 405-kyr eccentricity¹⁸ (Figs. 2c and 2e), several thousand Gt of carbon that degassed from the early Deccan magmas was sufficient to trigger the LMWE^6,20. The δ¹³C composition of volcanic CO₂ (~-5‰⁴⁶) is much more positive than that of other sources, such as peat (δ¹³C≈-25‰⁴⁷) and methane hydrates (δ¹³C ≈-60‰⁴⁸); therefore, the massive amount of volcanic CO₂ emitted through passive degassing of early Deccan magmas could have led to muted and prolonged δ¹³C negative excursions during the LMWE (Fig. 2), as well as disruption of the short-eccentricity cycle in the oceanic record (Fig. 3c, right). Although passive degassing triggered the LMWE, whether the amount of released CO₂ was sufficient to cause ~ 4°C of warming is still debated^19,20. Notably, 405-kyr long eccentricity cycles were significant during the LMWE according to both terrestrial and marine records (Fig. 3a), and even 100-kyr short eccentricity cycles were significant according to the terrestrial record (Fig. 3b, right), suggesting that Deccan CO₂ outgassing likely enhanced the climate sensitivity to orbital forcing, leading to a global warming of ~ 4°C^12,20. After the K-Pg boundary, with a decrease in CO₂ released from the Deccan magma^20,27, the carbon cycle was mainly controlled by orbitally driven carbon pools with more negative δ¹³C values, leading to larger CIEs (Fig. 2) and significant short-eccentricity cycles in both terrestrial and marine records (Figs. 3b left and 3c left). Due to the decrease in CO₂ emissions, pCO₂ and its warming effects were muted during the Dan-C2 event. The above speculated release scenario of Deccan volcanic CO₂ was further substantiated by long-term ocean-atmosphere-sediment carbon cycle reservoir (LOSCAR)⁴⁹ model simulations, which showed that prior to the K-Pg boundary, more CO₂ was released through passive degassing²⁰ (intrusive:extrusive = 5:1) or that half of the CO₂ was released (50:50 outgassing scenario) but with a higher emission rate⁶. The greater volume and higher rate of CO₂ release before compared to after the K-Pg boundary confirm that Deccan volcanism likely contributed to both the LMWE and Dan-C2 event but contributed more to the LMWE.

In conclusion, we provide a new terrestrial δ¹³C record of the LMWE and Dan-C2 event in low-latitude East Asia, which can be compared with marine records, further expanding the global distribution of these events. The inconsistency of marine records for the Dan-C2 event is related to the drastic ecosystem changes caused by mass extinction, especially in the heterogeneous ocean, while the specific mechanism remains to be revealed by additional studies. We hypothesize that Deccan volcanism and orbital forcing played different roles in the carbon cycles during the LMWE and Dan-C2 event. Deccan volcanic CO₂ triggered the LMWE through passive degassing, disturbed the carbon cycle and amplified the sensitivity of the climate to orbital forcing, whereas the Dan-C2 event was mainly controlled by orbital forcing, with weakening of the volcanic perturbation.

The Nanxiong Basin is in southeastern China (Fig. 1a); it is elongated, with its axis oriented from northeast to southwest. Continuous red fluvial–lacustrine clastics spanning from the Upper Cretaceous to the lower Palaeocene are preserved in the basin. Extensive chronostratigraphic, stratigraphic, palaeontological, and palaeoclimatic works have been carried out in the CGD-CGY section (also called the Datang section)^15,16,22,23, which is dominated by muddy siltstone and silty mudstone with interbedded sandstone and conglomerate (Fig. 1b). Many fossils have been preserved, for instance, dinosaurs, dinosaur footprints, dinosaur eggs, and mammals. Various palaeosol layers with pedogenic carbonate formed on the red clastics²². Previous studies have shown that the palaeoclimate in this basin across the K-Pg boundary was mainly hot and (semi)arid; moreover, the climatic evolution was consistent with marine records^4,15. A total of 274 fresh samples from the upper part of the Zhenshui Formation to the lower Xiahui part of the Shanghu Formation (CGD section; Fig. 1b) were collected at approximately 1 m intervals. The chronology of the CGD section follows that of Ma et al.⁴ and is based on the palaeomagnetic framework of Clyde et al.²² combined with a chronological control point on the K-Pg boundary.

The carbon isotopic composition of the bulk carbonate was determined using a Thermo Fisher Isotope ratio mass spectrometer (Mat 253) coupled with a GasBench II at the Laboratory for Stable Isotope Geochemistry, Institute of Geology and Geophysics, IGGCAS, through the production of CO₂ after reaction with phosphoric acid. Acid digestion was performed in a GasBench II in continuous flow mode at a temperature of 72 ± 0.1°C and a reaction time of 60 min, through which the generated CO₂ was transferred by high-purity (99.999%) He carrier gas into the mass spectrometer. The standard deviation of δ¹³C values was calculated from replicate analyses of an internal laboratory calcite standard, which is better than 0.15‰. The measured δ¹³C values are reported relative to those of the Vienna Pee Dee Belemnite (V-PDB).

To determine the contribution of orbital forcing to the carbon cycle, evolutionary power spectra were generated through the δ¹³C records of both the Nanxiong Basin and ODP 1262 using Acycle software⁵⁰. In preparation for this analysis, the δ¹³C records were interpolated linearly at 5-kyr and 2-kyr intervals for the Nanxiong Basin and ODP 1262 records, respectively, and detrended using local regression smoothing (LOWESS). The fast Fourier transform (FFT) method was selected, and the sliding windows were 500 kyr and 150 kyr for the complete and discrete time windows, respectively.

Competing interests:

The authors declare no competing interests.

Author contribution:

Mingming Ma and Xiuming Liu designed this study; Mingming Ma and Mengdi Wang carried them out; Huixin Huang conducted the time-series analyses; Mingming Ma prepared the manuscript with contributions from all co-authors.

Acknowledgement:

This research was supported by the National Natural Science Foundation of China (awards 42277440 and 42130507), and the IGCP Project 679.

Barnet, J. S. K. et al. A High-Fidelity Benthic Stable Isotope Record of Late Cretaceous–Early Eocene Climate Change and Carbon‐Cycling. Paleoceanogr. Paleoclimatology 34, 672–691 (2019).
Zachos, J. C., Mccarren, H., Murphy, B., Röhl, U. & Westerhold, T. Tempo and scale of late Paleocene and early Eocene carbon isotope cycles: Implications for the origin of hyperthermals. Earth Planet. Sci. Lett. 299, 242–249 (2010).
Li, L. & Keller, G. Abrupt deep-sea warming at the end of the Cretaceous. Geology 26, 995–998 (1998).
Ma, M. et al. Deccan Traps Volcanism Implicated in the Extinction of Non-Avian Dinosaurs in Southeastern China. Geophys. Res. Lett. 49, e2022GL100342 (2022).
Zhang, L. et al. Deccan volcanism caused coupled pCO2 and terrestrial temperature rises, and pre-impact extinctions in northern China. Geology 46, 271–274 (2018).
Hull, P. M. et al. On impact and volcanism across the Cretaceous-Paleogene boundary. Science 367, 266–272 (2020).
Arreguín-Rodríguez, G. J. et al. Benthic foraminiferal turnover across the Dan-C2 event in the eastern South Atlantic Ocean (ODP Site 1262). Palaeogeogr. Palaeoclimatol. Palaeoecol. 572, 110410 (2021).
Krahl, G. et al. Environmental changes occurred during the Early Danian at the Rio Grande Rise, South Atlantic Ocean. Glob. Planet. Change 191, 103197 (2020).
Quillévéré, F., Norris, R. D., Kroon, D. & Wilson, P. A. Transient ocean warming and shifts in carbon reservoirs during the early Danian. Earth Planet. Sci. Lett. 265, 600–615 (2008).
Coccioni, R. et al. The Dan-C2 hyperthermal event at Gubbio (Italy): Global implications, environmental effects, and cause(s). Earth Planet. Sci. Lett. 297, 298–305 (2010).
Gilabert, V., Batenburg, S. J., Arenillas, I. & Arz, J. A. Contribution of orbital forcing and Deccan volcanism to global climatic and biotic changes across the Cretaceous-Paleogene boundary at Zumaia, Spain. Geology 50, 21–25 (2022).
Gilabert, V., Arenillas, I., Arz, J. A., Batenburg, S. J. & Robinson, S. A. Multiproxy analysis of paleoenvironmental, paleoclimatic and paleoceanographic changes during the early Danian in the Caravaca section (Spain). Palaeogeogr. Palaeoclimatol. Palaeoecol. 576, 110513 (2021).
Gilmour, I. et al. A high-resolution nonmarine record of an early Danian hyperthermal event, Boltysh crater, Ukraine. Geology 41, 783–786 (2013).
Jolley, D. W., Gilmour, I., Gilmour, M., Kemp, D. B. & Kelley, S. P. Long-term resilience decline in plant ecosystems across the danian dan-c2 hyperthermal event, boltysh crater, ukraine. J. Geol. Soc. 172, 491–498 (2015).
Ma, M., He, M., Zhao, M., Peng, C. & Liu, X. Evolution of atmospheric circulation across the Cretaceous–Paleogene (K–Pg) boundary interval in low-latitude East Asia. Glob. Planet. Change 199, 103435 (2021).
Zhao, M., Ma, M., He, M., Qiu, Y. & Liu, X. Evaluation of the four potential Cretaceous-Paleogene (K-Pg) boundaries in the Nanxiong Basin based on evidences from volcanic activity and paleoclimatic evolution. Sci. China Earth Sci. 64, 631–641 (2021).
Westerhold, T., Röhl, U., Donner, B., McCarren, H. K. & Zachos, J. C. A complete high-resolution Paleocene benthic stable isotope record for the central Pacific (ODP Site 1209). Paleoceanography 26, PA002092 (2011).
Barnet, J. S. K., Littler, K., Kroon, D., Leng, M. J. & Zachos, J. C. A new high-resolution chronology for the late Maastrichtian warming event: Establishing robust temporal links with the onset of Deccan volcanism. Geology 46, 147–150 (2017).
Tobin, T. S., Bitz, C. M. & Archer, D. Modeling climatic effects of carbon dioxide emissions from Deccan Traps volcanic eruptions around the Cretaceous–Paleogene boundary. Palaeogeogr. Palaeoclimatol. Palaeoecol. 478, 139–148 (2017).
Nava, A. H. et al. Reconciling early Deccan Traps CO₂ outgassing and pre-KPB global climate. Proc. Natl. Acad. Sci. 118, e2007797118 (2021).
Cox, A. A. & Keller, C. B. A Bayesian inversion for emissions and export productivity across the end-Cretaceous boundary. Science 381, 1446–1451 (2023).
Clyde, W. C. et al. New Paleomagnetic and Stable-Isotope Results from the Nanxiong Basin, China: Implications for the K/T Boundary and the Timing of Paleocene Mammalian Turnover. J. Geol. 118, 131–143 (2010).
Wang Y. Paleoclimate changes of the late Cretaceous-early Paleocene in the Nanxiong basin, south China (in Chinese with English abstract). (Nanjing University, Nangjing, 2012).
Abels, H. A. et al. Terrestrial carbon isotope excursions and biotic change during Palaeogene hyperthermals. Nat. Geosci. 5, 326–329 (2012).
Chen, Z. et al. Structure of the carbon isotope excursion in a high-resolution lacustrine Paleocene–Eocene Thermal Maximum record from central China. Earth Planet. Sci. Lett. 408, 331–340 (2014).
Schoene, B. et al. U-Pb constraints on pulsed eruption of the Deccan Traps across the end-Cretaceous mass extinction. Science 363, 862–866 (2019).
Sprain, C. J. et al. The eruptive tempo of Deccan volcanism in relation to the Cretaceous-Paleogene boundary. Science 363, 866–870 (2019).
Laskar, J., Fienga, A., Gastineau, M. & Manche, H. La2010: a new orbital solution for the long-term motion of the Earth. Astron. Astrophys. 532, 784–785 (2011).
Kroon, D. & Zachos, J. C. Leg 208 synthesis: cenozoic climate cycles and excursions. Sci. Results 208, 1–55 (2007).
Artemieva, N., Morgan, J., & Expedition 364 Science Party. Quantifying the Release of Climate-Active Gases by Large Meteorite Impacts With a Case Study of Chicxulub. Geophys. Res. Lett. 44, 10180–10188 (2017).
Alegret, L., Thomas, E. & Lohmann, K. C. End-Cretaceous marine mass extinction not caused by productivity collapse. Proc. Natl. Acad. Sci. U.S.A. 109, 728–32 (2012).
Zachos, J. C. et al. Rapid Acidification of the Ocean During the Paleocene-Eocene Thermal Maximum. Science 308, 1611–1615 (2005).
McCarren, H., Thomas, E., Hasegawa, T., Röhl, U. & Zachos, J. C. Depth dependency of the Paleocene-Eocene carbon isotope excursion: Paired benthic and terrestrial biomarker records (Ocean Drilling Program Leg 208, Walvis Ridge). Geochem. Geophys. Geosyst. 9, Q10008 (2008).
Bowen, G. J., Beerling, D. J., Koch, P. L., Zachos, J. C. & Quattlebaum, T. A humid climate state during the Palaeocene/Eocene thermal maximum. Nature 432, 495–499 (2004).
Hilting, A. K., Kump, L. R. & Bralower, T. J. Variations in the oceanic vertical carbon isotope gradient and their implications for the Paleocene-Eocene biological pump. Paleoceanogr. Paleoclimatology 23, PA3222 (2008).
Hsü, K. J. et al. Environmental and evolutionary consequences of mass-mortality at the end of the Cretaceous. Science 216, 249–256 (1982).
D’Hondt, S. Consequences of the Cretaceous/Paleogene mass extinction for marine ecosystems. Annu. Rev. Ecol. Evol. Syst. 36, 295–317 (2005).
Hull, P. M. & Norris, R. D. Diverse patterns of ocean export productivity change across the Cretaceous-Paleogene boundary: New insights from biogenic barium. Paleoceanography 26, PA3205 (2011).
Esmeray-Senlet, S. et al. Evidence for reduced export productivity following the Cretaceous/Paleogene mass extinction. Paleoceanography 30, 718–738 (2015).
Jiang, S., Bralower, T. J., Patzkowsky, M. E., Kump, L. R. & Schueth, J. D. Geographic controls on nannoplankton extinction across the Cretaceous/Palaeogene boundary. Nat. Geosci. 3, 280–285 (2010).
Alegret, L. & Thomas, E. Cretaceous/Paleogene boundary bathyal paleo-environments in the central North Pacific (DSDP Site 465), the Northwestern Atlantic (ODP Site 1049), the Gulf of Mexico and the Tethys: The benthic foraminiferal record. Palaeogeogr. Palaeoclimatol. Palaeoecol. 224, 53–82 (2005).
Nauter-Alves, A. et al. Biotic turnover and carbon cycle dynamics in the early Danian event (Dan-C2): New insights from Blake Nose, North Atlantic. Glob. Planet. Change 221, 104046 (2023).
Faul, K. L., Anderson, L. D. & Delaney, M. L. Late Cretaceous and early Paleogene nutrient and paleoproductivity records from Blake Nose, western North Atlantic Ocean. Paleoceanography 18, PA000722 (2003).
Lunt, D. J. et al. A model for orbital pacing of methane hydrate destabilization during the Palaeogene. Nat. Geosci. 4, 775–778 (2011).
Nordt, L. C., Atchley, S. & Dworkin, S. Terrestrial evidence for two greenhouse events in the latest Cretaceous. Geol. Soc. Am. 13, 4–9 (2003).
Grard, A., Francois, L. M., Dessert, C., B. Dupré & Y. Goddéris. Basaltic volcanism and mass extinction at the Permo-Triassic boundary: Environmental impact and modeling of the global carbon cycle. Earth Planet. Sci. Lett. 234, 207–221 (2005).
DeConto, R. M. et al. Past extreme warming events linked to massive carbon release from thawing permafrost. Nature 484, 87–91 (2012).
Dickens, G. R., Castillo, M. M. & Walker, J. C. A blast of gas in the latest Paleocene: Simulating first-order effects of massive dissociation of oceanic methane hydrate. Geology 25, 259–262 (1997).
Zeebe, R. E. LOSCAR: Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir Model v2.0.4. Geosci. Model Dev. 5, 149–166 (2012).
Li, M., Hinnov, L. & Kump, L. Acycle: Time-series analysis software for paleoclimate research and education. Comput. Geosci. 127, 12–22 (2019).

There is NO Competing Interest.

SI.docx

Download PDF

Journal Publication

published 10 May, 2024

Read the published version in Communications Earth & Environment →

Version 1

posted

You are reading this latest preprint version

Perturbations of the global carbon cycle across the Cretaceous–Palaeogene boundary

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Results

Discussion

Was the Dan-C2 event a global hyperthermal event?

The roles of Deccan volcanism and orbital forcing in the carbon cycle

Materials and Methods

Declarations

Competing interests:

Author contribution:

Acknowledgement:

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1