The Askervein Hill Project: Mean wind variations at fixed heights above ground

This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.     

The behaviour of aerosol optical depth during clear summer days in Athens,Greece

Summary The integral aerosol optical depths (k ) at the hour of 08:20 Local Standard Time (LST), are compared with those calculated previously at 11:20 and 14:20 LST, for clear days during summer in Athens over the period 1962–1988. The mean values at 08:20 LST were consistently lower than the values at 11:20 and 14:20 LST. The influence of the vertical wind profile on the values ofk was also investigated. A comparison was made of the wind profiles at 02:00 and 14:00 LST, for days in which the 11:20 and 14:20 LST values ofk were 0.200 andk 0.350, respectively. The corresponding bulk wind shear s was also found for the period 1980–1988. The most significant results occurred with the first category of days. The resultant wind velocities from the surface to the 900 hPa level, in each hour were higher by 2–4 m·s^–1 with respect to the corresponding values for the second category. At 02:00 LST the bulk wind shear showed a considerable difference (1.8) between the two categories of days in the surface to 700 hPa layer at 02:00 LST. Finally, the associated weather conditions that appear to initiate a period of low values ofk (k 0.200) at 11:20 and 14:20 LST were examined for the period 1980–1988. Fifteen such cases were identified and it was found that they all occurred after the passage of weak cold fronts.With 6 Figures     

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     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.     

A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.     

Effect of finite sampling on atmospheric spectra

The effect of a finite averaging time on variances is well known, but its effect on power spectra is less clearly understood. We present numerical solutions for the spectral distortion arising from sampling over a finite time interval T and show that the commonly used filter function (1 – sinc²f T), valid for variances, is a reasonable approximation for power spectra only when T 10 _m , where f is the cyclic frequency, and _m is the dominant time scale of the process. Our results exhibit an increasingly steeper low-frequency roll-off as T decreases relative to _m , indicating that the measured spectrum is subject to a greater suppression of the lower frequencies (f > 1/T) than predicted by (1 – sinc²f T). This suppression is, in a sense, compensated by an overestimation of spectral estimates in the frequency range f 1/T.     

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.     

The effects of orography on precipitation

In this paper, a model simulating the effects of topography and altitude on precipitation is presented. Topography has its maximum effect on precipitation when the angle which the wind makes with the slope direction approaches zero and the inclination of the slope is near 45°. The smaller the angle , the greater the influence of slope on precipitation. When <45°, the larger the inclination, the greater the influence of slope on precipitation and the less the difference in precipitation between the windward and the leeward slopes. When <45°, the reverse holds. But for in the range of 0°–45° and in the range 30°–60°, differences in precipitation on both the windward and leeward slopes are not so well marked and can be neglected in general. In condition of uniform slope inclination, precipitation on the windward slope increases with altitude at first and then decreases after attaining a height (H _m ) of maximum precipitation; alsoH _m is greater, the drier the air mass. When the terrain on the windward side is stepped in shape, it is possible that more than one height of maximum precipitation will occur.     

A new theory of the energy spectrum

In a recent paper, the author introduced a new viscous boundary layer, called the mesolayer, in turbulent shear flow. Its importance stems from its location between the inner and outer regions which are controlled by the law of the wall and Reynolds number similarity, respectively. This intrusion prevents the classical overlap assumption which appears to be fundamental in the derivation of the classical logarithmic behavior. The mesolayer has a thickness proportional to Taylor's microscale . This, and the analogy between the energy equation for the spectrum function of isotropic turbulence and the momentum equation for shear flow, suggest the existence of a similar region in wavenumber space with wavenumber k ~ ^-1. This mesoregion separates the inner region k ~ k _s(where k _s^-1 and is the Kolmogorov length) and the outer region k k _e(where k _e ^-1 and l is the energy-containing eddy size) and again invalidates the overlap assumption which appears to be fundamental in the derivation of the classical k ^-5/3-behavior of the inertial subrange.Incorporation of the mesoregion into the argument leads to a new theory with k ^-5/3-behavior in two regions (^-1 k k _s) and (k _e k ^-1) although with two different coefficients of proportionality (Kolmogorov constants). This leads to a wandering of the spectrum curve about the classical k ^-5/3 line similar to a wandering in turbulent shear flow about the logarithmic curve. This is clearly indicated by the data for the variation of the Kolmogorov constant.Other data support the new theory. In particular, the location of the point k _mwhere the curve of the nonlinear energy-transfer function goes through zero shows agreement with the theory, i.e., k _m^-1.     

Similarity functions D for water vapor in the unstable atmospheric boundary layer

