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.     

Use of Streamflow Data to Estimate Base Flow/Ground‐Water Recharge For Wisconsin1

Abstract: The average annual base flow/recharge was determined for streamflow‐gaging stations throughout Wisconsin by base‐flow separation. A map of the State was prepared that shows the average annual base flow for the period 1970‐99 for watersheds at 118 gaging stations. Trend analysis was performed on 22 of the 118 streamflow‐gaging stations that had long‐term records, unregulated flow, and provided aerial coverage of the State. The analysis found that a statistically significant increasing trend was occurring for watersheds where the primary land use was agriculture. Most gaging stations where the land cover was forest had no significant trend. A method to estimate the average annual base flow at ungaged sites was developed by multiple‐regression analysis using basin characteristics. The equation with the lowest standard error of estimate, 9.5%, has drainage area, soil infiltration and base flow factor as independent variables. To determine the average annual base flow for smaller watersheds, estimates were made at low‐flow partial‐record stations in 3 of the 12 major river basins in Wisconsin. Regression equations were developed for each of the three major river basins using basin characteristics. Drainage area, soil infiltration, basin storage and base‐flow factor were the independent variables in the regression equations with the lowest standard error of estimate. The standard error of estimate ranged from 17% to 52% for the three river basins.     

Evaluation of Ground‐Water and Surface‐Water Exchanges Using Streamflow Difference Analyses1

Abstract: In the karstic lower Flint River Basin, limestone fracturing, jointing, and subsequent dissolution have resulted in the development of extensive secondary permeability and created a system of major conduits that facilitate the exchange of water between the Upper Floridan aquifer and Flint River. Historical streamflow data from U.S. Geological Survey gaging stations located in Albany and Newton, Georgia, were used to quantify ground‐water and surface‐water exchanges within a 55.3 km section of the Flint River. Using data from 2001, we compared estimates of ground‐water flux using a time adjustment method to a water balance equation and found that these independent approaches yielded similar results. The associated error was relatively large during high streamflow when unsteady conditions prevail, but much lower during droughts. Flow reversals were identified by negative streamflow differences and verified with in situ data from temperature sensors placed inside large spring conduits. Long‐term (13 years) analysis showed negative streamflow differentials (i.e., a losing stream condition) coincided with high river stages and indicated that streamflow intrusion into the aquifer could potentially exceed 150 m³/s. Although frequent negative flow differentials were evident, the Flint River was typically a gaining stream and showed a large net increase in flow between the two gages when examined over the period 1989‐2003. Ground‐water contributions to this stream section averaged 2‐42 m³/s with a mean of 13 m³/s. The highest rate of ground‐water discharge to the Flint River occurred during the spring when regional ground‐water levels peaked following heavy winter and spring rains and corresponding rates of evapotranspiration were low. During periods of extreme drought, ground‐water contributions to the Flint River declined.     

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.     

Allocation of Streamflow Depletion Impacts under Nonlinear Conditions

A method is proposed for the equitable allocation of impacts of groundwater pumping on streamflow. The method is intended for cases in which the pumping activity of multiple entities has impacts on streamflow and these impacts are computed by perturbation. It is shown that when the response of streamflow to pumping is nonlinear, simple methods for impact calculation can fail. The proposed method is developed for the case when there are four entities that impact streamflow. The method relies on the calculation of impacts by perturbation of the simulation model from different base pumping levels. When four entities are evaluated, 16 runs of the simulation model are required. It is shown the proposed method produces estimated impacts for each individual entity that are equitable because they meet the requirement that the impacts of each entity sum to the total impacts of all entities acting together and the impacts attributed to each entity do not depend on the order of calculation. A brief example demonstrates the approach.     

CONSUMPTIVE USE OF STREAMFLOW INCREASES IN THE COLORADO RIVER BASIN1

ABSTRACT: This study examined the disposition of streamflow increases that could be created by vegetation management on forest land along the upper reaches of the Colorado River. A network optimization model was used to simulate water flow, storage, consumptive use, and loss within the entire Colorado River Basin with and without the flow increases, according to various scenarios incorporating both current and future consumptive use levels as well as existing and potential institutional constraints. Results indicate that very little of the flow increases would be consumptively used at current use levels, or even at future use levels, if water allocation institutions remain unchanged. Given future use levels and economically based water allocation institutions, up to one-half of the flow increases could be consumptively used. The timing of streamflow increases, and the institutional constraints on water allocation, often limit the potential for consumptive use of flow increases.     

