Description: Two interglacial epochs are included in the suite of Paleoclimate Modeling Intercomparison Project (PMIP4) simulations in the Coupled Model Intercomparison Project (CMIP6). The experimental protocols for simulations of the mid-Holocene (midHolocene, 6000 years before present) and the Last Interglacial (lig127k, 127 000 years before present) are described here. These equilibrium simulations are designed to examine the impact of changes in orbital forcing at times when atmospheric greenhouse gas levels were similar to those of the preindustrial period and the continental configurations were almost identical to modern ones. These simulations test our understanding of the interplay between radiative forcing and atmospheric circulation, and the connections among large-scale and regional climate changes giving rise to phenomena such as land–sea contrast and high-latitude amplification in temperature changes, and responses of the monsoons, as compared to today. They also provide an opportunity, through carefully designed additional sensitivity experiments, to quantify the strength of atmosphere, ocean, cryosphere, and land-surface feedbacks. Sensitivity experiments are proposed to investigate the role of freshwater forcing in triggering abrupt climate changes within interglacial epochs. These feedback experiments naturally lead to a focus on climate evolution during interglacial periods, which will be examined through transient experiments. Analyses of the sensitivity simulations will also focus on interactions between extratropical and tropical circulation, and the relationship between changes in mean climate state and climate variability on annual to multi-decadal timescales. The comparative abundance of paleoenvironmental data and of quantitative climate reconstructions for the Holocene and Last Interglacial make these two epochs ideal candidates for systematic evaluation of model performance, and such comparisons will shed new light on the importance of external feedbacks (e.g., vegetation, dust) and the ability of state-of-the-art models to simulate climate changes realistically.