An analysis was performed of experimental data obtained at fixed ship stations during AMTEX 1974 and 1975. This allowed the calculation of the bulk transfer relationships for water vapor and sensible heat in the atmospheric boundary layer for different interpretations of the thickness scale of the boundary layer. It was found that scaling based on the observed thickness, which herein was taken as the height of the lowest value in the potential temperature profile under unstable conditions, produces least scatter in the calculations. The results obtained for the similarity function c( _i ) of the bulk heat transfer coefficient are in good agreement with the results of previous studies. As observed earlier (Brutsaert and Mawdsley, 1976; Mawdsley and Brutsaert, 1977), under unstable conditions the similarity functions D() of the bulk water vapor transfer coefficient are smaller than the corresponding C() functions for sensible heat. In the case of inversion height scaling, the results can be represented by d( _i ) = 0.65 c( _i ).     

An application of the E-ε turbulence model for studying coupled air-sea boundary-layer structure

An E- turbulence model is used to study air-sea interaction characteristics and turbulence structure using a coupled model for air-sea boundary layers. The E- turbulence model consists of equations for the turbulent kinetic energy, the energy-dissipation, and for the turbulent exchange coefficient expressed in terms of turbulent kinetic energy and energy-dissipation. The energy-dissipation equations for the air-sea interface are solved analytically to obtain boundary conditions for energy-dissipation at the interface. The air-sea interaction and turbulence characteristics of the E- model are compared with those of the mixing-length model and with available observations.The simulations demonstrate that the air-sea interaction parameters obtained by the E- model agree well with observations. The numerical studies also show that the E- turbulence model with appropriate constants can give good results in modeling coupled air-sea boundary-layer flows.     

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.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

Turbulent characteristics over a rough natural surface part I: Turbulent structures

A comprehensive experiment for turbulence and profiles has been carried out over a very rough natural mallee bushland. In Part I, the large inherent structures of the turbulence are presented. Short-period space-time correlations of wind and temperature indicate that periodic fluctuations pervade the lower surface layer above the canopy. These are associated with the appearance of large turbulent structures. Prominent peaks in the spectra are identified with the arrival of these structures, which show that they have a relatively large length scale over rough surfaces. These findings suggest a different mechanism for the enhancement of turbulent diffusivity over rough surfaces from that by the simple generation of turbulence by plant wakes.Analysis via conditional sampling of the turbulent components has revealed that bursts make almost the same contributions as gusts (at the measurement height of 8.4 m above ground) for either momentum or temperature under neutral conditions; for unstable conditions and high hole size H (see Equation (3)), gusts dominate for momentum while bursts dominate for heat transfer. This implies that the mechanisms for momentum and heat transfer are different.     

A Field Study of the Mean Pressure About a Windbreak

To provide additional field data for assessingwindbreak flow models, mean ground-level pressurehas been measured upstream and downstream from along porous fence (height H = 1.25 m, resistancecoefficient k _r = 2.4). Measurements were madeduring periods of near-neutral stability and near-normallyincident flow, with the fence standing on bare soil(roughness length, z ₀ 0.8 cm;H/z ₀ 160), or within a plant canopy. The mean pressure field,measured far from the ends of the fence, was foundto be quite insensitive to mean wind direction( , zero for perpendicular flow), for| | less than about 25°.In the absence of a canopy, during each measurementperiod the minimum pressure occurred at the closestsampling location to leeward of the windbreak, thepressure-gradient in most cases beingmaximally-adverse in the immediate lee, and decayingwith increasing downwind distance (x). On one day ofmeasurements, however, the pressure gradient over2 x/H 6 (H = windbreak height) resembled theleeward plateau identified by Wang and Taklein their numerical studies. Perhaps thisoccasional feature was only due to instrumenterror. Nevertheless a plateau of sorts wasindicated in similar measurements by Judd andPrendergast (with H = 1.92 m, z ₀ 1.2 cm;H/z ₀ 160, k _r 3). Therefore,existence of a leeward pressure plateau behind athin fence cannot be definitely ruled out.When the windbreak was placed in a canopy, minimumsurface pressure was displaced downwind. Thisagrees with the wind-tunnel study of Judd, Raupach and Finnigan,and is consistent with a simple simulation reported here.     

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.     

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.     

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.     

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.     

The carbon isotopic ratio of atmospheric carbon dioxide at Tsukuba,Japan

To find out the secular and seasonal trends of the ¹³C value and CO₂ concentration in the surface air and the determination of the ¹³C in the atmospheric CO₂ collected at Tsukuba Science City was carried out during the period from July 1981 to October 1983. The monthly average of the ¹³C value of CO₂ in the surface air collected at 1400 LMT ranged from -7.52 to \s-8.45 with an average of -7.96±0.25 and the CO₂ concentration in the air varied from 334.5 l 1^-1 to 359 l 1^-1 with an average of 347.2±6.3 l 1^-1. The ¹³C value is high in summer and low in winter and is negatively correlated with the CO₂ concentration. In general, the relationship between the ¹³C and the CO₂ concentration is explainable by a simple mixing model of two different constant carbon isotopic species but the relationship does not always follow the model. The correlation between the ¹³C value and the CO₂ concentration is low during the plant growth season and high at other times. The observed negative deviation of the ¹³C value from the simple mixing model in the plant growth season is partly due to the isotopic fractionation process which takes place in the land biota.