Vertical cross-spectra of horizontal velocity components at the Boulder observatory

Cross-spectra between horizontal wind components at different levels of the Boulder Atmospheric Observatory (BAO) tower lead to the following conclusions:

1. Davenport's hypothesis is satisfied that coherence decays exponentially with the ratio of vertical separation to horizontal wave length, at least to very small values of coherence.
2. The decay coefficients increase with z/L for z/L < 0.5. For larger stabilities, irregular fluctuations with periods of order 10–20 min have considerable vertical coherence. Results at BAO are quite consistent with those elsewhere.
3. Eddy slopes in vertical planes increase with wind shear up to a point where the slope (horizontal delay over vertical separation) is just above 2. Beyond that point, the systematic increase of slopes with shear ceases. Since wind shear decreases upward, slopes tend to decrease upward. Slopes for lateral components are significantly larger than those for u-components.

     

Wind characteristics at the Boulder Atmospheric Observatory

Wind speeds at the 300 m tower at the Boulder Atmospheric Observatory have been analyzed. This tower is located in slightly rolling farmland. The following conclusions have been reached:

1. For west winds, the terrain is sufficiently uniform for simple surface-layer theory to be adequate without modification even though the air has moved up a small slope to reach the tower. For south and southeast winds, ‘effective’ roughness lengths must be introduced, which are significantly larger than the ‘true’ roughness length.
2. Useful wind estimates up to 150 m can be made from winds at 10 m and stability information, provided the ‘effective’ roughness length is known.
3. The observations are consistent with a von Kármán constant of 0.35.

     

Atmospheric boundary layer research at Cabauw

At Cabauw, The Netherlands, a 213 m high mast specifically built for meteorological research has been operational since 1973. Its site, construction, instrumentation and observation programs are reviewed. Regarding analysis of the boundary layer at Cabauw, the following subjects are discussed:

- terrain roughness;

- Monin-Obukhov theory in practice;

- the structure of stable boundary layers;

- observed evolution of fog layers;

- inversion rise and early morning entrainment;

- use of the geostrophic wind as a predictor for wind profiles;

- height variation of wind climate statistics;

- air pollution applications: long range transport and short range dispersion;

- dependence of sound wave propagation on boundary-layer structure;

- testing of weather and climate models.

     

The internal boundary layer — A review

A review is given of relevant work on the internal boundary layer (IBL) associated with:

1. Small-scale flow in neutral conditions across an abrupt change in surface roughness,
2. Small-scale flow in non-neutral conditions across an abrupt change in surface roughness, temperature or heat/moisture flux,
3. Mesoscale flow, with emphasis on flow across the coastline for both convective and stably stratified conditions.

The major theme in all cases is on the downstream, modified profile form (wind and temperature), and on the growth relations for IBL depth.     

Multiple windbreaks: An aeolean ensemble

Near-neutral measurements of the turbulent wind field within and above a sequence of 15 parallel windbreaks on a flat pastoral site are presented. The windbreak fences each had a porosity of 60% and were equally-spaced at 6 times their height (h = 2 m). The following conclusions seem justified for wind directions within 10 ° of the normal to the array:

1. Above the windbreaks (2h), mean windspeeds first decreased and then increased asymptotically to a value in equilibrium with the new surface roughness. At 0.5h, windspeeds exhibited a slow increase down the entire array.
2. Reflecting differences in approach flows, the drag on the initial fence was almost twice that on barriers farther downstream. This reduction in momentum extraction per windbreak was associated with an elevation in the zero-plane displacement to a level equal to 0.8h.
3. At positions well-removed from the initial fences, mean windspeeds were reduced throughout the entire region below shelter height. In this region, the flow became increasingly dominated by downward moving air with velocities much greater than the local average. The zone of reduced turbulence was small, extending only 2h downstream of a barrier at a height of 0.25h. This corresponded with the region excluded from smoke trails released at the top of windbreaks.
4. An approximate TKE budget mid-way between windbreaks 7 and 8 suggests that shear and wake production peak near z = h and that production is balanced by dissipation and vertical transport components. Advective and inertial interaction terms are negligible at this midway position but are likely to be major sources of TKE closer to the windbreak. Local equilibrium is attained above z = 1.5h implying the existence of a constant-stress layer.

The measurements show the practical difficulty of simultaneously reducing both mean windspeeds and turbulence levels with repeated windbreaks at conventional spacings for horticultural applications.     

Limitations of an eddy-correlation technique for the determination of the carbon dioxide and sensible heat fluxes

A modified infrared CO₂ gas analyzer, a small thermocouple assembly, a heated-thermocouple anemometer for horizontal wind, and a propeller-type vertical wind sensor were used to measure the eddy fluxes of heat and CO₂ above a corn crop. Experimental results of these fluxes are discussed. The main sources of errors of the eddy fluxes using these instruments were estimated:

1. Sensors with a time constant of 0.5 s appear to be fast enough to detect most of the vertical CO₂ transfer as long as the sensors are located at least one meter above the crop surface.
2. The deviation from steady-state conditions for 10-min periods was found to have a significant effect on the eddy flux estimates.
3. Temperature fluctuations of the air sample passing through the CO₂ infrared gas analyzer were found to be non-negligible but could be easily corrected.
4. A 1° misalignment of the vertical anemometer affected these eddy fluxes by less than 10% under all circumstances studied.

     

Atmospheric turbulence within and above a Douglas-fir stand. Part I: Statistical properties of the velocity field

This is the first of two papers reporting the results of a study of the turbulence regimes and exchange processes within and above an extensive Douglas-fir stand. The experiment was conducted on Vancouver Island during a two-week rainless period in July and August 1990. The experimental site was located on a 5^o slope. The stand, which was planted in 1962, and thinned and pruned uniformly in 1988, had a (projected) leaf area index of 5.4 and a heighth=16.7 m. Two eddy correlation units were operated in the daytime to measure the fluctuations in the three velocity components, air temperature and water vapour density, with one mounted permanently at a height of 23.0m (z/h=1.38) and the other at various heights in the stand with two to three 8-hour periods of measurement at each level. Humidity and radiation regimes both above and beneath the overstory and profiles of wind speed and air temperature were also measured. The most important findings are:

1. A marked secondary maximum in the wind speed profile occurred in the middle of the trunk space (aroundz/h=0.12). The turbulence intensities for the longitudinal and lateral velocity components increased with decreasing height, but the intensity for the vertical velocity component had a maximum atz/h=0.60 (middle of the canopy layer). Magnitudes of the higher order moments (skewness and kurtosis) for the three velocity components were higher in the canopy layer than in the trunk space and above the stand.
2. There was a 20% reduction in Reynolds stress fromz/h=1.00 to 1.38. Negative Reynolds stress or upward momentum flux perisistently occurred atz/h=0.12 and 0.42 (base of the canopy), and was correlated with negative wind speed gradients at the two heights. The longitudinal pressure gradient due to the land-sea/upslope-downslope circulations was believed to be the main factor responsible for the negative Reynolds stress.
3. Momentum transfer was highly intermittent. Sweep and ejection events dominated the transfer atz/h=0.60, 1.00 and 1.38, with sweeps playing the more important role of the two atz/h=0.60 and 1.00 and the less important role atz/h=1.38. But interaction events were of greater magnitude than sweep and ejection events atz/h=0.12 and 0.42.

     

Bare-ground surface heat and water exchanges under dry conditions: Observations and parameterization

A simplified land-surface parameterization is tested against bare-soil data collected during the EFEDA experiment conducted in Spain in June 1991. A complete data set, made up of soil properties as well as hydrological and atmospheric measurements, is described and discussed. The 11-day data set is characterized by very dry conditions and high surface temperatures during the day. Large values of sensible and soil heat fluxes and small values of surface evaporation (≈1 mm/day) were observed. This data set was modelled, leading to the following conclusions:

1. In the model, the parameterization provides values of the soil thermal properties and subsequently of the predicted soil heat fluxes which are overestimated when compared with the observations.
2. Following the literature, a value of the ratio between the roughness lengths for momentumZ _oand heatZ _ohof close to 10 for fairly homogeneous areas of bare soil and vegetation is used. This value leads to a fair prediction of the surface temperature. If the roughness lengths were taken to be equal, as is often assumed in atmospheric modelling, a poorer prediction results.
3. Finally, the vapor phase transfer mode is found dominant close to the surface and a modified parameterization including this effect is proposed. It allows a fair prediction of both surface evaporation and near-surface water content.

     

Characteristics of the sea breeze in Attica,Greece

The characteristics of the sea breeze in the Attica region of Greece, in which Athens is located, have been studied for occasions of weak synoptic-scale pressure gradient. The analysis is based on synoptic observations from six meteorological stations, three on the coast and three inland. The three inland stations and one of the coastal stations lie almost in a straight line at different distances from the coast. For each meteorological station, the basic characteristics of the sea breeze were determined, i.e.,

1. The mean number of sea-breeze days for each calendar month.
2. The monthly mean wind speed for each synoptic hour.
3. The times of onset and cessation of the sea breeze.
4. The monthly vector mean wind, and its constancy ‘Constancy’ is defined as 100{itV{inr}/V{ins}}, where {itV{inr}} is the magnitude of the vector mean wind, and {itV{ins}} is the scalar mean wind speed. See Brooks and Carruthers (1953). (In this paper, the factor 100 is not used.) for each synoptic hour.
5. For days on which there was a sea breeze at Helliniko (the coastal reference station), the percentage number of days on which there was also a sea breeze at the given station.

