EVALUATION OF A STREAM AQUIFER ANALYSIS TEST USING ANALYTICAL SOLUTIONS AND FIELD DATA1

ABSTRACT: Considerable advancements have been made in the development of analytical solutions for predicting the effects of pumping wells on adjacent streams and rivers. However, these solutions have not been sufficiently evaluated against field data. The objective of this research is to evaluate the predictive performance of recently proposed analytical solutions for unsteady stream depletion using field data collected during a stream/aquifer analysis test at the Tamarack State Wildlife Area in eastern Colorado. Two primary stream/aquifer interactions exist at the Tamarack site: (1) between the South Platte River and the alluvial aquifer and (2) between a backwater stream and the alluvial aquifer. A pumping test is performed next to the backwater stream channel. Drawdown measured in observation wells is matched to predictions by recently proposed analytical solutions to derive estimates of aquifer and streambed parameters. These estimates are compared to documented aquifer properties and field measured streambed conductivity. The analytical solutions are capable of estimating reasonable values of both aquifer and streambed parameters with one solution capable of simultaneously estimating delayed aquifer yield and stream flow recharge. However, for long term water management, it is reasonable to use simplified analytical solutions not concerned with early‐time delayed yield effects. For this site, changes in the water level in the stream during the test and a varying water level profile at the beginning of the pumping test influence the application of the analytical solutions.     

STREAMFLOW DEPLETION: MODELING OF REDUCED BASEFLOW ANI INDUCED STREAM INFILTRATION FROM SEASONALLY PUMPED WELLS1

ABSTRACT: Numerical modeling techniques are used to analyze streamflow depletion for stream‐aquifer systems with baseflow. The analyses calculated two flow components generated by a pumping well located at a given distance from a river that is hydraulically connected to an unconfined aquifer. The two components are induced stream infiltration and reduced baseflow; both contribute to total streamflow depletion. Simulation results suggest that the induced infiltration, the volume of water discharged from the stream to the aquifer, has a shorter term impact on streamflow, while the reduced baseflow curves show a longer term effect. The peak impacts of the two hydrologic processes on streamflow occur separately. The separate analysis helps in understanding the hydrologic interactions between stream and aquifer. Practically, it provides useful information about contaminant transport from stream to aquifer when water quality is a concern, and for areas where water quantity is an issue, the separate analysis offers additional information to the development of water resource management plan.     

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.     

Effects of Rainfall and Ground‐Water Pumping on Streamflow in Mākaha,O’ahu,Hawai’i1

Abstract: Land‐use/land‐cover changes in Mākaha valley have included the development of agriculture, residential dwellings, golf courses, potable water supply facilities, and the introduction of alien species. The impact of these changes on surface water and ground water resources in the valley is of concern. In this study, streamflow, rainfall, and ground‐water pumping data for the upper part of the Mākaha valley watershed were evaluated to identify corresponding trends and relationships. The results of this study indicate that streamflow declined during the ground‐water pumping period. Mean and median annual streamflow have declined by 42% (135 mm) and 56% (175 mm), respectively, and the mean number of dry stream days per year has increased from 8 to 125. Rainfall across the study area appears to have also declined though it is not clear whether the reduction in rainfall is responsible for all or part of the observed streamflow decline. Mean annual rainfall at one location in the study area declined by 14% (179 mm) and increased by 2% (48 mm) at a second location. Further study is needed to assess the effect of ground‐water pumping and to characterize the hydrologic cycle with respect to rainfall, infiltration, ground‐water recharge and flow in the study area, and stream base flow and storm flow.     

WATER LOSSES ALONG A REACH OF THE PECOS RIVER IN NEW MEXICO1

ABSTRACT: A reach of the Pecos River, located in eastern New Mexico, was examined to evaluate losses of river flows due to evaporation, seepage, and transpiration. An accurate assessment of the water losses along this reach is critical for determining how water rights are adjudicated for water users in the Pecos basin and interstate compact accounting. Water losses significantly impact flows through critical habitat for species protected under the Endangered Species Act. Daily losses of river flows were analyzed for the study reach that extends from immediately below the Pecos River confluence with Taiban Creek to the United States Geological Survey (USGS) gage near Acme. The analysis was completed with consideration for other processes including flood wave travel times and attenuation along with stream bank storage and returns. The analysis was completed using daily stream flow data from USGS gages located along the study reach. Empirical seasonal functions were developed to relate flow loss to the flow rate in the river. The functions were ultimately developed to provide a method for comparing the effects of different river flows on the available water supply.     

