The PhanSST global database of Phanerozoic sea surface temperature proxy data

Paleotemperature proxy data form the cornerstone of paleoclimate research and are integral to understanding the evolution of the Earth system across the Phanerozoic Eon. Here, we present PhanSST, a database containing over 150,000 data points from five proxy systems that can be used to estimate past sea surface temperature. The geochemical data have a near-global spatial distribution and temporally span most of the Phanerozoic. Each proxy value is associated with consistent and queryable metadata fields, including information about the location, age, and taxonomy of the organism from which the data derive. To promote transparency and reproducibility, we include all available published data, regardless of interpreted preservation state or vital effects. However, we also provide expert-assigned diagenetic assessments, ecological and environmental flags, and other proxy-specific fields, which facilitate informed and responsible reuse of the database. The data are quality control checked and the foraminiferal taxonomy has been updated. PhanSST will serve as a valuable resource to the paleoclimate community and has myriad applications, including evolutionary, geochemical, diagenetic, and proxy calibration studies.

Ultra-depleted hydrogen isotopes in hydrated glass record Late Cretaceous glaciation in Antarctica

The Early Jurassic Butcher Ridge Igneous Complex (BRIC) in the Transantarctic Mountains contains abundant and variably hydrated silicic glass which has the potential to preserve a rich paleoclimate record. Here we present Fourier Transform Infrared Spectroscopic data that indicates BRIC glasses contain up to ~8 wt.% molecular water (H₂O_m), and low (<0.8 wt.%) hydroxyl (OH) component, interpreted as evidence for secondary hydration by meteoric water. BRIC glasses contain the most depleted hydrogen isotopes yet measured in terrestrial rocks, down to δD = −325 ‰. In situ ⁴⁰Ar/³⁹Ar geochronology of hydrated glasses with ultra-depleted δD values yield ages from 105 Ma to 72 Ma with a peak at c. 91.4 Ma. Combined, these data suggest hydration of BRIC glasses by polar glacial ice and melt water during the Late Cretaceous, contradicting paleoclimate reconstructions of this period that suggest Antarctica was ice-free and part of a global hot greenhouse.

Analysis of volcanic glass from the Transantarctic Mountains suggests that 90 Million years ago glaciation was widespread in Antarctica, a period in Earth’s history when the continent was considered to be ice-free and part of a global greenhouse.

Archaeal lipids trace ecology and evolution of marine ammonia-oxidizing archaea

Archaeal membrane lipids are widely used for paleotemperature reconstructions, yet these molecular fossils also bear rich information about ecology and evolution of marine ammonia-oxidizing archaea (AOA). Here we identified thermal and nonthermal behaviors of archaeal glycerol dialkyl glycerol tetraethers (GDGTs) by comparing the GDGT-based temperature index (TEX₈₆) to the ratio of GDGTs with two and three cyclopentane rings (GDGT-2/GDGT-3). Thermal-dependent biosynthesis should increase TEX₈₆ and decrease GDGT-2/GDGT-3 when the ambient temperature increases. This presumed temperature-dependent (PTD) trend is observed in GDGTs derived from cultures of thermophilic and mesophilic AOA. The distribution of GDGTs in suspended particulate matter (SPM) and sediments collected from above the pycnocline-shallow water samples-also follows the PTD trend. These similar GDGT distributions between AOA cultures and shallow water environmental samples reflect shallow ecotypes of marine AOA. While there are currently no cultures of deep AOA clades, GDGTs derived from deep water SPM and marine sediment samples exhibit nonthermal behavior deviating from the PTD trend. The presence of deep AOA increases the GDGT-2/GDGT-3 ratio and distorts the temperature-controlled correlation between GDGT-2/GDGT-3 and TEX₈₆. We then used Gaussian mixture models to statistically characterize these diagnostic patterns of modern AOA ecology from paleo-GDGT records to infer the evolution of marine AOA from the Mid-Mesozoic to the present. Long-term GDGT-2/GDGT-3 trends suggest a suppression of today's deep water marine AOA during the Mesozoic-early Cenozoic greenhouse climates. Our analysis provides invaluable insights into the evolutionary timeline and the expansion of AOA niches associated with major oceanographic and climate changes.

