Deliquescence induces eddy covariance and estimable dry deposition errors

The effects of deliquescent processes on eddy covariance estimates of dry deposition are considered, both theoretically and using some representative data. Turbulent and latent heat fluxes near the Earth's surface imply a vertical “saturation ratio flux” (Fairall, C.W., Atmospheric Environment, 18 (1984) 1329). Deliquescence, instantaneous relative to turbulent time-scales, responds to saturation ratio (or relative humidity) fluctuations and induces covariance between vertical winds and size-resolved number concentrations of hygroscopic aerosols. The induced covariance represents an error in the estimation of surface exchange by direct application of the eddy covariance technique. Under deliquescent conditions (high relative humidity), resulting errors in dry deposition estimates are shown to be often as large as typical deposition velocities reported for small particles, depending on the shape of the number distribution and the magnitudes of heat and vapor fluxes in the boundary layer.     

Particle fluxes above forests: observations, methodological considerations and method comparisons

This paper reports a study designed to test, evaluate and compare micro-meteorological methods for determining the particle number flux above forest canopies. Half-hour average particle number fluxes above a representative broad-leaved forest in Denmark derived using eddy covariance range from -7x10(7) m(-2) s(-1) (1st percentile) to 5x10(7) m(-2) s(-1) (99th percentile), and have a median value of -1.6x10(6) m(-2) s(-1). The statistical uncertainties associated with the particle number flux estimates are larger than those for momentum fluxes and imply that in this data set approximately half of the particle number fluxes are not statistically different to zero. Particle number fluxes from relaxed eddy accumulation (REA) and eddy covariance are highly correlated and of almost identical magnitude. Flux estimates from the co-spectral and dissipation methods are also correlated with those from eddy covariance but exhibit higher absolute magnitude of fluxes.     

Growing season total gaseous mercury (TGM) flux measurements over an Acer rubrum L. stand

Relaxed eddy accumulation (REA) measurements of the total gaseous mercury (TGM) flux measurements were taken over a deciduous forest predominantly composed of Red Maple (Acer rubrum L.) during the growing season of 2004 and the second half of the growing season of 2005. The magnitudes of the flux estimates were in the range of published results from other micrometeorological mercury fluxes taken above a tall canopy and larger than estimates from flux chambers. The magnitude and direction of the flux were not static during the growing season. There was a significant trend (p < 0.001), from net deposition of TGM in early summer to net evasion in the late summer and early fall before complete senescence. A growing season atmosphere-canopy total mercury (TGM) compensation point during unstable daytime conditions was estimated at background ambient concentrations (1.41 ng m^?3). The trend in the seasonal net TGM flux indicates that long term dry deposition monitoring is needed to accurately estimate mercury loading over a forest ecosystem.     

Measurements and modelling of cloudwater deposition to moorland and forests

The objective of the study was to measure the size dependence of cloudwater deposition and associated average ionic fluxes to vegetated surfaces. Measurements were made over a forest canopy at Dunslair heights in south Scotland and a moorland site, Great Dun Fell, in northern England. Measurements were made using the gradient and eddy correlation techniques. Eddy correlation measurements were made using an ultrasonic anemometer, a Knollenberg forward scattering spectrometer probe (to measure liquid water fluxes and fluxes of droplets in 1 microm size intervals) and a GSI particulate volume monitor (to measure liquid-water fluxes). Measurements were made at Great Dun Fell of the size dependence of droplet deposition velocity, using the gradient technique with two Knollenberg probes. Simultaneous gradient and eddy correlation measurements were also made at Great Dun Fell of average cloud-water fluxes, together with chemical analysis of cloud water composition, using a continuous analysis system. At Dunslair Heights, eddy correlation measurements were made using both the Knollenberg and Gerber Scientific Instruments (GSI) probes, while simultaneous gradient measurements using two GSI probes were also attempted. Samples of cloud water were collected at Dunslair Heights, using passive string collectors for chemical analysis by ion chromatography. The major findings of the study were: 1. The droplet deposition velocities measured by the two techniques were similar. 2. The deposition velocities were a strong function of droplet size. Considerable resistance to deposition was evident for droplets of less than 5 microm radius. Deposition velocities for particles from about 6 to 8 microm exceeded those for momentum. 3. Except when the droplets were very small or the winds very light, bulk cloud-water deposition velocities were about 80% or more of the momentum deposition velocities to forests.     

