Historical Groundwater Trends in Northern New England and Relations with Streamflow and Climatic Variables

Water‐level trends spanning 20, 30, 40, and 50 years were tested using month‐end groundwater levels in 26, 12, 10, and 3 wells in northern New England (Maine, New Hampshire, and Vermont), respectively. Groundwater levels for 77 wells were used in interannual correlations with meteorological and hydrologic variables related to groundwater. Trends in the contemporary groundwater record (20 and 30 years) indicate increases (rises) or no substantial change in groundwater levels in all months for most wells throughout northern New England. The highest percentage of increasing 20‐year trends was in February through March, May through August, and October through November. Forty‐year trend results were mixed, whereas 50‐year trends indicated increasing groundwater levels. Whereas most monthly groundwater levels correlate strongly with the previous month's level, monthly levels also correlate strongly with monthly streamflows in the same month; correlations of levels with monthly precipitation are less frequent and weaker than those with streamflow. Groundwater levels in May through August correlate strongly with annual (water year) streamflow. Correlations of groundwater levels with streamflow data and the relative richness of 50‐ to 100‐year historical streamflow data suggest useful proxies for quantifying historical groundwater levels in light of the relatively short and fragmented groundwater data records presently available.     

Surface Water Seepage Effects on Shallow Ground‐Water Quality Along the Rio Grande in Northern New Mexico1

Abstract: Interactions between surface irrigation water, shallow ground water, and river water may have effects on water quality that are important for both drinking water supplies and the ecological function of rivers and floodplains. We investigated water quality in surface water and ground water, and how water quality is influenced by surface water inputs from an unlined irrigation system in the Alcalde Valley of the Rio Grande in northern New Mexico. From August 2005 to July 2006, we sampled ground water and surface water monthly and analyzed for concentrations of major cations and anions, specific conductance, pH, dissolved oxygen, and water levels. Results indicate that irrigation ditch seepage caused an increase in ground water levels and that the Rio Grande is a gaining stream in this region. Temporal and spatial differences were found in ion concentrations in shallow ground water as it flowed from under the ditch toward the river. Ground‐water ion concentrations were higher when the ditch was not flowing compared with periods during peak irrigation season when the ditch was flowing. Ditch inputs diluted ion concentrations in shallow ground water at well positions near the ditch. Specifically, lower ion concentrations were detected in ground water at well positions located near the ditch and river compared with well positions located in the middle of an agricultural field. Results from this project showed that ditch inputs influenced ion concentrations and were associated with ground‐water recharge. In arid region river valleys, careful consideration should be given to management scenarios that change seepage from irrigation systems, because in some situations reduced seepage could negatively affect ground‐water recharge and water quality.     

Controls on Temperature in Salmonid‐Bearing Headwater Streams in Two Common Hydrogeologic Settings,Kenai Peninsula,Alaska

Headwater streams are the most numerous in terms of both number and length in the conterminous United States and play important roles as spawning and rearing grounds for numerous species of anadromous fish. Stream temperature is a controlling variable for many physical, chemical, and biological processes and plays a critical role in the overall health and integrity of a stream. We investigated the controls on stream temperature in salmon‐bearing headwater streams in two common hydrogeologic settings on the Kenai Peninsula, Alaska: (1) drainage‐ways, which are low‐gradient streams that flow through broad valleys; and (2) discharge‐slopes, which are high gradient streams that flow through narrow valleys. We hypothesize local geomorphology strongly influences surface‐water and groundwater interactions, which control streamflow at the network scale and stream temperatures at the reach scale. The results of this study showed significant differences in stream temperatures between the two hydrogeologic settings. Observed stream temperatures were higher in drainage‐way sites than in discharge‐slope sites, and showed strong correlations as a continuous function with the calculated topographic metric flow‐weighted slope. Additionally, modeling results indicated the potential for groundwater discharge to moderate stream temperature is not equal between the two hydrogeologic settings, with groundwater having a greater moderating effect on stream temperature at the drainage‐way sites.     

