Modelling the migration of contaminants through variably saturated dual-porosity, dual-permeability chalk

In the Hesbaye region in Belgium, tracer tests performed in variably saturated fissured chalk rocks presented very contrasting results in terms of transit times, according to artificially controlled water recharge conditions prevailing during the experiments. Under intense recharge conditions, tracers migrated across the partially or fully saturated fissure network, at high velocity in accordance with the high hydraulic conductivity and low effective porosity (fracture porosity). At the same time, a portion of the tracer was temporarily retarded in the almost immobile water located in the matrix. Under natural infiltration conditions, the fissure network remained inactive. Tracers migrated downward through the matrix, at low velocity in relation with the low hydraulic conductivity and the large porosity of the matrix. Based on these observations, Brouyère et al. (2004a) [Brouyère, S., Dassargues, A., Hallet, V., 2004a. Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation, J. Contam. Hydrol., 72 (1-4), 135-164, doi: 10.1016/j.conhyd.2003.10.009] proposed a conceptual model in order to explain the migration of solutes in variably saturated, dual-porosity, dual-permeability chalk. Here, mathematical and numerical modelling of tracer and contaminant migration in variably saturated fissured chalk is presented, considering the aforementioned conceptual model. A new mathematical formulation is proposed to represent the unsaturated properties of the fissured chalk in a more dynamic and appropriate way. At the same time, the rock water content is partitioned between mobile and immobile water phases, as a function of the water saturation of the chalk rock. The groundwater flow and contaminant transport in the variably saturated chalk is solved using the control volume finite element method. Modelling the field tracer experiments performed in the variably saturated chalk shows the adequacy and usefulness of the new conceptual, mathematical and numerical model.     

Chemical modifications of groundwater contaminated by recharge of treated sewage effluent

Long-term monitoring of the chemical composition of recharge sewage effluent and associated contaminated groundwater from the Dan Region Sewage Reclamation Project shows, after 16 years of recharge operation, the presence of a distinct saline plume (up to 400 mg/l Cl), extending 1600 m downgradient in the Coastal Plain aquifer of Israel. The recorded electrolyte composition of groundwater in the vicinity of the recharge area reflects the variations in the compositions of the sewage effluents, as well as water-rock interactions induced by the recharge of treated sewage effluents. The original sewage composition was modified, particularly during early stages of effluent migration in the unsaturated zone, by cation-exchange and adsorption reactions. Since the soil sorption capacity is finite these reactions caused only limited modifications, and once the system reached a steady state the inorganic composition of the contaminated groundwater became similar to that of the recharge water. Decomposition of organic matter in the unsaturated zone resulted in CO₂ generation and dissolution of CaCO₃ minerals in the aquifer. It was shown that chemical and/or bio-degradation of organic matter takes place mainly in the unsaturated (vadose) zone. Hence, monitoring the efficiency of the vadose zone to retain contaminants is essential for evaluating the quality of groundwater since it was shown that organic compounds behave almost conservatively once the effluents enter and flow within the saturated zone.     

非正规垃圾填埋场包气带介质的污染物阻滞能力研究

非正规垃圾填埋场已成为地下水的潜在污染源,给城市供水安全带来了极大风险。以北京某非正规垃圾填埋场的3种典型包气带介质为研究对象,以Br-为示踪剂,采取土柱实验和数值模拟相结合的方法,通过对场区3种典型包气带介质饱和与非饱和状态下的Br-迁移速率的比较,探讨包气带介质对污染物的阻滞能力。结果表明,Br-在不同包气带介质饱和状态下的迁移速率为中砂>细砂>砂质粉土,这可能与不同介质中粉粒和粘粒的含量有关,同时介质的有机质含量也会影响Br-的迁移速率,但前者因素占了主导作用;在Br-迁移过程中,不同包气带介质的含水量是影响其迁移速率的主要因素,若能合理控制场区的地下水位,使得包气带介质长期处于非饱和状态,将会大大延缓垃圾渗滤液中污染物向地下水迁移;针对特定的非正规垃圾填埋场,可采取抽水等措施,控制场区的细砂和砂质粉土层处于非饱和状态,充分利用包气带介质在非饱和状态时较强的污染物阻滞能力,阻隔或延缓垃圾渗滤液中的污染物进入地下水中。     

