Lateral coherence of longitudinal wind components in strong winds

A combination of lateral coherence measurements of wind speed at five locations suggests that the decay constant is a monotonically increasing function of the ratio of separation to height, under neutral conditions.     

Tropical atmospheric stability and wind structures as inferred from the Indian MST radar observations

Summary The MST (Mesosphere-Stratosphere-Tropospher) Radar Facility at Gadanki (13.47° N, 79.18° E), near Tirupati, Andhra Pradesh, India has been operated over seven diurnal cycles—three in November 1994, one in September 1995 and three in January–February 1996 with an objective to study the wind and stability characteristics in the troposphere and lower-stratosphere. The radar-measured height profiles of both zonal (EW) and meridional (NS) wind components and near-simultaneous radiosonde measurements from Madras (13.04° N, 80.7° E) and Bangalore (12.85° N, 77.58° E), the two stations close to either side of the radar site, have been compared and they are found to be in gross agreement within the limitations of the sensing techniques.The results of the study also indicated multiple stable and turbulent structures/stratification throughout the height region from about 4 to 30 km. It is noticed that the stable layers are well marked around the altitudes 4 km, 12 km and the tropopause while the turbulent layers exist a few kilometers below the tropopause. These stable and turbulent layer structures showed good correspondence with the radar-measured wind gradients and also with the radiosonde-derived temperature and wind distributions over Madras. The maximum positive gradient in the signal-tonoise ratio (SNR) which corresponds to radar tropopause is found to coincide with the greater potential temperature gradient and smaller wind gradient. The time evolution of atmospheric stability structure, derived from the SNR, spectral width and vertical wind revealed a diffused tropopause or tropopause weakening which is found to be associated with broader spectral width and larger gradients of winds. This feature is considered to be due either to the instability associated with large vertical gradients in horizontal winds (dynamical instability) or to the instability generated by the convection (convective instability).With 6 Figures     

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.     

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.     

Wind and temperature profile prediction for diabatic surface layers including strong inversion cases

The micrometeorological research program in Antarctica has provided extensive data on wind and temperature profile structure under strong to extreme inversion conditions (Dalrymple et al., 1966; Lettau et al., 1977). The basic similarity hypotheses and limiting conditions for prediction of diabatic surface layer profiles are summarized. The model by Businger et al. (1971) for dimensionless shear and temperature gradients is revised to conform with the new results for strong stability. A novel similarity hypothesis is introduced to complete the step from shear and gradient prediction to prediction of absolute wind speed, wind energy, and temperature on the basis of prescribed external factors of surface layer structure. The physics of interactions between predicted profile tilting and curving are discussed and used to explain several micrometeorological paradoxes, including that of the elevated minimum of air temperature observed occasionally near the active surface when the energy budget is of the nocturnal type.     

Spectra,variances and length scales in a marine stable boundary layer dominated by a low level jet

A flow situation over coastal waters of the Baltic Sea is studied. The boundary layer was characterized by stable stratification and the presence of a pronounced low level jet at very low height, 30–150 m, above the surface of the sea. The atmospheric surface layer was apparently extremely shallow; thus the non-dimensional wind gradients and temperature gradients derived from measurements at 8 m do not show adherence to Monin-Obukhov similarity, in sharp contrast to findings from the same site at similar stability conditions but with no low level jet. Instead these quantities are shown to be governed by scales characteristic of stable shear flow away from the surface. The height to the jet centre appears to be an important quantity. Thus, for the cases with the lowest jet height values (30–50 m), some turbulent characteristics of the flow (non-dimensional velocity standard deviations and the correlation between the longitudinal and vertical velocity) have values similar to those found for the zero pressure-gradient laboratory boundary layer over a flat plate (the so called canonical boundary layer) rather than the typical values found in atmospheric boundary-layer flow. It was inferred that the large scale fluctuations known as inactive turbulence, as well as gravity waves, were suppressed in this case.     

One-dimensional checkerboards and blending heights

