Estimating Travel Time in Bank Filtration Systems from a Numerical Model Based on DTS Measurements

An approach is presented to determine the seasonal variations in travel time in a bank filtration system using a passive heat tracer test. The temperature in the aquifer varies seasonally because of temperature variations of the infiltrating surface water and at the soil surface. Temperature was measured with distributed temperature sensing along fiber optic cables that were inserted vertically into the aquifer with direct push equipment. The approach was applied to a bank filtration system consisting of a sequence of alternating, elongated recharge basins and rows of recovery wells. A SEAWAT model was developed to simulate coupled flow and heat transport. The model of a two‐dimensional vertical cross section is able to simulate the temperature of the water at the well and the measured vertical temperature profiles reasonably well. MODPATH was used to compute flowpaths and the travel time distribution. At the study site, temporal variation of the pumping discharge was the dominant factor influencing the travel time distribution. For an equivalent system with a constant pumping rate, variations in the travel time distribution are caused by variations in the temperature‐dependent viscosity. As a result, travel times increase in the winter, when a larger fraction of the water travels through the warmer, lower part of the aquifer, and decrease in the summer, when the upper part of the aquifer is warmer.     

Relative Recovery of Thermal Energy and Fresh Water in Aquifer Storage and Recovery Systems

This paper explores the relationship between thermal energy and fresh water recoveries from an aquifer storage recovery (ASR) well in a brackish confined aquifer. It reveals the spatial and temporal distributions of temperature and conservative solutes between injected and recovered water. The evaluation is based on a review of processes affecting heat and solute transport in a homogeneous aquifer. In this simplified analysis, it is assumed that the aquifer is sufficiently anisotropic to inhibit density‐affected flow, flow is axisymmetric, and the analysis is limited to a single ASR cycle. Results show that the radial extent of fresh water at the end of injection is greater than that of the temperature change due to the heating or cooling of the geological matrix as well as the interstitial water. While solutes progress only marginally into low permeability aquitards by diffusion, conduction of heat into aquitards above and below is more substantial. Consequently, the heat recovery is less than the solute recovery when the volume of the recovered water is lower than the injection volume. When the full volume of injected water is recovered the temperature mixing ratio divided by the solute mixing ratio for recovered water ranges from 0.95 to 0.6 for ratios of maximum plume radius to aquifer thickness of 0.6 to 4.6. This work is intended to assist conceptual design for dual use of ASR for conjunctive storage of water and thermal energy to maximize the potential benefits.     

Geochemical processes during five years of aquifer storage recovery

A key factor in the long-term viability of aquifer storage recovery (ASR) is the extent of mineral solution interaction between two dissimilar water types and consequent impact on water quality and aquifer stability. We collected geochemical and isotopic data from three observation wells located 25, 65, and 325 m from an injection well at an experimental ASR site located in a karstic, confined carbonate aquifer in South Australia. The experiment involved five major injection cycles of a total of 2.5 x 10(5) m3 of storm water (total dissolved solids [TDS] approximately 150 mg/L) into the brackish (TDS approximately 2400 mg/L) aquifer. Approximately 60% of the mixture was pumped out during the fifth year of the experiment. The major effect on water quality within a 25 m radius of the injection well following injection of storm water was carbonate dissolution (35 +/- 6 g of CaCO3 dissolved/m3 of aquifer) and sulfide mineral oxidation (50 +/- 10 g as FeS2/m3 after one injection). < 0.005% of the total aquifer carbonate matrix was dissolved during each injection event, and approximately 0.2% of the total reduced sulfur. Increasing amounts of ambient ground water was entrained into the injected mixture during each of the storage periods. High 14C(DIC) activities and slightly more negative delta13C(DIC) values measured immediately after injection events show that substantial CO2(aq) is produced by oxidation of organic matter associated with injectant. There were no detectable geochemical reactions while pumping during the recovery phase in the fifth year of the experiment.     

Study of the feasibility of an aquifer storage and recovery system in a deep aquifer in Belgium

Abstract

