Infiltration of river water to a shallow aquifer investigated with H/He, noble gases and CFCs

Noble gas isotopes (³He, ⁴He, Ne, Ar, Kr, Xe), tritium (³H), chlorofluorocarbons (CFCs) and dissolved oxygen (O₂) were seasonally measured in a small groundwater system recharged by infiltration of river water at Linsental, northeastern Switzerland. All Groundwater samples contained an excess of atmospheric noble gases (‘excess air’) usually with an elemental composition equal to air. The concentrations of atmospheric noble gases in the groundwater were used to calculate the excess air component and the water temperature at recharge. The noble gas temperatures (NGTs) in the boreholes close to the river vary seasonally, however, the average NGT of all samples lies close to the mean annual temperature of the river water. Groundwater ages were calculated using the tritium/helium-3 (³H/³He) dating method. The water ages of the samples obtained near the river depend on the amount of recently infiltrated river water and are young during times of active river discharge. In contrast, the mean water age of about 3 years of the deep aquifer remained nearly constant over the sampling period. The observed CFC-11 (CFCl₃) and CFC-12 (CF₂Cl₂) concentrations are significantly higher than the atmospheric equilibrium concentrations and therefore CFCs do not provide any direct information on the residence time of the groundwater. Nevertheless, the CFC excess in the groundwater shows a linear increase with the ³H/³He age. Additionally, both accumulation of radiogenic He (⁴He_rad) and O₂ consumption are strongly correlated with residence time. All these correlations can be interpreted either in terms of mixing of recently infiltrated river water with older groundwater or in terms of accumulation/consumption rates.     

Ground water flow analysis of a mid-Atlantic outer coastal plain watershed, Virginia, U.S.A

Models for ground water flow (MODFLOW) and particle tracking (MODPATH) were used to determine ground water flow patterns, principal ground water discharge and recharge zones, and estimates of ground water travel times in an unconfined ground water system of an outer coastal plain watershed on the Delmarva Peninsula, Virginia. By coupling recharge and discharge zones within the watershed, flowpath analysis can provide a method to locate and implement specific management strategies within a watershed to reduce ground water nitrogen loading to surface water. A monitoring well network was installed in Eyreville Creek watershed, a first-order creek, to determine hydraulic conductivities and spatial and temporal variations in hydraulic heads for use in model calibration. Ground water flow patterns indicated the convergence of flow along the four surface water features of the watershed; primary discharge areas were in the nontidal portions of the watershed. Ground water recharge zones corresponded to the surface water features with minimal development of a regional ground water system. Predicted ground water velocities varied between < 0.01 to 0.24 m/day, with elevated values associated with discharge areas and areas of convergence along surface water features. Some ground water residence times exceeded 100 years, although average residence times ranged between 16 and 21 years; approximately 95% of the ground water resource would reflect land use activities within the last 50 years.     

Pesticide residues in ground water of the San Joaquin Valley, California

A regional assessment of non-point-source contamination of pesticide residues in ground water was made of the San Joaquin Valley, an intensively farmed and irrigated structural trough in central California. About 10% of the total pesticide use in the USA is in the San Joaquin Valley. Pesticides detected include atrazine, bromacil, 2.4-DP, diazinon, dibromochloropropane, 1,2-dibromoethane, dicamba, 1,2-dichloropropane, diuron, prometon, prometryn, propazine and simazine. All are soil applied except diazinon.

Pesticide leaching is dependent on use patterns, soil texture, total organic carbon in soil, pesticide half-life and depth to water table. Leaching is enhanced by flood-irrigation methods except where the pesticide is foliar applied such as diazinon. Soils in the western San Joaquin Valley are fine grained and are derived primarily from marine shales of the Coast Ranges. Although shallow ground water is present, the fewest number of pesticides were detected in this region. The fine-grained soil inhibits pesticide leaching because of either low vertical permeability or high surface area; both enhance adsorption on to solid phases. Soils of the valley floor tend to be fine grained and have low vertical permeability. Soils in the eastern part of the valley are coarse grained with low total organic carbon and are derived from Sierra Nevada granites. Most pesticide leaching is in these alluvial soils, particularly in areas where depth to ground water is less than 30m. The areas currently most susceptible to pesticide leaching are eastern Fresno and Tulare Counties.

