Impact of Watershed Subdivision and Soil Data Resolution on SWAT Model Calibration and Parameter Uncertainty1

Abstract: Impact of watershed subdivision and soil data resolution on Soil Water Assessment Tool (SWAT) model calibration and parameter uncertainty is investigated by creating 24 different watershed model configurations for two study areas in northern Indiana. SWAT autocalibration tool is used to calibrate 14 parameters for simulating seven years of daily streamflow records. Calibrated parameter sets are found to be different for all 24 watershed configurations, however in terms of calibrated model output, their effect is minimal. In some cases, autocalibration is followed by manual calibration to correct for low flows, which were underestimated during autocalibration. In addition to normal validation using four years of streamflow data for each calibrated parameter set, cross‐validation (using a calibrated parameter set from one model configuration to validate observations using another configuration) is performed to investigate the effect of different model configurations on streamflow prediction. Results show that streamflow output during cross‐validation is not affected, thus highlighting the non‐unique nature of calibrated parameters in hydrologic modeling. Finally, parameter uncertainty is investigated by extracting good parameter sets during the autocalibration process. Parameter uncertainty analysis suggests that significant parameters show very narrow range of uncertainty across different watershed configurations compared with nonsignificant parameters. Results from recalibration of some configurations using only six significant parameters were comparable to that from calibration using 14 parameters, suggesting that including fewer significant parameters could reduce the uncertainty arising from model parameters, and also expedite the calibration process.     

SWATmodel: A Multi‐Operating System,Multi‐Platform SWAT Model Package in R

The Soil and Water Assessment Tool (SWAT) model (Arnold et al., 1998) is a popular watershed management tool. Currently, the SWAT model, actively supported by the U.S. Department of Agriculture and Texas A&M, operates only on Microsoft^® Windows, which hinders modelers that use other operating systems (OS). This technical note introduces the Comprehensive R Archive Network (CRAN) distributed “SWATmodel” package which allows SWAT 2005 and 2012 to be widely distributed and run as a linear model‐like function on multiple OS and processor platforms. This allows researchers anywhere in the world using virtually any OS to run SWAT. In addition to simplifying the use of SWAT across computational platforms, the SWATmodel package allows SWAT modelers to utilize the analytical capabilities, statistical libraries, modeling tools, and programming flexibility inherent to R. The software allows watershed modelers to develop a simple hydrological watershed model conceptualization of the SWAT model and to obtain a first approximation of the minimum expected results a more complicated model should deliver. As a proof of concept, we test the SWAT model by initializing and calibrating 314 U.S. Geological Survey stream gages in the Chesapeake Bay watershed and present the results.     

SENSITIVITY ANALYSIS,CALIBRATION, AND VALIDATIONS FOR A MULTISITE AND MULTIVARIABLE SWAT MODEL1

The ability of a watershed model to mimic specified watershed processes is assessed through the calibration and validation process. The Soil and Water Assessment Tool (SWAT) watershed model was implemented in the Beaver Reservoir Watershed of Northwest Arkansas. The objectives were to: (1) provide detailed information on calibrating and applying a multisite and multivariable SWAT model; (2) conduct sensitivity analysis; and (3) perform calibration and validation at three different sites for flow, sediment, total phosphorus (TP), and nitrate‐nitrogen (NO₃‐N) plus nitrite‐nitrogen (NO₂‐N). Relative sensitivity analysis was conducted to identify parameters that most influenced predicted flow, sediment, and nutrient model outputs. A multi objective function was defined that consisted of optimizing three statistics: percent relative error (RE), Nash‐Sutcliffe Coefficient (R_NS²), and coefficient of determination (R²). This function was used to successfully calibrate and validate a SWAT model of Beaver Reservoir Watershed at multi‐sites while considering multivariables. Calibration and validation of the model is a key factor in reducing uncertainty and increasing user confidence in its predictive abilities, which makes the application of the model effective. Information on calibration and validation of multisite, multivariable SWAT models has been provided to assist watershed modelers in developing their models to achieve watershed management goals.     

