On the log-linear wind profile and the relationship between shear stress and stability characteristics over the sea

Wind and stability characteristics in the atmospheric surface boundary layer at a height,Z, less than 20 m above the sea were examined in nine oceanic investigations. The analysis lends further support to the utility of the log-linear wind-profile law in the stability region of –0.4Z/L0.9, whereL is the Monin-Obukhov length. However, it is also shown that, inasmuch as better than 90% of the measurements fall within the range of ¦Z/L¦ 0.25, and inasmuch as this correction to the drag coefficient under neutral conditions amounts to less than 10%, the familiar logarithmic wind law may be used rather than the log-linear form. A wind-stress drag coefficient,C _d (=1.2×10^–3 between 1.0 m Z 18.3 m), is thus recommended for general deepwater oceanic applications. The situation over shallow water, which is different, is discussed briefly.     

Simplified expressions for vegetation roughness length and zero-plane displacement as functions of canopy height and area index

Using a previous treatment of drag and drag partition on rough surfaces, simple analytic expressions are derived for the roughness length (z ₀) and zero-plane displacement (d) of vegetated surfaces, as functions of canopy height (h) and area index (). The resulting expressions provide a good fit to numerous field and wind tunnel data, and are suitable for applications such as surface parameterisations in atmospheric models.     

Intermittent wind close to the ground within a grass canopy

Wind speed was measured at a height of 1 cm above the ground and at several other heights in and above a canopy of tall fescue grass (Festuca arundinacea) using single hot-wire and triple hot-film anemometers. The plant area density in the canopy was concentrated close to the ground, with 75% of the plant area standing belowz=15 cm, wherez is height above the ground. The frequency distributions of horizontal wind speeds,s, were sharply skewed towards positive values at all measurement heights, but were most highly skewed near the ground where the coefficient of skewness ranged from 1.6 to 2.9. Above mid-canopy height, the frequency distribution ofs was described reasonably well by a Gumbel extreme value distribution. Average wind speed,S, decreased exponentially with depth into the canopy with an exponential scale length of abouth/2.8, whereh is the height of the canopy. Atz=1 cm, the value ofS was about 11% of the surface-layeru _*. The standard deviation of the fluctuations of the vertical and horizontal components of the wind speed also decreased exponentially with depth inside the canopy with a scale length of abouth/2.5.Inside the canopy, the Eulerian integral time scales for the vertical ( _w ) and horizontal ( _u ) components of wind speed were about 0.1 s and 1.0 s, respectively, and were approximately constant with height. Above the canopy, these time scales increased sharply and, atz=2.25h, _w and _u were approximately 1.0 and 3.0s, respectively. Turbulence length scales in the vertical and downwind directions, _u and _w ·U, respectively, were approximately 1 cm for heights between 1 to 10 cm above the ground inside the canopy, while atz=2.25h, they were about 55 cm and 277 cm. Relatively quiescent periods (lulls) in the air close to the ground were interrupted frequently by gusts. The frequency of occurrence of gusts appears to be correlated with the value of the local shear near the top of the canopy.     

The extent of the unstable Monin-Obukhov layer for temperature and humidity above complex hilly grassland

Potential temperature, specific humidity and wind profiles measured by radiosondes under unstable but windy conditions during FIFE in northeastern Kansas were analyzed within the framework of Monin-Obukhov similarity. Around 86% of these profiles were found to have a height range over which the similarity, formulated in terms of the Businger-Dyer functions, is valid and for which the resulting surface fluxes are in good agreement with independent measurements at ground stations. When scaled with the surface roughness z ₀ = 1.05 m and the displacement height d ₀ = 26.9 m, for the potential temperature this height range was 45 (±31) (z – d ₀ )/z ₀ 104 (±54) and the comparison of the profile-derived surface fluxes with the independent measurements gave a correlation coefficient of r = 0.96. For the specific humidity these values are 42 (±29) (z – d ₀ )/z ₀ 96 (±38) and r = 0.94. In terms of the height of the bottom of the inversion H _i , in the morning hours the upper limit of (z – d ₀ ) in the Monin-Obukhov layer is approximately 0.3H _i , whereas for a fully developed ABL it is closer to 0.1H _i . Probably, as a result of the short sampling times and perhaps also of the small gradients under the windy conditions, the exact height range of validity was difficult to establish from a mere inspection of these profiles.     

Aerodynamic Roughness Length of Fresh Snow

