Measurements of solar radiation and estimation of optical depth in the high Andes of Peru

Summary During an expedition to the high Andes of Southern Peru in June–July 1977, measurements of direct solar radiation in four spectral bands (0.270–0.530–0.630–0.695–2.900 ) were conducted at six sites in elevations ranging from sea level to 5645 m. These measurements were evaluated in Langley plots to determine total optical depths () and irradiances at the top of the atmosphere. In addition, water vapor optical depths (_wv) were calculated from the mean radiosounding over Lima during the expedition, and Rayleigh (_ray) and ozone (_oz) optical depths were obtained from published tabulations. Subtracting _ray, _oz, and _wv from yielded estimates of aerosol optical depth _aer. The components _ray and _oz decrease from the shorter towards the longer wavelength bands and from the lower towards the higher elevation sites; _aer also decreases towards the higher elevations. Particularly pronounced is the decrease of _aer and from the lowlands of the Pacific coast to the highlands of the interior, reflecting the effect of a persistent lower-tropospheric inversion and the contrast from the marine boundary layer to the clear atmosphere of the high Andes.With 4 Figures     

Direct Measurement Of Shear Stress During Snow Saltation

Wind-tunnel experiments were carried out to measure the shear stressduring snow saltation. Shear stress acting on the snow surface, measured directly with a newly developed drag meter system, revealedthat the shear stress increased with the development ofsaltation. This result supports Owen's hypothesis that the saltationlayer acts as an increased roughness to the flow above the saltationlayer, leading to an increase in surface shear stress. To investigate the contribution of the grain borne shear stress_g and the fluid shear stress _f to the increment of the total shear stress _total, _g was calculated from the loss of horizontal momentum of saltating snowparticles. Since _g is the largest contribution to theincrement of _total, the collision of thesaltating particles is dominant for the shear stressmodification. The results qualitatively support the numericalsimulation reported by McEwan and Willetts.     

Second-order closure study of the covariance between chemically reactive species in the surface layer

A second-order modelling technique is used to investigate the influence of turbulence on chemical reactions. The covariance and variance equations for the NO-O₃-NO₂ system are developed as a function of the ratio of the timescale of turbulence ( _t ) and the timescale of chemistry (_Ch): the first Damköhler number ( _t /_Ch). Special attention is given to the calculation of the covariance between NO and O₃ normalized by the product of their means, the so-called intensity of segregation (I _S ). This parameter quantifies the state of mixing of two chemical species.The intensity of segregation is calculated as a function of the flux of NO and the first Damköhler number. The model results presented illustrate the importance of taking the effect of turbulence on chemical reactions into account for higher values of the NO flux, for values of the ratio O₃/NO larger than 12.5 and for values of the ratio _t /_CH larger than 0.1. For such cases, the effective reaction rates are slower than if the chemical species are assumed to be uniformly mixed.     

Drag and drag partition on rough surfaces

An analytic treatment of drag and drag partition on rough surfaces is given. The aims are to provide simple predictive expressions for practical applications, and to rationalize existing laboratory and atmospheric data into a single framework. Using dimensional analysis and two physical hypotheses, theoretical predictions are developed for total stress (described by the square root of the canopy drag coefficient), stress partition (described by the ratio _S/ of the stress _s on the underlying ground surface to total stress ), zero-plane displacement and roughness length. The stress partition prediction is the simple equation _S/= 1/(1+), where = C_RC_S the ratio of element and surface drag coefficients. This prediction agrees very well with data and is free of adjustable constants. Other predictions also agree well with a range of laboratory and atmospheric data.     

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.     

Disagreements between gradient-diffusion and Lagrangian stochastic dispersion models,even for sources near the ground

It is well known that if turbulent mass convection is modelled as diffusion, errors result unless trajectories from the source (ath) to the point of observation (z _p ) comprise many statistically-independent segments (Taylor, 1921). We show that this is not guaranteed merely by the Lagrangian timescale () at the source being small (e.g., source at ground), but that a better criterion istmax[(h), (z _p )], wheret is a typical travel time toz _p .     

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     

Statistical properties of 49 years of rainfall rate events

