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.     

WATERSHED SCALE MODELING OF CRITICAL SOURCE AREAS OF RUNOFF GENERATION AND PHOSPHORUS TRANSPORT1

ABSTRACT: A curve number based model, Soil and Water Assessment Tool (SWAT), and a physically based model, Soil Moisture Distribution and Routing (SMDR), were applied in a headwater watershed in Pennsylvania to identify runoff generation areas, as runoff areas have been shown to be critical for phosphorus management. SWAT performed better than SMDR in simulating daily streamflows over the four‐year simulation period (Nash‐Sutcliffe coefficient: SWAT, 0.62; SMDR, 0.33). Both models varied streamflow simulations seasonally as precipitation and watershed conditions varied. However, levels of agreement between simulated and observed flows were not consistent over seasons. SMDR, a variable source area based model, needs further improvement in model formulations to simulate large peak flows as observed. SWAT simulations matched the majority of observed peak flow events. SMDR overpredicted annual flow volumes, while SWAT underpredicted the same. Neither model routes runoff over the landscape to water bodies, which is critical to surface transport of phosphorus. SMDR representation of the watershed as grids may allow targeted management of phosphorus sources. SWAT representation of fields as hydrologic response units (HRUs) does not allow such targeted management.     

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.     

Simulation of Daily Flow Pathways,Tile‐Drain Nitrate Concentrations,and Soil‐Nitrogen Dynamics Using SWAT

Tile drainage significantly alters flow and nutrient pathways and reliable simulation at this scale is needed for effective planning of nutrient reduction strategies. The Soil and Water Assessment Tool (SWAT) has been widely utilized for prediction of flow and nutrient loads, but few applications have evaluated the model's ability to simulate pathway‐specific flow components or nitrate‐nitrogen (NO₃‐N) concentrations in tile‐drained watersheds at the daily time step. The objectives of this study were to develop and calibrate SWAT models for small, tile‐drained watersheds, evaluate model performance for simulation of flow components and NO₃‐N concentration at daily intervals, and evaluate simulated soil‐nitrogen dynamics. Model evaluation revealed that it is possible to meet accepted performance criteria for simulation of monthly total flow, subsurface flow (SSF), and NO₃‐N loads while obtaining daily surface runoff (SURQ), SSF, and NO₃‐N concentrations that are not satisfactory. This limits model utility for simulating best management practices (BMPs) and compliance with water quality standards. Although SWAT simulates the soil N‐cycle and most predicted fluxes were within ranges reported in agronomic studies, improvements to algorithms for soil‐N processes are needed. Variability in N fluxes is extreme and better parameterization and constraint, through use of more detailed agronomic data, would also improve NO₃‐N simulation in SWAT. 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.     

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.     

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.     

Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho1

Stratton, Benjamin T., Venakataramana Sridhar, Molly M. Gribb, James P. McNamara, and Balaji Narasimhan, 2009. Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho. Journal of the American Water Resources Association (JAWRA) 45(6):1390‐1408. Abstract: The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed, near Boise Idaho to investigate its water balance components both temporally and spatially. Calibrating and validating SWAT is necessary to enable our understanding of the water balance components in this semiarid watershed. Daily streamflow data from four streamflow gages were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event‐based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain‐on‐snow (ROS) events during the validation period between 2005 and 2007. Especially, soil moisture was generally underestimated compared with observations from two monitoring pits. However, our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture, and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion‐related snowmelt or ROS weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting), and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.     

Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds1

Abstract: The watershed scale Soil and Water Assessment Tool (SWAT) model divides watersheds into smaller subwatersheds for simulation of rainfall‐runoff and sediment loading at the field level and routing through stream networks. Typically, the SWAT model first needs to be calibrated and validated for accurate estimation through adjustment of sensitive input parameters (i.e., Curve Number values, USLE P, slope and slope‐length, and so on). However, in some instances, SWAT‐simulated results are greatly affected by the watershed delineation and Digital Elevation Models (DEM) cell size. In this study, the SWAT ArcView GIS Patch II was developed for steep sloping watersheds, and its performance was evaluated for various threshold values and DEM cell size scenarios when delineating subwatersheds using the SWAT model. The SWAT ArcView GIS Patch II was developed using the ArcView GIS Avenue program and Spatial Analyst libraries. The SWAT ArcView GIS Patch II improves upon the SWAT ArcView GIS Patch I because it reflects the topographic factor in calculating the field slope‐length of Hydrologic Response Units in the SWAT model. The simulated sediment value for 321 subwatersheds (watershed delineation threshold value of 25 ha) is greater than that for 43 subwatersheds (watershed delineation threshold value of 200 ha) by 201% without applying the SWAT ArcView GIS Patch II. However, when the SWAT ArcView GIS Patch II was applied, the difference in simulated sediment yield decreases for the same scenario (i.e., difference in simulated sediment with 321 subwatersheds and 43 subwatersheds) was 12%. The simulated sediment value for DEM cell size of 50 m is greater than that for DEM cell size of 10 m by 19.8% without the SWAT ArcView GIS Patch II. However, the difference becomes smaller (3.4% difference) between 50 and 10 m with the SWAT ArcView GIS Patch II for the DEM scenarios. As shown in this study, the SWAT ArcView GIS Patch II can reduce differences in simulated sediment values for various watershed delineation and DEM cell size scenarios. Without the SWAT ArcView GIS Patch II, variations in the SWAT‐simulated results using various watershed delineation and DEM cell size scenarios could be greater than those from input parameter calibration. Thus, the results obtained in this study show that the SWAT ArcView GIS Patch II should be used when simulating hydrology and sediment yield for steep sloping watersheds (especially if average slope of the subwatershed is >25%) for more accurate simulation of hydrology and sediment using the SWAT model. The SWAT ArcView GIS Patch II is available at http://www.EnvSys.co.kr/~swat for free download.     

