Observational evidence for the Monin-Obukhov similarity under all stability conditions

Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the Monin-Obukhov similarity for describing the turbulent fluctu- ations of 3D winds under all stability conditions and to obtain the turbulence characteristics under different weather conditions. The results show that the dimensionless standard deviations of turbulent velocity com- ponents (σ u /u* , σ v /u* , σ w /u * ) and dimensionless turbulent kinetic energy (TKE) can be well described by "1/3" power law relationships under stable, neutral, and unstable conditions, with σ u /u * > σ v /u * > σ w /u* . Land use and land cover changes mainly impact dimensionless standard deviations of horizontal component fluctuations, but they have very little on those of the vertical component. The dimensionless standard devi- ations of wind components and dimensionless TKE are remarkably affected by different weather conditions; the deviations of horizontal wind component and dimensionless TKE present fog day > clear sky > overcast > cloudy; the trend of the vertical wind component is the reverse. The surface drag coefficient at a Nan- jing suburban measurement site during the observation period was obviously higher than at other reported plains and plateau areas, and was approximately one order larger in magnitude than the reported plains areas. Dimensionless standard deviation of temperature declined with increasing |z /L| with an approximate "-1/3" slope in unstable stratification and "-2/3" slope in stable stratification.     

Organised Motion in a Tall Spruce Canopy: Temporal Scales,Structure Spacing and Terrain Effects

This study investigates the organised motion near the canopy-atmosphere interface of a moderately dense spruce forest in heterogeneous, complex terrain. Wind direction is used to assess differences in topography and surface properties. Observations were obtained at several heights above and within the canopy using sonic anemometers and fast-response gas analysers over the course of several weeks. Analysed variables include the three-dimensional wind vector, the sonic temperature, and the concentration of carbon dioxide. Wavelet analysis was used to extract the organised motion from time series and to derive its temporal scales. Spectral Fourier analysis was deployed to compute power spectra and phase spectra. Profiles of temporal scales of ramp-like coherent structures in the vertical and longitudinal wind components showed a reversed variation with height and were of similar size within the canopy. Temporal scales of scalar fields were comparable to those of the longitudinal wind component suggesting that the lateral scalar transport dominates. The existence of a – 1 power law in the longitudinal power spectra was confirmed for a few cases only, with a majority showing a clear 5/3 decay. The variation of effective scales of organised motion in the longitudinal velocity and temperature were found to vary with atmospheric stability, suggesting that both Kelvin-Helmholtz instabilities and attached eddies dominate the flow with increasing convectional forcing. The canopy mixing-layer analogy was observed to be applicable for ramp-like coherent structures in the vertical wind component for selected wind directions only. Departures from the prediction of m = Λ _w L _s ⁻¹ = 8–10 (where Λ _w is the streamwise spacing of coherent structures in the vertical wind w and L _s is a canopy shear length scale) were caused by smaller shear length scales associated with large-scale changes in the terrain as well as the vertical structure of the canopy. The occurrence of linear gravity waves was related to a rise in local topography and can therefore be referred to as mountain-type gravity waves. Temporal scales of wave motion and ramp-like coherent structures were observed to be comparable.     

Aerosol light absorption, black carbon, and elemental carbon at the Fresno Supersite, California

