Field Testing the Riparian Ecosystem Management Model on a Riparian Buffer in the North Carolina Upper Coastal Plain

The riparian ecosystem management model (REMM) was field tested using five years (2005‐2009) of measured hydrologic and water quality data on a riparian buffer located in the Tar‐Pamlico River Basin, North Carolina. The buffer site received NO₃‐N loading from an agricultural field that was fertilized with inorganic fertilizer. Field results showed the buffer reduced groundwater NO₃‐N concentration moving to the stream over a five‐year period. REMM was calibrated hydrologically using daily field‐measured water table depths (WTDs), and with monthly NO₃‐N concentrations in groundwater wells. Results showed simulated WTDs and NO₃‐N concentrations in good agreement with measured values. The mean absolute error and Willmott's index of agreement for WTDs varied from 13‐45 cm and 0.72‐0.92, respectively, while the root mean square error and Willmott's index of agreement for NO₃‐N concentrations ranged from 1.04‐5.92 mg/l and 0.1‐0.86, respectively, over the five‐year period. REMM predicted plant nitrogen (N) uptake and denitrification were within ranges reported in other riparian buffer field studies. The calibrated and validated REMM was used to simulate 33 years of buffer performance at the site. Results showed that on average the buffer reduced NO₃‐N concentrations from 12 mg/l at the field edge to 0.7 mg/l at the stream edge over the simulation period, while the total N and NO₃‐N load reductions from the field edge to the stream were 77 and 82%, respectively.     

Estimation of Daily Streamflow of Southeastern Coastal Plain Watersheds by Combining Estimated Magnitude and Sequence

Commonly used methods to predict streamflow at ungauged watersheds implicitly predict streamflow magnitude and temporal sequence concurrently. An alternative approach that has not been fully explored is the conceptualization of streamflow as a composite of two separable components of magnitude and sequence, where each component is estimated separately and then combined. Magnitude is modeled using the flow duration curve (FDC), whereas sequence is modeled by transferring streamflow sequence of gauged watershed(s). This study tests the applicability of the approach on watersheds ranging in size from about 25‐7,226 km² in Southeastern Coastal Plain (U.S.) with substantial surface storage of wetlands. A 19‐point regionalized FDC is developed to estimate streamflow magnitude using the three most selected variables (drainage area, hydrologic soil index, and maximum 24‐h precipitation with a recurrence interval of 100 years) by a greedy‐heuristic search process. The results of validation on four watersheds (Trent River, North Carolina: 02092500; Satilla River, Georgia: 02226500; Black River, South Carolina: 02136000; and Coosawhatchie River, South Carolina: 02176500) yielded Nash‐Sutcliffe efficiency values of 0.86‐0.98 for the predicted magnitude and 0.09‐0.84 for the predicted daily streamflow over a simulation period of 1960‐2010. The prediction accuracy of the method on two headwater watersheds at Santee Experimental Forest in coastal South Carolina was weak, but comparable to simulations by MIKE‐SHE.     

Characterizing Uncertainty in Daily Streamflow Estimates at Ungauged Locations for the Massachusetts Sustainable Yield Estimator

Hydrologic characterization at ungauged locations is one of the quintessential challenges of hydrology. Beyond simulation of historical streamflows, it is similarly important to characterize the level of uncertainty in hydrologic estimates. In tandem with updates to Massachusetts Sustainable Yield Estimator, this work explores the application of global uncertainty estimates to daily streamflow simulations. Expanding on a method developed for deterministic modeling, this approach produces confidence intervals on daily streamflow developed through nonlinear spatial interpolation of daily streamflow using flow duration curves; the 95% confidence is examined. Archived cross‐validations of daily streamflows from 66 watersheds in and around Massachusetts are used to evaluate an approach to uncertainty characterization. Neighboring sites are treated as ungauged, producing relative errors that can be resampled and applied to target sites. The method, with some modification, is found to provide appropriately narrow confidence intervals that contain 95% of the observed streamflows in cross‐validation. Further characterizing uncertainty, multiday means of daily streamflow are evaluated. Working through cross‐validation in Massachusetts, two‐ to three‐month averages of daily streamflow show the best performance. These two approaches to uncertainty characterization inform how streamflow simulation produced for prediction in ungauged basins can be used for water resources management.     

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.     

