A Spatially Explicit Model for Mapping Headwater Streams

Headwater streams are the primary sources of water in a drainage network and serve as a critical hydrologic link between the surrounding landscape and larger, downstream surface waters. Many states, including North Carolina, regulate activity in and near headwater streams for the protection of water quality and aquatic resources. A fundamental tool for regulatory management is an accurate representation of streams on a map. Limited resources preclude field mapping every headwater stream and its origin across a large region. It is more practical to develop a model for headwater streams based on a sample of field data that can then be extrapolated to a larger area of interest. The North Carolina Division of Water Quality has developed a cost‐effective method for modeling and mapping the location, length, and flow classification (intermittent and perennial) of headwater streams. We used a multiple logistic regression approach that combined field data and terrain derivatives for watersheds located in the Triassic Basins ecoregion. Field data were collected using a standard methodology for identifying headwater streams and origins. Terrain derivatives were generated from digital elevation models interpolated from bare‐earth Light Detection and Range data. Model accuracies greater than 80% were achieved in classifying stream presence and absence, stream length and perennial stream length, but were not as consistent in predicting intermittent stream length.     

Bankfull Regional Curves for North and Northwest Florida Streams1

Abstract: Regional curves, which relate bankfull channel dimensions and discharge to watershed drainage area, are developed to aid in identifying the bankfull stage in ungaged watersheds, and estimating the bankfull discharge and dimensions for river studies and natural channel design applications. This study assessed 26 stable stream reaches in two hydro‐physiographic regions of the Florida Coastal Plain: the Northwest Florida Coastal Plain (NWFCP) and the North Florida Coastal Plain (NFCP). Data from stream reaches in Georgia and Alabama were also used to develop the Florida regional curves, since they are located in the same hydro‐physiographic region. Reaches were selected based on the presence of U.S. Geological Survey gage stations and indicators of limited watershed development (e.g., <10% impervious surface). Analyses were conducted to determine bankfull channel dimensions, bankfull discharge, average channel slope, and Rosgen stream classification. Based on these data, significant relationships were found between bankfull cross‐sectional area, width, mean depth, and discharge as a function of drainage area for both regions. Data from this study suggested that bankfull discharges and channel dimensions were larger from NWFCP streams than from Coastal Plain streams in North Carolina and Maryland. Bankfull discharges were similar between NFCP and Georgia coastal plain streams; therefore, the data were combined into one regional curve. In addition, the data were stratified by Rosgen stream type. This stratification strengthened the relationships of bankfull width and mean depth as a function of drainage area.     

HYDRAULIC GEOMETRY RELATIONSHIPS FOR URBAN STREAMS THROUGHOUT THE PIEDMONT OF NORTH CAROLINA1

ABSTRACT: Hydraulic geometry relationships, or regional curves, relate bankfull stream channel dimensions to watershed drainage area. Hydraulic geometry relationships for streams throughout North Carolina vary with hydrology, soils, and extent of development within a watershed. An urban curve that is the focus of this study shows the bankfull features of streams in urban and suburban watersheds throughout the North Carolina Piedmont. Seventeen streams were surveyed in watersheds that had greater than 10 percent impervious cover. The watersheds had been developed long enough for the streams to redevelop bankfull features, and they had no major impoundments. The drainage areas for the streams ranged from 0.4 to 110.3 square kilometers. Cross‐sectional and longitudinal surveys were conducted to determine the channel dimension, pattern, and profile of each stream and power functions were fitted to the data. Comparisons were made with regional curves developed previously for the rural Piedmont, and enlargement ratios were produced. These enlargement ratios indicated a substantial increase in the hydraulic geometry for the urban streams in comparison to the rural streams. A comparison of flood frequency indicates a slight decrease in the bankfull discharge return interval for the gaged urban streams as compared to the gaged rural streams. The study data were collected by North Carolina State University (NCSU), the University of North Carolina at Charlotte (UNC), and Charlotte Storm Water Services. Urban regional curves are useful tools for applying natural channel design in developed watersheds. They do not, however, replace the need for field calibration and verification of bankfull stream channel dimensions.     

