Air temperature and precipitation outputs from 10 CMIP6 GCMs were input to the Precipitation-Runoff Modeling System hydrologic model to evaluate water and energy responses in three headwater basins to projected climate change over the twenty-first century. The headwater basins (398–801 km 2 ) are located within the Mataura River basin in the South Island of New Zealand. CMIP6 datasets included two emission scenarios [shared socioeconomic pathways (SSPs) SSP2-4.5 and SSP5-8.5]. Half of the 10 GCMs selected in the study have equilibrium climate sensitivity (ECS) values above 4.5°C, which has been considered the upper end of equilibrium climate sensitivity. Modeling results included increased annual air temperature, evapotranspiration, and precipitation by the end of the twenty-first century for both SSP emissions scenarios, both high- and low-ECS GCMs, and all three headwater basins. Monthly precipitation and evapotranspiration totals also increased for all or most months. Monthly streamflow changes generally corresponded with monthly precipitation changes. Snowpack decreased significantly in depth and seasonal duration in all basins. However, streamflow increased for all SSP and ECS groups and basins because increased precipitation was consistently greater than increased evapotranspiration losses. Sources of uncertainty include the GCMs, climate sensitivity, downscaling, bias adjustment, emission scenarios, and the hydrologic model. Simulated hydrologic responses based on climate data from GCMs with ECS values of greater than 4.5°C could be less plausible since previous studies have suggested true ECS ranges from 1.5° to 4.5°C.

First posted September 1, 2021 For additional information, contact: Director,California Water Science CenterU.S. Geological Survey6000 J Street, Placer HallSacramento, California 95819 The U.S. Geological Survey, in cooperation with the California Department of Water Resources (DWR), has constructed a new spatially detailed Precipitation-Runoff Modeling System (PRMS) model for the Merced River Basin, California, which is a tributary of the San Joaquin River in California. Operated through an Object User Interface (OUI) with Ensemble Streamflow Prediction (ESP) and daily climate distribution preprocessing functionality, the model is calibrated primarily to simulate (and eventually, forecast) year-to-year variations of inflows to Lake McClure during the critical April–July snowmelt season. The model is intended to become part of a suite of methods used by DWR for estimating daily streamflow from the Merced River Basin, especially during the snowmelt season. This study describes the results of the application of an analysis tool that simulates responses to climate and land-use variations at a higher spatial resolution than previously available to DWR.A geographic information system was used to delineate the model domain, that is, areas draining to a single outlet at U.S. Geological Survey streamflow-gaging station 11270900, Merced River below Merced Falls Dam, near Snell, CA (also known as California Data Exchange Center station MRC), and subdrainage areas, including four draining to internal gages used as calibration targets. Using this delineation, three contiguous subbasins were recognized and, along with the model domain and nested calibration targets, are the simulation units evaluated in this report.An auto-calibration tool, LUCA (Let Us CAlibrate), was used for each calibration node, from headwaters to basin outlet, and then parameters were manually adjusted to complete the calibration. The main objective was to match April–July snowmelt seasonal discharge values of simulated streamflow to observed (measured or reconstructed) discharge values. Calibration or validation periods used site-specific streamflows—mostly from October 1, 1988, through September 30, 2013—but differed according to the period-of-record available for the measurements collected at internal gages or reconstructed flows for the single outlet.The accuracy of the Merced PRMS streamflow simulations varied seasonally, as compared to observed values. Based on statistical results, the Merced PRMS model satisfactorily simulated snowmelt seasonal streamflows. April–July calibrations for all areas had small negative bias (not greater than 7 percent) and low relative error (less than 8 percent). Less satisfactory performance for other seasons was attributed to several factors: (1) high uncertainty in low or zero flows in summer and fall, (2) lack of accounting for basin withdrawals and anthropogenic water use, (3) unavailability and (or) inaccuracy of observed (measured) meteorological input data, and (4) uncertainty in reconstructed streamflow data.With some additional refinement, the Merced PRMS model may be used for forecasting seasonal and longer-term streamflow variations; evaluating forecasted and past climate and land cover changes; providing water-resource managers with a consistent and documented method for estimating streamflow at ungaged sites within the basin; and aiding environmental studies, hydraulic design, water management, and water-quality projects in the Merced River Basin.

