基于GIS城市机动车尾气扩散迁移研究

中国城市大气污染日趋严重,汽车尾气污染所占比重越来越大。交通、道路、管制及气候条件等是影响城市机动车尾气排放的主要因素。本文根据机动车尾气研究的实际,研究了专题MGIS数据库结构, 并以广州为例,建立了基于GIS的城市机动车尾气扩散迁移预报系统,在此基础上实现了机动车尾气污染状况的可视化。     

Field study of TCE diffusion profiles below DNAPL to assess aquitard integrity

An area where a free-product accumulation of trichloroethylene (TCE) dense non-aqueous phase liquid (DNAPL) occurs at the bottom of a 10-m-thick surficial sand aquifer was studied to determine the integrity of the underlying, 20-m-thick, clayey silt aquitard formed of glaciolacustrine sediment. TCE concentration-versus-depth profiles determined from aquitard cores collected at five locations indicated penetration of detectable TCE 2.5 to 3.0 m into the aquitard. Two of the profiles show persistent DNAPL at the aquitard interface, while two others indicate that DNAPL, present initially, was completely dissolved away producing concentration declines at the aquitard interface. The fifth profile suggests shallow DNAPL penetration (<0.5 m) into the aquitard, however, this penetration, which was likely caused by cross-contamination during core collection or cone penetrometry (CPT) of the aquitard interface, did not increase the maximum depth of TCE penetration. Combining the field profiles with one-dimensional model simulations, downward migration of the aqueous TCE front, defined as the EPA MCL of 5 microg/l, which was below the analytical detection limit, was projected to a distance between 4 and 5 m below the top of the aquitard. Using a single set of estimated aquitard parameter values, simulations of aqueous TCE migration into the aquitard provided a good fit to four of the field profiles with a migration time of 35 to 45 years, consistent with the history of TCE use at the site. These simulations indicate aqueous TCE migration is diffusion-dominated with only small advective influence by the downward groundwater velocity of 2 to 3 cm/year or less in the aquitard due to pumping of the underlying aquifer to supply water to the facility in the past 50 years. The applicability of the parameter values was confirmed by in situ diffusion experiments of 1-year duration, in which stainless steel cylinders containing DNAPL were inserted into the aquitard. The diffusion-dominated nature of the profiles indicates that the aquitard provides long-term protection of the underlying aquifer from contamination from this DNAPL zone. Simulations of long-term migration of the TCE solute front indicate breakthrough to the lower aquifer at 1200 years for the no advection scenario and at 500 years if the strong downward hydraulic gradient persists. However, even after breakthrough, the mass flux through the aquitard to the underlying aquifer remains relatively low, and when considered in terms of potential impacts to pumping wells, concentrations are not expected to increase significantly above present-day MCLs. The use of contaminant profiles of different time and distance scales, in addition to hydraulic data, dramatically improves the ability to assess aquitard integrity, and provides improved transport parameter values for estimating contaminant arrival times and fluxes. The apparent lack of deep preferential pathways for TCE migration, such as open fractures, is probably due to the softness of the silty aquitard deposit and minimal physical or chemical weathering of the aquitard provides long-term protection of the underlying aquifer from contamination from this DNAPL zone. Simulations of long-term migration of the TCE solute front indicate breakthrough to the lower aquifer at 1200 years for the no advection scenario and at 500 years if the strong downward hydraulic gradient persists. However, even after the breakthrough, the mass flux through the aquitard to the underlying aquifer remains relatively low, and when considered in terms of potential impacts to pumping wells , concentrations are not expected to increase significantly above present-day MCLs. The use of contaminant profiles of different time and distance scales, in addition to hydraulic data, dramatically improves the ability to assess aquitard integrity, and provides improved transport parameter values for estimating contaminant arrival times and fluxes. The apparent lack of deep preferential pathways for TCE migration, such as open fractures, is probably due to the softness of the silty aquitard deposit and minimal physical or chemical weathering of the aquitard.     

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

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

A method for estimating distributions of mass transfer rate coefficients with application to purging and batch experiments

