Contaminant transport in dual-porosity media with dissolved organic matter and bacteria present as mobile colloids

In riverbank filtration, contaminant transport is affected by colloidal particles such as dissolved organic matter (DOM) and bacterial particles. In addition, the subsurface heterogeneity influences the behavior of contaminant transport in riverbank filtration. A mathematical model is developed to describe the contaminant transport in dual-porosity media in the presence of DOM and bacteria as mobile colloids. In the model development, a porous medium is divided into the mobile and immobile regions to consider the presence of ineffective micropores in physically heterogeneous riverbanks. We assume that the contaminant transport in the mobile region is controlled by the advection and dispersion while the contaminant transport in the immobile region occurs due to the molecular diffusion. The contaminant transfer between the mobile and immobile regions takes place by diffusive mass transfer. The mobile region is conceptualized as a four-phase system: two mobile colloidal phases, an aqueous phase, and a solid matrix. The complete set of governing equations is solved numerically with a fully implicit finite difference method. The model results show that in riverbank filtration, the contaminant can migrate further than expected due to the presence of DOM and bacteria. In addition, the contaminant mobility increases further in the presence of the immobile region in aquifers. A sensitivity analysis shows that in dual-porosity media, earlier breakthrough of the contaminant takes place as the volumetric fraction of the mobile region decreases. It is also demonstrated that as the contaminant mass transfer rate coefficient between the mobile and immobile regions increases, the contaminant concentration gradient between the two regions reverses at earlier pore volumes. The contaminant mass transfer coefficient between the mobile and immobile regions mainly controls the tailing effect of the contaminant breakthrough. The contaminant breakthrough curves are sensitive to changes in contaminant adsorption and desorption rate coefficients on DOM and bacteria. In situations where the contaminant is released in the presence of DOM and bacteria in dual-porosity media, the early breakthrough and tailing occur due to the colloidal facilitation and presence of immobile regions.     

Transport of a non-volatile contaminant in unsaturated porous media in the presence of colloids

Stable colloidal particles can travel long distances in subsurface environments and carry particle-reactive contaminants with them to locations further than predicted by the conventional advective-dispersive transport equation. When such carriers exist in a saturated porous medium, the system can be idealized as consisting of three phases: an aqueous phase, a carrier phase, and a stationary solid matrix phase. However, when colloids are present in an unsaturated porous medium, the system representation should include one more phase, i.e. the air phase. In the work reported, a mathematical model was developed to describe the transport and fate of the colloidal particles and a non-volatile contaminant in unsaturated porous media. The model is based on mass balance equations in a four-phase porous medium. Colloid mass transfer mechanisms among aqueous, solid matrix, and air phases, and contaminant mass transfer between aqueous and colloid phases are represented by kinetic expressions. Governing equations are non-dimensionalized and solved to investigate colloid and contaminant transport in an unsaturated porous medium. A sensitivity analysis of the transport model was utilized to assess the effects of several parameters on model behavior. The colloid transport model matches successfully with experimental data of Wan and Wilson. The presence of air-water interface retards the colloid transport significantly counterbalancing the facilitating effect of colloids. However, the retardation of contaminant transport by colloids is highly dependent on the properties of the contaminant and the colloidal surface.     

An independent prediction of the effect of dissolved organic matter on the transport of polycyclic aromatic hydrocarbons

The transport of polycyclic aromatic hydrocarbons (PAH) in porous media in the presence of dissolved organic matter (DOM) was predicted with a transport bicontinuum model using independently obtained relationships to derive transport parameters for describing the effect of PAH binding to the DOM. The sorption constants of PAHs to soil and their binding constants to DOM were derived from basic correlations with K(ow) (indicator of hydrophobicity). The kinetic (rate) constants were derived from previously published correlations with K(p) (sorption constant). The independently obtained sorption and rate constants were corrected for binding to DOM and were used to predict the breakthrough curves (BTC) of contaminants in the presence and the absence of DOM. Column results confirmed the independently predicted BTC of PAHs in the presence of DOM that did not sorb to the solid phase, as well as the effect of DOM on the rate of the sorption and desorption processes. These findings confirm the ability to quantitatively describe how DOM facilitates transport of contaminants in the subsurface using independently derived parameters.     

Transport of reactive colloids and contaminants in groundwater: effect of nonlinear kinetic interactions

