Developing Nutrient Criteria for Streams: An Evaluation of the Frequency Distribution Method1

Abstract: The U.S. Environmental Protection Agency recommends two statistical methods to States and Tribes for developing nutrient criteria. One establishes a criterion as the 75th percentile of a reference‐population frequency distribution, the other uses the 25th percentile of a general‐population distribution; the U.S. Environmental Protection Agency suggests either method results in similar criteria. To evaluate each method, the Montana Department of Environmental Quality (MT DEQ) assembled data from STORET and other sources to create a nutrient general population. MT DEQ’s reference‐stream project provided reference population data. Data were partitioned by ecoregions, and by seasons (winter, runoff, and growing) defined for the project. For each ecoregion and season, nutrient concentrations at the 75th percentile of the reference population were matched to their corresponding concentrations in the general population. Additionally, nutrient concentrations from five regional scientific studies were matched to their corresponding reference population concentrations; each study linked nutrients to impacts on water uses. Reference‐to‐general population matches were highly variable between ecoregions, as nutrients at the 75th percentile of reference corresponded to percentiles ranging from the 4th to the 97th of the general population. In contrast, case studies‐to‐reference matches were more consistent, matching on average to the 86th percentile of reference, with a coefficient of variation of 13%.     

A Regional‐Scale Habitat Suitability Model to Assess the Effects of Flow Reduction on Fish Assemblages in Michigan Streams1

Zorn, Troy G., Paul W. Seelbach, and Edward S. Rutherford, 2012. A Regional‐Scale Habitat Suitability Model to Assess the Effects of Flow Reduction on Fish Assemblages in Michigan Streams. Journal of the American Water Resources Association (JAWRA) 48(5): 871‐895. DOI: 10.1111/j.1752‐1688.2012.00656.x Abstract: In response to concerns over increased use and potential diversion of Michigan’s freshwater resources, and the resulting state legislative mandate, an advisory council created an integrated assessment model to determine the potential for water withdrawals to cause an adverse resource impact to fish assemblages in Michigan’s streams. As part of this effort, we developed a model to predict how fish assemblages characteristic of different stream types would change in response to decreased stream base flows. We describe model development and use in this case study. The model uses habitat suitability information (i.e., catchment size, base‐flow yield, and July mean water temperature) for over 40 fish species to predict assemblage structure in an individual river segment under a range of base‐flow reductions. By synthesizing model runs for individual fish species at representative segments for each of Michigan’s 11 ecological stream types, we developed curves describing how typical fish assemblages in each type respond to flow reduction. Each stream type‐specific, fish response curve was used to identify streamflow reduction levels resulting in adverse resource impacts to characteristic fish populations, the regulatory standard. Used together with a statewide map of stream types, our model provided a spatially comprehensive framework for evaluating impacts of flow withdrawals on biotic communities across a diverse regional landscape.     

A Multi‐Agency Nutrient Dataset Used to Estimate Loads,Improve Monitoring Design,and Calibrate Regional Nutrient SPARROW Models1

Saad, David A., Gregory E. Schwarz, Dale M. Robertson, and Nathaniel L. Booth, 2011. A Multi‐Agency Nutrient Dataset Used to Estimate Loads, Improve Monitoring Design, and Calibrate Regional Nutrient SPARROW Models. Journal of the American Water Resources Association (JAWRA) 47(5):933‐949. DOI: 10.1111/j.1752‐1688. 2011.00575.x Abstract: Stream‐loading information was compiled from federal, state, and local agencies, and selected universities as part of an effort to develop regional SPAtially Referenced Regressions On Watershed attributes (SPARROW) models to help describe the distribution, sources, and transport of nutrients in streams throughout much of the United States. After screening, 2,739 sites, sampled by 73 agencies, were identified as having suitable data for calculating long‐term mean annual nutrient loads required for SPARROW model calibration. These sites had a wide range in nutrient concentrations, loads, and yields, and environmental characteristics in their basins. An analysis of the accuracy in load estimates relative to site attributes indicated that accuracy in loads improve with increases in the number of observations, the proportion of uncensored data, and the variability in flow on observation days, whereas accuracy declines with increases in the root mean square error of the water‐quality model, the flow‐bias ratio, the number of days between samples, the variability in daily streamflow for the prediction period, and if the load estimate has been detrended. Based on compiled data, all areas of the country had recent declines in the number of sites with sufficient water‐quality data to compute accurate annual loads and support regional modeling analyses. These declines were caused by decreases in the number of sites being sampled and data not being entered in readily accessible databases.     

