On the similarity in atmospheric fluctuations of carbon dioxide,water vapor and temperature over vegetated fields

This paper describes the similarity between atmospheric fluctuations of carbon dioxide, water vapor and temperature using data which cover a wide range of instability (0.02 < < 10). The is the Monin-Obukhov stability parameter including the humidity effect.The spectral analysis shows that the coherency between fluctuations of carbon dioxide and water vapor or temperature is very close to unity, and the phase difference is basically out of phase for whole frequency ranges analyzed. The stability dependence of the normalized standard deviation of carbon dioxide is very similar to those of water vapor and temperature. The normalized standard deviation is about 2.5 under near neutral conditions, and it decreases with increasing instability following the -1/3; power law as (-)^-1/3. The skewness factors of carbon dioxide, water vapor and temperature show a systematic departure with increasing instabilities for 0.02 < s- < 1, and level off at high instabilities for 1 < -\s < 10. The stability dependence of the flatness factors is not so clear as that noted in the standrard deviation and skewness factors. Dissipation rates of carbon dioxide, water vapor and temperature variance are well related to the spectral peak wavelength. This seems to be real since the local production and local dissipation rates are the main terms, almost balancing one another in the variance budget equations for scalar entities.     

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.     

Properties of spatial temperature derivatives in the atmospheric surface layer

A pair of parallel cold wires separated in either the vertical or lateral direction was used to obtain the three components _x, _y, _z of the temperature derivative in the streamwise, lateral and vertical directions, respectively. The average absolute skewness values of _x and _z are nonzero and approximately equal, while the skewness of _y is approximately zero. These results appear to be consistent with the presence of a large, three-dimensional organised structure in the surface layer. There is an apparent low-frequency contamination in the spectral density of _y and _z due mainly to small errors in estimating the sensitivity of the cold wires. The temperature derivatives were high-pass filtered, the filter being set to remove possible contributions from the large structure and to minimise low-frequency sensitivity contamination. The filtered rms ratios \~_x/\~_y and \~_x/\~_z were in the range 0.7 to 0.9, a result in qualitative agreement with that obtained in the laboratory boundary layer by Sreenivasan et al. (1977). The skewness of filtered _x or _z is negligible, consistent with local isotropy of small-scale temperature fluctuations and in support of the high wavenumber spectral isotropy discussed in Antonia and Chambers (1978).     

The Kolmogorov constant for carbon dioxide in the atmospheric surface layer over a paddy field

From measurements in the atmospheric surface layer over a paddy field, the Kolmogorov constants for CO₂ and longitudinal wind velocity were obtained. In this study, the nondimensional dissipation rate _nc = (1–16 _v )^-1/2 for CO₂ variance and = (1–16 _v )^-1/4 – _v for turbulent energy were used, assuming the equality of the local production term and the local dissipation term, and neglecting the divergence flux term in the budget equation. The value of the constant for CO₂ was consistent with recent determinations for temperature and humidity. The constant for longitudinal wind velocity showed good agreement with other recent observations.     

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.     

Convective Profile Constants Revisited

