Estimating vertebrate,benthic macroinvertebrate,and diatom taxa richness in raftable Pacific Northwest rivers for bioassessment purposes

The number of sites sampled must be considered when determining the effort necessary for adequately assessing taxa richness in an ecosystem for bioassessment purposes; however, there have been few studies concerning the number of sites necessary for bioassessment of large rivers. We evaluated the effect of sample size (i.e., number of sites) necessary to collect vertebrate (fish and aquatic amphibians), macroinvertebrate, and diatom taxa from seven large rivers in Oregon and Washington, USA during the summers of 2006–2008. We used Monte Carlo simulation to determine the number of sites needed to collect 90–95% of the taxa 75–95% of the time from 20 randomly located sites on each river. The river wetted widths varied from 27.8 to 126.0 m, mean substrate size varied from 1 to 10 cm, and mainstem distances sampled varied from 87 to 254 km. We sampled vertebrates at each site (i.e., 50 times the mean wetted channel width) by nearshore-raft electrofishing. We sampled benthic macroinvertebrates nearshore through the use of a 500-μm mesh kick net at 11 systematic stations. From each site composite sample, we identified a target of 500 macroinvertebrate individuals to the lowest possible taxon, usually genus. We sampled benthic diatoms nearshore at the same 11 stations from a 12-cm² area. At each station, we sucked diatoms from soft substrate into a 60-ml syringe or brushed them off a rock and rinsed them with river water into the same jar. We counted a minimum of 600 valves at 1,000× magnification for each site. We collected 120–211 diatom taxa, 98–128 macroinvertebrate taxa, and 14–33 vertebrate species per river. To collect 90-95% of the taxa 75-95% of the time that were collected at 20 sites, it was necessary to sample 11–16 randomly distributed sites for vertebrates, 13–17 sites for macroinvertebrates, and 16–18 sites for diatoms. We conclude that 12–16 randomly distributed sites are needed for cost-efficient sampling of vertebrate richness in the main stems of our study rivers, but 20 sites markedly underestimates the species richness of benthic macroinvertebrates and diatoms in those rivers.     

Using Information on Spatial Variability of Small Estuaries in Designing Large-Scale Estuarine Monitoring Programs

In the early 1990s, EPA's Environmental Monitoring and Assessment Program (EMAP) documented the ecological condition of the overall population of small estuaries along the mid-Atlantic coast of the United States. However, the Program did not provide detailed information on the condition of individual estuaries less than 260 km² in surface area, a group of estuaries of concern to environmental managers. To address the needs of environmental managers, when EMAP returned to the region in summer 1997, it included a study of the spatial variability of ecological indicators within individual small estuaries. At 127 probability-based sites in 10 estuaries, EMAP measured a variety of parameters of water quality and sediments, including dissolved oxygen (DO), nutrients, grain size of sediments, contaminants in sediments, and community structures of benthic macroinvertebrates. From this information the ecological condition (e.g., percent area with DO concentrations below 5mg L^–1) for each estuary, along with 90% confidence interval, was determined. The width of the confidence interval was then recalculated for sample sizes ranging from two stations to the total number of stations sampled in that estuary. Confidence interval widths were then plotted against sample size. These plots can be useful in designing future regional monitoring programs with a goal of describing conditions in individual systems as well as broad geographic regions. Results illustrate that beyond five stations per estuary, the reduction in the width of the confidence interval with increasing sampling intensity is relatively small; however, individual program managers need to determine "how small is small enough."     

Mayfly (Insecta: Ephemeroptera) community structure as an indicator of the ecological status of a stream in the Niger Delta area of Nigeria

Ephemeroptera is an important group of insects used in the bioassessment and monitoring of freshwater bodies worldwide because of their relative abundance in a wide variety of substrates and their increasing chances of detecting pollution impacts. In this study, their faunistic composition and spatiotemporal variations in density and diversity in River Orogodo (Southern Nigeria) was investigated at five ecologically distinct stations over a 12-month period. The mayfly nymph community responses to environmental variables were evaluated by means of biological measures and multivariate analysis (redundancy analysis [RDA]). Thirteen morphologically distinct taxa belonging to six families were identified. The dominant taxa were Afrobaetodes pusillus (23.1%), Baetis sp. (13.7%), and Caenis cibaria (11.4%). The density of Ephemeroptera differed significantly (p?<?0.05) both in space and time. Diversity was influenced by substrate heterogeneity which in turn was influenced by catchment processes such as flooding and anthropogenic activities especially abattoir effluent. Based on the RDA ordination and relative abundance data, Baetis sp. dominated at impacted stations while a more equitable distribution of species were observed in less disturbed sites. Water velocity, canopy cover, nature of bottom sediments, and the amount of dissolved oxygen also accounted for the variations in Ephemeroptera densities at the different stations. Shannon diversity, taxa richness, and evenness were lowest in station 3 (the abattoir discharge site).     