This study presents the results from a series of wind-tunnel experiments designed to investigate the aerodynamic roughness length z ₀ of fresh snow under no-drift conditions. A two-component hot-film anemometer was employed to obtain vertical profiles of velocity statistics in a zero pressure gradient turbulent boundary layer for flow over naturally deposited snow surfaces. The roughness of these snow surfaces was measured by means of digital photography to capture characteristic length scales that can be related to z ₀. Our results show that, under aerodynamically rough conditions, the mean value of the roughness length for fresh snow is \({\langle{z}_{0}\rangle= 0.24}\) mm with a standard deviation σ(z ₀) = 0.05 mm. In this study, we show that variations in z ₀ are associated with variations in the roughness geometry. The roughness measurements suggest that the estimated values of z ₀ are consistent with the presence of irregular roughness structures that develop during snowfalls that mimic ballistic deposition processes.     

On the role of roughness lengths in flux parameterizations of boundary-layer models

Flux-profile relationships based on surface-layer similarity theory are used to derive relationships between the Monin-Obukhov stability parameter = z/L and the bulk Richardson number Ri _b . In contrast to previous studies, the roughness length for heat, z _0h ,is assumed unequal to the roughness length for momentum, z _0m .For the stable case, an analytic expression of in terms of Ri _b can be derived and in the unstable case, the solution is obtained through a simple iterative process.Errors introduced from the simplification of z _0h = z _0m are evaluated and are shown to be very significant in most cases. Thus, this error in many practical applications may invalidate the intended solution.     

A model of airflow above changes in surface heat flux,temperature and roughness for neutral and unstable conditions

A numerical model of airflow in the lowest 50–100 m of the atmosphere above changes in surface roughness and temperature or heat flux has been developed based on boundary layer approximations, the Businger-Dyer hypotheses for the non-dimensional wind shear and heat flux and a mixing length hypothesis.Results have been obtained for several situations, in particular, airflow with neutral upstream conditions encountering a step change in surface temperature or heat flux with no roughness change. In these cases large increases in shear stress at the outer edge of the internal boundary layer are predicted. The case of unstable upstream flow encountering a step change to zero heat flux is also considered.Two situations that may be encountered near the shores of the Great Lakes are considered.Notation B Businger-Dyer constant (= 16.0) in form for _M, _H - c _p Specific heat at constant pressure - g Acceleration due to gravity - H Upward vertical heat flux - H ₀ , H ₁ Surface heat fluxes for x < 0, x 0 - k von Kármán's constant ( = 0.4) - l Mixing length - L Monin-Obukhov length - L ₀ Upstream value of L - m Ratio of roughness lengths (= z ₁/z ₀) - RL ^* Non-dimensional parameter, see Equations (20, 22 and 24) - RL ₁ ^* Same as RL ^* but with z ₁ scaling (= mRL ^*) - T Scaled temperature - T ₀ (z) Upstream temperature profile - u ₀, u ₁(x) Surface friction velocities for x < 0, x 0 - U, W Horizontal and vertical mean velocities - U ₀ (z) Upstream velocity profile - x, z Horizontal and vertical coordinates - z _i Local roughness length     

Flux-gradient profiles for momentum and heat over an urban surface

Summary In this paper, we evaluate the applicability of flux-gradient relationships for momentum and heat for urban boundary layers within the Monin-Obukhov similarity (MOS) theory framework. Although the theory is widely used for smooth wall boundary layers, it is not known how well the theory works for urban layers. To address this problem, we measured the vertical profiles of wind velocity, air temperature, and fluxes of heat and momentum over a residential area and compared the results to theory. The measurements were done above an urban canopy whose mean height z_h is 7.3 m. 3-D sonic anemometers and fine wire thermocouples were installed at 4 heights in the region 1.5z_h < z < 4z_h. We found the following: (1) The non-dimensional horizontal wind speed has good agreement with the stratified logarithmic profile predicted using the semi-empirical Monin-Obukov similarity (MOS) function, when it was scaled by the surface friction velocity that is derived from the shear stress extrapolated to the roof-top level. (2) The scaled gradient of horizontal wind speed followed a conventional semi-empirical function for a flat surface at a level (z/z_h = 2.9), whereas, in the vicinity of the canopy height was larger than the commonly-used empirical relationship. (3) The potential temperature profile above the canopy shows dependency on the atmospheric stability and the scaled gradient of temperature is in good agreement with a conventional shear function for heat. In the case of heat, the dependency on height was not found. (4) The flux-gradient relationship for momentum and heat in the region 1.5z_h < z < 4z_h was rather similar to that for flat surfaces than that for vegetated canopies.     

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.     

Heat and mass transfer to and from surfaces with dense vegetation or similar permeable roughness

