SPARROW Modeling of Nitrogen Sources and Transport in Rivers and Streams of California and Adjacent States,U.S.

The SPARROW (SPAtially Referenced Regressions On Watershed attributes) model was used to evaluate the spatial distribution of total nitrogen (TN) sources, loads, watershed yields, and factors affecting transport and decay in the stream network of California and portions of adjacent states for the year 2002. The two major TN sources to local catchments on a mass basis were fertilizers and manure (51.7%) and wastewater discharge (15.9%). Other sources contributed < 12%. Fertilizer use is widespread in the Central Valley region of California, and also important in several other regions because of the diversity of California agriculture. Precipitation, sand content of surficial soils, wetlands, and tile drains were important for TN movement to stream reaches. Median streamflow in the study area is about 0.04 m³/s. Aquatic losses of nitrogen were found to be most important in intermittent and small to medium sized streams (0.2‐14 m³/s), while larger streams showed less loss, and therefore are important for TN transport. Nitrogen loss in reservoirs was found to be insignificant, possibly because most of the larger ones are located upstream of nitrogen sources. The model was used to show loadings, sources, and tributary inputs to several major rivers. The information provided by the SPARROW model is useful for determining both the major sources contributing nitrogen to streams and the specific tributaries that transport the load.     

Stream Discharge and Riparian Land Use Influence In-Stream Concentrations and Loads of Phosphorus from Central Plains Watersheds

Total annual nutrient loads are a function of both watershed characteristics and the magnitude of nutrient mobilizing events. We investigated linkages among land cover, discharge and total phosphorus (TP) concentrations, and loads in 25 Kansas streams. Stream monitoring locations were selected from the Kansas Department of Health and Environment stream chemistry long-term monitoring network sites at or near U.S. Geological Survey stream gauges. We linked each sample with concurrent discharge data to improve our ability to estimate TP concentrations and loads across the full range of possible flow conditions. Median TP concentration was strongly linked (R ² = 76%) to the presence of cropland in the riparian zones of the mostly perennial streams. At baseflow, discharge data did not improve prediction of TP, but at high flows discharge was strongly linked to concentration (a threshold response occurred). Our data suggest that on average 88% of the total load occurred during the 10% of the time with the greatest discharge. Modeled reductions in peak discharges, representing increased hydrologic retention, predicted greater decreases in total annual loads than reductions of ambient concentrations because high discharge and elevated phosphorus concentrations had multiplicative effects. No measure of land use provided significant predictive power for concentrations when discharge was elevated or for concentration rise rates under increasing discharge. These results suggest that reductions of baseflow concentrations of TP in streams without wastewater dischargers may be managed by reductions of cropland uses in the riparian corridor. Additional measures may be needed to manage TP annual loads, due to the large percentage of the TP load occurring during a few high-flow events each year.     

A Regional Modeling Framework of Phosphorus Sources and Transport in Streams of the Southeastern United States1

García, Ana María, Anne B. Hoos, and Silvia Terziotti, 2011. A Regional Modeling Framework of Phosphorus Sources and Transport in Streams of the Southeastern United States. Journal of the American Water Resources Association (JAWRA) 47(5):991‐1010. DOI: 10.1111/j.1752‐1688.2010.00517.x Abstract: We applied the SPARROW model to estimate phosphorus transport from catchments to stream reaches and subsequent delivery to major receiving water bodies in the Southeastern United States (U.S.). We show that six source variables and five land‐to‐water transport variables are significant (p < 0.05) in explaining 67% of the variability in long‐term log‐transformed mean annual phosphorus yields. Three land‐to‐water variables are a subset of landscape characteristics that have been used as transport factors in phosphorus indices developed by state agencies and are identified through experimental research as influencing land‐to‐water phosphorus transport at field and plot scales. Two land‐to‐water variables – soil organic matter and soil pH – are associated with phosphorus sorption, a significant finding given that most state‐developed phosphorus indices do not explicitly contain variables for sorption processes. Our findings for Southeastern U.S. streams emphasize the importance of accounting for phosphorus present in the soil profile to predict attainable instream water quality. Regional estimates of phosphorus associated with soil‐parent rock were highly significant in explaining instream phosphorus yield variability. Model predictions associate 31% of phosphorus delivered to receiving water bodies to geology and the highest total phosphorus yields in the Southeast were catchments with already high background levels that have been impacted by human activity.     