A method for determining volatile organic solvents in chalk pore waters (Southern and Eastern England) and its relevance to the evaluation of groundwater contamination

Conventional methods of sampling groundwater in Chalk aquifers sample the mobile water in the fissures only and not the non-mobile water held in the microporous matrix. However, the microporous matrix can play an important role in modifying groundwater quality by attenuating contaminant migration through diffusion exchange and, in the long term, by acting as a reservoir for contaminant retention.A method for extracting and analysing chalk pore waters for chlorinated solvents has been developed. Chalk pore water-depth profiles at the site of a spillage of tetrachloroethene demonstrated that the immiscible phase of the solvent had migrated to depths of 50 m, producing solvent concentrations in the matrix pore waters of up to 40,000 μg L⁻¹.In addition, the pore water-depth profile indicated that the immiscible phase had ponded upon a less permeable horizon within the Chalk aquifer at shallow depths.     

Sorption of hydrophobic organic pollutants in saturated soil systems

Sorption equilibria and rates were characterized for a matrix of four aquifer sands and two slightly to moderately hydrophobic organic solutes (nitrobenzene and lindane), and the effects of sorption on the behavior of these solutes in saturated systems of the soils were determined. Experimental data were used to test and evaluate a variety of mathematical models for predicting contaminant fate and transport in groundwater systems.Observed equilibrium relationships between soil and solution phase solute concentrations were found to be described best by the nonlinear Freundlich isotherm model. It was further determined that the sorption process in the systems tested is rate controlled, requiring several days to approach equilibrium in completely mixed batch reactors. Subsequent modeling of solute transport in continuous flow soil column reactors was found to be most successful when rate-controlled models were used, the best results were obtained with a dual-resistance model incorporating the coupled mass transport steps of boundary-layer and intraparticle diffusion.     

Impact of urbanization on hydrochemical evolution of groundwater and on unsaturated-zone gas composition in the coastal city of Tel Aviv, Israel

The coastal city of Tel Aviv was founded at the beginning of the 20th century. The number of its inhabitants and its water consumption increased rapidly. This study analyses a 15-year record (1934-1948) of pre-industrial development of groundwater chemistry in the urban area. Archive data on concentrations of major ions, dissolved gases (CO2 and O2), organic matter, and pH were available for each half-year during the period of 1934-1948. The major factors causing changes in the chemistry of groundwater flowing in three sandy sub-aquifers have been seawater encroachment due to overpumping, and infiltration of effluents from pit-latrine collectors. Influence of these factors decreases with depth. Landward-penetrating seawater passed through clayey coastal sediments, interbedded among sands and calcareous sandstones, and spread into the Kurkar Group aquifer. This has led to exchange of sodium (dominant in seawater) with calcium adsorbed on clay particles, enriching groundwater with calcium. Intensity of cation exchange decreases inland and with depth. Infiltration of pit-latrine effluents has introduced large amounts of ammonium into the unsaturated zone. Its rapid oxidation in unsaturated sediments has caused massive nitrate production, accompanied by pore-water acidification. This process induces dissolution of vadose carbonate, resulting in enrichment of groundwater recharge in calcium. Anthropogenically induced dissolution of calcite in the unsaturated zone has been the major factor for the increase of Ca2+ concentration in groundwater, accounting for about 80% of this increase. In the interface zone, an additional 20% of calcium has been supplied by cation exchange. Owing to pH increase caused by denitrification in the aquifer, Ca(2+)-rich waters supersaturated with calcite could be formed, especially in the capillary fringe of the uppermost sub-aquifer, which could induce calcite precipitation and ultimately lead to the cementation of sandy aquifers. Urban development has caused drastic changes in the gas content in the unsaturated zone and in groundwater. Carbon dioxide was intensively generated by nitrification-denitrification processes, by hydration of urea, to a lesser degree by oxidation of organic matter, and probably by anoxic biodegradation of organics. Between 1934 and 1948, concentrations of CO2 in unsaturated sediment air rose from 3.2% to 7.6%. In the unsaturated zone, oxygen consumption for oxidation of ammonium and organic matter lowered O2 concentrations in sediment air to unusually low values of 3.9-12.9%. Nitrification in the urban unsaturated zone could thus serve as a pump, sucking in atmospheric oxygen at a rate of about 0.3-0.5 g m-2 day-1. The extreme concentrations of CO2 and O2 in unsaturated sediments have been preserved due to production and consumption of gas under conditions of diminishing areas open to the atmosphere, uncovered by buildings and by roads.     