Using GIS to Delineate Headwater Stream Origins in the Appalachian Coalfields of Kentucky

Headwater streams have a significant nexus or physical, chemical, and/or biological connection to downstream reaches. Generally, defined as 1st‐3rd order with ephemeral, intermittent, or perennial flow regimes, these streams account for a substantial portion of the total stream network particularly in mountainous terrain. Due to their often remote locations, small size, and large numbers, conducting field inventories of headwater streams is challenging. A means of estimating headwater stream location and extent according to flow regime type using publicly available spatial data is needed to simplify this complex process. Using field‐collected headwater point of origin data from three control watersheds, streams were characterized according to a set of spatial parameters related to topography, geology, and soils. These parameters were (1) compared to field‐collected point of origin data listed in three nearby Jurisdictional Determinations, (2) used to develop a geographic information system (GIS)‐based stream network for identifying ephemeral, intermittent, and perennial streams, and (3) applied to a larger watershed and compared to values obtained using the high‐resolution National Hydrography Dataset (NHD). The parameters drainage area and local valley slope were the most reliable predictors of flow regime type. Results showed the high‐resolution NHD identified no ephemeral streams and 9 and 65% fewer intermittent and perennial streams, respectively, than the GIS model.     

Classification of Ephemeral,Intermittent, and Perennial Stream Reaches Using a TOPMODEL‐Based Approach

Whether a waterway is temporary or permanent influences regulatory protection guidelines, however, classification can be subjective due to a combination of factors, including time of year, antecedent moisture conditions, and previous experience of the field investigator. Our objective was to develop a standardized protocol using publically available spatial information to classify ephemeral, intermittent, and perennial streams. Our hypothesis was that field observations of flow along the stream channel could be compared to results from a hydrologic model, providing an objective method of how these stream reaches can be identified. Flow‐state sensors were placed at ephemeral, intermittent, and perennial stream reaches from May to December 2011 in the Appalachian coal basin of eastern Kentucky. This observed flow record was then used to calibrate the simulated saturation deficit in each channel reach based on the topographic wetness index used by TOPMODEL. Saturation deficit values were categorized as flow or no‐flow days, and the simulated record of streamflow was compared to the observed record. The hydrologic model was more accurate for simulating flow during the spring and fall seasons. However, the model effectively identified stream reaches as intermittent and perennial in each of the two basins.     

Channel Roughness in North Carolina Mountain Streams

Channel roughness, often described by Manning's n, is used to represent the amount of resistance that flow encounters, and has direct implications on velocity and discharge. Ideally, n is calculated from a long‐term record of channel discharge and hydraulic geometry. In the absence of these data, a combination of photo references and a validated qualitative method is preferable to simply choosing n arbitrarily or from a table. The purpose of this study was to use United States Geological Survey (USGS) streamflow data to calculate roughness coefficients for streams in the mountains of North Carolina. Five USGS gage stations were selected for this study, representing drainage areas between 71.5 and 337 km². Photo references of the study sites are presented. Measured discharges were combined with hydraulic geometry at a cross‐section to calculate roughness coefficients for flows of interest. At bankfull flow, n ranged between 0.039 and 0.064 for the five study sites. Roughness coefficients were not constant for all flows in a channel, and fluctuated over a large range. At all sites, roughness was highest during low‐flow conditions, then quickly decreased as flow increased, up to the bankfull elevation.     

AVERAGE DISCHARGE,PERENNIAL FLOW INITIATION,AND CHANNEL INITIATION ‐ SMALL SOUTHERN APPALACHIAN BASINS1

