Time series analyses of concentration and wind fluctuations

Analyses of concentration fluctuation (C) spectra from boundary-layer smoke plume experiments at six separate locations show that the spectra from these experiments generally exhibit an inertial subrange at high frequencies with a slope of -5/3 and indicate peak energy at a time period of about 50 to 100 s. These periods of peak energy are a factor of two to five less than those for the peak of the wind speed fluctuation (u or v) spectra. A general spectral formula fits normalized spectra from the U.S. and Australia, where the frequency, n, is made dimensionless by multiplying by the plume dispersion parameter, _y , and dividing by the wind speed, u. Peak energy occurs at a dimensionless frequency of n _y/u equal to about 0.15. The Kolmogorov constant in the inertial subrange is estimated from a set of averaged spectra. Cross-spectra indicate little relation between concentration and wind fluctuations. However, most of the correlation that exists is due to periods larger than about 10 or 20 s.     

The Interference of Higher-Order Statistics of the Concentration Field Produced by Two Point Sources According to a Generalized Fluctuating Plume Model

The higher-order correlation functions for the concentrationfluctuations arising from a two-point-source configuration have beencalculated analytically within the context of the phenomenology of afluctuating plume model (viz., a meandering plume model that explicitlyincorporates internal fluctuations). Explicit expressions for thesecond-, third-, and fourth-order correlationfunctions between the concentrationfluctuations produced by two point sources are given in terms of the sourceseparation d and the five physically based parameters that define thegeneralized fluctuating plume model: namely, the absolute plume dispersion,_a, which determines the outer plume length scale; the relative plume dispersion, _r, which determines the inner plume length scale; the fluctuation intensity, i_r, in relative coordinates, which determines the internal concentration fluctuation level; the correlation coefficient, r,between the positions of the centroids of the two interfering plumes; and,the correlation coefficient, r^*, between the concentration fluctuationsof the two plumes in relative coordinates, which determines the degree ofinternal mixing of the two scalars. Furthermore, the form of the totalconcentration probability density function arising from the interferenceproduced by two point sources is presented. Predictions for the second-ordercorrelation function, , and for the total concentration probabilitydensity function have been compared with some new experimental data fora two-point-source configuration in grid turbulence generated in awater-channel simulation. These results are in good agreement with the dataand suggest that the analytical model for the second-order correlationfunction and the total concentration probability density function canreproduce many qualitative trends in the interaction of plumes from twosources.     

Microwave vegetation optical depth and inverse modelling of soil emissivity using Nimbus/SMMR satellite observations

Summary A radiative transfer model has been used to determine the large scale effective 6.6 GHz and 37 GHz optical depths of the vegetation cover. Knowledge of the vegetation optical depth is important for satellite-based large scale soil moisture monitoring using microwave radiometry. The study is based on actual observed large scale surface soil moisture data and observed dual polarization 6.6 and 37 GHz Nimbus/SMMR brightness temperatures over a 3-year period. The derived optical depths have been compared with microwave polarization differences and polarization ratios in both frequencies and with Normalized Difference Vegetation Index (NDVI) values from NOAA/AVHRR. A synergistic approach to derive surface soil emissivity from satellite observed brightness temperatures by inverse modelling is described. This approach improves the relationship between satellite derived surface emissivity and large scale top soil moisture fromR ²=0.45 (no correction for vegetation) toR ²=0.72 (after correction for vegetation). This study also confirms the relationship between the microwave-based MPDI and NDVI earlier described and explained in the literature.List of Symbols f frequency [Hz] - f _i(p) fractional absorption at polarizationp - h surface roughness - h h cos² - H horizontal polarization - n _i complex index of refraction - p polarization (H orV) - R _s microwave surface reflectivity - T _B(p) brightness temperature at polarizationp - T ^* normalized brightness temperature - T polarization difference (T _v-T _H) - T _s temperature of soil surface - T _c temperature of canopy - T _max daily maximum air temperature - T _min daily minimum air temperature - V vertical polarization - soil moisture distribution factor; also used for the constant to partition the influence of bound and free water components to the dielectric constant of the mixture - empirical complex constant related to soil texture - microwave transmissivity of vegetation (=e ^–) - ^* effective transmissivity of vegetation (assuming =0) - microwave emissivity - _s emissivity of smooth soil surface - _rs emissivity of rough soil surface - _vs emissivity of vegetated surface - soil moisture content (% vol.) - K dielectric constant [F·m^–1] - K _fw dielectric constant of free water [F·m^–1] - K _ss dielectric constant of soil solids [F·m^–1] - K _m dielectric constant of mixture [F·m^–1] - K _o permittivity of free space [8.854·10^–12 F·m^–1] - high frequency limit ofK _wf [F·m^–1] - wavelength [m] - incidence angle [degrees from nadir] - polarization ratio (T _H/T _V) - b soil bulk density [gr·cm^–3] - s soil particle density [gr·cm^–3] - R surface reflectivity in red portion of spectrum - NIR surface reflectivity in near infrared portion of spectrum - _eff effective conductivity of soil extract [mS·cm^–1] - vegetation optical depth - _6.6 vegetation optical depth at 6.6 GHz - ₃₇ vegetation optical depth at 37 GHz - ^* effective vegetation optical depth (assuming =0) - single scattering albedo of vegetation With 12 Figures     

