Comparison of constant volume balloons, model trajectories and tracer transport during ETEX

In the field phases of the European Tracer EXperiment (ETEX), an inert tracer was released for 12 h into the atmosphere and samples taken at several locations downwind. During the same time, several Constant Volume Balloons (CVB) (10 and 6 for ETEX first and second release, respectively) were launched into different altitudes and followed as far as 21–188 km, to indicate the initial dispersion directions of the tracer puff. A model simulating the CVB behaviour in hydrostatic meso-scale model forecasts is applied to ETEX data to demonstrate its capability to predict the tracer puff mean axis over long distances (−2000 km). CVB model results are first compared to air parcels trajectories and 2D (i.e. isentropic, isobaric and isodensity) trajectories. Then they are compared to the measured CVB trajectories and finally to the tracer puff trajectories. As expected, the CVB model and isodensity model trajectories are found to be identical. The 16 CVBs calculated trajectories nearly overlap the real ones over 21–188 km with mean absolute horizontal transport deviations less than 20 km (average value of 8.2 km). The corresponding relative transport deviations are less than 45% with an average value of 20.6%. Better predictions are obtained for the ETEX second release. During the 60 h following ETEX’s first release start, the simulated CVBs are mainly found in the area of the maximum surface concentrations of the released tracer, up to 2000 km. Up to 36 h after ETEX second tracer release start, the simulated CVB trajectories predict well the mean axis of the tracer puff, but failed later.     

Evaluation of the effect of meteorological data resolution on Lagrangian particle dispersion simulations using the ETEX experiment

This paper presents results from a series of numerical experiments designed to evaluate operational long-range dispersion model simulations, and to investigate the effect of different temporal and spatial resolution of meteorological data from numerical weather prediction models on these simulations. Results of Lagrangian particle dispersion simulations of the first tracer release of the European Tracer Experiment (ETEX) are presented and compared with measured tracer concentrations. The use of analyzed data of higher resolution from the European Center for Medium-Range Weather Forecasts (ECMWF) model produced significantly better agreement between the concentrations predicted with the dispersion model and the ETEX measurements than the use of lower resolution Navy Operational Global Atmospheric Prediction System (NOGAPS) forecast data. Numerical experiments were performed in which the ECMWF model data with lower vertical resolution (4 instead of 7 levels below 500 mb), lower temporal resolution (12 h instead of 6 h intervals), and lower horizontal resolution (2.5° instead of 0.5°) were used. Degrading the horizontal or temporal resolution of the ECMWF data resulted in decreased accuracy of the dispersion simulations. These results indicate that flow features resolved by the numerical weather prediction model data at approximately 45 km horizontal grid spacing and 6 h time intervals, but not resolved at 225 km spacing and 12 h intervals, made an important contribution to the long-range dispersion.     

The field campaigns of the European Tracer Experiment (ETEX): overview and results

As part of the European Tracer Experiment (ETEX) two successful atmospheric experiments were carried out in October and November, 1994. Perfluorocarbon (PFC) tracers were released into the atmosphere in Monterfil, Brittany, and air samples were taken at 168 stations in 17 European countries for 72 h after the release. Upper air tracer measurements were made from three aircraft. During the first experiment a westerly air flow transported the tracer plume north-eastwards across Europe. During the second release the flow was eastwards. The results from the ground sampling network allowed the determination of the cloud evolution as far as Sweden, Poland and Bulgaria. This demonstrated that the PFT technique can be successfully applied in long-range tracer experiments up to 2000 km. Typical background concentrations of the tracer used are around 5–7 fl ?^-1 in ambient air. Concentrations in the plume ranged from 10 to above 200 fl/?^-1. The tracer release characteristics, the tracer concentrations at the ground and in upper air, the routine and additional meteorological observations at the ground level and in upper air, trajectories derived from constant-level balloons and the meteorological input fields for long-range transport models are assembled in the ETEX database. The ETEX database is accessible via the Internet. Here, an overview is given of the design of the experiment, the methods used and the data obtained.     

Meteorological conditions at the release site during the two tracer experiments

The state of the boundary layer at the release site during the two tracer experiments is described, based on measurements of atmospheric turbulence carried out by a sonic anemometer, profiles of horizontal and vertical wind by a SODAR as well as frequent radiosonde releases. The boundary layer height is derived from radiosoundings, modelled and discussed. The study is meant as background for discussions of the tracer behaviour near the release site, and possible influence of the local meteorology at the release site on the long-range dispersion. The difference in the initial meso-meteorological conditions between the two experiments is discussed and some of its consequences for the dispersion of the released tracer are clarified. It is concluded that the first experiment is well suited for validation of long-range transport and dispersion models, whereas the second experiment calls for further model development. It is concluded that the micrometeorological measurements are very important as background information for the understanding of plume behaviour on small as well as long scales.     

