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.     

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.     

Comparison of a new Eulerian model with a modified Lagrangian approach for particle distribution and deposition indoors

Understanding of aerosol dispersion characteristics has many scientific and engineering applications. It is recognized that Eulerian or Lagrangian approach has its own merits and limitations. A new Eulerian model has been developed and it adopts a simplified drift–flux methodology in which external forces can be incorporated straightforwardly. A new near-wall treatment is applied to take into account the anisotropic turbulence for the modified Lagrangian model. In the present work, we present and compare both Eulerian and Lagrangian models to simulate particle dispersion in a small chamber. Results reveal that the standard k–ε Lagrangian model over-predicts particle deposition compared to the present turbulence-corrected Lagrangian approach. Prediction by the Eulerian model agrees well with the modified Lagrangian model.     

Global backtracking of anthropogenic radionuclides by means of a receptor oriented ensemble dispersion modelling system in support of Nuclear-Test-Ban Treaty verification

In this paper, we introduce a methodology for quality assessment of backtracking models. We present results illustrating the level of agreement between the backtracking models, and the accuracy of each model and the ensemble model in resolving the geo-temporal reference of a single point source. Both assessments are based on an ensemble of 12 different Lagrangian particle dispersion modelling (LPDM) systems utilized in receptor oriented (adjoint) mode during an international numerical experiment dedicated to source region estimation.As major result, we can confirm that the findings of Galmarini et al. [2004b. Ensemble prediction forecasting—Part II: application and evaluation. Atmospheric Environment 38, 4619–4632] and Delle Monache and Stull [2003. An ensemble air-quality forecast over Europe during an ozone episode. Atmospheric Environment 37, 3469–3474], regarding the superiority of the ensemble dispersion forecast over a single forecast, do also apply to LPDM when utilized for backtracking purposes, in particular if only vague a priori knowledge of the source time is available. This, however, is a likely situation in the context of the global nuclear monitoring performed by the Provisional Technical Secretariat (PTS) of the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO), where quick but reliable source location identification is required. We introduce a simple methodology as a template for a future electronic emergency response system in the field of dispersion modelling.     

Assimilation of conventional and satellite wind observations in a mesoscale atmospheric model for studying atmospheric dispersion

A mesoscale atmospheric model PSU/NCAR MM5 is used to provide operational weather forecasts for a nuclear emergency response decision support system on the southeast coast of India. In this study the performance of the MM5 model with assimilation of conventional surface and upper-air observations along with satellite derived 2-d surface wind data from QuickSCAT sources is examined. Two numerical experiments with MM5 are conducted: one with static initialization using NCEP FNL data and second with dynamic initialization by assimilation of observations using four dimensional data assimilation (FDDA) analysis nudging for a pre-forecast period of 12 h. Dispersion simulations are conducted for a hypothetical source at Kalpakkam location with the HYSPLIT Lagrangian particle model using simulated wind field from the above experiments. The present paper brings out the differences in the atmospheric model predictions and the differences in dispersion model results from control and assimilation runs. An improvement is noted in the atmospheric fields from the assimilation experiment which has led to significant alteration in the trajectory positions, plume orientation and its distribution pattern. Sensitivity tests using different PBL and surface parameterizations indicated the simple first order closure schemes (Blackadar, MRF) coupled with the simple soil model have given better results for various atmospheric fields. The study illustrates the impact of the assimilation of the scatterometer wind and automated weather stations (AWS) observations on the meteorological model predictions and the dispersion results.     

An efficient Lagrangian stochastic model of vertical dispersion in the convective boundary layer

We consider the one-dimensional case of vertical dispersion in the convective boundary layer (CBL) assuming that the turbulence field is stationary and horizontally homogeneous. The dispersion process is simulated by following Lagrangian trajectories of many independent tracer particles in the turbulent flow field, leading to a prediction of the mean concentration. The particle acceleration is determined using a stochastic differential equation, assuming that the joint evolution of the particle velocity and position is a Markov process. The equation consists of a deterministic term and a random term. While the formulation is standard, attention has been focused in recent years on various ways of calculating the deterministic term using the well-mixed condition incorporating the Fokker–Planck equation. Here we propose a simple parameterisation for the deterministic acceleration term by approximating it as a quadratic function of velocity. Such a function is shown to represent well the acceleration under moderate velocity skewness conditions observed in the CBL. The coefficients in the quadratic form are determined in terms of given turbulence statistics by directly integrating the Fokker–Planck equation. An advantage of this approach is that, unlike in existing Lagrangian stochastic models for the CBL, the use of the turbulence statistics up to the fourth order can be made without assuming any predefined form for the probability distribution function (PDF) of the velocity. The main strength of the model, however, lies in its simplicity and computational efficiency. The dispersion results obtained from the new model are compared with existing laboratory data as well as with those obtained from a more complex Lagrangian model in which the deterministic acceleration term is based on a bi-Gaussian velocity PDF. The comparison shows that the new model performs well.     

