Bioassessment Data Research Articles

Perspectives from the Global South inform and affirm the contemporary river conservation paradigm: A commentary on Feio et al.

Perspectives from the Global South inform and affirm the contemporary river conservation paradigm: A commentary on Feio et al.

Eutrophication thresholds associated with protection of biological integrity in California wadeable streams

Eutrophication is one of the most pervasive stressors impacting streams, often leading to loss of biodiversity or change in natural functions. To protect against these adverse effects, managers can set targets for environmental indicators to limit eutrophication that are likely to maintain high biological integrity. To identify protective management targets, we evaluated the responses of three bioassessment indices (one for benthic macroinvertebrates and two for benthic algae assemblages) to five eutrophication indicators (total nitrogen [TN], total phosphorus [TP], benthic chlorophyll-a [chl-a], benthic ash-free dry mass [AFDM], and percent macroalgal cover [% cover] of the streambed). First, we used a bioassessment data set of 1249 sites in California to create logistic regression models of the likelihood of achieving several biointegrity goals for each index along increasing gradients of each indicator. Then, we evaluated eutrophication thresholds at concentrations corresponding to several relative probabilities (from 50 to 95 %), reflecting a range of policy makers’ potential tolerance for risk of failing to meet biointegrity goals. Finally, we validated the thresholds with relative risk assessment, and identified the lowest validated threshold across the three indices. All eutrophication indicators were significantly associated with increased risks to biointegrity, resulting in a set of validated thresholds for each biointegrity goal. For example, thresholds of 0.24 mg/L TN, 0.05 mg/L TP, 44 mg/m2 benthic chl-a, 2.5 mg/cm2 AFDM, and 26 % cover would achieve index scores above the 10th percentile of reference (a biointegrity goal that has been used in certain regulatory applications) with an 80 % probability. If these thresholds were applied ambient monitoring statewide dataset, TN was the most pervasively exceeded threshold (37 % of sites), particularly within the highly agricultural Central Valley (76 %), although exceedances of TN and AFDM were also common in the urbanized South Coast region (68 % and 72 % of sites, respectively).

Prioritizing Stream Protection, Restoration and Management Actions Using Landscape Modeling and Spatial Analysis

Watersheds are often degraded by human activities, reducing their ability to provide ecosystem functions and services. While governmental agencies have put forward plans for improving watershed health, resources are limited, and choices must be made as to which watersheds to prioritize and what actions to take. Prioritization tools with sufficient specificity, resolution, and automation are needed to guide decisions on restoration and management actions across large scales. To address this need, we developed a set of tools to support the protection of streams and associated riparian habitats across the state of California. We developed and tested watershed condition estimation models based on bioassessment data, used the EPA’s StreamCat dataset to identify stressors, incorporated environmental justice factors and developed reach-specific models to prioritize actions. We applied the prioritization tools statewide and were able to identify 18% of stream reaches that are in good condition but that are most vulnerable to existing stressors and an additional 19% of stream reaches that are degraded and are highest priority for restoration and management. The remaining 63% of stream reaches were prioritized for protection and periodic monitoring or minor remedial actions. The results of this project can help regional stakeholders and agencies prioritize hundreds of millions of dollars being spent to protect, acquire, and restore stream and riparian habitats. The methods are directly transferable by using any regional condition and stress data that can be readily obtained.

Application of Flow-Ecology Analysis to Inform Prioritization for Stream Restoration and Management Actions

A key challenge in managing flow alteration is determining the severity and pattern of alteration associated with the degradation of biological communities. Understanding these patterns helps managers prioritize locations for restoration and flow management actions. However, the choices made about how to use these flow-ecology relationships can have profound implications on management decisions (e.g., which biological endpoints, which thresholds, which seasonal flow components to use). We describe a process for using flow-ecology relationships to prioritize management actions that 1) Represents the most relevant components of the annual hydrograph, 2) Demonstrates an appropriate level of sensitivity in order to discriminate locations to inform decision making, 3) Aims to protect multiple biological assemblages, 4) Reduces misclassification of priority areas (i.e., error of omission). Our approach is based on the functional flows approach which uses multiple flow metrics that describe the frequency, timing, magnitude, duration, and rate of change of seasonal process-based components of the annual hydrograph. Using this approach, we performed a flow-ecology analysis of regional bioassessment data, through which we determined where flow alteration impacts biology and prioritized reaches for changes in flow management to protect aquatic resources in a highly urbanized region of southern California, where managing scarce water resources leads to difficult decisions about tradeoffs that require technical information. We identified three important functional flow metrics for each of two bioassessment indices, one based on benthic macroinvertebrates, and another based on benthic algae. Based on thresholds that describe levels of alteration as well as thresholds describing the probability of achieving a healthy biological condition, we compared nine biological threshold combinations for each index. We found instances of flow alteration that impact biological condition highly variable (0–100% of subbasins) between combinations and we present a method for finding the most appropriate combination for prioritizing locations for flow management. We apply the final thresholds to the study region and propose 16 subbasins of high priority for implementing flow management and restoration. Importantly, we show that focusing on a single biological group would result in biologically altered locations being effectively ignored.

