The Role of Climate and Human Influences in the Dry‐Up of the Jinci Springs,China1

Abstract: One of the largest karst springs in North China, the Jinci Springs, dried up and has remained dry since 1994. We develop a correlation analysis with time‐lag and a regression analysis with time‐lag to study the relation between spring flow and precipitation. This allows us to obtain a better understanding of karst hydrological processes by differentiating the contribution of variation in precipitation from anthropogenic impacts on the dry‐up of Jinci Springs. We divided the karstic hydrological processes into two phases: pre‐1961 and post‐1961. In the first phase (i.e., 1954‐1960) the groundwater recharge was affected by precipitation alone, and in the second phase (i.e., 1961‐1994) the groundwater recharge was influenced by both precipitation and human activities. Using precipitation and groundwater recharge data in the first phase, we set up a groundwater recharge model with time‐lags. By running the time‐lags model, we acquired the groundwater recharge likely to occur under the sole effect of precipitation in the second phase. Using a water‐balance calculation, we conclude that the groundwater recharge exhibited statistical stationarity, and the Jinci Springs dry‐up was the result of anthropogenic activities. At least three specific types of anthropogenic activities contributed to the drying‐up of Jinci Springs: (1) groundwater pumping accounts for 51%, (2) the dewatering from coal mining accounts for 33%, (3) and dam‐building 14%. The drying‐up of Jinci Springs meant that the groundwater drained from the aquifer’s fractures, and subsequently changed the structure of the karst aquifer. Although groundwater exploitation has been reduced, the flow at Jinci Springs has not reoccurred.     

Defining the Diurnal Pattern of Snowmelt Using a Beta Distribution Function

Snow is an important component of the hydrologic cycle for many regions worldwide. In addition to vital water resources, snowmelt can be important for forest ecosystem dynamics and flood risk. However, standard design events in the United States lack a design snowmelt event, including only precipitation events, though snowmelt has been shown to be larger than rainfall. In this article, we present a method using hourly snow water equivalent data to develop and test a function for representing the diurnal pattern of snowmelt. A two‐parameter beta distribution function is modified for the purposes of this study and found to fit the pattern of snowmelt well with a root mean squared error of 0.008. Soil moisture sensors were additionally utilized to assess the timing of the snowmelt water outflow from the base of the snowpack that supports the shape of the function, but suggests that the timing of losses recorded on snow pillows lag as much as 3 h. Further testing of the function showed the shape of the function to be accurate. The methods developed and tested in this paper can be applied for design purposes comparing snowmelt and rainfall events or to improve hydrological models investigating processes such as streamflow or groundwater recharge.     

Estimating Water Budgets and Vertical Leakages for Karst Lakes in North‐Central Florida (United States) Via Hydrological Modeling1

Lin, Zhulu, 2011. Estimating Water Budgets and Vertical Leakages for Karst Lakes in North‐Central Florida (United States) Via Hydrological Modeling. Journal of the American Water Resources Association (JAWRA) 1‐16. DOI: 10.1111/j.1752‐1688.2010.00513.x Abstract: Newnans, Lochloosa, and Orange Lakes are closely hydrologically connected karst lakes located in north‐central Florida, United States. The complex karst hydrology in this region poses a great challenge to the hydrological modeling that is essential to the development of Total Maximum Daily Loads for these lakes. We used a Hydrological Simulation Program – Fortran model coupled with the parallel Parameter ESTimation model calibration and uncertainty analysis software to estimate effectively the hydrological interactions between the lakes and the underlying upper Floridan aquifer and the water budgets for these three lakes. The net results of the lake‐groundwater interactions in Newnans and Orange Lakes are that both lakes recharge the underlying upper Floridan aquifer, with the recharge rate of the latter one magnitude greater than that of the former. However, for Lochloosa Lake, the net lake‐groundwater interaction is that the lake gains water from groundwater in a significant amount, approximately 40% of its total terrestrial water input. The annual average vertical leakages estimated for Newnans, Lochloosa, and Orange Lakes are 6.0 × 10⁶, ?8.9 × 10⁶, and 44.4 × 10⁶ m³, respectively. The average vertical hydraulic conductance (K_v/b) of the units between a lake bottom and the underlying upper Floridan aquifer in this region are also estimated to be from 1.26 × 10^?4 to 1.01 × 10^?3 day^?1.     

Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas1