Estimation of Daily Streamflow of Southeastern Coastal Plain Watersheds by Combining Estimated Magnitude and Sequence

Commonly used methods to predict streamflow at ungauged watersheds implicitly predict streamflow magnitude and temporal sequence concurrently. An alternative approach that has not been fully explored is the conceptualization of streamflow as a composite of two separable components of magnitude and sequence, where each component is estimated separately and then combined. Magnitude is modeled using the flow duration curve (FDC), whereas sequence is modeled by transferring streamflow sequence of gauged watershed(s). This study tests the applicability of the approach on watersheds ranging in size from about 25‐7,226 km² in Southeastern Coastal Plain (U.S.) with substantial surface storage of wetlands. A 19‐point regionalized FDC is developed to estimate streamflow magnitude using the three most selected variables (drainage area, hydrologic soil index, and maximum 24‐h precipitation with a recurrence interval of 100 years) by a greedy‐heuristic search process. The results of validation on four watersheds (Trent River, North Carolina: 02092500; Satilla River, Georgia: 02226500; Black River, South Carolina: 02136000; and Coosawhatchie River, South Carolina: 02176500) yielded Nash‐Sutcliffe efficiency values of 0.86‐0.98 for the predicted magnitude and 0.09‐0.84 for the predicted daily streamflow over a simulation period of 1960‐2010. The prediction accuracy of the method on two headwater watersheds at Santee Experimental Forest in coastal South Carolina was weak, but comparable to simulations by MIKE‐SHE.     

IDENTIFICATION OF CHANGES IN STREAMFLOW CHARACTERISTICS1

ABSTRACT: A set of procedures for identifying changes in selected streamflow characteristics at sites having long‐term continuous streamflow records is illustrated by using streamflow data from the Waccamaw River at Freeland, North Carolina for the 55‐year period of 1940–1994. Data were evaluated and compared to streamflow in the adjacent Lumber River Basin to determine if changes in streamflow characteristics in the Waccamaw River were localized and possibly the result of some human activity, or consistent with regional variations. Following 1963, droughts in the Waccamaw Basin seem to have been less severe than in the Lumber Basin, and the annual one‐, seven‐, and 30‐day low flows exhibited a slightly increasing trend in the Waccamaw River. Mean daily flows in the Waccamaw River at the 90 percent exceedance level (low flows) during 1985–194, a relatively dry period, were very nearly equal to flows at the same exceedance level for 1970–1979, which represents the 10‐year period between 1940 and 1994 with the highest flows. Prior to the 1980s, flows per unit drainage area in the Waccamaw Basin were generally less than those in the Lumber Basin, but after 1980, the opposite was true. The ratio of base flow to runoff in the Waccamaw River may have changed relative to that in the Lumber River in the late 1970s. There was greater variability in Waccamaw River streamflow than in Lumber River flow, and flow variability in the Waccamaw River may have increased slightly during 1985–1994.     

LONG-RANGE STREAMFLOW FORCASTING USING NONPARAMETRIC REGRESSION1

ABSTRACT: A class of nonparametric procedures is developed for producing long-range streamflow forecasts. The forecasting procedures, which are based solely on daily streamflow data, utilize nonparametric regression to relate a forecast variable to a covariate variable. The forecast variable is a function of future streamflow and can take a wide variety of forms. The covariate variable is a function of antecedent streamflow. The forecasting procedures are quite flexible, both in terms of the duration of the forecast period and the types of forecast variables that can be considered. The procedures are used to develop long-term (1–4 months) forecasts of minimum daily flow of the Potomac River at Washington, D.C. This forecast information is an integral component of water management activities for the Washington, D.C. metropolitan area.     

Towards Real‐Time Continental Scale Streamflow Simulation in Continuous and Discrete Space

