Search USGS⌕ Search

Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

At least 55 records · Page 3Linked to original sources

Conceptual and numerical groundwater flow model of the Cedar River alluvial aquifer system with simulation of drought stress on groundwater availability near Cedar Rapids, Iowa, for 2011 through 2013

Between July 2011 and February 2013, the City of Cedar Rapids observed water level declines in their horizontal collector wells approaching 11 meters. As a result, pumping from these production wells had to be halted, and questions were raised about the reliability of the alluvial aquifer under future drought conditions. The U.S. Geological Survey, in cooperation with the City of Cedar Rapids, completed a study to better understand the effects of drought stress on the Cedar River alluvial aquifer using a numerical groundwater flow model. Previously published groundwater flow models were combined with newly collected airborne, waterborne, down-hole, and land-based geophysical survey data and provided a detailed three-dimensional lithologic model of the Cedar River alluvial aquifer and surrounding area. An improved conceptual model for the groundwater flow system and a lithologic model were used to build and inform a numerical groundwater flow model capable of simulating water levels observed in the City of Cedar Rapids horizontal collector wells during the 2012 drought. Model performance was assessed primarily on the ability of the model to simulate water table elevation at six monitoring wells. Statistical tests were used to assess the numerical model during the calibration period, and results varied from satisfactory to unsatisfactory, likely because of stage changes in the Cedar River and production well withdrawal rates near monitoring wells. Simulated water levels during the 2012 drought indicated a depression near the horizontal collector wells, although simulated elevations at these locations and at monitoring wells were generally overestimated compared to measured values. The numerical groundwater flow model was modified to account for a decrease in seepage rate caused by low flow in the Cedar River and increased production. With seepage rate modification, model results improved; the simulated water table elevations were like those observed in horizontal collector and monitoring wells. Results demonstrated the ability of the model to simulate water levels observed in the horizontal collector wells during the 2012 drought when accounting for a decrease in infiltration from the Cedar River.

Iowa↗

Groundwater potentiometric-surface altitude in 2022 and groundwater-level changes between 1968, 1991, and 2022, in the alluvial aquifer in the Big Lost River Valley, south-central Idaho

The U.S. Geological Survey and the Idaho Department of Water Resources measured groundwater levels during spring 2022 and autumn 2022 to create detailed potentiometric-surface maps for the alluvial aquifer in the Big Lost River Valley in south-central Idaho. Wells were assigned to shallow, intermediate, and deep water-bearing units based on well depth, groundwater potentiometric-surface altitude, and hydrogeologic unit. Potentiometric-surface contours were created for each of the three water-bearing units for spring 2022 and autumn 2022. Groundwater flow generally follows topography down valley to the south. The groundwater-level data also were used to calculate changes in groundwater levels from spring to autumn 2022 and from historical measurement events in 1968 and 1991 to 2022. Groundwater levels declined at most wells from spring 1968 to spring 2022 and from spring 1991 to spring 2022. Although groundwater-level changes are sensitive to interannual wet and dry periods, long-term groundwater-level declines suggest that recharge and down-valley groundwater flows are insufficient to fully recover groundwater-level declines from pumping in some parts of the alluvial aquifer in the Big Lost River Valley.

Idaho↗

Assessment of groundwater availability in the Northern Atlantic Coastal Plain aquifer system From Long Island, New York, to North Carolina

Executive Summary The U.S. Geological Survey began a multiyear regional assessment of groundwater availability in the Northern Atlantic Coastal Plain (NACP) aquifer system in 2010 as part of its ongoing regional assessments of groundwater availability of the principal aquifers of the Nation. The goals of this national assessment are to document effects of human activities on water levels and groundwater storage, explore climate variability effects on the regional water budget, and provide consistent and integrated information that is useful to those who use and manage the groundwater resource. As part of this nationwide assessment, the USGS evaluated available groundwater resources within the NACP aquifer system from Long Island, New York, to northeastern North Carolina. The northern Atlantic Coastal Plain physiographic province depends heavily on groundwater to meet agricultural, industrial, and municipal needs. The groundwater assessment of the NACP aquifer system included an evaluation of how water use has changed over time; this evaluation primarily used groundwater budgets and development of a numerical modeling tool to assess system responses to stresses from future human uses and climate trends. This assessment focused on multiple spatial and temporal scales to examine changes in groundwater pumping, storage, and water levels. The regional scale provides a broad view of the sources and demands on the system with time. The sub-regional scale provides an evaluation of the differing response of the aquifer system across geographic areas allowing for closer examination of the interaction between different aquifers and confining units and the changes in these interactions under pumping and recharge conditions in 2013 and hydrologic stresses as much as 45 years in the future. By focusing on multiple scales, water-resource managers may utilize this study to understand system response to changes as they affect the system as a whole. The NACP aquifer system extends from Long Island to northeastern North Carolina, and includes aquifers primarily within New York, New Jersey, Delaware, Maryland, Virginia, and North Carolina. The seaward-dipping sedimentary wedge that underlies the northern Atlantic Coastal Plain physiographic province forms a complex groundwater system. Although the NACP aquifer system is recognized by the U.S. Geological Survey as one of the smallest of the 66 principal aquifer systems in the Nation, it ranks 13th overall in terms of total groundwater withdrawals and is 7th in population served. Despite abundant precipitation [about 45 inches per year (in/yr)], the supply of fresh surface water in this region is limited because many of the surface waters in this area are brackish estuaries, contributing to why many communities in the northern Atlantic Coastal Plain physiographic province rely heavily on groundwater to meet their water needs. Increases in population and changes in land use during the past 100 years have resulted in diverse increased demands for freshwater throughout the northern Atlantic Coastal Plain physiographic province with groundwater serving as a vital source of drinking water for the nearly 20 million people who live in the region. Total groundwater withdrawal in 2013 was estimated to be about 1,300 million gallons per day (Mgal/d) and accounts for about 40 percent of the drinking water supply with the densely populated areas tending to have the highest rates of withdrawals and, therefore, being most susceptible to effects from these withdrawals over time. Water levels in many of the confined aquifers are decreasing by as much as 2 feet per year (ft/yr) in response to extensive development and subsequent increased withdrawals throughout the region. Total water-level decreases (drawdowns) are more than 100 feet (ft) in some aquifers from their predevelopment (before 1900) levels. These drawdowns extend across state lines and under the Chesapeake and Delaware Bays, creating the potential for interstate aquifer management issues. Regional water-resources managers in the northern Atlantic Coastal Plain physiographic province face challenges beyond competing local domestic, industrial, agricultural, and environmental demands for water. Large changes in regional water use have made the State-level management of aquifer resources more difficult because of hydrologic effects that extend beyond State boundaries. The northern Atlantic Coastal Plain physiographic province is underlain by a wedge of unconsolidated to partially consolidated sediments that are typically thousands of feet thick along the coastline with a maximum thickness of about 10,000 ft near the edge of the continental shelf. The NACP aquifer system consists of nine confined aquifers and nine confining units capped by an unconfined surficial aquifer that is bounded laterally from the west by the contact between Coastal Plain sediments and the upland Piedmont bedrock. This aquifer system extends to the east to the limit of the Continental Shelf, however, the boundary between fresh and saline groundwater is considered to be much closer to the shoreline and varies vertically by aquifer. Precipitation over the region for average conditions from 2005 to 2009 is about 61,800 Mgal/d, but about 70 percent of it is lost to evapotranspiration resulting in an inflow of about 19,600 Mgal/d entering the groundwater system as aquifer recharge. Most of this recharge enters the aquifer system and flows through the shallow unconfined aquifer and either discharges to streams or directly to coastal waters without reaching the deep, confined aquifer system. In addition to recharge from precipitation, other sources of water include the return of wastewater from domestic septic systems of about 240 Mgal/d, about 60 Mgal/d of water released from storage in the confined system, and about 30 Mgal/d of lateral inflow at the boundary between freshwater and saltwater in response to pumping for conditions in 2013. The outflow needed to balance the inflows was subdivided between streamflow, discharge to tidal portions of streams, and coastal discharge. The hydrologic budget developed for current [2013] conditions determined that 93 percent of the total outflow was to surface waters with about 70 percent divided evenly between streamflow and shallow coastal discharge and 23 percent as discharge to tidal waters. The remaining 7 percent of the total outflow components include withdrawals from both the surficial and confined aquifers of the groundwater system. The groundwater availability assessment of the NACP aquifer system highlights the importance of analyses at both the regional and local scales to understand how changes in land use, water use, and climate have affected groundwater resources and how these resources may change in the future. The investigation included assessments of the regional changes in water levels and budgets across State lines, the importance of considering storage change in the confining units, the response of the aquifer system to a continuation of current [2013] hydrologic stresses into the future, and the potential effects of climate change and sea-level rise on the aquifer system. The Potomac aquifer group includes two of the most widely used aquifers in the NACP aquifer system, the Potomac-Patapsco and Potomac-Patuxent regional aquifers, providing about 24 percent of the total groundwater used in the region. Withdrawals from large pumping centers in this deep, confined aquifer group have resulted in substantial decreases in water-levels across state lines, particularly between southern Virginia and northeastern North Carolina as well as between southern New Jersey and northern Delaware where water levels in the Potomac-Patapsco aquifer have decreased by as much as 200 ft and 50 ft, respectively from predevelopment to current [2013] conditions. This response in water levels also is reflected in changes in water budgets where, for example, about 20 percent of the total response to pumping in Virginia is met by inducing flow from adjacent States. Understanding and quantifying these hydrologic effects that extend beyond State boundaries is critical for the State- and regional-level management of aquifer resources. The cumulative storage loss from the intervening confining units throughout the entire NACP aquifer system was about 35 percent of the total storage loss from predevelopment to current [2013] conditions. In geographic areas such as Delmarva Peninsula, Maryland, and New Jersey, the water released from storage in the confining units makes up the majority of the total storage release from the groundwater system and is becoming proportionally more important over time as the surficial aquifer approaches equilibrium with respect to pumping and recharge stresses as of 2013. Storage loss from the confining units is of particular concern because, unlike in the sands that comprise the confined aquifers, water removed from the clayey confining unit sediments cannot be replenished as these units gradually compress. This non-recoverable storage loss, if great enough, can result in land subsidence where these units are thick and the release from storage is relatively large and contributes to increased concerns for sea-level rise in areas such as the lower portion of the Chesapeake Bay. Groundwater usage increased dramatically in the NACP aquifer system during post-World War II era from the mid-1940s to early the 1980s, with withdrawals increasing from about 400 Mgal/d to more than 1,300 Mgal/d. Although groundwater withdrawals have been relatively constant since the early 1980s, about half of the total groundwater withdrawn from the NACP aquifer system since 1900 was withdrawn in the past 30 years. An analysis of the response of the groundwater system to a continuation of the current [2013] pumping for an additional 30 years into the future shows that the flow system continues to adjust in terms of changes in water budget components, water levels, and the boundary between freshwater and saltwater as it approaches equilibrium. The largest change in water budget components is the reduction in the amount of water released from storage. Across the entire NACP aquifer system, the reduction of storage release from 7 to 4 percent of the total water budget change is accounted for by reductions in groundwater discharge to streams and coastal waters. Locally, a similar response is calculated for each of the geographic areas except for Virginia where the amount of water released from storage accounts for about 25 percent of the total change in water budget. This finding suggests that the groundwater flow system in Virginia is not approaching equilibrium under the current [2013] stresses and, therefore, water levels will continue to decrease even if the pumping remains constant. An analysis of the change in water levels in the Potomac-Patapsco aquifer as pumping is continued 30 years into the future reveals that the largest decreases in water levels throughout the NACP aquifer system will occur in the southern Virginia and northeastern North Carolina parts of the study area. It is these areas that also see the greatest potential for increased lateral movement of saline groundwater in the deep, confined portion of the groundwater flow system in response to a continuation of the current [2013] pumping rates. The potential effects of long-term climate change and variability on the hydrologic system and availability of water resources in the NACP aquifer system continue to be of serious societal concern. These concerns include the effects of changes in aquifer recharge and in sea-level rise on the groundwater flow system. An assessment of the potential effects of a prolonged drought during current [2013] pumping conditions indicated that the reductions in recharge associated with droughts, including additional irrigation withdrawals required to meet increased crop water demand, have the greatest effects on water levels and streamflows in the surficial aquifer, and changes in water levels in the confined aquifers primarily resulted from the increased withdrawals associated with increased irrigation pumping; this response was most apparent in the Delmarva Peninsula. These results suggest that water levels may not be susceptible to the effects of droughts in the confined aquifers of the NACP aquifer system not used for irrigation, unlike in the unconfined surficial aquifer. A second analysis also was conducted to assess the effects of sea-level rise on the groundwater system throughout the northern Atlantic Coastal Plain physiographic province because recent analyses of the relative rates of sea-level rise along the Atlantic coast indicate that the Mid-Atlantic region represents a hot spot with anomalously higher rates of sea-level rise than observed elsewhere in the United States. Groundwater levels rose from 0 to 3 ft in response to a 3-ft simulated change in sea-level position, with the largest response occurring along the shoreline and away from non-tidal streams. About 37 percent (or 10,000 square miles) of the area of the northern Atlantic Coastal Plain physiographic province may experience about a 0.5-ft or more increase in water levels with the 3-ft increase in sea-level position, whereas about 18 percent (almost 5,000 square miles) of land of the northern Atlantic Coastal Plain physiographic province may experience a 2-ft or more increase in water levels with the 3-ft increase in sea-level position. These increases in the water table are of particular concern in low-lying areas where the unsaturated (vadose) zone is already thin, thus creating concerns for groundwater inundation of subsurface infrastructure, such as basements, septic systems, and subway systems. This increase in the water table also will likely alter the distribution of groundwater discharge to surface-water bodies thus increasing groundwater flow to streams that would have otherwise discharged directly to coastal waters. Throughout the NACP aquifer system, this redistribution of groundwater discharge results in an additional 2 percent of base flow in streams. Although the increases in groundwater discharge to streams (and corresponding decreases in discharge to coastal waters) calculated for the entire NACP aquifer system and its geographic areas represent only a small increase compared with current [2013] conditions, this redistribution of groundwater discharge from the coast to streams locally can alter the delivery of freshwater input to coastal receiving waters and have ecohydrological implications on the sensitive ecosystems which rely on a balance of groundwater discharge and surface-water flow.