Mittelstet, Aaron R., Michael D. Smolen, Garey A. Fox, and Damian C. Adams, 2011. Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas. Journal of the American Water Resources Association (JAWRA) 1‐8. DOI: 10.1111/j.1752‐1688.2011.00524.x Abstract: We compared two approaches to administration of groundwater law on a hydrologic model of the North Canadian River, an alluvial aquifer in northwestern Oklahoma. Oklahoma limits pumping rates to retain 50% aquifer saturated thickness after 20 years of groundwater use. The Texas Panhandle Groundwater Conservation District’s (GCD) rules limit pumping to a rate that consumes no more than 50% of saturated thickness in 50 years, with reevaluation and readjustment of permits every 5 years. Using a hydrologic model (MODFLOW), we simulated river‐groundwater interaction and aquifer dynamics under increasing levels of “development” (i.e., increasing groundwater withdrawals). Oklahoma’s approach initially would limit groundwater extraction more than the GCD approach, but the GCD approach would be more protective in the long run. Under Oklahoma rules more than half of aquifer storage would be depleted when development reaches 65%. Reevaluation of permits under the Texas Panhandle GCD approach would severely limit pumping as the 50% level is approached. Both Oklahoma and Texas Panhandle GCD approaches would deplete alluvial base flow at approximately 10% development. Results suggest periodic review of permits could protect aquifer storage and river base flow. Modeling total aquifer storage is more sensitive to recharge rate and aquifer hydraulic conductivity than to specific yield, while river leakage is most sensitive to aquifer hydraulic conductivity followed by specific yield.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

A Simple Process‐Based Snowmelt Routine to Model Spatially Distributed Snow Depth and Snowmelt in the SWAT Model1

Abstract: We present a method to integrate a process‐based (PB) snowmelt model that requires only daily temperature and elevation information into the Soil and Water Assessment Tool (SWAT) model. The model predicts the spatiotemporal snowpack distribution without adding additional complexity, and in fact reduces the number of calibrated parameters. To demonstrate the utility of the PB model, we calibrate the PB and temperature‐index (TI) SWAT models to optimize agreement with stream discharge on a 46‐km² watershed in northwestern Idaho, United States, for 10 individual years and use the calibrated parameters for the year with the best agreement to run the model for 15 remaining years. Stream discharge predictions by the PB and TI model were similar, although the PB model simulated snowmelt more accurately than the TI model for the remaining 15‐year period. Spatial snow distributions predicted by the PB model better matched observations from LandSat imagery and a SNOTEL station. Results for this watershed show that including PB snowmelt in watershed models is feasible, and calibration of TI‐based watershed models against discharge can incorrectly predict snow cover.     

The Effects of Irrigation and Climate on the High Plains Aquifer: A County‐Level Econometric Analysis

The High Plains Aquifer (HPA) underlies parts of eight states and 208 counties in the central area of the United States (U.S.). This region produces more than 9% of U.S. crops sales and relies on the aquifer for irrigation. However, these withdrawals have diminished the stock of water in the aquifer. In this paper, we investigate the aggregate county‐level effect on the HPA of groundwater withdrawal for irrigation, of climate variables, and of energy price changes. We merge economic theory and hydrological characteristics to jointly estimate equations describing irrigation behavior and a generalized water balance equation for the HPA. Our simple water balance model predicts, at average values for irrigation and precipitation, an HPA‐wide average decrease in the groundwater table of 0.47 feet per year, compared to 0.48 feet per year observed on average across the HPA during this 1985–2005 period. The observed distribution and predicted change across counties is in the (?3.22, 1.59) and (?2.24, 0.60) feet per year range, respectively. The estimated impact of irrigation is to decrease the water table by an average of 1.24 feet per year, whereas rainfall recharges the level by an average of 0.76 feet per year. Relative to the past several decades, if groundwater use is unconstrained, groundwater depletion would increase 50% in a scenario where precipitation falls by 25% and the number of degree days above 36°C doubles. Editor’s note : This paper is part of the featured series on Optimizing Ogallala Aquifer Water Use to Sustain Food Systems. See the February 2019 issue for the introduction and background to the series.     

