Search

Near-surface winds play a crucial role in the climate of Antarctica, but accurately quantifying and understanding their drivers is complex. They result from the contribution of two distinct families of drivers: the large-scale pressure gradient and surface-induced pressure gradients known as katabatic and thermal wind. The extrapolation of vertical potential temperature above the boundary layer down to the surface enables us to separate and quantify the contribution of these different pressure gradients in the momentum budget equations. Using this method applied to outputs of the regional atmospheric model MAR at a 3-hourly resolution, we find that the seasonal and spatial variability in near-surface winds in Adélie Land is dominated by surface processes. On the other hand, high-frequency temporal variability (3-hourly) is mainly controlled by large-scale variability everywhere in Antarctica, except on the coast. In coastal regions, although the katabatic acceleration surpasses all other accelerations in magnitude, none of the katabatic or large-scale accelerations can be identified as the single primary driver of near-surface wind variability. The angle between the large-scale acceleration and the surface slope is a key factor in explaining strong wind speed events: the highest-wind-speed events happen when the katabatic and large-scale forcing are aligned, although each acceleration, when acting alone, can also cause strong wind speed. This study underlines the complexity of the drivers of Antarctic surface winds and the value of the momentum budget decomposition to identify drivers at different spatial and temporal scales.

The Greenland Ice Sheet is a key contributor to sea level rise. By melting, the ice sheet thins, inducing higher surface melt due to lower surface elevations, accelerating the melt coming from global warming. This process is called the melt–elevation feedback and can be considered by using two types of models: either (1) atmospheric models, which can represent the surface mass balance (SMB), or SMB estimates resulting from simpler models such as positive degree day models or (2) ice sheet models representing the surface elevation evolution. The latter ones do not represent the surface mass balance explicitly as well as polar-oriented climate models. A new coupling between the MAR (Modèle Atmosphérique Régional) regional climate model and the PISM (Parallel Ice Sheet Model) ice sheet model is presented here following the CESM2 (Community Earth System Model; SSP5-8.5, Shared Socioeconomic Pathway) scenario until 2100 at the MAR lateral boundaries. The coupling is extended to 2200 with a stabilised climate (+7 ∘C compared to 1961–1990) by randomly sampling the last 10 years of CESM2 to force MAR and reaches a sea level rise contribution of 64 cm. The fully coupled simulation is compared to a one-way experiment where surface topography remains fixed in MAR. However, the surface mass balance is corrected for the melt–elevation feedback when interpolated on the PISM grid by using surface mass balance vertical gradients as a function of local elevation variations (offline correction). This method is often used to represent the melt–elevation feedback and prevents a coupling which is too expensive in computation time. In the fully coupled MAR simulation, the ice sheet morphology evolution (changing slope and reducing the orographic barrier) induces changes in local atmospheric patterns. More specifically, wind regimes are modified, as well as temperature lapse rates, influencing the melt rate through modification of sensible heat fluxes at the ice sheet margins. We highlight mitigation of the melt lapse rate on the margins by modifying the surface morphology. The lapse rates considered by the offline correction are no longer valid at the ice sheet margins. If used (one-way simulation), this correction implies an overestimation of the sea level rise contribution of 2.5 %. The mitigation of the melt lapse rate on the margins can only be corrected by using a full coupling between an ice sheet model and an atmospheric model.

We use an observationally calibrated ice-sheet model to investigate the future trajectory of the Antarctic ice sheet related to uncertainties in the future balance between sub-shelf melting and ice discharge, on the one hand, and the surface mass balance, on the other. Our ensemble of simulations, forced by a panel of climate models from the sixth phase of the Coupled Model Intercomparison Project (CMIP6), suggests that the ocean will be the primary driver of short-term Antarctic mass loss, initiating ice loss in West Antarctica already during this century. The atmosphere initially plays a mitigating role through increased snowfall, leading to an Antarctic contribution to global mean sea-level rise by 2100 of 6 (−8 to 15) cm under a low-emission scenario and 5.5 (−10 to 16) cm under a very high-emission scenario. However, under the very high-emission pathway, the influence of the atmosphere shifts beyond the end of the century, becoming an amplifying driver of mass loss as the ice sheet's surface mass balance decreases. We show that this transition occurs when Antarctic near-surface warming exceeds a critical threshold of +7.5 ∘C, at which the increase in surface runoff outweighs the increase in snow accumulation, a signal that is amplified by the melt–elevation feedback. Therefore, under the very high-emission scenario, oceanic and atmospheric drivers are projected to result in a complete collapse of the West Antarctic ice sheet along with significant grounding-line retreat in the marine basins of the East Antarctic ice sheet, leading to a median global mean sea-level rise of 2.75 (6.95) m by 2300 (3000). Under a more sustainable socio-economic pathway, we find that the Antarctic ice sheet may still contribute to a median global mean sea-level rise of 0.62 (1.85) m by 2300 (3000). However, the rate of sea-level rise is significantly reduced as mass loss is likely to remain confined to the Amundsen Sea Embayment, where present-day climate conditions seem sufficient to commit to a continuous retreat of Thwaites Glacier.

Antarctica, the coldest and driest continent, is home to the largest ice sheet, whose mass is predominantly recharged by snowfall. A common feature of polar regions is the warming associated with snowfall, as moist oceanic air and cloud cover increase the surface temperature. Consequently, snow that accumulates on the ice sheet is deposited under unusually warm conditions. Here we use a polar-oriented regional atmospheric model to study the statistical difference between average and snowfall-weighted temperatures. During snowfall, the warm anomaly scales with snowfall amount, with the strongest sensitivity occurring at low-accumulation sites. Heavier snowfall in winter helps to decrease the annual snowfall-weighted temperature, but this effect is overwritten by the event-scale warming associated with precipitating atmospheric systems, which particularly contrast with the extremely cold conditions that occur in winter. Consequently, the seasonal range of snowfall-weighted temperature is reduced by 20 %. On the other hand, the annual snowfall-weighted temperature shows 80 % more interannual variability than the annual temperature due to the irregularity of snowfall occurrence and its associated temperature anomaly. Disturbances of the apparent annual temperature cycle and interannual variability have important consequences for the interpretation of water isotopes in precipitation, which are deposited with snowfall and commonly used for paleotemperature reconstructions from ice cores.

Global warming has already substantially altered the Arctic cryosphere. Due to the Arctic warming amplification, the temperature is increasing more strongly, leading to pervasive changes in this area. Recent years were notably marked by melt records over the Greenland Ice Sheet, while other regions such as Svalbard seem to remain less influenced. This raises the question of the current state of the Greenland Ice Sheet and the various ice caps in the Arctic for which few studies are available. Here, we run the regional climate model (RCM) Modèle Atmosphérique Régional (MAR) at a resolution of 6 km over four different domains covering all Arctic land ice to produce a unified surface mass balance product from 1950 to the present day. We also compare our results to large-scale indices to better understand the heterogeneity of the evolutions across the Arctic and their links to recent climate change. We find a sharp decrease of surface mass balance (SMB) over the western Arctic (Canada and Greenland) in relationship with the atmospheric blocking situations that have become more frequent in summer, resulting in a 41 % increase of the melt rate since 1950. This increase is not seen over the Russian Arctic permanent ice areas, where melt rates have increased by only 3 % on average, illustrating a heterogeneity in the Arctic SMB response to global warming.

Both regional climate models (RCMs) and remote sensing (RS) data are essential tools in understanding the response of polar regions to climate change. RCMs can simulate how certain climate variables, such as surface melt, runoff and snowfall, are likely to change in response to different climate scenarios but are subject to biases and errors. RS data can assist in reducing and quantifying model uncertainties by providing indirect observations of the modeled variables on the present climate. In this work, we improve on an existing scheme to assimilate RS wet snow occurrence data with the “Modèle Atmosphérique Régional” (MAR) RCM and investigate the sensitivity of the RCM to the parameters of the scheme. The assimilation is performed by nudging the MAR snowpack temperature to match the presence of liquid water observed by satellites. The sensitivity of the assimilation method is tested by modifying parameters such as the depth to which the MAR snowpack is warmed or cooled, the quantity of water required to qualify a MAR pixel as “wet” (0.1 % or 0.2 % of the snowpack mass being water), and assimilating different RS datasets. Data assimilation is carried out on the Antarctic Peninsula for the 2019–2021 period. The results show an increase in meltwater production (+66.7 % on average, or +95 Gt), along with a small decrease in surface mass balance (SMB) (−4.5 % on average, or −20 Gt) for the 2019–2020 melt season after assimilation. The model is sensitive to the tested parameters, albeit with varying orders of magnitude. The prescribed warming depth has a larger impact on the resulting surface melt production than the liquid water content (LWC) threshold due to strong refreezing occurring within the top layers of the snowpack. The values tested for the LWC threshold are lower than the LWC for typical melt days (approximately 1.2 %) and impact results mainly at the beginning and end of the melting period. The assimilation method will allow for the estimation of uncertainty in MAR meltwater production and will enable the identification of potential issues in modeling near-surface snowpack processes, paving the way for more accurate simulations of snow processes in model projections.

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−1 between 1992 and 1996 to 372 Gt yr−1 between 2016 and 2020. In Greenland, the rate of mass loss is 169±9 Gt yr−1 between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86 Gt yr−1 in 2017 to 444 Gt 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±9 Gt yr−1) and, to a lesser extent, from the Antarctic Peninsula (13±5 Gt yr−1). East Antarctica remains close to a state of balance, with a small gain of 3±15 Gt yr−1, 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).

Regional climate models (RCMs) and reanalysis datasets provide valuable information for assessing the vulnerability of ice shelves to collapse over Antarctica, which is important for future global sea level rise estimates. Within this context, this paper examines variability in snowfall, near-surface air temperature and melt across products from the Met Office Unified Model (MetUM), Regional Atmospheric Climate Model (RACMO) and Modèle Atmosphérique Régional (MAR) RCMs, as well as the ERA-Interim and ERA5 reanalysis datasets. Seasonal and trend decomposition using LOESS (STL) is applied to split the monthly time series at each model grid cell into trend, seasonal and residual components. Significant systematic differences between outputs are shown for all variables in the mean and in the seasonal and residual standard deviations, occurring at both large and fine spatial scales across Antarctica. Results imply that differences in the atmospheric dynamics, parametrisation, tuning and surface schemes between models together contribute more significantly to large-scale variability than differences in the driving data, resolution, domain specification, ice sheet mask, digital elevation model and boundary conditions. Despite significant systematic differences, high temporal correlations are found for snowfall and near-surface air temperature across all products at fine spatial scales. For melt, only moderate correlation exists at fine spatial scales between different RCMs and low correlation between RCM and reanalysis outputs. Root mean square deviations (RMSDs) between all outputs in the monthly time series for each variable are shown to be significant at fine spatial scales relative to the magnitude of annual deviations. Correcting for systematic differences results in significant reductions in RMSDs, suggesting the importance of observations and further development of bias-correction techniques.

Recent warm atmospheric conditions have damaged the ice shelves of the Antarctic Peninsula through surface melt and hydrofracturing and could potentially initiate future collapse of other Antarctic ice shelves. However, model projections with similar greenhouse gas scenarios suggest large differences in cumulative 21st-century surface melting. So far it remains unclear whether these differences are due to variations in warming rates in individual models or whether local feedback mechanisms of the surface energy budget could also play a notable role. Here we use the polar-oriented regional climate model MAR (Modèle Atmosphérique Régional) to study the physical mechanisms that would control future surface melt over the Antarctic ice shelves in high-emission scenarios RCP8.5 and SSP5-8.5. We show that clouds enhance future surface melt by increasing the atmospheric emissivity and longwave radiation towards the surface. Furthermore, we highlight that differences in meltwater production for the same climate warming rate depend on cloud properties and particularly cloud phase. Clouds containing a larger amount of supercooled liquid water lead to stronger melt, subsequently favouring the absorption of solar radiation due to the snowmelt–albedo feedback. As liquid-containing clouds are projected to increase the melt spread associated with a given warming rate, they could be a major source of uncertainties in projections of the future Antarctic contribution to sea level rise.

Regional climate models compute ice sheet surface mass balance (SMB) over a mask that defines the area covered by glacier ice, but ice masks have not been harmonised between models. Intercomparison studies of modelled SMB therefore use a common ice mask. The SMB in areas outside the common ice mask, which are typically coastal and high-precipitation regions, is discarded. Ice mask differences change integrated SMB by between 40.5 and 140.6 Gt yr−1 (1.8 % to 6.0 % of ensemble mean SMB), equivalent to the entire Antarctic mass imbalance. We conclude there is a pressing need for a common ice mask protocol.

We introduce PARASO, a novel five-component fully coupled regional climate model over an Antarctic circumpolar domain covering the full Southern Ocean. The state-of-the-art models used are the fast Elementary Thermomechanical Ice Sheet model (f.ETISh) v1.7 (ice sheet), the Nucleus for European Modelling of the Ocean (NEMO) v3.6 (ocean), the Louvain-la-Neuve sea-ice model (LIM) v3.6 (sea ice), the COnsortium for Small-scale MOdeling (COSMO) model v5.0 (atmosphere) and its CLimate Mode (CLM) v4.5 (land), which are here run at a horizontal resolution close to 1/4∘. One key feature of this tool resides in a novel two-way coupling interface for representing ocean–ice-sheet interactions, through explicitly resolved ice-shelf cavities. The impact of atmospheric processes on the Antarctic ice sheet is also conveyed through computed COSMO-CLM–f.ETISh surface mass exchange. In this technical paper, we briefly introduce each model's configuration and document the developments that were carried out in order to establish PARASO. The new offline-based NEMO–f.ETISh coupling interface is thoroughly described. Our developments also include a new surface tiling approach to combine open-ocean and sea-ice-covered cells within COSMO, which was required to make this model relevant in the context of coupled simulations in polar regions. We present results from a 2000–2001 coupled 2-year experiment. PARASO is numerically stable and fully operational. The 2-year simulation conducted without fine tuning of the model reproduced the main expected features, although remaining systematic biases provide perspectives for further adjustment and development.

In order to understand the evolution of the climate of Antarctica, dominant processes that control surface and low-atmosphere meteorology need to be accurately captured in climate models. We used the regional climate model MAR (v3.11) at 10 km horizontal resolution, forced by ERA5 reanalysis over a 9-year period (2010–2018) to study the impact of drifting snow (designating here the wind-driven transport of snow particles below and above 2 m) on the near-surface atmosphere and surface in Adelie Land, East Antarctica. Two model runs were performed, one with and one without drifting snow, and compared to half-hourly in situ observations at D17, a coastal and windy location of Adelie Land. We show that sublimation of drifting-snow particles in the atmosphere drives the difference between model runs and is responsible for significant impacts on the near-surface atmosphere. By cooling the low atmosphere and increasing its relative humidity, drifting snow also reduces sensible and latent heat exchanges at the surface (−5.7 W m−2 on average). Moreover, large and dense drifting-snow layers act as near-surface cloud by interacting with incoming radiative fluxes, enhancing incoming longwave radiation and reducing incoming shortwave radiation in summer (net radiative forcing: 5.7 W m−2). Even if drifting snow modifies these processes involved in surface–atmosphere interactions, the total surface energy budget is only slightly modified by introducing drifting snow because of compensating effects in surface energy fluxes. The drifting-snow driven effects are not prominent near the surface but peak higher in the boundary layer (fourth vertical level, 12 m) where drifting-snow sublimation is the most pronounced. Accounting for drifting snow in MAR generally improves the comparison at D17, especially for the representation of relative humidity (mean bias reduced from −14.0 % to −0.7 %) and incoming longwave radiation (mean bias reduced from −20.4 W m−2 to −14.9 W m−2). Consequently, our results suggest that a detailed representation of drifting-snow processes is required in climate models to better capture the near-surface meteorology and surface–atmosphere interactions in coastal Adelie Land.

We compare the performance of five different regional climate models (RCMs) (COSMO-CLM2, 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−1 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−1. However, individual model estimates vary from 1961±70 to 2519±118 Gt yr−1. 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−1 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.

Drifting snow, or the wind-driven transport of snow particles originating from clouds and the surface below and above 2 m above ground and their concurrent sublimation, is a poorly documented process on the Antarctic ice sheet, which is inherently lacking in most climate models. Since drifting snow mostly results from erosion of surface particles, a comprehensive evaluation of this process in climate models requires a concurrent assessment of simulated drifting-snow transport and the surface mass balance (SMB). In this paper a new version of the drifting-snow scheme currently embedded in the regional climate model MAR (v3.11) is extensively described. Several important modifications relative to previous version have been implemented and include notably a parameterization for drifting-snow compaction of the uppermost snowpack layer, differentiated snow density at deposition between precipitation and drifting snow, and a rewrite of the threshold friction velocity above which snow erosion initiates. Model results at high resolution (10 km) over Adélie Land, East Antarctica, for the period 2004–2018 are presented and evaluated against available near-surface meteorological observations at half-hourly resolution and annual SMB estimates. The evaluation demonstrates that MAR resolves the local drifting-snow frequency and transport up to the scale of the drifting-snow event and captures the resulting observed climate and SMB variability, suggesting that this model version can be used for continent-wide applications.

The Greenland Ice Sheet (GrIS) will be losing mass at an accelerating pace throughout the 21st century, with a direct link between anthropogenic greenhouse gas emissions and the magnitude of Greenland mass loss. Currently, approximately 60 % of the mass loss contribution comes from surface melt and subsequent meltwater runoff, while 40 % are due to ice calving. In the ablation zone covered by bare ice in summer, most of the surface melt energy is provided by absorbed shortwave fluxes, which could be reduced by solar geoengineering measures. However, so far very little is known about the potential impacts of an artificial reduction in the incoming solar radiation on the GrIS surface energy budget and the subsequent change in meltwater production. By forcing the regional climate model MAR with the latest CMIP6 shared socioeconomic pathways (SSP) future emission scenarios (SSP245, SSP585) and associated G6solar experiment from the CNRM-ESM2-1 Earth system model, we estimate the local impact of a reduced solar constant on the projected GrIS surface mass balance (SMB) decrease. Overall, our results show that even in the case of a low-mitigation greenhouse gas emissions scenario (SSP585), the Greenland surface mass loss can be brought in line with the medium-mitigation emissions scenario (SSP245) by reducing the solar downward flux at the top of the atmosphere by ∼ 40 W/m2 or ∼ 1.5 % (using the G6solar experiment). In addition to reducing global warming in line with SSP245, G6solar also decreases the efficiency of surface meltwater production over the Greenland ice sheet by damping the well-known positive melt–albedo feedback. With respect to a MAR simulation where the solar constant remains unchanged, decreasing the solar constant according to G6solar in the MAR radiative scheme mitigates the projected Greenland ice sheet surface melt increase by 6 %. However, only more constraining geoengineering experiments than G6solar would allow us to maintain a positive SMB until the end of this century without any reduction in our greenhouse gas emissions.

The future surface mass balance (SMB) will influence the ice dynamics and the contribution of the Antarctic ice sheet (AIS) to the sea level rise. Most of recent Antarctic SMB projections were based on the fifth phase of the Coupled Model Intercomparison Project (CMIP5). However, new CMIP6 results have revealed a +1.3 ∘C higher mean Antarctic near-surface temperature than in CMIP5 at the end of the 21st century, enabling estimations of future SMB in warmer climates. Here, we investigate the AIS sensitivity to different warmings with an ensemble of four simulations performed with the polar regional climate model Modèle Atmosphérique Régional (MAR) forced by two CMIP5 and two CMIP6 models over 1981–2100. Statistical extrapolation enables us to expand our results to the whole CMIP5 and CMIP6 ensembles. Our results highlight a contrasting effect on the future grounded ice sheet and the ice shelves. The SMB over grounded ice is projected to increase as a response to stronger snowfall, only partly offset by enhanced meltwater run-off. This leads to a cumulated sea-level-rise mitigation (i.e. an increase in surface mass) of the grounded Antarctic surface by 5.1 ± 1.9 cm sea level equivalent (SLE) in CMIP5-RCP8.5 (Relative Concentration Pathway 8.5) and 6.3 ± 2.0 cm SLE in CMIP6-ssp585 (Shared Socioeconomic Pathways 585). Additionally, the CMIP6 low-emission ssp126 and intermediate-emission ssp245 scenarios project a stabilized surface mass gain, resulting in a lower mitigation to sea level rise than in ssp585. Over the ice shelves, the strong run-off increase associated with higher temperature is projected to decrease the SMB (more strongly in CMIP6-ssp585 compared to CMIP5-RCP8.5). Ice shelves are however predicted to have a close-to-present-equilibrium stable SMB under CMIP6 ssp126 and ssp245 scenarios. Future uncertainties are mainly due to the sensitivity to anthropogenic forcing and the timing of the projected warming. While ice shelves should remain at a close-to-equilibrium stable SMB under the Paris Agreement, MAR projects strong SMB decrease for an Antarctic near-surface warming above +2.5 ∘C compared to 1981–2010 mean temperature, limiting the warming range before potential irreversible damages on the ice shelves. Finally, our results reveal the existence of a potential threshold (+7.5 ∘C) that leads to a lower grounded-SMB increase. This however has to be confirmed in following studies using more extreme or longer future scenarios.

We present projections of West Antarctic surface mass balance (SMB) and surface melt to 2080–2100 under the RCP8.5 scenario and based on a regional model at 10 km resolution. Our projections are built by adding a CMIP5 (Coupled Model Intercomparison Project Phase 5) multi-model-mean seasonal climate-change anomaly to the present-day model boundary conditions. Using an anomaly has the advantage to reduce CMIP5 model biases, and a perfect-model test reveals that our approach captures most characteristics of future changes despite a 16 %–17 % underestimation of projected SMB and melt rates. SMB over the grounded ice sheet in the sector between Getz and Abbot increases from 336 Gt yr−1 in 1989–2009 to 455 Gt yr−1 in 2080–2100, which would reduce the global sea level changing rate by 0.33 mm yr−1. Snowfall indeed increases by 7.4 % ∘C−1 to 8.9 % ∘C−1 of near-surface warming due to increasing saturation water vapour pressure in warmer conditions, reduced sea-ice concentrations, and more marine air intrusion. Ice-shelf surface melt rates increase by an order of magnitude in the 21st century mostly due to higher downward radiation from increased humidity and to reduced albedo in the presence of melting. There is a net production of surface liquid water over eastern ice shelves (Abbot, Cosgrove, and Pine Island) but not over western ice shelves (Thwaites, Crosson, Dotson, and Getz). This is explained by the evolution of the melt-to-snowfall ratio: below a threshold of 0.60 to 0.85 in our simulations, firn air is not entirely depleted by melt water, while entire depletion and net production of surface liquid water occur for higher ratios. This suggests that western ice shelves might remain unaffected by hydrofracturing for more than a century under RCP8.5, while eastern ice shelves have a high potential for hydrofracturing before the end of this century.

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.

The ERA5 reanalysis, recently made available by the European Centre for Medium-Range Weather Forecasts (ECMWF), is a new reanalysis product at a high resolution replacing ERA-Interim and is considered to provide the best climate reanalysis over Greenland to date. However, so far little is known about the performance of ERA5 over the Greenland Ice Sheet (GrIS). In this study, we compare the near-surface climate from the new ERA5 reanalysis to ERA-Interim, the Arctic System Reanalysis (ASR) as well as to a state-of-the-art polar regional climate model (MAR). The results show (1) that ERA5 does not outperform ERA-Interim significantly when compared with near-surface climate observations over GrIS, but ASR better models the near-surface temperature than both ERA reanalyses. (2) Polar regional climate models (e.g., MAR) are still a useful tool to downscale the GrIS climate compared to ERA5, as in particular the near-surface temperature in summer has a key role for representing snow and ice processes such as the surface melt. However, assimilating satellite data and using a more recent radiative scheme enable both ERA and ASR reanalyses to represent more satisfactorily than MAR the downward solar and infrared fluxes. (3) MAR near-surface climate is not affected when forced at its lateral boundaries by either ERA5 or ERA-Interim. Therefore, forcing polar regional climate models with ERA5 starting from 1950 will enable long and homogeneous surface mass balance reconstructions.

Understanding the interannual variability of surface mass balance (SMB) and surface melting in Antarctica is key to quantify the signal-to-noise ratio in climate trends, identify opportunities for multi-year climate predictions and assess the ability of climate models to respond to climate variability. Here we simulate summer SMB and surface melting from 1979 to 2017 using the Regional Atmosphere Model (MAR) at 10 km resolution over the drainage basins of the Amundsen Sea glaciers in West Antarctica. Our simulations reproduce the mean present-day climate in terms of near-surface temperature (mean overestimation of 0.10 ∘C), near-surface wind speed (mean underestimation of 0.42 m s−1), and SMB (relative bias  % over Thwaites glacier). The simulated interannual variability of SMB and melting is also close to observation-based estimates. For all the Amundsen glacial drainage basins, the interannual variability of summer SMB and surface melting is driven by two distinct mechanisms: high summer SMB tends to occur when the Amundsen Sea Low (ASL) is shifted southward and westward, while high summer melt rates tend to occur when ASL is shallower (i.e. anticyclonic anomaly). Both mechanisms create a northerly flow anomaly that increases moisture convergence and cloud cover over the Amundsen Sea and therefore favors snowfall and downward longwave radiation over the ice sheet. The part of interannual summer SMB variance explained by the ASL longitudinal migrations increases westward and reaches 40 % for Getz. Interannual variation in the ASL relative central pressure is the largest driver of melt rate variability, with 11 % to 21 % of explained variance (increasing westward). While high summer SMB and melt rates are both favored by positive phases of El Niño–Southern Oscillation (ENSO), the Southern Oscillation Index (SOI) only explains 5 % to 16 % of SMB or melt rate interannual variance in our simulations, with moderate statistical significance. However, the part explained by SOI in the previous austral winter is greater, suggesting that at least a part of the ENSO–SMB and ENSO–melt relationships in summer is inherited from the previous austral winter. Possible mechanisms involve sea ice advection from the Ross Sea and intrusions of circumpolar deep water combined with melt-induced ocean overturning circulation in ice shelf cavities. Finally, we do not find any correlation with the Southern Annular Mode (SAM) in summer.

Sublimation of snow particles during transport has been recognized as an important ablation process on the Antarctic ice sheet. The resulting increase in moisture content and cooling of the ambient air are thermodynamic negative feedbacks that both contribute to increase the relative humidity of the air, inhibiting further sublimation when saturation is reached. This self-limiting effect and the associated development of saturated near-surface air layers in drifting snow conditions have mainly been described through modelling studies and a few field observations. A set of meteorological data, including drifting snow mass fluxes and vertical profiles of relative humidity, collected at site D17 in coastal Adélie Land (East Antarctica) during 2013 is used to study the relationship between saturation of the near-surface atmosphere and the occurrence of drifting snow in a katabatic wind region that is among the most prone to snow transport by wind. Atmospheric moistening by the sublimation of the windborne snow particles generally results in a strong increase in relative humidity with the magnitude of drifting snow and a decrease in its vertical gradient, suggesting that windborne-snow sublimation can be an important contributor to the local near-surface moisture budget. Despite a high incidence of drifting snow at the measurement location (60.1 % of the time), saturation, when attained, is however most often limited to a thin air layer below 1 m above ground. The development of a near-surface saturated air layer up to the highest measurement level of 5.5 m is observed in only 8.2 % of the drifting snow occurrences or 6.3 % of the time and mainly occurs in strong wind speed and drift conditions. This relatively rare occurrence of ambient saturation is explained by the likely existence of moisture-removal mechanisms inherent to the katabatic and turbulent nature of the boundary-layer flow that weaken the negative feedback of windborne-snow sublimation. Such mechanisms, potentially quite active in katabatic-generated windborne-snow layers all over Antarctica, may be very important in understanding the surface mass and atmospheric moisture budgets of the ice sheet by enhancing windborne-snow sublimation.

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−1 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.

Estimates for the recent period and projections of the Antarctic surface mass balance (SMB) often rely on high-resolution polar-oriented regional climate models (RCMs). However, RCMs require large-scale boundary forcing fields prescribed by reanalyses or general circulation models (GCMs). Since the recent variability of sea surface conditions (SSCs, namely sea ice concentration, SIC, and sea surface temperature, SST) over the Southern Ocean is not reproduced by most GCMs from the 5th phase of the Coupled Model Intercomparison Project (CMIP5), RCMs are then subject to potential biases. We investigate here the direct sensitivity of the Antarctic SMB to SSC perturbations around the Antarctic. With the RCM “Modèle Atmosphérique Régional” (MAR), different sensitivity experiments are performed over 1979–2015 by modifying the ERA-Interim SSCs with (i) homogeneous perturbations and (ii) mean anomalies estimated from all CMIP5 models and two extreme ones, while atmospheric lateral boundary conditions remained unchanged. Results show increased (decreased) precipitation due to perturbations inducing warmer, i.e. higher SST and lower SIC (colder, i.e. lower SST and higher SIC), SSCs than ERA-Interim, significantly affecting the SMB of coastal areas, as precipitation is mainly related to cyclones that do not penetrate far into the continent. At the continental scale, significant SMB anomalies (i.e greater than the interannual variability) are found for the largest combined SST/SIC perturbations. This is notably due to moisture anomalies above the ocean, reaching sufficiently high atmospheric levels to influence accumulation rates further inland. Sensitivity experiments with warmer SSCs based on the CMIP5 biases reveal integrated SMB anomalies (+5 % to +13 %) over the present climate (1979–2015) in the lower range of the SMB increase projected for the end of the 21st century.

Since the 2000s, a change in the atmospheric circulation over the North Atlantic resulting in more frequent blocking events has favoured warmer and sunnier weather conditions over the Greenland Ice Sheet (GrIS) in summer, enhancing the melt increase. This circulation change is not represented by general circulation models (GCMs) of the Coupled Model Intercomparison Project Phase 5 (CMIP5), which do not predict any circulation change for the next century over the North Atlantic. The goal of this study is to evaluate the impact of an atmospheric circulation change (as currently observed) on projections of the future GrIS surface mass balance (SMB). We compare GrIS SMB estimates simulated by the regional climate model MAR forced by perturbed reanalysis (ERA-Interim with a temperature correction of +1, +1.5, and +2 °C at the MAR lateral boundaries) over 1980–2016 to projections of the future GrIS SMB from MAR simulations forced by three GCMs over selected periods for which a similar temperature increase of +1, +1.5, and +2 °C is projected by the GCMs in comparison to 1980–1999. Mean SMB anomalies produced with perturbed reanalysis over the climatologically stable period 1980–1999 are similar to those produced with MAR forced by GCMs over future periods characterised by a similar warming over Greenland. However, over the 2 last decades (2000–2016) when an increase in the frequency of blocking events has been observed in summer, MAR forced by perturbed reanalysis suggests that the SMB decrease could be amplified by a factor of 2 if such atmospheric conditions persist compared to projections forced by GCMs for the same temperature increase but without any circulation change.

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−1) 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.