Articles | Volume 17, issue 1
Research article
11 Jan 2021
Research article |  | 11 Jan 2021

A multi-model CMIP6-PMIP4 study of Arctic sea ice at 127 ka: sea ice data compilation and model differences

Masa Kageyama, Louise C. Sime, Marie Sicard, Maria-Vittoria Guarino, Anne de Vernal, Ruediger Stein, David Schroeder, Irene Malmierca-Vallet, Ayako Abe-Ouchi, Cecilia Bitz, Pascale Braconnot, Esther C. Brady, Jian Cao, Matthew A. Chamberlain, Danny Feltham, Chuncheng Guo, Allegra N. LeGrande, Gerrit Lohmann, Katrin J. Meissner, Laurie Menviel, Polina Morozova, Kerim H. Nisancioglu, Bette L. Otto-Bliesner, Ryouta O'ishi, Silvana Ramos Buarque, David Salas y Melia, Sam Sherriff-Tadano, Julienne Stroeve, Xiaoxu Shi, Bo Sun, Robert A. Tomas, Evgeny Volodin, Nicholas K. H. Yeung, Qiong Zhang, Zhongshi Zhang, Weipeng Zheng, and Tilo Ziehn

The Last Interglacial period (LIG) is a period with increased summer insolation at high northern latitudes, which results in strong changes in the terrestrial and marine cryosphere. Understanding the mechanisms for this response via climate modelling and comparing the models' representation of climate reconstructions is one of the objectives set up by the Paleoclimate Modelling Intercomparison Project for its contribution to the sixth phase of the Coupled Model Intercomparison Project. Here we analyse the results from 16 climate models in terms of Arctic sea ice. The multi-model mean reduction in minimum sea ice area from the pre industrial period (PI) to the LIG reaches 50 % (multi-model mean LIG area is 3.20×106 km2, compared to 6.46×106 km2 for the PI). On the other hand, there is little change for the maximum sea ice area (which is 15–16×106 km2 for both the PI and the LIG. To evaluate the model results we synthesise LIG sea ice data from marine cores collected in the Arctic Ocean, Nordic Seas and northern North Atlantic. The reconstructions for the northern North Atlantic show year-round ice-free conditions, and most models yield results in agreement with these reconstructions. Model–data disagreement appear for the sites in the Nordic Seas close to Greenland and at the edge of the Arctic Ocean. The northernmost site with good chronology, for which a sea ice concentration larger than 75 % is reconstructed even in summer, discriminates those models which simulate too little sea ice. However, the remaining models appear to simulate too much sea ice over the two sites south of the northernmost one, for which the reconstructed sea ice cover is seasonal. Hence models either underestimate or overestimate sea ice cover for the LIG, and their bias does not appear to be related to their bias for the pre-industrial period. Drivers for the inter-model differences are different phasing of the up and down short-wave anomalies over the Arctic Ocean, which are associated with differences in model albedo; possible cloud property differences, in terms of optical depth; and LIG ocean circulation changes which occur for some, but not all, LIG simulations. Finally, we note that inter-comparisons between the LIG simulations and simulations for future climate with moderate (1 % yr−1) CO2 increase show a relationship between LIG sea ice and sea ice simulated under CO2 increase around the years of doubling CO2. The LIG may therefore yield insight into likely 21st century Arctic sea ice changes using these LIG simulations.

1 Introduction

The Last Interglacial period (LIG) was the last time global temperature was substantially higher than the pre-industrial period (PI) at high northern latitudes. It is important in helping us understand warm-climate sea ice and climate dynamics (Otto-Bliesner et al.2013, 2017; Capron et al.2017; Fischer et al.2018). Stronger LIG spring and summertime insolation contributed to this warmth, as well as feedbacks amplifying the initial insolation signal, in particular feedbacks related to the marine and land cryosphere. Previous climate model simulations of the LIG, forced by appropriate greenhouse gas (GHG) and orbital changes, have failed to capture the observed high temperatures at higher latitudes (Malmierca-Vallet et al.2018; Masson-Delmotte et al.2011; Otto-Bliesner et al.2013; Lunt et al.2013). Models used during the previous Coupled Model Intercomparison Project 5 (CMIP5) disagree on the magnitude of Arctic sea ice retreat during the LIG: the diversity of sea ice behaviour across models was linked to the spread in simulated surface temperatures and in the magnitude of the polar amplification (Otto-Bliesner et al.2013; Lunt et al.2013; IPCC2013). However it was difficult to compare some of the LIG simulations because they were not all run using identical protocol. These studies thus highlighted the need of a systematic approach to study the role of Arctic sea ice changes during the LIG.

Coupled Model Intercomparison Projects (CMIPs) coordinate and design climate model protocols for past, present and future climates and have become an indispensable tool to facilitate our understanding of climate change (IPCC2007, 2013; Eyring et al.2016). The Paleoclimate Model Intercomparison Project 4 (PMIP4) is one of the individual Model Intercomparison Projects that is taking part in CMIP6 (Kageyama et al.2018). Within this framework, a common experimental protocol for LIG climate simulation was developed by Otto-Bliesner et al. (2017). CMIP models differ among each other in their physical formulation, numerical discretisation and code implementation. However, this CMIP6-PMIP4 LIG standard protocol facilitates model inter-comparison work.

Alongside a previous lack of a common experimental protocol, our ability to evaluate CMIP models has previously been hindered by difficulties in determining LIG sea ice extent from marine core evidence (e.g. Otto-Bliesner et al.2013; Sime et al.2013; Malmierca-Vallet et al.2018; Stein et al.2017). Planktonic foraminifer assemblages that include a subpolar component suggest reduced sea ice in the Arctic Ocean (Nørgaard-Pedersen et al.2007; Adler et al.2009). Microfauna found in LIG marine sediments recovered from the Beaufort Sea Shelf, an area characterised by ice-free conditions during summers today, also support ice-free conditions during those times; this indicates that more saline Atlantic water was present on the Beaufort Shelf, suggesting reduced perennial Arctic sea ice during some part of the LIG (Brigham-Grette and Hopkins1995). On the other hand, a reconstruction of LIG Arctic sea ice changes based on sea ice biomarker proxies (see below for details) suggests that the central part of the LIG Arctic Ocean remained covered by ice throughout the year, while a significant reduction of LIG sea ice occurred across the Barents Sea continental margin (Stein et al.2017). On the modelling side, no previous coupled climate model has simulated an ice-free Arctic during the LIG (Otto-Bliesner et al.2006; Lunt et al.2013; Otto-Bliesner et al.2013; Stein et al.2017).

Here we address the question of LIG Arctic sea ice by providing a new marine core synthesis. Additionally, the CMIP6-PMIP4 LIG experimental protocol developed by Otto-Bliesner et al. (2017) provides the systematic framework to enable us to examine the question of the simulation of LIG Arctic sea ice using a multi-model approach. This is important given the current level of interest in the ability of climate models to accurately represent key Arctic climate processes during warm periods, including sea ice formation and melting. We compare the LIG Arctic sea ice simulated by each model against our new data synthesis and investigate why different models show different Arctic sea ice behaviour.

2 Materials and methods

2.1 Current Arctic sea ice

Our main objective is to investigate LIG sea ice. However, a quick assessment of the sea ice simulated in the reference state, i.e. the pre-industrial control experiment (referred to as piControl in the CMIP6 terminology, and PI in this paper) was necessary. In the absence of extensive sea ice data for the PI, we used data for a recent period before the current sea ice cover significant decrease. We use the NOAA Optimum Interpolation version 2 data (Reynolds et al.2002) for the period 1982 to 2001. The sea ice data in this dataset are obtained from different satellite and in situ observations. We have used the monthly time series at a resolution of 1. This dataset is termed “NOAA_OI_v2” in the rest of this paper.

2.2 Marine records of LIG Arctic sea ice

We focus here on records of sea ice from marine cores. Table 1 provides a summary of LIG sea ice information and data obtained from marine sediment cores collected in the Arctic Ocean, Nordic Seas and northern North Atlantic. South of 78 N, the records show ice-free conditions. Most of these sea ice records are derived from quantitative estimates of sea surface parameters based on dinoflagellate cysts (dinocysts). North of 78 N the sea-ice-related records are rare and different types of indicators were used. In addition to dinocysts, the records are based on biomarkers linked to phototrophic productivity in sea ice and on foraminifers and ostracods that both provide indication on water properties and indirectly on sea ice (de Vernal et al.2013b). Between 78 and 87 N, the faunal data have been interpreted as indicating seasonal sea ice cover conditions during the LIG.

Cronin et al. (2010)Stein et al. (2017)Cronin et al. (2010)Stein et al. (2017)Nørgaard-Pedersen et al. (2007)Kremer et al.2018bMatthiessen et al.2001Matthiessen and Knies (2001)Stein et al. (2017)Stein et al. (2017)Adler et al. (2009)Cronin et al. (2010)Kremer et al.2018aCronin et al. (2010)20112013Zhuravleva et al.2017201120082008Eynaud (1999)Eynaud et al.2004Hillaire-Marcel et al.200120082011Hodell et al. (2009)Hillaire-Marcel et al. (2011)Penaud et al. (2008)Eynaud et al.2000

Table 1Marine core records of Arctic sea ice from MIS5e. The references indicated for the dinocyst reconstructions are those for the initial core, the reconstruction itself follows de Vernal et al. (2013a, b, 2020) (see the main text for details). SIC stands for sea ice concentration. Question marks indicate that seasonal duration or annual mean SIC are uncertain and not quantitatively estimated. More general qualitative statements are still possible and are given in the “Qualitative sea ice state” column.

Download XLSX

Among sea ice cover indicators, dinocyst assemblages have been used as quantitative proxy based on the application of the modern analogue technique applied to a standardised reference modern data base developed from surface sediment samples collected at middle to high latitudes of the Northern Hemisphere (de Vernal et al.2005b, 2013b, 2020). The sea ice estimates from dinocysts used here are from different studies (see references in Table 1) and reconstructions based the new database, including 71 taxa and 1968 stations (de Vernal et al.2020). The reference sea ice data used for calibration are the monthly 1955–2012 average of the National Snow and Ice Data Center (NSIDC) (Walsh et al.2016). The results are expressed in term of annual mean of sea ice cover concentration or as the number of months with > 50 % of sea ice. The error of prediction for sea ice concentration is ±12 % and that of sea ice cover duration through the year is ±1.5 months yr−1. Such values are very close to the interannual variability in areas occupied by seasonal sea ice cover (see de Vernal et al.2013b).

