Filtration and transport of Bacillus subtilis spores and the F-RNA phage MS2 in a coarse alluvial gravel aquifer: implications in the estimation of setback distances

Filtration of Bacillus subtilis spores and the F-RNA phage MS2 (MS2) on a field scale in a coarse alluvial gravel aquifer was evaluated from the authors' previously published data. An advection-dispersion model that is coupled with first-order attachment kinetics was used in this study to interpret microbial concentration vs. time breakthrough curves (BTC) at sampling wells. Based on attachment rates (katt) that were determined by applying the model to the breakthrough data, filter factors (f) were calculated and compared with f values estimated from the slopes of log (cmax/co) vs. distance plots. These two independent approaches resulted in nearly identical filter factors, suggesting that both approaches are useful in determining reductions in microbial concentrations over transport distance. Applying the graphic approach to analyse spatial data, we have also estimated the f values for different aquifers using information provided by some other published field studies. The results show that values of f, in units of log (cmax/co) m(-1), are consistently in the order of 10(-2) for clean coarse gravel aquifers, 10(-3) for contaminated coarse gravel aquifers, and generally 10(-1) for sandy fine gravel aquifers and river and coastal sand aquifers. For each aquifer category, the f values for bacteriophages and bacteria are in the same order-of-magnitude. The f values estimated in this study indicate that for every one-log reduction in microbial concentration in groundwater, it requires a few tens of meters of travel in clean coarse gravel aquifers, but a few hundreds of meters in contaminated coarse gravel aquifers. In contrast, a one-log reduction generally only requires a few meters of travel in sandy fine gravel aquifers and sand aquifers. Considering the highest concentration in human effluent is in the order of 10(4) pfu/l for enteroviruses and 10(6) cfu/100 ml for faecal coliform bacteria, a 7-log reduction in microbial concentration would comply with the drinking water standards for the downgradient wells under natural gradient conditions. Based on the results of this study, a 7-log reduction would require 125-280 m travel in clean coarse gravel aquifers, 1.7-3.9 km travel in contaminated coarse gravel aquifers, 33-61 m travel in clean sandy fine gravel aquifers, 33-129 m travel in contaminated sandy fine gravel aquifers, and 37-44 m travel in contaminated river and coastal sand aquifers. These recommended setback distances are for a worst-case scenario, assuming direct discharge of raw effluent into the saturated zone of an aquifer. Filtration theory was applied to calculate collision efficiency (alpha) from model-derived attachment rates (katt), and the results are compared with those reported in the literature. The calculated alpha values vary by two orders-of-magnitude, depending on whether collision efficiency is estimated from the effective particle size (d10) or the mean particle size (d50). Collision efficiency values for MS-2 are similar to those previously reported in the literature (e.g. ) [DeBorde, D.C., Woessner, W.W., Kiley, QT., Ball, P., 1999. Rapid transport of viruses in a floodplain aquifer. Water Res. 33 (10), 2229-2238]. However, the collision efficiency values calculated for Bacillus subtilis spores were unrealistic, suggesting that filtration theory is not appropriate for theoretically estimating filtration capacity for poorly sorted coarse gravel aquifer media. This is not surprising, as filtration theory was developed for uniform sand filters and does not consider particle size distribution. Thus, we do not recommend the use of filtration theory to estimate the filter factor or setback distances. Either of the methods applied in this work (BTC or concentration vs. distance analyses), which takes into account aquifer heterogeneities and site-specific conditions, appear to be most useful in determining filter factors and setback distances.     

Use of vegetable oil in a pilot-scale denitrifying barrier