Simulating Concentration Fluctuation Time Series with Intermittent Zero Periods and Level Dependent Derivatives

Intermittent concentration fluctuation time series were produced with a stochastic numerical model derived from the assumption that the concentration fluctuations at a fixed receptor in a point-source plume can be modelled as a first order Markov process. The time derivative of concentration was assumed to be level-dependent and constrained by a stationary lognormal probability density function. The input parameters required to reconstruct the intermittent time series are the intermittency factor , the conditional fluctuation intensity i _p ² , and the time scale T _c . A clipped lognormal probability distribution was used to describe the fluctuation time series. Good agreement between the stochastic simulation and experimental water-channel data was demonstrated by comparing the time derivative of concentration and the upcrossing rates over a range of intermittency factors = 0.7 to 0.01 and fluctuation intensities i _w ² = 2.2 to 7.5.     

On characteristics of wind direction fluctuations in the surface layer

It is shown that the ratio of standard deviation of lateral velocity to the friction velocity, /u _*, and therefore wind direction fluctuations, are sensitive to mesoscale terrain properties. Under neutral conditions, /u _* is almost 40% larger in rolling terrain than over a horizontal surface. In the lee of a low mountain, the fluctuations may be 2.5 times as strong as over horizontal terrain. In contrast, vertical velocity fluctuations are little influenced by mesoscale terrain features.Now with Air Weather Service, Offutt AFB, Omaha, Nebraska.     

A statistical-dynamical downscaling procedure for global climate simulations

Summary A statistical-dynamical downscaling procedure for global climate simulations is described. The procedure is based on the assumption that any regional climate is associated with a specific frequency distribution of classified large-scale weather situations. The frequency distributions are derived from multi-year episodes of low resolution global climate simulations. Highly resolved regional distributions of wind and temperature are calculated with a regional model for each class of large-scale weather situation. They are statistically evaluated by weighting them with the according climate-specific frequency. The procedure is exemplarily applied to the Alpine region for a global climate simulation of the present January climate.List of Symbols west-east mesh size in geographic coordinates south-north mesh size in geographic coordinates N number of large-scale weather classes n number of regional-scale event classes p pressure P probability Ø large-scale event regional-scale event q _v specific humidity potential temperature u west-east wind component v south-north wind componentAbbreviations AGL above ground level - LT local time - UTC universal time coordinated With 13 Figures     

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.     

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.     

Temperature distribution and current system in a convectively mixed lake

