Distribution and Characteristics of Blackwater Rivers and Streams of the Contiguous United States.

Blackwater rivers and streams have stained or tea-colored water from tannins released by decaying plant matter. Natural conditions in these waters often differ from non-blackwater systems. For example, the pH and oxygen levels in waters can be very low, but completely natural. We examined an existing USEPA data set and found that blackwaters existed across the contiguous United States but were most common in the east. Water chemistry showed differences between blackwater and non-blackwater sites, but differences were not consistent across ecoregions making national scale generalizations difficult. Physical habitat data analysis did not show dramatic differences between blackwater and non-blackwater sites. Blackwater typically arises from streams that drain tannin-rich bogs/muskeg and wetlands, so as expected blackwater sites had a shorter Euclidean distance to wetlands than non-blackwater sites and existed in watersheds with a higher percentage of wetland habitat. Blackwaters in Northern and Temperate Plains tended to have higher acid neutralizing capacity, conductivity, and lower True Color; a visual color scale used for water purity. We posit that differences were because Color and Dissolved Organic Carbon at these sites were from buried wetland deposits rather than contemporary wetland habitats. Research needs that may increase our understanding and management of blackwaters include development of an operational definition that includes a classification framework and reference conditions for different blackwater types, identification of stressors and their associated dynamics that negatively impact blackwater systems, and development of data-driven, consistent, and repeatable assessment methods, including development of targets, protective of unique conditions in blackwater rivers and streams.

A national-scale vegetation multimetric index (VMMI) as an indicator of wetland condition across the conterminous United States.

In 2011, the US Environmental Protection Agency and its partners conducted the first National Wetland Condition Assessment at the continental-scale of the conterminous United States. A probability design for site selection was used to allow an unbiased assessment of wetland condition. We developed a vegetation multimetric index (VMMI) as a parsimonious biological indicator of ecological condition applicable to diverse wetland types at national and regional scales. Vegetation data (species presence and cover) were collected from 1138 sites that represented seven broad estuarine intertidal and inland wetland types. Using field collected data and plant species trait information, we developed 405 candidate metrics with potential for distinguishing least disturbed (reference) from most disturbed sites. Thirty-five of the metrics passed range, repeatability, and responsiveness screens and were considered as potential component metrics for the VMMI. A permutation approach was used to calculate thousands of randomly constructed potential national-scale VMMIs with 4, 6, 8, or 10 metrics. The best performing VMMI was identified based on limited redundancy among constituent metrics, sensitivity, repeatability, and precision. This final VMMI had four broadly applicable metrics (floristic quality index, relative importance of native species, richness of disturbance-tolerant species, and relative cover of native monocots). VMMI values and weights from the survey design for probability sites (n = 967) were used to estimate wetland area in good, fair, and poor condition, nationally and for each of 10 ecoregion by wetland type reporting groups. Strengths and limitations of the national VMMI for describing ecological condition are highlighted.

Characterizing nonnative plants in wetlands across the conterminous United States.

