Constraints on the onset duration of the Paleocene-Eocene Thermal Maximum.

The Paleocene-Eocene Thermal Maximum (PETM, approx. 56 Ma) provides a test case for investigating how the Earth system responds to rapid greenhouse gas-driven warming. However, current rates of carbon emissions are approximately 10 Pg C yr-1, whereas those proposed for the PETM span orders of magnitude-from ≪1 Pg C yr-1 to greater than the anthropogenic rate. Emissions rate estimates for the PETM are hampered by uncertainty over the total mass of PETM carbon released as well as the PETM onset duration. Here, I review constraints on the onset duration of the carbon isotope excursion (CIE) that is characteristic of the event with a focus on carbon cycle model-based attempts that forgo the need for a traditional sedimentary age model. I also review and compare existing PETM carbon input scenarios employing the Earth system model cGENIE and suggest another possibility-that abrupt input of an isotopically depleted carbon source combined with elevated volcanic outgassing over a longer interval can together account for key features of the PETM CIE.This article is part of a discussion meeting issue 'Hyperthermals: rapid and extreme global warming in our geological past'

Development of a novel empirical framework for interpreting geological carbon isotope excursions, with implications for the rate of carbon injection across the PETM

AbstractAs an episode of rapid global warming associated with the release of massive quantities of carbon to the atmosphere and oceans, the Paleocene–Eocene Thermal Maximum (PETM, ∼56 Ma) is considered a potential analog for modern anthropogenic carbon emissions. However, the prevailing order of magnitude uncertainty in the rate of carbon release during the PETM precludes any straightforward comparison between the paleo-record and the modern. Similar barriers exist to the interpretation of many other carbon isotope excursions in the geological record. Here we use the Earth system model cGENIE to quantify the consequences of differing carbon emissions rates on the isotopic record of different carbon reservoirs. We explore the consequences of a range of emissions scenarios – from durations of carbon input of years to millennia and constant versus pulsed emissions rates, and trace how the isotopic signal is imprinted on the different carbon reservoirs. From this, we identify a characteristic relationship between the difference in carbon isotope excursion sizes between atmospheric CO2 and dissolved inorganic carbon (DIC) and the duration of carbon emissions. To the extent that available isotopic data spanning the PETM constrain the size of the marine and atmospheric carbon isotopic excursions, applying this empirical relationship suggests the duration of the component of carbon emissions that dominates the isotopic signal could be less than 3000 yr. However, utilizing the ratio of excursion size in the atmosphere to ocean as a metric to constrain duration of carbon emissions highlights the necessity to strengthen estimates for these two measurements across the PETM. Our general interpretive framework could be equally applied in assessing rates of carbon emissions for other geological events

Earth System Model Analysis of How Astronomical Forcing Is Imprinted Onto the Marine Geological Record::The Role of the Inorganic (Carbonate) Carbon Cycle and Feedbacks

Astronomical cycles are strongly expressed in marine geological records, providing important insights into Earth system dynamics and an invaluable means of constructing age models. However, how various astronomical periods are filtered by the Earth system and the mechanisms by which carbon reservoirs and climate components respond, particularly in absence of dynamic ice sheets, is unclear. Using an Earth system model that includes feedbacks between climate, ocean circulation, and inorganic (carbonate) carbon cycling relevant to geological timescales, we systematically explore the impact of astronomically modulated insolation forcing and its expression in model variables most comparable to key paleoceanographic proxies (temperature, the δ13C of inorganic carbon, and sedimentary carbonate content). Temperature predominately responds to short and long eccentricity and is little influenced by the modeled carbon cycle feedbacks. In contrast, the cycling of nutrients and carbon in the ocean generates significant precession power in atmospheric CO2, benthic ocean δ13C, and sedimentary wt% CaCO3, while inclusion of marine sedimentary and weathering processes shifts power to the long eccentricity period. Our simulations produce reduced pCO2 and dissolved inorganic carbon (DIC) δ13C at long eccentricity maxima and, contrary to early Cenozoic marine records, CaCO3 preservation in the model is enhanced during eccentricity-modulated warmth. Additionally, the magnitude of δ13C variability simulated in our model underestimates marine proxy records. These model-data discrepancies hint at the possibility that the Paleogene silicate weathering feedback was weaker than modeled here and that additional organic carbon cycle feedbacks are necessary to explain the full response of the Earth system to astronomical forcing

