The extent of ice growth during the Pliocene M2 glaciation (~3.3 Ma) has been called into question, with benthic foraminiferal oxygen isotope records interpreted primarily as a cooling signal. Here we improve the benthic foraminiferal Mg/Ca paleothermometer, allowing bottom water temperature reconstructions with a precision of ±0.2-0.3°C (1 s.d.). Applying this approach to M2 implies a significant increase in ice volume (~55 m SLE) that was more tightly coupled to a drop in CO2 than to ocean temperature. We suggest that the M2 glaciation was driven by a reduction in northern hemisphere poleward heat transport, and amplified by a reduction in southern hemisphere poleward heat transport caused by restriction of the Indonesian Seaway. The cryosphere growth drove the atmospheric CO2 decrease, which likely contributed to the overall magnitude of ice growth. These results demonstrate the sensitivity of the cryosphere to changes in ocean heat transport in a similar to modern climate.

AbstractBoron isotope ratios (δ11B) are used across the Earth Sciences and are increasing analyzed by Multi‐Collector Inductively Coupled Plasma Mass Spectrometry (MC‐ICPMS). Accurate δ11B MC‐ICPMS analysis requires boron purification from the sample matrix using ion‐exchange column chromatography. However, the traditional gravity‐drip column method is time‐consuming and prone to airborne contamination due to its long duration and open resin surface. To address these issues, we designed a novel, simple, and reliable column chromatography technique called “peri‐columns.” This method uses a peristaltic pump to generate vacuum on a commonly used column set up. This method uses sealed collection beakers and does not require solutions to pass through pump tubing, minimizing contamination. The duration is reduced by eight‐fold, processing 12 samples in just 1.5 hr. It also yields low and consistent total procedural blanks, averaging 11 pg. The efficiency and efficacy of this method were tested by repeated boron purification from calcium carbonate and high‐sodium matrices with international and in‐house reference materials. The results matched those obtained using the gravity column method and fell within our laboratory long‐term and international certified values. The mean δ11B and 2SD (standard deviation) of repeatedly processed NIST 8301f were 14.57 ± 0.26‰ (n = 31), NIST 8301c was 24.19 ± 0.33‰ (n = 10), STAiG‐F1 was 16.20 ± 0.26‰ (n = 13), and seawater was 39.52 ± 0.32‰ (n = 10). All the components of our techniques are commercially available, and it is easily adaptable to other laboratories and isotope systems.

Electron backscatter diffraction (EBSD) analysis enables a unique perspective of the internal microstructure of foraminiferal calcite. Specifically, EBSD provides crystallographic data from within the test, highlighting the highly organised “mesocrystal” structure of crystallographically aligned domains throughout the test, formed by sequential deposits of microgranular calcite. We compared EBSD maps across the test walls of both poorly preserved and well-preserved specimens of the planktonic foraminifera species Globigerinoides ruber and Morozovella crater. The EBSD maps, paired with information about intra-test distributions of Mg/Ca ratios, allowed us to examine the effects of different diagenetic processes on the foraminifera test. In poorly preserved specimens EBSD data show extensive reorganisation of the biogenic crystal microstructure, indicating differing phases of dissolution, re-precipitation and overgrowth. The specimens with the greatest degree of microstructural reorganisation also show an absence of higher concentration magnesium bands, which are typical features of well-preserved specimens. These findings provide important insights into the extent of post-depositional changes, in both microstructure and geochemical signals that must be considered when utilising foraminifera to generate proxy archive data.

