Seasonal mixing from intermittent flow drives concentration-discharge behaviour in a stream affected by coal mine drainage

Abandoned mining operations continue to severely degrade many ecosystems worldwide by releasing acidic water and/or heavy metals into surface and groundwater. Contaminant concentrations in affected streams vary with discharge in patterns that reflect both geochemical reactions and variable mixing of contaminated and non-contaminated waters. However, controls on concentration-discharge (C-Q) patterns remain unclear, particularly for constituents that experience changing solubility across redox and pH gradients. Understanding the C-Q behaviour of contaminants aids in predicting both downstream transport and effects on aquatic life under variable flow. Here, we examined the C-Q behaviours of non-reactive (Na, K, Ca, Mg, Cl⁻) and reactive (Fe, Mn, Al, H⁺, SO₄²⁻) solutes in a stream contaminated with acid mine drainage in northeastern Ohio, USA. Concentration-discharge patterns at the watershed outlet primarily reflected mixing of contaminated baseflow with intermittent inputs of high pH water draining from a passive limestone treatment system into the stream. The treatment system acted as an ephemeral tributary that mitigated contamination in the stream by diluting solutes, raising pH, and driving metal precipitation, but only when flow was present during wet seasons. Consequently, AMD-derived reactive solutes (H⁺, Fe, Mn, Al) decreased with increasing stream discharge while relatively conservative solutes (e.g., Ca, Mg, K, Na) decreased only slightly or were chemostatic. This study highlights both the unique C-Q patterns of reactive solutes when compared to those of non-reactive solutes and the potential for intermittent streams to control C-Q behaviour in headwater catchments.     

Hydraulic subsurface measurements and hydrodynamic modelling as indicators for groundwater flow systems in the Rotondo granite,Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man‐made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high‐relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of Seasonal Pumping on Groundwater Sources and Flow System,Nagaoka Plain,Japan

Intensive groundwater development in the urban area of the Nagaoka Plain, Japan, has induced changes in the pH and saturation index of calcite in groundwater. To account for these chemical changes, it is important to determine seasonal variations of recharge and the groundwater flow system in the aquifer. This study identified the sources and flow system of groundwater in this urban area by a comprehensive method using stable isotope data and a numerical groundwater model of the Nagaoka Plain. Stable isotope evidence shows that the groundwater is recharged by meteoric water originating from low‐elevation areas rather than the mountains surrounding the plain. The water table in the study area is drawn down during the winter and recovers in the other seasons. Numerical modeling shows that discharge occurs primarily along the Shinano River during the recovery period, whereas discharge is centered in urbanized areas during the drawdown period, when a conical depression of the water table stimulates recharge from the immediate area. These results are indications of a local groundwater flow system, with its recharge area between the Shinano River and the urban areas, which is governed by intensive seasonal groundwater extraction.     

Spatial patterns of groundwater‐lake exchange – implications for acid neutralization processes in an acid mine lake

Exchange of groundwater and lake water with typically quite different chemical composition is an important driver for biogeochemical processes at the groundwater‐lake interface, which can affect the water quality of lakes. This is of particular relevance in mine lakes where anoxic and slightly acidic groundwater mixes with oxic and acidic lake water (pH < 3). To identify links between groundwater‐lake exchange rates and acid neutralization processes in the sediments, exchange rates were quantified and related to pore‐water pH, sulfate and iron concentrations as well as sulfate reduction rates within the sediment. Seepage rates measured with seepage meters (?2.5 to 5.8 L m^‐2 d^‐1) were in reasonable agreement with rates inverted from modeled chloride profiles (?1.8 to 8.1 L m^‐2 d^‐1). Large‐scale exchange patterns were defined by the (hydro)geologic setting but superimposed by smaller scale variations caused by variability in sediment texture. Sites characterized by groundwater upwelling (flow into the lake) and sites where flow alternated between upwelling and downwelling were identified. Observed chloride profiles at the alternating sites reflected the transient flow regime. Seepage direction, as well as seepage rate, were found to influence pH, sulfate and iron profiles and the associated sulfate reduction rates. Under alternating conditions proton‐consuming processes, for example, sulfate reduction, were slowed. In the uppermost layer of the sediment (max. 5 cm), sulfate reduction rates were significantly higher at upwelling (>330 nmol g^‐1 d^‐1) compared to alternating sites (<220 nmol g^‐1 d^‐1). Although differences in sulfate reduction rates could not be explained solely by different flux rates, they were clearly related to the prevailing groundwater‐lake exchange patterns and the associated pH conditions. Our findings strongly suggest that groundwater‐lake exchange has significant effects on the biogeochemical processes that are coupled to sulfate reduction such as acidity retention and precipitation of iron sulfides. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessment of groundwater resources of Quaternary aquifer system in concordance with avoidable negative impact on geoenvironment,south-eastern part of Lithuania

