A wind-driven snow redistribution module for Alpine3D v3.3.0: adaptations designed for downscaling ice sheet surface mass balance

Ice sheet surface mass balance describes the net snow accumulation at the ice sheet surface. On the Antarctic ice sheet, winds redistribute snow, resulting in a surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 inline-formulakm. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 inline-formula $M2inlinescrollmathmlunit\mathrm{normal\; mm}\mathrm{normal\; w}.\mathrm{normal\; e}.{\mathrm{normal\; yr}}^{-\mathrm{normal\; 1}}$ 61pt15ptsvg-formulamathimgababe4b727a34893f1cecbd88151fb23 gmd-16-3203-2023-ie00001.svg61pt15ptgmd-16-3203-2023-ie00001.png (13 %) and increases variance explained by 24 %. Despite these improvements, our results demonstrate that considerable uncertainty stems from the employed saltation model, confounding simulations of surface mass balance variability.