Chemical transport models

Background, aim, and scope  Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. Materials and methods  The upward mixing rate forming 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. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. Results  This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. Discussion  In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants calculated with chemical and air-quality models. In addition, this scheme was also compared with a commonly used, local, eddy-diffusivity scheme. Simulated concentrations of NO₂ by the COM scheme and new parameterization of the in-canopy resistance are closer to the observations when compared to those obtained from using the local eddy-diffusivity scheme. Conclusions  Concentrations calculated with the COM scheme and new parameterization of in-canopy resistance, are in general higher and closer to the observations than those obtained by the local, eddy-diffusivity scheme (on the order of 15–22%). Recommendations and perspectives  To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO₂) were compared for the years 1999 and 2002. The comparison was made for the entire domain used in simulations performed by the chemical European Monitoring and Evaluation Program Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.     

Atmospheric pollution transport: the parallelization of a transport & chemistry code

We discuss a model for computing the chemical reactions occurring among air pollutants and predict their dispersion in the atmosphere. We consider a parallel implementation on a CRAY T3D, using up to 128 processors. A regular subdomain decomposition is considered to optimize inter-processor communications, while load balancing – the main technical difficulty to obtain a good efficiency – is achieved by re-assigning temporarily pieces of work to other processors, according to a global strategy. The performance of our implementation is described by an analytical relation whose validity is checked on realistic data sets. A speedup of 100 is obtained with 128 processors.     

Southeastward transport of Asian dust: Source,transport and its contributions to Taiwan

Atmospheric Aluminum measured in northern Taiwan from 2003 to 2006 is used as a dust tracer, from which dust concentrations are derived, and major Asian dust events are determined. The source locations for the major dust events are traced back and identified, and the processes leading to the southeastward transport of Asian dust is investigated. The derived dust concentrations are compared to the local PM₁₀ (particle with size less than 10 μm) concentrations, and the impacts of Asian dust on the air quality of Taiwan are quantified.According to the backward trajectory and dust observation analyses, most of the southeastward transport of major Asian dust events originate from Mongolia and Inner Mongolia in northern China, and only one out of 16 events is generated from western China. Modeling studies and weather analyses of dust events suggest that the southeastward transport of Asian dust is usually generated behind a surface front and transported downwind behind the associated upper level trough. The associated upper level trough is usually deep, in which the northwesterly wind behind the trough favors the southeastward transport of dust to lower latitudes. Dust transported to Taipei generally occur during periods of large-scale subsidence.Asian dust contributes about 15 μg m^?3 of aerosol particles to northern Taiwan during winter monsoon, which accounts for about 24–30% of the PM₁₀ concentrations to the northern Taiwan. The contributions of Asian dust are raised pronouncedly to about 60–70% during major dust events. The impacts of Asian dust on Taiwan's air quality are most substantial in December. The Asian dust impacts decrease in other months, but still remain at around 30% in the late winter to early spring.     

The effect of coupled transport phenomena in the Opalinus Clay and implications for radionuclide transport

In this study, the potential effects of coupled transport phenomena on radionuclide transport in the vicinity of a repository for vitrified high-level radioactive waste (HLW) and spent nuclear fuel (SF) hosted by the Opalinus Clay in Switzerland, at times equal to or greater than the expected lifetime of the waste canisters (about 1000 years), are addressed. The solute fluxes associated with advection, chemical diffusion, thermal and chemical osmosis, hyperfiltration and thermal diffusion have been incorporated into a simple one-dimensional transport equation. The analytical solution of this equation, with appropriate parameters. shows that thermal osmosis is the only coupled transport mechanism that could, on its own, have a strong effect on repository performance. Based on the results from the analytical model, two-dimensional finite-difference models incorporating advection and thermal osmosis, and taking conservation of fluid mass into account, have been formulated. The results show that, under the conditions in the vicinity of the repository at the time scales of interest, and due to the constraints imposed by conservation of fluid mass, the advective component of flow will oppose and cancel the thermal-osmotic component. The overall conclusion is that coupled phenomena will only have a very minor impact on radionuclide transport in the Opalinus Clay, in terms of fluid and solute fluxes, at least under the conditions prevailing at times equal to or greater than the expected lifetime of the waste canisters (about 1000 years).     

