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.     

Three long-range transport models compared to the ETEX experiment: a performance study

For operational or research purposes (dispersion computations of radioactive effluents during nuclear emergency situations, simulations of chemical pollution in the vicinity of thermal power plants), different models of passive dispersion in the atmosphere have been developed at the Environment Department of EDF’s R and D Division. This report presents the comparison of the performances of three such models: DIFTRA (lagrangian puff model, with operational goal), DIFEUL (three dimensional eulerian) and DIFPAR (Monte Carlo particle model) for the simulation of the first ETEX release, an international tracer campaign during which a passive tracer cloud has been followed over Europe. The results obtained in this study give model vs. experience differences of the same order as the model vs. experience differences observed during an international model comparison experiment using data of the Chernobyl release, the ATMES exercise. In addition to the standard statistical scores used in the evaluation of the performances of the transport models two asymmetric scores (in contradistinction with the Figure of Merit in Space) are proposed: “efficiency” and “power”. Their aim is to separate the two manners in which a model may be wrong: by predicting presence of pollutant while none is measured or conversely predicting absence when pollutant is actually detected.     

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.     

Chemistry and aerosols in the marine boundary layer: 1-D modelling of the three ACE-2 Lagrangian experiments

Three Lagrangian experiments were conducted during IGAC's second aerosol characterization experiment (ACE-2) in the area between Portugal, Tenerife and Madeira in June/July 1997. During each Lagrangian experiment, a boundary layer air mass was followed for about 30 h, and the temporal evolution of its chemical and aerosol composition was documented by a series of vertical profiles and horizontal box pattern flown by the Meteorological Research Flight research aircraft Hercules C130. The wealth of observational data that has been collected during these three Lagrangian experiments is the basis for the development and testing of a one-dimensional Lagrangian boundary layer model with coupled gas, aqueous, and aerosol phase chemistry. The focus of this paper is on current model limitations and strengths. We show that the model is able to represent the dynamical and chemical evolution of the marine boundary layer, in some cases requiring adjustments of the subsidence velocity and of the surface heat fluxes. Entrainment of a layer rich in ozone and carbon monoxide from a residual continental boundary layer into the marine boundary layer as well as in-cloud oxidation of sulphur dioxide by hydrogen peroxide are simulated, and coherent results are obtained, concerning the evolution of the small, presumably sulphate–ammonia aerosol mode.     

Sensitivity of the solution of the Elder problem to density, velocity and numerical perturbations

In this paper the Elder problem is studied with the purpose of evaluating the inherent instabilities associated with the numerical solution of this problem. Our focus is first on the question of the existence of a unique numerical solution for this problem, and second on the grid density and fluid density requirements necessary for a unique numerical solution. In particular we have investigated the instability issues associated with the numerical solution of the Elder problem from the following perspectives: (i) physical instability issues associated with density differences; (ii) sensitivity of the numerical solution to idealization irregularities; and, (iii) the importance of a precise velocity field calculation and the association of this process with the grid density levels that is necessary to solve the Elder problem accurately. In the study discussed here we have used a finite element Galerkin model we have developed for solving density-dependent flow and transport problems, which will be identified as TechFlow. In our study, the numerical results of Frolkovic and de Schepper [Frolkovic, P. and H. de Schepper, 2001. Numerical modeling of convection dominated transport coupled with density-driven flow in porous media, Adv. Water Resour., 24, 63-72.] were replicated using the grid density employed in their work. We were also successful in duplicating the same result with a less dense grid but with more computational effort based on a global velocity estimation process we have adopted. Our results indicate that the global velocity estimation approach recommended by Yeh [Yeh, G.-T., 1981. On the computation of Darcian velocity and mass balance in finite element modelling of groundwater flow, Water Resour. Res., 17(5), 1529-1534.] allows the use of less dense grids while obtaining the same accuracy that can be achieved with denser grids. We have also observed that the regularity of the elements in the discretization of the solution domain does make a difference in obtaining a unique stationary solution for this problem. The results of our study also indicate that the density differences are critical in the solution of the Elder problem and that high density differences lead to the physical instability that is inherent with this problem. Other than the physical instability associated with the level of density differences used in the Elder problem, the following two points should be considered in solving the Elder problem in a consistent manner: (i) strict attention should be paid to the vertical grid Peclet number in developing the criteria for convergent grid selection; and, (ii) with a globally continuous velocity calculation stable solutions can be obtained at lower grid densities.     