Particle light absorption (b_ap), black carbon (BC), and elemental carbon (EC) measurements at the Fresno Supersite during the summer of 2005 were compared to examine the equivalency of current techniques, evaluate filter-based b_ap correction methods, and determine the EC mass absorption efficiency (σ_ap) and the spectral dependence of b_ap. The photoacoustic analyzer (PA) was used as a benchmark for in-situ b_ap. Most b_ap measurement techniques were well correlated (r ≥ 0.95). Unadjusted Aethalometer (AE) and Particle Soot Absorption Photometer (PSAP) b_ap were up to seven times higher than PA b_ap at similar wavelengths because of absorption enhancement by backscattering and multiple scattering. Applying published algorithms to correct for these effects reduced the differences to 24 and 17% for the AE and PSAP, respectively, at 532 nm. The Multi-Angle Absorption Photometer (MAAP), which accounts for backscattering effects, overestimated b_ap relative to the PA by 51%. BC concentrations determined by the AE, MAAP, and Sunset Laboratory semi-continuous carbon analyzer were also highly correlated (r ≥ 0.93) but differed by up to 57%. EC measured with the IMPROVE/STN thermal/optical protocols, and the French two-step thermal protocol agreed to within 29%. Absorption efficiencies determined from PA b_ap and EC measured with different analytical protocols averaged 7.9 ± 1.5, 5.4 ± 1.1, and 2.8 ± 0.6 m²/g at 532, 670, and 1047 nm, respectively. The Angström exponent (α) determined from adjusted AE and PA b_ap ranged from 1.19 to 1.46. The largest values of α occurred during the afternoon hours when the organic fraction of total carbon was highest. Significant biases associated with filter-based measurements of b_ap, BC, and EC are method-specific. Correcting for these biases must take into account differences in aerosol concentration, composition, and sources.     

Effects of Topographical Slope Angle and Atmospheric Stratification on Surface-Layer Turbulence

The effect of topographical slope angle and atmospheric stratification on turbulence intensities in the unstably stratified surface layer have been parameterized using observations obtained from a three-dimensional sonic anemometer installed at 8 m height above the ground at the Seoul National University (SNU) campus site in Korea for the years 1999–2001. Winds obtained from the sonic anemometer are analyzed according to the mean wind direction, since the topographical slope angle changes significantly along the azimuthal direction. The effects of the topographical slope angle and atmospheric stratification on surface-layer turbulence intensity are examined with these data. It is found that both the friction velocity and the variance for each component of wind normalized by the mean wind speed decrease with increase of the topographical slope angle, having a maximum decreasing rate at very unstable stratification. The decreasing rate of the normalized friction velocity (u _* /U) is found to be much larger than that of the turbulence intensity of each wind component due to the reduction of wind shear with increase in slope angle under unstable stratification. The decreasing rate of the w component of turbulence intensity (σ _w /U) is the smallest over the downslope surface whereas that of the u component (σ _u /U) has a minimum over the upslope surface. Consequently, σ _w /u _* has a maximum increasing rate with increase in slope angle for the downslope wind, whereas σ _u /u _* has its maximum for the upslope wind. The sloping terrain is found to reduce both the friction velocity and turbulence intensity compared with those on a flat surface. However, the reduction of the friction velocity over the sloping terrain is larger than that of the turbulence intensity, thereby enhancing the turbulence intensity normalized by the friction velocity over sloping terrain compared with that over a flat surface.     

Measurement of the photodissociation coefficient of NO2 in the atmosphere: II,stratospheric measurements

The photodissociation coefficient of NO₂, J _{NO 2}, has been measured from a balloon platform in the stratosphere. Results from two balloon flights are reported. High Sun values of J _{NO 2} measured were 10.5±0.3 and 10.3±0.3×10^-3 s^-1 at 24 and 32 km respectively. The decrease in J _{NO 2} at sunset was monitored in both flights. The measurements are found to be in good agreement with calculations of J _{NO 2} using a simplified isotropic multiple scattering computer routine.     

Variations in atmospheric CO2 at the Meteorological Research Institute, Tsukuba, Japan

Since April 1986, measurements of the CO₂ concentration in the surface air have been conducted at the Meteorological Research Institure (MRI, 36°04 N, 140°07 E, 25 m above sea level) in Tsukuba, located 50 km northeast of Tokyo, Japan. The CO₂ data measured over times between 11:00 Japan Standard Time (JST) and 16:00 JST (C _N ) were considered to be representative of the air (within a few ppmv) in the planetary boundary layer. To evaluate the representative CO₂ level on a spatial scale larger than that of the C _N record, the CO₂ data with hour-to-hour variation less than 1 ppmv were selected (C _P ). Comparison of these data with those of Ryori (39°02 N, 141°50 E), a continental station operated by the Japan Meteorological Agency, indicates that the C _P record provides a representative CO₂ level in the air on spatial scales of at least a few hundred kilometers.The C _N record allows an investigation of the internanual changes in photosynthesis/respiration against changes in climatological parameters. Within a small temperature anomaly (ca.±1 °C) respiration is sensitive to the temperature change, while photosynthesis is less sensitive. When the temperature anomaly is large, however, photosynthesis and respiration tend to be competitive.     

