Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Velocity Adjustment and Passive Scalar Diffusion in and Above an Urban Canopy in Response to Various Approach Flows

We used wind-tunnel experiments to investigate velocity-field adjustment and scalar diffusion behaviour in and above urban canopies located downwind of various roughness elements. Staggered arrays of rectangular blocks of various heights H and plan area ratios λ_p were used to model the urban canopies. The velocity field in the roughness sublayer (height \({z \lesssim 2H}\)) reached equilibrium at distances proportional to \({\sqrt{L_{\rm c}H}}\) where L _c is the canopy-drag length scale determined as a function of λ_p and the block side length L. A distance of about \({20\sqrt{L_{\rm c}H}}\) was required for adjustment at z = H/2 (in the canopy), and a distance of about \({10\sqrt{L_{\rm c}H}}\) was required at z = 2H (near the top of the roughness sublayer). Diffusion experiments from a ground emission source revealed that differences in upwind roughness conditions had negligible effects on the plume growth near the source (up to a few multiples of L from the source) if the source was located at a fetch F larger than about \({10\sqrt{L_{\rm c}H}}\) from the upwind edge of the canopy. However, at locations farther downwind (more than several multiples of L from the source), upwind conditions had considerable effects on the plume growth. For a representative urban canopy, it was shown that a much larger fetch than required for velocity-field adjustment in the roughness sublayer was necessary to eliminate the effects of upwind conditions on plume widths at 24L downwind from the source.     

Turbulent Structure of a Simplified Urban Fluid Flow Studied Through Stereoscopic Particle Image Velocimetry

Stereoscopic particle image velocimetry was used to provide a three-dimensional characterization of the flow around a simplified urban model defined by a 5 by 7 array of blocks, forming four parallel streets, perpendicular to the incoming wind direction corresponding to a zero angle of incidence. Channeling of the flow through the array under consideration was observed, and its effect increased as the incoming wind direction, or angle of incidence (AOI), was changed from \(0^{\circ }\) to \(15^{\circ }\), \(30^{\circ }\), and \(45^{\circ }\). The flow between blocks can be divided into two regions: a region of low turbulence kinetic energy (TKE) levels close to the leeward side of the upstream block, and a high TKE area close to the downstream block. The centre of the arch vortex is located in the low TKE area, and two regions of large streamwise velocity fluctuation bound the vortex in the spanwise direction. Moreover, a region of large spanwise velocity fluctuation on the downstream block is found between the vortex legs. Our results indicate that the reorientation of the arch vortex at increasing AOI is produced by the displacement of the different TKE regions and their interaction with the shear layers on the sides and top of the upstream and downstream blocks, respectively. There is also a close connection between the turbulent structure between the blocks and the wind gusts. The correlations among gust components were also studied, and it was found that in the near-wall region of the street the correlations between the streamwise and spanwise gusts \(R_{uv}\) were dominant for all four AOI cases. At higher wall-normal positions in the array, the \(R_{uw}\) correlation decreased with increasing AOI, whereas the \(R_{uv}\) coefficient increased as AOI increased, and at \({\textit{AOI}}=45^{\circ }\) all three correlations exhibited relatively high values of around 0.4.     

Edge-to-Stem Variability in Wet-Canopy Evaporation From an Urban Tree Row

Evaporation from wet-canopy (\(E_\mathrm{C}\)) and stem (\(E_\mathrm{S}\)) surfaces during rainfall represents a significant portion of municipal-to-global scale hydrologic cycles. For urban ecosystems, \(E_\mathrm{C}\) and \(E_\mathrm{S}\) dynamics play valuable roles in stormwater management. Despite this, canopy-interception loss studies typically ignore crown-scale variability in \(E_\mathrm{C}\) and assume (with few indirect data) that \(E_\mathrm{S}\) is generally \({<}2\%\) of total wet-canopy evaporation. We test these common assumptions for the first time with a spatially-distributed network of in-canopy meteorological monitoring and 45 surface temperature sensors in an urban Pinus elliottii tree row to estimate \(E_\mathrm{C}\) and \(E_\mathrm{S}\) under the assumption that crown surfaces behave as “wet bulbs”. From December 2015 through July 2016, 33 saturated crown periods (195 h of 5-min observations) were isolated from storms for determination of 5-min evaporation rates ranging from negligible to 0.67 \(\hbox {mm h}^{-1}\). Mean \(E_\mathrm{S}\) (0.10 \(\hbox {mm h}^{-1}\)) was significantly lower (\(p < 0.01\)) than mean \(E_\mathrm{C}\) (0.16 \(\hbox {mm h}^{-1}\)). But, \(E_\mathrm{S}\) values often equalled \(E_\mathrm{C}\) and, when scaled to trunk area using terrestrial lidar, accounted for 8–13% (inter-quartile range) of total wet-crown evaporation (\(E_\mathrm{S}+E_\mathrm{C}\) scaled to surface area). \(E_\mathrm{S}\) contributions to total wet-crown evaporation maximized at 33%, showing a general underestimate (by 2–17 times) of this quantity in the literature. Moreover, results suggest wet-crown evaporation from urban tree rows can be adequately estimated by simply assuming saturated tree surfaces behave as wet bulbs, avoiding problematic assumptions associated with other physically-based methods.     