A web-based tool for assessing the condition of benthic diatom assemblages in streams and rivers of the conterminous United States.

Benthic diatom assemblages are known to be indicative of water quality but have yet to be widely adopted in biological assessments in the United States due to several limitations. Our goal was to address some of these limitations by developing regional multi-metric indices (MMIs) that are robust to inter-laboratory taxonomic inconsistency, adjusted for natural covariates, and sensitive to a wide range of anthropogenic stressors. We aggregated bioassessment data from two national-scale federal programs and used a data-driven analysis in which all-possible combinations of 2-7 metrics were compared for three measures of performance. After ranking the best-performing MMIs, we selected the final MMIs by evaluating stress-response relations in independent regional datasets of diatom samples paired with measures of several water-quality stressors, including herbicides and streamflow flashiness. Each regional MMI performed well at calibration sites and represented diverse aspects of the structure and function of diatom communities. Most metrics included in the best MMIs were modeled to account for natural variation including climate, topography, soil characteristics, lithology, and groundwater influence on streamflow. MMI performance improved with higher numbers of component metrics, but this effect diminished beyond six metrics. Component metrics of MMIs were associated with a broad suite of measured stressors in every region, including salinity, nutrients, herbicides, and streamflow flashiness. We provide a web-based software application that allows users in the conterminous United States to apply our MMIs to their own datasets and compare MMI scores from their sites to a broader regional context.

Identifying Functional Flow Linkages Between Stream Alteration and Biological Stream Condition Indices Across California

Large state or regional environmental flow programs, such as the one based on the California Environmental Flows Framework, rely on broadly applicable relationships between flow and ecology to inform management decisions. California, despite having high flow and bioassessment data density, has not established relationships between specific elements of the annual hydrograph and biological stream condition. To address this, we spatially and temporally linked USGS gage stations and biological assessment sites in California to identify suitable paired sites for comparisons of streamflow alteration with biological condition at a statewide scale. Flows were assessed using a set of functional flow metrics that provide a comprehensive way to compare alteration and seasonal variation in streamflow across different locations. Biological response was evaluated using the California Stream Condition Index (CSCI) and Algal Stream Condition Index (ASCI), which quantify biological conditions by translating benthic invertebrate or algal resources and watershed-scale environmental data into an overall measure of stream health. These indices provide a consistent statewide standard for interpreting bioassessment data, and thus, a means of quantitatively comparing stream conditions throughout the state. The results indicate that indices of biological stream condition were most closely associated with flow alteration in seasonality and timing metrics, such as fall pulse timing, dry-season timing, and wet season timing. Magnitude metrics such as dry-season baseflow, wet season baseflow, and the fall pulse magnitude were also important in influencing biological stream conditions. Development of ecological flow needs in large-scale environmental programs should consider that alteration to any of the seasonal flow components (e.g., dry-season baseflow, fall pulse flow, wet-season baseflow, spring recession flow) may be important in restructuring biological communities.

Prioritizing management goals for stream biological integrity within the developed landscape context

