Wind structure near a city centre

Wind speed and temperature profiles to a height of 8 m were recorded for 30-, 60-, and 90-min averaging times over a striated snow surface at the geographic South Pole during the austral winter of 1975. A gradient Richardson number was calculated for each averaging time to determine conditions of neutral stability under which the logarithmic wind law would hold. A log-linear regression technique was used to determine values of aerodynamic roughness height (Z ₀) for those profile averages recorded in conditions of neutral stability. A plot of Z ₀ as a function of average wind direction revealed a variation in Z ₀ of almost three orders of magnitude, from 0.01 to 7 cm, over 120 deg of wind direction. A simple model is presented to justify the fact that aerodynamic roughness is a function of wind direction and erosion history.     

A theory for the scalar roughness and the scalar transfer coefficients over snow and sea ice

Although the bulk aerodynamic transfer coefficients for sensible (C _H ) and latent (C _E ) heat over snow and sea ice surfaces are necessary for accurately modeling the surface energy budget, they have been measured rarely. This paper, therefore, presents a theoretical model that predicts neutral-stability values of C _H and C _E as functions of the wind speed and a surface roughness parameter. The crux of the model is establishing the interfacial sublayer profiles of the scalars, temperature and water vapor, over aerodynamically smooth and rough surfaces on the basis of a surface-renewal model in which turbulent eddies continually scour the surface, transferring scalar contaminants across the interface by molecular diffusion. Matching these interfacial sublayer profiles with the semi-logarithmic inertial sublayer profiles yields the roughness lengths for temperature and water vapor. When coupled with a model for the drag coefficient over snow and sea ice based on actual measurements, these roughness lengths lead to the transfer coefficients. C _E is always a few percent larger than CH. Both decrease monotonically with increasing wind speed for speeds above 1 m s^–1, and both increase at all wind speeds as the surface gets rougher. Both, nevertheless, are almost always between 1.0 × 10^–3 and 1.5 × 10^–3.     

Critical Evaluation of Scalar Roughness Length Parametrizations Over a Melting Valley Glacier

We present a field investigation over a melting valley glacier on the Tibetan Plateau. In the ablation zone, aerodynamic roughness lengths (z _0M ) vary on the order of 10⁻⁴–10⁻² m, whose evolution corresponds to three melt phases with distinct surface cover and moisture exchange: snow (sublimation/evaporation), bare ice (deposition/condensation), and ice hummocks (sublimation/evaporation). Bowen-ratio similarity is validated in the stably stratified katabatic winds, which suggests a useful means for data quality check. A roughness sublayer is regarded as irrelevant to the present ablation season, because selected characteristics of scalar turbulence over smooth snow are quite similar to those over hummocky ice. We evaluate three parametrizations of the scalar roughness lengths (z _0T for temperature and z _0q for humidity), viz. key factors for the accurate estimation of sensible heat and latent heat fluxes using the bulk aerodynamic method. The first approach is based on surface-renewal models and has been widely applied in glaciated areas; the second has never received application over an ice/snow surface, despite its validity in (semi-)arid regions; the third, a derivative of the first, is proposed for use specifically over rough ice defined as z _0M > 10⁻³ m or so. This empirical z _0M threshold value is deemed of general relevance to glaciated areas (e.g. ice sheet/cap and valley/outlet glaciers), above which the first approach gives notably underestimated z _0T,q . The first and the third approaches tend to underestimate and overestimate turbulent heat/moisture exchange, respectively, frequently leading to relative errors higher than 30%. Comparatively, the second approach produces fairly low errors in energy flux estimates both in individual melt phases and over the whole ablation season; it thus emerges as a practically useful choice to parametrize z _0T,q in glaciated areas. Moreover, we find all three candidate parametrizations unable to predict diurnal variations in the excess resistances to humidity transfer, thus encouraging more efforts for improvement.     

Field And Wind-Tunnel Studies Of Aerodynamic Roughness Length

The aerodynamic roughness length (z₀) values of three Gobi desert surfaces were obtained by measurement of the boundary-layer wind profile in the field. To clarify the factors affecting the Gobi surface aerodynamic roughness length, a wind-tunnel experiment was conducted. The wind-tunnel simulation shows that z₀ values increase with increasingsize and coverage of roughness elements. Especially, the shape and height of roughnesselements are more important than other factors in affecting roughness length. The roughness length increases with decreasing values of the geometric parameter (the ratio of element horizontal surface area to height, ) of roughness elements. But at a higher free stream velocity, the height is more important than the shape in affecting roughness length.     

