SENSITIVITY OF GROUNDWATER FLOW MODELS TO VERTICAL VARIABILITY OF AQUIFER CONSTANTS1

ABSTRACT: The Ogallala aquifer in the Oklahoma Panhandle is in need of better management because of increased groundwater demand which has caused declines in static water levels at an alarming rate. A groundwater management computer model was developed for the Ogallala aquifer in the Texas Panhandle and treats the aquifer as a homogeneous system. In this study, the computer model has been modified in order to evaluate the effects of vertical layering on semi-static water level changes which occur during the dewatering of a single unconfined aquifer. The modified model was applied to a study area near Guymon, Oklahoma, using both the homogeneous and the multilayered cases. The aquifer is characterized by a saturated thickness of 400 feet. The accumulated drawdown values of the homogeneous and the multilayered cases demonstrate that an average difference of approximately 22% of the original saturated thickness occurs between the two cases before the base of the aquifer is encountered. Approximately 25% more time is required to dewater the layered aquifer. Thus, vertical variations of lithology in an aquifer such as the Ogallala should be considered when prediction is made relative to groundwater management.     

Assessment of Groundwater Depletion and Implications for Management in the Denver Basin Aquifer System

The Denver Basin Aquifer System (DBAS) is a critical groundwater resource along the Colorado Front Range. Groundwater depletion has been documented over the past few decades due to the increased water use among users, presenting long‐term sustainability challenges. A spatiotemporal geostatistical analysis is used to estimate potentiometric surfaces and evaluate groundwater storage changes between 1990 and 2016 in each of the four DBAS aquifers. Several key depletion patterns and spatial water‐level changes emerge in this work. Hydraulic head changes are the largest in the west‐central side of the DBAS and have decreased in some areas by up to 180 m since 1990, while areas to the northwest show increases in hydraulic head by over 30.5 m. The Denver and Arapahoe aquifers show the largest groundwater storage losses, with the highest rates occurring in the 2000s. The results highlight uncertainty in the volumetric predictions under various storage coefficient calculations and emphasize the importance of representative aquifer characterization. The observed groundwater storage depletions are due to a combination of factors, which include population growth increasing the demand for water, variable precipitation, and drought influencing recharge, and increased groundwater pumping. The methods applied in this study are transferable to other groundwater systems and provide a framework that can help assess groundwater depletion and inform management decisions at other locations.     

Mason Valley Groundwater Model: Linking Surface Water and Groundwater in the Walker River Basin,Nevada1

Carroll, Rosemary W.H., Greg Pohll, David McGraw, Chris Garner, Anna Knust, Doug Boyle, Tim Minor, Scott Bassett, and Karl Pohlmann, 2010. Mason Valley Groundwater Model: Linking Surface Water and Groundwater in the Walker River Basin, Nevada. Journal of the American Water Resources Association (JAWRA) 46(3):554-573. DOI: 10.1111/j.1752-1688.2010.00434.x Abstract: An integrated surface water and groundwater model of Mason Valley, Nevada is constructed to replicate the movement of water throughout the different components of the demand side of water resources in the Walker River system. The Mason Valley groundwater surface water model (MVGSM) couples the river/drain network with agricultural demand areas and the groundwater system using MODFLOW, MODFLOW’s streamflow routing package, as well as a surface water linking algorithm developed for the project. The MVGSM is capable of simulating complex feedback mechanisms between the groundwater and surface water system that is not dependent on linearity among the related variables. The spatial scale captures important hydrologic components while the monthly stress periods allow for seasonal evaluation. A simulation spanning an 11-year record shows the methodology is robust under diverse climatic conditions. The basin-wide modeling approach predicts a river system generally gaining during the summer irrigation period but losing during winter months and extended periods of drought. River losses to the groundwater system approach 25% of the river’s annual budget. Reducing diversions to hydrologic response units will increase river flows exiting the model domain, but also has the potential to increase losses from the river to groundwater storage.     

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.     

The Transboundary Nature of the Allende–Piedras Negras Aquifer Using a Numerical Model Approach

