Evaluating Pre‐ and Post‐Fire Peak Discharge Predictions across Western U.S. Watersheds

This study reviews five models commonly used in post‐fire hydrologic assessments: the Rowe Countryman and Storey (RCS), United States Geological Survey (USGS) Linear Regression Equations, USDA Windows Technical Release 55 (USDA TR‐55), Wildcat5, and U.S. Army Corps of Engineers (USACE) Hydrologic Modeling System (HEC‐HMS). The models are applied to eight diverse basins in the western United States (U.S.) (Arizona, California, Colorado, Montana, and Washington) affected by wildfires and assessed by input parameters, calibration methods, model constraints, and performance. No one model is versatile enough for application to all study sites. Results show inconsistency between model predictions for events across the sites and less confidence with larger return periods (25‐ and 50‐year events) and post‐fire predictions. The RCS method performs well, but application is limited to southern California. The USGS linear regression model has wider regional application, but performance is less reliable at the large recurrence intervals and post‐fire predictions are reliant on a subjective modifier. Of the three curve number‐based models, Wildcat5 performs best overall without calibration, whereas the calibrated TR‐55 and HEC‐HMS models show significant improvement in pre‐fire predictions. Results from our study provide information and guidance to ultimately improve model selection for post‐fire prediction and encourage uniform parameter acquisition and calibration across the western U.S.     

Modeling Postfire Response and Recovery using the Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS)1

Abstract: This paper investigates application of the Army Corps of Engineers’ Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS) to a burned watershed in San Bernardino County, California. We evaluate the HEC‐HMS’ ability to simulate discharge in prefire and postfire conditions in a semi arid watershed and the necessary parameterizations for modeling hydrologic response during the immediate, and subsequent recovery, period after a wildfire. The model is applied to City Creek watershed, which was 90% burned during the Old Fire of October 2003. An optimal spatial resolution for the HEC‐HMS model was chosen based on an initial sensitivity analysis of subbasin configurations and related model performance. Five prefire storms were calibrated for the selected model resolution, defining a set of parameters that reasonably simulate prefire conditions. Six postfire storms, two from each of the following rainy (winter) seasons were then selected to simulate postfire response and evaluate relative changes in parameter values and model behavior. There were clear trends in the postfire parameters [initial abstractions (I_a), curve number (CN), and lag time] that reveal significant (and expected) changes in watershed behavior. CN returns to prefire (baseline) values by the end of Year 2, while I_a approaches baseline by the end of the third rainy season. However, lag time remains significantly lower than prefire values throughout the three‐year study period. Our results indicate that recovery of soil conditions and related runoff response is not entirely evidenced by the end of the study period (three rainy seasons postfire). Understanding the evolution of the land surface and related hydrologic properties during the highly dynamic postfire period, and accounting for these changes in model parameterizations, will allow for more accurate and reliable discharge simulations in both the immediate, and subsequent, rainy seasons following fire.     

Effects of Watershed Impervious Cover on Dissolved Silica Loading in Storm Flow1

Abstract: Dissolved silica (DSi) availability is a factor that affects the composition of algal populations in aquatic ecosystems. DSi cycling is tightly linked to the hydrological cycle, which is affected by human alterations of the landscape. Development activities that increase impervious cover change watershed hydrology and may increase the discharge of DSi‐poor rainwater and decrease the discharge of DSi‐rich ground water into aquatic ecosystems, possibly shifting algal community composition toward less desirable assemblages. In this study, DSi loadings from two adjacent coastal watersheds with different percent impervious cover were compared during four rain and five nonrain events. Loadings in the more impervious watershed contained a significantly larger proportion of surface runoff than base flow (ground‐water discharge) and had lower [DSi] water during rain events than the less impervious watershed. Application of the Soil Conservation Service Curve Number (CN) method showed that the minimum rainfall height necessary to yield runoff was significantly lower for the more impervious watershed, implying that runoff volumes increase with impervious cover as well as the frequency of runoff‐yielding events. Empirical data collected during this study and estimates derived from the CN method suggest that impervious cover may be responsible for both short‐term DSi limitation during rain events as well as long‐term reduction of DSi inputs into aquatic ecosystems.     

GIS‐Based Predictive Models of Hillslope Runoff Generation Processes1

