Sensitivity of ocean circulation to warming during the Early Eocene greenhouse.

Multiple abrupt warming events ("hyperthermals") punctuated the Early Eocene and were associated with deep-sea temperature increases of 2 to 4 °C, seafloor carbonate dissolution, and negative carbon isotope (δ13C) excursions. Whether hyperthermals were associated with changes in the global ocean overturning circulation is important for understanding their driving mechanisms and feedbacks and for gaining insight into the circulation's sensitivity to climatic warming. Here, we present high-resolution benthic foraminiferal stable isotope records (δ13C and δ18O) throughout the Early Eocene Climate Optimum (~53.26 to 49.14 Ma) from the deep equatorial and North Atlantic. Combined with existing records from the South Atlantic and Pacific, these indicate consistently amplified δ13C excursion sizes during hyperthermals in the deep equatorial Atlantic. We compare these observations with results from an intermediate complexity Earth system model to demonstrate that this spatial pattern of δ13C excursion size is a predictable consequence of global warming-induced changes in ocean overturning circulation. In our model, transient warming drives the weakening of Southern Ocean-sourced overturning circulation, strengthens Atlantic meridional water mass aging gradients, and amplifies the magnitude of negative δ13C excursions in the equatorial to North Atlantic. Based on model-data consistency, we conclude that Eocene hyperthermals coincided with repeated weakening of the global overturning circulation. Not accounting for ocean circulation impacts on δ13C excursions will lead to incorrect estimates of the magnitude of carbon release driving hyperthermals. Our finding of weakening overturning in response to past transient climatic warming is consistent with predictions of declining Atlantic Ocean overturning strength in our warm future.

Toward a Cenozoic history of atmospheric CO2.

The geological record encodes the relationship between climate and atmospheric carbon dioxide (CO2) over long and short timescales, as well as potential drivers of evolutionary transitions. However, reconstructing CO2 beyond direct measurements requires the use of paleoproxies and herein lies the challenge, as proxies differ in their assumptions, degree of understanding, and even reconstructed values. In this study, we critically evaluated, categorized, and integrated available proxies to create a high-fidelity and transparently constructed atmospheric CO2 record spanning the past 66 million years. This newly constructed record provides clearer evidence for higher Earth system sensitivity in the past and for the role of CO2 thresholds in biological and cryosphere evolution.

Why the Early Paleozoic was intrinsically prone to marine extinction.

The geological record of marine animal biodiversity reflects the interplay between changing rates of speciation versus extinction. Compared to mass extinctions, background extinctions have received little attention. To disentangle the different contributions of global climate state, continental configuration, and atmospheric oxygen concentration (pO2) to variations in background extinction rates, we drive an animal physiological model with the environmental outputs from an Earth system model across intervals spanning the past 541 million years. We find that climate and continental configuration combined to make extinction susceptibility an order of magnitude higher during the Early Paleozoic than during the rest of the Phanerozoic, consistent with extinction rates derived from paleontological databases. The high extinction susceptibility arises in the model from the limited geographical range of marine organisms. It stands even when assuming present-day pO2, suggesting that increasing oxygenation through the Paleozoic is not necessary to explain why extinction rates apparently declined with time.

What the geological past can tell us about the future of the ocean's twilight zone.

Paleontological reconstructions of plankton community structure during warm periods of the Cenozoic (last 66 million years) reveal that deep-dwelling 'twilight zone' (200-1000 m) plankton were less abundant and diverse, and lived much closer to the surface, than in colder, more recent climates. We suggest that this is a consequence of temperature's role in controlling the rate that sinking organic matter is broken down and metabolized by bacteria, a process that occurs faster at warmer temperatures. In a warmer ocean, a smaller fraction of organic matter reaches the ocean interior, affecting food supply and dissolved oxygen availability at depth. Using an Earth system model that has been evaluated against paleo observations, we illustrate how anthropogenic warming may impact future carbon cycling and twilight zone ecology. Our findings suggest that significant changes are already underway, and without strong emissions mitigation, widespread ecological disruption in the twilight zone is likely by 2100, with effects spanning millennia thereafter.

Volcanic CO2 degassing postdates thermogenic carbon emission during the end-Permian mass extinction.