First posted August 6, 2019 For additional information, contact: Director, Oregon Water Science CenterU.S. Geological Survey2130 SW 5th AvenuePortland, Oregon 97201 Accurate forecasts of the streamflow expected during late spring and summer in the Upper Klamath River Basin in southern-central Oregon and northern California are used by water management agencies to balance water allocations for agriculture, aquatic habitat, and hydropower-production needs. Streamflow forecasts are also used by irrigation farmers for planning. The forecasts are typically made twice a month starting as early in the water year as December. Multiple regression equations relating real-time snowpack and precipitation conditions to seasonal streamflow volumes have been used for many years in forecasting. However, with warming temperature trends and lower snowpack, such forecasts based on historical data could become less reliable in the future. If the timing and relation of snowpack and precipitation are outside of the range of the historical data used to create the equations, the forecasts become extrapolations. Statistical forecast equations are also limited in their ability to forecast streamflow in groundwater-dominated basins having inter-annual lag. As an additional method for seasonal streamflow forecasting, a physical-process-based hydrologic model employing the Precipitation-Runoff Modeling System (PRMS) was developed in cooperation with the U.S. Bureau of Reclamation for the Upper Klamath Basin in this study. The model was calibrated for the portion of the basin draining into Upper Klamath Lake. PRMS is a deterministic, distributed-parameter, physical-process-based modeling system developed by the U.S. Geological Survey. It simulates daily streamflow, snow, solar radiation, evapotranspiration, surface-water, and groundwater processes within the basin. A model calibration and validation period for water years 2000–15 and water years 1984–99, respectively, was used. The model was calibrated and validated using measured streamflow, snowpack, evapotranspiration, and solar radiation data sets. Interpolated daily precipitation and air temperature data from 32 meteorological stations within and surrounding the Upper Klamath Basin were used as model input. Performance statistics, used to evaluate how well simulated daily streamflow matched with measured streamflow included percent bias, percent relative error, and root-mean-square error. The statistics were computed annually, monthly, for October–March, and for April–September. With the exception of the October–March period, percent bias statistics were all within plus or minus 5-percent for both the calibration and validation periods. Limitations to using the model are error in the precipitation and air temperature input time series data, which include measurement error and error in the spatial interpolation method. Other errors include measured daily streamflow data, which were adjusted for consumptive use losses to make them more closely resemble natural streamflow for calibration.The model developed for the Upper Klamath Basin can be used to forecast streamflow from the Sprague and Williamson River Basins and inflow to Upper Klamath Lake. Reliable forecasts at these locations are needed for managing water for irrigation, ecosystem health, and power production. Using the models in a forecast application requires assembling model input data sets of anticipated daily precipitation and minimum and maximum air temperature for the period after the date the forecast is made and the end of the forecasted period. These climate data sets can be based on historical or synthetic records, at the discretion of the forecaster. With the Ensemble Streamflow Prediction method, a suite of streamflow scenarios is simulated using multiple years of climate data as model input. The forecasted streamflow is determined from knowing the exceedance probabilities of the simulated streamflows. In this study, the model and the Ensemble Streamflow Prediction method were used to forecast the volume of inflow to Upper Klamath Lake for a 6-month period from April 1, 2015, to September 30, 2015, using a range of climate data sets based on El Niño Southern Oscillation (ENSO) criteria. Because 2015 was a warm phase ENSO period, climate data for 10 warm phase ENSO years from 1980 to 2010 were used as input to the model. The simulated April–September 2015 UKL inflow volume based on measured 2015 climate data was 482,000 acre-feet, which was very close to the 50th percent exceedance probability computed from 10 simulated scenarios that used warm phase ENSO climate input data from 1980–2010.

