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     

Wind and wind forces in a plantation spruce forest

Observations have been made of the tubulent structure within and above a dense (LAI=10.2) plantation spruce forest along with measurements of the movement of individual trees. The mean statistics of the turbulence and the turbulence spectra are compared with observations in other crops and complementary wind-tunnel studies using 175 scale plastic trees. The measurements show that momentum transport and the subsequent motion of the trees is dominated by intermittent sweep/ejection events associated with honami waves moving across the forest. The trees themselves act as forced damped harmonic oscillators and appear to short circuit the normal turbulent energy dissipation process by efficiently absorbing energy at their resonant frequencies. It is argued that understanding the nature and formation of honami waves over forests and crops is a crucial problem in agricultural and forest meteorology because of their important role both in turbulent transport and in causing wind damage.     

Shear convection

The paper discusses convection in the presence of wind shear, a condition analysed previously by Zilitinkevich (1971). This region (between the forced and fully developed convective layers) was also considered by Betchov and Yaglom (1971). In the present paper the author endeavors to develop a consistent analysis from the basic hypothesis of a very weak interaction between the vertical convective motions and mechanical turbulence, employing a new similarity model of the turbulent regime. Additional experimental data are introduced. Unlike the notation used in the references quoted above, this regime is termed shear convection rather than free convection. The latter is traditionally regarded as synonymous with the terms pure or fully developed convection.     

Charakteristische Zirkulationstypen in mittleren Breiten der nördlichen Hemisphäre

Zusammenfassung In der vorliegenden Arbeit werden neue Klassifikationsprinzipien für Großwetterlagen entwickelt. Bisher wurde bei Wetterlagenklassifikationen das Druckfeld zugrunde gelegt, wobei quasistationären Druckzentren eine nicht berechtigte Vorrangstellung eingeräumt wurde. In der hier versuchten Klassifikation wird vom Strömungsfeld ausgegangen, das in elementare Formen zerlegt wird. Eine zu diesem Zweck durchgeführte statistische Untersuchung ergab, daß alle im Strömungsfeld auftretenden Zirkulationstypen auf drei Grundformen zurückgeführt werden können. Diese Grundformen sind:Driften, Wellen undWirbel.Die Untersuchung ergab im einzelnen, daß in mittleren Breiten der nördlichen Hemisphäre bei 49% aller untersuchten Fälle Driften, bei 23% Wellen und bei 28% Wirbel auftraten.In der hier durchgeführten Klassifikation wird das Druckfeld durch das Strömungsfeld und der Begriff Großwetterlage durch den umfassenderen Begriff des Zirkulationstyps ersetzt. Damit wird der unberechtigte Vorrang der Druckformen bei der Wetterlagenklassifikation aufgegeben. Die Klassifizierung der Zirkulationstypen ergibt sich schließlich durch Kombination der drei Zirkulationselemente: Drift, Welle und Wirbel.

---

Summary New principles of classification for large-scale weather situations are outlined in this paper. Hitherto the pressure-field has been taken as a basis for such classifications by conceding a precedence of an unjustified position to semi-permanent centres of pressure. The new classification starts from the field of large-scale motions, which is dissected in elementary models. A statistical test yielded the possibility to reduce all types of atmospheric circulations in the following three elementary models:drifts, waves andeddies.In detail it was found out, that drifts occur in 49%, waves in 23% and eddies in 28% of all cases investigated.In the new classification the term pressure-field is substituted by field of motion and the expression large-scale weather situation by the more comprehensive conception type of circulation. By that the unjustified priority of pressure-centers in classifying weather situations is abolished. At last the classification of the types of circulation follows from a combination of the three elementary models: drift, wave and eddy.

---

