Simulation of Eocene extreme warmth and high climate sensitivity through cloud feedbacks

Where did they come from, where did they go? Niche conservatism in woody and herbaceous plants and implications for plant-based paleoclimatic reconstructions

PREMISE

The ecological conditions that constrain plants to an environmental niche are assumed to be constant through time. While the fossil record has been used previously to test for niche conservatism of woody flowering plants, additional studies are needed in other plant groups especially since they can provide insight with paleoclimatic reconstructions, high biodiversity in modern terrestrial ecosystems, and significant contributions to agriculture.

METHODS

We tested climatic niche conservatism across time by characterizing the climatic niches of living herbaceous ginger plants (Zingiberaceae) and woody dawn redwood (Metasequoia) against paleoniches reconstructed based on fossil distribution data and paleoclimatic models.

RESULTS

Despite few fossil Zingiberaceae occurrences in the latitudinal tropics, unlike living Zingiberaceae, extinct Zingiberaceae likely experienced paratropical conditions in the higher latitudes, especially in the Cretaceous and Paleogene. The living and fossil distributions of Metasequoia largely remain in the upper latitudes of the northern hemisphere. The Zingiberaceae shifted from an initial subtropical climatic paleoniche in the Cretaceous, toward a temperate regime in the late Cenozoic; Metasequoia occupied a more consistent climatic niche over the same time intervals.

CONCLUSIONS

Because of the inconsistent climatic niches of Zingiberaceae over geologic time, we are less confident of using them for taxonomic-based paleoclimatic reconstruction methods like nearest living relative, which assume a consistent climatic niche between extant and extinct relatives; we argue that the consistent climatic niche of Metasequoia is more appropriate for these reconstructions. Niche conservatism cannot be assumed between extant and extinct plants and should be tested further in groups used for paleoclimatic reconstructions.

A 485-million-year history of Earth's surface temperature

A long-term record of global mean surface temperature (GMST) provides critical insight into the dynamical limits of Earth's climate and the complex feedbacks between temperature and the broader Earth system. Here, we present PhanDA, a reconstruction of GMST over the past 485 million years, generated by statistically integrating proxy data with climate model simulations. PhanDA exhibits a large range of GMST, spanning 11° to 36°C. Partitioning the reconstruction into climate states indicates that more time was spent in warmer rather than colder climates and reveals consistent latitudinal temperature gradients within each state. There is a strong correlation between atmospheric carbon dioxide (CO2) concentrations and GMST, identifying CO2 as the dominant control on variations in Phanerozoic global climate and suggesting an apparent Earth system sensitivity of ~8°C.

DeepMIP-Eocene-p1: multi-model dataset and interactive web application for Eocene climate research

Paleoclimate model simulations provide reference data to help interpret the geological record and offer a unique opportunity to evaluate the performance of current models under diverse boundary conditions. Here, we present a dataset of 35 climate model simulations of the warm early Eocene Climatic Optimum (EECO; ~ 50 million years ago) and corresponding preindustrial reference experiments. To streamline the use of the data, we apply standardised naming conventions and quality checks across eight modelling groups that have carried out coordinated simulations as part of the Deep-Time Model Intercomparison Project (DeepMIP). Gridded model fields can be downloaded from an online repository or accessed through a new web application that provides interactive data exploration. Local model data can be extracted in CSV format or visualised online for streamlined model-data comparisons. Additionally, processing and visualisation code templates may serve as a starting point for advanced analysis. The dataset and online platform aim to simplify accessing and handling complex data, prevent common processing issues, and facilitate the sharing of climate model data across disciplines.

The radiative feedback continuum from Snowball Earth to an ice-free hothouse

Paleoclimate records have been used to estimate the modern equilibrium climate sensitivity. However, this requires understanding how the feedbacks governing the climate response vary with the climate itself. Here we warm and cool a state-of-the-art climate model to simulate a continuum of climates ranging from a nearly ice-covered Snowball Earth to a nearly ice-free hothouse. We find that the pre-industrial (PI) climate is near a stability optimum: warming leads to a less-stable (more-sensitive) climate, as does cooling of more than 2K. Physically interpreting the results, we find that the decrease in stability for climates colder than the PI occurs mainly due to the albedo and lapse-rate feedbacks, and the decrease in stability for warmer climates occurs mainly due to the cloud feedback. These results imply that paleoclimate records provide a stronger constraint than has been calculated in previous studies, suggesting a reduction in the uncertainty range of the climate sensitivity.