Specific Yield Functions for Estimating Evapotranspiration from Diurnal Surface Water Cycles

The White method has been routinely used to estimate evapotranspiration using diurnal variations in groundwater levels. Applications to surface water systems (e.g., wetlands) are less common. For applications to surface water systems, a stage‐dependent specific yield function must be defined. This is especially important for small wetlands formed in topographic depressions with bowl shaped bathymetries. Existing formulations of the specific yield function include weighting factors that impact the relative importance of the soil and open water specific yields on the composite value. Three formulations of the specific yield function from the literature were compared and found to produce varied results. Based on a comparison with empirical estimates of specific yield based on observed ratios of net precipitation to water level rise, one of the existing formulations is generalized and recommended for general use. The recommended function is dependent on wetland bathymetry, magnitude of the diurnal fluctuation, spatial extent of the equilibration area, and soil‐specific yield. A sensitivity analysis was conducted to examine the relative importance of these variables. The specific yield function is independent of wetland size and is strongly dependent on the basin profile coefficient (p), an indication of wetland shape. For most natural wetlands, bathymetry strongly influences specific yield.     

Water Balances of Two Piedmont Headwater Catchments: Implications for Regional Hydrologic Landscape Classification

In the Piedmont of North Carolina, a traditionally water‐rich region, reservoirs that serve over 1 million people are under increasing pressure due to naturally occurring droughts and increasing land development. Innovative development approaches aim to maintain hydrologic conditions of the undisturbed landscape, but are based on insufficient target information. This study uses the hydrologic landscape concept to evaluate reference hydrology in small headwater catchments surrounding Falls Lake, a reservoir serving Raleigh and the greater Triangle area. Researchers collected one year of detailed data on water balance components, including precipitation, evapotranspiration, streamflow, and shallow subsurface storage from two headwater catchments representative of two hydrologic landscapes defined by differences in soils and topographic characteristics. The two catchments are similar in size and lie within the same physiographic region, and during the study period they showed similar water balances of 26‐30% Q, ?4 to 5% ΔS, 59‐65% evapotranspiration, and 9‐10% G. However, the steeper, more elevated catchment exhibited perennial streamflow and nongrowing season runoff ratios (Q/P) of 33%, whereas the flat, low‐lying stream was drier during the growing season and exhibited Q/P ratios of 52% during the nongrowing season. A hydrologic landscape defined by topography and soil characteristics helps characterize local‐scale reference hydrology and may contribute to better land management decisions.     

Instream Restoration to Improve the Ecohydrologic Function of a Subalpine Meadow: Pre‐implementation Modeling with HEC‐RAS

Vegetation in subalpine meadows in the Sierra Nevada Mountains is particularly vulnerable to lowering of groundwater levels because wet meadow vegetation is reliant upon shallow groundwater during the dry summer growing season. These ecosystems are especially vulnerable to channel incision as meadow aquifers are hydrologically connected to tributaries, and many have not yet recovered from previous anthropogenic influences. While instream restoration projects have become a common approach, lack of postrestoration monitoring and communication often result in a trial‐and‐error approach. In this study we demonstrate that preimplementation modeling of possible instream restoration solutions, chosen to raise stream stage and subsequently groundwater levels, is a useful tool for evaluating and comparing potential channel modifications. Modeling allows us to identify strategic locations and specific methods. Results show additional sediment depth and roughness on tributaries along with introduced woody debris (simulated by high roughness) on the Tuolumne River are the most effective means of raising stream stage. Results demonstrate that restoration efforts are most efficient in tributary streams. Managers and planners can more efficiently direct resources while minimizing the potential for negative impacts or failed restoration projects by modeling the possible effects of multiple restoration scenarios before implementation.     

Streamflow Depletion Modeling: Methods for an Adaptable and Conjunctive Water Management Decision Support Tool

