Hydrologic and Water Quality Integration Tool: HydroWAMIT

A spatially distributed and continuous hydrologic model focusing on total maximum daily load (TMDL) projects was developed. Hydrologic models frequently used for TMDLs such as the hydrologic simulation program—FORTRAN (HSPF), soil and water assessment tool (SWAT), and generalized watershed loading function (GWLF) differ considerably in terms of spatial resolution, simulated processes, and linkage flexibility to external water quality models. The requirement of using an external water quality model for simulating specific processes is not uncommon. In addition, the scale of the watershed and water quality modeling, and the need for a robust and cost-effective modeling framework justify the development of alternative watershed modeling tools for TMDLs. The hydrologic and water quality integration tool (HydroWAMIT) is a spatially distributed and continuous time model that incorporates some of the features of GWLF and HSPF to provide a robust modeling structure for TMDL projects. HydroWAMIT operates within the WAMIT structure, developed by Omni Environmental LLC for the Passaic River TMDL in N. J. HydroWAMIT is divided into some basic components: the hydrologic component, responsible for the simulation of surface flow and baseflow from subwatersheds; the nonpoint-source (NPS) component, responsible for the calculation of the subwatershed NPS loads; and the linkage component, responsible for linking the flows and loads from HydroWAMIT to the water quality analysis simulation program (WASP). HydroWAMIT operates with the diffusion analogy flow model for flow routing. HydroWAMIT provides surface runoff, baseflow and associated loads as outputs for a daily timestep, and is relatively easy to calibrate compared to hydrologic models like HSPF. HydroWAMIT assumes that the soil profile is divided into saturated and unsaturated layers. The water available in the unsaturated layer directly affects the surface runoff from pervious areas. Surface runoff from impervious areas is calculated separately according to precipitation and the impervious fractions of the watershed. Baseflow is given by a linear function of the available water in the saturated zone. The utility of HydroWAMIT is illustrated for the North Branch and South Branch Raritan River Watershed (NSBRW) in New Jersey. The model was calibrated, validated, and linked to the WASP. The NPS component was tested for total dissolved solids. Available weather data and point-source discharges were used to prepare the meteorological and flow inputs for the model. Digital land use, soil type datasets, and digital elevation models were used for determining input data parameters and model segmentation. HydroWAMIT was successfully calibrated and validated for monthly and daily flows for the NSBRW outlet. The model statistics obtained using HydroWAMIT are comparable with statistics of HSPF and SWAT applications for medium and large drainage areas. The results show that HydroWAMIT is a feasible alternative to HSPF and SWAT, especially for large-scale TMDLs that require particular processes for water quality simulation and minor hydrologic model calibration effort.     

Evaluating Watershed Experiments through Recursive Residual Analysis

The question of watershed response to changes in land use and land cover has received a great deal of attention in the hydrologic literature. One of the primary tools for quantifying such changes has been the paired watershed approach. In general this approach applies a linear regression analysis to relate water yield from a control watershed to a physiographically similar watershed in close proximity, and to predict changes in yield following treatment. Although much research on paired watershed experiments has focused on quantification of the impacts of land-use and land cover changes on water yield, little attention has been brought to evaluating explicitly the methodology used to quantify such effects. An alternative method is proposed for examining treatment impacts and their duration, through the application of a cumulative recursive residual test of model stability. The results from a paired watershed study of ponderosa pine watersheds in north-central Arizona are analyzed using both a traditional linear regression analysis and a recursive residuals approach. The results suggest that the linear regression approach may fail to account for the true complexity of watershed response to vegetation treatments, by underestimating the interactions of treatment impacts, shifts in climatic drivers, and revegetation rates.     

Approach Using Active Groundwater Storage for Hydrologic Model Calibration in West-Central Florida