This paper focuses upon south-eastern Lithuanian Quaternary aquifer system groundwater resources formation modelling. Groundwater model calibration has been performed for a pre-development and transient flow conditions. The results demonstrate that there is an intense interaction between groundwater and surface water bodies which form groundwater resources and runoff. For Quaternary cover the majority of unconfined groundwater outflows to surface water streams, the remaining part discharges through the confined interglacial/interstadial aquifers and lateral outflows across the boundaries. Groundwater prognostic exploitable resources can be obtained without a significant negative impact on the geoenvironment. The main sources of exploitable resources formation are increase in groundwater recharge and lateral inflow (40.7%), decreased outflow via streams and lateral boundaries (41.2%) and drawing up additional flow from the hydraulically connected streams (14.5%).     

Comparison of two modeling approaches for groundwater–surface water interactions

An assessment of interactions between groundwater and surface water was carried out by applying two different modeling approaches to a small‐scale study area in the municipality of Havelock, Quebec. The first approach involved a commonly used sequential procedure that consists in determining the daily recharge rate using a quasi 2D infiltration model (HELP), applied in the next step as an imposed flux to a 3D finite‐element groundwater flow model. The flow model was calibrated under steady‐state and transient conditions against measured water levels. The second approach was based on a recently developed physically based, 3D fully coupled groundwater–surface water flow model (CATHY) applied to the entire flow domain in an integrated manner. Implementation, calibration, and results of the simulations for both approaches are presented and discussed. For equal annual precipitation (1038 mm/y) and evapotranspiration (556 mm/y), the second approach computed a recharge rate of 233 mm/y (8.9% higher than the first approach) and a net upward flow from the fractured aquifer (the first approach predicted a net downward flow to the rock). The simulated annual discharge was similar for the two approaches (9.6% difference). Both approaches were found to be useful in understanding the interactions between groundwater and surface water, although limitations are apparent in the sequential procedure's inability to account for surface–subsurface feedbacks, for instance near stream reaches where groundwater discharge is prevalent. The decoupled, two‐model approach provides disaggregated surface, vadose, and aquifer flows, and a simple aperçu at the different components of total discharge. The fully coupled model accounts for continuous water exchanges between the land surface, subsurface, and stream channel in a more complex manner, and produces a better match against observed data. Copyright © 2012 John Wiley & Sons, Ltd.     

Changes in Projected Spatial and Seasonal Groundwater Recharge in the Upper Colorado River Basin

The Colorado River is an important source of water in the western United States, supplying the needs of more than 38 million people in the United States and Mexico. Groundwater discharge to streams has been shown to be a critical component of streamflow in the Upper Colorado River Basin (UCRB), particularly during low‐flow periods. Understanding impacts on groundwater in the basin from projected climate change will assist water managers in the region in planning for potential changes in the river and groundwater system. A previous study on changes in basin‐wide groundwater recharge in the UCRB under projected climate change found substantial increases in temperature, moderate increases in precipitation, and mostly periods of stable or slight increases in simulated groundwater recharge through 2099. This study quantifies projected spatial and seasonal changes in groundwater recharge within the UCRB from recent historical (1950 to 2015) through future (2016 to 2099) time periods, using a distributed‐parameter groundwater recharge model with downscaled climate data from 97 Coupled Model Intercomparison Project Phase 5 (CMIP5) climate projections. Simulation results indicate that projected increases in basin‐wide recharge of up to 15% are not distributed uniformly within the basin or throughout the year. Northernmost subregions within the UCRB are projected an increase in groundwater recharge, while recharge in other mainly southern subregions will decline. Seasonal changes in recharge also are projected within the UCRB, with decreases of 50% or more in summer months and increases of 50% or more in winter months for all subregions, and increases of 10% or more in spring months for many subregions.     

