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.     

Non-linear response of wet deposition to emissions reduction: A model study

An Eulerian atmospheric model with complex chemistry (Acidic Deposition and Oxidant Model) and a Lagrangian model with linear chemistry (Ontario Ministry of the Environment Trajectory Model) were used to simulate the wet SO₄²⁻ deposition pattern over eastern North America for 16 days during April 1981.The two model results agree reasonably well with each other when the 16 day average values are compared. They also show reasonable agreement with observed data. Having established the ability of the models to predict deposition patterns for 1981 emissions, reduction scenarios with 50% SO_x and 50% SO_x and NO_x of the 1981 emissions were studied through the Eulerian model. Near the heavy emissions area, the reduction in SO₄²⁻ wet deposition is only about 30–40%. In this respect the linear Lagrangian model departs significantly from the Eulerian model. This non-linearity in response is attributed to the role of oxidants in controlling the conversion of SO₂ to SO₄²⁻.     

Examination of a subgrid-scale parameterization for the transport of pollutants in a nonprecipitating cumulus cloud ensemble

The mass flux based subgrid-scale parameterization technique of Gidel (1983, J. geophys. Res.88, 6587–6599) is re-examined for use in Eulerian long-range transport models. Specifically, the parameterization scheme is incorporated into the STEM-II Eulerian transport/transformation/removal model and the model is used to investigate pollutant transport in a nonprecipitating cumulus cloud ensemble. The effects of entrainment, detrainment, evaporation and the transport by subsidence, updrafts and turbulent diffusion are included in the analysis. Presented simulation results indicate that the parameterization is able to treat the rapid vertical transport by cloud updrafts, enables the calculation of reaction rates based on subgrid-scale concentrations, and is readily adopted by Eulerian models.     

Modelling annual pasture dynamics: Application to stomatal ozone deposition

Modelling ozone (O₃) deposition for impact risk assessment is still poorly developed for herbaceous vegetation, particularly for Mediterranean annual pastures. High inter-annual climatic variability in the Mediterranean area makes it difficult to develop models characterizing gas exchange behaviour and air pollutant absorption suitable for risk assessment. This paper presents a new model to estimate stomatal conductance (g_s) of Trifolium subterraneum, a characteristic species of dehesa pastures. The MEDPAS (MEDiterranean PAStures) model couples 3 modules estimating soil water content (SWC), vegetation growth and g_s. The g_s module is a reparameterized version of the stomatal component of the EMEP DO₃SE O₃ deposition model. The MEDPAS model was applied to two contrasting years representing typical dry and humid springs respectively and with different O₃ exposures. The MEDPAS model reproduced realistically the g_s seasonal and inter-annual variations observed in the field. SWC was identified as the major driver of differences across years. Despite the higher O₃ exposure in the dry year, meteorological conditions favoured 2.1 times higher g_s and 56 day longer growing season in the humid year compared to the dry year. This resulted in higher ozone fluxes absorbed by T. subterraneum in the humid year. High inter-family variability was found in gas exchange rates, therefore limiting the relevance of single species O₃ deposition flux modelling for dehesa pastures. Stomatal conductance dynamics at the canopy level need to be considered for more accurate O₃ flux modelling for present and future climate scenarios in the Mediterranean area.     

Modeling Atmospheric Mercury Deposition in the Vicinity of Power Plants

Abstract

Two mathematical models of the atmospheric fate and transport of mercury (Hg), an Eulerian grid–based model and a Gaussian plume model, are used to calculate the atmospheric deposition of Hg in the vicinity (i.e., within 50 km) of five coal–fired power plants. The former is applied using two different horizontal resolutions: coarse (84 km) and fine (16.7 km). More than 96% of the power plant Hg emissions are calculated with the plume model to be transported beyond 50 km from the plants. The grid–based model predicts a lower fraction to be transported beyond 50 km: >91% with a coarse resolution and >95% with a fine resolution. The contribution of the power plant emissions to total Hg deposition within a radius of 50 km from the plants is calculated to be <8% with the plume model, <14% with the Eulerian model with a coarse resolution, and <10% with the Eulerian model with a fine resolution. The Eulerian grid–based model predicts greater local impacts than the plume model because of artificially enhanced vertical dispersion; the former predicts about twice as much Hg deposition as the latter when the area considered is commensurate with the resolution of the grid–based model. If one compares the local impacts for an area that is significantly less than the grid–based model resolution, then the grid–based model may predict lower local deposition than the plume model, because two compensating errors affect the results obtained with the grid–based model: initial dilution of the power plant emissions within one or more grid cells and enhanced vertical mixing to the ground.     

