A high-resolution non-invasive approach to quantify oxygen transport across the capillary fringe and within the underlying groundwater

Oxygen transport across the capillary fringe is relevant for many biogeochemical processes. We present a non-invasive technique, based on optode technology, to measure high-resolution concentration profiles of oxygen across the unsaturated/saturated interface. By conducting a series of quasi two-dimensional flow-through laboratory experiments, we show that vertical hydrodynamic dispersion in the water-saturated part of the capillary fringe is the process limiting the mass transfer of oxygen. A number of experimental conditions were tested in order to investigate the influence of grain size and horizontal flow velocity on transverse vertical dispersion in the capillary fringe. In the same setup, analogous experiments were simultaneously carried out in the fully water-saturated zone, therefore allowing a direct comparison with oxygen transfer across the capillary fringe. The outcomes of the experiments under various conditions show that oxygen transport in the two zones of interest (i.e., the unsaturated/saturated interface and the saturated zone) is characterized by very similar transverse dispersion coefficients. An influence of the capillary fringe morphology on oxygen transport has not been observed. These results may be explained by the narrow grain size distribution used in the experiments, leading to a steep decline in water saturation at the unsaturated/saturated interface and to the absence of trapped gas in this transition zone. We also modeled flow (applying the van Genuchten and the Brooks-Corey relationships) and two-dimensional transport across the capillary fringe, obtaining simulated profiles of equivalent aqueous oxygen concentration that were in good agreement with the observations.     

Aerobic and anaerobic transformations of sulfide in a sewer system--field study and model simulations.

The formation and fate of sulfide in a force main and a downstream-located gravity sewer were investigated in an extensive field study. Sulfide formation in the force main was significant. However, during 14 minutes of transport in the gravity sewer, the sulfide concentration decreased 30%, on average. An application of a conceptual sewer process model for simulating the formation and fate of sulfide was demonstrated. Overall, the model predicted that approximately 90% of the decrease of the sulfide concentration in the gravity sewer was the result of sulfide oxidation and that only a small fraction entered the sewer atmosphere, causing odor and corrosion. Even so, the model predicted concrete corrosion rates of up to 1.2 mm/y in the gravity sewer section.     

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.     

Simulation of solute transport across low-permeability barrier walls

Low-permeability, non-reactive barrier walls are often used to contain contaminants in an aquifer. Rates of solute transport through such barriers are typically many orders of magnitude slower than rates through the aquifer. Nevertheless, the success of remedial actions may be sensitive to these low rates of transport. Two numerical simulation methods for representing low-permeability barriers in a finite-difference groundwater-flow and transport model were tested. In the first method, the hydraulic properties of the barrier were represented directly on grid cells and in the second method, the intercell hydraulic-conductance values were adjusted to approximate the reduction in horizontal flow, allowing use of a coarser and computationally efficient grid. The alternative methods were tested and evaluated on the basis of hypothetical test problems and a field case involving tetrachloroethylene (PCE) contamination at a Superfund site in New Hampshire. For all cases, advective transport across the barrier was negligible, but preexisting numerical approaches to calculate dispersion yielded dispersive fluxes that were greater than expected. A transport model (MODFLOW-GWT) was modified to (1) allow different dispersive and diffusive properties to be assigned to the barrier than the adjacent aquifer and (2) more accurately calculate dispersion from concentration gradients and solute fluxes near barriers. The new approach yields reasonable and accurate concentrations for the test cases.     

Modeling the formation and fate of odorous substances in collection systems.

A conceptual model that simulates the formation and fate of odorous substances in branched collection systems is presented. The model predicts the activity of the relevant biomass phenotypes under aerobic, anoxic, and anaerobic conditions in force mains and gravity sewers. The formation and fate of individual, malodorous substances in the bulk water, biofilms, and sediments are modeled. The release of odorous compounds from the bulk water to the sewer gas phase, their fate in the gas phase, and their subsequent release into the urban atmosphere is simulated. Examples of model application include the prediction of hydrogen sulfide and malodorous fermentation products from force mains and gravity sewers.     

A different approach for predicting reaeration rates in gravity sewers and completely mixed tanks.