Stream management goals for biological integrity may be difficult to achieve in developed landscapes where channel modification and other factors constrain in-stream conditions. To evaluate potential constraints on biological integrity, we developed a statewide landscape model for California that estimates ranges of likely scores for a macroinvertebrate-based index that are typical at a site with the observed level of landscape alteration. This context can support prioritization decisions for stream management, like identifying reaches for restoration or enhanced protection based on how observed scores relate to model predictions. Median scores were accurately predicted by the model for all sites in California with bioassessment data (Pearson correlation r = 0.75 between observed and predicted for calibration data, r = 0.72 for validation). The model also predicted that 15% of streams statewide are constrained for biological integrity within their present developed landscape, particularly for urban and agricultural areas in the South Coast, Central Valley, and Bay Area regions. We worked with a local stakeholder group from the San Gabriel River watershed (Los Angeles County, California) to evaluate how the statewide model could support local management decisions. To achieve this purpose, we created an interactive application, the Stream Classification and Priority Explorer, that compares observed scores with predictions from the landscape model to assign priorities. We observed model predictions consistent with the land-use gradient from the upper to lower watershed, where potential limits to achieving biological integrity were more common in the heavily-urbanized lower watershed. However, most of the sites in the lower watershed scored within their predicted ranges, and were therefore given a low priority for restoration. In contrast, 2 low-scoring sites in the undeveloped upper watershed were prioritized for causal assessment and possible future restoration, whereas 3 high-scoring sites were prioritized for protection. The availability of geospatial and bioassessment data at the national level suggests that these tools can easily be applied to inform management decisions at other locations where altered landscapes may limit biological integrity.

Empirically-based modeling and mapping to consider the co-occurrence of ecological receptors and stressors

Part of the ecological risk assessment process involves examining the potential for environmental stressors and ecological receptors to co-occur across a landscape. In this study, we introduce a Bayesian joint modeling framework for use in evaluating and mapping the co-occurrence of stressors and receptors using empirical data, open-source statistical software, and Geographic Information Systems tools and data. To illustrate the approach, we apply the framework to bioassessment data on stream fishes and nutrients collected from a watershed in southwestern Ohio. The results highlighted the joint model's ability to parse and exploit statistical dependencies in order to provide empirical insight into the potential environmental and ecotoxicological interactions influencing co-occurrence. We also demonstrate how probabilistic predictions can be generated and mapped to visualize spatial patterns in co-occurrences. For practitioners, we believe that this data-driven approach to modeling and mapping co-occurrence can lead to more quantitatively transparent and robust assessments of ecological risk.

Evaluation of suitability and comparability of stream assessment indices using macroinvertebrate data sets from the Northern Lakes and Forests Ecoregion

Researchers and managers within the Upper Midwest currently use a variety of sampling methodologies and biological indices to assess ecological condition of stream systems. With multiple entities collecting bioassessment data it is important that we determine the comparability of data and the indices derived from these data for effective assessment of natural systems. In this study we assessed the similarity of data collected by different agencies and we focused on data from one watershed to examine the outputs of different indices for stream assessment, and the temporal variation of index score within sites. We compared duplicate macroinvertebrate community data collected by the Little River Band of Ottawa Indians and the Michigan Department of Environmental Quality for overall community composition and index scores derived from these data. Duplicate samples were similar in composition index scores. Taxonomic resolution was addressed and indicated that genus level resolution gives a more favorable score when using indices. We also evaluated the utility of currently available macroinvertebrate indices of biotic integrity to assess data from the Big Manistee River watershed. The indices evaluated were the Hilsenhoff biotic index, the benthic community index for the Northern Lakes and Forests (NLFBCI), the Great Lakes Environmental Assessment Survey (GLEAS) procedure 51 for macroinvertebrates and a biological condition gradient model for the Upper Midwest. Outputs from the indices were moderately correlated (Spearman rank order correlation, r=0.35–0.698) though they indicated different assessments of overall site integrity. Compared with larger scale regional indices, locally calibrated indices generally classified sites as having better biological condition. Replicate samples collected within sites indicated the GLEAS had higher levels of variability (0–265%CV) within sites than the other indices (<10%CV). Data from long-term (10 year) monitoring stations were used to evaluate seasonal and long-term index performance. There were differences in index score classifications from spring and fall samples indicating that standardization of sampling time is necessary for comparative analysis. Temporal trends over 10 years reveal natural variation and set the baseline for evaluating the influence of anthropogenic effects. Overall, results indicate that choice of index can alter assessment of site condition. For bioassessment in the Big Manistee River watershed the NLFBCI performs well and accurately reflects site condition.

An approach for determining bioassessment performance and comparability