Nonnative plants are widely recognized as stressors to wetlands and other ecosystems. They may compete with native plant species or communities and alter ecosystem properties, which can affect ecological condition, posing challenges to resource managers. As part of the United States Environmental Protection Agency's National Wetland Condition Assessment (NWCA), we characterized the status of nonnative plants in wetlands across the conterminous United States (US). Our primary goals were to (1) document the composition of nonnative taxa at 1138 NWCA sites sampled in 2011 and (2) estimate the areal extent of wetland under stress from nonnative plants within the NWCA 2011 sampled population of ~ 25 million ha of wetland (represented by 967 sampled probability sites and the NWCA survey design). A total of 443 unique nonnative taxa were observed, encompassing a species pool adapted to diverse ecological conditions. For individual sites, the number of nonnative taxa ranged from 0 to 29, and total absolute cover of nonnatives ranged from 0 to 160%. We devised the nonnative plant indicator (NNPI) as a categorical indicator of stress (low to very high) from the collective set of nonnative plant taxa occurring at a particular location, based on a decision matrix of exceedance values for nonnative richness, relative frequency, and relative cover. Wetland area of the sampled population occurring in each NNPI category was estimated at the scale of the conterminous US and within five large ecoregions and four broad wetland types. Potential stress from nonnative plants, as indicated by the NNPI category, was low for approximately 61% (~ 15.3 million ha), moderate for about 20% (~ 5.2 million ha), high for about 10% (~ 2.48 million ha), and very high for about 9% (~ 2.2 million ha) of the wetland area in the entire sampled population. Percent of wetland area with high and very high NNPI varied by ecoregional subpopulations: greater within interior and western ecoregions (~ 29 to 87%) than within ecoregions in the eastern half of the nation (~ 11%). Among wetland type subpopulations, greater percent of wetland area with high and very high NNPI was observed for herbaceous vs. woody types and for inland vs. estuarine types. Estimates of wetland area by NNPI categories are expected to be useful to policy makers or resource managers for prioritizing management actions by identifying situations where stress from nonnative plants is most extensive. We also considered four exploratory analyses aimed at providing ecological information useful in interpreting NNPI extent results. We conducted three population-scale analyses examining ecoregional and wetland type population means for (1) the three NNPI metrics, (2) absolute cover of growth-habit groups of nonnative plants, and (3) metrics describing human-mediated disturbance. Finally, we examined ecological relationships with site-level NNPI status using a random forest (RF) analysis with NNPI as the response variable and predictor variables including ecoregion, wetland type, and a variety of characteristics describing natural vegetation structure, environment, and human-mediated disturbance.

Use of national-scale data to examine human-mediated additions of heavy metals to wetland soils of the US.

Soil concentrations of 12 heavy metals that have been linked to various anthropogenic activities were measured in samples collected from the uppermost horizon in approximately 1000 wetlands across the conterminous US as part of the 2011 National Wetland Condition Assessment (NWCA). The heavy metals were silver (Ag), cadmium (Cd), cobalt (Co), chromium (Cr), copper (Cu), nickel (Ni), lead (Pb), antimony (Sb), tin (Sn), vanadium (V), tungsten (W), and zinc (Zn). Using thresholds to distinguish natural background concentrations from human-mediated additions, we evaluated wetland soil heavy metal concentrations in the conterminous US and four regions using a Heavy Metal Index (HMI) that reflects human-mediated heavy metal loads based on the number of elements above expected background concentration. We also examined the individual elements to detect concentrations of heavy metals above expected background that frequently occur in wetland soils. Our data show that wetland soils of the conterminous US typically have low heavy metal loads, and that most of the measured elements occur nationally in concentrations below thresholds that relate to anthropogenic activities. However, we found that soil lead is more common in wetland soils than other measured elements, occurring nationally in 11.3% of the wetland area in concentrations above expected natural background (> 35 ppm). Our data show positive relationships between soil lead concentration and four individual landscape metrics: road density, percent impervious surface, housing unit density, and population density in a 1-km radius buffer area surrounding a site. These relationships, while evident on a national level, are strongest in the eastern US, where the highest road densities and greatest population densities occur. Because lead can be strongly bound to wetland soils in particular, maintenance of the good condition of our nation's wetlands is likely to minimize risk of lead mobilization.

Benthic macroinvertebrate field sampling effort required to produce a sample adequate for the assessment of rivers and streams of Neuquén Province, Argentina.

This multi-year pilot study evaluated a proposed field method for its effectiveness in the collection of a benthic macroinvertebrate sample adequate for use in the condition assessment of streams and rivers in the Neuquén Province, Argentina. A total of 13 sites, distributed across three rivers, were sampled. At each site, benthic macroinvertebrates were collected at 11 transects. Each sample was processed independently in the field and laboratory. Based on a literature review and resource considerations, the collection of 300 organisms (minimum) at each site was determined to be necessary to support a robust condition assessment, and therefore, selected as the criterion for judging the adequacy of the method. This targeted number of organisms was collected at all sites, at a minimum, when collections from all 11 transects were combined. Subsequent bootstrapping analysis of data was used to estimate whether collecting at fewer transects would reach the minimum target number of organisms for all sites. In a subset of sites, the total number of organisms frequently fell below the target when fewer than 11 transects collections were combined.Site conditions where < 300 organisms might be collected are discussed. These preliminary results suggest that the proposed field method results in a sample that is adequate for robust condition assessment of the rivers and streams of interest. When data become available from a broader range of sites, the adequacy of the field method should be reassessed.