KAPPA-G,A Non-Gaussian Plume Dispersion Model: Description and Evaluation Against Tracer Measurements

This paper presents a mixed methodology for the simulation of atmospheric disperson phenomena in which vertical diffusion is computed using an analytical solution of the K-theory equation, while horizontal diffusion is simulated by the Gaussian formula. This new formulation, while maintaining a simple analytical form for the concentration field, incorporates the effects of power-law vertical profiles of both wind speed and eddy diffusivity. The performance of this approach, which has been implemented into a full computer package (KAPPA-G), is evaluated by comparison with data from SF₆ tracer experiments.     

Measurements and modelling of the ozone deposition velocity to concrete tiles,including the effect of diffusion

A range of models were fitted to the experimental time-dependent curves for the deposition velocity of O₃ to concrete floor tile samples. The models included modified Langmuir isotherms assuming adsorption of O₃ on the material surfaces and models assuming direct reaction on and diffusion of O₃ into the material from the air. The best fit was obtained with a simple two-parameter model assuming direct reaction of O₃ with adsorbed surface water and direct diffusion of O₃ into the material from the air. However, models assuming an additional second-order reaction of O₃ with an adsorbed surface species with a given start mass gave improved curve fit in the first 500 min. Applying the best model to experimental data obtained over the whole air humidity range resulted in markedly lower equilibrium deposition velocities than those measured after 48 h. The modelling gave a deposition velocity minimum in the 50–70% relative air humidity range in agreement with observations. The deposition velocity minimum could be explained with a reduced constant for the reaction of O₃ with water or OH⁻ ions on the surface.     

Modeling atmospheric mercury deposition in the vicinity of power plants

Two mathematical models of the atmospheric fate and transport of mercury (Hg), an Eulerian grid-based model and a Gaussian plume model, are used to calculate the atmospheric deposition of Hg in the vicinity (i.e., within 50 km) of five coal-fired power plants. The former is applied using two different horizontal resolutions: coarse (84 km) and fine (16.7 km). More than 96% of the power plant Hg emissions are calculated with the plume model to be transported beyond 50 km from the plants. The grid-based model predicts a lower fraction to be transported beyond 50 km: >91% with a coarse resolution and >95% with a fine resolution. The contribution of the power plant emissions to total Hg deposition within a radius of 50 km from the plants is calculated to be <8% with the plume model, <14% with the Eulerian model with a coarse resolution, and <10% with the Eulerian model with a fine resolution. The Eulerian grid-based model predicts greater local impacts than the plume model because of artificially enhanced vertical dispersion; the former predicts about twice as much Hg deposition as the latter when the area considered is commensurate with the resolution of the grid-based model. If one compares the local impacts for an area that is significantly less than the grid-based model resolution, then the grid-based model may predict lower local deposition than the plume model, because two compensating errors affect the results obtained with the grid-based model: initial dilution of the power plant emissions within one or more grid cells and enhanced vertical mixing to the ground.     

Current Status of Air Quality Models for Particulate Matter

ABSTRACT

The current status of the mathematical modeling of atmospheric particulate matter (PM) is reviewed in this paper. Simulating PM requires treating various processes, including the formation of condensable species, the gas/ particle partitioning of condensable compounds, and in some cases, the evolution of the particle size distribution. The algorithms available to simulate these processes are reviewed and discussed. Eleven 3-dimensional (3-D) Eulerian air quality models for PM are reviewed in terms of their formulation and past applications. Results of past performance evaluations of 3-D Eulerian PM models are presented. Currently, 24-hr average PM_2.5 concentrations appear to be predicted within 50% for urban-scale domains. However, there are compensating errors among individual particulate species. The lowest errors tend to be associated with SO₄ ^2-, while NO₃ ^-, black carbon (BC), and organic carbon (OC) typically show larger errors due to uncertainties in emissions inventories and the prediction of the secondary OC fraction. Further improvements and performance evaluations are recommended.     

Atmospheric fate of non-volatile and ionizable compounds

