COMPARISON OF HSPF OUTPUTS USING FTABLES GENERATED WITH FIELD SURVEY AND DIGITAL DATA1

ABSTRACT: The Hydrological Simulation Program‐FORTRAN (HSPF) describes discharge from a stream reach based on function tables (FTABLES) that relate stream stage, surface area, volume, and discharge. For this study, five FTABLE scenarios were compared to assess their effect on daily discharge rates predicted using HSPF. Four “field‐based” FTABLE scenarios were developed using detailed cross section surveys collected at predefined intervals along 14 reaches in the study watershed. A fifth “digital‐based” scenario was developed using digital elevation models (DEMs) and Natural Resource Conservation Service (NRCS) Regional Hydraulic Geometry Curves. The Smirnov k‐sample test was used to compare average daily discharge rates simulated with HSPF using the five FTABLE scenarios. No significant difference in simulated stream discharge was found (p = 0.99) between the five FTABLE scenarios. Additional examination of the four field‐based scenarios revealed that the number of cross sections per stream reach used to generate FTABLES had little effect on the resulting stage discharge relationship. These findings suggest that FTABLES generated using digital data are a viable option when simulating stream discharge with HSPF and that if field data are used to generate FTABLES, using fewer cross sections will not adversely affect simulated discharge predictions.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.     

ENHANCING FECAL COLIFORM TOTAL MAXIMUM DAILY LOAD MODELS THROUGH BACTERIAL SOURCE TRACKING1

ABSTRACT: Surface water impairment by fecal coliform bacteria is a water quality issue of national scope and importance. In Virginia, more than 400 stream and river segments are on the Commonwealth's 2002 303(d) list because of fecal coliform impairment. Total maximum daily loads (TMDLs) will be developed for most of these listed streams and rivers. Information regarding the major fecal coliform sources that impair surface water quality would enhance the development of effective watershed models and improve TMDLs. Bacterial source tracking (BST) is a recently developed technology for identifying the sources of fecal coliform bacteria and it may be helpful in generating improved TMDLs. Bacterial source tracking was performed, watershed models were developed, and TMDLs were prepared for three streams (Accotink Creek, Christians Creek, and Blacks Run) on Virginia's 303(d) list of impaired waters. Quality assurance of the BST work suggests that these data adequately describe the bacteria sources that are impairing these streams. Initial comparison of simulated bacterial sources with the observed BST data indicated that the fecal coliform sources were represented inaccurately in the initial model simulation. Revised model simulations (based on BST data) appeared to provide a better representation of the sources of fecal coliform bacteria in these three streams. The coupled approach of incorporating BST data into the fecal coliform transport model appears to reduce model uncertainty and should result in an improved TMDL.     

MODELING FECAL COLIFORM CONTAMINATION IN THE RIO GRANDE1

This study examines sources of fecal coliform in Segment 2302 of the Rio Grande, located south of the International Falcon Reservoir in southern Texas. The watershed is unique because the contributing drainage areas lie in Texas and Mexico. Additionally, the watershed is mostly rural, with populated communities known as “colonias.” The colonias lack sewered systems and discharge sanitary water directly to the ground surface, thus posing an increased health hazard from coliform bacteria. Monitoring data confirm that Segment 2302 is not safe for contact recreation due to elevated fecal coliform levels. The goal of the study was to simulate the observed exceedences in Segment 2302 and evaluate potential strategies for their elimination. Fecal coliform contributions from ranching and colonia discharges were modeled using the Hydrologic Simulation Program‐Fortran (HSPF). Model results indicated that the regulatory 30‐day geometric mean fecal coliform concentration of 200 colony forming units (cfu) per 100 milliliters is exceeded approximately three times per year for a total of 30 days. Ongoing initiatives to improve wastewater facilities will reduce this to approximately once per year for 14 days. Best management practices will be necessary to reduce cattle access to streams and eliminate all exceedences. The developed model was limited by the relatively sparse flow and fecal coliform data.     

