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.     

Modeling Escherichia Coli and Its Sources in an Urban Bayou with Hydrologic Simulation Program—FORTRAN

Bacterial levels in Buffalo Bayou in Houston commonly exceed contact recreation standards. Potential sources of bacteria include wastewater treatment plants, sanitary sewer overflows, septic systems, wet and dry nonpoint-source discharges via direct runoff and pipes, direct deposition, and sediment. A water-quality model in the Hydrologic Simulation Program—FORTRAN (HSPF) was calibrated and validated for hydrology, sediment, and Escherichia coli and subsequently used to evaluate the impacts of the bacterial sources in the watershed. In addition, simple estimates of bacterial loads were calculated along with source evaluations from load duration curves. Load reductions based upon the simple estimates indicated that water-quality standards were met by reducing dry-weather indicator bacterial loads by 69% and wet-weather loads by 98%. When these load reductions were implemented in the HSPF model, however, standards were not met under dry-weather conditions. Residual nonpoint-source loading was found to cause the discrepancy between simple load estimate calculations and the developed water-quality model. This paper demonstrates that runoff can play a significant role in maintaining high levels of bacteria under all flow conditions and that understanding the temporal variations in bacterial source loading is critical to ensure that load reductions will achieve water-quality standards.     

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.     

Retention of Sediment and Nutrient Loads with Peak Runoff Control

The purpose of this research was to evaluate peak runoff control as a water protection method to reduce sediment and nutrient loads. Increased eutrophication of surface waters and risk of floods demands cost effective methods to reduce pollutant input and risks of flooding. With the peak runoff control it is possible to cut the main peaks and store the runoff water temporarily in ditches. The method decreases the suspended solids (SS) and nutrient loads by reducing flow velocities, and improving the settling of sediment particles. The method was tested in two heavily drained adjacent peat harvesting areas suffering considerable erosion. The peak flows were cut by 27–87%, the SS load by 61–94%, the total nitrogen (Ntot) by 45–91%, and the total phosphorus (Ptot) load by 47–88%. The peak runoff control method operated most effectively during extreme events when most of the SS load is transported. A detailed particle analysis of runoff water showed that water detention reduced the median particle size of SS load as the largest particles settle. The results clearly indicate that the peak runoff control is an effective method to control the sediment loads and peak flows from peatland drainage.     

Regional Regression Models of Annual Streamflow for the United States

Estimates of annual streamflow volumes are needed in many different types of hydrologic studies. Usually a streamgauge is unavailable at the location of interest, hence regional methods that relate streamflow to readily measured geomorphic and climate characteristics provide a practical solution. Hydrologic, geomorphic, and climatic characteristics of 1,553 undeveloped watersheds across the United States are used to develop regional regression equations that relate the first two moments of annual streamflow to readily measured basin and climate characteristics. These relations are summarized for each of 18 major U.S. water resource regions. The relationships are remarkably precise, with adjusted R2 values ranging from 90.2–99.8% and an average value of 96.2% across the continent. The usefulness of these relationships is evaluated by deriving their information content in terms of equivalent record length. These results indicate that regional models of annual streamflow, including runoff maps, are less accurate than suggested by traditional goodness-of-fit statistics. We also provide estimates of precipitation and temperature elasticity of streamflow, by region.     

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.     

Retrospective Comparison of Watershed Analysis Risk Management Framework and Hydrologic Simulation Program Fortran Applications to Mica Creek Watershed

As part of an ongoing watershed model comparison program for forested watersheds, Watershed Analysis Risk Management Framework (WARMF V5.18) and Hydrologic Simulation Program Fortran (HSPF V10) were independently applied to the Mica Creek Watershed in Idaho. A comprehensive model comparison was made in terms of watershed delineation, hydrologic formulations, model parameterization, meteorological data, hydrologic calibration, and hydrologic verification. Comparison was not made for water quality, which was not simulated in the HSPF application. It was concluded that WARMF is a mechanistic model structured to simulate the hydrologic processes, whereas HSPF is an empirical water budget model. The WARMF is suitable for application to forested watersheds. It successfully predicted stream flows comparable to measured values. The HSPF results were also good, if one ignores an unrealistic amount of water loss to inactive groundwater and an empirical treatment of rain-on-snow events.     

Improvement of Regression Simulation in Fluvial Sediment Loads

In the mid-Atlantic region of the United States, sediment loads from stream runoff generally change more rapidly in the rising limb than in the falling limb of a storm hydrograph. As a result, sediment load reaches its peak prior to flow peak, an observation known as clockwise hysteresis. This dynamic load–flow relationship is poorly reproduced by the existing multivariate linear regression models. This paper explores regressors that attempt to incorporate observed features in a statistical model and thus improve load estimates. These included inverse discharge and flow-change regressors. The load estimates using three regression models for eight rivers are compared, and recommended regression equations are proposed.     