Computing surface concentration fluxes of trace gases from a variational method using measured variance and single-level concentration data

Surface fluxes of O₃, CO₂ and SO₂ were estimated from a variational method by using measured concentrations and variances of these trace gases. The measurements were taken over a deciduous forest when it was fully leafed during the summer of 1988 and when it was leafless during the winter of 1990. A flux–variance relation and a flux–gradient relation were employed as constraints in a cost function which is minimized to find the optimal estimate of concentration fluxes of the gases under study. Fluxes of O₃, CO₂ and SO₂ from the variational method were compared with fluxes estimated by the flux–variance relation and measured using an eddy correlation technique. Results show that the variational method improves the estimates of fluxes.     

Aerosol properties,in-canopy gradients,turbulent fluxes and VOC concentrations at a pristine forest site in Amazonia

Aerosol physical and chemical properties were measured in a forest site in central Amazonia (Cuieiras reservation, 2.61S; 60.21W) during the dry season of 2004 (Aug–Oct). Aerosol light scattering and absorption, mass concentration, elemental composition and size distributions were measured at three tower levels (Ground: 2 m; Canopy: 28 m, and Top: 40 m). For the first time, simultaneous eddy covariance fluxes of fine mode particles and volatile organic compounds (VOC) were measured above the Amazonian forest canopy. Aerosol fluxes were measured by eddy covariance using a Condensation Particle Counter (CPC) and a sonic anemometer. VOC fluxes were measured by disjunct eddy covariance using a Proton Transfer Reaction Mass Spectrometer (PTR-MS). At nighttime, a strong vertical gradient of phosphorus and potassium in the aerosol coarse mode was observed, with higher concentrations at Ground level. This suggests a source of primary biogenic particles below the canopy. Equivalent black carbon measurements indicate the presence of light-absorbing aerosols from biogenic origin. Aerosol number size distributions typically consisted of superimposed Aitken (76 nm) and accumulation modes (144 nm), without clear events of new particle formation. Isoprene and monoterpene fluxes reached respectively 7.4 and 0.82 mg m^?2 s^?1 around noon. An average fine particle flux of 0.05 ± 0.10 10⁶ m^?2 s^?1 was calculated, denoting an equilibrium between emission and deposition fluxes of fine mode particles at daytime. No significant correlations were found between VOC and fine mode aerosol concentrations or fluxes.     

Apply appropriate models in the study of dry deposition fluxes of particulate matter and polycyclic aromatic hydrocarbons (PAHs) in Tsukuba, Japan

In this study, plates for downward flux and upward flux were used to measure atmospheric dry deposition fluxes for particulate mass and polycyclic aromatic hydrocarbons (PAHs) in TERC (Tsukuba), Japan. Ambient particles concentrations were also collected using a high-volume air sampler, and ambient particle size distributions between 0.01 μm and 13.1 μm were measured using a low-pressure cascade impactor to characterise the PAHs levels and dry deposition. The results indicated that the average cumulative fraction of dry deposition flux for particles and PAHs which attached with them was caused by the particle size of greater than 1.2-6.3 μm (97%).     

Measurement of the dry deposition flux of NH3 on to coniferous forest

The dry deposition flux of NH3 to coniferous forest was determined by the micrometeorological gradient method using a 36 m high tower. Aerodynamic characteristics of the site were studied, using a second tower erected in the forest 100 m from the first. Fluxes and gradients of heat and momentum measured on both towers indicated a fairly homogeneous turbulent flow field over the studied area of the forest. Site specific flux profile functions for heat were derived from continuous measurements of turbulent fluxes and gradients. These functions were used to derive fluxes from the observed gradients of NH3. In total, eighty 90-min NH3 flux runs were performed. The results indicate a strong nonstomatal uptake of NH3 by the forest. A representative dry deposition velocity for NH3 of 3.6 cm(-1) s was derived. The annual average flux was roughly estimated to be equivalent to 50 kg N ha(-1) yr, significantly higher than the critical load for coniferous forest.     

Atmospheric deposition of polychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans, and dioxin-like polychlorinated biphenyls in the Kanto Region, Japan

