Displaced-Beam Small Aperture Scintillometer Test. Part Ii: Cases-99 Stable Boundary-Layer Experiment

The performance of the Scintec displaced-beam small aperture scintillometer (DBSAS) in the stable boundary layer (SBL) is investigated using data gathered during the CASES-99 experiment in Kansas, U.S.A. The DBSAS is superior to the eddy-covariance method in determining vertical fluxes of sensible heat and momentumclose to the ground and/or over short (< 1 min) averaging intervals. Both aspects are of importance in the shallow and non-stationary SBL.The friction velocity, u_*, the temperature scale, _*, and from these the sensible heat flux, H, were calculated from the indirectly determined dissipation rate, , and the structure parameter of temperature, C_T ², by the DBSAS, which was operated over a path length of 112 m. All these variables are compared with eddy-covariance data for 10-minute time averages. Previously reported systematic errors in the DBSAS, overestimation of u_* for low u_* values and underestimation of u_* for high u_* values, have in part been dealt with by adjusting the beam displacement distance from 2.7 mm to 2.6 mm in the calculations. The latter adjustment is presented as a working hypothesis, not a general solution.     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.     

Turbulence characteristics in a near neutrally stratified urban atmosphere

Turbulence measurements from the city of Uppsala, Sweden, are analysed. Measurements were taken at two sites: one in the central area, ca. 6 m above roof level, the average building height being ca. 15 m; the other at ca. 8 and 50 m above the ground on a tower situated 100 m downwind of a sharp discontinuity between the densely built-up urban area and flat grass-covered land. The average stability was close to neutral, the range being -0.2 < z/L < 0.2. The main emphasis of the study is on the non-dimensional standard deviations of the velocity components _i /u _*t and on the corresponding non-dimensional energy spectra, u _*t being a local velocity scale defined as _i /( _l is the local momentum flux). Comparison with results obtained from surface-layer measurements at ideal sites (with u _*, being the ordinary friction velocity) shows good general agreement. The most complete agreement is found for the tower 50 m measurements, a result which is notable as this measurement point is found to be within a distinctly transitional zone between the urban and post-urban boundary layers. The results from the central city measurement point are also fairly close to the ideal results, the deviations found being small in view of the fact that the site is probably inside the layer in which the roughness elements (the buildings) have direct influence. The measurements at the tower 8 m level show certain distinct deviations from ideal results: all three _i /u _*l , are higher by ca. 10%, the excessive energy being found at the low frequency end of the spectrum. Arguments are presented for this feature to be due to a spectral lag effect.     

Standard deviation of vertical two-point longitudinal velocity differences in the atmospheric boundary layer

The standard deviation of vertical two-point longitudinal velocity fluctuation differences is analyzed experimentally with eleven sets of turbulence measurements obtained at the NASA 150-m ground-winds tower site at Cape Kennedy, Florida. It is concluded that /u _*0 is proportional to (fz/u _*0)^0.22, where the coefficient of proportionality is a function of fz/u _*0 and u _*0/fL ₀. The quantities f and L₀ denote the Coriolis parameter and the surface Monin-Obukhov stability length, respectively; u _*0 is the surface friction velocity; z is the vertical distance between the two points over which the velocity difference is calculated; and zz is the mean height of the mid-point of the interval z above natural grade. The results of the analysis are valid for 20<-u _*0/fL ₀<2000.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

Similarity forms for ground-source surface-layer diffusion

Past work on analyzing ground-source diffusion data in terms of surface-layer similarity theory is reviewed; these analyses assume that _z /L orh/L is a function of u ^* x/L (where h = Q/ dy). It is argued that an alternative scaling, h ^*/L versus x/L, is nearly as universal in that it is very weakly influenced by surface roughness, except for a modest influence in the free convective case (h ^* = Q/u ^* dy); the advantage of this scaling is that it eliminates the need to reassess as vertical diffusion progresses. The Prairie Grass data set is adjusted for the difference in source and sampling heights, and is plotted with this scaling. Simple analytic equations are suggested that fit the resultant data plots for stable and unstable conditions, and suggestions are made towards practical application of these results.On assignment from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce.     

Observation of lake breeze penetration and subsequent development of the thermal internal boundary layer for the Nanticoke II Shoreline Diffusion Experiment