Last Glacial Maximum pattern effects reduce climate sensitivity estimates

Here, we show that the Last Glacial Maximum (LGM) provides a stronger constraint on equilibrium climate sensitivity (ECS), the global warming from increasing greenhouse gases, after accounting for temperature patterns. Feedbacks governing ECS depend on spatial patterns of surface temperature ("pattern effects"); hence, using the LGM to constrain future warming requires quantifying how temperature patterns produce different feedbacks during LGM cooling versus modern-day warming. Combining data assimilation reconstructions with atmospheric models, we show that the climate is more sensitive to LGM forcing because ice sheets amplify extratropical cooling where feedbacks are destabilizing. Accounting for LGM pattern effects yields a median modern-day ECS of 2.4°C, 66% range 1.7° to 3.5°C (1.4° to 5.0°C, 5 to 95%), from LGM evidence alone. Combining the LGM with other lines of evidence, the best estimate becomes 2.9°C, 66% range 2.4° to 3.5°C (2.1° to 4.1°C, 5 to 95%), substantially narrowing uncertainty compared to recent assessments.

Climate system asymmetries drive eccentricity pacing of hydroclimate during the early Eocene greenhouse

The early Eocene Climatic Optimum (EECO) represents the peak of Earth's last sustained greenhouse climate interval. To investigate hydroclimate variability in western North America during the EECO, we developed an orbitally resolved leaf wax δ2H record from one of the most well-dated terrestrial paleoclimate archives, the Green River Formation. Our δ2Hwax results show ∼60‰ variation and evidence for eccentricity and precession forcing. iCESM simulations indicate that changes in the Earth's orbit drive large seasonal variations in precipitation and δ2H of precipitation at our study site, primarily during the summer season. Our findings suggest that the astronomical response in δ2Hwax is attributable to an asymmetrical climate response to the seasonal cycle, a "clipping" of precession forcing, and asymmetric carbon cycle dynamics, which further enhance the influence of eccentricity modulation on the hydrological cycle during the EECO. More broadly, our study provides an explanation for how and why eccentricity emerges as a dominant frequency in climate records from ice-free greenhouse worlds.

Persistent high latitude amplification of the Pacific Ocean over the past 10 million years

While high latitude amplification is seen in modern observations, paleoclimate records, and climate modeling, better constraints on the magnitude and pattern of amplification would provide insights into the mechanisms that drive it, which remain actively debated. Here we present multi-proxy multi-site paleotemperature records over the last 10 million years from the Western Pacific Warm Pool (WPWP) - the warmest endmember of the global ocean that is uniquely important in the global radiative feedback change. These sea surface temperature records, based on lipid biomarkers and seawater Mg/Ca-adjusted foraminiferal Mg/Ca, unequivocally show warmer WPWP in the past, and a secular cooling over the last 10 million years. Compiling these data with existing records reveals a persistent, nearly stationary, extratropical response pattern in the Pacific in which high latitude (~50°N) temperatures increase by ~2.4° for each degree of WPWP warming. This relative warming pattern is also evident in model outputs of millennium-long climate simulations with quadrupling atmospheric CO₂, therefore providing a strong constraint on the future equilibrium response of the Earth System.

Model spread in tropical low cloud feedback tied to overturning circulation response to warming

Among models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), here we show that the magnitude of the tropical low cloud feedback, which contributes considerably to uncertainty in estimates of climate sensitivity, is intimately linked to tropical deep convection and its effects on the tropical atmospheric overturning circulation. First, a reduction in tropical ascent area and an increased frequency of heavy precipitation result in high cloud reduction and upper-tropospheric drying, which increases longwave cooling and reduces subsidence weakening, favoring low cloud reduction (Radiation-Subsidence Pathway). Second, increased longwave cooling decreases tropospheric stability, which also reduces subsidence weakening and low cloudiness (Stability-Subsidence Pathway). In summary, greater high cloud reduction and upper-tropospheric drying (negative longwave feedback) lead to a more positive cloud feedback among CMIP6 models by contributing to a greater reduction in low cloudiness (positive shortwave feedback). Varying strengths of the two pathways contribute considerably to the intermodel spread in climate sensitivity.

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.