Assessing Parameter Uncertainty of a Semi‐Distributed Hydrology Model for a Shallow Aquifer Dominated Environmental System

This paper examines the performance of a semi‐distributed hydrology model (i.e., Soil and Water Assessment Tool [SWAT]) using Sequential Uncertainty FItting (SUFI‐2), generalized likelihood uncertainty estimation (GLUE), parameter solution (ParaSol), and particle swarm optimization (PSO). We applied SWAT to the Waccamaw watershed, a shallow aquifer dominated Coastal Plain watershed in the Southeastern United States (U.S.). The model was calibrated (2003‐2005) and validated (2006‐2007) at two U.S. Geological Survey gaging stations, using significant parameters related to surface hydrology, hydrogeology, hydraulics, and physical properties. SWAT performed best during intervals with wet and normal antecedent conditions with varying sensitivity to effluent channel shape and characteristics. In addition, the calibration of all algorithms depended mostly on Manning's n‐value for the tributary channels as the surface friction resistance factor to generate runoff. SUFI‐2 and PSO simulated the same relative probability distribution tails to those observed at an upstream outlet, while all methods (except ParaSol) exhibited longer tails at a downstream outlet. The ParaSol model exhibited large skewness suggesting a global search algorithm was less capable of characterizing parameter uncertainty. Our findings provide insights regarding parameter sensitivity and uncertainty as well as modeling diagnostic analysis that can improve hydrologic theory and prediction in complex watersheds. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Implications of Conceptual Channel Representation on SWAT Streamflow and Sediment Modeling

Hydrologic modeling outputs are influenced by how a watershed system is represented. Channel routing is a typical example of the mathematical conceptualization of watershed landscape and processes in hydrologic modeling. We investigated the sensitivity of accuracy, equifinality, and uncertainty of Soil and Water Assessment Tool (SWAT) modeling to channel dimensions to demonstrate how a conceptual representation of a watershed system affects streamflow and sediment modeling. Results showed the amount of uncertainty and equifinality strongly responded to channel dimensions. On the other hand, the model performance did not significantly vary with the changes in the channel representation due to the degree of freedom allowed by the conceptual nature of hydrologic modeling in the parameter calibration. Such findings demonstrated good modeling performance statistics do not necessarily mean small output uncertainty, and partial improvements in the watershed representation may neither increase modeling accuracy nor reduce uncertainty. We also showed the equifinality and uncertainty of hydrologic modeling are case‐dependent rather than specific to models or regions, suggesting great caution should be used when attempting to transfer uncertainty analysis results to other modeling studies, especially for ungauged watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Large‐Scale Fine‐Resolution Hydrological Modeling Using Parameter Regionalization in the Missouri River Basin

This study simulated crop and water yields in the Missouri River Basin (MRB; 1,371,000 km²), one of the largest river basins in the United States, using the Soil and Water Assessment Tool (SWAT) at a fine resolution of 12‐digit Hydrological Unit Codes (HUCs) using the regionalization calibration approach. Very few studies have simulated the entire MRB, and those that have developed were at a coarser resolution of 8‐digit HUCs and were minimally calibrated. The MRB was first divided into three subbasins and was further divided into eleven regions. A “head watershed” was selected in each region and was calibrated for crop and water yields. The parameters from the calibrated head watershed were extrapolated to other subwatersheds in the region to complete comprehensive spatial calibration. The simulated crop yields at the head watersheds were in close agreement with observed crop yields. Spatial validation of the aggregated crop yields resulted in reasonable predictions for all crops except dryland corn in a few regions. Simulated and observed water yields in head watersheds and also in the validation locations were in close agreement in naturalized streams and poor agreement in streams with high groundwater‐surface water interactions and/or reservoirs found upstream of the gauges. Overall, the SWAT model was able to reasonably capture the hydrological and crop growth dynamics occurring in the basin despite some limitations.     

