USGS Science⌕ Search

Geology topics

Marshall W. Gannett

Publications and source records attributed to Marshall W. Gannett.

At least 19 recordsLinked to original sources
USGS2021

Groundwater, biodiversity, and the role of flow system scale

Groundwater-dependent ecosystems and species (GDEs) are found throughout watersheds at locations of groundwater discharge, yet not all GDEs are the same, nor are the groundwater systems supporting them. Groundwater moves along a variety of flow paths of different lengths and with different contributing areas, ranging from shorter local flow paths with low discharge and large seasonal variability to streams, springs and wetlands to longer regional flow paths with potentially larger discharge and low seasonal variability, commonly at low basin elevations. How does this variation in physical hydrology affect the type and distribution of GDEs? Using data on hypsographic position, groundwater-dependent species distributions, groundwater pumping and streamflow from Oregon, USA, we provide a conceptual model and initial supporting evidence demonstrating that spatial variation in groundwater flow path scales, illustrated using basin hypsography, is a driver of non-random distribution of GDEs across watersheds. Further, we posit that the spatial variation in primary stressors to groundwater (e.g. pumping and climate change) will differentially affect GDEs depending on their hypsographic position. Furthermore, because of their use for irrigation and municipal water supply, regional groundwater systems and associated species are more likely to be studied and receive regulatory protection. Our initial data point to a disproportionate focus on larger discharge, lower elevation GDEs, which leads to a bias in our understanding of the full suite of biodiversity associated with groundwater discharge as well as their stressors and potential mechanisms for protection.

Ecohydrology↗
USGS2019

Landscape controls on the distribution and ecohydrology of central Oregon springs

Small springs in semiarid landscapes are essential for maintaining aquatic biodiversity and supporting livestock grazing operations. However, little is known about controls on the distribution and physical characteristics of small springs, the aquatic species they support, or their sensitivity to disturbance. We address this information gap in the Crooked River subbasin, a tributary of the Deschutes River in Oregon. We conducted spatial analyses on 2,519 mapped springs to investigate the influence of landscape controls (precipitation and bedrock permeability) on spring density in the Crooked River subbasin and the adjacent Upper Deschutes subbasin. Spring density was highest in areas of low bedrock permeability ( P < 0.0001) and high annual precipitation ( P < 0.0001). We suggest that the high density of small springs on low‐permeability bedrock indicates that these springs generally have short, shallow flow paths and thus may be susceptible to forecasted climate changes. A survey of 137 springs in the Crooked River subbasin revealed the hydrogeologic setting affects spring discharge type ( P = 0.017), temperature ( P = 0.011), and pH ( P = 0.026). We found a high frequency of anthropogenic impacts on springs: 95% of diffuse‐discharge springs and 79% of discrete‐discharge springs were disturbed by livestock grazing. Species inventories at 10 of the most intact surveyed springs confirm that small springs are biologically diverse, with 151 total species of plants and 135 total taxa of macroinvertebrates. Springs in the Crooked River subbasin are ecologically important habitats but require careful management to protect against livestock disturbance and development.

Oregon↗
USGS2017

Simulation of groundwater and surface-water flow in the upper Deschutes Basin, Oregon

This report describes a hydrologic model for the upper Deschutes Basin in central Oregon developed using the U.S. Geological Survey (USGS) integrated Groundwater and Surface-Water Flow model (GSFLOW). The upper Deschutes Basin, which drains much of the eastern side of the Cascade Range in Oregon, is underlain by large areas of permeable volcanic rock. That permeability, in combination with the large annual precipitation at high elevations, results in a substantial regional aquifer system and a stream system that is heavily groundwater dominated. The upper Deschutes Basin is also an area of expanding population and increasing water demand for public supply and agriculture. Surface water was largely developed for agricultural use by the mid-20th century, and is closed to additional appropriations. Consequently, water users look to groundwater to satisfy the growing demand. The well‑documented connection between groundwater and the stream system, and the institutional and legal restrictions on streamflow depletion by wells, resulted in the Oregon Water Resources Department (OWRD) instituting a process whereby additional groundwater pumping can be permitted only if the effects to streams are mitigated, for example, by reducing permitted surface-water diversions. Implementing such a program requires understanding of the spatial and temporal distribution of effects to streams from groundwater pumping. A groundwater model developed in the early 2000s by the USGS and OWRD has been used to provide insights into the distribution of streamflow depletion by wells, but lacks spatial resolution in sensitive headwaters and spring areas. The integrated model developed for this project, based largely on the earlier model, has a much finer grid spacing allowing resolution of sensitive headwater streams and important spring areas, and simulates a more complete set of surface processes as well as runoff and groundwater flow. In addition, the integrated model includes improved representation of subsurface geology and explicitly simulates the effects of hydrologically important fault zones not included in the previous model. The upper Deschutes Basin GSFLOW model was calibrated using an iterative trial and error approach using measured water-level elevations (water levels) from 800 wells, 144 of which have time series of 10 or more measurements. Streamflow was calibrated using data from 21 gage locations. At 14 locations where measured flows are heavily influenced by reservoir operations and irrigation diversions, so called “ naturalized ” flows, with the effects of reservoirs and diversion removed, developed by the Bureau of Reclamation, were used for calibration. Surface energy and moisture processes such as solar radiation, snow accumulation and melting, and evapotranspiration were calibrated using national datasets as well as data from long-term measurement sites in the basin. The calibrated Deschutes GSFLOW model requires daily precipitation, minimum and maximum air temperature data, and monthly data describing groundwater pumping and artificial recharge from leaking irrigation canals (which are a significant source of groundwater recharge). The calibrated model simulates the geographic distribution of hydraulic head over the 5,000 ft range measured in the basin, with a median absolute residual of about 53 ft. Temporal variations in head resulting from climate cycles, pumping, and canal leakage are well simulated over the model area. Simulated daily streamflow matches gaged flows or calculated naturalized flows for streams including the Crooked and Metolius Rivers, and lower parts of the mainstem Deschutes River. Seasonal patterns of runoff are less well fit in some upper basin streams. Annual water balances of streamflow are good over most of the model domain. Model fit and overall capabilities are appropriate for the objectives of the project. The integrated model results confirm findings from other studies and models indicating that most streamflow in the upper Deschutes Basin comes directly from groundwater discharge. The integrated model provides additional insights about the components of streamflow including direct groundwater discharge to streams, interflow, groundwater discharge to the land surface (Dunnian flow), and direct runoff (Hortonian flow). The new model provides improved capability for exploring the timing and distribution of streamflow capture by wells, and the hydrologic response to changes in other external stresses such as canal operation, irrigation, and drought. Because the model uses basic meteorological data as the primary input; and simulates surface energy and moisture balances, groundwater recharge and flow, and all components of streamflow; it is well suited for exploring the hydrologic response to climate change, although no such simulations are included in this report. The model was developed as a tool for future application; however, example simulations are provided in this report. In the example simulations, the model is used to explore the influence of well location and geologic structure on stream capture by pumping wells. Wells were simulated at three locations within a 12-mi area close to known groundwater discharge areas and crossed by a regional fault zone. Simulations indicate that the magnitude and timing of stream capture from pumping is largely controlled by the geographic location of the wells, but that faults can have a large influence on the propagation of pumping stresses.

