Trends in Western U.S. Snowpack and Related Upper Colorado River Basin Streamflow1

Miller, W. Paul and Thomas C. Piechota, 2011. Trends in Western U.S. Snowpack and Related Upper Colorado River Basin Streamflow. Journal of the American Water Resources Association (JAWRA) 47(6):1197–1210. DOI: 10.1111/j.1752‐1688.2011.00565.x Abstract: Water resource managers in the Western United States (U.S.) are currently faced with the challenge of adapting to unprecedented drought and uncertain impacts of climate change. Recent research has indicated increasing regional temperature and changes to precipitation and streamflow characteristics throughout the Western U.S. As such, there is increased uncertainty in hydroclimatological forecasts, which impact reservoir operations and water availability throughout the Western U.S., particularly in the Colorado River Basin. Previous research by the authors hypothesized a change in the character of precipitation (i.e., the frequency and amount of rainfall and snowfall events) throughout the Colorado River Basin. In the current study, 398 snowpack telemetry stations were investigated for trends in cumulative precipitation, snow water equivalent, and precipitation events. Observations of snow water equivalent characteristics were compared to observations in streamflow characteristics. Results indicate that the timing of the last day of the snow season corresponds well to the volume of runoff observed over the traditional peak flow season (April through July); conversely, the timing of the first day of the snow season does not correspond well to the volume of runoff observed over the peak flow season. This is significant to water resource managers and river forecasters, as snowpack characteristics may be indicative of a productive or unproductive runoff season.     

Relative Effects of Climate and Water Use on Base‐Flow Trends in the Lower Klamath Basin1

Abstract: Since the 1940s, snow water equivalent (SWE) has decreased throughout the Pacific Northwest, while water use has increased. Climate has been proposed as the primary cause of base‐flow decline in the Scott River, an important coho salmon rearing tributary in the Klamath Basin. We took a comparative‐basin approach to estimating the relative contributions of climatic and non‐climatic factors to this decline. We used permutation tests to compare discharge in 5 streams and 16 snow courses between “historic” (1942‐1976) and “modern” (1977‐2005) time periods, defined by cool and warm phases, respectively, of the Pacific Decadal Oscillation. April 1 SWE decreased significantly at most snow courses lower than 1,800 m in elevation and increased slightly at higher elevations. Correspondingly, base flow decreased significantly in the two streams with the lowest latitude‐adjusted elevation and increased slightly in two higher‐elevation streams. Base‐flow decline in the Scott River, the only study stream heavily utilized for irrigation, was larger than that in all other streams and larger than predicted by elevation. Based on comparison with a neighboring stream draining wilderness, we estimate that 39% of the observed 10 Mm³ decline in July 1‐October 22 discharge in the Scott River is explained by regional‐scale climatic factors. The remainder of the decline is attributable to local factors, which include an increase in irrigation withdrawal from 48 to 103 Mm³/year since the 1950s.     

If Stationarity is Dead,What Do We Do Now?1

Galloway, Gerald E., 2011. If Stationarity Is Dead, What Do We Do Now? Journal of the American Water Resources Association (JAWRA) 47(3):563‐570. DOI: 10.1111/j.1752‐1688.2011.00550.x Abstract: In January 2010, hydrologists, climatologists, engineers, and scientists met in Boulder, Colorado, to discuss the report of the death of hydrologic stationarity and the implications this might have on water resources planning and operations in the United States and abroad. For decades planners have relied on design guidance from the Interagency Advisory Committee on Water Data Bulletin 17B that was based upon the concept of stationarity. After 2½ days of discussion it became clear that the assembled community had yet to reach an agreement on whether or not to replace the assumption of stationarity with an assumption of nonstationarity or something else. Hydrologists were skeptical that data gathered to this point in the 21st Century point to any significant change in river parameters. Climatologists, on the other hand, point to climate change and the predicted shift away from current conditions to a more turbulent flood and drought filled future. Both groups are challenged to provide immediate guidance to those individuals in and outside the water community who today must commit funds and efforts on projects that will require the best estimates of future conditions. The workshop surfaced many approaches to dealing with these challenges. While there is good reason to support additional study of the death of stationarity, its implications, and new approaches, there is also a great need to provide those in the field the information they require now to plan, design, and operate today’s projects.     

