Modeling sewage leakage and transport in carbonate aquifer using carbamazepine as an indicator.

The Western Mountain Aquifer (Yarkon-Taninim) of Israel is one of the country's major water resources and partially flows through a karst system. During late winter 2013, maintenance actions were performed on a central sewage pipe that caused sewage to leak into the creek located in the study area. Carbamazepine (CBZ) was used as an indicator for the presence of sewage in the groundwater. The research goal was to develop a mathematical model for quantifying flow and contaminant transport processes in the karst/fractured-porous unsaturated zone and groundwater system. The model was used to simulate CBZ transport during and after an observed sewage leakage event. A quasi-3D dual permeability numerical model represents the 'vadose zone - aquifer' system, by a series of 1D vertical flow and transport equations solved in a variably-saturated zone and by 3D-saturated flow and transport equation in groundwater. The results of simulation showed that after the leakage stopped, significant amounts of CBZ were retained in the porous matrix of the unsaturated zone below the creek. Water redistribution and slow recharge during the dry summer season contributed to a continuous supply of CBZ to the groundwater in the vicinity of the creek and hundreds of meters downstream.

Rainfall-induced fecal indicator organisms transport from manured fields: model sensitivity analysis.

Microbial quality of surface waters attracts attention due to food- and waterborne disease outbreaks. Fecal indicator organisms (FIOs) are commonly used for the microbial pollution level evaluation. Models predicting the fate and transport of FIOs are required to design and evaluate best management practices that reduce the microbial pollution in ecosystems and water sources and thus help to predict the risk of food and waterborne diseases. In this study we performed a sensitivity analysis for the KINEROS/STWIR model developed to predict the FIOs transport out of manured fields to other fields and water bodies in order to identify input variables that control the transport uncertainty. The distributions of model input parameters were set to encompass values found from three-year experiments at the USDA-ARS OPE3 experimental site in Beltsville and publicly available information. Sobol' indices and complementary regression trees were used to perform the global sensitivity analysis of the model and to explore the interactions between model input parameters on the proportion of FIO removed from fields. Regression trees provided a useful visualization of the differences in sensitivity of the model output in different parts of the input variable domain. Environmental controls such as soil saturation, rainfall duration and rainfall intensity had the largest influence in the model behavior, whereas soil and manure properties ranked lower. The field length had only moderate effect on the model output sensitivity to the model inputs. Among the manure-related properties the parameter determining the shape of the FIO release kinetic curve had the largest influence on the removal of FIOs from the fields. That underscored the need to better characterize the FIO release kinetics. Since the most sensitive model inputs are available in soil and weather databases or can be obtained using soil water models, results indicate the opportunity of obtaining large-scale estimates of FIO transport from fields based on publicly available rather than site-specific information.

Effect of treated domestic wastewater on soil physicochemical and microbiological properties.

A main concern with reuse of treated domestic wastewater (DWW) in irrigation is its possible effect on the soil. Few studies have focused on DWW treated in on-site settings, which generally use low-tech systems that can be constructed and serviced locally. One such system is the recirculating vertical flow constructed wetland (RVFCW). The aim of this study was to assess short- to midterm effects of irrigation with DWW treated in the RVFCW. Four groups of plastic barrels, filled with a sandy loam soil, were irrigated for 36 mo with fresh water (FW), FW with added fertilizer, raw DWW, or DWW treated in the RVFCW followed by ultraviolet disinfection. Principal component analysis revealed that the soil irrigated with treated DWW had physicochemical properties similar to those irrigated with FW amended with fertilizer. Levels of surfactants in soil irrigated with treated DWW were identical to those expected from standard irrigation practices, abating concerns for possible changes in soil hydraulic properties. was not detected in the soil irrigated with treated DWW, demonstrating the importance of disinfection of treated effluents before reuse in irrigation. Furthermore, irrigation with treated DWW did not alter the bacterial community structure according to terminal restriction fragment analysis. This 3-yr study suggests that the practice of irrigation with RVFCW effluents is safe. Continuation of the experiment is required to determine whether longer-term irrigation might show a different pattern.

Advection dominated transport of polycyclic aromatic hydrocarbons in amended sediment caps.

