MOBILITY AND AQUATIC TOXICITY OF COPPER IN AN URBAN WATERSHED1

ABSTRACT: Abundant use of copper based products has resulted in increased violation of copper water quality criteria in runoff from urban storm water systems. The objectives of this work were to understand the mobility and toxicity of copper in an urban watershed and to apportion the amount of copper entering the freshwater receiving stream from different urban land covers using a mass balance approach. Sixteen rainfall events collected from the University of Connecticut study watershed between August 1998 and September 2000 were analyzed to assess copper flux in an urban storm water system. Mean flow weighted dissolved copper concentrations observed in the study for copper based architectural material runoff, pervious area runoff, impervious area runoff, and in the receiving stream were 1210 ± 840, 9 ± 3, 8 ± 2, and 14 ± 7 μg/L, respectively. Mean dissolved copper concentrations in the receiving stream exceeded Connecticut's water quality criteria. Despite exceeding the dissolved concentration based criteria, cupric ion concentrations at the system outlet remained below 0.05 μg/L for all storms analyzed, and no acute toxicity (using Daphnia pulex as the test organism) was measured in samples collected from the stream.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

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.     

STORMFLOW SIMULATION USING A GEOGRAPHICAL INFORMATION SYSTEM WITH A DISTRIBUTED APPROACH1

ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.     

EFFECTS OF MULTI‐SCALE ENVIRONMENTAL CHARACTERISTICS ON AGRICULTURAL STREAM BIOTA IN EASTERN WISCONSIN1

ABSTRACT: The U.S. Geological Survey examined 25 agricultural streams in eastern Wisconsin the determine relations between fish, invertebrate, and algal metrics and multiple spatial scales of land cover, geologic setting, hydrologic, aquatic habitat, and water chemistry data. Spearman correlation and redundancy analyses were used to examine relations among biotic metrics and environmental characteristics. Riparian vegetation, geologic, and hydrologic conditions affected the response of biotic metrics to watershed agricultural land cover but the relations were aquatic assemblage dependent. It was difficult to separate the interrelated effects of geologic setting, watershed and buffer land cover, and base flow. Watershed and buffer land cover, geologic setting, reach riparian vegetation width, and stream size affected the fish IBI, invertebrate diversity, diatom IBI, and number of algal taxa; however, the invertebrate FBI, percentage of EPT, and the diatom pollution index were more influenced by nutrient concentrations and flow variability. Fish IBI scores seemed most sensitive to land cover in the entire stream network buffer, more so than watershed‐scale land cover and segment or reach riparian vegetation width. All but one stream with more than approximately 10 percent buffer agriculture had fish IBI scores of fair or poor. In general, the invertebrate and algal metrics used in this study were not as sensitive to land cover effects as fish metrics. Some of the reach‐scale characteristics, such as width/depth ratios, velocity, and bank stability, could be related to watershed influences of both land cover and geologic setting. The Wisconsin habitat index was related to watershed geologic setting, watershed and buffer land cover, riparian vegetation width, and base flow, and appeared to be a good indicator of stream quality Results from this study emphasize the value of using more than one or two biotic metrics to assess water quality and the importance of environmental characteristics at multiple scales.     

Boosted Regression Tree Models to Explain Watershed Nutrient Concentrations and Biological Condition

Boosted regression tree (BRT) models were developed to quantify the nonlinear relationships between landscape variables and nutrient concentrations in a mesoscale mixed land cover watershed during base‐flow conditions. Factors that affect instream biological components, based on the Index of Biotic Integrity (IBI), were also analyzed. Seasonal BRT models at two spatial scales (watershed and riparian buffered area [RBA]) for nitrite‐nitrate (NO₂‐NO₃), total Kjeldahl nitrogen, and total phosphorus (TP) and annual models for the IBI score were developed. Two primary factors — location within the watershed (i.e., geographic position, stream order, and distance to a downstream confluence) and percentage of urban land cover (both scales) — emerged as important predictor variables. Latitude and longitude interacted with other factors to explain the variability in summer NO₂‐NO₃ concentrations and IBI scores. BRT results also suggested that location might be associated with indicators of sources (e.g., land cover), runoff potential (e.g., soil and topographic factors), and processes not easily represented by spatial data indicators. Runoff indicators (e.g., Hydrological Soil Group D and Topographic Wetness Indices) explained a substantial portion of the variability in nutrient concentrations as did point sources for TP in the summer months. The results from our BRT approach can help prioritize areas for nutrient management in mixed‐use and heavily impacted watersheds.     