ABSTRACT: Regional average evapotranspiration estimates developed by water balance techniques are frequently used to estimate average discharge in ungaged streams. However, the lower stream size range for the validity of these techniques has not been explored. Flow records were collected and evaluated for 16 small streams in the Southern Appalachians to test whether the relationship between average discharge and drainage area in streams draining less than 200 acres was consistent with that of larger basins in the size range (> 10 square miles) typically gaged by the U.S. Geological Survey (USGS). This study was designed to evaluate predictors of average discharge in small ungaged streams for regulatory purposes, since many stream regulations, as well as recommendations for best management practices, are based on measures of stream size, including average discharge. The average discharge/drainage area relationship determined from gages on large streams held true down to the perennial flow initiation point. For the southern Appalachians, basin size corresponding to perennial flow is approximately 19 acres, ranging from 11 to 32 acres. There was a strong linear relationship (R²= 0.85) between average discharge and drainage area for all streams draining between 16 and 200 acres, and the average discharge for these streams was consistent with that predicted by the USGS Unit Area Runoff Map for Georgia. Drainage area was deemed an accurate predictor of average discharge, even in very small streams. Channel morphological features, such as active channel width, cross‐sectional area, and bankfull flow predicted from Manning's equation, were not accurate predictors of average discharge. Monthly baseflow statistics also were poor predictors of average discharge.     

BANKFULL DISCHARGE RECURRENCE INTERVALS AND REGIONAL HYDRAULIC GEOMETRY RELATIONSHIPS: PATTERNS IN THE PACIFIC NORTHWEST,USA1

ABSTRACT: The model bankfull discharge recurrence interval (annual series) (T_a) in streams has been approximated at a 1.5‐year flow event. This study tests the linkage between regional factors (climate, physiography, and ecoregion) and the frequency of bank‐full discharge events in the Pacific Northwest (PNW). Patterns of Ta were found to be significant when stratified by EPA Ecoregion. The mean value for T_a in the PNW is 1.4 years; however, when the data is stratified by ecoregion, the humid areas of western Oregon and Washington have a mean value of 1.2 years, while the dryer areas of Idaho and eastern Oregon and Washington have a mean value of 1.4 to 1.5 years. Among the four factors evaluated, vegetation association and average annual precipitation are the primary factors related to channel form and T_a. Based on the results of the T_a analyses, regional hydraulic geometry relationships of streams were developed for the PNW, which relate variables, such as bank‐full cross‐sectional area, width, depth, and velocity, to bankfull discharge and drainage area. The verification of T_a values, combined with the development of regional hydraulic geometry relationships, provides geographically relevant information that will result in more accurate estimates of hydraulic geometry variables in the PNW.     

Southwestern Intermittent and Ephemeral Stream Connectivity

Ephemeral and intermittent streams are abundant in the arid and semiarid landscapes of the Western and Southwestern United States (U.S.). Connectivity of ephemeral and intermittent streams to the relatively few perennial reaches through runoff is a major driver of the ecohydrology of the region. These streams supply water, sediment, nutrients, and biota to downstream reaches and rivers. In addition, they provide runoff to recharge alluvial and regional groundwater aquifers that support baseflow in perennial mainstem stream reaches over extended periods when little or no precipitation occurs. Episodic runoff, as well as groundwater inflow to surface water in streams support limited naturally occurring riparian communities. This paper provides an overview and comprehensive examination of factors affecting the hydrologic, chemical, and ecological connectivity of ephemeral and intermittent streams on perennial or intermittent rivers in the arid and semiarid Southwestern U.S. Connectivity as influenced and moderated through the physical landscape, climate, and human impacts to downstream waters or rivers is presented first at the broader Southwestern scale, and secondly drawing on a specific and more detailed example of the San Pedro Basin due to its history of extensive observations and research in the basin. A wide array of evidence clearly illustrates hydrologic, chemical, and ecological connectivity of ephemeral and intermittent streams throughout stream networks.     

Estimating Contributions of Nitrate and Herbicides From Groundwater to Headwater Streams,Northern Atlantic Coastal Plain,United States1