The atmospheric bulk (dry and wet) deposition of dioxins was investigated at four locations (Tokyo, Yokohama, Tsukuba, and Tanzawa) in the Kanto region (in Japan) over one year using a stainless-steel pot. Annual average polychlorinated dibenzo-p-dioxins/polychlorinated dibenzofurans (PCDD/PCDF) deposition fluxes were estimated to be from 450 to 1300 ng/m2/yr, and the annual average TEQ fluxes from 5.7 to 17 ng-TEQ/m2/yr at the four locations. The PCDD/PCDF deposition flux was higher in winter than in summer. The deposition flux could be related to ambient temperature, particularly for less chlorinated PCDDs/PCDFs, while the deposition flux is not necessarily related to the amount of precipitation. The PCDD/PCDF deposition flux increased as the particle deposition flux increased, for the winter samples. Based on the ratio of the PCDD/PCDF deposition fluxes to the particle deposition fluxes, the contribution of the reentrainment of soil particles to the TEQ of PCDD/PCDF deposition was considered to be negligible in this region. Based on the air concentrations monitored near our deposition sampling points by the municipalities, the ratio of the annual deposition flux to the annual average air concentration was roughly estimated to be 0.082 cm/s. The range of deposition flux in the Kanto region was estimated to be from 1.5 to 31 (median: 9.8) ng-TEQ/m2/yr based on the range of air concentration data measured by the municipalities. The total annual deposition flux in the entire Kanto region was estimated to range from 50 to 900 g-TEQ/yr (median 320 g-TEQ/yr). This estimated flux was of the same order as the sum of estimated emissions from municipal solid waste incinerators and industrial waste incinerators in the Kanto region. The contributions of dioxin-like PCBs in Yokohama, Tsukuba, and Tanzawa depositions were less than 10% of the total TEQ; however, in Tokyo it was almost equal to or more than 50%.     

Particle deposition in low-speed,high-turbulence flows

The experimental and numerical study considers the concentration of airborne particulate contaminants, such as spores of spoilage fungi, and their deposition on a surface, in a Petri dish, and on a warm box-shaped product placed in a food-processing environment. Field measurements by standard, active and passive samplers provide typical values of airborne concentrations and specific deposition fluxes. Velocity and turbulence data from field studies are used as input in large eddy simulations of the process, and estimates of deposition fluxes are of the same order of magnitude as those deduced from field measurements. Particle deposition is shown to be associated with near-wall coherent structures. Flow reversal, simulated by impulsive start, is shown to give higher deposition rates than steady mean flows.     

Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

Elevational patterns of sulfur deposition at a site in the Catskill Mountains,New York

In this paper we report measurements of SO^2-₄ fluxes in throughfall and bulk deposition across an elevational transect from 800 to 1275 m on Slide Mountain in the Catskill Mountains of southeastern New York State. The net throughfall flux of SO^2-₄ (throughfall-bulk deposition), which we attribute to cloud and dry deposition, increased by roughly a factor of 13 across this elevational range. Part of the observed increase results from the year-round exposure of evergreen foliage at the high-elevation sites, compared to the lack of foliage in the dormant season in the deciduous canopies at low elevations. Comparison of the net throughfall flux with estimates of cloud deposition suggests that both cloud deposition and dry deposition increased with elevation. Dry deposition estimates from a nearby monitoring site fall within the measured range of net throughfall flux for SO^2-₄. The between-site variation in net throughfall flux was very high at the high-elevation sites, and less so at the lower sites, suggesting that studies of atmospheric deposition at high-elevations will be complicated by extreme spatial variability in deposition rates. Studies of atmospheric deposition in mountainous areas of the eastern U.S. have often emphasized cloud water deposition, but these results suggest that elevational increases in dry deposition may also be important.     

Dry deposition of particles to wave surfaces: II. Wind tunnel experiments

