# Biodiversity and Environmental Factors Structuring Diatom Assemblages of Mineral Saline Springs in the French Massif Central

Aude Beauger, Olivier Voldoire, Elisabeth Allain, Pierre Gosseaume,
Christelle Blavignac, Lory-Anne Baker, Carlos E Wetzel

► To cite this version:

Aude Beauger, Olivier Voldoire, Elisabeth Allain, Pierre Gosseaume, Christelle Blavignac, et al. Biodiversity and Environmental Factors Structuring Diatom Assemblages of Mineral Saline Springs in the French Massif Central. Diversity, 2023, 15 (2), pp.283. 10.3390/d15020283. hal-03992632

HAL Id: hal-03992632

https://hal.science/hal-03992632

Submitted on 16 Feb 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Biodiversity and Environmental Factors Structuring Diatom Assemblages of Mineral Saline Springs in the French Massif Central

Aude Beauger<sup>1,2,3,* ID</sup>, Olivier Voldoire<sup>1</sup>, Elisabeth Allain<sup>1</sup>, Pierre Gosseaume<sup>1</sup>, Christelle Blavignac<sup>4 ID</sup>, Lory-Anne Baker<sup>3,5</sup> and Carlos E. Wetzel<sup>6</sup>

1 Université Clermont Auvergne, CNRS, GEOLAB, F-63000 Clermont-Ferrand, France

2 LTSER "Zone Atelier Loire", F-37000 Tours, France

3 LTSER "Zone Atelier Territoires Uranifères", F-63000 Clermont-Ferrand, France

4 Centre Imagerie Cellulaire Santé, UCA PARTNER, F-63000 Clermont-Ferrand, France

5 Université Clermont Auvergne, CNRS, LMGE, F-63000 Clermont-Ferrand, France

6 Luxembourg Institute of Science and Technology (LIST), Environmental Research and Innovation Department (ERIN), Observatory for Climate, Environment and Biodiversity (OCEB), 41 rue du Brill, 4422 Belvaux, Luxembourg

* Correspondence: aude.beauger@uca.fr

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Beauger, A.; Voldoire, O.; Allain, E.; Gosseaume, P.; Blavignac, C.; Baker, L.-A.; Wetzel, C.E.

Biodiversity and Environmental Factors Structuring Diatom Assemblages of Mineral Saline Springs in the French Massif Central.

*Diversity* **2023**, *15*, 283.

https://doi.org/10.3390/d15020283

Academic Editors: Michael Wink, Luc Legal and Michel Baguette

Received: 20 January 2023

Revised: 9 February 2023

Accepted: 13 February 2023

Published: 15 February 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Springs are abundant and present worldwide and are among the most threatened ecosystems on Earth. The main aim of this study was to evaluate the biodiversity of seventy-nine (79) mineral saline springs situated in the French Massif Central, focusing on the species richness and the estimated richness and also on the diatom community composition. The influence of physical and chemical factors on both richness and communities was analyzed in order to have a better knowledge of the diatom ecological preferences and identify species typical of saline springs. Since December 2014, an on-going inventory of mineral springs has been in progress. For each spring, physical and chemical characteristics were measured, and benthic diatoms were sampled. The richness was the lowest in the springs presenting a man-made construction around the emergence. In the other springs, the highest richness was associated with the lowest lithium, sodium, total dissolved solid concentrations, and conductivity. Mineralization and some ions (bromine, calcium, chloride, fluoride, lithium, potassium, and sodium) were found to be the most critical drivers of diatom community composition. Some diatom species were typical of specific abiotic conditions, such as *Navicula sanctamargaritae*, which was associated with the highest potassium concentration. These species could appear as bio-indicators of these conditions.

**Keywords:** diatom communities; mineral saline springs; biodiversity; environmental factors

## 1. Introduction

Springs are abundant and present worldwide, "with estimated numbers in the tens of millions" [1]. Spring types have been classified based on hydrology, geology, hydrochemistry, water temperature, biological assemblages, and human use [1-3]. Thus, for example, the springs feeding streams that are crenal (springs situated at headwaters according to the longitudinal models of running waters) [4] can be distinguished from mineral springs that are characterized by high ion concentrations. For this last group, displaying high mineralization, their hydrochemistry reflects the lithology of the aquifer.

Mineral and thermal springs are formed by hydrogeological processes linked to volcanic and tectonic activities. All spring water comes from precipitation that seep into the ground under the combined action of gravity and capillary action to the hydrostatic level [3]. A significant portion of rainwater feeds shallow aquifers, while elsewhere it slowly infiltrates faults up to several kilometers deep. Indeed, mineral springs coincide with the fault zones [5]. Hot waters circulating in the earth's crust will dissolve minerals from the
rocks through which they pass [6]. Then, the water rises through hydrostatic or gaseous pressure with different major and trace elements linked to the deep warming, favoring their dissolution. The knowledge of aquifer temperature and the total electric charge of mobile elements allows us to predict these elemental concentrations in waters [7]. As springs are fed from thermally buffered underground sources, they show a slight variation in temperature [3] and other physical and chemical characteristics, such as the concentrations of sodium, calcium, and potassium ions [5,8].

Springs are ecotones that rely on groundwater to surface water [9], and they are known to be groundwater-dependent ecosystems [10] and isolated biotopes [9,11] capable of hosting specific biocenoses due to their discontinuous presence in the landscape [9]. The associated organisms and communities are often unique and have specific traits and distribution patterns [12–16]. These environments are generated by upwellings, and surface waters often host significant biodiversity and species adapted to the constraints imposed by these ecosystems [14,15].

Many species of algae appear to have developed adaptations to live in water characterized by high salinity and/or high concentrations of specific ions. Most alkaline hot springs contain between 1000 and 2000 mg of total dissolved solids (TDS) per liter; some have salinities higher than seawater [6,17]. Depending on the mineral elements present in the water, such as sodium, some springs create inland saline areas under the influence of saline waters [18]. According to the classification of Krieger [19], springs are moderately salty when TDS varies between 3000 ppm and 10,000 ppm. In this case, different ecosystems may develop where a halophilic flora grows [18,20], including benthic microscopic algae: the diatoms. These algae are characterized by a siliceous skeleton (frustule). They are found in almost every aquatic environment, including fresh and marine waters and terrestrial ecosystem. They are generally considered bioindicators of environmental quality because of their high sensitivity to elements of geogenic origin and the nutrient loading that can occur as a result of contamination [13]. In mineral springs, these microalgae are studied by diatomists all over the world (e.g., [21–31]). Many studies have been performed in Central Europe [21,22,24], others in other countries of Europe, e.g., France [14–16], Italy [25], and Portugal (Azores archipelago) [26]. In other parts of the world, investigations have to date been carried out in Turkey [23], Chile [27], the USA (California) [28], Africa [29], Thailand [30], Russia [31], Iceland [32,33], Canada [32], Japan [32], and New Zealand [33].

In the French Massif Central, numerous mineral saline springs are present and sodium bicarbonate waters form the main hydromineral province in France [5]. These waters have a dominant endogenous character. Their geographical distribution shows a relationship with recent volcanism. This province is made up of basins located on the major volcanogenic directions (N-S, N-W–S-E, N-E–S-W) or at the point of intersection of two of these directions [5]. In this part of France, few studies have been conducted on diatoms living in mineral springs [12,34–36] and on a limited number of mineral springs. Thus, these studies did not address the inventory of the biodiversity of these algae. As a consequence, since the survey initiated in 2014, new taxa have been recently discovered in these springs, such as *Navicula sanctamargaritae* Beauger [14], *Sellaphora labernardierei* Beauger, C. E. Wetzel & Ector [15], *Craticula lecohui* Beauger, C. E. Wetzel & Ector [16], *Pseudostaurosira bardii* Beauger, C. E. Wetzel & Ector [37], and *Chamaepinnularia salina* Beauger, C. E. Wetzel, Allain & Ector [38]. A new genus called *Fontina* Beauger, C. E. Wetzel & Ector in [39] was also observed and described in a mineral spring in that area. Moreover, preliminary research, focusing on diatom biodiversity from 17 springs in this area, was performed between November 2014 and April 2015 [40], and has provided some preliminary information on the ecology of some species encountered in these springs.

The main aim of the present study was to evaluate the biodiversity of the saline springs of the Massif Central related to the species richness and estimated richness of the diatom community. The influence of environmental factors on both diatom richness and communities was also analyzed in order (1) to increase our understanding on the
ecological preferences of the observed diatom species, and (2) to identify species typical of saline springs.

**2. Materials and Methods**

2.1. *Study Site*

Between December 2014 and April 2022, an on-going inventory of mineral springs, and the diatoms inhabiting these habitats, in the French Massif Central was carried out. The area under study covers a zone of approximately 85,000 km² and comprises 15% of the territory of the metropolitan France. It is influenced by oceanic, continental, and Mediterranean climates. Each spring was sampled once only.

Among the mineral springs sampled since 2014 in the French Massif Central, we examined the saline springs for diatom biodiversity with a TDS range varying between 3000 ppm and 10,000 ppm, thereby retaining some 79 emergences according to Krieger [19]. The springs retained are situated in different places of the Massif Central (Figure 1; Table S1). Most are located in the same area (Puy-de-Dôme department). Among the 79 springs, 47 present an artificial construction around the emergence sites, such as Tennis and Lefort springs (Figure 2). For some sites, two or three springs were sampled, such as at Poix springs and Châtel-Guyon springs. Poix springs are peculiar in that they are known to contain bitumen associated with salt water, methane, and traces of hydrogen sulfide [41].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Map of the sampling sites (mineral springs). The correspondence between the spring number and its name is presented in Table S1.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Photography of some studied springs: Tennis (a) and Lefort (b).</figcaption></figure>

## 2.2. Diatom Sampling and Physical and Chemical Analysis

