The influence of substrates and physicochemical factors on the composition of diatom assemblages in karst springs and their applicability in water-quality assessment

The methodological standards of monitoring programs recommend collecting material from stones, but this community is often absent in slow-flowing waters, whereas a great number of localities do offer bryophytic material for sampling. In order to investigate whether the kind of substrate (epilithic or bryophytic) influences the diatom-based biomonitoring results the structure of diatom assemblages from 13 anthropogenically altered springs were investigated. To examine differences between diatom assemblages, the species richness of diatom assemblages, H′ (Shannon-Wiener diversity index), Trophic Diatom Index (TDI) and Specific Pollution sensitivity Index (SPI) were compared. The TDI and SPI differed significantly (P = 0.002 and P = 0.004) between springs located in villages and springs away from villages, and did not differ significantly between epilithic and bryophytic samples (P > 0.05). No significant differences in species richness or the diversity index (H′) (P > 0.05) were found between the epilithic and bryophytic diatom assemblages. This study suggests that both the kinds of substrates can be successfully used for diatom-based water-quality biomonitoring in mesotrophic and eutrophic environments.


Introduction
Diatoms are one of the most commonly used organisms for assessing natural and human-related features of aquatic environments. Their value as bioindicators was appreciated long ago, and diatom analysis is still done for reconstruction of the paleoenvironment and for water quality/environmental health assessment. Their usefulness is due to the quick and precise response of the diatom assemblage structure to changing conditions (Stoermer & Smol, 1999). To standardize the methodology of using diatoms for water-quality assessment and to minimize substraterelated variability of diatom assemblage structure, it has been recommended to study communities developing on a single type of substrate; for example, in analyses of running waters the epilithon is recommended (e.g. Round, 1991Round, , 1993Kelly et al., 1998;Townsend & Gell, 2005).
Once the general regularities of how specific diatom assemblages form in response to pH, salinity, trophic state and several other physicochemical parameters were recognized, a number of detailed studies of habitat-diatom assemblage relationships were undertaken. The relationships between benthic communities and substrate type have been the subject of several investigations and extensive discussion (Douglas, 1958;Stevenson & Hashim, 1989;Jüttner et al., 1996;Rothfritz et al., 1997;Sabater et al., 1998;Potapova & Charles, 2005;Townsend & Gell, 2005;Cantonati & Spitale, 2009). As the relationships between substrate and the organisms inhabiting it are much more complicated in benthic and periphytic communities (Cox, 1990;Burkholder, 1996;Townsend & Gell, 2005) than in plankton, the effects of different factors on assemblage structure in running waters are much less well understood. The influence of biotic and abiotic substrates on diatom assemblages has been studied in rivers on broad geographical scales (e.g. Jüttner et al., 1996;Rothfritz et al., 1997;Potapova & Charles, 2005;Townsend & Gell, 2005) but has received much less attention in springs (Cantonati, 2001;Cantonati & Spitale, 2009). Springs provide excellent opportunities for such studies because the fluctuation of water properties is minimal (Sabater & Roca, 1990;Cantonati et al., 2006;Dumnicka et al., 2007;Cantonati & Spitale, 2009) and the input of allochthonous organisms from other environments is limited. The available data indicate that the geomorphological features of an area and the water characteristics are major factors governing the diatom distribution in springs (Werum, 2001;Werum & Lange-Bertalot, 2004;Cantonati et al., 2006).
Most of the data on diatoms in springs come from mountain areas in natural or near-natural conditions (Sabater & Roca, 1990, 1992Cantonati, 1998;Battegazzore et al., 2004;Werum & Lange-Bertalot, 2004), where the epilithon generally occurs. In lowland areas, the role of epilithic communities is smaller, because slow flow allows sediment to accumulate and become the predominant substrate in springs, together with submerged vascular and nonvascular macrophytes. Macrophytes such as bryophytes occur in most springs and are one of the most common habitats for a very diverse microflora (Cantonati & Spitale, 2009;Hájkova et al., 2011). A comparison of the structure of diatom assemblages in bryophytic and epilithic communities from non-mountain, mesotrophic springs can provide applicative data to use for water-quality biomonitoring; epilithic assemblages are recommended for waterquality assessment, but such assemblages are virtually absent in several aquatic lowland systems (e.g. Potapova & Charles, 2005). In order to investigate whether the kind of substrate (epilithic or bryophytic) commonly used for assessment of water quality diatom influences, e.g. Trophic Diatom Index (TDI; Kelly & Whitton, 1995) and Specific Pollution sensitivity Index (SPI; CEMAGREF, 1982), the diatom assemblages from 13 anthropogenically altered springs were investigated.

