Abstract

Bryological composition, water chemistry, and environmental factors were characterized on 67 Belgian travertines. We explore the relationship between these environmental factors and the community composition, species richness, or presence of individual species using Redundancy Analysis with Hellinger’s transformation (tb-RDA) or Generalized Linear Models (GLMs). The best variables explaining the community composition are slope, NO3, NH4, and PO4. The species richness is negatively related to canopy cover and PO4. Palustriella commutata tends to be more frequent when the slope is steeper and to a lesser degree when the canopy cover is lower. Eucladium verticillatum tends to be slightly more frequent when canopy cover and NH4 concentrations are lower. Cratoneuron filicinum is more frequent at higher Mg concentrations and Pellia endiviifolia is more frequent at lower PO4 concentrations and higher NO3 concentrations. Brachythecium rivulare showed wide ecological amplitude and almost none of the tested environmental factors seem to be related to its presence. The study identifies eutrophication as the main factor responsible for habitat deterioration. Practical indications on the best ways to maintain or to enhance the quality of these petrifying sources are given.

1. Introduction

“Travertine” or “tufa” deposits include a wide variety of calcareous substrates that are characterized not only by their mineral composition but also by their morphology and processes by which they form. There are many interpretations and definitions of the terms tufa and travertine (Symoens et al. [1], Couderc [2], Pentecost [3], De Zuttere [4], Viles and Goudie [5], Ford and Pedley [6], Janssen and Swennen [7], Merz-Preiß and Riding [8], Boch et al. [9], Franco et al. [10], and Brusa and Cerabolini [11]). The term “travertine” is often used to describe calcareous substrates formed during warm periods or in warm environment including thermal springs and that contain no or very few organic material other than bacteria [6, 7]. In contrast, the term “tufa” is generally used to describe substrates formed in colder environments and containing organic material [6]. It is possible for both warm and cold travertine formation to occur simultaneously, as evidenced by several Belgian locations including the Hoyoux river near Huy and on the “Ry de Matignolles” stream near Treignes [10].

Although the European definition of the habitat [12] uses the term “tufa” we will use the term “travertine” in accordance with Pentecost [13] as the term “tufa” relates more to a soft and poorly consolidated variety of travertine.

In the framework of the mapping of Natura 2000 sites in Wallonia, south Belgium, we had the opportunity to undertake research on the ecology of the priority habitat “Petrifying springs with tufa formation (Cratoneurion) 7220.” The code 7220 is specific to this priority habitat in the Interpretation Manual of European Union Habitats [12]. This priority habitat generally consists of small point or linear formations dominated by bryophytes (Cratoneurion commutati Koch, 1928, communities). The manual [12] provides a list of characteristic species for each habitat type; however, these lists have to be adapted at the national or regional scale. In Wallonia, two bryophyte species are considered strictly linked to the 7220 habitat, Palustriella commutata (Hedw.) Ochyra and Eucladium verticillatum (With.) Bruch & Schimp. Palustriella commutata (Hedw.) Ochyra has been split into two species by some authors [14]; however, in this study we consider both as Palustriella commutata.

The following are considered to be good companion species and are also included in the Walloon definition when assessing the habitat conservation status [15]: Aneura pinguis (L.) Dumort., Conocephalum conicum (L.) Dumort., Jungermannia atrovirens Dumort., Leiocolea badensis (Gottsche.) Jörg., Pellia endiviifolia (Dicks.) Dumort., Preissia quadrata (Scop.) Nees, Brachythecium rivulare Schimp., Bryum pseudotriquetrum (Hedw.) P. Gaertn. et al., Ctenidium molluscum (Hedw.) Mitt., Cratoneuron filicinum (Hedw.) Spruce, Dichodontium pellucidum (Hedw.) Schimp., Didymodon tophaceus (Brid.) Lisa, Fissidens adianthoides Hedw., Fissidens crassipes Wilson ex Bruch & Schimp., Gymnostomum calcareum Nees & Hornsch., Hymenostylium recurvirostrum (Hedw.) Dixon, Philonotis calcarea (Bruch & Schimp.) Schimp., Plagiomnium ellipticum (Brid.) T. J. Kop., Plagiomnium rostratum (Schrad.) T. J. Kop., Plagiomnium undulatum (Hedw.) T. J. Kop., and Seligeria trifaria (Brid.) Lindb. Cratoneuron filicinum is considered a good characteristic species in other regions and countries, for example, Flanders and the Brussels-Capital Region [16], Great Britain [17], and the Netherlands [18]. However, it was not included in the small group of characteristic species as it is far more widespread in Wallonia than the two other species [19], occurring in a variety of different habitats, including man-made calcareous forest tracks. Nonetheless Cratoneuron filicinum is considered to be a good companion species for the 7220 habitat in Wallonia [15].