Past terrestrial hydroclimate sensitivity controlled by Earth system feedbacks

Despite tectonic conditions and atmospheric CO₂ levels (pCO₂) similar to those of present-day, geological reconstructions from the mid-Pliocene (3.3-3.0 Ma) document high lake levels in the Sahel and mesic conditions in subtropical Eurasia, suggesting drastic reorganizations of subtropical terrestrial hydroclimate during this interval. Here, using a compilation of proxy data and multi-model paleoclimate simulations, we show that the mid-Pliocene hydroclimate state is not driven by direct CO₂ radiative forcing but by a loss of northern high-latitude ice sheets and continental greening. These ice sheet and vegetation changes are long-term Earth system feedbacks to elevated pCO₂. Further, the moist conditions in the Sahel and subtropical Eurasia during the mid-Pliocene are a product of enhanced tropospheric humidity and a stationary wave response to the surface warming pattern, which varies strongly with land cover changes. These findings highlight the potential for amplified terrestrial hydroclimate responses over long timescales to a sustained CO₂ forcing.

In contrast to future projections, paleoclimate records often find wetter subtropics in tandem with elevated CO₂. Here, a compilation of proxies and simulations are used to reveal the climate dynamics and feedbacks responsible for generating wet subtropics during the mid-Pliocene.

The latitudinal temperature gradient and its climate dependence as inferred from foraminiferal δ 18 O over the past 95 million years

The temperature difference between low and high latitudes is one measure of the efficiency of the global climate system in redistributing heat and is used to test the ability of models to represent the climate system through time. Here, we show that the latitudinal temperature gradient has exhibited a consistent inverse relationship with global mean sea-surface temperature for at least the past 95 million years. Our results help reduce conflicts between climate models and empirical estimates of temperature and argue for a fundamental consistency in the dynamics of heat transport and radiative transfer across vastly different background states.

The latitudinal temperature gradient is a fundamental state parameter of the climate system tied to the dynamics of heat transport and radiative transfer. Thus, it is a primary target for temperature proxy reconstructions and global climate models. However, reconstructing the latitudinal temperature gradient in past climates remains challenging due to the scarcity of appropriate proxy records and large proxy–model disagreements. Here, we develop methods leveraging an extensive compilation of planktonic foraminifera δ¹⁸O to reconstruct a continuous record of the latitudinal sea-surface temperature (SST) gradient over the last 95 million years (My). We find that latitudinal SST gradients ranged from 26.5 to 15.3 °C over a mean global SST range of 15.3 to 32.5 °C, with the highest gradients during the coldest intervals of time. From this relationship, we calculate a polar amplification factor (PAF; the ratio of change in >60° S SST to change in global mean SST) of 1.44 ± 0.15. Our results are closer to model predictions than previous proxy-based estimates, primarily because δ¹⁸O-based high-latitude SST estimates more closely track benthic temperatures, yielding higher gradients. The consistent covariance of δ¹⁸O values in low- and high-latitude planktonic foraminifera and in benthic foraminifera, across numerous climate states, suggests a fundamental constraint on multiple aspects of the climate system, linking deep-sea temperatures, the latitudinal SST gradient, and global mean SSTs across large changes in atmospheric CO₂, continental configuration, oceanic gateways, and the extent of continental ice sheets. This implies an important underlying, internally driven predictability of the climate system in vastly different background states.

Globally resolved surface temperatures since the Last Glacial Maximum

Climate changes across the past 24,000 years provide key insights into Earth system responses to external forcing. Climate model simulations^1,2 and proxy data^3-8 have independently allowed for study of this crucial interval; however, they have at times yielded disparate conclusions. Here, we leverage both types of information using paleoclimate data assimilation^9,10 to produce the first proxy-constrained, full-field reanalysis of surface temperature change spanning the Last Glacial Maximum to present at 200-year resolution. We demonstrate that temperature variability across the past 24 thousand years was linked to two primary climatic mechanisms: radiative forcing from ice sheets and greenhouse gases; and a superposition of changes in the ocean overturning circulation and seasonal insolation. In contrast with previous proxy-based reconstructions^6,7 our results show that global mean temperature has slightly but steadily warmed, by ~0.5 °C, since the early Holocene (around 9 thousand years ago). When compared with recent temperature changes¹¹, our reanalysis indicates that both the rate and magnitude of modern warming are unusual relative to the changes of the past 24 thousand years.