Conjunctive Water Use in Confined Basalt Aquifers: An Evaluation Using Geochemistry,a Numerical Model,and Historical Water Level

As withdrawals from deep compartmentalized aquifers increasingly exceed recharge throughout the western United States, conjunctive water use management alternatives have become an applied research priority. This study highlights both details and limitations of the role of irrigation canal seepage as groundwater recharge, revealing the regional limitations of canal seepage as a dependable source of recharge in overdrawn aquifers. A suite of geochemical indicators were used together with a numerical model to evaluate current and future management scenarios focused on recharge derived from seepage from a region‐wide irrigation canal system. Twenty‐five years of static groundwater level data were used to relate spatial trends determined using geochemistry and groundwater modeling with “on‐the‐ground” management practices, which vary based on acreage, crop, and irrigation scheduling. Increasing groundwater age determined using isotope analysis, and declines in potentiometric heads, each correlate with increasing distance from the canal reaches. Predictive modeling indicates that if pumping is gradually reduced, as has been suggested by management agencies, that recharge from canal seepage will be negligible by 2035 due to regional groundwater through‐flow and the pattern of potentiometric head recovery. Unfortunately, historic hydrographs suggest that under current groundwater development conditions most wells are not sustainable, irrespective of proximity to the canal.     

Classification of Ephemeral,Intermittent, and Perennial Stream Reaches Using a TOPMODEL‐Based Approach

Whether a waterway is temporary or permanent influences regulatory protection guidelines, however, classification can be subjective due to a combination of factors, including time of year, antecedent moisture conditions, and previous experience of the field investigator. Our objective was to develop a standardized protocol using publically available spatial information to classify ephemeral, intermittent, and perennial streams. Our hypothesis was that field observations of flow along the stream channel could be compared to results from a hydrologic model, providing an objective method of how these stream reaches can be identified. Flow‐state sensors were placed at ephemeral, intermittent, and perennial stream reaches from May to December 2011 in the Appalachian coal basin of eastern Kentucky. This observed flow record was then used to calibrate the simulated saturation deficit in each channel reach based on the topographic wetness index used by TOPMODEL. Saturation deficit values were categorized as flow or no‐flow days, and the simulated record of streamflow was compared to the observed record. The hydrologic model was more accurate for simulating flow during the spring and fall seasons. However, the model effectively identified stream reaches as intermittent and perennial in each of the two basins.     

Modeling the Potential Impact of Seasonal and Inactive Multi‐Aquifer Wells on Contaminant Movement to Public Water‐Supply Wells1

Johnson, R.L., B.R. Clark, M.K. Landon, L.J. Kauffman, and S.M. Eberts, 2011. Modeling the Potential Impact of Seasonal and Inactive Multi‐Aquifer Wells on Contaminant Movement to Public Water‐Supply Wells. Journal of the American Water Resources Association (JAWRA) 47(3):588‐596. DOI: 10.1111/j.1752‐1688.2011.00526.x Abstract: Wells screened across multiple aquifers can provide pathways for the movement of surprisingly large volumes of groundwater to confined aquifers used for public water supply (PWS). Using a simple numerical model, we examine the impact of several pumping scenarios on leakage from an unconfined aquifer to a confined aquifer and conclude that a single inactive multi‐aquifer well can contribute nearly 10% of total PWS well flow over a wide range of pumping rates. This leakage can occur even when the multi‐aquifer well is more than a kilometer from the PWS well. The contribution from multi‐aquifer wells may be greater under conditions where seasonal pumping (e.g., irrigation) creates large, widespread downward hydraulic gradients between aquifers. Under those conditions, water can continue to leak down a multi‐aquifer well from an unconfined aquifer to a confined aquifer even when those multi‐aquifer wells are actively pumped. An important implication is that, if an unconfined aquifer is contaminated, multi‐aquifer wells can increase the vulnerability of a confined‐aquifer PWS well.     

Flow‐Biology Relationships Based on Fish Habitat Guilds in North Carolina