Evaluation of a catalytic reduction technique for the measurement of total reactive odd-nitrogen NO y in the atmosphere

A catalytic reduction technique for the measurement of total reactive odd-nitrogen NO _y in the atmosphere was evaluated in laboratory and field tests. NO _y component species include NO, NO₂, NO₃, HNO₃, N₂O₅, CH₃COO₂NO₂(PAN), and particulate nitrate. The technique utilizes the reduction of the higher oxides to NO in reaction with CO on a metal catalyst and the subsequent detection of NO by chemiluminescence produced in reaction with O₃. The efficiency and linearity of the conversion of the principal NO _y species were examined for mixing ratios in the range of 0.1 to 100 parts per billion by volume (ppbv). Results of tests with Au, Ni, and stainless steel as the catalyst in the temperature range of 25–500°C showed Au to be the preferred catalyst. NH₃, HCN, N₂O, CH₄, and various chlorine and sulfur compounds were checked as possible sources of NO _y interference with the Au catalyst. The effects of pressure, O₃, and H₂O on NO _y conversion were also examined. The results of the checks and tests in the laboratory showed the technique to be suitable for initial NO _y measurements in the atmosphere. The technique was subsequently tested in ambient air at a remote ground-based field site located near Niwot Ridge, Colorado. The results of conversion and inlet tests made in the field and a summary of the NO _y data are included in the discussion.     

Boundary-layer resistances of artificial leaves in turbulent air: I. Leaves parallel to the mean flow

Boundary-layer resistance to heat transfer from plates was studied in a wind tunnel which produced turbulence with streamwise intensity in the range 3.5 to 25% and a longitudinal integral scale of the streamwise turbulence component (L _u,x) in the range from 8 to 100 mm. It was found that heat transfer enhancement occurred due to the turbulence, and that at any given intensity, this enhancement was determined by the ratio of L _u,x to the characteristic dimension of the plate.     

On the Time Evolution of the Turbulent Kinetic Energy Spectrum for Decaying Turbulence in the Convective Boundary Layer

Our focus is the time evolution of the turbulent kinetic energy for decaying turbulence in the convective boundary layer. The theoretical model with buoyancy and inertial transfer terms has been extended by a source term due to mechanical energy and validated against large-eddy simulation data. The mechanical effects in a boundary layer of height z _i at a convective surface-layer height z = 0.05z _i are significant in the time evolution of the vertical component of the spectrum, i.e. they enhance the decay time scale by more than an order of magnitude. Our findings suggest that shear effects seem to feedback to eddies with smaller wavenumbers, preserving the original shape of the spectrum, and preventing the spectrum from shifting towards shorter wavelengths. This occurs in the case where thermal effects only are considered.     

Daily, weekly and seasonal relationships among VOCs, NOx x and O3 in a semi-urban area near Barcelona

Daily, weekly, and seasonal patterns of O₃, NOx _x and VOCs and their relationship to meteorological conditions were studied in a semi-urban site near Barcelona by means of five-day long campaigns that included weekend and labor days in December, March, June, and October. The plant protection thresholds for ozone and NO₂ were exceeded, respectively, on all the studied days in summer and on all the studied days. Ozone formation was predominantly local and relied on photochemical processes with VOCs playing a controlling role. Formaldehyde, acetaldehyde, methanol, toluene, isoprene, and acetone (in this order) presented the highest O₃ formation potential during the studied periods. These results highlight the important role in O₃ formation played by VOC species such as acetaldehyde, methanol, and acetone, that all have a significant biogenic component. Thus, these VOCs must be taken into account in the discussion of any ozone abatement strategy.     

Spectral Short-circuiting and Wake Production within the Canopy Trunk Space of an Alpine Hardwood Forest

