Vertical Wellbore Flow Monitoring for Assessing Spatial and Temporal Flow Relationships with a Dynamic River Boundary

A useful tool for identifying the temporal and spatial ambient wellbore flow relationships near a dynamic river boundary is to monitor ambient vertical wellbore flow with an electromagnetic borehole flowmeter. This is important because the presence of the wellbore can result in significant mixing or exchange of groundwater vertically across the aquifer. Mixing or exchanging groundwater within the well-screen section can have significant impacts on the distribution of contaminants within the aquifer and adverse effects on the representativeness of groundwater samples collected from the monitoring well. Ambient monitoring data, collected from long screened wells at Hanford’s 300-Area Integrated Field Research Challenge site, located approximately 260 m from the Columbia River, demonstrate that vertical wellbore flow exhibits both a positive and inverse temporal relationship with periodic river-stage fluctuations that can change over short distances between wells. The spatial distribution of these vertical flows across the well field indicates two general regions of ambient wellbore flow behavior. The western region of the site is characterized by vertical flows that are positively related to river-stage fluctuations. In contrast, the eastern region of the site exhibits vertical flows that are inversely related to river-stage fluctuations. The cause of this opposite relationship is not completely understood; however, the positive relationships appear to be associated with high-energy Hanford formation flood deposits. These flood deposits have a well-defined northwest-southeast trend and are believed to coincide with a local paleochannel. The inverse relationships are attributed to an erosional, subsurface high in the Hanford/Ringold Formation contact between the site and the Columbia River. Under these complex hydrogeologic and hydrodynamic conditions, the behavior of ambient vertical wellbore flow in monitoring wells near a dynamic river boundary can have important implications for collecting groundwater-quality samples, for contributing to contaminant distribution within an aquifer system, and for implementing effective remediation strategies.     

Concentration Profiles in Screened Wells under Static and Pumped Conditions

Purge and pump samples from screened wells reflect concentration averaging and contaminant redistribution by wellbore flow. These issues were assessed in a screened well at the Hanford Site by investigating the vertical profile of a technetium-99 plume in a conventional well under static and pumped conditions. Specific conductance and technetium-99 concentrations were well correlated, and this enabled measurement of specific conductance to be used as a surrogate for technetium-99 concentration. Time-series measurements were collected during purging from three specific conductance probes installed in the well at 1.2, 3.1, and 4.9 m below the static water level in a 7.7-m-deep screened well. The vertical contaminant profile adjacent to the well in the aquifer was calculated using the concentration profile in the well during pumping, the pumping flow rate, and a wellbore flow and mixing model. The plume was found to be stratified in the aquifer—the highest concentrations occurred adjacent to the upper part of the screened interval. The purge and pump sample concentrations were 41% to 58% of the calculated peak concentration in the aquifer. Plume stratification in the aquifer adjacent to the well screen became more pronounced as pumping continued. Extended pumping may have partially reversed the effect of contaminant redistribution in the aquifer by wellbore flow and allowed the stratification of the plume to be more observable. It was also found that the vertical profile of contamination in the well under static (i.e., nonpumping conditions) was not representative of the profile in the aquifer. Thus, passive or micropurge sampling techniques, which sample the wellbore water at different depths, would not yield results representative of the aquifer in this well.     

Removal of river-stage fluctuations from well response using multiple regression

Many contaminated unconfined aquifers are located in proximity to river systems. In groundwater studies, the physical presence of a river is commonly represented as a transient-head boundary that imposes hydrologic responses within the intersected unconfined aquifer. The periodic fluctuation of river-stage height at the boundary produces associated responses within the adjacent aquifer system, the magnitude of which is a function of the existing well, aquifer, boundary conditions, and characteristics of river-stage fluctuations. The presence of well responses induced by the river stage can significantly limit characterization and monitoring of remedial activities within the stress-impacted area. This article demonstrates the use of a time-domain, multiple-regression, convolution (superposition) method to develop well/aquifer river response function (RRF) relationships. Following RRF development, a multiple-regression deconvolution correction approach can be applied to remove river-stage effects from well water-level responses. Corrected well responses can then be analyzed to improve local aquifer characterization activities in support of optimizing remedial actions, assessing the area-of-influence of remediation activities, and determining mean groundwater flow and contaminant flux to the river system.     

