Exponential-linear stability correction functions for weak to moderate instability near the ground

Stability correction functions which combine the exponent of z/L, and a linear term in z/L, are proposed for the unstable case. The functions provide a reasonably close fit to the _m and _h results of Dyer and Hicks (1970) for 0 < –z/L 1, but they cannot be extended to cases of strong instability. Attractive features are the ability to integrate the expressions directly in terms of z/L, and a particularly close fit of the integrated result to experimentally derived _m values.     

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     

Local Similarity Relationships In The Urban Boundary Layer

To investigate turbulent structures in an urban boundary layer (UBL) with many tallbuildings, a number of non-dimensional variable groups based on turbulent observationsfrom a 325-m meteorological tower in the urban area of Beijing, China, are analyzedin the framework of local similarity. The extension of surface-layer similarity to localsimilarity in the stable and unstable boundary layer is also discussed. According to localsimilarity, dimensionless quantities of variables: e.g., velocity and temperature standarddeviations _i/u_*l (i=u,v,w) and_T/T_*l,correlation coefficients of uw and wT covariance, gradients of wind and temperature_m and _h, and dissipation rates of turbulent kinetic energy (TKE) andtemperature variance and _N can be represented as a functiononly of a local stability parameter z/, where is the local Obukhovlength and z is the height above ground. The average dissipation rates of TKE andtemperature variance are computed by using the u spectrum, and the uw and wTcospectra in the inertial subrange. The functions above were found to be in a goodagreement with observational behaviour of turbulence under unstable conditions, butthere were obvious differences in the stable air.     

The spectral velocity tensor for homogeneous boundary-layer turbulence

We have postulated a simple model for the spectral tensor _ij (k) of an anisotropic, but homogeneous turbulent velocity field. It is a simple generalization of the spectral tensor _inf ^{ij piso}(k) for isotropic turbulence and we show how in the limit of isotropy, _ij (k) becomes equal to _inf ^{ij piso}(k). Whereas _inf ^{ij piso}(k) is determined entirely by one scalar function of k = ¦k¦, namely the energy spectrum, we need three independent scalar functions of k to specify _ij (k). We show how it is possible by means of the three stream-wise velocity component spectra to determine the three scalar functions in _ij (k) by solving two uncoupled, ordinary linear differential equations of first and second order. The analytic form of the component spectra each has a set of three parameters: the variance and the integral length scale of the velocity component and a dimensionless parameter, which governs the curvature of the spectrum in the transition domain from the inertial subrange towards lower wave numbers. When the three sets of parameters are the same, the three spectra correspond to isotropic turbulence and they are all interrelated and related to the energy spectrum. We show how it is possible to obtain these spectral forms in the neutral surface layer and in the convective boundary layer from data reported in the literature. The spectral tensor is used to predict the lateral coherences for all three velocity components and these predictions are compared with coherences obtained in two experiments, one using three masts at a horizontally homogeneous site in Denmark and one employing two aircraft flying in formation over eastern Colorado. Comparison shows reasonable agreement although with considerable experimental scatter.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Flux–Profile Relationships Over a Fetch Limited Beech Forest

The influence of an internal boundary layer and a roughness sublayer on flux–profile relationships for momentum and sensible heat have been investigated for a closed beech forest canopy with limited fetch conditions. The influence was quantified by derivation of local scaling functions for sensible heat flux and momentum (_h and _m) and analysed as a function of atmospheric stability and fetch. For heat, the influences of the roughness sublayer and the internal boundary layer were in agreement with previous studies. For momentum, the strong vertical gradient of the flow just above the canopy top for some wind sectors led to an increase in _m, a feature that has not previously been observed. For a fetch of 500 m over the beech forest during neutral atmospheric conditions, there is no height range at the site where profiles can be expected to be logarithmic with respect to the local surface. The different influence of the roughness sublayer on _h and _m is reflected in the aerodynamic resistance for the site. The aerodynamic resistance for sensible heat is considerably smaller than the corresponding value for momentum.     

Stability functions correct at the free convection limit and consistent for for both the surface and Ekman layers

For the thermal stability function _h used to calculate heat and moisture fluxes in the surface layer, we choose a formulation which has the theoretically correct free convection limit % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeikaiabgk% HiTGqaciaa-PhacaqGVaGaamitaiaabMcadaahaaWcbeqaaiabgkHi% TiaaigdacaGGVaGaaG4maaaaaaa!3DFE!\[{\rm{(}} - z{\rm{/}}L{\rm{)}}^{ - 1/3} \]. We then use the experimental result that z/L Ri to deduce a formulation with an exponent -1/6 for the momentum stability function _m. This formulation also resolves the matching problem at the interface between the surface and Ekman layers. The proposed functions are found to remain reasonably close to another formulation that is well supported by observations and has exponents -1/2 for _h and -1/4 for _m. The intent of the proposals is mainly to clarify and simplify the parameterization of the convective boundary layer in present day atmospheric models, without significantly altering the results.     