Groundwater pumping depletes streamflow, which can have significant ecological impacts depending on the magnitude of depletion relative to environmental flow needs. Streamflow depletion estimates from groundwater pumping have been quantified using both analytical and numerical methods, but are not routinely compared to environmental flow needs or used in practical water management tools. Decision support tools that incorporate groundwater dynamics are becoming increasingly necessary for water managers as groundwater regulations become more important in environmental policy, particularly concerning the preservation of environmental flow needs. We develop and apply methods for a web‐based decision support tool for conjunctive groundwater and surface water management, demonstrated using a case study watershed in British Columbia, Canada. Open‐source data are analyzed with a combination of spatial algorithms and previously developed analytical models, such that the tool can be applied to other regions. Streamflow depletion estimates are calculated in four regions in the largely undeveloped Bulkley Valley, British Columbia. Our transparent methodology has geographic and data input flexibility which is a significant improvement on currently existing water management tool methods.     

Air‐Water Temperature Relationships in the Trout Streams of Southeastern Minnesota's Carbonate‐Sandstone Landscape

Carbonate‐sandstone geology in southeastern Minnesota creates a heterogeneous landscape of springs, seeps, and sinkholes that supply groundwater into streams. Air temperatures are effective predictors of water temperature in surface‐water dominated streams. However, no published work investigates the relationship between air and water temperatures in groundwater‐fed streams (GWFS) across watersheds. We used simple linear regressions to examine weekly air‐water temperature relationships for 40 GWFS in southeastern Minnesota. A 40‐stream, composite linear regression model has a slope of 0.38, an intercept of 6.63, and R² of 0.83. The regression models for GWFS have lower slopes and higher intercepts in comparison to surface‐water dominated streams. Regression models for streams with high R² values offer promise for use as predictive tools for future climate conditions. Climate change is expected to alter the thermal regime of groundwater‐fed systems, but will do so at a slower rate than surface‐water dominated systems. A regression model of intercept vs. slope can be used to identify streams for which water temperatures are more meteorologically than groundwater controlled, and thus more vulnerable to climate change. Such relationships can be used to guide restoration vs. management strategies to protect trout streams.     

Variations in Base‐Flow Nitrate Flux in a First‐Order Stream and Riparian Zone1

Abstract: Nonpoint source pollution, which contributes to contamination of surface waters, is difficult to control. Some pollutants, particularly nitrate (), are predominantly transmitted through ground water. Riparian buffer zones have the potential to remove contaminants from ground water and reduce the amount of that enters surface water. This is a justification for setting aside vegetated buffer strips along waterways. Many riparian zone hydrologic models assume uniform ground‐water flow through organic‐rich soil under reducing conditions, leading to effective removal of ground‐water prior to discharge into a stream. However, in a small first‐order stream in the mid‐Atlantic coastal plain, base‐flow generation was highly variable (spatially and temporally). Average base‐flow loads were greater in winter than summer, and higher during a wetter year than in dryer years. Specific sections of the stream consistently received greater amounts of high ground water than others. Areas within the riparian zone responsible for most of the exported from the watershed are termed “critical areas.” Over this 5‐year study, most of the exported during base flow originated from a critical area comprising less than 10% of the total riparian zone land area. Allocation of resources to address and improve mitigation function in critical areas should be a priority for continued riparian zone research.     

Riparian Ground‐Water Flow Patterns Using Flownet Analysis: Evapotranspiration‐Induced Upwelling and Implications for N Removal1

Abstract: Ground‐water flow paths constrain the extent of nitrogen (N) sinks in deep, stratified soils of riparian wetlands. We examined ground‐water flow paths at four forested riparian wetlands in deep, low gradient, stratified deposits subjected to Southern New England’s temperate, humid climate. Mid‐day piezometric heads were recorded during the high water table period in April/May and again in late November at one site. Coupling field data with a two‐dimensional steady‐state ground‐water flow model, flow paths and fluxes were derived to 3 m depths. April/May evapotranspiration (ET) dominated total outflux (44‐100%) while flux to the stream was <10% of total outflux. ET exerted upward ground‐water flux through shallow carbon‐rich soils, increasing opportunities for N transformations and diverting flow from the stream. Dormant season results showed a marked increase in flux to the stream (27% of the total flux). Riparian sites with deep water tables (naturally or because of increased urbanization or other hydrologic modifications) or shallow root zones may not generate ground‐water upwelling to meet evaporative demand, thereby increasing the risk of N movement to streams. As water managers balance issues of water quality with water quantity, they will be faced with decisions regarding riparian management. Further work towards refining our understanding of ET mediation of N and water flux at the catchment scale will serve to inform these decisions.     