Development of bacteria and benthic total maximum daily loads: a case study, Linville Creek, Virginia

Two total maximum daily load (TMDL) studies were performed for Linville Creek in Rockingham County, Virginia, to address bacterial and benthic impairments. The TMDL program is an integrated watershed management approach required by the Clean Water Act. This paper describes the procedures used by the Center for TMDL and Watershed Studies at Virginia Tech to develop the Linville Creek TMDLs and discusses the key lessons learned from and the ramifications of the procedures used in these and other similar TMDL studies. The bacterial impairment TMDL was developed using the Hydrological Simulation Program-Fortran (HSPF). Fecal coliform loads were estimated through an intensive source characterization process. The benthic impairment TMDL was developed using the Generalized Watershed Loading Function (GWLF) model and the reference watershed approach. The bacterial TMDL allocation scenario requires a 100% reduction in cattle manure direct-deposits to the stream, a 96% reduction in nonpoint-source loadings to the land surface, and a 95% reduction in wildlife direct-deposits to the stream. Sediment was identified as the primary benthic stressor. The TMDL allocation scenario for the benthic impairment requires an overall reduction of 12.3% of the existing sediment loads. Despite the many drawbacks associated with using watershed-scale models like HSPF and GWLF to develop TMDLs, the detailed watershed and pollutant-source characterization required to use these and similar models creates information that stakeholders need to select appropriate corrective measures to address the cause of the water quality impairment when implementing the TMDL.     

CONSTRUCTING A PROBLEM SOLVING ENVIRONMENT TOOL FOR HYDROLOGIC ASSESSMENT OF LAND USE CHANGE1

ABSTRACT: This paper describes the integration of a comprehensive hydrological model known as the Hydrological Simulation Program Fortran (HSPF) into a problem solving environment (PSE) for watershed management. The original PSE concept was a structure providing web‐based access to a suite of models, including HSPF and other models of in‐stream hydrodynamics, biological impacts and economic effects, for the watershed‐wide assessment of alternative land use scenarios. The present paper describes only the HSPF integration into the PSE program. Example applications to the 148 square kilometer (57 square mile) Back Creek subwatershed in the upper Roanoke River system (1,479 square kilometers or 571 square miles) in southwest Virginia are used to illustrate important concepts and linkages between land development and hydrological change using hypothetical' what if'scenarios. The features of HSPF and its limitations in this context are discussed. The paper as such is a proof‐of‐concept paper and not a completion report. It is intended to describe the PSE tool building process rather than analysis of the many possible simulation outcomes. However, the dominance of raw imperviousness as a contributor to hydrograph response is apparent in all the PSE simulations described in this paper.     

ESCHERICHIA COLI SURVIVAL IN MANTLED KARST SPRINGS AND STREAMS,NORTHWEST ARKANSAS OZARKS,USA1

Recent studies indicate fecal coliform bacterial concentrations, including Escherichia coli (E. coli), characteristically vary by several orders of magnitude, depending on the hydrology of storm recharge and discharge. E. coli concentrations in spring water increase rapidly during the rising limb of a storm hydrograph, peak prior to or coincident with the peak of the storm pulse, and decline rapidly, well before the recession of the storm hydrograph. This suggests E. coli are associated with resuspension of sediment during the onset of turbulent flow, and indicates viable bacteria reside within the spring and stream sediments. E. coli inoculated chambers were placed in spring and stream environments within the mantled karst of northwest Arkansas to assess long term (> 75 days) E. coli viability. During the 75‐day study, a 4‐log die‐off of E. coli was observed for chambers placed in the Illinois River, and a 5‐log die‐off for chambers placed in Copperhead Spring. Extrapolation of the regression line for each environment indicates E. coli concentration would reach 1 most probable number (MPN)/100 g sediment at Copperhead Spring in about 105 days, and about 135 days in the Illinois River, based on a starting inoculation of 2.5 × 107 MPN E. coli/100 g of sediment. These in situ observations indicate it is possible for E. coli to survive in these environments for at least four months with no fresh external inputs.     

