Fishes of Indiana streams: current and historic assemblage structure

We examined Indiana fish assemblages using taxonomy and ecological categories to assess temporal shifts in community structure and recent environmental relationships. Historic (1945) and recent (1996–2007) presence/absence data were compiled by subbasin and analyzed with Nonmetric Multidimensional Scaling (NMS) ordination and by species richness. Canonical Correspondence Analysis (CCA) was used to test taxonomic identity and ecological category abundance data for explanation with recent (1996–2007) environmental variables. We found a decrease in assemblage heterogeneity for recent assemblages and an increase in the number of tolerant species per subbasin. Recent Indiana streams are dominated by tolerant fishes with generalist life history strategies and low functional variation. The use of ecological categories resulted in weaker relationships with environmental variables than analyses with taxonomic identities. Analyses using taxonomy resulted in strong assemblage explanation from stream size and flow variation, while analyses using ecological categories resulted in strong assemblage explanation from habitat variation in silt substrates and flow. Analyses of recent assemblage structure using ecological categories resulted in decreased assemblage variation among subbasins than in analyses using taxonomic identities. We found that fish assemblages of Indiana streams are structured primarily by habitat complexity and have been altered during the past 50 years through multiple disturbances including fragmentation, siltation, and species introductions.     

Influences of low-head dams on the fish assemblages in the headwater streams of the Qingyi watershed, China

The influences of low-head dams on the fish assemblages were examined in this study, using fish data collected in six treatment and five reference sites at three low-head dams in the headwater streams of the Qingyi watershed, China. Comparing with those in the reference sites, local habitat variables were significantly altered by low-head dams in the treatment sites, involving wider channel (only in the impoundment area), deeper water and slower flow. Fish species richness varied significantly across seasons, not across site categories, suggesting that these low-head dams did not alter species richness. However, significant decreases in fish abundance and density were observed in the impoundment areas immediately upstream of dams, but not in the plunge areas downstream. Fish assemblage structures kept relative stability across seasons, and their significant difference between-site was only observed between the impoundment areas and the sites far from dams upstream. This variation in assemblage structures was due to the differing relative abundance of some co-occurring species; more lentic but less lotic fish was observed in the impoundment areas. The spatial and temporal patterns of fish assemblages were correlated with local habitat in this study area. Wetted width had negative correlation with fish species richness, abundance and density, respectively. Water temperature also positively affected species richness. In addition, wetted width, water depth, current velocity and substrate were the important habitat variables influencing assemblage structures. Our results suggested that, by modifying local habitat characteristics, low-head dams altered fish abundance and density in the impoundment areas immediately upstream of dam, not in the plunge areas immediately downstream, and thereby influenced fish assemblage structures in these stream segments.     

Remotely-sensed indicators for monitoring the general condition of “natural habitat” in watersheds: an application for Delaware’s Nanticoke River watershed

Over the past two decades there has been increasing interest in developing indicators to monitor environmental change. Remote sensing techniques have been primarily used to generate information on land use/land cover changes. The US Fish and Wildlife Service has used this technology to monitor wetland trends and recently developed a set of remotely-sensed indicators to characterize and assess trends in the integrity of natural habitat in watersheds. The indices largely focus on the extent of “natural” cover throughout a given watershed, with an emphasis on locations important to fish, wildlife, and water quality. Six indices address natural habitat extent and four deal with human-caused disturbance. A composite index of natural habitat integrity combining the habitat extent and habitat disturbance indices may be formulated to provide an overall numeric value for a watershed or subbasin. These indices facilitate comparison between watersheds (and subbasins) and assesssment of trends useful for environmental monitoring. This paper describes the indices and presents an example of their application for characterizing and assessing conditions of subbasins within Delaware’s Nanticoke River watershed.     

Hierarchical large-scale to local-scale influence of abiotic factors in summer-fragmented Mediterranean rivers: structuring effects on fish distributions,assemblage composition and species richness