A new semiempirical approach is presented for predicting air-to-water oxygen transfer rates in mixed tanks and gravity sewers, using methods adopted from mixing theory. First, a flocculation unit was used to impart selected mean velocity gradients (G) into a completely mixed tank, from which oxygen was first removed, and dissolved oxygen concentrations were measured with time. Regression analysis was used to fit the rate of oxygen transfer equation against G. The reaeration rate in completely mixed reactors was found to be proportional to G2 (R2 = 0.987). Subsequently, G was linked to headloss in sewers, and the equation was calibrated using a slope-adjustable, 27-m-long, gravity-flow, experimental sewer (internal diameter, D = 0.16 m). Here, the reaeration rate was proportional to G1 (R2 = 0.981). The equation was compared with existing oxygen transfer models and validated against experimental data from the literature, to which the overall mass transfer coefficient for oxygen, K(L)a, derived by the new approach, conformed well.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

Experimental analysis of particle concentration heterogeneity in a ventilated scale chamber

The mixing processes of the aerosol particles from an outdoor environment in a ventilated scale chamber were experimentally studied. The particles were classified into five groups by size: 0.3–0.5 μm, 0.5–1.0 μm, 1.0–3.0 μm, 3.0–5.0 μm and 5.0–10.0 μm. The developing process for the concentration of each particle group was measured in different kinds of flow fields.The results show that the flow field configuration can effectively influence the dispersion time rate of the particles at certain positions. The increase in particle diameter can decrease the dispersion time rate. When the gas flow velocity is high, the particle dispersion time rate is independent of particle size; but when the gas flow velocity is low, particle size can significantly affect the particle dispersion time rate because the turbulent diffusion becomes important in the air and particle transport. The uniformity of the particle concentration for certain positions in steady state tends to be controlled by the inflow velocity, flow field configuration and the particle diameters.     

A large-scale experiment on mass transfer of trichloroethylene from the unsaturated zone of a sandy aquifer to its interfaces

A large-scale experiment was conducted to investigate the transport of trichloroethylene (TCE) vapors in the unsaturated zone and to determine the mass transfer to the groundwater and the atmosphere. The experiment involved injection of 5 1 of TCE in the unsaturated zone under controlled conditions, with multidepth sampling of gas and water through the unsaturated zone and across the capillary zone into underlying groundwater. The mass transfer of TCE vapors from the vadose zone to the atmosphere was quantified using a vertical flux chamber. A special soil water sampler was used to monitor transport across the capillary fringe. Experimental data indicated that TCE in the unsaturated zone was mainly transported to the atmosphere and this exchange reduced significantly the potential for groundwater pollution. The maximum measured TCE flux to the atmosphere was about 3 g/m(2)/day. Observed and calculated fluxes based on vertical TCE vapor concentration gradients and Fick's law were in good agreement. This confirms that TCE vapor transport under the experimental conditions was governed essentially by molecular diffusion. TCE vapors also caused a lower, but significant contamination of the underlying groundwater by dispersion across the capillary fringe with a corresponding maximum flux of about 0.1 g/m(2)/day. This mass transfer to groundwater is partly uncertain due to an inadvertent entry of some nonaqueous phase liquid (NAPL) from the source area into the saturated zone. Application of an analytical solution to estimate the TCE flux from the unsaturated zone to the groundwater indicated that this phenomenon is not only influenced by molecular diffusion but also by vertical dispersion. The mass balance indicates that, under the given experimental conditions (e.g. proximity of the source emplacement relative to the soil surface, relatively high permeable porous medium), nearly 95% of the initial TCE mass was transferred to the atmosphere.     

Simulation of boundary layer trajectory dispersion sensitivity to soil moisture conditions: MM5 and Noah-based investigation

