Assessment of Bioenergy Cropping Scenarios for the Boone River Watershed in North Central Iowa,United States

Several biofuel cropping scenarios were evaluated with an improved version of Soil and Water Assessment Tool (SWAT) as part of the CenUSA Bioenergy consortium for the Boone River Watershed (BRW), which drains about 2,370 km² in north central Iowa. The adoption of corn stover removal, switchgrass, and/or Miscanthus biofuel cropping systems was simulated to assess the impact of cellulosic biofuel production on pollutant losses. The stover removal results indicate removal of 20 or 50% of corn stover in the BRW would have negligible effects on streamflow and relatively minor or negligible effects on sediment and nutrient losses, even on higher sloped cropland. Complete cropland conversion into switchgrass or Miscanthus, resulted in reductions of streamflow, sediment, nitrate, and other pollutants ranging between 23‐99%. The predicted nitrate reductions due to Miscanthus adoption were over two times greater compared to switchgrass, with the largest impacts occurring for tile‐drained cropland. Targeting of switchgrass or Miscanthus on cropland ≥2% slope or ≥7% slope revealed a disproportionate amount of sediment and sediment‐bound nutrient reductions could be obtained by protecting these relatively small areas of higher sloped cropland. Overall, the results indicate that all biofuel cropping systems could be effectively implemented in the BRW, with the most robust approach being corn stover removal adopted on tile‐drained cropland in combination with a perennial biofuel crop on higher sloped landscapes. 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.     

Influence of Bioenergy Crop Production and Climate Change on Ecosystem Services

Land use change can significantly affect the provision of ecosystem services and the effects could be exacerbated by projected climate change. We quantify ecosystem services of bioenergy‐based land use change and estimate the potential changes of ecosystem services due to climate change projections. We considered 17 bioenergy‐based scenarios with Miscanthus, switchgrass, and corn stover as candidate bioenergy feedstock. Soil and Water Assessment Tool simulations of biomass/grain yield, hydrology, and water quality were used to quantify ecosystem services freshwater provision (FWPI), food (FPI) and fuel provision, erosion regulation (ERI), and flood regulation (FRI). Nine climate projections from Coupled Model Intercomparison Project phase‐3 were used to quantify the potential climate change variability. Overall, ecosystem services of heavily row cropped Wildcat Creek watershed were lower than St. Joseph River watershed which had more forested and perennial pasture lands. The provision of ecosystem services for both study watersheds were improved with bioenergy production scenarios. Miscanthus in marginal lands of Wildcat Creek (9% of total area) increased FWPI by 27% and ERI by 14% and decreased FPI by 12% from the baseline. For St. Joseph watershed, Miscanthus in marginal lands (18% of total area) improved FWPI by 87% and ERI by 23% while decreasing FPI by 46%. The relative impacts of land use change were considerably larger than climate change impacts in this paper. 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.     

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.     

Reducing Nitrogen Export from the Corn Belt to the Gulf of Mexico: Agricultural Strategies for Remediating Hypoxia

SPAtially Referenced Regression on Watershed models developed for the Upper Midwest were used to help evaluate the nitrogen‐load reductions likely to be achieved by a variety of agricultural conservation practices in the Upper Mississippi‐Ohio River Basin (UMORB) and to compare these reductions to the 45% nitrogen‐load reduction proposed to remediate hypoxia in the Gulf of Mexico (GoM). Our results indicate that nitrogen‐management practices (improved fertilizer management and cover crops) fall short of achieving this goal, even if adopted on all cropland in the region. The goal of a 45% decrease in loads to the GoM can only be achieved through the coupling of nitrogen‐management practices with innovative nitrogen‐removal practices such as tile‐drainage treatment wetlands, drainage–ditch enhancements, stream‐channel restoration, and floodplain reconnection. Combining nitrogen‐management practices with nitrogen‐removal practices can dramatically reduce nutrient export from agricultural landscapes while minimizing impacts to agricultural production. With this approach, it may be possible to meet the 45% nutrient reduction goal while converting less than 1% of cropland in the UMORB to nitrogen‐removal practices. Conservationists, policy makers, and agricultural producers seeking a workable strategy to reduce nitrogen export from the Corn Belt will need to consider a combination of nitrogen‐management practices at the field scale and diverse nitrogen‐removal practices at the landscape scale.     

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 novel framework to classify marginal land for sustainable biomass feedstock production

