Montana Bureau of Mines and Geology Open-File Report 638 2013
HYDROGEOLOGIC INVESTIGATION OF THE LOWER BEAVERHEAD STUDY AREA BEAVERHEAD COUNTY, MONTANA GROUNDWATER MODELING REPORT
Julie A. Butler and Ginette Abdo Montana Bureau of Mines and Geology Ground Water Investigations Program
HYDROGEOLOGIC INVESTIGATION OF THE LOWER BEAVERHEAD STUDY AREA BEAVERHEAD COUNTY, MONTANA
GROUNDWATER MODELING REPORT
Julie A. Butler and Ginette Abdo
Montana Bureau of Mines and Geology Ground Water Investigations Program
Montana Bureau of Mines and Geology Open-File Report 638
2013
Montana Bureau of Mines and Geology Open-File Report 638 TABLE OF CONTENTS
Introduction ...... 1 Report Purpose ...... 1 General Setting ...... 1 Climate ...... 1 Physiography ...... 3 Man-Made Hydrologic Features ...... 4 Model Objectives ...... 4 Conceptual Model ...... 4 Geologic Framework ...... 5 Hydrogeologic Units ...... 7 Alluvial Aquifer ...... 8 Tertiary Sediment Aquifer...... 10 Volcanic Rock Aquifer ...... 12 Groundwater Flow System ...... 13 Vertical Gradients ...... 16 Hydrologic Boundaries ...... 18 Aquifer Properties ...... 18 Tertiary Sediment Aquifer Test ...... 18 Volcanic Rock Aquifer Test ...... 19 Alluvial Aquifer Tests ...... 21 Summary of Hydraulic Properties ...... 21 Sources and Sinks ...... 21 Groundwater Budget ...... 21 Computer Code ...... 24 Groundwater Flow Model Construction ...... 25 Model Grid ...... 25 Hydraulic Parameters ...... 27 Boundary Conditions ...... 27 Model Borders ...... 27 Surface-Water Bodies ...... 27 Irrigation Field Recharge ...... 28 Well Withdrawals ...... 28 Calibration ...... 29 Selection of Calibration Targets ...... 29 Steady-State Calibration ...... 30 Transient Calibration ...... 33 Sensitivity Analysis ...... 36 Model Veri ication ...... 36 Predictive Simulations ...... 37 Scenario 1 ...... 38 Scenario 2 ...... 38 Scenario 3 ...... 39 Scenario 4 ...... 41 Canal Seepage Scenarios ...... 41 Scenarios 5 and 6 ...... 41 Scenario 7 ...... 47 Summary and Conclusions ...... 48
v Butler and Abdo, 2013 Assumptions and Limitations ...... 48 Model Predictions ...... 50 Recommendations ...... 52 References ...... 52 Appendix A. Model File Index ...... 55 Appendix B. Estimating Irrigation Recharge in the Groundwater Budget ...... 61 Appendix C. Estimating Well Withdrawals in the Groundwater Budget ...... 65 Appendix D. NRCS Irrigation Water Requirements—Crop Data Summaries ...... 69
FIGURES Figure 1. Site location and physiography ...... 2 Figure 2. Study area ...... 3 Figure 3. Deviation from annual average precipitation in Dillon, 1900–2010 ...... 4 Figure 4. Beaverhead River and tributaries of interest in the model ...... 5 Figure 5. Study area geology ...... 6 Figure 6. Cross section through the lower Beaverhead River basin ...... 8 Figure 7. Exposures of the primary hydrogeologic units in the study area ...... 9 Figure 8. Location of sedimentary bedrock outcrop, observed August 4, 2010 ...... 10 Figure 9. Groundwater monitoring sites, including MBMG well installations ...... 11 Figure 10. 3D cross sections through loodplain ...... 11 Figure 11. Estimation of the volcanic rock aquifer extent ...... 12 Figure 12. Groundwater low system—bench potentiometric surface ...... 13 Figure 13. Groundwater low system— loodplain (shallow alluvial) potentiometric surface ....14 Figure 14. Groundwater levels in luenced by irrigation practices ...... 15 Figure 15. Vertical gradient monitoring sites ...... 16 Figure 16. Site A and B hydrographs ...... Figure 17. Site F hydrographs ...... 18 Figure 18. Tertiary sediment aquifer test site ...... 19 Figure 19. Volcanic rock aquifer test site ...... 20 Figure 20. Active cell coverage in the Lower Beaverhead model ...... 25 Figure 21. Model cross section through Row 237 of the model ...... 26 Figure 22. Model boundary conditions ...... 28 Figure 23. Irrigation ield recharge in the model ...... 29 Figure 24. Hydraulic conductivity array generated by steady-state calibration of the model .....32 Figure 25. Potentiometric surface generated by steady-state calibration of the model ...... 33 Figure 26. Computed vs. observed hydraulic head in the steady-state model ...... 34 Figure 27. Computed vs. observed heads of selected wells from the 2010 transient model ...... 35 Figure 28. RMS error caused by changing K and recharge values ...... 37 Figure 29. Potentiometric surface and head residuals from a 50% decrease in K values ...... 38 Figure 30. Potentiometric surface and head residuals from a 50% increase in recharge values ...... 39 Figure 31. Well locations in pumping and canal-seepage scenarios ...... 41 Figure 32. Stream depletion in the inal year of scenario 2 ...... 42 Figure 33. Spatial distribution of stream depletion in Beaverhead River during year 20 of scenario 2 ...... 43 Figure 34. Stream depletion in the inal year of scenario 3 ...... 44 Figure 35. Stream depletion in Beaverhead River throughout scenario 3...... 45 Figure 36. Drawdown in wells OW-1 and OW-2 in scenario 3 ...... 45 Figure 37. Stream depletion in the inal year of scenario 4 ...... 46 vi Montana Bureau of Mines and Geology Open-File Report 638 Figure 38. Base low to Willard and Black Sloughs in scenarios 1, 3, and 5 ...... 47 Figure 39. Base low to Willard and Black Sloughs in scenarios 1, 3, and 6 ...... 48 Figure 40. Base low to Willard and Black Sloughs in scenarios 1, 3, and 7 ...... 49 Figure 41. Base low to Beaverhead River in scenarios 1, 3, and 7 ...... 50 Figure 42. Plan view of maximum drawdown in scenarios 3 and 7 ...... 51
TABLES Table 1. Ranges of aquifer property values based on aquifer test results in the study area ...... 21 Table 2. Hydraulic conductivity values used in other western Montana large-area groundwater models ...... 22 Table 3. Comparison of ield-estimated and model aquifer property values ...... 22 Table 4. Details of the model grid as listed in GMS ...... 26 Table 5. Base low comparison between ield estimates and model results ...... 34 Table 6. Stress periods and time steps in the 2010 transient model run ...... 35 Table 7. Water budget values in the steady-state and transient versions of the model ...... 36 Table 8. Summary of pumping and canal seepage scenarios ...... 40
vii
Montana Bureau of Mines and Geology Open-File Report 638 INTRODUCTION in a previous study (MBMG, 2008) and because several West Bench well-permit applications were The main economy in the lower Beaverhead awaiting approval for some time. The study in- River basin is agriculture, which is dependent on volved aquifer testing, groundwater and surface- groundwater and surface-water irrigation. The water monitoring, well installations, water budget basin was closed to new surface-water appropria- analyses, and an evaluation of groundwater-level tions by Legislative authority effective April 1, trends. Results of the investigation were intended 1993, as part of the Jefferson–Madison River Basin to provide a scienti ic basis in which to help land- closure [Montana Department of Natural Resources owners and governing agencies make informed (DNRC), 2003]. In a closed basin, the DNRC may management decisions. Other interest groups could not grant new surface-water rights except in re- also bene it from the study’s baseline data in their stricted circumstances. This closure, combined future work, such as watershed health improve- with increased irrigation demands, resulted in an ment projects. increased number of high-volume irrigation wells. However, a 2006 Montana Supreme Court decision Report Purpose recognized impacts to stream low by pre-stream This report provides documentation of the capture of tributary groundwater, which effectively procedures and assumptions inherent in the model closed the basin to new groundwater develop- and communicates the indings of the model; it is ment (Montana Supreme Court, 2006). In 2007, the intended to allow the model to be evaluated and Montana Legislature passed House Bill 831 in an used by others. All iles needed to operate the attempt to improve this situation; the bill required groundwater model are posted to the program a hydrogeologic assessment to determine whether website (http://www.mbmg.mtech.edu/gwip/ a new well would result in a “net depletion” of gwip.asp). The iles are intended to enable quali ied surface water and affect a prior appropriator. If individuals to use the model developed by GWIP to depletion of surface water was shown, the appli- test speci ic scenarios of interest, or to provide a cant would then have to submit a plan to mitigate starting point for site-speci ic analysis. the depletion. General Setting Applications for new well permits have led to con licts between senior surface-water rights hold- The Lower Beaverhead study area is located in ers and junior groundwater rights holders. A com- southwestern Montana between Dillon and Bea- mon objection to the well permits is that ground- verhead Rock ( igs. 1, 2) and covers approximately water withdrawals will reduce stream low and 110 square miles. The study area encompasses lower groundwater levels. Based on House Bill 831, the lower portion of the Beaverhead River and applications for new well permits must address its loodplain, as well as the East Bench and West stream depletion and groundwater drawdown in a Bench ( ig. 2). scienti ically sound manner. Climate The Lower Beaverhead groundwater investiga- The long-term average annual precipitation in tion (Abdo and others, 2013) was conducted to Dillon is 13.17 in based on a 111-year period of build on previous hydrologic and hydrogeologic record. The short-term annual average is 11.46 in studies in the area [Montana Bureau of Mines and based on the past 30 years [Western Regional Cli- Geology (MBMG), 2008; Weight and Snyder, 2007; mate Center (WRCC), 2011]. In general, precipita- U.S. Bureau of Reclamation, 2008; H. Sessoms and tion was above average from 1900 to 1930 ( ig. 3). J. Bauder, written commun., 2007]. The primary From the 1930s through 2007, most of the annual objectives were to determine if groundwater draw- precipitation was below the long-term average; down and stream depletion were occurring due to only 17 years of that 77-year period showed above- high-capacity irrigation well withdrawals, and to average precipitation. With the exception of the evaluate possible impacts to streams from future past 2 years, the last decade has seen below-aver- groundwater development. Particular attention age precipitation, and in 7 of those years the devia- was given to the West Bench area, in part because tion below normal was 3 in or more. Nearly half stream depletion was examined on the East Bench of the annual rainfall fell between April and July 1 Butler and Abdo, 2013
Figure 1. The lower Beaverhead River basin is located in southwest Montana, approximately 23 miles northeast of Clark Canyon Reservoir. The basin is bordered by the Blacktail, Pioneer, and Ruby Mountains. of each year. The average maximum temperature (WRCC, 2011). Almost 90% of this snow fell be- over the period of record occurred in July (83.3°F), tween November and April of each year. A SNOTEL and the average minimum temperature occurred in (SNOpack TELemetry) station (Mule Creek Station January (12.6°F). #656) is located about 20 miles northwest of Dillon The average annual snowfall in Dillon is 37.3 in, at an elevation of 8,300 ft. Snow-water equivalent based on the 116-year period of record from the data (31-year record) indicate the annual aver- University of Montana–Western weather station age maximum to be 17.36 in at this site (SNOTEL,
2 Montana Bureau of Mines and Geology Open-File Report 638
Figure 2. Dillon marks the southern tip of the Lower Beaverhead study area, and Beaverhead Rock marks its northern tip. Three highways pass through the area, as well as two major canals. Three of the Beaverhead River’s tributaries (Blacktail, Stone, and Spring Creeks) originate from outside the study area. 2011). The average date of the maximum snow- the Tendoy, Snowcrest, and Blacktail Ranges to the water equivalent is May 11th, and the average date south ( ig. 1). Major tributaries to the Beaverhead of snowpack disappearance is June 15th. Therefore, River include Grasshopper Creek, Blacktail Deer the snowpack takes an average of 35 days to melt Creek, and Rattlesnake Creek. Grasshopper Creek completely. lows toward the southeast and joins the Beaver- Physiography head River above Barretts Diversion. Blacktail Deer Creek lows to the northwest in a northwest–south- The Beaverhead River drainage encompasses an east-trending valley that is nearly at right angles area of about 2,895 square miles below the Clark to the Beaverhead River Valley; it joins the Beaver- Canyon Reservoir, 23 miles southwest of Dillon, head River near Dillon. Rattlesnake Creek lows to- Montana ( ig. 1). The reservoir receives water from wards the southeast and also joins the Beaverhead Red Rock River and Horse Prairie Creek. The Bea- River near Dillon. Within the study area, extending verhead River lows northeast through the Beaver- north of Dillon to Beaverhead Rock ( ig. 2), Stone head Canyon and into the Beaverhead River Valley Creek and Spring Creek low northwest into the for about 45 miles until its con luence with the Big Beaverhead River from the Ruby Mountains. Hole and Ruby Rivers near Twin Bridges. There it In the Dillon area, the valley is about 2 miles forms the headwaters of the Jefferson River, a tribu- wide and expands to a maximum width of about tary to the Missouri River. 3 miles in the central part of the study area. The The basin is bounded by the Pioneer Mountains loodplain is bounded on the east and west by to the west, the Ruby Mountains to the east, and thick sequences of sediments that form benches, 3 Butler and Abdo, 2013
Figure 3. Annual precipitation has been below the long-term average for most of the past 80 years. Since 1930 there have only been three times when the annual total precipitation has been above the long-term average for 2 consecutive years. Model Objectives which are referred to in this report as the East and The primary objective of groundwater modeling West Benches. The East Bench is about 130 to 260 was to investigate the degree to which groundwater ft above the loodplain, and the West Bench rises withdrawals from high-capacity wells on the West about 80 to 160 ft above the loodplain. At Beaver- Bench lead to stream depletion of the lower Bea- head Rock, the loodplain is constricted by bed- verhead River and its tributaries. It was also used rock and narrows to about 1,000 ft wide. The river to evaluate the degree to which additional canal valley ranges in elevation from 5,100 ft in Dillon to seepage would offset such stream depletion. The about 4,800 ft near Beaverhead Rock. model thus served as a predictive tool. Man-Made Hydrologic Features Although the model area encompasses the east and west sides of the river, modeling efforts were Two main irrigation canals divert water from focused on the West Bench to simulate pumping the Beaverhead River to the East and West Bench- scenarios in the volcanic rock aquifer and the Ter- es: the East Bench Canal, operated by the East tiary sediment aquifer. The effects of pumping were Bench Irrigation District, and the West Side Canal, evaluated for Black Slough, Willard Slough, Albers operated by the Clark Canyon Water Supply Com- Slough, and the Beaverhead River ( ig. 4). The East pany ( ig. 2). The 53-mile East Bench Canal was Bench was the focus of an MBMG investigation completed in 1964, and provides full irrigation (MBMG, 2008) for which groundwater modeling il- service to 21,800 acres and supplemental service lustrated the effects of pumping proximal and distal to 28,000 acres (U.S. Bureau of Reclamation, 2008) to the Beaverhead River; therefore, pumping on the on the East Bench. The canal diverts water at Bar- East Bench was not addressed in this study. retts Diversion Dam, 11 miles downstream from the Clark Canyon Dam. The full capacity of the canal CONCEPTUAL MODEL is 440 cfs, and it extends about 21 miles through the study area. The West Side Canal supplies water A conceptual model is an interpretation or to about 6,855 acres on the West Bench and has a working description of the characteristics and capacity of approximately 160 cfs. This canal is di- dynamics of the physical groundwater low system. verted in Dillon and runs about 14 miles north until It is based on the analysis of all available hydrogeo- it terminates within the study area. logic data for the study area. The conceptual model includes the system’s geologic framework, aquifer properties, groundwater low directions, loca- tions and rates of recharge and discharge, and the 4 Montana Bureau of Mines and Geology Open-File Report 638
Figure 4. The Beaverhead River and three of its tributaries (Albers Slough, Black Slough, and Willard Slough) were the focus of the stream depletion modeling scenarios. All three are located on the West Bench. locations and hydraulic characteristics of natural Pioneer, Ruby, Tendoy, Snowcrest, and Blacktail boundaries (ASTM, 2004; Mandle, 2002). Ranges that border the Beaverhead River Basin is Geologic Framework composed of crystalline metamorphic rock. The structural controls in the Beaverhead River Basin A sound understanding of the existing geologic include the northeast-trending Ruby Fault Zone units and their distribution in three dimensions is along the basin’s southeast side, and the northeast- critical to a proper understanding of groundwater trending faults within the river valley (Ruppel and low. Ruppel and others (1993) provided descrip- others, 1993). Near Beaverhead Rock, the north- tions of the geology in the Lower Beaverhead study west-trending McCartney Fault Zone cross-cuts the area. Additional descriptions of the geologic units lower end of the basin ( ig. 5). In this area, faulting in this area were provided by Vuke and others has brought the Madison Limestone (Mm) to the (2007; ig. 5). These data were supplemented with surface, constricting the loodplain. Mesozoic rocks hydrogeologic concepts presented in the MBMG consisting of mudstone, siltstone, and limestone 2008 Lower Beaverhead investigation (MBMG, are also exposed in this area. The July 2005 Dillon 2008), Weight and Snyder (2007), and the Upper earthquake and other recent seismic activity indi- Beaverhead Basin hydrogeologic study (Uthman cate that some of the faults in the basin are active and Beck, 1998). (M. Stickney, MBMG Earthquake Studies Director, Most of the bedrock associated with the oral commun., 2011). 5 Butler and Abdo, 2013
Figure 5. The study area geology ranges from Mississippian to Quaternary in age. Most of the bedrock surrounding the valley consists of crystalline metamorphic rock. The basin-fi ll consists of Tertiary deposits such as the Renova and Six Mile Creek Formations. The relatively young and shallow Quaternary deposits are composed of clays, silts, sands, and gravels. 6 Montana Bureau of Mines and Geology Open-File Report 638 The basin- ill between Dillon and Beaverhead Hydrogeologic Units Rock may be about 1,000 ft thick (R. Thomas, Pro- The Lower Beaverhead basin’s hydrogeologic fessor of Geology, University of Montana Western, units were discerned from the geologic framework, oral commun., 2011). The bulk of this basin- ill con- ield observations, and well logs. Well logs from the sists of Tertiary deposits, namely the Renova and MBMG Groundwater Information Center (GWIC) Six Mile Creek Formations of the Bozeman Group, include such information as the well location, which are the main Tertiary units in southwestern lithologic descriptions, and well-completion details. Montana. The Renova volcanic and volcanoclastic The GWIC logs were reviewed to help identify the sequence was deposited during the early to middle primary types and locations of aquifers in the study Miocene in a broad continuous wedge across the area. Lithologic descriptions in each log were com- basin over a low-relief loodplain. The depositional pared with those in surrounding well logs and with sequence was overwhelmed with volcanically geologic maps. derived sediment, which typically is ine-grained The three principal aquifers identi ied in the material with low permeability. Non-volcanic facies study area are: the alluvial aquifer that forms the include sandstone, carbonaceous shale, lignite, and sur icial deposit in the Beaverhead River Valley; the limestone deposited in lakes and streams (Alt and Tertiary sediment aquifer that underlies the alluvi- Hyndman, 1986). um and also comprises the East and West Benches; During the middle Miocene, the basin was seg- and the volcanic rock aquifer on the northern West mented into several grabens by basin-and-range Bench ( ig. 6). Outcrops of the alluvium, Tertiary style faulting. Sequences of non-volcanic and vol- sediments, and volcanic rock are shown in igure canic sediments known as the Six Mile Creek For- 7. The geologic map for the area ( ig. 5) indicates mation illed the Beaverhead and other grabens in Quaternary sands and gravels (Qgr) overlying the southwest Montana during the middle Miocene to Tertiary sediments on the West Bench. The Qua- late Pliocene. The Six Mile Creek Formation is gen- ternary/Tertiary contact on the West Bench is not erally coarser-grained than the underlying Renova well de ined, and sediments deposited during this Formation and consists of mudstone, siltstone, period in time were probably formed under similar conglomerate with local occurrences of limestone, conditions. For this reason, the Quaternary sedi- volcanic fallout ash, pyroclastic ash low tuffs, ments that blanket the West Bench are considered fallout tuffs, and basalt lows. The Six Mile Creek part of the Tertiary sediment aquifer within the Formation is generally thickest near the axis of the study area. valleys and thins in the uplands. Field observations and GWIC well logs have Also during the Tertiary Period, volcanic low suggested the presence of a fourth aquifer located deposits formed as magma intruded through older near the study area’s western boundary, composed rocks and lowed at the surface. The sur icial de- of sedimentary bedrock. Based on an outcrop posits inter ingered with and were overlain by the ( ig. 8) and well log descriptions, it appears to be Tertiary sediments. The unit consists mainly of well-consolidated sandstone, thus distinguishing rhyodacite, an extrusive volcanic rock, and ranges it from the surrounding unconsolidated Tertiary from purple-brown to gray-brown in color. Within sediments. Analyses of outcrops and thin sections the lower Beaverhead basin, it is found west of Dil- indicate the formation is likely a Tertiary conglom- lon and in northern portions of the West Bench. erate in a sand matrix, with Belt quartzite clasts During the Quaternary Period, clays, silts, sands, and chalcedony veinlets (R. Berg, geologist, MBMG, and gravels were deposited over the Tertiary for- oral and written commun., 2010). Groundwater mations of the Beaverhead basin. These deposits monitoring and GWIC logs of deep wells in this area are primarily found in the valley bottoms as allu- revealed lower groundwater levels and well yields, vium, and in the uplands as alluvial fans and land- which suggest the unit is not part of the volcanic slide deposits. rock aquifer or the Tertiary sediment aquifer and is disconnected from the more productive (and well-used) hydrogeologic units in the study area. Therefore, for the purposes of this investigation, the sedimentary bedrock unit was not considered 7 Butler and Abdo, 2013
Figure 6. This schematic of a geologic cross section illustrates the relation of the alluvial, Tertiary sediment, and volcanic rock aquifers to one other. The cross section line is shown in fi gure 5 (A–A’). to be one of the principal aquifers. Based on these ield observations and GWIC The following discussion characterizes the well logs, it was postulated that a con ining clay lay- three aquifers in terms of thickness and lithologic er may underlie the alluvial aquifer. Con ining units, composition. Quantitative hydrogeologic proper- such as clay layers, can signi icantly affect ground- ties, such as transmissivity and storativity, are water low paths. Understanding the continuity provided in the Aquifer Properties section. and extent of less permeable layers is important to predicting groundwater/surface-water interactions Alluvial Aquifer and groundwater drawdown from pumping. The Based on drill cuttings observed in the ield extent of the clay layer under the shallow alluvial and in other drillers’ logs, about 25 to 30 ft of silts, aquifer was investigated by mapping drillers’ logs sands, gravels, and cobbles underlie the surface from the GWIC database to create geologic cross of the valley loor. This thickness varies and can sections. This approach provided a means to visu- extend to greater depths. This coarse material con- ally identify the extent and continuity of the clay tains cobbles up to 3 in. in diameter. Depth to water across the loodplain. Recorded well locations were in this aquifer ranges from approximately 3 to 13 ft con irmed by comparing them to the landowners’ below ground surface. reported locations. Well logs in which lithologic During well installations in the northern lood- descriptions were missing or were determined to plain area ( ig. 9), a gray silty clay layer was identi- be inaccurate were eliminated from the evaluation. ied underlying the shallow alluvial aquifer (MBMG, The elevation for each well was obtained by using 2008). Beneath the clay layer was gray indurated online elevation- inder software with an accuracy silt identi ied as Tertiary sediments. In the central of ±5 ft (GPS Visualizer, 2011). Limitations to this loodplain ( ig. 9), light brown clay and silty clay approach included the inaccuracy associated with were noted in two intervals: at 22 to 50 ft and 65 to wells mislocated within landowner parcels, corre- 85 ft below ground surface. Indurated silt encoun- lation errors due to geologic discontinuities, a lack tered at 85 ft below ground surface was considered of wells that fully penetrate the shallow alluvium in to be the top of the Tertiary sediments. GWIC well some areas, and ambiguous geologic descriptions logs in the area also indicated the presence of clay provided in the drillers’ logs. or other less permeable units underlying the allu- The resulting cross sections illustrated that the vium. reported clayey materials underlying the alluvium 8 Montana Bureau of Mines and Geology Open-File Report 638
Figure 7. These photographs show the three principal aquifers in the study area. Note the fi ne-grained nature of the Tertiary sediments. Secondary porosity in the volcanic rock makes it the most prolifi c aquifer in the study area.
9 Butler and Abdo, 2013
Figure 8. A sedimentary bedrock outcrop, observed in August 2010, is located in the northern portion of the West Bench. appear to be continuous in some areas and less to be the top of the Tertiary. Farther south, (well so in others ( ig. 10). In the more continuous ar- 255492, T. 6 S., R. 8 W.) the Quaternary–Tertiary eas, these materials may effectively form localized contact was less distinct. Several clay layers were con ining beds separating the alluvial and Tertiary encountered, with the thickest being about 20 ft at sediment aquifers; in contrast, the discontinuities 65 to 85 ft below ground surface. Below 85 ft, indu- allow for a direct connection between the alluvium rated silt was encountered and considered Tertiary and Tertiary sediments. Because evidence of a dis- sediments. The thickness of the Tertiary sediments tinct, continuous clay unit was lacking, the clay was beneath the loodplain is uncertain. The deepest numerically modeled as part of the Tertiary sedi- well completed in these sediments is 460 ft deep, ment aquifer (see Model Grid section). and it does not fully penetrate the unit. Tertiary Sediment Aquifer The West Bench Tertiary sediments consist of ine to coarse sands, gravels, cemented gravels, The Tertiary sediments (Ts, ig. 5), which un- and interbedded clay and silt. The numerous clay derlie the alluvium in the loodplain and comprise and silt layers generally range from 1 to 12 ft thick. most of the East and West Benches, include a wide However, a 40-ft-thick silty clay (60 to 100 ft below range of lithologies. Figure 7 illustrates the over- ground surface) was identi ied during one well all ine-grained nature of the sediments. Depth to installation (well 254962, T. 6 S., R. 8 W.). Bedrock water in this aquifer ranges from approximately 3 noted in well logs indicates that Tertiary sediments to 35 ft below ground surface in the loodplain; 13 may only be about 60 ft thick in some areas of the to 127 ft on the East Bench; and 2 to 300 ft on the West Bench. The deepest wells on the West Bench West Bench. that are completed in Tertiary sediments are about During installation of well 242403 ( ig. 9, T. 5 S., 400 to 500 ft deep. R. 7 W.) in the loodplain, indurated silt and sand The East Bench Tertiary sediments are similar were observed below 65 ft and were considered to the West Bench, in that the lithology consists of 10 Montana Bureau of Mines and Geology Open-File Report 638
Figure 9. The groundwater monitoring network includes private water-supply wells and dedicated monitoring wells.
Figure 10. Geologic cross sections were generated from GWIC well logs to examine the extent of clay underlying the alluvium in the fl oodplain. 11 Butler and Abdo, 2013 interbedded gravel, sand, silt, and clay. Drill cut- Volcanic Rock Aquifer tings in the Spring Creek drainage ( ig. 9) indicate Through geologic mapping, drillers’ logs, ield that the lithology is dominated by indurated silty observations, and geophysical surveys, volcanic sand with less abundant layers of clay, sand, and rock locations have been approximated on the West gravel; clay layers ranged from about 10 to 30 ft Bench. Sur icial exposures were either mapped thick. The thickness of the Tertiary sediments is previously (Ruppel and others, 1993; Tv in ig. 5) generally much deeper than that of the West Bench, or observed in the ield during this study. Figure as noted in the drillers’ logs. The maximum thick- 11 illustrates the horizontal extent of the volcanic ness of the unit is uncertain. The deepest well rock aquifer based on these data. The rhyodacite completed on the East Bench was drilled to nearly is vesicular and exhibits secondary porosity as a 700 ft, and it did not fully penetrate the Tertiary result of fracturing and dissolution of phenocrysts. sediments.
Figure 11. Location of volcanic rock outcrops and wells in which volcanic rock was noted in the driller’s log. Well logs and drill cuttings indicate the rock is of variable thickness (3–200 ft) and interbeds with the Tertiary sediments. The full lateral extent and maximum thickness of the volcanic rock are unknown. 12 Montana Bureau of Mines and Geology Open-File Report 638 The volcanic rock of the West Bench is believed to 2013) suggest that groundwater low within the be older than or contemporaneous with the Ter- study area is irrigation-driven, with most (approxi- tiary sediments, and well logs indicate that it both mately 88%) of the groundwater recharge derived underlies and inter ingers with the sediments. The from canal seepage and water applied to irrigated volcanic low deposits appear to be variable in ields. Mapped potentiometric surfaces ( igs. 12, thickness, ranging from a few feet to 200 ft in well 13) show additional in low coming from the up- log descriptions; the maximum thickness is un- lands along the east and west edges of the study known. Depth to water in this aquifer ranges from area, and from the upper Beaverhead basin south of approximately 10 to 35 ft below ground surface. Dillon; local precipitation is also a groundwater re- Groundwater Flow System charge component. However, the dominant control on the low system is water originating from irriga- Water-balance calculations (Abdo and others, tion return lows, as evidenced in hydrographs of
Figure 12. Potentiometric surface map of the Tertiary sediment aquifer, based on data collected in August 2010. Groundwater fl ows from the East and West Benches towards the Beaverhead River Valley. 13 Butler and Abdo, 2013
Figure 13. Water-table map of the alluvial aquifer, based on data collected in August 2010. Groundwater fl ows to the northeast, from Dillon toward Beaverhead Rock. wells within the lower Beaverhead basin’s irrigated least one reach of the river might be losing just areas ( ig.14); they show an annual pattern of a north of Dillon, it is primarily a gaining stream steady rise from late spring through late summer, within the study area. Likewise, the river’s tribu- followed by a plateau in the fall, and then a decline taries are believed to be almost entirely gaining. through early spring. Another major groundwater discharge mechanism The majority of groundwater discharge from in the lower Beaverhead basin is evapotranspira- the study area occurs as surface-water out low to tion (ET), notably on irrigated lands. Discharge also the Beaverhead River and its tributaries. While occurs through irrigation well withdrawals, while the alluvial water-table map ( ig. 13) indicates at domestic and public water supply wells withdraw 14 Montana Bureau of Mines and Geology Open-File Report 638
Figure 14. Water-level hydrographs from wells on irrigated land in the Lower Beaverhead study area. The groundwater-level fl uctuations refl ect the infl uence of irrigation recharge. Well 109060 is located in sec. 25, T. 6 S., R. 9 W., and well 259540 is located in sec. 29, T. 6 S., R. 7 W
15 Butler and Abdo, 2013 relatively small amounts (Abdo and others, 2013). loodplain is due to the higher transmissivity and Lastly, groundwater low discharges at the down- geometry of the alluvial aquifer. The slightly steep- gradient pinch point of the study area, near Bea- er gradient of the East Bench vs. the West Bench verhead Rock; due to the shallow bedrock at this could be caused by an overall higher transmissivity location (Geologic Framework section), most of the of the aquifers on the West Bench, the higher leak- low is believed to move into streams and through age rate of the East Bench Canal as compared to the shallow alluvium. Volumetric estimates of wa- the West Side Canal (Groundwater Budget section), ter moving through the system are provided in the and/or topographic controls. Groundwater Budget section. Vertical Gradients The three units (alluvial, Tertiary sediment, and volcanic rock aquifers) are hydrogeologically con- Vertical groundwater gradients were evalu- nected. The volcanic rock aquifer is distinct in that ated in 2010 in three areas: between the alluvial groundwater moves through fractures rather than and Tertiary sediment aquifers in the loodplain; porous media. However, the fractures are so ex- within the Tertiary sediment aquifer on the East tensive that, when viewed at the scale of the study and West Benches; and between the volcanic rock area, they can be treated as equivalent porous and Tertiary sediment aquifers on the West Bench. media. As described below, both upward and downward Within the study area’s aquifer system, ground- gradients were observed. water lows from the benches to the loodplain ( ig. In the loodplain, Site A ( ig. 15) included a 12). The groundwater gradient on the East Bench shallow (23.5 ft) and a deep (108 ft) monitoring is about 0.01. The gradient is slightly gentler on well. During the 2010 irrigation season, an upward the West Bench, about 0.009, and lattens to about gradient occurred between the deeper Tertiary 0.006 in the loodplain. The lower gradient in the sediment aquifer and the alluvium ( ig. 16). A small
Figure 15. Seven sets of surface-water and groundwater monitoring sites were used to evaluate vertical gradients. Both upward and downward gradients were observed in the data. 16 Montana Bureau of Mines and Geology Open-File Report 638 upward gradient was also observed between the al- well pair showed a consistent downward gradient. luvium and the river during part of this time; in the Another well pair at Site F, located in the up- fall and winter the gradient decreased and even- per reach of Black Slough and near the West Side tually reversed to a downward gradient from the Canal, was completed in the Tertiary sediments at river to the groundwater. Site B ( ig. 15) featured a 77 ft deep (well 254839) and in the volcanic rock shallow (24 ft) and a deep (94 ft) monitoring well. at 158 ft deep (well 254815). The hydrographs for The water-level elevations in the shallower aquifer these wells ( ig. 17) indicate the gradient changed were about 4 ft higher than those of the deep moni- seasonally. During the 2010 irrigation season, a toring well, thus indicating a downward gradient downward gradient of about -0.011 occurred from between the aquifers. the Tertiary sediments to the volcanic rock. Outside Three sets of well pairs on the West Bench of the irrigation season in 2010, the gradient was (Sites C, D, and E; ig. 15) show a downward gra- most often upward (0.017) from the volcanic rock dient. However, these well pairs were located to the Tertiary sediments. between 85 and 715 ft apart, and these distances are great enough to lend uncertainty in geologic discontinuities that could affect the accuracy of the vertical gradients. On the East Bench at Site G, the
Figure 16. Water elevations in the Beaverhead River and nearby monitoring wells at Sites A and B (see fi g. 15 for locations). At Site A, water elevations indicate the river can gain or lose water from groundwater depending on the time of year. At Site B, the river consis- tently gained water from the alluvial aquifer. 17 Butler and Abdo, 2013
Figure 17. Hydrographs in adjacent monitoring wells 254839 (completed in Tertiary sediments) and 254815 (completed in volcanic rock) indicate a downward gradient from the Tertiary sediments to the volcanic rock during the irrigation season, and an upward gradient during the non-irrigating season of 2010. Hydrologic Boundaries Aquifer Properties The most signi icant hydrologic boundary Sources of data for the aquifer properties in within the study area is the Beaverhead River. It the vicinity of the study area include aquifer test lows through the center of the valley and receives data available from water-rights permit reports groundwater recharge from both the West and East obtained from the Montana DNRC, aquifer tests Benches, thus acting as a groundwater sink ( ig. conducted as part of this study, previous hydrogeo- 12). Both the East Bench and the West Side Canals logic reports for the lower Beaverhead River basin are hydrologic controls on the groundwater low (MBMG, 2008), and values from similar groundwa- system, in that canal seepage has a strong in luence ter studies and low models in western Montana. on the groundwater gradient and groundwater Aquifer properties typically estimated from aquifer recharge. tests are transmissivity, hydraulic conductivity (K), The pinch point near Beaverhead Rock is anoth- and storativity (S). The range of parameter values er important hydrologic feature; here, groundwater exhibited in the alluvial, Tertiary sediment, and low is constricted and forced to emerge as springs, volcanic rock aquifers were established using these as base low to streams, or as groundwater low available data. Aquifer tests conducted speci ically through the narrow outlet from the Lower Beaver- for this investigation are brie ly discussed below. head valley. Other features of signi icance are the Tertiary Sediment Aquifer Test Beaverhead River watershed boundaries, which are Irrigation well 220021 was pumped for 3 days assumed to act as groundwater divides. The north- (May 18–21, 2010) at 300 gpm. Water levels were ern border of the study area approximates one such monitored in two monitoring wells (wells 254962 watershed boundary ( ig. 1). Those located outside and 254963), two shallow domestic wells (wells of the study area include the ridgeline of the Pio- 108978 and 258390) and at two locations in Wil- neer Mountains to west, and the ridgeline of the lard Slough. All wells were completed in the Tertia- Ruby Mountains to the east. ry sediments. Figure 18 shows the locations of the monitoring sites. Maximum drawdown in wells 254962 and 18 Montana Bureau of Mines and Geology Open-File Report 638
Figure 18. Tertiary sediment aquifer test site. Monitoring during the aquifer test of irrigation well 220021 included both wells and surface-water sites.
254963 were about 42 ft and 34 ft, respectively. value indicates that the aquifer is con ined in the Drawdown in the pumping well was about 230 ft. vicinity of the test site. There was no observed drawdown in the shallow Volcanic Rock Aquifer Test domestic wells or any apparent in luence to sur- face-water low at either Willard Slough monitoring Irrigation well 220080, completed in the vol- location. canic rock aquifer, was pumped from October 15 Transmissivity estimates were derived us- to 18, 2010 at an average rate of 1,420 gpm. The ing the Cooper–Jacob (1946) composite plot and monitoring network included seven wells and one straight-line analyses. Estimates generated from surface-water site in Black Slough ( ig. 19). Five the two methods agree (412 ft2/day and 405–522 of the seven monitoring wells (204226, 224244, ft2/day, respectively). A storativity of 0.00098 was 254767, 254815, and 254840) were completed in calculated using the Cooper–Jacob composite plot the volcanic rock aquifer, one (254839) was com- method (Abdo and others, 2013). This storativity pleted in the Tertiary sediment aquifer, and one (259558) was a shallow (16 ft deep) piezometer 19 Butler and Abdo, 2013
Figure 19. Volcanic rock aquifer test site. Monitoring during the aquifer test of irrigation well 220280 included seven wells and one surface-water site. completed in the alluvial channel of Black Slough. aquifer test analysis software, AQTESOLVTM. The Drawdown due to pumping was observed in Cooper–Jacob composite plot, straight-line, and all monitoring wells, though the water level in the distance drawdown methods (Cooper and Jacob, shallow piezometer (259558) was not affected un- 1946) were used to estimate transmissivity and til after pumping ceased. The maximum drawdown storativity. Transmissivity estimates ranged from in the pumping well (220080) was 4 ft. 42,500 to 75,000 ft2/day, and storativity values The van der Kamp (1989) method was used ranged from 0.0026 to 0.016 (Abdo and others, to extend the drawdown data because water lev- 2013). The storativity values indicate the aquifer els in the pumped well did not reach steady-state. ranges from uncon ined to semi-con ined. The drawdown data were then analyzed using the 20 Montana Bureau of Mines and Geology Open-File Report 638 Alluvial Aquifer Tests Sources and Sinks While tests were not conducted within the allu- As noted in the Groundwater Flow System sec- vial aquifer as part of this investigation, tests were tion, sources of groundwater recharge within the performed previously by the MBMG (MBMG, 2008) Lower Beaverhead study area include seepage from and by consultants for well permit applications canals and ditches; water applied to irrigated ields; (Water Rights, Inc., written commun., 2007). Test local precipitation; leakage from the Beaverhead results indicated the transmissivity of the alluvial River and its tributaries; and groundwater in low aquifer ranges from about 18,000 to 37,000 ft2/day, from the upper Beaverhead River basin as well as and storativity ranges from 0.003 to 0.15 (table 1). the eastern and western uplands. The sinks (i.e., Summary of Hydraulic Properties groundwater discharge mechanisms) in the Lower Beaverhead study area include surface-water out- Table 1 provides a summary of aquifer property low to the Beaverhead River and its tributaries; estimates based on aquifer test results from this evapotranspiration; groundwater low exiting the study, previous studies, and water permit applica- valley at Beaverhead Rock; and well withdrawals, tions within the study area. Values used for various most notably from irrigation wells. The next section aquifer types in some other groundwater modeling provides estimates of groundwater lux for these efforts in western Montana are shown in table 2. sources and sinks. This table includes the model presented by Uth- man and Beck (1998), in which aquifer properties Groundwater Budget were based on aquifer tests conducted in the upper A groundwater budget quantitatively summa- Beaverhead basin’s alluvial and upper Tertiary rizes the processes within the conceptual model. aquifers. While there is inherent uncertainty associated with The ranges of aquifer property values were the calculations, a groundwater budget is useful for used in groundwater low calculations and ground- determining the relative importance of different water modeling; they are all within the broad processes. A groundwater budget combines water ranges of expected values as described in numer- entering and leaving the study area from sources ous groundwater textbooks for similar materials. and sinks, respectively. The idea of a water budget The groundwater modeling effort further assessed is the same as the more general law of mass bal- these aquifer property estimates. Table 3 compares ance; that is, matter cannot disappear or be cre- the ield-estimated K and S values in table 1 with ated spontaneously. Thus, the amount of water that those used in the model. enters over a period of time must be equal to the
7DEOH$TXLIHUSURSHUWLHVIRUWKHDOOXYLXP7HUWLDU\VHGLPHQWVDQGYROFDQLFURFNXVHGLQWKLVVWXG\
7\SLFDO 6DWXUDWHG 7UDQVPLVVLYLW\ +\GUDXOLF 7KLFNQHVV $TXLIHU7\SH IWGD\ 6WRUDWLYLW\ *UDGLHQW IW 6RXUFH $OOXYLXP ± ± GD\DTXLIHUWHVW
7HUWLDU\VHGLPHQWV ± GD\DTXLIHUWHVW :HVW%HQFK 7HUWLDU\VHGLPHQWV ± :DWHUULJKWVSHUPLWV (DVW%HQFK 7HUWLDU\VHGLPHQWV GD\DTXLIHUWHVWZDWHUULJKWV 15 IORRGSODLQ LQIRUPDWLRQ
9ROFDQLFURFN ± ± 1$ GD\DTXLIHUWHVW
3XPSLQJ:HOO 0%0* 3XPSLQJ:HOO 0%0*WKLVUHSRUW 8QSXEOLVKHGGDWD0RQWDQD3URYLVLRQDO:DWHU5LJKWV3HUPLW%ZULWWHQFRPPXQ:DWHU5LJKWV,QF0LVVRXOD 0RQWDQD 3XPSLQJ:HOO 0%0* XQSXEOLVKHGGDWDZULWWHQFRPPXQ:DWHU5LJKWV,QF0LVVRXOD0RQWDQD 3XPSLQJ:HOO 0%0*WKLVUHSRUW 21 Butler and Abdo, 2013 Table 2. Hydraulic conductivity values used in other western Montana large-area groundwater models.
8SSHU%HDYHUKHDG²8WKPDQ %HFN '15&WUDQVLHQWVLPXODWLRQ Layer 1 (Quaternary alluvium), 25 to 170 ft thick 10–1,800 ft/day Layer 2 (Tertiary basin-fill) 5–10 ft/day
+D\HV&UHHN²:DUHQ '15&/LPLWHG7UDQVLHQW Belt Argillite–Missoula Group, Mount Shields Fm., Member 3 0.1–0.75 ft/day
+HOHQD9DOOH\$TXLIHU²%ULDUDQG0DGLVRQ 86*6VWHDG\VWDWHVLPXODWLRQ Layer 1 Upper 35 ft thickness of aquifer 80 ft/day Layer 2 Next 75 ft thickness 40 ft/day Layer 3 170 to 1000 ft thickness beneath layers 1 and 2 40 ft/day
/RZHU%HDYHUKHDG²0%0* WUDQVLHQWVLPXODWLRQ Alluvium 75 ft/day Mesozoic bedrock 1–2 ft/day
*DOODWLQ9DOOH\DQG0DGLVRQ3ODWHDX²0%0* WUDQVLHQWVLPXODWLRQ Alluvium 82–131 ft/day
'UXPPRQG9DOOH\².DXIIPDQ 0RQWDQD6WDWH8QLYHUVLW\*UDGXDWH7KHVLVWUDQVLHQWVLPXODWLRQ Alluvium 26–45 ft/day
Table 3. Comparison of field-estimated and model aquifer property values.
Field K Range1 Model K Range Field Storativity Model Storativity (ft/day) (ft/day) Geometric Range (ft/day) Range* Mean2 Min Max Min Max Min Max Min Max Floodplain—Alluvium 600 1,233 110 1,800 764 0.003 0.15 0.15 Floodplain—Tertiary sediments 0.87 6.09 45.2 134 95.5 NR 0.08 0.10 East Bench—Tertiary sediments 3.53 6.43 3.74 8.92 5.74 0.004 0.05 0.08 West Bench—Tertiary sediments 0.87 1.13 1.50 16.5 6.47 0.001 0.05 0.08 Northern West Bench—Volcanic 142 252 433.5 433.5 0.003 0.018 0.15 bedrock Other bedrock — — 0.12 0.15 0.14 NR 0.01 0.05 1K values documented in Hydrogeologic Setting section of this report. 2This denotes the mean of the modeled PEST zones that represented the given aquifer type. *Storativity was not estimated using PEST. S values are discussed in the Transient Calibration section. Note. NR, not reported.
amount of water that leaves over that same time GI + CL + IR = BR + GD + WW ± S, (1) period, plus or minus any change in aquifer stor- age; in a groundwater system, changes in storage where: are directly related to changes in groundwater GI, groundwater in low along the study area levels. The general form of the mass balance equa- boundary; tion is: CL, irrigation canal leakage; Inputs = Outputs ± Changes in storage. IR, irrigation recharge; BR, base low to Beaverhead River and its The mass balance equation can be expanded tributaries (net loss to streams); for the Lower Beaverhead study area to include the GD, groundwater discharge through components below; in lows are presented on the Beaverhead Rock pinch point; left side of the equation and out lows are presented WW, withdrawal from wells; and on the right. S, changes in storage.
22 Montana Bureau of Mines and Geology Open-File Report 638 These components were based on the study area at the south edge of the model. The total in low hydrogeology and land uses that affect groundwa- from boundaries was an estimated 6,068 acre-ft/yr. ter. Certain components re lect the net result of a Leakage from irrigation canals (CL) provides a combination of sub components. For instance, the source of recharge to the groundwater. GIS analysis irrigation recharge component accounts for both of canals showed that 12.6 miles of the East Bench precipitation and evapotranspiration on irrigated Canal and 14.3 miles of the West Side Canal are lands. The following discussion summarizes the present in the modeled portion of the study area. water budget components used in the model as Leakage is quite evident in several hydrographs shown in equation 1. All component estimates have ( ig. 14). During the summer of 2010, a series of inherent uncertainty, some to a greater degree than low measurements were made along each canal, others. For this reason the model components were and diversion records for each day of the low not calibrated precisely to ield estimates; rather, measurements were obtained. Based on these data, the ield estimates were used as guides to which seepage rates were estimated for several reaches of the model components were approximated. each canal, and an average value per canal was cal- Groundwater in low (GI) enters the study area culated (Abdo and others, 2013). The average seep- from upgradient areas along the East Bench, West age rates used as initial model inputs for the East Bench, and in the loodplain at Dillon (i.e., the Bench Canal and West Side Canal were 2.8 cfs/mile upgradient border of the Lower Beaverhead Val- and 1.1 cfs/mile, respectively. Using this approach, ley). Along the study area’s East Bench boundary, a total of 18,100 acre-ft/year was estimated to in il- groundwater moves through the Tertiary sediment trate from the canals during the irrigation season. aquifer, which is relatively thick and continuous. When water (irrigation water plus precipita- In contrast, in low from the West Bench bound- tion) is applied to a ield in excess of crop demand ary is derived from the Tertiary sediment aquifer and evaporation, the excess must either runoff or as well as sedimentary rock and volcanic rock ( ig. in iltrate to the subsurface. On irrigated ields, the 5); the bedrock aquifers are believed to provide water that recharges groundwater is termed irriga- less recharge due to their relatively low hydraulic tion recharge (IR). A linear groundwater recharge conductivity values. Along the study area boundary rate was calculated for the irrigated land within the in the loodplain, groundwater lows from the Up- study area as follows: per Beaverhead Valley into the Lower Beaverhead Recharge Rate = P + R , (3) Valley; the majority of in low is derived from the in irr highly transmissive alluvial aquifer, while a smaller where Pin is total precipitation and Rirr is the amount enters through the underlying Tertiary sed- groundwater recharge from irrigation, which was iment aquifer. The estimate of groundwater enter- calculated by the NRCS IWR method (appendix B). ing the site was calculated using groundwater low The IWR method considers the recharge rates for nets and Darcy’s law: three irrigation types: lood, pivot, and sprinkler. Q = KiA, (2) In non-irrigated areas, groundwater recharge from where: precipitation was assumed to be negligible. Using Q, groundwater low; this approach, the irrigation recharge was calcu- K, hydraulic conductivity; lated to be approximately 14,000 acre-ft/yr on the i, hydraulic gradient; and West Bench, 16,000 acre-ft/yr on the East Bench, A, low tube area (tube width multiplied by and 18,000 acre-ft/yr in the loodplain, for a total saturated aquifer thickness). of 48,000 acre-ft/yr. Base low to the Beaverhead River and its tribu- Using this approach, groundwater in low from the taries (BR) is the primary means by which ground- study-area boundaries was calculated. About 142 water exits the aquifers of the Lower Beaverhead acre-ft/yr and 3,161 acre-ft/yr were estimated to study area. Although ield data (e.g., stage, low, and low in from the boundaries of the West Bench and alluvial groundwater-level measurements) show East Bench, respectively. Approximately 2,765 acre- that the Beaverhead River loses low to the alluvial ft/yr was estimated to low into the model from the aquifer in certain reaches, it is generally a gaining alluvium and underlying Tertiary sediment aquifer stream. To estimate base low to the river, surface- 23 Butler and Abdo, 2013 water out lows were subtracted from surface-water detail on this approach is provided in appendix C. in lows. The in lows consisted of the Beaverhead The resulting well withdrawal estimate was ap- River at Dillon and its tributaries that originate proximately 8,150 acre-ft/yr. from outside the study area and join the river with- Annual water budgets often ind a balance in the study area, including Blacktail Deer Creek between the amount of water that enters and exits and Stone Creek. Surface-water out lows consisted the study area; however, changes in storage (S) in of only the river at Beaverhead Rock. The resulting the aquifer system can cause the volumes of in low base low within the study area was estimated to be and out low to differ. 2010 was the 23rd wettest between 25,000 and 35,000 acre-ft/yr. This esti- year in the past 111 years of recorded precipita- mate was of a range rather than a single value due tion in the Dillon area. Examination of groundwater to the inherent uncertainty in the approach. hydrographs revealed a gain in storage within the In the Lower Beaverhead study area, ground- study area in 2010. The volume of the change in water discharges through the pinch point at Bea- storage was estimated by comparing water levels verhead Rock ( ig. 2). Groundwater discharge from January 2010 and January 2011 for the West (GD) was estimated using Darcy’s Law, as with the Bench, East Bench, and loodplain areas. An aver- groundwater in low component. Because water age change in groundwater level was calculated well logs indicate that bedrock underlies the allu- from numerous wells within each of the areas. The vium at the pinch point, most of the groundwater average water-level change was then multiplied likely lows through the alluvium. Using Darcy’s by the estimated effective porosity for each area Law, the alluvial aquifer thickness in this location and by the total acreage. The resulting estimate was approximated to be 30 ft and was assigned a K of groundwater held in storage during 2010 was of 1,200 ft/day; the underlying bedrock had a thick- about 17,000 acre-ft. Most of this increase in stored ness of 470 ft in the model and was assigned a K of groundwater occurred on the East Bench (72%; 0.1 ft. The valley width was measured using aerial Abdo and others, 2013). photography and is about 1,000 ft wide at this loca- The model water budget results in both the tion. Last, the groundwater gradient was estimated steady-state and transient versions of the model from potentiometric maps ( igs. 12,13). The result- are presented in the Calibration section because ing discharge estimate was about 1,500 acre-ft/yr. certain budget components were adjusted during Well withdrawals (WW) are another groundwa- calibration. Results showed that irrigation (i.e., ter sink in the Lower Beaverhead study area. Only canal leakage and applied irrigation water) was irrigation wells and public water supply (PWS) the predominant source of recharge to the model, wells were simulated in the model; withdrawals which demonstrates its considerable in luence on from domestic wells were not considered due to groundwater low within the Lower Beaverhead their relatively low water usage. The number of study area. The primary means of groundwater dis- such wells within the study area was estimated by charge from the model was base low to the Beaver- inventorying irrigation and PWS wells through the head River and its tributaries. GWIC database. Water-rights records and permit applications were also inventoried through the COMPUTER CODE DNRC Water Rights Bureau and through consul- Groundwater Modeling Systems (GMS) software tants who produced hydrogeologic reports for was used to develop a MODFLOW-2000 groundwa- Dillon-area well permit applications (PBS&J, writ- ter low model (Aquaveo, 2010). MODFLOW-2000 ten commun., 2005; WET, written commun., 2004; is a widely accepted groundwater low program Water Rights, Inc., written commun., 2007); these developed by the US Geological Survey (Harbaugh records were reviewed for such data as crop acre- and others, 2000). It simulates groundwater low age, crop type, and water use allotments. An annual numerically using a inite-difference method. The withdrawal volume was then calculated for each version of GMS used for this modeling was GMS well based on this information as well as the IWR 7.1.2, with a build date of April 16, 2010. The ver- per-acre crop requirements (appendix B). Sum- sion of MODFLOW-2000 operated in GMS 7.1.2 was mertime aerial photos and transducer data (where Version 1.18.01, compiled June 20, 2008. available) were also used to verify well use. Further 24 Montana Bureau of Mines and Geology Open-File Report 638 PEST is a general purpose parameter estimation grid in the X, Y, and Z dimensions respectively are utility developed by John Doherty of Watermark 76,646, 74,335, and 500 ft. This rectangular grid Numerical Computing (Doherty, 2010). PEST was frame encompasses the Lower Beaverhead study used for automated parameter estimation in cer- area; some cells within grid frame were inactivated tain model runs. The version of PEST operated in in order for the model domain to best correspond GMS 7.1.2 was PEST Version 12.0. with the study area ( ig. 20). Cells were 200 ft x 200 ft, and the model had two layers, 215 rows, and 230 GROUNDWATER FLOW MODEL columns. The model thickness was 500 ft to ap- CONSTRUCTION proximate the portion of the aquifer in which most irrigation wells are completed (based on reported Model Grid total depths of irrigation wells, GWIC), and the The GMS project was operated using the North saturated thickness ranged from about 200 to 500 American Datum (NAD) 1983 Montana State Plane ft. Table 4 provides additional numeric details of coordinates, with vertical units of U.S. Survey Feet. the model grid. The model grid was created using a uniform grid The top of layer 1 was de ined using data de- frame with an X origin of 1,140,727 ft, Y origin of rived from the U.S. Geologicial Survey 1/3-Arc Sec- 359,905 ft, and Z origin of 4700 ft. Lengths of the ond National Elevation Dataset (USGS, 2009). These
Figure 20. Delineation of the active cells in the Lower Beaverhead model domain. The model featured a uniform grid spacing of 200-ft by 200-ft cells. 25 Butler and Abdo, 2013 Table 4. Details of the model grid as listed in GMS. Grid Type Cell Centered X origin (ft) 1,140,727 Y origin (ft) 359,905 Z origin (ft) 4,700 Length in X (ft) 76,646 Length in Y (ft) 74,335 Length in Z (ft) 500 Rotation angle 0° Minimum scalar 4,827 Maximum scalar 5,155 Number of rows (i) 372 Number of columns (j) 383 Number of layers (k) 2 Number of nodes 358,080 Number of cells 284,952 Number of active cells 146,635 Number of inactive cells 138,317 data were converted into a scatter point dataset and imported into GMS as a text ile. This scatter point set is referred to here as the Digital Elevation Model (DEM) scatter point set. The DEM scatter point spacing is about 186 ft, which is similar to the cell size of 200 ft. The bottom of layer 1 was de- ined by a surface derived from two elements. The oodplain. irst element was a surface de ined by subtracting fl 30 ft from the elevation of the DEM scatter point set. This irst approach was used in the loodplain, where the upper 30 ft of the grid approximated the alluvial aquifer. The second element was a compos- ite of lat surfaces that changed elevation westward and eastward from the loodplain onto the benches. The shifts in elevation of the bottom surface cor- respond with large shifts in elevation of the top surface (i.e., the land surface). Layer 1 was about 250 ft thick toward the eastern and western extents of the grid, where the model represents the Tertia- ry sediments and volcanic bedrock. This thickness ensured that the maximum depth to groundwater would remain above the bottom of layer 1 and prevent cells from drying ( ig. 21). In the lood- plain, layer 2 represented the Tertiary sediments underlying the alluvium; as with layer 1, layer 2 represented the Tertiary sediments on the East Bench and a combination of Tertiary sediments The grid has a uniform thickness of 500 ft and features orientation in the model grid. Figure 21. Model cross section through Row 237, which runs in a West–East two layers, with the upper layer representing alluvium within and volcanic bedrock on the West Bench. Because the model thickness was held constant (500 ft), the
26 Montana Bureau of Mines and Geology Open-File Report 638 thickness of layer 2 was variable; it was relatively into four general categories: the model borders, thick in the loodplain (approximately 470 ft) and surface-water bodies, aerial recharge (precipitation thinner at the east and west edges of the benches and applied irrigation water), and well withdraw- (approximately 250 ft). als. Figure 22 shows the boundary conditions in the Prior to arriving at this layer con iguration, model with the exception of aerial recharge, which another version was tested. The initial version is shown in igure 23. of the model grid tightly discretized the vertical Model Borders hydrogeology vertically using four layers. The goal of the four-layer design was to reproduce site-scale Speci ied- lux boundaries were used to repre- observations in the volcanic rock aquifer on the sent in low along the East and West Benches, and northern West Bench (Aquifer Properties section). from the upper Beaverhead River basin into Dil- However, the degree of detail could not be ad- lon from the south. These boundaries were placed equately calibrated via PEST due to a lack of cali- along contours based on the potentiometric surface bration targets within certain layers (Steady-State map developed for this study, and they replicate Calibration section). the relatively stable groundwater setting observed in these areas. A speci ied- lux boundary simulates Hydraulic Parameters water in low or out low to the aquifer system as Prior to model calibration, K and S values were a user-de ined volumetric rate. Flux values were assigned to polygonal zones in the model based on estimated using a low net approach (Groundwater the aquifer property estimates from aquifer tests Budget section), with hydraulic conductivity (K) performed during this investigation and previous and gradient values based on aquifer property es- investigations (Aquifer Properties section). The timates. Along the East Bench border, the lux rate polygon extents were based on the hydrogeologic represented seepage from the East Bench Canal in units of the conceptual model (Geologic Framework addition to the upland in low component, as the section); the units included the loodplain alluvium, canal was located just outside of the model domain. the volcanic rock on the West Bench, and the Ter- At the portion of the model border near Beaver- tiary sediments. The initial K values were modi ied head Rock, the loodplain constricts as groundwa- during the steady-state model calibration process. ter low exits the model domain. The Drain Package Storativity was introduced in the transient mod- was assigned to this model border. This MODFLOW el. Values were assigned in polygonal zones based package allows groundwater low to exit the mod- on estimates presented in the Hydrogeologic Setting eled aquifer. No- low boundaries were set along section. The simulation was divided into monthly groundwater low lines at the remainder of the increments to best simulate seasonal changes and north and northeast model border; these no- low to use the available monthly water-level datasets in boundaries run parallel to groundwater low lines model calibration. During the transient model cali- on the potentiometric surface map. bration, parameter values were adjusted to render Surface-Water Bodies the observed water-level luctuations. A speci ied- lux boundary was used to repre- Boundary Conditions sent seepage from the West Side Canal; the lux rate The boundary conditions of a numerical was based on the average rate obtained from two groundwater model are assigned to all of the three- 2010 canal seepage runs (Abdo and others, 2013). dimensional boundary surfaces of the aquifer Because the Beaverhead River both contributes (re- system and to internal sources and sinks (ASTM, charges) water and drains (discharges) water from 2004). Boundary conditions represent the sources the aquifer system, the MODFLOW River Package of recharge and discharge to the groundwater low was used; this package allows groundwater low to system, and/or the hydraulic head at the edges of enter as well as exit the model. The larger sloughs the modeled domain. within the study area are believed to only drain The boundary conditions for the Lower Beaver- water from the aquifer system, and were simulated head model follow those discussed in the Hydrolog- as drains. It should be noted that the River Package ic Boundaries, Sources and Sinks, and Groundwater and Drain Package of MODFLOW do not calculate Budget sections of this report. They can be grouped stream discharge; the modules are used to calculate 27 Butler and Abdo, 2013
Figure 22. Boundary conditions in the groundwater fl ow model included specifi ed fl uxes, head-dependent fl uxes (such as the river and drains), and pumping wells. the gain or loss of groundwater to those features. tion ield recharge for much of the loodplain was Stream depletion from pumping is determined by eliminated during the calibration process, as dis- the change in such gains or losses. cussed in the Steady-State Calibration section. Irrigation Field Recharge Well Withdrawals A linear groundwater recharge rate was cal- Well withdrawals were estimated using the ap- culated for the irrigated areas within the model proach described in the Groundwater Budget sec- (equation 3, Groundwater Budget section). The tion and in appendix C. The wells were simulated approach calculated the cumulative recharge from using a speci ied- lux boundary (Well Package). precipitation and applied irrigation water, minus Only irrigation wells and PWS wells were simulated ET from crop consumption. Recharge rates were in the model. Domestic well withdrawals were not categorized into the three irrigation types within simulated due to their relatively minor effect on the the study area: lood, pivot, and sprinkler ( ig. 23). groundwater low system; they pump at relatively In non-irrigated areas, groundwater recharge from low rates and are spaced at a low density within precipitation was assumed to be negligible. Irriga- the study area. 28 Montana Bureau of Mines and Geology Open-File Report 638
Figure 23. Irrigation recharge was applied to the model in accordance with irrigation type. Irrigated fi elds from Montana Department of Revenue (2010). Some irrigated parcels in the fl oodplain were removed from the dataset and are not shown here (see text). CALIBRATION contoured using the default kriging method in Surfer Version 9 (Golden Software, Inc.; Abdo and Selection of Calibration Targets others, 2013). Because there was no major shift in Observed groundwater elevations were used as contours from month to month, the most robust calibration targets in the model. Groundwater level monthly dataset (July 2012) was used for the data were collected at selected wells monthly or steady-state calibration. bi-monthly during the project, beginning at various Prior to steady-state calibration, 16 of the 90 times during the fall of 2009 and winter of 2010, wells were removed from the July 2010 dataset and continuing until June 2011. About 90 well sites because they were located outside of the model provided adequate records of static water levels. domain. Another two wells were removed due to Without modi ication, the monthly 2010 data sets the non-static nature of their water levels. Dur- yielded quite similar potentiometric surfaces when ing calibration, three wells (207332, 152570, and 29 Butler and Abdo, 2013 237993) were deleted from the calibration data set included the residual mean, which should be close due to the inability of the model to calibrate suc- to zero in a well-calibrated model (i.e., the positive cessfully based on the polygon array. They are all and negative residuals balance one another); the located in the central West Bench area and are all mean of the absolute value of the residuals, which deeper wells (340 to 400 ft). They were removed is a measure of the average error in the model; and from the calibration data set due to their anoma- the root means square error, which is the square lously low water levels relative to the surrounding root of the average of the squared residuals. Last, wells (approximately 200 to 250 ft lower than the the residual standard deviation was divided by the surrounding water levels). The most likely cause overall range in observed values; this statistic is for the anomalies is that these wells are screened in a useful measure of how the spread of error com- the sedimentary bedrock aquifer discussed in the pares to the gradient across the entire model do- Geologic Framework section. This aquifer is found main. A value of less than 10% is generally accept- in portions of the central West Bench and appears able for calibration (Rumbaugh, 2007). Calibration to be disconnected from the shallower Tertiary data iles are provided with each set of groundwa- sediment and volcanic rock aquifers of the West ter model iles in appendix A. Note that the calibra- Bench. tion statistics were based only on data measured Also during calibration, control points (i.e., from the 69 actual observation wells and did not imaginary observations wells) were added to take the seven control points into account. better it heads to observed water levels. A total Steady-State Calibration of 14 points were added after preliminary PEST runs generated an unrealistic hydraulic gradient The steady-state version of the model simulates in certain areas due to a lack of observation data. average annual conditions for all components of The majority of control points were placed along recharge and discharge, and represents the sys- the model grid borders, where observation wells tem in equilibrium with a speci ied set of stresses. were lacking. As discussed in the Boundary Con- The steady-state model is useful for predicting the ditions section, the model borders were aligned ultimate impact to the groundwater low system with potentiometric surface contours (i.e., lines of from a new stress, such as a pumping well, and for uniform head); however, initial PEST K con igura- evaluating the overall groundwater budget. tions produced variable heads along these borders. The steady-state model was calibrated to ob- Consequently, the position of the potentiometric served values (i.e., calibration targets) through the surface was estimated in these areas and control use of manual trial and error as well as inverse points were entered to guide the model calcula- modeling, also known as automated parameter es- tions toward a realistic result. timation (PEST). Manual calibration was performed The removal of observation points and addition irst and involved adjusting input parameters of imaginary points resulted in a total of 83 cali- (hydraulic conductivity, recharge, and bed con- bration targets, 69 of which were real monitoring ductance values) until MODFLOW converged on a sites. The calibration tolerance interval was set at solution and produced reasonable head and water 5% of the range of observed groundwater eleva- budget values. Typically, only one parameter was tions within the model domain. This 5% criterion adjusted per model iteration in order to isolate its equated to a ±15 ft head residual; the residual is in luence relative to other input parameters. the difference between the modeled head value and PEST was used following the manual calibration the observed value. This value was selected based to further re ine the hydraulic conductivity assign- on the results of models of similar scale in Mon- ments in the model. In this approach, polygonal tana and Utah (Kauffman, 1999; Uthman and Beck, hydraulic conductivity zones were de ined, and 1998; Waren, 1998; Brooks and others, 2003). recharge values were held constant. Polygons were Along with the head residual criterion of ±15 ft, drawn to allow hydraulic conductivity to vary about error statistics were used during calibration. The the model. The polygons were based on known selected measure of error was the head residual, or suspected geologic boundaries, such as those which is the difference between the modeled head discussed in the Geologic Framework section. An value and the observed value. Error statistics effort was made to limit the number of K polygons 30 Montana Bureau of Mines and Geology Open-File Report 638 to reasonably sized zones that contained multiple and a more balanced steady-state budget. monitoring wells. The zones were assigned initial The model’s in low was further decreased by value ranges based on manual calibration results reducing the amount of aerial irrigation recharge and on the aquifer property estimates discussed in the loodplain. Aerial recharge to irrigated areas in the Aquifer Properties section (table 1). PEST mostly in the northern half of the loodplain was model runs were repeated and adjustments made removed, which was a mix of lood and pivot irriga- to the polygon con igurations and values in order to tion and covered approximately 48% of the total minimize the difference between computed heads loodplain area in the model. This reduction in re- and observed water levels. The generated hydraulic charge resulted in slightly lower head values in the conductivities and groundwater budget were evalu- northern loodplain and a more balanced average ated relative to the conceptual model, the results of annual budget. It may be that much of the excess aquifer tests, and the manually calculated ground- irrigation water applied in this portion of the study water budget to ensure they were reasonable. area returns to streams and ditches as surface- The modeled ranges of hydraulic conductivity water return lows. In this area of the model, the values ( ig. 24) were comparable to those estimated water table is shallow, so there is limited space in during the ield investigation (table 3), with the the aquifer to accommodate recharge. Excess water exception of the values assigned to the Tertiary may tend to be captured by the many smaller drain sediment aquifer underlying the alluvium in the ditches in the loodplain that were not modeled. loodplain. The modeled loodplain K values were This drain capture would explain why the model higher than aquifer test estimates, which was likely could not accommodate the excess waters without a result of bulk properties of the lower alluvium the addition of these ditches, for which there was and the upper Tertiary aquifer. Considering the limited information available. high value of the alluvial K and the limit of out low The third change aimed at lowering the north- at the Beaverhead Rock pinch point, the parameter ern loodplain heads involved the modeled ground- estimation process yielded a loodplain Tertiary water out low at the downstream end of the valley sediment aquifer K value that was higher than the (Beaverhead Rock). The layer 1 out low was in- East and West Bench K values. creased to re lect the higher alluvial K in the model The steady-state model was able to replicate (1,800 ft/day at Beaverhead Rock vs. the ield-esti- groundwater levels within a reasonable error ( ig. mated range of 600 to 1,233 ft/day). Furthermore, 25). Fifty-nine of the 69 computed heads were a low groundwater out low rate (480 acre-ft/yr) within the calibration target criterion of ±15 ft was set in layer 2; the initial conceptual model of (green bars in ig. 25). The 10 heads that did not the system assumed that out low is minimal be- meet this criterion had head values that ranged neath the alluvium at Beaverhead Rock. from 15.9 to 25.5 ft above observed levels (aver- The residuals for the calibrated model are age of 18.3 ft), and all fell within the north-central shown in igure 26. This graph displays the differ- loodplain area (yellow bars in ig. 25). Thus, this ence between observed water levels and the mod- area was the highest source of error in the model. el’s head results. The erroneously high heads in the The error most likely resulted from an imbalance of north-central loodplain were the points that plot water entering and exiting this portion of the lood- above the “perfect match” line. The error statistics plain. That is, more water entered than exited the shown in the igure are reasonable for a low model model area, which in turn resulted in heads higher of this scale; this is best indicated by the residual than the observed water levels. standard deviation over the range in head, which Though these loodplain heads did not meet the was 2.4%. This result met the calibration criterion calibration criterion of 15 ft, they did improve con- of 10% (Rumbaugh, 2007). siderably through certain calibration efforts. Spe- Stream base low was also used as a calibration ci ically, three main changes were made. First, the measure, because predicting changes in base low East Bench Canal seepage rate was reduced from was a primary modeling objective. Base low is the 2.8 cfs/mile (the upper range of the ield range esti- component of stream low supplied by groundwa- mate of 0.8 to 3.3 cfs/mile; Abdo and others, 2013) ter (Fetter, 2001). Calibration efforts were focused to 1.0 cfs/mile, which resulted in lower head values on the sloughs used in the pumping scenarios, 31 Butler and Abdo, 2013
Figure 24. The distribution of hydraulic conductivity values in the model was based on manual trial and error and on PEST.
32 Montana Bureau of Mines and Geology Open-File Report 638
Figure 25. Potentiometric surface of layer 1, generated in the steady-state model. The green bars represent cali- bration targets that were within ±15 ft; yellow bars represent targets that were were between 16 and 25.5 ft. namely Black, Willard, and Albers Sloughs. For each Transient Calibration slough, low measurements from an upstream and The transient version of the model used the downstream location measured during the non- aquifer properties from the steady-state model and irrigation season were used to determine a gain in added the element of time. The transient model low per mile. This base low amount was projected was used to simulate time-dependent stresses, over the stream distance featured in the model such as seasonal irrigation activities. Irrigation and used as a calibration target for each slough. A ield recharge, canal leakage, and well withdrawals difference of 15% or less between this value and were modi ied from the steady-state version of the the model results was considered reasonable. In model. Canal and irrigation recharge were applied the steady-state model, the differences between the from April 15 to October 15 to approximate the ir- ield and model results ranged from 9% to 14%, rigation season; well withdrawal rates were deter- which met this calibration criterion (table 5). mined based on available pumping data (appendix 33 Butler and Abdo, 2013
Figure 26. The residuals for the calibrated model display the difference between the observed water levels and those calculated by the model. The water levels that plot above the “perfect match” line were the high modeled levels that fell primarily in the fl oodplain. C). All other speci ied luxes, such as groundwater and relatively large seasonal luctuation was ob- in low, were held at constant steady-state values served, which was in the irrigated portions of the throughout the year, because the available data Tertiary sediments and alluvium. The hydraulic did not indicate a marked seasonal change in their conductivity values generated by calibration of rates. the steady-state model were used in the transient The model was calibrated to monthly data from modeling effort; minor adjustments were made to January 2010 through December 2010. This period some K values, after which the steady-state model of record was modeled as monthly stress periods, was rerun and re-calibrated. each of which had one time step (table 6). The Table 3 includes ield estimated storativity val- model was calibrated using the same observation ues and those assigned to the model. The bedrock wells (i.e., calibration targets) as in the steady-state units were assigned values ranging from 0.01 to calibration; however, the number of available water 0.05. Because bedrock observation well data were levels varied slightly from month to month, and a limited and showed little to no seasonal change different set of water-level values was input each through the period of record, calibration efforts month. The calibration tolerance interval was the were not focused on the bedrock system. Manual same as that of the steady-state Table 5. Baseflow comparison between field estimates and model results. model (±15 ft). Transient model calibration was conducted Baseflow: Baseflow: primarily by adjusting S values Flow Average Field Average until observed transient water- Measurement Measurement Model Result Percent level changes were reasonably Stream Events (cfs) (cfs) Difference replicated by the model. Efforts Willard Slough 4 2.0 2.3 14% were focused where a distinct Albers Slough 2 12.1 13.3 10% Black Slough 9 3.4 3.7 9% 34 Montana Bureau of Mines and Geology Open-File Report 638 7DEOH6WUHVVSHULRGVDQGWLPHVWHSVLQWKHWUDQVLHQW calibration resulted in storativity values ranging PRGHO from 0.05 to 0.08 in the portion of the model rep- 6WUHVVSHULRG 1RRIWLPH resenting the Tertiary sediment aquifer, and 0.15 'DWH7LPH OHQJWK GD\V VWHSV in the portions representing the alluvial and West $0 Bench volcanic rock aquifers. These values approxi- $0 mated the observed water-level changes caused by $0 seasonal changes in irrigation recharge. $0 Hydrographs comparing modeled to observed water levels are shown in igure 27; the well loca- $0 tions are shown in igure 25. The hydrographs illus- $0 trate that the modeled heads are in general agree- $0 ment with the observations, with the exception of $0 the northern loodplain area ( ig. 27, well 242404). $0 Like the steady-state model, this area exhibited $0 heads higher than the observed levels, though the $0 annual pattern in head was comparable. $0 Table 7 provides the water budget results in both the steady-state and transient versions of
Figure 27. The observed heads and computed heads were in general agreement. Well 242404, however, is an example of a northern fl oodplain well that showed considerable difference between the observed and computed water levels. 35 Butler and Abdo, 2013 the model. Comparison of the two model budgets model run while all other parameters remained at reveals that their total inputs and outputs are their calibrated values. The analysis was conducted similar; likewise, most of their individual budget on the steady-state model rather than the transient components are similar. The only major differences model due the relatively brief data record (1 year). are a lack of storage in the steady-state budget and R and K were adjusted four times (by +25%, -25%, a large discrepancy in base low. Due to the nature +50% and -50%), which resulted in a total of eight of a steady-state simulation, changes in storage simulations. The root mean square (RMS) error cannot be modeled; consequently, water that would was the calibration criterion used to evaluate the low into storage in a transient simulation must model’s sensitivity to each parameter change. be represented as another out low component in a Results showed the model output to be most steady-state simulation. In this particular steady- sensitive to the changes in K values, most notably state simulation, storage out low took the form of decreases in K. A 50% decrease in K values pro- base low to the Beaverhead River and its tributar- duced the largest RMS error ( ig. 28), which indi- ies; base low in the steady-state simulation was cates that the model is more sensitive to decreases 7,731 acre-ft/yr greater than that of the transient in K than in R. In contrast, increases in K and R simulation (table 7). The storage out low compo- generated comparable RMS errors. nent in the transient simulation agreed with the The effect on the spatial distribution of head 2010 water budget analysis (Water Budget section). residuals was also examined. A spatial sensitivity Model budget results re lect the fact that irriga- analysis helps to identify areas where con idence in tion (i.e., canal leakage and applied irrigation wa- parameter estimates is most important and, con- ter) was the predominant source of recharge to the versely, where accurate estimates are less impor- model, demonstrating its considerable in luence tant. Thus, the analysis can identify areas where on groundwater low within the Lower Beaverhead future data collection efforts should be focused. study area. The primary means of groundwater dis- This analysis indicated a considerably higher error charge from the model was base low to the Beaver- in the East and West Benches relative to the lood- head River and its tributaries. plain ( igs. 29, 30); this higher sensitivity is due in Sensitivity Analysis part to the relatively low K values in the benches. These results suggest that efforts to estimate K and A sensitivity analysis was performed to assess R values should be more focused on the Tertiary the sensitivity of the calibrated model to uncertain- sediment and volcanic rock aquifers rather than the ty in the aquifer parameter estimates. During the alluvial aquifer. Such data collection could also help analysis, recharge (R) and K were adjusted system- to independently estimate the parameters, thus atically, such that one parameter was changed per improving the model it Table 7. Water-budget values in the steady-state and transient versions of the model. and re ining the concep- Steady-State 2010 Transient tual understanding of Budget Values Budget Values the groundwater low (ac-ft/yr) (ac-ft/yr) system. Inputs Model Veriϐication East Bench inflow 3,161 3,161 Floodplain groundwater inflow 2,765 2,765 Model veri ication West Bench inflow 142 142 is a process in which West Side Canal leakage 5,921 5,925 calibrated parameters East Bench Canal leakage 4,641 4,594 and stresses are used to Infiltration from applied irrigation water 35,994 36,044 reproduce a second set TOTAL INPUT 52,624 52,631 of ield data; the process Outputs is intended to provide Irrigation + PWS well withdrawals 8,153 8,160 greater con idence in River (net gain at Beaverhead Rock, including sloughs) 41,631 33,918 the model (Anderson Outflow at Beaverhead Rock pinch point 2,840 2,840 and Woessner, 2002). Storage -- 7,705 Field data within the TOTAL OUTPUT 52,624 52,622 36 Montana Bureau of Mines and Geology Open-File Report 638
Figure 28. RMS error caused by changing K and recharge values by various percentages. A 50% decrease in K values produced the largest RMS error, which indicates that the model is more sensitive to decreases in K than in R. Errors generated by increases in K and R were comparable.
Lower Beaverhead study area were limited out- were intended to predict groundwater levels and side of the 2010–2011 study period. The available groundwater low under the hypothetical modeled water-level records did not constitute a second set conditions. In reality, future conditions will inevi- of ield data, so model veri ication was not possible. tably differ from the modeled conditions due to However, the model could be veri ied in the future changes in climate, land use, and other factors. if monitoring of water levels, irrigation practices, Scenarios were run to predict changes in the and climate continues in the study area. Veri ication groundwater low system, with particular atten- would test whether the model output (i.e., head tion given to the effects of varying pumping well values) agrees reasonably well with the observed rates and locations within the aquifer of the north- water levels in area wells. ern West Bench. The simulations were created by extending the model stress periods into the future PREDICTIVE SIMULATIONS and specifying the stresses to be tested. Speci ically, Prediction is one of the three main applications 240 1-month stress periods were set up to oper- of modeling (Anderson and Woessner, 2002). In the ate the model for 20 years. Pumping stresses were Lower Beaverhead groundwater investigation, the applied for 2 months of the irrigation season each modeling purpose was to predict the consequences year. Pumping rates and durations were based on of a proposed action; namely, the consequences groundwater usage estimates of existing irrigation of pumping new high-capacity irrigation wells on wells in the area. the West Bench. However, the predictive modeling The predictive simulations also included canal scenarios described below were not attempts to seepage scenarios, in which groundwater recharge predict the real future; in other words, a given 20- was increased by extending the canal low into year simulation was not intended to represent the the pre- and/or post-irrigation season. The canal groundwater levels and groundwater low that will seepage scenarios were compared to baseline and occur in the next 20 years. Rather, the scenarios pumping scenario results to examine how the ad- 37 Butler and Abdo, 2013
Figure 29. Head residuals resulting from a 50% decrease in K values indicated a considerably higher error in the East and West Benches relative to the fl oodplain; this higher sensitivity is due in part to the relatively low K values in the benches. those of the baseline scenario in order to predict ditional groundwater recharge offset stream deple- stream depletion and groundwater drawdown. tion. Stream depletion was examined in the Bea- Scenario 2 verhead River, Black Slough, Willard Slough, and Albers Slough. Scenario 2 involved well A ( ig. 31), a hypotheti- Table 8 provides a description of the modeling cal irrigation well completed in the volcanic rock scenarios, and igure 31 shows the pumping well aquifer, pumping at 1,500 gpm for 2 months of locations. The results from pumping scenarios 2 the irrigation season. The simulation resulted in a through 4 and all three canal seepage scenarios are maximum drawdown of 6.8 ft at well A, which oc- presented. The baseline scenario (scenario 1) is curred in August of the inal year of pumping (year brie ly presented for comparison purposes. 20). The increase in drawdown from year to year decreased over time, with an increase of 0.04 ft Scenario 1 between years 19 and 20. The irst scenario was the baseline scenario and The highest depletion rates also occurred in the featured only the pumping wells from the transient inal year of pumping ( ig. 32). Because the model model, which pumped seasonally at their assigned simulated the amount of groundwater lowing into 2010 rates throughout the simulation. The results each stream, stream depletion was estimated by of each subsequent scenario were compared to subtracting the base low during the pumping sce- 38 Montana Bureau of Mines and Geology Open-File Report 638
Figure 30. Head residuals resulting from a 50% increase in recharge values indicated higher errors in the East and West Benches than in the fl oodplain. A 50% increase in K values generated a similar degree and spatial distribution of model sensitivity. nario from that of the baseline scenario. The Beaverhead River was divided into nine When comparing the depletion in the three segments of equal 2-mile lengths to examine the sloughs, a distance–magnitude relationship can distribution of depletion along a stream reach. The be observed. The highest depletion was in Black depletion for each segment was extracted from the Slough, which was closest to the pumping well, model output and then calculated as a percentage followed by Willard and then Albers Slough. The of the overall depletion of the entire river reach. Beaverhead River had a depletion rate between The percentage for each segment is shown in igure that of Black Slough and Willard Slough. Figure 32 33. The greatest depletion occurred in the middle also illustrates a distance–time depletion relation- segments of the river, which were slightly upstream ship. For instance, a relatively immediate response of the segment closest to the pumping well. was observed in the closest stream to well A (Black Scenario 3 Slough) during and soon after the pumping period. In contrast, depletion in the farthest stream (the Scenario 3 featured pumping well A (scenario Beaverhead River) gradually increased during the 2) along with a second well (well B) completed pumping period and did not reach its maximum de- in the volcanic rock aquifer ( ig. 31). Both wells A pletion rate until 2 months after pumping ceased. and B were pumped at 1,500 gpm for 2 months of the irrigation season. The simulation results were
39 Butler and Abdo, 2013
7DEOH6XPPDU\RISXPSLQJDQGFDQDOVHHSDJHVFHQDULRV
3XPSLQJ6FHQDULRV 3XPSLQJ5DWH 3XPSLQJ 1RRI1HZ RI:HOO V 'XUDWLRQ 3XPSLQJ 6FHQDULR 'HVFULSWLRQ :HOOV :HOO1DPH V JSP 0RQWKV
FRQGLWLRQV 1$ 1$ 1$ 1$ ZHOOLQ YROFDQLFDTXLIHU :HOO$ -XQH$XJXVW ZHOOVLQYROFDQLF DTXLIHU :HOO$% ZHOO -XQH$XJXVW ZHOOVLQ7HUWLDU\ VHGLPHQWDTXLIHU :HOO&' ZHOO -XQH$XJXVW &DQDO6HHSDJH6FHQDULRV &DQDO 3XPSLQJ5DWH 3XPSLQJ 3UH3RVW RI:HOO V 'XUDWLRQ 3XPSLQJ 6FHQDULR 'HVFULSWLRQ 6HDVRQ3HULRG :HOO1DPH V JSP 0RQWKV