In the field, the sampling sites were georeferenced using a DGPS Geo7x (Trimble, Sunnyvale, CA, USA) in Lambert 93. For each spring, several parameters were measured in-situ: conductivity (µS cm⁻¹), pH (pH unit), and water temperature (°C), using a multi-parameter WTW probe FC Multi 340i, and also dissolved oxygen (% and mg L⁻¹) using a ProODO oxygen probe (Ysi, Yellow Springs, OH, USA). In addition, carbonate concentration (HCO₃⁻) was determined directly in the field using a HACH Digital Titrator, sulfuric acid (0.1600 N and 1.600 N), and Bromocresol Green-Methyl Red Indicator (Hach method 8203) (Hach, Loveland, CO, USA). Moreover, for each spring, a water sample was collected in clean, sterile, 180 mL polypropylene vials with polyethylene caps and was transported in a cooler to the GEOLAB laboratory in Clermont-Ferrand (CNRS UMR 6042, University Clermont-Auvergne), whereupon they were stored at 4 °C.

Diatoms were sampled in relation to the dominant substrate observed for each spring (Table S1). As underlined by Baker et al. [42], the substrate collected in springs was found to influence diatom community composition less than local environmental conditions. The sediments were recovered by scraping off the first few millimeters of mud directly within the vial, and the stones and travertine (a sedimentary rock) by brushing them with a toothbrush. Samples were fixed with 99% ethanol to a final concentration of 70%.

Once back in the GEOLAB laboratory, and following filtration of the samples, the concentrations of the various anions and cations present (as (mg L⁻¹): lithium, sodium, ammonium, potassium, magnesium, calcium, fluoride, chloride, nitrite, nitrate, phosphate, and sulphate) were analyzed within three days using high-pressure ion chromatography. For cations, a Thermo Scientific Dionex ICS1100 system was used, and for anions, a Thermo Scientific Dionex Aquion system (Thermo Fisher Scientific, Courtaboeuf, France) was used.

## 2.3. Slide Preparation and Microscopy