Oregon↗
USGS2017

Geothermal implications of a refined composition-age geologic map for the volcanic terrains of southeast Oregon, northeast California, and southwest Idaho, USA

Sufficient temperatures to generate steam likely exist under most of the dominantly volcanic terrains of southeast Oregon, northeast California, and southeast Idaho, USA, but finding sufficient permeability to allow efficient advective heat exchange is an outstanding challenge. A new thematic interpretation of existing state-level geologic maps provides an updated and refined distribution of the composition and age of geologic units for the purposes of assessing the implications for measurement and development of geothermal resources. This interpretation has been developed to better understand geothermal and hydrologic resources of the region. Comparison of the new geologic categories with available hydrologic data shows that younger volcanogenic terrains tend to have higher primary permeability than older terrains. Decrease in primary permeability with age is attributable to weathering and hydrothermal alteration of volcanogenic deposits to pore-filling clays and deposition of secondary deposits (e.g., zeolites). Spring density as a function of geology and precipitation can be used to infer groundwater flow path length within the upper aquifers. Beneath the upper aquifers, we postulate that, due to accelerated hydrothermal alteration at temperatures ~>30 °C, primary permeability at depths of geothermal interest will be limited, and that secondary permeability is a more viable target for hydrothermal fluid withdrawal. Because open fractures resulting from tensile stresses will affect all geologic layers, regions with a significant amount of groundwater flow through shallow, structurally controlled secondary permeability may overlay zones of deep secondary permeability. Regardless of whether the shallow permeability is connected with the deep permeability, shallow groundwater flow can mask the presence of deep hydrothermal flow, resulting in blind geothermal systems. Ideally, hydraulic connectivity between shallow and deep secondary permeability is limited, so that shallow groundwater does not cool potential geothermal reservoirs.

Conference Paper↗
USGS2015

Groundwater levels, trends, and relations to pumping in the Bureau of Reclamation Klamath Project, Oregon and California

The use of groundwater to supplement surface-water supplies for the Bureau of Reclamation Klamath Project in the upper Klamath Basin of Oregon and California markedly increased between 2000 and 2014. Pre-2001 groundwater pumping in the area where most of this increase occurred is estimated to have been about 28,600 acre-feet per year. Subsequent supplemental pumping rates have been as high as 128,740 acre-feet per year. During this period of increased pumping, groundwater levels in and around the Bureau of Reclamation Klamath Project have declined by about 20-25 feet. Water-level declines are largely due to the increased supplemental pumping, but other factors include increased pumping adjacent to the Klamath Project and drying climate conditions. This report summarizes the distribution and magnitude of supplemental groundwater pumping and groundwater-level declines, and characterizes the relation between the stress and response in subareas of the Klamath Project to aid decision makers in developing groundwater-management strategies.

California, Oregon↗
USGS2015

Effects of groundwater pumping on agricultural drains in the Tule Lake subbasin, Oregon and California

Since 2001, irrigators in the upper Klamath Basin have increasingly turned to groundwater to compensate for reductions in surface-water allocation caused by shifts from irrigation use to instream flows for Endangered Species Act listed fishes. The largest increase in groundwater pumping has been in and around the Bureau of Reclamation&rsquo;s Klamath Irrigation Project, which includes the Tule Lake subbasin in the southern part of the upper Klamath Basin. Agricultural drains on the Klamath Project are an important source of water for downstream irrigators and for the Tule Lake and Lower Klamath Lake National Wildlife Refuges. U.S. Geological Survey regional groundwater-flow model simulations and records of irrigation-return flow pumped from the Tule Lake subbasin into the adjacent Lower Klamath Lake subbasin have indicated that water-level declines from pumping may be causing decreased flow of shallow groundwater to agricultural drains. To better define the effect of increased pumping on drain flow and on the water balance of the groundwater system, the annual water volume pumped from drains in three subareas of the Tule Lake subbasin was estimated and a fine-grid, local groundwater model of the Tule Lake subbasin was constructed. Results of the agricultural-drain flow analysis indicate that groundwater discharge to drains has decreased such that flows in 2012 were approximately 32,400 acre-ft less than the 1997&ndash;2000 average flow. This decrease was concentrated in the northern and southeastern parts of the subbasin, which corresponds with the areas of greatest groundwater pumping. Model simulation results of the Tule Lake subbasin groundwater model indicate that increased supplemental pumping is the dominant stress to the groundwater system in the subbasin. Simulated supplemental pumping and decreased recharge from irrigation between 2000 and 2010 totaled 323,573 acre-ft, 234,800 acre-ft (73 percent) of which was from supplemental pumping. The response of the groundwater system to this change in stress included about 180,500 acre-ft (56 percent) of decreased groundwater discharge to drains and a 126,000 acre-ft (39 percent) reduction in aquifer storage. The remaining 5 percent came from reduced groundwater flow to other model boundaries, including the Lost River, the Tule Lake sumps, and interbasin flow.