Dual-domain solute transfer and transport processes: evaluation in batch and transport experiments

Soil macropore networks establish a dual-domain transport scenario in which water and solutes are preferentially channeled through soil macropores while slowly diffusing into and out of the bulk soil matrix. The influence of macropore networks on intra-ped solute diffusion and preferential transport in a soil typical of subsurface-drained croplands in the Midwestern United States was studied in batch- and column-scale experiments. In the batch diffusion studies with soil aggregates, the estimated diffusion radius (length) of the soil aggregates corresponded to the half-spacing of the aggregate fissures, suggesting that the intra-ped fissures reduced the diffusion impedance and preferentially allowed solutes to diffuse into the soil matrix. In the column-scale solute transport experiments, the average diffusion radius (estimated from HYDRUS-2D simulations and a first-order diffusive transfer term) was nearly double that of the batch-scale study. This increase may be attributed to a loss of pore continuity and a compounding of the small diffusion impedance through macropores at the larger scale. The column-scale solute transport experiments also suggest that two preferential networks exist in the soil. At and near soil saturation, a primary network of large macropores (possibly root channels and earthworm burrows) dominate advective transport, causing a high degree of physical and sorption nonequilibrium and simultaneous breakthrough of a nonreactive (bromide) and a reactive (alachlor) solute. As the saturation level decreases, the primary network drains, while transport through smaller macropores (possibly intra-ped features) continues, resulting in a reduced degree of nonequilibrium and separation in the breakthrough curves of bromide and alachlor.     

Gas phase transport in gravity sewers--A methodology for determination of horizontal gas transport and ventilation.

A method was developed for determination of horizontal gas transport and ventilation in gravity sewers. This was achieved by changing the composition of the sewer atmosphere by pulse injection of oxygen gas and subsequently measuring the oxygen concentration in a downstream manhole. Conventional tracer techniques may require sampling and may also affect the environment. The method developed is simple, based on direct monitoring and without environmental or toxic effects. The method was developed based on measurements in an intercepting gravity sewer. The horizontal gas transport processes were quantified by measuring the velocity and dispersion of the gas in the sewer atmosphere. Based on 54 measurements, the gas velocity was found to vary between 0.05 and 0.22 m/s. The coefficients of dispersion were calculated to be in the range 0.05 to 1.1 m2/s. Climatic conditions did not significantly influence the gas phase transport.     

Long-term scenarios for surface transport

Current trends in transport indicate that the system is moving away from a sustainable development (e.g., due to rising CO₂ emissions), and that major changes in technology and public policy, as well as in the behaviour of individuals, are necessary to make the transport system more compatible with environmental sustainability. This indicates the need to assess a set of future images for transport in relation to the environment. In this paper, reference and expert scenarios, which can act as a judgement framework for a sustainable transport system, are constructed on the basis of the recently developed 'spider model'. Based on a set of distinct characteristics of a transport system, represented by eight axes in the spatial, institutional, economic and social psychological field, an evaluation framework is constructed, which visualises the driving forces that largely influence the future of the transport system. There are several directions in which these factors may develop, and each of them will separately or in combination lead to entirely different transport systems. In this way, many scenarios can be constructed by connecting points on the successive axes. Such scenarios may range from market-oriented to regulatory pictures; the first may lead to a transport system in which individual, the second in which collective modes of transport dominate. Next, expected and desired scenarios are constructed from the opinions of Dutch transport experts, who responded to a nationwide survey. The expected scenario indicates that many current trends will continue, while the transport system is largely the same as the current one. The desired scenario, on the other hand, gives a more collective system, in which also many new modes are operating. The conclusion is that there are many roads for achieving a sustainable transport system, but that whichever one is chosen, it will be hard to follow.     

Colloid-facilitated transport of radionuclides through fractured media