Validation of a Lagrangian dispersion model implementing different kernel methods for density reconstruction

In this paper the inert version of a Lagrangian particle model named photochemical Lagrangian particle model (PLPM) is described and validated. PLPM implements four density reconstruction algorithms based on the kernel density estimator. All these methods are fully grid-free but they differ each other in considering local or global features of the particles distribution, in treating the Cartesian directions separately or together and in being based on receptors or particles positions in space. Each kernel has been shown to have both advantages and disadvantages, but the overall good performances of the model when compared with the well known Copenhagen and Kincaid data sets are very encouraging in view of its extension to fully chemically active simulations, currently under development.     

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.     

An atmospheric dispersion model for the environmental impact assessment of thermal power plants in Japan--a method for evaluating topographical effects

An atmospheric dispersion model was developed for the environmental impact assessment of thermal power plants in Japan, and a method for evaluating topographical effects using this model was proposed. The atmospheric dispersion model consists of an airflow model with a turbulence closure model based on the algebraic Reynolds stress model and a Lagrangian particle dispersion model (LPDM). The evaluation of the maximum concentration of air pollutants such as SO2, NOx, and suspended particulate matter is usually considered of primary importance for environmental impact assessment. Three indices were therefore estimated by the atmospheric dispersion model: the ratios (alpha and beta, respectively) of the maximum concentration and the distance of the point of the maximum concentration from the source over topography to the respective values over a flat plane, and the relative concentration distribution [gamma(x)] along the ground surface projection of the plume axis normalized by the maximum concentration over a flat plane. The atmospheric dispersion model was applied to the topography around a power plant with a maximum elevation of more than 1,000 m. The values of alpha and beta evaluated by the atmospheric dispersion model varied between 1 and 3 and between 1 and 0.4, respectively, depending on the topographical features. These results and the calculated distributions of y(x) were highly similar to the results of the wind tunnel experiment. Therefore, when the slope of a hill or mountain is similar to the topography considered in this study, it is possible to evaluate topographical effects on exhaust gas dispersion with reasonable accuracy using the atmospheric dispersion model as well as wind tunnel experiments.     

The effect of improved nowcasting of precipitation on air quality modeling

The predictive potential of air quality models and thus their value in emergency management and public health support are critically dependent on the quality of their meteorological inputs. The atmospheric flow is the primary cause of the dispersion of airborne substances. The scavenging of pollutants by cloud particles and precipitation is an important sink of atmospheric pollution and subsequently determines the spatial distribution of the deposition of pollutants. The long-standing problem of the spin-up of clouds and precipitation in numerical weather prediction models limits the accuracy of the prediction of short-range dispersion and deposition from local sources. The resulting errors in the atmospheric concentration of pollutants also affect the initial conditions for the calculation of the long-range transport of these pollutants. Customary the spin-up problem is avoided by only using NWP (Numerical Weather Prediction) forecasts with a lead time greater than the spin-up time of the model. Due to the increase of uncertainty with forecast range this reduces the quality of the associated forecasts of the atmospheric flow.In this article recent improvements through diabatic initialization in the spin-up of large-scale precipitation in the Hirlam NWP model are discussed. In a synthetic example using a puff dispersion model the effect is demonstrated of these improvements on the deposition and dispersion of pollutants with a high scavenging coefficient, such as sulphur, and a low scavenging coefficient, such as cesium-137. The analysis presented in this article leads to the conclusion that, at least for situations where large-scale precipitation dominates, the improved model has a limited spin-up so that its full forecast range can be used. The implication for dispersion modeling is that the improved model is particularly useful for short-range forecasts and the calculation of local deposition. The sensitivity of the hydrological processes to proper initialization implies that the spin-up problem may reoccur with changes in the model and increased model resolution. Spin-up should be an ongoing concern for atmospheric modelers.     

