Observed variations of σθ and σΦ in the nocturnal drainage flow in a deep valley

The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms^-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and     

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.     

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.     

Markov-chain simulation of particle dispersion in inhomogeneous flows: The mean drift velocity induced by a gradient in Eulerian velocity variance

The Langevin equation is used to derive the Markov equation for the vertical velocity of a fluid particle moving in turbulent flow. It is shown that if the Eulerian velocity variance _wE is not constant with height, there is an associated vertical pressure gradient which appears as a force-like term in the Markov equation. The correct form of the Markov equation is: w(t + t) = aw(t) + b _wE + (1 – a)T _L ( _wE ²)/z, where w(t) is the vertical velocity at time t, a random number from a Gaussian distribution with zero mean and unit variance, T _L the Lagrangian integral time scale for vertical velocity, a = exp(–t/T _L), and b = (1 – a ²)^1/2. This equation can be used for inhomogeneous turbulence in which the mean wind speed, _wE and T _L vary with height. A two-dimensional numerical simulation shows that when this equation is used, an initially uniform distribution of tracer remains uniform.     

On the memory of wind standard deviation for upstream roughness

Big eddies in the outer part of the atmospheric boundary layer contribute to the variance of the horizontal velocity fluctuations near the surface. Because of the slow adjustment of these eddies to new boundary conditions, they carry the roughness characteristics of a large upstream terrain. A scaling relation is proposed that accounts for the memory effects in the big eddies. It is concluded that the standard deviation of the horizontal wind ( _u ) measured at a given height is representative for the shear stress at greater height. This gives at least qualitative support to existing work where u is used for exposure correction of mean wind.     

Seasonal variation of Σθ with wind speed,direction and stability

For an airport site near Visakhapatnam, India, and based on 10 years of data for the months of January, April, August and October, values of are given as a function of wind speed, wind direction and Pasquill diffusion category.     

Estimates of diabatic wind speed profiles from near-surface weather observations

In this paper we analyse diabatic wind profiles observed at the 213 m meteorological tower at Cabauw, the Netherlands. It is shown that the wind speed profiles agree with the well-known similarity functions of the atmospheric surface layer, when we substitute an effective roughness length. For very unstable conditions, the agreement is good up to at least 200 m or z/L–7(z is height, L is Obukhov length scale). For stable conditions, the agreement is good up to z/L1. For stronger stability, a semi-empirical extension is given of the log-linear profile, which gives acceptable estimates up to ~ 100 m. A scheme is used for the derivation of the Obukhov length scale from single wind speed, total cloud cover and air temperature. With the latter scheme and the similarity functions, wind speed profiles can be estimated from near-surface weather data only. The results for wind speed depend on height and stability. Up to 80 m, the rms difference with observations is on average 1.1 m s^–1. At 200 m, 0.8 m s^–1 for very unstable conditions increasing to 2.1 m s^–1 for very stable conditions. The proposed methods simulate the diurnal variation of the 80 m wind speed very well. Also the simulated frequency distribution of the 80 m wind speed agrees well with the observed one. It is concluded that the proposed methods are applicable up to at least 100 m in generally level terrain.     

Characteristics of the Lyman-alpha humidiometer and measurements of turbulent humidity fluctuations over the ocean

The characteristics of a Lyman-alpha humidiometer have been carefully examined in an air-conditioned test chamber. The results confirm that when carefully used, this humidiometer is suitable for measurements of turbulent humidity fluctuations. Measurements with a Lyman-alpha humidiometer were carried out in the surface boundary layer over the ocean. The relation between turbulent intensity ( _a = a ^ov2) and the friction humidity (a _*) can be expressed as _a = l.6a _*. The spectrum of turbulent humidity for wind speeds larger than 3 m s ^–1 conforms to the similarity law in the surface boundary layer. The spectrum has two characteristic normalized frequencies, namely, a higher peak and a secondary peak (or a shoulder).     

Time series analyses of concentration and wind fluctuations

Analyses of concentration fluctuation (C) spectra from boundary-layer smoke plume experiments at six separate locations show that the spectra from these experiments generally exhibit an inertial subrange at high frequencies with a slope of -5/3 and indicate peak energy at a time period of about 50 to 100 s. These periods of peak energy are a factor of two to five less than those for the peak of the wind speed fluctuation (u or v) spectra. A general spectral formula fits normalized spectra from the U.S. and Australia, where the frequency, n, is made dimensionless by multiplying by the plume dispersion parameter, _y , and dividing by the wind speed, u. Peak energy occurs at a dimensionless frequency of n _y/u equal to about 0.15. The Kolmogorov constant in the inertial subrange is estimated from a set of averaged spectra. Cross-spectra indicate little relation between concentration and wind fluctuations. However, most of the correlation that exists is due to periods larger than about 10 or 20 s.     

A note on plume geometry from a tall stack

Measurements of the cross-wind standard deviation () of an instantaneous plume from a 380 m stack have been made under convective conditions in the boundary layer. Comparison of these measurements with those calculated from the buoyancy-induced diffusion formulation of Weil (1979) shows good agreement.     

A random-walk model for dispersion of heavy particles in turbulent air flow

A random-walk model is presented for calculating the dispersion of heavy particles in a turbulent air flow when only air turbulence statistics and the drag characteristics of the particle are known. Algebraic expressions for the modification of air velocity variance ² and Lagrangian autocorrelation tune-scale T _L,due to particle inertia effects, are derived. These expressions introduce only a very small computational overhead on the random-walk models for inertia-less particles of Wilson et al. (1983). Measurements of T _Land by Snyder and Lumley (1971) for four different particles are used to determine constants in the heavy-particle model. It is shown that the agreement between the model, for a single set of constants, and the dispersion measurements is good for the 47 m hollow glass, 87 m glass, and 47 m copper particles. The predictions for the 87 m corn pollen particles show less satisfactory agreement by underestimating dispersion measurements by 15% after 0.4s. Finally, some aspects of the model's application to spray dispersion in and above a crop canopy are considered.     

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.     

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.     