Summer rainfall shortage and low or no flows in Mediterranean rivers make it difficult for fish to attain their final preferenda and hinder the assembling of fish communities. We investigated in this study the degree to which large-scale features of the watersheds and reach-scale habitats correlated with a non-random structure of fish species' richness, composition and distribution under these conditions by intensively sampling a river system during summer and using Geographical Information Systems data. The results of multivariate ordination and association analyses showed that multi-scale factors were hierarchically related and that they correlated with non-random fish distributions and assemblage structure. They also highlighted large-scale factors shaping reach-scale water availability for fish during summer. Fish species' richness, composition and distribution varied along a gradient of water shortage and connectivity loss which were more severe with the smaller drainage area into the reaches upstream. Fish species' distributions along these gradients were strongly correlated with both life-history traits (adult size) and biogeographical origin (native vs. exotic). The results for the summer environmental segregation of exotic versus small endemic richness characteristics indicate strong species' interactions forcing the latter into the upstream, harsher seasonal habitats. Special fitness costs are to be expected in drought years, in regulated rivers, and in general, in a climate change scenario in Mediterranean river networks. The study emphasizes that the conservation of the threatened Mediterranean freshwater fish faunas requires operating at the proper spatial scale. In particular, reach-scale habitat improvements are adequate but not enough during the summer in fragmenting river networks, unless there is proper watershed-scale management of the flow regimes.     

Fish assemblages of the Casiquiare River, a corridor and zoogeographical filter for dispersal between the Orinoco and Amazon basins

Aim The aim of this study was to determine whether the Casiquiare River functions as a free dispersal corridor or as a partial barrier (i.e. filter) for the interchange of fish species of the Orinoco and Negro/Amazon basins using species assemblage patterns according to geographical location and environmental features. Location The Casiquiare, Upper Orinoco and Upper Negro rivers in southern Venezuela, South America. Methods Our study was based on an analysis of species presence/absence data and environmental information (11 habitat characteristics) collected by the authors and colleagues between the years 1984 and 1999. The data set consisted of 269 sampled sites and 452 fish species (> 50,000 specimens). A wide range of habitat types was included in the samples, and the collection sites were located at various points along the entire length of the Casiquiare main channel, at multiple sites on its tributary streams, as well as at various nearby sites outside the Casiquiare drainage, within the Upper Orinoco and Upper Rio Negro river systems. Most specimens and field data used in this analysis are archived in the Museo de Ciencias Naturales in Guanare, Venezuela. We performed canonical correspondence analysis (CCA) based on species presence/absence using two versions of the data set: one that eliminated sites having < 5 species and species occurring at < 5 sites; and another that eliminated sites having < 10 species and species occurring at < 10 sites. Cluster analysis was performed on sites based on species assemblage similarity, and a separate analysis was performed on species based on CCA loadings. Results The CCA results for the two versions of the data set were qualitatively the same. The dominant environmental axis contrasted assemblages and sites associated with blackwater vs. clearwater conditions. Longitudinal position on the Casiquiare River was correlated (r² = 0.33) with CCA axis‐1 scores, reflecting clearwater conditions nearer to its origin (bifurcation of the Orinoco) and blackwater conditions nearer to its mouth (junction with the Rio Negro). The second CCA axis was most strongly associated with habitat size and structural complexity. Species associations derived from the unweighted pair‐group average clustering method and pair‐wise squared Euclidean distances calculated from species loadings on CCA axes 1 and 2 showed seven ecological groupings. Cluster analysis of species assemblages according to watershed revealed a stronger influence of local environmental conditions than of geographical proximity. Main conclusions Fish assemblage composition is more consistently associated with local environmental conditions than with geographical position within the river drainages. Nonetheless, the results support the hypothesis that the mainstem Casiquiare represents a hydrochemical gradient between clearwaters at its origin and blackwaters at its mouth, and as such appears to function as a semi‐permeable barrier (environmental filter) to dispersal and faunal exchanges between the partially vicariant fish faunas of the Upper Orinoco and Upper Negro rivers.     

Genetic diversity and species diversity of stream fishes covary across a land-use gradient

Genetic diversity and species diversity are expected to covary according to area and isolation, but may not always covary with environmental heterogeneity. In this study, we examined how patterns of genetic and species diversity in stream fishes correspond to local and regional environmental conditions. To do so, we compared population size, genetic diversity and divergence in central stonerollers (Campostoma anomalum) to measures of species diversity and turnover in stream fish assemblages among similarly sized watersheds across an agriculture–forest land-use gradient in the Little Miami River basin (Ohio, USA). Significant correlations were found in many, but not all, pair-wise comparisons. Allelic richness and species richness were strongly correlated, for example, but diversity measures based on allele frequencies and assemblage structure were not. In-stream conditions related to agricultural land use were identified as significant predictors of genetic diversity and species diversity. Comparisons to population size indicate, however, that genetic diversity and species diversity are not necessarily independent and that variation also corresponds to watershed location and glaciation history in the drainage basin. Our findings demonstrate that genetic diversity and species diversity can covary in stream fish assemblages, and illustrate the potential importance of scaling observations to capture responses to hierarchical environmental variation. More comparisons according to life history variation could further improve understanding of conditions that give rise to parallel variation in genetic diversity and species diversity, which in turn could improve diagnosis of anthropogenic influences on aquatic ecosystems.     