Abstract: Successful nonpoint source pollution control using best management practice placement is a complex process that requires in‐depth knowledge of the locations of runoff source areas in a watershed. Currently, very few simulation tools are capable of identifying critical runoff source areas on hillslopes and those available are not directly applicable under all runoff conditions. In this paper, a comparison of two geographic information system (GIS)‐based approaches: a topographic index model and a likelihood indicator model is presented, in predicting likely locations of saturation excess and infiltration excess runoff source areas in a hillslope of the Savoy Experimental Watershed located in northwest Arkansas. Based on intensive data collected from a two‐year field study, the spatial distributions of hydrologic variables were processed using GIS software to develop the models. The likelihood indicator model was used to produce probability surfaces that indicated the likelihood of location of both saturation and infiltration excess runoff mechanisms on the hillslope. Overall accuracies of the likelihood indicator model predictions varied between 81 and 87% for the infiltration excess and saturation excess runoff locations respectively. On the basis of accuracy of prediction, the likelihood indicator models were found to be superior (accuracy 81‐87%) to the predications made by the topographic index model (accuracy 69.5%). By combining statistics with GIS, runoff source areas on a hillslope can be identified by incorporating easily determined hydrologic measurements (such as bulk density, porosity, slope, depth to bed rock, depth to water table) and could serve as a watershed management tool for identifying critical runoff source areas in locations where the topographic index or other similar methods do not provide reliable results.     

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.     

A LONG‐TERM,WATERSHED‐SCALE,EVALUATION OF THE IMPACTS OF ANIMAL WASTE BMPs ON INDICATOR BACTERIA CONCENTRATIONS1

ABSTRACT: Driven by increasing concerns about bacterial pollution from agricultural sources, states such as Virginia have initiated cost sharing programs that encourage the use of animal waste best management practices (BMPs) to control this pollution. Although a few studies have shown that waste management BMPs are effective at the field scale, their effectiveness at the watershed scale and over the long term is unknown. The focus of this research was to evaluate the effectiveness of BMPs in reducing bacterial pollution at the watershed scale and over the long term. To accomplish this goal, a 1,163 ha watershed located in the Piedmont region of Virginia was monitored over a ten‐year period. Fecal coliforms (FC) and fecal streptococci (FS) were measured as indicators of bacterial pollution. A pre‐BMP versus post‐BMP design was adopted. Major BMPs implemented were manure storage facilities, stream fencing, water troughs, and nutrient management. Seasonal Kendall trend analysis revealed a significant decreasing trend during the post‐BMP period for FC concentrations at the watershed outlet, but not at the subwatershed level. Implementation of BMPs also resulted in a significant reduction in the geometric mean of FS concentrations. FC concentrations in streamflow at the watershed outlet exceeded the Virginia primary standard 86 and 74 percent of the time during pre‐BMP and post‐BMP periods, respectively. Corresponding exceedances for the secondary standard were 50 and 41 percent. Violations decreased only slightly during the post‐BMP period. The findings of this study suggest that although BMP implementation can be expected to accomplish some improvement in water quality, BMP implementation alone may not ensure compliance with current water quality standards.     

THE IMPACT OF RUNOFF GENERATION MECHANISMS ON THE LOCATION OF CRITICAL SOURCE AREAS1

ABSTRACT: Identifying phosphorus (P) source areas and transport pathways is a key step in decreasing P loading to natural water systems. This study compared the effects of two modeled runoff generation processes ‐ saturation excess and infiltration excess ‐ on total phosphorus (TP) and soluble reactive phosphorus (SRP) concentrations in 10 catchment streams of a Catskill mountain watershed in southeastern New York. The spatial distribution of runoff from forested land and agricultural land was generated for both runoff processes; results of both distributions were consistent with Soil Conservation Service‐Curve Number (SCS‐CN) theory. These spatial runoff distributions were then used to simulate stream concentrations of TP and SRP through a simple equation derived from an observed relation between P concentration and land use; empirical results indicate that TP and SRP concentrations increased with increasing percentage of agricultural land. Simulated TP and SRP stream concentrations predicted for the 10 catchments were strongly affected by the assumed runoff mechanism. The modeled TP and SRP concentrations produced by saturation excess distribution averaged 31 percent higher and 42 percent higher, respectively, than those produced by the infiltration excess distribution. Misrepresenting the primary runoff mechanism could not only produce erroneous concentrations, it could fail to correctly locate critical source areas for implementation of best management practices. Thus, identification of the primary runoff mechanism is critical in selection of appropriate models in the mitigation of nonpoint source pollution. Correct representation of runoff processes is also critical in the future development of biogeochemical transport models, especially those that address nutrient fluxes.     

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.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

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.     