Massive carbon dioxide (CO2) emissions are widely assumed to be the driver of the end-Permian mass extinction (EPME). However, the rate of and total CO2 released, and whether the source changes with time, remain poorly understood, leaving a key question surrounding the trigger for the EPME unanswered. Here, we assimilate reconstructions of atmospheric Pco2 and carbonate δ13C in an Earth system model to unravel the history of carbon emissions and sources across the EPME. We infer a transition from a CO2 source with a thermogenic carbon isotopic signature associated with a slower emission rate to a heavier, more mantle-dominated volcanic source with an increased rate of emissions. This implies that the CO2 degassing style changed as the Siberian Traps emplacement evolved, which is consistent with geochemical proxy records. Carbon cycle feedbacks from terrestrial ecosystem disturbances may have further amplified the warming and the severity of marine extinctions.

Spatial patterns of climate change across the Paleocene-Eocene Thermal Maximum.

The Paleocene-Eocene Thermal Maximum (PETM; 56 Ma) is one of our best geological analogs for understanding climate dynamics in a "greenhouse" world. However, proxy data representing the event are only available from select marine and terrestrial sedimentary sequences that are unevenly distributed across Earth's surface, limiting our view of the spatial patterns of climate change. Here, we use paleoclimate data assimilation (DA) to combine climate model and proxy information and create a spatially complete reconstruction of the PETM and the climate state that precedes it ("PETM-DA"). Our data-constrained results support strong polar amplification, which in the absence of an extensive cryosphere, is related to temperature feedbacks and loss of seasonal snow on land. The response of the hydrological cycle to PETM warming consists of a narrowing of the Intertropical Convergence Zone, off-equatorial drying, and an intensification of seasonal monsoons and winter storm tracks. Many of these features are also seen in simulations of future climate change under increasing anthropogenic emissions. Since the PETM-DA yields a spatially complete estimate of surface air temperature, it yields a rigorous estimate of global mean temperature change (5.6 ∘C; 5.4 ∘C to 5.9 ∘C, 95% CI) that can be used to calculate equilibrium climate sensitivity (ECS). We find that PETM ECS was 6.5 ∘C (5.7 ∘C to 7.4 ∘C, 95% CI), which is much higher than the present-day range. This supports the view that climate sensitivity increases substantially when greenhouse gas concentrations are high.

Continental configuration controls ocean oxygenation during the Phanerozoic.

The early evolutionary and much of the extinction history of marine animals is thought to be driven by changes in dissolved oxygen concentrations ([O2]) in the ocean1-3. In turn, [O2] is widely assumed to be dominated by the geological history of atmospheric oxygen (pO2)4,5. Here, by contrast, we show by means of a series of Earth system model experiments how continental rearrangement during the Phanerozoic Eon drives profound variations in ocean oxygenation and induces a fundamental decoupling in time between upper-ocean and benthic [O2]. We further identify the presence of state transitions in the global ocean circulation, which lead to extensive deep-ocean anoxia developing in the early Phanerozoic even under modern pO2. Our finding that ocean oxygenation oscillates over stable thousand-year (kyr) periods also provides a causal mechanism that might explain elevated rates of metazoan radiation and extinction during the early Palaeozoic Era6. The absence, in our modelling, of any simple correlation between global climate and ocean ventilation, and the occurrence of profound variations in ocean oxygenation independent of atmospheric pO2, presents a challenge to the interpretation of marine redox proxies, but also points to a hitherto unrecognized role for continental configuration in the evolution of the biosphere.

Post-extinction recovery of the Phanerozoic oceans and biodiversity hotspots.