Atrazine in a Stream-Aquifer System: Estimation of Aquifer Properties from Atrazine Concentration Profiles

Lincolns municipal wellfield consists of 44 wells developed in an alluvial aquifer adjacent to the Platte River near Ashland, Nebraska Induced recharge from the river is the primary source of water for the wellfield. Wafer samples were collected on a periodic basis from the Platte River arid two transects of monitoring wells. These samples were analyzed for the herbicide atrazine, which was used as a tracer of induced recharge in this stream-aquifer system. Atrazine concentrations in the river and aquifer were much less than 1.0 ppb during late fall and winter, but increased to as high as 18.9 ppb during spring and summer, associated with runoff from upgradient agricultural lands. There was approximately a 21-day lag time from the first detection of increasing atrazine concentration in the river to the first detection in monitoring wells immediately adjacent to the river. This lag time was relatively constant throughout the year and from one year to the next, even with major fluctuations of river stage and wellfield production. This consistency of lag time indicated that the travel times from the river to the first set of monitoring wells immediately adjacent to the river were fairly constant.
Paths of preferential flow were identified in the aquifer at a depth of 25 to 35 feet below land surface. This aquifer zone appeared to play a significant role in movement of water from beneath the river into the wellfield.
Aquifer dispersivity was calculated using a method described by Hoehn and Santschi (1987). Macrodispersivity (A_L) was shown to increase linearly over the scale of the wellfield. Calculated values of A_L were within limits of other reported values for this type of aquifer material and agreed well with values reported by Hoehn and Santschi (1987); These findings will be extremely beneficial for planning and management of the municipal wellfield.     

Effect of rapidly changing river stage on uranium flux through the hyporheic zone

Measurement of ground water/surface water interaction within the hyporheic zone is increasingly recognized as an important aspect of subsurface contaminant fate and transport. Understanding the interaction between ground water and surface water is critical in developing a complete conceptual model of contaminant transport through the hyporheic zone. At the Hanford Site near Richland, Washington, ground water contaminated with uranium discharges to the Columbia River through the hyporheic zone. Ground water flux varies according to changes in hydraulic gradient caused by fluctuating river stage, which changes in response to operation of dams on the Columbia River. Piezometers and continuous water quality monitoring probes were installed in the hyporheic zone to provide long-term, high-frequency measurement of hydraulic gradient and estimated uranium concentrations. Subsequently, the flux of water and uranium was calculated for each half-hour time period over a 15-month study period. In addition, measurement of water levels in the near-shore unconfined aquifer enhanced the understanding of the relationship between river stage, aquifer elevation, and uranium flux. Changing river stage resulted in fluctuating hydraulic gradient within the hyporheic zone. Further, influx of river water caused lower uranium concentrations as a result of dilution. The methods employed in this study provide a better understanding of the interaction between surface and ground water in a situation with a dynamically varying vertical hydraulic gradient and illustrate how the combination of relatively standard methods can be used to derive an accurate estimation of water and contaminant flux through the hyporheic zone.     

Innovative Contaminant Mass Flux Monitoring in an Aquifer Subject to Tidal Effects