California, Oregon↗
USGS2015

Geologic and geomorphic controls on the occurrence of fens in the Oregon Cascades and implications for vulnerability and conservation

Montane fens are biologically diverse peat-forming wetlands that develop at points of groundwater discharge. To protect these ecosystems, it is critical to understand their locations on the landscape and the hydrogeologic systems that support them. The upper Deschutes Basin has a groundwater flow system that supports baseflow in many rivers, but little is known about the wetland types and groundwater dependence of the thousands of wetlands within the watershed. In 292 randomly selected wetlands, we quantified landscape metrics thought useful for discriminating montane fens from non-peat-forming wetlands. We inspected these wetlands and classified 67 of them as fens. Of the landscape metrics, only geology reliably differentiated fens from other types of wetlands. Nearly all fens develop in low-permeability glacial till found at approximately 1400&ndash;1800 m in elevation, and are concentrated in areas mantled by pumice deposits that originated primarily from the eruption of Mt. Mazama approximately 7700 years BP. Stratigraphic and hydrologic factors indicate the fens are supplied by perched aquifers in glacial till, instead of the deeper regional aquifer system. Their hydrogeologic setting makes the fens highly vulnerable to expected changes to recharge associated with climate change, but not to groundwater pumping from the regional aquifer.

Oregon↗
USGS2014

Evaluation of alternative groundwater-management strategies for the Bureau of Reclamation Klamath Project, Oregon and California

The water resources of the upper Klamath Basin, in southern Oregon and northern California, are managed to achieve various complex and interconnected purposes. Since 2001, irrigators in the Bureau of Reclamation Klamath Irrigation Project (Project) have been required to limit surface-water diversions to protect habitat for endangered freshwater and anadromous fishes. The reductions in irrigation diversions have led to an increased demand for groundwater by Project irrigators, particularly in drought years. The potential effects of sustained pumping on groundwater and surface-water resources have caused concern among Federal and state agencies, Indian tribes, wildlife groups, and groundwater users. To aid in the development of a viable groundwater-management strategy for the Project, the U.S. Geological Survey, in collaboration with the Klamath Water and Power Agency and the Oregon Water Resources Department, developed a groundwater-management model that links groundwater simulation with techniques of constrained optimization. The overall goal of the groundwater-management model is to determine the patterns of groundwater pumping that, to the extent possible, meet the supplemental groundwater demands of the Project. To ensure that groundwater development does not adversely affect groundwater and surface-water resources, the groundwater-management model includes constraints to (1) limit the effects of groundwater withdrawal on groundwater discharge to streams and lakes that support critical habitat for fish listed under the Endangered Species Act, (2) ensure that drawdowns do not exceed limits allowed by Oregon water law, and (3) ensure that groundwater withdrawal does not adversely affect agricultural drain flows that supply a substantial portion of water for irrigators and wildlife refuges in downslope areas of the Project. Groundwater-management alternatives were tested and designed within the framework of the Klamath Basin Restoration Agreement (currently [2013] awaiting authorizing Federal legislation), which would establish a permanent limit on the amount of surface water that can be diverted annually to the Project. Groundwater-management scenarios were evaluated for the period 1970•2004; supplemental groundwater demand by the Project was estimated as the part of irrigation demand that would not have been satisfied by the surface-water diversion allowed under the Klamath Basin Restoration Agreement. Over the 35-year management period, 22 years have supplemental groundwater demand, which ranges from a few thousand acre-feet (acre-ft) to about 100,000 acre-ft in the driest years. The results of the groundwater-management model indicate that supplemental groundwater pumping by the Project can be managed to avoid adverse effects to groundwater discharge that supports critical aquatic habitat. The existing configuration of wells in the Project would be able to meet groundwater-pumping goals in 14 of the 22 years with supplemental groundwater demand; however, substantial irrigation shortages can be expected during drought periods when the demand for supplemental groundwater is highest. The maximum irrigation-season withdrawal calculated by the groundwater-management model is about 60,000 acre-ft, the average withdrawal in drought years is about 54,000 acre-ft, and the amount of unmet groundwater demand reaches a maximum of about 45,000 acre-ft. A comparison of optimized groundwater withdrawals by geographic region shows that the highest annual withdrawals are associated with wells in the Tule Lake and Klamath Valley regions of the Project. The patterns of groundwater withdrawal also show that a substantial amount of the available pumping capacity is unused due to the restrictions imposed by drawdown constraints. Subsequent model applications were used to evaluate the sensitivity of optimization results to various factors. A sensitivity analysis quantified the changes in optimized groundwater withdrawals that result from changes in drawdown-constraint limits. The analysis showed the potential for substantial increases in withdrawals of groundwater with less restrictive drawdown limits at drawdown-control sites in the California part of the model. Systematic variation of the drains-constraint limit yielded a trade-off curve between optimized groundwater withdrawals and the allowable reduction in groundwater discharge to the Project drain system. Additional model applications were used to assess the value of increasing the pumping capacity of the network of wells serving the Project, and the relation between reduced off-Project groundwater pumping and increased pumping by Project irrigators.

California;Oregon↗
USGS2013

Spatial variability of the response to climate change in regional groundwater systems -- examples from simulations in the Deschutes Basin, Oregon