Using synchronous multi-level high frequency velocity measurements, the turbulence spectra within the trunk space of an alpine hardwood forest were analysed. The spectral short-circuiting of the energy cascade for each velocity component was well reproduced by a simplified spectral model that retained return-to-isotropy and component-wise work done by turbulence against the drag and wake production. However, the use of an anisotropic drag coefficient was necessary to reproduce these measured component-wise spectra. The degree of anisotropy in the vertical drag was shown to vary with the element Reynolds number. The wake production frequency in the measured spectra was shown to be consistent with the vortex shedding frequency at constant Strouhal number given by f _vs = 0.21ū/d, where d can be related to the stem diameter at breast height (dbh) and ū is the local mean velocity. The energetic scales, determined from the inflection point instability at the canopy–atmosphere interface, appear to persist into the trunk space when , where C _du is the longitudinal drag coefficient, a _cr is the crown-layer leaf area density, h _c is the canopy height, and β is the dimensionless momentum absorption at the canopy top.     

The Surface Energy Budget and Its Impact on the Freeze-thaw Processes of Active Layer in Permafrost Regions of the Qinghai-Tibetan Plateau

The surface energy budget is closely related to freeze-thaw processes and is also a key issue for land surface process research in permafrost regions.In this study,in situ data collected from 2005 to 2015 at the Tanggula site were used to analyze surface energy regimes,the interaction between surface energy budget and freeze-thaw processes.The results confirmed that surface energy flux in the permafrost region of the Qinghai-Tibetan Plateau exhibited obvious seasonal variations.Annual average net radiation(R_n)for 2010 was 86.5 W m^-2,with the largest being in July and smallest in November.Surface soil heat flux(G₀)was positive during warm seasons but negative in cold seasons with annual average value of 2.7 W m^-2.Variations in R_n and G₀ were closely related to freeze-thaw processes.Sensible heat flux(H)was the main energy budget component during cold seasons,whereas latent heat flux(LE)dominated surface energy distribution in warm seasons.Freeze-thaw processes,snow cover,precipitation,and surface conditions were important influence factors for surface energy flux.Albedo was strongly dependent on soil moisture content and ground surface state,increasing significantly when land surface was covered with deep snow,and exhibited negative correlation with surface soil moisture content.Energy variation was significantly related to active layer thaw depth.Soil heat balance coefficient K was>1 during the investigation time period,indicating the permafrost in the Tanggula area tended to degrade.     

Flux between soil and atmosphere, vertical concentration profiles in soil, and turnover of nitric oxide: 2. Experiments with naturally layered soil cores

Intact soils cores were taken with a stainless steel corer from a sandy podzol and a loamy luvisol, and used to measure the flux (J) of NO between soil and atmosphere and the vertical profile of the NO mixing ratios (m) in the soil atmosphere, both as function of the NO mixing ratio (m _a) in the atmosphere of the headspace. These measurements were repeated after stepwise excavation of the soil column from the top, e.g. by removing the upper 2 cm soil layer. The gaseous diffusion coefficients of NO in the soil cores were either computed from soil porosity or were determined from experiments using SF6. The NO fluxes (J) that were actually measured at the soil surface were compared to the fluxes which were calculated either from the vertical NO profiles (J _c ) or from the NO production and uptake rates (J _m ) determined in the excavated soil samples. In the podzol, the actually measured (J) and the calculated (J _m , J_m) NO fluxes agreed within a factor of 2. In the luvisol, the measured NO fluxes (J) and those calculated from the vertical NO profiles (J _c ) also agreed well, but in the upper 6 cm soil layer the NO fluxes (J _m ) calculated from NO production and uptake rates were up to 7 times higher than the measured NO fluxes. This poor agreement was probably due to the inhomogeneous distribution of NO production and consumption processes and the change of diffusivities within the top layers of the luvisol. Indeed, the luvisol showed a pronounced maximum of the NO mixing ratios at about 6 cm depth, whereas the podzol column exhibited a steady and exponential decrease of the NO mixing ratios with depth. The inhomogeneities in the luvisol were confirmed by incubation of the soil cores under anoxic conditions. This treatment resulted in production of NO at several depths indicating a zonation of increased potential activities within the luvisol profile which may have biased the modelling of the NO surface flux from turnover measurements in soil samples. Inhomogeneities could be achieved even in homogenized soil by fertilization with nitrate solution.     

