Use of Boron as a Tracer for Recycled Water in Brackish Aquifers

An investigation was conducted to develop a simple method for tracking the fate and transport of recycled water following recharge into a shallow brackish aquifer (caprock aquifer) in a coastal area on the Island of Oahu, Hawaii. Several naturally present chemical constituents including the boron isotopic signature (δ?11B) were used to characterize each of the caprock aquifer source waters and the recycled water. Because of the influence of seawater, only δ?11B could be used to clearly distinguish the recycled water from the source waters and the caprock aquifer water. Estimates of the mixing ratios of source waters in the aquifer were made and a method was developed to determine the fraction of wastewater present in a brackish water sample recovered from a monitoring well during recharge operations without addition of a tracer. This method can be adapted to monitor any other brackish aquifer subjected to wastewater recharge.     

Fate of Polybrominated Diphenyl Ethers, Nonylphenol, and Estrogenic Activity during Managed Infiltration of Wastewater Effluent

About a billion cubic meters of wastewater effluent are artificially recharged annually in the United States for maintenance of groundwater levels and prevention of seawater intrusion. There is concern that trace contaminants, including various endocrine disrupting compounds (EDCs), are not completely removed during infiltration, leading to deterioration of groundwater quality. In this work, we investigate the mechanisms and efficiency of EDC removal at the Sweetwater Recharge Facility, which is used to recharge secondary effluent from the Roger Road Wastewater Treatment Plant in Tucson, Ariz. Material was collected from the top meter of sediments in two infiltration basins and analyzed for extractable nonylphenol (NP), polybrominated diphenyl ethers (PBDE) and total estrogenic activity. The basins differed significantly in length of service (7 versus 15 years). Nevertheless, profiles of extractable contaminants and estrogenic activity were similar in the two basins. Results suggest that hydrophobic determinants of estrogenic activity are efficiently retained in surface sediments during soil-aquifer treatment. However, measurable levels of PBDEs, NP, and estrogenic activity are present in infiltrate that reaches the local unconfined aquifer at ～ 40?m below land surface.     

Water Table Fluctuation in the Presence of a Time-Varying Exponential Recharge and Depth-Dependent ET in a Two-Dimensional Aquifer System with an Inclined Base

An analytical solution is presented for water table fluctuation between ditch drains in presence of exponential recharge and depth-dependent evapotranspiration (ET) from groundwater table in a two-dimensional gently sloping aquifer. The groundwater head above the drain is small compared to the saturated thickness of the aquifer. A sound mathematical transformation is devised to transform the two-dimensional groundwater flow equation into a simple form, which makes possible to obtain an analytical solution. The transient midpoint water table variations from the proposed solution compare well with the already existing solutions for horizontal aquifer. A numerical example is used to illustrate the combined effect of depth-dependent ET coupled with a time-varying exponential recharge on the water table fluctuation. The inclusion of a depth-dependent ET in the solution results in water table decline at a faster rate as compared to the case when ET is not considered. With an increase in slope of the aquifer base, water table profiles become asymmetric and the water table divide shifts towards the lower drain. The height of the water table profiles increases on moving away from the boundary of the aquifer and the highest level of the ground water table is obtained in the central portion of the aquifer basin due to the presence of drainage ditches on the aquifer boundary. When the effect of ET is incorporated in combination with recharge, the analytical solution results in accurate and reliable estimates of water table fluctuations under situations subjected to a number of controlling factors. This study will be useful for alleviation of drainage problems of the aquifers receiving surface recharge and surrounded by streams.     

Wastewater Reuse Effects on Soil Hydraulic Conductivity