Penetration of herbicides to groundwater in an unconfined chalk aquifer following normal soil applications

The persistence and penetration of the herbicides isoproturon and chlorotoluron in an unconfined chalk aquifer has been monitored over a 4-year period through soil sampling, shallow coring and groundwater monitoring. Chlorotoluron was applied on plots as a marker compound, having never been used previously on that, or surrounding fields. The fieldsite had a 5 degree slope with soil depths of 0.5 to 1.5 m and a water table between 20 and 5 m from the soil surface. Where the water table was deepest (9-20 m below surface (mbs)) little or no positive herbicide detections were made. However, where the water table was at only 4-5 mbs, a regular pesticide signal of around 0.1 microg/l for isoproturon and chlorotoluron could be distinguished. Over the winter recharge period automatic borehole samplers revealed a series of short-lived peaks of isoproturon and chlorotoluron reaching up to 0.8 microg/l. This is consistent with a preferential flow mechanism operating at this particular part of the field. Such peaks were occurring over 2 years after the last application of these compounds. Shallow coring failed to uncover any significant pesticide pulse moving through the deep unsaturated zone matrix at the fieldsite.     

Distribution of metals in a polluted aquifer: A comparison of aquifer suspended material to fine sediments of the adjacent environment

Metal concentrations were determined for groundwater suspended matter from a site in the coastal aquifer of Israel which has been irrigated with secondary sewage effluents since the 1960's. Suspended matter was collected from the aquifer saturated zone by pumping and by a multi-layer sampler. Fine sediments were collected from both the unsaturated and saturated zones of the contaminated aquifer, as well as from an adjacent uncontaminated environment. Ag, Cu, Fe, Mn and Zn were leached from the samples in three sequential chemical extractions which are taken to represent the carbonate, organic and oxide phases. Comparison of the aquifer samples to those of the adjacent environment showed that Fe and Mn are primarily enriched in non-mobile fine sediments and not in suspended matter, whereas the concentrations of Zu, Cu and Ag show up to an order of magnitude enrichment in the mobile suspended matter in groundwater. The enrichment of these metals in the suspended matter indicates that metals from sewage effluents and agricultural activities have reached the groundwater.     

Effects of a reactive barrier and aquifer geology on metal distribution and mobility in a mine drainage impacted aquifer

The Nickel Rim aquifer has been impacted for five decades by a metal-rich plume generated from the Nickel Rim mine tailings impoundment. Metals released by the oxidation of pyrrhotite in the unsaturated zone of the tailings migrate into the downgradient aquifer, affecting both the groundwater and the aquifer solids. A reactive barrier has been installed in the aquifer to remove sulfate and metals from the groundwater. The effect of the reactive barrier on metal concentrations in the aquifer solids has not previously been studied. In this study, a series of selective extraction procedures was applied to cores of aquifer sediment, to ascertain the distribution of metals among various solid phases present in the aquifer. Extraction results were combined with groundwater chemistry, geochemical modelling and solid-phase microanalyses, to assess the potential mobility of metals under changing geochemical conditions. Reactions within the reactive barrier caused an increase in the solid-phase carbonate content downgradient from the barrier. The concentrations of poorly crystalline, oxidized phases of Mn and Fe, as well as concentrations of Cr(III) associated with oxidized Fe, and poorly crystalline Zn, are lower downgradient from the barrier, whereas total solid-phase metal concentrations remain constant. Iron and Mn accumulate as oxidized, easily extractable forms in a peat layer overlying the aquifer. Although these oxides may buffer reducing plumes, they also have the potential to release metals to the groundwater, should a reduced condition be imposed on the aquifer by remedial actions.     