To achieve food and energy security, sustainable bioenergy has become an important goal for many countries. The use of marginal lands to produce energy crops is one strategy for achieving this goal, but what is marginal land? Current definitions generally focus on a single criterion, primarily agroeconomic profitability. Herein, we present a framework that incorporates multiple criteria including profitability of current land use, soil health indicators (erosion, flooding, drainage, or high slopes), and environmental degradation resulting from contamination of surface water or groundwater resources. We tested this framework for classifying marginal land in the state of Nebraska and estimated the potential for using marginal land to produce biofuel crops. Our results indicate that approximately 1.6 million ha, or 4 million acres, of land (approximately 8% of total land area) could be classified as marginal on the basis of at least two criteria. Second-generation lignocellulosic bioenergy crops such as switchgrass ( Panicum virgatum L.), miscanthus (Miscanthus giganteus), native prairie grasses, and short-rotation woody crops could be grown on this land in redesigned landscapes that meet energy and environmental needs, without significant impacts on food or feed production. Calculating tradeoffs between the economics of redesigned landscapes and current practices at the field scale is the next step for determining functional designs for integrating biofuel feedstock production into current land management practices.     

The myth of nitrogen fertilization for soil carbon sequestration

Intensive use of N fertilizers in modern agriculture is motivated by the economic value of high grain yields and is generally perceived to sequester soil organic C by increasing the input of crop residues. This perception is at odds with a century of soil organic C data reported herein for the Morrow Plots, the world's oldest experimental site under continuous corn (Zea mays L.). After 40 to 50 yr of synthetic fertilization that exceeded grain N removal by 60 to 190%, a net decline occurred in soil C despite increasingly massive residue C incorporation, the decline being more extensive for a corn-soybean (Glycine max L. Merr.) or corn-oats (Avena sativa L.)-hay rotation than for continuous corn and of greater intensity for the profile (0-46 cm) than the surface soil. These findings implicate fertilizer N in promoting the decomposition of crop residues and soil organic matter and are consistent with data from numerous cropping experiments involving synthetic N fertilization in the USA Corn Belt and elsewhere, although not with the interpretation usually provided. There are important implications for soil C sequestration because the yield-based input of fertilizer N has commonly exceeded grain N removal for corn production on fertile soils since the 1960s. To mitigate the ongoing consequences of soil deterioration, atmospheric CO(2) enrichment, and NO(3)(-) pollution of ground and surface waters, N fertilization should be managed by site-specific assessment of soil N availability. Current fertilizer N management practices, if combined with corn stover removal for bioenergy production, exacerbate soil C loss.     

Modeling Water Quantity and Sulfate Concentrations in the Devils Lake Watershed Using Coupled SWAT and CE‐QUAL‐W2

Devils Lake is an endorheic lake in the Red River of the North basin in northeastern North Dakota. During the last two decades, the lake water level has risen by nearly 10 m, causing floods that have cost more than 1 billion USD in mitigation measures. Another increase of approximately 1.5 m in the lake water level would cause spillage into the Sheyenne River. To alleviate this potentially catastrophic spillage, two artificial outlets were constructed. However, the artificial drainage of water into the Sheyenne River raises water quality concerns because the Devils Lake water contains significantly higher concentrations of dissolved solids, particularly sulfate. In this study, the Soil and Water Assessment Tool (SWAT) was coupled with the CE‐QUAL‐W2 model to simulate both water balance and sulfate concentrations in the lake. The SWAT model performed well in simulating daily flow in tributaries with E_NS > 0.5 and |PBIAS| < 25%, and reproduced the lake water level with a root mean square error of 0.35 m for the study period from 1995 to 2014. The water temperature and sulfate concentrations simulated by CE‐QUAL‐W2 for the lake are in general agreement with the field observations. The model results show that the operation of the two outlets since August 2005 has lowered the lake level by 0.70 m. Furthermore, the models show pumping water from the two outlets raises sulfate concentrations in the Sheyenne River from ~100 to >500 mg/L. 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.     

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.     

Continental‐Scale River Flow Modeling of the Mississippi River Basin Using High‐Resolution NHDPlus Dataset

