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

Future sea level contribution from Antarcticainferred from CMIP5 model forcing and itsdependence on precipitation ansatz

Various observational estimates indicate growing mass loss at Antarctica's margins but also heavier precipitation across the continent. In the future, heavier precipitation fallen on Antarctica will counteract any stronger iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we use from nine CMIP5 models future projections, ranging from strong mitigation efforts to business-as-usual, to run an ensemble of ice-sheet simulations. We test, how the precipitation boundary condition determines Antarctica's sea-level contribution. The spatial and temporal varying climate forcings drive ice-sheet simulations. Hence, our ensemble inherits all spatial and temporal climate patterns, which is in contrast to a spatial mean forcing. Regardless of the applied boundary condition and forcing, some areas will lose ice in the future, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea. In general the simulated ice-sheet thickness grows in a broad marginal strip, where incoming storms deliver topographically controlled precipitation. This strip shows the largest ice thickness differences between the applied precipitation boundary conditions too. On average Antarctica's ice mass shrinks for all future scenarios if the precipitation is scaled by the spatial temperature anomalies coming from the CMIP5 models. In this approach, we use the relative precipitation increment per degree warming as invariant scaling constant. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies of the CMIP5 models directly. Here, the scaling factors show a distinct spatial pattern across Antarctica. Furthermore, the diagnosed mean scaling across all considered climate forcings is larger than the values deduced from ice cores. In general, the scaling is higher across the East Antarctic Ice Sheet, lower across the West Antarctic Ice Sheet, and lowest around the Siple Coast. The latter is located on the east side of the Ross Ice Shelf

DMI Report 21-17 Including a dynamic Greenland Ice Sheet in the EC-Earth global climate model

Recent observations have indicated rapidly increasing mass loss from the Greenland Ice Sheet. To explore the interactions and feedbacks of the ice sheets in the climate system, it is important to develop coupled climate-ice sheet models. The integration of an ice sheet model in a global model is challenging, and, currently, relatively few climate models include a two-way coupling to a dynamical ice sheet model. In this work package, we have continued developing the coupled ice sheet-climate model system comprising the global climate model EC-Earth and the Parallel Ice Sheet Model (PISM) for Greenland. The new model system, EC-Earth3-GrIS, is upgraded to include the recent model versions, EC-Earth3 and PISM version 1.2. In addition, a new module has been developed to handle the exchange of information between the ice sheet model and EC-Earth using the OASIS3- MCT software interface. The new module reads output from the ice sheet model and exchanges the fields with the relevant EC-Earth components. The ice sheet mask and topography are provided to the atmosphere and land surface components. The heat and freshwater fluxes from basal melt and ice discharge are provided to the ocean module via the runoff-mapper that routes surface runoff into the ocean. The new module also prepares the forcing fields for the ice sheet model, i.e., subsurface temperature and surface mass balance. These fields are calculated in EC- Earth3 using a land ice surface parameterization, developed explicitly for the Greenland ice sheet. The parameterization contains a responsive snow and ice albedo scheme and includes land ice characteristics in the calculation of heat and energy transfer at the surface. Experiments with and without the land ice surface parameterization have been carried out for preindustrial and present-day conditions to assess the influence of the surface parameterization on the calculated surface mass balance. The results show that the ice sheet responds stronger and more realistically to forcing changes when the new surface parameterization is used. Besides the model development, the results from experiments with the first model version, EC- Earth-PISM, have been analyzed. These results stress that a decent surface scheme with a responsive snow albedo scheme is necessary for investigating mass balance changes of the Greenland Ice Sheet. Overall, our results indicate that the feedbacks induced by the interactive ice sheet have a significant influence on Arctic climate change under warming conditions. In warm scenarios where the CO2 level is raised to four times the preindustrial level, the coupled model has a colder Arctic surface, a fresher ocean, and more sea-ice in winter

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

Brief communication: A submarine wall protecting the Amundsen Sea intensifies melting of neighboring ice shelves

Disintegration of ice shelves in the Amundsen Sea, in front of the West Antarctic Ice Sheet, has the potential to cause sea level rise by inducing an acceleration of ice discharge from upstream grounded ice. Moore et al. (2018) proposed that using a submarine wall to block the penetration of warm water into the subsurface cavities of these ice shelves could reduce this risk. We use a global sea ice–ocean model to show that a wall shielding the Amundsen Sea below 350 m depth successfully suppresses the inflow of warm water and reduces ice shelf melting. However, these warm water masses get redirected towards neighboring ice shelves, which reduces the net effectiveness of the wall. The ice loss is reduced by 10 %, integrated over the entire Antarctic continent

Analysis of the surface mass balance for deglacial climate simulations

