Evidence for a rapid release of carbon at the Paleocene-Eocene thermal maximum

Edited* by Wallace S. Broecker, Columbia University, Palisades, NY, and approved August 5, 2013 (received for review May 14, 2013)
September 16, 2013
110 (40) 15908-15913
Letter
Unsettled puzzle of the Marlboro clays
Peter Stassen, Robert P. Speijer, Ellen Thomas
Letter
Layering in the Paleocene/Eocene boundary of the Millville core is drilling disturbance
Paul N. Pearson, Christopher J. Nicholas
Letter
Onset of carbon isotope excursion at the Paleocene-Eocene thermal maximum took millennia, not 13 years
Richard E. Zeebe, Gerald R. Dickens [...] Ellen Thomas

Significance

Calcium carbonate and carbon isotope records from the rhythmically bedded Marlboro Clay, deposited during the onset of the PETM CIE, show that the massive release of isotopically light carbon was instantaneous, providing important constraints for the magnitude of carbon released and potential mechanisms.

Abstract

The Paleocene/Eocene thermal maximum (PETM) and associated carbon isotope excursion (CIE) are often touted as the best geologic analog for the current anthropogenic rise in pCO2. However, a causal mechanism for the PETM CIE remains unidentified because of large uncertainties in the duration of the CIE’s onset. Here, we report on a sequence of rhythmic sedimentary couplets comprising the Paleocene/Eocene Marlboro Clay (Salisbury Embayment). These couplets have corresponding δ18O cycles that imply a climatic origin. Seasonal insolation is the only regular climate cycle that can plausibly account for δ18O amplitudes and layer counts. High-resolution stable isotope records show 3.5‰ δ13C decrease over 13 couplets defining the CIE onset, which requires a large, instantaneous release of 13C-depleted carbon. During the CIE, a clear δ13C gradient developed on the shelf with the largest excursions in shallowest waters, indicating atmospheric δ13C decreased by ∼20‰. Our observations and revised release rate are consistent with an atmospheric perturbation of 3,000-gigatons of carbon (GtC).
Deep sea carbon isotope and CaCO3 records across the Paleocene/Eocene thermal maximum (PETM) and associated carbon isotope excursion (CIE) (55.8 Mya) require a massive addition of 13C-depleted carbon to the ocean–atmosphere system in a geologically short interval of time (13). Proposed mechanisms include the destabilization of the global methane reservoir by a thermal trigger (2, 4) or physical disturbance (5), production of thermogenic CH4 and CO2 during the emplacement of a large igneous province (6, 7), wildfires burning peatlands (8), desiccation of a large epicontinental sea (9), decomposition of terrestrial permafrost (10), and bolide impact (11, 12). The only consensus is that a precise chronology providing rates for the onset of the CIE is essential to distinguish among these mechanisms. Efforts to establish such a chronology have relied on average sedimentation rates based on integrated magneto- and biostratigraphic constraints (1), identifying cycles in deep sea cores and assigning an orbital periodicity (13, 14), or measuring the concentration of extraterrestrial 3He and applying an estimated flux to establish sedimentation rates (15, 16). The lack of a precise timescale for the δ13C excursion precludes further advances in identification of the source(s) of the light carbon and calculation of the release rates and magnitudes. Until these are better quantified, the relevance of the CIE as an analog for the Anthropocene remains speculative.
Deep sea bulk carbonate records show that the CIE began with an abrupt initial δ13C decrease of ∼1‰ followed by a more gradual decrease of similar magnitude (17). Bulk carbonate δ13C values from Ocean Drilling Program (ODP) Site 690 record the initial decrease in adjacent 1-cm samples. The cyclostratigraphic age model for Site 690 of Röhl et al. (13) constrains this initial decrease to <750 y, whereas the 3He method of Murphy et al. (16) places the duration of this initial decrease at closer to 30 kya. Herein lies the chronologic conundrum: the 750-y and 30 kya durations of the initial excursion are based on the same deep sea δ13C record and are essentially indistinguishable given the error in each methodology, yet have drastically different implications for the size of the carbon release necessary to produce the CIE. This uncertainty highlights the dire need for a precise chronometer for the onset of the carbon isotope excursion.
Here, we present high-resolution bulk stable isotope and %CaCO3 records from the northern Salisbury Embayment [35°N paleolatitude (18)] on the Atlantic coastal plain, using the Millville (ODP 174X) and newly recovered Wilson Lake B cores. Both of these cores contain the upper Paleocene-lower Eocene Marlboro Clay unit that has been correlated to deep sea CIE sections using carbon isotope and biostratigraphy (19, 20). These cores have two unique features that distinguish them as outstanding temporal archives for the onset of the CIE: (i) expanded sections of the Marlboro Clay in both the Wilson Lake B (15.5 m) and Millville (12.6 m) cores and (ii) distinct and rhythmic bedding of silty clays through the entire section containing the δ13C excursion.
The δ13C excursions in the Wilson Lake B and Millville cores have amplitudes of −6‰ and −4.5‰, respectively (Fig. S1), with ∼3.5‰ of the decrease representing the virtually instantaneous initial δ13C decrease observed in deep sea records (13, 17). In Wilson Lake B, this initial step in the δ13C decrease is recorded from 112.29 to 111.00 m, whereas the same decrease is found between 273.85 and 273.39 m in the Millville core. In contrast, the correlative δ13C decrease occurs over 1 cm at ODP Site 690 (13, 17). Therefore, our shallow marine records provide at least 45–130 times the temporal resolution afforded by the deep ocean record.
In the Wilson Lake B and Millville cores, the Marlboro Clay is characterized by rhythmic couplets of silty, kaolinitic clay (21, 22) distinguished by 1- to 2-mm layers of swelling smectite clays and micaceous silts, recurring every 1–3 cm through the entirety of the unit (mean, 1.9 ± 0.8 cm at 1σ) (Fig. 1). Harris et al. (19) noted similar layering in the Marlboro Clay interval in the Ancora core (ODP Leg 174X). The distinct layered beds are found in stratigraphically equivalent exposures in Medford, NJ, and in the South Dover Bridge core from Maryland (23). The prevalence and similarity of these couplets in at least five locations demonstrates that they are primary depositional features in the Marlboro Clay of the Salisbury Embayment. If these sedimentary cycles are demonstrably periodic, such layering offers the possibility of assigning a precise chronology to the onset of the CIE.
Fig. 1.
Examples of the remarkably cyclic sedimentary layering found throughout the Marlboro Clay in the Wilson Lake (A and B) and Millville (C) cores. (A) Core photograph of freshly split core from Wilson Lake B from 105.45 to 106.06 m (346–348 ft), ∼2-cm-scale swelling clay couplets. (B) Same interval after desiccation. Original swelling clays now appear as dark banding. (C) Whole round section of the Millville core (IODP 174AX) from 262.13 to 263.96 m (860–866 ft, white cards) showing the same cyclic clay couplets. Photo taken shortly after recovery, core was scraped but not split; however, the couplets are still visible.