Testing of the Modified Streambank Erosion and Instream Phosphorus Routines for the SWAT Model

In some watersheds, streambanks are a source of two major pollutants, phosphorus (P) and sediment. P originating from both uplands and streambanks can be transported and stored indefinitely on floodplains, streambanks, and in closed depressions near the stream. The objectives of this study were to (1) test the modified streambank erosion and instream P routines for the Soil and Water Assessment Tool (SWAT) model in the Barren Fork Creek watershed in northeast Oklahoma, (2) predict P in the watershed with and without streambank‐derived P, and (3) determine the significance of streambank erosion P relative to overland P sources. Measured streambank and channel parameters were incorporated into a flow‐calibrated SWAT model and used to estimate streambank erosion and P for the Barren Fork Creek using modified streambank erosion and instream P routines. The predicted reach‐weighted streambank erosion was 40 kg/m vs. the measured 42 kg/m. Streambank erosion contributed 47% of the total P to the Barren Fork Creek and improved P predictions compared to observed data, especially during the high‐flow events. Of the total P entering the stream system, approximately 65% was removed via the watershed outlet and 35% was stored in the floodplain and stream system. This study successfully applied the SWAT model's modified streambank erosion and instream P routines and demonstrated that streambank‐derived P can improve P modeling at the watershed scale. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Evaluation of Existing and Modified Wetland Equations in the SWAT Model

The ability to accurately simulate flow and nutrient removal in treatment wetlands within an agricultural, watershed‐scale model is needed to develop effective plans for meeting nutrient reduction goals associated with protection of drinking water supplies and reduction of the Gulf of Mexico hypoxic zone. The objectives of this study were to incorporate new equations for wetland hydrology and nutrient removal in Soil and Water Assessment Tool (SWAT), compare model performance using original and improved equations, and evaluate the ramifications of errors in watershed and tile drain simulation on prediction of NO₃‐N dynamics in wetlands. The modified equations produced Nash‐Sutcliffe Efficiency values of 0.88 to 0.99 for daily NO₃‐N load predictions, and percent bias values generally less than 6%. However, statistical improvement over the original equations was marginal and both old and new equations provided accurate simulations. The new equations reduce the model's dependence on detailed monitoring data and hydrologic calibration. Additionally, the modified equations increase SWAT's versatility by incorporating a weir equation and an irreducible nutrient concentration and temperature coefficient. Model improvements enhance the utility of SWAT for simulating flow and nutrients in wetlands and other impoundments, although performance is limited by the accuracy of inflow and NO₃‐N predictions from the contributing watershed. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

A Simple Process‐Based Snowmelt Routine to Model Spatially Distributed Snow Depth and Snowmelt in the SWAT Model1

Abstract: We present a method to integrate a process‐based (PB) snowmelt model that requires only daily temperature and elevation information into the Soil and Water Assessment Tool (SWAT) model. The model predicts the spatiotemporal snowpack distribution without adding additional complexity, and in fact reduces the number of calibrated parameters. To demonstrate the utility of the PB model, we calibrate the PB and temperature‐index (TI) SWAT models to optimize agreement with stream discharge on a 46‐km² watershed in northwestern Idaho, United States, for 10 individual years and use the calibrated parameters for the year with the best agreement to run the model for 15 remaining years. Stream discharge predictions by the PB and TI model were similar, although the PB model simulated snowmelt more accurately than the TI model for the remaining 15‐year period. Spatial snow distributions predicted by the PB model better matched observations from LandSat imagery and a SNOTEL station. Results for this watershed show that including PB snowmelt in watershed models is feasible, and calibration of TI‐based watershed models against discharge can incorrectly predict snow cover.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

Addressing data challenges in riverine nutrient load modeling of an intensively managed agro-industrial watershed

