Estimation of the Antarctic surface mass balance using MAR (1979–2015) and identification of dominant processes

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 surface mass balance 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 MERRA2, for the period 1979–2015, and we compare our 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 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 erosion-deposition processes not included in MAR, where the drifting snow module has been switched off, and probably underestimated in RACMO2 by a factor of three. As a consequence, the amount of drifting snow sublimating in the atmospheric boundary layer remains a potentially large mass sink needed to be better constrained. Moreover, MAR generally simulates larger SMB and snowfall amounts than RACMO2 inland, whereas snowfall rates are significantly lower in MAR than in RACMO2 at the ice sheet margins. This divergent behaviour at the margins results from differences in model parameterisations, as MAR explicitly advects precipitating particles through the atmospheric layers and sublimates snowflakes in the undersaturated katabatic layer, whereas in RACMO2 precipitation is added to the surface without advection through the atmosphere. Consequently, we corroborate a recent study concluding that sublimation of precipitation in the low-level atmospheric layers is a significant mass sink for the Antarctic SMB, as it may represent ∼ 240 ± 25 Gt yr-1 of difference in snowfall between RACMO2 and MAR for the period 1979–2015, which is 10 % of the simulated snowfall loaded on the ice sheet and more than twice the surface snow sublimation as currently simulated by MAR.