Hydrologic model calibration is always a challenging and tedious process especially for the calibration of complex models, which includes continuous hydrograph models, requires sophisticated calibration methods. The Hydrologic Simulation Program-FORTRAN (HSPF) is one of the popular and powerful time variable hydrologic models. However, in order to improve the assessment of hydrologic activities in shallow ground water settings, the model needs to be reliably calibrated for ground water contribution. Little guidance is provided in the literature concerning the manner of this contribution. In fact, the most common calibration of HSPF uses subjective parameter fitting and focuses on the attainment of statistical goodness of fit of runoff fluxes and water levels, ignoring ground water components. The goal of this research is using a different approach to calibrate HSPF with observed water table records. In this study, HSPF is applied on a small area in west-central Florida and calibrated by comparing active ground water storage to well elevation records in range land and forested land covers. The Nash-Sutcliffe efficiency and correlation coefficient computed using observed and simulated daily flows are 0.91 and 0.96 at Peace River, respectively, also with good fair results for other stations in the model domain. The study shows that improved calibration of the model can be achieved if active ground water storage and well records are compared for timing and magnitude of fluctuations.     

Real-Time Watershed Delineation System Using Web-GIS

A Web-based geographic information system (GIS) that can delineate watersheds in real time was developed for use in hydrologic analysis and to support hydrologic model operation on the Internet. The system integrates a watershed delineation (WD) engine, common gateway interfaces, a Web-GIS user interface, and supports digital spatial data including vector and grid formats. The WD program utilizes a double-seed array-replacement algorithm to obtain a watershed boundary from point coordinates. The system provides a user interface for the selection of an outlet point from a map display in the Web browser using MapServer Web-GIS capability. The WD and data extraction system has been implemented for all of Indiana, with extensive verification conducted in the 2,082.7 km2 Wildcat Creek watershed in Indiana. The time to obtain results and the quality of results are acceptable for use as a real-time system for WD via the Web.     

Modeling Nonpoint Source Pollutant Losses from a Small Watershed Using HSPF Model

Hydrologic models play an important role in the assessment of nonpoint source (NPS) pollution, which is essential for the environmental management of water resources. The present study has been undertaken to evaluate the applicability of a physically based continuous time scale, hydrological, and water quality computer model—Hydrologic Simulation Program-Fortran (HSPF)—in simulating runoff and sediment associated NPS pollutant losses from a small mixed type (land under agriculture, shrubs and forest, rocks, grasses) watershed of the Damodar Valley Corporation, Hazaribagh, India. Water soluble NO3–N, NH4–N, and P were considered as pollutants and their concentrations in the runoff were measured at the outlet of the watershed, randomly for 15 dates during the monsoon season (June–October) of 2000 and 2001. The model calibration and validation results reveal that the seasonal trend of HSPF simulated runoff, sediment yield, and NPS pollutants compared reasonably with their measured counterparts. Although the concentrations of pollutants were generally overpredicted for NO3–N and underpredicted for NH4–N and water-soluble P in the month of June when fertilizers releasing NH4–N and P are applied in rice fields, the differences in the mean concentration were not significantly different at a 95% level of confidence. Variation in the simulated losses of water soluble N and P species between the years occurred largely due to differences in the amount and distribution of rainfall. These results indicate that the HSPF model can be used as a tool for simulating runoff and sediment associated NPS pollution losses from the study area.     

Development of an Empirical Lag Time Equation

A regression analysis was performed on measured lag times from gauged watersheds to develop a lag time equation. The watersheds are part of the Agricultural Research Service’s database. They are located in several states and are comprised of varying terrain. The goal of the analysis was to develop a lag time equation that is useful in hydrologic modeling. The study included measurements from approximately 10,000 direct runoff events from 52 watersheds to determine which watershed parameters are best for predicting lag time. The lag time was found to correlate strongly with the longest hydraulic length of the watershed. Therefore an equation was developed that used only this parameter. The inclusion of any other watershed characteristics in the equation did not improve its ability to predict the lag time. Finally, the National Resource Conservation Service procedures for calculating watershed lag time were used to determine the lag times of the watersheds. These estimated lag times were then compared with the measured lag time of the watershed. It was found that the use of these methods generally underpredicted the true lag time of a watershed.     

