Modeling Changes to Streamflow,Sediment, and Nutrient Loading from Land Use Changes Due to Potential Natural Gas Development

Natural gas development using hydraulic fracturing has many potential environmental impacts, but among the most certain is the land disturbance required to build the well pads and other infrastructure required to drill and extract the gas. We used the Soil and Water Assessment Tool (SWAT) model to investigate how natural gas development could impact streamflow and sediment, total nitrogen (TN), and total phosphorous (TP) loadings in the upper Delaware River Basin (DRB), a relatively undeveloped watershed of 7,950 km² that lies above the Marcellus Shale formation. If gas development was permitted, our projections show the DRB could experience development of over 600 well pads to extract natural gas at build out, which, with supporting infrastructure (roads, gathering pipelines), could convert over 5,000 ha from existing land uses in the study area. In subbasins with development activity we found sediment, TN, and TP yields could increase by an average of 15, 0.08, and 0.03 kg/ha/yr, respectively (an increase of 2, 3, and 15%, respectively) for each one percent of subbasin land area converted into natural gas infrastructure. At the study area outlet on the Delaware River at Port Jervis, New York, we found increases in the annual average streamflow and sediment, nitrogen, and phosphorus loads of up to 0.01, 0.2, 0.2, and 1%, respectively, for a rapid development year, and 0.08, 1.3, 2.0, and 11%, respectively, for the full development scenario. 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.     

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.     

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.     

Evaluation of Dynamically Dimensioned Search Algorithm for Optimizing SWAT by Altering Sampling Distributions and Searching Range

The primary advantage of Dynamically Dimensioned Search (DDS) algorithm is that it outperforms other optimization techniques in both convergence speed and searching ability for parameter sets that satisfy statistical guidelines while requiring only one algorithm parameter (perturbation factor) in the optimization process. Conventionally, a default value of 0.2 is used as the perturbation factor, where a normal distribution is applied with mean sampling distribution of zero and variance of one. However, the perturbation factor sensitivity to the performance of DDS for watershed modeling is still unknown. The fixed‐form sampling distribution may result in finding parameters at the local scale rather than global in the sampling space. In this study, the efficiency of DDS was evaluated by altering the perturbation factor (from 0.05 to 1.00) and the selection of sampling distribution (normal and uniform) on hydrologic and water quality predictions in a lowland agricultural watershed in Texas, United States. Results show that the use of altered perturbation factor may cause variations in convergence speed or the ability to find better solutions. In addition, DDS results were found to be very sensitive to sampling distribution selections, where DDS‐N (normal distribution) outperformed DDS‐U (uniform distribution) in all case scenarios. The choice of sampling distributions could be the potential major factor to be attributed for the performance of auto‐calibration techniques for watershed simulation models.     

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

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

Modeling of Sediment Yield From Anjeni-Gauged Watershed,Ethiopia Using SWAT Model1

Setegn, Shimelis G., Bijan Dargahi, Ragahavan Srinivasan, and Assefa M. Melesse, 2010. Modeling of Sediment Yield From Anjeni-Gauged Watershed, Ethiopia Using SWAT Model. Journal of the American Water Resources Association (JAWRA) 46(3):514-526. DOI: 10.1111/j.1752-1688.2010.00431.x Abstract: The Soil and Water Assessment Tool (SWAT) was tested for prediction of sediment yield in Anjeni-gauged watershed, Ethiopia. Soil erosion and land degradation is a major problem on the Ethiopian highlands. The objectives of this study were to evaluate the performance and applicability of SWAT model in predicting monthly sediment yield and assess the impacts of subbasin delineation and slope discretization on the prediction of sediment yield. Ten years monthly meteorological, flow and sediment data were used for model calibration and validation. The annual average measured sediment yield was 24.6 tonnes/ha. The annual average simulated sediment yield was 27.8 and 29.5 tones/ha for calibration and validation periods, respectively. The study found that the observed values showed good agreement with the simulated sediment yield with Nash-Sutcliffe efficiency (NSE) = 0.81, percent bias (PBIAS) = 28%, RMSE-observations standard deviation ratio (RSR) = 0.23, and coefficient of determination (R²) = 0.86 for calibration and NSE = 0.79, PBIAS = 30%, RSR = 0.29, and R² = 0.84 for validation periods. The model can be used for further analysis of different management scenarios that could help different stakeholders to plan and implement appropriate soil and water conservation strategies.     

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.     

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.     

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.     

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

