HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

FORESTS AND OTHER FACTORS ASSOCIATED WITH STREAMFLOWS IN EAST TEXAS1

Effects of proportion of watersheds in forest and watershed physiographic factors on mean annual streamflow (1965-76), median flow, and 12 flood flow characteristics were regionally analyzed for 19 unregulated streams in East Texas. Annual streamflow increased with decreasing proportion of forest area. Differences in annual streamflow between full forest cover and bare watersheds could be as much as 200 mm. Other things being equal, the minimum watershed area required to generate 0.142 cm (5 cfs), a criterion used by the U.S. Corps of Engineering in regulating dredge and fill activity for water pollution abatement in East Texas streams, is 70 km² (27 mi²). Of the 31 physio-climatic parameters analyzed, watershed area, percent forest area, shape index, spring precipitation, and annual temperature were the most significant in affecting streamflow characteristics in East Texas. Using 2-3 of these five variables, all of the 14 streamflow characteristics can be estimated with accuracy ranging from acceptable to excellent levels.     

HYDROLOGIC IMPACTS OF A LARGE-SCALE MOUNTAIN PINE BEETLE (DENDROCTONUS PONDEROSAE HOPKINS) EPIDEMIC1

The Jack Creek watershed, a 133 km² (51.5 mi²) drainage in southwestern Montana, was impacted by a mountain pine beetle (Dendroctonus ponderosae Hopkins) epidemic in 1975–1977 which killed an estimated 35 percent of its total timber. Analyses of USGS streamflow data for four years prior to and five years after mortality suggest a 15 percent post-epidemic increase in annual water yield, a two-to three-week advance in the annual hydrograph, a 10 percent increase in low flows and little increase of peak runoff.     

MODELING DEVELOPED COASTAL WATERSHEDS WITH THE AGRICULTURAL NON-POINT SOURCE MODEL1

ABSTRACT: Many coastal states are facing increasing urban growth along their coast lines. The growth has caused urban non-point source nitrogen runoff to be a major contributor to coastal and estuarine enrichment. Water resource managers are responsible for evaluating the impacts from point and non-point sources in developed watersheds and developing strategies to manage future growth. Non-point source models provide an effective approach to these management challenges. The Agricultural Non-Point Source Model (AGNPS) permits the incorporation of important spatial information (soils, landuse, topography, hydrology) in simulating surface hydrology and nitrogen non-point source runoff. The AGNPS model was adapted for developed coastal watersheds by deriving urban coefficients that reflect urban landuse classes and the amount of impervious surface area. Popperdam Creek watershed was used for model parameter development and model calibration. Four additional watersheds were simulated to validate the model. The model predictions of the peak flow and total nitrogen concentrations were close to the field measurements for the five sub-basins simulated. Measured peak flow varied by 30 fold among the sub-basins. The average simulated peak flow was within 14 percent of the average measured peak flow. Measured total nitrogen loads varied over an order of magnitude among the sub-basins yet error between the measured and simulated loads for a given sub-basin averaged 5 percent. The AGNPS model provided better estimates of nitrogen loads than widely used regression methods. The spatial distribution of important watershed characteristics influenced the impacts of urban landuse and projecting future residential expansion on runoff, sediment and nitrogen yields. The AGNPS model provides a useful tool to incorporate these characteristics, evaluate their importance, and evaluate fieldscale to watershed-scale urban impacts.     

HYDROLOGIC CALIBRATION OF THE SWAT MODEL IN A WATERSHED CONTAINING FRAGIPAN SOILS1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) model, designed for use on rural ungaged basins and incorporating a GRASS GIS interface, was used to model the hydrologic response of the Ariel Creek watershed of northeastern Pennsylvania. Model evaluation of daily flow prior to calibration revealed a deviation of runoff volumes (D_v) of 68.3 percent and a Nash-Sutcliffe coefficient of-0.03. Model performance was affected by unusually large observed snowmelt events and the inability of the model to accurately simulate baseflow, which was influenced by the presence of fragipans. Seventy-five percent of the soils in the watershed contain fragipans. Model calibration yielded a D_v of 39.9 percent and a Nash-Sutcliffe coefficient of 0.04, when compared on a daily basis. Monthly comparisons yielded a Nash-Sutcliffe coefficient of 0.14. Snowmelt events in the springs of 1993 and 1994, which were unusually severe, were not adequately simulated. Neglecting these severe events, which produced the largest and third largest measured flows for the period of record, a D_v of 4.1 percent and Nash-Sutcliffe coefficient of 0.20 were calculated on a daily comparison, while on a monthly basis the Nash-Sutciffe coefficient was 0.55. These results suggest that the SWAT model is better suited to longer period simulations of hydrologic yields. Baseflow volumes were accurately simulated after calibration (D_v= -0.2 percent). Refinements made to the algorithms controlling subsurface hydrology and snowmelt, to better represent the presence of fragipans and snowmelt events, would likely improve model performance.     