Mass transfer between aquifer material and groundwater is often modeled as first-order rate-limited sorption or diffusive exchange between mobile zones and immobile zones with idealized geometries. Recent improvements in experimental techniques and advances in our understanding of pore-scale heterogeneity demonstrate that two (or even a few) rate coefficients are insufficient in many cases. Here, we investigate a piece-wise linear model for a continuous distribution of rate coefficients, that has several advantages over previously used ‘statistical' distribution models (with functional form from gamma or lognormal PDF's): (1) distributions of arbitrary, even bimodal, shapes can be represented; (2) linear estimation methods can be applied to determine the distribution from experimental data; (3) the uncertainty in the distribution can be determined for each of its sections; and (4) the relationship between the time scales of available data and those of estimatable mass transfer processes can be investigated. A statistical model refinement algorithm is presented that reduces the number of parameters (sections of the piece-wise linear model) to the admissible minimum. We show that purging experiments allow estimation of a wider zone of the rate distribution than do batch experiments, and hence will provide predictions that are accurate over a wider range of time scales. Finally, in an application to TCE gas-purging desorption data, the piece-wise linear rate-distribution model has a higher probability of being adequate than those using a gamma or lognormal distribution or a single rate coefficient.     

Analytical solution of two-dimensional solute transport in an aquifer–aquitard system

This study deals with two-dimensional solute transport in an aquifer–aquitard system by maintaining rigorous mass conservation at the aquifer–aquitard interface. Advection, longitudinal dispersion, and transverse vertical dispersion are considered in the aquifer. Vertical advection and diffusion are considered in the aquitards. The first-type and the third-type boundary conditions are considered in the aquifer. This study differs from the commonly used averaged approximation (AA) method that treats the mass flux between the aquifer and aquitard as an averaged volumetric source/sink term in the governing equation of transport in the aquifer. Analytical solutions of concentrations in the aquitards and aquifer and mass transported between the aquifer and upper or lower aquitard are obtained in the Laplace domain, and are subsequently inverted numerically to yield results in the real time domain (the Zhan method). The breakthrough curves (BTCs) and distribution profiles in the aquifer obtained in this study are drastically different from those obtained using the AA method. Comparison of the numerical simulation using the model MT3DMS and the Zhan method indicates that the numerical result differs from that of the Zhan method for an asymmetric case when aquitard advections are at the same direction. The AA method overestimates the mass transported into the upper aquitard when an upward advection exists in the upper aquitard. The mass transported between the aquifer and the aquitard is sensitive to the aquitard Peclet number, but less sensitive to the aquitard diffusion coefficient.     

Contributions of advective and diffusive oxygen transport through multilayer composite caps over mine waste

The relative contributions of four mechanisms of oxygen transport in multilayer composite (MLC) caps placed over oxygen-consuming mine waste were evaluated using numerical and analytical methods. MLC caps are defined here as caps consisting of earthen and geosynthetic (polymeric) components where a composite barrier layer consisting of a geomembrane (1-2 mm thick polymeric sheet) overlying a clay layer is the primary barrier to transport. The transport mechanisms that were considered are gas-phase advective transport, gas-phase diffusive transport, liquid-phase advective transport via infiltrating precipitation and liquid-phase diffusive transport. A numerical model was developed to simulate gas-phase advective-diffusive transport of oxygen through a multilayer cap containing seven layers. This model was also used to simulate oxygen diffusion in the liquid phase. An approximate analytical method was used to compute the advective flux of oxygen in the liquid phase. The numerical model was verified for limiting cases using an analytical solution. Comparisons were also made between model predictions and field data for earthen caps reported by others. Results of the analysis show that the dominant mechanism for oxygen transport through MLC caps is gas-phase diffusion. For the cases that were considered, the gas-phase diffusive flux typically comprises at least 99% of the total oxygen flux. Thus, designers of MLC caps should focus on design elements and features that will limit diffusion of gas-phase oxygen.     

Field-scale effective matrix diffusion coefficient for fractured rock: results from literature survey

