Search

Earth's hydrological cycle is expected to intensify in response to global warming, with a “wet-gets-wetter, dry-gets-drier” response anticipated over the ocean. Subtropical regions (∼15°–30°N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data-modeling approach to reconstruct global and zonal-mean rainfall patterns during the early Eocene (∼56–48 million years ago). The Deep-Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid- (30°–60°N/S) and high-latitudes (>60°N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0°–15°N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter-Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation-evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter-model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy-derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation-induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patterns.

Estimates of global mean near-surface air temperature (global SAT) for the Cenozoic era rely largely on paleo-proxy data of deep-sea temperature (DST), with the assumption that changes in global SAT covary with changes in the global mean deep-sea temperature (global DST) and global mean sea-surface temperature (global SST). We tested the validity of this assumption by analyzing the relationship between global SST, SAT, and DST using 25 different model simulations from the Deep-Time Model Intercomparison Project simulating the early Eocene Climatic Optimum (EECO) with varying CO2 levels. Similar to the modern situation, we find limited spatial variability in DST, indicating that local DST estimates can be regarded as a first order representative of global DST. In line with previously assumed relationships, linear regression analysis indicates that both global DST and SAT respond stronger to changes in atmospheric CO2 than global SST by a similar factor. Consequently, this model-based analysis validates the assumption that changes in global DST can be used to estimate changes in global SAT during the early Cenozoic. Paleo-proxy estimates of global DST, SST, and SAT during EECO show the best fit with model simulations with a 1,680 ppm atmospheric CO2 level. This matches paleo-proxies of EECO atmospheric CO2, indicating a good fit between models and proxy-data.

Earth's hydrological cycle is expected to intensify in response to global warming, with a wet-gets-wetter, dry-gets-drier response anticipated over the ocean. Subtropical regions (similar to 15 degrees-30 degrees N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data-modeling approach to reconstruct global and zonal-mean rainfall patterns during the early Eocene (similar to 56-48 million years ago). The Deep-Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid-(30 degrees-60 degrees N/S) and high-latitudes (>60 degrees N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0 degrees-15 degrees N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter-Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation-evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter-model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy-derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation-induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patte

The total meridional heat transport (MHT) is relatively stable across different climates. Nevertheless, the strength of individual processes contributing to the total transport are not stable. Here we investigate the MHT and its main components especially in the atmosphere, in five coupled climate model simulations from the Deep-Time Model Intercomparison Project (DeepMIP). These simulations target the early Eocene climatic optimum, a geological time period with high CO2 concentrations, analog to the upper range of end-of-century CO2 projections. Preindustrial and early Eocene simulations, at a range of CO2 levels are used to quantify the MHT changes in response to both CO2 and non-CO2 related forcings. We found that atmospheric poleward heat transport increases with CO2, while oceanic poleward heat transport decreases. The non-CO2 boundary conditions cause more MHT toward the South Pole, mainly through an increase in the southward oceanic heat transport. The changes in paleogeography increase the heat transport via transient eddies at the northern mid-latitudes in the Eocene. The Eocene Hadley cells do not transport more heat poleward, but due to the warmer atmosphere, especially the northern cell, circulate more heat in the tropics, than today. The monsoon systems' poleward latent heat transport increases with rising CO2 concentrations, but this change is counterweighted by the globally smaller Eocene monsoon area. Our results show that the changes in the monsoon systems' latent heat transport is a robust feature of CO2 warming, which is in line with the currently observed precipitation increase of present day monsoon systems.

Estimates of global mean near-surface air temperature (global SAT) for the Cenozoic era rely largely on paleo-proxy data of deep-sea temperature (DST), with the assumption that changes in global SAT covary with changes in the global mean deep-sea temperature (global DST) and global mean sea-surface temperature (global SST). We tested the validity of this assumption by analyzing the relationship between global SST, SAT, and DST using 25 different model simulations from the Deep-Time Model Intercomparison Project simulating the early Eocene Climatic Optimum (EECO) with varying CO2 levels. Similar to the modern situation, we find limited spatial variability in DST, indicating that local DST estimates can be regarded as a first order representative of global DST. In line with previously assumed relationships, linear regression analysis indicates that both global DST and SAT respond stronger to changes in atmospheric CO2 than global SST by a similar factor. Consequently, this model-based analysis validates the assumption that changes in global DST can be used to estimate changes in global SAT during the early Cenozoic. Paleo-proxy estimates of global DST, SST, and SAT during EECO show the best fit with model simulations with a 1,680 ppm atmospheric CO2 level. This matches paleo-proxies of EECO atmospheric CO2, indicating a good fit between models and proxy-data.