The sorption of radionuclides on natural colloids may significantly modify their transport behaviour through fractured media, since radionuclides bound to colloids may not be subject to the important retardation mechanisms of matrix diffusion and sorption onto pore surfaces. This paper reports on theoretical and experimental work aimed at assessing the relevance of colloid-facilitated transport to repository safety analyses, with specific reference to the Swiss case. Transport models are presented, developed in conjunction with field- and laboratory-based studies of deep groundwater in the crystalline basement of northern Switzerland, in which colloid size distributions, compositions and sorption properties have been measured. Various potential mechanisms giving rise to both reversible and irreversible sorption are discussed. In the first case, a simple approach is examined which is based on previously reported models of colloid transport and assumes reversible, linear sorption on colloids, for which experimental data have been obtained. It is shown that transport of radionuclides would not, in general, be significantly enhanced because of this process. A more recently developed and more complex model is then described incorporating irreversible sorption, in which case the transport of radionuclides tends to be strongly dependent on the extent of colloid-fracture wall interaction.     

Energy-dependent transport of cadmium by Stichococcus bacillaris

The investigations showed that Stichococcus bacillaris took up cadmium by means of adsorption and energy-dependent transport. Cadmium transport into alga cells was inhibited in the dark as well as under the influence of the uncouplers: 2,4-dinitrophenol (DNP), carbonyl cyanide m-chlorophenylhydrazone (CCCP), and ATPase inhibitor - N, N′ - dicyclohexylcarbodiimide (DCCD). It has also been stated that CCCP caused ATP level decrease in cells.     

Environmental transport of endogenous dairy manure estrogens

Although estrogens originating from dairy manure applied to agricultural soils as a fertilizer can potentially contaminate surface water and groundwater, the variables that control transport are poorly understood. Our objective was to assess the potential for off-site movement of endogenous dairy cattle estrogens when manure is applied on fields at agronomically relevant fertilization rates. Estrone (E1), 17α-estradiol (α-E2), and 17β-estradiol (β-E2) were used in laboratory sorption, desorption, and transformation incubations with both manure and an agriculturally relevant soil. Sorption on manure containing 44% organic carbon exceeded sorption on soil containing 0.8% organic carbon by 20 to 150 times, following the pattern of β-E2 > α-E2 > E1. Approximately 20% of E1 and 17% of α-E2 were desorbed from manure, whereas only about 4% of β-E2 was desorbed. Thirty to seventy percent of α-E2 and β-E2 were converted to E1 in soil and manure, making it imperative that transformation reactions be considered when predicting transport and potential biological effects in the environment. Overall results indicate that high organic carbon concentrations and relatively low amounts of desorption inhibit the potential for off-site transport of endogenous dairy manure estrogens.     

Second-order modeling of arsenite transport in soils

Rate limited processes including kinetic adsorption-desorption can greatly impact the fate and behavior of toxic arsenic compounds in heterogeneous soils. In this study, miscible displacement column experiments were carried out to investigate the extent of reactivity during transport of arsenite in soils. Arsenite breakthrough curves (BTCs) of Olivier and Windsor soils exhibited strong retardation with diffusive effluent fronts followed by slow release or tailing during leaching. Such behavior is indicative of the dominance of kinetic retention reactions for arsenite transport in the soil columns. Sharp decrease or increase in arsenite concentration in response to flow interruptions (stop-flow) further verified that non-equilibrium conditions are dominant. After some 40-60 pore volumes of continued leaching, 30-70% of the applied arsenite was retained by the soil in the columns. Furthermore, continued arsenite slow release for months was evident by the high levels of residual arsenite concentrations observed during leaching. In contrast, arsenite transport in a reference sand material exhibited no retention where complete mass recovery in the effluent solution was attained. A second-order model (SOM) which accounts for equilibrium, reversible, and irreversible retention mechanisms was utilized to describe arsenite transport results from the soil columns. Based on inverse and predictive modeling results, the SOM model successfully depicted arsenite BTCs from several soil columns. Based on inverse and predictive modeling results, a second-order model which accounts for kinetic reversible and irreversible reactions is recommended for describing arsenite transport in soils.     

Sorption nonequilibrium during solute transport