Abstract: Groundwater transport often complicates understanding of surface‐water contamination. We estimated the regional flux of nitrate and selected herbicides from groundwater to nontidal headwater streams of the Atlantic Coastal Plain (New Jersey through North Carolina) based on late‐winter or spring base‐flow samples from 174 streams. Sampled streams were selected randomly, and flux estimates are based on resulting population estimates rather than on empirical models, which have been used previously for similar estimates. Base‐flow flux in the estimated 8,834 headwater streams of the study area are an estimated 21,200 kg/day of nitrate (as N) and 5.83, 0.565, and 20.7 kg/day of alachlor, atrazine, and metolachlor (and selected degradates), respectively. Base‐flow flux of alachlor and metolachlor is <3% of the total base‐flow flux of those compounds plus degradates. Base‐flow flux of nitrate and herbicides as a percentage of applications is typically highest in well‐drained areas and lowest in areas with abundant poor drainage and anoxic conditions. In Coastal Plain watersheds of Albemarle and Pamlico Sounds, <2% of applied nitrogen reaches headwater streams as base flow. On the Delmarva Peninsula part of the Chesapeake Bay watershed, however, more than 10% of such applications are transported through groundwater to streams, and base‐flow nitrate flux represents 70% of total nitrogen flux in headwater streams.     

Morphology Characteristics of Southern Appalachian Wilderness Streams1

Zink, Jason M., Gregory D. Jennings, and G. Alexander Price, 2012. Morphology Characteristics of Southern Appalachian Wilderness Streams. Journal of the American Water Resources Association (JAWRA) 48(4): 762‐773. DOI: 10.1111/j.1752‐1688.2012.00647.x Abstract: Watersheds without urbanization or impacts from logging are rare in the southern Appalachian Mountains. The Joyce Kilmer/Slickrock Wilderness of North Carolina and Tennessee contains 24 km² of old‐growth forest, with the balance of the wilderness in a mature second‐growth forest. The watersheds of Little Santeetlah and Slickrock Creek are located within the wilderness. Morphological information, including channel dimensions and longitudinal profiles, was gathered from 14 alluvial stream reaches in these watersheds. The study sites had drainage areas from 0.25 to 41.6 km² and stream slopes from 0.014 to 0.104 m/m. Bankfull cross‐section dimensions of the study stream reaches were strongly correlated to drainage area across the observed range of slopes and bed morphology. Cross‐section area and width relationships for the streams in this study did not differ significantly from regional curves for the mountain physiographic region of North Carolina. Observations of these reaches did not suggest a definitive rule regarding the proportion of steps and riffles in streams. Pools occupied greater than 50% of the length in all stream reaches with slopes less than 0.07 m/m. Significant correlation existed between step height ratio and slope, suggesting that step height can be approximated as the product of channel width and slope. Riffle length and riffle slope ratios were also significantly correlated with slope, though pool spacing was not.     

Wet/Dry Mapping: Using Citizen Scientists to Monitor the Extent of Perennial Surface Flow in Dryland Regions

Wet/dry mapping provides a low-cost, comprehensive snapshot for monitoring flow conditions in rivers with interrupted perennial (spatially intermittent) surface flow. When used in conjunction with more traditional point-specific stream flow or groundwater measurements, it provides a better understanding of hydrologic systems at the broad landscape or watershed scale. Through use of trained volunteers, we mapped reaches with surface water during the driest time of year to track annual variation in length and location of perennial flow. Data from 12 years of wet/dry mapping on the San Pedro River in Arizona, USA, showed 62 reaches with surface flow in every year, totaling 32% of the river length through the San Pedro Riparian National Conservation Area. They also show areas with high year-to-year variation in flow length, which indicate changes in local groundwater conditions and may provide early warning of ecological changes. Data and maps from this project have been useful for a wide variety of conservation, management, and research efforts.     

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.     

Morphology of a Wetland Stream