Our biomarker approach for sea ice reconstruction is based on the determination of a highly branched isoprenoid (HBI) with 25 carbons (C25 HBI monoene = IP25) (Belt et al.2007). This biomarker is only biosynthesised by specific diatoms living in the Arctic sea ice (Brown et al.2014), meaning the presence of IP25 in the sediments is a direct proof for the presence of past Arctic sea ice. Meanwhile, this biomarker approach has been used successfully in numerous studies dealing with the reconstruction of past Arctic sea ice conditions during the late Miocene to Holocene (for a review, see Belt2018). By combining the sea ice proxy IP25 with (biomarker) proxies for open-water (phytoplankton productivity such as brassicasterol, dinosterol or a specific tri-unsaturated HBI, HBI-III), the so-called PIP25 index has been developed (Müller et al.2011; Belt et al.2015; Smik et al.2016). Based on a comparison (“calibration”) PIP25 data obtained from surface sediments with modern satellite-derived (spring) sea ice concentration maps (Müller et al.2011; Xiao et al.2015; Smik et al.2016), the PIP25 approach may allow a more semi-quantitative reconstruction of present and past Arctic Ocean sea ice conditions from marine sediments, i.e. estimates of spring sea ice concentration (or in the Central Arctic probably more the summer situation due to light limitations for algae growth in the other seasons). Based on these data, one may separate “permanent to extended sea ice cover” (> 0.75) and “seasonal sea ice cover”; (0.75–0.1), perhaps including the sub-groups “ice-edge” (0.75–0.5) and “less/reduced sea ice” (0.5–0.1), and “ice-free” (< 0.1). For pros and cons of this approach, we refer to a recent review by Belt (2018).

Figure 1Map showing the location of the LIG Arctic sediment cores listed in Table 1. Open symbols correspond to records with uncertain chronology, and filled symbols correspond to records with good chronology. The map background has been created using (last access: 1 January 2020).

Based on several IP25/PIP25 records obtained from central Arctic Ocean sediment cores (see Fig. 1 for core locations and Table 1 for data), perennial sea ice cover probably existed during the LIG in the Central Arctic, whereas along the Barents Sea continental margin, influenced by the inflow of warm Atlantic Water, sea ice was significantly reduced (Stein et al.2017). However, Stein et al. (2017) emphasise that the PIP25 records obtained from the central Arctic Ocean cores indicating a perennial sea ice cover have to be interpreted cautiously as the biomarker concentrations are very low to absent (see Belt2018 for further discussion). The productivity of algal material (ice and open water) must have been quite low, so that (almost) nothing reached the seafloor or is preserved in the sediments, and there must have been periods during the LIG when some open-water conditions occurred, since subpolar foraminifers and coccoliths were found in core PS51/038 and PS2200 (Stein et al.2017). It is however unclear whether these periods equate to more than 1 month yr−1 of open water (or seasonal ice conditions). This explains why some sites show both seasonal and perennial interpretations at the same site. The reader is referred to the original publications (Table 1) for more information on these data. Furthermore and importantly, a new revised 230Th chronology of late Quaternary sequences from the central Arctic Ocean (Hillaire-Marcel et al.2017) questions the age model of some of the data listed in Table 1. Thus, further verification of age control is still needed and the data from the central Arctic Ocean should be interpreted with caution. We have therefore marked the chronological control as “uncertain” for these cores, while the chronological control is good for cores outside the Central Arctic.

The information given by the different types of sea ice indicators shows that care should be taken when comparing them with model results. We have used the qualitative information given in Table 1, taking into account the threshold given in this table. Indeed, for instance, “perennial sea ice cover” does not automatically mean 100 % sea ice cover, or a sea ice concentration (SIC) of 1.0. It means rather that there is sea ice but not necessarily at a concentration of 100 % over the core site throughout the year (i.e. the summer season is not totally ice-free). Most qualitative reconstructions cite a threshold of 75 %, which we have therefore used in our model-data comparison. We have also used the quantitative mean annual sea ice reconstructions. Finally, similar to studies for future climate, we have considered the Arctic to be ice-free when, in any given month, the total area of sea ice is less than 1×106 km2. This means that some marine core sites could remain ice covered for the summer, but the Arctic would nevertheless remain technically ice-free.

2.3 CMIP6-PMIP4 models

The last Coupled Model Intercomparison Project Phase 5 (CMIP5) collected climate simulations performed with 60 different numerical models by 26 research institutes around the world (IPCC2013). The follow-on CMIP6 archive, to be completed in 2020, is expected to gather model outputs from over 30 research institutes. Of these, currently 15 models have run the CMIP6-PMIP4 LIG simulation (Table 2). We present results here from all these models.

Ziehn et al.Sidorenko et al. (2015)Sidorenko et al. (2019, 2015)Danabasoglu et al. (2020)Voldoire et al.Hazeleger et al. (2012)Li et al. (2020)Zheng et al.Kelley et al.Williams et al.Volodin et al.Boucher et al.Goosse et al.Hajima et al.20182020Cao et al.Guo et al.Seland et al.

Table 2Overview of models that have run the CMIP6-PMIP4 LIG simulation. For each model, denomination, physical core components, horizontal and vertical grid specifications, details on prescribed vs. interactive boundary conditions, reference publication, and LIG simulation length are shown.

Download XLSX

Table 2 provides an overview of the models used in this study. They are state-of-the-art coupled general circulation models (GCMs) and Earth System Models (ESMs) simulating the atmosphere, ocean, sea ice and land surface processes dynamics with varying degrees of complexity. These 15 CMIP6-PMIP4 models have been developed for several years by individual institutes across the world and, in the context of CMIP6, are used in the same configuration to seamlessly simulate past, present and future climate. We have added the results from the LOVECLIM Earth System Model of Intermediate Complexity, which can be used for longer simulations.

Table 2 shows the following qualities for each model: model denomination, physical core components, horizontal and vertical grid specifications, details on prescribed vs. interactive boundary conditions, relative publication for an in-depth model description, and LIG simulation length (spin-up and production runs).

2.4 PMIP4 LIG (lig127k) simulation protocol

Results shown here are from the main Tier 1 LIG simulation, from the standard CMIP6-PMIP4 LIG experimental protocol (Otto-Bliesner et al.2017). The prescribed LIG (lig127k) protocol differs from the CMIP6 pre-industrial (PI) simulation protocol in astronomical parameters and the atmospheric trace greenhouse gas concentrations (GHG). LIG astronomical parameters are prescribed according to Berger and Loutre (1991), and atmospheric trace GHG concentrations are based on ice core measurements. Table 3 from Otto-Bliesner et al. (2017) summarises the protocol. All models followed this protocol, except CNRM-CM6-1 for which the most important forcings for the LIG, i.e. the astronomical parameters, have been imposed at the recommended values, but the GHG have been kept at their pre-industrial values of 284.3170 ppm for CO2, 808.2490 ppb for CH4 and 273.0211 ppb for N2O. All other boundary conditions, including solar activity, ice sheets, aerosol emissions, etc., are identical to PI protocol. Both the Greenland and Antarctica ice sheets are known to have shrunk during the interglacial with different timings, and therefore taking PI characteristics for the lig127k protocol is an approximation, particularly for the Antarctic ice sheet, which was possibly smaller than in the PI at that time (Otto-Bliesner et al.2017). The Greenland ice sheet likely reached a minimum at around 120 ka and was probably still close to its PI size at 127 ka. Given the dating uncertainties and the difficulty for models to include the largest changes in ice sheets for 127 ka, i.e. changes in West Antarctica, the choice of the PMIIP4 working group on interglacials was to use the PI ice sheets as boundary conditions for the Tier 1 PMIP4-CMIP6 experiments presented here and to foster sensitivity experiments to ice sheet characteristics at a later stage. In terms of the Greenland ice sheet, the approximation is considered quite good and ideal for starting transient experiments through the whole interglacial period.

Table 3Astronomical parameters and atmospheric trace gas concentrations used to force LIG and PI simulations.

Download Print Version | Download XLSX

LIG simulations were initialised either from a previous LIG run, or from the standard CMIP6 protocol pre-industrial simulations, using constant 1850 GHGs, ozone, solar, tropospheric aerosol, stratospheric volcanic aerosol and land use forcing.

Although PI and LIG spin-ups vary between the models, most model groups aimed to allow the land and oceanic masses to attain approximate steady state i.e. to reach atmospheric equilibrium and to achieve an upper-oceanic equilibrium. LIG production runs are all between 100–200 years long, which is generally within the appropriate length for Arctic sea ice analysis (Guarino et al.2020).

The LIG orbital parameters result in modifications of the definitions of the months and seasons (in terms of start and end dates within a year). Since daily data was not available for all models to re-compute LIG-specific monthly averages, we have corrected these averages using the method of Bartlein and Shafer (2019). Unless otherwise specified, we use these results adjusted for the LIG calendar throughout this paper.

2.5 The CMIP6 1pctCO2 protocol

We compare the response to the lig127k forcings to idealised forcings for future climate. We have chosen to use the 1pctCO2 simulation from the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima Eyring et al.2016). These simulations start from the PI (piControl) experiment and the atmospheric CO2 concentration is gradually increased by 1 % yr−1 for at least 150 years, i.e. 10 years after atmospheric CO2 quadrupling.

3 Results: simulated Arctic sea ice

Since all LIG production runs are at least 100 years in length, all model results are averaged over at least 100 years. We refer to the multi-model mean throughout as the MMM. We consider both the sea ice area (SIA), defined as the sum, over all Northern Hemisphere ocean cells, of the sea ice concentration × the cell area and the sea ice extent (SIE), defined as the sum of the areas of ocean cells where the sea ice concentration is larger than 0.15. Both quantities are used in sea ice studies, SIE has been used widely in IPCC AR5 (Vaughan et al.2013), while SIA tends to be used more for CMIP6 analyses (e.g. SIMIP Community2020).

3.1 PI sea ice