Compensatory Mitigation for Streams Under the Clean Water Act: Reassessing Science and Redirecting Policy1

Doyle, Martin W. and F. Douglas Shields, 2012. Compensatory Mitigation for Streams Under the Clean Water Act: Reassessing Science and Redirecting Policy. Journal of the American Water Resources Association (JAWRA) 48(3): 494-509. DOI: 10.1111/j.1752-1688.2011.00631.x Abstract: Current stream restoration science is not adequate to assume high rates of success in recovering ecosystem functional integrity. The physical scale of most stream restoration projects is insufficient because watershed land use controls ambient water quality and hydrology, and land use surrounding many restoration projects at the time of their construction, or in the future, do not provide sufficient conditions for functional integrity recovery. Reach scale channel restoration or modification has limited benefits within the broader landscape context. Physical habitat variables are often the basis for indicating success, but are now increasingly seen as poor surrogates for actual biological function; the assumption “if you build it they will come” lacks support of empirical studies. If stream restoration is to play a continued role in compensatory mitigation under the United States Clean Water Act, then significant policy changes are needed to adapt to the limitations of restoration science and the social environment under which most projects are constructed. When used for compensatory mitigation, stream restoration should be held to effectiveness standards for actual and measurable physical, chemical, or biological functional improvement. To achieve improved mitigation results, greater flexibility may be required for the location and funding of restoration projects, the size of projects, and the restoration process itself.     

Concentration‐Duration‐Frequency Curves for Stream Turbidity: Possibilities for Assessing Biological Impairment1

Abstract: Siltation and subsequent biological impairment is a national problem prompting state regulatory agencies to develop sediment total maximum daily loads (TMDL) for many streams. To support TMDL targets for reduced sediment yield in disturbed watersheds, a critical need exists for stream assessments to identify threshold concentrations of suspended sediment that impact aquatic biota. Because of the episodic nature of stream sediment transport, thresholds should not only be a function of sediment concentration, but also of duration and dose frequency. Water quality sondes can collect voluminous amounts of turbidity data, a surrogate for suspended sediment, at intervals that can be used to characterize concentration, duration, and frequency of elevated turbidity events. To characterize turbidity sonde data in an ecologically relevant manner, a methodology for concentration‐duration‐frequency (CDF) curves was developed using turbidity doses that relate to different levels of biological impairment. To illustrate this methodology, turbidity CDF curves were generated for two sites on Little Pigeon River in the Great Smoky Mountains National Park, Tennessee, using over 30,000 sonde data measurements per site for a one‐year period. Utilizing a Poisson arrival approach, turbidity spikes were analyzed stochastically by observing the frequency and duration of recorded events over a turbidity level that relates to a biological dose response. An exponential equation was used to fit duration and frequency of a specified turbidity level to generate concentric‐shaped CDF curves, where at specific turbidities longer durations occurred less frequently and conversely shorter durations occurred more frequently. The significance of the equation fit to the data was accomplished with a Kolmogorov‐Smirnov goodness‐of‐fit test. Our findings showed that the CDF curves derived by an exponential function performed reasonable well, with most curves significant at a 95% confidence level. These CDF curves were then used to demonstrate how they could be used to assess biological impairment, and identify future research needs for improved development of sediment TMDLs.     

Tributary Plunging in an Urban Lake (Onondaga Lake): Drivers,Signatures, and Implications1

