Geology topics

Adam J. Stonewall

Publications and source records attributed to Adam J. Stonewall.

At least 19 recordsLinked to original sources

USGS2026

Assessing the state of hydrologic science in the Upper Klamath Basin—A comprehensive review of data, tools, and models

Water demand in the Upper Klamath Basin (UKB) from various stakeholders and ecological needs often outstrips available supply, leading to persistent management challenges. This study reviews the state of hydrologic science within the UKB as of 2025—specifically, the tools, data, and models available for assessing five key components of the water system: (1) surface water; (2) precipitation; (3) evapotranspiration; (4) groundwater; and (5) water use. The UKB water supply is critical for Native American communities, regional agriculture, and federally listed fishes and faces challenges from competing needs, climate variability, and operational/regulatory requirements. We assess existing datasets, regional and national models, and historical studies to understand the available resources and identify gaps that may hinder integrated water assessments and management. Our findings indicate areas where improvements in data collection and model precision could improve the accuracy of water-availability forecasts and support water-management practices. This review can inform near-term forecasting, assist in optimizing water-resource data collection and management strategies, and support regional water-availability assessments of the basin.

Oregon↗

USGS2023

Assessing the effects of chloride deicer applications on groundwater near the Siskiyou Pass, southwestern Oregon, July 2018–February 2021

The U.S. Geological Survey, in cooperation with the Oregon Department of Transportation (ODOT), evaluated the effects of cold-weather chloride deicers (road deicing chemicals) on groundwater quality, with a focus on chloride, near the Siskiyou Pass in southwestern Oregon. The study covered the period during July 2018 through February 2021. Between the years 2016 and 2020 ODOT applied up to 16,000 gallons per mile of chloride deicer and 143,000 pounds per mile of road salt along an 11-mile stretch of Interstate 5 (I-5) through the Siskiyou Pass. Despite the benefit of safer driving conditions, there are potentially negative environmental effects associated with the use of chloride-based deicers (such as magnesium chloride and sodium chloride). The results from this study are intended to help ODOT assess the water-quality effects from the application of chloride deicers at the Siskiyou Pass and inform decisions on how those chemicals are used. Dissolved chloride concentrations tended to be greater in groundwater downgradient from I-5 compared to groundwater upgradient from the interstate. Specific conductance was a good predictor of dissolved chloride concentration (R 2 = 0.905). Continuous monitoring showed that specific conductance measurements were greater at four downgradient spring-fed sites at the end of the study period compared with measurements at the beginning of the study. The study results indicate that chloride levels in shallow groundwater downgradient from I-5 are increasing, but dissolved chloride concentrations in domestic wells are not above the U.S. Environmental Protection Agency drinking water recommendations. The approach and methods used in this study, with modifications as site conditions warrant, can be applied in other areas of chloride deicer application to determine if groundwater is affected.

Oregon↗

USGS2022

Assessing the impact of chloride deicer application in the Siskiyou Pass, southern Oregon

Chloride deicers have been applied by the Oregon Department of Transportation (ODOT) to Interstate Route 5 (I–5) from the Oregon-California border north to mile marker 10 for several years in the high-elevation area known as the Siskiyou Pass. Magnesium chloride (MgCl 2 ) and sodium chloride (NaCl) are applied to keep the interstate highway safe for drivers and allow for efficient transport of goods and people through adverse weather conditions, particularly snow and ice. The U.S. Geological Survey entered into a cooperative agreement with ODOT to research the effects of chloride deicers in the Carter and Wall Creek watersheds that drain the vicinity of the Siskiyou Pass. The Stochastic Empirical Loading and Dilution Model (SELDM) was used to estimate combinations of prestorm-streamflow, stormflow, highway-runoff, and event mean constituent concentrations (EMCs), as well as stormwater-constituent loads at sites of interest. The study evaluated the effects of roadway application of chloride deicers on downstream and highway-runoff conditions (particularly EMCs), exceedance rates of criterion maximum concentrations, and concurrent runoff loads of stormwater constituents from a site of interest. SELDM was also used to evaluate the efficiency of hydrograph extension best management practices to reduce peak constituent concentrations. Several SELDM scenarios were developed as sensitivity analyses to evaluate the model benefit of collecting specific local sets of data, such as streamflow, precipitation, highway-runoff and riverine water-quality samples, and volumetric runoff coefficient statistics. Results of the study showed that for SELDM modeling in the Siskiyou Pass area, (1) the inclusion of local streamflow data is important for obtaining accurate downstream EMCs, (2) the inclusion of precipitation data is important for highway and concurrent runoff load calculations, and (3) water-quality constituent EMC data from highway runoff and upstream stormflows are the most important data to collect for highway runoff and upstream water-quality constituent concentration statistics.