An attempt was also made to determine further characteristics, such as the inland penetration of the sea breeze, its depth, the spatial and temporal variation of wind speed and direction, and the existence of the return flow. Finally, the properties of the land breeze are briefly outlined.     

Numerical study on the bulk heat transfer coefficient for a variety of vegetation types and densities

The relationship between the geometrical structure of a canopy layer and the bulk transfer coefficient was investigated using a numerical canopy model. The following results were obtained:

1. The bulk transfer coefficients for momentum and heat, C _M and C _H , change with non-dimensional canopy density C _* each has a maximum.
2. The value of C _M is always larger than the value of C _H for a canopy with c _m > c _h , c _m and c _h being the drag coefficient and the heat transfer coefficient of an individual canopy element, respectively.
3. The value of C _* at which C _H has its maximum value is larger than the value of C _* at which C _M has its maximum. Therefore, the reciprocal of the sublayer Stanton number b _h ^?1 ranges between 50 and 65 for C _* around 0.1 while it ranges between 0 and 30 for C _* < 10^?2 and C _* > 2 (when c _m = 0.5).
4. The value of B _H ^?1 in the present study is consistent with most available observations, except for canopies of medium density (when C _* is around 0.1) for which no observational value has been obtained.

     

Problems of atmospheric shear flows and their laboratory simulation

Shear flows generated by movement of the atmosphere near the earth's surface are accompanied by complexities not ordinarily encountered in the treatment of turbulent boundary layers. Problems arising from the following physical features are considered:

1. thermal stratification;
2. surface roughness in the form of forests and cities;
3. non-uniformity of surface roughness and/or temperature (leading to 3-dimensional turbulent boundary layers);
4. surface irregularities in the form of hilly and mountainous topography.

The complex nature of atmospheric shear flows has stimulated efforts to study their characteristics in the laboratory under controlled conditions. Accordingly, questions of similarity between the laboratory and the atmospheric flows for both mean and turbulent quantities arise. Similarity criteria, or appropriate scaling relationships, are discussed. Wind tunnels designed for investigations related to atmospheric shear flows are described. These facilities are shown to have a capability for simulating such flows for a wide range of the physical features listed above.     

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.     

Bulk characteristics of heat transfer in the unstable,baroclinic atmospheric boundary layer

Field data for the unstable, baroclinic, atmospheric boundary layer over land and over the sea are considered in the context of a general similarity theory of vertical heat transfer. The dependence of δθ/θ_* upon logarithmic functions of h _c z _T and stability (through the similarity function C) is clearly demonstrated in the data. The combined data support the conventional formulation for the heat transfer coefficient δθ/θ_* when,

1. the surface scaling length is z _T (« z ₀), the height at which the surface temperature over land is obtained by extrapolation of the temperature profile
2. the height scale is taken as the depth of convective mixing h _c
3. the temperature profile equivalent of the von Karman constant is taken as 0.41
4. areal average, rather than single point, values of δθ are employed in strongly baroclinic conditions. No significant effect of baroclinity or the height scale ratio as proposed in the general theory is found. Variations in C about a linear regression relation against stability are most probably due to uncertainties in the areal surface temperature and to experimental errors in general temperature measurements.

     

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

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

Lake Kinneret (Sea of Galilee) and its exceptional wind system

From 1973–1976, research was performed around the Sea of Galilee, aimed at examining the wind regime in the area and whether the area develops a land-sea breeze despite its particular topographical location.

The main conclusions were:
1. During the summer mornings a lake breeze develops, blowing towards the shores of the lake. It ceases at the peak of its development when a westerly wind, originating in the development of a breeze along the Israeli Mediterranean coast, plunges towards the lake.
2. Late at night, a wind flow develops from the land towards the lake, which combines with the katabatic winds that blow along the steep slopes surrounding the Kinneret.
3. The stations at the upper level, at a height of 400–500 m above the Kinneret, are not affected by the lake breeze during the day or by the land breeze at night.
4. In winter, the Kinneret lake breeze is almost as developed as in summer, because the westerly winds, originating in the Mediterranean sea breeze which hardly develops in this season, do not plunge into the Kinneret.

     

The scavenging of nitrate by clouds and precipitation

A model with spectral microphysics was developed to describe the scavenging of nitrate aerosol particles and HNO₃ gas. This model was incorporated into the dynamic framework of an entraining air parcel model with which we computed the uptake of nitrate by cloud drops whose size distribution changes with time because of condensation, collision-coalescence and break-up. Significant differences were found between the scavenging behavior of nitrate and our former results on the scavenging behavior of sulfate. These reflect the following chemical and microphysical differences between the two systems:

1. nitrate particles occur in a larger size range than sulfate particles.
2. HNO₃ has a much greater solubility than SO₂ and is taken up irreversibly inside the drops in contrast to SO₂.
3. nitric acid in the cloud water is formed directly on uptake of HNO₃ gas whereas on uptake of SO₂ sulfuric acid is formed only after the reaction with oxidizing agents such as e.g., H₂O₂ or O₃.
4. nitrate resulting from uptake of HNO₃ is confined mainly to small drops, whereas sulfate resulting from uptake of SO₂ is most concentrated in the largest, oldest drops, which have had the greatest time for reaction.

Sensitivity studies showed that the nitrate concentration of small drops is significantly affected by the mass accommodation coefficient.     

Surface Heterogeneity Effects on Regional-Scale Fluxes in the Stable Boundary Layer: Aerodynamic Roughness Length Transitions

The effects of abrupt streamwise transitions of the aerodynamic roughness length ( $z_\mathrm{o}$ z o ) on the stable atmospheric boundary layer are evaluated using a series of large-eddy simulations based on the first Global Energy and Water Cycle Experiment Atmospheric Boundary Layer intercomparison study (GABLS1). Four $z_\mathrm{o}$ z o values spanning three orders of magnitude are used to create all possible binary distributions with each arranged into patches of characteristic length scales equal to roughly one-half, one, and two times the equivalent homogeneous boundary-layer height. The impact of the heterogeneity on mean profiles of wind speed and temperature, on surface fluxes of heat and momentum, and on internal boundary-layer dynamics are considered. It is found that $z_\mathrm{o}$ z o transitions do not significantly alter the functional relationship between the average surface fluxes and the mean profiles of wind speed and potential temperature. Although this suggests that bulk similarity theory is applicable for modelling the stable boundary layer over $z_\mathrm{o}$ z o heterogeneity, effective surface parameters must still be specified. Existing models that solve for effective roughness lengths of momentum and heat are evaluated and compared to values derived from the simulation data. The existing models are unable to accurately reproduce both the values of the effective aerodynamic roughness lengths and their trends as functions of patch length scale and stability. A new model for the effective aerodynamic roughness length is developed to exploit the benefits of the other models tested. It accurately accounts for the effects of the heterogeneity and stratification on the blending height and effective aerodynamic roughness length. The new model provides improved average surface fluxes when used with bulk similarity.     

Surface influence upon vertical profiles in the nocturnal boundary layer

Near-surface wind profiles in the nocturnal boundary layer, depth h, above relatively flat, tree-covered terrain are described in the context of the analysis of Garratt (1980) for the unstable atmospheric boundary layer. The observations at two sites imply a surface-based transition layer, of depth z _*, within which the observed non-dimensional profiles Φ _M ⁰ are a modified form of the inertial sub-layer relation \(\Phi _M \left( {{z \mathord{\left/ {\vphantom {z L}} \right. \kern-0em} L}} \right) = \left( {{{1 + 5_Z } \mathord{\left/ {\vphantom {{1 + 5_Z } L}} \right. \kern-0em} L}} \right)\) according to $$\Phi _M^{\text{0}} \simeq \left( {{{1 + 5z} \mathord{\left/ {\vphantom {{1 + 5z} L}} \right. \kern-\nulldelimiterspace} L}} \right)\exp \left[ { - 0.7\left( {{{1 - z} \mathord{\left/ {\vphantom {{1 - z} z}} \right. \kern-\nulldelimiterspace} z}_ * } \right)} \right]$$ , where z is height above the zero-plane displacement and L is the Monin-Obukhov length. At both sites the depth z _* is significantly smaller than the appropriate neutral value (z _*N ) found from the previous analysis, as might be expected in the presence of a buoyant sink for turbulent kinetic energy.     

Aircraft measurements of oxides of nitrogen along the eastern rim of the Asian continent: Winter observations

Aircraft observations of oxides of nitrogen (NO _y ), measured with a ferrous sulfate converter, over the sea surrounding the Japanese islands (30–43° N, 131–141° E) were carried out in the winter of 1983 and 1984 at altitudes mostly between 3 and 8 km. NO _y defined here is the sum of NO, NO₂, and other unstable oxides of nitrogen that are converted to NO by ferrous sulfate. The main observations were:

1. Over the Pacific Ocean between the latitudes of 30–35° N, the observed NO _y mixing ratio between 3 and 8 km was a fairly constant 200 pptv. The NO mixing ratio increased with altitude from 15 pptv at 3 km to 35 pptv at 7 km.
2. Over the Sea of Japan, tropospheric NO _y mesured between 1 and 6 km started increasing with latitude North of 35° N and reached about 1000 pptv at 40° N.
3. NO _y was measured in an air mass transported from the stratosphere near a tropopause fold region. When the ozone mixing ratio was between 80 and 140 ppbv, the NO _y mixing ratio was about 200 pptv.