Validation of the operational emergency response model at the Swedish Meteorological and Hydrological Institute using data from etex and the chernobyl accident

The Eulerian atmospheric tracer transport model MATCH (Multiscale Atmospheric Transport and Chemistry model) has been extended with a Lagrangian particle model treating the initial dispersion of pollutants from point sources. The model has been implemented at the Swedish Meteorological and Hydrological Institute in an emergency response system for nuclear accidents and can be activated on short notice to provide forecast concentration and deposition fields.The model has been used to simulate the transport of the inert tracer released during the ETEX experiment and the transport and deposition of ¹³⁷Cs from the Chernobyl accident. Visual inspection of the results as well as statistical analysis shows that the extent, time of arrival and duration of the tracer cloud, is in good agreement with the observations for both cases, with a tendency towards over-prediction for the first ETEX release. For the Chernobyl case the simulated deposition pattern over Scandinavia and over Europe as a whole agrees with observations when observed precipitation is used in the simulation. When model calculated precipitation is used, the quality of the simulation is reduced significantly and the model fails to predict major features of the observed deposition field.     

Droplet trajectories and deposition with an airblast forestry sprayer

Abstract

The characteristics of the air pattern and spray droplet trajectories from the Algonquin airblast forestry sprayer were investigated. An experiment with water and Rhodamine tracer dye was conducted in a shelterwood location to evaluate droplet deposition on Kromekote^® card targets above and within a canopy of poplar saplings. Computer simulation of a free jet was used to model the airblast in order to investigate droplet trajectories. Measured deposits indicated that the effective swath width is 14 m from the blower outlet. Mean simulated droplet trajectory data were in good agreement with the droplet density data observed in the field. Both sets of results predicted a very low level of spray deposition beyond 14 m.     

Validation of the lagrangian particle dispersion model FLEXPART against large-scale tracer experiment data

A comprehensive validation of FLEXPART, a recently developed Lagrangian particle dispersion model based on meteorological data from the European Centre for Medium-Range Weather Forecasts, is described in this paper. Measurement data from three large-scale tracer experiments, the Cross-Appalachian Tracer Experiment (CAPTEX), the Across North America Tracer Experiment (ANATEX) and the European Tracer Experiment (ETEX) are used for this purpose. The evaluation is based entirely on comparisons of model results and measurements paired in space and time. It is found that some of the statistical parameters often used for model validation are extremely sensitive to small measurement errors and should not be used in future studies. 40 cases of tracer dispersion are studied, allowing a validation of the model performance under a variety of different meteorological conditions. The model usually performs very well under undisturbed meteorological conditions, but it is less skilful in the presence of fronts. The two ETEX cases reveal the full range of the model’s skill, with the first one being among the best cases studied, and the second one being, by far, the worst. The model performance in terms of the statistical parameters used stays rather constant with time over the periods (up to 117 h) studied here. It is shown that the method used to estimate the concentrations at the receptor locations has a significant effect on the evaluation results. The vertical wind component sometimes has a large influence on the model results, but on the average only a slight improvement over simulations which neglect the vertical wind can be demonstrated. Subgrid variability of mixing heights is important and must be accounted for.     

Measurements of wind velocities in a street canyon

Measurements of wind velocities in an urban street canyon with a height-to-width ratio of about 1.4 were made when ambient winds aloft were approximately perpendicular to the street. The mean wind velocities were determined by analysis of trajectories of tracer balloons that were released in the canyon and photographed in rapid sequence. The trajectories indicate the presence of a primary vortex cell within the canyon, provided the ambient wind velocity exceeds 1.5–2.0 m s⁻¹. Measurements with hot-wire anemometers suggest that vehicular traffic at the street level is significant in increasing turbulence up to heights of approximately 7 m.     

A replacement for simple back trajectory calculations in the interpretation of atmospheric trace substance measurements

