Search USGS⌕ Search

Geology topics

Derek W. Ryter

Publications and source records attributed to Derek W. Ryter.

16 recordsLinked to original sources
USGS2026

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the effects of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study, we constructed the Strawberry Creek integrated hydrological model to increase this understanding. The Strawberry Creek integrated hydrological model is an integrated surface runoff and groundwater model, GSFLOW to simulate surface runoff and infiltration and groundwater flow. The model has 263 rows and 176 columns in each of its three layers. The model area includes the Strawberry Creek and four adjacent watersheds. The precipitation-runoff modeling system was calibrated using two streamflow gaging stations, and the groundwater and surface-water flow model was calibrated to the reported spring diversion discharge and a sparse number of groundwater-level measurements. The Strawberry Creek integrated hydrological model was run with and without diversions active, and simulated streamflow was compared, indicating that in the headwaters of Strawberry Creek, about 35 percent of the diversion flow was captured from base flow.

California↗
USGS2026

Assessing future hydrologic extremes using an integrated hydrology and river operations model in the Russian River watershed

Study region The Russian River watershed, situated in coastal, northern California, experiences hydrologic extremes, including periodic droughts and flooding. Water managers are working to maintain sustainable water supplies and environmental flows, while mitigating flood risks. Study focus This paper introduces an integrated hydrology and river operations model for the Russian River watershed. This model is distinct from models in previous studies because it represents surface-groundwater interactions and uses climate forcings to estimate dynamic water use demands that are superimposed onto both reservoir operations and water supply constraints. The model was used to examine three historical (1990–2015) and eight future (2016–2099) water use and climate change scenarios. New hydrological insights for the region The direct connection between streams and aquifers facilitated both annual aquifer replenishment by high winter streamflows and streamflow depletion by groundwater wells (19 % of pumped groundwater in alluvial aquifers from stream leakage) during critical low flow periods. Simulated streamflow changes included 59 % longer and 54 % more severe streamflow droughts, 26 % lower seasonal low streamflows, and up to 125 % higher peak streamflows, averaged over future climate and water use scenarios, suggesting increased future flood and water availability risks. Results showed the importance of reservoir operations for mitigating the impacts of increased hydroclimatic volatility, despite a decrease in reservoir reliability at Lake Mendocino, suggesting that reservoir management may be used to decrease future risks.

California↗
USGS2025

Updating and recalibrating the integrated Santa Rosa Plain Hydrologic Model to assess stream depletion and to simulate future climate and management scenarios in Santa Rosa, Sonoma County, California

The Santa Rosa Plain Hydrologic Model (SRPHM) was developed and published in 2014 through a collaboration between the U.S. Geological Survey (USGS) and Sonoma Water to analyze the hydrologic system in the Santa Rosa Plain watershed, help meet the increasing demand for fresh water, and prepare for future uncertainties in water resources. The original model simulated hydrological conditions and water use from water years 1975 to 2010. Recently (2023), the USGS, in cooperation with Sonoma Water and the California State Water Resources Control Board, updated the SRPHM model to extend its simulation period to the end of the 2018 calendar year, incorporate new estimates of rural and agricultural water use, and use efficient input format for climate variables. The updated model was recalibrated, and evaluation of the new model calibration is included in this report. This report presents the results of comparing the hydraulic heads, streamflow, and groundwater budget simulated by the updated model with those generated by the original model and observed data. The main difference in the simulated budget between the original and updated SRPHM is the estimates of agricultural pumping, rural domestic pumping, and return flow generated from rural water use that was not simulated in the original model. The revised agricultural pumping is simulated using the agricultural package, which constrains pumping to available groundwater. The use of the agricultural package leads to a more realistic estimation of agricultural water use, with revised agricultural pumping being one-third less than that in the original model. The revised rural pumping is about half of the pumping in the original model because of using detailed parcel data to estimate population density in rural areas instead of coarse census tracts. Overall, average total inflows for water years 2006–10 simulated by the updated model were about 2 percent less than the original model, and the average total updated outflows were nearly 5 percent less than the original model. The updated model was then used to generate stream depletion maps, simulate climate change scenarios during 2019–99, and simulate water rights allocation using the Model for Decision Support in Integrated River Basin Management (MODSIM). The results from simulating eight future climate scenarios indicated either an increase in groundwater storage or no significant change in the next 80 years, along with an increase in recharge, an increase in actual evapotranspiration in six out of eight climate projections, and an increase in surface runoff. The increases in the simulated future groundwater storage, recharge, evapotranspiration, and runoff in most climate projections are mainly driven by the projected increase in precipitation in most of the future climate scenarios. The updated model also was used to test a pilot case study demonstrating water-resource allocation among different users with different water rights using the integrated MODSIM-Groundwater and Surface-Water Flow Model (GSFLOW) platform. The updated SRPHM serves as a valuable tool for analyzing historical and future hydrologic conditions in the Santa Rosa Plain watershed and preparing for future uncertainties.