Wind tunnel measurements of particle dry deposition to wavy and flat surfaces were made to estimate the enhancement of deposition rates due to waves on water surfaces. Measurements were made of 4.0 and 6.7 μm uranine particles at wind speeds of 5 and 10 m s⁻¹ to sinusoidal waves with height to length ratios 2a/λ=0.1 and 0.03 and to flat surfaces. Results showed that deposition was greatest to the upslope portion of the wave, accounting for 40–45% of the total mass, followed by the trough (30%), downslope (15%), and crest (10–15%). These results generally agreed within experimental variability with modeling predictions (Zufall et al., 1999). Deposition was enhanced at the upslope due to the effects of particle interception and impaction on the wave. Total deposition to the wave surfaces was greater than deposition to the flat surface for a large majority of the cases. The average increase in deposition to both wave surfaces for both particle sizes and wind speeds over deposition to the flat surface was 80%.     

Measuring eddy covariance fluxes of ozone with a slow-response analyser

Ozone (O₃) fluxes above a temperate mountain grassland were measured by means of the eddy covariance (EC) method using a slow-response O₃ analyser. The resultant flux loss was corrected for by a series of transfer functions which model the various sources of high- and, in particular, low-pass filtering. The resulting correction factors varied on average between 1.7 and 3.5 during night and daytime, respectively. A cospectral analysis confirmed the accuracy of this approach. O₃ fluxes were characterised by a comparatively large random uncertainty, which during daytime typically amounted to 60%. EC O₃ fluxes were compared against O₃ flux measurements made concurrently with the flux-gradient (FG) method. The two methods generally agreed well, except for a period between sunrise and early afternoon, when the FG method was suspected of being affected by the presence of photochemical sources/sinks. O₃ flux magnitudes and deposition velocities determined with the EC method compared nicely with the available literature from grassland studies. We conclude that our understanding of the causes and consequences of various sources of flux loss (associated with any EC system) has sufficiently matured so that also less-than-ideal instrumentation may be used in EC flux applications, albeit at the cost of relatively large empirical corrections.     

Conditional sampling for measuring mercury vapor fluxes

Surface–atmosphere mercury fluxes are difficult to measure accurately. Current techniques include dynamic flux chambers and micrometeorological gradient and aerodynamic approaches, all of which have a number of intrinsic problems associated with them. We have adapted conditional sampling (relaxed eddy accumulation), a micrometeorological technique commonly used to measure other trace gas fluxes, to measure surface–air mercury fluxes. Our initial flux measurement campaign over an agricultural soil consisted of two 1-week measurement periods, and was longer in duration than previously reported mercury flux measurement periods. Fluxes during both measurement periods ranged between 190.5 (evolution) and –91.7 ng m⁻² h⁻¹ (deposition) with an average evolution of 9.67 ng m⁻² h⁻¹. The data showed significant diurnal trends, weakly correlated with shallow soil temperatures and solar radiation. This initial trial run indicates that conditional sampling has much promise for the accurate quantification of both short and long-term mercury fluxes.     

Vertical transport of accumulation mode particles between two street canyons and the urban boundary layer

Concentrations and turbulent fluxes of accumulation mode particles were measured during the 2004–2005 ‘Canopy and Aerosol Particle Interaction in Toulouse Urban Layer’ project (CAPITOUL) at the top of two intersecting street canyons and in the urban boundary layer (UBL) in Toulouse, France. Particle numbers were strongly affected by boundary layer depth and showed limited sensitivity to local emissions. Differences in the diurnal patterns of particle numbers were observed between the finer fraction (0.3–0.4 μm) and coarser fraction (1.6–2.0 μm) of accumulation mode particles, indicating different processes of formation, evolution and transportation may be dominant. Highest particle numbers were observed in the narrow street canyon which had more limited local emissions and comparatively small particle fluxes. However, the improved ventilation rate in the wider canyon was also associated with the downward mixing of particles into the street canyon from the UBL. The results from this study clearly illustrate the temporal and spatial variability of particle numbers and fluxes in the urban atmosphere.     

Modeling the atmospheric transport and deposition of polychlorinated dibenzo-p-dioxins and dibenzofurans in North America