The fossil record of marine invertebrates has long fuelled the debate as to whether or not there are limits to global diversity in the sea1-5. Ecological theory states that, as diversity grows and ecological niches are filled, the strengthening of biological interactions imposes limits on diversity6,7. However, the extent to which biological interactions have constrained the growth of diversity over evolutionary time remains an open question1-5,8-11. Here we present a regional diversification model that reproduces the main Phanerozoic eon trends in the global diversity of marine invertebrates after imposing mass extinctions. We find that the dynamics of global diversity are best described by a diversification model that operates widely within the exponential growth regime of a logistic function. A spatially resolved analysis of the ratio of diversity to carrying capacity reveals that less than 2% of the global flooded continental area throughout the Phanerozoic exhibits diversity levels approaching ecological saturation. We attribute the overall increase in global diversity during the Late Mesozoic and Cenozoic eras to the development of diversity hotspots under prolonged conditions of Earth system stability and maximum continental fragmentation. We call this the 'diversity hotspots hypothesis', which we propose as a non-mutually exclusive alternative to the hypothesis that the Mesozoic marine revolution led this macroevolutionary trend12,13.

Exploring the impact of climate change on the global distribution of non-spinose planktonic foraminifera using a trait-based ecosystem model.

Planktonic foraminifera are one of the primary calcifiers in the modern ocean, contributing 23%-56% of total global pelagic carbonate production. However, a mechanistic understanding of how physiology and environmental conditions control their abundance and distribution is lacking, hindering the projection of the impact of future climate change. This understanding is important, not only for ecosystem dynamics, but also for marine carbon cycling because of foraminifera's key role in carbonate production. Here we present and apply a global trait-based ecosystem model of non-spinose planktonic foraminifera ('ForamEcoGEnIE') to assess their ecology and global distribution under future climate change. ForamEcoGEnIE considers the traits of calcium carbonate production, shell size, and foraging. It captures the main characteristic of biogeographical patterns of non-spinose species - with maximum biomass concentrations found in mid- to high-latitude waters and upwelling areas. The model also reproduces the magnitude of global carbonate production relatively well, although the foraminifera standing stock is systematically overestimated. In response to future scenarios of rising atmospheric CO2 (RCP6 and RCP8.5), on a regional scale, the modelled foraminifera biomass and export flux increases in the subpolar regions of the North Atlantic and the Southern Ocean while it decreases everywhere else. In the absence of adaptation, the biomass decline in the low-latitude South Pacific suggests extirpation. The model projects a global average loss in non-spinose foraminifera biomass between 8% (RCP6) and 11% (RCP8.5) by 2050 and between 14% and 18% by 2100 as a response to ocean warming and associated changes in primary production and ecological dynamics. Global calcium carbonate flux associated with non-spinose foraminifera declines by 13%-18% by 2100. That decline can slow down the ocean carbonate pump and create short-term positive feedback on rising atmospheric pCO2 .

Decreasing Phanerozoic extinction intensity as a consequence of Earth surface oxygenation and metazoan ecophysiology.

The decline in background extinction rates of marine animals through geologic time is an established but unexplained feature of the Phanerozoic fossil record. There is also growing consensus that the ocean and atmosphere did not become oxygenated to near-modern levels until the mid-Paleozoic, coinciding with the onset of generally lower extinction rates. Physiological theory provides us with a possible causal link between these two observations-predicting that the synergistic impacts of oxygen and temperature on aerobic respiration would have made marine animals more vulnerable to ocean warming events during periods of limited surface oxygenation. Here, we evaluate the hypothesis that changes in surface oxygenation exerted a first-order control on extinction rates through the Phanerozoic using a combined Earth system and ecophysiological modeling approach. We find that although continental configuration, the efficiency of the biological carbon pump in the ocean, and initial climate state all impact the magnitude of modeled biodiversity loss across simulated warming events, atmospheric oxygen is the dominant predictor of extinction vulnerability, with metabolic habitat viability and global ecophysiotype extinction exhibiting inflection points around 40% of present atmospheric oxygen. Given this is the broad upper limit for estimates of early Paleozoic oxygen levels, our results are consistent with the relative frequency of high-magnitude extinction events (particularly those not included in the canonical big five mass extinctions) early in the Phanerozoic being a direct consequence of limited early Paleozoic oxygenation and temperature-dependent hypoxia responses.

Paleocene/Eocene carbon feedbacks triggered by volcanic activity.