Data limitations often challenge the reliability of water quality models, especially in intensively managed watersheds. While numerous studies report successful hydrological model setup and calibration, few have addressed in detail the data challenges for multisite and multivariable model calibration to an intensively managed watershed. In this study, we address some of these challenges based on our reflective experience calibrating the Soil and Water Assessment Tool (SWAT) to the Upper Sangamon River Watershed in central Illinois based on daily flow, annual crop yield, and monthly sediment, nitrate, and total phosphorus loads. We highlight some challenges in SWAT calibration processes due to data errors and inconsistencies, and insufficient precipitation and water quality observations. Following, we demonstrate the merits of additional weather and water quality observations that could help reduce input uncertainties, and we provide suggestions for selecting appropriate observations for the model calibration. After dealing with the data issues, we show that the SWAT model could be calibrated with acceptable results for the case study watershed.     

Spatial and Temporal Evaluation of Hydrological Response to Climate and Land Use Change in Three South Dakota Watersheds

This study analyzed changes in hydrology between two recent decades (1980s and 2010s) with the Soil and Water Assessment Tool (SWAT) in three representative watersheds in South Dakota: Bad River, Skunk Creek, and Upper Big Sioux River watersheds. Two SWAT models were created over two discrete time periods (1981‐1990 and 2005‐2014) for each watershed. National Land Cover Datasets 1992 and 2011 were, respectively, ingested into 1981‐1990 and 2005‐2014 models, along with corresponding weather data, to enable comparison of annual and seasonal runoff, soil water content, evapotranspiration (ET), water yield, and percolation between these two decades. Simulation results based on the calibrated models showed that surface runoff, soil water content, water yield, and percolation increased in all three watersheds. Elevated ET was also apparent, except in Skunk Creek watershed. Differences in annual water balance components appeared to follow changes in land use more closely than variation in precipitation amounts, although seasonal variation in precipitation was reflected in seasonal surface runoff. Subbasin‐scale spatial analyses revealed noticeable increases in water balance components mostly in downstream parts of Bad River and Skunk Creek watersheds, and the western part of Upper Big Sioux River watershed. Results presented in this study provide some insight into recent changes in hydrological processes in South Dakota watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Automatic Calibration of Hydrologic Models With Multi‐Objective Evolutionary Algorithm and Pareto Optimization1

Abstract: In optimization problems with at least two conflicting objectives, a set of solutions rather than a unique one exists because of the trade‐offs between these objectives. A Pareto optimal solution set is achieved when a solution cannot be improved upon without degrading at least one of its objective criteria. This study investigated the application of multi‐objective evolutionary algorithm (MOEA) and Pareto ordering optimization in the automatic calibration of the Soil and Water Assessment Tool (SWAT), a process‐based, semi‐distributed, and continuous hydrologic model. The nondominated sorting genetic algorithm II (NSGA‐II), a fast and recent MOEA, and SWAT were called in FORTRAN from a parallel genetic algorithm library (PGAPACK) to determine the Pareto optimal set. A total of 139 parameter values were simultaneously and explicitly optimized in the calibration. The calibrated SWAT model simulated well the daily streamflow of the Calapooia watershed for a 3‐year period. The daily Nash‐Sutcliffe coefficients were 0.86 at calibration and 0.81 at validation. Automatic multi‐objective calibration of a complex watershed model was successfully implemented using Pareto ordering and MOEA. Future studies include simultaneous automatic calibration of water quality and quantity parameters and the application of Pareto optimization in decision and policy‐making problems related to conflicting objectives of economics and environmental quality.     

Application of SWAT with and without Variable Source Area Hydrology to a Large Watershed