We analyze the checkerboard problem of many alternating surfaces with different properties, on scales up to (say) 3,000 m. Power-law representations of the vertical profiles of mean wind speed and eddy diffusivity lead to solutions in terms of Kelvin and trigonometric functions.These solutions are used to determine blending heights (_*), where deviations from the mean of concentration, or of vertical flux density, fall to some small fraction, , of their value at the surface. Values of _*are important for regional and larger-scale meteorological models. In smaller scale micrometeorological studies, they may serve also as the top levels of surface boundary layers.An important result for both theoretical and experimental contexts is that deviations of flux persist with elevation much more strongly than those of concentration, so that, in general, _* should be based on flux rather than concentration. Representative values of _*, for = 0.05, are of order 5 and 30 m for surface pattern wavelengths of 10² and 10³ m, respectively. Values of _* are robust to changes in adopted power-law indices, and are independent of wind speed. Surface roughness has a mild but calculable effect.     

Determination of similarity functions of the resistance laws for the planetary boundary layer using surface-layer similarity functions

Functional forms of the universal similarity functions A, B (for wind components parallel and normal to the surface stress), and C (for potential temperature difference) are determined based on the generalized theory of the resistance laws for the Planetary Boundary Layer (PBL). The similarity-profile functions for the surface layer are matched with the velocity and temperature-defect profiles that are assumed to have shapes modified by certain powers of nondimensional height z/h, where h is the PBL height. The powers of the outer-layer profile functions are determined, so that the functions become negligible in the surface layer. To close the temperature defect law, an assumption that the temperature gradient across the top of the PBL is continuous with the stratification of the overlying atmosphere is used. The result of this assumption is that nondimensional momentum and temperature profiles in the PBL can be described in terms of four basic ratios: (1) roughness ratio = /h (2) scale-height ratio =|f|h/u_*, (3) ambient stratification parameter =h/_*, and (4) stability parameter =h/L, where L is the Monin-Obukhov length, z₀ is the surface roughness, is the upper-air stratification, u _* is the friction velocity, and _* is the temperature scale at the surface. For stable conditions, the scale-height ratio can be related to the atmospheric stability and the upperair stratification, and the generalized similarity and Rossby number similarity theories become identical. Under appropriate boundary conditions, function A is explicitly dependent on the stability parameter , while B is a function of scale-height ratio , which in turn depends on the stability. Function C is shown to be dependent on the stability and the upper-air stratification, due to the closure assumption used for the temperature profile.The suggested functional forms are compared with other empirical approximations by several authors. The general framework used to determine the functional forms needs to be tested against good boundary-layer measurements.     

On the streamline configuration associated with miles' theory of sea wave growth

The Miles theory of sea wave growth implies that, in a coordinate system moving with the wave phase velocity, closed streamline cats'-eyes are formed around the critical height where the mean wind speed equals the wave phase speed. From inspection of the equations of motion (the Rayleigh system), it is shown that these cats'-eyes may lie either over the wave crests or over the wave troughs, depending upon the behaviour of the mean wind profile. The validity of using Miles first approximation of the streamfunction' to approximate the actual streamline pattern is also discussed.     

Final results of the CONDORS convective diffusion experiment

TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF ₆ values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, _y and _z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations _vand _w 0.6w^* (the convective scale velocity), as measured atz>100m. Surface release _y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger _v and reduced wind speed measured near the surface. Decreased _y at small dimensionless average times is also noted. Finally, convectively scaled dy, C _y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC _y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.     

On the turbulence structure in the stable boundary layer over the Greenland ice sheet

Turbulence measurements from a 30 m tower in the stably stratifiedboundary layer over the Greenland ice sheet are analyzed. The observationsinclude profile and eddy-correlation measurements at various levels. Atfirst, the analysis of the turbulence data from the lowest level (2 m aboveground) shows that the linear form of the non-dimensional wind profile(m) is in good agreement with the observations for z/L <0.4, whereL represents the Obukhov length. A linear regression yieldsm=1+5.8z/L. The non-dimensional temperature profile (h) at the2m level shows no tendency to increase with increasing stability. The datafrom the upper levels of the tower are analyzed in terms of both localscaling and surface-layer scaling. The m and the h values show atendency to level off at large stability (z/>0.4) where represents the local Obukhov length. Hence, the linear form of the functions is no longer appropriate under such conditions. The bestcorrespondence to the data can be achieved when using the expression ofBeljaars and Holtslag for m and h. The vertical profiles of theturbulent fluxes, the wind velocity variances and temperature variance arealso determined. The momentum flux profile and the profiles of wind speedvariances are in general agreement with other observations if a welldeveloped low-level wind maximum occurs, and the height of this maximum isused as a height scale.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