Particle transport in unsaturated fractured chalk under arid conditions

A series of field and laboratory experiments were conducted to study the mechanisms of particle detachment and transport from fractures in vadose chalk. Experiments of intermittent flow events along fracture surfaces were carried out in the laboratory. In the field, water was percolated from land surface via a discrete fracture into a compartmental sampler installed inside a horizontal corehole located I m below the surface. The mass, size distribution, and composition of the particles drained from the fracture voids were examined along with flow rates and salt dissolution. Two boreholes penetrating the underlying saturated zone were sampled and analyzed for colloidal concentration and composition. Most of the particle and solute release at the drained effluents occurred during the first several hours of flow, but erratic pulses of particles were still observed after long periods of time. Most of the detached particles had a mean diameter of >2 microm, while the mobile colloidal phase in the groundwater had a mean diameter of approximately 1 microm. Mineralogical composition of the groundwater colloids and the particles detached from the upper vadose fracture were similar. Laboratory observations demonstrated the importance of the existence of a coating layer, made of weathered particles and salts, on particle detachment. The results of this study suggest that: (1) particle detachment causes flow-rate variability in the unsaturated fracture; (2) the mechanisms of particle detachment and salt dissolution within the fracture are linked: and (3) although most of the detached particles are large and likely to accumulate inside fractures, some colloidal particles also eroded from the fracture void and are likely to be transported to the groundwater.     

Estimating spatially-variable first-order rate constants in groundwater reactive transport systems

Numerical reactive transport models are often used as tools to assess aquifers contaminated with reactive groundwater solutes as well as investigating mitigation scenarios. The ability to accurately simulate the fate and transport of solutes, however, is often impeded by a lack of information regarding the parameters that define chemical reactions. In this study, we employ a steady-state Ensemble Kalman Filter (EnKF), a data assimilation algorithm, to provide improved estimates of a spatially-variable first-order rate constant λ through assimilation of solute concentration measurement data into reactive transport simulation results. The methodology is applied in a steady-state, synthetic aquifer system in which a contaminant is leached to the saturated zone and undergoes first-order decay. Multiple sources of uncertainty are investigated, including hydraulic conductivity of the aquifer and the statistical parameters that define the spatial structure of the parameter field. For the latter scenario, an iterative method is employed to identify the statistical mean of λ of the reference system. Results from all simulations show that the filter scheme is successful in conditioning the λ ensemble to the reference λ field. Sensitivity analyses demonstrate that the estimation of the λ values is dependent on the number of concentration measurements assimilated, the locations from which the measurement data are collected, the error assigned to the measurement values, and the correlation length of the λ fields.     

Source and persistence of pesticides in a semi-confined chalk aquifer of southeast England

Pesticide contamination in groundwater is an increasing problem that poses a significant long-term threat to water quality. Following the detection of elevated concentrations of diuron in boreholes in a semi-confined chalk aquifer from southeast England, a sampling programme was undertaken. Between 2003 and 2004 diuron was observed in 90% of groundwaters analysed. In 60% of groundwater samples metabolites of diuron were more prevalent than the parent compound. Longer-term (1989-2005) monitoring shows that pollution of the aquifer by atrazine, simazine, and more recently diuron, shows a positive correlation with periods of high groundwater levels. Results from groundwater residence time indicators suggest that the highest diuron concentrations are associated with waters containing the greatest proportion of recent recharge. There is some evidence to indicate that diuron occurrence can be spatially related to areas of urban and industrial development and is probably correlated with amenity usage.     

Competitive sorption of organic contaminants in chalk