Inconsistencies in the Eocene climates of East Asia have been revealed in both geological studies and simulations. Several earlier reconstructions showed an arid zonal band in mid-latitude China, but others showed a humid climate in the same region. Moreover, previous Eocene modeling studies have demonstrated that climate models can simulate both scenarios in China. Therefore, it is essential to investigate the cause of this model spread. We conducted a series of experiments using Norwegian Earth System Model 1-F and examined the impact of mountains in Southern China on the simulated Eocene climate. These mountains, including the Gangdese and Southeast Mountains, are located along the main path of water vapor transport to East Asia. Our results reveal that the Southeast Mountains play the dominant role in controlling the simulated precipitation in Eastern China during the Eocene. When the heights of the Southeast Mountains exceed similar to 2,000 m, an arid zonal band appears in mid-latitude China, whereas humid climates appear in Eastern China when the elevation of the Southeast Mountains is relatively low.

The early Eocene greenhouse climate maintained by high atmospheric CO2 concentrations serves as a testbed for future climate changes dominated by increasing CO2 forcing. In particular, the early Eocene Arctic region is important in the context of future CO2 driven climate warming in the northern polar region and associated shrinking Arctic sea ice. Here, we present early Eocene Arctic sea ice simulations carried out by six coupled climate models within the framework of the Deep-Time Model Intercomparison Project (DeepMIP). We find differences in sea ice responses to CO2 changes across the ensemble and compare the results with available proxy-based sea ice reconstructions from the Arctic Ocean. Most of the models simulate seasonal sea ice presence at high CO2 levels (≥ 840 ppmv = 3× pre-industrial (PI) level of 280 ppmv). However, the threshold when sea ice permanently disappears from the ocean varies considerably between the models (from <840 ppmv to >1680 ppmv). Based on a one-dimensional energy balance model analysis we find that the greenhouse effect likely caused by increased atmospheric water vapor concentration plays an important role in the inter-model spread in Arctic winter surface temperature changes in response to a CO2 rise from 1× to 3× the PI level. Furthermore, differences in simulated surface salinity in the Arctic Ocean play an important role in the control of local sea ice formation. These differences result from different implementations of river run-off between the models, but also from differences in the exchange of waters between a brackish Arctic and a more saline North Atlantic Ocean that are controlled by the width of the gateway between both basins. As there is no geological evidence for Arctic sea ice in the early Eocene, its presence in most of the simulations with 3× PI CO2 level indicates either a higher CO2 level and/or an overly weak polar sensitivity in these models.

The early Eocene (∼56–48 Myr ago) is characterized by high CO2 estimates (1,200–2,500 ppmv) and elevated global temperatures (∼10°C–16°C higher than modern). However, the response of the hydrological cycle during the early Eocene is poorly constrained, especially in regions with sparse data coverage (e.g., Africa). Here, we present a study of African hydroclimate during the early Eocene, as simulated by an ensemble of state-of-the-art climate models in the Deep-time Model Intercomparison Project (DeepMIP). A comparison between the DeepMIP pre-industrial simulations and modern observations suggests that model biases are model- and geographically dependent, however, these biases are reduced in the model ensemble mean. A comparison between the Eocene simulations and the pre-industrial suggests that there is no obvious wetting or drying trend as the CO2 increases. The results suggest that changes to the land sea mask (relative to modern) in the models may be responsible for the simulated increases in precipitation to the north of Eocene Africa. There is an increase in precipitation over equatorial and West Africa and associated drying over northern Africa as CO2 rises. There are also important dynamical changes, with evidence that anticyclonic low-level circulation is replaced by increased south-westerly flow at high CO2 levels. Lastly, a model-data comparison using newly compiled quantitative climate estimates from paleobotanical proxy data suggests a marginally better fit with the reconstructions at lower levels of CO2.