Optimization of Watershed Control Strategies for Reservoir Eutrophication Management

A vital key to the development of a reservoir eutrophication management strategy is to link the watershed-nutrient model to the model of reservoir water quality. To develop a cost-effective optimization model, a coupled watershed-reservoir model with an optimization model has been developed to design control strategies in the watershed in a planning time horizon. This methodology can help reduce the phosphorus concentration of a reservoir to the standard level. In this study, the weather data for the next 10 years was generated using downscaled GCM data to simulate the watershed phosphorus load using the SWAT model. Then an optimal model for selection and placement of best management practices (BMP) at watershed scale is developed by linking the coupled watershed and reservoir models with a genetic algorithm. This model is able to identify the minimum present cost design (type and location) of BMP structural alternatives. The objective of water quality is obtained using a system dynamic model for reservoir phosphorus concentration to determine a permissible phosphorus load as the main agent of eutrophication in a reservoir. Structural BMPs in this study include, filter strips, parallel terraces, grade stabilization structures, and detention ponds. The optimum solution was obtained through a trade-off curve between cost and exceedance magnitude from the standard of reservoir phosphorus concentration. The case study is the Aharchai River Watershed upstream of the Satarkhan Reservoir in the northwestern part of Iran.     

Simulating the Effect of Sulfate Addition on Methylmercury Output from a Wetland

The watershed analysis risk management framework (WARMF) model was applied to Wetland S6 of the Marcell Experimental Forest, using the data from a field experiment, conducted to investigate the effect of sulfate additions on mercury methylation in the wetland. The wetland was modeled as interconnected land catchments. Actual meteorology data and mercury and sulfate concentrations of precipitation were input to the model. To simulate the sulfate sprinkling, the experimental section of the bog was irrigated with sulfate water on the actual dates of sulfate additions. The model simulated wetland outflows that matched the measured outflows with an R-square of 0.856. WARMF also simulated other phenomena observed in the experiment: higher sulfate and MeHg levels at the wetland outlet after every sulfate addition, and higher sulfate and MeHg levels in the pore water of the bog after only the May addition, not the July and September additions. According to WARMF, the low groundwater table in May allowed the sprinkled sulfate to percolate to the soil stratum 10–30 cm below the ground level of the bog, where the pore water was sampled. In July and September, the sulfate could not reach that zone because the percolation was blocked by high groundwater tables. The sampled soil stratum was not the site of methylation that contributed MeHg to the wetland outlet. The saturated zone of the top 10 cm of bog was the site that produced MeHg, which was flushed to the outlet after all sulfate additions. WARMF predicted that quadrupling the sulfate deposition would increase the MeHg output by 216%, which might become lower with more data and better model calibration.     

Predictive Modeling of Storm-Water Runoff Quantity and Quality for a Large Urban Watershed

A predictive model for storm-water runoff was implemented on a GIS platform based on the unit area loading method and Browne’s empirical relation for soil characteristics for the Upper Ballona Creek Watershed in Los Angeles. The heterogeneity of the watershed was quantified by dividing it into many small subareas and applying lumped parameters for each. Characterization of total pollutant load by land-use types to total loads was achieved through zeroth-order regularization and limited memory Broyden–Fletcher–Goldfarb–Shanno bound constrained optimization techniques. Relative form was used in the objective function to compensate for strong contributions of high magnitude variables. Model predictions showed reasonable agreement with pollutant loadings, using Zn as an example, measured at the mass emission site at watershed mouth. The predicted runoff volumes using the developed quantity model were in good agreement with the data and had R2 of 0.86. The RMS error of the quality model was 9?kg, which is low compared to the mean discharge of 77?kg/event.     

Variations of Time of Concentration Estimates Using NRCS Velocity Method