Many organizations in the USA collect aquatic bioassessment data using different sampling and analysis methods, most of which have unknown performance in terms of data quality produced. Thus, the comparability of bioassessments produced by different organizations is often unknown, ultimately affecting our ability to make comprehensive assessments on large spatial scales. We evaluated a pilot approach for determining bioassessment performance using macroinvertebrate data obtained from several states in the Southeastern USA. Performance measures evaluated included precision, sensitivity, and responsiveness to a human disturbance gradient, defined in terms of a land disturbance index value for each site, combined with a value for specific conductance, and instream habitat quality. A key finding of this study is the need to harmonize ecoregional reference conditions among states so as to yield more comparable and consistent bioassessment results. Our approach was also capable of identifying potential areas for refinement such as reevaluation of less precise, sensitive, or responsive metrics that may result in suboptimal index performance. Higher performing bioassessments can yield information beyond "impaired" versus "unimpaired" condition. Acknowledging the limitations of this pilot study, we would recommend that performance evaluations use at least 50 sites, 10 of which are ecoregional reference sites. Efforts should be made to obtain data from the entire human disturbance gradient in an ecoregion to improve statistical confidence in performance measures. Having too few sites will result in an under-representation of certain parts of the disturbance gradient (e.g., too few poor quality sites), which may bias sensitivity and responsiveness estimates.

A study of the effect of physical and chemical stressors on biological integrity within the San Diego hydrologic region

Environmental agencies across the United States have searched for adequate methods to assess anthropogenic impacts on the environment. Biological assessments, which compare the taxonomic composition of an aquatic assemblage to relevant biocriteria, have surfaced as an effective method to assess the ecological integrity of US waterbodies. In this study, bioassessment data were collected and analyzed in conjunction with physical habitat and chemical stressor data for streams and rivers within the San Diego basin from 1998 through 2005. Physical stressors such as sediment loading, riparian destruction, and in-stream habitat homogenization affect many locations in the region. However, physical habitat measures alone were found to frequently overestimate the biological integrity of streams in the region. Many sites within the San Diego Basin, although unaffected by physical stressors, continue to exhibit low biological integrity scores. Sites with low biological integrity tend to possess higher specific conductance and salinity compared to sites with high biological integrity. We suggest that one possible reason for these differences is the source water used for municipal purposes.

Non-biting midges in biodiversity conservation and environmental assessment: Findings from boreal freshwater ecosystems

Non-biting midges (Diptera: Chironomidae) are frequently overlooked in freshwater biodiversity surveys and environmental assessment, yet they are commonly the most abundant and the most diverse taxon in freshwater ecosystems. We reviewed the diversity patterns and assemblage–environment relationships of non-dipterans and chironomids mainly in boreal freshwater ecosystems building on previously reported findings. Although generally the same environmental gradients are correlated with assemblage structure, their relative importance varies between chironomids and non-dipterans. Chironomid assemblage response to and recovery from human impacts are also likely to differ from that of other common benthic taxa. Thus, environmental assessments may be biased if chironomids are not included. Different surrogacy approaches have thus far shown little success in accounting for chironomid species richness and assemblage structure, and there appears to be no easy short-cut for the examination of chironomids as part of freshwater surveys. However, we show that genus-level identification of pupal exuviae provides a reliable and rapid way of estimating chironomid species richness at least in boreal freshwater ecosystems. In addition, we demonstrate that the inclusion of chironomids may increase the signal-to-noise ratio in bioassessment data sets, and that this information can be obtained with modest increases in costs.

Vulnerability of biological metrics and multimetric indices to effects of climate change

Aquatic ecosystems and their fauna are vulnerable to a variety of climate-related changes. Benthic macroinvertebrates are used frequently by water-quality agencies to monitor the status of aquatic resources. We used several regionally distributed state bioassessment data sets to analyze how climate change might influence metrics used to define ecological condition of streams. Many widely used, taxonomically based metrics were composed of both cold- and warm-water-preference taxa, and differing responses of these temperature-preference groups to climate-induced changes in stream temperatures could undermine assessment of stream condition. Climate responsiveness of these trait groups varied among states and ecoregions, but the groups generally were sensitive to changing temperature conditions. Temperature sensitivity of taxa and their sensitivity to organic pollution were moderately but significantly correlated. Therefore, metrics selected for condition assessments because taxa are sensitive to disturbance or to conventional pollutants also were sensitive to changes in temperature. We explored the feasibility of modifying metrics by partitioning components based on temperature sensitivity to reduce the likelihood that responses to climate change would confound responses to impairment from other causes and to facilitate tracking of climate-change-related taxon losses and replacements.

Evaluation of the Biocriteria of Streams in the San Diego Hydrologic Region