California↗
USGS2025

Simulation of the impacts of projected climate change on groundwater resources in the urban, semiarid Yucaipa Valley watershed, southern California using an integrated hydrologic model

Managing water resources in semiarid watersheds is challenging due to limited supply and uncertain future climate conditions. This paper examines the impact of future climate changes on an urban watershed in southern California using an integrated hydrologic model. GSFLOW modeling software is used to simulate the nonlinear relationships between climate trends and precipitation partitioning into ET, runoff, and subsurface storage. Four global circulation models (GCMs), each with two greenhouse-gas scenarios, RCP45 and RCP85 are used to project future climate conditions. GCMs include the CanESM2, CNRM-CM5, HadGEM2-ES, and MIROC5 models. The model's simulated hydrologic conditions are compared with historical data to assess changes in water budgets and groundwater supply. Results indicate decreased groundwater storage in most scenarios due to increased natural evapotranspiration, vegetation consumptive use, and streamflow out of the watershed. Only scenarios with substantially increased future precipitation show increased groundwater storage. The study also highlights increased future aridity despite the rise in precipitation and large precipitation events forecast by GCMs, which increase the risk of urban floods and decrease stream leakage and water available to vegetation.

California↗
USGS2025

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv↗
USGS2020

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↗
USGS2016

Summary of U.S. Geological Survey studies conducted in cooperation with the Citizen Potawatomi Nation, central Oklahoma, 2011–14

The U.S. Geological Survey conducted hydrologic studies and published three U.S. Geological Survey scientific investigations reports in cooperation with the Citizen Potawatomi Nation from 2011 to 2014 to characterize the quality and quantity of water resources. The study areas of those reports consisted of approximately 960 square miles in parts of three counties in central Oklahoma. This study area has multiple abundant sources of water, being underlain by three principal aquifers (alluvial/terrace, Central Oklahoma, and Vamoosa-Ada), being bordered by two major rivers (North Canadian and Canadian), and having several smaller drainages including the Little River in the central part of the study area and Salt Creek in the southeastern part of the study area. The Central Oklahoma aquifer (also referred to as the “Garber-Wellington aquifer”) underlies approximately 3,000 square miles in central Oklahoma in parts of Cleveland, Logan, Lincoln, Oklahoma, and Pottawatomie Counties and much of the study area. Water from these aquifers is used for municipal, industrial, commercial, agricultural, and domestic supplies. Much of the water in the study area is of good quality; however, in some parts of this area water quality was impaired by very hard surface water and groundwater, large chloride concentrations in some smaller streams, relatively large concentrations of nitrogen and phosphorus nutrients and large counts of fecal-indicator bacteria in the North Canadian River, and uranium concentrations that exceeded the U.S. Environmental Protection Agency Maximum Contaminant Level of 30 micrograms per liter for public water supplies in water samples collected from a small number of wells. Most stream-water samples collected from the Little River by the U.S. Geological Survey in 2012–13 had dissolved solids concentrations exceeding the U.S. Environmental Protection Agency Secondary Maximum Contaminant Level for public water supplies of 500 milligrams per liter. Larger numbers of organic compounds were measured in water samples collected from the North Canadian River than the Little River. Numerical groundwater-flow models were created to characterize flow systems in aquifers underlying this study area and areas of particular interest within the study area. Those models were used to estimate sustainable groundwater yields from parts of the North Canadian River alluvial aquifer, characterize groundwater/surface-water interactions, and estimate the effects of a 10-year simulated drought on streamflows and water levels in alluvial and bedrock aquifers. Pumping of wells at the Iron Horse Industrial Park was estimated to cause negligible infiltration of water from the adjoining North Canadian River. A 10-year simulated drought of 50 percent of normal recharge was tested for the period 1990–2000. For this period, the total amount of groundwater in storage was estimated to decrease by 8.6 percent in the North Canadian River alluvial aquifer and approximately 0.2 percent in the Central Oklahoma aquifer, and groundwater flow to streams was estimated to decrease by 28–37 percent. This volume of groundwater loss showed that the Central Oklahoma aquifer is a bedrock aquifer that has relatively low rates of recharge from the land surface. The simulated drought decreased simulated streamflow, composed of base flow, in the North Canadian River at Shawnee, Okla., which did not recover to predrought conditions until the relatively wet year of 2007 after the simulated drought period.