For the present-day we have satellite and in situ observations with which to evaluate the models. The use of present-day sea ice data implies that we might expect the simulated PI sea ice to be generally somewhat larger than the observed mean. Indeed the atmospheric CO2 levels for the years for which we chose the observation dataset (1982 to 2001) were between 340 and 370 ppm, compared to the PI level of 280 ppm. Figure 2 shows the mean seasonal cycle of the Arctic sea ice extent simulated for the PI and LIG alongside the observed Arctic sea ice extent.

Figure 2Mean seasonal cycle of the Arctic sea ice area (SIA, left-hand side) and sea ice extent (SIE, right-hand side), in 106 km2, simulated for the PI and LIG periods by the PMIP4 models. The top row shows the results for PI. The grey shading shows the monthly minimum and maximum in the SIA and SIE observed over the years 1982–2001, as given by the NOAA_OI_v2 dataset. The second and fourth row show the LIG results, with no calendar adjustment and with calendar adjustment, respectively. The third and bottom row show the corresponding LIG–PI anomalies, with no calendar adjustment and with calendar adjustment, respectively.


Table 4Sea ice area and extent (in 106 km2) for the PI and LIG simulations (calendar-adjusted values). MMM stands for the multi-model mean, SD for the multi-model standard deviation. NA stands for not available.

Download Print Version | Download XLSX

The summer minimum monthly MMM SIA for the PI is 6.46 ± 1.41×106 km2, compared to the observed 1981 to 2002 mean of 5.65×106 km2. In terms of SIE, the summer minimum for PI is 8.89 ± 1.41×106 km2, to be compared to the observed 7.73×106 km2. Interestingly this MMM PI area and extent is a little larger than the 1981–2002 area. The majority of the simulations show a realistic representation of the geographical extent for the summer minimum (Fig. 3, Table 4), with 9 out of 16 models showing a slightly smaller area compared to the present-day observations and 7 showing an overestimated area. LOVECLIM, EC-Earth, FGOALS-g3, GISS-E2-1-G and INM-CM4-8 clearly simulate too much ice (Table 4). The other models generally exhibit realistic PI summer minimum ice conditions. The detail of the geographical distribution of sea ice for the models, the MMM and the NOAA_OI_v2 datasets (Fig. 3) confirms the results in terms of Arctic sea ice extent. Overestimations appear to be due to too much sea ice being simulated in the Barents–Kara area (LOVECLIM, FGOALS-g3, GISS-E2-1-G), in the Nordic Seas (EC-Earth, FGOALS-g3) and in Baffin Bay (LOVECLIM, INM-CM4-8, EX-Earth). MIROC-ES2L performs rather poorly for the PI, with insufficient ice close to the continents. The other models generally match the 0.15 isoline from the NOAA_OI_v2 dataset in a realistic manner. The winter maximum monthly MMM areas show little difference between the present-day and PI simulated areas. The MMM PI area is 15.16 ± 1.90×106 km2, compared to the observed 1981 to 2002 mean of 14.44×106 km2. For both the summer and winter, the simulations and observations mostly agree on the month that the minimum and maximum are attained: August–September for the minimum and February–March for the maximum for every model.

Figure 3PI sea ice concentration for the month of minimum SIA as computed for Fig. 2. The magenta contour shows the 0.15 isocontour of the NOAA_OI_v2 observations (Reynolds et al.2002, see the Data availability section) averaged over the years 1982–2001.

Before we carry out the comparison between model results and sea ice cover reconstructions for the LIG period, we compare the results of the models for PI to the observations at the reconstruction sites (Fig. 4 for the comparison of annual mean sea ice concentrations and Fig. 5a and b for winter and summer). Models generally overestimate sea ice cover at the three northernmost sites in summer and in annual mean and over the seven northernmost sites for the winter season. Those sites are actually very close to the sea ice edge and the overestimation could correspond to the fact that the observations are for 1981 to 2002 period, which was already warmer than the pre-industrial one.

Figure 4Sea ice annual concentration simulated for PI, for the multi-model mean (MMM) and for each model. The colour filling of the symbols on the maps correspond to the observed values at each site, which are classified into three categories according to the NOAA_OI_v2 dataset: perennial cover (9 to 12 months), seasonal cover (3 to 9 months) and ice-free state (0 to 3 months). On the MMM panel, for each data site the colour of the symbol outline corresponds to the number of models simulating the observed ice cover. On the panels for individual models, the shape of the symbol depends on the observed result being below the simulated one (triangle down), above the simulated one (triangle up) or in the same category as the simulated one (circle). The number of sites for which reconstructions are equal to and above the number of months simulated by models are written at the bottom-right corner of each panel.

Figure 5Model–data comparison as a function of latitude and record site for PI (a, b) and LIG (c, d). For each LIG data site, the NOAA_OI_v2 observations (PI) or reconstructions (LIG) are shown in the first column of each plot, and the model results are in the columns to the right. Both the model results and the NOAA_OI_v2 observations are shown in terms of sea ice fraction averaged over the month of minimum (a, c) or maximum (b, d) Northern Hemisphere sea ice area and the previous and following months. For the LIG, the qualitative assessments (eighth column of Table 1) have been used for records with good chronological control. The letter next to the name of the site stands for indicator used for the reconstruction: dinocysts (“d”) or IP25/PIP25 (“i”). For summer conditions, dark blue shading is used for the “no sea ice” category, light blue shading is used for “summer probably ice-free” conditions, and white shading is used for “summer sea ice concentration > 75 %” and “perennial sea ice”. For winter, dark blue shading is used for “ice-free all year round” conditions and white shading is used for “seasonal sea ice conditions” and “perennial sea ice”. The model results are averaged as they are for PI and shown following the colour scale on the right-hand side of the plots.


3.2 LIG sea ice

The models show a minimum monthly MMM SIA for the LIG of 3.20 ± 1.50×106 km2, and a maximum MMM SIA of 15.95 ± 2.61×106 km2. In terms of SIE, the minimum MMM extent is 5.39 ± 2.13×106 km2, while the maximum MMM extent is equal to 18.38 ± 3.12×106 km2. Thus, compared to the PI results, there is a reduction of ca. 50 % in the MMM minimum (summer) monthly SIA in the LIG results, and of nearly 40 % in terms of SIE, but a slight increase in the winter monthly MMM SIA and SIE. Every model shows an often substantial reduction in summer sea ice between the PI and LIG.

There is a large amount of inter-model variability for the LIG SIA and SIE during the summer (Fig. 6 and Table 4). Out of the 16 models, 1 model, HadGEM3, shows a LIG Arctic Ocean free of sea ice in summer, i.e. with an SIE lower than 1×106 km2. CESM2 and NESM3 show low SIA values (slightly above 2×106 km2) in summer for the LIG simulation, but their minimum SIE values are below 3×106 km2. Both HadGEM3 and CESM2 realistically capture the PI Arctic sea ice seasonal cycle. On the other hand, NESM3 overestimates winter ice and the amplitude of the seasonal cycle in SIA and SIE, while simulating realistic PI values for both SIA and SIE (Cao et al.2018). This seasonal cycle is amplified in the LIG simulation, with an increase in SIA and SIE in winter and a decrease in summer, following the insolation forcing. Hence, the difference in the response of these models to LIG forcing in terms of sea ice does not appear to only depend on differences in PI sea ice representation.

Figure 6LIG sea ice concentration for the month of minimum SIA (computed with calendar adjustment) as computed for Fig. 2. The magenta contour shows the 0.15 isocontour of the corresponding PI simulation.

For the winter, only one model (EC-Earth) simulates a decrease in SIA and SIE of around 2 × 106 km2, two other models (ACCESS and INM-CM4-8) simulate a slight decrease in SIA and SIE, all other models simulate an increase in both SIA and SIE. All models therefore show a larger sea ice area amplitude for LIG than for PI, and the range of model results is larger for LIG than for PI. The summer season and the seasons of sea ice growth and decay are therefore key to understanding the behaviour of LIG sea ice and the inter-model differences, as will be confirmed in Sect. 4.

3.3 LIG model–data comparison

We limit our comparison to the sites for which the chronology is good. These cores mostly show ice-free conditions in summer, except for the northernmost site (core PS92/039-2), which is at least 75 % covered by ice in summer (Fig. 5c). Two other sites at high latitude (PS213861 and PS93/006-1, for which sea ice has been reconstructed based on dinocysts and IP25/PIP25), show summer conditions which are “probably ice-free”. Only four models simulate more than 75 % sea ice concentration over the northernmost site, but they also simulate more than 75 % sea ice concentration over the two following sites (in descending order of latitudes), and FGOALS-g3 simulates more than 75 % sea ice concentration for another four sites for which the reconstructions show no sea ice. On the other hand, 10 models simulate no sea ice concentration at all over the reconstruction sites in summer, and therefore probably overestimate the LIG summer sea ice reduction. From these reconstructions, we cannot distinguish the performance of the models simulating a strong reduction of sea ice from the model simulating a nearly total disappearance of summer sea ice in the Arctic. Apart from FGOALS-g3, which simulates extensive sea ice cover for both periods, there does not appear to be a strong relationship between the PI and LIG model results over the data sites: models which simulate sea ice cover over the three northernmost sites at the LIG do not necessarily simulate large sea ice concentrations over these sites for PI (e.g. LOVECLIM, AWIESM1 and AWIESM2).

For the winter season, the reconstructions show the four northernmost sites to be ice covered. The reconstructions for most other sites are qualitatively given as “nearly ice-free all year round” or “ice-free all year round”. Model results are generally in agreement with the reconstructions for the three to four northernmost sites (Fig. 5d). Most models simulate sea ice over some of the sites characterised by “nearly ice-free all year round” conditions, and only one model (IPSLCM6) simulates sea ice cover over a site for which the reconstructions show ice-free conditions. The model–data agreement is therefore quite good for the winter season. In this case, the model results for LIG appear to be related to their results for PI, with models simulating more sea ice for LIG being those simulating more sea ice for PI.

