Diatom biomonitoring of streams: Reliability of reference sites and the response of metrics to environmental variations across temporal scales

Benthic diatoms are widely used indicators of human impacts on stream ecosystems because they are very responsive to changing environmental conditions. However, little research has explicitly focused on their reliability with regards to temporal variation in assemblage structure and environmental conditions. We examined variability in diatom-environment relationships at bi-weekly, monthly, and yearly time scales from 7 reference, 7 agricultural, and 2 acid mine drainage (AMD)-impacted streams, and how nutrient and pH fluctuations may affect the interpretation of diatom metrics and the Diatom Model Affinity (DMA) index. Reference streams had less bi-weekly variability in NO₃-N concentrations than non-reference streams. The % eutraphentic diatoms and DMA scores were more strongly correlated with seasonal means of NO₃-N and PO₄-P concentrations than with same day concentrations. Most nutrient indicator metrics had strong correlations with watershed land use. All 14 non-AMD streams experienced substantial increases in NO₃-N and decreases in temperature from November to May, which were associated with high species turnover, substantial changes in community structure, reduced diversity and richness, increased relative abundances of high nutrient diatoms, and decreases in low nutrient diatoms and DMA scores. The % acidophilic diatoms and DMA scores were significantly correlated with increased pH associated with greater precipitation at AMD sites from December to April (r = ?0.77, r = 0.62, respectively; P < 0.01). Yearly, DMA scores for all reference streams were consistently in the minimally impaired category, whereas scores for non-reference streams varied among impairment categories. Reference sites serve as reliable benchmarks for diatom ecological integrity during the summer. In this region, June to October is a recommended time period for diatom sampling in monitoring programs because subsequent shifts in hydrologic regimes, nutrients, and diatom assemblages occurred, affecting all sites and masking among stream differences attributable to agricultural land uses.     

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

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

Diatom taxa and assemblages for establishing nutrient criteria of lakes with anthropogenic hydrologic alteration

Stressor-response models offer guidance for concentration-based nutrient criteria in lakes under human intervention. Diatom-based statistics from biological responses were incorporated to derive taxon-specific and community-level change points (thresholds) of phosphorous and nitrogen in 77 Yangtze floodplain lakes. Diatom metrics relating with conductivity were adopted as response variables, since conductivity explained the maximum variation (38.1%) in diatom assemblages via Bootstrapped regression trees. Nonparametric change-point analysis and Threshold Indicator Taxa ANalysis showed threshold responses of diatom community structure at 0.05–0.08 mg TP/L in connected lakes and 0.02–0.04 mg TP/L in isolated lakes. Distinct community change points of sensitive diatoms occurred at 0.96–1.63 mg TN/L in connected lakes and 0.52–0.63 mg TN/L in isolated lakes. Diatom community structures of tolerant taxa were substantially altered beyond 0.22–0.23 mg/L in connected lakes and 0.52–0.69 mg NO_x/L in isolated lakes. Hydrological river-lake connectivity differed significantly in ecological nutrient criteria with more TN/TP criteria and less NO_x criteria in connected lakes. Given the ecological significance and biological integrity, diatom-based statistics can provide more reliable change points (thresholds) for nutrient criteria than Chl a-nutrient relationships.     

Using algal metrics and biomass to evaluate multiple ways of defining concentration-based nutrient criteria in streams and their ecological relevance

We examined the utility of nutrient criteria derived solely from total phosphorus (TP) concentrations in streams (regression models and percentile distributions) and evaluated their ecological relevance to diatom and algal biomass responses. We used a variety of statistics to characterize ecological responses and to develop concentration-based nutrient criteria (derived from ecological effects) for streams in Connecticut, USA, where urbanization is the primary cause of watershed alteration. Mean background TP concentration in the absence of anthropogenic land cover was predicted to be 0.017 mg/l, which was similar to the 25th percentile of all study sites. Increased TP concentrations were significantly correlated with altered diatom community structure, decreased percent low P diatoms and diatoms sensitive to impervious cover, and increased percent high P diatoms, diatoms that increase with greater impervious cover, and chlorophyll a (P < 0.01). Variance partitioning models showed that shared effects of anthropogenic land cover and chemistry (i.e., chemistry affected by land cover) represented the majority of explained variation in diatom metrics and chlorophyll a. Bootstrapped regression trees, threshold indicator taxa analysis, and boosted regression trees identified TP concentrations at which strong responses of diatom metrics and communities occurred, but these values varied among analyses. When considering ecological responses, scientifically defensible and ecologically relevant TP criteria were identified at (1) 0.020 mg/l for designating highest quality streams and restoration targets, above which sensitive taxa steeply declined, tolerant taxa increased, and community structure changed, (2) 0.040 mg/l, at which community level change points began to occur and sensitive diatoms were greatly reduced, (3) 0.065 mg/l, above which most sensitive diatoms were lost and tolerant diatoms steeply increased to their maxima, and (4) 0.082 mg/l, which appeared to be a saturated threshold, beyond which substantially altered community structure was sustained. These criteria can inform anti-degradation policies for high quality streams, discharge permit decisions, and future strategies for watershed development and managment. Our results indicated that management practices and decisions at the watershed scale will likely be important for improving degraded streams and conserving high quality streams. Results also emphasized the importance of incorporating ecological responses and considering the body of evidence from multiple conceptual approaches and statistical analyses for developing nutrient criteria, because solely relying on one approach could lead to misdirected decisions and resources.     