Spatiotemporal Variability of Snow Depletion Curves Derived from SNODAS for the Conterminous United States, 2004‐2013

Assessment of water resources at a national scale is critical for understanding their vulnerability to future change in policy and climate. Representation of the spatiotemporal variability in snowmelt processes in continental‐scale hydrologic models is critical for assessment of water resource response to continued climate change. Continental‐extent hydrologic models such as the U.S. Geological Survey National Hydrologic Model (NHM) represent snowmelt processes through the application of snow depletion curves (SDCs). SDCs relate normalized snow water equivalent (SWE) to normalized snow covered area (SCA) over a snowmelt season for a given modeling unit. SDCs were derived using output from the operational Snow Data Assimilation System (SNODAS) snow model as daily 1‐km gridded SWE over the conterminous United States. Daily SNODAS output were aggregated to a predefined watershed‐scale geospatial fabric and used to also calculate SCA from October 1, 2004 to September 30, 2013. The spatiotemporal variability in SNODAS output at the watershed scale was evaluated through the spatial distribution of the median and standard deviation for the time period. Representative SDCs for each watershed‐scale modeling unit over the conterminous United States (n = 54,104) were selected using a consistent methodology and used to create categories of snowmelt based on SDC shape. The relation of SDC categories to the topographic and climatic variables allow for national‐scale categorization of snowmelt processes.     

Decentralized Groundwater Recharge Systems Using Roofwater and Stormwater Runoff1

Stephens, Daniel B., Mark Miller, Stephanie J. Moore, Todd Umstot, and Deborah J. Salvato, 2011. Decentralized Groundwater Recharge Systems Using Roofwater and Stormwater Runoff. Journal of the American Water Resources Association (JAWRA) 48(1): 134‐144. DOI: 10.1111/j.1752‐1688.2011.00600.x Abstract: Stormwater capture for groundwater recharge in urban areas is usually conducted at the regional level by water agencies. Field and modeling studies in New Mexico indicate that stormwater diverted to retention basins may recharge about 50% of precipitation that falls on the developed area, even in dry climates. Comparable volumes of recharge may be expected at homes, subdivisions, or commercial properties with low‐impact development (LID) technologies for stormwater control that promote recharge over evapotranspiration. Groundwater quality has not been significantly impacted at sites that have been recharging stormwater to aquifers for decades. Distributed recharge systems may be a good alternative to centralized regional facilities where there is limited land for constructing spreading basins or little funding for new infrastructure. LID technologies borrowed from stormwater managers are important tools for groundwater managers to consider to enhance recharge.     

iTree‐Hydro: Snow Hydrology Update For The Urban Forest Hydrology Model1

Yang, Yang, Theodore A. Endreny, and David J. Nowak, 2011. iTree‐Hydro: Snow Hydrology Update for the Urban Forest Hydrology Model. Journal of the American Water Resources Association (JAWRA) 47(6):1211–1218. DOI: 10.1111/j.1752‐1688.2011.00564.x Abstract: This article presents snow hydrology updates made to iTree‐Hydro, previously called the Urban Forest Effects—Hydrology model. iTree‐Hydro Version 1 was a warm climate model developed by the USDA Forest Service to provide a process‐based planning tool with robust water quantity and quality predictions given data limitations common to most urban areas. Cold climate hydrology routines presented in this update to iTree‐Hydro include: (1) snow interception to simulate the capture of snow by the vegetation canopy, (2) snow unloading to simulate the release of snow triggered by wind, (3) snowmelt to simulate the solid to liquid phase change using a heat budget, and (4) snow sublimation to simulate the solid to gas phase via evaporation. Cold climate hydrology routines were tested with research‐grade snow accumulation and weather data for the winter of 1996‐1997 at Umpqua National Forest, Oregon. The Nash‐Sutcliffe efficiency for open area snow accumulation was 0.77 and the Nash‐Sutcliffe efficiency for under canopy was 0.91. The USDA Forest Service offers iTree‐Hydro for urban forest hydrology simulation through http://www.iTreetools.org .     