Oregon↗

USGS2021

Historical streamflow and stage data compilation for the Lower Columbia River, Pacific Northwest

The U.S. Geological Survey mined data from a variety of national and state agencies including USGS, Oregon Water Resources Department, National Oceanic and Atmospheric Administration, Washington Department of Ecology, Pacific Northwest National Laboratory, Portland State University, and U.S. Army Corps of Engineers. A comprehensive dataset of streamflow, stage, and tidal elevations for the Lower Columbia River basin was compiled. Data were compiled from gaging stations in Oregon and Washington along the Columbia River from Astoria to The Dalles and along the Willamette River from Salem to Portland. Tidal gages along the Washington, Oregon, and California coasts were also compiled. Seasonal maximum values were calculated for both streamflow and stage for the winter (November–March) and spring (April–July) flow seasons, as well as for the full water year when underlying data were available. The aggregated datasets are available at https://doi.org/10.5066/P9R6RT0Z .

California, Oregon, Washington↗

USGS2020

Development of regional skew coefficients for selected flood durations in the Columbia River Basin, northwestern United States and British Columbia, Canada

Flood-frequency (hereinafter frequency) estimates provide information used to design, operate, and maintain hydraulic structures such as bridges and dams. Failures of these structures could cause catastrophic loss of property, life, or both. In addition to frequency estimates that use annual peak streamflow, frequency estimates of flood durations are required to safely and effectively operate the numerous dams in the Columbia River Basin of the northwestern United States, and British Columbia, Canada. Frequency studies rely on U.S. Geological Survey Guidelines for Determining Flood Flow Frequency (Bulletin 17C, published in 2018). A major consideration in estimating frequencies is the use of skew coefficients, which measure the asymmetry of flood flow distributions. Large uncertainties are associated with estimating the at-site skew coefficients directly from streamflow records, which are limited in length. Skew also is sensitive to extreme events for limited record lengths. Bulletin 17C recommends using regional skew coefficients to weight with the at-site skew estimate for more reliable frequency estimates. In this study, streamflow records from 313 unregulated U.S. Geological Survey streamgage sites and 97 regulated sites with naturalized streamflow records provided by the U.S. Army Corps of Engineers were used to develop regional skew models for the Columbia River Basin. The naturalized streamflow records were synthesized by removing regulatory components such as withdrawals and reservoir storage. Skew models were developed for 1-, 3-, 7-, 10-, 15-, 30-, and 60-day flood durations and used to estimate regional skew coefficients for the Columbia River Basin. This report used Bayesian statistical regression methods to develop and analyze regional skew models based on hydrologically important basin characteristics. After examining a suite of available basin characteristics, mean annual precipitation had the strongest correlation to skew across the flood durations. Regional skew regression models were fit using mean annual precipitation for selected subbasins in the Columbia River Basin.

Columbia River Basin↗

USGS2020

Evaluation of restoration alternatives using hydraulic models of lake outflow at Wapato Lake National Wildlife Refuge, northwestern Oregon