The National Weather Service (NWS) forecasts floods at approximately 3,600 locations across the United States (U.S.). However, the river network, as defined by the 1:100,000 scale National Hydrography Dataset‐Plus (NHDPlus) dataset, consists of 2.7 million river segments. Through the National Flood Interoperability Experiment, a continental scale streamflow simulation and forecast system was implemented and continuously operated through the summer of 2015. This system leveraged the WRF‐Hydro framework, initialized on a 3‐km grid, the Routing Application for the Parallel Computation of Discharge river routing model, operating on the NHDPlus, and real‐time atmospheric forcing to continuously forecast streamflow. Although this system produced forecasts, this paper presents a study of the three‐month nowcast to demonstrate the capacity to seamlessly predict reach scale streamflow at the continental scale. In addition, this paper evaluates the impact of reservoirs, through a case study in Texas. Validation of the uncalibrated model using observed hourly streamflow at 5,701 U.S. Geological Survey gages shows 26% demonstrate PBias ≤ |25%|, 11% demonstrate Nash‐Sutcliffe Efficiency (NSE) ≥ 0.25, and 6% demonstrate both PBias ≤ |25%| and NSE ≥ 0.25. When evaluating the impact of reservoirs, the analysis shows when reservoirs are included, NSE ≥ 0.25 for 56% of the gages downstream while NSE ≥ 0.25 for 11% when they are not. The results presented here provide a benchmark for the evolving hydrology program within the NWS and supports their efforts to develop a reach scale flood forecasting system for the country.     

EMERGENCY WATER SUPPLIES FROM GROUND WATER IN HUMID REGIONS1

ABSTRACT: This paper focuses on the development and testing of a mathematical model of an emergency ground water supply operated principally during periods of low streamflow. The process of ground water withdrawal and recharge is simulated taking account of streamflow, water demand, evapotranspiration, natural and artificial recharge and increased evapotranspiration due to artificial recharge, ground water pumpage, and streamflow contribution to pumped water. The model determines whether natural recharge is possible in less time than the return period of drought and also whether artificial recharge is needed. By simulating operation over a long period of time, the model can examine different droughts of short and long duration and can test the operating rules for ground water storage development in an area. Submodels analyze the components of the operating process including ground water flow into the stream, seepage losses, stream portion of well discharge due to induced infiltration and recharge from rainfall or water spreading. The model has been tested for areas in the humid northeastern United States.     

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.     

Analyzing the Effects of Groundwater Pumping on an Urban Stream‐Aquifer System

When the cone of influence of a pumping well reaches a nearby river, the resulting hydraulic gradient can induce enhanced seepage of streamflow into the aquifer. The rate of seepage is often modeled using analytical solutions that are simple to apply but may not reproduce field data due to mathematical assumptions not being met in the field. Furthermore, the appropriateness of such models has not been investigated in detail due to difficulty in measuring streamflow loss in the field. In this study, a field experiment was conducted on a reach of the South Platte River near Denver, Colorado to estimate pumping‐induced streamflow loss. A network of stream gauges, monitoring wells, and in situ measurements was used to observe streamflow rates, groundwater levels, and temperature to assess if pumping wells have a significant impact on streamflow, and to compare observed streamflow depletion against analytical solutions. Data collected suggest that pumping wells have a noticeable impact on streamflow. The analytical solutions proved accurate if streamflow was low and constant but performed poorly if streamflow was high and variable. Therefore, for this reach, the use of analytical solutions to predict streamflow may only be appropriate under low‐flow, constant‐flow conditions. Methods and results can be used to guide other streamflow depletion studies and to inform cases of pumping‐induced streamflow depletion, particularly in regard to water rights.     

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.     

Characterization of Drought in the South Atlantic,United States

Drought has been less extensively characterized in the humid South Atlantic compared to the arid western United States. Our objective was to characterize drought in the South Atlantic and to understand whether drought has become more severe in this region over time. Here we used monthly streamflow to characterize hydrological drought. Hydrological drought occurred when streamflow fell below the 20th percentile over three consecutive months and terminated once streamflow remained above the 20th percentile for three consecutive months. We characterized the frequency, duration, magnitude, and severity of events using the above definition. Significant changes in drought characteristics were tested with Mann‐Kendall over three periods: 1930‐2010, 1930‐1969, and 1970‐2010. We show that 71% of drought events were shorter than six months, while 7% were multiyear events. There was little evidence of trends in drought characteristics to support the claim of drought becoming more severe in the South Atlantic over the 20th Century. The one exception was a significant increase in the joint probability of nearby basins being simultaneously in drought conditions in the southern portion of the study area from 1970 to 2010. While drought characteristics have changed little through time, decreasing average streamflow in non drought periods coupled with increasing water demand provide the context within which recent multiyear drought events have produced significant stress on existing water infrastructure.     

Adaptation of Climate Model Projections of Streamflow to Account for Upstream Anthropogenic Impairments