Trajectory calculations are often used for the interpretation of atmospheric trace substances measurements. However, two important effects are normally neglected or not considered systematically: first, measurements of trace substances sample finite volumes of air, whereas a trajectory tracks the path of an infinitesimally small particle; second, turbulence and convection. Advection by the deformative synoptic-scale atmospheric flow is responsible for the fact that a compact measurement volume is in fact turned into filamentary structures at earlier times, which a single trajectory cannot represent. Turbulence and convection add to this by causing a growth of the volume (backward in time) where processes such as emissions can affect measured concentrations. In this paper, we show that both effects are substantial and may be the largest sources of error when using trajectory calculations to establish source–receptor relationships. We use backward simulations with a Lagrangian particle dispersion model (LPDM) and cluster analysis of the particle positions to derive more representative single “trajectories” (transport paths) and trajectory ensembles. This reduces errors caused by filamentation and backward growth of the measurement volume to a few percent as compared to using a single, mean-wind trajectory and also yields estimates of the spread of the region of influence. Thus, we recommend to replace simple back trajectory calculations for interpretation of atmospheric trace substances measurements in the future by backward simulations with LPDMs, possibly followed by the clustering of particle positions as introduced in this paper.     

A sulfur hexafluoride tracer study at the tracy power plant—I. Preliminary results

A plume tracer and flow visualization experiment has been conducted in the vicinity of the Tracy Power Station near Reno, Nevada during the period 7–20 November 1983. The experiment was undertaken as a preliminary project to help design a larger, full-scale experiment planned for the future. The major objective of the study was to determine the extent of plume impact on the ground, and to ascertain how plume impact is influenced by both terrain and meteorological factors. The tracer experiment involved the release of gaseous SF₆ through the stack of the Tracy Power Station, and the deployment of 53 syringe samplers distributed in an area within 10 km of the release point. Supplementing the tracer measurements were meteorological, remote sensing, and photographic measurements, all designed to visualize the dispersion of the plume. Ten tests were conducted, each lasting a period of up to 9 h. Hourly samples were collected and analyzed for SF₆. Five of the tests resulted in significant plume touchdown within the test area. Although at elevations above the surrounding terrain winds tend to be strong and blow predominantly from the west, SF₆ concentrations were frequently observed both east and west of the power plant. Concentration maxima were localized at discrete locations and altitudes depending on the local meteorological conditions.     

Comparison of line– and point–source releases of tracer gases

Field measurements were made of greenhouse gas emissions from a wastewater treatment system using open path monitoring with detection by FTIR spectroscopy. Emission rates were determined by the ratio technique using a sulfur hexafluoride tracer gas released from a line source. As a quality control check, a second tracer gas – ethylene – was released from various single point locations. This paper presents a comparison of the line-source and point-source tracer releases for approximating emissions from the area source. The two types of tracer release showed excellent agreement when both release points were two hundred meters from the FTIR beam path. Data for other release points also were comparable, once differences in vertical dispersion as a function of distance are taken into account.     

Macrodispersion by diverging radial flows in randomly heterogeneous porous media

Radial flow takes place in a heterogeneous porous formation where the transmissivity T is modelled as a stationary random space function (RSF). The steady flow is driven by a given rate, and the mean velocity is radial. A pulse-like of a tracer is injected in the porous formation, and the thin plume spreads due to the fluctuations of the velocity which results a RSF as well. Transport is characterized by the mean front, and by the second spatial moment of the plume. We are primarily interested in tracer macrodispersion modelling. With the neglect of pore-scale dispersion, macrodispersion coefficients are computed at the second order of approximation, without neglecting the head-gradient fluctuations. Although transport is non-ergodic at the source, it is shown that ergodicity is achieved at small distances from the source. This is due to the fact that close to the source local velocities are quite large, and therefore solute particles become uncorrelated very soon. Under ergodic conditions, we compare macrodispersion mechanism in radial flows with that occurring in mean uniform flows. At short distances the spreading effect is highly enhanced by the large variability of the flow field, whereas at large distances transport exhibits a lesser dispersion due to the reduction of velocities. This supports the explanation provided by Indelman and Dagan (1999) to justify why the macrodispersivity is found smaller than that pertaining to mean uniform flows. The model is tested against a tracer transport experiment (Fernàndez-Garcia et al., 2004) by comparing the theoretical and experimental breakthrough curves. The accordance with real data, that is achieved without any fitting to concentration values, strengthens the capability of the proposed model to grasp the main features of such an experiment, the approximations as well as experimental uncertainties notwithstanding.     

Analysis of colloid and tracer breakthrough curves

