The use of disjunct eddy sampling methods for the determination of ecosystem level fluxes of trace gases

The concept of disjunct eddy sampling (DES) for use in measuring ecosystem-level micrometeorological fluxes is re-examined. The governing equations are discussed as well as other practical considerations and guidelines concerning this sampling method as it is applied to either the disjunct eddy covariance (DEC) or disjunct eddy accumulation (DEA) techniques. A disjunct eddy sampling system was constructed that could either be combined with relatively slow sensors (response time of 2 to 40 s) to measure fluxes using DEC, or could also be used to accumulate samples in stable reservoirs for later laboratory analysis (DEA technique). Both the DEC and DEA modes of this sampler were tested against conventional eddy covariance (EC) for fluxes of either CO2 (DEC) or isoprene (DEA). Good agreement in both modes was observed relative to the EC systems. However, the uncertainty in a single DEA flux measurement was considerable (40%) due to both the reduced statistical sampling and the analytical precision of the concentration difference measurements. We have also re-investigated the effects of nonzero mean vertical wind velocity on accumulation techniques as it relates to our DEA measurements. Despite the higher uncertainty, disjunct eddy sampling can provide an alternative technique to eddy covariance for determining ecosystem-level fluxes for species where fast sensors do not currently exist

Some aspects of the energy balance closure problem

International audienceAfter briefly discussing several reasons for the energy balance closure problem in the surface layer, the paper focuses on the influence of the low frequency part of the turbulence spectrum on the residual. Changes in the turbulent fluxes in this part of the turbulence spectrum were found to have a significant influence on the changes of the residual. Using the ogive method, it was found that the eddy-covariance method underestimates turbulent fluxes in the case of ogives converging for measuring times longer than the typical averaging interval of 30 min. Additionally, the eddy-covariance method underestimates turbulent fluxes for maximal ogive functions within the averaging interval, both mainly due to advection and non-steady state conditions. This has a considerable influence on the use of the eddy-covariance method

Eddy covariance raw data processing for CO2 and energy fluxes calculation at ICOS ecosystem stations

open18siThe eddy covariance is a powerful technique to estimate the surface-Atmosphere exchange of different scalars at the ecosystem scale. The EC method is central to the ecosystem component of the Integrated Carbon Observation System, a monitoring network for greenhouse gases across the European Continent. The data processing sequence applied to the collected raw data is complex, and multiple robust options for the different steps are often available. For Integrated Carbon Observation System and similar networks, the standardisation of methods is essential to avoid methodological biases and improve comparability of the results. We introduce here the steps of the processing chain applied to the eddy covariance data of Integrated Carbon Observation System stations for the estimation of final CO2, water and energy fluxes, including the calculation of their uncertainties. The selected methods are discussed against valid alternative options in terms of suitability and respective drawbacks and advantages. The main challenge is to warrant standardised processing for all stations in spite of the large differences in e.g. ecosystem traits and site conditions. The main achievement of the Integrated Carbon Observation System eddy covariance data processing is making CO2 and energy flux results as comparable and reliable as possible, given the current micrometeorological understanding and the generally accepted state-of-The-Art processing methodsopenSabbatini, Simone; Mammarella, Ivan; Arriga, Nicola; Fratini, Gerardo; Graf, Alexander; Hörtnagl, Lukas; Ibrom, Andreas; Longdoz, Bernard; Mauder, Matthias; Merbold, Lutz; Metzger, Stefan; Montagnani, Leonardo; Pitacco, Andrea; Rebmann, Corinna; Sedlák, Pavel; Šigut, Ladislav; Vitale, Domenico; Papale, DarioSabbatini, Simone; Mammarella, Ivan; Arriga, Nicola; Fratini, Gerardo; Graf, Alexander; Hörtnagl, Lukas; Ibrom, Andreas; Longdoz, Bernard; Mauder, Matthias; Merbold, Lutz; Metzger, Stefan; Montagnani, Leonardo; Pitacco, Andrea; Rebmann, Corinna; Sedlák, Pavel; Šigut, Ladislav; Vitale, Domenico; Papale, Dari

Lag time determination in DEC measurements with PTR-MS

The disjunct eddy covariance (DEC) method has emerged as a popular technique for micrometeorological flux measurements of volatile organic compounds (VOCs). It has usually been combined with proton transfer reaction mass spectrometry (PTR-MS), an online technique for VOC concentration measurements. However, the determination of the lag time between wind and concentration measurements has remained an important challenge. To address this issue, we studied the effect of different lag time methods on DEC fluxes. The analysis was based on both actual DEC measurements with PTR-MS and simulated DEC data derived from high frequency H<sub>2</sub>O measurements with an infrared gas analyzer. Conventional eddy covariance fluxes of H<sub>2</sub>O served as a reference in the DEC simulation. The individual flux measurements with PTR-MS were rather sensitive to the lag time methods, but typically this effect averaged out when the median fluxes were considered. The DEC simulation revealed that the maximum covariance method was prone to overestimation of the absolute values of fluxes. The constant lag time methods, one based on a value calculated from the sampling flow and the sampling line dimensions and the other on a typical daytime value, had a tendency to underestimate. The visual assessment method and our new averaging approach utilizing running averaged covariance functions did not yield statistically significant errors and thus fared better than the habitual choice, the maximum covariance method. Given this feature and the potential for automatic flux calculation, we recommend using the averaging approach in DEC measurements with PTR-MS. It also seems well suited to conventional eddy covariance applications when measuring fluxes near the detection limit