résumé La présente étude développe de nouveaux principes de classification des situations météorologiques. Alors que jusqu'ici on s'est fondé sur le champ de pression ce qui conduisait à attribuer aux centres d'action quasi stationnaires un rôle trop important, l'auteur part ici du champ de mouvement décomposé en formes élémentaires. Un examen statistique lui a montré que tous les types de circulation peuvent se ramener à trois formes fondamentales:courants, ondulations ettourbillons.Aux latitudes moyennes de l'hémisphère Nord les courants représentent le 49%, les ondulations le 23% et les tourbillons le 28%.Au champ de pression se substitue donc le champ de courant, et les situations météorologiques se groupent en types de circulation ce qui supprime le rôle prépondérant des formes isobariques. Le classement final des types de circulation résulte de la combinaison des trois types mentionnés:courants, ondulations ettourbillons.

---

---

Mit 15 Textabbildungen     

Stable boundary layers: observations,models and variability part II: Data analysis and averaging effects

As argued in Part I (Derbyshire, 1995), variability is a key issue in stable boundary layers, and differences in variability between observations and idealized models may imply sytematic biases. Here we discuss how data analysis can be geared to allow for variability and thus consistency with models. Instrumental errors, smoothing methods and vertical discretization are considered. We then show how statistical averaging broadly improves the agreement of heterogeneous results in Part I with the Brost-Wyngaard closure. Recommendations are made for the information needed to analyze apparent differences between homogeneous and heterogeneous stable boundary layers.Part of UK Meteorological Office Atmospheric Processes Research Division     

Red dust rain within the Spanish Mediterranean area

Dust rain belongs to the climatological conditions of the Iberian Peninsula's Mediterranean seaboard. Traditionally known as muddy' or bloody rains, red dust rain has become a topical issue as a result of more frequent rainfall of this kind in recent years. In spite of the difficulties that studying this phenomenon involves, owing to the lack of systematically kept records at meteorological observatories, details are provided of the high frequency of such phenomena within the Spanish Mediterranean area, using records from other phenological sources. The study analyzes the chemical composition of red dust rain and the atmospheric factors giving rise to such phenomena. Whether there is a relationship betwen the higher frequency of dust rain and possible changes in global atmospheric circulation is as yet unclear.     

On the enhanced smoothing over topography in some mesometeorological models

An equation is derived for the components of the horizontal (turbulent) frictional force in the -coordinate system with special attention to mesometeorological flow models. The starting point is the horizontal equation of motion in its flux-form in the -system in which we replace (following Reynolds' procedure) the velocity components u,v and % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbai % aaaaa!37B8! \[ \dot \sigma \] aswell as other relevant quantities by terms of the form u = + u,..., = ± + % MathType!MTEF!2!1!+- % feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbai % Gbauaaaaa!37C3! \[ \dot \sigma ' \] , etc. ( = time average of u; u = fluctuating part of u.) Next, the equation is averaged with respect to time and terms which we believe are small in mesometeorological flows, are neglected. On expressing by an appropriate expression that involves w, the result shows the appearance of two new terms which, have not been considered previously in the published literature. While the expression earlier used in the literature involved the -derivative of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG3bGbauaaaaaaaa!380B!\[\overline {u'w'} \] alone, the new terms add the -derivatives of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaadaahaaWcbeqaaiaaikdaaaaaaaaa!37EC!\[\overline {u'^2 } \] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG2bGbauaaaaaaaa!380A!\[\overline {u'v'} \] for the x-component of the force, and the -derivatives of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG2bGbauaadaahaaWcbeqaaiaaikdaaaaaaaaa!37ED!\[\overline {v'^2 } \]} and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG1bGbauaaceWG2bGbauaaaaaaaa!380A!\[\overline {u'v'} \] for the y-component, where and are the slopes of the -surfaces in the x- and y-directions, respectively. Further, a few numerical simulations of the sea-breeze over topography are carried out with and without the correction terms. It is shown that when corrections terms are not included the effective smoothing is stronger above the sloping regions and may amount to as high as 50 percent of the convergence with slopes of ~.04. The ìnclusìon of the new terms does not lead to any special computational difficulties and for that reason there is no compelling reason to neglect them, all the more so because, as is shown, the addition of the new terms results in a consistent apportioning of the degree of horizontal diffusion.On leave from CIMMS, Norman, OK.Now visiting Dept. of Met., Helsinki, Finland.     