Validation of rapid assessment methods to determine streamflow duration classes in the Pacific Northwest, USA.

United States Supreme Court rulings have created uncertainty regarding U.S. Clean Water Act (CWA) authority over certain waters, and established new data and analytical requirements for determining CWA jurisdiction. Thus, rapid assessment methods are needed that can differentiate between ephemeral, intermittent, and perennial streams. We report on the validation of several methods. The first (Interim Method) was developed through best professional judgment (BPJ); an alternative (Revised Method) resulted from statistical analysis. We tested the Interim Method on 178 study reaches in Oregon, and constructed the Revised Method based on statistical analysis of the Oregon data. Next, we evaluated the regional applicability of the methods on 86 study reaches across a variety of hydrologic landscapes in Washington and Idaho. During the second phase, we also compared the Revised Method with a similar approach (Combined Method) based on combined field data from Oregon, Washington, and Idaho. We further compared field-based methods with a GIS-based approach (GIS Method) that used the National Hydrography Dataset and a synthetic stream network. Evaluations of all methods compared results with actual streamflow duration classes. The Revised Method correctly determined known streamflow duration 83.9% of the time, versus 62.3% accuracy of the Interim Method and 43.6% accuracy for the GIS-based approach. The Combined Method did not significantly outperform the Revised Method. Analysis showed biological indicators most accurately discriminate streamflow duration classes. While BPJ established a testable hypothesis, this study illustrates the importance of quantitative field testing of rapid assessment methods. Results support a consistent method applicable across the Pacific Northwest.

Evaluation of an alternate method for sampling benthic macroinvertebrates in low-gradient streams sampled as part of the National Rivers and Streams Assessment.

Benthic macroinvertebrates are sampled in streams and rivers as one of the assessment elements of the US Environmental Protection Agency's National Rivers and Streams Assessment. In a 2006 report, the recommendation was made that different yet comparable methods be evaluated for different types of streams (e.g., low gradient vs. high gradient). Consequently, a research element was added to the 2008-2009 National Rivers and Streams Assessment to conduct a side-by-side comparison of the standard macroinvertebrate sampling method with an alternate method specifically designed for low-gradient wadeable streams and rivers that focused more on stream edge habitat. Samples were collected using each method at 525 sites in five of nine aggregate ecoregions located in the conterminous USA. Methods were compared using the benthic macroinvertebrate multimetric index developed for the 2006 Wadeable Streams Assessment. Statistical analysis did not reveal any trends that would suggest the overall assessment of low-gradient streams on a regional or national scale would change if the alternate method was used rather than the standard sampling method, regardless of the gradient cutoff used to define low-gradient streams. Based on these results, the National Rivers and Streams Survey should continue to use the standard field method for sampling all streams.

Sediment toxicity in mid-continent great rivers (USA).

As part of the Environmental Monitoring and Assessment Program for Great River Ecosystems (EMAP-GRE), sediment samples were collected from 447 randomly selected littoral sites along the main channels of the Ohio, Missouri, and Upper Mississippi Rivers between 2004 and 2006. Toxicity of these sediment samples was measured using a 7-day Hyalella azteca survival and growth test. Sixty-five sites (14.5%) exhibited lethal toxicity, and 130 sites (29.1%) exhibited decreased growth. In the EMAP-GRE probabilistic sampling design, each sampled site had a weight associated with it that determined the length (and proportion) of the river represented by that sample point in the population. Weighted whole-river estimates indicated that of the 4721 river km sampled, sediment from 15.9 ± 3.0% of the river (752 ± 50 km) were lethally toxic, 27.4 ± 3.5% (1289 ± 57 km) were toxic by way of growth inhibition, and 40.0 ± 3.7% (1887 ± 68 km) exhibited either lethal or growth toxicity. Selected toxic samples were analyzed for 21 pesticides, 20 polychlorinated biphenyl congeners, and 6 polybrominated diphenyl ether congeners. For all of the samples tested, the concentration levels of these analytes were mostly lower than known toxicity thresholds, and neither unionized ammonia concentration nor osmotic stress (as measured by conductivity) could account for the toxicity found in sediments. The spatial pattern of sediment toxicity cannot be readily explained by urbanization or agricultural land use at the subcatchment scale. We speculate that the distribution of toxic sediment is more likely due to a combination of localized sources, including polluted tributaries, and the redistribution of contaminated sediments from upriver. The sediment toxicity results from this study will be used, in combination with other sediment, biologic, and habitat metrics and indicators collected in the EMAP-GRE study, to help interpret and assess the condition of the Ohio, Upper Mississippi, and Missouri Rivers.