Analytical investigations of horizontal meso-scale momentum equations with Newtonian nudging

Summary This study presents an analytical investigation of the local behaviour of the solution to a mesoscale model with Newtonian nudging when observed winds are time varying. The analysis examines each Fourier component of the time series of observed winds. Unlike the case with a constant observed wind, the nudged wind vector does not asymptotically approach the observed wind. In response to sinusoidal oscillation of the observed wind, the nudged wind vector is always on a half circle connecting the vector ends of the observed and un-nudged modelled winds. When nudging parameter 0, the nudged wind vector approaches the un-nudged wind; when , the nudged wind vector approaches the observed wind. For commonly used values of nudging parameter , the modelled wind field always carries errors.A target nudging scheme is devised in this study in order to ensure the model result is identical to observed winds with sinusoidal oscillation. Investigation shows that such a target wind exists for a finite value of , and the magnitude of the target-nudging term is about the same as that of a normal nudging term if f, wheref is the Coriolis parameter and is the frequency of the wind oscillation.With 7 Figures     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

A radiation and energy balance study of mature forest and clear-cut sites

Components of the radiation and energy balances were measured over a clear-cut area and a mature, mixed forest during the summer of 1981 at the Petawawa National Forestry Institute, Chalk River, Ontario. The work concentrated on the clear-cut site which supported a canopy layer composed primarily of bracken fern and logging remnants.Forty days of radiation data were collected at the clear-cut site. After the first four weeks of measurements (the green season), most of the ferns quickly died, and their foliage changed appearance from a green to brownish colour (the brown season). The daily mean reflection coefficient of solar radiation determined over the green season was 0.20 and decreased to 0.13 for the brown season. The corresponding value for the forest was 0.13, based on a limited amount of data. The clear-cut site received 11% and 21% less net radiation than the forest on a 24-hr and daylight-hours basis, respectively, as a consequence of the higher reflection coefficient and larger daytime longwave radiation emission.A reversing temperature difference measurement system (RTDMS), incorporating ten-junction thermopiles was employed at each site in order to determine Bowen ratios () via differential psychrometry. Both systems performed well, especially the RTDMS over the forest which was capable of resolving very small differences of temperature, typically less than 0.2 °C over a height of 3 m. The mean hourly Bowen ratio, calculated from values from 0800 to 1600 hr, varied from 0.2 to 1.0 for the forest and from 0.4 to 0.8 for the clear-cut site in the green season.A significant canopy heat storage component of the energy balance, Q _S , was found at the clear-cut site. In the early morning, a portion of the available energy was used to heat the biomass materials and air within the canopy layer. The stored heat within the canopy was released later in the day, increasing the available energy total.The daily mean value of the Priestley-Taylor coefficient (Priestley and Taylor, 1972) for the green season at the clear-cut site was 1.14, and individual values tended to increase during wet surface conditions and decrease when the surface dried. The daylight mean value during dry canopy conditions at the forest was 1.05, and much higher values occurred when the canopy was wet. The enhancement of for the wet forest was a result of the evaporation of intercepted rain (which is not limited by stomatal resistance) and the concomitant transfer of sensible heat to the forest.     

A SIMPLE FORMULA FOR ATTENUATION OF EDDY FLUXES MEASURED WITH FIRST-ORDER-RESPONSE SCALAR SENSORS

A simple formula, (1 + (2fmc))-1,is proposed to estimate the attenuation of a scalar flux measurement made by eddy-correlation using a fast-response anemometer and a linear, first-order-response scalar sensor with a characteristic time constant c.In this formula, =7/8 for neutral and unstable stratification within the surface-flux layer and =1 both within the convective boundary layer (CBL) and for stable stratification in the surface layer.fm is the frequency of the peak of the logarithmic cospectrum and can be estimated from fm = nm /z, where z is the measurement height and is thewind speed at that height. The dimensionless frequency at the cospectral maximum nm is estimated here from observations of its behavioras a function of atmospheric stability, z/L within the surface layeror z/zi within the CBL, where L is the Obukhov stability length and zi is the depth of the CBL. The predicted dependence of flux attenuation on measurement height is discussed.     