In Wallonia [15], the 7220 habitat is mainly associated with calcareous springs and small streams (width < 5 m) where incrustation processes occur. An atypical form of the habitat is represented in Wallonia by the “travertins” formed in the Hoyoux river consisting of calcareous deposits across the river (barrages) or its tributaries whose width can reach 10 to 20 m. The habitat in Wallonia [15] is characterized by Ca-rich waters (110–120 mg·L−1) and high pH values (7.5–8.5) but depending on the subregion the substrate can be mostly sandy with calcareous incrustations (sandy-loam region) or more compact calcareous rocks (Condroz, Famenne-Calestienne, Lorraine) where they are locally called “crons” or “crânières” [1, 2024].

From a phytosociological point of view, most of the springs and small streams of Wallonia and Brussels-Capital Region can be attributed to the Montio-fontanae-Cardaminetea-amarae Braun-Blanquet et Tüxen, 1943, Class [16, 25, 26]. Following Bardat and Hauguel [26] and depending on the local situations the bryophytes communities can then be linked to the Caricion remotae Kästner, 1941, Alliance (intraforest communities of oligo-mesotrophic soils), to the Pellion endiviifoliae Bardat, 1998, prov.nov. Alliance (neutro-alkaline small streams), or to the Riccardio-Eucladion verticillati Bardat, 1998, prov.nov. Alliance (Ca-rich and thermophilous soils). Following Zechmeister and Mucina [16, 25], Belgian communities can be linked to the Cratoneurion commutati Koch, 1928, Alliance (sunny springs with highly oxygen-saturated waters), to the Adiantion Br.-Bl. ex Horvatic, 1934, em. Hoc. Loco Alliance (waters with high Ca contents and high summer temperatures), or to the Caricion remotae Kästner, 1941, Alliance (intraforest communities of oligotrophic waters).

The bryological aspects of travertines and their phytosociological composition are well documented in Belgium [1, 4, 16, 21, 22, 27, 28], the Netherlands [18, 27], Italy [29], France [2, 26, 30], and in a wider European context [25]. However, relatively few studies have been undertaken on datasets that consider the driving environmental variables that could influence the distribution of particular species and the species richness of these habitats [11, 3135].

The aims of this research were (1) to learn more about the variables, both physical and chemical, that can affect bryophyte diversity in the surveyed sites; (2) to evidence the most significant variables that can affect the presence or absence of the more characteristic species; (3) to give practical indications on the best ways to maintain or to enhance the quality of this habitat both in Wallonia and the Brussels-Capital Region.

2. Materials and Methods

2.1. Data Collection

Most of the surveyed sites lie in Wallonia, south Belgium (Figure 1), where 62 sites were identified, representing the widest variety of conditions within the 7220 habitat. Wallonia is not a homogeneous region as it consists of 5 subregions characterized by geomorphological specificities (from north to south: sandy Loess Plateau, psammitic Condroz, calcareous Famenne-Calestienne, acidic siliceous Ardennes, and calcareous Lorraine). Petrifying sources are mainly small habitats; however, they are relatively widespread in all these subregions except the Ardennes where it is virtually absent due to the acidic siliceous substrate and in the sandy Loess Plateau where there were also fewer sites. To increase the coverage, we added 5 sites within the Loess region by including some of the Brussels-Capital Region localities.

