Assessing Retrospective National Water Model Streamflow with Respect to Droughts and Low Flows in the Colorado River Basin

Streamflow monitoring in the Colorado River Basin (CRB) is essential to ensure diverse needs are met, especially during periods of drought or low flow. Existing stream gage networks, however, provide a limited record of past and current streamflow. Modeled streamflow products with more complete spatial and temporal coverage (including the National Water Model [NWM]), have primarily focused on flooding, rather than sustained drought or low flow conditions. Objectives of this study are to (1) evaluate historical performance of the NWM streamflow estimates (particularly with respect to droughts and seasonal low flows) and (2) identify characteristics relevant to model inputs and suitability for future applications. Comparisons of retrospective flows from the NWM to observed flows from the United States Geological Survey stream gage network over 22 years in the CRB reveal a tendency for underestimating low flow frequency, locations with low flows, and the number of years with low flows. We found model performance to be more accurate for the Upper CRB and at sites with higher precipitation, snow percent, baseflow index, and elevations. Underestimation of low flows and variable model performance has important implications for future applications: inaccurate evaluations of historical low flows and droughts, and less reliable performance outside of specific watershed/stream conditions. This highlights characteristics on which to focus future model development efforts.     

A SIMPLE METHOD FOR ESTIMATING BASEFLOW AT UNGAGED LOCATIONS1

ABSTRACT: Baseflow, or water that enters a stream from slowly varying sources such as ground water, can be critical to humans and ecosystems. We evaluate a simple method for estimating base‐flow parameters at ungaged sites. The method uses one or more baseflow discharge measurements at the ungaged site and longterm streamflow data from a nearby gaged site. A given baseflow parameter, such as the median, is estimated as the product of the corresponding gage site parameter and the geometric mean of the ratios of the measured baseflow discharges and the concurrent discharges at the gage site. If baseflows at gaged and ungaged sites have a bivariate lognormal distribution with high correlation and nearly equal log variances, the estimated baseflow parameters are very accurate. We tested the proposed method using long‐term streamflow data from two watershed pairs in the Driftless Area of southwestern Wisconsin. For one watershed pair, the theoretical assumptions are well met; for the other the log‐variances are substantially different. In the first case, the method performs well for estimating both annual and long‐term baseflow parameters. In the second, the method performs remarkably well for estimating annual mean and annual median baseflow discharge, but less well for estimating the annual lower decile and the long‐term mean, median, and lower decile. In general, the use of four measurements in a year is not substantially better than the use of two.     

TREE RING RECONSTRUCTIONS OF STREAMFLOW FOR THREE CANADIAN PRAIRIE RIVERS1

ABSTRACT: Information regarding long term hydrological variability is critical for the effective management of surface water resources. In the Canadian Prairie region, growing dependence on major river systems for irrigation and other consumptive uses has resulted in an increasing vulnerability to hydrological drought and growing interprovincial tension. This study presents the first dendrochronological records of streamflow for Canadian Prairie rivers. We present 1,113‐year, 522‐year, and 325‐year reconstructions of total water year (October to September) streamflow for the North Saskatchewan, South Saskatchewan, and Saskatchewan Rivers, respectively. The reconstructions indicate relatively high flows during the 20th Century and provide evidence of past prolonged droughts. Low flows during the 1840s correspond with aridity that extended over much of the western United States. Similarly, an exceptional period of prolonged low flow conditions, approximately 900 A.D. to 1300 A.D., is coincident with evidence of sustained drought across central and western North America. The 16th Century megadrought of the western United States and Mexico, however, does not appear to have had a major impact on the Canadian rivers. The dendrohydrological records illustrate the risks involved if future water policy and infrastructure development in the Canadian Prairies are based solely on records of streamflow variability over the historical record.     

Use of observed hydroclimatic trends to constrain projections of snowmelt season runoff in the Rio Grande headwaters