Here, we compare the ocean overturning circulation of the early Eocene (47–56 Ma) in eight coupled climate model simulations from the Deep-Time Model Intercomparison Project (DeepMIP) and investigate the causes of the observed inter-model spread. The most common global meridional overturning circulation (MOC) feature of these simulations is the anticlockwise bottom cell, fed by sinking in the Southern Ocean. In the North Pacific, one model (GFDL) displays strong deepwater formation and one model (CESM) shows weak deepwater formation, while in the Atlantic two models show signs of weak intermediate water formation (MIROC and NorESM). The location of the Southern Ocean deepwater formation sites varies among models and relates to small differences in model geometry of the Southern Ocean gateways. Globally, convection occurs in the basins with smallest local freshwater gain from the atmosphere. The global MOC is insensitive to atmospheric CO2 concentrations from 1× (i.e., 280 ppm) to 3× (840 ppm) pre-industrial levels. Only two models have simulations with higher CO2 (i.e., CESM and GFDL) and these show divergent responses, with a collapsed and active MOC, respectively, possibly due to differences in spin-up conditions. Combining the multiple model results with available proxy data on abyssal ocean circulation highlights that strong Southern Hemisphere-driven overturning is the most likely feature of the early Eocene. In the North Atlantic, unlike the present day, neither model results nor proxy data suggest deepwater formation in the open ocean during the early Eocene, while the evidence for deepwater formation in the North Pacific remains inconclusive.

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

