With the aim of studying the recent Greenland ice sheet (GrIS) surface mass balance (SMB) decrease relative to the last century, we have forced the regional climate MAR (Modèle Atmosphérique Régional; version 3.5.2) model with the ERA-Interim (ECMWF Interim Re-Analysis; 1979–2015), ERA-40 (1958–2001), NCEP–NCARv1 (National Centers for Environmental Prediction–National Center for Atmospheric Research Reanalysis version 1; 1948–2015), NCEP–NCARv2 (1979–2015), JRA-55 (Japanese 55-year Reanalysis; 1958–2014), 20CRv2(c) (Twentieth Century Reanalysis version 2; 1900–2014) and ERA-20C (1900–2010) reanalyses. While all these forcing products are reanalyses that are assumed to represent the same climate, they produce significant differences in the MAR-simulated SMB over their common period. A temperature adjustment of +1 °C (respectively −1 °C) was, for example, needed at the MAR boundaries with ERA-20C (20CRv2) reanalysis, given that ERA-20C (20CRv2) is ∼ 1 °C colder (warmer) than ERA-Interim over Greenland during the period 1980–2010. Comparisons with daily PROMICE (Programme for Monitoring of the Greenland Ice Sheet) near-surface observations support these adjustments. Comparisons with SMB measurements, ice cores and satellite-derived melt extent reveal the most accurate forcing datasets for the simulation of the GrIS SMB to be ERA-Interim and NCEP–NCARv1. However, some biases remain in MAR, suggesting that some improvements are still needed in its cloudiness and radiative schemes as well as in the representation of the bare ice albedo. Results from all MAR simulations indicate that (i) the period 1961–1990, commonly chosen as a stable reference period for Greenland SMB and ice dynamics, is actually a period of anomalously positive SMB (∼ +40 Gt yr⁻¹) compared to 1900–2010; (ii) SMB has decreased significantly after this reference period due to increasing and unprecedented melt reaching the highest rates in the 120-year common period; (iii) before 1960, both ERA-20C and 20CRv2-forced MAR simulations suggest a significant precipitation increase over 1900–1950, but this increase could be the result of an artefact in the reanalyses that are not well-enough constrained by observations during this period and (iv) since the 1980s, snowfall is quite stable after having reached a maximum in the 1970s. These MAR-based SMB and accumulation reconstructions are, however, quite similar to those from Box (2013) after 1930 and confirm that SMB was quite stable from the 1940s to the 1990s. Finally, only the ERA-20C-forced simulation suggests that SMB during the 1920–1930 warm period over Greenland was comparable to the SMB of the 2000s, due to both higher melt and lower precipitation than normal.

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.9 mm to global mean sea level, with the rate of mass loss rising from 105 Gt yr ⁻¹ between 1992 and 1996 to 372 Gt yr ⁻¹ between 2016 and 2020. In Greenland, the rate of mass loss is 169±9 Gt yr ⁻¹ between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86 Gt yr ⁻¹ in 2017 to 444 Gt yr ⁻¹ 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±9 Gt yr ⁻¹ ) and, to a lesser extent, from the Antarctic Peninsula (13±5 Gt yr ⁻¹ ). East Antarctica remains close to a state of balance, with a small gain of 3±15 Gt yr ⁻¹ , but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at https://doi.org/10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021).

