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

Patterns in and predictors of stream and river macroinvertebrate genera and fish species richness across the conterminous USA

Both native and non-native taxa richness patterns are useful for evaluating areas of greatest conservation concern. To determine those patterns, we analyzed fish and macroinvertebrate taxa richness data obtained at 3475 sites collected by the USEPA's National Rivers and Streams Assessment. We also determined which natural and anthropogenic variables best explained patterns in regional richness. Macroinvertebrate and fish richness increased with the number of sites sampled per region. Therefore, we determined residual taxa richness from the deviation of observed richness from predicted richness given the number of sites per region. Regional richness markedly exceeded average site richness for both macroinvertebrates and fish. Predictors of macroinvertebrate-genus and fish-species residual-regional richness differed. Air temperature was an important predictor in both cases but was positive for fish and negative for macroinvertebrates. Both natural and land use variables were significant predictors of regional richness. This study is the first to determine mean site and regional richness of both fish and aquatic macroinvertebrates across the conterminous USA, and the key anthropogenic drivers of regional richness. Thus, it offers important insights into regional USA biodiversity hotspots.

Choice of field and laboratory methods affects the detection of anthropogenic disturbances using stream macroinvertebrate assemblages

Accurate and precise detection of anthropogenic impacts on stream ecosystems using macroinvertebrates as biological indicators depends on the use of appropriate field and laboratory methods. We assessed the responsiveness to anthropogenic disturbances of assemblage metrics and composition by comparing commonly employed alternative combinations of field sampling and individuals counting methods. Four datasets were derived by, in the field 1) conducting multihabitat sampling (MH) or 2) targeting samples in leaf packs (single-habitat sampling - SH) and, in the laboratory A) counting all individuals of the samples, or B) simulating subsampling of 300 individuals per sample. We collected our data from 39 headwater stream sites in a drainage basin located in the Brazilian Cerrado. We used a previously published quantitative integrated disturbance index (IDI), based on both local and catchment disturbance measurements, to characterize the intensity of anthropogenic alterations at each site. Family richness and % Ephemeroptera, Plecoptera and Trichoptera (% EPT) individuals obtained from each dataset were tested against the IDI through simple linear regressions, and the differences in assemblage composition between least- and most-disturbed sites was tested using Permutational Multivariate Analysis of Variance (PERMANOVA). When counting all individuals, differences in taxonomic richness and assemblage composition of macroinvertebrate assemblages between least- and most-disturbed sites were more pronounced in the MH than in the SH sampling method. Leaf packs seemed to concentrate high abundance and diversity of macroinvertebrates in highly disturbed sites, acting as 'biodiversity hotbeds' in these situations, which likely reduced the response of the assemblages to the disturbance gradient when this substrate was targeted. However, MH sampling produced weaker results than SH when subsampling was performed. The % EPT individuals responded better to the disturbance gradient when SH was employed, and its efficiency was not affected by the subsampling procedure. We conclude that no single method was the best in all situations, and the efficiency of a sampling protocol depends on the combination of field and laboratory methods being used. Although the total count of individuals with multihabitat sampling obtained the best results for most of the evaluated variables, the decision of which procedures to use depends on the amount of time and resources available, on the variables of interest, on the availability of habitat types in the sites sampled, and on the other methods being employed in the sampling protocol.

The role of physical habitat and sampling effort on estimates of benthic macroinvertebrate taxonomic richness at basin and site scales

Taxonomic richness is one of the most important measures of biological diversity in ecological studies, including those with stream macroinvertebrates. However, it is impractical to measure the true richness of any site directly by sampling. Our objective was to evaluate the effect of sampling effort on estimates of macroinvertebrate family and Ephemeroptera, Plecoptera, and Trichoptera (EPT) genera richness at two scales: basin and stream site. In addition, we tried to determine which environmental factors at the site scale most influenced the amount of sampling effort needed. We sampled 39 sites in the Cerrado biome (neotropical savanna). In each site, we obtained 11 equidistant samples of the benthic assemblage and multiple physical habitat measurements. The observed basin-scale richness achieved a consistent estimation from Chao 1, Jack 1, and Jack 2 richness estimators. However, at the site scale, there was a constant increase in the observed number of taxa with increased number of samples. Models that best explained the slope of site-scale sampling curves (representing the necessity of greater sampling effort) included metrics that describe habitat heterogeneity, habitat structure, anthropogenic disturbance, and water quality, for both macroinvertebrate family and EPT genera richness. Our results demonstrate the importance of considering basin- and site-scale sampling effort in ecological surveys and that taxa accumulation curves and richness estimators are good tools for assessing sampling efficiency. The physical habitat explained a significant amount of the sampling effort needed. Therefore, future studies should explore the possible implications of physical habitat characteristics when developing sampling objectives, study designs, and calculating the needed sampling effort.