Matrix diffusion is an important mechanism for solute transport in fractured rock. We recently conducted a literature survey on the effective matrix diffusion coefficient, D_m^e, a key parameter for describing matrix diffusion processes at the field scale. Forty field tracer tests at 15 fractured geologic sites were surveyed and selected for the study, based on data availability and quality. Field-scale D_m^e values were calculated, either directly using data reported in the literature, or by reanalyzing the corresponding field tracer tests. The reanalysis was conducted for the selected tracer tests using analytic or semi-analytic solutions for tracer transport in linear, radial, or interwell flow fields. Surveyed data show that the scale factor of the effective matrix diffusion coefficient (defined as the ratio of D_m^e to the lab-scale matrix diffusion coefficient, D_m, of the same tracer) is generally larger than one, indicating that the effective matrix diffusion coefficient in the field is comparatively larger than the matrix diffusion coefficient at the rock-core scale. This larger value can be attributed to the many mass-transfer processes at different scales in naturally heterogeneous, fractured rock systems.Furthermore, we observed a moderate, on average trend toward systematic increase in the scale factor with observation scale. This trend suggests that the effective matrix diffusion coefficient is likely to be statistically scale-dependent. The scale-factor value ranges from 0.5 to 884 for observation scales from 5 to 2000 m. At a given scale, the scale factor varies by two orders of magnitude, reflecting the influence of differing degrees of fractured rock heterogeneity at different geologic sites. In addition, the surveyed data indicate that field-scale longitudinal dispersivity generally increases with observation scale, which is consistent with previous studies. The scale-dependent field-scale matrix diffusion coefficient (and dispersivity) may have significant implications for assessing long-term, large-scale radionuclide and contaminant transport events in fractured rock, both for nuclear waste disposal and contaminant remediation.     

The diffusion and sorption of volatile organic compounds through kaolinitic clayey soils

Laboratory experiments to estimate the effective molecular diffusion coefficient (D(e)) and sorption coefficient (K(d)) for volatile organic compounds through natural clayey soils were conducted using diffusion testing apparatus. The compounds tested were methyl ethyl ketone (MEK), toluene and trichloroethylene (TCE). The D(e) and K(d) values were determined by a curve fitting procedure. The compound losses, and the effects of porous disks used in the apparatus were significant. The transport of MEK was faster than that of TCE and toluene because of the lower sorption to the soils. The D(e) values of all the compounds were of the order of 10(-10) m(2)/s and smaller than the diffusion coefficient in pure aqueous solution at infinite dilution (D(0)), due to the tortuosity of the samples. The effects of the sample thickness on the parameter determination were not significant. Comparison to the K(d) values estimated from batch sorption tests and from organic carbon content (f(oc))-based predictions showed that the diffusion test results were intermediate between those from the other two methods. The diffusion tests use compacted soil samples and should be more relevant to in situ conditions, but the reliability of the tests is affected by large compound losses that cause uncertainties in their interpretation. It is recommended that more than one method be used to assess K(d) values.     

Analysis of the long-term behavior of solute transport with nonlinear equilibrium sorption using breakthrough curves and temporal moments

We analyzed the long-term behavior of breakthrough curves (BTCs) and temporal moments of a solute subjected to Freundlich equilibrium sorption (s = kc(n)). For one-dimensional transport in a homogeneous porous medium, we derived a power-law relation between travel time, tau, and solute displacement, chi, with the exponent being equal to the Freundlich n exponent. The mean solute velocity, derived from the first time moment, was found to change as tau(n-1). For n values larger than 0.66, the second time moment could be related to c chi(2/n), where c is a constant. An approach based on the use of a critical concentration was developed to estimate the presence of the asymptotic regime in the tail of the BTC. This approach was tested successfully using numerical case studies. One-dimensional numerical simulations with varying values of k, n and initial mass were run to verify the closed form analytical expressions for the large time behavior of temporal moments and the tailing part of breakthrough curves. Good agreement between the slope of the tailing part of log-log transformed BTCs and the predicted slope using asymptotic theory was found. Asymptotic theory in general underestimated the magnitude of the concentration in the tail. The quality of the estimated concentrations in the tail improved for small values of the dispersivity. Experimental BTCs of uranin and benazolin were analyzed in combination with sorption/desorption batch experiments using asymptotic theory. A good agreement between the value of n parameter derived from desorption experiment with benazolin and the value of the n parameter derived from the tail of the BTC was found.     

Quantifying the physical and chemical mass transfer processes for the fate and transport of Co(II)EDTA in a partially-weathered limestone-shale saprolite