Evaluating Flow Diversion Impacts to Groundwater‐Dependent Riparian Vegetation with Flow Alteration and Groundwater Model Analysis

An approach for assessing the potential ecologic response of groundwater‐dependent riparian vegetation to flow alteration is developed, focusing on change to groundwater. Groundwater requirements for riparian vegetation are reviewed in conjunction with flow alteration statistics. Where flow alteration coincides with groundwater‐related vegetation sensitivities, scenarios are developed for groundwater simulation. Groundwater depths and recession rates in the riparian zone are simulated for baseline and altered stream hydrographs, with changes to river stage and width represented with a transient, flow‐dependent boundary condition. Potential flow diversion from the Upper Gila River in New Mexico is examined. Statistical flow alteration analysis, applying prospective diversions to a 76‐year record of daily flow, shows that flows in the winter‐spring months and within the high‐pulse to small flood range are subject to greatest potential change. Groundwater simulation scenarios are developed for these flow conditions in representative dry, near‐average, and wet years. Differences in groundwater elevations, generally less than 0.25 m during the flow alteration period, dissipate rapidly following cessation of diversion. Relating groundwater depth, recession rates and range of fluctuations to riparian vegetation needs, we find adverse ecological response is not expected from groundwater impacts for the flow alteration examined.     

Estimating Lag to Peak between Rainfall and Peak Streamflow with a Mixed‐Effects Model

We test the use of a mixed‐effects model for estimating lag to peak for small basins in Maine (drainage areas from 0.8 to 78 km²). Lag to peak is defined as the time between the center of volume of the excess rainfall during a storm event and the resulting peak streamflow. A mixed‐effects model allows for multiple observations at sites without violating model assumptions inherent in traditional ordinary least squares models, which assume each observation is independent. The mixed model includes basin drainage area and maximum 15‐min rainfall depth for individual storms as explanatory features. Based on a remove‐one‐site cross‐validation analysis, the prediction errors of this model ranged from ?42% to +73%. The mixed model substantially outperformed three published models for lag to peak and one published model for centroid lag for estimating lag to peak for small basins in Maine. Lag to peak estimates are a key input to rainfall–runoff models used to design hydraulic infrastructure. The improved accuracy and consistency with model assumptions indicates that mixed models may provide increased data utilization that could enhance models and estimates of lag to peak in other regions.     

Evaluation of Rio Grande Management Alternatives Using a Surface‐Water/Ground‐Water Model1

Abstract: Previous investigations observed significant seepage losses from the Rio Grande to the shallow aquifer between Socorro and San Antonio, New Mexico. High‐resolution telescopic modeling was used along a 10‐km reach of the Rio Grande and associated drains and canals to evaluate several management alternatives aimed at improving river conveyance efficiency. Observed data consisted of ground‐water and surface‐water elevations, seepage rates along the Rio Grande and associated canals and drains, and borehole geology. Model calibration was achieved by adjusting hydraulic conductivity and specific storage until the output matched observed data. Sensitivity analyses indicated that the system was responsive to changes in hydrogeologic properties, especially when such alterations increased vertical connectivity between layers. The calibrated model predicted that removal of the low flow conveyance channel, a major channel draining the valley, would not only decrease river seepage by 67%, but also decrease total flow through the reach by 75%. The decreased flow through the reach would result in increased water logging and an average increase in ground‐water elevations of 1.21 meter. Simulations of the system with reduced riparian evapotranspiration rates or a relocated river channel also predicted decreased river seepage, but to a much lesser degree.     