Temporal Origin of the Couplets Within the Marlboro Clay

The remarkably rhythmic layers are the most prominent feature of the Marlboro Clay on first visual inspection of the freshly split Wilson Lake B cores (Fig. 1). Recognizing the potential for a climatic control, we selected two intervals from the Wilson Lake B core for high-resolution sampling (2 mm): the first from 111.313 to 111.542 m spans 15 layers, recording deposition just after the initial δ13C excursion (Fig. 2). A second transect between 107.936 and 108.036 m, spanning nine layers within the interval of low δ13C values that characterize the high-frequency signal at the height of the CIE (Fig. S2; SI Text). The high-resolution δ18O records from both segments covary with the physical expression of the layers, where each of the clay couplets corresponds to a single cycle in δ18O. Amplitudes of the δ18O cycles have a mean of 1.1 ± 0.5‰, including two cycles between 111.313 and 111.351 m that have amplitudes of 2‰.
Fig. 2.
Bulk δ18O record from Wilson Lake B core with the starting depth of 111.313 m (365.2 ft) in the core. Samples were taken every 2 mm. Gray bars show the position of the lamina relative to stable isotope samples. The sample transect was taken down the middle of the core. The core photo is below. The couplet tops are defined as the occurrence of swelling smectite clays and correspond to increasing or maxima in δ18O values.
The variability in δ18O reflects changes in temperature, the δ18O of water (δ18Owater), or a combination of both. If ascribed solely to temperature, the 1.1 ± 0.5‰ range in δ18O indicates oscillations of 4.8 ± 2.2 °C (24). If the δ18O of bulk carbonate (δ18Occ) was driven only by variations in δ18Owater resulting from changes in fresh water flux, the δ18O cycles indicate a range of 7.8 ± 3.6 practical salinity units (psu) in salinity [using a freshwater end-member δ18O of −5 ± 1‰ from the modern Carolinas (25) and mixing over a salinity range of 35 psu (Figs. S3 and S4; SI Text)]. We maintain that temperature must be a significant component of the intracouplet δ18O variability, especially for couplets that record 2‰ cycles; the required changes in salinity are far greater than are observed on the modern mid-Atlantic shelf (26, 27) or even at sites at comparable water depths off the Amazon fan (28). If half of the amplitude of the largest cycles were due to freshwater, the responsible process must be capable of producing cyclic temperature variations of ∼4.5 °C, affecting the entire shelf surface–water system.
The clay couplets in the Millville core are strikingly similar in character to those observed in Wilson Lake B and show δ18O cycles associated with each (Fig. 1). We count ∼750 couplets at Millville, beginning just before the initial δ13C decrease and continuing to the top of the unit, where the Marlboro is truncated by an erosional surface (total of 12.5 m). Presumably, this rhythmic bedding can be related to some cyclic depositional process, and given the tight association with periodic changes in δ18O, the most likely driver is something inherent to the climate system (29). Aside from tidal forcing of cyclic sedimentary packages, the dominant forcing within the climate system is related to changes in insolation. On depositional timescales, long period changes in Earth’s orbital parameters have the greatest influence on insolation: those of eccentricity (95–125 and 413 kya), obliquity (41 kya), and precession (19 and 23 kya).
For the sake of discussion, if we consider a scenario where the layers are related to the ∼20-kya precession of the equinoxes, where each clay couplet comprises a single cycle, the entirety of the unit would represent ∼15 My of deposition (750 layers at 20 kya/layer). This extreme duration is not supported by biostratigraphic constraints through the interval (23, 30), and more importantly, the sediments of the Marlboro Clay have reversed magnetizations (11) and were deposited entirely during magnetochron C24r, precluding an ∼15-My duration. Below the ∼20-kya band, but above the seasonal cycle, there is no cyclic change in radiative forcing due to changes in earth’s orbital parameters; a phenomenon referred to as the “orbital gap” by Munk et al. (31). A suite of periods and interfering harmonics within the millennial band have been suggested (32, 33), including a 1,500-y cycle observed in the North Atlantic that is purportedly orbitally driven (32). Even if each couplet was a single 1,500-y cycle, the Marlboro Clay sequence would represent about ∼1.13 My of deposition, which is not supported by any bio- or cyclo-stratigraphic evidence (13, 19, 34). More importantly, the expected radiative forcing due to these millennial-scale (or shorter, e.g., sunspot) cycles is only on the order of a few watts per meter squared (35, 36), and hence, these mechanisms lack the requisite insolation forcing to result in the observed δ18O changes reported here.
The next cycle of any substantial radiative consequence below orbital precession is due to seasonality (particularly at mid- and high latitudes). At 35°N, the inferred paleo-latitude for the Millville and Wilson Lake locations (18, 20), annual insolation varies from ∼200 to 500 Wm−2 (37) (far exceeding that of anything in the millennial band by several orders or even that of the precession cycle). A simple heat flux calculation indicates that a 300-Wm−2 variation would impart a 6 °C seasonal temperature cycle on the shelf (38). We use the Levitus and Boyer (39) gridded monthly data to estimate the contributions of temperature and δ18Owater on calcite precipitated on the modern Carolina shelf and predict a seasonal range in δ18Ocalcite of 1.2 ± 0.2‰ for shelf locations at 50-m water depth at ∼35°N (Fig. S4; SI Text). We note that oscillations in oxygen isotopes driven by seasonal temperature and/or fresh water flux variations in modern shelf settings are observed in shell transects of bivalves (4042) and planktonic foraminifera from sediment traps (43, 44). The δ18O cycles accompanying the layering observed in the Marlboro Clay are wholly consistent with seasonal changes forced by insolation at midlatitude sites.
The inference of seasonally paced deposition requires sedimentation rates on the order of 2 cm/y within the Salisbury Embayment. Mud accumulation rates in excess of 10 cm/y are observed on the modern Amazon shelf (45) and clays in the East China Sea (46), indicating that the high shelf sedimentation rate inferred from the Salisbury Embayment during the CIE cannot be excluded, especially under an enhanced hydrologic cycle (20). In fact, rapid accumulation of muds in ∼50-m water depth for decades to centuries are common on shelves in regions with high suspended sediment load (45, 47). Therefore, we contend that the rhythmic bedding throughout the Marlboro Clay represents annual deposition, where each couplet corresponds to a seasonal cycle, resulting from highly turbid waters in the Salisbury Embayment (20).