The kernel method to compute the intensity of segregation for reactive pollutants: Mathematical formulation

It is well known that turbulent dispersion influences chemical reactions and that computation of reactant concentrations or mean chemical reaction rates can suffer of serious error when small-scale atmospheric processes' effects on chemical transformation are neglected. A quantity that gives a measure of the influence of turbulent dispersion on second-order chemical reaction rates is the intensity of segregation. A nonparametric estimator based on the kernel method aimed at measuring the intensity of segregation is proposed. Numerical benchmark tests, in the case of a Gaussian plume, are performed to study the suitability of this technique. The estimator works well, especially for small and moderate separation from the plume centreline and generally in the smooth parts of the estimated function. The effective reaction rate is computed and the percentage error emerges to be less than 5% in the best estimation intervals, and less than 40% in the worst. A method to reduce percentage error is introduced and improved performances are observed. The estimator proposed turns out to be particularly suitable for Lagrangian air quality modelling because it permits conservation of the grid independence.     

Quantifying the effects of fumarate on in situ reductive dechlorination rates

In situ methods are needed to evaluate the effectiveness of chemical amendments at enhancing reductive dechlorination rates in groundwater that is contaminated with the priority pollutant, trichloroethene (TCE). In this communication, a method that utilizes single-well, “push–pull” tests to quantify the effects of chemical amendments on in situ reductive dechlorination rates is presented and demonstrated. Five push–pull tests were conducted in each of five monitoring wells located in a TCE-contaminated aquifer at the site of a former chemical manufacturing facility. Rates for the reductive dechlorination of the fluorinated TCE-surrogate, trichlorofluoroethene (TCFE), were measured before (test 1) and after (test 5) three successive additions (tests 2–4) of fumarate. Fumarate was selected to stimulate the growth and activity of indigenous microorganisms with the metabolic capability to reduce TCFE and TCE. In three wells, first-order rate constants for the reductive dechlorination of TCFE increased by 8.2–92 times following fumarate additions. In two wells, reductive dechlorination of TCFE was observed after fumarate additions but not before. The transformation behavior of fumarate was also monitored following each fumarate addition. Correlations between the reductive dechlorination of TCFE and the reduction of fumarate to succinate were observed, indicating that these reactions were supported by similar biogeochemical conditions at this site.     

A tool for estimating the contribution of water-soluble organic compounds to the particle mass and condensational growth in the atmosphere

In this theoretical investigation, we elucidated several factors governing the ability of organic, water-soluble vapours to participate into either the secondary organic aerosol formation or particle condensational growth in the atmosphere. The corresponding requirements for physico-chemical properties of the vapour were estimated. These estimates were also compared with the properties of several difunctional carboxylic acids identified in the atmosphere. Our analysis indicates that while many of the considered compounds are expected to contribute to the total particle mass in the atmosphere, their role in particle condensational growth process remains uncertain. This uncertainty emerges from the fact that the saturation ratio of an organic vapour does not alone determine its ability to act as a “nonvolatile” compound. Instead, the hygroscopic and chemical properties of the vapour, as well as the particle pH and size, have also to be taken into account.     

Joint propagation of variability and imprecision in assessing the risk of groundwater contamination