AN EVALUATION OF DISPERSION COEFFICIENTS FOR USE IN AIR QUALITY MODELS

Standard deviations of concentration in horizontal andvertical directions i.e. y andz have been estimated by using fivedifferent schemes based on empirical(due to Pasquill and Briggs)schemes and sophisticated methods(due to Irwin, Draxler, Taylor, Hanna et al.). The fiveschemes are discussed atlength. The purpose of this study is to make use ofmeteorological observations whichare routinely available, to test all the above methods andintercompare the resultswith one another and observations so that the sensitivityof each schemeunder various atmospheric stability conditions could beassessed. It has beenfound that the existing schemes are good enough to providereasonable estimates ofdispersion coefficient (y) during highly unstableconditions (Pasquill stability classes A and B). However, thesame is not true for the case when the stability increasesfrom C to F and turbulencedecreases, specifically during stable atmospheric conditions,when the observedvalues were found to be much higher than all the existingschemes. This suggests thatwhile we continue to use the current methods of estimatingthe dispersion parameters,a rigorous search is required for methods which give predictionswhich are close-to-realityduring such conditions which are represented by lowlevels (in terms of magnitude)of atmospheric turbulence leading to higher levelsof pollution.As one of the sophisticated methods requiresthe use of v and w (standard deviationsof wind velocity fluctuation in y and z directions),these have been estimated andvalidated with observed data (field experiments conductedby EPRI at Kincaid).Statistical evaluation of v and wbased on performance measures indicate a goodperformance of the parameterisations adopted in thisstudy. In the case of w duringunstable conditions a comparison of three differentschemes with observations is made.     

Evidence of recent climatic anomalies at Trieste (Italy)

A recent noticeable climatic anomaly in some meteomarine quantities is analysed based on observations during the period March 1, 1971-February 28, 1993 taken in Trieste (Italy). Marked deviations from the normal values of the coupled quantities atmospheric pressure-solar irradiance and atmospheric pressure-sea level elevation, during the last 22 years, are analysed. The procedures used bring out a clear and persistent configuration of high pressure, high irradiance and low sea level from 1989 on, previously unmatched.     

A numerical study of the Chicago lake breeze

A two-dimensional numerical model is presented for the study of the Chicago area lake breeze. The hydrostatic approximation used by previous workers in this field is not employed. Instead, the complete z-momentum equation is solved and the equation of continuity is retained in its original form. The numerical model employs a staggered mesh, and Alternating Direction Implicit methods are used for the integration of the equations. Computational stability is achieved without external filters, upstream space differencing, or artificial horizontal diffusion terms.The results obtained for realistic initial conditions but simplified boundary conditions compare reasonably well with observations obtained by Lyons in his 1967 observational study of the Chicago lake breeze. Limiting forms of the lower boundary condition bracket the observed wind speeds and other kinematic features of the breeze. Among the various eddy diffusivity models, that with constant eddy diffusivity with height produces the most realistic results. The models with a decrease of eddy diffusivity with height show unrealistic flow patterns in the upper layers. The stability of the atmosphere at the onset of the breeze is shown to have a pronounced effect on the intensity of the circulation, but not on the thickness of the inflow layer, which appears quite insensitive to large changes in the parameters involved.     

Interactions between internal gravity waves and their traumatic effect on a continuous stratification

An investigation is made of the detailed kinematic conditions which characterize breaking of internal waves in a continuous density field. It has been observed that breaking is associated with localized, abruptly appearing intensifications of the density gradient. Turbulence subsequently develops from these intensifications or traumata.Interactions between two or more disparate waves produce distortions more extreme than would be deduced from linear combination. This is exploited experimentally by crossing two weak internal wave beams from separate sources.Lagrangian integration of the motion field defined by the beams gives extrema of isopycnic slope, Richardson number defined by shear across isopycnics and dynamic vertical density gradient, and dynamic density gradient, under conditions for the appearance of traumata.     