In situ hydraulic tests in the active fault survey tunnel, Kamioka Mine, excavated through the active Mozumi-Sukenobu Fault zone and their hydrogeological significance

Abstract   The spatial hydrogeological and structural character of the active Mozumi-Sukenobu Fault (MSF) was investigated along a survey tunnel excavated through the MSF in the Kamioka Mine, central Japan. Major groundwater conduits on both sides of the MSF are recognized. One is considered to be a subvertical conduit between the tunnel and the surface, and the other is estimated to be a major reservoir of old meteoric water alongside the MSF. Our studies indicate that part of the MSF is a sub-vertical continuous barrier that obstructs younger meteoric water observed in the south-eastern part of the Active Fault Survey Tunnel (AFST) and water recharge to the rock mass intersected by the north-western part of the AFST. The MSF might be a continuous barrier resulting in the storage of a large quantity of older groundwater to the northwest. The observations and results of in situ hydraulic tests indicate that the major reservoir is not the fault breccia associated with the northeast–southwest trending faults of the MSF, but the zone in which blocks of fractured rocks occur beside high angle faults corresponding to X shears whose tectonic stress field coincides with the present regional stress field and antithetic Riedel shears of the MSF. The results from borehole investigations in the AFST indicate that secondary porosity is developed in the major reservoir due to the destruction of filling minerals and fracture development beside these shears. The increase in hydraulic conductivity is not directly related to increased density of fractures around the MSF. Development of secondary porosity could cause the increase in hydraulic conductivity around the MSF. Our results suggest that minor conduits of the fracture network are sporadically distributed in the sedimentary rocks around the MSF in the AFST.     

Groundwater discharge mechanism in semi‐arid regions and the role of evapotranspiration

The present study examined groundwater recharge/discharge mechanisms in the regional Central Sudan Rift Basins (CSRB). Aquifers in CSRB constitute poorly sorted silisiclastics of sand, clay and gravels deposited in closed hydrologic systems of the Cretaceous–Pleistocene fluviolacustrine environments. CSRB are bounded to the north by the highlands of the Central African Shear Zone (CAZS) that represents the surface and groundwater divides. Sporadic recharge in the peripheries of the basins along the CASZ occurs subsequent to decadal and centennial storm events. Inflow from the Nile into the aquifers represents an additional source of recharge. Thus, groundwater resources cannot be labelled fossil nor can they be readily recharged. Closed hydrologic troughs located adjacent to the influent Nile system mark areas of main groundwater discharge characterized by lower hydraulic heads. This study has examined mechanisms that derive the discharge of the groundwater in these closed basins and concluded that only evapotranspirative discharge can provide a plausible explanation. Groundwater abstraction is mainly through deep‐rooted trees and effective evaporation. The increase of TDS along the flow indicates local recharge at the peripheries of basins and shows the influence of evaporation and rock/water interaction. The decline in groundwater level along a flow path was calculated using Darcy's law to estimate average recharge and evapotranspirative discharge, which are equal under natural equilibrium and make the only fluxes in CSRB. Steady‐state 2D flow modelling has demonstrated that an average recharge of 4–8 mm yr^?1 and evapotranspirative discharge of 1–22 mm yr^?1 will maintain natural equilibrium in CSRB. Sporadic storms provide recharge in the highlands to preserve the current hydraulic gradient and maintain aquifer dynamics. Simulated recharge from the Nile totals about 17·5 mm yr^?1 and is therefore a significant contributor to the water balance. Copyright © 2007 John Wiley & Sons, Ltd.     

Groundwater Renewal in the Middle Odra River Catchment Based on Flow Model and Isotopic Data – the Lubin Głogów Copper Region,Poland

