USGS Visual Identifier

Simulated Ground-Water Flow in the Ogallala and Arikaree Aquifers, Rosebud Indian Reservation Area, South Dakota

By Andrew J. Long, Larry D. Putnam, and Janet M. Carter

Water-Resources Investigations Report 03-4043

Prepared in cooperation with the Rosebud Sioux Tribe


The Ogallala and Arikaree aquifers are important water resources in the Rosebud Indian Reservation area and are used extensively for irrigation, municipal, and domestic water supplies. Continued or increased withdrawals from the Ogallala and Arikaree aquifers in the Rosebud Indian Reservation area have the potential to affect water levels in these aquifers. This report describes a conceptual model of ground-water flow in these aquifers and documents the development and calibration of a numerical model to simulate ground-water flow. Data for a twenty-year period (water years 1979 through 1998) were analyzed for the conceptual model and included in steady-state and transient numerical simulations of ground-water flow for the same 20-year period.

A three-dimensional ground-water flow model, with two layers, was used to simulate ground-water flow in the Ogallala and Arikaree aquifers. The upper layer represented the Ogallala aquifer, and the lower layer represented the Arikaree aquifer. The study area was divided into grid blocks 1,640 feet (500 meters) on a side, with 153 rows and 180 columns.

Areal recharge to the Ogallala and Arikaree aquifers occurs from precipitation on the outcrop areas. The recharge rate for the steady-state simulation was 3.3 inches per year for the Ogallala aquifer and 1.7 inches per year for the Arikaree aquifer for a total recharge rate of 266 cubic feet per second.

Discharge from the Ogallala and Arikaree aquifers occurs through evapotranspiration, discharge to streams, and well withdrawals. Discharge rates in cubic feet per second for the steady-state simulation were 184 for evapotranspiration, 46.8 and 19.7 for base flow to the Little White and Keya Paha Rivers, respectively, and 11.6 for well withdrawals from irrigation use. Estimated horizontal hydraulic conductivity used for the numerical model ranged from 0.2 to 120 feet per day in the Ogallala aquifer and 0.1 to 5.4 feet per day in the Arikaree aquifer. A uniform vertical hydraulic conductivity value of 6.6x10-4 feet per day was applied to the Ogallala aquifer. Vertical hydraulic conductivity was estimated for five zones in the Arikaree aquifer and ranged from 8.6x10-6 to 7.2x10-1 feet per day. Average rates of recharge, maximum evapotranspiration, and well withdrawals were included in the steady-state simulation, whereas the time-varying rates were included in the transient simulation.

Model calibration was accomplished by varying parameters within plausible ranges to produce the best fit between simulated and observed hydraulic heads and base-flow discharges from the Ogallala and Arikaree aquifers. For the steady-state simulation, the root mean square error for simulated hydraulic heads for all wells was 26.8 feet. Simulated hydraulic heads were within 50 feet of observed values for 95 percent of the wells. For the transient simulation, the difference between the simulated and observed means for hydrographs was within 40 feet for all observation wells. The potentiometric surfaces of the two aquifers calculated by the steady-state simulation established initial conditions for the transient simulation.

A sensitivity analysis was used to examine the response of the calibrated steady-state model to changes in model parameters including horizontal and vertical hydraulic conductivity, evapotranspiration, recharge, and riverbed conductance. The model was most sensitive to recharge and horizontal hydraulic conductivity.





Purpose and Scope


Description of Study Area

Physiography, Land Use, and Climate

Drainage Features and Streamflow


Hydrogeologic Setting

Simulated Ground-Water Flow

Conceptual Model

Potentiometric Surface

Hydraulic Properties




Discharge to Streams

Well Withdrawals

Numerical Model

Model Design

Grid and Boundary Conditions

Representation of Hydraulic Properties

Representation of Recharge and Discharge

Model Calibration

Steady-State Simulation

Transient Simulation

Model Limitations



Supplemental Information

This report is available online in Portable Document Format (PDF). If you do not have the Adobe Acrobat PDF Reader, it is available for free download from Adobe Systems Incorporated.

Download the full report (PDF, 6.1MB).

The report can be downloaded in parts for a faster download time:

Pages 0-35 (3.3MB)

Pages 36-69 (3.0MB)

Document Accessibility: Adobe Systems Incorporated has information about PDFs and the visually impaired. This information provides tools to help make PDF files accessible. These tools convert Adobe PDF documents into HTML or ASCII text, which then can be read by a number of common screen-reading programs that synthesize text as audible speech. In addition, an accessible version of Acrobat Reader 5.0 for Windows (English only), which contains support for screen readers, is available. These tools and the accessible reader may be obtained free from Adobe at Adobe Access.

Send questions or comments about this report to the author, Andrew J. Long (605) 355-4560 ext. 237.

For more information about USGS activities in South Dakota, visit the USGS South Dakota District home page.

U.S. Department of the Interior, U.S. Geological Survey
Persistent URL:
Page Contact Information: Contact USGS
Last modified: Wednesday, December 07 2016, 01:16:23 PM
FirstGov button  Take Pride in America button