Exposure from groundwater contamination to aquatic receptors residing in receiving surface water is dependent upon the rate of contaminated groundwater discharge. Characterization of groundwater fluxes is challenging, especially in coastal environments where tidal fluctuations result in transient groundwater flows towards these receptors. This can also be further complicated by the high spatial heterogeneity of subsurface deposits enhanced by anthropogenic influences such as the mixing of natural sediments and backfill materials, the presence of subsurface built structures such as sheet pile walls or even occurrence of other sources of contaminant discharge. In this study, the finite volume point dilution method (FVPDM) was successfully used to characterize highly transient groundwater flows and contaminant mass fluxes within a coastal groundwater flow system influenced by marked tides. FVPDM tests were undertaken continuously for more than 48 h at six groundwater monitoring wells, in order to evaluate groundwater flow dynamics during several tide cycles. Contaminant concentrations were measured simultaneously which allowed calculating contaminant mass fluxes. The study highlighted the importance of the aquifer heterogeneity, with groundwater fluxes ranging from 10⁻⁷ to 10⁻³ m/s. Groundwater flux monitoring enabled a significant refinement of the conceptual site model, including the fact that inversion of groundwater fluxes was not observed at high tide. Results indicated that contaminant mass fluxes were particularly higher at a specific monitoring well, by more than three orders of magnitude, than at other wells of the investigated aquifer. This study provided crucial information for optimizing further field investigations and risk mitigation measures.     

On‐Site Vapor‐Phase Analysis as a Novel Approach for Monitoring Groundwater Wells

The results of comprehensive field testing of on‐site vapor‐phase‐based groundwater monitoring methods are presented to demonstrate their utility as a robust and cost‐effective approach for rapidly obtaining volatile organic compounds (VOCs) concentration data from a monitoring well. These methods—which rely on sensitive, commercially available field equipment to analyze vapor in equilibrium with groundwater—proved easy to implement and can be tailored to site‐specific needs, including multilevel sampling. During field testing, low‐flow groundwater concentrations could be reasonably estimated using submerged passive vapor diffusion samplers or field equilibration of collected groundwater (R² = 0.85 to 0.96). These two methods are not as reliant on in‐well mixing to overcome vertical stratification within wells as simpler headspace methods. The importance of well and aquifer‐specific factors on concentration data (and therefore method selection) is highlighted, including the effect of changing in‐well patterns due to seasonal temperature gradients. Results indicated that vertical stratification was relatively limited within the set of wells included in these studies, resulting in similar performance for short depth‐discrete passive vapor diffusion samplers (constructed from 40‐mL vials) and longer samplers (2.5 to 5 feet in length) designed to cover a larger portion of the screened interval. A year‐long, multi‐event evaluation demonstrated that vapor‐phase‐based monitoring methods are no more variable than conventional groundwater monitoring methods, with both types subject to similar spatial and temporal variability that can be difficult to reduce. Vapor sampling methods represent a promising approach for estimation of groundwater concentrations by reducing the cost liabilities associated with monitoring while providing a more sustainable approach.     

Low-Flow Purging and Sampling of Ground Water Monitoring Wells with Dedicated Systems

A field study was conducted to assess purging requirements for dedicated sampling systems in conventional monitoring wells and for pumps encased in short screens and buried within a shallow sandy aquifer. Low-flow purging methods were used, and wells were purged until water quality indicator parameters (dissolved oxygen, specific conductance, turbidity) and contaminant concentrations (chromate, trichloroethylene, dichloroethylene) reached equilibrium. Eight wells, varying in depth from 4.6 to 15.2 m below ground surface, were studied. The data show that purge volumes were independent of well depth or casing volumes. Contaminant concentrations equilibrated with less than 7.5 I. of purge volume in all wells. Initial contaminant concentration values were generally within 20 percent of final values. Water quality parameters equilibrated in less than 10 L in all wells and were conservative measures for indicating the presence of adjacent formation water. Water quality parameters equilibrated faster in dedicated sampling systems than in portable systems and initial turbidity levels were lower.     

Coupled Aquifer-Borehole Simulation

A model coupling fluid hydraulics in a borehole with fluid flow in an aquifer is developed in this paper. Conservation of momentum is used to create a one-dimensional steady-state model of vertical flow in an open borehole combined with radially symmetric flow in an aquifer and with inflow to the well through the wellbore screen. Both laminar and turbulent wellbore conditions are treated. The influence of inflow through the wellbore screen on vertical flow in the wellbore is included, using a relation developed by Siwoń (1987) . The influence of inflow reduces the predicted vertical variation in head up to 15% compared to a calculation of head losses due to fluid acceleration and the conventional Colebrook-White formulation of friction losses in a circular pipe. The wellbore flow model is embedded into the MODFLOW-2000 ground water flow code. The nonlinear conservation of momentum equations are iteratively linearized to calculate the conductance terms for vertical flow in the wellbore. The resulting simulations agree favorably with previously published results when the model is adjusted to meet the assumptions of the previous coupled models.     