Late Cretaceous Paleoceanographic Evolution and the Onset of Cooling in the Santonian at Southern High Latitudes (IODP Site U1513, SE Indian Ocean)

The latest Cenomanian to Santonian sedimentary record recovered at IODP Expedition 369 Site U1513 in the Mentelle Basin (SE Indian Ocean, paleolatitude 60°S at 85 Ma) is studied to interpret the paleoceanographic evolution in the Southern Hemisphere. The planktonic foraminiferal assemblage changes, the depth ecology preferences of different species, and the surface and seafloor temperature inferred from the stable isotopic values measured on foraminiferal tests provide meaningful information to the understanding of the Late Cretaceous climate. The hothouse climate during the Turonian-Santonian, characterized by weak latitudinal temperature gradients and high atmospheric CO2 concentrations, is followed by a progressive cooling during the Campanian. At Site U1513 the beginning of this climatic transition is nicely recorded within the Santonian, as indicated by an ∼1‰ increase in δ18O values of planktonic foraminifera suggesting a decline in surface water paleotemperatures of 4°C. The onset of cooling is mirrored by changes in the planktonic foraminiferal assemblages including extinctions among surface and deep dwellers, appearances and diversification of newly evolving taxa, and changes from predominantly epifaunal oxic to infaunal dysoxic/suboxic taxa among co-occurring benthic foraminifera. Overall, the data presented here document an interval in the Santonian during which the rate of southern high latitude cooling increased. Both surface and bottom waters were affected, although the cooling signal is more evident in the data for surface waters. This pattern of cooling ascribes the deterioration of the Late Cretaceous climate to decreased CO2 in the atmosphere and changes in the oceanic circulation correlated with enhanced meridional circulation.

Early Cretaceous sea surface temperature evolution in subtropical shallow seas

Late Cretaceous sea surface temperatures (SST) are, amongst others, traditionally reconstructed by compiling oxygen isotope records of planktonic foraminifera obtained from globally distributed pelagic IODP drill cores. In contrast, the evolution of Early Cretaceous SSTs is essentially based on the organic TEX₈₆ palaeothermometer, as oxygen-isotope data derived from well-preserved 'glassy' foraminifer calcite are currently lacking. In order to evaluate the extraordinary warm TEX₈₆-derived SSTs of the Barremian to Aptian (130-123 Ma) subtropics, we present highly resolved sclerochemical profiles of pristine rudist bivalve shells from Tethyan and proto-North Atlantic shallow water carbonate platforms. An inverse correlation of seasonal ontogenetic variations in δ¹⁸O_rudist and Mg/Ca ratios demonstrates the fidelity of oxygen isotopes as palaeotemperature proxy. The new data shows moderate mean annual SSTs (22-26 °C) for large parts of the Barremian and Aptian and transient warm pulses for the so-called Mid-Barremian Event and Oceanic Anoxic Event 1a (reaching mean annual SSTs of 28 to 30 °C). A positive shift in mean annual oxygen-isotope values (δ¹⁸O: ≤ - 0.3‰) coupled with invariant Mg/Ca ratios at the Barremian-Aptian boundary points to a significant net loss of ¹⁶O in Tethyan shallow-marine settings. As the positive oxygen-isotope rudist shell values are recorded immediately beneath a major superregional hiatal surface, they are interpreted to be related to a major cooling phase and potential glacio-eustatic sea-level lowering. Our new sclerochemical findings are in clear contrast to open ocean SST records based on TEX₈₆, which indicate exceptionally warm Barremian to earliest Aptian subtropical oceans and weak meridional SST gradients.

A 313 plastome phylogenomic analysis of Pooideae: Exploring relationships among the largest subfamily of grasses