Tritium in water molecules is an indicator of aquifer recharge with water of recent origin. Pesticide residues transported as dissolved species were not detected in non-tritiated water. Although pesticides were not detected in all samples containing high tritium, these samples are indicative of the presence of recharge water that interacted with agricultural soils.     

Steady, annual, and monthly recharge implied by deep unconfined aquifer flow

We consider the response of a deep unconfined horizontal aquifer to steady, annual, and monthly recharge. A groundwater divide and a zero head reservoir constrain the aquifer, so that sinusoidal monthly and aperiodic annual recharge fluctuations create transient specific discharge near the reservoir and an unsteady water table elevation inland. One existing and two new long-term data sets from the Plymouth-Carver Aquifer in southeastern Massachusetts calibrate and confirm hydraulic properties in a set of analytical models. [Geohydrology and simulated groundwater flow, 1992] data and a new power law for tritiugenic helium to tritium ratios calibrate the steady recharge that drives the classical parabolic model of steady hydraulics [Applied Hydrogeology, 2001]. Observed water table and gradient fluctuations calibrate the transient recharge models. In the latter regard, monitoring wells within 1 km of Buttermilk Bay exhibit appreciable specific discharge and reduced water table fluctuations. We apply [Trans Am Geophys Union 32(1951)238] periodic model to the monthly hydraulics and a recharge convolution integral [J Hydrol 126(1991)315] to annual flow. An infiltration fraction of 0.79 and a consumptive use coefficient of 1.08×10⁻⁸ m/s °C relate recharge to precipitation and daylight weighted temperature across all three time scales. Errors associated with this recharge relation decrease with increasing time scale.     

The environmental tritium concentration of underground water and its hydrological interpretation

Two models which attempt to correlate the tritium concentration of water taken from aquifers to aquifer parameters are discussed. The first takes into account flow along individual streamlines and relates aquifer parameters to the observed tritium concentration at outflow. For the Gambier Plain unconfined aquifer in southern Australia, the calculated tritium concentration at outflow derived from known aquifer parameters and an environmental tritium input function for the area, is 0.7 T.U. The mean tritium concentration of several springs at outflow is also 0.7 T.U.Using the complete mixing model and the approximation that samples withdrawn from an aquifer on Eyre Peninsula, South Australia, are fully mixed, mean annual recharge for the area is estimated at 3 cm/year.     

Ground water discharge and nitrate flux to the Gulf of Mexico

Ground water samples (37 to 186 m depth) from Baldwin County, Alabama, are used to define the hydrogeology of Gulf coastal aquifers and calculate the subsurface discharge of nutrients to the Gulf of Mexico. The ground water flow and nitrate flux have been determined by linking ground water concentrations to 3H/3He and 4He age dates. The middle aquifer (A2) is an active flow system characterized by postnuclear tritium levels, moderate vertical velocities, and high nitrate concentrations. Ground water discharge could be an unaccounted source for nutrients in the coastal oceans. The aquifers annually discharge 1.1 +/- 0.01 x 10(8) moles of nitrate to the Gulf of Mexico, or 50% and 0.8% of the annual contributions from the Mobile-Alabama River System and the Mississippi River System, respectively. In southern Baldwin County, south of Loxley, increasing reliance on ground water in the deeper A3 aquifer requires accurate estimates of safe ground water withdrawal. This aquifer, partially confined by Pliocene clay above and Pensacola Clay below, is tritium dead and contains elevated 4He concentrations with no nitrate and estimated ground water ages from 100 to 7000 years. The isotopic composition and concentration of natural gas diffusing from the Pensacola Clay into the A3 aquifer aids in defining the deep ground water discharge. The highest 4He and CH4 concentrations are found only in the deepest sample (Gulf State Park), indicating that ground water flow into the Gulf of Mexico suppresses the natural gas plume. Using the shape of the CH4-He plume and the accumulation of 4He rate (2.2 +/- 0.8 microcc/kg/1000 years), we estimate the natural submarine discharge and the replenishment rate for the A3 aquifer.     

