Search

We combine consistently dated benthic carbon isotopic records distributed over the entire Atlantic Ocean with numerical simulations performed by a glacial configuration of the Norwegian Earth System Model with active ocean biogeochemistry in order to interpret the observed Cibicides δ13C changes at the stadial–interstadial transition corresponding to the end of Heinrich Stadial 4 (HS4) in terms of ocean circulation and remineralization changes. We show that the marked increase in Cibicides δ13C observed at the end of HS4 between ∼2000 and 4200 m in the Atlantic can be explained by changes in nutrient concentrations as simulated by the model in response to the halting of freshwater input in the high-latitude glacial North Atlantic. Our model results show that this Cibicides δ13C signal is associated with changes in the ratio of southern-sourced (SSW) versus northern-sourced (NSW) water masses at the core sites, whereby SSW is replaced by NSW as a consequence of the resumption of deep-water formation in the northern North Atlantic and Nordic Seas after the freshwater input is halted. Our results further suggest that the contribution of ocean circulation changes to this signal increases from ∼40 % at 2000 m to ∼80 % at 4000 m. Below ∼4200 m, the model shows little ocean circulation change but an increase in remineralization across the transition marking the end of HS4. The simulated lower remineralization during stadials compared to during interstadials is particularly pronounced in deep subantarctic sites, in agreement with the decrease in the export production of carbon to the deep Southern Ocean during stadials found in previous studies.

Insoluble particles in ice cores record signatures of past climate parameters like vegetation dynamics, volcanic activity, and aridity. For some of them, the analytical detection relies on intensive bench microscopy investigation and requires dedicated sample preparation steps. Both are laborious, require in-depth knowledge, and often restrict sampling strategies. To help overcome these limitations, we present a framework based on flow imaging microscopy coupled to a deep neural network for autonomous image classification of ice core particles. We train the network to classify seven commonly found classes, namely mineral dust, felsic and mafic (basaltic) volcanic ash grains (tephra), three species of pollen (Corylus avellana, Quercus robur, Quercus suber), and contamination particles that may be introduced onto the ice core surface during core handling operations. The trained network achieves 96.8 % classification accuracy at test time. We present the system's potential and its limitations with respect to the detection of mineral dust, pollen grains, and tephra shards, using both controlled materials and real ice core samples. The methodology requires little sample material, is non-destructive, fully reproducible, and does not require any sample preparation procedures. The presented framework can bolster research in the field by cutting down processing time, supporting human-operated microscopy, and further unlocking the paleoclimate potential of ice core records by providing the opportunity to identify an array of ice core particles. Suggestions for an improved system to be deployed within a continuous flow analysis workflow are also presented.

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period in which the atmospheric CO2 concentration was approximately equal to the concentration we measure today (ca. 400 ppm). Sea surface temperature (SST) proxies indicate above-average warming over the North Atlantic in the mid-Pliocene with respect to the pre-industrial period, which may be linked to an intensified Atlantic Meridional Overturning Circulation (AMOC). Earlier results from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) show that the ensemble simulates a stronger AMOC in the mid-Pliocene than in the pre-industrial. However, no consistent relationship between the stronger mid-Pliocene AMOC and either the Atlantic northward ocean heat transport (OHT) or average North Atlantic SSTs has been found. In this study, we look further into the drivers and consequences of a stronger AMOC in mid-Pliocene compared to pre-industrial simulations in PlioMIP2. We find that all model simulations with a closed Bering Strait and Canadian Archipelago show reduced freshwater transport from the Arctic Ocean into the North Atlantic. This contributes to an increase in salinity in the subpolar North Atlantic and Labrador Sea that can be linked to the stronger AMOC in the mid-Pliocene. To investigate the dynamics behind the ensemble's variable response of the total Atlantic OHT to the stronger AMOC, we separate the Atlantic OHT into two components associated with either the overturning circulation or the wind-driven gyre circulation. While the ensemble mean of the overturning component is increased significantly in magnitude in the mid-Pliocene, it is partly compensated by a reduction in the gyre component in the northern subtropical gyre region. This indicates that the lack of relationship between the total OHT and AMOC is due to changes in OHT by the subtropical gyre. The overturning and gyre components should therefore be considered separately to gain a more complete understanding of the OHT response to a stronger mid-Pliocene AMOC. In addition, we show that the AMOC exerts a stronger influence on North Atlantic SSTs in the mid-Pliocene than in the pre-industrial, providing a possible explanation for the improved agreement of the PlioMIP2 ensemble mean SSTs with reconstructions in the North Atlantic.