Energy balance closure for the LITFASS-2003 experiment

In the first part, this paper synthesises the main results from a series of previous studies on the closure of the local energy balance at low-vegetation sites during the LITFASS-2003 experiment. A residual of up to 25% of the available energy has been found which cannot be fully explained either by the measurement uncertainty of the single components of the surface energy balance or by the length of the flux-averaging period. In the second part, secondary circulations due to heterogeneities in the surface characteristics (roughness, thermal and moisture properties) are discussed as a possible cause for the observed energy balance non-closure. This hypothesis seems to be supported from the fluxes derived from area-averaging measurement techniques (scintillometers, aircraft)

Stochastic superparameterization in quasigeostrophic turbulence

In this article we expand and develop the authors' recent proposed methodology for efficient stochastic superparameterization (SP) algorithms for geophysical turbulence. Geophysical turbulence is characterized by significant intermittent cascades of energy from the unresolved to the resolved scales resulting in complex patterns of waves, jets, and vortices. Conventional SP simulates large scale dynamics on a coarse grid in a physical domain, and couples these dynamics to high-resolution simulations on periodic domains embedded in the coarse grid. Stochastic SP replaces the nonlinear, deterministic eddy equations on periodic embedded domains by quasilinear stochastic approximations on formally infinite embedded domains. The result is a seamless algorithm which never uses a small scale grid and is far cheaper than conventional SP, but with significant success in difficult test problems. Various design choices in the algorithm are investigated in detail here, including decoupling the timescale of evolution on the embedded domains from the length of the time step used on the coarse grid, and sensitivity to certain assumed properties of the eddies (e.g. the shape of the assumed eddy energy spectrum). The different design choices are compared and contrasted on a stringent test suite for quasigeostrophic turbulence involving two-layer dynamics on a beta-plane forced by an imposed background shear. The algorithms developed are expected to be particularly useful in providing accurate and efficient stochastic parameterizations for use in ensemble-based state estimation and prediction

Evening Methane Emission Pulses from a Boreal Wetland Correspond to Convective Mixing in Hollows

Spatial and temporal heterogeneity of methane flux from boreal wetlands makes prediction and up-scaling challenging, both within and among wetland systems. Drivers of methane production and emissions are also highly variable, making empirical model development difficult and leading to uncertainty in methane emissions estimates from wetlands. Previous studies have examined this problem using point-scale (static chamber method) and ecosystem-scale (flux tower methods) measurements, but few studies have investigated whether different processes are observed at these scales. We analyzed methane emissions from a boreal fen, measured by both techniques, using data from the Boreal Ecosystem-Atmosphere Study. We sought to identify driving processes associated with methane emissions at two scales and explain diurnal patterns in emissions measured by the tower. The mean methane emission rates from flux chambers were greater than the daytime, daily mean rates measured by the tower, but the nighttime, daily mean emissions from the tower were often an order of magnitude greater than emissions recorded during the daytime. Thus, daytime measurements from either the tower or chambers would lead to a biased estimate of total methane emissions from the wetland. We found that the timing of nighttime emission events was coincident with the cooling and convective mixing within hollows, which occurred regularly during the growing season. We propose that diurnal thermal stratification in shallow pools traps methane by limiting turbulent transport. This methane stored during daytime heating is later released during evening cooling due to convective turbulent mixing

Chemical sensor resolution requirements for near-surface measurements of turbulent fluxes

Businger and Delany (1990) presented an approach to estimate the sensor resolution required to limit the contribution of the uncertainty in the chemical concentration measurement to uncertainty in the flux measurement to 10 % for eddy covariance, gradient, and relaxed eddy accumulation flux measurement methods. We describe an improvement to their approach to estimate required sensor resolution for the covariance method, and include disjunct eddy covariance. In addition, we provide data to support selection of a form for the dimensionless scalar standard deviation similarity function based on observations of the variance of water vapor fluctuations from recent field experiments. We also redefine the atmospheric parameter of Businger and Delany in a more convenient, dimensionless form. We introduce a "chemical parameter" based on transfer velocity parameterizations. Finally, we provide examples in which the approach is applied to measurement of carbon dioxide, dimethylsulfide, and hexachlorobenzene fluxes over water. The information provided here will be useful to plan field measurements of atmosphere-surface exchange fluxes of trace gases