A method is developed for choosing 21st Century streamflow projections among widely varying results from a large ensemble of climate model-driven simulations. We quantify observed trends in climate–streamflow relationships in the Rio Grande headwaters, which has experienced warming temperature and declining snowpack since the mid-20th Century. Prominent trends in the snowmelt runoff season are used to assess corresponding statistics in downscaled global climate model projections. We define “Observationally Consistent (OC)” simulations as those that reproduce historical changes to linear statistics of diminished snowpack–streamflow coupling in the headwaters and an associated increase in the contribution of spring season (post-peak snowpack) precipitation to streamflow. Only a modest fraction of the ensemble of simulations meets these consistency metrics. The subset of OC simulations projects significant decreases in headwaters flow, whereas the simulations that poorly replicate historical trends exhibit a much wider range of projected changes. These results bolster confidence in model-based projections of declining runoff in the Rio Grande headwaters in the snowmelt runoff season and offer an example of a methodology for evaluating model-based projections in basins with similar hydroclimates that have experienced pronounced climate changes in the recent historical record.     

Millennial‐Length Records of Streamflow From Three Major Upper Colorado River Tributaries1

Gray, Stephen T., Jeffrey J. Lukas, and Connie A. Woodhouse, 2011. Millennial‐Length Records of Streamflow From Three Major Upper Colorado River Tributaries. Journal of the American Water Resources Association (JAWRA) 47(4):702‐712. DOI: 10.1111/j.1752‐1688.2011.00535.x Abstract: Drought, climate change, and shifting consumptive use are prompting a widespread reassessment of water availability in the Upper Colorado River basin. Here, we present millennial‐length records of water year (October‐September) streamflow for key Upper Colorado tributaries: the White, Yampa, and Little Snake Rivers. Based on tree rings, these records represent the first paleohydrological reconstructions from these subbasins to overlap with a series of Medieval droughts (∼ad 800 to 1300). The reconstructions show marked interannual variability imbedded in nonstationary behavior over decadal to multidecadal time scales. These reconstructions suggest that, even in a millennial context, gaged flows from a handful of years (e.g., 1977 and 2002) were extremely dry. However, droughts of much greater duration and magnitude than any in the instrumental record were regular features prior to 1900. Likewise these reconstructions point to the unusual wetness of the gage period, and the potential for recent observations to paint an overly optimistic picture of regional water supplies. The future of the Upper Colorado River will be determined by a combination of inherent hydroclimatic variability and a broad range of human‐induced changes. It is then essential that regional water managers, water users, and policy makers alike consider a broader range of hydroclimatic scenarios than is offered by the gage record alone.     

Filyos River Streamflow Reconstruction from Tree‐Ring Chronologies with Nonparametric Approaches

The article presents nonparametric methods based on K nearest neighbors (KNNs), modified KNNs, and local polynomial techniques to reconstruct streamflow ensembles from tree‐ring data in Filyos River region (Turkey). Three methods were tested using cross‐validation for the overlap period, 1963‐1997 for which the tree‐ring and streamflow data are available. It was found that for the study where the length of the overlap period was limited, a nonparametric method based on a local polynomial technique provides simulations that have a slightly better solution than the other methods. After verification using standard statistical techniques, these methods were utilized to develop streamflow reconstructions from tree‐ring data for the paleo‐hydrologic period (1657‐1963). These reconstructions of seasonal low and high flows were discussed with the obtained flood duration curve. They were also compared with the historical archives and other tree‐ring reconstructions data available in the same river. Overall, the utility and limitations of these methods and the resulting streamflow simulations were discussed to assess the long‐term discharge behavior of Filyos River and to evaluate water supply reliability.     

Continental‐Scale River Flow Modeling of the Mississippi River Basin Using High‐Resolution NHDPlus Dataset

As a key component of the National Flood Interoperability Experiment (NFIE), this article presents the continental scale river flow modeling of the Mississippi River Basin (MRB), using high‐resolution river data from NHDPlus. The Routing Application for Parallel computatIon of Discharge (RAPID) was applied to the MRB with more than 1.2 million river reaches for a 10‐year study (2005‐2014). Runoff data from the Variable Infiltration Capacity (VIC) model was used as input to RAPID. This article investigates the effect of topography on RAPID performance, the differences between the VIC‐RAPID streamflow simulations in the HUC‐2 regions of the MRB, and the impact of major dams on the streamflow simulations. The model performance improved when initial parameter values, especially the Muskingum K parameter, were estimated by taking topography into account. The statistical summary indicates the RAPID model performs better in the Ohio and Tennessee Regions and the Upper and Lower Mississippi River Regions in comparison to the western part of the MRB, due to the better performance of the VIC model. The model accuracy also increases when lakes and reservoirs are considered in the modeling framework. In general, results show the VIC‐RAPID streamflow simulation is satisfactory at the continental scale of the MRB.     