In this study, the sensitivity of trajectory paths to anomalous soil moisture was analyzed during three different synoptic episodes in June 2006. The MM5 and Noah land surface models were used to simulate the response of the planetary boundary layer. The HYSPLIT model was used for trajectory analysis. It was found that the response in horizontal lower-level wind field was larger at regions where vertical wind velocity changes were also large. In addition, the sensitivity to soil moisture changes was significant and localized where convective activity was well developed and synoptic effects did not dominate. A non-local effect was felt over the rest of the domain where convection was not present since the model atmosphere reacted as a whole to compensate for induced changes in vertical velocity. This finding was supported by the fact that domain averaged vertical velocities changes were of the order of 0.2 cm s^?1 or less at about 650 hPa and about 200 times smaller than modeled local vertical velocity changes. The largest change in horizontal wind field near the surface was found for weak synoptic events on June 11–12 and June 22–23 while the stronger synoptic event of June 17–18 showed smaller differences. These changes in wind field conditions impacted the transport and dispersion of pollutants. To quantify the sensitivity of air quality estimates to soil moisture uncertainty, we have used three well known measures of trajectory differences: the absolute horizontal transport deviation (AHTD), the relative horizontal transport deviation (RHTD) and the absolute vertical transport deviation (AVTD) for an ensemble of 98 trajectories departing from a region well within the computational domain. For the June 11–12 event it was found that for wet and dry soil moisture experiments, AHTD, RHTD, and AVGTD can reach values in the range 60–100 km, 10–20% and 500–900 m at 24 h run time, respectively. For the June 17–18 and June 22–23 events these values of trajectory differences were reduced more than half. These differences in behavior between time periods are largely attributed to the combined effects of synoptic forcing and the sensitivity of planetary boundary layer to soil moisture changes during well developed convection. The implication for air quality studies is that the soil moisture anomaly and related uncertainty in planetary boundary layer response needs to be incorporated in order to construct an ensemble of the most probable scenarios in which pollutants are released and transported throughout a given target region.     

Characteristics of nitrate horizontal transport in a paddy field of the Tai Lake region,China

The characteristics of nitrate horizontal transport in a major paddy soil, Wu Shan soil in the Tai Lake region, were studied. The concentration of nitrate during horizontal movement decreased with the increasing in distance from the tracer source, the change following a logarithmic function. The concentration of the nitrate was strongly correlated with the soil moisture content, as an exponential function. The horizontal transport velocity of nitrate was significantly correlated with the distance of the tracer source as power function. Therefore, the velocity of nitrate horizontal transport was controlled by the concentration gradient of nitrate, and soil water potential gradient from beginning to the 20 cm mark in the horizontal column. However, the velocity of nitrate horizontal was stable beyond 20 cm, where it was controlled by soil matric potential.     

The influence of wind speed on airflow and fine particle transport within different building layouts of an industrial city

In atmospheric environment, the layout difference of urban buildings has a powerful influence on accelerating or inhibiting the dispersion of particle matters (PM). In industrial cities, buildings of variable heights can obstruct the diffusion of PM from industrial stacks. In this study, PM dispersed within building groups was simulated by Reynolds-averaged Navier-Stokes equations coupled Lagrangian approach. Four typical street building arrangements were used: (a) a low-rise building block with Height/base H/b = 1 (b = 20 m); (b) step-up building layout (H/b = 1, 2, 3, 4); (c) step-down building layout (H/b = 4, 3, 2, 1); (d) high-rise building block (H/b = 5). Profiles of stream functions and turbulence intensity were used to examine the effect of various building layouts on atmospheric airflow. Here, concepts of particle suspension fraction and concentration distribution were used to evaluate the effect of wind speed on fine particle transport. These parameters showed that step-up building layouts accelerated top airflow and diffused more particles into street canyons, likely having adverse effects on resident health. In renewal old industry areas, the step-down building arrangement which can hinder PM dispersion from high-level stacks should be constructed preferentially. High turbulent intensity results in formation of a strong vortex that hinders particles into the street canyons. It is found that an increase in wind speed enhanced particle transport and reduced local particle concentrations, however, it did not affect the relative location of high particle concentration zones, which are related to building height and layout.

Implications: This study has demonstrated the height variation and layout of urban architecture affect the local concentration distribution of particulate matter (PM) in the atmosphere and for the first time that wind velocity has particular effects on PM transport in various building groups. The findings may have general implications in optimization the building layout based on particle transport characteristics during the renewal of industrial cities. For city planners, the results and conclusions are useful for improving the local air quality. The study method also can be used to calculate the explosion risk of industrial dust for people who live in industrial cities.     

Reaction-enhanced mass transfer and transport from non-aqueous phase liquid source zones

