A major impact of anthropogenic climate change is the crossing of tipping points, which may have severe consequences such as the complete mass loss of the Greenland ice sheet (GrIS). At present, the GrIS is losing mass at an accelerated rate, largely due to a steep decrease in its surface mass balance (SMB; the balance between snow accumulation and surface ablation from melt and associated runoff). Previous work on the magnitude and nature of a threshold for GrIS complete melt remains controversial. Here, we explore a potential SMB threshold for complete melt of the GrIS; the impact and interplay of surface melt and glacial isostatic adjustment (GIA) in determining this threshold; and whether the GrIS exhibits characteristics commonly associated with tipping points, such as sensitivity to external forcing. To this end, we force the Community Ice Sheet Model v.2 (CISM2) by cycling different SMB climatologies previously calculated at multiple elevation classes with the Community Earth System Model v.2 (CESM2) in a two-way coupled CESM2-CISM2 transient simulation of the global climate and GrIS under high CO₂ forcing. The SMB calculation in CESM2 has been evaluated with contemporary observations and high-resolution modelling and includes an advanced representation of surface melt and snow-firn processes. We find a positive SMB threshold for complete GrIS melt of 230 ± 84 Gtyr^-1, corresponding to a 60 % decrease in SMB and to a global mean warming of +3.4 K compared to pre-industrial CESM2-CISM2 simulated values. In our simulations, a small change in the initial SMB forcing (from 255 to 230 Gtyr^-1) and global mean warming above pre-industrial levels (from +3.2 to +3.4 K) causes an abrupt change in the GrIS final volume (from 50 % mass to nearly complete deglaciation). This nonlinear behaviour is caused by the SMB-elevation feedback, which responds to changes in surface topography due to surface melt and GIA. The GrIS tips from ∼ 50 % mass towards nearly complete melt when the impact of melt-induced surface lowering outweighs that of GIA-induced bedrock uplift and the (initially positive) SMB becomes and remains negative for at least a few thousand years. We also find that the GrIS tips towards nearly complete melt when the ice margin in the central west unpins from a coastal region with high topography and SMB. We show that if we keep the SMB fixed (i.e. no SMB-elevation feedback) in this relatively confined region, the ice sheet retreat is halted and nearly complete GrIS melt is prevented even though the initial SMB forcing is past the threshold. Based on the minimum GrIS configuration in previous paleo-ice-sheet modelling studies, we suggest that the surface topography in the central west might have played a role in preventing larger GrIS loss during the last interglacial period ∼ 130-115 kyrBP.

Iceberg A-68 separated from the Larsen C Ice Shelf in July 2017 and the impact of this event on the local ocean circulation has yet to be assessed. Here, we conduct numerical simulations of ocean dynamics near and below the ice shelf pre- and post-calving. Results agree with in situ and remote observations of the area as they indicate that basal melt is primarily controlled by wintertime sea-ice formation, which in turn produces High Salinity Shelf Water (HSSW). After the calving event, we simulate a 50% increase in HSSW intrusion under the ice shelf, enhancing ocean heat delivery by 30%. This results in doubling of the melt rate under Gipps Ice Rise, suggesting a positive feedback for further retreat that could destabilize the Larsen C Ice Shelf. Assessing the impact of ice-front retreat on the heat delivery under the ice is crucial to better understand ice-shelf dynamics in a warming environment.

The Arctic is the region on Earth that is warming the fastest. At the same time, Arctic sea ice is reducing while the Greenland ice sheet (GrIS) is losing mass at an accelerated pace. Here, we study the seasonal impact of reduced Arctic sea ice on GrIS surface mass balance (SMB), using the Community Earth System Model version 2.1 (CESM2), which features an advanced, interactive calculation of SMB. Addressing the impact of sea-ice reductions on the GrIS SMB from observations is difficult due to the short observational records. Also, signals detected using transient climate simulations may be aliases of other forcings. Here, we analyze dedicated simulations from the Polar Amplification Model Intercomparison Project with reduced Arctic sea ice and compare them with preindustrial sea ice simulations while keeping all other forcings constant. In response to reduced sea ice, the GrIS SMB increases in winter due to increased precipitation, driven by the more humid atmosphere and increasing cyclones. In summer, surface melt increases due to a warmer, more humid atmosphere providing increased energy transfer to the surface through the sensible and latent heat fluxes, which triggers the melt-albedo feedback. Further, warming occurs throughout the entire troposphere over Baffin Bay. This deep warming results in regional enhancement of the 500 hPa geopotential heights over the Baffin Bay and Greenland, increasing blocking and heat advection over the GrIS’ surface. This anomalous circulation pattern has been linked to recent increases in the surface melt of the GrIS.

Earth system/ice-sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea-level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy-balance-based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice-sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO₂ scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.

