USGS Science⌕ Search

USGS · sir20185136

Simulation of groundwater flow and analysis of projected water use for the Rush Springs aquifer, western Oklahoma

Abstract

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation and the Oklahoma Water Resources Board, (1) quantified the groundwater resources of the Rush Springs aquifer in western Oklahoma by developing a numerical groundwater-flow model, (2) evaluated the effects of estimated equal-proportionate-share (EPS) pumping rates on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future, (3) assessed the uncertainty in the EPS scenario results, and (4) evaluated the effects of (a) projected groundwater-use rates extended 50 years into the future and (b) sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage. The Rush Springs aquifer is an important source of water for municipal and irrigation use by many communities and agricultural users in the study area. The study area is composed of about 4,970 square miles (3,181,003 acres) of Rush Springs aquifer bedrock deposits located in 14 counties. The study area also includes the alluvium and terrace deposits of the Canadian and Washita Rivers, as well as alluvium along the Little Washita River, Deer Creek, and a number of smaller tributaries of the Washita River that overlie the bedrock. A numerical groundwater-flow model of the Rush Springs aquifer was constructed by using MODFLOW with the Newton solver. Groundwater flow was simulated for January 1979–December 2015 by using monthly stress periods, and an initial steady-state stress period was configured to represent mean annual inflows and outflows. The model was calibrated to groundwater-level observations at selected wells, monthly base flow at nine streamgages, stream seepage as estimated for the conceptual water budget, and Fort Cobb Reservoir stage. The EPS scenarios for the Rush Springs 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.82, 0.49, and 0.43 acre-foot per acre per year, respectively. Given the 2,954,545-acre aquifer area used for the EPS scenarios, the 20-year rate corresponds to an annual yield of about 2,422,727 acre-feet per year. Groundwater storage at the end of the 20-year EPS scenario was about 13,321,000 acre-feet, or about 31,516,437 acre-feet (70 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean groundwater-level decline of about 152 feet. Water availability under the EPS pumping rate was primarily from the western area of the model. Saturation was sustained though the entire EPS scenario where the aquifer was sufficiently thick or a shallow hydraulic gradient was present. Fort Cobb Reservoir stage was below the dead-pool stage after about 5 years of 20-year EPS pumping. An uncertainty analysis was conducted to assess the uncertainty in the EPS scenario results. An ensemble of 400 random sets of possible parameter values was performed for the uncertainty analysis by using a multivariate normal distribution centered on the calibrated parameter values. The parameter bounds for the uncertainty analysis were determined by using the posterior covariance matrix, which allows for the incorporation of knowledge gained during the calibration process as well as observation uncertainty and the correlation between estimated parameters. The uncertainty results indicate a 95-percent confidence interval for the 20-year EPS pumping rate between 0.73 and 0.95 acre-foot per acre per year. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Rush Springs aquifer. The effects of well withdrawals were evaluated by comparing changes in groundwater storage between four 50-year scenarios using (1) no groundwater use, (2) mean groundwater use for the study period (1979–2015), (3) increasing groundwater use, and (4) groundwater use at the 2015 rate. The increasing-use scenario assumed a 38-percent increase in pumping over 50 years on the basis of 2010–60 demand projections for western Oklahoma. Simulated groundwater storage changes ranged between an increase of 6.3 percent for the scenario with no groundwater use, and 0.9 percent for the scenario with 2015 groundwater-use rates. For the Fort Cobb Reservoir surface watershed, simulated groundwater storage changes ranged between an increase of 23.6 percent for the scenario with no groundwater use and a decrease of 4.0 percent for the increasing groundwater-use scenario. Groundwater-level changes were generally greater in areas with a large concentration of groundwater wells and groundwater use such as the Fort Cobb Reservoir surface watershed. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on the Rush Springs aquifer groundwater storage and Fort Cobb Reservoir stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model. To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–1992), and upstream inflows to the Canadian and Washita Rivers and associated tributaries were reduced by 37 percent. Groundwater storage at the end of the hypothetical drought period in December 1992 was about 42,983,000 acre-feet, or about 3,525,000 acre-feet (7.6 percent) less than the groundwater storage of the calibrated numerical model. This change in groundwater storage is equivalent to a mean groundwater-level decline of 15.8 feet. Simulated mean base-flow declines at the Canadian and Washita River streamgages were between 39 and 59 percent during the drought period. The minimum stage in Fort Cobb Reservoir at the end of the hypothetical drought period was 1,311 feet, indicating a storage capacity of only 10 percent of active conservation pool storage. The Fort Cobb Reservoir storage declines mostly resulted from reduced base flows in Cobb, Lake, and Willow Creeks upstream from the reservoir.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 34.5° to 36.5° latitude; -99.75° to -97.75° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J.H. Ellis. 2018-11-29. Simulation of groundwater flow and analysis of projected water use for the Rush Springs aquifer, western Oklahoma. https://doi.org/10.3133/sir20185136

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Baseline water-quality inventory and summary for the Paterson Great Falls National Historical Park, New Jersey