The turbulent structure of the lake breeze penetration and subsequent development of the thermal internal boundary layer (TIBL) was observed using a kytoon-mounted ultrasonic anemometer-thermometer. The lake breeze penetrated with an upward rolling motion associated with the upward flow near the lake breeze front. After the lake breeze front passed, the behaviors of the velocity and temperature at the top of the lake breeze layer were similar to those found in convective boundary layers (CBL). Comparing _gq/_*, _u /w _* and _w /w _* between the present observation of TIBL development after the passage of the lake breeze front and CBL data from the literature, the /_* values showed reasonable agreement; however, _u /w _* and _w /W_* had smaller values in the TIBL than in the CBL at higher altitudes. This is due to the differences in the mean velocity profiles. While the CBL has a uniform velocity profile, the TIBL has a peak at lower elevation due to the lake breeze penetration; the velocity then decreases with height.Present address: The Institute of Behavioral Science, 1-35-7 Yoyogi, Tokyo 151, Japan.     

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.     

Boundary-layer parameterizations for applied dispersion modeling over urban areas

Many applied dispersion models require the knowledge of boundary-layer parameters such as sensible heat flux,Q _H , friction velocity,u _*, and turbulent energy components, _w and _v . Formulas are suggested for calculating these parameters over a wide variety of types of ground surfaces, based on simple observations of wind speed near the ground and fractional cloud cover, and specification of constants such as roughness length, albedo, and soil moisture availability. Observations ofu _*,Q _H , _w , and _v during field experiments in St. Louis and Indianapolis are used to test the formulas for urban sites. Relative errors of about ±20% in the predictions are seen to occur whenu _*,Q _H , _w , and _v are large. However, when these quantities are small (e.g.,u _* < 0.2 m/s), the errors in the predictions are as large as the mean value of the quantity itself.In addition, it is concluded from studies of available field data and theories that the magnitude of _w is not well-known at elevations above about 100m during the late afternoon and night. Some simple parameterizations for _w . are suggested that are consistent with the observed steady decrease in ground-level concentration in the afternoon and the sudden increase in concentration that can occur a few hours after sunset due to wind shears associated with a low-level jet, for continuous plumes emitted from moderate to tall stacks.     

On characteristics of wind direction fluctuations in the surface layer

It is shown that the ratio of standard deviation of lateral velocity to the friction velocity, /u _*, and therefore wind direction fluctuations, are sensitive to mesoscale terrain properties. Under neutral conditions, /u _* is almost 40% larger in rolling terrain than over a horizontal surface. In the lee of a low mountain, the fluctuations may be 2.5 times as strong as over horizontal terrain. In contrast, vertical velocity fluctuations are little influenced by mesoscale terrain features.Now with Air Weather Service, Offutt AFB, Omaha, Nebraska.     

A note on the analogy between momentum transfer across a rough solid surface and the air-sea interface

The aerodynamic classification of the resistance laws above solid surfaces is based on the use of a so-called Reynolds roughness number Re _s =h _s u _*/, whereh _s is the effective roughness height, -viscosity,u _*-friction velocity. The recent experimental studies reported by Toba and Ebuchi (1991), demonstrated that the observed variability of the sea roughness cannot be explained only on the basis of the classification of aerodynamic conditions of the sea surface proposed by Kitaigorodskii and Volkov (1965) and Kitaigorodskii (1968) even though the latter approach gains some support from recent experimental studies (see for example Geernaertet al. 1986). In this paper, an attempt is made to explain some of the recently observed features of the variability of surface roughness (Toba and Ebuchi, 1991; Donelanet al., 1993). The fluctuating regime of the sea surface roughness is also described. It is shown that the contribution from the dissipation subrange to the variability of the sea surface can be very important and by itself can explain Charnock's (1955) regime.     

The use of σ T as a temperature scale in the atmospheric surface layer

The standard deviation of temperature _T is proposed as a temperature scale and as a velocity scale to describe the behaviour of turbulent flows in the Atmospheric Surface Layer (ASL), instead of _* andu _* of the Monin—Obukhov similarity theory, and ofT _f andU _f used for free convection stability conditions. On the basis of experimental evidence reported in the literature, it is shown that _T T _f andv _* U _f in the free convection region, and _{T *} andv _* U _* in nearneutral and stable conditions. This implies that the proposed scales can be applied for all stabilities. Furthermore, a new length scale is proposed and its relation with Obukhov length is given. Also, a simple semi-empirical expression is presented with which _T andv _* can be evaluated in a rather simple way. Some examples of practical applications are given, e.g., a stability classification for unstable conditions.     

A laboratory study of the mechanism of the oxygen airglow

Further laboratory studies of emission by O(¹ S) and by O₂ A ³ _u ⁺ ,A³ _u andc ¹ _u ^– in the oxygen afterglow lead to the conclusion that Barth's mechanism for the excitation of the auroral green line O ₂ ^* +O(³P=O₂+O(¹S)–(1) is correct and that levelsv=6 and 7 of O₂ A ³ _u ⁺ are Barth precursors. The value ofk ₁=7×10^–11 cm³ s^–1 deduced for these levels is shown to be in fair agreement with atmospheric measurements.     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