EVALUATING HSPF IN AN ARID,URBANIZED WATERSHED1

ABSTRACT: The Hydrologic Simulation Program‐FORTRAN (HSPF) is a powerful time variable hydrologic model that has rarely been applied in arid environments. Here, the performance of HSPF in southern California was assessed, testing its ability to predict annual volume, daily average flow, and hourly flow. The model was parameterized with eight land use categories and physical watershed characteristics. It was calibrated using rainfall and measured flow over a five‐year period in a predominantly undeveloped watershed and it was validated using a subsequent 4‐year period. The process was repeated in a separate, predominantly urbanized watershed over the same time span. Annual volume predictions correlated well with measured flow in both the undeveloped and developed watersheds. Daily flow predictions correlated well with measured flow following rain events, but predictions were poor during extended dry weather periods in the developed watershed. This modeling difficulty during dry‐weather periods reflects the large influence of, and the poor accounting in the model for, artificially introduced water from human activities, such as landscape overwatering, that can be important sources of water in urbanized arid environments. Hourly flow predictions mistimed peak flows, reflecting spatial and temporal heterogeneity of rainfall within the watershed. Model correlation increased considerably when predictions were averaged over longer time periods, reaching an asymptote after an 11‐hour averaging window.     

Improvement in HSPF's Low‐Flow Predictions by Implementation of a Power Law Groundwater Storage‐Discharge Relationship

We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km² drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.     

Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain1

Epps, Thomas H., Daniel R. Hitchcock, Anand D. Jayakaran, Drake R. Loflin, Thomas M. Williams, and Devendra M. Amatya, 2012. Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12000 Abstract: Hydrologic monitoring was conducted in two first‐order lower coastal plain watersheds in South Carolina, United States, a region with increasing growth and land use change. Storm events over a three‐year period were analyzed for direct runoff coefficients (ROC) and the total storm response (TSR) as percent rainfall. ROC calculations utilized an empirical hydrograph separation method that partitioned total streamflow into sustained base flow and direct runoff components. ROC ratios ranged from 0 to 0.32 on the Upper Debidue Creek (UDC) watershed and 0 to 0.57 on Watershed 80 (WS80); TSR results ranged from 0 to 0.93 at UDC and 0.01 to 0.74 at WS80. Variability in event runoff generation was attributed to seasonal trends in water table elevation fluctuation as regulated by evapotranspiration. Groundwater elevation breakpoints for each watershed were identified based on antecedent water table elevation, streamflow, ROCs, and TSRs. These thresholds represent the groundwater elevation above which event runoff generation increased sharply in response to rainfall. For effective coastal land use decision making, baseline watershed hydrology must be understood to serve as a benchmark for management goals, based on both seasonal and event‐based surface and groundwater interactions.     

FLOW ANALYSIS OF LANDSLIDE DAMMED LAKE WATERSHEDS: A CASE STUDY1

Abstract: The Chi-Chi earthquake, which occurred on September 21, 1999, and had a magnitude of 7.3 on the Richter scale, resulted in an extensive landslide that blocked the Ching-Shui Creek in Taiwan, forming a large lake with a storage volume of 40 million m³. This paper describes an analytical procedure used to perform flow analysis of the Tsao-Ling watershed, which includes the new landslide dammed lake. In this study, a digital elevation model was applied to obtain the watershed geomorphic factors and stage-area storage function of the landslide dammed lake. Satellite images were used to identify the landslide area and the land cover change that occurred as a result of the earthquake. Two topography-based runoff models were applied for long term and short term streamflow analyses of the watershed because the watershed upstream of the landslide dam was ungauged. The simulated daily flow and storm runoff were verified using limited available measured data in the watershed, and good agreement was obtained. The proposed analytical procedure for flow analysis is considered promising for application to other landslide dammed lake watersheds.     

