On the parameterisation of the effective roughness length for momentum transfer over heterogeneous terrain

We develop a parameterisation for the effective roughness length of terrain that consists of a repeating sequence of patches, in which each patch is composed of strips of two roughness types. A numerical model with second-order closure in the turbulent stress is developed and used to show that: (i) the normalised Reynolds stress develops as a self-similar profile; (ii) the mixing-length parameterisation is a good first-order approximation to the Reynolds stress. These findings are used to characterise the blending layer, where the stress adjusts smoothly from its local surface value to its 'effective' value aloft. Previous studies have assumed that this adjustment occurs abruptly at a single level, often called the blending height. The blending layer is shown to be characterised by height scales that arise naturally in linear models of surface layer flow over roughness changes, and calculations with the numerical model show that these height scales remain appropriate in the nonlinear regime. This concept of the blending layer allows the development of a new parameterisation of the effective roughness length, which gives values for the effective roughness length that are shown to compare well with both atmospheric measurements and values determined from the second-order model.