The wastewater total suspended solids (TSS) concentration effects on the saturated hydraulic conductivity, Ks, of a clay and a loam soil were investigated on laboratory repacked soil cores by a constant head permeameter. Both municipal wastewater (MW) and artificial wastewater (AW) with different TSS concentrations were used, with the aim to evaluate, by comparison, the effects of biological activity. The development of a surface sealed layer was investigated in loam soil columns supplied with AW and equipped with water manometers at different depths to detect the hydraulic head gradient changes. In the loam soil, Ks reduced to about 80% of the initial value after infiltration of 175?mm of MW with TSS = 57–68?mg?L?1. Reductions in Ks were more remarkable in the clay soil. An empirical relationship was proposed to predict the relative hydraulic conductivity, Kr, i.e., the ratio between actual and initial hydraulic conductivity versus the cumulative density loading of TSS. Hydraulic head gradients in the top layer (0–20?mm) of the soil columns increased during application of AW, as a consequence of the formation of a sealed layer, denoting that the surface pore sealing was the main mechanism responsible for the observed Ks reductions. Laboratory data were gathered in a numerical simulation code specifically created to assess the consequences of Ks reduction on water movement through the soil profile. Simulation of both ponded and sprinkler irrigation with MW resulted in reduced infiltration and increased surface ponding condition compared to the application of fresh water (FW).     

Effects of Stormwater Infiltration on Quality of Groundwater Beneath Retention and Detention Basins

Infiltration of storm water through detention and retention basins may increase the risk of groundwater contamination, especially in areas where the soil is sandy and the water table shallow, and contaminants may not have a chance to degrade or sorb onto soil particles before reaching the saturated zone. Groundwater from 16 monitoring wells installed in basins in southern New Jersey was compared to the quality of shallow groundwater from 30 wells in areas of new-urban land use. Basin groundwater contained much lower levels of dissolved oxygen, which affected concentrations of major ions. Patterns of volatile organic compound and pesticide occurrence in basin groundwater reflected the land use in the drainage areas served by the basins, and differed from patterns in background samples, exhibiting a greater occurrence of petroleum hydrocarbons and certain pesticides. Dilution effects and volatilization likely decrease the concentration and detection frequency of certain compounds commonly found in background groundwater. High recharge rates in storm water basins may cause loading factors to be substantial even when constituent concentrations in infiltrating storm water are relatively low.     

Biodegradation of Aqueous Organic Matter over Seasonal Changes: Bioreactor Experiments with Indigenous Lake Water Bacteria

Artificial groundwater recharge for drinking water production involves infiltration of surface water through sandy soil and its capture into a groundwater aquifer. The transformation of aqueous organic matter is one of the central issues in this process. The purpose of this work was to assess the potential of indigenous microorganisms in the source water to contribute in the aqueous organic matter biodegradation. For this purpose, microorganisms were enriched from the source water in a fluidized-bed reactor (FBR) and used for kinetic studies on biodegradation of organic matter at ambient temperature range. Lake water (total organic carbon 5.8?mg?L?1) was continuously fed to the FBR containing porous carrier material to support biomass retention. In the inlet and outlet water there were on average 21±6 and 13±5×105?cells?mL?1, respectively. Biofilm accumulation (as volatile solids) reached 13.1?mg?g?1 dw carrier. In the continuous-flow mode and the batch tests, the highest oxygen consumption rate appeared in the summer, followed by the fall, spring, and winter. At low temperatures, the biodegradation of aqueous organic matter was relatively rapid initially for labile fractions followed by a slower phase for refractory fractions. The average temperature coefficient (Q10) in the system was 2.3 illustrating a strong temperature dependency of oxygen consumption. The isotopic analysis of dissolved inorganic carbon δ13CDIC analysis revealed 27 and 69% mineralizations of dissolved organic carbon at 23 and 6°C over 65 and 630 min, respectively. These results can be used to construct additional input parameters in modeling applications of artificial groundwater recharge process. The biological component especially, i.e., the biodegradation, is difficult to predict for on-site applications without experimental proof and thus the interpretation in this study will help formulate design predictions for the process.     

Storm-Water Infiltration and Focused Recharge Modeling with Finite-Volume Two-Dimensional Richards Equation: Application to an Experimental Rain Garden

Rain gardens are infiltration systems that provide volume and water quality control, recharge enhancement, as well as landscape, ecological, and economic benefits. A model for application to rain gardens based on Richards equation coupled to a surface water balance was developed, using a two-dimensional finite-volume code. It allows for alternating upper boundary conditions, including ponding and overflow, and can simulate heterogeneous soil-layering or more complex geometries to estimate infiltration and recharge. The algorithm is conservative, and exhibits good performance compared to standard models for several test cases (less than 0.1% absolute mass balance error); simulations were also performed for an experimental rain garden and comparisons to collected data are presented. The model accurately simulated the matrix flow, soil water distribution, as well as deep percolation (potential recharge) for a natural rainfall event in the controlled experimental setup.     