Spatial and temporal variation of fish assemblages and their associations to habitat variables in a mountain stream of north Tiaoxi River,China

The spatial and temporal variations of the fish assemblages in mountain streams of China are poorly understood. The relationships between the fish assemblage and selected habitat features were examined in the North Tiaoxi River, one of headwaters of Taihu Lake. A total of 3,348 individuals belonging to 5 orders, 11 families, 25 genera and 34 species were collected including 33 native species and one invasive species. Among those, about 20 species were endemic to China. Non-metric Multidimensional Scaling (NMDS) was applied to compare fish assemblage structures from upstream to downstream during four seasons. Species assemblages differed along the stream continuum, but there was little apparent change associated with the seasons. Species richness and Shannon-Weaver index (H′) tended to increase along the stream continuum from the upstream to downstream and the proportion of invertivorous fish tended to significantly decrease along the continuum with a parallel significant increase in the percentage of omnivores. Fish assemblages were significantly related to both water quality and habitat structure variables. Canonical Correspondence Analysis ordinations (CCA) revealed that 6 of the 14 selected environmental variables had significant relationships with the fish assemblage such as distance to source, stream width, altitude, pH, water depth, and water velocity and different sampling sites were associated with different environmental variables in different seasons. The main differences in fish assemblage structure and diversity within the whole watercourse are probably related to large-scale factors such distance to source, altitude and stream width. Differences of instream characteristics are likely to be caused by natural variability of the ecosystems but also, in some case, by anthropogenic influence like human settlements, agriculture and river embankment and pollution from small factory.     

Influence of the macrophyte <Emphasis Type="Italic">Eichhornia azurea</Emphasis> on fish assemblage of the Upper Paraná River floodplain (Brazil)

The architecture of aquatic macrophytes adds structural complexity to the littoral region. This increased habitat heterogeneity prompts greater diversity and stability of biotic communities. However, there are few studies that explore the ecological role of macrophytes in the Neotropical region. Then, the aim of this study was to determine spatial and diel influence of near-shore Eichhornia azurea stands on fish assemblage attributes and structure, in the Upper Paraná River floodplain, considering a meso-habitat scale. To achieve this objective, fish were trapped in four stands of macrophytes, in three positions (middle and border of the stand and open area) with samplings taken at different times of the day (6 h, 12 h and 18 h). A total of 537 fish of 16 species were caught. Dissolved oxygen varied significantly among positions. Fish assemblage attributes (species richness, evenness and Shannon Diversity Index) also varied among positions. Fish assemblage structure (summarized by a detrended correspondence analysis) and size structure were strongly influenced by times of the day and positions. The patterns observed can be explained by factors such as oxygen concentration, food resource availability and predation, all influenced by macrophytes. We suggest that any strategy to manage fish diversity and fish stocks has to consider aquatic macrophytes.     

Spatial Variation in Fish Assemblages across small Mediterranean Drainages: Effects of Habitat and Landscape Context

Much debate about assemblage organization in stream fish may stem from analysing the effects of both local and large-scale processes on assemblage attributes over whole geographic regions. This study addresses this issue, by examining the contribution of local habitat attributes and landscape context to fish assemblage variation across small Mediterranean drainages in southern Portugal. Fish abundance and species composition was estimated in 28 sites, across 10 drainages, in both a dry year (1999) and a wet year (2001), and related to two sets of variables reflecting habitat and landscape characteristics. Fish showed responses to both sets of variables with variance partitioning indicating that landscape context had important effects on species richness whereas habitat attributes were the primary determinants of local fish abundance. In general, high species richness was associated with larger drainage area and higher rainfall variability, whereas variation in species abundances mostly reflected the influence of width, depth, conductivity, current velocity, substrate size and emergent vegetation. The relative contributions of both landscape context and habitat attributes to species richness and abundance were generally lower in 1999 than in 2001, with much less diversified species–habitat relationships being found in the former dry year. These results point to the dynamic nature of assemblage organization, emphasizing the importance of innovative, multi-scale approaches in advancing our knowledge of fish assemblage structure in Mediterranean streams.     