STREAM DISCHARGE PREDICTION VIA A GRID BASED SOIL WATER ROUTING WITH PADDY FIELDS1

ABSTRACT: A grid based daily hydrologic model for a watershed with paddy fields was developed to predict the stream discharge. ASCII formatted elevation, soil, and land use data supported by the GRASS Geographic Information System are used to generate distributed results such as surface runoff and subsurface flow, soil water content, and evapotranspiration. The model uses a single flow path algorithm and simulates a water balance at each grid element. A linear reservoir assumption was used to predict subsurface runoff components. The model was applied to a 75.6 km² watershed located in the middle of South Korea, and observed stream flow hydrographs from 1995 and 1996 were compared to model predictions. The stream flow predictions of 1995 and 1996 generally agreed with the observed flow, resulting in a Nash‐Sutcliffe efficiency R² of 0.60 and 0.62, respectively. The hydraulic conductivity for percolating water through the saturated layer affected baseflow generation. The levee height of the paddy influenced the time and magnitude of the surface runoff, depending on irrigation management. The model will be used for making low flow management decisions by evaluating the role of each land use to stream flow, especially in case of paddy decrease by gradual urbanization of a watershed.     

Estimating Potential E. coli Sources in a Watershed Using Spatially Explicit Modeling Techniques1

Riebschleager, K.J., R. Karthikeyan, R. Srinivasan, and K. McKee, 2012. Estimating Potential E. coli Sources in a Watershed Using Spatially Explicit Modeling Techniques. Journal of the American Water Resources Association (JAWRA) 48(4): 745‐761. DOI: 10.1111/j.1752‐1688.2012.00649.x Abstract: The Spatially Explicit Load Enrichment Calculation Tool (SELECT) was automated to characterize waste and the associated pathogens from various sources within a mixed land use watershed. Potential Escherichia coli loads in Lake Granbury watershed were estimated using spatially variable governing factors, such as land use, soil condition, and distance to streams. A new approach for characterizing E. coli loads resulting from malfunctioning on‐site wastewater treatment systems (OWTSs) was incorporated into SELECT along with the Pollutant Connectivity Factor (PCF) module. The PCF component was applied to identify areas contributing E. coli loads during runoff events by incorporating the influence of potential E. coli loading, runoff potential, and travel distance to waterbodies. Simulation results indicated livestock and wildlife are potential E. coli contributing sources in the watershed. The areas in which these sources are potentially contributing are not currently monitored for E. coli. The bacterial water quality violations seen around Lake Granbury are most likely the result of malfunctioning OWTSs and pet wastes. SELECT results demonstrate the need to evaluate each contributing source separately to effectively allocate site specific best management practices (BMPs) utilizing stakeholder inputs. It also serves as a powerful screening tool for determining areas where detailed investigation is merited.     

GIS‐Based Stream Classification in a Mountain Watershed for Jurisdictional Evaluation

This study analyzed stream characteristics in a mountain watershed in southwestern Colorado and developed a three‐level hierarchical classification scheme using national datasets to demonstrate jurisdictional evaluation as “waters of the United States (U.S.)” under U.S. Clean Water Act Section 404 at the watershed scale. The National Hydrography Dataset and USGS StreamStats were used with field observations to classify streams in the 53 km² Cement Creek Watershed based on flow duration (Level 1), stream order (Level 2), and other biophysical metrics (Level 3). Kruskal‐Wallis tests and discriminant analysis showed significant differences among Level 2 classes. Level 3 classification used cluster analysis for stream length, distance to the downstream traditional navigable water (TNW), and the ratio of mean annual flow from the source stream to the TNW. Results showed all perennial and intermittent streams are jurisdictional relatively permanent waters (RPWs), which include over a third of all streams, 64% are intermittent or ephemeral, and almost half are ephemeral first order. All ephemeral reaches are non‐RPWs requiring significant nexus evaluation to determine jurisdiction. These ephemeral first‐order streams can contribute 5% of the annual flow to the TNW at the confluence, while the Cement Creek main stem contributes 21% of the TNW flow. The study demonstrated that the classification provides key biophysical and regulatory information to aid jurisdictional evaluations in mountain watersheds.     

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.     

