The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

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.     

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.     

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.     

RELATION OF NITRATE CONCENTRATIONS TO BASEFLOW IN THE RACCOON RIVER,IOWA1

ABSTRACT: Excessive nitrate‐nitrogen (nitrate) export from the Raccoon River in west central Iowa is an environmental concern to downstream receptors. The 1972 to 2000 record of daily streamflow and the results from 981 nitrate measurements were examined to describe the relation of nitrate to streamflow in the Raccoon River. No long term trends in streamflow and nitrate concentrations were noted in the 28‐year record. Strong seasonal patterns were evident in nitrate concentrations, with higher concentrations occurring in spring and fall. Nitrate concentrations were linearly related to streamflow at daily, monthly, seasonal, and annual time scales. At all time scales evaluated, the relation was improved when baseflow was used as the discharge variable instead of total streamflow. Nitrate concentrations were found to be highly stratified according to flow, but there was little relation of nitrate to streamflow within each flow range. Simple linear regression models developed to predict monthly mean nitrate concentrations explained as much as 76 percent of the variability in the monthly nitrate concentration data for 2001. Extrapolation of current nitrate baseflow relations to historical conditions in the Raccoon River revealed that increasing baseflow over the 20th century could account for a measurable increase in nitrate concentrations.     

AUTOMATED METHODS FOR ESTIMATING BASEFLOW AND GROUND WATER RECHARGE FROM STREAMFLOW RECORDS1

ABSTRACT: To quantify and model the natural ground water recharge process, six sites located in the midwest and eastern United States where previous water balance observations had been made were compared to computerized techniques to estimate: (1) base flow and (2) ground water recharge. Results from an existing automated digital filter technique for separating baseflow from daily streamflow records were compared to baseflow estimates made in the six water balance studies. Previous validation of automated baseflow separation techniques consisted only of comparisons with manual techniques. In this study, the automated digital filter technique was found to compare well with measured field estimates yielding a monthly coefficient of determination of 0.86. The recharge algorithm developed in this study is an automated derivation of the Rorabaugh hydrograph recession curve displacement method that utilizes daily streamflow. Comparison of annual recharge from field water balance measurements to those computed with the automated recession curve displacement method had coefficients of determination of 0.76 and predictive efficiencies of 71 percent. Monthly estimates showed more variation and are not advocated for use with this method. These techniques appear to be fast, reproducible methods for estimating baseflow and annual recharge and should be useful in regional modeling efforts and as a quick check on mass balance techniques for shallow water table aquifers.     

Bankfull Regional Curves for North and Northwest Florida Streams1

Abstract: Regional curves, which relate bankfull channel dimensions and discharge to watershed drainage area, are developed to aid in identifying the bankfull stage in ungaged watersheds, and estimating the bankfull discharge and dimensions for river studies and natural channel design applications. This study assessed 26 stable stream reaches in two hydro‐physiographic regions of the Florida Coastal Plain: the Northwest Florida Coastal Plain (NWFCP) and the North Florida Coastal Plain (NFCP). Data from stream reaches in Georgia and Alabama were also used to develop the Florida regional curves, since they are located in the same hydro‐physiographic region. Reaches were selected based on the presence of U.S. Geological Survey gage stations and indicators of limited watershed development (e.g., <10% impervious surface). Analyses were conducted to determine bankfull channel dimensions, bankfull discharge, average channel slope, and Rosgen stream classification. Based on these data, significant relationships were found between bankfull cross‐sectional area, width, mean depth, and discharge as a function of drainage area for both regions. Data from this study suggested that bankfull discharges and channel dimensions were larger from NWFCP streams than from Coastal Plain streams in North Carolina and Maryland. Bankfull discharges were similar between NFCP and Georgia coastal plain streams; therefore, the data were combined into one regional curve. In addition, the data were stratified by Rosgen stream type. This stratification strengthened the relationships of bankfull width and mean depth as a function of drainage area.     

MONTHLY RUNOFF DISTRIBUTION RELATED TO BASIN GEOGRAPHY1

ABSTRACT: Mean monthly runoff from ungaged drainage basins that have significant snowpacks each year can be estimated quite well by assuming that the time duration between snowfall and snowmelt is the predominant factor in temporal runoff distribution. That time span is related to basin temperatures which are, in turn, functions of basin elevation and latitude. Regional hydrologic analyses of gaged basin data create regression equations for estimating runoff distribution by month. These equations then can be applied to ungaged basins. Basin latitude and mean elevation are two independent variables that can be used in estimating monthly runoff distributions.     