Water Quality Signals from Rural Land Use and Exurbanization in a Mountain Landscape: What's Clear and What's Confounded?

In mountainous landscapes with high climatic and geomorphic variability, how do rural land uses and exurbanization alter hydrology and water quality? We evaluated effects of rural land use and exurbanization on streamflows, suspended sediment concentrations and loads, specific conductance, and summer water temperatures in 12 streams and rivers within the Upper Little Tennessee River basin in the southern Appalachian Mountains. Eleven streams featured low levels of development (>61% forest cover) but differed in land use patterning, basin size, annual precipitation, and watershed morphology. One urban stream, located within the largest town in the basin, provided the high development comparative endpoint. Even low levels of rural development and exurbanization were associated with substantial increases in suspended sediment concentrations, sediment loads, and summer stream temperature daily maxima and diurnal variation. Observed summer temperature increases were much larger than would be expected due to global climate change over the next century. Specific conductance was idiosyncratic among the smaller streams. These water quality changes were not accompanied by streamflow changes that were discernible amid the high natural variation in precipitation and geomorphology. The water quality findings suggest the need for applying the best management practices, including riparian buffers, to even low levels of rural development.     

Phosphorus, sediment, and Escherichia coli loads in unfenced streams of the Georgia Piedmont, USA

Contamination of unfenced streams with P, sediments, and pathogenic bacteria from cattle (Bos taurus) activity may be affected by the availability of shade and alternative water sources. The objectives of this study were to evaluate water quality in two streams draining tall fescue (Festuca arundinacea Schreb.)-common bermudagrass (Cynodon dactylon L.) pastures with different shade distribution, and to quantify the effects of alternative water sources on stream water quality. For 3 yr, loads of dissolved reactive phosphorus (DRP), total phosphorus (TP), and total suspended solids (TSS) were measured during storm flow, and loads of DRP, TP, TSS, and Escherichia coli were measured every 14 d during base flow. We also used GPS collars to determine amount of time cattle spent in riparian areas. Our results showed that cattle-grazed pastures with unfenced streams contributed significant loads of DRP, TP, TSS, and E. coli to surface waters (p < 0.01). Time spent by cattle in riparian areas as well as storm flow loads of DRP, TP, and TSS were larger (p < 0.08) in the pasture with the smaller amount of nonriparian shade. Water trough availability decreased base flow loads of TSS and E. coli in both streams, and decreased time cattle spent in riparian areas in the pasture with the smaller amount of nonriparian shade (p < 0.08). Our results indicate that possible BMPs to reduce contamination from cattle-grazed pastures would be to develop or encourage nonriparian shade and to provide cattle with alternative water sources away from the stream.     

Landscape Planning for Agricultural Nonpoint Source Pollution Reduction III: Assessing Phosphorus and Sediment Reduction Potential

Riparian buffers have the potential to improve stream water quality in agricultural landscapes. This potential may vary in response to landscape characteristics such as soils, topography, land use, and human activities, including legacies of historical land management. We built a predictive model to estimate the sediment and phosphorus load reduction that should be achievable following the implementation of riparian buffers; then we estimated load reduction potential for a set of 1598 watersheds (average 54 km²) in Wisconsin. Our results indicate that land cover is generally the most important driver of constituent loads in Wisconsin streams, but its influence varies among pollutants and according to the scale at which it is measured. Physiographic (drainage density) variation also influenced sediment and phosphorus loads. The effect of historical land use on present-day channel erosion and variation in soil texture are the most important sources of phosphorus and sediment that riparian buffers cannot attenuate. However, in most watersheds, a large proportion (approximately 70%) of these pollutants can be eliminated from streams with buffers. Cumulative frequency distributions of load reduction potential indicate that targeting pollution reduction in the highest 10% of Wisconsin watersheds would reduce total phosphorus and sediment loads in the entire state by approximately 20%. These results support our approach of geographically targeting nonpoint source pollution reduction at multiple scales, including the watershed scale.     

Temporal and Spatial Trends in Nutrient and Sediment Loading to Lake Tahoe,California‐Nevada,USA

Since 1980, the Lake Tahoe Interagency Monitoring Program (LTIMP) has provided stream‐discharge and water quality data—nitrogen (N), phosphorus (P), and suspended sediment—at more than 20 stations in Lake Tahoe Basin streams. To characterize the temporal and spatial patterns in nutrient and sediment loading to the lake, and improve the usefulness of the program and the existing database, we have (1) identified and corrected for sources of bias in the water quality database; (2) generated synthetic datasets for sediments and nutrients, and resampled to compare the accuracy and precision of different load calculation models; (3) using the best models, recalculated total annual loads over the period of record; (4) regressed total loads against total annual and annual maximum daily discharge, and tested for time trends in the residuals; (5) compared loads for different forms of N and P; and (6) tested constituent loads against land use‐land cover (LULC) variables using multiple regression. The results show (1) N and P loads are dominated by organic N and particulate P; (2) there are significant long‐term downward trends in some constituent loads of some streams; and (3) anthropogenic impervious surface is the most important LULC variable influencing water quality in basin streams. Many of our recommendations for changes in water quality monitoring and load calculation methods have been adopted by the LTIMP.     

THE IMPACT OF RUNOFF GENERATION MECHANISMS ON THE LOCATION OF CRITICAL SOURCE AREAS1

ABSTRACT: Identifying phosphorus (P) source areas and transport pathways is a key step in decreasing P loading to natural water systems. This study compared the effects of two modeled runoff generation processes ‐ saturation excess and infiltration excess ‐ on total phosphorus (TP) and soluble reactive phosphorus (SRP) concentrations in 10 catchment streams of a Catskill mountain watershed in southeastern New York. The spatial distribution of runoff from forested land and agricultural land was generated for both runoff processes; results of both distributions were consistent with Soil Conservation Service‐Curve Number (SCS‐CN) theory. These spatial runoff distributions were then used to simulate stream concentrations of TP and SRP through a simple equation derived from an observed relation between P concentration and land use; empirical results indicate that TP and SRP concentrations increased with increasing percentage of agricultural land. Simulated TP and SRP stream concentrations predicted for the 10 catchments were strongly affected by the assumed runoff mechanism. The modeled TP and SRP concentrations produced by saturation excess distribution averaged 31 percent higher and 42 percent higher, respectively, than those produced by the infiltration excess distribution. Misrepresenting the primary runoff mechanism could not only produce erroneous concentrations, it could fail to correctly locate critical source areas for implementation of best management practices. Thus, identification of the primary runoff mechanism is critical in selection of appropriate models in the mitigation of nonpoint source pollution. Correct representation of runoff processes is also critical in the future development of biogeochemical transport models, especially those that address nutrient fluxes.     

Factors Affecting Stream Nutrient Loads: A Synthesis of Regional SPARROW Model Results for the Continental United States1

Preston, Stephen D., Richard B. Alexander, Gregory E. Schwarz, and Charles G. Crawford, 2011. Factors Affecting Stream Nutrient Loads: A Synthesis of Regional SPARROW Model Results for the Continental United States. Journal of the American Water Resources Association (JAWRA) 47(5):891‐915. DOI: 10.1111/j.1752‐1688.2011.00577.x Abstract: We compared the results of 12 recently calibrated regional SPARROW (SPAtially Referenced Regressions On Watershed attributes) models covering most of the continental United States to evaluate the consistency and regional differences in factors affecting stream nutrient loads. The models – 6 for total nitrogen and 6 for total phosphorus – all provide similar levels of prediction accuracy, but those for major river basins in the eastern half of the country were somewhat more accurate. The models simulate long‐term mean annual stream nutrient loads as a function of a wide range of known sources and climatic (precipitation, temperature), landscape (e.g., soils, geology), and aquatic factors affecting nutrient fate and transport. The results confirm the dominant effects of urban and agricultural sources on stream nutrient loads nationally and regionally, but reveal considerable spatial variability in the specific types of sources that control water quality. These include regional differences in the relative importance of different types of urban (municipal and industrial point vs. diffuse urban runoff) and agriculture (crop cultivation vs. animal waste) sources, as well as the effects of atmospheric deposition, mining, and background (e.g., soil phosphorus) sources on stream nutrients. Overall, we found that the SPARROW model results provide a consistent set of information for identifying the major sources and environmental factors affecting nutrient fate and transport in United States watersheds at regional and subregional scales.     

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

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

Dry Weather Water Quality Loadings in Arid,Urban Watersheds of the Los Angeles Basin,California, USA1

Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.     

DETERMINING NUTRIENT EXPORT COEFFICIENTS AND SOURCE LOADING UNCERTAINTY USING IN‐STREAM MONITORING DATA1

ABSTRACT: Over a three‐year period, flow and nutrients were monitored at 13 sites in the upper North Bosque River watershed in Texas. Drainage areas above sampling sites differed in percent of dairy waste application fields, forage fields, wood/range, and urban land area. A multiple regression approach was used to develop total phosphorus (TP) and total nitrogen (TN) export coefficients for the major land uses in these heterogeneous drainage areas. The largest export coefficients were associated with dairy waste application fields followed by urban, forage fields, and wood/range. An empirical model was then established to assess nutrient contribution by major sources using developed export coefficients and point source loadings from municipal wastewater treatment. This model was verified by comparison of estimated loadings to measured in‐stream data. Monte Carlo simulation techniques were applied to provide an uncertainty analysis for nutrient loads by source, based on the variance associated with each export coefficient. The largest sources of nutrients contributing to the upper North Bosque River were associated with dairy waste application fields and forage fields, while the greatest relative uncertainty in source contribution was associated with loadings from urban and wood/range land uses.     

Comparison of Two Parametric Methods to Estimate Pesticide Mass Loads in California’s Central Valley1

Saleh, Dina K., David L. Lorenz, and Joseph L. Domagalski, 2010. Comparison of Two Parametric Methods to Estimate Pesticide Mass Loads in California’s Central Valley. Journal of the American Water Resources Association (JAWRA) 00(0):1‐11. DOI: 10.1111/j.1752‐1688.2010.00506.x Abstract: Mass loadings were calculated for four pesticides in two watersheds with different land uses in the Central Valley, California, by using two parametric models: (1) the Seasonal Wave model (SeaWave), in which a pulse signal is used to describe the annual cycle of pesticide occurrence in a stream, and (2) the Sine Wave model, in which first‐order Fourier series sine and cosine terms are used to simulate seasonal mass loading patterns. The models were applied to data collected during water years 1997 through 2005. The pesticides modeled were carbaryl, diazinon, metolachlor, and molinate. Results from the two models show that the ability to capture seasonal variations in pesticide concentrations was affected by pesticide use patterns and the methods by which pesticides are transported to streams. Estimated seasonal loads compared well with results from previous studies for both models. Loads estimated by the two models did not differ significantly from each other, with the exceptions of carbaryl and molinate during the precipitation season, where loads were affected by application patterns and rainfall. However, in watersheds with variable and intermittent pesticide applications, the SeaWave model is more suitable for use on the basis of its robust capability of describing seasonal variation of pesticide concentrations.     

NUTRIENT CONTRIBUTION OF NONPOINT SOURCE RUNOFF IN THE LAS VEGAS VALLEY1