In this study, bioassessment data collected between 1998 and 2005 were synthesized and analyzed for streams and rivers throughout the San Diego Hydrologic Region to provide a spatial and temporal context for the results of several monitoring projects conducted between 1998 and 2005 and to ascertain the applicability of the Southern California benthic macroinvertebrate index of biological integrity (SoCal B-IBI) to the region’s streams. The water quality of the sites studied in the region, as reflected by temporal and spatial analyses of SoCal B-IBI scores, was found to be quite poor. When streams were analyzed individually most showed stable scores over the time frame of the study with some showing better scores in the fall. Spatially, scores were found to be better farther away from the coast in the upstream reaches of the watersheds. This study further explored the applicability of the SoCal B-IBI to a focused geographic region by demonstrating the necessity of each component metric to the assignment of biological condition. Although all component metric scores were deemed to be necessary, the percent intolerant individuals score had a more significant effect in driving impairment. The analysis of the component metrics of the SoCal B-IBI provides useful insights to the changes in scores among and between the sampled sites in the region’s watersheds. Based on this study, natural resource management agencies responsible for managing water quality should incorporate regular measures of biological integrity into their water quality programs to ascertain regional and temporal trends.

Relationship of fish and macroinvertebrate communities in the mid-Atlantic uplands: Implications for integrated assessments

Indices developed for stream bioassessment are typically based on either fish or macroinvertebrate assemblages. These indices consist of metrics which subsume attributes of various species into aggregate measures reflecting community-level ecological responses to disturbance. However, little is known about the relationship between fish and macroinvertebrate metrics, or about how ecological health assessments are affected by assemblage-specific responses to disturbance. We used principal component analysis (PCA) and regression analysis of existing fish (n = 371) and macroinvertebrate (n = 442) stream bioassessment data from a multi-source dataset to determine broad scale, within-assemblage metric patterns, and to examine the intercorrelation of fish and macroinvertebrate metrics (n = 246) and their response to watershed area and land use/land cover gradients. Fish and macroinvertebrate metrics expressed as principal components (PCs) accounted for 72.4 and 85.4% of dataset variance, respectively, with PC-metric patterns reflecting aspects of stream impairment including water and habitat quality. Model components predicting fish metric response differed among fish PCs, with watershed area and macroinvertebrate metric response strongly correlated with the first fish PC, and remaining fish PC models consisting of watershed area, land use, and macroinvertebrate PCs. Correlation between fish and macroinvertebrate PCs, and models relating fish and macroinvertebrate PCs generally explained less variation (13–27%) than metric response models of fish (25–34%) and macroinvertebrates (8–38%) to watershed area and land use/land cover variables. Best-response models integrating fish and macroinvertebrate PCs, watershed area, and land use/land cover variables accounted for the greatest variation in fish PCs (32–50%) across sites. Because fish and macroinvertebrate metrics provide different information on ecological condition, integrated use of information from multiple groups may be appropriate when developing monitoring programs.

Bioassessment of the ecological integrity of river ecosystems using aquatic macroinvertebrates: an overview with a focus on South Africa

An overview is given of the bioassessment of the ecological integrity of river ecosystems using aquatic macroinvertebrates, focussing on the South African situation within a global context. Biotic indices and their use in aquatic bioassessment are covered, and comparative descriptions of the more important and/or widely used macroinvertebrate-based biotic indices are given. For each biotic index described, information is provided on the biotopes sampled, the sampling equipment used, the sampling protocol followed, the level of taxonomic resolution, whether identifications are laboratory- or field-based, the range of the final index value, and its current usage. A discussion on using biotic indices and interpreting bioassessment data follows. The multimetric and multivariate approaches to bioassessment are compared, as are the regional and multivariate approaches to the derivation of reference conditions. The use of invertebrate habitat assessment in river bioassessment studies is briefly discussed and the South African spatial framework for analysing bioassessment data is described. In conclusion, a number of potential avenues for further research regarding the use of macroinvertebrates in the bioassessment of river ecosystems are identified.

Natural variation in macroinvertebrate assemblages and the development of a biological banding system for interpreting bioassessment data—a preliminary evaluation using data from upland sites in the south-western Cape, South Africa

