Investigating the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle

A 2.5-dimensional climate model of intermediate complexity, CLIMBER-2, fully coupled with the GREMLINS 3-D thermo-mechanical ice sheet model is used to simulate the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle and to investigate the ice sheets responses to both insolation and atmospheric CO<sub>2</sub> concentration. This model reproduces the main phases of advance and retreat of Northern Hemisphere ice sheets during the last glacial cycle, although the amplitude of these variations is less pronounced than those based on sea level reconstructions. At the last glacial maximum, the simulated ice volume is 52.5×10<sup>15</sup> m<sup>3</sup> and the spatial distribution of both the American and Eurasian ice complexes is in reasonable agreement with observations, with the exception of the marine parts of these former ice sheets. <br> A set of sensitivity studies has also been performed to assess the sensitivity of the Northern Hemisphere ice sheets to both insolation and atmospheric CO<sub>2</sub>. Our results suggest that the decrease of summer insolation is the main factor responsible for the early build up of the North American ice sheet around 120 kyr BP, in agreement with benthic foraminifera δ<sup>18</sup>O signals. In contrast, low insolation and low atmospheric CO<sub>2</sub> concentration are both necessary to trigger a long-lasting glaciation over Eurasia

Investigating similarities and differences of the penultimate and last glacial terminations with a coupled ice sheet–climate model

Glacial terminations are marked by a re-organisation of the different components of the climate system. In particular, rapid ice sheet disintegration leads to multiple complex feedback loops that are still poorly understood. To further investigate this aspect, we use here a fully coupled Northern Hemisphere ice sheet–climate model to perform numerical experiments of the last two glacial terminations. We show that even if the first-order climate trajectory is similar for the two terminations, the difference in terms of solar insolation leads to important changes for the ice sheet–climate system. Warmer temperatures during the penultimate termination are compatible with higher sea level during the last interglacial period with respect to the Holocene. We simulate a last interglacial Greenland contribution to sea level rise of about 2 m of sea level equivalent. We also simulate warmer subsurface Southern Ocean, compatible with an additional contribution from the Antarctic ice sheet. In addition, even without considering freshwater flux to the ocean resulting from ice sheet melting, the two terminations display different Atlantic overturning circulation sensitivity, this circulation being more prone to collapses during the penultimate termination. Finally, with additional sensitivity experiments we show that, for the two terminations, the Northern Hemisphere insolation is the main driver for the ice sheet retreat even if vegetation changes have also to be taken into account to simulate the full deglaciation. Conversely, even though it impacts the temperature, greenhouse gas concentration change alone does not explain the amplitude of ice sheet retreat and only modulates its timing.</p

Relative importance of the mechanisms triggering the Eurasian ice sheet deglaciation in the GRISLI2.0 ice sheet model

The last deglaciation (21 to 8 ka) of the Eurasian ice sheet (EIS) is thought to have been responsible for a sea level rise of about 20 m. While many studies have examined the timing and rate of the EIS retreat during this period, many questions remain about the key processes that triggered the EIS deglaciation 21 kyr ago. Due to its large marine-based parts in the Barents–Kara (BKIS) and British Isles sectors, the BKIS is often considered to be a potential analogue of the current West Antarctic ice sheet (WAIS). Identifying the mechanisms that drove the EIS evolution might provide a better understanding of the processes at play in the West Antarctic destabilization. To investigate the relative impact of key drivers on the EIS destabilization, we used the three-dimensional ice sheet model GRISLI (GRenoble Ice Shelf and Land Ice) (version 2.0) forced by climatic fields from five Paleoclimate Modelling Intercomparison Project phases 3 and 4 (PMIP3, PMIP4) Last Glacial Maximum (LGM) simulations. In this study, we performed sensitivity experiments to test the response of the simulated Eurasian ice sheets to surface climate, oceanic temperatures (and thus basal melting under floating ice tongues), and sea level perturbations. Our results highlight that the EIS retreat simulated with the GRISLI model is primarily triggered by atmospheric warming. Increased atmospheric temperatures further amplify the sensitivity of the ice sheets to sub-shelf melting. These results contradict those of previous modelling studies mentioning the central role of basal melting on the deglaciation of the marine-based Barents–Kara ice sheet. However, we argue that the differences with previous works are mainly related to differences in the methodology followed to generate the initial LGM ice sheet. Due to the strong sensitivity of EIS to the atmospheric forcing highlighted with the GRISLI model and the limited extent of the confined ice shelves during the LGM, we conclude by questioning the analogy between EIS and the current WAIS. However, because of the expected rise in atmospheric temperatures, the risk of hydrofracturing is increasing and could ultimately put the WAIS in a configuration similar to the past Eurasian ice sheet.</p

Sensitivity of a Greenland ice sheet model to atmospheric forcing fields

International audiencePredicting the climate for the future and how it will impact ice sheet evolution requires coupling ice sheet models with climate models. However, before we attempt to develop a realistic coupled setup, we propose, in this study, to first analyse the impact of a model simulated climate on an ice sheet. We undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary conditions to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyrs of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed one, there are considerable deviations among the ice sheets on regional scales. These deviations can be explained by biases in temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations between the climate models are also due to the differences in the atmospheric general circulation. To account for these differences in the context of coupling ice sheet models with climate models, we conclude that appropriate downscaling methods will be needed. In some cases, systematic corrections of the climatic variables at the interface may be required to obtain realistic results for the Greenland ice sheet (GIS)