A differential equation is obtained to describe the concentration of passive admixtures (water vapor, sensible heat, pollutants, CO₂, etc.) of turbulent flow inside a dense and uniform vegetational canopy. The profiles of eddy diffusivity, wind speed and shear stress are assumed to be exponential decay functions of depth below the top of the canopy. This equation is solved for the case of a vegetation with constant concentration of the admixture at the foliage surfaces. The solution is used to formulate bulk mass or heat transfer coefficients, which can be applied to practical problems involving surfaces covered with a vegetation or with similar porous or fibrous roughness elements. The results are shown to be consistent with experimental data presented by Chamberlain (1966), Garratt and Hicks (1973) and Garratt (1978). Calculations with the model illustrate that, as compared to its behavior over surfaces with bluff roughness elements, ln(z ₀/z _0c ) (wherez ₀ is the momentum roughness andz _0c , the scalar roughness) for permeable roughness elements is relatively insensitive tou _* and practically independent ofz ₀.     

Heat flux parameterization for sparse and dense grasslands with the analytical land-atmosphere radiometer model (alarm)

The dependence of radiometric surface temperature (_s) on view angle and the unclear definition of the aerodynamic temperature, which is the temperature that gives the correct sensible heat flux estimate at a specified roughness length, bring about a challenge in estimating sensible heat flux from _s. An analytical-land-atmosphere-radiometer model (ALARM) has been developed to convert _s taken at any zenith view angle to a clearly defined equivalent isothermal surface temperature, _i, at a defined scalar roughness length. ALARM is an analytical model based on K-theory that links the foliage temperature profile to the radiometric surface temperature and the temperature felt by the turbulent lower atmosphere. ALARM has previously been applied with slightly different values of its parameters to several grassland sites of varying canopy density. Our objective in this study was to apply ALARM to these and to one additional dataset with a single parameterization. When compared to the reference (measured) values of sensible heat flux H, ALARM estimates of H had root mean square errors of about 35 W m^-2. These results were comparable to those from two other simple canopy models also tested with these datasets.     

Momentum and scalar transfer coefficients over aerodynamically smooth antarctic surfaces

Vertical profiles of wind speed, temperature and humidity were used to estimate the roughness lengths for momentum (z ₀), heat (z _H ) and moisture (z _Q) over smooth ice and snow surfaces. The profile-measurements were performed in the vicinity of a blue ice field in Queen Maud Land, East Antarctica. The values ofz ₀ over ice (3·10^–6 m) seem to be the smallest ever obtained over permanent, natural surfaces. The settling of snow on the ice and the loss of momentum at saltating snow particles serve as momentum dissipating processes during snow-drift events, expressed as a strong dependence ofz ₀ on u_#.The scalar roughness lengths and surface temperature can be evaluated from the temperature and humidity profile measurements if the ratioz _H /z _Q is specified. This new method circumvents the difficult measurement of surface temperature. The scalar roughness lengths seem to be approximately equal toz0 for a large range of low roughness Reynolds numbers, despite the frequent occurrence of drifting snow. Possible reasons for this agreement with theory of non-saltating flow are discussed.     

Equations for the Drag Force and Aerodynamic Roughness Length of Urban Areas with Random Building Heights

We use a conceptual model to investigate how randomly varying building heights within a city affect the atmospheric drag forces and the aerodynamic roughness length of the city. The model is based on the assumptions regarding wake spreading and mutual sheltering effects proposed by Raupach (Boundary-Layer Meteorol 60:375?C395, 1992). It is applied both to canopies having uniform building heights and to those having the same building density and mean height, but with variability about the mean. For each simulated urban area, a correction is determined, due to height variability, to the shear stress predicted for the uniform building height case. It is found that u _*/u _*R , where u _* is the friction velocity and u _*R is the friction velocity from the uniform building height case, is expressed well as an algebraic function of ?? and ?? _h /h _m , where ?? is the frontal area index, ?? _h is the standard deviation of the building height, and h _m is the mean building height. The simulations also resulted in a simple algebraic relation for z ₀/z _0R as a function of ?? and ?? _h /h _m , where z ₀ is the aerodynamic roughness length and z _0R is z ₀ found from the original Raupach formulation for a uniform canopy. Model results are in keeping with those of several previous studies.     

Estimation of the sensible heat flux of a semi-arid area using surface radiative temperature measurements

In the framework of an international field program for the study of semi-arid areas, observations were done in the region called La Crau in southern France. In this paper, the use of the surface radiative temperature for the determination of the sensible heat flux is addressed. We found that, once proper values of the roughness length of momentum (z ₀) and heat (z ₀h) are set, the sensible heat flux can be reliably predicted with a one-layer resistance model using standard observations of wind speed and air temperature, together with the surface temperature. The latter quantity has to be known with a precision better than ±2°C. From our observations, the value of the parameterB ^–1k ^–1 In (z ₀ z ₀h) was found to be 9.2, which falls between values quoted by Brutsaert (1982) for grass and bluff bodies.     