Wapato Lake National Wildlife Refuge near the city of Gaston in northwestern Oregon was established in 2013, and planning is underway to restore a more natural lake and wetland system after more than 100 years of agricultural activity on the lakebed. Several water-management and restoration alternatives are under consideration, one of which involves opening and reconnecting Wapato Lake’s outlet to allow flow in and out of the lake to Wapato Creek and downstream to the Tualatin River. The effects of this and other alternatives are being evaluated, partly through a detailed examination of the lake’s water budget. The water budget for the lake during 2011–13 was quantified by the U.S. Geological Survey in partnership with U.S. Fish and Wildlife Service and others. Results were incorporated in a spreadsheet-based Water Management Scenario Tool (WMST) for Wapato Lake, which predicts the effects of various management actions on daily lake level and potential habitat areas for waterfowl or other target species. Incorporating the effects of a hypothetical open outlet between the lake and the downstream river network in the WMST was accomplished by using a hydraulic model to simulate the flow-exchange rate between Wapato Lake and Wapato Creek over a wide range of lake levels and downstream river conditions. A Hydraulic Engineering Center-River Analysis System (HEC-RAS) one-dimensional unsteady flow model was constructed and calibrated for Wapato Creek and part of the Tualatin River using data from October 2011 to April 2013, and then was used to simulate daily lake/creek exchange flows in water years 1992–2014 under hypothetically constant lake levels. Results were used to populate a table of lake/creek flow-exchange rates for use in the WMST; a dynamic link between the WMST and HEC-RAS was unrealistic because it would require hundreds of calls to HEC-RAS and result in long run times for a single water-year’s WMST calculations with daily time steps. Predictions of daily outlet flows from the WMST were checked against HEC-RAS simulated flows under daily varying lake levels to ensure that the timing and magnitude of lake/creek exchange flows used by the WMST were consistent with those of the hydraulic model. Two scenarios were tested with a hypothetical open lake outlet to show how the WMST could be used to inform restoration planning—one scenario used a year-round open lake outlet, and the other scenario closed that outlet for part of the high-water winter season. Results showed that flows in and out of a year-round open lake outlet would dominate the lake’s water budget and produce water depths during winter and through mid-summer that might be too deep to support waterbird species that require shallow water. Closing the lake outlet during large winter storms and high-water conditions in the downstream river network would isolate the lake from surrounding rivers, keep the lake level lower, and retain substantially more shallow-water areas. Because of the ease with which management alternatives can be evaluated, a water-budget spreadsheet tool such as the WMST has been a valuable part of an analysis of potential water-management and restoration alternatives for Wapato Lake National Wildlife Refuge.

Oregon↗

USGS2019

Assessing potential effects of highway and urban runoff on receiving streams in total maximum daily load watersheds in Oregon using the stochastic empirical loading and dilution model

The Stochastic Empirical Loading and Dilution Model (SELDM) was developed by the U.S. Geological Survey (USGS) in cooperation with the Federal Highway Administration to simulate stormwater quality. To assess the effects of runoff, SELDM uses a stochastic mass-balance approach to estimate combinations of pre-storm streamflow, stormflow, highway runoff, event mean concentrations (EMCs) and stormwater constituent loads from a site of interest. In addition, SELDM can be used to assess the effects of stormwater Best Management Practices (BMPs), which are designed to mitigate the adverse effects of runoff into a waterbody. Adverse effects of stormwater on receiving waters are one of the greatest unsolved water-quality problems Nationwide. State DOTs, municipalities, Federal facilities, and private property owners who manage impervious surfaces need information about the potential magnitude of their contributions and the potential effectiveness of methods to mitigate the adverse effects of runoff. Because the efficacy of at-site controls are limited, information about the potential effectiveness of alternative strategies is needed. The USGS, in cooperation with the Oregon Department of Transportation (ODOT), conducted a study to research methods in which SELDM can be used to enhance the efficiency of ODOT’s stormwater program, support the development of a stormwater banking program, and meet environmental goals. Results can be used to develop a strategic, systems-level approach to stormwater management by considering entire watersheds instead of individual road crossings. Two watersheds, Bear Creek and Mill Creek, in western Oregon were selected for analysis. Within each watershed, seven road crossings were selected for demonstrating the utility of SELDM in nested basins. Precipitation statistics, pre-storm streamflow, runoff coefficients, and hydrograph recession factors were calculated for each location and used in SELDM to simulate flow, water-quality concentrations, and constituent loads in the upstream basin, from the highway (or developed area), and downstream from the road crossing. Three water-quality constituents were selected for modeling: suspended-sediment concentration (SSC), total phosphorus (TP), and total copper (TCu). Using water-quality transport curves, the relations between streamflow and SSC and between streamflow and TP were simulated. Concentrations of TCu were simulated by configuring a linear relation between SSC and TCu. A generic BMP was simulated using the median treatment statistics for flow reductions, hydrograph extensions, concentration reductions, and minimum irreducible concentrations from nine BMP categories with data from the 2012 International BMP database. Five simulation scenarios were modeled for demonstrative purposes. These simulations were used to evaluate potential effects of different watershed properties, water-quality inputs, and stormwater mitigation measures. Instream EMCs were compared to hypothetical water-quality criteria for suspended sediment, total phosphorus, and total copper to demonstrate the concept of water-quality risk analysis. For all five scenarios, it was assumed that highway runoff concentrations were independent of location or average annual daily traffic. These five scenarios are as follows: • Simulation Scenario 1—Natural Conditions (hereafter Simulation Scenario 1) represents conditions in an undeveloped watershed. This scenario demonstrates that the strategic placement of a hypothetical road crossing within a watershed could be used to avoid exceeding water-quality standards of TP and SSC, but that no location choice results in meeting TCu standards. Implementation of BMP had the most pronounced effects on downstream water-quality constituent EMCs at road crossings with the highest ratio of highway catchment area to upstream drainage area, but the largest effect of BMP treatment on mean annual load is based on highway catchment area alone. • Simulation Scenario 2—Current Conditions (hereafter Simulation Scenario 2) represents current watershed conditions, where all developed area upstream from the road crossing was modeled as a highway and combined with the undeveloped part of the upstream drainage area (scenario 2A) and where the output from scenario 2A is used for the upstream area (developed area and the undeveloped area), and where the road crossing is added as usual (scenario 2B). Scenario 2 results indicate that attaining water-quality standards is more difficult with upstream developed areas. Specific road-crossing sites can be selected to achieve the fewest water-quality exceedances per year, but water-quality targets are not met without BMP implementation, and in some instances are not achievable even with BMP implementation. Results from this scenario also serve to quantify the upper limit of constituent reduction if funding were available to implement BMPs to large areas of development, and to quantify how much area would need BMP implementation to achieve water-quality targets. • Simulation Scenario 3—Alternative Road Layouts (hereafter Simulation Scenario 3) was designed to assess the sensitivity of SELDM to various road layouts. In this scenario, different highway configurations were superimposed at one road crossing. Results indicate that downstream waterquality constituent EMCs did not exhibit much variation, but annual water-quality constituent loads varied considerably. • Simulation Scenario 4—Varying Road Width (hereafter Simulation Scenario 4) was designed to assess the sensitivity of SELDM to road width. Similar to scenario 3, the results indicate little variation in downstream water-quality constituent EMCs, but annual water-quality constituent loads increased in proportion to road width. • Simulation scenario 5—Changes to Impervious Area (hereafter Simulation Scenario 5) was designed to investigate the effects of changing amounts of imperviousness upstream from the road crossing. Results indicate that the downstream water-quality constituent EMCs are highly correlated with the percentage of impervious area upstream.