We consider the dispersion and elution of colloids and dissolved nonsorbing tracers within saturated heterogeneous porous media. Since flow path geometry in natural systems is often ill-characterized macroscopic (mean) flow rates and dispersion tensors are utilized in order to account for the sub-model scale microscopic fluctuations in media structure (and the consequent hydrodynamic profile). Even for tracer migration and dispersal this issue is far from settled.Here we consider how colloid and tracer migration phenomena can be treated consistently. Theoretical calculations for model flow geometries yield two quantitative predictions for the transport of free (not yet captured) colloids with reference to a non-sorbing dissolved tracer within the same medium: the average migration velocity of the free colloids is higher than that of the tracer; and that the ratio of the equivalent hydrodynamic dispersion rates of colloids and tracer is dependent only upon properties of the colloids and the porous medium, it is independent of pathlengths and fluid flux, once length scales are large enough.The first of these is well known, since even in simple flow paths free colloids must stay more centre stream. The second, if validated suggests how solute and colloid dispersion may be dealt with consistently in macroscopic migration models. This is crucial since dispersion is usually ill-characterized and unaddressed by the experimental literature. In this paper we present evidence based upon an existing Drigg field injection test for the validity of these predictions.We show that starting from experimental data the fitted dispersion rates of both colloids and non-sorbing tracers increase with the measured elution rates (obeying slightly different rules for tracers and colloids); and that the ratio of colloid and nonsorbing tracer elution rates, and the ratio of colloid and nonsorbing tracer dispersion rates may be dependent upon properties of the colloids and the medium (not the flow regime).It is important to realize that even for unretarded species, an earlier peak in the breakthrough curve does not necessarily correspond to a faster mean elution rate, or vice versa. But rather that a colloid may elute faster but disperse less than an equivalent tracer. Hence its peak may be retarded compared to that of the tracer, even assuming no retardation. Hence one must consider a combination of mean elution rate and mean dispersion rate, and not rely on “peak times” to corroborate chromatographic effects. The importance of this lies in the fact that these processes are not independent and yet upscale differently. Thus realistic estimates of effective colloid dispersion rates should be upscaled in a way consistent with that adopted for tracers within the same system.     

Mesoscale influence on long-range transport — evidence from ETEX modelling and observations

During the first European Tracer Experiment (ETEX) tracer gas was released from a site in Brittany, France, and subsequently observed over a range of 2000 km. Hourly measurements were taken at the National Environmental Research Institute (NERI) located at Risø, Denmark, using two measurement techniques. At this location, the observed concentration time series shows a double-peak structure occurring between two and three days after the release. By using the Danish Emergency Response Model of the Atmosphere (DERMA), which is developed at the Danish Meteorological Institute (DMI), simulations of the dispersion of the tracer gas have been performed. Using numerical weather-prediction data from the European Centre for Medium-Range Weather Forecast (ECMWF) by DERMA, the arrival time of the tracer is quite well predicted, so also is the duration of the passage of the plume, but the double-peak structure is not reproduced. However, using higher-resolution data from the DMI version of the HIgh Resolution Limited Area Model (DMI-HIRLAM), DERMA reproduces the observed structure very well. The double-peak structure is caused by the influence of a mesoscale anti-cyclonic eddy on the tracer gas plume about one day earlier.     

Transport of coliphage PRD1 in a surface flow constructed wetland.

A tracer study was conducted in a 3-ha surface flow constructed wetland to analyze transport performance of PRD1, an enteric virus model. The convection-dispersion equation (CDE), including a first-order reaction model, adequately simulated transport performance of PRD1 in the wetland under an average hydraulic loading rate of 82 mm/d. Convective velocity (v) and longitudinal dispersion coefficient (D) were estimated by modeling a conservative tracer (bromide) pulse through the wetland. Both PRD1 and bromide were simultaneously added to the entering secondary treated wastewater effluent. The mass of bromide and PRD1 recovered was 76 and 16%, respectively. The PRD1 decay rate was calculated to be 0.3/day. The findings of this study suggest that the CDE model and analytical moment equations represent a suitable option to characterize virus transport performance in surface flow constructed wetlands.     

多沙大河的扩散示踪研究

污染物在水体中的扩散规律研究已引起越来越多人的关注。通过现场扩散示踪实验求得扩散系数，进而对污染物的扩散规律予以定量描述，是一种可靠易行的方法。对黄河托克托至河曲段的实际研究表明，以Ｃｒ２Ｏ或Ｆ－为示踪剂，运用Ｔａｙｌｏｒ模型进行多沙大河的污染物扩散规律研究可以取得满意的结果。同时，利用现场示踪实验结果与经验公式的结合可以预测不同情况下的纵向扩散系数，从而得出污染物的允许排放量，这将有利于污染物在水体中行为的深入研究。     

Characterization of organic aerosol in Big Bend National Park,Texas