Figure 7 shows a quantitative model–data comparison in terms of annual mean sea ice concentration, which is the variable for which we have the highest number of reconstructions (Table 1). From this, we see that it is more difficult for the models to realistically capture sea ice change over the core sites near Greenland close to the sea ice edge. If we cross-compare the observation–model match for each model for both the PI (Fig. 4) and the LIG (Fig. 7) then FGOALS-g3 and NESM3 have difficulties in accurately capturing sea ice cover at the core site locations in the Nordic Seas, whilst AWIESM1 and NORESM1-F best display sea ice cover close to the sea ice edge near Greenland and in the Nordic seas for both time periods. It is these Nordic Seas sea ice edge differences (over the core sites listed in Table 1) that make the difference between the simulation–data matches for each model.

Figure 7Sea ice annual concentration during the LIG (computed with calendar adjustment), for the multi-model mean (MMM) and for each model. The colour-filling of the symbols on the maps correspond to the reconstructed values, classified into three categories: perennial cover (9 to 12 months), seasonal cover (3 to 9 months) and ice-free state (0 to 3 months). On the MMM panel, for each data site, the colour of the symbol outline corresponds to the number of models simulating the reconstructed ice cover. On the panels for individual models, the shape of the symbol depends on the model result being below the reconstructed one (triangle down), above the reconstructed one (triangle up) or in the same category as the reconstructed one (circle). The number of data points which are above, equal and below the number of months simulated by models are written in the bottom-right corner of each panel.

4 Discussion of model differences

Whilst we cannot yet definitely establish the most likely Arctic sea ice conditions during the LIG, we can investigate sea ice differences across models when we have sufficient model data. We have first performed this analysis for the three models for which we had sufficient data: CESM2, HadGEM3, and IPSLM6. These models each represent a distinct sea ice response to the LIG forcing, i.e. summer sea ice concentration less than 0.15 everywhere (HadGEM3), significant summer sea ice retreat with concentration less than 0.8 in the Central Arctic (CESM2) and modest summer sea ice retreat with a small area with sea ice concentration close to 1 in the Central Arctic (IPSLCM6).

Sea ice formation and melting can be affected by a large number of factors inherent to the atmosphere and the ocean dynamics, alongside the representation of sea ice itself within the model (i.e. the type of sea ice scheme used). In coupled models it can be extremely difficult to identify the causes of essentially coupled model behaviour. Nevertheless, we discuss the short-wave (SW) surface energy balance, ocean, and atmosphere circulations and comment on cloudiness and albedo changes.

4.1 Atmospheric energy budget differences

The atmospheric energy budget LIG–PI anomaly (Fig. 8) is negative in winter and strongly positive in summer, following the imposed insolation anomaly. These anomalies in total heat budget are dominated by the SW budget contribution from May to August. We split the SW budget into the downward (SWdn) and upward (SWup) contributions. Both fluxes are defined to be positive when they are downward and negative when they are upward. Hence, the total SW budget (in black) is the sum of the SWdn contribution (in red) and the SWup contribution (in blue). In this figure, a positive SWup anomaly means that the SWup is less intense at LIG than at PI, hence contributing to an increase in the net SW flux.

Figure 8The main components of the atmospheric energy budget at the surface averaged over the Arctic (70–90 N) for HadGEM3, CESM2 and IPSLCM6. The LIG–PI anomalies as shown as a function of the month for the total energy budget (Ftot, black), the SW budget (SW, violet), and for the downward (SWdn, red) and upward SW (SWup, blue) fluxes. The sign convention for all fluxes is that fluxes pointing downward are positive and fluxes pointing upward are negative. Panels (b), (d) and (f) show results for which the LIG calendar has been taken into account (for the LIG simulations), while panels (a), (c) and (e) show the results averaged on the PI calendar both for PI and LIG.


For all the models, the total heat budget anomaly is due to (i) an increased downward short-wave flux in spring resulting from the insolation forcing and (ii) a decreased upward short-wave flux in summer, related to the decrease of the albedo due to the smaller sea ice cover. During summer, this decrease in upward short-wave flux more than compensates the decrease in SWdn, which is maximum in August.

Figure 9Albedo over the Arctic for PI (a, b), LIG (c, d) and LIG–PI (e, f) for HadGEM3, IPSL-CM6 and CESM2. The albedo has been recomputed from the SWup and SWdn fluxes. Panels (a), (c) and (e) show the results for the whole Arctic, while panels (b), (d) and (f) show the results for areas where the sea ice fraction is larger than 0.9. All LIG values have been calendar adjusted.


The summer anomaly reaches 80 W m−2 in June for HadGEM3, 60 W m−2 for IPSLCM6 and 50 W m−2 for CESM2. The differences between the model results are due to a different phasing of the SWdn and SWup anomalies for HadGEM3, compared to the other two models: for HadGEM3, the two fluxes peak in June, while for CESM2 and IPSLCM6, the SWdn flux peaks in May and the SWup signal peaks in July, and thus the anomaly in these fluxes partly compensate for each other. HadGEM3 shows a larger net SW increase despite a SWdn anomaly which is smaller than for the other two models. On the other hand, HadGEM3's SWup component is stronger and always positive, which is different to the other two models, which show a negative SWup contribution in April–May. These differences are associated with differences in albedo for the three models (Fig. 9). HadGEM3's sea ice and Arctic Ocean albedos are always smaller than those simulated by IPSLCM6 and CESM2 and the difference is larger for LIG than for PI. The albedo simulated by HadGEM3 in May and June is particularly low compared to the two other models, which explains why the SWup component peaks earlier. The albedo LIG–PI anomalies over the whole Arctic show that the sea ice albedo feedback is most effective in HadGEM3.

In terms of cloudiness, IPSLCM6 shows differences in the properties of clouds, in terms of optical depth, between PI and LIG, but this could not be investigated due to a lack of data (thus far) for the other models. Thus we cannot tell if LIG–PI anomalies in SWdn fluxes, i.e. differences between HadGEM3's and CESM2 flux, also have a contribution due to cloud changes.

Figure 10The main components of the atmospheric energy budget averaged over the Arctic (70–90 N) for HadGEM3, CESM2, CNRM-CM6-1, IPSLCM6, INM-CM4-8, MIROC-ES2L, NorESM1 and NorESM2. The LIG–PI anomalies are shown as a function of the month for the total energy budget (black), the SW budget (violet), and for the downward (red) and upward SW (blue) fluxes. The sign convention for all fluxes is the same as for Fig. 8.


The comparison to other model results (Fig. 10) confirm that the behaviour of HadGEM3 is unusual in terms of energy budget. It is the only model in which the anomalies in SWup and SWdn are exactly in phase and produce a much larger anomaly in total heat budget, while in other models those anomalies are not in phase and partly compensate each other.

4.2 Ocean and atmosphere circulation differences

Changes in Arctic sea ice related to ocean heat transport have been found for the CESM large ensemble (Auclair and Tremblay2018). The differences can then be amplified by the sea ice albedo feedbacks. We check this in our models by calculating long-term means of the maximum meridional stream function at 26 N for the PI and LIG simulations. These are 19.5 and 18.7 for CESM2, 15.6 and 15.8 Sv for HadGEM3, and 12.9 and 10.4 for IPSLCM6. Thus, the CESM2 and HadGEM3 models exhibit an Atlantic Meridional Overturning Circulation (AMOC) that is almost unchanged between PI and LIG, while in the IPSLCM6 model the AMOC weakens. This implies that a reduced northward oceanic heat transport could prevent sea ice loss in the Central Arctic in some but not all models (see also Stein et al.2017).

Some differences in the response of sea ice to LIG forcing therefore appear to be due either to differences in atmospheric response (HadGEM3 vs. IPSL-CM6 and CESM2), similar to mechanisms found for current sea ice decline (e.g. He et al.2019; Olonscheck et al.2019) or to changes in ocean heat transport (CESM2 vs. IPSLCM6). But while AMOC changes partially explain the differences found between IPSL (more sea ice in Central Arctic) and CESM2 and HadGEM3 (less sea ice in Central Arctic), they do not explain differences between ice-free and ice-covered conditions in HadGEM3 and CESM2.

Differences in atmospheric circulation changes could also explain difference in sea ice response to LIG forcings. We therefore investigate LIG–PI anomalies in sea level pressure (Fig. 11). Most models simulate a decrease in summer mean sea level pressure largely encompassing the Arctic Ocean and adjacent continents. This decrease is not as strong over the Nordic Seas as it is over the Arctic, and this local heterogeneity over the Nordic Seas is model dependent. However, the anomaly in atmospheric circulation is more zonal over the Nordic Seas and northern North Atlantic in HadGEM3 than in CESM2 or IPSLCM6, and therefore differences in atmospheric circulation are probably not causing more warm air to enter the Arctic for HadGEM3 and are thus not the cause of HadGEM3 being so warm over the Arctic. The mean sea level pressure winter anomaly is characterised by a deepening of the Icelandic low for all models except NESM3.

Figure 11Anomalies (LIG–PI, hPa, shading) in mean sea level pressure for DJF (top plots) and JJA (bottom plots). Contours indicating PI values are superimposed: values every 5 hPa, 1005 hPa isobar in white, black contours for lower values and grey contours for higher values.

Other factors that remain to be investigated include clouds and ocean heat uptake in the Arctic in the different models, e.g. as a function of stratification.

4.3 Transient CO2 forced responses: LIG vs. transient 1pctCO2

The LIG has higher insolation than PI at high northern latitudes during spring and summer and less significant changes in winter insolation. This is distinct from the increased GHG, which is the dominant forcing for future climates. However, since sea ice minimum occurs in summer, it is of interest to consider possible relationships between CMIP6 model responses for the LIG and those for the transient 1pctCO2 experiments. A total of 12 models have the LIG, PI and 1pctCO2 simulations available. These include models with large, small and intermediate responses in sea ice for the LIG.

Figure 12 suggests that there is indeed such a relationship between the summer sea ice concentration decreases for LIG and the averages from years 50 to 70 of the transient 1pctCO2 simulations: the models that respond strongly at the LIG also respond strongly for the 1pctCO2 forcing, and the model with the smallest response for the LIG (INMCM4-8) has the smallest response to the 1pctCO2 forcing. The relationship shown in Fig. 12 does not last for later periods in the 1pctCO2 runs, when the winter sea ice is also affected by the increased greenhouse gas forcing. This implies inter-comparisons between the LIG simulation and simulations with moderate CO2 increase (during the transition to high CO2 levels) should be investigated.