Proteomic analysis of a model fish species exposed to individual pesticides and a binary mixture.

Pesticides are nearly ubiquitous in surface waters of the United States, where they often are found as mixtures. The molecular mechanisms underlying the toxic effects of sub-lethal exposure to pesticides as both individual and mixtures are unclear. The current work aims to identify and compare differentially expressed proteins in brains of male fathead minnows (Pimephales promelas) exposed for 72 h to permethrin (7.5 µg/L), terbufos (57.5 µg/L) and a binary mixture of both. Twenty-four proteins were found to be differentially expressed among all three treatments relative to the control using an ANOVA followed by a Dunnett's post hoc test (p ≤0.05). One protein was found to be differentially expressed among all treatment groups and one protein was in common between the terbufos and the mixture group. Fifteen spots were successfully sequenced using LC-MS/MS sequencing. Proteins associated with the ubiquitin-proteasome system, glycolysis, the cytoskeleton and hypoxia were enriched. As a second objective, we attempted to establish protein expression signatures (PES) for individual permethrin and terbufos exposures. We were unable to generate a useable PES for terbufos; however, the permethrin PES was able to distinguish between control and permethrin-exposed individuals in an independent experiment with an accuracy of 87.5%. This PES also accurately classified permethrin exposed individuals when the exposure occurred as part of a mixture. The identification of proteins differentially expressed as a result of pesticide exposure represent a step forward in the understanding of mechanisms of toxicity of permethrin and terbufos. They also allow a comparison of molecular responses of the binary mixture to single exposures. The permethrin PES is the first step in establishing a method to determine exposures in real-world scenarios.

Mercury contamination in fish in midcontinent great rivers of the United States: importance of species traits and environmental factors.

We measured mercury (Hg) concentrations in whole fish from the Upper Mississippi, Missouri, and Ohio Rivers to characterize the extent and magnitude of Hg contamination and to identify environmental factors influencing Hg accumulation. Concentrations were generally lower (80% of values between 20-200 ng g(-1) wet weight) than those reported for other regions (e.g., upper Midwest and Northeast U.S.). Mercury exceeded the risk threshold for belted kingfisher (Ceryle alcyon, the most sensitive species considered) in 33-75% of river length and 1-7% of river length for humans. Concentrations were lower in the Missouri than in the Mississippi and Ohio Rivers, consistent with continental-scale patterns in atmospheric Hg deposition. Body size and trophic guild were the best predictors of Hg concentrations, which were highest in large-bodied top predators. Site geochemical and landscape properties were weakly related with fish Hg. Moreover, relationships often ran contrary to conventional wisdom, and the slopes of the relationships (positive or negative) were inconsistent among fish guilds and rivers. For example, sulfate is positively associated with fish Hg concentrations but was negatively correlated with Hg in five of six regression models of tissue concentrations. Variables such as pH, acid neutralizing capacity, and total phosphorus did not occur at levels associated with high fish Hg concentrations, partially explaining the relatively low Hg values we observed.

Persistent organic pollutants in fish tissue in the mid-continental great rivers of the United States.

Great rivers of the central United States (Upper Mississippi, Missouri, and Ohio rivers) are valuable economic and cultural resources, yet until recently their ecological condition has not been well quantified. In 2004-2005, as part of the Environmental Monitoring and Assessment Program for Great River Ecosystems (EMAP-GRE), we measured legacy organochlorines (OCs) (pesticides and polychlorinated biphenyls, PCBs) and emerging compounds (polybrominated diphenyl ethers, PBDEs) in whole fish to estimate human and wildlife exposure risks from fish consumption. PCBs, PBDEs, chlordane, dieldrin and dichlorodiphenyltrichloroethane (DDT) were detected in most samples across all rivers, and hexachlorobenzene was detected in most Ohio River samples. Concentrations were highest in the Ohio River, followed by the Mississippi and Missouri Rivers, respectively. Dieldrin and PCBs posed the greatest risk to humans. Their concentrations exceeded human screening values for cancer risk in 27-54% and 16-98% of river km, respectively. Chlordane exceeded wildlife risk values for kingfisher in 11-96% of river km. PBDE concentrations were highest in large fish in the Missouri and Ohio Rivers (mean>1000 ng g(-1) lipid), with congener 47 most prevalent. OC and PBDE concentrations were positively related to fish size, lipid content, trophic guild, and proximity to urban areas. Contamination of fishes by OCs is widespread among great rivers, although exposure risks appear to be more localized and limited in scope. As an indicator of ecological condition, fish tissue contamination contributes to the overall assessment of great river ecosystems in the U.S.