Retreat of marine outlet glaciers often initiates depletion of inland ice through dynamic adjustments of the upstream glacier. The local topography of a fjord may promote or inhibit such retreat, and therefore fjord geometry constitutes a critical control on ice sheet mass balance. To quantify the processes of ice–topography interactions and enhance the understanding of the dynamics involved, we analyze a multitude of topographic fjord settings and scenarios using the Ice-sheet and Sea-level System Model (ISSM). We systematically study glacier retreat through a variety of artificial fjord geometries and quantify the modeled dynamics directly in relation to topographic features. We find that retreat in an upstream-widening or upstream-deepening fjord does not necessarily promote retreat, as suggested by previous studies. Conversely, it may stabilize a glacier because converging ice flow towards a constriction enhances lateral and basal shear stress gradients. An upstream-narrowing or upstream-shoaling fjord, in turn, may promote retreat since fjord walls or bed provide little stability to the glacier where ice flow diverges. Furthermore, we identify distinct quantitative relationships directly linking grounding line discharge and retreat rate to fjord topography and transfer these results to a long-term study of the retreat of Jakobshavn Isbræ. These findings offer new perspectives on ice–topography interactions and give guidance to an ad hoc assessment of future topographically induced ice loss based on knowledge of the upstream fjord geometry.

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period during which atmospheric CO2 levels were similar to recent historical values (∼400 ppm). Several proxy reconstructions for the mid-Pliocene show highly reduced zonal sea surface temperature (SST) gradients in the tropical Pacific Ocean, indicating an El Niño-like mean state. However, past modelling studies do not show these highly reduced gradients. Efforts to understand mid-Pliocene climate dynamics have led to the Pliocene Model Intercomparison Project (PlioMIP). Results from the first phase (PlioMIP1) showed clear El Niño variability (albeit significantly reduced) and did not show the greatly reduced time-mean zonal SST gradient suggested by some of the proxies. In this work, we study El Niño–Southern Oscillation (ENSO) variability in the PlioMIP2 ensemble, which consists of additional global coupled climate models and updated boundary conditions compared to PlioMIP1. We quantify ENSO amplitude, period, spatial structure and “flavour”, as well as the tropical Pacific annual mean state in mid-Pliocene and pre-industrial simulations. Results show a reduced ENSO amplitude in the model-ensemble mean (−24 %) with respect to the pre-industrial, with 15 out of 17 individual models showing such a reduction. Furthermore, the spectral power of this variability considerably decreases in the 3–4-year band. The spatial structure of the dominant empirical orthogonal function shows no particular change in the patterns of tropical Pacific variability in the model-ensemble mean, compared to the pre-industrial. Although the time-mean zonal SST gradient in the equatorial Pacific decreases for 14 out of 17 models (0.2 ∘C reduction in the ensemble mean), there does not seem to be a correlation with the decrease in ENSO amplitude. The models showing the most “El Niño-like” mean state changes show a similar ENSO amplitude to that in the pre-industrial reference, while models showing more “La Niña-like” mean state changes generally show a large reduction in ENSO variability. The PlioMIP2 results show a reasonable agreement with both time-mean proxies indicating a reduced zonal SST gradient and reconstructions indicating a reduced, or similar, ENSO variability.

The mid-Pliocene (∼3 Ma) is one of the most recent warm periods with high CO2 concentrations in the atmosphere and resulting high temperatures, and it is often cited as an analog for near-term future climate change. Here, we apply a moisture budget analysis to investigate the response of the large-scale hydrological cycle at low latitudes within a 13-model ensemble from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The results show that increased atmospheric moisture content within the mid-Pliocene ensemble (due to the thermodynamic effect) results in wetter conditions over the deep tropics, i.e., the Pacific intertropical convergence zone (ITCZ) and the Maritime Continent, and drier conditions over the subtropics. Note that the dynamic effect plays a more important role than the thermodynamic effect in regional precipitation minus evaporation (PmE) changes (i.e., northward ITCZ shift and wetter northern Indian Ocean). The thermodynamic effect is offset to some extent by a dynamic effect involving a northward shift of the Hadley circulation that dries the deep tropics and moistens the subtropics in the Northern Hemisphere (i.e., the subtropical Pacific). From the perspective of Earth's energy budget, the enhanced southward cross-equatorial atmospheric transport (0.22 PW), induced by the hemispheric asymmetries of the atmospheric energy, favors an approximately 1∘ northward shift of the ITCZ. The shift of the ITCZ reorganizes atmospheric circulation, favoring a northward shift of the Hadley circulation. In addition, the Walker circulation consistently shifts westward within PlioMIP2 models, leading to wetter conditions over the northern Indian Ocean. The PlioMIP2 ensemble highlights that an imbalance of interhemispheric atmospheric energy during the mid-Pliocene could have led to changes in the dynamic effect, offsetting the thermodynamic effect and, hence, altering mid-Pliocene hydroclimate.