Summary The paper presents a detailed investigation of magnitudes and properties pertaining to the large population of rainfall events recorded during 49 years at the Fabra Observatory in Barcelona. The study includes a statistical analysis of event durationT and rainfall quantityQ together with the statistical rainfall rate parameters: , ²(R) and maximum Rm within an event. The decorrelation time is also analysed. It is found thatQ, T, and can be well modelled by a lognormal distribution, but , ² and Rm are only so for a limited range of precipitation ratesR. The regression analysis between pairs of logarithms of the magnitudes investigated is generally good and a regression coefficient is often better than 0.9. Comparison with published work is also carried out. An attempt is made to discriminate between heavy and non-heavy rainfall rates, and the 50 mm/h threshold is used for the study. The twelve-monthly running average indicates that the rainfall amount has a small increasing trend over the fifty year period. However, this trend is reversed when considering heavy rains. Finally, the return period in years to exceeding a thresholdR within an event is also investigated and the distribution of the population of annual extremes is found to be Gumbel II.With 6 Figures     

The atmospheric tide as a continuous spectrum: Lunar semidiurnal tide in surface pressure

Summary The standard equations for the theory of atmospheric tides are solved here by an integral representation on the continuous spectrum of free oscillations. The model profile of back-ground temperature is that of the U.S. Standard Atmosphere in the lower and middle atmosphere, and in the lower thermosphere, above which an isothermal top extends to arbitrarily great heights. The top is warm enough to bring both the Lamb and the Pekeris modes into the continuous spectrum.Computations are made for semidiurnal lunar tidal pressure at sea level at the equator, and the contributions are partitioned according to vertical as well as horizontal structure. Almost all the response is taken up by the Lamb and Pekeris modes of the slowest westward-propagating gravity wave. At sea level, the Lamb-mode response is direct and is relatively insensitive to details of the temperature profile. The Pekeris mode at sea level has an indirect response-in competition with the Lamb mode-and, as has been known since the time of its discovery, it is quite sensitive to the temperature profile, in particular to stratopause temperature. In the standard atmosphere the Lamb mode contributes about +0.078 mb to tidal surface pressure at the equator and the Pekeris mode about –0.048 mb.The aim of this investigation is to illustrate some consequences of representing the tide in terms of the structures of free oscillations. To simplify that task as much as possible, all modifying influences were omitted, such as background wind and ocean or earth tide. Perhaps the main defect of this paper's implementation of the free-oscillation spectrum is that, in contrast to the conventional expansion in the structures of forced oscillations, it does not include dissipation, either implicity or explicity, and thus does not satisfy causality. Dissipation could be added implicity by means of an impedance condition, for example, which would cause up-going energy flux to exceed downgoing flux at the base of the isothermal top layer. To achieve complete causality, however, the dissipation must be modeled explicity. Nevertheless, since the Lamb and Pekeris modes are strongly trapped in the lower and middle atmosphere, where dissipation is rather weak (except possibly in the surface boundary layer), more realistic modeling is not likely to change the broad features of the present results.Symbols a earth's mean radius; expansion coefficient in (5.3) - b recursion variable in (7.4); proximity to resonance in (9.2) - c sound speed in (2.2); specific heatc _p in (2.2) - f Coriolis parameter 2sin in (2.2) - g standard surface gravity - h equivalent depth - i ; discretization index in (7.3) - j index for horizontal structure - k index for horizontal structure; upward unit vectork in (2.2) - m wave number in longitude - n spherical-harmonic degree; number of grid layers in a model layer - p tidal pressure perturbation; background pressurep ₀ - q heating function (energy per mass per time) - r tidal state vector in (2.1) - s tidal entropy perturbation; background entropys ₀ - t time - u tidal horizontal velocityu - w tidal vertical component of velocity - x excitation vector defined in (2.3); vertical coordinate lnp _*/p ₀ [except in (3.8), where it is lnp /p ₀] - y vertical-structure function in (7.1) - z geopotential height - A constant defined in (6.2) - C spherical-harmonic expansion coefficient in (3.6) - D vertical cross section defined in (5.6) and (5.9) - E eigenstate vector - F vertical-structure function for eigenstate pressure in (3.2) [re-defined with WKB scaling in (7.2)] - G vertical-structure function for eigenstate vertical velocity in (3.2) [re-defined with WKB scaling in (7.2)] - H pressure-scale height - I mode intensity defined in (8.1) - K quadratic form defined in (4.4) - L quadratic form defined in (4.4); horizontal-structure magnification factor defined in (5.11) - M vertical-structure magnification factor defined in (4.6) - P eigenstate pressure in (3.2); tidal pressure in (6.2) - R tidal state vector in (5.1) - S eigenstate entropy in (3.2); spherical surface area, in differential dS - T background molecular-scale (NOAA, 1976) absolute temperatureT ₀ - U eigenstate horizontal velocityU in (3.2); coefficient in (7.3) - V horizontal-structure functionV for eigenstate horizontal velocity in (3.2); recursion variable in (7.3) - W eigenstate vertical velocity in (3.2) - X excitation vector in (5.1) - Y surface spherical harmonic in (3.7) - Z Hough function defined in (3.6) - +dH/dz - (1––)/2 - Kronecker delta; Dirac delta; correction operator in (7.6) - equilibrium tide elevation - (square-root of Hough-function eigenvalue) - ratio of specific gas constant to specific heat for air=2/7 - longitude - - - background density ₀ - eigenstate frequency in (3.1) - proxy for heating functionq =c _P/t - latitude - tide frequency - operator for the limitz - horizontal-structure function for eigenstate pressure in (3.2) - Hough function defined in (6.2) - earth's rotation speed - horizontal gradient operator - ()₀ background variable - ()_* surface value of background variable - () value at base of isothermal top layer - Õ state vector with zerow-component - , energy product defined in (2.4) - | | energy norm - ()^* complex conjugate With 10 Figures     