Changes in Patterns of Streamflow From Unregulated Watersheds in Idaho,Western Wyoming,and Northern Nevada1

Clark, Gregory M., 2010. Changes in Patterns of Streamflow From Unregulated Watersheds in Idaho, Western Wyoming, and Northern Nevada. Journal of the American Water Resources Association (JAWRA) 46(3):486-497. DOI: 10.1111/j.1752-1688.2009.00416.x Abstract: Recent studies have identified a pattern of earlier spring runoff across much of North America. Earlier spring runoff potentially poses numerous problems, including increased risk of flooding and reduced summer water supply for irrigation, power generation, and migratory fish passage. To identify changing runoff patterns in Idaho streams, streamflow records were analyzed for 26 U.S. Geological Survey gaging stations in Idaho, western Wyoming, and northern Nevada, each with a minimum of 41 years of record. The 26 stations are located on 23 unregulated and relatively pristine streams that drain areas ranging from 28 to >35,000 km². Four runoff parameters were trend tested at each station for both the period of historical record and from 1967 through 2007. Parameters tested were annual mean streamflow, annual minimum daily streamflow, and the dates of the 25th and 50th percentiles of the annual total streamflow. Results of a nonparametric Mann-Kendall trend test revealed a trend toward lower annual mean and annual minimum streamflows at a majority of the stations, as well as a trend toward earlier snowmelt runoff. Significant downward trends over the period of historical record were most prevalent for the annual minimum streamflow (12 stations) and the 50th percentile of streamflow (11 stations). At most stations, trends were more pronounced during the period from 1967 through 2007. A regional Kendall test for water years 1967 through 2007 revealed significant regional trends in the percent change in the annual mean and annual minimum streamflows (0.67% less per year and 0.62% less per year, respectively), the 25th percentile of streamflow (12.3 days earlier), and the 50th percentile of streamflow (11.5 days earlier).     

Streamflow Response to Climate as Influenced by Geology and Elevation1

Mayer, Timothy D. and Seth W. Naman, 2011. Streamflow Response to Climate as Influenced by Geology and Elevation. Journal of the American Water Resources Association (JAWRA) 47(4):724‐738. DOI: 10.1111/j.1752‐1688.2011.00537.x Abstract: This study examines the regional streamflow response in 25 predominately unregulated basins to warmer winter temperatures and snowpack reductions over the last half century in the Klamath Basin of California and Oregon. Geologic controls of streamflow in the region result in two general stream types: surface‐dominated and groundwater‐dominated basins. Surface‐dominated basins were further differentiated into rain basins and snowmelt basins on the basis of elevation and timing of winter runoff. Streamflow characteristics and response to climate vary with stream type, as discussed in the study. Warmer winter temperatures and snowpack reductions have caused significantly earlier runoff peaks in both snowmelt and groundwater basins in the region. In the groundwater basins, the streamflow response to changes in snowpack is smoothed and delayed and the effects are extended longer in the summer. Our results indicate that absolute decreases in July‐September base flows are significantly greater, by an order of magnitude, in groundwater basins compared to surface‐dominated basins. The declines are important because groundwater basins sustain Upper Klamath Lake inflows and mainstem river flows during the typically dry summers of the area. Upper Klamath Lake April‐September net inflows have decreased an estimated 16% or 84 thousand acre‐feet (103.6 Mm³) since 1961, with the summer months showing proportionately more decline. These changes will exacerbate water supply problems for agriculture and natural resources in the region.     

Models,Assumptions, and Stakeholders: Planning for Water Supply Variability in the Colorado River Basin1