Measuring Flow Rate in Crystalline Bedrock Wells Using the Dissolved Oxygen Alteration Method

Determination of vertical flow rates in a fractured bedrock well can aid in planning and implementing hydraulic tests, water quality sampling, and improving interpretations of water quality data. Although flowmeters are highly accurate in flow rate measurement, the high cost and logistics may be limiting. In this study the dissolved oxygen alteration method (DOAM) is expanded upon as a low‐cost alternative to determine vertical flow rates in crystalline bedrock wells. The method entails altering the dissolved oxygen content in the wellbore through bubbler aeration, and monitoring the vertical advective movement of the dissolved oxygen over time. Measurements were taken for upward and downward flows, and under ambient and pumping conditions. Vertical flow rates from 0.06 to 2.30 Lpm were measured. To validate the method, flow rates determined with the DOAM were compared to pump discharge rates and found to be in agreement within 2.5%.     

Aerobic respiration in riparian exchange zones of regulated river corridors

River stage fluctuations drive surface water-groundwater exchanges within river corridors. This study evaluates how repeated daily stage fluctuations, representative of hydropeaking conditions, influence aerobic respiration of river-sourced dissolved organic carbon (DOC) in the riparian exchange zone using reactive flow and transport simulations. Over 50 hypothetical scenarios were modelled to evaluate how the duration of the daily flood signal, river DOC concentration, aquifer hydraulic conductivity and ambient groundwater flow condition affect the fate and transport of DOC and DO in the riparian aquifer. Time series subsurface snapshots highlight how the various factors influence the subsurface distribution of DOC and DO. The total mass of DOC respired per meter of river had a wide range depending on the parameters, spanning from 1.4 to 71 g over 24-h, with high hydraulic conductivity and losing ambient groundwater flow conditions favouring the largest amount of DOC respired. The ratio of DOC mass entering the riparian zone with the mass returning to the river showed that as little as 5% to as much as 76% of the DOC that enters the bank during stage fluctuations returns to the river. This return ratio is dependent on river DOC concentration, hydraulic conductivity and ambient groundwater flow. The results illustrate that stage variations due to river regulation can be a significant control on aerobic respiration in riparian exchange zones.     

Surface water-groundwater exchange dynamics in buried-valley aquifer systems

Groundwater is a primary source of drinking water worldwide, but excess nutrients and emerging contaminants could compromise groundwater quality and limit its usage as a drinking water source. As such contaminants become increasingly prevalent in the biosphere, a fundamental understanding of their fate and transport in groundwater systems is necessary to implement successful remediation strategies. The dynamics of surface water-groundwater (hyporheic) exchange within a glacial, buried-valley aquifer system are examined in the context of their implications for the transport of nutrients and contaminants in riparian sediments. High conductivity facies act as preferential flow pathways which enhance nutrient and contaminant delivery, especially during storm events, but transport throughout the aquifer also depends on subsurface sedimentary architecture (e.g. interbedded high and low conductivity facies). Temperature and specific conductance measurements indicate extensive hyporheic mixing close to the river channel, but surface water influence was also observed far from the stream-aquifer interface. Measurements of river stage and hydraulic head indicate that significant flows during storms (i.e., hot moments) alter groundwater flow patterns, even between consecutive storm events, as riverbed conductivity and, more importantly, the hydraulic connectivity between the river and aquifer change. Given the similar mass transport characteristics among buried-valley aquifers, these findings are likely representative of glacial aquifer systems worldwide. Our results suggest that water resources management decisions based on average (base) flow conditions may inaccurately represent the system being evaluated, and could reduce the effectiveness of remediation strategies for nutrients and emerging contaminants.     

Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.     