On the determination of the height of the Ekman boundary layer

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

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.     

On the temperature coefficient of PCO 2 in seawater

The commonly reported temperature coefficient of P. the equilibrium partial pressure of CO₂, is (P/T) _A,C ,which is about 15 ppm/°C, or 5% of the atmospheric partial pressure of CO₂. This coefficient, however, applies only to deep water, not to surface water which can exchange CO₂ with the atmosphere. The coefficient (P/T) _A,C ,, where designates constancy of the sum of atmospheric and surface-ocean CO₂, is the appropriate value for air-sea exchange. Numerical values are mass-dependent because the depth of the exchanging ocean layer must be specified. For a 100-m surface layer, the value is ca. 1.5 ppm/°C, or 0.5% of ambient CO₂. Editor's Note:In view of the interdisciplinary importance of the carbon dioxide-climate problem, this note on seawater chemistry should be of interest to specialists beyond the discipline of ocean chemistry.     

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     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.     

A comparison of airborne eddy correlation and bulk aerodynamic methods for ocean-air turbulent fluxes during cold-air outbreaks

Four bulk schemes (LKB, FG, D and DB), with the flux-profile relationships of Liuet al. (1979), Francey and Garratt (1981), Dyer (1974), and Dyer and Bradley (1982), are derived from the viscous interfacial-sublayer model of Liuet al. These schemes, with stability-dependent transfer coefficients, are then tested against the eddy-correlation fluxes measured at the 50 m flight level above the western Atlantic Ocean during cold-air outbreaks. The bulk fluxes of momentum (), sensible heat (H), and latent heat (E) are found to increase with various von Kármán constants (k _M for k _H forH, andk _E forE). Except that the LKB scheme overestimates by 28% (46Wm^–2), on the average, the fluxes estimated by the four bulk schemes appear to be in fairly good agreement with those of the eddy correlation method (magnitudes of biases within 10% for , 17% forH, and 13% forE). The results suggest that the overall fluxes and surface-layer scaling parameters are best estimated by FG and thatk _H <k _E . On the average, the FG scheme underestimates by 10% (0.032N m^–2) andE by 4% (12Wm^–2), and overestimatesH by 0.3% (0.5W m^–2). The equivalent neutral transfer coefficients at 10 m height of the FG scheme compare well with some schemes of those tested by Blanc (1985).The relative importance of various von Kármán constants, dimensionless gradients and roughness lengths to the oceanic transfer coefficients is assessed. The dependence of transfer coefficients on wind speeds and roughness lengths is discussed. The transfer coefficients for andE agree excellently between LKB and FG. However, the ratio of the coefficient forH of LKB to that of FG, increasing with decreasing stability, is very sensitive to stability at low winds, but approaches the neutral value of 1.25 at high winds.     

Properties of spatial temperature derivatives in the atmospheric surface layer

A pair of parallel cold wires separated in either the vertical or lateral direction was used to obtain the three components _x, _y, _z of the temperature derivative in the streamwise, lateral and vertical directions, respectively. The average absolute skewness values of _x and _z are nonzero and approximately equal, while the skewness of _y is approximately zero. These results appear to be consistent with the presence of a large, three-dimensional organised structure in the surface layer. There is an apparent low-frequency contamination in the spectral density of _y and _z due mainly to small errors in estimating the sensitivity of the cold wires. The temperature derivatives were high-pass filtered, the filter being set to remove possible contributions from the large structure and to minimise low-frequency sensitivity contamination. The filtered rms ratios \~_x/\~_y and \~_x/\~_z were in the range 0.7 to 0.9, a result in qualitative agreement with that obtained in the laboratory boundary layer by Sreenivasan et al. (1977). The skewness of filtered _x or _z is negligible, consistent with local isotropy of small-scale temperature fluctuations and in support of the high wavenumber spectral isotropy discussed in Antonia and Chambers (1978).     