The feasibility of aquifer storage and recovery (ASR) was tested in a deep aquifer near Koksijde, Belgium. To achieve this, oxic drinking water was injected into a deep aquifer (the Tienen Formation) that contains anoxic brackish water. The hydraulic properties of the aquifer were determined using a step-drawdown test. Chemical processes caused by the injection of the water were studied by two push—pull tests. The step-drawdown test was interpreted by means of an inverse numerical model, resulting in a transmissivity of 3.38 m²/d and a well loss coefficient of 0.00038 d²/m⁵. The push—pull tests identified mixing between the injection and pristine waters, and cation exchange, as the major processes determining the quality of the recovered water. Mobilization of DOC, aerobic respiration, denitrification and mobilization of phosphate were also observed.     

Isotopes and geochemistry in a managed aquifer recharge scheme: a case study of fresh water injection at the Damascus University Campus,Syria

A combination of stable isotopes (¹⁸O and ²H) and hydrochemistry has been applied to investigate storage processes in relation to aquifer storage and recovery (ASR) of the shallow alluvial Quaternary aquifer in Damascus basin. The stored water, entirely taken from the Figeh springs during flood periods, was injected in a single well having a brackish groundwater. Water samples were collected from four observation wells drilled in the Damascus University Campus (DUC) site during a 3‐year period (2006–2008). The injectant water, which deviates in its chemical and isotopic signatures from that of the ambient groundwater, shows that the stored water plume remains within close proximity to the injection well (IW) (<≈ 100 m). Thus, only two wells (W13 and W14) located at a distance less than 80 m from the injection point were affected by this injection. The observation wells located at longer distances from the IW (≈145 m and ≈ 600 m for wells W15 and WHz, respectively) were completely unaffected by the injection. Although most of the chemical and isotopic parameters usefully reflected the mixing process that occurs between the injectant water and ambient groundwater, the stable isotope (¹⁸O) and chloride (Cl⁻) were the most sensitive parameters that quickly reflect this signature. Using a simple mass balance, the calculated proportion of injectant water reaching the well W13 was in the range of 50–90%. This proportion was even lower (30–55%) in the case of well W14. Although the drought event prevailing during this study did not much help to inject further amounts of water, higher than the injected volume (0·2416 M m³) and also not favourable to better evaluate the fate and subsurface hydrological processes, these findings offer encouragement to continue the ASR activities, as an alternative way for better management of water resources in this basin facing intensive problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Aquifer Storage and Recovery Using Saline Aquifers: Hydrogeological Controls and Opportunities

Aquifer storage and recovery (ASR) is a valuable tool for managing variations in the supply and demand of freshwater, but system performance is highly dependent upon system-specific hydrogeological conditions including the salinity of the storage-zone native groundwater. ASR systems using storage zones containing saline (>10,000 mg/L of total dissolved solids) groundwater tend to have relatively low recovery efficiencies (REs). However, the drawbacks of low REs may be offset by lesser treatment requirements and may be of secondary importance where the stored water (e.g., excess reclaimed, surface, and storm waters) would otherwise go to waste and pose disposal costs. Density-dependent, solute-transport modeling results demonstrate that the RE of ASR systems using a saline storage zone is most strongly controlled by parameters controlling free convection (e.g., horizontal hydraulic conductivity) and mixing of recharged and native groundwater (e.g., dispersivity and aquifer heterogeneity). Preferred storage zone conditions are moderate hydraulic conductivities (5 to 20 m/d), low degrees of aquifer heterogeneity, and primary porosity-dominated siliclastic and limestones lithologies with effective porosities greater than 5%. Where hydrogeological conditions are less favorable, operational options are available to improve RE, such as preferential recovery from the top of the storage zone. Injection of large volumes of excess water currently not needed into saline aquifers could create valuable water resources that could be tapped in the future during times of greater need.     

Effect of multilayered groundwater mounds on water dynamics beneath a recharge basin: Numerical simulation and assessment of surface injection