Evaluation of the Meteo-France response in ETEX release 1

During ETEX Meteo-France applied part of its emergency response system for critical events developped in the framework of the World Meteorological Organization environmental emergency response program. The atmospheric transport model used to forecast the evolution of a passive tracer is an eulerian model called MEDIA. In real time this model is driven by meteorological data from ARPEGE, the operational numerical weather prediction model available at the Meteo-France operation center. The overall evaluation of the results show that the model can reproduce the cloud displacement, but there exists a stretching in the transport direction. In the ATMES-II phase, the results are closer to the observations when meteorological data from the European Center for Medium range Weather Forecast are used. A simulation using analyzed meteorological data from ARPEGE every 6 h slightly improve the results comparing with the real-time experiment. All the simulations we performed reveal that the quality of the atmospheric transport model is strongly dependent on the quality of the driving numerical weather prediction model.     

A statistical approach for estimating uncertainty in dispersion modeling: An example of application in southwestern USA

A method based on a statistical approach of estimating uncertainty in simulating the transport and dispersion of atmospheric pollutants is developed using observations and modeling results from a tracer experiment in the complex terrain of the southwestern USA. The method takes into account the compensating nature of the error components by representing all terms, except dispersion error and variance of stochastic processes. Dispersion error and the variance of the stochastic error are estimated using the maximum likelihood estimation technique applied to the equation for the fractional error. Mesoscale Model 5 (MM5) and a Lagrangian random particle dispersion model with three optional turbulence parameterizations were used as a test bed for method application. Modeled concentrations compared well with the measurements (correlation coefficients on the order of 0.8). The effects of changing two structural components (the turbulence parameterization and the model grid vertical resolution) on the magnitude of the dispersion error also were examined. The expected normalized dispersion error appears to be quite large (up to a factor of three) among model runs with various turbulence schemes. Tests with increased vertical resolution of the atmospheric model (MM5) improved most of the dispersion model statistical performance measures, but to a lesser extent compared to selection of a turbulence parameterization. Method results confirm that structural components of the dispersion model, namely turbulence parameterizations, have the most influence on the expected dispersion error.     

Assimilating concentration observations for transport and dispersion modeling in a meandering wind field

Assimilating concentration data into an atmospheric transport and dispersion model can provide information to improve downwind concentration forecasts. The forecast model is typically a one-way coupled set of equations: the meteorological equations impact the concentration, but the concentration does not generally affect the meteorological field. Thus, indirect methods of using concentration data to influence the meteorological variables are required. The problem studied here involves a simple wind field forcing Gaussian dispersion. Two methods of assimilating concentration data to infer the wind direction are demonstrated. The first method is Lagrangian in nature and treats the puff as an entity using feature extraction coupled with nudging. The second method is an Eulerian field approach akin to traditional variational approaches, but minimizes the error by using a genetic algorithm (GA) to directly optimize the match between observations and predictions. Both methods show success at inferring the wind field. The GA-variational method, however, is more accurate but requires more computational time. Dynamic assimilation of a continuous release modeled by a Gaussian plume is also demonstrated using the genetic algorithm approach.     

Dispersion over and around a steep obstacle for varying thermal stratification—numerical simulations

A rather steep dump of an open-cast working area for lignite is situated close to the nuclear research installation of Juelich. This dump is more than 200 m high; under light wind conditions the air flow goes around this obstacle during night-time (stable stratification) and over it in the early morning hours when the stratification is destabilized. This air flow is simulated by the non-hydrostatic mesoscale numerical model FITNAH; for the time 0500–0900 LST the concentration field (of a point source) is simulated with a Lagrangian particle model. The computed flow behaviour causes a bifurcation of the plume during night while after 0900 LST dispersion is nearly unaffected by topography.     

Testing the importance of accurate meteorological input fields and parameterizations in atmospheric transport modelling using dream - validation against ETEX-1

