Recurrence statistics of concentration fluctuations in plumes within a near-neutral atmospheric surface layer

This study examines the statistical properties of the concentration derivative, , for a dispersing plume in a near-neutrally stratified atmospheric surface layer. Towards this goal, the probability density function (pdf) of , and the conditional pdf of given a fixed concentration level, , have been measured. These pdfs are found to be modeled well by a generalizedq-Gaussian (gqG) distribution with intermittency exponent,q, equal to 0.3 and 3/4, respectively. These results highlight the strong intermittency effect (patchiness) of the small-scale concentration eddy structures in the plume. The distribution of time intervals between successive high peaks in the squared derivative process, x², is found to be well approximated by a power-law distribution, implying that occurrences of these high peaks are much more clustered than would be predicted by a Poisson or shot-noise process. The results are used to improve models for the joint pdf of and , and for the expected number of upcrossings per unit time interval of a fixed concentration level that have been proposed by Kristensenet al. (1989). The predictions of the improved models are in accord with observations, and suggest that the intercorrelation between and must be explicitly incorporated if good estimates of the upcrossing intensity are to be obtained.     

Potential refractivity as a similarity variable

The validity of a common radiometeorological application of Monin-Obukhov (M-O) similarity theory to potential refractivity (), which is a nonlinear combination of and q, is determined by whether the properly nondimensionalized gradient is a universal function of z/L. We develop expressions for the flux of (and its scaling parameter, *) in terms of temperature and moisture fluxes, and an M-O similarity expression for the vertical gradient. Results show that even if and q are accepted as exactly following M-O similarity expressions, when the surface layer is stable, does not obey such an expression. That is, when properly nondimensionalized, the vertical gradient of does not collapse to a single universal function of z/L. The assumption that behaves as a similarity variable is approximately correct for well-mixed surface layers under unstable and near-neutral conditions.The gradient of is an important factor in determining microwave propagation conditions. We demonstrate the error induced in a simple algorithm when is assumed to obey M-O similarity theory. An alternative methodology, consistent with the application of similarity theory to and q, is then developed without requiring that itself satisfy similarity theory.     

Final results of the CONDORS convective diffusion experiment

TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF ₆ values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, _y and _z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations _vand _w 0.6w^* (the convective scale velocity), as measured atz>100m. Surface release _y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger _v and reduced wind speed measured near the surface. Decreased _y at small dimensionless average times is also noted. Finally, convectively scaled dy, C _y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC _y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.     

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.     

Principal oscillation pattern analysis of the tropical 30- to 60-day oscillation

In Part I (Storch and Xu 1990) the principal oscillation pattern (POP) analysis of 200 mb equatorial velocity potential leads to the definition of a bivariate (POP-) index of the tropical 30- to 60-day oscillation. Using the POP prediction scheme this index is predictable for a few days in advance. In Part 11, the prediction of the equatorial velocity potential field, made by the POP method and made by two GCMs, is investigated. The POP index forecast can incorporate skillful forecasts of the equatorial velocity potential () field. Its ensemble correlation skill score passes the 0.50 level at 7 days, whereas persistence passes after 3 days. If there is a strong 30- to 60-day oscillation signal in the initial state, useful forecasts of more than 20 days are sometimes possible; if the initial signal is weak, the POP forecast fails. Also, the forecast skill of two GCMs is considered. The NCAR T31 GCM appears to be quite skillful in predicting the equatorial -field, and in particular the 30- to 60-day oscillation. Its skill, however, is less than that of the POP scheme. The CNRM T42 GCM seems not to be able to predict the regular development associated with the tropical 30- to 60-day oscillation. The power of the POP index in explaining the equatorial x-field is a measure of the strength and dominance of the 30- to 60-day oscillation. This measure at day 0 is an a priori indicator of the NCAR T31 GCM's skill in predicting the equatorial velocity potential field.The National Center for Atmospheric Research is sponsored by the National Science Foundation     

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.     

An approximate analytical solution to the diffusion equation for short-range dispersion from a continuous ground-level source