Transport of reactive colloids in groundwater may enhance the transport of contaminants in groundwater. Often, the interpretation of results of transport experiments is not a simple task as both reactions of colloids with the solid matrix and reactions of contaminants with the solid matrix and mobile and immobile colloids may be time dependent and nonlinear. Further colloid transport properties may differ from solute transport properties. In this paper, a one-dimensional model for coupled and contaminant in a porous medium (COLTRAP) is presented together with simulation results. Calculated breakthrough curves (BTC's) during contamination and decontamination show systematically the effect of nonlinear and kinetic interactions on contaminant transport in the presence of reactive colloids, and the effect of colloid transport properties that differ from solute transport properties. It is shown that in case of linear kinetic reactions, the rate of exchange of mobile and immobile colloids have a large impact on the shape of BTC's even if the solid matrix is saturated with respect to colloids. BTC's during the contamination and decontamination phase have identical shapes in this case. Moreover, the slow reactions of contaminants and colloids may lead to unretarded breakthrough of contaminants. Independent of reaction rates, nonlinear reactions lead to BTC's that are steeper during contamination than in the linear case. A characteristic aspect of nonlinear sorption is that shapes of BTC's differ during the contamination and decontamination phase. It has been observed that shapes of some of the simulated adsorption and desorption curves are similar as shapes found in experiments reported in literature. This stresses the importance of incorporating both kinetics and nonlinearity in models for coupled colloid and contaminant transport and the capability of COLTRAP to interpret experimental results. Finally, to figure out whether nonlinear processes play a role, it is very important to consider both contamination and decontamination in transport experiments.     

Contaminant degradation in physically and chemically heterogeneous aquifers

This paper examines the importance of the correlation between hydraulic conductivity (K) and degradation rate constant (k) during the transport of reactive contaminants in heterogeneous aquifers. We simulated reactive transport in an ensemble of two-dimensional heterogeneous aquifers. Two sets of transport simulations were conducted: one in which a perfect positive correlation was assumed between ln(K) and ln(k), and one in which a perfect negative correlation was assumed. We found that the sign of the correlation has important consequences for the contaminant transport. Qualitatively, a negative correlation leads to significantly more pronounced "fingering" of the contaminant plume than does a positive correlation, with potentially important consequences for downgradient receptors. Quantitatively, the expected behavior (as quantified by the contaminant mass remaining in the aquifer) is statistically different between the positive and negative cases: on average, more contaminant mass persists when K and k are negatively correlated. Also, the negative correlation leads to more variability between realizations of the ensemble, whereas a positive correlation induces relatively little variability between realizations. We discuss the implications of these findings for the management of contaminated aquifers.     

Modeling of flow and contaminant transport in coupled stream-aquifer systems

This paper describes an integrated approach for modeling flow and contaminant transport in hydraulically connected stream-aquifer systems. The code, FTSTREAM, extended the capabilities of the ground-water model, FTWORK, to incorporate chemical fate and transport in streams. Flow in the stream network is modeled as an unsteady, spatially varying flow, while transport modeling is based on a one-dimensional advection-dispersion equation. In addition to sorption and decay during transport in ground water, the model incorporates volatilization, settling and decay during transport in surface water. The interaction between surface water and ground water is accommodated by a leakage term and is implemented in the model using an iterative Picard-type procedure to ensure mass conservation across the interface between the two systems. The modeling approach is used to simulate contaminant transport in the Mad River, Ohio, which is hydraulically connected to a buried valley aquifer of sand and gravel outwash. The river is a receiving stream in the upstream part of the modeled area. Downstream, heavy pumping from a municipal well field causes the river to become a loosing stream. Induced infiltration from the river is responsible for a considerable portion of the well yield. The flow and transport model, developed for this study, simulates coupling between flow in the aquifer and the river. Hypothetical sources of contamination are introduced at selected locations in the upstream portion of the aquifer. The model is then used to simulate the expected transport in both the aquifer and the stream. A series of simulations elucidates the role of the river in facilitating the transport of the hypothetical contaminants in ground water and surface water. Effect of sorption, retardation and volatilization on contaminant transport is also examined for the case of the volatile organic compounds.     

MTBE Transport Screening Analysis in Vadose Zone Risk and Liability Apportionment Assessment

The problem of allocating liability cleanup costs is an arduous task when more than one potentially responsible party has contributed to the groundwater plume. This problem is most likely to be encountered when dealing with methyl-tert -butyl-ether (MTBE) contamination, as MTBE is seen to travel large distances in underlying aquifers. There has been a significant effort in the recent past to develop liability allocation methodologies that incorporate fate and transport behavior and toxicological characteristics of the contaminants. The application of such methods often requires estimation of contaminant input from the vadose zone into the underlying aquifer. A screening level analysis is presented here to develop preliminary insights on relative mass contributions arising from different source types. The analysis illustrates how different vadose zone conceptualizations lead to vastly different contaminant loadings. Parametric studies indicate that the contaminant flux into the aquifer is very sensitive to changes in water infiltration rates. Hence, a reliable estimate of this parameter is critical for equitable allocation of remedial costs. Conceptual model formulation, should focus on identifying whether the fuel present in the aquifer can flow as a separate phase. It is also important to obtain reliable estimates for fluid saturations at the site.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

MTBE Transport Screening Analysis in Vadose Zone Risk and Liability Apportionment Assessment

The problem of allocating liability cleanup costs is an arduous task when more than one potentially responsible party has contributed to the groundwater plume. This problem is most likely to be encountered when dealing with methyl- tert -butyl-ether (MTBE) contamination, as MTBE is seen to travel large distances in underlying aquifers. There has been a signi®cant effort in the recent past to develop liability allocation methodologies that incorporate fate and transport behavior and toxicological characteristics of the contaminants. The application of such methods often requires estimation of contaminant input from the vadose zone into the underlying aquifer. A screening level analysis is presented here to develop preliminary insights on relative mass contributions arising from different source types. The analysis illustrates how different vadose zone conceptualizations lead to vastly different contaminant loadings. Parametric studies indicate that the contaminant flux into the aquifer is very sensitive to changes in water infiltration rates. Hence, a reliable estimate of this parameter is critical for equitable allocation of remedial costs. Conceptual model formulation, should focus on identifying whether the fuel present in the aquifer can flow as a separate phase. It is also important to obtain reliable estimates for fluid saturations at the site.     

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

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

Anthropogenic contaminants as tracers in an urbanizing karst aquifer

Karst aquifers are uniquely vulnerable to contamination. In the Barton Springs segment of the karstic Edwards aquifer (Texas, U.S.A.), urban contaminants such as pesticides and volatile organic compounds frequently are detected in spring base flow. To determine whether contaminant concentrations change in response to storms, and if they therefore might act as tracers of focused recharge, samples were collected from Barton Springs at closely spaced intervals following three storms. Two herbicides (atrazine and simazine), two insecticides (carbaryl and diazinon), and a solvent (tetrachloroethene) described breakthrough curves over a 1-week period following one or more storms. The breakthrough curves were decomposed into two to five log-normal subcurves, which were interpreted as representing pulses of contaminants moving through the aquifer. Each subcurve could be used in the same way as an artificial tracer to determine travel time to and recovery at the spring. The contaminants have several advantages over artificial tracers: they represent the actual compounds of interest, they are injected essentially simultaneously at several points, and they are injected under those conditions when transport is of the most interest, i.e., following storms. The response of storm discharge, specific conductance, and contaminant loading at the spring depended on initial aquifer flow conditions, which varied from very low (spring discharge of 0.48 m3/s) to high (spring discharge of 2.7 m3/s): concentrations and recovery were the highest when initial aquifer flow conditions were low. This behavior provides information about aquifer structure and the influence of aquifer flow condition on transport properties.     

Processes controlling the distribution and natural attenuation of dissolved phenolic compounds in a deep sandstone aquifer

Processes controlling the distribution and natural attenuation (NA) of phenol, cresols and xylenols released from a former coal-tar distillation plant in a deep Triassic sandstone aquifer are evaluated from vertical profiles along the plume centerline at 130 and 350 m from the site. Up to four groups of contaminants (phenols, mineral acids, NaOH, NaCl) form discrete and overlapping plumes in the aquifer. Their distribution reflects changing source history with releases of contaminants from different locations. Organic contaminant distribution in the aquifer is determined more by site source history than degradation. Contaminant degradation at total organic carbon (TOC) concentrations up to 6500 mg l(-1) (7500 mg l(-1) total phenolics) is occurring by aerobic respiration NO3-reduction, Mn(IV)-/Fe(III)-reduction, SO4-reduction, methanogenesis and fermentation, with the accumulation of inorganic carbon, organic metabolites (4-hydroxybenzaldehyde, 4-hydroxybenzoic acid), acetate, Mn(II), Fe(II), S(-II), CH4 and H2 in the plume. Aerobic and NO3-reducing processes are restricted to a 2-m-thick plume fringe but Mn(IV)-/Fe(II)-reduction, SO4-reduction, methanogenesis and fermentation occur concomitantly in the plume. Dissolved H2 concentrations in the plume vary from 0.7 to 110 nM and acetate concentrations reach 200 mg l(-1). The occurrence of a mixed redox system and concomitant terminal electron accepting processes (TEAPs) could be explained with a partial equilibrium model based on the potential in situ free energy (deltaGr) yield for oxidation of H2 by specific TEAPs. Respiratory processes rather than fermentation are rate limiting in determining the distribution of H2 and TEAPs and H2 dynamics in this system. Most (min. 90%) contaminant degradation has occurred by aerobic and NO3-reducing processes at the plume fringe. This potential is determined by the supply of aqueous O2 and NO3 from uncontaminated groundwater, as controlled by transverse mixing, which is limited in this aquifer by low dispersion. Consumption to date of mineral oxides and SO4 is, respectively, <0.15% and 0.4% of the available aquifer capacity, and degradation using these oxidants is <10%. Fermentation is a significant process in contaminant turnover, accounting for 21% of degradation products present in the plume, and indicating that microbial respiration rates are slow in comparison with fermentation. Under present conditions, the potential for degradation in the plume is very low due to inhibitory effects of the contaminant matrix. Degradation products correspond to <22% mass loss over the life of the plume, providing a first-order plume scale half-life >140 years. The phenolic compounds are biodegradable under the range of redox conditions in the aquifer and the aquifer is not oxidant limited, but the plume is likely to be long-lived and to expand. Degradation is likely to increase only after contaminant concentrations are reduced and aqueous oxidant inputs are increased by dispersion of the plume. The results imply that transport processes may exert a greater control on the natural attenuation of this plume than aquifer oxidant availability.     

A new concept: the use of neutrally-buoyant microemulsions for DNAPL remediation

Even in the absence of mobilization of dense nonaqueous phase liquid (DNAPL), the microemulsion that forms when the surfactant solubilizes a dense contaminant such as trichloroethylene will be more dense than water and tends to migrate downward. This paper addresses the issue of migration with a new concept: surfactant enhanced aquifer remediation at neutral buoyancy. Laboratory results of surfactant remediation in two-dimensional model aquifers show that downward migration of microemulsion containing solubilized dense contaminants can be reduced to an acceptable level, even in the absence of capillary barriers in the aquifer. One model experiment was designed to exhibit a small degree of vertical migration and full capture of the microemulsion at the extraction well. The second experiment was designed to demonstrate the effect of large buoyancy forces that lead to excessive downward migration of the microemulsion. Density measurements of aqueous solutions containing sodium dihexyl sulfosuccinate surfactant, isopropanol, trichloroethylene, and sodium chloride are presented. A companion paper presents the results of the flow and transport calculations needed for this approach to surfactant flooding.     

New methodology to investigate potential contaminant mass fluxes at the stream-aquifer interface by combining integral pumping tests and streambed temperatures

The spatial pattern and magnitude of mass fluxes at the stream-aquifer interface have important implications for the fate and transport of contaminants in river basins. Integral pumping tests were performed to quantify average concentrations of chlorinated benzenes in an unconfined aquifer partially penetrated by a stream. Four pumping wells were operated simultaneously for a time period of 5 days and sampled for contaminant concentrations. Streambed temperatures were mapped at multiple depths along a 60m long stream reach to identify the spatial patterns of groundwater discharge and to quantify water fluxes at the stream-aquifer interface. The combined interpretation of the results showed average potential contaminant mass fluxes from the aquifer to the stream of 272microgm(-2)d(-1) MCB and 71microgm(-2)d(-1) DCB, respectively. This methodology combines a large-scale assessment of aquifer contamination with a high-resolution survey of groundwater discharge zones to estimate contaminant mass fluxes between aquifer and stream.     

Soil mobility of sewage sludge-derived dissolved organic matter, copper, nickel and zinc

A soil (sandy loam) column leaching study aimed to determine the extent of mobility and co-mobility of Cu, Ni, Zn and dissolved organic matter (DOM) released from a surface-application (equivalent to 50 t ds ha(-1) of anaerobically-digested sewage sludge. Leaching of DOM through the soil column was found to be almost un-retarded. Decidedly similar behaviour was exhibited by Ni suggesting that it migrated as organic complexes. Whilst Cu was also found to be leached, significant retardation was evident. However, the importance of DOM in promoting the mobility of both Cu and Ni was evidenced by their lack of mobility when added to the soil column as inorganic forms. The presence of DOM did not prevent Zn from becoming completely adsorbed by the soil solid phase. In relation to WHO drinking water guidelines, only Ni concentrations showed potential environmental significance, due to the relatively poor retention of Ni by the sludge solid phase.     

Humic acid enhanced remediation of an emplaced diesel source in groundwater. 1. Laboratory-based pilot scale test

The enhanced solubility of petroleum-derived compounds in humic acid solutions is the basis for a new groundwater remediation technology. In this unique pilot-scale test, a stationary contaminant source consisting of diesel fuel was placed below the water table in a model sand aquifer (1.2 x 5.5 x 1.8-m deep) and flushed with water at a flow rate of 2 cm/h over 5 years. At 51 days, laboratory grade humic acid was added to the water and maintained at a level of approximately 0.8 g/l. The addition of humic acid had only a small impact on the aqueous transport of the BTEX components, which were rapidly dissolved from the diesel, but had a large effect on the flushing of PAHs, including methylated naphthalenes (MNs). Binding to aqueous humic acid enhanced the solubilization of MNs two- to tenfold. During aqueous transport, biodegradation of the BTEX and PAHs occurred, limiting the lateral and longitudinal extent of the diesel contaminant plume in the model aquifer. It appears that through enhanced solubilization, the overall biodegradation rate of the MNs was increased. As the various MNs were depleted from the diesel source, the MN plume shrank and then disappeared.     

Assessing the natural attenuation of organic contaminants in aquifers using plume-scale electron and carbon balances: model development with analysis of uncertainty and parameter sensitivity

A quantitative methodology is described for the field-scale performance assessment of natural attenuation using plume-scale electron and carbon balances. This provides a practical framework for the calculation of global mass balances for contaminant plumes, using mass inputs from the plume source, background groundwater and plume residuals in a simplified box model. Biodegradation processes and reactions included in the analysis are identified from electron acceptors, electron donors and degradation products present in these inputs. Parameter values used in the model are obtained from data acquired during typical site investigation and groundwater monitoring studies for natural attenuation schemes. The approach is evaluated for a UK Permo-Triassic Sandstone aquifer contaminated with a plume of phenolic compounds. Uncertainty in the model predictions and sensitivity to parameter values was assessed by probabilistic modelling using Monte Carlo methods. Sensitivity analyses were compared for different input parameter probability distributions and a base case using fixed parameter values, using an identical conceptual model and data set. Results show that consumption of oxidants by biodegradation is approximately balanced by the production of CH4 and total dissolved inorganic carbon (TDIC) which is conserved in the plume. Under this condition, either the plume electron or carbon balance can be used to determine contaminant mass loss, which is equivalent to only 4% of the estimated source term. This corresponds to a first order, plume-averaged, half-life of > 800 years. The electron balance is particularly sensitive to uncertainty in the source term and dispersive inputs. Reliable historical information on contaminant spillages and detailed site investigation are necessary to accurately characterise the source term. The dispersive influx is sensitive to variability in the plume mixing zone width. Consumption of aqueous oxidants greatly exceeds that of mineral oxidants in the plume, but electron acceptor supply is insufficient to meet the electron donor demand and the plume will grow. The aquifer potential for degradation of these contaminants is limited by high contaminant concentrations and the supply of bioavailable electron acceptors. Natural attenuation will increase only after increased transport and dilution.     

Simulating the fate and transport of TCE from groundwater to indoor air

This work provides an exploratory analysis on the relative importance of various factors controlling the fate and transport of volatile organic contaminants (in this case, TCE) from a DNAPL source zone located below the water table and into the indoor air. The analysis is conducted using the multi-phase compositional model CompFlow Bio, with the base scenario problem geometry reminiscent of a field experiment conducted by Rivett [Rivett, M.O., (1995), Soil–gas signatures from volatile chlorinated solvents: Borden field experiments. Groundwater, 33(1), 84–98.] at the Borden aquifer where groundwater was observed to transport a contaminant plume a substantial distance without vertical mass transport of the contaminant across the capillary fringe and into the vadose zone. Results for the base scenario model indicate that the structure of the permeability field was largely responsible for deflecting the groundwater plume upward towards the capillary fringe, permitting aqueous phase diffusion to transport the TCE into the vadose zone. Alternative permeability realizations, generated as part of a Monte Carlo simulation process, at times deflected the groundwater plume downwards causing the extended thickness of the saturated zone to insulate the vadose zone from exposure to the TCE by upward diffusive transport. Comparison of attenuation coefficients calculated using the CompFlow Bio and Johnson and Ettinger [Johnson, P.C. and Ettinger, R.A., (1991), Heuristic model for predicting the intrusion rate of contaminant vapors into buildings. Environmental Science and Technology, 25, 1445–1452.] heuristic model exhibited fortuitous agreement for the base scenario problem geometry, with this agreement diverging for the alternative permeability realizations as well as when parameters such as the foundation slab fracture aperture, the indoor air pressure drop, the capillary fringe thickness, and the infiltration rate were varied over typical ranges.     

Colloid facilitated transport of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) to the groundwater at Ma Da Area, Vietnam

PCDD/Fs are hydrophobic organic substances and strongly sorbing to soil particles. Once adsorbed to soil particles they are believed to be virtually immobile. However, research in the last decades confirmed that strong sorbing contaminants may reach the groundwater via colloid-facilitated transport. This pathway has not been investigated before in Vietnam. Ma Da area, 100 km north of Ho Chi Minh City, was repeatedly sprayed during the Vietnam War (1962–1971) with herbicides like Agent Orange containing, beside others, the teratogenic contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). 11 surface soil samples and 12 water samples were collected in Ma Da area for analysis of PCDD/Fs in solids. Soil TCDD concentrations ranged from 1–41 ppt with a mean of 8.8 ppt and a mean I-TEQ of 9.7 ppt. Two surface water samples showed colloid bound TCDD (7 and 19 ppt). Groundwater samples showed elevated colloid bound PCDD concentrations (mean 770 ng/kg), mainly octachlorodibenzo-p-dioxin. Groundwater colloids separated by filtration did not show any TCDD. The results support that TCDD/Fs can be relocated from the top soil to the groundwater by colloidal pathway. They did not provide evidence that the dioxins bound to groundwater colloids are leftovers from the Second Indochinese War. However, this study reinforces that the colloidal transport pathway has to be included investigating the relocation of strong sorbing organic contaminants.     