This chapter describes how the constructed basin hydrologic models can be used in water-quality analysis; using descriptions of the watershed and channel-hydraulic components of the models, specific streamflow-routing network applications, and a hypothetical application for future use in water quality modeling. Basin modeling is part of a larger effort to build a comprehensive picture of water quality in the Willamette River Basin. Precipitation-Runoff Modeling System is a physical-process, deterministic, distributed-parameter modeling system designed to analyze the effects of precipitation, climate, and land-use on streamflow and general basin hydrology. The land use, slope, soils, and geology coverages were merged to create basinwide hydrologic response unit coverage for each major basin and subbasin in the Willamette River Basin. A one-dimensional water-quality model based on the Lagrangian reference frame was developed for use in simulating the transport of conservative constituents and applying reactions among constituents for branched river systems, tidal canal systems, and deltaic channels.

The U.S. Geological Survey (USGS) has developed the PRObability of Streamflow PERmanence (PROSPER) model, a GIS raster-based empirical model that provides streamflow permanence probabilities (probabilistic predictions) of a stream channel having year-round flow for any unregulated and minimally-impaired stream channel in the Pacific Northwest region, U.S. The model provides annual predictions for 2004-2016 at a 30-m spatial resolution based on monthly or annually updated values of climatic conditions and static physiographic variables associated with the upstream basin. Predictions correspond to any pixel on the channel network consistent with the medium resolution National Hydrography Dataset channel network stream grid. Total annual precipitation and percent forest cover were consistently the most important predictor variables among global and most subregional models, which had error rates between 17 and 22%. Probabilities were converted to wet and dry streamflow permanence classes with an associated confidence. Wet and dry classifications were used to derive descriptors that characterize the statistical and spatial distribution of streamflow permanence in three focal basins. Predicted dry channel segments account for 52-92% of the stream network across the three focal basins; streamflow permanence decreased during climatically drier years. Predictions are publicly available through the USGS StreamStats platform. Results demonstrate the utility of the PROSPER model as a tool for identifying areas that may be resilient or sensitive to drought conditions, allowing for management efforts that target protecting critical reaches. Importantly, PROSPER's successful predictive performance can be improved with new datasets of streamflow permanence underscoring the importance of field observations. (C) 2018 The Authors. Published by Elsevier B.V.