/ Little attention has been paid to wetland stream morphology in the geomorphological and environmental literature, and in the recently expanding wetland reconstruction field, stream design has been based primarily on stream morphologies typical of nonwetland alluvial environments. Field investigation of a wetland reach of Roaring Brook, Stafford, Connecticut, USA, revealed several significant differences between the morphology of this stream and the typical morphology of nonwetland alluvial streams. Six morphological features of the study reach were examined: bankfull flow, meanders, pools and riffles, thalweg location, straight reaches, and cross-sectional shape. It was found that bankfull flow definitions originating from streams in nonwetland environments did not apply. Unusual features observed in the wetland reach include tight bends and a large axial wavelength to width ratio. A lengthy straight reach exists that exceeds what is typically found in nonwetland alluvial streams. The lack of convex bank point bars in the bends, a greater channel width at riffle locations, an unusual thalweg location, and small form ratios (a deep and narrow channel) were also differences identified. Further study is needed on wetland streams of various regions to determine if differences in morphology between alluvial and wetland environments can be applied in order to improve future designs of wetland channels.KEY WORDS: Stream morphology; Wetland restoration; Wetland creation; Bankfull; Pools and riffles; Meanders; Thalweg     

Biophysical‐Regulatory Classification and Profiling of Streams Across Management Units and Ecoregions1

Caruso, Brian S. and Joshua Haynes, 2011. Biophysical‐Regulatory Classification and Profiling of Streams Across Management Units and Ecoregions. Journal of the American Water Resources Association (JAWRA) 00(0):1‐22. DOI: 10.1111/j.1752‐1688.2010.00522.x Abstract: Aquatic resources management in the United States (U.S.) under Clean Water Act Section 404 has become more complex after recent Supreme Court decisions and U.S. Army Corps of Engineers and Environmental Protection Agency (USEPA) guidance. Many intermittent/ephemeral and headwater streams may not be jurisdictional if they lack a significant nexus with navigable waters. Streams in semiarid USEPA Region 8 were classified based on hydrologic permanence and stream order using National Hydrography Dataset (NHD) Plus and GIS to provide information across broad spatial scales to aid with jurisdictional determinations (JDs). Four classes were developed for profiling across management units and ecoregions. Based on medium‐resolution NHDPlus data, intermittent streams comprise >¾, and first order streams constitute >½ of the total stream length in Region 8. Mountain states and ecoregions have the largest percentage of perennial first order streams, whereas the Dakotas, plains, and desert ecoregions have the greatest percentages of intermittent first order and intermittent higher order streams. In the Upper Colorado River Basin, >50% of reaches are intermittent first order, and 9% are perennial first order. NHDPlus data can significantly underestimate the length of headwater and intermittent streams, but can still be a valuable tool to help develop stream classes and for regional JD planning and analysis. Refinement of the stream classes using high resolution NHD data and other key catchment parameters can improve their utility for JDs.     

Effect of Soils on Water Quantity and Quality in Piedmont Forested Headwater Watersheds of North Carolina1

Boggs, Johnny, Ge Sun, David Jones, and Steven G. McNulty, 2012. Effect of Soils on Water Quantity and Quality in Piedmont Forested Headwater Watersheds of North Carolina. Journal of the American Water Resources Association (JAWRA) 1‐19. DOI: 10.1111/jawr.12001 Abstract: Water quantity and quality data were compared from six headwater watersheds on two distinct soil formations, Carolina Slate Belt (CSB) and Triassic Basins (TB). CSB soils are generally thicker, less erodible, and contain less clay content than soils found in TB. TB generated significantly more discharge/precipitation ratio than CSB (0.33 vs. 0.24) in the 2009 dormant season. In the 2009 growing season, TB generated significantly less discharge/precipitation ratio than CSB (0.02 vs. 0.07). Over the entire monitoring period, differences in discharge/precipitation ratios between CSB and TB were not significantly different (0.17 vs. 0.20, respectively). Storm‐flow rates were significantly higher in TB than CSB in both dormant and growing season. Benthic macroinvertebrate biotic index scores were excellent for all streams. Nutrient concentrations and exports in CSB and TB were within background levels for forests. Low‐stream nitrate and ammonium concentrations and exports suggested that both CSB and TB were nitrogen limited. Soils appear to have had a significant influence on seasonal and storm‐flow generation, but not on long‐term total water yield and water quality under forested conditions. This study indicated that watersheds on TB soils might be more prone to storm‐flow generation than on CSB soils when converted from forest to urban. Future urban growth in the area should consider differences in baseline hydrology and effects of landuse change on water quantity and quality.     