A modified version of the Multimedia Activity Model for Ionics MAMI, including two-layered atmosphere, air-water interface partitioning, intermittent rainfall and variable cloud coverage was developed to simulate the atmospheric fate of ten low volatility or ionizable organic chemicals. Probabilistic simulations describing the uncertainty of substance and environmental input properties were run to evaluate the impact of atmospheric parameters, ionization and air-water (or air-ice) interface enrichment.The rate of degradation and the concentration of OH radicals, the duration of dry and wet periods, and the parameters describing air-water partitioning (K_AW and temperature) and ionization (pK_a and pH) are the key parameters determining the potential for long range transport. Wet deposition is an important removal process, but its efficiency is limited, primarily by the duration of the dry period between precipitation events.Given the underlying model assumptions, the presence of clouds contributes to the higher persistence in the troposphere because of the capacity of cloud water to accumulate and transport non-volatile (e.g. 2,4-D) and surface-active chemicals (e.g. PFOA). This limits the efficiency of wet deposition from the troposphere enhancing long-range transport.     

Nighttime variations in HO2 radical mixing ratios at Rishiri Island observed with elevated monoterpene mixing ratios

HO₂ radical concentrations were measured by a laser-induced fluorescence instrument for three nighttime periods during the intensive field campaign at Rishiri Island, Japan, in June 2000. The HO₂ mixing ratio had temporal variations around its average of 4.2±1.2 (1σ) pptv and showed a positive correlation with the summed mixing ratio of four monoterpene species, α-pinene, β-pinene, camphene, and limonene, that sometimes reached 1 ppbv. Our model calculations suggested that ozonolysis reactions of monoterpenes were the main source of nighttime radicals and they explained 58% of measured HO₂ concentration levels. The model roughly reproduced the dependence of the HO₂ mixing ratio on the square root of the radical production rate due to the ozonolysis reactions of the monoterpenes. However, the absolute HO₂ mixing ratio was significantly underpredicted by the model. We discuss possible reasons in terms of misunderstood RO₂ chemistry, RO₂ interference with HO₂ observations, unknown radical production process associated by high NO₂ mixing ratio, and the contribution of unmeasured olefinic species to radical production via their reactions with ozone.     

Estimating contribution of precipitation scavenging of atmospheric particulate mercury to mercury wet deposition in Japan

Mercury wet deposition is dependent on both the scavenging of divalent reactive gaseous mercury (RGM) and atmospheric particulate mercury (Hg(p)) by precipitation. Estimating the contribution of precipitation scavenging of RGM and Hg(p) is important for better understanding the causes of the regional and seasonal variations in mercury wet deposition. In this study, the contribution of Hg(p) scavenging was estimated on the basis of the scavenging ratios of other trace elements (i.e., Cd, Cu, Mn, Ni, Pb and V) existing entirely in particulate form. Their wet deposition fluxes and concentrations in air, which were measured concurrently from April 2004 to March 2005 at 10 sites in Japan, were used in this estimation. The monthly wet deposition flux of mercury at each site correlated with the amount of monthly precipitation, whereas the Hg(p) concentrations in air tended to decrease during summer. There was a significant correlation (P<0.001) among the calculated monthly average scavenging ratios of trace elements, and the values in each month at each site were similar. Therefore, it is assumed the monthly scavenging ratio of Hg(p) is equivalent to the mean value of other trace elements. Using this scavenging ratio (W), the wet deposition flux (F) due to Hg(p) scavenging in each month was calculated by F=WKP, where K and P are the Hg(p) concentration and amount of precipitation, respectively. Relatively large fluxes due to Hg(p) scavenging were observed at a highly industrial site and at sites on the Japan Sea coast, which are strongly affected by the local sources and the long-range transport from the Asian continent, respectively. However, on average, at the 10 sites, the contribution of Hg(p) scavenging to the annual mercury deposition flux was 26%, suggesting that mercury wet deposition in Japan is dominated by RGM scavenging. This RGM should originate mainly from the in situ oxidation of Hg⁰ in the atmosphere.     

Plutonium partitioning in three-phase systems with water,granite grains,and different colloids

Low-solubility contaminants with high affinity for colloid surfaces may form colloid-associated species. The mobile characteristics of this species are, however, ignored by the traditional sorption/distribution experiments in which colloidal species contributed to the immobile fraction of the contaminants retained on the solids as a result of centrifugation or ultrafiltration procedures. The mobility of the contaminants in subsurface environments might be underestimated accordingly. Our results show that colloidal species of ²³⁹Pu in three-phase systems remained the highest percentages in comparison to both the dissolved species and the immobile species retained on the granite grains (solid phase), although the relative fraction of these three species depended on the colloid types. The real solid/liquid distribution coefficients (K _s/d) experimentally determined were generally smaller than the traditional K _s/d (i.e., the K _s+c/d in this study) by ~1,000 mL/g for the three-phase systems with the mineral colloids (granite particle, soil colloid, or kaolinite colloid). For the humic acid system, the traditional K _s/d was 140 mL/g, whereas the real K _s/d was approximately zero. The deviations from the real solid/liquid K _s/d were caused by the artificially increased immobile fraction of Pu. One has to be cautious in using K _s/d-based transport models to predict the fate and transport of Pu in the environment.     