The objective of the research is to quantify the relative contributions of physical and chemical mass transfer to the movement of Co(II/III)EDTA (chelates of Cobalt and Ethylene Diamine Tetraacetic Acid or EDTA) through a limestone-shale saprolite soil. Saprolite is a collective term referring to partially-weathered bedrock. It exists extensively in the subsurface. Because the parent bedding structures are maintained during the weathering process, saprolite soils are characterized by intensive fractures and secondary deposits of minerals such as Al-, Fe- and Mn-oxides on the fracture surfaces. Movement of reactive species through the soils may be influenced by diffusion into the rock matrix, a physical mass transfer (PMT) process, and interfacial chemical reactions, a chemical mass transfer (CMT) process. The PMT and CMT processes are phenomenologically similar but mechanistically different. In this research, previous laboratory observations from a Br and Co(II)EDTA tracer injection into an undisturbed saprolite soil column were used. Mechanistic reactive transport models were formulated to quantify the PMT and CMT processes. The PMT process was independently characterized by using the non-reactive tracer Br. Model parameters thus obtained were subsequently used as constraints to quantify the CMT processes involving Co(II)EDTA and its oxidation product Co(III)EDTA. Our calculations indicated that the PMT rates of the less reactive Co(III)EDTA were comparable with their theoretical CMT rates. In contrast, for the more reactive species Co(II)EDTA, CMT rates are higher than PMT rates. Evaluations of alternative CMT process models further confirmed one of our hypotheses on the basis of previous experimental understandings. The hypothesis suggested that competition from Fe-oxide for Co(II)EDTA may account for the majority of the decrease of Co(III)EDTA effluent concentrations that resulted in the separation of total Co and Co(III)EDTA breakthrough curves. Because Co(III)EDTA is more mobile than Co(II)EDTA in the subsurface, the results of this research suggest independent quantifications of CoEDTA PMT and CMT processes if laboratory results are to be interpreted correctly and scaled up for field and predictive uses.     

Explicit and implicit coupling during solute transport through clay membrane barriers

Simulations of salt (KCl) flux through a 1-m-thick clay membrane barrier (CMB) based on coupled solute transport theory are compared to simulated fluxes based on traditional advective-dispersive transport theory. The simulations are based on measured values for the effective salt-diffusion coefficient (Ds*) and chemico-osmotic efficiency coefficient (omega) for a bentonite-based barrier material subjected to KCl solutions. The results indicate that the exit salt flux is reduced due to both explicit coupling (hyperfiltration and chemico-osmotic counter-advection) and an implicit coupling effect resulting from the decrease in Ds* due to a decrease in the apparent tortuosity factor, tau a, with an increase in omega. Implicit coupling is shown to be more significant than explicit coupling for reducing and retarding salt flux through a CMB under diffusion-dominated conditions. Failure to account for the implicit coupling effect may result in unrealistic results, such as the existence of salt flux through a perfect (ideal) clay membrane (i.e., omega=1).     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

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.     

A model of oscillatory transport in granular soils, with application to barometric pumping and earth tides

A simple algebraic model is proposed to estimate the transport of a volatile or soluble chemical caused by oscillatory flow of fluid in a porous medium. The model is applied to the barometric pumping of vapors in the vadose zone, and to the transport of dissolved species by earth tides in an aquifer. In the model, the fluid moves sinusoidally with time in the porosity of the soil. The chemical concentration in the mobile fluid is considered to equilibrate with the concentration in the surrounding matrix according to a characteristic time governed by diffusion, sorption, or other rate processes. The model provides a closed form solution, to which barometric pressure data are applied in an example of pore gas motion in the vadose zone. The model predicts that the additional diffusivity due barometric pumping in an unfractured vadose zone would be comparable to the diffusivity in stagnant pore gas if the equilibration time is 1 day or longer. Water motion due to the M2 lunar tide is examined as an example of oscillatory transport in an aquifer. It is shown that the tidal motion of the water in an aquifer might significantly increase the vertical diffusivity of dissolved species when compared to diffusion in an absolutely stagnant aquifer, but the hydrodynamic dispersivity due to tidal motion or gravitational flow would probably exceed the diffusivity due to oscillatory advection.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

The effect of a biofilm on solute diffusion in fractured porous media

At sites in fractured rock where contamination has been exposed to the rock matrix for extended periods of time, the amount of contaminant mass residing in the matrix can be considerable. Even though it may be possible to diminish concentrations by the advection of clean water through the fracture features, back diffusion from mass held in the matrix will lead to a continuing source of contamination. In such an event, the development of a biofilm (a thin film of microbial mass) on the wall of the fractures may act to limit or prevent the back diffusion process. The objective of this preliminary study is to explore the influence imparted by the presence of a biofilm on the process of matrix diffusion. The investigation was conducted using radial diffusion cells constructed from rock core in which biofilm growth was stimulated in a central reservoir. Once biofilms were developed, forward diffusion experiments were conducted in which a conservative solute migrated from the central reservoir into the intact rock sample. Diffusion experiments were performed in a total of 11 diffusion cell pairs where biofilm growth was stimulated in one member of the pair and inhibited in the other. The effect of the presence of a biofilm on tracer diffusion was determined by comparison of the diffusion curves produced by each cell pair. A semi-analytical model that accounts for the presence of a biofilm was used to investigate the effect of the biofilm on mass transfer due to changes in the effective porosity, effective diffusion coefficient, and the depth of penetration of the biofilm into the intact rock. The results show that the biofilm acted to plug the rock matrix, rather than forming a discrete layer on the reservoir surface. The reduction in effective porosity due to the biofilm ranged from 6% to 52% with the majority of the samples in the 30% to 50% range. Based on the present results, with more efficient biofilm stimulation, it is reasonable to assume that a more complete plugging of the microcrack porosity might be possible, leaving a much thicker and efficient barrier than could be achieved via a surface biofilm.     