STORM RUNOFF PREDICTION BASED ON A SPATIALLY DISTRIBUTED TRAVEL TIME METHOD UTILIZING REMOTE SENSING AND GIS1

ABSTRACT: In this study, remotely sensed data and geographic information system (GIS) tools were used to estimate storm runoff response for Simms Creek watershed in the Etonia basin in northeast Florida. Land cover information from digital orthophoto quarter quadrangles (DOQQ), and enhanced thematic mapper plus (ETM+) were analyzed for the years 1990, 1995, and 2000. The corresponding infiltration excess runoff response of the study area was estimated using the U.S. Department of Agriculture (USDA), Natural Resources Conservation Service Curve Number (NRCS‐CN) method. A digital elevation model (DEM)/GIS technique was developed to predict stream response to runoff events based on the travel time from each grid cell to the watershed outlet. A comparison of predicted to observed stream response shows that the model predicts the total runoff volume with an efficiency of 0.98, the peak flow rate at an efficiency of 0.85, and the full direct runoff hydrograph with an average efficiency of 0.65. The DEM/GIS travel time model can be used to predict the runoff response of ungaged watersheds and is useful for predicting runoff hydrographs resulting from proposed large scale changes in the land use.     

PARAMETER IDENTIFICATION IN A TWO‐MULTIPLIER SEDIMENT YIELD MODEL1

ABSTRACT: A process based, distributed runoff erosion model (KINEROS2) was used to examine problems of parameter identification of sediment entrainment equations for small watersheds. Two multipliers were used to reflect the distributed nature of the sediment entrainment parameters: one multiplier for a raindrop induced entrainment parameter, and one multiplier for a flow induced entrainment parameter. The study was conducted in three parts. First, parameter identification was studied for simulated error free data sets where the parameter values were known. Second, the number of data points in the simulated sedigraphs was reduced to reflect the effect of temporal sampling frequency on parameter identification. Finally, event data from a small range‐land watershed were used to examine parameter identifiability when the parameter values are unknown. Results demonstrated that whereas unique multiplier values can be obtained for simulated error free data, unique parameter values could not be obtained for some event data. Unique multiplier values for raindrop induced entrainment and flow induced entrainment were found for events with greater than a two‐year return period (～25 mm) that also had at least 10 mm of rain in ten minutes. It was also found that the three‐minute sampling frequency used for the sediment sampler might be inadequate to identify parameters in some cases.     

EFFECTS OF WATERSHED REPRESENTATION ON RUNOFF AND SEDIMENT YIELD MODELING1

ABSTRACT: This paper evaluates the effects of watershed geometric representation (i.e., plane and channel representation) on runoff and sediment yield simulations in a semiarid rangeland watershed. A process based, spatially distributed runoff erosion model (KINEROS2) was used to explore four spatial representations of a 4.4 ha experimental watershed. The most complex representation included all 96 channel elements identifiable in the field. The least complex representation contained only five channel elements. It was concluded that oversimplified watershed representations greatly influence runoff and sediment yield simulations by inducing excessive infiltration on hillslopes and distorting runoff patterns and sediment fluxes. Runoff and sediment yield decrease systematically with decreasing complexity in watershed representation. However, less complex representations had less impact on runoff and sediment‐yield simulations for small rainfall events. This study concludes that the selection of the appropriate level of watershed representation can have important theoretical and practical implications on runoff and sediment yield modeling in semiarid environments.     

Mechanistic Simulation of Tree Effects in an Urban Water Balance Model1