A comparison of octanol-water partitioning between organic chemicals and their metabolites in mammals

Bioaccumulation models take various elimination and uptake processes into account, estimating rates from chemical lipophilicity, expressed as the octanol-water partition ratio (K_ow). Here, we focussed on metabolism, which transforms parent compounds into usually more polar metabolites, thus enhancing elimination. The aim of this study was to quantify the change in lipophilicity of relevant organic pollutants undergoing various biotransformation reactions in mammals. We considered oxidation reactions catalyzed by three enzyme groups: cytochrome P450 (CYP), alcohol dehydrogenase (ADH), and aldehyde dehydrogenase (ALDH). Estimated log K_ow values of a selected dataset of parent compounds were compared with the log K_ow of their first metabolites. The log K_ow decreased by a factor that varies between 0 and −2, depending on the metabolic pathway. For reactions mediated by CYP, the decrease in K_ow was one order of magnitude for hydroxylated and epoxidated compounds and two orders of magnitude for dihydroxylated and sulphoxidated xenobiotics. On the other hand, no significant change in lipophilicity was observed for compounds N-hydroxylated by CYP and for alcohols and aldehydes metabolized by ADH and ALDH. These trends could be anticipated by the calculus method of log K_ow. Yet, they were validated using experimental log K_ow values, when available. These relationships estimate the extent to which the elimination of pollutants is increased by biotransformation. Thus, the quantification of the K_ow reduction can be considered as a first necessary step in an alternative approach to anticipate biotransformation rates, which are hard to estimate with existing methods.     

Peroxidase-mediated removal of endocrine disrupting compound mixtures from water

Several classes of oxidative enzymes have shown promise for efficient removal of endocrine disrupting compounds (EDCs) that are resistant to conventional wastewater treatments. Although the kinetics of reactions between individual EDCs and selected oxidative enzymes are well documented in the literature, there has been little investigation of reactions with EDC mixtures. This makes it impossible to predict how enzyme-mediated treatment systems will perform since wastewater effluents generally contain multiple EDCs. This paper reports pseudo-first order rate constants for a model oxidative enzyme, horseradish peroxidase (HRP), during single-substrate (k₁) and mixed-substrate (k_1-MIX) reactions. Measured values are compared with literature values of three Michaelis-Menten parameters: half-saturation constant (K_M), enzyme turnover number (k_CAT), and the ratio k_CAT/K_M. Published reports had suggested that each of these could be correlated with HRP reactivity towards EDCs in mixtures, and empirical results from this study show that K_M can be used to predict the sequence of EDC removal reactions within a particular mixture. We also observed that k_1-MIX values were generally greater than k₁ values and that compounds exhibiting greatest estrogenic toxicities reacted most rapidly in a given mixture. Finally, because K_M may be tedious to measure for every EDC of interest, we have constructed a quantitative structure-activity relationship (QSAR) model to predict these values. This model predicts K_M quite accurately (R² = 89%) based on two molecular characteristics: molecular volume and hydration energy. Its accuracy makes this QSAR a useful tool for predicting which EDCs will be removed most efficiently during enzyme treatment of EDC mixtures.     

Adsorptive removal of an acid dye by lignocellulosic waste biomass activated carbon: equilibrium and kinetic studies

Chemically prepared activated carbon material derived from palm flower was used as adsorbent for removal of Amido Black dye in aqueous solution. Batch adsorption studies were performed for the removal of Amido Black 10B (AB10B), a di-azo acid dye from aqueous solutions by varying the parameters like initial solution pH, adsorbent dosage, initial dye concentration and temperature with three different particle sizes such as 100 μm, 600 μm and 1000 μm. The zero point charge was pH 2.5 and the maximum adsorption occurred at the pH 2.3. Experimental data were analyzed by model equations such as Langmuir, Freundlich and Temkin isotherms and it was found that the Freundlich isotherm model best fitted the adsorption data and the Freundlich constants varied from (K_F) 1.214, 1.077 and 0.884 for the three mesh sizes. Thermodynamic parameters such as ΔG, ΔH and ΔS were also calculated for the adsorption processes and found that the adsorption process is feasible and it was the endothermic reaction. Adsorption kinetics was determined using pseudo first-order, pseudo second-order rate equations and also Elovich model and intraparticle diffusion models. The results clearly showed that the adsorption of AB10B onto lignocellulosic waste biomass from palm flower (LCBPF) followed pseudo second-order model, and the pseudo second-order rate constants varied from 0.059 to 0.006 (g mg⁻¹ min) by varying initial adsorbate concentration from 25 mg L⁻¹ to 100 mg L⁻¹. Analysis of the adsorption data confirmed that the adsorption process not only followed intraparticle diffusion but also by the film diffusion mechanism.     