Timing of the Onset of the CIE from High-Resolution Stable Isotopes in Millville

The Millville core was chosen for high-resolution study through the CIE because the δ13C decrease (starting at 273.96 m) is wholly contained in the layered Marlboro Clay and therefore provides a chronometer through the CIE onset. We sampled the Millville core every 2 mm from the basal layer at 273.96–273.37 m, fully capturing the initial δ13C decrease. Between 273.772 and 273.532 m, the δ13C of bulk carbonate decreases by 3.9‰, from 1.45‰ to −2.48‰, over a span of 13 couplets (Fig. 3). This decrease is a maximum range because the −2.48‰ minimum value is part of a higher-frequency signal superimposed onto the near linear change of −3.5‰ over the 13 couplets. In contrast, %CaCO3 shows a more abrupt decrease, from 6% to 1% within one layer, with most of the change (4.25% to 1%) occurring across 4 mm (Fig. 3).
Fig. 3.
High-resolution bulk stable isotope and percent carbonate data from the Millville core (IODP Leg 174AX) through the onset of the CIE. Depth scales are in meters (Upper) and feet (Lower; the drillers used 10-ft drill rods and feet units are preserved here to facilitate comparison with the original core descriptions). The δ13C (red), δ18O (blue), and CaCO3 (green) records were generated by sampling the slabbed Millville cores at 2-mm intervals. Intervals with no data reflect too little or no CaCO3 (<0.1%) are shown by the gray bars. The position of the clay couplets is denoted by the inverted triangles at top.
Bulk δ18O decreases from −1.5‰ to −2.7‰ during the initial decrease in δ13C (Fig. 3). As with δ13C, the δ18O records higher-frequency cycles superimposed on this decreasing trend, which correspond to the depositional layers. We measured the thickness of each layer through the CIE recorded in Millville, using lightly polished slabs of the now dry core, and found a mean layer thickness through the interval of 19 ± 7 mm (SD). When the high-resolution δ18O data from the CIE interval are filtered using this layer frequency, the tight correspondence (Fig. S5) indicates that the oscillations in δ18O are reasonably described by the frequency of the sedimentary layers. Given that the oscillations in δ18O have large amplitudes (0.73 ± 0.23‰), implying regular ≥3 °C changes in temperature, the cyclicity in δ18O (and corresponding rhythmic bedding) are best explained by the seasonal insolation cycle. Therefore, we interpret the 13 observed layers through the onset of the CIE at Millville as 13 annual cycles, and the 750 layers within the Marlboro clay at Millville as representing 750 annual cycles.

Surface Water CO2(aq) Concentrations and Carbon Isotopic Equilibrium