Hydrodynamics of floodplain wetlands in a chalk catchment: The River Lambourn, UK

Variations in floodplain channel water levels and valley floor groundwater levels (measured in piezometers and boreholes) are examined at selected points along the course of the River Lambourn, a chalk river in southern England. A local alluvial gravel aquifer in the valley bottom is associated with numerous small wetlands that extend over much of the river's perennial profile. Variations in hydraulic gradient between local borehole levels and/or floodplain channel water levels are described for three sites in the seasonal section of the channel at Bockhampton, East Garston and West Shefford. The results indicate that observed groundwater levels are closely associated with flows from discrete springs at the margins of the channel and floodplain. However, as the floodplain widens and the alluvial gravel aquifer increases in size, the gravel aquifer accounts for a substantial down-valley component of groundwater flow with a diffuse vertical water flux. In the lower catchment, the exchange of flows between the gravel aquifer and the river enables some attenuation of floodplain water-table variability, providing a stable hydrological regime for valley-bottom wetlands. Catchment controls upon the local, valley-bottom, wetland regime are demonstrated with the application of a simple groundwater model developed using MODFLOW. The model is used to simulate groundwater discharge to the river in the upper and lower catchment, in addition to the water level regime at selected points in the valley bottom in the lower catchment. The results demonstrate the importance of taking catchment-scale water flow into account when managing isolated wetlands in a permeable catchment.     

The Utility of Multiple-Completion Monitoring Wells for Describing a Solvent Plume

At a study site in the midwestern United States, multiple-completion wells demonstrated that a vertical hydraulic gradient was responsible for the contamination pattern exhibited by chlorinated solvent plumes. The typical pattern consisted of little or no contamination in the upper portion of the aquifer with concentrations increasing with depth. When ground water contamination was discovered in an unexpected portion of the site, water level elevations and contaminant distribution data obtained from multiple-completion wells resulted in identification of the source location. The well eventually determined to be located in the source area displayed contaminant levels much higher in the upper zone of the aquifer — the opposite contamination pattern of other on-site wells. Such results indicated that the spill had occurred near this location and that solvent residing along the capillary fringe was continuing to contaminate the aquifer.     

Single Well Tracer Tests in Aquifer Characterization

When the purpose of aquifer testing is to yield data for modeling aqueous mass transport, pumping tests and gradient measurement can only partially satisfy characterization requirements. Effective porosity, ground water flow velocity, and the vertical distribution of hydraulic conductivity within the aquifer are left as unknowns. Single well tracer methods, when added to the testing program, can be used to estimate these parameters. A drift, and pumpback test yields porosity and velocity, and point-dilution testing yields depth-discrete hydraulic information, A single emplacement of tracer into a test well is sufficient to conduct both tests. The tracer tests are facilitated by a simple method for injecting and evenly distributing the tracer solution into a wellbore, and by new ion-selective electrode instrumentation, specifically designed for submersible service, for monitoring the concentration of tracers such as bromide.     

Estimation of recharge from floods in disconnected stream-aquifer systems

Stream-aquifer interaction has been the subject of much research for cases of good hydraulic connection (continuous saturated zone) between a river and an aquifer. Under these conditions, floods do not represent a very large net input to the aquifer because most of the water that enters the aquifer during the flood returns to the river when its stage recedes. The situation is different in disconnected stream-aquifer systems, where the streambed lies above the water level in the aquifer, thus preventing return flow from the aquifer. Under these conditions, floods may represent large, but hard to quantify, water inputs. Here, we present a methodology to estimate recharge from floods for disconnected stream-aquifer systems. Recharge is estimated as the product of a flood time function (dependent on the streamflow) and an unknown factor, which is obtained from calibrating a ground water flow model to aquifer heads. The approach can also benefit from concentration data, which can be very informative when river water concentrations vary over time. This methodology is applied to a field situation where recharge from river flooding is found to amount to nearly 15 million m(3)/year on the average, which represents 40% of the total aquifer inputs. Recharge from flooding helps explain major head recoveries, suggesting that basin water management programs should allow some floods to occur.     