In the area of intense mining within the Lubin Głogów Copper Region the investigation of groundwater renewal was made according to two methods: a numerical flow model and isotopic analysis. The results of model simulations show that within the shallow aquifer the groundwater renewal is quick and mostly depends on infiltration within the recharge zones. An average time of groundwater flow between recharge and discharge areas in the representative part of the entire system has been obtained as 130 years. The content of isotopes was determined in samples collected as well from the mine excavations as from water intakes. The measurements were done on 17 samples. Results of isotopic studies indicate that the relationship of δD – δ18O in groundwater gathers along the established local meteoric water line (MWL) δD = 7.58 δ18O + 4.45. Deep groundwaters that create an inflow to the mines were formed by infiltration in Holocene with reference to two of total four mines and probably in Eopleistocene in another two.     

Hydrodynamic and isotopic investigations for evaluating the mechanisms and amount of groundwater seepage through a rockslide dam

This study addresses the influence of landslide dams on surface water drainage and groundwater flow. In the study area of Scanno Lake and Sagittario River (Central Italy), a limestone rockslide‐avalanche formed a lake, which has an outlet that is occasionally active, showing infiltration into the rockslide dam. Several springs are present at the lake's base and are partly fed by seepage through the rockslide debris. Piezometric surveys, discharge measurements, pumping tests and chemical analyses are tools used to build a conceptual model of the groundwater flow and to evaluate the flow through the rockslide debris. Seasonal water isotopic signatures validate the assumed model, showing a mixing of infiltration recharge and groundwater seepage throughout the rockslide debris. Various recharge areas have been found for springs, pointing out those directly fed by the rockslide debris aquifer. Hypotheses about seasonal groundwater mixing between the regional carbonate aquifer and the rockslide debris aquifer are supported by isotope results. Seasonal changes in groundwater table level due to recharge and surface losses from seasonal outlet have been correlated with isotopic groundwater composition from the rockslide debris aquifer and the downstream springs; this relationship highlights the role of the rockslide dam body on the hydrodynamics of the studied area. Relationships between surface waters and groundwater in the area have been completely understood on the basis of water isotopic fingerprinting, finally obtaining a complete evaluation of groundwater renewable resources and its regimen. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrogeochemical Evolution Along Groundwater Flow Paths in the Manas River Basin,Northwest China

The impacts of long-term pumping on groundwater chemistry remain unclear in the Manas River Basin, Northwest China. In this study, major ions within five surface water and 105 groundwater samples were analyzed to identify hydrogeochemical processes affecting groundwater composition and evolution along the regional-scale groundwater flow paths using the multivariate techniques of hierarchical cluster analysis (HCA) and principal components analysis (PCA) and traditional graphical methods for analyzing groundwater geochemistry. HCA classified the groundwater samples into four clusters (C1 to C4). PCA reduced the dimensionality of geochemical data into three PCs, which explained 86% of the total variance. The results of HCA and PCA were used to identify three zones: “recharge,” “transition,” and “discharge.” In the recharge zone the groundwater type is Ca-HCO₃-SO₄ and is primarily impacted by the dissolution of calcite and silicate weathering. In the transition zone the groundwater type is Ca-HCO₃-SO₄-Cl and is impacted by rock dissolution and reverse ion exchange. In the discharge zone the groundwater type is Na-Cl and is impacted by evaporation and reverse ion exchange. In addition, anthropogenic activities impact the groundwater chemistry in the study area. The groundwater type generally changes from Ca-HCO₃-SO₄ in the recharge area to Na-Cl in the discharge area along the regional-scale groundwater flow paths. This study provides a process-based knowledge for understanding the interaction of groundwater flow patterns and geochemical evolution within the Manas River Basin.     

Groundwater–surface water interaction in Lake Nasser,Southern Egypt