Nitrate in drinking water is a hazard to both humans and animals. Contaminated water can cause methemoglobinemia and may pose a cancer risk. Permeable barriers containing innocuous oils, which stimulate denitrification, can remove nitrate from flowing groundwater. For this study, a sand tank (1.1 x 2.0 x 0.085 m in size) containing sand was used as a one-dimensional open-top scale model of an aquifer. A meter-long area near the center of the tank contained sand coated with soybean oil. This region served as a permeable denitrifying barrier. Water containing 20 mg l(-1) nitrate-N was pumped through the barrier at a high flow rate, 1112 l week(-1), for 30 weeks. During the 30-week study, the barrier removed 39% of the total nitrate-N present in the water. The barrier was most efficient during the first 10 weeks of the study when almost all of the nitrate and nitrogen was removed. Efficiency declined with time so that by week 30 almost no nitrate was removed by the system. Nitrite levels in the effluent water remained low throughout the study. Barriers could be used to protect groundwater from nitrate contamination or for the in situ treatment of contaminated water. At the low flow rates that exist in most aquifers, such barriers should be effective at removing nitrate from groundwater for a much longer period of time.     

Modeling potential vadose-zone transport of nitrogen from onsite wastewater systems at the development scale

Water in the urban front-range corridor of Colorado has become an increasingly critical resource as the state faces both supply issues as well as anthropogenic degradation of water quality in several aquifers used for drinking water. A proposed development (up to 1100 homes over two quarter-quarter sections) at Todd Creek, Colorado, a suburb of Westminster located about 20 miles northeast of Denver, is considering use of onsite wastewater systems (OWS) to treat and remove domestic wastewater. Local health and environmental agencies have concerns for potential impacts to local water quality. Nitrogen treatment in the vadose zone and subsequent transport to ground water at a development scale is the focus of this investigation. The numerical model HYDRUS 1D was used, with input based on site-specific data and several transport parameters estimated from statistical distribution, to simulate nitrate concentrations reaching ground water. The model predictions were highly sensitive to mass-loading of nitrogen from OWS and the denitrification rate coefficient. The mass loading is relatively certain for the large number of proposed OWS. However, reasonable values for the denitrification rate coefficients vary over three orders of magnitude. Using the median value from a cumulative frequency distribution function, based on rates obtained from the literature, resulted in simulated output nitrate concentrations that were less than 1% of regulatory maximum concentrations. Reasonable rates at the lower end of the reported range, corresponding to lower 95% confidence interval estimates, result in simulated nitrate concentrations reaching groundwater above regulatory limits.     

Pesticide storage and release in unsaturated soil in Illinois, USA.

The chemical fate and movement of pesticides may be subject to transient storage in unsaturated soils during periods of light rainfall, and subsequent release into shallow groundwater by increased rainfall. The objective of this study was to conduct field-scale experiments to determine the relative importance of transient storage and subsequent release of agrichemicals from the vadose zone into potential aquifers. Two field-scale experiments were conducted under a rain exclusion shelter. In the 1x experiment, atrazine and chlorpyrifos were applied at application-rate equivalents (1.6 kg ha(-1) and 1.3 kg ha(-1), respectively). In the 4x experiment, atrazine was applied in an amount that was four times greater than that usually applied to fields (6.7 kg ha(-1)). Water was either applied to simulate rain or withheld to simulate dry periods. In the 1x experiment, atrazine was detected in the water samples whereas chlorpyrifos was not detected in the majority of the samples. The dry period imposed on the treatment plot did not appear to result in storage of the chemicals, whereas the wet period resulted in greater leaching of atrazine, although the concentrations remained less than the Maximum Contaminant Level of 3 microg L(-1). Both chemicals were detected in soil samples collected from a 20- to 30-cm depth, but it appeared that both chemicals dissipated before the field experiment was concluded. It appeared that the one-time application of atrazine and chlorpyrifos at the label rates did not result in a sufficient mass to be stored and flushed in significant concentrations to the saturated zone. When atrazine was applied at 4x and a longer drought period was imposed on the treatment plot, the resulting concentrations of dissolved atrazine were still less than 3 microg L(-1) . Atrazine was detected in only the near-surface (0 to 15 cm) soil samples and the herbicide dissipated before the onset of the dry period in the treatment plot. The results of this field study demonstrated that atrazine and chlorpyrifos were not sufficiently persistent to be stored and then released in significantly large concentrations to the saturated zone. The dissipation half-life of atrazine in the 4x application was about 44 days. This study, in addition to others, suggested that atrazine may be less persistent in surface soil than has been generally reported.     