Oregon↗

USGS2018

Assessing roadway contributions to stormwater flows, concentrations, and loads with the StreamStats application

The Oregon Department of Transportation (ODOT) and other state departments of transportation need quantitative information about the percentages of different land cover categories above any given stream crossing in the state to assess and address roadway contributions to water-quality impairments and resulting total maximum daily loads. The U.S. Geological Survey, in cooperation with ODOT and the FHWA, added roadway and land cover information to the online StreamStats application to facilitate analysis of stormwater runoff contributions from different land covers. Analysis of 25 delineated basins with drainage areas of about 100 mi2 indicates the diversity of land covers in the Willamette Valley, Oregon. On average, agricultural, developed, and undeveloped land covers comprise 15%, 2.3%, and 82% of these basin areas. On average, these basins contained about 10 mi of state highways and 222 mi of non-state roads. The Stochastic Empirical Loading and Dilution Model was used with available water-quality data to simulate long-term yields of total phosphorus from highways, non-highway roadways, and agricultural, developed, and undeveloped areas. These yields were applied to land cover areas obtained from StreamStats for the Willamette River above Wilsonville, Oregon. This analysis indicated that highway yields were larger than yields from other land covers because highway runoff concentrations were higher than other land covers and the highway is fully impervious. However, the total highway area was a fraction of the other land covers. Accordingly, highway runoff mitigation measures can be effective for managing water quality locally, they may have limited effect on achieving basin-wide stormwater reduction goals.

Transportation Research Record↗

USGS2018

Preliminary flood-duration frequency estimates using naturalized streamflow records for the Willamette River Basin, Oregon

In this study, “naturalized” daily streamflow records, created by the U.S. Army Corps of Engineers and the Bureau of Reclamation, were used to compute 1-, 3-, 7-, 10-, 15-, 30-, and 60-day annual maximum streamflow durations, which are running averages of daily streamflow for the number of days in each duration. Once the annual maximum durations were computed, the floodduration frequencies could be estimated. The estimated flood-duration frequencies correspond to the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent probabilities of their occurring or being exceeded each year. For this report, the focus was on the Willamette River Basin in Oregon, which is a subbasin of the Columbia River Basin. This study is part of a larger one encompassing the entire Columbia Basin.

Oregon↗

USGS2017

Hydrogeologic framework and selected components of the groundwater budget for the upper Umatilla River Basin, Oregon