Delaware, Maryland, New Jersey, New York, North Ca↗

Hydrogeology, numerical simulation of groundwater flow, and effects of future water use and drought for reach 1 of the Washita River alluvial aquifer, Roger Mills and Custer Counties, western Oklahoma, 1980–2015

The Washita River alluvial aquifer is a valley-fill and terrace alluvial aquifer along the valley of the Washita River in western Oklahoma that provides a productive source of groundwater for agricultural irrigation and water supply. The Oklahoma Water Resources Board (OWRB) has designated the westernmost section of the aquifer in Roger Mills and Custer Counties, Okla., as reach 1 of the Washita River alluvial aquifer; reach 1 is the focus of this report. The OWRB issued an order on November 13, 1990, that established the maximum annual yield (MAY; 120,320 acre-feet per year [acre-ft/yr]) and equal-proportionate-share (EPS) pumping rate (2.0 acre-feet per acre per year [(acre-ft/acre)/yr]) for reach 1 of the Washita River alluvial aquifer. The MAY and EPS were based on hydrologic investigations that evaluated the effects of potential groundwater withdrawals on groundwater availability in the Washita River alluvial aquifer. Every 20 years, the OWRB is statutorily required to update the hydrologic investigation on which the MAY and EPS were based. Because 30 years have elapsed since the last order was issued, the U.S. Geological Survey, in cooperation with the OWRB, conducted a new hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the Washita River alluvial aquifer. The Washita River is the primary source of inflow to Foss Reservoir, a Bureau of Reclamation reservoir constructed in 1961 for flood control, water supply, and recreation. Foss Reservoir provides water for Bessie, Clinton, New Cordell, and Hobart, Okla. Nearly 98 percent of the total groundwater use from the Washita River alluvial aquifer during 1967 to 2015 was for irrigation; other uses of groundwater in the study area include public supply, mining, and agriculture. A hydrogeologic framework was developed for the Washita River alluvial aquifer and included the physical characteristics of the aquifer, the geologic setting, the hydraulic properties of hydrogeologic units, the potentiometric surface (water table), and groundwater-flow directions at a scale that captures the regional controls on groundwater flow. The Washita River alluvial aquifer consists of alluvium and terrace deposits that were transported primarily by water and range from clay to gravel in size. The terrace includes windblown deposits of silt size and, in some cases, contains gravel laid down at several levels along former courses of present-day rivers. A conceptual flow model is a simplified description of the aquifer system that includes hydrologic boundaries, major inflow and outflow sources of the groundwater-flow system, and a conceptual water budget with the estimated mean flows between those hydrologic boundaries. During the study period 1980–2015, mean annual groundwater withdrawals, predominantly used for agricultural irrigation, totaled 5,502 acre-ft/yr, or 14 percent of aquifer outflows. When applied across the 132-square-mile aquifer area used for modeling purposes (84,366 acres), mean annual recharge of 3.15 inches per year corresponds to a mean annual recharge volume of 22,169 acre-ft/yr, or 56 percent of aquifer inflows. The annual saturated-zone evapotranspiration outflow was 11,828 acre-ft/yr for the Washita River alluvial aquifer, or about 30 percent of aquifer outflows. For the Washita River alluvial aquifer, lateral flow was 17,157 acre-ft/yr, or 44 percent of the aquifer inflows. The conceptual flow model and hydrogeologic framework were used to conceptualize, design, and build the numerical groundwater-flow model. A numerical groundwater-flow model of the Washita River alluvial aquifer was constructed by using MODFLOW-2005. The Washita River alluvial aquifer groundwater-model grid was spatially discretized into 350-foot (ft) cells and two layers. Layer 1 represented the undifferentiated alluvium and terrace deposits of Quaternary age, and layer 2 represented the bedrock of Permian age, which was given a uniform nominal thickness of 100 ft. The groundwater-simulation period was temporally discretized into 433 monthly transient stress periods, representing January 1980 to December 2015. An initial 365-day steady-state stress period was configured to represent mean annual inflows and outflows from the Washita River alluvial aquifer for the study period. The groundwater-flow model was calibrated manually and by automated adjustment of model inputs by using PEST++. Calibration targets for the Washita River alluvial aquifer model included groundwater-level observations and reservoir-stage observations, as well as base-flow and stream-seepage estimates. Three groundwater-availability scenarios were used in the calibrated groundwater model to (1) estimate the EPS pumping rate that retains the saturated thickness that meets the minimum 20-year life of the aquifer, (2) quantify the effects of projected pumping rates on groundwater storage over a 50-year period, and (3) evaluate how projected pumping rates extended 50 years into the future and sustained hypothetical drought conditions over a 10-year period affect base flow and groundwater in storage. The results of the groundwater-availability scenarios could be used by the OWRB to reevaluate the established MAY of groundwater from the Washita River alluvial aquifer. EPS scenarios for the Washita River alluvial aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 1.7, 1.6, and 1.6 (acre-ft/acre)/yr, respectively. Given the aquifer area used for modeling purposes (84,366 acres), these rates correspond to annual yields of 142,579, 134,986, and 134,986 acre-ft/yr, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 281,000 acre-feet (acre-ft), or about 306,000 acre-ft (52 percent) less than the starting storage. Considering the land-surface area of the Washita River alluvial aquifer and using a specific yield of 0.12, this decrease in storage was equivalent to a mean groundwater-level decline of about 30 ft. The Washita River downstream from Foss Reservoir and most of the streams in the study area were dry at the end of the 20-year EPS scenario. Foss Reservoir stage was below the dead-pool stage of 1,597 ft after about 7 years of pumping in the 20-year EPS scenario. Four projected 50-year groundwater-use scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage in the Washita River alluvial aquifer. These four scenarios used (1) no groundwater use, (2) groundwater use at the 2015 pumping rate, (3) mean groundwater use for the simulation period, and (4) increasing groundwater use. Groundwater storage after 50 years with no groundwater use was 545,249 acre-ft, or 693 acre-ft (0.1 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean groundwater-level increase of 0.1 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 543,831 acre-ft, or 723 acre-ft (0.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage after 50 years with the mean pumping rate for the study period was 543,202 acre-ft, or 1,349 acre-ft (0.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage at the end of the 50-year period with an increasing demand groundwater-pumping rate, which was 38 percent greater than the 2015 groundwater-pumping rate, was 542,584 acre-ft, or 1,967 acre-ft (0.4 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.2 ft. A hypothetical 10-year-drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage in the Washita River alluvial aquifer and Foss Reservoir stage and storage. To simulate the hypothetical drought, recharge in the calibrated model was reduced by 50 percent during the simulated drought period (1983–1992). Groundwater storage at the end of the drought period in December 1992 was 562,000 acre-ft, or 36,000 acre-ft (6 percent) less than the groundwater storage of the calibrated groundwater model (598,000 acre-ft). At the end of the hypothetical drought, the largest changes in saturated thickness (as great as 43.5 ft) were in the area upgradient from Foss Reservoir, particularly in the terrace at the model boundary. Substantial decreases in the Foss Reservoir stage began during the fall of 1985 in conjunction with base-flow decreases of up to 100 percent at U.S. Geological Survey streamgage 07324200 Washita River near Hammon, Okla. These lake-stage declines outpaced groundwater-level declines in the surrounding aquifer. The minimum Foss Reservoir storage simulated during the drought period was 77,954 acre-ft, which was a decrease of 46 percent from the nondrought storage.

Oklahoma↗

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967–2011, was divided into two sub-intervals because of varying water use, 1970–80 and 1981–2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981–2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967–2011; 20,724 acre-ft/yr during 1970–80, and 13,739 acre-ft/yr during 1981–2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967–2011; 19,963 acre-ft/yr during 1970–80; and 12,285 acre-ft/yr during 1981–2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967–2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980–2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980–2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7–279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4–279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981–2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994–2004) during the 1980–2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗

Quality of groundwater in the Denver Basin aquifer system, Colorado, 2003-5

Groundwater resources from alluvial and bedrock aquifers of the Denver Basin are critical for municipal, domestic, and agricultural uses in Colorado along the eastern front of the Rocky Mountains. Rapid and widespread urban development, primarily along the western boundary of the Denver Basin, has approximately doubled the population since about 1970, and much of the population depends on groundwater for water supply. As part of the National Water-Quality Assessment Program, the U.S. Geological Survey conducted groundwater-quality studies during 2003–5 in the Denver Basin aquifer system to characterize water quality of shallow groundwater at the water table and of the bedrock aquifers, which are important drinking-water resources. For the Denver Basin, water-quality constituents of concern for human health or because they might otherwise limit use of water include total dissolved solids, fluoride, sulfate, nitrate, iron, manganese, selenium, radon, uranium, arsenic, pesticides, and volatile organic compounds. For the water-table studies, two monitoring-well networks were installed and sampled beneath agricultural (31 wells) and urban (29 wells) land uses at or just below the water table in either alluvial material or near-surface bedrock. For the bedrock-aquifer studies, domestic- and municipal-supply wells completed in the bedrock aquifers were sampled. The bedrock aquifers, stratigraphically from youngest (shallowest) to oldest (deepest), are the Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers. The extensive dataset collected from wells completed in the bedrock aquifers (79 samples) provides the opportunity to evaluate factors and processes affecting water quality and to establish a baseline that can be used to characterize future changes in groundwater quality. Groundwater samples were analyzed for inorganic, organic, isotopic, and age-dating constituents and tracers. This report discusses spatial and statistical distributions of chemical constituents and evaluates natural and human-related processes that affect water quality. Findings are synthesized to assess the vulnerability of the Denver Basin aquifer system to groundwater contamination. The chemistry of groundwater samples collected from the water-table wells was generally different from that of samples collected from the bedrock-aquifer wells. Samples from the water-table wells tended to have higher concentrations of total dissolved solids and most major ions. Concentrations of several constituents with potential human-health concerns, including nitrate, selenium, uranium, and arsenic, decreased with depth and were highest in samples from the water-table wells. Exceedances of drinking-water standards and water-quality benchmarks were more frequently associated with shallow groundwater samples; concentrations of total dissolved solids and sulfate exceeded water-quality benchmarks for about half or more of samples from the water-table wells. The sediments and rocks of the Denver Basin are natural sources of the trace elements selenium, uranium, and arsenic, which affect their concentrations in groundwater. Detections of organic contaminants, which are typically indicative of human sources of contamination to groundwater, were more frequent in samples from the water-table wells. Pesticide compounds and volatile organic compounds were detected in 33 and 62 percent, respectively, of water-table well samples. Detected organic contaminant concentrations were much less than the associated drinking-water standards. Samples collected from the bedrock aquifers had lower concentrations of total dissolved solids than did samples collected from the water-table wells, although within the bedrock-aquifer samples, concentrations increased from the Dawson to Denver to Arapahoe to Laramie-Fox Hills aquifers. Concentrations of total dissolved solids and many constituents varied spatially and with depth in the bedrock aquifers, likely as a result of ion-exchange and oxidation-reduction reactions, which are important processes affecting water quality. Major-ion chemistry generally evolved from a calcium-bicarbonate to calcium-sulfate composition, with some sodium-bicarbonate and sodium-sulfate facies in the deeper bedrock aquifers, likely resulting from longer residence times and more extensive water-rock interaction. Oxidation-reduction conditions generally evolved from oxic at the water table to anoxic with increasing depth in the bedrock aquifers. Most samples from the bedrock aquifers were anoxic. Exceedances of drinking-water standards and water-quality benchmarks for the bedrock aquifers occurred in 1 percent or less of samples for nitrate, selenium, or arsenic; there were no exceedances for uranium. Exceedances for total dissolved solids, sulfate, manganese, and iron were generally between about 10 and 20 percent for the bedrock-aquifer samples. Radon concentrations, which were only measured in samples collected from two of the bedrock aquifers, exceeded the lower proposed drinking-water standard for more than 90 percent of samples but exceeded the higher alternative standard for less than 5 percent of samples. Pesticide compounds and volatile organic compounds were detected in 3 and 22 percent, respectively, of bedrock-aquifer samples, all at concentrations that were that were much less than drinking-water standards. Water-quality data were synthesized to evaluate factors that affect spatial and depth variability in water quality and to assess aquifer vulnerability to contaminants from geologic materials and those of human origin. The quality of shallow groundwater in the alluvial aquifer and shallow bedrock aquifer system has been adversely affected by development of agricultural and urban areas. Land use has altered the pattern and composition of recharge. Increased recharge from irrigation water has mobilized dissolved constituents and increased concentrations in the shallow groundwater. Concentrations of most constituents associated with poor or degraded water quality in shallow groundwater decreased with depth; many of these constituents are not geochemically conservative and are affected by geochemical reactions such as oxidation-reduction reactions. Groundwater age tracers provide additional insight into aquifer vulnerability and help determine if young groundwater of potentially poor quality has migrated to deeper parts of the bedrock aquifers used for drinking-water supply. Age-tracer results were used to group samples into categories of young, mixed, and old groundwater. Groundwater ages transitioned from mostly young in the water-table wells to mostly mixed in the shallowest bedrock aquifer, the Dawson aquifer, to mostly old in the deeper bedrock aquifers. Although the bedrock aquifers are mostly old groundwater of good water quality, several lines of evidence indicate that young, contaminant-bearing recharge has reached shallow to moderate depths in some areas of the bedrock aquifers. The Dawson aquifer is the most vulnerable of the bedrock aquifers to contamination, but results indicate that the older (deeper) bedrock aquifers are also vulnerable to groundwater contamination and that mixing with young recharge has occurred in some areas. Heavy pumping has caused water-level declines in the bedrock aquifers in some parts of the Denver Basin, which has the potential to enhance the transport of contaminants from overlying units. Results of this study are consistent with the existing conceptual understanding of aquifer processes and groundwater issues in the Denver Basin and add new insight into the vulnerability of the bedrock aquifers to groundwater contamination.

Colorado↗

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia↗

Hydrogeology and simulated groundwater flow and availability in the North Fork Red River aquifer, southwest Oklahoma, 1980–2013

On September 8, 1981, the Oklahoma Water Resources Board established regulatory limits on the maximum annual yield of groundwater (343,042 acre-feet per year) and equal-proportionate-share (EPS) pumping rate (1.0 acre-foot per acre per year) for the North Fork Red River aquifer. The maximum annual yield and EPS were based on a hydrologic investigation that used a numerical groundwater-flow model to evaluate the effects of potential groundwater withdrawals on groundwater availability in the North Fork Red River aquifer. The Oklahoma Water Resources Board is statutorily required (every 20 years) to update the hydrologic investigation on which the maximum annual yield and EPS were based. Because 20 years have elapsed since the final order was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the North Fork Red River aquifer in Oklahoma. This report describes a hydrologic investigation of the North Fork Red River aquifer that includes an updated summary of the aquifer hydrogeology. As part of this investigation, groundwater flow and availability were simulated by using a numerical groundwater-flow model. The North Fork Red River aquifer in Beckham, Greer, Jackson, Kiowa, and Roger Mills Counties in Oklahoma is composed of about 777 square miles (497,582 acres) of alluvium and terrace deposits along the North Fork Red River and tributaries, including Sweetwater Creek, Elk Creek, Otter Creek, and Elm Fork Red River. The North Fork Red River is the primary source of surface-water inflow to Lake Altus, which overlies the North Fork Red River aquifer. Lake Altus is a U.S. Bureau of Reclamation reservoir with the primary purpose of supplying irrigation water to the Lugert-Altus Irrigation District. A hydrogeologic framework was developed for the North Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of a numerical groundwater-flow model of the North Fork Red River aquifer described in this report. A conceptual model of aquifer inflows and outflows was developed for the North Fork Red River aquifer to constrain the construction and calibration of a numerical groundwater-flow model that reasonably represented the groundwater-flow system. The conceptual-model water budget estimated mean annual inflows to and outflows from the North Fork Red River aquifer for the period 1980–2013 and included a sub-accounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from Lake Altus. The numerical groundwater-flow model simulated the period 1980–2013 and was calibrated to water-table-altitude observations at selected wells, monthly base flow at selected streamgages, net streambed seepage as estimated for the conceptual model, and Lake Altus stage. Groundwater-availability scenarios were performed by using the calibrated numerical groundwater-flow model to (1) estimate the EPS pumping rate that guarantees a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical (10-year) drought on base flow and groundwater storage. The results of the groundwater-availability scenarios could be used by the Oklahoma Water Resources Board to reevaluate the maximum annual yield of groundwater from the North Fork Red River aquifer. EPS scenarios for the North Fork Red River aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 0.59, 0.52, and 0.52 acre-foot per acre per year, respectively. Given the 497,582-acre aquifer area, these rates correspond to annual yields of about 294,000, 259,000, and 259,000 acre-feet per year, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 951,000 acre-feet, or about 1,317,000 acre-feet (58 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean water-level decline of about 22 feet. Most areas of the active alluvium near the North Fork Red River, Elk Creek, and Elm Fork Red River remained partially saturated through the end of the EPS scenario because of streambed seepage. Lake Altus storage was reduced to zero after 6–7 years of EPS pumping in each scenario. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the North Fork Red River aquifer and base flows in the North Fork Red River upstream from Lake Altus. The effects of well withdrawals were evaluated by comparing changes in groundwater storage and base flow between four 50-year scenarios using (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2013), (3) 2013 pumping rates, and (4) increasing demand pumping rates. The increasing demand pumping rates assumed a 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwest Oklahoma. Groundwater storage after 50 years with no pumping was about 2,606,000 acre-feet, or 137,000 acre-feet (5.5 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 2.3 feet. Groundwater storage after 50 years with the mean pumping rate for the study period (1980–2013) was about 2,476,000 acre-feet, or about 7,000 acre-feet (0.3 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 0.1 foot. Groundwater storage at the end of the 50-year period with 2013 pumping rates was about 2,398,000 acre-feet, or about 70,000 acre-feet (2.8 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.2 feet. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was about 2,361,000 acre-feet, or about 107,000 acre-feet (4.3 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.8 feet. Mean annual base flow simulated at the Carter streamgage (07301500) on North Fork Red River increased by about 4,000 acre-feet (10 percent) after 50 years with no pumping and decreased by about 5,400 acre-feet (13 percent) after 50 years with increasing demand pumping rates. Mean annual base flow simulated at the North Fork Red River inflow to Lake Altus increased by about 7,400 acre-feet (15 percent) after 50 years with no pumping and decreased by about 5,800 acre-feet (12 percent) after 50 years with increasing demand pumping rates. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage and Lake Altus stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1984–1993). Groundwater storage at the end of the drought period was about 2,271,000 acre-feet, or about 426,000 acre-feet (15.8 percent) less than the groundwater storage of the calibrated numerical model. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 7.1 feet. At the end of the 10-year hypothetical drought period, base flows at the Sweetwater (07301420), Carter (07301500), Headrick (07305000), and Snyder (07307010) streamgages had decreased by about 37, 61, 44, and 45 percent, respectively. The minimum Lake Altus storage simulated during the drought period was 403 acre-feet, which was a decline of 92 percent from the nondrought storage. Reduced base flows in the North Fork Red River were the primary cause of Lake Altus storage declines.

Oklahoma↗

Groundwater conditions during 2009 and changes in groundwater levels from 1984 to 2009, Columbia Plateau Regional Aquifer System, Washington, Oregon, and Idaho

Groundwater elevations in three basalt units and one unconsolidated hydrogeologic unit in the Columbia Plateau Regional Aquifer System were measured and evaluated to provide a regional overview of groundwater conditions in spring 2009. Water levels for the Saddle Mountains unit, the Wanapum unit, the Grande Ronde unit, and for the overlying Overburden unit were measured in 1,752 wells during spring 2009 by the U.S. Geological Survey (USGS) and 10 other Federal, State, Tribal, and local agencies, including 66 wells located and measured by the USGS specifically for this study. These data were analyzed to determine the presence of spatial correlation of groundwater levels with distance and direction from each other. Groundwater flow in the Palouse Slope structural region showed evidence of being more continuous relative to groundwater flow in the Yakima Fold Belt, where the geologic complexity may contribute to compartmentalization of groundwater flow. This information was used to interpolate the generalized groundwater elevations for each of the basalt hydrogeologic units and to provide information on regional flow. Water-level change maps were constructed for the three basalt hydrogeologic units and the Overburden (unconsolidated) unit. Groundwater levels measured in spring 1984 and 2009 in 470 wells were compared. Small to moderate groundwater-level declines were measured in most wells, although declines greater than 100 ft and as great as 300 ft were measured in many wells. Essentially unchanged groundwater levels were measured in other wells. Of the wells measured in 1984 and 2009, water levels declined in 83 percent of the wells, and declines greater than 25 ft were measured in 29 percent of all wells. The groundwater-level changes were greatest in the deeper hydrogeologic units. Mean groundwater-level changes ranged from a 7 ft decline for the Overburden unit to a 51 ft decline for the Grande Ronde unit. The average annual rates of groundwater-level change for the 25-year period ranged from a 0.3 ft/yr decline for the Overburden unit to a 2.0 ft/yr decline for the Grande Ronde unit. Groundwater level declines were identified throughout the Columbia Plateau, but areas with large and widespread declines were located in the central northern part of the study area, in parts of the Yakima River basin in Washington, in the Pullman-Moscow area in Washington and Idaho, and in parts of the Umatilla River basin in Oregon. These declines are in areas known to rely heavily on groundwater for irrigation and other uses.

Idaho, Oregon, Washington↗

Groundwater quality and age of secondary bedrock aquifers in the glaciated portion of eastern Nebraska, 2016–18

The Eastern Nebraska Water Resources Assessment (ENWRA) project was initiated in 2006 to assist water managers by developing a hydrogeologic framework and water budget for the glaciated portion of eastern Nebraska. Within the ENWRA area, the primary groundwater sources for municipal, domestic, and irrigation water needs are provided by withdrawals from alluvial, buried paleovalley, and the High Plains aquifer (where present). Generally, other bedrock aquifers are considered a secondary water source. However, in some areas, such as parts of Sarpy and Nemaha Counties, these secondary bedrock aquifers are the only source of water within glaciated upland areas. To improve the understanding of the quality, geochemistry, and age of groundwater from bedrock aquifers, the U.S. Geological Survey (USGS), in cooperation with the ENWRA group, which includes the Lewis and Clark, Lower Elkhorn, Lower Platte North, Lower Platte South, Nemaha, and Papio-Missouri River Natural Resources Districts, designed a study to sample 31 wells completed in the secondary bedrock aquifers and analyze samples for major ions, physical properties, nutrients, stable isotopes, and selected age tracers. Of the 31 samples collected for this report, 22 samples were collected from the Dakota aquifer contained in the Dakota Sandstone, 3 from the Niobrara aquifer contained in the Niobrara Formation of Colorado Group, and 6 from Paleozoic aquifers contained in undifferentiated Paleozoic-age units. The results of this study indicate that major ion data collected from the Dakota aquifer can be used for assessing the quality, recharge source, and age of groundwater. Calcium bicarbonate dominant samples were characterized as modern or mixed, indicating that, in these areas, groundwater is unconfined and is recharged by precipitation and (or) surface water. If groundwater extraction rates exceed recharge rates, total dissolved solid concentrations may increase as a result of upwelling of groundwater from deeper units or formations, which can adversely affect groundwater quality. Sampling results presented in this report indicate water quality is good, but that groundwater in the Dakota aquifer with calcium bicarbonate water type may be vulnerable to surface contamination. In contrast, groundwater sampled from the Dakota aquifer, having a dominant water type other than calcium bicarbonate, generally has low dissolved oxygen and nitrate concentrations, and higher concentrations of total dissolved solids and trace elements, including iron and strontium. The geochemical characteristics of noncalcium bicarbonate samples from the Dakota aquifer indicated confining conditions and limited groundwater recharge from local precipitation. Apparent groundwater ages estimated from radiocarbon (carbon-14) sampling of noncalcium bicarbonate samples from the Dakota aquifer indicated that the time of groundwater recharge to the Dakota aquifer occurred during Pleistocene time. Depleted stable isotopes results indicate recharge during a colder climate. Groundwater under confined conditions is not easily recharged from precipitation or surface water. Future groundwater-level monitoring in locations where the Dakota aquifer appears to be confined could provide information to evaluate whether groundwater supplies remain sufficient to meet future municipal, domestic, and irrigation needs. For the Niobrara aquifer and Paleozoic aquifers, the dominant water type was not a diagnostic indicator of recharge source, age, and groundwater quality as with the Dakota aquifer. Most likely this is because the host formation was dominated by calcium-carbonate-rich rocks; however, few samples were collected from these aquifers to be able to confirm this interpretation. Samples collected from wells completed in the Niobrara aquifer and Paleozoic aquifers and characterized as calcium sulfate water type have statistically significantly higher concentrations of total dissolved solids compared to other samples from the Niobrara aquifer and Paleozoic aquifers characterized as calcium bicarbonate. Given that six of the nine of samples collected from the Niobrara and Paleozoic aquifers indicated modern recharge, these secondary bedrock aquifers are reliant on precipitation to sustain groundwater levels and may be vulnerable to a multiyear drought. Well yields of the Niobrara and Paleozoic aquifers are dependent on the presence of secondary porosity and these units offer little storage. Samples collected from wells completed in Paleozoic aquifers were the most isotopically enriched and similar to modern precipitation and had the highest concentrations of nitrate, indicating that groundwater is affected by agricultural activities. Future groundwater sampling would be beneficial to characterize groundwater-quality changes within the Niobrara and Paleozoic aquifers over time.

Nebraska↗

Hydrogeology and simulation of groundwater flow and analysis of projected water use for the Canadian River alluvial aquifer, western and central Oklahoma

This report describes a study of the hydrogeology and simulation of groundwater flow for the Canadian River alluvial aquifer in western and central Oklahoma conducted by the U.S. Geological Survey in cooperation with the Oklahoma Water Resources Board. The report (1) quantifies the groundwater resources of the Canadian River alluvial aquifer by developing a conceptual model, (2) summarizes the general water quality of the Canadian River alluvial aquifer groundwater by using data collected during August and September 2013, (3) evaluates the effects of estimated equal proportionate share (EPS) on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future by using numerical groundwater-flow models, and (4) evaluates the effects of present-day groundwater pumping over a 50-year period and sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage by using numerical flow models. The Canadian River alluvial aquifer is a Quaternary-age alluvial and terrace unit consisting of beds of clay, silt, sand, and fine gravel sediments unconformably overlying Tertiary-, Permian-, and Pennsylvanian-age sedimentary rocks. For groundwater-flow modeling purposes, the Canadian River was divided into Reach I, extending from the Texas border to the Canadian River at the Bridgeport, Okla., streamgage (07228500), and Reach II, extending downstream from the Canadian River at the Bridgeport, Okla., streamgage (07228500), to the confluence of the river with Eufaula Lake. The Canadian River alluvial aquifer spans multiple climate divisions, ranging from semiarid in the west to humid subtropical in the east. The average annual precipitation in the study area from 1896 to 2014 was 34.4 inches per year (in/yr). A hydrogeologic framework of the Canadian River alluvial aquifer was developed that includes the areal and vertical extent of the aquifer and the distribution, texture variability, and hydraulic properties of aquifer materials. The aquifer areal extent ranged from less than 0.2 to 8.5 miles wide. The maximum aquifer thickness was 120 feet (ft), and the average aquifer thickness was 50 ft. Average horizontal hydraulic conductivity for the Canadian River alluvial aquifer was calculated to be 39 feet per day, and the maximum horizontal hydraulic conductivity was calculated to be 100 feet per day. Recharge rates to the Canadian River alluvial aquifer were estimated by using a soil-water-balance code to estimate the spatial distribution of groundwater recharge and a water-table fluctuation method to estimate localized recharge rates. By using daily precipitation and temperature data from 39 climate stations, recharge was estimated to average 3.4 in/yr, which corresponds to 8.7 percent of precipitation as recharge for the Canadian River alluvial aquifer from 1981 to 2013. The water-table fluctuation method was used at one site where continuous water-level observation data were available to estimate the percentage of precipitation that becomes groundwater recharge. Estimated annual recharge at that site was 9.7 in/yr during 2014. Groundwater flow in the Canadian River alluvial aquifer was identified and quantified by a conceptual flow model for the period 1981–2013. Inflows to the Canadian River alluvial aquifer include recharge to the water table from precipitation, lateral flow from the surrounding bedrock, and flow from the Canadian River, whereas outflows include flow to the Canadian River (base-flow gain), evapotranspiration, and groundwater use. Total annual recharge inflows estimated by the soil-water-balance code were multiplied by the area of each reach and then averaged over the simulated period to produce an annual average of 28,919 acre-feet per year (acre-ft/yr) for Reach I and 82,006 acre-ft/yr for Reach II. Stream base flow to the Canadian River was estimated to be the largest outflow of groundwater from the aquifer, measured at four streamgages, along with evapotranspiration and groundwater use, which were relatively minor discharge components. Objectives for the numerical groundwater-flow models included simulating groundwater flow in the Canadian River alluvial aquifer from 1981 to 2013 to address groundwater use and drought scenarios, including calculation of the EPS pumping rates. The EPS for the alluvial and terrace aquifers is defined by the Oklahoma Water Resources Board as the amount of fresh water that each landowner is allowed per year per acre of owned land to maintain a saturated thickness of at least 5 ft in at least 50 percent of the overlying land of the groundwater basin for a minimum of 20 years. The groundwater-flow models were calibrated to water-table altitude observations, streamgage base flows, and base-flow gain to the Canadian River. The Reach I water-table altitude observation root-mean-square error was 6.1 ft, and 75 percent of residuals were within ±6.7 ft of observed measurements. The average simulated stream base-flow residual at the Bridgeport streamgage (07228500) was 8.8 cubic feet per second (ft 3 /s), and 75 percent of residuals were within ±30 ft 3 /s of observed measurements. Simulated base-flow gain in Reach I was 8.8 ft 3 /s lower than estimated base-flow gain. The Reach II water-table altitude observation root-mean-square error was 4 ft, and 75 percent of residuals were within ±4.3 ft of the observations. The average simulated stream base-flow residual in Reach II was between 35 and 132 ft 3 /s. The average simulated base-flow gain residual in Reach II was between 11.3 and 61.1 ft 3 /s. Several future predictive scenarios were run, including estimating the EPS pumping rate for 20-, 40-, and 50-year life of basin scenarios, determining the effects of current groundwater use over a 50-year period into the future, and evaluating the effects of a sustained drought on water availability for both reaches. The EPS pumping rate was determined to be 1.35 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 3.08 (acre-ft/acre)/yr in Reach II for a 20-year period. For the 40- and 50-year periods, the EPS rate was determined to be 1.34 (acre-ft/acre)/yr in Reach I and 3.08 (acre-ft/acre)/yr in Reach II. Storage changes decreased in tandem with simulated groundwater pumping and were minimal after the first 15 simulated years for Reach I and the first 8 simulated years for Reach II. Groundwater pumping at year 2013 rates for a period of 50 years resulted in a 0.2-percent decrease in groundwater-storage volumes in Reach I and a 0.6-percent decrease in the groundwater-storage volumes in Reach II. The small changes in storage are due to groundwater use by pumping, which composes a small percentage of the total groundwater-flow model budgets for Reaches I and II. A sustained drought scenario was used to evaluate the effects of a hypothetical 10-year drought on water availability. A 10-year period was chosen where the effects of drought conditions would be simulated by decreasing recharge by 75 percent. In Reach I, average simulated stream base flow at the Bridgeport streamgage (07228500) decreased by 58 percent during the hypothetical 10-year drought compared to average simulated stream base flow during the nondrought period. In Reach II, average simulated stream base flows at the Purcell streamgage (07229200) and Calvin streamgage (07231500) decreased by 64 percent and 54 percent, respectively. In Reach I, the groundwater-storage drought scenario resulted in a storage decline of 30 thousand acre-feet, or an average decline in the water table of 1.2 ft. In Reach II, the groundwater-storage drought scenario resulted in a storage decline of 71 thousand acre-feet, or an average decline in the water table of 2.0 ft.

Oklahoma↗

Aquifers of Arkansas: protection, management, and hydrologic and geochemical characteristics of groundwater resources in Arkansas

Sixteen aquifers in Arkansas that currently serve or have served as sources of water supply are described with respect to existing groundwater protection and management programs, geology, hydrologic characteristics, water use, water levels, deductive analysis, projections of hydrologic conditions, and water quality. State and Federal protection and management programs are described according to regulatory oversight, management strategies, and ambient groundwater-monitoring programs that currently (2013) are in place for assessing and protecting groundwater resources throughout the State. Physical attributes, groundwater geochemistry, and groundwater quality are described for each of the 16 aquifers of the State. Information in regard to the hydrology and geochemistry of each of the aquifers is summarized from about 550 historical and recent publications. Additionally, more than 8,000 sites with groundwater-quality data were obtained from the U.S. Geological Survey National Water Information System and the Arkansas Department of Environmental Quality databases and entered into a spatial database to investigate distribution and trends in chemical constituents for each of the aquifers. The 16 aquifers of the State were divided into two major physiographic regions of the State: the Coastal Plain Province (referred to as Coastal Plain) of eastern and southern Arkansas, which includes 11 of the 16 aquifers, and the Interior Highlands Division (referred to as Interior Highlands) of western Arkansas, which includes the remaining 5 aquifers. The 11 aquifers in the Coastal Plain consist of various geologic units that are Cenozoic in age and consist primarily of Cretaceous, Tertiary, and Quaternary sands, gravels, silts, and clays. Groundwater in the Coastal Plain represents one of the most valuable natural resources in the State, driving the economic engines of agriculture, while also supplying abundant water for commercial, industrial, and public-supply use. In terms of age from youngest to oldest, the aquifers of the Coastal Plain include Quaternary alluvial aquifers, including the Mississippi River Valley alluvial aquifer (the most important aquifer in Arkansas in terms of volume of use and economic benefits), the Jackson Group (a regional confining unit that served for decades as an important source of domestic supply), and the Cockfield, Sparta, Cane River, Carrizo, Wilcox, Nacatoch, Ozan, Tokio, and Trinity aquifers. The Mississippi River Valley alluvial aquifer accounts for approximately 94 percent of all groundwater used in the State, and the aquifer is used primarily for irrigation purposes. The Sparta aquifer is the second most important aquifer in terms of use, and the aquifer was used in the past dominantly as a source of public and industrial supply, although increasing irrigation use is occurring because of critically declining water levels in the Mississippi River Valley alluvial aquifer. Other aquifers of the Coastal Plain generally are used as important local sources of domestic, industrial, and public supply, in addition to other minor uses. Water quality generally is good for all aquifers of the Coastal Plain, except for elevated iron concentrations and localized areas of high salinity. The high salinity results from intrusion from underlying formations, evapotranspiration processes in areas of low recharge, and inadequate flushing in downgradient areas of residual salinity from deposition in marine environments. Trends in the spatial distribution of individual chemical constituents are related to position along the flow path for most aquifers of the Coastal Plain. These trends include elevated iron and nitrate concentrations with lower pH values and dissolved solids in groundwater from the outcrop areas, transitioning to lower iron and nitrate (related to changes in redox) and higher pH and dissolved solids (dominantly from the dissolution of carbonate minerals) in groundwater downgradient from outcrop areas. Groundwater generally trended from a calcium- to a sodium-bicarbonate water type with increasing cation exchange along the flow path. The Interior Highlands of western Arkansas has less reported groundwater use than other areas of the State, reflecting a combination of factors. These factors include prevalent and increasing use of surface water, less intensive agricultural uses, lower population and industry densities, lesser potential yield of the resource, and lack of detailed reporting. The overall low yields of aquifers of the Interior Highlands result in domestic supply as the dominant use, with minor industrial, public, and commercial-supply use. Where greater volumes are required for growth of population and industry, surface water is the greatest supplier of water needs in the Interior Highlands. The various aquifers of the Interior Highlands generally occur in shallow, fractured, well-indurated, structurally modified bedrock of this mountainous region of the State, as compared to the relatively flat-lying, unconsolidated sediments of the Coastal Plain. In terms of age from youngest to oldest, the aquifers of the Interior Highlands include: the Arkansas River Valley alluvial aquifer, the Ouachita Mountains aquifer, the Western Interior Plains confining system, the Springfield Plateau aquifer, and the Ozark aquifer. Spatial trends in groundwater geochemistry in the Interior Highlands differ greatly from trends noted for aquifers of the Coastal Plain. In the Coastal Plain, the prevalence of long regional flow paths results in regionally predictable and mappable geochemical changes along the flow paths. In the Interior Highlands, short, topographically controlled flow paths (from hilltops to valleys) within small watersheds represent the predominant groundwater-flow system. As such, dense data coverage from numerous wells would be required to effectively characterize these groundwater basins and define small-scale geochemical changes along any given flow path for aquifers of the Interior Highlands. Changes in geochemistry generally were related to rock type and residence time along individual flow paths. Dominant changes in geochemistry for the Ouachita Mountains aquifer and the Western Interior Plains confining system are attributed to rock/water interaction and changes in redox zonation along the flow path. In these areas, groundwater evolves along flow paths from a calcium- to a sodium-bicarbonate water type with increasing reducing conditions resulting in denitrification, elevated iron and manganese concentrations, and production of methane in the more geochemically evolved and strongest reducing conditions. In the Ozark and Springfield Plateau aquifers, rapid influx of surface-derived contaminants, especially nitrogen, coupled with few to no attenuation processes was attributed to the karst landscape developed on Mississippian- and Ordovician-age carbonate rocks of the Ozark Plateaus. Increasing nitrate concentrations are related to increasing agricultural land use, and areas of mature karst development result in higher nitrate concentrations than areas with less karst features.

Arkasas↗

Status and understanding of groundwater quality in the Monterey-Salinas Shallow Aquifer Study Unit, 2012–13: California GAMA Priority Basin Project

Groundwater quality in the approximately 7,820-square-kilometer (km 2 ) Monterey-Salinas Shallow Aquifer (MS-SA) study unit was investigated from October 2012 to May 2013 as part of the second phase of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is in the central coast region of California in the counties of Santa Cruz, Monterey, and San Luis Obispo. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in cooperation with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The MS-SA study was designed to provide a statistically robust assessment of untreated-groundwater quality in the shallow aquifer systems. The assessment was based on water-quality samples collected by the U.S. Geological Survey from 100 groundwater sites and 70 household tap sites, along with ancillary data such as land use and well-construction information. The shallow aquifer systems were defined by the depth interval of wells associated with domestic supply. The MS-SA study unit consisted of four study areas—Santa Cruz (210 km 2 ), Pajaro Valley (360 km 2 ), Salinas Valley (2,000 km 2 ), and Highlands (5,250 km 2 ). This study had two primary components: the status assessment and the understanding assessment . The first primary component of this study—the status assessment —assessed the quality of the groundwater resource indicated by data from samples analyzed for volatile organic compounds (VOCs), pesticides, and naturally present inorganic constituents, such as major ions and trace elements. The status assessment is intended to characterize the quality of groundwater resources in the shallow aquifer system of the MS-SA study unit, not the treated drinking water delivered to consumers by water purveyors. As opposed to the public wells, however, water from private wells, which often tap the shallow aquifer, is usually consumed without any treatment. The second component of this study—the understanding assessment —identified the natural and human factors that potentially affect groundwater quality by evaluating land-use characteristics, measures of location, geologic factors, groundwater age, and geochemical conditions of the shallow aquifer. An additional component of this study was a comparison of MS-SA water-quality results to those of the GAMA Monterey Bay and Salinas Valley Groundwater Basins study unit. This study unit covered much of the same areal extent as the MS-SA, but assessed the deeper, public drinking-water aquifer system. Relative concentrations (sample concentration divided by the benchmark concentration) were used to evaluate concentrations of constituents in groundwater samples relative to water-quality benchmarks for those constituents that have Federal or California benchmarks, such as maximum contaminant levels. For organic and special-interest constituents, relative concentrations were classified as high, greater than 1.0; moderate, greater than 0.1 and less than or equal to 1.0; or low, less than or equal to 0.1. For inorganic constituents, relative concentrations were classified as high, greater than 1.0; moderate, greater than 0.5 and less than or equal to 1.0; or low, less than or equal to 0.5. A relative concentration greater than 1.0 indicates that the concentration was greater than a benchmark. Aquifer-scale proportions were used to quantify regional-scale groundwater quality. The aquifer-scale proportions are the areal percentages of the shallow aquifer system where relative concentrations for a given constituent or class of constituents were high, moderate, or low. Inorganic constituents were measured at high and moderate relative concentrations more frequently than organic constituents. In the MS-SA study unit, inorganic constituents with benchmarks were detected at high relative concentrations in 51 percent of the study unit. The greatest proportions of high relative concentrations of trace elements and radioactive constituents were in the Highlands and Santa Cruz study areas, whereas high relative concentrations of nutrients were most often detected in the Salinas Valley and Pajaro Valley study areas and salinity indicators were most often detected in the Highlands and Salinas Valley study areas. The trace elements detected at high relative concentrations were arsenic, boron, iron, manganese, molybdenum, selenium, and strontium. The radioactive constituents detected at high relative concentrations were adjusted gross alpha radioactivity and uranium. The nutrient detected at high relative concentrations was nitrate plus nitrite. The salinity indicators detected at high relative concentrations were chloride, sulfate, and total dissolved solids. Organic constituents (VOCs and pesticides) were not detected at high relative concentrations in any of the study areas. The fumigant 1,2-dichloropropane was detected at moderate relative concentrations. The VOC chloroform and the pesticide simazine were the only organic constituents detected in more than 10 percent of samples. The constituents of special interest NDMA ( N -nitrosodimethylamine) and perchlorate were detected at high relative concentrations in the MS-SA study unit. Selected constituents were evaluated with explanatory factors to identify potential sources or processes that could explain their presence and distribution. Trace elements and radioactive constituents came from natural sources and were not elevated by anthropogenic sources or processes, except for selenium and the radioactive constituent uranium. Arsenic, manganese, iron, selenium, and uranium concentrations were all influenced by oxidation-reduction conditions. Unlike other trace elements and radioactive constituents, uranium and selenium can be affected by agricultural practices. Uranium and selenium can be released from aquifer sediments as a result of irrigation recharge water interacting with bicarbonate systems. Nitrate can be strongly affected by anthropogenic sources. Nitrate concentrations were significantly higher in modern groundwater, indicating recent inputs of nitrate to the shallow aquifer system. Nitrate was positively correlated with agricultural land use, indicating that irrigation-return water could be leaching nitrogen fertilizer and naturally present nitrate to elevate nitrate concentrations in shallow groundwater. The salinity indicators total dissolved solids, chloride, and sulfate all had natural sources in the MS-SA study unit, primarily marine sediments. Concentrations of the constituents were elevated as a result of evaporative concentration of irrigation water or precipitation. Sulfate concentrations were significantly correlated to agricultural land use, indicating that agricultural land-use practices are a contributing source of sulfate to groundwater. The samples with most of the detections of VOCs were from sites in the Pajaro Valley and northern part of the Salinas Valley. Most of the samples with pesticide detections were from sites in the Salinas Valley study area. The herbicide simazine was positively correlated to the percentage of agricultural land use, and its concentrations were higher in modern groundwater than in pre-modern groundwater. Perchlorate, similar to nitrate, has natural and anthropogenic sources. Correlations of perchlorate to dissolved oxygen, nitrate, and percentage of agricultural land use indicated that the irrigation-return water could be leaching naturally present perchlorate, as well as perchlorate from historical applications of Chilean nitrate fertilizer, to increase perchlorate concentrations in groundwater. The quality of the water in the shallow aquifer system from this study was compared with the quality of water in the public drinking-water aquifer in a previous GAMA (MS-PA) study in the same area. The shallow system was more oxic and had more sites with modern groundwater than the public drinking-water aquifer, which was more anoxic and had sites with more pre-modern groundwater. Arsenic and selenium were found at high relative concentrations in a greater proportion of the shallow system. Manganese and iron were found at high relative concentrations in a greater proportion of the public drinking-water aquifer. Uranium was found at higher relative concentrations in a greater proportion of the shallow system. Concentrations of arsenic, iron, manganese, and molybdenum are not likely to change much as groundwater percolates from the shallow system to the public drinking-water aquifer because there are no anthropogenic sources affecting these constituents. Uranium and selenium concentrations in the public drinking-water aquifer could be affected by the higher concentrations in the shallow system because of irrigation-return water, however. Nitrate and salinity indicators had concentrations that were much higher in the shallow system than the deeper public drinking-water aquifer. High concentrations of these constituents in the shallow system could lead to increased concentrations in the public drinking-water aquifer in parts of the study units because of land-use practices, such as irrigated agriculture. Organic constituents were detected more frequently in the public drinking-water aquifer than in the shallow system, possibly because more of the sites sampled in the public drinking-water aquifer were in urban areas compared to the sites sampled for the shallow system or because sources of contamination have decreased as a result of changes in use at the land surface.

California↗

Groundwater levels in the Denver Basin bedrock aquifers of Douglas County, Colorado, 2011–19

Municipal and domestic water users in Douglas County, Colorado, rely on groundwater from the bedrock aquifers in the Denver Basin aquifer system as part of their water supply. The four principal Denver Basin bedrock aquifers are, from shallowest to deepest, the Dawson aquifer (divided administratively into “upper” and “lower” Dawson aquifers in Douglas County), the Denver aquifer, the Arapahoe aquifer, and the Laramie-Fox Hills aquifer. Increased groundwater pumping in response to rapid population growth and development has led to declining groundwater levels in Douglas County, where groundwater is a primary water source for densely populated and rural communities. The U.S. Geological Survey, in cooperation with the Rural Water Authority of Douglas County, began a study in 2011 to assess the groundwater resources of the Denver Basin bedrock aquifers within the county. The primary purpose of this report is to present a summary of groundwater levels measured during the study period (2011–19) and present results from statistical analyses of changes in groundwater-level elevations, reported above the land-surface datum, North American Vertical Datum of 1988, through time. During the study period, January 2011 through June 2019, discrete groundwater levels were routinely measured at 36 wells producing from Denver Basin bedrock aquifers within Douglas County. Of the 36 wells, 15 are instrumented with pressure transducers that record groundwater-level measurements at hourly intervals, and these data were temporally aggregated into time-series records. During 2011, wells were added to the monitoring network in phases, so that the start dates of the well records are noncontemporaneous. To keep temporal analysis among wells consistent, the periods of record used in statistical analyses were from February 2012 through February 2019 for the discrete data and from January 2012 through June 2019 for the time-series data. The upper Dawson, lower Dawson, Denver, and Arapahoe aquifers had some wells with rises in calculated groundwater-level elevations, but most wells showed declines on the basis of statistically significant trends and the relative differences in static groundwater-level elevations between the February 2012 and February 2019 measurements. Neither of the two wells in the Laramie-Fox Hills aquifer showed significant trends in groundwater-level elevations, and these wells had few static discrete measurements, precluding a comparison between 2012 and 2019 static groundwater-level elevations. Of the 13 wells in the upper Dawson, lower Dawson, Denver, and Arapahoe aquifers with significant trends in discrete groundwater-level elevation measurements, the records of 12 wells demonstrated negative trends during the study period. The upper Dawson, lower Dawson, Denver, and Arapahoe aquifers had median significant trends of −0.23, −0.31, −0.92, and −2.26 feet per year, respectively. Although the Arapahoe aquifer had the greatest negative median trend, this median only represents one well with significant trends. Otherwise, the Denver aquifer had the next greatest negative trend, with a median trend of −0.92 foot per year. Significant trends in time-series groundwater-level elevations agreed with significant trends in discrete groundwater-level elevations; for all wells with statistically significant trends in discrete and in time-series groundwater-level elevation data, trend estimates from the two records were within 0.1 foot per year of each other. Potentiometric-surface maps of the upper Dawson, lower Dawson, and Denver aquifers, created using discrete static groundwater levels measured in February 2019, show that groundwater flow direction for the upper Dawson, lower Dawson, and Denver aquifers is generally from south to north. Results of this study could guide future groundwater monitoring in the county and aid in long-term planning of water resources.

Colorado↗

An update of hydrologic conditions and distribution of selected constituents in water, eastern Snake River aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2019–21

Since 1952, wastewater discharged to infiltration ponds (also called “percolation ponds”) and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain (ESRP) aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy (DOE), maintains groundwater-monitoring networks at the INL to determine hydrologic trends and to delineate the movement of radiochemical and chemical wastes in both the aquifer and perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from the ESRP aquifer and perched groundwater wells from the USGS groundwater monitoring networks during 2019–21. From March–May 2018 to March–May 2021, water levels in wells completed in the ESRP aquifer increased in the northern part of the INL and decreased in the southwestern part. Water-level increases ranged from 0.02 to 1.04 feet in the northern part and decreases ranged from 0.03 to 2.94 feet in the southwestern part of the INL. Detectable concentrations of radiochemical constituents in water samples from wells in the ESRP aquifer at the INL generally decreased or remained constant during 2019–21. Decreases in concentrations were attributed to radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In 2021, tritium was detected above reporting levels in water samples collected from 46 of 105 aquifer wells and ranged from 150±50 to 4,280±150 picocuries per liter (pCi/L). Tritium concentrations from eight wells completed in deep perched groundwater near the Advanced Test Reactor Complex (ATRC) generally were greater than or equal to the reporting level during at least one sampling event during 2019–21, and concentrations ranged from 160±50 to 2,097±107 pCi/L. Concentrations of strontium-90 in water from 12 of 45 aquifer wells sampled in 2021 exceeded the reporting level, and concentrations ranged from 2.5±0.7 to 299±6 pCi/L. During 2021, concentrations of strontium-90 from five wells completed in deep perched groundwater at the ATRC equaled or exceeded the reporting levels, and concentrations ranged from 3±0.9 pCi/L to 27.8±1.3 pCi/L. Concentrations of cesium-137 were less than the reporting level in all but one aquifer well, and concentrations of plutonium-238, plutonium-239, -240 (undivided), and americium-241 were less than the reporting level in water samples from all aquifer wells sampled during this study period. Dissolved chromium concentrations in water samples from 64 ESRP aquifer wells ranged from less than (<) 0.5 to 76.4 micrograms per liter (μg/L). During 2019–21, dissolved chromium was detected in water from wells completed in deep perched groundwater above the ESRP aquifer at the ATRC, and concentrations ranged from <1 to 82.1 μg/L. In 2021, concentrations of dissolved sodium in water from most ESRP aquifer wells in the southern part of the INL were greater than the western tributary groundwater background concentration of 8.3 milligrams per liter (mg/L). During 2021, dissolved sodium concentrations in water from 15 wells completed in deep perched groundwater ranged from 11.7 to 122.5 mg/L. Variations in sodium concentrations in aquifer wells and perched groundwater zones are attributed to either migration of remnant water from the former chemical-waste ponds or disposal volume and composition variability in percolation ponds installed in 2008. In 2021, concentrations of chloride in most water samples from ESRP aquifer wells south of the Idaho Nuclear Technology and Engineering Center (INTEC) and at the Central Facilities Area (CFA) exceeded background concentrations. Chloride concentrations in water from wells south of the INTEC have generally decreased because of discontinued chloride disposal to the legacy percolation ponds since 2002 when the discharge of wastewater was discontinued. During 2019–21, dissolved chloride concentrations in deep perched groundwater above the ESRP aquifer from 18 wells at the ATRC ranged from 8.15 to 231 mg/L. In 2021, sulfate concentrations in water samples from ESRP aquifer wells in the south-central part of the INL that exceeded the background concentration of sulfate, ranged from 21 to 141 mg/L. The greater-than-background concentrations in water from these wells are attributed to sulfate disposal at the ATRC infiltration ponds or the legacy INTEC percolation ponds. In 2021, sulfate concentrations in water samples from aquifer wells near the Radioactive Waste Management Complex (RWMC) were mostly greater than background concentrations. The maximum dissolved sulfate concentration in shallow perched groundwater near the ATRC was 575 mg/L in 2021. During 2021, dissolved sulfate concentrations in water from wells completed in deep perched groundwater near the cold waste ponds at the ATRC ranged from 22.3 to 519 mg/L. In 2021, concentrations of nitrate in water from most ESRP aquifer wells at and near the INTEC exceeded the western tributary groundwater background concentration of 0.655 mg/L. Concentrations of nitrate in aquifer wells southwest of INTEC and farther away from the influence of disposal areas and the Big Lost River, in intermittent source of surface water recharge to the aquifer, show a general decrease in nitrate concentration over time. Two aquifer wells south of INTEC show increasing trends that could result from wastewater beneath the INTEC tank farm being mobilized to the aquifer. During 2019–21, water samples from several ESRP aquifer wells were collected and analyzed for volatile organic compounds (VOCs). Twelve VOCs were detected, and 1–4 VOCs were detected in water samples from 10 wells. The most frequently detected VOCs include carbon tetrachloride (tetrachloromethane), trichloromethane, tetrachloroethene, 1,1,1-trichloroethane, and trichloroethene. In 2019–21, concentrations for all VOCs were less than their respective maximum contaminant levels (MCLs) for drinking water, except carbon tetrachloride in one well, trichloroethene in two wells, and vinyl chloride in one well. During 2019–21, variability and bias were evaluated from 34 replicate and 14 blank quality-assurance samples. Results from replicate analyses were investigated to evaluate sample variability. Constituents with acceptable reproducibility were major ions, trace elements, nutrients, and VOCs. All radiochemical constituents including gross alpha- and beta- radioactivity, strontium-90, cesium-137, and tritium, had acceptable reproducibility. Bias from sample contamination was evaluated from equipment, field, and source-solution blanks. Chloride and sulfate were detected slightly above their respective method detection limits in equipment and field blanks, but at concentrations well below the co-collected sample for that well. These chloride and sulfate detections in the field and equipment blanks were inconsequential because they weren’t detected above the analysis-specific variability for those constituents as determined by replicate sample result evaluation. None of the detections of nutrients and trace inorganic constituents were high enough to indicate environmental sample or analytical procedure bias.

Idaho↗

Geohydrology, water quality, and simulation of groundwater flow in the stratified-drift aquifer system in Virgil Creek and Dryden Lake Valleys, Town of Dryden, Tompkins County, New York

In 2002, the U.S. Geological Survey, in cooperation with the Tompkins County Planning Department and the Town of Dryden, New York, began a study of the stratified-drift aquifer system in the Virgil Creek and Dryden Lake Valleys in the Town of Dryden, Tompkins County. The study provided geohydrologic data needed by the town and county to develop a strategy to manage and protect their water resources. In this study area, three extensive confined sand and gravel aquifers (the upper, middle, and lower confined aquifers) compose the stratified-drift aquifer system. The Dryden Lake Valley is a glaciated valley oriented parallel to the direction of ice movement. Erosion by ice extensively widened and deepened the valley, truncated bedrock hillsides, and formed a nearly straight, U-shaped bedrock trough. The maximum thickness of the valley fill in the central part of the valley is about 400 feet (ft). The Virgil Creek Valley in the east part of the study area underwent less severe erosion by ice than the Dryden Lake Valley, and hence, it has a bedrock floor that is several hundred feet higher in altitude than that in the Dryden Lake Valley. The sources and amounts of recharge were difficult to identify in most areas because the confined aquifers are overlain by confining units. However, in the vicinity of the Virgil Creek Dam, the upper confined aquifer crops out at land surface in the floodplain of a gorge eroded by Virgil Creek, and this is where the aquifer receives large amounts of recharge from precipitation that directly falls over the aquifer and from seepage losses from Virgil Creek. The results of streamflow measurements made in Virgil Creek where it flows through the gorge indicated that the stream lost 1.2 cubic feet per second (ft 3 /s) or 0.78 million gallons per day (Mgal/d) of water in the reach extending from 220 ft downstream from the dam to 1,200 ft upstream from the dam. In the southern part of the study area, large amounts of recharge also replenish the stratified-drift aquifers at the Valley Heads Moraine, which consists of heterogeneous sediments including coarse-grained outwash and kame sediments, as well as zones containing till with a fine-grained matrix. In the southern part of the study area, the confining units are thin and likely to be discontinuous in some places, resulting in windows of permeable sediment, which can more readily transmit recharge from precipitation and from tributaries that lose water as they flow over the valley floor. In contrast, in the northern part of the study area, the confining units are thick, continuous, and comprise homogeneous fine-grained sediments that more effectively confine the aquifers than in the southern part of the study area. Most groundwater in the northern part of the study area discharges to the Village of Dryden municipal production wells, to the outlet to Dryden Lake, to Virgil Creek, and as groundwater underflow that exits the northern boundary of the study area. Most northward-flowing groundwater in the southern part of the study area discharges to Dryden Lake, to the inlet to Dryden Lake, and to homeowner, nonmunicipal community (a mobile home community and several apartments), and commercial wells. Most of this pumped water is returned to the groundwater system via septic systems. Most southward-flowing groundwater in the southern part of the study area discharges to the headwaters of Owego Creek and to agricultural wells; some flow also exits the southern boundary of the study area as groundwater underflow. The largest user of groundwater in the study area is the Village of Dryden. Water use in the village has approximately tripled between the early 1970s when withdrawals ranged between 18 and 30 million gallons per year (Mgal/yr) and from 2000 through 2008 when withdrawals ranged between 75 and 85 Mgal/yr. The estimated groundwater use by homeowners, nonmunicipal communities, and small commercial facilities outside the area supplied by the Village of Dryden municipal wells is estimated to be about 18.4 Mgal/yr. Most of this pumped water is returned to the groundwater system via septic systems. For this investigation, an aquifer test was conducted at the Village of Dryden production well TM 981 (finished in the middle confined aquifer at a well depth of 72 ft) at the Jay Street pumping station during June 19&ndash;21, 2007. The aquifer test consisted of pumping production well TM 981 at 104 gallons per minute over a 24-hour period. The drawdown in well TM 981 at the end of 24 hours of pumping was 19.2 ft. Results of the aquifer-test analysis for a partially penetrating well in a confined aquifer indicated that the transmissivity was 1,560 feet squared per day, and the horizontal hydraulic conductivity was 87 feet per day, based on a saturated thickness of 18 ft. During 2003&ndash;5, 14 surface-water samples were collected at 8 sites, including Virgil Creek, Dryden Lake outlet, and several tributaries. During 2003 through 2009, eight groundwater samples were collected from eight wells, including three municipal production wells, two test wells, and three domestic wells. Calcium dominates the cation composition, and bicarbonate dominates the anion composition in most groundwater and surface-water samples. None of the common inorganic constituents collected exceeded any Federal or State water-quality standards. Results from a three-dimensional, finite-difference groundwater-flow model were used to compute a water budget and to estimate the areal extent of the zone of groundwater contribution to the Village of Dryden municipal production wells. The model-computed water budget indicated that the sources of recharge to the confined aquifer system are precipitation that falls directly on the valley-fill sediments (40 percent of total recharge), stream leakage (35.5 percent), seepage from wetlands and ponds (12 percent), unchanneled runoff and groundwater inflow from the uplands (8.5 percent), and groundwater underflow into the eastern end of the model area (4 percent). Most groundwater discharges to surface-water bodies, including Dryden Lake (33 percent), streams (33 percent), and wetlands and ponds (10 percent of the total). In addition, some groundwater discharges as underflow out of the southern and northern ends of the model area (15 percent), to simulated pumping wells (4.5 percent), and to drains that represent seepage from the bluffs exposed in the gorge in the vicinity of the Virgil Creek Dam (4.5 percent). The areal extents of the zones of groundwater contribution for Village of Dryden municipal production wells TM 202 (Lake Road pump station, finished in the upper confined aquifer) and TM 981 (Jay Street pump station, finished in the middle confined aquifer) are 0.5 square mile (mi 2 ) and 0.9 mi 2 , respectively. The areal extent of the zone of contribution to production well TM 202 extends 2.2 miles (mi) southeast into the Virgil Creek Valley, whereas production well TM 981 extends 3.8 mi south in the Dryden Lake Valley. The areal extent of the zone of contribution to production well TM1046 (South Street pump station) is 1.4 mi 2 and extends 2.4 mi into Dryden Lake Valley and 0.5 mi into Virgil Creek Valley.

New York↗