Estimating risks of groundwater contamination often require schemes for representing and propagating uncertainties relative to model input parameters. The most popular method is the Monte Carlo method whereby cumulative probability distributions are randomly sampled in an iterative fashion. The shortcoming of the approach, however, arises when probability distributions are arbitrarily selected in situations where available information is incomplete or imprecise. In such situations, alternative modes of information representation can be used, for example the nested intervals known as “possibility distributions”. In practical situations of groundwater risk assessment, it is common that certain model parameters may be represented by single probability distributions (representing variability) because there are data to justify these distributions, while others are more faithfully represented by possibility distributions (representing imprecision) due to the partial nature of available information. This paper applies two recent methods, designed for the joint-propagation of variability and imprecision, to a groundwater contamination risk assessment. Results of the joint-propagation methods are compared to those obtained using both interval analysis and the Monte Carlo method with a hypothesis of stochastic independence between model parameters. The two joint-propagation methods provide results in the form of families of cumulative distributions of the probability of exceeding a certain value of groundwater concentration. These families are delimited by an upper cumulative distribution and a lower distribution respectively called Plausibility and Belief after evidence theory. Slight differences between the results of the two joint-propagation methods are explained by the different assumptions regarding parameter dependencies. Results highlight the point that non-conservative results may be obtained if single cumulative probability distributions are arbitrarily selected for model parameters in the face of imprecise information and the Monte Carlo method is used under the assumption of stochastic independence. The proposed joint-propagation methods provide upper and lower bounds for the probability of exceeding a tolerance threshold. As this may seem impractical in a risk-management context, it is proposed to introduce “a-posteriori subjectivity” (as opposed to the “a-priori subjectivity” introduced by the arbitrary selection of single probability distributions) by defining a single indicator of evidence as a weighted average of Plausibility and Belief, with weights to be defined according to the specific context.     

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.     

Comparison of solubility databases used for HLW performance assessment

Elemental solubility limits are employed in many performance assessment studies of the disposal of high-level radioactive waste. A comparison of solubility databases used in a number of recent assessments from Canada, Finland, Japan, Sweden, Switzerland and the U.S.A. is presented. Despite similarities in the chemical conditions under which such solubilities are defined, there are very large ranges in the values selected for most elements (several orders of magnitude). In some cases, differences in solubilities can be attributed to differences in the thermodynamic databases used but, due to generally incomplete documentation, both of such databases and the criteria for data selection, detailed analysis of the reason for such discrepancies is precluded. A key factor, however, is the fairly subjective choice of the solid phase considered to be “solubility-limiting”. Procedures are recommended which would increase the transparency of solubility database compilation in future analyses.     

Regional NOx emission inversion through a four-dimensional variational approach using SCIAMACHY tropospheric NO2 column observations

In this paper, the NOx emission scaling factors applied over the 2001 National Emissions Inventory (NEI) are estimated through a four-dimensional variational (4D-Var) approach using SCIAMACHY (Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY) tropospheric NO₂ columns measured during summer 2004. In the “top-down” approach, two-month average NO₂ columns are assimilated into a regional chemical transport model (CTM), STEM, using different assimilation setups. In a basic setup, NOx emissions are adjusted by assimilating the NO₂ columns. A more general setup of emission inversion allows the initial O₃ concentrations be adjusted along with the NOx emissions. A final case is set up to assimilate both the NO₂ columns and O₃ measurement from various platforms while allowing adjustments of both the NOx emissions and the initial O₃ concentrations. It is found that the addition of O₃ measurements did not improve the NOx emission inversion. With the NOx emission at surface and upper levels being adjusted separately, results from four cases show that the elevated NOx emission reduction ranges from 8.9% to 11.4%, and the surface NOx emission reduction is up to 6.6%. All the cases show NOx emission reduction in Ohio valley and Washington, District of Columbia areas.     

Chemical leasing: cooperative business models for sustainable chemicals management. Summary of research projects commissioned by the Austrian Federal Ministry of Agriculture, Forestry, Environment and Water Management