Dynamic modeling of nitrogen losses in river networks unravels the coupled effects of hydrological and biogeochemical processes

The importance of lotic systems as sinks for nitrogen inputs is well recognized. A fraction of nitrogen in streamflow is removed to the atmosphere via denitrification with the remainder exported in streamflow as nitrogen loads. At the watershed scale, there is a keen interest in understanding the factors that control the fate of nitrogen throughout the stream channel network, with particular attention to the processes that deliver large nitrogen loads to sensitive coastal ecosystems. We use a dynamic stream transport model to assess biogeochemical (nitrate loadings, concentration, temperature) and hydrological (discharge, depth, velocity) effects on reach-scale denitrification and nitrate removal in the river networks of two watersheds having widely differing levels of nitrate enrichment but nearly identical discharges. Stream denitrification is estimated by regression as a nonlinear function of nitrate concentration, streamflow, and temperature, using more than 300 published measurements from a variety of US streams. These relations are used in the stream transport model to characterize nitrate dynamics related to denitrification at a monthly time scale in the stream reaches of the two watersheds. Results indicate that the nitrate removal efficiency of streams, as measured by the percentage of the stream nitrate flux removed via denitrification per unit length of channel, is appreciably reduced during months with high discharge and nitrate flux and increases during months of low-discharge and flux. Biogeochemical factors, including land use, nitrate inputs, and stream concentrations, are a major control on reach-scale denitrification, evidenced by the disproportionately lower nitrate removal efficiency in streams of the highly nitrate-enriched watershed as compared with that in similarly sized streams in the less nitrate-enriched watershed. Sensitivity analyses reveal that these important biogeochemical factors and physical hydrological factors contribute nearly equally to seasonal and stream-size related variations in the percentage of the stream nitrate flux removed in each watershed.     

Fish assemblage of a Pampasic stream (Buenos Aires,Argentina): temporal variations and relationships with environmental variables

In the Argentine Pampa fertile soils were originally covered by grasslands, but at present are intensively cultivated. We assessed the specific composition of the fish assemblage of El Pescado stream and compared it with that recorded in 1991–1993, when land use in the watershed consisted in natural pastures. The persistence of the fish assemblage between the two studies was rather high: 0.76. Abundance, biomass and species richness were higher during drought periods. Connectivity with the huge Río de la Plata hydrographic system seems the most important contribution to the high and stable species richness of El Pescado stream.     

Fish-habitat relationship in a tropical river under anthropogenic influences

This study analyzes the interaction of fish assemblages with 14 physicochemical and hydrogeomorphological variables at 31 sampling stations along the watershed of the Meia Ponte River, Upper Paraná Basin, Central Brazil, during low and high water seasons in 2001. This river and its tributaries drain both urban and agricultural areas. Fish were caught with sieves along a 100 m stretch demarcated in every sampling site, where environmental variables were also measured. A total of 3508 individuals belonging to 31 species were collected. Fish abundance and environmental data matrices were submitted to a multivariate analysis of co-inertia. Two axes were retained for interpretation (total variance explained = 63.65%) indicating that pH, water temperature, conductivity, chemical dissolved oxygen, and turbidity, all have an influence on fish assemblage structure. The co-structure found (fish assemblages and physicochemical variables) is correlated in both of the axes considered (r = 0.73 and r = 0.68, respectively), and is statistically significant (Monte Carlo test, P < 0.001). This co-structure is regulated by seasonality, but is influenced by fish habitat preferences, spawning and available food, the extent and effects of anthropogenic activities (domestic sewage, agriculture, ranching, urban areas) and the position of sampling stations along the watershed. Handling editor: S. M. Thomaz     

Distribution and abundance of fish eggs and larvae in three tributaries of the Upper Uruguay River (Brazil)

The present study aimed to investigate the distribution and abundance of fish eggs and larvae in three important tributaries (Chapecó, Peixe, and Ligeiro rivers) of the Upper Uruguay River. The spatial and temporal distribution of fish eggs and larvae were studied as well as the correlation between environmental parameters and the abundance of ichthyoplankton species. The study was conducted between October 2005 and September 2006. Ichthyoplankton samples were collected at night with cylindroconical 0.5-mm mesh plankton nets every 5 days. Of the 591 samples collected, 170 contained ichthyoplankton organisms, resulting in the capture of 12,847 fish eggs and 962 fish larvae. Twenty-seven fish species were observed, of which 69% were Characiforms and 27% were Siluriforms. Among the fish species captured, the representatives were predominantly young forms of small and medium size fishes, with rheophilic species occurring infrequently. Eggs occurred exclusively between October and January, while the highest larvae occurrence was observed between November and December. Fish larvae assemblage structure was shown to be related to some environmental variables. There was a tendency of higher values of water temperature and velocity at the lower sampling sites than at the upper ones, as well an increase number of eggs and larvae. The study tributaries serve as reproduction sites and nursery areas for several fish species of the Upper Uruguay River; thus, the maintenance of their integrity is important for the preservation of diversity and enhancement of fisheries in the region.     

Assignment and evaluation of ground beetle (Coleoptera: Carabidae) assemblages to sites on different scales in a grassland landscape

Carabid beetles and environmental parameters were investigated in 52 grassland sites with three replicate pitfall traps in each site and in the valley of the River Eider in Schleswig–Holstein (northern Germany) with 61 pitfall traps. Environmental parameters included water content of soil, sand content, organic matter content and pH. Ground beetle assemblages were derived by detrended correspondence analyses (DCA) and characterised by the specific environmental conditions as means for each assemblage. On the regional scale including all investigated sites of Schleswig–Holstein, five assemblages were differentiated. On the local scale including the investigated sites in the valley of the River Eider, three assemblages were found corresponding well with those found on the regional scale. Environmental conditions at the sites of the five assemblages were correlated with land use data, soil types, and water level stages provided by three maps of a geographic information system (GIS). The GIS maps were combined to develop smaller areas with land use, soil type and water level stage information. The characteristic environmental conditions were assigned to each area to derive the spatial distribution of the five ground beetle assemblages. Spatial prediction was correct for 65% of investigated sites. The potential area of each assemblage was estimated for the valley. The different grassland areas were evaluated as potential habitats for ground beetle species comparing total species richness with the regional species richness of each assemblage. The comparison shows that species richness in the evaluated assemblages is relatively low compared to the regional potential.     

Multi-scale assessment of habitats and stressors influencing stream fish assemblages in the Lake Pontchartrain Basin,USA

This study determined if a spatial hierarchy existed with respect to reach-scale habitat, subwatershed-scale, and watershed-scale geomorphology and land use stressors on fish assemblages in southeastern U.S. coastal plain streams. During May–October 2009–2012, fishes were sampled by seine, and habitat was assessed at 50 reaches in the Lake Pontchartrain Basin (USA). Using partial redundancy analysis (pRDA) a variance decomposition procedure was used to partial out influences of confounding covariables at each spatial scale. Reach-scale habitat had the strongest association with the assemblage. Stream width, depth, aquatic vegetation and human debris cover, rapid habitat assessment score, and large woody debris volume were the most important variables. At subwatershed and watershed scales, natural and anthropogenic characteristics were important, including elevation, gradient, watershed area, wetland cover, stream density, road, dam and oil/gas well densities. Six species were associated most strongly with the watershed variables, compared to reach- and subwatershed-scale variables. These species had more “r-selected” life-history strategies (e.g., smaller eggs, shorter life spans, multiple broods, longer spawning season, and trophic generalists). In contrast, most species that were associated strongest with reach-scale variables exhibited more “k-selected” life-history traits (e.g., larger eggs, longer life spans, shorter spawning season, single brood, and trophic specialists).     

Patterns of fish species richness in the Seine River basin,France

Variation in fish assemblage structure was examined in different rivers of the Seine Basin. Factor analysis was used to identify similarities among samples and species, and to show ichtyological changes along an upstream-downstream gradient. Fish species richness was correlated with catchment area, distance from the source and number of individuals. Relationships between species richness and number of individuals in a sample were largely an artifact of sampling.Species richness increased with river size, reached a maximum in midsize rivers, then decreased in large rivers. This species richness pattern was consistent with the model of the River Continuum Concept. However, possible causes could include greater pollution effects in large rivers. Distribution of feeding guilds was related to river size. Species richness and proportion of omnivores and piscivores increased with river size, whereas species richness and proportion of invertivores declined downstream. These patterns suggest that trophic diversity of fish assemblages may be related to food availability.     