A Field Study of the Mean Pressure About a Windbreak

To provide additional field data for assessingwindbreak flow models, mean ground-level pressurehas been measured upstream and downstream from along porous fence (height H = 1.25 m, resistancecoefficient k _r = 2.4). Measurements were madeduring periods of near-neutral stability and near-normallyincident flow, with the fence standing on bare soil(roughness length, z ₀ 0.8 cm;H/z ₀ 160), or within a plant canopy. The mean pressure field,measured far from the ends of the fence, was foundto be quite insensitive to mean wind direction( , zero for perpendicular flow), for| | less than about 25°.In the absence of a canopy, during each measurementperiod the minimum pressure occurred at the closestsampling location to leeward of the windbreak, thepressure-gradient in most cases beingmaximally-adverse in the immediate lee, and decayingwith increasing downwind distance (x). On one day ofmeasurements, however, the pressure gradient over2 x/H 6 (H = windbreak height) resembled theleeward plateau identified by Wang and Taklein their numerical studies. Perhaps thisoccasional feature was only due to instrumenterror. Nevertheless a plateau of sorts wasindicated in similar measurements by Judd andPrendergast (with H = 1.92 m, z ₀ 1.2 cm;H/z ₀ 160, k _r 3). Therefore,existence of a leeward pressure plateau behind athin fence cannot be definitely ruled out.When the windbreak was placed in a canopy, minimumsurface pressure was displaced downwind. Thisagrees with the wind-tunnel study of Judd, Raupach and Finnigan,and is consistent with a simple simulation reported here.     

Determination of Roughness Lengths for Heat and Momentum Over Boreal Forests

The roughness length for momentum (z_0m), zero-plane displacementheight (d), and roughness length for heat (z_0h) are importantparameters used to estimate land-atmosphere energy exchange. Although many different approaches have been developed to parameterizemomentum and heat transfer, existing parameterizations generally utilizehighly simplified representations of vegetation structure. Further, a mismatch exists between the treatments used for momentum and heat exchange and those used for radiative energy exchanges. In this paper, parameterizations are developed to estimate z_0m, d, and z_0h for forested regimes using information related to tree crown density and structure. The parameterizations provide realistic representationfor the vertical distribution of foliage within canopies, and include explicit treatment for the effects of the canopy roughness sublayer and leaf drag on momentum exchange. The proposed parameterizationsare able to realistically account for site-to-site differences in roughness lengths that arise from canopy structural properties.Comparisons between model predictions and field measurements show good agreement, suggesting that the proposed parameterizations capture the most important factors influencing turbulent exchange of momentumand heat over forests.     

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.     

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.     

On the determination of zero-plane displacement and roughness length for flow over forest canopies

This paper discusses the importance of the aerodynamic characteristics of forest and other similar canopies to modelling of boundary-layer flow and to estimating the diffusivity coefficients of turbulence transfer mechanisms over such canopies.The hypothesis of Marunich (1971) reported by Tajchman (1981) that the zero-plane displacement, d, equals the upward displacement of the flow trajectory, is critically examined. It is concluded that Marunich's hypothesis is conceptually incorrect and that calculations of d based on Marunich's hypothesis are inherently in error.This paper presents a method based on the mass conservation principle and uses wind profiles in and above a forest canopy as the sole input for determining d, z ₀ and u _*.Sensitivities of calculated results to measurements errors of wind profile data are evaluated. It is found that an error of less than 1% in wind in the logarithmic regime above the canopy can introduce up to 100% errors in calculated values of d, z ₀ and u _*. It is also found that the high sensitivity to wind data accuracy, characteristic of the present method, can be used as a guide for the selection of high quality canopy wind data.     

Calculation of velocity skewness in real and artificial plant canopies

Wind velocities within a plant canopy are much more strongly skewed than those of the air flow above. We have examined the governing Eulerian equations for the velocity products u _i, u _j u_k using data from a wind tunnel study with an artificial canopy consisting of an array of 5 cm lengths of monofilament fishing line, and from measurements in corn (Zea mays L).Simple parameterizations for pressure-velocity correlations, and for the quadruple velocity products allowed reasonably accurate calculations of the third moments using measured profiles of the mean velocity, variance and covariance fields. Comparisons of individual terms in the rate equations for ovu _i, u _j u _krevealed that diffusion (from above) and mean shear were most important in creating large skewness in the canopy. A drag term also contributed but was of lesser importance. These terms were balanced by return-to-isotropy and a turbulence interaction term. A quasi-Gaussian approximation considerably underestimated the magnitude of the fourth moments within the canopy.