On each Walloon site a comprehensive sampling regime (random selection of all known locations) was undertaken during autumn, 2012; the bryophytes were identified in the laboratory (A. Sotiaux). Sampling of the 5 locations in the Brussels-Capital Region was undertaken during autumn, 2014, with bryophytes identified in the laboratory (J.-M. Couvreur). We collected bryophytes that are directly associated with areas of travertine and to avoid “noise” in the final dataset we also dropped before analysis species that were not directly related to the 7220 habitats. Moreover, as it is sometimes difficult to assess the rate of abundance or cover by these bryophytes especially in this kind of rather heterogeneous habitat, we decided to only mention the presence/absence of these species.

Distance from the spring was not measured but the physical and physicochemical variables, including slope, orientation, canopy cover, and water chemistry, were measured at the location of each bryophyte sample. The physical variables were measured first, including slope, orientation, and canopy cover. Slope was estimated visually and classified into three groups: low (0.5–7.5%), moderate (7.5–20%), and steep (>20%). The orientation was measured using a compass (the four cardinal points plus the four intermediate ones). The percentage of canopy cover (%) was measured visually. Water chemistry, the second group of variables, was measured at each site on the same day that the bryophytes were collected. The water sample was obtained by collecting one litre of water flushing through the site in a glass bottle. The glass bottles were kept cool by storing them in cool boxes with freeze packs, until they were transferred to a laboratory fridge. The pH was measured in the lab using a WTW197i multimeter; other physicochemical variables were measured using a DR3900 Hach Lange photometer: conductivity (μS/cm), nitrates (N-NO3 mg·L−1), ammonium (N-NH4 mg·L−1), soluble orthophosphates (P-PO4 mg·L−1), nitrites (N-NO2 mg·L−1), water hardness (dH degrees), calcium (Ca mg·L−1), and magnesium (Mg mg·L−1). Information including locality, date, location ( and Lambert 1972 coordinates), and biogeographical region (Atlantic or Continental mentioned in the variable “Zone”) was also recorded at each site.

2.2. Analysis

All of the analyses were performed with R 3.2.1 [36]. We removed one dataset from the analysis (Leffe6) for which NO2 and NO3 values were completely unrealistic, even in highly eutrophicated waters. Four missing values in three different explanatory variables (two for PO4, one for Ca, and one for Mg) were modelled because the AIC approaches we intended to use are incompatible with missing values. We built three Gaussian linear models, one for each variable with missing values, with the other environmental variables as predictors; however, the species data were not used. The missing data were replaced by the predicted value from these models.

The relationship between bryophyte community composition (presence/absence matrix) and the environmental factors was studied with a Hellinger transformation based Redundancy Analysis (tb-RDA [37]) performed with vegan 2.3-1 package [38]. We checked that the environmental variables were able to explain at least a part of the variance of the species matrix with a global permutation test (ANOVA function in vegan, 999 permutations). We then applied a simple forward selection procedure based on sequential permutation tests (ordistep function from vegan) to identify the environmental variables that are most related to species composition.

The relationship between species richness (total number of species) and the environmental predictors was analysed with a Poisson Generalized Linear Model (GLM). Binomial GLMs were used to characterize the presence/absence of a few individual species relative to the environmental predictors. These binomial GLMs were built only for species that were present on 20 of the 67 sites (Palustriella commutata, Eucladium verticillatum, Cratoneuron filicinum, Pellia endiviifolia, and Brachythecium rivulare).

For all GLMs, model conditions (linearity, distribution, and outliers) were checked with residuals plots and overdispersion was checked for the Poisson GLM. Multicollinearity between the explanatory variables was checked graphically and by computing Variance Inflation Factors (VIFs). If necessary, some of the explanatory variables were dropped in order to keep the VIFs < 5. We centred the pH explanatory variable on its mean value to avoid predictions at the intercept for pH = 0, that is, completely outside the range of observed values. The absence of spatial correlation of the models residuals was checked with splines correlograms.