Modeling Stream Channel Characteristics From Drainage‐Enforced DEMs in Puget Sound,Washington, USA1

Abstract: Mapping stream channels and their geomorphic attributes is an important step in many watershed research and management projects. Often insufficient field data exist to map hydromorphologic attributes across entire drainage basins, necessitating the application of hydrologic modeling tools to digital elevation models (DEMs) via a geographic information system (GIS). In this article, we demonstrate methods for deriving synthetic stream networks via GIS across large and diverse basins using drainage‐enforced DEMs, along with techniques for estimating channel widths and gradient on the reach scale. The two‐step drainage enforcement method we used produced synthetic stream networks that displayed a high degree of positional accuracy relative to the input streams. The accuracies of our estimated channel parameters were assessed with field data, and predictions of bankfull width, wetted width and gradient were strongly correlated with measured values (r² = 0.92, r² = 0.95, r² = 0.88, respectively). Classification accuracies of binned channel attributes were also high. Our methodology allows for the relatively rapid mapping of stream channels and associated morphological attributes across large geographic areas. Although initially developed to provide salmon recovery planners with important salmon habitat information, we suggest these methodologies are relevant to a variety of research and management questions.     

Bankfull Regional Curves for the Alleghany Plateau/Valley and Ridge,Piedmont, and Coastal Plain Regions of Maryland

Regional curves are empirical relationships that can help identify the bankfull stage in ungaged watersheds and aid in designing the riffle dimension in stream restoration projects. Bankfull regional curves were developed from gage stations with drainage areas less than 102 mi² (264.2 km²) for the Alleghany Plateau/Valley and Ridge (AP/VR), Piedmont, and Coastal Plain regions of Maryland. The AP/VR regions were combined into one region for this project. These curves relate bankfull discharge, cross‐sectional area, width, and mean depth to drainage area within the same hydro‐physiographic region (region with similar rainfall/runoff relationship). The bankfull discharge curve for the Coastal Plain region was further subdivided into the Western Coastal Plain (WCP) and Eastern Coastal Plain (ECP) region due to differences in topography and runoff. Results show that the Maryland Piedmont yields the highest bankfull discharge rate per unit drainage area, followed by the AP/VR, WCP, and ECP. Likewise, the Coastal Plain and AP/VR streams have less bankfull cross‐sectional area per unit drainage area than the Piedmont. The average bankfull discharge return interval across the three hydro‐physiographic regions was 1.4 years. The Maryland regional curves were compared to other curves in the eastern United States. The average bankfull discharge return interval for the other studies ranged from 1.1 to 1.8 years.     

Evaluating the Eco‐Geomorphological Condition of Restored Streams Using Visual Assessment and Macroinvertebrate Metrics

The Stream Performance Assessment (SPA), a new rapid assessment method, was applied to 93 restored, 21 impaired, 29 reference, and 13 reference streams with some incision throughout North Carolina. Principal component analysis (PCA) indicated restored streams align more closely with reference streams rather than impaired streams. Further, PCA‐based factor analysis revealed restored streams were similar to reference streams in terms of morphologic condition, but exhibited a greater range of scores relative to aquatic habitat and bedform. Macroinvertebrate sampling and GIS watershed analyses were conducted on 84 restored streams. SPA and watershed data were compared to Ephemeroptera, Plecoptera, and Trichoptera (EPT) taxa to determine which factors indicate stream health. SPA and watershed factors were used in least squares, ridge, and principal component regression (PCR) to develop a prediction model for EPT taxa. All three methods produced reasonable predictions for EPT taxa. Cross‐validation indicated ridge regression resulted in the lowest prediction error. The ridge model was then used to predict EPT taxa numbers for 21 impaired and 25 reference streams in addition to the 84 restored streams. Statistical comparisons of the predicted scores indicated urban streams (>10% impervious watershed cover) have lower expected numbers of EPT taxa. Rural restored streams have macroinvertebrate metric scores similar to those predicted for rural reference streams.