First posted May 31, 2017 For additional information, contact: Director, Oregon Water Science Center U.S. Geological Survey 2130 SW 5th Avenue Portland, Oregon 97201 This report presents a summary of the hydrogeology of the upper Umatilla River Basin, Oregon, based on characterization of the hydrogeologic framework, horizontal and vertical directions of groundwater flow, trends in groundwater levels, and components of the groundwater budget. The conceptual model of the groundwater flow system integrates available data and information on the groundwater resources of the upper Umatilla River Basin and provides insights regarding key hydrologic processes, such as the interaction between the groundwater and surface water systems and the hydrologic budget.The conceptual groundwater model developed for the study area divides the groundwater flow system into five hydrogeologic units: a sedimentary unit, three Columbia River basalt units, and a basement rock unit. The sedimentary unit, which is not widely used as a source of groundwater in the upper basin, is present primarily in the lowlands and consists of conglomerate, loess, silt and sand deposits, and recent alluvium. The Columbia River Basalt Group is a series of Miocene flood basalts that are present throughout the study area. The basalt is uplifted in the southeastern half of the study area, and either underlies the sedimentary unit, or is exposed at the surface. The interflow zones of the flood basalts are the primary aquifers in the study area. Beneath the flood basalts are basement rocks composed of Paleogene to Pre-Tertiary sedimentary, volcanic, igneous, and metamorphic rocks that are not used as a source of groundwater in the upper Umatilla River Basin.The major components of the groundwater budget in the upper Umatilla River Basin are (1) groundwater recharge, (2) groundwater discharge to surface water and wells, (3) subsurface flow into and out of the basin, and (4) changes in groundwater storage.Recharge from precipitation occurs primarily in the upland areas of the Blue Mountains. Mean annual recharge from infiltration of precipitation for the upper Umatilla River Basin during 1951–2010 is about 9.6 inches per year (in/yr). Annual recharge from precipitation for water year 2010 ranged from 3 in. in the lowland area to about 30 in. in the Blue Mountains. Using Kahle and others (2011) data and methods from the Columbia Plateau regional model, average annual recharge from irrigation is estimated to be about 2.2 in/yr for the 13 square miles of irrigated land in the upper Umatilla River Basin.Groundwater discharges to streams throughout the year and is a large component of annual streamflow in the upper Umatilla River Basin. Upward vertical hydraulic gradients near the Umatilla River indicate the potential for groundwater discharge. Groundwater discharge to the Umatilla River generally occurs in the upper part of the basin, upstream from the main stem.Groundwater development in the upper Umatilla River Basin began sometime after 1950 (Davies-Smith and others, 1988; Gonthier and Bolke, 1991). By water year 2010, groundwater use in the upper Umatilla River Basin was approximately 11,214 acre-feet (acre-ft). Total groundwater withdrawals for the study area were estimated at 7,575 acre-ft for irrigation, 3,173 acre-ft for municipal use, and 466 acre-ft for domestic use.Total groundwater flow into or from the study area depends locally on geology and hydraulic head distribution. Estimates of subsurface flow were calculated using the U.S. Geological Survey Columbia Plateau regional groundwater flow model. Net flux values range from 25,000 to 27,700 acre-ft per year and indicate that groundwater is moving out of the upper Umatilla River Basin into the lower Umatilla River Basin.Water level changes depend on storage changes within an aquifer, and storage changes depend on the storage properties of the aquifer, as well as recharge to or discharge from the aquifer. Groundwater level data in the upper Umatilla River Basin are mostly available from wells in Columbia River basalt units, which indicate areas of long-term water level declines in the Grande Ronde basalt unit near Pendleton and Athena, Oregon. Groundwater levels in the Wanapum basalt unit do not show long-term declines in the upper Umatilla River Basin. Because of pumping, some areas in the upper Umatilla River Basin have shown a decrease, or reversal, in the upward vertical head gradient.Key data needs are improvement of the spatial and temporal distribution of water-level data collection and continued monitoring of streamflow gaging sites. Additionally, refinement of recharge estimates would enhance understanding of the processes that provide the groundwater resources in the upper Umatilla River Basin.

First posted October 20, 2017 For additional information, contact: Director, Oregon Water Science Center U.S. Geological Survey 2130 SW 5th Avenue Portland, Oregon 97201 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.