A cross‐sectional model, based on the two dimensional groundwater flow equation of Edelman, was applied at seven transects distributed over four geological cross sections to estimate groundwater heads and recharge from/or groundwater discharge to Lake Nasser. The lake with a length of 500 km and an average width of 12 km was created over the period 1964–1970, the time for constructing the Aswan High Dam (AHD). The model, constrained by regional‐scale groundwater flow and groundwater head data in the vicinity of the lake, was successfully calibrated to timeseries of piezometeric heads collected at the cross sections in the period 1965–2004. Inverse modeling yielded high values for the horizontal hydraulic conductivity in the range of 6.0 to 31.1 m day^?1 and storage coefficient between 0.01 and 0.40. The results showed the existence of a strong vertical anisotropy of the aquifer. The calibrated horizontal permeability is systematically higher than the vertical permeability (≈1000:1). The calibrated model was used to explore the recharge from/or groundwater discharge to Lake Nasser at the seven transects for a 40‐year period, i.e. from 1965 to 2004. The analysis for the last 20‐year period, 1985–2004, revealed that recharge from Lake Nasser reduced by 37% compared to the estimates for the first 20‐year period, 1965–1984. In the period 1965–2004, seepage of Lake Nasser to the surrounding was estimated at 1.15 × 10⁹ m³ year^?1. This led to a significant rise of the groundwater table. Variance‐based sensitivity and uncertainty analysis on the Edelman results were conducted applying quasi‐Monte Carlo sequences (Latin Hypercube sampling). The maximum standard deviation of the total uncertainty on the groundwater table was 0.88 m at Toshka (west of the lake). The distance from the lake, followed by the storage coefficient and hydraulic conductivity, were identified as the most sensitive parameters. Copyright © 2012 John Wiley & Sons, Ltd.     

Minor improvement for intertidal seaweeds and invertebrates after acid mine drainage diversion at Britannia Beach, Pacific Canada

In December 2001, acid mine drainage (AMD) from an abandoned copper mine at Britannia Beach (British Columbia, Canada) was diverted to flow from Britannia Creek into an outfall at 30 m depth in Howe Sound. Britannia Beach was studied in early 2003 to determine whether AMD diversion resulted in improved conditions for intertidal organisms. Species number and abundance have increased at the intertidal zone since AMD diversion, although they were still lower than at an unpolluted control site nearby (Furry Creek). Survivorship and growth rates of transplanted Mytilus trossulus (mussel) have increased since AMD diversion, although they were still significantly lower than at the control site. Transplanted Fucus gardneri (seaweed) performed better than before the AMD diversion; at Britannia Beach the chlorophyll a concentration in tissues was not significantly different from that at the control site, although the concentration of chlorophyll c in tissues and the chlorophyll c to a ratio was lower than at the control site six weeks after transplantation. Britannia Beach is still subject to leaching of metals from surrounding soils, low levels of AMD coming down the creek, and AMD discharge from the deep outfall. Although there has been an improvement, the intertidal environment at Britannia Beach still seems unable to support normal growth and survival of organisms.     

Evidence of karst from behaviour of the Asmari limestone aquifer at the Khersan 3 Dam site,southern Iran

Abstract

The Asmari limestone formation is the major aquifer system at the Khersan 3 Dam site, Zagros, Iran. Characterization of the aquifer system and study of karst development are essential for forecasting leakage potential and to plan remediation works. The aquifer functioning and karst structure were evaluated by geology, well hydrodynamics and natural tracing studies, showing one unconfined and three artesian sub-aquifers, the last being characterized by rapid flow, with contributing old recharge water, and a recharge at higher elevation than the unconfined sub-aquifer. The anticline structure outcrops the artesian sub-aquifers downstream of the site. The confining layers disappear and the aquifer system discharges as a mix of all groundwater. Artesian groundwater is estimated to contribute about 80% (dry season) and 50% (wet season) of flow in the mixing zone. Very few karst features are observed in cores and galleries at the site, despite some karst landforms in the reservoir area. However the aquifer exhibits fast flow of karst type in the artesian sub-aquifers.     

Propagation of drought through groundwater—a new approach using linear reservoir theory