Hydrological regimes of alpine catchments are expected to be strongly affected by climate change, mostly due to their dependence on snow and ice dynamics. While seasonal changes have been studied extensively, studies on changes in the timing and magnitude of annual extremes remain rare. This study investigates the effects of climate change on runoff patterns in six contrasting Alpine catchments in Austria using a process-based, semi-distributed hydrological model and projections from 14 regional and global climate model combinations for two representative concentration pathways, namely RCP4.5 and RCP8.5. The study catchments represent a spectrum of different hydrological regimes, from pluvial-nival to nivo-glacial, as well as distinct topographies and land forms, characterizing different elevation zones across the eastern Alps to provide a comprehensive picture of future runoff changes. The climate projections are used to model river runoff in 2071-2100, which are then compared to the 1981-2010 reference period for all study catchments. Changes in the timing and magnitude of annual maximum and minimum flows, as well as in monthly runoff and snowmelt, are quantified and analyzed. Our results indicate a substantial shift to earlier occurrences in annual maximum flows by 9 to 31gd and an extension of the potential flood season by 1 to 3 months for high-elevation catchments. For low-elevation catchments, changes in the timing of annual maximum flows are less pronounced. Magnitudes of annual maximum flows are likely to increase by 2g%-18g% under RCP4.5, while no clear changes are projected for four catchments under RCP8.5. The latter is caused by a pronounced increase in evaporation and decrease in snowmelt contributions, which offset increases in precipitation. In the future, minimum annual runoff will occur 13-31gd earlier in the winter months for high-elevation catchments, whereas for low-elevation catchments a shift from winter to autumn by about 15-100gd is projected, with generally larger changes for RCP8.5. While all catchments show an increase in mean magnitude of minimum flows by 7-30% under RCP4.5, this is only the case for four catchments under RCP8.5. Our results suggest a relationship between the elevation of catchments and changes in the timing of annual maximum and minimum flows. For the magnitude of the extreme flows, a relationship is found between catchment elevation and annual minimum flows, whereas this relationship is lacking between elevation and annual maximum flow.

Future Greenland ice sheet (GrIS) melt projections are limited by the lack of explicit melt calculations within most global climate models and the high computational cost of dynamical downscaling with regional climate models (RCMs). Here, we train artificial neural networks (ANNs) to obtain relationships between quantities consistently available from global climate model simulations and annually integrated GrIS surface melt. To this end, we train the ANNs with model output from the Community Earth System Model 2.1 (CESM2), which features an interactive surface melt calculation based on a downscaled surface energy balance. We find that ANNs compare well with an independent CESM2 simulation and RCM simulations forced by a CMIP6 subset. The ANNs estimate a melt increase for 2,081–2,100 ranging from 414 (Formula presented.) 275 Gt (Formula presented.) (SSP1-2.6) to 1,378 (Formula presented.) 555 Gt (Formula presented.) (SSP5-8.5) for the full CMIP6 suite. The primary source of uncertainty throughout the 21st century is the spread of climate model sensitivity.

The Greenland ice sheet (GrIS) is now losing mass at a rate of 0.7 mm of sea level rise (SLR) per year. Here we explore future GrIS evolution and interactions with global and regional climate under high greenhouse gas forcing with the Community Earth System Model version 2.1 (CESM2.1), which includes an interactive ice sheet component (the Community Ice Sheet Model v2.1 [CISM2.1]) and an advanced energy balance-based calculation of surface melt. We run an idealized 350-year scenario in which atmospheric CO₂ concentration increases by 1% annually until reaching four times pre-industrial values at year 140, after which it is held fixed. The global mean temperature increases by 5.2 and 8.5 K by years 131–150 and 331–350, respectively. The projected GrIS contribution to global mean SLR is 107 mm by year 150 and 1,140 mm by year 350. The rate of SLR increases from 2 mm yr⁻¹ at year 150 to almost 7 mm yr⁻¹ by year 350. The accelerated mass loss is caused by rapidly increasing surface melt as the ablation area expands, with associated albedo feedback and increased sensible and latent heat fluxes. This acceleration occurs for a global warming of approximately 4.2 K with respect to pre-industrial and is in part explained by the quasi-parabolic shape of the ice sheet, which favors rapid expansion of the ablation area as it approaches the interior “plateau.”.

Recent research shows increasing decadal ice mass losses from the Greenland and Antarctic Ice Sheets and more generally from glaciers worldwide in the light of continued global warming. Here, in an update of our previous ISMASS paper (Hanna et al., 2013), we review recent observational estimates of ice sheet and glacier mass balance, and their related uncertainties, first briefly considering relevant monitoring methods. Focusing on the response to climate change during 1992–2018, and especially the post-IPCC AR5 period, we discuss recent changes in the relative contributions of ice sheets and glaciers to sea-level change. We assess recent advances in understanding of the relative importance of surface mass balance and ice dynamics in overall ice-sheet mass change. We also consider recent improvements in ice-sheet modelling, highlighting data-model linkages and the use of updated observational datasets in ice-sheet models. Finally, by identifying key deficiencies in the observations and models that hamper current understanding and limit reliability of future ice-sheet projections, we make recommendations to the research community for reducing these knowledge gaps. Our synthesis aims to provide a critical and timely review of the current state of the science in advance of the next Intergovernmental Panel on Climate Change Assessment Report that is due in 2021.

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin-up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet-enabled simulations with CESM2-CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).