Infiltration of Water into Soil with Cracks

This paper presents the physical basis of the FRACTURE submodel for simulating infiltration of precipitation∕irrigation water into relatively dry, cracked, fine-textured soils. The FRACTURE submodel forms part of the HYDRUS-ET variably saturated flow∕transport model. Infiltration into the soil matrix is formally divided into two components: (1) Vertical infiltration through the soil surface; and (2) lateral infiltration via soil cracks. The first component is described and solved using the 1D Richards' equation. Excess water that does not infiltrate through the soil surface is either considered to be runoff, if no soil cracks are present, or routed into soil cracks from where it may laterally infiltrate into the soil matrix. Horizontal infiltration from soil cracks into the soil matrix is calculated using the Green-Ampt approach and incorporated as a positive source∕sink term Sf in the Richards' equation describing flow in the matrix. In addition to the hydraulic properties of the soil matrix, the FRACTURE submodel requires parameters characterizing the soil cracks, notably the specific crack length per surface area lc and the relationship between crack porosity Pc and the gravimetric soil water content w. An example problem shows that infiltration from soil cracks can be an important process affecting the soil water regime of cracked soils. A comparison with the more traditional approach, involving surface infiltration only, indicates important differences in the soil water content distribution during a rainfall∕irrigation event. This extension of the classical approach to include crack infiltration significantly improves the identification and prediction of the soil water regime.     

Nitrogen Transformations during Soil–Aquifer Treatment of Wastewater Effluent—Oxygen Effects in Field Studies

Depth-dependent oxygen concentrations and aqueous-phase total ammonia and nitrate/nitrite ion concentrations were measured in the field during the infiltration of wastewater effluent. Measurements illustrated the dependence of nitrogen fate and transport on oxygen availability. Infiltration basins were operated by alternating wet (infiltration) and dry periods. During infiltration periods, ammonia was removed within the top few feet of sediments via adsorption. Biochemical activity rapidly eliminated residual molecular oxygen in the infiltrate, making the soil profile anoxic. During dry periods, oxygen reentered the basin profile and sorbed ammonia was converted to nitrate via nitrification. Oxygen penetrated to a depth of about 0.6?m?(2?ft) within the first few days of dry periods. At greater depths, oxygen levels increased more slowly due to a combination of slow transport kinetics and biochemical (nitrogenous) oxygen demand. During normal wet/dry basin cycles consisting of about 4 wet and 4 dry days, the local vadose zone remained anoxic at depths greater than about 1.5?m?(5?ft) below land surface. As a consequence, conditions for denitrification were satisfied in the deeper sediments. That is, the nitrate nitrogen produced in near surface sediments moved freely downward with infiltrating water where it encountered an extensive anoxic zone before reaching local monitoring or extraction wells. The relative importance of dissolved organics and sorbed ammonia as electron donors for denitrification reactions remains to be established.     

Water Table Rise in Sloping Aquifer due to Canal Seepage and Constant Recharge

Analytical solutions of the linearized Boussinesq equation and a fully implicit finite-difference numerical solution of the nonlinear Boussinesq equation were obtained to study transient and steady-state water table rise in a homogeneous, isotropic, and incompressible unconfined sloping aquifer. The rise was due to seepage from two canals located at different elevations above the sloping impermeable barrier and constant recharge from the land surface. Proposed analytical solutions were verified with existing analytical solutions for a horizontal aquifer and were found in close agreement. The effect of recharge and slope of the impermeable barrier on water table rise predicted by both analytical and numerical solutions was studied by considering a numerical example. The effect of the linearization of the Boussinesq equation on water table rise was also studied by comparing the water table heights predicted by the numerical solution with those computed from the analytical solution. The analytical solution overestimates water table elevations compared to those obtained from the numerical solution, and the difference in water table in the middle region decreases with increase in time.     

Controlled Field Experiment for Performance Evaluation of Septic Tank Effluent Treatment during Soil Infiltration

