Sub-kilometer Scale Snow Depth Distribution on Sea Ice of Different Ages and Thickness

Abstract. Accurately representing the snow depth (SND) distribution on sea ice is essential for sea ice thickness (SIT) retrievals, ecological studies, and climate modeling. Using co-located SND and SIT measurements from multiple Arctic and Antarctic campaigns, this study examines sub-kilometer-scale SND variability, considering both ice type and SIT, and identifies the most suitable statistical distributions to represent SND across different ice ages and thicknesses. First, we examine the statistical properties of SND and their dependence on SIT, finding a linear increase of SND with SIT for new and first-year ice, reflecting concurrent seasonal growth. The ratio between the standard deviation and the mean SND is referred to as the coefficient of variation (CV). A consistent CV ≈ 0.50 is observed to be independent of SIT, allowing variability to be estimated directly from the mean SND. Notably, flooded snow exhibits a lower CV. Furthermore, we investigate four probability density functions (Normal, Log-normal, Gamma, and Skew) and find that the best-fit distribution depends on ice ages, SIT, deformation, and meteorological events such as snow fall and drift. Finally, SND correlation lengths derived from semi-variograms show a positive relation with SIT and are enhanced by snow drift events. The results reveal substantial differences in SND distributions across ice types and SIT during winter and summer, underscoring the importance of ice-condition-dependent parameterizations for representing sub-kilometer SND variability. These findings support improved parameterizations of SND variability at sub-grid scale in remote sensing and climate models.