Abstract: A combination of long‐term fixed‐frequency and robotic monitoring information for a polluted urban lake, Onondaga Lake, New York, and two of its tributaries is used to resolve the propensity for, and occurrences of, tributary plunging. Cooler temperatures (T) and higher salinity (S) are primarily responsible for the elevated density and plunging of one of these polluted streams for the summer through early fall interval. In‐lake transport of this plunging tributary, which receives inputs from combined sewer overflows (CSOs), is tracked by its high S during dry weather, its high turbidity (T_n) with associated lower S (dilution with rainwater) following runoff events, and by its characteristic ionic composition. These signatures are documented extending from the creek mouth, through a connecting navigation channel, through the inflow zone of the lake, and into metalimnetic depths of pelagic portions of the lake. The entry of this polluted tributary below the depth interval(s) of primary production and contact recreation has important implications for the ongoing major rehabilitation program for this lake. The plunging phenomenon diminishes the benefits previously expected for related features of the lake’s water quality from ongoing management efforts to abate CSO inputs and reduce nonpoint nutrient loading from the tributary. Previously this tributary tended to instead enter the upper layers of the lake during the operation of an adjoining soda ash manufacturing facility (closure in 1986), as a result of high lake S caused by the industry’s ionic waste discharge.     

Nutrient Bioassays of Growth Parameters for Algae in the North Bosque River of Central Texas1

Abstract: Nutrient dose‐response bioassays were conducted using water from three sites along the North Bosque River. These bioassays provided support data for refinement of the Soil and Water Assessment Tool (SWAT) model used in the development of two phosphorus TMDLs for the North Bosque River. Test organisms were native phytoplanktonic algae and stock cultured Pseudokirchneriella subcapitata (Korshikov) Hindak. Growth was measured daily by in vivo fluorescence. Algal growth parameters for maximum growth (μ_max) and half‐saturation constants for nitrogen (K_N) or phosphorus (K_P) were determined by fitting maximum growth rates associated with each dose level to a Monod growth rate function. Growth parameters of native algae were compared between locations and to growth parameters of P. subcapitata and literature values. No significant differences in half‐saturation constants were indicated within nutrient treatment for site or algal type. Geometric mean K_N was 32 μg/l and for K_P 7 μg/l. A significant difference was detected in maximum growth rates between algae types but not between sites or nutrient treatments. Mean μ_max was 1.5/day for native algae and 1.2/day for stock algae. These results indicate that watershed‐specific maximum growth rates may need to be considered when modeling algal growth dynamics with regard to nutrients.     

Assessing Impervious Surface Connectivity and Applications for Watershed Management1

Abstract: Although total impervious area (TIA) is often used as an indicator of urban disturbance, recent studies suggest that the subset of impervious surfaces that route stormwater runoff directly to streams via stormwater pipes, called directly connected impervious area (DCIA), may be a better predictor of stream ecosystem alteration. We evaluated the differences between TIA and DCIA in the Shepherd Creek catchment, a small (1.85‐km²), suburban basin in Cincinnati, Ohio. Imperviousness determinations were calculated based on publicly available geographic information system (GIS) data and parcel‐scale field assessments, and these direct assessments were compared to DCIA calculated from published, empirical relationships. Impervious and semi‐impervious area comprised 13.1% of the catchment area, with 56.3% of the impervious area connected. When summarized by subcatchments (0.26‐1.85 km²), TIA measured in the field (11‐23%) was considerably higher than that calculated from the National Land Cover Data Imperviousness Layer (7‐18%). In contrast, TIA calculated based on aerial photos was similar to TIA calculated from field assessments, thus indicating that photo interpretation may be adequate for catchment‐scale (>25 ha) TIA determinations. While these GIS data sources can be used to calculate TIA, on‐site assessments were necessary to accurately determine DCIA within residential parcels. There was a wide variation in percent connectivity across parcels, and, subsequently, DCIA was not accurately predicted from empirical relationships with TIA. We discuss applications of DCIA data that highlight the importance of parcel‐scale field assessments for managing suburban watersheds.     