Nitrogen dynamics of a mountain forest on dolomitic limestone--a scenario-based risk assessment

The dominant nitrogen (N) fluxes were simulated in a mountain forest ecosystem on dolomitic bedrock in the Austrian Alps. Based on an existing small-scale climate model the simulation encompassed the present situation and a 50-yr projection. The investigated scenarios were current climate, current N deposition (SC1) and future climate (+2.5 degrees C and +10% annual precipitation) with three levels of N deposition (SC2, 3, 4). The microbially mediated N transformation, including the emission of nitrogen oxides, was calculated with PnET-N-DNDC. Soil hydrology was calculated with HYDRUS and was used to estimate the leaching of nitrate. The expected change of the forest ecosystem due to changes of the climate and the N availability was simulated with PICUS. The incentive for the project was the fact that forests on dolomitic limestone stock on shallow Rendzic Leptosols that are rich in soil organic matter are considered highly sensitive to the expected environmental changes. The simulation results showed a strong effect due to increased temperatures and to elevated levels of N deposition. The outflux of N, both as nitrate (6-25kg Nha(-1)yr(-1)) and nitrogen oxides (1-2kg Nha(-1)yr(-1)), from the forest ecosystem are expected to increase. Temperature exerts a stronger effect on the N(2)O emission than the increased rate of N deposition. The main part of the N emission will occur as N(2) (15kg Nha(-1)yr(-1)). The total N loss is partially offset by increased rates of N uptake in the biomass due to an increase in forest productivity.     

The addition of organic carbon and nitrate affects reactive transport of heavy metals in sandy aquifers

Organic carbon introduction in the soil to initiate remedial measures, nitrate infiltration due to agricultural practices or sulphate intrusion owing to industrial usage can influence the redox conditions and pH, thus affecting the mobility of heavy metals in soil and groundwater. This study reports the fate of Zn and Cd in sandy aquifers under a variety of plausible in-situ redox conditions that were induced by introduction of carbon and various electron acceptors in column experiments. Up to 100% Zn and Cd removal (from the liquid phase) was observed in all the four columns, however the mechanisms were different. Metal removal in column K1 (containing sulphate), was attributed to biological sulphate reduction and subsequent metal precipitation (as sulphides). In the presence of both nitrate and sulphate (K2), the former dominated the process, precipitating the heavy metals as hydroxides and/or carbonates. In the presence of sulphate, nitrate and supplemental iron (Fe(OH)(3)) (K3), metal removal was also due to precipitation as hydroxides and/or carbonates. In abiotic column, K4, (with supplemental iron (Fe(OH)(3)), but no nitrate), cation exchange with soil led to metal removal. The results obtained were modeled using the reactive transport model PHREEQC-2 to elucidate governing processes and to evaluate scenarios of organic carbon, sulphate and nitrate inputs.     

Trends in nitrate concentrations and determination of its origin using stable isotopes (18O and 15N) in groundwater of the Western Central Valley, Costa Rica

A study was conducted to evaluate long-term trends in nitrate concentrations and to try to identify the origin of nitrate using stable isotopes (15N(NO3-) and 18O(NO3-)) in the aquifers of the western Central Valley, Costa Rica, where more than 1 million people depend on groundwater to satisfy their daily needs. Data from 20 sites periodically sampled for 4 to 17 years indicate an increasing trend in nitrate concentrations at five sites, which in a period ranging from 10 to 40 years, will exceed recommended maximum concentrations. Results of isotopic analysis indicate a correspondence between land use patterns and the isotopic signature of nitrate in groundwater and suggest that urbanization processes without adequate waste disposal systems, followed by coffee fertilization practices, are threatening water quality in the region. We conclude that groundwater management in this area is not sustainable, and that land use substitution processes from agricultural activity to residential occupation that do not have proper sewage disposal systems may cause a significant increment in the nitrate contaminant load.     

A PAH fate model for San Francisco Bay

