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.     

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.     

Contaminant transport in groundwater in the presence of colloids and bacteria: Model development and verification

Colloids and bacteria (microorganisms) naturally exist in groundwater aquifers and can significantly impact contaminant migration rates. A conceptual model is first developed to account for the different physiochemical and biological processes, reaction kinetics, and different transport mechanisms of the combined system (contaminant–colloids–bacteria). All three constituents are assumed to be reactive with the reactions taking place between each constituent and the porous medium and also among the different constituents. A general linear kinetic reaction model is assumed for all reactive processes considered. The mathematical model is represented by fourteen coupled partial differential equations describing mass balance and reaction processes. Two of these equations describe colloid movement and reactions with the porous medium, four equations describe bacterial movement and reactions with colloids and the porous medium, and the remaining eight equations describe contaminant movement and its reactions with bacteria, colloids, and the porous medium. The mass balance equations are numerically solved for two-dimensional groundwater systems using a third-order, total variance-diminishing scheme (TVD) for the advection terms. Due to the complex coupling of the equations, they are solved iteratively each time step until a convergence criterion is met. The model is tested against experimental data and the results are favorable.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

A novel two-dimensional model for colloid transport in physically and geochemically heterogeneous porous media

A two-dimensional model for colloid transport in geochemically and physically heterogeneous porous media is presented. The model considers patchwise geochemical heterogeneity, which is suitable to describe the chemical variability of many surficial aquifers with ferric oxyhydroxide-coated porous matrix, as well as spatial variability of hydraulic conductivity, which results in heterogeneous flow field. The model is comprised of a transient fluid flow equation, a transient colloid transport equation, and an equation for the dynamics of colloid deposition and release. Numerical simulations were carried out with the model to investigate the colloid transport behavior in layered and randomly heterogeneous porous media. Results demonstrate that physical and geochemical heterogeneities markedly affect the colloid transport behavior. Layered physical or geochemical heterogeneity can result in distinct preferential flow paths of colloidal particles. Furthermore, the combined effect of layered physical and geochemical heterogeneity may result in enhanced or reduced preferential flow of colloids. Random distribution of physical heterogeneity (hydraulic conductivity) results in a random flow field and an irregularly distributed colloid concentration profile in the porous medium. Contrary to random physical heterogeneity, the effect of random patchwise geochemical heterogeneity on colloid transport behavior is not significant. It is mostly the mean value of geochemical heterogeneity rather than its distribution that governs the colloid transport behavior.     

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.     

Gel barrier formation in unsaturated porous media

The gel barrier formation by a gelling liquid (Colloidal Silica) injection in an unsaturated porous medium is investigated by developing a mathematical model and conducting numerical simulations. Gelation process is initiated by adding electrolytes such as NaCl, and the gel phase consisting of cross-linked colloidal silica particles grows as the gelation process proceeds. The mathematical model describing the transport and gelation of Colloidal Silica (CS) is based on coupled mass balance equations for the gel mixture (the sol phase plus the gel phase), gel phase (cross-linked colloidal silica particles plus water captured between cross-linked particles), and colloidal silica particles (discrete and cross-linked) and NaCl in the sol (suspension of discrete colloidal silica particles in water) and gel phases. The solutions in terms of volumetric fraction of the gel phase yield the gel mixture viscosity via the dependency on the volumetric fraction of gel phase. This dependency is determined from a kinetic gelation model with time-normalized viscosity curves. The proposed model is verified by comparing experimentally and numerically determined hydraulic conductivities of gel-treated soil columns at different CS injection volumes. The numerical experiments indicate that an impermeable gel layer is formed within the time period twice the gel-point in a one-dimensional flow system. At the same normalized time corresponding to twice the gel-point, the CS solutions with lower NaCl concentrations result in further migration and poor performance in plugging the pore space. The viscosity computation proposed in this study is compared with another method available in the literature. It is observed that the other method estimates the viscosity at the mixing zone higher than the one proposed by the authors. The proposed model can simulate realistic injection scenarios with various combinations of operating parameters such as NaCl concentration and NaCl mixing time, and thus providing guidelines in performing this technology on site.     

Effect of dissolved organic matter and bacteria on contaminant transport in riverbank filtration