Figure 12LIG vs. 1pctCO2 July–August–September sea ice areas (for sea ice concentrations larger than 0.15). The results for the 1pctCO2 simulations have been averaged for years 50 to 70.


5 Conclusions

The Last Interglacial period (LIG) was the last time global temperature was substantially higher than the pre-industrial period at high northern latitudes (Otto-Bliesner et al.2013; Capron et al.2017; Otto-Bliesner et al.2017; Fischer et al.2018; Otto-Bliesner et al.2021). To help understand the role of Arctic sea ice in these changes, we present a new synthesis of LIG sea ice information using marine core data collected in the Arctic Ocean, Nordic Seas and northern North Atlantic and compare this to PMIP4-LIG simulations.

Our synthesis shows that south of 79 N in the Atlantic and Nordic seas the LIG was definitely seasonally ice-free. These southern sea ice records provide quantitative estimates of sea surface parameters based on dinoflagellate cysts (dinocysts). North of 79 N the sea-ice-related records are more difficult to obtain and interpret. However, the core at 81.5 N brings evidence of summer being probably seasonally ice-free during the LIG from two indicators: dinocysts and IP25/PIP25. The northernmost core with good chronology is located at 81.9 N and shows evidence of substantial (> 75 %) sea ice concentration all year round. Other cores, with debated chronologies, have not been used for model–data comparisons in the present study.

Model results from 16 models show a multi model mean (MMM) summer SIA LIG of 3.20 ± 1.29×106 km2, and a winter monthly MMM area of 15.95 ± 1.21×106 km2. This is a reduction in SIA of 50 % for the minimum summer month between the PI and LIG but almost no change for the winter month MMM. Every model shows an often substantial reduction in summer sea ice between the PI and LIG. For the winter, only one of the 16 models shows a (small) winter reduction in sea ice between the PI and LIG. This reinforces that the key seasons for understanding LIG warming are the spring, summer and autumn.

We investigate reasons for inter-model differences in LIG Arctic sea ice simulations: we find that the LIG total heat budget anomaly in the Arctic is due to (i) an increased downward short-wave flux in spring, resulting from the insolation forcing, and (ii) a decreased upward short-wave flux in summer, related to the decrease of the albedo due to the smaller sea ice cover. During summer, this decrease in upward short-wave flux more than compensates the decrease in the SWdn, which is at a maximum in August. Differences between the model results are due to a different phasing of the up and down short-wave anomalies in the different models and are associated with the differences in model albedo.

Analysis of IPSLCM6 results shows differences in the properties of clouds, in terms of optical depth, between PI and LIG. Further work is required to identify if this is also important for other models. Changes in Arctic sea ice may also be related to ocean heat transport. Here, we have shown that ocean circulation changes occur for some (but not all) LIG simulations. Other factors that remain to be investigated include clouds and ocean heat uptake in the Arctic in the different models.

Most models agree with the reconstructed year-round ice-free northern North Atlantic. Model–data disagreement for the LIG occur over the Nordic Seas, close to Greenland and at the boundary with the Arctic Ocean, where many models overestimate annual mean sea ice concentration. This is not fully related to the model performance for summer. Indeed, 12 of 16 models simulate little sea ice cover over the northernmost site and 10 of the models simulate less than 25 % sea ice concentration over the site at 81.5 N. It is not possible, from the available data, to decide on the best models, in particular in terms of summer sea ice. The northernmost site appears to discriminate those models that simulate very little sea ice at this site. However, models which do simulate > 75 % summer sea ice concentration at this site also simulate > 75 % summer sea ice concentration for the two sites at 81.5 and 79.2 N, just south of the northernmost site, which is not realistic. More reconstructions with good chronology are needed in the Central Arctic to determine which model behaviour is more realistic, and in particular if the summer ice-free Arctic simulated by the HadGEM3 model alone, among the 16 models, is possible. This would be key in assessing ESMs used for future projections with respect to climates with much warmer summers than today. This means that it is all the more crucial that there appear to be a nearly linear relationship between the ESM simulations of summer sea ice for the near future (years 50 to 70 of transient 1pctCO2 simulations) and that simulated for the LIG: the models which respond strongly to the LIG forcing also respond strongly for the 1pctCO2 forcing. This implies inter-comparisons between the LIG simulation and simulations with a moderate CO2 increase (during the transition to high CO2 levels) may yield insight into likely 21st century Arctic sea ice changes, especially if we achieve a more extensive characterisation of LIG Arctic sea ice from marine cores.

Data availability

The NOAA Optimum Interpolation (OI) V2 dataset for sea ice concentration has been retrieved from (NOAA/OAR/ESRL PSL2019). The dataset used for the present study is the monthly dataset: (last access: 8 December 2019).

The original output data from the model simulations used in this study are available from the Earth System Grid Federation (, last access: 31 May 2020), the data repository for CMIP6 simulations, or on open repositories, as listed in this paper's companion paper by Otto-Bliesner et al. (2021).