Abstract: Declining reservoir storage has raised the specter of the first water shortage on the Lower Colorado River since the completion of Glen Canyon and Hoover Dams. This focusing event spurred modeling efforts to frame alternatives for managing the reservoir system during prolonged droughts. This paper addresses the management challenges that arise when using modeling tools to manage water scarcity under variable hydroclimatology, shifting use patterns, and institutional complexity. Assumptions specified in modeling simulations are an integral feature of public processes. The policymaking and management implications of assumptions are examined by analyzing four interacting sources of physical and institutional uncertainty: inflow (runoff), depletion (water use), operating rules, and initial reservoir conditions. A review of planning documents and model reports generated during two recent processes to plan for surplus and shortage in the Colorado River demonstrates that modeling tools become useful to stakeholders by clarifying the impacts of modeling assumptions at several temporal and spatial scales. A high reservoir storage‐to‐runoff ratio elevates the importance of assumptions regarding initial reservoir conditions over the three‐year outlook used to assess the likelihood of reaching surplus and shortage triggers. An ensemble of initial condition predictions can provide more robust initial conditions estimates. This paper concludes that water managers require model outputs that encompass a full range of future potential outcomes, including best and worst cases. Further research into methods of representing and communicating about hydrologic and institutional uncertainty in model outputs will help water managers and other stakeholders to assess tradeoffs when planning for water supply variability.     

Contrasting Lumped and Distributed Hydrology Models for Estimating Climate Change Impacts on California Watersheds1

Maurer, Edwin P., Levi D. Brekke, and Tom Pruitt, 2010. Contrasting Lumped and Distributed Hydrology Models for Estimating Climate Change Impacts on California Watersheds. Journal of the American Water Resources Association (JAWRA) 46(5):1024–1035. DOI: 10.1111/j.1752-1688.2010.00473.x Abstract: We compare the projected changes to streamflows for three Sierra Nevada rivers using statistically downscaled output from 22 global climate projections. The downscaled meteorological data are used to drive two hydrology models: the Sacramento Soil Moisture Accounting model and the variable infiltration capacity model. These two models differ in their spatial resolution, computational time step, and degree and objective of calibration, thus producing significantly different simulations of current and future streamflow. However, the projected percentage changes in monthly streamflows through mid-21st Century generally did not differ, with the exceptions of streamflow during low flow months, and extreme low flows. These findings suggest that for physically based hydrology models applied to snow-dominated basins in Mediterranean climate regimes like the Sierra Nevada, California, model formulation, resolution, and calibration are secondary factors for estimating projected changes in extreme flows (seasonal or daily). For low flows, hydrology model selection and calibration can be significant factors in assessing impacts of projected climate change.     

Evaluation of Short-to-Medium Range Streamflow Forecasts Obtained Using an Enhanced Version of SRM1

Harshburger, Brian J., Karen S. Humes, Von P. Walden, Brandon C. Moore, Troy R. Blandford, and Albert Rango, 2010. Evaluation of Short-to-Medium Range Streamflow Forecasts Obtained Using an Enhanced Version of SRM. Journal of the American Water Resources Association (JAWRA) 46(3):603-617. DOI: 10.1111/j.1752-1688.2010.00437.x Abstract: As demand for water continues to escalate in the western United States, so does the need for accurate streamflow forecasts. Here, we describe a methodology for generating short-to-medium range (1 to 15 days) streamflow forecasts using an enhanced version of the Snowmelt Runoff Model (SRM), snow-covered area data derived from MODIS products, data from Snow Telemetry stations, and meteorological forecasts. The methodology was tested on three mid-elevation, snowmelt-dominated basins ranging in size from 1,600 to 3,500 km². To optimize the model performance and aid in its operational implementation, two enhancements have been made to SRM: (1) the use of an antecedent temperature index method to track snowpack cold content, and (2) the use of both maximum and minimum critical temperatures to partition precipitation into rain, snow, or a mixture of rain and snow. The comparison of retrospective model simulations with observed streamflow shows that the enhancements significantly improve the model performance. Streamflow forecasts generated using the enhanced version of the model compare well with the observed streamflow for the earlier leadtimes; forecast performance diminishes with leadtime due to errors in the meteorological forecasts. The three basins modeled in this research are typical of many mid-elevation basins throughout the American West, thus there is potential for this methodology to be applied successfully to other mountainous basins.     