The Ecotoxicological Recovery of Ely Creek and Tributaries (Lee County, VA) after Remediation of Acid Mine Drainage

The Ely Creek watershed (Lee County, VA) was determined in 1995 to be the most negatively affected by acid mine drainage (AMD) within the Virginia coalfield. This determination led the US Army Corps of Engineers to design and build passive wetland remediation systems at two major AMD seeps affecting Ely Creek. This study was undertaken to determine if ecological recovery had occurred in Ely Creek. The results indicate that remediation had a positive effect on all monitoring sites downstream of the remediated AMD seeps. At the site most impacted by AMD, mean pH was 2.93 prior to remediation and improved to 7.14 in 2004. Benthic macroinvertebrate surveys revealed that one AMD influenced site had increased taxa richness from zero taxa in 1997 to 24 in 2004. While in situ testing of Asian clams resulted in zero survival at five of seven AMD influenced sites prior to remediation, some clams survived at all sites after. Clam survival was found to be significantly less than upstream references at only two sites, both downstream of un-mitigated AMD seeps in 2004. An ecotoxicological rating (ETR) system that combined ten biotic and abiotic parameters was developed as an indicator of the ecological status for each study site. A comparison of ETRs from before and after remediation demonstrated that all sites downstream of the remediation had experienced some level of recovery. Although the remediation has improved the ecological health of Ely Creek, un-mitigated AMD discharges are still negatively impacting the watershed.     

The Ecotoxicological Recovery of Ely Creek and Tributaries (Lee County, Va) After Remediation of Acid Mine Drainage

The Ely Creek watershed (Lee County, VA) was determined in 1995 to be the most negatively affected by acid mine drainage (AMD) within the Virginia coalfield. This determination led the US Army Corps of Engineers to design and build passive wetland remediation systems at two major AMD seeps affecting Ely Creek. This study was undertaken to determine if ecological recovery had occurred in Ely Creek. The results indicate that remediation had a positive effect on all monitoring sites downstream of the remediated AMD seeps. At the site most impacted by AMD, mean pH was 2.93 prior to remediation and improved to 7.14 in 2004. Benthic macroinvertebrate surveys revealed that one AMD influenced site had increased taxa richness from zero taxa in 1997 to 24 in 2004. While in situ testing of Asian clams resulted in zero survival at five of seven AMD influenced sites prior to remediation, some clams survived at all sites after. Clam survival was found to be significantly less than upstream references at only two sites, both downstream of un-mitigated AMD seeps in 2004. An ecotoxicological rating (ETR) system that combined ten biotic and abiotic parameters was developed as an indicator of the ecological status for each study site. A comparison of ETRs from before and after remediation demonstrated that all sites downstream of the remediation had experienced some level of recovery. Although the remediation has improved the ecological health of Ely Creek, un-mitigated AMD discharges are still negatively impacting the watershed.     

The problem of using fixed-area subsampling methods to estimate macroinvertebrate richness: a case study with Neotropical stream data

Subsampling has been widely applied in the laboratory to process freshwater macroinvertebrate samples. Currently, many governmental agencies and research groups apply the fixed-count approach, targeting a number of individuals per sample, and at the same time keeping track of the number of quadrats (fraction of the sample) processed. However, fixed-area methods are still in use. The objective of this paper was to evaluate the reliability of macroinvertebrate taxonomic richness estimates developed from processing a standard number of subsampling quadrats (i.e., fixed-area approaches). We used a dataset from 18 tropical stream sites experiencing three different levels of human disturbance (most-, intermediate-, and least-disturbed). With 12 quadrats processed (half the sample), the collection curves started to stabilize, and for more than half of the sites studied, it was possible to sample at least 80 % of the total taxonomic richness of the sample. However, we observed that the minimum number of quadrats to achieve 80 % of taxonomic richness was strongly negatively correlated with the number of individuals collected in each site: the fewer the individuals in a sample, the greater the processed proportion of that sample needed to represent it properly. Thus our results indicate that for any given areal subsampling effort (any fixed fraction of the sample), samples with different numbers of individuals will be represented differently in terms of the proportion of the total number of taxa of the whole samples, those with greater numbers being overestimated and those with fewer numbers being underestimated. Therefore, we do not recommend the use of fixed-area subsampling methods alone if the main purpose is to measure and analyze taxonomic richness; instead, we encourage researchers to use fixed-count approaches for this purpose.     