Spatial and Temporal Patterns of Low Streamflow and Precipitation Changes in the Chesapeake Bay Watershed

Spatial and temporal patterns in low streamflows were investigated for 183 streamgages located in the Chesapeake Bay Watershed for the period 1939–2013. Metrics that represent different aspects of the frequency and magnitude of low streamflows were examined for trends: (1) the annual time series of seven‐day average minimum streamflow, (2) the scaled average deficit at or below the 2% mean daily streamflow value relative to a base period between 1939 and 1970, and (3) the annual number of days below the 2% threshold. Trends in these statistics showed spatial cohesion, with increasing low streamflow volume at streamgages located in the northern uplands of the Chesapeake Bay Watershed and decreasing low streamflow volume at streamgages in the southern part of the watershed. For a small subset of streamgages (12%), conflicting trend patterns were observed between the seven‐day average minimum streamflow and the below‐threshold time series and these appear to be related to upstream diversions or the influence of reservoir‐influenced streamflows in their contributing watersheds. Using multivariate classification techniques, mean annual precipitation and fraction of precipitation falling as snow appear to be broad controls of increasing and decreasing low‐flow trends. Further investigation of seasonal precipitation patterns shows summer rainfall patterns, driven by the Atlantic Multidecadal Oscillation, as the main driver of low streamflows in the Chesapeake Bay Watershed.     

Streamflow Hydrology Estimate Using Machine Learning (SHEM)

Continuity and accuracy of near real‐time streamflow gauge (streamgage) data are critical for flood forecasting, assessing imminent risk, and implementing flood mitigation activities. Without these data, decision makers and first responders are limited in their ability to effectively allocate resources, implement evacuations to save lives, and reduce property losses. The Streamflow Hydrology Estimate using Machine Learning (SHEM) is a new predictive model for providing accurate and timely proxy streamflow data for inoperative streamgages. SHEM relies on machine learning (“training”) to process and interpret large volumes (“big data”) of historic complex hydrologic information. Continually updated with real‐time streamflow data, the model constructs a virtual dataset index of correlations and groups (clusters) of relationship correlations between selected streamgages in a watershed and under differing flow conditions. Using these datasets, SHEM interpolates estimated discharge and time data for any indexed streamgage that stops transmitting data. These estimates are continuously tested, scored, and revised using multiple regression analysis processes and methodologies. The SHEM model was tested in Idaho and Washington in four diverse watersheds, and the model's estimates were then compared to the actual recorded data for the same time period. Results from all watersheds revealed a high correlation, validating both the degree of accuracy and reliability of the model.     

REGIONAL APPLICATION OF AN APPROXIMATE STREAMFLOW PARTITIONING METHOD1

ABSTRACT: The approximate streamflow partitioning method which uses daily rainfall and streamfiow data was applied in Coastal Plain, Coastal Flatwoods, and Southern Piedmont physiographic regions for estimation of the surface and subsurface flow components of total streainflow. Sizes of the watersheds ranged from 9.6 km² to 1,030 km. Although the streamflow partitioning method was developed and tested on the Coastal Plain physiographic region, results indicate that the procedure can be applied to other physiographic regions where available data are limited to daily values. The effect of channelization on the partitioned flow components in the Coastal Plain and Coastal Flatwoods physiographic areas was also examined. While channelization was found to decrease the storm-time base, it had no significant effect on the relative percentages of the partitioned flow components.     

Effects of Impervious Area and BMP Implementation and Design on Storm Runoff and Water Quality in Eight Small Watersheds