Instream Bacteria Influences from Bird Habitation of Bridges

The representativeness of ambient water samples collected from bridge crossings has occasionally been challenged because critics contend birds nesting on bridges elevate fecal indicator bacteria concentrations over samples collected from river reaches not spanned by bridges. This study was designed to evaluate the influence, if any, of bridge‐dwelling bird colonies on instream bacteria concentrations. Three bridges in central Texas were sampled under dry‐weather conditions for instream Escherichia coli. Two bridges were inhabited by migratory cliff swallows and one was devoid of birds. Numerous samples were collected from locations upstream, at the upstream bridgeface, and downstream of each bridge to determine whether significant increases in E. coli occurred in a downstream direction when birds were present. E. coli values increased significantly at bridgeface and downstream locations compared to upstream locations throughout the nesting season. During peak bird activity in May, bacteria geometric mean concentrations at bridgeface and downstream locations jumped from background levels <50 to >190 colony forming units (CFU)/100 mL, well above the state geometric mean criterion of 126 CFU/100 mL for primary contact recreation use. Results confirmed that under dry‐weather conditions bird colonies can have a significant impact on bacteria concentrations in the vicinity of the bridges they inhabit and therefore, to avoid this impact, monitoring should occur upstream of bridges.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

Using Spatially Detailed Water‐Quality Data and Solute‐Transport Modeling to Support Total Maximum Daily Load Development1

Walton‐Day, Katherine, Robert L. Runkel, and Briant A. Kimball, 2012. Using Spatially Detailed Water‐Quality Data and Solute‐Transport Modeling to Support Total Maximum Daily Load Development. Journal of the American Water Resources Association (JAWRA) 48(5): 949‐969. DOI: 10.1111/j.1752‐1688.2012.00662.x Abstract: Spatially detailed mass‐loading studies and solute‐transport modeling using OTIS (One‐dimensional Transport with Inflow and Storage) demonstrate how natural attenuation and loading from distinct and diffuse sources control stream water quality and affect load reductions predicted in total maximum daily loads (TMDLs). Mass‐loading data collected during low‐flow from Cement Creek (a low‐pH, metal‐rich stream because of natural and mining sources, and subject to TMDL requirements) were used to calibrate OTIS and showed spatially variable effects of natural attenuation (instream reactions) and loading from diffuse (groundwater) and distinct sources. OTIS simulations of the possible effects of TMDL‐recommended remediation of mine sites showed less improvement to dissolved zinc load and concentration (14% decrease) than did the TMDL (53‐63% decrease). The TMDL (1) assumed conservative transport, (2) accounted for loads removed by remediation by subtracting them from total load at the stream mouth, and (3) did not include diffuse‐source loads. In OTIS, loads were reduced near their source; the resulting concentration was decreased by natural attenuation and increased by diffuse‐source loads during downstream transport. Thus, by not including natural attenuation and loading from diffuse sources, the TMDL overestimated remediation effects at low flow. Use of the techniques presented herein could improve TMDLs by incorporating these processes during TMDL development.     

Simulation of Surface Water for Un‐Gauged Areas with Storage‐Attenuation Wetlands1