Spatial and Temporal Patterns of Diversity and Distribution of the Upper Juruá River Fish Community (Brazilian Amazon)

We studied spatial and temporal patterns in fish species composition and diversity at the upper Juruá River located in the west Brazilian Amazon. We collected with gillnet 822 fishes belonging to 90 species in the main Juruá River, its tributaries and the floodplain lakes during wet and dry seasons. Fish abundance and species richness were greater in the dry season. During that season, fishes may be concentrated due to the low water level, being caught more easily by gillnets. There has been a trend towards a greater fish biomass caught in lakes. This might be associated with a greater environmental stability as lakes may be less subject to large variations in water level. The fish communities differed between the two seasons and between lakes and the lotic environments (main river and tributaries). Fish species from the family Curimatidae were most abundant in the lakes, while Pimelodus spp. and Hypostomus spp. predominated in the main Juruá River. Seasonal variations in fish communities may be related to differences in the migratory behavior among fish species. Such spatial and temporal patterns influencing fish community structure at the Upper Juruá Extractive Reserve must be accounted for in management and conservation strategies.     

Altitudinal gradients in stream fish diversity and the prevalence of diadromy in the Sixaola River basin,Costa Rica

Landscape-scale patterns of freshwater fish diversity and assemblage structure remain poorly documented in many areas of Central America, while aquatic ecosystems throughout the region have been impacted by habitat degradation and hydrologic alterations. Diadromous fishes may be especially vulnerable to these changes, but there is currently very little information available regarding their distribution and abundance in Central American river systems. We sampled small streams at 20 sites in the Sixaola River basin in southeastern Costa Rica to examine altitudinal variation in the diversity and species composition of stream fish assemblages, with a particular focus on diadromous species. A set of environmental variables was also measured in the study sites to evaluate how changes in fish assemblage structure were related to gradients in stream habitat. Overall, fish diversity and abundance declined steeply with increasing elevation, with very limited species turnover. The contribution of diadromous fishes to local species richness and abundance increased significantly with elevation, and diadromous species dominated assemblages at the highest elevation sites. Ordination of the sampling sites based on fish species composition generally arranged sites by elevation, but also showed some clustering based on geographic proximity. The dominant gradient in fish community structure was strongly correlated with an altitudinal habitat gradient identified through ordination of the environmental variables. The variation we observed in stream fish assemblages over relatively small spatial scales has significant conservation implications and highlights the ecological importance of longitudinal connectivity in Central American river systems.     

Impacts of urbanization on stream habitats and macroinvertebrate communities in the tributaries of Qiangtang River,China

The impacts of watershed urbanization on streams have been studied worldwide, but are rare in China. We examined relationships among watershed land uses and stream physicochemical and biological attributes, impacts of urbanization on overall stream conditions, and the response pattern of macroinvertebrate assemblage metrics to the percent of impervious area (PIA) of watersheds in the middle section of the Qiantang River, Zhejiang Province, China. Environmental variables and benthic macroinvertebrates of 60 stream sites with varied levels of watershed urban land use were sampled in April, 2010. Spearman correlation analysis showed watershed urbanization levels significantly correlated with increased stream depth, width, and values of conductivity, total nitrogen, ammonia, phosphate, calcium, magnesium, and chemical oxygen demand for the study streams. There was significant difference in total taxa richness, Empheroptera, Plecoptera, and Trichoptera (EPT) taxa richness, and Diptera taxa richness, percentages of individual abundances of EPT, Chironomidae, shredders, filterers, and scrapers, and Shannon–Wiener diversity index between reference streams and urban impacted streams. In contrast, percentages of individual abundances for collectors, oligochaeta, and tolerant taxa, and biotic index were significantly higher in urban impacted than reference streams. All the above metrics were significantly correlated with PIA. The response patterns of total taxa richness, EPT taxa richness, and Shannon–Wiener diversity index followed a drastic decrease at thresholds of 3.6, 3.7, and 5.5% of PIA, respectively. Our findings indicate that stream benthic macroinvertebrate metrics are effective indicators of impacts of watershed urban development, and the PIA-imperviousness thresholds we identified could potentially be used for setting benchmarks for watershed development planning and for prioritizing high valued stream systems for protection and rehabilitation.