ALBERT

Description: The ice–substrate interface is an important boundary condition for ice sheet modelling. The substrate affects the ice sheet by allowing sliding through sediment deformation and accommodating the storage and drainage of subglacial water. We present three datasets on a 1 : 5 000 000 scale with different geological parameters for the region that was covered by the ice sheets in North America, including Greenland and Iceland. The first dataset includes the distribution of surficial sediments, which is separated into continuous, discontinuous and predominantly rock categories. The second dataset includes sediment grain size properties, which is divided into three classes: clay, silt and sand, based on the dominant grain size of the fine fraction of the glacial sediments. The third dataset is the generalized bedrock geology. We demonstrate the utility of these datasets for governing ice sheet dynamics by using an ice sheet model with a simulation that extends through the last glacial cycle. In order to demonstrate the importance of the basal boundary conditions for ice sheet modelling, we changed the shear friction angle to account for a weaker substrate and found changes up to 40 % in ice thickness compared to a reference run. Although incorporation of the ice–bed boundary remains model dependent, our dataset provides an observational baseline for improving a critical weakness in current ice sheet modelling (https://doi.org/10.1594/PANGAEA.895889, Gowan et al., 2018b).

Description: The ice-substrate interface is an important boundary condition for ice sheet modelling. The substrate affects the ice sheet by allowing sliding through sediment deformation and accommodating the storage and drainage of subglacial water. We present three datasets with different geological parameters for the region that was covered by the ice sheets in North America, including Greenland and Iceland. The first dataset includes the distribution surficial sediments, which is separated into continuous, discontinuous and predominantly rock categories. The second dataset includes sediment grain size properties, which is divided into three classes: clay, silt and sand, based on the dominant grain size of the glacial sediments. The third dataset is the generalized bedrock geology. We demonstrate the utility of these datasets for governing ice sheet dynamics by using an ice sheet model with a simulation that extends through the last glacial cycle. Changes in ice thickness by up to 40 % relative to a reference simulation happened when the shear friction angle was reduced to account for a weaker substrate. These datasets provide a basis to improve the basal boundary conditions in ice sheet models. Gowan, E. J., Niu, L., Knorr, G., and Lohmann, G., 2018. Geology datasets of North America for use with ice sheet models, link to datafiles. PANGAEA, https://doi.pangaea.de/10.1594/PANGAEA.895889

Description: We present the Alfred Wegener Institute's contribution to the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) wherein we employ the Community Earth System Models (COSMOS) that include a dynamic vegetation scheme. This work builds on our contribution to Phase 1 of the Pliocene Model Intercomparison Project (PlioMIP1) wherein we employed the same model without dynamic vegetation. Our input to the PlioMIP2 special issue of Climate of the Past is twofold. In an accompanying paper we compare results derived with COSMOS in the framework of PlioMIP2 and PlioMIP1. With this paper we present details of our contribution with COSMOS to PlioMIP2. We provide a description of the model and of methods employed to transfer reconstructed mid-Pliocene geography, as provided by the Pliocene Reconstruction and Synoptic Mapping Initiative Phase 4 (PRISM4), to model boundary conditions. We describe the spin-up procedure for creating the COSMOS PlioMIP2 simulation ensemble and present large-scale climate patterns of the COSMOS PlioMIP2 mid-Pliocene core simulation. Furthermore, we quantify the contribution of individual components of PRISM4 boundary conditions to characteristics of simulated mid-Pliocene climate and discuss implications for anthropogenic warming. When exposed to PRISM4 boundary conditions, COSMOS provides insight into a mid-Pliocene climate that is characterised by increased rainfall (+0.17 mm d−1) and elevated surface temperature (+3.37 ∘C) in comparison to the pre-industrial (PI). About two-thirds of the mid-Pliocene core temperature anomaly can be directly attributed to carbon dioxide that is elevated with respect to PI. The contribution of topography and ice sheets to mid-Pliocene warmth is much smaller in contrast – about one-quarter and one-eighth, respectively, and nonlinearities are negligible. The simulated mid-Pliocene climate comprises pronounced polar amplification, a reduced meridional temperature gradient, a northwards-shifted tropical rain belt, an Arctic Ocean that is nearly free of sea ice during boreal summer, and muted seasonality at Northern Hemisphere high latitudes. Simulated mid-Pliocene precipitation patterns are defined by both carbon dioxide and PRISM4 paleogeography. Our COSMOS simulations confirm long-standing characteristics of the mid-Pliocene Earth system, among these increased meridional volume transport in the Atlantic Ocean, an extended and intensified equatorial warm pool, and pronounced poleward expansion of vegetation cover. By means of a comparison of our results to a reconstruction of the sea surface temperature (SST) of the mid-Pliocene we find that COSMOS reproduces reconstructed SST best if exposed to a carbon dioxide concentration of 400 ppmv. In the Atlantic to Arctic Ocean the simulated mid-Pliocene core climate state is too cold in comparison to the SST reconstruction. The discord can be mitigated to some extent by increasing carbon dioxide that causes increased mismatch between the model and reconstruction in other regions.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.

Description: We introduce a new hypothesis concerning the role of internal climate dynamics in the non-linear transitions from interglacial to glacial (IG-G) state since the Mid Pleistocene Transition (MPT). These transitions encompass large and abrupt changes in atmospheric CO2, ice volume, and temperature that we suggest involve critical interactions between insolation and high amplitude oscillations in ocean/atmosphere circulation patterns. Specifically, we highlight the large amplitude of millennial-scale climate oscillations across the transition from Marine Isotope Stage (MIS) 5 to 4, which we argue led to amplified cooling of the deep ocean and we demonstrate that analogous episodes of extreme cooling systematically preceded glacial periods of the last 800 kyr. We suggest that such cooling necessitates a reconfiguration of the deep ocean to avoid a density paradox between northern and southern-sourced deep waters (SSW), which could be accomplished by increasing the relative volume and or salinity of SSW, thus providing the necessary storage capacity for the subsequent (delayed) and relatively abrupt drawdown of CO2. We therefore explain the transient decoupling of Antarctic temperature from CO2 across MIS 5/4 as a direct consequence of millennial activity at that time. We further show that similar climatic decoupling typically occurred during times of low obliquity and was a ubiquitous feature of IG-G transitions over the past 800 kyr, producing the appearance of bimodality in records of CO2, benthic δ18O and others. Finally we argue that the apparent lack of bimodality in the pre-MPT record of benthic δ18O implies that the dynamics associated with IG-G transitions changed across the MPT

Description: We present the Alfred Wegener Institute's contribution to the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) wherein we employ the Community Earth System Models (COSMOS) that include a dynamic vegetation scheme. This work builds on our contribution to Phase 1 of the Pliocene Model Intercomparison Project (PlioMIP1) wherein we employed the same model without dynamic vegetation. Our input to the PlioMIP2 special issue of Climate of the Past is twofold. In an accompanying paper we compare results derived with COSMOS in the framework of PlioMIP2 and PlioMIP1. With this paper we present details of our contribution with COSMOS to PlioMIP2. We provide a description of the model and of methods employed to transfer reconstructed mid-Pliocene geography, as provided by the Pliocene Reconstruction and Synoptic Mapping Initiative Phase 4 (PRISM4), to model boundary conditions. We describe the spin-up procedure for creating the COSMOS PlioMIP2 simulation ensemble and present large-scale climate patterns of the COSMOS PlioMIP2 mid-Pliocene core simulation. Furthermore, we quantify the contribution of individual components of PRISM4 boundary conditions to characteristics of simulated mid-Pliocene climate and discuss implications for anthropogenic warming. When exposed to PRISM4 boundary conditions, COSMOS provides insight into a mid-Pliocene climate that is characterised by increased rainfall (+0.17 mm d−1) and elevated surface temperature (+3.37 ∘C) in comparison to the pre-industrial (PI). About two-thirds of the mid-Pliocene core temperature anomaly can be directly attributed to carbon dioxide that is elevated with respect to PI. The contribution of topography and ice sheets to mid-Pliocene warmth is much smaller in contrast – about one-quarter and one-eighth, respectively, and nonlinearities are negligible. The simulated mid-Pliocene climate comprises pronounced polar amplification, a reduced meridional temperature gradient, a northwards-shifted tropical rain belt, an Arctic Ocean that is nearly free of sea ice during boreal summer, and muted seasonality at Northern Hemisphere high latitudes. Simulated mid-Pliocene precipitation patterns are defined by both carbon dioxide and PRISM4 paleogeography. Our COSMOS simulations confirm long-standing characteristics of the mid-Pliocene Earth system, among these increased meridional volume transport in the Atlantic Ocean, an extended and intensified equatorial warm pool, and pronounced poleward expansion of vegetation cover. By means of a comparison of our results to a reconstruction of the sea surface temperature (SST) of the mid-Pliocene we find that COSMOS reproduces reconstructed SST best if exposed to a carbon dioxide concentration of 400 ppmv. In the Atlantic to Arctic Ocean the simulated mid-Pliocene core climate state is too cold in comparison to the SST reconstruction. The discord can be mitigated to some extent by increasing carbon dioxide that causes increased mismatch between the model and reconstruction in other regions.

Description: 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.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land ice sheet growth and sea level fall) the multimillennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future, we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.