Comparison of macroinvertebrate sampling methods for nonwadeable streams.

Bioassessment of nonwadeable streams in the United States is increasing, but methods for these systems are not as well-developed as for wadeable streams. In this study, we compared six macroinvertebrate field sampling methods for nonwadeable streams adapted from those used by three major programs: the U.S. Environmental Protection Agency's Environmental Monitoring and Assessment Program-Surface Waters, the U.S. Geological Survey's National Water Quality Assessment Program, and the Ohio Environmental Protection Agency, Division of Surface Water Biocriteria Program. We performed all six methods at 60 sites across four rivers and measured water chemistry and physical habitat at each site to assess abiotic conditon. Sites were divided into two groups: those influenced by navigational lock and dam structures (restricted flow, or RF) and those free-flowing or with lowhead dams (run-of-the-river, or ROR). Metrics based on passive Hester-Dendy artificial substrate samplers differed greatly from active sampling methods (i.e., using nets) but represented abiotic conditions well in both ROR and RF sites. Although metric values were similar across certain sampling methods, the metrics significantly correlated with abiotic variables varied among methods and between ROR and RF sites. These results emphasize that methods are not interchangeable, and the ability to detect certain stressors depends on sampling method.

Electrofishing in boatable rivers: does sampling design affect bioassessment metrics?

Data were collected from 60 boatable sites using an electrofishing design that permitted comparisons of the effects of designs and distances on fish assemblage metrics. Sites were classified a priori as Run-of-the-River (ROR) or Restricted Flow (RF). Data representing four different design options (i.e., 1000 and 2000 m for both single and paired banks) were extracted from the dataset and analyzed. Friedman tests comparing metric values among the designs detected significant differences for all richness metrics at both types of sites and for catch per unit effort and percent tolerant species at ROR sites. Richness metrics were generally higher for the two 2000-m designs than for the two 1000-m designs. When plotted against cumulative electrofishing distance, the percent change in metrics declined sharply within approximately 1000 m, after which metrics usually varied by less than 10%. These data demonstrate that designs electrofishing 1000 m of shoreline are sufficient for bioassessments on boatable rivers similar to those in this study, regardless of whether the shoreline is along a single bank or distributed equally among paired banks. However, at sites with depths greater than 4 m, it may be advisable to employ nighttime electrofishing or increase day electrofishing designs to 2000 m.

A performance comparison of metric scoring methods for a multimetric index for Mid-Atlantic highlands streams.

When biological metrics are combined into a multimetric index for bioassessment purposes, individual metrics must be scored as unitless numbers to be combined into a single index value. Among different multimetric indices, methods of scoring metrics may vary widely in the type of scaling used and the way in which metric expectations are established. These differences among scoring methods may influence the performance characteristics of the final index that is created by summing individual metric scores. The Macroinvertebrate Biotic Integrity Index (MBII), a multimetric index, was developed previously for first through third order streams in the Mid-Atlantic highlands of the United States. In this study, six metric scoring methods were evaluated for the MBII using measures related to site condition and index variability, including the degree of overlap between impaired and reference distributions, relationships to a stressor gradient, within-sample index variability, temporal variability, and the minimum detectable difference. Measures of index variability were affected to a greater degree than those of index responsiveness by both the type of scaling (discrete or continuous) and the method of setting expectations. A scoring method using continuous scaling and setting metric expectations using the 95th percentile of the entire distribution of sites performed the best overall for the MBII. These results showed that the method of scoring metrics affects the properties of the final index, particularly variability, and should be examined in developing a multimetric index because these properties can affect the number of condition classes (e.g., unimpaired, impaired) an index can distinguish.

