Spatially-resolved mean flow and turbulence help explain observed erosion and deposition patterns of snow over Antarctic sea ice

Sea ice-atmosphere interactions are major drivers of patterns of sea ice flows and deformations in Polar regions, and affect snow erosion and deposition at the surface. Here, we combine analyses of sea ice surface topography at very high-resolutions (~1-10 cm), and Large Eddy Simulations (LES) to study surface drag and snow erosion and deposition patterns from process scales to floe scales (~ 1 cm – 100 m). The snow/ice elevations were obtained using a Terrestrial Laser Scanner during the SIPEX II (Sea Ice Physics and Ecosystem eXperiment II) research voyage to East Antarctica (September-November 2012). LES are performed on a regular domain adopting a mixed pseudo-spectral/finite difference spatial discretization. A scale-dependent dynamic subgrid-scale model based on Lagrangian time averaging is adopted to determine the eddy-viscosity in the bulk of the flow. Effects of larger-scale features of the surface on wind flows (those features that can be resolved in the LES) are accounted for through an immersed boundary method. Conversely, drag forces caused by subgrid-scale features of the surface should be accounted for through a parameterization. However, the effective aerodynamic roughness parameter z0 for snow/ice is not known. Hence, a novel dynamic approach is utilized, in which z0 is determined using the constraint that the total momentum flux (drag) must be independent on grid-filter scale. We focus on three ice floe surfaces. The first of these surfaces (October 6, 2012) is used to test the performance of the model, validate the algorithm, and study the spatial distributed fields of resolved and modeled stress components. The following two surfaces, scanned at the same location before and after a snow storm event (October 20 and 23, 2012), are used to propose an application to study how spatially resolved mean flow and turbulence relates to observed patterns of snow erosion and deposition. We show how erosion and deposition patterns are correlated with the computed stresses, with modeled stresses having higher explanatory power. Deposition is mainly occurring in wake regions of specific ridges that strongly affect wind flow patterns. These larger ridges lock in place elongated streaks of relatively high speeds with axes along the stream-wise direction, and which are largely responsible for the observed erosion.