Oklahoma↗
USGS2016

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↗
USGS2015

Numerical simulation of groundwater flow, resource optimization, and potential effects of prolonged drought for the Citizen Potawatomi Nation Tribal Jurisdictional Area, central Oklahoma

A hydrogeological study including two numerical groundwater-flow models was completed for the Citizen Potawatomi Nation Tribal Jurisdictional Area of central Oklahoma. One numerical groundwater-flow model, the Citizen Potawatomi Nation model, encompassed the jurisdictional area and was based on the results of a regional-scale hydrogeological study and numerical groundwater flow model of the Central Oklahoma aquifer, which had a geographic extent that included the Citizen Potawatomi Nation Tribal Jurisdictional Area. The Citizen Potawatomi Nation numerical groundwater-flow model included alluvial aquifers not in the original model and improved calibration using automated parameter-estimation techniques. The Citizen Potawatomi Nation numerical groundwater-flow model was used to analyze the groundwater-flow system and the effects of drought on the volume of groundwater in storage and streamflow in the North Canadian River. A more detailed, local-scale inset model was constructed from the Citizen Potawatomi Nation model to estimate available groundwater resources for two Citizen Potawatomi Nation economic development zones near the North Canadian River, the geothermal supply area and the Iron Horse Industrial Park. Groundwater pumping rates at potential well locations were optimized using the most recent version of the U.S. Geological Survey Groundwater-Management Process for MODFLOW. The objectives of optimization were to determine if a total pumping rate of 500 gallons per minute could be pumped from 5 wells at the geothermal supply area and to maximize discharge from 16 wells at the Iron Horse Industrial Park without exceeding specified head drawdown constraints at the pumping wells and thus prevent groundwater depletion. The inset model was used to estimate North Canadian River streamflow depletion caused by optimized pumping at the Iron Horse Industrial Park because water quality was a concern, and the river may have degraded water quality compared to water in other parts of the alluvial aquifer. The fate of streamflow that infiltrates into groundwater because of pumping was not directly determined, but it was assumed that this water could end up in the well discharge, and was considered to be a maximum proportion of well discharge derived from the North Canadian River. The total optimized continuous pumping rate from five managed wells at the geothermal supply area was 638 gallons per minute, which exceeded the target pumping rate of 500 gallons per minute. The total continuous pumping rate from 16 wells at the Iron Horse Industrial Park was 1,472 gallons per minute, which induced stream infiltration of approximately 4.1 gallons per minute (approximately 0.3 percent of the total well discharge) from the North Canadian River. To estimate the effects of drought on water resources in the Citizen Potawatomi Nation Tribal Jurisdictional Area, a hypothetical 10-year drought during which precipitation would decrease by 50 percent was simulated by decreasing model groundwater recharge by the same proportion for the period 1990–2000 of the transient model. The effects of the drought were estimated by calculating the change in the volume of groundwater storage and groundwater flow to streams at the end of the drought period, and the change in simulated streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., during and after the drought. The hypothetical decrease in recharge during the simulated drought caused groundwater in storage over the entire model in the study area to decrease by 361,500 acre-feet (14,100 acre-feet in the North Canadian River alluvial aquifer and 347,400 acre-feet in the Central Oklahoma aquifer), or approximately 0.2 percent of the total groundwater in storage over the drought period. This small percentage of groundwater loss showed that the Central Oklahoma aquifer as a bedrock aquifer has relatively low rates of recharge from the surface relative to the approximate storage. The budget for base flow to the North Canadian River indicated that the change in groundwater flow to the North Canadian River decreased during the 10-year drought by 386,500 acre-feet, or 37 percent. In all other parts of the Citizen Potawatomi Nation Tribal Jurisdictional Area, base flow decreased by 292,000 acre-feet, or 28 percent. Streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., decreased during the hypothetical drought by as much as 28 percent, and the mean change in streamflow decreased as much as 16 percent. Streamflow at the Shawnee streamflow-gaging station did not recover to nondrought conditions until about 3 years after the simulated drought ended, during the relatively wet year of 2007.

Oklahoma↗
USGS2015

Simulation of groundwater flow and analysis of the effects of water-management options in the North Platte Natural Resources District, Nebraska