The effects of increases in effective impervious area (EIA) and the implementation of water quality protection designed detention pond best management practices (BMPs) on storm runoff and stormwater quality were assessed in Gwinnett County, Georgia, for the period 2001‐2008. Trends among eight small watersheds were compared, using a time trend study design. Significant trends were detected in three storm hydrologic metrics and in five water quality constituents that were adjusted for variability in storm characteristics and climate. Trends in EIA ranged from 0.10 to 1.35, and changes in EIA treated by BMPs ranged from 0.19 to 1.32; both expressed in units of percentage of drainage area per year. Trend relations indicated that for every 1% increase in watershed EIA, about 2.6, 1.1, and 1.5% increases in EIA treated by BMPs would be required to counteract the effects of EIA added to the watersheds on peak streamflow, stormwater yield, and storm streamflow runoff, respectively. Relations between trends in EIA, BMP implementation, and water quality were counterintuitive. This may be the result of (1) changes in constituent inputs in the watersheds, especially downstream of areas treated by BMPs; (2) BMPs may have increased the duration of stormflow that results in downstream channel erosion; and/or (3) spurious relationships between increases in EIA, BMP implementation, and constituent inputs with development rates.     

MODELING ACTUAL EVAPOTRANSPIRATION FROM FORESTED WATERSHEDS ACROSS THE SOUTHEASTERN UNITED STATES1

ABSTRACT: About 50 to 80 percent of precipitation in the southeastern United States returns to the atmosphere by evapotranspiration. As evapotranspiration is a major component in the forest water balances, accurately quantifying it is critical to predicting the effects of forest management and global change on water, sediment, and nutrient yield from forested watersheds. However, direct measurement of forest evapotranspiration on a large basin or a regional scale is not possible. The objectives of this study were to develop an empirical model to estimate long‐term annual actual evapotranspiration (ART) for forested watersheds and to quantify spatial AET patterns across the southeast. A geographic information system (GIS) database including land cover, daily streamflow, and climate was developed using long term experimental and monitoring data from 39 forested watersheds across the region. Using the stepwise selection method implemented in a statistical modeling package, a long term annual AET model was constructed. The final multivariate linear model includes four independent variables—annual precipitation, watershed latitude, watershed elevation, and percentage of forest coverage. The model has an adjusted R² of 0.794 and is sufficient to predict long term annual ART for forested watersheds across the southeastern United States. The model developed by this study may be used to examine the spatial variability of water availability, estimate annual water loss from mesoscale watersheds, and project potential water yield change due to forest cover change.     

IDENTIFICATION OF CHANGES IN STREAMFLOW CHARACTERISTICS1

ABSTRACT: A set of procedures for identifying changes in selected streamflow characteristics at sites having long‐term continuous streamflow records is illustrated by using streamflow data from the Waccamaw River at Freeland, North Carolina for the 55‐year period of 1940–1994. Data were evaluated and compared to streamflow in the adjacent Lumber River Basin to determine if changes in streamflow characteristics in the Waccamaw River were localized and possibly the result of some human activity, or consistent with regional variations. Following 1963, droughts in the Waccamaw Basin seem to have been less severe than in the Lumber Basin, and the annual one‐, seven‐, and 30‐day low flows exhibited a slightly increasing trend in the Waccamaw River. Mean daily flows in the Waccamaw River at the 90 percent exceedance level (low flows) during 1985–194, a relatively dry period, were very nearly equal to flows at the same exceedance level for 1970–1979, which represents the 10‐year period between 1940 and 1994 with the highest flows. Prior to the 1980s, flows per unit drainage area in the Waccamaw Basin were generally less than those in the Lumber Basin, but after 1980, the opposite was true. The ratio of base flow to runoff in the Waccamaw River may have changed relative to that in the Lumber River in the late 1970s. There was greater variability in Waccamaw River streamflow than in Lumber River flow, and flow variability in the Waccamaw River may have increased slightly during 1985–1994.     