A mass balance model was applied to simulate the long-term fate of PAHs in San Francisco Bay. The model treats the Bay as a single box with interacting water and sediment compartments, and includes loading, volatilization, outflow to the ocean, degradation, and burial in deep sediment. The estimated time required for loss of one-half of the mass in the Bay in the absence of loading ranged from 20 d for naphthalene to 5 yr for benzo(b)fluoranthene. Uncertainty analysis using Monte Carlo simulation indicated a high degree of influence and uncertainty for degradation rates, suggesting that improved estimates of degradation would significantly improve the predictive ability of the model. A comparison of model calculations to literature values suggested that external PAH loading to San Francisco Bay was at or above previous estimates of 3600 kgyr(-1), and that degradation in the Bay was within the range of commonly published estimates for high molecular weight PAHs (4.0 x 10(-5) to 4.0 x 10(-4)d(-1)).     

Development of the high spatial resolution EMICAT2000 emission model for air pollutants from the north-eastern Iberian Peninsula (Catalonia, Spain)

Ozone (O(3)) pollution episodes take place in Catalonia (NE of the Iberian Peninsula), mainly during summertime. The complex O(3) behaviour could be understood by using a Chemical Transport Model (CTM). Emission inventories provide the spatial and temporal emissions distribution of the O(3) precursors and other pollutants required by this approach. We developed the EMICAT2000 model with high spatial (cells of 1 km(2)) and temporal (1h) resolutions, to estimate the emissions during the year 2000 from Catalonia. Total annual emissions were 107 kt yr(-1) of NO(x), 137 kt yr(-1) of NMVOC, 267 kt yr(-1) of CO, 65 kt yr(-1) of SO(2), 24 kt yr(-1) of TSP and 32,175 kt yr(-1) of equivalent CO(2). Main NO(x) sources are on-road traffic (58%) and industries (38%). Main NMVOC sources are on-road traffic (36%), vegetation (34%) and use of solvents (13%). Speciation was established according to the Carbon Bond IV mechanism. EMICAT2000 generates directly the data files required for the third generation CTM Models-3/CMAQ.     

Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers. I. Model formulation and performance

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. The objectives of this study are to develop and assess the performance of a mechanistic flow and reactive transport model which couples the most relevant physical, geochemical and biochemical processes involved in wastewater plume evolution in sandy aquifers. The numerical model solves for variably saturated groundwater flow and reactive transport of multiple carbon- and nitrogen-containing species in a three-dimensional porous medium. The reactive transport equations are solved using the Strang splitting method which is shown to be accurate for Monod and first- and second-order kinetic reactions, and two to four times more efficient than sequential iterative splitting. The reaction system is formulated as a fully kinetic chemistry problem, which allows for the use of several special-purpose ordinary differential equation (ODE) solvers. For reaction systems containing both fast and slow kinetic reactions, such as the combined nitrogen-carbon system, it is found that a specialized stiff explicit solver fails to obtain a solution. An implicit solver is more robust and its computational performance is improved by scaling of the fastest reaction rates. The model is used to simulate wastewater migration in a 1-m-long unsaturated column and the results show significant oxidation of dissolved organic carbon (DOC), the generation of nitrate by nitrification, and a slight decrease in pH.     

The value of coastal wetlands for hurricane protection

Coastal wetlands reduce the damaging effects of hurricanes on coastal communities. A regression model using 34 major US hurricanes since 1980 with the natural log of damage per unit gross domestic product in the hurricane swath as the dependent variable and the natural logs of wind speed and wetland area in the swath as the independent variables was highly significant and explained 60% of the variation in relative damages. A loss of 1 ha of wetland in the model corresponded to an average USD 33,000 (median = USD 5000) increase in storm damage from specific storms. Using this relationship, and taking into account the annual probability of hits by hurricanes of varying intensities, we mapped the annual value of coastal wetlands by 1 km x 1 km pixel and by state. The annual value ranged from USD 250 to USD 51,000 ha(-1) yr(-1), with a mean of USD 8240 ha(-1) yr(-1) (median = USD 3230 ha(-1) yr(-1)) significantly larger than previous estimates. Coastal wetlands in the US were estimated to currently provide USD 23.2 billion yr(-1) in storm protection services. Coastal wetlands function as valuable, selfmaintaining "horizontal levees" for storm protection, and also provide a host of other ecosystem services that vertical levees do not. Their restoration and preservation is an extremely cost-effective strategy for society.     

