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.     

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     

The use of σ T as a temperature scale in the atmospheric surface layer

The standard deviation of temperature _T is proposed as a temperature scale and as a velocity scale to describe the behaviour of turbulent flows in the Atmospheric Surface Layer (ASL), instead of _* andu _* of the Monin—Obukhov similarity theory, and ofT _f andU _f used for free convection stability conditions. On the basis of experimental evidence reported in the literature, it is shown that _T T _f andv _* U _f in the free convection region, and _{T *} andv _* U _* in nearneutral and stable conditions. This implies that the proposed scales can be applied for all stabilities. Furthermore, a new length scale is proposed and its relation with Obukhov length is given. Also, a simple semi-empirical expression is presented with which _T andv _* can be evaluated in a rather simple way. Some examples of practical applications are given, e.g., a stability classification for unstable conditions.     

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     

On Modelling the One-Dimensional Atmospheric Boundary Layer

Several commonly used turbulence closure schemes forthe atmospheric boundarylayer (ABL) are applied to simulate neutral, nocturnal and diurnal cyclesituations in a one-dimensional ABL. Results obtained with the differentschemes, E-, E- and its modified versions, and twoversions ofthe q ² Level 2.5, are compared and discussed.     

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     

Flux Footprints Within and Over Forest Canopies

The characteristics of turbulence within a forest arespatially heterogeneous and distinct from thoseassociated with the surface boundary layer. Consequently, the size and probability distribution offlux footprints emanating from sources below aforest canopy have the potential to differ from thoseobserved above forests.A Lagrangian random walk model was used toinvestigate this problem since no analytical solutionof the diffusion equation exists. Model calculations suggest that spatialcharacteristics of flux footprint probabilitydistributions under forest canopies are muchcontracted, compared to those evaluated in the surfaceboundary layer. The key factors affecting thestatistical spread of the flux footprint, and theposition of the peak of its probability distribution,are horizontal wind velocity and the standarddeviations of vertical and horizontal velocityfluctuations. Consequently, canopies, which attenuatemean horizontal wind speed, or atmospheric conditions,which enhance vertical velocity fluctuations, willcontract flux footprint distributions mostly near thefloor of a forest. It was also found that theprobability distributions of the flux footprint arenarrower when horizontal wind velocity fluctuationsare considered, instead of the simpler case that considers only vertical velocity fluctuations and meanhorizontal wind velocity.     

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.     

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.     

On the similarity functions A and B as determined from the ‘VOVES’ experiment

Data collected during the french VOVES-1977 experiment are used for the determination of the functionsA() andB() of the stability parameter . It is found that, although pertaining to different geographic conditions, they lead to the same kind of dependency ofA andB upon as found previously considered data. The influence of baroclinicity onA andB is also studied and it is found in that, statistically, consideration of this parameter does not lead either to a decrease of the large scatter of experimental points or to a significant improvement of the similarity theory. It is also shown that a satisfactory estimation of the surface geostrophic wind can be made from the network of meteorological stations, to compute such functions.     

On the representation of frequency spectra in meteorology

The common representation of frequency spectra in meteorology and climatology is discussed. It is pointed out that this representation is misleading since spectral peaks and spectral gaps are obtained even when the spectrum density is monotonously decreasing in the whole frequency range. A plea is made for using the spectrum distribution function, F() which gives an unambiguous picture of the distribution of variance with frequency.     

Über die mittlere lineare sukzessive Differenz und ihre Anwendung auf die Untersuchung meteorologischer Beobachtungsreihen

Zusammenfassung Eine endliche Reihe (Sequenz) wird als eine der möglichen Permutationen ihrer Glieder aufgefaßt. Es wird gezeigt, daß die Summe der absoluten Differenzen der aufeinanderfolgenden Glieder gleich ist , wo die natürliche Zahlen sind und nur von der Rangordnung der Glieder der Reihe (von der Permutation) abhängen; die _j sind von der Reihenfolge unabhängig und werden durch die Dispersion der Reihenglieder bestimmt. Die _j und die _j werden separat untersucht; der Erwartungswert der erwähnten Differenzsumme wird abgeleitet. Verschiedene bereits bekannte und auch erstmalig hier vorgeschlagene Maßzahlen werden geprüft. An Reihen jährlicher Regenmengen wird die Rolle der _j und der _j und das Verhalten der besprochenen Maßzahlen veranschaulicht.

---