In the Negev desert, Israel, a chemical industrial complex is located over fractured Eocene chalk formations where transfer of water and solutes between fracture voids and matrix pores affects migration of contaminants in the fractures due to diffusion into the chalk matrix. This study tests sorption and sorption competition between contaminants in the chalk matrix to make it possible to evaluate the potential for contaminant attenuation during transport in fractures. Single solute sorption isotherms on chalk matrix material for five common contaminants (m-xylene, ametryn, 1,2-dichloroethane, phenanthrene, and 2,4,6-tribromophenol) were found to be nonlinear, as confirmed in plots of Kd versus initial solution concentration. Over the studied concentration ranges, m-xylene Kd varied by more than a factor of 100, ametryn Kd by a factor of 4, 1,2-dichloroethane Kd by more than a factor of 3, phenanthrene Kd by about a factor of 2, and 2,4,6-tribromophenol Kd by a factor of 10. It was earlier found that sorption is to the organic matter component of the chalk matrix and not to the mineral phases (Chemosphere 44 (2001) 1121). Nonlinear sorption isotherms indicate that there is at least some finite sorption domain. Bi-solute competition experiments with 2,4,6-tribromophenol as the competitor were designed to explore the nature of the finite sorption domain. All of the isotherms in the bi-solute experiments are more linear than in the single solute experiments, as confirmed by smaller variations in Kd as a function of initial solution concentration. For both m-xylene and ametryn, there is a small nonlinear component or domain that was apparently not susceptible to competition by 2,4,6-tribromophenol. The nonlinear sorption domain(s) is best expressed at low solution concentrations. Inert-solvent-normalized single and bi-solute sorption isotherms demonstrate that ametryn undergoes specific force interactions with the chalk sorbent. The volume percent of phenanthrene sorbed at the liquid solubility limit is calculated to be 13% v:v in both the single and bi-solute experiments. This value exceeds what may be reasonably interpreted as partitioning of phenanthrene into organic matter, despite the relative linearity of the phenanthrene sorption isotherm (compared with other compounds) in both single and bi-solute systems.     

Multiphase flow and transport in fractured clay/sand sequences.