Test of Scintillometer Saturation Correction Methods Using Field Experimental Data

Saturation of large aperture scintillometer (LAS) signals can result in sensible heat flux measurements that are biased low. A field study with LASs of different aperture sizes and path lengths was performed to investigate the onset of, and corrections for, signal saturation. Saturation already occurs at \({C_n^2 \approx 0.074 D^{5/3} \lambda^{1/3} L^{-8/3}}\), where \({C_n^2}\) is the structure parameter of the refractive index, D is the aperture size, λ is the wavelength, L is the transect length, which is smaller than theoretically derived saturation limits. At a transect length of 1 km, a height of 2.5 m, and aperture ≈0.15 m the correction factor exceeds 5% already at \({C_n^2=2\times 10^{-12}{\rm m}^{-2/3}}\), which will affect many practical applications of scintillometry. The Clifford correction method, which only depends on \({C_n^2}\) and the transect geometry, provides good saturation corrections over the range of conditions observed in our study. The saturation correction proposed by Ochs and Hill results in correction factors that are too small in large saturation regimes. An inner length scale dependence of the saturation correction factor was not observed. Thus for practical applications the Clifford correction method should be applied.     

On the Relevance of $$\ln (z_{0} /z_{0T})= kB^{-1}$$

The assumption that the roughness Reynolds number \(( Re_{*})\) can be used as a basis for quantifying the boundary-layer property \({ kB}^{-1} (= \ln (z_{0}/z_{0T}))\) as in some modern numerical models is questioned. While \({ Re}_{*}\) is a useful property in studies of pipe flow, it appears to have only marginal applicability in the case of treeless terrain, as studied in the two experimental situations presented here. For both the daytime and night-time cases there appears to be little correlation between \({ kB}^{-1}\) and \({ Re}_{*}\). For daytime, the present studies indicate that the assumption \({ kB}^{-1} \approx 2\) is acceptable, while for night-time, the scatter involved in relating \({ kB}^{-1}\) to \({ Re}_{*}\) suggests there is little reason to assume a direct relationship. However, while the scatter affecting all of the night-time results is large, there remains a significant correlation between the heat and momentum fluxes upon which an alternative methodology for describing bulk air–surface exchange at night could be constructed. The friction coefficient (\(C_{f}\)) and the turbulent Stanton number \(({ St}_{*})\) are discussed as possible alternatives for describing bulk properties of the air layer adjacent to the surface. While describing the surface roughness in terms of the friction coefficient provides an attractive simplification relative to the conventional methodologies based on roughness length and stability considerations, use of the Stanton number shares many of uncertainties that affect \({ kB}^{-1}\). The transitions at dawn and dusk remain demanding situations to address.     

Katabatic Flow: A Closed-Form Solution with Spatially-Varying Eddy Diffusivities