The interpretation of annual couplets makes a testable prediction: the response time of the surface water carbonate system should be much more rapid than the total equilibration time between the surface ocean carbon reservoir and the atmosphere. This difference should be reflected in the %CaCO3 and bulk δ13C records from the shelf. Paleo-water depths are estimated (Fig. S6; SI Text) to be ∼60 m at Millville at the start of the PETM; an atmospheric pCO2 increase would propagate through the water column to this depth in a matter of weeks (48), increasing [CO2(aq)/H2CO3] enough to shift the carbonate equilibrium away from CO3=. In contrast, the δ13C response depends on the rate of CO2 exchange between the surface ocean and atmosphere (modern exchange is ∼90 gigatons of carbon per year) (49), and radiocarbon measurements have been used to quantify this isotopic equilibration as occurring on the order of a decade (50, 51).
Our ability to resolve the differential response times in %CaCO3 (<1 couplet) and δ13C (∼13 couplets) at Millville (Fig. 3) requires a very high sedimentation rate. An annual forcing for the couplets matches the observed changes in the modern ocean with respect to the rates of CO2 invasion (<1 y) and carbon isotopic equilibration (∼1 decade). Therefore, the most parsimonious explanation for the difference in response times between %CaCO3 and δ13C is a large and rapid injection of CO2, supporting our interpretation that the couplets reflect annual depositional cycles. Any other interpretation for the origin of the couplets must account for the temporal responses of the carbonate chemistry and isotope change observed on the shelf during the PETM CIE. Invoking longer period forcing creates a conundrum, requiring as yet unspecified sources of carbon and additional feedbacks to explain these observations of the PETM CIE (52).

Reconciling the Marlboro Clay and Deep Ocean Chronologies

The rhythmic bedding within the Marlboro Clay implies a timescale for the onset of the CIE that on first inspection appears much different (by two to three orders of magnitude) than the chronology based on deep sea records (decades vs. 1–10 kya). We agree that published chronologies for the recovery interval derived from deep sea records are correct (e.g., cycle counting vs. 3He). However, we argue that the anatomy of the CIE onset and its attendant chronology can only be obtained from recorders that are largely insulated from the open ocean, and the timescales being compared here are wholly compatible because (i) the Marlboro Clay records only the initiation of the CIE, <1000 y, and (ii) bulk δ13C records from the open ocean are fatally imprinted by mixing/diffusion with older carbon from the deep ocean and have much lower sedimentation rates.
The δ13C from the various Marlboro Clay sites indicate that most if not all Marlboro Clay sections do not record the recovery phase of the CIE documented in deep ocean sites. At the top of the Marlboro Clay, bulk δ13C values remain −2.5‰ to −4.5‰ lower than preexcursion values (Fig. 4 A and B). On the contrary, deep ocean sites have δ13C values that are <2.5‰ lower than preexcursion values at the start and <1‰ by the end of the recovery phase. The recovery phase of the CIE is well established as being on the order of 150 kya based on modeling (5355), and these time scales have been validated by temporally well-constrained pCO2 events associated with large igneous province emplacement (56, 57). A pCO2 pulse introduced into the ocean–atmosphere system with an attendant δ13C anomaly decays exponentially during the recovery phase, as modeled. If a CO2 pulse is added directly to the atmosphere, the atmosphere and surface ocean (to a lesser extent) will initially experience much higher CO2 concentrations and have a much larger δ13C anomaly than the deep ocean. However, this initial change is both extremely rapid and transient, making it very difficult to observe without a high-resolution recorder that is in close communication with the atmosphere. Lack of this critical observation to date should not be construed as a flaw in our hypothesis; rather, a specific prediction that warrants significant further work (see below). After the CO2 perturbation, the atmosphere and surface ocean will come into near equilibrium with the deep ocean on the 200- to 2,000-y scale (i.e., ∼70% of the initial pulse has dissipated) (58; Fig. S7) and then follow the well-constrained decay of the anomaly over ∼150 kya in the absence of additional perturbation, as predicted by geochemical modeling (5355).
Fig. 4.
(A) Comparison of the PETM carbon isotope excursion from several sites along a dip-transect from shallow to deep on the mid-Atlantic shelf (water depths calculated using simple projection onto a dip line with a slope of 1:1,000; SI Text). From deepest to shallowest: Bass River (red) (11), Ancora (green) (11), South Dover Bridge (black) (23), Clayton (dk. blue)(11), Wilson Lake B (lt. blue), Medford (orange; SI Text). Note the deepest sites show the smallest excursions, whereas the shallowest show the largest excursions, including the Medford locality, which is based on both outcrop and core. All Marlboro Clay sections are capped by an erosional unconformity. All Marlboro Clay sections are capped by an erosional unconformity. Depths in each core are recalculated with the initiation of the CIE defined as 0 m. (B) Relationship between total carbon isotope excursion SI Text and its paleo-water depth for each site along the dip transect. A logarithmic best fit (R2 = 0.91) indicates that the very shallowest sites, those in unbuffered communication with the atmosphere, should show a total excursion of approximately −20‰ (SI Text). Data from the South Dover Bridge core not included due to large uncertainty in projection onto the dip transect.
A CO2 perturbation with 13C-depleted carbon has different expressions in the atmosphere and surface ocean because of the mixing and diffusion with the deep ocean (SI Text). Measurements of surface ocean 14C show that the CO2 pool in those waters is a mix of modern atmosphere and much older CO2 ventilated from the deep ocean (59). Bomb 14C production during the 1950s and 1960s demonstrates that the timing and magnitude of a carbon isotope perturbation differs between the atmosphere, surface ocean, and deep ocean. The atmosphere recorded a near instantaneous shift with a large Δ14C response, whereas peak surface ocean invasion was about a decade later with an attenuated Δ 14C signal (51). This response is nearly identical to our observations of the onset of the CIE at Millville: the precipitous drop in wt% CaCO3 (evidence of a rapid increase in [CO2]aq), followed by a slower, decadal scale invasion of 13C-depleted carbon to the surface water system (Fig. 3). Most of the modern deep ocean has yet to be influenced by bomb 14C (Fig. S4). This modern experiment implies that the large (up to 8.2‰) δ13C signal of the Marlboro Clay represents the initial phase of ocean invasion (<2,000 y). It also predicts that the atmospheric response was larger than the 7‰ δ13C decrease observed in soil carbonates (60).
Aside from visual similarity, there is no other a priori reason to assign deep sea chronologies to shelf records (thereby forcing them to artificially conform to a recovery phase). In the case of the Salisbury Embayment, this miscorrelation is simply a function of our poor understanding on the exact timing of the erosional truncation at the top of the Marlboro Clay, within the overall global-scale context of the CIE. Therefore, our interpretation of the shelf sites containing the Marlboro Clay is that they do not preserve any of the later sediments that would comprise the recovery phase and hence do not record the time period of chemical and isotopic reequilibration of the surface ocean/atmosphere system with the deep ocean. This perspective effectively divorces the shelf chronologies from the deep ocean records as they pertain to the onset of the CIE. The δ13C values in the Marlboro Clay are best interpreted as a record of the initial CO2 perturbation to the atmosphere and its infiltration into the surface ocean. In fact, this interpretation is allowable by biostratigraphic data, which have 100- to 200-kya resolution at best (23). The apparent δ13C recovery in the upper part of the Marlboro actually represents the transfer of CO2 (with its anomalous δ13C value) from the atmosphere/surface ocean reservoir into the deep ocean. Using our timescale, the ∼750-y timespan recorded is wholly consistent with the 200- to 2,000-y timescale of Archer et al.’s (58) synthesis of model predictions for the current anthropogenic increase. Thus, the shelf recovery is a separate event from the recovery observed in the deep ocean carbonate records, in both its root cause and its timescale, and comparison between the two is inappropriate and has led to significant miscorrelation (6163).