Incorporating Variable Source Area Hydrology into a Spatially Distributed Direct Runoff Model1

Buchanan, Brian, Zachary M. Easton, Rebecca Schneider, and M. Todd Walter, 2011. Incorporating Variable Source Area Hydrology Into a Spatially Distributed Direct Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(1): 43‐60. DOI: 10.1111/j.1752‐1688.2011.00594.x Abstract: Few hydrologic models simulate both variable source area (VSA) hydrology, and runoff‐routing at high enough spatial resolutions to capture fine‐scale hydrologic pathways connecting VSA to the stream network. This paper describes a geographic information system‐based operational model that simulates the spatio‐temporal dynamics of VSA runoff generation and distributed runoff‐routing, including through complex artificial drainage networks. The model combines the Natural Resource Conservation Service’s Curve Number (CN) equation for estimating storm runoff with the topographic index concept for predicting the locations of VSA and a runoff‐routing algorithm into a new spatially distributed direct hydrograph (SDDH) model (SDDH‐VSA). Using a small agricultural watershed in central New York, SDDH‐VSA results were compared to those from a SDDH model using the traditional land use assumptions for the CN (SDDH‐CN). The SDDH‐VSA model generally agreed better with observed discharge than the SDDH‐CN model (average, Nash‐Sutcliffe efficiency of 0.69 vs. 0.58, respectively) and resulted in more realistic spatial patterns of runoff‐generating areas. The SDDH approach did not correctly capture the timing of runoff from small storms in dry periods. Despite this type of limitation, SDDH‐VSA extends the applicability of the SDDH technique to VSA conditions, providing a basis for new tools to help identify critical management areas and assess water quality risks due to landscape alterations.     

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.     

SPATIALLY EXPLICIT HYDROLOGIC MODELING OF LAND USE CHANGE1

ABSTRACT: Using a Geographic Information System (GIS), a method is presented to develop a spatially explicit time series of land use in an urbanizing watershed. The method is prefaced on the existence of independent observations of land use at different times and data that describes the spatial‐temporal land use transition characteristics of the watershed between these two points in time. A method is then presented to generalize the TR‐55 graphical method, a common lumped hydrologic model for estimating peak discharge, for use in a spatially explicit scheme. This scheme predicts peak discharge throughout a watershed, rather than at a single selected watershed outlet. Coupling these two methods allows the engineer to model both the temporal and spatial evolution of peak discharge for the watershed. An illustrative watershed in a suburban area of Washington, DC is selected to demonstrate the methods. The model results from these analyses are presented graphically to highlight the complex features in peak discharge behavior that exist both spatially, as a function of position within the watershed drainage network, and temporally, as the watershed undergoes urbanization. These features are not commonly noted in most hydrologic analyses but are captured in these analyses because of the high spatial and temporal resolution of the methods presented. The physical implications of the modeled results are discussed in the context of the information content of a stream gauge located at the overall outlet of the illustrative watershed. This work shows that the common practice of transposition of gauge information to locations internal to the watershed would neglect internal variability in peak discharge behavior, and could potentially lead to the determination of inappropriate design discharges.     

An Optimization Method for Estimating Constituent Mean Concentrations in Base Flow‐Dominated Flow

