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

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

InitMIP-Antarctica:An ice sheet model initialization experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue