Using benefit transfer to estimate housing value increases from improved water clarity: A case study of lakes in Kosciusko County, Indiana
This study provides step-by-step guidance for practitioners and local stakeholders on how to use existing study results to conduct benefit transfer (BT), and ultimately make informed predictions of how improvements in lake water clarity may benefit surrounding communities. The procedures are demonstrated using a publicly available meta-dataset developed by the United States Environmental Protection Agency, and a subsequent meta-analysis that synthesizes the literature on how improvements in water clarity impact home values. The BT procedures are demonstrated using a case study of 14 large lakes in Kosciusko County, Indiana. Lake-specific average increases in home values, as well as the value of the housing stock in aggregate, are calculated for illustrative improvements in lake water clarity. This analysis provides a critical bridge to better connect high-quality, academic research with real-world policy analysis, and ultimately serves to better equip local governments and stakeholders to make more informed policy and land use decisions.
Spatial and Temporal Variability in Stream Thermal Regime Drivers for Three River Networks During the Summer Growing Season
Many cold-water dependent aquatic organisms are experiencing habitat and population declines from increasing water temperatures. Identifying mechanisms which drive local and regional stream thermal regimes facilitates restoration at ecologically relevant scales. Stream temperatures vary spatially and temporally both within and among river basins. We developed a modeling process to identify statistical relationships between drivers of stream temperature and covariates representing landscape, climate, and management-related processes. The modeling process was tested in 3 study areas of the Pacific Northwest USA during the growing season (May [start], August [warmest], September [end]). Across all months and study systems, covariates with the highest relative importance represented the physical landscape (elevation [1], catchment area [3], main channel slope [5]) and climate covariates (mean monthly air temperature [2] and discharge [4]). Two management covariates (ground water use [6] and riparian shade [7]) also had high relative importance. Across the growing season (for all basins) local reach slope had high relative importance in May, but transitioned to a regional main channel slope covariate in August and September. This modeling process identified regionally similar and locally unique relationships among drivers of stream temperature. High relative importance of management-related covariates suggested potential restoration actions for each system.
Random Forest models to estimate bankfull and low flow channel widths and depths across the conterminous United States
Channel dimensions (width and depth) at varying flows influence a host of instream ecological processes, as well as habitat and biotic features; they are a major consideration in stream habitat restoration and instream flow assessments. Models of widths and depths are often used to assess climate change vulnerability, develop endangered species recovery plans, and model water quality. However, development and application of such models require specific skillsets and resources. To facilitate acquisition of such estimates, we created a dataset of modeled channel dimensions for perennial stream segments across the conterminous U.S. We used random forest models to predict wetted width, thalweg depth, bankfull width, and bankfull depth from several thousand field measurements of the National Rivers and Streams Assessment. Observed channel widths varied from <5 m to >2000 m and depths varied from <2 m to >125 m. Metrics of watershed area, runoff, slope, land use, and more were used as model predictors. The models had high pseudo R-squared values (0.70 to 0.91) and median absolute errors within ±6% to ±21% of the interquartile range of measured values across ten stream orders. Predicted channel dimensions can be joined to 1.1 million stream segments of the 1:100K resolution National Hydrography Dataset Plus (version 2.1). These predictions, combined with a rapidly growing body of nationally available data, will further enhance our ability to study and protect aquatic resources.
Wetland Flowpaths Mediate Nitrogen and Phosphorus Concentrations across the Upper Mississippi River Basin
Eutrophication, harmful algal blooms, and human health impacts are critical environmental challenges resulting from excess nitrogen and phosphorus in surface waters. Yet we have limited information regarding how wetland characteristics mediate water quality across watershed scales. We developed a large, novel set of spatial variables characterizing hydrological flowpaths from wetlands to streams, that is, "wetland hydrological transport variables," to explore how wetlands statistically explain the variability in total nitrogen (TN) and total phosphorus (TP) concentrations across the Upper Mississippi River Basin (UMRB) in the United States. We found that wetland flowpath variables improved landscape-to-aquatic nutrient multilinear regression models (from = 0.89 to 0.91 for TN; = 0.53 to 0.84 for TP) and provided insights into potential processes governing how wetlands influence watershed-scale TN and TP concentrations. Specifically, flowpath variables describing flow-attenuating environments, for example, subsurface transport compared to overland flowpaths, were related to lower TN and TP concentrations. Frequent hydrological connections from wetlands to streams were also linked to low TP concentrations, which likely suggests a nutrient source limitation in some areas of the UMRB. Consideration of wetland flowpaths could inform management and conservation activities designed to reduce nutrient export to downstream waters.
Spatial analysis of future climate risk to stormwater infrastructure
Climate change is expected to result in more intense precipitation events that will affect the performance and design requirements of stormwater infrastructure. Such changes will vary spatially, and climate models provide a range of estimates of the effects on events of different intensities and recurrence. Infrastructure performance should be evaluated against the expected range of events, not just rare extremes. We present a national-scale, spatially detailed screening assessment of the potential effects of climatic change on precipitation, stormwater runoff, and associated design requirements. This is accomplished through adjustment relative to multiple future climate scenarios of precipitation intensity-duration-frequency analyses presented in NOAA Atlas 14, which are commonly used in infrastructure design. Future precipitation results are estimated for each Atlas 14 station (these currently omit the Pacific Northwest). Results are interpolated using a geographically conditioned regression kriging approach to provide information about potential climate change impacts in a format more directly useful to local stormwater managers. The intensity of 24-h events with 2-year or greater recurrence is likely to increase in most areas of the United States leading to increased runoff and potential need for increased storage volumes. Changes in more frequent events (e.g., the 90th percentile event) commonly used in design of green infrastructure are relatively less.
Considerations for Using Alternative Technologies in Nutrient Management on Cape Cod: Beyond Cost and Performance
Mitigating non-point source nitrogen in coastal estuaries is economically, environmentally, logistically, and socially challenging. On Cape Cod, Massachusetts, nitrogen management includes both traditional, centralized wastewater treatment and sewering as well as a number of alternative technologies. We conducted semi-structured interviews with 37 participants from governmental and non-governmental organizations as well as related industries to identify the barriers and opportunities for the use of alternative technologies to mitigate nitrogen pollution. The interviews were recorded, transcribed, and then analyzed using content analysis and rhetorical analysis. Cost and technical capacity to reduce nitrogen were the most discussed considerations. Beyond those, there were a slew of additional considerations that also impacted whether a technology would be installed, permitted, and socially accepted. These included: maintenance and monitoring logistics, comparisons to sewering, co-benefits, risk/uncertainty, community culture, extent of public engagement, permitting/regulatory challenges, and siting considerations. The insights about these additional considerations are valuable for transferring to other coastal areas managing nutrient impairments that may have not yet factored in these considerations when making decisions about how to meet water quality goals.
River Basin Export Reduction Optimization Support Tool; a tool to screen options for reducing nutrient loads while minimizing cost
Excess loading of nitrogen and phosphorus to river networks causes environmental harm, but reducing loads from large river basins is difficult and expensive. We develop a new tool, the River Basin Export Reduction Optimization Support Tool (RBEROST) to identify least-cost combinations of management practices that will reduce nutrient loading to target levels in downstream and mid-network waterbodies. We demonstrate the utility of the tool in a case study in the Upper Connecticut River basin in New England, USA. The total project cost of optimized lowest-cost plans ranged from $18.0 million to $41.0 million per year over 15 years depending on user specifications. Plans include both point source and non-point source management practices, and most costs are associated with urban stormwater practices. Adding a 2% margin of safety to loading targets improved estimated probability of success from 37.5% to 99%. The large spatial scale of RBEROST, and the consideration of both point and non-point source contributions of nutrients, makes it well suited as an initial screening tool in watershed planning.
Early Influence of the COVID-19 Pandemic on Volunteer Water Monitoring Programs in the United States and Canada
Volunteer water monitoring programs generate new scientific knowledge, contribute data to decision-making processes, and increase social networks, technical knowledge, and skills of participants. Declaration of the COVID-19 pandemic threatened the ability of these programs to continue to engage volunteers to achieve such outcomes. A national water monitoring network hosted a brainstorming webinar to facilitate communication across programs to identify potential solutions to pandemic-influenced challenges. Following that webinar, a survey of United States and Canadian volunteer monitoring programs that was conducted about 3 months into the pandemic revealed that 72% of 80 responding programs planned to carry on through the 2020 field season despite most having experienced delayed starts. Other common program modifications implemented in the first months of the pandemic included adding COVID-19 safety information to program guidance, changing field team composition, monitoring timing and logistics, and adopting new communications strategies. Most programs reported loss or anticipated loss in number of data observations (74%) and volunteers (66%), while 44% reported known or anticipated losses in funding. Seventeen percent of responding programs were able to swiftly develop distance learning tools to train participants, which led to increased program capacity to reach broader audiences.
Assessing Evidence of Phosphorus Concentration Trends in North American Fresh Waters
The U.S. EPA's National Aquatic Resource Surveys (NARS) documented evidence of widespread, unexplained total phosphorus (TP) concentration increases in lakes and streams across the United States during the 2000 - 2012 time period. To examine the robustness of evidence for this trend, we used additional monitoring datasets to calculate rates of TP change in thousands of individual waterbodies across the U.S. during the same time frame, and compared them against TP change rates calculated in the same manner for waterbodies that were resurveyed under NARS in different years. For the additional datasets, median rates of TP change were substantially lower than median rates calculated using NARS data. To further examine differences between NARS and non-NARS results in specific waterbodies, we assembled composite datasets for 52 predominantly northern lakes that by chance had been sampled under both NARS and other sampling programs during the same time frame. Using only NARS data, the median calculated TP change rate for this set of lakes was positive, and similar to that for the larger set of 401 resurveyed NARS lakes. However, when additional sample data were included, the median calculated TP change rate for these lakes was much lower. Results suggest that increasing TP concentrations in waterbodies may not have been as ubiquitous as suggested. They also illustrate a need to supplement randomized continental-scale monitoring with detailed, site-focused investigations.
Methods for Estimating Locations of Housing Units Served by Private Domestic Wells in the United States Applied to 2010
In 1990, the last time the decennial census included a question on domestic drinking water source, it was estimated that private domestic water wells (PDWs) supplied household water to about 15.1 million housing units (15% of the population) in the United States (U.S.). PDWs are not regulated by the Safe Drinking Water Act, and with few exceptions, are not subject to the water quality testing required of public water suppliers. We expanded two methods in estimating housing units reliant on PDWs from an Oklahoma pilot study (Weaver et al. 2017), nationally. Both use 1990 census data on drinking water sources as a baseline. The first method uses housing unit change and private well drilling logs for 20 states. This allows for the rate of well use to change between 1990 and 2010 in these states. The second, based solely on housing unit change, assumes a constant rate of well use. Ordinary least squares regression demonstrated ( = 0.78) that the methods yield similar estimates for nationwide well use. Using the housing unit change method, it is estimated that in 2010, 23 million housing units were reliant on PDWs (17% of the population). We provide these estimates at the census block group and census block resolution. This dataset will assist in a better understanding of the reliance on PDWs in the U.S., and position local, tribal, state, and national groups to better protect this water resource from contaminant sources.
Monitoring Turbidity in San Francisco Estuary and Sacramento-San Joaquin Delta Using Satellite Remote Sensing
This study utilizes satellite data to investigate water quality conditions in the San Francisco Estuary and its upstream delta, the Sacramento-San Joaquin River Delta. To do this, this study derives turbidity from the European Space Agency satellite Sentinel-2 acquired from September 2015 to June 2019 and conducts a rigorous validation with in situ measurements of turbidity from optical sensors at continuous monitoring stations. This validation includes 965 matchup comparisons between satellite and in situ sensor data across 22 stations, yielding = 0.63 and 0.75 for Nephelometric Turbidity Unit and Formazin Nephelometric Unit (FNU) stations, respectively. This study then applies remote sensing to evaluate patterns in turbidity during the Suisun Marsh Salinity Control Gates Action ("Gates action"), a pilot study designed to increase habitat access and quality for the endangered Delta Smelt. The basic strategy was to direct more freshwater into Suisun Marsh, creating more low salinity habitat that would then have higher (and more suitable) turbidity than upstream river channels. For all seven acquisitions considered from June 29 to September 27, 2018, turbidity conditions in Bays and Sloughs subregions were consistently higher (and more suitable) (26-47 FNU) than what was observed in the upstream River region (13-25 FNU). This overall pattern was observed when comparing images acquired during similar tidal stages and heights.
in an Urban World: New Contexts for Hydraulic Conductivity
Hydraulic conductivity () is a key hydrologic parameter widely recognized to be difficult to estimate and constrain, with little consistent assessment in disturbed, urbanized soils. To estimate , it is either measured, or simulated by pedotransfer functions, which relate to easily measured soil properties. We measured in urbanized soils by double-ring infiltrometer ( ), near-saturated tension infiltrometry ( ), and constant head borehole permeametry ( ), along with other soil properties across the major soil orders in 12 United States cities. We compared measured with that predicted from the pedotransfer function, ROSETTA. We found that regardless of soil texture, was consistently larger than ; with the latter having slightly less sample variance. was dependent upon specific subsurface conditions, and contrary to common expectations, did not always decrease with depth. Based on either soil textural class, or percent textural separates (sand, silt clay), ROSETTA did not accurately predict measured for surface nor subsurface soils. We go on to discuss how varies in urban landscapes, the role of measurement methods and artifacts in the perception of this metric, and implications for hydrologic modeling. Overall, we aim to inspire consistency and coherence when addressing -related challenges in sustainable urban water management.
A landscape assessment and associated dataset of stream confluences for the conterminous U.S
Stream confluences are important components of fluvial networks. Hydraulic forces meeting at stream confluences often produce changes in streambed morphology and sediment distribution. These changes often increase habitat heterogeneity relative to upstream and downstream locations, which have led some to identify them as biological hotspots. Despite their potential ecological importance, there are relatively few empirical studies documenting ecological patterns upstream and downstream of confluences. We have produced a publicly available dataset of stream confluences and associated watershed attributes for the conterminous USA. The dataset includes 1,085,629 stream confluences and 383 attributes for each confluence organized into 15 dataset tables for both tributary and mainstem upstream catchments and watersheds. Themes in the dataset include hydrology (e.g., stream order), land cover, land cover change, geology (e.g., calcium content of underlying lithosphere), physical condition (e.g., precipitation), measures of ecological integrity, and stressors (e.g., impaired streams). Additionally, we used measures of ecological integrity to assess the condition of the stream confluences. Aside from a generally positive east-to-west gradient in ecological condition, we found that approximately one-third of the confluences had markedly contrasting ecological conditions between mainstem and tributary, catchment and watershed, or both. The dataset should support many, multifaceted studies of stream confluence ecology.
When, where, and how to intervene? Trade-offs between time and costs in coastal nutrient management
Policies and regulations designed to address nutrient pollution in coastal waters are often complicated by delays in environmental and social systems. Social and political inertia may delay implementation of cleanup projects, and even after the best nutrient pollution management practices are developed and implemented, long groundwater travel times may delay the impact of inland or upstream interventions. These delays and the varying costs of nutrient removal alternatives used to meet water quality goals combine to create a complex dynamic decision problem with trade-offs about when, where, and how to intervene. We use multi-objective optimization to quantify the trade-offs between costs and minimizing the time to meet in-bay nutrient reduction goals represented as a Total Maximum Daily Load (TMDL). We calculate the impact of using in-bay (in-situ) nutrient removal through shellfish aquaculture relative to waiting for traditional source control to be implemented. We apply these methods to the Three Bays Watershed in Cape Cod, Massachusetts. In gross benefit terms, not accounting for any social costs, this equates to an average value of 37¢ (2035 TMDL target date) and 11¢ (2060 TMDL target date) per animal harvested over the plan implementation period. Our results encourage the consideration of alternative and in-situ approaches to tackle coastal pollution while traditional source control is implemented and its effects realized over time.
Lake Water Levels and Associated Hydrologic Characteristics in the Conterminous U.S
Establishing baseline hydrologic characteristics for lakes in the U.S. is critical to evaluate changes to lake hydrology. We used the U.S. EPA National Lakes Assessment 2007 and 2012 surveys to assess hydrologic characteristics of a population of ~45,000 lakes in the conterminous U.S. based on probability samples of ~1,000 lakes/yr distributed across nine ecoregions. Lake hydrologic study variables include water-level drawdown (i.e., vertical decline and horizontal littoral exposure) and two water stable isotope-derived parameters: evaporation-to-inflow (E:I) and water residence time. We present 1) national and regional distributions of the study variables for both natural and man-made lakes and 2) differences in these characteristics between 2007 and 2012. In 2007, 59% of the population of U.S. lakes had or drawdown relative to water levels in ecoregional reference lakes with minimal human disturbances; while in 2012, only 20% of lakes were significantly drawn down beyond normal ranges. Water isotope-derived variables did not differ significantly between survey years in contrast to drawdown. Median E:I was 20% indicating that flow-through processes dominated lake water regimes. For 75% of U.S. lakes, water residence time was < 1 year and was longer in natural vs. man-made lakes. Our study provides baseline ranges to assess local and regional lake hydrologic status and inform management decisions in changing environmental conditions.
Comparison and Evaluation of Gridded Precipitation Datasets in a Kansas Agricultural Watershed Using SWAT
Gridded precipitation datasets are becoming a convenient substitute for gauge measurements in hydrological modeling; however, these data have not been fully evaluated across a range of conditions. We compared four gridded datasets (Daily Surface Weather and Climatological Summaries [DAYMET], North American Land Data Assimilation System [NLDAS], Global Land Data Assimilation System [GLDAS], and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) as precipitation data sources and evaluated how they affected hydrologic model performance when compared with a gauged dataset, Global Historical Climatology Network-Daily (GHCN-D). Analyses were performed for the Delaware Watershed at Perry Lake in eastern Kansas. Precipitation indices for DAYMET and PRISM precipitation closely matched GHCN-D, whereas NLDAS and GLDAS showed weaker correlations. We also used these precipitation data as input to the Soil and Water Assessment Tool (SWAT) model that confirmed similar trends in streamflow simulation. For stations with complete data, GHCN-D based SWAT-simulated streamflow variability better than gridded precipitation data. During low flow periods we found PRISM performed better, whereas both DAYMET and NLDAS performed better in high flow years. Our results demonstrate that combining gridded precipitation sources with gauge-based measurements can improve hydrologic model performance, especially for extreme events.
Testing of the Storm Water Management Model Low Impact Development Modules
Stormwater infrastructure designers and operators rely heavily on the United States Environmental Protection Agency's Storm Water Management Model (SWMM) to simulate stormwater and wastewater infrastructure performance. Since its inception in the late 1970s, improvements and extensions have been tested and evaluated rigorously to verify the accuracy of the model. As a continuation of this progress, the main objective of this study was to quantify how accurately SWMM simulates the hydrologic activity of low impact development (LID) storm control measures. Model performance was evaluated by quantitatively comparing empirical data to model results using a multievent, multiobjective calibration method. The calibration methodology utilized the PEST software, a Parameter ESTimation tool, to determine unmeasured hydrologic parameters for SWMM's LID modules. The calibrated LID modules' Nash-Sutcliffe efficiencies averaged 0.81; average percent bias (PBIAS) -9%; average ratio of root mean square error to standard deviation of measured values 0.485; average index of agreement 0.94; and the average volume error, simulated vs. observed, was +9%. SWMM accurately predicted the timing of peak flows, but usually underestimated their magnitudes by 10%. The average volume reduction, measured outflow volume divided by inflow volume, was 48%. We had more difficulty in calibrating one study, an infiltration trench, which identified a significant limitation of the current version of the SWMM LID module; it cannot simulate lateral exfiltration of water out of the storage layers of a LID storm control measure. This limitation is especially severe for a deep LIDs, such as infiltration trenches. Nevertheless, SWMM satisfactorily simulated the hydrologic performance of eight of the nine LID practices.
Full Water-Cycle Monitoring in an Urban Catchment Reveals Unexpected Water Transfers (Detroit MI, USA)
A goal in urban water management is to reduce the volume of stormwater runoff in urban systems and the effect of combined sewer overflows into receiving waters. Effective management of stormwater runoff in urban systems requires an accounting of various components of the urban water balance. To that end, precipitation, evapotranspiration, sewer flow, and groundwater in a 3.40-hectare sewershed in Detroit, Michigan were monitored to capture the response of the sewershed to stormwater flow prior to implementation of stormwater control measures. Monitoring results indicate that stormflow in sewers was not initiated unless rain depth was 3.6 mm or greater. Evapotranspiration removed more than 40 percent of the precipitation in the sewershed whereas pipe flow accounted for 19 to 85 percent of the losses. Flows within the sewer that could not be associated with direct precipitation indicate an unexpected exchange of water between the leaky sewer and the groundwater system, pathways through abandoned or failing residential infrastructure, or a combination of both. Groundwater data indicate that groundwater flows into the leaky combined sewer rather than out. This research demonstrates that urban hydrologic fluxes can modulate the local water cycle in complex ways which affect the efficiency of the wastewater system, effectiveness of stormwater management, and, ultimately, public health.
Differing Modes of Biotic Connectivity within Freshwater Ecosystem Mosaics
We describe a collection of aquatic and wetland habitats in an inland landscape, and their occurrence within a terrestrial matrix, as a "freshwater ecosystem mosaic" (FEM). Aquatic and wetland habitats in any FEM can vary widely, from permanently ponded lakes, to ephemerally ponded wetlands, to groundwater-fed springs, to flowing rivers and streams. The terrestrial matrix can also vary, including in its influence on flows of energy, materials, and organisms among ecosystems. Biota occurring in a specific region are adapted to the unique opportunities and challenges presented by spatial and temporal patterns of habitat types inherent to each FEM. To persist in any given landscape, most species move to recolonize habitats and maintain mixtures of genetic materials. Species also connect habitats through time if they possess needed morphological, physiological, or behavioral traits to persist in a habitat through periods of unfavorable environmental conditions. By examining key spatial and temporal patterns underlying FEMs, and species-specific adaptations to these patterns, a better understanding of the structural and functional connectivity of a landscape can be obtained. Fully including aquatic, wetland, and terrestrial habitats in FEMs facilitates adoption of the next generation of individual-based models that integrate the principles of population, community, and ecosystem ecology.
Managing Uncertainty in Runoff Estimation with the U.S. Environmental Protection Agency National Stormwater Calculator
The U.S. Environmental Protection Agency National Stormwater Calculator (NSWC) simplifies the task of estimating runoff through a straightforward simulation process based on the EPA Stormwater Management Model. The NSWC accesses localized climate and soil hydrology data, and options to experiment with low-impact development (LID) features for parcels up to 5 ha in size. We discuss how the NSWC treats the urban hydrologic cycle and focus on the estimation uncertainty in soil hydrology and its impact on runoff simulation by comparing field-measured soil hydrologic data from 12 cities to corresponding NSWC estimates in three case studies. The default NSWC hydraulic conductivity is 10.1 mm/h, which underestimates conductivity measurements for New Orleans, Louisiana (95 ± 27 mm/h) and overestimates that for Omaha, Nebraska (3.0 ± 1.0 mm/h). Across all cities, the NSWC prediction, on average, underestimated hydraulic conductivity by 10.5 mm/h compared to corresponding measured values. In evaluating how LID interact with soil hydrology and runoff response, we found direct hydrologic interaction with pre-existing soil shows high sensitivity in runoff prediction, whereas LID isolated from soils show less impact. Simulations with LID on higher permeability soils indicate that nearly all of pre-LID runoff is treated; while features interacting with less-permeable soils treat only 50%. We highlight the NSWC as a screening-level tool for site runoff dynamics and its suitability in stormwater management.
Benthic algal (periphyton) growth rates in response to nitrogen and phosphorus: Parameter estimation for water quality models
Nitrogen (N) and phosphorus (P) are significant pollutants that can stimulate nuisance blooms of algae. Water quality models (e.g., WASP, CE-QUAL-R1, CE-QUAL-ICM, QUAL2k) are valuable and widely used management tools for algal accrual due to excess nutrients in the presence of other limiting factors. These models utilize the Monod and Droop equations to associate algal growth rate with dissolved nutrient concentration and intra-cellular nutrient content. Having accurate parameter values is essential to model performance, however published values for model parameterization are limited, particularly for benthic (periphyton) algae. We conducted a 10-day mesocosm experiment and measured diatom-dominated periphyton biomass accrual through time as chlorophyll (chl ) and ash-free dry mass (AFDM) in response to additions of N (range 5-11,995 μg NO-N/L) and P (range 0.89-59.51 μg SRP/L). Resulting half saturation coefficients and growth rates are similar to other published values, but minimum nutrient quotas are higher than those previously reported. Saturation concentration for N ranged from 150 to 2450 μg NO-N/L based on chl and from 8.5 to 60 μg NO-N/L when based on AFDM. Similarly, the saturation concentration for P ranged from 12 to 29 μg-P/L based on chl , and from 2.5 to 6.1 μg-P/L based on AFDM. These saturation concentrations provide an upper limit for streams where diatom growth can be expected to respond to nutrient levels and a benchmark for reducing nutrient concentrations to a point where benthic algal growth will be limited.