MONTHLY RUNOFF PREDICTIONS BASED ON RAINFALL FORECASTS IN A SMALL OKLAHOMA WATERSHED1

ABSTRACT: Conditions under which monthly rainfall forecasts translate into monthly runoff predictions that could support water resources planning and management activities were investigated on a small watershed in central Oklahoma. Runoff response to rainfall forecasts was simulated using the hydrologic model SWAT. Eighteen scenarios were examined that represented combinations of wet, average, and dry antecedent rainfall conditions, with wet, normal, and dry forecasted rainfall. Results suggest that for the climatic and physiographic conditions under consideration, rainfall forecasts could offer potential application opportunities in surface water resources but only under certain conditions. Pronounced wet and dry antecedent rainfall conditions were shown to have greater impact on runoff than forecasts, particularly in the first month of a forecast period. Large forecast impacts on runoff occurred under wet antecedent conditions, when the fraction of forecasted rainfall contributing to runoff was greatest. Under dry antecedent conditions, most of the forecasted rainfall was absorbed in the soil profile, with little immediate runoff response. Persistent three‐month forecasts produced stronger impacts on runoff than one‐month forecasts due to cumulative effects in the hydrologic system. Runoff response to antecedent conditions and forecasts suggest a highly asymmetric utility function for rainfall forecasts, with greatest decision‐support potential for persistent wet forecasts under wet antecedent conditions when the forecast signal is least dampened by soil‐storage effects. Under average and dry antecedent conditions, rainfall forecasts showed little potential value for practical applications in surface water resources assessments.     

SWMHMS - SMALL WATERSHED MONTHLY HYDROLOGIC MODELING SYSTEM1

ABSTRACT: SWMHMS is a conceptual computer modeling program developed to simulate monthly runoff from a small nonurban watershed. The input needed to run model simulations include daily precipitation, monthly data for evapotranspiration determination (average temperature, crop consumptive coefficients, and percent daylight hours), and six watershed parameter values. Evapotranspiration was calculated with the Blaney-Criddle equation while surface runoff was determined using the Soil Conservation Service curve number procedure. For watershed parameter evaluation, SWMHMS provides options for both optimization and sensitivity analysis. Observed runoff data are required along with the model input previously mentioned in order to conduct parameter optimization. SWMEIMS was tested with data from six watersheds located in different regions of the United States. Model accuracy was generally found to be very good except on watersheds having substantial snowfall accumulation. In having only six watershed parameters, SWMHMS is less complex to use than many other computer programs that calculate monthly runoff. Consequently, SWMHMS may find its greatest application as an educational tool for students learning principles of hydrologic modeling, such as parameter evaluation procedures and the impacts of input data uncertainty on model results.     

Evaluation of the MIKE SHE Model for Application in the Loess Plateau,China1

Abstract: Quantifying the hydrologic responses to land use/land cover change and climate variability is essential for integrated sustainable watershed management in water limited regions such as the Loess Plateau in Northwestern China where an adaptive watershed management approach is being implemented. Traditional empirical modeling approach to quantifying the accumulated hydrologic effects of watershed management is limited due to its complex nature of soil and water conservation practices (e.g., biological, structural, and agricultural measures) in the region. Therefore, the objective of this study was to evaluate the ability of the distributed hydrologic model, MIKE SHE to simulate basin runoff. Streamflow data measured from an overland flow‐dominant watershed (12 km²) in northwestern China were used for model evaluation. Model calibration and validation suggested that the model could capture the dominant runoff process of the small watershed. We found that the physically based model required calibration at appropriate scales and estimated model parameters were influenced by both temporal and spatial scales of input data. We concluded that the model was useful for understanding the rainfall‐runoff mechanisms. However, more measured data with higher temporal resolution are needed to further test the model for regional applications.     

VERIFICATION AND APPLICATION OF REGIONAL EQUATIONS FOR THE STORAGE FUNCTION RUNOFF MODEL1

