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.     

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.     

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.     

Introduction to SWAT+, A Completely Restructured Version of the Soil and Water Assessment Tool

SWAT+ is a completely restructured version of the Soil and Water Assessment Tool (SWAT) that was developed to face present and future challenges in water resources modeling and management and to meet the needs of the worldwide user community. It is expected to improve code development and maintenance; support data availability, analysis, and visualization; and enhance the model's capabilities in terms of the spatial representation of elements and processes within watersheds. The most important change is the implementation of landscape units and flow and pollutant routing across the landscape. Also, SWAT+ offers more flexibility than SWAT in defining management schedules, routing constituents, and connecting managed flow systems to the natural stream network. To test the basic hydrologic function of SWAT+, it was applied to the Little River Experimental Watershed (Georgia) without enhanced overland routing and compared with previous models. SWAT+ gave similar results and inaccuracies as these models did for streamflow and water balance. Taking full advantage of the new capabilities of SWAT+ regarding watershed discretization and landscape and river interactions is expected to improve simulations in future studies. While many capabilities of SWAT have already been enhanced in SWAT+ and new capabilities have been added, the model will continue to evolve in response to advancements in scientific knowledge and the demands of the growing worldwide user community. 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.     

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.     

EVALUATION OF THE APPLICABILITY OF THE SWAT MODEL FOR COASTAL WATERSHEDS IN SOUTHEASTERN LOUISIANA1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) has been used for hydrologic analyses at various watershed scales. However, little is known about the model's performance in coastal watersheds. In this study SWAT was evaluated for its applicability in three Louisiana coastal watersheds: the Amite, Tickfaw, and Tangipahoa River watersheds. The model was calibrated with daily discharge from 1976 to 1977 and validated from 1979 to 1999 for the Amite and Tangipahoa and with daily discharge from 1979 to 1989 for the Tickfaw. Deviation of mean discharge and the Nash‐Sutcliffe model efficiency were used to evaluate model behavior. The study found that Manning's roughness coefficient for the main channel, SCS curve number, and soil evaporation compensation factor were the most sensitive parameters for these coastal watersheds. The Manning's roughness coefficient showed the greatest effect on the response time of surface runoff, suggesting the critical role of channel routing in hydrologic modeling for lowland watersheds. The SWAT model demonstrated an excellent performance, with Nash‐Sutcliffe efficiencies of 0.935, 0.940, and 0.960 for calibrations of the Amite, Tickfaw, and Tangipahoa watersheds, respectively, and of 0.851, 0.811, and 0.867 for validations. The modeling results demonstrate that SWAT is capable of simulating hydrologic processes for medium scale to large scale coastal lowland watersheds in Louisiana.     

COMPARISON OF PROCESS‐BASED AND ARTIFICIAL NEURAL NETWORK APPROACHES FOR STREAMFLOW MODELING IN AN AGRICULTURAL WATERSHED1

ABSTRACT: The performance of the Soil and Water Assessment Tool (SWAT) and artificial neural network (ANN) models in simulating hydrologic response was assessed in an agricultural watershed in southeastern Pennsylvania. All of the performance evaluation measures including Nash‐Sutcliffe coefficient of efficiency (E) and coefficient of determination (R²) suggest that the ANN monthly predictions were closer to the observed flows than the monthly predictions from the SWAT model. More specifically, monthly streamflow E and R² were 0.54 and 0.57, respectively, for the SWAT model calibration period, and 0.71 and 0.75, respectively, for the ANN model training period. For the validation period, these values were ?0.17 and 0.34 for the SWAT and 0.43 and 0.45 for the ANN model. SWAT model performance was affected by snowmelt events during winter months and by the model's inability to adequately simulate base flows. Even though this and other studies using ANN models suggest that these models provide a viable alternative approach for hydrologic and water quality modeling, ANN models in their current form are not spatially distributed watershed modeling systems. However, considering the promising performance of the simple ANN model, this study suggests that the ANN approach warrants further development to explicitly address the spatial distribution of hydrologic/water quality processes within watersheds.     

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.     

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.     

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.     

Use of SWAT to Estimate Spatial Scaling of Phosphorus Export Coefficients and Load Reductions Due to Agricultural BMPS

