Exchange of sulfur pollution between Bulgaria and Greece

The EMAP model (Eulerian Model for Air Pollution) is applied for calculating the sulfur concentration and deposition fields for 1995 as based upon Bulgarian and Greek sources. The country-to-country budgets show that about 4% of the sulfur oxides emitted by Bulgaria are deposited over Greek territory, estimated at 28 kt S. Only 2% of sulfur compounds emitted by Greece are deposited over Bulgaria, estimated at 6.2 kt S for 1995. This data is in agreement with EMEP/MSC-W estimates and provides more details concerning time and space. The results can be used in decision-making, negotiating and the development of contamination strategies.     

Evaluating non-equilibrium solute transport in small soil columns

Displacement studies on leaching of bromide and two pesticides (atrazine and isoproturon) were conducted under unsaturated steady state flow conditions in 24 small undisturbed soil columns (5.7 cm in diameter and 10 cm long) each collected from two sites differing in soil structure and organic carbon content in North Germany. There were large and irregular variabilities in the characteristics of both soils, as well as in the shapes of breakthrough curves (BTCs) of different columns, including some with early breakthrough and increased tailing, qualitatively indicating the presence of preferential flow. It was estimated that one preferential flow column (PFC) at site A, and four at site B, contributed, respectively to 11% and 58% of the accumulated leached fraction and to more than 80% of the maximum observed standard deviation (SD) in the field-scale concentration and mass flux of pesticides at two sites. The bromide BTCs of two sites were analyzed with the equilibrium convection-dispersion equation (CDE) and a non-equilibrium two-region/mobile-immobile model. Transport parameters of these models for individual BTCs were determined using a curve fitting program, CXTFIT, and by the time moment method. For the CDE based equilibrium model, the mean values of retardation factor, R, considered separately for all columns, PFCs or non-preferential flow columns (NPFCs) were comparable for the two methods; significant differences were observed in the values of dispersion coefficients of two sites using the two estimation methods. It was inferred from the estimated parameters of non-equilibrium model that 5-12% of water at site A, and 12% at site B, was immobile during displacement in NPFCs. The corresponding values for PFCs of two sites were much larger, ranging from 25% to 51% by CXTFIT and from 24% to 72% by the moment method, suggesting the role of certain mechanisms other than immobile water in higher degrees of non-equilibrium in these columns. Peclet numbers in PFCs of both sites were consistently smaller than five, indicating the inadequacy of the non-equilibrium model to incorporate the effect of all forms of non-equilibrium in PFCs. Overall, the BTCs of individual NPFCs, PFCs and of field average concentration at the two sites were better reproduced with parameters obtained from CXTFIT than by the moment method. The moment method failed to capture the peak concentrations in PFCs, but tended to describe the desorption and tail branches of BTCs better than the curve fitting approach.     

Biodegradation during contaminant transport in porous media: 4. Impact of microbial lag and bacterial cell growth

Miscible-displacement experiments were conducted to examine the impact of microbial lag and bacterial cell growth on the transport of salicylate, a model hydrocarbon compound. The impacts of these processes were examined separately, as well as jointly, to determine their relative effects on biodegradation dynamics. For each experiment, a column was packed with porous medium that was first inoculated with bacteria that contained the NAH plasmid encoding genes for the degradation of naphthalene and salicylate, and then subjected to a step input of salicylate solution. The transport behavior of salicylate was non-steady for all cases examined, and was clearly influenced by a delay (lag) in the onset of biodegradation. This microbial lag, which was consistent with the results of batch experiments, is attributed to the induction and synthesis of the enzymes required for biodegradation of salicylate. The effect of microbial lag on salicylate transport was eliminated by exposing the column to two successive pulses of salicylate, thereby allowing the cells to acclimate to the carbon source during the first pulse. Elimination of microbial lag effects allowed the impact of bacterial growth on salicylate transport to be quantified, which was accomplished by determining a cell mass balance. Conversely, the impact of microbial lag was further investigated by performing a similar double-pulse experiment under no-growth conditions. Significant cell elution was observed and quantified for all conditions/systems. The results of these experiments allowed us to differentiate the effects associated with microbial lag and growth, two coupled processes whose impacts on the biodegradation and transport of contaminants can be difficult to distinguish.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

Ozone patterns for three metropolitan statistical areas in North Carolina, USA

As part of an effort by the state of North Carolina to develop a State Implementation Plan (SIP) for 1-h peak ozone control, a network of ozone stations was established to monitor surface ozone concentrations across the state. Between 19 and 23 ozone stations made continuous surface measurements between 1993 and 1995 surrounding three major metropolitan statistical areas (MSAs): Raleigh/Durham (RDU), Charlotte/Mecklenburg (CLT), and Greensboro/High Point/Winston-Salem (GSO). Statistical averages of the meteorological and ozone data were performed at each Metropolitan Statistical Area (MSA) to study trends and/or relationships on high ozone days (days in which one of the MSA sites measured an hourly ozone concentration90.0 ppbv). County emission maps of precursor gases, wind roses, total area averages of ozone, total downwind averages of ozone deviations, upwind averages of ozone, and a modified delta ozone analysis were all obtained and analyzed. The results of this study show a reduction in the delta ozone relative to an earlier study at RDU, but no average significant change at CLT (no comparison can be made for GSO). The statistical data analyses in this study are used to quantify the importance of local contributions and regional transport, to ozone air pollution in the MSAs.     