In order to determine the most important explanatory variables for these GLMs, we applied an AICc based model selection procedure as described by Burnham and Anderson [39] and shortly summarized hereafter. GLMs with all possible combinations of explanatory variables were computed. Their AICc and AICc model weights were calculated. The AICc model weight quantify the model selection uncertainty and can be interpreted as the probability for a given model to be selected as the best model (the one with lowest AICc) in a set of models if we could resample the data. Then for each explanatory variable we computed a variable weight as the sum of the AICc weights of the models in which this explanatory variable is present. This variable weight is a measure of the relative importance of the explanatory variables. We interpreted only the explanatory variables with an AICc weight > 0.6. The shrinkage model averaged coefficients were then computed along with their unconditional standard errors. These model averaged coefficients are the mean of the coefficient of all models weighted by the model quality (AICc model weights).

For the biological interpretation of the GLMs outputs we plotted graphs, based on model averaged coefficients, for the most important explanatory variables in the range of their observed values. The explanatory variables that are not drawn on the graphs are fixed to their mean value.

3. Results

3.1. Dataset

The results of the field campaign are presented in Table 1(a) (physicochemical variables) and Table 1(b) (species) as two matrices with the 67 sites as lines and the species or environmental factors as columns. The “Leffe6” site is mentioned; however, it was removed due to “out of range” values of NO2 and NO3. The 4 missing values were replaced by their predicted ones in the final matrix by using an explanatory GLM model based on the other variables. These values are, respectively, 0.055 and 0.058 mg·L−1 PO4 for sites MRT1 and MRT2, respectively, 51.06 mg·L−1 Ca for site GVR2 and 25.52 mg·L−1 Mg for site MRT4.

Three pairs of explanatory variables were highly correlated: conductivity-Ca (0.8), conductivity-water hardness (0.77), and water hardness-Ca (0.85). We decided to drop the variables conductivity and water hardness which were highly correlated to the Ca concentration and to keep the last one for the next step analysis.

3.2. Key Variables of Species and Sites Assemblages

The tb-RDA analysis (Table 2) using the list of previously retained quantitative and qualitative variables (pH, NO3, NH4, PO4, NO2, Ca, Mg, slope converted into numerical values from 1 to 3, canopy cover and orientation). The first unconstrained axis (PCA1) explains a rather high amount of variation, almost comparable to the first constrained axis, RDA1. This suggests that there is another environmental variable, not measured during this study, which could influence the species composition of the communities (see Section 4). Nevertheless the variables explain 36.5% of the variation and a Monte Carlo permutation test executed on this RDA delivers a highly significant value for the model () that confirms the global model is relevant.

Inside the constrained variance the first axis explained 44.9% (= eigenvalue RDA1 = 0.4282/0.9524) of this constrained variance and the second axis 18.1% (= eigenvalue RDA2 = 0.1721/0.9524). The use of the same dataset without the two outliers yields approximately the same results (38.3% of constrained variation, 42.7% for the first axis, and 15.4% for the second axis). As the parameter orientation is not retained by the forward selection model and because it produces a poorly readable figure, we also tested the same model without the orientation variable (Figure 2). The results are the following: 28.3% of constrained variation, 55.9% for the first axis, and 17.7% for the second axis.

The RDA analysis is potted both with and without the orientation variable (Figure 2). The scaling 2 option [37] allows the following comments and shows that slope, NO3, NH4, PO4, and canopy cover play an important role in the dispersion of the sites along the first axis.