Revealing the Diversity of Natural Hydrologic Regimes in California with Relevance for Environmental Flows Applications

Alterations to flow regimes for water management objectives have degraded river ecosystems worldwide. These alterations are particularly profound in Mediterranean climate regions such as California with strong climatic variability and riverine species highly adapted to the resulting flooding and drought disturbances. However, defining environmental flow targets for Mediterranean rivers is complicated by extreme hydrologic variability and often intensive water management legacies. Improved understanding of the diversity of natural streamflow patterns and their spatial arrangement across Mediterranean regions is needed to support the future development of effective flow targets at appropriate scales for management applications with minimal resource and data requirements. Our study addresses this need through the development of a spatially explicit reach‐scale hydrologic classification for California. Dominant hydrologic regimes and their physio‐climatic controls are revealed, using available unimpaired and naturalized streamflow time‐series and generally publicly available geospatial datasets. This methodology identifies eight natural flow classes representing distinct flow sources, hydrologic characteristics, and catchment controls over rainfall‐runoff response. The study provides a broad‐scale hydrologic framework upon which flow‐ecology relationships could subsequently be established towards reach‐scale environmental flows applications in a complex, highly altered Mediterranean region.     

DAM‐INDUCED MODIFICATIONS TO UPPER ALLEGHENY RIVER STREAMFLOW PATTERNS AND THEIR BIODIVERSITY IMPLICATIONS1

ABSTRACT: This study evaluates the streamflow characteristics of the upper Allegheny River during the periods preceding (1936 to 1965) and following (1966 to 1997) completion of the Kinzua Dam in northwestern Pennsylvania. Inter‐period trends in seasonal patterns of discharge and peak flow at three downstream sites are compared to those at two upstream sites to determine the influence of this large dam on surface water hydrology. Climatic records indicate that significant changes in annual total and seasonal precipitation occurred over the twentieth century. Increased runoff during the late summer through early winter led to increased discharge both upstream and downstream during these months, while slightly less early‐year rainfall produced minor reductions in spring flood peaks since 1966. The Kinzua Dam significantly enhanced these trends downstream, creating large reductions in peak flow, while greatly augmenting low flow during the growing season. This reduction in streamflow variability, coupled with other dam‐induced changes, has important biodiversity implications. The downstream riparian zone contains numerous threatened/endangered species, many of which are sensitive to the type of habitat modifications produced by the dam. Flood dynamics under the current post‐dam conditions are likely to compound the difficulties of maintaining their long‐term viability.     

A Fully Distributed Implementation of Mean Annual StreamFlow Regional Regression Equations1

Abstract: Estimates of mean annual streamflow are needed for a variety of hydrologic assessments. Away from gage locations, regional regression equations that are a function of upstream area, precipitation, and temperature are commonly used. Geographic information systems technology has facilitated their use for projects, but traditional approaches using the polygon overlay operator have been too inefficient for national scale applications. As an alternative, the Elevation Derivatives for National Applications (EDNA) database was used as a framework for a fully distributed implementation of mean annual streamflow regional regression equations. The raster “flow accumulation” operator was used to efficiently achieve spatially continuous parameterization of the equations for every 30 m grid cell of the conterminous United States (U.S.). Results were confirmed by comparing with measured flows at stations of the Hydro‐Climatic Data Network, and their applications value demonstrated in the development of a national geospatial hydropower assessment. Interactive tools at the EDNA website make possible the fast and efficient query of mean annual streamflow for any location in the conterminous U.S., providing a valuable complement to other national initiatives (StreamStats and the National Hydrography Dataset Plus).     

Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data1

