Next Article in Journal
The Role of Debris Cover in Catchment Runoff: A Case Study of the Hailuogou Catchment, South-Eastern Tibetan Plateau
Next Article in Special Issue
Conserving Rivers and Their Biodiversity in Tanzania
Previous Article in Journal
Reducing Hydro-Meteorological Risk by Nature-Based Solutions: What Do We Know about People’s Perceptions?
Previous Article in Special Issue
Ecological and Conservation Value of Small Standing-Water Ecosystems: A Systematic Review of Current Knowledge and Future Challenges
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Water
Volume 11
Issue 12
10.3390/w11122602
water-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diatom Biodiversity in Karst Springs of Mediterranean Geographic Areas with Contrasting Characteristics: Islands vs Mainland

by
Giuseppina G. Lai
1,*,
Sara Burato
2,
Bachisio M. Padedda
1,*,
Raffaella Zorza
2,
Elisabetta Pizzul
3,
Cristina Delgado
4,
Antonella Lugliè
1 and
Marco Cantonati
5
1
Department of Architecture, Design and Urban Planning, University of Sassari, via Piandanna 4, 07100 Sassari, Italy
2
ARPA FVG (Regional Agency for Environmental Protection Friuli Venezia-Giulia), Via Cairoli 14, 33057 Palmanova, Italy
3
Department of Life Sciences, University of Trieste, via Giorgieri 10, 34127 Trieste, Italy
4
Department of Ecology and Animal Biology, University of Vigo, Campus As Lagoas Marcosende, 36330 Vigo, Spain
5
MUSE–Museo delle Scienze, Limnology & Phycology Section, Corso del Lavoro e della Scienza 3, 38122 Trento, Italy
*
Authors to whom correspondence should be addressed.
Water 2019, 11(12), 2602; https://doi.org/10.3390/w11122602
Submission received: 30 October 2019 / Revised: 28 November 2019 / Accepted: 5 December 2019 / Published: 10 December 2019
(This article belongs to the Special Issue Multiplicity, Characteristics, Main Impacts, and Stewardship of Natural and Artificial Freshwater Environments: Consequences for Biodiversity Conservation)
Browse Figures
Review Reports Versions Notes

Abstract

:
Karst ecosystems are considered as priority environments for the protection of biodiversity on a global scale. This study provides a first comparative analysis of epilithic diatom flora from karst springs in two Mediterranean geographic areas (Spain and Italy) with contrasting characteristics (islands vs mainland). We investigated twenty-three springs with different anthropogenic impact levels once in the winter season between 2007 and 2017 (N = 23). A total of 176 diatom taxa (56 genera) were found of which 101 (44 genera) were observed in single sites. A general good biotic integrity was revealed by structural indices (species richness, diversity and evenness). However, crenophilous species were generally present and abundant in less impacted springs. Comparing islands and mainland, significant differences were found in species composition and diversity (H’) based on multivariate analyses (global R = 0.610; p = 0.001) and t-test (t = 2.304; p = 0.031). Discharge and Cl− were the most significant variables in determining diatom assemblages. Our results confirm the role of springs as multiple ecotones and refuges for rare species and suggest that the geographic insularity may be an important factor in maintaining diatom biodiversity.

1. Introduction

Karst springs are freshwater environments of great ecological value and strategic water resources in the Mediterranean area since they provide valuable ecosystem services, among which habitat for high biodiversity and drinking water supply [1]. In this geographic area, due to the strong climatic seasonality and lower water availability during dry periods (summer–autumn), several large springs have been captured since ancient times and they are almost the exclusive drinking water source for many urban centres [2,3].
Karst ecosystems have been indicated as priority environments for the protection of biodiversity on a global scale [4]. Unlike other freshwater ecosystems, several springs are still relatively natural environments [5] and can provide specific habitats for crenobiontic and crenophilous species [6,7]. As multiple ecotones, they include mosaic-like structures or patches of different aquatic, semi-aquatic and semi-terrestrial microhabitats (e.g., mosses and debris layers) which can host a large number of species and serve as refuges for rare and threatened organisms and least-impaired habitat relicts (LIHRe) [7,8].
Karst springs are intrinsically sensitive to variations in the hydrological cycle, water abstraction and pollution [9,10], and are more exposed to these impacts in the Mediterranean area because of the variable rainfall regime, climate change and increasing water use e.g., [11,12]. Hydromorphological alterations in these ecosystems can make some species vulnerable, especially those with limited dispersal, and cause an impoverishment of the biota, in particular, aquatic organisms. In fact, the quality and quantity of water and substrates are the main factors influencing the development of aquatic plant communities [13]. Further, habitat alteration is one of the major drivers of species loss [7]. Biodiversity, especially of microorganisms, play an important role in abiotic processes and maintenance of karst systems [14]. Recent studies highlighted the need for integrated multidisciplinary approaches, also based on biology and ecology, for sustainable use and proper management of karst springs [1,15].
Diatoms are an abundant group of microalgae in springs [16]. As primary producers they sustain the food webs in these ecosystems, significantly contributing to their functioning [17,18]. They are also useful indicators of important environmental features of springs, such as water quality (chemical composition and trophic conditions), water flow permanence, hydraulic regime, discharge, shading and lithology of the aquifers e.g., [19,20].
In the Mediterranean area, studies on diatoms from karst springs were carried out mainly in mountain regions in Slovenia [21], Croatia [22], Bosnia and Herzegovina [23], and Italy [18,24,25]. Only a few studies focused on low-altitude springs on islands: Majorca, the largest of the Balearic Islands (Spain) [26] and Sardinia (Italy), the second-largest island in the Mediterranean basin [27,28,29,30,31,32]. In general, these studies pointed out high levels of biodiversity and the presence of species new to science (e.g., Navicula veronesis Lange-Bertalot & Cantonati and Sellaphora gologonica G.G. Lai, Ector & C.E. Wetzel). They also provided information on the ecology of species and their relationships with environmental variables, considering both natural and anthropogenic disturbances, and suggested good prospects in the use of diatoms for the definition of ecological integrity and vulnerability aspects of karst springs. However, these studies were limited and focused on a single spring or groups of springs in specific areas, and comparisons on a wider geographic scale are not available in the literature.
This study compares the epilithic diatom flora from karst springs of four regions in two Mediterranean geographic areas (Spain and Italy). We hypothesized significant differences among diatom assemblages due to contrasting characteristics (islands vs mainland) and local factors. The main objectives were: (1) to examine the species composition and structure of diatom assemblages; (2) to evaluate similarities/dissimilarities among geographic areas; (3) to explore relationships between diatom taxa and environmental variables.

2. Materials and Methods

2.1. Study Sites

For this study, 23 springs were selected from four datasets built on a regional scale: Majorca Island (Spain), and Sardinia, Veneto and Friuli (Italy) (Figure 1, Table 1).
The three springs studied in Majorca Island are located in the Tramuntana mountain range, in the northern part of the island, and cover an altitudinal range of 56–471 m a.s.l. The Tramuntana Mountains extend parallel to the north-western coast and present numerous springs with different sizes and discharge [33]. The springs studied in this area are perennial and rheocrenic systems, and some of them are captured for drinking purposes.
In Sardinia, the eight springs studied are located in different areas of the island and distributed over an elevational range of 55–780 m a.s.l. They are mostly perennial and rheocrenic systems. Only two springs, SG and VA are classified as limno-rheocrenic and seasonal, respectively. The springs SG and DN are “Natural Monuments” according to Regional Law 31/1998 that defines the areas of noteworthy naturalistic value for conservation and protection in Sardinia. Five springs are located in Sites of Community Importance (SCI) for the Mediterranean biogeographic region: UL and FO (ITB 021107 “Monte Albo”) SG (ITB022212 “Supramonte di Oliena, Orgosolo e Urzulei–Su Sercone”), VA (ITB021156 “Monte Gonare”) and DN (ITB041111 “Monte Linas–Marganai”). Most of these springs are captured for drinking purposes.
In the Veneto Region, almost all springs (five out of six) are located in the south-eastern Alpine foothills, in a geographic area called Lessinia. One spring (PR) is located on the Monte Baldo mountain range (south-eastern Alps). They cover an elevation range of 183–675 m a.s.l. These springs are perennial systems and are virtually all captured (tapped) to some extent. Three springs (BR, BA, and TO) are classified as rheo-limnocrenic, and three springs (FC, LF and PR) as limno-rheocrenic systems. This is typically due to artificial (concrete) water collection basins that occupy part of the spring-head.
Finally, in Friuli, the six springs studied cover an altitudinal range of 454–1249 m a.s.l and are perennial and typical rheocrenic systems. The springs BM and GD are included in the list of geosites of regional and national interest, respectively [34]. The site GD is classified as “Alpine Cave” (habitat code SC1) [35]. The springs RR is captured for drinking purposes.

2.2. Sampling