The two variables NH4 and PO4 are very closely correlated, suggesting that they can act in a similar way on the species assemblages and between factors themselves; however, NO3 is negatively correlated with both NH4 and PO4. The two characteristic species Eucladium verticillatum and Palustriella commutata are close together suggesting they appear both on sites with higher slope and lower cover of canopy and lower NH4 and PO4 concentrations. Pellia endiviifolia shows a positive relationship with NO3 whereas Brachythecium rivulare and Oxyrrhynchium hians (Hedw.) Loeske are weakly correlated with PO4, NH4, and high pH values. Bryum pseudotriquetrum is negatively correlated with canopy cover. Most other species are clustered together away from these extremes. They show mostly shorter projections, indicating that they are either present in most of the sites or related to intermediate ecological conditions; however, they could also be related to some variables on third or higher axis which is not displayed.

The result of the forward model selection of the RDA retained a model with (sorted by decreasing importance) slope, Mg, Ca, and canopy cover (spe ~ slope + Mg + Ca + cover). The same analysis without the outliers (“Haumont5” and “Givry2”) generates a slightly different model with slope, canopy cover, Mg, and Ca as the most significant factors (spe ~ slope + cover + Mg + Ca).

3.3. Species Richness

The results show (Table 3) that only two explanatory variables (canopy cover and PO4) are strongly supported by the data (variable weight = 0.992 and 0.960, resp.) while two other ones are moderately supported by the data (NH4 and Mg) with a weight of, respectively, 0.709 and 0.650. All but the Mg variables are negatively correlated with species richness.

To check the potential influence of outliers we performed the same analysis after excluding Haumont5 and Givry2 sites which have very high PO4 and NH4 values. The model selection results are comparable to the first analysis including outliers; AIC retains only canopy cover and PO4 as the main explanatory variables (Table 3) (columns on the right for each parameter).

Figures 3 and 4 show the results of two predicted scenarios based on the complete model (2 outliers included) with the mean values of other variables. Figure 3 represents the relation between the PO4 concentration and the species richness and Figure 4 represents the relation between the canopy cover and the species richness. It can be seen that for the canopy cover the expected number of species rises by 2 or 3 species (going approximately from 2 to 4 or even 5) when the cover falls from 100 to 0%. In the case of the PO4 concentration the expected number of species falls to zero when the PO4 concentration rises from 0 to 0.5 mg·L−1 or more.

3.4. Presence/Absence of Characteristic and Companion Species

For each species we began by excluding the variables that were highly correlated using the VIF (Variance Inflation Factor) values. The only variable to exclude from the five species was orientation and thus the binomial GLM analysis included the following explanatory variables: zone, NO2, Mg, pH_mean, Ca, NH4, NO3, canopy cover, slope, and PO4.

3.4.1. Palustriella commutata

The most important variables delivered by the AIC procedure are “steep slope” and “moderate slope” and to a lesser extend “canopy cover” (Table 4). The frequency of Palustriella commutata tends to increase when the slope is steeper and decreases when canopy cover is higher (Figure 5).

3.4.2. Eucladium verticillatum

The most important variables delivered by the AIC procedure are canopy cover and NH4, both are negatively correlated with the probability of the presence of the species (Table 5). The frequency of Eucladium verticillatum tends to decrease when the canopy cover increases and for higher concentrations of NH4, however the size of the effect is relatively limited as shown on Figure 6.

3.4.3. Cratoneuron filicinum

The only important variable delivered by the AIC procedure is the Mg concentration (positive relationship) and to a lesser degree NH4 (negative relationship) (Table 6, Figure 7). The results remain similar after excluding the “Laarbeek (Kwel)” site whose outlying value (65.5 mg·L−1) could have unduly influenced the correlation.

3.4.4. Pellia endiviifolia

The most important variables delivered by the AIC procedure are PO4 (negative relationship) and NO3 (positive relationship) (Table 7; Figure 8).

3.4.5. Brachythecium rivulare

The only important variable delivered by the AIC procedure is the Ca concentration (positive relationship) (Table 8). However Figure 9 shows that this correlation is weak and is based almost entirely on the outlying value (251 mg·L−1) of the “Bois de Dieleghem” site. It would thus be sensible to not pay too much attention to this correlation.

4. Discussion