Effects of Rainfall and Ground‐Water Pumping on Streamflow in Mākaha,O’ahu,Hawai’i1

Abstract: Land‐use/land‐cover changes in Mākaha valley have included the development of agriculture, residential dwellings, golf courses, potable water supply facilities, and the introduction of alien species. The impact of these changes on surface water and ground water resources in the valley is of concern. In this study, streamflow, rainfall, and ground‐water pumping data for the upper part of the Mākaha valley watershed were evaluated to identify corresponding trends and relationships. The results of this study indicate that streamflow declined during the ground‐water pumping period. Mean and median annual streamflow have declined by 42% (135 mm) and 56% (175 mm), respectively, and the mean number of dry stream days per year has increased from 8 to 125. Rainfall across the study area appears to have also declined though it is not clear whether the reduction in rainfall is responsible for all or part of the observed streamflow decline. Mean annual rainfall at one location in the study area declined by 14% (179 mm) and increased by 2% (48 mm) at a second location. Further study is needed to assess the effect of ground‐water pumping and to characterize the hydrologic cycle with respect to rainfall, infiltration, ground‐water recharge and flow in the study area, and stream base flow and storm flow.     

USE OF HISTORIC FLOOD INFORMATION IN ESTIMATING FLOOD PEAKS ON UNGAGED WATERSHEDS1

Regional procedures to estimate flood magnitudes for ungaged watersheds typically ignore available site-specific historic flood information such as high water marks and the corresponding flow estimates, otherwise referred to as limited site-specific historic (LSSH) flood data. A procedure to construct flood frequency curves on the basis of LSSH flood observations is presented. Simple inverse variance weighting is employed to systematically combine flood estimates obtained from the LSSH data base with those from a regional procedure to obtain improved estimtes of flood peaks on the ungaged watershed. For the region studied, the variance weighted estimates of flow had a lower logarithmic standard error than either the regional or the LSSH flow estimates, when compared to the estimates determined by three standard distributions for gaged watersheds investigated in the development of the methodology. Use of the simple inverse variance weighting procedure is recommended when “reliable” estimates of LSSH floods for the ungaged site are available.     

BASE FLOW RECESSION RATES,LOW FLOWS,AND HYDROLOGIC FEATURES OF SMALL WATERSHEDS IN PENNSYLVANIA,USA1

This paper examines the relationships between measurable watershed hydrologic features, base flow recession rates, and the Q_7,10 low flow statistic (the annual minimum seven‐day average streamflow occurring once every 10 years on average). Base flow recession constants were determined by analyzing hydrograph recession data from 24 small (>130 km²), unregulated watersheds across five major physiographic provinces of Pennsylvania, providing a highly variable dataset. Geomorphic, hydrogeologic, and land use parameters were determined for each watershed. The base flow recession constant was found to be most strongly correlated to drainage density, geologic index, and ruggedness number (watershed slope); however, these three parameters are intercorrelated. Multiple regression models were developed for predicting the recession rate, and it was found that only two parameters, drainage density and hydrologic soil group, were required to obtain good estimates of the recession constant. Equations were also developed to relate the recession rates to Q_7,10 per unit area, and to the Q_7,10/Q₅₀ ratio. Using these equations, estimates of base flow recession rates, Q_7,10, and streamflow reduction under drought conditions can be made for small, ungaged basins across a wide range of physiography.     

Impact of Length of Dataset on Streamflow Calibration Parameters and Performance of APEX Model

Due to resource constraints, long‐term monitoring data for calibration and validation of hydrologic and water quality models are rare. As a result, most models are calibrated and, if possible, validated using limited measured data. However, little research has been done to determine the impact of length of available calibration data on model parameterization and performance. The main objective of this study was to evaluate the impact of length of calibration data (LCD) on parameterization and performance of the Agricultural Policy Environmental eXtender model for predicting daily, monthly, and annual streamflow. Long‐term (1984‐2015) measured daily streamflow data from Rock Creek watershed, an agricultural watershed in northern Ohio, were used for this study. Data were divided into five Short (5‐year), two Medium (15‐year), and one Long (25‐year) streamflow calibration data scenarios. All LCD scenarios were calibrated and validated at three time steps: daily, monthly, and annual. Results showed LCD affected the ability of the model to accurately capture temporal variability in simulated streamflow. However, overall average streamflow, water budgets, and crop yields were simulated reasonably well for all LCD scenarios.     