Impact of Turbidity Currents on Reservoir Sedimentation

All lakes created on natural rivers are subjected to reservoir sedimentation. The construction of a dam significantly modifies the flow conditions of natural streams inside and downstream of an artificial lake. The sediment concentration is often high during the flood season, and the entering flow shows a greater density than the ambient fluid. Suspended load can therefore be carried along the reservoir bottom to the dam in the form of turbidity currents. This paper presents research results that help to better understand this physical phenomenon, which contributes to reservoir sedimentation. It is based on in situ measurements, a laboratory scale model of turbidity currents and numerical flow simulations. The study of a thousand-year flood in the Luzzone Reservoir in the Swiss Alps using the developed computer model revealed the potential of such a tool. In particular, the impact on the sediment deposits was analyzed. A valuable evaluation of the incidence of such a turbidity flow is presented and its effects are compared to observations. Significant progress has been made in understanding the importance of turbidity currents in reservoir sedimentation.     

Effect of non-point source pollution on water quality of the Weihe River

Non-point source pollution is an important factor that affects the water quality of the Weihe River. To study the effect of the pollution on the water quality of the Weihe River, five flood events and three normal discharge events during non-flood period were monitored from July to December in 2006, at the Lin-tong section. In order to identify how sediment influenced the water quality of the river, raw and supernatant samples taken from the monitored events were analyzed. Supernatant samples were siphoned from 5 cm below the water surface of the raw samples at the lab afterwards the raw samples were shaken up and laid for 30 minutes in the beakers. The results indicated that: 1) The concentrations of SS, COD, TP and TN in the raw samples from flood events were higher than those from normal discharge events. The higher values of the COD, TP and TN in the raw samples resulted from natural humic matters in surface soil. 2) The load transport rate of analyzed water quality indexes increased gradually to its maximum and then drop down, matching that in the discharge hydrograph. The concentrations of SS, NH3-N, NO2-N, NO3-N, COD, and TP in the raw samples increased initially and then decreased, while the concentration of dissolved orthophosphate and TN in the raw samples decreased gradually and then increased. The peak time of concentration and load transport rate of SS as well as COD, TP and TN in the raw samples were close to or lagged behind the time of peak flow.Generally, the peak time of concentration and load transport rate of dissolved orthophosphate,dissolved total-phosphate, NH3-N, NO2-N, and NO3-N occurred prior to the time of peak flow. 3) The mean concentration method was used to calculate the NPS pollution load at Lin-tong section, and the results are credible. In 2006, the proportions of the NPS pollution load to the total load for COD, TP,TN and inorganic nitrogen were more than 30%.     

Assessment of Agricultural NonPoint Source Model for a Watershed in Tropical Environment

Very little work on the application of watershed modeling has been done in the tropical climatic conditions of Thailand to explore the nature of environmental problems arising from nonpoint source pollution due to agricultural activities, and to evaluate possible remedial measures and strategies. The present study attempts to verify the suitability of a nonpoint source pollution model, the Agricultural NonPoint Source model, for the Huai Nong Prong watershed in Southeastern Thailand. Extensive fieldwork was carried out to collect data and information needed for the model preparation and application. The study has revealed that simulated runoff volume, sediment, and nutrient yield from the watershed with mixed land use and relatively high slopes match favorably with observed data. For the ten rainfall events simulated, the coefficient of performance, a measure of model efficiency (equal to zero for a perfect match), was 0.09, 0.47, 0.09, and 0.03 for runoff volume, sediment yield, total nitrogen, and total phosphorus, respectively. The model, however, could not accurately simulate peak flow rates, suggesting the need for changes in the modeling approach or governing equations and relationships to calculate peak discharges in a tropical environment.     

Pesticide Runoff Loads from Lawns and Golf Courses

Pesticide runoff loads from grass surfaces were estimated for 29 chemicals commonly applied to U.S. lawns and golf courses. Data on pesticide properties and typical application rates and schedules were developed and summarized as input parameters for the TurfPQ runoff model. Weather data for each of 9 U.S. cities were generated by the USCLIMATE model and modified by the addition of growing season irrigation. Simulation runs were made for each chemical, grass surface (lawns, greens, fairways), and city, and the results were summarized as mean annual and 1-in-10 year annual maximum daily pesticide loads. These loads varied greatly with pesticide, grass surface, and city, ranging from 0 to 875?g/ha for mean loads and 0 to 818?g/ha for 1-in-10 year daily loads. Mean annual loads averaged over the 29 chemicals and 3 grass surfaces were found to be closely related to growing season precipitation. Variations among the nine cities were well-captured by three general climate categories: Humid, characterized by abundant precipitation and warm temperatures, represented by Atlanta and Houston; temperate, with moderate precipitation and temperature, as with Albany, Columbus, Madison, and Olympia; and dry, with sparse precipitation, represented by Bismarck, Fresno, and Roswell. Mean annual pesticide runoff was 37, 9, and 2?g/ha in the humid, temperate, and dry regions, respectively.     