A numerical model (Queen's University Multi-Phase Flow Simulator, QUMPFS) was used to assess the rate of trichloroethylene (TCE) dense, non-aqueous phase liquid (DNAPL) migration through fractured clay, with special attention focused on the influence of interbedded sand lenses. The presence of these sand lenses was found to increase the time required for the non-wetting phase to migrate through the full 30 m vertical extent of the clay sequence from a few days to several years. Applied vertical hydraulic gradients were found to be moderately influential in systems consisting solely of fractured clays, yet one of the dominant factors controlling speed of vertical migration when sand lenses were present. Larger displacement pressure of the sands relative to that of the fractures leads to slower DNAPL migration rates, due to the delays that occur during build-up of capillary pressures. Dissolution of DNAPL and subsequent matrix diffusion of the aqueous phase has little effect on the rate of DNAPL migration through systems consisting of fractured clay only, yet slows the rate of migration in systems containing sand lenses. In all cases examined, the rate of DNAPL loading to the lower aquifer far exceeded the rate of aqueous phase mass loading. It was also found that DNAPL reaches the lower aquifer at approximately the same time as the aqueous phase plumes even for systems experiencing downward groundwater flow due to the attenuation of the aqueous phase through matrix diffusion.     

地下水中轻质有机污染物(LNAPL)透镜体研究

在二维砂槽模型中模拟了轻质油在均匀多孔介质地下水非饱和区中的运移过程。模拟结果表明,地下水毛细区是轻质油污染的重点区,除了ＬＮＡＰＬ的残留以外,进入地下水饱和的ＬＮＡＰＬ终将被地下水顶托回到毛细区中,毛细区以上的约大多民将在重力作用下进入毛细区中,试验中观察到达稳定状态时ＬＮＡＰＬ透镜体的上边缘略微高出毛细区。利用多孔介质毛细管模型,建立了利用界面张力、接触角、介质特征孔隙直径等物理量估算不同位置     

Variability of aquifer sorption properties in a field experiment on groundwater transport of organic solutes: Methods and preliminary results

Methods to characterize the organic solute sorption distribution coefficient, organic carbon content, and specific surface area of aquifer solids from the site of a field experiment on solute transport in groundwater were refined for application to small subsamples of 10-cm depth increments taken from 5-cm diameter cores. Initial results indicate that the average sorption characteristics of the Borden aquifer do not vary appreciably along the trajectory of the solute plumes. However, the sorption distribution coefficient of tetrachloroethylene varied over nearly an order of magnitude among 10-cm depth increments in one core sample. Preliminary evidence suggests that the sorption distribution coefficients for four halogenated organic solutes vary proportionally among core strata. However, the distribution coefficients for sorption of tetrachloroethylene on various depth increments are not well correlated with either organic carbon content or specific surface area, suggesting that as yet unidentified mineral phases may play a significant role in sorption of such solutes by the sandy aquifer solids.     

Transport of Escherichia coli and solutes during waste water infiltration in an urban alluvial aquifer

Recharge of waste water in an unconsolidated poorly sorted alluvial aquifer is a complex process, both physically and hydrochemically. The aim of this paper is to analyse and conceptualise vertical transport mechanisms taking place in an urban area of extensive wastewater infiltration by analysing and combining the water balance, the microbial (Escherichia coli) mass balance, and the mass balance for dissolved solutes. For this, data on sediment characteristics (grain size, organic carbon, reactive iron, and calcite), groundwater levels, and concentrations of E. coli in groundwater and waste water were collected. In the laboratory, data on E. coli decay rate coefficients, and on bacteria retention characteristics of the sediment were collected via column experiments. The results indicated that shallow groundwater, at depths of 50 m below the surface, was contaminated with E. coli concentrations as high as 10(6) CFU/100 mL. In general, E. coli concentrations decreased only 3 log units from the point of infiltration to shallow groundwater. Concentrations were lower at greater depths in the aquifer. In laboratory columns of disturbed sediments, bacteria removal was 2-5 log units/0.5 cm column sediment. Because of the relatively high E. coli concentrations in the shallow aquifer, transport had likely taken place via a connected network of pores with a diameter large enough to allow bacterial transport instead of via the sediment matrix, which was inaccessible for bacteria, as was clear from the column experiments. The decay rate coefficient was determined from laboratory microcosms to be 0.15 d(-1). Assuming that decay in the aquifer was similar to decay in the laboratory, then the pore water flow velocity between the point of infiltration and shallow groundwater, coinciding with a concentration decrease of 3 log units, was 0.38 m/d, and therefore, transport in this connected network of pores was fast. According to the water balance of the alluvial aquifer, determined from transient groundwater modelling, groundwater flow in the aquifer was mainly in vertical downward direction, and therefore, the mass balance for dissolved solutes was simulated using a 1D transport model of a 200 m column of the Quaternary Alluvium aquifer. The model, constructed with PHREEQC, included dual porosity, and was able to adequately simulate removal of E. coli, cation-exchange, and nitrification. The added value of the use of E. coli in this study was the recognition of relatively fast transport velocities occurring in the aquifer, and the necessity to use the dual porosity concept to investigate vertical transport mechanisms. Therefore, in general and if possible, microbial mass balances should be considered more systematically as an integral part of transport studies.     

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.     

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.     

Numerical simulation of a fine-grained denitrification layer for removing septic system nitrate from shallow groundwater

One of the most common methods to dispose of domestic wastewater involves the release of septic effluent from drains located in the unsaturated zone. Nitrogen from such systems is currently of concern because of nitrate contamination of drinking water supplies and eutrophication of coastal waters. It has been proposed that adding labile carbon sources to septic distribution fields could enhance heterotrophic denitrification and thus reduce nitrate concentrations in shallow groundwater. In this study, a numerical model which solves for variably saturated flow and reactive transport of multiple species is employed to investigate the performance of a drain field design that incorporates a fine-grained denitrification layer. The hydrogeological scenario simulated is an unconfined sand aquifer. The model results suggest that the denitrification layer, supplemented with labile organic carbon, may be an effective means to eliminate nitrogen loading to shallow groundwater. It is also shown that in noncalcareous aquifers, the denitrification reaction may provide sufficient buffering capacity to maintain near neutral pH conditions beneath and down gradient of the drain field. Leaching of excess dissolved organic carbon (DOC) from the denitrification layer is problematic, and causes an anaerobic plume to develop in simulations where the water table is less than 5-6 m below ground surface; this anaerobic plume may lead to other down gradient changes in groundwater quality. A drain field and denitrification layer of smaller dimensions is shown to be just as effective for reducing nitrate, but has the benefit of reducing the excess DOC leached from the layer. This configuration will minimize the impact of wastewater disposal in areas where the water table is as shallow as 3.5 m.