Phosphorus export coefficients (kg/ha/yr) from selected land covers, also called phosphorus yields, tend to get smaller as contributing areas get larger because some of the phosphorus mobilized on local fields gets trapped during transport to regional watershed outlets. Phosphorus traps include floodplains, wetlands, and lakes, which can then become impaired by eutrophication. The Sunrise River watershed in east central Minnesota, United States, has numerous lakes impaired by excess phosphorus. The Sunrise is tributary to the St. Croix River, whose much larger watershed is terminated by Lake St. Croix, also impaired by excess phosphorus. To support management of these impairments at both local and regional scales, a Soil and Water Assessment Tool (SWAT) model of the Sunrise watershed was constructed to estimate load reductions due to selected best management practices (BMPs) and to determine how phosphorus export coefficients scaled with contributing area. In this study, agricultural BMPs, including vegetated filter strips, grassed waterways, and reduction of soil‐phosphorus concentrations reduced phosphorus loads by 4‐20%, with similar percentage reductions at field and watershed spatial scales. Phosphorus export coefficients from cropland in rotation with corn, soybeans, and alfalfa decreased as a negative power function of contributing area, from an average of 2.12 kg/ha/yr at the upland field scale (~0.6 km²) to 0.63 kg/ha/yr at the major river basin scale (20,000 km²). 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.     

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.     

Evaluation of Three Watershed‐Scale Pesticide Environmental Transport and Fate Models1

Abstract: The U.S. Environmental Protection Agency (USEPA) Office of Pesticide Programs (OPP) has completed an evaluation of three watershed‐scale simulation models for potential use in Food Quality Protection Act pesticide drinking water exposure assessments. The evaluation may also guide OPP in identifying computer simulation tools that can be used in performing aquatic ecological exposure assessments. Models selected for evaluation were the Soil Water Assessment Tool (SWAT), the Nonpoint Source Model (NPSM), a modified version of the Hydrologic Simulation Program‐Fortran (HSPF), and the Pesticide Root Zone Model‐Riverine Water Quality (PRZM‐RIVWQ) model. Simulated concentrations of the pesticides atrazine, metolachlor, and trifluralin in surface water were compared with field data monitored in the Sugar Creek watershed of Indiana’s White River basin by the National Water Quality Assessment (NAWQA) program. The evaluation not only provided USEPA with experience in using watershed models for estimating pesticide concentration in flowing water but also led to the development of improved statistical techniques for assessing model accuracy. Further, it demonstrated the difficulty of representing spatially and temporally variable soil, weather, and pesticide applications with relatively infrequent, spatially fixed, point estimates. It also demonstrated the value of using monitoring and modeling as mutually supporting tools and pointed to the need to design monitoring programs that support modeling.     

Development and Testing of a Physically Based Model of Streambank Erosion for Coupling with a Basin‐Scale Hydrologic Model SWAT

A comprehensive streambank erosion model based on excess shear stress has been developed and incorporated in the hydrological model Soil and Water Assessment Tool (SWAT). It takes into account processes such as weathering, vegetative cover, and channel meanders to adjust critical and effective stresses while estimating bank erosion. The streambank erosion model was tested for performance in the Cedar Creek watershed in north‐central Texas where streambank erosion rates are high. A Rapid Geomorphic field assessment (RAP‐M) of the Cedar Creek watershed was done adopting techniques developed by the Natural Resources Conservation Service (NRCS), and the stream segments were categorized into various severity classes. Based on the RAP‐M field assessment, erosion pin sites were established at seven locations within the severely eroding streambanks of the watershed. A Monte Carlo simulation was carried out to assess the sensitivity of different parameters that control streambank erosion such as critical shear stress, erodibility, weathering depth, and weathering duration. The sensitive parameters were adjusted and the model was calibrated based on the bank erosion severity category identified by the RAP‐M field assessment. The average observed erosion rates were in the range 25‐367 mm year^?1. The SWAT model was able to reasonably predict the bank erosion rates within the range of variability observed in the field (R² = 0.90; E = 0.78). 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.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

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.     

An Open Source GIS‐Based Decision Support System for Watershed Evaluation of Best Management Practices

Economic costs, water quantity/quality benefits, and cost effectiveness of agricultural best management practices (BMPs) at a watershed scale are increasingly examined using integrated economic‐hydrologic models. However, these models are typically complex and not user‐friendly for examining the effects of various BMP scenarios. In this study, an open source geographic information system (GIS)‐based decision support system (DSS), named the watershed evaluation of BMPs (WEBs), was developed for creating BMP scenarios and simulating economic costs and water quantity/quality benefits at farm field, subbasin, and watershed scales. This DSS or WEBs interface integrated a farm economic model, the Soil and Water Assessment Tool (SWAT), and an optimization model within Whitebox Geospatial Analysis Tools (GAT), an open source GIS software. The DSS was applied to the 14.3‐km² Gully Creek watershed, a coastal watershed in southern Ontario, Canada that drains directly into Lake Huron. BMPs that were evaluated included conservation tillage, nutrient management, cover crop, and water and sediment control basins. In addition to assessing economic costs, water quantity/quality benefits, and cost effectiveness of BMPs, the DSS can be also used to examine prioritized BMP types/locations and corresponding economic and water quantity/quality tradeoffs in the study watershed based on environmental targets or budget constraints. Further developments of the DSS including interface transfer to other watersheds are also discussed. 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 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).