In this study, we analyzed 313 plastid genomes (plastomes) of Poaceae with a focus on expanding our current knowledge of relationships among the subfamily Pooideae, which represented over half the dataset (164 representatives). In total, 47 plastomes were sequenced and assembled for this study. This is the largest study of its kind to include plastome-level data, to not only increase sampling at both the taxonomic and molecular levels with the aim of resolving complex and reticulate relationships, but also to analyze the effects of alignment gaps in large-scale analyses, as well as explore divergences in the subfamily with an expanded set of 14 accepted grass fossils for more accurate calibrations and dating. Incorporating broad systematic assessments of Pooideae taxa conducted by authors within the last five years, we produced a robust phylogenomic reconstruction for the subfamily, which included all but two supergeneric taxa (Calothecinae and Duthieeae). We further explored how including alignment gaps in plastome analyses oftentimes can produce incorrect or misinterpretations of complex or reticulate relationships among taxa of Pooideae. This presented itself as consistently changing relationships at specific nodes for different stripping thresholds (percentage-based removal of gaps per alignment column). Our summary recommendation for large-scale genomic plastome datasets is to strip alignment columns of all gaps to increase pairwise identity and reduce errant signal from poly A/T bias. To do this we used the "mask alignment" tool in Geneious software. Finally, we determined an overall divergence age for Pooideae of roughly 84.8 Mya, which is in line with, but slightly older than most recent estimates.

Paleoceanography of the Late Cretaceous northwestern Tethys Ocean: Seasonal upwelling or steady thermocline?

In this study we attempted to assess whether seasonal upwelling or a steady thermocline persisted at the western margin of the Tethys Ocean during the late Turonian–early Coniacian interval. For this scope, we employed novel and published stable oxygen isotope (δ¹⁸O) data of various organisms (bivalves, bivalves, brachiopods, fish and belemnites). New seasonally resolved temperature estimates were based on the δ¹⁸O record of sequentially sampled inoceramid (Inoceramus sp.) and rudist (Hippurites resectus) shells from the Scaglia Rossa and Gosau deposits of northern Italy and western Austria, respectively. Diagenetic screening was performed using reflected light, cathodoluminescence (CL), scanning electron microscopy (SEM) and stable isotope analysis. Originally preserved δ¹³C and δ¹⁸O values were used to characterize the lifestyle of the bivalves and detect vital effects that could have biased oxygen isotope-based temperature reconstructions. Inoceramid δ¹⁸O values provide–for the first time–information on temperatures of Tethyan benthic waters, which were, on average, 14.4 ± 0.6 °C and fluctuated seasonally within a range of less than 2 °C. Such a thermal regime is in line with the temperatures postulated for late Turonian boreal water masses and support the existence of a cold water supply from the North Atlantic to the Tethyan bottom. Bottom cooling, however, did not affect the shallow water environment. In fact, the rudist-based temperature estimates for shallow water environment revealed a mean annual range of 11 °C, between 24 and 35 °C (assuming a seasonally constant δ18O_w = 1.0 ‰), which are among the warmest temperatures recorded over the entire Late Cretaceous. Our findings, thus, suggest a strong thermal and food web decoupling between the two environments. The absence of a seasonal vertical homogenization of different water bodies suggests the existence of a steady thermocline and, therefore, contrasts with the presence of an active coastal upwelling in the region as hypothesized by previous authors.

Temperate rainforests near the South Pole during peak Cretaceous warmth

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years^1-5, driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume⁶. In the near absence of proximal geological records from south of the Antarctic Circle, it is disputed whether polar ice could exist under such environmental conditions. Here we use a sedimentary sequence recovered from the West Antarctic shelf-the southernmost Cretaceous record reported so far-and show that a temperate lowland rainforest environment existed at a palaeolatitude of about 82° S during the Turonian-Santonian age (92 to 83 million years ago). This record contains an intact 3-metre-long network of in situ fossil roots embedded in a mudstone matrix containing diverse pollen and spores. A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric carbon dioxide concentrations of 1,120-1,680 parts per million by volume and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo under high levels of atmospheric carbon dioxide.