The Allende–Piedras Negras (APN) aquifer is located between the states of Texas (United States [U.S.]) and Coahuila (Mexico). The Rio Grande crosses the aquifer, acting as a natural and political divide between the countries. However, it remains unclear whether the APN aquifer can be considered a truly transboundary aquifer flow system, which would potentially require joint management by two different administrative jurisdictions. The main purpose of this study was to evaluate the transboundary nature of this aquifer. This was achieved by developing a detailed hydrogeological model to analyze the direction of volumetric fluxes within the APN aquifer using Visual MODFLOW. The model simulated a spatially averaged cumulative drawdown of 0.76 m for the entire aquifer over an 18‐year modeling period (2000–2017). The flow convergence zone, previously located below the Rio Grande, has shifted to the U.S. side in most locations, driven by higher pumping rates of the wells located near the river. This shift of the convergence zone from one country to the other means that groundwater recharge from one side flows underneath the river to the other side. This qualifies the APN aquifer as a “transboundary groundwater flow system.” The procedure followed in this study may be applied to other aquifers that straddle the U.S.–Mexico border and may motivate future modeling studies on other poorly studied transboundary aquifers around the world and thereby enable bi‐national aquifer management.     

Using a Mass Balance Model to Evaluate Groundwater Budget of Seawater‐Intruded Island Aquifers1

Jang, Cheng‐Shin, Chen‐Wuing Liu, Shih‐Kai Chen, and Wen‐Sheng Lin, 2011. Using a Mass Balance Model to Evaluate Groundwater Budget of Seawater‐Intruded Island Aquifers. Journal of the American Water Resources Association (JAWRA) 48(1): 61‐73. DOI: 10.1111/j.1752‐1688.2011.00593.x Abstract: The study developed a mass balance model to evaluate the groundwater budget of seawater‐intruded island aquifers using limited available data. The Penghu islands were selected as a study area. As sparse observed data were available in the islands, methods of combining water and chloride balances were used to determine the amounts of groundwater pumping, seawater intrusion, aquifer storages, and safe yields in the shallow and deep aquifers. The groundwater budget shows that seawater intrusion to freshwater aquifers was 1.38 × 10⁶ and 0.29 × 10⁶ m³/year in the shallow and deep aquifers, respectively, indicating that the seawater intrusion is severe in the both aquifers. The safe yield of the shallow aquifer was 14.56 × 10⁶ m³/year in 2005 which was four times higher than that of the deep aquifer (3.70 × 10⁶ m³/year). However, the annual pumping amounts in the shallow and deep aquifers were 4.77 × 10⁶ and 3.63 × 10⁶ m³/year, respectively. Although the safe yield of the shallow aquifer is enough for all water resources demands, only 55% of exploitation amount was extracted from the shallow aquifer due to its poor water quality. Groundwater exploitation in the deep aquifer should be significantly reduced and regulated by a dynamic management of pumping scheme because the annual pumping amounts are close to the safe yield and seawater intrusion occurs continually. Additionally, to alleviate further aquifer salination, at least half of the current annual groundwater abstraction should be reduced.     

ASSESSMENT OF HYDRAULIC PROPERTIES IN AN UNCONFINED ALLUVIAL AQUIFER NEAR GRAND ISLAND,NEBRASKA1

ABSTRACT: A method is presented to analyze time-drawdown data from one or more observation wells for the calculation of four hydraulic parameters for unconfined aquifers: vertical hydraulic conductivity, horizontal hydraulic conductivity, storage coefficient, and specific yield. The hydraulic parameter results are further analyzed for reliability and the possible ranges of the actual parameter values. After verification using a theoretical example, the method was used to analyze pumping test data from 22 observation wells in an unconfined alluvial aquifer near Grand Island, Nebraska. Results indicate that this method can be used to efficiently calculate the four hydraulic parameters in this type of aquifer. The method can also identify the impact of measurement errors on the parameter estimates, and provide ranges of the actual parameter values. The parameter values calculated using this method were compared to those determined using other theories. It is found that this method is very useful for calculating the hydraulic properties from pumping test data and for analyzing the parameter reliability.     

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.     

MODELING CONCENTRATION VARIATIONS IN HIGH-CAPACITY WELLS: IMPLICATIONS FOR GROUNDWATER SAMPLING1