应用TOXI模型进行重金属在河流中迁移转化规律的研究

ＴＯＸＩ模型是美国环保局的ＷＡＳＰ（水质分析模拟程序）中的一部分。由于该模型考虑了污染物质在河流中溶解态，颗粒态和底泥态的分布关系，对于模拟重金属在河流中的浓度变化效果较好。为保证计算结果的可靠性，首先应对重金属浓度的衰减变化特性进行研究，本文重点讨论的是重金属在河流中迁移转化的过程，并如何将这一过程用数学公式来进行模拟。     

On the hydro-dispersive equivalence between multi-layered mineral barriers

In the context of municipal solid waste and hazardous waste disposal, the notion of "equivalence" between different barrier designs appears in regulatory documents from several industrialized countries. While in the past, equivalence has been thought of mainly in terms of contaminant travel times, in recent years it has been defined more in terms of the magnitude of a disposal site's potential impact on groundwater resources. This paper presents some original analytical solutions to the problem of contaminant migration through a multi-layered mineral barrier. The solutions account for the two major mechanisms of subsurface contaminant migration, namely, advection and diffusion-dispersion. An example application using the proposed solutions and a numerical model illustrates how one multi-layered mineral barrier can be considered superior to another from a strictly hydro-dispersive viewpoint. The influence of partial saturation of the mineral barrier is investigated using a numerical solution to the Richards equation for unsaturated flow. It is emphasized that conclusions relative to the superiority of one multi-layered barrier, with respect to another, should not only consider hydro-dispersive aspects, but also other processes such as the mechanical and chemical evolutions of the different barrier components. Although such phenomena are poorly addressed by existing models, failure to take them into account, at least in a qualitative fashion, may lead to unconservative conclusions with respect to barrier equivalence.     

Intercomparison of reactive transport models applied to UO2 oxidative dissolution and uranium migration

Oxidative dissolution of uranium dioxide (UO₂) and the subsequent migration of uranium in a subsurface environment and an underground waste disposal have been simulated with reactive transport models. In these systems, hydrogeological and chemical processes are closely entangled and their interdependency has been analyzed in detail, notably with respect to redox reactions, kinetics of mineralogical evolution and hydrodynamic migration of species of interest.Different codes, where among CASTEM, CHEMTRAP and HYTEC, have been used as an intercomparison and verification exercise. Although the agreement between codes is satisfactory, it is shown that the discretization method of the transport equation (i.e. finite elements (FE) versus mixed-hybrid FE and finite differences) and the sequential coupling scheme may lead to systematic discrepancies.     

Transport modeling of copper and cadmium with linear and nonlinear retardation factors

The predictive accuracy of using the one-dimensional advection–dispersion equation to evaluate the fate and transport of solute in a soil column is usually dependent on the proper determination of chemical retardation factors. Typically, the distribution coefficient (K_d) obtained by fitting the linear sorption isotherm has been extensively used to consider general geochemical reactions on solute transport in a low-concentration range. However, the linear distribution coefficient cannot be adequately utilized to describe the solute fate at a higher concentration level. This study employed the nonlinear equilibrium-controlled sorption parameters to determine the retardation factor used in column leaching experiments. Copper and cadmium transportation in a lateritic silty-clay soil column was examined. Through the explicit finite-difference calculations with a third-order total-variation-diminishing (TVD) numerical solution scheme, all results of the theoretical copper and cadmium breakthrough curves (BTCs) simulated by using the Freundlich nonlinear retardation factors revealed good agreement with the experimental observations.     

污染场地地下水中汞污染反应运移模拟

汞（Hg）是人们持续关注的全球环境污染物之一,其对地下水的污染严重威胁着与地下水相关的生态环境系统.汞在地下水系统中的物理与地球化学反应过程的准确刻画是研究汞迁移转化规律的重点和难点.基于某工业场地汞污染数据,首先采用PHREEQC研究地下水中无机二价汞的存在形态,然后利用PHT3D程序建立汞污染物反应性溶质运移二维剖面模型.该模型考虑了汞污染物在地下水系统中的对流、弥散过程及地球化学反应过程（包括水相络合作用、表面络合吸附作用及受动力学控制的氧化还原作用）.结果表明,无机二价汞的存在形态以HgCl₂和Hg（OH）Cl占主导地位,氧化还原作用是影响地下水中汞污染反应性运移的主要控制因素;另一方面,水合氧化铁HFO对汞迁移的阻滞影响较小,而溶解性有机质对汞较强的络合作用不能忽视.本文研究成果可为预测与评估特定污染场地地下水汞污染的变化趋势及制定相应的修复策略提供科学依据.