A realistic simulation of the surface mass balance (SMB) is essential for simulating past and future ice-sheet changes. As most state-of-the-art Earth system models (ESMs) are not capable of realistically representing processes determining the SMB, most studies of the SMB are limited to observations and regional climate models and cover the last century and near future only. Using transient simulations with the Max Planck Institute ESM in combination with an energy balance model (EBM), we extend previous research and study changes in the SMB and equilibrium line altitude (ELA) for the Northern Hemisphere ice sheets throughout the last deglaciation. The EBM is used to calculate and downscale the SMB onto a higher spatial resolution than the native ESM grid and allows for the resolution of SMB variations due to topographic gradients not resolved by the ESM. An evaluation for historical climate conditions (1980–2010) shows that derived SMBs compare well with SMBs from regional modeling. Throughout the deglaciation, changes in insolation dominate the Greenland SMB. The increase in insolation and associated warming early in the deglaciation result in an ELA and SMB increase. The SMB increase is caused by compensating effects of melt and accumulation: the warming of the atmosphere leads to an increase in melt at low elevations along the ice-sheet margins, while it results in an increase in accumulation at higher levels as a warmer atmosphere precipitates more. After 13 ka, the increase in melt begins to dominate, and the SMB decreases. The decline in Northern Hemisphere summer insolation after 9 ka leads to an increasing SMB and decreasing ELA. Superimposed on these long-term changes are centennial-scale episodes of abrupt SMB and ELA decreases related to slowdowns of the Atlantic meridional overturning circulation (AMOC) that lead to a cooling over most of the Northern Hemisphere

On the reduced sensitivity of the Atlantic overturning to Greenland ice sheet melting in projections: a multi-model assessment

Large uncertainties exist concerning the impact of Greenland ice sheet melting on the Atlantic meridional overturning circulation (AMOC) in the future, partly due to different sensitivity of the AMOC to freshwater input in the North Atlantic among climate models. Here we analyse five projections from different coupled ocean–atmosphere models with an additional 0.1 Sv (1 Sv = 10 6 m3/s) of freshwater released around Greenland between 2050 and 2089. We find on average a further weakening of the AMOC at 26°N of 1.1 ± 0.6 Sv representing a 27 ± 14% supplementary weakening in 2080–2089, as compared to the weakening relative to 2006–2015 due to the effect of the external forcing only. This weakening is lower than what has been found with the same ensemble of models in an identical experimen - tal set-up but under recent historical climate conditions. This lower sensitivity in a warmer world is explained by two main factors. First, a tendency of decoupling is detected between the surface and the deep ocean caused by an increased thermal stratification in the North Atlantic under the effect of global warming. This induces a shoaling of ocean deep ventilation through convection hence ventilating only intermediate levels. The second important effect concerns the so-called Canary Current freshwater leakage; a process by which additionally released fresh water in the North Atlantic leaks along the Canary Current and escapes the convection zones towards the subtropical area. This leakage is increasing in a warming climate, which is a consequence of decreasing gyres asymmetry due to changes in Ekman rumping. We suggest that these modifications are related with the northward shift of the jet stream in a warmer world. For these two reasons the AMOC is less susceptible to freshwater perturbations (near the deep water formation sides) in the North Atlantic as compared to the recent historical climate conditions. Finally, we propose a bilinear model that accounts for the two former processes to give a conceptual explanation about the decreasing AMOC sensitivity due to freshwater input. Within the limit of this bilinear model, we find that 62 ± 8% of the reduction in sensitivity is related with the changes in gyre asymmetry and freshwater leakage and 38 ± 8% is due to the reduction in deep ocean ventilation associated with the increased stratification in the North Atlantic

Reaching the 1.5 degree limit: what does it mean for West Antarctica and the global mean sea level?

What are the benefits of limiting the global warming to 1.5 degree with respect to pre-industrial conditions for the vulnerable region of West Antarctica which might be prone to positive feedback mechanisms between ocean circulation, melting of shelf ice and instabilities of the ice sheet? There are indications that West Antarctic ice sheet instabilities have occurred in the Last Interglacial around 125.000 years ago. At that time the polar surface temperature was about 2K warmer than today. The question under which circumstances a tipping point may be reached and if this may happen again is therefore highly relevant, especially since a disintegration of the West Antarctic ice sheet could cause a global sea level rise between 3 and 5 m. Here we address this question with variable resolution, global coupled ice sheet - shelf ice - ocean - atmosphere multi-century simulations. With our innovative ocean modelling approach in the Finite Element Sea-ice Ocean Model FESOM it is possible to refine the ocean resolution to up to 3 km in the Amundsen Sea and 10 km around the whole Antarctica while keeping it relatively coarse in the order of a couple of hundred km in dynamically not very active regions such as the subtropical regions. This means that we can simulate the feedback between ocean and ice in the relevant regions highly resolved given that the ice sheet model runs at a resolution of 5 to 10 km. Three different emission scenarios are applied up to 2100, two of them limiting the global mean temperature increase to 1.5 ◦ C and 2 ◦ C respectively and one of them assuming business-as-usual conditions (IPCC SRES RCP8.5 scenario). The simulations are extended to 2400 with the greenhouse gas and aerosol concentrations kept constant at 2100 levels, respectively, to be able to simulate the long-term implications of different global warming levels