Typical sand caps used for sediment remediation have little sorption capacity to retard the migration of hydrophobic contaminants such as PAHs that can be mobilized by significant groundwater flow. Laboratory column experiments were performed using contaminated sediments and capping materials from a creosote contaminated USEPA Superfund site. Azoic laboratory column experiments demonstrated rapid breakthrough of lower molecular weight PAHs when groundwater seepage was simulated through a column packed with coarse sand capping material. After eight pore volumes of flow, most PAHs measured showed at least 50% of initial source pore water concentrations at the surface of 65 cm capping material. PAH concentration in the cap solids was low and comparable to background levels typically seen in urban depositional sediment, but the pore water concentrations were high. Column experiments with a peat amendment delayed PAH breakthrough. The most dramatic result was observed for caps amended with activated carbon at a dose of 2% by dry weight. PAH concentrations in the pore water of the activated carbon amended caps were 3-4 orders of magnitude lower (0.04 ± 0.02 µg/L for pyrene) than concentrations in the pore water of the source sediments (26.2 ± 5.6 µg/L for pyrene) even after several hundred pore volumes of flow. Enhancing the sorption capacity of caps with activated carbon amendment even at a lower dose of 0.2% demonstrated a significant impact on contaminant retardation suggesting consideration of active capping for field sites prone to groundwater upwelling or where thin caps are desired to minimize change in bathymetry and impacts to aquatic habitats.

Colloid transport in porous media: impact of hyper-saline solutions.

The transport of colloids suspended in natural saline solutions with a wide range of ionic strengths, up to that of Dead Sea brines (10(0.9) M) was explored. Migration of microspheres through saturated sand columns of different sizes was studied in laboratory experiments and simulated with mathematical models. Colloid transport was found to be related to the solution salinity as expected. The relative concentration of colloids at the columns outlet decreased (after 2-3 pore volumes) as the solution ionic strength increased until a critical value was reached (ionic strength > 10(-1.8) M) and then remained constant above this level of salinity. The colloids were found to be mobile even in the extremely saline brines of the Dead Sea. At such high ionic strength no energetic barrier to colloid attachment was presumed to exist and colloid deposition was expected to be a favorable process. However, even at these salinity levels, colloid attachment was not complete and the transport of ∼ 30% of the colloids through the 30-cm long columns was detected. To further explore the deposition of colloids on sand surfaces in Dead Sea brines, transport was studied using 7-cm long columns through which hundreds of pore volumes were introduced. The resulting breakthrough curves exhibited a bimodal shape whereby the relative concentration (C/C(0)) of colloids at the outlet rose to a value of 0.8, and it remained relatively constant (for the ∼ 18 pore volumes during which the colloid suspension was flushed through the column) and then the relative concentration increased to a value of one. The bimodal nature of the breakthrough suggests different rates of colloid attachment. Colloid transport processes were successfully modeled using the limited entrapment model, which assumes that the colloid attachment rate is dependent on the concentration of the attached colloids. Application of this model provided confirmation of the colloid aggregation and their accelerated attachment during transport through soil in high salinity solution.

Line-source multi-tracer test for assessing high groundwater velocity.

Segmented line-source multi-tracer injection is suggested as an effective method for assessing groundwater velocities and flow directions in subsurfaces characterized by high water flux. Modifying the common techniques of injecting a tracer into a well became necessary after point-source natural and forced gradient tracer tests ended with no reliable information on the local groundwater flow. The tracer's line-source increases the likelihood of success of the test and could provide additional information regarding the lateral heterogeneity of the aquifer. In a field experiment conducted in the northwestern part on the Dead Sea coast, tracers were injected into an 8-m-long line injection system perpendicular to the assumed flow direction. The injection system was divided into four separate segments with four different tracers. An array of five boreholes located within a 10 × 10 m area downstream was used for monitoring the tracers' transport. Two dye tracers (uranine and Na naphthionate) were injected in a long pulse of several hours into two of the injection pipe segments. Two other tracers (Rhenium oxide and Gd-DTPA) were instantaneously injected into the other two segments. The tracers were detected 0.7 to 2.3 h after injection in four of the five observation wells, located 2.3 to 10 m away from the injection system. The groundwater velocity was determined to be ∼80 to 170 m/d, based on the recoveries of the tracers. The groundwater flow direction was derived based on the arrival of the tracers and was found to be quite consistent with the apparent direction of the hydraulic gradient.