This study aimed to identify the driving variables that explain the species richness and the distribution of the bryophyte assemblages within the 7220 habitat. The environmental driving factors discussed within this study apply only to the 7220 habitat of travertine forming springs and watercourses, and it would not be advisable to apply these findings outside of this habitat such as (less calcareous) types of water courses in Wallonia and Brussels-Capital Region. The springs and watercourses in this study are characterized by medium to high levels of Ca (60–300 mg·L−1) [40] and high pH values (7.5–8.5). The extreme values in our samples range from 76.5 to 251 mg·L−1 for Ca (the 27.8 value seeming doubtful) and pH from 6.2 to 8.2. In comparison Brusa and Cerabolini [11] recorded Ca values between 33.2 and 60.2 mg·L−1 and pH values between 7.5 and 8.4.

4.1. Global Analysis

The results of the tb-RDA showed that only 36.5% of the total variance in the dataset can be explained by the environmental variables, considered within this study. This suggests that other variables not measured here could explain better the observed species assemblages. Following other authors [41, 42], we suggest that factors such as the substrate itself (sand, loam, and rock), the micro topography [11, 16, 18, 33, 43], the water turbulence [2, 8, 44], the distance from the spring [8, 13], the presence and characteristics of travertine-productive microorganisms like cyanobacteria [8, 16, 45], and the local relative humidity could play an important role in the species assemblages even if they would probably be difficult to quantify.

The results show that the most important variables explaining the rest of the variance are the slope, NO3, NH4, and PO4 concentrations (Figure 2). The slope has been identified by many authors as a key factor for the presence of the travertine communities [18] mainly because a steep slope generates more turbulence that enhances CO2 degassing and thus calcite precipitation that accounts for more than 80% of the precipitation in travertine [13, 45]. The negative relationship between the PO4 concentration in the spring water and the presence of mosses specific to petrifying sources should be further explored but could possibly be explained by the fact that PO4 has an inhibitive role in the formation of calcite [3, 46, 47]. Pentecost [13] also highlights the negative role of phosphate pollution on the deposits formation and on the travertine communities. These results are partly confirmed by the forward selection that has similarly retained slope and canopy cover but also Mg and Ca (as well as NH4 when the 2 possible outliers are dropped). That means that the relief, the canopy cover, and the degree of eutrophication play an important role among the variables. It is not surprising that NH4 and PO4 variables act in the same way (Figure 2). More curiously though is the opposite relation between these two variables and the NO3 concentration though a similar observation has been made by Denys and Oosterlynck [48] in Flanders. The second axis showed that the canopy cover also plays an important role in the dispersion of the species assemblages. This result is expected since we know that the exposure and the degree of relative humidity play a key role in the ecology of most of bryophytes species [49].

The analysis also indicates that some species are more specifically linked to some variables considered in this study. Eucladium verticillatum and Palustriella commutata appear mostly in the same kind of locations characterized by a moderate or steep slope and with no or very low levels of eutrophication. This result is similar to analysis in Italy, using the same multivariate analysis [11]. Canopy cover is also negatively correlated with these sites indicating that these species occur more frequently at enlightened sites. Pellia endiviifolia shows a tendency to grow in NO3-rich sites, which is confirmed by other observation that the species can grow in oligotrophic and eutrophic waters [40]. Finally Brachythecium rivulare and Oxyrrhynchium hians are more likely to be found in more eutrophicated sites (determined by PO4 and NH4 values).

4.2. Species Richness

The Poisson GLM analysis and the AIC have retained two major significant variables concerning the species richness, that is, canopy cover and PO4 concentration. However other variables that have not been included in this analysis, such as soil type and water turbulence, should be considered in future studies. For example, Pentecost [13] has stated that two other variables, site size and heterogeneity, are also important factors affecting species richness. It can be seen in Figure 4 that all other registered variables being equal there could be a gain of only 2 or 3 species when clearing completely a theoretical site. This gain does not predict if the gained species would be characteristic species or other species.