Description: The Last Glacial Maximum (LGM, 21 000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: Little is known about the climate evolution on the Kamchatka Peninsula during the last glacial–interglacial transition as existing climate records do not reach beyond 12 ka BP. In this study, a summer temperature record for the past 20 kyr is presented. Branched glycerol dialkyl glycerol tetraethers, terrigenous biomarkers suitable for continental air temperature reconstructions, were analyzed in a sediment core from the western continental margin off Kamchatka in the marginal northwest Pacific (NW Pacific). The record suggests that summer temperatures on Kamchatka during the Last Glacial Maximum (LGM) equaled modern temperatures. We suggest that strong southerly winds associated with a pronounced North Pacific High pressure system over the subarctic NW Pacific accounted for the warm conditions. A comparison with an Earth system model reveals discrepancies between model and proxy-based reconstructions for the LGM temperature and atmospheric circulation in the NW Pacific realm. The deglacial temperature development is characterized by abrupt millennial-scale temperature oscillations. The Bølling–Allerød warm phase and the Younger Dryas cold spell are pronounced events, suggesting a connection to North Atlantic climate variability.

Description: We present here the first results, for the preindustrial and mid-Holocene climatological periods, of the newly developed isotope-enhanced version of the fully coupled Earth system model MPI-ESM, called hereafter MPI-ESM-wiso. The water stable isotopes H162O, H182O and HDO have been implemented into all components of the coupled model setup. The mid-Holocene provides the opportunity to evaluate the model response to changes in the seasonal and latitudinal distribution of insolation induced by different orbital forcing conditions. The results of our equilibrium simulations allow us to evaluate the performance of the isotopic model in simulating the spatial and temporal variations of water isotopes in the different compartments of the hydrological system for warm climates. For the preindustrial climate, MPI-ESM-wiso reproduces very well the observed spatial distribution of the isotopic content in precipitation linked to the spatial variations in temperature and precipitation rate. We also find a good model–data agreement with the observed distribution of isotopic composition in surface seawater but a bias with the presence of surface seawater that is too 18O-depleted in the Arctic Ocean. All these results are improved compared to the previous model version ECHAM5/MPIOM. The spatial relationships of water isotopic composition with temperature, precipitation rate and salinity are consistent with observational data. For the preindustrial climate, the interannual relationships of water isotopes with temperature and salinity are globally lower than the spatial ones, consistent with previous studies. Simulated results under mid-Holocene conditions are in fair agreement with the isotopic measurements from ice cores and continental speleothems. MPI-ESM-wiso simulates a decrease in the isotopic composition of precipitation from North Africa to the Tibetan Plateau via India due to the enhanced monsoons during the mid-Holocene. Over Greenland, our simulation indicates a higher isotopic composition of precipitation linked to higher summer temperature and a reduction in sea ice, shown by positive isotope–temperature gradient. For the Antarctic continent, the model simulates lower isotopic values over the East Antarctic plateau, linked to the lower temperatures during the mid-Holocene period, while similar or higher isotopic values are modeled over the rest of the continent. While variations of isotopic contents in precipitation over West Antarctica between mid-Holocene and preindustrial periods are partly controlled by changes in temperature, the transport of relatively 18O-rich water vapor near the coast to the western ice core sites could play a role in the final isotopic composition. So, more caution has to be taken about the reconstruction of past temperature variations during warm periods over this area. The coupling of such a model with an ice sheet model or the use of a zoomed grid centered on this region could help to better describe the role of the water vapor transport and sea ice around West Antarctica. The reconstruction of past salinity through isotopic content in sea surface waters can be complicated for regions with strong ocean dynamics, variations in sea ice regimes or significant changes in freshwater budget, giving an extremely variable relationship between the isotopic content and salinity of ocean surface waters over small spatial scales. These complicating factors demonstrate the complexity of interpreting water isotopes as past climate signals of warm periods like the mid-Holocene. A systematic isotope model intercomparison study for further insights on the model dependency of these results would be beneficial.

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 propose a surface melt scheme for glaciated land surfaces, which only requires monthly mean short-wave radiation and temperature as inputs, yet implicitly accounts for the diurnal cycle of short-wave radiation. The scheme is deduced from the energy balance of a daily melt period, which is defined by a minimum solar elevation angle. The scheme yields a better spatial representation of melting than common empirical schemes when applied to the Greenland Ice Sheet, using a 1948–2016 regional climate and snowpack simulation as a reference. The scheme is physically constrained and can be adapted to other regions or time periods.

Description: The Last Glacial Maximum (LGM, ~ 21,000 years ago) has been a major focus for evaluating how well state-of-the-art climate models simulate climate changes as large as those expected in the future using paleoclimate reconstructions. A new generation of climate models have been used to generate LGM simulations as part of the Palaeoclimate Modelling Intercomparison Project (PMIP) contribution to the Coupled Model Intercomparison Project (CMIP). Here we provide a preliminary analysis and evaluation of the results of these LGM experiments (PMIP4-CMIP6) and compare them with the previous generation of simulations (PMIP3-CMIP5). We show that the PMIP4-CMIP6 are globally less cold and less dry than the PMIP3-CMIP5 simulations, most probably because of the use of a more realistic specification of the northern hemisphere ice sheets in the latest simulations although changes in model configuration may also contribute to this. There are important differences in both atmospheric and ocean circulation between the two sets of experiments, with the northern and southern jet streams being more poleward and the changes in the Atlantic Meridional Overturning Circulation being less pronounced in the PMIP4-CMIP6 simulations than in the PMIP3-CMIP5 simulations. Changes in simulated precipitation patterns are influenced by both temperature and circulation changes. Differences in simulated climate between individual models remain large so, although there are differences in the average behaviour across the two ensembles, the new simulation results are not fundamentally different from the PMIP3-CMIP5 results. Evaluation of large-scale climate features, such as land-sea contrast and polar amplification, confirms that the models capture these well and within the uncertainty of the palaeoclimate reconstructions. Nevertheless, regional climate changes are less well simulated: the models underestimate extratropical cooling, particularly in winter, and precipitation changes. The spatial patterns of increased precipitation associated with changes in the jet streams are also poorly captured. However, changes in the tropics are more realistic, particularly the changes in tropical temperatures over the oceans. Although these results are preliminary in nature, because of the limited number of LGM simulations currently available, they nevertheless point to the utility of using paleoclimate simulations to understand the mechanisms of climate change and evaluate model performance.

Description: Holocene sea surface temperature trends and variability are underestimated in models compared to paleoclimate data. The idea is presented that the local trends and variability are related, which is elaborated in a conceptual framework of the stochastic climate model. The relation is a consequence of the fluctuation–dissipation theorem, connecting the linear response of a system to its statistical fluctuations. Consequently, the spectrum can be used to estimate the timescale-dependent climate response. The non-normality in the propagation operator introduces enhanced long-term variability related to nonequilibrium and/or Earth system sensitivity. The simple model can guide us to analyze comprehensive models' behavior.

Description: The Greenland Ice Sheet (GrIS) mass loss has been accelerating at a rate of about 20 ± 10 Gt/yr2 since the end of the 1990's, with around 60 % of this mass loss directly attributed to enhanced surface meltwater runoff. However, in the climate and glaciology communities, different approaches exist on how to model the different surface mass balance (SMB) components using: (1) complex physically-based climate models which are computationally expensive; (2) intermediate complexity energy balance models; (3) simple and fast positive degree day models which base their inferences on statistical principles and are computationally highly efficient. Additionally, many of these models compute the SMB components based on different spatial and temporal resolutions, with different forcing fields as well as different ice sheet topographies and extents, making inter-comparison difficult. In the GrIS SMB model intercomparison project (GrSMBMIP) we address these issues by forcing each model with the same data (i.e., the ERA-Interim reanalysis) except for two global models for which this forcing is limited to the oceanic conditions, and at the same time by interpolating all modelled results onto a common ice sheet mask at 1 km horizontal resolution for the common period 1980–2012. The SMB outputs from 13 models are then compared over the GrIS to (1) SMB estimates using a combination of gravimetric remote sensing data from GRACE and measured ice discharge, (2) ice cores, snow pits, in-situ SMB observations, and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Our results reveal that the mean GrIS SMB of all 13 models has been positive between 1980 and 2012 with an average of 340 ± Gt/yr, but has decreased at an average rate of −7.3 Gt/yr2 (with a significance of 96 %), mainly driven by an increase of 8.0 Gt/yr2 (with a significance of 98 %) in meltwater runoff. Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting the need for accurate representation of the GrIS ablation zone extent and processes driving the surface melt. In addition, a higher density of in-situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 mWE/yr due to large discrepancies in modelled snowfall accumulation. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of same order than RCMs with observations and remain then useful tools for long-term simulations. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present day SMB relative to observations, suggesting that biases are not systematic among models.

Description: During the Last Interglacial (LIG, ∼130–115kiloyears (kyr) before present (BP)), the northern high latitudes were characterized by higher temperatures than those of the late Holocene and a lower Greenland Ice Sheet (GIS). However, the impact of a reduced GIS on the global climate has not yet been well constrained. In this study, we quantify the contribution of the GIS to LIG warmth by performing various sensitivity studies based on equilibrium simulations, employing the Community Earth System Models (COSMOS), with a focus on height and extent of the GIS. We present the first study on the effects of a reduction in the GIS on the surface temperature (TS) on a global scale and separate the contribution of astronomical forcing and changes in GIS to LIG warmth. The strong Northern Hemisphere summer warming of approximately 2°C (with respect to pre-industrial) is mainly caused by increased summer insolation. Reducing the height by  ∼ 1300m and the extent of the GIS does not have a strong influence during summer, leading to an additional global warming of only +0.24°C compared to the purely insolation-driven LIG. The effect of a reduction in the GIS is, however, strongest during local winter, with up to +5°C regional warming and with an increase in global average temperature of +0.48°C. In order to evaluate the performance of our LIG simulations, we additionally compare the simulated TS anomalies with marine and terrestrial proxy-based LIG temperature anomalies derived from three different proxy data compilations. Our model results are in good agreement with proxy records with respect to the warming pattern but underestimate the magnitude of temperature change when compared to reconstructions, suggesting a potential misinterpretation of the proxy records or deficits in our model. However, we are able to partly reduce the mismatch between model and data by additionally taking into account the potential seasonal bias of the proxy record and/or the uncertainties in the dating of the proxy records for the LIG thermal maximum. The seasonal bias and the uncertainty of the timing are estimated from new transient model simulations covering the whole LIG. The model–data comparison improves for proxies that represent annual mean temperatures when the GIS is reduced and when we take the local thermal maximum during the LIG (130–120kyrBP) into account. For proxy data that represent summer temperatures, changes in the GIS are of minor importance for sea surface temperatures. However, the annual mean and summer temperature change over Greenland in the reduced GIS simulations seems to be overestimated as compared to the local ice core data, which could be related to the interpretation of the recorder system and/or the assumptions of GIS reduction. Thus, the question regarding the real size of the GIS during the LIG has yet to be answered.