Migration of Injected Wastewater with High Levels of Ammonia in a Saline Aquifer in South Florida

Treated wastewater with high levels of ammonia has been injected, since March 1983 into the deep saline units of the Lower Floridan aquifer (LFA) from a treatment plant near the east coast of Miami-Dade County in southeastern Florida. Monitoring wells in the plant recorded ammonia concentrations above ambient levels at hydrogeologic units located about 1000 ft (304.8 m) above injection depths between 2500 and 2800 ft (762 and 853 m) below sea level. A solute-transport model was developed to assess the horizontal and vertical extent of the injected ammonia, with ammonia moving from the injected zone into the overlying units: the upper semiconfining unit, the uppermost permeable zone of the LFA, and the middle semiconfining units of the Avon Park Formation. Ammonia is assumed to be transported under the effects of local heterogeneity in a porous limestone aquifer with high-salinity ambient groundwater and via upward migration through quasi-vertical pathways. A flow model of the migration of the injected ammonia was calibrated with PEST using head, salinity, and ammonia concentration data measured from 1983 to 2013. Borehole geophysical data support the high permeability of the uppermost permeable zone in the LFA. Average simulated head, normalized salinity, and ammonia concentration residuals over all monitoring wells were −1.37 ft, 0.01, and −0.67 mg/L, respectively. Model results are consistent with undetectable ammonia concentrations in the Upper Floridan aquifer.     

Patterns and Rates of Ground-Water Flow on Long Island, New York

Increased ground-water contamination from human activities on Long Island has prompted studies to define the pattern and rate of ground-water movement. A two-dimensional, fine-mesh, finite-element model consisting of 11,969 nodes and 22,880 elements was constructed to represent ground-water flow along a north-south section through central Long Island. The model represents average hydrologic conditions within a corridor approximately 15 miles wide. The model solves discrete approximations of both the potential and stream functions. The resulting flownet depicts flow paths and defines the vertical distribution of flow within the section. Ground-water flow rates decrease with depth. Sixty-two percent of the water flows no deeper than the upper glacial (water-table) aquifer, 38 percent enters the underlying Magothy aquifer, and only 3.1 percent enters the Lloyd aquifer. The limiting streamlines for flow to the Magothy and Lloyd aquifers indicate that aquifer recharge areas are narrow east-west bands through the center of the island. The recharge area of the Magothy aquifer is only 5.4 miles wide; that of the Lloyd aquifer is less than 0.5 miles. The distribution of ground-water traveltime and a flownet are calculated from model results; both are useful in the investigation of contaminant transport or the chemical evolution of ground water within the flow system. A major discontinuity in traveltime occurs across the streamline which separates the flow subsystems of the two confined aquifers. Water that reaches the Lloyd aquifer attains traveltimes as high as 10,000 years, whereas water that has not penetrated deeper than the Magothy aquifer attains traveltimes of only 2,000 years. The finite-element approach used in this study is particularly suited to ground-water systems that have complex hydrostratigraphy and cross-sectional symmetry.     

A Method for Conducting Simultaneous Convergent Tracer Tests in Multilayered Aquifers

Forced gradient tracer tests between two boreholes can be used to study contaminant transport processes at the small field scale or investigate the transport properties of an aquifer. Full depth tests, in which tracer samples are collected just from the discharge of the abstraction borehole, often give rise to breakthrough curves with multiple peaks that are usually attributed to different flow paths through the aquifer that can rarely be identified from the test results alone. Tests in selected levels of the aquifer, such as those between packer‐isolated sections of the boreholes, are time consuming, expensive; and the identification of major transport pathways is not guaranteed. We present a method for simultaneously conducting multiple tracer tests covering the full depth of the boreholes, in which tracer sampling and monitoring is carried out by a novel multilevel sampling system allowing high frequency and cumulative sampling options. The method is applied to a tracer test using fluorescein conducted in the multilayered sandstone aquifer beneath the city of Birmingham, UK, producing six well‐defined tracer breakthrough curves.