Establishing expectations for pan-European diatom based ecological status assessments

The European Union (EU) Water Framework Directive depends, for effective implementation, on Member States (MSs) agreeing to a concept of the unimpacted “reference” state, which will then provide the “expected” value in Ecological Quality Ratio (EQR) calculations. Reference assemblages of organism groups will, in turn, vary, due to geological, hydrological, climatic, physicochemical and biological factors. Member States tackle this by establishing “types” which share common characteristics. However, for the purposes of ensuring consistent application, broad transboundary types were also established within five Geographical Intercalibration Groups (GIGs, referred to here as “regions”) as part of the EU's intercalibration exercise. In this paper, we evaluate these types using river diatom assemblages and also provide reference threshold values for two common metrics used in pan-European diatom assessments. A database was assembled, representing 14 EU Member States from Ireland and Portugal in the West, to Estonia and Cyprus in the East, in order to explore biogeographical patterns in assemblages unaffected by anthropogenic pressures. Multivariate analyses were used to examine this pattern and its relationship with geographic, typological and abiotic parameters. After taxonomic harmonisation, NMDS ordination of samples indicated weak differences in assemblage composition among regions. ANOSIM analyses, in turn, indicated that MS was the best factor to group similar samples whereas alkalinity, recognised as the primary environmental variable structuring diatom communities, although significant, explained less variability in the dataset. This, we believe, reflects the importance of methodological factors other than taxonomy (e.g. counting protocol, sample season) that may be constant within a MS but which vary between MSs. When two diatom metrics, the TI and IPS, were applied to the data, differences in the distribution of the metric scores between MS were generally not statistically significant even though some differences between regions were apparent. A trend of increasing values of TI (decreasing values of IPS) was observed in the sequence: Nordic < Alpine < Mediterranean < Central-Baltic < Eastern Continental regions. Additionally, some differences were observed among types within the Mediterranean and Nordic regions, though not for other regions. The data used in this exercise provides us with a region and, in some cases, a type specific benchmark dataset against which national reference data can be compared.     

Diatoms are better indicators of urban stream conditions: A case study in Beijing,China

Urbanization dramatically affects hydrology, water quality and aquatic ecosystem composition. Here we characterized changes in diatom assemblages along an urban-to-rural gradient to assess impacts of urbanization on stream conditions in Beijing, China. Diatoms, water chemistry, and physical variables were measured at 22 urban (6 in upstream and 16 in downstream) and 7 rural reference stream sites during July and August of 2013. One-way ANOVA showed that water physical and chemical variables were significantly different (p < 0.05) between urban downstream and both reference and urban upstream sites, but not between reference and urban upstream sites (p > 0.05). Similarly, structural metrics, including species richness (S), Shannon diversity (H′), species evenness (J′) and Simpson diversity (D′), were significantly different (p < 0.05) between urban downstream and both reference and urban upstream sites, but not (p > 0.05) between reference and urban upstream sites. However, diatom assemblages were very different among all sites. Achnanthidium minutissima was a consistent dominant species in reference sites; Staurosira construens var. venter and Pseudostaurosira brevistriata were the dominant species in urban upstream sites; and Nitzschia palea was the dominant species in urban downstream sites. Clustering analyses based on the relative abundance of diatom species, showed all the samples fit into three groups: reference sites, urban upstream sites, and urban downstream sites. Canonical correspondence analysis (CCA) and Monte Carlo permutation tests showed that concentration of K⁺, EC, TN, Cl⁻ and pH were positively correlated with relative abundance of dominant diatom species in urban downstream samples; WT and F⁻ were correlated with reference and urban stream diatom composition. Our results demonstrate that the composition of diatom species was more sensitive to urbanization than the water physical and chemical parameters, and that diatom assemblage structure metrics more accurately assessed water quality. Some species, such as Amphora pediculus and Cocconeis placentula were among the dominant species in low nutrients stream sites; however, they were considered to be high nutrient indicators in some streams in USA. We suggest using caution in applying indicator indices based on species composition from other regions. It is necessary to build a complete set of diatom species data and their co-ordinate environment data for specific regions.     