Chromium transport, oxidation, and adsorption in manganese-coated sand

We examine how the processes of advection, dispersion, oxidation-reduction, and adsorption combine to affect the transport of chromium through columns packed with pyrolusite (beta-MnO2)-coated sand. We find that beta-MnO2 effectively oxidizes Cr(III) to Cr(VI) and that the extent of oxidation is sensitive to changes in pH, pore water velocity, and influent concentrations of Cr(III). Cr(III) oxidation rates, although initially high, decline well before the supply of beta-MnO2 is depleted, suggesting that a reaction product inhibits the conversion of Cr(III) to Cr(VI). Rate-limited reactions govern the weak adsorption of each chromium species, with Cr(III) adsorption varying directly with pH and Cr(VI) adsorption varying inversely with pH. The breakthrough data on chromium transport can be matched closely by calculations of a simple model that accounts for (1) advective-dispersive transport of Cr(III), Cr(VI), and dissolved oxygen, (2) first-order kinetics adsorption of the reduced and oxidized chromium species, and (3) nonlinear rate-limited oxidation of Cr(III) to Cr(VI). Our work supplements the limited database on the transport of redox-sensitive metals in porous media and provides a means for quantifying the coupled processes that contribute to this transport.     

A laboratory technique for investigation of diffusion and transformation of volatile organic compounds in low permeability media

A laboratory diffusion cell technique that permits spatial and temporal estimates of porewater concentrations over short intervals suitable for estimation of effective diffusion coefficients (De) and degradation rate constants (k) of volatile organic compounds (VOCs) in saturated low permeability media is presented. The diffusion cell is a sealed cylinder containing vapour reservoirs for sampling, including a vapour reservoir source and an array of vapour-filled "mini-boreholes", which are maintained water- and sediment-free by slightly negative porewater pressures. The vapour reservoirs were sampled by Solid Phase Micro-Extraction (SPME), resulting in minimal disturbance to the experimental system. Porewater concentrations are estimated from the measured vapour concentrations. Experiments were conducted using a non-reactive medium and five VOCs with a range in partitioning properties. Calibration experiments showed that equilibrium partition coefficients could be used for calculating concentrations in the vapour reservoir source from concentrations in the SPME coating after a 1-min microextraction and that the reservoir concentration was insignificantly affected by sampling. However, equilibrium was not reached during the one-min extraction of the boreholes; the microextraction reduced the borehole vapour concentrations, leading to diffusion of VOCs from porewater into the vapour-filled borehole. Thus, empirical partitioning coefficients were required for the determination of porewater VOC concentrations. The experimental data and numerical modelling indicate masses extracted by SPME extraction are relatively small, with minimal perturbation on processes studied in diffusion experiments. This technique shows promise for laboratory investigation of diffusion and transformation processes in low permeability media.     

Transport modeling of copper and cadmium with linear and nonlinear retardation factors

The predictive accuracy of using the one-dimensional advection–dispersion equation to evaluate the fate and transport of solute in a soil column is usually dependent on the proper determination of chemical retardation factors. Typically, the distribution coefficient (K_d) obtained by fitting the linear sorption isotherm has been extensively used to consider general geochemical reactions on solute transport in a low-concentration range. However, the linear distribution coefficient cannot be adequately utilized to describe the solute fate at a higher concentration level. This study employed the nonlinear equilibrium-controlled sorption parameters to determine the retardation factor used in column leaching experiments. Copper and cadmium transportation in a lateritic silty-clay soil column was examined. Through the explicit finite-difference calculations with a third-order total-variation-diminishing (TVD) numerical solution scheme, all results of the theoretical copper and cadmium breakthrough curves (BTCs) simulated by using the Freundlich nonlinear retardation factors revealed good agreement with the experimental observations.     

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.