skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The 405 kyr and 2.4 Myr eccentricity components in Cenozoic carbon isotope records
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
Award ID(s):
1658023
PAR ID:
10213257
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Climate of the Past
Volume:
15
Issue:
1
ISSN:
1814-9332
Page Range / eLocation ID:
91 to 104
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Climate variability is typically amplified towards polar regions. The underlying causes, notably albedo and humidity changes, are challenging to accurately quantify with observations or models, thus hampering projections of future polar amplification. Polar amplification reconstructions from the ice-free early Eocene (∼56–48 Ma) can exclude ice albedo effects, but the required tropical temperature records for resolving timescales shorter than multi-million years are lacking. Here, we reconstruct early Eocene tropical sea surface temperature variability by presenting an up to ∼4 kyr resolution biomarker-based temperature record from Ocean Drilling Program (ODP) Site 959, located in the tropical Atlantic Ocean. This record shows warming across multiple orbitally paced carbon cycle perturbations, coeval with high-latitude-derived deep-ocean bottom waters, showing that these events represent transient global warming events (hyperthermals). This implies that orbital forcing caused global temperature variability through carbon cycle feedbacks. Importantly, deep-ocean temperature variability was amplified by a factor of 1.7–2.3 compared to the tropical surface ocean, corroborating available long-term estimates. This implies that fast atmospheric feedback processes controlled meridional temperature gradients on multi-million year, as well as orbital, timescales during the early Eocene. Our combined records have several other implications. First, our amplification factor is somewhat larger than the same metric in fully coupled simulations of the early Eocene (1.1–1.3), suggesting that models slightly underestimate the non-ice-related – notably hydrological – feedbacks that cause polar amplification of climate change. Second, even outside the hyperthermals, we find synchronous eccentricity-forced temperature variability in the tropics and deep ocean that represent global mean sea surface temperature variability of up to 0.7 °C, which requires significant variability in atmospheric pCO2. We hypothesize that the responsible carbon cycle feedbacks that are independent of ice, snow, and frost-related processes might play an important role in Phanerozoic orbital-scale climate variability throughout geological time, including Pleistocene glacial–interglacial climate variability. 
    more » « less
  2. Abstract Astronomical (or Milanković) forcing of the Earth system is key to understanding rhythmic climate change on time scales ≳104 y. Paleoceanographic and paleoclimatological applications concerned with past astronomical forcing rely on astronomical calculations (solutions), which represent the backbone of cyclostratigraphy and astrochronology. Here we present state‐of‐the‐art astronomical solutions over the past 3.5 Gyr. Our goal is to provide tuning targets and templates for interpreting deep‐time cyclostratigraphic records and designing external forcing functions in climate models. Our approach yields internally consistent orbital and precession‐tilt solutions, including fundamental solar system frequencies, orbital eccentricity and inclination, lunar distance, luni‐solar precession rate, Earth's obliquity, and climatic precession. Contrary to expectations, we find that the long eccentricity cycle (LEC) (previously assumed stable and labeled “metronome,” recent period ∼405 kyr), can become unstable on long time scales. Our results reveal episodes during which the LEC is very weak or absent and Earth's orbital eccentricity and climate‐forcing spectrum are unrecognizable compared to the recent past. For the ratio of eccentricity‐to‐inclination amplitude modulation (recent individual periods of ~2.4 and ~1.2 Myr, frequently observable in paleorecords) we find a wide distribution around the recent 2:1 ratio, that is, the system is not restricted to a 2:1 or 1:1 resonance state. Our computations show that Earth's obliquity was lower and its amplitude (variation around the mean) significantly reduced in the past. We therefore predict weaker climate forcing at obliquity frequencies in deep time and a trend toward reduced obliquity power with age in stratigraphic records. For deep‐time stratigraphic and modeling applications, the orbital parameters of our 3.5‐Gyr integrations are made available at 400‐year resolution. 
    more » « less
  3. Abstract. Paleoceanographic interpretations of Plio-Pleistocene climate variability over the past 5 million years rely on the evaluation of event timing of proxy changes in sparse records across multiple ocean basins. In turn, orbital-scale chronostratigraphic controls for these records are often built from stratigraphic alignment of benthic foraminiferal stable oxygen isotope (δ18O) records to a preferred dated target stack or composite. This chronostratigraphic age model approach yields age model uncertainties associated with alignment method, target selection, the assumption that the undated record and target experienced synchronous changes in benthic foraminiferal δ18O values, and the assumption that any possible stratigraphic discontinuities within the undated record have been appropriately identified. However, these age model uncertainties and their impact on paleoceanographic interpretations are seldom reported or discussed. Here, we investigate and discuss these uncertainties for conventional manual and automated tuning techniques based on benthic foraminiferal δ18O records and evaluate their impact on sedimentary age models over the past 3.5 Myr using three sedimentary benthic foraminiferal δ18O records as case studies. In one case study, we present a new benthic foraminiferal δ18O record for International Ocean Discovery Program (IODP) Site U1541 (54°13′ S, 125°25′ W), recently recovered from the South Pacific on IODP Expedition 383. The other two case studies examine published benthic foraminiferal δ18O records of Ocean Drilling Program (ODP) Site 1090 and the ODP Site 980/981 composite. Our analysis suggests average age uncertainties of 3 to 5 kyr associated with manually derived versus automated alignment, 1 to 3 kyr associated with automated probabilistic alignment itself, and 2 to 6 kyr associated with the choice of tuning target. Age uncertainties are higher near stratigraphic segment ends and where local benthic foraminiferal δ18O stratigraphy differs from the tuning target. We conclude with recommendations for community best practices for the development and characterization of age uncertainty of sediment core chronostratigraphies based on benthic foraminiferal δ18O records. 
    more » « less
  4. 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
  5. The late Paleocene and early Eocene (LPEE) are characterized by long-term (million years, Myr) global warming and by transient, abrupt (kiloyears, kyr) warming events, termed hyperthermals. Although both have been attributed to greenhouse (CO2) forcing, the longer-term trend in climate was likely influenced by additional forcing factors (i.e., tectonics) and the extent to which warming was driven by atmospheric CO2remains unclear. Here, we use a suite of new and existing observations from planktic foraminifera collected at Pacific Ocean Drilling Program Sites 1209 and 1210 and inversion of a multiproxy Bayesian hierarchical model to quantify sea surface temperature (SST) and atmospheric CO2over a 6-Myr interval. Our reconstructions span the initiation of long-term LPEE warming (~58 Ma), and the two largest Paleogene hyperthermals, the Paleocene–Eocene Thermal Maximum (PETM, ~56 Ma) and Eocene Thermal Maximum 2 (ETM-2, ~54 Ma). Our results show strong coupling between CO2and temperature over the long- (LPEE) and short-term (PETM and ETM-2) but differing Pacific climate sensitivities over the two timescales. Combined CO2and carbon isotope trends imply the carbon source driving CO2increase was likely methanogenic, organic, or mixed for the PETM and organic for ETM-2, whereas a source with higher δ13C values (e.g., volcanic degassing) is associated with the long-term LPEE. Reconstructed emissions for the PETM (5,800 Gt C) and ETM-2 (3,800 Gt C) are comparable in mass to future emission scenarios, reinforcing the value of these events as analogs of anthropogenic change. 
    more » « less