Generation of Ensemble Streamflow Forecasts Using an Enhanced Version of the Snowmelt Runoff Model1

Harshburger, Brian J., Von P. Walden, Karen S. Humes, Brandon C. Moore, Troy R. Blandford, and Albert Rango, 2012. Generation of Ensemble Streamflow Forecasts Using an Enhanced Version of the Snowmelt Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(4): 643‐655. DOI: 10.1111/j.1752‐1688.2012.00642.x Abstract: As water demand increases in the western United States, so does the need for accurate streamflow forecasts. We describe a method for generating ensemble streamflow forecasts (1‐15 days) using an enhanced version of the snowmelt runoff model (SRM). Forecasts are produced for three snowmelt‐dominated basins in Idaho. Model inputs are derived from meteorological forecasts, snow cover imagery, and surface observations from Snowpack Telemetry stations. The model performed well at lead times up to 7 days, but has significant predictability out to 15 days. The timing of peak flow and the streamflow volume are captured well by the model, but the peak‐flow value is typically low. The model performance was assessed by computing the coefficient of determination (R²), percentage of volume difference (Dv%), and a skill score that quantifies the usefulness of the forecasts relative to climatology. The average R² value for the mean ensemble is >0.8 for all three basins for lead times up to seven days. The Dv% is fairly unbiased (within ±10%) out to seven days in two of the basins, but the model underpredicts Dv% in the third. The average skill scores for all basins are >0.6 for lead times up to seven days, indicating that the ensemble model outperforms climatology. These results validate the usefulness of the ensemble forecasting approach for basins of this type, suggesting that the ensemble version of SRM might be applied successfully to other basins in the Intermountain West.     

Getting From Here to Where? Flood Frequency Analysis and Climate1

Stedinger, Jery R. and Veronica W. Griffis, 2011. Getting From Here to Where? Flood Frequency Analysis and Climate. Journal of the American Water Resources Association (JAWRA) 47(3):506‐513. DOI: 10.1111/j.1752‐1688.2011.00545.x Abstract: Modeling variations in flood risk due to climate change and climate variability are a challenge to our profession. Flood‐risk computations by United States (U.S.) federal agencies follow guidelines in Bulletin 17 for which the latest update 17B was published in 1982. Efforts are underway to update that remarkable document. Additional guidance in the Bulletin as to how to address variation in flood risk over time would be welcome. Extensions of the log‐Pearson type 3 model to include changes in flood risk over time would be relatively easy mathematically. Here an example of the use of a sea surface temperature anomaly to anticipate changes in flood risk from year to year in the U.S. illustrates this opportunity. Efforts to project the trend in the Mississippi River flood series beg the question as to whether an observed trend will continue unabated, has reached its maximum, or is really nothing other than climate variability. We are challenged with the question raised by Milly and others: Is stationarity dead? Overall, we do not know the present flood risk at a site because of limited flood records. If we allow for historical climate variability and climate change, we know even less. But the issue is not whether stationarity is dead – the issue is how to use all the information available to reliably forecast flood risk in the future: “Where do we go from here?”     

Stationarity of Streamflow Records and Their Influence on Bankfull Regional Curves1