Executive Summary 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.

Oregon↗

USGS2017

Developing flood-inundation maps for Johnson Creek, Portland, Oregon

Digital flood-inundation maps were created for a 12.9‑mile reach of Johnson Creek by the U.S. Geological Survey (USGS). The flood-inundation maps depict estimates of water depth and areal extent of flooding from the mouth of Johnson Creek to just upstream of Southeast 174th Avenue in Portland, Oregon. Each flood-inundation map is based on a specific water level and associated streamflow at the USGS streamgage, Johnson Creek at Sycamore, Oregon (14211500), which is located near the upstream boundary of the maps. The maps produced by the USGS, and the forecasted flood hydrographs produced by National Weather Service River Forecast Center can be accessed through the USGS Flood Inundation Mapper Web site ( http://wimcloud.usgs.gov/apps/FIM/FloodInundationMapper.html ) . Water-surface elevations were computed for Johnson Creek using a combined one-dimensional and two‑dimensional unsteady hydraulic flow model. The model was calibrated using data collected from the flood of December 2015 (including the calculated streamflows at two USGS streamgages on Johnson Creek) and validated with data from the flood of January 2009. Results were typically within 0.6 foot (ft) of recorded or measured water-surface elevations from the December 2015 flood, and within 0.8 ft from the January 2009 flood. Output from the hydraulic model was used to create eight flood inundation maps ranging in stage from 9 to 16 ft. Boundary condition hydrographs were identical in shape to those from the December 2015 flood event, but were scaled up or down to produce the amount of streamflow corresponding to a specific water-surface elevation at the Sycamore streamgage (14211500). Sensitivity analyses using other hydrograph shapes, and a version of the model in which the peak flow is maintained for an extended period of time, showed minimal variation, except for overbank areas near the Foster Floodplain Natural Area. Simulated water-surface profiles were combined with light detection and ranging (lidar) data collected in 2014 to delineate water-surface extents for each of the eight modeled stages. The availability of flood-inundation maps in conjunction with real-time data from the USGS streamgages along Johnson Creek and forecasted hydrographs from the National Weather Service Northwest River Forecast Center will provide residents of the watershed and emergency management personnel with valuable information that may aid in flood response, including potential evacuations, road closures, and mitigation efforts. In addition, these maps may be used for post-flood recovery efforts.

Oregon↗

USGS2016

Evaluation of flood inundation in Crystal Springs Creek, Portland, Oregon

Efforts to improve fish passage have resulted in the replacement of six culverts in Crystal Springs Creek in Portland, Oregon. Two more culverts are scheduled to be replaced at Glenwood Street and Bybee Boulevard (Glenwood/Bybee project) in 2016. Recently acquired data have allowed for a more comprehensive understanding of the hydrology of the creek and the topography of the watershed. To evaluate the impact of the culvert replacements and recent hydrologic data, a Hydrologic Engineering Center-River Analysis System hydraulic model was developed to estimate water-surface elevations during high-flow events. Longitudinal surface-water profiles were modeled to evaluate current conditions and future conditions using the design plans for the culverts to be installed in 2016. Additional profiles were created to compare with the results from the most recent flood model approved by the Federal Emergency Management Agency for Crystal Springs Creek and to evaluate model sensitivity. Model simulation results show that water-surface elevations during high-flow events will be lower than estimates from previous models, primarily due to lower estimates of streamflow associated with the 0.01 and 0.002 annual exceedance probability (AEP) events. Additionally, recent culvert replacements have resulted in less ponding behind crossings. Similarly, model simulation results show that the proposed replacement culverts at Glenwood Street and Bybee Boulevard will result in lower water-surface elevations during high-flow events upstream of the proposed project. Wider culverts will allow more water to pass through crossings, resulting in slightly higher water-surface elevations downstream of the project during high-flows than water-surface elevations that would occur under current conditions. For the 0.01 AEP event, the water-surface elevations downstream of the Glenwood/Bybee project will be an average of 0.05 ft and a maximum of 0.07 ft higher than current conditions. Similarly, for the 0.002 AEP event, the water-surface elevations will be an average of 0.04 ft and a maximum of 0.19 ft higher than current conditions.

Oregon↗

USGS2016

Development of a CE-QUAL-W2 temperature model for Crystal Springs Lake, Portland, Oregon

During summer 2014, lake level, streamflow, and water temperature in and around Crystal Springs Lake in Portland, Oregon, were measured by the U.S. Geological Survey and the City of Portland Bureau of Environmental Services to better understand the effect of the lake on Crystal Springs Creek and Johnson Creek downstream. Johnson Creek is listed as an impaired water body for temperature by the Oregon Department of Environmental Quality (ODEQ), as required by section 303(d) of the Clean Water Act. A temperature total maximum daily load applies to all streams in the Johnson Creek watershed, including Crystal Springs Creek. Summer water temperatures downstream of Crystal Springs Lake and the Golf Pond regularly exceed the ODEQ numeric criterion of 64.4 °F (18.0 °C) for salmonid rearing and migration. To better understand temperature contributions of this system, the U.S. Geological Survey developed two-dimensional hydrodynamic water temperature models of Crystal Springs Lake and the Golf Pond. Model grids were developed to closely resemble the bathymetry of the lake and pond using data from a 2014 survey. The calibrated models simulated surface water elevations to within 0.06 foot (0.02 meter) and outflow water temperature to within 1.08 °F (0.60 °C). Streamflow, water temperature, and lake elevation data collected during summer 2014 supplied the boundary and reference conditions for the model. Measured discrepancies between outflow and inflow from the lake, assumed to be mostly from unknown and diffuse springs under the lake, accounted for about 46 percent of the total inflow to the lake. Model simulations (scenarios) were run with lower water surface elevations in Crystal Springs Lake and increased shading to the lake to assess the relative effect the lake and pond characteristics have on water temperature. The Golf Pond was unaltered in all scenarios. The models estimated that lower lake elevations would result in cooler water downstream of the Golf Pond and shorter residence times in the lake. Increased shading to the lake would also provide substantial cooling. Most management scenarios resulted in a decrease in 7-day average of daily maximum values by about 2.0– 4.7 °F (1.1 –2.6 °C) for outflow from Crystal Springs Lake during the period of interest. Outflows from the Golf Pond showed a net temperature reduction of 0.5–2.7 °F (0.3–1.5 °C) compared to measured values in 2014 because of solar heating and downstream warming in the Golf Pond resulting from mixing with inflow from Reed Lake.

Oregon↗

USGS2015

Development of a HEC-RAS temperature model for the North Santiam River, northwestern Oregon

A one-dimensional, unsteady streamflow and temperature model (HEC-RAS) of the North Santiam and Santiam Rivers was developed by the U.S. Geological Survey to be used in conjunction with previously developed two-dimensional hydrodynamic water-quality models (CE-QUAL-W2) of Detroit and Big Cliff Lakes upstream of the study area. In conjunction with the output from the previously developed models, the HEC-RAS model can simulate streamflows and temperatures within acceptable limits (mean error [bias] near zero; typical streamflow errors less than 5 percent; typical water temperature errors less than 1.0 °C) for the length of the North Santiam River downstream of Big Cliff Dam under a series of potential future conditions in which dam structures and/or dam operations are modified to improve temperature conditions for threatened and endangered fish. Although a two-dimensional (longitudinal, vertical) CE-QUAL-W2 model for the North Santiam and Santiam Rivers downstream of Big Cliff Dam exists, that model proved unstable under highly variable flow conditions. The one-dimensional HEC-RAS model documented in this report can better simulate cross-sectional-averaged stream temperatures under a wide range of flow conditions. The model was calibrated using 2011 streamflow and temperature data. Measured data were used as boundary conditions when possible, although several lateral inflows and their associated water temperatures, including the South Santiam River, were estimated using statistical models. Streamflow results showed high accuracy during low-flow periods, but predictions were biased low during large storm events when unmodeled ephemeral tributaries contributed to the actual streamflow. Temperature results showed low annual bias against measured data at two locations on the North Santiam River and one location on the Santiam River. Mean absolute errors using 2011 hourly data ranged from 0.4 to 0.7 °C. Model results were checked against 2012 data and showed a positive bias at the Santiam River station (+0.6 ˚C). Annual mean absolute errors using 2012 hourly data ranged from 0.4 to 0.8 °C. Much of the error in temperature predictions resulted from the model’s inability to accurately simulate the full range of diurnal fluctuations during the warmest months. Future iterations of the model could be improved by the collection and inclusion of additional streamflow and temperature data, especially near the mouth of the South Santiam River. Presently, the model is able to predict hourly and daily water temperatures under a wide variety of conditions with a typical error of 0.8 and 0.7 °C, respectively.

Oregon↗

USGS2015

Simulations of a hypothetical temperature control structure at Detroit Dam on the North Santiam River, northwestern Oregon

Water temperature models of Detroit Lake, Big Cliff Lake, and the North Santiam River in northwestern Oregon were used to assess the potential for a hypothetical structure with variable intake elevations and an internal connection to power turbines at Detroit Dam (scenario SlidingWeir ) to release more natural, pre-dam temperatures year round. This hypothetical structure improved outflow temperature control from Detroit Dam while meeting minimum dry-season release rates and lake levels specified by the rule curve specified for Detroit Lake. A water temperature target based on long-term, without-dams temperature estimates was developed and used to guide the Detroit Lake model to blend releases from the user-defined outlets at Detroit Dam. Simulations that included warm surface water releases during the spring and summer, and cool, deep hypolimnetic water releases later during autumn typically met the temperature target. Immediately downstream of Detroit Dam, these simulations resulted in temperatures within the range of the without-dams temperature estimates for most of the year until about November. The minimum release rates of flow imposed at Detroit Dam during late summer and early autumn exceeded unregulated, without-dams flow estimates. This higher flow led to temperatures near the low end of the without-dams temperature range 46.3 river miles downstream at Greens Bridge from July to September; the high flows released from Detroit Dam were less susceptible to downstream warming than the low unregulated flows. Simulations that blended warm and cool water from different outlets at Detroit Dam resulted in less daily temperature variation compared to the without-dams scenarios as far downstream as Greens Bridge. Estimated egg-emergence days for endangered Upper Willamette River Chinook salmon ( Oncorhynchus tshawytscha ) and Upper Willamette River winter steelhead ( Oncorhynchus mykiss ) were assessed for all scenarios. Estimated spring Chinook fry emergence under SlidingWeir scenarios was 9 days later immediately downstream of Big Cliff Dam, and 4 days later at Greens Bridge compared with existing structural scenarios at Detroit Dam. Despite the inclusion of a hypothetical sliding weir at Detroit Dam, temperatures exceeded without-dams temperatures during November and December. These late-autumn exceedances likely represent the residual thermal effect of Detroit Lake operated to meet minimum dry-season release rates (supporting instream habitat and irrigation requirements) and lake levels specified by the current (2014) operating rules (supporting recreation and flood mitigation).

Oregon↗

USGS2014

Comparison of historical streamflows to 2013 Streamflows in the Williamson, Sprague, and Wood Rivers, Upper Klamath Lake Basin, Oregon

In 2013, the Upper Klamath Lake Basin, Oregon, experienced a dry spring, resulting in an executive order declaring a state of drought emergency in Klamath County. The 2013 drought limited the water supply and led to a near-total cessation of surface-water diversions for irrigation above Upper Klamath Lake once regulation was implemented. These conditions presented a unique opportunity to understand the effects of water right regulation on streamflows. The effects of regulation of diversions were evaluated by comparing measured 2013 streamflow with data from hydrologically similar years. Years with spring streamflow similar to that in 2013 measured at the Sprague River gage at Chiloquin from water years 1973 to 2012 were used to define a Composite Index Year (CIY; with diversions) for comparison to measured 2013 streamflows (no diversions). The best-fit 6 years (1977, 1981, 1990, 1991, 1994, and 2001) were used to determine the CIY. Two streams account for most of the streamflow into Upper Klamath Lake: the Williamson and Wood Rivers. Most streamflow into the lake is from the Williamson River Basin, which includes the Sprague River. Because most of the diversion regulation affecting the streamflow of the Williamson River occurred in the Sprague River Basin, and because of uncertainties about historical flows in a major diversion above the Williamson River gage, streamflow data from the Sprague River were used to estimate the change in streamflow from regulation of diversions for the Williamson River Basin. Changes in streamflow outside of the Sprague River Basin were likely minor relative to total streamflow. The effect of diversion regulation was evaluated using the “Baseflow Method,” which compared 2013 baseflow to baseflow of the CIY. The Baseflow Method reduces the potential effects of summer precipitation events on the calculations. A similar method using streamflow produced similar results, however, despite at least one summer precipitation event. The result of the analysis estimates that streamflow from the Williamson River Basin to Upper Klamath Lake increased by approximately 14,100 acre-feet between July 1 and September 30 relative to prior dry years as a result of regulation of surface-water diversions in 2013. Quantifying the change in streamflow from regulation of diversion for the Wood River Basin was likely less accurate due to a lack of long-term streamflow data. An increase in streamflow from regulation of diversions in the Wood River Basin of roughly 5,500 acre-feet was estimated by comparing the average August and September streamflow in 2013 with historical August and September streamflow. Summing the results of the estimated streamflow gain of the Williamson River Basin (14,100 acre-feet) and Wood River (5,500 acre-feet) gives a total estimated increase in streamflow into Upper Klamath Lake resulting from the July 1–September 2013 regulation of diversions of approximately 19,600 acre-feet.

Oregon↗

USGS2013

Development of CE-QUAL-W2 models for the Middle Fork Willamette and South Santiam Rivers, Oregon

Hydrodynamic (CE-QUAL-W2) models of Hills Creek Lake (HCL), Lookout Point Lake (LOP), and Dexter Lake (DEX) on the Middle Fork Willamette River (MFWR), and models of Green Peter Lake and Foster Lake on the South Santiam River systems in western Oregon were updated and recalibrated for a wide range of flow and meteorological conditions. These CE-QUAL-W2 models originally were developed by West Consultants, Inc., for the U.S. Army Corps of Engineers. This study by the U.S. Geological Survey included a reassessment of the models’ calibration in more recent years—2002, 2006, 2008, and 2011—categorized respectively as low, normal, high, and extremely high flow calendar years. These years incorporated current dam-operation practices and more available data than the time period used in the original calibration. Modeled water temperatures downstream of both HCL and LOP-DEX on the MFWR were within an average of 0.68 degree Celsius (°C) of measured values; modeled temperatures downstream of Foster Dam on the South Santiam River were within an average of 0.65°C of measured values. A new CE-QUAL-W2 model was developed and calibrated for the riverine MFWR reach between Hills Creek Dam and the head of LOP, allowing an evaluation of the flow and temperature conditions in the entire MFWR system from HCL to Dexter Dam. The complex bathymetry and long residence time of HCL, combined with the relatively deep location of the power and regulating outlet structures at Hills Creek Dam, led to a HCL model that was highly sensitive to several outlet and geometric parameters related to dam structures (STR TOP, STR BOT, STR WIDTH). Release temperatures from HCL were important and often persisted downstream as they were incorporated in the MFWR model and the LOP-DEX model (downstream of MFWR). The models tended to underpredict the measured temperature of water releases from Dexter Dam during the late-September-through-December drawdown period in 2002, and again (to a lesser extent) in 2011, but simulations were much more accurate in 2006 and 2008. This episodic model bias may have been a result of hot, dry conditions; lower lake elevations; and earlier drawdown at both HCL and LOP in 2002. These dry conditions in 2002 may have contradicted assumptions inherent in the estimation of certain model inputs, such as unmeasured inflows and water temperatures, which may respond differently during dry years than during normal and wet years. This report documents the development and calibration of new and revised flow and water-temperature models for riverine and reservoir reaches in the Middle Fork Willamette River and South Santiam River systems. Methods and model parameter values were established for the accurate simulation of flows and temperatures in these systems under current conditions. By extension, these models should be able to accurately simulate flows and temperatures under potential future conditions in which dam operations and dam outlet structures may be changed as part of a strategy to improve habitat, fish passage, and temperature conditions for endangered fish.

Oregon↗

USGS2012

Suspended-sediment characteristics for the Johnson Creek basin, Oregon, water years 2007-10

Significant Findings An analysis of suspended-sediment transport in the Johnson Creek basin, Oregon, during water years 2007–10 indicated that: Streamflow characteristics for the 4 years of study were not extremely dry or wet, and represented near-average conditions. Computed average annual suspended-sediment loads were 1,890 and 4,640 tons at the Gresham and Milwaukie stations, respectively. More than 70 percent of suspended-sediment transport in the watershed occurred during the high-flow months of November, December, and January. Less than 10 percent of suspended-sediment transport in the watershed occurred during April–October. About 50 percent of all suspended-sediment load is transported during the highest 1 percent of streamflows. The January 2009 streamflow peak was the third highest in the 70-year record for Johnson Creek. About 50 percent of suspended-sediment transport in water year 2009 occurred in January. The drainage area upstream of the Gresham streamflow-gaging station constitutes about 30 percent of the drainage area at the Milwaukie station, but accounted for about 40 percent of the suspended sediment and 45 percent of the streamflow at the Milwaukie station. On an annual basis, most of the higher sediment yield at the Gresham station, relative to the Milwaukie station, can be explained by the higher streamflow yield at the Gresham station rather than by higher suspended-sediment concentration.

Oregon↗