The Paleocene-Eocene Thermal Maximum (PETM) was a period of geologically-rapid carbon release and global warming ~56 million years ago. Although modelling, outcrop and proxy records suggest volcanic carbon release occurred, it has not yet been possible to identify the PETM trigger, or if multiple reservoirs of carbon were involved. Here we report elevated levels of mercury relative to organic carbon-a proxy for volcanism-directly preceding and within the early PETM from two North Sea sedimentary cores, signifying pulsed volcanism from the North Atlantic Igneous Province likely provided the trigger and subsequently sustained elevated CO2. However, the PETM onset coincides with a mercury low, suggesting at least one other carbon reservoir released significant greenhouse gases in response to initial warming. Our results support the existence of 'tipping points' in the Earth system, which can trigger release of additional carbon reservoirs and drive Earth's climate into a hotter state.

Ecosystem function after the K/Pg extinction: decoupling of marine carbon pump and diversity.

The ocean biological pump is the mechanism by which carbon and nutrients are transported to depth. As such, the biological pump is critical in the partitioning of carbon dioxide between the ocean and atmosphere, and the rate at which that carbon can be sequestered through burial in marine sediments. How the structure and function of planktic ecosystems in the ocean govern the strength and efficiency of the biological pump and its resilience to disruption are poorly understood. The aftermath of the impact at the Cretaceous/Palaeogene (K/Pg) boundary provides an ideal opportunity to address these questions as both the biological pump and marine plankton size and diversity were fundamentally disrupted. The excellent fossil record of planktic foraminifera as indicators of pelagic-biotic recovery combined with carbon isotope records tracing biological pump behaviour, show that the recovery of ecological traits (diversity, size and photosymbiosis) occurred much later (approx. 4.3 Ma) than biological pump recovery (approx. 1.8 Ma). We interpret this decoupling of diversity and the biological pump as an indication that ecosystem function had sufficiently recovered to drive an effective biological pump, at least regionally in the South Atlantic.

Temperature controls carbon cycling and biological evolution in the ocean twilight zone.

Theory suggests that the ocean's biological carbon pump, the process by which organic matter is produced at the surface and transferred to the deep ocean, is sensitive to temperature because temperature controls photosynthesis and respiration rates. We applied a combined data-modeling approach to investigate carbon and nutrient recycling rates across the world ocean over the past 15 million years of global cooling. We found that the efficiency of the biological carbon pump increased with ocean cooling as the result of a temperature-dependent reduction in the rate of remineralization (degradation) of sinking organic matter. Increased food delivery at depth prompted the development of new deep-water niches, triggering deep plankton evolution and the expansion of the mesopelagic "twilight zone" ecosystem.

Orographic evolution of northern Tibet shaped vegetation and plant diversity in eastern Asia.

The growth of the Tibetan Plateau throughout the past 66 million years has profoundly affected the Asian climate, but how this unparalleled orogenesis might have driven vegetation and plant diversity changes in eastern Asia is poorly understood. We approach this question by integrating modeling results and fossil data. We show that growth of north and northeastern Tibet affects vegetation and, crucially, plant diversity in eastern Asia by altering the monsoon system. This northern Tibetan orographic change induces a precipitation increase, especially in the dry (winter) season, resulting in a transition from deciduous broadleaf vegetation to evergreen broadleaf vegetation and plant diversity increases across southeastern Asia. Further quantifying the complexity of Tibetan orographic change is critical for understanding the finer details of Asian vegetation and plant diversity evolution.

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.

Bistability in the redox chemistry of sediments and oceans.

For most of Earth's history, the ocean's interior was pervasively anoxic and showed occasional shifts in ocean redox chemistry between iron-buffered and sulfide-buffered states. These redox transitions are most often explained by large changes in external inputs, such as a strongly altered delivery of iron and sulfate to the ocean, or major shifts in marine productivity. Here, we propose that redox shifts can also arise from small perturbations that are amplified by nonlinear positive feedbacks within the internal iron and sulfur cycling of the ocean. Combining observational evidence with biogeochemical modeling, we show that both sedimentary and aquatic systems display intrinsic iron-sulfur bistability, which is tightly linked to the formation of reduced iron-sulfide minerals. The possibility of tipping points in the redox state of sediments and oceans, which allow large and nonreversible geochemical shifts to arise from relatively small changes in organic carbon input, has important implications for the interpretation of the geological rock record and the causes and consequences of major evolutionary transitions in the history of Earth's biosphere.