We examine the spatial variability of the response of aquifer systems to climate change in and adjacent to the Cascade Range volcanic arc in the Deschutes Basin, Oregon using downscaled global climate model projections to drive surface hydrologic process and groundwater flow models. Projected warming over the 21st century is anticipated to shift the phase of precipitation toward more rain and less snow in mountainous areas in the Pacific Northwest, resulting in smaller winter snowpack and in a shift in the timing of runoff to earlier in the year. This will be accompanied by spatially variable changes in the timing of groundwater recharge. Analysis of historic climate and hydrologic data and modeling studies show that groundwater plays a key role in determining the response of stream systems to climate change. The spatial variability in the response of groundwater systems to climate change, particularly with regard to flow-system scale, however, has generally not been addressed in the literature. Here we simulate the hydrologic response to projected future climate to show that the response of groundwater systems can vary depending on the location and spatial scale of the flow systems and their aquifer characteristics. Mean annual recharge averaged over the basin does not change significantly between the 1980s and 2080s climate periods given the ensemble of global climate models and emission scenarios evaluated. There are, however, changes in the seasonality of groundwater recharge within the basin. Simulation results show that short-flow-path groundwater systems, such as those providing baseflow to many headwater streams, will likely have substantial changes in the timing of discharge in response changes in seasonality of recharge. Regional-scale aquifer systems with flow paths on the order of many tens of kilometers, in contrast, are much less affected by changes in seasonality of recharge. Flow systems at all spatial scales, however, are likely to reflect interannual changes in total recharge. These results provide insights into the possible impacts of climate change to other regional aquifer systems, and the streams they support, where discharge points represent a range of flow system scales.

Oregon↗
USGS2013

Evapotranspiration from marsh and open-water sites at Upper Klamath Lake, Oregon, 2008--2010