WATER QUALITY MODEL OUTPUT UNCERTAINTY AS AFFECTED BY SPATIAL RESOLUTION OF INPUT DATA1

ABSTRACT: Resolution of the input GIS data used to parameterize distributed‐parameter hydrologic/water quality models may affect uncertainty in model outputs and impact the subsequent application of model results in watershed management. In this study we evaluated the impact of varying spatial resolutions of DEM, land use, and soil data (30 × 30 m, 100 × 100 m, 150 × 150 m, 200 × 200 m, 300 × 300 m, 500 × 500 m, and 1,000 × 1,000 m) on the uncertainty of SWAT predicted flow, sediment, NO₃‐N, and TP transport. Inputs included measured hydrologic, meteorological, and watershed characteristics as well as water quality data from the Moores Creek watershed in Washington County, Arkansas. The SWAT model output was most affected by input DEM data resolution. A coarser DEM data resolution resulted in decreased representation of watershed area and slope and increased slope length. Distribution of pasture, forest, and urban areas within the watershed was significantly affected at coarser resolution of land use and resulted in significant uncertainty in predicted sediment, NO₃‐N, and TP output. Soils data resolution had no significant effect on flow and NO₃‐N predictions; however, sediment was overpredicted by 26 percent, and TP was underpredicted by 26 percent at 1,000 m resolution. This may be due to change in relative distribution of various hydrologic soils groups (HSGs) in the watershed. Minimum resolution for input GIS data to achieve less than 10 percent model output error depended upon the output variable of interest. For flow, sediment, NO₃‐N, and TP predictions, minimum DEM data resolution should range from 30 to 300 m, whereas minimum land use and soils data resolution should range from 300 to 500 m.     

Estimating Evapotranspiration for Dryland Cropping Systems in the Semiarid Texas High Plains Using SWAT

The Soil and Water Assessment Tool (SWAT) is one of the most widely used watershed models for simulating hydrology in response to agricultural management practices. However, limited studies have been performed to evaluate the SWAT model's ability to estimate daily and monthly evapotranspiration (ET) in semiarid regions. ET values were simulated using ArcSWAT 2012 for a lysimeter field managed under dryland conditions at the USDA‐ARS Conservation and Production Research Laboratory at Bushland, Texas, and compared with measured lysimeter values from 2000 to 2010. Two scenarios were performed to compare SWAT's performance: (1) use of default plant leaf area index (LAI) values in the embedded plant database and (2) adjusted LAI values. Scenario 1 resulted in an “unsatisfactory” Nash‐Sutcliffe efficiency (NSE) of 0.42 and 0.38 for the calibration and validation periods, respectively. Scenario 2 resulted in a “satisfactory” NSE value for the calibration period while achieving a “good” NSE of 0.70 for the validation period. SWAT generally underestimated ET at both the daily and monthly levels. Overestimation during fallow years may be due to the limitations of the pothole function used to simulate furrow diking. Users should be aware of potential errors associated with using default LAI parameters. Inaccuracies in ET estimation may also stem from errors in the plant stress functions, particularly when evaluating water management practices for dryland watersheds.     

EFFECTS OF STATSGO AND SSURGO AS INPUTS ON SWAT MODEL'S SNOWMELT SIMULATION1