As a key component of the National Flood Interoperability Experiment (NFIE), this article presents the continental scale river flow modeling of the Mississippi River Basin (MRB), using high‐resolution river data from NHDPlus. The Routing Application for Parallel computatIon of Discharge (RAPID) was applied to the MRB with more than 1.2 million river reaches for a 10‐year study (2005‐2014). Runoff data from the Variable Infiltration Capacity (VIC) model was used as input to RAPID. This article investigates the effect of topography on RAPID performance, the differences between the VIC‐RAPID streamflow simulations in the HUC‐2 regions of the MRB, and the impact of major dams on the streamflow simulations. The model performance improved when initial parameter values, especially the Muskingum K parameter, were estimated by taking topography into account. The statistical summary indicates the RAPID model performs better in the Ohio and Tennessee Regions and the Upper and Lower Mississippi River Regions in comparison to the western part of the MRB, due to the better performance of the VIC model. The model accuracy also increases when lakes and reservoirs are considered in the modeling framework. In general, results show the VIC‐RAPID streamflow simulation is satisfactory at the continental scale of the MRB.     

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.     

An Agent-Based Modeling Approach for Determining Corn Stover Removal Rate and Transboundary Effects

Bioenergy production involves different agents with potentially different objectives, and an agent’s decision often has transboundary impacts on other agents along the bioenergy value chain. Understanding and estimating the transboundary impacts is essential to portraying the interactions among the different agents and in the search for the optimal configuration of the bioenergy value chain. We develop an agent-based model to mimic the decision making by feedstock producers and feedstock-to-biofuel conversion plant operators and propose multipliers (i.e., ratios of economic values accruing to different segments and associated agents in the value chain) for assessing the transboundary impacts. Our approach is generic and thus applicable to a variety of bioenergy production systems at different sites and geographic scales. We apply it to the case of producing ethanol using corn stover in Iowa, USA. The results from the case study indicate that stover removal rate is site specific and varies considerably with soil type, as well as other factors, such as stover price and harvesting cost. In addition, ethanol production using corn stover in the study region would have strong positive ripple effects, with the values of multipliers varying with greenhouse gas price and national energy security premium. The relatively high multiplier values suggest that a large portion of the value associated with corn stover ethanol production would accrue to the downstream end of the value chain instead of stover producers.     

Denitrification in the shallow ground water of a tile-drained, agricultural watershed

Nonpoint-source pollution of surface water by N is considered a major cause of hypoxia. Because Corn Belt watersheds have been identified as major sources of N in the Mississippi River basin, the fate and transport of N from midwestern agricultural watersheds have received considerable interest. The fate and transport of N in the shallow ground water of these watersheds still needs additional research. Our purpose was to estimate denitrification in the shallow ground water of a tile-drained, Corn Belt watershed with fine-grained soils. Over a 3-yr period, N was monitored in the surface and ground water of an agricultural watershed in central Illinois. A significant amount of N was transported past the tile drains and into shallow ground water. The ground water nitrate was isotopically heavier than tile drain nitrate, which can be explained by denitrification in the subsurface. Denitrifying bacteria were found at depths to 10 m throughout the watershed. Laboratory and push-pull tests showed that a significant fraction of nitrate could be denitrified rapidly. We estimated that the N denitrified in shallow ground water was equivalent to 0.3 to 6.4% of the applied N or 9 to 27% of N exported via surface water. These estimates varied by water year and peaked in a year of normal precipitation after 2 yr of below average precipitation. Three years of monitoring data indicate that shallow ground water in watersheds with fine-grained soils may be a significant N sink compared with N exported via surface water.     

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.     

AutoRAPID: A Model for Prompt Streamflow Estimation and Flood Inundation Mapping over Regional to Continental Extents

This article couples two existing models to quickly generate flow and flood‐inundation estimates at high resolutions over large spatial extents for use in emergency response situations. Input data are gridded runoff values from a climate model, which are used by the Routing Application for Parallel computatIon of Discharge (RAPID) model to simulate flow rates within a vector river network. Peak flows in each river reach are then supplied to the AutoRoute model, which produces raster flood inundation maps. The coupled tool (AutoRAPID) is tested for the June 2008 floods in the Midwest and the April‐June 2011 floods in the Mississippi Delta. RAPID was implemented from 2005 to 2014 for the entire Mississippi River Basin (1.2 million river reaches) in approximately 45 min. Discretizing a 230,000‐km² area in the Midwest and a 109,500‐km² area in the Mississippi Delta into thirty‐nine 1° by 1° tiles, AutoRoute simulated a high‐resolution (~10 m) flood inundation map in 20 min for each tile. The hydrographs simulated by RAPID are found to perform better in reaches without influences from unrepresented dams and without backwater effects. Flood inundation maps using the RAPID peak flows vary in accuracy with F‐statistic values between 38.1 and 90.9%. Better performance is observed in regions with more accurate peak flows from RAPID and moderate to high topographic relief.     