Interactions between groundwater mounds caused by a geologic layer contrast affect the efficiency of managed aquifer recharge in arid areas. However, research has rarely examined the roles of groundwater mounding size variations on soil water dynamics in a stratified vadose zone in response to a sustained infiltration source. Numerical experiments were conducted on a two-dimensional vertical-section domain using HYDRUS software to simulate the behaviours of two adjacent (upper and lower) groundwater mounds underlying an infiltration basin subjected to clay loam and sandy alternately-layered soil profiles. The model successfully predicted the volume and extent of perched water and approximated vertical travel times during events generating downward fluxes from the surface injection. The response time of the mounding width (lateral extension) to the surface injection was delayed as compared to that of the mounding height (vertical extension), especially for the lower water mound. The mounding heights and widths show a strongly positive correlation with the infiltration rates of both high- and low-permeability layers where the injected water mounded, while the water storage amounts in the high- and low-permeability layers were governed by the mounding height and width, respectively. Exploratory simulations were then employed to assess the dependence of groundwater mounding behaviours and recharge performances on surface injection strategies. Results suggest that, by reducing injection rate or shortening injection duration, the near-term fraction of the surface injection converted to deep recharge is likely to be increased due to the narrowed groundwater mounding size, which would be limited by the water-retarding effect of layer contrasts. This study has important implications for predicting and understanding multilayered groundwater mounding behaviours and associated water mass balance under the geologic stratification, and is expected to aid in optimizing the infiltration basin operation for aquifer recharge.     

An assessment of aquifer storage recovery using ground water flow models

Owing to increased demands on ground water accompanied by increased drawdowns, technologies that use recharge options, such as aquifer storage recovery (ASR), are being used to optimize available water resources and reduce adverse effects of pumping. In this paper, three representative ground water flow models were created to assess the impact of hydrogeologic and operational parameters/factors on recovery efficiency of ASR systems. Flow/particle tracking and solute transport models were used to track the movement of water during injection, storage, and recovery. Results from particle tracking models consistently produced higher recovery efficiency than the solute transport models for the parameters/properties examined because the particle tracking models neglected mixing of the injected and ambient water. Mixing between injected and ambient water affected recovery efficiency. Results from this study demonstrate the interactions between hydrogeologic and operational parameters on predictions of recovery efficiency. These interactions are best simulated using coupled numerical ground water flow and transport models that include the effects of mixing of injected water and ambient ground water.     

Small‐Scale ASR Between Flow Barriers in a Saline Aquifer

Regular aquifer storage recovery, ASR, is often not feasible for small‐scale storage in brackish or saline aquifers because fresh water floats to the top of the aquifer where it is unrecoverable. Flow barriers that partially penetrate a brackish or saline aquifer prevent a stored volume of fresh water from expanding sideways, thus increasing the recovery efficiency. In this paper, the groundwater flow and mixing is studied during injection, storage, and recovery of fresh water in a brackish or saline aquifer in a flow‐tank experiment and by numerical modeling to investigate the effect of density difference, hydraulic conductivity, pumping rate, cyclic operation, and flow barrier settings. Two injection and recovery methods are investigated: constant flux and constant head. Fresh water recovery rates on the order of 65% in the first cycle climbing to as much as 90% in the following cycles were achievable for the studied configurations with constant flux whereas the recovery efficiency was somewhat lower for constant head. The spatial variation in flow velocity over the width of the storage zone influences the recovery efficiency, because it induces leakage of fresh water underneath the barriers during injection and upconing of salt water during recovery.     

Integrated assessment of lateral flow, density effects and dispersion in aquifer storage and recovery

Aquifer storage and recovery (ASR) involves the injection of freshwater into an aquifer for later recovery and use. This paper investigates three major factors leading to reduction in performance of ASR systems in brackish or saline aquifers: lateral flow, density-driven flow and dispersive mixing. Previous analyses of aquifer storage and recovery (ASR) have considered at most two of the above processes, but never all three together, and none have considered lateral flow and density effects together. In this analysis, four dimensionless parameters are defined to give an approximate characterisation of lateral flow, dispersive mixing, mixed convection (density effects during pumping) and free convection (density effects during storage). An extensive set of numerical models spanning a wide parameter range is then used to develop a predictive framework using the dimensionless numbers. If the sum of the four dimensionless numbers (denoted R_ASR) exceeds 10, the ASR operation is likely to fail with no recoverable freshwater, while if R_ASR < 0.1, the ASR operation is likely to provide at least some recovery of freshwater. The predictive framework is tested using limited data available from ASR field sites, broadly lending support to the framework. This study has several important implications. Firstly, the lack of completeness of field data sets in the literature must be rectified if we are to properly characterise mixed-convective flow processes in ASR operations. Once data are available, the dimensionless numbers can be used to identify suitable ASR sites and the desirable operational conditions that maximise recovery efficiencies.     

