Development of a nonlocal convective mixing scheme with varying upward mixing rates for use in air quality and chemical transport models

BACKGROUND, AIM, AND SCOPE: Asymmetrical convective non-local scheme (CON) with varying upward mixing rates is developed for simulation of vertical turbulent mixing in the convective boundary layer in air quality and chemical transport models. MATERIALS AND METHODS: The upward mixing rate form the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. RESULTS: This scheme provides a less rapid mass transport out of surface layer into other layers than other asymmetrical convective mixing schemes. DISCUSSION: In this paper, we studied the performance of a nonlocal convective mixing scheme with varying upward mixing in the atmospheric boundary layer and its impact on the concentration of pollutants calculated with chemical and air-quality models. This scheme was additionally compared versus a local eddy-diffusivity scheme (KSC). Simulated concentrations of NO(2) and the nitrate wet deposition by the CON scheme are closer to the observations when compared to those obtained from using the KSC scheme. CONCLUSIONS: Concentrations calculated with the CON scheme are in general higher and closer to the observations than those obtained by the KSC scheme (of the order of 15-20%). Nitrate wet deposition calculated with the CON scheme are in general higher and closer to the observations than those obtained by the KSC scheme. RECOMMENDATIONS AND PERSPECTIVES: To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO(2)) and nitrate wet deposition was compared for the year 2002. The comparison was made for the whole domain used in simulations performed by the chemical European Monitoring and Evaluation Programme Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.     

Development of a transport phenomena model for accidental releases of heavy gases in an industrial environment

The present paper discusses the development of a transport phenomena model for the dispersion of heavy gases. It describes the three-dimensional transient dispersion processes of an accidentally released contaminant. The model accounts for heavy gas effects and allows one to explicitly resolve the effect of twodimensional buildings in the ambient flow field. The results of the simulation are compared with those of other calculation procedures for conditions where the different approaches are applicable. In addition, comparisons are made with results from field experiments on heavy gas dispersion.     

Development of an in-line system for the analysis of 4,4'-DDT in water

An in-line system for trace persistent organic pollutants (POPs) in water was developed by using a laboratory-made hollow fiber membrane (HFM) unit connected with a high-resolution gas chromatograph-mass spectrometer (HRGC-MS). The semivolatile organic compound, 4,4'-Dichlorodiphenyl trichloroethane (4,4'-DDT), was chosen as a representative of a persistent organic compound. The synthetic water contaminated with 4,4'-DDT was passed through the HFM unit, the extraction occurred by the analyte pervaporated and permeated, then stripped into HRGC-MS. Several factors were investigated for the high extraction efficiency. The best performance was obtained at sample and stripping gas flow rates of 6 and 9 mLmin-1, respectively, and desorption temperature of 60 degrees C. At this temperature, the diffusion rate was enhanced by 15 times over 25 degrees C. A wide linear dynamic range was obtained, i.e., 0.10-1.0 mgL-1, with a limit of detection (LOD) of 90 microgL-1. The extraction efficiency of 4,4'-DDT in real water samples was in the range of 83-94%. Real water samples were analyzed and 0.6 microgL-1 of 4,4'-DDT was found in unregistered bottled water and 7.0 microgL-1 in tap water.     

Evaluation of an aerosol optical scheme in the chemistry-transport model CHIMERE

This paper presents an aerosol optical scheme developed in the chemistry-transport model CHIMERE dedicated to calculate optical properties of particles. Such developments are very helpful as they complement the usual validation with PM (Particulate Matter) ground-based measurements by using surface (AERONET/PHOTONS network) and satellite (MODIS) remote sensing observations. To reach this goal, Aerosol Optical Thickness (AOT), column-averaged Single Scattering Albedo (SSA) and asymmetry parameter (g) are calculated at 440 nm, 675 nm, 870 nm and 1020 nm (AERONET wavelengths) under three hypotheses on the particle mixing state (external, internally homogeneous and core-shell). Furthermore and in addition to optical calculations, an original development has been made to estimate column volume size distributions in CHIMERE, directly comparable with AERONET retrievals. Comparisons between simulations and observations are made over Western Europe for the year 2003 but also for one specific case focused on ammonium nitrate aerosols. Observed AOT display a seasonal cycle (with highest values during summer) rather well reproduced by the model but biases with observational data have been found depending on seasons. In fall, winter and early spring, modeled AOT values agree well with AERONET retrievals with small negative biases. Focus on a pollution episode of ammonium nitrate origin during March 2003 reveals that CHIMERE is able to well reproduce the fine mode volume size distribution retrieved by AERONET, leading to good agreements between modeled and observed AOT. In late spring and summer, AERONET AOT values are underpredicted by the model, which could be due to uncertainties in modeling secondary species.     