The objective of the following research is to theoretically quantify the enhancement of interphase mass transfer of dissolved non-aqueous phase liquid (NAPL) compounds from the non-aqueous phase to the aqueous phase and the enhancement of dispersive mass transport from a NAPL zone due to destruction of dissolved NAPL compounds. For relatively slow reaction rates, such as for permanganate and perchloroethene (PCE), local-scale mass transfer enhancement is expected to be small. Dispersive mass transport with reaction from a horizontal NAPL zone can be quantified using equations derived for a mathematically equivalent falling film reactor system. In contrast to local-scale interphase mass transfer, dispersive mass transport from NAPL zones may be significantly increased by reaction. Enhancement factors due to destruction of the NAPL compound(s) are mainly dependent on NAPL solubility and oxidant concentration and to a lesser extent on reaction rate, stoichiometry, and transverse dispersion coefficients. Higher NAPL solubility and/or lower oxidant concentration reduces the maximum expected enhancement factor. Reaction enhancement factors for mass transport from NAPL zones are expected to be in the range of 5-50 for permanganate and chlorinated solvents. Theoretical results suggest that assuming instantaneous reaction rates may be appropriate for dispersive mass transport from NAPL zones.     

Simulation and analysis of solute transport in 2D fracture/pipe networks: the SOLFRAC program

The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL (labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm). It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments performed in single fractures.     

Effective pollutant emission heights for atmospheric transport modelling based on real-world information

Emission data needed as input for the operation of atmospheric models should not only be spatially and temporally resolved. Another important feature is the effective emission height which significantly influences modelled concentration values. Unfortunately this information, which is especially relevant for large point sources, is usually not available and simple assumptions are often used in atmospheric models. As a contribution to improve knowledge on emission heights this paper provides typical default values for the driving parameters stack height and flue gas temperature, velocity and flow rate for different industrial sources. The results were derived from an analysis of the probably most comprehensive database of real-world stack information existing in Europe based on German industrial data. A bottom-up calculation of effective emission heights applying equations used for Gaussian dispersion models shows significant differences depending on source and air pollutant and compared to approaches currently used for atmospheric transport modelling.     

High-resolution simulation of volcanic sulfur dioxide dispersion over the Miyake Island

After severe eruptions of the volcano at Miyake Island in August 2000, a large amount of volcanic gas was released into the atmosphere. To simulate flows and dispersion of sulfur dioxide (SO₂) over Miyake Island, a set of numerical models was developed. The multi-nesting method was adopted to reflect a realistic meteorological field and to sufficiently resolve the flow over the island with a diameter of 8 km. The outermost model was the Regional Spectral Model (RSM) of the Japan Meteorological Agency (JMA) with a horizontal grid size of 10 km. Finer atmospheric structure was simulated with the nonhydrostatic model jointly developed by the Meteorological Research Institute and the Numerical Prediction Division of JMA (MRI/NPD-NHM) with grid intervals of 2 km, 400 m and 100 m. Realistic topography of the island was represented in the innermost model. The Lagrangian particle method was applied to the dispersion model, which is driven by the meteorological field of the 100 m grid MRI/NPD-NHM. The random walk procedure was used to represent the turbulent diffusion. The model was verified in four cases. Simulated SO₂ concentrations agreed well with observed concentrations at a monitoring station including temporal variation. Under a large synoptic change, however, accurate prediction became difficult. Further numerical experiments have been done to investigate characteristics of the flow and the distribution of SO₂. Steady inflows, classified according to the surface wind speed and direction, were assumed. Simulated SO₂ distribution on the ground apparently depends on the surface wind. Under relatively weak inflow, there is a large diurnal change in SO₂ distribution, affected by the thermally induced flow. SO₂ gas is widely spread downstream in the nighttime but hardly reaches the coastal area in the daytime. On the other hand, SO₂ gas steadily reached the downstream coast with little diurnal variation under the stronger inflow. Ground temperature, as well as the static stability of the inflow, also influences downstream wind, turbulent diffusivity and SO₂ distribution.     

A tracer study of headspace ventilation in a collector sewer