Simulation of submarine ground water discharge to a marine estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only approximately 10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.     

Ground water recharge and flow characterization using multiple isotopes

Stable isotopes of delta(18)O, delta(2)H, and (13)C, radiogenic isotopes of (14)C and (3)H, and ground water chemical compositions were used to distinguish ground water, recharge areas, and possible recharge processes in an arid zone, fault-bounded alluvial aquifer. Recharge mainly occurs through exposed stream channel beds as opposed to subsurface inflow along mountain fronts. This recharge distribution pattern may also occur in other fault-bounded aquifers, with important implications for conceptualization of ground water flow systems, development of ground water models, and ground water resource management. Ground water along the mountain front near the basin margins contains low delta(18)O, (14)C (percent modern carbon [pmC]), and (3)H (tritium units [TU]), suggesting older recharge. In addition, water levels lie at greater depths, and basin-bounding faults that locally act as a flow barrier may further reduce subsurface inflow into the aquifer along the mountain front. Chemical differences in ground water composition, attributed to varying aquifer mineralogy and recharge processes, further discriminate the basin-margin and the basin-center water. Direct recharge through the indurated sandstones and mudstones in the basin center is minimal. Modern recharge in the aquifer is mainly through the broad, exposed stream channel beds containing coarse sand and gravel where ground water contains higher delta(18)O, (14)C (pmC), and (3)H (TU). Spatial differences in delta(18)O, (14)C (pmC), and (3)H (TU) and occurrences of extensive mudstones in the basin center suggest sluggish ground water movement, including local compartmentalization of the flow system.     

Spatio-temporal analysis of potential aquifer recharge: Application to the Basin of Mexico

Regional estimates of aquifer recharge are needed in data-scarce regions such as the Basin of Mexico, where nearly 20 million people are located and where the Basin’s aquifer system represents the main water source. In order to develop the spatio-temporal estimates of aquifer recharge and to analyze to what extent urban growth has affected aquifer recharge, this work presents a daily soil water balance which uses different vegetation and soil types as well as the effect of topography on climatological variables and evapotranspiration. The soil water balance was applied on a daily time step in the Basin of Mexico for the period 1975–1986, obtaining an annually-lumped potential recharge flow of 10.9–23.8 m³/s (35.9–78.1 mm) in the entire Basin, while the monthly values for the year with the largest lumped recharge value (1981 = 78.1 mm) range from 1 m³/s (0.3 mm) in December to 87.9 m³/s (23.7 mm) in June. As aquifer recharge in the Basin mainly occurs by subsurface flow from its enclosing mountains as Mountain Block Recharge, urban growth has had a minimal impact on aquifer recharge, although it has diminished recharge in the alluvial plain.     

Ground water/surface water interaction in a fractured rock aquifer

In a recent field study of ground water/surface water interaction between a bedrock stream and an underlying fractured rock aquifer, it was determined that the majority of ground water discharge occurred through sparsely located vertical fractures. In this paper, the dominant mechanisms governing ground water/surface water exchange in such an environment are investigated using a numerical model. The study was conducted using several conceptual models based on the field study results. Although the field results provided the motivation for the modeling study, it was not intended to match modeling and field results directly. In addition, the extent of capture zones for discharging or recharging fractures was explored. The results of this study are intended to provide a better understanding of contaminant migration in the vicinity of bedrock streams. Based on the numerical results, the rate of ground water discharge (or recharge) was found to depend on the aperture size of the discharging feature, and on the distribution of hydraulic head with depth within the fracture network. It was determined that the extent of both the capture zone and reverse capture zone for an individual fracture can be extremely large, and will be determined by the height of the stream stage, the fracture apertures of the network, and the hydraulic-head distribution within the network. Because both the stream stage and the hydraulic-head distribution are transient, the size of the capture zone and/or the reverse capture zone for an individual fracture may change significantly over time. As a result, the migration path for contaminants within the fracture network and between the surface and subsurface will also vary significantly with time.     