Development and application of a three-dimensional aerosol chemical transport model,PMCAMx

A three-dimensional chemical transport model (PMCAMx) is used to simulate PM mass and composition in the eastern United States for a July 2001 pollution episode. The performance of the model in this region is evaluated, taking advantage of the highly time and size-resolved PM and gas-phase data collected during the Pittsburgh Air Quality Study (PAQS). PMCAMx uses the framework of CAMx and detailed aerosol modules to simulate inorganic aerosol growth, aqueous-phase chemistry, secondary organic aerosol formation, nucleation, and coagulation. The model predictions are compared to hourly measurements of PM_2.5 mass and composition at Pittsburgh, as well as to measurements from the AIRS and IMPROVE networks. The performance of the model for the major PM_2.5 components (sulfate, ammonium, and organic carbon) is encouraging (fractional errors are in general smaller than 50%). Additional improvements are possible if the rainfall measurements are used instead of the meteorological model predictions. The modest errors in ammonium predictions and the lack of bias for the total (gas and particulate) ammonium suggest that the improved ammonia inventory used is reasonable. The significant errors in aerosol nitrate predictions are mainly due to difficulties in simulating the nighttime formation of nitric acid. The concentrations of elemental carbon (EC) in the urban areas are significantly overpredicted. This is a problem related to both the emission inventory but also the different EC measurement methods that have been used in the two measurement networks (AIRS and IMPROVE) and the actual development of the inventory. While the ability of the model to reproduce OC levels is encouraging, additional work is necessary to confirm that that this is due to the right reasons and not offsetting errors in the primary emissions and the secondary formation. The model performance against the semi-continuous measurements in Pittsburgh appears to be quite similar to its performance against daily average measurements in a wide range of stations across the Eastern US. This suggests that the skill of the model to reproduce the diurnal variability of PM_2.5 and its major components is as good as its ability to reproduce the daily average values and also the significant value of high temporal resolution measurements for model evaluation.     

Numerical analysis of water and solute transport in variably-saturated fractured clayey till

This study numerically investigates the influence of initial water content and rain intensities on the preferential migration of two fluorescent tracers, Acid Yellow 7 (AY7) and Sulforhodamine B (SB), through variably-saturated fractured clayey till. The simulations are based on the numerical model HydroGeoSphere, which solves 3D variably-saturated flow and solute transport in discretely-fractured porous media. Using detailed knowledge of the matrix, fracture, and biopore properties, the numerical model is calibrated and validated against experimental high-resolution tracer images/data collected under dry and wet soil conditions and for three different rain events. The model could reproduce reasonably well the observed preferential migration of AY7 and SB through the fractured till, although it did not capture the exact depth of migration and the negligible impact of the dead-end biopores in a near-saturated matrix. A sensitivity analysis suggests fast flow mechanisms and dynamic surface coating in the biopores, and the presence of a plough pan in the till.     

An analysis of the sulfur budget and interstate sulfur transport for Colorado

The sulfur budget for Colorado is analyzed in order to understand the importance of interstate transport and natural sources to wet sulfur deposition in the state. A climatological box model is used to separate the various contributions to sulfur deposition. We demonstrate that sulfur dioxide molecules emitted in Colorado have a low probability of wet deposition in Colorado. Therefore, as much as 80 per cent or more of the wet sulfur deposition in the state must be ascribed to external anthropogenic sulfur sources, such as non-ferrous metal smelters and electric powerplants unless some large unspecified natural sulfur source exists.     

Development of an improved optical transmission technique for black carbon (BC) analysis

A new optical transmission technique for black carbon (BC) analysis was developed to minimize interferences due to scattering effects in filter samples. A standard thermal analysis method (VDI, 1999) is used to link light attenuation by the filter samples to elemental carbon (EC) concentration. Scattering effects are minimized by immersion of the filters in oil of a similar refractive index, as is often done for microscopy purposes. Light attenuation was measured using both a white light source and a red LED of 650 nm. The usual increase in overestimation of BC concentrations with decreasing BC amount in filter samples was found considerably reduced. Some effects of BC properties (e.g. fractal dimension, microstructure and size distribution) on the specific attenuation coefficient B_ATN, however, are still present for the treated samples. B_ATN was found close to 1 m² g⁻¹ for dry-dispersed industrial BC and 7 m² g⁻¹ for nebulized BC. Good agreement was found between the oil immersion, integrating sphere and a polar photometer technique and Mie calculations. The average specific attenuation coefficient of ambient samples in oil varied between 7 and 11 m² g⁻¹ for white light and 6 and 9 m² g⁻¹ for red light (LED). B_ATN was found to have much less site variation for the treated than for the untreated samples. The oil immersion technique improved also the correlation with thermally analyzed EC. This new immersion technique therefore presents a considerable improvement over conventional optical transmission techniques and may therefore serve as a simple, fast and cost-effective alternative to thermal methods.     

