Eocene sedimentary calcium carbonate contents and stable isotope composition of benthic foraminifera

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

Collection of data on particulate and dissolved soil phosphorus in the Jena Experiment (Main Experiment, time series since 2002)

This data set contains two time series of measurements of dissolved phosphorus (organic, inorganic and total with a biweekly resolution) and dissolved inorganic phosphorus with a seasonal resolution. In addition, data on phosphorus from soil samples measured in 2007 and fractionated by different acid-extrations (Hedley fractions) are provided. All data measured at the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. 1. Dissolved phosphorus in soil solution: Suction plates installed on the field site in 10, 20, 30 and 60 cm depth were used to sample soil pore water. Cumulatively extracted soil solution was collected every two weeks from October 2002 to May 2006. The biweekly samples from 2002, 2003 and 2004 were analyzed for dissolved organic phosphorus (DOP), dissolved inorganic phosphorus (PO4P) and dissolved total phosphorus (TDP) by Continuous Flow Analyzer (CFA SAN ++, SKALAR [Breda, The Netherlands]). 2. Seasonal values of dissolved inorganic phosphorus in soil solution were calculated as volume-weighted mean values of the biweekly measurements (spring = March to May, summer = June to August, fall = September to November, winter = December to February). 3. Phosphorus fractions in soil: Five independent soil samples per plot were taken in a depth of 0-15 cm using a soil corer with an inner diameter of 1 cm. The five samples per plot were combined to one composite sample per plot. A four-step sequential P fractionation (Hedley fractions) was applied and concentrations of P fractions in soil were measured photometrically (molybdenum blue-reactive P) with a Continuous Flow Analyzer (Bran&Luebbe, Germany).

Dissolved phosphorus in soil solution of the Jena Experiment (Main Experiment, year 2003)

This data set contains measurements of dissolved phosphorus (total dissolved nitrogen: TDP, dissolved inorganic phosphorus: PO4P and dissolved organic phosphorus: DOP) in samples of soil water collected in 2003 from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 mm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. Manual soil matric potential measurements were used to regulate the vacuum system. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled bi-weekly, in 2003 at the 07.03.2003; 24.03.2003; 07.04.2003; 22.04.2003; 07.05.2003; 20.05.2003; 03.06.2003; 28.07.2003; 12.09.2003; 22.09.2003; 07.10.2003; and 21.10.2003, and analyzed for dissolved inorganic P (PO4P) and total dissolved phosphorus (TDP). Inorganic phosphorus concentrations in the soil solution were measured photometrically with a continuous flow analyzer (CFA SAN++, Skalar [Breda, The Netherlands]). Ammonium molybdate catalyzed by antimony tartrate reacts in an acidic medium with phosphate and forms a phospho-molybdic acid complex. Ascorbic acid reduces this complex to an intensely blue-colored complex. Total dissolved P in soil solution was analyzed by irradiation with UV and oxidation with K2S2O8 followed by reaction with ammonium molybdate (Skalar catnr. 503-553w/r). As the molybdic complex forms under strongly acidic conditions, we could not exclude the hydrolysis of labile organic P compounds in our samples. Furthermore, the molybdate reaction is not sensitive for condensed phosphates. The detection limits of both TDP and PO4P were 0.02 mg P l-1 (CFA, Skalar). Dissolved organic P (DOP) in soil solution was calculated as the difference between TDP and PO4P. In a low number of samples, TDP was equal to or smaller than PO4P; in these cases, DOP was assumed to be zero.

Dissolved phosphorus in soil solution of the Jena Experiment (Main Experiment, year 2004)

This data set contains measurements of dissolved phosphorus (total dissolved nitrogen: TDP, dissolved inorganic phosphorus: PO4P and dissolved organic phosphorus: DOP) in samples of soil water collected in 2004 from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. Manual soil matric potential measurements were used to regulate the vacuum system. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled bi-weekly, in 2004 at the 15.01.2004; 30.01.2004; 12.02.2004; 27.02.2004; 09.03.2004; 25.03.2004; 21.04.2004; 07.05.2004; and 24.05.2004, and analyzed for dissolved inorganic P (PO4P) and total dissolved phosphorus (TDP). Inorganic phosphorus concentrations in the soil solution were measured photometrically with a continuous flow analyzer (for samples collected until spring 2004: CFA SAN++, Skalar [Breda, The Netherlands]; for samples collected later: CFA Autoanalyzer [Bran&Luebbe, Norderstedt, Germany]). Ammonium molybdate catalyzed by antimony tartrate reacts in an acidic medium with phosphate and forms a phospho-molybdic acid complex. Ascorbic acid reduces this complex to an intensely blue-colored complex. Total dissolved P in soil solution was analyzed by irradiation with UV and oxidation with K2S2O8 followed by reaction with ammonium molybdate (Skalar catnr. 503-553w/r). As the molybdic complex forms under strongly acidic conditions, we could not exclude the hydrolysis of labile organic P compounds in our samples. Furthermore, the molybdate reaction is not sensitive for condensed phosphates. The detection limits of both TDP and PO4P were 0.02 mg P l-1 (CFA, Skalar) and 0.04 mg P l-1 (Autoanalyzer, Bran&Luebbe). Dissolved organic P (DOP) in soil solution was calculated as the difference between TDP and PO4P. In a low number of samples, TDP was equal to or smaller than PO4P; in these cases, DOP was assumed to be zero.