The Paterson Great Falls National Historical Park (PAGR), established in 2011 in the City of Paterson, New Jersey, contains many cultural, historical, and natural resources. Because of the PAGR’s young age and limited resources, a long record of natural resource data is not available for use to guide future park management and planning decisions. A detailed natural resource inventory was identified as an immediate need in the “Foundation Document—Paterson Great Falls National Historical Park,” given the history of extensive industrial activity throughout the PAGR and the Passaic River Basin. The U.S. Geological Survey, in cooperation with the National Park Service and the Urban Waters Federal Partnership for the Lower Passaic River, inventoried water-quality monitoring efforts and compiled historical water-quality data from secondary data sources to define baseline conditions in the PAGR and the surrounding study area. This report documents the methods used to harmonize water-quality data aggregated from multiple sources for the PAGR study area. Procedures implemented during data harmonization include compiling data from multiple sources, assessing the dataset for completeness, addressing metadata concerns, and creating a consistent and quality-assured dataset for future analysis. The harmonized dataset contains 14 aggregated surface-water sites, 1,103 individual sampling days, and 32,000 water-quality records throughout the period of record defined for this study (1952–2021). Summary statistics and comparison to water-quality criteria are provided for focus parameters that can support PAGR management, including field properties, suspended material, total dissolved solids, specific conductance, chloride, sulfate, Escherichia coli, nutrients, contaminants of concern related to fish consumption advisories, and litter and floatable debris. A spatial data gap was identified for the recreational area near the base of the Great Falls and immediately downstream from a known combined sewer overflow discharge location. Analytic data gaps were identified for continuous streamflow and water-quality parameters, floatable debris characterization, and harmful algal blooms identification. Continuous collection of contaminants of concern was identified as a temporal data gap, as these parameters were not consistently collected within the study area. The water-quality data synthesis and summary presented in this report can assist PAGR staff in understanding the complex baseline water resources within the study area and can serve as a baseline to design future water-quality monitoring efforts and evaluate changes in water quality over time. The harmonized dataset, as well as future New Jersey Department of Environmental Protection and New Jersey Harbor Discharge Group monitoring data and assessments, can also be used to examine the status and trends in surface-water quality, support compliance with the Clean Water Act (33 U.S.C. 1251 et seq.), and further investigate changes to park water resources that can affect ecological and human health.

New Jersey↗

Slip tendency and dilation tendency of Quaternary faults in the Great Basin region, Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming

This report presents the results of slip tendency and dilation tendency analysis of Quaternary faults in the Great Basin region, in the states of Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming. This effort was undertaken as part of two studies: the U.S. Department of Energy Geothermal Technologies Office-funded Innovative Geothermal Exploration through Novel Investigations of Undiscovered Systems (INGENIOUS) project, which is a multidisciplinary, multi-institution effort to develop new methodologies and best practices to accelerate the discovery of new, commercially viable hydrothermal resources; and the U.S. Geological Survey Geothermal Resources Investigations Project. The goal of this work is to identify faults and fault segments that are appropriately oriented to be stress-loaded for slip or to dilate under the ambient stress conditions. Both conditions may make such faults more likely to host permeability and fluid flow relative to other faults in the region. Slip tendency and dilation tendency were calculated for all fault segments in the study area from the U.S. Geological Survey Quaternary Fault and Fold Database. Results demonstrate that faults with high slip tendency and dilation tendency are common in the Great Basin region. The areas containing many high slip tendency and dilation tendency fault segments are associated with high strain rates. These areas also host a relatively high proportion of the known hydrothermal systems. This indicates that faults with high slip tendency and dilation tendency in areas with high strain-rate are potentially most prospective for hosting hydrothermal systems. The database published in coordination with this report contains the slip tendency and dilation tendency results data, including the fault geometry and stress data used to calculate these quantities, for both the Great Basin study area and the smaller INGENIOUS study area.

Arizona, California, Idaho, Nevada, Oregon, Utah, ↗

A spatial analysis of the groundwater emergence flood hazard in Long Island, New York, and coastal areas surrounding Long Island Sound in New York, Connecticut, and Rhode Island

Long Island, New York, and nearby coastal areas surrounding Long Island Sound are densely populated and, like other coastal areas, are susceptible to flooding from several potential sources, including stormwater from precipitation events, tidal flooding and storm surge, and groundwater emergence flooding. Groundwater emergence fooding refers to the intersection of a rising water table with land surface or critical infrastructure. Many studies of flood drivers either neglect or only briefly discuss how shallow groundwater conditions may contribute to or exacerbate flood conditions. As a part of a comprehensive study of compound flood hazards in the coastal areas of Long Island and Long Island Sound, a spatial analysis was completed by the U.S. Geological Survey, in cooperation with the Environmental Protection Agency’s Long Island Sound Partnership, formerly called the Long Island Sound Study, using available regional datasets to characterize the potential hazard for groundwater emergence flooding. The approximately 3,100 square mile study area was subdivided into 11,407 grid cells that were 900 meters by 900 meters (approximately 3,000 feet by 3,000 feet), for the purpose of integrating the spatial datasets to calculate and map the groundwater emergence flood hazard. The depth to the water table, hydrologic soil groups, and National Land Cover Database were harmonized to the common grid. A groundwater emergence flood hazard score was calculated for each grid cell for current average conditions following a set of rules accounting for the depth to the water table and the percentage of area within each cell with slow infiltrating soils. A scenario with a higher sea-level position was also calculated for the Long Island part of the study area. The calculated groundwater emergence flood hazard score was reviewed in conjunction with the National Land Cover Database to identify developed areas and associated infrastructure that may be at risk of potential groundwater emergence flooding. Study results indicate that the groundwater emergence flood hazard is highest in coastal areas and near surface water where the water table is close to ground surface. Inland areas away from surface waterbodies are not likely to be exposed to groundwater emergence flooding. For Long Island, in the scenario with higher sea-level position, a greater groundwater emergence flood hazard is calculated for some locations closer to the coast and where land is submerged. Away from the coast and surface-water drainage, the groundwater emergence flood hazard is similar between the current average sea-level condition and the higher sea-level position scenario.

Connecticut, New York, Rhode Island↗