The negative relation with PO4, and with NH4 when keeping the outliers in the dataset, indicates that the more oligotrophic the conditions the more rich the community. Other studies indicate that this kind of relationship is not always attributed to eutrophication, and human impacts are sometimes positively correlated to species richness [48]. These studies may address a wider range of conditions than in the present study and suggest that this relationship is scale-dependent. The majority of the richest sites were characterized by PO4 concentration < 0.1 mg·L−1 which is very close to the 0.05 mg·L−1 threshold determined by Galoux et al. [41] for the reference conditions of Walloon watercourses. Brusa and Cerabolini [11] in Italy using multivariate analysis found that the 3 major factors explaining their species assemblages were the permanency/nonpermanency of the water flow, the slope, and the canopy cover. In this study the slope is not retained in the species richness analysis; however, it is an important variable for Palustriella commutata.

4.3. Characteristic and Companion Species
4.3.1. Palustriella commutata

Our analysis shows that slope and to a lesser extend canopy cover are important variables explaining the presence of Palustriella commutata in the surveyed sites. The absence of canopy combined with a steep slope increases the probability of presence of Palustriella commutata (Figure 5). The slope has been detected as the only significant variable by Brusa and Cerabolini [11]. The fact that slope seems to play an important role in the probability of presence of the species seems contradictory with the well-known occurrence of the species in calcareous meadows in Belgium [16, 28]. It is possible that the properties of the water including highly mineralized or calcareous waters and permanently wet surfaces could be even more important than slope [11, 19, 28, 50]. It has been mentioned by many authors [13, 43, 51, 52] that there are two subspecies or species [14] Palustriella commutata and P. falcata whose ecological preferences are slightly different. According to these authors P. commutata occurs relatively more often on boggy sites and around small springs and seepages with a higher pH and higher Ca concentrations, and it tolerates higher rates of nitrates and grows almost exclusively on limestone. P. falcata seems more restricted to open habitat of turbulent water at higher altitude and in less bases-rich waters. Concerning the canopy cover other authors [1, 13, 18, 50] mention the fact that this factor must play a role even if the species tolerates much shade. The species seems also to tolerate a certain amount of eutrophication [18, 34, 50] even if it grows mostly in oligotrophic circumstances.

4.3.2. Eucladium verticillatum

Although canopy cover and NH4 are both negatively linked with the probability of presence of the species these variables seem not to play a major role. It is probable that other variables could also influence species richness including less permanently wet substrate by capillarity or because of being splashed. Indeed we found the species abundant around a spring on a steep sandy substrate in the Rouge-Cloître in the Brussels-Capital Region that is not splashed at all but is permanently wet due to capillarity and high relative humidity. The relative weak negative relationship with canopy cover is mentioned by some authors [2, 24, 50] that found the species in much shaded habitats.

4.3.3. Cratoneuron filicinum

This species is not recognized in Wallonia as a characteristic species of the 7220 habitat as it is widely distributed [19], occupying calcareous meadows, calcareous springs, and wet rocks, as well as forest tracks where dolomite gravels have been added to harden the soil. In the Walloon watercourses Cratoneuron filicinum has been identified as one of the few indicative species of more or less calcareous streams [41]. In our sites the species has a relatively wide ecological amplitude as regards canopy cover and light to moderate eutrophication. The tolerance to a certain level of eutrophication is also mentioned by Bailly et al. [50] and by Sossey-Alaoui and Rosillon [40]. Mg proves to be a key-element in our sites and this observation should be related to the fact that the species is especially common on dolomite substrates, such as forests tracks that are rich in Mg, as well as on marl substrates in Lorraine [19].

4.3.4. Pellia endiviifolia

This is a widespread species in Wallonia and Brussels [19, 53] commonly associated with neutral to basic waters. It is also known to tolerate shade as well as moderately eutrophic waters [50]. It is therefore not surprising that none of these variables show clear relationships with the probability of presence in our sites. The apparent contradictory relationships with PO4 and NO3 are hard to clarify. Sossey-Alaoui and Rosillon [40] indicate that the species can be an indicator species both in oligotrophic streams and in nitrates-rich streams, providing that the waters are mineral-rich (high conductivity/high Ca concentration).

