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.     

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.     

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     

On the inner-layer scale height of boundary-layer flow over low hills

Jackson and Hunt's (1975) equation for the depth of the inner layer of flows over low hills does not depend on any closure assumption as contrarily supposed in literature. This equation contains a constant which can arbitrarily be specified. It is suggested that this inner-layer constant should be determined from experimental data. A preliminary check with some data from the Askervein experiment suggests that Jackson and Hunt's equation fits these data almost as well as Jensen's equation provided that fitted inner-layer constants are used.     

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.     

Third- and fourth-order mixed moments of turbulent velocity and temperature fluctuations in the atmospheric surface layer

Fourth-order mixed moments of velocity and temperature fluctuations, measured within the atmospheric surface layer, are compared with results obtained by assuming the quasi-Gaussian approximation. Standard deviations of the products uw, u and w(u and w are the longitudinal and vertical velocity fluctuations; is the temperature fluctuation) are in good agreement with those obtained using the quasi-Gaussian assumption. Good agreement is also obtained between measured and Gaussian estimates of fourth-order moments including all three fluctuations u, w, Schwarz inequalities, commonly used in the clipping approximation in turbulence modelling, are found to provide bounds for third-order moments of w, that are too conservative. More reasonable, tighter, bounds for these moments are given by inequalities obtained by Lumley.     

Nitrate in the atmospheric boundary layer of the tropical South Pacific: Implications regarding sources and transport

Weekly bulk aerosol samples collected at Funafuti, Tuvalu (8°30S, 179°12E), American Samoa (14°15S, 170°35W), and Rarotonga (21°15S, 159°45W), from 1983 through most of 1987 have been analyzed for nitrate and other constituents. The mean nitrate concentration is about 0.11 g m^–3 at each of these stations: 0.107±0.011 g m^–3 at Funafuti; 0.116±0.008 at American Samoa; and 0.117±0.010 at Rarotonga. Previous measurements of mineral aerosol and trace metal concentrations at American Samoa are among the lowest ever recorded for the near-surface troposphere and indicate that this region is minimally affected by transport of soil material and pollutants from the continents. Consequently, the nitrate concentration of 0.11 g m^–3 can be regarded as the natural level for the remote marine boundary layer of the tropical South Pacific Ocean. In contrast, over the tropical North Pacific which is significantly impacted by the transport of material from Asia and North America, the mean nitrate concentrations are about three times higher, 0.29 and 0.36 g m^–3 at Midway and Oahu, respectively. The major sources of the nitrate over the tropical South Pacific are still very uncertain. A very significant correlation between the nitrate concentrations at American Samoa and the concentrations of ²¹⁰Pb suggests that transport from continental sources might be important. This continental source could be lightning, which occurs most frequently over the tropical continents. A near-zero correlation with ⁷Be indicates that the stratosphere and upper troposphere are probably not the major sources. A significant biogenic source would be consistent with the higher mean nitrate concentrations, 0.16 to 0.17 g m^–3, found over the equatorial Pacific at Fanning Island (3°55N, 159°20W) and Nauru (0°32S, 166°57E). The lack of correlation between nitrate and nss sulfate at American Samoa does not necessarily preclude an important role for marine biogenic sources.     

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.     

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.     

Trajectory Curvature As A Selection Criterion For valid Lagrangian Stochastic Dispersion Models

Among well-mixed multi-dimensional Lagrangian stochastic (LS) dispersion models, we observe that those in poorest agreement with observations produce spiralling trajectories, with an associated reduction in dispersion. We therefore investigate statistics of increments d ' to the orientation '= arctan(W'/U') of the Lagrangian velocity-fluctuation vector – as a possible means to distinguish the better LS models within the well-mixed class. Zero-spin models, having d' = 0, are found to provide best agreement with observations. It is not clear however, whether imposition of the zero-spin property selects (in conjunction with the well-mixed condition) a unique model.     

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.     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

Eddy energy dissipation rate and puff diffusion during calms

The paper considers a puff diffusion in its inertial stage when particle separation obeys the laws of the inertial subrange and depends only on eddy energy dissipation rate . The can be determined in the surface layer by the turbulent kinetic energy equation. Similarity equations connect with diffusion measure .A simple analytical model has been deduced to estimate pollutants diffusion during calms.     

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.     

Regular and seasonal behaviour of atmospheric radio noise field strength (ARNFS) over low latitude station calcutta

Summary The electromagnetic radiation of cloud discharge known as atmospheric radio noise field strength (ARNFS) shows a gradual fall from a frequency of 9 kHz to 81 kHz as studied over a period of two years at Calcutta, very close to Bay of Bengal. The main characteristic features of ARNFS at Calcutta are that-(i) ARNFS shows that midday median value is smaller than midnight median value in all months, (ii) level of daily minimum is higher in February and monsoon compared to other seasons, (iii) sunrise effect and sunset effect are well correlated with local sunrise and sunset times, (iv) the magnitude of sunrise fade and sunrise fade rate are maximum in April and lowest during winter period, (v) the magnitude of sunset fade is higher in premonsoon and postmonsoon while it is lowest in monsoon, (vi) number of occurrence of both sunrise effect and sunset effect is remark-ably smaller in monsoon. The positions of the sun and of atmospheric sources are jointly the causes of seasonal and diurnal variations. The missing of sunrise effect and sunset effect are due to local cloud activity and variation of electron density during geomagnetic storms.With 7 Figures     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

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.     

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     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

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.