Abstract: A nontraditional application of the Hydrological Simulation Program – FORTRAN (HSPF) model to simulate freshwater discharge to upper Charlotte Harbor along Florida’s west coast was performed. This application was different from traditional HSPF applications in three ways. First, the domain of the model was defined based on the hydraulic characteristics of the landforms using small distributed parameter discretization. Second, broad wetland land forms, representing more than 20% of this area, were simulated as reaches with storage‐attenuation characteristics and not as pervious land segments (PERLNDs). Finally, the reach flow‐tables (F‐Tables) were configured in a unique way to be calibrated representing the uncertainty of the storage‐attenuation process. Characterizing wetlands as hydrography elements allows flow from the wetlands to be treated as a stage‐dependent flux. The study was conducted for the un‐gauged portion of the Peace and Myakka rivers in west‐central Florida. Due to low gradient tidal influences, a large portion of the basin is un‐gauged. The objective of this study was to simulate stream flow discharges and to estimate freshwater inflow from these un‐gauged areas to upper Charlotte Harbor. Two local gauging stations were located within the model domain and were used for calibration. Another gauge with a shorter period of record was used for verification. A set of global hydrologic parameters were selected and tested using the parameter optimization software (PEST) during the calibration. Model results were evaluated using PEST and well‐known statistical indices. The correlation coefficients were very high (0.899 and 0.825) for the two calibration stations. Further testing of this approach appears warranted for watersheds with significant wetlands coverage.     

DINOSAUR NATIONAL MONUMENT: THE EVOLUTION OF A FEDERAL RESERVED WATER RIGHT1

ABSTRACT: Dinosaur National Monument, in northwestern Colorado, has become a test case in the establishment of a federal reserved water right to instream flows. For the first time, the Interior Department was forced to rigorously defend its claims in a watershed where the federal government did not control the upstream reaches. Inadequate quantification of minimum flow requirements, court orders, and an apparent Congressional ban on the spending of Water Resources Program funds by the Park Service to quantify its water rights have already placed the Service in a difficult position to protect instream flows for maintaining the ecological integrity of the Monument. As late as 1983, administrators of the Park Service were divided over their legal strategy, many wanting to pursue a policy of claiming “natural, historic” flows rather than “minimum” flows. The conditional right to instream flows panted to the Park Service in 1978 was subject to quantification within five years. That deadline has been extended, but it is not likely that the case will reach final settlement this decade. Until the design and conduct of federal water rights quantifications better integrate public policy and law with science, the principle lesson from Dinosaur may have to be repeated.     

TMDL Balance: A Model for Coastal Water Pollutant Loadings

Bacterial contamination accounts for more than 60% of the impairments included on the 2008 Texas 303(d) List. Many of these bacterial impairments are along the Texas Gulf Coast because coastal waters often are regulated for oyster harvesting, which have strict water quality standards. Under the Clean Water Act, each one of these impaired waterbodies requires a total maximum daily load (TMDL) study to be performed. A recent, statewide study recommended the development and application of simple modeling approaches to address the majority of Texas's bacteria TMDLs, including “… simple load duration curve, GIS [geographic information systems], and/or mass balance models.” We developed the TMDL Balance model in response to this recommendation. TMDL Balance is a steady state, mass balance, GIS‐based model for simulating pollutant loads and concentrations in coastal systems. The model uses plug‐flow reactor and continuously‐stirred tank reactor equations to route spatially distributed point and nonpoint source loads through a watershed via overland flow, non‐tidal flow, and tidal flow, decaying the loads via first‐order kinetics. In this paper, we explain the development of the watershed loading portion of the TMDL Balance model, demonstrating the methodology through a case study: computing bacterial loads in the Copano Bay watershed of southeast Texas. The application highlights an example of distributing bacterial sources spatially based on land use data.     

ADDRESSING TOTAL PHOSPHORUS IMPAIRMENTS WITH WATER QUALITY TRADING1