Sampling was carried out once in the winter season between 2007 and 2017 (N = 23). The geographic position and altitude of the springs were recorded by a GPS. Water temperature, pH and conductivity were measured in situ using digital multi-parametric probes [24,26,32].
Anthropogenic disturbances were evaluated in an ordinal scale with disturbance scores: 1 = low; 2 = medium; 3 = high [24]. Shading conditions and current velocity were assessed using a five-score scale [36]. In most springs, discharge was measured by the volumetric method using a graduated bucket and a chronometer. Repeated measures were made close to the diatom sampling points. Only in some large springs (TR, BM and AL), discharge was obtained by estimates based on current velocity/area ratio. The discharge measured during sampling was reported using a seven-range scale (<0.5; 0.5–1; 1–5; 5–10; 10–50; >50; >100 L sec-1).
Diatom and water samples for physical and chemical analyses were collected simultaneously. Diatoms were collected with a toothbrush from hard natural substrata (five stones) along a longitudinal transect from the emergence point up to a max distance of 5 m following the European Standard [37]. All samples were preserved with a formaldehyde solution (4% v/v) immediately after sampling. Some typically selected sampling points in each geographic region are reported in Figure 2.
Water samples were collected at the emergence point using 1-L polypropylene and polyethylene bottles and were transported in cold and dark conditions for the laboratory analyses. Further details for the methods used are given in previous studies [24,26,32].

2.3. Laboratory Analyses

Water samples were analysed in a laboratory using standard methods [38,39,40]. Diatom samples were treated to obtain a suspension of clean frustules. Organic matter was eliminated by an oxidation process on a heating plate using hydrogen peroxide (H2O2) and hydrochloric acid (HCl) was added to remove carbonates [41]. After rinsing with distilled water, permanent slides were mounted using Styrax (refractive index = 1.59) and Naphrax (refractive index = 1.73) mounting medium. Diatom observations and counts (~400 valves/frustules for each slide) were performed using a light microscope equipped with phase contrast and micrometric scale at ×1000 magnification. Diatom species were identified according to classic European floristic books [20,42,43,44,45,46,47], and recent taxonomic literature e.g., [48,49,50,51].

2.4. Data Processing and Statistical Analyses

The structure of diatom assemblages was examined by species richness (R), Shannon-Wiener diversity index (H’) [52] and Pielou’s evenness index (J’) [53], calculated using OMNIDIA 6.0 software [54]. Differences in structural indices between islands and the mainland were evaluated by a paired t-test [55,56] and after verification of normal distribution with the Shapiro-Wilktest [57] using R 3.6.1 [58]. For this test probability (p) <0.05 was considered significant.
All diatom data were converted in relative abundances (RA) for the statistical analyses. To minimize the influence of rare species, only the species with RA ≥ 5% in at least two sites are included in the analyses.
Ecological preferences of the most abundant diatom species (RA ≥ 5%) for pH, trophic state, salinity and moisture were attributed according to Van Dam et al. [59]. The taxa with a wider ecological range were equally placed among the respective autecological levels to obtain a simplified framework.
The similarity among diatom assemblages of the four geographic regions was analysed using a non-metric multidimensional scaling ordination (nMDS) [60]. The ordination was performed on a Bray-Curtis similarity matrix of species data log (x + 1) transformed [61]. The significance of the differences was tested by a one-way analysis of similarities (ANOSIM) [62]. For this analysis, p < 0.05 was considered significant. The ANOSIM pairwise test was also performed for each pair of geographic regions. nMDS and ANOSIM were performed using the software PRIMER 6 [63].
Relationships between diatom species and environmental variables were explored by a Canonical Correspondence Analysis (CCA) [64] after the previous assessment of the length of the gradient (>4) by means of a Detrended Correspondence Analysis (DCA) [65] of diatom data using Canoco 4.5 [66]. The two matrices, for physical and chemical data (except for pH), and for diatom data, were log (x+1) transformed. All canonical axes were used to assess the significant variables through analyses by means of a Monte Carlo test (499 permutations).
The indicator value method (IndVal) [67] was used to identify the most characteristic species of three groups defined by the CCA. Indicator values were calculated using abundance and frequency of occurrence of each species in each group as the input data. The species with p < 0.05 based on 499 permutations from the Monte Carlo test were considered to be indicator species. The IndVal was performed by means of the R package indicspecies (ver. 1.7.1).

3. Results

3.1. Environmental Variables

The values of the environmental variables measured and analysed for all springs studied are presented in Table 2. The water temperature ranged from 6.7 to 16.4 °C. The lowest values (6.7–8.7 °C) were measured at BM, AL, UC and TR in Friuli and the highest values (15.1–16.4 °C) at FO, PB and MS in Sardinia and SV and FL in Majorca. pH was slightly alkaline (>7) at several sites and circumneutral (slightly <7) at FP, SV in Majorca, and at VA in Sardinia. More alkaline values of pH (8.24–8.71) were recorded for springs in Friuli and for only one spring (DN) in Sardinia. Conductivity was low-intermediate in all sites (131–1140 µS cm−1). Impacts values (score > 6) were recorded at several springs, mainly in Veneto and Sardinia. Higher shading values (score ≥ 4) were recorded at five springs: FL and SV in Majorca, VA in Sardinia and FC and BR in Veneto. Discharge and current velocity showed very different values among sites. Higher discharge (up to >100 L s−1) and current velocity (score = 4) were generally recorded in the springs of Friuli. The P-PO43− values, available for all sites except for Veneto, were higher at DR (161 µg L−1) and MS (66 µg L−1) in Sardinia. The N-NO3− concentrations were higher for springs in Veneto (1580–5869 µg L−1) and Friuli (2300–6200 µg L−1). Cl− ranged from 0.4 to 195 mg L−1 with lower values and narrower range for springs in Friuli (0.4–4.1 mg L−1) and higher values and a wider range in Majorca (15.7–45.0 mg L−1) and Sardinia (21.3–195 mg L−1). Ca2+ and Mg2+ showed higher values in Veneto (respectively, 52.4–111.5 mg L−1 and 13.6–52.7 mg L−1) and Majorca (respectively, 75.3–97.5 mg L−1 and 13.0–35.1 mg L−1).

3.2. Species Composition and Structure of Diatom Assemblages

The complete list of species was reported in the supplementary table. A total of 176 taxa from 56 genera were found of which 5 (5 genera) were centric and 171 (51 genera) pennate. Ellerbeckia arenaria was the most abundant centric species with a maximum RA = 5.3%. The most rich-species genera were Gomphonema (20), Nitzschia (17), and Navicula (14), followed by Achnanthidium (10), Diploneis (9), and Sellaphora (6). Overall, the species common to samples collected in the islands were 22 and those common to samples collected in the mainland were 11.
The most abundant taxa (RA ≥ 5% in at least one site) were 46 from 25 genera (26% of the total species) (Table 3). Among these, the most frequent species (i.e., occurring in more than six sites) were 13: Achnanthidium lineare, A. minutissimum, A. pyrenaicum, Amphora pediculus, Caloneis fontinalis, Cocconeis euglypta, C. pseudolineata, Denticula tenuis, Humidophila contenta, Meridion circulare, Planothidium frequentissimum, P. lanceolatum and Sellaphora nigri. Considering the most abundant taxa, ecological information was available for 24 taxa (pH and trophic state), 26 taxa (salinity) and 21 taxa (moisture). These species are mainly alkaliphilous (63%), oligohalobous (65%), and linked to the aquatic environment (86%). The trophic requirements highlighted a prevalence of species characteristic of eutrophic and mesotrophic environments, respectively 46% and 29%, in respect of species characteristic of oligotrophic environments (25%).
The taxa found in single sites were 101 from 44 genera (57% of the total species). Among these, the most abundant were Cymbella tridentina, Eunotia arcubus, Fallacia insociabilis, F. muraloides, Gomphonema aff. cymbelliclinum, Hannaea arcus, Nitzschia frustulum, Platessa conspicua and Psammothidium grischunum. The number of taxa found in singles sites was higher in the islands (66) than the mainland (35). The number of rare taxa found in single sites was also higher in the islands (62) than the mainland (24). Distribution, the number of taxa found in single sites and the number of rare taxa (RA < 5%) found in singles sites are presented in Figure 3 and Figure 4.
The results of the structural indices are reported in Table 4. Species richness varied from 8 at VA to 54 at MS, both in Sardinia. The minimum and maximum values of the Shannon-Wiener diversity index (H’) were recorded at FP (1.40) in Majorca and at MS (4.48) in Sardinia. The values of Pielou’s evenness index (J’) ranged from 0.31 at UC, in Friuli to 0.81 at DR, in Sardinia and BA, in Veneto.

3.3. Comparison among Diatom Assemblages

