Onset and Energetics of Double-Diffusive Convection in the Ice Shelf/Ocean Boundary Layer

The interaction between Ice-Shelves and the Ocean is an important component of the response of ice sheets to future warming oceans. Observational data in the ocean boundary layer beneath ice shelves is limited and the turbulent flow in the boundary layer is not well characterised. Our work uses small scale (9m depth) direct numerical simulations (DNS) of the Ice-Shelf-Ocean Boundary Layer, inspired by field observations made beneath the George VI Ice Shelf. Here, warm water has been observed directly beneath the ice shelf, and yet the observed melt rates are modest. To study this scenario, we simulate a forced turbulent flow underlying an ice shelf where the ice base is represented by a dynamic melting boundary condition. As the ice melts, a pool of relatively cold, fresh water develops below the ice base. Thermal diffusion causes the underlying water to cool and can drive turbulent convection. At the same time, the salinity gradient in the halocline is stabilising, but develops over a longer time scale. As a result, two flow regimes exist: one with active turbulent convection driven by double-diffusion of heat and salt, and the other with stratified turbulence leading to mixing of the halocline. By varying control parameters, we identify the transition between the flow regimes in terms of the temperature contrast (thermal driving) and the level of turbulence in the far field. We form a ratio of the turbulent velocities associated with the peak in density compared to those associated with the mechanically forced turbulence. The time-dependent behaviour of this ratio is discussed, and it is shown to predict well the regime transition. We explain the process of transition by considering the energetic exchange between the kinetic energy and available potential energy. Finally, we discuss the relevance of the findings to an alternative simulation set up with a shear forcing as well as the application to observations.