A comparison of techniques to sample salamander assemblages along highland streams of Maryland

Amphibians may be useful indicators of biological condition in small streams so determining which sampling technique maximizes encounters at the least cost and at the optimal time of year is important. Area constrained surveys (ACS), used by the Maryland Biological Stream Survey, were tested against cover board surveys, drift fences with pitfall and funnel traps, quadrat leaf litter searches, and leaf litter bags. Sixteen, 100 m-long sites were established in headwater streams in the Savage River State Forest in Garrett County, Maryland. Each technique was randomly assigned to a 25 m stream section within each overall sampling site, and sites were sampled once each month from May to October (2005) with additional sampling in March and April (2006). Area constrained surveys yielded means of 2.7 taxa and 14.9 total individuals per sampling visit, which was significantly higher than the yield of all other methods in all months except October and March, when yields were low for all techniques. Area constrained surveys were also significantly more cost-effective per taxon and per individual compared to all other methods. September produced the most taxa and individuals, October and March produced the least, and yields for April through August were similar to September. We employed removal sampling at four sites in April 2006, but abundance could not be estimated because a significant linear decrease in the accumulated catch versus catch per unit effort did not occur for three of the sites.     

Benthic Macroinvertebrate Responses to Increasing Levels of Cattle Grazing in Blue Ridge Mountain Streams, Virginia, USA

The relationship between benthic macroinvertebrate assemblages and cattle density was assessed from fall 2002 through spring 2004 in five small streams that represented a gradient of cattle grazing intensity. All study stream reaches were in pasture with no woody riparian vegetation, but varied in the intensity of cattle grazing (0 cattle ha⁻¹ at site 1 to 2.85 cattle ha⁻¹ at site 5). Regression analysis indicated highly significant and strong macroinvertebrate metric responses to cattle density during most sampling periods. The majority of metrics responded negatively to increased grazing, while a few (total taxa richness, number of sensitive taxa, and % collector filterers) increased along the gradient before declining at the most heavily grazed sites. Total number of sensitive taxa and % Coleoptera had the strongest relationship with cattle density throughout the study period. During some sampling periods, nearly 80% of the variation in these metrics was explained by cattle density. The elmid beetle, Oulimnius, had a particularly strong negative response to the grazing gradient. Study site groupings based on taxa composition, using detrended correspondence analysis (DCA), indicated that benthic samples collected from the reference site and light rotational grazing site were more similar in macroinvertebrate taxa composition than samples collected from the intermediate grazing and heavy grazing sites. Our findings demonstrate that biological integrity, as measured by benthic macroinvertebrate metrics and assemblage composition, is highly related to cattle density in small streams in the Blue Ridge mountains, Virginia, USA. This suggests that the degree of agricultural intensity should be given consideration in stream assessments, as well as land use planning and regulatory decisions.     

Linking Monitoring and Effects Research: EMAP's Intensive Site Network Program

The EMAP program has been organized into three primary elements: Multi-Tier Design, Indicators, and Index Sites. The Index Sites program (DISPro - Demonstration Intensive Site Project) is the primary activity within the Index Sites element of EMAP. This project represents an inter-agency effort between EPA/ORD and DOI/NPS to develop a demonstration of an intensive site network of monitoring and research locations throughout the United States, utilizing the Nation's parklands as "outdoor laboratories." Twelve parks were selected to establish this demonstration. These 12 parks were selected because they are readily accessible, have a history of monitoring environmental information, and represent a broad spectrum of ecological communities. EMAP, through DISPro, is examining whether a "network" of sites existing within the parks can be used to address monitoring issues for global-scale environmental stressors (e.g., air pollution) as well as locale-specific stressors (e.g., air deposition, water-borne) and coordinated with cause-effect, issue-based research related to these environmental stressors. As a first activity, EPA will provide each of the sites with the instrumentation to monitor UV-B. The intent of the program is to initiate a consistent air monitoring program at each site to be followed by consistent monitoring within other media. The project will initiate research projects at all the sites (eventually) to examine the effects of environmental stressors of importance at each of the sites.     

