Modelling sea ice and melt ponds evolution

We present a mathematical model describing the evolution of sea ice and meltwater during summer. The system is described by two coupled partial differential equations for the ice thickness h(x,t) and pond depth w(x,t) fields. The model is similar, in principle, to the one put forward by Luthije et al. (2006), but it features i) a modified melting term, ii) a non-uniform seepage rate of meltwater through the porous ice medium and a minimal coupling with the atmosphere via a surface wind shear term, ?s (Scagliarini et al. 2020). We test, in particular, the sensitivity of the model to variations of parameters controlling fluid- dynamic processes at the pond level, namely the variation of turbulent heat flux with pond depth and the lateral melting of ice enclosing a pond. We observe that different heat flux scalings determine different rates of total surface ablations, while the system is relatively robust in terms of probability distributions of pond surface areas. Finally, we study pond morphology in terms of fractal dimensions, showing that the role of lateral melting is minor, whereas there is evidence of an impact from the initial sea ice topography.