A tracer model, the DREAM, which is based on a combination of a near-range Lagrangian model and a long-range Eulerian model, has been developed. The meteorological meso-scale model, MM5V1, is implemented as a meteorological driver for the tracer model. The model system is used for studying transport and dispersion of air pollutants caused by a single but strong source as, e.g. an accidental release from a nuclear power plant. The model system including the coupling of the Lagrangian model with the Eulerian model are described. Various simple and comprehensive parameterizations of the mixing height, the vertical dispersion, and different meterological input data have been implemented in the combined tracer model, and the model results have been validated against measurements from the ETEX-1 release. Several different statistical parameters have been used to estimate the differences between the parameterizations and meterological input data in order to find the best performing solution.     

Norwegian Meteorological Institute’s real-time dispersion model snap (Severe Nuclear Accident Program): Runs for ETEX and ATMES II experiments with different meteorological input

The Norwegian Meteorological Institute (DNMI) has developed and implemented for operational use a real-time dispersion model Severe Nuclear Accident Program (SNAP) with capability for predicting concentrations and depositions of the radioactive debris from large accidental releases. SNAP has been closely linked to DNMI’s operational numerical weather prediction (NWP) models.How good are these predictions? Participation in ETEX has partly answered this question. DNMI used SNAP with LAM50S giving meteorological input for these real-time dispersion calculations. LAM50S Limited Area Model with 50 km grid squareswas DNMI’s operational NWP model in 1994 when ETEX took place.In this article we report on how SNAP performed in the first of the ETEX releases in near-real-time mode, using LAM50S—and in hindcast mode for ATMES II, using “ECMWF 1995: ETEX Data set (ATMES II)”as meteorological input data. These two input data sets came from NWP models with quite different characteristics but with similar resolution in time and space.The results from these dispersion simulations matched closely. Deviations early in the simulation period shrank to insignificant differences later on. Since both input data sets were based on “weather analysis” and had similar resolution in space and time, SNAP described the dispersion of the released material very similar in these two simulations.     

An Atmospheric Dispersion Model for the Environmental Impact Assessment of Thermal Power Plants in Japan—A Method for Evaluating Topographical Effects

Abstract

An atmospheric dispersion model was developed for the environmental impact assessment of thermal power plants in Japan, and a method for evaluating topographical effects using this model was proposed. The atmospheric dispersion model consists of an airflow model with a turbulence closure model based on the algebraic Reynolds stress model and a Lagrangian particle dispersion model (LPDM). The evaluation of the maximum concentration of air pollutants such as SO₂, NO_x, and suspended particulate matter is usually considered of primary importance for environmental impact assessment. Three indices were therefore estimated by the atmospheric dispersion model: the ratios (α and β, respectively) of the maximum concentration and the distance of the point of the maximum concentration from the source over topography to the respective values over a flat plane, and the relative concentration distribution [γ(x)] along the ground surface projection of the plume axis normalized by the maximum concentration over a flat plane. The atmospheric dispersion model was applied to the topography around a power plant with a maximum elevation of more than 1000 m. The values of α and β evaluated by the atmospheric dispersion model varied between 1 and 3 and between 1 and 0.4, respectively, depending on the topographical features. These results and the calculated distributions of γ(x) were highly similar to the results of the wind tunnel experiment. Therefore, when the slope of a hill or mountain is similar to the topography considered in this study, it is possible to evaluate topographical effects on exhaust gas dispersion with reasonable accuracy using the atmospheric dispersion model as well as wind tunnel experiments.     

Development of the ClearSky smoke dispersion forecast system for agricultural field burning in the Pacific Northwest