A Simple Mixing Length Formulation for the Eddy-Diffusivity Parameterization of Dry Convection

In this paper a simple mixing length formulation for the eddy-diffusivityparameterization of dry convection is suggested. The new formulation relates the mixinglength to the square root of the turbulent kinetic energy (e) and a time scale ( ):l = e. To close the parameterization the time scale is calculated as a functionof the boundary-layer height (h) and the convective velocity scale (w_*), h/w_*. Thesimpler approach of a constant time scale is also studied. The simulation of a case of dry atmosphericconvection with a one-dimensional boundary-layer model shows that the model with the new formulationreproduces quite well the main properties of the convective boundary layer. In particular,the entrainment is realistically represented by the new mixing length, which has the advantage of naturallydecreasing with the turbulent kinetic energy. Sensitivity studies to the surface flux and the lapserate, in the context of a simplified situation, show the robustness of the new formulation.     

Microwave vegetation optical depth and inverse modelling of soil emissivity using Nimbus/SMMR satellite observations

Summary A radiative transfer model has been used to determine the large scale effective 6.6 GHz and 37 GHz optical depths of the vegetation cover. Knowledge of the vegetation optical depth is important for satellite-based large scale soil moisture monitoring using microwave radiometry. The study is based on actual observed large scale surface soil moisture data and observed dual polarization 6.6 and 37 GHz Nimbus/SMMR brightness temperatures over a 3-year period. The derived optical depths have been compared with microwave polarization differences and polarization ratios in both frequencies and with Normalized Difference Vegetation Index (NDVI) values from NOAA/AVHRR. A synergistic approach to derive surface soil emissivity from satellite observed brightness temperatures by inverse modelling is described. This approach improves the relationship between satellite derived surface emissivity and large scale top soil moisture fromR ²=0.45 (no correction for vegetation) toR ²=0.72 (after correction for vegetation). This study also confirms the relationship between the microwave-based MPDI and NDVI earlier described and explained in the literature.List of Symbols f frequency [Hz] - f _i(p) fractional absorption at polarizationp - h surface roughness - h h cos² - H horizontal polarization - n _i complex index of refraction - p polarization (H orV) - R _s microwave surface reflectivity - T _B(p) brightness temperature at polarizationp - T ^* normalized brightness temperature - T polarization difference (T _v-T _H) - T _s temperature of soil surface - T _c temperature of canopy - T _max daily maximum air temperature - T _min daily minimum air temperature - V vertical polarization - soil moisture distribution factor; also used for the constant to partition the influence of bound and free water components to the dielectric constant of the mixture - empirical complex constant related to soil texture - microwave transmissivity of vegetation (=e ^–) - ^* effective transmissivity of vegetation (assuming =0) - microwave emissivity - _s emissivity of smooth soil surface - _rs emissivity of rough soil surface - _vs emissivity of vegetated surface - soil moisture content (% vol.) - K dielectric constant [F·m^–1] - K _fw dielectric constant of free water [F·m^–1] - K _ss dielectric constant of soil solids [F·m^–1] - K _m dielectric constant of mixture [F·m^–1] - K _o permittivity of free space [8.854·10^–12 F·m^–1] - high frequency limit ofK _wf [F·m^–1] - wavelength [m] - incidence angle [degrees from nadir] - polarization ratio (T _H/T _V) - b soil bulk density [gr·cm^–3] - s soil particle density [gr·cm^–3] - R surface reflectivity in red portion of spectrum - NIR surface reflectivity in near infrared portion of spectrum - _eff effective conductivity of soil extract [mS·cm^–1] - vegetation optical depth - _6.6 vegetation optical depth at 6.6 GHz - ₃₇ vegetation optical depth at 37 GHz - ^* effective vegetation optical depth (assuming =0) - single scattering albedo of vegetation With 12 Figures     

The Kolmogorov constant for carbon dioxide in the atmospheric surface layer over a paddy field

From measurements in the atmospheric surface layer over a paddy field, the Kolmogorov constants for CO₂ and longitudinal wind velocity were obtained. In this study, the nondimensional dissipation rate _nc = (1–16 _v )^-1/2 for CO₂ variance and = (1–16 _v )^-1/4 – _v for turbulent energy were used, assuming the equality of the local production term and the local dissipation term, and neglecting the divergence flux term in the budget equation. The value of the constant for CO₂ was consistent with recent determinations for temperature and humidity. The constant for longitudinal wind velocity showed good agreement with other recent observations.     

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.