Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009

Kyle W. Davis; Larry D. Putnam

doi:10.3133/sir20135069

Home/Archive/USGS

USGS · sir20135069

Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009

Kyle W. Davis·Larry D. Putnam

View original source ↗Read source document ↗Export metadata ↓

Abstract

The Ogallala aquifer is an important water resource for the Rosebud Sioux Tribe in Gregory and Tripp Counties in south-central South Dakota and is used for irrigation, public supply, domestic, and stock water supplies. To better understand groundwater flow in the Ogallala aquifer, conceptual and numerical models of groundwater flow were developed for the aquifer. A conceptual model of the Ogallala aquifer was used to analyze groundwater flow and develop a numerical model to simulate groundwater flow in the aquifer. The MODFLOW–NWT model was used to simulate transient groundwater conditions for water years 1985–2009. The model was calibrated using statistical parameter estimation techniques. Potential future scenarios were simulated using the input parameters from the calibrated model for simulations of potential future drought and future increased pumping. Transient simulations were completed with the numerical model. A 200-year transient initialization period was used to establish starting conditions for the subsequent 25-year simulation of water years 1985–2009. The 25-year simulation was discretized into three seasonal stress periods per year and used to simulate transient conditions. A single-layer model was used to simulate flow and mass balance in the Ogallala aquifer with a grid of 133 rows and 282 columns and a uniform spacing of 500 meters (1,640 feet). Regional inflow and outflow were simulated along the western and southern boundaries using specified-head cells. All other boundaries were simulated using no-flow cells. Recharge to the aquifer occurs through precipitation on the outcrop area. Model calibration was accomplished using the Parameter Estimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. This program was designed to estimate parameter values that are statistically the most likely set of values to result in the smallest differences between simulated and observed values, within a given set of constraints. The potentiometric surface of the aquifer calculated during the 200-year initialization period established initial conditions for the transient simulation. Water levels for 38 observation wells were used to calibrate the 25-year simulation. Simulated hydraulic heads for the transient simulation were within plus or minus 20 feet of observed values for 95 percent of observation wells, and the mean absolute difference was 5.1 feet. Calibrated hydraulic conductivity ranged from 0.9 to 227 feet per day (ft/d). The annual recharge rates for the transient simulation (water years 1985–2009) ranged from 0.60 to 6.96 inches, with a mean of 3.68 inches for the Ogallala aquifer. This represents a mean recharge rate of 280.5 ft 3 /s for the model area. Discharge from the aquifer occurs through evapotranspiration, discharge to streams through river leakage and flow from springs and seeps, and well withdrawals. Water is withdrawn from wells for irrigation, public supply, domestic, and stock uses. Simulated mean discharge rates for water years 1985–2009 were about 185 cubic feet per second (ft 3 /s) for evapotranspiration, 66.7 ft 3 /s for discharge to streams, and 5.48 ft 3 /s for well withdrawals. Simulated annual evapotranspiration rates ranged from about 128 to 254 ft 3 /s, and outflow to streams ranged from 52.2 to 79.9 ft 3 /s. A sensitivity analysis was used to examine the response of the calibrated model to changes in model parameters for horizontal hydraulic conductivity, recharge, evapotranspiration, and spring and riverbed conductance. The model was most sensitive to recharge and maximum potential evapotranspiration and least sensitive to riverbed and spring conductances. Two potential future scenarios were simulated: a potential drought scenario and a potential increased pumping scenario. To simulate a potential drought scenario, a synthetic drought record was created, the mean of which was equal to 60 percent of the mean estimated recharge rate for the 25-year simulation period. Compared with the results of the calibrated model (non-drought simulation), the simulation representing a potential drought scenario resulted in water-level decreases of as much as 30 feet for the Ogallala aquifer. To simulate the effects of potential future increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 25-year simulation period. Compared with the results of the calibrated model, the simulation representing an increased pumping scenario resulted in water-level decreases of as much as 26 feet for the Ogallala aquifer. Groundwater budgets for the potential future scenario simulations were compared with the transient simulation representing water years 1985–2009. The simulation representing a potential drought scenario resulted in lower aquifer recharge from precipitation and decreased discharge from streams, springs, seeps, and evapotranspiration. The simulation representing a potential increased pumping scenario was similar to results from the transient simulation, with a slight increase in well withdrawals and a slight decrease in discharge from river leakage and evapotranspiration. This numerical model is suitable as a tool that could be used to better understand the flow system of the Ogallala aquifer, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the study area. The model also is useful to help assess the response of the aquifer to additional stresses, including potential drought conditions and increased well withdrawals.

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 42.9979° to 43.7619° latitude; -100.2333° to -98.4985° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Kyle W. Davis, Larry D. Putnam. 2013. Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009. https://doi.org/10.3133/sir20135069

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

KEEP EXPLORING

Related USGS reports

USGS2026

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↗

USGS2026

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

USGS2026

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↗