A nitrogen budget of the Changjiang river catchment

Based on 1997-1998 field investigations in the Changjiang river mouth, rain sampling from the river's upper reaches to the mouth, historical data, and relevant literature, the various sources of Total Nitrogen (TN) and Dissolved Inorganic Nitrogen (DIN) in the Changjiang river catchment and N transport in the Changjiang river mouth were estimated. The export fluxes of various form of N were mainly controlled by the river runoff, and the export fluxes of NO3-N, DIN and TN in 1998 (an especially heavy flood year) were 1438 10(3) tonnes (t) yr(-1) or 795.1 kg km(-2) yr(-1), 1746 10(3) t yr(-1) or 965.4 kg km(-2) yr(-1) and 2849 10(3) t yr(-1) or 1575.3 kg km(-2) yr(-1), respectively. The TN and DIN in the Changjiang river came mainly from precipitation, agricultural nonpoint sources, N lost from fertilizer and soil, and point sources of industrial waste and residential sewage discharge, which were about 56.2% and 62.3%, 15.4% and 18.5%, 17.1% and 14.4%, respectively, of the N outflow at the Changjiang river mouth; maximum transport being in the middle reaches.     

Fate of pyrene in contaminated soil amended with alternate electron acceptors

Creosote-contaminated soil samples from the Libby Ground Water Contamination Superfund Site in Libby, MT, were amended with the potential alternate electron acceptors (AEA) nitrate (KNO3), manganese oxide (MnO2), and amorphous iron oxyhydroxide (FeOOH) and incubated at low oxygen tensions (0-6% O2). The fate of 14C-pyrene was evaluated with respect to the different soil amendments. The fate of 14C from the radiolabeled pyrene with regard to mineralization and bound residue formation within soil humic fractions was not significantly different from controls for the iron and manganese amended soils. Nitrate amendments appeared to stimulate 14C-pyrene mineralization at a level of 170 mg NO3-N kg(-1), and inhibit mineralization at 340 mg NO3-N kg(-1). The stimulatory effect did not appear to be the result of nitrate serving as an electron acceptor. Although AEA amendments did not significantly affect the rate or extent of 14C-pyrene mineralization, results of oxygen-deprived incubations (purged with N2) indicate that AEA may be utilized by the microbial community in the unsaturated contaminated soil system.     

Mixing-controlled biodegradation in a toluene plume--results from two-dimensional laboratory experiments

Various abiotic and biotic processes such as sorption, dilution, and degradation are known to affect the fate of organic contaminants, such as petroleum hydrocarbons in saturated porous media. Reactive transport modeling of such plumes indicates that the biodegradation of organic pollutants is, in many cases, controlled by mixing and therefore occurs locally at the plume's fringes, where electron donors and electron-acceptors mix. Herein, we aim to test whether this hypothesis can be verified by experimental results obtained from aerobic and anaerobic degradation experiments in two-dimensional sediment microcosms. Toluene was selected as a model compound for oxidizable contaminants. The two-dimensional microcosm was filled with quartz sand and operated under controlled flow conditions simulating a contaminant plume in otherwise uncontaminated groundwater. Aerobic degradation of toluene by Pseudomonas putida mt-2 reduced a continuous 8.7 mg L(-1) toluene concentration by 35% over a transport distance of 78 cm in 15.5 h. In comparison, under similar conditions Aromatoleum aromaticum strain EbN1 degraded 98% of the toluene infiltrated using nitrate (68.5+/-6.2 mg L(-1)) as electron acceptor. A major part of the biodegradation activity was located at the plume fringes and the slope of the electron-acceptor gradient was steeper during periods of active biodegradation. The distribution of toluene and the significant overlap of nitrate at the plume's fringe indicate that biokinetic and/or microscale transport processes may constitute additional limiting factors. Experimental data is corroborated with results from a reactive transport model using double Monod kinetics. The outcome of the study shows that in order to simulate degradation in contaminant plumes, detailed data sets are required to test the applicability of models. These will have to deal with the incorporation of existing parameters coding for substrate conversion kinetics and microbial growth.     