Estimation of the Monin-Obukhov similarity functions from a spectral model

From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function _m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - C_i constants i=u, v, w, (see Equation (5) - k_me,k_mT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - k_mi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - k_sb, k_c characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - _m nondimensional mean shear - _H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - _D nondimensional divergence of vertical transports of turbulent linetic energy     

Analytical solutions for the Ekman layer

The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U _*z (0.4,U _*=friction velocity,z=elevation) to the geostrophic asymptoteK _GU _*d forzd is solved through an expansion in fd/U _*1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU _*d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U _*d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U _* (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U _* and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U _*d[1–exp(–z/d)] are developed in appendices.     

Mass transfer to bean leaves

A vapour of radio-lead (²¹²Pb) has been used to measure the Sherwood number, Sh, of model leaves at various angles of incidence,, to the airstream in a wind tunnel. The results for=0 are compared with Pohlhausen's formula and the results for 0, with Powell's experiments. The local values of Sh on the upwind and downwind sides of discs have been obtained. For leaves in the canopy, Sh was found to be about 25% greater than would be predicted by applying Pohlhausen's equation without correction for orientation.     

An observational study of heat fluxes and their relationships with net radiation

Experimental evidence indicates that the diurnal behaviour of the fluxes of heat into the ground and into the atmosphere versus the net flux of radiation can be modelled by closed curves, the hourly values folowing one another in either a clockwise or counter clockfashion. A general formulation to express the different heat fluxes as a function of net radiation is proposed. This relationship between the different heat fluxes and can be expressed as a sum of three terms: the first indicates a direct proportionality to , the second gives the deviation from linearity and depends on /t, and the third gives the value of the flux when = 0. The formulae are then expressed as a function of time and the ratios between the heat fluxes and are evaluated. A comparison with the approximations generally used shows that the latter may be considered as particular cases of the more general equations proposed here.     

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.     

Turbulence Reynolds number and the turbulent kinetic energy balance in the atmospheric surface layer

The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .     

Non-dimensional wind and temperature profiles in the atmospheric surface layer: A re-evaluation

Previous results of non-dimensional wind and temperature profiles as functions of ( = z/L) show systematic deviations between different experiments. These discrepancies are generally believed not to reflect real differences but rather instrumental shortcomings. In particular, it is clear that flow distortion has not been adequately treated in most previous experiments. In the present paper, results are presented from a surface-layer field experiment where great care was taken to remove any effects from this kind of error and also to minimize other measuring errors. Data from about 90 30-min runs with turbulence measurements at three levels (3, 6, and 14 m) and simultaneous profile data have been analysed to yield information on flux-gradient relationships for wind and temperature.The flux measurements themselves show that the fluxes of momentum and sensible heat are constant within ± 7% on average for the entire 14 m layer in daytime conditions and when the stratification is slightly stable. For more stable conditions, the flux starts to decrease systematically somewhere in the layer 6 to 14 m. From a large body of data for near-neutral conditions (¦¦ 0.1), values are derived for von Kármán's constant: 0.40 ± 0.01 and for _h at neutrally, 0.95 ± 0.04. The range of uncertainty indicated here is meant to include statistical uncertainty as well as the effect of possible systematic errors.Data for _m and _h for an extended stability range (1 > > – 3) are presented. Several formulas for _m and _h appearing in the literature have been used in a comparative study. But first all the formulas have been modified in accordance with the following assumptions: = 0.40 and ( _h ) _{= 0} = 0.95; deviations from this result in the various studies are due to incomplete correction for flow distortion. After new corrections are introduced, the various formulas were compared with the present measurements and with each other. It is found that after this modification, the most generally used formulas for _m and _h for unstable conditions, i.e., those of Businger et al. (1971) and Dyer (1974) agree with each other to within ± 10% and with the present data. For stable conditions, the various formulas still disagree to some extent. The conclusion in relation to the present data is not as clear as for the unstable runs, because of increased scatter. It is, however, found that the modified curve of Businger et al. (1971) for _h fits the data well, whereas for _m , Dyer's (1974) curve appears to give slightly better agreement.     

A dynamical analysis of the drag conditions at the sea surface

When applied to a sea surface, shortcomings are noted for the ordinary classification of drag conditions at rigid underlying surfaces according to the Reynolds roughness number Re _s . It is shown that in the case of mobile underlying surfaces, it would be more natural to use the dynamical classification of drag conditions according to the order of magnitude of the ratio ( = /) of the momentum flux toward the waves ( _w) to the viscous momentum flux through the surface ( _w). The relevant estimates of for the main stages of development of the wind waves indicate that the observed values of the drag coefficient of the sea surface correspond to the case of underdeveloped roughness.     

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.     

Modelling radiation quantities and photolysis frequencies in the troposphere

STAR (System for Transfer of Atmospheric Radiation) was developed to calculate accurately and efficiently the irradiance, the actinic flux, and the radiance in the troposphere. Additionally a very efficient calculation scheme to computer photolysis frequencies for 21 different gases was evolved. STAR includes representative data bases for atmospheric constituents, especially aerosol particles. With this model package a sensitivity study of the influence of different parameter on photolysis frequencies in particular of O₃ to Singlet D oxygen atoms, of NO₂, and of HCHO was performed. The results show the quantitative effects of the influence of the solar zenith angle, the ozone concentration and vertical profile, the aerosol particles, the surface albedo, the temperature, the pressure, the concentration of NO₂, and different types of clouds on the photolysis frequencies.Notation I _A(, ) actinic flux - I _H(, ) irradiance - L(, , , ) radiance - wavelength - azimuth angle - cosine of zenith angle - _s cosine of solar zenith angle - optical depth - _s scattering coefficient - _c extinction coefficient - _o single scattering albedo - p _mix mixed phase function - g _mix mixed asymmetry factor - J _gas photolysis frequency     

The impact of turbulence closure schemes on predictions of the Mixed Spectral Finite-Difference model for flow over topography

Impacts of different closure schemes in the Mixed Spectral Finite-Difference model (Beljaarset al., 1987) for neutrally stratified atmospheric surface-layer flow over complex terrain are studied. Six different closure schemes, (Z+z ₀), Mixing Length,E–(Z+z ₀),E–,E–– andq ² l are compared. Model results for flow over an infinite series of sinusoidal ridges are examined in the context of the inner and outer layers defined by Jackson and Hunt (1975). Results are compared with rapid distortion estimates of the changes in normal stresses. The effects of streamline curvature are also examined in a qualitative sense.     

AN E – ε – ℓ TURBULENCE CLOSURE SCHEME FOR PLANETARY BOUNDARY-LAYER MODELS: THE NEUTRALLY STRATIFIED CASE

The standard E – model generates aplanetary boundary layerthat appears to be much too deep. The cause of theproblem is traced to the equation for the dissipationrate () of turbulent kinetic energy (E), specifically theparameterization of dissipation production anddestruction. In the context of atmosphericboundary-layer modelling, we argue that a part of thedissipation production should be modelled as the inputto the spectral cascade from the energy-containingpart of the spectrum, with a characteristic length , while the equilibrium imbalancebetween local production and destruction ofdissipation is modelled as proportional toE2/E, as in the standard model. Wepropose an E – – turbulence closurescheme, in which both the mixing length, m, and are prescribed. The importance ofthe equation is diminished, though itstill determines the dissipation rate in the Eequation.     

The response of a propeller anemometer to turbulent flow with the mean wind vector perpendicular to the axis of rotation

The system transfer function ¦H(v)¦² at frequencyv (units of Hz) for a vertical velocity propeller anemometer in a statistically stationary and horizontally homogeneous turbulent flow is determined from: (1) experimental estimates of propeller velocity spectra; and (2) estimates of Eulerian vertical velocity spectra based on the hypothesis that degradation of the input vertical velocity Fourier components occurs in the inertial subrange. The experimental estimates of ¦H(v)¦² were adequately summarized with the mathematical expression for the system transfer function of a first-order system with parameterT which has units of time and is analogous to the time constant of a horizontal velocity propeller anemometer. Dimensional analysis techniques and the Monin-Obukhov similarity hypothesis were used to construct a model for the system parameterT which yielded the result that _w /D ₁ ( _w /)^1/3, where _w , andD ₁ denote the standard deviation of the input vertical velocity fluctuations, the horizontal mean wind speed, and the diameter of the propeller, respectively. The system parameterT is interpreted in terms of the time required for the propeller velocity statistics to become asymptotically independent of time upon being released from rest in a statistically stationary turbulent flow.Currently on leave of absence from the Indian Institute of Technology, New Delhi, India.     

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