Defining quantitative stream disturbance gradients and the additive role of habitat variation to explain macroinvertebrate taxa richness

Most studies dealing with the use of ecological indicators and other applied ecological research rely on some definition or concept of what constitutes least-, intermediate- and most-disturbed condition. Currently, most rigorous methodologies designed to define those conditions are suited to large spatial extents (nations, ecoregions) and many sites (hundreds to thousands). The objective of this study was to describe a methodology to quantitatively define a disturbance gradient for 40 sites in each of two small southeastern Brazil river basins. The assessment of anthropogenic disturbance experienced by each site was based solely on measurements strictly related to the intensity and extent of anthropogenic pressures. We calculated two indices: one concerned site-scale pressures and the other catchment-scale pressures. We combined those two indices into a single integrated disturbance index (IDI) because disturbances operating at both scales affect stream biota. The local- and catchment-scale disturbance indices were weakly correlated in the two basins (r = 0.21 and 0.35) and both significantly (p < 0.05) reduced site EPT (insect orders Ephemeroptera, Plecoptera, Trichoptera) richness. The IDI also performed well in explaining EPT richness in the basin that presented the stronger disturbance gradient (R² = 0.39, p < 0.001). Natural habitat variability was assessed as a second source of variation in EPT richness. Stream size and microhabitats were the key habitat characteristics not related to disturbances that enhanced the explanation of EPT richness over that attributed to the IDI. In both basins the IDI plus habitat metrics together explained around 50% of EPT richness variation. In the basin with the weaker disturbance gradient, natural habitat explained more variation in EPT richness than did the IDI, a result that has implications for biomonitoring studies. We conclude that quantitatively defined disturbance gradients offer a reliable and comprehensive characterization of anthropogenic pressure that integrates data from different spatial scales.     

Diatoms as indicators of isolated herbaceous wetland condition in Florida,USA