The variability of macroinvertebrate assemblages was investigated at 27 upland reference sites in the south-western Cape, South Africa. Multivariate analyses showed that sites did not group on the basis of geomorphological zonation, i.e. mountain stream and foothill-cobble bed. When separate analyses were undertaken for mountain stream (n = 21) and foothill-cobble bed sites (n = 6), assemblages formed three and two groups, respectively. Similarity amongst groups ranged from 47% to 52%, while within-group similarity was between 54% and 67%. Environmental variables shown to contribute to this variability included distance from source, cation ratio ([Na+]+[K+]/([Na+]+[K+]+[Ca2+]+[Mg2+]), pH, longitude and stream width. Whilst overall variability in the metrics of the biotic index, SASS (South African Scoring System), is high at reference sites, the interpretation of monitoring-site data using biological bands derived from a range of reference sites, ensured that variability was taken into account and that detection of disturbance at a monitoring site was not impeded. A biological banding system has been developed for upland sites in the south-western Cape, together with a list of reference or expected SASS-taxa. This list includes details pertaining to seasonality and biotope preferences. The ability to define reference conditions that take intrinsic variability amongst reference sites into account is important for the accurate interpretation of bioassessment data.

Estimation and application of indicator values for common macroinvertebrate genera and families of the United States

Abstract Tolerance of macroinvertebrate taxa to chemical and physical stressors is widely used in the analysis and interpretation of bioassessment data, but many estimates lack empirical bases. Our main objective was to estimate genus- and family-level indicator values (IVs) from a data set of macroinvertebrate communities, chemical, and physical stressors collected in a consistent manner throughout the United States. We then demonstrated an application of these IVs to detect alterations in benthic macroinvertebrate assemblages along gradients of urbanization in New England and Alabama. Principal components analysis (PCA) was used to create synthetic gradients of chemical stressors, for which genus- and family-level weighted averages (WAs) were calculated. Based on results of PCA, WAs were calculated for three synthetic gradients (ionic concentration, nutrient concentration, and dissolved oxygen/water temperature) and two uncorrelated physical variables (suspended sediment concentration and percent fines). Indicator values for each stress gradient were subsequently created by transforming WAs into ten ordinal ranks based on percentiles of values across all taxa. Mean IVs of genera and families were highly correlated to road density in Alabama and New England, and supported the conclusions of independent assessments of the chemical and physical stressors acting in each geographic area. Family IVs were nearly as responsive to urbanization as genus IVs. The limitations of widespread use of these IVs are discussed.

Seasonal variability of macroinvertebrate assemblages in two regions of South Africa: implications for aquatic bioassessment

Seasonal variability of macroinvertebrate assemblages was examined in two regions of South Africa: Western Cape and Mpumalanga. Sampling was undertaken at reference sites using the rapid bioassessment method SASS4 (South African Scoring System, Version 4). This study examined the influence of sampling season on the occurrence of taxa, macroinvertebrate assemblages and SASS scores. In the Western Cape a few taxa were more common at certain times of the year, and multivariate analysis of macroinvertebrate assemblages showed that assemblages grouped by season. The number of taxa and the Average Score per Taxon (ASPT) were significantly different among seasons in the Western Cape, with fewer taxa recorded in winter compared to summer, and significantly higher ASPT values recorded in winter and spring than summer and autumn. Whereas more taxa were recorded in autumn than in spring, a higher proportion of sensitive and high-scoring taxa were recorded in spring. In Mpumalanga macroinvertebrate assemblages did not group by season, and differences in SASS scores were not significant, although more taxa were recorded in winter than in spring or autumn, and ASPTs were slightly higher in winter than in spring. Results are discussed in relation to ecological reference conditions and the interpretation of bioassessment data.

Modification of an ecotoxicological rating to bioassess small acid mine drainage‐impacted watersheds exclusive of benthic macroinvertebrate analysis

The modified ecotoxicological rating (METR) approach to synthesizing integrative bioassessment data into a single number ranging from 0 to 100 was found to differentiate between multiple levels of impacts from acid mine drainage (AMD). Our objective was to develop a more cost-effective and time-efficient bioassessment technique than previously used in other large scale ecotoxicological ratings (ETRs) by minimizing the number of parameters required to rank stations to only those most descriptive of the benthic macroinvertebrate community responses to AMD. Nineteen physical, chemical, toxicological, and ecological measurements were made at 38 stations in two adjacent watersheds. The most descriptive parameters were selected through multiple linear regression analysis, bivariate correlation analysis, and one-way analysis of variance. We found that habitat assessment, 30-d in situ Asian clam survival, mean conductivity, and mean total water column concentration of aluminum and manganese were the most descriptive parameters. The METR constructed from these parameters was equally effective at differentiating stations as were two previous published ETRs that incorporated up to 10 parameters, including benthic macroinvertebrate indices. When the METR was applied to a new watershed, the scores were significantly correlated with benthic macroinvertebrate indices for those stations.