Thomas F. Cuffney

Responses of benthic macroinvertebrates to environmental changes associated with urbanization in nine metropolitan areas

Responses of benthic macroinvertebrates along gradients of urban intensity were investigated in nine metropolitan areas across the United States. Invertebrate assemblages in metropolitan areas where forests or shrublands were being converted to urban land were strongly related to urban intensity. In metropolitan areas where agriculture and grazing lands were being converted to urban land, invertebrate assemblages showed much weaker or nonsignificant relations with urban intensity because sites with low urban intensity were already degraded by agriculture. Ordination scores, the number of EPT taxa, and the mean pollution-tolerance value of organisms at a site were the best indicators of changes in assemblage condition. Diversity indices, functional groups, behavior, and dominance metrics were not good indicators of urbanization. Richness metrics were better indicators of urban effects than were abundance metrics, and qualitative samples collected from multiple habitats gave similar results to those of single habitat quantitative samples (riffles or woody snags) in all metropolitan areas. Changes in urban intensity were strongly correlated with a set of landscape variables that was consistent across all metropolitan areas. In contrast, the instream environmental variables that were strongly correlated with urbanization and invertebrate responses varied among metropolitan areas. The natural environmental setting determined the biological, chemical, and physical instream conditions upon which urbanization acts and dictated the differences in responses to urbanization among metropolitan areas. Threshold analysis showed little evidence for an initial period of resistance to urbanization. Instead, assemblages were degraded at very low levels of urbanization, and response rates were either similar across the gradient or higher at low levels of urbanization. Levels of impervious cover that have been suggested as protective of streams (5-10%) were associated with significant assemblage degradation and were not protective.

On the application of multilevel modeling in environmental and ecological studies

This paper illustrates the advantages of a multilevel/hierarchical approach for predictive modeling, including flexibility of model formulation, explicitly accounting for hierarchical structure in the data, and the ability to predict the outcome of new cases. As a generalization of the classical approach, the multilevel modeling approach explicitly models the hierarchical structure in the data by considering both the within- and between-group variances leading to a partial pooling of data across all levels in the hierarchy. The modeling framework provides means for incorporating variables at different spatiotemporal scales. The examples used in this paper illustrate the iterative process of model fitting and evaluation, a process that can lead to improved understanding of the system being studied.

Responses of Physical, Chemical, and Biological Indicators of Water Quality to a Gradient of Agricultural Land Use in the Yakima River Basin, Washington

The condition of 25 stream sites in the Yakima River Basin, Washington, were assessed by the U.S. Geological Surveys National Water-Quality Assessment Program. Multimetric condition indices were developed and used to rank sites on the basis of physical, chemical, and biological characteristics. These indices showed that sites in the Cascades and Eastern Cascades ecoregions were largely unimpaired. In contrast, all but two sites in the Columbia Basin ecoregion were impaired, some severely. Agriculture (nutrients and pesticides) was the primary factor associated with impairment and all impaired sites were characterized by multiple indicators of impairment. All indices of biological condition (fish, invertebrates, and algae) declined as agricultural intensity increased. The response exhibited by invertebrates and algae suggested a threshold response with conditions declining precipitously at relatively low levels of agricultural intensity and little response at moderate to high levels of agricultural intensity. This pattern of response suggests that the success of mitigation will vary depending upon where on the response curve the mitigation is undertaken. Because the form of the community condition response is critical to effective water-quality management, the National Water-Quality Assessment Program is conducting studies to examine the response of biota to gradients of land-use intensity and the relevance of these responses to water-quality management. These land-use gradient pilot studies will be conducted in several urban areas starting in 1999.

Ambiguous taxa: effects on the characterization and interpretation of invertebrate assemblages

Abstract Damaged and immature specimens often result in macroinvertebrate data that contain ambiguous parent–child pairs (i.e., abundances associated with multiple related levels of the taxonomic hierarchy such as Baetis pluto and the associated ambiguous parent Baetis sp.). The choice of method used to resolve ambiguous parent–child pairs may have a very large effect on the characterization of invertebrate assemblages and the interpretation of responses to environmental change because very large proportions of taxa richness (73–78%) and abundance (79–91%) can be associated with ambiguouparents. To address this issue, we examined 16 variations of 4 basic methods for resolving ambiguous taxa: RPKC (remove parent, keep child), MCWP (merge child with parent), RPMC (remove parent or merge child with parent depending on their abundances), and DPAC (distribute parents among children). The choice of method strongly affected assemblage structure, assemblage characteristics (e.g., metrics), and the ability to detect responses along environmental (urbanization) gradients. All methods except MCWP produced acceptable results when used consistently within a study. However, the assemblage characteristics (e.g., values of assemblage metrics) differed widely depending on the method used, and data should not be combined unless the methods used to resolve ambiguous taxa are well documented and are known to be comparable. The suitability of the methods was evaluated and compared on the basis of 13 criteria that considered conservation of taxa richness and abundance, consistency among samples, methods, and studies, and effects on the interpretation of the data. Methods RPMC and DPAC had the highest suitability scores regardless of whether ambiguous taxa were resolved for each sample separately or for a group of samples. Method MCWP gave consistently poor results. Methods MCWP and DPAC approximate the use of family-level identifications and operational taxonomic units (OTU), respectively. Our results suggest that restricting identifications to the family level is not a good method of resolving ambiguous taxa, whereas generating OTUs works well provided that documentation issues are addressed.