The atmospheric fate of polychlorinated dibenzo-p-dioxins (PCDDs) and polychlorinated dibenzofurans (PCDFs) was simulated for the year 2000 in North America using a SMOKE/CMAQ-based chemical transport model that was modified for this purpose. The 1999 USEPA emission inventories of PCDD/Fs and criteria pollutants were used. The 1995 Canadian emission inventory of criteria pollutants and the 1995 Canadian area source emissions for PCDD/Fs were used with the 2000 Canadian point source emissions. Modifications to CMAQ involved coupling it with dual organic matter (OM) absorption and black carbon (BC) adsorption models to calculate PCDD/F gas–particle partitioning. The model satisfactorily reproduced the particle bound fractions at all rural sites for which there were measured data and across the whole domain, the modeled vs. measured differences in particle bound fractions were less than 20% for nearly all congeners. The model predicted ambient air PCDD/F concentrations were also consistent with measurements. Simulated deposition fluxes were within 58% of direct measurements. PCDD/F atmospheric depositions to each of the Great Lakes were estimated for the year 2000. The results indicate that approximately 76% of the total deposition of PCDD/Fs to the Great Lakes (in W-TEQ, or toxic equivalent units as defined by the World Health Organization) is attributed to PCDD/Fs absorbed into OM in aerosol. For all of the lakes, more than 92% of all deposition is particle phase wet deposition and only 5–8% is particle phase dry deposition. Wet deposition from the gas phase is negligible. Of the 17 toxic PCDD/F congeners, the Cl_4–5DD/F compounds contribute approximately 70% to the total atmospheric deposition to the Great Lakes. The seasonal changes in the PCDD/F deposition flux track variations in ambient temperature.     

Spatial variability of throughfall fluxes in a spruce forest

The spatial variability of throughfall deposition of H(+), Ca(2+), Mg(2+), Na(+), K(+), Cl(-), NO(3)(-), NH(4)(+), O(4)(2-) to a Norway spruce (Picea abies (L.) Karst.) forest was intensively examined during the period October 1986 to October 1987. Large systematic spatial variability of the atmospheric deposition within the forest was observed. The flux of throughfall water was higher away from the trunk compared to the flux close to the trunk. In contrast to this, the deposition of all substances was considerably higher close to the trunk compared to the deposition at the periphery of the canopy. A linear decrease in deposition as a function of the distance from the nearest tree trunk was found. Further, the deposition varied quite dramatically between trees according to their size. The observed spatial variability in throughfall may be due to variabilities in the processes taking part in altering the distribution and composition of the precipitated water as it moves through the canopy. The influence of these processes of precipitation, wash-off, dry deposition and canopy exchange is discussed, and it is found that both increased dry deposition and canopy exchange in the tree tops contribute to the higher solute fluxes found close to the tree trunk.     

Remote sampling of a CO2 point source in an urban setting

High frequency CO₂ and wind speed measurements were used to examine the urban baseline eddy covariance CO₂ flux and analyse the CO₂ rich plume from a local power station. A reliable relationship between high frequency CO₂ maxima and the rate of CO₂ emission at the power station was established. This relationship was shown to be highly dependant on wind speed. The ensemble mean plume was found to be Gaussian in horizontal profile with a width dependant on wind speed. The relationship between peak CO₂ mixing ratio and averaging time was shown to be a simple power law with a time exponent of approximately 0.5. The large, short pulses in CO₂ mixing ratio in the power plant plume were found to have an approximately Lorentzian shape. These pulses generated negative vertical eddy flux measurements so data from the plume sector were necessarily excluded from the flux baseline results. The plume-excluded flux had a similar magnitude and variability to those reported in other urban CO₂ flux studies despite this site not being ideal due to the proximity of roughness elements to the measurement point.     

Limitations of long-path averaging instruments

Long-path averaging instruments measure the average velocity or concentration of a substance or substances over an averaging path. These measurements are then often used for calculation of the average concentration and mass flow rate of the substance. The purpose of this paper is to describe some of the limitations of these instruments and to suggest ways in which these limitations can be minimized. Two limitations were examined: measuring concentration in a single dimension (e.g., ignoring the variation in concentration over the width of the sample plane), and deriving an average concentration without considering velocity effects. The resultant errors will be application-specific. Estimates of the second source of error can be obtained from the covariance of concentration and velocity profiles over the path length. Unfortunately, suitable field data were not available, and to illustrate the method, estimates of the error were obtained for a range of possible concentration and velocity profiles. Errors of 50% or greater in the mass flow were incurred for the concentration and velocity profiles considered. This error was reduced to a negligible level by segmenting the averaging path length. It is recommended that velocity and concentration profiles be obtained for a broad range of applications to enable the importance of covariance errors to be better assessed.