Atmospheric Response to PETM CIE

A requirement of the near instantaneous release of 13C-depleted CO2 is a much larger excursion in the atmosphere than has been measured to date. We show δ13C excursions of approximately −8.2‰ at the shallowest shelf sites (Medford, 30-m water depth) to −3.5‰ at the deepest (Bass River, ∼73 m; Fig. 4A). At the peak of the CIE, a clear δ13C gradient had developed on the shelf, with the lowest values and largest excursions recorded at the shallowest sites (Medford) and the smallest excursions at the deepest sites (Bass River). When the total magnitude of each excursion is compared with paleo-water depth of the corresponding site, an unambiguous trend emerges that can be extrapolated to the origin, yielding an expected atmospheric perturbation of approximately −20‰ (Fig. 4B; SI Text). Because the magnitude of the δ13C anomalies are a function of the extent to which mixing with the deep ocean overprints each record (as evidenced by modern bomb 14C), we hypothesize that the total atmospheric excursion is much larger than the measured 8‰ excursion at Medford and has simply yet to be identified. Furthermore, as an atmospheric perturbation integrates over the successively larger and slower reacting reservoirs, the total magnitude of the excursion is attenuated; at the opposite extreme of the gradient, deep sea δ13C excursions (based on benthic foraminiferal records) are −2.7 ‰ in the South Atlantic and Southern Ocean (1) and approximately −2‰ in the Indo-Pacific (64).
We note that some shelf sections have been used to argue for a more protracted release of carbon (14, 65). Most notably, Cui et al. (14) identified a δ13C excursion of −4.5‰ in bulk organic matter from a section near Spitsbergen. The total excursion there is slightly less than predicted from our observations in the Salisbury Embayment (unless the site is at a paleo-water depth of ∼70 m; Fig. 4B), an effect that is probably related to the incorporation time and residence time of organic matter on the shelf, which is known to have a lag on the order of 1,000 y (66) and is differentially diachronous (67). More importantly, the slow release rate is a function of assigning precessional forcing to cyclicity in manganese and iron in constructing their timescale, in part based on the assumption of synchrony with the deep sea, which we argue above is inappropriate for the δ13C onset in shelf localities.

Implications for the Rate of Carbon Release and Sequence of Events at the PETM

The single greatest hurdle hindering understanding of the PETM CIE has been uncertainty in the timing of the carbon added to the ocean–atmosphere system: the release schedule greatly affects the amount of 13C-depleted carbon necessary to produce the globally observed CIE at any given isotope composition, due to the differential reaction time of Earth’s exchangeable carbon reservoirs (68, 69). The second unknown has been the size of the atmospheric response. Our high-resolution stable isotope records from the Marlboro Clay provide constraints for both. We demonstrate that the initial release was rapid, if not instantaneous. A best fit of the relationship between the total CIE magnitude and paleo-water depth at each site (Fig. 4B) predicts an atmospheric excursion of −20‰ (R2 = 0.91; SI Text). Assuming a pre-CIE atmospheric reservoir of 2,000 GtC (with a δ13C of −6‰) (70) and an instantaneous release, a mass balance calculation gives an estimate of the amount of carbon necessary to produce the ∼20‰ atmospheric excursion. No realistic amount of organic carbon (approximately −26‰) can produce a −20‰ atmospheric change (>100,000 GtC is needed). Thermogenic (−40‰) and biogenic methane (−60‰) sources would require 2,900 and 1,200 GtC, respectively, to produce the −20‰ atmospheric excursion. Given the rapidity of the onset, magnitude of the δ13C excursion, and that the observed calcite compensation depth shoaling in deep ocean requires ∼3,000 GtC (3), two mechanisms meet these criteria: large igneous province-produced thermogenic methane (6, 7) and cometary carbon (11, 12). The latter is consistent with the recent discovery of a substantial accumulation of nonbiogenic magnetic nanoparticles in the Marlboro clay, whose origin is best ascribed to impact condensate (71). If released as CO2, this would be consistent with observations of an ∼5 °C global warming, although we note that the radiative effect of a methane release, while short lived, is substantially greater than CO2. Finally, the revised timescale for the rate of carbon release at the onset of the PETM limits its usefulness as an analog for our current anthropogenic release.

