Observations and models agree that the Greenland Ice Sheet (GrIS) surface mass balance (SMB) has decreased since the end of the 1990s due to an increase in meltwater runoff and that this trend will accelerate in the future. However, large uncertainties remain, partly due to different approaches for modelling the GrIS SMB, which have to weigh physical complexity or low computing time, different spatial and temporal resolutions, different forcing fields, and different ice sheet topographies and extents, which collectively make an inter-comparison difficult. Our GrIS SMB model intercomparison project (GrSMBMIP) aims to refine these uncertainties by intercomparing 13 models of four types which were forced with the same ERA-Interim reanalysis forcing fields, except for two global models. We interpolate all modelled SMB fields onto a common ice sheet mask at 1 km horizontal resolution for the period 1980-2012 and score the outputs against (1) SMB estimates from a combination of gravimetric remote sensing data from GRACE and measured ice discharge; (2) ice cores, snow pits and in situ SMB observations; and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting model deficiencies in an accurate representation of the GrIS ablation zone extent and processes related to surface melt and runoff. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of the same order as RCMs compared with observations and therefore remain useful tools for long-term simulations or coupling with ice sheet models. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present-day SMB relative to observations, suggesting that biases are not systematic among models and that this ensemble estimate can be used as a reference for current climate when carrying out future model developments. However, a higher density of in situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 m w.e. yr-1 due to large discrepancies in modelled snowfall accumulation. @en

Surface meltwater ponding has been implicated as a major driver for recent ice shelf collapse as well as the speedup of tributary glaciers in the northeast Antarctic Peninsula. Surface melt on the NAP is impacted by the strength and frequency of westerly winds, which result in sporadic foehn flow. We estimate changes in the frequency of foehn flow and the associated impact on snow melt, density, and the percolation depth of meltwater over the period 1982–2017 using a regional climate model and passive microwave data. The first of two methods extracts spatial patterns of melt occurrence using empirical orthogonal function analysis. The second method applies the Foehn Index, introduced here to capture foehn occurrence over the full study domain. Both methods show substantial foehn-induced melt late in the melt season since 2015, resulting in compounded densification of the near-surface snow, with potential implications for future ice shelf stability.@en

Ice loss from the Greenland ice sheet is one of the largest sources of contemporary sea-level rise (SLR). While process-based models place timescales on Greenland’s deglaciation, their confidence is obscured by model shortcomings including imprecise atmospheric and oceanic couplings. Here, we present a complementary approach resolving ice sheet disequilibrium with climate constrained by satellite-derived bare-ice extent, tidewater sector ice flow discharge and surface mass balance data. We find that Greenland ice imbalance with the recent (2000–2019) climate commits at least 274 ± 68 mm SLR from 59 ± 15 × 103 km2 ice retreat, equivalent to 3.3 ± 0.9% volume loss, regardless of twenty-first-century climate pathways. This is a result of increasing mass turnover from precipitation, ice flow discharge and meltwater run-off. The high-melt year of 2012 applied in perpetuity yields an ice loss commitment of 782 ± 135 mm SLR, serving as an ominous prognosis for Greenland’s trajectory through a twenty-first century of warming.@en

Between 2003-2016, the Greenland ice sheet (GrIS) was one of the largest contributors to sea level rise, as it lost about 255 Gt of ice per year. This mass loss slowed in 2017 and 2018 to about 100 Gt yr−1. Here we examine further changes in rate of GrIS mass loss, by analyzing data from the GRACE-FO (Gravity Recovery and Climate Experiment – Follow On) satellite mission, launched in May 2018. Using simulations with regional climate models we show that the mass losses observed in 2017 and 2018 by the GRACE and GRACE-FO missions are lower than in any other two year period between 2003 and 2019, the combined period of the two missions. We find that this reduced ice loss results from two anomalous cold summers in western Greenland, compounded by snow-rich autumn and winter conditions in the east. For 2019, GRACE-FO reveals a return to high melt rates leading to a mass loss of 223 ± 12 Gt month−1 during the month of July alone, and a record annual mass loss of 532 ± 58 Gt yr−1.@en

This study presents a data set of daily, 1 km resolution Greenland ice sheet (GrIS) surface mass balance (SMB) covering the period 1958–2015. Applying corrections for elevation, bare ice albedo and accumulation bias, the high-resolution product is statistically downscaled from the native daily output of the polar regional climate model RACMO2.3 at 11 km. The data set includes all individual SMB components projected to a down-sampled version of the Greenland Ice Mapping Project (GIMP) digital elevation model and ice mask. The 1 km mask better resolves narrow ablation zones, valley glaciers, fjords and disconnected ice caps. Relative to the 11 km product, the more detailed representation of isolated glaciated areas leads to increased precipitation over the southeastern GrIS. In addition, the downscaled product shows a significant increase in runoff owing to better resolved low-lying marginal glaciated regions. The combined corrections for elevation and bare ice albedo markedly improve model agreement with a newly compiled data set of ablation measurements.@en

Glaciers in the Arctic respond sensitively to climate change, recording the polar amplification of global warming with increasing mass loss. Here, we use glacier mass balances in Svalbard and northern Arctic Canada to categorize tropospheric variability and the associated summer circulation over the Arctic. We establish a link between annual glacier mass balances and their respective atmospheric forcings since 1950 using GRACE/GRACE-FO satellite data (2002–2021), as well as regional climate models and reanalysis data (1950–2019). We find that asynchronous behaviour of mass balance between the regions has become very likely since the early 2000s, exceeding the range of previous decadal variability. Related tropospheric circulation exhibits more meridional patterns, a greater influence of meridional heat advection and a wavier summer circulation. The traceable impact on glacier mass balances emphasizes the importance of dynamic next to thermodynamic climate changes for the future of glacier mass loss, Arctic ecology and societal impacts.@en

Mass loss from the Greenland Ice Sheet (GrIS) has accelerated over the past two decades, coincident with rapid Arctic warming and increasing moisture transport over Greenland by atmospheric rivers (ARs). Summer ARs affecting western Greenland trigger GrIS melt events, but the physical mechanisms through which ARs induce melt are not well understood. This study elucidates the coupled surface-atmosphere processes by which ARs force GrIS melt through analysis of the surface energy balance (SEB), cloud properties, and local- to synoptic-scale atmospheric conditions during strong summer AR events affecting western Greenland. ARs are identified in MERRA-2 reanalysis (1980-2017) and classified by integrated water vapor transport (IVT) intensity. SEB, cloud, and atmospheric data from regional climate model, observational, reanalysis, and satellite-based datasets are used to analyze melt-inducing physical processes during strong,.90th percentile ''AR901'' events. Near AR ''landfall,'' AR901 days feature increased cloud cover that reduces net shortwave radiation and increases net longwave radiation. As these oppositely signed radiative anomalies partly cancel during AR901 events, increased melt energy in the ablation zone is primarily provided by turbulent heat fluxes, particularly sensible heat flux. These turbulent heat fluxes are driven by enhanced barrier winds generated by a stronger synoptic pressure gradient combined with an enhanced local temperature contrast between cool over-ice air and the anomalously warm surrounding atmosphere. During AR901 events in northwest Greenland, anomalous melt is forced remotely through a clear-sky foehn regime produced by downslope flow in eastern Greenland.@en