Development and evaluation of the Lake Macroinvertebrate Integrity Index (LMII) for New Jersey lakes and reservoirs.

In response to the recent focus by the U.S. Environmental Protection Agency on bioassessment of lakes, a multimetric index was developed for New Jersey lakes and reservoirs using benthic macroinvertebrates. Benthic samples were collected from reference and impaired lakes with muck and intermediate sediments in central and northern New Jersey during summer 1997. We used a stepwise process to evaluate properties of candidate metrics and selected five for the Lake Macroinvertebrate Integrity Index (LMII): Hilsenhoff Biotic Index (HBI), percent chironomid individuals, percent collector-gatherer taxa, percent oligochaetes/leeches, and number of Diptera taxa. We scored metrics as the fraction of the best expected value (based on all sites) achieved at a site and summed them into the LMII. Evaluation of the LMII showed that it discriminated well between reference and impaired lakes and was strongly related to several potential stressors. Chemical and physical gradients distinguished between reference and impaired lakes, and the LMII summarized these gradients well. The LMII corresponded strongly with land use, but some lakes with more urban land use still achieved high scores. Based on a power analysis, the ability of the LMII to detect differences in condition was sensitive to the number of samples from each lake.

Methods development and use of macroinvertebrates as indicators of ecological conditions for streams in the Mid-Atlantic Highlands Region.

The Mid-Atlantic Highlands Assessment (MAHA) included the sampling of macroinvertebrates from 424 wadeable stream sites to determine status and trends, biological conditions, and water quality in first through third order streams in the Mid-Atlantic Highlands Region (MAHR) of the United States in 1993-1995. We identified reference and impaired sites using water chemistry and habitat criteria and evaluated a set of candidate macroinvertebrate metrics using a stepwise process. This process examined several metric characteristics, including ability of metrics to discriminate reference and impaired sites, relative scope of impairment, correlations with chemical and habitat indicators of stream disturbance, redundancy with other metrics, and within-year variability. Metrics that performed well were compared with metrics currently being used by three states in the region: Pennsylvania, Virginia, and West Virginia. Some of the metrics used by these states did not perform well when evaluated using regional data, while other metrics used by all three states in some form, specifically number of taxa, number of EPT taxa, and Hilsenhoff Biotic Index, performed well overall. Reasons for discrepancies between state and regional evaluations of metrics are explored. We also provide a set of metrics that, when used in combination, may provide a useful assessment of stream conditions in the MAHR.

Development and evaluation of a Macroinvertebrate Biotic Integrity Index (MBII) for regionally assessing Mid-Atlantic Highlands Streams.

The Macroinvertebrate Biotic Integrity Index (MBII) was developed from data collected at 574 wadeable stream reaches in the Mid-Atlantic Highlands region (MAHR) by the U.S. Environmental Protection Agency's (USEPA) Environmental Monitoring and Assessment Program (EMAP). Over 100 candidate metrics were evaluated for range, precision, responsiveness to various disturbances, relationship to catchment area, and redundancy. Seven metrics were selected, representing taxa richness (Ephemeroptera richness, Plecoptera richness, Trichoptera richness), assemblage composition (percent non-insect individuals, percent 5 dominant taxa), pollution tolerance [Macroinvertebrate Tolerance Index (MTI)], and one functional feeding group (collector-filterer richness). We scored metrics and summed them, then ranked the resulting index through use of independently evaluated reference stream reaches. Although sites were classified into lowland and upland ecoregional groups, we did not need to develop separate scoring criteria for each ecoregional group. We were able to use the same metrics for pool and riffle composite samples, but we had to score them differently. Using the EMAP probability design, we inferred the results, with known confidence bounds, to the 167,797 kilometers of wadeable streams in the Mid-Atlantic Highlands. We classified 17% of the target stream length in the MAHR as good, 57% as fair, and 26% as poor. Pool-dominated reaches were relatively rare in the MAHR, and the usefulness of the MBII was more difficult to assess in these reaches. The process used for developing the MBII is widely applicable and resulted in an index effective in evaluating region-wide conditions and distinguishing good and impaired reaches among both upland and lowland streams dominated by riffle habitat.