The long-term environmental behaviour of strontium and barium released from former mine workings in the granites of the Sunart region of Scotland, UK

The concentrations of strontium and barium have been measured in water, sediment and the shells of mussels (Mytilus edulis) from a river system in the Sunart region of Scotland, UK. The aim was to establish the fate and mobility of these elements, which are slowly being released from old mine workings on the Strontian granites. Enhanced strontium (1500-2000 microg l(-1) and 250-290 microg l(-1)) and barium concentrations (316 microg l(-1) and 83 microg l(-1)) were found in the waters originating from the two mine drains studied. Both element were also found at significant levels in the river sediments taken from the vicinity of each drainage site (Sr: 225 microg g(-1) and 120-125 microg g(-1); Ba: 1380 microg g(-1) and 126-170 microg g(-1)). The data suggests that the sediments are acting as a reservoir for these group II cations from where they become distributed throughout the river system. Strontium is found to be incorporated into the shells (3.16-3.46 microg g(-1)) and pearls (3.57 microg g(-1)) of the blue mussel, located at the estuarine margin some 10 km downstream, at values close to the maximum expected (3.3% by weight of the calcium content). The study presents a view of the fate of barium and strontium in a river system over a prolonged period of time. As such it provides valuable information for studies that seek to model the impact of the accidental release of barium and strontium (including the important radionuclide 90Sr) into the environment.     

Comparison of source apportionments of fine particulate matter at two San Jose speciation trends network sites

In this study, the chemical composition of fine particulate matter samples collected at U.S. Environmental Protection Agency Speciation Trends Network sites in San Jose, CA, from February 2000 to February 2005 were analyzed. A San Jose site was initially established at 4th Street and then subsequently moved to Jackson Street in mid-2002. These sites are approximately 1 km apart. There were no known major changes in the nature of the sources in the area over this period. The study used positive matrix factorization model to extract the source profiles and their mass contributions and to compare the results for the congruence of the source apportionments between these two nearby sites. In the case of the 4th Street site, the average mass was apportioned to wood combustion (32.1 +/- 2.5%), secondary nitrate (22.3 +/- 2%), secondary sulfate (10.7 +/- 0.6%), fresh sea salt (7.7 +/- 0.9%), gasoline vehicles (7.3 +/- 0.5%), aged sea salt (6.8 +/- 0.4%), road dust (6.7 +/- 0.7%), diesel emissions (3.9 +/- 0.3%), and a Ni-related industrial source (2.5 +/- 0.4%). At the Jackson Street site, the average mass was apportioned to wood combustion (33.6 +/- 2.6%), secondary nitrate (20.3 +/- 1.9%), secondary sulfate (13.9 +/- 0.9%), aged sea salt (12.4 +/- 0.7%), gasoline vehicle (8.3 +/- 0.6%), fresh sea salt (5.3 +/- 0.5%), diesel emission (3.2 +/- 0.3%), road dust (1.9 +/- 0.1%), and Ni-related industrial source (1.3 +/- 0.1%). Conditional probability function analysis was used to help identify local sources. These results suggested that moving the sampling site a short distance had little effect on the nature of the resolved source types although some differences in their quantitative impacts were obtained in the positive matrix factorization analyses.     

Effects of increased deposition of atmospheric nitrogen on an upland moor: leaching of N species and soil solution chemistry