Impacts of Climate Change on Extreme Precipitation Events Over Flamingo Tropicana Watershed1

Abstract: Climate change, particularly the projected changes to precipitation patterns, is likely to affect runoff both regionally and temporally. Extreme rainfall events are expected to become more intense in the future in arid urban areas and this will likely lead to higher streamflow. Through hydrological modeling, this article simulates an urban basin response to the most intense storm under anthropogenic climate change conditions. This study performs an event‐based simulation for shorter duration storms in the Flamingo Tropicana (FT) watershed in Las Vegas, Nevada. An extreme storm, defined as a 100‐year return period storm, is selected from historical records and perturbed to future climatic conditions with respect to multimodel multiscenario (A1B, A2, B1) bias corrected and spatially disaggregated data from the World Climate Research Programme's (WCRP's) database. The cumulative annual precipitation for each 30‐year period shows a continuous decrease from 2011 to 2099; however, the summer convective storms, which are considered as extreme storms for the study area, are expected to be more intense in future. Extreme storm events show larger changes in streamflow under different climate scenarios and time periods. The simulated peak streamflow and total runoff volume shows an increase from 40% to more than 150% (during 2041‐2099) for different climate scenarios. This type of analysis can help evaluate the vulnerability of existing flood control system and flood control policies.     

Modeling Streamflow and Water Quality Sensitivity to Climate Change and Urban Development in 20 U.S. Watersheds

Watershed modeling in 20 large, United States (U.S.) watersheds addresses gaps in our knowledge of streamflow, nutrient (nitrogen and phosphorus), and sediment loading sensitivity to mid‐21st Century climate change and urban/residential development scenarios. Use of a consistent methodology facilitates regional scale comparisons across the study watersheds. Simulations use the Soil and Water Assessment Tool. Climate change scenarios are from the North American Regional Climate Change Assessment Program dynamically downscaled climate model output. Urban and residential development scenarios are from U.S. Environmental Protection Agency's Integrated Climate and Land Use Scenarios project. Simulations provide a plausible set of streamflow and water quality responses to mid‐21st Century climate change across the U.S. Simulated changes show a general pattern of decreasing streamflow volume in the central Rockies and Southwest, and increases on the East Coast and Northern Plains. Changes in pollutant loads follow a similar pattern but with increased variability. Ensemble mean results suggest that by the mid‐21st Century, statistically significant changes in streamflow and total suspended solids loads (relative to baseline conditions) are possible in roughly 30‐40% of study watersheds. These proportions increase to around 60% for total phosphorus and total nitrogen loads. Projected urban/residential development, and watershed responses to development, are small at the large spatial scale of modeling in this study.     

DOCUMENTED AUTUMNAL STREAMFLOW INCREASE WITHOUT MEASURABLE PRECIPITATION1

ABSTRACT: In the south-central interior of British Columbia, the summer and fall of 1987 were very dry. Streamflow was extremely low throughout the region, and an extended rainless period coincided with the onset of the winter dormant period. During a five-day rainless period from October 9 to October 13, which should have been a period of declining streamflow in all watersheds, 19 of 31 small forested watersheds having suitable natural flow record experienced increasing discharge. The data is reviewed and analyzed with an initial attempt made to explain the streamflow changes among the various watersheds quantitatively. The only explanation for the widespread increases in streamflow in the region during this dry period is that greatly reduced transpiration, as the vegetation went dormant, permitted increased volumes of water in the soil to drain into surface streams rapidly enough to be recorded as increased discharge in 61 percent of the study streams.     

EFFECTS OF BASIN-SCALE TIMBER HARVEST ON WATER YIELD AND PEAK STREAMFLOW1