Variability of surface roughness and turbulence intensities at a coastal site in India

Surface-layer features with different prevailing wind directions for two distinct seasons (Southwest Monsoon and Northeast Monsoon) on the west coast of India are studied using data obtained from tower-based sensors at a site located about 500 m from the coast. Only daytime runs have been used for the present analysis. The surface boundary-layer fluxes have been estimated using the eddy correlation method. The surface roughnessz ₀ obtained using the stability-corrected wind profiles (Paulson, 1970) has been found to be low for the Southwest monsson season. For the other season,z ₀ is relatively high. The drag coefficientC _D varies with height in the NE monsoon season but not in the season with lowz ₀. This aspect is reflected in the wind profiles for the two seasons and is discussed in detail. The scaling behaviour of friction velocityu _* and the turbulence intensity of longitudinal, lateral and vertical winds _u, _v and _w, respectively) are further examined to study their dependence on fetch. Our study shows that for the non-dimensional case, _u/u_* and _v/u_* do not show any surface roughness dependence in either season. On the other hand, for _w/u_* for the season with lowz ₀, the values are seen to agree well with that of Panofskyet al. (1977) for homogeneous terrain whereas for the other season with highz ₀, the results seem to conform more to the values observed by Smedman and Högström (1983) for coastal terrain. The results are discussed in the light of observations by other investigators.     

A Richardson Number criterion for the air-sea interface

It is shown that the observationally determined roughness relation z ₀ = u _* ²/g in which g is the acceleration of gravity, u _*, is the friction velocity in air, and = 0.0185 (Wu, 1982) for the wind profile over the sea surface relative to the surface current, is consistent with the existence of a Richardson Number criterion at the air-sea interface in which the critical Richardson Number, Ri_c = 1, such that all the shear energy is converted into potential energy.     

Convective Profile Constants Revisited

This paper examines the interpolation betweenBusinger–Dyer (Kansas-type) formulae,_u = (1 -1 6 )^-1/4 and_t = (1 - 16 )^-1/2, and free convection forms. Based on matching constraints, the constants, a_u and a_t, in the convective flux-gradient relations, _u = (1 - a_u )^-1/3 and _t = (1 - a_t )^-1/3, are determined. It isshown that a_u and a_t cannot be completely independent if convective forms are blended with theKansas formulae. In other words, these relationships already carryinformation about a_u and a_t. This follows because the Kansas relations cover a wide stability range (up to = - 2), which includes a lower part of the convective sublayer (about 0.1 < - < 2). Thus, there is a subrange where both Kansas and convective formulae are valid. Matching Kansas formulae and free convection relations within thesubrange 0.1 < - < 2 and independently smoothing ofthe blending function are used to determine a_u and a_t. The values a_u = 10 for velocity and a_t = 34for scalars (temperature and humidity) give a good fit. This new approacheliminates the need for additional independent model constants and yields a`smooth' blending between Kansas and free-convection profileforms in the COARE bulk algorithm.     

On the determination of the height of the Ekman boundary layer

The heighth of the Ekman turbulent boundary layer determined by the momentum flux profile is estimated with the aid of considerations of similarity and an analysis of the dynamic equations. Asymptotic formulae have been obtained showing that, with increasing instability,h increases as ¦¦^1/2 (where is the non-dimensional stratification parameter); with increasing stability, on the other hand,h decreases as ^–1/2. For comparison, a simple estimate of the boundary-layer heighth _u determined by the velocity profile is given. As is shown, in unstable stratification,h _u behaves asymptotically as ¦¦^–1, i.e., in a manner entirely different from that ofh .     

Radiation balance in a tropical city — Delhi (India)

The impact of urbanization is assessed by comparing values of the radiation parameters at an urban location with those of a rural site. Urban Delhi was divided according to land-use and the effects of urbanization was studied on incoming short-wave (K), albedo, incoming longwave (L, outgoing longwave (L), and net radiation (Q ^*), were individually studied at four representative sites (Rural, Commercial, Residential and Industrial). MaximumK was observed in the rural and commercial areas whereas highL was observed in the commercial and industrial locations. High depletion ofK of the order of 13% was observed for the industrial location in the winter season. An increase ofL in, the industrial location is of the same order as that of the commercial location, i.e. 20%. The residential location shows quite moderate (4.6%) depletion ofK in comparison with other sites.Nomenclature U Urban - R Rural - K Incoming Short-wave radiation - L Incoming Long-wave radiation - L Outgoing Long-wave radiation - Q ^* Net Radiation - Albedo - K ^* Net Short-wave radiation - L ^* Net Long-wave radiation