- Sustainable chemistry - Section editors: Klaus Günter Steinhäuser, Steffi Richter, Petra Greiner, Jutta Penning, Michael AngrickBackground, Aim and Scope Chemicals play a vital role in the day-to-day life of industrialised societies. Their use is not restricted to the chemical enterprises per se, but is a crucial part of production processes in a lot of industrial sectors. Traditional instruments of environmental policy (such as bans, restrictions) can only deal with the most hazardous substances. The Johannesburg Implementation Plan of 2002 calls for more sustainable patterns of production and consumption, and sets the year of 2020 as a goal to use chemicals in a way that human health and the environment are not endangered. Political instruments should not only gather more knowledge about the properties of chemicals, but should also stimulate the environmentally sound use of chemicals. Existing business models should therefore be reviewed in relation to this strategic approach to encourage marketing options with respect to the environmental focus.Main Features Business models were examined for their effects on the consumption of chemicals and amount of waste emissions in relation to their economic potential. Different possibilities for cooperation of supplier, user and disposal companies were elaborated and examined with a view to the specific situation in Austria.Results and Discussion A range of cooperative models – summarised under the term ‘chemical leasing’ - was identified, which can contribute to a more efficient use of resources. 12 main possible application areas (cleaning, lubrication, paint stripping and others) have been identified in Austria. If chemical leasing models were applied in these areas, the amounts of chemicals currently used could be reduced by one third (53,000 tonnes per year). Cost reductions of up to 15 % can be expected.Conclusion The application of chemical leasing models can contribute considerably to achieving more sustainable and resource-efficient patterns of production. The Austrian Ministry for Environment has therefore decided to subsidise the further practical implementation of these new service-oriented business models. Pilot projects in 4 enterprises, which are supervised by consulting companies, are currently being carried out. Recommendation and Outlook The experiences of the pilot projects will serve as valuable building blocks for the wider use of chemical leasing models. Furthermore, the UNIDO Cleaner Production Centres have expressed their clear interest and will examine the possibility to use chemical leasing as a part of their Cleaner Production Programmes.     

Modelling tritium and phosphorus transport by preferential flow in structured soil

Subsurface solute transport through structured soil is studied by model interpretation of experimental breakthrough curves from tritium and phosphorus tracer tests in three intact soil monoliths. Similar geochemical conditions, with nearly neutral pH, were maintained in all the experiments. Observed transport differences for the same tracer are thus mainly due to differences in the physical transport process between the different monoliths. The modelling is based on a probabilistic Lagrangian approach that decouples physical and chemical mass transfer and transformation processes from pure and stochastic advection. Thereby, it enables explicit quantification of the physical transport process through preferential flow paths, honouring all independently available experimental information. Modelling of the tritium breakthrough curves yields a probability density function of non-reactive solute travel time that is coupled with a reaction model for linear, non-equilibrium sorption–desorption to describe the phosphorus transport. The tritium model results indicate that significant preferential flow occurs in all the experimental soil monoliths, ranging from 60–100% of the total water flow moving through only 25–40% of the total water content. In agreement with the fact that geochemical conditions were similar in all experiments, phosphorus model results yield consistent first-order kinetic parameter values for the sorption–desorption process in two of the three soil monoliths; phosphorus transport through the third monolith cannot be modelled because the apparent mean transport rate of phosphorus is anomalously rapid relative to the non-adsorptive tritium transport. The occurrence of preferential flow alters the whole shape of the phosphorus breakthrough curve, not least the peak mass flux and concentration values, and increases the transported phosphorus mass by 2–3 times relative to the estimated mass transport without preferential flow in the two modelled monoliths.     

Laboratory models to evaluate phytotoxicity of sulphadimethoxine on terrestrial plants

Antimicrobial drugs soil contamination is a possible side effect of spreading animal wastes from intensive farming on arable land: waste often contain persistent drugs (like sulphonamides) that pollute soil. In laboratory tests different crop plants and weeds demonstrated, on both synthetic medium and soil, toxicity and sulphadimethoxine bioaccumulation. These data suggest potential adverse implications for organisms higher up the food net, including humans. The increasingly used technique of'soil top dressing will impose the monitoring of antimicrobial contamination in arable lands and their crops, to protect natural ecosystems and consumer. A possible tool for this monitoring is the use of weeds that can constitute a “mesh” where antimicrobial contamination can be detected.     

Levels, profile and pattern of PCDDs and PCDFs in samples related to the production and use of chlorine

Analyses of sludge from graphite electrodes used in the chloralkali process show total levels of PCDFs as high as 650000 pg/g sludge. The levels of tetra-, penta-and hexaCDFs were found to be approximately the same. The levels for the corresponding PCDDs were below the detection level. The dominating congeners within each group are the toxic 2,3,7,8-substituted PCDFs at levels of 340000 pg/g, This typical pattern, called the “chloralkali pattern” can also be found in soil samples taken at another chloralkali plant. The typical pattern can also be identified in a solution of ferric chloride. Traces of “chloralkali pattern” are identified in a sludge sample from the drinking water purification plant, the source for PCDFs in the sludge is unknown.     

