Meltwater transport and mixing-layer growth near the ice–ocean interface
Abstract
Ice melting into saline water plays a fundamental role in the dynamics near the ice–ocean interface in polar oceans. The physics of meltwater transport involves a non-trivial interplay between the thermodynamics at the interface, hydrodynamic transport in the bulk and the properties of the ambient ocean. The key control parameters are the density ratio Rρ, which is proportional to the ambient salinity and measures the balance between the temperature and salinity effects on density, together with the Lewis number L = κT/κS, which compares thermal and solutal diffusivities. In quiescent horizontal configurations, increasing the salinity is known to slow down melting, with the melt rate transitioning from subdiffusive to diffusive as upper Rρ increases. Here, we assess the role of turbulence in this transition, using highly resolved numerical simulations of the two-dimensional Boussinesq equations with a slowly melting upper boundary. We analyse the non-stationary growth of the thermal and solutal mixing layers, varying the Lewis number and the density ratio. While meltwater transport is continuously driven by convection within the bulk, we identify a transition from convection to diffusion close to the interface. This transition is reflected by the formation of an interfacial boundary layer that regulates the flux of meltwater pouring into the turbulent bulk. The boundary layer has little effect on thermal transport, but strongly suppresses solutal transport for upper Rρ ≳ 10, only allowing a fraction ∝ Rρ−1 of the input flux to reach the bulk. Using mixing-layer diagnostics based on solutal-concentration thresholds, we observe that the turbulent layer grows super-diffusively ∝ t , while the interfacial boundary layer expands diffusively ∝ t0.5 but with a non-universal prefactor depending on the Lewis number and density ratio. The superdiffusive growth of the mixing layer challenges the commonly assumed picture of a fully diffusive regime at high salinity. Overall, our results indicate that double-diffusive effects are here confined to the interface, and highlight potential limitations of diagnostics based on fixed concentration thresholds in oceanographic applications.
Related
- Influence of irradiance and mixing layer height on the vertical trace matter distribution in the lower planetary boundary layer – drone-based investigation
- Precession and meltwater control of greenhouse gas emissions from terrestrial sources during the last glacial period
- Vertical transport and confinement of weakly buoyant particles in a convective ocean mixed-layer model
- From prototypical regimes to observed complex dynamics and CO ₂ transport in the Amazonian nocturnal boundary layer
- Experimental mineralogical replacement to investigate the role of texture on interface morphology, reaction kinetics and pore-scale transport during metasomatism
BibTeX
@article{allende2026meltwater,
title = {{Meltwater transport and mixing-layer growth near the ice–ocean interface}},
author = {Sofia Allende and Louis-Alexandre Couston and Simon Thalabard and Benjamin Favier},
journal = {Journal of Fluid Mechanics},
year = {2026},
month = aug,
doi = {10.1017/jfm.2026.11887},
eprint = {2601.18674},
archivePrefix = {arXiv},
url = {https://hal.science/hal-05752582v1},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Meltwater transport and mixing-layer growth near the ice–ocean interface
AU - Sofia Allende
AU - Louis-Alexandre Couston
AU - Simon Thalabard
AU - Benjamin Favier
PY - 2026
DA - 2026/08/07
JO - Journal of Fluid Mechanics
DO - 10.1017/jfm.2026.11887
UR - https://hal.science/hal-05752582v1
AB - Ice melting into saline water plays a fundamental role in the dynamics near the ice–ocean interface in polar oceans. The physics of meltwater transport involves a non-trivial interplay between the thermodynamics at the interface, hydrodynamic transport in the bulk and the properties of the ambient ocean. The key control parameters are the density ratio Rρ, which is proportional to the ambient salinity and measures the balance between the temperature and salinity effects on density, together with the Lewis number L = κT/κS, which compares thermal and solutal diffusivities. In quiescent horizontal configurations, increasing the salinity is known to slow down melting, with the melt rate transitioning from subdiffusive to diffusive as upper Rρ increases. Here, we assess the role of turbulence in this transition, using highly resolved numerical simulations of the two-dimensional Boussinesq equations with a slowly melting upper boundary. We analyse the non-stationary growth of the thermal and solutal mixing layers, varying the Lewis number and the density ratio. While meltwater transport is continuously driven by convection within the bulk, we identify a transition from convection to diffusion close to the interface. This transition is reflected by the formation of an interfacial boundary layer that regulates the flux of meltwater pouring into the turbulent bulk. The boundary layer has little effect on thermal transport, but strongly suppresses solutal transport for upper Rρ ≳ 10, only allowing a fraction ∝ R_ρ^−1 of the input flux to reach the bulk. Using mixing-layer diagnostics based on solutal-concentration thresholds, we observe that the turbulent layer grows super-diffusively ∝ t , while the interfacial boundary layer expands diffusively ∝ t0.5 but with a non-universal prefactor depending on the Lewis number and density ratio. The superdiffusive growth of the mixing layer challenges the commonly assumed picture of a fully diffusive regime at high salinity. Overall, our results indicate that double-diffusive effects are here confined to the interface, and highlight potential limitations of diagnostics based on fixed concentration thresholds in oceanographic applications.
SN - 0022-1120
SN - 1469-7645
ER -