The Nieuwstadt closed-form solution for the stationary Ekman layer is generalized for katabatic flows within the conceptual framework of the Prandtl model. The proposed solution is valid for spatially-varying eddy viscosity and diffusivity (O’Brien type) and constant Prandtl number (Pr). Variations in the velocity and buoyancy profiles are discussed as a function of the dimensionless model parameters \(z_0 \equiv \hat{z}_0 \hat{N}^2 Pr \sin {(\alpha )} |\hat{b}_\mathrm{s} |^{-1}\) and \(\lambda \equiv \hat{u}_{\mathrm{ref}}\hat{N} \sqrt{Pr} |\hat{b}_\mathrm{s} |^{-1}\), where \(\hat{z}_0\) is the hydrodynamic roughness length, \(\hat{N}\) is the Brunt-Väisälä frequency, \(\alpha \) is the surface sloping angle, \(\hat{b}_\mathrm{s}\) is the imposed surface buoyancy, and \(\hat{u}_{\mathrm{ref}}\) is a reference velocity scale used to define eddy diffusivities. Velocity and buoyancy profiles show significant variations in both phase and amplitude of extrema with respect to the classic constant \(\textit{K}\) model and with respect to a recent approximate analytic solution based on the Wentzel-Kramers-Brillouin theory. Near-wall regions are characterized by relatively stronger surface momentum and buoyancy gradients, whose magnitude is proportional to \(z_0\) and to \(\lambda \). In addition, slope-parallel momentum and buoyancy fluxes are reduced, the low-level jet is further displaced toward the wall, and its peak velocity depends on both \(z_0\) and \(\lambda \).     

Observations of Atmospheric Methane and Carbon Dioxide Mixing Ratios: Tall-Tower or Mountain-Top Stations?

Mountain-top observations of greenhouse gas mixing ratios may be an alternative to tall-tower measurements for regional scale source and sink estimation. To investigate the equivalence or limitations of a mountain-top site as compared to a tall-tower site, we used the unique opportunity of comparing in situ measurements of methane (\(\hbox {CH}_{4}\)) and carbon dioxide (\(\hbox {CO}_{2}\)) mixing ratios at a mountain top (986 m above sea level, a.s.l.) with measurements from a nearby (distance 28.4 km) tall tower, sampled at almost the same elevation (1009 m a.s.l.). Special attention was given to, (i) how local wind statistics and greenhouse gas sources and sinks at the mountain top influence the observations, and (ii) whether mountain-top observations can be used as for those from a tall tower for constraining regional greenhouse gas emissions. Wind statistics at the mountain-top site are clearly more influenced by local flow systems than those at the tall-tower site. Differences in temporal patterns of the greenhouse gas mixing ratios observed at the two sites are mostly related to the influence of local sources and sinks at the mountain-top site. Major influences of local sources can be removed by applying a statistical filter (\(5{\mathrm{th}}\) percentile) or a filter that removes periods with unfavourable flow conditions. In the best case, the bias in mixing ratios between the mountain-top and the tall-tower sites after the application of the wind filter was \({-}0.0005\pm 0.0010\) ppm for methane (September, 0000–0400 UTC) and \(0.11\pm 0.18\) ppm for \(\hbox {CO}_{2}\) (February, 1200–1600 UTC). Temporal fluctuations of atmospheric \(\hbox {CH}_{4}\) and \(\hbox {CO}_{2}\) mixing ratios at both stations also showed good agreement (apart from \(\hbox {CO}_{2}\) during summertime) as determined by moving bi-weekly Pearson correlation coefficients (up to 0.96 for \(\hbox {CO}_{2}\) and 0.97 for \(\hbox {CH}_{4}\)). When only comparing mixing ratios minimally influenced by local sources (low bias and high correlation coefficients), our measurements indicate that mountain-top observations are comparable to tall-tower observations.     