Nishat, Bushra and S.M. Mahbubur Rahman, 2009. Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data. Journal of the American Water Resources Association (JAWRA) 45(6):1313‐1327. Abstract: Large‐scale water resources modeling can provide useful insights on future water availability scenarios for downstream nations in anticipation of proposed upstream water resources projects in large international river basins (IRBs). However, model set up can be challenging due to the large amounts of data requirement on both static states (soils, vegetation, topography, drainage network, etc.) and dynamic variables (rainfall, streamflow, soil moisture, evapotranspiration, etc.) over the basin from multiple nations and data collection agencies. Under such circumstances, satellite remote sensing provides a more pragmatic and convenient alternative because of the vantage of space and easy availability from a single data platform. In this paper, we demonstrate a modeling effort to set up a water resources management model, MIKE BASIN, over the Ganges, Brahmaputra, and Meghna (GBM) river basins. The model is set up with the objective of providing Bangladesh, the lowermost riparian nation in the GBM basins, a framework for assessing proposed water diversion scenarios in the upstream transboundary regions of India and deriving quantitative impacts on water availability. Using an array of satellite remote sensing data on topography, vegetation, and rainfall from the transboundary regions, we demonstrate that it is possible to calibrate MIKE BASIN to a satisfactory level and predict streamflow in the Ganges and Brahmaputra rivers at the entry points of Bangladesh at relevant scales of water resources management. Simulated runoff for the Ganges and Brahmaputra rivers follow the trends in the rated discharge for the calibration period. However, monthly flow volume differs from the actual rated flow by (?) 8% to (+) 20% in the Ganges basin, by (?) 15 to (+) 12% in the Brahmaputra basin, and by (?) 15 to (+) 19% in the Meghna basin. Our large‐scale modeling initiative is generic enough for other downstream nations in IRBs to adopt for their own modeling needs.     

AutoRAPID: A Model for Prompt Streamflow Estimation and Flood Inundation Mapping over Regional to Continental Extents

This article couples two existing models to quickly generate flow and flood‐inundation estimates at high resolutions over large spatial extents for use in emergency response situations. Input data are gridded runoff values from a climate model, which are used by the Routing Application for Parallel computatIon of Discharge (RAPID) model to simulate flow rates within a vector river network. Peak flows in each river reach are then supplied to the AutoRoute model, which produces raster flood inundation maps. The coupled tool (AutoRAPID) is tested for the June 2008 floods in the Midwest and the April‐June 2011 floods in the Mississippi Delta. RAPID was implemented from 2005 to 2014 for the entire Mississippi River Basin (1.2 million river reaches) in approximately 45 min. Discretizing a 230,000‐km² area in the Midwest and a 109,500‐km² area in the Mississippi Delta into thirty‐nine 1° by 1° tiles, AutoRoute simulated a high‐resolution (~10 m) flood inundation map in 20 min for each tile. The hydrographs simulated by RAPID are found to perform better in reaches without influences from unrepresented dams and without backwater effects. Flood inundation maps using the RAPID peak flows vary in accuracy with F‐statistic values between 38.1 and 90.9%. Better performance is observed in regions with more accurate peak flows from RAPID and moderate to high topographic relief.     

COMPARISON OF REGRESSION AND TIME-SERIES METHODS FOR SYNTHESIZING MISSING STREAMFLOW RECORDS1

ABSTRACT: Regression and time-series techniques have been used to synthesize and predict the stream flow at the Foresta Bridge gage from information at the upstream Pohono Bridge gage on the Merced River near Yosemite National Park. Using the available data from two time periods (calendar year 1979 and water year 1986), we evaluated the two techniques in their ability to model the variation in the observed flows and in their ability to predict stream flow at the Foresta Bridge gage for the 1979 time period with data from the 1986 time period. Both techniques produced reasonably good estimates and forecasts of the flow at the downstream gage. However, the regression model was found to have a significant amount of autocorrelation in the residuals, which the time-series model was able to eliminate. The time-series technique presented can be of great assistance in arriving at reasonable estimates of flow in data sets that have large missing portions of data.     

A STOCHASTIC DYNAMIC PROGRAMMING BASED APPROACH TO THE OPERATION OF A MULTI-RESERVOIR SYSTEM1