Decentralized systems are responsible for treating approximately 25% of the wastewater generated in the United States. The most common decentralized system involves onsite treatment using a septic tank unit followed by dispersal to a subsurface soil infiltration unit where percolation to groundwater occurs. To evaluate the hydraulic and purification processes occurring during soil treatment of septic tank effluent (STE), a field experiment was initiated in the Spring of 2003 with continued operation and monitoring for 2 years. A replicated factorial design (22) was employed to evaluate three infiltrative surface architectures (ISAs) (open, stone, and synthetic) and two daily hydraulic loading rates (HLRs) (4 and 8?cm/day). Pilot-scale test cells were established in native sandy loam soils at the Mines Park Test Site located on the Colorado School of Mines campus in Golden, Colo. STE was obtained from a nearby multifamily apartment building and applied to the test cells daily. Field monitoring included baseline characterization of soil and site properties, routine characterization of the STE applied, observations of STE ponding on the infiltrative surface, periodic measurement of constant-head infiltration rates, and periodic sampling and analyses of the soil pore water at 60- or 120-cm depths below the infiltrative surface. Monitoring revealed that the ISA and HLR influenced the rate and extent of hydraulic capacity loss during soil treatment. For example, an open horizontal infiltrative surface maintained an infiltration capacity that was 40–80% higher than one covered with either washed stones or synthetic aggregate. Purification of STE during infiltration and percolation through the sandy loam soil was very high. The cumulative mass removed during 2 years of operation for dissolved organic carbon, total nitrogen, and total phosphorus averaged 94, 42, and 99%, respectively. While there was no significant difference in the purification performance based on ISA or HLR, an increase in the vadose zone depth slightly increased purification.     

Storage Volume and Overflow Risk for Infiltration Basin Design

This study presents a risk-based approach for the design of infiltration basins. The design parameters include basin storage volume, drain time, and overflow risk. At a basin site, the storm-water detention storage volume is determined by design runoff capture volume, tributary watershed area, and runoff coefficient. The basin geometry is dictated by the water volume balance between the surface storage volume in the basin and the subsurface storage capacity in soil pores. The drain time at a basin site is found to be a function of initial soil water content, soil porosity, infiltration rate, and distance to the ground-water table. After knowing the basin geometry and size, the basin's performance can be evaluated by the overflow risk analysis using the local average event rainfall depth and interevent time. In practice, a sensitivity test on overflow risk can be conducted with a range of basin storage volumes. The risk-based approach presented in this study provides an algorithm to calculate the long-term runoff capture percentage for a basin size. The diminishing return on runoff capture percentage can serve as a basis to select the proper basin storage volume at the site.     

Integrated Water Management for the 21st Century: Problems and Solutions

Most of the projected global population increases will take place in third world countries that already suffer from water, food, and health problems. Increasingly, the various water uses (municipal, industrial, and agricultural) must be coordinated with, and integrated into, the overall water management of the region. Sustainability, public health, environmental protection, and economics are key factors. More storage of water behind dams and especially in aquifers via artificial recharge is necessary to save water in times of water surplus for use in times of water shortage. Municipal wastewater can be an important water resource but its use must be carefully planned and regulated to prevent adverse health effects and, in the case of irrigation, undue contamination of groundwater. While almost all liquid fresh water of the planet occurs underground as groundwater, its long-term suitability as a source of water is threatened by nonpoint source pollution from agriculture and other sources and by aquifer depletion due to groundwater withdrawals in excess of groundwater recharge. In irrigated areas, groundwater levels may have to be controlled with drainage or pumped well systems to prevent waterlogging and salinization of soil. Salty drainage waters must then be handled in an ecologically responsible way. Water short countries can save water by importing most of their food and electric power from other countries with more water, so that in essence they also get the water that was necessary to produce these commodities and, hence, is virtually embedded in the commodities. This “virtual” water tends to be a lot cheaper for the receiving country than developing its own water resources. Local water can then be used for purposes with higher social, ecological, or economic returns or saved for the future. Climate changes in response to global warming caused by carbon dioxide emission are difficult to predict in space and time. Resulting uncertainties require flexible and integrated water management to handle water surpluses, water shortages, and weather extremes. Long-term storage behind dams and in aquifers may be required. Rising sea levels will present problems in coastal areas.     