The nMDS ordination showed a clear separation of the diatom assemblages in the springs of Friuli from those of all other regions (Figure 5). The springs RR in Friuli, VA in Sardinia and TO in Veneto were also quite separated from their respective groups. Significant differences among assemblages were confirmed by the ANOSIM test (global R = 0.610; p = 0.001) and pairwise tests (Table 5).
The t-test, performed on the structural indices, indicated significant differences only in diversity as expressed by Shannon-Wiener’s index (t = 2.304; p = 0.031).

3.4. Relationships with Environmental Variables

In the CCA analysis (Figure 6), the first two axes accounted for 60.5% of the total variance of diatom species and environmental data (axis 1: 44.8% and axis 2: 15.7%). Ordination data distinguished three main groups of sites and species. The first group, situated in the right side of the plot includes the sites of the Friuli, with higher discharge and more alkaline waters. This group is composed of Achnanthidium lineare, A. pyrenaicum, Denticula tenuis, Gomphonema elegantissimum, and G. tergestinum. The second group, situated at the upper left part of the plot contains the sites of the islands, characterized by higher temperature and mineralization of water (Cl− and conductivity). This group is composed of Amphora pediculus, Caloneis fontinalis, Cocconeis euglypta, C. pseudolineata, Meridion circulare, Humidophila contenta, Kolbesia gessneri, Planothidium frequentissimum, and Pseudostaurosira alvareziae. The third group, situated at the bottom left part of the plot, include the sites of the Veneto, characterized by higher shading and impacts. This group is composed of Amphora inariensis, Sellaphora nigri, and S. seminulum. Overall, the significant variables for diatom assemblages were discharge (p = 0.018) and Cl− (p = 0.002).
The IndVAL analysis identified a total of eight indicator species for the three distinct groups of sites based on the CCA analysis (Figure 7, Table 6).

4. Discussion

4.1. Main Characteristics of the Diatom Assemblages

Overall, the taxa abundant and common to more than six sites were Achnanthidium lineare, A. minutissimum, A. pyrenaicum, Amphora pediculus, Caloneis fontinalis, Cocconeis euglypta, C. pseudolineata, Denticula tenuis, Humidophila contenta, Meridion circulare, Planothidium frequentissimum, P. lanceolatum, and Sellaphora nigri. Most of these taxa were found in karst springs of different geographic areas in Europe, such as France [68], Poland [69], Austria [70], Slovenia [21], and Bosnia and Herzegovina [23]. By contrast, a large number of taxa occurred in single sites and most of them were rare (RA < 5%). Our observations are supported by previous studies on springs in Italy [32,36] and France [68]. The high number of taxa found in single sites suggests a general marked heterogeneity of assemblages but also a specificity of each site at a regional level, since these species often have high affinity with specific abiotic conditions [32]. This aspect also underlines the role of the springs as refuges for many different species as recently reported by Taxböck et al. [7]. Several species with a strong affinity to spring environments or crenophilous species were found: Caloneis fontinalis, Cymbella tridentina, Diploneis krammeri Lange-Bertalot & Reichardt, Eunotia arcubus, Eunotia minor (Kützing) Grunow, Geissleria gereckei Cantonati & Lange-Bertalot, Gomphonema elegantissimum, Gomphonema lateripunctatum Reichardt & Lange-Bertalot, Meridion circulare, and Odontium mesodon. Most of them were recorded in karst springs both in high and low-altitude, with different disturbance levels e.g., [18,71,72]. Further, Meridion circulare and Odontium mesodon constitute a common association of species of karst springs of different geographic areas. In our study, the crenophilous species were generally more abundant in springs with low impact: C. fontinalis at SG and BR (score = 3–5), C. tridentina at AL (score = 2), E. arcubus at SV (score = 3) G. elegantissimum at BM (score = 2), M. circulare at UL, VA and RR (score = 3) and O. mesodon at RR (score = 3). Cymbella tridentina, found only at AL (1249 m a.s.l.) in Friuli, seems to be the only crenophilous species with a distribution restricted to high-altitude springs in good accordance with previous studies that reported this species from the uppermost sections of carbonate spring-fed streams of the Alps at elevations >1200 m a.s.l. [18,73]. Many taxa are not closely linked to the aquatic environment according to Van Dam et al. [59]. Examples of the pseudaerial/euaerial species, belonging to categories 4+5, (i.e., species occurring in temporarily dry places or outside water bodies) were Achnanthes coarctata (Brébisson ex W. Smith) Grunow, Cosmioneis pusilla (W. Smith) D.G. Mann & Stickle, Denticula subtilis Grunow, Humidophila gallica (W. Smith) Lowe, Kociolek, Q. You, Q. Wang & Stepanek, Diploneis oblongella (Nägeli ex Kützing) Cleve-Euler, Ellerbeckia arenaria, Fallacia insociabilis, Eunotia arcubus, Humidophila contenta and Nitzschia valdestriata Aleem & Hustedt. Whilst being absent in the springs of Friuli, characterized by high discharge, these species were observed mainly in springs with low discharge in Sardinia and Veneto. In these sites, these species can find suitable parts of substrata, temporarily uncovered by the water film due to the low discharge and water abstraction. Their presence underlines the role of springs as multiple ecotones, in particular, the land-water transition.
The diatom flora of the springs studied was composed by a small group of centric species rarely observed in spring environments: Aulacoseira sp., Cyclotella meneghiniana Kützing, Ellerbeckia arenaria, Melosira varians C. Agardh, and Pleurosira laevis (Ehrenberg) Compère. These centric diatoms were found only in the Sardinian springs, generally with very low abundance, and could be favoured by low current velocity. Similar findings were reported in previous studies on thermo-mineral springs of the island [74,75].

4.2. Structure of the Diatom Assemblages

In general, the springs studied hosted rich and diversified diatom communities with a balanced distribution of the species as expressed by species richness (R), diversity (H’), and evenness (J’) indices. High species richness (≥25 taxa) was found at several sites, in Sardinia and Friuli, both with higher naturalness and strongly modified by water abstraction systems. In fact, the coexistence of numerous species seems to depend on a complex combination of local abiotic factors and contrastingly disturbed microhabitats, besides the degree of naturalness of the site [8,32]. The values of species richness were comparable among spring types, i.e., R-L vs LR, except for the spring SG, which showed a higher number of species. This is in agreement with an extensive study carried out on different spring types in the southeastern Alps [76].
The t-test performed on structural indices revealed significant differences between islands and the mainland only in diversity as expressed by the Shannon index (H’).

4.3. Comparison of Diatom Assemblages from Islands and Mainland

The nMDS ordination and ANOSIM test performed on the most abundant taxa did not show a real distinctiveness among diatom assemblages from islands and the mainland. However, the diatom assemblages from the sites of Friuli formed a well-separated cluster and revealed high dissimilarity with those from all other regions. Differences in this group of springs seem to be attributable to specific local factors, such as hydrological features (high discharge and current velocity) and lower mineralization level of the water. In fact, these factors are considered among the most important drivers for the growth and distribution of diatom species [17]. In the nMDS plot, the springs RR, VA and TO were also quite separated from their groups, respectively in Friuli, Sardinia, and Veneto. In Sardinia, VA was the only spring with a seasonal regime. The water flow permanence was an important driver for diatom flora in this site as reported in a previous study [32]. Further, this factor was indicated as a relevant local hydrological factor for diatom assemblages in springs [8,77]. By contrast, differences for the springs RR and TO seem less clear and are probably due to a combination of different factors.
Considering the whole data set, we found further differences in diatom assemblages from the islands and the mainland. For example, the number of taxa found in single sites and the number of rare taxa (also found in single sites) were higher in springs of the islands than those of the mainland. This suggests that the role of springs as refuges for rare species may be greater in the islands and that the condition of geographic insularity may be important in maintaining a high level of biodiversity. In fact, a limited number of sites and low abundance were indicated as characteristics of rare diatom species [78]. In addition, springs are considered in themselves water islands with specific and well-differentiated biocoenoses due to their scattered distribution [8,9,10,11,12,13,14,15,16,17,18,19] and the presence of species with restricted distributions are expected in more isolated freshwater systems because of geographic, physical, and chemical barriers [8].

4.4. Relationships Diatoms-Environmental Variables