Evaluating the reliability of time series analysis to estimate variable riparian travel times by numerical groundwater modelling

The transition zones between rivers and adjacent riparian aquifers are locations of high biogeochemical activities that contribute to a removal of potentially hazardous substances in the aquatic system. The potential of the removal processes depends highly on subsurface water travel times, which can be determined by using the propagation of electrical conductivity (EC) signal from the river into the riparian aquifer. Although this method has been applied and verified in many studies, we observe possible limitations for the usage of EC fluctuation analysis. Our findings are based on EC time series analyses during storm events and artificial hydropeaks induced by watermill operations. Travel times derived by cross‐correlation analysis were compared with travel times calculated based on backward particle tracking of a calibrated transient numerical groundwater flow model. The cross‐correlation method produced only reasonable travel times for the artificial hydropeaks. In contrast, cross‐correlation analysis of the EC data during natural storm events resulted in implausibly negative or unrealistically low travel times for the bulk of the data sets. We conclude that the reason for this behaviour is, first, the low EC contrast between river and groundwater in connection with a strong damping of the infiltrating river EC signal into the subsurface during storm events. Second, the existence of old and less‐mineralized riparian water between the river and the monitoring well resulted in bank‐storage‐driven EC breakthrough curves with earlier arrival times and the subsequent estimation of implausible riparian travel times.     

Dampening of fluctuations in groundwater temperature by heat exchange between the aquifer and the adjacent layers

Heat storage in aquifers is attractive from the point of view of energy conservation. However, the influence on the chemical and microbiological composition of groundwater is insufficiently known. The nature and the extent of this influence can best be studied in the field. An understanding of the transport of heat in the aquifer will be indispensable for the planning of field experiments. For this purpose several aspects of non-steady and periodic heat transport are discussed, with special attention to the heat exchange between an aquifer and its adjacent layers. As in the case of steady heat transport, the changes of the temperature around an injection well can be expressed as a function of the retention time of the injected water. Hence, the heat flow pattern can be found by the calculation of groundwater streamlines and retention times.     

Thermal Impact of Residential Ground-Water Heat Pumps

A computer simulation study was conducted to quantify the potential thermal impact of residential water-source heat pump usage on ground-water aquifers. In a first phase of the study, weather data for nine locations throughout the country were used to estimate the energy requirements for heating and air conditioning a typical residence. These energy requirements were then translated into the volumetric water demands for a selected heat pump at each location. A representative model aquifer was then defined and its characteristics used, along with the heat pump water requirements and design ΔT's (difference between inlet and outlet water temperature) to identify the important parameters that contribute to heat transfer and to model the movement of the thermal front resulting from injection of heat pump discharge water at the nine locations. The major factor that determines the heat pump thermal impact was found to be the net amount of heat injected into, or removed from an aquifer. Other significant factors included well design, heat pump design ΔT, and physical properties of the aquifer such as thickness, porosity and dispersivity. The study showed that, in climates where winter heating demand is very nearly equal to summer cooling demands, the injection of heat pump discharge water did not cause any significant modification of the ambient model aquifer temperature. However, in hot or cold climates where air conditioning or heating demand dominates, measurable thermal changes occurred in the model aquifer. In most cases, the maximum temperature     

Facilitating Open Pit Mine Closure with Managed Aquifer Recharge