Time of concentration (Tc) is the time required for runoff to travel from the hydraulically most distant point to the outlet of a watershed. The Natural Resources Conservation Service (NRCS) velocity method commonly is used to estimate Tc for hydrologic analysis and design. The NRCS velocity method applies the physical concept that travel time is a function of runoff flow length and flow velocity. Time of concentration for 96 Texas watersheds is independently estimated by three research teams using the NRCS velocity method. Drainage areas of the 96 watersheds considered in the study are approximately 0.8–440.3?km2 (0.3–170?mi2). Digital elevation models having a grid size of 30?m were used to derive watershed physical characteristics using ArcGIS or HEC-GeoHMS. Average channel width was estimated from 1?m or 1?ft digital orthoimagery quarter quadrangle or aerial photography. Each team made independent decisions to estimate parameters needed for different flow segments for the NRCS velocity method. Estimates of time of concentration made by three research teams are compared, and both graphic comparison and statistical summary demonstrate that time of concentration estimated using the NRCS velocity method is subject to large variation, dependent on the analyst-derived parameters used to estimate flow velocity. Because of the propensity for different analysts to arrive at different results, caution is required in application of the NRCS velocity method to estimate Tc.     

Rational Coefficients for Steeply Sloped Watersheds

When computing peak discharges for the design of drainage systems using the rational method, it is important to have an accurate value for the rational coefficient (C). For steeply sloped watersheds the origin of values of the rational coefficient are unknown and lack even modeling verification. A model that shows the relationship between the rational coefficient and watershed slope was developed for steeply sloped watersheds. Using Horton’s infiltration equation, Manning’s equation, the velocity method for computing times of concentration, and generalized intensity-duration-frequency curves, a model was developed to test the effect of variation of several watershed characteristics on the relationship between slope and the rational coefficient. Analyses with the model showed that both Manning’s coefficient and land use had the greatest effect on the relationship between C and slope. A mathematical function was then developed from data generated from the Horton–Manning model. This model allows C to be estimated for a given slope and a value of Manning’s coefficient for the land cover. A rational coefficient at a 6% slope is also required input. The model was tested using several watersheds with moderate to steep slopes. This relationship should be used to better estimate values of C on steep slopes, and thereby, lead to more accurately hydrologic designs.     

Upland Erosion Modeling in a Semihumid Environment via the Water Erosion Prediction Project Model

The major water quality impairment in the midwest United States is sediment eroded from agricultural lands. Yet, few understand the spatial and temporal variability of erosion, or soil erosion dynamics, in relation to precipitation, topography, land management, and severe events. The objectives of this paper are to (1) develop a methodology for estimating long-term spatial soil erosion and water runoff losses and (2) explore issues in applying an established physical-based process model, Water Erosion Prediction Project (WEPP), to a large area by establishing a prototype system for the state of Iowa. This study for the first time provides a comparison of the model predictions against long-term measurements of the sediment delivery ratio (SDR) in the South Amana Catchment of the Clear Creek Watershed (CCW), a heavily instrumented watershed that is roughly 10 times the maximum WEPP fold size. To further examine the performance of WEPP in a semihumid environment, such as the CCW, where runoff and raindrop impact to erosion may be significant, the SDR was plotted as a function of the runoff coefficient, defined as the runoff/rainfall ratio. In addition, the WEPP predictions are compared against the statistical relation of SDR vs. runoff coefficient developed by Piest et al. in 1975) for watersheds in Iowa. It is shown that WEPP follows the trend shown by Piest et al. quite closely and performs well for continuous simulations extended up to 300?years.     

Case Study: Evaluation of Streamflow Partitioning Methods