Acknowledgments

This work benefitted greatly from discussions with Ken Miller, who noted that other cores from the New Jersey coastal plain displayed layering in the Marlboro, which prompted us to look at the Millville core in detail. Dennis Kent and Robert Kopp provided valuable discussions. Richard Mortlock, Nicole Abdul, and Elizabeth Miller helped with sample preparation and stable isotope analyses. We thank Jeff Severinghaus and Gerald Haug who provided insightful reviews that greatly benefitted the final manuscript. We also acknowledge three anonymous reviewers who provided helpful reviews on an earlier version of the manuscript. This work was funded by National Science Foundation Grant 0958867.

Supporting Information

Supporting Information (PDF)
Supporting Information
sd01.xlsx
sd02.xlsx
sd03.xlsx
sd04.xlsx
sd05.xlsx

References

1
JP Kennett, LD Stott, Abrupt deep-sea warming, palaeoceanographic changes and benthic extinctions at the end of the Palaeocene. Nature 353, 225–229 (1991).
2
GR Dickens, JR O'Neil, DK Rea, RM Owen, Dissociation of oceanic methane hydrate as a cause of the carbon isotope excursion at the end of the Paleocene. Paleoceanography 10, 965–971 (1995).
3
JC Zachos, et al., Rapid acidification of the ocean during the Paleocene-Eocene thermal maximum. Science 308, 1611–1615 (2005).
4
R Matsumoto, 13C anomalies of carbonates and a new paradigm 'gas-hydrate hypothesis'. J Geol Soc Jap 101, 902–924 (1995).
5
ME Katz, BS Cramer, GS Mountain, S Katz, KG Miller, Uncorking the bottle: What triggered the Paleocene/Eocene thermal maximum methane release? Paleoceanography 16, 549–562 (2001).
6
H Svensen, et al., Release of methane from a volcanic basin as a mechanism for initial Eocene global warming. Nature 429, 542–545 (2004).
7
M Storey, RA Duncan, CC Swisher, Paleocene-Eocene thermal maximum and the opening of the Northeast Atlantic. Science 316, 587–589 (2007).
8
AC Kurtz, LR Kump, MA Arthur, JC Zachos, A Paytan, Early Cenozoic decoupling of the global carbon and sulfur cycles. Paleoceanography 18, 1090 (2003).
9
JA Higgins, DP Schrag, Beyond methane: Towards a theory for the Paleocene-Eocene Thermal Maximum. Earth Planet Sci Lett 245, 523–537 (2006).
10
RM DeConto, et al., Past extreme warming events linked to massive carbon release from thawing permafrost. Nature 484, 87–91 (2012).
11
DV Kent, et al., A case for a comet impact trigger for the Paleocene/Eocene thermal maximum and carbon isotope excursion. Earth Planet Sci Lett 211, 13–26 (2003).
12
BS Cramer, DV Kent, Bolide summer: The Paleocene/Eocene thermal maximum as a response to an extraterrestrial trigger. Palaeogeogr Palaeoclimatol Palaeoecol 224, 144–166 (2005).
13
U Röhl, T Westerhold, TJ Bralower, JC Zachos, On the duration of the Paleocene–Eocene thermal maximum (PETM). Geochem Geophys Geosyst 8, 1–13 (2007).
14
Y Cui, et al., Slow release of fossil carbon during the Palaeocene-Eocene thermal maximum. Nat Geosci 4, 481–485 (2011).
15
KA Farley, SF Eltgroth, An alternative age model for the Paleocene-Eocene thermal maximum using extraterrestrial 3He. Earth Planet Sci Lett 208, 135–148 (2003).
16
BH Murphy, KA Farley, JC Zachos, An extraterrestrial 3He-based timescale for the Paleocene-Eocene thermal maximum (PETM) from Walvis Ridge, IODP Site 1266. Geochim Cosmochim Acta 74, 5098–5108 (2010).
17
S Bains, RM Corfield, RD Norris, Mechanisms of climate warming at the end of the paleocene. Science 285, 724–727 (1999).
18
RD Müller, M Sdrolias, C Gaina, B Steinberger, C Heine, Long-term sea-level fluctuations driven by ocean basin dynamics. Science 319, 1357–1362 (2008).
19
AD Harris, et al., Integrated stratigraphic studies of Paleocene–lowermost Eocene sequences, New Jersey Coastal Plain: Evidence for glacioeustatic control. Paleoceanography 25, PA3211 (2010).
20
RE Kopp, et al., An Appalachian Amazon? Magnetofossil evidence for the development of a tropical river-like system in the mid-Atlantic United States during the Paleocene-Eocene thermal maximum. Paleoceanography 24, PA4211 (2009).
21
NH Darton, The Marlboro clay. Econ Geol 43, 154–155 (1948).
22
T Gibson, L Bybell, D Mason, Stratigraphic and climatic implications of clay mineral changes around the Paleocene/Eocene boundary of the northeastern US margin. Sediment Geol 134, 65–92 (2000).
23
JM Self-Trail, DS Powars, DK Watkins, GA Wandless, Calcareous nannofossil assemblage changes across the Paleocene-Eocene thermal maximum: Evidence from a shelf setting. Mar Micropaleontol 92-93, 61–80 (2012).
24
S Epstein, R Buchsbaum, HA Lowenstam, HC Urey, Revised carbonate-water isotopic temperature scale. Geol Soc Am Bull 64, 1315–1326 (1953).
25
C Kendall, TB Coplen, Distribution of oxygen‐18 and deuterium in river waters across the United States. Hydrol Processes 15, 1363–1393 (2001).
26
Chang G, et al. (2002) Nearshore physical processes and bio-optical properties in the New York Bight. J. Geophys. Res 107(C9):3133 (10.10292001).
27
LP Atkinson, J Blanton, W Chandler, T Lee, Climatology of the southeastern United States continental shelf waters. J Geophys Res 88, 4705–4718 (1983).
28
SJ Lentz, R Limeburner, The Amazon River Plume during AMASSEDS: Spatial characteristics and salinity variability. J Geophys Res 100, 2355–2375 (1995).
29
FB Van Houten, Cyclic sedimentation and the origin of analcime-rich upper Triassic Lockatong Formation, west-central New Jersey and Adjacent Pennsylvania. Am J Sci 260, 561–576 (1962).
30
A Kahn, MP Aubry, Provincialism associated with the Paleocene/Eocene thermal maximum: Temporal constraint. Mar Micropaleontol 52, 117–131 (2004).
31
W Munk, M Dzieciuch, S Jayne, Millennial climate variability: Is there a tidal connection? J Clim 15, 370–385 (2002).
32
G Bond, et al., Persistent solar influence on North Atlantic climate during the Holocene. Science 294, 2130–2136 (2001).
33
CD Keeling, TP Whorf, The 1,800-year oceanic tidal cycle: A possible cause of rapid climate change. Proc Natl Acad Sci USA 97, 3814–3819 (2000).
34
Sugarman PJ, Miller KG, Browning JV, et al. (2005) Millville Site. Proceedings of the Ocean Drilling Program, Initial reports, eds Miller KG, Sugarman PJ, Browning JV, et al. (Ocean Drilling Program, College Station, TX), Volume 174AX (Suppl.), 1–94.
35
P Foukal, C Fröhlich, H Spruit, TM Wigley, Variations in solar luminosity and their effect on the Earth’s climate. Nature 443, 161–166 (2006).
36
GC Reid, Solar total irradiance variations and the global sea surface temperature record. J Geophys Res 96, 2835–2844 (1991).
37
P Huybers, Early Pleistocene glacial cycles and the integrated summer insolation forcing. Science 313, 508–511 (2006).
38
Visbeck M, Chassignet EP, Curry RG, Delworth TL, Dickson RR, Krahmann G (2013) The ocean's response to North Atlantic oscillation variability, in The North Atlantic oscillation: Climatic significance and environmental impact, eds Hurrell JW, Kushnir Y, Ottersen G, Visbeck M, (American Geophysical Union, Washington, DC).
39
Levitus S, Boyer TP (1994) World Ocean Atlas 1994. Volume 2. Oxygen. (National Environmental Satellite, Data, and Information Service, Washington, DC).
40
DE Krantz, AT Kronick, DF Williams, A model for interpreting continental-shelf hydrographic processes from the stable isotope and cadmium: Calcium profiles of scallop shells. Palaeogeogr Palaeoclimatol Palaeoecol 64, 123–140 (1988).
41
DS Jones, DF Williams, MA Arthur, Growth history and ecology of the Atlantic surf clam, Spisula solidissima (Dillwyn), as revealed by stable isotopes and annual shell increments. J Exp Mar Biol Ecol 73, 225–242 (1983).
42
BR Schöne, et al., Climate records from a bivalved Methuselah (Arctica islandica, Mollusca; Iceland). Palaeogeogr Palaeoclimatol Palaeoecol 228, 130–148 (2005).
43
DF Williams, AW Bé, RG Fairbanks, Seasonal stable isotopic variations in living planktonic foraminifera from Bermuda plankton tows. Palaeogeogr Palaeoclimatol Palaeoecol 33, 71–102 (1981).
44
R Thunell, E Tappa, C Pride, E Kincaid, Sea-surface temperature anomalies associated with the 1997–1998 El Nino recorded in the oxygen isotope composition of planktonic foraminifera. Geology 27, 843–846 (1999).
45
SA Kuehl, DJ DeMaster, CA Nittrouer, Nature of sediment accumulation on the Amazon continental shelf. Cont Shelf Res 6, 209–225 (1986).
46
DJ DeMaster, BA McKee, CA Nittrouer, Q Jiangchu, C Guodong, Rates of sediment accumulation and particle reworking based on radiochemical measurements from continental shelf deposits in the East China Sea. Cont Shelf Res 4, 143–158 (1985).
47
CA Nittrouer, et al., The geological record preserved by Amazon shelf sedimentation. Cont Shelf Res 16, 817–841 (1996).
48
WS Broecker, TH Peng Tracers in the Sea (Eldigio Press Lamont Doherty Geological Observatory, Palisades, NY), pp. 690 (1982).
49
U Siegenthaler, JL Sarmiento, Atmospheric carbon dioxide and the ocean. Nature 365, 119–125 (1993).
50
R Revelle, HE Suess, Carbon dioxide exchange between atmosphere and ocean and the question of an increase of atmospheric CO2 during the past decades. Tellus 9, 18–27 (1957).
51
T-H Peng, RM Key, HG Östlund, Temporal variations of bomb radiocarbon inventory in the Pacific Ocean. Mar Chem 60, 3–13 (1998).
52
Zeebe RE, Zachos JC (2013) Long-term legacy of massive carbon input to the Earth system: Anthropocene vs. Eocene. Phil Trans R Soc A 371(2001):20120006.
53
K Caldeira, MR Rampino, Carbon dioxide emissions from Deccan volcanism and a K/T boundary greenhouse effect. Geophys Res Lett 17, 1299–1302 (1990).
54
C Dessert, et al., Erosion of Deccan Traps determined by river geochemistry: Impact on the global climate and the Sr-87/Sr-86 ratio of seawater. Earth Planet Sci Lett 188, 459–474 (2001).
55
GR Dickens, MM Castillo, JCG Walker, A blast of gas in the latest Paleocene: Simulating first-order effects of massive dissociation of oceanic methane hydrate. Geology 25, 259–262 (1997).
56
MF Schaller, JD Wright, DV Kent, Atmospheric PCO₂ perturbations associated with the Central Atlantic Magmatic Province. Science 331, 1404–1409 (2011).
57
MF Schaller, JD Wright, DV Kent, PE Olsen, Rapid emplacement of the Central Atlantic Magmatic Province as a net sink for CO2. Earth Planet Sci Lett 323–324, 27–39 (2012).
58
D Archer, et al., Atmospheric lifetime of fossil fuel carbon dioxide. Annu Rev Earth Planet Sci 37, 117–134 (2009).
59
M Stuiver, GW Pearson, TF Braziunas, Radiocarbon age calibration of marine samples back to 9000 cal yr BP. Radiocarbon 28, 980–1021 (2006).
60
GJ Bowen, BB Bowen, Mechanisms of PETM global change constrained by a new record from central Utah. Geology 36, 379–382 (2008).
61
A Sluijs, et al., Environmental precursors to rapid light carbon injection at the Palaeocene/Eocene boundary. Nature 450, 1218–1221 (2007).
62
JC Zachos, et al., Extreme warming of mid-latitude coastal ocean during the Paleocene-Eocene thermal maximum: Inferences from TEX86 and isotope data. Geology 34, 737–740 (2006).
63
P Stassen, E Thomas, RP Speijer, Integrated stratigraphy of the Paleocene‐Eocene thermal maximum in the New Jersey Coastal Plain: Toward understanding the effects of global warming in a shelf environment. Paleoceanography 27, PA4210 (2012).
64
F Nunes, RD Norris, Abrupt reversal in ocean overturning during the Palaeocene/Eocene warm period. Nature 439, 60–63 (2006).
65
MJ Nicolo, GR Dickens, CJ Hollis, South Pacific intermediate water oxygen depletion at the onset of the Paleocene-Eocene thermal maximum as depicted in New Zealand margin sections. Paleoceanography 25, PA4210 (2010).
66
G Mollenhauer, et al., An evaluation of 14C age relationships between co-occurring foraminifera, alkenones, and total organic carbon in continental margin sediments. Paleoceanography 20, PA1016 (2005).
67
A Schneider-Mor, GJ Bowen, Coupled and decoupled responses of continental and marine organic-sedimentary systems through the Paleocene-Eocene thermal maximum, New Jersey margin. Paleoceanography 28, 105–115 (2013).
68
K Panchuk, A Ridgwell, L Kump, Sedimentary response to Paleocene-Eocene thermal maximum carbon release: A model-data comparison. Geology 36, 315–318 (2008).
69
RE Zeebe, JC Zachos, GR Dickens, Carbon dioxide forcing alone insufficient to explain Palaeocene-Eocene thermal maximum warming. Nat Geosci 2, 576–580 (2009).
70
B Hönisch, et al., The geological record of ocean acidification. Science 335, 1058–1063 (2012).
71
H Wang, DV Kent, MJ Jackson, Evidence for abundant isolated magnetic nanoparticles at the Paleocene-Eocene boundary. Proc Natl Acad Sci USA 110, 425–430 (2013).