Model Based Nitrate TMDLs for Two Agricultural Watersheds of Southeastern Minnesota1

Abstract: In this study, a set of nitrogen reduction strategies were modeled to evaluate the feasibility of improving water quality to meet total maximum daily loads (TMDLs) in two agricultural watersheds. For this purpose, a spatial‐process model was calibrated and used to predict monthly nitrate losses (1994‐96) from Sand and Bevens Creek watersheds located in south‐central Minnesota. Statistical comparison of predicted and observed flow and nitrate losses gave r² coefficients of 0.75 and 0.70 for Sand Creek watershed and 0.72 and 0.67 for Bevens Creek watershed, respectively. Modeled alternative agricultural management scenarios included: six different N application rates over three application timings and three different percentages of crop land with subsurface drainage. Predicted annual nitrate losses were then compared with nitrate TMDLs assuming a 30% reduction in observed nitrate losses is required. Reductions of about 33 (8.6 to 5.8 kg/ha) and 35% (23 to 15 kg/ha) in existing annual nitrate losses are possible for Sand and Bevens Creek watersheds, respectively, by switching the timing of fertilizer application from fall to spring. Trends towards increases in tile‐drained crop land imply that attaining nitrate TMDLs in future may require other alternative management practices in addition to fertilizer management such as partial conversion of crop land to pasture.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

Urban Streamflow Response to Imported Water and Water Conservation Policies in Los Angeles,California

Los Angeles has a long history of importing water; however, drought, climate change, and environmental mitigation have forced the City to focus on developing more local water sources (target of 50% local supply by 2035). This study aims to improve understanding of water cycling in Los Angeles, including the impacts of imported water and water conservation policies. We evaluate the influence of local water restrictions on discharge records for 12 years in the Ballona Creek (urban) and Topanga Creek (natural) watersheds. Results show imported water has significantly altered the timing and volume of streamflow in the urban Ballona watershed, resulting in runoff ratios above one (more streamflow than precipitation). Further analysis comparing pre‐ vs. during‐mandatory water conservation periods shows there is a significant decrease in dry season streamflow during‐conservation in Ballona, indicating that prior to conservation efforts, heavy irrigation and other outdoor water use practices were contributing to streamflow. The difference between summer streamflow pre‐ vs. during‐conservation is enough to serve 160,000 customers in Los Angeles. If Los Angeles returns to more watering days, educating the public on proper irrigation rates is critical for ensuring efficient irrigation and conserving water; however, if water restrictions remain in place, the City must take the new flow volumes into account for complying with water quality standards in the region.     

SIMULATING HYDROLOGIC AND WATER QUALITY IMPACTS IN AN URBANIZING WATERSHED1

ABSTRACT: The Hydrologic Simulation Program‐Fortran (HSPF) was calibrated and used to assess the future effects of various land development scenarios on water quality in the Polecat Creek watershed in Caroline County, Virginia. Model parameters related to hydrology and water quality were calibrated and validated using observed stream flow and water quality data collected at the watershed outlet and the outlets of two sub water sheds. Using the county's Comprehensive Plan, land use scenarios were developed by taking into account the trends and spatial distributions of future development. The simulation results for the various land use scenarios indicate that runoff volume and peak rate increased as urban areas increased. Urbanization also increased sediment loads mainly due to increases in channel erosion. Constituent loads of total Kjeldal nitrogen, orthophosphorus, and total phosphorous for Polecat Creek watershed slightly decreased under future development scenarios. These reductions are due to increases in urban areas that typically contribute smaller quantities of nitrogen and phosphorous, as compared to agricultural areas. However, nitrate loads increased for the future land use scenarios, as compared to the existing land use. The increases in nitrate loads may result from increases in residential land and associated fertilizer use and concurrent decreases in forested land. The procedures used in this paper could assist local, state, and regional policy makers in developing land management strategies that minimize environmental impacts while allowing for future development.     