Subfilter-scale Fluxes over a Surface Roughness Transition. Part I: Measured Fluxes and Energy Transfer Rates

A wind-tunnel experiment was designed and carried out to study the effect of a surface roughness transition on subfilter-scale (SFS) physics in a turbulent boundary layer. Specifically, subfilter-scale stresses are evaluated that require parameterizations and are key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. The surface transition considered in this study consists of a sharp change from a rough, wire-mesh covered surface to a smooth surface. The resulting magnitude jump in aerodynamic roughnesses, M = ln(z ₀₁/z ₀₂), where z ₀₁ and z ₀₂ are the upwind and downwind aerodynamic surface roughnesses respectively, is similar to that of past experimental studies in the atmospheric boundary layer. The two-dimensional velocity fields used in this study are measured using particle image velocimetry and are acquired at several positions downwind of the roughness transition as well as over a homogeneous smooth surface. Results show that the SFS stress, resolved strain rate and SFS transfer rate of resolved kinetic energy are dependent on the position within the boundary layer relative to the surface roughness transition. A mismatch is found in the downwind trend of the SFS stress and resolved strain rate with distance from the transition. This difference of behaviour may not be captured by some eddy-viscosity type models that parameterize the SFS stress tensor as proportional to the resolved strain rate tensor. These results can be used as a benchmark to test the ability of existing and new SFS models to capture the spatial variability SFS physics associated with surface roughness heterogeneities.     

Incorporation of a high-roughness lower boundary into a mesoscale model for studies of dry deposition over complex terrain

For flow over natural surfaces, there exists a roughness sublayer within the atmospheric surface layer near the boundary. In this sublayer (typically 50z ₀ deep in unstable conditions), the Monin-Obukhov (M-O) flux profile relations for homogeneous surfaces cannot be applied. We have incorporated a modified form of the M-O stability functions (Garratt, 1978, 1980, 1983) in a mesoscale model to take account of this roughness sublayer and examined the diurnal variation of the boundary-layer wind and temperature profiles with and without these modifications. We have also investigated the effect of the modified M-O functions on the aerodynamic and laminar-sublayer resistances associated with the transfer of trace gases to vegetation. Our results show that when an observation height or the lowest level in a model is within the roughness sublayer, neglect of the flux-profile modifications leads to an underestimate of resistances by 7% at the most.     

The Parameterisation of Scalar Transfer over Rough Ice

We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z _s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z ₀). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z ₀ below a threshold value of approximately 10^?3 m, the model performs well and in accord with earlier studies. However, with growing hummock size, z ₀ increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z ₀ = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z _s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z ₀  >  10^?3 m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z _s /z ₀)  =  1.5  ? 0.2 ln (Re _*)  ? 0.11(ln (Re _*))².     

Aerodynamic Scaling for Estimating the Mean Height of Dense Canopies

We used an aerodynamic method to objectively determine a representative canopy height, using standard meteorological measurements. The canopy height may change if the tree height is used to represent the actual canopy, but little work to date has focused on creating a standard for determining the representative canopy height. Here we propose the ‘aerodynamic canopy height’ h _a as the most effective means of resolving the representative canopy height for all forests. We determined h _a by simple linear regression between zero-plane displacement d and roughness length z ₀, without the need for stand inventory data. The applicability of h _a was confirmed in five different forests, including a forest with a complex canopy structure. Comparison with stand inventory data showed that h _a was almost equivalent to the representative height of trees composing the crown surface if the forest had a simple structure, or to the representative height of taller trees composing the upper canopy in forests with a complex canopy structure. The linear relationship between d and z ₀ was explained by assuming that the logarithmic wind profile above the canopy and the exponential wind profile within the canopy were continuous and smooth at canopy height. This was supported by observations, which showed that h _a was essentially the same as the height defined by the inflection point of the vertical profile of wind speed. The applicability of h _a was also verified using data from several previous studies.     

Flux–Profile Relationships Over a Fetch Limited Beech Forest

