We developed and tested a species-level Bray–Curtis (BC) similarity to reference index and a genus-level diatom model affinity (DMA) index to quantify agricultural and acid mine drainage (AMD) impacts on streams in the Western Allegheny Plateau of southeastern Ohio. Decreased similarity to references sites was found in impaired streams, and diatom metrics further indicated how assemblages were impacted. Sites identified by index scores as impaired had significantly greater conductivity, PO4-P, Ca2 , Mg2 , Na , Cl−, % pasture in the upstream watershed, and forest fragmentation than minimally impaired sites (p < 0.05). Percent acidophilic diatoms significantly increased with reduced alkalinity caused by mining impacts (p < 0.01). Relative abundances of species indicating low P and N status decreased with increased pasture and row crops in upstream watersheds (p < 0.01), whereas abundance of motile diatoms increased (p < 0.01). In a case study, diatom, macroinvertebrate, and fish assemblages were compared among 18 sites in a single watershed (6 sites along an AMD impact gradient, 6 NaOH-treated AMD sites, and 6 sites with no mining impacts). All assemblages indicated severe impairment by AMD, but fish were least useful because they were absent from 5 of the 6 sites. Macroinvertebrates indicated unimpaired conditions at treated sites, but fish signaled potential problems (lower than expected biomass and presence of deformities) despite high species richness and index of biotic integrity scores. Diatom indices indicated significant impairment at AMD and treated sites (p < 0.01). Diatom metrics and indices very effectively signaled agricultural and mining impacts at the regional scale and were especially useful in the case study because they provided finer resolution of AMD effects and additional information of ecological importance missed by or in conflict between macroinvertebrate and fish assemblages. The DMA and BC indices responded similarly to stressors at the regional scale, but BC more effectively signaled impacts at the watershed scale. However, DMA was an effective assessment tool that might be used more easily than BC by novice phycologists and watershed groups because it requires less taxonomic expertise. If implemented in future watershed or regional studies, diatom indices could greatly benefit policy, management plans, and current monitoring efforts.
Human activities, such as agriculture, urbanization, and coal mining, threaten or impair thousands of stream kilometers throughout Appalachian states (Virginia, West Virginia, Ohio, and Pennsylvania) (USEPA 2000) and around the world (Grimm et al. 2008). Agriculture and coal mining affect streams and biological assemblages in the Western Allegheny Plateau (WAP) Omernik level III ecoregion of southeastern Ohio by contributing to sedimentation, nutrient loading, and acidic pH with high concentrations of metals (Omernik 1987, OEPA 2000). Biological communities damaged by these stressors lose their resiliency, ability to handle anthropogenic pressures, and ultimately, ecosystem functions, and biodiversity (Baron et al. 2002). Assessments of stream biology can characterize and quantify impacts and are needed to develop management plans, policy, or discharge permits.
In the WAP, acid mine drainage (AMD) has acidic pH and high concentrations of metals resulting from oxidation of FeS in coal mining waste products. AMD impairs diversity (Verb and Vis 2000), productivity (Niyogi et al. 2002), and ecosystem functions, such as metabolism and organic matter breakdown (Niyogi et al. 2001), in streams. AMD has devastated >7200 stream kilometers throughout Appalachia (USEPA 2000). In extremely acidic waters, except for a few macroinvertebrates (Monteith and Evans 2005), algae might be the only biomonitoring organisms available to complement chemical measurements. Thus, the algal community is a valuable and potentially informative tool for assessing AMD impacts.
Depending on land use and conservation practices, agriculture can increase nutrient and sediment loads to streams (Jones et al. 2001), leading to nuisance algal growths (Dodds et al. 1997, Biggs 2000) and loss of habitat heterogeneity (Roth et al. 1996), diversity, and ecological integrity (Allan 2004). Organic matter processing, stream metabolism, and other ecosystem services also can be negatively affected by agricultural activities (Young et al. 2008). Monitoring methods must provide accurate distinction among and quantification of stressors to be effective tools for restoring, managing, and conserving stream ecological integrity (Karr and Yoder 2004).
Diatoms are especially useful for characterizing impacts because they are predominant primary producers in streams and have a wide range of ecological tolerances (Patrick 1973, Lowe 1974, van Dam et al. 1994). Indices of biotic integrity (IBIs), which assess structural and functional components of assemblages, are effective tools for assessing anthropogenic impacts and are used to inform restoration and conservation policy (Karr 1991, Karr and Yoder 2004). Fish and macroinvertebrate indices are used in the WAP region of Ohio, but diatom indices have not been developed or tested on a regional scale (OEPA 2000). Different assemblages respond to different stressors that operate at multiple scales and, thus, can give conflicting diagnoses. Therefore, multiple lines of evidence from various taxonomic groups are important in stream assessments (Johnson and Hering 2009). Diatoms might be more responsive than other assemblages to certain stressors, and diatom-based monitoring tools that can be used to assess regional stressors are needed to complement or supplement results from other assemblages and help refine management efforts.
Indices and metrics based on biological data from least-impaired reference sites in a region help set realistic restoration and conservation goals and improve our understanding of how pollution and land use affect ecological integrity (Karr and Chu 2000). Least-impaired reference sites reflect how biological communities would exist, function, and persist with minimal anthropogenic impacts and are useful benchmarks from which to judge impairment (Stoddard et al. 2006). Diatom assemblages at least-impaired reference sites converge on a similar structure at the regional scale, and a test site's loss of similarity to reference sites can indicate anthropogenic stressors (Passy and Bode 2004, Wang et al. 2005). Changes in specific metrics, such as representation of taxonomic groups, nutrient optima, and morphologies, can be used to identify specific stressors that might lead to the overall loss of resemblance (Barbour et al. 1999). Metrics and similarity measures are standardized and can be used to communicate stressor effects in streams throughout the country (Barbour et al. 1999).
Our goal was to develop diatom monitoring tools for agricultural and AMD impacts. We created 2 indices that scored sites based on their mean similarity to reference-site diatom assemblages: 1) a species-level index based on mean Bray–Curtis (BC) similarity, and 2) a genus-level index based on % similarity of taxonomic groups. A genus-level index, if effective, might be more useful than a species-level index to novice phycologists and watershed groups wanting to integrate diatoms into their biological monitoring programs. We documented the effects of AMD and agriculture on the ecological similarity of study-site diatom assemblages to reference-site assemblages and correlations of changes with individual metrics. We developed and tested the indices/metrics at the regional scale and further tested them in a watershed-scale case study. We also compared the performances of indices based on diatoms, fish, and macroinvertebrates.
The geology of the WAP is predominantly sandstone, shale, limestone, and coal (Omernik 1987). We sampled 60 regional sites (10 reference sites, 50 test sites) throughout the WAP during summer 2005 and 2006 (Fig. 1). We chose regional reference and test sites randomly from an Ohio Environmental Protection Agency (OEPA) database of all sites historically sampled for fish, macroinvertebrates, or chemistry. From these 60 sites, we selected 7 reference sites that currently meet OEPA nutrient (PO4-P < 0.06 mg/L, NO3-N < 0.47 mg/L), ionic composition (<500 µS/cm, Ca2+ < 69, Mg2+ < 18, Na+ < 17, Cl− < 30 mg/L), and physical habitat criteria (OEPA 1999). We added the criteria that reference sites have >66% forest in the upstream watershed, no mining activity, and no point sources of pollution. We selected 3 additional reference sites from Omernik level IV subecoregions (Pittsburgh Low Plateau, Unglaciated Upper Muskingum Basin, and Monongahela Transition Zone; Omernik 1987). We included these subecoregions to increase the ability of our indices to represent least-impaired conditions in the WAP. However, these subecoregions have more human disturbance and agriculture and less forested area than other subecoregions in the WAP, and no site in these subecoregions met all reference criteria. In each subecoregion, we chose the site that failed the fewest criteria (each site marginally failed 1 requirement). Fifteen percent (10/60) of the sites sampled were reference sites. Reference sites were on 2nd- to 5th-order wadeable streams with watersheds that ranged from 6.8 to 248 km2. The 50 test sites were variously affected by agriculture, forest fragmentation, AMD, elevated conductivity, and nutrients.
Leading Creek sites
The Leading Creek watershed drains 388 km2 of steep hills with narrow valley floors (Fig. 1). AMD from coal mining and agriculture heavily affect the watershed, which is ∼67% forested and 26% pasture and hay (Kennedy et al. 2004, Rankin 2005). We sampled 18 sites in the Leading Creek watershed during August 2006. Six sites represented an AMD impact gradient (AMD sites), 6 were affected by AMD discharge treated with NaOH (treated sites), and 6 had no mining impacts but some agricultural influences (non-AMD sites). Treated sites were 1.12 to 6.92 km from the next nearest site.
Sampling and laboratory techniques
We used the National Land Cover Database (USGS 1992) to determine % forest, % pasture, and % row crops in the watershed upstream of each sampling site. We determined connectivity, a measure of forest fragmentation, from US Forest Service data (USFS 2005) at 5 spatial scales (0.023, 0.073, 0.656, 5.905, and 53.144 km2). We calculated connectivity as the probability that 2 adjacent 30 × 30 m pixels are forested, given that the first is forested. We report only the values calculated at the largest spatial scale because they were most strongly associated with indices, metrics, and water chemistry.
We collected samples for water chemistry from each site and filtered them on location through 0.45-µm-pore size filters (Millipore®, Billerica, Massachusetts). We measured PO4-P, NO3-N, SO42−, Cl−, alkalinity, and acidity with US Environmental Protection Agency (EPA)-approved protocols (APHA 1995). We analyzed Ca2+, Mg2+, K+, Na+, Fe2+, Mn2+, and Al3+ with an inductively coupled plasma atomic emission spectrometer (Varian Vista, Palo Alta, California). We measured conductivity, pH, and temperature in situ with handheld probes (Waterproof ECTestr® and pHTestr 30®; Oakton, Vernon Hills, Illinois).
At each stream, we proportionately sampled the available benthic habitats (riffles and pools) and substrata (rock, cobble, silt, and sand) in a 50-m reach to produce a composite sample of periphyton from 10 sampling locations (Stevenson and Bahls 1999). We sampled rocks by scraping them with a firm-bristled toothbrush and soft sediments by suction with a pipette and a 7.1-cm2 O-ring. We preserved the composite samples with 2.5% CaCO3-buffered glutaraldehyde. We cleaned a 10-mL subsample with 30% H2O2 and 50% HNO3 to remove organic matter and allowed the material to settle on round cover slips in a chamber to disperse diatoms randomly and evenly (Battarbee 1973). We mounted cover slips on slides with NAPHRAXTM (The Biology Shop, Hazelbrook, New South Wales, Australia) and counted diatoms at 1000× with a light microscope (Olympus BX40TM, Center Valley, Pennsylvania). We used the keys of Krammer and Lange-Bertalot (1986, 1988, 1991a, b) to identify 500 valves per sample along a transect. We calculated relative abundances of species, Shannon index of diversity, and species richness for each sample.
Diatom indices and metrics
We used the predominant genera found at reference sites to create a model diatom assemblage (diatom model affinity [DMA] index; Passy and Bode 2004). We defined predominant genera as those with >1% relative abundance in >1 reference site. Recent taxonomic splitting of genera (e.g., Round et al. 1990, Lange-Bertalot 2001) was not used because congeneric diatoms frequently represent similar conditions (van Dam et al. 1994, Hill et al. 2001). For this index, Achnanthidium included Planothidium species; Cymbella included Encyonema, Encyonopsis, and Reimeria species; Nitzschia included Tryblionella species; Navicula included Geissleria, Hippodonta, Luticola, Placoneis, and Sellophora species; and Fragilaria included Synedra, Staurosirella, and Staurosira species. We grouped genera with similar ecological attributes to increase robustness and ease of model use. We grouped Achnanthidium and Cymbella because they are typical of waters with high dissolved O2 and low trophic states (van Dam et al. 1994, Wang et al. 2005). We grouped Amphora and Cocconeis because they are typical of high pH, saprobity, and nutrients, whereas we grouped Fragilaria and Gomphonema because the species collected represented high nutrient conditions (van Dam et al. 1994).
We used the mean of each taxonomic group to set the maximum attainable score for each category of the index within 1 standard deviation (SD) as: Achnanthidium + Cymbella: 50%, Nitzschia: 20%, Navicula: 15%, Amphora + Cocconeis: 10%, and Fragilaria + Gomphonema: 5% (Table 1). We set the maximum attainable score for each taxonomic group close to its respective mean, except for Achnanthidium + Cymbella, which had a large SD because the 3 subecoregion reference sites had lower-than-expected abundances of these genera (17–20%). The mean relative abundance of Achnanthidium + Cymbella for the 7 reference sites that met all reference criteria was 44%. We set the value for this group at 50%, which was within 1 SD of the mean of all 10 reference sites, because these genera typically represent less-impaired conditions (KYDOW 1993, Wang et al. 2005). Scoring sites in this manner set the index on a scale of 0–100 for ease of interpretation.
We scored sites based on % similarity to the DMA index. High similarity indicated less impaired conditions, and low similarity indicated more impaired conditions. Percent similarity was calculated as (Whittaker and Fairbanks 1958):Table 1, DMAMAX for Nitzschia = 20%) and r is the relative abundance of the same taxon at a study site. Scores were determined by summing the smaller of m or r for each taxonomic group. For example, if a study site had 20% Achnanthidium + Cymbella, 30% Nitzschia, 5% Navicula, 8% Amphora + Cocconeis, and 37% Fragilaria + Gomphonema, it would be scored as 20 + 20 + 5 + 8 + 5 = 58% similarity. Thus, index scores are reduced if test sites have relative abundances greater or less than DMA values.
BC similarity is a resemblance coefficient that often is used to analyze species composition patterns with abundance data (Bray and Curtis 1957, Legendre and Legendre 1998). We used a species-based BC similarity index to calculate the similarity between a study site and each reference site as:Legendre and Legendre 1998). We calculated scores as the mean BC similarity between assemblages at study sites and at the reference sites.
We used correlations of water chemistry and landuse variables associated with human impacts with 15 diatom assemblage metrics to assess how specific components of diatom communities contributed to loss of resemblance between test- and reference-site assemblages in response to stressors. We tested metrics associated with taxonomic composition (species richness, Shannon diversity, % dominant diatom species), morphology (prostrate, erect, stalked, unattached, motile, very motile), and autecology (trophic state, N and P optima, pH optima, % eutraphentic, % acidophilic). Morphology reflects the physical structure of a diatom assemblage (Hill et al. 2000, Fore and Grafe 2002, Wang et al. 2005). In particular, motile diatoms are useful for assessing siltation (Stevenson and Bahls 1999), agricultural, and anthropogenic impacts (Fore and Grafe 2002). Autecology reflects general ecological conditions and specific stressors. We used diatom species indicators of high and low nutrient (P and N) status developed for the Eastern Highlands (synonymous with Appalachian Forests located in the Temperate Eastern Forests, Omernik level II ecoregion), which includes the WAP, because they were developed and tested specifically for the region and were more powerful indicators of trophic state than were indicators developed from European studies (Potapova and Charles 2007). We also tested a % eutraphentic diatom metric (Hill et al. 2000) that uses genus attributes and relative abundance of acidophilic diatoms, because of the frequent occurrence of AMD in the region. We report only metrics that were significantly correlated with anthropogenic stressors.
We used NCSS2004 to conduct exploratory univariate and multivariate statistical tests (Number Cruncher Statistical Systems, Kaysville, Utah). We screened data for normality, and either log(x)- or √(x)-transformed data to correct skewness. We used principal components analysis (PCA) of a correlation matrix to summarize the structure of chemical variables in the data set. We used Pearson correlations to indicate relationships among index/metric scores, environmental, and landscape variables.
Canonical correspondence analysis (CCA)
A detrended correspondence analysis (DCA) of √(x)-transformed species relative abundance data had a gradient of 2.3 SD, which indicated a unimodal relationship among most species and environmental variables. Therefore, we conducted CCA with CANOCO version 4.5 (Microcomputer Power, Ithaca, New York; ter Braak and Šmilauer 2002) to explore species–environmental variable relationships. CCA ordinated sites based on relative abundances of diatom species while constraining the solution based on regressions with environmental variables. We eliminated species from CCA if they were not present in >1 site with >1% relative abundance. This step reduced the number of diatom species in the regional data set from 342 to 116. We used a stepwise variable-selection procedure and included only significant (999 Monte Carlo permutations; p < 0.05) and noncollinear variables in the CCA. We used 999 Monte Carlo permutations to test significance of axes. We used Pearson correlations of indices, metrics, common genera, water chemistry, and landuse variables with the first 3 CCA axes to indicate their relationships to stressors and diatom assemblages.
Defining levels of impairment
Narrative categories help communicate results to the public and policy makers and can indicate benchmarks of ecological condition (Herlihy et al. 2008). We used DMA scores to define levels of dissimilarity because changes in genera typically are more conservative than changes in species (Hill et al. 2001). We used the 25th percentile score of reference sites to define the boundary between unimpaired and moderately impaired conditions and the 50th percentile of test sites to define the boundary between moderate and severe impairment. We used this arbitrary categorization of impaired sites only for the purpose of determining if diatom similarity measures could indicate different levels of stressors present at the 60 study sites. The 50th percentile probably would change if sites were added to the data set in the future, and therefore, this score should not be used as a level of attainment status in future studies. However, our use of the 25th percentile score to indicate unimpaired conditions was not arbitrary because it was based on the distribution of values at reference sites, and it is a common criterion used in stream biomonitoring (Herlihy et al. 2008).
We tested for differences in environmental variables associated with human activity and correlated with indices/metrics (Ca2+, Mg2+, Na+, PO4-P, Cl−) among impairment categories (factor variable) with a multivariate analysis of variance (MANOVA; Wilks' Lamda). We used 1-way analyses of variance (ANOVAs) with Tukey–Kramer tests for multiple comparisons to determine significant differences among impairment categories for each variable (Zar 1999). We also used ANOVAs with Tukey–Kramer tests for multiple comparisons to test for differences in metric values among impairment categories.
Leading Creek watershed
We used 18 sites in the Leading Creek watershed as a case study to evaluate the ability of the diatom indices to detect impairment from agriculture and AMD relative to existing fish and macroinvertebrate IBIs. We chose 6 AMD sites, 6 treated sites, and 6 non-AMD sites based on recent studies of water chemistry, macroinvertebrates, and fish (Kennedy et al. 2004, Rankin 2005). We used fish and macroinvertebrate data from Rankin (2005) to calculate species richness, macroinvertebrate impairment scores, and fish IBI scores for each site. Macroinvertebrate scores ranged from 1 to 4, with 4 being the least impaired. Higher fish IBI scores indicated less impaired conditions. We used DCA of √(x)-transformed diatom relative abundances to explore patterns among groups of sites along the main axis of variation. We used a MANOVA (Wilks' Lambda) to test for differences in the response variables of pH, SO42−, and conductivity (chemical variables indicating mining impairment) among groups of sites. We used 1-way ANOVAs followed by Tukey–Kramer tests to determine significant differences in pH, conductivity, SO42−, and assemblage structure (fish, macroinvertebrate, and diatom richness and index scores) among groups of sites. When needed, we transformed variables used in the ANOVAs to meet assumptions of normality and homogeneity of variances.
BC and DMA
A total of 342 taxa were identified from the 60 sites. Shannon diversity ranged from 2.0 to 3.8, and species richness ranged from 44 to 80. DMA scores ranged from 30.7 to 86.4, and mean BC similarity ranged from 10.8 to 50.1% for the 60 sites. BC and DMA scores were strongly correlated (r = 0.75, p < 0.001).
PCA of chemistry data identified conductivity and nutrients as environmental variables contributing most to the diatom data structure (axis 1: 50% of variation, axis 2: 19% of variation). Factor loadings on axes 1 and 2 were greatest for Ca2+ (0.47, 0.46, respectively), Mg2+ (0.48, 0.37), Na+ (0.46, −0.30), Cl− (−0.26, −0.22), PO4-P (0.41, −0.20), and K+ (0.34, −0.69). Percent forested area, forest connectivity, DMA scores, and BC scores were negatively correlated with these water chemistry variables and conductivity (r = −0.29 to −0.65; Table 2). DMA and BC were positively correlated with % forested area (r = 0.42, r = 0.37) and negatively correlated with % pasture (r = −0.34, r = −0.29) and row crops (r = −0.33, r = −0.29) (Table 2). DMA and BC were strongly correlated with forest connectivity at the 53-km2 scale (r = 0.63 and 0.60, p < 0.001; Table 2).
Species richness was correlated with Ca2+ (r = −0.46), Mg2+ (r = −0.32), conductivity (r = −0.37), and forest connectivity (r = 0.29). Shannon diversity was correlated only with Ca2+ (r = −0.35) and pH (r = 0.29). Relative abundances of diatoms indicating low P and low N were positively correlated with forest connectivity (r = 0.48, r = 0.59) and % forest (r = 0.37, r = 0.54) and negatively correlated with % pasture (r = −0.35, r = −0.46, respectively) (Table 2). Relative abundance of diatoms indicating low N also was negatively correlated with % row crops (r = −0.41, p < 0.01). Diatoms indicating high P and N were negatively correlated with forest connectivity (r = −0.36, r = −0.34) and % forest (r = −0.36, r = −0.46, respectively) and positively correlated with % pasture (r = 0.32, r = 0.35). Relative abundance of eutraphentic diatoms was not correlated with any land use, but was correlated with Ca2+, conductivity, and alkalinity (Table 2). Acidophilic diatom relative abundances were negatively correlated with pH and alkalinity (r = −0.44, r = −0.68). Percent motile diatoms was the only morphological metric correlated with measures of human disturbance and was negatively correlated with forest connectivity (r = −0.41) and % forest (r = −0.46) and positively correlated with % pasture (r = 0.51) and % row crops (r = 0.35).
Diatom community analysis
Alkalinity, PO4-P, Ca2+, Mg2+, Na+, pH, SO42−, and watershed area (p < 0.05) were kept in the stepwise CCA with diatom species and habitat data (Fig. 2). These variables explained 31.2% of the total variation in the species data (axis 1: 12%, axis 2: 5.7%, axis 3: 4.4%, all other axes: 9.1%). BC, DMA, % forest, % pasture, and % row crops were used as passive variables to display their relationships with environmental gradients on the CCA biplot (Fig. 2). All 10 reference sites were located in the lower left quadrant of the CCA biplot. Measures of taxonomic composition, genera, index scores, metrics, land use, and water chemistry were significantly correlated with axes 1, 2, and 3 (Table 3).
Axis 1 represented an alkalinity gradient. Sites with positive axis-1 scores were impacted by AMD, indicated by strong correlations of alkalinity, pH, cations (associated with leaching by H+), and conductivity (Table 3). Percent acidophilic diatoms also was strongly correlated with axis 1 scores (r = 0.77, p < 0.01). Relative abundances of Cocconeis and Amphora species, which are typical of alkaline conditions, and Cymbella species were negatively correlated with axis 1 scores (r = −0.59, r = −0.63, r = −0.37, respectively).
Axis 2 represented an agricultural impact gradient and was significantly correlated with % forest, forest connectivity, % pasture, % row crops, diatom metrics indicating high and low P or N concentrations, and % motile diatoms (Table 3). Sites with greater relative abundances of diatoms indicative of high P and high N conditions had more positive axis 2 scores (r = 0.37, r = 0.36, respectively), whereas sites with greater diatom abundances indicative of low P and low N conditions had more negative axis 2 scores (r = −0.63, r = −0.68). Relative abundances of Achnanthidium and Cymbella were negatively correlated with axis 2 scores (r = −0.73, r = −0.37), and were associated with lower % pasture and higher forest connectivity. Relative abundance of eutraphentic diatoms was positively correlated with axis 2 scores (r = 0.36), as were relative abundances of Navicula, Nitzschia, and Cocconeis (r = 0.40, r = 0.28, r = 0.41). Percent pasture and % row crops were positively correlated with PO4-P, Ca2+, Mg2+, and Na+ (Fig. 2). Sites with high similarity (BC and DMA) to reference sites (i.e., less impairment) had more negative CCA axis 1 and 2 scores (Table 3).
Twenty-three sites had scores ≥25th percentile of DMA reference scores (scores > 67) and were categorized as unimpaired. The remaining 37 sites were categorized as moderately impaired (scores = 55–66) or severely impaired (scores < 55) as set by the 50th percentile of impaired sites in this study. BC scores differed significantly among sites in different impairment categories (ANOVA, p < 0.01). BC scores were highest for unimpaired sites, intermediate for moderately impaired sites, and lowest for severely impaired sites (p < 0.01). Ca2+, Mg2+, PO4-P, Na+, and Cl−, important indicators of anthropogenic stress, differed significantly among sites in different impairment categories (MANOVA, F10,104 = 4.57, p < 0.001). Conductivity, Ca2+, Mg2+, Na+, Cl− and PO4-P were significantly higher at moderately and severely impaired sites than at unimpaired sites (ANOVA; Fig. 3A–F). Unimpaired sites had higher % forest, lower % pasture (Fig. 3G, H), and higher forest connectivity than did moderately or severely impaired sites (p < 0.01). Relative abundances of diatoms indicative of low P and N conditions were highest and relative abundances of motile and eutraphentic diatoms were lowest at unimpaired sites (ANOVA; Fig. 4A–D). Species richness and Shannon diversity tended to be higher at moderately impaired sites than at unimpaired sites, but were significantly lowest at severely impaired sites (ANOVA, p < 0.01; Fig. 4E, F).
A total of 186 diatom taxa were identified from the 18 sites. Shannon diversity ranged from 1.56 to 3.90, and species richness ranged from 12 to 85. The 3 groups of sites differed distinctly in community composition (Fig. 5). Non-AMD sites had scores in the middle of DCA axis 1. Treated sites were on the far left and AMD sites were on the far right of axis 1.
Shannon diversity was significantly greater at non-AMD sites (mean ± SD; 3.26 ± 0.48) and treated sites (3.00 ± 0.27) than at AMD sites (2.13 ± 0.44). Diatom assemblages at severely impaired AMD sites (pH < 4) were dominated by Eunotia exigua (Brébisson ex Kützing) Rabenhorst and Frustulia rhomboides v. saxonica (Rabenhorst) De Toni. At pH > 5.8, Achnanthidium minutissimum (Kützing) Czarnecki was dominant at AMD sites despite high SO42− (>158 mg/L) and conductivity (980 µS/cm). Conductivity at treated sites ranged from 1550 to 9940 µS/cm, and diatom assemblages were dominated by brackish-water species, such as Navicula recens (Lange-Bertalot) Lange-Bertalot, Pleurosira laevis (Ehrenberg) Compère, Nitzschia reversa W. Smith, and Thalassiosira lacustris (Grunow) Hasle.
Conductivity, pH, and SO4 (important indicators of AMD) differed among sites in the 3 impairment categories (MANOVA, F6,24 = 15.54, p < 0.001). pH at treated sites was comparable to pH at non-AMD sites (Fig. 6A), but SO42− and conductivity were significantly higher at treated than at non-AMD sites (ANOVA, p < 0.001; Fig. 6B, C). Conductivity was significantly greater at treated sites than at AMD sites because NaOH treatments increased the already high conductivity associated with AMD. AMD sites had significantly lower pH than did sites in the other 2 groups (p < 0.001), higher conductivity and SO42− than non-AMD sites (p < 0.001), and similar SO42− concentrations to treated sites (Fig. 6A–C).
DMA scores and BC % similarity at sites were strongly correlated (r = 0.89, p < 0.001). With the exception of 1 outlier site (score = 39) with an assemblage dominated by Navicula perminuta Grunow, non-AMD sites had DMA scores that ranged from 58 to 79 (Fig. 6D). This site clustered with other non-AMD sites in the DCA, but was separated from the cluster along axis 2 (Fig. 5). Elevated nutrient concentrations (NO3-N > 1 mg/L, PO4-P > 0.25 mg/L) at non-AMD sites with low DMA and BC scores indicated that these sites were affected by agriculture. DMA scores reflected the gradient of AMD impairment, but because of the high abundance of A. minutissimum in moderately impacted streams (pH > 5.8), scores were higher than expected (>65) (Fig. 6D). DMA scores ranged from 17 to 39 at AMD streams with pH < 5.8 and from 38 to 51 at treated sites. These scores indicated severe impairment based on the regional data set. BC % similarity at non-AMD sites ranged from 35 to 42%, with the exception of the outlier site (18.8%), and were significantly greater than BC % similarity at AMD and treated sites (p < 0.01; Fig. 6E). BC % similarity ranged from 1.4 to 27.7% at AMD sites and from 14.6 to 26.2% at treated sites. Diatom richness was significantly higher at non-AMD sites than at AMD sites (Fig. 6F).
Macroinvertebrate impairment scores and species richness and fish IBI scores were significantly lower at AMD sites than at other sites and indicated impairment, but these metrics were similar between treated and non-AMD sites (Fig. 6G–I). Fish species richness was significantly different among all 3 groups (p < 0.01; Fig. 6J). Fish were absent from 5 of the 6 AMD impacted sites, and fish species richness was greatest at treated sites, intermediate at non-AMD sites, and lowest at AMD sites. Fish with deformities were present in 4 of the 6 treated sites, but deformed fish were not observed in any other sites. Fish biomass was lower than expected at treated sites when compared to similar regional reference sites (E. T. Rankin, Environmental Program Manager, Institute for Local Government Administration and Rural Development, Ohio University, personal communication).
In general, all 3 assemblages indicated severe impairment at AMD sites, but macroinvertebrates and fish did not detect clear differences between treated and non-AMD sites (Fig. 6G–J). Fish data indicated potential problems (low biomass, deformities) at treated sites, but these results were in conflict with high species richness and IBI scores. The diatom assemblage was the only one to signal impaired conditions clearly via assemblage structure (Fig. 5) and loss of similarity to regional reference streams (Fig. 6D–F). Fish and macroinvertebrates showed clear improvements in treated sites compared to AMD sites, a result that indicated that NaOH treatment improved conditions, even though other stressors associated with the discharge remained ecologically problematic.
BC, DMA, and diatom metrics very effectively represented components of the diatom community that were useful for assessing, interpreting, and communicating how anthropogenic activities impacted streams. Human impacts on landscapes strongly influence habitat conditions and biological communities of streams (Allan 2004). The indices and metrics developed and tested in our study signaled gradients of human impacts, particularly AMD, agriculture, and forest fragmentation. Water chemistry variables, such as conductivity, cation concentrations, and PO4-P, were associated with landuse patterns and further linked diatom community responses with human impacts. Sites with index and metric scores indicating moderate or severe impairment were in areas with low forestation, increased agricultural and mining activities, and highly fragmented forests (Figs 2–Figure 3.4). Loss of cations and P from disturbed watersheds, especially as a result of leaching from agricultural (Thornton and Dise 1998, Rhodes et al. 2001) or deforested lands (Swank et al. 2001), can increase stress on aquatic communities.
Ionic composition of stream water is extremely important, and Ca2+, alkalinity, pH, and HCO3− strongly influence diatom communities (Biggs 1995, Potapova and Charles 2003, Soininen 2007). Agriculture can elevate conductivity and cation concentrations as a result of fertilizing, liming, animal wastes, and erosion or particulate deposition to streams (Johnson et al. 1997, Carpenter and Waite 2000, Dow and Zampella 2000). AMD also can increase leaching of cations because of elevated H+ concentrations. Percent acidophilic diatoms was an excellent indicator of AMD with low alkalinity but occasionally high cation concentrations, indicated by strong significant correlations with alkalinity and CCA axis 1 (Fig. 2). Impaired sites had significantly greater Ca2+, Mg2+, and Na+ concentrations than did unimpaired sites (Fig. 3A–C) because of AMD impacts at some sites (CCA axis 1) and agricultural impacts at others (CCA Axis 2) (Fig. 2). The ratio of the sum of base cations to alkalinity indicates mineral weathering (1:1) or probable anthropogenic sources of cations (>1:1) (Rhodes et al. 2001). Most severely and moderately impacted sites had ratios substantially >1:1, whereas this ratio rarely was exceeded (and then only slightly) at unimpaired sites. The lack of correlation between Ca2+ and alkalinity (r = 0.19, p = 0.14) (Fig. 2) further indicated that both agricultural land use and AMD contribute to increased cations because Ca2+ and Mg2+ concentrations would be correlated with Ca- or Mg-bicarbonates if natural chemical weathering were the main source of cations (Leland and Porter 2000).
Diatom assemblages were exceptional indicators of AMD and land use in the watershed and adjacent to streams because they were strongly structured by the chemical habitat, which is a result of the water source, percolation through the soils, and runoff in the upstream watershed. Forest fragmentation and deforestation can greatly increase leaching and runoff of PO4-P and cations into streams and can increase sedimentation, especially if the fragmentation is caused by clearing of land for pasture or crops (Owens et al. 2003, Kayser and Isselstein 2005). In our study, diatom indices were more strongly correlated with forest connectivity than with % forest. Thus, processes mediating nutrient inputs to streams might be affected more by conditions in the area immediately adjacent to sites than by processes at the scale of the watershed, and nutrient input might be high even if the watershed is highly forested in other areas away from the stream. Impaired sites had significantly greater cation, PO4-P, and Cl− concentrations, which were more strongly correlated with forest connectivity than % forest in the upstream watershed (Table 2), than did unimpaired sites. Diatoms indicating low P and N status were much more abundant in watersheds with high % forest and connectivity than in watersheds with low % forest or low connectivity (Figs 3, 4, Tables 2, 3). Diatoms indicating low P and N were significantly more abundant at unimpaired than impaired sites (Fig. 4) and were more strongly correlated with chemistry and landuse variables than were high P and N diatoms (Tables 2, 3). Together, these results suggest that sensitive species might be more responsive than tolerant species to increased nutrients. Genus-level % eutraphentic diatoms was less informative than species-level indicators of P or N status as a metric indicating nutrient status (Tables 2, 3). Motile diatoms increased as % pasture increased and % forest decreased, probably because of increased erosion and sedimentation in streams affected by agricultural land use. Motile diatoms also increased along CCA axis 1, which represented an AMD gradient, probably because metal oxide precipitates can smother benthic habitats.
Water chemistry, BC, DMA, and diatom assemblages indicated AMD and NaOH treated AMD impacts (Figs 5, 6). Severely impaired AMD sites had low BC and DMA scores and low species diversity dominated by acidophilic genera, such as Eunotia and Frustulia. Fish were less informative of the AMD gradient because they were absent from 5 of the 6 sites. More macroinvertebrate species (6–24 species) than fish species were present along the AMD gradient, but diatoms were most species rich (12–59) and, therefore, have potential to be informative of AMD impact gradients.
NaOH-treated AMD sites had extremely high conductivities, which were reflected in diatom assemblages characterized by high abundances of brackish-water diatoms. Nonnative diatoms also thrive in these sites (Smucker et al. 2008). Assemblages at these sites have experienced chronic Cl−, Na+, and SO42− toxicity from the discharge (Kennedy et al. 2003). Treated discharge limits fish biomass and growth (E. T. Rankin, personal communication). Deformed fish were found in ⅔ of sites downstream of the discharge. Sensitive aquatic fauna were impaired at sites with >3700 µS/cm in Leading Creek (Kennedy et al. 2003), but the fish IBI and fish richness did not differ between non-AMD and treated AMD sites (Rankin 2005; Fig. 6I, J). Macroinvertebrate data also indicated minimal impairment in these sites (Rankin 2005; Fig. 6D, E). All 6 treated AMD sites were scored as severely impaired by both the DMA and BC indices, but diatom diversity and species richness was not significantly different from non-AMD sites because of the high diversity of brackish-water diatoms. This result underscores the potential problem of using diversity as a metric of human disturbance in some situations.
Diatom communities downstream of NaOH-treated AMD discharge were significantly impaired, and understanding this impairment might indicate mechanisms by which higher trophic levels were impacted. Discharge of NaOH-treated AMD waste is permitted in this watershed, in part because of the conflicting results of macroinvertebrate (no impact) and fish (moderate impact based on the metrics used) monitoring. However, that 3rd-order Leading Creek cannot effectively function and process its effluent is clearly indicated by its diatom assemblage. Including diatoms in the monitoring program might provide the evidence of an ecological problem at these sites that is needed to change the permit.
Species (BC)- vs genus (DMA)-level taxonomic resolution
As with macroinvertebrates, genus-level diatom taxonomy is adequate and informative for some stream assessments (Hill et al. 2001, Wang et al. 2005). Many species within genera have similar growth forms and responses to environmental conditions (Hill et al. 2001). In our study, the DMA and BC indices performed similarly at the regional scale and signaled gradients of AMD, agriculture, and forest fragmentation. Moreover, DMA and BC scores at sites were strongly correlated, probably because numerous genera were represented by a single or few dominant species (e.g., Achnanthidium by A. minutissimum and A. deflexum [C. W. Reimer] J. C. Kingston). These results indicate that DMA could be an effective tool for novice phycologists in biomonitoring programs. However, species-based nutrient metrics more effectively indicated anthropogenic impairments than did the genus-based % eutraphentic diatoms metric.
BC similarity provided finer resolution of environmental impacts than did the DMA index in the watershed case study. Higher-than-expected DMA scores were found in AMD sites because one species of an important DMA genus, A. minutissimum, was able to dominate assemblages in slightly acidic (Verb and Vis 2000) and moderately disturbed streams (Stevenson and Bahls 1999). Moderate AMD impairment is difficult to detect and score properly (Verb and Vis 2000, Hamsher et al. 2004). Moderately impaired sites were scored correctly as impaired with the BC index, whereas 2 moderately impaired AMD sites were scored as unimpaired even though conductivity and SO42− levels were high at these sites. Thus, BC similarity was a more effective index of impairment than DMA.
In conclusion, the DMA and BC indices, along with diatom metrics, effectively diagnosed anthropogenic impairments (agriculture and AMD) of streams. Diatom assemblages were responsive to landscape level and local environmental conditions and can be used to inform management practices and conservation of streams. Diatom monitoring can improve biomonitoring efforts by contributing information missed by or in conflict between other assemblages.
Funding was provided by USEPA Science to Achieve Results grant (R831365). We thank Emily Hollingsworth and Jason Zalack for assistance with sampling and Dina Lopez and Prosper Gbolo for assistance with water chemistry data. We thank James Dyer for compiling GIS data, Michael Hughes for assistance with creating the sampling map, and Ed Rankin for site selection and discussions of this topic. Discussions with Leska Fore and feedback from 2 anonymous reviewers improved the manuscript. Comments from Jared DeForest, Kelly Johnson, and Brian McCarthy are also appreciated.
Formulation of the diatom model affinity (DMA), which is based on genus-level taxonomic composition from 10 reference sites. DMAMAX is the maximum possible score for each category and the total index. Scores for each taxonomic category are calculated as the lesser value of the % counted in the sample or the % in the model (DMAMAX).
Pearson correlations of diatom index scores, metrics associated with nutrients and acid mine drainage, chemical, and landuse variables. DMA = diatom model affinity, BC = Bray–Curtis similarity to reference sites; low P, high P, low N, high N are relative abundances of diatoms indicating low and high nutrient concentrations (Potapova and Charles 2007) for the Eastern Highlands, which includes the Western Allegheny Plateau; EC = electrical conductivity; eutraphentic = relative abundance of diatoms indicating eutrophic conditions (Hill et al. 2000); acid = relative abundance of acidophilic diatoms; motile = relative abundance of motile diatoms; Pff243 = forest connectivity at the 53-km2 scale (see Methods for explanation). Values in bold are statistically significant (p < 0.05), * = p < 0.01, ** = p < 0.001.
Correlations of index scores, taxonomic composition, chemistry, and landuse variables with the first 3 canonical correspondence analysis (CCA) axes. DMA = diatom model affinity, BC = Bray–Curtis similarity to reference sites. Variable names are as in Table 1. Values in bold are statistically significant (p < 0.05).
 1 E-mail addresses: email@example.com