Among the environmental variables measured and analysed in this study, discharge and Cl−, respectively associated with the sites of the islands and Friuli, explained the significant amount of variance in diatom assemblages according to CCA analysis. Our results seem to reflect well the characteristics of these sites. For example, water mineralization was reported as an important feature of the springs in Majorca Island [26,79]. In addition, the low discharge was indicated among the most important environmental variables determining the diatom assemblages in low-elevation carbonate springs [72].
The IndVal analysis revealed indicator species quite consistent with the two major gradients highlighted by the CCA analysis. The indicator taxa for the springs of Friuli (group 1) were Achnanthidium pyrenaicum, A. lineare, Gomphonema elegantissimum and G. tergestinum. The first two species of Achnanthidium genus were observed in all sites of Friuli. These species are reported as well-adapted at high discharge and current velocity in the literature e.g., [45,80]. Gomphonema tergestinum is also reported as a species present in fast-flowing running waters [45] while G. elegantissimum seems to prefer moderate current velocity [80] but in these springs may be favoured by the availability of nutrients. For the springs of the islands (group 2) Cocconeis euglypta was the only indicator species. This species seems to be typical of flowing waters with medium-high mineralization [81,82]. However, habitat and ecology of this species are not yet well known because it was not consistently identified over time e.g., [45]. The springs of Veneto (group 3) were characterized by Amphora inariensis, Sellaphora seminulum and S. nigri. Amphora inariensis was found in undisturbed and sometimes moderately impacted freshwater environments [45]. Sellaphora seminulum was reported in springs and intermittently-drying small water bodies but the ecology of this species is not still precisely known [45]. Instead, S. nigri seems to be related to conditions of greater anthropogenic impact and was indicated as a well-known indicator of nutrient enrichment [51]. This species was reported in low-elevation carbonate springs affected by anthropogenic disturbance [72]. However, NGS data suggest great genetic diversity, and the consequent likely necessity of future revisions of its taxonomy and ecology [45].

5. Conclusions

This study provides the first comparison of diatom biodiversity in karst springs with different anthropogenic impact levels in two Mediterranean geographic areas with contrasting characteristics (islands and mainland). In general, the biotic integrity level was good both in sites with greater naturalness and in those strongly modified by water abstraction systems. However, crenophilous species were generally present and abundant in less impacted springs. The diatom flora was very heterogeneous and significant differences were found in the diversity of assemblages between islands and the mainland. The springs in the islands also recorded the highest number of species observed in single sites, as well as of rare species in single sites. Based on the most abundant species, significant differences were observed according to our initial hypothesis. Discharge and Cl− were two important variables in the springs studied and significant factors influencing diatom assemblages. Our results confirm the role of springs as refuges for rare species and suggest that the geographic insularity may be an important factor in maintaining diatom biodiversity. However, Mediterranean karst springs are fragile ecosystems increasingly exposed to natural and anthropogenic pressures. Morphological alterations and uncontrolled water abstraction as a result of inadequate management may cause a reduction of their ecological integrity and a loss of biodiversity in the future.

Supplementary Materials

The supplementary materials are available online at https://www.mdpi.com/2073-4441/11/12/2602/s1.

Author Contributions

G.G.L. and M.C. conceived and designed the paper. G.G.L., S.B., C.D., and M.C. carried out sampling and identification of the diatom species. R.Z. and E.P. provided environmental data. B.M.P. performed all statistical analyses and figures. G.G.L. wrote the paper and prepared tables and figures. G.G.L., B.M.P., C.D., A.L. and M.C. revised and edited the manuscript.

Funding

This research received no external funding. The activities of Prof. Antonella Lugliè were supported by the research fund of the University of Sassari (fondo di Ateneo per la ricerca 2019).

Acknowledgments