Integrated Hydrologic Modeling and GIS in Water Resources Management

The integration of a physically-based distributed model with a geographic information system (GIS) in watershed-based water resources management is presented, and an example watershed is chosen to demonstrate the spatial database and modeling system developed in this study. The spatial data is first processed by GIS. The model is then used to simulate runoff hydrographs. It operates at a daily time step on 1 × 1 km grid squares and simulates important hydrologic processes including evapotranspiration, snowmelt, infiltration, aquifer recharge, ground-water flow, and overland and channel runoff. Finally, the model result is displayed by using GIS. This study demonstrates that the integration of a physically based distributed model and GIS may successfully and efficiently implement the watershed-based water resources management. Not only does this process facilitate examination of a wider range of alternatives that would be impossible by using conventional methods, but it also provides a living management that could be modified and updated by water managers once the watershed condition is changed.     

Biotreatment of Surfactant Flush Wastewater from Wood Treatment Soil

The high cost of surfactant enhanced aquifer remediation (SEAR) could be reduced if contaminants in the extracted surfactant flushing water could be biodegraded, and the surfactant reinjected. To test this concept, ex situ biotreatment using immobilized bacteria was simulated in laboratory columns. Nonionic surfactant solutions of 1,000–5,000?mg/L Tergitol NP-10 (TNP10) were flushed through contaminated soil collected from a wood-treatment site. The highest TNP10 dose increased the effluent concentrations of tetrachlorophenol (TeCP) and tetrachlorohydroquinone (TCHQ) by about 3 and 16 times, respectively. These wash solutions were then treated by the aerobic bacterium Sphingomonas chlorophenolica RA2 immobilized in polyurethane foam. The immobilized bacteria were capable of degrading pentachlorophenol, TeCP and TCHQ in soil flushing wastewater containing up to 4,900?mg/L TNP10. Surfactant sorption to the biotreatment columns occurred, but these losses decreased over time as the sorption capacity of the foam was exhausted. The results suggest that SEAR wastewater could be biotreated, thus enabling reinjection of the surfactant.     

Analytical Model of Surge Flow in Nonprismatic Permeable Channels and Its Application in Arid Regions

A surge running down a dry wadi bed as a consequence of a controlled water release from a reservoir—e.g., for artificial groundwater recharge—represents a free boundary problem. After some time, when aiming for groundwater recharge, the infiltration equals inflow and thus forms a kind of “standing” wave. The numerical solution of such phenomena generally involves considerable problems. For avoiding the numerical inconvenience resulting from the complex interacting surface/subsurface flow, we present an analytical solution of the slightly modified zero-inertia (ZI) equations. The development introduces a momentum-representative cross section for portraying the transient development of momentum and refers to a channel with constant slope, irregular geometry, and a permeable channel bed with significant infiltration. Due to the structure of the solution, any arbitrary infiltration model can be used for quantifying the infiltration losses. For both synthetic prismatic and nonprismatic test channels, the robust and easy-to-use analytical ZI model shows an excellent match with the results of a comparative numerical simulation. Finally, the ZI model is employed for simulating a surge flow downstream of the Wadi Ahin groundwater recharge dam (Oman), in order to perform a scenario for artificial groundwater recharge in a natural wadi channel reach. This realistic application illustrates the potential of the new approach by even computing an almost standing wave and shows its applicability for an accurate and robust evaluation of release strategies.     

Temperature Effects on the Infiltration Rate through an Infiltration Basin BMP