An early Cambrian greenhouse climate

The oceans of the early Cambrian (~541 to 509 million years ago) were the setting for a marked diversification of animal life. However, sea temperatures-a key component of the early Cambrian marine environment-remain unconstrained, in part because of a substantial time gap in the stable oxygen isotope (δ¹⁸O) record before the evolution of euconodonts. We show that previously overlooked sources of fossil biogenic phosphate have the potential to fill this gap. Pristine phosphatic microfossils from the Comley Limestones, UK, yield a robust δ¹⁸O signature, suggesting sea surface temperatures of 20° to 25°C at high southern paleolatitudes (~65°S to 70°S) between ~514 and 509 million years ago. These sea temperatures are consistent with the distribution of coeval evaporite and calcrete deposits, peak continental weathering rates, and also our climate model simulations for this interval. Our results support an early Cambrian greenhouse climate comparable to those of the late Mesozoic and early Cenozoic, offering a framework for exploring the interplay between biotic and environmental controls on Cambrian animal diversification.

Oxygen Isotopes in Foraminifera: Overview and Historical Review

Foraminiferal tests are a common component of many marine sediments. The oxygen isotope ratio (δ18O) of test calcite is frequently used to reconstruct aspects of their life environment. The δ18O depends mainly on the isotope ratio of the water it is precipitated from, the temperature of calcification, and, to a lesser extent, the carbonate ion concentration. Foraminifera and other organisms can potentially preserve their original isotope ratio for many millions of years, although diagenetic processes can alter the ratios. Work on oxygen isotope ratios of foraminifera was instrumental in the discovery of the orbital theory of the ice ages and continues to be widely used in the study of rapid climate change. Compilations of deep sea benthic foraminifer oxygen isotopes have revealed the long history of global climate change over the past 100 million years. Planktonic foraminifer oxygen isotopes are used to investigate the history of past sea surface temperatures, revealing the extent of past ‘greenhouse’ warming and global sea surface temperatures.

Cenozoic mean greenhouse gases and temperature changes with reference to the Anthropocene

Cenozoic greenhouse gases (GHG) variations and warming periods underscore the extreme rates of current climate change, with major implications for the adaptability and survivability of terrestrial and marine habitats. Current rise rate of greenhouse gases, reaching 3.3 ppm CO₂ per year during March 2015-2016, is the fastest recorded since the Paleocene-Eocene Thermal Event (PETM) when carbon release to the atmosphere was about an order of magnitude less than at present. The ice core evidence of concentration of (GHG) and temperatures in the atmosphere/ocean/cryosphere system over the last 740 kyr suggests that the rate of rise in GHG over the last ~260 years, CO₂ rates rising from 0.94 ppm yr^-1 in 1959 (315.97 ppm) to 1.62 ppm yr^-1 in 2000 (369.52 ppm) to 3.05 ppm yr^-1 in 2015 (400.83 ppm), constitutes a unique spike in the history of the atmosphere. The reliance of pre-740 kyr paleoclimate estimates on multiple proxies, including benthic and plankton fossils, fossil plants, residual organic matter, major and trace elements in fossils, sediments and soils, place limits on the resolution of pre-upper Pleistocene paleoclimate estimates, rendering it likely recorded mean Cenozoic paleoclimate trends may conceal abrupt short-term climate fluctuations. However, as exemplified by the Paleocene-Eocene thermal maximum (PETM) and earlier GHG and temperature spikes associated with major volcanic and asteroid impact events, the long-term residence time of CO₂ in the atmosphere extends the signatures of abrupt warming events to within detection limits of multiple paleoproxies. The mean post-1750 temperature rise rate (approximately ~0.0034 °C per yr, or ~0.008 °C per yr where temperature is not masked by sulfur aerosols) exceeds those of the PETM (approximately ~0.0008-0.0015 °C per yr) by an order of magnitude and mean glacial termination warming rates (last glacial termination [LGT] ~ 0.00039; Eemian ~0.0004 °C per yr) by near to an order of magnitude. Consistent with previous interglacial peaks an increasing likelihood of collapse of the Atlantic Meridional Ocean Circulation is threatening a severe stadial event.