This paper examines the interpolation betweenBusinger–Dyer (Kansas-type) formulae,_u = (1 -1 6 )^-1/4 and_t = (1 - 16 )^-1/2, and free convection forms. Based on matching constraints, the constants, a_u and a_t, in the convective flux-gradient relations, _u = (1 - a_u )^-1/3 and _t = (1 - a_t )^-1/3, are determined. It isshown that a_u and a_t cannot be completely independent if convective forms are blended with theKansas formulae. In other words, these relationships already carryinformation about a_u and a_t. This follows because the Kansas relations cover a wide stability range (up to = - 2), which includes a lower part of the convective sublayer (about 0.1 < - < 2). Thus, there is a subrange where both Kansas and convective formulae are valid. Matching Kansas formulae and free convection relations within thesubrange 0.1 < - < 2 and independently smoothing ofthe blending function are used to determine a_u and a_t. The values a_u = 10 for velocity and a_t = 34for scalars (temperature and humidity) give a good fit. This new approacheliminates the need for additional independent model constants and yields a`smooth' blending between Kansas and free-convection profileforms in the COARE bulk algorithm.     

Circulation Variability Reflected in Ice Core and Lake Records of the Southern Tropical Andes

The circulation mechanisms of climate anomalies in the southern tropical Andes are of particular interest for the January–February core of the precipitation season. With this focus, we evaluate in context upper-air and surface analyses, water level measurements of Lake Titicaca, and records of net balance and ¹⁸O from ice cores. Precipitation is more abundant with enhanced and southward expanded easterlies through a deep layer of the troposphere over the southern tropical Andes. Concomitant with this is a southward displaced circulation system over the equatorial Atlantic, entailing reduced interhemispheric gradient of sea surface temperature (SST; cold/warm anomalies in the North/South), more southerly position of the surface wind confluence and Intertropical Convergence Zone, and thus more abundant rainfall in Northeast Brazil. Such ensemble of circulation departures in boreal winter is common to the high phase of the Southern Oscillation.¹⁸O in the ice cores from Peru's Quelccaya Icecap, as wellas the cores from Sajama and Ilimani in Bolivia is more negative with more abundant precipitation, both in the same annual cycle and on interannual timescales. The large-scale circulation departures associated with the more negative ¹⁸O are in the sense as for anomalously abundant precipitation activity over the southern tropical Andes. The variability of ¹⁸O seasonally and interannually appears to be controlled mainly by the fate of the water vapor along its trajectory and over the Andes, rather than by the SST of the South Atlantic source region.     

Characteristics of the Lyman-alpha humidiometer and measurements of turbulent humidity fluctuations over the ocean

The characteristics of a Lyman-alpha humidiometer have been carefully examined in an air-conditioned test chamber. The results confirm that when carefully used, this humidiometer is suitable for measurements of turbulent humidity fluctuations. Measurements with a Lyman-alpha humidiometer were carried out in the surface boundary layer over the ocean. The relation between turbulent intensity ( _a = a ^ov2) and the friction humidity (a _*) can be expressed as _a = l.6a _*. The spectrum of turbulent humidity for wind speeds larger than 3 m s ^–1 conforms to the similarity law in the surface boundary layer. The spectrum has two characteristic normalized frequencies, namely, a higher peak and a secondary peak (or a shoulder).     

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.     

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     

Measurements of fluctuations of vertical wind velocity,temperature and radio-refractive index above the sea

Mean spectra of vertical wind velocity, temperature, and humidity and co-spectra of vertical turbulent heat and moisture fluxes are reported, normalized in terms of the similarity theory of Monin-Obukhov. The measurements were made in April 1967 at light-house Alte Weser at a height of 30 m above the German Bight. Ten spectra measured under conditions of moderate instability (z/L –0.08) have been used for the evaluation of a mean spectral curve. The humidity fluctuations have been determined from simultaneous turbulence measurements of temperature and radiorefractive-index, which were obtained by a microwave refractometer.Similar to the BOMEX measurements discussed by Phelps and Pond (1971) our results also show a dissimilarity between the fluctuations of temperature and humidity. Accordingly, in the range of lower and higher natural frequencies, significant differences were present between the co-spectra of the vertical heat and moisture fluxes, although the mean peak frequencies of both the fluxes nearly correspond to that of the mean vertical wind-velocity spectrum (f _m 0.35).     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.     

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.     

Effect of land management on ecosystem carbon fluxes at a subalpine grassland site in the Swiss Alps

Summary The influence of agricultural management on the CO₂ budget of a typical subalpine grassland was investigated at the Swiss CARBOMONT site at Rigi-Seebodenalp (1025m a.s.l.) in Central Switzerland. Eddy covariance flux measurements obtained during the first growing season from the mid of spring until the first snow fall (17 Mai to 25 September 2002) are reported. With respect to the 10-year average 1992–2001, we found that this growing season had started 10 days earlier than normal, but was close to average temperature with above-normal precipitation (100–255% depending on month). Using a footprint model we found that a simple approach using wind direction sectors was adequate to classify our CO₂ fluxes as being controlled by either meadow or pasture. Two significantly different light response curves could be determined: one for periods with external interventions (grass cutting, cattle grazing) and the other for periods without external interventions. Other than this, meadow and pasture were similar, with a net carbon gain of –128±17g Cm^–2 on the undisturbed meadow, and a net carbon loss of 79±17g Cm^–2 on the managed meadow, and 270±24g Cm^–2 on the pasture during 131 days of the growing season, respectively. The grass cut in June reduced the gross CO₂ uptake of the meadow by 50±2% until regrowth of the vegetation. Cattle grazing reduced gross uptake over the whole vegetation period (37±2%), but left respiration at a similar level as observed in the meadow.     

Occurrence of roll circulations in a shallow boundary layer

Meteorological measurements taken at the Näsudden wind turbine site during slightly unstable conditions have been analyzed. The height of the convective boundary layer (CBL) was rather low, varying between 60 and 300 m. Turbulence statistics near the ground followed Monin-Obukhov similarity, whereas the remaining part of the boundary layer can be regarded as a near neutral upper layer. In 55% of the runs, horizontal roll vortices were found. Those were the most unstable runs, with -z _i/L > 5. Spectra and co-spectra are used to identify the structures. Three roll indicators were identified: (i) a low frequency peak in the spectrum of the lateral component at low level; (ii) a corresponding increase in the vertical component at mid-CBL; (iii) a positive covariance {ovvw} together with positive wind shear in the lateral direction (V/z) in the CBL. By applying these indicators, it is possible to show that horizontal roll circulations are likely to be a common phenomenon over the Baltic during late summer and early winter.     

Estimating product yields of carbon-containing products from the atmospheric photooxidation of ambient air alkenes

The yields of products have been calculated for the reactions of hydroxyl radicals and ozone with 19 of the two-through-six carbon anthropogenic alkenes. Based on their rate of reaction, mechanisms of reactions and the ambient air distribution for these alkenes their seasonal ambient air yields have been estimated.Aldehydes predominate as products irrespective of season, with smaller yields of several ketones. Other minor products include carboxylic acids, carbon monoxide, carbon dioxide, and alkenes. About a two-fold increase is estimated in the yields of hot biradicals and their products from summer to winter.One sensitivity analysis was made by recomputing yields at a different OH radical to O₃ concentration than assumed most likely in the calculations discussed above. In addition, the sensitivity of product yields to an estimated range of seasonally averaged sunset-to-sunrise NO₃ radical concentrations was calculated. The effects of free radical reactions are discussed, but these are believed to make a relatively minor contribution within the NO _x -rich atmospheres that contain anthropogenic alkenes.The uncertainties in product yields associated with the range of NO₃ radical concentrations assumed present is relatively small for aldehydes, as is the decrease in yield of the one carbon hot biradical. Larger uncertainties occur for ketones. Significant decreases in yields occur for larger hot biradicals, especially the branched-chain hot radicals in the presence of NO₃ radicals.     

Numerical simulation of particle trajectories in inhomogeneous turbulence,II: Systems with variable turbulent velocity scale

It is shown that for the purpose of trajectory simulation, the vertical velocityw _L (t) of a fluid element, which is moving in a system (such as a forest canopy, or the unstably stratified atmospheric surface layer) whose turbulent velocity scale _w is height-dependent, must be chosen from a frequency-distribution which is asymmetric aboutw _L = 0. If the gradient _w /z varies only slowly with height, correct trajectories may be obtained by adding a bias (where _L is the length scale) to a fluctuating velocity chosen from a symmetric distribution with variance _w ²(z).     

The atmospheric tide as a continuous spectrum: Lunar semidiurnal tide in surface pressure

Summary The standard equations for the theory of atmospheric tides are solved here by an integral representation on the continuous spectrum of free oscillations. The model profile of back-ground temperature is that of the U.S. Standard Atmosphere in the lower and middle atmosphere, and in the lower thermosphere, above which an isothermal top extends to arbitrarily great heights. The top is warm enough to bring both the Lamb and the Pekeris modes into the continuous spectrum.Computations are made for semidiurnal lunar tidal pressure at sea level at the equator, and the contributions are partitioned according to vertical as well as horizontal structure. Almost all the response is taken up by the Lamb and Pekeris modes of the slowest westward-propagating gravity wave. At sea level, the Lamb-mode response is direct and is relatively insensitive to details of the temperature profile. The Pekeris mode at sea level has an indirect response-in competition with the Lamb mode-and, as has been known since the time of its discovery, it is quite sensitive to the temperature profile, in particular to stratopause temperature. In the standard atmosphere the Lamb mode contributes about +0.078 mb to tidal surface pressure at the equator and the Pekeris mode about –0.048 mb.The aim of this investigation is to illustrate some consequences of representing the tide in terms of the structures of free oscillations. To simplify that task as much as possible, all modifying influences were omitted, such as background wind and ocean or earth tide. Perhaps the main defect of this paper's implementation of the free-oscillation spectrum is that, in contrast to the conventional expansion in the structures of forced oscillations, it does not include dissipation, either implicity or explicity, and thus does not satisfy causality. Dissipation could be added implicity by means of an impedance condition, for example, which would cause up-going energy flux to exceed downgoing flux at the base of the isothermal top layer. To achieve complete causality, however, the dissipation must be modeled explicity. Nevertheless, since the Lamb and Pekeris modes are strongly trapped in the lower and middle atmosphere, where dissipation is rather weak (except possibly in the surface boundary layer), more realistic modeling is not likely to change the broad features of the present results.Symbols a earth's mean radius; expansion coefficient in (5.3) - b recursion variable in (7.4); proximity to resonance in (9.2) - c sound speed in (2.2); specific heatc _p in (2.2) - f Coriolis parameter 2sin in (2.2) - g standard surface gravity - h equivalent depth - i ; discretization index in (7.3) - j index for horizontal structure - k index for horizontal structure; upward unit vectork in (2.2) - m wave number in longitude - n spherical-harmonic degree; number of grid layers in a model layer - p tidal pressure perturbation; background pressurep ₀ - q heating function (energy per mass per time) - r tidal state vector in (2.1) - s tidal entropy perturbation; background entropys ₀ - t time - u tidal horizontal velocityu - w tidal vertical component of velocity - x excitation vector defined in (2.3); vertical coordinate lnp _*/p ₀ [except in (3.8), where it is lnp /p ₀] - y vertical-structure function in (7.1) - z geopotential height - A constant defined in (6.2) - C spherical-harmonic expansion coefficient in (3.6) - D vertical cross section defined in (5.6) and (5.9) - E eigenstate vector - F vertical-structure function for eigenstate pressure in (3.2) [re-defined with WKB scaling in (7.2)] - G vertical-structure function for eigenstate vertical velocity in (3.2) [re-defined with WKB scaling in (7.2)] - H pressure-scale height - I mode intensity defined in (8.1) - K quadratic form defined in (4.4) - L quadratic form defined in (4.4); horizontal-structure magnification factor defined in (5.11) - M vertical-structure magnification factor defined in (4.6) - P eigenstate pressure in (3.2); tidal pressure in (6.2) - R tidal state vector in (5.1) - S eigenstate entropy in (3.2); spherical surface area, in differential dS - T background molecular-scale (NOAA, 1976) absolute temperatureT ₀ - U eigenstate horizontal velocityU in (3.2); coefficient in (7.3) - V horizontal-structure functionV for eigenstate horizontal velocity in (3.2); recursion variable in (7.3) - W eigenstate vertical velocity in (3.2) - X excitation vector in (5.1) - Y surface spherical harmonic in (3.7) - Z Hough function defined in (3.6) - +dH/dz - (1––)/2 - Kronecker delta; Dirac delta; correction operator in (7.6) - equilibrium tide elevation - (square-root of Hough-function eigenvalue) - ratio of specific gas constant to specific heat for air=2/7 - longitude - - - background density ₀ - eigenstate frequency in (3.1) - proxy for heating functionq =c _P/t - latitude - tide frequency - operator for the limitz - horizontal-structure function for eigenstate pressure in (3.2) - Hough function defined in (6.2) - earth's rotation speed - horizontal gradient operator - ()₀ background variable - ()_* surface value of background variable - () value at base of isothermal top layer - Õ state vector with zerow-component - , energy product defined in (2.4) - | | energy norm - ()^* complex conjugate With 10 Figures     

Free convection frequency spectra of atmospheric turbulence over the sea

Frequency spectra of atmospheric turbulenceS (f) in the inertial subrange are considered in the free convection regime over the sea surface in a case of motionless instrument measurements (Eulerian frequency spectra). The frequency spectra formulaef _* S (f)/ ² =c (f _*/f)^5/3 for wind velocity (=1–3), temperature (=t) and humidity (=e) fluctuations are derived on the basis of similarity theory and the –5/3 law. These relations also can be derived from a consideration of convective large-scale advection of small eddies. The frequency scalef _* = (N ^{1 2}/)^1/2 (H/z ²)^1/3 is the lower bound of the inertial subrange and it is of order 10^–2 Hz.The spectra formulae are compared with direct measurements of atmospheric turbulence from the fixed research tower in the coastal zone of the Black Sea in calm weather. It is shown that these formulae are realized at least over two to three decades of the frequency range (approximately from 10^–2 to 10 Hz) and values of the numerical coefficients are found. The derived formulae can be used for calculations of sensible and latent heat fluxes by measuring the high-frequency range of spectra at a fixed point at low wind speeds when the conventional inertial dissipation method is not applicable.     

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.