We present results from an ensemble of eight climate models, each of which has carried out simulations of the early Eocene climate optimum (EECO, similar to 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 degrees 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 co

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

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

The mid-Cretaceous was one of the warmest intervals of the past 140 million years (Myr) driven by atmospheric CO2 levels around 1000 ppmv. In the near absence of proximal geological records from south of the Antarctic Circle, it remains disputed whether polar ice could exist under such environmental conditions. Here we present results from a unique sedimentary sequence recovered from the West Antarctic shelf. This by far southernmost Cretaceous record contains an intact ~3 m-long network of in-situ fossil roots. The roots are embedded in a mudstone matrix bearing diverse pollen and spores, indicative of a temperate lowland rainforest environment at a palaeolatitude of ~82°S during the Turonian–Santonian (93–83 Myr). A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric CO2 contents of 1120–1680 ppmv and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo in high-CO2 climate worlds.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years, 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.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years1,2,3,4,5, driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume6. 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.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years, 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.

The mid-Cretaceous was one of the warmest intervals of the past 140 million years (Myr) driven by atmospheric CO2 levels around 1000 ppmv. In the near absence of proximal geological records from south of the Antarctic Circle, it remains disputed whether polar ice could exist under such environmental conditions. Here we present results from a unique sedimentary sequence recovered from the West Antarctic shelf. This by far southernmost Cretaceous record contains an intact ~3 m-long network of in-situ fossil roots. The roots are embedded in a mudstone matrix bearing diverse pollen and spores, indicative of a temperate lowland rainforest environment at a palaeolatitude of ~82°S during the Turonian–Santonian (93–83 Myr). A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric CO2 contents of 1120–1680 ppmv and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo in high-CO2 climate worlds.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years1,2,3,4,5, driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume6. 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.

The mid-Cretaceous was one of the warmest intervals of the past 140 million years (Myr) driven by atmospheric CO2 levels around 1000 ppmv. In the near absence of proximal geological records from south of the Antarctic Circle, it remains disputed whether polar ice could exist under such environmental conditions. Here we present results from a unique sedimentary sequence recovered from the West Antarctic shelf. This by far southernmost Cretaceous record contains an intact ~3 m-long network of in-situ fossil roots. The roots are embedded in a mudstone matrix bearing diverse pollen and spores, indicative of a temperate lowland rainforest environment at a palaeolatitude of ~82°S during the Turonian–Santonian (92–83 Myr). A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric CO2 contents of 1120–1680 ppmv and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo in high-CO2 climate worlds.

The mid-Cretaceous was one of the warmest intervals of the past 140 million years (Myr) driven by atmospheric CO2 levels around 1000 ppmv. In the near absence of proximal geological records from south of the Antarctic Circle, it remains disputed whether polar ice could exist under such environmental conditions. Here we present results from a unique sedimentary sequence recovered from the West Antarctic shelf. This by far southernmost Cretaceous record contains an intact ~3 m-long network of in-situ fossil roots. The roots are embedded in a mudstone matrix bearing diverse pollen and spores, indicative of a temperate lowland rainforest environment at a palaeolatitude of ~82°S during the Turonian–Santonian (92–83 Myr). A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric CO2 contents of 1120–1680 ppmv and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo in high-CO2 climate worlds.

We investigate the impact of different CO2 levels and different subarctic gateway configurations on the surface temperatures during the latest Cretaceous using the Earth System Model COSMOS. The simulated temperatures are compared with the surface temperature reconstructions based on a recent compilation of the latest Cretaceous proxies. In our numerical experiments, the CO2 level ranges from 1 to 6 times the preindustrial (PI) CO2 level of 280 ppm. On a global scale, the most reasonable match between modeling and proxy data is obtained for the experiments with 3 to 5 × PI CO2 concentrations. However, the simulated low- (high-) latitude temperatures are too high (low) as compared to the proxy data. The moderate CO2 levels scenarios might be more realistic, if we take into account proxy data and the dead zone effect criterion. Furthermore, we test if the model-data discrepancies can be caused by too simplistic proxy-data interpretations. This is distinctly seen at high latitudes, where most proxies are biased toward summer temperatures. Additional sensitivity experiments with different ocean gateway configurations and constant CO2 level indicate only minor surface temperatures changes (<~1°C) on a global scale, with higher values (up to ~8°C) on a regional scale. These findings imply that modeled and reconstructed temperature gradients are to a large degree only qualitatively comparable, providing challenges for the interpretation of proxy data and/or model sensitivity. With respect to the latter, our results suggest that an assessment of greenhouse worlds is best constrained by temperatures in the midlatitudes.

We investigate the impact of different CO2 levels and different subarctic gateway configurations on the surface temperatures during the latest Cretaceous using the Earth System Model COSMOS. The simulated temperatures are compared with the surface temperature reconstructions based on a recent compilation of the latest Cretaceous proxies. In our numerical experiments, the CO2 level ranges from 1 to 6 times the preindustrial (PI) CO2 level of 280 ppm. On a global scale, the most reasonable match between modeling and proxy data is obtained for the experiments with 3 to 5 × PI CO2 concentrations. However, the simulated low- (high-) latitude temperatures are too high (low) as compared to the proxy data. The moderate CO2 levels scenarios might be more realistic, if we take into account proxy data and the dead zone effect criterion. Furthermore, we test if the model-data discrepancies can be caused by too simplistic proxy-data interpretations. This is distinctly seen at high latitudes, where most proxies are biased toward summer temperatures. Additional sensitivity experiments with different ocean gateway configurations and constant CO2 level indicate only minor surface temperatures changes (<~1°C) on a global scale, with higher values (up to ~8°C) on a regional scale. These findings imply that modeled and reconstructed temperature gradients are to a large degree only qualitatively comparable, providing challenges for the interpretation of proxy data and/or model sensitivity. With respect to the latter, our results suggest that an assessment of greenhouse worlds is best constrained by temperatures in the midlatitudes.

CO2-induced global warming will affect seasonal to decadal temperature patterns. Expected changes will be particularly strong in extratropical regions where temperatures will increase at faster rates than at lower latitudes. Despite that, it is still poorly constrained how precisely short-term climate dynamics will change in a generally warmer world, particularly in nearshore surface waters in the extratropics, i.e., the ecologically most productive regions of the ocean on which many human societies depend. Specifically, a detailed knowledge of the relationship between pCO2 and seasonal SST is crucial to understand interactions between the ocean and the atmosphere. In the present investigation, we have studied for the first time how rising atmospheric pCO2 levels forced surface temperature changes in Central Europe (paleolatitude ~45 °N) during the mid-Oligocene (fromca. 31 to 25Ma), a time interval of Earth history during which global conditions were comparable to those predicted for the next few centuries. For this purpose, we computed numerical climate models for the Oligocene (winter, summer, annual average) assuming an atmospheric carbon dioxide rise from 400 to 560 ppm (current level to two times pre-industrial levels, PAL) and from 400 to 840 ppm (= three times PAL), respectively. These models were compared to seasonally resolved sea surface temperatures (SST) reconstructed from δ18O values of fossil bivalve shells (Glycymeris planicostalis, G. obovata, Palliolum pictum, Arctica islandica and Isognomon maxillata sandbergeri) and shark teeth (Carcharias cuspidata, C. acutissimaand Physogaleus latus) collected fromthe shallow water deposits of the Mainz and Kassel Basins (Germany). Multi-taxon oxygen isotope-based reconstructions suggest a gradual rise of temperatures in surface waters (upper 30 to 40m), on average, by asmuch as 4 °C during the Rupelian stage followed by a 4 °C cooling during the Chattian stage. Seasonal temperature amplitudes increased by ca.