ABSTRACT: A Geographic Information System (GIS) based non‐point source runoff model is developed for the Las Vegas Valley, Nevada, to estimate the nutrient loads during the years 2000 and 2001. The estimated nonpoint source loads are compared with current wastewater treatment facilities loads to determine the non‐point source contribution of total phosphorus (TP), total nitrogen (TN), and total suspended solids (TSS) on a monthly and annual time scale. An innovative calibration procedure is used to estimate the pollutant concentrations for different land uses based on available water quality data at the outlet. Results indicate that the pollutant concentrations are higher for the Las Vegas Valley than previous published values for semi‐arid and arid regions. The total TP and TN loads from nonpoint sources are approximately 15 percent and 4 percent, respectively, of the total load to the receiving water body, Lake Mead. The TP loads during wet periods approach the permitted loads from the wastewater treatment plants that discharge into Las Vegas Wash. In addition, the GIS model is used to track pollutant loads in the stream channels for one of the subwatersheds. This is useful for planning the location of Best Management Practices to control nonpoint pollutant loads.     

Linking Landscape Characteristics and High Stream Nitrogen in the Oregon Coast Range: Red Alder Complicates Use of Nutrient Criteria

Red alder (Alnus rubra), a nitrogen(N)‐fixing deciduous broadleaf tree, can strongly influence N concentrations in western Oregon and Washington. We compiled a database of stream N and GIS‐derived landscape characteristics in order to examine geographic variation in N across the Oregon Coast Range. Basal area of alder, expressed as a percent of watershed area, accounted for 37% and 38% of the variation in summer nitrate and total N (TN) concentrations, respectively. Relationships between alder and nitrate were strongest in winter when streamflow and landscape connections are highest. Distance to the coast and latitude, potential surrogates for sea salt inputs, and watershed area were also related to nitrate concentrations in an all‐subsets regression analysis, which accounted for 46% of the variation in summer nitrate concentrations. The model with the lowest Akaike's Information Criterion did not include developed or agricultural land cover, probably because few watersheds in our database had substantial levels of these land cover classes. Our results provide evidence, at a regional scale, that background sources and processes cause many Coast Range streams to exceed proposed nutrient criteria, and that the prevalence of a single tree species (N‐fixing red alder) exerts a dominant control over stream N concentrations across this region.     

The Influence of Two‐Stage Ditches with Constructed Floodplains on Water Column Nutrients and Sediments in Agricultural Streams

The two‐stage ditch is a novel management practice originally implemented to increase bank stability through floodplain restoration in channelized agricultural streams. To determine the effects of two‐stage construction on sediment and nutrient loads, we monitored turbidity, and also measured total suspended solids (TSS), dissolved inorganic nitrogen (N) species, and phosphorus (P) after two‐stage ditch construction in reference and manipulated reaches of four streams. Turbidity decreased during floodplain inundation at all sites, but TSS and P, soluble reactive phosphorus (SRP) and total phosphorus (TP) decreased only in the two‐stage ditches with longer duration of inundation. Both TSS and TP were positively correlated within individual streams, but neither were correlated with turbidity. Phosphorus was elevated in the stream to which manure was applied adjacent to the two‐stage reach, but not the reference reach, suggesting that landscape nutrient management plans could restrict nutrient transport to the stream, ultimately determining the efficacy of instream management practices. In addition, ammonium and nitrate decreased in two‐stage reaches with lower initial N concentrations. Overall, results suggest that turbidity, TSS, and TP were reduced during floodplain inundation, but the two‐stage alone may not be effective for managing high inorganic N loads.     

THE EFFECTS OF CLIMATE CHANGE ON STREAM FLOW AND NUTRIENT LOADING'

ABSTRACT: This study assesses the potential impact of climate change on stream flow and nutrient loading in six watersheds of the Susquehanna River Basin using the Generalized Watershed Loading Function (GWLF). The model was used to simulate changes in stream flow and nutrient loads under a transient climate change scenario for each watershed. Under an assumption of no change in land cover and land management, the model was used to predict monthly changes in stream flow and nutrient loads for future climate conditions. Mean annual stream flow and nutrient loads increased for most watersheds, but decreased in one watershed that was intensively cultivated. Nutrient loading slightly decreased in April and late summer for several watersheds as a result of early snowmelt and increasing evapotranspiration. Spatial and temporal variability of stream flow and nutrient loads under the transient climate scenario indicates that different approaches for future water resource management may be useful.     