The effect of drought on groundwater heads and discharge is often complex and poorly understood. Therefore the propagation of a drought from groundwater recharge to discharge and the influence of aquifer characteristics on the propagation was analysed by tracking a drought in recharge through a linear reservoir. The recharge was defined as a sinusoid function with a period of 1 year. The decrease in recharge owing to drought was simulated by multiplying the recharge during 1 year with a drought fraction between 0 and 1, which represents a decrease in the recharge of 100 to 0%, respectively. The droughts were identified using the threshold level approach, with a threshold that is constant in time. For this case analytical formulations were derived, which express the drought duration and deficit in the groundwater discharge in terms of the decrease in recharge, the reservoir coefficient that characterizes aquifer properties and the height of the threshold level. The results showed that the delay in the groundwater system caused a shift of the main part of the decrease in recharge from the high‐flow to the low‐flow period. This resulted in an increase in drought deficit for discharge compared with the drought deficit for recharge. Also the development of multiyear droughts caused an increase in drought deficit. The attenuation in the groundwater system caused a decrease in drought deficit. In most cases the net effect of these processes was an increase of drought deficit as a result of the propagation through groundwater. Only for small droughts the deficit decreased from recharge to discharge. The amount of increase or decrease depends on the reservoir coefficient and the severity of the drought. Under most conditions a maximum in the drought deficit occurred for a reservoir coefficient of around 200 days. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Effects of urbanization on streamflow in the Atlanta area (Georgia,USA): a comparative hydrological approach

For the period from 1958 to 1996, streamflow characteristics of a highly urbanized watershed were compared with less‐urbanized and non‐urbanized watersheds within a 20 000 km² region in the vicinity of Atlanta, Georgia: in the Piedmont and Blue Ridge physiographic provinces of the southeastern USA. Water levels in several wells completed in surficial and crystalline‐rock aquifers were also evaluated. Data were analysed for seven US Geological Survey (USGS) stream gauges, 17 National Weather Service rain gauges, and five USGS monitoring wells. Annual runoff coefficients (RCs; runoff as a fractional percentage of precipitation) for the urban stream (Peachtree Creek) were not significantly greater than for the less‐urbanized watersheds. The RCs for some streams were similar to others and the similar streams were grouped according to location. The RCs decreased from the higher elevation and higher relief watersheds to the lower elevation and lower relief watersheds: values were 0·54 for the two Blue Ridge streams, 0·37 for the four middle Piedmont streams (near Atlanta), and 0·28 for a southern Piedmont stream. For the 25 largest stormflows, the peak flows for Peachtree Creek were 30% to 100% greater than peak flows for the other streams. The storm recession period for the urban stream was 1–2 days less than that for the other streams and the recession was characterized by a 2‐day storm recession constant that was, on average, 40 to 100% greater, i.e. streamflow decreased more rapidly than for the other streams. Baseflow recession constants ranged from 35 to 40% lower for Peachtree Creek than for the other streams; this is attributed to lower evapotranspiration losses, which result in a smaller change in groundwater storage than in the less‐urbanized watersheds. Low flow of Peachtree Creek ranged from 25 to 35% less than the other streams, possibly the result of decreased infiltration caused by the more efficient routing of stormwater and the paving of groundwater recharge areas. The timing of daily or monthly groundwater‐level fluctuations was similar annually in each well, reflecting the seasonal recharge. Although water‐level monitoring only began in the 1980s for the two urban wells, water levels displayed a notable decline compared with non‐urban wells since then; this is attributed to decreased groundwater recharge in the urban watersheds due to increased imperviousness and related rapid storm runoff. Copyright © 2001 John Wiley & Sons, Ltd.     

Quantifying changes of effective springshed area and net recharge through recession analysis of spring flow

The average flow of Silver Springs, one of the largest magnitude springs in Central Florida, declined 32% from 2000 to 2012. The average groundwater head in the springshed declined 0.14 m, and the spring pool altitude increased 0.24 m during the same period. This paper presents a novel explanation of the spring flow recession curve for Silver Springs using the Torricelli model, which uses the groundwater head at a sentinel well, the spring pool altitude and the net recharge to groundwater. The effective springshed area and net recharge (defined as recharge minus groundwater pumping and evapotranspiration) were estimated based on the observed recession slopes for spring flow, groundwater head and spring pool altitude. The results indicate that the effective springshed area continuously declined since 1989 and the net recharge declined since the 1970s with a significant drop in 2002. Subsequent to 2002, the net recharge increased modestly but not to the levels prior to the 1990s. The reduction in net recharge was caused by changes in hydroclimatic conditions including precipitation and air temperature, along with groundwater withdrawals, which contributed to the declined spring flow. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.