Water allocation in the Upper Klamath Basin has become difficult in recent years due to the increase in occurrence of drought coupled with continued high water demand. Upper Klamath Lake is a central component of water distribution, supplying water downstream to the Klamath River, supplying water for irrigation diversions, and providing habitat for various species within the lake and surrounding wetlands. Evapotranspiration (ET) is a major component of the hydrologic budget of the lake and wetlands, and yet estimates of ET have been elusive—quantified only as part of a lumped term including other substantial water-budget components. To improve understanding of ET losses from the lake and wetlands, measurements of ET were made from May 2008 through September 2010. The eddy-covariance method was used to monitor ET at two wetland sites continuously during this study period and the Bowen-ratio energy-balance method was used to monitor open-water lake evaporation at two sites during the warmer months of the 3 study years. Vegetation at one wetland site (the bulrush site) consists of a virtual monoculture of hardstem bulrush (formerly Scirpus acutus, now Schoenoplectus acutus), and at the other site (the mixed site) consists of a mix of about 70 percent bulrush, 15 percent cattail (Typha latifolia), and 15 percent wocus (Nuphar polysepalum). Measured ET at these two sites was very similar (means were ±2.5 percent) and mean wetland ET is computed as a 70 to 30 percent weighted average of the bulrush and mixed sites, respectively, based on community-type distribution estimated from satellite imagery. Biweekly means of wetland ET typically vary from maximum values of around 6 to 7 millimeters per day during midsummer, to minimum values of less than 1 mm/d during midwinter. This strong annual signal primarily reflects life-cycle changes in the wetland vegetation, and the annual variation of radiative input to the surface and resulting temperature. The perennial vegetation begins each growing season submerged, emerges from the dead litter mat around late May or early June, reaches a maximum height of about 2.2 meters (m) during summer, senesces in October, and subsequently lodges over, contributing to the dead litter mat from previous years. Hydroperiods last about 5 to 6 months, typically beginning in January or February and ending in July or August, and have a minor influence on the annual ET cycle. These hydroperiods result from lake levels that typically vary about 1.3 m, from around 0.6 to 0.9 m above the wetland surface, to around 0.4 to 0.7 m below the wetland surface. An estimate of 3-year annual wetland ET, made by substituting early- and late-season data measured during 2009 for the missing periods in early 2008 and late 2010, is 0.938 meter per year (m/yr). Daily values of alfalfa-based reference ET (ET r ) were retrieved from the Bureau of Reclamation AgriMet Web site (http://www.usbr.gov/pn/agrimet/index.html) and are aggregated into biweekly, annual, and 3-year values (for consistency, the 3-year values are also computed using substitute data from 2009 for early 2008 and late 2010). These ET r values are computed from weather data measured at the nearby Agency Lake weather station (AGKO), and are based on the assumption that the alfalfa crop is green and vigorous year-round. The 3-year value of ET r is 1.145 m/yr, about 22 percent greater than wetland ET. A comparison of 2008–2010 alfalfa and pasture growing season actual ET with wetland ET is made using data from the more distant Klamath Falls AgriMet weather station (KFLO) because actual alfalfa and pasture ET are not computed for the AGKO site. During the 190-day average alfalfa growing season, wetland ET (0.779 m) is about 7 percent less than alfalfa ET (0.838 m). During the 195-day average pasture growing season, wetland ET (0.789 m) is about 18 percent greater than pasture ET (0.671 m). Assuming alfalfa and pasture ET are equal to wetland ET during the non-growing season, annual estimates become 0.997 m, 0.938 m, and 0.820 m from alfalfa, wetland, and pasture, respectively. Wetland crop coefficients (K c =ET/ET r ) are computed at daily, biweekly, and annual time steps. Approximate formulas are given to estimate daily values of growing season Kc, thereby allowing computation of daily growing season ET using ET r from the AGKO weather station. Biweekly values of growing season Kc are computed from ensemble average values of ET and ET r during the 3 study period growing seasons, and a single, mean Kc is computed for the non-growing season. Together, these provide relatively accurate estimates of biweekly ET during the study (RMSE=0.396 and 0.347 mm/d, r2 = 62 and 0.971 at the bulrush and mixed sites, respectively). A fourth-order polynomial fit of the biweekly growing season values to day of year provides a more automated form of ET computation. Measured ET at the bulrush wetland site during the current study compares very closely with growing-season ET estimated during a study in 1997 at nearly the same location. During the earlier study, ET was measured four times, using eddy covariance for 1- to 2-day periods, and was estimated between measurement periods using a Penman-Monteith model, calibrated to the measurements. Differences between time series of ET from the two studies are similar to interannual differences within the current study. Compared to the 1997 study, the current study measured larger ET rates in early summer and smaller rates in late summer, resulting in very similar growing-season totals. A study conducted in 2000 estimated ET from nearby fallowed cropland, using the Bowen-ratio energy balance method supplemented with Priestley-Taylor and crop-coefficient ET modeling. Seasonal timing of ET from three different crop types varied considerably, but growing-season totals were remarkably similar, at 0.435 ± 0.009 m. Wetland ET measured during the current study, evaluated over the same growing season was 0.718 m, or about 65 percent greater than the fallowed cropland ET. Open-water evaporation from Upper Klamath Lake was measured at two locations during the warmer months of 2008–2010 using the Bowen-ratio energy balance method. Measured rates were in general agreement with those measured in 2003 using the same method. Open-water evaporation and wetland ET were nearly equal during late June through early August, when wetland vegetation was green and abundant. As expected, open-water evaporation consistently exceeded wetland ET during late summer, as wetland ET responded to vegetation senescence while open water evaporation responded to extra available energy in the form of heat previously stored in the lake. Overall, open-water evaporation was 20 percent greater than wetland ET during the same period.

Oregon↗
USGS2013

Analysis of 1997–2008 groundwater level changes in the upper Deschutes Basin, Central Oregon

Groundwater-level monitoring in the upper Deschutes Basin of central Oregon from 1997 to 2008 shows water-level declines in some places that are larger than might be expected from climate variations alone, raising questions regarding the influence of groundwater pumping, canal lining (which decreases recharge), and other human influences. Between the mid-1990s and mid-2000s, water levels in the central part of the basin near Redmond steadily declined as much as 14 feet. Water levels in the Cascade Range, in contrast, rose more than 20 feet from the mid-1990s to about 2000, and then declined into the mid-2000s, with little or no net change. An existing U.S. Geological Survey regional groundwater-flow model was used to gain insights into groundwater-level changes from 1997 to 2008, and to determine the relative influence of climate, groundwater pumping, and irrigation canal lining on observed water-level trends. To utilize the model, input datasets had to be extended to include post-1997 changes in groundwater pumping, changes in recharge from precipitation, irrigation canal leakage, and deep percolation of applied irrigation water (also known as on-farm loss). Mean annual groundwater recharge from precipitation during the 1999–2008 period was 25 percent less than during the 1979–88 period because of drying climate conditions. This decrease in groundwater recharge is consistent with measured decreases in streamflow and discharge to springs. For example, the mean annual discharge of Fall River, which is a spring-fed stream, decreased 12 percent between the 1979–88 and 1999–2008 periods. Between the mid-1990s and late 2000s, groundwater pumping for public-supply and irrigation uses increased from about 32,500 to 52,000 acre-feet per year, partially because of population growth. Between 1997 and 2008, the rate of recharge from leaking irrigation canals decreased by about 58,000 acre-feet per year as a result of lining and piping of canals. Decreases in recharge from on-farm losses over the past decade were relatively small, approaching an estimated 1,000 acre-feet per year by the late 2000s. All these changes in the hydrologic budget contributed to declines in groundwater levels. Groundwater flow model simulations indicate that climate variations have the largest influence on groundwater levels throughout the upper Deschutes Basin, and that impacts from pumping and canal lining also contribute but are largely restricted to the central part of the basin that extends north from near Benham Falls to Lower Bridge, and east from Sisters to the community of Powell Butte. Outside of this central area, the water-level response from changes in pumping and irrigation canal leakage cannot be discerned from the larger response to climate-driven changes in recharge. Within this central area, where measured water-level declines have generally ranged from about 5 to 14 feet since the mid-1990s, climate variations are still the dominant factor influencing groundwater levels, accounting for approximately 60–70 percent of the measured declines. Post-1994 increases in groundwater pumping account for about 20–30 percent of the measured declines in the central part of the basin, depending on location, and decreases in recharge due to canal lining account for about 10 percent of the measured declines. Decreases in recharge from on-farm losses were simulated, but the effects were negligible compared to climate influences, groundwater pumping, and the effects of canal lining and piping. Observation well data and model simulation results indicate that water levels in the Cascade Range rose and declined tens of feet in response to wet and dry climate cycles over the past two decades. Water levels in the central part of the basin, in contrast, steadily declined during the same period, with the rate of decline lessening during wet periods. This difference is because the water-level response from recharge is damped as water moves (diffuses) from the principal recharge area in the Cascade Range to discharge points along the main stems of the Deschutes, Crooked, and Metolius Rivers in the central part of the basin. Water levels in the central part of the basin respond more to multi-decadal climate trends than shorter term changes. Groundwater-flow simulations show that the effects from increased pumping and decreased irrigation canal leakage extend south into the Bend area. However, the only wells presently monitored in the Bend area are heavily influenced by the Deschutes River, which dampens any response of water levels to external stresses such as groundwater pumping, changes in canal leakage, or climate variations.

Oregon↗
USGS2012

Groundwater simulation and management models for the upper Klamath Basin, Oregon and California

The upper Klamath Basin encompasses about 8,000 square miles, extending from the Cascade Range east to the Basin and Range geologic province in south-central Oregon and northern California. The geography of the basin is dominated by forested volcanic uplands separated by broad interior basins. Most of the interior basins once held broad shallow lakes and extensive wetlands, but most of these areas have been drained or otherwise modified and are now cultivated. Major parts of the interior basins are managed as wildlife refuges, primarily for migratory waterfowl. The permeable volcanic bedrock of the upper Klamath Basin hosts a substantial regional groundwater system that provides much of the flow to major streams and lakes that, in turn, provide water for wildlife habitat and are the principal source of irrigation water for the basin's agricultural economy. Increased allocation of surface water for endangered species in the past decade has resulted in increased groundwater pumping and growing interest in the use of groundwater for irrigation. The potential effects of increased groundwater pumping on groundwater levels and discharge to springs and streams has caused concern among groundwater users, wildlife and Tribal interests, and State and Federal resource managers. To provide information on the potential impacts of increased groundwater development and to aid in the development of a groundwater management strategy, the U.S. Geological Survey, in collaboration with the Oregon Water Resources Department and the Bureau of Reclamation, has developed a groundwater model that can simulate the response of the hydrologic system to these new stresses. The groundwater model was developed using the U.S. Geological Survey MODFLOW finite-difference modeling code and calibrated using inverse methods to transient conditions from 1989 through 2004 with quarterly stress periods. Groundwater recharge and agricultural and municipal pumping are specified for each stress period. All major streams and most major tributaries for which a substantial part of the flow comes from groundwater discharge are included in the model. Groundwater discharge to agricultural drains, evapotranspiration from aquifers in areas of shallow groundwater, and groundwater flow to and from adjacent basins also are simulated in key areas. The model has the capability to calculate the effects of pumping and other external stresses on groundwater levels, discharge to streams, and other boundary fluxes, such as discharge to drains. Historical data indicate that the groundwater system in the upper Klamath Basin fluctuates in response to decadal climate cycles, with groundwater levels and spring flows rising and declining in response to wet and dry periods. Data also show that groundwater levels fluctuate seasonally and interannually in response to groundwater pumping. The most prominent response is to the marked increase in groundwater pumping starting in 2001. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. Pumping from wells close (within a few miles) to groundwater discharge features, such as springs, drains, and certain streams, can affect those features within weeks or months of the onset of pumping, and the impacts can be essentially fully manifested in several years. Simulations indicate that seasonal variations in pumping rates are buffered by the groundwater system, and peak impacts are closer to mean annual pumping rates than to instantaneous rates. Thus, pumping effects are, to a large degree, spread out over the entire year. When pumping locations are distant (more than several miles) from discharge features, the effects take many years or decades to fully impact those features, and much of the pumped water comes from groundwater storage over a broad geographic area even after two decades. Moreover, because the effects are spread out over a broad area, the impacts to individual features are much smaller than in the case of nearby pumping. Simulations show that the discharge features most affected by pumping in the area of the Bureau of Reclamation's Klamath Irrigation Project are agricultural drains, and impacts to other surface-water features are small in comparison. A groundwater management model was developed that uses techniques of constrained optimization along with the groundwater flow model to identify the optimal strategy to meet water user needs while not violating defined constraints on impacts to groundwater levels and streamflows. The coupled groundwater simulation-optimization models were formulated to help identify strategies to meet water demand in the upper Klamath Basin. The models maximize groundwater pumping while simultaneously keeping the detrimental impacts of pumping on groundwater levels and groundwater discharge within prescribed limits. Total groundwater withdrawals were calculated under alternative constraints for drawdown, reductions in groundwater discharge to surface water, and water demand to understand the potential benefits and limitations for groundwater development in the upper Klamath Basin. The simulation-optimization model for the upper Klamath Basin provides an improved understanding of how the groundwater and surface-water system responds to sustained groundwater pumping within the Bureau of Reclamation's Klamath Project. Optimization model results demonstrate that a certain amount of supplemental groundwater pumping can occur without exceeding defined limits on drawdown and stream capture. The results of the different applications of the model demonstrate the importance of identifying constraint limits in order to better define the amount and distribution of groundwater withdrawal that is sustainable.

Oregon;California↗
USGS2007

Ground-Water Hydrology of the Upper Klamath Basin, Oregon and California

The upper Klamath Basin spans the California-Oregon border from the flank of the Cascade Range eastward to the Basin and Range Province, and encompasses the Klamath River drainage basin above Iron Gate Dam. Most of the basin is semiarid, but the Cascade Range and uplands in the interior and eastern parts of the basin receive on average more than 30 inches of precipitation per year. The basin has several perennial streams with mean annual discharges of hundreds of cubic feet per second, and the Klamath River at Iron Gate Dam, which represents drainage from the entire upper basin, has a mean annual discharge of about 2,100 cubic feet per second. The basin once contained three large lakes: Upper and Lower Klamath Lakes and Tule Lake, each of which covered areas of 100 to 150 square miles, including extensive marginal wetlands. Lower Klamath Lake and Tule Lake have been mostly drained, and the former lake beds are now cultivated. Upper Klamath Lake remains, and is an important source of irrigation water. Much of the wetland surrounding Upper Klamath Lake has been diked and drained, although efforts are underway to restore large areas. Upper Klamath Lake and the remaining parts of Lower Klamath and Tule Lakes provide important wildlife habitat, and parts of each are included in the Klamath Basin National Wildlife Refuges Complex. The upper Klamath Basin has a substantial regional ground-water flow system. The late Tertiary to Quaternary volcanic rocks that underlie the region are generally permeable, with transmissivity estimates ranging from 1,000 to 100,000 feet squared per day, and compose a system of variously interconnected aquifers. Interbedded with the volcanic rocks are late Tertiary sedimentary rocks composed primarily of fine-grained lake sediments and basin-filling deposits. These sedimentary deposits have generally low permeability, are not good aquifers, and probably restrict ground-water movement in some areas. The regional ground-water system is underlain and bounded on the east and west by older Tertiary volcanic and sedimentary rocks that have generally low permeability. Eight regional-scale hydrogeologic units are defined in the upper Klamath Basin on the basis of surficial geology and subsurface data. Ground water flows from recharge areas in the Cascade Range and upland areas in the basin interior and eastern margins toward stream valleys and interior subbasins. Ground water discharge to streams throughout the basin, and most streams have some component of ground water (baseflow). Some streams, however, are predominantly ground-water fed and have relatively constant flows throughout the year. Large amounts of ground water discharges in the Wood River subbasin, the lower Williamson River area, and along the margin of the Cascade Range. Much of the inflow to Upper Klamath Lake can be attributed to ground-water discharge to streams and major spring complexes within a dozen or so miles from the lake. This large component of ground water buffers the lake somewhat from climate cycles. There are also ground-water discharge areas in the eastern parts of the basin, for example in the upper Williamson and Sprague River subbasins and in the Lost River subbasin at Bonanza Springs. Irrigated agriculture is an integral part of the economy of the upper Klamath Basin. Although estimates vary somewhat, roughly 500,000 acres are irrigated in the upper Klamath Basin, about 190,000 acres of which are part of the Bureau of Reclamation Klamath Project. Most of this land is irrigated with surface water. Ground water has been used for many decades to irrigate areas where surface water is not available, for example outside of irrigation districts and stream valleys. Ground water has also been used as a supplemental source of water in areas where surface water supplies are limited and during droughts. Ground water use for irrigation has increased in recent years due to drought and shifts in surface-water allocation from irrigati

Scientific Investigations Report↗
USGS2006

An evaluation and review of water-use estimates and flow data for the Lower Klamath and Tule Lake National Wildlife Refuges, Oregon and California

The Lower Klamath and Tule Lake National Wildlife Refuges, located in the upper Klamath Basin of Oregon and California, encompass approximately 46,700 and 39,100 acres, respectively. Demand for water in the semiarid upper Klamath Basin has increased in recent years, resulting in the need to better quantify water availability and use in the refuges. This report presents an evaluation of water-use estimates for both refuges derived on the basis of two approaches. One approach used evaporation and evapotranspiration estimates and the other used measured inflow and outflow data. The quality of the inflow and outflow data also was assessed. Annual water use in the refuges, using evapotranspiration estimates, was computed with the use of different rates for each of four land-use categories. Annual water-use rates for grain fields, seasonal wetlands, permanently flooded wetlands with emergent vegetation, and open-water bodies were 2.5, 2.9, 2.63, and 4.07 feet per year, respectively. Total water use was estimated as the sum of the products of each rate and the number of acres in its associated land-use category. Mean annual (2003-2005) water use for the Lower Klamath and Tule Lake refuges was approximately 124,000 and 95,900 acre-feet, respectively. To estimate water deliveries needed for each refuge, first, annual precipitation for 2003-2005 was subtracted from the annual water use for those years. Then, an adjusted total was obtained by adding 20 percent to the difference to account for salinity flushing. Resulting estimated mean annual adjusted needed water deliveries in 2003-2005 for the Lower Klamath and Tule Lake refuges were 107,000 and 82,800 acre-feet, respectively. Mean annual net inflow to the refuges for 2003-2005 was computed by subtracting estimated and measured surface-water outflows from inflows. Mean annual net inflow during the 3-year period for the Lower Klamath refuge, calculated for a subsection of the refuge, was approximately 73,700 acre-feet. The adjusted needed water delivery for this section of the refuge, calculated from evapotranspiration estimates, was approximately 77,600 acre-feet. For the Tule Lake refuge, mean annual net inflow during the 3-year period was approximately 76,100 acre-feet, which is comparable to the estimated annual needed water delivery for the refuge of 82,800 acre-feet. For 1962-2005, mean annual net inflow to the Lower Klamath refuge was approximately 49,800 acre-feet, about 23,900 acre-feet less than for 2003-2005. Although mean April-September net inflows for 1962-2005 and 2003-2005 have remained fairly constant, annual net inflow has increased for October-March, which accounts for the difference. Consistently higher autumn and winter flow deliveries since the mid-1980s reflect a significant change in refuge management. More sections of the refuge are currently managed as seasonal wetlands than were in the 1960s and 1970s. Flow records for the Ady Canal at State Line Road, Klamath Straits Drain at State Line Road, and D Pumping Plant were evaluated for their data quality. On the basis of USGS flow-record criteria, all three flow records were rated as 'poor.' By definition, 95 percent of the daily flows in a record having this rating could be in error by more than 15 percent.

Scientific Investigations Report↗
USGS2004

Simulation of regional ground-water flow in the Upper Deschutes Basin, Oregon

This report describes a numerical model that simulates regional ground-water flow in the upper Deschutes Basin of central Oregon. Ground water and surface water are intimately connected in the upper Deschutes Basin and most of the flow of the Deschutes River is supplied by ground water. Because of this connection, ground-water pumping and reduction of artificial recharge by lining leaking irrigation canals can reduce the amount of ground water discharging to streams and, consequently, streamflow. The model described in this report is intended to help water-management agencies and the public evaluate how the regional ground-water system and streamflow will respond to ground-water pumping, canal lining, drought, and other stresses. Ground-water flow is simulated in the model by the finite-difference method using MODFLOW and MODFLOWP. The finite-difference grid consists of 8 layers, 127 rows, and 87 columns. All major streams and most principal tributaries in the upper Deschutes Basin are included. Ground-water recharge from precipitation was estimated using a daily water-balance approach. Artificial recharge from leaking irrigation canals and on-farm losses was estimated from diversion and delivery records, seepage studies, and crop data. Ground-water pumpage for irrigation and public water supplies, and evapotranspiration are also included in the model. The model was calibrated to mean annual (1993-95) steady-state conditions using parameter-estimation techniques employing nonlinear regression. Fourteen hydraulic-conductivity parameters and two vertical conductance parameters were determined using nonlinear regression. Final parameter values are all within expected ranges. The general shape and slope of the simulated water-table surface and overall hydraulic-head distribution match the geometry determined from field measurements. The fitted standard deviation for hydraulic head is about 76 feet. The general magnitude and distribution of ground-water discharge to streams is also well simulated throughout the model. Ground-water discharge to streams in the area of the confluence of the Deschutes, Crooked, and Metolius Rivers is closely matched. The model was also calibrated to transient conditions from 1978 to 1997 using traditional trial-and-error methods. Climatic cycles during this period provided an excellent regional hydrologic signal for calibration. Climate-driven water-level fluctuations are simulated with reasonable accuracy over most of the model area. The timing and magnitude of simulated water-level fluctuations caused by annual pulses of recharge from precipitation match those observed reasonably well, given the limitations of the time discretization in the model. Water-level fluctuations caused by annual canal leakage are simulated very well over most of the area where such fluctuations occur. The transient model also simulates the volumetric distribution and temporal variations in ground-water discharge reasonably well. The match between simulated and measured volume of and variations in ground-water discharge is, however, somewhat dependent on geographic scale. The rates of and variations in ground-water discharge are matched best at regional scales. Example simulations were made to demonstrate the utility of the model for evaluating the effects of ground-water pumping or canal lining. Pumping simulations show that pumped water comes largely from aquifer storage when pumping begins, but as the water table stabilizes, the pumping increasingly diminishes the discharge to streams and, hence, streamflow. The time it takes for pumping to affect streamflow varies spatially depending, in general, on the location of pumping relative to the discharge areas. Canal-lining simulations show similar effects.

Water-Resources Investigations Report↗
USGS2002

Geologic framework of the regional ground-water flow system in the Upper Deschutes Basin, Oregon

Ground water is increasingly relied upon to satisfy the needs of a growing population in the upper Deschutes Basin, Oregon. Hydrogeologic studies are being undertaken to aid in management of the ground-water resource. An understanding of the geologic factors influencing ground-water flow is basic to those investigations. The geology of the area has a direct effect on the occurrence and movement of ground water. The permeability and storage properties of rock material are influenced by the proportion, size, and degree of interconnection of open spaces the rocks contain. These properties are the result of primary geologic processes such as volcanism and sedimentation, as well as subsequent processes such as faulting, weathering, or hydrothermal alteration. The geologic landscape in the study area evolved during about 30 million years of volcanic activity related to a north-south trending volcanic arc, the current manifestation of which are today&rsquo;s Cascade Range volcanoes. The oldest rock unit in the upper Deschutes Basin study area, the John Day Formation, is a sequence of upper Eocene to lower Miocene volcanic and sedimentary rocks. Weathering and alteration of the rocks has resulted in very low permeability; consequently, the unit forms the hydrologic basement for the regional ground- water flow system throughout much of the area. The Deschutes Formation and age-equivalent deposits that overlie the John Day Formation, in contrast, are highly permeable and are the most widely used ground-water-bearing units in the study area. The Deschutes Formation consists of a variety of volcanic and sedimentary deposits ranging in age from late Miocene to Pliocene (approximately 7.5 to 4.0 million years). Three distinct depositional environments previously described for the formation provide useful hydrogeologic subdivisions. The ancestral Deschutes River deposits and some units within the arc- adjacent alluvial-plain region are among the highest yielding units within the Deschutes Formation, with some wells producing up to a few thousand gallons per minute. Opal Springs basalt, Pelton basalt, and the rhyodacite dome complex near Steelhead Falls are particularly productive subunits within the Deschutes Formation and provide tens to hundreds of cubic feet per second of ground-water discharge to the Deschutes and Crooked Rivers, upstream of Round Butte Dam. Most ground-water recharge in the upper Deschutes Basin occurs in Quaternary deposits of the Cascade Range and Newberry Volcano. These deposits are highly permeable, and the fractured character of the lava flows facilitates rapid infiltration of precipitation and snowmelt, as well as movement of ground water to lower elevations. Additional recharge from canal leakage occurs along sections of unlined canals near Bend, constructed on lava flows from Newberry Volcano. Hydrothermal alteration and secondary mineralization at depth beneath the Cascade Range and Newberry Volcano has drastically reduced the permeability of the material in those regions, effectively restricting most ground water to the strata above the altered rocks. The top of the hydrothermally altered region is considered the base of the regional ground-water system beneath the Cascade Range and Newberry Volcano. Structural features influence ground-water flow within the upper Deschutes Basin mainly by juxtaposing materials with contrasting permeability. This juxtaposition can be caused by fault movement or by the influence of a fault on subsequent deposition. Several depositional centers have formed along the base of fault-line scarps or in grabens within the study area, and the infilling sedimentary deposits have permeability that differs from the surrounding rocks. The effects of faults on ground-water flow may be masked in some areas. For example, the water-table gradient changes slope in the vicinity of the Sisters fault zone, but the slope change also corresponds with a major precipitation gradient change; therefore, any influence of the fault zone is unclear. Geologic units in the Deschutes Basin were divided into several distinct hydrogeologic units. In some instances the units correspond to existing stratigraphic divisions. In other instances, hydrogeologic units correspond to different facies within a single stratigraphic unit or formation. The hydrogeologic units include Quaternary sediment, deposits of the Cascade Range and Newberry Volcano, four zones within the Deschutes Formation and age-equivalent rocks that roughly correspond with depositional environments, and pre-Deschutes-age strata.

Water-Resources Investigations Report↗