The mid-Pliocene warm period (mPWP; ∼3.2 million years ago) is seen as the most recent time period characterized by a warm climate state, with similar to modern geography and ∼400 ppmv atmospheric CO2 concentration, and is therefore often considered an interesting analogue for near-future climate projections. Paleoenvironmental reconstructions indicate higher surface temperatures, decreasing tropical deserts, and a more humid climate in West Africa characterized by a strengthened West African Monsoon (WAM). Using model results from the second phase of the Pliocene Modelling Intercomparison Project (PlioMIP2) ensemble, we analyse changes of the WAM rainfall during the mPWP by comparing them with the control simulations for the pre-industrial period. The ensemble shows a robust increase in the summer rainfall over West Africa and the Sahara region, with an average increase of 2.5 mm/d, contrasted by a rainfall decrease over the equatorial Atlantic. An anomalous warming of the Sahara and deepening of the Saharan Heat Low, seen in >90 % of the models, leads to a strengthening of the WAM and an increased monsoonal flow into the continent. A similar warming of the Sahara is seen in future projections using both phase 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5). Though previous studies of future projections indicate a west–east drying–wetting contrast over the Sahel, PlioMIP2 simulations indicate a uniform rainfall increase in that region in warm climates characterized by increasing greenhouse gas forcing. We note that this effect will further depend on the long-term response of the vegetation to the CO2 forcing.

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ∼ 8–12 ∘C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

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.

The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to the Coupled Model Intercomparison Project (CMIP6) is the Tier 1 Last Interglacial experiment for 127 000 years ago (lig127k), designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models as for the future and following a common experimental protocol. Here we present a first analysis of a multi-model ensemble of 17 climate models, all of which have completed the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima) experiments. The equilibrium climate sensitivity (ECS) of these models varies from 1.8 to 5.6 ∘C. The seasonal character of the insolation anomalies results in strong summer warming over the Northern Hemisphere continents in the lig127k ensemble as compared to the CMIP6 piControl and much-reduced minimum sea ice in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than the CMIP6 midHolocene simulations as expected from the larger insolation anomalies at 127 than 6 ka. New synthesis for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and for precipitation over the Northern Hemisphere continents. The model–data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6–Paleoclimate Modeling Intercomparison Project (PMIP4) lig127k simulations, in combination with the proxy record, improve our confidence in future projections of monsoons, surface temperature, and Arctic sea ice, thus providing a key target for model evaluation and optimization.

This study presents simulations of Greenland surface melt for the Eemian interglacial period (∼130 000 to 115 000 years ago) derived from regional climate simulations with a coupled surface energy balance model. Surface melt is of high relevance due to its potential effect on ice core observations, e.g., lowering the preserved total air content (TAC) used to infer past surface elevation. An investigation of surface melt is particularly interesting for warm periods with high surface melt, such as the Eemian interglacial period. Furthermore, Eemian ice is the deepest and most compressed ice preserved on Greenland, resulting in our inability to identify melt layers visually. Therefore, simulating Eemian melt rates and associated melt layers is beneficial to improve the reconstruction of past surface elevation. Estimated TAC, based on simulated melt during the Eemian, could explain the lower TAC observations. The simulations show Eemian surface melt at all deep Greenland ice core locations and an average of up to ∼30 melt days per year at Dye-3, corresponding to more than 600 mm water equivalent (w.e.) of annual melt. For higher ice sheet locations, between 60 and 150 mmw.e.yr-1 on average are simulated. At the summit of Greenland, this yields a refreezing ratio of more than 25 % of the annual accumulation. As a consequence, high melt rates during warm periods should be considered when interpreting Greenland TAC fluctuations as surface elevation changes. In addition to estimating the influence of melt on past TAC in ice cores, the simulated surface melt could potentially be used to identify coring locations where Greenland ice is best preserved.