Quality comparison of tap water vs. bottled water in the industrial city of Yanbu (Saudi Arabia)

The objective of this study was to determine the relative influence of reach-specific habitat variables and geographic location on benthic invertebrate assemblages within six ecoregions across the Western USA. This study included 417 sites from six ecoregions. A total of 301 taxa were collected with the highest richness associated with ecoregions dominated by streams with coarse substrate (19-29 taxa per site). Lowest richness (seven to eight taxa per site) was associated with ecoregions dominated by fine-grain substrate. Principle component analysis (PCA) on reach-scale habitat separated the six ecoregions into those in high-gradient mountainous areas (Coast Range, Cascades, and Southern Rockies) and those in lower-gradient ecoregions (Central Great Plains and Central California Valley). Nonmetric multidimensional scaling (NMS) models performed best in ecoregions dominated by coarse-grain substrate and high taxa richness, along with coarse-grain substrates sites combined from multiple ecoregions regardless of location. In contrast, ecoregions or site combinations dominated by fine-grain substrate had poor model performance (high stress). Four NMS models showed that geographic location (i.e. latitude and longitude) was important for: (1) all ecoregions combined, (2) all sites dominated by coarse-grain sub strate combined, (3) Cascades Ecoregion, and (4) Columbia Ecoregion. Local factors (i.e. substrate or water temperature) seem to be overriding factors controlling invertebrate composition across the West, regardless of geographic location.     

Assessment of Benthic Infaunal Condition on the Mainland Shelf of Southern California

Benthic infauna were sampled from 251 Southern California Bight (SCB) mainland shelf sites in the summer of 1994. Sample sites were selected using a stratified random design, with the primary strata being depth zone, geography, and proximity to point and non-point discharges. Benthic infaunal condition was assessed using the Benthic Response Index (BRI), and by comparing dominant taxa and community parameters (e.g., number of taxa) among strata. Ninety-one percent of sediments in the SCB were found to contain healthy benthic communities. Most stations with altered benthos were located near river mouths, in Santa Monica Bay, or on the Palos Verdes Shelf. Deviations at sites with altered benthic communities were mostly limited to minor changes in species composition, rather than to large declines in diversity or abundance.     

A Bird Community Index of Biotic Integrity for the Mid-Atlantic Highlands

We report on the development and preliminary application of a songbird community-based index of biotic integrity. The bird community index (BCI) sorts bird species found at sample sites into a series of values representing the proportional species richness of 20 behavioral and physiological response guilds. Relative proportions of specialist and generalist guilds are used to assign a composite score to each site. Scores from multiple sites indicate the overall biotic integrity of the study area. The BCI is intended to function as a landscape-scale indicator of biotic integrity, integrating conditions across large sample sites containing diverse ecological resources and intensities of human use. We developed the BCI with data from a 1994 pilot study in central Pennsylvania, then applied our preliminary index in 1995 and 1996 to independent samples of sites across the Mid-Atlantic Highlands Assessment area (MAHA). The 1995 and 1996 sample sites were selected using the probability-based sampling design of the Environmental Monitoring and Assessment Program (EMAP), and therefore represent the total land area in MAHA. Our preliminary assessment indicates that MAHA exhibits six categories of biotic integrity, and that more than 40% of the land area supports the two highest biotic integrity categories. Pending BCI refinement and incorporation of landscape and vegetation explanatory variables, the BCI will be included in a suite of indicators designed to provide an assessment of overall ecological condition in MAHA.     

The Ecological Condition of Gulf of Mexico Resources from Perdido Key to Port St. Joe, Florida, USA: Part I. Coastal Beach Resources

Using the approach established by EPA's Environmental Monitoring and Assessment Program (EMAP), a shoreline monitoring survey was conducted in August and September 1999, encompassing the Florida Panhandle from Perdido Key, Florida to Port St. Joe, Florida. The objective of this survey was to demonstrate the use of a probabilistic survey for monitoring and estimating the condition of swimmable beach areas. Thirty stations were sampled using a probabilistic sampling design. Hydrographic data were collected in addition to samples for water chemistry. Bacterial indicators, enterococci and fecal coliforms, were enumerated from the water according to the EPA Beaches Environmental Assessment Closure and Health (BEACH) Program and Florida state guidelines. Additional criteria for site condition included the presence or absence of primary and secondary dunes, anthropogenic debris and vegetation. Based on EMAP evaluation guidelines and Florida state criteria, a baseline assessment of the condition of the Gulf of Mexico beach resources surveyed is presented.     