Importance of river water recharge to the San Joaquin Valley groundwater system

Groundwater is not a sustainable resource, unless abstraction is balanced by recharge. Identifying the sources of recharge in a groundwater basin is critical for sustainable groundwater management. We studied the importance of river water recharge to groundwater in the south‐eastern San Joaquin Valley (24,000 km², population 4 million). We combined dissolved noble gas concentrations, stable isotopes, tritium, and carbon‐14 analyses to analyse the sources, mechanisms, and timescales of groundwater recharge. Area‐representative groundwater sampling and numerical model input data enabled a stable isotope mass balance and quantitative estimates of river and local recharge. River recharge, identified by a lighter stable isotope signature, represents 47 ± 4% of modern groundwater in the San Joaquin Valley (recharged after 1950) but only 26 ± 4% of premodern groundwater (recharged before 1950). This implies that the importance of river water recharge in the San Joaquin valley has nearly doubled and is likely the result of a 40% increase in total recharge, caused by river water irrigation return flows and increased stream depletion and river recharge due to groundwater pumping. Compared with the large and long‐duration capacity for water storage in the subsurface, storage of water in rivers is limited in time and volume, as evidenced by cold river recharge temperatures resulting from fast infiltration and recharge. Groundwater banking of seasonal surface water flows and expansion of managed aquifer recharge practices therefore appear to be a natural and promising method for increasing the resilience of the San Joaquin Valley water supply system.     

Utilization of decadal tritium variation for assessing the residence time of base flow

An iterative algorithm is presented that allows the user to model the subsurface residence time of shallow ground water comprising stream base flow based on decadal scale variation of tritium concentrations. The algorithm accounts for the effects of radioactive decay, the shallow subsurface mixing of ground water with precipitation, and ground water flux. The inverse of the best-fitting modeled flux through the saturated zone is equivalent to the residence time. The data required for this model include at least two measurements of tritium in base flow for a given stream location made at least a decade apart and the long-term tritium input in precipitation for the region of interest. The model is sensitive to relatively small changes in tritium concentrations and is limited by analytic uncertainties to an accuracy of approximately +/-5 years. The algorithm was applied to stream base flow for several basins in the Piedmont Province of Georgia in which tritium concentrations were measured during the early 1990s and again in the 2000s. The model results produced highly concordant residence times for three hydrogeologically similar basins in the Upper Ocmulgee Basin in North Central Georgia. The best estimate of the average residence time for ground water comprising base flow in this Piedmont basin using this new method is between approximately 14 and 18 years. These results are generally consistent with calculations made in previous studies, and these relatively long residence times can be attributed to the storage of water in the clay soils that dominate Piedmont Province watersheds.     

Inferring source waters from measurements of carbonate spring response to storms

We infer information about the nature of groundwater flow within a karst aquifer from the physical and chemical response of a spring to storm events. The spring discharges from the Maynardville Limestone in Bear Creek Valley, Tennessee. Initially, spring discharge peaks approximately 1–2 h from the midpoint of summer storms. The initial peak is likely due to surface loading, which pressurizes the aquifer and results in water moving out of storage. All of the storms monitored exhibited recessions that follow a master recession curve very closely, indicating that storm response is fairly consistent and repeatable, independent of the time between storms and the configuration of the rain event itself. Electrical conductivity initially increases for 0.5–2.9 days (longer for smaller storms), the result of moving older water out of storage. This is followed by a 2.1–2.5 day decrease in conductivity, resulting from an increasing portion of low conductivity recharge water entering the spring. Stable carbon isotope data and the calcite saturation index of the spring water also support this conceptual model. Spring flow is likely controlled by displaced water from the aquifer rather than by direct recharge through the soil zone.     

Hydrochemical and stable isotope evidence for the extent and nature of the effective Chalk aquifer of north Norfolk, UK

