Correction to: Nature Reviews Earth & Environment https://doi.org/10.1038/s43017-023-00507-9, published online 23 January 2024.

In the version of the article initially published, in Fig. 5, under “Radar altimeter”, “O(16–160 m)” previously read “O(16–160 km)”. This has now been corrected in the HTML and PDF versions of the article.

Most of the Greenland and Antarctic ice sheets are covered with firn — the transitional material between snow and glacial ice. Firn is vital for understanding ice-sheet mass balance and hydrology, and palaeoclimate. In this Review, we synthesize knowledge of firn, including its formation, observation, modelling and relevance to ice sheets. The refreezing of meltwater in the pore space of firn currently prevents 50% of meltwater in Greenland from running off into the ocean and protects Antarctic ice shelves from catastrophic collapse. Continued atmospheric warming could inhibit future protection against mass loss. For example, warming in Greenland has already contributed to a 5% reduction in firn pore space since 1980. All projections of future firn change suggest that surface meltwater will have an increasing impact on firn, with melt occurring tens to hundreds of kilometres further inland in Greenland, and more extensively on Antarctic ice shelves. Although progress in observation and modelling techniques has led to a well-established understanding of firn, the large uncertainties associated with meltwater percolation processes (refreezing, ice-layer formation and storage) must be reduced further. A tighter integration of modelling components (firn, atmosphere and ice-sheet models) will also be needed to better simulate ice-sheet responses to anthropogenic warming and to quantify future sea-level rise.

Ice losses from the Greenland and Antarctic ice sheets have accelerated since the 1990s, accounting for a significant increase in the global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume, and in Earth's gravity field. Between 1992 and 2020, the ice sheets contributed 21.0±1.9g€¯mm to global mean sea level, with the rate of mass loss rising from 105g€¯Gtg€¯yr-1 between 1992 and 1996 to 372g€¯Gtg€¯yr-1 between 2016 and 2020. In Greenland, the rate of mass loss is 169±9g€¯Gtg€¯yr-1 between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86g€¯Gtg€¯yr-1 in 2017 to 444g€¯Gtg€¯yr-1 in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (82±9g€¯Gtg€¯yr-1) and, to a lesser extent, from the Antarctic Peninsula (13±5g€¯Gtg€¯yr-1). East Antarctica remains close to a state of balance, with a small gain of 3±15g€¯Gtg€¯yr-1, but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at 10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021).

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.

We evaluate modelled Antarctic ice sheet (AIS) near-surface climate, surface mass balance (SMB) and surface energy balance (SEB) from the updated polar version of the regional atmospheric climate model, RACMO2 (1979-2016). The updated model, referred to as RACMO2.3p2, incorporates upper-air relaxation, a revised topography, tuned parameters in the cloud scheme to generate more precipitation towards the AIS interior and modified snow properties reducing drifting snow sublimation and increasing surface snowmelt.

Comparisons of RACMO2 model output with several independent observational data show that the existing biases in AIS temperature, radiative fluxes and SMB components are further reduced with respect to the previous model version. The model-integrated annual average SMB for the ice sheet including ice shelves (minus the Antarctic Peninsula, AP) now amounts to 2229ĝ€Gtĝ€yĝ'1, with an interannual variability of 109ĝ€Gtĝ€yĝ'1. The largest improvement is found in modelled surface snowmelt, which now compares well with satellite and weather station observations. For the high-resolution ( ĝ1/4 ĝ€5.5ĝ€km) AP simulation, results remain comparable to earlier studies.

The updated model provides a new, high-resolution data set of the contemporary near-surface climate and SMB of the AIS; this model version will be used for future climate scenario projections in a forthcoming study.