Climate tipping elements are large-scale subsystems of the Earth that may transgress critical thresholds (tipping points) under ongoing global warming, with substantial impacts on the biosphere and human societies. Frequently studied examples of such tipping elements include the Greenland Ice Sheet, the Atlantic Meridional Overturning Circulation (AMOC), permafrost, monsoon systems, and the Amazon rainforest. While recent scientific efforts have improved our knowledge about individual tipping elements, the interactions between them are less well understood. Also, the potential of individual tipping events to induce additional tipping elsewhere or stabilize other tipping elements is largely unknown. Here, we map out the current state of the literature on the interactions between climate tipping elements and review the influences between them. To do so, we gathered evidence from model simulations, observations, and conceptual understanding, as well as examples of paleoclimate reconstructions where multi-component or spatially propagating transitions were potentially at play. While uncertainties are large, we find indications that many of the interactions between tipping elements are destabilizing. Therefore, we conclude that tipping elements should not only be studied in isolation, but also more emphasis has to be put on potential interactions. This means that tipping cascades cannot be ruled out on centennial to millennial timescales at global warming levels between 1.5 and 2.0 ∘C or on shorter timescales if global warming surpassed 2.0 ∘C. At these higher levels of global warming, tipping cascades may then include fast tipping elements such as the AMOC or the Amazon rainforest. To address crucial knowledge gaps in tipping element interactions, we propose four strategies combining observation-based approaches, Earth system modeling expertise, computational advances, and expert knowledge.

The geological record encodes the relationship between climate and atmospheric carbon dioxide (CO 2 ) over long and short timescales, as well as potential drivers of evolutionary transitions. However, reconstructing CO 2 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 CO 2 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 CO 2 thresholds in biological and cryosphere evolution.

The Mg/Ca ratio of foraminiferal calcite is a widely used empirical proxy for ocean temperature. Foraminiferal Mg/Ca-temperature relationships are based on extant species and are species-specific, introducing uncertainty when applying them to the fossil tests of extinct groups. Many modern species show remarkable heterogeneity in their intra-test Mg distributions, typically due to the presence of high Mg bands, which have a biological origin. Importantly, banding patterns differ between species, which could affect Mg/Ca-temperature relationships. Few studies have looked at intra-test variability in Mg/Ca ratios in extinct species of foraminifera, despite the obvious implications for paleothermometry. We used electron probe microanalysis (EPMA) to investigate intra-test Mg distributions in the fossil tests of two species of planktonic foraminifera from the extinct muricate mixed-layer-dwelling genus Morozovella, commonly used in Paleogene sea surface temperature reconstructions. Both M. aragonensis and M. crater show striking Mg banding patterns with multiple high and low Mg/Ca band pairs throughout the test wall in all chambers. The intra-test Mg variability in M. aragonensis and M. crater is similar to that in modern species widely used in paleoclimate reconstructions and banding patterns are consistent with published growth models for modern forms, albeit with subtle differences. The presence of Mg bands supports the application of Mg/Ca-palaeothermometry in extinct Morozovella species as well as the utility of EPMA for examining preservation of foraminifera tests in paleoclimatological studies. However, we emphasize the importance of rigorous assessments of inter-and intra-test Mg variability when using microanalytical techniques for foraminiferal Mg/Ca paleothermometry.

The timing and modalities of calcareous phytoplankton community and evolutionary responses to the Eocene-Oligocene transition (EOT, similar to 34 Ma) are still under-investigated. In order to better constrain the dynamics of these pelagic primary producers during the climate transition, we conducted high resolution assemblage analysis on calcareous nannofossils across a similar to 19 m-thick interval of nannofossil ooze at Ocean Drilling Program (ODP) Site 756 (Ninetyeast Ridge, Indian Ocean; Peirce et al. 1989) (paleolatitude similar to 43 degrees S; Zachos et al. 1992). We explored the diversity patterns against a new integrated planktonic foraminifera and calcareous nannofossil biostratigraphy produced for the site, as well as new benthic foraminifera and bulk sediment stable isotope (C, O) records, which documents ocean-climate changes, and provides independent chemostratigraphy. The study section spans nannofossil Zones NP20-NP23 (equivalent to CNE20-CNO4) and lasts 5.5 Myr. The results show that the hankeninid extinction falls within the similar to 4.5 m-thick EOT isotopic interval (0.67 m below the base of the second positive delta O-18 shift - EOIS), which is consistent with previous studies, making Hole 756C one of a few sites globally boasting both the familiar stepped delta O-18 and delta C-13 structure of the EOT and the primary biostratigraphic marker defining the base of the Oligocene. A series of potentially useful new calcareous nannofossil bioevents were identified that could help improve dating and correlation of this crucial interval. In this context, changes in calcareous nannofossil assemblages observed across EOT are interpreted in terms of modifications of paleoecological parameters that typically control the abundance and distribution of different taxa. Variations in sea surface temperature and nutrient availability are considered to be the most likely triggers for the calcareous phytoplankton changes observed across EOT. Specifically, our data suggest that increased nutrients in the mixed layer played a key role in shaping the late Eocene - early Oligocene calcareous nannofossil assemblages.