Mathematical model for analysis of recirculating vertical flow constructed wetlands.

The recirculating vertical flow constructed wetland (RVFCW) was developed for the treatment of domestic wastewater (DWW). In this system, DWW is applied to a vertical flow bed through which it trickles into a reservoir located beneath the bed. It is then recirculated back to the root zone of the bed. In this study, a compartmental model was developed to simulate the RVFCW. The model, which addresses transport and removal kinetics of total suspended solids, 5-day biological oxygen demand and nitrogen, was fitted to kinetical results obtained from pilot field setups and a local sensitivity analysis was performed on the model parameters and operational conditions. This analysis showed that after 5h of treatment water quality is affected more by stochastic events than by the model parameter values, emphasizing the stability of the RVFCW system to large variations in operational conditions. Effluent quality after 1h of treatment, when the sensitivity analysis showed the parameter impacts to be largest, was compared to model predictions. The removal rate was found to be dependent on the recirculation rate. The predictions correlated well with experimental observations, leading to the conclusion that the proposed model is a satisfactory tool for studying RVFCWs.

Colloid transport in a heterogeneous partially saturated sand column.

Colloid transport was studied in heterogeneous sand columns under unsaturated steady-state conditions, using two sizes of acid-cleaned sand to pack the column. Heterogeneity was created by placing three continuous tubes of fine sand (3.6% of the total volume) within a column of coarse sand (mean grain diameters 0.36 and 1.2 mm, respectively). Experiments were performed under three flow rates (0.1, 0.2, and 0.4 cm/ min) applied by a rain simulator atthe top of the column. Constant water-content profile in the coarse sand was achieved by applying corresponding suction at the column bottom. Three sizes of latex microspheres (1, 0.2, and 0.02 microm) and soluble tracers (LiBr), diluted in a weak base (pH 7.3, ionic strength 0.0023 M) solution, were used simultaneously. Introduction of preferential pathways reduced front-arrival time about 2-fold and increased colloid recovery which, at the 0.2 cm/min flow rate, was higher than at 0.4 and 0.1 cm/min. Maximum solution flux from coarse to fine sand (due to differences in matric pressure) at 0.2 cm/min, verified by hydrodynamic modeling, could explain this phenomenon. Results suggest that in heterogeneous soil, maximum colloid recovery does not necessarily occur at maximum water content. This has clear implications for colloid transport in natural soils, many of which are heterogeneous.

Dilution tests in a low-permeability fractured aquifer: matrix diffusion effect.

A point dilution test is commonly used in single-borehole tracer experiments designed to determine the Darcy velocity of a formation. This method is based on the concept that, in a borehole, a tracer's concentration declines as a consequence of the water flux. Based on theoretical simulations and field observations, this study indicates that for low-permeability, yet highly porous fractured formations, the common practice of excluding the effect of diffusive mass flux between the dissolved tracer within the borehole and the surrounding matrix may lead to significant errors in the assessment of the Darcy velocity. This conclusion was confirmed by a model adapted to simulate experimental data collected from a tracer test performed in a vertical, large-diameter (25-cm) borehole drilled along a subvertical fracture intersecting a chalk formation.

Evaluation of soil flushing potential for clean-up of desert soil contaminated by industrial wastewater.

The flushing potential of a desert loess soil contaminated by the flame retardant Tetrabromobisphenol A (TBBPA), chloride (Cl(-)) and bromide (Br(-)) was studied in undisturbed laboratory column experiments (20 cm diameter, 45 cm long) and a small field plot (2 x 2 m). While the soluble inorganic ions (Cl(-) and Br(-)) were efficiently flushed from the soil profile after less than three pore volumes (PV) of water, about 50% of the initial amount of TBBPA in the soil was also flushed, despite its hydrophobic nature. TBBPA leaching was made possible due to a significant increase in the pH of the soil solution from 7.5 to 9, which increased TBBPA aqueous solubility. The remaining TBBPA mass in the soil was not mobilized from its initial location in the topsoil due to the decrease in pH at this horizon. In situ soil flushing demonstrated that this method is a feasible treatment for reducing soil contamination at this site.