Interannual variability of the simulated hydrology in a climatic model-implications for drought

A 10-year simulation of the atmosphere has been performed with a two-level global general circulation model. The model takes account of realistic topography and the hydrological cycle, and computes its own cloud cover, snowfall and sea ice distributions. The seasonally varying sea surface temperatures are specified from climatology. The results presented in this paper are restricted to hydrological aspects of the integration, particularly the soil moisture variations. The model exhibits considerable temporal and regional variability in its simulation of the hydrological cycle, generally reproducing quite fair agreement with observations. Some geographical regions were found to have very noticeable interannual variability, with the occurrence of annual or multi-annual drought a feature at a number of locations. The implications of such naturally occurring drought are discussed.     

Numerical studies of heat island circulations

A two dimensional model has been set up to investigate the circulation induced by an urban heat island in the absence of synoptic winds. The boundary conditions need to be formulated carefully and due to difficulties arising here, we restrict our attention to cases of initially stable thermal stratification. Heat island circulations are allowed to develop from rest and prior to the appearance of the final symmetric double cell pattern, a transitional multi-cell pattern is observed in some cases. The influence on the steady state circulation of various parameters is studied, among which are eddy transfer coefficients, the heat island intensity, the initial temperature stratification and the heat island size. Some results are presented for a case in which differential surface cooling beneath an initially stable atmosphere produces a circulation and an unstable layer capped by an elevated inversion over the city. It is hoped that this case is vaguely representative of the night-time heat island with no geostrophic wind.Notation c_p Specific heat at constant pressure - g Acceleration due to gravity - H Top of integration region - K_z Vertical eddy transfer coefficient - K_x, K_xH, K_xm Horizontal eddy transfer coefficients for heat and momentum - l ixing length - p Pressure - p₀ Reference surface pressure (1000 mb) - P_H (x, t) Pressure at z = H - R Specific gas constant for dry air - t Time - u, w Horizontal and vertical velocities - x, z Horizontal and vertical coordinates - x₁, x₂ Positions of discontinuities in surface temperature field (see Figure 2) - x_a Heat island half-width - x_b Boundary of integration region - Parameter in formula for eddy coefficients (variable-K case) = 18.0 - _s Intensity of heat island - Potential temperature field - Reference absolute temperature (variable-K case) - _r Reference temperature (° C) - _s Surface temperature - Q Air density     

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.     

A note on the analogy between momentum transfer across a rough solid surface and the air-sea interface

The aerodynamic classification of the resistance laws above solid surfaces is based on the use of a so-called Reynolds roughness number Re _s =h _s u _*/, whereh _s is the effective roughness height, -viscosity,u _*-friction velocity. The recent experimental studies reported by Toba and Ebuchi (1991), demonstrated that the observed variability of the sea roughness cannot be explained only on the basis of the classification of aerodynamic conditions of the sea surface proposed by Kitaigorodskii and Volkov (1965) and Kitaigorodskii (1968) even though the latter approach gains some support from recent experimental studies (see for example Geernaertet al. 1986). In this paper, an attempt is made to explain some of the recently observed features of the variability of surface roughness (Toba and Ebuchi, 1991; Donelanet al., 1993). The fluctuating regime of the sea surface roughness is also described. It is shown that the contribution from the dissipation subrange to the variability of the sea surface can be very important and by itself can explain Charnock's (1955) regime.     

Interannual variation of Gulf Stream heat transport in a high-resolution model forced by reanalysis data

The variability present in a 1/6th degree Atlantic ocean simulation forced by analysed wind stress and heat flux over a 20-year period is investigated by means of heat transport diagnostics. A section is defined which follows the Gulf Stream and its seaward extension, and transport of heat across this section is analysed to reveal the physical mechanisms responsible for intergyre heat exchanges on a variety of time scales. Heat transport across another section that crosses the Gulf Stream is also diagnosed to reveal the temporal behaviour of the gyre circulation. The Ekman response to wind stress variations accounts for the annual cycle and much of the interannual variability in both measures. For the intergyre heat transports, cancellation by transient-mean flow terms leads to a very weak annual cycle. Transient eddies account for approximately half the total intergyre transport of 0.7 Petawatts. They also account for a significant fraction of the interannual variability, but separate experiments with repeated-annual-cycle forcing indicate that the transient eddy component of the heat transport variability is internally generated. Links between the intergyre transport, the wind-driven gyre circulation, the surface heat budget and the atmospheric North Atlantic Oscillation are discussed.     

Carbon dioxide emissions in a methane economy

Increasing reliance on natural gas (methane) to meet global energy demands holds implications for atmospheric CO₂ concentrations. Analysis of these implications is presented, based on a logistic substitution model viewing energy technologies like biological species invading an econiche and substituting in case of superiority for existing species. This model suggests gas will become the dominant energy source and remain so for 50 years, peaking near 70 percent of world supply. Two scenarios of energy demand are explored, one holding per capita consumption at current levels, the second raising the global average in the year 2100 to the current U.S. level. In the first (efficiency) scenario concentrations peak about 450 ppm, while in the second (long wave) they near 600 ppm. Although projected CO₂ concentrations in a methane economy are low in relation to other scenarios, the projections confirm that global climate warming is likely to be a major planetary concern throughout the twenty-first century. A second finding is that data on past growth of world per capita energy consumption group neatly into two pulses consistent with longwave theories in economics.     

A generalization of the Kolmogorov-von Kármán relationship and some further implications on the values of the constants

A general relation between the Kolmogorov (₁) and von Kármán (k) constants appropriate to flows in the surface boundary layer is obtained by including stability and divergence dependent terms in the development. The result is used to show that the author's earlier inference of the valuek = 0.36 drawn from the observed value ,k ^4/3=0.141 cannot have been underestimated due to the neglect of stability or divergence effects. It is concluded that if k is to be 0.4, then ₁ must be 0.5.Work performed under the auspices of the U.S. Energy Research and Development Administration.     

The Effect of Atmospheric Stability on the Surface-Layer Characteristics in a Low-Wind Area of Tropical West Africa

This study details the observed effects ofatmospheric stability on characteristics of thesurface layer in a low wind speed (U = 1.5 m s-1)regime of tropical West Africa. Theaerodynamic roughness length, z0, anddisplacement height, d, obtained from profilewind-speed data at our bush land site (height 2 m)have values of 0.24 ± 0.10 m and 1.54 ± 0.04 mrespectively. In the unstable range (-2.5 < Ri < -0.1; Riis gradient Richardson number), thestandard deviation in wind speed fluctuations, u, increased from 0.57 ± 0.19 m s-1 toa maximum of 0.7 ± 0.2 m s-1 in near-neutralconditions, and in the stable range, the parameterdecreased rapidly to 0.41 ± 0.15 m s-1 at Ri 0.2.In the same stability range, the horizontal winddispersion, , decreased withincreasing stability from 19 ± 8 deg. to 13 ± 5 deg.The surface-layer integral quantity, u/u*, when plottedas a function of stability, is in agreement with theempirical results. The ratio ofsensible heat flux (estimated) to the net radiationranged between 0.1 and 0.2 at nighttime,increasing to about 0.5 during the daytime, and showeda strong dependency on season.     

Stable boundary layers: Observations,models and variability part I: Modelling and measurements

Selected field measurements of evening stable boundary layers are presented in detail comparable with published Large Eddy Simulation results. Such models appear to match idealized theories more closely than do some boundary-layer observations. Any attempt to compare detailed observations with idealized models therefore highlights the variability of the real boundary layer.Here direct turbulence measurements across the stable boundary layer from a heterogeneous and an ideal site are contrasted. Recommendations are made for the information needed to distinguish heterogeneous and ideal cases.The companion paper (Part II) discusses further the issues of data, analysis in the presence of variability, and the effects of averaging over heterogeneous terrain.Part of UK Meteorological Office Atmospheric Process Research Division.