During the Miocene (23.0–5.3 Ma) North Africa experienced both humid and arid intervals, but the underlying cause of these transitions is unknown. Earth's climate was characterized by a unipolar icehouse with a dynamic Antarctic ice sheet, which may have influenced regional hydrology through atmospheric teleconnections. However, the Miocene also witnessed the restriction of the Mesopotamian Seaway, which may have had significant climatic impacts. The Maltese il‐Blata section (Central Mediterranean) comprises Late Oligocene to Early Miocene marine deposits previously used to constrain the timing of the Mesopotamian Seaway restriction using the εNd tracer. The location of this section also makes it sensitive to climatic changes in the North African region, and biogeochemical changes in the central Mediterranean. Here, we present lithological and geochemical records of the il‐Blata section. We find a marked shift in lithology and an increase in sedimentation rate coeval with the Early Miocene (∼19–20 Ma) restriction of the Mesopotamian Seaway. Concomitant changes in bulk sediment CaCO 3 , Sr/Ca, K/Al, Ti/Al, Zr/Al, and Si/Ti support a major humid climate transition and associated intensification of river systems over western North Africa. We propose that these changes in North African hydroclimate reflect either a tipping point effect in a gradually warming global climate, or are the result of the initial restriction of the Mesopotamian Seaway, perhaps through consequent changes in Atlantic Meridional Overturning Circulation and the West African Monsoon. We also suggest the restriction of the Mesopotamian Seaway inhibited phosphorite deposition at low latitudes.

Trace element and δ18O values of foraminifera are widely used to reconstruct oceanic temperatures throughout the Cenozoic and beyond. Previous work evaluating the geochemistry of foraminifera with differing degrees of physical preservation have shown that Mg/Ca and δ18O paleothermometers give discrepant values in recrystallized tests, with planktonic oxygen isotopes often yielding significantly lower temperatures than Mg/Ca ratios. To study the mobility of elements during diagenesis, we performed microspatial trace element analyses in Eocene Morozovella. Element maps show that trace element banding is readily identifiable and preserved, to an extent, in texturally recrystallized tests. A reaction-diffusion model was used to test whether the preservation of Mg-banding and the decoupling of δ18O and Mg/Ca values could be the result of diffusively limited “closed-system” recrystallization. Results show that, in a closed system, internal features (such as Mg-banding) will dissipate prior to changes in bulk Mg/Ca composition, while the bulk δ18O value will typically change faster than Mg/Ca. This is observed regardless of what partitioning coefficient is used for Mg and demonstrates that the planktonic Mg/Ca proxy is more diagenetically robust than the δ18O proxy. Thus, this model can explain the observed decoupling of these two proxies. Furthermore, the preservation of intra-test Mg-banding shows potential for use in evaluating the preservation of primary Mg/Ca values and hence the accuracy of paleotemperature reconstructions.