ABSTRACT: Streamflow changes resulting from clearcut harvest of lodgepole pine (Pinus contorta) on a 2145 hectare drainage basin are evaluated by the paired watershed technique. Thirty years of continuous daily streamflow records were used in the analysis, including 10 pre-harvest and 20 post-harvest years of data. Regression analysis was used to estimate the effects of timber harvest on annual water yield and annual peak discharge. Removal of 14 million board feet of lodgepole pine (Pinus contorta) from about 526 hectares (25 percent of the basin) produced an average of 14.7 cm additional water yield per year, or an increase of 52 percent. Mean annual daily maximum discharge also increased by 1.6 cubic meters per second or 66 percent. Increases occurred primarily during the period of May through August with little or no change in wintertime streamflows. Results suggest that clearcutting conifers in relatively large watersheds (> 2000 ha) may produce significant increases in water yield and flooding. Implications of altered streamflow regimes are important for assessing the future ecological integrity of stream ecosystems subject to large-scale timber harvest and other disturbances that remove a substantial proportion of the forest cover.     

Impact of Length of Dataset on Streamflow Calibration Parameters and Performance of APEX Model

Due to resource constraints, long‐term monitoring data for calibration and validation of hydrologic and water quality models are rare. As a result, most models are calibrated and, if possible, validated using limited measured data. However, little research has been done to determine the impact of length of available calibration data on model parameterization and performance. The main objective of this study was to evaluate the impact of length of calibration data (LCD) on parameterization and performance of the Agricultural Policy Environmental eXtender model for predicting daily, monthly, and annual streamflow. Long‐term (1984‐2015) measured daily streamflow data from Rock Creek watershed, an agricultural watershed in northern Ohio, were used for this study. Data were divided into five Short (5‐year), two Medium (15‐year), and one Long (25‐year) streamflow calibration data scenarios. All LCD scenarios were calibrated and validated at three time steps: daily, monthly, and annual. Results showed LCD affected the ability of the model to accurately capture temporal variability in simulated streamflow. However, overall average streamflow, water budgets, and crop yields were simulated reasonably well for all LCD scenarios.     

Temporal Changes in Streamflow and Attribution of Changes to Climate and Landuse in Wisconsin Watersheds

Previous historic trends analyses on 21st Century hydrologic data in the United States generally focus on annual flow statistics and have continued to use USGS hydro‐climatic data network (HCDN) stations, although post‐1988 diversions and runoff regulations are not reflected in the HCDN. Using a more recent dataset, Geospatial Attributes of Gages for Evaluating Streamflow, version II (GAGES II), compiled by Falcone (2012), which includes more watersheds with reference conditions, a comprehensive analysis of changes in seasonal, and annual streamflow in Wisconsin watersheds is demonstrated. Given the pronounced influence of seasonal hydrology in Wisconsin watersheds, the objective of this study is to elucidate the nature of temporal (annual, seasonal, and monthly) changes in runoff. Considerable temporal and regional variability was found in annual and seasonal streamflow changes between the two historic periods 1951‐1980 and 1981‐2010 considered in the study. For example, the northern watersheds show relatively small changes in streamflow discharge ranging from ?6.0 to 4.2%, while the southern watersheds show relatively large increases in streamflow discharge ranging from 13.1 to 18.2%. To apportion streamflow changes to climate and nonclimatic factors, a method based on potential evapotranspiration changes is demonstrated. Results show that nonclimatic factors account for more than 60% of changes in annual runoff in Wisconsin watersheds considered in the study.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

Factors Affecting Sediment Oxygen Demand Dynamics in Blackwater Streams of Georgia’s Coastal Plain1

