Benthic Diatom Communities of a Large Mediterranean River under the Influence of a Thermal Effluent ()
1. Introduction
Diatoms are unicellular algae with a wide spectrum of responses to seasonal and environmental variation, and with optimum ranges of temperature to grow [1] [2] [3] . Each species has different tolerances and preferences, and some have therefore been used as indicators of environmental changes and conditions [4] [5] [6] .
Temperature has a significant role in all biochemical and physiological functions of organisms and influences the morphology, physiology, behavior, growth, reproduction, and distribution of species [7] [8] , and it has been noted as a main factor influencing primary production [9] . In fact, the rate of photosynthesis depends directly on temperature because it is an enzyme controlled process [10] . It has been reported that warming generally increases the primary production [8] .
The importance of temperature in rivers has been widely recognized [9] [11] [12] , and the effects of its alteration on aquatic species cover a wide spectrum of direct and indirect effects that range from minor importance to lethal effects [13] . Changes in community structure as response to thermal disturbances have been detected even with a temperature alteration of a few degrees [14] and depend on the preference and tolerance of species to different temperatures as well as on the level of heating.
To generate thermal power, nuclear power stations use nuclear fission to heat water and drive steam turbines that then produce electricity; but this process requires large volumes of water for its cooling system in order to remove the waste heat produced. The increase in river water temperature caused by these thermal discharges has been shown to alter biological and ecological components of aquatic ecosystems [11] [15] [16] , but the effects are variable, and depend on the levels and quantity of heated discharge and on the biological features of the environment [17] [18] . Depending on the design and the operating units of the power plants, water temperature in effluent sites can increase by as much as 8˚C [19] . However, in Europe, legislation requires that the temperature downstream of the effluent should not increase by more than 3˚C [20] .
Many authors have studied the ecological effects of temperature in aquatic environments [1] [13] [15] , and several such studies have been based on diatoms [2] [7] [10] [21] . Furthermore, diatoms have been also used as biological proxy to assess the effects of climate change [5] [22] [23] [24] . The impacts of thermal effluents on benthic diatom communities have been studied mostly in estuarine and coastal regions [17] [18] [25] [26] [27] , and to a lesser extent in lakes, rivers and streams [28] [29] [30] [31] [32] . However, literature dealing with the effects of thermal pollution on benthic diatom communities of Mediterranean rivers is absent, even though this type of alteration is frequent in the watersheds of the Mediterranean basin.
This study aimed to assess changes in the community structure of benthic diatoms caused by the thermal pollution produced by the cooling system of a nuclear power station (Ascó nuclear power station). This is one of the main anthropogenic factors exerting pressure on the lower Ebro River and has been subjecting the river to a sustained heating during the last 30 years, therefore providing an excellent opportunity for assessing the long-term effects of water warming on benthic communities. For this purpose, surveys at sites located before and after the effluent were conducted, and to minimize the potential influence of substrate heterogeneity, artificial substrates deployed over the same temperature gradient than natural surfaces were also analyzed.
2. Materials and Methods
2.1. Study Area
The Ebro River, located in the NE of the Iberian Peninsula (Figure 1(a)) has a
(a) (b)
Figure 1. (a) Map of the lower Ebro River showing the study area; (b) Location of sampling sites.
length of 928 km; its basin has a surface of 85,534 km2 being one of the most important tributaries to the Mediterranean Sea. Over 180 dams regulate the river flow and the lower part is regulated by two large reservoirs (Mequinensa with a capacity of 1534 hm3 and Riba-roja with a capacity of 207 hm3) built in 1964 and 1969, respectively for hydropower purposes. The last downstream dam is located at Flix, a small reservoir with a capacity of 11.4 hm3.
The Ascó nuclear power station is located at the right margin of the lower Ebro River, 10 km downstream the Flix dam, between Ascó and Flix towns, and at about 110 km from the river delta (Figure 1(b)). It was built in 1984 and has two reactors with a gross electrical power output of about 2050 MWe and a thermal reactor power of about 5900 MWt. (data available at http://www.anav.es). The power station has granted a concession of 72.3 m3/s of the Ebro’s flow for its cooling system, and a weir has been built to collect the river water to the condensers. After its use the water is returned to the river with an average thermal increase of 3˚C [33] .
The study area has a total length of 2 km that comprise 1 km before and after the nuclear power station (41˚12'0"N, 0˚34'10"E).
2.2. Diatom Sampling and Preparation
In order to compare benthic community features of a site unimpacted by the heated effluent with those under its influence, three sampling sites were selected: a control site (C), located upstream the nuclear power station, and two impacted sites (I1 and I2) covering the thermal plume, located downstream of the effluent outlet, on the right and left river margins, respectively (Figure 1(b)).
Three sampling campaigns were conducted in August, October and December of 2013. In every occasion, three replicates were collected at each site from both natural substrata (pebbles) and artificial substrata (fired clay bricks placed with a colonization period of 6 weeks). During the summer campaign the artificial substrates placed on site I1 were not recovered due to vandalism.
For every sampling site and occasion, physicochemical data were recorded. A YSI 556 multi-parameter probe was used to measure dissolved oxygen (mg/l), oxygen saturation (%), pH, salinity (ppt) and conductivity (mS/cm); current velocity at 60% of total water depth was recorded with a Braystoke BFM 001 current meter; total dissolved nitrogen (TDN), total nitrogen (TN), soluble reactive phosphorus (SRP), total phosphorus (TP) were measured according to Koroleff (1977, 1983) [34] [35] ; and the total chlorophyll concentration was calculated using the colorimetric method [36] . In each sampling site, water temperature (˚C) was monitored at intervals of 30 minutes during all the study period with a TCtemp1000 data logger Madgetech.
Benthic diatom samples were collected according to the recommendations of Kelly et al. (1998) [37] . The suspension was fixed in 4% formaldehyde solution. At the laboratory, benthic samples were oxidized with H2O2 30% v/v for several hours in order to remove the organic matter. HCl 37% v/v was added to evaporate the carbonates from the samples, as described in Renberg (1990) [38] . Clean valves were permanently mounted with Naphrax© (refractive index 1.74). The permanent slides were examined using a LEICA DMI 3000 B light microscope equipped with differential interference contrast (DIC) with a 100 times oil immersion objective (n.a = 1.40).
Identification of diatoms was done to species level mainly following Krammer and Lange-Bertalot (1986-1991) [39] but other taxonomic and floristic works were also used when needed. A minimum of 400 valves were counted each time.
2.3. Data Analysis
Water temperature values recorded over the study period were analyzed to identify variations and trends, the difference of temperature between control and impacted sites was calculated (Diff_T) and the temperature variability at each site was represented by the standard deviation values (TempSD).
Differences in values of environmental variables between sites were tested with analysis of variance (ANOVA) with Tukey post hoc test performed using software SPSS 19 (SPSS Inc., Chicago, IL, USA).
Diatom abundance is presented as relative percentages and it was square-root transformed in order to reduce the effect of highly variable population densities on ordination scores. All environmental variables that expressed concentration were logarithmically transformed before analysis to avoid skewed distributions.
Descriptive community parameters were calculated: Richness (S), Shannon-Wiener’s diversity index (H’) H’ = −SUM (Pi*Loge(Pi)), and Pielou’s evenness index (J’) J’ = H/log(S).
Sites were ordered in relation to their species composition using Non-metric Multidimensional Scaling (NMDS) and significant differences were identified using 1-way Analysis of Similarities test (ANOSIM) [40] , that hypothesizes for differences between groups of samples (defined a priori) through randomization methods on a resemblance matrix; ANOSIM provides an R statistic value that reflects the amount of dissimilarity associated with each group, R values close to one indicate very different composition, whereas values near to zero indicate little difference. Then, in order to identify resemblances between sample groups and to identify taxa that contributed to dissimilarity among sites, a Similarity Percentage Analysis (SIMPER) was performed.
Finally, relationship between diatom assemblages and environmental variables was assessed with a multivariate distanced-based linear regression model (DISTLM) [41] and a set of explanatory variables was identified. The model was visualized through a distance-based redundancy analysis (dbRDA) performed using PRIMER V6 software [42] with the add-on package PERMANOVA+ [43] .
3. Results
3.1. Environmental Characteristics
The average values for physicochemical parameters measured at each sampling site are shown in Table 1. Water temperature showed constantly higher values at impacted sites as consequence of the water heating produced by the cooling system of the nuclear power station (Figure 2), and was significantly different between control and impacted sites (ANOVA p = 0.008) (C ≠ I1, C ≠ I2, I1 = I2). The mean values recorded over the study period were 20.54˚C (C), 23.04˚C (I1) and 22.98˚C (I2); while the mean difference of temperature recorded between C and I1 was 2.39˚C and 2.33˚C between C and I2. Water velocity showed mean values of 0.26 m/s at control site, and 0.13 m/s and 0.11 m/s at I1and I2 respectively; significant differences between control and impacted sites were found (ANOVA p = 0.000) (C ≠ I1, C ≠ I2, I1 = I2). The other measured environmental variables (dissolved oxygen, pH, conductivity, soluble reactive phosphorus,
Table 1. Values of physicochemical parameters measured at each sampling site. (T = temperature, Diff. T = temperature difference, Temp SD = temperature variability, DO = dissolved oxygen, Cond = conductivity, SPR = soluble reactive phosphate, TP = total phosphorus, TDN = total dissolved nitrogen, TN = total nitrogen, Chl a = chlorophyll a).
Figure 2. Water temperature recorded over the study period at control (C) and impacted (I1, I2) sites.
total phosphorus, total dissolved nitrogen, total nitrogen and depth) only showed minor variation and did not present significant differences between sites (ANOVA p > 0.05). Dissolved oxygen at control site showed highest values than impacted sites in winter (10.23 mg/l) and in summer (8.46 mg/l) but in autumn I2 showed the highest value (7.71 mg/l); pH showed variation from 7.8 to 8.1; conductivity varied from 0.84 mS/cm in summer (C) to 1.31 mS/cm in winter (I2); soluble reactive phosphorus showed a minimum of 29.6 μg/l (C, winter) and a maximum of 53 μg/l (I1, summer); total phosphorus ranged from 111.5 μg/l (C, winter) to 598 μg/l (I2, summer); total dissolved nitrogen varied from 1319.4 μg/l (I1, autumn) to 1712.2 μg/l (I1, winter); total nitrogen showed a minimum of 1999.9 μg/l (I1, autumn) and a maximum of 3120.4 μg/l (I1, winter); and the mean depth values recorded were 0.74 m (C), 0.74 m (I1) and 0.80 m (I2).
3.2. Diatom Assemblages
During the study period a total of 85 species in natural substrate and 78 species in artificial substrate were found and are listed with their taxon authors and relative abundances in Table A1. Seasonal changes were observed in the diatom community along the study period. In natural substrate assemblages, Amphora pediculus and Nitzschia inconspicua were the dominant species, sharing this dominance with Navicula capitatoradiata in summer, with Amphora copulata in autumn, and with Reimeria uniseriata in winter. Artificial substrate assemblages were dominated in summer by Nitzschia inconspicua, N. palea and Cocconeis placentula var. lineata; in autumn by Amphora pediculus, Cocconeis placentula var. euglypta and Nitzschia inconspicua; and in winter by Amphora pediculus, Cocconeis placentula var. lineata and C. placentula var. trilineata.
Concerning diatom diversity (Table A2), there were no significant differences (ANOVA p > 0.05) between control and impacted sites. However, when mean annual values were analyzed, slightly higher values of species richness and diversity indices were found at impacted sites (Table A3).
The NMDS ordination (Figure 3) displays the spatial distribution of the control and impacted sites and the stress obtained was 0.18 and 0.17 for natural and artificial substrates, respectively. For both types of substrate, the assemblage composition was analyzed with ANOSIM and showed significant differences between C and I1, and between C and I2; but not between I1 and I2 (Table 2).
The Simper analysis (Table A3) for natural substrate assemblages showed that the mean dissimilarity between control and impacted sites was 42.22% and Amphora pediculus, Nitzschia inconspicua and Navicula capitatoradiata were the species with highest percentage of contribution to dissimilarity between groups. While for artificial substrate, the mean dissimilarly was 39.97% and the species with the highest contribution were Amphora pediculus, Nitzschia inconspicua and Rhoicosphenia abbreviata.
The dbRDA analysis performed for natural substrate (Figure 4), revealed that
Table 2. Values of R statistic and significance level of differences between C, I1 and I2 groups, obtained by ANOSIM test for diatom communities of natural and artificial substrate.
Significance: *p ≤ 0.05;**p ≤ 0.01; ***p ≤ 0.001.
Figure 3. Two dimensional NMDS plots based on Bray-Curtis similarities of square-root transformed diatom abundance data. (a) Natural substrate ordination; (b) Artificial substrate ordination.
(a)(b)
Figure 4. Distance-based redundancy analysis (dbRDA) ordination of natural substrate data: (a) Samples displayed by site and season and vectors showing correlation between explaining variables and dbRDA axes; (b) Samples displayed by season and site and vectors showing correlation between the five species with highest contribution to the dissimilarity between control and impacted sites and dbRDA axes. (NCPR = Navicula capitatoradiata, APED = Amphora pediculus, NINC = Nitzschia inconspicua, RUNI = Reimeria uniseriata, CPLT9 = Cocconeis placentula var. trilineata).
the set of variables selected by the DISTLM (T˚, total nitrogen, T˚ difference, chlorophyll and T˚ variability) explained 57.27% of fitted variation and 26.81% of total variation in the two first axes; while the dbRDA performed on artificial substrate (Figure 5), revealed that the set of variables selected by the DISTLM (dissolved oxygen, T˚ difference, total phosphorus, pH and chlorophyll) explained 67.24% of fitted variation and 39.55% of total variation in the first two axes. Water velocity was not selected by the DISTLM as part of the explanatory variables set.
The first axis of the dbRDA plots of both natural (Figure 4) and artificial (Figure 5) substrates distinguished samples from control and impacted sites and in both cases the axis was strongly correlated with the difference in water temperature caused by Ascó nuclear power station. The second dbRDA axis, also of both natural (Figure 4) and artificial (Figure 5) substrates, basically distinguished autumn and winter samples from those of summer and was strongly correlated with a gradient of temperature and nutrient levels (total nitrogen for natural substrate and total phosphorus for artificial substrate) associated with the seasonal variation in the fluvial system. The five species with the highest contribution to the dissimilarity between control and impacted sites are represented in the dbRDA plots (Figure 4 and Figure 5).
4. Discussion
The presence of the nuclear power station has been producing a sustained increase of water temperature occurring over the last 30 years at the Ebro River.
The values recorded for environmental variables and the distribution of samples indicated a seasonal variation explained by the natural fluctuation of conductivity, pH, dissolved oxygen and nutrients related with changes in temperature and flow as consequence of the annual cycle and hydrodynamics of the river. However, this process has been altered by a thermal increase consistently greater than 2˚C in the river, and the results obtained in this study show the existence of two different diatom assemblages inhabiting in sites before and after the nuclear power station.
Most of environmental variables measured shared the same values or showed a little variation between control and impacted sites; therefore, the differences detected in diatoms assemblages could be mostly attributable to the warming effect, either by its direct influence or by its interaction with other functional processes. The sensitivity of diatoms to changes in water temperature is widely recognized [44] [45] [46] . Increases in temperature have complex effects, for instance affecting the diffusion rates of chemicals and reducing the amount of oxygen that water may maintain; these changes in the environmental conditions will very likely affect the reproductive rates and metabolism of the algae [47] [48] and therefore lead to changes in community structure.
Diatom assemblages were significantly different between control and impacted sites; these differences were mainly due to variation in community composition expressed as species abundances rather than species presence or absence.
(a)(b)
Figure 5. Distance-based Redundancy Analysis (dbRDA) ordination of artificial substrate data: (a) Samples displayed by site and season and vectors showing correlation between explaining variables and dbRDA axes; (b) Samples displayed by season and site and vectors showing correlation between the five species with highest contribution to the dissimilarity between control and impacted sites and dbRDA axes. (APED = Amphora pediculus, NINC = Nitzschia inconspicua, NPAL = Nitzschia palea, RABB = Rhoicosphenia abbreviata, CPLT9 = Cocconeis placentula var. trilineata).
These changes in abundance could be related to specific physiological responses of species to their optimal temperature ranges, but may also be related to shifts as consequence of interspecies interactions as competition or due to the influence of other environmental variables. Our data do not allow to attribute the observed changes in community structure solely to the temperature alteration, but evidence that warming is a determinant factor influencing or enhancing other factors on the structure of communities.
In this study, the species pool did not show significant variation; we found slightly higher values in species richness and diversity indices at the impacted sites. For algae, it has been documented that diversity increases from 0˚C to 25˚C and starts to decrease at temperatures above 30˚C [3] [8] [9] , though changes in community structure are usually more evident at temperatures from 25˚C to 30˚C rather than <25˚C [7] . During the study period, temperature never exceeded 25˚C and changes detected in species composition were minor and due to species with low relative abundances (<5%). Similar results were reported in a study including benthic epilithic communities under thermal influence, where Hillebrand et al. (2010) [26] found that elevating the water temperature increased temporal beta-diversity and decreased compositional stability of communities; and instead of changes in species richness, it was observed a change of proportion of species from the same pool. Changes in diatom community structure as consequence of thermal alteration were documented by Squires et al. (1979) [31] who found that the algal flora was significantly affected at the section immediately below the discharge point of a power station; and Vinson and Rushforth (1989) [32] noted that diversity and species richness increased with temperature and maximum values were reached between 25˚C and 30˚C, beyond this temperature species diversity decreased; parallel results were also found by Patrick (1969) [2] .
The colonization on artificial substrate seemed to be dominated by opportunistic diatom species with fast growth rates such as Amphora pediculus, Nitzschia inconspicua and Rhoicosphenia abbreviata, which can quickly form large blooms and compete with other algal species with slower growth rates, as has been previously highlighted by Snoeijs and Prentice (1989) [27] . These species also showed high abundance on natural substrate and in both cases (natural and artificial substrate), dominance was shared with Cocconeis spp., which did not show a clear preference between control and impacted sites. These results are opposite to those found by Stevenson et al. (1996) [48] who detected a shift to dominance of Cocconeis in warmer waters but this discrepancy could be explained, perhaps, by the results of De Nicola (1996) [7] , who noted that Cocconeis tended to be more abundant in waters above 25˚C and as mentioned before, we did not record values exceeding that temperature. Interestingly, although there are some community differences between natural and artificial substrate, both provided essentially the same picture of thermal influence. This agrees with some previous works where again it was found that benthic diatom communities tend to be much more affected by the environmental conditions than by substrate type [49] [50] [51] .
Diatom communities proved to be sensitive to water warming even though this alteration did not exceed 3˚C. The factors that seemed to have most effect on the benthic assemblages inhabiting the area influenced by the nuclear power station, were the seasonal variation and thermal alteration caused by the heated effluent.
Nowadays, aquatic ecosystems are threatened as consequence of greater water demands and climate change, and to ensure their adequate management, it is evident the need to better understand the response of biota to thermal alterations. This is especially important for Mediterranean ecosystems, since this region of the world is going to be among the most impacted ones by climate change, and in particular by global warming [52] . By using the thermal gradient in the nuclear power station plume it is possible to cover part of the range of future scenarios of temperature and therefore, our results could be of interest to predict changes in benthic communities under global warming scenarios. However, it is important to note that this work focused on a local species pool along a period of time when the regional species pool did not change, and it has been pointed out that global warming will lead to turnover also in the regional species pools by actions as emigration or adaptation of species from other regions [53] ; thus, changes in local species could also be influenced by changes on regional scales if the temperature increase affect larger areas over longer periods [26] .
We think that the information generated here will contribute to a better understanding of the effects of increasing temperature on the benthic diatom communities of Mediterranean rivers and hence will provide useful baseline data for predicting the effects of global warming under future projected scenarios.
Acknowledgements
This study was funded by the Secretaría de Educación Superior, Ciencia, Tecnología e Innovación (SENESCYT) of Ecuador, which provided a doctoral research fellowship to the first author. Special thanks are given to David Mateu and Lluís Jornet (IRTA) for their invaluable help in the sampling campaigns.
Appendix
Table A1. List of diatom taxa found and their relative abundances (%) over the study period at control (C), and impacted sites (I1 and I2).
Table A2. Diatom community descriptive parameters for over the study period at control (C) and impacted sites. Richness (S), Shannon-Wiener’s diversity index (H’, as log2) and Pielou’s evenness index (J’).
(a) (b)
Table A3. Similarity percentages analysis (SIMPER) of diatom taxa showing mean dissimilarity between control and impacted sites and percentages of taxa contribution until reach 50%. (a) Natural substrate; (b) Artificial substrate.