ABSTRACT: High-capacity wells are used as a convenient and economical means of sampling groundwater quality. Although the inherent limitations of using these wells are generally recognized, little has been done to investigate how these wells actually sample groundwater. A semi-analytical particle tracking model is used to illustrate the influence of variable vertical contaminant distributions and aquifer heterogeneity on the composition of water samples from these wells during short pumping periods. The hypothetical pumping well used in the simulations is located in an unconfined, alluvial aquifer with a shallow water table and concentration gradients of nitrate-nitrogen contamination. This is a typical setting for many irrigated areas in the United States. The main conclusions are: (1) high-capacity wells underestimate the average amount of contamination within an aquifer; (2) shapes of concentration-time curves for high-capacity wells appear to be governed by the distribution of the contaminant and travel times to the well; (3) variables such as well construction, pumping rate, and hydrogeologic properties contribute to the magnitude of the concentration-time curves at individual high-capacity wells; and (4) a sampling strategy using concentration-time curves based on the behavioral characteristics of the well rather than individual samples will provide a much better framework for interpreting spatial contaminant distributions.     

Is Standardized Precipitation Index (SPI) a Useful Indicator to Forecast Groundwater Droughts? — Insights from a Karst Aquifer

The study sought to understand the relationships between meteorological and groundwater droughts on water levels and spring discharges in Edwards Aquifer, Texas. Standardized Precipitation Index (SPI)‐styled Standardized Groundwater Index (SGI) was used to quantify groundwater droughts. SGI time series signal was delayed and damped, while SPI was volatile. SGI values correlated well with SPI values that were observed five to eight months ago. Dynamic regression models with lagged SPI terms and autoregressive integrated moving average errors indicated a statistically significant yet weak relationship between Lag‐1 SPI and SGI. The utility of SPI for groundwater drought forecasting was minimal in this aquifer. Nonseasonal and seasonal autoregressive terms played an important role in forecasting SGI and highlighted the need for long‐term, high‐resolution monitoring to properly characterize groundwater droughts. Spring flows exhibited stronger and quicker responses to meteorological droughts than changes in storage. In aquifers with spring discharges, groundwater monitoring programs must make efforts to inventory and monitor them. Groundwater drought contingency measures can be initiated using SPI but this indicator is perhaps inappropriate to remove groundwater drought restrictions.     

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.     

ANALYSIS OF PUMPING TESTS: SIGNIFICANCE OF WELL DIAMETER,PARTIAL PENETRATION,AND NOISE1

ABSTRACT: The nonlinear least squares (NLS) method was applied to pumping and recovery aquifer test data in confined and unconfined aquifers with finite diameter and partially penetrating pumping wells, and with partially penetrating piezometers or observation wells. It was demonstrated that noiseless and moderately noisy drawdown data from observation points located less than two saturated thicknesses of the aquifer from the pumping well produced an exact or acceptable set of parameters when the diameter of the pumping well was included in the analysis. The accuracy of the estimated parameters, particularly that of specific storage, decreased with increases in the noise level in the observed drawdown data. With consideration of the well radii, the noiseless drawdown data from the pumping well in an unconfined aquifer produced good estimates of horizontal and vertical hydraulic conductivities and specific yield, but the estimated specific storage was unacceptable. When noisy data from the pumping well were used, an acceptable set of parameters was not obtained. Further experiments with noisy drawdown data in an unconfined aquifer revealed that when the well diameter was included in the analysis, hydraulic conductivity, specific yield and vertical hydraulic conductivity may be estimated rather effectively from piezometers located over a range of distances from the pumping well. Estimation of specific storage became less reliable for piezometers located at distances greater than the initial saturated thickness of the aquifer. Application of the NLS to field pumping and recovery data from a confined aquifer showed that the estimated parameters from the two tests were in good agreement only when the well diameter was included in the analysis. Without consideration of well radii, the estimated values of hydraulic conductivity from the pumping and recovery tests were off by a factor of four.     

Effective Transboundary Aquifer Areas: An Approach for Transboundary Groundwater Management