Use of PFU protozoan community structural and functional characteristics in assessment of water quality in a large, highly polluted freshwater lake in China

Structural and functional parameters of protozoan communities colonizing on PFU (polyurethane foam unit) artificial substrate were assessed as indicators of water quality in the Chaohu Lake, a large, shallow and highly polluted freshwater lake in China. Protozoan communities were sampled 1, 3, 6, 9 and 14 days after exposure of PFU artificial substrate in the lake during October 2003. Four study stations with the different water quality gradient changes along the lake were distinguishable in terms of differences in the community's structural (species richness, individual abundance, etc.) and functional parameters (protozoan colonization rates on PFU). The concentrations of TP, TN, COD and BOD as the main chemical indicators of pollution at the four sampling sites were also obtained each year during 2002-2003 for comparison with biological parameters. The results showed that the species richness and PFU colonization rate decreased as pollution intensity increased and that the Margalef diversity index values calculated at four sampling sites also related to water quality. The three functional parameters based on the PFU colonization process, that is, S(eq), G and T90%, were strongly related to the pollution status of the water. The number of protozoan species colonizing on PFU after exposure of 1 to 3 days was found to give a clear comparative indication of the water quality at the four sampling stations. The research provides further evidence that the protozoan community may be utilized effectively in the assessment of water quality and that the PFU method furnishes rapid, cost-effective and reliable information that may be useful for measuring responses to pollution stress in aquatic ecosystems.     

Relating Benthic Infaunal Community Structure to Environmental Variables in Estuaries using Nonmetric Multidimensional Scaling and Similarity Analysis

In 1994, 19 stations were sampled (2 replicates/station) with Young grabs in association with the EMAP-Estuaries Carolian Province Base Monitoring in Florida. A total of 295 unique benthic infauna taxa and 9647 individuals were identified and enumerated. Environmental data (bottom-water quality, sediment grain size, sediment metals, and organics) and benthic community data were analyzed using hierarchical agglomerative cluster analysis and ordination via nonmetric multidimensional scaling. Bray-Curtis similarities and Euclidean distance were used as the distance measures for biotic and abiotic data, respectively. Multivariate analyses were complemented by examining incremental contributions of benthic taxa to similarity values using a percentage similarity technique. A low-salinity site in a tributary to the St. Johns River had benthic communities uniquely different from those of moderate-to high-salinity sites. A diverse assemblage of polychaetes, gastropods, bivalves, amphipods, sipunculans, and phoronids was consistently associated with relatively unimpacted sites in the Indian River Lagoon. Infaunal community structure in the northern portion of the study area was influenced by the nearby Atlantic Ocean. Community shifts in association with latitudinal gradients and concentrations of sediment metals and organics were apparent. The nonparametric multivariate techniques used in this study were particularly effective at delineating and defining fine-scale community differences.     

Sediment Quality of Estuaries in the Southeastern U.S.

A study was conducted to assess the condition of estuaries in the EMAP Carolinian Province (Cape Henry, VA - St. Lucie Inlet, FL). Synoptic measures of sediment contamination, toxicity, and macroinfaunal condition were made at 82 and 86 stations in 1994 and 1995, respectively, in accordance with a probabilistic sampling design. These data were used to estimate percentages of degraded vs. undegraded estuarine area from the perspective of sediment quality. Each year a sizable portion of the province (36% in 1994, 51 % in 1995) showed some evidence of either degraded benthic assemblages, contaminated sediments in excess of bioeffect guidelines, or significant sediment toxicity (based on Ampelisca abdita and Microtox® assays). However, co-occurrences of a degraded benthos and adverse exposure conditions (sediment contamination and/or toxicity) were much less extensive – 17% of the province in 1994 and 25% in 1995. Each year only four sites, representing 5% of the province in 1994 and 8% in 1995, had degraded infauna accompanied by both sediment contamination and toxicity, suggesting that strong contaminant-induced effects on the benthos (based on such combined weight-of-evidence) were limited to a fairly small percentage of estuarine area province-wide. PCBs and pesticides (lindane, dieldrin, DDT and derivatives) were the most dominant contaminants over the two-year period. The broad-scale sampling design of EMAP was not intended to support detailed characterizations of potential pollutant impacts within individual estuaries. Thus, some estuaries classified as undegraded may include additional degraded portions outside the immediate vicinity of randomly sampled sites. Such localized impacts (not accounted for in the above estimates) were detected in this study at additional nonrandom supplemental sites near potential contaminant sources.     