The North Platte Natural Resources District (NPNRD) has been actively collecting data and studying groundwater resources because of concerns about the future availability of the highly inter-connected surface-water and groundwater resources. This report, prepared by the U.S. Geological Survey in cooperation with the North Platte Natural Resources District, describes a groundwater-flow model of the North Platte River valley from Bridgeport, Nebraska, extending west to 6 miles into Wyoming. The model was built to improve the understanding of the interaction of surface-water and groundwater resources, and as an optimization tool, the model is able to analyze the effects of water-management options on the simulated stream base flow of the North Platte River. The groundwater system and related sources and sinks of water were simulated using a newton formulation of the U.S. Geological Survey modular three-dimensional groundwater model, referred to as MODFLOW–NWT, which provided an improved ability to solve nonlinear unconfined aquifer simulations with wetting and drying of cells. Using previously published aquifer-base-altitude contours in conjunction with newer test-hole and geophysical data, a new base-of-aquifer altitude map was generated because of the strong effect of the aquifer-base topography on groundwater-flow direction and magnitude. The largest inflow to groundwater is recharge originating from water leaking from canals, which is much larger than recharge originating from infiltration of precipitation. The largest component of groundwater discharge from the study area is to the North Platte River and its tributaries, with smaller amounts of discharge to evapotranspiration and groundwater withdrawals for irrigation. Recharge from infiltration of precipitation was estimated with a daily soil-water-balance model. Annual recharge from canal seepage was estimated using available records from the Bureau of Reclamation and then modified with canal-seepage potentials estimated using geophysical data. Groundwater withdrawals were estimated using land-cover data, precipitation data, and published crop water-use data. For fields irrigated with surface water and groundwater, surface-water deliveries were subtracted from the estimated net irrigation requirement, and groundwater withdrawal was assumed to be equal to any demand unmet by surface water. The groundwater-flow model was calibrated to measured groundwater levels and stream base flows estimated using the base-flow index method. The model was calibrated through automated adjustments using statistical techniques through parameter estimation using the parameter estimation suite of software (PEST). PEST was used to adjust 273 parameters, grouped as hydraulic conductivity of the aquifer, spatial multipliers to recharge, temporal multipliers to recharge, and two specific recharge parameters. Base flow of the North Platte River at Bridgeport, Nebraska, streamgage near the eastern, downstream end of the model was one of the primary calibration targets. Simulated base flow reasonably matched estimated base flow for this streamgage during 1950–2008, with an average difference of 15 percent. Overall, 1950–2008 simulated base flow followed the trend of the estimated base flow reasonably well, in cases with generally increasing or decreasing base flow from the start of the simulation to the end. Simulated base flow also matched estimated base flow reasonably well for most of the North Platte River tributaries with estimated base flow. Average simulated groundwater budgets during 1989–2008 were nearly three times larger for irrigation seasons than for non-irrigation seasons. The calibrated groundwater-flow model was used with the Groundwater-Management Process for the 2005 version of the U.S. Geological Survey modular three-dimensional groundwater model, MODFLOW–2005, to provide a tool for the NPNRD to better understand how water-management decisions could affect stream base flows of the North Platte River at Bridgeport, Nebr., streamgage in a future period from 2008 to 2019 under varying climatic conditions. The simulation-optimization model was constructed to analyze the maximum increase in simulated stream base flow that could be obtained with the minimum amount of reductions in groundwater withdrawals for irrigation. A second analysis extended the first to analyze the simulated base-flow benefit of groundwater withdrawals along with application of intentional recharge, that is, water from canals being released into rangeland areas with sandy soils. With optimized groundwater withdrawals and intentional recharge, the maximum simulated stream base flow was 15–23 cubic feet per second (ft 3 /s) greater than with no management at all, or 10–15 ft 3 /s larger than with managed groundwater withdrawals only. These results indicate not only the amount that simulated stream base flow can be increased by these management options, but also the locations where the management options provide the most or least benefit to the simulated stream base flow. For the analyses in this report, simulated base flow was best optimized by reductions in groundwater withdrawals north of the North Platte River and in the western half of the area. Intentional recharge sites selected by the optimization had a complex distribution but were more likely to be closer to the North Platte River or its tributaries. Future users of the simulation-optimization model will be able to modify the input files as to type, location, and timing of constraints, decision variables of groundwater withdrawals by zone, and other variables to explore other feasible management scenarios that may yield different increases in simulated future base flow of the North Platte River.

Nebraska↗
USGS2014

Hydrogeology and simulation of groundwater flow in the Central Oklahoma (Garber-Wellington) Aquifer, Oklahoma, 1987 to 2009, and simulation of available water in storage, 2010–2059

