Skip Links

USGS - science for a changing world

Scientific Investigations Report 2011–5155

Prepared in cooperation with the Bureau of Reclamation, Washington State Department of Ecology, and the Yakama Nation

Numerical Simulation of Groundwater Flow for the Yakima River Basin Aquifer System, Washington

By D.M. Ely, M.P. Bachmann, and J.J. Vaccaro

Thumbnail of and link to report PDF (15.6 MB)Abstract

A regional, three-dimensional, transient numerical model of groundwater flow was constructed for the Yakima River basin aquifer system to better understand the groundwater-flow system and its relation to surface-water resources. The model described in this report can be used as a tool by water-management agencies and other stakeholders to quantitatively evaluate proposed alternative management strategies that consider the interrelation between groundwater availability and surface-water resources.

The model was constructed using the U.S. Geological Survey finite-difference model MODFLOW. The model uses 1,000-foot grid cells that subdivide the model domain by 600 rows and 600 columns. Forty-eight hydrogeologic units in the model are included in 24 model layers. The Yakima River, all major tributaries, and major agricultural drains are included in the model as either drain cells or streamflow-routing cells. Recharge was estimated from previous work using physical process models. Groundwater pumpage specified in the model is derived from monthly pumpage values previously estimated from another component of this study. The pumpage values include estimates for wells with standby/reserve rights that are used in drought years.

The model was calibrated to the transient conditions for October 1959 to September 2001. Calibration was completed by using traditional trial-and-error methods and automated parameter-estimation techniques. The model simulates the shape and slope of the water table that generally is consistent with mapped water levels. At well observation points, the average difference between simulated and measured hydraulic heads is -49 feet with a root-mean-square error divided by the total difference in water levels of 4 percent. Simulated river streamflow was compared to measured streamflow at seven sites. Annual differences between measured and simulated streamflow for the sites ranged from 1 to 9 percent. Calibrated model output includes a 42-year estimate of a monthly water budget for the aquifer system.

Five applications (scenarios) of the model were completed to obtain a better understanding of the relation between pumpage and surface-water resources and groundwater levels. For the first three scenarios, the calibrated transient model was used to simulate conditions without: (1) pumpage from all hydrogeologic units, (2) pumpage from basalt hydrogeologic units, and (3) exempt-well pumpage. The simulation results indicated potential streamflow capture by the existing pumpage from 1960 through 2001. The quantity of streamflow capture generally was inversely related to the total quantity of pumpage eliminated in the model scenarios. For the fourth scenario, the model simulated 1994 through 2001 under existing conditions with additional pumpage estimated for pending groundwater applications. The differences between the calibrated model streamflow and this scenario indicated additional decreases in streamflow of 91 cubic feet per second in the model domain. Existing conditions representing 1994 through 2001 were projected through 2025 for the fifth scenario and indicated additional streamflow decreases of 38 cubic feet per second and groundwater-level declines.

First posted September 16, 2011

For additional information contact:
Director, Washington Water Science Center
U.S. Geological Survey
934 Broadway, Suite 300
Tacoma, Washington 98402

Part or all of this report is presented in Portable Document Format (PDF); the latest version of Adobe Reader or similar software is required to view it. Download the latest version of Adobe Reader, free of charge.

Suggested citation:

Ely, D.M., Bachmann, M.P., and Vaccaro, J.J., 2011, Numerical simulation of groundwater flow for the Yakima River basin aquifer system, Washington: U.S. Geological Survey Scientific Investigations Report 2011-5155, 90 p.




Purpose and Scope

Description of Study Area

Well Data and Their Use

Hydrogeologic Units

Numerical Simulation of Groundwater Flow

Model Calibration and Sensitivity

Model Uncertainty and Limitations

Appropriate Uses of the Groundwater Model

Model-Derived Groundwater Budget

Model Applications

Summary and Conclusions


References Cited

Appendix A. Assignment of Groundwater Pumpage to Model Layers

Appendix B. Estimates of Pumpage for Pending Groundwater Applications

Accessibility FOIA Privacy Policies and Notices logo U.S. Department of the Interior | U.S. Geological Survey
Page Contact Information: Contact USGS
Page Last Modified: Thursday, January 10, 2013, 08:10:20 PM