Muricate planktonic foraminifera comprise an extinct clade that was diverse and abundant in the Paleogene oceans and are widely used in palaeoclimate research as geochemical proxy carriers for the upper oceans. Their characteristic wall texture has surface projections called “muricae” formed by upward deflection and mounding of successive layers of the test wall. The group is generally considered to have lacked “true spines”: that is, acicular calcite crystals embedded in and projecting from the test surface such as occur in many modern and some Paleogene groups. Here we present evidence from polished sections, surface wall scanning electron microscope images and test dissections, showing that radially orientated crystalline spine-like structures occur in the centre of muricae in various species of Acarinina and Morozovella and projected from the test wall in life. Their morphology and placement in the wall suggest that they evolved independently of true spines. Nevertheless, they may have served a similar range of functions as spines in modern species, including aiding buoyancy and predation and especially harbouring algal photosymbionts, the function for which we suggest they probably first evolved. Our observations strengthen the analogy between Paleogene mixed-layer-dwelling planktonic foraminifera and their modern spinose counterparts.

Geology is the science of how the Earth functions and has evolved and, as such, it can contribute to our understanding of the climate system and how it responds to the addition of carbon dioxide (CO 2 ) to the atmosphere and oceans.Observations from the geological record show that atmospheric CO 2 concentrations are now at their highest levels in at least the past 3 million years.Furthermore, the current speed of human-induced CO 2 change and warming is nearly without precedent in the entire geological record, with the only known exception being the instantaneous, meteorite-induced event that caused the extinction of non-bird-like dinosaurs 66 million years ago.In short, whilst atmospheric CO 2 concentrations have varied dramatically during the geological past due to natural processes, and have often been higher than today, the current rate of CO 2 (and therefore temperature) change is unprecedented in almost the entire geological past.The geological record shows that changes in temperature and greenhouse gas concentrations have direct impacts on sea-level, the hydrological cycle, marine and terrestrial ecosystems, and the acidification and oxygen depletion of the oceans.Important climate phenomena, such as the El Niño-Southern Oscillation (ENSO) and the monsoons, which today affect the socio-economic stability and food and water security of billions of people, have varied markedly with past changes in climate.Climate reconstructions from around the globe show that climate change is not globally uniform, but tends to exhibit a consistent pattern, with changes at the poles larger than elsewhere.This polar amplification is seen in ancient warmer-than-modern time intervals like the Eocene epoch, about 50 million years ago and, more recently, in the Pliocene, about 3 million years ago.The warmest intervals of the Pliocene saw the disappearance of summer sea ice from the Arctic.The loss of ice cover during the Pliocene was one of the many rapid climate changes observed in the record, which are often called climate tipping points.The geological record can be used to calculate a quantity called Equilibrium Climate Sensitivity, which is the amount of warming caused by a doubling of atmospheric CO 2 , after various processes in the climate system have reached equilibrium.Recent estimates suggest that global mean climate warms between 2.6 and 3.9°C per doubling of CO 2 once all slow Earth system processes have reached equilibrium.The geological record provides powerful evidence that atmospheric CO 2 concentrations drive climate change, and supports multiple lines of evidence that greenhouse gases emitted by human activities are altering the Earth's climate.Moreover, the amount of anthropogenic greenhouse gases already in the atmosphere means that Earth is committed to a certain degree of warming.As the Earth's climate changes due to the burning of fossil fuels and changes in land-use, the planet we live on will experience further changes that will have increasingly drastic effects on human societies.An assessment of past climate changes helps to inform policy decisions regarding future climate change.Earth scientists will also have an important role to play in the delivery of any policies aimed at limiting future climate change.