Development and validation of a macroinvertebrate index of biotic integrity (IBI) for assessing urban impacts to Northern California freshwater wetlands

Despite California policies requiring assessment of ambient wetland condition and compensatory wetland mitigations, no intensive monitoring tools have been developed to evaluate freshwater wetlands within the state. Therefore, we developed standardized, wadeable field methods to sample macroinvertebrate communities and evaluated 40 wetlands across Northern California to develop a macroinvertebrate index of biotic integrity (IBI). A priori reference sites were selected with minimal urban impacts, representing a best-attainable condition. We screened 56 macroinvertebrate metrics for inclusion in the IBI based on responsiveness to percent urbanization. Eight final metrics were selected for inclusion in the IBI: percent three dominant taxa; scraper richness; percent Ephemeroptera, Odonata, and Trichoptera (EOT); EOT richness; percent Tanypodinae/Chironomidae; Oligochaeta richness; percent Coleoptera; and predator richness. The IBI (potential range 0–100) demonstrated significant discriminatory power between the reference (mean = 69) and impacted wetlands (mean = 28). It also declined with increasing percent urbanization (R ² = 0.53, p < 0.005) among wetlands in an independent validation dataset (n = 14). The IBI was robust in showing no significant bias with environmental gradients. This IBI is a functional tool to determine the ecological condition at urban (stormwater and flood control ponds), as well as rural freshwater wetlands (stockponds, seasonal wetlands, and natural ponds). Biological differences between perennial and non-perennial wetlands suggest that developing separate indicators for these wetland types may improve applicability, although the existing data set was not sufficient for exploring this option.     

松花江流域大型底栖动物生物完整性指数构建及其适用性

利用生物完整性指数评价河流健康状态,对于水环境管理决策具有重要的实践意义。基于大型底栖动物构建生物完整性指数(B-IBI),并评价松花江流域的水生态系统健康状况。在松花江主要干支流设定37个采样点,分别于2016年6、9月进行环境因子和大型底栖动物调查研究。最终从28个候选参数中确定了种类总数、摇蚊种类数、敏感种百分比、Hilsenhoff指数、Marglef指数作为核心参数构建B-IBI。通过0～10赋分法,计算得到了松花江流域全部采样点的生物完整性评价得分。结果显示,松花江流域内60%区域生物状态存在不同程度的损害。另外,B-IBI能够综合反映松花江大型底栖动物群落多样性、生境质量、理化水质等,具有一定的适用性。     

Evaluation of the Impact of Acid Mine Drainage on the Chemistry and the Macrobenthos in the Carolina Stream (San Luis – Argentina)

The purpose of this study was to evaluate the impact of acid mine drainage on the chemistry and the macrobenthos of the Carolina stream (San Luis – Argentina). Samples were obtained in the years 1997–1998 at two sites: site C₁, located 200 m upstream of the drainage, and site C₂, located 800 m downstream. The system buffer capacity was evaluated in the non – contaminated site by means of the buffer index calculation. The physico – chemical changes observed as a consequence of the contribution of acid mine drainage (AMD) were: a decreasing of pH and an increase in the ionic concentration, especially sulfate and Fe coming from the oxidation produced by chemiolithotrophic bacteria. The values obtained indicated a low buffer capacity and a high intrinsic vulnerability of the system to resist the impact originated by the AMD, producing a remarkable decreasing of pH of the receiving stream. These changes caused modifications in the original benthic community that was replaced by organisms more tolerant to the acid stress. A reduction in the abundance and in the taxonomic richness of the benthic macroinvertebrates was observed when compared with the reference station. An increase in the proportion of Chironomidae and of Acari and a decrease in the proportion of the remaining taxa were observed. The most sensitive groups were Ephemeroptera, Trichoptera and Mollusca. The community was mostly affected by the following variables: pH, conductivity, sulfate and dissolved total Fe.