A field-scale tracer test was conducted to evaluate in-situ ventilation rates in a major collector sewer. The sewer under study was approximately 11 km long and ranged from 0.61 to 2.1 m in diameter. For the purposes of the tracer testing, the collector was divided into four reaches, each of which was tested individually. The tracer test involved injecting a measured volume of CO gas into a manhole over a short time period. CO concentrations were then measured in the collector headspace at selected manholes along the length of the reach. The technique employed successfully measured average headspace velocities over extended lengths of the collector. In a section that had a relatively stagnant headspace, approximately 1.1 km of sewer could be evaluated, with substantial tracer loss attributed to losses to manholes. In a section of the sewer with elevated headspace velocities, a section approximately 7.0 km long was successfully tested with one injection of tracer gas. The velocities observed in the collector varied substantially with time and location in the collector. The lowest velocities measured were in the upstream sections, with a minimum observed value of 3.8 m/min. The highest velocities were observed in the downstream sections, with a maximum value of 31.5 m/min. The presence of a substantial drop structure appeared to reduce the headspace velocity in the upstream reach. In general, there was an increasing trend in gas-phase flows with distance along the length of the collector. Flows at the discharge end of the collector were almost 2 orders of magnitude greater than those at the beginning.     

A Tracer Study of Headspace Ventilation in a Collector Sewer

ABSTRACT

A field-scale tracer test was conducted to evaluate in-situ ventilation rates in a major collector sewer. The sewer under study was ~11 km long and ranged from 0.61 to 2.1 m in diameter. For the purposes of the tracer testing, the collector was divided into four reaches, each of which was tested individually. The tracer test involved injecting a measured volume of CO gas into a manhole over a short time period. CO concentrations were then measured in the collector headspace at selected manholes along the length of the reach.

The technique employed successfully measured average headspace velocities over extended lengths of the collector. In a section that had a relatively stagnant headspace, ~1.1 km of sewer could be evaluated, with substantial tracer loss attributed to losses to manholes. In a section of the sewer with elevated headspace velocities, a section ~7.0 km long was successfully tested with one injection of tracer gas. The velocities observed in the collector varied substantially with time and location in the collector. The lowest velocities measured were in the upstream sections, with a minimum observed value of 3.8 m/min. The highest velocities were observed in the downstream sections, with a maximum value of 31.5 m/min. The presence of a substantial drop structure appeared to reduce the headspace velocity in the upstream reach. In general, there was an increasing trend in gas-phase flows with distance along the length of the collector. Flows at the discharge end of the collector were almost 2 orders of magnitude greater than those at the beginning.     

Stochastic modeling of chemical oxygen demand transformations in gravity sewers.

Wastewater quality characteristics in terms of biomass, its substrates, and the corresponding kinetic and stoichiometric parameters were determined based on 109 wastewater samples originating from five different campaigns in four different sewer networks. Quality parameters were determined by model calibration of measured wastewater oxygen uptake rates applying a model that describes the aerobic breakdown of wastewater organic matter. Thereafter, the distributions of the parameters were analyzed. Two of the five datasets were obtained at the upstream end of a five-km-long, intercepting gravity sewer. For each of these upstream wastewater samples, downstream samples were collected with a delay corresponding to the residence time. The upstream distributions of the wastewater composition were used as boundary conditions for a Monte Carlo simulation. The calculated downstream distributions were compared to the measured downstream distributions and good agreement was observed.     

An evaluation of flare combustion efficiency using open-path Fourier transform infrared technology

Open-path Fourier transform infrared (OP-FTIR) technology was used to evaluate the combustion efficiency of a flare for comparison to several combustion models. Most flares have been considered an effective method for destroying organic compounds and anything that burns. As the Btu content of the flare gas is reduced, the combustion efficiency may also be reduced. Recent studies have suggested that lower Btu flares may have efficiencies as low as 65%. In addition, models have been developed to predict the effect of wind speed and stack discharge velocity on the combustion efficiency. This study was conducted on a low-Btu flare gas that is primarily CO. While the models would predict efficiencies as low as 30%, the sampling using OP-FTIR showed most combustion efficiencies well above 90%. Three methods were used to track combustion efficiency: monitoring the ratio of CO to CO2, monitoring the ratio of CO to tracer gas, and dispersion modeling. This study was complicated by the presence of two flare stacks, thus two tracer gases were used--SF6 and CF4. A method was developed for distinguishing between the two stacks and quantifying the efficiency in each stack.