Earth System Model Analysis of How Astronomical Forcing Is Imprinted Onto the Marine Geological Record::The Role of the Inorganic (Carbonate) Carbon Cycle and Feedbacks

Astronomical cycles are strongly expressed in marine geological records, providing important insights into Earth system dynamics and an invaluable means of constructing age models. However, how various astronomical periods are filtered by the Earth system and the mechanisms by which carbon reservoirs and climate components respond, particularly in absence of dynamic ice sheets, is unclear. Using an Earth system model that includes feedbacks between climate, ocean circulation, and inorganic (carbonate) carbon cycling relevant to geological timescales, we systematically explore the impact of astronomically modulated insolation forcing and its expression in model variables most comparable to key paleoceanographic proxies (temperature, the δ13C of inorganic carbon, and sedimentary carbonate content). Temperature predominately responds to short and long eccentricity and is little influenced by the modeled carbon cycle feedbacks. In contrast, the cycling of nutrients and carbon in the ocean generates significant precession power in atmospheric CO2, benthic ocean δ13C, and sedimentary wt% CaCO3, while inclusion of marine sedimentary and weathering processes shifts power to the long eccentricity period. Our simulations produce reduced pCO2 and dissolved inorganic carbon (DIC) δ13C at long eccentricity maxima and, contrary to early Cenozoic marine records, CaCO3 preservation in the model is enhanced during eccentricity-modulated warmth. Additionally, the magnitude of δ13C variability simulated in our model underestimates marine proxy records. These model-data discrepancies hint at the possibility that the Paleogene silicate weathering feedback was weaker than modeled here and that additional organic carbon cycle feedbacks are necessary to explain the full response of the Earth system to astronomical forcing

Astronomically paced changes in overturning circulation in the Western North Atlantic during the middle Eocene

North Atlantic Deep Water (NADW) currently redistributes heat and salt between Earth’s ocean basins, and plays a vital role in the ocean-atmosphere CO2 exchange. Despite its crucial role in today’s climate system, vigorous debate remains as to when deep-water formation in the North Atlantic started. Here, we present datasets from carbonate-rich middle Eocene sediments from the Newfoundland Ridge, revealing a unique archive of paleoceanographic change from the progressively cooling climate of the middle Eocene. Well-defined lithologic alternations between calcareous ooze and clay-rich intervals occur at the ∼41-kyr beat of axial obliquity. Hence, we identify obliquity as the driver of middle Eocene (43.5–46 Ma) Northern Component Water (NCW, the predecessor of modern NADW) variability. High-resolution benthic foraminiferal δ18O and δ13C suggest that obliquity minima correspond to cold, nutrient-depleted, western North Atlantic deep waters. We thus link stronger NCW formation with obliquity minima. In contrast, during obliquity maxima, Deep Western Boundary Currents were weaker and warmer, while abyssal nutrients were more abundant. These aspects reflect a more sluggish NCW formation. This obliquity-paced paleoceanographic regime is in excellent agreement with results from an Earth system model, in which obliquity minima configurations enhance NCW formation

System Architecture and Planetary Obliquity: Implications for Long-term Habitability