This paper provides the first chronology for the deep ice core from the East Greenland Ice-core Project (EGRIP) over the Holocene and the late last glacial period. We rely mainly on volcanic events and common peak patterns recorded by dielectric profiling (DEP) and electrical conductivity measurement (ECM) for the synchronization between the EGRIP, North Greenland Eemian Ice Drilling (NEEM) and North Greenland Ice Core Project (NGRIP) ice cores in Greenland. We transfer the annual-layer-counted Greenland Ice Core Chronology 2005 (GICC05) from the NGRIP core to the EGRIP ice core by means of 381 match points, typically spaced less than 50 years apart. The NEEM ice core has previously been dated in a similar way and is only included to support the match-point identification. We name our EGRIP timescale GICC05-EGRIP-1. Over the uppermost 1383.84 m, we establish a depth–age relationship dating back to 14 967 years b2k (years before the year 2000 CE). Tephra horizons provide an independent validation of our match points. In addition, we compare the ratio of the annual layer thickness between ice cores in between the match points to assess our results in view of the different ice-flow patterns and accumulation regimes of the different periods and geographical regions. For the next years, this initial timescale will be the basis for climatic reconstructions from EGRIP high-resolution proxy data sets, e.g. stable water isotopes, chemical impurity or dust records.

Palaeoclimate simulations improve our understanding of the climate, inform us about the performance of climate models in a different climate scenario, and help to identify robust features of the climate system. Here, we analyse Arctic warming in an ensemble of 16 simulations of the mid-Pliocene Warm Period (mPWP), derived from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The PlioMIP2 ensemble simulates Arctic (60–90∘ N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 ∘C compared to the pre-industrial period, with a multi-model mean (MMM) increase of 7.2 ∘C. The Arctic warming amplification ratio relative to global SAT anomalies in the ensemble ranges from 1.8 to 3.1 (MMM is 2.3). Sea ice extent anomalies range from −3.0 to -10.4×106 km2, with a MMM anomaly of -5.6×106 km2, which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer sea-ice-free conditions (≤1×106 km2) in their mPWP simulation. The ensemble tends to underestimate SAT in the Arctic when compared to available reconstructions, although the degree of underestimation varies strongly between the simulations. The simulations with the highest Arctic SAT anomalies tend to match the proxy dataset in its current form better. The ensemble shows some agreement with reconstructions of sea ice, particularly with regard to seasonal sea ice. Large uncertainties limit the confidence that can be placed in the findings and the compatibility of the different proxy datasets. We show that while reducing uncertainties in the reconstructions could decrease the SAT data–model discord substantially, further improvements are likely to be found in enhanced boundary conditions or model physics. Lastly, we compare the Arctic warming in the mPWP to projections of future Arctic warming and find that the PlioMIP2 ensemble simulates greater Arctic amplification than CMIP5 future climate simulations and an increase instead of a decrease in Atlantic Meridional Overturning Circulation (AMOC) strength compared to pre-industrial period. The results highlight the importance of slow feedbacks in equilibrium climate simulations, and that caution must be taken when using simulations of the mPWP as an analogue for future climate change.

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ∼400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.7 and 5.2 ∘C relative to the pre-industrial era with a multi-model mean value of 3.2 ∘C. Annual mean total precipitation rates increase by 7 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases by 4.3 ∘C over land and 2.8 ∘C over the oceans. There is a clear pattern of polar amplification with warming polewards of 60∘ N and 60∘ S exceeding the global mean warming by a factor of 2.3. In the Atlantic and Pacific oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. There is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (equilibrium climate sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble Earth system response to a doubling of CO2 (including ice sheet feedbacks) is 67 % greater than ECS; this is larger than the increase of 47 % obtained from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea surface temperatures are used to assess model estimates of ECS and give an ECS range of 2.6–4.8 ∘C. This result is in general accord with the ECS range presented by previous Intergovernmental Panel on Climate Change (IPCC) Assessment Reports.

A range of future climate scenarios are projected for high atmospheric CO2 concentrations, given uncertainties over future human actions as well as potential environmental and climatic feedbacks. The geological record offers an opportunity to understand climate system response to a range of forcings and feedbacks which operate over multiple temporal and spatial scales. Here, we examine a single interglacial during the late Pliocene (KM5c, ca. 3.205±0.01 Ma) when atmospheric CO2 exceeded pre-industrial concentrations, but were similar to today and to the lowest emission scenarios for this century. As orbital forcing and continental configurations were almost identical to today, we are able to focus on equilibrium climate system response to modern and near-future CO2. Using proxy data from 32 sites, we demonstrate that global mean sea-surface temperatures were warmer than pre-industrial values, by ∼2.3 ∘C for the combined proxy data (foraminifera Mg∕Ca and alkenones), or by ∼3.2–3.4 ∘C (alkenones only). Compared to the pre-industrial period, reduced meridional gradients and enhanced warming in the North Atlantic are consistently reconstructed. There is broad agreement between data and models at the global scale, with regional differences reflecting ocean circulation and/or proxy signals. An uneven distribution of proxy data in time and space does, however, add uncertainty to our anomaly calculations. The reconstructed global mean sea-surface temperature anomaly for KM5c is warmer than all but three of the PlioMIP2 model outputs, and the reconstructed North Atlantic data tend to align with the warmest KM5c model values. Our results demonstrate that even under low-CO2 emission scenarios, surface ocean warming may be expected to exceed model projections and will be accentuated in the higher latitudes.

The Greenland ice sheet contributes increasingly to global sea level rise. Its history during past warm intervals is a valuable reference for future sea level projections. We present ice sheet simulations for the Eemian interglacial period (∼130 000 to 115 000 years ago), a period with warmer-than-present summer climate over Greenland. The evolution of the Eemian Greenland ice sheet is simulated with a 3-D higher-order ice sheet model, forced with a surface mass balance derived from regional climate simulations. Sensitivity experiments with various surface mass balances, basal friction, and ice flow approximations are discussed. The surface mass balance forcing is identified as the controlling factor setting the minimum in Eemian ice volume, emphasizing the importance of a reliable surface mass balance model. Furthermore, the results indicate that the surface mass balance forcing is more important than the representation of ice flow for simulating the large-scale ice sheet evolution. This implies that modeling of the future contribution of the Greenland ice sheet to sea level rise highly depends on an accurate surface mass balance.

An equilibrium simulation of Marine Isotope Stage 3 (MIS3) climate with boundary conditions characteristic of Greenland Interstadial 8 (GI-8; 38 kyr BP) is carried out with the Norwegian Earth System Model (NorESM). A computationally efficient configuration of the model enables long integrations at relatively high resolution, with the simulations reaching a quasi-equilibrium state after 2500 years. We assess the characteristics of the simulated large-scale atmosphere and ocean circulation, precipitation, ocean hydrography, sea ice distribution, and internal variability. The simulated MIS3 interstadial near-surface air temperature is 2.9 ∘C cooler than the pre-industrial (PI). The Atlantic meridional overturning circulation (AMOC) is deeper and intensified by ∼13 %. There is a decrease in the volume of Antarctic Bottom Water (AABW) reaching the Atlantic. At the same time, there is an increase in ventilation of the Southern Ocean, associated with a significant expansion of Antarctic sea ice and concomitant intensified brine rejection, invigorating ocean convection. In the central Arctic, sea ice is ∼2 m thicker, with an expansion of sea ice in the Nordic Seas during winter. Attempts at triggering a non-linear transition to a cold stadial climate state, by varying atmospheric CO2 concentrations and Laurentide Ice Sheet height, suggest that the simulated MIS3 interstadial state in the NorESM is relatively stable, thus underscoring the role of model dependency, and questioning the existence of unforced abrupt transitions in Greenland climate in the absence of interactive ice sheet–meltwater dynamics.

Understanding the behavior of the Greenland ice sheet in a warmer climate, and particularly its surface mass balance (SMB), is important for assessing Greenland's potential contribution to future sea level rise. The Eemian interglacial period, the most recent warmer-than-present period in Earth's history approximately 125 000 years ago, provides an analogue for a warm summer climate over Greenland. The Eemian is characterized by a positive Northern Hemisphere summer insolation anomaly, which complicates Eemian SMB calculations based on positive degree day estimates. In this study, we use Eemian global and regional climate simulations in combination with three types of SMB models – a simple positive degree day, an intermediate complexity, and a full surface energy balance model – to evaluate the importance of regional climate and model complexity for estimates of Greenland's SMB. We find that all SMB models perform well under the relatively cool pre-industrial and late Eemian. For the relatively warm early Eemian, the differences between SMB models are large, which is associated with whether insolation is included in the respective models. For all simulated time slices, there is a systematic difference between globally and regionally forced SMB models, due to the different representation of the regional climate over Greenland. We conclude that both the resolution of the simulated climate as well as the method used to estimate the SMB are important for an accurate simulation of Greenland's SMB. Whether model resolution or the SMB method is most important depends on the climate state and in particular the prevailing insolation pattern. We suggest that future Eemian climate model intercomparison studies should include SMB estimates and a scheme to capture SMB uncertainties.

Rapid retreat of Greenland's marine-terminating glaciers coincides with regional warming trends, which have broadly been used to explain these rapid changes. However, outlet glaciers within similar climate regimes experience widely contrasting retreat patterns, suggesting that the local fjord geometry could be an important additional factor. To assess the relative role of climate and fjord geometry, we use the retreat history of Jakobshavn Isbræ, West Greenland, since the Little Ice Age (LIA) maximum in 1850 as a baseline for the parameterization of a depth- and width-integrated ice flow model. The impact of fjord geometry is isolated by using a linearly increasing climate forcing since the LIA and testing a range of simplified geometries.We find that the total length of retreat is determined by external factors – such as hydrofracturing, submarine melt and buttressing by sea ice – whereas the retreat pattern is governed by the fjord geometry. Narrow and shallow areas provide pinning points and cause delayed but rapid retreat without additional climate warming, after decades of grounding line stability. We suggest that these geometric pinning points may be used to locate potential sites for moraine formation and to predict the long-term response of the glacier. As a consequence, to assess the impact of climate on the retreat history of a glacier, each system has to be analyzed with knowledge of its historic retreat and the local fjord geometry.

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.

Understanding of long-term dynamics of glaciers and ice caps is vital to assess their recent and future changes, yet few long-term reconstructions using ice flow models exist. Here we present simulations of the maritime Hardangerjøkulen ice cap in Norway from the mid-Holocene through the Little Ice Age (LIA) to the present day, using a numerical ice flow model combined with glacier and climate reconstructions. In our simulation, under a linear climate forcing, we find that Hardangerjøkulen grows from ice-free conditions in the mid-Holocene to its maximum extent during the LIA in a nonlinear, spatially asynchronous fashion. During its fastest stage of growth (2300–1300 BP), the ice cap triples its volume in less than 1000 years. The modeled ice cap extent and outlet glacier length changes from the LIA until today agree well with available observations. Volume and area for Hardangerjøkulen and several of its outlet glaciers vary out-of-phase for several centuries during the Holocene. This volume–area disequilibrium varies in time and from one outlet glacier to the next, illustrating that linear relations between ice extent, volume and glacier proxy records, as generally used in paleoclimatic reconstructions, have only limited validity. We also show that the present-day ice cap is highly sensitive to surface mass balance changes and that the effect of the ice cap hypsometry on the mass balance–altitude feedback is essential to this sensitivity. A mass balance shift by +0.5 m w.e. relative to the mass balance from the last decades almost doubles ice volume, while a decrease of 0.2 m w.e. or more induces a strong mass balance–altitude feedback and makes Hardangerjøkulen disappear entirely. Furthermore, once disappeared, an additional +0.1 m w.e. relative to the present mass balance is needed to regrow the ice cap to its present-day extent. We expect that other ice caps with comparable geometry in, for example, Norway, Iceland, Patagonia and peripheral Greenland may behave similarly, making them particularly vulnerable to climate change.

The last interglacial (LIG, ~130–116 ka, ka = 1000 yr ago) is characterized by high-latitude warming and is therefore often considered as a possible analogue for future warming. However, in contrast to predicted future greenhouse warming, the LIG climate is largely governed by variations in insolation. Greenhouse gas (GHG) concentrations were relatively stable and similar to pre-industrial values, with the exception of the early LIG when, on average, GHGs were slightly lower. We performed six time-slice simulations with the low-resolution version of the Norwegian Earth System Model covering the LIG. In four simulations only the orbital forcing was changed. In two other simulations, representing the early LIG, additionally the GHG forcing was reduced. With these simulations we investigate (1) the different effects of GHG versus insolation forcing on the temperatures during the LIG; (2) whether reduced GHGs can explain the low temperatures reconstructed for the North Atlantic; and (3) the timing of the observed LIG peak warmth. Our simulations show that the insolation forcing results in seasonal and hemispheric differences in temperature. In contrast, a reduction in the GHG forcing causes a global and seasonal-independent cooling. Furthermore, we compare modelled temperatures with proxy-based LIG sea-surface temperatures along a transect in the North Atlantic. The modelled North Atlantic summer sea-surface temperatures capture the general trend of the reconstructed summer temperatures, with low values in the early LIG, a peak around 125 ka, and a steady decrease towards the end of the LIG. Simulations with reduced GHG forcing improve the model–data fit as they show lower temperatures in the early LIG. Furthermore we show that the timing of maximum summer and winter surface temperatures is in line with the local summer and winter insolation maximum at most latitudes. Two regions where the maximum local insolation and temperature do not occur at the same time are Antarctica and the Southern Ocean. The austral summer insolation has a late maximum at ~115 ka. In contrast the austral summer temperatures in Antarctica show maxima at both ~130 ka and ~115 ka, and the Southern Ocean temperatures peak only at ~130 ka. This is probably due to the integrating effect of the ocean, storing heat from other seasons and resulting in relatively warm austral summer temperatures. Reducing the GHG concentrations in the early LIG (125 and 130 ka) results in a similar timing of peak warmth, except over Antarctica. There, the lower austral summer temperatures at 130 ka shift the maximum warmth to a single peak at 115 ka.

In the Pliocene Model Intercomparison Project (PlioMIP), eight state-of-the-art coupled climate models have simulated the mid-Pliocene warm period (mPWP, 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), northward ocean heat transport and ocean stratification simulated with these models. None of the models participating in PlioMIP simulates a strong mid-Pliocene AMOC as suggested by earlier proxy studies. Rather, there is no consistent increase in AMOC maximum among the PlioMIP models. The only consistent change in AMOC is a shoaling of the overturning cell in the Atlantic, and a reduced influence of North Atlantic Deep Water (NADW) at depth in the basin. Furthermore, the simulated mid-Pliocene Atlantic northward heat transport is similar to the pre-industrial. These simulations demonstrate that the reconstructed high-latitude mid-Pliocene warming can not be explained as a direct response to an intensification of AMOC and concomitant increase in northward ocean heat transport by the Atlantic.

The last interglaciation (~130 to 116 ka) is a time period with a strong astronomically induced seasonal forcing of insolation compared to the present. Proxy records indicate a significantly different climate to that of the modern, in particular Arctic summer warming and higher eustatic sea level. Because the forcings are relatively well constrained, it provides an opportunity to test numerical models which are used for future climate prediction. In this paper we compile a set of climate model simulations of the early last interglaciation (130 to 125 ka), encompassing a range of model complexities. We compare the simulations to each other and to a recently published compilation of last interglacial temperature estimates. We show that the annual mean response of the models is rather small, with no clear signal in many regions. However, the seasonal response is more robust, and there is significant agreement amongst models as to the regions of warming vs cooling. However, the quantitative agreement of the model simulations with data is poor, with the models in general underestimating the magnitude of response seen in the proxies. Taking possible seasonal biases in the proxies into account improves the agreement, but only marginally. However, a lack of uncertainty estimates in the data does not allow us to draw firm conclusions. Instead, this paper points to several ways in which both modelling and data could be improved, to allow a more robust model–data comparison.

Using simulated climate data from the comprehensive coupled climate model IPSL CM4, we simulate the Greenland ice sheet (GrIS) during the Eemian interglaciation with the three-dimensional ice sheet model SICOPOLIS. The Eemian is a period 126 000 yr before present (126 ka) with Arctic temperatures comparable to projections for the end of this century. In our simulation, the northeastern part of the GrIS is unstable and retreats significantly, despite moderate melt rates. This result is found to be robust to perturbations within a wide parameter space of key parameters of the ice sheet model, the choice of initial ice temperature, and has been reproduced with climate forcing from a second coupled climate model, the CCSM3. It is shown that the northeast GrIS is the most vulnerable. Even a small increase in melt removes many years of ice accumulation, giving a large mass imbalance and triggering the strong ice-elevation feedback. Unlike the south and west, melting in the northeast is not compensated by high accumulation. The analogy with modern warming suggests that in coming decades, positive feedbacks could increase the rate of mass loss of the northeastern GrIS, exceeding the recent observed thinning rates in the south.

The Last Glacial Maximum (LGM; 21 000 yr before present) was a period of low atmospheric greenhouse gas concentrations, when vast ice sheets covered large parts of North America and Europe. Paleoclimate reconstructions and modeling studies suggest that the atmospheric circulation was substantially altered compared to today, both in terms of its mean state and its variability. Here we present a suite of coupled model simulations designed to investigate both the separate and combined influences of the main LGM boundary condition changes (greenhouse gases, ice sheet topography and ice sheet albedo) on the mean state and variability of the atmospheric circulation as represented by sea level pressure (SLP) and 200-hPa zonal wind in the North Atlantic sector. We find that ice sheet topography accounts for most of the simulated changes during the LGM. Greenhouse gases and ice sheet albedo affect the SLP gradient in the North Atlantic, but the overall placement of high and low pressure centers is controlled by topography. Additional analysis shows that North Atlantic sea surface temperatures and sea ice edge position do not substantially influence the pattern of the climatological-mean SLP field, SLP variability or the position of the North Atlantic jet in the LGM.

Following the Early Eocene climatic optimum (EECO, ~55–50 Ma), climate deteriorated and gradually changed the earth from a greenhouse into an icehouse, with major cooling events at the Eocene-Oligocene boundary (∼34 Ma) and the Middle Miocene (∼15 Ma). It is believed that the opening of the Drake Passage had a marked impact on the cooling at the Eocene-Oligocene boundary. Based on an Early Eocene simulation, we study the sensitivity of climate and ocean circulation to tectonic events such as the closing of the West Siberian Seaway, the deepening of the Arctic-Atlantic Seaway, the opening of the Drake Passage, and the constriction of the Tethys and Central American seaways. The opening of the Drake Passage, together with the closing of the West Siberian Seaway and the deepening of the Arctic-Atlantic Seaway, weakened the Southern Ocean Deep Water (SODW) dominated ocean circulation and led to a weak cooling at high latitudes, thus contributing to the observed Early Cenozoic cooling. However, the later constriction of the Tethys and Central American Seaways is shown to give a strong cooling at southern high latitudes. This cooling was related to the transition of ocean circulation from a SODW-dominated mode to the modern-like ocean circulation dominated by North Atlantic Deep Water (NADW).

We simulate the last glacial inception, 115 000 years ago, with a three dimensional thermomechanical ice sheet model of the Northern Hemisphere, forced by a comprehensive coupled climate model. High oceanic heat transport into the Nordic Seas prevents large scale ice growth over Scandinavia. Glacial inception in the region starts on the highest mountains in the south when sea surface temperatures in the Nordic Seas are reduced by at least 3 °C. Ice growth in Northern Scandinavia requires a cooling by at least 4 °C. This is in good agreement with marine proxy data from the Nordic Seas and North Atlantic as well as available terrestrial data. This study thus provides a physical understanding and revised timing of the first glacier advance over Scandinavia.

Using four different climate models, we investigate sea level pressure variability in the extratropical North Atlantic in the preindustrial climate (1750 AD) and at the Last Glacial Maximum (LGM, 21 kyrs before present) in order to understand how changes in atmospheric circulation can affect signals recorded in climate proxies. In general, the models exhibit a significant reduction in interannual variance of sea level pressure at the LGM compared to pre-industrial simulations and this reduction is concentrated in winter. For the preindustrial climate, all models feature a similar leading mode of sea level pressure variability that resembles the leading mode of variability in the instrumental record: the North Atlantic Oscillation (NAO). In contrast, the leading mode of sea level pressure variability at the LGM is model dependent, but in each model different from that in the preindustrial climate. In each model, the leading (NAO-like) mode of variability explains a smaller fraction of the variance and also less absolute variance at the LGM than in the preindustrial climate. The models show that the relationship between atmospheric variability and surface climate (temperature and precipitation) variability change in different climates. Results are model-specific, but indicate that proxy signals at the LGM may be misinterpreted if changes in the spatial pattern and seasonality of surface climate variability are not taken into account.