Middle Eocene greenhouse warming facilitated by diminished weathering feedback
Robin van der Ploeg ORCID: orcid.org/0000-0001-5096-74961,
David Selby2,3,
Margot J. Cramwinckel ORCID: orcid.org/0000-0002-6063-836X1,
Yang Li ORCID: orcid.org/0000-0001-7441-77124,5,
Steven M. Bohaty ORCID: orcid.org/0000-0002-1193-73986,
Jack J. Middelburg ORCID: orcid.org/0000-0003-3601-90721 &
Appy Sluijs1
Nature Communications volume 9, Article number: 2877 (2018) Cite this article
Climate and Earth system modelling
Element cycles
Palaeoceanography
The Middle Eocene Climatic Optimum (MECO) represents a ~500-kyr period of global warming ~40 million years ago and is associated with a rise in atmospheric CO2 concentrations, but the cause of this CO2 rise remains enigmatic. Here we show, based on osmium isotope ratios (187Os/188Os) of marine sediments and published records of the carbonate compensation depth (CCD), that the continental silicate weathering response to the inferred CO2 rise and warming was strongly diminished during the MECO—in contrast to expectations from the silicate weathering thermostat hypothesis. We surmise that global early and middle Eocene warmth gradually diminished the weatherability of continental rocks and hence the strength of the silicate weathering feedback, allowing for the prolonged accumulation of volcanic CO2 in the oceans and atmosphere during the MECO. These results are supported by carbon cycle modeling simulations, which highlight the fundamental importance of a variable weathering feedback strength in climate and carbon cycle interactions in Earth's history.
The chemical weathering of silicate rocks represents a negative feedback mechanism that is generally considered to modulate atmospheric CO2 levels and Earth's climate on geological timescales1,2. This phenomenon has been studied for various carbon cycle perturbations and episodes of global warming in the geological past, including Pleistocene deglaciations, the Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma), and the Cretaceous and Jurassic Oceanic Anoxic Events (OAEs), mainly through the application of isotope ratios of marine sediments that are sensitive to shifts in weathering fluxes or compositions on the appropriate timescales3,4,5. For many of these phases, it is now relatively well established that enhanced continental weathering contributed to CO2 drawdown and climatic recovery4,6,7. However, the available data spanning the Middle Eocene Climatic Optimum (MECO; ~40 Ma) pose questions regarding the functioning of the weathering feedback8. Over a period of ~500 kyr, global ocean temperatures rose gradually by up to ~5 °C in association with an increase in atmospheric CO2 concentrations, sourced from a reservoir with a stable carbon isotopic composition (δ13C) close to that of the ocean9,10,11,12,13. Importantly, the inferred rise in atmospheric CO2 and temperature over ~500 kyr during the MECO should have led to increased weathering and alkalinity supply to the oceans, but reconstructions show that the oceans acidified8,10. Therefore, reconstructing the global weathering response during the MECO is instrumental to improving our fundamental understanding of carbon cycle dynamics on such intermediate timescales of ~500 kyr8.
A promising proxy to reconstruct changes in continental weathering during the MECO is the osmium isotope ratio of marine sediments at the time of deposition (187Os/188Osinitial, or Osi)14,15. The 187Os/188Os ratio of the global ocean is governed by the relative input of radiogenic Os (187Os/188Os = ~1.4) through continental weathering of ancient crustal rocks, and the relative input of unradiogenic Os (187Os/188Os = 0.13) through hydrothermal activity at mid-ocean ridges and weathering of fresh mantle-derived rocks, with additional contributions from extraterrestrial sources14. Osmium is a quasi-conservative element that is well-mixed in the ocean, and has a short oceanic residence time (generally ~104 yr in the open ocean, but residence times of ~103 yr have been inferred for very restricted settings)14,16. Variations in the 187Os/188Os ratio of seawater are thus indicative of changes in continental weathering relative to the other sources on timescales shorter than, or similar to, climate and carbon cycle processes such as greenhouse warming, ocean acidification, and carbonate compensation14,15. Seawater Os is incorporated in the metalliferous and organic phases of marine sediments without isotopic fractionation, and remains a closed isotopic system from the time of deposition17,18,19. As such, Osi values are calculated on the basis that radiogenic 187Os ingrowth is derived solely from post-depositional 187Re (rhenium) decay. Shifts to higher (radiogenic) Osi values, which are attributed to a global increase in continental silicate weathering rates, have been recorded for carbon cycle perturbations such as the Toarcian OAE and the PETM and Eocene Thermal Maximum 2 (ETM2) transient global warming events5,15,20.
A second parameter that is often used to constrain changes in continental weathering is the carbonate compensation depth (CCD). The CCD is the depth in the oceans at which carbonate delivery is balanced by carbonate dissolution, and is modulated by the interplay of volcanic CO2 degassing, the weathering of silicate rocks and organic-rich sediments on land, and the burial of marine carbonates and organic carbon21. As such, changes in the position of the CCD as reflected in sediments play a crucial role in reconstructions of carbon cycle change, both on multi-million year timescales and during transient perturbations such as the MECO22.
In this study, we present Osi records of marine sediments from three locations in different ocean basins in combination with a compilation of published CCD records8 to reconstruct global changes in continental weathering during the MECO. Rather than an Osi increase expected from globally enhanced weathering, we document a modest global Osi decrease during the MECO that may represent an episode of enhanced volcanism and/or associated basalt weathering. In fact, prolonged CCD shoaling precludes an increase in total continental weathering rates in response to CO2 rise and greenhouse warming. We employ a series of simulations with the carbon cycle model LOSCAR23 together with an independent osmium cycle model to demonstrate that this combination of observations can only be successfully reconciled on MECO timescales by invoking enhanced volcanism together with a diminished continental weathering feedback. Finally, we surmise that such a reduced silicate weathering feedback may have resulted from a progressive decrease in the weatherability of the continents during the Eocene. A variable silicate weathering feedback strength may have been important for other enigmatic climate and carbon cycle perturbations in Earth's history.
Middle Eocene osmium isotope records
We present Re-Os data and Osi values for middle Eocene sediments from Ocean Drilling Program (ODP) Site 959 in the equatorial Atlantic along the African continental margin, ODP Site 1263 on the Walvis Ridge in the South Atlantic, and Integrated Ocean Drilling Program (IODP) Site U1333 in the equatorial Pacific (Fig. 1; Supplementary Data 1; Supplementary Figs. 1–3). The Re and Os abundances are significantly enriched in the relatively organic-rich, siliceous sediments of Site 959 (Re = 10–60 ppb, Os = 100–300 ppt) relative to the carbonate-rich pelagic sediments of Sites 1263 and U1333 (Re = 0.02–0.2 ppb, Os = 10–40 ppt). The abundances of 192Os, the Os isotope best representing the amount of hydrogenous Os chelated by organic matter at the time of deposition24, increase slightly over the study interval at Site 959, but are essentially stable at the other two sites (Fig. 1). We calculate Osi values of 0.46 to 0.60 at all study sites (Fig. 1), which is in good agreement with previously published middle Eocene Osi values from Site 959 sediments25,26 and with Osi values from ferromanganese crusts that document a progressive increase in the 187Os/188Os composition of seawater during the Cenozoic27,28,29 (Fig. 2).
Osi values (in blue) and 192Os concentrations (in red) for the analyzed middle Eocene sediments from the three different sites. a ODP Site 959; b ODP Site 1263; c IODP Site U1333. The MECO interval is defined based on TEX86 values for Site 959 (in black; Cramwinckel et al.13) and bulk carbonate stable oxygen isotope ratios (δ18O) for Site 1263 (in black; Bohaty et al.10). The MECO is characterized by low carbonate content at Site U1333 (in grey; Westerhold et al.84). The error bars indicate fully propagated analytical uncertainties (2σ)
Comparison of Osi records from the MECO with Osi records from the PETM and ETM2, shown against the overall Osi evolution of the Cenozoic and the relative weathering feedback strength of the Cenozoic. a MECO data from Site 959 (in red), Site 1263 (in blue) and Site U1333 (in green) plotted against age (GTS2012)74. See Methods for discussion of the age models for the study sites. b MECO data from Sites 959, 1263, and U1333 (this study); PETM and ETM2 data from DSDP Site 549 (in purple) as published in Peucker-Ehrenbrink & Ravizza15; Cenozoic data from ferromanganese crusts D11 and CD29 (in black) as published in Klemm et al.27 and Burton28, respectively, based on the updated age model of Nielsen et al.29. c Model estimates of the relative continental weathering feedback strength of the Cenozoic as published in Caves et al.57, based on their CO2 scenario 1 and a logarithmic expression for the weathering feedback
At Site 959, the Osi values range between approximately 0.56 and 0.60 for most of the middle Eocene study interval, with the exception of a decrease to 0.51 during the MECO at ~580 mbsf (Fig. 1). Importantly, the lack of an increase in the Osi values during the MECO implies that weathering rates of felsic silicate rocks did not increase in response to CO2 rise and accompanied warming, while such an increase would be expected from theory and published Osi records from analogous carbon cycle perturbations3,7,15 (Fig. 2b). Furthermore, the relative invariability of both the Osi records and the 192Os abundances—which scale to organic matter content—implies that the balance of Os fluxes to the oceans and uptake of Os in sedimentary organic matter did not appreciably change during the MECO.
Although the magnitude of the negative Osi shift at Site 959 is small (~0.05), it exceeds the maximum analytical uncertainty (2σ = 0.01) by a factor of 5. The shift starts at the onset of MECO warming and is also present at Sites 1263 and U1333, where it is similar in magnitude (Figs. 1, 2). Interestingly, the Osi profile of Site U1333 is characterized by two separate excursions to lower, less radiogenic values rather than the gradual and continuous decrease that is observed at Site 959. The Osi profile at Site 1263 shows trends intermediate to Sites 959 and U1333. Nevertheless, the lowest Osi values in all three records occur toward the end of the MECO, which is coincident with the peak warming phase10. In addition, the return towards pre-MECO values is synchronous with the termination of the MECO at all three sites, implying that the Osi shift lasted for the entire duration of the event (~500 kyr). The absolute Osi values differ slightly between sites, likely because of differences in coastal proximity and oceanographic setting30,31. However, the general timing and magnitude of the Osi shift are reproduced at all sites, indicating that the Osi shift records a change in the 187Os/188Os composition of the global ocean. The global character and synchroneity of the Osi shift at the end of the MECO also indicate that osmium isotope stratigraphy is a promising tool for correlation of the event between sites in future studies (Fig. 2a).
In principle, the modest negative Osi shift during the MECO may be caused by an increase in the unradiogenic Os flux from hydrothermal and/or extraterrestrial sources, a decrease in the radiogenic Os flux from continental weathering, or a decrease in the 187Os/188Os composition of the continental weathering flux through a transient change in the exposure of different rock types, such as basalts7. There is no evidence for an extraterrestrial impact during the MECO. Furthermore, a reduction in continental silicate weathering rates during an episode of greenhouse warming seems paradoxical and unlikely, even though our Osi records clearly show no evidence of the expected increase in continental weathering. It is difficult to exclude a warming-induced change in regional climates and precipitation patterns—which could have affected the contributions of rock types with different 187Os/188Os compositions to the continental weathering flux3,32—as a cause for the Osi shift. However, this would still require a different cause for MECO warming.
Finally, the Osi shift could reflect a short-lived increase in mid-ocean ridge hydrothermal activity or an episode of increased volcanism and associated weathering of mafic silicate rocks24,33,34. Mass balance calculations with a progressive two-component mixing model that involves seawater and basalts (see Methods; Supplementary Data 2) show that the Osi shift across the MECO may correspond to a 10–15% increase in the contribution of the mantle-derived Os flux relative to the continental Os flux. Although there is no indication for the emplacement of a large igneous province during the middle Eocene8, an episode of volcanic activity at mid-ocean ridges or on land could have increased the Os flux from basalts, and consequently resulted in a decrease of the 187Os/188Os composition of the oceans that is consistent with our Osi records. Moreover, enhanced volcanism would provide a mechanism for the atmospheric CO2 rise that has been inferred for the MECO8,11, perhaps similar to the Late Cretaceous episode of greenhouse warming associated with volcanic eruptions from the Deccan Traps33,35,36. Potential events that have been dated at approximately the right age in the middle Eocene include (1) a pulse of metamorphic decarbonation associated with Himalayan uplift and metamorphism37,38, (2) increased arc volcanism around the Pacific rim39 and especially in the Caribbean, related to an ignimbrite flare-up in the Sierra Madre Occidental of Mexico40,41,42, (3) an episode of magmatism in the East African Rift zone43, in particular in Southern Ethiopia and Northern Kenya44,45, and/or (4) mid-ocean ridge volcanism in the North Atlantic, due to rifting in East Greenland and activity of the Iceland hotspot46,47,48. However, the timing and magnitude of these events are at present not sufficiently well resolved to establish a direct causal link with the MECO. Additionally, it is unclear if increased Himalayan uplift and metamorphism would be compatible with the observed negative Osi shift, as the Himalayas are generally considered to contribute relatively radiogenic Os to the continental weathering flux49,50. Yet, the effects of Himalayan uplift and subsequent weathering on the Cenozoic Osi record are likely small51,52.
Carbon and osmium cycle modeling
Enhanced volcanism and/or hydrothermal activity may represent the most parsimonious scenario to explain the modest negative Osi shift and atmospheric CO2 rise during the MECO. However, a strong silicate weathering response to greenhouse warming through focused weathering of fresh basalts is in disagreement with the extensive carbonate dissolution observed in deep ocean basins8,10. Therefore, total continental weathering fluxes must have remained approximately constant during the event. Collectively, the available data indicate that CO2 was added to the ocean-atmosphere system through enhanced volcanism, leading to warming, but was not neutralized through the silicate weathering feedback, leading to sustained ocean acidification.
To test the plausibility of scenarios involving enhanced volcanism and/or diminished continental weathering during the MECO, we performed a series of carbon cycle simulations with the box model LOSCAR23 by prescribing fluxes with the transient shift that is inferred from our Osi records (see Methods; Fig. 3; Supplementary Figs. 4–9). For consistency, we have also modeled the 187Os/188Os composition of the global ocean by applying the same LOSCAR carbon cycle fluxes as forcing to a box model of the Os cycle (see Methods; Supplementary Software 1). In addition to a ~0.05 decrease in the 187Os/188Os ratio of seawater, our target scenario for the MECO involves a rise in atmospheric CO2 concentrations, a slight increase in the δ13C of dissolved inorganic carbon in the deep ocean, and a shoaling of the CCD over ~500 kyr8. Since there are no high-resolution pCO2 records available for the MECO, the target scenario includes an approximate doubling of atmospheric CO2 concentrations relative to middle Eocene background values of 500–1000 ppmv11,53. Furthermore, the magnitude of CCD change during the event possibly varied between the different ocean basins10, so we incorporate a conservative estimate of at least 500 m shoaling in the Atlantic and Pacific in our target scenario.
LOSCAR and Os cycle model simulations of the most likely MECO scenario. a Forcing for two scenarios involving a gradual, transient 20% increase in the volcanic CO2 flux over ~500 kyr. The solid lines represent a scenario in which the silicate and carbonate weathering fluxes are allowed to vary in response to CO2 forcing (normal weathering feedback), while the dashed lines represent a scenario in which these weathering fluxes are kept constant (diminished weathering feedback). Only the latter scenario corresponds to all observations. b Model response in the 187Os/188Os composition of the global ocean, shown against smoothed fits to the MECO Osi records from the study sites. c Model CCD response of different ocean basins, shown against carbonate content (wt%) records for different depths in the Atlantic, Indian and Pacific oceans as compiled by Sluijs et al8. d Model atmospheric pCO2 response and pH response for the surface Atlantic and Pacific oceans. e Model δ13C response for the DIC of the deep Atlantic and Pacific oceans. For a full description of the LOSCAR model, see Zeebe23
All model simulations result in a decrease in the 187Os/188Os ratio of seawater (Fig. 3; Supplementary Figs. 4–9). Although a gradual, linear increase in volcanism of 10–20% over ~500 kyr is sufficient to cause CO2 accumulation in the ocean-atmosphere system, and hence global warming on MECO timescales, this scenario results in a deepening of the CCD instead of the observed shoaling (Fig. 3; Supplementary Fig. 4)8,10. A similar behavior of the CCD is observed in previous LOSCAR simulations of the MECO8 and the Late Cretaceous warming episode36. Crucially, the model is only able to reproduce CO2 rise in conjunction with shoaling of the CCD on these timescales if we invoke enhanced volcanism together with a diminished weathering feedback by maintaining the silicate and carbonate weathering fluxes constant (Fig. 3; Supplementary Fig. 5). Although the magnitude of this modeled CCD response is smaller than the shoaling inferred from deep-sea carbonate records8,10, we note that the model may underestimate CCD changes because it does not account for the additional effects of biological carbonate compensation54. Such a reduction in net carbonate production resulting from ocean acidification might amplify the CCD response for a given CO2 rise. In all scenarios, the model reproduces the modest increase in deep-sea benthic foraminifer δ13C values during the event9,10 because of a resulting decrease in carbonate versus constant organic carbon burial. Finally, the relatively rapid termination of the MECO is forced by a recovery of silicate weathering. We note that this does not need to represent a pronounced tectonic event, as the absolute magnitude of the flux imbalances is relatively small, but regionally enhanced weathering in the aftermath of the MECO would be consistent with observations from the Tethys region55.
To reconcile our Osi records and model results with global warming and atmospheric CO2 rise on MECO timescales8,9,10,11, we hypothesize that a long-term reduction in the strength of the negative silicate weathering feedback occurred in the Eocene56,57, due to a progressive reduction in the weatherability of the continents—the sum of all factors affecting chemical weathering other than climate58,59. For millions of years prior to the MECO, the Earth was generally characterized by high atmospheric CO2 levels53 and very warm climates60 (Supplementary Fig. 10), as well as flat continental relief. Thick, cation-depleted soils developed and transport-limited weathering regimes prevailed61,62, and consequently the weatherability of Earth's surface may have gradually decreased over the course of the Eocene. Indeed, such a progressive reduction in weathering feedback strength during the Eocene has been inferred from inverse modeling simulations of weathering fluxes based on Cenozoic pCO2 and δ13C records57 (Fig. 2c). With the strength of the weathering feedback strongly diminished, a small increase in volcanism or hydrothermal activity would lead to the accumulation of large amounts of CO2 in the ocean-atmosphere system, resulting in prolonged warming and ocean acidification during the MECO.
Changes in weatherability have also been suggested to explain other episodes of apparent decoupling between silicate weathering and climate59, for example during major glaciations in the Paleozoic and Neoproterozoic63,64,65. Our interpretations of a limited weathering response during the MECO suggest that a variable silicate weathering feedback strength57 can indeed act as a driver for sustained global warming on geological timescales, with potential importance to other enigmatic phases of carbon cycle change in Earth's history. Moreover, a variable weathering feedback strength governed by the interplay between tectonics, climate and the weatherability of the continents fundamentally challenges the parameterization of the silicate weathering feedback in carbon cycle models, especially those used to model transient perturbations such as the OAEs and the PETM. We therefore argue that future studies of these events should focus on exploring changes in temperature, atmospheric CO2, and the CCD in conjunction with the strength of the weathering feedback.
The samples used in this study were derived from middle Eocene sedimentary units from three different sites: organic-rich sediments from ODP Site 959 in the equatorial Atlantic along the African continental margin, carbonate-rich pelagic sediments from ODP Site 1263 on the Walvis Ridge in the south Atlantic, and carbonate-rich pelagic sediments from IODP Site U1333 in the equatorial Pacific (Supplementary Fig. 1). The total organic carbon (TOC) contents of these middle Eocene sediments range between 0.1 and 2%, with the highest TOC abundances occurring at Site 95966. Rock samples of 20–40 g were selected across the middle Eocene interval between approximately 42 and 38 Ma, with the highest sampling resolution across the MECO.
Bulk samples were freeze-dried or oven-dried at 50 °C and subsequently powdered using a ceramic pestle and mortar, in order to homogenize the Re and Os within the samples. Contact with metal surfaces was avoided so as not to contaminate the sample set. All Re and Os isotope analyses were performed at the Laboratory for Source Rock and Sulfide Geochronology and Geochemistry, and the Arthur Holmes Laboratory at the Durham Geochemistry Centre, Durham University (UK). Samples were digested in a CrO3-H2SO4 solution (0.25 g/g CrO3 in 8 mL of 4 N H2SO4) following the well-established methods of Selby & Creaser67, which have been shown to significantly limit the contribution of detrital Re and Os to the hydrogenous fraction bound to organic matter.
Powdered samples of approximately 0.3–1 g were loaded into Carius tubes with a known amount of 185Re + 190Os tracer solution (spike) and dissolved in 8 mL of CrO3-H2SO4 solution. The Carius tubes were then sealed and heated in an oven at 220 °C for 48 h. Osmium was isolated from the CrO3-H2SO4 sample solution by using solvent extraction with chloroform (CHCl3), and then back extracted by hydrobromic acid (HBr). The Os was further purified through micro distillation. Rhenium was isolated by evaporating 1 mL of the CrO3–H2SO4 sample solution to dryness, followed by solvent extraction involving sodium hydroxide (NaOH) and acetone (C3H6O)68,69. The Re was further purified by anion chromatography.
Following purification, the Re and Os fractions were loaded onto Ni and Pt filaments, respectively, together with 0.5 μL BaNO3 and BaOH activator solutions, respectively67. Rhenium and osmium isotope ratios were determined by negative thermal ionization mass spectrometry, using Faraday cups for Re and a Secondary Electron Multiplier for Os in peak-hopping mode.
Re and Os isotope ratios were corrected for instrumental mass fractionation, as well as spike and blank contributions. Procedural blanks for Re and Os in this study were 12 ± 3 pg/g and 0.07 ± 0.05 fg/g, respectively, with an 187Os/188Os value of 0.25 ± 0.15 (n = 3). The 187Re/188Os and 187Os/188Os uncertainties (2σ) include full propagation of uncertainties in weighing, mass spectrometer measurements, spike calibrations, blank corrections, and reproducibility of standards.
The 187Os/188Osinitial ratios (Osi) were calculated by correcting for post-depositional 187Re decay over time with the following equation:
$${}^{187}\mathrm{Os}/{}^{188}\mathrm{Os}_{{\mathrm{initial}}}\left( {\mathrm{Os}_{{i}}} \right) = {}^{187}\mathrm{Os}/{}^{188}\mathrm{Os}_{{\mathrm{measured}}} - {}^{187}\mathrm{Re}/{}^{188}\mathrm{Os}_{{\mathrm{measured}}} \ast \left( {{\mathrm{e}}^{\lambda t} - 1} \right)$$
where λ is the 187Re decay constant (1.666 · 10−11 yr−1)70 and t is the age of the rock. Given the high Re abundances in the organic-rich sediments from Site 959, we have used best estimates for the depositional ages of each of these samples. An age of 40 Ma was used for all samples from Sites 1263 and U1333, because improved age estimates would result in variations in Osi values of 0.1% or less on average. All results are listed in Supplementary Data 1. The Re–Os isotopic system is expected to have remained closed for the sample set, given that the cores were all fresh, unweathered, and showed no evidence of post-depositional events (e.g., veining, etc.). Further, where the Re–Os data has sufficient spread in isochron plot space to yield statistically robust isochrons, a geologically reasonable Re–Os isochron age is obtained (e.g., Site 959; see below for details).
Evaluation of Re and Os data
Although the studied samples were collected for evaluating changes in Osi rather than establishing isochrons, the Re–Os data of the sediments from Site 959 show a positive correlation between 187Re/187Os and 187Os/188Os, which results in an isochron age that is in good agreement with the age of the MECO between 40.5 and 40.0 Ma (Supplementary Figs. 2, 3). In contrast, the 187Re/187Os and 187Os/188Os data for Sites 1263 and U1333 do not have sufficient spread in isochron plot space, and hence cannot yield statistically geologically meaningful age estimates.
Age models
We adopt the age model of Cramwinckel et al.13 for Site 959 (Supplementary Fig. 11). This is based on initial71 and recently improved13 calcareous nannofossil biostratigraphy. The model also uses the long-term 187Os/188Os minimum at 34.65 Ma recorded at this site26, and TEX86 data that mark the MECO warming13. Moreover, we use the highest TEX86 value during the MECO peak warming and the lowest TEX86 value at the onset of the MECO as reported by Cramwinckel et al.13 to tentatively correlate to minima and maxima in the δ18O records of Bohaty et al.10, which were assigned ages of 40.06 and 40.52 Ma, respectively. Better age models are available for the other two sites. For Site 1263, we use a published age model10 based on magnetostratigraphy and bulk carbonate δ18O and δ13C chemostratigraphy. For Site U1333, an astronomically calibrated magnetostratigraphic age model72 was used in combination with calcareous nannofossil events73. All ages were adjusted to the framework of the GTS 201274 and tie points for the age models are listed in Supplementary Tables 1, 2 and 3.
Calculating changes in Os fluxes across the MECO
The 187Os/188Os composition of seawater is controlled by the balance between input fluxes from continental, mantle-derived, and extraterrestrial sources. However, the flux of extraterrestrial Os is generally assumed to be negligible and constant75,76, so our Osi records can be used to directly infer changes in relative contributions of the continental and mantle-derived Os sources across the MECO. To evaluate an increase in the mantle-derived Os flux, we developed a progressive, two-component mixing model for the release of Os from mantle-derived basalts that incorporates both the Os abundance and 187Os/188Os composition of seawater and basalts. This model is an adaptation of the two-component mixing model for strontium (Sr) isotopes of Faure (1986, Equations (9.2) and (9.10))77, with modifications to consider the larger range of Os isotope variations in comparison to Sr isotope variations.
From the relative molar concentrations of natural Os isotopes, we know:
$$\frac{{\left[ {{\mathrm{Os}}} \right] - \left[ {\;{}^{187}{\mathrm{Os}}} \right]}}{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]}} = 7.4$$
where [Os] represents the molar concentration (in mol/kg) of total Os (i.e., 186Os + 187Os + 188Os + 189Os + 190Os + 192Os), and [187Os] and [188Os] represent the molar concentrations (in mol/kg) of 187Os and 188Os, respectively78.
Equation (2) can be rewritten as:
$$\left[ {\;^{187}{\mathrm{Os}}} \right] = \frac{R}{{7.4 + R}}\left[ {{\mathrm{Os}}} \right]$$
$$\left[ {\;{}^{188}{\mathrm{Os}}} \right] = \frac{1}{{7.4 + R}}\left[ {{\mathrm{Os}}} \right]$$
where R = [187Os]/[188Os].
Two-component mixing between seawater and basalts can then be expressed for both 187Os and 188Os as:
$$\left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{mix}}} = \frac{{\left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{sw}}} \ast M_{{\mathrm{sw}}} + \left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{bas}}} \ast M_{{\mathrm{bas}}}}}{{M_{{\mathrm{sw}}} + M_{{\mathrm{bas}}}}}$$
$$\left[ {\;{}^{188}{\mathrm{Os}}} \right]_{\mathrm{mix}} = \frac{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]_{{\mathrm{sw}}} \ast M_{{\mathrm{sw}}} + \left[ {\;{}^{188}{\mathrm{Os}}} \right]_{{\mathrm{bas}}} \ast M_{{\mathrm{bas}}}}}{{M_{\mathrm{sw}} + M_{{\mathrm{bas}}}}}$$
where M represents the mass of a component (in kg) and the subscripts sw, bas and mix represent seawater, basalts and the eventual mix between the two, respectively.
We now define:
$$\Delta M_{{\mathrm{bas}}} = \frac{{M_{{\mathrm{bas}}}}}{{M_{{\mathrm{sw}},{\mathrm{initial}}}}}$$
$$f = \frac{{\Delta M_{{\mathrm{bas}}}}}{{M_{{\mathrm{sw}},{\mathrm{progressive}}} + \Delta M_{{\mathrm{bas}}}}} = \frac{{\Delta M_{{\mathrm{bas}}}}}{{M_{{\mathrm{mix}}}}}$$
where ΔMbas is an infinitesimal representing the mass of basalts added during a mixing step relative to the mass of seawater initially present, and f represents the amount of basalts added during a mixing step relative to the total amount of seawater and basalts present during progressive mixing (Mmix).
Equations (3)–(8) can then be combined as follows:
$$\begin{array}{l}\left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{mix}}} = f \ast \left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{bas}}} + \left( {1 - f} \right) \ast \left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{sw}}}\\ \hskip 70pt= f \ast \frac{{R_{{\mathrm{bas}}}}}{{7.4 + R_{{\mathrm{bas}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{bas}}} + \left( {1 - f} \right) \ast \frac{{R_{{\mathrm{sw}}}}}{{7.4 + R_{{\mathrm{sw}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{sw}}}\end{array}$$
$$\begin{array}{l}\left[ {\;{}^{188}{\mathrm{Os}}} \right]_{{\mathrm{mix}}} = f \ast \left[ {\;{}^{188}{\mathrm{Os}}} \right]_{{\mathrm{bas}}} + \left( {1 - f} \right) \ast \left[ {\;{}^{188}{\mathrm{Os}}} \right]_{{\mathrm{sw}}}\\ \hskip 70pt= f \ast \frac{1}{{7.4 + R_{{\mathrm{bas}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{bas}}} + \left( {1 - f} \right) \ast \frac{1}{{7.4 + R_{{\mathrm{sw}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{sw}}}\end{array}$$
Finally, dividing equation (9) by equation (10) yields:
$$\begin{array}{l}R_{{\mathrm{mix}}} = \frac{{\left[ {\;{}^{187}{\mathrm{Os}}} \right]_{{\mathrm{mix}}}}}{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]_{{\mathrm{mix}}}}} = \frac{{f \ast \frac{{R_{{\mathrm{bas}}}}}{{7.4 + R_{{\mathrm{bas}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{bas}}} + \left( {1 - f} \right) \ast \frac{{R_{{\mathrm{sw}}}}}{{7.4 + R_{{\mathrm{sw}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{sw}}}}}{{f \ast \frac{1}{{7.4 + R_{{\mathrm{bas}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{bas}}} + \left( {1 - f} \right) \ast \frac{1}{{7.4 + R_{{\mathrm{sw}}}}}\left[ {{\mathrm{Os}}} \right]_{{\mathrm{sw}}}}}\end{array}$$
where R is the 187Os/188Os composition of the corresponding components (i.e., seawater, basalts, and the eventual mix between the two). Equations (7)–(11) can then be used to estimate the extent of mixing between seawater and basalts during the MECO by progressively calculating Rmix until our observed Osi shift is reproduced (see Supplementary Data 2). We assumed the pre-MECO 187Os/188Os ratio of seawater to be ~0.55 based on an average of pre-MECO Osi values recorded for the three sites and the Os concentration of seawater to be 10 ppq (~53 fmol/kg, similar to present-day values)14. Furthermore, we used an 187Os/188Os ratio of 0.13 for the mantle and mantle-derived basalts79,80, as well as an Os abundance of 1 ppt (~5.3 pmol/kg) for basalts80. Finally, we assumed that the maximum amount of basalt that can theoretically be added to seawater represents ~1% of the total mass of the ocean, as estimated for OAE231,81, and used increments of 0.01% for the value of ΔMbas.
Based on an Osi shift of 0.05 from the pre-MECO value of ~0.55 to a peak MECO value of ~0.50, we calculated a relative increase in the mantle-derived Os flux of ~13% across the event, which would equal the addition of Os from basalts with a mass of ~0.13% relative to the total mass of the ocean (Supplementary Data 2). Similar results are obtained if we estimate the relative increase in the 188Os flux, rather than the total Os flux. It is important to note that mantle-derived Os could also have been released to seawater through direct addition from magmatic degassing or hydrothermal inputs instead of basalt dissolution, but regardless of the mechanism, a ~13% increase in the mantle-derived Os flux during the MECO would be sufficient to reproduce our observed Osi shift and would correspond to the cumulative release of ~9.4 · 106 mol of mantle-derived Os. We also performed our calculations with the Osi values of each individual site: for Site 959, an Osi shift from 0.560 to 0.505 would yield a relative increase in the mantle-derived Os flux of ~14%; for Site 1263, an Osi shift from 0.530 to 0.485 would yield an increase of ~12%; for Site U1333, an Osi shift from 0.515 to 0.460 would yield an increase of ~16%. These differences are most likely to be attributed to the resolution of our records. To accommodate for this range of flux estimates, we adopted a best estimate of 10–15% for the increase in the mantle-derived Os flux during the MECO, but also explored the effects of an increase of up to 20% because we are unlikely to have sampled the lowest Osi values in any of our records due to the relatively low resolution of our dataset.
LOSCAR and Os cycle modeling
Carbon cycle simulations were performed using the Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir (LOSCAR) model23. In this box model, modified from Walker and Kasting82, carbon and several other biogeochemical tracers (e.g., alkalinity, phosphate, oxygen) are cycled through atmospheric and oceanic reservoirs. The model ocean is coupled to a sediment module and consists of surface-water, intermediate-water, and deep-water boxes of the four main Paleogene ocean basins (Atlantic, Indian, Pacific and Tethys). The model is designed to simulate the PETM at 56 Ma, but the minor changes in paleogeography compared to the middle Eocene at 40 Ma are not of relevance to the simple LOSCAR model. In these simulations, we use default parameter settings for the Paleogene setup. Equilibrium pCO2 is set at 750 ppm, consistent with pCO2 estimates based on planktic foraminifer boron isotope ratios (δ11B)53, and by default, silicate and carbonate weathering are implemented in the model as a feedback response to atmospheric CO2 concentrations. The CCD definition follows the default LOSCAR setup and is taken as the sediment depth level at which sedimentary CaCO3 contents fall below 10 wt%.
We explored the effects of changes in volcanism and/or continental weathering with the constraints from our Osi records to assess which scenario is able to reproduce a more realistic MECO target. We first simulated several scenarios with a gradual, linear increase in the volcanic CO2 flux (+10, +15, and +20%) over ~500 kyr, either while allowing the silicate and carbonate weathering fluxes to vary in response to CO2 forcing (Supplementary Fig. 4), or while maintaining these weathering fluxes at constant values (Supplementary Fig. 5). Subsequently, we performed several simulations invoking silicate weathering as a forcing rather than a feedback, by prescribing a gradual, linear decrease in the silicate weathering flux (−10, −15, and −20%) over ~500 kyr, while keeping the volcanic CO2 flux and the carbonate weathering flux at constant values (Supplementary Fig. 6). Finally, we tested the effect of an increase in volcanism (+5%) combined with a decrease in silicate weathering (−5%) (Supplementary Fig. 7); the effect of a combined decrease in silicate and carbonate weathering (both −10%) (Supplementary Fig. 8); and the effect of a decrease in silicate weathering (−10%) while maintaining a carbonate weathering feedback (Supplementary Fig. 9). For an overview of all model scenarios, see Supplementary Table 4.
In order to demonstrate that our LOSCAR model simulations are consistent with the Osi records, the scenarios outlined above were also applied to a separate box model of the Os cycle. This Os cycle model is inspired by the work of Richter & Turekian83 and many subsequent studies, including Peucker-Ehrenbrink & Ravizza14. We fully derive the equations used to model the Os cycle in the ocean below.
We first define N as the total molar inventory of Os (including all Os isotopes) in seawater, and 187N and 188N as the molar inventories of 187Os and 188Os in seawater, respectively. The 187Os/188Os composition of seawater (Rsw) is thus expressed as:
$$R_{{\mathrm{sw}}} = \frac{{\;{}^{187}N}}{{\;{}^{188}N}}$$
Subsequently, changes in Rsw over time can be written as:
$$\begin{array}{l}\frac{{{\mathrm{d}}R_{{\mathrm{sw}}}}}{{{\mathrm{d}}t}} = \frac{{\mathrm{d}}}{{{\mathrm{d}}t}}\left( {\frac{{\,^{187}N}}{{\,^{188}N}}} \right) = \frac{{\,^{188}N\frac{{{\mathrm{d}}^{187}N}}{{{\mathrm{d}}t}} \,- \,^{187}N\frac{{{\mathrm{d}}^{188}N}}{{{\mathrm{d}}t}}}}{{\left( {\,^{188}N} \right)^2}} = \frac{1}{{\,^{188}N}}\left[ {\frac{{{\mathrm{d}}^{187}N}}{{{\mathrm{d}}t}} - R_{{\mathrm{sw}}}\frac{{{\mathrm{d}}^{188}N}}{{{\mathrm{d}}t}}} \right]\end{array}$$
Multiplying equation (13) by 188N gives:
$$\;{}^{188}N\frac{{{\mathrm{d}}R_{{\mathrm{sw}}}}}{{{\mathrm{d}}t}} = \frac{{{\mathrm{d}}{}^{187}N}}{{{\mathrm{d}}t}} - R_{{\mathrm{sw}}}\frac{{{\mathrm{d}}{}^{188}N}}{{{\mathrm{d}}t}}$$
Changes in N, 187N and 188N over time can then be written as follows:
$$\frac{{{\mathrm{d}}N}}{{{\mathrm{d}}t}} = F_{{\mathrm{riv}}} + F_{{\mathrm{hyd}}} + F_{{\mathrm{ext}}} - F_{{\mathrm{sed}}}$$
$$\begin{array}{l}\frac{{{\mathrm{d}}{}^{187}N}}{{{\mathrm{d}}t}} = F_{{\mathrm{riv}}}\left( {\frac{{\left[ {\;{}^{187}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{riv}}} + \,F_{{\mathrm{hyd}}}\left( {\frac{{\left[ {\;{}^{187}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{hyd}}}+ \,F_{{\mathrm{ext}}}\left( {\frac{{\left[ {\;{}^{187}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{ext}}} - \,F_{{\mathrm{sed}}}\left( {\frac{{\left[ {\;{}^{187}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{sed}}}\end{array}$$
$$\begin{array}{l}\frac{{{\mathrm{d}}{}^{188}N}}{{{\mathrm{d}}t}} = F_{{\mathrm{riv}}}\left( {\frac{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{riv}}} + \,F_{{\mathrm{hyd}}}\left( {\frac{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{hyd}}}\,+ \,F_{{\mathrm{ext}}}\left( {\frac{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{ext}}} - \,F_{{\mathrm{sed}}}\left( {\frac{{\left[ {\;{}^{188}{\mathrm{Os}}} \right]}}{{\left[ {{\mathrm{Os}}} \right]}}} \right)_{{\mathrm{sed}}}\end{array}$$
where F represents the fluxes of Os (in mol/yr) from and to various reservoirs and the subscripts sw, riv, hyd, ext and sed represent seawater, riverine, hydrothermal, extraterrestrial and sediment reservoirs, respectively14,83.
Substituting equations (3) and (4) into equations (16) and (17), respectively, yields:
$$\begin{array}{l}\frac{{{\mathrm{d}}{}^{187}N}}{{{\mathrm{d}}t}} = F_{{\mathrm{riv}}}\frac{{R_{{\mathrm{riv}}}}}{{7.4 + R_{{\mathrm{riv}}}}} + F_{{\mathrm{hyd}}}\frac{{R_{{\mathrm{hyd}}}}}{{7.4 + R_{{\mathrm{hyd}}}}}+\, F_{{\mathrm{ext}}}\frac{{R_{{\mathrm{ext}}}}}{{7.4 + R_{{\mathrm{ext}}}}} - F_{{\mathrm{sed}}}\frac{{R_{{\mathrm{sed}}}}}{{7.4 + R_{{\mathrm{sed}}}}}\end{array}$$
$$\begin{array}{l}\frac{{{\mathrm{d}}{}^{188}N}}{{{\mathrm{d}}t}} = F_{{\mathrm{riv}}}\frac{1}{{7.4 + R_{{\mathrm{riv}}}}} + F_{{\mathrm{hyd}}}\frac{1}{{7.4 + R_{{\mathrm{hyd}}}}}+ F_{{\mathrm{ext}}}\frac{1}{{7.4 + R_{{\mathrm{ext}}}}} - F_{{\mathrm{sed}}}\frac{1}{{7.4 + R_{{\mathrm{sed}}}}}\end{array}$$
Finally, substituting equations (18) and (19) into equation (14) and combining with equation (4) results in:
$$\begin{array}{l}\frac{N}{{7.4 + R_{{\mathrm{sw}}}}}\frac{{{\mathrm{d}}R_{{\mathrm{sw}}}}}{{{\mathrm{d}}t}} = F_{{\mathrm{riv}}}\frac{{R_{{\mathrm{riv}}} - R_{{\mathrm{sw}}}}}{{7.4 + R_{{\mathrm{riv}}}}} + F_{{\mathrm{hyd}}}\frac{{R_{{\mathrm{hyd}}} - R_{{\mathrm{sw}}}}}{{7.4 + R_{{\mathrm{hyd}}}}}+ F_{{\mathrm{ext}}}\frac{{R_{{\mathrm{ext}}} - R_{{\mathrm{sw}}}}}{{7.4 + R_{{\mathrm{ext}}}}} - F_{{\mathrm{sed}}}\frac{{R_{{\mathrm{sed}}} - R_{{\mathrm{sw}}}}}{{7.4 + R_{{\mathrm{sed}}}}}\end{array}$$
which relates changes in Rsw over time to the fluxes of total Os (F), the 187Os/188Os compositions of these fluxes (R) and the amount of total Os in the ocean (N). Because there is no isotopic fractionation associated with Os burial (i.e., Rsed = Rsw), the net effect of the sedimentary Os flux (Fsed) in equation (20) is zero.
Together, equations (15) and (20) can be used to simulate any transient perturbation of the Os cycle. We first constructed a steady state model based on flux estimates and 187Os/188Os values for the present-day Os cycle with a 187Os/188Os ratio of seawater of 1.06 (see Supplementary Table 5). For the middle Eocene Os cycle, we assumed that the total Os inventory and the total input and output fluxes of Os are similar to present-day values, and recalculated the steady state riverine and hydrothermal Os fluxes for the pre-MECO 187Os/188Os ratio of seawater of 0.55 by assuming that the 187Os/188Os composition of these fluxes has remained unchanged. Subsequently, we used scaled silicate weathering and volcanic degassing fluxes from the LOSCAR model simulations to force our model of the Os cycle. The modeled changes in the 187Os/188Os ratio of seawater are included in the respective figures of all model scenarios (Fig. 3 of the main text and Supplementary Figs. 4–9). The full code used to perform the Os cycle model simulations is included as an R script in Supplementary Software 1.
The authors declare that all data supporting the results of this study are available in the Supplementary Information files associated with this manuscript.
Berner, R. A., Lasaga, A. C. & Garrels, R. M. The carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years. Am. J. Sci. 283, 641–683 (1983).
Article ADS CAS Google Scholar
Walker, J. C. G., Hays, P. B. & Kasting, J. F. A negative feedback mechanism for the long-term stabilization of Earth's surface temperature. J. Geophys. Res. 86, 9776–9782 (1981).
Ravizza, G. E., Norris, R. N., Blusztajn, J. & Aubry, M.-P. An osmium isotope excursion associated with the late Paleocene thermal maximum: evidence of intensified chemical weathering. Paleoceanography 16, 155–163 (2001).
Article ADS Google Scholar
Pogge von Strandmann, P. A. E., Jenkyns, H. C. & Woodfine, R. G. Lithium isotope evidence for enhanced weathering during Oceanic Anoxic Event 2. Nat. Geosci. 6, 668–672 (2013).
Cohen, A. S., Coe, A. L., Harding, S. M. & Schwark, L. Osmium isotope evidence for the regulation of atmospheric CO2 by continental weathering. Geology 32, 157–160 (2004).
Vance, D., Teagle, D. A. H. & Foster, G. L. Variable quaternary chemical weathering fluxes and imbalances in marine geochemical budgets. Nature 458, 493–496 (2009).
Article ADS CAS PubMed Google Scholar
Dickson, A. J. et al. Evidence for weathering and volcanism during the PETM from Arctic Ocean and Peri-Tethys osmium isotope records. Palaeogeogr. Palaeoclimatol. Palaeoecol. 438, 300–307 (2015).
Sluijs, A., Zeebe, R. E., Bijl, P. K. & Bohaty, S. M. A middle Eocene carbon cycle conundrum. Nat. Geosci. 6, 429–434 (2013).
Bohaty, S. M. & Zachos, J. C. Significant Southern Ocean warming event in the late middle Eocene. Geology 31, 1017 (2003).
Bohaty, S. M., & Zachos, J. C., & Florindo, F., & Delaney, M. L. Coupled greenhouse warming and deep-sea acidification in the middle Eocene. Paleoceanography 24, PA2207 (2009).
Bijl, P. K. et al. Transient Middle Eocene atmospheric CO2 and temperature variations. Science 330, 819–821 (2010).
Boscolo Galazzo, F. et al. The middle Eocene climatic optimum: a multi-proxy record of paleoceanographic changes in the South Atlantic (ODP Site 1263). Paleoceanography 29, 1143–1161 (2014).
Cramwinckel, M. J. et al. Synchronous tropical and polar temperature evolution in the Eocene. Nature doi:10.1038/s41586-018-0272-2 (2018).
Peucker-Ehrenbrink, B. & Ravizza, G. E. The marine osmium isotope record. Terra Nov. 12, 205–219 (2000).
Peucker-Ehrenbrink, B ., & Ravizza, G. E . in Chapter 8: osmium isotope stratigraphy (eds F.M., Gradstein, J.G., Ogg, M., Schmitz, & G., Ogg) In The Geologic Time Scale 2012 2-Volume Set. 145–166 Elsevier 2012: Boston.
Rooney, A. D. et al. Tracking millennial-scale Holocene glacial advance and retreat using osmium isotopes: insights from the Greenland ice sheet. Quat. Sci. Rev. 138, 49–61 (2016).
Ravizza, G. E. & Turekian, K. K. The osmium isotopic composition of organic-rich marine sediments. Earth. Planet. Sci. Lett. 110, 1–6 (1992).
Cohen, A. S., Coe, A. L., Bartlett, J. M. & Hawkesworth, C. J. Precise Re-Os ages of organic-rich mudrocks and the Os isotope composition of Jurassic seawater. Earth. Planet. Sci. Lett. 167, 159–173 (1999).
Ravizza, G. E. & Zachos, J. C. Records of Cenozoic ocean chemistry. Treatise Geochem. 6, 551–581 (2003).
Them, T. R. et al. Evidence for rapid weathering response to climatic warming during the Toarcian Oceanic Anoxic Event. Sci. Rep. 7, 5003 (2017).
Article PubMed Central ADS CAS PubMed Google Scholar
Ridgwell, A. & Zeebe, R. E. The role of the global carbonate cycle in the regulation and evolution of the Earth system. Earth. Planet. Sci. Lett. 234, 299–315 (2005).
Pälike, H. et al. A Cenozoic record of the equatorial Pacific carbonate compensation depth. Nature 488, 609–614 (2012).
Zeebe, R. E. LOSCAR: Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir model v2.0.4. Geosci. Model Dev. 5, 149–166 (2012).
Cohen, A. S. & Coe, A. L. New geochemical evidence for the onset of volcanism in the Central Atlantic magmatic province and environmental change at the Triassic-Jurassic boundary. Geology 30, 267–270 (2002).
Ravizza, G. E. Osmium-isotope geochemistry of Site 959: implications for Re-Os sedimentary geochronology and reconstruction of past variations in the Os-isotopic composition of seawater. Proc. Ocean Drill. Program, Sci. Results 159, 181–186 (1998).
Ravizza, G. E., & Paquay, F. S. Os isotope chemostratigraphy applied to organic-rich marine sediments from the Eocene-Oligocene transition on the West African margin (ODP Site 959). Paleoceanography 23, PA2204 (2008).
Klemm, V., Levasseur, S., Frank, M., Hein, J. R. & Halliday, A. N. Osmium isotope stratigraphy of a marine ferromanganese crust. Earth. Planet. Sci. Lett. 238, 42–48 (2005).
Burton, K. W. Global weathering variations inferred from marine radiogenic isotope records. J. Geochem. Explor. 88, 262–265 (2006).
Nielsen, S. G. et al. Thallium isotope evidence for a permanent increase in marine organic carbon export in the early Eocene. Earth. Planet. Sci. Lett. 278, 297–307 (2009).
Paquay, F. S. & Ravizza, G. E. Heterogeneous seawater 187Os/188Os during the late Pleistocene glaciations. Earth. Planet. Sci. Lett. 349–350, 126–138 (2012).
Du Vivier, A. D. C. et al. Marine 187Os/188Os isotope stratigraphy reveals the interaction of volcanism and ocean circulation during Oceanic Anoxic Event 2. Earth. Planet. Sci. Lett. 389, 23–33 (2014).
Elsworth, G., & Galbraith, E., & Halverson, G., & Yang, S. Enhanced weathering and CO2 drawdown caused by latest Eocene strengthening of the Atlantic meridional overturning circulation. Nat. Geosci. 10, 213–216 (2017).
Ravizza, G. E. & Peucker-Ehrenbrink, B. Chemostratigraphic evidence of Deccan volcanism from the marine osmium isotope record. Science 302, 1392–1395 (2003).
Turgeon, S. C. & Creaser, R. A. Cretaceous oceanic anoxic event 2 triggered by a massive magmatic episode. Nature 454, 323–326 (2008).
Robinson, N., Ravizza, G. E., Coccioni, R., Peucker-Ehrenbrink, B. & Norris, R. D. A high-resolution marine 187Os/188Os record for the late Maastrichtian: distinguishing the chemical fingerprints of Deccan volcanism and the KP impact event. Earth. Planet. Sci. Lett. 281, 159–168 (2009).
Henehan, M. J., Hull, P. M., Penman, D. E., Rae, J. W. B. & Schmidt, D. N. Biogeochemical significance of pelagic ecosystem function: an end-Cretaceous case study. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150510 (2016).
Kerrick, D. M. & Caldeira, K. Paleoatmospheric consequences of CO2 released during early Cenozoic regional metamorphism in the Tethyan orogen. Chem. Geol. 108, 201–230 (1993).
Kerrick, D. M. & Caldeira, K. Was the Himalayan orogen a climatically significant coupled source and sink for atmospheric CO2 during the Cenozoic? Earth. Planet. Sci. Lett. 173, 195–203 (1999).
Cambray, H. & Cadet, J.-P. Synchronisme de l'activité volcanique d'arc: mythe ou réalité? Comptes rendus l'Académie des Sci. Série 2. Sci. la Terre Des. planètes 322, 237–244 (1996).
McDowell, F. W. & Mauger, R. L. K-Ar and U-Pb zircon chronology of late Cretaceous and Tertiary magmatism in central Chihuahua State, Mexico. Geol. Soc. Am. Bull. 106, 118–132 (1994).
Sigurdsson, H. et al. History of circum-Caribbean explosive volcanism: 40Ar/39Ar dating of tephra layers. Proc. Ocean Drill. Program, Sci. Results 165, 299–314 (2000).
Aguirre-Díaz, G. J. & Labarthe-Hernández, G. Fissure ignimbrites: fissure-source origin for voluminous ignimbrites of the Sierra Madre Occidental and its relationship with Basin and Range faulting. Geology 31, 773–776 (2003).
Bailey, D. K. Episodic alkaline igneous activity across. Afr.: Implic. causes Cont. Break-Up. Geol. Soc. Lond., Spec. Publ. 68, 91–98 (1992).
ADS Google Scholar
George, R., Rogers, N. & Kelley, S. Earliest magmatism in Ethiopia: evidence for two mantle plumes in one flood basalt province. Geology 26, 923 (1998).
Rooney, T. O. The Cenozoic magmatism of East-Africa: Part I — flood basalts and pulsed magmatism. Lithos 286–287, 264–301 (2017).
Torsvik, T. H. & Cocks, L. R. M. Norway in space and time: a Centennial cavalcade. Nor. Geol. Tidsskr. 85, 73–86 (2005).
Mjelde, R. et al. Magmatic and tectonic evolution of the North Atlantic. J. Geol. Soc. Lond. 165, 31–42 (2008).
Torsvik, T. H. et al. Continental crust beneath southeast Iceland. Proc. Natl. Acad. Sci. 112, E1818–E1827 (2015).
Singh, S. K., Trivedi, J. R. & Krishnaswami, S. Re-Os isotope systematics in black shales from the Lesser Himalaya: their chronology and role in the 187Os/188Os evolution of seawater. Geochim. Cosmochim. Acta 63, 2381–2392 (1999).
Pierson-Wickmann, A. C., Reisberg, L. & France-Lanord, C. The Os isotopic composition of Himalayan river bedloads and bedrocks: importance of black shales. Earth. Planet. Sci. Lett. 176, 203–218 (2000).
Sharma, M., Wasserburg, G. J., Hofmann, A. W. & Chakrapani, G. J. Himalayan uplift and osmium isotopes in oceans and rivers. Geochim. Cosmochim. Acta 63, 4005–4012 (1999).
Paul, M. et al. Dissolved osmium in Bengal plain groundwater: implications for the marine Os budget. Geochim. Cosmochim. Acta 74, 3432–3448 (2010).
Anagnostou, E. et al. Changing atmospheric CO2 concentration was the primary driver of early Cenozoic climate. Nature 533, 380–384 (2016).
Luo, Y., Boudreau, B. P., Dickens, G. R., Sluijs, A. & Middelburg, J. J. An alternative model for CaCO3 over-shooting during the PETM: biological carbonate compensation. Earth. Planet. Sci. Lett. 453, 223–233 (2016).
Spofforth, D. J. A. et al. Organic carbon burial following the middle Eocene climatic optimum in the central western Tethys. Paleoceanography 25, (2010).
Maher, K. & Chamberlain, C. P. Hydrologic regulation of chemical weathering and the geologic carbon cycle. Science 343, 1502–1504 (2014).
Caves, J. K., Jost, A. B., Lau, K. V. & Maher, K. Cenozoic carbon cycle imbalances and a variable weathering feedback. Earth. Planet. Sci. Lett. 450, 152–163 (2016).
Kump, L. R ., & Arthur, M. A . in Global chemical erosion during the Cenozoic: weatherability balances the budget In: W, Ruddiman (ed.) in Tectonic Uplift and Climate Change. 399–426 Plenum Press: New York.
Kump, L. R., Brantley, S. L. & Arthur, M. A. Chemical weathering, atmospheric CO2, and climate. Annu. Rev. Earth. Planet. Sci. 28, 611–667 (2000).
Zachos, J. C., Dickens, G. R. & Zeebe, R. E. An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature 451, 279–283 (2008).
West, A. J., Galy, A. & Bickle, M. J. Tectonic and climatic controls on silicate weathering. Earth. Planet. Sci. Lett. 235, 211–228 (2005).
Froelich, P. N. & Misra, S. Was the late Paleocene-early Eocene hot because earth was flat? An ocean lithium isotope view of mountain building, continental weathering, carbon dioxide, and Earth's Cenozoic climate. Oceanography 27, 36–49 (2014).
Kump, L. R. et al. A weathering hypothesis for glaciation at high atmospheric pCO2 during the late Ordovician. Palaeogeogr. Palaeoclimatol. Palaeoecol. 152, 173–187 (1999).
Mills, B., Watson, A. J., Goldblatt, C., Boyle, R. & Lenton, T. M. Timing of Neoproterozoic glaciations linked to transport-limited global weathering. Nat. Geosci. 4, 861–864 (2011).
Goddéris, Y. et al. Onset and ending of the late Palaeozoic ice age triggered by tectonically paced rock weathering. Nat. Geosci. 10, 382–386 (2017).
Wagner, T. Late Cretaceous to early Quaternary organic sedimentation in the eastern Equatorial Atlantic. Palaeogeogr. Palaeoclimatol. Palaeoecol. 179, 113–147 (2002).
Selby, D. & Creaser, R. A. Re–Os geochronology of organic rich sediments: an evaluation of organic matter analysis methods. Chem. Geol. 200, 225–240 (2003).
Li, C., Qu, W. J., Du, A. D. & Sun, W. Comprehensive study on extraction of rhenium with acetone in Re-Os isotopic dating. Rock. Miner. Anal. 28, 233–238 (2009).
Cumming, V. M., Selby, D. & Lillis, P. G. Re-Os geochronology of the lacustrine Green River formation: Insights into direct depositional dating of lacustrine successions, Re-Os systematics and paleocontinental weathering. Earth. Planet. Sci. Lett. 359–360, 194–205 (2012).
Smoliar, M. I., Walker, R. J. & Morgan, J. W. Re-Os ages of group IIA, IIIA, IVA, and IVB iron meteorites. Science 271, 1099–1102 (1996).
Shafik, S., Watkins, D. K. & Shin, I. C. Calcareous nannofossil Paleogene biostratigraphy, Côte d'Ivoire-Ghana Marginal Ridge, eastern Equatorial Atlantic. Proc. Ocean Drill. Program, Sci. Results 159, 413–431 (1998).
Pälike, H. et al. The heartbeat of the Oligocene climate system. Science 314, 1894–1898 (2006).
Toffanin, F., Agnini, C., Rio, D., Acton, G. & Westerhold, T. Middle Eocene to early Oligocene calcareous nannofossil biostratigraphy at IODP Site U1333 (equatorial Pacific). Micropaleontology 59, 69–82 (2013).
Gradstein, F. M., & Ogg, J. G., & Schmitz, M., & Ogg, G. The Geologic Time Scale 2012 2-Volume Set. (Elsevier: Amsterdam).
Pegram, W. J. & Turekian, K. K. The osmium isotopic composition change of Cenozoic sea water as inferred from a deep-sea core corrected for meteoritic contributions. Geochim. Cosmochim. Acta 63, 4053–4058 (1999).
Peucker-Ehrenbrink, B. Accretion of extraterrestrial matter during the last 80 million years and its effect on the marine osmium isotope record. Geochim. Cosmochim. Acta 60, 3187–3196 (1996).
Faure, G. Principles of Isotope Geology. (Wiley: New York, 1986)
Li, G. & Elderfield, H. Evolution of carbon cycle over the past 100 million years. Geochim. Cosmochim. Acta 103, 11–25 (2013).
Meisel, T., Walker, R. J., Irving, A. J. & Lorand, J. P. Osmium isotopic compositions of mantle xenoliths: a global perspective. Geochim. Cosmochim. Acta 65, 1311–1323 (2001).
Martin, C. E. Osmium isotopic characteristics of mantle-derived rocks. Geochim. Cosmochim. Acta 55, 1421–1434 (1991).
Hauff, F., Hoernle, K., Tilton, G., Graham, D. W. & Kerr, A. C. Large volume recycling of oceanic lithosphere over short time scales: geochemical constraints from the Caribbean Large Igneous Province. Earth. Planet. Sci. Lett. 174, 247–263 (2000).
Walker, J. C. G. & Kasting, J. F. Effects of fuel and forest conservation on future levels of atmospheric carbon dioxide. Palaeogeogr. Palaeoclimatol. Palaeoecol. 97, 151–189 (1992).
Richter, F. M. & Turekian, K. K. Simple models for the geochemical response of the ocean to climatic and tectonic forcing. Earth. Planet. Sci. Lett. 119, 121–131 (1993).
Westerhold, T. et al. Orbitally tuned timescale and astronomical forcing in the middle Eocene to early Oligocene. Clim. Past. 10, 955–973 (2014).
This research used samples provided by the Ocean Drilling Program (ODP) and Integrated Ocean Drilling Program (IODP). This work was carried out under the program of the Netherlands Earth System Science Centre (NESSC), which is financially supported by the Ministry of Education, Culture and Science (OCW) of the Netherlands. D.S. acknowledges support of the Total endowment fund and the Dida Scholarship of CUG Wuhan. Y.L. thanks Alan Rooney (Yale University) for intellectual freedom and a grant (SKL-K201706) from the Institute of Geology and Geophysics, Chinese Academy of Sciences. We thank Richard Zeebe (University of Hawai'i) for discussions on LOSCAR modeling, Antonia Hofmann, Geoff Nowell and Chris Ottley (Durham University) and Natasja Welters and Arnold van Dijk (Utrecht Geolab) for laboratory assistance.
Department of Earth Sciences, Faculty of Geosciences, Utrecht University, Utrecht, 3584 CB, The Netherlands
Robin van der Ploeg, Margot J. Cramwinckel, Jack J. Middelburg & Appy Sluijs
Department of Earth Sciences, Durham University, Durham, DH1 3LE, UK
David Selby
State Key Laboratory of Geological Processes and Mineral Resources, School of Earth Resources, China University of Geosciences, Wuhan, 430074, Hubei, China
Department of Geology and Geophysics, Yale University, New Haven, Connecticut, 06511, USA
Yang Li
State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, 10029, China
Ocean and Earth Science, National Oceanography Centre, University of Southampton Waterfront Campus, Southampton, SO14 3ZH, UK
Steven M. Bohaty
Robin van der Ploeg
Margot J. Cramwinckel
Jack J. Middelburg
Appy Sluijs
R.v.d.P., S.M.B., J.J.M., and A.S. designed the study. R.v.d.P. and D.S. generated the osmium isotope records, M.J.C. performed the carbon cycle simulations, R.v.d.P. and Y.L. conducted the osmium cycle modeling, and all authors contributed to data interpretation. R.v.d.P. wrote the manuscript with input from all authors.
Correspondence to Robin van der Ploeg.
The authors declare no competing interests.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Description of Additional Supplementary Files
Supplementary Data 1
Supplementary Software 1
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
van der Ploeg, R., Selby, D., Cramwinckel, M.J. et al. Middle Eocene greenhouse warming facilitated by diminished weathering feedback. Nat Commun 9, 2877 (2018). https://doi.org/10.1038/s41467-018-05104-9
Proxy evidence for state-dependence of climate sensitivity in the Eocene greenhouse
E. Anagnostou
E. H. John
G. L. Foster
Nature Communications (2020)
Fish proliferation and rare-earth deposition by topographically induced upwelling at the late Eocene cooling event
Junichiro Ohta
Kazutaka Yasukawa
Yasuhiro Kato
Scientific Reports (2020)
Carbon cycle instability and orbital forcing during the Middle Eocene Climatic Optimum
Martino Giorgioni
Luigi Jovane
Ercan Özcan
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Editors' Highlights
Show results from All journals This journal
Nature Communications (Nat Commun) ISSN 2041-1723 (online)