This study was designed to investigate the leaching response of an upland moorland to long-term (10 yr) ammonium nitrate additions of 40, 80 and 120 kg N ha(-1) yr(-1) and to relate this response to other indications of potential system damage, such as acidification and cation displacement. Results showed increases in nitrate leaching only in response to high rates of N input, in excess of 96 and 136 kg total N input ha(-1) yr(-1) for the organic Oh horizon and mineral Eag horizon, respectively. Individual N additions did not alter ammonium leaching from either horizon and ammonium was completely retained by the mineral horizon. Leaching of dissolved organic nitrogen (DON) from the Oh horizon was increased by the addition of 40 kg N ha(-1) yr(-1), but in spite of increases, retention of total dissolved nitrogen reached a maximum of 92% and 95% of 80 kg added N ha(-1) yr(-1) in the Oh and Eag horizons, respectively. Calcium concentrations and calcium/aluminium ratios were decreased in the Eag horizon solution with significant acidification mainly in the Oh horizon leachate. Nitrate leaching is currently regarded as an early indication of N saturation in forest systems. Litter C:N ratios were significantly lowered but values remained above a threshold predicted to increase leaching of N in forests.     

Effect of water-table fluctuation on dissolution and biodegradation of a multi-component, light nonaqueous-phase liquid

Light nonaqueous-phase liquids (LNAPLs) such as gasoline and diesel fuel are among the most common causes of soil and groundwater contamination. Dissolution and subsequent advective transport of LNAPL components can negatively impact water supplies, while biodegradation is thought to be an important sink for this class of contaminants. We present a laboratory investigation of the effect of a water-table fluctuation on dissolution and biodegradation of a multi-component LNAPL (85% hexadecane, 5% toluene, 5% ethylbenzene, and 5% 2-methylnapthalene on a molar basis) in a pair of similar model aquifers (80 cm x 50 cm x 3 cm), one of which was subjected to a water-table fluctuation. Water-table fluctuation resulted in LNAPL and air entrapment below the water table, an increase in the vertical extent of the LNAPL source zone (by factor 6.7), and an increase in the volume of water passing through the source zone (by factor ~18). Effluent concentrations of dissolved LNAPL components were substantially higher and those of dissolved nitrate lower in the model aquifer where a fluctuation had been induced. Thus, water-table fluctuation led to enhanced biodegradation activity (28.3 mmol of nitrate consumed compared to 16.3 mmol in the model without fluctuation) as well as enhanced dissolution of LNAPL components. Despite the increased biodegradation, fluctuation led to increased elution of dissolved LNAPL components from the system (by factors 10-20). Hence, water-table fluctuations in LNAPL-contaminated aquifers might be expected to result in increased exposure of downgradient receptors to LNAPL components. Accordingly, water-table fluctuations in contaminated aquifers are probably undesirable unless the LNAPL is of minimal solubility or the dissolved-phase plume is not expected to reach a receptor due to distance or the presence of some form of containment.     

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.     

Nutrient dynamics as indicators of karst processes: comparison of the Chalk aquifer (Normandy, France) and the Edwards aquifer (Texas, U.S.A.)

Karst aquifers display a range of geologic and geomorphic characteristics in a wide range of climatic and land-use settings; identification of transport dynamics representative of karst aquifers in general could help advance our understanding of these complex systems. To this end, nutrient, turbidity, and major ion dynamics in response to storms were compared at multiple sites in two karst aquifers with contrasting characteristics and settings: the Chalk aquifer (Eure Department, Normandy, France) and the Barton Springs segment of the Edwards Aquifer (Texas, U.S.A.). The Chalk aquifer is typified by high matrix porosity, thick surficial deposits (up to 30 m thick), and agricultural land use; the Barton Springs segment is typified by low matrix porosity, outcropping limestone, and urban land use. Following one to three storms, from 5 to 16 samples from springs and wells were analyzed for major ions, and specific conductance and turbidity were monitored continuously. Comparison of the chemographs indicated some generalized responses, including an increase in turbidity and potassium concentrations and a decrease in major ion and nitrate concentrations with infiltrating storm runoff. Factor analysis of major ions and turbidity revealed strikingly similar behavior of the chemical variables for the two aquifers: The first two factors, explaining more than 75% of the variability, illustrate that dynamics of most major ions (including nitrate) are opposed to those of turbidity and of potassium. The results demonstrate that potassium and nitrate are effective tracers of infiltrating storm runoff and resident ground water, respectively, and the similar results for these two highly contrasting aquifers suggest that the dynamics identified might be applicable to karst systems in general.