Karst Hydrological Processes and Grey System Model1

Hao, Yonghong, Jiaojuan Zhao, Huamin Li, Bibo Cao, Zhongtang Li, and Tian‐Chyi J. Yeh, 2012. Karst Hydrological Processes and Grey System Model. Journal of the American Water Resources Association (JAWRA) 48(4): 656‐666. DOI: 10.1111/j.1752‐1688.2012.00640.x Abstract: The karst hydrological processes are the response of karst groundwater system to precipitation. This study provided a concept model of karst hydrological processes. The hydraulic response time of spring discharge to precipitation includes the time that precipitation penetrates through the vadose zone, and the subsequent groundwater pressure wave propagates to a spring outlet. Due to heterogeneities in karst aquifers, the hydraulic response time is different in different areas. By using grey system theory, we proposed a karst hydrological model that offers a calculation of hydraulic response time, and a response model of spring discharge to precipitation. Then, we applied the models to the Liulin Springs Basin, China. In the south part of the Liulin Springs Basin, where large fields of carbonate rocks outcrop with intensive karstification, the hydraulic response time is one year. In the north, where the Ordovician karst aquifer is covered by Quaternary loess sediments, the response time is seven years. The grey system GM(1,3) response model of spring discharge to precipitation was applied in consideration of the hydraulic response time. The model calibration showed that the average error was 6.55%, and validation showed that the average error was 12.19%.     

An Aquifer Storage and Recovery system with reclaimed wastewater to preserve native groundwater resources in El Paso, Texas

The traditional concept of Aquifer Storage and Recovery (ASR) has been emphasized and extensively applied for water resources conservation in arid and semi-arid regions using groundwater systems as introduced in Pyne's book titled Groundwater Recharge and Wells. This paper extends the ASR concept to an integrated level in which either treated or untreated surface water or reclaimed wastewater is stored in a suitable aquifer through a system of spreading basins, infiltration galleries and recharge wells; and part or all of the stored water is recovered through production wells, dual function recharge wells, or by streams receiving increased discharge from the surrounding recharged aquifer as needed. In this paper, the author uses the El Paso Water Utilities (EPWU) ASR system for injection of reclaimed wastewater into the Hueco Bolson aquifer as an example to address challenges and resolutions faced during the design and operation of an ASR system under a new ASR system definition. This new ASR system concept consists of four subsystems: source water, storage space-aquifer, recharge facilities and recovery facilities. Even though facing challenges, this system has successfully recharged approximately 74.7 million cubic meters (19.7 billion gallons) of reclaimed wastewater into the Hueco Bolson aquifer through 10 recharge wells in the last 18 years. This ASR system has served dual purposes: reuse of reclaimed wastewater to preserve native groundwater, and restoration of groundwater by artificial recharge of reclaimed wastewater into the Hueco Bolson aquifer.     

Streamflow Response to Climate as Influenced by Geology and Elevation1