Abstract: A semidistributed, physical‐based Urban Forest Effects – Hydrology (UFORE‐Hydro) model was created to simulate and study tree effects on urban hydrology and guide management of urban runoff at the catchment scale. The model simulates hydrological processes of precipitation, interception, evaporation, infiltration, and runoff using data inputs of weather, elevation, and land cover along with nine channel, soil, and vegetation parameters. Weather data are pre‐processed by UFORE using Penman‐Monteith equations to provide potential evaporation terms for open water and vegetation. Canopy interception algorithms modified established routines to better account for variable density urban trees, short vegetation, and seasonal growth phenology. Actual evaporation algorithms allocate potential energy between leaf surface storage and transpiration from soil storage. Infiltration algorithms use a variable rain rate Green‐Ampt formulation and handle both infiltration excess and saturation excess ponding and runoff. Stream discharge is the sum of surface runoff and TOPMODEL‐based subsurface flow equations. Automated calibration routines that use observed discharge has been coupled to the model. Once calibrated, the model can examine how alternative tree management schemes impact urban runoff. UFORE‐Hydro model testing in the urban Dead Run catchment of Baltimore, Maryland, illustrated how trees significantly reduce runoff for low intensity and short duration precipitation events.     

Spatial and Temporal Evaluation of Hydrological Response to Climate and Land Use Change in Three South Dakota Watersheds

This study analyzed changes in hydrology between two recent decades (1980s and 2010s) with the Soil and Water Assessment Tool (SWAT) in three representative watersheds in South Dakota: Bad River, Skunk Creek, and Upper Big Sioux River watersheds. Two SWAT models were created over two discrete time periods (1981‐1990 and 2005‐2014) for each watershed. National Land Cover Datasets 1992 and 2011 were, respectively, ingested into 1981‐1990 and 2005‐2014 models, along with corresponding weather data, to enable comparison of annual and seasonal runoff, soil water content, evapotranspiration (ET), water yield, and percolation between these two decades. Simulation results based on the calibrated models showed that surface runoff, soil water content, water yield, and percolation increased in all three watersheds. Elevated ET was also apparent, except in Skunk Creek watershed. Differences in annual water balance components appeared to follow changes in land use more closely than variation in precipitation amounts, although seasonal variation in precipitation was reflected in seasonal surface runoff. Subbasin‐scale spatial analyses revealed noticeable increases in water balance components mostly in downstream parts of Bad River and Skunk Creek watersheds, and the western part of Upper Big Sioux River watershed. Results presented in this study provide some insight into recent changes in hydrological processes in South Dakota watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Discussion

Epps et al. (2013) derived Curve Number (CN) values for two forested headwater watersheds in the Lower Coastal Plain of South Carolina during the 2008‐2011 period from rainfall‐runoff data, resulting in 23 events for the Upper Debidue Creek (UDC) watershed and in 20 events for Watershed 80 (WS80). D'Asaro and Grillone analyzed the P, CN data of the UDC watershed finding an evident “complacent” behavior, characterized by a declining CN with increasing P but without approaching a stable value at large storms. In this case, the traditional runoff CN equation does not fit well with the rainfall‐runoff data that indicate a partial source area watershed behavior and are more aptly modeled by the equation introduced by D'Asaro and Grillone (2012), who introduced a C parameter in the well‐known runoff CN formula. The C value, that represents the source area (fraction of drainage area) of the basin that produces runoff with a fixed CN < 100, is here considered to vary with the water table elevation (WTE), to which it is well correlated.     

Impact of Changing Climate and Land Cover on Flood Magnitudes in the Delaware River Basin,USA

Changing climate and land cover are expected to impact flood hydrology in the Delaware River Basin over the 21st Century. HEC‐HMS models (U.S. Army Corps of Engineers Hydrologic Engineering Center‐Hydrologic Modeling System) were developed for five case study watersheds selected to represent a range of scale, soil types, climate, and land cover. Model results indicate that climate change alone could affect peak flood discharges by ?6% to +58% a wide range that reflects regional variation in projected rainfall and snowmelt and local watershed conditions. Land cover changes could increase peak flood discharges up to 10% in four of the five watersheds. In those watersheds, the combination of climate and land cover change increase modeled peak flood discharges by up to 66% and runoff volumes by up to 44%. Precipitation projections are a key source of uncertainty, but there is a high likelihood of greater precipitation falling on a more urbanized landscape that produces larger floods. The influence of climate and land cover changes on flood hydrology for the modeled watersheds varies according to future time period, climate scenario, watershed land cover and soil conditions, and flood frequency. The impacts of climate change alone are typically greater than land cover change but there is substantial geographic variation, with urbanization the greater influence on some small, developing watersheds.     