Ice sheet model dependency of the simulated Greenland Ice Sheet in the mid-Pliocene

The understanding of the nature and behavior of ice sheets in past warm periods is important for constraining the potential impacts of future climate change. The Pliocene warm period (between 3.264 and 3.025 Ma) saw global temperatures similar to those projected for future climates; nevertheless, Pliocene ice locations and extents are still poorly constrained. We present results from the efforts to simulate mid-Pliocene Greenland Ice Sheets by means of the international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP). We compare the performance of existing numerical ice sheet models in simulating modern control and mid-Pliocene ice sheets with a suite of sensitivity experiments guided by available proxy records. We quantify equilibrated ice sheet volume on Greenland, identifying a potential range in sea level contributions from warm Pliocene scenarios. A series of statistical measures are performed to quantify the confidence of simulations with focus on inter-model and inter-scenario differences. We find that Pliocene Greenland Ice Sheets are less sensitive to differences in ice sheet model configurations and internal physical quantities than to changes in imposed climate forcing. We conclude that Pliocene ice was most likely to be limited to the highest elevations in eastern and southern Greenland as simulated with the highest confidence and by synthesizing available regional proxies; however, the extent of those ice caps needs to be further constrained by using a range of general circulation model (GCM) climate forcings

The future sea-level contribution of the Greenland ice sheet: A multi-model ensemble study of ISMIP6

The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6).We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90-50 and 32-17mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean. © Author(s) 2020

The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6

The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90±50 and 32±17 mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean

The PMIP4 Last Glacial Maximum experiments: preliminary results and comparison with the PMIP3 simulations

The Last Glacial Maximum (LGM, ∼ 21 000 years ago) has been a major focus for evaluating how well state-of-the-art climate models simulate climate changes as large as those expected in the future using paleoclimate reconstructions. A new generation of climate models has been used to generate LGM simulations as part of the Paleoclimate Modelling Intercomparison Project (PMIP) contribution to the Coupled Model Intercomparison Project (CMIP). Here, we provide a preliminary analysis and evaluation of the results of these LGM experiments (PMIP4, most of which are PMIP4-CMIP6) and compare them with the previous generation of simulations (PMIP3, most of which are PMIP3-CMIP5). We show that the global averages of the PMIP4 simulations span a larger range in terms of mean annual surface air temperature and mean annual precipitation compared to the PMIP3-CMIP5 simulations, with some PMIP4 simulations reaching a globally colder and drier state. However, the multi-model global cooling average is similar for the PMIP4 and PMIP3 ensembles, while the multi-model PMIP4 mean annual precipitation average is drier than the PMIP3 one. There are important differences in both atmospheric and oceanic circulations between the two sets of experiments, with the northern and southern jet streams being more poleward and the changes in the Atlantic Meridional Overturning Circulation being less pronounced in the PMIP4-CMIP6 simulations than in the PMIP3-CMIP5 simulations. Changes in simulated precipitation patterns are influenced by both temperature and circulation changes. Differences in simulated climate between individual models remain large. Therefore, although there are differences in the average behaviour across the two ensembles, the new simulation results are not fundamentally different from the PMIP3-CMIP5 results. Evaluation of large-scale climate features, such as land–sea contrast and polar amplification, confirms that the models capture these well and within the uncertainty of the paleoclimate reconstructions. Nevertheless, regional climate changes are less well simulated: the models underestimate extratropical cooling, particularly in winter, and precipitation changes. These results point to the utility of using paleoclimate simulations to understand the mechanisms of climate change and evaluate model performance

Methane and carbon dioxide fluxes and their regional scalability for the European Arctic wetlands during the MAMM project in summer 2012

Airborne and ground-based measurements of methane (CH4), carbon dioxide (CO2) and boundary layer thermodynamics were recorded over the Fennoscandian landscape (67–69.5° N, 20–28° E) in July 2012 as part of the MAMM (Methane and other greenhouse gases in the Arctic: Measurements, process studies and Modelling) field campaign. Employing these airborne measurements and a simple boundary layer box model, net regional-scale (~ 100 km) fluxes were calculated to be 1.2 ± 0.5 mg CH4 h−1 m−2 and −350 ± 143 mg CO2 h−1 m−2. These airborne fluxes were found to be relatively consistent with seasonally averaged surface chamber (1.3 ± 1.0 mg CH4 h−1 m−2) and eddy covariance (1.3 ± 0.3 mg CH4 h−1 m−2 and −309 ± 306 mg CO2 h−1 m−2) flux measurements in the local area. The internal consistency of the aircraft-derived fluxes across a wide swath of Fennoscandia coupled with an excellent statistical comparison with local seasonally averaged ground-based measurements demonstrates the potential scalability of such localised measurements to regional-scale representativeness. Comparisons were also made to longer-term regional CH4 climatologies from the JULES (Joint UK Land Environment Simulator) and HYBRID8 land surface models within the area of the MAMM campaign. The average hourly emission flux output for the summer period (July–August) for the year 2012 was 0.084 mg CH4 h−1 m−2 (minimum 0.0 and maximum 0.21 mg CH4 h−1 m−2) for the JULES model and 0.088 mg CH4 h−1 m−2 (minimum 0.0008 and maximum 1.53 mg CH4 h−1 m−2) for HYBRID8. Based on these observations both models were found to significantly underestimate the CH4 emission flux in this region, which was linked to the under-prediction of the wetland extents generated by the models