An easily-evaluated expression for the dimensionless concentration profile (z/z ₀,/z ₀, z ₀/L) = = cu _*/kQ (or z ₀cu_*/kQ) downwind of a continuous ground-level area (or line) source in the stable surface layer is obtained by integrating the diffusion equation using the Shwetz approximation method (c = concentration, Q = source strength, k = von Kárman's constant). The analytical solution compares closely with concentration profiles obtained using a trajectory-simulation model over a useful range of heights, the important discrepancies occurring at the upper edge of the plume. The analytical solution is used to generate predictions of ground-level concentration for the Project Prairie Grass experiments; good agreement with the observations is obtained at all downwind distances (50 to 800 m).     

A multifractal analysis of lidar measured water vapour

Accurate and fast-response measurements of space-time observations of specific humidity were made above a drying land surface at the University of California at Davis, using the Los Alamos water Raman-lidar. In an attempt to quantify the space-time intermittency features of turbulent flows in the lower atmosphere, a multifractal analysis of these water vapour measurements was performed. The structure of the specific humidity, (x, t), was analyzed quantifying a scalar gradient measure both in time and space, for all possible one-dimensional cuts, i.e. and . The results confirm the multifractal nature of this scalar gradient measure (a type of scalar dissipation rate) and show that humidity measurements at fixed times (x) are more intermittent (e.g. have less entropy dimension) than those at fixed locations in space (t). Similar multifractal behaviour of the spatial data, with and without a transformation from the observed wind velocities, supports the validity of Taylor's hypothesis for the studied fields.     

The budget of carbon dioxide variance in the surface layer over vegetated fields

The budget equation for carbon dioxide variance can be represented by production, dissipation and flux divergence terms. Each term is measured under near neutral to moderately unstable conditions over vegetated fields. The flux divergence term is about an order of magnitude smaller than production and dissipation terms, though it shows a loss for 0.006 < _v < 1 and a gain for 1 < - _v < 10. Here, _v is the Monin-Obukhov stability parameter including humidity effect. As expected from a closure of the budget, the nondimensional production and dissipation terms are basically identical and represented by the same functional form: (1–16 _v )^–1/2.     

Similarity forms for ground-source surface-layer diffusion

Past work on analyzing ground-source diffusion data in terms of surface-layer similarity theory is reviewed; these analyses assume that _z /L orh/L is a function of u ^* x/L (where h = Q/ dy). It is argued that an alternative scaling, h ^*/L versus x/L, is nearly as universal in that it is very weakly influenced by surface roughness, except for a modest influence in the free convective case (h ^* = Q/u ^* dy); the advantage of this scaling is that it eliminates the need to reassess as vertical diffusion progresses. The Prairie Grass data set is adjusted for the difference in source and sampling heights, and is plotted with this scaling. Simple analytic equations are suggested that fit the resultant data plots for stable and unstable conditions, and suggestions are made towards practical application of these results.On assignment from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce.     

Estimates of diabatic wind speed profiles from near-surface weather observations

In this paper we analyse diabatic wind profiles observed at the 213 m meteorological tower at Cabauw, the Netherlands. It is shown that the wind speed profiles agree with the well-known similarity functions of the atmospheric surface layer, when we substitute an effective roughness length. For very unstable conditions, the agreement is good up to at least 200 m or z/L–7(z is height, L is Obukhov length scale). For stable conditions, the agreement is good up to z/L1. For stronger stability, a semi-empirical extension is given of the log-linear profile, which gives acceptable estimates up to ~ 100 m. A scheme is used for the derivation of the Obukhov length scale from single wind speed, total cloud cover and air temperature. With the latter scheme and the similarity functions, wind speed profiles can be estimated from near-surface weather data only. The results for wind speed depend on height and stability. Up to 80 m, the rms difference with observations is on average 1.1 m s^–1. At 200 m, 0.8 m s^–1 for very unstable conditions increasing to 2.1 m s^–1 for very stable conditions. The proposed methods simulate the diurnal variation of the 80 m wind speed very well. Also the simulated frequency distribution of the 80 m wind speed agrees well with the observed one. It is concluded that the proposed methods are applicable up to at least 100 m in generally level terrain.     

Size-Resolved Characterisation of Soluble Ions in the Particles in the Tropospheric Plume of Masaya Volcano, Nicaragua: Origins and Plume Processing

We present the first application of a multi-stage impactor to study volcanic particle emissions to the troposphere from Masaya volcano, Nicaragua. Concentrations of soluble SO₄ ^2–,Cl^–, F^–, NO₃ ^–, K⁺, Na⁺,NH₄ ⁺, Ca²⁺ and Mg²⁺ were determined in 11 size bins from 0.07 m to >25.5 m. The near-source size distributions showed major modes at 0.5m (SO₄ ^2–, H⁺,NH₄ ⁺); 0.2 m and 5.0 m (Cl^–) and 2.0–5.0 m(F^–). K⁺ and Na⁺ mirrored the SO₄ ^2– size-resolvedconcentrations closely, suggesting that these were transported primarily asK₂SO₄ and Na₂SO₄ in acidic solution, while Mg²⁺ andCa²⁺ presented modes in both <1 m and >1 m particles. Changes in relative humidity were studied by comparing daytime (transparent plume) and night-time (condensed plume) results. Enhanced particle growth rates were observed in the night-time plume as well as preferential scavenging of soluble gases, such as HCl, by condensed water. Neutralisation of the acidic aerosol by background ammonia was observed at the crater rim and to a greater extent approximately 15 km downwind of the active crater. We report measurements of re-suspended near-source volcanic dust, which may form a component of the plume downwind. Elevated levels ofSO₄ ^2–, Cl^–, F^–,H⁺, Na⁺, K⁺ and Mg²⁺ were observed around the 10 m particle diameter in this dust. The volcanic SO₄ ^2– flux leaving the craterwas 0.07 kg s^–1.     

Increasing cloud in a warming world

Cloud amount records for the U.S.A. have been analyzed in the context of the warming world analogue model described by Lough et al. (1983). Cloud amount increases over practically the entire U.S.A. in all seasons. This result considerably strengthens the more tentative conclusion of Henderson-Sellers (1986) that cloud amount increases over Europe in the same warming world scenario. These results are in contrast to the few numerical model predictions of cloud changes in warming world experiments. A possible, rather tantalizing, conclusion is that current GCM cloud prediction schemes tend to enhance temperature increases through cloud-climate feedback whereas the historical data could suggest a negative feedback. Part, possibly all, of this difference may be the result of the fundamental distinction between the two experimental scenarios: the equilibrium change modelled by GCMs as compared to the smaller transient change represented by the historical analogue. On the other hand the current real-world experiment is, itself, a transient change in boundary and atmospheric conditions. At the least, surface-observed cloudiness seems to offer a useful and complementary data source with which to examine one aspect of the performance of numerical climate models.     

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.     

On the log-linear wind profile and the relationship between shear stress and stability characteristics over the sea

Wind and stability characteristics in the atmospheric surface boundary layer at a height,Z, less than 20 m above the sea were examined in nine oceanic investigations. The analysis lends further support to the utility of the log-linear wind-profile law in the stability region of –0.4Z/L0.9, whereL is the Monin-Obukhov length. However, it is also shown that, inasmuch as better than 90% of the measurements fall within the range of ¦Z/L¦ 0.25, and inasmuch as this correction to the drag coefficient under neutral conditions amounts to less than 10%, the familiar logarithmic wind law may be used rather than the log-linear form. A wind-stress drag coefficient,C _d (=1.2×10^–3 between 1.0 m Z 18.3 m), is thus recommended for general deepwater oceanic applications. The situation over shallow water, which is different, is discussed briefly.     

Turbulent dispersion in a stable layer with a quadratic exchange coefficient

Dispersion from a line source into a stable boundary layer of thickness l is analysed through solution of the diffusion equation assuming an exchange coefficient K(z) (1 – z/l) ² and wind profiles u(z) z ⁿ, with n = 0,1. Estimates of ground-level concentrations are made by developing analytic formulae where this is possible. A general method of solution using Laplace transformation and Green's function techniques is developed as an alternative to the eigenfunction expansion method discussed previously (Robson, 1983).     

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     

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.     

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.     

On the temperature coefficient of PCO 2 in seawater

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