The post-harvest burning of agricultural fields is commonly used to dispose of crop residue and provide other desired services such as pest control. Despite careful regulation of burning, smoke plumes from field burning in the Pacific Northwest commonly degrade air quality, particularly for rural populations. In this paper, ClearSky, a numerical smoke dispersion forecast system for agricultural field burning that was developed to support smoke management in the Inland Pacific Northwest, is described. ClearSky began operation during the summer through fall burn season of 2002 and continues to the present. ClearSky utilizes Mesoscale Meteorological Model version 5 (MM5v3) forecasts from the University of Washington, data on agricultural fields, a web-based user interface for defining burn scenarios, the Lagrangian CALPUFF dispersion model and web-served animations of plume forecasts. The ClearSky system employs a unique hybrid source configuration, which treats the flaming portion of a field as a buoyant line source and the smoldering portion of the field as a buoyant area source. Limited field observations show that this hybrid approach yields reasonable plume rise estimates using source parameters derived from recent field burning emission field studies. The performance of this modeling system was evaluated for 2003 by comparing forecast meteorology against meteorological observations, and comparing model-predicted hourly averaged PM_2.5 concentrations against observations. Examples from this evaluation illustrate that while the ClearSky system can accurately predict PM_2.5 surface concentrations due to field burning, the overall model performance depends strongly on meteorological forecast error. Statistical evaluation of the meteorological forecast at seven surface stations indicates a strong relationship between topographical complexity near the station and absolute wind direction error with wind direction errors increasing from approximately 20° for sites in open areas to 70° or more for sites in very complex terrain. The analysis also showed some days with good forecast meteorology with absolute mean error in wind direction less than 30° when ClearSky correctly predicted PM_2.5 surface concentrations at receptors affected by field burns. On several other days with similar levels of wind direction error the model did not predict apparent plume impacts. In most of these cases, there were no reported burns in the vicinity of the monitor and, thus, it appeared that other, non-reported burns were responsible for the apparent plume impact at the monitoring site. These cases do not provide information on the performance of the model, but rather indicate that further work is needed to identify all burns and to improve burn reports in an accurate and timely manner. There were also a number of days with wind direction errors exceeding 70° when the forecast system did not correctly predict plume behavior.     

High-resolution simulation of volcanic sulfur dioxide dispersion over the Miyake Island

After severe eruptions of the volcano at Miyake Island in August 2000, a large amount of volcanic gas was released into the atmosphere. To simulate flows and dispersion of sulfur dioxide (SO₂) over Miyake Island, a set of numerical models was developed. The multi-nesting method was adopted to reflect a realistic meteorological field and to sufficiently resolve the flow over the island with a diameter of 8 km. The outermost model was the Regional Spectral Model (RSM) of the Japan Meteorological Agency (JMA) with a horizontal grid size of 10 km. Finer atmospheric structure was simulated with the nonhydrostatic model jointly developed by the Meteorological Research Institute and the Numerical Prediction Division of JMA (MRI/NPD-NHM) with grid intervals of 2 km, 400 m and 100 m. Realistic topography of the island was represented in the innermost model. The Lagrangian particle method was applied to the dispersion model, which is driven by the meteorological field of the 100 m grid MRI/NPD-NHM. The random walk procedure was used to represent the turbulent diffusion. The model was verified in four cases. Simulated SO₂ concentrations agreed well with observed concentrations at a monitoring station including temporal variation. Under a large synoptic change, however, accurate prediction became difficult. Further numerical experiments have been done to investigate characteristics of the flow and the distribution of SO₂. Steady inflows, classified according to the surface wind speed and direction, were assumed. Simulated SO₂ distribution on the ground apparently depends on the surface wind. Under relatively weak inflow, there is a large diurnal change in SO₂ distribution, affected by the thermally induced flow. SO₂ gas is widely spread downstream in the nighttime but hardly reaches the coastal area in the daytime. On the other hand, SO₂ gas steadily reached the downstream coast with little diurnal variation under the stronger inflow. Ground temperature, as well as the static stability of the inflow, also influences downstream wind, turbulent diffusivity and SO₂ distribution.     

Evaluating the use of outputs from comprehensive meteorological models in air quality modeling applications

Currently used dispersion models, such as the AMS/EPA Regulatory Model (AERMOD), process routinely available meteorological observations to construct model inputs. Thus, model estimates of concentrations depend on the availability and quality of meteorological observations, as well as the specification of surface characteristics at the observing site. We can be less reliant on these meteorological observations by using outputs from prognostic models, which are routinely run by the National Oceanic and Atmospheric Administration (NOAA). The forecast fields are available daily over a grid system that covers all of the United States. These model outputs can be readily accessed and used for dispersion applications to construct model inputs with little processing. This study examines the usefulness of these outputs through the relative performance of a dispersion model that has input requirements similar to those of AERMOD. The dispersion model was used to simulate observed tracer concentrations from a Tracer Field Study conducted in Wilmington, California in 2004 using four different sources of inputs: (1) onsite measurements; (2) National Weather Service measurements from a nearby airport; (3) readily available forecast model outputs from the Eta Model; and (4) readily available and more spatially resolved forecast model outputs from the MM5 prognostic model. The comparison of the results from these simulations indicate that comprehensive models, such as MM5 and Eta, have the potential of providing adequate meteorological inputs for currently used short-range dispersion models such as AERMOD.