A statistical procedure is developed to adjust natural streamflows simulated by dynamical models in downstream reaches, to account for anthropogenic impairments to flow that are not considered in the model. The resulting normalized downstream flows are appropriate for use in assessments of future anthropogenically impaired flows in downstream reaches. The normalization is applied to assess the potential effects of climate change on future water availability on the Rio Grande at a gage just above the major storage reservoir on the river. Model‐simulated streamflow values were normalized using a statistical parameterization based on two constants that relate observed and simulated flows over a 50‐year historical baseline period (1964–2013). The first normalization constant is a ratio of the means, and the second constant is the ratio of interannual standard deviations between annual gaged and simulated flows. This procedure forces the gaged and simulated flows to have the same mean and variance over the baseline period. The normalization constants can be kept fixed for future flows, which effectively assumes that upstream water management does not change in the future, or projected management changes can be parameterized by adjusting the constants. At the gage considered in this study, the effect of the normalization is to reduce simulated historical flow values by an average of 72% over an ensemble of simulations, indicative of the large fraction of natural flow diverted from the river upstream from the gage. A weak tendency for declining flow emerges upon averaging over a large ensemble, with tremendous variability among the simulations. By the end of the 21st Century the higher‐emission scenarios show more pronounced declines in streamflow.     

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.     

Spatiotemporal Evaluation of Simulated Evapotranspiration and Streamflow over Texas Using the WRF‐Hydro‐RAPID Modeling Framework

This study assesses a large‐scale hydrologic modeling framework (WRF‐Hydro‐RAPID) in terms of its high‐resolution simulation of evapotranspiration (ET) and streamflow over Texas (drainage area: 464,135 km²). The reference observations used include eight‐day ET data from MODIS and FLUXNET, and daily river discharge data from 271 U.S. Geological Survey gauges located across a climate gradient. A recursive digital filter is applied to decompose the river discharge into surface runoff and base flow for comparison with the model counterparts. While the routing component of the model is pre‐calibrated, the land component is uncalibrated. Results show the model performance for ET and runoff is aridity‐dependent. ET is better predicted in a wet year than in a dry year. Streamflow is better predicted in wet regions with the highest efficiency ~0.7. In comparison, streamflow is most poorly predicted in dry regions with a large positive bias. Modeled ET bias is more strongly correlated with the base flow bias than surface runoff bias. These results complement previous evaluations by incorporating more spatial details. They also help identify potential processes for future model improvements. Indeed, improving the dry region streamflow simulation would require synergistic enhancements of ET, soil moisture and groundwater parameterizations in the current model configuration. Our assessments are important preliminary steps towards accurate large‐scale hydrologic forecasts.     

PREDICTION OF MASTER RECESSION CURVES AND BASEFLOW RECESSIONS IN THE LUQUILLO MOUNTAINS OF PUERTO RICO1

ABSTRACT: An analysis of hydrograph recessions and rainfall data was performed to estimate the recession constants for two watersheds in the Luquillo mountains of Puerto Rico. To account for seasonal rainfall patterns, the data were grouped into dry and wet seasons. Sets of three Master Recession Curves (MRC) per season for each watershed were developed: one using the Matching Strip Method (MS) and two using variations of the Correlation Method (CM). These variations were the envelope line (CME) and the least squares regression (CMR). Other regression based analytical expressions that consider the streamflow recession as an autore‐gressive or an integrated moving average process were also applied. The regression based methods performed consistently better than the graphical ones and they proved to be faster, easier, and less subjective. The recession constants from these methods were then used to estimate the time it would take the streamflow to reach the critical Q₉₉ flow duration. Based on this study, once the streamflow reaches Q₉₀, water managers have 6 to 12 days warning before streamflow reaches critical levels.     

USE OF STREAMFLOW INCREASES FROM VEGETATION MANAGEMENT IN THE VERDE RIVER BASIN1

ABSTRACT: Although the effects of vegetation management on streamflow have been studied in many locations, the effects of augmented streamflow on downstream water users have not been carefully analyzed. This study examines the routing of streamflow increases that could be produced in the Verde River Basin of Arizona. Reservoir management and water routing to users in the Salt River Valley around Phoenix were carefully modeled. Simulation of water routing with and without vegetation modification indicates that, under current institutional conditions, less than one-half of the streamflow increase would reach consumptive users as surface water. Most of the remainder would accumulate in storage until a year of unusually heavy runoff, when it would add to reservoir spills. Under alternative scenarios, from 39 to 58 percent of the streamflow increase was delivered to consumptive users.