The effects of pore-water velocity, solute hydrophobicity, and sorbent organic-carbon content on sorption nonequilibrium during solute transport were evaluated. Nonequilibrium transport was observed to increase with pore-water velocity, solute hydrophobicity, and sorbent organic-carbon content. Nonequilibrium transport of neutral organic compounds was not detected with low organic-carbon (TOC = 0.33 g kg⁻¹) aquifer material, but was detected on higher organic sorbents from the unsaturated zone (TOC = 2.6 g kg⁻¹) and the soil surface (TOC = 6.9 g kg⁻¹). For solute-sorbent combinations yielding retardation factors > 2, nonequilibrium during transport was observed. After experimentally accounting for slow solute diffusion in the aqueous phase and isotherm nonlinearity as potential contributors to nonequilibrium solute transport, sorption nonequilibrium was attributed to slow solute diffusion within the organic-carbon matrix.     

Synoptic-scale transport of ozone into Southern Ontario

This study focuses on synoptic-scale transport of ozone as it affects Southern Ontario. This process has been analyzed for the summer in 2001, as an example period of a frequent event that usually occurs during summer in this region. The work was carried out using the mesoscale modeling system generation 5 (MM5)/sparse matrix operator kernel emission modeling system (SMOKE)/community multiscale air quality (CMAQ) regional air quality modeling system, together with observational data from monitoring stations located throughout the modeling domain. Other different analyses have been carried out to supply more information apart from that obtained by the modeling system. A back-trajectory cluster methodology was used to evaluate the magnitude of the effects studied and an analysis of wind direction and cloud cover revealed a significant correlation with ozone concentration (R²=0.5–0.6). Synoptic sea-surface level pressure (SLP) patterns were also analyzed to examine other meteorological aspects. The contribution of natural background ozone to the total amount within the region was compared with that from synoptic-scale transport. The influence of emission of pollutants from selected areas on ozone concentrations in Southern Ontario was also analyzed. As relevant results of these analyses, the model predicts that background ozone is the largest contribution to the ground-level ozone concentration during days in which low values were recorded. However, when smog episodes occurred, the model predicts that around 60% of the ozone formed by anthropogenic emissions of pollutants is due to releases from nearby US states.     

Global-scale environmental transport of persistent organic pollutants

In order to realistically simulate both chemistry and transport of atmospheric organic pollutants, it is indispensable that the applied models explicitly include coupling between different components of the global environment such as atmosphere, hydrosphere, cryosphere and soil system. A model with such properties is presented.

The atmospheric part of the model is based on the equations in a general contravariant form which permits easy changes of the coordinate system by redefining the metric tensor of a specifically employed coordinate system. Considering a need to include explicitly the terrain effects, the terrain following spherical coordinate system is chosen from among many possible coordinate systems. This particular system is a combination of the Gal-Chen coordinates, commonly employed in mesoscale meteorological models, and the spherical coordinates, typical for global atmospheric models.

In addition to atmospheric transport, the model also simulates the exchange between air and different types of underlying surfaces such as water, soil, snow, and ice. This approach permits a realistic representation of absorption and delayed re-emission of pollutants from the surface to the atmosphere and, consequently, allows to capture hysteresis-like effects of the exchange between the atmosphere and the other components of the system. In this model, the most comprehensive numerical representation of the exchange is that for soil. In particular, the model includes a realistic soil module which simulates both diffusion and convection of a tracer driven by evaporation from the soil, precipitation, and gravity.

The model is applied to a long-term simulation of the transport of pesticides (hexachlorocyclohexanes in particular). Emission fluxes from the soil are rigorously computed on the basis of the realistic data of the agricultural application. All four modelled systems, i.e. atmosphere, soil, hydrosphere and cryosphere, are driven by objectively analysed meteorological data supplemented, when necessary, by climatological information. Therefore, the verification against the observed data is possible. The comparison of the model results and the observations taken at remote stations in the Arctic indicates that the presented global modelling system is able to capture both trends and short-term components in the observed time series of the concentrations, and therefore, provides a useful tool for the evaluation of the source–receptor relationships.     

Evaluating the sensitivity of a subsurface multicomponent reactive transport model with respect to transport and reaction parameters.