The Miocene climate was dynamic, oscillating between major glaciation events and greenhouse conditions (the so-called Miocene Climatic Optimum or MCO). However, forcing factors responsible for climatic transitions from one state to another are not fully understood, partly because palaeoclimatological records from northern mid to high latitudes are scarce. To better resolve climatic changes of the Miocene epoch in the northern middle latitudes we studied a unique, nearly complete sedimentary record (Sdr. Vium borehole) spanning the upper Aquitanian to the Tortonian of the North Sea Basin. Newly obtained sea surface temperatures (SSTs) from our Miocene core revealed that the North Sea Basin was up to 20°C warmer than today, reaching the temperature maximum during the worldwide MCO (Herbert et al. 2020). Our high-resolution δ13C, TOC and C/N records, as well as elemental detrital ratios (Si/Al, Zr/Rb, Zr/Al) derived from XRF reveal important changes in the source of organic matter and detrital coarse fraction of the sediment. During the Miocene the location of the Sdr. Vium borehole was situated in a proximal setting, with water depths varying between 0 and ~200 m, partly due to advancing and retreating delta lobes and partly due to relative sea level changes. We observe that the depositional environment had a large impact on our record. By far the most important of these changes is a condensed interval associated with phosphatization, pyritization, and glauconite, associated with a major shift from a dark brown, organic-rich, bioturbated silty clay with thin sand lenses (the Hodde Formation) towards a green and brown clay with high concentrations of green glaucony pellets of fine sand grade (the Ørnhøj Formation). This shift is related to the subsidence of the North Sea Basin and marks the onset of a sediment-starvation in the basin.

The Miocene epoch, spanning 23.03–5.33 Ma, was a dynamic climate of sustained, polar amplified warmth. Miocene atmospheric CO2 concentrations are typically reconstructed between 300 and 600 ppm and were potentially higher during the Miocene Climatic Optimum (16.75–14.5 Ma). With surface temperature reconstructions pointing to substantial midlatitude and polar warmth, it is unclear what processes maintained the much weaker‐than‐modern equator‐to‐pole temperature difference. Here, we synthesize several Miocene climate modeling efforts together with available terrestrial and ocean surface temperature reconstructions. We evaluate the range of model‐data agreement, highlight robust mechanisms operating across Miocene modeling efforts and regions where differences across experiments result in a large spread in warming responses. Prescribed CO2 is the primary factor controlling global warming across the ensemble. On average, elements other than CO2, such as Miocene paleogeography and ice sheets, raise global mean temperature by ∼2°C, with the spread in warming under a given CO2 concentration (due to a combination of the spread in imposed boundary conditions and climate feedback strengths) equivalent to ∼1.2 times a CO2 doubling. This study uses an ensemble of opportunity: models, boundary conditions, and reference data sets represent the state‐of‐art for the Miocene, but are inhomogeneous and not ideal for a formal intermodel comparison effort. Acknowledging this caveat, this study is nevertheless the first Miocene multi‐model, multi‐proxy comparison attempted so far. This study serves to take stock of the current progress toward simulating Miocene warmth while isolating remaining challenges that may be well served by community‐led efforts to coordinate modeling and data activities within a common analytical framework.

The Eocene–Oligocene transition (EOT) was a climate shift from a largely ice-free greenhouse world to an icehouse climate, involving the first major glaciation of Antarctica and global cooling occurring ∼34 million years ago (Ma) and lasting ∼790 kyr. The change is marked by a global shift in deep-sea δ18O representing a combination of deep-ocean cooling and growth in land ice volume. At the same time, multiple independent proxies for ocean temperature indicate sea surface cooling, and major changes in global fauna and flora record a shift toward more cold-climate-adapted species. The two principal suggested explanations of this transition are a decline in atmospheric CO2 and changes to ocean gateways, while orbital forcing likely influenced the precise timing of the glaciation. Here we review and synthesise proxy evidence of palaeogeography, temperature, ice sheets, ocean circulation and CO2 change from the marine and terrestrial realms. Furthermore, we quantitatively compare proxy records of change to an ensemble of climate model simulations of temperature change across the EOT. The simulations compare three forcing mechanisms across the EOT: CO2 decrease, palaeogeographic changes and ice sheet growth. Our model ensemble results demonstrate the need for a global cooling mechanism beyond the imposition of an ice sheet or palaeogeographic changes. We find that CO2 forcing involving a large decrease in CO2 of ca. 40 % (∼325 ppm drop) provides the best fit to the available proxy evidence, with ice sheet and palaeogeographic changes playing a secondary role. While this large decrease is consistent with some CO2 proxy records (the extreme endmember of decrease), the positive feedback mechanisms on ice growth are so strong that a modest CO2 decrease beyond a critical threshold for ice sheet initiation is well capable of triggering rapid ice sheet growth. Thus, the amplitude of CO2 decrease signalled by our data–model comparison should be considered an upper estimate and perhaps artificially large, not least because the current generation of climate models do not include dynamic ice sheets and in some cases may be under-sensitive to CO2 forcing. The model ensemble also cannot exclude the possibility that palaeogeographic changes could have triggered a reduction in CO2.