In eastern England the Chalk aquifer is covered by extensive Pleistocene deposits which influence the hydraulic conditions and hydrochemical nature of the underlying aquifer. In this study, the results of geophysical borehole logging of groundwater temperature and electrical conductivity and depth sampling for major ion concentrations and stable isotope compositions (δ¹⁸O and δ²H) are interpreted to reveal the extent and nature of the effective Chalk aquifer of north Norfolk. It is found that the Chalk aquifer can be divided into an upper region of fresh groundwater, with a Cl concentration of typically less than 100 mg l⁻¹, and a lower region of increasingly saline water. The transition between the two regions is approximately 50 m below sea-level, and results in an effective aquifer thickness of 50–60 m in the west of the area, but less than 25 m where the Eocene London Clay boundary is met in the east of the area. Hydrochemical variations in the effective aquifer are related to different hydraulic conditions developed in the Chalk. Where the Chalk is confined by low-permeability Chalky Boulder Clay, isotopically depleted groundwater (δ¹⁸O less than −7.5‰) is present, in contrast to those areas of unconfined Chalk where glacial deposits are thin or absent (δ¹⁸O about −7.0‰). The isotopically depleted groundwater is evidence for groundwater recharge during the late Pleistocene under conditions when mean surface air temperatures are estimated to have been 4.5°C cooler than at the present day, and suggests long groundwater residence times in the confined aquifer. Elevated molar Mg:Ca ratios of more than 0.2 resulting from progressive rock-water interaction in the confined aquifer also indicate long residence times. A conceptual hydrochemical model for the present situation proposes that isotopically depleted groundwater, occupying areas where confined groundwater dates from the late Pleistocene, is being slowly modified by both diffusion and downward infiltration of modem meteoric water and diffusive mixing from below with an old saline water body.     

Management of an aquifer with artificial recharge using water balance

Abstract

This study considers the management of an aquifer with artificial recharge whose water resources are mainly used for irrigation and where the irrigation return flows and natural discharge are considerable. These two variables are dependent on the use and increase of water resources. The interaction between some inflow and outflow variables is taken into account in the computation of the actual volume of water available, there being a pre-established artificial recharge of the aquifer. Conversely, the artificial recharge for a given water demand may be calculated. The most important inflow and outflow variables are used in an equation system (the water balance for each period of time) for these interactions and the total inflow and outflow. An iterative procedure is developed that solves the equation and with it the computation of the actual availability of water and storage.     

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.     

Radiocarbon and δ13C Values Related to Ground-Water Recharge and Mixing

Ground-water levels in the Upper Floridan aquifer beneath the southeastern coast of South Carolina have undergone pumpage-induced declines approaching 20 ft below sea level at the southern end of Hilton Head Island. This scenario suggests the potential exists for the inducement of recharge to the Upper Floridan aquifer across the island, which could affect the quality of water being pumped by wells. However, low radiocarbon concentrations in ground-water samples (0.5 to 1.4 ± 0.1 PMC) indicate that most of the water is relict ground water reflecting prepumpage ground-water flow conditions in the Upper Floridan aquifer. The isotopic data indicate long residence times and water-chemistry evolution more characteristic of ground-water recharge occurring farther inland prior to the commencement of pumpage in the late 1800s. Radiocarbon concentrations (as Percent Modern Carbon) and stable carbon isotope ratios (as δ¹³C in dissolved inorganic carbon) determined during this study and reported in other studies on and around Hilton Head Island varied in a systematic manner. Heavier δ¹³C values (–2.8 to –1.6 per mil) in ground water beneath southern Hilton Head Island reflect ground-water discharge from prepumpage flowpaths originating over 100 miles away, hence a depletion in radiocarbon concentration with corrected ground-water ages no younger than 16,000 yrs BP. In contrast, lighter δ¹³C values (–13.9 to –8.67 per mil) beneath the northern part of the island indicate recent recharge as a result of water-level declines, and recharge in areas off the island that have not changed as a result of pumpage (evidenced by enrichment in radiocarbon with corrected ground-water ages no older than 4,000 yrs BP). This suggests that the δ¹³C composition of ground water in the Upper Floridan aquifer is a useful indicator of mixing between ground waters from different sources, and can be used to delineate recharge-discharge patterns. This approach may be applicable to other aquifers of highly evolved ground-water chemistry in regional carbonate aquifer systems that may be receiving recent recharge. Moreover, this approach could prove useful in delineating the contribution of recent water being captured by pumped wells as part of wellhead protection programs designed to assess aquifer vulnerability from surficial contaminant sources.     