The storage function model is a nonlinear rainfall-runoff model that has been developed for and applied to flood runoff analysis in Japan. This paper extends the model applicability by developing practical equations for estimating model parameters which are appropriate on a regional basis, i.e., so-called regional equations. Previously, the parameters were computed from historical data for a specific basin or from relationships that do not account for land use and topography. To develop the regionalized equations, model parameters were identified for 91 flood events from 22 watersheds in Japan by applying a mathematical optimization technique. Results from 39 of these events were statistically compared and regional relationships were determined as a function of land use, basin area and rainfall intensity. The utility of the estimated equations were tested by computing runoff hydrographs for lumped basins. The estimated parameters were also applied in a distributed watershed model formulation. Both applications showed acceptable results that validate the use of the regionalized relationships.     

MODEL ACCURACY IN SNOWMELT-RUNOFF FORECASTS EXTENDING FROM 1 TO 20 DAYS1

ABSTRACT: This paper examines the performance of snowmelt-runoff models in conditions approximating real-time forecast situations. These tests are one part of an intercomparison of models recently conducted by the World Meteorological Organization (WMO). Daily runoff from the Canadian snowmelt basin Illecille. waet (1155 km², 509–3150 m a.s.l.) was forecast for 1 to 20 days ahead. The performance of models was better than in a previous WMO project, which dealt with runoff simulations from historical data, for the following reasons: (1) conditions for models were more favorable than a real-time forecast situation because measured input data and not meteorological forecast inputs were distributed to the modelers; (2) the selected test basin was relatively easy to handle and familiar from the previous WMO project; and (3) all kinds of updating were allowed so that some models even improved their accuracy towards longer forecast times. Based on this experience, a more realistic follow-up project can be imagined which would include temperature forecasts and quantitative precipitation forecasts instead of measured data.     

EVALUATION OF INFLOW TO MIRROR LAKE,NEW HAMPSHIRE1

ABSTRACT: Measured stream discharge plus calculated ground water discharge (total measured runoff) were compared with runoff calculated by the unit-runoff method for the two largest watersheds of Mirror Lake for 1981–1983. Runoff calculated by the unit-runoff method, using Hubbard Brook watershed 3 as the index watershed, was greater than the total measured runoff into Mirror Lake during periods of high flow and slightly less than the total measured runoff into Mirror Lake during periods of low flow. Annual calculated unit runoff was 17 to 37 percent greater than total measured runoff. Differences in monthly runoff are far greater, ranging from 0 to greater than 100 percent. For high flows the calculated unit runoff is about 2 times greater than total measured runoff. For low flows the northwest basin of Mirror Lake has the greatest ground water contribution compared to the other two basins. In contrast, Hubbard Brook watershed 3 has the least ground water contribution.     

UNCERTAINTY ANALYSIS OF RUNOFF ESTIMATES FROM A RUNOFF CONTOUR MAP1

ABSTRACT: The U.S. Environmental Protection Agency (EPA) in cooperation with the U.S. Geological Survey (USGS) conducted an analysis to quantify the uncertainty associated with interpolating runoff to specific sites using a runoff contour map. We interpolated runoff to 93 gaged watersheds from a runoff contour map using (1) hand interpolation to the watershed outlet, (2) a computer interpolation to the watershed outlet, and (3) hand interpolation to the watershed centroid. We compared the interpolated values to the actual gaged values and found that there was a bias in the average interpolated value for runoff estimated at basin outlets, with interpolated values being less than the actual. We found no significant difference between the hand interpolation method and the computer interpolation method except that the computer method tended to have higher variability due to factors inherent to the software used. There were no strong spatial correlations or regional patterns in the runoff interpolations, which indicates that there are no regional biases introduced in the development of the contour map. We determined that we could estimate runoff, on the average, within approximately 8.9 cm (3.5 in; 15 percent) of the measured value using the three methods. The results of this work indicate that runoff contour maps can he used in regional studies to estimate runoff to ungaged systems with quantifiable uncertainty.     

STOCHASTIC ANALYSIS OF TIME-AGGREGATED HYDROLOGIC DATA1