Evaluation of Ground‐Water and Surface‐Water Exchanges Using Streamflow Difference Analyses1

Abstract: In the karstic lower Flint River Basin, limestone fracturing, jointing, and subsequent dissolution have resulted in the development of extensive secondary permeability and created a system of major conduits that facilitate the exchange of water between the Upper Floridan aquifer and Flint River. Historical streamflow data from U.S. Geological Survey gaging stations located in Albany and Newton, Georgia, were used to quantify ground‐water and surface‐water exchanges within a 55.3 km section of the Flint River. Using data from 2001, we compared estimates of ground‐water flux using a time adjustment method to a water balance equation and found that these independent approaches yielded similar results. The associated error was relatively large during high streamflow when unsteady conditions prevail, but much lower during droughts. Flow reversals were identified by negative streamflow differences and verified with in situ data from temperature sensors placed inside large spring conduits. Long‐term (13 years) analysis showed negative streamflow differentials (i.e., a losing stream condition) coincided with high river stages and indicated that streamflow intrusion into the aquifer could potentially exceed 150 m³/s. Although frequent negative flow differentials were evident, the Flint River was typically a gaining stream and showed a large net increase in flow between the two gages when examined over the period 1989‐2003. Ground‐water contributions to this stream section averaged 2‐42 m³/s with a mean of 13 m³/s. The highest rate of ground‐water discharge to the Flint River occurred during the spring when regional ground‐water levels peaked following heavy winter and spring rains and corresponding rates of evapotranspiration were low. During periods of extreme drought, ground‐water contributions to the Flint River declined.     

Incorporating Variable Source Area Hydrology into a Spatially Distributed Direct Runoff Model1

Buchanan, Brian, Zachary M. Easton, Rebecca Schneider, and M. Todd Walter, 2011. Incorporating Variable Source Area Hydrology Into a Spatially Distributed Direct Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(1): 43‐60. DOI: 10.1111/j.1752‐1688.2011.00594.x Abstract: Few hydrologic models simulate both variable source area (VSA) hydrology, and runoff‐routing at high enough spatial resolutions to capture fine‐scale hydrologic pathways connecting VSA to the stream network. This paper describes a geographic information system‐based operational model that simulates the spatio‐temporal dynamics of VSA runoff generation and distributed runoff‐routing, including through complex artificial drainage networks. The model combines the Natural Resource Conservation Service’s Curve Number (CN) equation for estimating storm runoff with the topographic index concept for predicting the locations of VSA and a runoff‐routing algorithm into a new spatially distributed direct hydrograph (SDDH) model (SDDH‐VSA). Using a small agricultural watershed in central New York, SDDH‐VSA results were compared to those from a SDDH model using the traditional land use assumptions for the CN (SDDH‐CN). The SDDH‐VSA model generally agreed better with observed discharge than the SDDH‐CN model (average, Nash‐Sutcliffe efficiency of 0.69 vs. 0.58, respectively) and resulted in more realistic spatial patterns of runoff‐generating areas. The SDDH approach did not correctly capture the timing of runoff from small storms in dry periods. Despite this type of limitation, SDDH‐VSA extends the applicability of the SDDH technique to VSA conditions, providing a basis for new tools to help identify critical management areas and assess water quality risks due to landscape alterations.     

The Role of Ground Water in Generating Streamflow in Headwater Areas and in Maintaining Base Flow1

Abstract: The volume and sustainability of streamflow from headwaters to downstream reaches commonly depend on contributions from ground water. Streams that begin in extensive aquifers generally have a stable point of origin and substantial discharge in their headwaters. In contrast, streams that begin as discharge from rocks or sediments having low permeability have a point of origin that moves up and down the channel seasonally, have small incipient discharge, and commonly go dry. Nearly all streams need to have some contribution from ground water in order to provide reliable habitat for aquatic organisms. Natural processes and human activities can have a substantial effect on the flow of streams between their headwaters and downstream reaches. Streams lose water to ground water when and where their head is higher than the contiguous water table. Although very common in arid regions, loss of stream water to ground water also is relatively common in humid regions. Evaporation, as well as transpiration from riparian vegetation, causing ground‐water levels to decline also can cause loss of stream water. Human withdrawal of ground water commonly causes streamflow to decline, and in some regions has caused streams to cease flowing.     