Understanding water flow and its relative quantities through different pathways is vital for watershed management. Like many problems in hydrology, numbers of methods have been proposed for streamflow partitioning. Five methods were identified as being the most relevant and least input intensive. This study tested performance of these methods against separately measured surface and subsurface flow data from the coastal plain physiographic region of the southeastern United States. Separately measured surface and subsurface flow were collected for 12 years (1970–1981) in a field scale watershed by the Southeast Watershed Research Laboratory of the USDA-Agricultural Research Service. Results of comparative analysis indicated that Method IV performed the best. Results also indicated that accuracy of this method is highly dependent upon the proper estimation of the “fraction coefficient” that is based on many physical and hydrologic characteristics of the watershed. This study concluded that deterministic/empirical methods such as Boughton’s Method IV, require proper parameter value for increased accuracy.     

Modeling of Vegetation-Erosion Dynamics in Watershed Systems

Vegetation and erosion are a pair of competing and interactive factors that affect the quality of watershed ecosystems. The objective of this study is to develop an innovative approach for conceptualizing and simulating the vegetation-erosion dynamics. Differential equations of vegetation-erosion dynamics have been developed to describe the relevant vegetation processes, with the relevant solution methods being provided. Based on the developed model, a vegetation-erosion chart can be produced for predicting the tendencies of vegetation and erosion under different land-use conditions. Thus decision supports in terms of desired measures to improve the system conditions can be provided. In general, vegetation of a watershed may exist in three states, including (1) vegetation-developing and erosion-reducing; (2) vegetation-deteriorating and erosion-increasing; and (3) transitional state between states (1) and (2). Humans may change a watershed system from one state into another. The effort needed for such a change depends on the distance between the present position and the destination one as shown on the vegetation-erosion chart. The developed model has been applied to three regions, including the Xiaojiang, Heishui, and Shengou Watersheds in China. The results demonstrate that the proposed vegetation-erosion dynamics is a powerful tool for simulating and predicting vegetation evolutions in the watersheds. Generally, reforestation and erosion-control measures would improve vegetation coverage slowly in the first 10 years, but become much faster in the second 10 years; this implies that a long-term strategy is needed. The results also indicate that, for revegetating hilly areas, erosion control is critical; merely planting trees and shrubs is insufficient for greening the exposed land.     

Environmental Index for Assessing Spatial Bias in Watershed Sampling Networks

Bias in the design of stream-sampling networks can be a major cause of inaccurate characterization of ambient water quality. At the state level, sampling bias can impact a state’s ability to produce an accurate assessment of the water quality of all state waters. At the federal level, this bias has hindered attempts by the U.S. Environmental Protection Agency (U.S. EPA) to produce nationwide assessments of water quality. Three types of bias commonly occur in water-quality assessments: design, analytical, and statistical. This paper focuses on “design bias,” especially spatial-design bias, in stream monitoring networks. A geographic information system method is described to develop an environmental index to help recognize spatial design bias, and to prioritize areas (subwatersheds) for sampling. The environmental index, which is developed using data that are generally available from federal and state agencies, provides a means of differentiating the component parts of a watershed, its subwatersheds, in terms of two sets of features: natural landscape features and anthropogenic features or “stressors.” Together, these features largely determine the variability of the quantity and quality of water discharged from watersheds. Subwatersheds with higher environmental-index values are expected to have more variable water quality over the course of a year than subwatersheds with lower index values. They are also expected to exert a greater influence on basinwide water quality. These properties of the index make it useful for (1) identifying the possible presence of spatial bias in existing watershed-sampling networks, (2) making sampling design decisions, and (3) aggregating data from subwatersheds into basinwide measures of water quality. Development and use of the index is demonstrated in a watershed in southern New Hampshire.     

Integrated Storm-Water Management for Watershed Sustainability