Haucke, Jessica and Katherine A. Clancy, 2011. Stationarity of Streamflow Records and Their Influence on Bankfull Regional Curves. Journal of the American Water Resources Association (JAWRA) 47(6):1338–1347. DOI: 10.1111/j.1752‐1688.2011.00590.x Abstract: Bankfull regional curves, which are curves that establish relationships among channel morphology, discharge, drainage area, are used extensively for stream restoration. These curves are developed upon the assumption that streamflows maintain stationarity over the entire record. We examined this assumption in the Driftless Area of southwestern Wisconsin where agricultural soil retention practices have changed, and precipitation has increased since the 1970s. We developed a bankfull regional curve for this area using field surveys of bankfull channel performed during 2008‐2009 and annual series of peak streamflows for 10 rivers with streamflow records ranging from the 1930s to 2009. We found bankfull flows to correlate to a 1.1 return period. To evaluate gage data statistics, we used the sign test to compare our channel morphology to historic 1.5 return period discharge (Q_1.5) for five time periods: 1959‐1972, 1973‐1992, 1993‐2008, 1999‐2008, and the 1959‐2008 period of record. Analysis of the historic gage data indicated that there has been a more than 30% decline in Q_1.5 since 1959. Our research suggests that land conservation practices may have a larger impact on gaging station stationarity than annual precipitation changes do. Additionally, historic peak flow data from gages, which have records that span land conservation changes, may need to be truncated to represent current flow regimes.     

Does More Storage Give California More Water?

Increasing reservoir storage is commonly proposed to mitigate increasing water demand and provide drought reserves, especially in semiarid regions such as California. This paper examines the value of expanding surface reservoir capacity in California using hydroeconomic modeling for historical conditions, a future warm‐dry climate, and California's recently adopted policy to end groundwater overdraft. Results show expanding surface storage capacity rarely provides sizable economic value in most of California. On average, expanding facilities north of California's Delta provides some benefit in 92% of 82 years modeled under historical conditions and in 61% of years modeled in a warm‐dry climate. South of California's Delta, expanding storage capacity provides no benefits in 14% of years modeled under historical conditions and 99% of years modeled with a warm‐dry climate. Results vary across facilities between and within regions. The limited benefit of surface storage capacity expansion to statewide water supply should be considered in planning California's water infrastructure.     

A Streamflow Forecasting Framework using Multiple Climate and Hydrological Models1

Abstract: Water resources planning and management efficacy is subject to capturing inherent uncertainties stemming from climatic and hydrological inputs and models. Streamflow forecasts, critical in reservoir operation and water allocation decision making, fundamentally contain uncertainties arising from assumed initial conditions, model structure, and modeled processes. Accounting for these propagating uncertainties remains a formidable challenge. Recent enhancements in climate forecasting skill and hydrological modeling serve as an impetus for further pursuing models and model combinations capable of delivering improved streamflow forecasts. However, little consideration has been given to methodologies that include coupling both multiple climate and multiple hydrological models, increasing the pool of streamflow forecast ensemble members and accounting for cumulative sources of uncertainty. The framework presented here proposes integration and offline coupling of global climate models (GCMs), multiple regional climate models, and numerous water balance models to improve streamflow forecasting through generation of ensemble forecasts. For demonstration purposes, the framework is imposed on the Jaguaribe basin in northeastern Brazil for a hindcast of 1974‐1996 monthly streamflow. The ECHAM 4.5 and the NCEP/MRF9 GCMs and regional models, including dynamical and statistical models, are integrated with the ABCD and Soil Moisture Accounting Procedure water balance models. Precipitation hindcasts from the GCMs are downscaled via the regional models and fed into the water balance models, producing streamflow hindcasts. Multi‐model ensemble combination techniques include pooling, linear regression weighting, and a kernel density estimator to evaluate streamflow hindcasts; the latter technique exhibits superior skill compared with any single coupled model ensemble hindcast.     

Projections of 21st Century Sierra Nevada Local Hydrologic Flow Components Using an Ensemble of General Circulation Models1