Dewatering of open pit mines can lower the regional water table for distances of several kilometers from the pit. When the mine is closed, dewatering operations usually cease, and the water table near the pit begins to rise. If the pit is backfilled, the water table will eventually recover, but this recovery may take several hundred years. However, if the extracted water is re-injected into the subsurface, then this may accelerate recovery of the water table. We show that there is an optimal distance for re-injection, which is sufficiently far from the mine to minimize the amount of groundwater that flows back to the pit during mine operations (and hence necessitate additional pumping) but is still close enough to speed up the water table recovery post-mine closure. The optimal injection distance increases with the aquifer hydraulic diffusivity and the mine life (duration of dewatering and injection), and typically ranges between about two and nine times the radius of the mine pit. Where the mine pit is not backfilled, the relative reduction in drawdown due to injecting all the pumped water at the optimal distance is between approximately 10% and 50% after a recovery time equal to the mining period, increasing to 30% to 90% after a recovery time five times the mining period. The relative drawdown reduction due to managed aquifer recharge will be even greater for a pit which is backfilled when mining ceases.     

Fluctuations of electrical conductivity as a natural tracer for bank filtration in a losing stream

A key parameter used in the assessment of bank filtration is the travel time of the infiltrated river water during the passage through groundwater. We analyze time series of electrical conductivity (EC) in the river and adjacent groundwater observation wells to investigate travel times of young hyporheic groundwater in adjoining channelized and restored sections of River Thur in North-East Switzerland. To quantify mixing ratios and mean residence times we perform cross-correlation analysis and non-parametric deconvolution of the EC time series. Measurements of radon-222 in the groundwater samples validate the calculated residence times. A simple relationship between travel time and distance to the river has not been observed. Therefore, we speculate that the lateral position and depth of the thalweg as well as the type of bank stabilization might control the infiltration processes in losing rivers. Diurnal oscillations of EC observed in the river and in nearby observation wells facilitate analyzing the temporal variation of infiltration. The diurnal oscillations are particularly pronounced in low flow situations, while the overall EC signal is dominated by individual high-flow events. Differences in travel times derived from diurnal and overall EC signals thus reflect different infiltration regimes.     

Geophysical and Hydrochemical Identification of Flow Paths with Implications for Water Quality at an ARR Site

Two key challenges regarding the design and operation of aquifer recharge and recovery (ARR) systems are evaluating aquifer heterogeneity and understanding hydrochemical interactions. Uncertainty in this respect can impact the volume of recoverable water and the improvement in water quality. The objective of this research is to leverage the advantages of geophysical measurements and hydrochemical sampling to reveal the properties of an ARR site to inform current ARR system operations and future design decisions. Electrical resistivity tomography was used to image the subsurface below two key infiltration/extraction areas of an ARR site in Colorado, USA. Hydrochemical measurements on transects intersecting the geophysical measurements resolved bulk parameters (i.e., total organic carbon, nitrate, and major cations and anions) and trace organic chemicals (e.g., pharmaceuticals, personal care products). Conservative tracers were also used to estimate degrees of mixing and water travel times and to better assess the performance of the ARR site regarding water quality changes and water recovery. The electrical resistivity measurements suggest that certain areas of the infiltration basins have hydraulic connections to the extraction wells through preferential flow paths, compared with other infiltration basins that are separated by fine‐grained materials from their respective extraction wells. The hydrochemical results indicate that consistent improvements in water quality can be achieved in these preferential flow paths within relatively short travel times (<5 d) at this ARR site.     

Optimizing Multiwell Aquifer Storage and Recovery Systems for Energy Use and Recovery Efficiency

As more aquifer storage and recovery (ASR) systems are employed for management of water resources, the skillful operation of multiwell ASR systems has become very important to improve their performance. In this study, we developed MODFLOW and MT3DMS models to simulate a multiwell ASR system in a synthetic aquifer to assess effects of hydrogeological and operational factors on the performance of the multiwell ASR system. We evaluated a simplified (dual well) ASR system in comparison with complex system (three-, four-, five-, and seven-well systems). Recovery and energy efficiencies were calculated using the model simulations. Factors such as higher hydraulic conductivity and longitudinal dispersivity significantly reduced the recovery and energy efficiencies of the system. In contrast, increasing the volume of recharged water increased the recovery efficiency; however, the energy efficiency was reduced. Recovery and energy efficiencies also plummet when there is an increase in the underlying regional gradient and the designed storage duration. Operating the system multiple times can yield higher volume of potable water, but the energy efficiency may not vary significantly after the second operating cycle. Single-well systems and multiwell systems exhibit similar responses to changes in physical factors, although operational factors have a more pronounced effect on the multiwell systems. One of the major findings was that fewer wells in a multiwell ASR system can yield higher volume of potable water and better output with respect to the electrical power being consumed. The results provide design engineers with guidelines for optimizing performance of the multiwell ASR systems.     