Regional scale flood modeling using NEXRAD rainfall, GIS, and HEC-HMS/RAS: a case study for the San Antonio River Basin Summer 2002 storm event

This paper develops a framework for regional scale flood modeling that integrates NEXRAD Level III rainfall, GIS, and a hydrological model (HEC-HMS/RAS). The San Antonio River Basin (about 4000 square miles, 10,000 km2) in Central Texas, USA, is the domain of the study because it is a region subject to frequent occurrences of severe flash flooding. A major flood in the summer of 2002 is chosen as a case to examine the modeling framework. The model consists of a rainfall-runoff model (HEC-HMS) that converts precipitation excess to overland flow and channel runoff, as well as a hydraulic model (HEC-RAS) that models unsteady state flow through the river channel network based on the HEC-HMS-derived hydrographs. HEC-HMS is run on a 4 x 4 km grid in the domain, a resolution consistent with the resolution of NEXRAD rainfall taken from the local river authority. Watershed parameters are calibrated manually to produce a good simulation of discharge at 12 subbasins. With the calibrated discharge, HEC-RAS is capable of producing floodplain polygons that are comparable to the satellite imagery. The modeling framework presented in this study incorporates a portion of the recently developed GIS tool named Map to Map that has been created on a local scale and extends it to a regional scale. The results of this research will benefit future modeling efforts by providing a tool for hydrological forecasts of flooding on a regional scale. While designed for the San Antonio River Basin, this regional scale model may be used as a prototype for model applications in other areas of the country.     

Assessing Parameter Uncertainty of a Semi‐Distributed Hydrology Model for a Shallow Aquifer Dominated Environmental System

This paper examines the performance of a semi‐distributed hydrology model (i.e., Soil and Water Assessment Tool [SWAT]) using Sequential Uncertainty FItting (SUFI‐2), generalized likelihood uncertainty estimation (GLUE), parameter solution (ParaSol), and particle swarm optimization (PSO). We applied SWAT to the Waccamaw watershed, a shallow aquifer dominated Coastal Plain watershed in the Southeastern United States (U.S.). The model was calibrated (2003‐2005) and validated (2006‐2007) at two U.S. Geological Survey gaging stations, using significant parameters related to surface hydrology, hydrogeology, hydraulics, and physical properties. SWAT performed best during intervals with wet and normal antecedent conditions with varying sensitivity to effluent channel shape and characteristics. In addition, the calibration of all algorithms depended mostly on Manning's n‐value for the tributary channels as the surface friction resistance factor to generate runoff. SUFI‐2 and PSO simulated the same relative probability distribution tails to those observed at an upstream outlet, while all methods (except ParaSol) exhibited longer tails at a downstream outlet. The ParaSol model exhibited large skewness suggesting a global search algorithm was less capable of characterizing parameter uncertainty. Our findings provide insights regarding parameter sensitivity and uncertainty as well as modeling diagnostic analysis that can improve hydrologic theory and prediction in complex watersheds. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

SENSITIVITY CONSIDERATIONS WHEN MODELING HYDROLOGIC PROCESSES WITH DIGITAL ELEVATION MODEL1

ABSTRACT: The purpose of this paper is to investigate the sensitivity of a hydrologic models to the type of DEM used. This was done while modeling basin water quality with 1:24,000 and 1:250,000 U.S. Geological Survey DEMs as input to model hydro‐logical processes. The manner in which the model results were sensitive to the choice of raster cell size (scale) is investigated in this study. The Broadhead watershed, located in New Jersey, USA, was chosen as a study area. Curve numbers were estimated by a trial and error to match simulated and observed total discharge. Monthly runoff for the watershed was used in the calibration process. Higher runoff volumes were simulated by the model when the 1:24,000 DEM were used as input data, probably due to the finer resolution which simulated increased average slope and hence higher estimated runoff from the watershed. As the simulated slope of the watershed is flatten with the 1:250,000 DEM, the response of stream flow was delayed and simulated less runoff volume.