A mathematical model for the transport of hydrophobic organic contaminants in an aquifer under simplistic riverbank filtration conditions is developed. The model considers a situation where contaminants are present together with dissolved organic matter (DOM) and bacteria. The aquifer is conceptualized as a four-phase system: two mobile colloidal phases, an aqueous phase, and a stationary solid phase. An equilibrium approach is used to describe the interactions of contaminants with DOM, bacteria, and solid matrix. The model is composed of bacterial transport equation and contaminant transport equation. Numerical simulations are performed to examine the contaminant transport behavior in the presence of DOM and bacteria. The simulation results illustrate that contaminant transport is enhanced markedly in the presence of DOM and bacteria, and the impact of DOM on contaminant mobility is greater than that of bacteria under examined conditions. Sensitivity analysis demonstrates that the model is sensitive to changes of three lumped parameters: K+1 (total affinity of stationary solid phase to contaminants), K+2 (total affinity of DOM to contaminants), and K+3 (total affinity of bacteria to contaminants). In a situation where contaminants exist simultaneously with DOM and bacteria, contaminant transport is mainly affected by a ratio of K+1/K+2/K+3, which can vary with changes of equilibrium distribution coefficient of contaminants and/or colloidal concentrations. In riverbank filtration, the influence of DOM and bacteria on the transport behavior of contaminants should be accounted to accurately predict the contaminant mobility.     

Preliminary 3-D site-scale studies of radioactive colloid transport in the unsaturated zone at Yucca Mountain,Nevada

The U.S. Department of Energy (DOE) is actively investigating the technical feasibility of permanent disposal of high-level nuclear waste in a repository to be situated in the unsaturated zone (UZ) at Yucca Mountain (YM), Nevada. In this study we investigate, by means of numerical simulation, the transport of radioactive colloids under ambient conditions from the potential repository horizon to the water table. The site hydrology and the effects of the spatial distribution of hydraulic and transport properties in the Yucca Mountain subsurface are considered. The study of migration and retardation of colloids accounts for the complex processes in the unsaturated zone of Yucca Mountain, and includes advection, diffusion, hydrodynamic dispersion, kinetic colloid filtration, colloid straining, and radioactive decay. The results of the study indicate that the most important factors affecting colloid transport are the subsurface geology and site hydrology, i.e., the presence of faults (they dominate and control transport), fractures (the main migration pathways), and the relative distribution of zeolitic and vitric tuffs. The transport of colloids is strongly influenced by their size (as it affects diffusion into the matrix, straining at hydrogeologic unit interfaces, and transport velocity) and by the parameters of the kinetic-filtration model used for the simulations. Arrival times at the water table decrease with an increasing colloid size because of smaller diffusion, increased straining, and higher transport velocities. The importance of diffusion as a retardation mechanism increases with a decreasing colloid size, but appears to be minimal in large colloids.     

Colloid transport in unsaturated porous media: the role of water content and ionic strength on particle straining

Packed column and mathematical modeling studies were conducted to explore the influence of water saturation, pore-water ionic strength, and grain size on the transport of latex microspheres (1.1 microm) in porous media. Experiments were carried out under chemically unfavorable conditions for colloid attachment to both solid-water interfaces (SWI) and air-water interfaces (AWI) using negatively charged and hydrophilic colloids and modifying the solution chemistry with a bicarbonate buffer to pH 10. Interaction energy calculations and complementary batch experiments were conducted and demonstrated that partitioning of colloids to the SWI and AWI was insignificant across the range of the ionic strengths considered. The breakthrough curve and final deposition profile were measured in each experiment indicating colloid retention was highly dependent on the suspension ionic strength, water content, and sand grain size. In contrast to conventional filtration theory, most colloids were found deposited close to the column inlet, and hyper-exponential deposition profiles were observed. A mathematical model, accounting for time- and depth-dependent straining, produced a reasonably good fit for both the breakthrough curves and final deposition profiles. Experimental and modeling results suggest that straining--the retention of colloids in low velocity regions of porous media such as grain junctions--was the primary mechanism of colloid retention under both saturated and unsaturated conditions. The extent of stagnant regions of flow within the pore structure is enhanced with decreasing water content, leading to a greater amount of retention. Ionic strength also contributes to straining, because the number of colloids that are held in the secondary energy minimum increases with ionic strength. These weakly associated colloids are prone to be translated to stagnation regions formed at grain-grain junctions, the solid-water-air triple point, and dead-end pores and then becoming trapped.     

Modeling radium and radon transport through soil and vegetation

A one-dimensional flow and transport model was developed to describe the movement of two fluid phases, gas and water, within a porous medium and the transport of 226Ra and 222Rn within and between these two phases. Included in this model is the vegetative uptake of water and aqueous 226Ra and 222Rn that can be extracted from the soil via the transpiration stream. The mathematical model is formulated through a set of phase balance equations and a set of species balance equations. Mass exchange, sink terms and the dependence of physical properties upon phase composition couple the two sets of equations. Numerical solution of each set, with iteration between the sets, is carried out leading to a set-iterative compositional model. The Petrov-Galerkin finite element approach is used to allow for upstream weighting if required for a given simulation. Mass lumping improves solution convergence and stability behavior. The resulting numerical model was applied to four problems and was found to produce accurate, mass conservative solutions when compared to published experimental and numerical results and theoretical column experiments. Preliminary results suggest that the model can be used as an investigative tool to determine the feasibility of phytoremediating radium and radon-contaminated soil.     