Past climates inform our future.

As the world warms, there is a profound need to improve projections of climate change. Although the latest Earth system models offer an unprecedented number of features, fundamental uncertainties continue to cloud our view of the future. Past climates provide the only opportunity to observe how the Earth system responds to high carbon dioxide, underlining a fundamental role for paleoclimatology in constraining future climate change. Here, we review the relevancy of paleoclimate information for climate prediction and discuss the prospects for emerging methodologies to further insights gained from past climates. Advances in proxy methods and interpretations pave the way for the use of past climates for model evaluation-a practice that we argue should be widely adopted.

Algal plankton turn to hunting to survive and recover from end-Cretaceous impact darkness.

The end-Cretaceous bolide impact triggered the devastation of marine ecosystems. However, the specific kill mechanism(s) are still debated, and how primary production subsequently recovered remains elusive. We used marine plankton microfossils and eco-evolutionary modeling to determine strategies for survival and recovery, finding that widespread phagotrophy (prey ingestion) was fundamental to plankton surviving the impact and also for the subsequent reestablishment of primary production. Ecological selectivity points to extreme post-impact light inhibition as the principal kill mechanism, with the marine food chain temporarily reset to a bacteria-dominated state. Subsequently, in a sunlit ocean inhabited by only rare survivor grazers but abundant small prey, it was mixotrophic nutrition (autotrophy and heterotrophy) and increasing cell sizes that enabled the eventual reestablishment of marine food webs some 2 million years later.

The impact of marine nutrient abundance on early eukaryotic ecosystems.

The rise of eukaryotes to ecological prominence represents one of the most dramatic shifts in the history of Earth's biosphere. However, there is an enigmatic temporal lag between the emergence of eukaryotic organisms in the fossil record and their much later ecological expansion. In parallel, there is evidence for a secular increase in the availability of the key macronutrient phosphorus (P) in Earth's oceans. Here, we use an Earth system model equipped with a size-structured marine ecosystem to explore relationships between plankton size, trophic complexity, and the availability of marine nutrients. We find a strong dependence of planktonic ecosystem structure on ocean nutrient abundance, with a larger ocean nutrient inventory leading to greater overall biomass, broader size spectra, and increasing abundance of large Zooplankton. If existing estimates of Proterozoic marine nutrient levels are correct, our results suggest that increases in the ecological impact of eukaryotic algae and trophic complexity in eukaryotic ecosystems were directly linked to restructuring of the global P cycle associated with the protracted rise of surface oxygen levels. Our results thus suggest an indirect but potentially important mechanism by which ocean oxygenation may have acted to shape marine ecological function during late Proterozoic time.

Fundamentally different global marine nitrogen cycling in response to severe ocean deoxygenation.

The present-day marine nitrogen (N) cycle is strongly regulated by biology. Deficiencies in the availability of fixed and readily bioavailable nitrogen relative to phosphate (P) in the surface ocean are largely corrected by the activity of diazotrophs. This feedback system, termed the "nitrostat," is thought to have provided close regulation of fixed-N speciation and inventory relative to P since the Proterozoic. In contrast, during intervals of intense deoxygenation such as Cretaceous ocean anoxic event (OAE) 2, a few regional sedimentary δ15N records hint at the existence of a different mode of marine N cycling in which ammonium plays a major role in regulating export production. However, the global-scale dynamics during this time remain unknown. Here, using an Earth System model and taking the example of OAE 2, we provide insights into the global marine nitrogen cycle under severe ocean deoxygenation. Specifically, we find that the ocean can exhibit fundamental transitions in the species of nitrogen dominating the fixed-N inventory--from nitrate (NO3-) to ammonium (NH4+)--and that as this transition occurs, the inventory can partially collapse relative to P due to progressive spatial decoupling between the loci of NH4+ oxidation, NO3- reduction, and nitrogen fixation. This finding is relatively independent of the specific state of ocean circulation and is consistent with nitrogen isotope and redox proxy data. The substantive reduction in the ocean fixed-N inventory at an intermediate state of deoxygenation may represent a biogeochemical vulnerability with potential implications for past and future (warmer) oceans.