Channel Evolution Models as Predictors of Sediment Yield

This paper recounts our predictions of channel evolution of the Black Vermillion River (BVR) and sediment yields associated with the evolutionary sequence. Channel design parameters allowed for the prediction of stable channel form and coincident sediment yields. Measured erosion rates and basin‐specific bank erosion curves aided in prediction of the stream channel succession time frame. This understanding is critical in determining how and when to mitigate a myriad of instability consequences. The BVR drains approximately 1,062 km² in the glaciated region of Northeast Kansas. Once tallgrass prairie, the basin has been modified extensively for agricultural production. As such, channelization has shortened the river by nearly 26 km from pre‐European dimensions; shortening combined with the construction of numerous flow‐through structures have produced dramatic impacts on discharge and sediment dynamics. Nine stream reaches were established within three main tributaries of the BVR in 2007. Reaches averaged 490 m in length, were surveyed, and assessed for channel stability, while resurveys were conducted annually through 2010 to monitor change. This work illustrates the association of current stream state, in‐channel sediment contributions, and prediction of future erosion rates based on stream evolution informed by multiple models. Our findings suggest greater and more rapid sedimentation of a federal reservoir than has been predicted using standard sediment prediction methods.     

Assessing the Impacts of E. coli Laden Streambed Sediment on E. coli Loads over a Range of Flows and Sediment Characteristics

Understanding sediment Escherichia coli levels (i.e., pathogen indicators) and their contribution to the water column during resuspension is critical for predicting in‐stream E. coli levels and the potential risk to human health. The U.S. Environmental Protection Agency's current water quality testing strategies, however, rely on water borne E. coli concentrations to assess stream E. coli levels and identify impaired waters. In this work, we conducted a scenario analysis using a range of flows, sediment/water bacteria fractions, and particle sizes to which E. coli attach to assess the impact of E. coli in streambed sediments on water column E. coli levels. We used simple sediment transport theory to calculate the potential total E. coli concentrations in a stream with and without the resuspension process. Results clearly indicate that inclusion of resuspending sediment attached E. coli is necessary for watershed assessments and data on sediment attached E. coli concentrations is much needed. When neglecting the streambed sediment E. coli concentrations, the model predicted average E. coli loads of 10⁷ Colony Forming Units (CFU)/s; however, when streambed sediment E. coli concentrations were included in the model, the predictions ranged from 10¹⁰ to 10¹⁴ CFU/s. To evaluate the predictions, E. coli data in the streambed sediment and the water column were monitored in Squaw Creek, Iowa. Comparisons between measured and predicted E. coli loads yielded an R²‐value of 0.85.     

Monitored and Modeled Correlations of Sediment and Nutrients with Chesapeake Bay Water Clarity

This article analyzes the correlations of the observed and modeled light attenuation coefficient, Kd, with in situ total suspended solids (TSS) and chlorophyll‐a concentrations in Chesapeake Bay (CB) tidal waters, and with sediment and nutrient loads from the Chesapeake watershed. Light attenuation is closely related to in situ TSS and chlorophyll‐a concentrations, however, the strength of the correlation differs among the CB segments. There are distinct differences between saline and tidal fresh segments in the main Bay, but less distinction among saline and tidal fresh segments in the tidal tributaries. The correlation between Kd with sediment and nutrient loads is complicated by the lag times of TSS and the chlorophyll‐a responses to reductions in nutrient and sediment loads from the watershed, and also due to the diverse load sources. Three sets of model sensitivity scenarios were performed with: (1) differential sediment and nutrient loads; (2) selective sediment source types; and (3) geographically isolated inputs. The model results yield similar findings as those based on observed data and provide information regarding the effect of sediment on specific water bodies. Based on the model results a method was developed to determine sediment and nutrient load reductions needed to achieve the water clarity standards of the CB segments.     