Multilevel regression models describing regional patterns of invertebrate and algal responses to urbanization across the USA

Abstract Multilevel hierarchical regression was used to examine regional patterns in the responses of benthic macroinvertebrates and algae to urbanization across 9 metropolitan areas of the conterminous USA. Linear regressions established that responses (intercepts and slopes) to urbanization of invertebrates and algae varied among metropolitan areas. Multilevel hierarchical regression models were able to explain these differences on the basis of region-scale predictors. Regional differences in the type of land cover (agriculture or forest) being converted to urban and climatic factors (precipitation and air temperature) accounted for the differences in the response of macroinvertebrates to urbanization based on ordination scores, total richness, Ephemeroptera, Plecoptera, Trichoptera richness, and average tolerance. Regional differences in climate and antecedent agriculture also accounted for differences in the responses of salt-tolerant diatoms, but differences in the responses of other diatom metrics (% eutraphenic, % sensitive, and % silt tolerant) were best explained by regional differences in soils (mean % clay soils). The effects of urbanization were most readily detected in regions where forest lands were being converted to urban land because agricultural development significantly degraded assemblages before urbanization and made detection of urban effects difficult. The effects of climatic factors (temperature, precipitation) on background conditions (biogeographic differences) and rates of response to urbanization were most apparent after accounting for the effects of agricultural development. The effects of climate and land cover on responses to urbanization provide strong evidence that monitoring, mitigation, and restoration efforts must be tailored for specific regions and that attainment goals (background conditions) may not be possible in regions with high levels of prior disturbance (e.g., agricultural development).

Stream macroinvertebrate response models for bioassessment metrics: addressing the issue of spatial scale

We developed independent predictive disturbance models for a full regional data set and four individual ecoregions (Full Region vs. Individual Ecoregion models) to evaluate effects of spatial scale on the assessment of human landscape modification, on predicted response of stream biota, and the effect of other possible confounding factors, such as watershed size and elevation, on model performance. We selected macroinvertebrate sampling sites for model development (n = 591) and validation (n = 467) that met strict screening criteria from four proximal ecoregions in the northeastern U.S.: North Central Appalachians, Ridge and Valley, Northeastern Highlands, and Northern Piedmont. Models were developed using boosted regression tree (BRT) techniques for four macroinvertebrate metrics; results were compared among ecoregions and metrics. Comparing within a region but across the four macroinvertebrate metrics, the average richness of tolerant taxa (RichTOL) had the highest R2 for BRT models. Across the four metrics, final BRT models had between four and seven explanatory variables and always included a variable related to urbanization (e.g., population density, percent urban, or percent manmade channels), and either a measure of hydrologic runoff (e.g., minimum April, average December, or maximum monthly runoff) and(or) a natural landscape factor (e.g., riparian slope, precipitation, and elevation), or a measure of riparian disturbance. Contrary to our expectations, Full Region models explained nearly as much variance in the macroinvertebrate data as Individual Ecoregion models, and taking into account watershed size or elevation did not appear to improve model performance. As a result, it may be advantageous for bioassessment programs to develop large regional models as a preliminary assessment of overall disturbance conditions as long as the range in natural landscape variability is not excessive.

Fish and Invertebrate Flow-Biology Relationships to Support the Determination of Ecological Flows for North Carolina†

A method was developed to characterize fish and invertebrate responses to flow alteration in the state of North Carolina. This method involved using 80th percentile linear quantile regressions to relate six flow metrics to the diversity of riffle-run fish and benthic Ephemeroptera, Plecoptera, and Trichoptera (EPT) richness. All twelve flow-biology relationships were found to be significant, with both benthos and fish showing negative responses to ecodeficits and reductions in flow. The responses of benthic richness to reduced flows were consistent and generally greater than that of fish diversity. However, the riffle-run fish guild showed the greatest reductions in diversity in response to summer ecodeficits. The directional consistency and differential seasonal sensitivities of fish and invertebrates to reductions in flow highlight the need to consider seasonality when managing flows. In addition, all relationships were linear, and therefore do not provide clear thresholds to support ecological flow determinations and flow prescriptions to prevent the degradation of fish and invertebrate communities in North Carolina rivers and streams. A method of setting ecological flows based on the magnitude of change in biological condition that is acceptable to society is explored.

Aquatic Ecosystems as Indicators of Status and Trends in Water Quality

This chapter provides an introduction and overview of the use of aquatic ecosystems as indicators of water quality. The monitoring of biological communities (fish, algae, and invertebrates) as a means of assessing water quality conditions is emphasized along with ecotoxicological studies and measures of ecosystem function. Issues related to the design, implementation, and analysis of monitoring programs are reviewed and discussed with the objective of continually improving water quality conditions through adaptive management.