Two-dimensional distribution of microbial activity and flow patterns within naturally fractured chalk.

The two-dimensional distribution of flow patterns and their dynamic change due to microbial activity were investigated in naturally fractured chalk cores. Long-term biodegradation experiments were conducted in two cores ( approximately 20 cm diameter, 31 and 44 cm long), intersected by a natural fracture. 2,4,6-tribromophenol (TBP) was used as a model contaminant and as the sole carbon source for aerobic microbial activity. The transmissivity of the fractures was continuously reduced due to biomass accumulation in the fracture concurrent with TBP biodegradation. From multi-tracer experiments conducted prior to and following the microbial activity, it was found that biomass accumulation causes redistribution of the preferential flow channels. Zones of slow flow near the fracture inlet were clogged, thus further diverting the flow through zones of fast flow, which were also partially clogged. Quantitative evaluation of biodegradation and bacterial counts supported the results of the multi-tracer tests, indicating that most of the bacterial activity occurs close to the inlet. The changing flow patterns, which control the nutrient supply, resulted in variations in the concentrations of the chemical constituents (TBP, bromide and oxygen), used as indicators of biodegradation.

Biodegradation of 2,4,6-tribromophenol during transport in fractured chalk.

The effect of physicochemical conditions (residence time, oxygen concentrations, and chalk characteristics) on the biodegradation of 2,4,6-tribromophenol (TBP) during transport was investigated in low-permeability fractured-chalk cores. Long-term (approximately 600 d) biodegradation experiments were conducted in two cores (approximately 21 cm diameter, 31 and 44 cm long, respectively), intersected by a natural fracture. TBP was used as a model contaminant and as the sole carbon source for aerobic microbial activity. Bacterial isolates were recovered and identified by both Biolog identification kit and 16S rDNA sequences from batch enrichment cultures. One of the strains, with 98% similarity (based on the 16S rDNA data) to Achromobacter xylosoxidans, was shown to have the ability to degrade TBP in the presence of chalk. The decrease in TBP concentration along the fracture due to biodegradation was not affected by reducing the residence time from 49 to 8 min. In contrast, adding oxygen to the water at the inlet and increasing the flow rates improved TBP removal. Although the matrix pore-size distribution limits microbial activity to the fracture void, the chalk appears to provide an excellent environment for biodegradation activity. Approximately 90% of TBP removal occurred within 10 cm of the TBP source, indicating that in-situ bioremediation can be used to remove organic contaminants in low-permeability fractured rocks if nutrient-delivery pathways within the aquifer are secured.

Impact of microbial activity on the hydraulic properties of fractured chalk.

The impact of microbial activity on fractured chalk transmissivity was investigated on a laboratory scale. Long-term experiments were conducted on six fractured chalk cores (20 cm diameter, 23-44 cm long) containing a single natural fracture embedded in a porous matrix. Biodegradation experiments were conducted under various conditions, including several substrate and oxygen concentrations and flow rates. 2,4,6-Tribromophenol (TBP) was used as a model contaminant (substrate). TBP biodegradation efficiency depended mainly on the amount of oxygen. However, under constant oxygen concentration at the core inlet, elevating the flow rates increased the removal rate of TBP. Transmissivity reduction was clearly related to TBP removal rate, following an initial slow decline and a further sharp decrease with time. The fracture's transmissivity was reduced by as much as 97% relative to the initial value, with no leveling off of the clogging process. For the most extreme cases, reductions of 262 and 157 microm in the equivalent hydraulic apertures were recorded for fractures with initial apertures of 495 and 207 microm, respectively. The reductions in fracture transmissivity occurred primarily because of clogging by bacterial cells and extracellular polymeric substances (EPS) produced by the bacteria. Most of the biodegradation activity was concentrated near the fracture inlet, where the most suitable biodegradation conditions (nutrients and oxygen) prevailed, suggesting that the clogging had occurred in that vicinity. The clogging must have changed the structure of the fracture void, thereby reducing the active volume participating in flow and transport processes. This phenomenon caused accelerated transport of non-reactive tracers and doubled the fracture's dispersivity under constant flow rates.