The health of freshwater biota is dependent on streamflow, yet identification of the flow regimes required to maintain ecological integrity remains challenging to states in the United States seeking to establish ecological flows. We tested the relationship between decreases in streamflow and Shannon‐Weaver diversity index of fish species for four flow‐based habitat guilds: riffle, riffle‐run, pool‐run, and pool in North Carolina. We found species that prefer shallow habitats, such as riffles and riffle‐runs were the most sensitive to decreases in streamflow; whereas no significant relationships were found for pool or pool‐run species. The sensitivity to decreases in streamflow was greatest during summer and fall, when streams are naturally lower. When all fish habitat guilds were included in the assessment of flow‐biology relationships, there were no significant relationships to decreases in streamflow. As the sensitivity of fish to reductions in streamflow is not constant across habitat guilds, combining all fish species together for flow‐biology analyses may greatly underestimate the response of fish species to decreases in flow and should be acknowledged when establishing ecological flows.     

EVALUATION OF AN INDUCED INFILTRATION MODEL AS APPLIED TO GLACIAL AQUIFER SYSTEMS1

ABSTRACT: Numerical models were used to examine the limitations of the assumptions used in an analytical induced infiltration model. The assumptions tested included negligible streambed effects, negligible areal recharge, two-dimensional ground water flow, fully penetrating rivers and wells, and constant surface water stage. For situations that deviate from the underlying assumptions, the analytical model becomes a less reliable predictor of induced infiltration. The numerical experiments show that streambed effects cannot be neglected if the streambed conductivity is more than one order of magnitude lower than the aquifer hydraulic conductivity. Areal recharge cannot be neglected if the ground water basin receives more than 5 in/yr of areal recharge. Three-dimensional flow effects due to well partial penetration cannot be neglected if the ratio of horizontal hydraulic conductivity to vertical hydraulic conductivity (K_h/K_u) is greater than 10. Surface water elevation fluctuations can significantly influence the induced infiltration rate, depending on the degree of fluctuations and the ground water hydraulic gradient.     

Estimating Hydraulic Properties of Volcanic Aquifers Using Constant‐Rate and Variable‐Rate Aquifer Tests1

Abstract: In recent years the ground‐water demand of the population of the island of Maui, Hawaii, has significantly increased. To ensure prudent management of the ground‐water resources, an improved understanding of ground‐water flow systems is needed. At present, large‐scale estimations of aquifer properties are lacking for Maui. Seven analytical methods using constant‐rate and variable‐rate withdrawals for single wells provide an estimate of hydraulic conductivity and transmissivity for 103 wells in central Maui. Methods based on constant‐rate tests, although not widely used on Maui, offer reasonable estimates. Step‐drawdown tests, which are more abundantly used than other tests, provide similar estimates as constant‐rate tests. A numerical model validates the suitability of analytical solutions for step‐drawdown tests and additionally provides an estimate of storage parameters. The results show that hydraulic conductivity is log‐normally distributed and that for dike‐free volcanic rocks it ranges over several orders of magnitude from 1 to 2,500 m/d. The arithmetic mean, geometric mean, and median values of hydraulic conductivity are respectively 520, 280, and 370 m/d for basalt and 80, 50, and 30 m/d for sediment. A geostatistical approach using ordinary kriging yields a prediction of hydraulic conductivity on a larger scale. Overall, the results are in agreement with values published for other Hawaiian islands.     

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.     

Clogging of an Effluent Dominated Semiarid River: A Conceptual Model of Stream‐Aquifer Interactions1

Abstract: Water managers in arid and semiarid regions increasingly view treated wastewater (effluent) as an important water resource. Artificial recharge basins allow effluent to seep into the ground relieving stressed aquifers, however these basins frequently clog due to physical, chemical, and biological processes. Likewise effluent is increasingly used to maintain perennial base flow for dry streambeds, however, little is known about the impact of effluent on streambed hydraulic conductivity and stream‐aquifer interactions. We address this issue by investigating: if a clogging layer forms, how the formation of a clogging layer alters stream‐aquifer connections, and what hydrologic factors control the formation and removal of clogging layers. We focused on the Upper Santa Cruz River, Arizona where effluent from the Nogales International Waste Water Treatment Plant sustains perennial flow. Monthly sampling, along a 30 km river reach, was done with two foci: physical streambed transformations and water source identification using chemical composition. Historical dataset were included to provide a larger context for the work. Results show that localized clogging occurs in the Upper Santa Cruz River. The clogging layers perch the stream and shallow streambed causing desaturation below the streambed. With these results, a conceptual model of clogging is established in the context of a semiarid hydrologic cycle: formation during the hot premonsoon months when flow is nearly constant and removal by large flood flows (>10 m³/s) during the monsoon season. However, if the intensity of flooding during the semiarid hydrologic cycle is lessened, the dependent riparian area can experience a die off. This conceptual model leads us to the conclusion that effluent dominated riparian systems are inherently unstable due to the clogging process. Further understanding of this process could lead to improved ecosystem restoration and management.     