The response of the Greenland Ice Sheet (GrIS) to a warmer climate is uncertain on long time scales. Climate models, such as those participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), are used to assess this uncertainty. The Community Earth System Model version 2.1 (CESM2) is a CMIP6 model capable of running climate simulations with either one-way coupling (fixed ice sheet geometry) or two-way coupling (dynamic geometry) to the GrIS. The model features prognostic snow albedo, online downscaling using elevation classes, and a firn pack to refreeze percolating melt water. Here we evaluate the representation of the GrIS surface energy balance and surface mass balance in CESM2 at 1° resolution with fixed GrIS geometry. CESM2 agrees closely with ERA-Interim reanalysis data for key controls on GrIS SMB: surface pressure, sea ice extent, 500 hPa geopotential height, wind speed, and 700 hPa air temperature. Cloudsat-CALIPSO data show that supercooled liquid-containing clouds are adequately represented, whereas comparisons to Moderate Resolution Imaging Spectroradiometer and CM SAF Cloud, Albedo, and Surface Radiation data set from Advanced Very High Resolution Radiometer data second edition data suggest that CESM2 underestimates surface albedo. The seasonal cycle and spatial patterns of surface energy balance and surface mass balance components in CESM2 agree well with regional climate model RACMO2.3p2, with GrIS-integrated melt, refreezing, and runoff bracketed by RACMO2 counterparts at 11 and 1 km. Time series of melt, runoff, and SMB show a break point around 1990, similar to RACMO2. These results suggest that GrIS SMB is realistic in CESM2, which adds confidence to coupled ice sheet-climate experiments that aim to assess the GrIS contribution to future sea level rise.

The Community Earth System Model version 2.1 (CESM2.1) is used to investigate the evolution of the Greenland ice sheet (GrIS) surface mass balance (SMB) under an idealized CO₂ forcing scenario of 1% increase until stabilization at 4× pre-industrial at model year 140. In this simulation, the SMB calculation is coupled with the atmospheric model, using a physically based surface energy balance scheme for melt, explicit calculation of snow albedo, and a realistic treatment of polar snow and firn compaction. By the end of the simulation (years 131–150), the SMB decreases with 994 Gt yr⁻¹ with respect to the pre-industrial SMB, which represents a sea-level rise contribution of 2.8 mm yr⁻¹. For a threshold of 2.7-K global temperature increase with respect to pre-industrial, the rate of expansion of the ablation area increases, the mass loss accelerates due to loss of refreezing capacity and accelerated melt, and the SMB becomes negative 6 years later. Before acceleration, longwave radiation is the most important contributor to increasing energy for melt. After acceleration, the large expansion of the ablation area strongly reduces surface albedo. This and much increased turbulent heat fluxes as the GrIS-integrated summer surface temperature approaches melt point become the major sources of energy for melt.

The Greenland Ice Sheet (GrIS) mass balance is examined with an Earth system/ice sheet model that interactively couples the GrIS to the broader Earth system. The simulation runs from 1850 to 2100, with historical and SSP5-8.5 forcing. By the mid-21st century, the cumulative GrIS contribution to global mean sea level rise (SLR) is 23 mm. During the second half of the 21st century, the surface mass balance becomes negative in all drainage basins, with an additional SLR contribution of 86 mm. The annual mean GrIS mass loss in the last two decades is 2.7-mm sea level equivalent (SLE) year⁻¹. The increased SLR contribution from the surface mass balance (3.1 mm SLE year⁻¹) is partly offset by reduced ice discharge from thinning and retreat of outlet glaciers. The southern GrIS drainage basins contribute 73% of the mass loss in mid-century but 55% by 2100, as surface runoff increases in the northern basins.

A collection of scientific analyses, metrics, and visualizations for robust validation of ice sheet models is presented using the Land Ice Verification and Validation toolkit (LIVVkit), version 2.1, and the LIVVkit Extensions repository (LEX), version 0.1. This software collection targets stand-alone ice sheet or coupled Earth system models, and handles datasets and analyses that require high-performance computing and storage. LIVVkit aims to enable efficient and fully reproducible workflows for postprocessing, analysis, and visualization of observational and model-derived datasets in a shareable format, whereby all data, methodologies, and output are distributed to users for evaluation. Extending from the initial LIVVkit software framework, we demonstrate Greenland ice sheet simulation validation metrics using the coupled Community Earth System Model (CESM) as well as an idealized stand-alone high-resolution Community Ice Sheet Model, version 2 (CISM2), coupled to the Albany/FELIX velocity solver (CISM-Albany or CISM-A). As one example of the capability, LIVVkit analyzes the degree to which models capture the surface mass balance (SMB) and identifies potential sources of bias, using recently available in situ and remotely sensed data as comparison. Related fields within atmosphere and land surface models, e.g., surface temperature, radiation, and cloud cover, are also diagnosed. Applied to the CESM1.0, LIVVkit identifies a positive SMB bias that is focused largely around Greenland's southwest region that is due to insufficient ablation.