ABSTRACT: Water quality trading is a voluntary economic process that provides an opportunity for dischargers to reduce the costs associated with meeting a discharge limitation. Trading can provide a cost effective solution for point sources (i.e., wastewater treatment plants) to meet strict effluent limitations set in response to total maximum daily loads (TMDLs). A successful trading program often depends on first determining the trading suitability of a pollutant for a particular watershed. A simple technical approach has been developed to identify sub‐watersheds within the Raritan River Basin, New Jersey, where water quality trading could provide a cost effective and scientifically feasible method for addressing total phosphorus impairments. The methodology presented will serve as a model to conduct similar analyses in other watersheds. The Raritan River Basin was divided into 12 subwatershed‐based study areas. Point‐nonpoint source trading opportunities were examined for each study area by examining the point and nonpoint source total phosphorus loading to impaired water bodies. Of the 12 subwatersheds examined, four had a high potential for implementing a successful trading program. Since instream phosphorus concentrations are closely related to soil erosion, an additional analysis was performed to examine soil erodibility. Recommendations are presented for conducting an economic analysis following the feasibility study.     

REACH SCALE HYDRAULIC ASSESSMENT OF INSTREAM SALMONID HABITAT RESTORATION1

ABSTRACT: This study investigates the use of a two‐dimensional hydrodynamic model (River2D) for an assessment of the effects of instream large woody debris and rock groyne habitat structures. The bathymetry of a study reach (a side channel of the Chilliwack River located in southwestern British Columbia) was surveyed after the installation of 11 instream restoration structures. A digital elevation model was developed and used with a hydrodynamic model to predict local velocity, depth, scour, and habitat characteristics. The channel was resurveyed after the fall high‐flow season during which a bankfull event occurred. Pre‐flood and post‐flood bathymetry pool distributions were compared. Measured scour was compared to predicted shear and pre‐flood and post‐flood fish habitat indices for coho salmon (Oncorhynchus kisutch) and steelhead trout (O. mykiss) were compared. Two‐dimensional flow model velocity and depth predictions compare favorably to measured field values with mean standard errors of 24 percent and 6 percent, respectively, while areas of predicted high shear coincide with the newly formed pool locations. At high flows, the fish habitat index used (weighted usable area) increased by 150 percent to 210 percent. The application of the hydrodynamic model indicated a net habitat benefit from the restoration activities and provides a means of assessing and optimizing planned works.     

Sensitivity of Thermal Habitat of a Trout Stream to Potential Climate Change,Wisconsin, United States1

Abstract: Airborne thermal remote sensing from four flights on a single day from a single‐engine airplane was used to collect thermal infrared data of a 10.47‐km reach of the upper East Branch Pecatonica River in southwest Wisconsin. The study uses a one‐dimensional stream temperature model calibrated with the longitudinal profiles of stream temperature created from the four thermal imaging flights and validated with three days of continuous stream temperature data from instream data loggers on the days surrounding the thermal remote‐sensing campaign. Model simulations were used to quantify the sensitivity of stream thermal habitat to increases in air and groundwater temperature and changes in base flow. The simulations indicate that stream temperatures may reach critical maximum thresholds for brook trout (Salvelinus fontinalis) and brown trout (Salmo trutta) mortality, particularly if both air temperature increases and base flow declines. The approach demonstrates that thermal infrared data can greatly assist stream temperature model validation due to its high spatial resolution, and that this spatially continuous stream temperature data can be used to pinpoint spatial heterogeneity in groundwater inflow to streams. With this spatially distributed data on thermal heterogeneity and base‐flow accretion, stream temperature models considering various climate change scenarios are able to identify thermal refugia that will be critical for fisheries management under a changing climate.     

Improvement in HSPF's Low‐Flow Predictions by Implementation of a Power Law Groundwater Storage‐Discharge Relationship

We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km² drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.     

Characterizing a Major Urban Stream Restoration Project: Nine Mile Run (Pittsburgh,Pennsylvania, USA)

