Coupled ice sheet–climate modeling under glacial and pre-industrial boundary conditions

In the standard Paleoclimate Modelling Intercomparison Project (PMIP) experiments, the Last Glacial Maximum (LGM) is modeled in quasi-equilibrium with atmosphere–ocean–vegetation general circulation models (AOVGCMs) with prescribed ice sheets. This can lead to inconsistencies between the modeled climate and ice sheets. One way to avoid this problem would be to model the ice sheets explicitly. Here, we present the first results from coupled ice sheet–climate simulations for the pre-industrial times and the LGM. Our setup consists of the AOVGCM ECHAM5/MPIOM/LPJ bidirectionally coupled with the Parallel Ice Sheet Model (PISM) covering the Northern Hemisphere. The results of the pre-industrial and LGM simulations agree reasonably well with reconstructions and observations. This shows that the model system adequately represents large, non-linear climate perturbations. A large part of the drainage of the ice sheets occurs in ice streams. Most modeled ice stream systems show recurring surges as internal oscillations. The Hudson Strait Ice Stream surges with an ice volume equivalent to about 5 m sea level and a recurrence interval of about 7000 yr. This is in agreement with basic expectations for Heinrich events. Under LGM boundary conditions, different ice sheet configurations imply different locations of deep water formation

Coupled ice sheet–climate modeling under glacial and pre-industrial boundary conditions

In the standard Paleoclimate Modelling Intercomparison Project (PMIP) experiments, the Last Glacial Maximum (LGM) is modeled in quasi-equilibrium with atmosphere–ocean–vegetation general circulation models (AOVGCMs) with prescribed ice sheets. This can lead to inconsistencies between the modeled climate and ice sheets. One way to avoid this problem would be to model the ice sheets explicitly. Here, we present the first results from coupled ice sheet–climate simulations for the pre-industrial times and the LGM. Our setup consists of the AOVGCM ECHAM5/MPIOM/LPJ bidirectionally coupled with the Parallel Ice Sheet Model (PISM) covering the Northern Hemisphere. The results of the pre-industrial and LGM simulations agree reasonably well with reconstructions and observations. This shows that the model system adequately represents large, non-linear climate perturbations. A large part of the drainage of the ice sheets occurs in ice streams. Most modeled ice stream systems show recurring surges as internal oscillations. The Hudson Strait Ice Stream surges with an ice volume equivalent to about 5 m sea level and a recurrence interval of about 7000 yr. This is in agreement with basic expectations for Heinrich events. Under LGM boundary conditions, different ice sheet configurations imply different locations of deep water formation

Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison

Earlier large-scale Greenland ice sheet sea-level projections (e.g. those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of this initMIP-Greenland intercomparison exercise is to compare, evaluate, and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project Phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of (1) the initial present-day state of the ice sheet and (2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly); they should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions, and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises

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

Modeling the spreading of glacial meltwater from the Amundsen and Bellingshausen Seas

It has been suggested that an increased melting of continental ice in the Amundsen Sea (AS) and Bellingshausen Sea (BS) is a likely source of the observed freshening of Ross Sea (RS) water. To test this hypothesis, we simulate the spreading of glacial meltwater using the Finite Element Sea Ice/Ice Shelf/Ocean Model. Based on the spatial distribution of simulated passive tracers, most of the basal meltwater from AS ice shelves flows toward the RS with more than half of the melt originating from the Getz Ice Shelf. Further, the model results show that a slight increase of the basal mass loss can substantially intensify the transport of meltwater into the RS due to a strengthening of the melt-driven shelf circulation and the westward flowing coastal current. This supports the idea that the basal melting of AS and BS ice shelves is one of the main sources for the RS freshening

Future sea level contribution from Antarctica inferred from CMIP5 model forcing and its dependence on precipitation ansatz

Abstract. Various observational estimates indicate growing mass loss at Antarctica's margins as well as heavier precipitation across the continent. Simulated future projections reveal that heavier precipitation, falling on Antarctica, may counteract amplified iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we test how the ansatz (implementation in a mathematical framework) of the precipitation boundary condition shapes Antarctica's sea level contribution in an ensemble of ice sheet simulations. We test two precipitation conditions: we either apply the precipitation anomalies from CMIP5 models directly or scale the precipitation by the air temperature anomalies from the CMIP5 models. In the scaling approach, it is common to use a relative precipitation increment per degree warming as an invariant scaling constant. We use future climate projections from nine CMIP5 models, ranging from strong mitigation efforts to business-as-usual scenarios, to perform simulations from 1850 to 5000. We take advantage of individual climate projections by exploiting their full temporal and spatial structure. The CMIP5 projections beyond 2100 are prolonged with reiterated forcing that includes decadal variability; hence, our study may underestimate ice loss after 2100. In contrast to various former studies that apply an evolving temporal forcing that is spatially averaged across the entire Antarctic Ice Sheet, our simulations consider the spatial structure in the forcing stemming from various climate patterns. This fundamental difference reproduces regions of decreasing precipitation despite general warming. Regardless of the boundary and forcing conditions applied, our ensemble study suggests that some areas, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea, will lose ice in the future. In general, the simulated ice sheet thickness grows along the coast, where incoming storms deliver topographically controlled precipitation. In this region, the ice thickness differences are largest between the applied precipitation methods. On average, Antarctica shrinks for all future scenarios if the air temperature anomalies scale the precipitation. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies directly. The analysis reveals that the mean scaling inferred from climate models is larger than the commonly used values deduced from ice cores; moreover, it varies spatially: the highest scaling is across the East Antarctic Ice Sheet, and the lowest scaling is around the Siple Coast, east of the Ross Ice Shelf. The discrepancies in response to both precipitation ansatzes illustrate the principal uncertainty in projections of Antarctica's sea level contribution

The role of an interactive Greenland ice sheet in the coupled climate-ice sheet model EC-Earth-PISM

AbstractIce sheet processes are often simplified in global climate models as changes in ice sheets have been assumed to occur over long time scales compared to ocean and atmospheric changes. However, numerous observations show an increasing rate of mass loss from the Greenland Ice Sheet and call for comprehensive process-based models to explore its role in climate change. Here, we present a new model system, EC-Earth-PISM, that includes an interactive Greenland Ice Sheet. The model is based on the EC-Earth v2.3 global climate model in which ice sheet surface processes are introduced. This model interacts with the Parallel Ice Sheet Model (PISM) without anomaly or flux corrections. Under pre-industrial climate conditions, the modeled climate and ice sheet are stable while keeping a realistic interannual variability. In model simulations forced into a warmer climate of four times the pre-industrial CO2 concentration, the total surface mass balance decreases and the ice sheet loses mass at a rate of about 500 Gt/year. In the climate warming experiments, the resulting freshwater flux from the Greenland Ice Sheet increases 55% more in the experiments with the interactive ice sheet and the climate response is significantly different: the Arctic near-surface air temperature is lower, substantially more winter sea ice covers the northern hemisphere, and the ocean circulation is weaker. Our results indicate that the melt-albedo feedback plays a key role for the response of the ice sheet and its influence on the changing climate in the Arctic. This emphasizes the importance of including interactive ice sheets in climate change projections.</jats:p