We thank Giuseppe Cossu, Salvatore Sulis (Sardinia Forestry Corps) and EN.A.S. (Ente Acque della Sardegna) for the field inspections, Bruno Giordanu and Maria Antonietta Mariani for their support in sampling and the colleagues of the University of Sassari (Laboratory of Aquatic Ecology) for the analysis of water samples. We are also grateful to the colleagues of ARPA FVG (Regional Agency for Environmental Protection of Friuli-Venezia Giulia), in particular Ivan Martinuzzi for the analysis of water samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Parise, M.; Gabrovsek, F.; Kaufmann, G.; Ravbar, N. Recent advances in karst research: From theory to fieldwork and applications. In Advances in Karst Research: Theory, Fieldwork and Application; Parise, M., Gabrovsek, F., Kaufmann, G., Ravbar, N., Eds.; Geological Society: London, UK, 2018; Special Publications 466; pp. 1–24. [Google Scholar]
  2. Fiorillo, F.; Stevanović, Z. Introductory editorial thematic issue: Mediterranean karst hydrogeology. Environ. Earth Sci. 2015, 74, 1–3. [Google Scholar] [CrossRef] [Green Version]
  3. Bakalowicz, M.; Fleury, P.; Dörfliger, N.; Seidel, J.L. Coastal karst aquifers in Mediterranean regions. A valuable ground water resource in complex aquifers. In Proceedings of the Technologia de la Intrusion de Agua de Mar en Acuiferos Costeros: Paises Mediterraneos (TIAC ’03), Alicante, Spain, 11–14 March 2003; López Geta, J.A., de Dios Gómez, J., de la Orden, V., Ramos, G., Rodríguez, L., Eds.; Instituto Geológico y Minero de España: Madrid, Spain, 2003; pp. 125–128. [Google Scholar]
  4. Bonacci, O. Karst landscape ecohydrology. In Proceedings of the International Symposium on Water Management and Hydraulic Engineering, Ohrid, Macedonia, 1–5 September 2009; pp. 781–790. [Google Scholar]
  5. Zelnik, I.; Balanč, T.; Toman, M.J. Diversity and Structure of the tychoplankton diatom community in the limnocrene spring Zelenci (Slovenia) in relation to environmental factors. Water 2018, 10, 361. [Google Scholar] [CrossRef] [Green Version]
  6. Maiolini, B.; Carolli, M.; Silveri, L. Ephemeroptera, Plecoptera and Trichoptera in springs in Trentino (south-eastern Alps). J. Limnol. 2011, 70 (Suppl. 1), 122–133. [Google Scholar] [CrossRef]
  7. Taxböck, L.; Linder, H.P.; Cantonati, M. To what extent are Swiss springs refugial habitats for sensitive and endangered diatom taxa? Water 2017, 9, 967. [Google Scholar] [CrossRef] [Green Version]
  8. Cantonati, M.; Füreder, L.; Gerecke, R.; Jüttner, I.; Cox, E.J. Crenic habitats, hotspots for freshwater biodiversity conservation: Toward an understanding of their ecology. Freshw. Sci. 2012, 312, 463–480. [Google Scholar] [CrossRef]
  9. Meyer, A.; Meyer, E.I.; Meyer, C. Lotic communities of two small temporary karstic stream systems (East Westphalia, Germany) along a longitudinal gradient of hydrological intermittency. Limnologica 2003, 33, 271–279. [Google Scholar] [CrossRef] [Green Version]
  10. Daly, D.; Dassargues, A.; Drew, D.; Dunne, S.; Goldscheider, N.; Neale, S.; Popescu, I.C.; Zwahlen, F. Main concepts of the “European approach” to karst-groundwater-vulnerability assessment and mapping. Hydrogeol. J. 2002, 10, 340–345. [Google Scholar] [CrossRef] [Green Version]
  11. Bolle, H.J. Mediterranean Climate: Variability and Trends, 1st ed.; Springer: Berlin, Germany, 2003; pp. 1–372. [Google Scholar]
  12. Intergovernmental Panel on Climate Change—IPCC. Climate Change 2007: The Physical Science Basis; Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate, Change; Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K.B., Tignor, M., Miller, H.L., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2007; pp. 1–996. [Google Scholar]
  13. Spaɫek, K.; Procków, J. New, rare and constant habitats for endangered aquatic plant communities: The importance of microhabitats for global biodiversity. Diversity 2019, 11, 96. [Google Scholar] [CrossRef] [Green Version]
  14. Herman, J.W.; Culver, D.C.; Salzman, J. Groundwater ecosystems and the service of water purification. Stan. Environ. Law J. 2001, 20, 479–495. [Google Scholar]
  15. Goldscheider, N. A holistic approach to groundwater protection and ecosystem services in karst terrains. AQUA Mundi 2012, Am06046, 117–124. [Google Scholar] [CrossRef]
  16. Poulíčková, A.; Hájková, P.; Křenková, P.; Hájek, M. Distribution of diatoms and bryophytes on linear transects through spring fens. Nova Hedwigia 2004, 78, 411–424. [Google Scholar] [CrossRef]
  17. Cantonati, M.; Gerecke, R.; Bertuzzi, E. Springs of the Alps–sensitive ecosystems to environmental change: From biodiversity assessments to long-term studies. Hydrobiologia 2006, 562, 59–96. [Google Scholar] [CrossRef]
  18. Mogna, M.; Cantonati, M.; Andreucci, F.; Angeli, N.; Berta, G.; Miserere, L. Diatom communities and vegetation of springs in the south-western Alps. Acta Bot. Croat. 2015, 742, 265–285. [Google Scholar] [CrossRef] [Green Version]
  19. Werum, M. Die Kieselalgengesellschaften in Quellen. Abhängigkeit von Geologie und anthropogener Beeinflussung in Hessen (Bundesrepublik Deutschland); Hessisches Landesamt fur Umwelt und Geologie: Wiesbaden, Germany, 2001; p. 273.
  20. Rott, E.; Pfister, P.; Pipp, E.; Pall, K.; Binder, N.; Ortler, K. Indikationslisten für Aufwuchsalgen in Fließgewässern. Teil 2; Bundesministerium für Land-u. Forstwirtschaft, Wasserwirtschaftskataster: Wien, Austria, 1999; p. 248. [Google Scholar]
  21. Menegalija, T.; Kosi, G. Distribution of diatoms in springs in the Julian Alps (NW Slovenia). Nat. Slov. 2008, 101, 21–37. [Google Scholar]
  22. Žutinić, P.; Petrić, I.; Gottstein, S.; Udovič, M.G.; Borojević, K.K.; Kamberović, J.; Kolda, A.; Plenković-Moraj, A.; Ternjej, I. Microbial mats as shelter microhabitat for amphipods in an intermittent karstic spring. Knowl. Manag. Aquat. Ecosyst. 2018, 419, 7. [Google Scholar] [CrossRef]
  23. Dedić, A.; Plenković-Moraj, A.; Borojević, K.K.; Hafner, D. The first report on periphytic diatoms on artificial and natural substrate in the karstic spring Bunica, Bosnia and Herzegovina. Acta Bot. Croat. 2015, 742, 393–406. [Google Scholar] [CrossRef] [Green Version]
  24. Angeli, N.; Cantonati, M.; Spitale, D.; Lange-Bertalot, H. A comparison between diatom assemblages in two groups of carbonate, low-altitude springs with different levels of anthropogenic disturbances. Fottea 2010, 101, 115–128. [Google Scholar] [CrossRef] [Green Version]
  25. Cantonati, M.; Bertuzzi, E.; Scalfi, A. CRENODAT (biodiversity assessment and integrity evaluation of springs of Trentino (Italian Alps) and long-term ecological research): Project design and preliminary results. In Integrated Watershed Management: Perspectives and Problems; Beheim, E., Rajwar, G.S., Haigh, M., Krecek., J., Eds.; Springer, Capital Publishing Company: New Delhi, India, 2010; pp. 121–132. [Google Scholar]
  26. Delgado, C.; Ector, L.; Novais, M.H.; Blanco, S.; Hoffmann, L.; Pardo, I. Epilithic diatoms of springs and spring–fed streams in Majorca Island (Spain) with the description of a new diatom species Cymbopleura margalefii sp. nov. Fottea 2013, 132, 87–104. [Google Scholar] [CrossRef] [Green Version]
  27. Dell’Uomo, A. Chorological analysis of the algal flora of “Su Gologone” spring in Sardinia. In Proceedings of the International Symposium on Biogeographical Aspects of Insularity, Rome, Italy, 18–22 May 1987; Accademia Nazionale dei Lincei: Rome, Italy, 1990; pp. 417–427. [Google Scholar]
  28. Lai, G.G.; Padedda, B.M.; Wetzel, C.E.; Lugliè, A.; Sechi, N.; Ector, L. Valutazione della qualità ambientale della sorgente carsica Su Gologone (Sardegna centro-orientale) mediante l’impiego delle diatomee epilitiche: Primi risultati. Biol. Ambient. 2015, 291, 37–45. [Google Scholar]
  29. Lai, G.G.; Padedda, B.M.; Wetzel, C.E.; Lugliè, A.; Sechi, N.; Ector, L. Epilithic diatom assemblages and environmental quality of the Su Gologone karst spring (central-eastern Sardinia, Italy). Acta Bot. Croat. 2016, 751, 129–143. [Google Scholar] [CrossRef] [Green Version]
  30. Lai, G.G.; Ector, L.; Lugliè, A.; Sechi, N.; Wetzel, C.E. Sellaphora gologonica sp. nov. (Bacillariophyta, Sellaphoraceae), a new diatom species from a Mediterranean karst spring (Sardinia, Italy). Phytotaxa 2018, 3562, 145–157. [Google Scholar] [CrossRef]
  31. Lai, G.G.; Ector, L.; Wetzel, C.E.; Sechi, N.; Lugliè, A.; Padedda, B.M. Periphytic diatoms of the Mediterranean karst spring Sa Vena (Su Gologone system, Sardinia, Italy): Relationships with environmental variables and effects of an extreme flash flood. Inland Waters 2018, 83, 284–293. [Google Scholar] [CrossRef]
  32. Lai, G.G.; Padedda, B.M.; Ector, L.; Wetzel, C.E.; Lugliè, A.; Cantonati, M. Mediterranean karst springs: Diatom biodiversity hotspots under the pressure of hydrological fluctuation and nutrient enrichment. Plant Biosyst. 2019. [Google Scholar] [CrossRef]
  33. Llobera, M.; Ferriol, A. Seasonal fluctuations in the chemistry of limestone springs from the Tramuntana Mountain Range (Mallorca, Spain). Verh. Int. Ver. Limnol. 1994, 25, 1445–1448. [Google Scholar] [CrossRef]
  34. Available online: http://sgi.isprambiente.it/GeositiWeb/default.aspx (accessed on 10 October 2019).
  35. Poldini, L.; Oriolo, G.; Vidali, M.; Tomasella, M.; Stoch, F.; Orel, G. Manuale Degli Habitat del Friuli Venezia Giulia. Strumento a Supporto Della Valutazione D’impatto Ambientale (VIA), Ambientale Strategica (VAS) e di Incidenza Ecologica (VIEc). Regione Autonoma Friuli Venezia Giulia—Direzione Centrale Ambiente e Lavori Pubblici—Serv. VIA, Univ. Degli STUDI di Trieste—Dip. Biologia. Available online: http://www.regione.fvg.it/ambiente/manuale/home.htm (accessed on 30 July 2019).
  36. Cantonati, M.; Spitale, D. The role of environmental variables in structuring epiphytic and epilithic diatom assemblages in springs and streams of the Dolomiti Bellunesi National Park (south-eastern Alps). Fundam. Appl. Limnol. Arch. Hydrobiol. 2009, 1742, 117–133. [Google Scholar] [CrossRef]
  37. CEN. Water Quality—Guidance Standard for the Routine Sampling and Pretreatment of Benthic Diatoms from Rivers; European Standard EN 13946; European Committee for Standardization: Brussels, Belgium, 2003; pp. 1–14. [Google Scholar]
  38. Agenzia per la Protezione Dell’ambiente e per i Servizi Tecnici/Istituto di Ricerca Sulle Acque—Consiglio Nazionale delle Ricerche. Metodi Analitici per le Acque; Rapporti 29/2003; I.G.E.R. srl. APAT: Rome, Italy, 2003; ISBN 88-448-0083-7. Available online: http://www.isprambiente.gov.it/it/pubblicazioni/manuali-e-linee-guida/metodi-analitici-per-le-acque (accessed on 7 October 2017).
  39. APHA. Standard Methods for the Examination of Water and Wastewater, 20th ed.; APHA, AWWA & WEF, American Public Health Association: Washington, DC, USA, 2000; pp. 1–1220. [Google Scholar]
  40. Strickland, J.D.H.; Parsons, T.R. A Practical Handbook of Seawater Analysis, Bulletin 167, 2nd ed.; Fisheries Research Board of Canada: Ottawa, ON, Canada, 1972; pp. 1–328. [Google Scholar]
  41. Istituto Superiore per la Protezione e la Ricerca Ambientale. Metodi Biologici per le Acque Superficali Interne; Manuali e Linee Guida 111/2014; Istituto Superiore per la Protezione e la Ricerca Ambientale: Rome, Italy, 2003. Available online: http://admin.isprambiente.gov.it/it/pubblicazioni/manuali-e-linee-guida/metodi-biologici-per-le-acque-superficiali-interne (accessed on 10 June 2017).
  42. Krammer, K.; Lange-Bertalot, H. Bacillariophyceae. Teil 1–4; Gustav Fischer Verlag: Stuttgart, Germany; New York, NY, USA, 1991. [Google Scholar]
  43. Krammer, K.; Lange-Bertalot, H. Bacillariophyceae. English and French Translation of the Keys Teil 5; Spektrum Akademischer Verlag: Berlin/Heidelberg, Germany, 2000; pp. 1–311. [Google Scholar]
  44. Lange-Bertalot, H.; Cavacini, P.; Tagliaventi, N.; Alfinito, S. Diatoms of Sardinia. Rare and 76 New Species in Rock Pools and Other Ephemeral Waters; Iconographia Diatomologica 12; A.R.G. Gantner Verlag, KG: Ruggell, Germany, 2003; pp. 1–438. [Google Scholar]
  45. Lange-Bertalot, H.; Hofmann, G.; Werum, M.; Cantonati, M. Freshwater Benthic Diatoms of Central Europe: Over 800 Common Species Used in Ecological Assessment; English edition with updated taxonomy and added species; Koeltz Botanical Books: Schmitten, Oberreifenberg, Germany, 2017; pp. 1–942. [Google Scholar]
  46. Wojtal, A.Z. Species Composition and Distribution of Diatom Assemblages in Spring Waters from Various Geological Formations in Southern Poland; Bibliotheca Diatomologica 59; J. Cramer: Stuttgart, Germany, 2013; pp. 1–436. [Google Scholar]
  47. Żelazna-Wieczorek, J. Diatom Flora in Springs of Lódź Hills (Central Poland). Biodiversity, Taxonomy, and Temporal Changes of Episammic Diatom Assemblages in Springs Affected by Human Impact, Diatom Monographs 13; A.R.G. Gantner: Ruggell, Liechtenstein, Germany, 2011; pp. 1–419. [Google Scholar]
  48. Cantonati, M.; Lange-Bertalot, H. Geissleria gereckei sp. nov. (Bacillariophyta) from leaf-litter covered stones of very shaded carbonate mountain springs with extremely low discharge. Phycol. Res. 2009, 57, 171–177. [Google Scholar] [CrossRef]
  49. Novais, M.H.; Blanco, S.; Hlúbiková, D.; Falasco, E.; Gomà, J.; Delgado, C.; Ivanov, P.; Àcs, E.; Morais, M.; Hoffmann, L.; et al. Taxonomic revision and biogeography of the Gomphonema rosenstockianum and G. tergestinum species complex. Fottea 2009, 9, 257–274. [Google Scholar] [CrossRef] [Green Version]
  50. Cejudo-Figueiras, C.; Morales, E.A.; Wetzel, C.E.; Blanco, S.; Hoffmann, L.; Ector, L. Analysis of the type of Fragilaria construens var. subsalina (Bacillariophyceae) and description of two morphologically related taxa from Europe and the United States. Phycologia 2011, 501, 67–77. [Google Scholar] [CrossRef]
  51. Wetzel, C.E.; Ector, L.; Van de Vijver, B.; Compère, P.; Mann, D.G. Morphology, typification and critical analysis of some ecologically important small naviculoid species (Bacillariophyta). Fottea 2015, 152, 203–234. [Google Scholar] [CrossRef] [Green Version]
  52. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication, 1st ed.; University of Illinois Press: Urbana, IL, USA, 1949; pp. 1–117. [Google Scholar]
  53. Pielou, E.C. Ecological Diversity; Wiley: New York, NY, USA, 1975; pp. 1–165. [Google Scholar]
  54. Lecointe, C.; Coste, M.; Prygiel, J. “Omnidia”: Software for taxonomy, calculation of diatom indices and inventories management. Hydrobiologia 1993, 269, 509–513. [Google Scholar] [CrossRef]
  55. Student. The probable error of a mean. Biometrika 1908, 61, 1–25. [Google Scholar]
  56. Fisher, R.A. Applications of “Student’s’ distribution”. Metron 1925, 5, 90–104. [Google Scholar]
  57. Shapiro, S.S.; Wilk, M.B. An analysis of variance test for normality (complete samples). Biometrika 1965, 52, 591–611. [Google Scholar] [CrossRef]
  58. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: http://www.r-project.org/index.html (accessed on 20 September 2019).
  59. Van Dam, H.; Mertens, A.; Sinkeldam, J. A coded checklist and ecological indicator values of freshwater diatoms from The Netherlands. Netherland J. Aquat. Ecol. 1994, 28, 117–133. [Google Scholar]
  60. Kruskal, J.B. Nonmetric multidimensional scaling: A numerical method. Psychometrika 1964, 29, 115–129. [Google Scholar] [CrossRef]
  61. Bray, J.R.; Curtis, J.T. An ordination of the upland forest communities of Southern Wisconsin. Ecol. Monogr. 1957, 27, 325–349. [Google Scholar] [CrossRef]
  62. Clarke, K.R.; Warwick, R.M. Change in Marine Communities: An Approach to Statistical Analysis and Interpretation, 1st ed.; Plymouth Marine Laboratory: Plymouth, UK, 1994; pp. 1–144. [Google Scholar]
  63. Clarke, K.R.; Gorley, R.N. PRIMER v6: User Manual/Tutorial (Plymouth Routines in Multivariate Ecological Research); Plymouth Marine Laboratory: Plymouth, UK, 2006; pp. 1–190. [Google Scholar]
  64. Ter Braak, C.J.F. Canonical correspondence analysis: A new eigenvector technique for multivariate direct gradient analysis. Ecology 1986, 67, 1167–1179. [Google Scholar] [CrossRef] [Green Version]
  65. Hill, M.O.; Gauch, H.G., Jr. Detrended correspondence analysis: An improved ordination technique. Vegetatio 1980, 42, 47–58. [Google Scholar] [CrossRef]
  66. Ter Braak, C.J.F.; Šmilauer, P. CANOCO Reference Manual and Cano Draw for Windows User’s Guide: Software for Canonical Community Ordination (Version 4.5); Microcomputer Power: Ithaca, NY, USA, 2002; pp. 1–500. [Google Scholar]
  67. Dufrêne, M.; Legendre, P. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol. Monogr. 1997, 673, 345–366. [Google Scholar] [CrossRef]
  68. Bertrand, J.; Renon, J.-P.; Monnier, O. Les diatomées des sources du rebord karstique de la Beauce de la région orléanaise. Symbioses 1999, 1, 3–14. [Google Scholar]
  69. Wojtal, A.Z.; Sobczyk, Ł. The influence of substrates and physicochemical factors on the composition of diatom assemblages in karst springs and their applicability in water-quality assessment. Hydrobiologia 2012, 6951, 97–108. [Google Scholar] [CrossRef] [Green Version]
  70. Niedermayr, R.; Schagerl, M. Structuring factors of the phytobenthos community along a mountain headwater (Kalkalpen National Park, Austria). Fundam. Appl. Limnol. 2010, 1772, 93–104. [Google Scholar] [CrossRef] [Green Version]
  71. Cantonati, M.; Lange-Bertalot, H. Diatom biodiversity of springs in the Berchtesgaden National Park (northern Alps, Germany), with the ecological and morphological characterization of two species new to science. Diatom Res. 2010, 252, 251–280. [Google Scholar] [CrossRef]
  72. Cantonati, M.; Angeli, N.; Spitale, D.; Lange–Bertalot, H. A new Navicula (Bacillariophyta) species from low–elevation carbonate springs affected by anthropogenic disturbance. Fottea 2016, 162, 255–265. [Google Scholar] [CrossRef]
  73. Cantonati, M.; Lange-Bertalot, H.; Scalfi, A.; Angeli, N. Cymbella tridentina sp. nov. (Bacillariophyta), a crenophilous diatom from carbonate springs of the Alps. J. N. Am. Benthol. Soc. 2010, 293, 775–788. [Google Scholar] [CrossRef]
  74. Lai, G.G.; Padedda, B.M.; Wetzel, C.E.; Cantonati, M.; Sechi, N.; Lugliè, A.; Ector, L. Diatom assemblages from different substrates of the Casteldoria thermo-mineral spring (Northern Sardinia, Italy). Bot. Lett. 2019, 1661, 14–31. [Google Scholar] [CrossRef]
  75. Lai, G.G.; Beauger, A.; Wetzel, C.E.; Padedda, B.M.; Voldoire, O.; Lugliè, A.; Allain, E.; Ector, L. Diversity, ecology and distribution of benthic diatoms in thermo-mineral springs in Auvergne (France) and Sardinia (Italy). PeerJ 2019, 7, e7238. [Google Scholar] [CrossRef]
  76. Cantonati, M.; Rott, E.; Spitale, D.; Angeli, N.; Komárek, J. Are benthic algae related to spring types? Freshw. Sci. 2012, 31, 481–498. [Google Scholar] [CrossRef]
  77. Gerecke, R.; Cantonati, M.; Spitale, D.; Stur, E.; Wiedenbrug, S. The challenges of long-term ecological research in springs in the northern and southern Alps: Indicator groups, habitat diversity, and medium-term change. J. Limnol. 2011, 70 (Suppl. 1), 168–187. [Google Scholar] [CrossRef]
  78. Noga, T.; Rybak, M. First record of Pinnularia subinterrupta Krammer & Schroeter in Poland—A rare species in Europe. Biodivers. Res. Conserv. 2017, 451, 17–21. [Google Scholar]
  79. Moyá, G.; Ferriol, A.; Llobera, M. Chlorophyll contents in the freshwater vegetation of Serra de Tramuntana springs (Mallorca, Spain). Verh. Int. Ver. Limnol. 1991, 24, 2002–2006. [Google Scholar]
  80. Falasco, E.; Piano, E.; Bona, F. Guida al riconoscimento e all’ecologia delle principali diatomee fluviali dell’Italia nord occidentale. Biol. Ambient. 2013, 27, 1–292. [Google Scholar]
  81. Monnier, O.; Rimet, F.; Bey, M.; Chavaux, R.; Ector, L. Sur l’identité de Cocconeis euglypta Ehrenberg 1854 et C. lineata Ehrenberg 1843. Un approche par le sources historiques. Diatomania 2007, 11, 30–45. [Google Scholar]
  82. Blanco, S.; Cejudo-Figueiras, C.; Álvarez-Blanco, I.; Becares, E.; Hoffmann, L.; Ector, L. Atlas de las Diatomeas de la Cuenca del Duero, 1st ed.; Área de Publicaciones, Universidad de León: Leon, Spain, 2010; pp. 1–382. [Google Scholar]