Application of continuous time random walk theory to nonequilibrium transport in soil

Continuous time random walk (CTRW) formulations have been demonstrated to provide a general and effective approach that quantifies the behavior of solute transport in heterogeneous media in field, laboratory, and numerical experiments. In this paper we first apply the CTRW approach to describe the sorbing solute transport in soils under chemical (or) and physical nonequilibrium conditions by curve-fitting. Results show that the theoretical solutions are in a good agreement with the experimental measurements. In case that CTRW parameters cannot be determined directly or easily, an alternative method is then proposed for estimating such parameters independently of the breakthrough curve data to be simulated. We conduct numerical experiments with artificial data sets generated by the HYDRUS-1D model for a wide range of pore water velocities (υ) and retardation factors (R) to investigate the relationship between CTRW parameters for a sorbing solute and these two quantities (υ, R) that can be directly measured in independent experiments. A series of best-fitting regression equations are then developed from the artificial data sets, which can be easily used as an estimation or prediction model to assess the transport of sorbing solutes under steady flow conditions through soil. Several literature data sets of pesticides are used to validate these relationships. The results show reasonable performance in most cases, thus indicating that our method could provide an alternative way to effectively predict sorbing solute transport in soils. While the regression relationships presented are obtained under certain flow and sorption conditions, the methodology of our study is general and may be extended to predict solute transport in soils under different flow and sorption conditions.     

Levels of PCDDs and PCDFs in products and effluent from the Swedish pulp and paper industry and chloralkali process

Within the Swedish Dioxin Survey various samples from the pulp and paper industry and the chloralkali process have been analyzed by congener specific analytical methods. In addition to the generally discussed “bleaching pattern” of the tetrachlorinated congeners, these samples also contained higher chlorinated congeners like hexa-CDDs, hepta-CDFs, octa-CDD and octa-CDF. Consequently it is recommended that samples from the pulp industry should be analyzed for all PCDDs and PCDFs (tetra- through octa-). Counted as Nordic Toxic Equivalents (NTEQ), the recycled pulp samples had the highest contamination level followed by TMP, unbleached sulfite and bleached softwood and hardwood. In addition to the bleaching process, various chemicals used in the pulping, bleaching and wastewater treatment can contribute to the contamination.     

A hybrid modeling approach to resolve pollutant concentrations in an urban area

A modeling tool that can resolve contributions from individual sources to the urban environment is critical for air-toxics exposure assessments. Air toxics are often chemically reactive and may have background concentrations originated from distant sources. Grid models are the best-suited tools to handle the regional features of these chemicals. However, these models are not designed to resolve pollutant concentrations on local scales. Moreover, for many species of interest, having reaction time scales that are longer than the travel time across an urban area, chemical reactions can be ignored in describing local dispersion from strong individual sources making Lagrangian and plume-dispersion models practical. In this study, we test the feasibility of developing an urban hybrid simulation system. In this combination, the Community Multi-scale Air Quality model (CMAQ) provides the regional background concentrations and urban-scale photochemistry, and local models such as Hybrid Single Particle Lagrangian Integrated Trajectory model (HYSPLIT) and AMS/EPA Regulatory Model (AERMOD) provide the more spatially resolved concentrations due to local emission sources. In the initial application, the HYSPLIT, AERMOD, and CMAQ models are used in combination to calculate high-resolution benzene concentrations in the Houston area. The study period is from 18 August to 4 September of 2000. The Mesoscale Model 5 (MM5) is used to create meteorological fields with a horizontal resolution of 1×1 km². In another variation to this approach, multiple HYSPLIT simulations are used to create a concentration ensemble to estimate the contribution to the concentration variability from point sources. HYSPLIT simulations are used to model two sources of concentration variability; one due to variability created by different particle trajectory pathways in the turbulent atmosphere and the other due to different flow regimes that might be introduced when using gridded data to represent meteorological data fields. The ensemble mean concentrations determined by HYSPLIT plus the concentrations estimated by AERMOD are added to the CMAQ calculated background to estimate the total mean benzene concentration. These estimated hourly mean concentrations are also compared with available field measurements.