Stochastic models fitted to hydrologic data of different time scales are interrelated because the higher time scale data (aggregated data) are derived from those of lower time scale. Relationships between the statistical properties and parameters of models of aggregated data and of original data are examined in this paper. It is also shown that the aggregated data can be more accurately predicted by using a valid model of the original data than by using a valid model of the aggregated data. This property is particularly important in forecasting annual values because only a few annual values are usually available and the resulting forecasts are relatively inaccurate if models based only on annual data are used. The relationships and forecasting equations are developed for general aggregation time and can be used for hourly and daily, daily and monthly or monthly and yearly data. The method is illustrated by using monthly and yearly streamflow data. The results indicate that various statistical characteristics and parameters of the model of annual data can be accurately estimated by using the monthly data and forecasts of annual data by using monthly models have smaller one step ahead mean square error than those obtained by using annual data models.     

GENERALIZED WATERSHED LOADING FUNCTIONS FOR STREAM FLOW NUTRIENTS1

ABSTRACT: Loading functions are proposed as a general model for estimating monthly nitrogen and phosphorus fluxes in stream flow. The functions have a simple mathematical structure, describe a wide range of rural and urban nonpoint sources, and couple surface runoff and ground water discharge. Rural runoff loads are computed from daily runoff and erosion and monthly sediment yield calculations. Urban runoff loads are based on daily nutrient accumulation rates and exponential wash off functions. Ground water discharge is determined by lumped parameter unsaturated and saturated zone soil moisture balances. Default values for model chemical parameters were estimated from literature values. Validation studies over a three-year period for an 850 km² watershed showed that the loading functions explained at least 90 percent of the observed monthly variation in dissolved and total nitrogen and phosphorus fluxes in stream flow. Errors in model predictions of mean monthly fluxes were: dissolved phosphorus - 4 percent; total phosphorus - 2 percent; dissolved nitrogen - 18 percent; and total nitrogen - 28 percent. These results were obtained without model calibration.     

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.     

WATERSHED CONFIGURATION AND GEOGRAPHIC INFORMATION SYSTEM PARAMETERIZATION FOR SPUR MODEL HYDROLOGIC SIMULATIONS1

ABSTRACT: A grid cell geographic information system (GIS) is used to parameterize SPUR, a quasi-physically based surface runoff model in which a watershed is configured as a set of stream segments and contributing areas. GIS analysis techniques produce various watershed configurations by progressive simplification of a stream network delineated from digital elevation models (DEM). We used three watershed configurations: ≥ 2nd, ≥ 4th, and ≥ 13th Shreve order networks, where the watershed contains 28, 15, and 1 channel segments with 66, 37, and 3 contributing areas, respectively. Watershed configuration controls simulated daily and monthly sums of runoff volumes. For the climatic and topographic setting in southeastern Arizona the ≥ 4th order configuration of the stream network and contributing areas produces results that are typically as good as the ≥ 2nd order network. However both are consistently better than the ≥ 13th order configuration. Due to the degree of parameterization in SPUR, model simulations cannot be significantly improved by increasing watershed configuration beyond the ≥ 4th order network. However, a range of Soil Conservation Service curve numbers derived from rainfall/runoff data can affect model simulations. Higher curve numbers yield better results for the ≥ 2nd order network while lower curve numbers yield better results for the ≥ 4th order network.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

SENSITIVITY CONSIDERATIONS WHEN MODELING HYDROLOGIC PROCESSES WITH DIGITAL ELEVATION MODEL1

ABSTRACT: The purpose of this paper is to investigate the sensitivity of a hydrologic models to the type of DEM used. This was done while modeling basin water quality with 1:24,000 and 1:250,000 U.S. Geological Survey DEMs as input to model hydro‐logical processes. The manner in which the model results were sensitive to the choice of raster cell size (scale) is investigated in this study. The Broadhead watershed, located in New Jersey, USA, was chosen as a study area. Curve numbers were estimated by a trial and error to match simulated and observed total discharge. Monthly runoff for the watershed was used in the calibration process. Higher runoff volumes were simulated by the model when the 1:24,000 DEM were used as input data, probably due to the finer resolution which simulated increased average slope and hence higher estimated runoff from the watershed. As the simulated slope of the watershed is flatten with the 1:250,000 DEM, the response of stream flow was delayed and simulated less runoff volume.     

FLOW PERSISTENCE AS A MODIFIER OF SEASONAL RUNOFF PATTERNS1