The influence of an internal boundary layer and a roughness sublayer on flux–profile relationships for momentum and sensible heat have been investigated for a closed beech forest canopy with limited fetch conditions. The influence was quantified by derivation of local scaling functions for sensible heat flux and momentum (_h and _m) and analysed as a function of atmospheric stability and fetch. For heat, the influences of the roughness sublayer and the internal boundary layer were in agreement with previous studies. For momentum, the strong vertical gradient of the flow just above the canopy top for some wind sectors led to an increase in _m, a feature that has not previously been observed. For a fetch of 500 m over the beech forest during neutral atmospheric conditions, there is no height range at the site where profiles can be expected to be logarithmic with respect to the local surface. The different influence of the roughness sublayer on _h and _m is reflected in the aerodynamic resistance for the site. The aerodynamic resistance for sensible heat is considerably smaller than the corresponding value for momentum.     

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet

To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.     

Derivation of the relationship between the Obukhov stability parameter and the bulk Richardson number for flux-profile studies

Using the relationship between the bulk Richardson numberR _z and the Obukhov stability parameterz/L (L is the Obukhov length), formally obtained from the flux-profile relationships, methods to estimatez/L are discussed. Generally,z/L can not be uniquely solved analytically from flux-profile relationships, and it may be defined using routine observations only by iteration. In this paper, relationships ofz/L in terms ofR _z obtained semianalytically were corrected for variable aerodynamic roughnessz ₀ and for aerodynamic-to-temperature roughness ratiosz ₀/z _T, using the flux-profile iteration procedure. Assuming the so-called log-linear profiles to be valid for the nearneutral and moderately stable region (z/L<1), a simple relationship is obtained. For the extension to strong stability, a simple series expansion, based on utilisation of specified universal functions, is derived.For the unstable region, a simple form based on utilisation of the Businger-Dyer type universal functions, is derived. The formulae yield good estimates for surfaces having an aerodynamic roughness of 10^–5 to 10^–1 m, and an aerodynamic-to-temperature roughness ratio ofz ₀/z _T=0.5 to 7.3. When applied to the universal functions, the formulae yield transfer coefficients and fluxes which are almost identical with those from the iteration procedure.     

On the impact of thermal stability on some rough flow effects over mobile surfaces

Calculations are made of the effects of thermal stability under a range of conditions, over the sea and land, on the physical factors (including the critical wind speed) affecting dust-storm generation, snow drift, and rough sea conditions. The computational procedure involves the surface friction velocity, u _*, and its relation with the aerodynamic roughness over aerodynamically rough, mobile surfaces. The results indicated that even at relatively high wind speeds, thermal effects under extreme advection situations may be significant, particularly for those properties of the agitated surface dependent on u _* ³ and u _* ⁴.     

The roughness of wind waves

In this paper, a new dynamic roughness equation for airflow over wind waves is proposed, based on relationships of dimensionless parameters and the 1986 HEXMAX field data. The equation can be considered as a modification of the Charnock formula. The data are also compared with the parameterizations of Toba and Koga, and of Hsu, and the consequences of the new equation for aerodynamic drag are discussed. In the new equation, the drag coefficient is expressed as a function of both wind speed and wave age c _p/u_*. This testifies to the strong wave-age dependence of the drag, found in several experimental studies. It is also in good agreement with results of a theoretical model of the airflow-seawave interaction, proposed by Janssen (1989).Also: Dept. of Applied Physics, Techn. Univ. Delft, the Netherlands.     

Aerodynamic Parameters of Regular Arrays of Rectangular Blocks with Various Geometries