Evaluation of substrate removal kinetics for UASB reactors treating chlorinated ethanes

Purpose

Lack of focus on the treatment of wastewaters bearing potentially hazardous pollutants like 1,1,2 trichloroethane and 1,1,2,2 tetrachloroethane in anaerobic reactors has provided an impetus to undertake this study. The objective of this exercise was to quantify the behavior of upflow anaerobic sludge blanket reactors and predict their performance based on the overall organic substrate removal.

Methods

The reactors (wastewater-bearing TCA (R2), and wastewater-bearing TeCA (R3)) were operated at different hydraulic retention times (HRTs), i.e., 36, 30, 24, 18, and 12?h corresponding to food-to-mass ratios varying in the range of 0.2?C0.7?mg chemical oxygen demand (COD) mg^?1 volatile suspended solids day^?1. The process kinetics of substrate utilization was evaluated on the basis of experimental results, by applying three mathematical models namely first order, Grau second order, and Michaelis-Menten type kinetics.

Results

The results showed that the lowering of HRT below 24?h resulted in reduced COD removal efficiencies and higher effluent pollutant concentrations in the reactors. The Grau second-order model was successfully applied to obtain the substrate utilization kinetics with high value of R ² (>0.95). The Grau second-order substrate removal constant (K ₂) was calculated as 1.12 and 7.53?day^?1 for reactors R2 and R3, respectively.

Conclusion

This study demonstrated the suitability of Grau second-order kinetic model over other models, for predicting the performance of reactors R2 and R3, in treating wastewaters containing chlorinated ethanes under different organic and hydraulic loading conditions.     

Simulation of aerosols and gas-phase species over Europe with the Polyphemus system. Part II: Model sensitivity analysis for 2001

This paper presents a multi-pollutant sensitivity study of an air quality model over Europe with a focus on aerosols. Following the evaluation presented in the companion paper, the aim here is to study the sensitivity of the model to input data, mathematical parameterizations and numerical approximations. To that end, 30 configurations are derived from a reference configuration of the model by changing one input data set, one parameterization or one numerical approximation at a time. Each of these configurations is compared to the same reference simulation over two time periods of the year 2001, one in summer and one in winter. The sensitivity of the model to the different configurations is evaluated through a statistical comparison between the simulation results and through comparisons to available measurements. The species studied are ozone (O₃), nitrogen dioxide (NO₂), sulfur dioxide (SO₂), ammonia (NH₃), coarse and fine aerosol particles (PM_c and PM_2.5), sulfate, nitrate, ammonium, chloride and sodium.For all species, the modeled concentrations are very sensitive to the parameterization used for vertical turbulent diffusion and to the number of vertical levels. For the other configurations considered in this work, the sensitivity of the modeled concentration to configuration choice varies with the species and the period of the year. O₃ is impacted by options related to boundary conditions. PM_c is sensitive to sea-salt related options, to options influencing deposition and to options related to mass transfer between gas and particulate phases. PM_2.5 is sensitive to a larger number of options than PM_c: sea-salt, boundary conditions, heterogeneous reactions, aqueous chemistry and gas/particle mass transfer. NO₂ is strongly influenced by heterogeneous reactions. Nitrate shows the highest variability of all species studied. As with NO₂, nitrate is strongly sensitive to heterogeneous reactions but also to mass transfer, thermodynamic related options, aqueous chemistry and computation of the wet particle diameter. While SO₂ is mostly sensitive to aqueous chemistry, sulfate is also sensitive to boundary conditions and, to a lesser extent, to heterogeneous reactions. As with nitrate, ammonium is largely impacted by the different configuration choices, although the sensitivity is slightly lower than for nitrate. NH₃ is sensitive to aqueous chemistry, mass transfer and heterogeneous reactions. Chloride and sodium are impacted by sea-salt related options, by options influencing deposition and by options concerning the aqueous-phase module.