Diatom samples were prepared for light microscopic (LM) observation following the method of Prygiel and Coste [43] by cleaning a small sub-sample of epipelic and epilithic raw material with hydrogen peroxide (H₂O₂, 35% v/v; Caldic France, Cournon d'Auvergne, France) and hydrochloric acid (HCl 37% v/v; Fisher Scientific, Loughborough, UK). The cleaned material was then rinsed several times and subsequently diluted with distilled water to avoid excessive concentrations of diatom valves on the slides. Lastly, a drop of the diluted cleaned material was dried on coverslips and mounted in Naphrax®. LM observations and morphometric measurements were performed using a DM2700M microscope (Leica, Wetzlar, Germany) with a 100× oil immersion objective using differential interference contrast. Various taxonomic keys were used [44–58]. For each spring, ≈400 valves were enumerated on random transects, and the count data were converted to percentage relative abundance of the total count.
## 2.4. Data Analysis

### 2.4.1. Statistical Analyses on the Physical and Chemical Variables

The Piper diagram was created using the data derived from the 79 mineral spring samples [59]. It represents the anion and cation concentrations on two specific triangles, whose sides show the relative contents of each of the major ions compared to the total of these ions (cations on the left triangle, anions on the right triangle). The relative position of an analytical result on each of these two triangles makes it possible to specify the dominance of the variance anionic or cationic species present. For both triangles, a diamond is situated, on which is plotted the intersection of the two lines from the points identified on each triangle. This point of intersection represents the global analysis of the sample; its relative position allows one to specify the type of water concerned, such as bicarbonate-rich with sodium and potassium. For this, we used DIAGRAMMES software (version 6.59) [60].

The physical and chemical variables of the dataset were studied using principal component analyses (PCA). Analysis of $NO_2^-$ ions was not performed as several concentrations were below the detection limit; nor was it performed on dissolved oxygen because of missing values due to a technical problem. Kruskal-Wallis tests were conducted on the physical and chemical variables to test whether the groups discriminated by the PCA were significant.

### 2.4.2. Diatoms Data Statistical Analyses

Because the identification of 300–400 individuals did not provide an exhaustive list of the species richness (S) of a spring [61], we extrapolated the richness values using the Chao1 index [62], which estimates the number of unobserved species based on the number of species observed once or twice in the sampling. Statistical analyses of local diversity were conducted using the species richness, S, and the richness estimated by Chao1, and also the Shannon (H) and Piélou (J) indices. Kruskal-Wallis tests were applied to test for significant differences between taxonomic metrics (S, estimated richness (Chao1) and Shannon-Wiener diversity and Piélou indices) among the groups of springs observed in the PCA.

A Canonical Correspondence Analysis was performed on the species with a relative abundance ≥ 1% of the community. The data used for CCA were log-transformed, after being tested with the Detrended Correspondence Analyses (DCA) [63]. According to DCA analyses, gradient lengths of three or higher values supported the suitability of the data for CCA application [64,65].

The environmental variables with a correlation above 0.5 in the CCA axes were transformed, using the k-means-based discretization method, to obtain classes of concentration used in the Multi-Response Permutation Procedures (MRPP) analysis [66]. MRPP tests were performed to assess differences between lithium, sodium, conductivity, and pH. MRPP is a nonparametric, multivariate procedure that tests whether there is significant difference between two or more groups. The method compares dissimilarities within and between groups. The output of MRPP includes the test statistic T, which indicates separation among groups (a large negative value of T indicates a high relative discrimination of groups), whilst the A statistics estimates the within-group homogeneity. For the A-statistic, if it approaches 1, the groups are completely different; if the A statistic approaches 0, the heterogeneity within groups equals what would be expected by chance; if the A statistic approaches -1, the groups are homogeneous [66]. The *p*-value describes the extent to which the separation between groups is probably due only to chance [66]. For the MRPP, the following options were used: Sorensen Distance Measure, n/sum (n) weighting of groups, and pairwise comparison was made. Lastly, Indicator Species Analysis was used to identify the most characteristic species of each group. The statistical significance of each species indicator value was tested using a Monte-Carlo permutation test (999 permutations, *p* < 0.05).

The analyses were performed using XLSTAT 2016 1.1 (Addinsoft, Paris, France) and PC-ORD 7 softwares (Gleneden Beach, OR, USA) [67,68].
# 3. Results

## 3.1. *Physical and Chemical Characteristics*

Based on the Piper diagram, several springs were bicarbonate-rich with sodium and potassium (blue diamond) (Figure 3). On the opposite side, the orange diamond represents the springs characterized by chlorinated and sulphated calcium waters, such as Font Salado 1 and 2 (numbers 50 and 51 on the map), Pont des Eaux (73), St Coust (46), and Châtel-Guyon 1 and 2 (47, 48). The mineral springs characterized by sodium and potassium chloride or sodium sulphate waters (pink diamond) were, for example, Poix 1 and 3 (52, 53), Croizat (15), and Félix (16). Some springs are intermediate.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Piper diagram highlighting the main water types of the 79 mineral springs sampled. Several springs were bicarbonate-rich with sodium and potassium (blue diamond), others were characterized by chlorinated and sulphated calcium waters (orange diamond). At last, some mineral springs were characterized by sodium and potassium chloride or sodium sulphate waters (pink diamond).</figcaption></figure>

Considering all the springs studied, the conductivity was high and varied between 2112 and 12,370 µS cm⁻¹, while TDS values varied between 1665 and 10,493 mg L⁻¹ (Table 1). The water of the different springs was acidic. Cation concentrations were also high, such as lithium (between 1.04 and 25.73 mg L⁻¹) and sodium (between 185.77 and 2323.05 mg L⁻¹) (Table 1). Chloride concentration varied between 24.88 and 351 mg L⁻¹ (Table 1). For the nitrates, concentrations were low, and the highest value was measured at the spring Poix 3 (24.5 mg L⁻¹). Without this value, the mean was low, at around 0.60 mg L⁻¹ (±1.12).
<table><caption>**Table 1.** Environmental variables measured at the study sites: parameters measured in situ, TDS and cation and anion concentrations (in mg L⁻¹).</caption><thead><tr><th>Environmental Variables</th><th>Mean (±SD)</th></tr></thead><tbody><tr><td>Conductivity (µS cm<sup>-1</sup>)</td><td>6073.85 ± 2112.00</td></tr><tr><td>pH (pH unit)</td><td>6.87 ± 0.47</td></tr><tr><td>Water temperature (°C)</td><td>16.82 ± 8.24</td></tr><tr><td>TDS (mg L<sup>-1</sup>)</td><td>5025.80 ± 1665.47</td></tr><tr><td>Li<sup>+</sup></td><td>6.03 ± 4.08</td></tr><tr><td>Na<sup>+</sup></td><td>1065.66 ± 541.40</td></tr><tr><td>NH<sub>4</sub><sup>+</sup></td><td>1.12 ± 0.78</td></tr><tr><td>K<sup>+</sup></td><td>125.52 ± 63.71</td></tr><tr><td>Mg<sup>2+</sup></td><td>124.63 ± 81.93</td></tr><tr><td>Ca<sup>2+</sup></td><td>209.73 ± 147.31</td></tr><tr><td>F<sup>-</sup></td><td>1.18 ± 1.33</td></tr><tr><td>Cl<sup>-</sup></td><td>996.88 ± 774.70</td></tr><tr><td>Br<sup>-</sup></td><td>2.27 ± 1.88</td></tr><tr><td>NO<sub>3</sub><sup>-</sup></td><td>0.91 ± 2.92</td></tr><tr><td>PO<sub>4</sub><sup>3-</sup></td><td>0.15 ± 0.59</td></tr><tr><td>SO<sub>4</sub><sup>2-</sup></td><td>176.16 ± 243.02</td></tr><tr><td>HCO<sub>3</sub><sup>-</sup></td><td>2312.70 ± 749.37</td></tr></tbody></table>

Figure 4 shows the results of a Principal Component Analysis using the spring physical and chemical variables. PCA axes 1 and 2 explain a total of 38.2 and 16.1%, respectively, of the variance in the environmental data. In the left part of the first factorial plane, a first group (group 1) appears, with the springs highly mineralized with low nitrate and phosphate concentrations. These springs were also carbonated, with a high concentration of lithium, sodium, chloride, and bromine ions present. This group includes springs such as Font de Vie (number 65 on the map), Croizat (15), and the springs situated in the city of St-Nectaire (numbers 26 to 34) and in the area of the former thermal station of Ste Marguerite (numbers 1 to 13). The conductivity of this group was approximately 8000 µS cm⁻¹ on average. Lys, the geyser situated in the city of Vals-les-Bains and Adrienne (14, 60, and 62) was associated with a high concentration of fluoride (5.8 mg L⁻¹ on average) and carbonates (4482 mg L⁻¹ on average). In the right part of the first factorial plane (group 2), there are less mineralized mineral springs (mean conductivity was approximately 4500 µS cm⁻¹). Lastly, Poix 3 (53) was apart and associated with high concentrations of magnesium, calcium, nitrates, and phosphates.

Kruskal-Wallis tests were performed on the physical and chemical variables to determine whether the two groups discriminated using PCA are significant. (Table 2). The two groups were significant for all abiotic variables, magnesium, calcium, fluoride, nitrates, and sulphates excepted.

## 3.2. Diatom Composition and Diversity: Influence of the Environmental Variables

The species richness (S) (mean ± SD = 12.7 ± 10.6) ranged from one species at Tennis spring (number 1) and 67 at Bonnefont spring (number 71), the estimated richness (Chao1) (mean ± SD = 12.7 ± 10.8) from one (at Tennis spring) et 137.8 (at Bonnefont spring) underlining the differences between the studied springs (Table S1). It appeared that for many springs that have a man-made artificial construction around the emergence sites, specific richness and estimated richness (Chao1) was lower than for the natural springs (Table 3). The same result was obtained using the Shannon and Piélou indexes.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Ordination plot (F1 X F2) of a Principal Correspondence Analysis of the samples based on the physical and chemical variables. Green triangles represent the mineral springs studied and the spring number is presented in Table S1. The two groups of springs are presented.</figcaption></figure>

The Shannon–Wiener diversity index ranged from 0 (spring Tennis) to 3.53 (Tête de Lion spring (number 23)), and the Piélou index ranged from 0 (spring Tennis) to 0.269 (Croizat spring (number 15)). To test whether differences existed between the two groups of springs observed in the PCA, Kruskal–Wallis tests were also performed. The species richness (*p* = 0.020) and the Shannon–Wiener diversity index (*p* = 0.043) were statistically significant. Springs in group 2 had a higher biodiversity than those in group 1 (Figure 5).
<table><caption>**Table 2.** Results of the Kruskal–Wallis tests performed on the physical and chemical variables using the two groups created with the PCA. *p*-values in bold are statistically significant.</caption><thead><tr><th>Variable</th><th>p-Value</th></tr></thead><tbody><tr><td>Li<sup>+</sup></td><td><0.0001</td></tr><tr><td>Na<sup>+</sup></td><td><0.0001</td></tr><tr><td>NH<sub>4</sub><sup>+</sup></td><td>0.000</td></tr><tr><td>K<sup>+</sup></td><td><0.0001</td></tr><tr><td>Mg<sup>2+</sup></td><td>0.106</td></tr><tr><td>Ca<sup>2+</sup></td><td>0.896</td></tr><tr><td>F<sup>-</sup></td><td>0.052</td></tr><tr><td>Cl<sup>-</sup></td><td><0.0001</td></tr><tr><td>Br<sup>-</sup></td><td><0.0001</td></tr><tr><td>NO<sub>3</sub><sup>-</sup></td><td>0.848</td></tr><tr><td>PO<sub>4</sub><sup>3-</sup></td><td><0.009</td></tr><tr><td>SO<sub>4</sub><sup>2-</sup></td><td>0.229</td></tr><tr><td>HCO<sup>3-</sup></td><td><0.0001</td></tr><tr><td>Conductivity</td><td><0.0001</td></tr><tr><td>pH</td><td><0.0001</td></tr><tr><td>Water temperature</td><td><0.0001</td></tr><tr><td>TDS</td><td><0.0001</td></tr></tbody></table>

<table><caption>**Table 3.** Mean of the different diversity characteristics (mean ± SD) of the natural versus “artificial” springs.</caption><thead><tr><th>Variables</th><th>Natural Springs</th><th>“Artificial” Springs</th></tr></thead><tbody><tr><td>S</td><td>18.94 ± 12.66</td><td>8.60 ± 6.21</td></tr><tr><td>S.chao1</td><td>26.90 ± 24.02</td><td>10.73 ± 8.54</td></tr><tr><td>Shannon index</td><td>1.57 ± 0.76</td><td>0.91 ± 0.51</td></tr><tr><td>Pielou index</td><td>0.09 ± 0.03</td><td>0.12 ± 0.06</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Box plots of the significant parameters (using Kruskal–Wallis tests) performed on the studied springs. The red crosses represent the mean of data and the black bullets, the outlier points. * represent the statistically significant *p*-value.</figcaption></figure>

Overall, 207 different species were observed (Table S2). The species, with a relative abundance ≥ 1% of the community, encountered in the studied spring included: *Crenotia angustior* (Grunow) Wojtal [24] (Figure 6a), *C. thermalis* (Rabenhorst) Wojtal [24] (Figure 6b–d), *Chamaepinnularia salina* Beauger, C. E. Wetzel, Allain & Ector [38] (Figure 7g), *Fragilaria famelica* (Kützing) Lange-Bertalot [69] (Figure 8a,b), *Halamphora coffeaeformis* (Agardh) Levkov [57] (Figure 8c,d), *Navicula sanctamargaritae* Beauger [14] (Figure 8e), *Navicula*
*veneta* Kützing [70] (Figure 8f), *Nitzschia communis* Rabenhorst [71] (Figure 9a,b), *N. valde-costata* Lange-Bertalot et Simonsen [72] (Figure 9c,d), *Pinnularia kuetzingii* Krammer [73] (Figure 7h), *Planothidium frequentissimum* (Lange-Bertalot) Lange-Bertalot [74] (Figure 7a–d), *Sellaphora labernardierei* Beauger, C. E. Wetzel et Ector [15] (Figure 6), and *Surirella patella* Kützing [70] (Figure 9e,f). Among the species, with a relative abundance < 1% of the whole community, included: *Achnanthidium minutissimum* (Kützing) Czarnecki [75] and *Ad-lafia bryophila* (Petersen) Moser Lange-Bertalot & Metzeltin [76], *Fallacia pygmaea* (Kützing) Stickle et D. G. Mann [77], *Gomphonema lippertii* Reichardt & Lange-Bertalot [78], *Halamphora normanii* (Rabenhorst) Levkov [57], *Mastogloia elliptica* (C. A. Agardh) Cleve [79], *Navicula salinarum* Grunow in Cleve & Grunow [80], and *Pseudostaurosira cubonii* D. M. Williams & C. E. Wetzel [81].

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Scanning Electron Microscopy illustrations of some individuals observed in the studied springs. *Crenotia angustior* (a). *Crenotia thermalis* (b–d).</figcaption></figure>

The eigenvalues of the first CCA axis (0.257) and the second CCA axis (0.114) were both significant (*p* = 0.01 with Monte Carlo permutations test with 999 permutations), and the total explained variance in the species data was 18.5% (Figure 10). The diatom–environment correlations (Pearson) for CCA axis-1 (0.875) and CCA axis-2 (0.681) were high, indicating a relatively strong relationship between diatom species and measured environmental variables. *Pinnularia kuetzingii*, *Fragilaria famelica*, and *Navicula veneta* were associated with the highest values of pH, conductivity, and TDS, and also with the highest concentrations in lithium, potassium, chloride, and bromine ion concentrations. *Crenotia angustior*, *Navicula*
*sanctamargaritae*, *Sellaphora labernardierei*, and *Surirella patella* were linked with the highest concentrations of magnesium and calcium ions, and *Crenotia thermalis* and *Nitzschia valdecostata* with the highest fluoride ion concentration.

<figure><img src="image_10.png" /></figure>

**Figure 7.** Scanning Electron Microscopy illustrations of some individuals observed in the studied springs. *Planothidium frequentissimum* (**a–d**). *Sellaphora labernardierei* (**e,f**). *Chamaepinnularia salina* (**g**). *Pinnularia kuetzingii* (**h**).
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Scanning Electron Microscopy illustrations of some individuals observed in the studied springs. *Fragilaria famelica* (**a,b**). *Halamphora coffeaeformis* (**c,d**). *Navicula sanctamargaritae* (**e**). *Navicula veneta* (**f**).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Scanning Electron Microscopy illustrations of some individuals observed in the studied springs. *Nitzschia valdecostata* (**a,b**). *Nitzschia communis* (**c,d**). *Surirella patella* (**e,f**).</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Canonical Correspondence Analysis (CCA) diagram showing the ordination of diatom species (in blue) and sampling site (i.e., springs, green triangles) along the first two axes and their correlation with environmental variables (in red). The spring numbers are presented in Table S1. TDS: Total Dissolved Solids. CRTH: *Crenotia thermalis*, CRNA: *Crenotia angustior*, CHAMSAL: *Chamaepinnularia salina*, FFAM: *Fragilaria famelica*, HACO: *Halamphora coffeaeformis*, NYSG: *Navicula sanctamargaritae*, NVEN: *Navicula veneta*, NCOM: *Nitzschia communis*, NVLC: *Nitzschia valdecostata*, PKUT: *Pinnularia kuetzingii*, PLFR: *Planothidium frequentissimum*, SLAB: *Sellaphora labernardierei* and SPAT: *Surirella patella*.</figcaption></figure>

The variables with a correlation above 0.5 on the axes of the CCA were lithium, sodium, potassium, calcium, chloride, bromine, fluoride, TDS, conductivity, and pH. Then, they were transformed in classes according to concentration (using a k-means-based discretization method) used in the Kruskal–Wallis and the MRPP analyses (Table 4).

Differences in species richness (S) were statistically significant using Kruskal–Wallis tests between the different classes of lithium, sodium, TDS, and conductivity (Table 5). The highest richness was observed in class 1 of these four variables (Figure 11). Differences in estimated species richness were also statistically significant between the different classes of TDS. For the Shannon index, differences were observed only between the classes of lithium and conductivity. Finally, the Piélou index appears to be different for each of the classes considered.
<table><caption>**Table 4.** Range of environmental variables for each class.</caption><thead><tr><th></th><th>Class 1</th><th>Class 2</th><th>Class 3</th></tr></thead><tbody><tr><td>Li<sup>+</sup> (mg L<sup>-1</sup>)</td><td>[1.04–5.80)</td><td>[5.8–11.30)</td><td>≥11.3</td></tr><tr><td>Na<sup>+</sup> (mg L<sup>-1</sup>)</td><td>[185.77–1178.70)</td><td>[1178.7–1833.50)</td><td>≥1833.5</td></tr><tr><td>K<sup>+</sup> (mg L<sup>-1</sup>)</td><td>[10.46–94.07)</td><td>[94.07–165.44)</td><td>≥164.44</td></tr><tr><td>Ca<sup>2+</sup> (mg L<sup>-1</sup>)</td><td>[11.58–207.50)</td><td>[207.5–410.38)</td><td>≥410.38</td></tr><tr><td>Cl<sup>-</sup> (mg L<sup>-1</sup>)</td><td>[24.88–1118.89)</td><td>[1118.89–2066.06)</td><td>≥2066.06</td></tr><tr><td>Br<sup>-</sup> (mg L<sup>-1</sup>)</td><td>[0–1.44)</td><td>[1.44–3.81)</td><td>≥3.81</td></tr><tr><td>F<sup>-</sup> (mg L<sup>-1</sup>)</td><td>[0.092–3.49)</td><td>≥3.49</td><td></td></tr><tr><td>TDS</td><td>[3018–5175)</td><td>[5175–7195)</td><td>≥7195</td></tr><tr><td>Conductivity (µS cm<sup>-1</sup>)</td><td>[3340–6350)</td><td>[6350–9220)</td><td>≥9220</td></tr><tr><td>pH</td><td>[6.12–6.82)</td><td>[6.82–7.56)</td><td>≥7.56</td></tr></tbody></table>

<table><caption>**Table 5.** Results of the Kruskal–Wallis tests performed on the physical and chemical variables using the classes performed on the lithium, sodium, potassium, calcium, chloride, bromine, fluoride, TDS, conductivity, and pH. *p*-values in bold are statistically significant.</caption><thead><tr><th>Variables</th><th>Species Richness<br>(p-Value)</th><th>Estimated Richness<br>(p-Value)</th><th>Shannon Index<br>(p-Value)</th><th>Pielou Index<br>(p-Value)</th></tr></thead><tbody><tr><td>Li<sup>+</sup></td><td><strong>0.048</strong></td><td>0.118</td><td><strong>0.036</strong></td><td><strong>0.003</strong></td></tr><tr><td>Na<sup>+</sup></td><td><strong>0.019</strong></td><td>0.055</td><td>0.055</td><td><strong>0.010</strong></td></tr><tr><td>K<sup>+</sup></td><td>0.181</td><td>0.384</td><td>0.640</td><td>0.063</td></tr><tr><td>Ca<sup>2+</sup></td><td>0.578</td><td>0.525</td><td>0.335</td><td>0.617</td></tr><tr><td>Cl<sup>-</sup></td><td>0.160</td><td>0.265</td><td>0.309</td><td>0.171</td></tr><tr><td>Br<sup>-</sup></td><td>0.103</td><td>0.084</td><td>0.303</td><td>0.626</td></tr><tr><td>F<sup>-</sup></td><td>0.139</td><td>0.089</td><td>0.094</td><td>0.911</td></tr><tr><td>TDS</td><td><strong>0.001</strong></td><td><strong>0.003</strong></td><td>0.054</td><td><strong>0.035</strong></td></tr><tr><td>Conductivity</td><td><strong>0.020</strong></td><td>0.054</td><td><strong>0.045</strong></td><td><strong>0.005</strong></td></tr><tr><td>pH</td><td>0.172</td><td>0.394</td><td>0.061</td><td>0.159</td></tr></tbody></table>

MRPPs were used to test whether the lithium, sodium, potassium, calcium, chloride, bromine, fluoride, TDS, conductivity, and pH classes induced differences in the diatom communities. Diatom composition was significantly different between the three lithium classes (T = -8.3, A = 0.05, and *p* = 0.00000062); sodium classes (T = -6.0, A = 0.04, and *p* = 0.00000807), potassium classes (T = -4.9, A = 0.03, and *p* = 0.00034706), calcium classes (T = -4.5, A = 0.03, and *p* = 0.00061100), chloride (T = -4.5, A = 0.02, and *p* = 0.00062031), bromine (T = -4.7, A = 0.03, and *p* = 0.00053149), TDS classes (T = -4.8, A = 0.03, and *p* = 0.00044247), the three conductivity classes (T = -6.3, A = 0.04, and *p* = 0.00002596) and the three different pH classes (T = -8.9, A = 0.05, and *p* = 0.00000015). No differences were observed between the two classes of fluoride concentration.

Using Indicator Species Analysis, *Crenotia angustior* was associated with the lowest concentrations of lithium and sodium, the lowest conductivity and pH, but the highest concentration of calcium (such as *Sellaphora labernardierei*) (Table 6). *Planothidium frequentis-simum* lived in waters with the lowest lithium, sodium, bromine ion concentrations (such as *Sellaphora labernardierei*) and the lowest TDS, conductivity, and pH and also in the medium concentration of chloride. In contrast, *Halamphora coffeaeformis* was associated with the highest concentrations of lithium ions (such as *Fragilaria famelica* and *Nitzschia commu-nis*), sodium (such as *Nitzschia valdecostata*), potassium (such as *Navicula sanctamargaritae*), chloride, bromine, and the highest TDS, conductivity, and pH (such as *Navicula veneta*). *Surirella patella* lived in waters characterized by intermediate concentrations of bromine, and *Crenotia thermalis* in waters characterized by intermediate pH.
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Box plots of the species richness made using the classes of the significant variables (using Kruskal-Wallis tests) of the studied springs. The red crosses represent the mean of data and the black bullets, the outlier points.</figcaption></figure>

<table><caption>**Table 6.** List of indicator species for concentration classes with an observed indicator value (IV) greater than 20%, along with mean standard deviation and significant *p*-values (* < 0.05). Different diatoms were linked to each group.</caption><thead><tr><th>Indicator Species</th><th>Observed Indicator Value (IV)</th><th>Mean</th><th>S. dev</th><th>p*</th></tr></thead><tbody><tr><td><strong>Lithium (mg L<sup>-1</sup>)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>[1.04–5.8)</td><td></td><td></td><td></td><td></td></tr><tr><td>Crenotia angustior</td><td>32.2</td><td>14.8</td><td>5.63</td><td>0.0150</td></tr><tr><td>Planothidium frequentissimum</td><td>44.2</td><td>29.3</td><td>4.81</td><td>0.0067</td></tr><tr><td>≥11.3</td><td></td><td></td><td></td><td></td></tr><tr><td>Fragilaria famelica</td><td>27.3</td><td>17.9</td><td>5.51</td><td>0.0701</td></tr><tr><td>Halamphora coffeaeformis</td><td>61.6</td><td>11.6</td><td>5.04</td><td>0.0017</td></tr><tr><td>Nitzschia communis</td><td>28.7</td><td>14.0</td><td>5.36</td><td>0.0234</td></tr><tr><td><strong>Sodium (mg L<sup>-1</sup>)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>[185.77–1178.7)</td><td></td><td></td><td></td><td></td></tr><tr><td>Crenotia angustior</td><td>35.3</td><td>14.8</td><td>5.55</td><td>0.0117</td></tr></tbody></table>
<table><caption>Table 6. *Cont.*</caption><thead><tr><th>Indicator Species</th><th>Observed Indicator Value (IV)</th><th>Mean</th><th>S. dev</th><th>p*</th></tr></thead><tbody><tr><td>Planothidium frequentissimum<br>≥1833.5</td><td>49.4</td><td>29.6</td><td>5.11</td><td>0.0050</td></tr><tr><td>Halamphora coffeaeformis</td><td>47.6</td><td>12.4</td><td>5.61</td><td>0.0017</td></tr><tr><td>Nitzschia valdecostata</td><td>25.4</td><td>13.8</td><td>5.49</td><td>0.0484</td></tr><tr><td><strong>Potassium (mg L<sup>-1</sup>)</strong><br>≥164.44</td><td></td><td></td><td></td><td></td></tr><tr><td>Halamphora coffeaeformis</td><td>36.4</td><td>10.6</td><td>3.81</td><td>0.0010</td></tr><tr><td>Navicula sanctamargaritae</td><td>41.5</td><td>31.7</td><td>3.07</td><td>0.0080</td></tr><tr><td><strong>Calcium (mg L<sup>-1</sup>)</strong><br>≥410.38</td><td></td><td></td><td></td><td></td></tr><tr><td>Crenotia angustior</td><td>45.3</td><td>15.3</td><td>6.01</td><td>0.0030</td></tr><tr><td>Sellaphora labernardierei</td><td>33.1</td><td>15.0</td><td>5.63</td><td>0.0140</td></tr><tr><td><strong>Chloride (mg L<sup>-1</sup>)</strong><br>[1118.89–2066.06]</td><td></td><td></td><td></td><td></td></tr><tr><td>Planothidium frequentissimum<br>≥2066.06</td><td>44.1</td><td>29.6</td><td>4.79</td><td>0.0140</td></tr><tr><td>Halamphora coffeaeformis</td><td>42.4</td><td>12.1</td><td>5.23</td><td>0.0020</td></tr><tr><td><strong>Bromine (mg L<sup>-1</sup>)</strong><br>[0–1.44]</td><td></td><td></td><td></td><td></td></tr><tr><td>Planothidium frequentissimum</td><td>44.0</td><td>28.0</td><td>3.77</td><td>0.0010</td></tr><tr><td>Sellaphora labernardierei<br>[1.44–3.81]</td><td>30.5</td><td>13.7</td><td>4.13</td><td>0.0030</td></tr><tr><td>Surirella patella<br>≥3.81 mg L<sup>-1</sup></td><td>45.8</td><td>21.1</td><td>4.28</td><td>0.0010</td></tr><tr><td>Halamphora coffeaeformis</td><td>33.5</td><td>10.5</td><td>3.85</td><td>0.0020</td></tr><tr><td><strong>TDS (mg L<sup>-1</sup>)</strong><br>[3018–5175]</td><td></td><td></td><td></td><td></td></tr><tr><td>Planothidium frequentissimum<br>≥7195</td><td>55.9</td><td>28.9</td><td>4.62</td><td>0.0010</td></tr><tr><td>Halamphora coffeaeformis</td><td>33.1</td><td>11.6</td><td>4.99</td><td>0.0050</td></tr><tr><td><strong>Conductivity (µS cm<sup>-1</sup>)</strong><br>[3340–6350]</td><td></td><td></td><td></td><td></td></tr><tr><td>Crenotia angustior</td><td>28.8</td><td>15.5</td><td>6.36</td><td>0.0451</td></tr><tr><td>Planothidium frequentissimum<br>≥9220 µS cm<sup>-1</sup></td><td>55.1</td><td>30.0</td><td>5.47</td><td>0.0017</td></tr><tr><td>Halamphora coffeaeformis</td><td>59.5</td><td>12.6</td><td>6.46</td><td>0.0017</td></tr><tr><td><strong>pH</strong><br>[6.12–6.82]</td><td></td><td></td><td></td><td></td></tr><tr><td>Crenotia angustior</td><td>41.2</td><td>15.4</td><td>6.58</td><td>0.0117</td></tr><tr><td>Planothidium frequentissimum<br>[6.82–7.56]</td><td>44.0</td><td>30.1</td><td>5.70</td><td>0.0351</td></tr><tr><td>Crenotia thermalis<br>≥7.56</td><td>44.8</td><td>33.4</td><td>4.03</td><td>0.0017</td></tr><tr><td>Halamphora coffeaeformis</td><td>31.7</td><td>12.8</td><td>6.37</td><td>0.0250</td></tr><tr><td>Navicula veneta</td><td>56.1</td><td>25.5</td><td>6.63</td><td>0.0033</td></tr></tbody></table>

# 4. Discussion

## 4.1. *Physical and Chemical Characteristics of the Saline Springs*

Geothermal activity stimulates water penetration into the rocks, even though the Auvergne volcanoes have been extinct for thousands of years [82]. These superficial infiltrations and the characteristic volcanic imprint reflect and indeed are the cause of the high-water temperature, high mineralization, and water enriched with sodium and chloride observed.
Considering the 79 studied saline springs with a TDS between 3000 ppm and 10,000 ppm, numerous springs were bicarbonate-rich, with sodium and potassium and sodium and potassium chloride or sodium sulphate waters, confirming the saline character of the emergences. Some springs were found to be intermediate in these respects. Although all springs were characterized by high TDS, two groups of springs were observed with highly mineralized emergences (with a high concentration of carbonate, lithium, sodium, chloride, and bromine) with low nitrate and phosphate concentrations in contrast to the others, which were less mineralized. As diatom species are primary producers, they are more directly affected by physical and chemical factors. Thus, these differences in water characteristics could also induce differences in diatom composition. In addition, the Poix 3 spring was seen to differ from the other springs examined in the study and was associated with high concentrations of magnesium, calcium, nitrates, and phosphates due to the geological particularity of the site [83].

## 4.2. Diatom Structure and Composition

These abiotic characteristics are at the origin of a halophilic flora rarely encountered in the continental environment. Thus, organisms such as diatoms associated with springs are often unique [13], and these microalgae have been observed in all studied springs, including those characterized by high water temperature and high mineralization, which are factors generally considered to limit diatom growth (e.g., [84]). The species present are considered cosmopolitan elements, underlining easy dispersal mechanisms [85]. In contrast, some taxa—despite being widespread—are confined to specialized environments such as hot springs [86]. According to Owen et al. [33], these “subcosmopolitan taxa occur wherever suitable habitats are present and consequently develop a mosaic-like geographic distribution”. “Such phycogeographic variability has been reported for diatoms as well as many other phycological groups” [33,87,88].

A total of 207 different species were observed. Diatom communities from saline springs studied in the French Massif Central share many species with other communities from thermal springs in Europe, including Poland [24] and Italy [89,90], but also with other parts of the world, such as Chile [27] and Canada [32]. Because diatoms are known for their high dispersal ability [91], it is perhaps not that surprising for the species with a relative abundance ≥ 1% of the community to be found in different specific regions across Europe, or even the world, with the exception of *Sellaphora labernardierei*, which so far seems to be an endemic species of the Massif Central.

Some rare species were observed, such as *Pseudostaurosira cubonii*, described from a spring in Italia (Sicilia) [81]. It is the first occurrence of this species in France. Other species were also present in freshwater, such as *Achnanthidium minutissimum* and *Adlafia bryophila* [55]. New species were also observed, such as *Navicula* sp. and *Pseudostaurosira* sp., at the spring called “Saurier 3 grotte” and for which a taxonomic description is still in progress. Moreover, the richness was the lowest in the springs that present a man-made construction around the emergence sites, such as Croizat, Lefort, and Tennis springs. For the others, the richness varied between 10 and 137 (estimated richness). These results support those of other studies, such as Cantonati et al. [9] and Lai et al. [90], both of whom observed high species-richness values and suggest the influence of the ecotonal nature of springs and of different local abiotic factors. Differences in species richness were observed between the two groups of springs generated using PCA. The highest richness was associated with the lowest concentrations. This observation was confirmed by the differences observed between the classes of lithium, sodium, TDS, and conductivity, as the highest richness was associated with the lowest concentrations for each of them. Lai et al. [90] observed higher species richness in medium-mineral springs (conductivity 590–1193 µS cm⁻¹), whilst in this study it is associated with a higher conductivity range (3340–6349 µS cm⁻¹).

Among the species with a relative abundance ≥ 1% of the community, *Crenotia thermalis* was present in the mineral springs, except in 13 of them, such as Tennis, Giraudon, and Poix 2 and 3. This species is occasionally encountered in running waters and is charac-
teristic of electrolyte-rich inland habitats, particularly thermal and mineral springs [24,55]. This was also observed, for example, in Italy [89], Poland [24], and Chile [27]. In the case of *C. angustior*, the species was observed in 18 studied springs. According to Wojtal [24], *C. angustior* was generally abundant in waters with very high conductivity, chlorides, sulphates, and other ion concentrations (such as Mg²⁺, Ca²⁺, Mg²⁺, Na⁺, K⁺), but low in nitrate and dissolved oxygen. Only Tennis, Croizat, Giraudon, Poix 2 and 3, and Camuse did not host one of the two *Crenotia* species.

*Chamaepinnularia salina* was present in seven of the springs studied, such as Sail (where it was described), Bard 2, Adrienne, and Vallée Couze d'Ardes. This species is mostly found in mineral springs of Na-K-HCO₃ water type [38].

*Fragilaria famelica* was encountered in 24 studied springs. This species, also encountered in Poland and Portugal, is typical of brackish waters [92–94] and, according to Wojtal [24], is associated with saline springs with a high concentration of chloride and sulphate, but very low dissolved oxygen and nitrate content. This was confirmed by the results of the CCA. Indeed, this species was associated with the highest concentration in chloride but also bromine and potassium.

*Halamphora coffeaeformis* was observed in 12 studied springs. It was mainly the springs in the area of Ste Marguerite (springs 1 to 13 on the map) and of the city of St-Nectaire (springs 26 to 34 on the map). Known to live in water naturally rich in salt and electrolytes [55], the species is here described for the first time in La Montagne spring (Châteldon) and is therefore a denizen of mineral springs [58]. *H. coffeaeformis* has also been sampled in Canada [33], Chile [27], and Tajikistan [95].

*Navicula sanctamargaritae* was present in almost all springs (14 springs excepted). This species is associated with high mineralization levels of water [14,89], as confirmed here by CCA. This species had previously only been encountered in Italy [89]. *Navicula veneta*, a more cosmopolitan species common in electrolyte-rich to brackish waters [51] was also well represented (in 39 studied springs). This species was hitherto observed in various countries, including Chile [27], Tajikistan [95], Poland [24], and Italy [89].

Regarding the two *Nitzschia* species that were often dominant, *N. communis* was present in 15 spring. It has been considered a significant element of the diatom assemblages of thermal springs in volcanic environments in São Miguel Island, Azores [94]. In addition, this species is found in calcium-rich waters of moderate-very high specific conductivity, very low dissolved oxygen, and very low nitrate concentration in Polish springs [24]. This species was associated with high fluoride concentration in the CCA. It has previously been encountered in springs from many countries around the world, including Kenya [33], Iceland [32], Poland [24], Greece [96], and Slovakia [97]. *N. valdecostata* is a saline species found in Chile [27] and also in a thermal spring in California (USA) [28].

*Pinnularia kuetzingii* was present in 27 studied springs. It is known to occur in waters with medium to higher electrolyte content, particularly in thermal springs [48], and has also been observed in Bosnia and Herzegovina and in petrifying *Cratoneurion* springs (petrifying springs with tufa formation) in Belgium [98]. In the CCA, this species was associated with high conductivity and TDS, and also high lithium, sodium, potassium, chloride, and bromine concentrations.

*Planothidium frequentissimum* was also present in most of the springs studied (with the exception of 27, such as Laplace and Ceyssat). This species has earlier been observed in other springs of Europe, such as Italy, except in those with higher temperature and conductivity [90]. This species is widespread and encountered in many springs from different countries around the world, e.g., Belgium [98], Chile [27], Italy [89], Poland [24], Spain [26], and Iran [99].

*Sellaphora labernardierei* was present in 18 mineral springs, such as Sail, Zagat, Etoile, Daguillon, and Saute Bedel, and it seems to be associated with the presence of nitrates [15,76]. In the CCA, it was associated with high calcium concentration.
Lastly, *Surirella patella*, observed in 34 springs studied, is also known to be a brackish water species [88] and was therefore associated with high calcium concentration in the CCA. Previously, it has been encountered in mineral springs in Georgia [100].

### 4.3. *Environmental Variables Associated with Diatom Assemblages*

Mineralization was an important factor in the separation of the saline mineral springs with high TDS and conductivity and high bromine, calcium, chloride, fluoride, lithium, potassium, and sodium ion concentrations. Our results confirmed those of previous studies carried out in other countries where diatom assemblage composition is clearly influenced by conductivity [27,29,30,33,101] and also by high concentrations of calcium, chloride, fluoride, lithium, potassium, and sodium ions [27,90]. In these studies, magnesium, carbonate, and sulphates are also significant variables. Using the MRPP and ISA analyses, *C. angustior* was associated with the lowest concentrations in lithium and sodium, the lowest conductivity and pH but the highest concentration of calcium, while *F. famelica* was linked with the highest lithium concentration. For *H. coffeaeformis*, TDS, conductivity, lithium, sodium, potassium, chloride, and bromine ions were the variables positively correlated with this species, showing it to be an indicator of highly mineralized water. On the contrary, *Planothidium frequentissimum* was found to live in waters with the lowest concentrations of lithium, sodium, and bromine ions; medium chloride ion concentration; and the lowest TDS and conductivity ([3340–6350] µS cm⁻¹). This species was seen to be an indicator of less mineralized water. Lai et al. [90] showed this species to be abundant in springs with conductivity ranging from 5270 to 6270 µS cm⁻¹.

*Navicula sanctamargaritae* was only associated with concentration of potasismus ≥ 164.44 mg L⁻¹. Thus, it could be an indicator of the springs with a high concentration of this cation. *Nitzschia communis* was linked with the highest concentrations in lithium, whilst *Nitzschia valdecostata* was linked with the highest concentrations in sodium. According to Owen et al. [33], in New Zealand, *N. communis* is linked with high carbonate alkalinity. Regarding *Sellaphora labernardierei*, it was found by us to be linked to the highest calcium concentrations and the lowest bromine concentrations, whilst *Surirella patella* was associated with waters characterized by intermediate bromine concentration.

pH was also a significant variable, confirming the results obtained in other countries such as New Zealand, Thailand, and Italy [30,33,89]. In the studied saline springs, pH was acidic (mean ± SD: 6.87 ± 0.47). Using MRPP and ISA analyses, *C. angustior* and *P. frequentissimum* were associated with acidic waters, *C. thermalis* to waters with a pH that varied between 6.82 and 7.55, and *H. coffeaeformis* and *N. veneta* with basic pH.

Water temperature was not found to be an important driver of diatom community composition, contrary to other, previous studies [26,30,31,101].

## 5. Conclusions

The study evaluated and synthesized diatom biodiversity from seventy-nine (79) mineral saline springs located at the Massif Central region of France, along with the influence of environmental factors in the ecology of these communities. As springs are among the most threatened ecosystems on Earth [1,102], it is crucial to better understand their abiotic characteristics and the diatom communities that inhabit these extreme environments, both in order to protect as well as, if necessary, restore them. As found, diatom species richness was lowest in the springs presenting an artificial construction around the emergence sites, thereby underlining the need to restore such sites. In the other springs studied, the highest richness was associated with the lowest lithium, sodium, and TDS concentrations and conductivity. Mineralization and specific ions were the most critical driver of diatom community composition. Some species were typical of abiotic conditions, such as *N. sanctamargaritae*, which was associated with the highest potassium ion concentration; hence, this species would appear to be a good bio-indicator of such conditions. It will be interesting in future studies to investigate further the influence of the abiotic characteristics of all mineral springs analyzed, more especially those with TDS values < 3000 ppm.
**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d15020283/s1, Table S1: List of the studied mineral springs and species richness (S), estimated richness (S.chao1), Shannon and Piélou indexes; Table S2: List of the diatom species encountered.

**Author Contributions:** A.B. conceived the research; A.B., E.A., O.V., L.-A.B. and P.G. performed fieldwork; A.B., E.A., O.V., C.B., P.G. and C.E.W. analyzed the data and wrote and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partly funded by the Direction Régionale de l'Environnement, de l'Aménagement et du Logement Auvergne-Rhône-Alpes, by Clermont Auvergne Métropole, by the "Fédération des Recherches en Environnement" of Clermont-Ferrand (UBP/CNRS FR 3467/INRA), by the OBLA (Réseau d'Observation de la Biodiversité de la Loire et de ses Affluents) of the ZAL (LTSER Zone Atelier Loire) and partly by a grant from the Centre National de la Recherche Scientifique (CNRS, NEEDS, project Acronym: DISCOVER) in the contexte of the Zone Atelier Territoires Uranifères. Co-funding was provided in the framework of the DIATOMS project (LIST—Luxembourg Institute of Science and Technology).

**Institutional Review Board Statement:** Not applicable.

**Data Availability Statement:** The datasets used or analyzed during this study are available from the corresponding author on request.

**Acknowledgments:** The authors thank the two anonymous referees for their valuable and constructive comments and thank Hugh D. Loxdale for his English editing.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Stevens, L.E.; Schenk, E.R.; Springer, A.E. Springs ecosystem classification. *Ecol. Appl.* **2021**, *31*, e2218. [CrossRef] [PubMed]

2. Glazier, D.S. *Encyclopedia of Inland Waters*; Likens, G.E., Ed.; Academic Press Elsevier: Oxford, UK, 2009; Volume 1, pp. 734–735. ISBN 978-0-1237-0626-3.

3. Glazier, D.S. *Reference Module in Earth Systems and Environmental Sciences*; Elsevier: Amsterdam, The Netherlands, 2014. [CrossRef]

4. Ward, J.V. The four-dimensional nature of lotic ecosystems. *J. N. Am. Benthol. Soc.* **1989**, *8*, 2–8. [CrossRef]

5. Risler, J.J. Description et classification géologique des sources minérales et thermales du Massif Central. *BRGM Rep.* **1974**, *74*, 24.

6. Jonker, C.Z.; Van Ginkel, C.; Olivier, J. Association between Physical and Geochemical Characteristics of Thermal Springs and Algal Diversity in Limpopo Province, South Africa. *Water SA* **2013**, *39*, 95–104. [CrossRef]

7. Michard, G. Behaviour of major elements and some trace elements (Li, Rb, Cs, St, Fe, Mn, W, F) in deep hot waters from granitic areas. *Chem. Geol.* **1990**, *89*, 117–134. [CrossRef]

8. Van der Kamp, G. The hydrogeology of springs in relation to the biodiversity of spring fauna: A review. *J. Kans. Entomol. Soc.* **1995**, *68*, 4–17.

9. 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**, *31*, 463–480. [CrossRef]

10. Kløve, B.; Ala-Aho, P.; Bertrand, G.; Boukalova, Z.; Ertürk, A.; Goldscheider, N.; Ilmonen, J.; Karakaya, N.; Kupfersberger, H.; Kvoerner, J.; et al. Groundwater dependent ecosystems. Part I: Hydroecological status and trends. *Environ. Sci. Policy* **2011**, *14*, 770–781. [CrossRef]

11. Werum, M. *Die Kieselalgen Gesellschaften in Quellen: Abhängigkeit von Geologie und Anthropogener Beeinfl Ussung in Hessen (Bunderepublik Deutschland)*; Hessische Landesamt für Umwelt und Geologie: Wiesbaden, Germany, 2001; ISBN 3-89026-331-3.

12. Chaouite, J. Contribution à L'étude des Protistes des Eaux Minérales et Thermo-Minérales en Auvergne. Ph.D. Thesis, Université Blaise Pascal, Clermont-Ferrand, France, 1987.

13. Segadelli, S.; Cantonati, M.; Bertoni, E.; Spitale, D.; Angeli, N.; Borsato, A. Can reference spring diatom communities be predicted from simple aquifer and emergence-site characteristics? In *Abstract Book 9th Use of Algae for Monitoring Rivers and Comparable Habitats (UAMR Ich, June 15th–17th 2015) and International Workshop on Benthic Algae Taxonomy (InBAT, June 17th–19th 2015)*; Cantonati, M., Kelly, M.G., Rott, E., Sabater, S., Stevenson, J.R., Whitton, B.A., Schneider, S., Shubert, L.E., Van de Vijver, B., Vis, M.L., et al., Eds.; Museo Delle Scienze—MUSE: Trento, Italy, 2015; p. 37.

14. Beauger, A.; Voldoire, O.; Mertens, A.; Le Cohu, R.; Van de Vijver, B. Two new *Navicula* species (Bacillariophyceae) from Western Europe. *Phytotaxa* **2015**, *230*, 172–182. [CrossRef]

15. Beauger, A.; Wetzel, C.E.; Voldoire, O.; Garreau, A.; Ector, L. *Sellaphora labernardierei* (Sellaphoraceae, Bacillariophyta), a new epilithic species from French spring and four new combinations within the genus Sellaphora. *Phytotaxa* **2016**, *260*, 235–246. [CrossRef]
16. Beauger, A.; Wetzel, C.E.; Voldoire, O.; Garreau, A.; Ector, L. Morphology and ecology of *Craticula lecohui* sp. nov. (Bacillariophyceae) from hydrothermal springs (Puy-de-Dôme, Massif Central, France) and comparison with similar *Craticula* species. *Nova Hedwig. Beih.* **2017**, *146*, 7–22. [CrossRef] [PubMed]

17. Castenholz, R.W. The thermophilic cyanophytes of Iceland and the upper temperature limit. *J. Phycol.* **1969**, *5*, 360–368. [CrossRef] [PubMed]

18. Reimold, R.J.; Queen, W.H. *Ecology of Halophytes*; Academic Press: New York, NY, USA; London, UK, 1974; p. 607. [CrossRef]

19. Krieger, R.A. The chemistry of saline waters. *Groundwater* **1963**, *4*, 7–12. [CrossRef]

20. Poljakoff-Mayber, A.; Gale, J. *Plants in Saline Environments*; Springer: New York, NY, USA, 1975; Volume 15, p. 219. [CrossRef]

21. Werum, M.; Lange-Bertalot, H. Diatoms in springs from Central Europe and elsewhere under the influence of hydrogeology and anthropogenic impacts. *Iconogr. Diatomol.* **2004**, *13*, 3–417.

22. Żelazna-Wieczorek, J. Diatom flora in springs of Łódź Hills (Central Poland). Biodiversity, taxonomy, and temporal changes of epipsammic diatom assemblages in springs affected by human impact. *Diatom Monogr.* **2011**, *13*, 1–419.

23. Solak, C.N.; Wojtal, A.Z. Diatoms in springs and streams of Türkmen Mt. (Sakarya river basin) common in Turkish inland waters. *Pol. Bot. J.* **2012**, *57*, 375–425.

24. Wojtal, A.Z. *Species Composition and Distribution of Diatom Assemblages in Spring Waters from Various Geological Formations in Southern Poland*; Bibliotheca Diatomologica: Stuttgart, German, 2013; Volume 59, pp. 1–436, ISBN 978-3-443-57050-7.

25. 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**, *356*, 145–157. [CrossRef]

26. Delgado, C.; Gonçalves, V.; Blanco, S.; Almeida, S.F. A new diatom (Bacillariophyceae) species from a thermal spring in Azores archipelago (São Miguel island, Atlantic Ocean). *Bot. Sci.* **2020**, *99*, 169–181. [CrossRef]

27. Angel, A.; Vila, I.; Díaz, C.; Molina, X.; Sepúlveda, P. Geothermal diatoms: Seasonal variability in the El Tatio geothermal field (Altiplano, Chile). *Adv. Microbiol.* **2018**, *8*, 211–234. [CrossRef]

28. Ekins, L.; Rushforth, S.R. Diatom flora of Cowboy Hot Spring, Mono County, California. *Great Basin Nat.* **1986**, *46*, 612–624.

29. Mpawenayo, B.; Cocquyt, C.; Nindorera, A. Les diatomées (Bacillariophyta) et autres algues des sources thermales du Burundi (Afrique Centrale) en relation avec les caractéristiques physiques et chimiques des eaux. *Belg. J. Bot.* **2005**, *138*, 152–164.

30. Pumas, C.; Pruetiworanan, S.; Peerapornpisal, Y. Diatom diversity in some hot springs of northern Thailand. *Botanica* **2018**, *24*, 69–86. [CrossRef]

31. Nikulina, T.V.; Kociolek, J.P. Diatoms from Hot Springs from Kuril and Sakhalin Islands (Far East, Russia). In *The Diatom World. Cellular Origin, Life in Extreme Habitats and Astrobiology*; Seckbach, J., Kociolek, P., Eds.; Springer: Dordrecht, The Netherlands, 2011; Volume 19, ISBN 978-94-007-1326-0.

32. Villeneuve, V.; Pienitz, R. Composition diatomifère de quatre sources thermales au Canada, en Islande et au Japon. *Diatom Res.* **1998**, *13*, 149–175. [CrossRef]

33. Owen, R.B.; Renaut, R.W.; Jones, B. Geothermal diatoms: A comparative study of floras in hot spring systems of Iceland, New Zealand, and Kenya. *Hydrobiologia* **2008**, *610*, 175–192. [CrossRef]

34. Héribaud, J. *Les Diatomées d'Auvergne*; Librairie des Sciences Naturelles: Clermont-Ferrand, France; Paris, France, 1893; p. 255.

35. Chaouite, J.; Romagoux, J.C. La flore diatomique des eaux minérales de Vichy (Allier) et de Chaudes-Aigues (Cantal). *Rev. Des Sci. Nat. D'auvergne* **1989**, *55*, 7–23.

36. Tudesque, L. Etude Hydrobiologique de 3 Sites Halophiles: Bard, Saint-Nectaire et Les Saladis. Master's Thesis, (Diplôme Universitaire Supérieur) (mention Hydrobiol.), Université Blaise Pascal, Clermont-Ferrand, France, 1996; p. 71.

37. Beauger, A.; Wetzel, C.E.; Voldoire, O.; Ector, L. *Pseudostaurosira bardii* (Fragilariaceae, Bacillariophyta), a new species from a saline hydrothermal spring of the Massif Central (France). *Bot. Lett.* **2019**, *166*, 3–13. [CrossRef]

38. Beauger, A.; Wetzel, C.E.; Allain, E.; Bertin, C.; Voldoire, O.; Breton, V.; Baker, L.A.; Kolovi, S.; Biron, D.; Ector, L. *Chamaepinnularia salina* (Bacillariophyta), a new diatom species from French mineral springs (Massif Central). *Phytotaxa* **2022**, *538*, 55–73. [CrossRef]

39. Beauger, A.; Wetzel, C.E.; Voldoire, O.; Allain, E.; Breton, V.; Miallier, D.; Ector, L. *Fontina* Gen. nov. (Bacillariophyta): A new diatom genus from a thermo-mineral spring of the French Massif Central (France). *Diatom Res.* **2022**, *37*, 51–61. [CrossRef]

40. Beauger, A.; Voldoire, O.; Wetzel, C.E.; Allain, E.; Millan, F.; Breton, V.; Kolovi, S.; Ector, L. Biodiversité et écologie des diatomées dans les sources minérales de la région de Sainte Marguerite (Saint-Maurice-ès-Allier, Massif central, France). *BIOM* **2020**, *1*, 21–34. [CrossRef]

41. Baker, L.A.; Beauger, A.; Wetzel, C.E.; Voldoire, O.; Blavignac, C.; Allain, E.; Ector, L.; Biron, D. Brackish diatom species observed at the bituminous Puy de la Poix spring: An island of curiosity. *BIOM* **2022**, *3*, 40–51. [CrossRef]

42. Baker, L.A.; Biron, D.; Millan, F.; Voldoire, O.; Breton, V.; Allain, E.; Wetzel, C.E.; Ector, L.; Beauger, A. The substrate, a key factor or not, to explain the species diversity of diatom communities in mineral springs. *Bot. Lett.* **2022**, *169*, 155–165. [CrossRef]

43. Prygiel, J.; Coste, M. *Guide Méthodologique Pour la Mise en Œuvre de L'indice Biologique Diatomées NF T 90-354*; Etude Agences de l'Eau-Cemagref Bordeaux, Agences de l'Eau: Douai, France, 2000; p. 134.

44. Krammer, K.; Lange-Bertalot, H. *Bacillariophyceae*. 1. Teil: Naviculaceae. In *Süsswasserflora von Mitteleuropa*, 2nd ed.; Ettl, H., Gerloff, J., Heynig, H., Mollenhauer, D., Eds.; Gustav Fischer Verlag: Stuttgart, Germany, 1997; Volume 2, p. 876.

45. Krammer, K.; Lange-Bertalot, H. *Bacillariophyceae*. 2. Teil: Bacillariaceae, Epithemiaceae, Surirellaceae. In *Süsswasserflora von Mitteleuropa*, 2nd ed.; Ettl, H., Gerloff, J., Heynig, H., Mollenhauer, D., Eds.; Gustav Fischer Verlag: Stuttgart, Germany, 1997; Volume 2, p. 611.
46. Krammer, K.; Lange-Bertalot, H. *Bacillariophyceae*. 3. Teil: *Centrales, Fragilariaceae, Eunotiae*. In *Süsswasserflora von Mitteleuropa*, 2nd ed.; Ettl, H., Gerloff, J., Heynig, H., Mollenhauer, D., Eds.; Gustav Fischer Verlag: Stuttgart, Germany, 2000; Volume 2, p. 576.

47. Krammer, K.; Lange-Bertalot, H. *Bacillariophyceae*. 5. Teil: English and French translation of the keys. In *Süsswasserflora von Mitteleuropa*, 2nd ed.; Büdel, B., Gärtner, G., Krienitz, L., Lokhorst, G.M., Eds.; Spektrum Akademischer Verlag: Heidelberg, Germany, 2000; Volume 2, p. 599.

48. Krammer, K.; Lange-Bertalot, H. *Bacillariophyceae*. 4. Teil: *Achnanthaceae*, kritische Ergänzungen zu *Navicula* (Lineoatae) und *Gomphonema*. In *Süsswasserflora von Mitteleuropa*, 2nd ed.; Ettl, H., Gärtner, G., Heynig, H., Mollenhauer, D., Eds.; Spektrum Akademischer Verlag: Heidelberg, Germany, 2004; Volume 2, p. 468.

49. Krammer, K. The genus *Pinnularia*. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2000; Volume 1, p. 703.

50. Krammer, K. Cymbella. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2002; Volume 3, p. 584.

51. Krammer, K. Cymbopleura, Delicata, Navicymbula, Gomphocymbellopsis, Afrocymbella. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2003; Volume 4, p. 530.

52. Lange-Bertalot, H.; Navicula sensu stricto, 10 Genera separated from Navicula sensu lato Frustulia. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2001; Volume 2, p. 526.

53. Kulikovskiy, M.S.; Lange-Bertalot, H.; Witkowski, A.; Dorofeyuk, N.I.; Genkal, S.I. *Diatom Assemblages from Sphagnum Bogs of the World. Part I: Nur bog in Northern Mongolia*; Bibliotheca Diatomologica: Stuttgart, German, 2010; Volume 55, pp. 1–326, ISBN 978-3-443-57046-0.

54. Lange-Bertalot, H.; Bak, M.; Witkowski, A.; Tagliaventi, N. Eunotia and some related genera. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2011; Volume 6, p. 747.

55. 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; Cantonati, M., Kelly, M.G., Lange-Bertalot, H., Eds.; Koeltz Botanical Books: Schmitten-Oberreifenberg, Germany, 2017; p. 842.

56. 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**, *15*, 203–234. [CrossRef]

57. Levkov, Z.; Mitic-Kopanja, D.; Reichardt, E. The diatom genus *Gomphonema* from the Republic of Macedonia. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2016; Volume 8, p. 552.

58. Levkov, Z. Amphora sensu lato. In *Diatoms of Europe*; Lange-Bertalot, H., Ed.; A.R.G. Gantner Verlag Kommanditgesellschaft: Ruggell, Liechtenstein, 2009; Volume 5, p. 916.

59. Piper, A.M. A Graphical Procedure in the Geochemical Interpretation of Water Analysis. *Eos Trans. Am. Geophys. Union* **1944**, *25*, 914–928. [CrossRef]

60. Simler, R. Software "Diagrammes"; Laboratoire d'Hydrologie d'Avignon, Université d'Avignon et Pays du Vaucluse: Avignon, France, 2018; Available online: http://www.lha.univ-avignon.fr (accessed on 18 January 2023).

61. Soininen, J.; Heino, J.; Kokocinski, M.; Muotka, T. Local-regional diversity relationship varies with spatial scale in lotic diatoms. *J. Biogeogr.* **2009**, *36*, 720–727. [CrossRef]

62. Chao, A. Nonparametric Estimation of the Number of the classes in a Population. *Scand. J. Stat.* **1984**, *11*, 265–270. [CrossRef]

63. Hill, M.O.; Gauch, H.G.J. Detrended correspondence analysis: An improved ordination technique. *Vegetatio* **1980**, *42*, 47–58. [CrossRef]

64. Ter Braak, C.J.F.; Prentice, I.C. A theory of gradient analysis. *Adv. Ecol. Res.* **1988**, *18*, 271–317. [CrossRef]

65. Birks, H.J.B. Quantitative palaeoenvironmental reconstructions. In *Statistical Modelling of Quaternary Science Data*; Maddy, D., Brew, J.J., Eds.; Quaternary Research Association: Cambridge, UK, 1995; pp. 161–254.

66. McCune, B.; Grace, J.B. *Analysis of Ecological Communities*; MJM Software Design: Gleneden Beach, OR, USA, 2002; p. 304.

67. Addinsoft. XLSTAT Statistical and Data Analysis Solution. Paris, France, 2022. Available online: https://www.xlstat.com (accessed on 18 January 2023).

68. McCune, B.; Mefford, M.J. PC-ORD. In *Multivariate Analysis of Ecological Data*. Version 7; MJM Software Design: Gleneden Beach, OR, USA, 2016.

69. Lange-Bertalot, H. Zur systematischen Bewertung der bandförmigen Kolonien bei *Navicula* und *Fragilaria*. Kriterien für die Vereinigung von *Synedra* (subgen. *Synedra*) Ehrenberg mit *Fragilaria* Lyngbye. In *Nova Hedwigia*; 1980; Volume 33, pp. 723–787.

70. Kützing, F.T. *Die Kieselschaligen Bacillarien Oder Diatomeen*; Zu finden bei W. Köhne: Nordhausen, Germany, 1844; p. 152.

71. Rabenhorst, L. *Die Algen Sachsens. Resp. Mittel-Europa's Gesammelt und Herausgegeben von Dr. L. Rabenhorst, Dresden, Germay, 1848-1860, Dec. 1-100. No. 1-1000*; Harvard University: Cambridge, MA, USA, 1848.

72. Lange-Bertalot, H.; Simonsen, R. A taxonomic revision of the *Nitzschiae lanceolatae* Grunow. 2. European and related extra European freshwater and brackish water taxa. *Bacillaria* **1978**, *1*, 11–111.

73. Krammer, K. *Pinnularia: Eine Monographie der Europäischen Taxa*; Bibliotheca Diatomologica: Berlin, Germany; Stuttgart, Germany, 1992; Volume 26, pp. 1–353.

74. Lange-Bertalot, H. Neue Kombinationen von Taxa aus *Achnanthes* Bory (sensu lato). *Iconogr. Diatomol.* **1999**, *6*, 270–283.
75. Czarnecki, D.B. The freshwater diatoms culture collection at Loras College, Dubuque, Iowa. In Proceedings of the 11th International Diatom Symposium, San Francisco, CA, USA, 12–17 August 1990; Memoirs of the California Academy of Sciences, 1994, Kockiolek, J.P., Eds.; Volume 17, pp. 155–174.

76. Moser, G.; Lange-Bertalot, H.; Metzeltin, D. Insel der Endemiten. Geobotanisches Phänomen Neukaledonien. Island of Endemics New Caledonia—A Geobotanical Phenomenon. *Bibliotheca Diatomologica: Berlin, Germany; Stuttgart, Germany*, 1998; Volume 38, pp. 1–464.

77. Round, F.E.; Crawford, R.M.; Mann, D.G. *The Diatoms Biology and Morphology of the Genera*; Cambridge University Press: Cambridge, UK, 1990; p. 747.

78. Reichardt, E. Zur Revision der Gattung *Gomphonema*. Die Arten um G. *affine/insigne*, G. *angustatum/micropus*, G. *acuminatum* sowie gomphonemoide Diatomeen aus dem Oberligozän in Böhmen. In *Iconographia Diatomologica*; A.R.G. Gantner: Ruggell, Liechtenstein, 1999; Volume 8, pp. 1–203.

79. Schmidt, A. *Atlas der Diatomaceen-Kunde*; Series IV; Heft 47; O.R. Reisland: Leipzig, Germany, 1893; pp. 185–188.

80. Cleve, P.T.; Grunow, A. Beiträge zur Kenntniss der arctischen Diatomeen. *Kongliga Sven. Vetensk.-Akad. Handl.* **1880**, *17*, 1–121.

81. Williams, D.M.; Wetzel, C.E. Description of a new *Pseudostaurosira* based on "*Fragilaria virescens* f. *parva*" from Erbario Crittogamico Italiano. *Bot. Lett.* **2020**, *167*, 86–94. [CrossRef]

82. Boineau, R.; Maisonneuve, J. Les sources minérales du Massif central français et leur cadre géologique. *Rapp. BRGM 72-SGN-151-MCE* **1972**, 32. Available online: http://infoterre.brgm.fr/rapports/72-SGN-151-MCE.pdf (accessed on 18 January 2023).

83. Gonnard, F. Sur le quartz du calcaire bitumineux de la Limagne. *Bull. Société Française Minéralogie* **1906**, 29, 362–365. [CrossRef]

84. Żelazna-Wieczorek, J.; Olszyński, R.M.; Nowicka-Krawczyk, P. Half a century of research on diatoms in athalassic habitats in central Poland. *Oceanol. Hydrobiol. Stud.* **2015**, *44*, 51–67. [CrossRef]

85. Hoffman, L. Geographic distribution of blue-green algae. *Hydrobiologia* **1996**, *336*, 33–40. [CrossRef]

86. Padisák, J. Phytoplankton. In *The Lakes Handbook: Limnology and Limnetic Ecology*; O'Sullivan, P.E., Reynolds, C.S., Eds.; Blackwell Science Ltd.: Singapore, 2003; Volume 1, pp. 251–308, ISBN 978-0-6320-4797-0.

87. Coesel, P.F.M. Biogeography of desmids. *Hydrobiologia* **1996**, *336*, 41–53. [CrossRef]

88. Tyler, P.A. Endemism in freshwater algae: With special reference to the Australian region. *Hydrobiologia* **1996**, *336*, 127–135. [CrossRef]

89. 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. [CrossRef]

90. Lai, G.G.; Ector, L.; Wetzel, C.E.; Lugliè, A.; Padedda, B.M. Environmental factors structuring diatom assemblages in thermo-mineral springs of Sardinia, Italy. *Freshw. Sci.* **2022**, *41*, 45–61. [CrossRef]

91. Kristiansen, J. Dispersal of freshwater algae—A review. *Hydrobiologia* **1996**, *336*, 151–157. [CrossRef]

92. Van Dam, H.; Mertens, A.; Sinkeldam, J. A coded checklist and ecological indicator values of freshwater diatoms from The Netherlands. *Neth. J. Aquat. Ecol.* **1994**, *28*, 117–133. [CrossRef]

93. Rakowska, B. Diatom communities in a salt spring at Pełczyska (Central Poland). *Biologia* **1997**, *52*, 489–493.

94. Quintela, A.; Almeida, S.F.P.; Terroso, D.; Ferreira da Silva, E.; Forjaz, V.; Rocha, F. Diatom assemblages of thermal and mineral waters from volcanic environments in São Miguel Island, Azores. *Diatom Res.* **2013**, *28*, 407–417. [CrossRef]

95. Niyatbekov, T.; Barinova, S. Diatom Species Richness in Algal Flora of Pamir, Tajikistan. *Eur. Sci. J.* **J. 2018**, *14*, 301–323. [CrossRef]

96. Economou-Amili, A. *On Diatoms from Thermal Springs of Greece*; Institute of Systematic Botany, University of Athens: Athens, Greece, 1976.

97. Rehakova, Z. Diatoms from thermal waters and mud in Piestany Spa (Slovakia). *Arch. Hydrobiol.* **1976**, *14*, 1–175.

98. Denys, L.; Oosterlynck, P. Diatom assemblages of non-living substrates in petrifying *Cratoneurion* springs from lower Belgium. *Fottea* **2015**, *15*, 123–138. [CrossRef]

99. Musaabad, L.A.; Mirzahasanlou, J.P.; Mahmoodlu, M.G.; Bahlakeh, A. Diatom flora in three Springs of Golestan Province. *J. Phycol.* **2019**, *3*, 432–442. [CrossRef]

100. Barinova, S.; Kukhaleishvili, L.; Nevo, E.; Janelidze, Z. Diversity and ecology of algae in the Algeti National Park as a part of the Georgian system of protected areas. *Turk. J. Bot.* **2011**, *35*, 729–774. [CrossRef]

101. Leira, M.; Meijide-Failde, R.; Torres, E. Diatom communities in thermo-mineral springs of Galicia (NW Spain). *Diatom Res.* **2017**, *32*, 29–42. [CrossRef]

102. Cantonati, M.; Fensham, R.J.; Stevens, L.E.; Gerecke, R.; Glazier, D.S.; Goldscheider, N.; Knight, R.L.; Richardson, J.S.; Springer, A.E.; Tockner, K. Urgent plea for global protection of springs. *Conserv. Biol.* **2020**, *35*, 378–382. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.