Information & Authors

Information

Published in

Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 110 | No. 40
October 1, 2013
PubMed: 24043840

Classifications

Submission history

Published online: September 16, 2013
Published in issue: October 1, 2013

Keywords

  1. carbon cycle
  2. climate change

Acknowledgments

This work benefitted greatly from discussions with Ken Miller, who noted that other cores from the New Jersey coastal plain displayed layering in the Marlboro, which prompted us to look at the Millville core in detail. Dennis Kent and Robert Kopp provided valuable discussions. Richard Mortlock, Nicole Abdul, and Elizabeth Miller helped with sample preparation and stable isotope analyses. We thank Jeff Severinghaus and Gerald Haug who provided insightful reviews that greatly benefitted the final manuscript. We also acknowledge three anonymous reviewers who provided helpful reviews on an earlier version of the manuscript. This work was funded by National Science Foundation Grant 0958867.

Notes

*This Direct Submission article had a prearranged editor.

Authors

Affiliations

James D. Wright1 [email protected]
Department of Earth and Planetary Sciences, Rutgers University, Piscataway, NJ 08854; and
Morgan F. Schaller
Department of Earth and Planetary Sciences, Rutgers University, Piscataway, NJ 08854; and
Department of Geological Sciences, Brown University, Providence, RI 02912

Notes

1
To whom correspondence should be addressed. E-mail: [email protected].
Author contributions: J.D.W. and M.F.S. designed research, performed research, analyzed data, and wrote the paper.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Evidence for a rapid release of carbon at the Paleocene-Eocene thermal maximum
    Proceedings of the National Academy of Sciences
    • Vol. 110
    • No. 40
    • pp. 15851-16283

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media