Pollutant coefficients have been widely used to assess runoff nonpoint source pollution from individual land uses (e.g., agricultural, residential) of a watershed. Pollutant coefficients, known as event mean concentrations (EMCs), were developed by the U.S. Environmental Protection Agency's Nationwide Urban Runoff Program (NURP) to serve as a national measure for characterizing pollutant loading in a receiving water body. The term “baseflow pollutant coefficient (BPC)” is used in this study as a surrogate for EMC to describe mean concentration of pollutants in base flow‐dominated flow. A method for characterizing base flow quantity and quality for different land uses was explored using inverse modeling with two optimization techniques (a least square method and a genetic algorithm [GA] optimization), land use information, and streamflow quantity and quality data. The inverse model was formulated as a constrained minimization problem and demonstrated with data for 15 watersheds in Indiana. Results showed that estimated pollutant coefficients are comparable to the published literature. This indicates that the proposed method has the potential to effectively estimate constituent mean concentrations for pollutant load determination in gauged and ungauged watersheds, albeit more analysis with larger and more robust datasets is desirable to further refine and validate the accuracy of the approach.     

Spatial Calibration and Temporal Validation of Flow for Regional Scale Hydrologic Modeling1

Abstract: Physically based regional scale hydrologic modeling is gaining importance for planning and management of water resources. Calibration and validation of such regional scale model is necessary before applying it for scenario assessment. However, in most regional scale hydrologic modeling, flow validation is performed at the river basin outlet without accounting for spatial variations in hydrological parameters within the subunits. In this study, we calibrated the model to capture the spatial variations in runoff at subwatershed level to assure local water balance, and validated the streamflow at key gaging stations along the river to assure temporal variability. Ohio and Arkansas‐White‐Red River Basins of the United States were modeled using Soil and Water Assessment Tool (SWAT) for the period from 1961 to 1990. R² values of average annual runoff at subwatersheds were 0.78 and 0.99 for the Ohio and Arkansas Basins. Observed and simulated annual and monthly streamflow from 1961 to 1990 is used for temporal validation at the gages. R² values estimated were greater than 0.6. In summary, spatially distributed calibration at subwatersheds and temporal validation at the stream gages accounted for the spatial and temporal hydrological patterns reasonably well in the two river basins. This study highlights the importance of spatially distributed calibration and validation in large river basins.     

Evaluation of the MIKE SHE Model for Application in the Loess Plateau,China1

Abstract: Quantifying the hydrologic responses to land use/land cover change and climate variability is essential for integrated sustainable watershed management in water limited regions such as the Loess Plateau in Northwestern China where an adaptive watershed management approach is being implemented. Traditional empirical modeling approach to quantifying the accumulated hydrologic effects of watershed management is limited due to its complex nature of soil and water conservation practices (e.g., biological, structural, and agricultural measures) in the region. Therefore, the objective of this study was to evaluate the ability of the distributed hydrologic model, MIKE SHE to simulate basin runoff. Streamflow data measured from an overland flow‐dominant watershed (12 km²) in northwestern China were used for model evaluation. Model calibration and validation suggested that the model could capture the dominant runoff process of the small watershed. We found that the physically based model required calibration at appropriate scales and estimated model parameters were influenced by both temporal and spatial scales of input data. We concluded that the model was useful for understanding the rainfall‐runoff mechanisms. However, more measured data with higher temporal resolution are needed to further test the model for regional applications.     

Herbicide Transport in Goodwater Creek ExperimentalWatershed: I. Long‐Term Research on Atrazine1

Lerch, R.N., E.J. Sadler, K.A. Sudduth, C. Baffaut, and N.R. Kitchen, 2010. Herbicide Transport in Goodwater Creek Experimental Watershed: I. Long‐Term Research on Atrazine. Journal of the American Water Resources Association (JAWRA) 1‐15. DOI: 10.1111/j.1752‐1688.2010.00503.x Abstract: Atrazine continues to be the herbicide of greatest concern relative to contamination of surface waters in the United States (U.S.). The objectives of this study were to analyze trends in atrazine concentration and load in Goodwater Creek Experimental Watershed (GCEW) from 1992 to 2006, and to conduct a retrospective assessment of the potential aquatic ecosystem impacts caused by atrazine contamination. Located within the Central Claypan Region of northeastern Missouri, GCEW encompasses 72.5 km² of predominantly agricultural land uses, with an average of 21% of the watershed in corn and sorghum. Flow‐weighted runoff and weekly base‐flow grab samples were collected at the outlet to GCEW and analyzed for atrazine. Cumulative frequency diagrams and linear regression analyses generally showed no significant time trends for atrazine concentration or load. Relative annual loads varied from 0.56 to 14% of the applied atrazine, with a median of 5.9%. A cumulative vulnerability index, which takes into account the interactions between herbicide application, surface runoff events, and atrazine dissipation kinetics, explained 63% of the variation in annual atrazine loads. Based on criteria established by the U.S. Environmental Protection Agency, atrazine reached concentrations considered harmful to aquatic ecosystems in 10 of 15 years. Because of its vulnerability, atrazine registrants will be required to work with farmers in GCEW to implement practices that reduce atrazine transport.     