Comparison between well-preserved and texturally recrystallized foraminifera suggests that δ 18 O values can change significantly, while trace-element (e.g. Mg/Ca) ratios may remain largely unchanged. EPMA analyses show moderate preservation of fine-scale features in trace element abundance such as Mg-banding, even in frosty planktonic foraminifera whose δ 18 O values are more similar to benthic temperatures; this process occurred in the absence of partial dissolution. This behavior has been modeled with a reaction-diffusion model, which simulates recrystallization with diffusively limited exchange with external water. This model assumes that recrystallization of the test is largely internal, and facilitated by fluid films, which also facilitate the transport of ions into the fossil. Due to the abundance of oxygen, and relative scarcity of Mg 2+ , Ca 2+ and Sr 2+ in seawater, δ 18 O values will change more rapidly in the carbonate phase. This “internal buffering” model can be applied to make testable predictions for other element and isotope ratios (such as Ca, Mg or Sr isotopes) which could be used to provide additional constraint on the effects of these diagenetic processes. These results show that micro-spatial analyses of foraminifera can be used to screen for numerous diagenetic processes, which can in turn be used to improve confidence in the use of different geochemical paleothermometers.

We present results from an ensemble of eight climate models, each of which has carried out simulations of the early Eocene climate optimum (EECO, ∼ 50 million years ago). These simulations have been carried out in the framework of the Deep-Time Model Intercomparison Project (DeepMIP; http://www.deepmip.org, last access: 10 January 2021); thus, all models have been configured with the same paleogeographic and vegetation boundary conditions. The results indicate that these non-CO2 boundary conditions contribute between 3 and 5 C to Eocene warmth. Compared with results from previous studies, the DeepMIP simulations generally show a reduced spread of the global mean surface temperature response across the ensemble for a given atmospheric CO2 concentration as well as an increased climate sensitivity on average. An energy balance analysis of the model ensemble indicates that global mean warming in the Eocene compared with the preindustrial period mostly arises from decreases in emissivity due to the elevated CO2 concentration (and associated water vapour and long-wave cloud feedbacks), whereas the reduction in the Eocene in terms of the meridional temperature gradient is primarily due to emissivity and albedo changes owing to the non-CO2 boundary conditions (i.e. the removal of the Antarctic ice sheet and changes in vegetation). Three of the models (the Community Earth System Model, CESM; the Geophysical Fluid Dynamics Laboratory, GFDL, model; and the Norwegian Earth System Model, NorESM) show results that are consistent with the proxies in terms of the global mean temperature, meridional SST gradient, and CO2, without prescribing changes to model parameters. In addition, many of the models agree well with the first-order spatial patterns in the SST proxies. However, at a more regional scale, the models lack skill. In particular, the modelled anomalies are substantially lower than those indicated by the proxies in the southwest Pacific; here, modelled continental surface air temperature anomalies are more consistent with surface air temperature proxies, implying a possible inconsistency between marine and terrestrial temperatures in either the proxies or models in this region. Our aim is that the documentation of the large-scale features and model–data comparison presented herein will pave the way to further studies that explore aspects of the model simulations in more detail, for example the ocean circulation, hydrological cycle, and modes of variability, and encourage sensitivity studies to aspects such as paleogeography, orbital configuration, and aerosols.