The input variables for a numerical model of reactive solute transport in groundwater include both transport parameters, such as hydraulic conductivity and infiltration, and reaction parameters that describe the important chemical and biological processes in the system. These parameters are subject to uncertainty due to measurement error and due to the spatial variability of properties in the subsurface environment. This paper compares the relative effects of uncertainty in the transport and reaction parameters on the results of a solute transport model. This question is addressed by comparing the magnitudes of the local sensitivity coefficients for transport and reaction parameters. General sensitivity equations are presented for transport parameters, reaction parameters, and the initial (background) concentrations in the problem domain. Parameter sensitivity coefficients are then calculated for an example problem in which uranium(VI) hydrolysis species are transported through a two-dimensional domain with a spatially variable pattern of surface complexation sites. In this example, the reaction model includes equilibrium speciation reactions and mass transfer-limited non-electrostatic surface complexation reactions. The set of parameters to which the model is most sensitive includes the initial concentration of one of the surface sites, the formation constant (Kf) of one of the surface complexes and the hydraulic conductivity within the reactive zone. For this example problem, the sensitivity analysis demonstrates that transport and reaction parameters are equally important in terms of how their variability affects the model results.     

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.     

Illuminating reactive microbial transport in saturated porous media: demonstration of a visualization method and conceptual transport model

We demonstrate a method to study reactive microbial transport in saturated translucent porous media using the bacteria Pseudomonas fluorescens 5RL genetically engineered to carry a plasmid with bioluminescence genes inducible by salicylate. Induced bacteria were injected into a cryolite grain filled chamber saturated with a sterile non-growth-promoting (phosphorus limited) chemical mixture containing salicylate as an aromatic hydrocarbon analogue. The amount of light produced by the bacteria serves as an estimator of the relative efficiency of aerobic biodegradation since bioluminescence is dependent on both salicylate and oxygen but only consumes oxygen. Bioluminescence was captured with a digital camera and analyzed to study the evolving spatial pattern of the bulk oxygen consuming reactions. As fluid flow transported the bacteria through the chamber, bioluminescence was observed to initially increase until an oxygen depletion zone developed behind the advective front. Bacterial transport was modeled with the advection dispersion equation and oxygen concentration was modeled assuming bacterial consumption via Monod kinetics with consideration of additional effects of rate-limited mass transfer from residual gas bubbles. Consistent with previous measurements, bioluminescence was considered proportional to oxygen consumed. Using the observed bioluminescence, model parameters were fit that were consistent with literature values and produced results in good agreement with the experimental data. These findings demonstrate potential for using this method to investigate the complex spatial and temporal dynamics of reactive microbial transport in saturated porous media.     

Modeling colloid transport for performance assessment

The natural system is expected to contribute to isolation at the proposed high-level nuclear waste (HLW) geologic repository at Yucca Mountain, NV (YM). In developing performance assessment (PA) computer models to simulate long-term behavior at YM, colloidal transport of radionuclides has been proposed as a critical factor because of the possible reduced interaction with the geologic media. Site-specific information on the chemistry and natural colloid concentration of saturated zone groundwaters in the vicinity of YM is combined with a surface complexation sorption model to evaluate the impact of natural colloids on calculated retardation factors (RF) for several radioelements of concern in PA. Inclusion of colloids into the conceptual model can reduce the calculated effective retardation significantly. Strongly sorbed radionuclides such as americium and thorium are most affected by pseudocolloid formation and transport, with a potential reduction in RF of several orders of magnitude. Radioelements that are less strongly sorbed under YM conditions, such as uranium and neptunium, are not affected significantly by colloid transport, and transport of plutonium in the valence state is only moderately enhanced. Model results showed no increase in the peak mean annual total effective dose equivalent (TEDE) within a compliance period of 10,000 years, although this is strongly dependent on container life in the base case scenario. At longer times, simulated container failures increase and the TEDE from the colloidal models increased by a factor of 60 from the base case. By using mechanistic models and sensitivity analyses to determine what parameters and transport processes affect the TEDE, colloidal transport in future versions of the TPA code can be represented more accurately.     

Acid rain: chemistry and transport

This review describes the more important features of the emission, chemistry, transport and deposition of pollutants involved in acid deposition. Global emissions, both natural and man-made, of sulphur and nitrogen oxides are discussed and examples of spatial distributions and trends over the last century presented. The more significant chemical and physical processes involved in the transformation of the primary emissions into their acidic end products are described, including a summary of the approximate timescales of the processes involved. Measurements and modelled calculations of spatial and temporal patterns in the deposition of acidic pollutants by both wet and dry pathways are presented.