The Antarctic ice sheet mass balance is a major component of the sea level budget and results from the difference of two fluxes of a similar magnitude: ice flow discharging in the ocean and net snow accumulation on the ice sheet surface, i.e. the surface mass balance (SMB). Separately modelling ice dynamics and SMB is the only way to project future trends. In addition, mass balance studies frequently use regional climate models (RCMs) outputs as an alternative to observed fields because SMB observations are particularly scarce on the ice sheet. Here we evaluate new simulations of the polar RCM MAR forced by three reanalyses, ERA-Interim, JRA-55, and MERRA-2, for the period 1979–2015, and we compare MAR results to the last outputs of the RCM RACMO2 forced by ERA-Interim. We show that MAR and RACMO2 perform similarly well in simulating coast-to-plateau SMB gradients, and we find no significant differences in their simulated SMB when integrated over the ice sheet or its major basins. More importantly, we outline and quantify missing or underestimated processes in both RCMs. Along stake transects, we show that both models accumulate too much snow on crests, and not enough snow in valleys, as a result of drifting snow transport fluxes not included in MAR and probably underestimated in RACMO2 by a factor of 3. Our results tend to confirm that drifting snow transport and sublimation fluxes are much larger than previous model-based estimates and need to be better resolved and constrained in climate models. Sublimation of precipitating particles in low-level atmospheric layers is responsible for the significantly lower snowfall rates in MAR than in RACMO2 in katabatic channels at the ice sheet margins. Atmospheric sublimation in MAR represents 363 Gt yr⁻¹ over the grounded ice sheet for the year 2015, which is 16 % of the simulated snowfall loaded at the ground. This estimate is consistent with a recent study based on precipitation radar observations and is more than twice as much as simulated in RACMO2 because of different time residence of precipitating particles in the atmosphere. The remaining spatial differences in snowfall between MAR and RACMO2 are attributed to differences in advection of precipitation with snowfall particles being likely advected too far inland in MAR.

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⁻¹ due to large discrepancies in modelled snowfall accumulation.

We compare the performance of five different regional climate models (RCMs) (COSMO-CLM², HIRHAM5, MAR3.10, MetUM, and RACMO2.3p2), forced by ERA-Interim reanalysis, in simulating the near-surface climate and surface mass balance (SMB) of Antarctica. All models simulate Antarctic climate well when compared with daily observed temperature and pressure, with nudged models matching daily observations slightly better than free-running models. The ensemble mean annual SMB over the Antarctic ice sheet (AIS) including ice shelves is 2329±94 Gt yr⁻¹ over the common 1987–2015 period covered by all models. There is large interannual variability, consistent between models due to variability in the driving ERA-Interim reanalysis. Mean annual SMB is sensitive to the chosen period; over our 30-year climatological mean period (1980 to 2010), the ensemble mean is 2483 Gt yr⁻¹. However, individual model estimates vary from 1961±70 to 2519±118 Gt yr⁻¹. The largest spatial differences between model SMB estimates are in West Antarctica, the Antarctic Peninsula, and around the Transantarctic Mountains. We find no significant trend in Antarctic SMB over either period. Antarctic ice sheet (AIS) mass loss is currently equivalent to around 0.5 mm yr⁻¹ of global mean sea level rise (Shepherd et al., 2020), but our results indicate some uncertainty in the SMB contribution based on RCMs. We compare modelled SMB with a large dataset of observations, which, though biased by undersampling, indicates that many of the biases in SMB are common between models. A drifting-snow scheme improves modelled SMB on ice sheet surface slopes with an elevation between 1000 and 2000 m, where strong katabatic winds form. Different ice masks have a substantial impact on the integrated total SMB and along with model resolution are factored into our analysis. Targeting undersampled regions with high precipitation for observational campaigns will be key to improving future estimates of SMB in Antarctica.

The Antarctic ice sheet (AIS) is sensitive to short‐term extreme meteorological events that can leave long‐term impacts on the continent's surface mass balance (SMB). We investigate the impacts of atmospheric rivers (ARs) on the AIS precipitation budget using an AR detection algorithm and a regional climate model (Modèle Atmosphérique Régional) from 1980 to 2018. While ARs and their associated extreme vapor transport are relatively rare events over Antarctic coastal regions (∼3 days per year), they have a significant impact on the precipitation climatology. ARs are responsible for at least 10% of total accumulated snowfall across East Antarctica (localized areas reaching 20%) and a majority of extreme precipitation events. Trends in AR annual frequency since 1980 are observed across parts of AIS, most notably an increasing trend in Dronning Maud Land; however, interannual variability in AR frequency is much larger. This AR behavior appears to drive a significant portion of annual snowfall trends across East Antarctica, while controlling the interannual variability of precipitation across most of the AIS. AR landfalls are most likely when the circumpolar jet is highly amplified during blocking conditions in the Southern Ocean. There is a fingerprint of the Southern Annular Mode (SAM) on AR variability in West Antarctica with SAM+ (SAM−) favoring increased AR frequency in the Antarctic Peninsula (Amundsen‐Ross Sea coastline). Given the relatively large influence ARs have on precipitation across the continent, it is advantageous for future studies of moisture transport to Antarctica to consider an AR framework especially when considering future SMB changes.