Mayer, Timothy D. and Seth W. Naman, 2011. Streamflow Response to Climate as Influenced by Geology and Elevation. Journal of the American Water Resources Association (JAWRA) 47(4):724‐738. DOI: 10.1111/j.1752‐1688.2011.00537.x Abstract: This study examines the regional streamflow response in 25 predominately unregulated basins to warmer winter temperatures and snowpack reductions over the last half century in the Klamath Basin of California and Oregon. Geologic controls of streamflow in the region result in two general stream types: surface‐dominated and groundwater‐dominated basins. Surface‐dominated basins were further differentiated into rain basins and snowmelt basins on the basis of elevation and timing of winter runoff. Streamflow characteristics and response to climate vary with stream type, as discussed in the study. Warmer winter temperatures and snowpack reductions have caused significantly earlier runoff peaks in both snowmelt and groundwater basins in the region. In the groundwater basins, the streamflow response to changes in snowpack is smoothed and delayed and the effects are extended longer in the summer. Our results indicate that absolute decreases in July‐September base flows are significantly greater, by an order of magnitude, in groundwater basins compared to surface‐dominated basins. The declines are important because groundwater basins sustain Upper Klamath Lake inflows and mainstem river flows during the typically dry summers of the area. Upper Klamath Lake April‐September net inflows have decreased an estimated 16% or 84 thousand acre‐feet (103.6 Mm³) since 1961, with the summer months showing proportionately more decline. These changes will exacerbate water supply problems for agriculture and natural resources in the region.     

Long‐Term Trends in Streamflow and Precipitation in Northwest California and Southwest Oregon, 1953‐2012

Using nonparametric Mann‐Kendall tests, we assessed long‐term (1953‐2012) trends in streamflow and precipitation in Northern California and Southern Oregon at 26 sites regulated by dams and 41 “unregulated” sites. Few (9%) sites had significant decreasing trends in annual precipitation, but September precipitation declined at 70% of sites. Site characteristics such as runoff type (groundwater, snow, or rain) and dam regulation influenced streamflow trends. Decreasing streamflow trends outnumbered increasing trends for most months except at regulated sites for May‐September. Summer (July‐September) streamflow declined at many sites, including 73% of unregulated sites in September. Applying a LOESS regression model of antecedent precipitation vs. average monthly streamflow, we evaluated the underlying streamflow trend caused by factors other than precipitation. Decreasing trends in precipitation‐adjusted streamflow substantially outnumbered increasing trends for most months. As with streamflow, groundwater‐dominated sites had a greater percent of declining trends in precipitation‐adjusted streamflow than other runoff types. The most pristine surface‐runoff‐dominated watersheds within the study area showed no decreases in precipitation‐adjusted streamflow during the summer months. These results suggest that streamflow decreases at other sites were likely due to more increased human withdrawals and vegetation changes than to climate factors other than precipitation quantity.     

Observed Changes in Climate and Streamflow in the Upper Rio Grande Basin

Observed streamflow and climate data are used to test the hypothesis that climate change is already affecting Rio Grande streamflow volume derived from snowmelt runoff in ways consistent with model‐based projections of 21st‐Century streamflow. Annual and monthly changes in streamflow volume and surface climate variables on the Upper Rio Grande, near its headwaters in southern Colorado, are assessed for water years 1958–2015. Results indicate winter and spring season temperatures in the basin have increased significantly, April 1 snow water equivalent (SWE) has decreased by approximately 25%, and streamflow has declined slightly in the April–July snowmelt runoff season. Small increases in precipitation have reduced the impact of declining snowpack on trends in streamflow. Changes in the snowpack–runoff relationship are noticeable in hydrographs of mean monthly streamflow, but are most apparent in the changing ratios of precipitation (rain + snow, and SWE) to streamflow and in the declining fraction of runoff attributable to snowpack or winter precipitation. The observed changes provide observational confirmation for model projections of decreasing runoff attributable to snowpack, and demonstrate the decreasing utility of snowpack for predicting subsequent streamflow on a seasonal basis in the Upper Rio Grande Basin.     

Effect of Snow Cover Conditions on the Hydrologic Regime: Case Study in a Pluvial‐Nival Watershed,Japan1