A determination of the CH4, NO x and CO2 emissions from the Prudhoe Bay,Alaska oil development

In this paper we quantify the CH₄, CO₂ and NO _x emissions during routine operations at a major oil and gas production facility, Prudhoe Bay, Alaska, using the concentrations of combustion by products measured at the NOAA-CMDL observatory at Barrow, Alaska and fuel consumption data from Prudhoe Bay. During the 1989 and 1990 measurement campaigns, 10 periods (called events) were unambiguously identified where surface winds carry the Prudhoe Bay emissions to Barrow (approximately 300 km). The events ranged in duration from 8–48 h and bring ambient air masses containing substantially elevated concentrations of CH₄, CO₂ and NO _y to Barrow. Using the slope of the observed CH₄ vs CO₂ concentrations during the events and the CO₂ emissions based on reported fuel consumption data, we calculate annual CH₄ emissions of (24+/–8)×10³ metric tons from the facility. In a similar manner, the annual NO _x emissions are calculated to be (12+/–4)×10³ metric tons, which is in agreement with an independently determined value. The calculated CH₄ emissions represent the amount released during routine operations including leakage. However this quantity would not include CH₄ released during non-routine operations, such as from venting or gas flaring.     

A test of isotropy and Taylor's hypothesis in the atmospheric boundary layer

This study attempts to determine the scales of turbulence in a high Reynolds number shear flow near which transition to isotropy occurs and the scales for which Taylor's hypothesis is applicable. The flow studied was the wind near height x ₃ = 2 m above a flat land surface. Four hot-wire anemometers were mounted in a three-dimensional array with equal separations between 1.8 m and 2 cm in three different directions. Theoretical cross-spectra were computed from the observed spectra of downwind velocity fluctuations assuming isotropy and Taylor's hypothesis. Comparison between these and the observed cross-spectra revealed that the turbulence in the flow studied was consistent with both assumptions provided k _1x3&> 20, where k ₁ is the radian wavenumber; this was the lower bound to which no departure from isotropy could be detected by the experiment. For 4 k _1x3 20, the observations are consistent with symmetry of the turbulence about the downstream direction. That part of Taylor's hypothesis relating observed frequency at a stationary sensor to the downstream wavenumber component appears to be justified within experimental error for k _1x3& > 3.     

Size-resolved aerosol chemical concentrations at rural and urban sites in Central California, USA

Aerosol size distributions were measured with Micro Orifice Uniform Deposit Impactor (MOUDI) cascade impactors at the rural Angiola and urban Fresno Supersites in California's San Joaquin Valley during the California Regional PM₁₀/PM_2.5 Air Quality Study (CRPAQS) winter campaign from December 15, 2000 to February 3, 2001. PM_2.5 filter samples were collected concurrently at both sites with Sequential Filter Samplers (SFS). MOUDI nitrate (NO₃⁻) concentrations reached 66 μg/m³ on January 6, 2001 during the 1000–1600 PST (GMT-8) period. Pair-wise comparisons between PM_2.5 MOUDI and SFS concentrations revealed high correlations at the Angiola site (r > 0.93) but more variability (r < 0.85) at the Fresno site for NO₃⁻, sulfate (SO₄⁼), and ammonium (NH₄⁺). Correlations were higher at Fresno (r > 0.87) than at Angiola (r < 0.7) for organic carbon (OC), elemental carbon (EC), and total carbon (TC). NO₃⁻ and SO₄⁼ size distributions in Fresno were multi-modal and wider than the uni-modal distributions observed at Angiola. Geometric mean diameters (GMD) were smaller for OC and EC than for NO₃⁻ and SO₄⁼ at both sites. OC and EC were more concentrated on the lowest MOUDI stage (0.056 µm) at Angiola than at Fresno. The NO₃⁻ GMD increased from 0.97 to 1.02 µm as the NO₃⁻ concentration at Angiola increased from 43 to 66 µg m^− 3 during a PM_2.5 episode from January 4–7, 2001. There was a direct relationship between GMD and NO₃⁻ and SO₄⁼ concentrations at Angiola but no such relationships for OC or EC. This demonstrates that secondary aerosol formation increases both concentration and particle size for the rural California environment.     