Allocation of Streamflow Depletion Impacts under Nonlinear Conditions

A method is proposed for the equitable allocation of impacts of groundwater pumping on streamflow. The method is intended for cases in which the pumping activity of multiple entities has impacts on streamflow and these impacts are computed by perturbation. It is shown that when the response of streamflow to pumping is nonlinear, simple methods for impact calculation can fail. The proposed method is developed for the case when there are four entities that impact streamflow. The method relies on the calculation of impacts by perturbation of the simulation model from different base pumping levels. When four entities are evaluated, 16 runs of the simulation model are required. It is shown the proposed method produces estimated impacts for each individual entity that are equitable because they meet the requirement that the impacts of each entity sum to the total impacts of all entities acting together and the impacts attributed to each entity do not depend on the order of calculation. A brief example demonstrates the approach.     

ANALYSIS OF POND SEEPAGE FOR APPLICATION IN FISHERIES AND AQUACULTURE1

ABSTRACT: The purpose of this research was to examine through modeling and experimentation if seepage out of a pond through stratified soil can be predicted, and effectively collected and managed to augment streamflow during a low precipitation period extending three months or more. The 55 m² experimental pond with sandy/loamy banks was excavated to hardpan, and its bottom was approximately 0.7 above the water table. Output from a mathematical model containing both bottom and bank seepage elements agreed with experimental data, and showed that as compared to bottom seepage, the bank seepage contributed approximately 25 percent of the total seepage. Seepage collection (as measured from a circumscribing ditch) linearly varied with stage (r² < 0.99). There was an 8 to 22 percent over‐collection at the lower pond stages, and a 9 to 45 percent under‐collection at the highest stage. As an example of its utility, the model was applied to estimate the pond size and shape needed to supply a hypothetical stream and maintain fish stocks during a three‐month low‐precipitation period. Future work will focus on nutrient transport and removal.     

Revealing the Diversity of Natural Hydrologic Regimes in California with Relevance for Environmental Flows Applications

Alterations to flow regimes for water management objectives have degraded river ecosystems worldwide. These alterations are particularly profound in Mediterranean climate regions such as California with strong climatic variability and riverine species highly adapted to the resulting flooding and drought disturbances. However, defining environmental flow targets for Mediterranean rivers is complicated by extreme hydrologic variability and often intensive water management legacies. Improved understanding of the diversity of natural streamflow patterns and their spatial arrangement across Mediterranean regions is needed to support the future development of effective flow targets at appropriate scales for management applications with minimal resource and data requirements. Our study addresses this need through the development of a spatially explicit reach‐scale hydrologic classification for California. Dominant hydrologic regimes and their physio‐climatic controls are revealed, using available unimpaired and naturalized streamflow time‐series and generally publicly available geospatial datasets. This methodology identifies eight natural flow classes representing distinct flow sources, hydrologic characteristics, and catchment controls over rainfall‐runoff response. The study provides a broad‐scale hydrologic framework upon which flow‐ecology relationships could subsequently be established towards reach‐scale environmental flows applications in a complex, highly altered Mediterranean region.     

Automated Algorithms for Heuristic Base‐Flow Separation1

