Let be the thickness of the thermal boundary layer immediately below , let be the kinematic viscosity, let be the thermal diffusivity, and let be the critical local Rayleigh number. The density contrast driving the boundary layer follows from the density anomaly of water:
A Local Rayleigh-number closure sets
and therefore
By Fourier's law, the upward heat flux through this layer is
where is the water's thermal conductivity. This expression applies while the quantity inside braces is positive.
Vigorous convection keeps the interior nearly uniform, so the upper thermal boundary layer carries a temperature drop of order . A boundary layer becomes unstable when its local Rayleigh number reaches an order-fixed critical value :
This is marginal-stability thermal boundary-layer scaling using the Local Rayleigh-number closure: boundary-layer buoyancy overcomes viscous and thermal-diffusive stabilization. Equivalently, a buoyancy-driven velocity scale makes order . Thus
There is no in this closure, as required for a thin boundary layer controlled by its own stability. The exponent is fixed by scaling; its absolute coefficient is not. Here represents the geometry and boundary conditions, together with any order-one boundary-layer coefficient. The question supplies no numerical value for it. The assumption requires , positive , and properties sufficiently constant over the temperature contrast.
In an ice-covered freshwater lake, a stagnant layer colder than the density-maximum temperature can cap a warmer turbulent region. Maximizing the boundary-layer heat flux obtained from a Local Rayleigh-number closure gives an interface temperature satisfying
and a maximum upward heat flux proportional to .