Water 11 02602 g001 550
Figure 1. Geographic position of the four regions studied in the two Mediterranean areas: (A) Majorca Island (Spain); (B) Sardinia, (C) Veneto, (D) Friuli (Italy).
Figure 1. Geographic position of the four regions studied in the two Mediterranean areas: (A) Majorca Island (Spain); (B) Sardinia, (C) Veneto, (D) Friuli (Italy).
Water 11 02602 g001
Water 11 02602 g002 550
Figure 2. Some typically selected sampling points: (A) Font des Pí; (B) Son Vich; (C) Su Gologone-SVM; (D) S’Ulidone; (E) Buso delle Anguane; (F) La Ferrara; (G) Goriuda; (H) Barman.
Figure 2. Some typically selected sampling points: (A) Font des Pí; (B) Son Vich; (C) Su Gologone-SVM; (D) S’Ulidone; (E) Buso delle Anguane; (F) La Ferrara; (G) Goriuda; (H) Barman.
Water 11 02602 g002
Water 11 02602 g003 550
Figure 3. Distribution and number of taxa found in single sites. All codes of the sites are reported in Table 1.
Figure 3. Distribution and number of taxa found in single sites. All codes of the sites are reported in Table 1.
Water 11 02602 g003
Water 11 02602 g004 550
Figure 4. Distribution and number of rare taxa (RA < 5%) found in single sites. All codes of the sites are reported in Table 1.
Figure 4. Distribution and number of rare taxa (RA < 5%) found in single sites. All codes of the sites are reported in Table 1.
Water 11 02602 g004
Water 11 02602 g005 550
Figure 5. Non-metric multidimensional scaling (nMDS) ordination plot for diatom assemblages from the springs studied. All codes of the sites are reported in Table 1.
Figure 5. Non-metric multidimensional scaling (nMDS) ordination plot for diatom assemblages from the springs studied. All codes of the sites are reported in Table 1.
Water 11 02602 g005
Water 11 02602 g006 550
Figure 6. Canonical Correspondence Analysis (CCA) ordination plot for the diatom assemblages. Vectors = environmental variables; circles = spring sites, triangles = diatom species. All codes of sites are reported in Table 1. Acronyms of the species: ACLI = Achnanthidium lineare; ADMI = Achnanthidium minutissimum; ADPY = Achnanthidium pyrenaicum; AINA = Amphora inariensis; APED = Amphora pediculus; CFON = Caloneis fontinalis; CEUG = Cocconeis euglypta; COPL = Cocconeis pseudolineata; DTEN = Denticula tenuis; GELG = Gomphonema elegantissimum; GTER = Gomphonema tergestinum; HUCO = Humidophila contenta; KGES = Kolbesia gessneri; MCIR = Meridion circulare; PLFR = Planothidium frequentissimum; PTLA = Planothidium lanceolatum; PALV = Pseudostaurosira alvareziae; SNIG = Sellaphora nigri; SSEM = Sellaphora seminulum.
Figure 6. Canonical Correspondence Analysis (CCA) ordination plot for the diatom assemblages. Vectors = environmental variables; circles = spring sites, triangles = diatom species. All codes of sites are reported in Table 1. Acronyms of the species: ACLI = Achnanthidium lineare; ADMI = Achnanthidium minutissimum; ADPY = Achnanthidium pyrenaicum; AINA = Amphora inariensis; APED = Amphora pediculus; CFON = Caloneis fontinalis; CEUG = Cocconeis euglypta; COPL = Cocconeis pseudolineata; DTEN = Denticula tenuis; GELG = Gomphonema elegantissimum; GTER = Gomphonema tergestinum; HUCO = Humidophila contenta; KGES = Kolbesia gessneri; MCIR = Meridion circulare; PLFR = Planothidium frequentissimum; PTLA = Planothidium lanceolatum; PALV = Pseudostaurosira alvareziae; SNIG = Sellaphora nigri; SSEM = Sellaphora seminulum.
Water 11 02602 g006
Water 11 02602 g007 550
Figure 7. 1–28. Light microscopy (LM). 1–5: Achnanthidium lineare W. Smith; 6–10: Achnanthidium pyrenaicum (Hustedt) Kobayasi; 11: Amphora inariensis Krammer; 12–13: Cocconeis euglypta Ehrenberg; 14–20: Gomphonema elegantissimum Reichardt & Lange-Bertalot; 21–26: Gomphonema tergestinum (Grunow) Fricke; 27: Sellaphora nigri (De Notaris) C.E. Wetzel & Ector; 28: Sellaphora seminulum (Grunow) D.G. Mann. Scale bars = 10 μm.
Figure 7. 1–28. Light microscopy (LM). 1–5: Achnanthidium lineare W. Smith; 6–10: Achnanthidium pyrenaicum (Hustedt) Kobayasi; 11: Amphora inariensis Krammer; 12–13: Cocconeis euglypta Ehrenberg; 14–20: Gomphonema elegantissimum Reichardt & Lange-Bertalot; 21–26: Gomphonema tergestinum (Grunow) Fricke; 27: Sellaphora nigri (De Notaris) C.E. Wetzel & Ector; 28: Sellaphora seminulum (Grunow) D.G. Mann. Scale bars = 10 μm.
Water 11 02602 g007
Table
Table 1. Main characteristics of the springs studied. Type: R = Rheocrenic; R-L = Rheo-limnocrenic; L-R = Limno- rheocrenic. Regime: P = Perennial; S = Seasonal.
Table 1. Main characteristics of the springs studied. Type: R = Rheocrenic; R-L = Rheo-limnocrenic; L-R = Limno- rheocrenic. Regime: P = Perennial; S = Seasonal.
Spring NameCodeAltitude (m.a.s.l.)Latitude (N)Longitude (E)TypeRegimeSpring-Tapping
Majorca
Font des PíFP47139°46′14″2°48′26″RPX
Font de s′ OllaFL5639°45′23″ 2°42′42″RPX
Son VichSV35639°39′12″ 2°32′08″RPX
Sardinia
Fruncu ‘e OcheFO5540°33′84″ 009°40′67″RPX
PubusinuPB21239°24′34″008°31′39″RPX
San Giovanni-DNDN19039°20′13″008°37′37″RPX
San Giovanni-DRDR16840°19′14″009°36′55″RPX
Sa VarvaVA78040°13′39″009°12′31″RS
Sa Vena Manna-SMS27640°50′25″008°48′20″RPX
Su Gologone-SVMSG11640°17′34″009°29′80″L-RP
S′UlidoneUL60640°34′30″009°38′35″RPX
Veneto
BrusaferriBR67545°35′57″11°12′59″R-LPX
Buso delle AnguaneBA37045°33′04″11°17′07″R-LPX
Fonte del CoppoFC18345°28′39″11°06′51″L-RPX
La FerraraLF64045°32′59″11°11′08″L-RPX
PrealbaPR44545°37′34″10°49′47″L-RPX
TorricelleTO22045°28′09″11°00′16″R-LPX
Friuli
TorreTR52846°18′27″13°16′21″RP
UcceaUC1.12246°20′02″13°19′30″RP
BarmanBM65246°20′21″13°17′18″RP
Rio ResearticoRR45446°22′22″13°13′01″RPX
GoriudaGD70146°23′44″13°26′8″RP
AlbaAL1.24946°28′12″13°13′48″RP
Table
Table 2. Values of the environmental variables measured and analysed in all springs studied.
Table 2. Values of the environmental variables measured and analysed in all springs studied.
MajorcaSardiniaVenetoFriuli
Variable/SpringFLFPSVFOPBDNDRVAMSSGULFCBAPRTOLFBRTRUCBMRRGDAL
T [°C]15.712.515.115.115.914.213.413.216.413.09.51312.211.91212.310.38.78.36.710.813.17.3
pH [units]7.256.656.917.627.808.107.746.987.497.527.207.837.97.717.67.27.28.328.248.718.388.478.40
Conductivity [μS cm−1]6124167963395714784397601140326357424333438436498578154170136148131229
Impacts [score]41387893933105103105752332
Shading [score]52432324133522234113131
Discharge [L s−1]10–500.5–15–100.5–1<0.5<0.5<0.50.5–1<0.51–5<0.5<0.5<0.50.5–1<0.5<0.50.02>100 10–50>1005–100.5–1>50
Current velocity [score]42332223233322232444414
P-PO43− [µg L−1]1111478161866109------010010200
N-NO3− [µg L−1]22513423137146030920931549363480649225713542257564358691580300023002300230062003000
Cl− [mg L−1]24.115.745.028.456.760.339.042.5195.021.346.07.815.09.911.715.06.01.00.50.40.64.10.5
Ca2+ [mg L−1]97.575.385.830.042.032.048.056.054.032.018.0-52.4--106.2111.524.435.823.721.325.422.0
Mg2+ [mg L−1]22.513.035.118.227.921.918.251.048.619.418.2-13.6--52.715.77.41.04.68.92.212.8
Table
Table 3. Floristic list and relative abundance of the most abundant diatom taxa (RA ≥ 5% in at least one site).
Table 3. Floristic list and relative abundance of the most abundant diatom taxa (RA ≥ 5% in at least one site).
MajorcaSardiniaVenetoFriuli
Taxa FPFLSVFOPBDNDRVAMSSGULBRBAFCLFPRTOTRUCBMRRGDAL
Achnanthidium exile (Kützing) Heiberg9.2
Achnanthidium lineare W. Smith 1.7 0.9 0.6 0.8 7.72.40.50.55.47.0
Achnanthidium minutissimum (Kützing) Czarnecki 76.67.964.616.435.131.25.34.30.93.61.043.39.831.469.215.7 4.4
Achnanthidium pyrenaicum (Hustedt) Kobayasi 6.1 2.5 7.2 52.674.930.06.363.422.9
Amphora inariensis Krammer 0.6 7.812.340.21.61.6
Amphora indistincta Levkov 1.0 0.53.0 0.711.12.9
Amphora pediculus (Kützing) Grunow 2.127.01.26.9 1.517.60.715.321.717.910.711.97.93.9 0.21.00.50.5
Caloneis fontinalis (Grunow) Cleve-Euler 0.7 0.7 3.96.0 5.6 1.0 0.8
Cocconeis euglypta Ehrenberg 0.914.3 7.76.17.5 12.111.11.7 0.2
Cocconeis pseudolineata (Geitler) Lange-Bertalot 1.646.2 0.57.20.59.4 0.51.0
Cymbella affinis Kützing 1.4 1.59.8
Cymbella tridentina Lange-Bertalot, M. Cantonati & A. Scalfi 32.3
Derticula tenuis Kützing0.2 8.0 0.57.00.52.01.5
Diploneis separanda Lange-Bertalot0.2 6.3 0.52.1 0.4
Ellerbeckia arenaria (Moore ex Ralfs) R.M. Crawford 5.3 0.71.5
Encyonema minutum (Hilse) D.G. Mann 0.3 1.51.57.3
Encyonema ventricosum (Agardh) Grunow 0.511.1 2.04.4
Eunotia arcubus Nörpel & Lange-Berlalot 15.3
Fallacia insociabilis (Krasske) D.G. Mann 0.2 20.0
Fallacia muraloides (Hustedt) D.G. Mann 5.7
Fallacia sublucidula (Hustedt) D.G. Mann 0.25.3
Gomphonema aff. cymbelliclinum Reichardt & Lange-Bertalot 30.2
Gomphonema elegantissimum Reichardt & Lange-Bertalot 26.85.847.5 3.0
Gomphonema micropus Kützing 0.56.2 1.2 1.5 0.2
Gomphonema rosenstockianum Lange-Bertalot & Reichardt 0.9 3.7 18.6
Gomphonema sp.2 12.0
Gomphonema subclavatum Grunow 0.5 0.2 1.517.4
Gomphonema tergestinum (Grunow) Fricke 7.21.9 6.3
Gomphosphenia lingulatiformis (Lange-Bertalot & Reichardt) Lange-Bertalot 32.70.6
Hannaea arcus (Ehrenberg) R.M. Patrick 10.2
Humidophila contenta (Grunow) Lowe, Kociolek et al. 3.7 0.2 0.2 5.9 3.6 2.0 1.657.3
Kolbesia gessneri (Hustedt) Aboal 9.6 12.31.9
Meridion circulare (Greville) C. Agardh 0.25.7 1.70.261.2 16.2 0.6 41.5 3.0
Nitzschia frustulum (Kützing) Grunow 7.4
Nitzschia inconspicua Grunow 0.7 17.40.90.20.5 2.1
Nitzschia recta Hantzsch ex Rabenhorst 0.2 5.2
Odontidium mesodon (Ehrenberg) Kützing 0.2 0.50.5 18.00.50.5
Planothidium frequentissimum (Lange-Bertalot) Lange-Bertalot 9.73.24.528.61.22.5 1.8 16.01.30.71.30.68.02.3
Planothidium lanceolatum (Brébisson ex Kützing) Lange-Bertalot 12.9 0.7 0.96.013.64.40.4 1.50.01.9 0.50.5
Platessa conspicua (A. Mayer) Lange-Bertalot 11.9
Psammothidium grischunum (Wuthrich) Bukhtiyarova & Round 7.7
Pseudostaurosira alvareziae Cejudo-Figueras, E.A. Morales & Ector 12.3 5.7
Rhoicosphenia abbreviata (C. Agardh) Lange-Bertalot 1.8 0.2 5.90.9
Sellaphora atomoides (Grunow) C.E. Wetzel & Van de Vijver 1.7 1.420.8
Sellaphora nigri (De Notaris) C.E. Wetzel & Ector0.213.2 10.20.7 5.5 11.60.99.00.026.31.0
Sellaphora seminulum (Grunow) D.G. Mann 2.3 1.38.20.60.021.7
Table
Table 4. Values of species richness (R), Shannon-Wiener diversity index (H’) and Pielou’s evenness index (J’) for diatom assemblages of the springs studied.
Table 4. Values of species richness (R), Shannon-Wiener diversity index (H’) and Pielou’s evenness index (J’) for diatom assemblages of the springs studied.
SpringRH’J’
Majorca
FP141.400.37
FL193.190.75
SV161.940.49
Sardinia
FO273.770.79
PB222.930.66
DN182.210.53
DR434.370.81
VA81.440.48
MS544.480.78
SG374.050.78
UL253.220.69
Veneto
BR162.510.76
BA212.900.81
FC102.000.77
LF201.490.45
PR132.470.78
TO161.830.58
Friuli
TR201.940.47
UC261.430.31
BM351.910.39
RR192.800.67
GD342.090.43
AL222.760.62
Table
Table 5. Results of a one-way analysis of similarities (ANOSIM) test evaluating global and pairwise differences of the diatom assemblages in the four geographic regions. Significant values are reported in bold.
Table 5. Results of a one-way analysis of similarities (ANOSIM) test evaluating global and pairwise differences of the diatom assemblages in the four geographic regions. Significant values are reported in bold.
GroupsR Statisticp Significance Level
Global effect 0.6100.001
Pairwise tests:
Veneto vs Sardinia0.2780.008
Veneto vs Majorca0.0680.333
Veneto vs Friuli0.8980.002
Sardinia vs Majorca0.0520.358
Sardinia vs Friuli0.9280.001
Majorca vs Friuli0.9630.012
Table
Table 6. Indicator taxa from the three groups of sites visualized after the CCA according to indicator species analysis (IndVal).
Table 6. Indicator taxa from the three groups of sites visualized after the CCA according to indicator species analysis (IndVal).
Groups of SitesSpeciesCodeIndValp-Value
1 (Friuli)Achnanthidium pyrenaicumADPY0.975 0.001
Achnanthidium lineareACLI0.950 0.002
Gomphonema elegantissimumGELG0.8160.005
Gomphonema tergestinumGTER0.7070.018
2 (Majorca + Sardinia)Cocconeis euglyptaCEUG0.8500.003
3 (Veneto)Amphora inariensisAINA0.9110.001
Sellaphora seminulumSSEM0.8010.007
Sellaphora nigriSNIG0.7910.027

Share and Cite

MDPI and ACS Style

Lai, G.G.; Burato, S.; Padedda, B.M.; Zorza, R.; Pizzul, E.; Delgado, C.; Lugliè, A.; Cantonati, M. Diatom Biodiversity in Karst Springs of Mediterranean Geographic Areas with Contrasting Characteristics: Islands vs Mainland. Water 2019, 11, 2602. https://doi.org/10.3390/w11122602

AMA Style

Lai GG, Burato S, Padedda BM, Zorza R, Pizzul E, Delgado C, Lugliè A, Cantonati M. Diatom Biodiversity in Karst Springs of Mediterranean Geographic Areas with Contrasting Characteristics: Islands vs Mainland. Water. 2019; 11(12):2602. https://doi.org/10.3390/w11122602

Chicago/Turabian Style

Lai, Giuseppina G., Sara Burato, Bachisio M. Padedda, Raffaella Zorza, Elisabetta Pizzul, Cristina Delgado, Antonella Lugliè, and Marco Cantonati. 2019. "Diatom Biodiversity in Karst Springs of Mediterranean Geographic Areas with Contrasting Characteristics: Islands vs Mainland" Water 11, no. 12: 2602. https://doi.org/10.3390/w11122602

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop