Abstract. Blowing snow impacts Antarctic ice sheet surface mass balance by snow redistribution and sublimation. However, numerical models poorly represent blowing snow processes, while direct observations are limited in space and time. Satellite retrieval of blowing snow is hindered by clouds and only the strongest events are considered. Here, we develop a blowing snow detection (BSD) algorithm for ground-based remote-sensing ceilometers in polar regions and apply it to ceilometers at Neumayer III and Princess Elisabeth (PE) stations, East Antarctica. The algorithm is able to detect (heavy) blowing snow layers reaching 30 m height. Results show that 78 % of the detected events are in agreement with visual observations at Neumayer III station. The BSD algorithm detects heavy blowing snow 36 % of the time at Neumayer (2011Neumayer ( -2015 and 13 % at PE station (2010)(2011)(2012)(2013)(2014)(2015)(2016). Blowing snow occurrence peaks during the austral winter and shows around 5 % interannual variability. The BSD algorithm is capable of detecting blowing snow both lifted from the ground and occurring during precipitation, which is an added value since results indicate that 92 % of the blowing snow is during synoptic events, often combined with precipitation. Analysis of atmospheric meteorological variables shows that blowing snow occurrence strongly depends on fresh snow availability in addition to wind speed. This finding challenges the commonly used parametrizations, where the threshold for snow particles to be lifted is a function of wind speed only. Blowing snow occurs predominantly during storms and overcast conditions, shortly after precipitation events, and can reach up to 1300 m a.g.l. in the case of heavy mixed events (precipitation and blowing snow together). These results suggest that synoptic conditions play an important role in generating blowing snow events and that fresh snow availability should be considered in determining the blowing snow onset.

In this study, we evaluate output of near-surface atmospheric variables over the Antarctic Ice Sheet from four reanalyses: the new European Centre for Medium-Range Weather Forecasts ERA-5 and its predecessor ERA-Interim, the Climate Forecast System Reanalysis (CFSR), and the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). The near-surface temperature, wind speed, and relative humidity are compared with datasets of in situ observations, together with an assessment of the simulated surface mass balance (approximated by precipitation minus evaporation). No reanalysis clearly stands out as the best performing for all areas, seasons, and variables, and each of the reanalyses displays different biases. CFSR strongly overestimates the relative humidity during all seasons whereas ERA-5 and MERRA-2 (and, to a lesser extent, ERA-Interim) strongly underestimate relative humidity during winter. ERA-5 captures the seasonal cycle of near-surface temperature best and shows the smallest bias relative to the observations. The other reanalyses show a general temperature underestimation during the winter months in the Antarctic interior and overestimation in the coastal areas. All reanalyses underestimate the mean near-surface winds in the interior (except MERRA-2) and along the coast during the entire year. The winds at the Antarctic Peninsula are overestimated by all reanalyses except MERRA-2. All models are able to capture snowfall patterns related to atmospheric rivers, with varying accuracy. Accumulation is best represented by ERA-5, although it underestimates observed surface mass balance and there is some variability in the accumulation over the different elevation classes, for all reanalyses.

Abstract. 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.

Abstract. Antarctic ice sheet mass loss is currently equivalent to around 1 mm year−1 of global mean sea level rise. Most mass is lost due to sub-ice shelf melting and calving of icebergs. Ice sheet models of the Antarctic ice sheet have thus largely concentrated on parameterising sub-shelf and calving processes. However, surface mass balance (SMB) is also of crucial importance in controlling the stability and evolution of the vast Antarctic ice sheet. In this paper we compare the performance of five different regional climate models (COSMO-CLM2, HIRHAM5, MAR3.10, MetUM and RACMO2.3p2) in simulating the near surface climate and SMB of Antarctica. Our results show that, when regional climate models (RCMs) are forced by the ERA-Interim reanalysis, the integrated Antarctic ice sheet ensemble mean annual SMB is 2329 ± 94 Gigatonnes (Gt) year−1 over the common 1987 to 2015 period. However, individual model estimates vary from 1961 ± 70 to 2519 ± 118 Gt year−1. The large differences are mostly explained by different SMB estimates in West Antarctica and the peninsula as well as around the Transantarctic mountains. The calculated annual average SMB is very sensitive to the period chosen but over the climatological mean period of 1980 to 2010 the ensemble mean is 2486 Gt year−1. The interannual variability in SMB is consistent between the models and dominated by variability in the driving ERA-Interim reanalysis. The declining trend in Antarctic SMB reported in other studies is also very sensitive to period chosen and models disagree on the sign and magnitude of the trend in Antarctic SMB over the ERA-Interim period. Evaluation of models shows that they simulate Antarctic climate well when compared with daily observed temperature (Pearson correlation of 0.85 and higher) and pressure (bias ranges from −0.39 hPa in HIRHAM5 to −6.01 hPa in MAR with a mean of −3.49 hPa over all models) and nudged models, constrained within the domain as well as at lateral boundaries, perform better than un-nudged models. We compare modelled surface mass balance with a large dataset of observations which, though biased by undersampling in some regions, indicates that many of the biases in modelled SMB are common between models. The inclusion of drifting snow schemes improves modelled SMB on ice sheet slopes between 1000 and 2000 m where strong katabatic winds form but other regions where precipitation rates are high lack observations needed for the evaluation of different SMB estimates. Different ice masks have a substantial impact on the integrated total SMB and along with model resolution is therefore factored into our analysis. The majority of the different values for continental SMB are due to differences in modelled precipitation at relatively few grid points in coastal areas. Our analysis suggests that targeting coastal areas for observational campaigns will be key to improving and refining estimates of the total surface mass balance of Antarctica.

Snowfall rate (SR) estimates over Antarctica are sparse and characterised by large uncertainties. Yet, observations by precipitation radar offer the potential to get better insight in Antarctic SR. Relations between radar reflectivity (Ze) and snowfall rate (Ze-SR relations) are however not available over Antarctica. Here, we analyse observations from the first Micro Rain Radar (MRR) in Antarctica together with an optical disdrometer (Precipitation Imaging Package; PIP), deployed at the Princess Elisabeth station. The relation Ze = A*SR Bwas derived using PIP observations and its uncertainty was quantified using a bootstrapping approach, randomly sampling within the range of uncertainty.This uncertainty was used to assess the uncertainty in snowfall rates derived