Subsurface Injection of Treated Sewage into a Saline-Water Aquifer at St. Petersburg, Florida — Water-Quality Changes and Potential for Recovery of Injected Sewage

The city of St. Petersburg is testing subsurface injection of treated sewage into the Floridan aquifer as a means of eliminating discharge of sewage to surface waters and as a means of storing treated sewage for future nonpotable reuse. The injection zone at the test site at the start of injection contained saline water with chloride concentrations ranging from 14,000 to 20,000 milligrams per liter (mg/1). Treated sewage with a mean chloride concentration of 170 mg/1 was injected through a single well for 12 months at a mean rate of 4.7 × 10⁵ cubic feet per day. The volume of water injected during the year was 1.7 × 10⁸ cubic feet. Dissolved oxygen was contained in the sewage prior to injection. Water removed from the injection zone during injection was essentially free of oxygen. Probable growth of denitrifying bacteria and, thus, microbial denitri-fication, was suggested by bacterial counts in water from two observation wells that were close to the injection well. The volume fraction of treated sewage in water from wells located 35 feet and 733 feet from the injection well and open to the upper part of the injection zone stabilized at about 0.9 and 0.75, respectively. Chloride concentrations stabilized at about 1,900 mg/1 in water from the well that was 35 feet from the injection well and stabilized at about 4,000 mg/1 in water from the well that was 733 feet from the injection well. These and other data suggest that very little near injection-quality treated sewage would be recoverable from storage in the injection zone.     

Analysis of benzene transport in a two‐dimensional aquifer model

In this study, the fate and transport of aqueous benzene was investigated in a laboratory‐scale homogeneous aquifer by conducting a two‐dimensional plume test. Benzene solution was introduced as a pulse type along the width of the aquifer model through a recharge zone situated at the upper‐left part of the model and followed by a steady state flow. Solution samples were collected at various locations on the front side of the model to capture two‐dimensional plumes at discrete time intervals. The benzene plumes showed a moderate retardation relative to chloride plumes observed from the previous study conducted for the same aquifer model. The retardation factor was obtained from the ratio of travel distances of benzene peaks to chloride peaks from the injection point, computed using a line integral method. Mass recovery of aqueous benzene revealed that there was a significant reduction of benzene mass, indicating the occurrence of volatilization and/or irreversible sorption during transport. Thus, retardation along with volatilization and/or irreversible sorption may be important processes affecting the fate and transport of aqueous benzene in the aquifer model. Copyright © 2005 John Wiley & Sons, Ltd.     

Influence of variable salinity conditions in a tidal creek on riparian groundwater flow and salinity dynamics

While recent studies have revealed that tidal fluctuations in an estuary significantly affect groundwater flows and salt transport in the riparian zone, only seawater salinity in the estuary has been considered. A numerical study is conducted to investigate the influence of estuarine salinity variations on the groundwater flow and salt dynamics in the adjacent aquifer to extend our understanding of these complex and dynamic systems. Tidal salinity fluctuations (synchronous with estuary stage) were found to alter the magnitude and distribution of groundwater discharge to the estuary, which subsequently impacted on groundwater salinity patterns and residence times, especially in the riparian zone. The effects of salinity fluctuations were not fully captured by adopting a constant, time-averaged estuarine salinity. The modelling analysis also included an assessment of the impact of a seasonal freshwater flush in the estuary, similar to that expected in tropical climates (e.g. mean estuary level during flood significantly greater than average), on adjacent groundwater flow and salinity conditions. The three-month freshwater flushing event temporarily disrupted the salt distribution and re-circulation patterns predicted to occur under conditions of constant salinity and tidal water level fluctuations in the estuary. The results indicate that the salinity variations in tidal estuaries impact significantly on estuary–aquifer interaction and need to be accounted for to properly assess salinity and flow dynamics and groundwater residence times of riparian zones.