A three-dimensional integrated groundwater/surface-water model was developed for the upper Deschutes Basin in central Oregon to better understand groundwater and surface-water flow and interconnections, and to provide the capability for simulating the hydrologic response to external stresses including climate change. The model was calibrated using groundwater and streamflow observation data collected from 1980 through 2013. This data release contains all of the input and output files for simulation of historic conditions (the base run) plus three hypothetical pumping scenarios described in the associated USGS report: Simulation of Groundwater and Surface-Water Flow in the Upper Deschutes Basin, Oregon by Marshall W. Gannett, Kenneth E. Lite, Jr., John C. Risley, Esther M. Pischel, and Jonathan L. La Marche, U.S. Geological Survey Scientific Investigations Report 2017-5097 (https://doi.org/10.3133/sir20175097). This data release includes the GSFLOW source code.

A previously calibrated two-dimensional hydrodynamic and water-quality model (CE-QUAL-W2) of Detroit Lake in western Oregon was used in conjunction with inflows derived from Precipitation Runoff Modeling System (PRMS) hydrologic models to examine in-lake and downstream water temperature effects under future climate conditions. Current and hypothetical operations and structures at Detroit Dam were imposed on boundary conditions derived from downscaled General Circulation Models in base (1990-1999) and future (2059-2068) periods. Compared with the base period, future air temperatures were about 2 degrees C warmer year-round. Higher air temperature and lower precipitation under the future period resulted in a 23% reduction in mean annual PRMS-simulated discharge and a 1 degrees C increase in mean annual estimated stream temperatures flowing into the lake compared to the base period. Simulations incorporating current operational rules and minimum release rates at Detroit Dam to support downstream habitat, irrigation, and water supply during key times of year resulted in lower future lake levels. That scenario results in a lake level that is above the dam's spillway crest only about half as many days in the future compared to historical frequencies. Managing temperature downstream of Detroit Dam depends on the ability to blend warmer water from the lake's surface with cooler water from deep in the lake, and the spillway is an important release point near the lake's surface. Annual average in-lake and release temperatures from Detroit Lake warmed 1.1 degrees C and 1.5 degrees C from base to future periods under present-day dam operational rules and fill schedules. Simulated dam operations such as beginning refill of the lake 30 days earlier or reducing minimum release rates (to keep more water in the lake to retain the use of the spillway) mitigated future warming to 0.4 and 0.9 degrees C below existing operational scenarios during the critical autumn spawning period for endangered salmonids. A hypothetical floating surface withdrawal at Detroit Dam improved temperature control in summer and autumn (0.6 degrees C warmer in summer, 0.6 degrees C cooler in autumn compared to existing structures) without altering release rates or lake level management rules. Published by Elsevier B.V.

First posted March 31, 2016 For additional information, contact: Director, Western Fisheries Research Center U.S. Geological Survey 6505 NE 65th Street Seattle, Washington 98115 http://wfrc.usgs.gov/ We constructed a one-dimensional daily averaged water-temperature model to simulate Trinity River temperatures for 1980–2013. The purpose of this model is to assess effects of water-management actions on water temperature and to provide water temperature inputs for a salmon population dynamics model. Simulated meteorological data, observed streamflow data, and observed water temperatures were used as model inputs to simulate a continuous 34-year time series of historical daily mean water temperature at eight locations along 112.2 river miles from Lewiston Dam near Weaverville, California, downstream to the Klamath River confluence. To demonstrate the utility of the model to inform management actions, we simulated three management alternatives to assess the effects of bypass flow augmentation in a drought year, 1994, and compared those results to the simulated historical baseline, referred to as the “No Action” alternative scenario. Augmentation flows from the Lewiston Dam bypass consist of temperature-controlled releases capable of cooling downstream water temperatures in hot times of the year, which can reduce the probability of disease outbreaks in fish populations. Outputs from the Trinity River water-temperature model were then used as inputs to an existing water-temperature model of the Klamath River to evaluate the effect of augmentation flow releases on water temperatures in the lower Klamath River. We structured the Trinity River water-temperature model in River Basin Model-10 (RBM10), which uses a simple equilibrium flow model, assuming discharge in each river segment on each day is transmitted downstream instantaneously. The model uses a heat-budget formulation to quantify heat flux at the air-water interface. Inputs for the heat budget are calculated from daily mean meteorological data, including net shortwave solar radiation, net longwave atmospheric radiation, air temperature, wind speed, vapor pressure, and a psychrometric constant needed to calculate the Bowen ratio. The modeling domain was divided into eight reaches ranging in length from 8.8 to 20.6 miles, which were calibrated and validated separately with observed water temperature data collected irregularly from 1980 to 2013. Root mean square errors of observed and simulated water temperatures for the eight reaches ranged from 0.25 to 1.12 degrees Celsius (°C). Mean absolute errors ranged from 0.18 to 0.89 °C. For model validation, a k-fold cross-validation technique was used. Validation root mean square error and mean absolute error for the eight reaches ranged from 0.24 to 1.11 °C and from 0.18 to 0.89 °C, respectively. Augmentation scenarios were based on historical hydrological and meteorological data, combined with prescribed flow and temperature releases from Lewiston Dam provided by the Bureau of Reclamation. Water releases were scheduled to achieve targeted flows of 2,500, 2,800, and 3,200 cubic feet per second in the lower Klamath River from mid-August through late September, coinciding with the upstream migration of adult fall-run Chinook salmon (Oncorhynchus tshawytscha). Water temperatures simulated at river mile 5.7 on the Klamath River showed a 5 °C decrease from the No Action historical baseline, which was near or greater than 23 °C when augmentation began in mid-August. Thereafter, an approximate 1 °C difference among augmentation scenarios emerged, with the decrease in water temperature commensurate to the level of augmentation. All augmentation scenarios simulated water temperatures equal to or less than 21 °C from mid-August through late September. Water temperatures equal to or greater than 23 °C are of particular interest because of a thermal threshold known to inhibit upstream migration of salmon. When temperatures exceed this approximate 23 °C threshold, Chinook salmon are known to congregate in high densities in thermal refugias and show extended residence times, which can potentially trigger epizootic outbreaks such as of Ichthyophthirius multifiliis (“Ich”) and Flavobacterium columnare (“Columnaris”) that were the causative factors of the Klamath River fish kill in 2002. A model with the ability to simulate water temperatures in response to management actions at the basin scale is a valuable asset for water managers who must make decisions about how best to use limited water resources, which directly affect the state of fisheries in the Klamath Basin.

In 2012, the U.S. Geological Survey and the Oregon Department of Transportation began a cooperative study to demonstrate use of the Stochastic Empirical Loading and Dilution Model (SELDM) for runoff-quality analyses in Oregon. SELDM can be used to estimate stormflows, constituent concentrations, and loads from the area upstream of a stormflow discharge site, from the site of interest and in the receiving waters downstream of the discharge. SELDM also can be used to assess the potential effectiveness of best management practices (BMP) for mitigating potential effects of runoff in receiving waters. Nominally, SELDM is a highway-runoff model, but it is well suited for analysis of runoff from other land uses as well. This report provides case studies and examples to demonstrate stochastic-runoff modeling concepts and to demonstrate application of the model. Basin characteristics from six Oregon highway study sites were used to demonstrate various applications of the model. The highway catchment and upstream basin drainage areas of these study sites ranged from 3.85 to 11.83 acres and from 0.16 to 6.56 square miles, respectively. The upstream basins of two sites are urbanized, and the remaining four sites are less than 5 percent impervious. SELDM facilitates analysis by providing precipitation, pre-storm streamflow, and other variables by region or from hydrologically similar sites. In Oregon, there can be large variations in precipitation and streamflow among nearby sites. Therefore, spatially interpolated geographic information system data layers containing storm-event precipitation and pre-storm streamflow statistics specific to Oregon were created for the study using Kriging techniques. Concentrations and loads of cadmium, chloride, chromium, copper, iron, lead, nickel, phosphorus, and zinc were simulated at the six Oregon highway study sites by using statistics from sites in other areas of the country. Water‑quality datasets measured at hydrologically similar basins in the vicinity of the study sites in Oregon were selected and compiled to estimate stormflow-quality statistics for the upstream basins. The quality of highway runoff and some upstream stormflow constituents were simulated by using statistical moments (average, standard deviation, and skew) of the logarithms of data. Some upstream stormflow constituents were simulated by using transport curves, which are relations between stormflow and constituent concentrations. Stochastic analyses were done by using SELDM to demonstrate use of the model and to illustrate the types of information that stochastic analyses may provide: 1. An analysis was done to demonstrate use of dilution factors as an initial reconnaissance tool for comparing relative risk among sites. 2. An analysis of hardness-dependent, water-quality criteria was done to illustrate the effects of variations in hardness and flow on the application and interpretation of such criteria. This analysis shows that hardness-dependent criteria can vary by an order of magnitude among storm events because hardness is diluted by stormflows. 3. An analysis of uncertainties in input and output values was done to demonstrate that properly selected robust datasets are needed to represent conditions at a site of interest. This analysis shows that the rate of water-quality exceedances that are measured or simulated may depend on sample size and the luck of the draw. 4. An analysis was done to demonstrate that SELDM and other Monte Carlo models may generate extreme values from input statistics, which may or may not be feasible based on physicochemical or hydrological limits. 5. An analysis of BMP modeling methods was done to demonstrate use of the model for estimating treatment requirements for meeting water-quality objectives. 6. An analysis of the use of grab sampling and nonstochastic upstream modeling methods was done to evaluate the potential effects on modeling outcomes. Additional analyses using surrogate water-quality datasets for the upstream basin and highway catchment were provided for six Oregon study sites to illustrate the risk-based information that SELDM will produce. These analyses show that the potential effects of highway runoff on receiving-water quality downstream of the outfall depends on the ratio of drainage areas (dilution), the quality of the receiving water upstream of the highway, and the concentration of the criteria of the constituent of interest. These analyses also show that the probability of exceeding a water-quality criterion may depend on the input statistics used, thus careful selection of representative values is important.

An understanding of the role of hydro-climatic and geographic regimes on regional actual evapotranspiration (AET) change is essential to improving our knowledge on predicting water availability in a changing climate. This study investigates the relationship between AET change for a 60 year period (1951–2010) and the runoff sensitivity in 255 undisturbed catchments over the US. The runoff sensitivity to climate change is simply defined as the relative magnitude between runoff and precipitation changes with time. Runoff sensitivity can readily explain the conflicting directions of AET changes under similar precipitation change. Under increasing precipitation, AET decreases when runoff is increasing more rapidly than precipitation based on the water balance. Conversely, AET increases when runoff is decreasing more rapidly than precipitation. This result indicates that runoff sensitivity to climate change is a key factor for understanding regional water availability change at the catchment scale. In addition, a stepwise multiple regression analysis and a geographically weighted regression analysis show that the portion of evergreen forest and the mean elevation of a catchment may play a secondary role in the spatial pattern of the AET change, and the relative importance of such explanatory variables may change over space.

Prior to operational changes in 2007, Detroit Dam on the North Santiam River in western Oregon had a well-documented effect on downstream water temperature that was problematic for endangered salmonid fish species. In this U.S. Geological Survey study, done in cooperation with the U.S. Army Corps of Engineers, an existing calibrated CE-QUAL-W2 model of Detroit Lake (the impounded waterbody behind Detroit Dam) was used to determine how changes in dam operation or changes to the structural release points of Detroit Dam might affect downstream water temperatures under a range of historical hydrologic and meteorological conditions. Many combinations of environmental, operational, and structural options were explored with the model. Two downstream temperature targets were used along with three sets of environmental forcing conditions representing normal, hot/dry, and cool/wet conditions. Three structural options were modeled, including the use of existing outlets, one hypothetical variable-elevation outlet such as a sliding gate, and a hypothetical combination of a floating outlet and a fixed-elevation outlet. Finally, four sets of operational guidelines were explored to gain an understanding of the effects of imposing different downstream minimum streamflows or managing the level of the lake with different timelines in autumn. Several conclusions can be made from these interim model scenarios: * Temperature targets just downstream of Detroit Dam can be met through a combination of new dam outlets or a delayed drawdown of the lake in autumn. * Spring and summer dam operations greatly affect the available release temperatures and operational flexibility later in the autumn. Releasing warm water during mid-summer tends to keep more cool water available for release in autumn. * The ability to meet downstream temperature targets during spring depends on the characteristics of the available outlets. Under existing conditions, for example, although warm water sometimes is present at the lake surface, such water may not be available for release if the lake level is either well below or well above the spillway crest in spring and early summer. * Managing lake releases to meet downstream temperature targets depends on having outlet structures that can access both (warm) lake surface water and (cold) deeper lake water throughout the year. The existing outlets at Detroit Dam do not allow near-surface waters to be released during times when the lake surface level is below the spillway (spring and autumn). * Model simulations indicate that delayed drawdown of Detroit Lake in autumn would result in better control over release temperatures. * Compared to the existing outlets at Detroit Dam, floating or sliding-gate outlet structures can provide greater control over release temperatures because they provide better access to warm water at the lake surface and cooler water at depth. This report provides interim study results to the U.S. Army Corps of Engineers. The full study will be completed in 2012.

A study by the U.S. Geological Survey (USGS) evaluated the hydrologic response to different projected carbon emission scenarios of the 21st century using a hydrologic simulation model. This study involved five major steps: (1) setup, calibrate and evaluated the Precipitation Runoff Modeling System (PRMS) model in 14 basins across the United States by local USGS personnel; (2) acquire selected simulated carbon emission scenarios from the World Climate Research Programme's Coupled Model Intercomparison Project; (3) statistical downscaling of these scenarios to create PRMS input files which reflect the future climatic conditions of these scenarios; (4) generate PRMS projections for the carbon emission scenarios for the 14 basins; and (5) analyze the modeled hydrologic response. This report presents an overview of this study, details of the methodology, results from the 14 basin simulations, and interpretation of these results. A key finding is that the hydrological response of the different geographical regions of the United States to potential climate change may be different, depending on the dominant physical processes of that particular region. Also considered is the tremendous amount of uncertainty present in the carbon emission scenarios and how this uncertainty propagates through the hydrologic simulations.

General Circulation Model simulations of future climate through 2099 project a wide range of possible scenarios. To determine the sensitivity and potential effect of long-term climate change on the freshwater resources of the United States, the U.S. Geological Survey Global Change study, "An integrated watershed scale response to global change in selected basins across the United States" was started in 2008. The long-term goal of this national study is to provide the foundation for hydrologically based climate change studies across the nation. Fourteen basins for which the Precipitation Runoff Modeling System has been calibrated and evaluated were selected as study sites. Precipitation Runoff Modeling System is a deterministic, distributed parameter watershed model developed to evaluate the effects of various combinations of precipitation, temperature, and land use on streamflow and general basin hydrology. Output from five General Circulation Model simulations and four emission scenarios were used to develop an ensemble of climate-change scenarios for each basin. These ensembles were simulated with the corresponding Precipitation Runoff Modeling System model. This fact sheet summarizes the hydrologic effect and sensitivity of the Precipitation Runoff Modeling System simulations to climate change for the Sprague River Basin near Chiloquin, Oregon.

Computer model simulations were run to determine the effects of dam removal on water temperatures along the Klamath River, located in south-central Oregon and northern California, using flow requirements defined in the 2010 Biological Opinion of the National Marine Fisheries Service. A one-dimensional, daily averaged water temperature model (River Basin Model-10) developed by the U.S. Environmental Protection Agency Region 10, Seattle, Washington, was used in the analysis. This model had earlier been configured and calibrated for the Klamath River by the U.S. Geological Survey for the U.S. Department of the Interior, Klamath Secretarial Determination to simulate the effects of dam removal on water temperatures for current (2011) and future climate change scenarios. The analysis for this report was performed outside of the scope of the Klamath Secretarial Determination process at the request of the Bureau of Reclamation Technical Services Office, Denver, Colorado.For this analysis, two dam scenarios were simulated: “dams in” and “dams out.” In the “dams in” scenario, existing dams in the Klamath River were kept in place. In the “dams out” scenario, the river was modeled as a natural stream, without the J.C. Boyle, Copco1, Copco2, and Iron Gate Dams, for the entire simulation period. Output from the two dam scenario simulations included daily water temperatures simulated at 29 locations for a 50-year period along the Klamath River between river mile 253 (downstream of Link River Dam) and the Pacific Ocean. Both simulations used identical flow requirements, formulated in the 2010 Biological Opinion, and identical climate conditions based on the period 1961–2009.Simulated water temperatures from January through June at almost all locations between J.C. Boyle Reservoir and the Pacific Ocean were higher for the “dams out” scenario than for the “dams in” scenario. The simulated mean monthly water temperature increase was highest [1.7–2.2 degrees Celsius (°C)] in May downstream of Iron Gate Dam. However, from August to December, dam removal generally cooled water temperatures. During these months, water temperatures decreased 1°C or more between Copco Lake and locations 50 miles or more downstream. The greatest mean monthly temperature decrease was 4°C in October just downstream of Iron Gate Dam. Near the ocean, the effects of dam removal were small (less than 0.2°C) for most months. However, the mean November temperature near the ocean was almost 0.5°C cooler with dam removal.