The interaction between deep oceanic currents and an ice base is critical to accurately predict global ice melting rates, yet predictions are often affected by inaccuracies due to inadequate dynamical modelling of the ice–water interface morphology. To improve current predictive models, we numerically investigate the evolution of the ice–water interface under a subsurface turbulent shear-dominated flow, focusing on the time and length scales that govern both global and local morphological features. Based on our previous work (Perissutti, Marchioli & Soldati 2024 Intl J. Multiphase Flow 181, 105007), where we confirmed the existence of a threshold Reynolds number below which only streamwise-oriented topography forms and above which a larger-scale spanwise topography emerges and coexists with the streamwise structures, we explore three orders of magnitude for the Stefan number (the ratio of sensible heat to latent heat). We examine its impact on ice melting and its role in shaping the interface across the two distinct morphodynamic regimes. We identify characteristic time scales of ice melting and demonstrate that the key features of ice morphodynamics scale consistently with the Stefan number and the Péclet number (the ratio of heat advection to diffusion) in both regimes. These scaling relationships can be leveraged to infer the main morphodynamic characteristics of the ice–water interface from direct numerical simulation datasets generated at computationally feasible values of Péclet and Stefan numbers, enabling the incorporation of morphodynamics into geophysical melting models and thereby enhancing their predictive accuracy.