The aerodynamic effects of various configurations of an urban array were investigated in a wind-tunnel experiment. Three aerodynamic parameters characterising arrays—the drag coefficient (C _d ), roughness length (z _o) and displacement height (d)—are used for analysis. C _d is based on the direct measurement of the total surface shear using a floating element, and the other two parameters are estimated by logarithmic fitting of the measured wind profile and predetermined total drag force. The configurations of 63 arrays used for measurement were designed to estimate the effects of layout, wind direction and the height variability of the blocks on these parameters for various roughness packing densities. The results are summarised as follows: (1) The estimated C _d and z _o of the staggered arrays peak against the plan area index (λ _p ) and frontal area index (λ _f ), in contrast with values for the square arrays, which are less sensitive to λ _p and λ _f . In addition, the square arrays with a wind direction of 45° have a considerably larger C _d , and the wind direction increases z _o/H by up to a factor of 2. (2) The effect of the non-uniformity of roughness height on z _o is more remarkable when λ _f exceeds 20%, and the discrepancy in z _o is particularly remarkable and exceeds 200%. (3) The effect of the layout of tall blocks on C _d is stronger than that of short blocks. These results indicate that the effects of both wind direction and the non-uniformity of the heights of buildings on urban aerodynamic parameters vary greatly with λ _p and λ _f ; hence, these effects should be taken into account by considering the roughness packing density.     

An analytical method for the determination of the roughness parameters over complex regions

The values of roughness length for momentum z ₀and zero-plane displacement d ₀over a hilly rough complex region with vegetation were evaluated without any assumption concerning z ₀and d ₀.It was found that for widely scattered profile data, the method of least squares will not give a reasonable result in determining the roughness parameters. For this purpose, the method of maximum correlation was introduced instead. This method gave a fair result for captive balloon observations conducted in hilly terrain mainly covered with forest in the northwestern part of the Kanto Plain, Japan.     

北京城市下垫面动力学参数的计算与分析

采用北京325 m铁塔2008—2012年的单层超声观测资料,基于莫宁-奥布霍夫相似理论(Monin-Obukhov similarity theory)和前人提出的最小误差分析方法,计算了铁塔周边下垫面的零平面位移高度和动力粗糙度长度。结果表明,由于铁塔位于北京市区,其周边下垫面呈现极其复杂的非均匀性,所以对应铁塔周边不同的扇区,零平面位移高度和动力粗糙度长度各有不同。平均而言,在2008—2012年间,铁塔周边下垫面的零平面位移高度为34.4 m,动力粗糙度长度为1.16 m。此外,综合前人的计算结果发现,铁塔周边的零平面位移高度和动力粗糙度长度在2001年之前呈显著增加的趋势,而在2001年以后并未增长,这一现象与铁塔周边的城市化进程相对应。     

Momentum Absorption in Rough-Wall Boundary Layers with Sparse Roughness Elements in Random and Clustered Distributions

A wind-tunnel experiment has been used to investigate momentum absorption by rough surfaces with sparse random and clustered distributions of roughness elements. An unusual (though longstanding) method was used to measure the boundary-layer depth δ and friction velocity u _* and thence to infer the functional relationship z ₀/h = f(λ) between the normalised roughness length z ₀/ h and the roughness density λ (where z ₀ is the roughness length and h the mean height of the roughness elements). The method for finding u _* is based on fitting the velocity defect in the outer layer to a functional form for the dimensionless velocity-defect profile in a canonical zero-pressure-gradient boundary layer. For the conditions investigated here, involving boundary layers over sparse roughness with strong local heterogeneity, this velocity-defect-law method is found to be more robust than several alternative methods for finding u _* (uw covariance, momentum integral and slope of the logarithmic velocity profile).The experimental results show that, (1) there is general agreement in the relationship z ₀/h = f(λ) between the present experiment with random arrays and other wind-tunnel experiments with regular arrays; (2) the main effect of clustering is to increase the scatter in the z ₀/h = f(λ) relationship, through increased local horizontal heterogeneity; (3) this scatter obscures any trend in the z ₀/h = f(λ) relationship in response to clustering; and (4) the agreement between the body of wind-tunnel data (taken as a whole) and field data is good, though with scatter for which it is likely that a major contribution stems from local horizontal heterogeneity in the field.     

Roughness length for heat over an urban canopy

The roughness length for heat z_T was evaluated over an urban canopy, using the measured sensible heat flux and radiometric temperature. To overcome thermal heterogeneity in the urban area, the measured radiometric temperature was transformed into the equivalent temperature of an upward longwave radiation flux. The equivalent temperature was found to provide an effective parameterization of the radiometric temperature. The daytime average of the resulting ln(z_T/z₀) was 10, where z₀ is the aerodynamic roughness length. This result generally agrees with previous studies; however, the anthropogenic heat is a large uncertainty, which could cause an error at least 240% in z_T.