The effect of local-scale physical heterogeneity and nonlinear,rate-limited sorption/desorption on contaminant transport in porous media

Nonideal transport of contaminants in porous media has often been observed in laboratory characterization studies. It has long been recognized that multiple processes associated with both physical and chemical factors can contribute to this nonideal transport behavior. To fully understand system behavior, it is important to determine the relative contributions of these multiple factors when conducting contaminant transport and fate studies. In this study, the relative contribution of physical-heterogeneity-related processes versus those of nonlinear, rate-limited sorption/desorption to the observed nonideal transport of trichloroethene in an undisturbed aquifer core was determined through a series of miscible-displacement experiments. The results of experiments conducted using the undisturbed core, collected from a Superfund site in Tucson, AZ, were compared to those obtained from experiments conducted using the same aquifer material packed homogeneously. The results indicate that both physical and chemical factors, specifically preferential flow and associated rate-limited diffusive mass-transfer and rate-limited sorption/desorption, respectively, contributed to the nonideal behavior observed for trichloroethene transport in the undisturbed core. A successful prediction of trichloroethene transport in the undisturbed core was made employing a mathematical model incorporating multiple sources of nonideal transport, using independently determined model parameters to account for the multiple factors contributing to the nonideal transport behavior. The simulation results indicate that local-scale physical heterogeneity controlled the nonideal transport behavior of trichloroethene in the undisturbed core, and that nonlinear, rate-limited sorption/desorption were of secondary importance.