Abstract: Sierra Nevada snowmelt and runoff is a key source of water for many of California’s 38 million residents and nearly the entire population of western Nevada. The purpose of this study was to assess the impacts of expected 21st Century climatic changes in the Sierra Nevada at the subwatershed scale, for all hydrologic flow components, and for a suite of 16 General Circulation Models (GCMs) with two emission scenarios. The Soil and Water Assessment Tool (SWAT) was calibrated and validated at 35 unimpaired streamflow sites. Results show that temperatures are projected to increase throughout the Sierra Nevada, whereas precipitation projections vary between GCMs. These climatic changes drive a decrease in average annual streamflow and an advance of snowmelt and runoff by several weeks. The largest streamflow reductions were found in the mid‐range elevations due to less snow accumulation, whereas the higher elevation watersheds were more resilient due to colder temperatures. Simulation results showed that decreases in snowmelt affects not only streamflow, but evapotranspiration, surface, and subsurface flows, such that less water is available in spring and summer, thus potentially affecting aquatic and terrestrial ecosystems. Declining spring and summer flows did not equally affect all subwatersheds in the region, and the subwatershed perspective allowed for identification for the most sensitive basins throughout the Sierra Nevada.     

Evaluation of the Temporal Transferability of a Model Describing Dissolved Solids in Streams of the Upper Colorado River Basin1

Kenney, Terry A. and Susan G. Buto, 2012. Evaluation of the Temporal Transferability of a Model Describing Dissolved Solids in Streams of the Upper Colorado River Basin. Journal of the American Water Resources Association (JAWRA) 48(5): 1041‐1053. DOI: 10.1111/j.1752‐1688.2012.00667.x Abstract: The application of a nonlinear least‐squares regression model describing the sources and transport of dissolved solids in streams of the Upper Colorado River Basin, and that was calibrated using data from water year 1991, was evaluated for use in predicting annual dissolved‐solids loads for the years 1974 through 1998. Simulations for each water year were run using annual climate data. To evaluate how well the model captures the observed annual variability across the basin, differences in predicted annual dissolved‐solids loads for each simulated year and 1991 were compared with differences in monitored annual loads. The temporal trend of the differences between predicted annual loads for the simulated years and the load for 1991 generally followed the trend of the monitored loads. The model appears to underpredict the largest annual loads and overpredict some of the smaller annual loads. An underprediction bias for wetter years was evident in the residuals as was an overprediction bias, to a lesser degree, for drier years. A regression analysis on the residuals suggests that the underprediction bias is associated with precipitation differences from 1991 and with previously defined downward trends in dissolved‐solids concentrations in the basin. In general, given the representative climatic conditions, the model adequately performs throughout the period examined. However, the model is most transferable to years with climatic conditions similar to 1991.     

Simulated Impacts of Three Decadal Climate Variability Phenomena on Water Yields in the Missouri River Basin1

Mehta, Vikram M., Norman J. Rosenberg, and Katherin Mendoza, 2011. Simulated Impacts of Three Decadal Climate Variability Phenomena on Water Yields in the Missouri River Basin. Journal of the American Water Resources Association (JAWRA) 47(1):126‐135. DOI: 10.1111/j.1752‐1688.2010.00496.x Abstract: The Missouri River Basin (MRB) is the largest river basin in the United States (U.S.), and is one of the most important crop and livestock‐producing regions in the world. In a previous study of associations between decadal climate variability (DCV) phenomena and hydro‐meteorological (HM) variability in the MRB, it was found that positive and negative phases of the Pacific Decadal Oscillation (PDO), the tropical Atlantic sea‐surface temperature gradient variability (TAG), and the west Pacific warm pool (WPWP) temperature variability were significantly associated with decadal variability in precipitation and 2‐meter air temperature in the MRB, with combinations of various phases of these DCV phenomena associated with drought, flood, or neutral HM conditions. Here, we report on a methodology developed and applied to assess whether the aforementioned DCVs directly affect the hydrology of the MRB. The Hydrologic Unit Model of the U.S. (HUMUS) was used to simulate water yields in response to realistic values of the PDO, TAG, and WPWP at 75 widely distributed, eight‐digit hydrologic unit areas within the MRB. HUMUS driven by HM anomalies in both the positive and negative phases of the PDO and TAG resulted in major impacts on water yields, as much as ±20% of average water yield in some locations. Impacts of the WPWP were smaller. The combined and cumulative effects of these DCV phenomena on the MRB HM and water availability can be dramatic with important consequences for the MRB.