Modeling of Hyperconcentrated Sediment-Laden Floods in Lower Yellow River

This paper presents a rapid forecast model for simulating hyperconcentrated sediment-laden floods in the Lower Yellow River. The model is a hybrid of a conventional one-dimensional mathematical model for unsteady sediment-laden flow and an artificial neural networks model for encapsulation of numerical results. The former provides detailed river flood routing information under typical scenarios, whereas the latter extracts modeling outputs from the former and establishes a station-specific model for efficient flood forecasting. Three typical floods that occurred in the Lower Yellow River in 1977, 1982, and 1996 are simulated. Not only the hybrid model predictions are found to be in close agreement with measured data, but also the computational speed is significantly enhanced. It is found that sediment transport is of significance with regard to the flooding behavior of hyperconcentrated flows. Therefore, the model presented herein is of particular use for rivers with high sediment concentration.     

Depth-Averaged Two-Dimensional Numerical Modeling of Unsteady Flow and Nonuniform Sediment Transport in Open Channels

A depth-averaged two-dimensional (2D) numerical model for unsteady flow and nonuniform sediment transport in open channels is established using the finite volume method on a nonstaggered, curvilinear grid. The 2D shallow water equations are solved by the SIMPLE(C) algorithms with the Rhie and Chow’s momentum interpolation technique. The proposed sediment transport model adopts a nonequilibrium approach for nonuniform total-load sediment transport. The bed load and suspended load are calculated separately or jointly according to sediment transport mode. The sediment transport capacity is determined by four formulas which are capable of accounting for the hiding and exposure effects among different size classes. An empirical formula is proposed to consider the effects of the gravity on the sediment transport capacity and the bed-load movement direction in channels with steep slopes. Flow and sediment transport are simulated in a decoupled manner, but the sediment module adopts a coupling procedure for the computations of sediment transport, bed change, and bed material sorting. The model has been tested against several experimental and field cases, showing good agreement between the simulated results and measured data.     

Evaluation of Surface Water Runoff from Fly Ash–Stabilized and Nonstabilized Soil Surfaces

This study evaluated the constituent make up of simulated rainwater runoff from Class C fly ash–stabilized and nonstabilized clay soil using laboratory test pads to assess the potential for impacts to surface water from the use of uncovered fly ash–stabilized soils as potential roadbed material. Recirculated runoff from test pads was sampled and tested during three simulated rainfall events over an 84-day trial period. All samples were analyzed for trace metals. Analytical results from the simulated runoff were screened to identify five indicator parameters in the runoff that were used as the basis for assessing potential environmental effects to surface waters. Runoff water results from fly ash–stabilized test pads for these indicator parameters were compared to water quality benchmarks. Based on the low concentrations measured in runoff relative to applicable criteria, and on the conservative nature of the experimental methods relative to typical field conditions, we concluded that surface runoff from fly ash–stabilized soil would not present significant adverse effects to surface water if used uncovered on low traffic exposed surfaces.     

Case Study: Hydraulic Modeling of Runoff Processes in Ghanaian Inland Valleys

The inland valleys of West Africa are strategic in terms of food security and poverty alleviation, but scientific studies on hydrologic processes happening in these environments have not been well documented. Modeling approaches presented in this paper are an attempt to better comprehend hydraulic phenomena occurring in inland valleys. An inland valley situated in the Northern Region of Ghana is set as the study site. The inland valley comprises well-drained uplands and hydromorphic valley bottoms. There are several earthen dams across the valley bottoms, which are at the same time seasonal wetlands cultivated to rice during the rainy season. A finite volume model for the shallow water equations is developed to numerically simulate surface runoff flows in the valley bottoms during flood events. Innovation is necessitated to handle a series of different hydraulic phenomena. Flux-splitting and data reconstruction techniques are used to achieve stable computation in the complex topography of the valley bottoms. Standard problems of oblique hydraulic jump and dam break flows are used to test the accuracy of the numerical model. The Manning’s roughness coefficient is determined from calibration in another Ghanaian watershed located in the Eastern Region. Using actually observed time series data of rainfall intensity, surface flows during the rainfall events are simulated in the computational domain representing the valley bottoms of the study area. Observed data of water levels in the dams are compared to predictions, and discrepancies between them are examined from the hydrological point of view. In the case of a hypothetical flood event, cascading collapses of the dams and flooding of cultivated fields are reproduced.     

Engineered Levee Breaches for Flood Mitigation

Inland flooding from excessive precipitation and surface runoff remains the cause of extensive damage, loss of property, and human suffering worldwide. In this project, an engineered levee breach was considered as a measure to mitigate floods, and design attributes of the optimal levee breach were examined. The engineered levee breach creates depression waves that destructively interfere with the flood wave to crop peak flood stages. Furthermore, the engineered levee breach permits control over flood stage at a point along a channel where flood stage reduction is most beneficial in terms of hazard mitigation. The roles of flood plain storage, breach size, flood discharge, flood duration, and breach timing in the optimal design of engineered levee breaches were examined using a shallow-water model, and a strategy to design levee breaches is obtained by scaling the model results. It is shown that substantial flood stage reduction can be achieved with an engineered breach. This approach to flood fighting could prove useful in the mitigation of relatively short, extreme floods.     

New Approach to Evaluate Pollutant Removal by Storm-Water Treatment Devices

A methodology was developed to monitor and evaluate the removal of solids and associated constituents by a nutrient separating baffle box (NSBB) storm-water treatment device treating runoff from a 4.3 ha (10.6 acre) residential watershed discharging into the Indian River Lagoon, Florida. The NSBB was monitored over a 359-day time period using autosamplers to quantify water column removal during runoff events, and by quantifying and analyzing solids that accumulated within the NSBB. Flow composited influent and effluent samples were collected to represent water column performance. Event mean concentration (EMC) reduction was moderate (mean: 17%) and variable (range: ?39 to 68%) for suspended solids, and negative for nitrogen, phosphorus, fecal coliforms chromium, and copper. The mass of solids that accumulated in bottom chambers and in a strainer screen was quantified and analyzed for nitrogen, phosphorus, heavy metals, and polycyclic aromatic hydrocarbons. A quantitative evaluative framework was devised to estimate the total pollutant mass removal by NSBB, which consisted of the summation of the separately calculated mass removals for water column, bottom chamber material, and strainer screen material. The water column accounted for only 4% of total solids that accumulated in the NSBB, which was equally divided between bottom chamber and strainer screen. Removal of nitrogen, phosphorus, and metals could be accounted for only by considering mass accumulations. Results suggest that overall assessment of pollutant removal by NSBB must be cognizant of the materials not captured by typical autosamplers: larger size sediment particles, large floating and suspended matter, and the pollutants associated with these materials. Using water column EMCs as the sole measure of performance significantly underestimated loading reduction of storm-water constituents by the NSBB. The monitoring and evaluative methodology applied to the NSBB may be applicable to load reduction evaluations for other storm-water treatment devices with a similar function.     

Method for Estimating Sediment Transport in Ice-Covered Channels

A method is proposed for estimating rates of sediment transport in ice-covered alluvial channels. The method extends existing, open-water procedures for estimating rates of sediment transport to conditions of ice-covered flow. A key aspect of the method is the assessment of flow resistance attributable to bed-surface drag. That assessment is used to estimate rates of bed load and suspended load, and thereby total bed-sediment transport rate. Estimation of ice-covered suspended load additionally entails an approximation whereby open-water suspended load is scaled in proportion to the ratio of a reference sediment concentration for ice-covered flow relative to that for open-water flow. The reference concentration is calculated in terms of bed-load rate and shear velocity attributed to bed-surface drag. Flume data are used to develop the method and tentatively verify it. Field verification of the method presently is hampered by the absence of field data on bed sediment transport in ice-covered channels.     

3D Numerical Modeling of Flow and Sediment Transport in Laboratory Channel Bends

The development of a fully three-dimensional finite volume morphodynamic model, for simulating fluid and sediment transport in curved open channels with rigid walls, is described. For flow field simulation, the Reynolds-averaged Navier–Stokes equations are solved numerically, without reliance on the assumption of hydrostatic pressure distribution, in a curvilinear nonorthogonal coordinate system. Turbulence closure is provided by either a low-Reynolds number k?ω turbulence model or the standard k?ε turbulence model, both of which apply a Boussinesq eddy viscosity. The sediment concentration distribution is obtained using the convection-diffusion equation and the sediment continuity equation is applied to calculate channel bed evolution, based on consideration of both bed load and suspended sediment load. The governing equations are solved in a collocated grid system. Experimental data obtained from a laboratory study of flow in an S-shaped channel are utilized to check the accuracy of the model’s hydrodynamic computations. Also, data from a different laboratory study, of equilibrium bed morphology associated with flow through 90° and 135° channel bends, are used to validate the model’s simulated bed evolution. The numerically-modeled fluid and sediment transportation show generally good agreement with the measured data. The calculated results with both turbulence models show that the low-Reynolds k?ω model better predicts flow and sediment transport through channel bends than the standard k?ε model.