Human Responses and Adaptation in a Changing Climate: A Framework Integrating Biological, Psychological, and Behavioural Aspects

Climate change is one of the biggest challenges of our times. Its impact on human populations is not yet completely understood. Many studies have focused on single aspects with contradictory observations. However, climate change is a complex phenomenon that cannot be adequately addressed from a single discipline's perspective. Hence, we propose a comprehensive conceptual framework on the relationships between climate change and human responses. This framework includes biological, psychological, and behavioural aspects and provides a multidisciplinary overview and critical information for focused interventions. The role of tipping points and regime shifts is explored, and a historical perspective is presented to describe the relationship between climate evolution and socio-cultural crisis. Vulnerability, resilience, and adaptation are analysed from an individual and a community point of view. Finally, emergent behaviours and mass effect phenomena are examined that account for mental maladjustment and conflicts.

Investigating Mesozoic Climate Trends and Sensitivities With a Large Ensemble of Climate Model Simulations

The Mesozoic era (∼252 to 66 million years ago) was a key interval in Earth's evolution toward its modern state, witnessing the breakup of the supercontinent Pangaea and significant biotic innovations like the early evolution of mammals. Plate tectonic dynamics drove a fundamental climatic transition from the early Mesozoic supercontinent toward the Late Cretaceous fragmented continental configuration. Here, key aspects of Mesozoic long-term environmental changes are assessed in a climate model ensemble framework. We analyze so far the most extended ensemble of equilibrium climate states simulated for evolving Mesozoic boundary conditions covering the period from 255 to 60 Ma in 5 Myr timesteps. Global mean temperatures are generally found to be elevated above the present and exhibit a baseline warming trend driven by rising sea levels and increasing solar luminosity. Warm (Triassic and mid-Cretaceous) and cool (Jurassic and end-Cretaceous) anomalies result from pCO2 changes indicated by different reconstructions. Seasonal and zonal temperature contrasts as well as continental aridity show an overall decrease from the Late Triassic-Early Jurassic to the Late Cretaceous. Meridional temperature gradients are reduced at higher global temperatures and less land area in the high latitudes. With systematic sensitivity experiments, the influence of paleogeography, sea level, vegetation patterns, pCO2, solar luminosity, and orbital configuration on these trends is investigated. For example, long-term seasonality trends are driven by paleogeography, but orbital cycles could have had similar-scale effects on shorter timescales. Global mean temperatures, continental humidity, and meridional temperature gradients are, however, also strongly affected by pCO2.

Probing the Ecology and Climate of the Eocene Southern Ocean With Sand Tiger Sharks Striatolamia macrota

Many explanations for Eocene climate change focus on the Southern Ocean-where tectonics influenced oceanic gateways, ocean circulation reduced heat transport, and greenhouse gas declines prompted glaciation. To date, few studies focus on marine vertebrates at high latitudes to discern paleoecological and paleoenvironmental impacts of this climate transition. The Tertiary Eocene La Meseta (TELM) Formation has a rich fossil assemblage to characterize these impacts; Striatolamia macrota, an extinct (†) sand tiger shark, is abundant throughout the La Meseta Formation. Body size is often tracked to characterize and integrate across multiple ecological dimensions. †S. macrota body size distributions indicate limited changes during TELMs 2-5 based on anterior tooth crown height (n = 450, mean = 19.6 ± 6.4 mm). Similarly, environmental conditions remained stable through this period based on δ18OPO4 values from tooth enameloid (n = 42; 21.5 ± 1.6‰), which corresponds to a mean temperature of 22.0 ± 4.0°C. Our preliminary ε Nd (n = 4) results indicate an early Drake Passage opening with Pacific inputs during TELM 2-3 (45-43 Ma) based on single unit variation with an overall radiogenic trend. Two possible hypotheses to explain these observations are (1) †S. macrota modified its migration behavior to ameliorate environmental changes related to the Drake Passage opening, or (2) the local climate change was small and gateway opening had little impact. While we cannot rule out an ecological explanation, a comparison with climate model results suggests that increased CO2 produces warm conditions that also parsimoniously explain the observations.

The enigma of Oligocene climate and global surface temperature evolution