CALIBRATION OF STORM LOADS IN THE SOUTH PRONG WATERSHED,FLORIDA, USING BASINS/HSPF1

ABSTRACT: The South Prong watershed is a major tributary system of the Sebastian River and adjacent Indian River Lagoon. Continued urbanization of the Sebastian River drainage basin and other watersheds of the Indian River Lagoon is expected to increase runoff and nonpoint source pollutant loads. The St. Johns River Water Management District developed watershed simulation models to estimate potential impacts on the ecological systems of receiving waters and to assist planners in devising strategies to prevent further degradation of water resources. In the South Prong system, a storm water sampling program was carried out to calibrate the water quality components of the watershed model for total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN). During the period of May to November 1999, water quality and flow data were collected at three locations within the watershed. Two of the sampling stations were located at the downstream end of major watercourses. The third station was located at the watershed outlet. Five storm events were sampled and measured at each station. Sampling was conducted at appropriate intervals to represent the rising limb, peak, and recession limb of each storm event. The simulations were handled by HSPF (Hydrologic Simulation Program‐Fortran). Results include calibration of the hydrology and calibration of the individual storm loads. The hydrologic calibration was continuous over the period 1994 through 1999. Simulated storm runoff, storm loads, and event mean concentrations were compared with their corresponding observed values. The hydrologic calibration showed good results. The outcome of the individual storm calibrations was mixed. Overall, however, the simulated storm loads agreed reasonably well with measured loads for a majority of the storms.     

Chemical transport from paired agricultural and restored prairie watersheds

A five-year record of streamflow and chemical sampling data was evaluated to assess the effects of large-scale prairie restoration on transport of NO3-N, Cl, and SO4 loads from paired 5,000-ha watersheds located in Jasper County, Iowa. Water quality conditions monitored during land use conversion from row crop agriculture to native prairie in the Walnut Creek watershed were compared with a highly agricultural control watershed (Squaw Creek). Combining hydrograph separation with a load estimation program, baseflow and stormflow loads of NO3-N, Cl, and SO4 were estimated at upstream and downstream sites on Walnut Creek and a downstream site on Squaw Creek. Chemical export in both watersheds was found to occur primarily with baseflow, with baseflow transport greatest during the late summer and fall. Lower Walnut Creek watershed, which contained the restored prairie areas, exported less NO3-N and Cl compared with upper Walnut Creek and Squaw Creek watersheds. Average flow-weighted concentrations of NO3-N exceeded 10 mg/L in upper Walnut Creek and Squaw Creek, but were estimated to be 6.6 mg/L in lower Walnut Creek. Study results demonstrate the utility of partitioning loads into baseflow and stormflow components to identify sources of pollutant loading to streams.     

MODELING SEDIMENT AND PHOSPHORUS LOSSES IN AN AGRICULTURAL WATERSHED TO MEET TMDLs1

ABSTRACT: This paper studies the effectiveness of alternative farm management strategies at improving water quality to meet Total Maximum Daily Loads (TMDLs) in agricultural watersheds. A spatial process model was calibrated using monthly flow, sediment, and phosphorus (P) losses (1994 to 1996) from Sand Creek watershed in south‐central Minnesota. Statistical evaluation of predicted and observed data gave r² coefficients of 0.75, 0.69, and 0.49 for flow (average 4.1 m³/s), sediment load (average 0.44 ton/ha), and phosphorus load (average 0.97 kg/ha), respectively. The calibrated model was used to evaluate the effects of conservation tillage, conversion of crop land to pasture, and changes in phosphorus fertilizer application rate on pollutant loads. TMDLs were developed for sediment and P losses based on existing water quality standards and guidelines. Observed annual sediment and P losses exceeded these TMDLs by 59 percent and 83 percent, respectively. A combination of increased conservation tillage, reduced application rates of phosphorus fertilizer, and conversion of crop land to pasture could reduce sediment and phosphorus loads by 23 percent and 20 percent of existing loads, respectively. These reductions are much less than needed to meet TMDLs, suggesting that control of sediment using buffer strips and control of point sources of phosphorus are needed for the remaining reductions.     