Water-Channel Estimation of Eulerian and Lagrangian Time Scales of the Turbulence in Idealized Two-Dimensional Urban Canopies

Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian \((T^{\mathrm{E}})\) and Lagrangian \((T^{\mathrm{L}})\) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width (W) to the height (H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, \(T_\mathrm{u}^\mathrm{L} \) and \(T_\mathrm{w}^\mathrm{L} \), follow Raupach’s linear law within the constant-flux layer. The same holds true for \(T_\mathrm{w}^\mathrm{L} \) in both the canopies analyzed \((AR= 1\) and \(AR= 2\)) and also for \(T_\mathrm{u}^\mathrm{L} \) when \(AR = 1\). In contrast, for \(AR = 2\), \(T_\mathrm{u}^\mathrm{L} \) follows Raupach’s law only above \(z=2H\). Below that level, \(T_\mathrm{u}^\mathrm{L} \) is nearly constant with height, showing at \(z=H\) a value approximately one order of magnitude greater than that found for \(AR = 1\). It is shown that the assumption usually adopted for flat terrain, that \(\beta =T^{\mathrm{L}}/T^{\mathrm{E}}\) is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, \(\gamma /i_\mathrm{u} \) fits well \(\beta _\mathrm{u} =T_\mathrm{u}^\mathrm{L} /T_\mathrm{u}^\mathrm{E} \) in both the configurations by choosing \(\gamma \) to be 0.35 (here, \(i_\mathrm{u} =\sigma _\mathrm{u} / \bar{u} \), where \(\bar{u} \) and \(\sigma _\mathrm{u} \) are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, \(\beta _\mathrm{w} =T_\mathrm{w}^\mathrm{L} /T_\mathrm{w}^\mathrm{E} \) follows approximately \(\gamma /i_\mathrm{w} =0.65/\left( {\sigma _\mathrm{w} /\bar{u} } \right) \) for \(z > 2H\), irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum \((K_\mathrm{{T}})\) and the Kolmogorov constant \((C_0)\). It is found that \(C_0\) depends appreciably on the velocity component both for the flat terrain and canopy flow, even though for the latter case it is insensitive to AR values. In all the three experimental configurations analyzed here, \(K_\mathrm{{T}}\) shows an overall linear growth with height in agreement with the linear trend predicted by Prandtl’s theory.     

Spatially-Averaged Temperature Structure Parameter Over a Heterogeneous Surface Measured by an Unmanned Aerial Vehicle

The structure parameter of temperature, \({C_{T}^{2}}\) , in the lower convective boundary layer was measured using the unmanned mini aerial vehicle M²AV. The measurements were carried out on two hot summer days in July 2010 over a heterogeneous land surface around the boundary-layer field site of the Lindenberg Meteorological Observatory—Richard-Aßmann-Observatory of the German Meteorological Service. The spatial series of \({C_{T}^{2}}\) showed considerable variability along the flight path that was caused by both temporal variations and surface heterogeneity. Comparison of the aircraft data with \({C_{T}^{2}}\) values derived from tower-based in situ turbulence measurements showed good agreement with respect to the diurnal variability. The decrease of \({C_{T}^{2}}\) with height as predicted by free-convection scaling could be confirmed for the morning and afternoon flights while the flights around noon suggest a different behaviour.     

Scale Properties of Anisotropic and Isotropic Turbulence in the Urban Surface Layer

The scale properties of anisotropic and isotropic turbulence in the urban surface layer are investigated. A dimensionless anisotropic tensor is introduced and the turbulent tensor anisotropic coefficient, defined as C, where \(C = 3d_{3}\,+\,1 (d_{3}\) is the minimum eigenvalue of the tensor) is used to characterize the turbulence anisotropy or isotropy. Turbulence is isotropic when \(C \approx 1\), and anisotropic when \(C \ll 1\). Three-dimensional velocity data collected using a sonic anemometer are analyzed to obtain the anisotropic characteristics of atmospheric turbulence in the urban surface layer, and the tensor anisotropic coefficient of turbulent eddies at different spatial scales calculated. The analysis shows that C is strongly dependent on atmospheric stability \(\xi = (z-z_{\mathrm{d}})/L_{{\textit{MO}}}\), where z is the measurement height, \(z_{\mathrm{d}}\) is the displacement height, and \(L_{{\textit{MO}}}\) is the Obukhov length. The turbulence at a specific scale in unstable conditions (i.e., \(\xi < 0\)) is closer to isotropic than that at the same scale under stable conditions. The maximum isotropic scale of turbulence is determined based on the characteristics of the power spectrum in three directions. Turbulence does not behave isotropically when the eddy scale is greater than the maximum isotropic scale, whereas it is horizontally isotropic at relatively large scales. The maximum isotropic scale of turbulence is compared to the outer scale of temperature, which is obtained by fitting the temperature fluctuation spectrum using the von Karman turbulent model. The results show that the outer scale of temperature is greater than the maximum isotropic scale of turbulence.     

Observations of the Temperature and Humidity Structure Parameter Over Heterogeneous Terrain by Airborne Measurements During the LITFASS-2003 Campaign

We used numerical simulations to investigate the general relationship between urban morphology and the intensity of wind gusts in built-up areas at the pedestrian level. The simulated urban boundary layer developed over a 19.2 km (length) \(\times \) 4.8 km (width) \(\times \) 1.0 km (height) simulation domain, with 2-m resolution in all directions, to explicitly resolve the detailed shapes of buildings and the flow at the pedestrian level. This complex computation was accomplished using the lattice Boltzmann method and by implementing a large-eddy simulation model. To generalize the results, a new parameter that expresses the intensity of gusts (the gust index, \({\tilde{U}}_{ max})\) was defined as the local maximum wind speed divided by the freestream velocity. In addition, this parameter was decomposed into the mean wind-speed ratio, \({\tilde{U}} \) and turbulent gust ratio, \({\tilde{U}}^{{\prime }}\) to evaluate the qualities of gusts. These parameters were useful for quantitatively comparing the gust intensities within urban canopies at different locations or even among different experiments. In addition, the entire horizontal domain was subdivided into homogeneous square patches, in which both the simulated gust parameters and the morphological characteristics of building geometries were averaged. This procedure masked the detailed structure of individual buildings but retained the bulk characteristics of the urban morphology. At the pedestrian level, the gust index decreased with increasing building cover. Compared to \({\tilde{U}} \), the quantity \({\tilde{U}}^{{\prime }}\) notably contributed to the index throughout the range of plan area index \((\lambda _p)\) values. The dependences of all normalized wind-speed ratios transiently changed at \(\lambda _p =~0.28\). In cases where \(\lambda _p < 0.28, {\tilde{U}} \) decreased with increasing \(\lambda _p \), although \({\tilde{U}}^{{\prime }}\) was almost constant. In cases where \(\lambda _p > 0.28, {\tilde{U}}\) was almost constant and \({\tilde{U}}^{{\prime }}\) decreased with increasing \(\lambda _p \). This was explained by the change in flow regimes within the building canyon. At a higher elevation above the canopy layer, \(\lambda _p \) becomes less relevant to normalized wind-speed ratios, and instead the aerodynamic roughness length became important.     

Evaluation of Urban Local-Scale Aerodynamic Parameters: Implications for the Vertical Profile of Wind Speed and for Source Areas

Nine methods to determine local-scale aerodynamic roughness length \((z_{0})\) and zero-plane displacement \((z_{d})\) are compared at three sites (within 60 m of each other) in London, UK. Methods include three anemometric (single-level high frequency observations), six morphometric (surface geometry) and one reference-based approach (look-up tables). A footprint model is used with the morphometric methods in an iterative procedure. The results are insensitive to the initial \(z_{d}\) and \(z_{0}\) estimates. Across the three sites, \(z_{d}\) varies between 5 and 45 m depending upon the method used. Morphometric methods that incorporate roughness-element height variability agree better with anemometric methods, indicating \(z_{d}\) is consistently greater than the local mean building height. Depending upon method and wind direction, \(z_{0}\) varies between 0.1 and 5 m with morphometric \(z_{0}\) consistently being 2–3 m larger than the anemometric \(z_{0}\). No morphometric method consistently resembles the anemometric methods. Wind-speed profiles observed with Doppler lidar provide additional data with which to assess the methods. Locally determined roughness parameters are used to extrapolate wind-speed profiles to a height roughly 200 m above the canopy. Wind-speed profiles extrapolated based on morphometric methods that account for roughness-element height variability are most similar to observations. The extent of the modelled source area for measurements varies by up to a factor of three, depending upon the morphometric method used to determine \(z_{d}\) and \(z_{0}\).     

Momentum- and Heat-Flux Parametrization at Dome C,Antarctica: A Sensitivity Study

An extensive meteorological observational dataset at Dome C, East Antarctic Plateau, enabled estimation of the sensitivity of surface momentum and sensible heat fluxes to aerodynamic roughness length and atmospheric stability in this region. Our study reveals that (1) because of the preferential orientation of snow micro-reliefs (sastrugi), the aerodynamic roughness length \(z_{0}\) varies by more than two orders of magnitude depending on the wind direction; consequently, estimating the turbulent fluxes with a realistic but constant \(z_{0}\) of 1 mm leads to a mean friction velocity bias of \(24\,\%\) in near-neutral conditions; (2) the dependence of the ratio of the roughness length for heat \(z_{0t}\) to \(z_{0}\) on the roughness Reynolds number is shown to be in reasonable agreement with previous models; (3) the wide range of atmospheric stability at Dome C makes the flux very sensitive to the choice of the stability functions; stability function models presumed to be suitable for stable conditions were evaluated and shown to generally underestimate the dimensionless vertical temperature gradient; as these models differ increasingly with increases in the stability parameter z / L, heat flux and friction velocity relative differences reached \(100\,\%\) when \(z/L > 1\); (4) the shallowness of the stable boundary layer is responsible for significant sensitivity to the height of the observed temperature and wind data used to estimate the fluxes. Consistent flux results were obtained with atmospheric measurements at heights up to 2 m. Our sensitivity study revealed the need to include a dynamical parametrization of roughness length over Antarctica in climate models and to develop new parametrizations of the surface fluxes in very stable conditions, accounting, for instance, for the divergence in both radiative and turbulent fluxes in the first few metres of the boundary layer.     

Aerodynamic Properties of Rough Surfaces with High Aspect-Ratio Roughness Elements: Effect of Aspect Ratio and Arrangements

We examine the effect of varying roughness-element aspect ratio on the mean velocity distributions of turbulent flow over arrays of rectangular-prism-shaped elements. Large-eddy simulations (LES) in conjunction with a sharp-interface immersed boundary method are used to simulate spatially-growing turbulent boundary layers over these rough surfaces. Arrays of aligned and staggered rectangular roughness elements with aspect ratio >1 are considered. First the temporally- and spatially-averaged velocity profiles are used to illustrate the aspect-ratio effects. For aligned prisms, the roughness length (\(z_\mathrm{o}\)) and the friction velocity (\(u_*\)) increase initially with an increase in the roughness-element aspect ratio, until the values reach a plateau at a particular aspect ratio. The exact value of this aspect ratio depends on the coverage density. Further increase in the aspect ratio changes neither \(z_\mathrm{o}\), \(u_*\) nor the bulk flow above the roughness elements. For the staggered cases, \(z_\mathrm{o}\) and \(u_*\) continue to increase for the surface coverage density and the aspect ratios investigated. To model the flow response to variations in roughness aspect ratio, we turn to a previously developed phenomenological volumetric sheltering model (Yang et al., in J Fluid Mech 789:127–165, 2016), which was intended for low to moderate aspect-ratio roughness elements. Here, we extend this model to account for high aspect-ratio roughness elements. We find that for aligned cases, the model predicts strong mutual sheltering among the roughness elements, while the effect is much weaker for staggered cases. The model-predicted \(z_\mathrm{o}\) and \(u_*\) agree well with the LES results. Results show that the model, which takes explicit account of the mutual sheltering effects, provides a rapid and reliable prediction method of roughness effects in turbulent boundary-layer flows over arrays of rectangular-prism roughness elements.     

Optimizing the Determination of Roughness Parameters for Model Urban Canopies

We present an objective optimization procedure to determine the roughness parameters for very rough boundary-layer flow over model urban canopies. For neutral stratification the mean velocity profile above a model urban canopy is described by the logarithmic law together with the set of roughness parameters of displacement height d, roughness length \(z_0\), and friction velocity \(u_*\). Traditionally, values of these roughness parameters are obtained by fitting the logarithmic law through (all) the data points comprising the velocity profile. The new procedure generates unique velocity profiles from subsets or combinations of the data points of the original velocity profile, after which all possible profiles are examined. Each of the generated profiles is fitted to the logarithmic law for a sequence of values of d, with the representative value of d obtained from the minima of the summed least-squares errors for all the generated profiles. The representative values of \(z_0\) and \(u_*\) are identified by the peak in the bivariate histogram of \(z_0\) and \(u_*\). The methodology has been verified against laboratory datasets of flow above model urban canopies.     

Understanding Diurnality and Inter-Seasonality of a Sub-tropical Urban Heat Island

We quantify the spatial and temporal aspects of the urban heat-island (UHI) effect for Kanpur, a major city in the humid sub-tropical monsoon climate of the Gangetic basin. Fixed station measurements are used to investigate the diurnality and inter-seasonality in the urban–rural differences in surface temperature (\({\Delta } T_\mathrm{s}\)) and air temperature (\({\Delta } T_\mathrm{c}\)) separately. The extent of the spatial variations of the nighttime \({\Delta } T_\mathrm{c}\) and \({\Delta } T_\mathrm{s}\) is investigated through mobile campaigns and satellite remote sensing respectively. Nighttime \({\Delta } T_\mathrm{c}\) values dominate during both the pre-monsoon (maximum of 3.6 \(^\circ \hbox {C}\)) and the monsoon (maximum of 2.0 \(^\circ \hbox {C}\)). However, the diurnality in \({\Delta } T_\mathrm{s}\) is different, with higher daytime values during the pre-monsoon, but very little diurnality during the monsoon. The nighttime \({\Delta } T_\mathrm{s}\) value is mainly associated with differences in the urban–rural incoming longwave radiative flux (\(r^{2}=0.33\) during the pre-monsoon; 0.65 during the monsoon), which, in turn, causes a difference in the outgoing longwave radiative flux. This difference may modulate the nighttime \({\Delta } T_\mathrm{c}\) value as suggested by significant correlations (\(r^{2}=0.68\) for the pre-monsoon; 0.50 for the monsoon). The magnitude of \({\Delta } T_\mathrm{c}\) may also be modulated by advection, as it is inversely related with the urban wind speed. A combination of in situ, remotely sensed, and model simulation data were used to show that the inter-seasonality in \({\Delta } T_\mathrm{s}\), and, to a lesser extent, in \({\Delta } T_\mathrm{c}\), may be related to the change in the land use of the rural site between the pre-monsoon and the monsoon periods. Results suggest that the degree of coupling of \({\Delta } T_\mathrm{s}\) and \({\Delta } T_\mathrm{c}\) may be a strong function of land use and land cover.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Turbulent Flow Over Large Roughness Elements: Effect of Frontal and Plan Solidity on Turbulence Statistics and Structure

Wind-tunnel experiments were carried out on fully-rough boundary layers with large roughness (\(\delta /h \approx 10\), where h is the height of the roughness elements and \(\delta \) is the boundary-layer thickness). Twelve different surface conditions were created by using LEGO? bricks of uniform height. Six cases are tested for a fixed plan solidity (\(\lambda _\mathrm{P}\)) with variations in frontal density (\(\lambda _\mathrm{F}\)), while the other six cases have varying \(\lambda _\mathrm{P}\) for fixed \(\lambda _\mathrm{F}\). Particle image velocimetry and floating-element drag-balance measurements were performed. The current results complement those contained in Placidi and Ganapathisubramani (J Fluid Mech 782:541–566, 2015), extending the previous analysis to the turbulence statistics and spatial structure. Results indicate that mean velocity profiles in defect form agree with Townsend’s similarity hypothesis with varying \(\lambda _\mathrm{F}\), however, the agreement is worse for cases with varying \(\lambda _\mathrm{P}\). The streamwise and wall-normal turbulent stresses, as well as the Reynolds shear stresses, show a lack of similarity across most examined cases. This suggests that the critical height of the roughness for which outer-layer similarity holds depends not only on the height of the roughness, but also on the local wall morphology. A new criterion based on shelter solidity, defined as the sheltered plan area per unit wall-parallel area, which is similar to the ‘effective shelter area’ in Raupach and Shaw (Boundary-Layer Meteorol 22:79–90, 1982), is found to capture the departure of the turbulence statistics from outer-layer similarity. Despite this lack of similarity reported in the turbulence statistics, proper orthogonal decomposition analysis, as well as two-point spatial correlations, show that some form of universal flow structure is present, as all cases exhibit virtually identical proper orthogonal decomposition mode shapes and correlation fields. Finally, reduced models based on proper orthogonal decomposition reveal that the small scales of the turbulence play a significant role in assessing outer-layer similarity.     

Simplification and Validation of a Spectral-Tensor Model for Turbulence Including Atmospheric Stability

A spectral-tensor model of non-neutral, atmospheric-boundary-layer turbulence is evaluated using Eulerian statistics from single-point measurements of the wind speed and temperature at heights up to 100 m, assuming constant vertical gradients of mean wind speed and temperature. The model has been previously described in terms of the dissipation rate \(\epsilon \), the length scale of energy-containing eddies \(\mathcal {L}\), a turbulence anisotropy parameter \(\varGamma \), the Richardson number Ri, and the normalized rate of destruction of temperature variance \(\eta _\theta \equiv \epsilon _\theta /\epsilon \). Here, the latter two parameters are collapsed into a single atmospheric stability parameter z / L using Monin–Obukhov similarity theory, where z is the height above the Earth’s surface, and L is the Obukhov length corresponding to \(\{Ri,\eta _\theta \}\). Model outputs of the one-dimensional velocity spectra, as well as cospectra of the streamwise and/or vertical velocity components, and/or temperature, and cross-spectra for the spatial separation of all three velocity components and temperature, are compared with measurements. As a function of the four model parameters, spectra and cospectra are reproduced quite well, but horizontal temperature fluxes are slightly underestimated in stable conditions. In moderately unstable stratification, our model reproduces spectra only up to a scale \(\sim \) 1 km. The model also overestimates coherences for vertical separations, but is less severe in unstable than in stable cases.