4.3.5. Brachythecium rivulare

Together with Pellia endiviifolia this species is relatively widespread in Wallonia and Brussels-Capital Region [19, 53]. It occurs along water courses and in marshes, though it is less linked to calcareous waters. In Wallonia it has been demonstrated that it has very wide ecological amplitude that is common in every natural region independent of the substrate [41]. It has also wide ecological amplitude occurring in a range of trophic levels and tolerating shade but preferring open sites [50]. This high ecological amplitude mirrors in the quasi-absence of significant variables explaining the presence of the species in our dataset.

4.4. Management Measures

One of the main objectives of this study was to define management objectives that could be implemented to maintain or improve the conservation status of the habitat 7220. As only relatively small proportion of the observed variability can be explained by the surveyed variables; the inclusion of further variables is recommended for future studies. Additionally among the surveyed variables some cannot be managed at all like the slope and Mg concentrations. Nevertheless, our results provide information on key variables that should be considered when managing these sites. The management measures may be different depending on their intention and may vary if the aim is to enhance biodiversity or favouring particular species. Generally speaking, we have seen that one of the most prominent results is that the species assemblage characteristic of the 7220 habitat (Palustriella commutata and Eucladium verticillatum) favours very low PO4 concentration (PO4 < 0.2 mg·L−1) and low canopy cover. This is also the case but to a lesser degree with NH4. Graham and Farr [34] have stated that these two species are also mostly present when PO4 concentration is less than 0.05 mg·L−1. Brachythecium rivulare, Oxyrrhynchium hians, and even Pellia endiviifolia can occur in eutrophic conditions; however, in these conditions there is a lower probability of presence of the two characteristic species.

The more open the site the more diverse the moss assemblages will be even if the expected number of additional species is only 2 or 3 (Figure 4); however, the favoured species would not automatically be characteristic species of the 7220 habitat. If canopy cover plays a role in favouring species richness, it would be advisable not to cut down all existing trees and bushes but to undertake thinning to avoid possible radical negative changes in relative humidity and temperature on the bryophyte communities.

We recommend that regular monitoring of the water chemistry is the best way to detect abnormal changes especially eutrophication (PO4, NH4, NO3). Monitoring regimes are already in place for most of the Walloon and Brussels sites but should be extended to all 7220 known sites.

4.5. Suggestions for Future Studies

To try to intercept other significant factors in the distribution and the presence/absence of the most characteristic species of the 7220 habitat, further studies should seek to take into account the following variables: type of substrate (sandy, calcareous, silty, etc.), the micro topography, the water turbulence, the distance from the spring, the relative humidity, and the presence of algae/cyanobacteria.

5. Conclusions

In this study, we highlighted the variables explaining the species richness and the distribution of moss assemblages of the 7220 habitat in the southern part of Belgium. The driving factors highlighted in this study can only be applied to the 7220 habitat and would probably not be valuable for studies at sites dominated by less calcareous waters in Wallonia and Brussels-Capital Region. Our study recommends that monitoring of water chemistry (eutrophication) should be applied to all 7220 habitats, and this information should be used to support practical management actions to maintain or achieve a good status of conservation.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to thank their colleagues Fabrice Etienne and Etienne Peiffer who collected the water samples and recoded the site variables. Other colleagues helped them when measuring the parameter, Daniel Galoux and Jonathan Devriese. Christine Keulen provided them with useful information about the thresholds of the different water courses in Wallonia. Mathias Engelbeen and Wim Van Den Eynden gave them access to the sites of the Brussels-Capital Region. They are also very grateful to David Zeleny who kindly accepted to check their multivariate analysis. They thank Gareth Farr (British Geological Survey) for undertaking a review of the English manuscript. Finally, they would like to thank Pierre Gerard who allowed them to carry out this study.