Impact of Climate and Land Use Change on Streamflow and Sediment Yield in a Snow‐Dominated Semiarid Mountainous Watershed

This study investigates the impact of climate and land use change on the magnitude and timing of streamflow and sediment yield in a snow‐dominated mountainous watershed in Salt Lake County, Utah using a scenario approach and the Hydrological Simulation Program — FORTRAN model for the 2040s (year 2035–2044) and 2090s (year 2085–2094). The climate scenarios were statistically and dynamically downscaled from global climate models. Land use and land cover (LULC) changes were estimated in two ways — from a regional planning scenario and from a deterministic model. Results indicate the mean daily streamflow in the Jordan River watershed will increase by an amount ranging from 11.2% to 14.5% in the 2040s and from 6.8% to 15.3% in the 2090s. The respective increases in sediment load in the 2040s and 2090s is projected to be 6.7% and 39.7% in the canyons and about 7.4% to 14.2% in the Jordan valley. The historical 50th percentile timing of streamflow and sediment load is projected to be shifted earlier by three to four weeks by mid‐century and four to eight weeks by late‐century. The projected streamflow and sediment load results establish a nonlinear relationship with each other and are highly sensitive to projected climate change. The predicted changes in streamflow and sediment yield will have implications for water supply, flood control and stormwater management.     

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.     

Regional Blue and Green Water Balances and Use by Selected Crops in the U.S.

The availability of freshwater is a prerequisite for municipal development and agricultural production, especially in the arid and semiarid portions of the western United States (U.S.). Agriculture is the leading user of water in the U.S. Agricultural water use can be partitioned into green (derived from rainfall) and blue water (irrigation). Blue water can be further subdivided by source. In this research, we develop a hydrologic balance by 8‐Digit Hydrologic Unit Code using a combination of Soil and Water Assessment Tool simulations and available human water use estimates. These data are used to partition agricultural groundwater usage by sustainability and surface water usage by local source or importation. These predictions coupled with reported agricultural yield data are used to predict the virtual water contained in each ton of corn, wheat, sorghum, and soybeans produced and its source. We estimate that these four crops consume 480 km³ of green water annually and 23 km³ of blue water, 12 km³ of which is from groundwater withdrawal. Regional trends in blue water use from groundwater depletion highlight heavy usage in the High Plains, and small pockets throughout the western U.S. This information is presented to inform water resources debate by estimating the cost of agricultural production in terms of water regionally. This research illustrates the variable water content of the crops we consume and export, and the source of that water.     

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.     

Hydrological Variability and Uncertainty of Lower Missouri River Basin Under Changing Climate

The lower Missouri River Basin has experienced increasing streamflow and flooding events, with higher risk of extreme hydrologic impacts under changing climate. The newly available North American Regional Climate Change Assessment Program (NARCCAP) climate projections were used as atmospheric forcing for Soil and Water Assessment Tool (SWAT) model which runs with varying potential evapotranspiration (PET) methods to assess the hydrological change and uncertainty of 2040‐2069 over 1968‐1997. The NARCCAP temperature and precipitation predictions were refined using a bias correction method. The results show that, following the seasonal variability of precipitation, various water fluxes would increase in most seasons except the summer. Expected precipitation tends to increase in intensity with little change in frequency, triggering faster surface water concentration to form floods. The greatest streamflow increase would occur from November to February, increasing by around 10% on average. An increase of 3% occurs in the other months except for July and August in which river discharge decreases by around 2%. The climate predictions contribute more uncertainty annually, but PET algorithms gain more influence in winter or when other weather factors such as wind play a relatively more important role on evapotranspiration flux. This study predicts an even wetter environment compared to the historically very wet period, with the possibility of more flooding.