Urban stream restoration continues to be used as an ecological management tool, despite uncertainty about the long‐term sustainability and resilience of restored systems. Evaluations of restoration success often focus on specific instream indicators, with limited attention to the wider basin or parallel hydrologic and geomorphic process. A comprehensive understanding of urban stream restoration progress is particularly important for comparisons with nonurban sites as urban streams can provide substantial secondary benefits to urban residents. Here, we utilize a wide range of indicators to retrospectively examine the restoration of Nine Mile Run, a multi‐million dollar stream restoration project in eastern Pittsburgh (Pennsylvania, USA). Examination of available continuous hydrological data illustrates the high cost of failures to incorporate the data into planning and adaptive management. For example, persistent extreme flows drive geomorphic degradation threatening to reverse hydrologic connections created by the restoration and impact the improved instream biotic communities. In addition, human activities associated with restoration efforts suggest a positive feedback as the stream restoration has focused effort on the basin beyond the reach. Ultimately, urban stream restoration remains a potentially useful management tool, but continued improvements in post‐project assessment should include examination of a wider range of indicators.     

Effect of Flow Depth and Velocity on Nitrate Loss Rates in Natural Channels1

Carleton, James N. and Yusuf M. Mohamoud, 2012. Effect of Flow Depth and Velocity on Nitrate Loss Rates in Natural Channels. Journal of the American Water Resources Association (JAWRA) 1‐12. DOI: 10.1111/jawr.12007 Abstract: Loss rates of nitrate from streams and rivers are governed by movement of the ion from water column to anoxic bed sediments. Quantitative representations of nitrate in streams and rivers have often treated such losses as governed by first‐order mechanisms that are invariant with respect to potential modulating factors other than temperature. Results of studies in recent years, however, suggest that rates of water column‐sediment mass transfer are influenced by stream geometry and associated hydraulics. We develop expressions for the instream nitrate loss rate coefficient, k, as a function of water velocity and depth, using hydraulic geometry to empirically relate velocity to depth for two cases: (1) variability in mean conditions among reaches; and (2) temporal variability in conditions at a single reach, under changing flow. The result is expressions for k as functions of water column depth. Measured stream k values reported in the literature are shown to be well represented by expressions developed for the first case, and the potential for application to probabilistic analysis is briefly examined. We explore the latter case using the Hydrologic Simulation Program – FORTRAN (HSPF) model, modified to incorporate the dependence of k on instantaneous stream depth. In example simulations of two nitrate‐exporting watersheds, the incorporation of depth‐dependence of k produces improvement in the model’s ability to match observed stream nitrate concentrations.     

Critical Evaluation of How the Rosgen Classification and Associated “Natural Channel Design” Methods Fail to Integrate and Quantify Fluvial Processes and Channel Response1

Abstract: Over the past 10 years the Rosgen classification system and its associated methods of “natural channel design” have become synonymous to some with the term “stream restoration” and the science of fluvial geomorphology. Since the mid 1990s, this classification approach has become widely adopted by governmental agencies, particularly those funding restoration projects. The purposes of this article are to present a critical review, highlight inconsistencies and identify technical problems of Rosgen’s “natural channel design” approach to stream restoration. This paper’s primary thesis is that alluvial streams are open systems that adjust to altered inputs of energy and materials, and that a form‐based system largely ignores this critical component. Problems with the use of the classification are encountered with identifying bankfull dimensions, particularly in incising channels and with the mixing of bed and bank sediment into a single population. Its use for engineering design and restoration may be flawed by ignoring some processes governed by force and resistance, and the imbalance between sediment supply and transporting power in unstable systems. An example of how C5 channels composed of different bank sediments adjust differently and to different equilibrium morphologies in response to an identical disturbance is shown. This contradicts the fundamental underpinning of “natural channel design” and the “reference‐reach approach.” The Rosgen classification is probably best applied as a communication tool to describe channel form but, in combination with “natural channel design” techniques, are not diagnostic of how to mitigate channel instability or predict equilibrium morphologies. For this, physically based, mechanistic approaches that rely on quantifying the driving and resisting forces that control active processes and ultimate channel morphology are better suited as the physics of erosion, transport, and deposition are the same regardless of the hydro‐physiographic province or stream type because of the uniformity of physical laws.