Benthic, epiphytic, and phytoplanktonic diatoms, as well as soil and water physical–chemical parameters, were sampled from 70 small (average 0.86 ha) isolated depressional herbaceous wetlands located along a gradient of human disturbance in peninsular Florida to (1) compare diatom assemblage structure between algal types; (2) develop biological indicators of wetland condition; (3) examine synecological relationships between diatom structure and environmental variables, with the ultimate goal of developing an index of biological integrity using a single assemblage. Collected diatom samples were enumerated to 250 valves and identified to species or subspecies. An assessment of wetland condition was made using a landscape-scale human disturbance score (Landscape Development Intensity index, LDI), calculated for each site using land use maps and GIS.Assemblages from both impaired and reference sites were compared using blocked multi-response permutation procedures, the percent similarity index, and visually examined using non-metric multidimensional scaling (NMDS). No ecologically significant compositional differences were found within sites. Mantel's test (Mantel's r = 0.29, p < 0.0001) and NMDS (stress: 14.52, variance: 78.5%) identified epiphytic diatoms as the most responsive to human disturbance. Strong significant correlations (|r_s| > 0.50, p < 0.05) were found between epiphytic NMDS site scores and soil pH, specific conductivity, water total phosphorous, and LDI, while soil pH, water color, soil TP, and turbidity were also significantly correlated (p < 0.05).Metrics to assess wetland condition were developed using epiphytic abundance data. Epiphytic taxa sensitive or tolerant to human landscape modification were identified using Indicator Species Analysis, and autecological indices relating diatom sensitivity to nutrients, pH, dissolved oxygen levels, saprobity, salinity, and trophic status were calculated. Fourteen final metrics were identified, scored on an ordinal scale, and combined into the Diatom Index of Wetland Condition (DIWC). The DIWC was highly correlated with the disturbance score (Spearman's r_s = −0.71, p < 0.0001), although the results need to be validated.     

Multivariate geostatistical methods for analysis of relationships between ecological indicators and environmental factors at multiple spatial scales

As all biodiversity-related variables, ecological indicators are influenced by environmental factors working at different spatial scales. However, assessing the relationship between environmental factors and ecological indicators is limited to a set of spatial scales determined a priori. This a priori assumption can hide important relationships, especially for ecological indicators with a complex spatial structure that can be driven, for example, by the influence of multiple pollutants with different dispersion ranges or by the influence of local and regional factors such as land-cover and climate. To relate ecological indicators and environmental factors without assuming a priori spatial scales of analysis, we used a Linear Model of Coregionalization. This method has been used in literature to analyze the joint distribution of biodiversity variables. Here we show that it can be used to gain insight into spatial patterns of relationships between ecological indicators and underlying environmental factors. We applied this method to a region of south-west Europe, relating data from land-cover, altitude and climate with an ecological indicator, the abundance of fruticose lichen species, known to be very sensitive to multiple environmental factors. Based on variogram analysis we identified distinct spatial scales of relationships between the ecological indicator and environmental factors. For each spatial scale we described relationships using Principal Component Analysis applied to the coregionalization matrices. This way we could assess how strong the relationship between each environmental factor and ecological indicator at each spatial scale was: at medium scales (c. 15 km) open spaces areas (a proxy for particle emissions) were more important; at larger scales (c. 45 km) open spaces, artificial areas (a proxy for gaseous pollutants) and also climate were preponderant. Thus, multivariate geostatistics provided a tool to improve knowledge on relationships between ecological indicators and environmental factors at multiple spatial scales without setting a priori spatial scales of analysis.     

Choice of biota in stream assessment and monitoring programs in tropical streams: A comparison of diatoms,macroinvertebrates and fish

The objective of this study was to compare the responses of diatoms, macroinvertebrates and fish to agriculture, urbanization and mining in the Manyame River Basin. Water quality sampling and benthic diatom, macroinvertebrate and fish community data were collected in April (end of the rain season) and September (dry season) 2013 at 44 sampling stations spread out across three land-use categories: commercial agricultural, communal agricultural and urban-mining areas. Commercial agricultural areas were relatively pristine as they were characterized by mature deciduous riparian forest strips which acted as riparian buffers thus protecting water resources from nonpoint source pollution. In communal agricultural areas a combination of poor agricultural practices (stream bank cultivation, overgrazing, soil erosions) and high human population densities had negative effects on water quality of streams draining these areas. Streams in urban-mining areas were highly stressed, being impacted primarily by physical habitat degradation and both point and nonpoint sources of pollution. A suite of environmental variables that varied with land-use pattern was assessed to find the combination of variables that best explained patterns of biota community composition. Community metrics i.e. the Trophic Diatom Index (TDI) based on diatoms, the South African Scoring system version 5 (SASS 5) based on macroinvertebrates and the Fish Assemblage Integrity Index (FAII) were used to determine the ecological status of study streams in relation to human-induced stressors. Data were also subjected to multivariate statistical techniques; canonical correspondence analysis (CCA), mantel test and cluster analysis to determine environmental gradients along which the diatom, macroinvertebrate and fish assemblages were distributed as well as to elucidate hypothesized differences in response to stressors among communities per land-use type. Using CCA, we assessed the individual importance of a set of environmental variables on each biotic community structure. ANOVA, showed a significant difference (p < 0.05) in physical and chemical variables among commercial agricultural, communal agricultural and urban-mining sampling stations with no significant differences (p > 0.05) between the 2 sampling periods. Based on CCAs carried out using individual variables, the strengths of relationships between diatoms and macroinvertebrates was generally high for nutrient levels, organic and metal pollution and other variables. However, fish assemblages showed a relatively low association with all water quality variables in the study; this might be explained by the high abundance of omnivores and air breathers which are able to tolerate a variety of environmental conditions. These patterns were also confirmed by the mantel test as well as the other CCAs carried out to investigate the simultaneous effects of environmental variables. These findings indicate that diatoms are more powerful indicators in accessing ecological stream/river quality and have potential for application in routine monitoring programs in tropical streams.     

Community congruence of plants,invertebrates and birds in natural and constructed shallow open-water wetlands: Do we need to monitor multiple assemblages?

Biomonitoring is a common means of evaluating wetlands. It is based on the premise that the community composition of one taxonomic group is indicative of overall biology and the underlying environmental conditions at a wetland. To be a good bioindicator, there must be adequate concordance between the indicator group and other biotic assemblages. Otherwise, multi-assemblage monitoring is necessary to glean a complete picture of wetland condition. In 32 sites ranging from reference wetlands to stormwater retention ponds, we evaluated concordance in community composition among the six most commonly monitored wetland assemblages: waterfowl, wetland dependent songbirds, aquatic macroinvertebrates, and plants in the wet meadow, emergent, and open-water vegetation zones. We also assessed agreement in environmental correlates among these six assemblages and investigated the impact of human disturbance on cross-assemblage concordance. We found that cross-assemblage concordance was positive (p < 0.03 in 14 of 15 pair-wise comparisons, p = 0.06 in 15th case), but relatively low (Mantel R values 0.11–0.37), suggesting that the assemblages are mediocre surrogates for one another. Yet, we found very strong agreement among environmental correlates of the six assemblages, especially along the first axis of assemblage-specific ordinations (mean Spearman rho = 0.923), indicating that despite low concordance, the six assemblages are likely responding to the same environmental gradients. Thus, while a single assemblage may not serve as a surrogate for the other assemblages, it should yield an adequate estimate of underlying environmental conditions and the degree of disturbance. Most important among the environmental correlates were sediment and water nutrient levels, shoreline slope, and the size of wet meadow and emergent vegetation zones. Perhaps most interestingly, we found that the strength of cross-assemblage concordance was greatest in reference wetlands and was lower (p ≤ 0.05) in constructed wetlands. This implies that cross-assemblage concordance present in undisturbed sites may not persist in disturbed wetlands where several of these cross-assemblage relationships deteriorate. Furthermore, a general change in cross-assemblage concordance may itself be indicative of human disturbance in wetlands.     

Reference diatom assemblage response to restoration of an acid mine drainage stream

A prevalent legacy of coal mining within Appalachia and elsewhere is acid mine drainage (AMD), which drastically alters both the chemical and biological components of the receiving waters. Hewett Fork is one such affected stream. Although AMD treatment has reduced acidity considerably downstream, the ability of this stream to sustain a biological community compared to those found in reference conditions remains unclear. To assess this, tiles colonized with diatom assemblages from a reference stream were transplanted into Hewett Fork in 5 locations along a 6.9 km stream length and sampled after one, three, and six weeks. Diatom assemblage structure metrics, including species evenness (J’), species richness (S), relative abundance of dominant taxon, and Shannon diversity (H′), as well as chlorophyll a concentrations, Bray–Curtis dissimilarities, and Acid Mine Drainage Diatom Index of Biotic Integrity (AMD-DIBI) scores were calculated for each site and sampling time. One-way ANOVAs of structural metrics showed significant differences (P ≤ 0.001) between the reference site and the 2nd and 3rd most upstream sites within the study reach for the duration of the study, with the exception of the relative abundance of dominant taxa at an intermediate site during the third week. Conversely, the most downstream Hewett Fork assemblage, located 11.6 km from the primary AMD input, did not differ significantly (P > 0.05) from that of the reference assemblage for any structural metrics after the initial sampling period. Throughout the study, only three sites obtained “good” AMD-DIBI narrative class:the reference site (weeks 1, 3, and 6), the most downstream site (weeks 1, 3, and 6; 11.6 km downstream of primary AMD input) and the uppermost site (weeks 1 and 6; 4.7 km downstream of primary AMD input). Results suggested that after an initial one-week acclimation period, assemblages at the uppermost and most downstream sites along the reach were relatively similar to those found in reference conditions, while sites within the middle region continued to show signs of impairment, although the factor(s) causing this impairment remain unknown. These findings suggest that while treatment has been effective on a site-specific basis, the expected linear-response to treatment was not achieved due to underlying factors that are inhibiting reference-like biological communities from reestablishing within the affected stream reach.     

Development of a diatom index of biotic integrity for acid mine drainage impacted streams

Acid mine drainage (AMD) resulting from extensive coal mining throughout Appalachia since the 1800s has caused a legacy of severe acid and dissolved metal loads to thousands of stream miles, which has critically impacted aquatic life and ecological attributes. Relationships of diatoms and macroinvertebrates with AMD have been established, but no index specifically designed to quantify AMD impacts using diatoms has been created, nor have the response of multiple organism groups been compared for their utility as indices assessing AMD severity.For the purpose of developing an effective assessment and management strategy for AMD impacted streams, this study created and tested a multi-metric AMD-diatom index of biotic integrity (AMD-DIBI) and compared its response to AMD severity with an already established multi-metric macroinvertebrate community index (ICI). In 2006, 41 sites in southeast Ohio were sampled that represented an AMD impact gradient and non-AMD impacted reference sites. Metrics comprising the AMD-DIBI were selected based on their responsiveness to AMD and nutrient impacts. In the following year, the AMD-DIBI and its metrics were tested on a validation dataset consisting of 18 sites in an AMD impacted watershed. Results indicated a significant correlation between AMD-DIBI and ICI scores, and both indices and all metrics were strongly correlated with water chemistry variables indicative of AMD pollution (P < 0.05). Stepwise multiple regression selected alkalinity and conductivity as most influential to AMD-DIBI (adjusted r² = 0.70) and ICI scores (adjusted r² = 0.66). Narrative classes (e.g., Poor, Fair, Good, and Excellent) defined by index scores provided effective classifications of AMD severity. When tested on the watershed scale, AMD-DIBI and its metrics very successfully quantified AMD gradients and coal mining impacts as indicated by canonical correspondence analysis. This newly developed AMD-DIBI will be very useful for assessing impairment, sensitivity, and recovery of diatom communities in streams damaged or threatened by coal mining activities. In addition, because the AMD-DIBI was very responsive to a gradient of AMD pollution, it could be used in future studies measuring the long-term status of streams and effectiveness of various remediation methods. This study highlights the responsive power of diatom-based metrics.     

Effects of urban and non-urban land cover on nitrogen and phosphorus runoff to Chesapeake Bay

The aim of this study was to determine the effects of catchment and riparian stream buffer-wide urban and non-urban land cover/land use (LC/LU) on total nitrogen (TN) and total phosphorus (TP) runoff to the Chesapeake Bay. The effects of the composition and configuration of LC/LU patches were explored in particular. A hybrid-statistical-process model, the SPAtially Referenced Regression On Watershed attributes (SPARROW), was calibrated with year 1997 watershed-wide, average annual TN and TP discharges to Chesapeake Bay. Two variables were predicted: (1) yield per unit watershed area and (2) mass delivered to the upper estuary. The 166,534 km² watershed was divided into 2339 catchments averaging 71 km². LC/LU was described using 16 classes applied to both the catchments and also to riparian stream buffers alone. Seven distinct landscape metrics were evaluated. In all, 167 (TN) and 168 (TP) LC/LU class metric combinations were tested in each model calibration run. Runs were made with LC/LU in six fixed riparian buffer widths (31, 62, 125, 250, 500, and 1000 meters (m)) and entire catchments. The significance of the non-point source type (land cover, manure and fertilizer application, and atmospheric deposition) and factors affecting land-to-water delivery (physiographic province and natural or artificial land surfaces) was assessed. The model with a 31 m riparian stream buffer width accounted for the highest variance of mean annual TN (r² = 0.9366) and TP (r² = 0.7503) yield (mass for a specified time normalized by drainage area). TN and TP loadings (mass for a specified time) entering the Chesapeake Bay were estimated to be 1.449 × 10⁸ and 5.367 × 10⁶ kg/yr, respectively. Five of the 167 TN and three of the 168 TP landscape metrics were shown to be significant (p-value ≤ 0.05) either for non-point sources or land-to-water delivery variables. This is the first demonstration of the significance of riparian LC/LU and landscape metrics on water quality simulation in a watershed as large as the Chesapeake Bay. Land cover metrics can therefore be expected to improve the precision of estimated TN and TP annual loadings to the Chesapeake Bay and may also suggest changes in land management that may be beneficial in control of nutrient runoff to the Chesapeake Bay and similar watersheds elsewhere.     

Multi-scale spatial and partitioning analyses of the reef-fish community composition of the Yucatan fringing reef system

Coral reef and reef-fish communities are complex systems that have spatial structures that are influenced by multiple factors and processes that interact at several scales. We analysed composition data of reef-fish communities in the Yucatan Peninsula fringing reef system using a hierarchical multi-scale survey design to elucidate the origins of the species distribution patterns through multivariate canonical ordination and partitioning analysis. Twelve reef sites were surveyed along the fringing reef system, nested in three sectors, each with distinctive human interests (tourism, fishing, protection). Line transects were placed on reef lagoon, reef front, reef slope and reef terrace habitat for a total of 480 transects. The communities were composed of many rare and some dominant species (e.g., Thalasoma bifasciatum, Acanthurus bahianus, A. coeruleus, Sparisoma aurofrenatum). Redundancy analyses revealed that environmental variables explained between 11% and 22% of the variation in fish community data according to the habitat analysed. Spatial variables, represented by the geographic coordinates of transects, explained a low percentage (2–5%) of variation in the fish community data. Among the environmental variables that significantly explained (P < 0.05) the patterns of variation on the fish community were depth, topographical complexity, and the percentage of the substratum covered by algae, calcareous floor, rock and rubble. Although the total explained variation and proportion of the variation explained by each set of variables were consistently low, most of the effects on the fish community composition data were significant. Explanations are proposed concerning those results. Evincing the distribution patterns on the biological systems of multi-species as well as the causes that shape them is key for conservation planning and essential in highly threatened regions, such as the Caribbean coast of the Yucatan Peninsula.     

Development,calibration, and validation of a littoral zone plant index of biotic integrity (PIBI) for lacustrine wetlands

We examine lacustrine wetland plant assemblages in the Central Corn Belt Plain portion of the Lake Michigan basin and developed a multimetric plant index of biotic integrity (PIBI). Our objectives were to determine the structural and functional attributes of littoral zone plant assemblages of least-impacted lacustrine wetlands, establish and test candidate metrics, statistically test and calibrate metrics, and finally validate a PIBI along a disturbance gradient. Of 35 candidate metrics, we chose 11 metrics that were grouped into four categories: species richness and composition, species tolerance, guild structure, and vegetation abundance. Based on Spearman correlations, we identified a suite of metrics, particularly those related to species richness and tolerance that had a strong response to human-induced habitat change. The overall PIBI correlated strongly with independent measures of habitat quality (p < 0.001) using a qualitative habitat index developed for lacustrine habitats. We validated the lacustrine PIBI by comparing index response to various landuse, landcover, and management types. Least impacted lakes and lakes classified as recreational or undergoing ecological restoration were not statistically separable and received the highest index scores, while the lowest scores were associated with industrial and residential land use. Least-impacted sites differ significantly (p < 0.001) from both industrial and residential lakes.     

Validating Landsat-based landscape metrics with fine-grained land cover data

ContextModerate-grained data may not always represent landscape structure in adequate detail which could cause misleading results. Certain metrics have been shown to be predictable with changes in scale; however, no studies have verified such predictions using independent fine-grained data.ObjectivesOur objective was to use independently derived land cover datasets to assess relationships between metrics based on fine- and moderate-grained data for a range of analysis extents. We focus on metrics that previous literature has shown to have predictable relationships across scales.MethodsThe study area was located in eastern Connecticut. We compared a 1 m land cover dataset to a 30 m resampled dataset, derived from the 1 m data, as well as two Landsat-based datasets. We examined 11 metrics which included cover areas and patch metrics. Metrics were analyzed using analysis extents ranging from 100 to 1400 m in radius.ResultsThe resampled data had very strong linear relationships to the 1 m data, from which it was derived, for all metrics regardless of the analysis extent size. Landsat-based data had strong correlations for most cover area metrics but had little or no correlation for patch metrics. Increasing analysis areas improved correlations.ConclusionsRelationships between coarse- and fine-grained data tend to be much weaker when comparing independent land cover datasets. Thus, trends across scales that are found by resampling land cover are likely to be unsuitable for predicting the effects of finer-scale elements in the landscape. Nevertheless, coarser data shows promise in predicting fine-grained for cover area metrics provided the analysis area used is sufficiently large.     

A review of indicators of estuarine tidal wetland condition

This review critically evaluates indicators of tidal wetland condition based on 36 indicator development studies and indicators developed as part of U.S. state tidal wetland monitoring programs. Individual metrics were evaluated based on relative scores on two sets of evaluation factors. A rigor score evaluated metric development based on conceptual relevance, indicator development method, degree of independent validation, and temporal and spatial extent tested. An applicability score evaluated metrics based on cost of data collection, probable spatial extent of applicability, technical complexity, and indicator responsiveness. The majority of indicators could be classified as biotic condition indicators (81%), with vegetation (37%) and macroinvertebrate (28%) metrics composing the largest proportion. Most metrics provided a conceptual model or scientific justification (97%), were developed by correlation to environmental gradients (46%), were tested over multiple seasons or years (49%) and at multiple sites (88%). Few were independently validated (18%). Average rigor score was 10 (on a scale of 0–25) and ranged between 1 and 21. Highest rigor scores were for trematode community metrics (community similarity index, species richness) and metrics of grass shrimp (Palaemonetes pugio) individuals (gene expression, relative fecundity, embryo hatching success, larval survival). Most metrics had a high cost of data collection (63%), required field and laboratory processing (84%), would be applicable across the U.S. (72%), and were responsive to the variable of interest (44%). Mean applicability score was 4.9 (range: 2–8). Highest scores were found for metrics that only required field collection of data using simple or no instrumentation. Lowest scoring metrics required expensive equipment, specialized taxonomic knowledge, complex laboratory analysis, and/or culturing of organisms. Scores for individual metrics were grouped by indicator, then averaged and rescaled between 0 and 100 to provide a composite evaluation of the indicator they measured. Among major indicator types, biotic indicators had the highest rigor scores (mean = 44, range 20–79), followed by indicators of chemical/physical characteristics (mean = 36, range 16–56), landscape condition (mean = 31, range 24–37), and hydrology/geomorphology indicators (mean = 21, range 4–52). In contrast, biotic indicators scored lowest for applicability (mean = 58, range 25–100) and indicators of landscape condition scored highest. The results of this review suggest that the development and selection of tidal wetland indicators could be vastly improved by employing a standardized development methodology that provides uniform information about each indicator. In addition, tidal wetland indicator research should focus on the development of indicators of ecological processes and disturbance regimes.     

The reliability of a composite biodiversity indicator in predicting bird species richness at different spatial scales

Several biodiversity features can be linked to landscape heterogeneity, that, in turn, can be informative for management and conservation purposes. Usually, the more the landscape is complex the more the biodiversity increases. Biodiversity indicators can be a useful tool to assess biodiversity status, in function of landscape heterogeneity. In this study, we developed a biodiversity indicator, based on Shannon diversity index and built from distribution maps of protected species. With such an approach, we seek to evaluate the feasibility of using a combination of target species as a surrogate for assessing the status of the whole bird community. Our approach was spread over multiple spatial scales, to determine which was the most informative. We selected four species protected by European regulation and generated a presence-absence map from species distribution modelling. We, therefore, used the FRAGSTATS biodiversity metric to calculate Shannon index for the overlapped presence-absence maps, at two spatial scales (500 m and 1000 m). Then, the relationships with the whole community was assessed through generalised least square models, at the spatial scale of 4 ha, 9 ha and 25 ha. Results showed that the higher rate of variability of community was explained by the biodiversity indicator at 1000 m scale. Indeed, the more informative spatial scale for the whole bird community was 9 ha. In addition, a pattern emerged about the relationships between biodiversity indicator and community richness, that is worth of further research. Our study demonstrates that the usefulness of surrogate species for biodiversity and community assessment can become clear only at a certain spatial scales. Indeed, they can be highly predictive of the whole community, and highly informative for conservation planning. Moreover, their use can optimize biodiversity monitoring and conservation, focusing on a small number of noteworthy species.     

Growth responses of Zostera capricorni to estuarine sediment conditions

This study comprised (1) a field survey of intertidal seagrass (Zostera capricorni) biomass, cover and photosynthetic potential and sediment characteristics at a range of contrasting sites in three New Zealand harbours, and (2) a microcosm experiment comparing plant responses to sediments from extant versus historical seagrass sites. The field survey showed that the sediment physico-chemical characteristics were generally consistent with the limited previous reports for Zostera environments, although the total P concentration range was higher (0.08–0.72 mg P g⁻¹). Overall, 52% of variation in seagrass cover was explained by sediment water content (R = 0.54) and organic content (R = −0.56). Twenty-two percent of variation in seagrass biomass was explained by sediment total P and redox potential (both R = −0.35). Intra-harbour seagrass–sediment relationships were more significant (explaining up to 82% of plant variation) but harbour-specific. In the microcosm experiment, threefold higher Z. capricorni biomass was maintained on extant than historical sediments but not conclusively linked to measure sediment characteristics. Overall, the results of this study demonstrate that significant relations can exist between estuarine sediment conditions and Z. capricorni growth responses, and suggest that detrimental change in sediment conditions may be a contributing factor in seagrass decline.