The Big Bend Regional Aerosol and Visibility Observational (BRAVO) Study was conducted in Big Bend National Park, Texas, July through October 1999. Daily PM_2.5 organic aerosol samples were collected on pre-fired quartz fiber filters. Daily concentrations were too low for detailed organic analysis by gas chromatography-mass spectrometry (GC-MS) and were grouped based on their air mass trajectories. A total of 12 composites, each containing 3–10 daily samples, were analyzed. Alkane carbon preference indices suggest primary biogenic emissions were small contributors to primary PM_2.5 organic matter (OM) during the first 3 months, while in October air masses advecting from the north and south were more strongly influenced by biogenic sources. A series of trace organic compounds previously shown to serve as particle phase tracers for various carbonaceous aerosol source types were examined. Molecular tracer species were generally at or below detection limits, except for the wood smoke tracer levoglucosan in one composite, so maximum possible source influences were calculated using the detection limit as an upper bound to the tracer concentration. Wood smoke was found not to contribute significantly to PM_2.5 OM, with contributions for most samples at <1% of the total organic particulate matter. Vehicular exhaust also appeared to make only minor contributions, with maximum possible influences calculated to be 1–4% of PM_2.5 OM. Several factors indicate that secondary organic aerosol formation was important throughout the study, and may have significantly altered the molecular composition of the aerosol during transport.     

Sensitivity of the DERMA long-range Gaussian dispersion model to meteorological input and diffusion parameters

The Danish Emergency Response Model of the Atmosphere (DERMA) is described and applied to the first ETEX experiment. By using analysed low-resolution numerical weather-prediction data from the global model of the European Centre for Medium-range Weather Forecast (ECMWF) as well as higher-resolution data from two versions of the High Resolution Limited Area Model (HIRLAM), which are operational at the Danish Meteorological Institute (DMI), the sensitivity of DERMA to the resolution of meteorological data is analysed by comparing DERMA results with concentration measurements. Furthermore, the sensitivity to boundary-layer height and diffusion parameters is studied. These parameters include the critical bulk Richardson number, which is used to estimate the atmospheric boundary-layer height, the horizontal eddy diffusivity and the Lagrangian turbulence time scale. The parameters, which provide the best performance of DERMA, are 0.25 for the critical bulk Richardson number, 6×10³ m² s^-1 for the horizontal eddy diffusivity, and 3 h for the Lagrangian time scale. DERMA is much more sensitive to boundary-layer parameters when using high-resolution DMI-HIRLAM data than when using data of lower resolution from the ECMWF. Finally, the bulk Richardson number method of boundary-layer height calculation applied to DMI-HIRLAM data is verified directly against routine radiosondes released under the tracer gas plume. The boundary-layer height estimates based on analysed NWP model data agree well with observations, and the agreement deteriorates as a function of forecast length.     

Modeling the ETEX plume dispersion with the Canadian emergency response model

The CANadian Emergency Response Model (CANERM) was used to simulate the dispersion resulting from the ETEX release of 23 October 1994. Dispersion simulations were done using three different data sets as meteorological input: the ECMWF/ETEX Data Set, data from the CMC Global Data Assimilation System, and results from a diagnostic execution of the Global Environmental Multiscale (GEM) model. Comparisons of the dispersion simulations are made with observed surface concentration data provided by the Joint Research Centre (JRC) of the European Commission. It is found that CANERM can simulate fairly well the main features of the cloud dispersion. The spatial and temporal evolution of the simulated cloud appear quite plausible, but a tendency to overestimate surface concentrations is apparent. The simulations provide a credible explanation for the two peaks observed at station NL01; the first peak appears to be associated with the passage of the head portion of the plume, while the second seems to be associated with the tail part. Verification scores indicate that the simulations using the ECMWF/ETEX data set and CMC global data are of equivalent quality. However, the simulations obtained using the GEM diagnostic fields are significantly better.     

Online update of model state and parameters of a Monte Carlo atmospheric dispersion model by using ensemble Kalman filter

For an atmospheric dispersion model designed for the assessment of nuclear accident consequences, some uncertain model parameters, such as source term and weather conditions, may influence the reliability of model predictions. In this respect, good estimations of both model state and uncertain parameters are required. In this paper, an ensemble Kalman filter (EnKF) based method for simultaneous state and parameter estimation, using off-site radiation monitoring data, is presented. This method is based on a stochastic state space model, which resembles the parameter errors with stochastic quantities. Three imperfect parameters, including the source release rate, wind direction and turbulence intensity were perturbed simultaneously, and multiple parameter estimation were performed. Having been tested against both simulated and real radiation monitoring data, the method was found to be able to realistically reconstruct the real scene of dispersion, as well as the uncertain parameters. The estimated parameters given by EnKF nicely converge to the true values, and the method also tracks the temporal variation of those parameters.