A comparison of a Eulerian and a Lagrangian time scale for over-water atmospheric flows during stable conditions

Lagrangian integral time scales were calculated from crosswind concentration distributions of oil-fog smoke released from a continuous point source over the ocean during stable atmospheric conditions assuming an exponential correlation function. Variance of the lateral velocity fluctuations, _v ², and the energy dissipation rate, , were obtained from simultaneous Eulerian measurements at the beach. An Eulerian energy dissipation scale defined as _v ²/ was then computed. The ratio of the Lagrangian integral scale to the Eulerian energy dissipation scale was found to be close to 1. This ratio was also estimated to be 1 based on physical and dimensional considerations regarding the cascade of energy. Length scales for longitudinal, lateral and vertical directions were interpreted with a model based on similarity considerations applicable for over-water atmospheric flows.Authored under contract EY-76-C-02-0016 with the U.S. Department of Energy. Accordingly, the U.S. Government retains a nonexclusive, royalty-free license to publish or reproduce the published form of this contribution, or allow others to do so, for U.S. Government purposes.     

Free convection scaling over the ocean derived from buoy measurements during gate

The structure of atmospheric turbulence in the surface layer over the open ocean is examined under conditions of local free convection. The raw data consist of profile and fluctuation measurements of wind and temperature as obtained from a meteorological buoy. For near neutral conditions and for waves running approximately along the wind direction, wave-induced wind fluctuations can be described by a simplified linear theory based on Miles (1957). In this case, the spectrum of wind velocity is given as the sum of two parts; for the turbulent part, the parameterization as obtained by Kaimal et al. (1972) applies, while the wave-induced part is parameterized using a simplification of Miles' linear theory. For cases of local free convection, the measurements of the vertical component of the wind velocity are well described by similarity theory; as expected, _w /(-uw)^1/2 is proportional to (- z/L)^1/3. In order to scale the longitudinal wind velocity component, it seems to be reasonable to extend the list of relevant parameters by the height of the mixed layer z _i. We obtain _u /(- uw)^1/2 (z/z _i)^1/3(- z/L)^1/3 with only a poor correlation coefficient of r = 0.6. Overall, the results of local free convection scaling obtained from direct measurements show good agreement with those obtained from profile measurements. A comparison between direct and indirect determination of turbulent fluxes of momentum shows an unexplained difference of about 20%. This discrepancy is mainly due to a gap in the uw-cospectrum at the swell frequency.     

Turbulence measurements during mistral winds with a 1-dimensional sonic anemometer

Fluctuations in the vertical wind velocity and air temperature were measured with a 1-dimensional sonic anemometer and fine thermocouple over a flat agricultural site in the Rhone Valley, France. Strong Mistral winds with speeds up to 20 m s^–1 kept atmospheric conditions very close to neutral and ensured stationarity. Friction velocities estimated both by eddy correlation (sonic plus Gill Bivane) and inertialdissipation (sonic only) methods agreed within 1 and 5 % respectively of traditional profile measurements over the measured range of 0.2 to 1.2 m s^–1. The coefficient of eddy transport for heat exceeded that of momentum by a factor of 1.38 (± 0.05), a result almost identical to that obtained in the Kansas experiment (Businger et al., 1971). For - 0.15 >= z/L >= 0.05, the ratio _w /u _* was 1.69 and 1.34 for unstable and stable conditions, respectively. For ¦z/L¦ >= 0.05, the ratio /T _* was 1.40 independent of whether neutrality was approached from either stable or unstable conditions.     

Similarity scaling applied to sodar observations of the convective boundary layer above an irregular hill

Nine profiles of the temperature structure parameter C _T ² and the standard deviation of vertical velocity fluctuations ( _w) in the convective boundary layer (CBL) were obtained with a monostatic Doppler sodar during the second intensive field campaign of the First ISLSCP Field Experiment in 1987. The results were analyzed by using local similarity theory. Local similarity curves depend on four parameters: the height of the mixed layer (z _i ), the depth of the interfacial layer (), and the temperature fluxes at the top of the mixed layer (Q _i ) and the surface (Q _o). Values of these parameters were inferred from sodar data by using the similarity curve for C _T ² and observations at three points in its profile. The effects of entrainment processes on the profiles of C _T ² and _wnear the top of the CBL appeared to be described well by local similarity theory. Inferred estimates of surface temperature flux, however, were underestimated in comparison to fluxes measured by eddy correlation. The measured values of _wappeared to be slightly smaller than estimates based on available parmeterizations. These discrepancies might have been caused by experimental error or, more likely, by the distortion of turbulence structure above the site by flow over the nonuniform terrain at the observation site.     

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.     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

Calculation of particle trajectories in the presence of a gradient in turbulent-velocity variance

Recent papers by Wilson et al. (1981b) and Legg and Raupach (1982) give methods for the calculation of particle trajectories in turbulence with a gradient in vertical velocity variance ² _w. However the two methods seem contradictory.This paper demonstrates that in systems in which l(d _w /dz) (where / is the length scale) varies only slowly with height z, the two methods give similar predictions, and indicates why this is the case. For a particular system in which the restriction on l(d _w /dz) is not satisfied, it is shown that neither method is correct but that a simple modification of the method of Wilson et al. (1981b) gives reasonable predictions.