[1] Recirculating well systems provide an engine for the in situ treatment of subsurface contaminants. Although numerous recirculating wells have been installed in the field, for such systems, there is a paucity of comprehensive monitoring data and models constrained to data appearing in the research literature. Here we present an extensive data set combined with detailed inverse and simulation analyses for a two-well groundwater recirculation system used for in situ bioremediation at Edwards Air Force Base in southern California. The ''conveyor belt'' flow system, which was established for in situ treatment of trichloroethylene (TCE) in two bioactive zones, was created by pumping water upward in one well and downward in another well, each well being screened in both the upper and lower aquifers. A bromide tracer test was conducted and extensively monitored for 60 days. Combined inverse analysis was conducted on hydraulic heads from 38 monitoring wells, 32 bromide concentration histories, and a constraint on the degree of recirculation that was based on TCE concentration data. Four different formulations involving alternative weighting schemes used in a nonlinear weighted least squares simulation-regression analysis were explored. The best formulation provided parameter estimates with tight bounds on estimated covariances, suggesting that the model provides a reasonable description of the hydrogeologic system. Our investigation indicates the geometry of the recirculation zone and the degree of recirculation under two different sets of operating conditions. Surprisingly, our analysis suggests that the effects of aquifer heterogeneity are not significant at this site under the conditions of forced recirculation. Furthermore, anomalous flow through an open monitoring well created significant vertical short-circuiting between the generally insulated aquifers. Flow through this small open conduit was equivalent to as much as 33% of the flow through the pumping wells. Using the model as a guide, we treated the aquifer system and bioactive zones as an equivalent mixed reactor to develop simple expressions relating effluent concentrations to influent concentrations. We demonstrate how these expressions are useful in predicting the removal of TCE that had undergone in situ bioremediation in the recirculatory treatment well system. The finite element model developed in this work serves as the foundation for a reactive transport simulator that we developed to analyze bioremediation which occurred during a 444 day experiment [Gandhi et al., 2002].

[1] We present an analysis of an extensively monitored full-scale field demonstration of in situ treatment of trichloroethylene (TCE) contamination by aerobic cometabolic biodegradation. The demonstration was conducted at Edwards Air Force Base in southern California. There are two TCE-contaminated aquifers at the site, separated from one another by a clay aquitard. The treatment system consisted of two recirculating wells located 10 m apart. Each well was screened in both of the contaminated aquifers. Toluene, oxygen, and hydrogen peroxide were added to the water in both wells. At one well, water was pumped from the upper aquifer to the lower aquifer. In the other well, pumping was from the lower to the upper aquifer. This resulted in a ''conveyor belt'' flow system with recirculation between the two aquifers. The treatment system was successfully operated for a 410 day period. We explore how well a finite element reactive transport model can describe the key processes in an engineered field system. Our model simulates TCE, toluene, oxygen, hydrogen peroxide, and microbial growth/death. Simulated processes include advective-dispersive transport, biodegradation, the inhibitory effect of hydrogen peroxide on biomass growth, and oxygen degassing. Several parameter values were fixed to laboratory values or values from previous modeling studies. The remaining six parameter values were obtained by calibrating the model to 7213 TCE concentration data and 6997 dissolved oxygen concentration data collected during the demonstration using a simulation-regression procedure. In this complex flow field involving reactive transport, TCE and dissolved oxygen concentration histories are matched very well by the calibrated model. Both simulated and observed toluene concentrations display similar high-frequency oscillations due to pulsed toluene injection approximately one half hour during each 8 hour period. Simulation results indicate that over the course of the demonstration, 6.9 kg of TCE was degraded and that in the upper aquifer a region 40 m wide extending 25 m down gradient of the treatment system was cleaned up to less than 100 mg L À1 from initial concentrations of approximately 700 mg L À1 . A smaller region was cleaned up to less than 30 mg L À1 . Simulations indicate that the cleaned up area in the upper aquifer would continue to expand for as long as treatment was continued.

Two technologies in combination, cometabolic bioremediation and in-well vapor stripping, were applied to reduce trichloroethylene (TCE) concentrations in groundwater at a contaminant source area without the need to pump contaminated groundwater to the surface for treatment. The vapor-stripping well reduced source TCE concentrations (as high as 6-9 mg/L) by over 95%. Effluent from the well then flowed to two bioremediation wells, where additional reductions of approximately 60% were achieved. TCE removal was extensively monitored (for research and not regulatory purposes) using an automated system that collected samples about every 45 min at 55 locations over an area of approximately 50 x 60 m2. During 4.5 months of system operation, total TCE mass removal was 8.1 kg, 7.1 kg of which resulted from in-well vapor stripping and 1.0 kg from biotreatment. The system reduced the average TCE concentration of about 3000 microg/L in the source-zone groundwater to about 250 microg/L in water leaving the treatment zone, effecting greater than 92% TCE removal. A 6 month rebound study after system operation ceased found TCE concentrations then increased significantly in the treatment zone due to diffusion from the fractured rock below and perhaps other processes, with mass increases of about 1.5 kg in the lower aquifer and 0.3 kg in the upper aquifer.