One aspect of integrated watershed management evaluates the impact of development on the local hydrologic cycle and, in particular, drinking water, wastewater, and storm-water infrastructure. Sustainable storm-water management focuses on selecting storm-water controls based on an understanding of the problems in local receiving waters that result from runoff discharges. For example, long-term problems associated with accumulations of pollutants in water bodies include sedimentation in conveyance systems and receiving waters, nuisance algal growths, inedible fish, undrinkable water, and shifts to less sensitive aquatic organisms. Short-term problems associated with high pollutant concentrations or frequent high flows (event-related) include swimming beach closures, water quality violations, property damage from increased flooding, and habitat destruction. A wide variety of individual storm-water controls usually must be combined to form a comprehensive wet weather management strategy. Unfortunately, combinations of controls are difficult to analyze. This will require new modeling techniques that can effectively evaluate a wide variety of control practices and land uses, while at the same time ensure that the flood-control objectives also are met. The results of these new models and novel techniques used for storm-water control then can be incorporated into an evaluation of the urban water cycle for a specific service area to determine whether storm-water controls can provide additional benefits such as reduction of potable water use and reduction of sanitary sewer overflow events.     

Retrofit Storm Water Retention Volume for Low Impact Development

The concept of low impact development (LID) applies decentralized on-site runoff source control to storm water management. LID is an integration of bioretention and vegetated landscapes to maintain a catchment’s hydrologic and ecological functions. In current practice, the LID implementation is quantified for the specified watershed development. During the dynamic development process, the existing LID facilities have to be improved according to the incremental changes in the watershed. This technical note presents an on-site hydrologic approach to relate the required incremental storm water retention volume to the alteration of surface imperviousness in the tributary area. This approach allows the storm water retention volume to be tailored according to the stage of the watershed development. Cumulatively, the total storage volume can be achieved though multiple stages of the watershed development. The incremental retention volume is found to be related to the local average event rainfall depth. Design charts were produced and normalized by the local average rainfall event depth for generalized applicability.     

Integrated Hydrologic Modeling and GIS in Water Resources Management

The integration of a physically-based distributed model with a geographic information system (GIS) in watershed-based water resources management is presented, and an example watershed is chosen to demonstrate the spatial database and modeling system developed in this study. The spatial data is first processed by GIS. The model is then used to simulate runoff hydrographs. It operates at a daily time step on 1 × 1 km grid squares and simulates important hydrologic processes including evapotranspiration, snowmelt, infiltration, aquifer recharge, ground-water flow, and overland and channel runoff. Finally, the model result is displayed by using GIS. This study demonstrates that the integration of a physically based distributed model and GIS may successfully and efficiently implement the watershed-based water resources management. Not only does this process facilitate examination of a wider range of alternatives that would be impossible by using conventional methods, but it also provides a living management that could be modified and updated by water managers once the watershed condition is changed.     

Event and Continuous Hydrologic Modeling with HEC-HMS

Event hydrologic modeling reveals how a basin responds to an individual rainfall event (e.g., quantity of surface runoff, peak, timing of the peak, detention). In contrast, continuous hydrologic modeling synthesizes hydrologic processes and phenomena (i.e., synthetic responses of the basin to a number of rain events and their cumulative effects) over a longer time period that includes both wet and dry conditions. Thus, fine-scale event hydrologic modeling is particularly useful for understanding detailed hydrologic processes and identifying the relevant parameters that can be further used for coarse-scale continuous modeling, especially when long-term intensive monitoring data are not available or the data are incomplete. Joint event and continuous hydrologic modeling with the Hydrologic Engineering Center’s Hydrologic Modeling System (HEC-HMS) is discussed in this technical note and an application to the Mona Lake watershed in west Michigan is presented. Specifically, four rainfall events were selected for calibrating/verifying the event model and identifying model parameters. The calibrated parameters were then used in the continuous hydrologic model. The Soil Conservation Service curve number and soil moisture accounting methods in HEC-HMS were used for simulating surface runoff in the event and continuous models, respectively, and the relationship between the two rainfall-runoff models was analyzed. The simulations provided hydrologic details about quantity, variability, and sources of runoff in the watershed. The model output suggests that the fine-scale (5?min time step) event hydrologic modeling, supported by intensive field data, is useful for improving the coarse-scale (hourly time step) continuous modeling by providing more accurate and well-calibrated parameters.