In the search for life beyond our solar system, attention should be focused on those planets that have the potential to maintain habitable conditions over the prolonged periods of time needed for the emergence and expansion of life as we know it. The observable planetary architecture is one of the determinants for long-term habitability as it controls the orbital evolution and ultimately the stellar fluxes received by the planet. With an ensemble of n-body simulations and obliquity models of hypothetical planetary systems, we demonstrate that the amplitude and period of the eccentricity, obliquity, and precession cycles of an Earth-like planet are sensitive to the orbital characteristics of a giant companion planet. A series of transient, ocean-coupled climate simulations show how these characteristics of astronomical cycles are decisive for the evolving surface conditions and long-term fractional habitability relative to the modern Earth. The habitability of Earth-like planets increases with the eccentricity of a Jupiter-like companion, provided that the mean obliquity is sufficiently low to maintain temperate temperatures over large parts of its surface throughout the orbital year. A giant companion closer in results in shorter eccentricity cycles of an Earth-like planet but longer, high-amplitude, obliquity cycles. The period and amplitude of obliquity cycles can be estimated to first order from the orbital pathways calculated by the n-body simulations. In the majority of simulations, the obliquity amplitude relates directly to the orbital inclination whereas the period of the obliquity cycle is a function of the nodal precession and the proximity of the giant companion

Assessing the impact of bioturbation on sedimentary isotopic records through numerical models

The disturbance of seafloor sediments by the activities of bottom-dwelling organisms, known as bioturbation, significantly alters the marine paleorecord by redistributing particles in the upper sediment layers. Consequently, ‘proxy’ signals recorded in these sediment particles, such as the size, abundance, or isotopic composition of plankton shells, are distorted by particle mixing. Accordingly, bioturbation can alter the apparent timing, duration, and magnitude of recorded events by smoothing climatic and oceanographic signals. In an extreme scenario, biological mixing can significantly obscure our view of the past by homogenizing the bioturbated layer, destroying sediment layering, and distorting the relative timing and intensity of past climatological events. Here we explore how bioturbation distorts proxy records of environmental events from a modeling perspective. First, we provide an overview and comparison of different numerical models created for simulating the movement and structural alteration of sediment by bioturbation. Next, we use an updated particle resolving model – iTURBO2 – to illustrate how various modes and intensities of bioturbation distort the signature of past climatological events, considering a range of conceptual shapes of vertical proxy profiles. Finally, we demonstrate how sampled proxy records can differ due to the combined effects of particle mixing and differential abundance changes that often concur with environmental transitions. We make the iTURBO2 MATLAB code openly available to facilitate further exploration of proxy biases due to bioturbation to aid the interpretation of the climatological record preserved in marine sediments

North Atlantic Drift Sediments Constrain Eocene Tidal Dissipation and the Evolution of the Earth-Moon System

Cyclostratigraphy and astrochronology are now at the forefront of geologic timekeeping. While this technique heavily relies on the accuracy of astronomical calculations, solar system chaos limits how far back astronomical calculations can be performed with confidence. High-resolution paleoclimate records with Milankovitch imprints now allow reversing the traditional cyclostratigraphic approach: Middle Eocene drift sediments from Newfoundland Ridge are well-suited for this purpose, due to high sedimentation rates and distinct lithological cycles. Per contra, the stratigraphies of Integrated Ocean Drilling Program Sites U1408–U1410 are highly complex with several hiatuses. Here, we built a two-site composite and constructed a conservative age-depth model to provide a reliable chronology for this rhythmic, highly resolved (<1 kyr) sedimentary archive. Astronomical components (g-terms and precession constant) are extracted from proxy time-series using two different techniques, producing consistent results. We find astronomical frequencies up to 4% lower than reported in astronomical solution La04. This solution, however, was smoothed over 20-Myr intervals, and our results therefore provide constraints on g-term variability on shorter, million-year timescales. We also report first evidence that the g4–g3 “grand eccentricity cycle” may have had a 1.2-Myr period around 41 Ma, contrary to its 2.4-Myr periodicity today. Our median precession constant estimate (51.28 ± 0.56″/year) confirms earlier indicators of a relatively low rate of tidal dissipation in the Paleogene. Newfoundland Ridge drift sediments thus enable a reliable reconstruction of astronomical components at the limit of validity of current astronomical calculations, extracted from geologic data, providing a new target for the next generation of astronomical calculations