Louisiana Waterthrush and Benthic Macroinvertebrate Response to Shale Gas Development

Because shale gas development is occurring over large landscapes and consequently is affecting many headwater streams, an understanding of its effects on headwater-stream faunal communities is needed. We examined effects of shale gas development (well pads and associated infrastructure) on Louisiana waterthrush Parkesia motacilla and benthic macroinvertebrate communities in 12 West Virginia headwater streams in 2011. Streams were classed as impacted (n = 6) or unimpacted (n = 6) by shale gas development. We quantified waterthrush demography (nest success, clutch size, number of fledglings, territory density), a waterthrush Habitat Suitability Index, a Rapid Bioassessment Protocol habitat index, and benthic macroinvertebrate metrics including a genus-level stream-quality index for each stream. We compared each benthic metric between impacted and unimpacted streams with a Student's t-test that incorporated adjustments for normalizing data. Impacted streams had lower genus-level stream-quality index scores; lower overall and Ephemeroptera, Plecoptera, and Trichoptera richness; fewer intolerant taxa, more tolerant taxa, and greater density of 0–3-mm individuals (P ≤ 0.10). We then used Pearson correlation to relate waterthrush metrics to benthic metrics across the 12 streams. Territory density (no. of territories/km of stream) was greater on streams with higher genus-level stream-quality index scores; greater density of all taxa and Ephemeroptera, Plecoptera, and Trichoptera taxa; and greater biomass. Clutch size was greater on streams with higher genus-level stream-quality index scores. Nest survival analyses (n = 43 nests) completed with Program MARK suggested minimal influence of benthic metrics compared with nest stage and Habitat Suitability Index score. Although our study spanned only one season, our results suggest that shale gas development affected waterthrush and benthic communities in the headwater streams we studied. Thus, these ecological effects of shale gas development warrant closer examination.

Land use change in the Atlantic Forest affects carbon and nitrogen sources of streams as revealed by the isotopic composition of terrestrial invertebrates

Terrestrial invertebrates link terrestrial systems to aquatic ones, making vegetal material produced in the watershed available to aquatic food webs. In this study, using carbon and nitrogen stable isotopes, we evaluated the importance of introduced C4 grasses as a source of carbon in aquatic food webs of headwater streams of the coastal Atlantic Forest located on the north coast of the State of São Paulo, in the southeastern region of Brazil. Terrestrial invertebrates were collected in two streams: one where the main land cover was pristine montane Atlantic Forest (forest stream) and another where the main land cover was introduced C4 forage grasses for livestock (pasture stream). The average δ13C of terrestrial invertebrates collected in the forest stream (−26.3±2.1‰) was significantly (p<0.01) smaller than the average δ13C of terrestrial invertebrates collected in the pasture stream (−15.7±4.7‰), denoting a larger contribution of C4grasses to terrestrial invertebrates of the pasture stream. The average δ15N of terrestrial invertebrates of the forest stream (4.1±2.4‰) was significantly (p<0.01) lower than the average δ15N of terrestrial invertebrates of the pasture stream (9.5±2.7‰). The relative contribution of C3 and C4 plants to terrestrial invertebrates was estimated using SIAR. In the forest stream, the C3 contribution was on average 0.75 (0.72 minimum to 0.79 maximum), and the C4 contribution was on average 0.25 (0.21 minimum to 0.28 maximum). In the pasture stream, the C3contribution decreased to 0.20 (0.14 minimum to 0.26 maximum), and the C4 contribution increased to 0.80 (0.74 minimum to 0.86 maximum). These results have several implications for the ecosystem functioning as well as for recent changes in environmental policies of Brazil. The lower nutritional value of C4 grasses may not only decrease invertebrate performance, but also alter the stoichiometry of several components of the aquatic food webs with potential consequence for the whole ecosystem functioning. On the public policy side, recent changes in the Brazilian Forest Act, a series of laws that regulate land cover at the property level, reduced the width of the forested riparian area with potentially dangerous consequences for aquatic ecosystems.