Abstract: The subjective nature of graphical base‐flow separation combined with the many applications of base‐flow time series derived from continuous streamflow data, motivates the development and application of automated algorithms for heuristic base‐flow separation. Base‐flow time series derived from gauged streamflow support diverse applications in engineering hydrology, catchment analysis, hydrogeologic investigations, regional low‐flow analysis, and recharge estimation. Whether based on graphical procedures for recession analysis or analytical expressions derived from fundamental equations of ground‐water flow, the variety of base‐flow separation algorithms belies the array of base‐flow definitions and interpretations that variously refer to dominant process, source, flow path, and characteristic response time. Algorithms that are invariant in their consistent – though heuristic – characterization of base‐flow response are particularly useful for interbasin comparisons of low‐flow characteristics and hydrologic regionalization. More adaptable algorithms provide application‐specific flexibility in allocating flow components like interflow to either quickflow or slowflow. Four widely used algorithms that produce consistent base‐flow time series using only gauged streamflow records are compared and contrasted with a complementary heuristic algorithm that incorporates hydrologic judgment explicitly, through manual parameterization. The utility of these inherently subjective algorithms is illustrated through a simple example of flow phase separation in a two‐component end‐member mixing model of dissolved chlorides in the Cuyahoga River.     

PREDICTION OF MASTER RECESSION CURVES AND BASEFLOW RECESSIONS IN THE LUQUILLO MOUNTAINS OF PUERTO RICO1

ABSTRACT: An analysis of hydrograph recessions and rainfall data was performed to estimate the recession constants for two watersheds in the Luquillo mountains of Puerto Rico. To account for seasonal rainfall patterns, the data were grouped into dry and wet seasons. Sets of three Master Recession Curves (MRC) per season for each watershed were developed: one using the Matching Strip Method (MS) and two using variations of the Correlation Method (CM). These variations were the envelope line (CME) and the least squares regression (CMR). Other regression based analytical expressions that consider the streamflow recession as an autore‐gressive or an integrated moving average process were also applied. The regression based methods performed consistently better than the graphical ones and they proved to be faster, easier, and less subjective. The recession constants from these methods were then used to estimate the time it would take the streamflow to reach the critical Q₉₉ flow duration. Based on this study, once the streamflow reaches Q₉₀, water managers have 6 to 12 days warning before streamflow reaches critical levels.     

IDENTIFICATION OF CHANGES IN STREAMFLOW CHARACTERISTICS1

ABSTRACT: A set of procedures for identifying changes in selected streamflow characteristics at sites having long‐term continuous streamflow records is illustrated by using streamflow data from the Waccamaw River at Freeland, North Carolina for the 55‐year period of 1940–1994. Data were evaluated and compared to streamflow in the adjacent Lumber River Basin to determine if changes in streamflow characteristics in the Waccamaw River were localized and possibly the result of some human activity, or consistent with regional variations. Following 1963, droughts in the Waccamaw Basin seem to have been less severe than in the Lumber Basin, and the annual one‐, seven‐, and 30‐day low flows exhibited a slightly increasing trend in the Waccamaw River. Mean daily flows in the Waccamaw River at the 90 percent exceedance level (low flows) during 1985–194, a relatively dry period, were very nearly equal to flows at the same exceedance level for 1970–1979, which represents the 10‐year period between 1940 and 1994 with the highest flows. Prior to the 1980s, flows per unit drainage area in the Waccamaw Basin were generally less than those in the Lumber Basin, but after 1980, the opposite was true. The ratio of base flow to runoff in the Waccamaw River may have changed relative to that in the Lumber River in the late 1970s. There was greater variability in Waccamaw River streamflow than in Lumber River flow, and flow variability in the Waccamaw River may have increased slightly during 1985–1994.     

Impact of Length of Dataset on Streamflow Calibration Parameters and Performance of APEX Model

Due to resource constraints, long‐term monitoring data for calibration and validation of hydrologic and water quality models are rare. As a result, most models are calibrated and, if possible, validated using limited measured data. However, little research has been done to determine the impact of length of available calibration data on model parameterization and performance. The main objective of this study was to evaluate the impact of length of calibration data (LCD) on parameterization and performance of the Agricultural Policy Environmental eXtender model for predicting daily, monthly, and annual streamflow. Long‐term (1984‐2015) measured daily streamflow data from Rock Creek watershed, an agricultural watershed in northern Ohio, were used for this study. Data were divided into five Short (5‐year), two Medium (15‐year), and one Long (25‐year) streamflow calibration data scenarios. All LCD scenarios were calibrated and validated at three time steps: daily, monthly, and annual. Results showed LCD affected the ability of the model to accurately capture temporal variability in simulated streamflow. However, overall average streamflow, water budgets, and crop yields were simulated reasonably well for all LCD scenarios.