A size-segregated particle dry deposition scheme for an atmospheric aerosol module

A parameterization of particle dry deposition has been developed for the Canadian Aerosol Module (CAM).This parameterization calculates particle dry deposition velocities as a function of particle size and density as well as relevant meteorological variables. It includes deposition processes, such as, turbulent transfer, Brownian diffusion, impaction, interception, gravitational settling and particle rebound. Particle growth under humid conditions is also considered. Sensitivity tests show that the parameterization provides deposition velocities comparable with recent field observations, especially for sub-micron particles. The present parameterization has also been evaluated using two empirical bulk resistance models, which were originally developed from field observations. The present parameterization has been implemented in CAM, with meteorological input provided by the Canadian Regional Climate Model (RCM) to the eastern North America. A comparison of the modelled dry deposition velocities to a variety of recent measurements that have been reported in the literature demonstrated that the current parameterization produces reasonable results. The main improvement of the current parameterization compared to earlier size-dependent particle dry deposition models is that the current one produces more realistic deposition velocities for sub-micron particles and agrees better with recently published field measurements.     

Study of contaminant transport at an open-tipping waste disposal site

Field and laboratory studies were conducted to estimate concentration of potential contaminants from landfill in the underlying groundwater, leachate, and surface water. Samples collected in the vicinity of the landfill were analyzed for physiochemical parameters, organic contaminants, and toxic heavy metals. Water quality results obtained were compared from published data and reports. The results indicate serious groundwater and surface water contamination in and around the waste disposal site. Analysis of the organic samples revealed that the site contains polychlorinated biphenyls and other organo-chlorine chemicals, principally chloro-benzenes. Although the amount of PCB concentration discovered was not extreme, their presence indicates a potentially serious environmental threat. Elevated concentrations of lead, copper, nickel, manganese, cadmium, and cobalt at the downgradient indicate that the contamination plume migrated further from the site, and the distribution of metals and metals containing wastes in the site is nonhomogeneous. These results clearly indicate that materials are poorly contained and are at risk of entering the environment. Therefore, full characterization of the dump contents and the integrity of the site are necessary to evaluate the scope of the problem and to identify suitable remediation options.     

土壤中有机污染物运移模型的对比分析

分别建立了考虑非平衡吸附和生物降解的单区域溶质运移动力学模型(简称单区域模型)和两区域溶质运移动力学模型(简称两区域模型),基于数值模拟方法,对两种模型进行了定量的对比分析.分析结果表明,参与流动的水分的含量、不可动水区的持留释放能力是影响污染物在土壤-水环境中分配规律以及两种模型模拟结果之间差异的重要因素;较大的可动水区和不可动水区之间的质量转移系数,可以通过增加持留释放能力而增强不可动水区对溶质运移的影响能力,因此在模拟过程中不应单纯依据不可动水的含量大小来确定采用何种模型进行模拟,而应同时考虑可动水区和不可动水区之间的质量转移系数对模拟效果的影响.     

Analytical solutions for reactive transport under an infiltration-redistribution cycle