Site description
The springs in question are situated in the Wy_ zyna Krakowsko-Częstochowska Upland in southern Poland (Table 1). This upland, covering 2,650 km 2 , extends from the Carpathian foothills in the vicinity of Kraków northward to the town of Częstochowa. It is ca. 80-km long and ca. 20-km wide, and rises to 350 m a.s.l. on average. Geologically, it is karst formed in Upper Jurassic calcareous rock up to 200-m thick and further altered since the Cretaceous (Ró _ zkowski, 1996), partly covered by Quaternary loess in the southern part and by postglacial sands in the northern part (Dynowska, 1983). Most of the spring water belongs to the bi-ionic HCO 3 -Ca 2? category. Local hydrogeology is the main natural factor determining the spatial diversity of spring water chemistry (Siwek, 2004). The springs from the southern part of the upland are supplied from the Jurassic aquifer and are typically small, with discharge between 0.4 and 15 l s -1 (Dumnicka et al., 2007). Rapid circulation in a shallow system of large karstic channels limits bedrock weathering there. The northern part of the study area belongs to the Cretaceous aquifer, and the springs supplied by it have greater discharge. The range of spring water chemistry is much broader there, due to diverse geology, several local water tables, and higher mineral dissolution of the bedrock (Siwek, 2004). The springs of the Krakowsko-Częstochowska Upland are very vulnerable to pollution because surface water flows rapidly into groundwater systems through natural conduits and is under significant anthropogenic pressure. Most of the springs in this study are located in highly populated areas. The selected springs differ in discharge, type of surroundings, degree of artificial modification and ionic composition (Tables 1, 2).