During spring and autumn, many lakes in temperate latitudes experience intensive convective mixing in the vertical, which leads to almost isothermal conditions with depth. Thus the regime of turbulence appears to be similar with that characteristic of convective boundary layers in the atmosphere. In the present paper a simple analytical approach, based on boundary-layer theory, is applied to convective conditions in lakes. The aims of the paper are firstly to analyze in detail the temperature distribution during these periods, and secondly to investigate the current system, created by the horizontal temperature gradient and wind action. For these purposes, simple analytical solutions for the current velocities are derived under the assumption of depth-constant temperatures. The density-induced current velocities are shown to be small, in the order of a few mm/sec. The analytical model of wind-driven currents is compared with field data. The solution is in good qualitative agreement with observed current velocities under the condition that the wind field is steady for a relatively long time and that residual effects from former wind events are negligible.The effect of the current system on an approximately depth-constant temperature distribution is then checked by using the obtained current velocity fields in the heat transfer equation and deriving an analytical solution for the corrected temperature field. These temperature corrections are shown to be small, which indicates that it is reasonable to describe the temperature distribution with vertical isotherms.Notation T temperature - t time - x, y, z cartesian coordinates - molecular viscosity - _h , _v horizontal and vertical turbulent viscosity - K _h ,K _v horizontal and vertical turbulent conductivity - Q heat flux through the water surface - D depth - u, v, w average current velocity components inx, y andz directions - f Coriolis parameter - p pressure - density - g gravity acceleration - a constant in the freshwater state equation - h _s deviation from the average water surface elevation - L _*,H _* length and depth scale - U _*,W _* horizontal and vertical velocity scale - T temperature difference scale - bottom slope - u _* friction velocity at the water surface - von Karman constant - L Monin-Obukhov length scale - buoyancy parameter - l turbulence length scale - C ₁,C ₂,C ₃ dimensionless constants in the expressions for the vertical turbulent viscosity - , dimensionless vertical coordinate and dimensionless local depth - angle between surface stress direction andx-axis - T _bx ,T _by bottom stress components - c bottom drag coefficient     

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.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

Turbulence mechanisms at an agricultural site

An extensive set of turbulence data from the 3- and 12-m heights taken over an agricultural site (Marsta, Sweden) are analyzed and compared with data from ideal sites.In unstable air, Monin-Obukhov similarity is found to be valid for the non-dimensional gradients of wind, _m , temperature, _h , and humidity, _e , for (only a few data), for _T /|T _*|,/ _E /|E _*| and for the non-dimensionalized inertial subrange spectra of temperature and humidity. Where comparison is possible, the unstable data also agree with those found in the Kansas study, with one remarkable exception, the inertial subrange constant of the temperature spectrum, ₁, being only 0.39, compared to the value 0.80 found at the Kansas site.On the stable side, most similarity predictions break down, with most of the data differing systematically from the corresponding Kansas results, the only exception being . The inertial subrange constants for temperature, ₁, and for humidity, ₁ are found to have the same values, 0.39 and 0.30, respectively, as they do on the unstable side. Remarkable similarity is found for the shape of the stable u- and - and e-spectra. In addition, this shape is found to be identical with that found in Kansas. The peak wavelength of the stable u-, and -spectra is found to be about four times larger than it is for the corresponding Kansas spectra. This is interpreted to be a result of the increased macro-roughness at the Marsta site as compared with that at the Kansas site. A possible explanation for the low ₁-value is discussed, suggesting that ₁ is not a universal constant, but instead dependent on the turbulent structure.     

Bistability of the thermohaline circulation identified through comprehensive 2-parameter sweeps of an efficient climate model

The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (F_a) and atmospheric moisture diffusivity (K_q). In a large ensemble of spinup experiments, the values of F_a and K_q are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of F_a and K_q, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.     

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     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

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.     

Estimation of the Monin-Obukhov similarity functions from a spectral model

From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function _m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - C_i constants i=u, v, w, (see Equation (5) - k_me,k_mT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - k_mi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - k_sb, k_c characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - _m nondimensional mean shear - _H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - _D nondimensional divergence of vertical transports of turbulent linetic energy     

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.     

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     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.