Infiltration Best Management Practices (BMPs) are becoming more readily acceptable as a means of reducing postdevelopment runoff volumes and peak flow rates to pre-construction levels, while simultaneously increasing recharge. However, the design, construction, and operation of infiltration basins to this point have not been standardized due to a lack of understanding of the infiltration processes that occur in these structures. Sizing infiltration BMPs to hold and store a predetermined volume of runoff, typically called the Water Quality Volume, has become a widely accepted practice. This method of sizing BMPs does not account for the infiltration that is occurring in the BMP during the storm event; which could result in significantly oversized BMPs. The objective of this study was to develop a methodology to simulate varying infiltration rates observed from a large scale rock infiltration basin BMP. The results should aid in improved design of such structures. This methodology is required to predict the performance of these sites using single event and continuous flow models. The study site is a Pervious Concrete Infiltration Basin BMP built in 2002 in a common area at Villanova University. The system consists of three infiltration beds filled with coarse aggregate, lined with geotextile filter fabric, overlain with pervious concrete and underlain by undisturbed silty sand. The BMP is extensively instrumented to facilitate water quantity and quality research. The infiltration performance of the site is the focus of the study. Recorded data indicates a wide variation of linear infiltration rates for smaller storm events. A model was developed using the Green–Ampt formula to characterize the infiltration occurring in the basin for small storm events characterized by an accumulated depth of water in the infiltration bed of less than 10?cm. The effectiveness and accuracy of the model were measured by comparing the model outputs with observed bed water elevation data recorded from instrumentation on site. Results show that for bed depths of <10?cm, hydraulic conductivity is the most sensitive parameter, and that the storm event measured infiltration rate is substantially less then the measured saturated hydraulic conductivity of the soil. The governing factor affecting hydraulic conductivity, and subsequently, infiltration rate is temperature; with higher rates occurring during warmer periods, affecting the infiltration rate by as much as 56%.     

Agroecological Impacts from Salinization and Waterlogging in an Irrigated River Valley

Extensive field data and calibrated flow and salt-transport models characterize the spatial and temporal patterns of salinity and waterlogging in an irrigated western river valley. Over three irrigation seasons, average seasonal aquifer recharge from irrigated fields in a 50,600?ha study area ranges from 0.59?to?0.99?m, including contribution from precipitation. The salinity of irrigation water varies from 618?to?1,090?mg/L. The water table is shallow, with 16 to 33% of irrigated land underlaid by an average water table less than 2?m deep. Average water table salinity ranges from 2,680?to?3,015?mg/L, and average soil salinity from 2,490?to?3,860?mg/L. Crop yield reductions from salinity and waterlogging range from 0 to 89% on fields, with regional averages ranging from 11 to 19%. Annual salt loading to the river from subsurface return flows, generated in large part by dissolution from irrigation recharge, averages about 533?kg/irrigated?ha?per?km. Upflux from shallow water tables under fallow ground contributes to about 65?million?m3 (52,600?acre-ft) per year of nonbeneficial consumption. Beyond problem identification, the developed database and models provide a basis for effectively addressing these problems through a systematic and comparative assessment of alternative solutions.     

Development of Simplified Solutions for Modeling Recession in Basins

In irrigation basins the decrease in the gradient of the water-surface elevation following inflow cutoff often leads to reduced rate of convergence, increased computational time, and reduced robustness of the numerical solutions of the recession phase. As the water surface levels off, the underlying physical problem simplifies, thus allowing the use of highly accurate yet simple alternate solutions to the full-numerical solution of the zero-inertia equations. For level basins, the simplification involves treating the stream as a static pool, in which water level only falls in response to infiltration. Graded basins may require partitioning the stream into a flowing and static pool, before water-surface eventually levels off over the entire stream length. Implementation of these solutions enhances computational efficiency and robustness of surface irrigation models without a concomitant loss of accuracy. This paper discusses numerical problems related to the recession phase computation in basins and proposes simplified and robust, yet highly accurate solutions. A comparison of the recession trajectories and final infiltration profiles predicted by the full-numerical solution of the zero-inertia equations, obtained by using double-precision floating-point arithmetic, and the simplified alternate solutions, which is robust enough to be implemented over a range of hardware–software capabilities, show that the two approaches yield essentially identical results. Finally, the general validity of the proposed solutions is tested by comparing predictions of recession trajectories and infiltration profiles with those obtained using a surface irrigation hydraulic model, SRFR.     

淮南顾桂地区某矿区水文地质特征研究

水文地质调查是矿产资源开发与利用不可缺少的环节,对井下安全生产至关重要。本文以淮南顾桂地区某矿山为例,分析了矿区水文地质条件。研究表明:新生界松散含水层距离矿层较远,对矿床充水影响较小;二叠系砂岩裂隙含水层和石灰岩岩溶裂隙含水层受裂隙、岩溶发育程度而导致富水性能不均匀,补给水源较少,出水时间较短,容易疏干,常见的井下出水形式表现为淋滴水。