Occurrence of roll circulations in a shallow boundary layer

Meteorological measurements taken at the Näsudden wind turbine site during slightly unstable conditions have been analyzed. The height of the convective boundary layer (CBL) was rather low, varying between 60 and 300 m. Turbulence statistics near the ground followed Monin-Obukhov similarity, whereas the remaining part of the boundary layer can be regarded as a near neutral upper layer. In 55% of the runs, horizontal roll vortices were found. Those were the most unstable runs, with -z _i/L > 5. Spectra and co-spectra are used to identify the structures. Three roll indicators were identified: (i) a low frequency peak in the spectrum of the lateral component at low level; (ii) a corresponding increase in the vertical component at mid-CBL; (iii) a positive covariance {ovvw} together with positive wind shear in the lateral direction (V/z) in the CBL. By applying these indicators, it is possible to show that horizontal roll circulations are likely to be a common phenomenon over the Baltic during late summer and early winter.     

An improved calibration for tethered balloon wind measurements

A previously published technique for using tethered spherical balloons as anemometers for measuring light low-level winds has been further developed. Earlier data on the relationship between the aerodynamic drag coefficient and the Reynolds number of spherical rubber balloons were combined with a large number of new data and re-analysed; and the errors in the relationship were estimated. The results allowed a more accurate calculation of wind speed from the deflection of a tethered balloon from the vertical. When combined with a new technique for calculating the effects of the tether, this enabled light to moderate low-level winds at fixed heights up to 600 m or more to be measured with simple, cheap, and readily mobile equipment; and a slight modification of the technique allowed measurement of winds in and above fog. Wind speeds measured by the ballon technique showed reasonably good agreement with measurements by an anemometer carried beneath the balloon.Glossary of Symbols a, b, c Coefficients in the relationship between lnC _d and lnR - A Quantity under square root in solution for lnV whena0 - C _d Wind drag coefficient for balloon - C _dc Value ofC _d given by calibration curve of Table I - D Dynamic wind pressure force on balloon - F Buoyant free lift of balloon with load - Re Reynold's number of balloon (sphere) - R = Re/10⁵ - r Radius of sphere - T Tension in tether - V Wind speed - ₈₃() =(lnC _dc -lnC _d ) when 83° , or 0 for other - Error in lnC _d - Elevation of tether where attached to balloon - Elevation of balloon from ground tether point - Molecular viscosity of air - Ratio of circumference to diameter of circle - Density of air     

Application of the Priestley-Taylor evaporation model to assess the influence of soil moisture on the evaporation from a large weighing lysimeter and class a pan

The influence of soil moisture on evaporation from a 6-m grass-covered lysimeter and from Class A pans was assessed for one summer using the -parameter of the Priestley-Taylor evaporation model appropriate for the individual surfaces computed on a daily basis. Net radiation over the pan was estimated from above-grass measurements using a correlation established between the two, using measurements made in the previous two summers. Changes in heat storage of the water were considered in the derivation of for the pan. A unique relationship for the particular conditions of the site was determined between the for the lysimeter and soil moisture, approaching 1.29 at soil moisture near field capacity, but decreasing to as low as 0.5 for dry soil. The corresponding relationship for the pan showed more scatter, but this was improved by using 5-day running means of evaporation and stratifying the data in terms of wind speed to yield a family of curves. Values for at wet soil conditions varied from 1.07 for 100 km day^–1 wind run to 1.17 for 250 km day^–1 wind run. For each curve, values of increased by about 20%; as the soil dried. The relationships may be used to reduce observed Class A pan evaporation to equivalent values for wet-soil conditions and to estimate near-surface soil moisture and actual evapotranspiration for this particular site. Extension of the technique to other areas requires derivation of similar relationships appropriate for those other locations     

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.