Influence of Riparian Seepage Zones on Nitrate Variability in Two Agricultural Headwater Streams

Riparian seeps have been recognized for their contributions to stream flow in headwater catchments, but there is limited data on how seeps affect stream water quality. The objective of this study was to examine the effect of seeps on the variability of stream NO₃‐N concentrations in FD36 and RS, two agricultural catchments in Pennsylvania. Stream samples were collected at 10‐m intervals over reaches of 550 (FD36) and 490 m (RS) on 21 occasions between April 2009 and January 2012. Semi‐variogram analysis was used to quantify longitudinal patterns in stream NO₃‐N concentration. Seep water was collected at 14 sites in FD36 and 7 in RS, but the number of flowing seeps depended on antecedent conditions. Seep NO₃‐N concentrations were variable (0.1‐29.5 mg/l) and were often greater downslope of cropped fields compared to other land uses. During base flow, longitudinal variability in stream NO₃‐N concentrations increased as the number of flowing seeps increased. The influence of seeps on the variability of stream NO₃‐N concentrations was less during storm flow compared to the variability of base flow NO₃‐N concentrations. However, 24 h after a storm in FD36, an increase in the number of flowing seeps and decreasing streamflow resulted in the greatest longitudinal variability in stream NO₃‐N concentrations recorded. Results indicate seeps are important areas of NO₃‐N delivery to streams where targeted adoption of mitigation measures may substantially improve stream water quality.     

WATERSHED SCALING EFFECT ON BASE FLOW NITRATE,VALLEY AND RIDGE PHYSIOGRAPHIC PROVINCE1

ABSTRACT: A study of stream base flow and NO₃‐N concentration was conducted simultaneously in 51 subwatersheds within the 116‐square‐kilometer watershed of East Mahantango Creek near Klingerstown, Pennsylvania. The study was designed to test whether measurable results of processes and observations within the smaller watersheds were similar to or transferable to a larger scale. Ancillary data on land use were available for the small and large watersheds. Although the source of land‐use data was different for the small and large watersheds, comparisons showed that the differences in the two land‐use data sources were minimal. A land use‐based water‐quality model developed for the small‐scale 7.3‐square‐kilometer watershed for a previous study accurately predicted NO₃‐N concentrations from sampling in the same watershed. The water‐quality model was modified and, using the imagery‐based land use, was found to accurately predict NO₃‐N concentrations in the subwatersheds of the large‐scale 116‐square‐kilometer watershed as well. Because the model accurately predicts NO₃‐N concentrations at small and large scales, it is likely that in second‐order streams and higher, discharge of water and NO₃‐N is dominated by flow from smaller first‐order streams, and the contribution of ground‐water discharge to higher order streams is minimal at the large scale.     

EFFECTS OF ADJACENT AGRICULTURAL ACTIVITIES AND WATERSHED CHARACTERISTICS ON STREAM MACROINVERTEBRATE COMMUNITIES1

Benthic macroinvertebrate communities in streams adjacent to cornfields, streams where cows had unrestricted access, and reference locations without agriculture were compared to examine the effects of local land use and land use/land cover in the watershed. At each local site, macroinvertebrates and a variety of habitat parameters were measured upstream, adjacent, downstream, and farther downstream of the local land use. A geographic information system (GIS) was used to calculate drainage basin area, land use/land cover percentages in each basin, and the distance from sample sites to the stream source. Three‐way analysis of covariance (ANCOVA) tests with date, site type, and sampling location as main effects were used to explore differences in macroinvertebrate metrics using median substrate size, percent hay/pasture area, and stream depth as covariates. The covariates significantly improved model fit and showed that multiple contributing factors influence community composition. Local impacts were greatest at sites where cows had access, probably because of sedimentation and embeddedness in the substrate. Differences between the upstream and the adjacent and downstream locations were not as great as expected, perhaps because upstream recolonization was reduced by agricultural impacts or because of differences in the intensity or proximity of agriculture to riparian areas in the watershed. The results underscore the importance of both local and watershed factors in controlling stream community composition.