ABSTRACT: An heuristic iterative technique based upon stochastic dynamic programming is presented for the analysis of the operation of a three reservoir ‘Y’ shaped hydroelectric system. The technique is initiated using historical inflow data for the downstream reservoir. At each iteration the optimal policies for the downstream hydroelectric generating unit are used to provide relative weightings or targets for operation of upstream reservoirs. New input inflows to the downstream reservoir are then obtained by running the historical streamflow record through the optimal policies for the upstream reservoirs. These flows are then used to develop a new operating policy for the downstream reservoir and hence new targets for the upstream reservoirs. The process is continued until the operating policies for each reservoir provide the same overall system benefit for two successive iterations. Results obtained from the procedure are compared to the results obtained by historical operation of the system. The procedure is shown to develop operating policies which give benefits which are as close to the historical benefits as can be expected given the choice of the number of storage state variables.     

Reconstructions of Columbia River Streamflow from Tree‐Ring Chronologies in the Pacific Northwest,USA

We developed Columbia River streamflow reconstructions using a network of existing, new, and updated tree‐ring records sensitive to the main climatic factors governing discharge. Reconstruction quality is enhanced by incorporating tree‐ring chronologies where high snowpack limits growth, which better represent the contribution of cool‐season precipitation to flow than chronologies from trees positively sensitive to hydroclimate alone. The best performing reconstruction (back to 1609 CE) explains 59% of the historical variability and the longest reconstruction (back to 1502 CE) explains 52% of the variability. Droughts similar to the high‐intensity, long‐duration low flows observed during the 1920s and 1940s are rare, but occurred in the early 1500s and 1630s‐1640s. The lowest Columbia flow events appear to be reflected in chronologies both positively and negatively related to streamflow, implying low snowpack and possibly low warm‐season precipitation. High flows of magnitudes observed in the instrumental record appear to have been relatively common, and high flows from the 1680s to 1740s exceeded the magnitude and duration of observed wet periods in the late‐19th and 20th Century. Comparisons between the Columbia River reconstructions and future projections of streamflow derived from global climate and hydrologic models show the potential for increased hydrologic variability, which could present challenges for managing water in the face of competing demands.     

WATERSHED SCALE MODELING OF CRITICAL SOURCE AREAS OF RUNOFF GENERATION AND PHOSPHORUS TRANSPORT1

ABSTRACT: A curve number based model, Soil and Water Assessment Tool (SWAT), and a physically based model, Soil Moisture Distribution and Routing (SMDR), were applied in a headwater watershed in Pennsylvania to identify runoff generation areas, as runoff areas have been shown to be critical for phosphorus management. SWAT performed better than SMDR in simulating daily streamflows over the four‐year simulation period (Nash‐Sutcliffe coefficient: SWAT, 0.62; SMDR, 0.33). Both models varied streamflow simulations seasonally as precipitation and watershed conditions varied. However, levels of agreement between simulated and observed flows were not consistent over seasons. SMDR, a variable source area based model, needs further improvement in model formulations to simulate large peak flows as observed. SWAT simulations matched the majority of observed peak flow events. SMDR overpredicted annual flow volumes, while SWAT underpredicted the same. Neither model routes runoff over the landscape to water bodies, which is critical to surface transport of phosphorus. SMDR representation of the watershed as grids may allow targeted management of phosphorus sources. SWAT representation of fields as hydrologic response units (HRUs) does not allow such targeted management.     

Adaptive Water Resource Planning in the South Saskatchewan River Basin: Use of Scenarios of Hydroclimatic Variability and Extremes

The South Saskatchewan River Basin is one of Canada's most threatened watersheds, with water supplies in most subbasins over‐allocated. In 2013, stakeholders representing irrigation districts, the environment, and municipalities collaborated with researchers and consultants to explore opportunities to improve the resiliency of the management of the Oldman and South Saskatchewan River subbasins. Streamflow scenarios for 2025‐2054 were constructed by the novel approach of regressing historical river flows against indices of large‐scale ocean‐atmosphere climate oscillations to derive statistical streamflow models, which were then run using projected climate indices from global climate models. The impacts of some of the most extreme scenarios were simulated using the hydrologic mass‐balance model Operational Analysis and Simulation of Integrated Systems (OASIS). Based on stakeholder observations, the project participants proposed and evaluated potential risk management and adaption strategies, e.g., modifying existing infrastructure, building new infrastructure, changing operations to supplement environmental flows, reducing demand, and sharing supply. The OASIS model was applied interactively at live modeling sessions with stakeholders to explore practical adaptation strategies. Our results, which serve as recommendations for policy makers, showed that forecast‐based rationing together with new expanded storage could dramatically reduce water shortages.     