Palaeogeographic regulation of glacial events during the Cretaceous supergreenhouse

The historical view of a uniformly warm Cretaceous is being increasingly challenged by the accumulation of new data hinting at the possibility of glacial events, even during the Cenomanian–Turonian (∼95 Myr ago), the warmest interval of the Cretaceous. Here we show that the palaeogeography typifying the Cenomanian–Turonian renders the Earth System resilient to glaciation with no perennial ice accumulation occurring under prescribed CO₂ levels as low as 420 p.p.m. Conversely, late Aptian (∼115 Myr ago) and Maastrichtian (∼70 Myr ago) continental configurations set the stage for cooler climatic conditions, favouring possible inception of Antarctic ice sheets under CO₂ concentrations, respectively, about 400 and 300 p.p.m. higher than for the Cenomanian–Turonian. Our simulations notably emphasize that palaeogeography can crucially impact global climate by modulating the CO₂ threshold for ice sheet inception and make the possibility of glacial events during the Cenomanian–Turonian unlikely.

Indirect evidence indicates the surprising occurrence of glacial events during the peak warmth of the Cretaceous world. Here, based on coupled climate-ice sheet model simulations, the authors show that such events were likely thwarted by palaeogeographic reorganisations and complex ocean-atmosphere feedbacks.

Eccentricity and obliquity paced carbon cycling in the Early Triassic and implications for post-extinction ecosystem recovery

The timing of marine ecosystem recovery following the End Permian Mass Extinction (EPME) remains poorly constrained given the lack of radiometric ages. Here we develop a high-resolution carbonate carbon isotope (δ¹³C_carb) record for 3.20 million years of the Olenekian in South China that defines the astronomical time-scale for the critical interval of major evolutionary and oceanic events in the Spathian. δ¹³C_carb documents eccentricity modulation of carbon cycling through the period and a strong obliquity signal. A shift in phasing between short and long eccentricity modulation, and amplification of obliquity, is nearly coincident with a 2% decrease in seawater δ¹³C_DIC, the last of a longer-term stepped decrease through the Spathian. The mid-Spathian shift in seawater δ¹³C_DIC to typical thermocline values is interpreted to record a major oceanic reorganization with global climate amelioration. Coincidence of the phasing shift with the first occurrence of marine reptiles (248.81 Ma), suggests that their invasion into the sea and the onset of a complex ecosystem were facilitated by restoration of deep ocean ventilation linked mechanistically to a change in the response of the oceanic carbon reservoir to astronomical forcing. Together these records place the first constraints on the duration of the post-extinction recovery to 3.35 myr.

Evidence for global cooling in the Late Cretaceous

The Late Cretaceous ‘greenhouse’ world witnessed a transition from one of the warmest climates of the past 140 million years to cooler conditions, yet still without significant continental ice. Low-latitude sea surface temperature (SST) records are a vital piece of evidence required to unravel the cause of Late Cretaceous cooling, but high-quality data remain illusive. Here, using an organic geochemical palaeothermometer (TEX₈₆), we present a record of SSTs for the Campanian–Maastrichtian interval (~83–66 Ma) from hemipelagic sediments deposited on the western North Atlantic shelf. Our record reveals that the North Atlantic at 35 °N was relatively warm in the earliest Campanian, with maximum SSTs of ~35 °C, but experienced significant cooling (~7 °C) after this to <~28 °C during the Maastrichtian. The overall stratigraphic trend is remarkably similar to records of high-latitude SSTs and bottom-water temperatures, suggesting that the cooling pattern was global rather than regional and, therefore, driven predominantly by declining atmospheric pCO₂ levels.