Abstract: Sediment oxygen demand (SOD) is believed to be an important process affecting dissolved oxygen (DO) concentrations in blackwater streams of the southeastern coastal plain. Because very few data on SOD are available, it is common for modelers to take SOD values from the literature for use with DO models. In this study, SOD was measured in seven blackwater streams of the Suwannee River Basin within the Georgia coastal plain for between August 2004 and April 2005. SOD was measured using four in situ chambers and was found to vary on average between 0.1 and 2.3 g O₂/m/day across the seven study sites throughout the study period. SOD was found to vary significantly between the watersheds within the Suwannee River Basin. However, land use was not found to be the driving force behind SOD values. Statistical analyses did find significant interaction between land use and watersheds suggesting that an intrinsically different factor in each of the watersheds may be affecting SOD and the low DO concentrations. Further research is needed to identify the factors driving SOD dynamics in the blackwater streams of Georgia’s coastal plain. Results from this study will be used by the Georgia Department of Natural Resources – Environmental Protection Division as model input data for the development and evaluation of DO total maximum daily loads in the Georgia coastal plain.     

WATER QUALITY IN AGRICULTURAL,URBAN, AND MIXED LAND USE WATERSHEDS1

ABSTRACT: Water quality and nonpoint source (NPS) pollution are important issues in many areas of the world, including the Inner Bluegrass Region of Kentucky where urban development is changing formerly rural watersheds into urban and mixed use watersheds. In watersheds where land use is mixed, the relative contributions of NPS pollution from rural and urban land uses can be difficult to separate. To better understand NPS pollution sources in mixed use watersheds, surface water samples were taken at three sites that varied in land use to examine the effect of land use on water quality. Within the group of three watersheds, one was predominately agriculture (Agricultural), one was predominately urban (Urban), and a third had relatively equal representation of both types of land uses (Mixed). Nitrogen (N), phosphorus (P), total suspended solids (TSS), turbidity, pH, temperature, and streamflow were measured for one year. Comparisons are made among watersheds for concentration and fluxes of water quality parameters. Nitrate and orthophosphate concentrations were found to be significantly higher in the Agricultural watershed. Total suspended solids, turbidity, temperature, and pH, were found to be generally higher in the Urban and Mixed watersheds. No differences were found for streamflow (per unit area), total phosphorus, and ammonium concentrations among watersheds. Fluxes of orthophosphate were greater in the Agricultural watershed that in the Urban watershed while fluxes of TSS were greater in the Mixed watershed when compared to the Agricultural watershed. Fluxes of nitrate, ammonium, and total phosphorus did not vary among watersheds. It is apparent from the data that Agricultural land uses are generally a greater source of nutrients than the Urban land uses while Urban land uses are generally a greater source of suspended sediment.     

HYDROLOGIC MODELING OF TWO GLACIATED WATERSHEDS IN NORTHEAST PENNSYLVANIA1

ABSTRACT: A hydrologic modeling study, using the Hydrologic Simulation Program - FORTRAN (HSPF), was conducted in two glaciated watersheds, Purdy Creek and Ariel Creek in northeastern Pennsylvania. Both watersheds have wetlands and poorly drained soils due to low hydraulic conductivity and presence of fragipans. The HSPF model was calibrated in the Purdy Creek watershed and verified in the Ariel Creek watershed for June 1992 to December 1993 period. In Purdy Creek, the total volume of observed stream-flow during the entire simulation period was 13.36 × 10⁶ m³ and the simulated streamflow volume was 13.82 × 10⁶ m³ (5 percent difference). For the verification simulation in Ariel Creek, the difference between the total observed and simulated flow volumes was 17 percent. Simulated peak flow discharges were within two hours of the observed for 30 of 46 peak flow events (discharge greater than 0.1 m³/sec) in Purdy Creek and 27 of 53 events in Ariel Creek. For 22 of the 46 events in Purdy Creek and 24 of 53 in Ariel Creek, the differences between the observed and simulated peak discharge rates were less than 30 percent. These 22 events accounted for 63 percent of total volume of streamflow observed during the selected 46 peak flow events in Purdy Creek. In Ariel Creek, these 24 peak flow events accounted for 62 percent of the total flow observed during all peak flow events. Differences in observed and simulated peak flow rates and volumes (on a percent basis) were greater during the snowmelt runoff events and summer periods than for other times.