ABSTRACT: Soil data comprise a basic input of SWAT (Soil and Water Assessment Tool) for a watershed application. For watersheds where site specific soil data are unavailable, the two U.S. Department of Agriculture (USDA) soil databases, the State Soil Geographic (STATSGO) and Soil Survey Geographic (SSURGO) databases, may be the best alternatives. Although it has been noted that SWAT models using the STATSGO and SSURGO data may give different simulation results for water, sediment, and agricultural chemical yields, information is scarce on the effects of using these two databases in predicting streamflows that are predominantly generated from melting snow in spring. The objective of this study was to assess the effects of using STATSGO versus SSURGO as an input for the SWAT model's simulation of the streamflows in the upper 45 percent of the Elm River watershed in eastern North Dakota. Designating the model as SWAT‐STATSGO when the STATSGO data were used and SWAT‐SSURGO when the SSURGO data were used, SWAT‐STATSGO and SWAT‐SSURGO were separately calibrated and validated using the observed daily streamflows. The results indicated that SWAT‐SSURGO provided an overall better prediction of the discharges than SWAT‐STATSGO, although both did a good and comparable job of predicting the high streamflows. However, SWAT‐STATSGO predicted the low streamflows more accurately and had a slightly better performance during the validation period. In addition, the discrepancies between the discharges predicted by these two SWAT models tended to be larger at upstream locations than at those farther downstream within the study area.     

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

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

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.     

Does the Temporal Resolution of Precipitation Input Influence the Simulated Hydrological Components Employing the SWAT Model?

This study aimed to evaluate the influence of sub‐daily precipitation time steps on model performance and hydrological components by applying the Green and Ampt infiltration method using the Soil and Water Assessment Tool (SWAT). Precipitation was measured at a resolution of 0.1 mm and aggregated to 5‐, 15‐, 30‐, and 60‐min time steps. Daily discharge data over a 10‐year period were used to calibrate and validate the model. Following a global sensitivity analysis, relevant parameters were optimized through an automatic calibration procedure using SWAT‐CUP for each time step. Daily performance statistics were almost equal among all four time steps (NSE ≈ 0.47). Discharge mainly consisted of groundwater flow (55%) and tile flow (42%), in reasonable proportions for the investigated catchment. In conclusion, model outputs were almost identical, showing simulations responded nearly independently of the chosen precipitation time step. This held true for (1) the selection of sensitive parameters, (2) performance statistics, (3) the shape of the hydrographs, and (4) flow components. However, a scenario analysis revealed that the precipitation time step becomes important when saturated hydraulic conductivities are low and curve numbers are high. The study suggests that there is no need in using precipitation time steps <1 h for lowland catchments dominated by soils with a low surface runoff potential if daily flow values are being considered. 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.     

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.     

Reservoir Operation with Feedback in a Coupled Land Surface and Hydrologic Model: A Case Study of the Huai River Basin,China

Regarding emerging large‐scale reservoir operation models, reports of reservoir operation feedback for hydrologic modeling are rare, and little attention has been paid to flood control. An operation scheme considering multilevel flood control (MLFC) was first proposed in this study, but more reservoir information was needed. Thus, an alternative scheme was proposed that consisted of a modified version of the reservoir operation scheme in the Soil and Water Assessment Tool Model (MSWAT scheme). These schemes were coupled to a land surface and hydrologic model system with feedback, i.e., a system in which reservoir operation can affect the subsequent simulation, and were investigated in the Huai River Basin. The results show reservoir storage and peak flow were generally overestimated by the original SWAT reservoir scheme (SWAT scheme). Compared with the SWAT scheme, the MSWAT scheme successfully reduced the simulated storage and peak flow at the reservoir stations. For the downstream stations, the streamflow simulations were improved at a significance level of 5%. The performances of the MSWAT and MLFC schemes at the reservoir stations were nearly equivalent. Importantly, reservoir operation feedback to hydrologic modeling was necessary because the reservoir operation effects could not be transferred downstream without it. The streamflow simulation of a reservoir station located on a flat plain was less sensitive to feedback than that of a mountain reservoir station.     

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%.     

Estimating Potential Climate Change Effects on the Upper Neuse Watershed Water Balance Using the SWAT Model

Climate change poses water resource challenges for many already water stressed watersheds throughout the world. One such watershed is the Upper Neuse Watershed in North Carolina, which serves as a water source for the large and growing Research Triangle Park region. The aim of this study was to quantify possible changes in the watershed’s water balance due to climate change. To do this, we used the Soil and Water Assessment Tool (SWAT) model forced with different climate scenarios for baseline, mid‐century, and end‐century time periods using five different downscaled General Circulation Models. Before running these scenarios, the SWAT model was calibrated and validated using daily streamflow records within the watershed. The study results suggest that, even under a mitigation scenario, precipitation will increase by 7.7% from the baseline to mid‐century time period and by 9.8% between the baseline and end‐century time period. Over the same periods, evapotranspiration (ET) would decrease by 5.5 and 7.6%, water yield would increase by 25.1% and 33.2%, and soil water would increase by 1.4% and 1.9%. Perhaps most importantly, the model results show, under a high emission scenario, large seasonal differences with ET estimated to decrease by up to 42% and water yield to increase by up to 157% in late summer and fall. Planning for the wetter predicted future and corresponding seasonal changes will be critical for mitigating the impacts of climate change on water resources.