The Late Cretaceous experienced significant cooling, yet a lack of low-latitude records mean the regional extent of this cooling is poorly constrained. Linnert et al. present a TEX₈₆ sea surface temperature record from a palaeolatitude of ~35 °N and show that Late Cretaceous cooling was global in nature.

Geobiological constraints on Earth system sensitivity to CO₂ during the Cretaceous and Cenozoic

Earth system climate sensitivity (ESS) is the long-term (>10³ year) response of global surface temperature to doubled CO₂ that integrates fast and slow climate feedbacks. ESS has energy policy implications because global temperatures are not expected to decline appreciably for at least 10³ year, even if anthropogenic greenhouse gas emissions drop to zero. We report provisional ESS estimates of 3 °C or higher for some of the Cretaceous and Cenozoic based on paleo-reconstructions of CO₂ and temperature. These estimates are generally higher than climate sensitivities simulated from global climate models for the same ancient periods (approximately 3 °C). Climate models probably do not capture the full suite of positive climate feedbacks that amplify global temperatures during some globally warm periods, as well as other characteristic features of warm climates such as low meridional temperature gradients. These absent feedbacks may be related to clouds, trace greenhouse gases (GHGs), seasonal snow cover, and/or vegetation, especially in polar regions. Better characterization and quantification of these feedbacks is a priority given the current accumulation of atmospheric GHGs.

Spatial variations in archaeal lipids of surface water and core-top sediments in the South china sea and their implications for paleoclimate studies

The South China Sea (SCS) is the largest marginal sea of the western Pacific Ocean, yet little is known about archaeal distributions and TEX₈₆-based temperatures in this unique oceanic setting. Here we report findings of abundances in both core lipids (CL) and intact polar lipids (IPL) of Archaea from surface water (CL only) and core-top sediments from different regions of the SCS. TEX₈₆-derived temperatures were also calculated for these samples. The surface water had extremely low abundances of CL (average of 0.05 ± 0.13 ng/liter; n = 75), with higher values present in regions where upwelling is known to occur. The core-top sediments had CL values of 0.1 to 0.9 μg/g, which are on the low end of CL concentrations reported for other marine sediments and may reflect the oligotrophic nature of the open SCS. The IPL of Archaea accounted for 6 to 36.4% of total lipids (CL plus IPL), indicating that the majority of archaeal lipids in core-top sediments were derived from nonliving cells. The TEX₈₆-based temperatures of surface water were overall lower than satellite-based sea surface temperatures or CTD-measured in situ temperatures. The core-top sediment samples, however, had TEX₈₆ temperatures very close to the mean annual sea surface temperatures, except for samples with water depths of less than 100 m. Our results demonstrated low and heterogeneous distributions of archaeal lipids in surface water and core-top sediments of the SCS, which may reflect local or regional differences in productivity of Archaea. While TEX₈₆-based temperatures for core-top marine sediments at deep water depths (>100 m) generally reflected mean annual sea surface temperatures, TEX₈₆ temperatures in surface water varied basin wide and underestimated sea surface temperatures in most locations for the season when surface water samples were collected.

The magnitude and duration of Late Ordovician-Early Silurian glaciation

Understanding ancient climate changes is hampered by the inability to disentangle trends in ocean temperature from trends in continental ice volume. We used carbonate "clumped" isotope paleothermometry to constrain ocean temperatures, and thereby estimate ice volumes, through the Late Ordovician-Early Silurian glaciation. We find tropical ocean temperatures of 32° to 37°C except for short-lived cooling by ~5°C during the final Ordovician stage. Evidence for ice sheets spans much of the study interval, but the cooling pulse coincided with a glacial maximum during which ice volumes likely equaled or exceeded those of the last (Pleistocene) glacial maximum. This cooling also coincided with a large perturbation of the carbon cycle and the Late Ordovician mass extinction.