Abstract The note presents a rational approach to modelling the source/sink due to vegetation or buoyancy effects that appear in the turbulent kinetic energy, E, equation and a supplementary equation for a length-scale determining variable, φ, when two-equation closure is applied to canopy and atmospheric boundary-layer flows. The approach implements only standard model coefficients C φ1 and C φ2 in the production and destruction terms of the φ equation, respectively. Numerical tests illustrate the practical applicability of the method, where, for example, simulations with the E–ω model (where $$φ=ω=ε/E$$ is the specific dissipation and $$ε$$ is the dissipation rate of E) properly reproduce both the surface-layer wind profile estimated from the Monin-Obukhov similarity theory and the mixing-height evolution observed above forested terrain in Southern Finland.