Abstract: Hydrologic monitoring in a small forested and mountainous headwater basin in Niigata Prefecture has been undertaken since 2000. An important characteristic of the basin is that the hydrologic regime contains pluvial elements year‐round, including rain‐on‐snow, in addition to spring snowmelt. We evaluated the effect of different snow cover conditions on the hydrologic regime by analyzing observed data in conjunction with model simulations of the snowpack. A degree‐day snow model is presented and applied to the study basin to enable estimation of the basin average snow water equivalent using air temperature at three representative elevations. Analysis of hydrological time series data and master recession curves showed that flow during the snowmelt season was generated by a combination of ground water flow having a recession constant of 0.018/day and diurnal melt water flow having a recession constant of 0.015/hour. Daily flows during the winter/snowmelt season showed greater persistence than daily flows during the warm season. The seasonal water balance indicated that the ratio of runoff to precipitation during the cold season (December to May) was about 90% every year. Seasonal snowpack plays an important role in defining the hydrologic regime, with winter precipitation and snowmelt runoff contributing about 65% of the annual runoff. The timing of the snowmelt season, indicated by the date of occurrence of the first significant snowmelt event, was correlated with the occurrence of low flow events. Model simulations showed that basin average snow water equivalent reached a peak around mid‐February to mid‐March, although further validation of the model is required at high elevation sites.     

SIMULATING THREE-DIMENSIONAL GROUND WATER RESPONSE IN A SMALL MOUNTAINOUS WATERSHED1

ABSTRACT: Snowmelt from deep mountainous snowpacks is seldom rapid enough to exceed infiltration rates; thus, the source of streamflow in many mountainous watersheds is snowmelt recharge through shallow ground water systems. The hydrologic response and interaction between surface and sub-surface flow processes in these watersheds, which is controlled by basin structure, the spatial distribution of snowmelt, and the hydrogeology of the subsurface, are not well understood. The purpose of this study was to test a three-dimensional ground water model using simulated snowmelt input to simulate ground water response to spatially distributed snowmelt on the Upper Sheep Creek Watershed located within the Reynolds Creek Experimental Watershed in Southwestern Idaho. The model was used to characterize the mountainous aquifer and to delineate the subsurface flow mechanisms. Difficulty in finding a reasonable combination of grid spacing and time stepping within the model was encountered due to convergence problems with the Picard solution to the non-linear variably saturated ground water flow equations. Simulation results indicated that flow may be either unconfined or confined depending on inflow rate and hydrogeologic conditions in the watershed. The flow mechanism had a much faster response time when confined flow occurred. Response to snowmelt from a snow drift approximately 90 m away took only a few hours when flow was confined. Simulated results showed good agreement with piezometer measurements both in magnitude and timing; however, convergence problems with the Picard solution limited applicability of the model.     

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.     

CHARACTERIZING GROUND WATER FLOW IN THE MUNICIPAL WELL FIELDS OF CEDAR RAPIDS,IOWA, WITH SELECTED ENVIRONMENTAL TRACERS1

ABSTRACT: Cedar Rapids obtains its municipal water supply from a shallow alluvial aquifer along the Cedar River in east-central Iowa. Water samples were collected and analyzed for selected isotopes and chlorofluorocarbons to characterize the ground-water flow system near the municipal well fields. Analyses of deuterium and oxygen-18 indicate that water in the alluvial aquifer and in the underlying carbonate bedrock aquifer was recharged from precipitation during modern climatic conditions. Analyses of tritium indicate modern, post-1952, water in the alluvial aquifer and older, pre-1952, water in the bedrock aquifer. Mixing of the modern and older waters occurs in areas where (1) the confining layer between the two aquifers is discontinuous, (2) the bedrock aquifer is fractured, or (3) pumping of supply wells induces the flow of water between aquifers. Analyses of chlorofluorocarbons were used to determine the date of recharge of water samples. Water in the bedrock aquifer likely was recharged prior to the 1950s. Water in the alluvial aquifer likely was recharged from the 1960s to 1990s. Biodegradation or sorption probably affected some of the ground water analyzed for chlorofluorocarbons. These processes reduce the concentrations of CFCs, which results in older than actual calculated dates of recharge.