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.     

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.     

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.     

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

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

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.     

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

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

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.     

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.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

ABSTRACT: Significant land cover changes have occurred in the watersheds that contribute runoff to the upper San Pedro River in Sonora, Mexico, and southeast Arizona. These changes, observed using a series of remotely sensed images taken in the 1970s, 1980s, and 1990s, have been implicated in the alteration of the basin hydrologic response. The Cannonsville subwatershed, located in the Catskill/Delaware watershed complex that delivers water to New York City, provides a contrast in land cover change. In this region, the Cannonsville watershed condition has improved over a comparable time period. A landscape assessment tool using a geographic information system (GIS) has been developed that automates the parameterization of the Soil and Water Assessment Tool (SWAT) and KINEmatic Runoff and EROSion (KINEROS) hydrologic models. The Automated Geospatial Watershed Assessment (AGWA) tool was used to prepare parameter input files for the Upper San Pedro Basin, a subwatershed within the San Pedro undergoing significant changes, and the Cannonsville watershed using historical land cover data. Runoff and sediment yield were simulated using these models. In the Cannonsville watershed, land cover change had a beneficial impact on modeled watershed response due to the transition from agriculture to forest land cover. Simulation results for the San Pedro indicate that increasing urban and agricultural areas and the simultaneous invasion of woody plants and decline of grasslands resulted in increased annual and event runoff volumes, flashier flood response, and decreased water quality due to sediment loading. These results demonstrate the usefulness of integrating remote sensing and distributed hydrologic models through the use of GIS for assessing watershed condition and the relative impacts of land cover transitions on hydrologic response.     

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.     

Augmenting Watershed Model Calibration with Incorporation of Ancillary Data Sources and Qualitative Soft Data Sources

Watershed simulation models such as the Soil & Water Assessment Tool (SWAT) can be calibrated using “hard data” such as temporal streamflow observations; however, users may find upon examination of model outputs, that the calibrated models may not reflect actual watershed behavior. Thus, it is often advantageous to use “soft data” (i.e., qualitative knowledge such as expected denitrification rates that observed time series do not typically exist) to ensure that the calibrated model is representative of the real world. The primary objective of this study is to evaluate the efficacy of coupling SWAT‐Check (a post‐evaluation framework for SWAT outputs) and IPEAT‐SD (Integrated Parameter Estimation and Uncertainty Analysis Tool‐Soft & hard Data evaluation) to constrain the bounds of soft data during SWAT auto‐calibration. IPEAT‐SD integrates 59 soft data variables to ensure SWAT does not violate physical processes known to occur in watersheds. IPEAT‐SD was evaluated for two case studies where soft data such as denitrification rate, nitrate attributed from subsurface flow to total discharge ratio, and total sediment loading were used to conduct model calibration. Results indicated that SWAT model outputs may not satisfy reasonable soft data responses without providing pre‐defined bounds. IPEAT‐SD provides an efficient and rigorous framework for users to conduct future studies while considering both soft data and traditional hard information measures in watershed modeling.     

Targeted Application of Seasonal Load Duration Curves Using Multivariate Analysis in Two Watersheds Flowing Into Lake Houston1

Teague, Aarin, Philip B. Bedient, and Birnur Guven, 2011. Targeted Application of Seasonal Load Duration Curves Using Multivariate Analysis in Two Watersheds Flowing Into Lake Houston. Journal of the American Water Resources Association (JAWRA) 47(3):620‐634. DOI: 10.1111/j.1752‐1688.2011.00529.x Abstract: Water quality is a problem in Lake Houston, the primary source of drinking water for the City of Houston, Texas, due to pollutant loads coming from the influent watersheds, including Spring Creek and Cypress Creek. Statistical analysis of the historic water quality data was developed to understand the source characterization and seasonality of the watershed. Multivariate analysis including principal component, cluster, and discriminant analysis provided a custom seasonal assessment of the watersheds so that loading curves may be targeted for season specific pollutant source characterization. The load duration curves have been analyzed using data collected by the U.S. Geologic Survey with corresponding City of Houston water quality data at the sites to characterize the behavior of the pollutant sources and watersheds. Custom seasons were determined for Spring Creek and Cypress Creek watersheds and pollutant source characterization compared between the seasons and watersheds.     

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.     

Calibration and Verification of SWMM for Low Impact Development

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

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.     

Investigation of the Curve Number Method For Surface Runoff Estimation In Tropical Regions

This study tests the applicability of the curve number (CN) method within the Soil and Water Assessment Tool (SWAT) to estimate surface runoff at the watershed scale in tropical regions. To do this, surface runoff simulated using the CN method was compared with observed runoff in numerous rainfall‐runoff events in three small tropical watersheds located in the Upper Blue Nile basin, Ethiopia. The CN method generally performed well in simulating surface runoff in the studied watersheds (Nash‐Sutcliff efficiency [NSE] > 0.7; percent bias [PBIAS] < 32%). Moreover, there was no difference in the performance of the CN method in simulating surface runoff under low and high antecedent rainfall (PBIAS for both antecedent conditions: ~30%; modified NSE: ~0.4). It was also found that the method accurately estimated surface runoff at high rainfall intensity (e.g., PBIAS < 15%); however, at low rainfall intensity, the CN method repeatedly underestimated surface runoff (e.g., PBIAS > 60%). This was possibly due to low infiltrability and valley bottom saturated areas typical of many tropical soils, indicating that there is scope for further improvements in the parameterization/representation of tropical soils in the CN method for runoff estimation, to capture low rainfall‐intensity events. In this study the retention parameter was linked to the soil moisture content, which seems to be an appropriate approach to account for antecedent wetness conditions in the tropics.     

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.     

The Influence of Forest Regrowth on the Stream Discharge in the North Carolina Piedmont Watersheds

This study focuses on the relationships of watershed runoff with historical land use/land cover (LULC) and climate trends. Over the 20th Century, LULC in the Southeast United States, particularly the North Carolina Piedmont, has evolved from an agriculture dominated to an extensively forested landscape with more recent localized urbanization. The regrowth of forest has an important influence on the hydrology of the region as it enhances ecosystem interaction with recent climate change. During 1920‐2009, the amount of precipitation in some parts of the North Carolina Piedmont forest regrowth area showed increasing trends without corresponding increments in runoff. We employed the Soil and Water Assessment Tool (SWAT) to backcast long‐term hydrologic behavior of watersheds in North Carolina with different LULC conditions: (1) LULC conversion from agricultural to forested area and (2) long‐term stable forested area. Comparing U.S. Geological Survey‐measured stream discharge with SWAT‐simulated stream discharge under the assumption of constant 2006 LULC, we found significant stream discharge underprediction by SWAT in two LULC conversion watersheds during the early simulation period (1920s) with differences gradually decreasing by the mid‐1970s. This model bias suggests that forest regrowth on abandoned agricultural land was a key factor contributing to mitigate the impact of increased precipitation on runoff due to increasing water consumption driven by changes in vegetation.