Modeling Biochemical Oxygen Demand Through the Middle and Lower Savannah River1

Abstract: In order to improve modeling accuracy and general understanding of lotic biochemical oxygen demand (BOD), this study characterized river metabolism with the current Georgia Environmental Protection Division method for the middle and lower Savannah River basin (MLSRB) and several alternative methods developed with 120‐day, long‐term biochemical oxygen demand (LTBOD) data from the MLSRB. The data were a subset of a larger two‐year LTBOD study to characterize and understand BOD in the MLSRB, located approximately between Augusta, Georgia, and Savannah, Georgia, along the border of Georgia and South Carolina. The LTBOD data included total oxygen loss and nitrogen speciation for separately quantifying nitrification. Results support the following insights and opportunities for modeling methods: (1) it is important to modeling accuracy that residuals be checked for even dispersion to avoid areas of over‐ and underprediction; (2) modeling with bounded, yet unfixed, rates is a sufficiently simple alternative to fixed‐rate modeling that can eliminate the need for manual adjustments and provide additional system understanding to inform regulation; (3) if fixed rates modeling is desired, model quality for this system might be improved through revising the current low rate (along with the associated f‐ratio updates) from 0.02/day rate to 0.006/day and potentially adding a new rate at 1.0/day in some cases; and (4) the current 57/43 ratio of slow/fast BOD is reasonable based on the 52/45/3 slow/fast/faster BOD proportions of this study.     

Assessment Tools for Urban Catchments: Developing Stressor Gradients1

Abstract: This is the first in a series of three articles designed to establish empirically defined biological indicators and thresholds for impairment for urbanized catchments, and to describe a process by which the biological condition of waterbodies in urbanized catchments can be applied. This article describes alternative gradients of urbanization for assessing and selecting a nationally applicable biological index (article 2 – Purcell et al., this issue ) and defining the potential of biological communities within a gradient of cumulative stressors (article 3 – Paul et al. this issue ). Gradients were designed to represent the most prominent mosaic of stressors found in urban settings. A primary urban gradient was assembled based on readily obtained information of urbanization to include three broad‐scale parameters: percent urban land use/land cover, population density, and road density. This gradient was used as the standard by which alternative urban gradients, which included fine‐scale instream chemical and hydrologic parameters, were assessed. Five alternative gradients were developed to provide numerous environmental management options based on availability of data from water program resources. The urban gradients were developed with the intent that they be applied throughout the country; therefore, data from three different regions of the United States (Mid‐Atlantic, Midwest, and Pacific Coast) were used to validate the urban gradient model. Our study showed that a relatively straightforward stressor gradient consisting of human population density, road density, and urban land use is useful in providing a framework for developing relevant biological indicators and evaluating the potential of biological communities as a basis for assessing attainment of designated aquatic life use.     

Nutrient Sources and Transport in the Missouri River Basin,with Emphasis on the Effects of Irrigation and Reservoirs1

Brown, Juliane B., Lori A. Sprague, and Jean A. Dupree, 2011. Nutrient Sources and Transport in the Missouri River Basin, With Emphasis on the Effects of Irrigation and Reservoirs. Journal of the American Water Resources Association (JAWRA) 47(5):1034‐1060. DOI: 10.1111/j.1752‐1688.2011.00584.x Abstract: SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were used to relate instream nutrient loads to sources and factors influencing the transport of nutrients in the Missouri River Basin. Agricultural inputs from fertilizer and manure were the largest nutrient sources throughout a large part of the basin, although atmospheric and urban inputs were important sources in some areas. Sediment mobilized from stream channels was a source of phosphorus in medium and larger streams. Irrigation on agricultural land was estimated to decrease the nitrogen load reaching the Mississippi River by as much as 17%, likely as a result of increased anoxia and denitrification in the soil zone. Approximately 16% of the nitrogen load and 33% of the phosphorus load that would have otherwise reached the Mississippi River was retained in reservoirs and lakes throughout the basin. Nearly half of the total attenuation occurred in the eight largest water bodies. Unlike the other major tributary basins, nearly the entire instream nutrient load leaving the outlet of the Platte and Kansas River subbasins reached the Mississippi River. Most of the larger reservoirs and lakes in the Platte River subbasin are upstream of the major sources, whereas in the Kansas River subbasin, most of the source inputs are in the southeast part of the subbasin where characteristics of the area and proximity to the Missouri River facilitate delivery of nutrients to the Mississippi River.     