Summary A series ofn members can be considered as one of the possible permutations of its members. It is shown that the sum of the linear successive differences is equal to the expression , where the _j are positive integers, dependent only upon the rank-order (the permutation) of the members, while the _j are independent of the order of the succession and are determined by the dispersion of the members of the series. The factors _j and _j are separately investigated; the expected value of the sum of the linear successive differences is established. Various related statistical measures, already in usage and new ones suggested here, are discussed. Series of yearly rainfall amounts are used to show the effects of the _j and _j and to discuss the behaviour of the various measures.

---

Résumé Une série, constituée parn valeurs, est regardée comme une des possibles permutations de ces valeurs. L'auteur montre que la somme des différences absolues, qui se présentent entre les valeurs consécutives de la série, est égale à l'expression . Les _j sont des nombres entiers positifs et ne dépendent que de l'ordre des membres de la série, tandis que les _j, indépendants de l'ordre, sont déterminés par la dispersion des membres. Les facteurs _j et _j sont étudiés séparément; l'espérance mathématique de la somme mentionnée est dérivée. Des paramètres statistiques déjà connus ou proposés ici pour la première fois, sont discutés. Le rôle des _j et des _j et le comportement des divers paramètres sont montrés à l'aide de séries de totaux annuels de pluies.

     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

MCC and gull flight behavior

Wintertime observations of mesoscale cellular convection (MCC) over the East China Sea have resulted in criteria that have a remarkable similarity to those reported by Woodcock (1975) in the study of thermals and gull flight behavior. It has been determined that the surface wind speed (V) and the air-sea temperature difference (T) prescribe unique and compatible conditons for both the occurrence of MCC and soaring by sea gulls. Specifically, the onset of MCC when V is between 5 and 9 m s^–1 is inversely proportional to T in the range 5 to 7 °C. Elsewhere, the onset of MCC occurs under conditions of direct proportionality between V and T. Necessary conditions for the occurrence of MCC due to heating from below are T 5 °C and V 5 m s^–1. The boundaries of the convective regime for MCC are discussed and interpreted in accordance with the regime for sea-gull soaring and similarity concepts.     

Diurnal variation of horizontal wind direction fluctuations in complex terrain at Geysers,Cal.

Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km² Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s^-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.     

The influence of surface cover and climate on energy partitioning and evaporation in a subarctic wetland

Energy partitioning and evaporation were measured over three wetland surfaces in a subarctic coastal marsh during pre-growing and growing periods. These surfaces included an alder/willow woodland, a sedge marsh and a raised backshore sedge meadow. A combination model analysis was used to assess the relative importance of surface resistance and meteorological conditions on the magnitude of the Bowen ratio, , during the growing period.Overall, the three surfaces experienced important site-to-site and seasonal differences in and evaporation, Q _E. During the non-foliated period, Q _E was largest and was smallest for the open water marsh, while the dry backshore site experienced the smallest Q _E and largest . The non-foliated woodland assumed intermediate values of and Q _E. After the vegetation covers were established, the woodland assumed the smallest and largest Q _E flux. It was also found that at the marsh site increased with the presence of a vegetation cover.Wind direction was always an important factor in determining Q _E and at all sites. was substantially larger and Q _E was smaller for onshore winds (i.e., originating from James Bay) than for offshore winds. The combination model analysis showed that canopy resistance at all sites was largest during warm offshore winds, which were associated with large saturation deficits. However, the effect of increased canopy resistance on during offshore winds was offset by a large climatological resistance, resulting in small values and large Q _E. When winds originated from James Bay, canopy resistance was smaller than for offshore winds, but the climatological resistance also was much smaller, resulting in larger and small Q _E. The results have important implications for changes in land cover and climate on the regional water balance.     

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.     

Equilibrium evaporation beneath a growing convective boundary layer

Expressions for the equilibrium surface Bowen ratio ( _s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer _i and the entrainment parameter A _R . If A_R is put equal to zero, the solution for _s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA _R and _i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.     

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.     

Field measurements of snow-drift threshold and mass fluxes,and related model simulations

Field measurements were carried out to calculate the threshold friction velocity for snow saltation, and mass fluxes during snow drift. The wind was measured in three components by an ultrasonic anemometer, and the mass fluxes were determined using an optical sensor (snow particle counter), acoustic sensors (Flowcapt) and mechanical traps. The threshold friction velocity was found to be correlated to the grain size (R²=0.75). The mass flux measurements were compared with numerical simulations of snow drift, and it was demonstrated that the maximum snow transport takes place at shear stress values of roughly two times the average shear stress over 20 min. By implementing a probability distribution for the shear stress the mass flux was simulated with only the mean measured value of the shear stress as input. This procedure enables the future use of the numerical model for operational applications.