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

Description: There is an increasing need to understand the pre-Quaternary warm climates, how climate–vegetation interactions functioned in the past, and how we can use this information to understand the present. Here we report vegetation modelling results for the Late Miocene (11–7 Ma) to study the mechanisms of vegetation dynamics and the role of different forcing factors that influence the spatial patterns of vegetation coverage. One of the key uncertainties is the atmospheric concentration of CO2 during past climates. Estimates for the last 20 million years range from 280 to 500 ppm. We simulated Late Miocene vegetation using two plausible CO2 concentrations, 280 ppm CO2 and 450 ppm CO2, with a dynamic global vegetation model (LPJ-GUESS) driven by climate input from a coupled AOGCM (Atmosphere-Ocean General Circulation Model). The simulated vegetation was compared to existing plant fossil data for the whole Northern Hemisphere. For the comparison we developed a novel approach that uses information of the relative dominance of different plant functional types (PFTs) in the palaeobotanical data to provide a quantitative estimate of the agreement between the simulated and reconstructed vegetation. Based on this quantitative assessment we find that pre-industrial CO2 levels are largely consistent with the presence of seasonal temperate forests in Europe (suggested by fossil data) and open vegetation in North America (suggested by multiple lines of evidence). This suggests that during the Late Miocene the CO2 levels have been relatively low, or that other factors that are not included in the models maintained the seasonal temperate forests and open vegetation.

Description: In order to account for coupled climate–soil processes, we have developed a soil scheme which is asynchronously coupled to a comprehensive climate model with dynamic vegetation. This scheme considers vegetation as the primary control of changes in physical soil characteristics. We test the scheme for a warmer (mid-Holocene) and colder (Last Glacial Maximum) climate relative to the preindustrial climate. We find that the computed changes in physical soil characteristics lead to significant amplification of global climate anomalies, representing a positive feedback. The inclusion of the soil feedback yields an extra surface warming of 0.24 °C for the mid-Holocene and an additional global cooling of 1.07 °C for the Last Glacial Maximum. Transition zones such as desert–savannah and taiga–tundra exhibit a pronounced response in the model version with dynamic soil properties. Energy balance model analyses reveal that our soil scheme amplifies the temperature anomalies in the mid-to-high northern latitudes via changes in the planetary albedo and the effective longwave emissivity. As a result of the modified soil treatment and the positive feedback to climate, part of the underestimated mid-Holocene temperature response to orbital forcing can be reconciled in the model.

Description: The last deglaciation is one of the best constrained global-scale climate changes documented by climate archives. Nevertheless, understanding of the underlying dynamics is still limited, especially with respect to abrupt climate shifts and associated changes in the Atlantic meridional overturning circulation (AMOC) during glacial and deglacial periods. A fundamental issue is how to obtain an appropriate climate state at the Last Glacial Maximum (LGM, 21 000 yr before present, 21 ka BP) that can be used as an initial condition for deglaciation. With the aid of a comprehensive climate model, we found that initial ocean states play an important role on the equilibrium timescale of the simulated glacial ocean. Independent of the initialization, the climatological surface characteristics are similar and quasi-stationary, even when trends in the deep ocean are still significant, which provides an explanation for the large spread of simulated LGM ocean states among the Paleoclimate Modeling Intercomparison Project phase 2 (PMIP2) models. Accordingly, we emphasize that caution must be taken when alleged quasi-stationary states, inferred on the basis of surface properties, are used as a reference for both model inter-comparison and data model comparison. The simulated ocean state with the most realistic AMOC is characterized by a pronounced vertical stratification, in line with reconstructions. Hosing experiments further suggest that the response of the glacial ocean is dependent on the ocean background state, i.e. only the state with robust stratification shows an overshoot behavior in the North Atlantic. We propose that the salinity stratification represents a key control on the AMOC pattern and its transient response to perturbations. Furthermore, additional experiments suggest that the stratified deep ocean formed prior to the LGM during a time of minimum obliquity (~ 27 ka BP). This indicates that changes in the glacial deep ocean already occur before the last deglaciation. In combination, these findings represent a new paradigm for the LGM and the last deglaciation, which challenges the conventional evaluation of glacial and deglacial AMOC changes based on an ocean state derived from 21 ka BP boundary conditions.