Arsenic Control During Aquifer Storage Recovery Cycle Tests in the Floridan Aquifer

Implementation of aquifer storage recovery (ASR) for water resource management in Florida is impeded by arsenic mobilization. Arsenic, released by pyrite oxidation during the recharge phase, sometimes results in groundwater concentrations that exceed the 10 µg/L criterion defined in the Safe Drinking Water Act. ASR was proposed as a major storage component for the Comprehensive Everglades Restoration Plan (CERP), in which excess surface water is stored during the wet season, and then distributed during the dry season for ecosystem restoration. To evaluate ASR system performance for CERP goals, three cycle tests were conducted, with extensive water‐quality monitoring in the Upper Floridan Aquifer (UFA) at the Kissimmee River ASR (KRASR) pilot system. During each cycle test, redox evolution from sub‐oxic to sulfate‐reducing conditions occurs in the UFA storage zone, as indicated by decreasing Fe²⁺/H₂S mass ratios. Arsenic, released by pyrite oxidation during recharge, is sequestered during storage and recovery by co‐precipitation with iron sulfide. Mineral saturation indices indicate that amorphous iron oxide (a sorption surface for arsenic) is stable only during oxic and sub‐oxic conditions of the recharge phase, but iron sulfide (which co‐precipitates arsenic) is stable during the sulfate‐reducing conditions of the storage and recovery phases. Resultant arsenic concentrations in recovered water are below the 10 µg/L regulatory criterion during cycle tests 2 and 3. The arsenic sequestration process is appropriate for other ASR systems that recharge treated surface water into a sulfate‐reducing aquifer.     

The water budget myth revisited: why hydrogeologists model

Within the ground water community, the idea persists that if one can estimate the recharge to a ground water system, one then can determine the size of a sustainable development. Theis addressed this idea in 1940 and showed it to be wrong-yet the myth continues. The size of a sustainable ground water development usually depends on how much of the discharge from the system can be "captured" by the development. Capture is independent of the recharge; it depends on the dynamic response of the aquifer system to the development. Ground water models were created to study the response dynamics of ground water systems; it is one of the principal reasons hydrogeologists model.     

Ground water flow parameterization of an Appalachian coal mine complex

We examined a large (240 km2) northern Appalachian bituminous coal basin (Irwin Syncline, Westmoreland County, Pennsylvania) comprising 27 mine complexes with nine major (> 2.5 x 10(3) L/min) discharges. The synclinal basin was divided into seven subbasins based on equilibrium hydraulic relationships established during the past 25 years. Recharge rates, mine pool velocity, and residence times respond to hydraulic changes in the overburden induced by mine subsidence. The estimated maximum depth for subsidence fractures is 60 m (30 times mined thickness) with recharge rates decreasing significantly in subbasins with thicker overburden (> 75 m). Calculated subbasin recharge rates range from 2 to 6 x 10(-4) L/min/m2 and are significantly lower than the previously used rate for the basin. Residence time of ground water in the Irwin subbasins calculated using average linear velocity ranged from one to five years and were more consistent with field observations than estimates obtained using discharge and basin volume area. A positive correlation (r2 = 0.80) exists between net alkalinity of the mine water-impacted discharges and residence time in the mine pools. Our results for the Irwin coal basin suggest that use of a subbasin approach incorporating overburden depth, mining methodology, and the extent of postmining inundation will lead to improved determination of ground water flow parameters in mined watersheds in northern Appalachia and elsewhere.