Decoupling Streamflow Responses to Climate Variability and Land Use/Cover Changes in a Watershed in Northern China

Restored annual streamflow (Q_r) and measured daily streamflow of the Chaohe watershed located in northern China and associated long‐term climate and land use/cover data were used to explore the effects of land use/cover change and climate variability on the streamflow during 1961‐2009. There were no significant changes in annual precipitation (P) and potential evapotranspiration, whereas Q_r decreased significantly by 0.81 mm/yr (p < 0.001) over the study period with a change point in 1999. We used 1961‐1998 as the baseline period (BP) and 1999‐2009 the change period (CP). The mean Q_r during the CP decreased by 39.4 mm compared with that in the BP. From 1979 to 2009, the grassland area declined by 69.6%, and the forest and shrublands increased by 105.4 and 73.1%, respectively. The land use/cover change and climate variability contributed for 58.4 and 41.6% reduction in mean annual Q_r, respectively. Compared with the BP, median and high flows in the CP decreased by 38.8 and up to 75.5%, respectively. The study concludes that large‐scale ecological restoration and watershed management in northern China has greatly decreased water yield and reduced high flows due to the improved land cover by afforestation leading to higher water loss through evapotranspiration. At a large watershed scale, land use/cover change could play as much of an important role as climate variability on water resources.     

LEGAL CHALLENGE OF THE UNGAGED WATERSHED1

The lack of uniform techniques for estimating design discharges in ungaged areas is a source of growing concern in the courts now faced with challenges to floodplain boundaries and culvert design for highway crossings. This paper summarizes a court case in which calculations of the design discharge and the hydraulic backwater effects of a major highway culvert were contested by the plaintiff. Emphasis in the paper is placed on the variation in computed flows and the interpretation of the court in the face of diverse hydrologic methods for the ungaged watershed. The results of a preliminary evaluation of ungaged watershed methods applicable to Pennsylvania are also reported in terms of standard error and bias.     

SIMULATION OF FRESHWATER DISCHARGES FROM UNGAGED AREAS TO THE SEBASTIAN RIVER,FLORIDA1

ABSTRACT: Due to alterations in the natural drainage system over the past several decades and intensified agricultural practices, freshwater discharges to the Sebastian River, Florida, have increased substantially. As a result, salinity patterns in the Sebastian River and adjacent Indian River lagoon have been disrupted and the influx of nutrients has increased. Recently, the St. Johns River Water Management District has developed a 3‐D hydrodynamic and salinity model for the Sebastian River and adjacent Indian River to study the effects of freshwater inflows, and to set guidelines for management of future freshwater discharges. Freshwater inflows to the Sebastian River are part of the input data of the hydrodynamic model. Except for the downstream drainage areas, inflows are gaged, and the data were used for calibration of the hydrologic simulations. Collectively, the downstream ungaged areas constitute about 16 percent of the total drainage area. Because of the significant contribution to the total drainage area, reliable estimates of freshwater discharges from the ungaged areas to the Sebastian River are needed. This case study illustrates the development of a set of model parameters, reflecting the hydrologic and physiographic characteristics of the entire region. In this context region applies to the watersheds located in the coastal area along the Indian River from Titusville in the north to Vero Beach in the south. The parameter set was first tested on a number of gaged drainage basins in the region, and was then applied to the ungaged areas.     

ALAWAT: A spatially allocated watershed model for approximating stream,sediment, and pollutant flows in Hawaii,USA

The Ala Wai Canal Watershed Model (ALAWAT) is a planning-level watershed model for approximating direct runoff, streamflow, sediment loads, and loads for up to five pollutants. ALAWAT uses raster GIS data layers including land use, SCS soil hydrologic groups, annual rainfall, and subwatershed delineations as direct model parameter inputs and can use daily total rainfall from up to ten rain gauges and streamflow from up to ten stream gauges. ALAWAT uses a daily time step and can simulate flows for up to ten-year periods and for up to 50 subwatersheds. Pollutant loads are approximated using a user-defined combination of rating curve relationships, mean event concentrations, and loading/washoff parameters for specific subwatersheds, land uses, and times of year. Using ALAWAT, annual average streamflow and baseflow relationships and urban suspended sediment loads were approximated for the Ala Wai Canal watershed (about 10,400 acres) on the island of Oahu, Hawaii. Annual average urban suspended sediments were approximated using two methods: mean event concentrations and pollutant loading and washoff. Parameters for the pollutant loading and washoff method were then modified to simulate the effect of various street sweeping intervals on sediment loads.     