Measurements of the humidity structure function parameters,C q 2and C Tq,over the ocean

Atmospheric turbulence measurements, including temperature and humidity fluctuations, were made from the R/V Acania off the coast of California in June, 1979. The purpose of the experiment was to investigate the scaling properties of the humidity structure function parameter (C _q ²) and temperaturehumidity cospectrum structure parameter (C _Tq) in the marine surface layer. The bulk parameterization method was used to obtain Monin-Obukhov Similarity (MOS) scaling parameters u _*, T _*, q _*and L. Assuming a neutral stability humidity drag coefficient c _qn = 1.3 × 10^-3the dimensionless humidity structure function parameter C _q ²Z^2/3/q_* ²was found to be 18% lower than the corresponding temperature function obtained by Wyngaard et al. (1971). Furthermore, the measurements indicate that the temperature-humidity fluctuations are highly coherent well into the inertial subrange. The results have direct application to turbulent scattering of waves propagating in the atmosphere (particularly microwaves) and methods of estimating air-sea surface fluxes.     

Vertical profiles of the structure parameter of temperature in the stable,nocturnal boundary layer

Vertical profiles of the structure parameter of temperature C _infT ^sup2 in the stable, nocturnal boundary layer (NBL) have been obtained with the analytic models described by Nieuwstadt (1984, 1985) and Sorbjan (1986) and the numerical model of Duynkerke and Driedonks (1987). These theoretical profiles are compared with observed profiles from the meteorological mast at Cabauw, The Netherlands. From the observations, it is found that C _infT ^sup2 is large in the surface layer and small at the top of the NBL. Observations during nights with moderate geostrophic winds or during the first few hours of nights with a high geostrophic wind show a continuous decrease of C _infT ^sup2 from the surface layer to the top of the NBL. Observations made later on nights with a high geostrophic wind show the development of a maximum of C _infT ^sup2 at about three quarters of the NBL. From the comparison with the models, we conclude that the observed profiles are most satisfactorily described by the model of Duynkerke and Driedonks.     

An approximate analytical solution for the deposition of heavy particles released from an elevated line source

This work re-examines and further develops an analytical solution for the deposition swath of heavy particles released in the atmosphere from an elevated source over uniform terrain, correcting the particle diffusivity for the crossing trajectory effect. The revised (approximate) analytical solution proves to be accurate within 20% over a wide range of micrometeorological conditions and particle size, despite its neglect of the turbulence component of the deposition flux. It compares very satisfactorily with experimental data and with the simulations of a Lagrangian stochastic model, provided the variable U(H)/w _g ≤7 (ratio of the mean horizontal wind speed at source height to the particle settling velocity). In this domain of validity, simple formulae relating the statistics of the deposition swath to U(H)/w _g are derived.     

Dependence of geostrophic drag on intensity of convection,baroclinicity, and acceleration

This paper examines the practical importance of stability, baroclinicity, and acceleration in the bulk ABL similarity formulations, in light of the random errors inherent in field measurements. This is done by propagating the measurement uncertainties through a theoretical model for the bulk ABL similarity functionsA ₀ andB ₀, under a range of assumed (but always unstable) conditions. It is shown that random measurement errors and acceleration effects may overwhelm most effects of baroclinicity and stability, once conditions are at least slightly unstable. Because of this, it is hard to discern a clear functional dependence ofA ₀ andB ₀ on degree of instability. Thus, for a given value ofh _i/z₀, whereh _i is the inversion height andz ₀ is the surface roughness, the geostrophic drag coefficient, which depends onA ₀ andB ₀, and weakly onh _i/z₀, may also be taken to be nearly independent of degree of instability.