Land Use,Geology, Enrichment,and Stream Biota in the Eastern Ridge and Valley Ecoregion: Implications for Nutrient Criteria Development1

Abstract: The eastern panhandle region of West Virginia is entirely within the Appalachian Ridge and Valley ecoregion. It is underlain by limestone in the eastern part and by shale and sandstone in the western part. Agricultural and urban development has affected the condition of the streams of this region. We examined samples from 165 stations in the Ridge and Valley, collected from 1998 to 2004. Land use, geological characteristics, physical and chemical parameters, and algal and macroinvertebrate assemblages were used to identify potential stressors that affect streams in the region. Our analyses indicated that both human land uses and ecoregional differences led to elevated nutrient concentrations in streams of the study areas. Multiple regression analyses indicated that both agricultural and urban land use in the watershed were associated with high nutrient concentrations (NO₂₊₃, total nitrogen, and total phosphorus) in streams. These elevated nutrient concentrations have led to increased algal biomass, increased trophic state, and degradation of macroinvertebrate community in the streams. Values of the West Virginia Stream Condition Index, as well as several other benthic macroinvertebrate metrics, decreased with increased nutrient concentrations and conductivity, especially in the limestone region. When regional differences were partitioned out in the analysis, nutrient concentrations became the strongest stressor in the limestone region while conductivity exhibited less of an effect on macroinvertebrate metrics. Meanwhile, periphyton diagnostic metrics also responded to increased nutrient concentrations, suggesting nutrients could be a cause of biological degradation in the Eastern Ridge and Valley region. Multiple approaches and multiple lines of evidence (reference approach and stressor‐response approach) were applied to develop nutrient benchmarks for different geological regions in the study watershed.     

Observations and Modeling of Stream Plunging in an Urban Lake1

Owens, Emmet M., Steven W. Effler, Anthony R. Prestigiacomo, David A. Matthews, and Susan M. O’Donnell, 2012. Observations and Modeling of Stream Plunging in an Urban Lake. Journal of the American Water Resources Association (JAWRA) 48(4): 707‐721. DOI: 10.1111/j.1752‐1688.2012.00646.x Abstract: The plunging behavior of two tributaries in Onondaga Lake, New York, is quantified based on a program of monitoring, process studies, and modeling. The dynamics of buoyancy of the tributaries are resolved with hourly measurements of temperature (T), specific conductance (SC), and turbidity (Tn) at the mouths, and observations every 6 h in the lake. Negative buoyancy of the tributaries is found to diminish and change rapidly during runoff events compared to dry periods. In‐lake patterns of the transport of plunging inflow are resolved for dry weather conditions using a dye tracer, and following a runoff event through measurements of T, SC, and Tn. The hydrodynamic/transport model ELCOM (Estuary Lake and Coastal Ocean Model) is demonstrated to perform well in simulating these patterns. Analyses conducted with the model establish the importance of diurnal effects and in‐lake mixing mediated by wind, the need for temporally detailed measurements during runoff events, and modifications of the plunging behavior of the urban tributary as it passes through a harbor. The model provides critical information to support rehabilitation programs for the lake by quantifying the transport of the two largest tributaries, particularly the distribution of the loads between the upper waters and stratified layers. The model predicts that 10% of the urban tributary load enters the upper waters of the lake within 24 h for a dry weather period; this portion increases to 30% for a runoff event.