MAPPING REGIONAL PRECIPITATION INTENSITY DURATION FREQUENCY ESTIMATES1

ABSTRACT: For regional precipitation frequency analyses, methods are needed to spatially interpolate or smooth point intensity duration frequency (IDF) estimates at gage sites for the purposes of visualization and estimation at ungaged sites. In this study to update IDF estimates for Michigan, the assumption is made that for practical purposes, the entire state may be treated as a homogeneous region in which annual maximum precipitation is identically distributed at each site apart from a site‐specific scaling factor, commonly known as the index flood. Several interpolation and smoothing techniques are evaluated for IDF estimation at ungaged sites, including trend surface analysis, thin plate splines, inverse distance weighting, and several kriging algorithms. Ordinary block kriging is recommended as a practical and objective method for smoothing the variability in the index flood values and developing isopluvial maps.     

Hydrologic Landscape Classification to Estimate Bristol Bay,Alaska Watershed Hydrology

Hydrologic landscapes (HLs) have proven to be a useful tool for broad scale assessment and classification of landscapes across the United States as they help organize larger geographical areas into areas of similar hydrologic characteristics. We developed a HL classification for the Bristol Bay watershed of southwest Alaska that incorporates indices of annual climate and seasonality, terrain, geology, and the influences of large lakes and glaciers. A HL classification is particularly useful in this large watershed because of its hydrologic and landscape variability, important salmon fishery, variety of environmental and potential anthropogenic stressors, and lack of widespread hydrologic data. Following creation of Bristol Bay basin‐wide HL classes, we compared the HL distributions within watersheds grouped by two calculated runoff parameters derived from available long‐term streamflow records and found HL distributions within these groups provided predictive insight on hydrologic behavior. Using these developed runoff groups, we estimated expected hydrologic behavior in watersheds across the larger Bristol Bay watershed that lacked gauged streamflow records. The HL approach provides a scientific basis for estimating the first‐order hydrologic behavior of watersheds and landscapes that lack detailed hydrologic information.     

STREAMFLOW DEPLETION: MODELING OF REDUCED BASEFLOW ANI INDUCED STREAM INFILTRATION FROM SEASONALLY PUMPED WELLS1

ABSTRACT: Numerical modeling techniques are used to analyze streamflow depletion for stream‐aquifer systems with baseflow. The analyses calculated two flow components generated by a pumping well located at a given distance from a river that is hydraulically connected to an unconfined aquifer. The two components are induced stream infiltration and reduced baseflow; both contribute to total streamflow depletion. Simulation results suggest that the induced infiltration, the volume of water discharged from the stream to the aquifer, has a shorter term impact on streamflow, while the reduced baseflow curves show a longer term effect. The peak impacts of the two hydrologic processes on streamflow occur separately. The separate analysis helps in understanding the hydrologic interactions between stream and aquifer. Practically, it provides useful information about contaminant transport from stream to aquifer when water quality is a concern, and for areas where water quantity is an issue, the separate analysis offers additional information to the development of water resource management plan.     

GENERALIZED LOW-FLOW FREQUENCY RELATIONSHIPS FOR UNGAGED SITES IN MASSACHUSETTS1

ABSTRACT: Regional hydrologic procedures such as generalized least squares regression and streamflow record augmentation have been advocated for obtaining estimates of both flood-flow and low-flow statistics at ungaged sites. While such procedures are extremely useful in regional flood-flow studies, no evaluation of their merit in regional low-flow estimation has been made using actual streamflow data. This study develops generalized regional regression equations for estimating the d-day, T-year low-flow discharge, Q_{d, t}, at ungaged sites in Massachusetts where d = 3, 7, 14, and 30 days. A two-parameter lognormal distribution is fit to sequences of annual minimum d-day low-flows and the estimated parameters of the lognormal distribution are then related to two drainage basin characteristics: drainage area and relief. The resulting models are general, simple to use, and about as precise as most previous models that only provide estimates of a single statistic such as Q_7,10. Comparisons are provided of the impact of using ordinary least squares regression, generalized least squares regression, and streamflow record augmentation procedures to fit regional low-flow frequency models in Massachusetts.