Development and Operational Testing of a Super‐Ensemble Artificial Intelligence Flood‐Forecast Model for a Pacific Northwest River

Coastal catchments in British Columbia, Canada, experience a complex mixture of rainfall‐ and snowmelt‐driven contributions to flood events. Few operational flood‐forecast models are available in the region. Here, we integrated a number of proven technologies in a novel way to produce a super‐ensemble forecast system for the Englishman River, a flood‐prone stream on Vancouver Island. This three‐day‐ahead modeling system utilizes up to 42 numerical weather prediction model outputs from the North American Ensemble Forecast System, combined with six artificial neural network‐based streamflow models representing various slightly different system conceptualizations, all of which were trained exclusively on historical high‐flow data. As such, the system combines relatively low model development times and costs with the generation of fully probabilistic forecasts reflecting uncertainty in the simulation of both atmospheric and terrestrial hydrologic dynamics. Results from operational testing by British Columbia's flood forecasting agency during the 2013‐2014 storm season suggest that the prediction system is operationally useful and robust.     

Estimating Evapotranspiration and Seepage for a Sinkhole Wetland From Diurnal Surface‐Water Cycles1

Abstract: This study used measured diurnal surface‐water cycles to estimate daily evapotranspiration (ET) and seepage for a seasonally flooded sinkhole wetland. Diurnal surface‐water cycles were classified into five categories based on the relationship between the surface‐water body and the surrounding ground‐water system (i.e., recharge/discharge). Only one class of diurnal cycles was found to be suitable for application of this method. This subset of diurnal cycles was used to estimate ET and seepage and the relative importance of each transfer process to the overall water budget. The method has limited utility for wetlands with erratic hydrologic regimes (e.g., wetlands in urban environments). This is due to violation of the critical assumption that the inflow/outflow rate remains constant throughout the day. For application to surface‐water systems, the method is typically applied with an assumed specific yield of 1.0. This assumption was found to be invalid for application to surface‐water systems with a noncylindrical pond geometry. An overestimation of ET by as much as 60% was found to occur under conditions of low pond stage and high water loss. The results demonstrate the high ET rates that can occur in isolated wetlands due to contrasting roughness and moisture conditions (oasis and clothesline effects). Estimated ET rates ranged from 4.1 to 18.7 mm/day during the growing season. Despite these large ET rates, seepage (recharge) was found to be the dominant water loss mechanism for the wetland.     

Salt Flushing,Salt Storage,and Controls on Selenium and Uranium: A 31‐Year Mass‐Balance Analysis of an Irrigated,Semiarid Valley

Salinity, selenium, and uranium pose water‐quality challenges for the Arkansas River in southeastern Colorado and other rivers that support irrigation in semiarid regions. This study used 31 years of continuous discharge and specific conductance (SC) monitoring data to assess interannual patterns in water quality using mass balance on a 120‐km reach of river. Discrete sampling data were used to link the SC records to salinity, selenium, and uranium. Several important patterns emerged. Consumptive use reduced discharge by a median value of 33% and drove corresponding increases in salinity and uranium concentrations. Increased water availability for irrigation from rainfall and upstream snowpack in 1995–1999 flushed additional salinity and uranium into the river in 1996–2000; average annual total dissolved solids (salinity) concentrations increased 25%, and loads increased 131%. Smaller flushing events have occurred, sometimes lagging an increase in water availability by about one year. The pattern indicates flushing of salts temporarily stored, evaporatively concentrated, or of geologic origin. Mobilization of selenium from the reach was minor compared to salinity and uranium, and net selenium removal from the river was suggested in some years. Several processes related to irrigation could be removing selenium. The results provide context for efforts to improve water quality in the Arkansas River and rivers in other semiarid regions.