Spatial Optimization of Six Conservation Practices Using Swat in Tile‐Drained Agricultural Watersheds

Targeting of agricultural conservation practices to the most effective locations in a watershed can promote wise use of conservation funds to protect surface waters from agricultural nonpoint source pollution. A spatial optimization procedure using the Soil and Water Assessment Tool was used to target six widely used conservation practices, namely no‐tillage, cereal rye cover crops (CC), filter strips (FS), grassed waterways (GW), created wetlands, and restored prairie habitats, in two west‐central Indiana watersheds. These watersheds were small, fairly flat, extensively agricultural, and heavily subsurface tile‐drained. The targeting approach was also used to evaluate the model's representation of conservation practices in cost and water quality improvement, defined as export of total nitrogen, total phosphorus, and sediment from cropped fields. FS, GW, and habitats were the most effective at improving water quality, while CC and wetlands made the greatest water quality improvement in lands with multiple existing conservation practices. Spatial optimization resulted in similar cost‐environmental benefit tradeoff curves for each watershed, with the greatest possible water quality improvement being a reduction in total pollutant loads by approximately 60%, with nitrogen reduced by 20‐30%, phosphorus by 70%, and sediment by 80‐90%.     

SWATMOD‐Prep: Graphical User Interface for Preparing Coupled SWAT‐MODFLOW Simulations

This article presents SWATMOD‐Prep, a graphical user interface that couples a SWAT watershed model with a MODFLOW groundwater flow model. The interface is based on a recently published SWAT‐MODFLOW code that couples the models via mapping schemes. The spatial layout of SWATMOD‐Prep guides the user through the process of importing shape files (sub‐basins, hydrologic response units [HRUs], river network) from an existing SWAT model, creating a grid, performing necessary geo‐processing operations to link the models, writing out SWAT‐MODFLOW files, and running the simulation. The option of creating a new single‐layer MODFLOW model for near‐surface alluvial aquifers is available, with the user prompted to provide groundwater surface elevation (through a digital elevation model), aquifer thickness, and necessary aquifer parameter values. The option of simulating nitrate transport in the aquifer also is available, using the reactive transport model RT3D. The interface is in the public domain. It is programmed in Python, with various software packages used for geo‐processing operations (e.g., selection, intersection of rasters) and inputting/outputting data, and is written for Windows. The use of SWATMOD‐Prep is demonstrated for the Little River Experimental Watershed, Georgia. SWATMOD‐Prep and SWAT‐MODFLOW executables are available with an accompanying user's manual at: http://swat.tamu.edu/software/swat-modflow/ . The user's manual also accompanies this article as Supporting Information.     

Tillage and nutrient source effects on water quality and corn grain yield from a flat landscape

Beneficial effects of leaving residue at the soil surface are well documented for steep lands, but not for flat lands that are drained with surface inlets and tile lines. This study quantified the effects of tillage and nutrient source on tile line and surface inlet water quality under continuous corn (Zea mays L.) from relatively flat lands (<3%). Tillage treatments were either fall chisel or moldboard plow. Nutrient sources were either fall injected liquid hog manure or spring incorporated urea. The experiment was on a Webster-Canisteo clay loam (Typic Endoaquolls) at Lamberton, MN. Surface inlet runoff was analyzed for flow, total solids, NO(3)-N, NH(4)-N, dissolved P, and total P. Tile line effluent was analyzed for flow, NO(3)-N, and NH(4)-N. In four years of rainstorm and snowmelt events there were few significant differences (p < 0.10) in water quality of surface inlet or tile drainage between treatments. Residue cover minimally reduced soil erosion during both snowmelt and rainfall runoff events. There was a slight reduction in mineral N losses via surface inlets from manure treatments. There was also a slight decrease (p = 0.025) in corn grain yield from chisel-plow plots (9.7 Mg ha(-1)) compared with moldboard-plow plots (10.1 Mg ha(-1)). Chisel plowing (approximately 30% residue cover) alone is not sufficient to reduce nonpoint source sediment pollution from these poorly drained flat lands to the extent (40% reduction) desired by regulatory agencies.     

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.