The Central Oklahoma (Garber-Wellington) aquifer underlies about 3,000 square miles of central Oklahoma. The study area for this investigation was the extent of the Central Oklahoma aquifer. Water from the Central Oklahoma aquifer is used for public, industrial, commercial, agricultural, and domestic supply. With the exception of Oklahoma City, all of the major communities in central Oklahoma rely either solely or partly on groundwater from this aquifer. The Oklahoma City metropolitan area, incorporating parts of Canadian, Cleveland, Grady, Lincoln, Logan, McClain, and Oklahoma Counties, has a population of approximately 1.2 million people. As areas are developed for groundwater supply, increased groundwater withdrawals may result in decreases in long-term aquifer storage. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, investigated the hydrogeology and simulated groundwater flow in the aquifer using a numerical groundwater-flow model. The purpose of this report is to describe an investigation of the Central Oklahoma aquifer that included analyses of the hydrogeology, hydrogeologic framework of the aquifer, and construction of a numerical groundwater-flow model. The groundwater-flow model was used to simulate groundwater levels and for water-budget analysis. A calibrated transient model was used to evaluate changes in groundwater storage associated with increased future water demands.

Oklahoma↗
USGS2014

Hydrogeologic framework of the uppermost principal aquifer systems in the Williston and Powder River structural basins, United States and Canada

The glacial, lower Tertiary, and Upper Cretaceous aquifer systems in the Williston and Powder River structural basins within the United States and Canada are the uppermost principal aquifer systems and most accessible sources of groundwater for these energy-producing basins. The glacial aquifer system covers the northeastern part of the Williston structural basin. The lower Tertiary and Upper Cretaceous aquifer systems are present in about 91,300 square miles (mi 2 ) of the Williston structural basin and about 25,500 mi 2 of the Powder River structural basin. Directly under these aquifer systems are 800 to more than 3,000 feet (ft) of relatively impermeable marine shale that serves as a basal confining unit. The aquifer systems in the Williston structural basin have a shallow (less than 2,900 ft deep), wide, and generally symmetrical bowl shape. The aquifer systems in the Powder River structural basin have a very deep (as much as 8,500 ft deep), narrow, and asymmetrical shape. The Williston structural basin has been an important oil and natural gas producing region since the 1950s, and production has increased substantially since the mid-2000s due to improved drilling and hydraulic fracturing methods from deep formations, such as the Bakken and Three Forks Formations. These improved methods require considerable volumes of freshwater mostly from shallow aquifers or surface water. Coal, lignite, and coal-bed natural gas are additional sources of energy in both basins that can affect the quality and quantity of shallow aquifers through strip mining and groundwater depletion. In 2011, the U.S. Geological Survey initiated a regional study of the glacial, lower Tertiary, and Upper Cretaceous aquifer systems in the Williston and Powder River structural basins with the goal to quantify groundwater availability. This report, together with a companion report of the conceptual flow model, provides an improved understanding of the groundwater flow systems and a basis for a numerical, regional groundwater-flow model. This study combines the lithostratigraphic units of the glacial, lower Tertiary, and Upper Cretaceous aquifer systems in the United States and Canada into 7 regional hydrogeologic units—glacial deposits, 4 bedrock aquifers, and 2 bedrock confining units—using general hydraulic properties. The glacial deposits are composed of till and glacial outwash sands and gravels with areas of cobbles and boulders. The four bedrock aquifers are the upper Fort Union, lower Fort Union, lower Hell Creek, and Fox Hills aquifers and are contained primarily in sandstone layers. The two confining units are the middle Fort Union hydrogeologic unit (shale) and upper Hell Creek hydrogeologic unit (contains less sandstone than the underlying lower Hell Creek aquifer). Water from hydrogeologic units in these three aquifer systems is relatively fresh and potable, whereas withdrawals seldom occur from units below the basal confining unit because of great depths (greater than 800 ft) and poor water quality. Analysis of about 300 electric (resistivity) and lithologic logs in the Williston structural basin and numerous existing publications for the Powder River structural basin were used to develop a three-dimensional hydrogeologic framework for both basins. Interpolated thicknesses of the glacial deposits, the lower Tertiary aquifer system, and the Upper Cretaceous aquifer system in the Williston structural basin are less than about 750; 2,250; and 1,050 ft, respectively. Interpolated thicknesses of the lower Tertiary aquifer system and the Upper Cretaceous aquifer system in the Powder River structural basin are less than about 7,180 and 5,070 ft, respectively. Interpolated horizontal hydraulic conductivity values for the Williston structural basin were as much as 25 feet per day (ft/d) in the glacial deposits and had smaller ranges in the lower Tertiary aquifer system (0.01–9.8 ft/d) and in the Upper Cretaceous aquifer system (0.06–5.5 ft/d). In the Powder River structural basin, the lower Tertiary aquifer system had a greater range of interpolated horizontal hydraulic conductivity values (0.10–11 ft/d) than the Upper Cretaceous aquifer system (0.02–5.7 ft/d). Transmissivity is greatest in the gravel zones of the glacial deposits (2,120 feet squared per day) and generally decreases with depth into the bedrock units. Regionally, water in the lower Tertiary and Upper Cretaceous aquifer systems flows in a northerly or northeasterly direction from the Powder River structural basin to the Williston structural basin. Groundwater flow in the Williston structural basin generally is easterly or northeasterly. Flow in the uppermost hydrogeologic units generally is more local and controlled by topography where unglaciated in the Williston structural basin than is flow in the glaciated part and in underlying aquifers. Groundwater flow in the Powder River structural basin generally is northerly with local variations greatest in the uppermost aquifers. Groundwater is confined, and flow is regional in the underlying aquifers.