Stream biomonitoring using macroinvertebrates around the globe: a comparison of large-scale programs

Water quality agencies and scientists are increasingly adopting standardized sampling methodologies because of the challenges associated with interpreting data derived from dissimilar protocols. Here, we compare 13 protocols for monitoring streams from different regions and countries around the globe. Despite the spatially diverse range of countries assessed, many aspects of bioassessment structure and protocols were similar, thereby providing evidence of key characteristics that might be incorporated in a global sampling methodology. Similarities were found regarding sampler type, mesh size, sampling period, subsampling methods, and taxonomic resolution. Consistent field and laboratory methods are essential for merging data sets collected by multiple institutions to enable large-scale comparisons. We discuss the similarities and differences among protocols and present current trends and future recommendations for monitoring programs, especially for regions where large-scale protocols do not yet exist. We summarize the current state in one of these regions, Latin America, and comment on the possible development path for these techniques in this region. We conclude that several aspects of stream biomonitoring need additional performance evaluation (accuracy, precision, discriminatory power, relative costs), particularly when comparing targeted habitat (only the commonest habitat type) versus site-wide sampling (multiple habitat types), appropriate levels of sampling and processing effort, and standardized indicators to resolve dissimilarities among biomonitoring methods. Global issues such as climate change are creating an environment where there is an increasing need to have universally consistent data collection, processing and storage to enable large-scale trend analysis. Biomonitoring programs following standardized methods could aid international data sharing and interpretation.

Long-term impacts on macroinvertebrates downstream of reclaimed mountaintop mining valley fills in Central Appalachia

Recent studies have documented adverse effects to biological communities downstream of mountaintop coal mining and valley fills (VF), but few data exist on the longevity of these impacts. We sampled 15 headwater streams with VFs reclaimed 11-33 years prior to 2011 and sampled seven local reference sites that had no VFs. We collected chemical, habitat, and benthic macroinvertebrate data in April 2011; additional chemical samples were collected in September 2011. To assess ecological condition, we compared VF and reference abiotic and biotic data using: (1) ordination to detect multivariate differences, (2) benthic indices (a multimetric index and an observed/expected predictive model) calibrated to state reference conditions to detect impairment, and (3) correlation and regression analysis to detect relationships between biotic and abiotic data. Although VF sites had good instream habitat, nearly 90 % of these streams exhibited biological impairment. VF sites with higher index scores were co-located near unaffected tributaries; we suggest that these tributaries were sources of sensitive taxa as drifting colonists. There were clear losses of expected taxa across most VF sites and two functional feeding groups (% scrapers and %shredders) were significantly altered. Percent VF and forested area were related to biological quality but varied more than individual ions and specific conductance. Within the subset of VF sites, other descriptors (e.g., VF age, site distance from VF, the presence of impoundments, % forest) had no detectable relationships with biological condition. Although these VFs were constructed pursuant to permits and regulatory programs that have as their stated goals that (1) mined land be reclaimed and restored to its original use or a use of higher value, and (2) mining does not cause or contribute to violations of water quality standards, we found sustained ecological damage in headwaters streams draining VFs long after reclamation was completed.

Towards a protocol for stream macroinvertebrate sampling in China

Standard protocols are critical for maximizing data comparability and aggregation in national monitoring programs, and taxa richness is a common indicator of site condition and biological diversity. There are two general approaches for sampling stream macroinvertebrate assemblages: targeted richest habitat and site wide. At seven sites, we compared three methods: Ontario Benthic Biomonitoring Network (OBBN), Environmental Monitoring and Assessment Program (EMAP), and Rapid Bioassessment Protocol (RBP). The OBBN method produced a biased sample at a site with a single small riffle, the RBP method produced the most total taxa, and the EMAP method produced the most taxa at four sites and the most individuals at six sites. The RBP method produced asymptotes for percent tolerant individuals, percent chironomid individuals, and Hilsenhoff Biotic Index score after five to ten stations. The EMAP method produced asymptotes for those metrics after 10 to 20 stations per site. The EMAP method typically required half the number of stations as the RBP method to obtain 70-90% of true taxa richness as estimated by the Jaccard coefficient. We conclude that the EMAP method is preferable because of its greater precision in taxa richness estimates.