ABSTRACT: The effect of flow persistence on seasonal patterns of watershed runoff was modeled by using runoff of the immediate antecedent month as an index. Monthly runoff was expressed as a function of monthly rainfall, season of the year, and runoff of the antecedent month. The three independent variables were expressed functionally as sliding polynomials, thus producing a piece-wise, form-free, three-dimensional causative structure. A model form allowing complete interactivity of the three independent variables could not be optimized because of insufficient data with high values of both antecedent runoff and monthly rainfall. A model with reduced interactivity was successfully optimized. Data sets from five watersheds ranging from 0.14 to 398 square miles were analyzed. Results were presented as a series of contour maps that showed contours of monthly runoff in the data space of season and monthly rain. In the series of maps, the patterns of the runoff contours changed with changing values of antecedent runoff. During the wet season of the year the contours changed significantly with antecedent runoff, but changes in the dry season were minimal. The quantitative change of runoff was more readily portrayed with cross-sections through the contoured surfaces.     

AI-based runoff simulation based on remote sensing observations: A case study of two river basins in the United States and Canada

Data-driven techniques are used extensively for hydrologic time-series prediction. We created various data-driven models (DDMs) based on machine learning: long short-term memory (LSTM), support vector regression (SVR), extreme learning machines, and an artificial neural network with backpropagation, to define the optimal approach to predicting streamflow time series in the Carson River (California, USA) and Montmorency (Canada) catchments. The moderate resolution imaging spectroradiometer (MODIS) snow-coverage dataset was applied to improve the streamflow estimate. In addition to the DDMs, the conceptual snowmelt runoff model was applied to simulate and forecast daily streamflow. The four main predictor variables, namely snow-coverage (S-C), precipitation (P), maximum temperature (T_max), and minimum temperature (T_min), and their corresponding values for each river basin, were obtained from National Climatic Data Center and National Snow and Ice Data Center to develop the model. The most relevant predictor variable was chosen using the support vector machine-recursive feature elimination feature selection approach. The results show that incorporating the MODIS snow-coverage dataset improves the models' prediction accuracies in the snowmelt-dominated basin. SVR and LSTM exhibited the best performances (root mean square error = 8.63 and 9.80) using monthly and daily snowmelt time series, respectively. In summary, machine learning is a reliable method to forecast runoff as it can be employed in global climate forecasts that require high-volume data processing.     

A New Open‐Access HUC‐8 Based Downscaled CMIP‐5 Climate Model Forecast Dataset for the Conterminous United States

Watershed‐scale hydrologic simulation models generally require climate data inputs including precipitation and temperature. These climate inputs can be derived from downscaled global climate simulations which have the potential to drive runoff forecasts at the scale of local watersheds. While a simulation designed to drive a local watershed model would ideally be constructed at an appropriate scale, global climate simulations are, by definition, arbitrarily determined large rectangular spatial grids. This paper addresses the technical challenge of making climate simulation model results readily available in the form of downscaled datasets that can be used for watershed scale models. Specifically, we present the development and deployment of a new Coupled Model Intercomparison Project phase 5 (CMIP5) based database which has been prepared through a scaling and weighted averaging process for use at the level of U.S. Geological Survey (USGS) Hydrologic Unit Code (HUC)‐8 watersheds. The resulting dataset includes 2,106 virtual observation sites (watershed centroids) each with 698 associated time series datasets representing average monthly temperature and precipitation between 1950 and 2099 based on 234 unique climate model simulations. The new dataset is deployed on a HydroServer and distributed using WaterOneFlow web services in the WaterML format. These methods can be adapted for downscaled General Circulation Model (GCM) results for specific drainage areas smaller than HUC‐8. Two example use cases for the dataset also are presented.     

APPLICATION OF THE SNOWMELT RUNOFF MODEL USING MULTIPLE-PARAMETER LANDSCAPE ZONES ON THE TOWANDA CREEK BASIN,PENNSYLVANIA1