The natural complexity, heterogeneity, and extent of transboundary aquifers around the world, have led to controversy over which method or criteria should be used to identify and delineate their boundaries. Currently, there is no standard methodology that aquifer‐sharing countries can use to delineate the area of a transboundary aquifer. In the case of Mexico and Texas, Mexico uses administrative boundaries, whereas Texas uses geological boundaries. This paper proposes a method for delineation and prioritization of aquifers (or aquifer areas) called effective transboundary aquifer areas (ETAAs), which uses a combination of physical criteria (geological boundaries, topography, and hydrography) and the location and density of active water wells in the borderland between Mexico and Texas. This method identifies the area of priority (productivity area) in the aquifer using pumping patterns or hot spots regardless of the aquifer’s surficial geological limits, therefore offering a more effective, local and practical management option at the transboundary level. Different geological features or pumping patterns will have different sizes and locations of ETAAs within the same aquifer. In West Texas, which is dominated by bolsons, the method produces limited options for ETAAs, whereas in South Texas in the easternmost border the identified ETAAs are more significant.     

Spatial and temporal trends of nitrate,chloride, and bromide concentration in an alluvial aquifer,North‐Central Texas,USA

The Seymour aquifer consists of unconfined outcrops of sand and gravel in a semiarid, agricultural region of north‐central Texas in the United States of America. Most water samples collected from the aquifer in 2015 had nitrate concentrations above the drinking water standard of 44.3 milligrams per liter (mg/L). Generally, areas with high nitrate concentration in 2010 remained high in 2015, although the median dropped by 3.9 mg/L. The largest decreases in nitrate concentration—up to 97 mg/L (60%)—were observed in wells with depths less than the median of 13.1 meters (m). However, other wells, including depths above and below the median, showed increases in nitrate concentration of up to 40 mg/L (42%). In 2015, chloride concentrations in six wells exceeded the secondary contaminant level of 250 mg/L, and one well had a chloride concentration of 1,810 mg/L. Past and ongoing agricultural practices, including cultivation of native grassland, application of fertilizer, and irrigation with nitrate‐contaminated groundwater, help sustain overall high nitrate concentrations within the aquifer. Local conditions governing nitrogen inputs and dilution result in significant improvement or worsening of the nitrate problem over relatively short timeframes. The pumping of groundwater from the aquifer may facilitate mixing with groundwater of increased salinity that has been affected by the dissolution of evaporites in underlying Permian bedrock.     

Evidence Supporting Cap and Trade as a Groundwater Policy Option for Reducing Irrigation Consumptive Use1

Thompson, Christopher L., Raymond J. Supalla, Derrel L. Martin, and Brian P. McMullen, 2009. Evidence Supporting Cap and Trade as a Groundwater Policy Option for Reducing Irrigation Consumptive Use. Journal of the American Water Resources Association (JAWRA) 45(6):1508‐1518. Abstract: In the American West water is becoming an increasingly scarce resource. Obligations to bordering states, endangered species protection, and long‐term resource sustainability objectives have created a need for most western states to reduce the consumption of irrigation water. In Nebraska specifically, the Nebraska Department of Natural Resources (NDNR) and local Natural Resource Districts (NRDs) are meeting a large part of this need by using a regulatory approach, commonly called groundwater allocation. The cost of allocation, which occurs in the form of reduced economic returns to irrigation, could be greatly reduced by using an integrated cap and trade approach. Much like environmental cap and trade programs which are used to reduce the cost of limiting environmental pollution, the trading of capped groundwater allocations can reduce the cost of limiting water use. In an analysis of a typical case in the Nebraska Republican Basin, we found that the impact of a water market to trade groundwater allocations depended on the size of the allocation and on the characteristics of the land and irrigation systems involved in the trade. Potential economic benefits from trade ranged from US$0 to US$120 per 1,000 cubic meters traded, from US$25 to US$250 per 1,000 cubic meters of reduction in consumptive use, and from US$16 to US$50 per hectare of irrigated land in the region. The highest benefits occurred at relatively high allocations, which capped withdrawals at 65‐75% of the expected unrestricted pumping level. These gains from trade would be split between buyers and sellers based on the negotiated selling price.     

Groundwaters at Risk: Wetland Loss Changes Sources,Lengthens Pathways,and Decelerates Rejuvenation of Groundwater Resources