Calibration and Verification of SWMM for Low Impact Development

The Storm Water Management Model was used to simulate runoff and nutrient export from a low impact development (LID) watershed and a watershed using traditional runoff controls. Predictions were compared to observed values. Uncalibrated simulations underpredicted weekly runoff volume and average peak flow rates from the multiple subcatchment LID watershed by over 80%; the single subcatchment traditional watershed had better predictions. Saturated hydraulic conductivity, Manning's n for swales, and initial soil moisture deficit were sensitive parameters. After calibration, prediction of total weekly runoff volume for the LID and traditional watersheds improved to within 12 and 5% of observed values, respectively. For the validation period, predicted total weekly runoff volumes for the LID and traditional watersheds were within 6 and 2% of observed values, respectively. Water quality simulation was less successful, Nash–Sutcliffe coefficients >0.5 for both calibration and validation periods were only achieved for prediction of total nitrogen export from the LID watershed. Simulation of a 100‐year, 24‐h storm resulted in a runoff coefficient of 0.46 for the LID watershed and 0.59 for the traditional watershed. Results suggest either calibration is needed to improve predictions for LID watersheds or expanded look‐up tables for Green–Ampt infiltration parameter values that account for compaction of urban soil and antecedent conditions are needed.     

SIMULATING THREE-DIMENSIONAL GROUND WATER RESPONSE IN A SMALL MOUNTAINOUS WATERSHED1

ABSTRACT: Snowmelt from deep mountainous snowpacks is seldom rapid enough to exceed infiltration rates; thus, the source of streamflow in many mountainous watersheds is snowmelt recharge through shallow ground water systems. The hydrologic response and interaction between surface and sub-surface flow processes in these watersheds, which is controlled by basin structure, the spatial distribution of snowmelt, and the hydrogeology of the subsurface, are not well understood. The purpose of this study was to test a three-dimensional ground water model using simulated snowmelt input to simulate ground water response to spatially distributed snowmelt on the Upper Sheep Creek Watershed located within the Reynolds Creek Experimental Watershed in Southwestern Idaho. The model was used to characterize the mountainous aquifer and to delineate the subsurface flow mechanisms. Difficulty in finding a reasonable combination of grid spacing and time stepping within the model was encountered due to convergence problems with the Picard solution to the non-linear variably saturated ground water flow equations. Simulation results indicated that flow may be either unconfined or confined depending on inflow rate and hydrogeologic conditions in the watershed. The flow mechanism had a much faster response time when confined flow occurred. Response to snowmelt from a snow drift approximately 90 m away took only a few hours when flow was confined. Simulated results showed good agreement with piezometer measurements both in magnitude and timing; however, convergence problems with the Picard solution limited applicability of the model.     

CALIBRATION OF A WATER-BALANCE MODEL FOR SMALL WATERSHEDS IN EASTERN OREGON1

ABSTRACT: The BURP water-balance model was calibrated for 13 small (0.46 to 7.00 mi²), forested watersheds in the Blue Mountains of eastern Oregon where snowmelt is the dominant source of runoff. BURP is the model name and is not an acronym. Six of the 16 parameters in BURP were calibrated. The subsurface recession coefficient and three subsurface water-storage parameters were most sensitive for simulating monthly flow. Calibrated subsurface recession coefficients ranged from 0.988 to 0.998. The subsurface-water storage parameters were calibrated at between 20 to 120 percent of their initial values obtained from a category III soil survey. That reconnaissance-level survey was apparently too broad to accurately reflect subsurface-water storage in small watersheds. Tests of model performance showed BURP is capable of producing accurate simulations of monthly flow for mountainous, snow-dominated watersheds with shallow (< 4 ft) soils when calibrated with 2 to 4 years of streamflow data. A regression of observed versus simulated monthly flows with data from all watersheds combined showed that BURP accounted for 85 percent of the variability in observed flows, which ranged from 0.01 to 20.8 inches, but underpredicted high flow months, with a slope of 1.15 that is significantly different from 1.0 (p = 0.05). Without prior calibration, subsurface-water storage parameters appeared to be the greatest source of potential error.