Methods
Epilithic, bryophytic samples and water for physical and chemical analyses were collected in September-November 2004. Composite samples of at least four subsamples from submerged stones were collected by brushing the upper surface of submerged stones and squeezing from submerged bryophytes. Although previous studies did not show essential differences in the structure of epilithic diatom associations on stones of differing size in slow-flowing waters (Wojtal & Sobczyk, 2006), we collected the epilithon samples from stones of similar size. Some of the environmental variables were measured at the time of diatom sampling. Water temperature was taken and specific conductivity and pH were measured using an Elmetron CC-102 conductivity metre and a CC-104 pH metre. Other chemistry data were obtained by ion chromatography (Dionex 100), atomic absorption spectrophotometry and flame AAS in the Institute of Botany laboratory in Szarów, Poland.
To remove carbonates, 10% HCl was added to the samples in 250-ml glass beakers and boiled for 15 min after 24 h. Then the material was washed at least five times with distilled water, with 24 h settling time between each decanting. Then the samples were boiled in 30% H 2 O 2 with small amounts of KClO 3 , or using concentrated ([30%) acids (H 2 SO 4 , HNO 3 ) to remove any leftover organic matter not digested by peroxide. After washing several times with distilled water the material was air dried on coverglasses and mounted in Naphrax Ò resin. The diatoms were observed with a Nikon Eclipse 600 light microscope with Differential Interference Contrast.
The relative abundance of particular taxa and the species richness of the assemblages were estimated on the basis of at least 300 diatom valves per sample. To examine differences between diatom assemblages, we compared the species richness of epilithic and bryophytic assemblages, and the following indices for assessing water quality were computed with Omnidia 4.1 software (Lecointe et al., 1993): H 0 (Shannon-Wiener diversity index), TDI and SPI.
Differences in Shannon-Wiener diversity index, TDI and SPI between bryophytic samples and epilithic     samples and between springs under anthropogenic pressure (in or near villages) and springs away from villages were analyzed using ANOVA with three factors: substrate, spring and location type and springs nested in location type. Taxa accounting for [1% of the diatom assemblages at more than three of the studied springs, after logarithmic transformation (natural logarithm) were included in the ordination analyses. To explain variability in the species composition the data were analyzed by principal correspondence analysis (PCA) using correlation matrix. To explain species-environment relationship the species data were analyzed by redundancy analysis (RDA) (CANOCO v4.5 software; ter Braak & Š milauer, 2002;Lepš & Šmilauer, 2003). The analyses were performed using the whole data set (bryophytic and epilithic diatom assemblages together) and then separately for bryophytic diatom assemblages only and epilithic diatom assemblages only.

Results
The pH range was \1 U (pH 7.0-7.7) for the whole data set. All analyzed localities showed comparable calcium concentration, alkalinity and temperature (Table 2). There was a relatively high variation of specific conductivity in the springs, ranging from 210 (Harcerza spring) to 727 lS cm -1 (spring in Węgrzynów), and considerable variation of nitrate concentration from 10.7 (Ruskie spring) to 38.2 mg l -1 (spring in Łany Wielkie). The diatom assemblages were dominated by neutrophilous and alkaliphilous, mesotraphentic and eutraphentic, aand b-mesosaprobous diatoms (Table 3).
We identified a total of 136 diatom taxa, the vast majority of which were typically epilithic or typically epiphytic forms. Only 17 were common (i.e. present at [10% relative abundance in at least one sample); of those, 10 reached[20% relative abundance in at least one sample. Average species richness in epilithic samples was 18 diatom taxa, ranging from 5 to 34; the corresponding figures for bryophytic samples are 18.5 and 12-29. Average species richness differed between southern (away from villages; less eutrophicated) and northern (in or near villages; more eutrophicated) waters: 24 and 12, respectively, in the epilithon, and 19 and 18 in bryophytic samples. The Shannon-Wiener diversity index (H 0 ) did not differ significantly Table 2 continued Locality no. .0% when all taxa having over 10% relative abundance are included. Twenty taxa accounting for[1% of the diatom assemblages at more than three of the springs were found in the epilithon, and 13 such species in the epibryon; those species were included in the statistical analyses. The species-locality relationship analysis (PCA) of all the samples revealed that the first two principal components accounted for 18.9% (axis 1) and 16.1% (axis 2) of the total variation of species composition (Fig. 1). The first axis separates samples collected from springs situated in or near villages (mainly areas covered by postglacial sands) from those situated in localities less affected by human impacts (mainly loess-covered areas) (Table 1). Axis 2 separates samples collected from epilithic and bryophytic substrates. The main positive contributions to the results were from Caloneis fontinalis, Meridion circulare, Planothidium lanceolatum, Gomphonema cf. pumilum, Diatoma mesodon and Diploneis krammeri Lange-Bertalot & E. Reichardt, which were more abundant in bryophytic samples, and Denticula tenuis, Eolimna minima (Grunow) Lange-Bertalot, Staurosirella pinnata, Staurosira venter, Planothidium reichardtii, Achnanthidium    pyrenaicum and Gomphonema micropus, which were more numerous in epilithic samples. RDA based on combined data from epilithic and bryophytic assemblages showed that the first two axes accounted for 24.9% (axis 1) and 13.5% (axis 2) of the total variability in the species composition (Fig. 2). More than half of the total variability of the relationship between species composition and environmental variables was explained by the first two axes: 35.0% by axis 1 and 19.9% by axis 2, but these axes were not significant (P [ 0.05).
In an additional PCA restricted to the bryophytic diatom assemblages the first two principal components accounted for 31.8% (axis 1) and 21.9% (axis 2) of the total variation of species composition. RDA for these assemblages showed that the first two axes explain 29.7% (axis 1) and 9.4% (axis 2) of variance in species composition and 52.5% (axis 1) and 17.6% (axis 2) the relationships between species composition and environmental variables. The first axix (RDA) was significant (P = 0.008) and correlated positively with dissolved oxygen concentration and negatively with specific conductivity and nitrate concentration. The first axis again separated the samples collected from springs with higher dissolved oxygen concentration and lower specific conductivity (Harcerza, Siamoszyce and Ruskie springs), with predominance of Gomphonema cf. pumilum, Caloneis fontinalis, Cocconeis placentula and Diatoma mesodon, which were replaced by Planothidium lanceolatum and Eolimna minima in springs with higher nutrient concentration and higher specific conductivity (in Wegrzynów and Łany Wielkie springs). Higher relative abundance of Staurosirella pinnata was related to springs with sandy bottom and high discharge (e.g. Imbramowice spring); the occurrence of Achnanthidium minutissimum, Meridion circulare and Amphora pediculus was typical of small springs with gravel and silt bottom. PCA analysis of epilithic diatom assemblages only revealed that the first two principal components accounted for 30.2% (axis 1) and 15.4% (axis 2) of the total variation of species composition. RDA for these assemblages showed that the first axes explain 23.3% (axis 1) and 11.4% (axis 2) of variation in a structure of diatom assemblages and 37.5% (axis 1) and 18.3% (axis 2) in the relationship between species composition and environmental factors, but both axes were not significant (P [ 0.05). The first axis separated samples collected in springs with elevated specific conductivity, high nitrates, sulphates, calcium and orthophosphates, from those less eutrophicated, with higher dissolved oxygen concentration. The largest contribution to separation of the more eutrophicated springs was from Eolimna minima, Planothidium lanceolatum, Staurosirella pinnata and Staurosira venter. The rest of the springs with higher dissolved oxygen and lower specific conductivity were mainly colonized by Achnanthidium minutissimum, Amphora pediculus, Planothidium reichardtii, Cocconeis placentula, Gomphonema cf. pumilum, Encyonema minutum (Hilse) D.G. Mann, Nitzschia fonticola Grunow, Denticula tenuis and Achnanthidium pyrenaicum.

Discussion
The minimal fluctuation of chemical parameters in springs (Sabater & Roca, 1990;Cantonati et al., 2006;Dumnicka et al., 2007), as compared to other freshwater systems and especially lakes, is related to the continuous renewal of the aquatic environment in springs. In karst areas, where pollutants are easily transported into the aquifer, spring water chemistry reflects changes in the whole area of the aquifer feeding the springs (Siwek, 2004). The composition of natural spring water is related to natural factors such as the weathering resistance of carbonate rocks, reflected mainly in total mineral content. High nitrate concentrations have been observed in several springs of this upland for decades, and are a permanent feature (Baścik et al., 2001;Siwek, 2004;Wojtal et al., 2009). The increased specific conductivity and Na ? , Cl -, K ? and SO 4 2concentrations are anthropogenic in nature (Ró _ zkowski, 1996;Baścik et al., 2001;Siwek, 2004), attributable to animal husbandry, agriculture and the absence of sewage treatment systems in villages (Galas, 2005). The differences in water quality (Tables 2, 3) were clearly reflected in the composition of epilithic and bryophytic diatom assemblages. The high bacterial contamination in springs in villages (Table 2) indicates pollution from wastewater.
Stones are one of the most common substrates in fastflowing waters, but are much rarer in lowland areas. Boulders, cobbles and gravel were virtually absent in 9 of the 36 springs checked preliminary to the study, and at 3 more localities there were no epilithic diatom assemblages present, presumably as a consequence of insufficient light availability. Those springs were excluded from further study. In most springs, however, there were submerged bryophytes. Despite negative interaction-physical and sometimes chemicalbetween epiphytic algae and host macrophytes, epiphytic diatoms do secure a location for growth, improved access to light and shelter from scouring (Burkholder 1996). According to Burkholder (1996), loosely attached epiphytes generally do not exhibit substrate preferences. Moreover, the preferences of epiphytic algae for specific host plants decrease with increasing nutrient concentration (Eminson & Moss, 1980;Hall & Smol, 1999). The clearest relationships between hosts and diatom assemblages are known from oligotrophic waters, where nutrient limitation is one of the most important factors governing diatom distributions. In these conditions, macrophytes are an additional nutrient source (Wetzel, 1983) or even the main one (Eminson & Moss, 1980). Unfortunately, the data on the diatom floras of spring bryophytes are sparse (Cantonati & Spitale, 2009). Most epilithic and bryophytic diatoms share similar characteristics enabling them to adapt to these habitats, including small size and adaptations enabling adhesion. Some species seemed to prefer certain substrates (e.g. Caloneis fontinalis in submerged bryophytes, Achnanthidium pyrenaicum on stones) but the general response of diatom assemblages as a whole to water chemistry was very similar. The results are in agreement with those from other studies of running waters (Jüttner et al., 1996;Rothfritz et al., 1997;Rott et al., 1998;Potapova & Charles, 2005) and springs (Cantonati & Spitale, 2009).
The relatively low species richness in the studied springs corresponds well with results from springs in the Pyrenees (Sabater & Roca, 1990, 1992 and Alps (Cantonati, 1998;Cantonati & Spitale, 2009), but could also be related to anthropogenic alteration of spring water quality in that upland, where some factor or set of factors limits the growth of the most sensitive taxa. In some springs, the epilithic diatom assemblages were dominated by the neutrophilous Achnanthidium minutissimum, whereas bryophytic samples (Fontinalis antipyretica L.) were dominated by the alkaliphilous Planothidium lanceolatum. The differences might be related to host (bryophytes) photosynthetic processes, as host biological processes drive a rapid increase of pH (Wetzel, 1983) and precipitation of calcium carbonate (Burkholder, 1996). In addition, Planothidium lanceolatum prefers substrates that provide protection from abrasion (Biggs, 1996), and such conditions can be expected in submerged bryophytes (Cantonati & Spitale, 2009). Note that although the observed assemblages consisted of common diatoms having up to 96% relative abundance, of which several probably are habitat generalists; they represented only one-third of the total number of taxa identified in diatom counts. Those findings are in agreement with Townsend & Gell (2005) who observed the occurrence of common taxa in large numbers on all substrates, whereas taxa unique to a single substrate attained 0.1-2% relative abundance. In our study, the rest of diatoms (85 taxa) were rare (\1%), and their presence varied greatly among springs. Sometimes the more common diatoms also reflected the individual character of the assemblages in certain springs. For example, the spring in Węgrzynów was the only one where Fragilariforma nitzschioides was observed, reaching over 60% relative abundance. Unlike in springs in the Alps (Cantonati & Spitale, 2009), we did not find significant differences in species richness or the diversity index (H 0 ) between the epilithic and the bryophytic diatom assemblages.
The general variability of the aquatic environment and especially of total ionic strength is known to play a very important part in determining the distribution of diatoms (Sabater & Roca, 1990;Soininen, 2007;Cantonati & Spitale, 2009;Angeli et al., 2010). Potapova & Charles (2005) demonstrated that the effect of ionic strength is of much greater significance than the type of substrate colonized. According to Hall & Smol (1999), the diatom assemblages on various substrates become uniform in meso-and eutrophic waters, where nutrient limitation is weak. Our results support this. Changes in diatom assemblages were associated mostly with nutrient concentration, which was higher in springs in or near villages, as reflected in the TDI and SPI indices. Specific conductivity, dissolved oxygen and correlated factors had a strong effect on the structure of the diatom assemblages studied.
The methodological standards of monitoring programs recommend collecting material from stones (Round, 1991(Round, , 1993Kelly et al., 1998), but this community is often absent in slow-flowing waters, and that limits the method to localities from which such material can be sampled. A great number of localities do offer bryophytic material for sampling. In our study, the data from both the epilithic and the bryophytic samples yielded diatom indices (TP, SPI) that reflected differences in water quality between springs, demonstrating that the two substrate types can be sampled for water quality assessment in slowflowing mesotrophic and eutrophic waters.