Transport of reactive solute in unsaturated soils under an infiltration-redistribution cycle is investigated. The study is based on the model of vertical flow and transport in the unsaturated zone proposed by Indelman et al. [J. Contam. Hydrol. 32 (1998) 77], and generalizes it by accounting for linear nonequilibrium kinetics. An exact analytical solution is derived for an irreversible desorption reaction. The transport of solute obeying linear kinetics is modeled by assuming equilibrium during the redistribution stage. The model which accounts for nonequilibrium during the infiltration and assumes equilibrium at the redistribution stage is termed partial equilibrium infiltration-redistribution model (PEIRM). It allows to derive approximate closed form solutions for transport in one-dimensional homogeneous soils. These solutions are further applied to computing the field-scale concentration by adopting the Dagan and Bresler [Soil Sci. Soc. Am. J. 43 (1979) 461] column model. The effect of soil heterogeneity on the solute spread is investigated by modeling the hydraulic saturated conductivity as a random function of horizontal coordinates. The quality of the PEIRM is illustrated by calculating the critical values of the Damk?hler number which provide the achievable accuracy in estimating the solute mass in the mobile phase. The distinguishing feature of transport during the infiltration-redistribution cycle as compared to that of infiltration only is the finite depth of solute penetration. For irreversible desorption, the maximum solute penetration W/theta(r) is determined by the amount of applied water W and the residual water content theta(r). For sorption-desorption kinetics, the maximum depth of penetration z(r)(e, infinity ) also depends on the ratio between the rate of application and the column-saturated conductivity. It is shown that z(r)(e, infinity ) is bounded between the depths W/(theta(r)+K(d)) and W/theta(r) corresponding to the maximum solute penetration for equilibrium transport and for irreversible desorption, respectively. This feature of solute penetration explains the unusual phenomena of plume contraction after an initial period of spreading [Lessoff, S.C., Indelman, P., Dagan, G., 2002. Solute transport in infiltration-redistribution cycles in heterogeneous soils. In Raats, P.A.C., Smiles, D.,Warrick, A.W. (Eds), Environmental Mechanics: Water, Mass and Energy Transport in the Biosphere. American Geophysical Union, pp. 133-144]. Unlike transport under equilibrium conditions, when the solute is completely concentrated at the front, the solute under nonequilibrium conditions is spread out behind the front. Heterogeneity leads to additional spreading of the plume.     

Simulation of dispersion of a power plant plume using an adaptive grid algorithm

A new dynamic adaptive grid algorithm has been developed for use in air quality modeling. This algorithm uses a higher order numerical scheme—the piecewise parabolic method (PPM)—for computing advective solution fields; a weight function capable of promoting grid node clustering by moving grid nodes; and a conservative interpolation equation using PPM for redistributing the solution field after movement of grid nodes. Applications of the algorithm to a model problem, in which emissions from a point source disperse through the atmosphere in time, reflect that the algorithm is able to capture not only the regional ozone plume distribution, but also the small-scale plume structure near the source. In contrast, the small-scale plume structure was not captured in the corresponding static grid solution. Performance achieved in model problem simulations indicates that the algorithm has the potential to provide accurate air quality modeling solutions at costs that may be significantly less than those incurred in obtaining equivalent static grid solutions.     

Multi-process herbicide transport in structured soil columns: experiments and model analysis

Model predictions of pesticide transport in structured soils are complicated by multiple processes acting concurrently. In this study, the hydraulic, physical, and chemical nonequilibrium (HNE, PNE, and CNE, respectively) processes governing herbicide transport under variably saturated flow conditions were studied. Bromide (Br-), isoproturon (IPU, 3-(4-isoprpylphenyl)-1,1-dimethylurea) and terbuthylazine (TER, N2-tert-butyl-6-chloro-N4-ethyl-1,3,5-triazine-2,4-diamine) were applied to two soil columns. An aggregated Ap soil column and a macroporous, aggregated Ah soil column were irrigated at a rate of 1 cm h(-1) for 3 h. Two more irrigations at the same rate and duration followed in weekly intervals. Nonlinear (Freundlich) equilibrium and two-site kinetic sorption parameters were determined for IPU and TER using batch experiments. The observed water flow and Br- transport were inversely simulated using mobile-immobile (MIM), dual-permeability (DPM), and combined triple-porosity (DP-MIM) numerical models implemented in HYDRUS-1D, with improving correspondence between empirical data and model results. Using the estimated HNE and PNE parameters together with batch-test derived equilibrium sorption parameters, the preferential breakthrough of the weakly adsorbed IPU in the Ah soil could be reasonably well predicted with the DPM approach, whereas leaching of the strongly adsorbed TER was predicted less well. The transport of IPU and TER through the aggregated Ap soil could be described consistently only when HNE, PNE, and CNE were simultaneously accounted for using the DPM. Inverse parameter estimation suggested that two-site kinetic sorption in inter-aggregate flow paths was reduced as compared to within aggregates, and that large values for the first-order degradation rate were an artifact caused by irreversible sorption. Overall, our results should be helpful to enhance the understanding and modeling of multi-process pesticide transport through structured soils during variably saturated water flow.     

Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test