During the Eocene, high-latitude regions were much warmer than today and substantial polar ice sheets were lacking. Indeed, the initiation of significant polar ice sheets near the end of the Eocene has been closely linked to global cooling. Here, we examine the relationship between global temperatures and continental-scale polar ice sheets following the establishment of ice sheets on Antarctica ∼34 million years ago, using records of surface temperatures from around the world. We find that high-latitude temperatures were almost as warm after the initiation of Antarctic glaciation as before, challenging our basic understanding of how climate works, and of the development of climate and ice volume through time.

Falling atmospheric CO₂ levels led to cooling through the Eocene and the expansion of Antarctic ice sheets close to their modern size near the beginning of the Oligocene, a period of poorly documented climate. Here, we present a record of climate evolution across the entire Oligocene (33.9 to 23.0 Ma) based on TEX₈₆ sea surface temperature (SST) estimates from southwestern Atlantic Deep Sea Drilling Project Site 516 (paleolatitude ∼36°S) and western equatorial Atlantic Ocean Drilling Project Site 929 (paleolatitude ∼0°), combined with a compilation of existing SST records and climate modeling. In this relatively low CO₂ Oligocene world (∼300 to 700 ppm), warm climates similar to those of the late Eocene continued with only brief interruptions, while the Antarctic ice sheet waxed and waned. SSTs are spatially heterogenous, but generally support late Oligocene warming coincident with declining atmospheric CO₂. This Oligocene warmth, especially at high latitudes, belies a simple relationship between climate and atmospheric CO₂ and/or ocean gateways, and is only partially explained by current climate models. Although the dominant climate drivers of this enigmatic Oligocene world remain unclear, our results help fill a gap in understanding past Cenozoic climates and the way long-term climate sensitivity responded to varying background climate states.

Polar amplification of Pliocene climate by elevated trace gas radiative forcing

Warm periods in Earth's history offer opportunities to understand the dynamics of the Earth system under conditions that are similar to those expected in the near future. The Middle Pliocene warm period (MPWP), from 3.3 to 3.0 My B.P, is the most recent time when atmospheric CO₂ levels were as high as today. However, climate model simulations of the Pliocene underestimate high-latitude warming that has been reconstructed from fossil pollen samples and other geological archives. One possible reason for this is that enhanced non-CO₂ trace gas radiative forcing during the Pliocene, including from methane (CH₄), has not been included in modeling. We use a suite of terrestrial biogeochemistry models forced with MPWP climate model simulations from four different climate models to produce a comprehensive reconstruction of the MPWP CH₄ cycle, including uncertainty. We simulate an atmospheric CH₄ mixing ratio of 1,000 to 1,200 ppbv, which in combination with estimates of radiative forcing from N₂O and O₃, contributes a non-CO₂ radiative forcing of 0.9 [Formula: see text] (range 0.6 to 1.1), which is 43% (range 36 to 56%) of the CO₂ radiative forcing used in MPWP climate simulations. This additional forcing would cause a global surface temperature increase of 0.6 to 1.0 °C, with amplified changes at high latitudes, improving agreement with geological evidence of Middle Pliocene climate. We conclude that natural trace gas feedbacks are critical for interpreting climate warmth during the Pliocene and potentially many other warm phases of the Cenezoic. These results also imply that using Pliocene CO₂ and temperature reconstructions alone may lead to overestimates of the fast or Charney climate sensitivity.

Proxy evidence for state-dependence of climate sensitivity in the Eocene greenhouse

Despite recent advances, the link between the evolution of atmospheric CO2 and climate during the Eocene greenhouse remains uncertain. In particular, modelling studies suggest that in order to achieve the global warmth that characterised the early Eocene, warmer climates must be more sensitive to CO2 forcing than colder climates. Here, we test this assertion in the geological record by combining a new high-resolution boron isotope-based CO2 record with novel estimates of Global Mean Temperature. We find that Equilibrium Climate Sensitivity (ECS) was indeed higher during the warmest intervals of the Eocene, agreeing well with recent model simulations, and declined through the Eocene as global climate cooled. These observations indicate that the canonical IPCC range of ECS (1.5 to 4.5 °C per doubling) is unlikely to be appropriate for high-CO2 warm climates of the past, and the state dependency of ECS may play an increasingly important role in determining the state of future climate as the Earth continues to warm.