Powder River Basin;Williston Basin↗
USGS2013

Geodatabase compilation of hydrogeologic, remote sensing, and water-budget-component data for the High Plains aquifer, 2011

The High Plains aquifer underlies almost 112 million acres in the central United States. It is one of the largest aquifers in the Nation in terms of annual groundwater withdrawals and provides drinking water for 2.3 million people. The High Plains aquifer has gained national and international attention as a highly stressed groundwater supply primarily because it has been appreciably depleted in some areas. The U.S. Geological Survey has an active program to monitor the changes in groundwater levels for the High Plains aquifer and has documented substantial water-level changes since predevelopment: the High Plains Groundwater Availability Study is part of a series of regional groundwater availability studies conducted to evaluate the availability and sustainability of major aquifers across the Nation. The goals of the regional groundwater studies are to quantify current groundwater resources in an aquifer system, evaluate how these resources have changed over time, and provide tools to better understand a systems response to future demands and environmental stresses. The purpose of this report is to present selected data developed and synthesized for the High Plains aquifer as part of the High Plains Groundwater Availability Study. The High Plains Groundwater Availability Study includes the development of a water-budget-component analysis for the High Plains completed in 2011 and development of a groundwater-flow model for the northern High Plains aquifer. Both of these tasks require large amounts of data about the High Plains aquifer. Data pertaining to the High Plains aquifer were collected, synthesized, and then organized into digital data containers called geodatabases. There are 8 geodatabases, 1 file geodatabase and 7 personal geodatabases, that have been grouped in three categories: hydrogeologic data, remote sensing data, and water-budget-component data. The hydrogeologic data pertaining to the northern High Plains aquifer is included in three separate geodatabases: (1) base data from a groundwater-flow model; (2) hydrogeology and hydraulic properties data; and (3) groundwater-flow model data to be used as calibration targets. The remote sensing data for this study were developed by the U. S. Geological Survey Earth Resources Observation and Science Center and include historical and predicted land-use/land-cover data and actual evapotranspiration data by using remotely sensed temperature data. The water-budget-component data contains selected raster data from maps in the “Selected Approaches to Estimate Water-Budget Components of the High Plains, 1940 Through 1949 and 2000 Through 2009” report completed in 2011 ( http://pubs.usgs.gov/sir/2011/5183/ ). Federal Geographic Data Committee compliant metadata were created for each spatial and tabular data layer in the geodatabases.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2012

Effects of linking a soil-water-balance model with a groundwater-flow model

A previously published regional groundwater-flow model in north-central Nebraska was sequentially linked with the recently developed soil-water-balance (SWB) model to analyze effects to groundwater-flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater-level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root-mean-squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB-generated recharge, the RMS difference between simulated and estimated base-flow target values for the groundwater-flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater-level and base-flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.

Nebraska↗
USGS2012

Altitude, age, and quality of groundwater, Papio-Missouri River Natural Resources District, eastern Nebraska, 1992 to 2009

The U.S. Geological Survey, in cooperation with the Papio-Missouri River Natural Resources District (PMRNRD), conducted this study to map the water-level altitude of 2009 within the Elkhorn River Valley, Missouri River Valley, and Platte River Valley alluvial aquifers; to present the predevelopment potentiometric-surface altitude within the Dakota aquifer; and to describe the age and quality of groundwater in the five principal aquifers of the PMRNRD in eastern Nebraska using data collected from 1992 to 2009. In addition, implications of alternatives to the current PMRNRD groundwater-quality monitoring approach are discussed. In the PMRNRD, groundwater altitude, relative to National Geodetic Vertical Datum of 1929, ranged from about 1,080 feet (ft) to 1,180 ft in the Elkhorn River Valley alluvial aquifer and from about 960 ft to 1,080 ft in the Missouri River Valley and Platte River Valley alluvial aquifers. In the PMRNRD, the estimated altitude of the potentiometric surface of the Dakota aquifer, predevelopment, ranged from about 1,100 ft to 1,200 ft. To assess groundwater age and quality, groundwater samples were collected from a total of 217 wells from 1992 to 2009 for analysis of various analytes. Groundwater samples collected in the PMRNRD from 1992 to 2009 and interpreted in this report were analyzed for age-dating analytes (chlorofluorocarbons), dissolved gases, major ions, trace elements, nutrients, stable isotope ratios, pesticides and pesticide degradates, volatile organic compounds, explosives, and 222radon. Apparent groundwater age was estimated from concentrations of chlorofluorocarbons measured in samples collected in 2000. Apparent groundwater-recharge dates ranged from older than 1940 in samples from wells screened in the Missouri River Valley alluvial aquifer to the early 1980s in samples from wells screened in the Dakota aquifer. Concentrations of major ions in the most recent sample per well collected from 1992 to 2009 indicate that the predominant water type was calcium bicarbonate. Samples from 4 wells exceeded the U.S. Environmental Protection Agency (USEPA) Secondary Drinking Water Regulation (SDWR) for sulfate [250 milligrams per liter (mg/L)], and samples from 4 wells exceeded the USEPA Drinking Water Advisory Table for sodium (30-60 mg/L). Eighteen of the 21 trace elements analyzed in samples from PMRNRD wells have USEPA drinking-water standards. Sixteen of the trace elements with USEPA standards were detected in the selected samples. In the samples selected for trace-element analysis, the only trace-element concentration that exceeded an enforceable USEPA drinking-water standard, the Maximum Contaminant Level (MCL), was for arsenic; arsenic concentration exceeded the USEPA MCL of 10 micrograms per liter (μg/L) in 4 percent of the samples. Trace-element concentrations that exceeded the USEPA SDWR or Lifetime Health Advisory level were iron (46 percent of the samples were greater than USEPA SDWR of 300 μg/L), manganese (70 percent of the samples were greater than USEPA SDWR of 50 μg/L), and strontium (4 percent of the samples were greater than USEPA Lifetime Health Advisory level of 4,000 μ/L). The concentration of nitrate plus nitrite as nitrogen (nitrate-N) in the most recent nutrient samples collected from the network wells and from one randomly selected well in the well nests from 1992 to 2009 for most wells (80 percent) ranged from less than 0.06 to 8.55 mg/L, with a median value of 0.12 mg/L. Concentrations of nitrate-N in 13 (7 percent) nutrient samples, 1992 to 2009, were greater than or equal to the USEPA MCL and Nebraska Title-118 standard of 10 mg/L, and concentrations of nitrate-N in 35 (18 percent) nutrient samples, 1992 to 2009, were greater than or equal to 5 mg/L, which is the PMRNRD action level for possible management implementation to reduce nitrate concentrations in groundwater. Of the 61 pesticides or pesticide degradates analyzed from 2007 to 2009, 21 were detected. Three of the 21 pesticides detected (alachlor, atrazine, and metolachlor) have established health-based criteria; all detections of these compounds were at concentrations less than their USEPA standards. From 2007 to 2009, 1 or more pesticide compounds were detected in 16 of the 82 network wells and in 18 of the 26 wells in well nests. From 2007 to 2009, the individual pesticide compounds that were detected most frequently were alachlor ethane sulfonic acid, a degradate of alachlor; deethylcyanazine acid, a degradate of cyanazine; and atrazine. Analytes with concentrations that exceeded 30 percent of the applicable Nebraska Title-118 standard were identified so that the PMRNRD can plan to monitor groundwater in the area and consider possible actions should the analyte concentrations continue to rise. The analytical results from the most recent samples collected in the network wells and all the wells in well nests from 1992 to 2009 indicate that, in at least 1 sample, there was a concentration that exceeded 30 percent of the Nebraska Title-118 standard for at least 1 of 3 major ions (chloride, fluoride, and sulfate), 1 nutrient (nitrate-N), 1 pesticide (atrazine), or 3 trace elements (arsenic, iron, and manganese). In addition, 30 percent of the USEPA MCL or Nebraska Title-118 standard for gross alpha activity likely was exceeded in samples from three wells screened in the Dakota aquifer. Study findings indicate that some alternatives to the current PMRNRD groundwater-sampling approach that could be considered are to collect fewer samples for nutrient analysis and to collect samples periodically for determining concentrations of additional analytes, particularly the analytes with concentrations that were at least 30 percent or more than the Nebraska Title-118 standard.

Nebraska↗
USGS2011

Selected approaches to estimate water-budget components of the High Plains, 1940 through 1949 and 2000 through 2009

The High Plains aquifer, underlying almost 112 million acres in the central United States, is one of the largest aquifers in the Nation. It is the primary water supply for drinking water, irrigation, animal production, and industry in the region. Expansion of irrigated agriculture throughout the past 60 years has helped make the High Plains one of the most productive agricultural regions in the Nation. Extensive withdrawals of groundwater for irrigation have caused water-level declines in many parts of the aquifer and increased concerns about the long-term sustainability of the aquifer. Quantification of water-budget components is a prerequisite for effective water-resources management. Components analyzed as part of this study were precipitation, evapotranspiration, recharge, surface runoff, groundwater discharge to streams, groundwater fluxes to and from adjacent geologic units, irrigation, and groundwater in storage. These components were assessed for 1940 through 1949 (representing conditions prior to substantial groundwater development and referred to as "pregroundwater development" throughout this report) and 2000 through 2009. Because no single method can perfectly quantify the magnitude of any part of a water budget at a regional scale, results from several methods and previously published work were compiled and compared for this study when feasible. Results varied among the several methods applied, as indicated by the range of average annual volumes given for each component listed in the following paragraphs. Precipitation was derived from three sources: the Parameter-Elevation Regressions on Independent Slopes Model, data developed using Next Generation Weather Radar and measured precipitation from weather stations by the Office of Hydrologic Development at the National Weather Service for the Sacramento-Soil Moisture Accounting model, and precipitation measured at weather stations and spatially distributed using an inverse-distance-weighted interpolation method. Precipitation estimates using these sources, as a 10-year average annual total volume for the High Plains, ranged from 192 to 199 million acre-feet (acre-ft) for 1940 through 1949 and from 185 to 199 million acre-ft for 2000 through 2009. Evapotranspiration was obtained from three sources: the National Weather Service Sacramento-Soil Moisture Accounting model, the Simplified-Surface-Energy-Balance model using remotely sensed data, and the Soil-Water-Balance model. Average annual total evapotranspiration estimated using these sources was 148 million acre-ft for 1940 through 1949 and ranged from 154 to 193 million acre-ft for 2000 through 2009. The maximum amount of shallow groundwater lost to evapotranspiration was approximated for areas where the water table was within 5 feet of land surface. The average annual total volume of evapotranspiration from shallow groundwater was 9.0 million acre-ft for 1940 through 1949 and ranged from 9.6 to 12.6 million acre-ft for 2000 through 2009. Recharge was estimated using two soil-water-balance models as well as previously published studies for various locations across the High Plains region. Average annual total recharge ranged from 8.3 to 13.2 million acre-ft for 1940 through 1949 and from 15.9 to 35.0 million acre-ft for 2000 through 2009. Surface runoff and groundwater discharge to streams were determined using discharge records from streamflow-gaging stations near the edges of the High Plains and the Base-Flow Index program. For 1940 through 1949, the average annual net surface runoff leaving the High Plains was 1.9 million acre-ft, and the net loss from the High Plains aquifer by groundwater discharge to streams was 3.1 million acre-ft. For 2000 through 2009, the average annual net surface runoff leaving the High Plains region was 1.3 million acre-ft and the net loss by groundwater discharge to streams was 3.9 million acre-ft. For 2000 through 2009, the average annual total estimated groundwater pumpage volume from two soil-water-balance models ranged from 8.7 to 16.2 million acre-ft. Average annual irrigation application rates for the High Plains ranged from 8.4 to 16.2 inches per year. The USGS Water-Use Program published estimated total annual pumpage from the High Plains aquifer for 2000 and 2005. Those volumes were greater than those estimated from the two soil-water-balance models. Total groundwater in storage in the High Plains aquifer was estimated as 3,173 million acre-ft prior to groundwater development and 2,907 million acre-ft in 2007. The average annual decrease of groundwater in storage between 2000 and 2007 was 10 million acre-ft per year.

High Plains Aquifer↗