In this study, two different versions of the Soil and Water Assessment Tool (SWAT) model were used to simulate the hydrology and biogeochemical response of the Cannonsville Reservoir watershed, in New York. The first version distributes overland flow in ways that are consistent with variable source area (VSA) hydrology driven by saturation excess runoff, whereas the second version is the standard version of SWAT. These two models were each calibrated for streamflow (Flow), particulate phosphorus (PP), total dissolved phosphorus (TDP), and sediment (Sed) against measured data from the 1,200 km² Cannonsville watershed. The standard version of the model yielded an r² between the measured and simulated data of 0.85, 0.73, 0.70, and 0.72 for Flow, Sed, TDP, and PP, respectively. The VSA version yielded an r² of 0.84, 0.69, 0.72, and 0.53 for Flow, Sed, TDP, and PP, respectively. The two models were then used to determine the maximum upper bound on the reduction in phosphorus loading by removing all of the corn in the watershed. The average reductions between the two models were 65 and 37% for PP and TDP, respectively. The VSA version was also used to estimate the effect of moving corn land in the watershed from the wettest, most runoff prone areas to the driest, least runoff prone areas, which cannot be done directly with the standard SWAT model.     

Runoff Simulation of the Headwaters of the Yellow River Using The SWAT Model With Three Snowmelt Algorithms1

Abstract: The Soil and Water Assessment Tool (SWAT) model combined with different snowmelt algorithms was evaluated for runoff simulation of an 114,345 km² mountainous river basin (the headwaters of the Yellow River), where snowmelt is a significant process. The three snowmelt algorithms incorporated into SWAT were as follows: (1) the temperature‐index, (2) the temperature‐index plus elevation band, and (3) the energy budget based SNOW17. The SNOW17 is more complex than the temperature‐based snowmelt algorithms, and requires more detailed meteorological and topographical inputs. In order to apply the SNOW17 in the SWAT framework, SWAT was modified to operate at the pixel scale rather than the normal Hydrologic Response Unit scale. The three snowmelt algorithms were evaluated under two parameter scenarios, the default and the calibrated parameters scenarios. Under the default parameters scenario, the parameter values were determined based on a review of the current literature. The purpose of this type of evaluation was to assess the applicability of SWAT in ungauged basins, where there is little observed data available for calibration. Under the calibrated parameters scenario, the parameters were calibrated using an automatic calibration program, the Shuffled Complex Evolution (SCE‐UA). The purpose of this type of evaluation was to assess the applicability of SWAT in gauged basins. Two time periods (1975‐1985 and 1986‐1990) of monthly runoff data were used in this study to evaluate the performance of SWAT with different snowmelt algorithms. Under the default parameters scenario, the SWAT model with complex energy budget based SNOW17 performed the best for both time periods. Under the calibrated parameters scenario, the parameters were calibrated using monthly runoff from 1975‐1985 and validated using monthly runoff from 1986‐1990. After parameter calibration, the performance of SWAT with the three snowmelt algorithms was improved from the default parameters scenario. Further, the SWAT model with temperature‐index plus elevation band performed as well as the SWAT model with SNOW17. The SWAT model with temperature‐index algorithm performed the poorest for both time periods under both scenarios. Therefore, it is suggested that the SNOW17 model be used for modeling ungauged basins; however, for gauged basins, the SNOW17 and simple temperature‐index plus elevation band models could provide almost equally good runoff simulation results.     

Evaluation of SWAT Impoundment Modeling Methods in Water and Sediment Simulations

Worldwide studies show 80%–90% of all sediments eroded from watersheds is trapped within river networks such as reservoirs, ponds, and wetlands. To represent the impact of impoundments on sediment routing in watershed modeling, Soil and Water Assessment Tool (SWAT) developers recommend to model reservoirs, ponds, and wetlands using impoundment tools (ITs). This study evaluates performance of SWAT ITs in the modeling of a small, agricultural watershed dominated by lakes and wetlands. The study demonstrates how to incorporate impoundments into the SWAT model, and discusses and evaluates involved parameters. The study then recommends an appropriate calibration sequence, i.e., landscape parameters calibration, followed by pond/wetlands calibration, then channel parameter calibrations, and lastly, reservoir parameter calibration. Results of this study demonstrate not following SWAT recommendation regarding modeling water land use as an impoundment depreciates SWAT performance, and may lead to misplaced calibration efforts and model over‐calibration. Further, the chosen method to model impoundments’ outflow significantly impacts sediment loads in the watershed, while streamflow simulation is not very sensitive. This study also allowed calculation of mass accumulation rates in modeled impoundments where the annual mass accumulation rate in wetlands (2.3 T/ha/yr) was 39% higher than mass accumulation rate in reservoirs (1.4 T/ha/yr).     