Wetland loss alters the hydrology of wetlandscapes in poorly understood ways. To quantify the effects of wetland loss on subsurface hydrology, a physically based hydrologic model that simulates the timing and pathways of subsurface hydrologic connections was coupled with wetland inventories over a 50‐year period during which substantial wetland loss occurred. The model revealed, based on vertical variations in saturated hydraulic conductivities, wetland loss of different degrees led to a contraction of catchment contributing areas to local surface waters but an expansion of contributing areas to the regional surface water body. This shift in groundwater contributing areas reflected (1) a decrease in baseflow contribution to the local surface water bodies, and (2) an increase in the transit time and length of subsurface hydrologic connections with an associated increase in the magnitude and age of baseflow discharging to the regional surface water body. The model also showed regions with thick permeable aquifers were particularly sensitive to the loss of wetlands. Our ability to predict these changes in hydrology of the watershed provides important support for designing science‐based policies to promote sustainable water resource management.     

Analyzing the Effects of Groundwater Pumping on an Urban Stream‐Aquifer System

When the cone of influence of a pumping well reaches a nearby river, the resulting hydraulic gradient can induce enhanced seepage of streamflow into the aquifer. The rate of seepage is often modeled using analytical solutions that are simple to apply but may not reproduce field data due to mathematical assumptions not being met in the field. Furthermore, the appropriateness of such models has not been investigated in detail due to difficulty in measuring streamflow loss in the field. In this study, a field experiment was conducted on a reach of the South Platte River near Denver, Colorado to estimate pumping‐induced streamflow loss. A network of stream gauges, monitoring wells, and in situ measurements was used to observe streamflow rates, groundwater levels, and temperature to assess if pumping wells have a significant impact on streamflow, and to compare observed streamflow depletion against analytical solutions. Data collected suggest that pumping wells have a noticeable impact on streamflow. The analytical solutions proved accurate if streamflow was low and constant but performed poorly if streamflow was high and variable. Therefore, for this reach, the use of analytical solutions to predict streamflow may only be appropriate under low‐flow, constant‐flow conditions. Methods and results can be used to guide other streamflow depletion studies and to inform cases of pumping‐induced streamflow depletion, particularly in regard to water rights.     

Coupling of the Water Cycle with Patterns of Urban Growth in the Baltimore Metropolitan Region,United States

Regional municipal water plans typically do not recognize complex coupling patterns or that increased withdrawals in one location can result in changes in water availability in others. We investigated the interaction between urban growth and water availability in the Baltimore metropolitan region where urban growth has occurred beyond the reaches of municipal water systems into areas that rely on wells in low‐productivity Piedmont aquifers. We used the urban growth model SLEUTH and the hydrologic model ParFlow.CLM to evaluate this interaction with urban growth scenarios in 2007 and 2030. We found decreasing groundwater availability outside of the municipal water service area. Within the municipal service area we found zones of increasing storage resulting from increased urban growth, where reduced vegetation cover dominated the effect of urbanization on the hydrologic cycle. We also found areas of decreasing storage, where expanding impervious surfaces played a larger role. Although the magnitude of urban growth and change in water availability for the simulation period were generally small, there was considerable spatial heterogeneity of changes in subsurface storage. This suggests that there are locally concentrated areas of groundwater sensitivity to urban growth where water shortages could occur or where drying up of headwater streams would be more likely. The simulation approach presented here could be used to identify early warning indicators of future risk.     

Seasonal and Flood‐Induced Variations in Groundwater–Surface Water Exchange in a Northern Coldwater Fishery

Groundwater upwelling is important to coldwater fisheries survival. This study used stable isotopes to identify upwelling zones within a watershed, then combined isotope analyses with reach‐scale monitoring to measure surface water–groundwater exchange over time. Research focused on Amity Creek, Minnesota, a basin that exemplifies conditions limiting coldwater species survival along Lake Superior's North Shore where shallow bedrock limits groundwater capacity, lowering baseflows and increasing temperatures. Groundwater‐fed reaches were identified through synoptic isotope sampling, with results highlighting the importance of isolated shallow surficial aquifers (glacially derived sands and gravels) for providing cold baseflow waters. In an alluvial reach, monitoring well results show groundwater was stored in two reservoirs: one that reacts quickly to changes in stream levels, and one that remained isotopically isolated under most flow conditions, but which helps sustain summer baseflows for weeks to months. A 500‐year flood demonstrated the capacity of high‐flow events to alter surface water–groundwater connectivity. The previously isolated reservoir was exchanged or mixed during the flood pulse, while incision lowered the water table for years. The results here provide insight for streams that lack substantial groundwater inputs yet maintain coldwater species at risk in a warming climate and an approach for managers seeking to protect cold baseflow sources.