Removal of NAPLs from the unsaturated zone using steam: prevention of downward migration by injecting mixtures of steam and air

Steam injection for remediation of porous media contaminated by nonaqueous phase liquids has been shown to be a potentially efficient technology. There is, however, concern that the technique may lead to downward migration of separate phase contaminant. In this work, a modification of the steam injection technology is presented, where a mixture of steam and air was injected. In two-dimensional experiments with unsaturated porous medium contaminated with nonaqueous phase liquids, it was demonstrated how injection of pure steam lead to severe downward migration. Similar experiments, where steam and air were injected simultaneously, resulted in practically no downward migration and still rapid cleanup was achieved. The processes responsible for the prevention of downward migration when injecting steam-air mixtures were analyzed using a nonisothermal multiphase flow and transport model. Hereby, three mechanisms were identified and it was demonstrated how the effectiveness of these mechanisms depended on the air-to-steam mixing ratio.     

Influence of rate-limited sorption on the cleanup of layered soils by vapor extraction

The conditions under which rate-limited sorption is important for cleanup of layered soils by vapor extraction are investigated. The investigation includes two steps: (a) the cleanup time is estimated for a number of scenario cases by means of a numerical model and (b) the numerical results are approximated using analytical solutions derived for simplified models. In this way, equations are derived, which give insight into the influence of different parameters characterizing the properties of the soil, the geometry of the formation, the mass transfer mechanisms in it, and the distribution of the contaminant mass in the different phases (gas phase, water phase and solid phase). The numerical model used is based on the advection-dispersion differential equations for Darcian isothermal airflow, local equilibrium contaminant mass transfer between gas phase and soil water and first-order kinetics for mass transfer between soil water and solid phase. The numerical results are approximated combining an analytical solution to estimate cleanup time in layered formations for local equilibrium sorption, which has been presented in a previous work (J. Contam. Hydrol., 36 (1999) 105). with an analytical solution based on the well-mixed reservoir model under consideration of rate-limited sorption. The analytical approximation of the cleanup time is in reasonable agreement with the numerical results and allows its estimation with small computational effort.     

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.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Transport of volatile compounds in porous media in the presence of a trapped gas phase

The presence of an immobile gaseous phase in an otherwise-saturated porous medium affects the transport of volatile compounds. The linear theory of partitioning tracers suggests that a volatile tracer introduced into such a system should be retarded with a constant retardation factor. Using high concentrations, however, the saturation of the gaseous phase will change as an effect of the tracer test itself. Competitive gas transfer among all volatile compounds and the change of saturation may lead to tracer concentrations that are temporarily higher than those injected. We analyze the system in the framework of the coherence theory by Helfferich [Soc. Pet. Eng. J. 21 (1) (1981) 51]. The governing equations are formulated as functions of total concentration, i.e., the mass of solutes in all phases per pore volume. Neglecting dispersion and mass-transfer kinetics, we derive the characteristic form of the resulting system of hyperbolic equations. In a system with N volatile compounds, a variation of the concentrations splits up into N waves, each traveling with its own characteristic velocity. If the presence of a gaseous phase is sustained, one wave will be a standing one. We perform numerical model calculations for tracers with various Henry's law coefficients and show that the results agree with the semi-analytical solution obtained by coherence theory.     

Anomalous transport of colloids and solutes in a shear zone

Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters.     

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

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

Modeling in situ ozonation for the remediation of nonvolatile PAH-contaminated unsaturated soils

Mathematical models were developed to investigate the characteristics of gaseous ozone transport under various soil conditions and the feasibility of in situ ozone venting for the remediation of unsaturated soils contaminated with phenanthrene. On the basis of assumptions for the mass transfer and reactions of ozone, three approaches were considered: equilibrium, kinetic, and lump models. Water-saturation-dependent reactions of gaseous ozone with soil organic matter (SOM) and phenanthrene were employed. The models were solved numerically by using the finite-difference method, and the model parameters were determined by using the experimental data of Hsu [The use of gaseous ozone to remediate the organic contaminants in the unsaturated soils, PhD Thesis, Michigan State Univ., East Lansing, MI, 1995]. The transport of gas-phase ozone is significantly retarded by ozone consumption due to reactions with SOM and phenanthrene, in addition to dissolution. An operation time of 156 h was required to completely remove phenanthrene in a 5-m natural soil column. In actual situations, however, the operation time is likely to be longer than the ideal time because of unknown factors including heterogeneity of the porous medium and the distribution of SOM and contaminant. The ozone transport front length was found to be very limited (< 1 m). The sensitivity analysis indicated that SOM is the single most important factor affecting in situ ozonation for the remediation of unsaturated soil contaminated with phenanthrene. Models were found to be insensitive to the reaction mechanisms of phenathrene with either gas-phase ozone or dissolved ozone. More study is required to quantify the effect of OH* formation on the removal of contaminant and on ozone transport in the subsurface.