ABSTRACT: The Snowmelt Runoff Model (SRM) is designed to compute daily stream discharge using satellite snow cover data for a basin divided into elevation zones. For the Towanda Creek basin, a Pennsylvania watershed with relatively little relief, analysis of snow cover images revealed that both elevation and land use affected snow accumulation and melt on the landscape. The distribution of slope and aspect on the watershed was also considered; however, these landscape features were not well correlated with the available snow cover data. SRM streamflow predictions for 1990, 1993 and 1994 snowmelt seasons for the Towanda Creek basin using a combination of elevation and land use zones yielded more precise streamflow estimates than the use of standard elevation zones alone. The use of multiple-parameter zones worked best in non-rain-on-snow conditions such as in 1990 and 1994 seasons where melt was primarily driven by differences in solar radiation. For seasons with major rain-on-snow events such as 1993, only modest improvements were shown since melt was dominated by rainfall energy inputs, condensation and sensible heat convection. Availability of GIS coverages containing satellite snow cover data and other landscape attributes should permit similar reformulation of multiple-parameter watershed zones and improved SRM streamflow predictions on other basins.     

NETWORK DESIGN FOR WATER SUPPLY FORECASTING IN THE WEST1

ABSTRACT: The objective is to develop techniques to evaluate how changes in basic data networks can improve accuracy of water supply forecasts for mountainous areas. The approach used was to first quantify how additional data would improve our knowledge of winter precipitation, and second to estimate how this knowledge translates, quantitatively, into improvement in forecast accuracy. A software system called DATANET was developed to analyze each specific gage network alternative. This system sets up a fine mesh of grid points over the basin. The long-term winter mean precipitation at each grid point is estimated using a simple atmospheric model of the orographic precipitation process. The mean runoff at each grid point is computed from the long-term mean precipitation estimate. The basic runoff model is calibrated to produce the observed long-term runoff. The error analysis is accomplished by comparing the error in forecasts based on the best possible estimate of precipitation using all available data with the error in the forecasts based on the best possible estimate of winter precipitation using only the gaged data. Different data network configurations of gage sites can be compared in terms of forecast errors.     

ASSESSING FLOOD MITIGATION ALTERNATIVES IN SHIJR AREA IN METROPOLITAN TAIPEI1

The Keelung River Basin in northern Taiwan lies immediately upstream of the Taipei metropolitan area. The Shijr area is in the lower basin and is subject to frequent flooding. This work applies micromanagement and source control, including widely distributed infiltration and detention/ retention runoff retarding measures, in the Wudu watershed above Shijr. A method is also developed that combines a genetic algorithm and a rainfall runoff model to optimize the spatial distribution of runoff retarding facilities. Downstream of Wudu in the Shijr area, five dredging schemes are considered. If 10‐year flood flows cannot be confined in the channel, then a levee embankment that corresponds to the respective runoff retarding scheme will be required. The minimum total cost is considered in the rule to select from the regional flood mitigation alternatives. The results of this study reveal that runoff retarding facilities installed in the upper and middle parts of the watershed are most effective in reducing the flood peak. Moreover, as the cost of acquiring land for the levee embankment increases, installing runoff retarding measures in the upper portion of the watershed becomes more economical.     

Developing Subseasonal to Seasonal Climate Forecast Products for Hydrology and Water Management

We describe a new effort to enhance climate forecast relevance and usability through the development of a system for evaluating and displaying real‐time subseasonal to seasonal (S2S) climate forecasts on a watershed scale. Water managers may not use climate forecasts to their full potential due to perceived low skill, mismatched spatial and temporal resolutions, or lack of knowledge or tools to ingest data. Most forecasts are disseminated as large‐domain maps or gridded datasets and may be systematically biased relative to watershed climatologies. Forecasts presented on a watershed scale allow water managers to view forecasts for their specific basins, thereby increasing the usability and relevance of climate forecasts. This paper describes the formulation of S2S climate forecast products based on the Climate Forecast System version 2 (CFSv2) and the North American Multi‐Model Ensemble (NMME). Forecast products include bi‐weekly CFSv2 forecasts, and monthly and seasonal NMME forecasts. Precipitation and temperature forecasts are aggregated spatially to a United States Geological Survey (USGS) hydrologic unit code 4 (HUC‐4) watershed scale. Forecast verification reveals appreciable skill in the first two bi‐weekly periods (Weeks 1–2 and 2–3) from CFSv2, and usable skill in NMME Month 1 forecast with varying skills at longer lead times dependent on the season. Application of a bias‐correction technique (quantile mapping) eliminates forecast bias in the CFSv2 reforecasts, without adding significantly to correlation skill.