Dissolved nitrogen in soil solution of the Jena Experiment (Main Experiment, year 2004)

This data set contains measurements of dissolved nitrogen (total dissolved nitrogen: TDN, dissolved organic nitrogen: DON, dissolved ammonium: NH4+, and dissolved nitrate: NO3-) in samples of soil water collected from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. In April 2002 glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for nitrate (NO3-) and ammonium (NH4+) concentrations with a continuous flow analyzer (CFA, Skalar, Breda, The Netherlands). Nitrate was analyzed photometrically after reduction to NO2- and reaction with sulfanilamide and naphthylethylenediamine-dihydrochloride to an azo-dye. Our NO3- concentrations contained an unknown contribution of NO2- that is expected to be small. Simultaneously to the NO3- analysis, NH4+ was determined photometrically as 5-aminosalicylate after a modified Berthelot reaction. The detection limits of NO3- and NH4+ were 0.02 and 0.03 mg N L-1, respectively. Total dissolved N in soil solution was analyzed by oxidation with K2S2O8 followed by reduction to NO2- as described above for NO3-. Dissolved organic N (DON) concentrations in soil solution were calculated as the difference between TDN and the sum of mineral N (NO3- + NH4+). In 5% of the samples, TDN was equal to or smaller than mineral N. In these cases, DON was assumed to be zero.

Analysis of organic matter from DSDP Leg 81 Holes

Organic geochemical and visual kerogen analyses were carried out on approximately 50 samples from Leg 81 (Rockall Plateau, North Atlantic). The sediments are from four sites (Sites 552-555), Pleistocene to Paleocene in age, and represent significantly different depositional environments and sources of organic matter. The Pleistocene glacial-interglacial cycles show differences in sedimentary organic matter based on Rock-Eval pyrolysis, organic phosphorus, and pyrolysis/mass-spectrometry analyses. Glacial samples contain more organic carbon, with a larger proportion of reworked organic matter. This probably reflects increased erosion of continental and shelf areas as a result of low sea level stands. Inter glacial samples contain a larger proportion of marine organic matter as determined by organic phosphorus and pyrolysis analyses. This immature, highly oxidized marine organic matter may be associated with the skeletal organic matrix of calcareous organisms. In addition, Rock-Eval data indicate no significant inorganic-carbonate contribution to the S3 pyrolysis peak. The Pliocene-Miocene sediments consist of pelagic, biogenic carbonates. The organic matter is similar to that of the Pleistocene interglacial periods; a mixture of oxidized marine organic matter and reworked, terrestrial detritus. The Paleocene-Oligocene organic matter reflects variations in source and depositional factors associated with the isolation of Rockall from Greenland. Paleocene sediments contain primarily terrestrial organic matter with evidence of in situ thermal stress resulting from interbedded lava flows. Late Paleocene and early Eocene organic matter suggests a highly oxidized marine environment, with major periods of deposition of terrestrially derived organic matter. These fluctuations in organic-matter type are probably the result of episodic shallowing and deepening of Rockall Basins. The final stage of Eocene/Oligocene sedimentation records the accelerating subsidence of Rockall and its isolation from terrestrial sources (Rockall and Greenland). This is shown by the increasingly marine character of the organic matter. The petroleum potential of sediments containing more than 0.5% organic carbon is poor because of their thermal immaturity and their highly oxidized and terrestrial organic-matter composition.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Oxygen concentration at time series station DYNAPROC

Continuous measurements between 0 and 200 m depth were performed every 2 h over two separate periods of four days at a station in the open northwestern Mediterranean Sea (Dyfamed Station) during the Dynaproc cruise in May 1995. Estimates of the daily variations in profiles of temperature, partial pressure of CO2, oxygen, chlorophyll a and nutrients were obtained. The distributions of the various physical and chemical properties were clearly different during the two time series, which were separated by a period of 11 days during which a wind event occurred. The mean daily utilization or production due to biological processes of dissolved inorganic carbon (DIC), nitrate+nitrite and oxygen were calculated along isopycnals using a vertical diffusion model. Between the surface and about 20 m depth, DIC was consumed and O2 released during the two time series while the nitrate+nitrite concentrations as well as supplies were zero. After the wind event, the O2 : C : N ratios of consumption (or production) were, on average, near the Redfield ratios, but during the first time series, the C : N utilization ratio between 20 and 35 m was two to three times that of Redfield stoichiometry and the oxygen release was low. The integrated net community production (NCP) in terms of carbon was equivalent during the two time series, whereas the chlorophyll a biomass was twice as high, on average, during the first time series but did decrease. These results imply that the production systems were different during the two periods. The first time series corresponds to a period at the end of production, due to the nutrient depletion in the euphotic layer. The formation of degradation products of the living material in dissolved organic form is probably important as indicated by the high C : N utilization ratios. The second time series corresponds to a reactivation of the primary production due to the upward shift of nutrients after the wind event.