Nonetheless, the exact data shown in Figs. 2 to 12 are also provided as a Supplement to this paper: the numbers are given as text files for Figs. 2, 8, 9, 10 and 12. NetCDF files are provided with the data shown on Figs. 3, 4, 6, 7 and 11. For each model, there is one netCDF file ( with sea ice variables (Figs. 3, 4, 6 and 7) and one netCDF file ( with the mean sea level pressure (Fig. 11).


The supplement related to this article is available online at:

Author contributions

MK, LCS, MS, MVG, RS and AdV are joint first authors for this paper. MK and LCS planned the study with the other QUIGS members. MK and MS analysed all model simulations and produced all model figures. LCS wrote the paper. MVG contributed substantially to the first draft and compiled all model information. AdV, IMV, RS and LCS compiled the sea ice dataset, and IMV produced the dataset map. DS co-planned some of the model analysis. MK, LCS, MS, MVG, AAO, PB, ECB, JC, MAC, CG, ANL, GL, KJM, LM, PM, KHN, BLOB, RO, SRB, DSyM, SST, XS, BS, RAT, EV, NKHY, WP, QZ, ZZ, TZ contributed model data. DS, DF, CB, and JS provide sea ice modelling advice. MK supervised the corrections with contributions from LCS, MS, DS, CG and BLOB. All authors read the draft and commented on the text.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4) (CP/GMD inter-journal SI)”. It is not associated with a conference.


We acknowledge the QUIGS (Quaternary Interglacials working group endorsed by PAGES and PMIP) for making this comparison possible, thanks in particular to the workshop organised by this group in Cambridge, UK, in July 2019. We are grateful to the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP6. We thank the climate modelling groups for producing and making their model output available, the Earth System Grid Federation (ESGF) for archiving the data and providing access to them, and the multiple funding agencies who support CMIP6 and ESGF. The Paleoclimate Modelling Intercomparison Project is thanked for coordinating the lig127k protocol and making the model–model and model–data comparisons possible within CMIP6. PMIP is endorsed by WCRP and CLIVAR. We also acknowledge NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, for their optimally interpolated sea ice product, downloaded from their website at Qiong Zhang acknowledges the HPC resources provided by the Swedish National Infrastructure for Computing (SNIC) at the National Supercomputer Centre (NSC). Bette L. Otto-Bliesner, Esther C. Brady and Robert A. Tomas acknowledge the CESM project, which is supported primarily by the National Science Foundation (NSF). This material is based upon work supported by the National Center for Atmospheric Research (NCAR). Computing and data storage resources, including the Cheyenne supercomputer (, Computational and Information Systems Laboratory2019), were provided by the Computational and Information Systems Laboratory (CISL) at NCAR. The Last Interglacial studies of Anne de Vernal have been supported by the Natural Sciences and Engineering Research Council of Canada and the “Fonds de recherche du Québec – Nature et technologies”. The NorESM simulations benefitted from resources provided by UNINETT Sigma2 – the National Infrastructure for High Performance Computing and Data Storage in Norway. The ACCESS-ESM 1.5 experiments were performed on Raijin at the NCI National Facility at the Australian National University, through awards under the National Computational Merit Allocation Scheme, the Intersect allocation scheme, and the UNSW HPC at NCI Scheme. LOVECLIM experiments were performed on UNSW HPC Katana. Ayako Abe-Ouchi, Ryouta O'ishi, and Sam Sherriff-Tadano thank JAMSTEC for use of the Earth Simulator supercomputer.

Financial support

This research has been supported by the PAlaeo-Constraints on Monsoon Evolution and Dynamics (PACMEDY) Belmont Forum project (grand no. 01LP1607A), German Federal Ministry of Education and Science (BMBF) PalMod II WP 3.3 (grant no. 01LP1924B), NERC (projects NE/P013279/1 and NE/P009271/1), European Union's Horizon 2020 research and innovation programme (grant agreement no. 820970), “Convention des Services Climatiques” from IPSL, Russian state (assignment project 0148-2019-0009), RSF (grant no. 20-17-00190), NSF (cooperative agreement no. 1852977), Australian Research Council (grant nos. FT180100606 and DP180100048), Arctic Challenge for Sustainability (ArCS) Project (grant no. JPMXD1300000000), Arctic Challenge for Sustainability II (ArCS II) Project (grant no. JPMXD1420318865), JSPS KAKENHI (grant no. 17H06104), MEXT KAKENHI (grant no. 17H06323), Natural Sciences and Engineering Research Council of Canada, and the “Fonds de recherche du Québec – Nature et technologies”.

Review statement

This paper was edited by Marie-France Loutre and reviewed by Julie Brigham-Grette and one anonymous referee.


Adler, R. E., Polyak, L., Ortiz, J. D., Kaufman, D. S., Channell, J. E. T., Xuan, C., Grottoli, A. G., Selln, E., and Crawford, K. A.: Sediment record from the western Arctic Ocean with an improved Late Quaternary age resolution: HOTRAX core HLY0503-8JPC, Mendeleev Ridge, Global Planet. Change, 68, 18–29,, 2009. a, b

Auclair, G. and Tremblay, L. B.: The role of ocean heat transport in rapid sea ice declines in the Community Earth System Model Large Ensemble., J. Geophys. Res.-Oceans, 123, 8941–8957,, 2018. a

Bartlein, P. J. and Shafer, S. L.: Paleo calendar-effect adjustments in time-slice and transient climate-model simulations (PaleoCalAdjust v1.0): impact and strategies for data analysis, Geosci. Model Dev., 12, 3889–3913,, 2019. a

Belt, S.: Source-specific biomarkers as proxies for Arctic and Antarctic sea ice, Org. Geochem., 125, 277–298,, 2018. a, b, c

Belt, S. T., Massé, G., Rowland, S. J., Poulin, M., Michel, C., and LeBlanc, B.: A novel chemical fossil of palaeo sea ice: IP25, Org. Geochem., 38, 16–27,, 2007. a

Belt, S. T., Cabedo-Sanz, P., Smik, L., Navarro-Rodriguez, A., Berben, S. M., Knies, J., and Husum, K.: Identification of paleo Arctic winter sea ice limits and the marginal ice zone: Optimised biomarker-based reconstructions of late Quaternary Arctic sea ice, Earth Planet. Sc. Lett., 431, 127–139,, 2015. a

Berger, A. and Loutre, M.-F.: Insolation values for the climate of the last 10 million years, Quaternary Sci. Rev., 10, 297–317, 1991. a

Boucher, O., Servonnat, J., Albright, A. L., Aumont, O., Balkanski, Y., Bastrikov, V., Bekki, S., Bonnet, R., Bony, S., Bopp, L., Braconnot, P., Brockmann, P., Cadule, P., Caubel, A., Cheruy, F., Cozic, A., Cugnet, D., D'Andrea, F., Davini, P., de Lavergne, C., Denvil, S., Dupont, E., Deshayes, J., Devilliers, M., Ducharne, A., Dufresne, J.-L., Ethé, C., Fairhead, L., Falletti, L., Foujols, M.-A., Gardoll, S., Gastineau, G., Ghattas, J., Grandpeix, J.-Y., Guenet, B., Guez, L., Guilyardi, E., Guimberteau, M., Hauglustaine, D., Hourdin, F., Idelkadi, A., Joussaume, S., Kageyama, M., Khadre-Traoré, A., Khodri, M., Krinner, G., Lebas, N., Levavasseur, G., Lévy, C., Lott, F., Lurton, T., Luyssaert, S., Madec, G., Madeleine, J.-B., Maignan, F., Marchand, M., Marti, O., Mellul, L., Meurdesoif, Y., Mignot, J., Musat, I., Ottlé, C., Peylin, P., Planton, Y., Polcher, J., Rio, C., Rousset, C., Sepulchre, P., Sima, A., Swingedouw, D., Thiéblemont, R., Vancoppenolle, M., Vial, J., Vialard, J., Viovy, N., and Vuichard, N.: Presentation and evaluation of the IPSL-CM6A-LR climate model, J. Adv. Model. Earth Sy., 12, e2019MS002010,, 2020. a

Brigham-Grette, J. and Hopkins, D. M.: Emergent marine record and paleoclimate of the last interglaciation along the northwest Alaskan coast, Quaternary Res., 43, 159–173,, 1995. a

Brown, T., Belt, S., Tatarek, A., and Mundy, C. J.: Source identification of the Arctic sea ice proxy IP25, Nat. Commun., 5, 4197,, 2014. a

Cao, J., Wang, B., Yang, Y.-M., Ma, L., Li, J., Sun, B., Bao, Y., He, J., Zhou, X., and Wu, L.: The NUIST Earth System Model (NESM) version 3: description and preliminary evaluation, Geosci. Model Dev., 11, 2975–2993,, 2018. a, b

Capron, E., Govin, A., Feng, R., Otto-Bliesner, B. L., and Wolff, E. W.: Critical evaluation of climate syntheses to benchmark CMIP6/PMIP4 127 ka Last Interglacial simulations in the high-latitude regions, Quaternary Sci. Rev., 168, 137–150,, 2017. a, b

Computational and Information Systems Laboratory: Cheyenne: HPE/SGI ICE XA System (Climate Simulation Laboratory), National Center for Atmospheric Research, Boulder, CO, USA,, 2019. a

Cronin, T., Gemery, L., Briggs Jr., W., Jakobsson, M., Polyak, L., and Brouwers, E.: Quaternary Sea-ice history in the Arctic Ocean based on a new Ostracode sea-ice proxy, Quaternary Sci. Rev., 29, 3415–3429,, 2010. a, b, c, d

Danabasoglu, G., Lamarque, J. F., Bachmeister, J., Bailey, D. A., DuVivier, A. K., Edwards, J., Emmons, L. K., Fasullo, J., Garcia, R., Gettelman, A., Hannay, C., Holland, M. M., Large, W. G., Lauritzen, P. H., Lawrence, D. M., Lenaerts, J. T. M., Lindsay, K., Lipscomb, W. H., Mills, M. J., Neale, R., Oleson, K. W., Otto-Bliesner, B., Phillips, A. S., Sacks, W., Tilmes, S., van Kampenhout, L., Vertenstein, M., Bertini, A., Dennis, J., Deser, C., Fischer, C., Fox-Kemper, B., Kay, J. E., Kinnison, D., Kushner, P. J., Larson, V. E., Long, M. C., Mickelson, S., Moore, J. K., Nienhouse, E., Polvani, L., Rasch, P. J., and Strand, W. G.: The Community Earth System Model version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001916., 2020. a

de Vernal, A. and Hillaire-Marcel, C.: Natural Variability of Greenland Climate, Vegetation and Ice Volume during the Last Million Years, Science, 320, 1622–1625,, 2008. a

de Vernal, A., Eynaud, F., Henry, M., Hillaire-Marcel, C., Londeix, L., Mangin, S., Matthiessen, J., Marret, F., Radi, T., Rochon, A., Solignac, S., and Turon, J.-L.: Reconstruction of sea-surface conditions at middle to high latitudes of the Northern Hemisphere, PANGAEA,, 2005a. a

de Vernal, A., Hillaire-Marcel, C., and Darby, D. A.: Variability of sea ice cover in the Chukchi Sea (western Arctic Ocean) during the Holocene, Paleoceanography, 20, PA4018,, 2005b. a

de Vernal, A., Gersonde, R., Goosse, H., Seidenkrantz, M.-S., and Wolff, E. W.: Sea ice in the paleoclimate system: the challenge of reconstructing sea ice from proxies – an introduction, Quaternary Sci. Rev., 79, 1–8,, 2013a. a

de Vernal, A., Rochon, A., Fréchette, B., Henry, M., Radi, T., and Solignac, S.: Reconstructing past sea ice cover of the Northern Hemisphere from dinocyst assemblages: status of the approach, Quaternary Sci. Rev., 79, 122–134,, 2013b. a, b, c, d

de Vernal, A., Radi, T., Zaragosi, S., Nieuwenhove, N. V., Rochon, A., Allana, E., Scheppere, S. D., Eynaud, F., Head, M. J., Limoges, A., Londeix, L., Marret, F., Matthiessen, J., Penaud, A., Pospelova, V., Price, A., and Richerol, T.: Distribution of common modern dinoflagellate cyst taxa in surface sediments of the Northern Hemisphere in relation to environmental parameters: The new n=1968 database, Mar. Micropaleontol., 159, 101796,, 2020. a, b, c

Eynaud, F.: Kystes de Dinoflagellés et Evolution paléoclimatique et paléohydrologique de l'Atlantique Nord au cours du Dernier Cycle Climatique du Quaternaire, PhD thesis, Bordeaux 1 University, France, 1999. a

Eynaud, F., Turon, J., Sánchez-Goñi, M., and Gendreau, S.: Dinoflagellate cyst evidence of `Heinrich-like events' off Portugal during the Marine Isotopic Stage 5, Mar. Micropaleontol., 40, 9–21,, 2000. a, b

Eynaud, F., Turon, J., and Duprat, J.: Comparison of the Holocene and Eemian palaeoenvironments in the South Icelandic Basin: dinoflagellate cysts as proxies for the North Atlantic surface circulation, Rev. Palaeobot. Palyno., 128, 55–79,, 2004. a, b

Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958,, 2016. a, b

Fischer, H., Meissner, K. J., Mix, A. C., Abram, N. J., Austermann, J., Brovkin, V., Capron, E., Colombaroli, D., Daniau, A.-L., Dyez, K. A., Felis, T., Finkelstein, S. A., Jaccard, S. L., McClymont, E. L., Rovere, A., Sutter, J., Wolff, E. W., Affolter, S. , Bakker, P., Ballesteros-Canovas, J. A., Barbante, C., Caley, T., Carlson, A. E., Churakova (Sidorova), O., Cortese, G., Cumming, B. F., Davis, B. A. S., de Vernal, A., Emile-Geay, J., Fritz, S. C., Gierz, P., Gottschalk, J., Holloway, M. D., Joos, F., Kucera, M., Loutre, M.-F., Lunt, D. J., Marcisz, K., Marlon, J. R., Martinez, P., Masson-Delmotte, V., Nehrbass-Ahles, C., Otto-Bliesner, B. L., Raible, C. C., Risebrobakken, B., Sanchez-Goni, M.-F. , Arrigo, J. S., Sarnthein, M., Sjolte, J., Stocker, T. F., Alvarez, P. A. V., Tinner, W., Valdes, P. J., Vogel, H., Wanner, H., Yan, Q., Yu, Z., Ziegler, M., and Zhou, L.: Palaeoclimate constraints on the impact of 2 C anthropogenic warming and beyond, Nat. Geosci., 11, 474–485,, 2018. a, b

Goosse, H., Brovkin, V., Fichefet, T., Haarsma, R., Huybrechts, P., Jongma, J., Mouchet, A., Selten, F., Barriat, P.-Y., Campin, J.-M., Deleersnijder, E., Driesschaert, E., Goelzer, H., Janssens, I., Loutre, M.-F., Morales Maqueda, M. A., Opsteegh, T., Mathieu, P.-P., Munhoven, G., Pettersson, E. J., Renssen, H., Roche, D. M., Schaeffer, M., Tartinville, B., Timmermann, A., and Weber, S. L.: Description of the Earth system model of intermediate complexity LOVECLIM version 1.2, Geosci. Model Dev., 3, 603–633,, 2010. a

Guarino, M.-V., Sime, L. C., Schroeder, D., Lister, G. M. S., and Hatcher, R.: Machine dependence and reproducibility for coupled climate simulations: the HadGEM3-GC3.1 CMIP Preindustrial simulation, Geosci. Model Dev., 13, 139–154,, 2020. a

Guo, C., Bentsen, M., Bethke, I., Ilicak, M., Tjiputra, J., Toniazzo, T., Schwinger, J., and Otterå, O. H.: Description and evaluation of NorESM1-F: a fast version of the Norwegian Earth System Model (NorESM), Geosci. Model Dev., 12, 343–362,, 2019. a

Hajima, T., Watanabe, M., Yamamoto, A., Tatebe, H., Noguchi, M. A., Abe, M., Ohgaito, R., Ito, A., Yamazaki, D., Okajima, H., Ito, A., Takata, K., Ogochi, K., Watanabe, S., and Kawamiya, M.: Development of the MIROC-ES2L Earth system model and the evaluation of biogeochemical processes and feedbacks, Geosci. Model Dev., 13, 2197–2244,, 2020. a

Hazeleger, W., Wang, X., Severijns, C., Ştefănescu, S., Bintanja, R., Sterl, A., Wyser, K., Semmler, T., Yang, S., Van den Hurk, B., van Noije, T., van der Linden, E., and van der Wiel, K.: EC-Earth V2. 2: description and validation of a new seamless earth system prediction model, Clim. Dynam., 39, 2611–2629, 2012. a

He, M., Hu, Y., Chen, N., Wang, D., Huang, J., and Stamnes, K.: High cloud coverage over melted areas dominates the impact of clouds on the albedo feedback in the Arctic., Sci. Rep.-UK, 9, 9529,, 2019. a

Hillaire-Marcel, C., de Vernal, A., Bilodeau, G., and Weaver, A. J.: Absence of deep-water formation in the Labrador Sea during the last interglacial period, Nature, 410, 1073–1077, 2001. a, b

Hillaire-Marcel, C., de Vernal, A., and McKay, J.: Foraminifer isotope study of the Pleistocene Labrador Sea, northwest North Atlantic (IODP Sites 1302/03 and 1305), with emphasis on paleoceanographical differences between its “inner” and “outer” basins, Mar. Geol., 279, 188–198,, 2011. a

Hillaire-Marcel, C., Ghaleb, B., de Vernal, A., Maccali, J., Cuny, K., Jacobel, A., Duc, C. L., and McManus, J.: A new chronology of late quaternary sequences from the central Arctic Ocean based on “extinction ages” of their excesses in 231Pa and 230Th, Geochem., Geophy., Geosy., 18, 4573–4585,, 2017. a

Hodell, D., Minth, E., Curtis, J., Mccave, I., Hall, I., Channell, J., and Xuan, C.: Surface and deep-water hydrography on Gardar Drift (Iceland Basin) during the Last Interglacial period, Earth Planet. Sci. Lett., 288, 10–19,, 2009. a

IPCC: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, Cambridge, United Kingdom and New York, USA, 2007. a

IPCC: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y. Bex, V., and Midgley, P. M., Cambridge, United Kingdom and New York, NY, USA,, 2013. a, b, c

Kageyama, M., Braconnot, P., Harrison, S. P., Haywood, A. M., Jungclaus, J. H., Otto-Bliesner, B. L., Peterschmitt, J.-Y., Abe-Ouchi, A., Albani, S., Bartlein, P. J., Brierley, C., Crucifix, M., Dolan, A., Fernandez-Donado, L., Fischer, H., Hopcroft, P. O., Ivanovic, R. F., Lambert, F., Lunt, D. J., Mahowald, N. M., Peltier, W. R., Phipps, S. J., Roche, D. M., Schmidt, G. A., Tarasov, L., Valdes, P. J., Zhang, Q., and Zhou, T.: The PMIP4 contribution to CMIP6 – Part 1: Overview and over-arching analysis plan, Geosci. Model Dev., 11, 1033–1057,, 2018. a

Kelley, M., Schmidt, G. A., Nazarenko, L., Bauer, S. E., Ruedy, R., Russell, G. L., Ackerman, A. S., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Cesana, G., Cheng, Y., Clune, T. L., Cook, B. I., Cruz, C. A., Del Genio, A. D., Elsaesser, G. S., Faluvegi, G., Kiang, N. Y., Kim, D., Lacis, A. A., Leboissetier, A., LeGrande, A. N., Lo, K. K., Marshall, J., Matthews, E. E., McDermid, S., Mezuman, K., Miller, R. L., Murray, L. T., Oinas, V., Orbe, C., Pérez García-Pando, C., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A., Shindell, D. T., Sun, S., Tausnev, N., Tsigaridis, K., Tselioudis, G., Weng, E., Wu, J., and Yao, M.-S.: GISS-E2.1: Configurations and climatology, J. Adv. Model. Earth Sy., 12, e2019MS002 025,, 2020. a

Kremer, A., Stein, R., Fahl, K., Bauch, H., Mackensen, A., and Niessen, F.: A 190-ka biomarker record revealing interactions between sea ice, Atlantic Water inflow and ice sheet activity in eastern Fram Strait, Arktos, 4, 1–17,, 2018a. a, b

Kremer, A., Stein, R., Fahl, K., Ji, Z., Yang, Z., Wiers, S., Matthiessen, J., Forwick, M., Löwemark, L., O'Regan, M., Chen, J., and Snowball, I.: Changes in sea ice cover and ice sheet extent at the Yermak Plateau during the last 160 ka – Reconstructions from biomarker records, Quaternary Sci. Rev., 182, 93–108,, 2018b. a, b

Li, L., Yu, Y., Tang, Y., Lin, P., Xie, J., Song, M., Dong, L., Zhou, T., Liu, L., Wang, L., Pu, Y., Chen, X., Chen, L., Xie, Z., Liu, H., Zhang, L., Huang, X., Feng, T., Zheng, W., Xia, K., Liu, H., Liu, J., Wang, Y., Wang, L., Jia, B., Xie, F., Wang, B., Zhao, S., Yu, Z., Zhao, B., and Wei, J.: The flexible global ocean‐atmosphere‐land system model grid‐point version 3 (fgoals‐g3): description and evaluation, J. Adv. Model. Earth Syst., 12, e2019MS002012,, 2020. a

Lunt, D. J., Abe-Ouchi, A., Bakker, P., Berger, A., Braconnot, P., Charbit, S., Fischer, N., Herold, N., Jungclaus, J. H., Khon, V. C., Krebs-Kanzow, U., Langebroek, P. M., Lohmann, G., Nisancioglu, K. H., Otto-Bliesner, B. L., Park, W., Pfeiffer, M., Phipps, S. J., Prange, M., Rachmayani, R., Renssen, H., Rosenbloom, N., Schneider, B., Stone, E. J., Takahashi, K., Wei, W., Yin, Q., and Zhang, Z. S.: A multi-model assessment of last interglacial temperatures, Clim. Past, 9, 699–717,, 2013. a, b, c

Malmierca-Vallet, I., Sime, L. C., Valdes, P. J., Capron, E., Vinther, B. M., and Holloway, M. D.: Simulating the Last Interglacial Greenland stable water isotope peak: The role of Arctic sea ice changes, Quaternary Sci. Rev., 198, 1–14,, 2018. a, b

Masson-Delmotte, V., Braconnot, P., Hoffmann, G., Jouzel, J., Kageyama, M., Landais, A., Lejeune, Q., Risi, C., Sime, L., Sjolte, J., Swingedouw, D., and Vinther, B.: Sensitivity of interglacial Greenland temperature and δ18O: ice core data, orbital and increased CO2 climate simulations, Clim. Past, 7, 1041–1059,, 2011. a

Matthiessen, J. and Knies, J.: Dinoflagellate cyst evidence for warm interglacial conditions at the northern Barents Sea margin, during marine isotope stage 5, J. Quaternary Sci., 16, 727–737, 2001. a

Matthiessen, J., Knies, J., Nowaczyk, N. R., and Ruediger, S.: Age determination of sediment core PS2138-1, PANGAEA,, 2001. a, b

Müller, J., Wagner, A., Fahl, K., Stein, R., Prange, M., and Lohmann, G.: Towards quantitative sea ice reconstructions in the northern North Atlantic: A combined biomarker and numerical modelling approach, Earth Planet. Sci. Lett., 306, 137–148,, 2011. a, b

NOAA/OAR/ESRL PSL: NOAA Optimum Interpolation (OI) Sea Surface Temperature (SST) V2, Boulder, Colorado, USA, available at:, last access: 8 December 2019. a

Nørgaard-Pedersen, N., Mikkelsen, N., Lassen, S. J., Kristoffersena, Y., and Sheldon, E.: Reduced sea ice concentrations in the Arctic Ocean during the last interglacial period revealed by sediment cores off northern Greenland, Paleoceanography, 22, PA1218,, 2007. a, b

O'ishi, R., Chan, W.-L., Abe-Ouchi, A., Sherriff-Tadano, S., and Ohgaito, R.: PMIP4/CMIP6 Last Interglacial simulations using different versions of MIROC, with and without vegetation feedback, Clim. Past Discuss.,, in review, 2020. a

Olonscheck, D., Mauritsen, T., and Notz, D.: Arctic sea-ice variability is primarily driven by atmospheric temperature fluctuations., Nat. Geosci., 12, 430–434,, 2019. a

Otto-Bliesner, B. L., Marshall, S. J., Overpeck, J. T., Miller, G. H., Hu, A., and CAPE Last Interglacial Project members: Simulating Arctic Climate Warmth and Icefield Retreat in the Last Interglaciation, Science, 311, 1751–1753,, 2006. a

Otto-Bliesner, B. L., Rosenbloom, N., Stone, E. J., McKay, N. P., Lunt, D. J., Brady, E. C., and Overpeck, J. T.: How warm was the last interglacial? New model-data comparisons, Philos. T. Roy. Soc. A, 371, 20130097,, 2013. a, b, c, d, e, f

Otto-Bliesner, B. L., Braconnot, P., Harrison, S. P., Lunt, D. J., Abe-Ouchi, A., Albani, S., Bartlein, P. J., Capron, E., Carlson, A. E., Dutton, A., Fischer, H., Goelzer, H., Govin, A., Haywood, A., Joos, F., LeGrande, A. N., Lipscomb, W. H., Lohmann, G., Mahowald, N., Nehrbass-Ahles, C., Pausata, F. S. R., Peterschmitt, J.-Y., Phipps, S. J., Renssen, H., and Zhang, Q.: The PMIP4 contribution to CMIP6 – Part 2: Two interglacials, scientific objective and experimental design for Holocene and Last Interglacial simulations, Geosci. Model Dev., 10, 3979–4003,, 2017. a, b, c, d, e, f, g

Otto-Bliesner, B. L., Brady, E. C., Zhao, A., Brierley, C. M., Axford, Y., Capron, E., Govin, A., Hoffman, J. S., Isaacs, E., Kageyama, M., Scussolini, P., Tzedakis, P. C., Williams, C. J. R., Wolff, E., Abe-Ouchi, A., Braconnot, P., Ramos Buarque, S., Cao, J., de Vernal, A., Guarino, M. V., Guo, C., LeGrande, A. N., Lohmann, G., Meissner, K. J., Menviel, L., Morozova, P. A., Nisancioglu, K. H., O'ishi, R., Salas y Melia, D., Shi, X., Sicard, M., Sime, L., Stepanek, C., Tomas, R., Volodin, E., Yeung, N. K. H., Zhang, Q., Zhang, Z., and Zheng, W.: Large-scale features of Last Interglacial climate: results from evaluating the lig127k simulations for the Coupled Model Intercomparison Project (CMIP6)–Paleoclimate Modeling Intercomparison Project (PMIP4), Clim. Past, 17, 63–94,, 2021. a, b

Penaud, A., Eynaud, F., Turon, J., Zaragosi, S., Marret, F., and Bourillet, J.: Interglacial variability (MIS 5 and MIS 7) and dinoflagellate cyst assemblages in the Bay of Biscay (North Atlantic), Mar. Micropaleontol., 68, 136–155,, 2008. a

Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C., and Wang, W.: An improved in situ and satellite SST analysis for climate, J. Climate, 15, 1609–1625, 2002. a, b

Seland, Ø., Bentsen, M., Seland Graff, L., Olivié, D., Toniazzo, T., Gjermundsen, A., Debernard, J. B., Gupta, A. K., He, Y., Kirkevåg, A., Schwinger, J., Tjiputra, J., Schancke Aas, K., Bethke, I., Fan, Y., Griesfeller, J., Grini, A., Guo, C., Ilicak, M., Hafsahl Karset, I. H., Landgren, O., Liakka, J., Onsum Moseid, K., Nummelin, A., Spensberger, C., Tang, H., Zhang, Z., Heinze, C., Iverson, T., and Schulz, M.: The Norwegian Earth System Model, NorESM2 – Evaluation of theCMIP6 DECK and historical simulations, Geosci. Model Dev. Discuss.,, in review, 2020. a

Sidorenko, D., Rackow, T., Jung, T., Semmler, T., Barbi, D., Danilov, S., Dethloff, K., Dorn, W., Fieg, K., Gößling, H. F., Handorf, D., Harig, S., Hiller, W., Juricke, S., Losch, M., Schröter, J., Sein D. V., and Wang, Q.: Towards multi-resolution global climate modeling with ECHAM6-FESOM. Part I: model formulation and mean climate, Clim. Dynam., 44, 757–780, 2015. a, b

Sidorenko, D., Goessling, H., Koldunov, N., Scholz, P., Danilov, S., Barbi, D., Cabos, W., Gurses, O., Harig, S., Hinrichs, C., Juricke, S., Lohmann, G., Losch, M., Mu, L., Rackow, T., Rakowsky, N., Sein, D., Semmler, T., Shi, X., Stepanek, C., Streffing, J., Wang, Q., Wekerle, C., Yang, H., and Jung, T.: Evaluation of FESOM2.0 Coupled to ECHAM6.3: Preindustrial and HighResMIP Simulations, J. Adv. Model. Earth Sy., 11, 3794–3815,, 2019. a

Sime, L. C., Risi, C., Tindall, J. C., Sjolte, J., Wolff, E. W., Masson-Delmotte, V., and Capron, E.: Warm climate isotopic simulations: what do we learn about interglacial signals in Greenland ice cores?, Quaternary Sci. Rev., 67, 59–80,, 2013. a

SIMIP Community: Arctic Sea Ice in CMIP6, Geophys. Res. Lett., 47, e2019GL086749,, 2020. a

Smik, L., Belt, S. T., Lieser, J. L., Armand, L. K., and Leventer, A.: Distributions of highly branched isoprenoid alkenes and other algal lipids in surface waters from East Antarctica: Further insights for biomarker-based paleo sea-ice reconstruction, Org. Geochem., 95, 71–80,, 2016. a, b

Stein, R., Fahl, K., Gierz, P., Niessen, F., and Lohmann, G.: Arctic Ocean sea ice cover during the penultimate glacial and the last interglacial, Nat. commun., 8, 373,, 2017. a, b, c, d, e, f, g, h, i, j, k

Tatebe, H., Tanaka, Y., Komuro, Y., and Hasumi, H.: Impact of deep ocean mixing on the climatic mean state in the Southern Ocean, Sci. Rep., 8, 14479,, 2018. a

Van Nieuwenhove, N. and Bauch, H. A.: Last interglacial (MIS 5e) surface water conditions at the Voring Plateau (Norwegian Sea), based on dinoflagellate cysts, Polar Res., 27, 175–186,, 2008. a

Van Nieuwenhove, N., Bauch, H. A., and Matthiessen, J.: Last interglacial surface water conditions in the eastern Nordic Seas inferred from dinocyst and foraminiferal assemblages, Mar. Micropaleontol., 66, 247–263,, 2008. a

Van Nieuwenhove, N., Bauch, H. A., Eynaud, F., Kandiano, E., Cortijo, E., and Turon, J.-L.: Evidence for delayed poleward expansion of North Atlantic surface waters during the last interglacial (MIS 5e), Quaternay Sci. Rev., 30, 934–946,, 2011. a, b, c

Van Nieuwenhove, N., Bauch, H. A., and Andruleit, H.: Multiproxy fossil comparison reveals contrasting surface ocean conditions in the western Iceland Sea for the last two interglacials, Palaeogeogr. Palaeocl., 370, 247–259,, 2013. a

Vaughan, D., Comiso, J., Allison, I., Carrasco, J., Kaser, G., Kwok, R., Mote, P., Murray, T., Paul, F., Ren, J., Rignot, E., Solomina, O., Steffen, K., and Zhang, T.: Observations: Cryosphere, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P., Cambridge University Press, Cambridge, UK, 2013. a

Voldoire, A., Saint‐Martin, D., Sénési, S., Decharme, B., Alias, A., Chevallier, M., Colin, J., Guérémy, J.‐F., Michou, M., Moine, M.‐P., Nabat, P., Roehrig, R., Salas, D., Séférian, R., Valcke, S., Beau, I., Belamari, S., Berthet, S., Cassou, C., Cattiaux, J., Deshayes, J., Douville, H., Ethé, C., Franchistéguy, L., Geoffroy, O., Lévy, C., Madec, G., Meurdesoif, Y., Msadek, R., Ribes, A., Sanchez‐Gomez, E., Terray, L., and Waldman, R.: Evaluation of CMIP6 DECK Experiments With CNRM-CM6.1, J. Adv. Model. Earth Sy., 11, 2177–2213,, 2019. a

Volodin, E. M., Mortikov, E. V., Kostrykin, S. V., Galin, V. Y., Lykossov, V. N., Gritsun, A. S., Diansky, N. A., Gusev, A. V., Iakovlev, N. G., Shestakova, A. A., and Emelina, S. V.: Simulation of the modern climate using the INM-CM48 climate model, Russ. J. Numer. Anal. M., 33, 367–374, 2018. a

Walsh, J. E., Chapman, W. L., and Fetterer, F.: Gridded Monthly Sea Ice Extent and Concentration, 1850 Onward, Version 1, National Snow and Ice Data Center,, 2016.  a

Williams, K., Copsey, D., Blockley, E., Bodas-Salcedo, A., Calvert, D., Comer, R., Davis, P., Graham, T., Hewitt, H., Hill, R., Hyder, P., Ineson, S., Johns, T. C., Keen, A. B., Lee, R. W., Megann, A., Milton, S. F., Rae, J. G. L., Roberts, M. J., Scaife, A. A., Schiemann, R., Storkey, D., Thorpe, L., Watterson, I., G., Walters, D. N., West, A., Wood, R. A., Woollings, T., and Xavier, P. K.: The Met Office global coupled model 3.0 and 3.1 (GC3.0 and GC3.1) configurations, J. Adv. Model. Earth Sy., 10, 357–380, 2018. a

Xiao, X., Stein, R., and Fahl, K.: MIS 3 to MIS 1 temporal and LGM spatial variability in Arctic Ocean sea ice cover: Reconstruction from biomarkers, Paleoceanography, 30, 969–983, 2015. a

Zheng, W., Yu, Y., Luan, Y., Zhao, S., He, B., Dong, L., Song, M., Lin, P., and Liu, H.: CAS-FGOALS Datasets for the Two Interglacial Epochs of the Holocene and the Last Interglacial in PMIP4, Adv. Atmos. Sci., 37, 1034–1044,, 2020. a

Zhuravleva, A., Bauch, H. A., and Van Nieuwenhove, N.: Last Interglacial (MIS5e) hydrographic shifts linked to meltwater discharges from the East Greenland margin, Quaternary Sci. Rev., 164, 95–109,, 2017. a, b

Ziehn, T., Lenton, A., Law, R. M., Matear, R. J., and Chamberlain, M. A.: The carbon cycle in the Australian Community Climate and Earth System Simulator (ACCESS-ESM1) – Part 2: Historical simulations, Geosci. Model Dev., 10, 2591–2614,, 2017. a

Short summary
The Last interglacial (ca. 127 000 years ago) is a period with increased summer insolation at high northern latitudes, resulting in a strong reduction in Arctic sea ice. The latest PMIP4-CMIP6 models all simulate this decrease, consistent with reconstructions. However, neither the models nor the reconstructions agree on the possibility of a seasonally ice-free Arctic. Work to clarify the reasons for this model divergence and the conflicting interpretations of the records will thus be needed.