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

Analysis of groundwater geochemistry in the Lansing, Michigan, area

Groundwater resources in the Lansing, Michigan, area are protected by a wellhead protection plan, but previous studies have indicated that groundwater in this area may be susceptible to surface contamination. The U.S. Geological Survey, in cooperation with the Lansing Board of Water & Light, analyzed geochemical and aquifer characteristics to assess the potential for contamination of groundwater resources. Recharge age categories were determined from tritium concentrations for some wells, and by using ternary diagrams and principal component analysis for other wells. These age classifications, combined with hydrogeologic and geochemical data, were analyzed by using statistical methods to identify trends that may indicate areas of increased groundwater susceptibility. Wells with modern recharge tended to have higher sodium, chloride, and chloride:bromide ratios than did wells with mixed recharge, suggesting greater surface connectivity and potential influence from anthropogenic sources such as road salt. In contrast, wells in the mixed age category exhibited more stable geochemistry and lower constituent concentrations. Classification and regression tree models identified aquifer geometry—particularly the thickness of the uppermost bedrock unit and presence of confining layers—as key predictors of recharge age and geochemical indicators. The trends identified from these analyses could help to reveal locations in the Lansing area that are among the most susceptible to contamination.

Scientific Investigations Report↗

USGS2026

Evaluating the number of realizations (train-test splits) necessary to ensure robust estimates of conventional hydrothermal resource favorability using Monte Carlo cross-validation, the Great Basin, Western USA

Supervised machine learning approaches used to produce favorability maps in resource assessments commonly partition an input dataset into training and testing subsets (that is, through train-test splits). However, if a single split is not representative of underlying resource conditions, the resulting favorability map may be biased, particularly with sparse or heterogeneous datasets. To mitigate the potential for unrepresentative favorability maps, recent works by U.S. Geological Survey researchers have used Monte Carlo cross-validation to implement a large number ( n) of randomized train-test splits to produce n realizations. The realizations are then aggregated using cell-wise median predictions to produce the resulting final favorability map, recognizing that the median will remove extreme values that might have been estimated from a comparatively few unrepresentative realizations. Theoretically, the median prediction at every cell converges to a single estimate as n approaches infinity, and the question addressed herein is: How many realizations are necessary to get a robust estimate (that is, an estimate unaffected by additional realizations)? Therefore, convergence metrics are described to identify the number of train-test splits and realizations needed to make a reliable map using the data and supervised machine learning workflow for conventional hydrothermal resource assessment of the Great Basin of the Western United States. For this study, robust prediction occurs with n approximately equal to 600 realizations, and smaller ensembles (for example, n is approximately equal to 200 or even 100) may be suitable during preliminary or exploratory analyses. Although the resulting threshold numbers will vary depending on the data and modeling approach used, the convergence metrics proposed and applied herein are applicable to other supervised machine learning studies that use Monte Carlo cross-validation.

California, Idaho, Nevada, Oregon, Utah↗

USGS2026

Groundwater budget for the Portneuf River Basin, southeastern Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, developed the first groundwater budget for the upper Portneuf River Basin in southeastern Idaho, a tributary of the Snake River. A groundwater budget for the Portneuf River Basin supports water resource and administrative management decisions within the basin and for downstream users of the eastern Snake River Plain aquifer. The analysis quantifies inflows, outflows, and residuals for an aquifer system defined by the Idaho Department of Water Resources under three recent precipitation conditions: an average year (2016), a dry year (2021), and a wet year (2023). The budget was calculated for the basin and for two subbasins delineated by the U.S. Geological Survey streamgage Portneuf River at Topaz, Idaho (13073000): the above Topaz subbasin, which includes the Portneuf-Gem valley, and the below Topaz subbasin, which includes Marsh Valley. Aquifer inflows include mountain-front recharge; areal recharge; recharge from applied water, including recharge from irrigation, canal seepage, and effluent returns; groundwater flow into the basin; and reservoir seepage. Outflows include groundwater pumping for agricultural, municipal, and domestic uses; subsurface outflow through the Portneuf gap; and discharge to streams as base flow. Estimated annual inflows for the full basin ranged from 116,880 acre-feet in 2021 to 154,880 acre-feet in 2023, and outflows ranged from 132,190 acre-feet in 2016 to 158,490 acre-feet in 2023. Basinwide residuals, the change in storage plus uncertainty, were −2,410 acre-feet in 2016, −18,410 acre-feet in 2021, and −3,610 acre-feet in 2023, and each subbasin had a mixture of negative and positive residuals. These basinwide results, interpreted as change in storage plus water-budget uncertainty, indicate minor net losses that increase with dry conditions, consistent with well depth-to-water measurements indicating gradual multi-decadal declines. Patterns in subbasin residuals are driven by irrigation practices and surface-water management; for example, the above Topaz subbasin relies more on groundwater for irrigation, resulting in more negative residuals during dry years. In contrast, the below Topaz subbasin depends more on surface water, producing positive residuals during dry years because of irrigation losses and canal seepage. The analyses for both subbasins indicate decades-long groundwater declines, which agree with the overall negative residuals calculated. The largest sources of uncertainty include climate inputs (precipitation and evapotranspiration) and poorly constrained components such as groundwater inflow and tributary stream infiltration. These findings provide insight into how climatic variables and irrigation practices affect groundwater inflows and outflows within the Portneuf River Basin.

Idaho↗