Evaluation of a Watershed Model for Estimating Daily Flow Using Limited Flow Measurements1

Abstract: The Soil and Water Assessment Tool (SWAT) model was evaluated for estimation of continuous daily flow based on limited flow measurements in the Upper Oyster Creek (UOC) watershed. SWAT was calibrated against limited measured flow data and then validated. The Nash‐Sutcliffe model Efficiency (NSE) and mean relative error values of daily flow estimations were 0.66 and 15% for calibration, and 0.56 and 4% for validation, respectively. Also, further evaluation of the model’s estimation of flow at multiple locations was conducted with parametric paired t‐test and nonparametric sign test at a 95% confidence level. Among the five main stem stations, four stations were statistically shown to have good agreement between predicted and measured flows. SWAT underestimated the flow of the fifth main stem station possibly because of the existence of complex flood control measures near to the station. SWAT estimated the daily flow at one tributary station well, but with relatively large errors for the other two tributaries. The spatial pattern of predicted flows matched the measured ones well. Overall, it was concluded from the graphical comparisons and statistical analyses of the model results that SWAT was capable of reproducing continuous daily flows based on limited flow data as is the case in the UOC watershed.     

ARCGIS‐SWAT: A GEODATA MODEL AND GIS INTERFACE FOR SWAT1

This paper presents ArcGIS‐SWAT, a geodata model and geographic information system (GIS) interface for the Soil and Water Assessment Tool (SWAT). The ArcGIS‐SWAT data model is a system of geodatabases that store SWAT geographic, numeric, and text input data and results in an organized fashion. Thus, it is proposed that a single and comprehensive geodatabase be used as the repository of a SWAT simulation. The ArcGIS‐SWAT interface uses programming objects that conform to the Component Object Model (COM) design standard, which facilitate the use of functionality of other Windows‐based applications within ArcGIS‐SWAT. In particular, the use of MS Excel and MATLAB functionality for data analysis and visualization of results is demonstrated. Likewise, it is proposed to conduct hydrologic model integration through the sharing of information with a not‐model‐specific hub data model where information common to different models can be stored and from which it can be retrieved. As an example, it is demonstrated how the Hydrologic Modeling System (HMS) ‐ a computer application for flood analysis ‐ can use information originally developed by ArcGIS‐SWAT for SWAT. The application of ArcGIS‐SWAT to the Seco Creek watershed in Texas is presented.     

HYDROLOGIC PARAMETERIZATION OF WATERSHEDS FOR RUNOFF PREDICTION USING SWAT1

ABSTRACT: The use of continuous time, distributed parameter hydrologic models like SWAT (Soil and Water Assessment Tool) has opened several opportunities to improve watershed modeling accuracy. However, it has also placed a heavy burden on users with respect to the amount of work involved in parameterizing the watershed in general and in adequately representing the spatial variability of the watershed in particular. Recent developments in Geographical Information Systems (GIS) have alleviated some of the difficulties associated with managing spatial data. However, the user must still choose among various parameterization approaches that are available within the model. This paper describes the important parameterization issues involved when modeling watershed hydrology for runoff prediction using SWAT with emphasis on how to improve model performance without resorting to tedious and arbitrary parameter by parameter calibration. Synthetic and actual watersheds in Indiana and Mississippi were used to illustrate the sensitivity of runoff prediction to spatial variability, watershed decomposition, and spatial and temporal adjustment of curve numbers and return flow contribution. SWAT was also used to predict stream runoff from actual watersheds in Indiana that have extensive subsurface drainage. The results of this study provide useful information for improving SWAT performance in terms of stream runoff prediction in a manner that is particularly useful for modeling ungaged watersheds wherein observed data for calibration is not available.     

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.