SIMULATING HYDROLOGIC AND WATER QUALITY IMPACTS IN AN URBANIZING WATERSHED1

ABSTRACT: The Hydrologic Simulation Program‐Fortran (HSPF) was calibrated and used to assess the future effects of various land development scenarios on water quality in the Polecat Creek watershed in Caroline County, Virginia. Model parameters related to hydrology and water quality were calibrated and validated using observed stream flow and water quality data collected at the watershed outlet and the outlets of two sub water sheds. Using the county's Comprehensive Plan, land use scenarios were developed by taking into account the trends and spatial distributions of future development. The simulation results for the various land use scenarios indicate that runoff volume and peak rate increased as urban areas increased. Urbanization also increased sediment loads mainly due to increases in channel erosion. Constituent loads of total Kjeldal nitrogen, orthophosphorus, and total phosphorous for Polecat Creek watershed slightly decreased under future development scenarios. These reductions are due to increases in urban areas that typically contribute smaller quantities of nitrogen and phosphorous, as compared to agricultural areas. However, nitrate loads increased for the future land use scenarios, as compared to the existing land use. The increases in nitrate loads may result from increases in residential land and associated fertilizer use and concurrent decreases in forested land. The procedures used in this paper could assist local, state, and regional policy makers in developing land management strategies that minimize environmental impacts while allowing for future development.     

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

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

Linking Hydrologic Alteration to Biological Impairment in Urbanizing Streams of the Puget Lowland,Washington, USA1

Abstract: We used a retrospective approach to identify hydrologic metrics with the greatest potential for ecological relevance for use as resource management tools (i.e., hydrologic indicators) in rapidly urbanizing basins of the Puget Lowland. We proposed four criteria for identifying useful hydrologic indicators: (1) sensitive to urbanization consistent with expected hydrologic response, (2) demonstrate statistically significant trends in urbanizing basins (and not in undeveloped basins), (3) be correlated with measures of biological response to urbanization, and (4) be relatively insensitive to potentially confounding variables like basin area. Data utilized in the analysis included gauged flow and benthic macroinvertebrate data collected at 16 locations in 11 King County stream basins. Fifteen hydrologic metrics were calculated from daily average flow data and the Pacific Northwest Benthic Index of Biological Integrity (B‐IBI) was used to represent the gradient of response of stream macroinvertebrates to urbanization. Urbanization was represented by percent Total Impervious Area (%TIA) and percent urban land cover (%Urban). We found eight hydrologic metrics that were significantly correlated with B‐IBI scores (Low Pulse Count and Duration; High Pulse Count, Duration, and Range; Flow Reversals, T_Qmean, and R‐B Index). Although there appeared to be a great deal of redundancy among these metrics with respect to their response to urbanization, only two of the metrics tested – High Pulse Count and High Pulse Range – best met all four criteria we established for selecting hydrologic indicators. The increase in these high pulse metrics with respect to urbanization is the result of an increase in winter high pulses and the occurrence of high pulse events during summer (increasing the frequency and range of high pulses), when practically none would have occurred prior to development. We performed an initial evaluation of the usefulness of our hydrologic indicators by calculating and comparing hydrologic metrics derived from continuous hydrologic simulations of selected basin management alternatives for Miller Creek, one of the most highly urbanized basins used in our study. We found that the preferred basin management alternative appeared to be effective in restoring some flow metrics close to simulated fully forested conditions (e.g., T_Qmean), but less effective in restoring other metrics such as High Pulse Count and Range. If future research continues to support our hypothesis that the flow regime, particularly High Pulse Count and Range, is an important control of biotic integrity in Puget Lowland streams, it would have significant implications for stormwater management.