Realistic models of contaminant transport in groundwater demand detailed characterization of the spatial distribution of subsurface hydraulic properties, while at the same time programmatic constraints may limit collection of pertinent hydraulic data. Fortunately, alternate forms of data can be used to improve characterization of spatial variability. We utilize a methodology that augments sparse hydraulic information (hard data) with more widely available hydrogeologic information to generate equiprobable maps of hydrogeologic properties that incorporate patterns of connected permeable zones. Geophysical and lithologic logs are used to identify hydrogeologic categories and to condition stochastic simulations using Sequential Indicator Simulation (SIS). The resulting maps are populated with hydraulic conductivity values using field data and Sequential Gaussian Simulation (SGS). Maps of subsurface hydrogeologic heterogeneity are generated for the purpose of examining groundwater flow and transport processes at the Faultless underground nuclear test, Central Nevada Test Area (CNTA), through large-scale, three-dimensional numerical modeling. The maps provide the basis for simulation of groundwater flow, while transport of radionuclides from the nuclear cavity is modeled using particle tracking methods. Sensitivity analyses focus on model parameters that are most likely to reduce the long travel times observed in the base case. The methods employed in this study have improved our understanding of the spatial distribution of preferential flowpaths at this site and provided the critical foundation on which to build models of groundwater flow and transport. The results emphasize that the impacts of uncertainty in hydraulic and chemical parameters are dependent on the radioactive decay of specific species, with rapid decay magnifying the effects of parameters that change travel time.     

Electrical resistivity tomography as monitoring tool for unsaturated zone transport: an example of preferential transport of deicing chemicals

Non-invasive spatially resolved monitoring techniques may hold the key to observe heterogeneous flow and transport behavior of contaminants in soils. In this study, time-lapse electrical resistivity tomography (ERT) was employed during an infiltration experiment with deicing chemical in a small field lysimeter. Deicing chemicals like potassium formate, which frequently impact soils on airport sites, were infiltrated during snow melt. Chemical composition of seepage water and the electrical response was recorded over the spring period 2010. Time-lapse electrical resistivity tomographs are able to show the infiltration of the melt water loaded with ionic constituents of deicing chemicals and their degradation product hydrogen carbonate. The tomographs indicate early breakthrough behavior in parts of the profile. Groundtruthing with pore fluid conductivity and water content variations shows disagreement between expected and observed bulk conductivity. This was attributed to the different sampling volume of traditional methods and ERT due to a considerable fraction of immobile water in the soil. The results show that ERT can be used as a soil monitoring tool on airport sites if assisted by common soil monitoring techniques.     

Possible evidence for transport of an iron cyanide complex by plants

Barley (Hordeum vulgare L.), oat (Avena sativa L.), and wild cane (Sorghum bicolor L.), were exposed to 15N-labeled ferrocyanide to determine whether these plant species can transport this iron cyanide complex. Plants were treated with ferrocyanide in a nutrient solution that simulated iron cyanide contaminated groundwater and soil solutions. This nutrient solution has been shown to maintain ferrocyanide speciation with minimal dissociation to free cyanide. Following treatment, all three plants showed dramatic enrichments in roots (delta 15N per thousand =1000-1500) and shoots (delta 15N per thousand =500). Barley and oat showed enrichment primarily in roots while wild cane showed a near equal enrichment in root and shoot tissues. Nitrogen-deficient barley plants treated with ferrocyanide showed a significantly greater 15N enrichment as compared to nitrogen-sufficient plants. While the results are suggestive of ferrocyanide transport by these plant species, additional study will be required to verify these results.     

Speciation and transport of arsenic in an acid sulfate soil-dominated catchment, eastern Australia

Factors controlling the transport of geogenically-derived arsenic from a coastal acid sulfate soil into downstream sediments are identified in this study with both solid-phase associations and aqueous speciation clearly critical to the mobility and toxicity of arsenic. The data from both sequential extractions and X-ray adsorption spectroscopy indicate that arsenic in the unoxidised Holocene acid sulfate soils is essentially non-labile in the absence of prolonged oxidation, existing primarily as arsenopyrite or as an arsenopyrite-like species, likely arsenian pyrite. Anthropogenically-accelerated pedogenic processes, which have oxidised this material over time, have greatly enhanced the potential bioavailability of arsenic, with solid-phase arsenic almost solely present as As(V) associated with secondary Fe(III) minerals present. Analyses of downstream sediments reveal that a portion of the arsenic is retained as a mixed As(III)/As(V) solid-phase, though not at levels considered to be environmentally deleterious. Determination of arsenic speciation in pore waters using high performance liquid chromatography/Inductively Coupled Plasma-Mass Spectrometry shows a dominance of As(III) in upstream pore waters whilst an unidentified As species reaches comparative levels within the downstream, estuarine locations. Pore water As(V) was detected at trace concentrations only. The results demonstrate the importance of landscape processes to arsenic transport and availability within acid sulfate soil environments.