skip to main content


Title: Atmospheric CO 2 over the Past 66 Million Years from Marine Archives
Throughout Earth's history, CO 2 is thought to have exerted a fundamental control on environmental change. Here we review and revise CO 2 reconstructions from boron isotopes in carbonates and carbon isotopes in organic matter over the Cenozoic—the past 66 million years. We find close coupling between CO 2 and climate throughout the Cenozoic, with peak CO 2 levels of ∼1,500 ppm in the Eocene greenhouse, decreasing to ∼500 ppm in the Miocene, and falling further into the ice age world of the Plio–Pleistocene. Around two-thirds of Cenozoic CO 2 drawdown is explained by an increase in the ratio of ocean alkalinity to dissolved inorganic carbon, likely linked to a change in the balance of weathering to outgassing, with the remaining one-third due to changing ocean temperature and major ion composition. Earth system climate sensitivity is explored and may vary between different time intervals. The Cenozoic CO 2 record highlights the truly geological scale of anthropogenic CO 2 change: Current CO 2 levels were last seen around 3 million years ago, and major cuts in emissions are required to prevent a return to the CO 2 levels of the Miocene or Eocene in the coming century. ▪  CO 2 reconstructions over the past 66 Myr from boron isotopes and alkenones are reviewed and re-evaluated. ▪  CO 2 estimates from the different proxies show close agreement, yielding a consistent picture of the evolution of the ocean-atmosphere CO 2 system over the Cenozoic. ▪  CO 2 and climate are coupled throughout the past 66 Myr, providing broad constraints on Earth system climate sensitivity. ▪  Twenty-first-century carbon emissions have the potential to return CO 2 to levels not seen since the much warmer climates of Earth's distant past.  more » « less
Award ID(s):
1806015
NSF-PAR ID:
10232685
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Annual Review of Earth and Planetary Sciences
Volume:
49
Issue:
1
ISSN:
0084-6597
Page Range / eLocation ID:
609 to 641
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. As atmospheric carbon dioxide (CO2) and temperatures increase with modern climate change, ancient hothouse periods become a focal point for understanding ecosystem function under similar conditions. The early Eocene exhibited high temperatures, high CO2 levels, and similar tectonic plate configuration as today, so it has been invoked as an analog to modern climate change. During the early Eocene, the greater Green River Basin (GGRB) of southwestern Wyoming was covered by an ancient hypersaline lake (Lake Gosiute; Green River Formation) and associated fluvial and floodplain systems (Wasatch and Bridger formations). The volcaniclastic Bridger Formation was deposited by an inland delta that drained from the northwest into freshwater Lake Gosiute and is known for its vast paleontological assemblages. Using this well-preserved basin deposited during a period of tectonic and paleoclimatic interest, we employ multiple proxies to study trends in provenance, parent material, weathering, and climate throughout 1 million years. The Blue Rim escarpment exposes approximately 100 m of the lower Bridger Formation, which includes plant and mammal fossils, solitary paleosol profiles, and organic remains suitable for geochemical analyses, as well as ash beds and volcaniclastic sandstone beds suitable for radioisotopic dating. New 40Ar/39Ar ages from the middle and top of the Blue Rim escarpment constrain the age of its strata to ∼ 49.5–48.5 Myr ago during the “falling limb” of the early Eocene Climatic Optimum. We used several geochemical tools to study provenance and parent material in both the paleosols and the associated sediments and found no change in sediment input source despite significant variation in sedimentary facies and organic carbon burial. We also reconstructed environmental conditions, including temperature, precipitation (both from paleosols), and the isotopic composition of atmospheric CO2 from plants found in the floral assemblages. Results from paleosol-based reconstructions were compared to semi-co-temporal reconstructions made using leaf physiognomic techniques and marine proxies. The paleosol-based reconstructions (near the base of the section) of precipitation (608–1167 mm yr−1) and temperature (10.4 to 12.0 ∘C) were within error of, although lower than, those based on floral assemblages, which were stratigraphically higher in the section and represented a highly preserved event later in time. Geochemistry and detrital feldspar geochronology indicate a consistent provenance for Blue Rim sediments, sourcing predominantly from the Idaho paleoriver, which drained the active Challis volcanic field. Thus, because there was neither significant climatic change nor significant provenance change, variation in sedimentary facies and organic carbon burial likely reflected localized geomorphic controls and the relative height of the water table. The ecosystem can be characterized as a wet, subtropical-like forest (i.e., paratropical) throughout the interval based upon the floral humidity province and Holdridge life zone schemes. Given the mid-paleolatitude position of the Blue Rim escarpment, those results are consistent with marine proxies that indicate that globally warm climatic conditions continued beyond the peak warm conditions of the early Eocene Climatic Optimum. The reconstructed atmospheric δ13C value (−5.3 ‰ to −5.8 ‰) closely matches the independently reconstructed value from marine microfossils (−5.4 ‰), which provides confidence in this reconstruction. Likewise, the isotopic composition reconstructed matches the mantle most closely (−5.4 ‰), agreeing with other postulations that warming was maintained by volcanic outgassing rather than a much more isotopically depleted source, such as methane hydrates. 
    more » « less
  2. Abstract. Since the middle Miocene (15 Ma, million years ago), the Earth's climate has undergone a long-term cooling trend, characterised by a reduction in ocean temperatures of up to 7–8 ∘C. The causes of this cooling are primarily thought to be due to tectonic plate movements driving changes in large-scale ocean circulation patterns, and hence heat redistribution, in conjunction with a drop in atmospheric greenhouse gas forcing (and attendant ice-sheet growth and feedback). In this study, we assess the potential to constrain the evolving patterns of global ocean circulation and cooling over the last 15 Ma by assimilating a variety of marine sediment proxy data in an Earth system model. We do this by first compiling surface and benthic ocean temperature and benthic carbon-13 (δ13C) data in a series of seven time slices spaced at approximately 2.5 Myr intervals. We then pair this with a corresponding series of tectonic and climate boundary condition reconstructions in the cGENIE (“muffin” release) Earth system model, including alternative possibilities for an open vs. closed Central American Seaway (CAS) from 10 Ma onwards. In the cGENIE model, we explore uncertainty in greenhouse gas forcing and the magnitude of North Pacific to North Atlantic salinity flux adjustment required in the model to create an Atlantic Meridional Overturning Circulation (AMOC) of a specific strength, via a series of 12 (one for each tectonic reconstruction) 2D parameter ensembles. Each ensemble member is then tested against the observed global temperature and benthic δ13C patterns. We identify that a relatively high CO2 equivalent forcing of 1120 ppm is required at 15 Ma in cGENIE to reproduce proxy temperature estimates in the model, noting that this CO2 forcing is dependent on the cGENIE model's climate sensitivity and that it incorporates the effects of all greenhouse gases. We find that reproducing the observed long-term cooling trend requires a progressively declining greenhouse gas forcing in the model. In parallel to this, the strength of the AMOC increases with time despite a reduction in the salinity of the surface North Atlantic over the cooling period, attributable to falling intensity of the hydrological cycle and to lowering polar temperatures, both caused by CO2-driven global cooling. We also find that a closed CAS from 10 Ma to present shows better agreement between benthic δ13C patterns and our particular series of model configurations and data. A final outcome of our analysis is a pronounced ca. 1.5 ‰ decline occurring in atmospheric (and ca. 1 ‰ ocean surface) δ13C that could be used to inform future δ13C-based proxy reconstructions. 
    more » « less
  3. Beaufort, Luc (Ed.)
    Abstract. The evolution of the Cenozoic cryosphere from unipolar to bipolar over the past 30 million years (Myr) is broadly known. Highly resolved records of carbonate (CaCO3) content provide insight into the evolution of regional and global climate, cryosphere, and carbon cycle dynamics. Here, we generate the first Southeast Atlantic CaCO3 content record spanning the last 30 Myr, derived from X-ray fluorescence (XRF) ln(Ca/Fe) data collected at Ocean Drilling Program Site 1264 (Walvis Ridge, SE Atlantic Ocean). We present a comprehensive and continuous depth and age model for the entirety of Site 1264 (~316 m; 30 Myr). This constitutes a key reference framework for future palaeoclimatic and palaeoceanographic studies at this location. We identify three phases with distinctly different orbital controls on Southeast Atlantic CaCO3 deposition, corresponding to major developments in climate, the cryosphere and the carbon cycle: (1) strong ~110 kyr eccentricity pacing prevails during Oligocene–Miocene global warmth (~30–13 Ma), (2) increased eccentricity-modulated precession pacing appears after the middle Miocene ClimateTransition (mMCT) (~14–8 Ma), and (3) pervasive obliquity pacing appears in the late Miocene (~7.7–3.3 Ma) following greater importance of high-latitude processes, such as increased glacial activity and high-latitude cooling. The lowest CaCO3 content (92 %–94 %) occurs between 18.5 and 14.5 Ma, potentially reflecting dissolution caused by widespread early Miocene warmth and preceding Antarctic deglaciation across the Miocene Climatic Optimum (~17–14.5 Ma) by 1.5 Myr. The emergence of precession pacing of CaCO3 deposition at Site 1264 after ~14 Ma could signal a reorganisation of surface and/or deep-water circulation in this region following Antarctic reglaciation at the mMCT. The increased sensitivity to precession at Site 1264 between 14 and 13 Ma is associated with an increase in mass accumulation rates (MARs) and reflects increased regional CaCO3 productivity and/or recurrent influxes of cooler, less corrosive deep waters. The highest carbonate content (%CaCO3) and MARs indicate that the late Miocene–early PlioceneBiogenic Bloom (LMBB) occurs between ~7.8 and 3.3Ma at Site 1264; broadly contemporaneous with the LMBB in the equatorial Pacific Ocean. At Site 1264, the onset of the LMBB roughly coincides with appearance of strong obliquity pacing of %CaCO3, reflecting increased high-latitude forcing. The global expression of the LMBB may reflect increased nutrient input into the global ocean resulting from enhanced aeolian dust and/or glacial/chemical weathering fluxes, due to enhanced glacial activity and increased meridional temperature gradients. Regional variability in the timing and amplitude of the LMBB may be driven by regional differences in cooling, continental aridification and/or changes in ocean circulation in the late Miocene. 
    more » « less
  4. null (Ed.)
    Prior to 2004, geological sampling in the Arctic Ocean was mainly restricted to near-surface Quaternary sediments. Thus, the long-term pre-Quaternary geological history is still poorly known. With the successful completion of the Arctic Coring Expedition (ACEX) (Integrated Ocean Drilling Program Expedition 302) in 2004, a new era in Arctic research began. Employing a novel multivessel approach, the first mission-specific platform (MSP) expedition of the Integrated Ocean Drilling Program proved that drilling in permanently ice-covered regions is possible. During ACEX, 428 m of Quaternary, Neogene, Paleogene, and Campanian sediment on Lomonosov Ridge were penetrated, providing new and unique insights into Cenozoic Arctic paleoceanographic and climate history. Although it was highly successful, ACEX also had three important limitations. The ACEX sequence contains either a large hiatus spanning the time interval from late Eocene to middle Miocene (based on the original biostratigraphic age model) or an interval of strongly reduced sedimentation rates (based on a more recent Os-Re-isotope-based age model). This is a critical time interval, spanning the time when prominent changes in global climate took place during the transition from the early Cenozoic Greenhouse Earth to the late Cenozoic Icehouse Earth. Furthermore, generally poor recovery during ACEX prevented detailed and continuous reconstruction of Cenozoic climate history. Finally, a higher-resolution reconstruction of Arctic rapid climate change during Neogene and Pleistocene times could not be achieved during ACEX. Therefore, Expedition 377 (Arctic Ocean Paleoceanography [ArcOP]) will return to the Lomonosov Ridge for a second MSP-type drilling campaign with the International Ocean Discovery Program to fill these major gaps in our knowledge on Arctic Ocean paleoenvironmental history through Cenozoic times and its relationship to global climate history. The overall goal of this drilling campaign is to recover a complete stratigraphic sedimentary record of the southern Lomonosov Ridge to meet our highest priority paleoceanographic objective, the continuous long-term Cenozoic climate history of the central Arctic Ocean. Furthermore, sedimentation rates two to four times higher than those of ACEX permit higher-resolution studies of Arctic climate change. The expedition goal can be achieved through careful site selection, the use of appropriate drilling technology and ice management, and by applying multiproxy approaches to paleoceanographic, paleoclimatic, and age-model reconstructions. The expedition will complete one primary deep drill hole (proposed Site LR-11B) to 900 meters below seafloor (mbsf), supplemented by a short drill site (LR-10B) to 50 mbsf, to recover an undisturbed uppermost (Quaternary) sedimentary section. This plan should ensure complete recovery so scientists can construct a composite section that spans the full age range through the Cenozoic. 
    more » « less
  5. null (Ed.)
    Abstract. Cenozoic stable carbon (δ13C) and oxygen (δ18O)isotope ratios of deep-sea foraminiferal calcite co-vary with the 405 kyreccentricity cycle, suggesting a link between orbital forcing, the climatesystem, and the carbon cycle. Variations in δ18O are partlyforced by ice-volume changes that have mostly occurred since the Oligocene.The cyclic δ13C–δ18O co-variation is found inboth ice-free and glaciated climate states, however. Consequently, thereshould be a mechanism that forces the δ13C cyclesindependently of ice dynamics. In search of this mechanism, we simulate theresponse of several key components of the carbon cycle to orbital forcing inthe Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir model(LOSCAR). We force the model by changing the burial of organic carbon in theocean with various astronomical solutions and noise and study the responseof the main carbon cycle tracers. Consistent with previous work, thesimulations reveal that low-frequency oscillations in the forcing arepreferentially amplified relative to higher frequencies. However, whileoceanic δ13C mainly varies with a 405 kyr period in themodel, the dynamics of dissolved inorganic carbon in the oceans and ofatmospheric CO2 are dominated by the 2.4 Myr cycle of eccentricity.This implies that the total ocean and atmosphere carbon inventory is stronglyinfluenced by carbon cycle variability that exceeds the timescale of the405 kyr period (such as silicate weathering). To test the applicability ofthe model results, we assemble a long (∼22 Myr) δ13C andδ18O composite record spanning the Eocene to Miocene(34–12 Ma) and perform spectral analysis to assess the presence of the2.4 Myr cycle. We find that, while the 2.4 Myr cycle appears to beovershadowed by long-term changes in the composite record, it is present asan amplitude modulator of the 405 and 100 kyr eccentricity cycles. 
    more » « less