Streamflow Responses to Climate Change: Analysis of Hydrologic Indicators in a New York City Water Supply Watershed

Recent works have indicated that climate change in the northeastern United States is already being observed in the form of shorter winters, higher annual average air temperature, and more frequent extreme heat and precipitation events. These changes could have profound effects on aquatic ecosystems, and the implications of such changes are less understood. The objective of this study was to examine how future changes in precipitation and temperature translate into changes in streamflow using a physically based semidistributed model, and subsequently how changes in streamflow could potentially impact stream ecology. Streamflow parameters were examined in a New York City water supply watershed for changes from model‐simulated baseline conditions to future climate scenarios (2081‐2100) for ecologically relevant factors of streamflow using the Indicators of Hydrologic Alterations tool. Results indicate that earlier snowmelt and reduced snowpack advance the timing and increase the magnitude of discharge in the winter and early spring (November‐March) and greatly decrease monthly streamflow later in the spring in April. Both the rise and fall rates of the hydrograph will increase resulting in increased flashiness and flow reversals primarily due to increased pulses during winter seasons. These shifts in timing of peak flows, changes in seasonal flow regimes, and changes in the magnitudes of low flow can all influence aquatic organisms and have the potential to impact stream ecology.     

ESTIMATING STREAMFLOW AND ASSOCIATED HYDRAULIC GEOMETRY,THE MID‐ATLANTIC REGION,USA1

ABSTRACT: Methods to estimate streamflow and channel hydraulic geometry were developed for unpaged streams in the Mid‐Atlantic Region. Observed mean annual streamflow and associated hydraulic geometry data from 75 gaging stations in the Appalachian Plateau, the Ridge and Valley, and the Piedmont Physiographic Provinces of the Mid‐Atlantic Region were used to develop a set of power functions that relate streamflow to drainage area and hydraulic geometry to streamflow. For all three physiographic provinces, drainage area explained 95 to 98 percent of the variance in mean annual streamflow. Relationships between mean annual streamflow and water surface width and mean flow depth had coefficients of determination that ranged from R²= 0.55 to R²= 0.91, but the coefficient of determination between mean flow velocity and mean annual streamflow was lower (R²= 0.44 to R²= 0.54). The advantages of using the regional regression models to estimate streamflow over a conceptual model or a water balance model are its ease of application and reduced input data needs. The prediction of the regression equations were tested with data collected as part of the U.S. Environmental Protection Agency (USEPA) Environmental Monitoring and Assessment Program (EMAP). In addition, equations to transfer streamflow from gaged to ungaged streams are presented.     

RESERVOIR OPERATION ANI EVALUATION OF DOWNSTREAM FLOW AUGMENTATION1

ABSTRACT: Operation of a storage‐based reservoir modifies the downstream flow usually to a value higher than that of natural flow in dry season. This could be important for irrigation, water supply, or power production as it is like an additional downstream benefit without any additional investment. This study addresses the operation of two proposed reservoirs and the downstream flow augmentation at an irrigation project located at the outlet of the Gandaki River basin in Nepal. The optimal operating policies of the reservoirs were determined using a Stochastic Dynamic Programming (SDP) model considering the maximization of power production. The modified flows downstream of the reservoirs were simulated by a simulation model using the optimal operating policy (for power maximization) and a synthetic long‐term inflow series. Comparing the existing flow (flow in river without reservoir operation) and the modified flow (flow after reservoir operation) at the irrigation project, the additional amount of flow was calculated. The reliability analysis indicated that the supply of irrigation could be increased by 25 to 100 percent of the existing supply over the dry season (January to April) with a reliability of more than 80 percent.